SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
SSPACE: a new stand-alone scaffolding tool for small and large genomes boetsie Bioinformatics 251 06-01-2016 08:46 AM
SVDetect - No links have been found with the selected type of structural variation tarunmishra.rsg Bioinformatics 10 10-31-2012 08:24 AM
SVDetect : unable to get links.filtered.unique.defined tarunmishra.rsg Bioinformatics 2 10-10-2012 10:02 AM
Is SSPACE good for Abyss assemblies? pmiguel Bioinformatics 25 11-21-2011 07:26 AM
Breast cancer research papers or articles links to bioinformatics... Patrick Bioinformatics 4 09-23-2009 07:30 AM

Reply
 
Thread Tools
Old 05-17-2013, 08:23 AM   #1
swe
Junior Member
 
Location: Sweden

Join Date: Sep 2010
Posts: 2
Default Scaffolding despite multiple links in SSPACE

I have denovo assembled a 1.8 Mbp genome with a PE 100 bp 360 insert library using abyss resulting in 69 contigs. After that I scaffolded these contigs with a MP 100 bp 3300 insert library using SSPACE-Premium resulting in 7 scaffolds.

My problem is that SSPACE makes connections between contigs although the contig link to multiple other contigs. This is the command I used for running SSPACE.

Code:
SSPACE_Premium_v2.2.pl -l libSSPACE.txt -s contigs.fna -k 20 -a 0.9 -x0 -b outdir -T 12 -p1

libraryfile libSSPACE:
lib4000 bwa rawseq_MP_1.txt.gz rawseq_MP_2.txt.gz 3356 0.2 RF

When running this I get a connection between contig 37 and contig 15 although contig 37 has links to 4 different contigs.

r37 has 1118 links with r14 and gap of 351 bases
r37 has 1051 links with r36 and gap of 488 bases
r37 has 595 links with f12 and gap of 334 bases
f37 has 542 links with f15 and gap of -242 bases

I guess that this depends on the building scaffolds part with the two different two methods/ratios for generating reliable scaffolds.

I do not trust these results and how can I set the parameters to get reliable scaffolds with my long insert library?
swe is offline   Reply With Quote
Old 05-28-2013, 07:27 AM   #2
seb.lees
Member
 
Location: France, Poitiers

Join Date: Sep 2012
Posts: 12
Default

Hi swe,

SSPACE choose the 'best' contig connection according to the number of links between contigs (defined as the number of pairs connecting the two contigs). If your r37 have more links to r15 than to the other contigs, r37-r15 will be chosen as the preferred connection.

The -a parameter controls the threshold difference. To my understanding, the lower the more stringent. For instance, a 0.9 value means that if 10 links connect r37 and r15, and 9 links connect r37 and other contigs, r37 and r15 will be scaffolded because you have a ratio of 0.9. So you could reduce your -a parameter to avoid connections you don't trust into.

See the Suplementary Figure S2 from the SSPACE paper: http://bioinformatics.oxfordjournals...4/578.abstract

However in your case, I don't know whether the ratio will be computed as: links r37-r15 vs links r37-all other contigs concurrently, or as: links r37-r15 vs links r37-each other contig independently.

Seb.

Last edited by seb.lees; 05-28-2013 at 07:51 AM.
seb.lees is offline   Reply With Quote
Reply

Tags
long insert, mate pair, scaffolding, sspace

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 08:34 AM.


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