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.

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.

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.

Ependymoma with RELA fusion has been defined as a novel entity of the revised World Health Organization 2016 classification of tumors of the central nervous system (CNS), characterized by fusion transcripts of the RELA gene and consequent pathological activation of the NFkB pathway. These tumors represent the majority of supratentorial ependymomas in children. The validation of diagnostic tools to identify this clinically relevant ependymoma entity is essential. Here, we have used interphase fluorescent in situ hybridization (FISH) for C11orf95 and RELA, immunohistochemistry (IHC) for p65-RelA and the recently developed DNA methylation-based classification besides conventional histopathology, and compared the precision of the methods in 40 supratentorial pediatric brain tumors diagnosed as ependymomas in the past years. Reverse transcription PCR (RT-PCR) and RNA sequencing were performed to explore discordant cases. Furthermore, we integrated imaging and clinical features as additional layers of information. The concordance between nuclear RelA expression by IHC and RELA FISH was 100%. Concordance between IHC and DNA methylation profiling, and between FISH and DNA methylation profiling was also high (96.4% and 95.2%, respectively). Thirty-four out of 40 (85%) cases were confirmed by integrated diagnoses as ependymal tumors, including 22 RELA-fused ependymomas (71% of ependymal tumors), two YAP1-fused ependymomas (6%), six non-RELA/non-YAP1 ependymomas (18%) and four ependymal/subependymal mixed tumors (12%). Ependymal/subependymal mixed tumors had an excellent clinical outcome despite the presence of histopathological signs of malignancy, suggesting that these tumors should not be diagnosed as classic ependymomas. DNA methylation profiling helped in the differential diagnosis of RELA-fused ependymomas. IHC and FISH, which are available in the majority of pathology laboratories, are valuable tools to identify RELA-fused ependymomas.

The current consensus recognizes four main medulloblastoma subgroups (wingless, Sonic hedgehog, group 3 and group 4). While medulloblastoma subgroups have been characterized extensively at the (epi-)genomic and transcriptomic levels, the proteome and phosphoproteome landscape remain to be comprehensively elucidated. Using quantitative (phospho)-proteomics in primary human medulloblastomas, we unravel distinct posttranscriptional regulation leading to highly divergent oncogenic signaling and kinase activity profiles in groups 3 and 4 medulloblastomas. Specifically, proteomic and phosphoproteomic analyses identify aberrant ERBB4-SRC signaling in group 4. Hence, enforced expression of an activated SRC combined with p53 inactivation induces murine tumors that resemble group 4 medulloblastoma. Therefore, our integrative proteogenomics approach unveils an oncogenic pathway and potential therapeutic vulnerability in the most common medulloblastoma subgroup.

Human neural stem cell cultures provide progenitor cells that are potential cells of origin for brain cancers. However, the extent to which genetic predisposition to tumor formation can be faithfully captured in stem cell lines is uncertain. Here, we evaluated neuroepithelial stem (NES) cells, representative of cerebellar progenitors. We transduced NES cells with MYCN, observing medulloblastoma upon orthotopic implantation in mice. Significantly, transcriptomes and patterns of DNA methylation from xenograft tumors were globally more representative of human medulloblastoma compared to a MYCN-driven genetically engineered mouse model. Orthotopic transplantation of NES cells generated from Gorlin syndrome patients, who are predisposed to medulloblastoma due to germline-mutated PTCH1, also generated medulloblastoma. We engineered candidate cooperating mutations in Gorlin NES cells, with mutation of DDX3X or loss of GSE1 both accelerating tumorigenesis. These findings demonstrate that human NES cells provide a potent experimental resource for dissecting genetic causation in medulloblastoma.

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.

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 invasion is the primary cause of treatment failure of glioblastoma (GBM). Previous studies on GBM invasion have long been forced to use the resected tumor mass cells. Here, a strategy to reliably isolate matching pairs of invasive (GBMINV ) and tumor core (GBMTC ) cells from the brains of 6 highly invasive patient-derived orthotopic models is described. Direct comparison of these GBMINV and GBMTC cells reveals a significantly elevated invasion capacity in GBMINV cells, detects 23/768 miRNAs over-expressed in the GBMINV cells (miRNAINV ) and 22/768 in the GBMTC cells (miRNATC ), respectively. Silencing the top 3 miRNAsINV (miR-126, miR-369-5p, miR-487b) successfully blocks invasion of GBMINV cells in vitro and in mouse brains. Integrated analysis with mRNA expression identifies miRNAINV target genes and discovers KCNA1 as the sole common computational target gene of which 3 inhibitors significantly suppress invasion in vitro. Furthermore, in vivo treatment with 4-aminopyridine (4-AP) effectively eliminates GBM invasion and significantly prolongs animal survival times (P = 0.035). The results highlight the power of spatial dissection of functionally accurate GBMINV and GBMTC cells in identifying novel drivers of GBM invasion and provide strong rationale to support the use of biologically accurate starting materials in understanding cancer invasion and metastasis.

