Musashi1 (Msi1) is a highly conserved RNA-binding protein that is required during the development of the nervous system. Msi1 has been characterized as a stem cell marker, controlling the balance between self-renewal and differentiation, and has also been implicated in tumorigenesis, being highly expressed in multiple tumor types. We analyzed Msi1 expression in a large cohort of medulloblastoma samples and found that Msi1 is highly expressed in tumor tissue compared with normal cerebellum. Notably, high Msi1 expression levels proved to be a sign of poor prognosis. Msi1 expression was determined to be particularly high in molecular subgroups 3 and 4 of medulloblastoma. We determined that Msi1 is required for tumorigenesis because inhibition of Msi1 expression by small-interfering RNAs reduced the growth of Daoy medulloblastoma cells in xenografts. To characterize the participation of Msi1 in medulloblastoma, we conducted different high-throughput analyses. Ribonucleoprotein immunoprecipitation followed by microarray analysis (RIP-chip) was used to identify mRNA species preferentially associated with Msi1 protein in Daoy cells. We also used cluster analysis to identify genes with similar or opposite expression patterns to Msi1 in our medulloblastoma cohort. A network study identified RAC1, CTGF, SDCBP, SRC, PRL, and SHC1 as major nodes of an Msi1-associated network. Our results suggest that Msi1 functions as a regulator of multiple processes in medulloblastoma formation and could become an important therapeutic target.

Glioblastoma (GBM) is a typically lethal type of brain tumor with a median survival of 15 months postdiagnosis. This negative prognosis prompted the exploration of alternative treatment options. In particular, the reliance of GBM on angiogenesis triggered the development of anti-VEGF (vascular endothelial growth factor) blocking antibodies such as bevacizumab. Although its application in human GBM only increased progression-free periods but did not improve overall survival, physicians and researchers still utilize this treatment option due to the lack of adequate alternatives. In an attempt to improve the efficacy of anti-VEGF treatment, we explored the role of the egfl7 gene in malignant glioma. We found that the encoded extracellular matrix protein epidermal growth factor-like protein 7 (EGFL7) was secreted by glioma blood vessels but not glioma cells themselves, while no major role could be assigned to the parasitic miRNAs miR-126/126*. EGFL7 expression promoted glioma growth in experimental glioma models in vivo and stimulated tumor vascularization. Mechanistically, this was mediated by an upregulation of integrin α5β1 on the cellular surface of endothelial cells, which enhanced fibronectin-induced angiogenic sprouting. Glioma blood vessels that formed in vivo were more mature as determined by pericyte and smooth muscle cell coverage. Furthermore, these vessels were less leaky as measured by magnetic resonance imaging of extravasating contrast agent. EGFL7-inhibition using a specific blocking antibody reduced the vascularization of experimental gliomas and increased the life span of treated animals, in particular in combination with anti-VEGF and the chemotherapeutic agent temozolomide. Data allow for the conclusion that this combinatorial regimen may serve as a novel treatment option for GBM.

Triple-negative breast cancer (TNBC) is an aggressive subtype characterized by extensive intratumoral heterogeneity. To investigate the underlying biology, we conducted single-cell RNA-sequencing (scRNA-seq) of >1500 cells from six primary TNBC. Here, we show that intercellular heterogeneity of gene expression programs within each tumor is variable and largely correlates with clonality of inferred genomic copy number changes, suggesting that genotype drives the gene expression phenotype of individual subpopulations. Clustering of gene expression profiles identified distinct subgroups of malignant cells shared by multiple tumors, including a single subpopulation associated with multiple signatures of treatment resistance and metastasis, and characterized functionally by activation of glycosphingolipid metabolism and associated innate immunity pathways. A novel signature defining this subpopulation predicts long-term outcomes for TNBC patients in a large cohort. Collectively, this analysis reveals the functional heterogeneity and its association with genomic evolution in TNBC, and uncovers unanticipated biological principles dictating poor outcomes in this disease.

