Glioblastoma Multiforme (GBM), the most common and lethal primary human brain tumor, exhibits multiple molecular aberrations. We report that loss of the transcription factor GATA4, a negative regulator of normal astrocyte proliferation, is a driver in glioma formation and fulfills the hallmarks of a tumor suppressor gene (TSG). Although GATA4 was expressed in normal brain, loss of GATA4 was observed in 94/163 GBM operative samples and was a negative survival prognostic marker. GATA4 loss occurred through promoter hypermethylation or novel somatic mutations. Loss of GATA4 in normal human astrocytes promoted high-grade astrocytoma formation, in cooperation with other relevant genetic alterations such as activated Ras or loss of TP53. Loss of GATA4 with activated Ras in normal astrocytes promoted a progenitor-like phenotype, formation of neurospheres, and the ability to differentiate into astrocytes, neurons, and oligodendrocytes. Re-expression of GATA4 in human GBM cell lines, primary cultures, and brain tumor-initiating cells suppressed tumor growth in vitro and in vivo through direct activation of the cell cycle inhibitor P21(CIP1), independent of TP53. Re-expression of GATA4 also conferred sensitivity of GBM cells to temozolomide, a DNA alkylating agent currently used in GBM therapy. This sensitivity was independent of MGMT (O-6-methylguanine-DNA-methyltransferase), the DNA repair enzyme which is often implicated in temozolomide resistance. Instead, GATA4 reduced expression of APNG (alkylpurine-DNA-N-glycosylase), a DNA repair enzyme which is poorly characterized in GBM-mediated temozolomide resistance. Identification and validation of GATA4 as a TSG and its downstream targets in GBM may yield promising novel therapeutic strategies.

Telomerase reverse transcriptase (TERT) promoter mutations were recently shown to drive telomerase activity in various cancer types, including medulloblastoma. However, the clinical and biological implications of TERT mutations in medulloblastoma have not been described. Hence, we sought to describe these mutations and their impact in a subgroup-specific manner. We analyzed the TERT promoter by direct sequencing and genotyping in 466 medulloblastomas. The mutational distributions were determined according to subgroup affiliation, demographics, and clinical, prognostic, and molecular features. Integrated genomics approaches were used to identify specific somatic copy number alterations in TERT promoter-mutated and wild-type tumors. Overall, TERT promoter mutations were identified in 21 % of medulloblastomas. Strikingly, the highest frequencies of TERT mutations were observed in SHH (83 %; 55/66) and WNT (31 %; 4/13) medulloblastomas derived from adult patients. Group 3 and Group 4 harbored this alteration in <5 % of cases and showed no association with increased patient age. The prognostic implications of these mutations were highly subgroup-specific. TERT mutations identified a subset with good and poor prognosis in SHH and Group 4 tumors, respectively. Monosomy 6 was mostly restricted to WNT tumors without TERT mutations. Hallmark SHH focal copy number aberrations and chromosome 10q deletion were mutually exclusive with TERT mutations within SHH tumors. TERT promoter mutations are the most common recurrent somatic point mutation in medulloblastoma, and are very highly enriched in adult SHH and WNT tumors. TERT mutations define a subset of SHH medulloblastoma with distinct demographics, cytogenetics, and outcomes.

Amplification of the C19MC oncogenic miRNA cluster and high LIN28 expression has been linked to a distinctly aggressive group of cerebral CNS-PNETs (group 1 CNS-PNETs) arising in young children. In this study, we sought to evaluate the diagnostic specificity of C19MC and LIN28, and the clinical and biological spectra of C19MC amplified and/or LIN28+ CNS-PNETs. We interrogated 450 pediatric brain tumors using FISH and IHC analyses and demonstrate that C19MC alteration is restricted to a sub-group of CNS-PNETs with high LIN28 expression; however, LIN28 immunopositivity was not exclusive to CNS-PNETs but was also detected in a proportion of other malignant pediatric brain tumors including rhabdoid brain tumors and malignant gliomas. C19MC amplified/LIN28+ group 1 CNS-PNETs arose predominantly in children <4 years old; a majority arose in the cerebrum but 24 % (13/54) of tumors had extra-cerebral origins. Notably, group 1 CNS-PNETs encompassed several histologic classes including embryonal tumor with abundant neuropil and true rosettes (ETANTR), medulloepithelioma, ependymoblastoma and CNS-PNETs with variable differentiation. Strikingly, gene expression and methylation profiling analyses revealed a common molecular signature enriched for primitive neural features, high LIN28/LIN28B and DNMT3B expression for all group 1 CNS-PNETs regardless of location or tumor histology. Our collective findings suggest that current known histologic categories of CNS-PNETs which include ETANTRs, medulloepitheliomas, ependymoblastomas in various CNS locations, comprise a common molecular and diagnostic entity and identify inhibitors of the LIN28/let7/PI3K/mTOR axis and DNMT3B as promising therapeutics for this distinct histogenetic entity.

