Question: MACS2 Peak Calling : Different number of peaks for each replicate
2
gravatar for Apoorva
4.8 years ago by
Apoorva260
United States
Apoorva260 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 ?

Thanks

ADD COMMENTlink modified 4.4 years ago by Biostar ♦♦ 20 • written 4.8 years ago by Apoorva260
4
gravatar for Ian
4.8 years ago by
Ian5.7k
University of Manchester, UK
Ian5.7k 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 modified 10 months ago by _r_am31k • written 4.8 years ago by Ian5.7k
1
gravatar for igor
4.8 years ago by
igor11k
United States
igor11k 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 4.8 years ago by igor11k
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: 1065 users visited in the last hour