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.

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.

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.

T cell exhaustion limits anti-tumor immunity and responses to immunotherapy. Here, we explored the microenvironmental signals regulating T cell exhaustion using a model of chronic lymphocytic leukemia (CLL). Single-cell analyses identified a subset of PD-1hi, functionally impaired CD8+ T cells that accumulated in secondary lymphoid organs during disease progression and a functionally competent PD-1int subset. Frequencies of PD-1int TCF-1+ CD8+ T cells decreased upon Il10rb or Stat3 deletion, leading to accumulation of PD-1hi cells and accelerated tumor progression. Mechanistically, inhibition of IL-10R signaling altered chromatin accessibility and disrupted cooperativity between the transcription factors NFAT and AP-1, promoting a distinct NFAT-associated program. Low IL10 expression or loss of IL-10R-STAT3 signaling correlated with increased frequencies of exhausted CD8+ T cells and poor survival in CLL and in breast cancer patients. Thus, balance between PD-1hi, exhausted CD8+ T cells and functional PD-1int TCF-1+ CD8+ T cells is regulated by cell-intrinsic IL-10R signaling, with implications for immunotherapy.

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.

The maternal mode of mitochondrial DNA (mtDNA) inheritance is central to human genetics. Recently, evidence for bi-parental inheritance of mtDNA was claimed for individuals of three pedigrees that suffered mitochondrial disorders. We sequenced mtDNA using both direct Sanger and Massively Parallel Sequencing in several tissues of eleven maternally related and other affiliated healthy individuals of a family pedigree and observed mixed mitotypes in eight individuals. Cells without nuclear DNA, i.e. thrombocytes and hair shafts, only showed the mitotype of haplogroup (hg) V. Skin biopsies were prepared to generate ρ° cells void of mtDNA, sequencing of which resulted in a hg U4c1 mitotype. The position of the Mega-NUMT sequence was determined by fluorescence in situ hybridization and two different quantitative PCR assays were used to determine the number of contributing mtDNA copies. Thus, evidence for the presence of repetitive, full mitogenome Mega-NUMTs matching haplogroup U4c1 in various tissues of eight maternally related individuals was provided. Multi-copy Mega-NUMTs mimic mixtures of mtDNA that cannot be experimentally avoided and thus may appear in diverse fields of mtDNA research and diagnostics. We demonstrate that hair shaft mtDNA sequencing provides a simple but reliable approach to exclude NUMTs as source of misleading results.

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.

Clonal evolution is involved in the progression of chronic lymphocytic leukemia (CLL). In order to link evolutionary patterns to different disease courses, we performed a long-term longitudinal mutation profiling study of CLL patients. Tracking somatic mutations and their changes in allele frequency over time and assessing the underlying cancer cell fraction revealed highly distinct evolutionary patterns. Surprisingly, in long-term stable disease and in relapse after long-lasting clinical response to treatment, clonal shifts are minor. In contrast, in refractory disease major clonal shifts occur although there is little impact on leukemia cell counts. As this striking pattern in refractory cases is not linked to a strong contribution of known CLL driver genes, the evolution is mostly driven by treatment-induced selection of sub-clones, underlining the need for novel, non-genotoxic treatment regimens.

The identity of the cell of origin is a key determinant of cancer subtype, progression, and prognosis. Group 3 medulloblastoma (MB) is a malignant childhood brain cancer with poor prognosis and few candidates as putative cell of origin. We overexpressed the group 3 MB genetic drivers MYC and Gfi1 in different candidate cells of origin in the postnatal mouse cerebellum. We found that S100b+ cells are competent to initiate group 3 MB, and we observed that S100b+ cells have higher levels of Notch1 pathway activity compared to Math1+ cells. We found that additional activation of Notch1 in Math1+ and Sox2+ cells was sufficient to induce group 3 MB upon MYC/Gfi1 expression. Together, our data suggest that the Notch1 pathway plays a critical role in group 3 MB initiation.

In chronic lymphocytic leukemia (CLL), a diverse set of genetic mutations is embedded in a deregulated epigenetic landscape that drives cancerogenesis. To elucidate the role of aberrant chromatin features, we mapped DNA methylation, seven histone modifications, nucleosome positions, chromatin accessibility, binding of EBF1 and CTCF, as well as the transcriptome of B cells from CLL patients and healthy donors. A globally increased histone deacetylase activity was detected and half of the genome comprised transcriptionally downregulated partially DNA methylated domains demarcated by CTCF CLL samples displayed a H3K4me3 redistribution and nucleosome gain at promoters as well as changes of enhancer activity and enhancer linkage to target genes. A DNA binding motif analysis identified transcription factors that gained or lost binding in CLL at sites with aberrant chromatin features. These findings were integrated into a gene regulatory enhancer containing network enriched for B-cell receptor signaling pathway components. Our study predicts novel molecular links to targets of CLL therapies and provides a valuable resource for further studies on the epigenetic contribution to the disease.

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.