Search
Question: Chip-seq analysis with input
0
2.3 years ago by
damian.kao0 wrote:

I read through the other posts about using RNA-seq packages to analyze Chip-seq data (https://support.bioconductor.org/p/72098/). It seems like the general consensus is to just ignore the input control samples or build a black-list and look at differential binding between IP samples.

If I do want to incorporate input control into my differential binding, is it valid to include that data as another factor in the design matrix and perform a difference of difference?

So for every library, I would have a "IP" factor with two levels (IP/Input) and also a "sample" factor with two levels for treatment and control.

I am not very well versed in R, would it be possible to make this kind of contrast? And would this type of contrast even be valid?

modified 2.3 years ago by Aaron Lun21k • written 2.3 years ago by damian.kao0
0
2.3 years ago by
Aaron Lun21k
Cambridge, United Kingdom
Aaron Lun21k wrote:

Can you do it? Yes, because GLMs are very flexible and can easily test for differences of differences.

Should you do it? That's harder to answer, but in most cases, probably not. See A: csaw with negative controls for an explanation of why looking for differences of differences tends to be counterproductive.

Thanks. I am trying to understand what you meant when you wrote:

However, this becomes problematic when changes in chromatin state (that affect input coverage) coincide with changes in binding....

Are you saying that for a given genomic region that has signal in the input control, that signal is composed of various factors (mappability, gc content, open chromatin...). One of the factors of input control signal could also be the protein of interest that we are IP-ing. To then try to determine a difference between the protein of interest and the input control (which includes signal from the protein of interest) would not be productive. You would essentially be equalizing out the sample vs input signal.