While molecular subgrouping has revolutionized medulloblastoma classification, the extent of heterogeneity within subgroups is unknown. Similarity network fusion (SNF) applied to genome-wide DNA methylation and gene expression data across 763 primary samples identifies very homogeneous clusters of patients, supporting the presence of medulloblastoma subtypes. After integration of somatic copy-number alterations, and clinical features specific to each cluster, we identify 12 different subtypes of medulloblastoma. Integrative analysis using SNF further delineates group 3 from group 4 medulloblastoma, which is not as readily apparent through analyses of individual data types. Two clear subtypes of infants with Sonic Hedgehog medulloblastoma with disparate outcomes and biology are identified. Medulloblastoma subtypes identified through integrative clustering have important implications for stratification of future clinical trials.

Glioblastoma therapies have remained elusive due to limitations in understanding mechanisms of growth and survival of the tumorigenic population. Using CRISPR-Cas9 approaches in patient-derived GBM stem cells (GSCs) to interrogate function of the coding genome, we identify actionable pathways responsible for growth, which reveal the gene-essential circuitry of GBM stemness and proliferation. In particular, we characterize members of the SOX transcription factor family, SOCS3, USP8, and DOT1L, and protein ufmylation as important for GSC growth. Additionally, we reveal mechanisms of temozolomide resistance that could lead to combination strategies. By reaching beyond static genome analysis of bulk tumors, with a genome-wide functional approach, we reveal genetic dependencies within a broad range of biological processes to provide increased understanding of GBM growth and treatment resistance.

We investigated the role of 3D genome architecture in instructing functional properties of glioblastoma stem cells (GSCs) by generating sub-5-kb resolution 3D genome maps by in situ Hi-C. Contact maps at sub-5-kb resolution allow identification of individual DNA loops, domain organization, and large-scale genome compartmentalization. We observed differences in looping architectures among GSCs from different patients, suggesting that 3D genome architecture is a further layer of inter-patient heterogeneity for glioblastoma. Integration of DNA contact maps with chromatin and transcriptional profiles identified specific mechanisms of gene regulation, including the convergence of multiple super enhancers to individual stemness genes within individual cells. We show that the number of loops contacting a gene correlates with elevated transcription. These results indicate that stemness genes are hubs of interaction between multiple regulatory regions, likely to ensure their sustained expression. Regions of open chromatin common among the GSCs tested were poised for expression of immune-related genes, including CD276 We demonstrate that this gene is co-expressed with stemness genes in GSCs and that CD276 can be targeted with an antibody-drug conjugate to eliminate self-renewing cells. Our results demonstrate that integrated structural genomics data sets can be employed to rationally identify therapeutic vulnerabilities in self-renewing cells.

Childhood brain tumors have suspected prenatal origins. To identify vulnerable developmental states, we generated a single-cell transcriptome atlas of >65,000 cells from embryonal pons and forebrain, two major tumor locations. We derived signatures for 191 distinct cell populations and defined the regional cellular diversity and differentiation dynamics. Projection of bulk tumor transcriptomes onto this dataset shows that WNT medulloblastomas match the rhombic lip-derived mossy fiber neuronal lineage and embryonal tumors with multilayered rosettes fully recapitulate a neuronal lineage, while group 2a/b atypical teratoid/rhabdoid tumors may originate outside the neuroectoderm. Importantly, single-cell tumor profiles reveal highly defined cell hierarchies that mirror transcriptional programs of the corresponding normal lineages. Our findings identify impaired differentiation of specific neural progenitors as a common mechanism underlying these pediatric cancers and provide a rational framework for future modeling and therapeutic interventions.

