Question: Downsampling ChIP-seq BAM files with spike in normalization factor before feeding into Diffbind
gravatar for shikhtechai
16 months ago by
shikhtechai0 wrote:

I am planning to use exogenous chromatin as a spike in control with my actual sample from mouse to perform ChIP-seq for peak calling and differential binding analysis for histone modifications. this involves down-sampling the uniquely mapped read files to the calculated normalization factor from the spike in. This is generally helpful for peak calling and visualizing in IGV browser. But I have not found any reference on whether it is considered in differential binding analysis. Also, nothing is mentioned in the Diffbind vignette. Could you please explain how this strategy might affect using the diffbind package?

I greatly appreciate your time to read and answer to my question! Thank you in advance!


chipseq diffbind chipqc • 523 views
ADD COMMENTlink modified 14 months ago by Rory Stark2.8k • written 16 months ago by shikhtechai0
Answer: Downsampling ChIP-seq BAM files with spike in normalization factor before feedin
gravatar for Rory Stark
14 months ago by
Rory Stark2.8k
CRUK, Cambridge, UK
Rory Stark2.8k wrote:

In you downsampled the reads for visualization and peak calling, you should use the down-sampled files for differential binding analysis. Normalization is even more important for differential analysis than for peak calling etc. Be aware that DiffBind will still normalize the reads to the relative number of reads in each library, so this may negate some of what you are trying to do. You can bypass DiffBind's normalization by setting the library sizes to appear all the same, as discussed in this thread: MAnorm output as DiffBind input..

I'm not exactly which normalization procedure you are using. We have looked at some methods where we use the control reads to calculate normalization factors rather than actually downsampling reads (hate to throw away good data!). See for example this paper:

Novel Quantitative ChIP-seq Methods Measure Absolute Fold-Change in ER Binding Upon Fulvestrant Treatment

a version of which is to be published in NAR soon.

ADD COMMENTlink written 14 months ago by Rory Stark2.8k
Please log in to add an answer.


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