MPprimer: a program for multiplex PCR primer design with high reliability

Wubin Qu , Zhiyong Shen and Chenggang Zhang .

About this program

MPprimer employs the widely used primer design program Primer3 [Rozen, et al. 2000] and primer specificity evaluation program MFEprimer [Qu, et al. 2009] to design and evaluate the candidate primers, followed with rigorous primer dimer examination to exclude unsuitable primers. A graph expanding algorithm derived from greedy algorithm was used to determine the optimal primer set combinations (PSCs) for multiplex PCR. In addition, MPprimer provides a virtual electrophotogram to help users choose the best PSC. In conclusion, MPprimer is a valuable tool for designing specific, no dimer formation, amplicons size constrained PSCs for multiplex PCR experiments.

Example 1: multiplex PCR primer design for five mouse genes using the MPprimer program ( Click here to run this example )

To examine the performance of MPprimer, we designed a primer set combination (10 primers) for multiplex PCR analysis of five mouse genes (beta-actin, B2m, Pgk1, Rpl13a and GAPDH). Firstly, input the five mouse genes sequences in the input box. Secondly, specify the product size for each sequences. Here, we specify the 851-1000, 650-720, 400-520, 101-200 and 201-300 for the five mouse genes ( Fig. 1 ). Thirdly, set the "primer size" as following: Min: 20 Opt: 25 Max: 30. Other settings use the default values. At last, click "Pick primers" to submit. Some seconds later, MPprimer will output the results as Fig. 2 and Fig. 3 shows. Fig. 2 shows the summary results of 15 PSCs. We select the fifth PSC for experimental validation. Fig. 3 is the detailed information of PSC 5, including primer sequence, Tm, product size, etc. Fig. 4 is the virtual electrophotogram of the PSC 5. Fig. 5 is the real electrophotogram of our experimental results. Fig. 4 perfectly predict the mobility of amplicons of PSC5 as Fig. 5 showed.

Fig. 1 Specify the product size for each sequence


Fig. 2 PSC 5 was selected for experimental validation

Fig. 3 Primer details of PSC 5 . Select all the primer sets for virtual electrophoresis

Fig. 4 Virtual electrophoresis for five primer sets of PSC 5

Fig. 5 Experimental validation of PSC 5. The electrophotogram was matched perfectly with the virtual electrophotogram. Lane 0 and Lane 9 show the DNA marker of DL2000 (TaKaRa); PCR results of DL2000 size marker, Lane 1 to 5 shows the results of standard PCR with one primer set (not primer set combination) in the PSC list; Lane 6 shows the multiplex PCR results with equal mixture of the five primer set combination. Lane 7 shows the results of multiplex PCR with partial combination (two primer sets) of the five primer sets and Lane 8 shows the results of multiplex PCR with three primer sets combinations. All products had the expected length, and there was no visible unspecific amplicons. The conditions for all of the PCR analyses are the same as: first denaturing at 94 °C for 5 min; 36 cycles for denaturing at 94 °C for 30 s; annealing at 58 °C for 45 s; extension at 72 °C for 50 s; final extension at 72 °C for 5 min.

Example 1: multiplex PCR primer design for detecing exon deletions of 79 exons of DMD gene [Run examples: tube 1, tube 2, tube 3, tube 4, tube 5].

Designing the multiplex PCR primers for the DMD gene is a good example to evaluate the performance of MPprimer. For the 79 exons of DMD gene, we have successfully used MPprimer to design the multiplex PCR primers for 20×, 20×, 20×, 14×, and 5× plex PCR reactions in five tubes in one reaction. For DMD exon deletion test, users need to utilize high-resolution techniques such as SDS-PAGE or HPLC to separate the amplicons of nearly close size. In this case, the parameter MinBS (for minimal binding space) is not the critical factor anymore for consideration. Therefore, for designing multiplex PCR primers for the exons of DMD gene, we use the MinBS value of “10 bp” to replace the default “2 mm” mobility distance which was usually used in the agarose gel electrophoresis. To simplify the primer design in MPprimer, we have now added the MinBS option in the home page of MPprimer web service . We thank the review for this good suggestion because passing this examination has made the MPprimer service more flexible and useful.

Additionally, the settings of “ PRODUCT SIZE RANGE ” are very important in multiplex PCR primer design in MPprimer, especially for this kind of 20× plex PCR primers design. In other words, if the user set the suitable “ PRODUCT SIZE RANGE ”, they can successfully design reliable multiplex PCR primers with MPprimer. Although MPprimer can automatically assign “ PRODUCT SIZE RANGE ” for each of the template sequences, it still needs users' modification in some cases. We also add this information in the “ User guide ” page of MPprimer web service.

Following is the primer design process of 79 exons of the DMD gene:

  1. Download the exon sequence from Ensembl web site ( http://www.ensembl.org/Homo_sapiens/Transcript/Exons?db=core;g=ENSG00000198947;r=X:31132808-33357558;t=ENST00000357033 ).
  2. Preparing DNA template sequences of the DMD exons. To detect the deletion of DMD exons, we need to pick up the primers located at the flanking sequence of the exon. Therefore, both of the intron sequences of upstream (+500 bp) and downstream (-500 bp) of each DMD exon are prepared ( Fig. 1 ). For the introns which sequence length is shorter than 500 bp such as intron 14-15 (107 bp) , 31-32 (396 bp) , 35-36 (309 bp), we directly used the whole intron sequence.
  3. Set the parameters as shown in Fig. 2 .
  4. Specific the parameters “Target region” and “ PRODUCT SIZE RANGE ” as shown in Fig. 3 .
  5. Run MPprimer and get the PSCs. Fig. 4 -7 are the output pages for designing multiplex PCR primers for exons of the DMD gene.

 

Fig. 1 . Preparation of DNA template sequence of exons of the DMD gene for MPprimer design

Fig. 2 Set the parameters for DNA template sequences of exons of the DMD gene

Fig. 3 Specific the parameters “Target region” and “ PRODUCT SIZE RANGE ”

Fig. 4 . MPprimer output for exon templates of DMD gene (1)

Fig. 5 MPprimer output for exon templates of DMD gene (2)

Fig. 6 MPprimer output for exon templates of DMD gene (3)

Fig. 7 MPprimer output for exon templates of DMD gene (4)