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.

Posterior fossa type A (PFA) ependymomas exhibit very low H3K27 methylation and express high levels of EZHIP (Enhancer of Zeste Homologs Inhibitory Protein, also termed CXORF67). Here we find that a conserved sequence in EZHIP is necessary and sufficient to inhibit PRC2 catalytic activity in vitro and in vivo. EZHIP directly contacts the active site of the EZH2 subunit in a mechanism similar to the H3 K27M oncohistone. Furthermore, expression of H3 K27M or EZHIP in cells promotes similar chromatin profiles: loss of broad H3K27me3 domains, but retention of H3K27me3 at CpG islands. We find that H3K27me3-mediated allosteric activation of PRC2 substantially increases the inhibition potential of EZHIP and H3 K27M, providing a mechanism to explain the observed loss of H3K27me3 spreading in tumors. Our data indicate that PFA ependymoma and DIPG are driven in part by the action of peptidyl PRC2 inhibitors, the K27M oncohistone and the EZHIP 'oncohistone-mimic', that dysregulate gene silencing to promote tumorigenesis.

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.

Diffuse Intrinsic Pontine Gliomas (DIPGs) are deadly paediatric brain tumours where needle biopsies help guide diagnosis and targeted therapies. To address spatial heterogeneity, here we analyse 134 specimens from various neuroanatomical structures of whole autopsy brains from nine DIPG patients. Evolutionary reconstruction indicates histone 3 (H3) K27M--including H3.2K27M--mutations potentially arise first and are invariably associated with specific, high-fidelity obligate partners throughout the tumour and its spread, from diagnosis to end-stage disease, suggesting mutual need for tumorigenesis. These H3K27M ubiquitously-associated mutations involve alterations in TP53 cell-cycle (TP53/PPM1D) or specific growth factor pathways (ACVR1/PIK3R1). Later oncogenic alterations arise in sub-clones and often affect the PI3K pathway. Our findings are consistent with early tumour spread outside the brainstem including the cerebrum. The spatial and temporal homogeneity of main driver mutations in DIPG implies they will be captured by limited biopsies and emphasizes the need to develop therapies specifically targeting obligate oncohistone partnerships.

Our ability to manage acute myeloid leukemia (AML) is limited by our incomplete understanding of the epigenetic disruption central to leukemogenesis, including improper histone methylation. Here we examine 16 histone H3 genes in 434 primary AML samples and identify Q69H, A26P, R2Q, R8H and K27M/I mutations (1.6%), with higher incidence in secondary AML (9%). These mutations occur in pre-leukemic hematopoietic stem cells (HSCs) and exist in the major leukemic clones in patients. They increase the frequency of functional HSCs, alter differentiation, and amplify leukemic aggressiveness. These effects are dependent on the specific mutation. H3K27 mutation increases the expression of genes involved in erythrocyte and myeloid differentiation with altered H3K27 tri-methylation and K27 acetylation. The functional impact of histone mutations is independent of RUNX1 mutation, although they at times co-occur. This study establishes that H3 mutations are drivers of human pre-cancerous stem cell expansion and important early events in leukemogenesis.

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.

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.

Giant cell lesions of the jaw (GCLJ) are debilitating tumors of unknown origin with limited available therapies. Here, we analyze 58 sporadic samples using next generation or targeted sequencing and report somatic, heterozygous, gain-of-function mutations in KRAS, FGFR1, and p.M713V/I-TRPV4 in 72% (42/58) of GCLJ. TRPV4 p.M713V/I mutations are exclusive to central GCLJ and occur at a critical position adjacent to the cation permeable pore of the channel. Expression of TRPV4 mutants in HEK293 cells leads to increased cell death, as well as increased constitutive and stimulated channel activity, both of which can be prevented using TRPV4 antagonists. Furthermore, these mutations induce sustained activation of ERK1/2, indicating that their effects converge with that of KRAS and FGFR1 mutations on the activation of the MAPK pathway in GCLJ. Our data extend the spectrum of TRPV4 channelopathies and provide rationale for the use of TRPV4 and RAS/MAPK antagonists at the bedside in GCLJ.

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.

