Question: Using Bed File With Cuffdiff To Determine Expression Of Specific List Of Genes.
0
gravatar for nepalis
5.3 years ago by
nepalis10
nepalis10 wrote:

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!

rnaseq bed cuffdiff • 2.1k views
ADD COMMENTlink modified 5.2 years ago by Devon Ryan86k • written 5.3 years ago by nepalis10
0
gravatar for Devon Ryan
5.2 years ago by
Devon Ryan86k
Freiburg, Germany
Devon Ryan86k wrote:

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 COMMENTlink written 5.2 years ago by Devon Ryan86k
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: 1210 users visited in the last hour