Question: Does differential binding analysis across different ChIP-Seq samples makes sense ?
gravatar for Ar
4.1 years ago by
United States
Ar990 wrote:

Hello All,

I am interested to know if the differential binding analysis across different biological conditions makes sense or not. For example, if there are 2 different samples say A and B, where A is IPed using antibody P1 and B is IPed using antibody P2. and have their respective inputs. One does the peak calling using the inputs and get the set of peak for both the samples. Now, you find that a peak at certain genomic loci is present in both A and B; however, looking at the IGV plot you find that the peak in A has more reads than peak in B (considering that all the samples are normalized). When you do the differential analysis you find that the enrichment of the peak is significant in A more than B. Based on the above scenario, I have few questions:

1) I would like to know if such differential analysis makes sense or not ? 2) How can a peak be quantified same as the peak in RNA-Seq ? In RNA-Seq, you quantify the peaks because biologically, one knows there be can multiple copies of same mRNA

Ideally, I think a peak that is present should be considered as 1 and 0 if not. Please correct me if I wrong.

chip-seq • 1.5k views
ADD COMMENTlink modified 4.1 years ago by Devon Ryan96k • written 4.1 years ago by Ar990
gravatar for Devon Ryan
4.1 years ago by
Devon Ryan96k
Freiburg, Germany
Devon Ryan96k wrote:

While possible, such an analysis shouldn't be performed naively. For example, if IPs with P1 are more efficient than those with P2 then one would expect all P1 peaks to be higher...but in a biologically meaningless way. Assuming you can get around that issue, just use diffbind for the differential binding analysis, it shows in its tutorial how to quantify peaks.

BTW, you should not be comparing presence/absence of peaks using a binary label, that's like saying a p-value of 0.49999999 and 0.5 are significantly different.

ADD COMMENTlink written 4.1 years ago by Devon Ryan96k

Thanks that was helpful!

ADD REPLYlink written 4.1 years ago by Ar990

This post is rather old, but I stumbled across it and thought it was an interesting question. You can actually deal with the issue of different IP efficiencies if you use spike-in chromatin from another organism that's somewhat closely related. This paper describes the method and we're currently using it in my lab.

ADD REPLYlink written 2.3 years ago by goodez480
Please log in to add an answer.


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