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.

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.

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.

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.

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.

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.

Cerebral cavernous malformations are clusters of aberrant vessels that can lead to severe neurological complications. Pathogenic loss-of-function variants in the CCM1, CCM2, or CCM3 gene are associated with the autosomal dominant form of the disease. While interpretation of variants in protein-coding regions of the genes is relatively straightforward, functional analyses are often required to evaluate the impact of non-coding variants. Because of multiple alternatively spliced transcripts and different transcription start points, interpretation of variants in the 5' untranslated and upstream regions of CCM1 is particularly challenging. Here, we identified a novel deletion of the non-coding exon 1 of CCM1 in a proband with multiple CCMs which was initially classified as a variant of unknown clinical significance. Using CRISPR/Cas9 genome editing in human iPSCs, we show that the deletion leads to loss of CCM1 protein and deregulation of KLF2, THBS1, NOS3, and HEY2 expression in iPSC-derived endothelial cells. Based on these results, the variant could be reclassified as likely pathogenic. Taken together, variants in regulatory regions need to be considered in genetic CCM analyses. Our study also demonstrates that modeling variants of unknown clinical significance in an iPSC-based system can help to come to a final diagnosis.

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.

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.

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.

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.

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.

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.

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.