SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
Velvet and Oases - choice of k value? MikeyG Bioinformatics 5 05-08-2014 09:48 PM
de novo assembly using Trinity versus Velvet-Oases Nol De novo discovery 8 10-26-2013 11:56 AM
RNA-Seq: Oases: Robust de novo RNA-seq assembly across the dynamic range of expressio Newsbot! Literature Watch 0 03-01-2012 07:10 AM
Velvet-Oases Memory Problem dacotahm Bioinformatics 0 10-10-2011 09:08 AM
Velvet/Oases transcript expression level? simon.jarman@aad.gov.au RNA Sequencing 3 08-11-2010 09:47 AM

Reply
 
Thread Tools
Old 11-06-2012, 12:28 PM   #1
AdrianP
Senior Member
 
Location: Ottawa

Join Date: Apr 2011
Posts: 130
Default Denovo RNA-Seq assembly using Velvet/Oases

Hello dear colleagues,

I feel like I should know the answer to what I am about to ask but I could not find a paper with a clear cut answer to what protocol I should apply. I am getting down to it:

I have illumina pair-end reads from RNA-Seq of a fish that is very distant from anything sequenced. I have about 100 mil reads in both directions.

I used Velvet Optimiser to perform denovo assemblies on has lengths 67 77 83 and 93 in order to look at the results and chose the best assembly. However, this is what I encountered (velvet):
Hash 67 has about 250k contigs
Hash 77 has about 125k contigs
Hash 83 has about 78k contigs
Hash 95 has about 15k contigs

So not that big of a suprise, lower hashes more contigs, higher hashes are stricter because of their length, so less contigs.

I polished these using Oasis, and the numbers did not change a lot.
250k 149k 93k 16k
n50s: 867 897 889 827

I tried the -merge function suggested in the Oasis manual, and got 324k contigs and a bit higher n50 about 1000.

Here comes the question, now what? which one is the best assembly? they can't all be as good. Which one should I annotate?

Any help? I am very stuck.
AdrianP is offline   Reply With Quote
Old 11-11-2012, 03:55 PM   #2
AdrianP
Senior Member
 
Location: Ottawa

Join Date: Apr 2011
Posts: 130
Default

bump. if anyone has any input i would love to hear it.
AdrianP is offline   Reply With Quote
Old 11-16-2012, 12:18 AM   #3
CPCantalapiedra
Member
 
Location: Zaragoza (Spain)

Join Date: Sep 2011
Posts: 38
Default

just asking, have you tried other software, like trinity to compare and get more clues for data interpretation?
Maybe you could try too to compare some of the contigs with very highly conserved proteins among vertebrates?
CPCantalapiedra is offline   Reply With Quote
Old 11-27-2012, 06:07 AM   #4
vbiaudet
Member
 
Location: Paris

Join Date: Apr 2011
Posts: 13
Default

You could consider each kmer, I used velvet and oases with each kmer in descending kmer:
I began with kmer=95
--> take the contigs and us it (-long parameter) in the next velvet/oases with kmer=83
--> continue this until the small kmer and you obtain one list of contigs

After you have to use an other assembly of contigs as TGICL or iAssembler

I applied this method to Arabidopsis transcripts and find 90% of contigs that are good compare to gene model...

Last edited by vbiaudet; 11-27-2012 at 07:34 AM.
vbiaudet is offline   Reply With Quote
Old 11-27-2012, 06:09 AM   #5
AdrianP
Senior Member
 
Location: Ottawa

Join Date: Apr 2011
Posts: 130
Default

how far down do I go?

Can anyone else explain to me what confidence is? I understand that it is when there are alternatively spliced transcripts, but what does a very low confidence mean?
AdrianP is offline   Reply With Quote
Old 01-03-2013, 02:04 AM   #6
2seq
Junior Member
 
Location: Israel

Join Date: Oct 2012
Posts: 9
Default

bump...would like to hear more about this. I also like vbiaudet 's idea very much for getting one list of contigs. However, how does this message differ from the oasis merge function?
2seq 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:10 AM.


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