Saturation analysis of ChIP-seq data for reproducible identification of binding peaks

peaksat: an R package for ChIP-seq peak saturation analysis

BACKGROUND

Epigenomic profiling assays such as ChIP-seq have been widely used to map the genome-wide enrichment profiles of chromatin-associated proteins and posttranslational histone modifications. Sequencing depth is a key parameter in experimental design and quality control. However, due to variable sequencing depth requirements across experimental conditions, it can be challenging to determine optimal sequencing depth, particularly for projects involving multiple targets or cell types.

RESULTS

We developed the peaksat R package to provide target read depth estimates for epigenomic experiments based on the analysis of peak saturation curves. We applied peaksat to establish the distinctive read depth requirements for ChIP-seq studies of histone modifications in different cell lines. Using peaksat, we were able to estimate the target read depth required per library to obtain high-quality peak calls for downstream analysis. In addition, peaksat was applied to other sequence-enrichment methods including CUT&RUN and ATAC-seq.

CONCLUSION

peaksat addresses a need for researchers to make informed decisions about whether their sequencing data has been generated to an adequate depth and subsequently sufficient meaningful peaks, and failing that, how many more reads would be required per library. peaksat is applicable to other sequence-based methods that include calling peaks in their analysis.

MITF Contributes to the Body Color Differentiation of Sea Cucumbers Apostichopus japonicus through Expression Differences and Regulation of Downstream Genes

Melanin, which is a pigment produced in melanocytes, is an important contributor to sea cucumber body color. MITF is one of the most critical genes in melanocyte development and melanin synthesis pathways. However, how MITF regulates body color and differentiation in sea cucumbers is poorly understood. In this study, we analyzed the expression level and location of MITF in white, purple, and green sea cucumbers and identified the genes regulated by MITF using chromatin immunoprecipitation followed by sequencing. The mRNA and protein expression levels of MITF were all highest in purple morphs and lowest in white morphs. In situ hybridization indicated that MITF mRNA were mainly expressed in the epidermis. We also identified 984, 732, and 1191 peaks of MITF binding in green, purple, and white sea cucumbers, which were associated with 727, 557, and 887 genes, respectively. Our findings suggested that MITF contributed to the body color differentiation of green, purple, and white sea cucumbers through expression differences and regulation of downstream genes. These results provided a basis for future studies to determine the mechanisms underlying body color formation and provided insights into gene regulation in sea cucumbers.

Pervasive and CpG-dependent promoter-like characteristics of transcribed enhancers

The temporal and spatial expression of genes is controlled by promoters and enhancers. Findings obtained over the last decade that not only promoters but also enhancers are characterized by bidirectional, divergent transcription have challenged the traditional notion that promoters and enhancers represent distinct classes of regulatory elements. Over half of human promoters are associated with CpG islands (CGIs), relatively CpG-rich stretches of generally several hundred nucleotides that are often associated with housekeeping genes. Only about 6% of transcribed enhancers defined by CAGE-tag analysis are associated with CGIs. Here, we present an analysis of enhancer and promoter characteristics and relate them to the presence or absence of CGIs. We show that transcribed enhancers share a number of CGI-dependent characteristics with promoters, including statistically significant local overrepresentation of core promoter elements. CGI-associated enhancers are longer, display higher directionality of transcription, greater expression, a lesser degree of tissue specificity, and a higher frequency of transcription-factor binding events than non-CGI-associated enhancers. Genes putatively regulated by CGI-associated enhancers are enriched for transcription regulator activity. Our findings show that CGI-associated transcribed enhancers display a series of characteristics related to sequence, expression and function that distinguish them from enhancers not associated with CGIs.

7C: Computational Chromosome Conformation Capture by Correlation of ChIP-seq at CTCF motifs

BACKGROUND

Knowledge of the three-dimensional structure of the genome is necessary to understand how gene expression is regulated. Recent experimental techniques such as Hi-C or ChIA-PET measure long-range chromatin interactions genome-wide but are experimentally elaborate, have limited resolution and such data is only available for a limited number of cell types and tissues.

RESULTS

While ChIP-seq was not designed to detect chromatin interactions, the formaldehyde treatment in the ChIP-seq protocol cross-links proteins with each other and with DNA. Consequently, also regions that are not directly bound by the targeted TF but interact with the binding site via chromatin looping are co-immunoprecipitated and sequenced. This produces minor ChIP-seq signals at loop anchor regions close to the directly bound site. We use the position and shape of ChIP-seq signals around CTCF motif pairs to predict whether they interact or not. We implemented this approach in a prediction method, termed Computational Chromosome Conformation Capture by Correlation of ChIP-seq at CTCF motifs (7C). We applied 7C to all CTCF motif pairs within 1 Mb in the human genome and validated predicted interactions with high-resolution Hi-C and ChIA-PET. A single ChIP-seq experiment from known architectural proteins (CTCF, Rad21, Znf143) but also from other TFs (like TRIM22 or RUNX3) predicts loops accurately. Importantly, 7C predicts loops in cell types and for TF ChIP-seq datasets not used in training.

