# Designing multiple degenerate primers via consecutive pairwise alignments

- Hamed Shateri Najafabadi
^{1, 2}, - Noorossadat Torabi
^{1, 3}and - Mahmood Chamankhah
^{4}Email author

**9**:55

**DOI: **10.1186/1471-2105-9-55

© Najafabadi et al; licensee BioMed Central Ltd. 2008

**Received: **04 April 2007

**Accepted: **27 January 2008

**Published: **27 January 2008

## Abstract

### Background

Different algorithms have been proposed to solve various versions of degenerate primer design problem. For one of the most general cases, multiple degenerate primer design problem, very few algorithms exist, none of them satisfying the criterion of designing low number of primers that cover high number of sequences. Besides, the present algorithms require high computation capacity and running time.

### Results

PAMPS, the method presented in this work, usually results in a 30% reduction in the number of degenerate primers required to cover all sequences, compared to the previous algorithms. In addition, PAMPS runs up to 3500 times faster.

### Conclusion

Due to small running time, using PAMPS allows designing degenerate primers for huge numbers of sequences. In addition, it results in fewer primers which reduces the synthesis costs and improves the amplification sensitivity.

## Background

Polymerase Chain Reaction, or PCR [1], is a ubiquitous technique which amplifies a specific region of DNA, so that enough copies of that region is available to be adequately tested, sequenced or manipulated in other fashions. In order to use PCR, one must know the exact sequences which lie on either side of the DNA region of interest. These sequences are used to design two synthetic DNA oligonucleotides, or *primers*, one complementary to each strand of the DNA double-helix and lying on opposite sides of the target region. The primers are typically 20–30 nucleotides long.

*S*=

*x*

_{1}

*x*

_{2}...

*x*

_{ l }, where

*x*

_{ i }⊆ ∑,

*x*

_{ i }≠ Ø and

*l*is the length of

*S*. A sequence is called

*degenerate*if some of its positions have several possible bases [3]. For example, in the primer

*P** = {G}{G}{C,G}{A}{T,C,G}{A} the third position is C or G and the fifth is C, T or G. The IUPAC illustration of

*P** will be GGSABA (Figure 1). The

*degeneracy*of a sequence is the number of unique sequence combinations it contains, which can be calculated as

*d*(

*S*) = Π

^{ l }

_{i=1}|

*x*

_{ i }|. For example,

*d*(

*P**) = 1 × 1 × 2 × 1 × 3 × 1 = 6. Degenerate primers are useful for amplifying several related genomic or cDNA sequences, and have been exploited in various applications such as amplifying DNA sequences of homologous genes or genes from a particular protein family and analysis of species diversity [4–6].

Traditionally, degenerate primers were designed manually by examining multiple alignments of the target sequences. However, several programs are now available for designing degenerate primers for aligned sequences. CODEHOP [7] and DePiCt [8] are programs for designing degenerate primers for aligned protein sequences in order to identify new members of protein families. For each given multiple sequence alignment, CODEHOP constructs a pair of primers. Each primer consists of a degenerate 3' core region, typically with degeneracy of at most 128, and a 5' consensus sequence that stabilizes annealing. It works well for small sets of proteins, taking into account the codon usage of the target genome as well as the desired annealing temperature. However, it is inappropriate for constructing primers with high degeneracy on large sets of long genomic sequences. DePiCt clusters the sequences using a simple similarity score and then designs a pair of primers for each cluster by translating conserved blocks of amino acids into nucleotides.

*P*=

*p*

_{1}

*p*

_{2}...

*p*

_{ l }covers the sequence

*S*if there is a substring

*F*of length

*l*in

*S*where for each character

*f*

_{ i }in

*F*,

*f*

_{ i }⊆

*p*

_{ i }). On the other hand, in order to reduce the probability of amplifying unrelated sequences, the degeneracy must be bounded. This contradictory nature of the degenerate primer design (DPD) problem has led to definition of several variants of this problem, all of which are NP-complete:

- 1.
Maximum Coverage Degenerate Primer Design (MC-DPD) tries to find a primer of length

*l*and degeneracy at most*d*_{ max }that covers a maximum number of strings (sequences) of a given input set, each of length*l*. HYDEN [9], an algorithm based on a heuristic approach, basically addresses this variant of DPD problem and was first used to design degenerate primers for a set of genomic sequences in order to find new human olfactory receptor genes [9, 10]. - 2.
Minimum Degeneracy Degenerate Primer Design (MD-DPD) addresses the problem of finding a primer of length