Diffuse midline glioma (DMG), H3 K27M-mutant, is a new entity in the updated WHO classification grouping together diffuse intrinsic pontine gliomas and infiltrating glial neoplasms of the midline harboring the same canonical mutation at the Lysine 27 of the histones H3 tail.Two hundred and fifteen patients younger than 18 years old with centrally-reviewed pediatric high-grade gliomas (pHGG) were included in this study. Comprehensive transcriptomic (n = 140) and methylation (n = 80) profiling was performed depending on the material available, in order to assess the biological uniqueness of this new entity compared to other midline and hemispheric pHGG.Tumor classification based on gene expression (GE) data highlighted the similarity of K27M DMG independently of their location along the midline. T-distributed Stochastic Neighbor Embedding (tSNE) analysis of methylation profiling confirms the discrimination of DMG from other well defined supratentorial tumor subgroups. Patients with diffuse intrinsic pontine gliomas (DIPG) and thalamic DMG exhibited a similarly poor prognosis (11.1 and 10.8 months median overall survival, respectively). Interestingly, H3.1-K27M and H3.3-K27M primary tumor samples could be distinguished based both on their GE and DNA methylation profiles, suggesting that they might arise from a different precursor or from a different epigenetic reorganization.These differences in DNA methylation profiles were conserved in glioma stem-like cell culture models of DIPG which mimicked their corresponding primary tumor. ChIP-seq profiling of H3K27me3 in these models indicate that H3.3-K27M mutated DIPG stem cells exhibit higher levels of H3K27 trimethylation which are correlated with fewer genes expressed by RNAseq. When considering the global distribution of the H3K27me3 mark, we observed that intergenic regions were more trimethylated in the H3.3-K27M mutated cells compared to the H3.1-K27M mutated ones.H3 K27M-mutant DMG represent a homogenous group of neoplasms compared to other pediatric gliomas that could be further separated based on the type of histone H3 variant mutated and their respective epigenetic landscapes. As these characteristics drive different phenotypes, these findings may have important implication for the design of future trials in these specific types of neoplasms.

Accurate pathological diagnosis is crucial for optimal management of patients with cancer. For the approximately 100 known tumour types of the central nervous system, standardization of the diagnostic process has been shown to be particularly challenging-with substantial inter-observer variability in the histopathological diagnosis of many tumour types. Here we present a comprehensive approach for the DNA methylation-based classification of central nervous system tumours across all entities and age groups, and demonstrate its application in a routine diagnostic setting. We show that the availability of this method may have a substantial impact on diagnostic precision compared to standard methods, resulting in a change of diagnosis in up to 12% of prospective cases. For broader accessibility, we have designed a free online classifier tool, the use of which does not require any additional onsite data processing. Our results provide a blueprint for the generation of machine-learning-based tumour classifiers across other cancer entities, with the potential to fundamentally transform tumour pathology.

Diverse genetic, epigenetic, and developmental programs drive glioblastoma, an incurable and poorly understood tumor, but their precise characterization remains challenging. Here, we use an integrative approach spanning single-cell RNA-sequencing of 28 tumors, bulk genetic and expression analysis of 401 specimens from the The Cancer Genome Atlas (TCGA), functional approaches, and single-cell lineage tracing to derive a unified model of cellular states and genetic diversity in glioblastoma. We find that malignant cells in glioblastoma exist in four main cellular states that recapitulate distinct neural cell types, are influenced by the tumor microenvironment, and exhibit plasticity. The relative frequency of cells in each state varies between glioblastoma samples and is influenced by copy number amplifications of the CDK4, EGFR, and PDGFRA loci and by mutations in the NF1 locus, which each favor a defined state. Our work provides a blueprint for glioblastoma, integrating the malignant cell programs, their plasticity, and their modulation by genetic drivers.

