An integrated ChIP-seq analysis platform with customizable workflows

ChEAP: ChIP-exo analysis pipeline and the investigation of Escherichia coli RpoN protein-DNA interactions

Genome-scale studies of the bacterial regulatory network have been leveraged by declining sequencing cost and advances in ChIP (chromatin immunoprecipitation) methods. Of which, ChIP-exo has proven competent with its near-single base-pair resolution. While several algorithms and programs have been developed for different analytical steps in ChIP-exo data processing, there is a lack of effort in incorporating them into a convenient bioinformatics pipeline that is intuitive and publicly available. In this paper, we developed ChIP-exo Analysis Pipeline (ChEAP) that executes the one-step process, starting from trimming and aligning raw sequencing reads to visualization of ChIP-exo results. The pipeline was implemented on the interactive web-based Python development environment - Jupyter Notebook, which is compatible with the Google Colab cloud platform to facilitate the sharing of codes and collaboration among researchers. Additionally, users could exploit the free GPU and CPU resources allocated by Colab to carry out computing tasks regardless of the performance of their local machines. The utility of ChEAP was demonstrated with the ChIP-exo datasets of RpoN sigma factor in E. coli K-12 MG1655. To analyze two raw data files, ChEAP runtime was 2 min and 25 s. Subsequent analyses identified 113 RpoN binding sites showing a conserved RpoN binding pattern in the motif search. ChEAP application in ChIP-exo data analysis is extensive and flexible for the parallel processing of data from various organisms.

Super-enhancer hypermutation alters oncogene expression in B cell lymphoma

Diffuse large B cell lymphoma (DLBCL) is the most common B cell non-Hodgkin lymphoma and remains incurable in around 40% of patients. Efforts to sequence the coding genome identified several genes and pathways that are altered in this disease, including potential therapeutic targets1-5. However, the non-coding genome of DLBCL remains largely unexplored. Here we show that active super-enhancers are highly and specifically hypermutated in 92% of samples from individuals with DLBCL, display signatures of activation-induced cytidine deaminase activity, and are linked to genes that encode B cell developmental regulators and oncogenes. As evidence of oncogenic relevance, we show that the hypermutated super-enhancers linked to the BCL6, BCL2 and CXCR4 proto-oncogenes prevent the binding and transcriptional downregulation of the corresponding target gene by transcriptional repressors, including BLIMP1 (targeting BCL6) and the steroid receptor NR3C1 (targeting BCL2 and CXCR4). Genetic correction of selected mutations restored repressor DNA binding, downregulated target gene expression and led to the counter-selection of cells containing corrected alleles, indicating an oncogenic dependency on the super-enhancer mutations. This pervasive super-enhancer mutational mechanism reveals a major set of genetic lesions deregulating gene expression, which expands the involvement of known oncogenes in DLBCL pathogenesis and identifies new deregulated gene targets of therapeutic relevance.

Scaffold association factor B (SAFB) is required for expression of prenyltransferases and RAS membrane association

Inhibiting membrane association of RAS has long been considered a rational approach to anticancer therapy, which led to the development of farnesyltransferase inhibitors (FTIs). However, FTIs proved ineffective against KRAS-driven tumors. To reveal alternative therapeutic strategies, we carried out a genome-wide CRISPR-Cas9 screen designed to identify genes required for KRAS4B membrane association. We identified five enzymes in the prenylation pathway and SAFB, a nuclear protein with both DNA and RNA binding domains. Silencing SAFB led to marked mislocalization of all RAS isoforms as well as RAP1A but not RAB7A, a pattern that phenocopied silencing FNTA, the prenyltransferase α subunit shared by farnesyltransferase and geranylgeranyltransferase type I. We found that SAFB promoted RAS membrane association by controlling FNTA expression. SAFB knockdown decreased GTP loading of RAS, abrogated alternative prenylation, and sensitized RAS-mutant cells to growth inhibition by FTI. Our work establishes the prenylation pathway as paramount in KRAS membrane association, reveals a regulator of prenyltransferase expression, and suggests that reduction in FNTA expression may enhance the efficacy of FTIs.

Unique and Shared Epigenetic Programs of the CREBBP and EP300 Acetyltransferases in Germinal Center B Cells Reveal Targetable Dependencies in Lymphoma

Inactivating mutations of the CREBBP and EP300 acetyltransferases are among the most common genetic alterations in diffuse large B cell lymphoma (DLBCL) and follicular lymphoma (FL). Here, we examined the relationship between these two enzymes in germinal center (GC) B cells, the normal counterpart of FL and DLBCL, and in lymphomagenesis by using conditional GC-directed deletion mouse models targeting Crebbp or Ep300. We found that CREBBP and EP300 modulate common as well as distinct transcriptional programs implicated in separate anatomic and functional GC compartments. Consistently, deletion of Ep300 but not Crebbp impaired the fitness of GC B cells in vivo. Combined loss of Crebbp and Ep300 completely abrogated GC formation, suggesting that these proteins partially compensate for each other through common transcriptional targets. This synthetic lethal interaction was retained in CREBBP-mutant DLBCL cells and could be pharmacologically targeted with selective small molecule inhibitors of CREBBP and EP300 function. These data provide proof-of-principle for the clinical development of EP300-specific inhibitors in FL and DLBCL.

Crunch: integrated processing and modeling of ChIP-seq data in terms of regulatory motifs

Although ChIP-seq has become a routine experimental approach for quantitatively characterizing the genome-wide binding of transcription factors (TFs), computational analysis procedures remain far from standardized, making it difficult to compare ChIP-seq results across experiments. In addition, although genome-wide binding patterns must ultimately be determined by local constellations of DNA-binding sites, current analysis is typically limited to identifying enriched motifs in ChIP-seq peaks. Here we present Crunch, a completely automated computational method that performs all ChIP-seq analysis from quality control through read mapping and peak detecting and that integrates comprehensive modeling of the ChIP signal in terms of known and novel binding motifs, quantifying the contribution of each motif and annotating which combinations of motifs explain each binding peak. By applying Crunch to 128 data sets from the ENCODE Project, we show that Crunch outperforms current peak finders and find that TFs naturally separate into "solitary TFs," for which a single motif explains the ChIP-peaks, and "cobinding TFs," for which multiple motifs co-occur within peaks. Moreover, for most data sets, the motifs that Crunch identified de novo outperform known motifs, and both the set of cobinding motifs and the top motif of solitary TFs are consistent across experiments and cell lines. Crunch is implemented as a web server, enabling standardized analysis of any collection of ChIP-seq data sets by simply uploading raw sequencing data. Results are provided both in a graphical web interface and as downloadable files.

Rational Targeting of Cooperating Layers of the Epigenome Yields Enhanced Therapeutic Efficacy against AML

Disruption of epigenetic regulation is a hallmark of acute myeloid leukemia (AML), but epigenetic therapy is complicated by the complexity of the epigenome. Herein, we developed a long-term primary AML ex vivo platform to determine whether targeting different epigenetic layers with 5-azacytidine and LSD1 inhibitors would yield improved efficacy. This combination was most effective in TET2 ^mut AML, where it extinguished leukemia stem cells and particularly induced genes with both LSD1-bound enhancers and cytosine-methylated promoters. Functional studies indicated that derepression of genes such as GATA2 contributes to drug efficacy. Mechanistically, combination therapy increased enhancer-promoter looping and chromatin-activating marks at the GATA2 locus. CRISPRi of the LSD1-bound enhancer in patient-derived TET2 ^mut AML was associated with dampening of therapeutic GATA2 induction. TET2 knockdown in human hematopoietic stem/progenitor cells induced loss of enhancer 5-hydroxymethylation and facilitated LSD1-mediated enhancer inactivation. Our data provide a basis for rational targeting of cooperating aberrant promoter and enhancer epigenetic marks driven by mutant epigenetic modifiers. SIGNIFICANCE: Somatic mutations of genes encoding epigenetic modifiers are a hallmark of AML and potentially disrupt many components of the epigenome. Our study targets two different epigenetic layers at promoters and enhancers that cooperate to aberrant gene silencing, downstream of the actions of a mutant epigenetic regulator.This article is highlighted in the In This Issue feature, p. 813.

Drug-Induced Expression-Based Computational Repurposing of Small Molecules Affecting Transcription Factor Activity

