Question: Quantification And Comparison Of Chip-Seq Peaks
3
gravatar for Ying W
8.5 years ago by
Ying W4.0k
South San Francisco, CA
Ying W4.0k wrote:

I have two chip-seq experiments under two conditions. I was wondering if it was possible to say that a peak is 'bigger' in one condition versus the other and if there is software that will quantify the difference. Current chip-seq peak callers seem to just say if there is a peak or not.

Thanks

differential chip-seq • 8.2k views
ADD COMMENTlink modified 6.1 years ago by dnaseiseq200 • written 8.5 years ago by Ying W4.0k
3
gravatar for 2184687-1231-83-
8.5 years ago by
2184687-1231-83-5.0k wrote:

One option is to use DESeq to analyse count data from the Chip-seq summits in your two conditions and test for differential intensity:

http://www-huber.embl.de/users/anders/DESeq/

ADD COMMENTlink modified 4 months ago by RamRS25k • written 8.5 years ago by 2184687-1231-83-5.0k
1

@avilella +1 I have been using DESeq for ChIP-seq data (comparing read counts over two replicates of two time points [4 samples]). Apparently there is a GLM model that will allow the input control read counts to be factored in... Have you looked at this?

ADD REPLYlink written 8.5 years ago by Ian5.6k
1
gravatar for Ian
8.5 years ago by
Ian5.6k
University of Manchester, UK
Ian5.6k wrote:

You could also try DIME or EdgeR.

I am currently trying DESeq, which seems to be producing sensible results. I am going to try DIME next as a comparison.

ADD COMMENTlink written 8.5 years ago by Ian5.6k
0
gravatar for Ying W
8.5 years ago by
Ying W4.0k
South San Francisco, CA
Ying W4.0k wrote:

Thanks, I will give DESeq a try,

I have found a paper that seems to address this issue but approaches it in an entirely different way

"Identifying dispersed epigenomic domains from ChIP-Seq data." http://www.ncbi.nlm.nih.gov/pubmed/21325299

Their software RSEG can be found here: http://smithlab.cmb.usc.edu/histone/rseg/

ADD COMMENTlink written 8.5 years ago by Ying W4.0k
0
gravatar for Ziegfeld
8.5 years ago by
Ziegfeld0
US
Ziegfeld0 wrote:

Also SICER. Rseg runs slower than SICER. You may want to use ChIPDiff which uses HMM and could be fed the output -island.bed files. You can tune confidence% and "fold change" parameters to solve your problem.

ADD COMMENTlink written 8.5 years ago by Ziegfeld0
0
gravatar for dnaseiseq
6.1 years ago by
dnaseiseq200
United Kingdom
dnaseiseq200 wrote:

Bailey T, Krajewski P, Ladunga I, Lefebvre C, Li Q, et al. (2013) Practical Guidelines for the Comprehensive Analysis of ChIP-seq Data. PLoS Comput Biol 9(11): e1003326. doi:10.1371/journal.pcbi.1003326

http://www.ploscompbiol.org/article/info:doi/10.1371/journal.pcbi.1003326

ADD COMMENTlink written 6.1 years ago by dnaseiseq200
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1702 users visited in the last hour