Bioinformatics Advance Access originally published online on December 14, 2004
Bioinformatics 2005 21(8):1701-1702; doi:10.1093/bioinformatics/bti219
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
MultiPLX: automatic grouping and evaluation of PCR primers
1Department of Bioinformatics, University of Tartu Tartu, Estonia
2Estonian Biocentre Riia str. 23, Tartu 51010, Estonia
*To whom correspondence should be addressed.
| Abstract |
|---|
|
|
|---|
Summary: MultiPLX is a new program for automatic grouping of PCR primers. It can use many different parameters to estimate the compatibility of primers, such as primerprimer interactions, primerproduct interactions, difference in melting temperatures, difference in product length and the risk of generating alternative products from the template. A unique feature of the MultiPLX is the ability to perform automatic grouping of large number (thousands) of primer pairs.
Availability: Binaries for Windows, Linux and Solaris are available from http://bioinfo.ebc.ee/download/. A graphical version with limited capabilities can be used through a web interface at http://bioinfo.ebc.ee/multiplx/. The source code of the program is available on request for academic users.
Contact: maido.remm{at}ut.ee
| INTRODUCTION |
|---|
|
|
|---|
Polymerase chain reaction (PCR) is widely used in different areas of science and biotechnology. Grouping of PCR primer sets for amplification in a single tube (multiplexing) provides substantial savings in terms of time, chemicals and, most importantly sample materials. Thus, multiplexing is a powerful way for optimizing the cost of genetic analysis and some procedures of PCR are patented (Piepenbrock et al., 2004), which further demonstrates its importance. While a lot of effort has been put into the experimental optimization of multiplex PCR conditions (Henegariu et al., 1997; Zangenberg et al., 1999), little attention has been paid to the prediction of primer compatibility in multiplex reactions. Smaller multiplex groups are often composed manually and verified by trial and error, but this kind of approach is suboptimal for large datasets. Although there are software programs available for testing the compatibility of PCR primer pairs for multiplexing (Rychlik, 1995; Singh et al., 1998; Vallone and Butler, 2004), these do not perform automatic grouping and are unable to handle large datasets automatically. We have developed a program, MultiPLX, for automatic primer compatibility testing and grouping. It can also be used for the compatibility evaluation of the existing primer groups. The program is written in standardized C++ for maximum performance and compatibility. We have built binaries for Linux (x86 compatible processors), Windows and Solaris operating systems. Most of the functionality can be used through the web-based interface also.
| ALGORITHM AND IMPLEMENTATION |
|---|
|
|
|---|
The workflow of MultiPLX is divided into two main tasks: (1) calculation of compatibility scores for existing PCR primers and (2) grouping of PCR primers based on the cut-off scores. By scores, we are referring to all measurable parameters that may affect the compatibility of different PCR primer sets. Currently, MultiPLX is capable of calculating eight different pre-defined score types:
- maximum binding energy (deltaG) of two primers including 3' ends of both primers
- maximum binding energy of 3' end of one primer with any region of another primer
- maximum binding energy of any region of different primers
- maximum binding energy of 3' end of one primer with any region of PCR product
- maximum binding energy of any region of primer with any region of PCR product
- maximum product length difference between compared PCR primer sets
- minimum product length difference between compared PCR primer sets
- maximum difference in primer melting temperatures between compared PCR primer sets.
Compatibility scores are calculated for all possible pairwise combinations of PCR primer sets. MultiPLX allows in-depth examination of all default scores, including the listing of alignments between primers and their products from different PCR sets. Each score type has a default cut-off value to determine whether two given PCR primer sets are compatible and thus can be amplified in the same group. Compatibility cut-off values for all scores can be adjusted by the user either as generic stringency criteria (low, normal, high) or as exact numerical values. Default compatibility cut-offs are generating approximately N3 plex (high stringency), 610 plex (normal stringency) or 1520 plex (low stringency) groups of 100 PCR primer sets.
Automatic grouping of primers is based on compatibility scores and corresponding cut-off values, mentioned above. Two primer sets are considered incompatible if at least one of the compatibility scores exceeds the predefined cut-off value. Two grouping algorithms are implemented.
Algorithm 1. Number of friends (default).
- Sort PCR primer sets by the number of compatible primer pairs they have.
- In a sorted list of primers, try to insert each primer set into existing groups. If this is not possible, create a new group and move the primer set there.
- Repeat until all PCR primer sets have been grouped.
- Report the resultlist of primer sets associated with their group numbers.
Algorithm 2. Random grouping.
The second algorithm is based on random grouping.
- Reorder PCR primer sets randomly.
- Try to group each primer set into existing groups. If this is not possible, create a new group and move the primer set there.
- Repeat until all PCR primer sets are grouped.
- Record the number of groups and group number for each primer.
- Repeat the steps 14 N times (default value of N is 10 000).
- Report the result that generated the lowest number of groups.
The random grouping approach typically generates groups of very different sizes. To make the number of primers in each group more uniform an extra optimization step can be invoked. Optimization tries to move PCR sets from larger groups into smaller ones and swap elements randomly between the groups, until all group sizes differ by no more than one primer pair. None of the grouping algorithms guarantees the optimal solution (lowest possible number of groups). However, the first algorithm typically gives smaller number of groups with shorter computation time. Randomized grouping may be desirable in situations where non-deterministic grouping solution is preferred.
MultiPLX can also be used for the evaluation of existing PCR groups against the calculated score values. Groups where some primer pairs break cut-off rules will be listed alongside the problematic primers. Using this list it is possible to examine, which primers do not fit together and need to be replaced. This option can be helpful for the evaluation of the existing PCR groups with problematic results.
The computation time required by the program for different numbers of primers is shown in Table 1.
|
| Acknowledgments |
|---|
This work was supported by the Estonian Ministry of Education and Research grant no. 0182649s04.
Received on October 11, 2004; revised on November 25, 2004; accepted on December 10, 2004
| REFERENCES |
|---|
|
|
|---|
Henegariu, O., Heerema, N.A., Dlouhy, S.R., Vance, G.H., Vogt, P.H. (1997) Multiplex PCR: critical parameters and step-by-step protocol. Biotechniques, 23, 504511[Web of Science][Medline].
Kaderali, L. and Schliep, A. (2002) Selecting signature oligonucleotides to identify organisms using DNA arrays. Bioinformatics, 18, 13401349
Piepenbrock, C., Adorjan, P., Rujan, T., Schmitt, A. (2004) Method for amplification of nucleic acids of low complexity. European patent WO2004015139.
Rychlik, W. (1995) Selection of primers for polymerase chain reaction. Mol. Biotechnol., 3, 129134[Web of Science][Medline].
Singh, V.K., Mangalam, A.K., Dwivedi, S., Naik, S. (1998) Primer premier: program for design of degenerate primers from a protein sequence. Biotechniques, 24, 318319[Medline].
Vallone, P.M. and Butler, J.M. (2004) AutoDimer: a screening tool for primer-dimer and hairpin structures. Biotechniques, 37, 226231[Medline].
Zangenberg, G., Saiki, R.K., Reynolds, R. (1999) Multiplex PCR: Optimization guidelines. In Innis, M., Gelfand, G.H., Sninsky, J.J. (Eds.). PCR Application: Protocols for Functional Genomics, , London Academic Press, pp. 7394.
This article has been cited by other articles:
![]() |
R. Andreson, T. Mols, and M. Remm Predicting failure rate of PCR in large genomes Nucleic Acids Res., June 1, 2008; 36(11): e66 - e66. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