CONCLUSION

7C predicts chromatin loops which can help to associate TF binding sites to regulated genes. Furthermore, profiling of hundreds of ChIP-seq datasets results in novel candidate factors functionally involved in chromatin looping. Our method is available as an R/Bioconductor package: http://bioconductor.org/packages/sevenC .

Sensitive and robust assessment of ChIP-seq read distribution using a strand-shift profile

Motivation

Chromatin immunoprecipitation followed by sequencing (ChIP-seq) can detect read-enriched DNA loci for point-source (e.g. transcription factor binding) and broad-source factors (e.g. various histone modifications). Although numerous quality metrics for ChIP-seq data have been developed, the 'peaks' thus obtained are still difficult to assess with respect to signal-to-noise ratio (S/N) and the percentage of false positives.

Results

We developed a quality-assessment tool for ChIP-seq data, strand-shift profile (SSP), which quantifies S/N and peak reliability without peak calling. We validated SSP in-depth using  ≥ 1000 publicly available ChIP-seq datasets along with virtual data to demonstrate that SSP provides a quantifiable and sensitive score to different S/Ns for both point- and broad-source factors, which can be standardized across diverse cell types and read depths. SSP also provides an effective criterion to judge whether a specific normalization or a rejection is required for each sample, which cannot be estimated by quality metrics currently available. Finally, we show that 'hidden-duplicate reads' cause aberrantly high S/Ns, and SSP provides an additional metric to avoid them, which can also contribute to estimation of peak mode (point- or broad-source) of samples.

Availability and implementation

SSP is open source software written in C++ and can be downloaded at https://github.com/rnakato/SSP.

Supplementary information

Supplementary data are available at Bioinformatics online.

Rapid Response and Slow Recovery of the H3K4me3 Epigenomic Marker in the Liver after Light-mediated Phase Advances of the Circadian Clock

Mammalian tissues display circadian rhythms in transcription, translation, and histone modifications. Here we asked how an advance of the light-dark cycle alters daily rhythms in the liver epigenome at the H3K4me3 (trimethylation of lysine 4 on histone 3) modification, which is found at active and poised gene promoters. H3K4me3 levels were first measured at 4 time points (zeitgeber time [ZT] 3, 8, 15, and 20) during a normal 12L:12D light-dark cycle. Peak levels were observed during the early dark phase at ZT15 and dropped to low levels around lights-on (ZT0) between ZT20 and ZT3. A 6-h phase advance at ZT18 (new lights-on after only 6 h of darkness) led to a transient extension of peak H3K4me3 levels. Although locomotor activity reentrained within a week after the phase advance, H3K4me3 rhythms failed to do so, with peak levels remaining in the light phase at the 1-week recovery time point. Eight weekly phase advances, with 1-week recovery times between each phase advance, further disrupted the H3K4me3 rhythms. Finally, we used the mPer2^Luc knockin mouse to determine whether the phase advance also disrupted Per2 protein expression. Similar to the results from the histone work, we found both a rapid response to the phase advance and a delayed recovery, the latter in sync with H3K4me3 levels. A model to explain these results is offered.

Analysis of the androgen receptor-regulated lncRNA landscape identifies a role for ARLNC1 in prostate cancer progression

The androgen receptor (AR) plays a critical role in the development of the normal prostate as well as prostate cancer. Using an integrative transcriptomic analysis of prostate cancer cell lines and tissues, we identified ARLNC1 (AR-regulated long non-coding RNA 1) as an important long non-coding RNA that is strongly associated with AR signaling in prostate cancer progression. Not only was ARLNC1 induced by AR protein, ARLNC1 stabilized the AR transcript via RNA-RNA interaction. ARLNC1 knockdown suppressed AR expression, global AR signaling, and prostate cancer growth in vitro and in vivo. Taken together, these data support a role for ARLNC1 in maintaining a positive feedback loop that potentiates AR signaling during prostate cancer progression, and identifies ARLNC1 as a novel therapeutic target.

Ritornello: high fidelity control-free chromatin immunoprecipitation peak calling

