Question: Summit Extension in ATAC-seq
gravatar for shankasal
2.3 years 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.6k views
ADD COMMENTlink modified 2.3 years ago by igor11k • written 2.3 years ago by shankasal10

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

ADD REPLYlink written 2.3 years ago by genomax87k

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 2.3 years ago by shankasal10
gravatar for igor
2.3 years ago by
United States
igor11k 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 2.3 years ago by igor11k
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: 1385 users visited in the last hour