Display options
Share it on

Front Genet. 2014 Nov 14;5:399. doi: 10.3389/fgene.2014.00399. eCollection 2014.

Quantifying the impact of inter-site heterogeneity on the distribution of ChIP-seq data.

Frontiers in genetics

Jonathan Cairns, Andy G Lynch, Simon Tavaré

Affiliations

  1. Nuclear Dynamics Group, The Babraham Institute Cambridge, UK ; Cancer Research UK Cambridge Institute, University of Cambridge Cambridge, UK.
  2. Cancer Research UK Cambridge Institute, University of Cambridge Cambridge, UK.

PMID: 25452765 PMCID: PMC4231950 DOI: 10.3389/fgene.2014.00399

Abstract

Chromatin Immunoprecipitation followed by sequencing (ChIP-seq) is a valuable tool for epigenetic studies. Analysis of the data arising from ChIP-seq experiments often requires implicit or explicit statistical modeling of the read counts. The simple Poisson model is attractive, but does not provide a good fit to observed ChIP-seq data. Researchers therefore often either extend to a more general model (e.g., the Negative Binomial), and/or exclude regions of the genome that do not conform to the model. Since many modeling strategies employed for ChIP-seq data reduce to fitting a mixture of Poisson distributions, we explore the problem of inferring the optimal mixing distribution. We apply the Constrained Newton Method (CNM), which suggests the Negative Binomial - Negative Binomial (NB-NB) mixture model as a candidate for modeling ChIP-seq data. We illustrate fitting the NB-NB model with an accelerated EM algorithm on four data sets from three species. Zero-inflated models have been suggested as an approach to improve model fit for ChIP-seq data. We show that the NB-NB mixture model requires no zero-inflation and suggest that in some cases the need for zero inflation is driven by the model's inability to cope with both artifactual large read counts and the frequently observed very low read counts. We see that the CNM-based approach is a useful diagnostic for the assessment of model fit and inference in ChIP-seq data and beyond. Use of the suggested NB-NB mixture model will be of value not only when calling peaks or otherwise modeling ChIP-seq data, but also when simulating data or constructing blacklists de novo.

Keywords: ChIP-seq; Negative Binomial; high-throughput sequencing; mixture model; zero-inflation

References

  1. Biometrics. 2005 Mar;61(1):179-85 - PubMed
  2. BMC Bioinformatics. 2006 Mar 20;7:157 - PubMed
  3. Genome Res. 2008 Sep;18(9):1509-17 - PubMed
  4. PLoS Comput Biol. 2008 Aug 22;4(8):e1000158 - PubMed
  5. Genome Biol. 2008;9(9):R137 - PubMed
  6. Nat Rev Genet. 2009 Oct;10(10):669-80 - PubMed
  7. BMC Bioinformatics. 2009 Sep 21;10:299 - PubMed
  8. Science. 2010 May 21;328(5981):1036-40 - PubMed
  9. Theor Biol Med Model. 2010 Jun 03;7:18 - PubMed
  10. Bioinformatics. 2011 Mar 1;27(5):713-4 - PubMed
  11. Bioinformatics. 2011 Mar 15;27(6):870-1 - PubMed
  12. PLoS Biol. 2011 Apr;9(4):e1001046 - PubMed
  13. Genome Biol. 2011 Jul 25;12(7):R67 - PubMed
  14. Nature. 2012 Jan 04;481(7381):389-93 - PubMed
  15. Stat Appl Genet Mol Biol. 2012 Mar 31;11(3):Article 9 - PubMed
  16. Genome Res. 2012 Dec;22(12):2489-96 - PubMed
  17. Genome Res. 2012 Sep;22(9):1813-31 - PubMed
  18. Stat Med. 2014 Jun 15;33(13):2238-50 - PubMed
  19. Front Genet. 2014 Apr 10;5:75 - PubMed

Publication Types