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.

Medulloblastoma, which is the most common malignant paediatric brain tumour, has a 70% survival rate, but standard treatments often lead to devastating life-long side effects and recurrence is fatal. One of the emerging strategies in the search for treatments is to determine the roles of tumour microenvironment cells in the growth and maintenance of tumours. The most attractive target is tumour-associated macrophages (TAMs), which are abundantly present in the Sonic Hedgehog (SHH) subgroup of medulloblastoma. Here, we report an unexpected beneficial role of TAMs in SHH medulloblastoma. In human patients, decreased macrophage number is correlated with significantly poorer outcome. We confirm macrophage anti-tumoural behaviour in both ex vivo and in vivo murine models of SHH medulloblastoma. Taken together, our findings suggest that macrophages play a positive role by impairing tumour growth in medulloblastoma, in contrast to the pro-tumoural role played by TAMs in glioblastoma, another common brain tumour.

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.

Kidney tumours are among the most common solid tumours in children, comprising distinct subtypes differing in many aspects, including cell-of-origin, genetics, and pathology. Pre-clinical cell models capturing the disease heterogeneity are currently lacking. Here, we describe the first paediatric cancer organoid biobank. It contains tumour and matching normal kidney organoids from over 50 children with different subtypes of kidney cancer, including Wilms tumours, malignant rhabdoid tumours, renal cell carcinomas, and congenital mesoblastic nephromas. Paediatric kidney tumour organoids retain key properties of native tumours, useful for revealing patient-specific drug sensitivities. Using single cell RNA-sequencing and high resolution 3D imaging, we further demonstrate that organoid cultures derived from Wilms tumours consist of multiple different cell types, including epithelial, stromal and blastemal-like cells. Our organoid biobank captures the heterogeneity of paediatric kidney tumours, providing a representative collection of well-characterised models for basic cancer research, drug-screening and personalised medicine.

Pediatric high-grade glioma (pHGG) is a major contributor to cancer-related death in children. In vitro and in vivo disease models reflecting the intimate connection between developmental context and pathogenesis of pHGG are essential to advance understanding and identify therapeutic vulnerabilities. Here we report establishment of 21 patient-derived pHGG orthotopic xenograft (PDOX) models and eight matched cell lines from diverse groups of pHGG. These models recapitulate histopathology, DNA methylation signatures, mutations and gene expression patterns of the patient tumors from which they were derived, and include rare subgroups not well-represented by existing models. We deploy 16 new and existing cell lines for high-throughput screening (HTS). In vitro HTS results predict variable in vivo response to PI3K/mTOR and MEK pathway inhibitors. These unique new models and an online interactive data portal for exploration of associated detailed molecular characterization and HTS chemical sensitivity data provide a rich resource for pediatric brain tumor research.

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.

Richter syndrome (RS) is the transformation of chronic lymphocytic leukemia (CLL) into aggressive lymphoma, most commonly diffuse large B-cell lymphoma (DLBCL). We characterize 58 primary human RS samples by genome-wide DNA methylation and whole-transcriptome profiling. Our comprehensive approach determines RS DNA methylation profile and unravels a CLL epigenetic imprint, allowing CLL-RS clonal relationship assessment without the need of the initial CLL tumor DNA. DNA methylation- and transcriptomic-based classifiers were developed, and testing on landmark DLBCL datasets identifies a poor-prognosis, activated B-cell-like DLBCL subset in 111/1772 samples. The classification robustly identifies phenotypes very similar to RS with a specific genomic profile, accounting for 4.3-8.3% of de novo DLBCLs. In this work, RS multi-omics characterization determines oncogenic mechanisms, establishes a surrogate marker for CLL-RS clonal relationship, and provides a clinically relevant classifier for a subset of primary "RS-type DLBCL" with unfavorable prognosis.

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.

