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.

In vivo functional investigation of oncogenes using somatic gene transfer has been successfully exploited to validate their role in tumorigenesis. For tumour suppressor genes this has proven more challenging due to technical aspects. To provide a flexible and effective method for investigating somatic loss-of-function alterations and their influence on tumorigenesis, we have established CRISPR/Cas9-mediated somatic gene disruption, allowing for in vivo targeting of TSGs. Here we demonstrate the utility of this approach by deleting single (Ptch1) or multiple genes (Trp53, Pten, Nf1) in the mouse brain, resulting in the development of medulloblastoma and glioblastoma, respectively. Using whole-genome sequencing (WGS) we characterized the medulloblastoma-driving Ptch1 deletions in detail and show that no off-targets were detected in these tumours. This method provides a fast and convenient system for validating the emerging wealth of novel candidate tumour suppressor genes and the generation of faithful animal models of human cancer.

The microscopic analysis of telomere features provides a wealth of information on the mechanism by which tumor cells maintain their unlimited proliferative potential. Accordingly, the analysis of telomeres in tissue sections of patient tumor samples can be exploited to obtain diagnostic information and to define tumor subgroups. In many instances, however, analysis of the image data is conducted by manual inspection of 2D images at relatively low resolution for only a small part of the sample. As the telomere feature signal distribution is frequently heterogeneous, this approach is prone to a biased selection of the information present in the image and lacks subcellular details. Here we address these issues by using an automated high-resolution imaging and analysis workflow that quantifies individual telomere features on tissue sections for a large number of cells. The approach is particularly suited to assess telomere heterogeneity and low abundant cellular subpopulations with distinct telomere characteristics in a reproducible manner. It comprises the integration of multi-color fluorescence in situ hybridization, immunofluorescence and DNA staining with targeted automated 3D fluorescence microscopy and image analysis. We apply our method to telomeres in glioblastoma and prostate cancer samples, and describe how the imaging data can be used to derive statistically reliable information on telomere length distribution or colocalization with PML nuclear bodies. We anticipate that relating this approach to clinical outcome data will prove to be valuable for pretherapeutic patient stratification.

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.

Dysembryoplastic neuroepithelial tumors (DNT) share V600E mutation in the BRAF gene with other low grade neuroepithelial tumors (LGNTs). FGFR1 internal tandem duplication of the tyrosine-kinase domain (FGFR1-ITD), another genetic alteration that also leads to MAP kinase pathway alteration, has been previously reported in LGNTs by whole-genome sequencing. In the present study we searched for FGFR1-ITD by droplet digital PCR (DDPCR™) and for FGFR1 point mutations by HRM-sequencing in a series of formalin-fixed paraffin-embedded (FFPE) LGNTs including 12 DNT, 2 oligodendrogliomas lacking IDH mutation and 1p/19q co- deletion (pediatric-type oligodendrogliomas; PTOs), 3 pediatric diffuse astrocytomas (PDAs), 14 gangliogliomas (GGs) and 5 pilocytic astrocytomas (PAs). We showed by DDPCR™ that 5/12 DNT, but none of the other LGNTs, demonstrated FGFR1-ITD. In addition, these cases also accumulated phosphorylated-FGFR1 protein as shown by immunohistochemistry. FGFR1G539R point mutation was only recorded in one DNT that also showed FGFR1-ITD. Interestingly, these FGFR1 alterations were mutually exclusive from BRAFV600E mutation that was recorded in 13 LGNTs (3 DNTs, 1 PTO, 2 PDAs, 5 GGs and 2 PAs). Therefore, FGFR1 alteration mainly represented by FGFR1-ITD is a frequent event in DNT. DDPCR™ is an easy and alternative method than whole-genome sequencing to detect FGFR1-ITD in FFPE brain tumors, in routine practice.

Epileptogenic tumors affecting children and young adults are a morphologically diverse collection of neuroepithelial neoplasms that, as a group, exhibit varying levels of glial and/or neuronal differentiation. Recent advances in molecular profiling technology, including comprehensive DNA sequencing and methylation analysis, have enabled the application of more precise and biologically relevant classification schemes to these tumors. In this report, we describe a morphologically and molecularly distinct epileptogenic neoplasm, the polymorphous low-grade neuroepithelial tumor of the young (PLNTY), which likely accounts for a sizable portion of oligodendroglioma-like tumors affecting the pediatric population. Characteristic microscopic findings most notably include infiltrative growth, the invariable presence of oligodendroglioma-like cellular components, and intense immunolabeling for cluster of differentiation 34 (CD34). Moreover, integrative molecular profiling reveals a distinct DNA methylation signature for PLNTYs, along with frequent genetic abnormalities involving either B-Raf proto-oncogene (BRAF) or fibroblast growth factor receptors 2 and 3 (FGFR2, FGFR3). These findings suggest that PLNTY represents a distinct biological entity within the larger spectrum of pediatric, low-grade neuroepithelial tumors.

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.

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.

