Accurate pathological diagnosis is crucial for optimal management of patients with cancer. For the approximately 100 known tumour types of the central nervous system, standardization of the diagnostic process has been shown to be particularly challenging-with substantial inter-observer variability in the histopathological diagnosis of many tumour types. Here we present a comprehensive approach for the DNA methylation-based classification of central nervous system tumours across all entities and age groups, and demonstrate its application in a routine diagnostic setting. We show that the availability of this method may have a substantial impact on diagnostic precision compared to standard methods, resulting in a change of diagnosis in up to 12% of prospective cases. For broader accessibility, we have designed a free online classifier tool, the use of which does not require any additional onsite data processing. Our results provide a blueprint for the generation of machine-learning-based tumour classifiers across other cancer entities, with the potential to fundamentally transform tumour pathology.

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.

MYCN is a master regulator controlling many processes necessary for tumor cell survival. Here, we unravel a microRNA network that causes tumor suppressive effects in MYCN-amplified neuroblastoma cells. In profiling studies, histone deacetylase (HDAC) inhibitor treatment most strongly induced miR-183. Enforced miR-183 expression triggered apoptosis, and inhibited anchorage-independent colony formation in vitro and xenograft growth in mice. Furthermore, the mechanism of miR-183 induction was found to contribute to the cell death phenotype induced by HDAC inhibitors. Experiments to identify the HDAC(s) involved in miR-183 transcriptional regulation showed that HDAC2 depletion induced miR-183. HDAC2 overexpression reduced miR-183 levels and counteracted the induction caused by HDAC2 depletion or HDAC inhibitor treatment. MYCN was found to recruit HDAC2 in the same complexes to the miR-183 promoter, and HDAC2 depletion enhanced promoter-associated histone H4 pan-acetylation, suggesting epigenetic changes preceded transcriptional activation. These data reveal miR-183 tumor suppressive properties in neuroblastoma that are jointly repressed by MYCN and HDAC2, and suggest a novel way to bypass MYCN function.

Pilocytic astrocytoma (PA) is the most frequent pediatric brain tumor. Activation of the MAPK pathway is well established as the oncogenic driver of the disease. It is most frequently caused by KIAA1549:BRAF fusions, and leads to oncogene induced senescence (OIS). OIS is thought to be a major reason for growth arrest of PA cells in vitro and in vivo, preventing establishment of PA cultures. Hence, valid preclinical models are currently very limited, but preclinical testing of new compounds is urgently needed. We transduced the PA short-term culture DKFZ-BT66 derived from the PA of a 2-year old patient with a doxycycline-inducible system coding for Simian Vacuolating Virus 40 Large T Antigen (SV40-TAg). SV40-TAg inhibits TP53/CDKN1A and CDKN2A/RB1, two pathways critical for OIS induction and maintenance. DNA methylation array and KIAA1549:BRAF fusion analysis confirmed pilocytic astrocytoma identity of DKFZ-BT66 cells after establishment. Readouts were analyzed in proliferating as well as senescent states, including cell counts, viability, cell cycle analysis, expression of SV40-Tag, CDKN2A (p16), CDKN1A (p21), and TP53 (p53) protein, and gene-expression profiling. Selected MAPK inhibitors (MAPKi) including clinically available MEK inhibitors (MEKi) were tested in vitro. Expression of SV40-TAg enabled the cells to bypass OIS and to resume proliferation with a mean doubling time of 45h allowing for propagation and long-term culture. Withdrawal of doxycycline led to an immediate decrease of SV40-TAg expression, appearance of senescent morphology, upregulation of CDKI proteins and a subsequent G1 growth arrest in line with the re-induction of senescence. DKFZ-BT66 cells still underwent replicative senescence that was overcome by TERT expression. Testing of a set of MAPKi revealed differential responses in DKFZ-BT66. MEKi efficiently inhibited MAPK signaling at clinically achievable concentrations, while BRAF V600E- and RAF Type II inhibitors showed paradoxical activation. Taken together, we have established the first patient-derived long term expandable PA cell line expressing the KIAA1549:BRAF-fusion suitable for preclinical drug testing.

A hallmark of pediatric low-grade glioma (pLGG) is aberrant signaling of the mitogen activated protein kinase (MAPK) pathway. Hence, inhibition of MAPK signaling using small molecule inhibitors such as MEK inhibitors (MEKi) may be a promising strategy.

