Genome-wide enrichment of methylated DNA followed by sequencing (MeDIP-seq) offers a reasonable compromise between experimental costs and genomic coverage. However, the computational analysis of these experiments is complex, and quantification of the enrichment signals in terms of absolute levels of methylation requires specific transformation. In this work, we present QSEA, Quantitative Sequence Enrichment Analysis, a comprehensive workflow for the modelling and subsequent quantification of MeDIP-seq data. As the central part of the workflow we have developed a Bayesian statistical model that transforms the enrichment read counts to absolute levels of methylation and, thus, enhances interpretability and facilitates comparison with other methylation assays. We suggest several calibration strategies for the critical parameters of the model, either using additional data or fairly general assumptions. By comparing the results with bisulfite sequencing (BS) validation data, we show the improvement of QSEA over existing methods. Additionally, we generated a clinically relevant benchmark data set consisting of methylation enrichment experiments (MeDIP-seq), BS-based validation experiments (Methyl-seq) as well as gene expression experiments (RNA-seq) derived from non-small cell lung cancer patients, and show that the workflow retrieves well-known lung tumour methylation markers that are causative for gene expression changes, demonstrating the applicability of QSEA for clinical studies. QSEA is implemented in R and available from the Bioconductor repository 3.4 (www.bioconductor.org/packages/qsea).

We present a generative model, Lux, to quantify DNA methylation modifications from any combination of bisulfite sequencing approaches, including reduced, oxidative, TET-assisted, chemical-modification assisted, and methylase-assisted bisulfite sequencing data. Lux models all cytosine modifications (C, 5mC, 5hmC, 5fC, and 5caC) simultaneously together with experimental parameters, including bisulfite conversion and oxidation efficiencies, as well as various chemical labeling and protection steps. We show that Lux improves the quantification and comparison of cytosine modification levels and that Lux can process any oxidized methylcytosine sequencing data sets to quantify all cytosine modifications. Analysis of targeted data from Tet2-knockdown embryonic stem cells and T cells during development demonstrates DNA modification quantification at unprecedented detail, quantifies active demethylation pathways and reveals 5hmC localization in putative regulatory regions.

Medulloblastoma is associated with rare hereditary cancer predisposition syndromes; however, consensus medulloblastoma predisposition genes have not been defined and screening guidelines for genetic counselling and testing for paediatric patients are not available. We aimed to assess and define these genes to provide evidence for future screening guidelines.

Current therapies for medulloblastoma, a highly malignant childhood brain tumour, impose debilitating effects on the developing child, and highlight the need for molecularly targeted treatments with reduced toxicity. Previous studies have been unable to identify the full spectrum of driver genes and molecular processes that operate in medulloblastoma subgroups. Here we analyse the somatic landscape across 491 sequenced medulloblastoma samples and the molecular heterogeneity among 1,256 epigenetically analysed cases, and identify subgroup-specific driver alterations that include previously undiscovered actionable targets. Driver mutations were confidently assigned to most patients belonging to Group 3 and Group 4 medulloblastoma subgroups, greatly enhancing previous knowledge. New molecular subtypes were differentially enriched for specific driver events, including hotspot in-frame insertions that target KBTBD4 and 'enhancer hijacking' events that activate PRDM6. Thus, the application of integrative genomics to an extensive cohort of clinical samples derived from a single childhood cancer entity revealed a series of cancer genes and biologically relevant subtype diversity that represent attractive therapeutic targets for the treatment of patients with medulloblastoma.

Illumina Infinium DNA methylation arrays are a cost-effective technology to measure DNA methylation at CpG sites genome-wide and across cohorts of normal and cancer samples. While copy number alterations are commonly inferred from array-CGH, SNP arrays, or whole-genome DNA sequencing, Illumina Infinium DNA methylation arrays have been shown to detect copy number alterations at comparable sensitivity. Here we present an accessible, interactive GenePattern notebook for the analysis of copy number variation using Illumina Infinium DNA methylation arrays. The notebook provides a graphical user interface to a workflow using the R/Bioconductor packages minfi and conumee. The environment allows analysis to be performed without the installation of the R software environment, the packages and dependencies, and without the need to write or manipulate code.

