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.

Ependymomas in children can arise throughout all compartments of the central nervous system (CNS). Highly malignant paediatric ependymoma subtypes are Group A tumours of the posterior fossa (PF-EPN-A) and RELA-fusion positive (ST-EPN-RELA) tumours in the supratentorial compartment. It was repeatedly reported in smaller series that accumulation of p53 is frequently observed in ependymomas and that immunohistochemical staining correlates with poor clinical outcome, while TP53 mutations are rare. Our TP53 mutation analysis of 130 primary ependymomas identified a mutation rate of only 3%. Immunohistochemical analysis of 398 ependymomas confirmed previous results correlating the accumulation of p53 with inferior outcome. Among the p53-positive ependymomas, the vast majority exhibited a RELA fusion leading to the hypothesis that p53 inactivation might be linked to RELA positivity.In order to assess the potential of p53 reactivation through MDM2 inhibition in ependymoma, we evaluated the effects of Actinomycin-D and Nutlin-3 treatment in two preclinical ependymoma models representing the high-risk subtypes PF-EPN-A and ST-EPN-RELA. The IC-50 of the agent as determined by metabolic activity assays was in the lower nano-molar range (0.2-0.7 nM). Transcriptome analyses of high-dose (100 nM), low-dose (5 nM) and non-treated cells revealed re-expression of p53 dependent genes including p53 upregulated modulator of apoptosis (PUMA) after low-dose treatment. At the protein level, we validated the Actinomycin-D induced upregulation of PUMA, and of p53 interaction partners MDM2 and p21. Proapoptotic effects of low-dose application of the agent were confirmed by flow cytometry. Thus, Actinomycin-D could constitute a promising therapeutic option for ST-EPN-RELA ependymoma patients, whose tumours frequently exhibit p53 inactivation.

Gain-of-function mutations in histone 3 (H3) variants are found in a substantial proportion of pediatric high-grade gliomas (pHGG), often in association with TP53 loss and platelet-derived growth factor receptor alpha (PDGFRA) amplification. Here, we describe a somatic mouse model wherein H3.3K27M and Trp53 loss alone are sufficient for neoplastic transformation if introduced in utero. H3.3K27M-driven lesions are clonal, H3K27me3 depleted, Olig2 positive, highly proliferative, and diffusely spreading, thus recapitulating hallmark molecular and histopathological features of pHGG. Addition of wild-type PDGFRA decreases latency and increases tumor invasion, while ATRX knockdown is associated with more circumscribed tumors. H3.3K27M-tumor cells serially engraft in recipient mice, and preliminary drug screening reveals mutation-specific vulnerabilities. Overall, we provide a faithful H3.3K27M-pHGG model which enables insights into oncohistone pathogenesis and investigation of future therapies.

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.

Ten-Eleven Translocation (TET) proteins have recently come to light as important epigenetic regulators conserved in multicellular organisms. TET knockdown studies in rodents have highlighted the critical role of these proteins for proper brain development and function. Mutations in mammalian mTET proteins and mTET2 specifically are frequent and deregulated in leukaemia and glioma respectively. Accordingly, we examined the role of mTET2 in tumorigenesis in larval haemocytes and adult heads in Drosophila melanogaster. Our findings showed that expression of mutant and wild type mTET2 resulted in general phenotypic defects in adult flies and accumulation of abdominal melanotic masses. Notably, flies with mTET2-R43G mutation at the N-terminus of mTET2 exhibited locomotor and circadian behavioural deficits, as well as reduced lifespan. Flies with mTET2-R1261C mutation in the catalytic domain, a common mutation in acute myeloid leukaemia (AML), displayed alterations affecting haemocyte haemostasis. Using transcriptomic approach, we identified upregulated immune genes in fly heads that were not exclusive to TET2 mutants but also found in wild type mTET2 flies. Furthermore, inhibiting expression of genes that were found to be deregulated in mTET2 mutants, such as those involved in immune pathways, autophagy, and transcriptional regulation, led to a rescue in fly survival, behaviour, and hemocyte number. This study identifies the transcriptomic profile of wild type mTET2 versus mTET2 mutants (catalytic versus non-catalytic) with indications of TET2 role in normal central nervous system (CNS) function and innate immunity.

