Chemotherapies active in preclinical studies frequently fail in the clinic due to lack of efficacy, which limits progress for rare cancers since only small numbers of patients are available for clinical trials. Thus, a preclinical drug development pipeline was developed to prioritize potentially active regimens for pediatric brain tumors spanning from in vitro drug screening, through intracranial and intra-tumoral pharmacokinetics to in vivo efficacy studies. Here, as an example of the pipeline, data are presented for the combination of 5-fluoro-2'-deoxycytidine and tetrahydrouridine in three pediatric brain tumor models. The in vitro activity of nine novel therapies was tested against tumor spheres derived from faithful mouse models of Group 3 medulloblastoma, ependymoma, and choroid plexus carcinoma. Agents with the greatest in vitro potency were then subjected to a comprehensive series of in vivo pharmacokinetic (PK) and pharmacodynamic (PD) studies culminating in preclinical efficacy trials in mice harboring brain tumors. The nucleoside analog 5-fluoro-2'-deoxycytidine (FdCyd) markedly reduced the proliferation in vitro of all three brain tumor cell types at nanomolar concentrations. Detailed intracranial PK studies confirmed that systemically administered FdCyd exceeded concentrations in brain tumors necessary to inhibit tumor cell proliferation, but no tumor displayed a significant in vivo therapeutic response. Despite promising in vitro activity and in vivo PK properties, FdCyd is unlikely to be an effective treatment of pediatric brain tumors, and therefore was deprioritized for the clinic. Our comprehensive and integrated preclinical drug development pipeline should reduce the attrition of drugs in clinical trials.

TET-family dioxygenases oxidize 5-methylcytosine (5mC) in DNA, and exert tumour suppressor activity in many types of cancers. Even in the absence of TET coding region mutations, TET loss-of-function is strongly associated with cancer. Here we show that acute elimination of TET function induces the rapid development of an aggressive, fully-penetrant and cell-autonomous myeloid leukaemia in mice, pointing to a causative role for TET loss-of-function in this myeloid malignancy. Phenotypic and transcriptional profiling shows aberrant differentiation of haematopoietic stem/progenitor cells, impaired erythroid and lymphoid differentiation and strong skewing to the myeloid lineage, with only a mild relation to changes in DNA modification. We also observe progressive accumulation of phospho-H2AX and strong impairment of DNA damage repair pathways, suggesting a key role for TET proteins in maintaining genome integrity.

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.

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.

Gains and losses in DNA methylation are prominent features of mammalian cell types. To gain insight into the mechanisms that promote shifts in DNA methylation and contribute to changes in cell fate, including malignant transformation, we performed genome-wide mapping of 5-methylcytosine and 5-hydroxymethylcytosine in purified mouse hematopoietic stem cells. We discovered extended regions of low methylation (canyons) that span conserved domains frequently containing transcription factors and are distinct from CpG islands and shores. About half of the genes in these methylation canyons are coated with repressive histone marks, whereas the remainder are covered by activating histone marks and are highly expressed in hematopoietic stem cells (HSCs). Canyon borders are demarked by 5-hydroxymethylcytosine and become eroded in the absence of DNA methyltransferase 3a (Dnmt3a). Genes dysregulated in human leukemias are enriched for canyon-associated genes. The new epigenetic landscape we describe may provide a mechanism for the regulation of hematopoiesis and may contribute to leukemia development.

DNA enrichment followed by sequencing is a versatile tool in molecular biology, with a wide variety of applications including genome-wide analysis of epigenetic marks and mechanisms. A common requirement of these diverse applications is a comparison of read coverage between experimental conditions. The amount of samples generated for such comparisons ranges from few replicates to hundreds of samples per condition for epigenome-wide association studies. Consequently, there is an urgent need for software that allows for fast and simple processing and comparison of sequencing data derived from enriched DNA.

