Comprehensive analysis of the REST transcription factor regulatory networks in IDH mutant and IDH wild-type glioma cell lines and tumors.

The RE1-silencing transcription factor (REST) acts either as a repressor or activator of transcription depending on the genomic and cellular context. REST is a key player in brain cell differentiation by inducing chromatin modifications, including DNA methylation, in a proximity of its binding sites. Its dysfunction may contribute to oncogenesis. Mutations in IDH1/2 significantly change the epigenome contributing to blockade of cell differentiation and glioma development. We aimed at defining how REST modulates gene activation and repression in the context of the IDH mutation-related phenotype in gliomas. We studied the effects of REST knockdown, genome wide occurrence of REST binding sites, and DNA methylation of REST motifs in IDH wild type and IDH mutant gliomas. We found that REST target genes, REST binding patterns, and TF motif occurrence proximal to REST binding sites differed in IDH wild-type and mutant gliomas. Among differentially expressed REST targets were genes involved in glial cell differentiation and extracellular matrix organization, some of which were differentially methylated at promoters or gene bodies. REST knockdown differently impacted invasion of the parental or IDH1 mutant glioma cells. The canonical REST-repressed gene targets showed significant correlation with the GBM NPC-like cellular state. Interestingly, results of REST or KAISO silencing suggested the interplay between these TFs in regulation of REST-activated and repressed targets. The identified gene regulatory networks and putative REST cooperativity with other TFs, such as KAISO, show distinct REST target regulatory networks in IDH-WT and IDH-MUT gliomas, without concomitant DNA methylation changes. We conclude that REST could be an important therapeutic target in gliomas.

SorLA restricts TNFα release from microglia to shape a glioma-supportive brain microenvironment.

SorLA, encoded by the gene SORL1, is an intracellular sorting receptor of the VPS10P domain receptor gene family. Although SorLA is best recognized for its ability to shuttle target proteins between intracellular compartments in neurons, recent data suggest that also its microglial expression can be of high relevance for the pathogenesis of brain diseases, including glioblastoma (GBM). Here, we interrogated the impact of SorLA on the functional properties of glioma-associated microglia and macrophages (GAMs). In the GBM microenvironment, GAMs are re-programmed and lose the ability to elicit anti-tumor responses. Instead, they acquire a glioma-supporting phenotype, which is a key mechanism promoting glioma progression. Our re-analysis of published scRNA-seq data from GBM patients revealed that functional phenotypes of GAMs are linked to the level of SORL1 expression, which was further confirmed using in vitro models. Moreover, we demonstrate that SorLA restrains secretion of TNFα from microglia to restrict the inflammatory potential of these cells. Finally, we show that loss of SorLA exacerbates the pro-inflammatory response of microglia in the murine model of glioma and suppresses tumor growth.

Mapping chromatin accessibility and active regulatory elements reveals pathological mechanisms in human gliomas.

Chromatin structure and accessibility, and combinatorial binding of transcription factors to regulatory elements in genomic DNA control transcription. Genetic variations in genes encoding histones, epigenetics-related enzymes or modifiers affect chromatin structure/dynamics and result in alterations in gene expression contributing to cancer development or progression. Gliomas are brain tumors frequently associated with epigenetics-related gene deregulation. We perform whole-genome mapping of chromatin accessibility, histone modifications, DNA methylation patterns and transcriptome analysis simultaneously in multiple tumor samples to unravel epigenetic dysfunctions driving gliomagenesis. Based on the results of the integrative analysis of the acquired profiles, we create an atlas of active enhancers and promoters in benign and malignant gliomas. We explore these elements and intersect with Hi-C data to uncover molecular mechanisms instructing gene expression in gliomas.

Unveiling new interdependencies between significant DNA methylation sites, gene expression profiles and glioma patients survival.

