Recurrent medulloblastoma is a therapeutic challenge because it is almost always fatal. Studies have confirmed that medulloblastoma consists of at least four distinct subgroups. We sought to delineate subgroup-specific differences in medulloblastoma recurrence patterns.

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.

Cranial radiotherapy improves survival of the most common childhood cancers, including brain tumors and leukemia. Unfortunately, long-term survivors are faced with consequences of secondary neoplasia, including radiation-induced meningiomas (RIMs). We characterized 31 RIMs with exome/NF2 intronic sequencing, RNA sequencing and methylation profiling, and found NF2 gene rearrangements in 12/31 of RIMs, an observation previously unreported in sporadic meningioma (SM). Additionally, known recurrent mutations characteristic of SM, including AKT1, KLF4, TRAF7 and SMO, were not observed in RIMs. Combined losses of chromosomes 1p and 22q were common in RIMs (16/18 cases) and overall, chromosomal aberrations were more complex than that observed in SM. Patterns of DNA methylation profiling supported similar cell of origin between RIMs and SMs. The findings indicate that the mutational landscape of RIMs is distinct from SMs, and have significant therapeutic implications for survivors of childhood cranial radiation and the elucidation of the molecular pathogenesis of meningiomas.Radiation-induced meningiomas are often more aggressive than sporadic ones. In this study, the authors perform an exome, methylation and RNA-seq analysis of 31 cases of radiation-induced meningioma and show NF2 rearrangement, an observation previously unreported in the sporadic tumors.

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 present an extensive assessment of mutation burden through sequencing analysis of >81,000 tumors from pediatric and adult patients, including tumors with hypermutation caused by chemotherapy, carcinogens, or germline alterations. Hypermutation was detected in tumor types not previously associated with high mutation burden. Replication repair deficiency was a major contributing factor. We uncovered new driver mutations in the replication-repair-associated DNA polymerases and a distinct impact of microsatellite instability and replication repair deficiency on the scale of mutation load. Unbiased clustering, based on mutational context, revealed clinically relevant subgroups regardless of the tumors' tissue of origin, highlighting similarities in evolutionary dynamics leading to hypermutation. Mutagens, such as UV light, were implicated in unexpected cancers, including sarcomas and lung tumors. The order of mutational signatures identified previous treatment and germline replication repair deficiency, which improved management of patients and families. These data will inform tumor classification, genetic testing, and clinical trial design.

The emergence of high-throughput, next-generation sequencing technologies has dramatically altered the way we assess genomes in population genetics and in cancer genomics. Currently, there are four commonly used whole-genome sequencing platforms on the market: Illumina's HiSeq2000, Life Technologies' SOLiD 4 and its completely redesigned 5500xl SOLiD, and Complete Genomics' technology. A number of earlier studies have compared a subset of those sequencing platforms or compared those platforms with Sanger sequencing, which is prohibitively expensive for whole genome studies. Here we present a detailed comparison of the performance of all currently available whole genome sequencing platforms, especially regarding their ability to call SNVs and to evenly cover the genome and specific genomic regions. Unlike earlier studies, we base our comparison on four different samples, allowing us to assess the between-sample variation of the platforms. We find a pronounced GC bias in GC-rich regions for Life Technologies' platforms, with Complete Genomics performing best here, while we see the least bias in GC-poor regions for HiSeq2000 and 5500xl. HiSeq2000 gives the most uniform coverage and displays the least sample-to-sample variation. In contrast, Complete Genomics exhibits by far the smallest fraction of bases not covered, while the SOLiD platforms reveal remarkable shortcomings, especially in covering CpG islands. When comparing the performance of the four platforms for calling SNPs, HiSeq2000 and Complete Genomics achieve the highest sensitivity, while the SOLiD platforms show the lowest false positive rate. Finally, we find that integrating sequencing data from different platforms offers the potential to combine the strengths of different technologies. In summary, our results detail the strengths and weaknesses of all four whole-genome sequencing platforms. It indicates application areas that call for a specific sequencing platform and disallow other platforms. This helps to identify the proper sequencing platform for whole genome studies with different application scopes.

