IMP PCR primers detect single nucleotide polymorphisms for Anopheles gambiae species identification, Mopti and Savanna rDNA types, and resistance to dieldrin in Anopheles arabiensis
Malaria Journal volume 5, Article number: 125 (2006)
Polymerase chain reactions to distinguish single-nucleotide polymorphisms are commonly used for mosquito identification and identifying insecticide resistance alleles. However, the existing methods used for primer design often result in analyses that are not robust or require additional steps.
Utilizing oligonucleotides that are unique in having an intentional mismatch to both templates three bases from the SNP at the 3-prime end, three new PCR assays that distinguish SNP targets using standard gel electrophoresis of undigested DNA fragments were developed and tested. These were applied to: (1) an alternative ribosomal DNA PCR assay to distinguish five members of the Anopheles gambiae complex; (2) detection of the Mopti and Savanna rDNA types; and (3) an assay to distinguish resistance to dieldrin (Rdl) alleles in Anopheles arabiensis.
Reproducible specific amplification of the target alleles was observed in all three assays. The results were consistent with existing analyses but proved simpler and the results more distinct in our hands.
The simplicity and effectiveness of the method should be utilized in these and other PCR analyses to increase their specificity and simplicity. These results have the potential to be extended not only to mosquito analyses but also to parasite and human polymorphisms.
Single nucleotide polymorphisms (SNPs) have been utilized as markers for species differentiation  and are implicated in insecticide resistance in many anopheline mosquitoes . While various polymerase chain reaction (PCR) approaches to discriminate SNPs exist, specific primers for the wild-type and/or mutant alleles are often designed to perfectly match the sequence of one allele whereas the second varies only by matching the alternative SNP sequence at the terminal 3' end. Therefore, it is critical that the biochemical specificity of the two primers is sufficient to determine the success of PCR amplification from a particular template. This PCR-specificity depends on two characteristics: the ability of any mismatch to prevent extension of the primer by the polymerase – many polymerases can synthesize template over a single terminal mismatch  – and the difference between the annealing temperature (Tm) of the perfect and mismatched primer . For both of these reasons, a single 3' mismatch has limited discrimination. In this report, intentional mismatches were used to increase the discrimination of PCR based primarily on the polymerase synthesis characteristic.
Three PCR analyses were created by the method described here:
The Anopheles gambiae complex is comprised of six genetically and behaviourally distinct species that are morphologically almost identical . Members of this complex have been identified as major vectors of human malaria parasites, and as many as four species may be sympatric. Several methods for identifying these species have been developed, such as gas chromatography of cuticular hydrocarbons , polytene chromosome arrangements (reviewed in ) and allozyme analyses . However, none of these has obtained as widespread use as PCR amplification of the ribosomal DNA (rDNA). This has become the standard method for species discrimination since rDNA polymorphism detected by Southern hybridization  spurred development of a PCR assay . While the Scott et al. PCR method is the most popular, non-specific bands and high rates of PCR amplification failure were observed in our hands.
In addition to the cryptic species of the An. gambiae complex described above, karyotype analysis of the female ovarian chromosomes has revealed further subdivision into chromosomal forms called Forrest, Mopti (M) and Savanna (S) . The observed frequencies of the inversions that define these types are correlated with different ecological zones suggesting that they may have an adaptive character. Because karyotype analysis is tedious and can only be performed on semi-gravid females, several alternative PCR-based methods have been developed. PCR amplification of the rDNA intergenic spacer IGS [10–12] and internal transcribed spacers ITS1 and ITS2  of An. gambiae s.s. have provided a molecular surrogate for karyotypes of Mopti and Savanna types in some but not all areas . The novel method reported here was compared with the diagnostic PCR-RFLP method of Fanello  which amplifies a portion of the IGS followed by digestion of the products with Hha I to produce diagnostic fragments of either rDNA type.
Resistance to dieldrin (Rdl) in An. gambiae s. l. is due to a single nucleotide mutation within the M2 transmembrane subunit of γ-aminobutyric acid (GABA) receptor . In Anopheles arabiensis, resistance is conferred by an alanine to serine substitution within the Rdl locus. Resistance to dieldrin per se (a cyclodiene) is of little immediate concern since it has been banned from use in public health and agriculture. However, recently it has been shown that cross-resistance to the phenylpyrazole fipronil in Anopheles stephensi and An. gambiae was also due to the presence of an Rdl resistance allele [15, 16].
In this report, intentional mismatch primers – IMPs- a method previously utilized by Papp et al.  – were used to introduce mismatches at the third nucleotide from the 3' end. This simple modification of intentionally adding a mismatch to the SNP-detection primers counter-intuitively increases the specificity of the primers, which is demonstrated here for three anopheline SNP detection analyses.
Mosquitoes were obtained from the Malaria Research and Reference Reagent Resource Center (MR4) holdings at the CDC in Atlanta, GA USA (Table 1). Larvae were reared at 27°C using a standard method  except that larva were fed Aquaricare™ Koi Floating Blend from the L2 stage to pupation.
Samples were prepared for PCR by the method of Rafferty et al. . PCR products were observed by separation on 0.5× TBE agarose gels run in 0.5× TBE buffer at 12 v/cm and fragment sizes were estimated using a 1 kb ladder marker (Invitrogen®). Thermal cycling for all analyses was performed in a Bio-Rad iCycler®. Taq DNA polymerase and the manufacturer's (Promega®) recommended buffer at 1× concentration was used for all reactions. PCR reactions consisted of 1 U of Taq polymerase, 0.3 mM MgCl2, all primers at 1 μM except QD-3T at 2 μM, 0.08 mM dNTPs, and buffer in 25 μl total volume. All primers used in these studies except those of Scott et al. are listed in Table 2.
An. gambiae species identification PCR cycling consisted of melting at 95°C for 5 min followed by 30 cycles of 95°C for 30 seconds, 58°C for 30 seconds, and 72°C for 30 seconds, followed by one cycle of 72°C for 5 min. Mopti and Savanna rDNA were analysed using PCR-amplified rDNA and Hha I digests by the method of Fanello . Hha I digests were performed overnight at 37°C. The locations of novel primers were selected at SNP sites using published sequence data [Genbank: AF470116 and AF470112]. PCR thermal cycling for the IMP primers consisted of melting at 95°C for 5 min followed by 30 cycles of 95°C for 30 seconds, 54°C for 30 seconds, and 72°C for 30 seconds, followed by one cycle of 72°C for 5 min.
An. arabiensis Rdl-specific primers were designed to be used with the previously reported An. gambiae internal positive control primers RDLF and RDLR . Thermal cycling consisted of melting at 95°C for 5 min followed by 30 amplification cycles of 95°C for 45 seconds, 53°C for 45 seconds, and 72°C for 45 seconds, followed by one cycle of 72°C for 10 min. Heterozygous F1 individuals were produced for analysis by crossing KGB (homozygous susceptible) females to SENN (homozygous resistant) males en masse.
Results and discussion
The primer strategy for An. gambiae species identification was similar to that of Scott et al. : one universal forward primer was used for all species in a cocktail with four species-specific reverse primers. However, neither the species-specific nor universal annealing sites were the same (Figure 1a). The species-specific DNA fragments that we obtain using the IMP primer set for An. gambiae s.l. species identification (Table 2) compare favorably to the products we obtain with the Scott et al. assay (Figure 2). While in the latter assay, non-specific artifact bands were observed, the IMP-protocol products are clear and consistent. Prior to developing these assays, primers were designed for the species identification except incorporating the deliberate mismatch at the penultimate base. These did not provide the desired level of specificity with the mismatch at this site. Additional primers were designed with a mismatch and the SNP in the middle in an attempt to thermodynamically alter the melting temperature sufficiently for discrimination. Neither of these approaches was successful.
In this study, only laboratory colonies of mosquitoes were tested for species ID, and validation on a larger scale as has been performed previously for Anopheles quadrimaculatus s.l.  will be necessary to confirm their utility. In order to ensure that the IMP primer sets would function correctly in strains other than those we analysed and for which sequence data was available, reported sequences were examined [1, 12, 21–23] and no polymorphisms were found that would interfere with the success of the protocol.
An alternative to the Fanello protocol for analysis of the IGS was developed that is in principle like that of Favia et al.  in that PCR, but not restriction digestion, is required. In the Fanello analysis, the expected resulting fragments of the M (367 + 23 bp) and S (257 + 110 + 23 bp) types include one only slightly smaller than the undigested 390 bp PCR product, therefore, high percentage agarose gels are used to distinguish products. The IMP primer design includes single base mismatches between the M and S rDNA in the target regions. They were designed to amplify distinct fragments of 426 or 335 bp in the M and S forms respectively (Figure. 3). The results with these primers are consistent with those of the Fanello protocol and additional evidence of their consistency was obtained in a larger unrelated study.
The Rdl dieldrin resistance PCR assay contains two perfect-match primers (RDLF & RDLR, ) that amplify the region containing the SNP and one forward and one reverse IMP primer that specifically bind the SNP found in An. arabiensis. The SNP site in An. arabiensis is located asymmetrically from the middle of the region between the RDLF and RDLR primers so we were able to design a forward and a reverse primer at that site (Figure 1b), the forward primer being specific for the resistant form and the reverse specific for the susceptible form. The four primers are used in a cocktail together and are expected to yield a control band and a resistant or susceptible band in homozygotes depending on the templates. The result of using the four-primer cocktail with a heterozygote yielded three bands as expected (Figure 4).
The novel method, upon which these assays is based, is that all SNP-specific primers were designed to contain an intentional mismatch for both target alleles at the third nucleotide from the 3' end and with the terminal base falling on the site of the SNP . These primers were tested over a wide range of annealing temperatures using a gradient cycler and were found to amplify their targets consistently and clearly over annealing temperatures of 50–65°C. This bodes favorably for reproducibility of the results obtained here on different machines. Successful amplification over a wide range of temperatures also suggests that the primer design is successful not due to a thermodynamic change leading to a lower Tm but rather due to the inability of a standard Taq polymerase to extend over a series with alternating mismatches at the 3' end.
Successful amplification using IMPs was achieved with all mosquito stages tested including larvae, pupae, and adult mosquitoes without extensive DNA purification. Therefore, these assays are useful for studies in the field where DNA preparation is not possible or time-prohibitive. The ease of design, technical simplicity, high specificity, and sensitivity recommend these as improved specific assays and the general method as a means to produce assays for both multicopy and single-copy SNPs for laboratory and field analyses. There is no apparent reason that similar primer designs would not be useful for analysis of human and Plasmodium SNPs and any application in which rapid sensitive SNP discrimination is necessary.
Scott JA, Brogdon WG, Collins FH: Identification of single specimens of the Anopheles gambiae complex by the polymerase chain reaction. Am J Trop Med Hyg. 1993, 49: 520-529.
Hemingway J, Hawkes NJ, McCarroll L, Ranson H: The molecular basis of insecticide resistance in mosquitoes. Insect Biochem Mol Biol. 2004, 34: 653-665. 10.1016/j.ibmb.2004.03.018.
Zhang J, Li K, Liao D, Pardinas JR, Chen L, Zhang X: Different applications of polymerases with and without proofreading activity in single-nucleotide polymorphism analysis. Lab Invest. 2003, 83: 1147-1154. 10.1097/01.LAB.0000081589.91390.DF.
Guo Z, Liu Q, Smith LM: Enhanced discrimination of single nucleotide polymorphisms by artificial mismatch hybridization. Nat Biotechnol. 1997, 15: 331-335. 10.1038/nbt0497-331.
Davidson G, Paterson HE, Coluzzi M, Mason GF, Micks DW: The Anopheles gambiae Complex. Genetics of Insect Vectors of Disease. Edited by: Wright JW, Pal R. 1967, Elsevier, 211-250.
Carlson DA, Service MW: Identification of mosquitoes of Anopheles gambiae species complex A and B by analysis of cuticular components. Science. 1980, 207: 1089-1091. 10.1126/science.7355276.
Coluzzi M, Sabatini A, Petrarca V, Di Deco MA: Chromosomal differentiation and adaptation to human environments in the Anopheles gambiae complex. Trans R Soc Trop Med Hyg. 1979, 73: 483-497. 10.1016/0035-9203(79)90036-1.
Mahon RJ, Green CA, Hunt RH: Diagnostic allozymes for routine identification of adults of the Anopheles gambiae complex (Diptera, Culicidae). Bull Entomol Res. 1976, 66: 25-31.
Collins FH, Mendez MA, Rasmussen MO, Mehaffey PC, Besansky NJ, Finnerty V: A ribosomal RNA gene probe differentiates member species of the Anopheles gambiae complex. Am J Trop Med Hyg. 1987, 37: 37-41.
McLain DK, Collins FH, Brandling-Bennett AD, Were JBO: Microgeographic variation in rDNA inergenic spacers of Anopheles gambiae in western Kenya. Heredity. 1989, 62: 257-264.
Fanello C, Santolamazza F, della Torre A: Simultaneous identification of species and molecular forms of Anopheles gambiae complex by PCR-RFLP. Med Vet Entomol. 2002, 16: 461-464. 10.1046/j.1365-2915.2002.00393.x.
Favia G, Lanfrancotti A, Spanos L, Siden-Kiamos I, Louis C: Molecular characterization of ribosomal DNA polymorphisms discriminating among chromosomal forms of Anopheles gambiae s.s. Insect Mol Biol. 2001, 10: 19-23. 10.1046/j.1365-2583.2001.00236.x.
Gentile G, Slotman M, Ketmaier V, Powell JR, Caccone A: Attempts to molecularly distinguish cryptic taxa in Anopheles gambiae s.s. Insect Mol Biol. 2001, 10: 25-32. 10.1046/j.1365-2583.2001.00237.x.
Du W, Awolola TS, Howell P, Koekemoer LL, Brooke BD, Benedict MQ, Coetzee M, Zheng L: Independent mutations in the Rdl locus confer dieldrin resistance to Anopheles gambiae and An. arabiensis. Insect Mol Biol. 2005, 14: 179-183. 10.1111/j.1365-2583.2005.00544.x.
Brooke BD, Hunt RH, Coetzee M: Resistance to dieldrin + fipronil assorts with chromosome inversion 2La in the malaria vector Anopheles gambiae. Med Vet Entomol. 2000, 14: 190-194. 10.1046/j.1365-2915.2000.00222.x.
Kolaczinski J, Curtis C: Laboratory evaluation of fipronil, a phenylpyrazole insecticide, against adult Anopheles (Diptera: Culicidae) and investigation of its possible cross-resistance with dieldrin in Anopheles stephensi. Pest Manag Sci. 2001, 57: 41-45. 10.1002/1526-4998(200101)57:1<41::AID-PS260>3.0.CO;2-N.
Papp A, Pinsonneault J, Cooke G, Sadee W: Single nucleotide polymorphism genotyping using allele-specific PCR and fluorescence melting curves. Biotechniques. 2003, 34: 1068-1072.
Benedict MQ: Care and maintenance of anopheline mosquito colonies. The Molecular Biology of Insect Disease Vectors. Edited by: Crampton JM, Beard CB, Louis C. 1997, New York: Chapman & Hall, 2-12.
Rafferty C, Campbell S, Wirtz R, Benedict M: Polymerase chain reaction-based identification and genotyping of Anopheles mosquitoes with a 96-pin bacterial replicator. Am J Trop Med Hyg. 2002, 66: 234-237.
Rutledge CR, Cornel AJ, Meek CL, Collins FH: Validation of a ribosomal DNA-polymerase chain reaction species diagnostic assay for the common malaria mosquito (Diptera:Culicidae) sibling species complex. J Med Entomol. 1996, 33: 952-954.
Gentile G, Della Torre A, Maegga B, Powell JR, Caccone A: Genetic differentiation in the African malaria vector, Anopheles gambiae s.s., and the problem of taxonomic status. Genetics. 2002, 161: 1561-1578.
Besansky NJ, Powell JR, Caccone A, Hamm DM, Scott JA, Collins FH: Molecular phylogeny of the Anopheles gambiae complex suggests genetic introgression between principal malaria vectors. Proc Natl Acad Sci USA. 1994, 91: 6885-6888. 10.1073/pnas.91.15.6885.
Cornel AJ, Collins FH: PCR of the ribosomal DNA intergenic spacer regions as a method for identifying mosquitoes in the Anopheles gambiae complex. Methods in Molecular Biology: Species Diagnostic Protocols: PCR and other nucleic acid methods. Edited by: Clapp JP. 1996, Totowa NJ: Humana Press, 50: 321-332.
We appreciate the American Type Culture Collection and the NIAID-funded contract N01-AI-85355 for their financial and moral support for the operation of the MR4 and their cooperation in developing and implementing the subcontract to the CDC Foundation that makes this research possible. Special thanks to all the donors of mosquito stocks without whose collegial generosity these experiments would not have been possible: Maureen Coetzee (KGB and SKUQUA), Vincent Corbel (KISUMU1), Mustafa Dukeen (SENN), Greg Lanzaro (MOPTI), and Rajendra Maharaj (OPHANSI). Thanks to Lynn Huynh for the careful editing and discussion of this manuscript.
EW carried out the primer design research, primer design and performed the first draft of the manuscript. PH participated in background research. Both EW and PH conducted molecular analysis, mosquito culture, and reviewed numerous revisions of the manuscript. MB participated in the design and coordination of the studies and helped to draft and performed the final review of the manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Wilkins, E.E., Howell, P.I. & Benedict, M.Q. IMP PCR primers detect single nucleotide polymorphisms for Anopheles gambiae species identification, Mopti and Savanna rDNA types, and resistance to dieldrin in Anopheles arabiensis. Malar J 5, 125 (2006). https://doi.org/10.1186/1475-2875-5-125
- Polymerase Chain Reaction
- Polymerase Chain Reaction Amplification
- Polymerase Chain Reaction Assay
- Anopheles Arabiensis