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.

Intrinsic malignant brain tumors, such as glioblastomas are frequently resistant to immune checkpoint blockade (ICB) with few hypermutated glioblastomas showing response. Modeling patient-individual resistance is challenging due to the lack of predictive biomarkers and limited accessibility of tissue for serial biopsies. Here, we investigate resistance mechanisms to anti-PD-1 and anti-CTLA-4 therapy in syngeneic hypermutated experimental gliomas and show a clear dichotomy and acquired immune heterogeneity in ICB-responder and non-responder tumors. We made use of this dichotomy to establish a radiomic signature predicting tumor regression after pseudoprogression induced by ICB therapy based on serial magnetic resonance imaging. We provide evidence that macrophage-driven ICB resistance is established by CD4 T cell suppression and Treg expansion in the tumor microenvironment via the PD-L1/PD-1/CD80 axis. These findings uncover an unexpected heterogeneity of response to ICB in strictly syngeneic tumors and provide a rationale for targeting PD-L1-expressing tumor-associated macrophages to overcome resistance to ICB.

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.

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.

Concurrent malignant brain tumors in patients with multiple sclerosis (MS) constitute a rare but paradigmatic phenomenon for studying neuroimmunological mechanisms from both molecular and clinical perspectives.

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.

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.

Medulloblastomas are the most common malignant brain tumours in children. Identifying and understanding the genetic events that drive these tumours is critical for the development of more effective diagnostic, prognostic and therapeutic strategies. Recently, our group and others described distinct molecular subtypes of medulloblastoma on the basis of transcriptional and copy number profiles. Here we use whole-exome hybrid capture and deep sequencing to identify somatic mutations across the coding regions of 92 primary medulloblastoma/normal pairs. Overall, medulloblastomas have low mutation rates consistent with other paediatric tumours, with a median of 0.35 non-silent mutations per megabase. We identified twelve genes mutated at statistically significant frequencies, including previously known mutated genes in medulloblastoma such as CTNNB1, PTCH1, MLL2, SMARCA4 and TP53. Recurrent somatic mutations were newly identified in an RNA helicase gene, DDX3X, often concurrent with CTNNB1 mutations, and in the nuclear co-repressor (N-CoR) complex genes GPS2, BCOR and LDB1. We show that mutant DDX3X potentiates transactivation of a TCF promoter and enhances cell viability in combination with mutant, but not wild-type, β-catenin. Together, our study reveals the alteration of WNT, hedgehog, histone methyltransferase and now N-CoR pathways across medulloblastomas and within specific subtypes of this disease, and nominates the RNA helicase DDX3X as a component of pathogenic β-catenin signalling in medulloblastoma.

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.

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.

Pediatric patients with relapsed or refractory disease represent a population with a desperate medical need. The aim of the INFORM (INdividualized Therapy FOr Relapsed Malignancies in Childhood) program is to translate next generation molecular diagnostics into a biomarker driven treatment strategy. The program consists of two major foundations: the INFORM registry providing a molecular screening platform and the INFORM2 series of biomarker driven phase I/II trials. The INFORM2 NivEnt trial aims to determine the recommended phase 2 dose (RP2D) of the combination treatment of nivolumab and entinostat (phase I) and to evaluate activity and safety (phase II).

Diffuse IDH-mutant astrocytoma mostly occurs in adults and carries a favorable prognosis compared to IDH-wildtype malignant gliomas. Acquired mismatch repair deficiency is known to occur in recurrent IDH-mutant gliomas as resistance mechanism towards alkylating chemotherapy. In this multi-institutional study, we report a novel epigenetic group of 32 IDH-mutant gliomas with proven or suspected hereditary mismatch repair deficiency. None of the tumors exhibited a combined 1p/19q deletion. These primary mismatch repair-deficient IDH-mutant astrocytomas (PMMRDIA) were histologically high-grade and were mainly found in children, adolescents and young adults (median age 14 years). Mismatch repair deficiency syndromes (Lynch or Constitutional Mismatch Repair Deficiency Syndrom (CMMRD)) were clinically diagnosed and/or germline mutations in DNA mismatch repair genes (MLH1, MSH6, MSH2) were found in all cases, except one case with a family and personal history of colon cancer and another case with MSH6-deficiency available only as recurrent tumor. Loss of at least one of the mismatch repair proteins was detected via immunohistochemistry in all, but one case analyzed. Tumors displayed a hypermutant genotype and microsatellite instability was present in more than half of the sequenced cases. Integrated somatic mutational and chromosomal copy number analyses showed frequent inactivation of TP53, RB1 and activation of RTK/PI3K/AKT pathways. In contrast to the majority of IDH-mutant gliomas, more than 60% of the samples in our cohort presented with an unmethylated MGMT promoter. While the rate of immuno-histochemical ATRX loss was reduced, variants of unknown significance were more frequently detected possibly indicating a higher frequency of ATRX inactivation by protein malfunction. Compared to reference cohorts of other IDH-mutant gliomas, primary mismatch repair-deficient IDH-mutant astrocytomas have by far the worst clinical outcome with a median survival of only 15 months irrespective of histological or molecular features. The findings reveal a so far unknown entity of IDH-mutant astrocytoma with high prognostic relevance. Diagnosis can be established by aligning with the characteristic DNA methylation profile, by DNA-sequencing-based proof of mismatch repair deficiency or immunohistochemically demonstrating loss-of-mismatch repair proteins.

