WACS: improving ChIP-seq peak calling by optimally weighting controls
Open Access
- 15 February 2021
- journal article
- research article
- Published by Springer Science and Business Media LLC in BMC Bioinformatics
- Vol. 22 (1), 1-21
- https://doi.org/10.1186/s12859-020-03927-2
Abstract
Chromatin immunoprecipitation followed by high throughput sequencing (ChIP-seq), initially introduced more than a decade ago, is widely used by the scientific community to detect protein/DNA binding and histone modifications across the genome. Every experiment is prone to noise and bias, and ChIP-seq experiments are no exception. To alleviate bias, the incorporation of control datasets in ChIP-seq analysis is an essential step. The controls are used to account for the background signal, while the remainder of the ChIP-seq signal captures true binding or histone modification. However, a recurrent issue is different types of bias in different ChIP-seq experiments. Depending on which controls are used, different aspects of ChIP-seq bias are better or worse accounted for, and peak calling can produce different results for the same ChIP-seq experiment. Consequently, generating “smart” controls, which model the non-signal effect for a specific ChIP-seq experiment, could enhance contrast and increase the reliability and reproducibility of the results. We propose a peak calling algorithm, Weighted Analysis of ChIP-seq (WACS), which is an extension of the well-known peak caller MACS2. There are two main steps in WACS: First, weights are estimated for each control using non-negative least squares regression. The goal is to customize controls to model the noise distribution for each ChIP-seq experiment. This is then followed by peak calling. We demonstrate that WACS significantly outperforms MACS2 and AIControl, another recent algorithm for generating smart controls, in the detection of enriched regions along the genome, in terms of motif enrichment and reproducibility analyses. This ultimately improves our understanding of ChIP-seq controls and their biases, and shows that WACS results in a better approximation of the noise distribution in controls.Keywords
Funding Information
- NSERC Discovery Grant (328154-2014)
- Queen Elizabeth II Graduate Scholarship in Science & Technology
- Compute Canada Resources-for-Research-Groups
This publication has 33 references indexed in Scilit:
- An integrated encyclopedia of DNA elements in the human genomeNature, 2012
- Using MACS to Identify Peaks from ChIP‐Seq DataCurrent Protocols in Bioinformatics, 2011
- FIMO: scanning for occurrences of a given motifBioinformatics, 2011
- BEDTools: a flexible suite of utilities for comparing genomic featuresBioinformatics, 2010
- Computation for ChIP-seq and RNA-seq studiesNature Methods, 2009
- A clustering approach for identification of enriched domains from histone modification ChIP-Seq dataBioinformatics, 2009
- PeakSeq enables systematic scoring of ChIP-seq experiments relative to controlsNature Biotechnology, 2009
- FindPeaks 3.1: a tool for identifying areas of enrichment from massively parallel short-read sequencing technologyBioinformatics, 2008
- Mapping and quantifying mammalian transcriptomes by RNA-SeqNature Methods, 2008
- High-Resolution Profiling of Histone Methylations in the Human GenomeCell, 2007