Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • 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!

  • #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


    • #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


      • #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:
        http://www.masterbioinformatica.com/...Z_FERNANDO.pdf

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Exploring the Dynamics of the Tumor Microenvironment
          by seqadmin




          The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
          07-08-2024, 03:19 PM
        • seqadmin
          Exploring Human Diversity Through Large-Scale Omics
          by seqadmin


          In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
          06-25-2024, 06:43 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, Yesterday, 11:09 AM
        0 responses
        16 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 07-19-2024, 07:20 AM
        0 responses
        148 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 07-16-2024, 05:49 AM
        0 responses
        124 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 07-15-2024, 06:53 AM
        0 responses
        111 views
        0 likes
        Last Post seqadmin  
        Working...
        X