Question: Summit Extension in ATAC-seq
1
gravatar for shankasal
12 months ago by
shankasal10
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 • 832 views
ADD COMMENTlink modified 12 months ago by igor7.6k • written 12 months ago by shankasal10
1

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

ADD REPLYlink written 12 months ago by genomax65k

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 12 months ago by shankasal10
0
gravatar for igor
12 months ago by
igor7.6k
United States
igor7.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 12 months ago by igor7.6k
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: 1689 users visited in the last hour