High-grade gliomas defined by histone 3 K27M driver mutations exhibit global loss of H3K27 trimethylation and reciprocal gain of H3K27 acetylation, respectively shaping repressive and active chromatin landscapes. We generated tumor-derived isogenic models bearing this mutation and show that it leads to pervasive H3K27ac deposition across the genome. In turn, active enhancers and promoters are not created de novo and instead reflect the epigenomic landscape of the cell of origin. H3K27ac is enriched at repeat elements, resulting in their increased expression, which in turn can be further amplified by DNA demethylation and histone deacetylase inhibitors providing an exquisite therapeutic vulnerability. These agents may therefore modulate anti-tumor immune responses as a therapeutic modality for this untreatable disease.

Functional heterogeneity within tumors presents a significant therapeutic challenge. Here we show that quiescent, therapy-resistant Sox2(+) cells propagate sonic hedgehog subgroup medulloblastoma by a mechanism that mirrors a neurogenic program. Rare Sox2(+) cells produce rapidly cycling doublecortin(+) progenitors that, together with their postmitotic progeny expressing NeuN, comprise tumor bulk. Sox2(+) cells are enriched following anti-mitotic chemotherapy and Smoothened inhibition, creating a reservoir for tumor regrowth. Lineage traces from Sox2(+) cells increase following treatment, suggesting that this population is responsible for relapse. Targeting Sox2(+) cells with the antineoplastic mithramycin abrogated tumor growth. Addressing functional heterogeneity and eliminating Sox2(+) cells presents a promising therapeutic paradigm for treatment of sonic hedgehog subgroup medulloblastoma.

Developmental signal transduction pathways act diversely, with context-dependent roles across systems and disease types. Glioblastomas (GBMs), which are the poorest prognosis primary brain cancers, strongly resemble developmental systems, but these growth processes have not been exploited therapeutically, likely in part due to the extreme cellular and genetic heterogeneity observed in these tumors. The role of Wnt/βcatenin signaling in GBM stem cell (GSC) renewal and fate decisions remains controversial. Here, we report context-specific actions of Wnt/βcatenin signaling in directing cellular fate specification and renewal. A subset of primary GBM-derived stem cells requires Wnt proteins for self-renewal, and this subset specifically relies on Wnt/βcatenin signaling for enhanced tumor burden in xenograft models. In an orthotopic Wnt reporter model, Wnthi GBM cells (which exhibit high levels of βcatenin signaling) are a faster-cycling, highly self-renewing stem cell pool. In contrast, Wntlo cells (with low levels of signaling) are slower cycling and have decreased self-renewing potential. Dual inhibition of Wnt/βcatenin and Notch signaling in GSCs that express high levels of the proneural transcription factor ASCL1 leads to robust neuronal differentiation and inhibits clonogenic potential. Our work identifies new contexts for Wnt modulation for targeting stem cell differentiation and self-renewal in GBM heterogeneity, which deserve further exploration therapeutically.

Alteration of tissue mechanical properties is a physical hallmark of solid tumors including gliomas. How tumor cells sense and regulate tissue mechanics is largely unknown. Here, we show that mechanosensitive ion channel Piezo regulates mitosis and tissue stiffness of Drosophila gliomas, but not non-transformed brains. PIEZO1 is overexpressed in aggressive human gliomas and its expression inversely correlates with patient survival. Deleting PIEZO1 suppresses the growth of glioblastoma stem cells, inhibits tumor development, and prolongs mouse survival. Focal mechanical force activates prominent PIEZO1-dependent currents from glioma cell processes, but not soma. PIEZO1 localizes at focal adhesions to activate integrin-FAK signaling, regulate extracellular matrix, and reinforce tissue stiffening. In turn, a stiffer mechanical microenvironment elevates PIEZO1 expression to promote glioma aggression. Therefore, glioma cells are mechanosensory in a PIEZO1-dependent manner, and targeting PIEZO1 represents a strategy to break the reciprocal, disease-aggravating feedforward circuit between tumor cell mechanotransduction and the aberrant tissue mechanics. VIDEO ABSTRACT.

Transcription factors (TFs) are spatially and temporally regulated during gut organ specification. Although accumulating evidence shows aberrant reactivation of developmental programs in cancer, little is known about how TFs drive lineage specification in development and cancer. We first defined gastrointestinal tissue-specific chromatin accessibility and gene expression during development, identifying the dynamic epigenetic regulation of SOX family of TFs. We revealed that Sox2 is not only essential for gastric specification, by maintaining chromatin accessibility at forestomach lineage loci, but also sufficient to promote forestomach/esophageal transformation upon Cdx2 deletion. By comparing our gastrointestinal lineage-specific transcriptome to human gastrointestinal cancer data, we found that stomach and intestinal lineage-specific programs are reactivated in Sox2high /Sox9high and Cdx2high cancers, respectively. By analyzing mice deleted for both Sox2 and Sox9, we revealed their potentially redundant roles in both gastric development and cancer, highlighting the importance of developmental lineage programs reactivated by gastrointestinal TFs in cancer.