Transcription factors are proteins that regulate gene expression by binding to cis-regulatory sequences such as promoters and enhancers. In embryonic stem (ES) cells, binding of the transcription factors OCT4, SOX2 and NANOG is essential to maintain the capacity of the cells to differentiate into any cell type of the developing embryo. It is known that transcription factors interact to regulate gene expression. In this study we show that combinatorial binding is strongly associated with co-localization of the transcriptional co-activator Mediator, H3K27ac and increased expression of nearby genes in embryonic stem cells. We observe that the same loci bound by Oct4, Nanog and Sox2 in ES cells frequently drive expression in early embryonic development. Comparison of mouse and human ES cells shows that less than 5% of individual binding events for OCT4, SOX2 and NANOG are shared between species. In contrast, about 15% of combinatorial binding events and even between 53% and 63% of combinatorial binding events at enhancers active in early development are conserved. Our analysis suggests that the combination of OCT4, SOX2 and NANOG binding is critical for transcription in ES cells and likely plays an important role for embryogenesis by binding at conserved early developmental enhancers. Our data suggests that the fast evolutionary rewiring of regulatory networks mainly affects individual binding events, whereas "gene regulatory hotspots" which are bound by multiple factors and active in multiple tissues throughout early development are under stronger evolutionary constraints.

Choroid plexus carcinoma (CPC) is a rare and aggressive tumor of infancy without a clear treatment strategy. This study describes the outcomes of children with CPC treated on the multi-institutional phase 2 SJYC07 trial and reports on the significance of clinical and molecular characteristics.

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.

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.

Medulloblastoma is the most common malignant brain tumor of childhood and is considered a tumor with low mutational burden (~1 Mut/Mb). Therefore, though the medulloblastoma genomes have been extensively characterized in literature, reports on potential hypermutations and underlying mutagenic processes in medulloblastomas are limited.

Taking a genome-wide association study approach, we identified inherited genetic variations in ACYP2 associated with cisplatin-related ototoxicity (rs1872328: P = 3.9 × 10(-8), hazard ratio = 4.5) in 238 children with newly diagnosed brain tumors, with independent replication in 68 similarly treated children. The ACYP2 risk variant strongly predisposed these patients to precipitous hearing loss and was related to ototoxicity severity. These results point to new biology underlying the ototoxic effects of platinum agents.

TET/JBP enzymes oxidize 5-methylpyrimidines in DNA. In mammals, the oxidized methylcytosines (oxi-mCs) function as epigenetic marks and likely intermediates in DNA demethylation. Here we present a method based on diglucosylation of 5-hydroxymethylcytosine (5hmC) to simultaneously map 5hmC, 5-formylcytosine, and 5-carboxylcytosine at near-base-pair resolution. We have used the method to map the distribution of oxi-mC across the genome of Coprinopsis cinerea, a basidiomycete that encodes 47 TET/JBP paralogs in a previously unidentified class of DNA transposons. Like 5-methylcytosine residues from which they are derived, oxi-mC modifications are enriched at centromeres, TET/JBP transposons, and multicopy paralogous genes that are not expressed, but rarely mark genes whose expression changes between two developmental stages. Our study provides evidence for the emergence of an epigenetic regulatory system through recruitment of selfish elements in a eukaryotic lineage, and describes a method to map all three different species of oxi-mCs simultaneously.

The transcription factor OCT4 is highly expressed in pluripotent embryonic stem cells which are derived from the inner cell mass of mammalian blastocysts. Pluripotency and self renewal are controlled by a transcription regulatory network governed by the transcription factors OCT4, SOX2 and NANOG. Recent studies on reprogramming somatic cells to induced pluripotent stem cells highlight OCT4 as a key regulator of pluripotency.

Genome-wide association studies (GWASs) identify single nucleotide polymorphisms (SNPs) that are enriched in individuals suffering from a given disease. Most disease-associated SNPs fall into non-coding regions, so that it is not straightforward to infer phenotype or function; moreover, many SNPs are in tight genetic linkage, so that a SNP identified as associated with a particular disease may not itself be causal, but rather signify the presence of a linked SNP that is functionally relevant to disease pathogenesis. Here, we present an analysis method that takes advantage of the recent rapid accumulation of epigenomics data to address these problems for some SNPs. Using asthma as a prototypic example; we show that non-coding disease-associated SNPs are enriched in genomic regions that function as regulators of transcription, such as enhancers and promoters. Identifying enhancers based on the presence of the histone modification marks such as H3K4me1 in different cell types, we show that the location of enhancers is highly cell-type specific. We use these findings to predict which SNPs are likely to be directly contributing to disease based on their presence in regulatory regions, and in which cell types their effect is expected to be detectable. Moreover, we can also predict which cell types contribute to a disease based on overlap of the disease-associated SNPs with the locations of enhancers present in a given cell type. Finally, we suggest that it will be possible to re-analyze GWAS studies with much higher power by limiting the SNPs considered to those in coding or regulatory regions of cell types relevant to a given disease.

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.

