Go Back   SEQanswers > Bioinformatics > Bioinformatics

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

Thread Tools
Old 10-07-2010, 02:28 PM   #1
Junior Member
Location: Uni of Miami, FL

Join Date: Apr 2010
Posts: 1
Default TopHat & Cufflinks failing to assemble full length transcripts


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?

jlhaner is offline   Reply With Quote
Old 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.
Cole Trapnell is offline   Reply With Quote
Old 10-11-2010, 11:10 PM   #3
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.
litc is offline   Reply With Quote
Old 10-13-2010, 11:46 AM   #4
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.
lshen is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 07:59 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO