Systematic integration of m6A regulators and autophagy-related genes in combination with long non-coding RNAs predicts survival in glioblastoma multiforme.

Glioblastoma multiforme (GBM) is probably the only tumor in which a unique epigenetic alteration, namely methylation of the MGMT gene, possesses direct clinical relevance. Now with the emergence of aberrant N6 methyladenosine (m6A) modifications (the most common epigenetic modification of mRNA, closely linked to the autophagy process) in cancer, the epi-transcriptomic landscape of GBM pathobiology has been expanded. Considering this, herein, we systematically analyzed m6A regulators, assessed their correlation with autophagy-related genes (ATG), and established a long non-coding RNAs (lncRNA)-dependent prognostic signature (m6A-autophagy-lncRNAs) for GBM. Our analysis identified a novel signature of five long non-coding RNAs (lncRNAs: ITGA6-AS1, AC124248.1, NFYC-AS1, AC025171.1, and AC005229.3) associated with survival of GBM patients, and four among them clearly showed cancer-associated potential. We further validated and confirmed the altered expression of two lncRNAs (AC124248.1, AC005229.3) in GBM associated clinical samples using RT-PCR. Concerning the prognostic ability, the obtained signature determined high-/low-risk groups in GBM patients and showed sensitivity to anticancer drugs. Collectively, the m6A-autophagy-lncRNAs signature presented in the study is clinically relevant and is the first attempt to systematically predict the potential interaction between the three key determinants (m6A, autophagy, lncRNA) in cancer, particularly in GBM.

ZBTB18 inhibits SREBP-dependent lipid synthesis by halting CTBPs and LSD1 activity in glioblastoma.

Enhanced fatty acid synthesis is a hallmark of tumors, including glioblastoma. SREBF1/2 regulate the expression of enzymes involved in fatty acid and cholesterol synthesis. Yet, little is known about the precise mechanism regulating SREBP gene expression in glioblastoma. Here, we show that a novel interaction between the co-activator/co-repressor CTBP and the tumor suppressor ZBTB18 regulates the expression of SREBP genes. In line with our findings, metabolic assays and glucose tracing analysis confirm the reduction in several phospholipid species upon ZBTB18 expression. Our study identifies CTBP1/2 and LSD1 as co-activators of SREBP genes and indicates that the functional activity of the CTBP-LSD1 complex is altered by ZBTB18. ZBTB18 binding to the SREBP gene promoters is associated with reduced LSD1 demethylase activity of H3K4me2 and H3K9me2 marks. Concomitantly, the interaction between LSD1, CTBP, and ZNF217 is increased, suggesting that ZBTB18 promotes LSD1 scaffolding function. Our results outline a new epigenetic mechanism enrolled by ZBTB18 and its co-factors to regulate fatty acid synthesis that could be targeted to treat glioblastoma patients.

T-cell dysfunction in the glioblastoma microenvironment is mediated by myeloid cells releasing interleukin-10.

Despite recent advances in cancer immunotherapy, certain tumor types, such as Glioblastomas, are highly resistant due to their tumor microenvironment disabling the anti-tumor immune response. Here we show, by applying an in-silico multidimensional model integrating spatially resolved and single-cell gene expression data of 45,615 immune cells from 12 tumor samples, that a subset of Interleukin-10-releasing HMOX1+ myeloid cells, spatially localizing to mesenchymal-like tumor regions, drive T-cell exhaustion and thus contribute to the immunosuppressive tumor microenvironment. These findings are validated using a human ex-vivo neocortical glioblastoma model inoculated with patient derived peripheral T-cells to simulate the immune compartment. This model recapitulates the dysfunctional transformation of tumor infiltrating T-cells. Inhibition of the JAK/STAT pathway rescues T-cell functionality both in our model and in-vivo, providing further evidence of IL-10 release being an important driving force of tumor immune escape. Our results thus show that integrative modelling of single cell and spatial transcriptomics data is a valuable tool to interrogate the tumor immune microenvironment and might contribute to the development of successful immunotherapies.

NF1 regulates mesenchymal glioblastoma plasticity and aggressiveness through the AP-1 transcription factor FOSL1.

The molecular basis underlying glioblastoma (GBM) heterogeneity and plasticity is not fully understood. Using transcriptomic data of human patient-derived brain tumor stem cell lines (BTSCs), classified based on GBM-intrinsic signatures, we identify the AP-1 transcription factor FOSL1 as a key regulator of the mesenchymal (MES) subtype. We provide a mechanistic basis to the role of the neurofibromatosis type 1 gene (NF1), a negative regulator of the RAS/MAPK pathway, in GBM mesenchymal transformation through the modulation of FOSL1 expression. Depletion of FOSL1 in NF1-mutant human BTSCs and Kras-mutant mouse neural stem cells results in loss of the mesenchymal gene signature and reduction in stem cell properties and in vivo tumorigenic potential. Our data demonstrate that FOSL1 controls GBM plasticity and aggressiveness in response to NF1 alterations.

The aberrant splicing of BAF45d links splicing regulation and transcription in glioblastoma.

Background: Glioblastoma, the most aggressive primary brain tumor, is genetically heterogeneous. Alternative splicing (AS) plays a key role in numerous pathologies, including cancer. The objectives of our study were to determine whether aberrant AS could play a role in the malignant phenotype of glioma and to understand the mechanism underlying its aberrant regulation. Methods: We obtained surgical samples from patients with glioblastoma who underwent 5-aminolevulinic fluorescence-guided surgery. Biopsies were taken from the tumor center as well as from adjacent normal-appearing tissue. We used a global splicing array to identify candidate genes aberrantly spliced in these glioblastoma samples. Mechanistic and functional studies were performed to elucidate the role of our top candidate splice variant, BAF45d, in glioblastoma. Results: BAF45d is part of the switch/sucrose nonfermentable complex and plays a key role in the development of the CNS. The BAF45d/6A isoform is present in 85% of over 200 glioma samples that have been analyzed and contributes to the malignant glioma phenotype through the maintenance of an undifferentiated cellular state. We demonstrate that BAF45d splicing is mediated by polypyrimidine tract-binding protein 1 (PTBP1) and that BAF45d regulates PTBP1, uncovering a reciprocal interplay between RNA splicing regulation and transcription. Conclusions: Our data indicate that AS is a mechanism that contributes to the malignant phenotype of glioblastoma. Understanding the consequences of this biological process will uncover new therapeutic targets for this devastating disease.

Microenvironment-Derived Regulation of HIF Signaling Drives Transcriptional Heterogeneity in Glioblastoma Multiforme.

The evolving and highly heterogeneous nature of malignant brain tumors underlies their limited response to therapy and poor prognosis. In addition to genetic alterations, highly dynamic processes, such as transcriptional and metabolic reprogramming, play an important role in the development of tumor heterogeneity. The current study reports an adaptive mechanism in which the metabolic environment of malignant glioma drives transcriptional reprogramming. Multiregional analysis of a glioblastoma patient biopsy revealed a metabolic landscape marked by varying stages of hypoxia and creatine enrichment. Creatine treatment and metabolism was further shown to promote a synergistic effect through upregulation of the glycine cleavage system and chemical regulation of prolyl-hydroxylase domain. Consequently, creatine maintained a reduction of reactive oxygen species and change of the α-ketoglutarate/succinate ratio, leading to an inhibition of HIF signaling in primary tumor cell lines. These effects shifted the transcriptional pattern toward a proneural subtype and reduced the rate of cell migration and invasion in vitroImplications: Transcriptional subclasses of glioblastoma multiforme are heterogeneously distributed within the same tumor. This study uncovered a regulatory function of the tumor microenvironment by metabolism-driven transcriptional reprogramming in infiltrating glioma cells. Mol Cancer Res; 16(4); 655-68. ©2018 AACR.

EGFR feedback-inhibition by Ran-binding protein 6 is disrupted in cancer.

Transport of macromolecules through the nuclear pore by importins and exportins plays a critical role in the spatial regulation of protein activity. How cancer cells co-opt this process to promote tumorigenesis remains unclear. The epidermal growth factor receptor (EGFR) plays a critical role in normal development and in human cancer. Here we describe a mechanism of EGFR regulation through the importin ß family member RAN-binding protein 6 (RanBP6), a protein of hitherto unknown functions. We show that RanBP6 silencing impairs nuclear translocation of signal transducer and activator of transcription 3 (STAT3), reduces STAT3 binding to the EGFR promoter, results in transcriptional derepression of EGFR, and increased EGFR pathway output. Focal deletions of the RanBP6 locus on chromosome 9p were found in a subset of glioblastoma (GBM) and silencing of RanBP6 promoted glioma growth in vivo. Our results provide an example of EGFR deregulation in cancer through silencing of components of the nuclear import pathway.

Expression differences of programmed death ligand 1 in de-novo and recurrent glioblastoma multiforme.

The biology of recurrent glioblastoma multiforme (GBM) is a dynamic process influenced by selection pressure induced by different antitumoural therapies. The poor clinical outcome of tumours in the recurrent stage necessitates the development of effective therapeutic strategies. Checkpoint-inhibition (PD1/PD-L1 Inhibition) is a hallmark of immunotherapy being investigated in ongoing clinical trials. The purpose of this study was to analyse the PD-L1 expression in de-novo and recurrent glioblastoma multiforme and to explore associated genetic alterations and clinical traits. We show that PD-L1 expression was reduced in recurrent GBM in comparison to de-novo GBM. Additionally, patients who received an extended dose of temozolomide (TMZ) chemotherapy showed a significantly reduced level of PD-L1 expression in the recurrence stage compared to the corresponding de-novo tumour. Our findings may provide an explanation for potentially lower response to immunotherapy in the recurrent stage due to the reduced expression of the therapeutic target PD-L1.

Comprehensive analysis of PD-L1 expression in glioblastoma multiforme.

Glioblastoma multiforme are highly malignant brain tumours with frequent genetic and epigenetic alterations. The poor clinical outcome of these tumours necessitates the development of new treatment options. Immunotherapies for glioblastoma multiforme including PD1/PD-L1 inhibition are currently tested in ongoing clinical trials. The purpose of this study was to investigate the molecular background of PD-L1 expression in glioblastoma multiforme and to find associated pathway activation and genetic alterations. We show that PD-L1 is up-regulated in IDH1/2 wildtype glioblastoma multiforme compared to lower-grade gliomas. In addition, a strong association of PD-L1 with the mesenchymal expression subgroup was observed. Consistent with that, NF1 mutation and corresponding activation of the MAPK pathway was strongly connected to PD-L1 expression. Our findings may explain different response to PD-L1 inhibition of patients in ongoing trials and may help to select patients that may profit of immunotherapy in the future.

Integrative Modeling Reveals Annexin A2-mediated Epigenetic Control of Mesenchymal Glioblastoma.

Glioblastomas are characterized by transcriptionally distinct subtypes, but despite possible clinical relevance, their regulation remains poorly understood. The commonly used molecular classification systems for GBM all identify a subtype with high expression of mesenchymal marker transcripts, strongly associated with invasive growth. We used a comprehensive data-driven network modeling technique (augmented sparse inverse covariance selection, aSICS) to define separate genomic, epigenetic, and transcriptional regulators of glioblastoma subtypes. Our model identified Annexin A2 (ANXA2) as a novel methylation-controlled positive regulator of the mesenchymal subtype. Subsequent evaluation in two independent cohorts established ANXA2 expression as a prognostic factor that is dependent on ANXA2 promoter methylation. ANXA2 knockdown in primary glioblastoma stem cell-like cultures suppressed known mesenchymal master regulators, and abrogated cell proliferation and invasion. Our results place ANXA2 at the apex of a regulatory cascade that determines glioblastoma mesenchymal transformation and validate aSICS as a general methodology to uncover regulators of cancer subtypes.

Integrative Network-based Analysis of Magnetic Resonance Spectroscopy and Genome Wide Expression in Glioblastoma multiforme.

The goal of this study was to identify correlations between metabolites from proton MR spectroscopy and genetic pathway activity in glioblastoma multiforme (GBM). Twenty patients with primary GBM were analysed by short echo-time chemical shift imaging and genome-wide expression analyses. Weighed Gene Co-Expression Analysis was used for an integrative analysis of imaging and genetic data. N-acetylaspartate, normalised to the contralateral healthy side (nNAA), was significantly correlated to oligodendrocytic and neural development. For normalised creatine (nCr), a group with low nCr was linked to the mesenchymal subtype, while high nCr could be assigned to the proneural subtype. Moreover, clustering of normalised glutamine and glutamate (nGlx) revealed two groups, one with high nGlx being attributed to the neural subtype, and one with low nGlx associated with the classical subtype. Hence, the metabolites nNAA, nCr, and nGlx correlate with a specific gene expression pattern reflecting the previously described subtypes of GBM. Moreover high nNAA was associated with better clinical prognosis, whereas patients with lower nNAA revealed a shorter progression-free survival (PFS).

Functional role of CLIC1 ion channel in glioblastoma-derived stem/progenitor cells.

BACKGROUND: Chloride channels are physiologically involved in cell division and motility. Chloride intracellular channel 1 (CLIC1) is overexpressed in a variety of human solid tumors compared with normal tissues, suggesting a potential involvement of CLIC1 in the regulation of tumorigenesis. This led us to investigate the role of CLIC1 in gliomagenesis. METHODS: We used the neurosphere system to isolate stem/progenitor cells from human glioblastomas (GBMs). CLIC1 targeting in GBM neurospheres was achieved by both lentiviral-mediated short-hairpin RNA transduction and CLIC1 antibody treatment, and its effect on stem-like properties was analyzed in vitro by proliferation and clonogenic assays and in vivo by orthotopic injection in immunocompromised mice. Channel activity was studied by perforated patch clamp technique. Differences in expression were analyzed by analysis of variance with Tamhane's multiple comparison test. Kaplan-Meier analyses and log-rank test were used to assess survival. All statistical tests were two-sided. RESULTS: CLIC1 was statistically significantly overexpressed in GBMs compared with normal brain tissues (P < .001) with a better survival of patients with CLIC1 low-expressing tumors (CLIC1(low) vs CLIC1(high) survival: χ(2) = 74.35; degrees of freedom = 1; log-rank P < .001). CLIC1 was variably expressed in patient-derived GBM neurospheres and was found enriched in the stem/progenitor compartment. CLIC1 silencing reduced proliferative (P < .01), clonogenic (P < .01), and tumorigenic capacity (P < .05) of stem/progenitor cells. The reduction of CLIC1 chloride currents with a specific CLIC1 antibody mirrored the biological effects of CLIC1 silencing in GBM patient-derived neurospheres. CONCLUSIONS: Reduced gliomagenesis after CLIC1 targeting in tumoral stem/progenitor cells and the finding that CLIC1 expression is inversely associated with patient survival suggest CLIC1 as a potential target and prognostic biomarker.

Resistance to hypoxia-induced, BNIP3-mediated cell death contributes to an increase in a CD133-positive cell population in human glioblastomas in vitro.

In addition to intrinsic regulatory mechanisms, brain tumor stemlike cells (BTSCs), a small subpopulation of malignant glial tumor-derived cells, are influenced by environmental factors. Previous reports showed that lowering oxygen tension induced an increase of BTSCs expressing CD133 and other stem cell-related genes and more pronounced clonogenic capacity in vitro. We investigated the mechanisms responsible for hypoxia-dependent induction of CD133-positive BTSCs in glioblastomas. We confirmed that cultures exposed to lowered oxygen levels showed a severalfold increase of CD133-positive BTSCs. Both the increase of CD133-positive cells and deceleration of the growth kinetics were reversible after transfer to normoxic conditions. Exposure to hypoxia induced BNIP3 (BCL2/adenovirus E1B 19-kDa protein-interacting protein 3)-dependent apoptosis preferentially in CD133-negative cells. In contrast, CD133-positive cells proved to be more resistant to hypoxia-induced programmed cell death. Application of the demethylating agent 5'-azacitidine resulted in an increase of BNIP3 expression levels in CD133-positive cells. Thus, epigenetic modifications led to their better survival in lowered oxygen tension. Moreover, the, hypoxia-induced increase of CD133-positive cells was inhibited after 5'-azacitidine treatment. These results suggest the possible efficacy of a novel therapy for glioblastoma focused on eradication of BTSCs by modifications of epigenetic regulation of gene expression.

The transcriptional network for mesenchymal transformation of brain tumours.

