TipMT: Identification of PCR-based taxon-specific markers
© The Author(s). 2017
Received: 16 June 2016
Accepted: 11 January 2017
Published: 11 February 2017
Molecular genetic markers are one of the most informative and widely used genome features in clinical and environmental diagnostic studies. A polymerase chain reaction (PCR)-based molecular marker is very attractive because it is suitable to high throughput automation and confers high specificity. However, the design of taxon-specific primers may be difficult and time consuming due to the need to identify appropriate genomic regions for annealing primers and to evaluate primer specificity.
Here, we report the development of a Tool for Identification of Primers for Multiple Taxa (TipMT), which is a web application to search and design primers for genotyping based on genomic data. The tool identifies and targets single sequence repeats (SSR) or orthologous/taxa-specific genes for genotyping using Multiplex PCR. This pipeline was applied to the genomes of four species of Leishmania (L. amazonensis, L. braziliensis, L. infantum and L. major) and validated by PCR using artificial genomic DNA mixtures of the Leishmania species as templates. This experimental validation demonstrates the reliability of TipMT because amplification profiles showed discrimination of genomic DNA samples from Leishmania species.
The TipMT web tool allows for large-scale identification and design of taxon-specific primers and is freely available to the scientific community at http://188.8.131.52/tipMT/.
KeywordsMolecular marker Specific primers PCR PCR Multiplex Web application
Polymerase chain reaction (PCR)-based typing methods are molecular diagnostic techniques widely used in biological and biomedical studies. The level of discriminatory power of PCR-based typing depends upon the molecular marker targeted. Therefore, identifying appropriate DNA target regions for primer annealing is a crucial step because these regions must be conserved within the target taxa but must vary among related taxa [1, 2].
Recent advances in next-generation sequencing technology are enabling genome sequencing projects at a significantly lower cost, even for non-model organisms. The resulting increase in the amount of genomic data available, combined with bioinformatics tools, have led to the identification of highly informative markers, such as microsatellites and orthologous or taxa-specific genes .
Microsatellites or single sequence repeats (SSR) are tandem repeated stretches of short nucleotide motifs, usually ranging from 1 to 6 bp, ubiquitously distributed in the genomes of eukaryotic organisms. These regions are more prone to genetic variation and the differences in the length of individual SSR loci can be easily screened by PCR. In fact, this technique has been useful for several studies including strain typing and population genetics [4, 5]. The conventional method of SSR discovery is time consuming and costly. Therefore, in silico mining analysis has been used to improve marker identification [6, 7].
Orthologs are homologous proteins in different species that evolved from a single ancestral sequence and are related by speciation events. These sequences tend to show more functional similarity than other homologs. The identification of orthologous genes is useful in a wide range of contexts, such as inference of gene function, comparative genomics, evolutionary conservation and sequence variability . Due to their importance, many tools have been developed to predict ortholog groups, including the widely used software OrthoMCL .
The demand is increasing for bioinformatic tools that automate analysis of genomic data generated by next-generation sequencing technology . An example is the development of automated procedures to facilitate species-specific primer design for diagnostic methods . Several web-based tools for facilitating primer design are available , but many of them are written mainly to assist in the primer design process and are not meant to search for targets and analyze primer specificity. The use of fully automated methods to search for molecular markers and the availability of genomic data for a growing number of taxa would increase the efficiency of PCR-based genotyping applications. Moreover, this strategy might save time and resources because the in silico evaluation of the candidate primers against the target genomic sequence are performed prior to testing them in the laboratory . Thus, there is a need for a tool to search for appropriate genomic target regions and then design specific primers towards the selected markers.
In this context, we have developed TipMT to meet the growing demand for easy-to-use software that facilitates the design of primers that target molecular markers to distinguish the genomic sequences of different taxa. This program only requires genomic sequences of a target species and offers, as an output, specific primers for a given taxa.
The aim of the software pipeline is to provide a set of primer pairs flanking polymorphic sequences to identify taxa among related species using PCR, given their genomic sequences. By taking advantage of sequence data from related species, the pipeline identifies orthologous and singleton (taxa-specific) genes or SSR regions as target sequences that are likely to identify unique taxon or taxa. Because primer specificity is a key step in a PCR reaction, the pipeline identifies all potential annealing sites for the primers selected based on alignment and thermodynamics. Then, the TipMT evaluates compatibility among specific primers for designing multiplex PCR reactions. Finally, the program generates a virtual gel with the result of a simulated standard or multiplex PCR assay, where taxa can be identified by the size variation of the predicted PCR products.
User Input Data
TipMT is flexible because it accepts three different input data as target sequences: genomes, predicted genes and user defined target sequences. The user may provide sequences by uploading files in FASTA format or entering Nucleotide Accession numbers, and then, corresponding sequences will be downloaded from the NCBI RefSeq database. There are two types of genomic sequences: target taxa and cross-reaction taxa. The former sequences will be used in the pipeline as templates for primer design and to check the specificity of the primers. The latter are sequences of species that should not cross-react during the PCR assays. Regions in the target taxa that have similarities to sequences in the cross-reaction taxa are not targeted during the primer-designing step. The number of taxa analyzed in TipMT is not limited, but the processing time is substantially increased for each taxon added. For example, in the test using Leishmania braziliensis and Leishmania infantum genomes as templates on “Ortholog target” mode run on a Intel(R) Xeon(R) X3430 2.40GHz, 4 CPUs and 8 GB RAM machine, the software spent 12 h, while the execution with Leishmania amazonensis, L. braziliensis and L. infantum genomes in the same mode run and machine spent 22 h.
The next step is the definition of the target sequences by the user, and each approach offers advantages. SSRs are highly polymorphic and tend to be conserved between closely related species . These features enhance the success rate in the search for taxon specific primers. Moreover, SSR targets only require genomic sequences. Orthologs sequences are less error-prone than repetitive regions during genome assembly, thus the amplification failure tends to be lower compared to SSRs. Additionally, conserved primers among closely related taxa are easily found using this approach and singleton sequences enable the search for taxon-specific primers. In the case orthologs target sequences, the user should provide the sequences of the predicted genes in the genome. Also, users can have their custom sequences of interest be targets to design specific primers.
If the user chooses an SSR as a target, the pipeline will search SSR regions in the target genomes, and these regions will be the template sequences in the primer design step. On the other hand, if the user selects ‘ortholog target’, orthologous and singleton genes will be identified in the predicted coding regions provided by the user. Finally, the user may choose to provide target sequences instead of using the TipMT target search mechanism.
Lastly, the following primer design constraint parameters are defined: 1) mismatch and gap tolerance (both parameters are known to effect PCR specificity ); 2) PCR product sizes; 3) the number of primers per target (high values increase the chance to find specific primers but also increase the processing time); and 4) minimal difference (this value increases the PCR product size between taxa).
Target search. One way to improve sensitivity in PCR is to find the most appropriate template region for primer design. Ortholog and singleton sequences are identified in the predicted coding sequences provided by the user using OrthoMCL with default parameters values. ProGeRF searches for SSR regions, without degeneration or gaps in the sequence (perfect repeats).
Mask similar regions. Conservation of the flanking regions of the target sequences is essential for a high quality PCR assay because a high number of primer annealing sites can cause failure of the PCR assay . Similar regions between target sequences and cross-reaction genomes are identified using a MEGABLAST search with default parameters. MEGABLAST was chosen due to its speed and its ability to handle slight differences in genomic sequences. Next, regions with more than 95% of identity are masked with lowercase nucleotides (initially, all sequences are set as uppercase in the database).
Primer design. Candidate primers are generated for each of the DNA template sequences using primer3 2.3.5 with default parameters values and an option that rejects the primer candidates with lowercase letters in the first 3’ end position. Because the high number of annealing sites influences primer specificity, this procedure decreases the rate of low-success specific primers, avoiding primer design in regions with similarity between species .
Specificity check. All candidate pairs of primers generated are evaluated for specificity using the alignment based e-PCR algorithm and by thermodynamic properties using MFEprimer software, with mismatch and gap tolerance chosen by the user. If both tools predict PCR products using the same pair of primers with same length, the pair of primers is selected for the next step.
Primer classification. The program recovers all potentially useful pairs of primers for the differentiation of taxa. If a pair of primers has only one amplification product in only one target genome, it is defined as ‘specific’. If a pair of primers has one amplification product in at least one other genome, it is defined as ‘multiple’. If it amplifies the same size products in all genomes, it is named as ‘conserved’, but if the PCR products have different sizes in all genomes, then it is a ‘single’ primer. ‘Single’ pairs of primers are capable of distinguishing all taxa in a simple PCR reaction because each has different sizes of amplification products in each genome.
Compatibility check. Specific primers are clustered into compatible groups for multiplexing PCR using MultiPLX, which tests all primer pairs for interactions, including dimer formation and differences in their melting temperatures.
Gel visualization. After a set of primers is chosen, the relative electrophoretic migration distances are calculated based on the expected length of the amplification products. Then, a virtual electrophoresis gel is generated showing the expected amplification profile as a result of a standard PCR assay using a mixture of target genomes as the template. Another virtual gel is generated for the amplification profile of the multiplex PCR reaction, to check for interactions among primers that generate undesired alternative products.
List of specific (G1) and single (S1) pairs of primers for in vitro validation
Sequence Forward (5’ - > 3’)
(5’ - > 3’)
L. infantum, L. braziliensis and L. major
Results and discussion
Selected pairs of primers are classified in one of the four categories: specific, multiple, single and conserved. After the user chooses the categories, information regarding the primer characteristics is reported, such as primer sequence, melting temperature range (°C) and PCR product length (bp). Additionally, users can choose a set of primers or save the primer information in a text file or visualize a virtual gel electrophoresis. The text file shows the following amplicon properties: 1) name; 2) forward primer sequence; 3) reverse primer sequence; and 4) primer melting temperature and amplicon size in each target genome. A list of the compatible pairs of primers that are optimal for multiplex PCR is also available to download.
In the e-GEL function, the visual output is a simulated conventional PCR assay, where each lane is a reaction with one selected primer and a mixture of target genomes as a template. The e-MPX function generates another virtual gel with the amplification profile of the multiplex PCR assay, where each lane is a mixture of all selected primers and one target genome or all target genomes as the template.
Multiplex PCR is a variant of PCR, which simultaneously amplifies many loci of interest in one single reaction by using more than one pair of primers. Setting up a multiplex PCR with consistent quality is not trivial; therefore, TipMT generates a file with groups of specific primers that are compatible in a multiplex assay, based on primer–primer interactions and differences in the melting temperatures.
Leishmania is a genus of flagellate protozoan that cause a broad spectrum of diseases, ranging from self-limiting localized cutaneous lesions to visceral leishmaniasis. More than 20 species of Leishmania cause infection in humans . Despite the wide taxonomic complexity of this genus, the gold standard for diagnosing Leishmania infections, parasitological assays, only discriminates genus, not species. The reference method for species identification is multilocus enzyme electrophoresis (MLEE). However, this method has several limitations, including the relatively small number of characterized loci and the requirement of a parasite culture that potentially biases the results . Therefore, the development of new molecular diagnostic methods could allow for a rapid and accurate diagnosis of Leishmania species in infection. This is particularly important in the context of leishmaniasis, since the distinction between Leishmania species is important to provide the appropriate treatment and design the most effective control measures. This specific identification is especially relevant in areas where different species occur simultaneously causing human disease. Finally, robust molecular markers might contribute to the characterization of parasite-specific features, such as virulence or drug resistance .
We thus applied TipMT pipeline to the genomes of different species of Leishmania to generated sets of primers for genotyping using Multiplex PCR. First, we tested the pipeline using L. amazonensis, L. braziliensis and L. infantum genomes as templates on “Ortholog target” mode and among all pair of primers generated, we chose the G1 set for in vitro validation (Table 1, Fig. 2). Because in this search no pair of primer was classified as single, we re-ran the pipeline using L. braziliensis, L. major and L. infantum on “Ortholog target” mode and we obtained the pair of primers S1 which was also used in in vitro experiments (Table 1, Fig. 2c). We also generated pairs of primers using L. braziliensis and L. infantum genomes on SSR (R1, R2 and R3 groups) or Ortholog (O1, O2 and O3 groups) target modes and selected 12 pairs of primer for in vitro validation (Additional file 1: Table S1, Additional file 2: Figure S1).
The PCR gel electrophoresis profile revealed that the result of the real experiment was the same as that predicted by the virtual electrophoresis gel analysis (Fig. 2 and Additional file 2: Figure S1).
Comparison to other primer design applications
Main features of similar web-based primer design tools
Search for Target sequences
Specificity check method
multiple sequences (up to 500)
(limited and pre-defined list of sequences)
primers information and virtual gel
multiple sequences from multiple taxa
(SSR, orthologs, singletons)
alignment and thermodynamic
primers information and virtual gel
MPprimer is a web-based tool that designs specific multiplex PCR primers, uses thermodynamic theories to estimate the stability of the primers, and has functions for predicting a group of compatible multiplex primers and generating a virtual electrophoretic gel for each group. The main limitations of MPprimer are a limited and pre-defined list of genome databases (model organisms) provided by the tool that are used for checking primer specificity and the lack of a mechanism to search for target sequences, which should then be pre-defined by the user. BatchPrimer3 allows users to design several types of primer, including generic primers, hybridization oligos, primers for SSR regions, SNP genotyping primers and DNA sequencing primers. The main drawbacks of the tool are: limited number of sequences (up to 500) that can be used for primer design; no step for checking primer specificity; and the tool does not design primers for multiplex PCR. jPCR (FastPCR online) provides primers for most PCR applications and tests the specificity using a quick local alignment screen between the reference database (user’s database sequences) and input sequence. However, the application is platform dependent (Java Runtime Environment), requires computational resources of the user and more computing power for large databases. Also, jPCR only searches for SSR as target sequences.
TipMT offers a combination of features that are not present in any other available web applications. TipMT receives multiple sequences as input and identifies target regions automatically. Moreover, primer specificity is tested by both alignment- and thermodynamic-based properties. Furthermore, TipMT provides functions to generate a virtual electrophoresisgel for conventional or multiplex PCR assays. This output gives users a visual result before performing a real PCR reaction. Finally, the identification of taxa specific pairs of primers from multiple genomic sequences is a straightforward analysis only in TipMT, since it is the only available tool that can find pairs of primers for multiple taxa in a single run using genomic data without prior definition of the target region.
The emergence of large-scale DNA sequencing projects in recent years has produced large amounts of data, opening many opportunities for genomic analyses. Here, we focus our attention on identifying molecular markers and designing efficient primers for taxa differentiation. The ideal pair of primers should be capable of distinguishing the target taxon and should not cross-react with other closely related species. Toward this aim, TipMT receives genomic sequences as input and integrates the process of primer design, from the search for target sequences to the evaluation of primer specificity. As an output, the web-application generates a plain text file with general information on the pairs of primers, based on taxa-specificity. The output also includes an image showing the result of a simulated PCR assay with selected pairs of primers. Finally, experimental validation shows the effectiveness of the proposed tool in finding a taxon-specific pair of primers.
The pairs of primers generated by TipMT are suitable for use in conventional or multiplex PCR assays, as determined by the parameter settings during the primer design step. Furthermore, primer design principles for conventional PCR and real-time quantitative PCR are quite similar. Thus, our tool could be used for designing primers for both methodologies by adjusting some parameters, such as PCR product size.
Future versions of TipMT will consider multi-copy genes as targets to improve PCR sensitivity and will also receive raw sequencing reads as input. Additional improvements may also be performed, for example, automatic ranking of ideal primers for multiplex PCR.
The application is platform independent, freely available and has a simple and user-friendly interface that allows for designing primers in a high-throughput manner, even for novice users. Furthermore, TipMT web page has a “Manual” section with a tutorial and running examples. TipMT can be applied to a broad spectrum of research topics including both molecular diagnostic and evolutionary studies.
Availability and requirements
Project name: TipMT, Tool for Identification of Primers for Multiple Taxa
Project home page: http://184.108.40.206/tipMT/
Operating system(s): Platform independent
Other requirements: Web browser (supported browsers: Firefox, Chrome)
Any restrictions to use by non-academics: no license needed
Basic local alignment search tool
Hypertext markup language
PHP hypertext preprocessor
We thank Michele Silva de Matos for her technical assistance.
This study was funded by Fundação de Amparo a Pesquisa do Estado de Minas Gerais (FAPEMIG), Instituto Nacional de Ciência e Tecnologia de Vacinas (INCTV), Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) and Pró-Reitoria de Pesquisa (PRPq), Universidade Federal de Minas Gerias (UFMG). DCB, CG, RTF are CNPq research fellows. GFRL, HOV and EVA received scholarships from CAPES and MSC received a scholarship from CNPq. The funding body has no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Availability of data and material
Data sharing is not applicable to this article as no genomic data were generated during the current study.
GFRL carried out the code development, implementation, and drafted the manuscript. MCS, HOV and EVA validated SSR and Orthologs primers. RSL implemented the TipMT website construction. CMFG, TSR and RTF provided scientific advice and the resources to develop the software. DCB conceived of the study, participated in its design and coordination, and manuscript writing. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Lucchi NW, Oberstaller J, Kissinger JC, Udhayakumar V. Malaria diagnostics and surveillance in the post-genomic era. Public Health Genomics. 2013;16:37–43.View ArticlePubMedPubMed CentralGoogle Scholar
- Wong SS, Fung KS, Chau S, Poon RW, Wong SC, Yuen K-Y. Molecular diagnosis in clinical parasitology: When and why? Exp Biol Med (Maywood). 2014;239:1443–60.
- Li J, Zhang Y, Liu S, Hong L, Sullivan M, McCutchan TF, Carlton JM, Su XZ. Hundreds of microsatellites for genotyping Plasmodium yoelii parasites. Mol Biochem Parasitol. 2009;166:153–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Ellegren H. Microsatellites: simple sequences with complex evolution. Nat Rev Genet. 2004;5:435–45.View ArticlePubMedGoogle Scholar
- Guichoux E, Lagache L, Wagner S, Chaumeil P, Léger P, Lepais O, Lepoittevin C, Malausa T, Revardel E, Salin F, Petit RJ. Current trends in microsatellite genotyping. Mol Ecol Resour. 2011;11:591–611.View ArticlePubMedGoogle Scholar
- Sharma PC, Grover A, Kahl G. Mining microsatellites in eukaryotic genomes. Trends Biotechnol. 2007;25:490–8.View ArticlePubMedGoogle Scholar
- Duran C, Appleby N, Edwards D, Batley J. Molecular genetic markers: discovery, applications, data storage and visualisation. 2009, 61:16–27.
- Powell S, Szklarczyk D, Trachana K, Roth A, Kuhn M, Muller J, Arnold R, Rattei T, Letunic I, Doerks T, Jensen LJ, von Mering C, Bork P. eggNOG v3.0: orthologous groups covering 1133 organisms at 41 different taxonomic ranges. Nucleic Acids Res. 2012;40(Database issue):D284–9.View ArticlePubMedGoogle Scholar
- Li L, Stoeckert CJ, Roos DS. OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003;13:2178–89.View ArticlePubMedPubMed CentralGoogle Scholar
- Magi A, Benelli M, Gozzini A, Girolami F, Torricelli F, Brandi ML. Bioinformatics for next generation sequencing data. Genes (Basel). 2010;1:294–307.Google Scholar
- Abd-elsalam K a. Bioinformatic tools and guideline for PCR primer design. African J Biotechnol. 2003;2:91–5.View ArticleGoogle Scholar
- Cao Y, Wang L, Xu K, Kou C, Zhang Y, Wei G, He J, Wang Y, Zhao L. Information theory-based algorithm for in silico prediction of PCR products with whole genomic sequences as templates. BMC Bioinformatics. 2005;6:190.View ArticlePubMedPubMed CentralGoogle Scholar
- Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, Rozen SG. Primer3--new capabilities and interfaces. Nucleic Acids Res. 2012;40:e115.View ArticlePubMedPubMed CentralGoogle Scholar
- Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26:841–2.View ArticlePubMedPubMed CentralGoogle Scholar
- Stajich JE, Block D, Boulez K, Brenner SE, Chervitz S a, Dagdigian C, Fuellen G, Gilbert JGR, Korf I, Lapp H, Lehväslaiho H, Matsalla C, Mungall CJ, Osborne BI, Pocock MR, Schattner P, Senger M, Stein LD, Stupka E, Wilkinson MD, Birney E. The Bioperl toolkit: Perl modules for the life sciences. Genome Res. 2002;12:1611–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.View ArticlePubMedGoogle Scholar
- Rice P, Longden I, Bleasby A. EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet. 2000;16:276–7.View ArticlePubMedGoogle Scholar
- Rotmistrovsky K, Jang W, Schuler GD. A web server for performing electronic PCR. Nucleic Acids Res. 2004;32(Web Server issue):W108–12.View ArticlePubMedPubMed CentralGoogle Scholar
- Qu W, Zhou Y, Zhang Y, Lu Y, Wang X, Zhao D, Yang Y, Zhang C. MFEprimer-2.0: a fast thermodynamics-based program for checking PCR primer specificity. Nucleic Acids Res. 2012;40(Web Server issue):gks552–.Google Scholar
- Kaplinski L, Andreson R, Puurand T, Remm M. MultiPLX: automatic grouping and evaluation of PCR primers. Bioinformatics. 2005;21:1701–2.View ArticlePubMedGoogle Scholar
- Lopes R da S, Moraes JWL, Rodrigues TDS, Bartholomeu DC: ProGeRF: Proteome and Genome Repeat Finder Utilizing a Fast Parallel Hash Function. 2015;2015:1–9.
- Andreson R, Möls T, Remm M. Predicting failure rate of PCR in large genomes. Nucleic Acids Res. 2008;36:e66.View ArticlePubMedPubMed CentralGoogle Scholar
- Remm M, Ants K, Metspalu A. Primer Design for Large-Scale Multiplex PCR and Arrayed Primer Extension (APEX). In: PCR technology: Current innovations. 2nd ed. 2004. p. 131–40.Google Scholar
- Alvar J, Vélez ID, Bern C, Herrero M, Desjeux P, Cano J, Jannin J, den Boer M, WHO Leishmaniasis Control Team. Leishmaniasis worldwide and global estimates of its incidence. PLoS One. 2012;7:e35671.View ArticlePubMedPubMed CentralGoogle Scholar
- Hernández C, Ramírez JD. Molecular diagnosis of vector-borne parasitic diseases. 2013, 2:1–10.
- Reithinger R, Dujardin J-C. Molecular diagnosis of leishmaniasis: current status and future applications. J Clin Microbiol. 2007;45:21–5.View ArticlePubMedGoogle Scholar
- You FM, Huo N, Gu YQ, Luo M-C, Ma Y, Hane D, Lazo GR, Dvorak J, Anderson OD. BatchPrimer3: a high throughput web application for PCR and sequencing primer design. BMC Bioinformatics. 2008;9:253.View ArticlePubMedPubMed CentralGoogle Scholar
- Kalendar R, Lee D, Schulman A. FastPCR software for PCR, in silico PCR, and oligonucleotide assembly and analysis. DNA Cloning Assem Methods. 2014;1116:271–302.
- Shen Z, Qu W, Wang W, Lu Y, Wu Y, Li Z, Hang X, Wang X, Zhao D, Zhang C. MPprimer: a program for reliable multiplex PCR primer design. BMC Bioinformatics. 2010;11:143.View ArticlePubMedPubMed CentralGoogle Scholar