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

  1. Peter N. Robinson1,2,3,6
  1. 1Institute for Medical and Human Genetics, Charité–Universitätsmedizin Berlin, 13353 Berlin, Germany;
  2. 2Berlin Brandenburg Center for Regenerative Therapies (BCRT), Charité–Universitätsmedizin Berlin, 13353 Berlin, Germany;
  3. 3Max Planck Institute for Molecular Genetics, 14195 Berlin, Germany;
  4. 4Department of Biostatistics, Clinical Research Unit, Berlin Institute of Health, Charité–Universitätsmedizin Berlin, 13353 Berlin, Germany;
  5. 5Labor für Pädiatrische Molekularbiologie, Charité–Universitätsmedizin Berlin, 10117, Berlin, Germany;
  6. 6Institute for Bioinformatics, Department of Mathematics and Computer Science, Freie Universität Berlin, 14195 Berlin, Germany
  1. Corresponding author: peter.robinson{at}charite.de

Abstract

Chromatin immunoprecipitation coupled with next-generation sequencing (ChIP-seq) is a powerful technology to identify the genome-wide locations of transcription factors and other DNA binding proteins. Computational ChIP-seq peak calling infers the location of protein–DNA interactions based on various measures of enrichment of sequence reads. In this work, we introduce an algorithm, Q, that uses an assessment of the quadratic enrichment of reads to center candidate peaks followed by statistical analysis of saturation of candidate peaks by 5′ ends of reads. We show that our method not only is substantially faster than several competing methods but also demonstrates statistically significant advantages with respect to reproducibility of results and in its ability to identify peaks with reproducible binding site motifs. We show that Q has superior performance in the delineation of double RNAPII and H3K4me3 peaks surrounding transcription start sites related to a better ability to resolve individual peaks. The method is implemented in C+l+ and is freely available under an open source license.

Footnotes

  • Received January 21, 2015.
  • Accepted July 6, 2015.

This article is distributed exclusively by Cold Spring Harbor Laboratory Press for the first six months after the full-issue publication date (see http://genome.cshlp.org/site/misc/terms.xhtml). After six months, it is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), as described at http://creativecommons.org/licenses/by-nc/4.0/.

| Table of Contents

Preprint Server