Large-scale molecular profiling studies in recent years have shown that central nervous system (CNS) tumors display a much greater heterogeneity in terms of molecularly distinct entities, cellular origins and genetic drivers than anticipated from histological assessment. DNA methylation profiling has emerged as a useful tool for robust tumor classification, providing new insights into these heterogeneous molecular classes. This is particularly true for rare CNS tumors with a broad morphological spectrum, which are not possible to assign as separate entities based on histological similarity alone. Here, we describe a molecularly distinct subset of predominantly pediatric CNS neoplasms (n = 60) that harbor PATZ1 fusions. The original histological diagnoses of these tumors covered a wide spectrum of tumor types and malignancy grades. While the single most common diagnosis was glioblastoma (GBM), clinical data of the PATZ1-fused tumors showed a better prognosis than typical GBM, despite frequent relapses. RNA sequencing revealed recurrent MN1:PATZ1 or EWSR1:PATZ1 fusions related to (often extensive) copy number variations on chromosome 22, where PATZ1 and the two fusion partners are located. These fusions have individually been reported in a number of glial/glioneuronal tumors, as well as extracranial sarcomas. We show here that they are more common than previously acknowledged, and together define a biologically distinct CNS tumor type with high expression of neural development markers such as PAX2, GATA2 and IGF2. Drug screening performed on the MN1:PATZ1 fusion-bearing KS-1 brain tumor cell line revealed preliminary candidates for further study. In summary, PATZ1 fusions define a molecular class of histologically polyphenotypic neuroepithelial tumors, which show an intermediate prognosis under current treatment regimens.

Ependymomas encompass a heterogeneous group of central nervous system (CNS) neoplasms that occur along the entire neuroaxis. In recent years, extensive (epi-)genomic profiling efforts have identified several molecular groups of ependymoma that are characterized by distinct molecular alterations and/or patterns. Based on unsupervised visualization of a large cohort of genome-wide DNA methylation data, we identified a highly distinct group of pediatric-type tumors (n = 40) forming a cluster separate from all established CNS tumor types, of which a high proportion were histopathologically diagnosed as ependymoma. RNA sequencing revealed recurrent fusions involving the pleomorphic adenoma gene-like 1 (PLAGL1) gene in 19 of 20 of the samples analyzed, with the most common fusion being EWSR1:PLAGL1 (n = 13). Five tumors showed a PLAGL1:FOXO1 fusion and one a PLAGL1:EP300 fusion. High transcript levels of PLAGL1 were noted in these tumors, with concurrent overexpression of the imprinted genes H19 and IGF2, which are regulated by PLAGL1. Histopathological review of cases with sufficient material (n = 16) demonstrated a broad morphological spectrum of tumors with predominant ependymoma-like features. Immunohistochemically, tumors were GFAP positive and OLIG2- and SOX10 negative. In 3/16 of the cases, a dot-like positivity for EMA was detected. All tumors in our series were located in the supratentorial compartment. Median age of the patients at the time of diagnosis was 6.2 years. Median progression-free survival was 35 months (for 11 patients with data available). In summary, our findings suggest the existence of a novel group of supratentorial neuroepithelial tumors that are characterized by recurrent PLAGL1 fusions and enriched for pediatric patients.

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.

