Diffuse intrinsic pontine glioma (DIPG) is an incurable brain tumor of childhood characterized by histone mutations at lysine 27, which results in epigenomic dysregulation. There has been a failure to develop effective treatment for this tumor. Using a combined RNAi and chemical screen targeting epigenomic regulators, we identify the polycomb repressive complex 1 (PRC1) component BMI1 as a critical factor for DIPG tumor maintenance in vivo. BMI1 chromatin occupancy is enriched at genes associated with differentiation and tumor suppressors in DIPG cells. Inhibition of BMI1 decreases cell self-renewal and attenuates tumor growth due to induction of senescence. Prolonged BMI1 inhibition induces a senescence-associated secretory phenotype, which promotes tumor recurrence. Clearance of senescent cells using BH3 protein mimetics co-operates with BMI1 inhibition to enhance tumor cell killing in vivo.

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.

High-grade gliomas defined by histone 3 K27M driver mutations exhibit global loss of H3K27 trimethylation and reciprocal gain of H3K27 acetylation, respectively shaping repressive and active chromatin landscapes. We generated tumor-derived isogenic models bearing this mutation and show that it leads to pervasive H3K27ac deposition across the genome. In turn, active enhancers and promoters are not created de novo and instead reflect the epigenomic landscape of the cell of origin. H3K27ac is enriched at repeat elements, resulting in their increased expression, which in turn can be further amplified by DNA demethylation and histone deacetylase inhibitors providing an exquisite therapeutic vulnerability. These agents may therefore modulate anti-tumor immune responses as a therapeutic modality for this untreatable disease.

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.

Histological grading of choroid plexus tumors (CPTs) remains the best prognostic tool to distinguish between aggressive choroid plexus carcinoma (CPC) and the more benign choroid plexus papilloma (CPP) or atypical choroid plexus papilloma (aCPP); however, these distinctions can be challenging. Standard treatment of CPC is very aggressive and often leads to severe damage to the young child's brain. Therefore, it is crucial to distinguish between CPC and less aggressive entities (CPP or aCPP) to avoid unnecessary exposure of the young patient to neurotoxic therapy. To better stratify CPTs, we utilized DNA methylation (DNAm) to identify prognostic epigenetic biomarkers for CPCs.

Spinal ependymal tumors form a histologically and molecularly heterogeneous group of tumors with generally good prognosis. However, their treatment can be challenging if infiltration of the spinal cord or dissemination throughout the central nervous system (CNS) occurs and, in these cases, clinical outcome remains poor. Here, we describe a new and relatively rare subgroup of spinal ependymal tumors identified using DNA methylation profiling that is distinct from other molecular subgroups of ependymoma. Copy number variation plots derived from DNA methylation arrays showed MYCN amplification as a characteristic genetic alteration in all cases of our cohort (n = 13), which was subsequently validated using fluorescence in situ hybridization. The histological diagnosis was anaplastic ependymoma (WHO Grade III) in ten cases and classic ependymoma (WHO Grade II) in three cases. Histological re-evaluation in five primary tumors and seven relapses showed characteristic histological features of ependymoma, namely pseudorosettes, GFAP- and EMA positivity. Electron microscopy revealed cilia, complex intercellular junctions and intermediate filaments in a representative sample. Taking these findings into account, we suggest to designate this molecular subgroup spinal ependymoma with MYCN amplification, SP-EPN-MYCN. SP-EPN-MYCN tumors showed distinct growth patterns with intradural, extramedullary localization mostly within the thoracic and cervical spine, diffuse leptomeningeal spread throughout the whole CNS and infiltrative invasion of the spinal cord. Dissemination was observed in 100% of cases. Despite high-intensity treatment, SP-EPN-MYCN showed significantly worse median progression free survival (PFS) (17 months) and median overall survival (OS) (87 months) than all other previously described molecular spinal ependymoma subgroups. OS and PFS were similar to supratentorial ependymoma with RELA-fusion (ST-EPN-RELA) and posterior fossa ependymoma A (PF-EPN-A), further highlighting the aggressiveness of this distinct new subgroup. We, therefore, propose to establish SP-EPN-MYCN as a new molecular subgroup in ependymoma and advocate for testing newly diagnosed spinal ependymal tumors for MYCN amplification.

