Question: Annotate A Gtf File (Cufflinks)
0
gravatar for Nicolas Rosewick
6.5 years ago by
Belgium, Brussels
Nicolas Rosewick8.1k wrote:

Hi,

I'm working on RNA-Seq data. Here's my basic workflow :

STAR -> cufflinks -> cuffmerge

So after that, I've a merged transcriptome gtf file. Now I want to annotate it to furthyer use it in differential expression analysis. So here's my idea. First extract the transcritp sequence using the gtf coordinates. Then, blast each sequence against ensembl transcript to annotate them.

What do you think ? Is there maybe another tool ?

thanks,

N.

gtf cufflinks • 3.0k views
ADD COMMENTlink modified 6.5 years ago by biorepine1.4k • written 6.5 years ago by Nicolas Rosewick8.1k

Hi, I am in similar situation. I have draft genome but no annotations. I need to do differential splicing analysis between conditions. How should I go about this ? Should I do tophat2 -> cufflinks -> cuffmerge -> Annotate GTF ?

If yes, then what would be the best way to annotate the GTF generated by cufflinks ?

Is there any alternative way for this ?

ADD REPLYlink written 4.5 years ago by geek_y9.9k
0
gravatar for biorepine
6.5 years ago by
biorepine1.4k
Spain
biorepine1.4k wrote:

Convert cufflinks GTF into BED12 format using the following script and overlap the genomic-c0ordinates with known reference genome (RefSeq/UCSC/ENSEMBLE) by using BED tools. It gives you proper annotation. Here no need to to do BLAST and you can save lot of time.

For new transcripts that are not in the above databases, you need to add them. In your case, you can obtain all the gencode and human body map discovered lncRNAs (human-hg19) and add them to the reference genome BED file. That should work.

ADD COMMENTlink modified 6.5 years ago • written 6.5 years ago by biorepine1.4k

The problem is that the annotation of the species is not very good and a lot of things are missing (linke lncRNAs)

ADD REPLYlink written 6.5 years ago by Nicolas Rosewick8.1k

i updated my answer. Hope it helps.

ADD REPLYlink written 6.5 years ago by biorepine1.4k
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: 896 users visited in the last hour