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.

iTHER is a Dutch prospective national precision oncology program aiming to define tumour molecular profiles in children and adolescents with primary very high-risk, relapsed, or refractory paediatric tumours. Between April 2017 and April 2021, 302 samples from 253 patients were included. Comprehensive molecular profiling including low-coverage whole genome sequencing (lcWGS), whole exome sequencing (WES), RNA sequencing (RNA-seq), Affymetrix, and/or 850k methylation profiling was successfully performed for 226 samples with at least 20% tumour content. Germline pathogenic variants were identified in 16% of patients (35/219), of which 22 variants were judged causative for a cancer predisposition syndrome. At least one somatic alteration was detected in 204 (90.3%), and 185 (81.9%) were considered druggable, with clinical priority very high (6.1%), high (21.3%), moderate (26.0%), intermediate (36.1%), and borderline (10.5%) priority. iTHER led to revision or refinement of diagnosis in 8 patients (3.5%). Temporal heterogeneity was observed in paired samples of 15 patients, indicating the value of sequential analyses. Of 137 patients with follow-up beyond twelve months, 21 molecularly matched treatments were applied in 19 patients (13.9%), with clinical benefit in few. Most relevant barriers to not applying targeted therapies included poor performance status, as well as limited access to drugs within clinical trial. iTHER demonstrates the feasibility of comprehensive molecular profiling across all ages, tumour types and stages in paediatric cancers, informing of diagnostic, prognostic, and targetable alterations as well as reportable germline variants. Therefore, WES and RNA-seq is nowadays standard clinical care at the Princess Máxima Center for all children with cancer, including patients at primary diagnosis. Improved access to innovative treatments within biology-driven combination trials is required to ultimately improve survival.

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.

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.

YAP1 fusion-positive supratentorial ependymomas predominantly occur in infants, but the molecular mechanisms of oncogenesis are unknown. Here we show YAP1-MAMLD1 fusions are sufficient to drive malignant transformation in mice, and the resulting tumors share histo-molecular characteristics of human ependymomas. Nuclear localization of YAP1-MAMLD1 protein is mediated by MAMLD1 and independent of YAP1-Ser127 phosphorylation. Chromatin immunoprecipitation-sequencing analyses of human YAP1-MAMLD1-positive ependymoma reveal enrichment of NFI and TEAD transcription factor binding site motifs in YAP1-bound regulatory elements, suggesting a role for these transcription factors in YAP1-MAMLD1-driven tumorigenesis. Mutation of the TEAD binding site in the YAP1 fusion or repression of NFI targets prevents tumor induction in mice. Together, these results demonstrate that the YAP1-MAMLD1 fusion functions as an oncogenic driver of ependymoma through recruitment of TEADs and NFIs, indicating a rationale for preclinical studies to block the interaction between YAP1 fusions and NFI and TEAD transcription factors.

Extra-cranial malignant rhabdoid tumors (MRTs) and cranial atypical teratoid RTs (ATRTs) are heterogeneous pediatric cancers driven primarily by SMARCB1 loss. To understand the genome-wide molecular relationships between MRTs and ATRTs, we analyze multi-omics data from 140 MRTs and 161 ATRTs. We detect similarities between the MYC subgroup of ATRTs (ATRT-MYC) and extra-cranial MRTs, including global DNA hypomethylation and overexpression of HOX genes and genes involved in mesenchymal development, distinguishing them from other ATRT subgroups that express neural-like features. We identify five DNA methylation subgroups associated with anatomical sites and SMARCB1 mutation patterns. Groups 1, 3, and 4 exhibit cytotoxic T cell infiltration and expression of immune checkpoint regulators, consistent with a potential role for immunotherapy in rhabdoid tumor patients.

Telomere maintenance by telomerase activation or alternative lengthening of telomeres (ALT) is a major determinant of poor outcome in neuroblastoma. Here, we screen for ALT in primary and relapsed neuroblastomas (n = 760) and characterize its features using multi-omics profiling. ALT-positive tumors are molecularly distinct from other neuroblastoma subtypes and enriched in a population-based clinical sequencing study cohort for relapsed cases. They display reduced ATRX/DAXX complex abundance, due to either ATRX mutations (55%) or low protein expression. The heterochromatic histone mark H3K9me3 recognized by ATRX is enriched at the telomeres of ALT-positive tumors. Notably, we find a high frequency of telomeric repeat loci with a neuroblastoma ALT-specific hotspot on chr1q42.2 and loss of the adjacent chromosomal segment forming a neo-telomere. ALT-positive neuroblastomas proliferate slowly, which is reflected by a protracted clinical course of disease. Nevertheless, children with an ALT-positive neuroblastoma have dismal outcome.

