Pooling replicates for CHIP-seq peak calling
0
0
Entering edit mode
4.0 years ago
Pappu ★ 2.1k

I am wondering what is the difference between pooling 2 replicates and 2 IGG controls for ChIP-seq peak calling compared to running the peak calling separately and doing IDR to integrate the peaks.

ChIP-Seq • 2.2k views
ADD COMMENT
1
Entering edit mode

IDR will give you those peaks somewhat consistent between replicates (and by this more reliable, probably preferred when looking for "true" binding locations) while merging just "takes it all" without the chance to make statements about reproducibility.

ADD REPLY
0
Entering edit mode

I agree. Maybe if you were having problems getting peaks due to being borderline on your adequate sequencing depth of the samples, or bad quality, etc. you could consider pooling the replicates to explore your results...

ADD REPLY
1
Entering edit mode

On top of the per-sample quality, what people often forget is that IDR requires that peaks are being called with relaxed settings, e.g. a p-value threshold (not q-value) of like 0.05 so very non-stringent. Reason is that the model requires both signal and noise and calling to stringently will remove much noise, making the model inadequate. You also have to sort data by the ranking metric, e.g. the -log10(p) field. This tutorial summarizes the essentials: https://hbctraining.github.io/Intro-to-ChIPseq/lessons/07_handling-replicates-idr.html

ADD REPLY
0
Entering edit mode

I agree. I saw pooling of 2 control replicates for peak calling in MACS2.

ADD REPLY

Login before adding your answer.

Traffic: 2292 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