We present Model-based Analysis of ChIP-Seq data, MACS, which analyzes data generated by short read sequencers such as Solexa's Genome Analyzer. MACS empirically models the shift size of ChIP-Seq tags, and uses it to improve the spatial resolution of predicted binding sites. MACS also uses a dynamic Poisson distribution to effectively capture local biases in the genome, allowing for more robust predictions. MACS compares favorably to existing ChIP-Seq peak-finding algorithms, and is freely available.

1-methyl-D-tryptophan (1-D-MT) is currently being used in clinical trials in patients with relapsed or refractory solid tumors with the aim of inhibiting indoleamine-2,3-dioxygenase (IDO)-mediated tumor immune escape. IDO is expressed in tumors and tumor-draining lymph nodes and degrades tryptophan (trp) to create an immunsuppressive micromilieu both by depleting trp and by accumulating immunosuppressive metabolites of the kynurenine (kyn) pathway. Here we show that proliferation of alloreactive T-cells cocultured with IDO1-positive human cancer cells paradoxically was inhibited by 1-D-MT. Surprisingly incubation with 1-D-MT increased kyn production of human cancer cells. Cell-free assays revealed that 1-D-MT did not alter IDO1 enzymatic activity. Instead, 1-D-MT induced IDO1 mRNA and protein expression through pathways involving p38 MAPK and JNK signalling. Treatment of cancer patients with 1-D-MT has transcriptional effects that may promote rather than suppress anti-tumor immune escape by increasing IDO1 in the cancer cells. These off-target effects should be carefully analyzed in the ongoing clinical trials with 1-D-MT.

Chromatin profiling has emerged as a powerful means of genome annotation and detection of regulatory activity. The approach is especially well suited to the characterization of non-coding portions of the genome, which critically contribute to cellular phenotypes yet remain largely uncharted. Here we map nine chromatin marks across nine cell types to systematically characterize regulatory elements, their cell-type specificities and their functional interactions. Focusing on cell-type-specific patterns of promoters and enhancers, we define multicell activity profiles for chromatin state, gene expression, regulatory motif enrichment and regulator expression. We use correlations between these profiles to link enhancers to putative target genes, and predict the cell-type-specific activators and repressors that modulate them. The resulting annotations and regulatory predictions have implications for the interpretation of genome-wide association studies. Top-scoring disease single nucleotide polymorphisms are frequently positioned within enhancer elements specifically active in relevant cell types, and in some cases affect a motif instance for a predicted regulator, thus suggesting a mechanism for the association. Our study presents a general framework for deciphering cis-regulatory connections and their roles in disease.

Despite rapid progress in characterizing transcription factor-driven reprogramming of somatic cells to an induced pluripotent stem cell (iPSC) state, many mechanistic questions still remain. To gain insight into the earliest events in the reprogramming process, we systematically analyzed the transcriptional and epigenetic changes that occur during early factor induction after discrete numbers of divisions. We observed rapid, genome-wide changes in the euchromatic histone modification, H3K4me2, at more than a thousand loci including large subsets of pluripotency-related or developmentally regulated gene promoters and enhancers. In contrast, patterns of the repressive H3K27me3 modification remained largely unchanged except for focused depletion specifically at positions where H3K4 methylation is gained. These chromatin regulatory events precede transcriptional changes within the corresponding loci. Our data provide evidence for an early, organized, and population-wide epigenetic response to ectopic reprogramming factors that clarify the temporal order through which somatic identity is reset during reprogramming.

In embryonic stem (ES) cells, bivalent chromatin domains with overlapping repressive (H3 lysine 27 tri-methylation) and activating (H3 lysine 4 tri-methylation) histone modifications mark the promoters of more than 2,000 genes. To gain insight into the structure and function of bivalent domains, we mapped key histone modifications and subunits of Polycomb-repressive complexes 1 and 2 (PRC1 and PRC2) genomewide in human and mouse ES cells by chromatin immunoprecipitation, followed by ultra high-throughput sequencing. We find that bivalent domains can be segregated into two classes -- the first occupied by both PRC2 and PRC1 (PRC1-positive) and the second specifically bound by PRC2 (PRC2-only). PRC1-positive bivalent domains appear functionally distinct as they more efficiently retain lysine 27 tri-methylation upon differentiation, show stringent conservation of chromatin state, and associate with an overwhelming number of developmental regulator gene promoters. We also used computational genomics to search for sequence determinants of Polycomb binding. This analysis revealed that the genomewide locations of PRC2 and PRC1 can be largely predicted from the locations, sizes, and underlying motif contents of CpG islands. We propose that large CpG islands depleted of activating motifs confer epigenetic memory by recruiting the full repertoire of Polycomb complexes in pluripotent cells.