Tumors of the same histologic type often comprise clinically and molecularly distinct subgroups; however, the etiology of these subgroups is unknown. Here, we report that histologically identical, but genetically distinct, ependymomas exhibit patterns of gene expression that recapitulate those of radial glia cells in the corresponding region of the central nervous system. Cancer stem cells isolated from ependymomas displayed a radial glia phenotype and formed tumors when orthotopically transplanted in mice. These findings identify restricted populations of radial glia cells as candidate stem cells of the different subgroups of ependymoma, and they support a general hypothesis that subgroups of the same histologic tumor type are generated by different populations of progenitor cells in the tissues of origin.

Serial Analysis of Gene Expression (SAGE) and microarrays have found a widespread application, but much ambiguity exists regarding the evaluation of these technologies. Cross-platform utilization of gene expression data from the SAGE and microarray technology could reduce the need for duplicate experiments and facilitate a more extensive exchange of data within the research community. This requires a measure for the correspondence of the different gene expression platforms. To date, a number of cross-platform evaluations (including a few studies using SAGE and Affymetrix GeneChips) have been conducted showing a variable, but overall low, concordance. This study evaluates these overall measures and introduces the between-ratio difference as a concordance measure pergene.

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.

Medulloblastomas comprise a heterogeneous group of tumours and can be subdivided into four molecular subgroups (WNT, SHH, Group 3 and Group 4) with distinct prognosis, biological behaviour and implications for targeted therapies. Few experimental models exist of the aggressive and poorly characterized Group 3 tumours. In order to establish a reproducible transplantable Group 3 medulloblastoma model for preclinical therapeutic studies, we acquired a patient-derived tumour sphere culture and inoculated low-passage spheres into the cerebellums of NOD-scid mice. Mice developed symptoms of brain tumours with a latency of 17-18 weeks. Neurosphere cultures were re-established and serially transplanted for 3 generations, with a negative correlation between tumour latency and numbers of injected cells. Xenografts replicated the phenotype of the primary tumour, including high degree of clustering in DNA methylation analysis, high proliferation, expression of tumour markers, MYC amplification and elevated MYC expression, and sensitivity to the MYC inhibitor JQ1. Xenografts maintained maintained expression of tumour-derived VEGFA and stromal-derived COX-2. VEGFA, COX-2 and c-Myc are highly expressed in Group 3 compared to other medulloblastoma subgroups, suggesting that these molecules are relevant therapeutic targets in Group 3 medulloblastoma.

Human cells have twenty-three pairs of chromosomes. In cancer, however, genes can be amplified in chromosomes or in circular extrachromosomal DNA (ecDNA), although the frequency and functional importance of ecDNA are not understood. We performed whole-genome sequencing, structural modelling and cytogenetic analyses of 17 different cancer types, including analysis of the structure and function of chromosomes during metaphase of 2,572 dividing cells, and developed a software package called ECdetect to conduct unbiased, integrated ecDNA detection and analysis. Here we show that ecDNA was found in nearly half of human cancers; its frequency varied by tumour type, but it was almost never found in normal cells. Driver oncogenes were amplified most commonly in ecDNA, thereby increasing transcript level. Mathematical modelling predicted that ecDNA amplification would increase oncogene copy number and intratumoural heterogeneity more effectively than chromosomal amplification. We validated these predictions by quantitative analyses of cancer samples. The results presented here suggest that ecDNA contributes to accelerated evolution in cancer.

