SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Oases: De novo transcriptome assembly of very short reads lcollado De novo discovery 58 02-07-2017 08:48 AM
query regarding singletons obtained in the assembly with 454 reads icebreaker Bioinformatics 0 01-03-2012 12:44 AM
Which assembler for de-novo Illumina transcriptome assembly with relatively few reads kmkocot Bioinformatics 1 05-17-2011 03:13 AM
PubMed: Local De Novo Assembly of RAD Paired-End Contigs Using Short Sequencing Reads Newsbot! Literature Watch 0 05-05-2011 11:40 PM
454 and Illumina data; classifying reads/contigs poisson200 Bioinformatics 0 10-29-2010 02:53 AM

Reply
 
Thread Tools
Old 03-01-2011, 06:44 AM   #1
Vickenstein
Junior Member
 
Location: Boston

Join Date: Mar 2011
Posts: 2
Default Assembling De Novo 454 Transcriptome Contigs and Singletons with Illumina Short Reads

Background: I have assembled 4 million None Normalized 454 Titanium de novo transcriptome Reads using Newbler 2.5, and have gotten amazing results, 88000 contigs with N50 of 951 bp, 35000 isotigs with N50 1500 longest being 7900, and 276564 singletons. Now I am about to sequence the same transcriptome using Illumina 75+ bp platfrom.
Problem: I am contemplating the assembly software that I should be using.
I really like newber's isoform prediction, and am not sure if it possible to merge both sets of raws reads together in an good assembly.
I haven't seen any software out that are able to utilize a transcriptome reference for assembling new reads
Vickenstein is offline   Reply With Quote
Old 03-02-2011, 05:58 AM   #2
aparna
Member
 
Location: USA

Join Date: Feb 2009
Posts: 15
Default

I worked on Newbler on transcriptome data,while the concept of isotigs is good- 88k contigs with N50 951 is horrible.
With Illumina data I would advise you to look at Oasis/velvet and scripture.
aparna is offline   Reply With Quote
Old 03-02-2011, 08:43 AM   #3
cram
Member
 
Location: Canada

Join Date: Nov 2008
Posts: 16
Default

Quote:
I worked on Newbler on transcriptome data,while the concept of isotigs is good- 88k contigs with N50 951 is horrible.
Actually, an N50 of 951 sounds very good to me if you're doing a de-novo transcript assembly of some reasonably complex eukaryote. Remember too, Newbler uses a different definition of contig than most other assemblers, and the isotig N50 is probably a better value to use when comparing to other tools.
cram is offline   Reply With Quote
Old 03-03-2011, 05:23 AM   #4
flxlex
Moderator
 
Location: Oslo, Norway

Join Date: Nov 2008
Posts: 415
Default

Try adding the Illumina reads to newbler 2.5! See (my blog): http://contig.wordpress.com/2011/01/...her-platforms/
flxlex is offline   Reply With Quote
Old 03-03-2011, 06:31 AM   #5
Vickenstein
Junior Member
 
Location: Boston

Join Date: Mar 2011
Posts: 2
Default

Thanks for the reply. From the look of it I might have to reassemble the raw reads from 454 and the new reads from Illumina using both Oasis/velvet and Newbler. I will compare the results between these two methods.
Vickenstein is offline   Reply With Quote
Old 03-04-2011, 04:43 AM   #6
BaCh
Member
 
Location: Germany

Join Date: May 2008
Posts: 79
Default

Quote:
Originally Posted by Vickenstein View Post
[...] Now I am about to sequence the same transcriptome using Illumina 75+ bp platfrom. [...]
I am contemplating the assembly software that I should be using.
I use the current development version of MIRA (V3.2.1.8) and just went through a RNASeq denovo 100bp with 22m reads.

Quote:
Originally Posted by Vickenstein View Post
I really like newber's isoform prediction, and am not sure if it possible to merge both sets of raws reads together in an good assembly.
It is. I regularly use MIRA for genome de-novo with 454 and Illumina (ranging from 36 to 100mers). Should also work with mixed transcriptome.

Quote:
Originally Posted by Vickenstein View Post
I haven't seen any software out that are able to utilize a transcriptome reference for assembling new reads
MIRA. No problem if your reference is a transcriptome, but stay away from trying to map RNASeq to a genome, that will fail miserably at intron/exon boundaries.

B.

Disclaimer 1: I'm the author of MIRA, your mileage may vary (but then I'd like to hear about it)
Disclaimer 2: for data sets with more than 40m reads you probably want to wait for a next version.
BaCh is offline   Reply With Quote
Old 03-04-2011, 12:22 PM   #7
sklages
Senior Member
 
Location: Berlin, DE

Join Date: May 2008
Posts: 628
Default

Quote:
Originally Posted by Vickenstein View Post
Thanks for the reply. From the look of it I might have to reassemble the raw reads from 454 and the new reads from Illumina using both Oasis/velvet and Newbler. I will compare the results between these two methods.
Depending on the library itself, I suspect that coverage issues might also prevent a "good assembly" with non-normalized data. I had a 3mio-titanium-reads-human-non-normalzied-cDNA dataset;one fourth of the whole library representing the same gene ... Newbler even failed to map and assemble this set (crashed during consensus calculation). Non-normalized libraries are not the best option to assemble (denovo) with NGS data .. but as I said, it depends on the library.

Sven
sklages is offline   Reply With Quote
Old 03-05-2011, 12:43 AM   #8
BaCh
Member
 
Location: Germany

Join Date: May 2008
Posts: 79
Default

Quote:
Originally Posted by sklages View Post
I had a 3mio-titanium-reads-human-non-normalzied-cDNA dataset;one fourth of the whole library representing the same gene ... Newbler even failed to map and assemble this set (crashed during consensus calculation).
*cough* 750k times the same gene? I would not expect any program to really assemble that de-novo if not specially primed. In mapping too, the coverage might be somewhat on the unexpected side.

B.
BaCh 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:55 AM.


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