Find consensus overlap between many bed files

Question: Find consensus overlap between many bed files

5 months ago by

goodez • 420

United States

goodez • 420 wrote:

I usually use bedtools intersect to find overlapping regions of bed files, but it seems like this tool can only output overlap between a pair of files. I need something that can do the following from many bed files and only report regions contained in all of the bed files.

Example input from 4 separate bed files:

chr1    50    100
chr1    60    120
chr1    30    90
chr1    50    90

Desired output:

chr1    60    90

Any tools for this? Maybe I should just cat all the bed files together and merge them?

overlap intersect bedtools genome • 496 views

ADD COMMENT • link •

modified 5 months ago by Alex Reynolds ♦ 26k • written 5 months ago by goodez • 420

5 months ago by

venu ♦ 5.7k

Germany

venu ♦ 5.7k wrote:

Maybe I should just cat all the bed files together and merge them?

You can do that. But there is also multiIntersect. Check here

ADD COMMENT • link written 5 months ago by venu ♦ 5.7k

Thanks! My servers at work had multiIntersectBed already installed, which I guess is also known as bedtools multiinter as genomax suggested. This is just the functionality I was looking for, nothing fancy, just common regions between many bed files.

ADD REPLY • link written 5 months ago by goodez • 420

5 months ago by

genomax ♦ 59k

United States

genomax ♦ 59k wrote:

Intersect multiple BED files (bedops option)

There is also

 bedtools multiinter

Tool:    bedtools multiinter (aka multiIntersectBed)
Version: v2.26.0
Summary: Identifies common intervals among multiple
         BED/GFF/VCF files.

Usage:   bedtools multiinter [OPTIONS] -i FILE1 FILE2 .. FILEn
         Requires that each interval file is sorted by chrom/start.

ADD COMMENT • link modified 5 months ago • written 5 months ago by genomax ♦ 59k

Thanks to you and venu for pointing me in the right direction. +1 for both

ADD REPLY • link written 5 months ago by goodez • 420

5 months ago by

Alex Reynolds ♦ 26k

Seattle, WA USA

Alex Reynolds ♦ 26k wrote:

To intersect intervals from 1 to N files, simply use bedops --intersect:

$ bedops --intersect A.bed B.bed ... N.bed > answer.bed

More details from the documentation: http://bedops.readthedocs.io/en/latest/content/reference/set-operations/bedops.html#intersect-i-intersect

This requires sorted BED files. You can use BEDOPS sort-bed to do this quickly.

I think the other toolkit suggestions now require sorted BED files, as well.

ADD COMMENT • link modified 5 months ago • written 5 months ago by Alex Reynolds ♦ 26k

I was able to find all overlapping regions of my BEDs using bedtools multiinter, although I had to filter the output to only keep regions found in all BEDs.

I'll try out bedops --intersect too since the output might be easier to keep in my pipeline. To sort the BEDs, I found bedtools sort worked great.

ADD REPLY • link written 5 months ago by goodez • 420

bedops --intersect will give you the answer you want without additional steps.

ADD REPLY • link written 5 months ago by Alex Reynolds ♦ 26k

Thanks it was just one pipe to awk using bedtools, but I'm going to switch to bedops. +1 for every answer lol they're all good

ADD REPLY • link written 5 months ago by goodez • 420

Please log in to add an answer.

Similar posts • Search »