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.

High tumor mutational burden (TMB) is an emerging biomarker of sensitivity to immune checkpoint inhibitors and has been shown to be more significantly associated with response to PD-1 and PD-L1 blockade immunotherapy than PD-1 or PD-L1 expression, as measured by immunohistochemistry (IHC). The distribution of TMB and the subset of patients with high TMB has not been well characterized in the majority of cancer types.

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.

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.

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.

The transcription factor eomesodermin (EOMES) promotes interleukin (IL)-10 expression in CD4+ T cells, which has been linked to immunosuppressive and cytotoxic activities. We detected cytotoxic, programmed cell death protein-1 (PD-1) and EOMES co-expressing CD4+ T cells in lymph nodes (LNs) of patients with chronic lymphocytic leukemia (CLL) or diffuse large B-cell lymphoma. Transcriptome and flow cytometry analyses revealed that EOMES does not only drive IL-10 expression, but rather controls a unique transcriptional signature in CD4+ T cells, that is enriched in genes typical for T regulatory type 1 (TR1) cells. The TR1 cell identity of these CD4+ T cells was supported by their expression of interferon gamma and IL-10, as well as inhibitory receptors including PD-1. TR1 cells with cytotoxic capacity accumulate also in Eµ-TCL1 mice that develop CLL-like disease. Whereas wild-type CD4+ T cells control TCL1 leukemia development after adoptive transfer in leukopenic Rag2-/- mice, EOMES-deficient CD4+ T cells failed to do so. We further show that TR1 cell-mediated control of TCL1 leukemia requires IL-10 receptor (IL-10R) signaling, as Il10rb-deficient CD4+ T cells showed impaired antileukemia activity. Altogether, our data demonstrate that EOMES is indispensable for the development of IL-10-expressing, cytotoxic TR1 cells, which accumulate in LNs of CLL patients and control TCL1 leukemia in mice in an IL-10R-dependent manner.

Sarcomas are malignant soft tissue and bone tumours affecting adults, adolescents and children. They represent a morphologically heterogeneous class of tumours and some entities lack defining histopathological features. Therefore, the diagnosis of sarcomas is burdened with a high inter-observer variability and misclassification rate. Here, we demonstrate classification of soft tissue and bone tumours using a machine learning classifier algorithm based on array-generated DNA methylation data. This sarcoma classifier is trained using a dataset of 1077 methylation profiles from comprehensively pre-characterized cases comprising 62 tumour methylation classes constituting a broad range of soft tissue and bone sarcoma subtypes across the entire age spectrum. The performance is validated in a cohort of 428 sarcomatous tumours, of which 322 cases were classified by the sarcoma classifier. Our results demonstrate the potential of the DNA methylation-based sarcoma classification for research and future diagnostic applications.

Two distinct genetically defined entities of ependymoma arising in the supratentorial compartment are characterized by the presence of either a C11orf95-RELA or a YAP-MAMLD1 fusion, respectively. There is growing evidence that supratentorial ependymomas without these genetic features exist. In this study, we report on 18 pediatric non-RELA/non-YAP supratentorial ependymomas that were systematically characterized by means of their histology, immunophenotype, genetics, and epigenomics. Comprehensive molecular analyses included high-resolution copy number analysis, methylation profiling, analysis of fusion transcripts by Nanostring technology, and RNA sequencing. Based upon histological and immunohistochemical features two main patterns were identified-RELA-like (n = 9) and tanycytic ependymomas (n = 6). In the RELA-like group histologically assigned to WHO grade III and resembling RELA-fused ependymomas, tumors lacked nuclear expression of p65-RelA as a surrogate marker for a pathological activation of the NF-κB pathway. Three tumors showed alternative C11orf95 fusions to MAML2 or NCOA1. A methylation-based brain tumor classifier assigned two RELA-like tumors to the methylation class "EP, RELA-fusion"; the others demonstrated no significant similarity score. Of the tanycytic group, 5/6 tumors were assigned a WHO grade II. No gene fusions were detected. Methylation profiling did not show any association with an established methylation class. We additionally identified two astroblastoma-like tumors that both presented with chromothripsis of chromosome 22 but lacked MN1 breaks according to FISH analysis. They revealed novel fusion events involving genes in chromosome 22. One further tumor with polyploid cytogenetics was interpreted as PFB ependymoma by the brain tumor methylation classifier but had no relation to the posterior fossa. Clinical follow-up was available for 16/18 patients. Patients with tanycytic and astroblastoma-like tumors had no relapse, while 2 patients with RELA-like ependymomas died. Our data indicate that in addition to ependymomas discovered so far, at least two more supratentorial ependymoma types (RELA-like and tanycytic) exist.

