Comparison of Four ChIP-Seq Analytical Algorithms Using Rice Endosperm H3K27 Trimethylation Profiling Data
Author Information
Author(s): Brandon M. Malone, Feng Tan, Susan M. Bridges, Zhaohua Peng
Primary Institution: Mississippi State University
Hypothesis
How do different peak calling algorithms compare in identifying H3K27me3 enrichment sites in rice endosperm?
Conclusion
All four peak calling algorithms produced different peaks, but they reached similar conclusions about the effect of H3K27me3 on gene expression.
Supporting Evidence
- Different peak calling algorithms produced varying peak sizes, numbers, and positions relative to genes.
- Despite differences in peaks, all programs indicated that H3K27me3 is associated with gene repression.
- GO analysis showed that genes associated with H3K27me3 were involved in multicellular organism development and response to stimuli.
Takeaway
Scientists used four different computer programs to find important spots in rice DNA that help control how genes are turned on or off, and they found that while the programs showed different results, they all agreed on the main idea.
Methodology
The study used ChIP-Seq to profile H3K27me3 enrichment sites in rice endosperm and compared four peak calling algorithms: FindPeaks, PeakSeq, USeq, and MACS.
Potential Biases
Different algorithms may introduce biases in peak identification due to their varying methodologies.
Limitations
The algorithms produced varying results, and some peaks identified by PeakSeq were likely false positives.
Participant Demographics
Rice strain Oryza sativa ssp japonica cv Nipponbare was used in the study.
Statistical Information
P-Value
p<2.2e-16
Statistical Significance
p<0.05
Digital Object Identifier (DOI)
Want to read the original?
Access the complete publication on the publisher's website