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.

We reconstructed the natural history and temporal evolution of the most common childhood brain malignancy, medulloblastoma, by single-cell whole-genome sequencing (sc-WGS) of tumours representing its major molecular sub-classes and clinical risk groups. Favourable-risk disease sub-types assessed (MBWNT and infant desmoplastic/nodular MBSHH) typically comprised a single clone with no evidence of further evolution. In contrast, highest risk sub-classes (MYC-amplified MBGroup3 and TP53-mutated MBSHH) were most clonally diverse and displayed gradual evolutionary trajectories. Clinically adopted biomarkers (e.g. chromosome 6/17 aberrations; CTNNB1/TP53 mutations) were typically early-clonal/initiating events, exploitable as targets for early-disease detection; in analyses of spatially distinct tumour regions, a single biopsy was sufficient to assess their status. Importantly, sc-WGS revealed novel events which arise later and/or sub-clonally and more commonly display spatial diversity; their clinical significance and role in disease evolution post-diagnosis now require establishment. These findings reveal diverse modes of tumour initiation and evolution in the major medulloblastoma sub-classes, with pathogenic relevance and clinical potential.

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.

Here, we present a protocol for deriving a continuum score for group 3 and 4 medulloblastoma tumor samples analyzed via RNA-sequencing or DNA methylation microarray. We describe steps for utilizing NMF-defined group 3/group 4 metagenes to calculate a continuum score between 0 and 1 that can be projected onto new sample data analyzed via RNA-sequencing. We then detail procedures for reverse engineering a continuum score for samples analyzed via DNA methylation microarray using a random forest classifier.

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.

Activation of the RAS-RAF-MEK-ERK signaling pathway is thought to be the key driver of pediatric low-grade astrocytoma (PLGA) growth. Sorafenib is a multikinase inhibitor targeting BRAF, VEGFR, PDGFR, and c-kit. This multicenter phase II study was conducted to determine the response rate to sorafenib in patients with recurrent or progressive PLGA.

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.

High morbidity and mortality are common traits of malignant tumours and identification of the cells responsible is a focus of on-going research. Many studies are now reporting the use of antibodies specific to Clusters of Differentiation (CD) cell surface antigens to identify tumour-initiating cell (TIC) populations in neural tumours. Medulloblastoma is one of the most common malignant brain tumours in children and despite a considerable amount of research investigating this tumour, the identity of the TICs, and the means by which such cells can be targeted remain largely unknown. Current prognostication and stratification of medulloblastoma using clinical factors, histology and genetic profiling have classified this tumour into four main subgroups: WNT, Sonic hedgehog (SHH), Group 3 and Group 4. Of these subgroups, SHH remains one of the most studied tumour groups due to the ability to model medulloblastoma formation through targeted deletion of the Shh pathway inhibitor Patched1 (Ptch1). Here we sought to utilise CD antibody expression to identify and isolate TIC populations in Ptch1 deleted medulloblastoma, and determine if these antibodies can help classify the identity of human medulloblastoma subgroups. Using a fluorescence-activated cell sorted (FACS) CD antibody panel, we identified CD24 as a marker of TICs in Ptch1 deleted medulloblastoma. CD24 expression was not correlated with markers of astrocytes or oligodendrocytes, but co-labelled with markers of neural progenitor cells. In conjunction with CD15, proliferating CD24+/CD15+ granule cell precursors (GCPs) were identified as a TIC population in Ptch1 deleted medulloblastoma. On human medulloblastoma, CD24 was found to be highly expressed on Group 3, Group 4 and SHH subgroups compared with the WNT subgroup, which was predominantly positive for CD15, suggesting CD24 is an important marker of non-WNT medulloblastoma initiating cells and a potential therapeutic target in human medulloblastoma. This study reports the use of CD24 and CD15 to isolate a GCP-like TIC population in Ptch1 deleted medulloblastoma, and suggests CD24 expression as a marker to help stratify human WNT tumours from other medulloblastoma subgroups.