The primary cilium, a signaling organelle projecting from the surface of a cell, controls cellular physiology and behavior. The presence or absence of primary cilia is a distinctive feature of a given tumor type; however, whether and how the primary cilium contributes to tumorigenesis are unknown for most tumors. Medulloblastoma (MB) is a common pediatric brain cancer comprising four groups: SHH, WNT, group 3 (G3), and group 4 (G4). From 111 cases of MB, we show that primary cilia are abundant in SHH and WNT MBs but rare in G3 and G4 MBs. Using WNT and G3 MB mouse models, we show that primary cilia promote WNT MB by facilitating translation of mRNA encoding β-catenin, a major oncoprotein driving WNT MB, whereas cilium loss promotes G3 MB by disrupting cell cycle control and destabilizing the genome. Our findings reveal tumor type-specific ciliary functions and underlying molecular mechanisms. Moreover, we expand the function of primary cilia to translation control and reveal a molecular mechanism by which cilia regulate cell cycle progression, thereby providing new frameworks for studying cilium function in normal and pathologic conditions.

Differentiating tumor recurrence or progression from pseudoprogression during surveillance of pediatric high-grade gliomas (PHGGs) using MRI, the primary imaging modality for evaluation of brain tumors, can be challenging. The aim of this study was to evaluate whether 11C-methionine PET, a molecular imaging technique that detects functionally active tumors, is useful for further evaluating MRI changes concerning for tumor recurrence during routine surveillance. Methods: Using 11C-methionine PET during follow-up visits, we evaluated 27 lesions in 26 patients with new or worsening MRI abnormalities for whom tumor recurrence was of concern. We performed quantitative and qualitative assessments of both 11C-methionine PET and MRI data to predict the presence of tumor recurrence. Further, to assess for an association with overall survival (OS), we plotted the time from development of the imaging changes against survival. Results: Qualitative evaluation of 11C-methionine PET achieved 100% sensitivity, 60% specificity, and 93% accuracy to correctly predict the presence of tumors in 27 new or worsening MRI abnormalities. Qualitative MRI evaluation achieved sensitivity ranging from 86% to 95%, specificity ranging from 40% to 60%, and accuracy ranging from 85% to 89%. The interobserver agreement for 11C-methionine PET assessment was 100%, whereas the interobserver agreement was only 50% for MRI (P < 0.01). Quantitative MRI and 11C-methionine PET evaluation using receiver-operating characteristics demonstrated higher specificity (80%) than did qualitative evaluations (40%-60%). Postcontrast enhancement volume, metabolic tumor volume, tumor-to-brain ratio, and presence of tumor as determined by consensus MRI assessment were inversely associated with OS. Conclusion:11C-methionine PET has slightly higher sensitivity and accuracy for correctly predicting tumor recurrence, with excellent interobserver agreement, than does MRI. Quantitative 11C-methionine PET can also predict OS. These findings suggest that 11C-methionine PET can be useful for further evaluation of MRI changes during surveillance of previously treated PHGGs.

In response to acute infection, naive CD8+ T cells expand, differentiate into effector cells, and then contract to a long-lived pool of memory cells after pathogen clearance. During chronic infections or in tumors, CD8+ T cells acquire an "exhausted" phenotype. Here we present genome-wide comparisons of chromatin accessibility and gene expression from endogenous CD8+ T cells responding to acute and chronic viral infection using ATAC-seq and RNA-seq techniques. Acquisition of effector, memory, or exhausted phenotypes was associated with stable changes in chromatin accessibility away from the naive T cell state. Regions differentially accessible between functional subsets in vivo were enriched for binding sites of transcription factors known to regulate these subsets, including E2A, BATF, IRF4, T-bet, and TCF1. Exhaustion-specific accessible regions were enriched for consensus binding sites for NFAT and Nr4a family members, indicating that chronic stimulation confers a unique accessibility profile on exhausted cells.