*l*and minimum degeneracy that covers all the input strings, each of which having a length equal to or greater than*l*. - 3.
Minimum Primers Degenerate Primer Design (MP-DPD) is applied when a set of strings of length

*l*is given, and finds a minimum number of primers of length*l*and degeneracy at most*d*_{ max }, so that each input string is covered by at least one primer.

MP-DPD has the constraint that all input sequences are of the same length as the primers, which is not the case for most real situations. Removing this constraint, i.e. allowing the strings to have arbitrary lengths, results in a more general problem, Multiple Degenerate Primer Design (MDPD) [2]. MDPD is to find a minimum number of primers of length at least *l*_{
min
}and degeneracy at most *d*_{
max
}, given a set of *n* strings of various lengths (equal to or greater than *l*_{
min
}), so that each input string is covered by at least one primer. A currently available algorithm for designing multiple degenerate primers, called PT-MIPS [2], has been developed in the context of SNP genotyping. It uses an iterative beam-search technique to construct progressively a set of primers until all sequences are covered.

In this work, we introduce a new algorithm for solving MDPD problems which consecutively uses an *ad hoc* p airwise a lignment for m ultiple p rimer s election – hence called PAMPS. We will show that PAMPS performs better than previous algorithms on different sets of input strings, i.e. results in smaller number of primers in a considerably less computation time.

## Results and Discussion

*l*

_{ min }= 15 and

*d*

_{ max }= 10

^{4}. Almost always PAMPS resulted in smaller primer sets. Only in few cases both PAMPS and PT-MIPS produced primer sets with equal sizes. To compare PAMPS and PT-MIPS quantitatively, we defined efficiency of PAMPS as

*m*

_{ MIPS }and

*m*

_{ PAMPS }represent the number of primers designed by PT-MIPS and PAMPS, respectively. Figure 3 illustrates the values of

*f*

_{ PAMPS }for different numbers of sequences and different sequence lengths. Obviously PAMPS outperforms PT-MIPS, especially when smaller sets of sequences or long sequences are used. In most situations, PAMPS decreases the number of final primers by 30%–35%. PAMPS outperforms PT-MIPS for a wide range of primer sizes and maximum degeneracy values (Figure 4).

^{®}CPU): solving MDPD problem for 100 input sequences of length 50 nucleotides is about 3380 times faster using PAMPS compared to PT-MIPS (Figure 5). This allows PAMPS to be used to design highly degenerate primers for thousands of input sequences each hundreds of nucleotides long. Hence, even though the number of designed primers using PAMPS and PT-MIPS may converge as the number of input sequences increases, considering computation time strongly encourages using PAMPS; for an input set of 10

^{4}random sequences of length 2000, PAMPS needs an average time of 228 seconds to complete the computations on a 2.0 GHz Intel

^{®}Core™ 2 CPU. We should mention that PT-MIPS did not yield in any results after three days of running the same job as PAMPS. Based on previous comparisons of PAMPS and PT-MIPS, we can estimate that for PT-MIPS it takes more than nine days to finish a job like this.

PT-MIPS [2] is previously compared with HYDEN [9]. Though HYDEN is basically designed to solve MC-DPD problems, it can be used iteratively to approximate MDPD problems, i.e. once a primer of length *l*_{
min
}and degeneracy at most *d*_{
max
}is found that covers the maximum number of input sequences, the sequences which are covered by this primer are subtracted from input set and HYDEN is run again on the remaining sequences. By repeating this procedure, eventually a set of primers is obtained which covers all sequences. Since it has been shown that PT-MIPS outperforms HYDEN [2] and as PAMPS outperforms PT-MIPS, we avoided the direct comparison of PAMPS and HYDEN.

*l*

_{ min }(Figure 6). Therefore, any subsequences of length

*l*

_{ min }from each output can be selected to be used for PCR amplification. If the longest possible PCR product is desired, then the very upstream subsequence should be used. However, for most PCR reactions it is important to have primers with similar annealing temperatures if a mixture of primers is used. Since different combinations of primers can be chosen, it is possible to select the primers that have similar annealing temperatures. PAMPS is accompanied by a simple iterative algorithm provided in a separate software that chooses the best combination of primers in order to achieve the minimum variance among primer annealing temperatures. Primer annealing temperatures are estimated as [11].

## Conclusion

