Question: MACS2 Peak Calling : Different number of peaks for each replicate
gravatar for ApoorvaB
3.9 years ago by
United States
ApoorvaB230 wrote:

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 ?


ADD COMMENTlink modified 3.4 years ago by Biostar ♦♦ 20 • written 3.9 years ago by ApoorvaB230
gravatar for Ian
3.9 years ago by
University of Manchester, UK
Ian5.6k wrote:

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 COMMENTlink written 3.9 years ago by Ian5.6k
gravatar for igor
3.9 years ago by
United States
igor8.9k wrote:

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 COMMENTlink written 3.9 years ago by igor8.9k
Please log in to add an answer.


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