Eukaryotic translation initiation factor 4E binding protein 1 (EIF4EBP1) encodes the 4EBP1 protein, a negative regulator of mRNA translation and a substrate of the mechanistic target of rapamycin (mTOR), whose function and relevance in cancer is still under debate. Here, we analyzed EIF4EBP1 expression in different glioma patient cohorts and investigated its mode of transcriptional regulation in glioblastoma cells. We verified that EIF4EBP1 mRNA is overexpressed in malignant gliomas, including isocitrate dehydrogenase (IDH)-wildtype glioblastomas, relative to non-neoplastic brain tissue in multiple publically available datasets. Our analyses revealed that EIF4EBP1 overexpression in malignant gliomas is neither due to gene amplification nor to altered DNA methylation, but rather results from aberrant transcriptional activation by distinct transcription factors. We found seven transcription factor candidates co-expressed with EIF4EBP1 in gliomas and bound to the EIF4EBP1 promoter, as revealed by chromatin immunoprecipitation (ChIP)-sequencing data. We investigated the ability of these candidates to activate the EIF4EBP1 promoter using luciferase reporter assays, which supported four transcription factors as candidate EIF4EBP1 regulators, namely MYBL2, ETS1, HIF-1A, and E2F6. Finally, by employing transient knock-down experiments to repress either of these transcription factors, we identified MYBL2 and ETS1 as the relevant transcriptional drivers of enhanced EIF4EBP1 expression in malignant glioma cells. Taken together, our findings confirm enhanced expression of EIF4EBP1 in malignant gliomas relative to non-neoplastic brain tissue and characterize the underlying molecular pathomechanisms.

N-Myc is a transcription factor that is aberrantly expressed in many tumor types and is often correlated with poor patient prognosis. Recently, several lines of evidence pointed to the fact that oncogenic activation of Myc family proteins is concomitant with reprogramming of tumor cells to cope with an enhanced need for metabolites during cell growth. These adaptions are driven by the ability of Myc proteins to act as transcriptional amplifiers in a tissue-of-origin specific manner. Here, we describe the effects of N-Myc overexpression on metabolic reprogramming in neuroblastoma cells. Ectopic expression of N-Myc induced a glycolytic switch that was concomitant with enhanced sensitivity towards 2-deoxyglucose, an inhibitor of glycolysis. Moreover, global metabolic profiling revealed extensive alterations in the cellular metabolome resulting from overexpression of N-Myc. Limited supply with either of the two main carbon sources, glucose or glutamine, resulted in distinct shifts in steady-state metabolite levels and significant changes in glutathione metabolism. Interestingly, interference with glutamine-glutamate conversion preferentially blocked proliferation of N-Myc overexpressing cells, when glutamine levels were reduced. Thus, our study uncovered N-Myc induction and nutrient levels as important metabolic master switches in neuroblastoma cells and identified critical nodes that restrict tumor cell proliferation.

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.

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.

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.

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.

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.

Primary brain tumors are a leading cause of cancer-related death in children, and medulloblastoma is the most common malignant pediatric brain tumor. The current molecular characterization of medulloblastoma is mainly based on protein-coding genes, while little is known about the involvement of long non-coding RNAs (lncRNAs). This study aimed to elucidate the role of the lncRNA OTX2-AS1 in medulloblastoma.

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.

The introduction of the classification of brain tumours based on their DNA methylation profile has significantly changed the diagnostic approach for cases with ambiguous histology, non-informative or contradictory molecular profiles or for entities where methylation profiling provides useful information for patient risk stratification, for example in medulloblastoma and ependymoma. We present our experience that combines a conventional molecular diagnostic approach with the complementary use of a DNA methylation-based classification tool, for adult brain tumours originating from local as well as national referrals. We report the frequency of IDH mutations in a large cohort of nearly 1550 patients, EGFR amplifications in almost 1900 IDH-wildtype glioblastomas, and histone mutations in 70 adult gliomas. We demonstrate how additional methylation-based classification has changed and improved our diagnostic approach. Of the 325 cases referred for methylome testing, 179 (56%) had a calibrated score of 0.84 and higher and were included in the evaluation. In these 179 samples, the diagnosis was changed in 45 (25%), refined in 86 (48%) and confirmed in 44 cases (25%). In addition, the methylation arrays contain copy number information that usefully complements the methylation profile. For example, EGFR amplification which is 95% concordant with our Real-Time PCR-based copy number assays. We propose here a diagnostic algorithm that integrates histology, conventional molecular tests and methylation arrays.

Biallelic inactivation of SMARCB1, encoding a member of the SWI/SNF chromatin remodeling complex, is the hallmark genetic aberration of atypical teratoid rhabdoid tumors (ATRT). Here, we report how loss of SMARCB1 affects the epigenome in these tumors. Using chromatin immunoprecipitation sequencing (ChIP-seq) on primary tumors for a series of active and repressive histone marks, we identified the chromatin states differentially represented in ATRTs compared with other brain tumors and non-neoplastic brain. Re-expression of SMARCB1 in ATRT cell lines enabled confirmation of our genome-wide findings for the chromatin states. Additional generation of ChIP-seq data for SWI/SNF and Polycomb group proteins and the transcriptional repressor protein REST determined differential dependencies of SWI/SNF and Polycomb complexes in regulation of diverse gene sets in ATRTs.

Resistance to radiation treatment remains a major clinical problem for patients with brain cancer. Medulloblastoma is the most common malignant brain tumor of childhood, and occurs in the cerebellum. Though radiation treatment has been critical in increasing survival rates in recent decades, the presence of resistant cells in a substantial number of medulloblastoma patients leads to relapse and death.