Liquid biopsies come of age offering unexploited potential to monitor and react to tumor evolution. We developed a cost-effective assay to non-invasively determine the immune status of glioblastoma (GBM) patients. Employing newly developed printed peptide microarrays we assessed the B-cell response against tumor-associated antigens (TAAs) in 214 patients. Firstly, sera of long-term (36+ months, LTS, n=10) and short-term (6-10 months, STS, n=14) surviving patients were screened for prognostic antibodies against 1745 13-mer peptides covering known TAAs (TNC, EGFR, GLEA2, PHF3, FABP5, MAGEA3). Next, survival associations were investigated in two retrospective independent multicenter validation sets (n=61, n=129, all IDH1-wildtype). Reliability of measurements was tested using a second array technology (spotted arrays). LTS/STS screening analyses identified 106 differential antibody responses. Evaluating the Top30 peptides in validation set 1 revealed three prognostic peptides. Prediction of TNC peptide VCEDGFTGPDCAE was confirmed in a second set (p=0.043, HR=0.66 [0.44-0.99]) and was unrelated to TNC protein expression. Median signals of printed arrays correlated with pre-synthesized spotted microarrays (p<0.0002, R=0.33). Multiple survival analysis revealed independence of age, gender, KPI and MGMT status. We present a novel peptide microarray immune assay that identified increased anti-TNC VCEDGFTGPDCAE serum antibody titer as a promising non-invasive biomarker for prolonged survival.

Mammalian gene regulation is dependent on tissue-specific enhancers that can act across large distances to influence transcriptional activity. Mapping experiments have identified hundreds of thousands of putative enhancers whose functionality is supported by cell type-specific chromatin signatures and striking enrichments for disease-associated sequence variants. However, these studies did not address the in vivo functions of the putative elements or their chromatin states and did not determine which genes, if any, a given enhancer regulates. Here we present a strategy to investigate endogenous regulatory elements by selectively altering their chromatin state using programmable reagents. Transcription activator-like (TAL) effector repeat domains fused to the LSD1 histone demethylase efficiently remove enhancer-associated chromatin modifications from target loci, without affecting control regions. We find that inactivation of enhancer chromatin by these fusion proteins frequently causes downregulation of proximal genes, revealing enhancer target genes. Our study demonstrates the potential of epigenome editing tools to characterize an important class of functional genomic elements.

Mutation is a fundamental process in tumorigenesis. However, the degree to which the rate of somatic mutation varies across the human genome and the mechanistic basis underlying this variation remain to be fully elucidated. Here, we performed a cross-cancer comparison of 402 whole genomes comprising a diverse set of childhood and adult tumors, including both solid and hematopoietic malignancies. Surprisingly, we found that the inactive X chromosome of many female cancer genomes accumulates on average twice and up to four times as many somatic mutations per megabase, as compared to the individual autosomes. Whole-genome sequencing of clonally expanded hematopoietic stem/progenitor cells (HSPCs) from healthy individuals and a premalignant myelodysplastic syndrome (MDS) sample revealed no X chromosome hypermutation. Our data suggest that hypermutation of the inactive X chromosome is an early and frequent feature of tumorigenesis resulting from DNA replication stress in aberrantly proliferating cells.

Recurrent medulloblastoma is a therapeutic challenge because it is almost always fatal. Studies have confirmed that medulloblastoma consists of at least four distinct subgroups. We sought to delineate subgroup-specific differences in medulloblastoma recurrence patterns.