Medulloblastoma comprises four main subgroups (WNT, SHH, Group 3 and Group 4) originally defined by transcriptional profiling. In primary medulloblastoma tissues, these groups are thought to be distinguishable using the immunohistochemical detection of β-catenin, filamin A, GAB1 and YAP1 protein markers. To investigate the utility of these markers for in vitro studies using medulloblastoma cell lines, immunoblotting and indirect immunofluorescence were employed for the detection of β-catenin, filamin A, GAB1 and YAP1 in both DAOY and D283 Med reference cell lines and the panel of six medulloblastoma cell lines derived in our laboratory from the primary tumor tissues of known molecular subgroups. Immunohistochemical detection of these markers was performed on formalin-fixed paraffin-embedded tissue of the matching primary tumors. The results revealed substantial divergences between the primary tumor tissues and matching cell lines in the immunoreactivity pattern of medulloblastoma-subgroup-specific protein markers. Regardless of the molecular subgroup of the primary tumor, all six patient-derived medulloblastoma cell lines exhibited a uniform phenotype: immunofluorescence showed the nuclear localization of YAP1, accompanied by strong cytoplasmic positivity for β-catenin and filamin A, as well as weak positivity for GAB1. The same immunoreactivity pattern was also found in both DAOY and D283 Med reference medulloblastoma cell lines. Therefore, we can conclude that various medulloblastoma cell lines tend to exhibit the same characteristics of protein marker expression under standard in vitro conditions. Such a finding emphasizes the importance of the analyses of primary tumors in clinically oriented medulloblastoma research and the urgent need to develop in vitro models of improved clinical relevance, such as 3D cultures and organotypic slice cultures.

Ependymoma with RELA fusion has been defined as a novel entity of the revised World Health Organization 2016 classification of tumors of the central nervous system (CNS), characterized by fusion transcripts of the RELA gene and consequent pathological activation of the NFkB pathway. These tumors represent the majority of supratentorial ependymomas in children. The validation of diagnostic tools to identify this clinically relevant ependymoma entity is essential. Here, we have used interphase fluorescent in situ hybridization (FISH) for C11orf95 and RELA, immunohistochemistry (IHC) for p65-RelA and the recently developed DNA methylation-based classification besides conventional histopathology, and compared the precision of the methods in 40 supratentorial pediatric brain tumors diagnosed as ependymomas in the past years. Reverse transcription PCR (RT-PCR) and RNA sequencing were performed to explore discordant cases. Furthermore, we integrated imaging and clinical features as additional layers of information. The concordance between nuclear RelA expression by IHC and RELA FISH was 100%. Concordance between IHC and DNA methylation profiling, and between FISH and DNA methylation profiling was also high (96.4% and 95.2%, respectively). Thirty-four out of 40 (85%) cases were confirmed by integrated diagnoses as ependymal tumors, including 22 RELA-fused ependymomas (71% of ependymal tumors), two YAP1-fused ependymomas (6%), six non-RELA/non-YAP1 ependymomas (18%) and four ependymal/subependymal mixed tumors (12%). Ependymal/subependymal mixed tumors had an excellent clinical outcome despite the presence of histopathological signs of malignancy, suggesting that these tumors should not be diagnosed as classic ependymomas. DNA methylation profiling helped in the differential diagnosis of RELA-fused ependymomas. IHC and FISH, which are available in the majority of pathology laboratories, are valuable tools to identify RELA-fused ependymomas.

Over the past decade, wingless-activated (WNT) medulloblastoma has been identified as a candidate for therapy de-escalation based on excellent survival; however, a paucity of relapses has precluded additional analyses of markers of relapse. To address this gap in knowledge, an international cohort of 93 molecularly confirmed WNT MB was assembled, where 5-year progression-free survival is 0.84 (95%, 0.763-0.925) with 15 relapsed individuals identified. Maintenance chemotherapy is identified as a strong predictor of relapse, with individuals receiving high doses of cyclophosphamide or ifosphamide having only one very late molecularly confirmed relapse (p = 0.032). The anatomical location of recurrence is metastatic in 12 of 15 relapses, with 8 of 12 metastatic relapses in the lateral ventricles. Maintenance chemotherapy, specifically cumulative cyclophosphamide doses, is a significant predictor of relapse across WNT MB. Future efforts to de-escalate therapy need to carefully consider not only the radiation dose but also the chemotherapy regimen and the propensity for metastatic relapses.

