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.

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.

Pilocytic astrocytoma (PA), the most common pediatric brain tumor, is driven by aberrant mitogen-activated protein kinase signaling most commonly caused by BRAF gene fusions or activating mutations. While 5-year overall survival rates exceed 95%, tumor recurrence or progression constitutes a major clinical challenge in incompletely resected tumors. Here, we used similarity network fusion (SNF) analysis in an integrative multi-omics approach employing RNA transcriptomic and mass spectrometry-based proteomic profiling to molecularly characterize PA tissue samples from 62 patients. Thereby, we uncovered that PAs segregated into two molecularly distinct groups, namely, Group 1 and Group 2, which were validated in three non-overlapping cohorts. Patients with Group 1 tumors were significantly younger and showed worse progression-free survival compared to patients with group 2 tumors. Ingenuity pathways analysis (IPA) and gene set enrichment analysis (GSEA) revealed that Group 1 tumors were enriched for immune response pathways, such as interferon signaling, while Group 2 tumors showed enrichment for action potential and neurotransmitter signaling pathways. Analysis of immune cell-related gene signatures showed an enrichment of infiltrating T Cells in Group 1 versus Group 2 tumors. Taken together, integrative multi-omics of PA identified biologically distinct and prognostically relevant tumor groups that may improve risk stratification of this single pathway driven tumor type.

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.

Ten-eleven translocation (TET) proteins are crucial epigenetic regulators highly conserved in multicellular organisms. TETs' enzymatic function in demethylating 5-methyl cytosine in DNA is required for proper development and TETs are frequently mutated in cancer. Recently, Drosophila melanogaster Tet (dTet) was shown to be highly expressed in developing fly brains and discovered to play an important role in brain and muscle development as well as fly behavior. Furthermore, dTet was shown to have different substrate specificity compared with mammals. However, the exact role dTet plays in glial cells and how ectopic TET expression in glial cells contributes to tumorigenesis and glioma is still not clear. Here, we report a novel role for dTet specifically in glial cell organization and number. We show that loss of dTet affects the organization of a specific glia population in the optic lobe, the "optic chiasm" glia. Additionally, we find irregularities in axon patterns in the ventral nerve cord (VNC) both, in the midline and longitudinal axons. These morphologic glia and axonal defects were accompanied by locomotor defects in developing larvae escalating to immobility in adult flies. Furthermore, glia homeostasis was disturbed in dTet-deficient brains manifesting in gain of glial cell numbers and increased proliferation. Finally, we establish a Drosophila model to understand the impact of human TET3 in glia and find that ectopic expression of hTET3 in dTet-expressing cells causes glia expansion in larval brains and affects sleep/rest behavior and the circadian clock in adult flies.

Ependymomas (EPN) are central nervous system tumors comprising both aggressive and more benign molecular subtypes. However, therapy of the high-risk subtypes posterior fossa group A (PF-A) and supratentorial RELA-fusion positive (ST-RELA) is limited to gross total resection and radiotherapy, as effective systemic treatment concepts are still lacking. We have recently described fibroblast growth factor receptors 1 and 3 (FGFR1/FGFR3) as oncogenic drivers of EPN. However, the underlying molecular mechanisms and their potential as therapeutic targets have not yet been investigated in detail. Making use of transcriptomic data across 467 EPN tissues, we found that FGFR1 and FGFR3 were both widely expressed across all molecular groups. FGFR3 mRNA levels were enriched in ST-RELA showing the highest expression among EPN as well as other brain tumors. We further identified high expression levels of fibroblast growth factor 1 and 2 (FGF1, FGF2) across all EPN subtypes while FGF9 was elevated in ST-EPN. Interrogation of our EPN single-cell RNA-sequencing data revealed that FGFR3 was further enriched in cycling and progenitor-like cell populations. Corroboratively, we found FGFR3 to be predominantly expressed in radial glia cells in both mouse embryonal and human brain datasets. Moreover, we detected alternative splicing of the FGFR1/3-IIIc variant, which is known to enhance ligand affinity and FGFR signaling. Dominant-negative interruption of FGFR1/3 activation in PF-A and ST-RELA cell models demonstrated inhibition of key oncogenic pathways leading to reduced cell growth and stem cell characteristics. To explore the feasibility of therapeutically targeting FGFR, we tested a panel of FGFR inhibitors in 12 patient-derived EPN cell models revealing sensitivity in the low-micromolar to nano-molar range. Finally, we gain the first clinical evidence for the activity of the FGFR inhibitor nintedanib in the treatment of a patient with recurrent ST-RELA. Together, these preclinical and clinical data suggest FGFR inhibition as a novel and feasible approach to combat aggressive EPN.