In this work we presented a new algorithm, called PAMPS, for solving MDPD problems. PAMPS exploits an altered pairwise alignment to select the subsequences which may be merged into degenerate primers. PAMPS was shown to run significantly faster than a previously developed software, PT-MIPS [2] and also gives better results (i.e. smaller sets of primers), reducing the synthesis costs of primers. Besides, when the number of mixed primers that are used in a PCR reaction are decreased, the concentration of the reacting primer increases, which usually improves the sensitivity of amplification. PAMPS, in contrast to previous algorithms, does not restrict the output to the exact primer length that was given; instead, it may result in primers longer than the requested length which allows selecting an appropriate primer in terms of annealing temperature. PAMPS can be used to design degenerate primers for amplification of genes with uncertain sequences, such as new members of gene families or libraries of antibody variable fragments. An implementation of PAMPS is provided in the Additional file 1.

## Methods

### Merging two aligned sequences

*S*

_{1}=

*x*

_{1}

*x*

_{2}...

*x*

_{ l }and

*S*

_{2}=

*y*

_{1}

*y*

_{2}...

*y*

_{ l }results in

*S*

_{1,2}= (

*x*

_{1}∪

*y*

_{1})(

*x*

_{2}∪

*y*

_{2})...(

*x*

_{ l }∪

*y*

_{ l }) (Figure 7). Obviously, the regions of each sequence that are located in a gap are of no value in designing a degenerate primer that can cover both sequences. Therefore, these regions should be removed and the two regions surrounding each gap should be joined, at a point that is referred to as a "split point" through this article. Obviously, a degenerate primer that covers both sequences is located between two split points.

*d*

_{ max }and those having lengths smaller than

*l*

_{ min }. To achieve this goal, we only retain those nucleotides that are located within at least one window with length

*l*

_{ min }and degeneracy at most

*d*

_{ max }. Obviously, this window cannot have a split point within. If no such a window could be found for a nucleotide, that nucleotide should be removed. This results in the removal of all nucleotides between two split points that are closer than

*l*

_{ min }. The remaining regions are joined together with a new split point (Figure 8).

### Alignment

The alignment algorithm that is used by PAMPS is very similar to the conventional global alignment [12]. However, the scoring methods differ in some details. Since the purpose is to achieve an alignment that results in a merged sequence with low degeneracy, we defined the score of each match/mismatch as*M*(*x*, *y*) = 2 - log_{2} |*x* ∪ *y*|.

*x*,

*y*⊆ ∑. The two sequences that are being aligned may contain some split points as they may themselves have been resulted from merging other sequences. In this case, passing over a split point has the same penalty as gap opening (see Figure 9). In this work, we set the penalty of gap opening to -10.0 and gap extension to 0.0, since for our purpose it is of no importance how long a gap is. Merging two split sequences causes all split points to be copied into the relevant positions in the merged sequence. After two aligned sequences are merged and refined (Figure 8), the alignment score is recalculated, since some of the portions that are scored in the original alignment may be removed in the refined sequence.

### Designing degenerate primers

In order to design degenerate primers, pairs of sequences should be aligned and merged consecutively until no more sequences could be merged (i.e. merging any more pairs of sequences results in primers either with lengths less than *l*_{
min
}or with degeneracy more than *d*_{
max
}). However, there are different combinations in which sequences can be merged, each of which may result in a different set of primers. The optimum set is the one that contains the least number of primers. PAMPS uses a procedure similar to MIPS [2] to search for the optimum set of primers:

*P*

_{1},

*P*

_{2}, ...,

*P*

_{ m }} covers the set S = {

*S*

_{1},

*S*

_{2}, ...,

*S*

_{ n }} (P covers S if for each

*S*

_{ j }∈ S, 1 ≤

*j*≤

*n*there is a

*P*

_{ i }∈ P, 1 ≤

*i*≤

*m*which covers

*S*

_{ j }). For the

*S*

_{n+1}to be covered by P,

*m*+1 "actions" are possible: merging

*S*

_{n+1}with

*P*

_{ i }(1 ≤

*i*≤

*m*), or adding a new primer (

*P*

_{m+1}, which is the same as

*S*

_{n+1}) to P. Thus, PAMPS starts with P = {

*P*

_{1}},

*P*

_{1}=

*S*

_{1}, and either merges

*S*

_{2}with

*P*

_{1}or adds it to P as

*P*

_{2}. P is expanded until it covers all sequences. If in any step one of the requirements of MDPD is not fulfilled, i.e. the length of a primer becomes less that

*l*

_{ min }or the degeneracy of a primer exceeds

*d*