Pineoblastoma is a rare and highly aggressive brain cancer of childhood, histologically belonging to the spectrum of primitive neuroectodermal tumors. Patients with germline mutations in DICER1, a ribonuclease involved in microRNA processing, have increased risk of pineoblastoma, but genetic drivers of sporadic pineoblastoma remain unknown. Here, we analyzed pediatric and adult pineoblastoma samples (n = 23) using a combination of genome-wide DNA methylation profiling and whole-exome sequencing or whole-genome sequencing. Pediatric and adult pineoblastomas showed distinct methylation profiles, the latter clustering with lower-grade pineal tumors and normal pineal gland. Recurrent variants were found in genes involved in PKA- and NF-κB signaling, as well as in chromatin remodeling genes. We identified recurrent homozygous deletions of DROSHA, acting upstream of DICER1 in microRNA processing, and a novel microduplication involving chromosomal region 1q21 containing PDE4DIP (myomegalin), comprising the ancient DUF1220 protein domain. Expresion of PDE4DIP and DUF1220 proteins was present exclusively in pineoblastoma with PDE4DIP gain.

Medulloblastoma is curable in approximately 70% of patients. Over the past decade, progress in improving survival using conventional therapies has stalled, resulting in reduced quality of life due to treatment-related side effects, which are a major concern in survivors. The vast amount of genomic and molecular data generated over the last 5-10 years encourages optimism that improved risk stratification and new molecular targets will improve outcomes. It is now clear that medulloblastoma is not a single-disease entity, but instead consists of at least four distinct molecular subgroups: WNT/Wingless, Sonic Hedgehog, Group 3, and Group 4. The Medulloblastoma Down Under 2013 meeting, which convened at Bunker Bay, Australia, brought together 50 leading clinicians and scientists. The 2-day agenda included focused sessions on pathology and molecular stratification, genomics and mouse models, high-throughput drug screening, and clinical trial design. The meeting established a global action plan to translate novel biologic insights and drug targeting into treatment regimens to improve outcomes. A consensus was reached in several key areas, with the most important being that a novel classification scheme for medulloblastoma based on the four molecular subgroups, as well as histopathologic features, should be presented for consideration in the upcoming fifth edition of the World Health Organization's classification of tumours of the central nervous system. Three other notable areas of agreement were as follows: (1) to establish a central repository of annotated mouse models that are readily accessible and freely available to the international research community; (2) to institute common eligibility criteria between the Children's Oncology Group and the International Society of Paediatric Oncology Europe and initiate joint or parallel clinical trials; (3) to share preliminary high-throughput screening data across discovery labs to hasten the development of novel therapeutics. Medulloblastoma Down Under 2013 was an effective forum for meaningful discussion, which resulted in enhancing international collaborative clinical and translational research of this rare disease. This template could be applied to other fields to devise global action plans addressing all aspects of a disease, from improved disease classification, treatment stratification, and drug targeting to superior treatment regimens to be assessed in cooperative international clinical trials.

We devised a high-throughput, cell-based assay to identify compounds to treat Group3 medulloblastoma (G3 MB). Mouse G3 MBs neurospheres were screened against a library of approximately 7,000 compounds including US Food and Drug Administration-approved drugs. We found that pemetrexed and gemcitabine preferentially inhibited G3 MB proliferation in vitro compared to control neurospheres and substantially inhibited G3 MB proliferation in vivo. When combined, these two drugs significantly increased survival of mice bearing cortical implants of mouse and human G3 MBs that overexpress MYC compared to each agent alone, while having little effect on mouse MBs of the sonic hedgehog subgroup. Our findings strongly suggest that combination therapy with pemetrexed and gemcitabine is a promising treatment for G3 MBs.

