Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • PacBio, SMART, assembly problem of de novo small genome

    Hey,

    I am very new here and I don't have mach experience with bioinfo and sequencing.

    I need your help because i recently send the DNA of one bacteria to be sequence by PacBio technology and i asked the sequencing company that I want a de novo genome assembly.

    At the beginning they send a report in witch they imply that my DNA was contaminated before sequencing. As i told them that priory to send them the DNA I dis an 16S sequencing with another company and all was good, they replay that is their mistake and they will re-assembly the genome again.

    Now I just got their new report and is a mess. (I am attaching the report).

    I don't even know what to ask from them anymore. Insted of one circular genome I have one non-circular chromosome and 7!!!!! plasmids.

    Any of your advise will be very welcomed.Click image for larger version

Name:	Results 1.jpg
Views:	1
Size:	87.9 KB
ID:	309555

    Click image for larger version

Name:	Results 2 (1).jpg
Views:	1
Size:	83.9 KB
ID:	309556

  • #2
    Nothing about the read metrics raises warning flags - they could certainly have generated longer library molecules though since they did not multiplex the library.

    Is the genome size about what you expected?

    I would try additional assemblies with subsets of the data for examples only the longest reads. 50x coverage would be plenty. I would also try Canu as an alternative assembler.

    Comment


    • #3
      Thx for your suggestion.

      The genome size should be bigger, around 5.5 Mb.

      I'm thinking that some of the "plasmids" are actually part of the chromosome. Or is also possible to have 2 chromosomes. This future can be found in this bacteria genus.

      Comment


      • #4
        Originally posted by iuliachiciudean View Post
        Thx for your suggestion.

        The genome size should be bigger, around 5.5 Mb.

        I'm thinking that some of the "plasmids" are actually part of the chromosome. Or is also possible to have 2 chromosomes. This future can be found in this bacteria genus.
        Judging from the total size of the contigs that would be my suspicion, because it looks like you have the whole 5.5Mb there in some form - just not in a single contig (which you can't guarantee, even with PacBio, all the time).

        Comment


        • #5
          Originally posted by Bukowski View Post
          Judging from the total size of the contigs that would be my suspicion, because it looks like you have the whole 5.5Mb there in some form - just not in a single contig (which you can't guarantee, even with PacBio, all the time).
          So do you think would be possible to solve this only by doing a new assembly with different parameters (which ones?)?

          Or should we do more sequencing?

          Comment


          • #6
            I think the confusion is that they are calling all the contigs except for the largest "plasmid" which would be premature to do so without other evidence, especially if they are called as non-circular. Other than that, it is an OK PacBio assembly. I'd say a majority of the time a bacterial genome assembles as a single contig but a significant minority is fragmented into multiple contigs.

            I'd blast the resulting contigs. Usually it is obvious if it is a plasmid from the blast results. If you really really want a better assembly then it might be worth to re-sequence now that newer kits have been released with longer read lengths. The short fragment size may have been due to DNA quality, though.

            I would try Canu assembly with a higher expected coverage since you have the read depth to do so. The only other thing I can think of is that some libraries have high rates of palindromic reads due to un-repaired ends folding back and acting as a barbell adapter. This can mess up assembly. The newest version of the SMRT assembler chooses a single subread from each well to avoid some of the problems with the palindromes (I think) so I'd check if they assembled with the System 6 update. We filter palindromes before Canu to avoid the problem.
            Providing nextRAD genotyping and PacBio sequencing services. http://snpsaurus.com

            Comment

            Latest Articles

            Collapse

            • 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
            • seqadmin
              The Impact of AI in Genomic Medicine
              by seqadmin



              Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
              02-26-2024, 02:07 PM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, 03-14-2024, 06:13 AM
            0 responses
            34 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-08-2024, 08:03 AM
            0 responses
            72 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-07-2024, 08:13 AM
            0 responses
            82 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-06-2024, 09:51 AM
            0 responses
            68 views
            0 likes
            Last Post seqadmin  
            Working...
            X