Presence of the knockdown resistance mutation, Vgsc-1014F in Anopheles gambiae and An. arabiensis in western Kenya
Parasites & Vectors volume 8, Article number: 616 (2015)
The voltage gated sodium channel mutation Vgsc-1014S (kdr-east) was first reported in Kenya in 2000 and has since been observed to occur at high frequencies in the local Anopheles gambiae s.s. population. The mutation Vgsc-1014F has never been reported from An. gambiae Complex complex mosquitoes in Kenya.
Molecularly confirmed An. gambiae s.s. (hereafter An. gambiae) and An. arabiensis collected from 4 different parts of western Kenya were genotyped for kdr from 2011 to 2013. Vgsc-1014F was observed to have emerged, apparently, simultaneously in both An. gambiae and An. arabiensis in 2012. A portion of the samples were submitted for sequencing in order to confirm the Vgsc-1014F genotyping results. The resulting sequence data were deposited in GenBank (Accession numbers: KR867642-KR867651, KT758295-KT758303). A single Vgsc-1014F haplotype was observed suggesting, a common origin in both species.
This is the first report of Vgsc-1014F in Kenya. Based on our samples, the mutation is present in low frequencies in both An. gambiae and An. arabiensis. It is important that we start monitoring relative frequencies of the two kdr genes so that we can determine their relative importance in an area of high insecticide treated net ownership.
The two most widely applied vector control tools, insecticide treated nets (ITNs) and indoor residual spraying (IRS) have contributed greatly to the decline in global malaria rates [1, 2]. Pyrethroids are the most commonly used insecticides in control programs due to their low human toxicity and high efficacy against vectors [3, 4]. Previously, DDT, an organochlorine, was the most widely used insecticide for vector control with its use spread out over multiple countries for malaria control [5, 6]. The widespread use of these insecticides has likely contributed to the selection of resistance across sub-Saharan Africa  (http://www.irmapper.com/).
Increased resistance to pyrethroids is particularly troubling since this is the only class of insecticides approved by WHO for use on ITNs . If ITNs are rendered ineffective, a surge in malaria transmission could follow . Resistance to pyrethroids has been reported from multiple sites in western Kenya [9, 10] with both target site and metabolic resistance mechanisms implicated [9–15]. DDT and pyrethroids function by binding to the voltage gated sodium channels (Vgsc) on the mosquito’s neurons delaying the closing of the sodium channel; prolonging the action potential and causing repetitive neuron firing, ultimately resulting in paralysis and death [8, 16].
In Anopheles gambiae s.l., knock down resistance (kdr) is commonly caused by mutations in the Vgsc- either from leucine (TTA) to phenylalanine (TTT) or leucine to serine (TCA) [11, 17] at codon 1014. Vgsc-1014S (formerly kdr-east) was first reported in Kenya in 2000 and has been observed to occur at high frequencies in the local An. gambiae populations [10, 11]. Thus far, there has been no report of the existence of Vgsc-1014F (formerly kdr-west) in Kenya but has been reported in Uganda and Tanzania in the recent past [18, 19]. Our work demonstrates the emergence of Vgsc-1014F in western Kenya in two principal malaria vectors, An. gambiae and An. arabiensis.
Material and methods
This study was conducted in four malaria endemic districts in western Kenya with two distinct Vector control interventions: Rachuonyo and Nyando where IRS (Deltamethrin in 2011 and lambdacyhalothrin in 2012) was combined with ITNs (treated with permethrin or deltamethrin); and in Bondo and Teso where only ITNs are deployed . Mosquito collections were performed annually between June and September in 2011, 2012 and 2013. Mosquito sampling, rearing and bioassays of emergent adults were conducted as described in Ochomo et al. .
Species identification & Vgsc genotyping
DNA was extracted from whole specimens and a PCR assay  was used to distinguish between An. gambiae and An. arabiensis. DNA samples were genotyped to identify the kdr genotype at amino acid position 1014 of the Vgsc using a modification of the protocol by Bass et al.,  as described in Mathias et al., .
Exon sequencing of Vgsc
Previous studies in western Kenya have only reported the presence of Vgsc-1014S mutation. Therefore, in order to confirm the presence of Vgsc-1014F and to determine if was a de novo origin, a subsample of the Vgsc-1014F carriers were Sanger sequenced. Prior to sequencing, conventional PCR was used to amplify the exon 20  which contains the 1014 locus]. Samples were sequenced at Centre for Genomic Research, University of Liverpool, UK and resulting sequences aligned using CodonCode aligner (http://www.codoncode.com/aligner/).
Analysis for the origin of Vgsc-1014F Mutation
Gene sequences obtained from the sequencing exons 20 and 27 were aligned using Codon Code aligner (http://www.codoncode.com/) and the contigs transferred to DnaSP (http://www.ub.edu/dnasp/) as FASTA files. The files were concatenated and then run using the PHASE algorithm in DnaSP . The phased files were exported as a Phylip file to TCS, a statistical parsimony software for phylogenetic network estimation (http://darwin.uvigo.es/software/tcs.html).
Frequency of Vgsc-1014S and Vgsc-1014F in the study sites from 2011 to 2013
We observed low frequencies of Vgsc-1014S in An. arabiensis,even though we had high frequencies of the same allele in An. gambiae, they were much lower than has been reported previously . We saw a simultaneous appearance of Vgsc-1014F inboth An. gambiae and An. arabiensis in 2012 in all four study sites (Table 1) and thereafter compared the mean frequencies of the genes among the three years (Table 2). Of these, 19 samples (12 An. gambiae and 7 An. arabiensis) were sequenced. 3 An. gambiae and 3 An. arabiensis were confirmed to be homozygous for Vgsc-1014F with one An. arabiensis heterozygote detected. The sequences were deposited in GenBank (Accession numbers: KR867642- KR867651, KT758295-KT758303).
Only a single 1014F haplotype was observed (Fig. 1), suggesting a common origin in the species and subsequent interspecific transfer. However it should be noted that our ability to resolve different haplotypes was constrained by the low levels of diversity at this locus and our small amplicon length (478bp).
This is the first report of Vgsc-1014F in Kenya, which appears to have emerged in both An. gambiae and An. arabiensis around 2012 and is confirmed via DNA sequencing in multiple samples. The gene has previously been observed in Uganda , then much later in Tanzania  and now in Kenya. We have developed this report to alert researchers and programmatic staff to the presence of Vgsc-1014F mutation in these two important Anopheles vectors so that they can modify their resistance marker screening procedures. It is important therefore that we start monitoring allele and genotype frequencies so that we can assess their impact in an area of high bednet ownership.
Pluess B, Tanser FC, Lengeler C, Sharp BL. Indoor residual spraying for preventing malaria. Cochrane Database Syst Rev. 2010;4:CD006657.
Killeen GF, Smith TA, Ferguson HM, Mshinda H, Abdulla S, Lengeler C, et al. Preventing childhood malaria in Africa by protecting adults from mosquitoes with insecticide-treated nets. PLoS Med. 2007;4(7):e229.
Zaim M, Aitio A, Nakashima N. Safety of pyrethroid treated mosquito nets. Med Vet Entomol. 2000;14:1–5.
Smolen M, Sang M, Lirff R. Hazards and Exposures Associated with DDT and Synthetic Pyrethroids Used for Vector Control. Geneva: World Wildlife Fund; 1999. p. 1–46.
Ratovonjato J, Randrianarivelojosia M, Rakotondrainibe ME, Raharimanga V, Andrianaivolambo L, Le Goff G, et al. Entomological and parasitological impacts of indoor residual spraying with DDT, alphacypermethrin and deltamethrin in the western foothill area of Madagascar. Malar J. 2014;13:21.
Tangena JA, Adiamoh M, D'Alessandro U, Jarju L, Jawara M, Jeffries D, et al. Alternative treatments for indoor residual spraying for malaria control in a village with pyrethroid- and DDT-resistant vectors in the Gambia. PLoS One. 2013;8(9):e74351.
Knox TB, Juma EO, Ochomo EO, Pates Jamet H, Ndungo L, Chege P, et al. An online tool for mapping insecticide resistance in major Anopheles vectors of human malaria parasites and review of resistance status for the Afrotropical region. Parasit Vectors. 2014;7:76.
Ranson H, N’Guessan R, Lines J, Moiroux N, Nkuni Z, Corbel V. Pyrethroid resistance in African Anopheline mosquitoes: what are the implications for malaria control? Trends Parasitol. 2011;27(2):91–8.
Ochomo E, Bayoh NM, Kamau L, Atieli F, Vulule J, Ouma C, et al. Pyrethroid susceptibility of malaria vectors in four Districts of western Kenya. Parasit Vectors. 2014;7:310.
Mathias D, Ochomo E, Atieli F, Ombok M, Bayoh N, Olang G, et al. Spatial and temporal variation in the kdr allele L1014S in Anopheles gambiae s.s. and phenotypic variability in susceptibility to insecticides in Western Kenya. Malar J. 2011;10:10.
Ranson H, Jensen B, Vulule J, Wang X, Hemingway J, Collins F. Identification of a point mutation in the voltage-gated sodium channel gene of Kenyan Anopheles gambiae associated with resistance to DDT and pyrethroids. Insect Mol Biol. 2000;9:491–7.
Stump A, Atieli F, Vulule J, Besansky N. Dynamics of the pyrethroid knockdown resistance allele in western Kenyan populations of Anopheles gambiae in response to insecticide-treated bed net trials. Am J Trop Med Hyg. 2004;70:591–6.
Kamau L, Agai D, Matoke D, Wachira L, Gikandi G, Vulule JM. Status of insecticide susceptibility in Anopheles gambiae sensu lato and Anopheles funestus mosquitoes from Western Kenya. J Insect Sci. 2007;8:11.
Ochomo E, Bayoh MN, Brogdon WG, Brogdon JE, Gimnig JE, Ouma C, et al. Pyrethroid resistance in Anopheles gambiae s.s. and Anopheles arabiensis in western Kenya: phenotypic, metabolic and target site characterizations of three populations. Med Vet Entomol. 2013;27:156–64.
Vulule J, Beach R, Atieli F, Roberts J, Mount D, Mwangi R. Reduced susceptibility of Anopheles gambiae to permethrin associated with the use of permethrin-impregnated bednets and curtains in Kenya. Med Vet Entomol. 1994;8:71–5.
Davies TG, Field LM, Usherwood PN, Williamson MS. DDT, pyrethrins, pyrethroids and insect sodium channels. IUBMB Life. 2007;59(3):151–62.
Martinez-Torres D, Chandre F, Williamson M, Darriet F, Berge J, Devonshire A, et al. Molecular characterization of pyrethroid knockdown resistance (kdr) in themajor malaria vector Anopheles gambiae s.s. Insect Mol Biol. 1998;7:179–84.
Verhaeghen K, Van Bortel W, Roelants P, Backeljau T, Coosemans M. Detection of the East andWest African kdr mutation in Anopheles gambiae and Anopheles arabiensis from Uganda using a new assay based on FRET/Melt Curve analysis. Malar J. 2006;5:16.
Kabula B, Kisinza W, Tungu P, Ndege C, Batengana B, Kollo D, et al. Co-occurrence and distribution of East (L1014S) and West (L1014F) African knock- down Magesa resistance in Anopheles gambiae sensu lato population of Tanzania. Trop Med Int Health. 2014;19(3):331–41.
Scott JA, Brogdon WG, Collins FH. Identification of single specimens of the Anopheles gambiae complex by polymerase chain reaction. Am J Trop Med Hyg. 1993;49:520–9.
Bass C, Nikou D, Donnelly MJ, Williamson MS, Ranson H, Ball A, et al. Detection of knockdown resistance (kdr) mutations in Anopheles gambiae: a comparison of two new high-throughput assays with existing methods. Malar J. 2007;6:111.
Pinto J, Lynd A, Elissa N, Donnelly MJ, Costa C, Gentile G. Caccone A, do Rosario VE: Co-occurrence of East and West African kdr mutations suggests high levels of resistance to pyrethroid insecticides in Anopheles gambiae from Libreville. Gabon Med Vet Entomol. 2006;20(1):27–32.
Stephens M, Smith NJ, Donnelly P. A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001;68(4):978–89.
The authors acknowledge the support of Dr. David Weetman of the Vector Group, Liverpool School of Tropical Medicine, Judith Wandera, Richard Amito, Edward Esalimba, Gideon Nyansikera and the technical support of the KEMRI/CDC Malaria Branch, Centre for Global Health Research and Centre for Biotechnology, Research and Development, KEMRI and the National Malaria Control Program. We also acknowledge the support of Dr. Tessa Knox, Dr. Abraham Mnavaza and Dr. Evan Mathenge of WHO. The research is funded by the Bill and Melinda Gates Foundation through WHO (Award #54497 awarded to Dr. Charles Mbogo). We are grateful to the Director, KEMRI for the permission to publish this data.
The authors declare that they have no competing interests.
EO, KS, NMB, LK, FA, JMV, CO, JG, MJD and CM designed and developed the study. EO, KS, BK, ER, NMB, LK, FA, JV, CO, KS, MJD and CM contributed to development of the protocol. EO, KS, ER and MJD performed gene sequencing and data analysis. EO, KS and BK performed the laboratory analysis of the samples. All authors took part in manuscript preparation, read and approved the final manuscript.
About this article
Cite this article
Ochomo, E., Subramaniam, K., Kemei, B. et al. Presence of the knockdown resistance mutation, Vgsc-1014F in Anopheles gambiae and An. arabiensis in western Kenya. Parasites Vectors 8, 616 (2015). https://doi.org/10.1186/s13071-015-1223-5
- Insecticide resistance
- Anopheles gambiae