In order to find clinically useful prognostic markers for glioma patients' survival, we employed Monte Carlo Feature Selection and Interdependencies Discovery (MCFS-ID) algorithm on DNA methylation (HumanMethylation450 platform) and RNA-seq datasets from The Cancer Genome Atlas (TCGA) for 88 patients observed until death. The input features were ranked according to their importance in predicting patients' longer (400+ days) or shorter (≤400 days) survival without prior classification of the patients. Interestingly, out of the 65 most important features found, 63 are methylation sites, and only two mRNAs. Moreover, 61 out of the 63 methylation sites are among those detected by the 450 k array technology, while being absent in the HumanMethylation27. The most important methylation feature (cg15072976) overlaps with the RE1 Silencing Transcription Factor (REST) binding site, and was confirmed to intersect with the REST binding motif in human U87 glioma cells. Six additional methylation sites from the top 63 overlap with REST sites. We found that the methylation status of the cg15072976 site affects transcription factor binding in U87 cells in gel shift assay. The cg15072976 methylation status discriminates ≤400 and 400+ patients in an independent dataset from TCGA and shows positive association with survival time as evidenced by Kaplan-Meier plots.

Computational Analysis of Molecular Interaction Networks Underlying Change of HIV-1 Resistance to Selected Reverse Transcriptase Inhibitors.

MOTIVATION: Despite more than two decades of research, HIV resistance to drugs remains a serious obstacle in developing efficient AIDS treatments. Several computational methods have been developed to predict resistance level from the sequence of viral proteins such as reverse transcriptase (RT) or protease. These methods, while powerful and accurate, give very little insight into the molecular interactions that underly acquisition of drug resistance/hypersusceptibility. Here, we attempt at filling this gap by using our Monte Carlo feature selection and interdependency discovery method (MCFS-ID) to elucidate molecular interaction networks that characterize viral strains with altered drug resistance levels. RESULTS: We analyzed a number of HIV-1 RT sequences annotated with drug resistance level using the MCFS-ID method. This let us expound interdependency networks that characterize change of drug resistance to six selected RT inhibitors: Abacavir, Lamivudine, Stavudine, Zidovudine, Tenofovir and Nevirapine. The networks consider interdependencies at the level of physicochemical properties of mutating amino acids, eg,: polarity. We mapped each network on the 3D structure of RT in attempt to understand the molecular meaning of interacting pairs. The discovered interactions describe several known drug resistance mechanisms and, importantly, some previously unidentified ones. Our approach can be easily applied to a whole range of problems from the domain of protein engineering. AVAILABILITY: A portable Java implementation of our MCFS-ID method is freely available for academic users and can be obtained at: http://www.ipipan.eu/staff/m.draminski/software.htm.

A Rough Set-Based Model of HIV-1 Reverse Transcriptase Resistome.

Reverse transcriptase (RT) is a viral enzyme crucial for HIV-1 replication. Currently, 12 drugs are targeted against the RT. The low fidelity of the RT-mediated transcription leads to the quick accumulation of drug-resistance mutations. The sequence-resistance relationship remains only partially understood. Using publicly available data collected from over 15 years of HIV proteome research, we have created a general and predictive rule-based model of HIV-1 resistance to eight RT inhibitors. Our rough set-based model considers changes in the physicochemical properties of a mutated sequence as compared to the wild-type strain. Thanks to the application of the Monte Carlo feature selection method, the model takes into account only the properties that significantly contribute to the resistance phenomenon. The obtained results show that drug-resistance is determined in more complex way than believed. We confirmed the importance of many resistance-associated sites, found some sites to be less relevant than formerly postulated and-more importantly-identified several previously neglected sites as potentially relevant. By mapping some of the newly discovered sites on the 3D structure of the RT, we were able to suggest possible molecular-mechanisms of drug-resistance. Importantly, our model has the ability to generalize predictions to the previously unseen cases. The study is an example of how computational biology methods can increase our understanding of the HIV-1 resistome.

Monte Carlo feature selection for supervised classification.

MOTIVATION: Pre-selection of informative features for supervised classification is a crucial, albeit delicate, task. It is desirable that feature selection provides the features that contribute most to the classification task per se and which should therefore be used by any classifier later used to produce classification rules. In this article, a conceptually simple but computer-intensive approach to this task is proposed. The reliability of the approach rests on multiple construction of a tree classifier for many training sets randomly chosen from the original sample set, where samples in each training set consist of only a fraction of all of the observed features. RESULTS: The resulting ranking of features may then be used to advantage for classification via a classifier of any type. The approach was validated using Golub et al. leukemia data and the Alizadeh et al. lymphoma data. Not surprisingly, we obtained a significantly different list of genes. Biological interpretation of the genes selected by our method showed that several of them are involved in precursors to different types of leukemia and lymphoma rather than being genes that are common to several forms of cancers, which is the case for the other methods. AVAILABILITY: Prototype available upon request.