peaks overlap: ATAC & CHIP seq
3
0
Entering edit mode
7.9 years ago

Hi guys. I have two peak files. One is from ATAC seq and the other is from H3K27ac chip seq. I am trying to merge these two peaks files together and get a file with overlapping peaks. There are several questions that I have here:

  1. what scripts or programs are good for that application? (bedops and bedmap?)

  2. I used MACS2 to call peaks: regular peaks calling for ATAC seq and broad peaks calling for H3K27ac chip seq. I am concerned whether the different peak calling methods will influence the overlapping results. Thank you!

ChIP-Seq • 7.0k views
ADD COMMENT
3
Entering edit mode
7.9 years ago

You can use bedops to do the merging and bedmap to get a list of peaks that overlap the merged regions:

$ bedops --merge A.bed B.bed C.bed ... > merge.bed
$ bedops --everything A.bed B.bed C.bed ... > union.bed
$ bedmap --echo --echo-map merge.bed union.bed > answer.bed

Or as a one-liner:

$ bedmap --echo --echo-map <(bedops --merge A.bed B.bed C.bed ...) <(bedops --everything A.bed B.bed C.bed ...) > answer.bed
ADD COMMENT
0
Entering edit mode

so can i combine peaks from different condition as well? or only replicates

ADD REPLY
1
Entering edit mode
7.9 years ago
Satyajeet Khare ★ 1.6k
  1. You can use bedtools intersect
  2. It should not matter much.
ADD COMMENT
0
Entering edit mode
7.9 years ago

You may want to consider using Fseq to call ATAC seq peaks instead of MACS2. H3K27ac gives fairly punctate peaks, so you should probably be fine, but any time you are doing intersections you should consider the size of the peak regions from both sets, and decide if how you want to define the peak boundaries for the intersected set.

ADD COMMENT

Login before adding your answer.

Traffic: 1179 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6