Molecular groups of medulloblastoma (MB) are well established. Novel risk stratification parameters include Group 3/4 (non-WNT/non-SHH) methylation subgroups I-VIII or whole-chromosomal aberration (WCA) phenotypes. This study investigates the integration of clinical and molecular parameters to improve risk stratification of non-WNT/non-SHH MB. Non-WNT/non-SHH MB from the HIT2000 study and the HIT-MED registries were selected based on availability of DNA-methylation profiling data. MYC or MYCN amplification and WCA of chromosomes 7, 8, and 11 were inferred from methylation array-based copy number profiles. In total, 403 non-WNT/non-SHH MB were identified, 346/403 (86%) had a methylation class family Group 3/4 methylation score (classifier v11b6) ≥ 0.9, and 294/346 (73%) were included in the risk stratification modeling based on Group 3 or 4 score (v11b6) ≥ 0.8 and subgroup I-VIII score (mb_g34) ≥ 0.8. Group 3 MB (5y-PFS, survival estimation ± standard deviation: 41.4 ± 4.6%; 5y-OS: 48.8 ± 5.0%) showed poorer survival compared to Group 4 (5y-PFS: 68.2 ± 3.7%; 5y-OS: 84.8 ± 2.8%). Subgroups II (5y-PFS: 27.6 ± 8.2%) and III (5y-PFS: 37.5 ± 7.9%) showed the poorest and subgroup VI (5y-PFS: 76.6 ± 7.9%), VII (5y-PFS: 75.9 ± 7.2%), and VIII (5y-PFS: 66.6 ± 5.8%) the best survival. Multivariate analysis revealed subgroup in combination with WCA phenotype to best predict risk of progression and death. The integration of clinical (age, M and R status) and molecular (MYC/N, subgroup, WCA phenotype) variables identified a low-risk stratum with a 5y-PFS of 94 ± 5.7 and a very high-risk stratum with a 5y-PFS of 29 ± 6.1%. Validation in an international MB cohort confirmed the combined stratification scheme with 82.1 ± 6.0% 5y-PFS in the low and 47.5 ± 4.1% in very high-risk groups, and outperformed the clinical model. These newly identified clinico-molecular low-risk and very high-risk strata, accounting for 6%, and 21% of non-WNT/non-SHH MB patients, respectively, may improve future treatment stratification.

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.

Long-term complications such as radiation-induced second malignancies occur in a subset of patients following radiation-therapy, particularly relevant in pediatric patients due to the long follow-up period in case of survival. Radiation-induced gliomas (RIGs) have been reported in patients after treatment with cranial irradiation for various primary malignancies such as acute lymphoblastic leukemia (ALL) and medulloblastoma (MB). We perform comprehensive (epi-) genetic and expression profiling of RIGs arising after cranial irradiation for MB (n = 23) and ALL (n = 9). Our study reveals a unifying molecular signature for the majority of RIGs, with recurrent PDGFRA amplification and loss of CDKN2A/B and an absence of somatic hotspot mutations in genes encoding histone 3 variants or IDH1/2, uncovering diagnostic markers and potentially actionable targets.

Large-scale molecular profiling studies in recent years have shown that central nervous system (CNS) tumors display a much greater heterogeneity in terms of molecularly distinct entities, cellular origins and genetic drivers than anticipated from histological assessment. DNA methylation profiling has emerged as a useful tool for robust tumor classification, providing new insights into these heterogeneous molecular classes. This is particularly true for rare CNS tumors with a broad morphological spectrum, which are not possible to assign as separate entities based on histological similarity alone. Here, we describe a molecularly distinct subset of predominantly pediatric CNS neoplasms (n = 60) that harbor PATZ1 fusions. The original histological diagnoses of these tumors covered a wide spectrum of tumor types and malignancy grades. While the single most common diagnosis was glioblastoma (GBM), clinical data of the PATZ1-fused tumors showed a better prognosis than typical GBM, despite frequent relapses. RNA sequencing revealed recurrent MN1:PATZ1 or EWSR1:PATZ1 fusions related to (often extensive) copy number variations on chromosome 22, where PATZ1 and the two fusion partners are located. These fusions have individually been reported in a number of glial/glioneuronal tumors, as well as extracranial sarcomas. We show here that they are more common than previously acknowledged, and together define a biologically distinct CNS tumor type with high expression of neural development markers such as PAX2, GATA2 and IGF2. Drug screening performed on the MN1:PATZ1 fusion-bearing KS-1 brain tumor cell line revealed preliminary candidates for further study. In summary, PATZ1 fusions define a molecular class of histologically polyphenotypic neuroepithelial tumors, which show an intermediate prognosis under current treatment regimens.

