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
                      Latest Developments in Precision Medicine
                      by seqadmin



                      Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

                      Somatic Genomics
                      “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
                      05-24-2024, 01:16 PM
                    • seqadmin
                      Recent Advances in Sequencing Analysis Tools
                      by seqadmin


                      The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
                      05-06-2024, 07:48 AM

                    ad_right_rmr

                    Collapse

                    News

                    Collapse

                    Topics Statistics Last Post
                    Started by seqadmin, 05-24-2024, 07:15 AM
                    0 responses
                    195 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 05-23-2024, 10:28 AM
                    0 responses
                    217 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 05-23-2024, 07:35 AM
                    0 responses
                    218 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 05-22-2024, 02:06 PM
                    0 responses
                    12 views
                    0 likes
                    Last Post seqadmin  
                    Working...
                    X