The aim of this study is to present a national surveillance report on pediatric central nervous system (CNS) tumors in Canada during the period between 2001 and 2015.

High-intensity therapy effectively treats most TP53 wild-type (TP53-WT) Sonic Hedgehog-subgroup medulloblastomas (SHH-MBs), but often cause long-term deleterious neurotoxicities in children. Recent clinical trials investigating reduction/de-escalation of therapy for TP53-WT SHH-MBs caused poor overall survival. Here, we investigated whether reduced levels of p53-pathway activation by low-intensity therapy potentially contribute to diminished therapeutic efficacy.

Childhood brain tumors have suspected prenatal origins. To identify vulnerable developmental states, we generated a single-cell transcriptome atlas of >65,000 cells from embryonal pons and forebrain, two major tumor locations. We derived signatures for 191 distinct cell populations and defined the regional cellular diversity and differentiation dynamics. Projection of bulk tumor transcriptomes onto this dataset shows that WNT medulloblastomas match the rhombic lip-derived mossy fiber neuronal lineage and embryonal tumors with multilayered rosettes fully recapitulate a neuronal lineage, while group 2a/b atypical teratoid/rhabdoid tumors may originate outside the neuroectoderm. Importantly, single-cell tumor profiles reveal highly defined cell hierarchies that mirror transcriptional programs of the corresponding normal lineages. Our findings identify impaired differentiation of specific neural progenitors as a common mechanism underlying these pediatric cancers and provide a rational framework for future modeling and therapeutic interventions.

Retinoblastoma binding protein 4 (Rbbp4) is a component of transcription regulatory complexes that control cell cycle gene expression. Previous work indicated that Rbbp4 cooperates with the Rb tumor suppressor to block cell cycle entry. Here, we use genetic analysis to examine the interactions of Rbbp4, Rb, and Tp53 in zebrafish neural progenitor cell cycle regulation and survival.

Histone H3.3 lysine-to-methionine substitutions K27M and K36M impair the deposition of opposing chromatin marks, H3K27me3/me2 and H3K36me3/me2. We show that these mutations induce hypotrophic and disorganized eyes in Drosophila eye primordia. Restriction of H3K27me3 spread in H3.3K27M and its redistribution in H3.3K36M result in transcriptional deregulation of PRC2-targeted eye development and of piRNA biogenesis genes, including krimp. Notably, both mutants promote redistribution of H3K36me2 away from repetitive regions into active genes, which associate with retrotransposon de-repression in eye discs. Aberrant expression of krimp represses LINE retrotransposons but does not contribute to the eye phenotype. Depletion of H3K36me2 methyltransferase ash1 in H3.3K27M, and of PRC2 component E(z) in H3.3K36M, restores the expression of eye developmental genes and normal eye growth, showing that redistribution of antagonistic marks contributes to K-to-M pathogenesis. Our results implicate a novel function for H3K36me2 and showcase convergent downstream effects of oncohistones that target opposing epigenetic marks.

