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.

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.

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.

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.

Both the perivascular niche (PVN) and the integration into multicellular networks by tumor microtubes (TMs) have been associated with progression and resistance to therapies in glioblastoma, but their specific contribution remained unknown. By long-term tracking of tumor cell fate and dynamics in the live mouse brain, differential therapeutic responses in both niches are determined. Both the PVN, a preferential location of long-term quiescent glioma cells, and network integration facilitate resistance against cytotoxic effects of radiotherapy and chemotherapy-independently of each other, but with additive effects. Perivascular glioblastoma cells are particularly able to actively repair damage to tumor regions. Population of the PVN and resistance in it depend on proficient NOTCH1 expression. In turn, NOTCH1 downregulation induces resistant multicellular networks by TM extension. Our findings identify NOTCH1 as a central switch between the PVN and network niche in glioma, and demonstrate robust cross-compensation when only one niche is targeted.

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.

In vivo functional investigation of oncogenes using somatic gene transfer has been successfully exploited to validate their role in tumorigenesis. For tumour suppressor genes this has proven more challenging due to technical aspects. To provide a flexible and effective method for investigating somatic loss-of-function alterations and their influence on tumorigenesis, we have established CRISPR/Cas9-mediated somatic gene disruption, allowing for in vivo targeting of TSGs. Here we demonstrate the utility of this approach by deleting single (Ptch1) or multiple genes (Trp53, Pten, Nf1) in the mouse brain, resulting in the development of medulloblastoma and glioblastoma, respectively. Using whole-genome sequencing (WGS) we characterized the medulloblastoma-driving Ptch1 deletions in detail and show that no off-targets were detected in these tumours. This method provides a fast and convenient system for validating the emerging wealth of novel candidate tumour suppressor genes and the generation of faithful animal models of human cancer.

Intrinsic malignant brain tumors, such as glioblastomas are frequently resistant to immune checkpoint blockade (ICB) with few hypermutated glioblastomas showing response. Modeling patient-individual resistance is challenging due to the lack of predictive biomarkers and limited accessibility of tissue for serial biopsies. Here, we investigate resistance mechanisms to anti-PD-1 and anti-CTLA-4 therapy in syngeneic hypermutated experimental gliomas and show a clear dichotomy and acquired immune heterogeneity in ICB-responder and non-responder tumors. We made use of this dichotomy to establish a radiomic signature predicting tumor regression after pseudoprogression induced by ICB therapy based on serial magnetic resonance imaging. We provide evidence that macrophage-driven ICB resistance is established by CD4 T cell suppression and Treg expansion in the tumor microenvironment via the PD-L1/PD-1/CD80 axis. These findings uncover an unexpected heterogeneity of response to ICB in strictly syngeneic tumors and provide a rationale for targeting PD-L1-expressing tumor-associated macrophages to overcome resistance to ICB.

Mutations in chromatin modifier genes are frequently associated with neurodevelopmental diseases. We herein demonstrate that the chromodomain helicase DNA-binding protein 7 (Chd7), frequently associated with CHARGE syndrome, is indispensable for normal cerebellar development. Genetic inactivation of Chd7 in cerebellar granule neuron progenitors leads to cerebellar hypoplasia in mice, due to the impairment of granule neuron differentiation, induction of apoptosis and abnormal localization of Purkinje cells, which closely recapitulates known clinical features in the cerebella of CHARGE patients. Combinatory molecular analyses reveal that Chd7 is required for the maintenance of open chromatin and thus activation of genes essential for granule neuron differentiation. We further demonstrate that both Chd7 and Top2b are necessary for the transcription of a set of long neuronal genes in cerebellar granule neurons. Altogether, our comprehensive analyses reveal a mechanism with chromatin remodellers governing brain development via controlling a core transcriptional programme for cell-specific differentiation.

Although machine learning (ML) has shown promise across disciplines, out-of-sample generalizability is concerning. This is currently addressed by sharing multi-site data, but such centralization is challenging/infeasible to scale due to various limitations. Federated ML (FL) provides an alternative paradigm for accurate and generalizable ML, by only sharing numerical model updates. Here we present the largest FL study to-date, involving data from 71 sites across 6 continents, to generate an automatic tumor boundary detector for the rare disease of glioblastoma, reporting the largest such dataset in the literature (n = 6, 314). We demonstrate a 33% delineation improvement for the surgically targetable tumor, and 23% for the complete tumor extent, over a publicly trained model. We anticipate our study to: 1) enable more healthcare studies informed by large diverse data, ensuring meaningful results for rare diseases and underrepresented populations, 2) facilitate further analyses for glioblastoma by releasing our consensus model, and 3) demonstrate the FL effectiveness at such scale and task-complexity as a paradigm shift for multi-site collaborations, alleviating the need for data-sharing.

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.

DICER1 syndrome is a tumor predisposition syndrome that is associated with up to 30 different neoplastic lesions, usually affecting children and adolescents. Here we identify a group of mesenchymal tumors which is highly associated with DICER1 syndrome, and molecularly distinct from other DICER1-associated tumors. This group of DICER1-associated mesenchymal tumors encompasses multiple well-established clinicopathological tumor entities and can be further divided into three clinically meaningful classes designated "low-grade mesenchymal tumor with DICER1 alteration" (LGMT DICER1), "sarcoma with DICER1 alteration" (SARC DICER1), and primary intracranial sarcoma with DICER1 alteration (PIS DICER1). Our study not only provides a combined approach to classify DICER1-associated neoplasms for improved clinical management but also suggests a role for global hypomethylation and other recurrent molecular events in sarcomatous differentiation in mesenchymal tumors with DICER1 alteration. Our results will facilitate future investigations into prognostication and therapeutic approaches for affected patients.

