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.

Malignant rhabdoid tumors (MRTs) represent one of the most aggressive childhood malignancies. No effective treatment options are available, and prognosis is, therefore, dismal. Previous studies have demonstrated that tumor organoids capture the heterogeneity of patient tumors and can be used to predict patient response to therapy. Here, we perform drug screening on patient-derived normal and tumor organoids to identify MRT-specific therapeutic vulnerabilities. We identify neddylation inhibitor MLN4924 as a potential therapeutic agent. Mechanistically, we find increased neddylation in MRT organoids and tissues and show that MLN4924 induces a cytotoxic response via upregulation of the unfolded protein response. Lastly, we demonstrate in vivo efficacy in an MRT PDX mouse model, in which single-agent MLN4924 treatment significantly extends survival. Our study demonstrates that organoids can be used to find drugs selectively targeting tumor cells while leaving healthy cells unharmed and proposes neddylation inhibition as a therapeutic strategy in MRT.

Medulloblastomas with chromothripsis developing in children with Li-Fraumeni Syndrome (germline TP53 mutations) are highly aggressive brain tumors with dismal prognosis. Conventional photon radiotherapy and DNA-damaging chemotherapy are not successful for these patients and raise the risk of secondary malignancies. We hypothesized that the pronounced homologous recombination deficiency in these tumors might offer vulnerabilities that can be therapeutically utilized in combination with high linear energy transfer carbon ion radiotherapy.

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.

Medulloblastoma, the most common malignant childhood brain tumor, exhibits distinct molecular subtypes and cellular origins. Genetic alterations driving medulloblastoma initiation and progression remain poorly understood. Herein, we identify GNAS, encoding the G protein Gαs, as a potent tumor suppressor gene that, when expressed at low levels, defines a subset of aggressive Sonic hedgehog (SHH)-driven human medulloblastomas. Ablation of the single Gnas gene in anatomically distinct progenitors in mice is sufficient to induce Shh-associated medulloblastomas, which recapitulate their human counterparts. Gαs is highly enriched at the primary cilium of granule neuron precursors and suppresses Shh signaling by regulating both the cAMP-dependent pathway and ciliary trafficking of Hedgehog pathway components. Elevation in levels of a Gαs effector, cAMP, effectively inhibits tumor cell proliferation and progression in Gnas-ablated mice. Thus, our gain- and loss-of-function studies identify a previously unrecognized tumor suppressor function for Gαs that can be found consistently across Shh-group medulloblastomas of disparate cellular and anatomical origins, highlighting G protein modulation as a potential therapeutic avenue.

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.

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.

Genomic gain of the proto-oncogene transcription factor gene MYCN is associated with poor prognosis in several childhood cancers. Here we present a comprehensive copy number analysis of MYCN in Wilms tumour (WT), demonstrating that gain of this gene is associated with anaplasia and with poorer relapse-free and overall survival, independent of histology. Using whole exome and gene-specific sequencing, together with methylation and expression profiling, we show that MYCN is targeted by other mechanisms, including a recurrent somatic mutation, P44L, and specific DNA hypomethylation events associated with MYCN overexpression in tumours with high risk histologies. We describe parallel evolution of genomic copy number gain and point mutation of MYCN in the contralateral tumours of a remarkable bilateral case in which independent contralateral mutations of TP53 also evolve over time. We report a second bilateral case in which MYCN gain is a germline aberration. Our results suggest a significant role for MYCN dysregulation in the molecular biology of Wilms tumour. We conclude that MYCN gain is prognostically significant, and suggest that the novel P44L somatic variant is likely to be an activating mutation.

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.

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.

