MCMC-ODPR: Primer design optimization using Markov Chain Monte Carlo sampling

Article English OPEN
Kitchen, James L ; Moore, Jonathan D ; Palmer, Sarah A ; Allaby, Robin G (2012)
  • Publisher: BioMed Central
  • Journal: BMC Bioinformatics, volume 13, pages 287-287 (issn: 1471-2105, eissn: 1471-2105)
  • Related identifiers: pmc: PMC3561117, doi: 10.1186/1471-2105-13-287
  • Subject: Molecular Biology | QA | R858-859.7 | Computer applications to medicine. Medical informatics | Software | QH301 | Biochemistry | Computer Science Applications | Biology (General) | QH301-705.5

<p>Abstract</p> <p>Background</p> <p>Next generation sequencing technologies often require numerous primer designs that require good target coverage that can be financially costly. We aimed to develop a system that would implement primer reuse to design degenerate primers that could be designed around SNPs, thus find the fewest necessary primers and the lowest cost whilst maintaining an acceptable coverage and provide a cost effective solution. We have implemented Metropolis-Hastings Markov Chain Monte Carlo for optimizing primer reuse. We call it the Markov Chain Monte Carlo Optimized Degenerate Primer Reuse (MCMC-ODPR) algorithm.</p> <p>Results</p> <p>After repeating the program 1020 times to assess the variance, an average of 17.14% fewer primers were found to be necessary using MCMC-ODPR for an equivalent coverage without implementing primer reuse. The algorithm was able to reuse primers up to five times. We compared MCMC-ODPR with single sequence primer design programs Primer3 and Primer-BLAST and achieved a lower primer cost per amplicon base covered of 0.21 and 0.19 and 0.18 primer nucleotides on three separate gene sequences, respectively. With multiple sequences, MCMC-ODPR achieved a lower cost per base covered of 0.19 than programs BatchPrimer3 and PAMPS, which achieved 0.25 and 0.64 primer nucleotides, respectively.</p> <p>Conclusions</p> <p>MCMC-ODPR is a useful tool for designing primers at various melting temperatures at good target coverage. By combining degeneracy with optimal primer reuse the user may increase coverage of sequences amplified by the designed primers at significantly lower costs. Our analyses showed that overall MCMC-ODPR outperformed the other primer-design programs in our study in terms of cost per covered base.</p>
  • References (13)
    13 references, page 1 of 2

    1. Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol 2000, 132:365-386.

    2. You FM, Huo N, Gu YQ, Luo MC, Ma Y, Hane D, Lazo GR, Dvorak J, Anderson OD: BatchPrimer3: a high throughput web application for PCR and sequencing primer design. BMC Bioinforma 2008, 9:253-253.

    3. Ye J, Coulouris G, Zaretskaya I, Cutcutache I, Rozen S, Madden T: PrimerBLAST: A tool to design target-specific primers for polymerase chain reaction. BMC Bioinforma 2012, 13:134.

    4. Najafabadi HS, Torabi N, Chamankhah M: Designing multiple degenerate primers via consecutive pairwise alignments. BMC Bioinforma 2008, 9:55.

    5. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol 1990, 215(3):403-410.

    6. Linhart C, Shamir R: The degenerate primer design problem. Bioinformatics 2002, 18(Suppl 1):172-181.

    7. Souvenir R, Buhler J, Stormo G, Zhang W: An iterative method for selecting degenerate multiplex PCR primers. Methods Mol Biol 2007, 402:245-268.

    8. Doi K, Imai H: A Greedy Algorithm for Minimizing the Number of Primers in Multiple PCR Experiments. Genome Inform Ser Workshop Genome Inform 1999, 10:73-82.

    9. Rachlin J, Ding C, Cantor C, Kasif S: MuPlex: multi-objective multiplex PCR assay design. Nucleic Acids Res 2005, 33(Web server issue):W544-W547.

    10. Liu YT, Carson DA: A novel approach for determining cancer genomic breakpoints in the presence of normal DNA. PLoS One 2007, 2(4):e380. doi:10.1371/journal.pone.0000380.

  • Metrics
    No metrics available
Share - Bookmark