Comparison of Four ChIP-Seq Analytical Algorithms Using Rice Endosperm H3K27 Trimethylation Profiling Data
Open Access
- 30 September 2011
- journal article
- research article
- Published by Public Library of Science (PLoS) in PLOS ONE
- Vol. 6 (9), e25260
- https://doi.org/10.1371/journal.pone.0025260
Abstract
Chromatin immunoprecipitation coupled with high throughput DNA Sequencing (ChIP-Seq) has emerged as a powerful tool for genome wide profiling of the binding sites of proteins associated with DNA such as histones and transcription factors. However, no peak calling program has gained consensus acceptance by the scientific community as the preferred tool for ChIP-Seq data analysis. Analyzing the large data sets generated by ChIP-Seq studies remains highly challenging for most molecular biology laboratories. Here we profile H3K27me3 enrichment sites in rice young endosperm using the ChIP-Seq approach and analyze the data using four peak calling algorithms (FindPeaks, PeakSeq, USeq, and MACS). Comparison of the four algorithms reveals that these programs produce very different peaks in terms of peak size, number, and position relative to genes. We verify the peak predictions using ChIP-PCR to evaluate the accuracy of peak prediction of the four algorithms. We discuss the approach of each algorithm and compare similarities and differences in the results. Despite their differences in the peaks identified, all of the programs reach similar conclusions about the effect of H3K27me3 on gene expression. Its presence either upstream or downstream of a gene is predominately associated with repression of the gene. Additionally, GO analysis finds that a substantially higher ratio of genes associated with H3K27me3 were involved in multicellular organism development, signal transduction, response to external and endogenous stimuli, and secondary metabolic pathways than the rest of the rice genome.This publication has 31 references indexed in Scilit:
- Computation for ChIP-seq and RNA-seq studiesNature Methods, 2009
- PeakSeq enables systematic scoring of ChIP-seq experiments relative to controlsNature Biotechnology, 2009
- Design and analysis of ChIP-seq experiments for DNA-binding proteinsNature Biotechnology, 2008
- Endosperm gene imprinting and seed developmentCurrent Opinion in Genetics & Development, 2007
- Polycomb group and trithorax group proteins in ArabidopsisBiochimica et Biophysica Acta (BBA) - Gene Structure and Expression, 2007
- High-Resolution Profiling of Histone Methylations in the Human GenomeCell, 2007
- Different Polycomb group complexes regulate common target genes in ArabidopsisEMBO Reports, 2006
- Control of Developmental Regulators by Polycomb in Human Embryonic Stem CellsCell, 2006
- Polycomb complexes repress developmental regulators in murine embryonic stem cellsNature, 2006
- The map-based sequence of the rice genomeNature, 2005