On the Primer Selection Problem in Polymerase Chain Reaction Experiments

Abstract In this paper we address the problem of primer selection in polymerase chain reaction experiments. We prove that the problem of minimizing the number of primers required to amplify a set of DNA sequences is NP -complete. Moreover, we show that it is also intractable to approximate solutions to this problem to within a constant times optimal. We develop a branch-and-bound algorithm that solves the primers minimization problem within reasonable time for typical instances. Next, we present an efficient approximation scheme for this problem, and prove that our heuristic always produces solutions with cost no worse than a logarithmic factor times optimal. Finally, we analyze a weighted variant, where both the number of primers as well as the sum of their costs are optimized simultaneously. We conclude by addressing the empirical performance of our methods on biological data.

[1]  P. Green,et al.  OSP: a computer program for choosing PCR and DNA sequencing primers. , 1991, PCR methods and applications.

[2]  Vasek Chvátal,et al.  A Greedy Heuristic for the Set-Covering Problem , 1979, Math. Oper. Res..

[3]  David S. Johnson,et al.  Approximation algorithms for combinatorial problems , 1973, STOC.

[4]  Carsten Lund,et al.  On the hardness of approximating minimization problems , 1994, JACM.

[5]  William R. Pearson,et al.  A New Approach to Primer Selection in Polymerase Chain Reaction Experiments , 1995, ISMB.

[6]  David S. Johnson,et al.  Computers and Intractability: A Guide to the Theory of NP-Completeness , 1978 .

[7]  C. Dieffenbach,et al.  A computer program for selection of oligonucleotide primers for polymerase chain reactions. , 1990, Nucleic acids research.

[8]  W R Pearson,et al.  Molecular characterization of alpha 1- and alpha 2-adrenoceptors. , 1991, Trends in pharmacological sciences.

[9]  K. Lucas,et al.  An improved microcomputer program for finding gene- or gene family-specific oligonucleotides suitable as primers for polymerase chain reactions or as probes , 1991, Comput. Appl. Biosci..

[10]  William R. Pearson,et al.  Molecular characterization of α1- and α2-adrenoceptors , 1991 .

[11]  J. Sninsky,et al.  Recent advances in the polymerase chain reaction , 1991, Science.

[12]  László Lovász,et al.  On the ratio of optimal integral and fractional covers , 1975, Discret. Math..

[13]  N. Arnheim,et al.  PCR analysis of DNA sequences in single cells: single sperm gene mapping and genetic disease diagnosis. , 1990, Genomics.

[14]  Reuven Bar-Yehuda,et al.  A Linear-Time Approximation Algorithm for the Weighted Vertex Cover Problem , 1981, J. Algorithms.