Spinal ependymal tumors form a histologically and molecularly heterogeneous group of tumors with generally good prognosis. However, their treatment can be challenging if infiltration of the spinal cord or dissemination throughout the central nervous system (CNS) occurs and, in these cases, clinical outcome remains poor. Here, we describe a new and relatively rare subgroup of spinal ependymal tumors identified using DNA methylation profiling that is distinct from other molecular subgroups of ependymoma. Copy number variation plots derived from DNA methylation arrays showed MYCN amplification as a characteristic genetic alteration in all cases of our cohort (n = 13), which was subsequently validated using fluorescence in situ hybridization. The histological diagnosis was anaplastic ependymoma (WHO Grade III) in ten cases and classic ependymoma (WHO Grade II) in three cases. Histological re-evaluation in five primary tumors and seven relapses showed characteristic histological features of ependymoma, namely pseudorosettes, GFAP- and EMA positivity. Electron microscopy revealed cilia, complex intercellular junctions and intermediate filaments in a representative sample. Taking these findings into account, we suggest to designate this molecular subgroup spinal ependymoma with MYCN amplification, SP-EPN-MYCN. SP-EPN-MYCN tumors showed distinct growth patterns with intradural, extramedullary localization mostly within the thoracic and cervical spine, diffuse leptomeningeal spread throughout the whole CNS and infiltrative invasion of the spinal cord. Dissemination was observed in 100% of cases. Despite high-intensity treatment, SP-EPN-MYCN showed significantly worse median progression free survival (PFS) (17 months) and median overall survival (OS) (87 months) than all other previously described molecular spinal ependymoma subgroups. OS and PFS were similar to supratentorial ependymoma with RELA-fusion (ST-EPN-RELA) and posterior fossa ependymoma A (PF-EPN-A), further highlighting the aggressiveness of this distinct new subgroup. We, therefore, propose to establish SP-EPN-MYCN as a new molecular subgroup in ependymoma and advocate for testing newly diagnosed spinal ependymal tumors for MYCN amplification.

YAP1 fusion-positive supratentorial ependymomas predominantly occur in infants, but the molecular mechanisms of oncogenesis are unknown. Here we show YAP1-MAMLD1 fusions are sufficient to drive malignant transformation in mice, and the resulting tumors share histo-molecular characteristics of human ependymomas. Nuclear localization of YAP1-MAMLD1 protein is mediated by MAMLD1 and independent of YAP1-Ser127 phosphorylation. Chromatin immunoprecipitation-sequencing analyses of human YAP1-MAMLD1-positive ependymoma reveal enrichment of NFI and TEAD transcription factor binding site motifs in YAP1-bound regulatory elements, suggesting a role for these transcription factors in YAP1-MAMLD1-driven tumorigenesis. Mutation of the TEAD binding site in the YAP1 fusion or repression of NFI targets prevents tumor induction in mice. Together, these results demonstrate that the YAP1-MAMLD1 fusion functions as an oncogenic driver of ependymoma through recruitment of TEADs and NFIs, indicating a rationale for preclinical studies to block the interaction between YAP1 fusions and NFI and TEAD transcription factors.

Extra-cranial malignant rhabdoid tumors (MRTs) and cranial atypical teratoid RTs (ATRTs) are heterogeneous pediatric cancers driven primarily by SMARCB1 loss. To understand the genome-wide molecular relationships between MRTs and ATRTs, we analyze multi-omics data from 140 MRTs and 161 ATRTs. We detect similarities between the MYC subgroup of ATRTs (ATRT-MYC) and extra-cranial MRTs, including global DNA hypomethylation and overexpression of HOX genes and genes involved in mesenchymal development, distinguishing them from other ATRT subgroups that express neural-like features. We identify five DNA methylation subgroups associated with anatomical sites and SMARCB1 mutation patterns. Groups 1, 3, and 4 exhibit cytotoxic T cell infiltration and expression of immune checkpoint regulators, consistent with a potential role for immunotherapy in rhabdoid tumor patients.