_{ max }, PAMPS backtracks to a previous P and continues with another "action" (Figure 10). PAMPS searches for all P's each of which covering all sequences, and chooses the one with the minimum |P| (i.e. chooses the P that covers all sequences with the minimum possible number of primers); however, the minimum |P| is guaranteed only if the following conditions are met; (1) no heuristic approach is employed; (2) no gap is allowed in alignment of sequences, i.e. penalty of gap opening is -8; (3) length of each sequence is equal to required primer length which turns the problem into MP-DPD. Finding the minimum |P| is simply a result of searching all combination of actions, which is obviously not possible for large sets of input sequences; hence the need for a heuristic approach is emerging.

PAMPS uses a similar heuristic approach as MIPS [2] to reduce the search space. Assume P_{1} is a previously found set of primers that contains *m* primers and covers *n* sequences, and P_{2} is a newly found set that also contains *m* primers and covers *n* sequences. P_{2} is only expanded if the sum of scores of its primers (see section Alignment) exceeds that of P_{1} (Figure 10).

## Declarations

### Acknowledgements

Authors are grateful to Manely Rashedan for her everlasting supports. This project is granted by Avesina Research Institute, Tehran. HSN and NT are also supported by University of Tehran.

## Authors’ Affiliations

## References

- Mullis K, Faloona F, Scharf S, Saiki R, Horn G, Erlich H: Specific enzymatic amplification of DNA
*in vitro*: The polymerase chain reaction.*Cold Spring Harbor Symp Quant Biol*1986, 51: 263–273.View ArticlePubMedGoogle Scholar - Souvenir R, Buhler J, Stormo G, Zhang W: Selecting degenerate multiplex PCR primers.
*Proceedings of the 3rd Workshop on Algorithms in Bioinformatics (WABI 2003)*2003, 512–526.View ArticleGoogle Scholar - Kwok S, Chang S, Sninsky J, Wang A: A guide to the design and use of mismatched and degenerate primers.
*PCR Methods Appl*1994, 3: S39-S47.View ArticlePubMedGoogle Scholar - Fuchs T, Malecova B, Linhart C, Sharan R, Khen M, Herwig R, Shmulevich D, Elkon R, Steinfath M, O'Brien JK, Radelof U, Lehrach H, Lancet D, Shamir R: DEFOG: A Practical Scheme for Deciphering Families of Genes.
*Genomics*2002, 80(3):1–8. 10.1006/geno.2002.6830View ArticleGoogle Scholar - Jarman SN: Amplicon: software for designing PCR primers on aligned DNA sequences.
*Bioinformatics*2007, 20(10):1644–1645. 10.1093/bioinformatics/bth121View ArticleGoogle Scholar - Jarman SN, Deagle BE, Gales NJ: Group-specific PCR for DNA-based analysis of species diversity and identity in dietary samples.
*Mol Ecol*2003, 13: 1313–1322. 10.1111/j.1365-294X.2004.02109.xView ArticleGoogle Scholar - Rose T, Schultz E, Henikoff J, Pietrokovski S, McCallum C, Henikoff S: Consensus-degenerate hybrid oligonucleotide primers for amplification of distantly related sequences.
*Nucleic Acids Res*1998, 26: 1628–1635. 10.1093/nar/26.7.1628PubMed CentralView ArticlePubMedGoogle Scholar - Wei X, Kuhn D, Narasimhan G: Degenerate primer design via clustering.
*Proceedings of the 2nd IEEE Computer Society Bioinformatics Conference (CSB 2003)*2003, 75–83.View ArticleGoogle Scholar - Linhart C, Shamir R: The degenerate primer design problem.
*Bioinformatics*2002, 18: S172-S180.View ArticlePubMedGoogle Scholar - Linhart C, Shamir R: The Degenerate Primer Design Problem: Theory and Applications.
*J Comput Biol*2005, 12: 431–456. 10.1089/cmb.2005.12.431View ArticlePubMedGoogle Scholar - Breslauer KJ, Frank R, Blocker H, Marky LA: Predicting DNA duplex stability from the base sequence.
*Proc Natl Acad Sci USA*1986, 83: 3746–3750. 10.1073/pnas.83.11.3746PubMed CentralView ArticlePubMedGoogle Scholar - Needleman SB, Wunsch CD: A general method applicable to the search for similarities in the amino acid sequence of two proteins.
*J Mol Biol*1970, 48: 443–453. 10.1016/0022-2836(70)90057-4View ArticlePubMedGoogle Scholar - Cornish-Bowden A: IUPAC-IUB symbols for nucleotide nomenclature.
*Nucleic Acids Res*1985, 13: 3021–3030. 10.1093/nar/13.9.3021PubMed CentralView ArticlePubMedGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.