Suppressor of Fused (SuFu), a tumour suppressor mutated in medulloblastoma, is a central player of Hh signalling, a pathway crucial for development and deregulated in cancer. Although the control of Gli transcription factors by SuFu is critical in Hh signalling, our understanding of the mechanism regulating this key event remains limited. Here, we show that the Itch/β-arrestin2 complex binds SuFu and induces its Lys63-linked polyubiquitylation without affecting its stability. This process increases the association of SuFu with Gli3, promoting the conversion of Gli3 into a repressor, which keeps Hh signalling off. Activation of Hh signalling antagonises the Itch-dependent polyubiquitylation of SuFu. Notably, different SuFu mutations occurring in medulloblastoma patients are insensitive to Itch activity, thus leading to deregulated Hh signalling and enhancing medulloblastoma cell growth. Our findings uncover mechanisms controlling the tumour suppressive functions of SuFu and reveal that their alterations are implicated in medulloblastoma tumorigenesis.

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.

Glioblastomas strongly invade the brain by infiltrating into the white matter along myelinated nerve fiber tracts even though the myelin protein Nogo-A prevents cell migration by activating inhibitory RhoA signaling. The mechanisms behind this long-known phenomenon remained elusive so far, precluding a targeted therapeutic intervention. This study demonstrates that the prevalent activation of AKT in gliomas increases the ER protein-folding capacity and enables tumor cells to utilize a side effect of RhoA activation: the perturbation of the IRE1α-mediated decay of SPARC mRNA. Once translation is initiated, glioblastoma cells rapidly secrete SPARC to block Nogo-A from inhibiting migration via RhoA. By advanced ultramicroscopy for studying single-cell invasion in whole, undissected mouse brains, we show that gliomas require SPARC for invading into white matter structures. SPARC depletion reduces tumor dissemination that significantly prolongs survival and improves response to cytostatic therapy. Our finding of a novel RhoA-IRE1 axis provides a druggable target for interfering with SPARC production and underscores its therapeutic value.

Persistent mortality rates of medulloblastoma (MB) and severe side effects of the current therapies require the definition of the molecular mechanisms that contribute to tumor progression. Using cultured MB cancer stem cells and xenograft tumors generated in mice, we show that low expression of miR-326 and its host gene β-arrestin1 (ARRB1) promotes tumor growth enhancing the E2F1 pro-survival function. Our models revealed that miR-326 and ARRB1 are controlled by a bivalent domain, since the H3K27me3 repressive mark is found at their regulatory region together with the activation-associated H3K4me3 mark. High levels of EZH2, a feature of MB, are responsible for the presence of H3K27me3. Ectopic expression of miR-326 and ARRB1 provides hints into how their low levels regulate E2F1 activity. MiR-326 targets E2F1 mRNA, thereby reducing its protein levels; ARRB1, triggering E2F1 acetylation, reverses its function into pro-apoptotic activity. Similar to miR-326 and ARRB1 overexpression, we also show that EZH2 inhibition restores miR-326/ARRB1 expression, limiting E2F1 pro-proliferative activity. Our results reveal a new regulatory molecular axis critical for MB progression.

Medulloblastoma is the most common malignant brain tumor in children. Here we describe a medulloblastoma model using Induced pluripotent stem (iPS) cell-derived human neuroepithelial stem (NES) cells generated from a Gorlin syndrome patient carrying a germline mutation in the sonic hedgehog (SHH) receptor PTCH1. We found that Gorlin NES cells formed tumors in mouse cerebellum mimicking human medulloblastoma. Retransplantation of tumor-isolated NES (tNES) cells resulted in accelerated tumor formation, cells with reduced growth factor dependency, enhanced neurosphere formation in vitro, and increased sensitivity to Vismodegib. Using our model, we identified LGALS1 to be a GLI target gene that is up-regulated in both Gorlin tNES cells and SHH-subgroup of medulloblastoma patients. Taken together, we demonstrate that NES cells derived from Gorlin patients can be used as a resource to model medulloblastoma initiation and progression and to identify putative targets.

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.

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.

We and others have demonstrated that MYC-amplified medulloblastoma (MB) cells are susceptible to class I histone deacetylase inhibitor (HDACi) treatment. However, single drug treatment with HDACi has shown limited clinical efficacy. We hypothesized that addition of a second compound acting synergistically with HDACi may enhance efficacy.