The rare benign giant cell tumour of bone (GCTB) is defined by an almost unique mutation in the H3.3 family of histone genes H3-3A or H3-3B; however, the same mutation is occasionally found in primary malignant bone tumours which share many features with the benign variant. Moreover, lung metastases can occur despite the absence of malignant histological features in either the primary or metastatic lesions. Herein we investigated the genetic events of 17 GCTBs including benign and malignant variants and the methylation profiles of 122 bone tumour samples including GCTBs. Benign GCTBs possessed few somatic alterations and no other known drivers besides the H3.3 mutation, whereas all malignant tumours harboured at least one additional driver mutation and exhibited genomic features resembling osteosarcomas, including high mutational burden, additional driver event(s), and a high degree of aneuploidy. The H3.3 mutation was found to predate the development of aneuploidy. In contrast to osteosarcomas, malignant H3.3-mutated tumours were enriched for a variety of alterations involving TERT, other than amplification, suggesting telomere dysfunction in the transformation of benign to malignant GCTB. DNA sequencing of the benign metastasising GCTB revealed no additional driver alterations; polyclonal seeding in the lung was identified, implying that the metastatic lesions represent an embolic event. Unsupervised clustering of DNA methylation profiles revealed that malignant H3.3-mutated tumours are distinct from their benign counterpart, and other bone tumours. Differential methylation analysis identified CCND1, encoding cyclin D1, as a plausible cancer driver gene in these tumours because hypermethylation of the CCND1 promoter was specific for GCTBs. We report here the genomic and methylation patterns underlying the rare clinical phenomena of benign metastasising and malignant transformation of GCTB and show how the combination of genomic and epigenomic findings could potentially distinguish benign from malignant GCTBs, thereby predicting aggressive behaviour in challenging diagnostic cases. © 2020 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.

Sarcomas are malignant soft tissue and bone tumours affecting adults, adolescents and children. They represent a morphologically heterogeneous class of tumours and some entities lack defining histopathological features. Therefore, the diagnosis of sarcomas is burdened with a high inter-observer variability and misclassification rate. Here, we demonstrate classification of soft tissue and bone tumours using a machine learning classifier algorithm based on array-generated DNA methylation data. This sarcoma classifier is trained using a dataset of 1077 methylation profiles from comprehensively pre-characterized cases comprising 62 tumour methylation classes constituting a broad range of soft tissue and bone sarcoma subtypes across the entire age spectrum. The performance is validated in a cohort of 428 sarcomatous tumours, of which 322 cases were classified by the sarcoma classifier. Our results demonstrate the potential of the DNA methylation-based sarcoma classification for research and future diagnostic applications.

Clear cell meningioma represents an uncommon variant of meningioma that typically affects children and young adults. Although an enrichment of loss-of-function mutations in the SMARCE1 gene has been reported for this subtype, comprehensive molecular investigations are lacking. Here we describe a molecularly distinct subset of tumors (n = 31), initially identified through genome-wide DNA methylation screening among a cohort of 3093 meningiomas, of which most were diagnosed histologically as clear cell meningioma. This cohort was further supplemented by an additional 11 histologically diagnosed clear cell meningiomas for analysis (n = 42). Targeted DNA sequencing revealed SMARCE1 mutations in 33/34 analyzed samples, accompanied by a nuclear loss of expression determined via immunohistochemistry and a decreased SMARCE1 transcript expression in the tumor cells. Analysis of time to progression or recurrence of patients within the clear cell meningioma group (n = 14) in comparison to those with meningioma WHO grade 2 (n = 220) revealed a similar outcome and support the assignment of WHO grade 2 to these tumors. Our findings indicate the existence of a highly distinct epigenetic signature of clear cell meningiomas, separate from all other variants of meningiomas, with recurrent mutations in the SMARCE1 gene. This suggests that these tumors may arise from a different precursor cell population than the broad spectrum of the other meningioma subtypes.