Glioblastomas exhibit a hierarchical cellular organization, suggesting that they are driven by neoplastic stem cells that retain partial yet abnormal differentiation potential. Here, we show that a large subset of patient-derived glioblastoma stem cells (GSCs) express high levels of Achaete-scute homolog 1 (ASCL1), a proneural transcription factor involved in normal neurogenesis. ASCL1hi GSCs exhibit a latent capacity for terminal neuronal differentiation in response to inhibition of Notch signaling, whereas ASCL1lo GSCs do not. Increasing ASCL1 levels in ASCL1lo GSCs restores neuronal lineage potential, promotes terminal differentiation, and attenuates tumorigenicity. ASCL1 mediates these effects by functioning as a pioneer factor at closed chromatin, opening new sites to activate a neurogenic gene expression program. Directing GSCs toward terminal differentiation may provide therapeutic applications for a subset of GBM patients and strongly supports efforts to restore differentiation potential in GBM and other cancers.

Long non-coding RNAs (lncRNAs) are increasingly recognized as functional units in cancer and powerful biomarkers; however, most remain uncharacterized. Here, we analyze 5,592 prognostic lncRNAs in 9,446 cancers of 30 types using machine learning. We identify 166 lncRNAs whose expression correlates with survival and improves the accuracy of common clinical variables, molecular features, and cancer subtypes. Prognostic lncRNAs are often characterized by switch-like expression patterns. In low-grade gliomas, HOXA10-AS activation is a robust marker of poor prognosis that complements IDH1/2 mutations, as validated in another retrospective cohort, and correlates with developmental pathways in tumor transcriptomes. Loss- and gain-of-function studies in patient-derived glioma cells, organoids, and xenograft models identify HOXA10-AS as a potent onco-lncRNA that regulates cell proliferation, contact inhibition, invasion, Hippo signaling, and mitotic and neuro-developmental pathways. Our study underscores the pan-cancer potential of the non-coding transcriptome for identifying biomarkers and regulators of cancer progression.

Glioblastomas are hierarchically organised tumours driven by glioma stem cells that retain partial differentiation potential. Glioma stem cells are maintained in specialised microenvironments, but whether, or how, they undergo lineage progression outside of these niches remains unclear. Here we identify the white matter as a differentiative niche for glioblastomas with oligodendrocyte lineage competency. Tumour cells in contact with white matter acquire pre-oligodendrocyte fate, resulting in decreased proliferation and invasion. Differentiation is a response to white matter injury, which is caused by tumour infiltration itself in a tumoursuppressive feedback loop. Mechanistically, tumour cell differentiation is driven by selective white matter upregulation of SOX10, a master regulator of normal oligodendrogenesis. SOX10 overexpression or treatment with myelination-promoting agents that upregulate endogenous SOX10, mimic this response, leading to niche-independent pre-oligodendrocyte differentiation and tumour suppression in vivo. Thus, glioblastoma recapitulates an injury response and exploiting this latent programme may offer treatment opportunities for a subset of patients.

Across the animal kingdom, adult tissue homeostasis is regulated by adult stem cell activity, which is commonly dysregulated in human cancers. However, identifying key regulators of stem cells in the milieu of thousands of genes dysregulated in a given cancer is challenging. Here, using a comparative genomics approach between planarian adult stem cells and patient-derived glioblastoma stem cells (GSCs), we identify and demonstrate the role of DEAD-box helicase DDX56 in regulating aspects of stemness in four stem cell systems: planarians, mouse neural stem cells, human GSCs, and a fly model of glioblastoma. In a human GSC line, DDX56 localizes to the nucleolus, and using planarians, when DDX56 is lost, stem cells dysregulate expression of ribosomal RNAs and lose nucleolar integrity prior to stem cell death. Together, a comparative genomic approach can be used to uncover conserved stemness regulators that are functional in both normal and cancer stem cells.

