MACS2 Peak Calling : Different number of peaks for each replicate
2
2
Entering edit mode
6.0 years ago
Apoorva ▴ 330

Hi All,

I am a newbie and I am doing a ChIP-Seq analysis on a small dataset of 6 samples (Wild-type and knockout having 3 replicates each). I ran MACS2 for peak calling  but the number of peaks in each replicate is very different. For example, WT-Rep1 has 1500 peaks and WT-Rep2 has only 34 peaks. Why is there such a variation ? Did I do something wrong or is the data bad ?

Thanks

ChIP-Seq MACS2 biological replicates • 3.6k views
ADD COMMENT
4
Entering edit mode
6.0 years ago
Ian 5.8k

Given that you have three replicates for each condition you might want to consider differential binding analysis, e.g. diffReps. It will use the replicates to estimate variance and return windows of the genome that have significantly more binding in one sample compared to the other. The input is the mapped reads, and does not require called peaks, which gets round calling peaks based on arbitrary thresholds.

ADD COMMENT
1
Entering edit mode
6.0 years ago
igor 12k

Those are actually not very different, but both are a little low. There can be a number of reasons why they like that. You can start by looking at the total number of reads, alignment rate, and duplication rate. That may explain why one is better than the other. You can also look at some of the actual peaks in a genome browser and see how the different replicates compare.

ADD COMMENT

Login before adding your answer.

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