MYC-driven Group 3 medulloblastoma (MB) (subtype II) is a highly aggressive childhood brain tumor. Sensitivity of MYC-driven MB to class I histone deacetylase inhibitors (HDACi) has been previously demonstrated in vitro and in vivo. In this study we characterize the transcriptional effects of class I HDACi in MYC-driven MB and explore beneficial drug combinations.

Patients with MYC-amplified Group 3 medulloblastoma (MB) (subtype II) show poor progression-free survival rates. Class I histone deacetylase inhibitors (HDACi) are highly effective for the treatment of MYC-amplified MB in vitro and in vivo. Drug combination regimens including class I HDACi may represent an urgently needed novel treatment approach for this high risk disease.

Pilocytic astrocytoma is the most common childhood brain tumor, characterized by constitutive MAPK activation. MAPK signaling induces oncogene-induced senescence (OIS), which may cause unpredictable growth behavior of pilocytic astrocytomas. The senescence-associated secretory phenotype (SASP) has been shown to regulate OIS, but its role in pilocytic astrocytoma remains unknown.Experimental Design: The patient-derived pilocytic astrocytoma cell culture model, DKFZ-BT66, was used to demonstrate presence of the SASP and analyze its impact on OIS in pilocytic astrocytoma. The model allows for doxycycline-inducible switching between proliferation and OIS. Both states were studied using gene expression profiling (GEP), Western blot, ELISA, and cell viability testing. Primary pilocytic astrocytoma tumors were analyzed by GEP and multiplex assay.

Ependymomas in children can arise throughout all compartments of the central nervous system (CNS). Highly malignant paediatric ependymoma subtypes are Group A tumours of the posterior fossa (PF-EPN-A) and RELA-fusion positive (ST-EPN-RELA) tumours in the supratentorial compartment. It was repeatedly reported in smaller series that accumulation of p53 is frequently observed in ependymomas and that immunohistochemical staining correlates with poor clinical outcome, while TP53 mutations are rare. Our TP53 mutation analysis of 130 primary ependymomas identified a mutation rate of only 3%. Immunohistochemical analysis of 398 ependymomas confirmed previous results correlating the accumulation of p53 with inferior outcome. Among the p53-positive ependymomas, the vast majority exhibited a RELA fusion leading to the hypothesis that p53 inactivation might be linked to RELA positivity.In order to assess the potential of p53 reactivation through MDM2 inhibition in ependymoma, we evaluated the effects of Actinomycin-D and Nutlin-3 treatment in two preclinical ependymoma models representing the high-risk subtypes PF-EPN-A and ST-EPN-RELA. The IC-50 of the agent as determined by metabolic activity assays was in the lower nano-molar range (0.2-0.7 nM). Transcriptome analyses of high-dose (100 nM), low-dose (5 nM) and non-treated cells revealed re-expression of p53 dependent genes including p53 upregulated modulator of apoptosis (PUMA) after low-dose treatment. At the protein level, we validated the Actinomycin-D induced upregulation of PUMA, and of p53 interaction partners MDM2 and p21. Proapoptotic effects of low-dose application of the agent were confirmed by flow cytometry. Thus, Actinomycin-D could constitute a promising therapeutic option for ST-EPN-RELA ependymoma patients, whose tumours frequently exhibit p53 inactivation.

