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.

We describe molecular convergence between BMI1 and CHD7 in the initiation of medulloblastoma. Identified in a functional genomic screen in mouse models, a BMI1High;CHD7Low expression signature within medulloblastoma characterizes patients with poor overall survival. We show that BMI1-mediated repression of the ERK1/2 pathway leads to increased proliferation and tumor burden in primary human MB cells and in a xenograft model, respectively. We provide evidence that repression of the ERK inhibitor DUSP4 by BMI1 is dependent on a more accessible chromatin configuration in G4 MB cells with low CHD7 expression. These findings extend current knowledge of the role of BMI1 and CHD7 in medulloblastoma pathogenesis, and they raise the possibility that pharmacological targeting of BMI1 or ERK may be particularly indicated in a subgroup of MB with low expression levels of CHD7.

Study of the origin and development of cerebellar tumours has been hampered by the complexity and heterogeneity of cerebellar cells that change over the course of development. Here we use single-cell transcriptomics to study more than 60,000 cells from the developing mouse cerebellum and show that different molecular subgroups of childhood cerebellar tumours mirror the transcription of cells from distinct, temporally restricted cerebellar lineages. The Sonic Hedgehog medulloblastoma subgroup transcriptionally mirrors the granule cell hierarchy as expected, while group 3 medulloblastoma resembles Nestin+ stem cells, group 4 medulloblastoma resembles unipolar brush cells, and PFA/PFB ependymoma and cerebellar pilocytic astrocytoma resemble the prenatal gliogenic progenitor cells. Furthermore, single-cell transcriptomics of human childhood cerebellar tumours demonstrates that many bulk tumours contain a mixed population of cells with divergent differentiation. Our data highlight cerebellar tumours as a disorder of early brain development and provide a proximate explanation for the peak incidence of cerebellar tumours in early childhood.

Although brain functionality is often remarkably robust to lesions and other insults, it may be fragile when these take place in specific locations. Previous attempts to quantify robustness and fragility sought to understand how the functional connectivity of brain networks is affected by structural changes, using either model-based predictions or empirical studies of the effects of lesions. We advance a geometric viewpoint relying on a notion of network curvature, the so-called Ollivier-Ricci curvature. This approach has been proposed to assess financial market robustness and to differentiate biological networks of cancer cells from healthy ones. Here, we apply curvature-based measures to brain structural networks to identify robust and fragile brain regions in healthy subjects. We show that curvature can also be used to track changes in brain connectivity related to age and autism spectrum disorder (ASD), and we obtain results that are in agreement with previous MRI studies.

In cancer, recurrent somatic single-nucleotide variants-which are rare in most paediatric cancers-are confined largely to protein-coding genes1-3. Here we report highly recurrent hotspot mutations (r.3A>G) of U1 spliceosomal small nuclear RNAs (snRNAs) in about 50% of Sonic hedgehog (SHH) medulloblastomas. These mutations were not present across other subgroups of medulloblastoma, and we identified these hotspot mutations in U1 snRNA in only <0.1% of 2,442 cancers, across 36 other tumour types. The mutations occur in 97% of adults (subtype SHHδ) and 25% of adolescents (subtype SHHα) with SHH medulloblastoma, but are largely absent from SHH medulloblastoma in infants. The U1 snRNA mutations occur in the 5' splice-site binding region, and snRNA-mutant tumours have significantly disrupted RNA splicing and an excess of 5' cryptic splicing events. Alternative splicing mediated by mutant U1 snRNA inactivates tumour-suppressor genes (PTCH1) and activates oncogenes (GLI2 and CCND2), and represents a target for therapy. These U1 snRNA mutations provide an example of highly recurrent and tissue-specific mutations of a non-protein-coding gene in cancer.