Telomerase reverse transcriptase (TERT) promoter mutations were recently shown to drive telomerase activity in various cancer types, including medulloblastoma. However, the clinical and biological implications of TERT mutations in medulloblastoma have not been described. Hence, we sought to describe these mutations and their impact in a subgroup-specific manner. We analyzed the TERT promoter by direct sequencing and genotyping in 466 medulloblastomas. The mutational distributions were determined according to subgroup affiliation, demographics, and clinical, prognostic, and molecular features. Integrated genomics approaches were used to identify specific somatic copy number alterations in TERT promoter-mutated and wild-type tumors. Overall, TERT promoter mutations were identified in 21 % of medulloblastomas. Strikingly, the highest frequencies of TERT mutations were observed in SHH (83 %; 55/66) and WNT (31 %; 4/13) medulloblastomas derived from adult patients. Group 3 and Group 4 harbored this alteration in <5 % of cases and showed no association with increased patient age. The prognostic implications of these mutations were highly subgroup-specific. TERT mutations identified a subset with good and poor prognosis in SHH and Group 4 tumors, respectively. Monosomy 6 was mostly restricted to WNT tumors without TERT mutations. Hallmark SHH focal copy number aberrations and chromosome 10q deletion were mutually exclusive with TERT mutations within SHH tumors. TERT promoter mutations are the most common recurrent somatic point mutation in medulloblastoma, and are very highly enriched in adult SHH and WNT tumors. TERT mutations define a subset of SHH medulloblastoma with distinct demographics, cytogenetics, and outcomes.

Two recurrent mutations, K27M and G34R/V, within histone variant H3.3 were recently identified in ∼50% of pHGGs. Both mutations define clinically and biologically distinct subgroups of pHGGs. Here, we provide further insight about the dominant-negative effect of K27M mutant H3.3, leading to a global reduction of the repressive histone mark H3K27me3. We demonstrate that this is caused by aberrant recruitment of the PRC2 complex to K27M mutant H3.3 and enzymatic inhibition of the H3K27me3-establishing methyltransferase EZH2. By performing chromatin immunoprecipitation followed by next-generation sequencing and whole-genome bisulfite sequencing in primary pHGGs, we show that reduced H3K27me3 levels and DNA hypomethylation act in concert to activate gene expression in K27M mutant pHGGs.

Pilocytic astrocytoma (PA) is the most frequent pediatric brain tumor. Activation of the MAPK pathway is well established as the oncogenic driver of the disease. It is most frequently caused by KIAA1549:BRAF fusions, and leads to oncogene induced senescence (OIS). OIS is thought to be a major reason for growth arrest of PA cells in vitro and in vivo, preventing establishment of PA cultures. Hence, valid preclinical models are currently very limited, but preclinical testing of new compounds is urgently needed. We transduced the PA short-term culture DKFZ-BT66 derived from the PA of a 2-year old patient with a doxycycline-inducible system coding for Simian Vacuolating Virus 40 Large T Antigen (SV40-TAg). SV40-TAg inhibits TP53/CDKN1A and CDKN2A/RB1, two pathways critical for OIS induction and maintenance. DNA methylation array and KIAA1549:BRAF fusion analysis confirmed pilocytic astrocytoma identity of DKFZ-BT66 cells after establishment. Readouts were analyzed in proliferating as well as senescent states, including cell counts, viability, cell cycle analysis, expression of SV40-Tag, CDKN2A (p16), CDKN1A (p21), and TP53 (p53) protein, and gene-expression profiling. Selected MAPK inhibitors (MAPKi) including clinically available MEK inhibitors (MEKi) were tested in vitro. Expression of SV40-TAg enabled the cells to bypass OIS and to resume proliferation with a mean doubling time of 45h allowing for propagation and long-term culture. Withdrawal of doxycycline led to an immediate decrease of SV40-TAg expression, appearance of senescent morphology, upregulation of CDKI proteins and a subsequent G1 growth arrest in line with the re-induction of senescence. DKFZ-BT66 cells still underwent replicative senescence that was overcome by TERT expression. Testing of a set of MAPKi revealed differential responses in DKFZ-BT66. MEKi efficiently inhibited MAPK signaling at clinically achievable concentrations, while BRAF V600E- and RAF Type II inhibitors showed paradoxical activation. Taken together, we have established the first patient-derived long term expandable PA cell line expressing the KIAA1549:BRAF-fusion suitable for preclinical drug testing.

