Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Demultiplexing a MiSeq library based on Amplicon Sequence divergence?

    Hi All,

    My lab works on a lot of diverse taxa. We want to sequence the mtDNA genomes of each using our MiSeq. We generated long-range PCR amplicons for each species, then pooled 10 taxonomically divergent species into each library (total Miseq run was 10 libraries containing 10 species in each library). Each library is labelled with a Nextera ID. We planned on demultiplexing each library based on the non-similarity in sequences. The run was successful, and now it's time to demultiplex the libraries. We had a guy that was going to do this for us, but he's super busy so we are going to try it ourselves...but we don't know how we should go about this. Any suggestions regarding programs or how to go about it? Since each library contains only amplicon sequence can we simple ask MIRA to build contigs? Any help is appreciated.

    Cheers,

  • #2
    All you want is mtDNA right?

    Just assemble it w/o demultiplexing, you should get mtDNA contigs for each of your taxa.

    How long is your fragment size? With MiSeq + your situation I recommend MIRA.

    Comment


    • #3
      Thanks Adrian. The mtDNA genomes are between 15-18 kb, however we had multiple PCR's of a smaller size (~5-10 kb) to generate a full genome, then we pooled it together. Basically, my pooled library contains mtDNA amplicons of 10 distantly related species (i.e., a bat, a fish, a bird, etc). The 2 closest species are both fish, but still very distantly related. So you think MIRA can just assemble them denovo by asking for contigs to be built? No need for reference genomes (some of which we have, some not)? Thanks for answering!

      Comment


      • #4
        As long as your PCR product are overlapping, you should be fine. Indeed you can just denovo assemble it by feeding it the entire library. Assemble each library individually.

        As a side thing, you might want to check the quality of your data. You said you have some reference genomes. Use bwa to map a library containing the mtDNA of the reference genome to the reference genome sequence, then visualize read coverage with tablet, it will tell you if the entire mtDNA is covered. It's like quality control. If some regions are poorly covered or not at all, then shucks, you might not get your circular contigs.

        Comment


        • #5
          Thanks Adrian, I'll give that a try. Have you ever tried this before? Or heard of it begin done? Would love to hear some 'proof of concept' testimonial!
          Cheers,

          Comment


          • #6
            I have never heard of your initial idea of demultiplexing based on similarity. That's essentially what an assembler does, separates reads based on which reads form contigs...

            There is no proof of concept here, it's plain and simple, you have metagenomics, multiple genomes in one library. Assemblers handle that well. It really depends on the depth that each mtDNA got sequenced at.

            Comment


            • #7
              The assemblers should be OK when the species are not closely related, if the species are too closely related you may get chimeric contigs.

              You may have problems due to uneven coverage, even if you try to pool equimolar amounts of each amplicon, in practice it never comes out quite equimolar, and some amplicons will be present at higher concentrations than others. This causes problems for some assemblers which assume even coverage.

              Do adapter trimming/quality trimming of the reads before assembly.

              As mentioned above, MIRA is a good assembler to try.

              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
              10 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, Yesterday, 06:07 PM
              0 responses
              9 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
              67 views
              0 likes
              Last Post seqadmin  
              Working...
              X