Background: Fusion genes play an important role in the tumorigenesis of many cancers. Next-generation sequencing (NGS) technologies have been successfully applied in fusion gene detection for the last several years, and a number of NGS-based tools have been developed for identifying fusion genes during this period. Most fusion gene detection tools based on RNA-seq data report a large number of candidates (mostly false positives), making it hard to prioritize candidates for experimental validation and further analysis. Selection of reliable fusion genes for downstream analysis becomes very important in cancer research. We therefore developed confFuse, a scoring algorithm to reliably select high-confidence fusion genes which are likely to be biologically relevant. Results: confFuse takes multiple parameters into account in order to assign each fusion candidate a confidence score, of which score ≥8 indicates high-confidence fusion gene predictions. These parameters were manually curated based on our experience and on certain structural motifs of fusion genes. Compared with alternative tools, based on 96 published RNA-seq samples from different tumor entities, our method can significantly reduce the number of fusion candidates (301 high-confidence from 8,083 total predicted fusion genes) and keep high detection accuracy (recovery rate 85.7%). Validation of 18 novel, high-confidence fusions detected in three breast tumor samples resulted in a 100% validation rate. Conclusions: confFuse is a novel downstream filtering method that allows selection of highly reliable fusion gene candidates for further downstream analysis and experimental validations. confFuse is available at https://github.com/Zhiqin-HUANG/confFuse.

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.

Although osteosarcoma (OS) is a rare cancer, it is the most common primary malignant bone tumor in children and adolescents. BRCAness is a phenotypical trait in tumors with a defect in homologous recombination repair, resembling tumors with inactivation of BRCA1/2, rendering these tumors sensitive to poly (ADP)-ribose polymerase inhibitors (PARPi). Recently, OS was shown to exhibit molecular features of BRCAness. Our goal was to develop a method complementing existing genomic methods to aid clinical decision making on administering PARPi in OS patients. OS samples with DNA-methylation data were divided to BRCAness-positive and negative groups based on the degree of their genomic instability (n = 41). Methylation probes were ranked according to decreasing variance difference between two groups. The top 2000 probes were selected for training and cross-validation of the random forest algorithm. Two-thirds of available OS RNA-Seq samples (n = 17) from the top and bottom of the sample list ranked according to genome instability score were subjected to differential expression and, subsequently, to gene set enrichment analysis (GSEA). The combined accuracy of trained random forest was 85% and the average area under the ROC curve (AUC) was 0.95. There were 449 upregulated and 1,079 downregulated genes in the BRCAness-positive group (fdr < 0.05). GSEA of upregulated genes detected enrichment of DNA replication and mismatch repair and homologous recombination signatures (FWER < 0.05). Validation of the BRCAness classifier with an independent OS set (n = 20) collected later in the course of study showed AUC of 0.87 with an accuracy of 90%. GSEA signatures computed for this test set were matching the ones observed in the training set enrichment analysis. In conclusion, we developed a new classifier based on DNA-methylation patterns that detects BRCAness in OS samples with high accuracy. GSEA identified genome instability signatures. Machine-learning and gene expression approaches add new epigenomic and transcriptomic aspects to already established genomic methods for evaluation of BRCAness in osteosarcoma and can be extended to cancers characterized by genome instability.

Pleomorphic xanthoastrocytoma (PXA) in its classic manifestation exhibits distinct morphological features and is assigned to CNS WHO grade 2 or grade 3. Distinction from glioblastoma variants and lower grade glial and glioneuronal tumors is a common diagnostic challenge. We compared a morphologically defined set of PXA (histPXA) with an independent set, defined by DNA methylation analysis (mcPXA). HistPXA encompassed 144 tumors all subjected to DNA methylation array analysis. Sixty-two histPXA matched to the methylation class mcPXA. These were combined with the cases that showed the mcPXA signature but had received a histopathological diagnosis other than PXA. This cohort constituted a set of 220 mcPXA. Molecular and clinical parameters were analyzed in these groups. Morphological parameters were analyzed in a subset of tumors with FFPE tissue available. HistPXA revealed considerable heterogeneity in regard to methylation classes, with methylation classes glioblastoma and ganglioglioma being the most frequent mismatches. Similarly, the mcPXA cohort contained tumors of diverse histological diagnoses, with glioblastoma constituting the most frequent mismatch. Subsequent analyses demonstrated the presence of canonical pTERT mutations to be associated with unfavorable prognosis among mcPXA. Based on these data, we consider the tumor type PXA to be histologically more varied than previously assumed. Histological approach to diagnosis will predominantly identify cases with the established archetypical morphology. DNA methylation analysis includes additional tumors in the tumor class PXA that share similar DNA methylation profile but lack the typical morphology of a PXA. DNA methylation analysis also assist in separating other tumor types with morphologic overlap to PXA. Our data suggest the presence of canonical pTERT mutations as a robust indicator for poor prognosis in methylation class PXA.

