Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • 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?

  • #2
    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, 06:51 AM.

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Strategies for Sequencing Challenging Samples
      by seqadmin


      Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
      03-22-2024, 06:39 AM
    • seqadmin
      Techniques and Challenges in Conservation Genomics
      by seqadmin



      The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

      Avian Conservation
      Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
      03-08-2024, 10:41 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Yesterday, 06:37 PM
    0 responses
    11 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, Yesterday, 06:07 PM
    0 responses
    10 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-22-2024, 10:03 AM
    0 responses
    51 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-21-2024, 07:32 AM
    0 responses
    68 views
    0 likes
    Last Post seqadmin  
    Working...
    X