Current preclinical models in tumor biology are limited in their ability to recapitulate relevant (patho-) physiological processes, including autophagy. Three-dimensional (3D) growth cultures have frequently been proposed to overcome the lack of correlation between two-dimensional (2D) monolayer cell cultures and human tumors in preclinical drug testing. Besides 3D growth, it is also advantageous to simulate shear stress, compound flux and removal of metabolites, e.g., via bioreactor systems, through which culture medium is constantly pumped at a flow rate reflecting physiological conditions. Here we show that both static 3D growth and 3D growth within a bioreactor system modulate key hallmarks of cancer cells, including proliferation and cell death as well as macroautophagy, a recycling pathway often activated by highly proliferative tumors to cope with metabolic stress. The autophagy-related gene expression profiles of 2D-grown cells are substantially different from those of 3D-grown cells and tumor tissue. Autophagy-controlling transcription factors, such as TFEB and FOXO3, are upregulated in tumors, and 3D-grown cells have increased expression compared with cells grown in 2D conditions. Three-dimensional cultures depleted of the autophagy mediators BECN1, ATG5 or ATG7 or the transcription factor FOXO3, are more sensitive to cytotoxic treatment. Accordingly, combining cytotoxic treatment with compounds affecting late autophagic flux, such as chloroquine, renders the 3D-grown cells more susceptible to therapy. Altogether, 3D cultures are a valuable tool to study drug response of tumor cells, as these models more closely mimic tumor (patho-)physiology, including the upregulation of tumor relevant pathways, such as autophagy.

Pleomorphic xanthoastrocytoma (PXA) in its classic manifestation exhibits distinct morphological features and is assigned to CNS WHO grade 2 or grade 3. Distinction from glioblastoma variants and lower grade glial and glioneuronal tumors is a common diagnostic challenge. We compared a morphologically defined set of PXA (histPXA) with an independent set, defined by DNA methylation analysis (mcPXA). HistPXA encompassed 144 tumors all subjected to DNA methylation array analysis. Sixty-two histPXA matched to the methylation class mcPXA. These were combined with the cases that showed the mcPXA signature but had received a histopathological diagnosis other than PXA. This cohort constituted a set of 220 mcPXA. Molecular and clinical parameters were analyzed in these groups. Morphological parameters were analyzed in a subset of tumors with FFPE tissue available. HistPXA revealed considerable heterogeneity in regard to methylation classes, with methylation classes glioblastoma and ganglioglioma being the most frequent mismatches. Similarly, the mcPXA cohort contained tumors of diverse histological diagnoses, with glioblastoma constituting the most frequent mismatch. Subsequent analyses demonstrated the presence of canonical pTERT mutations to be associated with unfavorable prognosis among mcPXA. Based on these data, we consider the tumor type PXA to be histologically more varied than previously assumed. Histological approach to diagnosis will predominantly identify cases with the established archetypical morphology. DNA methylation analysis includes additional tumors in the tumor class PXA that share similar DNA methylation profile but lack the typical morphology of a PXA. DNA methylation analysis also assist in separating other tumor types with morphologic overlap to PXA. Our data suggest the presence of canonical pTERT mutations as a robust indicator for poor prognosis in methylation class PXA.

The survival rate among children with relapsed neuroblastomas continues to be poor, and thus new therapeutic approaches identified by reliable preclinical drug testing models are urgently needed. Zebrafish are a powerful vertebrate model in preclinical cancer research. Here, we describe a zebrafish neuroblastoma yolk sac model to evaluate efficacy and toxicity of histone deacetylase (HDAC) inhibitor treatments. Larvae were engrafted with fluorescently labeled, genetically diverse, established cell lines and short-term cultures of patient-derived primary cells. Engrafted tumors progressed locally and disseminated remotely in an intact environment. Combination treatments involving the standard chemotherapy doxorubicin and HDAC inhibitors substantially reduced tumor volume, induced tumor cell death, and inhibited tumor cell dissemination to the tail region. Hence, this model allows for fast, cost-efficient, and reliable in vivo evaluation of toxicity and response of the primary and metastatic tumor sites to drug combinations.

Glioneuronal tumors are a heterogenous group of CNS neoplasms that can be challenging to accurately diagnose. Molecular methods are highly useful in classifying these tumors-distinguishing precise classes from their histological mimics and identifying previously unrecognized types of tumors. Using an unsupervised visualization approach of DNA methylation data, we identified a novel group of tumors (n = 20) that formed a cluster separate from all established CNS tumor types. Molecular analyses revealed ATRX alterations (in 16/16 cases by DNA sequencing and/or immunohistochemistry) as well as potentially targetable gene fusions involving receptor tyrosine-kinases (RTK; mostly NTRK1-3) in all of these tumors (16/16; 100%). In addition, copy number profiling showed homozygous deletions of CDKN2A/B in 55% of cases. Histological and immunohistochemical investigations revealed glioneuronal tumors with isomorphic, round and often condensed nuclei, perinuclear clearing, high mitotic activity and microvascular proliferation. Tumors were mainly located supratentorially (84%) and occurred in patients with a median age of 19 years. Survival data were limited (n = 18) but point towards a more aggressive biology as compared to other glioneuronal tumors (median progression-free survival 12.5 months). Given their molecular characteristics in addition to anaplastic features, we suggest the term glioneuronal tumor with ATRX alteration, kinase fusion and anaplastic features (GTAKA) to describe these tumors. In summary, our findings highlight a novel type of glioneuronal tumor driven by different RTK fusions accompanied by recurrent alterations in ATRX and homozygous deletions of CDKN2A/B. Targeted approaches such as NTRK inhibition might represent a therapeutic option for patients suffering from these tumors.