Primitive neuroectodermal tumors of the central nervous system (CNS-PNETs) are highly malignant neoplasms posing diagnostic challenge due to a lack of defining molecular markers. CNS neuroblastoma with forkhead box R2 (FOXR2) activation (CNS_NBL) emerged as a distinct pediatric brain tumor entity from a pool previously diagnosed as primitive neuroectodermal tumors of the central nervous system (CNS-PNETs). Current standard of identifying CNS_NBL relies on molecular analysis. We set out to establish immunohistochemical markers allowing safely distinguishing CNS_NBL from morphological mimics. To this aim we analyzed a series of 84 brain tumors institutionally diagnosed as CNS-PNET. As expected, epigenetic analysis revealed different methylation groups corresponding to the (1) CNS-NBL (24%), (2) glioblastoma IDH wild-type subclass H3.3 G34 (26%), (3) glioblastoma IDH wild-type subclass MYCN (21%) and (4) ependymoma with RELA_C11orf95 fusion (29%) entities. Transcriptome analysis of this series revealed a set of differentially expressed genes distinguishing CNS_NBL from its mimics. Based on RNA-sequencing data we established SOX10 and ANKRD55 expression as genes discriminating CNS_NBL from other tumors exhibiting CNS-PNET. Immunohistochemical detection of combined expression of SOX10 and ANKRD55 clearly identifies CNS_NBL discriminating them to other hemispheric CNS neoplasms harboring "PNET-like" microscopic appearance. Owing the rarity of CNS_NBL, a confirmation of the elaborated diagnostic IHC algorithm will be necessary in prospective patient series.

Atypical teratoid/rhabdoid tumor (AT/RT) is a highly malignant brain tumor in infants that is characterized by loss of nuclear expression of SMARCB1 or SMARCA4 proteins. Recent studies show that AT/RTs comprise three molecular subgroups, namely AT/RT-TYR, AT/RT-MYC and AT/RT-SHH. The subgroups show distinct expression patterns of genes involved in ciliogenesis, however, little is known about the functional roles of primary cilia in the biology of AT/RT. Here, we show that primary cilia are present across all AT/RT subgroups with specific enrichment in AT/RT-TYR patient samples. Furthermore, we demonstrate that primary ciliogenesis contributes to AT/RT biology in vitro and in vivo. Specifically, we observed a significant decrease in proliferation and clonogenicity following disruption of primary ciliogenesis in AT/RT cell line models. Additionally, apoptosis was significantly increased via the induction of STAT1 and DR5 signaling, as detected by proteogenomic profiling. In a Drosophila model of SMARCB1 deficiency, concomitant knockdown of several cilia-associated genes resulted in a substantial shift of the lethal phenotype with more than 20% of flies reaching adulthood. We also found significantly extended survival in an orthotopic xenograft mouse model of AT/RT upon disruption of primary ciliogenesis. Taken together, our findings indicate that primary ciliogenesis or its downstream signaling contributes to the aggressiveness of AT/RT and, therefore, may constitute a novel therapeutic target.