A molecular diagnostic method that incorporates information about the transcriptional status of all genes across multiple tissue types can strengthen confidence in cancer diagnosis.

Spinal ependymal tumors form a histologically and molecularly heterogeneous group of tumors with generally good prognosis. However, their treatment can be challenging if infiltration of the spinal cord or dissemination throughout the central nervous system (CNS) occurs and, in these cases, clinical outcome remains poor. Here, we describe a new and relatively rare subgroup of spinal ependymal tumors identified using DNA methylation profiling that is distinct from other molecular subgroups of ependymoma. Copy number variation plots derived from DNA methylation arrays showed MYCN amplification as a characteristic genetic alteration in all cases of our cohort (n = 13), which was subsequently validated using fluorescence in situ hybridization. The histological diagnosis was anaplastic ependymoma (WHO Grade III) in ten cases and classic ependymoma (WHO Grade II) in three cases. Histological re-evaluation in five primary tumors and seven relapses showed characteristic histological features of ependymoma, namely pseudorosettes, GFAP- and EMA positivity. Electron microscopy revealed cilia, complex intercellular junctions and intermediate filaments in a representative sample. Taking these findings into account, we suggest to designate this molecular subgroup spinal ependymoma with MYCN amplification, SP-EPN-MYCN. SP-EPN-MYCN tumors showed distinct growth patterns with intradural, extramedullary localization mostly within the thoracic and cervical spine, diffuse leptomeningeal spread throughout the whole CNS and infiltrative invasion of the spinal cord. Dissemination was observed in 100% of cases. Despite high-intensity treatment, SP-EPN-MYCN showed significantly worse median progression free survival (PFS) (17 months) and median overall survival (OS) (87 months) than all other previously described molecular spinal ependymoma subgroups. OS and PFS were similar to supratentorial ependymoma with RELA-fusion (ST-EPN-RELA) and posterior fossa ependymoma A (PF-EPN-A), further highlighting the aggressiveness of this distinct new subgroup. We, therefore, propose to establish SP-EPN-MYCN as a new molecular subgroup in ependymoma and advocate for testing newly diagnosed spinal ependymal tumors for MYCN amplification.

DEAD-Box Helicase 3 X-Linked (DDX3X) is frequently mutated in the Wingless (WNT) and Sonic hedghog (SHH) subtypes of medulloblastoma-the commonest malignant childhood brain tumor, but whether DDX3X functions as a medulloblastoma oncogene or tumor suppressor gene is not known. Here, we show that Ddx3x regulates hindbrain patterning and development by controlling Hox gene expression and cell stress signaling. In mice predisposed to Wnt- or Shh medulloblastoma, Ddx3x sensed oncogenic stress and suppressed tumor formation. WNT and SHH medulloblastomas normally arise only in the lower and upper rhombic lips, respectively. Deletion of Ddx3x removed this lineage restriction, enabling both medulloblastoma subtypes to arise in either germinal zone. Thus, DDX3X is a medulloblastoma tumor suppressor that regulates hindbrain development and restricts the competence of cell lineages to form medulloblastoma subtypes.

High-intensity therapy effectively treats most TP53 wild-type (TP53-WT) Sonic Hedgehog-subgroup medulloblastomas (SHH-MBs), but often cause long-term deleterious neurotoxicities in children. Recent clinical trials investigating reduction/de-escalation of therapy for TP53-WT SHH-MBs caused poor overall survival. Here, we investigated whether reduced levels of p53-pathway activation by low-intensity therapy potentially contribute to diminished therapeutic efficacy.

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.

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.