Glioblastoma frequently exhibits therapy-associated subtype transitions to mesenchymal phenotypes with adverse prognosis. Here, we perform multi-omic profiling of 60 glioblastoma primary tumours and use orthogonal analysis of chromatin and RNA-derived gene regulatory networks to identify 38 subtype master regulators, whose cell population-specific activities we further map in published single-cell RNA sequencing data. These analyses identify the oligodendrocyte precursor marker and chromatin modifier SOX10 as a master regulator in RTK I-subtype tumours. In vitro functional studies demonstrate that SOX10 loss causes a subtype switch analogous to the proneural-mesenchymal transition observed in patients at the transcriptomic, epigenetic and phenotypic levels. SOX10 repression in an in vivo syngeneic graft glioblastoma mouse model results in increased tumour invasion, immune cell infiltration and significantly reduced survival, reminiscent of progressive human glioblastoma. These results identify SOX10 as a bona fide master regulator of the RTK I subtype, with both tumour cell-intrinsic and microenvironmental effects.

The maternal mode of mitochondrial DNA (mtDNA) inheritance is central to human genetics. Recently, evidence for bi-parental inheritance of mtDNA was claimed for individuals of three pedigrees that suffered mitochondrial disorders. We sequenced mtDNA using both direct Sanger and Massively Parallel Sequencing in several tissues of eleven maternally related and other affiliated healthy individuals of a family pedigree and observed mixed mitotypes in eight individuals. Cells without nuclear DNA, i.e. thrombocytes and hair shafts, only showed the mitotype of haplogroup (hg) V. Skin biopsies were prepared to generate ρ° cells void of mtDNA, sequencing of which resulted in a hg U4c1 mitotype. The position of the Mega-NUMT sequence was determined by fluorescence in situ hybridization and two different quantitative PCR assays were used to determine the number of contributing mtDNA copies. Thus, evidence for the presence of repetitive, full mitogenome Mega-NUMTs matching haplogroup U4c1 in various tissues of eight maternally related individuals was provided. Multi-copy Mega-NUMTs mimic mixtures of mtDNA that cannot be experimentally avoided and thus may appear in diverse fields of mtDNA research and diagnostics. We demonstrate that hair shaft mtDNA sequencing provides a simple but reliable approach to exclude NUMTs as source of misleading results.

Medulloblastoma is one of the most common malignant brain tumor types in children, with an overall survival of 70%. Mortality is associated with metastatic relapsed tumors. Rho-associated kinases (ROCKs), important for epithelial-mesenchymal transition (EMT) and proper nervous system development, have previously been identified as a promising drug target to inhibit cancer growth and metastatic spread. Here, we show that ROCKs are expressed in medulloblastoma, with higher ROCK2 mRNA expression in metastatic compared to non-metastatic tumors. By evaluating three ROCK inhibitors in a panel of medulloblastoma cell lines we demonstrated that medulloblastoma cells were sensitive for pharmacological ROCK inhibition. The specific ROCK inhibitor RKI-1447 inhibited the tumorigenicity in medulloblastoma cells as well as impeded cell migration and invasion. Differential gene expression analysis suggested that ROCK inhibition was associated with the downregulation of signaling pathways important in proliferation and metastasis e.g., TNFα via NFκβ, TGFβ, and EMT. Expression of key proteins in these pathways such as RHOA, RHOB, JUN, and vimentin was downregulated in ROCK inhibited cells. Finally, we showed that ROCK inhibition by RKI-1447 suppressed medulloblastoma growth and proliferation in vivo. Collectively, our results suggest that ROCK inhibition presents a potential new therapeutic option in medulloblastoma, especially for children with metastatic disease.

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-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.

Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium, for which whole-genome and-for a subset-whole-transcriptome sequencing data from 2,658 cancers across 38 tumor types was aggregated, we systematically investigated potential viral pathogens using a consensus approach that integrated three independent pipelines. Viruses were detected in 382 genome and 68 transcriptome datasets. We found a high prevalence of known tumor-associated viruses such as Epstein-Barr virus (EBV), hepatitis B virus (HBV) and human papilloma virus (HPV; for example, HPV16 or HPV18). The study revealed significant exclusivity of HPV and driver mutations in head-and-neck cancer and the association of HPV with APOBEC mutational signatures, which suggests that impaired antiviral defense is a driving force in cervical, bladder and head-and-neck carcinoma. For HBV, HPV16, HPV18 and adeno-associated virus-2 (AAV2), viral integration was associated with local variations in genomic copy numbers. Integrations at the TERT promoter were associated with high telomerase expression evidently activating this tumor-driving process. High levels of endogenous retrovirus (ERV1) expression were linked to a worse survival outcome in patients with kidney cancer.