Ewing sarcoma (EwS) is characterized by EWSR1-ETS fusion transcription factors converting polymorphic GGAA microsatellites (mSats) into potent neo-enhancers. Although the paucity of additional mutations makes EwS a genuine model to study principles of cooperation between dominant fusion oncogenes and neo-enhancers, this is impeded by the limited number of well-characterized models. Here we present the Ewing Sarcoma Cell Line Atlas (ESCLA), comprising whole-genome, DNA methylation, transcriptome, proteome, and chromatin immunoprecipitation sequencing (ChIP-seq) data of 18 cell lines with inducible EWSR1-ETS knockdown. The ESCLA shows hundreds of EWSR1-ETS-targets, the nature of EWSR1-ETS-preferred GGAA mSats, and putative indirect modes of EWSR1-ETS-mediated gene regulation, converging in the duality of a specific but plastic EwS signature. We identify heterogeneously regulated EWSR1-ETS-targets as potential prognostic EwS biomarkers. Our freely available ESCLA (http://r2platform.com/escla/) is a rich resource for EwS research and highlights the power of comprehensive datasets to unravel principles of heterogeneous gene regulation by chimeric transcription factors.

Analysis of mutational signatures can reveal underlying molecular mechanisms of the processes that have imprinted the somatic mutations found in cancer genomes. Here, we analyze single base substitutions and small insertions and deletions in pediatric cancers encompassing 785 whole-genome sequenced tumors from 27 molecularly defined cancer subtypes. We identified only a small number of mutational signatures active in pediatric cancers, compared with previously analyzed adult cancers. Further, we report a significant difference in the proportion of pediatric tumors showing homologous recombination repair defect signatures compared with previous analyses. In pediatric leukemias, we identified an indel signature, not previously reported, characterized by long insertions in nonrepeat regions, affecting mainly intronic and intergenic regions, but also exons of known cancer genes. We provide a systematic overview of COSMIC v.3 mutational signatures active across pediatric cancers, which is highly relevant for understanding tumor biology and enabling future research in defining biomarkers of treatment response.

Glioneuronal tumors are a heterogenous group of CNS neoplasms that can be challenging to accurately diagnose. Molecular methods are highly useful in classifying these tumors-distinguishing precise classes from their histological mimics and identifying previously unrecognized types of tumors. Using an unsupervised visualization approach of DNA methylation data, we identified a novel group of tumors (n = 20) that formed a cluster separate from all established CNS tumor types. Molecular analyses revealed ATRX alterations (in 16/16 cases by DNA sequencing and/or immunohistochemistry) as well as potentially targetable gene fusions involving receptor tyrosine-kinases (RTK; mostly NTRK1-3) in all of these tumors (16/16; 100%). In addition, copy number profiling showed homozygous deletions of CDKN2A/B in 55% of cases. Histological and immunohistochemical investigations revealed glioneuronal tumors with isomorphic, round and often condensed nuclei, perinuclear clearing, high mitotic activity and microvascular proliferation. Tumors were mainly located supratentorially (84%) and occurred in patients with a median age of 19 years. Survival data were limited (n = 18) but point towards a more aggressive biology as compared to other glioneuronal tumors (median progression-free survival 12.5 months). Given their molecular characteristics in addition to anaplastic features, we suggest the term glioneuronal tumor with ATRX alteration, kinase fusion and anaplastic features (GTAKA) to describe these tumors. In summary, our findings highlight a novel type of glioneuronal tumor driven by different RTK fusions accompanied by recurrent alterations in ATRX and homozygous deletions of CDKN2A/B. Targeted approaches such as NTRK inhibition might represent a therapeutic option for patients suffering from these tumors.

The mechanisms underlying T-ALL relapse remain essentially unknown. Multilevel-omics in 38 matched pairs of initial and relapsed T-ALL revealed 18 (47%) type-1 (defined by being derived from the major ancestral clone) and 20 (53%) type-2 relapses (derived from a minor ancestral clone). In both types of relapse, we observed known and novel drivers of multidrug resistance including MDR1 and MVP, NT5C2 and JAK-STAT activators. Patients with type-1 relapses were specifically characterized by IL7R upregulation. In remarkable contrast, type-2 relapses demonstrated (1) enrichment of constitutional cancer predisposition gene mutations, (2) divergent genetic and epigenetic remodeling, and (3) enrichment of somatic hypermutator phenotypes, related to BLM, BUB1B/PMS2 and TP53 mutations. T-ALLs that later progressed to type-2 relapses exhibited a complex subclonal architecture, unexpectedly, already at the time of initial diagnosis. Deconvolution analysis of ATAC-Seq profiles showed that T-ALLs later developing into type-1 relapses resembled a predominant immature thymic T-cell population, whereas T-ALLs developing into type-2 relapses resembled a mixture of normal T-cell precursors. In sum, our analyses revealed fundamentally different mechanisms driving either type-1 or type-2 T-ALL relapse and indicate that differential capacities of disease evolution are already inherent to the molecular setup of the initial leukemia.

Actinomycin D is a potent cytotoxic drug against pediatric (and other) tumors that is thought to barely cross the blood-brain barrier. To evaluate its potential applicability for the treatment of patients with central nervous system (CNS) tumors, we established a cerebral microdialysis model in freely moving mice and investigated its CNS disposition by quantifying actinomycin D in cerebral microdialysate, brain tissue homogenate, and plasma. For this purpose, we developed and validated an ultraperformance liquid chromatography-tandem mass spectrometry assay suitable for ultra-sensitive quantification of actinomycin D in the pertinent biological matrices in micro-samples of only 20 µL, with a lower limit of quantification of 0.05 ng/mL. In parallel, we confirmed actinomycin D as a substrate of P-glycoprotein (P-gp) in in vitro experiments. Two hours after intravenous administration of 0.5 mg/kg, actinomycin D reached total brain tissue concentrations of 4.1 ± 0.7 ng/g corresponding to a brain-to-plasma ratio of 0.18 ± 0.03, while it was not detectable in intracerebral microdialysate. This tissue concentration exceeds the concentrations of actinomycin D that have been shown to be effective in in vitro experiments. Elimination of the drug from brain tissue was substantially slower than from plasma, as shown in a brain-to-plasma ratio of approximately 0.53 after 22 h. Because actinomycin D reached potentially effective concentrations in brain tissue in our experiments, the drug should be further investigated as a therapeutic agent in potentially susceptible CNS malignancies, such as ependymoma.

Cerebellar granule cells (GCs) are cells which comprise over 50% of the neurons in the entire nervous system. GCs enable the cerebellum to properly regulate motor coordination, learning, and consolidation, in addition to cognition, emotion and language. During GC development, maternal GC progenitors (GCPs) divide to produce not only postmitotic GCs but also sister GCPs. However, the molecular machinery for regulating the proportional production of distinct sister cell types from seemingly uniform GCPs is not yet fully understood. Here we report that Notch signaling creates a distinction between GCPs and leads to their proportional differentiation in mice. Among Notch-related molecules, Notch1, Notch2, Jag1, and Hes1 are prominently expressed in GCPs. In vivo monitoring of Hes1-promoter activities showed the presence of two types of GCPs, Notch-signaling ON and OFF, in the external granule layer (EGL). Single-cell RNA sequencing (scRNA-seq) and in silico analyses indicate that ON-GCPs have more proliferative and immature properties, while OFF-GCPs have opposite characteristics. Overexpression as well as knock-down (KD) experiments using in vivo electroporation showed that NOTCH2 and HES1 are involved cell-autonomously to suppress GCP differentiation by inhibiting NEUROD1 expression. In contrast, JAG1-expressing cells non-autonomously upregulated Notch signaling activities via NOTCH2-HES1 in surrounding GCPs, eventually suppressing their differentiation. These findings suggest that Notch signaling results in the proportional generation of two types of cells, immature and differentiating GCPs, which contributes to the well-organized differentiation of GCs.

More than 60% of supratentorial ependymomas harbor a ZFTA-RELA (ZRfus) gene fusion (formerly C11orf95-RELA). To study the biology of ZRfus, we developed an autochthonous mouse tumor model using in utero electroporation (IUE) of the embryonic mouse brain. Integrative epigenomic and transcriptomic mapping was performed on IUE-driven ZRfus tumors by CUT&RUN, chromatin immunoprecipitation sequencing, assay for transposase-accessible chromatin sequencing, and RNA sequencing and compared with human ZRfus-driven ependymoma. In addition to direct canonical NFκB pathway activation, ZRfus dictates a neoplastic transcriptional program and binds to thousands of unique sites across the genome that are enriched with PLAGL family transcription factor (TF) motifs. ZRfus activates gene expression programs through recruitment of transcriptional coactivators (Brd4, Ep300, Cbp, Pol2) that are amenable to pharmacologic inhibition. Downstream ZRfus target genes converge on developmental programs marked by PLAGL TF proteins, and activate neoplastic programs enriched in Mapk, focal adhesion, and gene imprinting networks. SIGNIFICANCE: Ependymomas are aggressive brain tumors. Although drivers of supratentorial ependymoma (ZFTA- and YAP1-associated gene fusions) have been discovered, their functions remain unclear. Our study investigates the biology of ZFTA-RELA-driven ependymoma, specifically mechanisms of transcriptional deregulation and direct downstream gene networks that may be leveraged for potential therapeutic testing.This article is highlighted in the In This Issue feature, p. 2113.

Matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS) has rarely been used in the field of therapeutic drug monitoring, partly because of the complexity of the ionization processes between the compounds to be quantified and the many MALDI matrices available. The development of a viable MALDI-MS method that meets regulatory guidelines for bioanalytical method validation requires prior knowledge of the suitability of (i) the MALDI matrix with the analyte class and properties for ionization, (ii) the crystallization properties of the MALDI matrix with automation features, and (iii) the MS instrumentation used to achieve sensitive and specific measurements in order to determine low pharmacological drug concentrations in biological matrices. In the present hybrid article/white paper, we review the developments required for the establishment of MALDI-MS assays for the quantification of drugs in tissues and plasma, illustrated with concrete results for the different steps. We summarize the necessary parameters that need to be controlled for the successful development of fully validated MALDI-MS methods according to regulatory authorities, as well as currently unsolved problems and promising ways to address them. Finally, we propose an expert opinion on future perspectives and needs in order to establish MALDI-MS as a universal method for therapeutic drug monitoring.

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.