Pediatric high-grade glioma (pHGG) is a major contributor to cancer-related death in children. In vitro and in vivo disease models reflecting the intimate connection between developmental context and pathogenesis of pHGG are essential to advance understanding and identify therapeutic vulnerabilities. Here we report establishment of 21 patient-derived pHGG orthotopic xenograft (PDOX) models and eight matched cell lines from diverse groups of pHGG. These models recapitulate histopathology, DNA methylation signatures, mutations and gene expression patterns of the patient tumors from which they were derived, and include rare subgroups not well-represented by existing models. We deploy 16 new and existing cell lines for high-throughput screening (HTS). In vitro HTS results predict variable in vivo response to PI3K/mTOR and MEK pathway inhibitors. These unique new models and an online interactive data portal for exploration of associated detailed molecular characterization and HTS chemical sensitivity data provide a rich resource for pediatric brain tumor research.

Cerebellar granule cells (GCs) are cells which comprise over 50% of the neurons in the entire nervous system. GCs enable the cerebellum to properly regulate motor coordination, learning, and consolidation, in addition to cognition, emotion and language. During GC development, maternal GC progenitors (GCPs) divide to produce not only postmitotic GCs but also sister GCPs. However, the molecular machinery for regulating the proportional production of distinct sister cell types from seemingly uniform GCPs is not yet fully understood. Here we report that Notch signaling creates a distinction between GCPs and leads to their proportional differentiation in mice. Among Notch-related molecules, Notch1, Notch2, Jag1, and Hes1 are prominently expressed in GCPs. In vivo monitoring of Hes1-promoter activities showed the presence of two types of GCPs, Notch-signaling ON and OFF, in the external granule layer (EGL). Single-cell RNA sequencing (scRNA-seq) and in silico analyses indicate that ON-GCPs have more proliferative and immature properties, while OFF-GCPs have opposite characteristics. Overexpression as well as knock-down (KD) experiments using in vivo electroporation showed that NOTCH2 and HES1 are involved cell-autonomously to suppress GCP differentiation by inhibiting NEUROD1 expression. In contrast, JAG1-expressing cells non-autonomously upregulated Notch signaling activities via NOTCH2-HES1 in surrounding GCPs, eventually suppressing their differentiation. These findings suggest that Notch signaling results in the proportional generation of two types of cells, immature and differentiating GCPs, which contributes to the well-organized differentiation of GCs.

In 2012, an international consensus paper reported that medulloblastoma comprises four molecular subgroups (WNT, SHH, Group 3, and Group 4), each associated with distinct genomic features and clinical behavior. Independently, multiple recent reports have defined further intra-subgroup heterogeneity in the form of biologically and clinically relevant subtypes. However, owing to differences in patient cohorts and analytical methods, estimates of subtype number and definition have been inconsistent, especially within Group 3 and Group 4. Herein, we aimed to reconcile the definition of Group 3/Group 4 MB subtypes through the analysis of a series of 1501 medulloblastomas with DNA-methylation profiling data, including 852 with matched transcriptome data. Using multiple complementary bioinformatic approaches, we compared the concordance of subtype calls between published cohorts and analytical methods, including assessments of class-definition confidence and reproducibility. While the lowest complexity solutions continued to support the original consensus subgroups of Group 3 and Group 4, our analysis most strongly supported a definition comprising eight robust Group 3/Group 4 subtypes (types I-VIII). Subtype II was consistently identified across all component studies, while all others were supported by multiple class-definition methods. Regardless of analytical technique, increasing cohort size did not further increase the number of identified Group 3/Group 4 subtypes. Summarizing the molecular and clinico-pathological features of these eight subtypes indicated enrichment of specific driver gene alterations and cytogenetic events amongst subtypes, and identified highly disparate survival outcomes, further supporting their biological and clinical relevance. Collectively, this study provides continued support for consensus Groups 3 and 4 while enabling robust derivation of, and categorical accounting for, the extensive intertumoral heterogeneity within Groups 3 and 4, revealed by recent high-resolution subclassification approaches. Furthermore, these findings provide a basis for application of emerging methods (e.g., proteomics/single-cell approaches) which may additionally inform medulloblastoma subclassification. Outputs from this study will help shape definition of the next generation of medulloblastoma clinical protocols and facilitate the application of enhanced molecularly guided risk stratification to improve outcomes and quality of life for patients and their families.

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.