The large diversity of central nervous system (CNS) tumor types in children and adolescents results in disparate patient outcomes and renders accurate diagnosis challenging. In this study, we prospectively integrated DNA methylation profiling and targeted gene panel sequencing with blinded neuropathological reference diagnostics for a population-based cohort of more than 1,200 newly diagnosed pediatric patients with CNS tumors, to assess their utility in routine neuropathology. We show that the multi-omic integration increased diagnostic accuracy in a substantial proportion of patients through annotation to a refining DNA methylation class (50%), detection of diagnostic or therapeutically relevant genetic alterations (47%) or identification of cancer predisposition syndromes (10%). Discrepant results by neuropathological WHO-based and DNA methylation-based classification (30%) were enriched in histological high-grade gliomas, implicating relevance for current clinical patient management in 5% of all patients. Follow-up (median 2.5 years) suggests improved survival for patients with histological high-grade gliomas displaying lower-grade molecular profiles. These results provide preliminary evidence of the utility of integrating multi-omics in neuropathology for pediatric neuro-oncology.

Poorly differentiated chordoma (PDC) is a recently recognized subtype of chordoma characterized by expression of the embryonic transcription factor, brachyury, and loss of INI1. PDC primarily affects children and is associated with a poor prognosis and limited treatment options. Here we describe the molecular and immune tumour microenvironment profiles of two paediatric PDCs produced using whole-genome, transcriptome and whole-genome bisulfite sequencing (WGBS) and multiplex immunohistochemistry. Our analyses revealed the presence of tumour-associated immune cells, including CD8+ T cells, and expression of the immune checkpoint protein, PD-L1, in both patient samples. Molecular profiling provided the rationale for immune checkpoint inhibitor (ICI) therapy, which resulted in a clinical and radiographic response. A dominant T cell receptor (TCR) clone specific for a brachyury peptide-MHC complex was identified from bulk RNA sequencing, suggesting that targeting of the brachyury tumour antigen by tumour-associated T cells may underlie this clinical response to ICI. Correlative analysis with rhabdoid tumours, another INI1-deficient paediatric malignancy, suggests that a subset of tumours may share common immune phenotypes, indicating the potential for a therapeutically targetable subgroup of challenging paediatric cancers.

Myxopapillary ependymoma (MPE) is a heterogeneous disease regarding histopathology and outcome. The underlying molecular biology is poorly understood, and markers that reliably predict the patients' clinical course are unknown.

Atypical teratoid/rhabdoid tumor (ATRT) is an aggressive central nervous system tumor characterized by loss of SMARCB1/INI1 protein expression and comprises three distinct molecular groups, ATRT-TYR, ATRT-MYC and ATRT-SHH. ATRT-SHH represents the largest molecular group and is heterogeneous with regard to age, tumor location and epigenetic profile. We, therefore, aimed to investigate if heterogeneity within ATRT-SHH might also have biological and clinical importance. Consensus clustering of DNA methylation profiles and confirmatory t-SNE analysis of 65 ATRT-SHH yielded three robust molecular subgroups, i.e., SHH-1A, SHH-1B and SHH-2. These subgroups differed by median age of onset (SHH-1A: 18 months, SHH-1B: 107 months, SHH-2: 13 months) and tumor location (SHH-1A: 88% supratentorial; SHH-1B: 85% supratentorial; SHH-2: 93% infratentorial, often extending to the pineal region). Subgroups showed comparable SMARCB1 mutational profiles, but pathogenic/likely pathogenic SMARCB1 germline variants were over-represented in SHH-2 (63%) as compared to SHH-1A (20%) and SHH-1B (0%). Protein expression of proneural marker ASCL1 (enriched in SHH-1B) and glial markers OLIG2 and GFAP (absent in SHH-2) as well as global mRNA expression patterns differed, but all subgroups were characterized by overexpression of SHH as well as Notch pathway members. In a Drosophila model, knockdown of Snr1 (the fly homologue of SMARCB1) in hedgehog activated cells not only altered hedgehog signaling, but also caused aberrant Notch signaling and formation of tumor-like structures. Finally, on survival analysis, molecular subgroup and age of onset (but not ASCL1 staining status) were independently associated with overall survival, older patients (> 3 years) harboring SHH-1B experiencing relatively favorable outcome. In conclusion, ATRT-SHH comprises three subgroups characterized by SHH and Notch pathway activation, but divergent molecular and clinical features. Our data suggest that molecular subgrouping of ATRT-SHH has prognostic relevance and might aid to stratify patients within future clinical trials.

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.

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.

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.