Question: How to compare peaks (called by MACS2) within a specific genomic area for two samples?
gravatar for Amlan
3.4 years ago by
United States
Amlan20 wrote:


I have two ATAC-Seq peak samples for the same subject. I want to quantify the ATAC-Seq peaks within a defined area around TSS for each of those samples. Then I want to compare the scores to find out which of the two samples has more open chromatin region within that area around TSS?

So, is there any way to find a score for all peaks within a specific genomic area?

Thanks in advance

chip-seq • 2.1k views
ADD COMMENTlink modified 3.4 years ago by Alex Reynolds31k • written 3.4 years ago by Amlan20
gravatar for Alex Reynolds
3.4 years ago by
Alex Reynolds31k
Seattle, WA USA
Alex Reynolds31k wrote:

Given a set of intervals and a set of peaks with some score value, you could use BEDOPS bedmap --max-element or bedmap --max to get the maximum scoring peak over the interval, or just the maximum score, respectively:

$ bedmap --echo --max-element --delim '\t' intervals.bed peaks.bed > answer.bed

If you're looking for measured size or other peak features, there are other operators like --echo-map-size, --count, etc. that offer mapping of those results. Feel free to clarify what you're after, or take a look at the documentation or bedmap --help.

ADD COMMENTlink modified 3.4 years ago • written 3.4 years ago by Alex Reynolds31k
gravatar for YaGalbi
3.4 years ago by
Biocomputing, MRC Harwell Institute, Oxford, UK
YaGalbi1.5k wrote:

Im working on ATAC-seq data at the moment.

I've used ChipPeakAnno to annotate the peaks in the narrowPeaks file (MACS2 output) for genomic features (TSS etc). You could then do a simple count of the number of peaks annotated for each feature.

Have you considered using DESEQ2 for getting differential peaks first? That is a far more sophisticated package for comparing samples.

ADD COMMENTlink written 3.4 years ago by YaGalbi1.5k
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: 880 users visited in the last hour