Question: Summit Extension in ATAC-seq
gravatar for shankasal
18 months ago by
shankasal10 wrote:

Something I've seen with ATAC-seq analysis is the extension of peak summits (ie. "Peak summits were extended +/- 250 bp and filtered for nonoverlapping, maximally significant 500-bp peaks"). Even DiffBind does the extension with the summits parameter in dba.counts.

I was hoping someone had a good explanation as to why extended summit peaks are used instead of the actual peak sizes called by the peak caller (MACS2 in my case). Moreover, what is the justification for choosing a specific extension size. I've seen 250 bp used frequently, but I've also seen smaller extensions of 100 bp.

Lastly, I'm not sure of the best way to extend peaks summits, as in, which tool to use and whether to use the summits file output from MACS2 or the summits position in the narrowPeaks file.

sequencing atac-seq • 1.1k views
ADD COMMENTlink modified 18 months ago by igor8.6k • written 18 months ago by shankasal10

Previous thread from @igor may be of interest: ATAC-seq peak calling with MACS

ADD REPLYlink written 18 months ago by genomax73k

Thanks, I've had a look at that before to decide on the shift and extsize I wanted to use. At this point though I'm looking more at generating a non-redundant peakset for count comparison which is usually where I see the summit extension implemented.

ADD REPLYlink written 18 months ago by shankasal10
gravatar for igor
18 months ago by
United States
igor8.6k wrote:

I noticed that as well. I assume it is to keep all the peaks the same size and to ignore the potential nucleosome laddering issues.

ADD COMMENTlink written 18 months ago by igor8.6k
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: 990 users visited in the last hour