SEQanswers

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

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

Join Date: Apr 2010
Posts: 1
Default 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
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
Default

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
litc
Member
 
Location: China

Join Date: Oct 2010
Posts: 24
Default

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
lshen
Member
 
Location: Toronto

Join Date: Jan 2008
Posts: 30
Default

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
Reply

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 06:09 PM.


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