With the advent of next generation high-throughput DNA sequencing technologies, omics experiments have become the mainstay for studying diverse biological effects on a genome wide scale. Chromatin immunoprecipitation (ChIP-seq) is the omics technique that enables genome wide localization of transcription factor (TF) binding or epigenetic modification events. Since the inception of ChIP-seq in 2007, many methods have been developed to infer ChIP-target binding loci from the resultant reads after mapping them to a reference genome. However, interpreting these data has proven challenging, and as such these algorithms have several shortcomings, including susceptibility to false positives due to artifactual peaks, poor localization of binding sites and the requirement for a total DNA input control which increases the cost of performing these experiments. We present Ritornello, a new approach for finding TF-binding sites in ChIP-seq, with roots in digital signal processing that addresses all of these problems. We show that Ritornello generally performs equally or better than the peak callers tested and recommended by the ENCODE consortium, but in contrast, Ritornello does not require a matched total DNA input control to avoid false positives, effectively decreasing the sequencing cost to perform ChIP-seq. Ritornello is freely available at https://github.com/KlugerLab/Ritornello.

Recent advances in ChIP-seq analysis: from quality management to whole-genome annotation

Chromatin immunoprecipitation followed by sequencing (ChIP-seq) analysis can detect protein/DNA-binding and histone-modification sites across an entire genome. Recent advances in sequencing technologies and analyses enable us to compare hundreds of samples simultaneously; such large-scale analysis has potential to reveal the high-dimensional interrelationship level for regulatory elements and annotate novel functional genomic regions de novo. Because many experimental considerations are relevant to the choice of a method in a ChIP-seq analysis, the overall design and quality management of the experiment are of critical importance. This review offers guiding principles of computation and sample preparation for ChIP-seq analyses, highlighting the validity and limitations of the state-of-the-art procedures at each step. We also discuss the latest challenges of single-cell analysis that will encourage a new era in this field.

Bioinformatics Tools for Genome-Wide Epigenetic Research

Epigenetics play a central role in the regulation of many important cellular processes, and dysregulations at the epigenetic level could be the source of serious pathologies, such as neurological disorders affecting brain development, neurodegeneration, and intellectual disability. Despite significant technological advances for epigenetic profiling, there is still a need for a systematic understanding of how epigenetics shapes cellular circuitry, and disease pathogenesis. The development of accurate computational approaches for analyzing complex epigenetic profiles is essential for disentangling the mechanisms underlying cellular development, and the intricate interaction networks determining and sensing chromatin modifications and DNA methylation to control gene expression. In this chapter, we review the recent advances in the field of "computational epigenetics," including computational methods for processing different types of epigenetic data, prediction of chromatin states, and study of protein dynamics. We also discuss how "computational epigenetics" has complemented the fast growth in the generation of epigenetic data for uncovering the main differences and similarities at the epigenetic level between individuals and the mechanisms underlying disease onset and progression.

The SeqAn C++ template library for efficient sequence analysis: A resource for programmers

BACKGROUND

The use of novel algorithmic techniques is pivotal to many important problems in life science. For example the sequencing of the human genome (Venter et al., 2001) would not have been possible without advanced assembly algorithms and the development of practical BWT based read mappers have been instrumental for NGS analysis. However, owing to the high speed of technological progress and the urgent need for bioinformatics tools, there was a widening gap between state-of-the-art algorithmic techniques and the actual algorithmic components of tools that are in widespread use. We previously addressed this by introducing the SeqAn library of efficient data types and algorithms in 2008 (Döring et al., 2008).

RESULTS

The SeqAn library has matured considerably since its first publication 9 years ago. In this article we review its status as an established resource for programmers in the field of sequence analysis and its contributions to many analysis tools.

CONCLUSIONS

We anticipate that SeqAn will continue to be a valuable resource, especially since it started to actively support various hardware acceleration techniques in a systematic manner.

The haploinsufficient tumor suppressor, CUX1, acts as an analog transcriptional regulator that controls target genes through distal enhancers that loop to target promoters

One third of tumor suppressors are haploinsufficient transcriptional regulators, yet it remains unknown how a 50% reduction of a transcription factor is translated at the cis-regulatory level into a malignant transcriptional program. We studied CUX1, a haploinsufficient transcription factor that is recurrently mutated in hematopoietic and solid tumors. We determined CUX1 DNA-binding and target gene regulation in the wildtype and haploinsufficient states. CUX1 binds with transcriptional activators and cohesin at distal enhancers across three different human cell types. Haploinsufficiency of CUX1 altered the expression of a large number of genes, including cell cycle regulators, with concomitant increased cellular proliferation. Surprisingly, CUX1 occupancy decreased genome-wide in the haploinsufficient state, and binding site affinity did not correlate with differential gene expression. Instead, differentially expressed genes had multiple, low-affinity CUX1 binding sites, features of analog gene regulation. A machine-learning algorithm determined that chromatin accessibility, enhancer activity, and distance to the transcription start site are features of dose-sensitive CUX1 transcriptional regulation. Moreover, CUX1 is enriched at sites of DNA looping, as determined by Hi-C analysis, and these loops connect CUX1 to the promoters of regulated genes. We propose an analog model for haploinsufficient transcriptional deregulation mediated by higher order genome architecture.