This study aimed to re-evaluate the prognostic impact of TP53 mutations and to identify specific chromosomal aberrations as possible prognostic markers in WNT-activated medulloblastoma (WNT-MB). In a cohort of 191 patients with WNT-MBs, mutations in CTNNB1, APC, and TP53 were analyzed by DNA sequencing. Chromosomal copy-number aberrations were assessed by molecular inversion probe technology (MIP), SNP6, or 850k methylation array hybridization. Prognostic impact was evaluated in 120 patients with follow-up data from the HIT2000 medulloblastoma trial or HIT registries. CTNNB1 mutations were present in 92.2%, and APC mutations in 6.8% of samples. One CTNNB1 wild-type tumor gained WNT activation due to homozygous FBXW7 deletion. Monosomy 6 was present in 78.6%, and more frequent in children than adults. 16.1% of tumor samples showed TP53 mutations, of those 60% with nuclear positivity for the p53 protein. Loss of heterozygosity at the TP53 locus (chromosome 17p13.1) was found in 40.7% (11/27) of TP53 mutant tumor samples and in 12.6% of TP53 wild-type cases (13/103). Patients with tumors harboring TP53 mutations showed significant worse progression-free survival (PFS; 5-year-PFS 68% versus 93%, p = 0.001), and were enriched for chromosomes 17p (p = 0.001), 10, and 13 losses. Gains of OTX2 (14q22.3) occurred in 38.9% of samples and were associated with poor PFS and OS (5-year-PFS 72% versus 93%, p = 0.017 resp. 5-year-OS 83% versus 97%, p = 0.006). Multivariable Cox regression analysis for PFS/OS identified both genetic alterations as independent prognostic markers. Our data suggest that patients with WNT-MB carrying TP53 mutations or OTX2 gains (58.1%) are at higher risk of relapse. Eligibility of these patients for therapy de-escalation trials needs to be debated.

Current therapies for medulloblastoma, a highly malignant childhood brain tumour, impose debilitating effects on the developing child, and highlight the need for molecularly targeted treatments with reduced toxicity. Previous studies have been unable to identify the full spectrum of driver genes and molecular processes that operate in medulloblastoma subgroups. Here we analyse the somatic landscape across 491 sequenced medulloblastoma samples and the molecular heterogeneity among 1,256 epigenetically analysed cases, and identify subgroup-specific driver alterations that include previously undiscovered actionable targets. Driver mutations were confidently assigned to most patients belonging to Group 3 and Group 4 medulloblastoma subgroups, greatly enhancing previous knowledge. New molecular subtypes were differentially enriched for specific driver events, including hotspot in-frame insertions that target KBTBD4 and 'enhancer hijacking' events that activate PRDM6. Thus, the application of integrative genomics to an extensive cohort of clinical samples derived from a single childhood cancer entity revealed a series of cancer genes and biologically relevant subtype diversity that represent attractive therapeutic targets for the treatment of patients with medulloblastoma.

Medulloblastoma (MB) is the most common malignant primary pediatric brain cancer. Among the most aggressive subtypes, Group 3 and Group 4 originate from stem/progenitor cells, frequently metastasize, and often display the worst prognosis, yet we know the least about the molecular mechanisms driving their progression. Here, we show that the transcription factor orthodenticle homeobox 2 (OTX2) promotes self-renewal while inhibiting differentiation in vitro and increases tumor initiation from MB stem/progenitor cells in vivo. To determine how OTX2 contributes to these processes, we employed complementary bioinformatic approaches to characterize the OTX2 regulatory network and identified novel relationships between OTX2 and genes associated with neuronal differentiation and axon guidance signaling in Group 3 and Group 4 MB stem/progenitor cells. In particular, OTX2 levels were negatively correlated with semaphorin (SEMA) signaling, as expression of 9 SEMA pathway genes is upregulated following OTX2 knockdown with some being potential direct OTX2 targets. Importantly, this negative correlation was also observed in patient samples, with lower expression of SEMA4D associated with poor outcome specifically in Group 4 tumors. Functional proof-of-principle studies demonstrated that increased levels of select SEMA pathway genes are associated with decreased self-renewal and growth in vitro and in vivo and that RHO signaling, known to mediate the effects of SEMA genes, is contributing to the OTX2 KD phenotype. Our study provides mechanistic insight into the networks controlled by OTX2 in MB stem/progenitor cells and reveals novel roles for axon guidance genes and their downstream effectors as putative tumor suppressors in MB.