Spatial heterogeneity of transcriptional and genetic markers between physically isolated biopsies of a single tumor poses major barriers to the identification of biomarkers and the development of targeted therapies that will be effective against the entire tumor. We analyzed the spatial heterogeneity of multiregional biopsies from 35 patients, using a combination of transcriptomic and genomic profiles. Medulloblastomas (MBs), but not high-grade gliomas (HGGs), demonstrated spatially homogeneous transcriptomes, which allowed for accurate subgrouping of tumors from a single biopsy. Conversely, somatic mutations that affect genes suitable for targeted therapeutics demonstrated high levels of spatial heterogeneity in MB, malignant glioma, and renal cell carcinoma (RCC). Actionable targets found in a single MB biopsy were seldom clonal across the entire tumor, which brings the efficacy of monotherapies against a single target into question. Clinical trials of targeted therapies for MB should first ensure the spatially ubiquitous nature of the target mutation.

The robustness of brain structural networks, estimated from diffusion MRI data, may be relevant to cognition. We investigate whether measures of network robustness, such as Ollivier-Ricci curvature, can explain cognitive impairment in multiple sclerosis (MS). We assessed whether local (i.e., cortical area) and/or global (i.e., whole brain) robustness, differs between cognitively impaired (MSCI) and non-impaired (MSNI) MS patients. Fifty patients, with Expanded Disability Status Scale mean (m): 3.2, disease duration m: 12 years, and age m: 40 years, were enrolled. Cognitive impairment scores were estimated from the Minimal Assessment of Cognitive Function in Multiple Sclerosis. Images were obtained in a 3T MRI using a diffusion protocol with a 2 min acquisition time. Brain structural networks were created using 333 cortical areas. Local and global robustness was estimated for each individual, and comparisons were performed between MSCI and MSNI patients. 31 MSCI and 10 MSNI patients were included in the analyses. Brain structural network robustness and centrality showed significant correlations with cognitive impairment. Measures of network robustness and centrality identified specific cortical areas relevant to MS-related cognitive impairment. These measures can be obtained on clinical scanners and are succinct yet accurate potential biomarkers of cognitive impairment.

Diffusion MRI (dMRI) reveals microstructural features of the brain white matter by quantifying the anisotropic diffusion of water molecules within axonal bundles. Yet, identifying features such as axonal orientation dispersion, density, diameter, etc., in complex white matter fiber configurations (e.g. crossings) has proved challenging. Besides optimized data acquisition and advanced biophysical models, computational procedures to fit such models to the data are critical. However, these procedures have been largely overlooked by the dMRI microstructure community and new, more versatile, approaches are needed to solve complex biophysical model fitting problems. Existing methods are limited to models assuming single fiber orientation, relevant to limited brain areas like the corpus callosum, or multiple orientations but without the ability to extract detailed microstructural features. Here, we introduce a new and versatile optimization technique (MIX), which enables microstructure imaging of crossing white matter fibers. We provide a MATLAB implementation of MIX, and demonstrate its applicability to general microstructure models in fiber crossings using synthetic as well as ex-vivo and in-vivo brain data.

Recurrent medulloblastoma and ependymoma are universally lethal, with no approved targeted therapies and few candidates presently under clinical evaluation. Nearly all recurrent medulloblastomas and posterior fossa group A (PFA) ependymomas are located adjacent to and bathed by the cerebrospinal fluid, presenting an opportunity for locoregional therapy, bypassing the blood-brain barrier. We identify three cell-surface targets, EPHA2, HER2 and interleukin 13 receptor α2, expressed on medulloblastomas and ependymomas, but not expressed in the normal developing brain. We validate intrathecal delivery of EPHA2, HER2 and interleukin 13 receptor α2 chimeric antigen receptor T cells as an effective treatment for primary, metastatic and recurrent group 3 medulloblastoma and PFA ependymoma xenografts in mouse models. Finally, we demonstrate that administration of these chimeric antigen receptor T cells into the cerebrospinal fluid, alone or in combination with azacytidine, is a highly effective therapy for multiple metastatic mouse models of group 3 medulloblastoma and PFA ependymoma, thereby providing a rationale for clinical trials of these approaches in humans.