Accurate identification of brain tumor molecular subgroups is increasingly important. We aimed to establish the most accurate and reproducible ependymoma subgroup biomarker detection techniques, across 147 cases from International Society of Pediatric Oncology (SIOP) Ependymoma II trial participants, enrolled in the pan-European "Biomarkers of Ependymoma in Children and Adolescents (BIOMECA)" study.

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.

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.

Loss of nuclear SMARCB1 (INI1/hSNF5/BAF47) protein expression due to biallelic mutations of the SMARCB1 tumor suppressor gene is a hallmark of atypical teratoid/rhabdoid tumors (ATRT), but the presence of cytoplasmic SMARCB1 protein in these tumors has not yet been described. In a series of 102 primary ATRT, distinct cytoplasmic SMARCB1 staining on immunohistochemistry was encountered in 19 cases (19%) and was highly over-represented in cases showing pathogenic sequence variants leading to truncation or mutation of the C-terminal part of SMARCB1 (15/19 vs. 4/83; Chi-square: 56.04, p = 1.0E-10) and, related to this, in tumors of the molecular subgroup ATRT-TYR (16/36 vs. 3/66; Chi-square: 24.47, p = 7.6E-7). Previous reports have indicated that while SMARCB1 lacks a bona fide nuclear localization signal, it harbors a masked nuclear export signal (NES) and that truncation of the C-terminal region results in unmasking of this NES leading to cytoplasmic localization. To determine if cytoplasmic localization found in ATRT is due to unmasking of NES, we generated GFP fusions of one of the SMARCB1 truncating mutations (p.Q318X) found in the tumors along with a p.L266A mutation, which was shown to disrupt the interaction of SMARCB1-NES with exportin-1. We found that while the GFP-SMARCB1(Q318X) mutant localized to the cytoplasm, the double mutant GFP-SMARCB1(Q318X;L266A) localized to the nucleus, confirming NES requirement for cytoplasmic localization. Furthermore, cytoplasmic SMARCB1(Q318X) was unable to cause senescence as determined by morphological observations and by senescence-associated β-galactosidase assay, while nuclear SMARCB1(Q318X;L266A) mutant regained this function. Selinexor, a selective exportin-1 inhibitor, was effective in inhibiting the nuclear export of SMARCB1(Q318X) and caused rapid cell death in rhabdoid tumor cells. In conclusion, inhibition of nuclear export restores nuclear localization and residual tumor suppressor function of truncated SMARCB1. Therapies aimed at preventing nuclear export of mutant SMARCB1 protein may represent a promising targeted therapy in ATRT harboring truncating C-terminal SMARCB1 mutations.

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.

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.

Polycomb Repressive Complex 2 (PRC2)-mediated histone H3K27 tri-methylation (H3K27me3) recruits canonical PRC1 (cPRC1) to maintain heterochromatin. In early development, polycomb-regulated genes are connected through long-range 3D interactions which resolve upon differentiation. Here, we report that polycomb looping is controlled by H3K27me3 spreading and regulates target gene silencing and cell fate specification. Using glioma-derived H3 Lys-27-Met (H3K27M) mutations as tools to restrict H3K27me3 deposition, we show that H3K27me3 confinement concentrates the chromatin pool of cPRC1, resulting in heightened 3D interactions mirroring chromatin architecture of pluripotency, and stringent gene repression that maintains cells in progenitor states to facilitate tumor development. Conversely, H3K27me3 spread in pluripotent stem cells, following neural differentiation or loss of the H3K36 methyltransferase NSD1, dilutes cPRC1 concentration and dissolves polycomb loops. These results identify the regulatory principles and disease implications of polycomb looping and nominate histone modification-guided distribution of reader complexes as an important mechanism for nuclear compartment organization.

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.