Medulloblastoma is associated with rare hereditary cancer predisposition syndromes; however, consensus medulloblastoma predisposition genes have not been defined and screening guidelines for genetic counselling and testing for paediatric patients are not available. We aimed to assess and define these genes to provide evidence for future screening guidelines.

Giant cell lesions of the jaw (GCLJ) are debilitating tumors of unknown origin with limited available therapies. Here, we analyze 58 sporadic samples using next generation or targeted sequencing and report somatic, heterozygous, gain-of-function mutations in KRAS, FGFR1, and p.M713V/I-TRPV4 in 72% (42/58) of GCLJ. TRPV4 p.M713V/I mutations are exclusive to central GCLJ and occur at a critical position adjacent to the cation permeable pore of the channel. Expression of TRPV4 mutants in HEK293 cells leads to increased cell death, as well as increased constitutive and stimulated channel activity, both of which can be prevented using TRPV4 antagonists. Furthermore, these mutations induce sustained activation of ERK1/2, indicating that their effects converge with that of KRAS and FGFR1 mutations on the activation of the MAPK pathway in GCLJ. Our data extend the spectrum of TRPV4 channelopathies and provide rationale for the use of TRPV4 and RAS/MAPK antagonists at the bedside in GCLJ.

Histone H3.3 glycine 34 to arginine/valine (G34R/V) mutations drive deadly gliomas and show exquisite regional and temporal specificity, suggesting a developmental context permissive to their effects. Here we show that 50% of G34R/V tumors (n = 95) bear activating PDGFRA mutations that display strong selection pressure at recurrence. Although considered gliomas, G34R/V tumors actually arise in GSX2/DLX-expressing interneuron progenitors, where G34R/V mutations impair neuronal differentiation. The lineage of origin may facilitate PDGFRA co-option through a chromatin loop connecting PDGFRA to GSX2 regulatory elements, promoting PDGFRA overexpression and mutation. At the single-cell level, G34R/V tumors harbor dual neuronal/astroglial identity and lack oligodendroglial programs, actively repressed by GSX2/DLX-mediated cell fate specification. G34R/V may become dispensable for tumor maintenance, whereas mutant-PDGFRA is potently oncogenic. Collectively, our results open novel research avenues in deadly tumors. G34R/V gliomas are neuronal malignancies where interneuron progenitors are stalled in differentiation by G34R/V mutations and malignant gliogenesis is promoted by co-option of a potentially targetable pathway, PDGFRA signaling.

Ependymoma is a heterogeneous entity of central nervous system tumors with well-established molecular groups. Here, we apply single-cell RNA sequencing to analyze ependymomas across molecular groups and anatomic locations to investigate their intratumoral heterogeneity and developmental origins. Ependymomas are composed of a cellular hierarchy initiating from undifferentiated populations, which undergo impaired differentiation toward three lineages of neuronal-glial fate specification. While prognostically favorable groups of ependymoma predominantly harbor differentiated cells, aggressive groups are enriched for undifferentiated cell populations. The delineated transcriptomic signatures correlate with patient survival and define molecular dependencies for targeted treatment approaches. Taken together, our analyses reveal a developmental hierarchy underlying ependymomas relevant to biological and clinical behavior.

YAP1 is a transcriptional coactivator and the principal effector of the Hippo signaling pathway, which is causally implicated in human cancer. Several YAP1 gene fusions have been identified in various human cancers and identifying the essential components of this family of gene fusions has significant therapeutic value. Here, we show that the YAP1 gene fusions YAP1-MAMLD1, YAP1-FAM118B, YAP1-TFE3, and YAP1-SS18 are oncogenic in mice. Using reporter assays, RNA-seq, ChIP-seq, and loss-of-function mutations, we can show that all of these YAP1 fusion proteins exert TEAD-dependent YAP activity, while some also exert activity of the C'-terminal fusion partner. The YAP activity of the different YAP1 fusions is resistant to negative Hippo pathway regulation due to constitutive nuclear localization and resistance to degradation of the YAP1 fusion proteins. Genetic disruption of the TEAD-binding domain of these oncogenic YAP1 fusions is sufficient to inhibit tumor formation in vivo, while pharmacological inhibition of the YAP1-TEAD interaction inhibits the growth of YAP1 fusion-expressing cell lines in vitro. These results highlight TEAD-dependent YAP activity found in these gene fusions as critical for oncogenesis and implicate these YAP functions as potential therapeutic targets in YAP1 fusion-positive tumors.