Low-grade gliomas (LGGs) account for about a third of all brain tumours in children. We conducted a detailed study of DNA methylation and gene expression to improve our understanding of the biology of pilocytic and diffuse astrocytomas. Pilocytic astrocytomas were found to have a distinctive signature at 315 CpG sites, of which 312 were hypomethylated and 3 were hypermethylated. Genomic analysis revealed that 182 of these sites are within annotated enhancers. The signature was not present in diffuse astrocytomas, or in published profiles of other brain tumours and normal brain tissue. The AP-1 transcription factor was predicted to bind within 200 bp of a subset of the 315 differentially methylated CpG sites; the AP-1 factors, FOS and FOSL1 were found to be up-regulated in pilocytic astrocytomas. We also analysed splice variants of the AP-1 target gene, CCND1, which encodes cell cycle regulator cyclin D1. CCND1a was found to be highly expressed in both pilocytic and diffuse astrocytomas, but diffuse astrocytomas have far higher expression of the oncogenic variant, CCND1b. These findings highlight novel genetic and epigenetic differences between pilocytic and diffuse astrocytoma, in addition to well-described alterations involving BRAF, MYB and FGFR1.

Blastemal histology in chemotherapy-treated pediatric Wilms tumors (nephroblastoma) is associated with adverse prognosis. To uncover the underlying tumor biology and find therapeutic leads for this subgroup, we analyzed 58 blastemal type Wilms tumors by exome and transcriptome sequencing and validated our findings in a large replication cohort. Recurrent mutations included a hotspot mutation (Q177R) in the homeo-domain of SIX1 and SIX2 in tumors with high proliferative potential (18.1% of blastemal cases); mutations in the DROSHA/DGCR8 microprocessor genes (18.2% of blastemal cases); mutations in DICER1 and DIS3L2; and alterations in IGF2, MYCN, and TP53, the latter being strongly associated with dismal outcome. DROSHA and DGCR8 mutations strongly altered miRNA expression patterns in tumors, which was functionally validated in cell lines expressing mutant DROSHA.

Metabolic adaptation is essential for cell survival during nutrient deprivation. We report that eukaryotic elongation factor 2 kinase (eEF2K), which is activated by AMP-kinase (AMPK), confers cell survival under acute nutrient depletion by blocking translation elongation. Tumor cells exploit this pathway to adapt to nutrient deprivation by reactivating the AMPK-eEF2K axis. Adaptation of transformed cells to nutrient withdrawal is severely compromised in cells lacking eEF2K. Moreover, eEF2K knockdown restored sensitivity to acute nutrient deprivation in highly resistant human tumor cell lines. In vivo, overexpression of eEF2K rendered murine tumors remarkably resistant to caloric restriction. Expression of eEF2K strongly correlated with overall survival in human medulloblastoma and glioblastoma multiforme. Finally, C. elegans strains deficient in efk-1, the eEF2K ortholog, were severely compromised in their response to nutrient depletion. Our data highlight a conserved role for eEF2K in protecting cells from nutrient deprivation and in conferring tumor cell adaptation to metabolic stress. PAPERCLIP:

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.

Spatial heterogeneity of transcriptional and genetic markers between physically isolated biopsies of a single tumor poses major barriers to the identification of biomarkers and the development of targeted therapies that will be effective against the entire tumor. We analyzed the spatial heterogeneity of multiregional biopsies from 35 patients, using a combination of transcriptomic and genomic profiles. Medulloblastomas (MBs), but not high-grade gliomas (HGGs), demonstrated spatially homogeneous transcriptomes, which allowed for accurate subgrouping of tumors from a single biopsy. Conversely, somatic mutations that affect genes suitable for targeted therapeutics demonstrated high levels of spatial heterogeneity in MB, malignant glioma, and renal cell carcinoma (RCC). Actionable targets found in a single MB biopsy were seldom clonal across the entire tumor, which brings the efficacy of monotherapies against a single target into question. Clinical trials of targeted therapies for MB should first ensure the spatially ubiquitous nature of the target mutation.

The HERBY trial was a phase II open-label, randomized, multicenter trial evaluating bevacizumab (BEV) in addition to temozolomide/radiotherapy in patients with newly diagnosed non-brainstem high-grade glioma (HGG) between the ages of 3 and 18 years. We carried out comprehensive molecular analysis integrated with pathology, radiology, and immune profiling. In post-hoc subgroup analysis, hypermutator tumors (mismatch repair deficiency and somatic POLE/POLD1 mutations) and those biologically resembling pleomorphic xanthoastrocytoma ([PXA]-like, driven by BRAF_V600E or NF1 mutation) had significantly more CD8+ tumor-infiltrating lymphocytes, and longer survival with the addition of BEV. Histone H3 subgroups (hemispheric G34R/V and midline K27M) had a worse outcome and were immune cold. Future clinical trials will need to take into account the diversity represented by the term "HGG" in the pediatric population.

