Question: Normalization using spike-in chromatin
gravatar for goodez
8 days ago by
Kansas City, MO
goodez160 wrote:

I have many human ChIP-seq samples which I need to compare. All of them received the same 10% spike-in of mouse chromatin. At first I was going to simply normalize each sample to the total number of spike-in reads. However I've been asked to do something a little more sophisticated now and need some advice.

I have used a method called Irreproducible Discovery Rate (IDR) to find ChIP-seq peaks (called by MACS2) that are very reproducible and high confidence. I need to only use these confident peak regions for my normalization. As in, I have counted up total mouse spike-in reads that fall into the confident peak regions for each sample. So for example for one antibody I have this scenario:

Sample                     Reads within confident peaks                         
Wildtype rep. 1                 1.02 million
Wildtype rep. 2                 0.78 million

Mutant rep. 1                     1.01 million
Mutant rep. 2                     0.60 million

What is the best way to normalize the samples in this case? And should I still be normalizing to RPM prior to the spike-in scaling factor? I am thinking probably I should skip the RPM step and just use the spike-in. Hopefully this question makes sense to someone. Thanks

ADD COMMENTlink written 8 days ago by goodez160
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: 1317 users visited in the last hour