Inhibition of oncogenes and reactivation of tumor suppressors are well-established goals in anticancer drug development. Unfortunately many oncogenes and tumor suppressors are not classically druggable, in that they lack a targetable enzymatic activity and associated binding pockets that small molecule drugs can be directed to. This is especially relevant for transcription factors, which have long been thought to be undruggable. To address this gap, we have developed and described CRAFTT, a broadly applicable computational drug-repositioning approach for targeting transcription factors. CRAFTT combines transcription factor target gene sets with drug-induced expression profiling to identify small molecules that can perturb transcription factor activity. Network analysis is then used to derive a modulation index (MI) and prioritize predictions.

MEF2B Instructs Germinal Center Development and Acts as an Oncogene in B Cell Lymphomagenesis

The gene encoding the MEF2B transcription factor is mutated in germinal center (GC)-derived B cell lymphomas, but its role in GC development and lymphomagenesis is unknown. We demonstrate that Mef2b deletion reduces GC formation in mice and identify MEF2B transcriptional targets in GC, with roles in cell proliferation, apoptosis, GC confinement, and differentiation. The most common lymphoma-associated MEF2B mutant (MEF2B^D83V) is hypomorphic, yet escapes binding and negative regulation by components of the HUCA complex and class IIa HDACs. Mef2b^D83V expression in mice leads to GC enlargement and lymphoma development, a phenotype that becomes fully penetrant in combination with BCL2 de-regulation, an event associated with human MEF2B mutations. These results identify MEF2B as a critical GC regulator and a driver oncogene in lymphomagenesis.

IFN-γ Induces Histone 3 Lysine 27 Trimethylation in a Small Subset of Promoters to Stably Silence Gene Expression in Human Macrophages

The mechanisms by which IFN-γ activates expression of interferon-stimulated genes that have inflammatory and host defense functions are well understood. In contrast, little is known about how IFN-γ represses gene expression. By using transcriptomic and epigenomic analysis, we found that stable repression of a small group of genes by IFN-γ is associated with recruitment of the histone methyltransferase EZH2 and deposition of the negative mark histone 3 lysine 27 trimethylation (H3K27me3) at their promoters. Repressed genes included MERTK, PPARG, and RANK, which have anti-inflammatory functions and promote osteoclast differentiation. Gene repression and H3K27me3 persisted after IFN-γ signaling was terminated, and these silenced genes were no longer responsive to glucocorticoids, IL-4, and M-CSF. These results identify cytokine-induced H3K27 trimethylation as a mechanism that stabilizes gene silencing in macrophages. IFN-γ-induced macrophage activation is thus reinforced by a chromatin-based mechanism that blocks anti-inflammatory and opposing pathways.

Integrative analysis of methylomic and transcriptomic data in fetal sheep muscle tissues in response to maternal diet during pregnancy

BACKGROUND

Numerous studies have established a link between maternal diet and the physiological and metabolic phenotypes of their offspring. In previous studies in sheep, we demonstrated that different maternal diets altered the transcriptome of fetal tissues. However, the mechanisms underlying transcriptomic changes are poorly understood. DNA methylation is an epigenetic mark regulating transcription and is largely influenced by dietary components of the one-carbon cycle that generate the methyl group donor, SAM. Therefore, in the present study, we tested whether different maternal diets during pregnancy would alter the DNA methylation and gene expression patterns in fetal tissues.

RESULTS

Pregnant ewes were randomly divided into two groups which received either hay or corn diet from mid-gestation (day 67 ± 5) until day 131 ± 1 when fetuses were collected by necropsy. A total of 1516 fetal longissimus dorsi (LD) tissues were used for DNA methylation analysis and gene expression profiling. Whole genome DNA methylation using methyl-binding domain enrichment analysis revealed 60 differentially methylated regions (DMRs) between hay and corn fetuses with 39 DMRs more highly methylated in the hay fetuses vs. 21 DMRs more highly methylated in the corn fetuses. Three DMRs (LPAR3, PLIN5-PLIN4, and the differential methylation of a novel lincRNA) were validated using bisulfite sequencing. These DMRs were associated with differential gene expression. Additionally, significant DNA methylation differences were found at the single CpG level. Integrative methylome and transcriptome analysis revealed an association between gene expression and inter-/intragenic methylated regions. Furthermore, intragenic DMRs were found to be associated with expression of neighboring genes.

CONCLUSIONS

The findings of this study imply that maternal diet from mid- to late-gestation can shape the epigenome and the transcriptome of fetal tissues, and putatively affect phenotypes of the lambs.

Analyzing DNA Methylation Patterns During Tumor Evolution

Epigenetic modifications play a key role in cellular development and tumorigenesis. Recent large-scale genomic studies have shown that mutations in players of the epigenetic machinery and concomitant perturbation of epigenomic patterning are frequent events in tumors. Among epigenetic marks, DNA methylation is one of the best studied. Hyper- and hypo-methylation events of specific regulatory elements (such as promoters and enhancers) are sometimes thought to be correlated with expression of nearby genes. High-throughput bisulfite converted sequencing is currently the technology of choice for studying DNA methylation in base-pair resolution and on whole-genome scale. Such broad and high-resolution coverage investigations of the epigenome provide unprecedented opportunities to analyze DNA methylation patterns, which are correlated with tumorigenesis, tumor evolution, and tumor progression. However, few computational pipelines are available to the public to perform systematic DNA methylation analysis. Utilizing open source tools, we here describe a comprehensive computational methodology to thoroughly analyze DNA methylation patterns during tumor evolution based on bisulfite converted sequencing data, including intra-tumor methylation heterogeneity.

AICDA drives epigenetic heterogeneity and accelerates germinal center-derived lymphomagenesis

Epigenetic heterogeneity is emerging as a feature of tumors. In diffuse large B-cell lymphoma (DLBCL), increased cytosine methylation heterogeneity is associated with poor clinical outcome, yet the underlying mechanisms remain unclear. Activation-induced cytidine deaminase (AICDA), an enzyme that mediates affinity maturation and facilitates DNA demethylation in germinal center (GC) B cells, is required for DLBCL pathogenesis and linked to inferior outcome. Here we show that AICDA overexpression causes more aggressive disease in BCL2-driven murine lymphomas. This phenotype is associated with increased cytosine methylation heterogeneity, but not with increased AICDA-mediated somatic mutation burden. Reciprocally, the cytosine methylation heterogeneity characteristic of normal GC B cells is lost upon AICDA depletion. These observations are relevant to human patients, since DLBCLs with high AICDA expression manifest increased methylation heterogeneity vs. AICDA-low DLBCLs. Our results identify AICDA as a driver of epigenetic heterogeneity in B-cell lymphomas with potential significance for other tumors with aberrant expression of cytidine deaminases.

In diffuse large B-cell lymphoma (DLBCL) increased epigenetic heterogeneity in the form of cytosine methylation is known to link to a poor clinical outcome. Here, the authors show that AICDA, an enzyme required for DLBCL pathogenesis, increases cytosine methylation heterogeneity.

MuSERA: Multiple Sample Enriched Region Assessment

Enriched region (ER) identification is a fundamental step in several next-generation sequencing (NGS) experiment types. Yet, although NGS experimental protocols recommend producing replicate samples for each evaluated condition and their consistency is usually assessed, typically pipelines for ER identification do not consider available NGS replicates. This may alter genome-wide descriptions of ERs, hinder significance of subsequent analyses on detected ERs and eventually preclude biological discoveries that evidence in replicate could support. MuSERA is a broadly useful stand-alone tool for both interactive and batch analysis of combined evidence from ERs in multiple ChIP-seq or DNase-seq replicates. Besides rigorously combining sample replicates to increase statistical significance of detected ERs, it also provides quantitative evaluations and graphical features to assess the biological relevance of each determined ER set within its genomic context; they include genomic annotation of determined ERs, nearest ER distance distribution, global correlation assessment of ERs and an integrated genome browser. We review MuSERA rationale and implementation, and illustrate how sets of significant ERs are expanded by applying MuSERA on replicates for several types of NGS data, including ChIP-seq of transcription factors or histone marks and DNase-seq hypersensitive sites. We show that MuSERA can determine a new, enhanced set of ERs for each sample by locally combining evidence on replicates, and prove how the easy-to-use interactive graphical displays and quantitative evaluations that MuSERA provides effectively support thorough inspection of obtained results and evaluation of their biological content, facilitating their understanding and biological interpretations. MuSERA is freely available at http://www.bioinformatics.deib.polimi.it/MuSERA/.