SMARCB1 (also known as SNF5, INI1, and BAF47), a core subunit of the SWI/SNF (BAF) chromatin-remodeling complex, is inactivated in nearly all pediatric rhabdoid tumors. These aggressive cancers are among the most genomically stable, suggesting an epigenetic mechanism by which SMARCB1 loss drives transformation. Here we show that, despite having indistinguishable mutational landscapes, human rhabdoid tumors exhibit distinct enhancer H3K27ac signatures, which identify remnants of differentiation programs. We show that SMARCB1 is required for the integrity of SWI/SNF complexes and that its loss alters enhancer targeting-markedly impairing SWI/SNF binding to typical enhancers, particularly those required for differentiation, while maintaining SWI/SNF binding at super-enhancers. We show that these retained super-enhancers are essential for rhabdoid tumor survival, including some that are shared by all subtypes, such as SPRY1, and other lineage-specific super-enhancers, such as SOX2 in brain-derived rhabdoid tumors. Taken together, our findings identify a new chromatin-based epigenetic mechanism underlying the tumor-suppressive activity of SMARCB1.

Glioblastoma, the most common and aggressive malignant brain tumor, is propagated by stem-like cancer cells refractory to existing therapies. Understanding the molecular mechanisms that control glioblastoma stem cell (GSC) proliferation and drug resistance may reveal opportunities for therapeutic interventions. Here we show that GSCs can reversibly transition to a slow-cycling, persistent state in response to targeted kinase inhibitors. In this state, GSCs upregulate primitive developmental programs and are dependent upon Notch signaling. This transition is accompanied by widespread redistribution of repressive histone methylation. Accordingly, persister GSCs upregulate, and are dependent on, the histone demethylases KDM6A/B. Slow-cycling cells with high Notch activity and histone demethylase expression are present in primary glioblastomas before treatment, potentially contributing to relapse. Our findings illustrate how cancer cells may hijack aspects of native developmental programs for deranged proliferation, adaptation, and tolerance. They also suggest strategies for eliminating refractory tumor cells by targeting epigenetic and developmental pathways.

Based on extensive pre-clinical studies, the oncolytic parvovirus H-1 (H-1PV) is currently applied to patients with recurrent glioblastoma in a phase I/IIa clinical trial (ParvOryx01, NCT01301430). Cure rates of about 40% in pediatric high-risk medulloblastoma (MB) patients also indicate the need of new therapeutic approaches. In order to prepare a future application of oncolytic parvovirotherapy to MB, the present study preclinically evaluates the cytotoxic efficacy of H-1PV on MB cells in vitro and characterizes cellular target genes involved in this effect. Six MB cell lines were analyzed by whole genome oligonucleotide microarrays after treatment and the results were matched to known molecular and cytogenetic risk factors. In contrast to non-transformed infant astrocytes and neurons, in five out of six MB cell lines lytic H-1PV infection and efficient viral replication could be demonstrated. The cytotoxic effects induced by H-1PV were observed at LD50s below 0.05 p. f. u. per cell indicating high susceptibility. Gene expression patterns in the responsive MB cell lines allowed the identification of candidate target genes mediating the cytotoxic effects of H-1PV. H-1PV induced down-regulation of key regulators of early neurogenesis shown to confer poor prognosis in MB such as ZIC1, FOXG1B, MYC, and NFIA. In MB cell lines with genomic amplification of MYC, expression of MYC was the single gene most significantly repressed after H-1PV infection. H-1PV virotherapy may be a promising treatment approach for MB since it targets genes of functional relevance and induces cell death at very low titers of input virus.

