PCR-based methods for the detection of L1014 kdr mutation in Anopheles culicifacies sensu lato
© Singh et al. 2009
Received: 09 March 2009
Accepted: 14 July 2009
Published: 14 July 2009
Skip to main content
© Singh et al. 2009
Received: 09 March 2009
Accepted: 14 July 2009
Published: 14 July 2009
Anopheles culicifacies s.l., a major malaria vector in India, has developed widespread resistance to DDT and is becoming resistant to pyrethroids–the only insecticide class recommended for the impregnation of bed nets. Knock-down resistance due to a point mutation in the voltage gated sodium channel at L1014 residue (kdr) is a common mechanism of resistance to DDT and pyrethroids. The selection of this resistance may pose a serious threat to the success of the pyrethroid-impregnated bed net programme. This study reports the presence of kdr mutation (L1014F) in a field population of An. culicifacies s.l. and three new PCR-based methods for kdr genotyping.
The IIS4-IIS5 linker to IIS6 segments of the para type voltage gated sodium channel gene of DDT and pyrethroid resistant An. culicifacies s.l. population from the Surat district of India was sequenced. This revealed the presence of an A-to-T substitution at position 1014 leading to a leucine-phenylalanine mutation (L1014F) in a few individuals. Three molecular methods viz. Allele Specific PCR (AS-PCR), an Amplification Refractory Mutation System (ARMS) and Primer Introduced Restriction Analysis-PCR (PIRA-PCR) were developed and tested for kdr genotyping. The specificity of the three assays was validated following DNA sequencing of the samples genotyped.
The genotyping of this An. culicifacies s.l. population by the three PCR based assays provided consistent result and were in agreement with DNA sequencing result. A low frequency of the kdr allele mostly in heterozygous condition was observed in the resistant population. Frequencies of the different genotypes were in Hardy-Weinberg equilibrium.
The Leu-Phe mutation, which generates the kdr phenotype in many insects, was detected in a pyrethroid and DDT resistant An. culicifacies s.l. population. Three PCR-based methods were developed for kdr genotyping. All the three assays were specific. The ARMS method was refractory to non-specific amplification in non-stringent amplification conditions. The PIRA-PCR assay is able to detect both the codons for the phenylalanine mutation at kdr locus, i.e., TTT and TTC, in a single assay, although the latter codon was not found in the population genotyped.
Anopheles culicifacies s.l. is the main malaria vector in the Indian subcontinent, affecting mainly rural areas, and contributes 60–65% of malaria cases in India . As this is an endophilic vector, indoor residual spraying (IRS) of DDT is the main strategy used for its control. This species is resistant to DDT in most parts of India. Although DDT is banned in many countries, the recent endorsement by World Health Organization for the use of DDT for IRS for malaria vector control  has renewed interest in this insecticide. Pyrethroids are the most commonly used insecticides for IRS and the only insecticide class recommended for impregnation of bed nets due to their relatively low mammalian toxicity and rapid knock down effect on insects. In India, the use of pyrethroids was initiated in 1990s to control malaria epidemics in areas where An. culicifacies s.l was resistant to DDT and malathion. Pyrethroid resistance in An. culicifacies s.l. was detected in Surat district of Gujarat state, western India, soon after its introduction in vector control programme .
DDT and pyrethroids are neurotoxins that act on the voltage-gated sodium channels by modifying their gating kinetics, resulting in the prolonged opening of individual channels leading to paralysis and death of the insect. One of the mechanisms of pyrethroid resistance in insects is referred to as knock-down resistance (kdr) caused by reduced target site sensitivity. The phenotype is commonly conferred by a single point mutation (L1014F/S/H) in the IIS6 segment of voltage gated sodium channel [4, 5]. Other mutations in different regions of the gene also confer knock-down resistance in some insects [4, 6], but among anophelines this is the only locus where point mutations have been reported to date conferring resistance. Only two points mutations have been reported in anophelines at this locus–L1014F in Anopheles gambiae  (West African kdr), Anopheles arabiensis  and Anopheles stephensi , and L1014S in An. gambiae  (East African kdr).
The Leu-Phe mutation at the kdr locus in An. culicifacies s.l. was reported by Hoti et al  using an allele-specific PCR assay, whose external primers Agd1 and Agd2 were based on An. gambiae sequences. However, authors failed to amplify the DNA region of interest using the primers Agd1 and Agd2 (see Additional file 1) due to mismatching with template DNA. As the evident cross resistance between DDT and pyrethroids in this resistant strain strongly suggested a kdr type phenotype, an attempt was made to confirm the presence of kdr -based insecticide resistance in An. culicifacies s.l. by DNA sequencing and developing alternative high throughput methods for kdr genotyping. The part of the voltage gated sodium channel spanning IIS4-IIS5 linker to IIS6 segments of An. culicifacies species A, B and C was first sequenced which contains at least five residues where mutations have been reported in other insects, namely, Met918 in the IIS4-IIS5 linker, Leu925, Thr929 and Leu932 in IIS5 and Leu1014 in IIS6 . Based on the sequences obtained, three PCR based assays for kdr genotyping were developed, i.e., Allele Specific-PCR (AS-PCR), Amplification Refractory Mutation System (ARMS) and Primer Introduced Restriction Analysis-PCR (PIRA-PCR). These assays were tested on a DDT/pyrethroid resistant An. culicifacies population from Surat district and validated through DNA sequencing. The ARMS and PIRA-PCR were developed as alternatives of AS-PCR, which is often alleged to be non-specific [12–14]. Additionally, the PIRA-PCR was developed as a method by which both the codons of Phe (TTT and TTC) at kdr locus can be detected by a single assay.
Anopheles culicifacies s.l were collected from the villages of Surat district of Gujarat state (21–22°N latitude and 73–74°E longitude) located within three Primary Health Centres (PHC), Amladam, Bareda and Dadheda, where a significant level of resistance to DDT and pyrethroids (7–13% and 60–69% mortalities recorded in adults against 4% DDT and 0.05% Deltamethrin, respectively, in WHO's standard insecticide susceptibility tests) have been reported . This An. culicifacies s.l. population is comprised of species B and C . Indoor resting mosquitoes were collected using hand aspirators between 0600 and 0800 AM. The mosquitoes were preserved in isopropanol and DNA from individual mosquitoes was isolated.
The total RNA of An. culicifacies species A, B and C was isolated using TRIzol Reagent (Invitrogen Corporation) according to the manufacturer's instructions. cDNA synthesis and RT-PCR were carried out using AccessQuick RT-PCR System (Promega Corporation). The first strand cDNA was synthesized by avian myoblastosis virus reverse transcriptase (AMV-RT) at 45°C for 45 min, followed by denaturation of RNA/DNA hybrid and inactivation of AMV-RT at 95°C for 2 min. Second strand cDNA synthesis and subsequent DNA amplification were done by Tfl DNA polymerase at the following conditions: 35 cycles of denaturation at 95°C for 30 sec, annealing at 45°C for 30 sec, and extension at 72°C for 1 min followed by a final extension of 5 min at 72°C. The primers used for RT-PCR were Kdr1F and Kdr1R. The product was re-amplified by semi-nested PCR using primers Kdr1F and KdrR, and thereafter sequencing was performed using BigDye Terminator v3.1.
The DNA sequences obtained from cDNA and genomic DNA analyses were deposited in GenBank (accession numbers: GQ268331–GQ268333 and GQ279245–GQ279247).
Three molecular methods for kdr genotyping were developed based on sequences obtained in this study, viz., Allele Specific PCR (AS-PCR), Amplification Refractory Mutation System (ARMS) and Primer Introduced Restriction Analysis PCR (PIRA-PCR) (Figure 1). In all the PCR-based assays, unless otherwise stated, uniform concentrations of reagents, except primers, were taken, which were: 1× Buffer, 1.5 mM of MgCl2, 200 μM of each dNTP and 0.5 unit AmpliTaq Gold® per 15 μl of reaction mixture. All the PCR reactions were carried out in a GeneAmp® PCR System 9700 (Applied Biosystem) thermal cycler. The three assays were tested on field collected An. culicifacies s.l for kdr genotyping and results were validated following DNA sequencing.
To check the specificity of ARMS, the assay was also performed with increased concentration of MgCl2 (up to 3.0 mM) and AmpliTaq® Taq DNA Polymerase (up to 2.5 units per 15 μl reaction mixtures) at an annealing temperature of 45°C.
PIRA-PCR is based on incorporation of a deliberate mismatch in one primer so that a new restriction enzyme recognition site is created in the PCR amplicon in the presence of a specific allele in the target DNA. A PIRA primer KdrP (forward) with sequence 5'-TCC TGG CTA CAG TAG TGA TAG GAA AaT T-3' was designed, where a T-to-A mismatch was incorporated at the 3rd base position from 3' end terminus (shown in small letter in primer sequence). The terminal AAaTT-3' sequence of the PIRA primer (KdrP) creates a recognition site for the restriction enzyme Apo I or Xap I, i.e., AAATTY, in the amplified PCR product with the primer KdrR (reverse) in the presence of the Phe codon (TTT or TTC) in the template DNA at the kdr locus. Hence the PIRA-PCR is capable of detecting both possible codons of phenylalanine (TTT or TTC). With the Leu (TTA) allele, no Apo I restriction site is created.
The predicted size of PCR product amplified with primers KdrP and KdrR is 134 base pairs, which when digested with ApoI or XapI, will be cleaved into 106 and 24 bp fragments (with four bases 5'-overhang) in the presence of Phe allele (TTT/TTC). No cleavage is expected to take place in presence of Leu (TTA) allele. The 134 and 106 bp fragments can be easily resolved on 3% agarose gel whereas 24 bp fragment is normally not expected to be visualized or difficult to discriminate from primer-dimer. Therefore, 134 bp band was scored as S and 106 bp as R.
All the mosquito samples genotyped as RS (n = 18) and RR (n = 1) were sequenced alongside 23 SS in order to validate assays and to check for other mutations in addition to the Leu-Phe change. The primers used for PCR amplification and sequencing were KdrF and KdrR. The sequence data were aligned and analysed using Mega 4 alignment utilities .
The Exact Test of Hardy-Weinberg equilibrium using a Markov chain on three genotypes, i.e., SS, RS and RR was done using software Arlequin ver 3.11 
DNA sequencing revealed the presence of a Leu (TTA) to Phe (TTT) mutation at residue L1014 in a total of 19 samples, 18 of which were heterozygote and one homozygous resistant individual. No other mutations were recorded, including L995H mutation reported in Iranian An. culicifacies s.l. .
Testing of AS-PCR and ARMS assays on An. culicifacies s.l. samples from Surat revealed presence of discrete bands of expected sizes, as visualized on agarose gel. With PIRA-PCR also all the expected bands were seen, except 24 bp fragment resulting from the cleavage of PCR product (134 bp). Therefore, the presence of 106 bp cleaved product was scored as R. A total of 186 samples were genotyped using AS-PCR, of which 167 were assigned as homozygous susceptible (SS), one as homozygous resistant (RR) and 18 as heterozygote (SR) based on the criteria mentioned in section 'Methods'. All of the RS and RR, and 80 of the SS samples, as genotyped by AS-PCR, were also assayed with ARMS and PIRA-PCR assays, with similar results to the AS-PCR.
The result of DNA sequencing of 42 samples genotyped by all the three assays (RR = 1, SR = 18, SS = 23) were in agreement with PCR based assays. Thus all three assays are specific and can be used for genotyping.
The AS-PCR, based on a single mismatch at the 3' end of the allele specific primer, was specific and no mispriming was observed. The ARMS was specific even with increased concentrations of MgCl2 (3.0 mM) and taq polymerase (3 Units of AmpliTaq Gold per 25 μl reaction).
The observed and expected heterozygote frequencies were 0.09677 and 0.10202 respectively, which did not differ significantly (p = 0.41220) from Hardy-Weinberg equilibrium.
Malaria control in India relies on control of the major vector An. culicifacies. Knock-down resistance in An. culicifacies may be a serious threat to the continued use of pyrethroids, which are currently the insecticide class of choice for control of multiple insecticide resistant vector population by IRS and the only insecticide class available for the impregnation of bed nets. The present study provides the first direct evidence of the presence of the L1014F kdr mutation in Indian An. culicifacies.
The point mutation at residue L1014 leading to Leu-to-Phe change in voltage gated sodium channel is a most common mechanism associated with knock-down resistance in insects against DDT and pyrethroids. Variant replacements of this residue, L1014H and L1014S, have also been reported . In anophelines, only two mutations, L1014F and L1014S, are reported conferring resistance [7–10]. No mutation at other location has been reported in anophelines except a novel mutation, L995H, reported in a single specimen of An. culicifacies s.l. from Iran . During this study, screening of limited numbers of An. culicifacies s.l. samples from Surat didn't revealed presence of L995H mutation. The implication of such mutation in DDT/pyrethroid resistance and frequency in resistant population is yet to be investigated.
For monitoring of kdr resistance in vector population, there is a need of simple and reliable methods. AS-PCR is a convenient and inexpensive method used for kdr detection. This assay format has widely been used with the malaria vector An. gambiae in Africa [7, 10]. For An. culicifacies s.l., an AS-PCR was designed by Hoti et al  using the external primers Agd1 and Agd2 originally designed for An. gambiae. These primers do not work with An. culicifacies s.l. from Surat district (see Additional file 1) due to mismatch with template DNA which was revealed after sequencing presented in this study and, therefore, not reliable. Further, a significant departure from HW equilibrium, resulting from deficiency of heterozygotes, had been noticed by the authors, while re-examining genotyping data of kdr in An. culicifacies s.l. by Hoti et al , which may be due to discrepancy in genotyping result. This necessitated the development of culicifacies-specific molecular assays for kdr genotyping.
The AS-PCR designed for SNP detection is based on single base primer-template mismatch at 3' of the primers. These primers are generally specific but in some instances the mismatched base does not prevent the extension process [12, 13]. Mispriming with the primer specific for susceptible (kds) allele has also been reported in An. gambiae , where the mismatched base is at penultimate position from the 3' terminus. To avoid such mispriming, an ARMS method was developed, where an additional deliberate single nucleotide mismatch is incorporated near the 3' end of the allele-specific primers, which increases the specificity of the ARMS primers .
PIRA-PCR is a very reliable method for SNP detection that reduces the possibility of genotyping error due to mispriming. The PIRA-PCR was initially used for genotyping of kdr in An. gambiae by Janeira et al , where the mismatch in the PIRA primer was incorporated at the penultimate 3' base in order to create a restriction site. In the PIRA primer designed here for An. culicifacies s.l., the mismatch was incorporated at the third base position from 3' end, which seems to be advantageous, being less refractory to amplification in comparison to those where mismatch is at penultimate base. One of the most significant advantages of current PIRA-PCR over previously described methods is its capability to detect both possible codons of phenylalanine, i.e., TTT and TTC, whereas other PCR based assays can detect only one codon. Although in the single mosquito population tested only the A-to-T substitution was encountered, which is the most common substitution in mosquitoes , the possibility an A-to-C substitution in other mosquitoes cannot be ruled out and has recently been seen in Culex quinquefasciatus where both substitutions (A-to-T/C) exist . It is possible that the A-to-C substitution is being overlooked in other species. The PIRA-PCR strategy described here could be used for screening of Leu-Phe based kdr for other mosquito species with little modification in primer sequences. The PIRA-PCR is relatively costly compared to AS-PCR or ARMS, due to involvement of a restriction enzyme step, and it is also slightly more time consuming. The PIRA-PCR may be misleading in case of incomplete digestion of amplified product where RR samples may be scored as RS.
Recently several other molecular methods in addition to PIRA-PCR have been used for kdr genotyping as alternative of AS-PCR such as Heated Oligonucleotide Ligation Assays (HOLA) , Sequence-Specific Oligonucleotide Probes (SSOP) , Fluorescence Resonance Energy Transfer/Melt Curve Analysis (FRET/MCA)  and fluorescent based genotyping by capillary electrophoresis , but most of these methods are expensive or require more sophisticated instruments. PCR-based methods on the other hand are easy, inexpensive and can be performed in the most of the laboratories equipped with Thermal Cycler.
The Leu-Phe mutation (L1014F) in the voltage gated sodium channel, which generates the kdr phenotype in many insects, was detected in a pyrethroid and DDT resistant An. culicifacies s.l. population from Surat, western India. Three PCR-based methods (AS-PCR, ARMS and PIRA-PCR) were developed for kdr genotyping based on DNA sequences obtained from An. culicifacies species A, B and C. All the three assays were specific. The ARMS method was refractory to non-specific amplification in non-stringent amplification conditions. The PIRA-PCR assay is able to detect both the codons for the phenylalanine mutation at kdr locus, i.e., TTT and TTC, in a single assay, although the latter codon was not found in the population genotyped.
This collaboration was funded by the Wellcome Trust. Authors are thankful to Mr Uday Kumar and Mr Shri Bhagwan for excellent technical assistance.
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.