B-cell fate determination requires the action of transcription factors that operate in a regulatory network to activate B-lineage genes and repress lineage-inappropriate genes. However, the dynamics and hierarchy of events in B-cell programming remain obscure. To uncouple the dynamics of transcription factor expression from functional consequences, we generated induction systems in developmentally arrested Ebf1-/- pre-pro-B cells to allow precise experimental control of EBF1 expression in the genomic context of progenitor cells. Consistent with the described role of EBF1 as a pioneer transcription factor, we show in a time-resolved analysis that EBF1 occupancy coincides with EBF1 expression and precedes the formation of chromatin accessibility. We observed dynamic patterns of EBF1 target gene expression and sequential up-regulation of transcription factors that expand the regulatory network at the pro-B-cell stage. A continuous EBF1 function was found to be required for Cd79a promoter activity and for the maintenance of an accessible chromatin domain that is permissive for binding of other transcription factors. Notably, transient EBF1 occupancy was detected at lineage-inappropriate genes prior to their silencing in pro-B cells. Thus, persistent and transient functions of EBF1 allow for an ordered sequence of epigenetic and transcriptional events in B-cell programming.

Genomic sequencing has driven precision-based oncology therapy; however, the genetic drivers of many malignancies remain unknown or non-targetable, so alternative approaches to the identification of therapeutic leads are necessary. Ependymomas are chemotherapy-resistant brain tumours, which, despite genomic sequencing, lack effective molecular targets. Intracranial ependymomas are segregated on the basis of anatomical location (supratentorial region or posterior fossa) and further divided into distinct molecular subgroups that reflect differences in the age of onset, gender predominance and response to therapy. The most common and aggressive subgroup, posterior fossa ependymoma group A (PF-EPN-A), occurs in young children and appears to lack recurrent somatic mutations. Conversely, posterior fossa ependymoma group B (PF-EPN-B) tumours display frequent large-scale copy number gains and losses but have favourable clinical outcomes. More than 70% of supratentorial ependymomas are defined by highly recurrent gene fusions in the NF-κB subunit gene RELA (ST-EPN-RELA), and a smaller number involve fusion of the gene encoding the transcriptional activator YAP1 (ST-EPN-YAP1). Subependymomas, a distinct histologic variant, can also be found within the supratetorial and posterior fossa compartments, and account for the majority of tumours in the molecular subgroups ST-EPN-SE and PF-EPN-SE. Here we describe mapping of active chromatin landscapes in 42 primary ependymomas in two non-overlapping primary ependymoma cohorts, with the goal of identifying essential super-enhancer-associated genes on which tumour cells depend. Enhancer regions revealed putative oncogenes, molecular targets and pathways; inhibition of these targets with small molecule inhibitors or short hairpin RNA diminished the proliferation of patient-derived neurospheres and increased survival in mouse models of ependymomas. Through profiling of transcriptional enhancers, our study provides a framework for target and drug discovery in other cancers that lack known genetic drivers and are therefore difficult to treat.

Mechanistic target of rapamycin (MTOR) cooperates with Hedgehog (HH) signaling, but the underlying mechanisms are incompletely understood. Here we provide genetic, biochemical, and pharmacologic evidence that MTOR complex 1 (mTORC1)-dependent translation is a prerequisite for HH signaling. The genetic loss of mTORC1 function inhibited HH signaling-driven growth of the cerebellum and medulloblastoma. Inhibiting translation or mTORC1 blocked HH signaling. Depleting 4EBP1, an mTORC1 target that inhibits translation, alleviated the dependence of HH signaling on mTORC1. Consistent with this, phosphorylated 4EBP1 levels were elevated in HH signaling-driven medulloblastomas in mice and humans. In mice, an mTORC1 inhibitor suppressed medulloblastoma driven by a mutant SMO that is inherently resistant to existing SMO inhibitors, prolonging the survival of the mice. Our study reveals that mTORC1-mediated translation is a key component of HH signaling and an important target for treating medulloblastoma and other cancers driven by HH signaling.

