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


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 ?


rnaseq bed cuffdiff • 2.2k views
ADD COMMENTlink modified 5.5 years ago by Devon Ryan88k • written 5.6 years ago by nepalis10
gravatar for Devon Ryan
5.5 years ago by
Devon Ryan88k
Freiburg, Germany
Devon Ryan88k 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.5 years ago by Devon Ryan88k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 781 users visited in the last hour