The inference of transcriptional networks that regulate transitions into physiological or pathological cellular states remains a central challenge in systems biology. A mesenchymal phenotype is the hallmark of tumour aggressiveness in human malignant glioma, but the regulatory programs responsible for implementing the associated molecular signature are largely unknown. Here we show that reverse-engineering and an unbiased interrogation of a glioma-specific regulatory network reveal the transcriptional module that activates expression of mesenchymal genes in malignant glioma. Two transcription factors (C/EBPbeta and STAT3) emerge as synergistic initiators and master regulators of mesenchymal transformation. Ectopic co-expression of C/EBPbeta and STAT3 reprograms neural stem cells along the aberrant mesenchymal lineage, whereas elimination of the two factors in glioma cells leads to collapse of the mesenchymal signature and reduces tumour aggressiveness. In human glioma, expression of C/EBPbeta and STAT3 correlates with mesenchymal differentiation and predicts poor clinical outcome. These results show that the activation of a small regulatory module is necessary and sufficient to initiate and maintain an aberrant phenotypic state in cancer cells.

The N-Myc-DLL3 cascade is suppressed by the ubiquitin ligase Huwe1 to inhibit proliferation and promote neurogenesis in the developing brain.

Self-renewal and proliferation of neural stem cells and the decision to initiate neurogenesis are crucial events directing brain development. Here we show that the ubiquitin ligase Huwe1 operates upstream of the N-Myc-DLL3-Notch pathway to control neural stem cell activity and promote neurogenesis. Conditional inactivation of the Huwe1 gene in the mouse brain caused neonatal lethality associated with disorganization of the laminar patterning of the cortex. These defects stemmed from severe impairment of neurogenesis associated with uncontrolled expansion of the neural stem cell compartment. Loss- and gain-of-function experiments in the mouse cortex demonstrated that Huwe1 restrains proliferation and enables neuronal differentiation by suppressing the N-Myc-DLL3 cascade. Notably, human high-grade gliomas carry focal hemizygous deletions of the X-linked Huwe1 gene in association with amplification of the N-myc locus. Our results indicate that Huwe1 balances proliferation and neurogenesis in the developing brain and that this pathway is subverted in malignant brain tumors.

Localizing hotspots of antisense transcription.

Analysis of the transcriptome by computational and experimental methods has established that sense-antisense transcriptional units are a common phenomenon. Although the regulatory potential of antisense transcripts has been experimentally verified in a number of studies, the biological importance of sense-antisense regulation of gene expression is still a matter of debate. Here, we report the identification of sequence features that are associated with antisense transcription. We show that the sequence composition of the first exon and the 5'end of the first intron of many human genes is similar to the sequence composition observed in promoter regions as measured by the density of known transcription regulatory motifs. Cloned intron-derived fragments were found to possess bidirectional promoter activity. In agreement with the reported abundance of antisense transcripts overlapping the 5'UTR, mapping of the 5'ends of antisense transcripts to the corresponding sense transcripts revealed that the first exon and the 5'end of the first intron are hotspots of antisense transcription as measured by the number of antisense transcription start sites per unit sequence. CpG dinucleotide suppression that is typically weak in non-methylated promoter regions is similarly weakened upstream as well as downstream of the first exon. In support of antisense transcripts playing a regulatory role, we find that 5'UTRs and first exons of genes with overlapping antisense transcripts are significantly longer than the genomic average. Interestingly, a similar size distribution of 5'UTRs and first exons is observed for genes silenced by CpG island methylation in human cancer.

DEK Expression is controlled by E2F and deregulated in diverse tumor types.

Deregulation of the retinoblastoma (pRB) tumor suppressor pathway associated with aberrant activity of E2F transcription factors is frequently observed in human cancer. Microarray based analyses have revealed a large number of potential downstream mediators of the tumor suppressing activity of pRB, including DEK, a fusion partner of CAN found in a subset of acute myeloid leukaemia (AML) patients carrying a (6; 9) translocation. Here we report that the expression of DEK is under direct control of E2F transcription factors. Chromatin immunoprecipitation assays show that the DEK promoter is bound by endogenous E2F in vivo. The DEK promoter is transactivated by E2F and mutation of E2F binding sites eliminates this effect. Expression levels of DEK in human tumors have been investigated by tissue micro array analysis. We find that DEK is overexpressed in many solid tumors such as colon cancer, larynx cancer, bladder cancer, and melanoma.