It is essential to understand the network of transcription factors controlling self-renewal of human embryonic stem cells (ESCs) and human embryonal carcinoma cells (ECs) if we are to exploit these cells in regenerative medicine regimes. Correlating gene expression levels after RNAi-based ablation of OCT4 function with its downstream targets enables a better prediction of motif-specific driven expression modules pertinent for self-renewal and differentiation of embryonic stem cells and induced pluripotent stem cells.We initially identified putative direct downstream targets of OCT4 by employing CHIP-on-chip analysis. A comparison of three peak analysis programs revealed a refined list of OCT4 targets in the human EC cell line NCCIT, this list was then compared to previously published OCT4 CHIP-on-chip datasets derived from both ES and EC cells. We have verified an enriched POU-motif, discovered by a de novo approach, thus enabling us to define six distinct modules of OCT4 binding and regulation of its target genes.A selection of these targets has been validated, like NANOG, which harbours the evolutionarily conserved OCT4-SOX2 binding motif within its proximal promoter. Other validated targets, which do not harbour the classical HMG motif are USP44 and GADD45G, a key regulator of the cell cycle. Over-expression of GADD45G in NCCIT cells resulted in an enrichment and up-regulation of genes associated with the cell cycle (CDKN1B, CDKN1C, CDK6 and MAPK4) and developmental processes (BMP4, HAND1, EOMES, ID2, GATA4, GATA5, ISL1 and MSX1). A comparison of positively regulated OCT4 targets common to EC and ES cells identified genes such as NANOG, PHC1, USP44, SOX2, PHF17 and OCT4, thus further confirming their universal role in maintaining self-renewal in both cell types. Finally we have created a user-friendly database (http://biit.cs.ut.ee/escd/), integrating all OCT4 and stem cell related datasets in both human and mouse ES and EC cells.In the current era of systems biology driven research, we envisage that our integrated embryonic stem cell database will prove beneficial to the booming field of ES, iPS and cancer research.

Glioblastoma (GBM) is the most lethal primary brain cancer characterized by therapeutic resistance, which is promoted by GBM stem cells (GSC). Here, we interrogated gene expression and whole-genome CRISPR/Cas9 screening in a large panel of patient-derived GSCs, differentiated GBM cells (DGC), and neural stem cells (NSC) to identify master regulators of GSC stemness, revealing an essential transcription state with increased RNA polymerase II-mediated transcription. The YY1 and transcriptional CDK9 complex was essential for GSC survival and maintenance in vitro and in vivo. YY1 interacted with CDK9 to regulate transcription elongation in GSCs. Genetic or pharmacologic targeting of the YY1-CDK9 complex elicited RNA m6A modification-dependent interferon responses, reduced regulatory T-cell infiltration, and augmented efficacy of immune checkpoint therapy in GBM. Collectively, these results suggest that YY1-CDK9 transcription elongation complex defines a targetable cell state with active transcription, suppressed interferon responses, and immunotherapy resistance in GBM. SIGNIFICANCE: Effective strategies to rewire immunosuppressive microenvironment and enhance immunotherapy response are still lacking in GBM. YY1-driven transcriptional elongation machinery represents a druggable target to activate interferon response and enhance anti-PD-1 response through regulating the m6A modification program, linking epigenetic regulation to immunomodulatory function in GBM.This article is highlighted in the In This Issue feature, p. 275.

Group3 (G3) medulloblastoma (MB) is one of the deadliest forms of the disease for which novel treatment is desperately needed. Here we evaluate ribociclib, a highly selective CDK4/6 inhibitor, with gemcitabine in mouse and human G3MBs. Ribociclib central nervous system (CNS) penetration was assessed by in vivo microdialysis and by IHC and gene expression studies and found to be CNS-penetrant. Tumors from mice treated with short term oral ribociclib displayed inhibited RB phosphorylation, downregulated E2F target genes, and decreased proliferation. Survival studies to determine the efficacy of ribociclib and gemcitabine combination were performed on mice intracranially implanted with luciferase-labeled mouse and human G3MBs. Treatment of mice with the combination of ribociclib and gemcitabine was well tolerated, slowed tumor progression and metastatic spread, and increased survival. Expression-based gene activity and cell state analysis investigated the effects of the combination after short- and long-term treatments. Molecular analysis of treated versus untreated tumors showed a significant decrease in the activity and expression of genes involved in cell-cycle progression and DNA damage response, and an increase in the activity and expression of genes implicated in neuronal identity and neuronal differentiation. Our findings in both mouse and human patient-derived orthotopic xenograft models suggest that ribociclib and gemcitabine combination therapy warrants further investigation as a treatment strategy for children with G3MB.

Ependymoma is a tumor of the brain or spinal cord. The two most common and aggressive molecular groups of ependymoma are the supratentorial ZFTA-fusion associated and the posterior fossa ependymoma group A. In both groups, tumors occur mainly in young children and frequently recur after treatment. Although molecular mechanisms underlying these diseases have recently been uncovered, they remain difficult to target and innovative therapeutic approaches are urgently needed. Here, we use genome-wide chromosome conformation capture (Hi-C), complemented with CTCF and H3K27ac ChIP-seq, as well as gene expression and DNA methylation analysis in primary and relapsed ependymoma tumors, to identify chromosomal conformations and regulatory mechanisms associated with aberrant gene expression. In particular, we observe the formation of new topologically associating domains ('neo-TADs') caused by structural variants, group-specific 3D chromatin loops, and the replacement of CTCF insulators by DNA hyper-methylation. Through inhibition experiments, we validate that genes implicated by these 3D genome conformations are essential for the survival of patient-derived ependymoma models in a group-specific manner. Thus, this study extends our ability to reveal tumor-dependency genes by 3D genome conformations even in tumors that lack targetable genetic alterations.

Aberrant CpG methylation is a universal epigenetic trait of cancer cell genomes. However, human cancer samples or cell lines preclude the investigation of epigenetic changes occurring early during tumour development. Here, we have used MeDIP-seq to analyse the DNA methylome of APC(Min) adenoma as a model for intestinal cancer initiation, and we present a list of more than 13,000 recurring differentially methylated regions (DMRs) characterizing intestinal adenoma of the mouse. We show that Polycomb Repressive Complex (PRC) targets are strongly enriched among hypermethylated DMRs, and several PRC2 components and DNA methyltransferases were up-regulated in adenoma. We further demonstrate by bisulfite pyrosequencing of purified cell populations that the DMR signature arises de novo in adenoma cells rather than by expansion of a pre-existing pattern in intestinal stem cells or undifferentiated crypt cells. We found that epigenetic silencing of tumour suppressors, which occurs frequently in colon cancer, was rare in adenoma. Quite strikingly, we identified a core set of DMRs, which is conserved between mouse adenoma and human colon cancer, thus possibly revealing a global panel of epigenetically modified genes for intestinal tumours. Our data allow a distinction between early conserved epigenetic alterations occurring in intestinal adenoma and late stochastic events promoting colon cancer progression, and may facilitate the selection of more specific clinical epigenetic biomarkers.

Asthma and autoimmune disease susceptibility has been strongly linked to genetic variants in the 17q21 haploblock that alter the expression of ORMDL3; however, the molecular mechanisms by which these variants perturb gene expression and the cell types in which this effect is most prominent are unclear. We found several 17q21 variants overlapped enhancers present mainly in primary immune cell types. CD4+ T cells showed the greatest increase (threefold) in ORMDL3 expression in individuals carrying the asthma-risk alleles, where ORMDL3 negatively regulated interleukin-2 production. The asthma-risk variants rs4065275 and rs12936231 switched CTCF-binding sites in the 17q21 locus, and 4C-Seq assays showed that several distal cis-regulatory elements upstream of the disrupted ZPBP2 CTCF-binding site interacted with the ORMDL3 promoter region in CD4+ T cells exclusively from subjects carrying asthma-risk alleles. Overall, our results suggested that T cells are one of the most prominent cell types affected by 17q21 variants.