Medulloblastomas (MB) molecularly designated as Group 3 (Grp 3) MB represent a more clinically aggressive tumor variant which, as a group, displays heterogeneous molecular characteristics and disease outcomes. Reliable risk stratification of Grp 3 MB would allow for appropriate assignment of patients to aggressive treatment protocols and, vice versa, for sparing adverse effects of high-dose radio-chemotherapy in patients with standard or low-risk tumors. Here we performed RNA-based analysis on an international cohort of 179 molecularly designated Grp 3 MB treated with HIT protocols. We analyzed the clinical significance of differentially expressed genes, thereby developing optimal prognostic subdivision of this MB molecular group. We compared the transcriptome profiles of two Grp 3 MB subsets with various outcomes (76 died within the first 60 months vs. 103 survived this period) and identified 224 differentially expressed genes (DEG) between these two clinical groups (Limma R algorithm, adjusted p-value < 0.05). We selected the top six DEG overexpressed in the unfavorable cohort for further survival analysis and found that expression of all six genes strongly correlated with poor outcomes. However, only high expression of KIRREL2 was identified as an independent molecular prognostic indicator of poor patients' survival. Based on clinical and molecular patterns, four risk categories were outlined for Grp 3 MB patients: i. low-risk: M0-1/MYC non-amplified/KIRREL2 low (n = 48; 5-year OS-95%); ii. standard-risk: M0-1/MYC non-amplified/KIRREL2 high or M2-3/MYC non-amplified/KIRREL2 low (n = 65; 5-year OS-70%); iii. high-risk: M2-3/MYC non-amplified/KIRREL2 high (n = 36; 5-year OS-30%); iv. very high risk-all MYC amplified tumors (n = 30; 5-year OS-0%). Cross-validated survival models incorporating KIRREL2 expression with clinical features allowed for the reclassification of up to 50% of Grp 3 MB patients into a more appropriate risk category. Finally, KIRREL2 immunopositivity was also identified as a predictive indicator of Grp 3 MB poor survival, thus suggesting its application as a possible prognostic marker in routine clinical settings. Our results indicate that integration of KIRREL2 expression in risk stratification models may improve Grp 3 MB outcome prediction. Therefore, simple gene and/or protein expression analyses for this molecular marker could be easily adopted for Grp 3 MB prognostication and may help in assigning patients to optimal therapeutic approaches in prospective clinical trials.

Chemotherapy resistance is a persistent clinical problem in relapsed high-risk neuroblastomas. We tested a panel of 15 drugs for sensitization of neuroblastoma cells to the conventional chemotherapeutic vincristine, identifying tariquidar, an inhibitor of the transmembrane pump P-glycoprotein (P-gp/ABCB1), and the ERBB family inhibitor afatinib as the top resistance breakers. Both compounds were efficient in sensitizing neuroblastoma cells to vincristine in trypan blue exclusion assays and in inducing apoptotic cell death. The evaluation of ERBB signaling revealed no functional inhibition, that is, dephosphorylation of the downstream pathways upon afatinib treatment but direct off-target interference with P-gp function. Depletion of ABCB1, but not ERRB4, sensitized cells to vincristine treatment. P-gp inhibition substantially broke vincristine resistance in vitro and in vivo (zebrafish embryo xenograft). The analysis of gene expression datasets of more than 50 different neuroblastoma cell lines (primary and relapsed) and more than 160 neuroblastoma patient samples from the pediatric precision medicine platform INFORM (Individualized Therapy For Relapsed Malignancies in Childhood) confirmed a pivotal role of P-gp specifically in neuroblastoma resistance at relapse, while the ERBB family appears to play a minor part.