Abstract:Background RNA polymerase II (PolII) is essential in gene transcription and ChIP-seq experiments have been used to study PolII binding patterns over the entire genome. However, since PolII enriched regions in the genome can be very long, existing peak finding algorithms for ChIP-seq data are not adequate for identifying such long regions. Methods Here we propose an enriched region detection method for ChIP-seq data to identify long enriched regions by combining a signal denoising algorithm with a false discovery rate (FDR) approach. The binned ChIP-seq data for PolII are first processed using a non-local means (NL-means) algorithm for purposes of denoising. Then, a FDR approach is developed to determine the threshold for marking enriched regions in the binned histogram. Results We first test our method using a public PolII ChIP-seq dataset and compare our results with published results obtained using the published algorithm HPeak. Our results show a high consistency with the published results (80-100%). Then, we apply our proposed method on PolII ChIP-seq data generated in our own study on the effects of hormone on the breast cancer cell line MCF7. The results demonstrate that our method can effectively identify long enriched regions in ChIP-seq datasets. Specifically, pertaining to MCF7 control samples we identified 5,911 segments with length of at least 4 Kbp (maximum 233,000 bp); and in MCF7 treated with E2 samples, we identified 6,200 such segments (maximum 325,000 bp). Conclusions We demonstrated the effectiveness of this method in studying binding patterns of PolII in cancer cells which enables further deep analysis in transcription regulation and epigenetics. Our method complements existing peak detection algorithms for ChIP-seq experiments.

Increased peak detection accuracy in over-dispersed ChIP-seq data with supervised segmentation models

Visual annotations and a supervised learning approach for evaluating and calibrating ChIP-seq peak detectors

Diffreps: Detecting Differential Chromatin Modification Sites from ChIP-seq Data with Biological Replicates.

Analyzing ChIP-seq Data Based on Multiple Knowledge Sources for Histone Modification.

Picking ChIP-seq peak detectors for analyzing chromatin modification experiments

Systematic Evaluation of Factors Influencing ChIP-seq Fidelity

Joint modelling of ChIP-seq data via a Markov random field model

Probabilistic Peak Calling And Controlling False Discovery Rate Estimations In Transcription Factor Binding Site Mapping From Chip-Seq

A signal processing approach for enriched region detection in RNA polymerase II ChIP-seq data

Identification of regions differentially modified by histone modification based on the ChIP-seq data

High-Throughput Computational Approaches to Analyzing Histone Modification Next-Generation Sequencing Data

CNN-Peaks: ChIP-Seq peak detection pipeline using convolutional neural networks that imitate human visual inspection

A Generalized Linear Model for Peak Calling in ChIP-Seq Data

Computational Methodology for ChIP-seq Analysis

Unified Analysis of Multiple ChIP-Seq Datasets

Knowledge Integration for Analyzing ChIP-seq

ML-Peaks: CHIP-Seq Peak Detection Pipeline using Machine Learning Techniques

Efficiently Identifying Genome-Wide Changes with Next-Generation Sequencing Data

DiffChIPL: a differential peak analysis method for high-throughput sequencing data with biological replicates based on limma

Qb-13006-xsl 54..70

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