×

Spatio-temporal model for multiple ChIP-seq experiments. (English) Zbl 1311.62165

Summary: The increasing availability of ChIP-seq data demands for advanced statistical tools to analyze the results of such experiments. The inherent features of high-throughput sequencing output call for a modelling framework that can account for the spatial dependency between neighboring regions of the genome and the temporal dimension that arises from observing the protein binding process at progressing time points; also, multiple biological/technical replicates of the experiment are usually produced and methods to jointly account for them are needed. Furthermore, the antibodies used in the experiment lead to potentially different immunoprecipitation efficiencies, which can affect the capability of distinguishing between the true signal in the data and the background noise. The statistical procedure proposed consist of a discrete mixture model with an underlying latent Markov random field: the novelty of the model is to allow both spatial and temporal dependency to play a role in determining the latent state of genomic regions involved in the protein binding process, while combining all the information of the replicates available instead of treating them separately. It is also possible to take into account the different antibodies used, in order to obtain better insights of the process and exploit all the biological information available.

MSC:

62M40 Random fields; image analysis
92B15 General biostatistics

Software:

COUNT; HPeak; BayesPeak
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Bao, Y., V. Vinciotti, E. Wit and P. T. Hoen (2014): “Joint modelling of ChIP-seq data via a Markov Random Field model,” Biostatistics, 15(2), 296-310.;
[2] Bardet, A. F., Q. He, J. Zeitlinger and A. Stark (2012): “A computational pipeline for comparative ChIP-seq analyses,” Nature Protocols, 7(1), 45-61.;
[3] Benjamini, Yoav and Yosef Hochberg (1995): “Controlling the false discovery rate: a practical and powerful approach to multiple testing,” Journal of the Royal Statistical Society. Series B (Methodological), 57(1), 289-300.; · Zbl 0809.62014
[4] Hilbe, J. M. (2011): Negative binomial regression, Cambridge University Press, Cambridge, England, UK.; · Zbl 1269.62063
[5] Kharchenko, Peter V., Michael Y. Tolstorukov and Michael Y. Park (2008): “Design and analysis of ChIP-seq experiments for DNA-binding proteins.” Nature Biotechnology, 26, 1351-1359.;
[6] Kindermann, Ross and J. Laurie Snell (1980): “Markov Random Fields and Their Applications,” American Mathematical Society, Providence, USA.; · Zbl 1229.60003
[7] Kuan, Pei Fen, Dongjun Chung, Guangjin Pan, James A. Thomson, Ron Stewart and Sündüz Keleş (2001): “A statistical framework for the analysis of ChIP-Seq data,” Journal of the American Statistical Association, 106(495), 891-903.; · Zbl 1229.62146
[8] Lauritzen, Steffen L. (1996): Graphical models. Oxford University Press, Oxford, England, UK.; · Zbl 0907.62001
[9] Qin, Zhaohui S., Jianjun Yu, Jincheng Shen, Christopher A. Maher, Ming Hu, Shanker Kalyana-Sundaram, Jindan Yu and Arul M. Chinnaiyan (2010): “HPeak: an HMM-based algorithm for defining read-enriched regions in ChIP-Seq data,” BMC Bioinformatics, 11(1), 369.;
[10] Ramos, Yolande F. M., Matthew S. Hestand, Matty Verlaan, Elise Krabbendam, Yavuz Ariyurek, Michiel van Galen, Hans van Dam, Gert-Jan B. van Ommen, Johan T. den Dunnen, Alt Zantema and Peter A. C. ’t Hoen (2010): “Genome-wide assessment of differential roles for p300 and CBP in transcription regulation,” Nucleic Acids Research, 38(16), 5396-5408.;
[11] Spyrou, C., R. Stark, A. G. Lynch and S. Tavar (2009): “BayesPeak: Bayesian analysis of ChIP-seq data,” BMC Bioinformatics, 10(1), 299.;
[12] Zeng, Xin, Rajendran Sanalkumar, Emery H. Bresnick, Hongda Li, Qiang Chang and Sündüz Keles (2013): “jMOSAiCS: joint analysis of multiple ChIP-seq datasets,” Genome Biology, 14(4), R38.;
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.