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.

Recent research identified histone H3 K27M mutations to be associated with a dismal prognosis in pediatric diffuse midline glioma (pDMG); however, data on detailed MRI characteristics with respect to H3 K27 mutation status and molecular subgroups (H3.1 and H3.3 K27M mutations) are limited.

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.

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.

Atypical teratoid/rhabdoid tumors (AT/RT) are the most common malignant brain tumors manifesting in infancy. They split into four molecular types. The major three (AT/RT-SHH, AT/RT-TYR, and AT/RT-MYC) all carry mutations in SMARCB1, the fourth quantitatively smaller type is characterized by SMARCA4 mutations (AT/RT-SMARCA4). Molecular characteristics of disease recurrence or metastatic spread, which go along with a particularly dismal outcome, are currently unclear. Here, we investigated tumor tissue from 26 patients affected by AT/RT to identify signatures of recurrences in comparison with matched primary tumor samples. Microscopically, AT/RT recurrences demonstrated a loss of architecture and significantly enhanced mitotic activity as compared to their related primary tumors. Based on DNA methylation profiling, primary tumor and related recurrence were grossly similar, but three out of 26 tumors belonged to a different molecular type or subtype after second surgery compared to related primary lesions. Copy number variations (CNVs) differed in six cases, showing novel gains on chromosome 1q or losses of chromosome 10 in recurrences as the most frequent alterations. To consolidate these observations, our cohort was combined with a data set of unmatched primary and recurrent AT/RT, which demonstrated chromosome 1q gain and 10 loss in 18% (n = 7) and 11% (n = 4) of the recurrences (n = 38) as compared to 7% (n = 3) and 0% (n = 0) in the primary tumors (n = 44), respectively. Similar to the observations made by DNA methylation profiling, RNA sequencing of our cohort revealed AT/RT primary tumors and matched recurrences clustering closely together. However, a number of genes showed significantly altered expression in AT/RT-SHH recurrences. Many of them are known tumor driving growth factors, involved in embryonal development and tumorigenesis, or are cell-cycle-associated. Overall, our work identifies subtle molecular changes that occur in the course of the disease and that may help define novel therapeutic targets for AT/RT recurrences.

Pediatric cancer is the leading cause of disease-related death in children and the need for better therapeutic options remains urgent. Due to the limited number of patients, target and drug development for pediatrics is often supplemented by data from studies focused on adult cancers. Recent evidence shows that pediatric cancers possess different vulnerabilities that should be explored independently from adult cancers.

Posterior fossa type A (PF-EPN-A, PFA) ependymoma are aggressive tumors that mainly affect children and have a poor prognosis. Histopathology shows significant intratumoral heterogeneity, ranging from loose tissue to often sharply demarcated, extremely cell-dense tumor areas. To determine molecular differences in morphologically different areas and to understand their clinical significance, we analyzed 113 PF-EPN-A samples, including 40 corresponding relapse samples. Cell-dense areas ranged from 0 to 100% of the tumor area and displayed a higher proportion of proliferating tumor cells (p < 0.01). Clinically, cell density was associated with poor progression-free and overall survival (pPFS = 0.0026, pOS < 0.01). Molecularly, tumor areas with low and high cell density showed diverging DNA methylation profiles regarding their similarity to distinct previously discovered PF-EPN-A subtypes in 9/21 cases. Prognostically relevant chromosomal changes at 1q and 6q showed spatial heterogeneity within single tumors and were significantly enriched in cell-dense tumor areas as shown by single-cell RNA (scRNA)-sequencing as well as copy number profiling and fluorescence in situ hybridization (FISH) analyses of different tumor areas. Finally, spatial transcriptomics revealed cell-dense areas of different tumors to be more similar than various different areas of the same tumor. High-density areas distinctly overexpressed genes encoding histone proteins, WNT5A, TGFB1, or IGF2. Relapsing tumors displayed a higher proportion of cell-dense areas (p = 0.036), a change in PF-EPN-A methylation subtypes (13/32 patients), and novel chromosome 1q gains and 6q losses (12/32 cases) compared to corresponding primary tumors. Our data suggest that PF-EPN-A ependymomas habor a previously unrecognized intratumoral heterogeneity with clinical implications, which has to be accounted for when selecting diagnostic material, inter alia, by histological evaluation of the proportion of cell-dense areas.