Gene-editing nucleases enable targeted modification of DNA sequences in living cells, thereby facilitating efficient knockout and precise editing of endogenous loci. Engineered nucleases also have the potential to introduce mutations at off-target sites of action. Such unintended alterations can confound interpretation of experiments and can have implications for development of therapeutic applications. Recently, two improved methods for identifying the off-target effects of zinc finger nucleases (ZFNs) were described-one using an in vitro cleavage site selection method and the other exploiting the insertion of integration-defective lentiviruses into nuclease-induced double-stranded DNA breaks. However, application of these two methods to a ZFN pair targeted to the human CCR5 gene led to identification of largely non-overlapping off-target sites, raising the possibility that additional off-target sites might exist. Here, we show that in silico abstraction of ZFN cleavage profiles obtained from in vitro cleavage site selections can greatly enhance the ability to identify potential off-target sites in human cells. Our improved method should enable more comprehensive profiling of ZFN specificities.

Differentiation of human embryonic stem cells (hESCs) provides a unique opportunity to study the regulatory mechanisms that facilitate cellular transitions in a human context. To that end, we performed comprehensive transcriptional and epigenetic profiling of populations derived through directed differentiation of hESCs representing each of the three embryonic germ layers. Integration of whole-genome bisulfite sequencing, chromatin immunoprecipitation sequencing, and RNA sequencing reveals unique events associated with specification toward each lineage. Lineage-specific dynamic alterations in DNA methylation and H3K4me1 are evident at putative distal regulatory elements that are frequently bound by pluripotency factors in the undifferentiated hESCs. In addition, we identified germ-layer-specific H3K27me3 enrichment at sites exhibiting high DNA methylation in the undifferentiated state. A better understanding of these initial specification events will facilitate identification of deficiencies in current approaches, leading to more faithful differentiation strategies as well as providing insights into the rewiring of human regulatory programs during cellular transitions.

A remarkable observation emerging from recent cancer genome analyses is the identification of chromothripsis as a one-off genomic catastrophe, resulting in massive somatic DNA structural rearrangements (SRs). Largely due to lack of suitable model systems, the mechanistic basis of chromothripsis has remained elusive. We developed an integrative method termed "complex alterations after selection and transformation (CAST)," enabling efficient in vitro generation of complex DNA rearrangements including chromothripsis, using cell perturbations coupled with a strong selection barrier followed by massively parallel sequencing. We employed this methodology to characterize catastrophic SR formation processes, their temporal sequence, and their impact on gene expression and cell division. Our in vitro system uncovered a propensity of chromothripsis to occur in cells with damaged telomeres, and in particular in hyperploid cells. Analysis of primary medulloblastoma cancer genomes verified the link between hyperploidy and chromothripsis in vivo. CAST provides the foundation for mechanistic dissection of complex DNA rearrangement processes.

Posterior fossa ependymoma comprises two distinct molecular variants termed EPN_PFA and EPN_PFB that have a distinct biology and natural history. The therapeutic value of cytoreductive surgery and radiation therapy for posterior fossa ependymoma after accounting for molecular subgroup is not known.

DDX3X encodes a DEAD-box family RNA helicase (DDX3) commonly mutated in medulloblastoma, a highly aggressive cerebellar tumor affecting both children and adults. Despite being implicated in several facets of RNA metabolism, the nature and scope of DDX3's interactions with RNA remain unclear. Here, we show DDX3 collaborates extensively with the translation initiation machinery through direct binding to 5'UTRs of nearly all coding RNAs, specific sites on the 18S rRNA, and multiple components of the translation initiation complex. Impairment of translation initiation is also evident in primary medulloblastomas harboring mutations in DDX3X, further highlighting DDX3's role in this process. Arsenite-induced stress shifts DDX3 binding from the 5'UTR into the coding region of mRNAs concomitant with a general reduction of translation, and both the shift of DDX3 on mRNA and decreased translation are blunted by expression of a catalytically-impaired, medulloblastoma-associated DDX3R534H variant. Furthermore, despite the global repression of translation induced by arsenite, translation is preserved on select genes involved in chromatin organization in DDX3R534H-expressing cells. Thus, DDX3 interacts extensively with RNA and ribosomal machinery to help remodel the translation landscape in response to stress, while cancer-related DDX3 variants adapt this response to selectively preserve translation.