Characterization of hundreds of regulatory landscapes in developing limbs reveals two regimes of chromatin folding

Complex regulatory landscapes control the pleiotropic transcriptional activities of developmental genes. For most genes, the number, location, and dynamics of their associated regulatory elements are unknown. In this work, we characterized the three-dimensional chromatin microarchitecture and regulatory landscape of 446 limb-associated gene loci in mouse using Capture-C, ChIP-seq, and RNA-seq in forelimb, hindlimb at three developmental stages, and midbrain. The fine mapping of chromatin interactions revealed a strong preference for functional genomic regions such as repressed or active domains. By combining chromatin marks and interaction peaks, we annotated more than 1000 putative limb enhancers and their associated genes. Moreover, the analysis of chromatin interactions revealed two regimes of chromatin folding, one producing interactions stable across tissues and stages and another one associated with tissue and/or stage-specific interactions. Whereas stable interactions associate strongly with CTCF/RAD21 binding, the intensity of variable interactions correlates with changes in underlying chromatin modifications, specifically at the viewpoint and at the interaction site. In conclusion, this comprehensive data set provides a resource for the characterization of hundreds of limb-associated regulatory landscapes and a framework to interpret the chromatin folding dynamics observed during embryogenesis.

Q-nexus: a comprehensive and efficient analysis pipeline designed for ChIP-nexus

Background

ChIP-nexus, an extension of the ChIP-exo protocol, can be used to map the borders of protein-bound DNA sequences at nucleotide resolution, requires less input DNA and enables selective PCR duplicate removal using random barcodes. However, the use of random barcodes requires additional preprocessing of the mapping data, which complicates the computational analysis. To date, only a very limited number of software packages are available for the analysis of ChIP-exo data, which have not yet been systematically tested and compared on ChIP-nexus data.

Results

Here, we present a comprehensive software package for ChIP-nexus data that exploits the random barcodes for selective removal of PCR duplicates and for quality control. Furthermore, we developed bespoke methods to estimate the width of the protected region resulting from protein-DNA binding and to infer binding positions from ChIP-nexus data. Finally, we applied our peak calling method as well as the two other methods MACE and MACS2 to the available ChIP-nexus data.

Conclusions

The Q-nexus software is efficient and easy to use. Novel statistics about duplication rates in consideration of random barcodes are calculated. Our method for the estimation of the width of the protected region yields unbiased signatures that are highly reproducible for biological replicates and at the same time very specific for the respective factors analyzed. As judged by the irreproducible discovery rate (IDR), our peak calling algorithm shows a substantially better reproducibility. An implementation of Q-nexus is available at http://charite.github.io/Q/.

Electronic supplementary material

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

Transcriptional Regulation by ATOH1 and its Target SPDEF in the Intestine

BACKGROUND & AIMS

The transcription factor atonal homolog 1 (ATOH1) controls the fate of intestinal progenitors downstream of the Notch signaling pathway. Intestinal progenitors that escape Notch activation express high levels of ATOH1 and commit to a secretory lineage fate, implicating ATOH1 as a gatekeeper for differentiation of intestinal epithelial cells. Although some transcription factors downstream of ATOH1, such as SPDEF, have been identified to specify differentiation and maturation of specific cell types, the bona fide transcriptional targets of ATOH1 still largely are unknown. Here, we aimed to identify ATOH1 targets and to identify transcription factors that are likely to co-regulate gene expression with ATOH1.

METHODS

We used a combination of chromatin immunoprecipitation and messenger RNA-based high-throughput sequencing (ChIP-seq and RNA-seq), together with cell sorting and transgenic mice, to identify direct targets of ATOH1, and establish the epistatic relationship between ATOH1 and SPDEF.

RESULTS

By using unbiased genome-wide approaches, we identified more than 700 genes as ATOH1 transcriptional targets in adult small intestine and colon. Ontology analysis indicated that ATOH1 directly regulates genes involved in specification and function of secretory cells. De novo motif analysis of ATOH1 targets identified SPDEF as a putative transcriptional co-regulator of ATOH1. Functional epistasis experiments in transgenic mice show that SPDEF amplifies ATOH1-dependent transcription but cannot independently initiate transcription of ATOH1 target genes.

CONCLUSIONS

This study unveils the direct targets of ATOH1 in the adult intestines and illuminates the transcriptional events that initiate the specification and function of intestinal secretory lineages.