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

Hello,

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 • 1.9k views
ADD COMMENTlink modified 3.0 years ago by Alex Reynolds30k • written 3.0 years ago by Amlan20
0
gravatar for Alex Reynolds
3.0 years ago by
Alex Reynolds30k
Seattle, WA USA
Alex Reynolds30k 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.0 years ago • written 3.0 years ago by Alex Reynolds30k
0
gravatar for YaGalbi
3.0 years ago by
YaGalbi1.5k
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.0 years ago by YaGalbi1.5k
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: 1275 users visited in the last hour