BET-bromodomain inhibition (BETi) has shown pre-clinical promise for MYC-amplified medulloblastoma. However, the mechanisms for its action, and ultimately for resistance, have not been fully defined. Here, using a combination of expression profiling, genome-scale CRISPR/Cas9-mediated loss of function and ORF/cDNA driven rescue screens, and cell-based models of spontaneous resistance, we identify bHLH/homeobox transcription factors and cell-cycle regulators as key genes mediating BETi's response and resistance. Cells that acquire drug tolerance exhibit a more neuronally differentiated cell-state and expression of lineage-specific bHLH/homeobox transcription factors. However, they do not terminally differentiate, maintain expression of CCND2, and continue to cycle through S-phase. Moreover, CDK4/CDK6 inhibition delays acquisition of resistance. Therefore, our data provide insights about the mechanisms underlying BETi effects and the appearance of resistance and support the therapeutic use of combined cell-cycle inhibitors with BETi in MYC-amplified medulloblastoma.

Diffuse intrinsic pontine glioma (DIPG) is an incurable malignant childhood brain tumor, with no active systemic therapies and a 5-year survival of less than 1%. Polyamines are small organic polycations that are essential for DNA replication, translation and cell proliferation. Ornithine decarboxylase 1 (ODC1), the rate-limiting enzyme in polyamine synthesis, is irreversibly inhibited by difluoromethylornithine (DFMO). Herein we show that polyamine synthesis is upregulated in DIPG, leading to sensitivity to DFMO. DIPG cells compensate for ODC1 inhibition by upregulation of the polyamine transporter SLC3A2. Treatment with the polyamine transporter inhibitor AMXT 1501 reduces uptake of polyamines in DIPG cells, and co-administration of AMXT 1501 and DFMO leads to potent in vitro activity, and significant extension of survival in three aggressive DIPG orthotopic animal models. Collectively, these results demonstrate the potential of dual targeting of polyamine synthesis and uptake as a therapeutic strategy for incurable DIPG.

Sonic hedgehog medulloblastoma encompasses a clinically and molecularly diverse group of cancers of the developing central nervous system. Here, we use unbiased sequencing of the transcriptome across a large cohort of 250 tumors to reveal differences among molecular subtypes of the disease, and demonstrate the previously unappreciated importance of non-coding RNA transcripts. We identify alterations within the cAMP dependent pathway (GNAS, PRKAR1A) which converge on GLI2 activity and show that 18% of tumors have a genetic event that directly targets the abundance and/or stability of MYCN. Furthermore, we discover an extensive network of fusions in focally amplified regions encompassing GLI2, and several loss-of-function fusions in tumor suppressor genes PTCH1, SUFU and NCOR1. Molecular convergence on a subset of genes by nucleotide variants, copy number aberrations, and gene fusions highlight the key roles of specific pathways in the pathogenesis of Sonic hedgehog medulloblastoma and open up opportunities for therapeutic intervention.

Ependymoma is a tumor of the brain or spinal cord. The two most common and aggressive molecular groups of ependymoma are the supratentorial ZFTA-fusion associated and the posterior fossa ependymoma group A. In both groups, tumors occur mainly in young children and frequently recur after treatment. Although molecular mechanisms underlying these diseases have recently been uncovered, they remain difficult to target and innovative therapeutic approaches are urgently needed. Here, we use genome-wide chromosome conformation capture (Hi-C), complemented with CTCF and H3K27ac ChIP-seq, as well as gene expression and DNA methylation analysis in primary and relapsed ependymoma tumors, to identify chromosomal conformations and regulatory mechanisms associated with aberrant gene expression. In particular, we observe the formation of new topologically associating domains ('neo-TADs') caused by structural variants, group-specific 3D chromatin loops, and the replacement of CTCF insulators by DNA hyper-methylation. Through inhibition experiments, we validate that genes implicated by these 3D genome conformations are essential for the survival of patient-derived ependymoma models in a group-specific manner. Thus, this study extends our ability to reveal tumor-dependency genes by 3D genome conformations even in tumors that lack targetable genetic alterations.