Pituicytoma (PITUI), granular cell tumor (GCT), and spindle cell oncocytoma (SCO) are rare tumors of the posterior pituitary. Histologically, they may be challenging to distinguish and have been proposed to represent a histological spectrum of a single entity. We performed targeted next-generation sequencing, DNA methylation profiling, and copy number analysis on 47 tumors (14 PITUI; 12 GCT; 21 SCO) to investigate molecular features and explore possibilities of clinically meaningful tumor subclassification. We detected two main epigenomic subgroups by unsupervised clustering of DNA methylation data, though the overall methylation differences were subtle. The largest group (n = 23) contained most PITUIs and a subset of SCOs and was enriched for pathogenic mutations within genes in the MAPK/PI3K pathways (12/17 [71%] of sequenced tumors: FGFR1 (3), HRAS (3), BRAF (2), NF1 (2), CBL (1), MAP2K2 (1), PTEN (1)) and two with accompanying TERT promoter mutation. The second group (n = 16) contained most GCTs and a subset of SCOs, all of which mostly lacked identifiable genetic drivers. Outcome analysis demonstrated that the presence of chromosomal imbalances was significantly associated with reduced progression-free survival especially within the combined PITUI and SCO group (p = 0.031). In summary, we observed only subtle DNA methylation differences between posterior pituitary tumors, indicating that these tumors may be best classified as subtypes of a single entity. Nevertheless, our data indicate differences in mutation patterns and clinical outcome. For a clinically meaningful subclassification, we propose a combined histo-molecular approach into three subtypes: one subtype is defined by granular cell histology, scarcity of identifiable oncogenic mutations, and favorable outcome. The other two subtypes have either SCO or PITUI histology but are segregated by chromosomal copy number profile into a favorable group (no copy number changes) and a less favorable group (copy number imbalances present). Both of the latter groups have recurrent MAPK/PI3K genetic alterations that represent potential therapeutic targets.