TET proteins safeguard bivalent promoters from de novo methylation in human embryonic stem cells

The TET enzymes oxidize 5-methylcytosine to 5-hydroxymethylcytosine, which can lead to DNA demethylation. However, direct connections between TET-mediated DNA demethylation and transcriptional output are difficult to establish due to challenges of distinguishing global versus locus-specific effects. Here we show that TET1/2/3 triple knockout (TKO) human embryonic stem cells (hESCs) exhibit prominent bivalent promoter hypermethylation without an overall corresponding gene expression decrease in the undifferentiated state. Focusing on the bivalent PAX6 locus, we find increased DNMT3B binding is associated with promoter hypermethylation, which precipitates a neural differentiation defect and failure of PAX6 induction during differentiation. dCas9-mediated locus-specific demethylation and global inactivation of DNMT3B in TKO hESCs partially reverses the hypermethylation at the PAX6 promoter and improves differentiation to neuroectoderm. Taken together with further genome-wide methylation and TET1 and DNMT3B ChIP-Seq analysis, we conclude that the TET proteins safeguard bivalent promoters from de novo methylation to ensure robust lineage-specific transcription upon differentiation.

High-intensity UV laser ChIP-seq for the study of protein-DNA interactions in living cells

Genome-wide mapping of transcription factor binding is generally performed by chemical protein–DNA crosslinking, followed by chromatin immunoprecipitation and deep sequencing (ChIP-seq). Here we present the ChIP-seq technique based on photochemical crosslinking of protein–DNA interactions by high-intensity ultraviolet (UV) laser irradiation in living mammalian cells (UV-ChIP-seq). UV laser irradiation induces an efficient and instant formation of covalent “zero-length” crosslinks exclusively between nucleic acids and proteins that are in immediate contact, thus resulting in a “snapshot” of direct protein–DNA interactions in their natural environment. Here we show that UV-ChIP-seq, applied for genome-wide profiling of the sequence-specific transcriptional repressor B-cell lymphoma 6 (BCL6) in human diffuse large B-cell lymphoma (DLBCL) cells, produces sensitive and precise protein–DNA binding profiles, highly enriched with canonical BCL6 DNA sequence motifs. Using this technique, we also found numerous previously undetectable direct BCL6 binding sites, particularly in condensed, inaccessible areas of chromatin.

Chromatin immunoprecipitation followed by DNA sequencing (ChIP-seq) can map transcription factor binding to a given genome. Here, Steube and colleagues devise a ChIP-seq technique based on photochemical crosslinking by high-intensity ultraviolet laser irradiation, and describe its utility.

ORIO (Online Resource for Integrative Omics): a web-based platform for rapid integration of next generation sequencing data