Myxopapillary ependymoma (MPE) is a heterogeneous disease regarding histopathology and outcome. The underlying molecular biology is poorly understood, and markers that reliably predict the patients' clinical course are unknown.

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.

Ependymomas encompass a heterogeneous group of central nervous system (CNS) neoplasms that occur along the entire neuroaxis. In recent years, extensive (epi-)genomic profiling efforts have identified several molecular groups of ependymoma that are characterized by distinct molecular alterations and/or patterns. Based on unsupervised visualization of a large cohort of genome-wide DNA methylation data, we identified a highly distinct group of pediatric-type tumors (n = 40) forming a cluster separate from all established CNS tumor types, of which a high proportion were histopathologically diagnosed as ependymoma. RNA sequencing revealed recurrent fusions involving the pleomorphic adenoma gene-like 1 (PLAGL1) gene in 19 of 20 of the samples analyzed, with the most common fusion being EWSR1:PLAGL1 (n = 13). Five tumors showed a PLAGL1:FOXO1 fusion and one a PLAGL1:EP300 fusion. High transcript levels of PLAGL1 were noted in these tumors, with concurrent overexpression of the imprinted genes H19 and IGF2, which are regulated by PLAGL1. Histopathological review of cases with sufficient material (n = 16) demonstrated a broad morphological spectrum of tumors with predominant ependymoma-like features. Immunohistochemically, tumors were GFAP positive and OLIG2- and SOX10 negative. In 3/16 of the cases, a dot-like positivity for EMA was detected. All tumors in our series were located in the supratentorial compartment. Median age of the patients at the time of diagnosis was 6.2 years. Median progression-free survival was 35 months (for 11 patients with data available). In summary, our findings suggest the existence of a novel group of supratentorial neuroepithelial tumors that are characterized by recurrent PLAGL1 fusions and enriched for pediatric patients.

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.

Glioblastoma IDH-wildtype presents with a wide histological spectrum. Some features are so distinctive that they are considered as separate histological variants or patterns for the purpose of classification. However, these usually lack defined (epi-)genetic alterations or profiles correlating with this histology. Here, we describe a molecular subtype with overlap to the unique histological pattern of glioblastoma with primitive neuronal component. Our cohort consists of 63 IDH-wildtype glioblastomas that harbor a characteristic DNA methylation profile. Median age at diagnosis was 59.5 years. Copy-number variations and genetic sequencing revealed frequent alterations in TP53, RB1 and PTEN, with fewer gains of chromosome 7 and homozygous CDKN2A/B deletions than usually described for IDH-wildtype glioblastoma. Gains of chromosome 1 were detected in more than half of the cases. A poorly differentiated phenotype with frequent absence of GFAP expression, high proliferation index and strong staining for p53 and TTF1 often caused misleading histological classification as carcinoma metastasis or primitive neuroectodermal tumor. Clinically, many patients presented with leptomeningeal dissemination and spinal metastasis. Outcome was poor with a median overall survival of only 12 months. Overall, we describe a new molecular subtype of IDH-wildtype glioblastoma with a distinct histological appearance and genetic signature.

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.

As the progression of low-grade diffuse astrocytomas into grade 4 tumors significantly impacts patient prognosis, a better understanding of this process is of paramount importance for improved patient care. In this project, we analyzed matched IDH-mutant astrocytomas before and after progression to grade 4 from six patients (discovery cohort) with genome-wide sequencing, 21 additional patients with targeted sequencing, and 33 patients from Glioma Longitudinal AnalySiS cohort for validation. The Cancer Genome Atlas data from 595 diffuse gliomas provided supportive information. All patients in our discovery cohort received radiation, all but one underwent chemotherapy, and no patient received temozolomide (TMZ) before progression to grade 4 disease. One case in the discovery cohort exhibited a hypermutation signature associated with the inactivation of the MSH2 and DNMT3A genes. In other patients, the number of chromosomal rearrangements and deletions increased in grade 4 tumors. The cell cycle checkpoint gene CDKN2A, or less frequently RB1, was most commonly inactivated after receiving both chemo- and radiotherapy when compared to other treatment groups. Concomitant activating PDGFRA/MET alterations were detected in tumors that acquired a homozygous CDKN2A deletion. NRG3 gene was significantly downregulated and recurrently altered in progressed tumors. Its decreased expression was associated with poorer overall survival in both univariate and multivariate analysis. We also detected progression-related alterations in RAD51B and other DNA repair pathway genes associated with the promotion of error-prone DNA repair, potentially facilitating tumor progression. In our retrospective analysis of patient treatment and survival timelines (n = 75), the combination of postoperative radiation and chemotherapy (mainly TMZ) outperformed radiation, especially in the grade 3 tumor cohort, in which it was typically given after primary surgery. Our results provide further insight into the contribution of treatment and genetic alterations in cell cycle, growth factor signaling, and DNA repair-related genes to tumor evolution and progression.

