Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • 20kb Plasmid Assembly

    Hello all,

    I've very new to this type of analysis and was hoping I could get some help. I am trying to assemble the sequences of a 20kb plasmid that I had sequenced. They are short single reads that when all added up have over 4000x coverage. What I have done so far is I have trimmed off the first 4 and last 5 bases of the reads, since they were the lowest values when I quality checked them. I am trying to reduce the number of overall reads by quality filtering then down to get somewhere around 50x coverage.

    My plan is to use velvet to assemble, but I am a bit confused about what to do with the output. I have tested it a couple of times and can get the output files, but I have no idea what the next step should be. How do I decide what contigs are good and what are bad? I know that my plasmid is about 20kb, so should I just dismiss anything larger? The contigs file has a lot of different sequences, and I am not sure how to narrow it down from there. Any help would be greatly appreciated!

  • #2
    You might want to try using MUSKET or one of the other k-mer based tools out there to correct errors. This may be more effective than your trimming, and you might want to try assembling both trimmed & untrimmed data

    How long are your reads? Longer is better -- which is why paired end is far better than single end.

    Do you know anything about your plasmid backbone? What host was the plasmid prepared in? Screening out contigs which match the host (e.g. E.coli DH10B) would be a valuable next step. Screening out contigs corresponding to the center of the backbone may be useful -- as well as identifying the vector-insert junctions.

    Comment


    • #3
      I'll look into MUSKET for sure. What I'm starting to wonder is if I should not trim to reduce my coverage, but take a percentage of the reads. The problem is I have not a clue how to do that.

      My reads are 100bp. I trimmed them to 90bp.

      I can't find much on the plasmid backbone, that is the point of the sequencing. I'm trying to find the backbone of a plasmid I am developing. I do know about 3kb of the sequence, since I have inserted it myself. I'm estimating the plasmid to be about 20kb, as that is what makes sense from some restriction digests. I can take portions of the 3kb I inserted and find them in different contigs that velvet spits out, so at least I know that my plasmid is in there.

      Comment


      • #4
        If you have an additional few hundred dollars to commit to the project, why not just run PacBio sequencing? Since the reads are O(length of the plasmid), it becomes MSA rather than assembly...


        Originally posted by gnatnog View Post
        Hello all,

        I've very new to this type of analysis and was hoping I could get some help. I am trying to assemble the sequences of a 20kb plasmid that I had sequenced. They are short single reads that when all added up have over 4000x coverage. What I have done so far is I have trimmed off the first 4 and last 5 bases of the reads, since they were the lowest values when I quality checked them. I am trying to reduce the number of overall reads by quality filtering then down to get somewhere around 50x coverage.

        My plan is to use velvet to assemble, but I am a bit confused about what to do with the output. I have tested it a couple of times and can get the output files, but I have no idea what the next step should be. How do I decide what contigs are good and what are bad? I know that my plasmid is about 20kb, so should I just dismiss anything larger? The contigs file has a lot of different sequences, and I am not sure how to narrow it down from there. Any help would be greatly appreciated!

        Comment


        • #5
          Originally posted by mchaisso View Post
          If you have an additional few hundred dollars to commit to the project, why not just run PacBio sequencing? Since the reads are O(length of the plasmid), it becomes MSA rather than assembly...

          That's a bit of a stretch; you won't have many reads quite that long, and they probably won't survive error correction.

          On the other hand, even one SMRT cell on such a small genome -- barring rampant host contamination -- should give a number of high quality long reads to easily assemble the genome

          Comment


          • #6
            Originally posted by gnatnog View Post
            I can't find much on the plasmid backbone, that is the point of the sequencing. I'm trying to find the backbone of a plasmid I am developing. I do know about 3kb of the sequence, since I have inserted it myself. I'm estimating the plasmid to be about 20kb, as that is what makes sense from some restriction digests. I can take portions of the 3kb I inserted and find them in different contigs that velvet spits out, so at least I know that my plasmid is in there.
            Try mapping the reads back to your backbone with Bowtie2. That would give you an estimate of what coverage you actually achieved for your plasmid. Also worth taking your longest contigs and BLASTing them against all known sequences; sometimes that can be an eye opener.

            I'm not a big fan of trimming, but the truth is I haven't played with it much. One challenge is that many assemblers have logic to deal with these issues, but not always clear how much -- so there is a complex interaction between various pre-processing tools & assemblers, and only empirically can you really work out the best strategy.

            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
            8 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, Yesterday, 06:07 PM
            0 responses
            8 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-22-2024, 10:03 AM
            0 responses
            49 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-21-2024, 07:32 AM
            0 responses
            66 views
            0 likes
            Last Post seqadmin  
            Working...
            X