Recurrent mutations affecting the histone H3.3 residues Lys27 or indirectly Lys36 are frequent drivers of pediatric high-grade gliomas (over 30% of HGGs). To identify additional driver mutations in HGGs, we investigated a cohort of 60 pediatric HGGs using whole-exome sequencing (WES) and compared them to 543 exomes from non-cancer control samples. We identified mutations in SETD2, a H3K36 trimethyltransferase, in 15% of pediatric HGGs, a result that was genome-wide significant (FDR = 0.029). Most SETD2 alterations were truncating mutations. Sequencing the gene in this cohort and another validation cohort (123 gliomas from all ages and grades) showed SETD2 mutations to be specific to high-grade tumors affecting 15% of pediatric HGGs (11/73) and 8% of adult HGGs (5/65) while no SETD2 mutations were identified in low-grade diffuse gliomas (0/45). Furthermore, SETD2 mutations were mutually exclusive with H3F3A mutations in HGGs (P = 0.0492) while they partly overlapped with IDH1 mutations (4/14), and SETD2-mutant tumors were found exclusively in the cerebral hemispheres (P = 0.0055). SETD2 is the only H3K36 trimethyltransferase in humans, and SETD2-mutant tumors showed a substantial decrease in H3K36me3 levels (P < 0.001), indicating that the mutations are loss-of-function. These data suggest that loss-of-function SETD2 mutations occur in older children and young adults and are specific to HGG of the cerebral cortex, similar to the H3.3 G34R/V and IDH mutations. Taken together, our results suggest that mutations disrupting the histone code at H3K36, including H3.3 G34R/V, IDH1 and/or SETD2 mutations, are central to the genesis of hemispheric HGGs in older children and young adults.

Medulloblastoma is the most common malignant brain tumor in children. Despite recent improvements in cure rates, prediction of disease outcome remains a major challenge and survivors suffer from serious therapy-related side-effects. Recent data showed that patients with WNT-activated tumors have a favorable prognosis, suggesting that these patients could be treated less intensively, thereby reducing the side-effects. This illustrates the potential benefits of a robust classification of medulloblastoma patients and a detailed knowledge of associated biological mechanisms.

In cancer, recurrent somatic single-nucleotide variants-which are rare in most paediatric cancers-are confined largely to protein-coding genes1-3. Here we report highly recurrent hotspot mutations (r.3A>G) of U1 spliceosomal small nuclear RNAs (snRNAs) in about 50% of Sonic hedgehog (SHH) medulloblastomas. These mutations were not present across other subgroups of medulloblastoma, and we identified these hotspot mutations in U1 snRNA in only <0.1% of 2,442 cancers, across 36 other tumour types. The mutations occur in 97% of adults (subtype SHHδ) and 25% of adolescents (subtype SHHα) with SHH medulloblastoma, but are largely absent from SHH medulloblastoma in infants. The U1 snRNA mutations occur in the 5' splice-site binding region, and snRNA-mutant tumours have significantly disrupted RNA splicing and an excess of 5' cryptic splicing events. Alternative splicing mediated by mutant U1 snRNA inactivates tumour-suppressor genes (PTCH1) and activates oncogenes (GLI2 and CCND2), and represents a target for therapy. These U1 snRNA mutations provide an example of highly recurrent and tissue-specific mutations of a non-protein-coding gene in cancer.

The third-generation tyrosine kinase inhibitor (TKI), osimertinib, has revolutionized the treatment of patients with non-small cell lung carcinoma with epidermal growth factor receptor (EGFR)-activating mutation, and resistant to first- and second-generation TKIs. Osimertinib is now also proposed as a first-line therapy, thus extending the scope of applications in lung oncology. Personalized medicine approaches are still necessary to monitor if patients are exposed to adequate concentrations of osimertinib during their treatment. It would also help to understand the appearance of new resistances in patients after several months of dosing with osimertinib. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) is currently the gold standard for the quantification of drugs in plasma enabling pharmacokinetic analyses and patient monitoring. In the present study, we propose an alternative to LC-MS/MS methods for the rapid and sensitive quantification of osimertinib in plasma using matrix-assisted laser desorption/ionization (MALDI) -MS. The presented assay requires only 3 min per sample for their preparation, analysis, and data extraction, and less than 3 h for quantification. A lower limit of quantification (LLOQ) of 5 ng/mL in plasma was retrieved. The method was fully validated, following the guidelines of the US Food and Drug Administration (FDA) and the European Medicines Agency (EMA) for bioanalytical method validation. The present developments prove the importance to consider alternative MS assays for time-efficient quantification of small molecule inhibitors in plasma in the context of personalized medicine for targeted therapies.