Established and emerging next generation sequencing (NGS)-based technologies allow for genome-wide interrogation of diverse biological processes. However, accessibility of NGS data remains a problem, and few user-friendly resources exist for integrative analysis of NGS data from different sources and experimental techniques. Here, we present Online Resource for Integrative Omics (ORIO; https://orio.niehs.nih.gov/), a web-based resource with an intuitive user interface for rapid analysis and integration of NGS data. To use ORIO, the user specifies NGS data of interest along with a list of genomic coordinates. Genomic coordinates may be biologically relevant features from a variety of sources, such as ChIP-seq peaks for a given protein or transcription start sites from known gene models. ORIO first iteratively finds read coverage values at each genomic feature for each NGS dataset. Data are then integrated using clustering-based approaches, giving hierarchical relationships across NGS datasets and separating individual genomic features into groups. In focusing its analysis on read coverage, ORIO makes limited assumptions about the analyzed data; this allows the tool to be applied across data from a variety of experiments and techniques. Results from analysis are presented in dynamic displays alongside user-controlled statistical tests, supporting rapid statistical validation of observed results. We emphasize the versatility of ORIO through diverse examples, ranging from NGS data quality control to characterization of enhancer regions and integration of gene expression information. Easily accessible on a public web server, we anticipate wide use of ORIO in genome-wide investigations by life scientists.

Genome-wide analysis of DNA methylation in hypothalamus and ovary of Capra hircus

BACKGROUND

DNA methylation is a frequently studied epigenetic modification due to its role in regulating gene expression and hence in biological processes and in determining phenotypic plasticity in organisms. Rudimentary DNA methylation patterns for some livestock species are publically available: among these, goat methylome deserves to be further explored.

RESULTS

Genome-wide DNA methylation maps of the hypothalamus and ovary from Saanen goats were generated using Methyl-CpG binding domain protein sequencing (MBD-seq). Analysis of DNA methylation patterns indicate that the majority of methylation peaks found within genes are located gene body regions, for both organs. Analysis of the distribution of methylated sites per chromosome showed that chromosome X had the lowest number of methylation peaks. The X chromosome has one of the highest percentages of methylated CpG islands in both organs, and approximately 50% of the CpG islands in the goat epigenome are methylated in hypothalamus and ovary. Organ-specific Differentially Methylated Genes (DMGs) were correlated with the expression levels.

CONCLUSIONS

The comparison between transcriptome and methylome in hypothalamus and ovary showed that a higher level of methylation is not accompanied by a higher gene suppression. The genome-wide DNA methylation map for two goat organs produced here is a valuable starting point for studying the involvement of epigenetic modifications in regulating goat reproduction performance.

Epigenetic Identity in AML Depends on Disruption of Nonpromoter Regulatory Elements and Is Affected by Antagonistic Effects of Mutations in Epigenetic Modifiers

We performed cytosine methylation sequencing on genetically diverse patients with acute myeloid leukemia (AML) and found leukemic DNA methylation patterning is primarily driven by nonpromoter regulatory elements and CpG shores. Enhancers displayed stronger differential methylation than promoters, consisting predominantly of hypomethylation. AMLs with dominant hypermethylation featured greater epigenetic disruption of promoters, whereas those with dominant hypomethylation displayed greater disruption of distal and intronic regions. Mutations in IDH and DNMT3A had opposing and mutually exclusive effects on the epigenome. Notably, co-occurrence of both mutations resulted in epigenetic antagonism, with most CpGs affected by either mutation alone no longer affected in double-mutant AMLs. Importantly, this epigenetic antagonism precedes malignant transformation and can be observed in preleukemic LSK cells from Idh2^R140Q or Dnmt3a^R882H single-mutant and Idh2^R140Q/Dnmt3a^R882H double-mutant mice. Notably, IDH/DNMT3A double-mutant AMLs manifested upregulation of a RAS signaling signature and displayed unique sensitivity to MEK inhibition ex vivo as compared with AMLs with either single mutation.Significance: AML is biologically heterogeneous with subtypes characterized by specific genetic and epigenetic abnormalities. Comprehensive DNA methylation profiling revealed that differential methylation of nonpromoter regulatory elements is a driver of epigenetic identity, that gene mutations can be context-dependent, and that co-occurrence of mutations in epigenetic modifiers can result in epigenetic antagonism. Cancer Discov; 7(8); 868-83. ©2017 AACR.This article is highlighted in the In This Issue feature, p. 783.

Common nonmutational NOTCH1 activation in chronic lymphocytic leukemia

Activating mutations of NOTCH1 (a well-known oncogene in T-cell acute lymphoblastic leukemia) are present in ∼4-13% of chronic lymphocytic leukemia (CLL) cases, where they are associated with disease progression and chemorefractoriness. However, the specific role of NOTCH1 in leukemogenesis remains to be established. Here, we report that the active intracellular portion of NOTCH1 (ICN1) is detectable in ∼50% of peripheral blood CLL cases lacking gene mutations. We identify a "NOTCH1 gene-expression signature" in CLL cells, and show that this signature is significantly enriched in primary CLL cases expressing ICN1, independent of NOTCH1 mutation. NOTCH1 target genes include key regulators of B-cell proliferation, survival, and signal transduction. In particular, we show that NOTCH1 transactivates MYC via binding to B-cell-specific regulatory elements, thus implicating this oncogene in CLL development. These results significantly extend the role of NOTCH1 in CLL pathogenesis, and have direct implications for specific therapeutic targeting.

The CREBBP Acetyltransferase Is a Haploinsufficient Tumor Suppressor in B-cell Lymphoma

Inactivating mutations of the CREBBP acetyltransferase are highly frequent in diffuse large B-cell lymphoma (DLBCL) and follicular lymphoma (FL), the two most common germinal center (GC)-derived cancers. However, the role of CREBBP inactivation in lymphomagenesis remains unclear. Here, we show that CREBBP regulates enhancer/super-enhancer networks with central roles in GC/post-GC cell fate decisions, including genes involved in signal transduction by the B-cell receptor and CD40 receptor, transcriptional control of GC and plasma cell development, and antigen presentation. Consistently, Crebbp-deficient B cells exhibit enhanced response to mitogenic stimuli and perturbed plasma cell differentiation. Although GC-specific loss of Crebbp was insufficient to initiate malignant transformation, compound Crebbp-haploinsufficient/BCL2-transgenic mice, mimicking the genetics of FL and DLBCL, develop clonal lymphomas recapitulating the features of the human diseases. These findings establish CREBBP as a haploinsufficient tumor-suppressor gene in GC B cells and provide insights into the mechanisms by which its loss contributes to lymphomagenesis.Significance: Loss-of-function mutations of CREBBP are common and early lesions in FL and DLBCL, suggesting a prominent role in lymphoma initiation. Our studies identify the cellular program by which reduced CREBBP dosage facilitates malignant transformation, and have direct implications for targeted lymphoma therapy based on drugs affecting CREBBP-mediated chromatin acetylation. Cancer Discov; 7(3); 322-37. ©2017 AACR.This article is highlighted in the In This Issue feature, p. 235.

CREBBP Inactivation Promotes the Development of HDAC3-Dependent Lymphomas

Somatic mutations in CREBBP occur frequently in B-cell lymphoma. Here, we show that loss of CREBBP facilitates the development of germinal center (GC)-derived lymphomas in mice. In both human and murine lymphomas, CREBBP loss-of-function resulted in focal depletion of enhancer H3K27 acetylation and aberrant transcriptional silencing of genes that regulate B-cell signaling and immune responses, including class II MHC. Mechanistically, CREBBP-regulated enhancers are counter-regulated by the BCL6 transcriptional repressor in a complex with SMRT and HDAC3, which we found to bind extensively to MHC class II loci. HDAC3 loss-of-function rescued repression of these enhancers and corresponding genes, including MHC class II, and more profoundly suppressed CREBBP-mutant lymphomas in vitro and in vivo Hence, CREBBP loss-of-function contributes to lymphomagenesis by enabling unopposed suppression of enhancers by BCL6/SMRT/HDAC3 complexes, suggesting HDAC3-targeted therapy as a precision approach for CREBBP-mutant lymphomas.

SIGNIFICANCE

Our findings establish the tumor suppressor function of CREBBP in GC lymphomas in which CREBBP mutations disable acetylation and result in unopposed deacetylation by BCL6/SMRT/HDAC3 complexes at enhancers of B-cell signaling and immune response genes. Hence, inhibition of HDAC3 can restore the enhancer histone acetylation and may serve as a targeted therapy for CREBBP-mutant lymphomas. Cancer Discov; 7(1); 38-53. ©2016 AACR.See related commentary by Höpken, p. 14This article is highlighted in the In This Issue feature, p. 1.

Systematic Discovery of Chromatin-Bound Protein Complexes from ChIP-seq Datasets

Chromatin immunoprecipitation followed by sequencing is an invaluable assay for identifying the genomic binding sites of transcription factors. However, transcription factors rarely bind chromatin alone but often bind together with other cofactors, forming protein complexes. Here, we describe a computational method that integrates multiple ChIP-seq and RNA-seq datasets to discover protein complexes and determine their role as activators or repressors. This chapter outlines a detailed computational pipeline for discovering and predicting binding partners from ChIP-seq data and inferring their role in regulating gene expression. This work aims at developing hypotheses about gene regulation via binding partners and deciphering the combinatorial nature of DNA-binding proteins.

ChiLin: a comprehensive ChIP-seq and DNase-seq quality control and analysis pipeline

BACKGROUND

Transcription factor binding, histone modification, and chromatin accessibility studies are important approaches to understanding the biology of gene regulation. ChIP-seq and DNase-seq have become the standard techniques for studying protein-DNA interactions and chromatin accessibility respectively, and comprehensive quality control (QC) and analysis tools are critical to extracting the most value from these assay types. Although many analysis and QC tools have been reported, few combine ChIP-seq and DNase-seq data analysis and quality control in a unified framework with a comprehensive and unbiased reference of data quality metrics.

RESULTS

ChiLin is a computational pipeline that automates the quality control and data analyses of ChIP-seq and DNase-seq data. It is developed using a flexible and modular software framework that can be easily extended and modified. ChiLin is ideal for batch processing of many datasets and is well suited for large collaborative projects involving ChIP-seq and DNase-seq from different designs. ChiLin generates comprehensive quality control reports that include comparisons with historical data derived from over 23,677 public ChIP-seq and DNase-seq samples (11,265 datasets) from eight literature-based classified categories. To the best of our knowledge, this atlas represents the most comprehensive ChIP-seq and DNase-seq related quality metric resource currently available. These historical metrics provide useful heuristic quality references for experiment across all commonly used assay types. Using representative datasets, we demonstrate the versatility of the pipeline by applying it to different assay types of ChIP-seq data. The pipeline software is available open source at https://github.com/cfce/chilin .

CONCLUSION

ChiLin is a scalable and powerful tool to process large batches of ChIP-seq and DNase-seq datasets. The analysis output and quality metrics have been structured into user-friendly directories and reports. We have successfully compiled 23,677 profiles into a comprehensive quality atlas with fine classification for users.

EpiMINE, a computational program for mining epigenomic data

Background

In epigenetic research, both the increasing ease of high-throughput sequencing and a greater interest in genome-wide studies have resulted in an exponential flooding of epigenetic-related data in public domain. This creates an opportunity for exploring data outside the limits of any specific query-centred study. Such data have to undergo standard primary analyses that are accessible with multiple well-stabilized programs. Further downstream analyses, such as genome-wide comparative, correlative and quantitative analyses, are critical in deciphering key biological features. However, these analyses are only accessible for computational researchers and completely lack platforms capable of handling, analysing and linking multiple interdisciplinary datasets with efficient analytical methods.

Results

Here, we present EpiMINE, a program for mining epigenomic data. It is a user-friendly, stand-alone computational program designed to support multiple datasets, for performing genome-wide correlative and quantitative analysis of ChIP-seq and RNA-seq data. Using data available from the ENCODE project, we illustrated several features of EpiMINE through different biological scenarios to show how easy some known observations can be verified. These results highlight how these approaches can be helpful in identifying novel biological features.

Conclusions

EpiMINE performs different kinds of genome-wide quantitative and correlative analyses, using ChIP-seq- and RNA-seq-related datasets. Its framework enables it to be used by both experimental and computational researchers. EpiMINE can be downloaded from https://sourceforge.net/projects/epimine/.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0095-z) contains supplementary material, which is available to authorized users.

Multi-tiered Reorganization of the Genome during B Cell Affinity Maturation Anchored by a Germinal Center-Specific Locus Control Region

During the humoral immune response, B cells undergo a dramatic change in phenotype to enable antibody affinity maturation in germinal centers (GCs). Using genome-wide chromosomal conformation capture (Hi-C), we found that GC B cells undergo massive reorganization of the genomic architecture that encodes the GC B cell transcriptome. Coordinate expression of genes that specify the GC B cell phenotype-most prominently BCL6-was achieved through a multilayered chromatin reorganization process involving (1) increased promoter connectivity, (2) formation of enhancer networks, (3) 5' to 3' gene looping, and (4) merging of gene neighborhoods that share active epigenetic marks. BCL6 was an anchor point for the formation of GC-specific gene and enhancer loops on chromosome 3. Deletion of a GC-specific, highly interactive locus control region upstream of Bcl6 abrogated GC formation in mice. Thus, large-scale and multi-tiered genomic three-dimensional reorganization is required for coordinate expression of phenotype-driving gene sets that determine the unique characteristics of GC B cells.

fluff: exploratory analysis and visualization of high-throughput sequencing data

Summary. In this article we describe fluff, a software package that allows for simple exploration, clustering and visualization of high-throughput sequencing data mapped to a reference genome. The package contains three command-line tools to generate publication-quality figures in an uncomplicated manner using sensible defaults. Genome-wide data can be aggregated, clustered and visualized in a heatmap, according to different clustering methods. This includes a predefined setting to identify dynamic clusters between different conditions or developmental stages. Alternatively, clustered data can be visualized in a bandplot. Finally, fluff includes a tool to generate genomic profiles. As command-line tools, the fluff programs can easily be integrated into standard analysis pipelines. The installation is straightforward and documentation is available at http://fluff.readthedocs.org.

Availability. fluff is implemented in Python and runs on Linux. The source code is freely available for download at https://github.com/simonvh/fluff.

Regulation of the BMP Signaling-Responsive Transcriptional Network in the Drosophila Embryo

The BMP signaling pathway has a conserved role in dorsal-ventral axis patterning during embryonic development. In Drosophila, graded BMP signaling is transduced by the Mad transcription factor and opposed by the Brinker repressor. In this study, using the Drosophila embryo as a model, we combine RNA-seq with Mad and Brinker ChIP-seq to decipher the BMP-responsive transcriptional network underpinning differentiation of the dorsal ectoderm during dorsal-ventral axis patterning. We identify multiple new BMP target genes, including positive and negative regulators of EGF signaling. Manipulation of EGF signaling levels by loss- and gain-of-function studies reveals that EGF signaling negatively regulates embryonic BMP-responsive transcription. Therefore, the BMP gene network has a self-regulating property in that it establishes a balance between its activity and that of the antagonistic EGF signaling pathway to facilitate correct patterning. In terms of BMP-dependent transcription, we identify key roles for the Zelda and Zerknüllt transcription factors in establishing the resulting expression domain, and find widespread binding of insulator proteins to the Mad and Brinker-bound genomic regions. Analysis of embryos lacking the BEAF-32 insulator protein shows reduced transcription of a peak BMP target gene and a reduction in the number of amnioserosa cells, the fate specified by peak BMP signaling. We incorporate our findings into a model for Mad-dependent activation, and discuss its relevance to BMP signal interpretation in vertebrates.

Embryogenesis involves the patterning of many different cell fates by a limited number of types of signals. One way that these signals promote a particular cell fate is through the induction of a complex, yet highly reproducible, gene expression programme that instructs changes in the cell. For example, there is a conserved role for BMP signals in specifying cell fates during dorsal-ventral axis patterning. Here, we have used genomics approaches to identify the gene expression programme implemented in response to BMP signaling during axis patterning in the Drosophila embryo. Part of the gene network downstream of BMP signaling includes members of the EGF signaling pathway, with our data highlighting reciprocal interactions between these two pathways. We have also determined genome-wide binding of BMP-responsive transcription factors to gain new insights into how the BMP gene network is activated. Our data reveal roles for specific transcription factors and insulator binding proteins, with the latter traditionally associated with the separation of transcriptional domains. Overall, our data will provide a platform for exploiting the tractability of the Drosophila embryo to determine which features of the network are critical drivers of BMP-induced cell fate changes during embryogenesis.

The transcriptional repressor Hes1 attenuates inflammation by regulating transcription elongation

Most of the known regulatory mechanisms that curb inflammatory gene expression target pre-transcription-initiation steps, and evidence for post-initiation regulation of inflammatory gene expression remains scarce. We found that the transcriptional repressor Hes1 suppressed production of CXCL1, a chemokine that is crucial for recruiting neutrophils. Hes1 negatively regulated neutrophil recruitment in vivo in a manner that was dependent on macrophage-produced CXCL1, and it attenuated the severity of inflammatory arthritis. Mechanistically, inhibition of Cxcl1 expression by Hes1 did not involve modification of transcription initiation. Instead, Hes1 inhibited signal-induced recruitment of the positive transcription-elongation complex P-TEFb and thereby prevented phosphorylation of RNA polymerase II at Ser2 and productive elongation. Thus, our results identify Hes1 as a homeostatic suppressor of inflammatory responses that exerts its suppressive function by regulating transcription elongation.

Muscleblind-like 1 suppresses breast cancer metastatic colonization and stabilizes metastasis suppressor transcripts

Post-transcriptional deregulation is a defining feature of metastatic cancer. While many microRNAs have been implicated as regulators of metastatic progression, less is known about the roles and mechanisms of RNA-binding proteins in this process. We identified muscleblind-like 1 (MBNL1), a gene implicated in myotonic dystrophy, as a robust suppressor of multiorgan breast cancer metastasis. MBNL1 binds the 3' untranslated regions (UTRs) of DBNL (drebrin-like protein) and TACC1 (transforming acidic coiled-coil containing protein 1)-two genes that we implicate as metastasis suppressors. By enhancing the stability of these genes' transcripts, MBNL1 suppresses cell invasiveness. Consistent with these findings, elevated MBNL1 expression in human breast tumors is associated with reduced metastatic relapse likelihood. Our findings delineate a post-transcriptional network that governs breast cancer metastasis through RNA-binding protein-mediated transcript stabilization.

Epigenetic profiles signify cell fate plasticity in unipotent spermatogonial stem and progenitor cells

Spermatogonial stem and progenitor cells (SSCs) generate adult male gametes. During in vitro expansion, these unipotent murine cells spontaneously convert to multipotent adult spermatogonial-derived stem cells (MASCs). Here we investigate this conversion process through integrative transcriptomic and epigenomic analyses. We find in SSCs that promoters essential to maintenance and differentiation of embryonic stem cells (ESCs) are enriched with histone H3-lysine4 and -lysine 27 trimethylations. These bivalent modifications are maintained at most somatic promoters after conversion, bestowing MASCs an ESC-like promoter chromatin. At enhancers, the core pluripotency circuitry is activated partially in SSCs and completely in MASCs, concomitant with loss of germ cell-specific gene expression and initiation of embryonic-like programs. Furthermore, SSCs in vitro maintain the epigenomic characteristics of germ cells in vivo. Our observations suggest that SSCs encode innate plasticity through the epigenome and that both conversion of promoter chromatin states and activation of cell type-specific enhancers are prominent features of reprogramming.

Spermatogonial stem cells (SSCs) spontaneously convert to multipotent adult spermatogonial-derived stem cells (MASCs). Here, the authors reveal the dynamics of bivalent histone H3-lysine4 and -lysine27 methylation signatures at somatic gene promoters in SSCs and ESC-like promoter chromatin states in MASCs.

Disruption of KMT2D perturbs germinal center B cell development and promotes lymphomagenesis

Mutations in the gene encoding the KMT2D (or MLL2) methyltransferase are highly recurrent and occur early during tumorigenesis in diffuse large B cell lymphoma (DLBCL) and follicular lymphoma (FL). However, the functional consequences of these mutations and their role in lymphomagenesis are unknown. Here we show that FL- and DLBCL-associated KMT2D mutations impair KMT2D enzymatic activity, leading to diminished global H3K4 methylation in germinal-center (GC) B cells and DLBCL cells. Conditional deletion of Kmt2d early during B cell development, but not after initiation of the GC reaction, results in an increase in GC B cells and enhances B cell proliferation in mice. Moreover, genetic ablation of Kmt2d in mice overexpressing Bcl2 increases the incidence of GC-derived lymphomas resembling human tumors. These findings suggest that KMT2D acts as a tumor suppressor gene whose early loss facilitates lymphomagenesis by remodeling the epigenetic landscape of the cancer precursor cells. Eradication of KMT2D-deficient cells may thus represent a rational therapeutic approach for targeting early tumorigenic events.

HNRNPA2B1 Is a Mediator of m(6)A-Dependent Nuclear RNA Processing Events

N(6)-methyladenosine (m(6)A) is the most abundant internal modification of messenger RNA. While the presence of m(6)A on transcripts can impact nuclear RNA fates, a reader of this mark that mediates processing of nuclear transcripts has not been identified. We find that the RNA-binding protein HNRNPA2B1 binds m(6)A-bearing RNAs in vivo and in vitro and its biochemical footprint matches the m(6)A consensus motif. HNRNPA2B1 directly binds a set of nuclear transcripts and elicits similar alternative splicing effects as the m(6)A writer METTL3. Moreover, HNRNPA2B1 binds to m(6)A marks in a subset of primary miRNA transcripts, interacts with the microRNA Microprocessor complex protein DGCR8, and promotes primary miRNA processing. Also, HNRNPA2B1 loss and METTL3 depletion cause similar processing defects for these pri-miRNA precursors. We propose HNRNPA2B1 to be a nuclear reader of the m(6)A mark and to mediate, in part, this mark's effects on primary microRNA processing and alternative splicing. PAPERCLIP.

Multiplexing of ChIP-Seq Samples in an Optimized Experimental Condition Has Minimal Impact on Peak Detection

Multiplexing samples in sequencing experiments is a common approach to maximize information yield while minimizing cost. In most cases the number of samples that are multiplexed is determined by financial consideration or experimental convenience, with limited understanding on the effects on the experimental results. Here we set to examine the impact of multiplexing ChIP-seq experiments on the ability to identify a specific epigenetic modification. We performed peak detection analyses to determine the effects of multiplexing. These include false discovery rates, size, position and statistical significance of peak detection, and changes in gene annotation. We found that, for histone marker H3K4me3, one can multiplex up to 8 samples (7 IP + 1 input) at ~21 million single-end reads each and still detect over 90% of all peaks found when using a full lane for sample (~181 million reads). Furthermore, there are no variations introduced by indexing or lane batch effects and importantly there is no significant reduction in the number of genes with neighboring H3K4me3 peaks. We conclude that, for a well characterized antibody and, therefore, model IP condition, multiplexing 8 samples per lane is sufficient to capture most of the biological signal.

The chromatin Remodeler CHD8 is required for activation of progesterone receptor-dependent enhancers

While the importance of gene enhancers in transcriptional regulation is well established, the mechanisms and the protein factors that determine enhancers activity have only recently begun to be unravelled. Recent studies have shown that progesterone receptor (PR) binds regions that display typical features of gene enhancers. Here, we show by ChIP-seq experiments that the chromatin remodeler CHD8 mostly binds promoters under proliferation conditions. However, upon progestin stimulation, CHD8 re-localizes to PR enhancers also enriched in p300 and H3K4me1. Consistently, CHD8 depletion severely impairs progestin-dependent gene regulation. CHD8 binding is PR-dependent but independent of the pioneering factor FOXA1. The SWI/SNF chromatin-remodelling complex is required for PR-dependent gene activation. Interestingly, we show that CHD8 interacts with the SWI/SNF complex and that depletion of BRG1 and BRM, the ATPases of SWI/SNF complex, impairs CHD8 recruitment. We also show that CHD8 is not required for H3K27 acetylation, but contributes to increase accessibility of the enhancer to DNaseI. Furthermore, CHD8 was required for RNAPII recruiting to the enhancers and for transcription of enhancer-derived RNAs (eRNAs). Taken together our data demonstrate that CHD8 is involved in late stages of PR enhancers activation.

A lot of research has been devoted during the last decades to understand the mechanisms that control gene promoters activity, however, much less is known about enhancers. Only recently, the use of genome-wide chromatin immunoprecipitation techniques has revealed the existence of more than 400,000 enhancers in the human genome. We are starting to understand the importance of these regulatory elements and how they are activated or repressed. In this work we discover that the chromatin remodeler CHD8 is recruited to Progesteron Receptor-dependent enhancers upon hormone treatment. CHD8 is required for late steps in the activation of these enhancers, including transcription of the enhancers and synthesis of eRNA (long noncoding RNAs derived form the enhancers).

Epigenomic evolution in diffuse large B-cell lymphomas

The contribution of epigenomic alterations to tumour progression and relapse is not well characterized. Here we characterize an association between disease progression and DNA methylation in diffuse large B-cell lymphoma (DLBCL). By profiling genome-wide DNA methylation at single-base pair resolution in thirteen DLBCL diagnosis–relapse sample pairs, we show that DLBCL patients exhibit heterogeneous evolution of tumour methylomes during relapse. We identify differentially methylated regulatory elements and determine a relapse-associated methylation signature converging on key pathways such as transforming growth factor-β (TGF-β) receptor activity. We also observe decreased intra-tumour methylation heterogeneity from diagnosis to relapsed tumour samples. Relapse-free patients display lower intra-tumour methylation heterogeneity at diagnosis compared with relapsed patients in an independent validation cohort. Furthermore, intra-tumour methylation heterogeneity is predictive of time to relapse. Therefore, we propose that epigenomic heterogeneity may support or drive the relapse phenotype and can be used to predict DLBCL relapse.

The contribution of epigenomic alterations to tumour progression and relapse is not well characterized. Here the authors characterize epigenetic evolution in aggressive B-cell lymphoma and find that epigenomic heterogeneity may not only support and drive the relapse phenotype but also be used to predict lymphoma relapse.

BCL6 orchestrates Tfh cell differentiation via multiple distinct mechanisms

Follicular helper T cells (Tfh cells) are required for T cell help to B cells, and BCL6 is the defining transcription factor of Tfh cells. However, the functions of BCL6 in Tfh cells have largely remained unclear. Here we defined the BCL6 cistrome in primary human germinal center Tfh cells to assess mechanisms of BCL6 regulation of CD4 T cells, comparing and contrasting BCL6 function in T and B cells. BCL6 primarily acts as a repressor in Tfh cells, and BCL6 binding was associated with control of Tfh cell migration and repression of alternative cell fates. Interestingly, although some BCL6-bound genes possessed BCL6 DNA-binding motifs, many BCL6-bound loci were instead characterized by the presence of DNA motifs for AP1 or STAT. AP1 complexes are key positive downstream mediators of TCR signaling and external stimuli. We show that BCL6 can directly bind AP1, and BCL6 depends on AP1 for recruitment to BCL6-binding sites with AP1 motifs, suggesting that BCL6 subverts AP1 activity. These findings reveal that BCL6 has broad and multifaceted effects on Tfh biology and provide insight into how this master regulator mediates distinct cell context-dependent phenotypes.

N6-methyladenosine marks primary microRNAs for processing

The first step in the biogenesis of microRNAs is the processing of primary microRNAs (pri-miRNAs) by the microprocessor complex, composed of the RNA-binding protein DGCR8 and the type III RNase DROSHA. This initial event requires recognition of the junction between the stem and the flanking single-stranded RNA of the pri-miRNA hairpin by DGCR8 followed by recruitment of DROSHA, which cleaves the RNA duplex to yield the pre-miRNA product. While the mechanisms underlying pri-miRNA processing have been determined, the mechanism by which DGCR8 recognizes and binds pri-miRNAs, as opposed to other secondary structures present in transcripts, is not understood. Here we find in mammalian cells that methyltransferase-like 3 (METTL3) methylates pri-miRNAs, marking them for recognition and processing by DGCR8. Consistent with this, METTL3 depletion reduced the binding of DGCR8 to pri-miRNAs and resulted in the global reduction of mature miRNAs and concomitant accumulation of unprocessed pri-miRNAs. In vitro processing reactions confirmed the sufficiency of the N(6)-methyladenosine (m(6)A) mark in promoting pri-miRNA processing. Finally, gain-of-function experiments revealed that METTL3 is sufficient to enhance miRNA maturation in a global and non-cell-type-specific manner. Our findings reveal that the m(6)A mark acts as a key post-transcriptional modification that promotes the initiation of miRNA biogenesis.

RNA-Seq and ChIP-Seq reveal SQSTM1/p62 as a key mediator of JunB suppression of NF-κB-dependent inflammation

Mice with epidermal deletion of JunB transcription factor displayed a psoriasis-like inflammation. The relevance of these findings to humans and the mechanisms mediating JunB function are not fully understood. Here we demonstrate that impaired JunB function via gene silencing or overexpression of a dominant negative mutant increased human keratinocyte cell proliferation but decreased cell barrier function. RNA-seq revealed over 500 genes affected by JunB loss of function, which included the upregulation of an array of proinflammatory molecules relevant to psoriasis. Among these were tumor necrosis factor α (TNFα), CCL2, CXCL10, IL6R, and SQSTM1, an adaptor protein involved in nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) activation. Chromatin immunoprecipitation (ChIP)-Seq and gene reporter analyses showed that JunB directly suppressed SQSTM1 by binding to a consensus AP-1 cis element located around 2 kb upstream of SQSTM1-transcription start site. Similar to JunB loss of function, SQSTM1-overexpression induced TNFα, CCL2, and CXCL10. Conversely, NF-κB inhibition genetically with a mutant IκBα or pharmacologically with pyrrolidine dithiocarbamate (PDTC) prevented cytokine, but not IL6R, induction by JunB deficiency. Taken together, our findings indicate that JunB controls epidermal growth, barrier formation, and proinflammatory responses through direct and indirect mechanisms, pinpointing SQSTM1 as a key mediator of JunB suppression of NF-κB-dependent inflammation.