Targeting oncogenic pathways holds promise for brain tumor treatment, but inhibition of Sonic Hedgehog (SHH) signaling has failed in SHH-driven medulloblastoma. Cellular diversity within tumors and reduced lineage commitment can undermine targeted therapy by increasing the probability of treatment-resistant populations. Using single-cell RNA-seq and lineage tracing, we analyzed cellular diversity in medulloblastomas in transgenic, medulloblastoma-prone mice, and responses to the SHH-pathway inhibitor vismodegib. In untreated tumors, we find expected stromal cells and tumor-derived cells showing either a spectrum of neural progenitor-differentiation states or glial and stem cell markers. Vismodegib reduces the proliferative population and increases differentiation. However, specific cell types in vismodegib-treated tumors remain proliferative, showing either persistent SHH-pathway activation or stem cell characteristics. Our data show that even in tumors with a single pathway-activating mutation, diverse mechanisms drive tumor growth. This diversity confers early resistance to targeted inhibitor therapy, demonstrating the need to target multiple pathways simultaneously.

Burkitt lymphoma (BL) is the most common B-cell lymphoma in children. Within the International Cancer Genome Consortium (ICGC), we performed whole genome and transcriptome sequencing of 39 sporadic BL. Here, we unravel interaction of structural, mutational, and transcriptional changes, which contribute to MYC oncogene dysregulation together with the pathognomonic IG-MYC translocation. Moreover, by mapping IGH translocation breakpoints, we provide evidence that the precursor of at least a subset of BL is a B-cell poised to express IGHA. We describe the landscape of mutations, structural variants, and mutational processes, and identified a series of driver genes in the pathogenesis of BL, which can be targeted by various mechanisms, including IG-non MYC translocations, germline and somatic mutations, fusion transcripts, and alternative splicing.

Suppressor of Fused (SuFu), a tumour suppressor mutated in medulloblastoma, is a central player of Hh signalling, a pathway crucial for development and deregulated in cancer. Although the control of Gli transcription factors by SuFu is critical in Hh signalling, our understanding of the mechanism regulating this key event remains limited. Here, we show that the Itch/β-arrestin2 complex binds SuFu and induces its Lys63-linked polyubiquitylation without affecting its stability. This process increases the association of SuFu with Gli3, promoting the conversion of Gli3 into a repressor, which keeps Hh signalling off. Activation of Hh signalling antagonises the Itch-dependent polyubiquitylation of SuFu. Notably, different SuFu mutations occurring in medulloblastoma patients are insensitive to Itch activity, thus leading to deregulated Hh signalling and enhancing medulloblastoma cell growth. Our findings uncover mechanisms controlling the tumour suppressive functions of SuFu and reveal that their alterations are implicated in medulloblastoma tumorigenesis.

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.

The DNA repair enzyme O6-methylguanine-DNA methyltransferase (MGMT) is commonly overexpressed in cancers and is implicated in the development of chemoresistance. The use of drugs inhibiting MGMT has been hindered by their haematologic toxicity and inefficiency. As a different strategy to inhibit MGMT we investigated cellular regulators of MGMT expression in multiple cancers. Here we show a significant correlation between Wnt signalling and MGMT expression in cancers with different origin and confirm the findings by bioinformatic analysis and immunofluorescence. We demonstrate Wnt-dependent MGMT gene expression and cellular co-localization between active β-catenin and MGMT. Pharmacological or genetic inhibition of Wnt activity downregulates MGMT expression and restores chemosensitivity of DNA-alkylating drugs in mouse models. These findings have potential therapeutic implications for chemoresistant cancers, especially of brain tumours where the use of temozolomide is frequently used in treatment.

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.

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.

Malignant rhabdoid tumour (MRT) is an often lethal childhood cancer that, like many paediatric tumours, is thought to arise from aberrant fetal development. The embryonic root and differentiation pathways underpinning MRT are not firmly established. Here, we study the origin of MRT by combining phylogenetic analyses and single-cell mRNA studies in patient-derived organoids. Comparison of somatic mutations shared between cancer and surrounding normal tissues places MRT in a lineage with neural crest-derived Schwann cells. Single-cell mRNA readouts of MRT differentiation, which we examine by reverting the genetic driver mutation underpinning MRT, SMARCB1 loss, suggest that cells are blocked en route to differentiating into mesenchyme. Quantitative transcriptional predictions indicate that combined HDAC and mTOR inhibition mimic MRT differentiation, which we confirm experimentally. Our study defines the developmental block of MRT and reveals potential differentiation therapies.