Ependymoma with RELA fusion has been defined as a novel entity of the revised World Health Organization 2016 classification of tumors of the central nervous system (CNS), characterized by fusion transcripts of the RELA gene and consequent pathological activation of the NFkB pathway. These tumors represent the majority of supratentorial ependymomas in children. The validation of diagnostic tools to identify this clinically relevant ependymoma entity is essential. Here, we have used interphase fluorescent in situ hybridization (FISH) for C11orf95 and RELA, immunohistochemistry (IHC) for p65-RelA and the recently developed DNA methylation-based classification besides conventional histopathology, and compared the precision of the methods in 40 supratentorial pediatric brain tumors diagnosed as ependymomas in the past years. Reverse transcription PCR (RT-PCR) and RNA sequencing were performed to explore discordant cases. Furthermore, we integrated imaging and clinical features as additional layers of information. The concordance between nuclear RelA expression by IHC and RELA FISH was 100%. Concordance between IHC and DNA methylation profiling, and between FISH and DNA methylation profiling was also high (96.4% and 95.2%, respectively). Thirty-four out of 40 (85%) cases were confirmed by integrated diagnoses as ependymal tumors, including 22 RELA-fused ependymomas (71% of ependymal tumors), two YAP1-fused ependymomas (6%), six non-RELA/non-YAP1 ependymomas (18%) and four ependymal/subependymal mixed tumors (12%). Ependymal/subependymal mixed tumors had an excellent clinical outcome despite the presence of histopathological signs of malignancy, suggesting that these tumors should not be diagnosed as classic ependymomas. DNA methylation profiling helped in the differential diagnosis of RELA-fused ependymomas. IHC and FISH, which are available in the majority of pathology laboratories, are valuable tools to identify RELA-fused ependymomas.

The current consensus recognizes four main medulloblastoma subgroups (wingless, Sonic hedgehog, group 3 and group 4). While medulloblastoma subgroups have been characterized extensively at the (epi-)genomic and transcriptomic levels, the proteome and phosphoproteome landscape remain to be comprehensively elucidated. Using quantitative (phospho)-proteomics in primary human medulloblastomas, we unravel distinct posttranscriptional regulation leading to highly divergent oncogenic signaling and kinase activity profiles in groups 3 and 4 medulloblastomas. Specifically, proteomic and phosphoproteomic analyses identify aberrant ERBB4-SRC signaling in group 4. Hence, enforced expression of an activated SRC combined with p53 inactivation induces murine tumors that resemble group 4 medulloblastoma. Therefore, our integrative proteogenomics approach unveils an oncogenic pathway and potential therapeutic vulnerability in the most common medulloblastoma subgroup.

Genomic characterization has begun to redefine diagnostic classifications of cancers. However, it remains a challenge to infer disease phenotypes from genomic alterations alone. To help realize the promise of genomics, we have performed a quantitative proteomics investigation using Stable Isotope Labeling by Amino Acids in Cell Culture (SILAC) and 41 tissue samples spanning the 4 genomically based subgroups of medulloblastoma and control cerebellum. We have identified and quantitated thousands of proteins across these groups and find that we are able to recapitulate the genomic subgroups based upon subgroup restricted and differentially abundant proteins while also identifying subgroup specific protein isoforms. Integrating our proteomic measurements with genomic data, we calculate a poor correlation between mRNA and protein abundance. Using EPIC 850 k methylation array data on the same tissues, we also investigate the influence of copy number alterations and DNA methylation on the proteome in an attempt to characterize the impact of these genetic features on the proteome. Reciprocally, we are able to use the proteome to identify which genomic alterations result in altered protein abundance and thus are most likely to impact biology. Finally, we are able to assemble protein-based pathways yielding potential avenues for clinical intervention. From these, we validate the EIF4F cap-dependent translation pathway as a novel druggable pathway in medulloblastoma. Thus, quantitative proteomics complements genomic platforms to yield a more complete understanding of functional tumor biology and identify novel therapeutic targets for medulloblastoma.

