Optimization Methods for Reliable Genomic Based Pathogen Detection

Optimization Methods for Reliable Genomic. Based Pathogen Detection Systems K. M. Konwar, I. I. Mandoiu, A. C. Russell, and A. A. Shvartsman Computer Science & Engineering Department University of Connecticut, Storrs, CT 06269

Abstract Recent advances in genomic technologies have opened the way for the development of Genomic-based Pathogen Detection Systems (GPDSs) that can provide early warning in case of rapidly proliferating outbreaks of new natural pathogens such as the SARS corona-virus or bio-terrorist attacks. A critical step of all GPDS architectures proposed to date is DNA amplification by Multiplexed Polymerase Chain Reaction (MP-PCR). In this poster we present ongoing theoretical and practical research on the minimum primer set selection problem for MPPCR. We give algorithms with improved approximation guarantees for this problem and report results of empirical experiments on both synthetic and public genomic database test cases showing that our algorithms are highly scalable and produce better results compared to previous heuristics.

GPDS Components and Requirements m Key GPDS components: q q m Selection of distinguishing DNA oligonucleotides based on available genomic sequences for the pathogens Selective amplification of collected genetic material Hybridization-based detection of present distinguishers Pathogen identification by comparison with stored signatures/barcodes of known pathogens GPDS design requirements q q High specificity and sensitivity of detection Discrimination between pathogens and non-pathogenic organisms Ability to work with trace amounts of genetic material, and to detect multiple pathogens at the same time Fully automated operation (should require minimal human intervention

Pathogen Detection System Architecture PCR Machine Set of (degenerate) primers Multiplex PCR Fluorescent nucleotides Set of probes combining barcode distinguishers with tags Multiplex PCR Mixture of (degenerate) primers Multiplex PCR Sample possibly containing minute traces of several pathogens Mixture of (degenerate) primers … … Amplified DNA sequences from sample Universal DNA Tag Array Distinguishing signature of pathogens present in sample

The Polymerase Chain Reaction Target Sequence Primers 5’ 3’ 3’ Polymerase 5’ 3’ r 1 Primer 2 Prime 3’ 5’ 5’ Repeat 20 -30 cycles Invented by Mullis in 1980 s, PCR uses short oligonucleotide primers and the DNA polymerase in a cyclic reaction to produce millions of copies of a target sequence of DNA

Primer Pair Selection Problem 5' 3' Reverse primer L Forward primer 3' 5' amplification locus • Given: • Genomic sequence around amplification locus • Primer length k • Amplification upperbound L • Find: Forward and reverse primers of length k that hybridize within a distance of L of each other and optimize amplification efficiency (melting temperatures, secondary structure, cross hybridization, etc. )

Multiplex PCR m Multiplex PCR (MP-PCR) q q q m Multiple DNA fragments amplified simultaneously Boundaries of each amplification fragment still defined by two oligonucleotide primers A primer may participate in the amplification of multiple targets Primer set selection q q Typically done by time-consuming trial and error An important objective is to minimize the total number of primers Ø Reduced assay cost Ø Higher effective concentration of primers higher amplification efficiency Ø Reduced unintended amplification Original DNA Template PCR Products

Other Applications of Multiplex PCR m Spotted microarray synthesis [Fernandes&Skiena’ 02] q q m Need unique pair for each one of the n amplification product, but primers can be used multiple times Potential to reduce #primers from O(n) to O(n 1/2) SNP Genotyping q q Thousands of SNPs that must genotyped using hybridization based methods (e. g. , single-base extension) Selective PCR amplification needed to improve accuracy of detection steps (whole-genome amplification less appropriate) No need for unique amplification! Primer minimization is critical Ø Reduced cost Ø Fewer multiplex PCR reactions

Primer Set Selection Problem • Given: • Genomic sequences around each amplification locus • Primer length k • Amplification upper bound L • Find: • Minimum size set of primers S of length k such that, for each amplification locus, there are two primers in S hybridizing to the forward and reverse sequences within a distance of L of each other • For some applications: S should contain a unique pair of primers amplifying each locus

Previous Work on Primer Set Selection m All previous works, e. g. , [Pearson et al. 96][Linhart & Shamir’ 02] [Souvenir et al. ’ 03], use problem formulations that decouple selection of forward and reverse primers, and hence cannot directly enforce constraints on amplification product length q To enforce bound of L on amplification length, select only primers that hybridize within L/2 bases of desired target Ø Ignores half of the feasible primer pairs! q m In worst case, this method can increase the number of primers by a factor of O(n) compared to the optimum Greedy set cover algorithm gives O(ln n) approximation factor for the “decoupled” formulation q Cannot find better approximation unless P=NP

Previous Work (contd. ) m [Fernandes&Skiena’ 02] model primer selection as a minimum multicolored subgraph problem: q q There is an edge colored by color i between primers u and v if they hybridize to i-th forward and reverse sequences within a distance of L q m Vertices of the graph correspond to candidate primers Goal is to find minimum size set of vertices inducing edges of all colors No non-trivial approximation factor known previously

Selection w/o Uniqueness Constraints • Can be seen as a “simultaneous set covering” problem: - The ground set is partitioned into n disjoint sets, each with 2 L elements - The goal is to select a minimum number of sets (== primers) that cover at least half of the elements in each partition • Naïve modifications of the greedy set cover algo do not work • Key idea: use potential function to measure progress towards fasibility. For primer selection, potential function counts the total number of elements that remain to be covered • Initially, = n. L • For feasible solutions, = 0

Greedy Approximation Algorithm Potential-Function Driven Greedy Algorithm m Select a primer that decreases potential function by the largest amount (breaking ties arbitrarily) m Repeat until feasibility is achieved m m Theorem: The greedy algorithm in returns a feasible primer set whose size is at most 1+ln ∆ times larger than the optimum, where ∆ is the maximum potential value decrease caused by a single primer For primer selection ∆ is equal to n. L in the worst case, and is much smaller in practice q The number of primers selected by the greedy algorithm is at most ln(n. L) larger than the optimum

Selection w/ Uniqueness Constraints m Can be modeled as minimum multicolored subgraph problem: add edge colored by color i between two primers if they amplify i-th target but do not amplify any other genomic sequence m Trivial approximation algorithm: select 2 primers for each amplification target q m O(n 1/2) approximation since at least n 1/2 primers required by every feasible solution No non-trivial approximation known previously

Integer Program Formulation q Variable xv for every graph node (candidate primer) v V; xv set to 1 if v is selected, and to 0 otherwise q Variable ye for every graph edge e E; ye set to 1 if corresponding primer pair selected to amplify one of the targets,

LP-Rounding Approximation LP-Rounding Algorithm • Solve linear programming relaxation • Select node u with probability xu m Theorem: With probability of at least 1/3, the number of selected nodes is within a factor of O(m 1/2 lnn) of the optimum, where m is the maximum number of edges sharing the same color and n is the number of nodes (candidate primers). m For primer selection, m L 2 approximation factor is O(Llnn)

Experimental Setting m SNP genotyping datasets q q m m Extracted from NCBI databases Randomly generated using uniform distribution C/C++ code, 2. 8 GHz Dell Power. Edge running Linux Compared algorithms q q G-FIX: greedy primer cover algorithm of Pearson et al. Ø Primers restricted to be within L/2 bases of amplified SNPs G-VAR: naïve modification of G-FIX Ø For each SNP, first selected primer can be up to L bases away from SNP Ø If first selected primer is L 1 bases away from the SNP, opposite sequence is truncated to a length of L- L 1 MIPS-PT: iterative beam-search heuristic of Souvenir et al. G-POT: potential function driven greedy algorithm

Experimental Results, NCBI tests G-FIX Targets MIPS-PT G-POT (Pearson et al. ) # G-VAR (G-FIX with dynamic truncation) (Souvenir et al. ) (Potentialfunction greedy) k #Primers 8 20 10 12 8 50 10 12 8 100 10 12 CPU sec #Primers CPU sec 7 9 14 13 23 31 17 37 53 0. 04 0. 03 0. 04 0. 13 0. 22 0. 14 0. 49 0. 37 0. 59 7 10 13 15 24 32 20 37 48 0. 08 0. 30 0. 36 0. 30 0. 89 0. 72 0. 84 8 13 18 21 30 41 32 50 75 10 15 26 48 150 246 226 844 2601 6 9 13 10 18 29 14 31 42 0. 10 0. 08 0. 11 0. 32 0. 33 0. 28 0. 58 0. 75 0. 61

Experimental Results, k=8

Experimental Results, k=12

Runtime, k=10

Ongoing Work on Primer Selection m Extending the greedy algorithm to degenerate primer selection q q m Incorporating improved hybridization models q m m Huge number of feasible candidate primers impractical to find primer with largest reduction in potential function The greedy algorithm remains provably good if only near-optimal choices are made in each step Allow hybridization with mismatches, enforce constraints on melting temperature, secondary structure, cross hybridization, etc. Closing gap between O(lnn) inapproximability bound and O(m 1/2 lnn) approximation factor for the minimum multi-colored subgraph problem Finding approximation algorithms and practical heuristics for partitioning into multiplexed PCR reactions (Aumann et al. WABI’ 03)

The String Barcoding Problem String barcoding is a pathogen identification technique recently proposed by Rash and Gusfield, and Bornemann et al. In this technique, a number of short oligos called distinguishers are spotted or synthesized on a microarray and hybridized with the fluorescently labeled DNA of unknown pathogens. The hybridization pattern can be viewed as a string of 0's and 1’s. The unknown pathogen can be identified by comparing this 0/1 pattern (its ``barcode'') with a set of pre-computed patterns for the pathogens. The main objective is to minimize the number of distinguishers needed to uniquely identify the pathogens. Given: Genomic sequences g 1, …, gn Find: Minimum number of strings t 1, …, tk Such that: For every gi gj, there exists a string tl which is the Watson. Crick complement for a substring of gi or gj, but not of both

Ongoing Work on String Barcoding m m m The greedy setcover algorithm, in which pairs of pathogens are viewed as elements to be covered, and candidate distinguishers are viewed as sets, is known to guarantee an approximation factor of 2 lnn An “information content” greedy algorithm was recently shown by Berman et al. to have an approximation factor of 1+lnn In ongoing work we explore heuristics for the following important extensions of the string barcoding problem: q q Probe mixtures as distinguishers. In spotted microarrays, it is feasible to spot a mixture consisting of a limited number of probes at any given array location. Using probe mixtures can reduce the number of spots on the array - hence barcode length - close to the information theoretical lower-bound of log 2 n Robust barcodes. Practical application of string barcoding is complicated by imperfect hybridization, experimental errors, and variability in pathogen genomic sequence. We are exploring robust barcodes using redundant distinguishers and error correcting schemes

References m m m R. J. Fernandes and S. S. Skiena. Microarray synthesis through multiple-use PCR primer design. Bioinformatics, 18: S 128–S 135, 2002. M. T. Hajiaghayi, K. Jain, K. M. Konwar, L. C. Lau, I. I. Mandoiu, A. C. Russell, A. A. Shvartsman, and V. V. Vazirani. The Minimum k-colored subgraph problem in haplotyping and DNA primer selection, submitted to ACM Symp. on Discrete Algorithms. K. M. Konwar, I. I. Mandoiu, A. C. Russell, and A. A. Shvartsman, Improved Algorithms for Minimum PCR Primer Set Selection with Amplification Length Constraints, submitted to 3 rd Asia Pacific Bioinformatics Conference. K. M. Konwar, I. I. Mandoiu, A. C. Russell, and A. A. Shvartsman, Approximation algorithms for minimum PCR primer set selection with amplification length and uniqueness constraints. ACM Computing Research Repository, Technical Report cs. DS/0406053, 2004. W. R. Pearson, G. Robins, D. E. Wrege, and T. Zhang. On the primer selection problem for polymerase chain reaction experiments. Discrete and Applied Mathematics, 71: 231– 246, 1996. R. Souvenir, J. Buhler, G. Stormo, and W. Zhang. Selecting degenerate multiplex PCR primers. In Proc. 3 rd Intl. Workshop on Algorithms in Bioinformatics (WABI), pages 512– 526, 2003.