Gliomas demonstrate epigenetic dysregulation exemplified by the Glioma CpG Island Methylator Phenotype (G-CIMP) seen in IDH1 mutant tumors. 5-Hydroxymethylcytosine (5hmC) is implicated in glioma pathogenesis; however, its role in IDH1 mutant gliomas is incompletely understood. To characterize 5hmC in IDH1 mutant gliomas further, we examine 5hmC in a cohort of IDH1 mutant and wild-type high-grade gliomas (HGG) using a quantitative locus-specific approach. Regions demonstrating high 5hmC abundance and differentially hydroxymethylated regions (DHMR) enrich for enhancers implicated in glioma pathogenesis. Among these regions, IDH1 mutant tumors possess greater 5hmC compared to wild type. 5hmC contributes to overall methylation status of G-CIMP genes. 5hmC targeting gene body regions correlates significantly with increased gene expression. In particular, a strong correlation between increased 5hmC and increased gene expression is identified for genes highly expressed in the IDH1 mutant cohort. Overall, locus-specific gain of 5hmC targeting regulatory regions and associated with overexpressed genes suggests a significant role for 5hmC in IDH1 mutant HGG.

While the preponderance of morbidity and mortality in medulloblastoma patients are due to metastatic disease, most research focuses on the primary tumor due to a dearth of metastatic tissue samples and model systems. Medulloblastoma metastases are found almost exclusively on the leptomeningeal surface of the brain and spinal cord; dissemination is therefore thought to occur through shedding of primary tumor cells into the cerebrospinal fluid followed by distal re-implantation on the leptomeninges. We present evidence for medulloblastoma circulating tumor cells (CTCs) in therapy-naive patients and demonstrate in vivo, through flank xenografting and parabiosis, that medulloblastoma CTCs can spread through the blood to the leptomeningeal space to form leptomeningeal metastases. Medulloblastoma leptomeningeal metastases express high levels of the chemokine CCL2, and expression of CCL2 in medulloblastoma in vivo is sufficient to drive leptomeningeal dissemination. Hematogenous dissemination of medulloblastoma offers a new opportunity to diagnose and treat lethal disseminated medulloblastoma.

Genomic sequencing has driven precision-based oncology therapy; however, the genetic drivers of many malignancies remain unknown or non-targetable, so alternative approaches to the identification of therapeutic leads are necessary. Ependymomas are chemotherapy-resistant brain tumours, which, despite genomic sequencing, lack effective molecular targets. Intracranial ependymomas are segregated on the basis of anatomical location (supratentorial region or posterior fossa) and further divided into distinct molecular subgroups that reflect differences in the age of onset, gender predominance and response to therapy. The most common and aggressive subgroup, posterior fossa ependymoma group A (PF-EPN-A), occurs in young children and appears to lack recurrent somatic mutations. Conversely, posterior fossa ependymoma group B (PF-EPN-B) tumours display frequent large-scale copy number gains and losses but have favourable clinical outcomes. More than 70% of supratentorial ependymomas are defined by highly recurrent gene fusions in the NF-κB subunit gene RELA (ST-EPN-RELA), and a smaller number involve fusion of the gene encoding the transcriptional activator YAP1 (ST-EPN-YAP1). Subependymomas, a distinct histologic variant, can also be found within the supratetorial and posterior fossa compartments, and account for the majority of tumours in the molecular subgroups ST-EPN-SE and PF-EPN-SE. Here we describe mapping of active chromatin landscapes in 42 primary ependymomas in two non-overlapping primary ependymoma cohorts, with the goal of identifying essential super-enhancer-associated genes on which tumour cells depend. Enhancer regions revealed putative oncogenes, molecular targets and pathways; inhibition of these targets with small molecule inhibitors or short hairpin RNA diminished the proliferation of patient-derived neurospheres and increased survival in mouse models of ependymomas. Through profiling of transcriptional enhancers, our study provides a framework for target and drug discovery in other cancers that lack known genetic drivers and are therefore difficult to treat.