The large diversity of central nervous system (CNS) tumor types in children and adolescents results in disparate patient outcomes and renders accurate diagnosis challenging. In this study, we prospectively integrated DNA methylation profiling and targeted gene panel sequencing with blinded neuropathological reference diagnostics for a population-based cohort of more than 1,200 newly diagnosed pediatric patients with CNS tumors, to assess their utility in routine neuropathology. We show that the multi-omic integration increased diagnostic accuracy in a substantial proportion of patients through annotation to a refining DNA methylation class (50%), detection of diagnostic or therapeutically relevant genetic alterations (47%) or identification of cancer predisposition syndromes (10%). Discrepant results by neuropathological WHO-based and DNA methylation-based classification (30%) were enriched in histological high-grade gliomas, implicating relevance for current clinical patient management in 5% of all patients. Follow-up (median 2.5 years) suggests improved survival for patients with histological high-grade gliomas displaying lower-grade molecular profiles. These results provide preliminary evidence of the utility of integrating multi-omics in neuropathology for pediatric neuro-oncology.

BRAF genomic alterations are the most common oncogenic drivers in pediatric low-grade glioma (pLGG). Arm 1 (n = 77) of the ongoing phase 2 FIREFLY-1 (PNOC026) trial investigated the efficacy of the oral, selective, central nervous system-penetrant, type II RAF inhibitor tovorafenib (420 mg m-2 once weekly; 600 mg maximum) in patients with BRAF-altered, relapsed/refractory pLGG. Arm 2 (n = 60) is an extension cohort, which provided treatment access for patients with RAF-altered pLGG after arm 1 closure. Based on independent review, according to Response Assessment in Neuro-Oncology High-Grade Glioma (RANO-HGG) criteria, the overall response rate (ORR) of 67% met the arm 1 prespecified primary endpoint; median duration of response (DOR) was 16.6 months; and median time to response (TTR) was 3.0 months (secondary endpoints). Other select arm 1 secondary endpoints included ORR, DOR and TTR as assessed by Response Assessment in Pediatric Neuro-Oncology Low-Grade Glioma (RAPNO) criteria and safety (assessed in all treated patients and the primary endpoint for arm 2, n = 137). The ORR according to RAPNO criteria (including minor responses) was 51%; median DOR was 13.8 months; and median TTR was 5.3 months. The most common treatment-related adverse events (TRAEs) were hair color changes (76%), elevated creatine phosphokinase (56%) and anemia (49%). Grade ≥3 TRAEs occurred in 42% of patients. Nine (7%) patients had TRAEs leading to discontinuation of tovorafenib. These data indicate that tovorafenib could be an effective therapy for BRAF-altered, relapsed/refractory pLGG. ClinicalTrials.gov registration: NCT04775485 .

Posterior fossa type A (PF-EPN-A, PFA) ependymoma are aggressive tumors that mainly affect children and have a poor prognosis. Histopathology shows significant intratumoral heterogeneity, ranging from loose tissue to often sharply demarcated, extremely cell-dense tumor areas. To determine molecular differences in morphologically different areas and to understand their clinical significance, we analyzed 113 PF-EPN-A samples, including 40 corresponding relapse samples. Cell-dense areas ranged from 0 to 100% of the tumor area and displayed a higher proportion of proliferating tumor cells (p < 0.01). Clinically, cell density was associated with poor progression-free and overall survival (pPFS = 0.0026, pOS < 0.01). Molecularly, tumor areas with low and high cell density showed diverging DNA methylation profiles regarding their similarity to distinct previously discovered PF-EPN-A subtypes in 9/21 cases. Prognostically relevant chromosomal changes at 1q and 6q showed spatial heterogeneity within single tumors and were significantly enriched in cell-dense tumor areas as shown by single-cell RNA (scRNA)-sequencing as well as copy number profiling and fluorescence in situ hybridization (FISH) analyses of different tumor areas. Finally, spatial transcriptomics revealed cell-dense areas of different tumors to be more similar than various different areas of the same tumor. High-density areas distinctly overexpressed genes encoding histone proteins, WNT5A, TGFB1, or IGF2. Relapsing tumors displayed a higher proportion of cell-dense areas (p = 0.036), a change in PF-EPN-A methylation subtypes (13/32 patients), and novel chromosome 1q gains and 6q losses (12/32 cases) compared to corresponding primary tumors. Our data suggest that PF-EPN-A ependymomas habor a previously unrecognized intratumoral heterogeneity with clinical implications, which has to be accounted for when selecting diagnostic material, inter alia, by histological evaluation of the proportion of cell-dense areas.

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.

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.