TopHat & Cufflinks failing to assemble full length transcripts

jlhaner · 10-07-2010, 02:28 PM

Hi,

First post on SeqAnswers. The discussions here are very useful.

We are using Tophat (v1.0.13) and Cufflinks (0.8.3) without a reference GTF and then use Cuffcompare to identify the assembled transcripts.

We are finding many transcripts reported as novel isoforms that we suspect are actually just the main transcript being divided into 2 fragments, or are leaving out a few exons at the beginning which are clearly covered by reads.

For example, if the gene has 34 exons, novel isoform j1 is identified as the first 23 exons, and j2 is made of the last 11 exons. Another example, a novel transcript is reported which begins from exon2, however there are an equal # of reads covering the first exon.

Examination of the .wig file shows coverage of the complete transcript but for some reason the full length transcript is not being assembled. We've changed the # of bp on either side of the splice junction, with no avail. We also run the butterfly search, and no change with that option either.

Does anyone have suggestions for us?

Thanks,
Jessica

Cole Trapnell · 10-11-2010, 06:31 AM

Are you certain there are spliced reads connecting those exons as well? You want to visualize the read alignments in IGV to ensure that you have reads spanning all the junctions.

litc · 10-11-2010, 11:10 PM

I'm also interested in mRNA-seq reads assembly. It seems that there is not a good soft do that work because of the alternative splice.

lshen · 10-13-2010, 11:46 AM

Try new version instead of this Cufflinks (0.8.3). I got greatly improved transcript assembling results with my single read data. Lot of them "full" length when compared with existing annotations.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
High performance computing & Tophat & Cufflinks	paolo.kunder	Bioinformatics	1	01-24-2012 09:02 AM
cufflinks not able to assemble transcripts in high read coverage area	apratap	Bioinformatics	2	11-03-2011 08:21 PM
some questions about running tophat & cufflinks	songyj	RNA Sequencing	0	10-18-2011 06:07 PM
TopHat 1.1 failing on colorspace SE reads	krobison	Bioinformatics	22	09-22-2011 12:25 PM
Reads map in Tophat, but Cufflinks doesn't find transcripts	janec	Bioinformatics	0	07-07-2011 02:42 AM

10-07-2010, 02:28 PM	#1
jlhaner Junior Member Location: Uni of Miami, FL Join Date: Apr 2010 Posts: 1	TopHat & Cufflinks failing to assemble full length transcripts Hi, First post on SeqAnswers. The discussions here are very useful. We are using Tophat (v1.0.13) and Cufflinks (0.8.3) without a reference GTF and then use Cuffcompare to identify the assembled transcripts. We are finding many transcripts reported as novel isoforms that we suspect are actually just the main transcript being divided into 2 fragments, or are leaving out a few exons at the beginning which are clearly covered by reads. For example, if the gene has 34 exons, novel isoform j1 is identified as the first 23 exons, and j2 is made of the last 11 exons. Another example, a novel transcript is reported which begins from exon2, however there are an equal # of reads covering the first exon. Examination of the .wig file shows coverage of the complete transcript but for some reason the full length transcript is not being assembled. We've changed the # of bp on either side of the splice junction, with no avail. We also run the butterfly search, and no change with that option either. Does anyone have suggestions for us? Thanks, Jessica

10-11-2010, 06:31 AM	#2
Cole Trapnell Senior Member Location: Boston, MA Join Date: Nov 2008 Posts: 212	Are you certain there are spliced reads connecting those exons as well? You want to visualize the read alignments in IGV to ensure that you have reads spanning all the junctions.

10-11-2010, 11:10 PM	#3
litc Member Location: China Join Date: Oct 2010 Posts: 24	I'm also interested in mRNA-seq reads assembly. It seems that there is not a good soft do that work because of the alternative splice.

10-13-2010, 11:46 AM	#4
lshen Member Location: Toronto Join Date: Jan 2008 Posts: 30	Try new version instead of this Cufflinks (0.8.3). I got greatly improved transcript assembling results with my single read data. Lot of them "full" length when compared with existing annotations.