Pediatric low-grade gliomas (pLGG) show heterogeneous responses to MAPK inhibitors (MAPKi) in clinical trials. Thus, more complex stratification biomarkers are needed to identify patients likely to benefit from MAPKi therapy. Here, we identify MAPK-related genes enriched in MAPKi-sensitive cell lines using the GDSC dataset and apply them to calculate class-specific MAPKi sensitivity scores (MSSs) via single-sample gene set enrichment analysis. The MSSs discriminate MAPKi-sensitive and non-sensitive cells in the GDSC dataset and significantly correlate with response to MAPKi in an independent PDX dataset. The MSSs discern gliomas with varying MAPK alterations and are higher in pLGG compared to other pediatric CNS tumors. Heterogenous MSSs within pLGGs with the same MAPK alteration identify proportions of potentially sensitive patients. The MEKi MSS predicts treatment response in a small set of pLGG patients treated with trametinib. High MSSs correlate with a higher immune cell infiltration, with high expression in the microglia compartment in single-cell RNA sequencing data, while low MSSs correlate with low immune infiltration and increased neuronal score. The MSSs represent predictive tools for the stratification of pLGG patients and should be prospectively validated in clinical trials. Our data supports a role for microglia in the response to MAPKi.

Pilocytic astrocytoma (PA), the most common pediatric brain tumor, is driven by aberrant mitogen-activated protein kinase signaling most commonly caused by BRAF gene fusions or activating mutations. While 5-year overall survival rates exceed 95%, tumor recurrence or progression constitutes a major clinical challenge in incompletely resected tumors. Here, we used similarity network fusion (SNF) analysis in an integrative multi-omics approach employing RNA transcriptomic and mass spectrometry-based proteomic profiling to molecularly characterize PA tissue samples from 62 patients. Thereby, we uncovered that PAs segregated into two molecularly distinct groups, namely, Group 1 and Group 2, which were validated in three non-overlapping cohorts. Patients with Group 1 tumors were significantly younger and showed worse progression-free survival compared to patients with group 2 tumors. Ingenuity pathways analysis (IPA) and gene set enrichment analysis (GSEA) revealed that Group 1 tumors were enriched for immune response pathways, such as interferon signaling, while Group 2 tumors showed enrichment for action potential and neurotransmitter signaling pathways. Analysis of immune cell-related gene signatures showed an enrichment of infiltrating T Cells in Group 1 versus Group 2 tumors. Taken together, integrative multi-omics of PA identified biologically distinct and prognostically relevant tumor groups that may improve risk stratification of this single pathway driven tumor type.

Current therapies for medulloblastoma, a highly malignant childhood brain tumour, impose debilitating effects on the developing child, and highlight the need for molecularly targeted treatments with reduced toxicity. Previous studies have been unable to identify the full spectrum of driver genes and molecular processes that operate in medulloblastoma subgroups. Here we analyse the somatic landscape across 491 sequenced medulloblastoma samples and the molecular heterogeneity among 1,256 epigenetically analysed cases, and identify subgroup-specific driver alterations that include previously undiscovered actionable targets. Driver mutations were confidently assigned to most patients belonging to Group 3 and Group 4 medulloblastoma subgroups, greatly enhancing previous knowledge. New molecular subtypes were differentially enriched for specific driver events, including hotspot in-frame insertions that target KBTBD4 and 'enhancer hijacking' events that activate PRDM6. Thus, the application of integrative genomics to an extensive cohort of clinical samples derived from a single childhood cancer entity revealed a series of cancer genes and biologically relevant subtype diversity that represent attractive therapeutic targets for the treatment of patients with medulloblastoma.