Medulloblastoma is the most common malignant brain tumor of childhood. Group 3 medulloblastoma, the most aggressive molecular subtype, frequently disseminates through the leptomeningeal cerebral spinal fluid (CSF) spaces in the brain and spinal cord. The mechanism of dissemination through the CSF remains poorly understood, and the molecular pathways involved in medulloblastoma metastasis and self-renewal are largely unknown. Here we show that NOTCH1 signaling pathway regulates both the initiation of metastasis and the self-renewal of medulloblastoma. We identify a mechanism in which NOTCH1 activates BMI1 through the activation of TWIST1. NOTCH1 expression and activity are directly related to medulloblastoma metastasis and decreased survival rate of tumor-bearing mice. Finally, medulloblastoma-bearing mice intrathecally treated with anti-NRR1, a NOTCH1 blocking antibody, present lower frequency of spinal metastasis and higher survival rate. These findings identify NOTCH1 as a pivotal driver of Group 3 medulloblastoma metastasis and self-renewal, supporting the development of therapies targeting this pathway.

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.

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.

Genomic characterization has begun to redefine diagnostic classifications of cancers. However, it remains a challenge to infer disease phenotypes from genomic alterations alone. To help realize the promise of genomics, we have performed a quantitative proteomics investigation using Stable Isotope Labeling by Amino Acids in Cell Culture (SILAC) and 41 tissue samples spanning the 4 genomically based subgroups of medulloblastoma and control cerebellum. We have identified and quantitated thousands of proteins across these groups and find that we are able to recapitulate the genomic subgroups based upon subgroup restricted and differentially abundant proteins while also identifying subgroup specific protein isoforms. Integrating our proteomic measurements with genomic data, we calculate a poor correlation between mRNA and protein abundance. Using EPIC 850 k methylation array data on the same tissues, we also investigate the influence of copy number alterations and DNA methylation on the proteome in an attempt to characterize the impact of these genetic features on the proteome. Reciprocally, we are able to use the proteome to identify which genomic alterations result in altered protein abundance and thus are most likely to impact biology. Finally, we are able to assemble protein-based pathways yielding potential avenues for clinical intervention. From these, we validate the EIF4F cap-dependent translation pathway as a novel druggable pathway in medulloblastoma. Thus, quantitative proteomics complements genomic platforms to yield a more complete understanding of functional tumor biology and identify novel therapeutic targets for medulloblastoma.

Human neural stem cell cultures provide progenitor cells that are potential cells of origin for brain cancers. However, the extent to which genetic predisposition to tumor formation can be faithfully captured in stem cell lines is uncertain. Here, we evaluated neuroepithelial stem (NES) cells, representative of cerebellar progenitors. We transduced NES cells with MYCN, observing medulloblastoma upon orthotopic implantation in mice. Significantly, transcriptomes and patterns of DNA methylation from xenograft tumors were globally more representative of human medulloblastoma compared to a MYCN-driven genetically engineered mouse model. Orthotopic transplantation of NES cells generated from Gorlin syndrome patients, who are predisposed to medulloblastoma due to germline-mutated PTCH1, also generated medulloblastoma. We engineered candidate cooperating mutations in Gorlin NES cells, with mutation of DDX3X or loss of GSE1 both accelerating tumorigenesis. These findings demonstrate that human NES cells provide a potent experimental resource for dissecting genetic causation in medulloblastoma.

Extensive molecular analyses of ependymal tumors have revealed that supratentorial and posterior fossa ependymomas have distinct molecular profiles and are likely to be different diseases. The presence of C11orf95-RELA fusion genes in a subset of supratentorial ependymomas (ST-EPN) indicated the existence of molecular subgroups. However, the pathogenesis of RELA fusion-negative ependymomas remains elusive. To investigate the molecular pathogenesis of these tumors and validate the molecular classification of ependymal tumors, we conducted thorough molecular analyses of 113 locally diagnosed ependymal tumors from 107 patients in the Japan Pediatric Molecular Neuro-Oncology Group. All tumors were histopathologically reviewed and 12 tumors were re-classified as non-ependymomas. A combination of RT-PCR, FISH, and RNA sequencing identified RELA fusion in 19 of 29 histologically verified ST-EPN cases, whereas another case was diagnosed as ependymoma RELA fusion-positive via the methylation classifier (68.9%). Among the 9 RELA fusion-negative ST-EPN cases, either the YAP1 fusion, BCOR tandem duplication, EP300-BCORL1 fusion, or FOXO1-STK24 fusion was detected in single cases. Methylation classification did not identify a consistent molecular class within this group. Genome-wide methylation profiling successfully sub-classified posterior fossa ependymoma (PF-EPN) into PF-EPN-A (PFA) and PF-EPN-B (PFB). A multivariate analysis using Cox regression confirmed that PFA was the sole molecular marker which was independently associated with patient survival. A clinically applicable pyrosequencing assay was developed to determine the PFB subgroup with 100% specificity using the methylation status of 3 genes, CRIP1, DRD4 and LBX2. Our results emphasized the significance of molecular classification in the diagnosis of ependymomas. RELA fusion-negative ST-EPN appear to be a heterogeneous group of tumors that do not fall into any of the existing molecular subgroups and are unlikely to form a single category.