In 2012, an international consensus paper reported that medulloblastoma comprises four molecular subgroups (WNT, SHH, Group 3, and Group 4), each associated with distinct genomic features and clinical behavior. Independently, multiple recent reports have defined further intra-subgroup heterogeneity in the form of biologically and clinically relevant subtypes. However, owing to differences in patient cohorts and analytical methods, estimates of subtype number and definition have been inconsistent, especially within Group 3 and Group 4. Herein, we aimed to reconcile the definition of Group 3/Group 4 MB subtypes through the analysis of a series of 1501 medulloblastomas with DNA-methylation profiling data, including 852 with matched transcriptome data. Using multiple complementary bioinformatic approaches, we compared the concordance of subtype calls between published cohorts and analytical methods, including assessments of class-definition confidence and reproducibility. While the lowest complexity solutions continued to support the original consensus subgroups of Group 3 and Group 4, our analysis most strongly supported a definition comprising eight robust Group 3/Group 4 subtypes (types I-VIII). Subtype II was consistently identified across all component studies, while all others were supported by multiple class-definition methods. Regardless of analytical technique, increasing cohort size did not further increase the number of identified Group 3/Group 4 subtypes. Summarizing the molecular and clinico-pathological features of these eight subtypes indicated enrichment of specific driver gene alterations and cytogenetic events amongst subtypes, and identified highly disparate survival outcomes, further supporting their biological and clinical relevance. Collectively, this study provides continued support for consensus Groups 3 and 4 while enabling robust derivation of, and categorical accounting for, the extensive intertumoral heterogeneity within Groups 3 and 4, revealed by recent high-resolution subclassification approaches. Furthermore, these findings provide a basis for application of emerging methods (e.g., proteomics/single-cell approaches) which may additionally inform medulloblastoma subclassification. Outputs from this study will help shape definition of the next generation of medulloblastoma clinical protocols and facilitate the application of enhanced molecularly guided risk stratification to improve outcomes and quality of life for patients and their families.

Functional heterogeneity within tumors presents a significant therapeutic challenge. Here we show that quiescent, therapy-resistant Sox2(+) cells propagate sonic hedgehog subgroup medulloblastoma by a mechanism that mirrors a neurogenic program. Rare Sox2(+) cells produce rapidly cycling doublecortin(+) progenitors that, together with their postmitotic progeny expressing NeuN, comprise tumor bulk. Sox2(+) cells are enriched following anti-mitotic chemotherapy and Smoothened inhibition, creating a reservoir for tumor regrowth. Lineage traces from Sox2(+) cells increase following treatment, suggesting that this population is responsible for relapse. Targeting Sox2(+) cells with the antineoplastic mithramycin abrogated tumor growth. Addressing functional heterogeneity and eliminating Sox2(+) cells presents a promising therapeutic paradigm for treatment of sonic hedgehog subgroup medulloblastoma.

One key advantage of the CRISPR/Cas9 system in comparison with other gene editing approaches lies in its potential for multiplexing. Here, we describe an elaborate procedure that allows the assembly of multiple gRNA expression cassettes into a vector of choice within a single step, termed ASAP(Adaptable System for Assembly of multiplexed Plasmids)-cloning. We demonstrate the utility of ASAP-cloning for multiple CRISPR-mediated applications, including efficient multiplex gene editing, robust transcription activation and convenient analysis of Cas9 activity in the presence of multiple gRNAs.

Pilocytic astrocytoma is the most common childhood brain tumor, characterized by constitutive MAPK activation. MAPK signaling induces oncogene-induced senescence (OIS), which may cause unpredictable growth behavior of pilocytic astrocytomas. The senescence-associated secretory phenotype (SASP) has been shown to regulate OIS, but its role in pilocytic astrocytoma remains unknown.Experimental Design: The patient-derived pilocytic astrocytoma cell culture model, DKFZ-BT66, was used to demonstrate presence of the SASP and analyze its impact on OIS in pilocytic astrocytoma. The model allows for doxycycline-inducible switching between proliferation and OIS. Both states were studied using gene expression profiling (GEP), Western blot, ELISA, and cell viability testing. Primary pilocytic astrocytoma tumors were analyzed by GEP and multiplex assay.

The TCF4 gene encodes for the basic helix-loop-helix transcription factor 4 (TCF4), which plays an important role in the development of the central nervous system (CNS). Haploinsufficiency of TCF4 was found to cause Pitt-Hopkins syndrome (PTHS), a severe neurodevelopmental disorder. Recently, the screening of a large cohort of medulloblastoma (MB), a highly aggressive embryonal brain tumor, revealed almost 20% of adult patients with MB of the Sonic hedgehog (SHH) subtype carrying somatic TCF4 mutations. Interestingly, many of these mutations have previously been detected as germline mutations in patients with PTHS. We show here that overexpression of wild-type TCF4 in vitro significantly suppresses cell proliferation in MB cells, whereas mutant TCF4 proteins do not to the same extent. Furthermore, RNA sequencing revealed significant upregulation of multiple well-known tumor suppressors upon expression of wild-type TCF4. In vivo, a prenatal knockout of Tcf4 in mice caused a significant increase in apoptosis accompanied by a decreased proliferation and failed migration of cerebellar granule neuron precursor cells (CGNP), which are thought to be the cells of origin for SHH MB. In contrast, postnatal in vitro and in vivo knockouts of Tcf4 with and without an additional constitutive activation of the SHH pathway led to significantly increased proliferation of CGNP or MB cells. Finally, publicly available data from human MB show that relatively low expression levels of TCF4 significantly correlate with a worse clinical outcome. These results not only point to time-specific roles of Tcf4 during cerebellar development but also suggest a functional linkage between TCF4 mutations and the formation of SHH MB, proposing that TCF4 acts as a tumor suppressor during postnatal stages of cerebellar development.