DNA hydroxymethylation profiling reveals that WT1 mutations result in loss of TET2 function in acute myeloid leukemia

Somatic mutations in IDH1/IDH2 and TET2 result in impaired TET2-mediated conversion of 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC). The observation that WT1 inactivating mutations anticorrelate with TET2/IDH1/IDH2 mutations in acute myeloid leukemia (AML) led us to hypothesize that WT1 mutations may impact TET2 function. WT1 mutant AML patients have reduced 5hmC levels similar to TET2/IDH1/IDH2 mutant AML. These mutations are characterized by convergent, site-specific alterations in DNA hydroxymethylation, which drive differential gene expression more than alterations in DNA promoter methylation. WT1 overexpression increases global levels of 5hmC, and WT1 silencing reduced 5hmC levels. WT1 physically interacts with TET2 and TET3, and WT1 loss of function results in a similar hematopoietic differentiation phenotype as observed with TET2 deficiency. These data provide a role for WT1 in regulating DNA hydroxymethylation and suggest that TET2 IDH1/IDH2 and WT1 mutations define an AML subtype defined by dysregulated DNA hydroxymethylation.

The oestrogen receptor alpha-regulated lncRNA NEAT1 is a critical modulator of prostate cancer

The androgen receptor (AR) plays a central role in establishing an oncogenic cascade that drives prostate cancer progression. Some prostate cancers escape androgen dependence and are often associated with an aggressive phenotype. The oestrogen receptor alpha (ERα) is expressed in prostate cancers, independent of AR status. However, the role of ERα remains elusive. Using a combination of chromatin immunoprecipitation (ChIP) and RNA-sequencing data, we identified an ERα-specific non-coding transcriptome signature. Among putatively ERα-regulated intergenic long non-coding RNAs (lncRNAs), we identified nuclear enriched abundant transcript 1 (NEAT1) as the most significantly overexpressed lncRNA in prostate cancer. Analysis of two large clinical cohorts also revealed that NEAT1 expression is associated with prostate cancer progression. Prostate cancer cells expressing high levels of NEAT1 were recalcitrant to androgen or AR antagonists. Finally, we provide evidence that NEAT1 drives oncogenic growth by altering the epigenetic landscape of target gene promoters to favour transcription.