Revealing the 3D composition of intact tissue specimens is essential for understanding cell and organ biology in health and disease. State-of-the-art 3D microscopy techniques aim to capture tissue volumes on an ever-increasing scale, while also retaining sufficient resolution for single-cell analysis. Furthermore, spatial profiling through multi-marker imaging is fast developing, providing more context and better distinction between cell types. Following these lines of technological advance, we here present a protocol based on FUnGI (fructose, urea and glycerol clearing solution for imaging) optical clearing of tissue before multispectral large-scale single-cell resolution 3D (mLSR-3D) imaging, which implements 'on-the-fly' linear unmixing of up to eight fluorophores during a single acquisition. Our protocol removes the need for repetitive illumination, thereby allowing larger volumes to be scanned with better image quality in less time, also reducing photo-bleaching and file size. To aid in the design of multiplex antibody panels, we provide a fast and manageable intensity equalization assay with automated analysis to design a combination of markers with balanced intensities suitable for mLSR-3D. We demonstrate effective mLSR-3D imaging of various tissues, including patient-derived organoids and xenografted tumors, and, furthermore, describe an optimized workflow for mLSR-3D imaging of formalin-fixed paraffin-embedded samples. Finally, we provide essential steps for 3D image data processing, including shading correction that does not require pre-acquired shading references and 3D inhomogeneity correction to correct fluorescence artefacts often afflicting 3D datasets. Together, this provides a one-week protocol for eight-fluorescent-marker 3D visualization and exploration of intact tissue of various origins at single-cell resolution.

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.

Study of the origin and development of cerebellar tumours has been hampered by the complexity and heterogeneity of cerebellar cells that change over the course of development. Here we use single-cell transcriptomics to study more than 60,000 cells from the developing mouse cerebellum and show that different molecular subgroups of childhood cerebellar tumours mirror the transcription of cells from distinct, temporally restricted cerebellar lineages. The Sonic Hedgehog medulloblastoma subgroup transcriptionally mirrors the granule cell hierarchy as expected, while group 3 medulloblastoma resembles Nestin+ stem cells, group 4 medulloblastoma resembles unipolar brush cells, and PFA/PFB ependymoma and cerebellar pilocytic astrocytoma resemble the prenatal gliogenic progenitor cells. Furthermore, single-cell transcriptomics of human childhood cerebellar tumours demonstrates that many bulk tumours contain a mixed population of cells with divergent differentiation. Our data highlight cerebellar tumours as a disorder of early brain development and provide a proximate explanation for the peak incidence of cerebellar tumours in early childhood.

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

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.

Enzymes that catalyse CpG methylation in DNA, including the DNA methyltransferases 1 (DNMT1), 3A (DNMT3A) and 3B (DNMT3B), are indispensable for mammalian tissue development and homeostasis1-4. They are also implicated in human developmental disorders and cancers5-8, supporting the critical role of DNA methylation in the specification and maintenance of cell fate. Previous studies have suggested that post-translational modifications of histones are involved in specifying patterns of DNA methyltransferase localization and DNA methylation at promoters and actively transcribed gene bodies9-11. However, the mechanisms that control the establishment and maintenance of intergenic DNA methylation remain poorly understood. Tatton-Brown-Rahman syndrome (TBRS) is a childhood overgrowth disorder that is defined by germline mutations in DNMT3A. TBRS shares clinical features with Sotos syndrome (which is caused by haploinsufficiency of NSD1, a histone methyltransferase that catalyses the dimethylation of histone H3 at K36 (H3K36me2)8,12,13), which suggests that there is a mechanistic link between these two diseases. Here we report that NSD1-mediated H3K36me2 is required for the recruitment of DNMT3A and maintenance of DNA methylation at intergenic regions. Genome-wide analysis shows that the binding and activity of DNMT3A colocalize with H3K36me2 at non-coding regions of euchromatin. Genetic ablation of Nsd1 and its paralogue Nsd2 in mouse cells results in a redistribution of DNMT3A to H3K36me3-modified gene bodies and a reduction in the methylation of intergenic DNA. Blood samples from patients with Sotos syndrome and NSD1-mutant tumours also exhibit hypomethylation of intergenic DNA. The PWWP domain of DNMT3A shows dual recognition of H3K36me2 and H3K36me3 in vitro, with a higher binding affinity towards H3K36me2 that is abrogated by TBRS-derived missense mutations. Together, our study reveals a trans-chromatin regulatory pathway that connects aberrant intergenic CpG methylation to human neoplastic and developmental overgrowth.

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.