SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Own a website which hosts plasmids? Rutia Vendor Forum 0 01-26-2016 06:19 AM
SPAdes contigs vpi Bioinformatics 3 12-02-2015 04:49 PM
Identifying plasmids olsb Bioinformatics 1 10-24-2015 10:17 AM
SPAdes trusted contigs aimc Bioinformatics 0 01-12-2015 10:22 AM
Plasmids part of core genome Mona Bioinformatics 2 06-16-2014 06:32 AM

Reply
 
Thread Tools
Old 07-26-2016, 08:48 AM   #1
ronaldrcutler
Member
 
Location: Virginia

Join Date: May 2016
Posts: 80
Default SPAdes plasmids output

Hello all. I have done a plasmidSPAdes run on my bacterial reads to look for any plasmids. In the final output directory, I have the contigs.fasta file which contains the plasmids with the suffix componenet_X denoting each respective plasmid.

When looking at the different kmer value folders the K21-K99 folders seem to contain assemblies of the bacterial contigs in the final_contigs.fasta file, as a normal SPAdes run would. However, the K127 folder final_contigs.fasta is instead similar to the contigs.fasta file mentioned above - it seems to only contain plasmids and not the bacterial contigs.

After comparing the K127 final_contigs.fasta and contigs.fasta files with diff, they have differences. My question is why are these two files different? Is the contigs.fasta file different because SPAdes ran the mismatch corrector? And also why are plasmids only output in the K127 folder final_contigs.fasta and not the other kmer value folders?

Thanks

Last edited by ronaldrcutler; 07-26-2016 at 08:51 PM.
ronaldrcutler is offline   Reply With Quote
Old 07-26-2016, 01:19 PM   #2
ronaldrcutler
Member
 
Location: Virginia

Join Date: May 2016
Posts: 80
Default

Update:

In the K127 folder there is also before_chromosome_removal.fasta file. When I try to see if contig in the final_contigs.fasta is in the aformentioned file, I do not find them. So are the sequences that are supposedly plasmid sequences from the final_contigs.fasta in the K127 folder and contigs.fasta in the main SPAdes output directory in the final_contigs.fasta files in the other folders? What exactly is the difference between them all?
ronaldrcutler is offline   Reply With Quote
Old 07-26-2016, 05:14 PM   #3
HESmith
Senior Member
 
Location: Bethesda MD

Join Date: Oct 2009
Posts: 509
Default

Quote:
And also why are plasmids only output in the K127 folder final_contigs.fasta and not the other kmer value folders?
Increasing the kmer length decreases the depth of coverage for each kmer. At some point, the depth becomes insufficient to assemble the de Bruijn graph accurately. Since plasmid sequences are present at much higher read depth, assembly of those reads can tolerate a higher kmer than the host cell sequences. It appears that setting k=127 is a threshold for assembling plasmid but not host cell reads.

Also, if you're only interested in the plasmid sequences, you may want to use plasmidSPAdes (described here).

Last edited by HESmith; 07-28-2016 at 09:30 AM.
HESmith is offline   Reply With Quote
Old 07-26-2016, 08:56 PM   #4
ronaldrcutler
Member
 
Location: Virginia

Join Date: May 2016
Posts: 80
Default

I was using plasmidSPAdes. If I understand correctly, the plasmid contigs are able to assemble at higher kmer values because they are present at a higher read depth. Why are they more present at higher read depths than the host cell? What about coverage?
ronaldrcutler is offline   Reply With Quote
Old 07-27-2016, 03:27 AM   #5
HESmith
Senior Member
 
Location: Bethesda MD

Join Date: Oct 2009
Posts: 509
Default

Plasmids are present in multiple copies per cell, whereas the host genome is (typically) a single copy. More copies = higher read depth.
HESmith is offline   Reply With Quote
Old 07-27-2016, 12:18 PM   #6
ronaldrcutler
Member
 
Location: Virginia

Join Date: May 2016
Posts: 80
Default

Okay, thanks for the clarification. So since there are no plasmids specified in the other final_contigs.fasta files from using kmer values less than 127, I can confidently say there are no plasmids in these and were only found when using a kmer value of 127?
ronaldrcutler is offline   Reply With Quote
Reply

Tags
bacteria genome, de novo assembly, genome assembly, spades

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 12:35 PM.


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