The current consensus recognizes four main medulloblastoma subgroups (wingless, Sonic hedgehog, group 3 and group 4). While medulloblastoma subgroups have been characterized extensively at the (epi-)genomic and transcriptomic levels, the proteome and phosphoproteome landscape remain to be comprehensively elucidated. Using quantitative (phospho)-proteomics in primary human medulloblastomas, we unravel distinct posttranscriptional regulation leading to highly divergent oncogenic signaling and kinase activity profiles in groups 3 and 4 medulloblastomas. Specifically, proteomic and phosphoproteomic analyses identify aberrant ERBB4-SRC signaling in group 4. Hence, enforced expression of an activated SRC combined with p53 inactivation induces murine tumors that resemble group 4 medulloblastoma. Therefore, our integrative proteogenomics approach unveils an oncogenic pathway and potential therapeutic vulnerability in the most common medulloblastoma subgroup.

Two recurrent mutations, K27M and G34R/V, within histone variant H3.3 were recently identified in ∼50% of pHGGs. Both mutations define clinically and biologically distinct subgroups of pHGGs. Here, we provide further insight about the dominant-negative effect of K27M mutant H3.3, leading to a global reduction of the repressive histone mark H3K27me3. We demonstrate that this is caused by aberrant recruitment of the PRC2 complex to K27M mutant H3.3 and enzymatic inhibition of the H3K27me3-establishing methyltransferase EZH2. By performing chromatin immunoprecipitation followed by next-generation sequencing and whole-genome bisulfite sequencing in primary pHGGs, we show that reduced H3K27me3 levels and DNA hypomethylation act in concert to activate gene expression in K27M mutant pHGGs.

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.

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.

Blastemal histology in chemotherapy-treated pediatric Wilms tumors (nephroblastoma) is associated with adverse prognosis. To uncover the underlying tumor biology and find therapeutic leads for this subgroup, we analyzed 58 blastemal type Wilms tumors by exome and transcriptome sequencing and validated our findings in a large replication cohort. Recurrent mutations included a hotspot mutation (Q177R) in the homeo-domain of SIX1 and SIX2 in tumors with high proliferative potential (18.1% of blastemal cases); mutations in the DROSHA/DGCR8 microprocessor genes (18.2% of blastemal cases); mutations in DICER1 and DIS3L2; and alterations in IGF2, MYCN, and TP53, the latter being strongly associated with dismal outcome. DROSHA and DGCR8 mutations strongly altered miRNA expression patterns in tumors, which was functionally validated in cell lines expressing mutant DROSHA.

Biallelic inactivation of SMARCB1, encoding a member of the SWI/SNF chromatin remodeling complex, is the hallmark genetic aberration of atypical teratoid rhabdoid tumors (ATRT). Here, we report how loss of SMARCB1 affects the epigenome in these tumors. Using chromatin immunoprecipitation sequencing (ChIP-seq) on primary tumors for a series of active and repressive histone marks, we identified the chromatin states differentially represented in ATRTs compared with other brain tumors and non-neoplastic brain. Re-expression of SMARCB1 in ATRT cell lines enabled confirmation of our genome-wide findings for the chromatin states. Additional generation of ChIP-seq data for SWI/SNF and Polycomb group proteins and the transcriptional repressor protein REST determined differential dependencies of SWI/SNF and Polycomb complexes in regulation of diverse gene sets in ATRTs.

Ependymal tumors across age groups are currently classified and graded solely by histopathology. It is, however, commonly accepted that this classification scheme has limited clinical utility based on its lack of reproducibility in predicting patients' outcome. We aimed at establishing a uniform molecular classification using DNA methylation profiling. Nine molecular subgroups were identified in a large cohort of 500 tumors, 3 in each anatomical compartment of the CNS, spine, posterior fossa, supratentorial. Two supratentorial subgroups are characterized by prototypic fusion genes involving RELA and YAP1, respectively. Regarding clinical associations, the molecular classification proposed herein outperforms the current histopathological classification and thus might serve as a basis for the next World Health Organization classification of CNS tumors.

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.

Ependymoma is a heterogeneous entity of central nervous system tumors with well-established molecular groups. Here, we apply single-cell RNA sequencing to analyze ependymomas across molecular groups and anatomic locations to investigate their intratumoral heterogeneity and developmental origins. Ependymomas are composed of a cellular hierarchy initiating from undifferentiated populations, which undergo impaired differentiation toward three lineages of neuronal-glial fate specification. While prognostically favorable groups of ependymoma predominantly harbor differentiated cells, aggressive groups are enriched for undifferentiated cell populations. The delineated transcriptomic signatures correlate with patient survival and define molecular dependencies for targeted treatment approaches. Taken together, our analyses reveal a developmental hierarchy underlying ependymomas relevant to biological and clinical behavior.

Glioblastoma (GBM) is a brain tumor that carries a dismal prognosis and displays considerable heterogeneity. We have recently identified recurrent H3F3A mutations affecting two critical amino acids (K27 and G34) of histone H3.3 in one-third of pediatric GBM. Here, we show that each H3F3A mutation defines an epigenetic subgroup of GBM with a distinct global methylation pattern, and that they are mutually exclusive with IDH1 mutations, which characterize a third mutation-defined subgroup. Three further epigenetic subgroups were enriched for hallmark genetic events of adult GBM and/or established transcriptomic signatures. We also demonstrate that the two H3F3A mutations give rise to GBMs in separate anatomic compartments, with differential regulation of transcription factors OLIG1, OLIG2, and FOXG1, possibly reflecting different cellular origins.

Medulloblastoma (MB) is a highly malignant pediatric brain tumor. Despite aggressive therapy, many patients succumb to the disease, and survivors experience severe side effects from treatment. MYC-driven MB has a particularly poor prognosis and would greatly benefit from more effective therapies. We used an animal model of MYC-driven MB to screen for drugs that decrease viability of tumor cells. Among the most effective compounds were histone deacetylase inhibitors (HDACIs). HDACIs potently inhibit survival of MYC-driven MB cells in vitro, in part by inducing expression of the FOXO1 tumor suppressor gene. HDACIs also synergize with phosphatidylinositol 3-kinase inhibitors to inhibit tumor growth in vivo. These studies identify an effective combination therapy for the most aggressive form of MB.

There is a pressing need to identify therapeutic targets in tumors with low mutation rates such as the malignant pediatric brain tumor medulloblastoma. To address this challenge, we quantitatively profiled global proteomes and phospho-proteomes of 45 medulloblastoma samples. Integrated analyses revealed that tumors with similar RNA expression vary extensively at the post-transcriptional and post-translational levels. We identified distinct pathways associated with two subsets of SHH tumors, and found post-translational modifications of MYC that are associated with poor outcomes in group 3 tumors. We found kinases associated with subtypes and showed that inhibiting PRKDC sensitizes MYC-driven cells to radiation. Our study shows that proteomics enables a more comprehensive, functional readout, providing a foundation for future therapeutic strategies.