jChIP: a graphical environment for exploratory ChIP-Seq data analysis

Background

Chromatin immunoprecipitation coupled with next-generation sequencing (ChIP-Seq) provides a powerful tool for discovering protein-DNA interactions. Still, the computational analysis of the great amount of ChIP-Seq data generated, involving mapping of raw data to reference genome, has been a bottle neck for most of researchers in the transcriptional and epigenetic fields. Thus, user-friendly ChIP-Seq processing method sare much needed to enable greater community of computational and bench biologists to exploit the power of ChIP-Seq technology .

Findings

jChIP is a graphical tool that was developed to analyze and display ChIP-Seq data. It matches reads to the corresponding loci downloaded from Ensembl Genes or Ensembl Regulation databases. jChIP provides a friendly interface for exploratory analysis of mapped reads as well as peak calling data. The built-in functions for graphical display of reads distribution allows to evaluate the quality and meaning of ChIP-Seq data.

Conclusion

jChIP is a user-friendly GUI-based software for the analysis of ChIP-Seq data within context of known genomic features. Further, jChIP provides tools for discovering new and refining known genome-wide protein binding patterns.

Histone methyltransferase MMSET/NSD2 alters EZH2 binding and reprograms the myeloma epigenome through global and focal changes in H3K36 and H3K27 methylation