Human neural stem cell cultures provide progenitor cells that are potential cells of origin for brain cancers. However, the extent to which genetic predisposition to tumor formation can be faithfully captured in stem cell lines is uncertain. Here, we evaluated neuroepithelial stem (NES) cells, representative of cerebellar progenitors. We transduced NES cells with MYCN, observing medulloblastoma upon orthotopic implantation in mice. Significantly, transcriptomes and patterns of DNA methylation from xenograft tumors were globally more representative of human medulloblastoma compared to a MYCN-driven genetically engineered mouse model. Orthotopic transplantation of NES cells generated from Gorlin syndrome patients, who are predisposed to medulloblastoma due to germline-mutated PTCH1, also generated medulloblastoma. We engineered candidate cooperating mutations in Gorlin NES cells, with mutation of DDX3X or loss of GSE1 both accelerating tumorigenesis. These findings demonstrate that human NES cells provide a potent experimental resource for dissecting genetic causation in medulloblastoma.

The human and mouse genomes contain instructions that specify RNAs and proteins and govern the timing, magnitude, and cellular context of their production. To better delineate these elements, phase III of the Encyclopedia of DNA Elements (ENCODE) Project has expanded analysis of the cell and tissue repertoires of RNA transcription, chromatin structure and modification, DNA methylation, chromatin looping, and occupancy by transcription factors and RNA-binding proteins. Here we summarize these efforts, which have produced 5,992 new experimental datasets, including systematic determinations across mouse fetal development. All data are available through the ENCODE data portal (https://www.encodeproject.org), including phase II ENCODE1 and Roadmap Epigenomics2 data. We have developed a registry of 926,535 human and 339,815 mouse candidate cis-regulatory elements, covering 7.9 and 3.4% of their respective genomes, by integrating selected datatypes associated with gene regulation, and constructed a web-based server (SCREEN; http://screen.encodeproject.org) to provide flexible, user-defined access to this resource. Collectively, the ENCODE data and registry provide an expansive resource for the scientific community to build a better understanding of the organization and function of the human and mouse genomes.

A hallmark of pediatric low-grade glioma (pLGG) is aberrant signaling of the mitogen activated protein kinase (MAPK) pathway. Hence, inhibition of MAPK signaling using small molecule inhibitors such as MEK inhibitors (MEKi) may be a promising strategy.

Telomere maintenance by telomerase activation or alternative lengthening of telomeres (ALT) is a major determinant of poor outcome in neuroblastoma. Here, we screen for ALT in primary and relapsed neuroblastomas (n = 760) and characterize its features using multi-omics profiling. ALT-positive tumors are molecularly distinct from other neuroblastoma subtypes and enriched in a population-based clinical sequencing study cohort for relapsed cases. They display reduced ATRX/DAXX complex abundance, due to either ATRX mutations (55%) or low protein expression. The heterochromatic histone mark H3K9me3 recognized by ATRX is enriched at the telomeres of ALT-positive tumors. Notably, we find a high frequency of telomeric repeat loci with a neuroblastoma ALT-specific hotspot on chr1q42.2 and loss of the adjacent chromosomal segment forming a neo-telomere. ALT-positive neuroblastomas proliferate slowly, which is reflected by a protracted clinical course of disease. Nevertheless, children with an ALT-positive neuroblastoma have dismal outcome.

The identity of the cell of origin is a key determinant of cancer subtype, progression, and prognosis. Group 3 medulloblastoma (MB) is a malignant childhood brain cancer with poor prognosis and few candidates as putative cell of origin. We overexpressed the group 3 MB genetic drivers MYC and Gfi1 in different candidate cells of origin in the postnatal mouse cerebellum. We found that S100b+ cells are competent to initiate group 3 MB, and we observed that S100b+ cells have higher levels of Notch1 pathway activity compared to Math1+ cells. We found that additional activation of Notch1 in Math1+ and Sox2+ cells was sufficient to induce group 3 MB upon MYC/Gfi1 expression. Together, our data suggest that the Notch1 pathway plays a critical role in group 3 MB initiation.