Forkhead box R2 (FOXR2) is a forkhead transcription factor located on the X chromosome whose expression is normally restricted to the testis. In this study, we performed a pan-cancer analysis of FOXR2 activation across more than 10,000 adult and pediatric cancer samples and found FOXR2 to be aberrantly upregulated in 70% of all cancer types and 8% of all individual tumors. The majority of tumors (78%) aberrantly expressed FOXR2 through a previously undescribed epigenetic mechanism that involves hypomethylation of a novel promoter, which was functionally validated as necessary for FOXR2 expression and proliferation in FOXR2-expressing cancer cells. FOXR2 promoted tumor growth across multiple cancer lineages and co-opted ETS family transcription circuits across cancers. Taken together, this study identifies FOXR2 as a potent and ubiquitous oncogene that is epigenetically activated across the majority of human cancers. The identification of hijacking of ETS transcription circuits by FOXR2 extends the mechanisms known to active ETS transcription factors and highlights how transcription factor families cooperate to enhance tumorigenesis.

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.

Ependymal tumors across age groups are currently classified and graded solely by histopathology. It is, however, commonly accepted that this classification scheme has limited clinical utility based on its lack of reproducibility in predicting patients' outcome. We aimed at establishing a uniform molecular classification using DNA methylation profiling. Nine molecular subgroups were identified in a large cohort of 500 tumors, 3 in each anatomical compartment of the CNS, spine, posterior fossa, supratentorial. Two supratentorial subgroups are characterized by prototypic fusion genes involving RELA and YAP1, respectively. Regarding clinical associations, the molecular classification proposed herein outperforms the current histopathological classification and thus might serve as a basis for the next World Health Organization classification of CNS tumors.