Pediatric ependymomas are highly recurrent tumors resistant to conventional chemotherapy. Telomerase, a ribonucleoprotein critical in permitting limitless replication, has been found to be critically important for the maintenance of tumor-initiating cells (TICs). These TICs are chemoresistant, repopulate the tumor from which they are identified, and are drivers of recurrence in numerous cancers. In this study, telomerase enzymatic activity was directly measured and inhibited to assess the therapeutic potential of targeting telomerase. Telomerase repeat amplification protocol (TRAP) (n = 36) and C-circle assay/telomere FISH/ATRX staining (n = 76) were performed on primary ependymomas to determine the prevalence and prognostic potential of telomerase activity or alternative lengthening of telomeres (ALT) as telomere maintenance mechanisms, respectively. Imetelstat, a phase 2 telomerase inhibitor, was used to elucidate the effect of telomerase inhibition on proliferation and tumorigenicity in established cell lines (BXD-1425EPN, R254), a primary TIC line (E520) and xenograft models of pediatric ependymoma. Over 60 % of pediatric ependymomas were found to rely on telomerase activity to maintain telomeres, while no ependymomas showed evidence of ALT. Children with telomerase-active tumors had reduced 5-year progression-free survival (29 ± 11 vs 64 ± 18 %; p = 0.03) and overall survival (58 ± 12 vs 83 ± 15 %; p = 0.05) rates compared to those with tumors lacking telomerase activity. Imetelstat inhibited proliferation and self-renewal by shortening telomeres and inducing senescence in vitro. In vivo, Imetelstat significantly reduced subcutaneous xenograft growth by 40 % (p = 0.03) and completely abolished the tumorigenicity of pediatric ependymoma TICs in an orthotopic xenograft model. Telomerase inhibition represents a promising therapeutic approach for telomerase-active pediatric ependymomas found to characterize high-risk ependymomas.

We recently reported that atypical teratoid rhabdoid tumors (ATRTs) comprise at least two transcriptional subtypes with different clinical outcomes; however, the mechanisms underlying therapeutic heterogeneity remained unclear. In this study, we analyzed 191 primary ATRTs and 10 ATRT cell lines to define the genomic and epigenomic landscape of ATRTs and identify subgroup-specific therapeutic targets. We found ATRTs segregated into three epigenetic subgroups with distinct genomic profiles, SMARCB1 genotypes, and chromatin landscape that correlated with differential cellular responses to a panel of signaling and epigenetic inhibitors. Significantly, we discovered that differential methylation of a PDGFRB-associated enhancer confers specific sensitivity of group 2 ATRT cells to dasatinib and nilotinib, and suggest that these are promising therapies for this highly lethal ATRT subtype.

Glioblastoma (GBM) is a deadly cancer in which cancer stem cells (CSCs) sustain tumor growth and contribute to therapeutic resistance. Protein arginine methyltransferase 5 (PRMT5) has recently emerged as a promising target in GBM. Using two orthogonal-acting inhibitors of PRMT5 (GSK591 or LLY-283), we show that pharmacological inhibition of PRMT5 suppresses the growth of a cohort of 46 patient-derived GBM stem cell cultures, with the proneural subtype showing greater sensitivity. We show that PRMT5 inhibition causes widespread disruption of splicing across the transcriptome, particularly affecting cell cycle gene products. We identify a GBM splicing signature that correlates with the degree of response to PRMT5 inhibition. Importantly, we demonstrate that LLY-283 is brain-penetrant and significantly prolongs the survival of mice with orthotopic patient-derived xenografts. Collectively, our findings provide a rationale for the clinical development of brain penetrant PRMT5 inhibitors as treatment for GBM.

Glioblastoma (GBM) is a cancer comprised of morphologically, genetically, and phenotypically diverse cells. However, an understanding of the functional significance of intratumoral heterogeneity is lacking. We devised a method to isolate and functionally profile tumorigenic clones from patient glioblastoma samples. Individual clones demonstrated unique proliferation and differentiation abilities. Importantly, naïve patient tumors included clones that were temozolomide resistant, indicating that resistance to conventional GBM therapy can preexist in untreated tumors at a clonal level. Further, candidate therapies for resistant clones were detected with clone-specific drug screening. Genomic analyses revealed genes and pathways that associate with specific functional behavior of single clones. Our results suggest that functional clonal profiling used to identify tumorigenic and drug-resistant tumor clones will lead to the discovery of new GBM clone-specific treatment strategies.