Sonic hedgehog medulloblastoma encompasses a clinically and molecularly diverse group of cancers of the developing central nervous system. Here, we use unbiased sequencing of the transcriptome across a large cohort of 250 tumors to reveal differences among molecular subtypes of the disease, and demonstrate the previously unappreciated importance of non-coding RNA transcripts. We identify alterations within the cAMP dependent pathway (GNAS, PRKAR1A) which converge on GLI2 activity and show that 18% of tumors have a genetic event that directly targets the abundance and/or stability of MYCN. Furthermore, we discover an extensive network of fusions in focally amplified regions encompassing GLI2, and several loss-of-function fusions in tumor suppressor genes PTCH1, SUFU and NCOR1. Molecular convergence on a subset of genes by nucleotide variants, copy number aberrations, and gene fusions highlight the key roles of specific pathways in the pathogenesis of Sonic hedgehog medulloblastoma and open up opportunities for therapeutic intervention.

Posterior fossa A (PFA) ependymomas are lethal malignancies of the hindbrain in infants and toddlers. Lacking highly recurrent somatic mutations, PFA ependymomas are proposed to be epigenetically driven tumors for which model systems are lacking. Here we demonstrate that PFA ependymomas are maintained under hypoxia, associated with restricted availability of specific metabolites to diminish histone methylation, and increase histone demethylation and acetylation at histone 3 lysine 27 (H3K27). PFA ependymomas initiate from a cell lineage in the first trimester of human development that resides in restricted oxygen. Unlike other ependymomas, transient exposure of PFA cells to ambient oxygen induces irreversible cellular toxicity. PFA tumors exhibit a low basal level of H3K27me3, and, paradoxically, inhibition of H3K27 methylation specifically disrupts PFA tumor growth. Targeting metabolism and/or the epigenome presents a unique opportunity for rational therapy for infants with PFA ependymoma.

SUFU alterations are common in human Sonic Hedgehog (SHH) subgroup medulloblastoma (MB). However, its tumorigenic mechanisms have remained elusive. Here, we report that loss of Sufu alone is unable to induce MB formation in mice, due to insufficient Gli2 activation. Simultaneous loss of Spop, an E3 ubiquitin ligase targeting Gli2, restores robust Gli2 activation and induces rapid MB formation in Sufu knockout background. We also demonstrated a tumor-promoting role of Sufu in Smo-activated MB (∼60% of human SHH MB) by maintaining robust Gli activity. Having established Gli2 activation as a key driver of SHH MB, we report a comprehensive analysis of its targetome. Furthermore, we identified Atoh1 as a target and molecular accomplice of Gli2 that activates core SHH MB signature genes in a synergistic manner. Overall, our work establishes the dual role of SUFU in SHH MB and provides mechanistic insights into transcriptional regulation underlying Gli2-mediated SHH MB tumorigenesis.

Medulloblastoma is a malignant pediatric tumor that arises from neural progenitors in the cerebellum. Despite a five-year survival rate of ~70%, nearly all patients incur adverse side effects from current treatment strategies that drastically impact quality of life. Roughly one-third of medulloblastoma are driven by aberrant activation of the Sonic Hedgehog (SHH) signaling pathway. However, the scarcity of genetic mutations in medulloblastoma has led to investigation of other mechanisms contributing to cancer pathogenicity including epigenetic regulation of gene expression. Here, we show that Helicase, Lymphoid Specific (HELLS), a chromatin remodeler with epigenetic functions including DNA methylation and histone modification, is induced by Sonic Hedgehog (SHH) in SHH-dependent cerebellar progenitor cells and the developing murine cerebella. HELLS is also up-regulated in mouse and human SHH medulloblastoma. Others have shown that HELLS activity generally results in a repressive chromatin state. Our results demonstrate that increased expression of HELLS in our experimental systems is regulated by the oncogenic transcriptional regulator YAP1 downstream of Smoothened, the positive transducer of SHH signaling. Elucidation of HELLS as one of the downstream effectors of the SHH pathway may lead to novel targets for precision therapeutics with the promise of better outcomes for SHH medulloblastoma patients.