Overexpression of the histone methyltransferase MMSET in t(4;14)+ multiple myeloma patients is believed to be the driving factor in the pathogenesis of this subtype of myeloma. MMSET catalyzes dimethylation of lysine 36 on histone H3 (H3K36me2), and its overexpression causes a global increase in H3K36me2, redistributing this mark in a broad, elevated level across the genome. Here, we demonstrate that an increased level of MMSET also induces a global reduction of lysine 27 trimethylation on histone H3 (H3K27me3). Despite the net decrease in H3K27 methylation, specific genomic loci exhibit enhanced recruitment of the EZH2 histone methyltransferase and become hypermethylated on this residue. These effects likely contribute to the myeloma phenotype since MMSET-overexpressing cells displayed increased sensitivity to EZH2 inhibition. Furthermore, we demonstrate that such MMSET-mediated epigenetic changes require a number of functional domains within the protein, including PHD domains that mediate MMSET recruitment to chromatin. In vivo, targeting of MMSET by an inducible shRNA reversed histone methylation changes and led to regression of established tumors in athymic mice. Together, our work elucidates previously unrecognized interplay between MMSET and EZH2 in myeloma oncogenesis and identifies domains to be considered when designing inhibitors of MMSET function.

Zfx facilitates tumorigenesis caused by activation of the Hedgehog pathway

The Hedgehog (Hh) signaling pathway regulates normal development and cell proliferation in metazoan organisms, but its aberrant activation can promote tumorigenesis. Hh-induced tumors arise from various tissues and they may be indolent or aggressive, as is the case with skin basal cell carcinoma (BCC) or cerebellar medulloblastoma, respectively. Little is known about common cell-intrinsic factors that control the development of such diverse Hh-dependent tumors. Transcription factor Zfx is required for the self-renewal of hematopoietic and embryonic stem cells, as well as for the propagation of acute myeloid and T-lymphoblastic leukemias. We report here that Zfx facilitates the development of experimental BCC and medulloblastoma in mice initiated by deletion of the Hh inhibitory receptor Ptch1. Simultaneous deletion of Zfx along with Ptch1 prevented BCC formation and delayed medulloblastoma development. In contrast, Zfx was dispensable for tumorigenesis in a mouse model of glioblastoma. We used genome-wide expression and chromatin-binding analysis in a human medulloblastoma cell line to characterize direct, evolutionarily conserved targets of Zfx, identifying Dis3L and Ube2j1 as two targets required for the growth of the human medulloblastoma cells. Our results establish Zfx as a common cell-intrinsic regulator of diverse Hh-induced tumors, with implications for the definition of new therapeutic targets in these malignancies.

Metastasis-suppressor transcript destabilization through TARBP2 binding of mRNA hairpins

Aberrant regulation of RNA stability has an important role in many disease states. Deregulated post-transcriptional modulation, such as that governed by microRNAs targeting linear sequence elements in messenger RNAs, has been implicated in the progression of many cancer types. A defining feature of RNA is its ability to fold into structures. However, the roles of structural mRNA elements in cancer progression remain unexplored. Here we performed an unbiased search for post-transcriptional modulators of mRNA stability in breast cancer by conducting whole-genome transcript stability measurements in poorly and highly metastatic isogenic human breast cancer lines. Using a computational framework that searches RNA sequence and structure space, we discovered a family of GC-rich structural cis-regulatory RNA elements, termed sRSEs for structural RNA stability elements, which are significantly overrepresented in transcripts displaying reduced stability in highly metastatic cells. By integrating computational and biochemical approaches, we identified TARBP2, a double-stranded RNA-binding protein implicated in microRNA processing, as the trans factor that binds the sRSE family and similar structural elements--collectively termed TARBP2-binding structural elements (TBSEs)--in transcripts. TARBP2 is overexpressed in metastatic cells and metastatic human breast tumours and destabilizes transcripts containing TBSEs. Endogenous TARBP2 promotes metastatic cell invasion and colonization by destabilizing amyloid precursor protein (APP) and ZNF395 transcripts, two genes previously associated with Alzheimer's and Huntington's disease, respectively. We reveal these genes to be novel metastasis suppressor genes in breast cancer. The cleavage product of APP, extracellular amyloid-α peptide, directly suppresses invasion while ZNF395 transcriptionally represses a pro-metastatic gene expression program. The expression levels of TARBP2, APP and ZNF395 in human breast carcinomas support their experimentally uncovered roles in metastasis. Our findings establish a non-canonical and direct role for TARBP2 in mammalian gene expression regulation and reveal that regulated RNA destabilization through protein-mediated binding of mRNA structural elements can govern cancer progression.

