Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • ddRADSeq Libraries

    We seem to be having some difficulties creating reproducible ddRADSeq libraries.

    We are noticing that after running the first of the processing steps (process radtags) - which removes the barcodes, filters for the rad sites, and checks for sequence quality - we have some evidence of perhaps incomplete restriction digestion.
    For instance - we have:
    rad cut site (i.e, CATG from NlaIII) - followed by target sequence - followed by the second cut site (i.e., MluCI - AATT) - followed by more sequence.

    Shouldn't the sequence end at the second cut site? Has anyone tried a sequential digest? We have been running the test libraries at 150bp SE on a MiSeq.

    Any help is greatly appreciated.

  • #2
    A single digest is far better than sequential digest if both enzymes are active in the same temperature and buffer (yours are).

    In ideal reaction condition incomplete digestion should be minimal. I wonder if the sequences beyond MluCI site are untrimmed adapter sequences.

    Comment


    • #3
      Thanks for the reply nucacidhunter. After looking over the data again there is a minimal amount of adapter attached after the cut site. The majority of reads are partial digestion by the MluCI restriction enzyme. I'm not sure why this is. We are starting with 200 ng of gDNA and adding 10 units total (1 ul) of each enzyme to a 50 ul reaction. It is then incubating in a thermalcycler for 3 hours at 37C.

      I'm wondering if we should try adding more of the MluCI RE compared to the NlaIII RE or extending the digest to 6 hours.

      Comment


      • #4
        Could you be getting chimeras of adapter - NlaIII - fragment 1 - MluCI - fragment 2 - MluCI - adapter?

        You would see different sequences after the first MluCI for a given locus if so.

        These are both 4-cutter enzymes... are you trying to sample an enormous number of loci?
        Providing nextRAD genotyping and PacBio sequencing services. http://snpsaurus.com

        Comment


        • #5
          You can add more RE provided that the final combined volume does not exceed 10%. It is better to do a test digest on MluCI to rule out bad batch of RE.

          The issue mentioned by SNPsaurus could be result of fragments ligating to each other rather than to adapters which could be overcome by increasing adapter quantity in reaction. It also can result from overcycling.
          Last edited by nucacidhunter; 11-08-2017, 01:47 PM. Reason: Omitting unnecessary line.

          Comment


          • #6
            I was thinking about this some more, and I think the use of two 4-cutters may be creating conditions that enrich for this sort of artifact.

            What is your size selection range? Most people do ddRAD with 300 bp fragments for PE 2x150 reads. If that is the case, you will have very few fragments available at 300 bp, since the genome is being chopped up into 125 bp fragments by the use of two 4-cutters. Even if chimeras happen at a low rate, these chimeras will be now long enough to be in the size selection range, and could be present at that size range as an appreciable fraction compared to the unlikely desired NlaIII-MluCI fragment of that size.
            Providing nextRAD genotyping and PacBio sequencing services. http://snpsaurus.com

            Comment


            • #7
              Wow. Thanks so much for the insight. Right now we are supposedly adding 5 fold adapter concentrations but we could up that. The calculations are based off of the enzyme trial step prior to the protocol, which you would think also has this incomplete cutting.

              Would using a 6-cutter and a 4-cutter enzyme in the double digest be more advantageous? Perhaps SphI (6 cutter) and NlaII (4 cutter).

              We were also thinking about selecting a larger range of basepairs in size selection step. Last round we did 300-400bp on our Pippen Prep. This time we were planning on doing 400-600bp, since this will be SE sequencing.

              We will be receiving a positive control from NEB, to test our enzyme function. Our current project is using deep-sea fish genomic DNA. There is a slight chance that the problem could be with this kind of genomic DNA.

              Comment


              • #8
                I'm curious what your goals are for the project. There might be several million of each site in the genome (if it is 1 Gb). Even if you took the largest 5% of fragments, that could be several hundred thousand loci being amplified and sequenced. Getting good read depth on that number of loci will require lots of sequencing per individual (10M reads or more). Plus the artifacts from SNPs creating new 4-cutter sites in your fragments would mean a huge rate of missing data in any kind of genetically diverse population.

                If you are just looking at population structure, then I would pair up a 6-cutter or 8-cutter with a 4-cutter. You could estimate the number of 6 or 8-cutter sites in the genome, then figure out from your size selection what portion would be kept (so keeping 400-600 bp fragments will mean keeping fragments bigger than the expected ~250 bp (the usual rate of a 4-cutter cutting)... only 10% of fragments will be 500 bp, for example.

                How many samples do you have? What's your budget for sequencing? I'd take the # loci x read depth desired and multiply by 3 or more to get how much sequencing will be needed. Loci amplify at different rates, and samples vary in how many reads they get so you need to pack in some buffer to get good depth at most loci in most samples.
                Providing nextRAD genotyping and PacBio sequencing services. http://snpsaurus.com

                Comment


                • #9
                  Using a 4 and 6 cutter combination and 100 bp window cut should give over 10k tags. I would suggest to use 5 ng of digested-ligated DNA in 10-12 PCR cycle and check the size distribution. The profile should give hints on size selection window and also can help to avoid repeat regions tags.

                  Comment


                  • #10
                    Thank you nucacidhunter and SNPsaurus! We are going to use all of this information to remake the libraries and run another test run on the MiSeq. I will follow-up with results. Thanks again you have been more thank helpful!

                    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
                    67 views
                    0 likes
                    Last Post seqadmin  
                    Working...
                    X