Medulloblastomas with extensive nodularity are cerebellar tumors characterized by two distinct compartments and variable disease progression. The mechanisms governing the balance between proliferation and differentiation in MBEN remain poorly understood. Here, we employ a multi-modal single cell transcriptome analysis to dissect this process. In the internodular compartment, we identify proliferating cerebellar granular neuronal precursor-like malignant cells, along with stromal, vascular, and immune cells. In contrast, the nodular compartment comprises postmitotic, neuronally differentiated malignant cells. Both compartments are connected through an intermediate cell stage resembling actively migrating CGNPs. Notably, we also discover astrocytic-like malignant cells, found in proximity to migrating and differentiated cells at the transition zone between the two compartments. Our study sheds light on the spatial tissue organization and its link to the developmental trajectory, resulting in a more benign tumor phenotype. This integrative approach holds promise to explore intercompartmental interactions in other cancers with varying histology.

Many genes that drive normal cellular development also contribute to oncogenesis. Medulloblastoma (MB) tumors likely arise from neuronal progenitors in the cerebellum, and we hypothesized that the heterogeneity observed in MBs with sonic hedgehog (SHH) activation could be due to differences in developmental pathways. To investigate this question, here we perform single-nucleus RNA sequencing on highly differentiated SHH MBs with extensively nodular histology and observed malignant cells resembling each stage of canonical granule neuron development. Through innovative computational approaches, we connect these results to published datasets and find that some established molecular subtypes of SHH MB appear arrested at different developmental stages. Additionally, using multiplexed proteomic imaging and MALDI imaging mass spectrometry, we identify distinct histological and metabolic profiles for highly differentiated tumors. Our approaches are applicable to understanding the interplay between heterogeneity and differentiation in other cancers and can provide important insights for the design of targeted therapies.

Tumor microtubes (TMs) connect glioma cells to a network with considerable relevance for tumor progression and therapy resistance. However, the determination of TM-interconnectivity in individual tumors is challenging and the impact on patient survival unresolved. Here, we establish a connectivity signature from single-cell RNA-sequenced (scRNA-Seq) xenografted primary glioblastoma (GB) cells using a dye uptake methodology, and validate it with recording of cellular calcium epochs and clinical correlations. Astrocyte-like and mesenchymal-like GB cells have the highest connectivity signature scores in scRNA-sequenced patient-derived xenografts and patient samples. In large GB cohorts, TM-network connectivity correlates with the mesenchymal subtype and dismal patient survival. CHI3L1 gene expression serves as a robust molecular marker of connectivity and functionally influences TM networks. The connectivity signature allows insights into brain tumor biology, provides a proof-of-principle that tumor cell TM-connectivity is relevant for patients' prognosis, and serves as a robust prognostic biomarker.

As whole-genome sequencing for cancer genome analysis becomes a clinical tool, a full understanding of the variables affecting sequencing analysis output is required. Here using tumour-normal sample pairs from two different types of cancer, chronic lymphocytic leukaemia and medulloblastoma, we conduct a benchmarking exercise within the context of the International Cancer Genome Consortium. We compare sequencing methods, analysis pipelines and validation methods. We show that using PCR-free methods and increasing sequencing depth to ∼ 100 × shows benefits, as long as the tumour:control coverage ratio remains balanced. We observe widely varying mutation call rates and low concordance among analysis pipelines, reflecting the artefact-prone nature of the raw data and lack of standards for dealing with the artefacts. However, we show that, using the benchmark mutation set we have created, many issues are in fact easy to remedy and have an immediate positive impact on mutation detection accuracy.

Chromothripsis and chromoanasynthesis are catastrophic events leading to clustered genomic rearrangements. Whole-genome sequencing revealed frequent complex genomic rearrangements (n = 16/26) in brain tumors developing in mice deficient for factors involved in homologous-recombination-repair or non-homologous-end-joining. Catastrophic events were tightly linked to Myc/Mycn amplification, with increased DNA damage and inefficient apoptotic response already observable at early postnatal stages. Inhibition of repair processes and comparison of the mouse tumors with human medulloblastomas (n = 68) and glioblastomas (n = 32) identified chromothripsis as associated with MYC/MYCN gains and with DNA repair deficiencies, pointing towards therapeutic opportunities to target DNA repair defects in tumors with complex genomic rearrangements.

Pediatric high-grade gliomas of the subclass MYCN (HGG-MYCN) are highly aggressive tumors frequently carrying MYCN amplifications, TP53 mutations, or both alterations. Due to their rarity, such tumors have only recently been identified as a distinct entity, and biological as well as clinical characteristics have not been addressed specifically. To gain insights into tumorigenesis and molecular profiles of these tumors, and to ultimately suggest alternative treatment options, we generated a genetically engineered mouse model by breeding hGFAP-cre::Trp53Fl/Fl::lsl-MYCN mice. All mice developed aggressive forebrain tumors early in their lifetime that mimic human HGG-MYCN regarding histology, DNA methylation, and gene expression. Single-cell RNA sequencing revealed a high intratumoral heterogeneity with neuronal and oligodendroglial lineage signatures. High-throughput drug screening using both mouse and human tumor cells finally indicated high efficacy of Doxorubicin, Irinotecan, and Etoposide as possible therapy options that children with HGG-MYCN might benefit from.