Pediatric embryonal brain tumors with multilayered rosettes demonstrate a unique oncogenic amplification of the chromosome 19 miRNA cluster, C19MC. Because oncogenic lesions often cause deregulation of vascular effectors, including procoagulant tissue factor (TF), this study explores whether there is a link between C19MC oncogenic miRNAs (oncomirs) and the coagulant properties of cancer cells, a question previously not studied. In a pediatric embryonal brain tumor tissue microarray, we observed an association between C19MC amplification and reduced fibrin content and TF expression, indicative of reduced procoagulant activity. In medulloblastoma cell lines (DAOY and UW228) engineered to express miR-520g, a biologically active constituent of the C19MC cluster, we observed reduced TF expression, procoagulant and TF signaling activities (responses to factor VIIa stimulation), and diminished TF emission as cargo of extracellular vesicles. Antimir and luciferase reporter assays revealed a specific and direct effect of miR-520g on the TF 3' untranslated region. Although the endogenous MIR520G locus is methylated in differentiated cells, exposure of DAOY cells to 5-aza-2'-deoxycytidine or their growth as stem cell-like spheres up-regulated endogenous miR-520g with a coincident reduction in TF expression. We propose that the properties of tumors harboring oncomirs may include unique alterations of the vascular microenvironment, including deregulation of TF, with a possible impact on the biology, therapy, and hemostatic adverse effects of both disease progression and treatment.

Angiocentric gliomas are pediatric low-grade gliomas (PLGGs) without known recurrent genetic drivers. We performed genomic analysis of new and published data from 249 PLGGs, including 19 angiocentric gliomas. We identified MYB-QKI fusions as a specific and single candidate driver event in angiocentric gliomas. In vitro and in vivo functional studies show that MYB-QKI rearrangements promote tumorigenesis through three mechanisms: MYB activation by truncation, enhancer translocation driving aberrant MYB-QKI expression and hemizygous loss of the tumor suppressor QKI. To our knowledge, this represents the first example of a single driver rearrangement simultaneously transforming cells via three genetic and epigenetic mechanisms in a tumor.

In neuroblastoma, the most common solid tumor of childhood, excellent prognosis is associated with extensive Schwann cell (SC) content and high-level expression of the neurotrophin receptor, NTRK1/TrkA, which is known to mediate neuroblastoma cell differentiation. We hypothesized that both stromal composition and neuroblastic differentiation are based on bidirectional neuroblastoma-SC interaction. Reanalysis of microarray data from human SY5Y neuroblastoma cells stably transfected with either NTRK1 or NTRK2 revealed upregulation of the mRNA for the SC growth factor, NRG1, in NTRK1-positive cells. Media conditioned by NTRK1-expressing neuroblastoma cells induced SC proliferation and migration, while antibody-based NRG1 neutralization significantly decreased these effects. Vice versa, NRG1-stimulated SC secreted the NTRK1-specific ligand, NGF. SC-conditioned medium activated the NTRK1 receptor in a neuroblastoma cell culture model conditionally expressing NTRK1 and induced differentiation markers in NTRK1-expressing cells. NTRK1 induction in neuroblastoma xenografts mixed with primary SC also significantly reduced tumor growth in vivo. We propose a model for NTRK1-mediated and NRG1-dependent attraction of adjacent SC, which in turn induce neuroblastic differentiation by secretion of the NTRK1-specific ligand, NGF. These findings have implications for understanding the mature and less malignant neuroblastoma phenotype associated with NTRK1 expression, and could assist the development of new therapeutic strategies for neuroblastoma differentiation.

We recently reported that atypical teratoid rhabdoid tumors (ATRTs) comprise at least two transcriptional subtypes with different clinical outcomes; however, the mechanisms underlying therapeutic heterogeneity remained unclear. In this study, we analyzed 191 primary ATRTs and 10 ATRT cell lines to define the genomic and epigenomic landscape of ATRTs and identify subgroup-specific therapeutic targets. We found ATRTs segregated into three epigenetic subgroups with distinct genomic profiles, SMARCB1 genotypes, and chromatin landscape that correlated with differential cellular responses to a panel of signaling and epigenetic inhibitors. Significantly, we discovered that differential methylation of a PDGFRB-associated enhancer confers specific sensitivity of group 2 ATRT cells to dasatinib and nilotinib, and suggest that these are promising therapies for this highly lethal ATRT subtype.