Question

Interpret Cufflinks output to detect novel isoforms

2

Entering edit mode

9.5 years ago

Ana ▴ 20

Good evening,

I followed the tuxedo protocol as described here. With the "normal" outputs that cufflinks provide (I'm searching on cuffdiff output files, but tell me if I'm wrong), how could I determine which of isoforms detected are actually novel ones?

I'm pretty new on this field, and I've found in many posts that using Cufflinks is the answer to find novel isoforms, or differential splicing too. Now that I followed up the "instructions" I'm stuck on trying to interpret the results, after searching all over I could. Searching on cufflinks manual didn't help me either. :/

Can someone explain me or provide me some tutorials/papers/etc, that might help me on this topic please?

Thanks in advance.

RNA-Seq cuffdiff cufflinks novel-isoforms • 6.5k views

ADD COMMENT • link updated 23 months ago by Ram 44k • written 9.5 years ago by Ana ▴ 20

0

Entering edit mode

I have done this by running cufflinks with no annotation, then using bedtools to basically remove all the regions that intersect with known genes... Or if you're interested in antisense transcription, look for transcripts on the opposite strand from known genes. Then apply some filters on length and rpkm and come up with a set of putative novel transcripts. Hopefully followed up by some kind of lab based verification, because I feel like this is all pretty hazy and based on thresholds. This may not be the right way to do it, but it seemed to make sense to me at the time.

ADD REPLY • link 9.5 years ago by Madelaine Gogol 5.3k

0

Entering edit mode

7.5 years ago

jila-hessani • 0

hi get cufflink manual from enter link description here

ADD COMMENT • link 7.5 years ago by jila-hessani • 0

Ram · Accepted Answer · 2015-06-07

Here are two tutorials on RNA-seq, both that involve discussion of Cufflinks that you might find generally useful:

With regard to your specific question. You might find the Cuffcompare and Cuffmerge commands and their documentation useful.

Cuffcompare has this description:

After assembling a transcriptome from one or more samples, you'll probably want to compare your assembly to known transcripts. Even if there is no "reference" transcriptome for the organism you're studying, you may want to compare the transcriptomes assembled from different RNA-Seq libraries. Cuffcompare helps you perform these comparisons and assess the quality of your assembly.

You can run cuffcompare directly but in recent versions of cufflinks, the cuffmerge command will do this for you. You should have '.tracking' output files. Refer to the 'class code' column of this file. The definition of "transfrag class codes" is provided at the bottom of the Cuffcompare manual page. For example, class "j" is described as "Potentially novel isoform (fragment): at least one splice junction is shared with a reference transcript"