Using Bed File With Cuffdiff To Determine Expression Of Specific List Of Genes.
1
0
Entering edit mode
8.1 years ago
nepalis ▴ 10

Hi,

I have a list of genes in bed format. I need to determine the expression levels of this genes in the experiment. I have aligned the file using tophat and i have cufflinks results too from older runs. Is there a way to determine the expression level for the genes in the bed format for tophat alignment and cufflinks output ? Is converting the bed file to gtf file and rerunning the rnaseq pipeline again only way ?

thanks!

cuffdiff bed rnaseq • 2.9k views
ADD COMMENT
0
Entering edit mode
8.0 years ago

It's a bit unusual to have genes in bed format, but I suppose that could be used. If you had the genes in GTF or GFF format, the normal reply would be to simply use htseq-count, but that won't work with your bed file. Your best bet is to load the bed file into R and create a GRangesList from it. Then you can open your SAM/BAM files and use countOverlaps to get per-gene counts. I suspect that creating a properly formatted GRangesList from a BED file will prove non-trivial, but it should be doable.

ADD COMMENT

Login before adding your answer.

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