Patterns of insecticide resistance and knock down resistance (kdr) in malaria vectors An. arabiensis, An. coluzzii and An. gambiae from sympatric areas in Senegal
Parasites & Vectorsvolume 9, Article number: 71 (2016)
Malaria vector control in Africa relies on insecticides targeting adult mosquito vectors via insecticide treated nets or indoor residual spraying. Despite the proven efficacy of these strategies, the emergence and rapid rise in insecticide resistance in malaria vectors raises many concerns about their sustainability. Therefore, the monitoring of insecticide resistance is essential for resistance management strategies implementation. We investigated the kdr mutation frequencies in 20 sympatric sites of An. arabiensis Patton, An. coluzzii Coetzee & Wilkerson and An. gambiae Giles and its importance in malaria vector control by evaluating the susceptibility to insecticides in four representative sites in Senegal.
Sibling species identification and kdr mutation detection were determined using polymerase chain reaction on mosquitoes collected using pyrethrum sprays collection in 20 sites belonging to two transects with differential insecticide selection pressure. The World Health Organization (WHO) tube test was used to determine phenotypic resistance of An. gambiae s.l. to DDT, deltamethrin, lambdacyholothrin, permethrin, bendiocarb and malathion in four representative sites.
The L1014F kdr mutation was widely distributed and was predominant in An. gambiae in comparison to An. arabiensis and An. coluzzii. The bioassay tests showed a general trend with a resistance to DDT and pyrethroids and a susceptibility to organophosphate and carbamate according to WHO thresholds. For deltamethrin and permethrin, the two most used insecticides, no significant difference were observed either between the two transects or between mortality rates suggesting no differential selection pressures on malaria vectors. The study of the KD times showed similar trends as comparable levels of resistance were observed, the effect being more pronounced for permethrin.
Our study showed a widespread resistance of malaria vectors to DDT and pyrethroids and a widespread distribution of the 1014F kdr allele. These combined observations could suggest the involvement of the kdr mutation. The existence of other resistance mechanisms could not be ruled out as a proportion of mosquitoes did not harbour the kdr allele whereas the populations were fully resistant. The susceptibility to carbamate and organophosphate could be exploited as alternative for insecticide resistance management.
Malaria vector control in Africa is based predominately on the use of residual insecticides through indoor residual spray and insecticide treated nets . Both methods have shown to be very effective against Anopheles mosquitoes [2, 3]. Pyrethroids are considered most suitable for bednets impregnation due to their insecticidal effect, relative safety for human and other mammals and their quick knock-down effect on mosquitoes, whereas other insecticide classes (organophosphates, carbamates and organochlorines) are mainly used for Indoor Residual Spraying . In some contexts, the use of both methods has significantly improved the prevention and control of malaria . However, the development of insecticide resistance has become a serious threat to the effectiveness of these control measures. One of the mechanisms involved in pyrethroid resistance in Anopheles gambiae Giles is caused by target-site insensitivity through a knock-down resistance (kdr) produced by two different points mutation at amino acid position 1014 of the voltage gated sodium channel gene. The first leads to leucine-to-phenylalanine substitution and is widely distributed in West Africa  whereas the second was described in East Africa and involves a leucine-to-serine change . To date, the latter is mainly found with the 1014F kdr allele and is suspected to be less involved in pyrethroid resistance than the 1014F allele [7, 8]. Within the An. gambiae complex, a sympatric ecological diversification is in progress and has leaded to the emergence of at least two incipient species (the M and S molecular forms). These forms repeatedly showed heterogeneous levels of divergence in most parts of Africa [9, 10], and are now recognised as separate species - An. coluzzii Coetzee & Wilkerson and An. gambiae respectively for M and S forms . In their “far west” African distribution, genetic isolation was observed along the Senegambian coasts  as well as in inland areas of south-eastern Senegal, where substantial reproductive isolation was evident and further supports the ongoing process of speciation in inland areas .
Previous studies have revealed the presence of the kdr mutation in An. gambiae and its absence in An. coluzzii populations even in sympatric areas where the mosquitoes are presumably subjected to similar selection pressures. These observations supported the possibilities of different resistance mechanism in An. coluzzii populations or a restriction in gene flow between the two species [14–17] that could therefore affect the spread of the kdr mutation. Subsequent studies conducted thereafter showed the presence of this mutation in An. coluzzii populations in several sites in Africa. This was first reported in Benin and was attributed to the existence of gene flow from An. gambiae to An. coluzzii rather than the same mutation that occurred in both species  and further in Cameroon .
Based on the presence of An. arabiensis Patton, An. coluzzii and An. gambiae in sympatric areas in south-eastern Senegal, we assessed the kdr frequencies in 20 sites and further studied its importance in terms of vector control by evaluating the susceptibility to insecticides in four representative sites.
Study design and study sites
The study was conducted in the Tambacounda region of southeastern Senegal (Fig. 1) during the rainy seasons, 2010 and 2011, which last from June to October, with a peak in August–September, and an annual rainfall of 500 mm.
In 2010, 20 sites were selected from two distinct transects. The first is situated along the Gambia River an agricultural areas growing rice and cash crops of banana with high pesticide usage. The second transect is along the National Road 7, in a less humid and mostly arid area with a little or no application of pesticides. In each site, indoor resting mosquitoes were collected during the daytime using pyrethrum sprays. Upon collection, mosquitoes were counted and identified as An. gambiae s.l. morphologically using the keys of Gillies & de Meillon . All the mosquito samples were stored individually after identification in numbered vials containing desiccant until laboratory processing.
In 2011, Anopheles larvae and pupae were collected from natural breeding sites in four selected sites (two in the first transect namely Koar and Sankagne and two in the second transect namely Djnkore and Wassadou). Upon collection, they were kept in separate labelled bottles, transported to the insectary and maintained at a relative humidity of 75 ± 5 % and a temperature of 28 ± 3 °C.
kdr molecular genotyping and species identification
For molecular identification, genomic DNA was extracted from the wings or legs of individual mosquitoes as described by Collins et al. . For each collection of mosquitoes randomly sampled from each village, 30–100 % of females belonging to the An. gambiae s.l. were identified to species using the molecular methods of Favia et al.  and Fanello et al. . The presence of L1014F mutation was confirmed using the method described by Martinez-Torres et al. .
Insecticide susceptibility tests
Insecticide susceptibility tests were conducted on unfed adult females aged from 2 to 5 days. Bioassays were carried out using WHO test kits for adult mosquitoes. Insecticide-impregnated papers were provided by The Vector Control Research Unit, School of Biological Sciences (Universiti Sains Malaysia), a WHO Collaborating Centre. The following diagnostic concentrations of insecticides were tested: 0.05 % deltamethrin, 4 % DDT, 0.75 % permethrin, 0.05 % lambdacyhalothrin, 5 % malathion and 0.1 % bendiocarb. For each insecticide, four replicates were exposed for 60 min. The mosquitoes were then transferred in tubes with untreated papers and kept under observation for 24 h.
The kdr allele frequency was estimated for each site and period of collection as the proportion of specimens found with the L1014F kdr alleles. The frequencies as well as the conformity to Hardy–Weinberg equilibrium were assessed using Genepop v.3.2. .
For each insecticide, the resistance status was studied using WHO criteria . For DDT and pyrethroid insecticides, knock-down rates (%) were assessed at 10, 15, 20, 30, 40, 50 and 60 min.
All statistical analyses were performed using R software (version 3.0.2).
Frequency and distribution of 1014F allele
In this study, the kdr mutation was confirmed in An. gambiae s.l. populations in all of the 20 sites selected (frequency range 0.05–0.33 %) with significant variations (χ2 = 76.3, df = 19, p < 0.001). Examination of 189 An. arabiensis, 115 An. coluzzii and 814 An. gambiae showed that only An. gambiae s.s. was carrying the kdr mutation in all the 20 sites (frequency range: 0.06–0.38 %). The kdr mutation was observed in 17 sites (frequency range in observed sites: 0.04–0.40 %) and 12 sites (frequency range in observed sites: 0.06–0.50 %) in An. arabiensis and An. coluzzii, respectively (Table 1, Fig. 2). The mean kdr frequencies by species were 14.38 ± 2.42 (An. arabiensis), 21.22 ± 3.85 (An. gambiae) and 14.69 ± 1.95 (An. coluzzii). No significant differences were observed between these means (one-way ANOVA, F = 1.83, p = 0.17).
These frequencies were significantly different only between An. gambiae and An. arabiensis (χ2 = 10.92, df = 1, p = 0.0009), irrespective of the transects. When considering each species within each different transect, the prevalence of the 1014F kdr allele was similar between all three species in each transect. In An. coluzzii, the 1014F kdr allele was the only allele observed in two of the nine localities prospected in transect two (Madina Dian and Bira) .
In all the bioassays, controls using wild An. gambiae s.l. populations from each site showed mortality rate less than 5 %, thus, no corrections were required in the test sample data.
All populations exhibited resistance to DDT (Fig. 3). The highest mortality rate was observed for Djnkore populations (mortality rate = 74.07 %, IC95 = 64.75–82.03 %). For pyrethroids, resistance was observed for both types: Type I (permethrin) and Type II (deltamethrin, lambdacyhalothrin) for all populations except in Sankagne (deltamethrin) and probably in Wassadou (permethrin) and Koar (lambdacyhalothrin) where a probable resistance is observed with the mortality rates of 94 % (IC95 = 87.4–97.77 %) and 94.12 % (IC95 = 87.64–97.81 %) for the two latter villages. Full susceptibility was observed for the organophosphate malathion (Koar and Sankagne) and for the carbamate bendiocarb (Koar, Sankagne and Djnkore). Wassadou populations exhibited resistance to bendiocarb with a mortality rate of 95.87 % (IC95 = 90.62–98.64 %).
Although a statistically significant difference was observed between mortality rates for sites showing confirmed resistance for DDT (χ2 = 22.3, df = 3, p < 0.001), no significant differences were observed for deltamethrin and permethrin, whereas mortality rates to lambdacyhalothrin were significantly different between Sankagne and Djnkore (χ2 = 4.4, df = 1, p = 0.03).
Knock down effects of pyrethroids
Despite the resistant status observed, the study of the knock down effects of DDT and pyrethroids showed different patterns. For DDT, the resistance observed was associated to a low effect on the four populations tested (Fig. 4). Indeed after 60 min exposure the highest knock down effect was below 75 % (population from Djnkore). For the other populations, the highest values were respectively 57.01, 53.43 and 31.35 % at 60 mn for Koar, Sankagne and Wassadou (Fig. 4). For these populations, KD50 and KD95 increased markedly with a near complete loss of knock down effect for the Wassadou population (Table 2). Despite the resistance status observed for Koar, a similar trend of the KD rates was observed for deltamethrin in Sankagne (p < 0.05) with 100 % knock down at 20 min (Fig. 5). For lambdacyhalothrin, the KD dynamics were quite similar for Koar, Sankagne and Djnkore populations (Fig. 6). While the KD50 times were similar, a slight increase of the KD95 was observed in the resistant population of Djnkore, but no significant difference was observed with the populations of Koar and Sankagne (p < 0.05 after Bonferroni correction). For permethrin, the pattern of KD effect is similar to the resistance status (Fig. 7). Whatever the population pairs considered, the difference was statistically significant (p < 0.05).
The present study compared the patterns of insecticide resistance and kdr resistance in three sympatric species of the An. gambiae complex. The study of the distribution of the L1014F kdr mutation showed its widespread presence in the study sites particularly in An. gambiae in comparison to An. arabiensis and An. coluzzii as observed elsewhere in Africa [7, 8, 25]. In Senegal, previous reports have documented the presence of the L1014F kdr allele in An. gambiae populations from the southeastern part of the country with frequencies ranging from 14.1 to 18.6 %  and in An. arabiensis populations from two suburbs of Dakar . Compared to the present study, these observed frequencies are low with respect to some of our sites herein, where the L1014F kdr mutation frequencies were up to 38 % in An. gambiae, 50 % in An. coluzzii and 40 % in An. arabiensis. However, these frequencies seem to be low when compared with other sites in Benin where frequencies of up to 80 % were recently observed . The comparison of the observed frequencies with those expected under Hardy-Weinberg equilibrium indicated in some cases (mainly in An. gambiae and An. arabiensis populations), a deficit or excess of heterozygotes. While the origin of the L1014F mutation in An. arabiensis is suspected to be a new and independent mutation, its presence in An. coluzzii has been suggested to occur by introgression from An. gambiae . Therefore, it is undoubtedly important to carry out additional studies to unravel the origin of this mutation in An. arabiensis and An. coluzzii in our context as other mechanisms could be involved.
Globally, the results of the bioassay tests showed that the populations were resistant to DDT and pyrethroids and susceptible to organophosphate and carbamate. This situation is relatively common in many sites in Africa . In West Africa, pyrethroid resistance is high, widespread and predominant in An. gambiae compared to An. arabiensis . Concerning carbamate and organophospahe insecticides, our results contrast with recent findings that showed high levels of resistance to these insecticides in the urban site of Dakar .
The resistance levels observed to pyrethroids varied greatly between the four sites studied. For deltamethrin and permethrin the lack of difference in mortality rates could reflect no differential selection pressures on malaria vectors. This could be due to the fact that these insecticides are the most widely used in insecticides-treated nets distributed by the National Malaria Control Program in Senegal. Indeed, Senegal is on the short list of African countries that have reached the RBM target of 80 % of households owning at least one insecticide-treated nets . Despite the resistance observed, the study of the KD dynamic showed a similar effect of deltamethrin on the resistant population of Koar compared to the susceptible population of Sankagne where a similar trend was observed. In contrast, the study of the effect of permethrin showed profiles similar to observed resistance. Except the fact that deltamethrin and permethrin are different chemical molecules (type I vs type II), this observation may express a more recent use of deltamethrin in the area in comparison with permethrin. Indeed, impregnated-bednets (Olyset Net©) were introduced in the study area in the 1990’s  and could have exerted strong insecticidal pressures on mosquitoes giving rise to the high levels of insecticide resistance observed. However, for lambdacyhalothrin, the difference observed between Koar and Djnkore reflects variation in resistance selection pressures, which results from differential use in agriculture representing the main human activities.
During this study a widespread resistance to DDT and pyrethroids and a widespread 1014F kdr allele distribution was observed. The observed resistance to DDT and pyrethroids combined with the high KD times could suggest the involvement of the 1014F kdr mutation, as observed elsewhere [34–37]. However, the involvement of other resistance mechanisms including the 1014S kdr mutation recently observed in West Africa and Senegal [31, 38, 39], could not be ruled out, as at least 50 % of the populations did not harbour the 1014F allele whereas the populations were fully resistant. This needs further investigations. The susceptibility to carbamate and organophosphate could be an alternative for insecticide resistance management.
WHO. Test procedures for insecticide resistance monitoring in malaria vector mosquitoes. Geneva: World Health Organization; 2013.
Curtis CF, Mnzava AE, Misra S, Rowland M. Bednets or spraying? Summary of the presentations and the discussion. Trans R Soc Trop Med Hyg. 1999;93:460.
Nyarango PM, Gebremeskel T, Mebrahtu G, Mufunda J, Abdulmumini U, Ogbamariam A, et al. A steep decline of malaria morbidity and mortality trends in Eritrea between 2000 and 2004: the effect of combination of control methods. Malar J. 2006;5:33.
WHO. Indoor residual spraying. An operational manual for indoor residual spraying (IRS) for malaria transmission control and elimination. Geneva: World Health Organization; 2013.
Martinez-Torres D, Chandre F, Williamson MS, Darriet F, Berge JB, Devonshire AL, et al. Molecular characterization of pyrethroid knockdown resistance (kdr) in the major malaria vector Anopheles gambiae s.s. Insect Mol Biol. 1998;7:179–84.
Ranson H, Jensen B, Vulule JM, Wang X, Hemingway J, Collins FH. 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.
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:91–8.
Corbel V, N’Guessan R. Distribution, mechanisms, impact and management of insecticide resistance in malaria vectors. A pragmatic review. Edited by Manguin S. Rijeka, Croatia: In Tech Publisher; 2013. p. 579–633.
Neafsey DE, Lawniczak MK, Park DJ, Redmond SN, Coulibaly MB, Traore SF, et al. SNP genotypying defines complex gene-flow boundaries among African malaria vector mosquitoes. Science. 2010;330:514–7.
Reidenbach KR, Neafsey DE, Costantini C, Sagnon NF, Simard F, Ragland GJ, et al. Patterns of genomic differentiation between ecologically differentiated M and S forms of An. gambiae in West and Central Africa. Genome Biol Evol. 2012;4:1202–12.
Coetzee M, Hunt RH, Wilkerson R, della Torre A, Coulibaly MB, Besansky NJ. Anopheles coluzzii and Anopheles amharicus, new members of the Anopheles gambiae complex. Zootaxa. 2013;3:246–74.
Nwakanma DC, Neafsey DE, Jawara M, Adiamoh M, Lund E, Rodrigues A, et al. Breakdown in the Process of Incipient Speciation in Anopheles gambiae. Genetics. 2013;193:1221–31.
Niang EA, Konaté L, Diallo M, Faye O, Dia I. Reproductive Isolation among Sympatric Molecular Forms of An. gambiae from Inland Areas of South-Eastern Senegal. PLoS One. 2014;9:e104622.
della Torre A, Fanello C, Akogbeto M, Dossou-Yovo J, Favia G, Petrarca V, et al. Molecular evidence of incipient speciation within Anopheles gambiae s.s. in West Africa. Insect Mol Biol. 2001;10:9–18.
Wondji CS, Simard F, Fontenille D. Evidence for genetic differentiation between the molecular forms M and S within the Forest chromosomal form of Anopheles gambiae in an area of sympatry. Insect Mol Biol. 2002;11:11–9.
Fanello C, Petrarca V, della Torre A, Santolamazza F, Dolo G, Coulibaly M, et al. The pyrethroid knock-down resistance in the Anopheles gambiae complex in Mali and further indication of incipient speciation within An. gambiae s.s. Insect Mol Biol. 2003;12:241–5.
Awolola TS, Oyewole IO, Amajoh CN, Idowu ET, Ajayi MB, Oduola A, et al. Distribution of the molecular forms of Anopheles gambiae and pyrethroid knock down resistance gene in Nigeria. Acta Trop. 2005;95:204–9.
Weill M, Chandre F, Brengues C, Manguin S, Akogbeto M, Pasteur N, et al. The kdr mutation occurs in the Mopti form of Anopheles gambiae s.s. through introgression. Insect Mol Biol. 2000;9:451–5.
Etang J, Vicente JL, Nwane P, Chouaibou M, Morlais I, Do Rosario VE, et al. Polymorphism of intron-1 in the voltage-gated sodium channel gene of Anopheles gambiae s.s. populations from Cameroon with emphasis on insecticide knockdown resistance mutations. Mol Ecol. 2009;18:3076–86.
Gillies MT, de Meillon B. The anophelinae of Africa, south of the Sahara. Johannesburg The South African Institute for Medical Research. 1068; 54: 343 pages.
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.
Favia G, della Torre A, Bagayoko M, Lanfrancotti A, Sagnon N, Toure YT. Molecular identification of sympatric chromosomal forms of Anopheles gambiae and further evidence of their reproductive isolation. Insect Mol Biol. 1997;6:377–83.
Fanello C, Santolamazza F, della Torre A. Simultaneous identification of species and molecular forms of the Anopheles gambiae complex by PCR-RFLP. Med Vet Entomol. 2002;16:461–4.
Raymond M, Rousset F. GENEPOP Version 1.2. A population genetics software for exact tests and ecumenicism. J Hered. 1995;86:248–9.
Dabire RK, Namountougou M, Diabate A, Soma DD, Bado J, Hyacinthe K, et al. Distribution and Frequency of kdr Mutations within Anopheles gambiae s.l. Populations and First Report of the Ace.1G119S Mutation in Anopheles arabiensis from Burkina Faso (West Africa). PLoS One. 2014;9:e101484.
Santolamazza F, Calzetta M, Etang J, Barrese E, Dia I, Caccone A, et al. Distribution of the knock-down resistance mutations in Anopheles gambiae molecular forms in west and west-central Africa. Malaria J. 2008;7:74.
Pagès F, Texier G, Pradines B, Gadiaga L, Machault V, Jarjaval F, et al. Malaria transmission in Dakar: A two-year survey. Malaria J. 2008;7:178.
Gnanguenon V, Agossa FR, Badirou K, Govoetchan R, Anagonou R, Oke-Agbo F, et al. Malaria vectors resistance to insecticides in Benin: current trends and mechanisms involved. Parasit Vectors. 2015;8:223.
Diabate A, Brengues C, Baldet T, Dabire KR, Hougard JM, Akogbeto M, et al. The spread of the Leu-Phe kdr mutation through Anopheles gambiae complex in Burkina Faso: genetic introgression and de novo phenomena. Trop Med Int Health. 2004;9:1267–73.
Chandre F, Darriet F, Manga L, Akogbeto M, Faye O, Mouchet J, et al. Status of pyrethroid resistance in Anopheles gambiae sensu lato. Bull World Health Organ. 1999;77:230–4.
Ndiath MO, Cailleau A, Orlandi-Pradines E, Bessell P, Pagès F, Trape JF, et al. Emerging knock-down resistance in Anopheles arabiensis populations of Dakar, Senegal: first evidence of a high prevalence of kdr-e mutation in West African urban area. Malar J. 2015;14(1):364.
Thwing JI, Perry RT, Townes DA, Diouf MB, Ndiaye S, Thior M. Success of Senegal’s first nationwide distribution of long-lasting insecticide-treated nets to children under five - contribution toward universal coverage. Malaria J. 2011;10:86.
Faye O, Konaté L, Gaye O, Fontenille D, Sy N, Diop A, et al. Impact de l’utilisation des moustiquaires pré-imprégnées de perméthrine sur la transmission du paludisme dans un village hyperendémique du Sénégal. Med Trop. 1998;58:355–60.
Tia E, Akogbeto M, Koffi A, Toure M, Adja AM, Koné M, et al. Situation de la résistance d’Anopheles gambiae s.s. (Diptera : Culicidae) aux pyréthrinoïdes et au DDT dans cinq écosystèmes agricoles de Côte-d’Ivoire. Bull Soc Pathol Exot. 2006;99:278–82.
Dabiré KR, Diabaté A, Namountougou M, Toé KH, Ouari A, Kengne P, et al. Distribution of pyrethroid and DDT resistance and the L1014F kdr mutation in Anopheles gambiae s.l. from Burkina Faso (West Africa). Trans R Soc Trop Med Hyg. 2009;103:1113–20.
Kwiatkowska RM, Platt N, Poupardin R, Irving H, Dabire RK, Mitchell S, et al. Dissecting the mechanisms responsible for the multiple insecticide resistance phenotype in Anopheles gambiae s.s. M form, from Vallée du Kou. Burkina Faso. Gene. 2013;519:98–106.
Okorie PN, Ademowo OG, Irving H, Kelly-Hope LA, Wondji CS. Insecticide susceptibility of Anopheles coluzzii and Anopheles gambiae mosquitoes in Ibadan, Southwest Nigeria. Med Vet Entomol. 2015;29:44–50.
Djegbe I, Boussari O, Sidick A, Martin T, Ranson H, Chandre F, et al. Dynamics of insecticide resistance in malaria vectors in Benin: first evidence of the presence of L1014S kdr mutation in Anopheles gambiae from West Africa. Malar J. 2011;10:261.
Badolo A, Traore A, Jones CM, Sanou A, Flood L, Guelbeogo WM, et al. Three years of insecticide resistance monitoring in Anopheles gambiae in Burkina Faso: resistance on the rise? Malar J. 2012;11:232.
Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution. 1984;38:1358–70.
This study was funded by the EU project INFRAVEC (Research capacity for the implementation of genetic control of Mosquitoes) and the European Commission’s Seventh Framework Research Programme under the grant agreement 228421).
The authors declare that they have no competing interests.
ID, EAN, MD conceived and designed the study. EAN and ID carried out the field collections, performed the experiments and analysed the data. EAN, LK, OF, MD and ID drafted the manuscript and read, critically revised and approved the final manuscript.