Question: Chip-seq on pooled replicates
0
gravatar for RT
2.4 years ago by
RT340
European Union
RT340 wrote:

Dear All,

Can someone please help me to understand why for ChIP-seq we should not merge the biological replicates before peak calling. I read in several papers, peaks calling is done on individual replicates. I am struggling to understand this. Ideally if I merge replicates I should have more power, means more coverage in peak regions.

My results also indicates that if I call peaks on individual replicates, I get around 10,000 peaks in each rep. I have 3 replicates for chip and 3 for control. If I merge these together and run macs2, I get very few peaks < 100. Can anyone please explain this.

Thanks. RT

pooled replicates chip-seq • 1.3k views
ADD COMMENTlink modified 2.4 years ago by EagleEye6.4k • written 2.4 years ago by RT340

I think sometimes replicates do just get pooled so it's not out of the question, and many fewer programs exist for properly handling chip-seq replicates separately. the IDR pipeline from encode tried to tackle the problem ("Irreproducibility Discovery Rate") but it is probably just good to manually inspect the data like @EagleEye recommends to see if there are problems with the enrichment or something.

ADD REPLYlink written 2.4 years ago by cmdcolin1.2k
0
gravatar for EagleEye
2.4 years ago by
EagleEye6.4k
Sweden
EagleEye6.4k wrote:

Before merging the replicates please check how the replicate samples/experiments correlated (example, PCA, deeptools etc.,).

Compare replicates to check their correlation

ADD COMMENTlink modified 2.4 years ago • written 2.4 years ago by EagleEye6.4k

Yeah it seems that the noise in a sample could be adding to this.

ADD REPLYlink written 2.4 years ago by datascientist28390
Please log in to add an answer.

Help
Access

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