[ChIP-seq. signal plot] How to merge ChIP-seq. biological replicates for signal profile plotting?
1
0
Entering edit mode
6.3 years ago
chiefcat ▴ 180

Dear all,

In ChIP-seq. analysis, what is the proper way to combine biological replicates for signal profile plotting?

At the moment, I'm using deepTools plotProfile to visualize the signal density around TSS. After plotting the replicates separately, I can see that they are highly similar. Then, I would like to pool the replicates and presenting them as one single line in the profile. Should I achieve this by merging the bigWig files (e.g. bigWigMerge from UCSC tools as suggested by Devon in this post?) generated from bamCompare? (Each of the replicates has been compared to their own input in bamCompare & normalized to 1x Reads Per Genomic Content (RPGC).

Thanks!

deepTools ChIP-Seq biological replicates • 3.5k views
ADD COMMENT
1
Entering edit mode
6.3 years ago

I'd normally suggest keeping them separate. While you can merge them (and then scale accordingly, since bigWigMerge gives you an easy-to-manipulate text file), I'm not of the opinion that you're actually gaining much by doing so. Keeping samples separate allows people to easily judge biological variability, which is what they really want to get out of figures.

ADD COMMENT
0
Entering edit mode

Hi Devon, thanks for your reply. Should I scale them according to no. of reads?

ADD REPLY
1
Entering edit mode

They're already 1x normalized, so merge and divide by the number of replicates.

ADD REPLY

Login before adding your answer.

Traffic: 2926 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6