Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • AP38
    Member
    • May 2016
    • 14

    Estimating the number of ddRAD fragments

    Hi all,

    I am wondering if there is an easy way to estimate the number of fragments resulting from a double digest reaction (in my case ddRAD). For a single digest, one can do the following:

    ((proportion of cut sites)^(enzyme cut frequency)) x (genome size) = number of RAD fragments

    For instance, with ecoRI (6 cutter) and a proportion of 0.25 across a genome of 2.5Gbs:

    (0.25)^6 x 2,500,000,000 = 610,351.6 RAD fragments


    Is there a similar way for calculate this for a double digest reaction? One could use simRAD in R or even use grep or cat commands on the reference fasta file (cat ref.fasta | tr -d "\n" | grep -o -E “enz1sequence.{200,400}enzyme2sequence” | wc -l).

    However, I would like to know if there is an easier solution to this such as the one described above for single digest.

    Thanks!
  • nucacidhunter
    Jafar Jabbari
    • Jan 2013
    • 1250

    #2
    Calculation method that you have described is incorrect because it assumes that nucleotide sequences in genomes are random and ignores the repetitive content. Two correct methods:

    1- Bioinformatically if there is a reference as you have mentioned
    2- Empirical (most accurate)

    Comment

    • AP38
      Member
      • May 2016
      • 14

      #3
      Thanks nucacidhunter,

      You are right about the calculation method. It is not the best because of the point you made but should still give a rough estimate. The empirical method is definitely the most accurate way but I would prefer to have a precise in silico estimate prior to jumping into lab work.

      Would you have a good bioinformatic method to propose?

      Thanks!

      Comment

      • nucacidhunter
        Jafar Jabbari
        • Jan 2013
        • 1250

        #4
        I know that there are some open source software that gives the number and size of the in silico digest fragments. Hopefully someone in the forum can give the names. I found the following thesis interesting in this regard:

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Pathogen Surveillance with Advanced Genomic Tools
          by seqadmin




          The COVID-19 pandemic highlighted the need for proactive pathogen surveillance systems. As ongoing threats like avian influenza and newly emerging infections continue to pose risks, researchers are working to improve how quickly and accurately pathogens can be identified and tracked. In a recent SEQanswers webinar, two experts discussed how next-generation sequencing (NGS) and machine learning are shaping efforts to monitor viral variation and trace the origins of infectious...
          03-24-2025, 11:48 AM
        • seqadmin
          New Genomics Tools and Methods Shared at AGBT 2025
          by seqadmin


          This year’s Advances in Genome Biology and Technology (AGBT) General Meeting commemorated the 25th anniversary of the event at its original venue on Marco Island, Florida. While this year’s event didn’t include high-profile musical performances, the industry announcements and cutting-edge research still drew the attention of leading scientists.

          The Headliner
          The biggest announcement was Roche stepping back into the sequencing platform market. In the years since...
          03-03-2025, 01:39 PM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 03-20-2025, 05:03 AM
        0 responses
        49 views
        0 reactions
        Last Post seqadmin  
        Started by seqadmin, 03-19-2025, 07:27 AM
        0 responses
        57 views
        0 reactions
        Last Post seqadmin  
        Started by seqadmin, 03-18-2025, 12:50 PM
        0 responses
        50 views
        0 reactions
        Last Post seqadmin  
        Started by seqadmin, 03-03-2025, 01:15 PM
        0 responses
        200 views
        0 reactions
        Last Post seqadmin  
        Working...