Probabilistic partitioning methods to find significant patterns in ChIP-Seq data

MOTIVATION

We have witnessed an enormous increase in ChIP-Seq data for histone modifications in the past few years. Discovering significant patterns in these data is an important problem for understanding biological mechanisms.

RESULTS

We propose probabilistic partitioning methods to discover significant patterns in ChIP-Seq data. Our methods take into account signal magnitude, shape, strand orientation and shifts. We compare our methods with some current methods and demonstrate significant improvements, especially with sparse data. Besides pattern discovery and classification, probabilistic partitioning can serve other purposes in ChIP-Seq data analysis. Specifically, we exemplify its merits in the context of peak finding and partitioning of nucleosome positioning patterns in human promoters.

AVAILABILITY AND IMPLEMENTATION

The software and code are available in the supplementary material.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

High-resolution chromatin immunoprecipitation (ChIP) sequencing reveals novel binding targets and prognostic role for SOX11 in mantle cell lymphoma

Sex determining region Y-box 11 (SOX11) expression is specific for mantle cell lymphoma (MCL) as compared with other non-Hodgkin's lymphomas. However, the function and direct-binding targets of SOX11 in MCL are largely unknown. We used high-resolution chromatin immunoprecipitation sequencing to identify the direct target genes of SOX11 in a genome-wide, unbiased manner and elucidate its functional significance. Pathway analysis identified WNT, PKA and TGF-beta signaling pathways as significantly enriched by SOX11-target genes. Quantitative chromatin immunoprecipitation sequencing and promoter reporter assays confirmed that SOX11 directly binds to individual genes and modulates their transcription activities in these pathways in MCL. Functional studies using RNA interference demonstrate that SOX11 directly regulates WNT in MCL. We analyzed SOX11 expression in three independent well-annotated tissue microarrays from the University of Wisconsin (UW), Karolinska Institute and British Columbia Cancer Agency. Our findings suggest that high SOX11 expression is associated with improved survival in a subset of MCL patients, particularly those treated with intensive chemotherapy. Transcriptional regulation of WNT and other biological pathways affected by SOX11-target genes may help explain the impact of SOX11 expression on patient outcomes.

The transcriptional modulator BCL6 as a molecular target for breast cancer therapy

Inappropriate expression or activation of transcription factors can drive patterns of gene expression, leading to the malignant behavior of breast cancer cells. We have found that the transcriptional repressor BCL6 is highly expressed in breast cancer cell lines, and its locus is amplified in about half of primary breast cancers. To understand how BCL6 regulates gene expression in breast cancer cells, we used chromatin immunoprecipitation followed by deep sequencing to identify the BCL6 binding sites on a genomic scale. This revealed that BCL6 regulates a unique cohort of genes in breast cancer cell lines compared with B-cell lymphomas. Furthermore, BCL6 expression promotes the survival of breast cancer cells, and targeting BCL6 with a peptidomimetic inhibitor leads to apoptosis of these cells. Finally, combining a BCL6 inhibitor and a signal transducer and activator of transcription3 inhibitor provided enhanced cell killing in triple-negative breast cancer cell lines, suggesting that combination therapy may be particularly useful. Thus, targeting BCL6 alone or in conjunction with other signaling pathways may be a useful therapeutic strategy for treating breast cancer.

C26 cancer-induced muscle wasting is IKKβ-dependent and NF-kappaB-independent

Existing data suggest that NF-kappaB signaling is a key regulator of cancer-induced skeletal muscle wasting. However, identification of the components of this signaling pathway and of the NF-κB transcription factors that regulate wasting is far from complete. In muscles of C26 tumor bearing mice, overexpression of dominant negative (d.n.) IKKβ blocked muscle wasting by 69% and the IκBα-super repressor blocked wasting by 41%. In contrast, overexpression of d.n. IKKα or d.n. NIK did not block C26-induced wasting. Surprisingly, overexpression of d.n. p65 or d.n. c-Rel did not significantly affect muscle wasting. Genome-wide mRNA expression arrays showed upregulation of many genes previously implicated in muscle atrophy. To test if these upregulated genes were direct targets of NF-κB transcription factors, we compared genome-wide p65 binding to DNA in control and cachectic muscle using ChIP-sequencing. Bioinformatic analysis of ChIP-sequencing data from control and C26 muscles showed very little p65 binding to genes in cachexia and little to suggest that upregulated p65 binding influences the gene expression associated with muscle based cachexia. The p65 ChIP-seq data are consistent with our finding of no significant change in protein binding to an NF-κB oligonucleotide in a gel shift assay, no activation of a NF-κB-dependent reporter, and no effect of d.n.p65 overexpression in muscles of tumor bearing mice. Taken together, these data support the idea that although inhibition of IκBα, and particularly IKKβ, blocks cancer-induced wasting, the alternative NF-κB signaling pathway is not required. In addition, the downstream NF-κB transcription factors, p65 and c-Rel do not appear to regulate the transcriptional changes induced by the C26 tumor. These data are consistent with the growing body of literature showing that there are NF-κB-independent substrates of IKKβ and IκBα that regulate physiological processes.

Detection and classification of peaks in 5' cap RNA sequencing data

BACKGROUND

The large-scale sequencing of 5' cap enriched cDNA promises to reveal the diversity of transcription initiation across entire genomes. The process of transcription is noisy, and there is often no single, exact start site. This creates the need for a fast and simple method of identifying transcription start peaks based on this type of data. Due to both biological and technical noise, many of the peaks seen are not real transcription initiation events. Classification of the observed peaks is an essential filtering step in the discovery of genuine initiation locations.

RESULTS

We develop a two-stage approach consisting of a fast and simple algorithm based on a sliding window with Poisson null distribution for detecting the genomic locations of peaks, followed by a linear support vector machine classifier to distinguish between peaks which represent the initiation of transcription and peaks that do not. Comparison of classification performance to the best existing method based on whole genome segmentation showed comparable precision and improved recall. Internal features, which are intrinsic to the data and require no further experiments, had high precision and recall rates. Addition of pooled external data or matched RNA sequencing data resulted in gains of recall with equivalent precision.

CONCLUSIONS

The Poisson sliding window model is an effective and fast way of taking the peak neighbourhood into account, and finding statistically significant peaks over a range of transcript expression values. It is orders of magnitude faster than doing whole genome segmentation. The support vector classification scheme has better precision and recall than existing methods. Integrating additional datasets is shown to provide minor gains in recall, in comparison to using only the cap-sequencing data.

Nencki Genomics Database--Ensembl funcgen enhanced with intersections, user data and genome-wide TFBS motifs

We present the Nencki Genomics Database, which extends the functionality of Ensembl Regulatory Build (funcgen) for the three species: human, mouse and rat. The key enhancements over Ensembl funcgen include the following: (i) a user can add private data, analyze them alongside the public data and manage access rights; (ii) inside the database, we provide efficient procedures for computing intersections between regulatory features and for mapping them to the genes. To Ensembl funcgen-derived data, which include data from ENCODE, we add information on conserved non-coding (putative regulatory) sequences, and on genome-wide occurrence of transcription factor binding site motifs from the current versions of two major motif libraries, namely, Jaspar and Transfac. The intersections and mapping to the genes are pre-computed for the public data, and the result of any procedure run on the data added by the users is stored back into the database, thus incrementally increasing the body of pre-computed data. As the Ensembl funcgen schema for the rat is currently not populated, our database is the first database of regulatory features for this frequently used laboratory animal. The database is accessible without registration using the mysql client: mysql –h database.nencki-genomics.org –u public. Registration is required only to add or access private data. A WSDL webservice provides access to the database from any SOAP client, including the Taverna Workbench with a graphical user interface.

Database URL:http://www.nencki-genomics.org.