Getting a gene list fom eCLIP bigBED file
1
0
Entering edit mode
4.6 years ago
o.hickman ▴ 10

I'm using eCLIP .bigbed files downloaded from the ENCODE project. I'm try to get a gene list - which genes have tags.

I've been trying to use the "Operate on genomic intervals - Join" function in Galaxy. This I've used before for .bed files. But it doesn't seem to work on .bigbed files.

Converting them to .bed files appears to be non-trivial as well.

Any help getting this to work, or an alternative method for getting a gene list for these data would be most appreciated!!!

Thanks in advance.

eCLIP bed bigbed galaxy encode • 1.6k views
ADD COMMENT
2
Entering edit mode
4.6 years ago

I don't know what the format of an eCLIP file is, but you can convert BigBed to BED with the bigBedToBed tool in the Kent toolkit.

If you're on OS X, you can install them via Homebrew:

$ brew install homebrew/science/kent-tools

For other platforms, you can install binaries (including bigBedToBed) one by one via http://hgdownload.cse.ucsc.edu/admin/exe/

Once you have this tool (or toolkit) installed, you run:

$ bigBedToBed input.bb output.bed

Once you have the BED file, you can do set operations via the BEDOPS kit with bedops, bedmap, etc., e.g.:

$ bedmap --echo --echo-map-id-uniq eCLIP.bed genes.bed > gene_IDs_overlapping_eCLIP_regions.bed
ADD COMMENT
0
Entering edit mode

Hi Thanks Alex, Thanks for your answer. I'm using windows. What platform do you suggest I use to install these tools? Oliver

ADD REPLY
0
Entering edit mode

Run a version of Linux within VirtualBox. https://en.m.wikipedia.org/wiki/VirtualBox

ADD REPLY

Login before adding your answer.

Traffic: 2087 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