Question: MACS2 with multiple replicates
gravatar for ariel.balter
3.3 years ago by
ariel.balter140 wrote:

When you use macs2 with multiple replicates is it just pooling the replicates or doing something more sophisticated?

chip-seq macs2 • 4.9k views
ADD COMMENTlink modified 2.2 years ago by bwd40 • written 3.3 years ago by ariel.balter140
gravatar for bwd
2.2 years ago by
bwd40 wrote:

Pooling replicates may not be the most appropriate method, e.g., it does not account for differences across replicates. Instead use either of the followings:

  • call peaks on each of the replicates separately using MACS, then process the replicates using MSPC, which keeps only reproducible peaks (combines p-values and then applies a threshold on their right-tail probability) and generates a set of "consensus" peaks.
  • use comparative peak callers such as jMOSAiCS, which is specifically developed to call peaks considering multiple replicates.
  • last option could be IDR, AFAIK it does not account for replicate type and does not consider strong peaks as "reproducible" if they are not overlapping with other peaks. This can lead to discarding true-binding sites who have significant p-value but not overlapping with peaks from other replicates; one can argue that such peaks reflect expected discrepancy between biological replicates (as opposed to technical replicates). Therefore, discarding them may not be the best practice. However, MSPC can "confirm" strong peaks in bio reps even if they do not overlap with any other peak (this is a configurable as you can adjust it with its -c parameter).
ADD COMMENTlink modified 8 weeks ago • written 2.2 years ago by bwd40
gravatar for Praneet Chaturvedi
3.3 years ago by
Cincinnati Children's Hospital and Medical Center
Praneet Chaturvedi110 wrote:

You can use ENCODE IDR pipeline with peak caller as MACS2 to call peaks which are reproducible across replicates. Here is the link :

Cheers !!


ADD COMMENTlink written 3.3 years ago by Praneet Chaturvedi110

I just want to know what method macs2 is using. Calling multiple replicates with macs2 is part of the encode pipeline.

ADD REPLYlink written 3.3 years ago by ariel.balter140
gravatar for James Ashmore
3.3 years ago by
James Ashmore3.0k
UK/Edinburgh/MRC Centre for Regenerative Medicine
James Ashmore3.0k wrote:

As far as I'm aware it just pools the replicates.

ADD COMMENTlink modified 3.3 years ago • written 3.3 years ago by James Ashmore3.0k

Yeah I think you are probably correct based on the output I got.

ADD REPLYlink written 3.3 years ago by ariel.balter140
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: 776 users visited in the last hour