- Open Access
Status of insecticide resistance and its biochemical and molecular mechanisms in Anopheles stephensi (Diptera: Culicidae) from Afghanistan
Malaria Journal volume 18, Article number: 249 (2019)
Insecticide resistance of Anopheles stephensi, the main malaria vector in eastern Afghanistan, has been reported previously. This study describes the biochemical and molecular mechanisms of resistance to facilitate effective vector control and insecticide resistance management.
Mosquito larvae were collected from the provinces of Kunar, Laghman and Nangarhar from 2014 to 2017. The susceptibility of the reared 3–4 days old adults was tested with deltamethrin 0.05%, bendiocarb 0.1%, malathion 5%, permethrin 0.75% and DDT 4%. Cytochrome P450 content and general esterase, glutathione S-transferase (GST) and acetylcholinesterase (AChE) activities were measured in the three field populations and the results were compared with those of the laboratory susceptible An. stephensi Beech strain. Two separate allele-specific PCR assays were used to identify L1014, L1014F and L1014S mutations in the voltage gated sodium channel gene of An. stephensi. Probit analysis, ANOVA and Hardy–Weinberg equilibrium were used to analyse bioassay, biochemical assay and gene frequency data respectively.
The population of An. stephensi from Kunar was susceptible to bendiocarb, apart from this, all populations were resistant to all the other insecticides tested. The differences between all values for cytochrome P450s, general esterases, GSTs and AChE inhibition rates in the Kunar, Laghman and Nangarhar populations were statistically significant when compared to the Beech strain, excluding GST activities between Kunar and Beech due to the high standard deviation in Kunar. The three different sodium channel alleles [L1014 (wild type), L1014F (kdr west) and L1014S (kdr east)] were all segregated in the Afghan populations. The frequencies of kdr east mutation were 22.9%, 32.7% and 35% in Kunar, Laghman and Nangarhar populations respectively. Kdr west was at the lowest frequency of 4.44%.
Resistance to different groups of insecticides in the field populations of An. stephensi from Kunar, Laghman and Nangarhar Provinces of Afghanistan is caused by a range of metabolic and site insensitivity mechanisms, including esterases, cytochrome P450s and GSTs combined with AChE and sodium channel target site insensitivity. The intensity and frequency of these mechanisms are increasing in these populations, calling for urgent reorientation of vector control programmes and implementation of insecticide resistance management strategies.
Malaria is a major endemic vector borne disease in Afghanistan. The 2018 World Malaria Report states that 27, 50 and 23% of the total population of 35,530,083 Afghans are at high, low and no risk of malaria, respectively . Although there are 6 potential malaria vector species, Anopheles stephensi is the major vector in the eastern provinces of Kunar, Laghman and Nangarhar [2,3,4]. This species is common in the Middle East and the Indian subcontinent extending to South China and Myanmar .
In Afghanistan, the malaria control is reliant on indoor residual spraying (IRS) with deltamethrin, or more recently bendiocarb, combined with distribution of LNs mostly PermaNet [6, 7]. During 2007–2016 over 2 million deltamethrin-treated LNs were distributed in Kunar, Laghman and Nangarhar , with a top up distribution of 45,000 LNs in Kunar and Laghman Provinces in 2017 .
Insecticide resistance in An. stephensi to all four classes of insecticides including DDT, malathion, bendiocarb, permethrin and deltamethrin in Kunar, Laghman and Nangarhar Provinces has been reported, although it remained susceptible to bendiocarb until 2014 in Afghanistan [8, 9]. Resistance to DDT, dieldrin, malathion and more recently pyrethroids has also been reported in An. stephensi from the Middle East and the Indian subcontinent [6, 9,10,11,12,13,14,15,16]. Selection of resistance to bendiocarb in An. stephensi was recently reported from Afghanistan .
Several mechanisms, including metabolic resistance and site insensitivity can cause insecticide resistance [12, 15, 17,18,19,20,21,22,23,24,25,26,27]. In a previous study on An. stephensi from Afghanistan, general esterases (GES), glutathione S-transferases (GSTs), cytochrome P450s and insensitive acetylcholinesterase (iAChEs) were implicated in insecticide resistance . Pyrethroid insecticide resistance in An. stephensi from India and Iran was associated with increased activity of GES and GSTs [12, 17, 20, 21, 29]. Involvement of GSTs in insecticide resistance is evident in many insects, including mosquitoes [17, 18, 24].
Knockdown resistance (kdr) mutation is widespread in Anopheles species in Africa especially Anopheles gambiae [30,31,32,33,34,35]. Originally, the L1014F mutation (later known as kdr west) was detected. In 2000, a second kdr mutation (kdr east) was detected in Kenyan An. gambiae , subsequently both mutations have been detected in other Anopheles conferring varying degrees of phenotypic pyrethroid resistance [5, 10, 37,38,39,40,41,42,43].
The first report of a kdr L1014F resistance mechanism in An. stephensi was in the DUB-S strain in 2003 . Recently, kdr east and kdr west mutations have been detected in An. stephensi from India . In 2014, target site insensitivity for pyrethroid insecticides was studied in An. stephensi from Kunar and Nangarhar Provinces of Afghanistan . The wild type susceptible L1014 allele in the sodium channel gene was the most prevalent followed by L1014S (kdr east, 21.4%) and L1014F (kdr west, 1.4%), no kdr homozygotes were collected.
The World Health Organization (WHO) recommends that insecticide susceptibility status of malaria vectors should be monitored annually [44, 45]. When insecticide resistance is detected, its intensity and the biochemical and molecular mechanisms should also be investigated [44, 45]. Accurate information on the underlying resistance mechanisms and their intensity or frequency in malaria vectors can then inform vector control programmes and ensure timely management of insecticide resistance. Following the WHO Global Plan for Insecticide Resistance Management (GPIRM) , in this study the insecticide resistance status and its underlying mechanisms were investigated in An. stephensi from Kunar, Laghman and Nangarhar Provinces in eastern Afghanistan.
The study areas were the provinces of Kunar (34.8466°N, 71.0973°E), Laghman (34.6898°N, 70.1456°E) and Nangarhar (34.1718°N, 70.6217°E) in East Afghanistan (Fig. 1), sample site details are given in Table 1.
Larval collection and mosquito rearing
Thousands of larvae were collected from multiple breeding sites in the provinces of Kunar, Laghman and Nangarhar from 2014 to 2017. Breeding sites with the highest larval density were used for sampling to obtain enough specimens for susceptibility tests and biochemical and molecular analyses.
Field collected larvae were reared to the adult at 25 ± 2 °C temperature and 75 ± 10% relative humidity. The adults were identified to species using Glick’s identification keys . Sugar-fed 3–4 days old adult An. stephensi mosquito specimens were used for bioassays and biochemical and molecular analyses. The susceptible Beech strain of An. stephensi was provided by the Department of Medical Entomology, School of Public Health, Tehran University of Medical Sciences, Iran.
Insecticide susceptibility tests
Insecticide susceptibility tests were carried out according to the standard WHO procedure for mosquito adults . The WHO supplied insecticide-impregnated papers of DDT 4%, malathion 5%, bendiocarb 0.1%, permethrin 0.75% and deltamethrin 0.05%, were used for bioassays. Mosquitoes were divided into batches of 25 before being exposed to the insecticide-treated papers for 1 h. Experiments were conducted under insectary conditions with a minimum of four replicates per bioassay. For control replicates, silicone-treated papers were used. The results of the bioassay were discarded if the mortality of the control replicates was over 20%. Abbott’s formula was used to correct the mortality if it was between 5 and 20% .
The biochemical assays were performed according to the protocol of WHO/WHOPES . The enzyme activity of GSTs and GES as well as the cytochrome P450s content and inhibition rates of AChE (using propoxur) were measured.
DNA from individual mosquitoes was extracted using the Livak buffer extraction method  with some modifications. Livak buffer contained 80 mM NaCl, 1.57% Tris Base, 0.5% sodium dodecyl sulfate (SDS), 5.5% sucrose and 50 mM EDTA. Individual mosquitoes were homogenized in 100 µl pre-heated (65 °C) Livak buffer in 1.5 ml Eppendorf tubes using a plastic pestle. Homogenates were incubated at 65 °C for 30 min. Potassium acetate was added to each tube to a final concentration of 1 M before incubating the mixture on ice for 30 min. The tubes were centrifuged at 12,000g for 15 min at 4 °C. Supernatants were transferred to clean tubes and mixed with 200 µl ice-cold ethanol, followed by centrifugation at 12,000g for 15 min at 4 °C. Pellets were rinsed in 100 µl 70% ice-cold ethanol, spun at 12,000g for 5 min at 4 °C and re-suspended in 50 µl pre-heated Tris–EDTA (TE) buffer or nuclease free water.
Primer design and PCR
To obtain a fragment of voltage gated sodium channel (vgsc) encompassing the kdr locus, a DNA section from IIS6 segment of An. stephensi was amplified using primers kdrF (5′-GGA CCA YGA TTT GCC AAG AT-3′) and kdrR (5′-TGG TGC AGA CAA GGA TGA AG-3′) in a reaction mixture (25 μl) that contained 1× buffer, 1.5 mM of MgCl2, 200 μM of each dNTP, 0.5 μM of each primer and 0.625 unit of Taq DNA polymerase. The conditions of PCR were: an initial denaturation at 95 °C for 5 min, followed by 35 cycles at 95 °C for 30 S, 48 °C for 30 S and 72 °C for 45 S, and a final extension step at 72 °C for 7 min. The PCR products were purified with Takapozist PCR purification kit (Takapozist, Iran) and were sequenced in both directions by Macrogen Inc, Korea using BigDye (Applied Biosystem Chemistry). A schematic diagram showing the stretch of the IIS6 of the sodium channel gene where the primers for cloning and screening sit is given in Fig. 2.
Positive controls and cloning
Positive control for kdr east mutation (L1014S) was the courtesy of Dr. OM Singh, National Malaria Research Center, India. Positive control for the kdr west mutation (L1014F) was from initial PCR screening of some specimens that survived bioassay using kdrF and kdrR primers. Wild type positive control (L1014) was the product of PCR of two specimens from the susceptible Beech strain of An. stephensi with aforementioned primers. All positive control fragments were confirmed by sequencing and then cloned in pTG19-T system (Sinaclon, Iran) according to the manufacturer’s instruction and used as positive control in genotyping reactions.
Genotyping kdr alleles
The screening method developed by Singh et al. was used in this study with some modifications . Two PCR reactions were carried out for each DNA sample from individual mosquito specimens. In both PCRs, the forward primer used was kdrF. In the first PCR (hereafter called PCR-1), the allele 1014F is discriminated from the wild type (L1014) and in the second PCR (hereafter called PCR-2), 1014S and wild type alleles are discriminated. The allele-specific primers were the same used by Singh et al. . They are St-L/SR (5′-GCG GGC AGG GCG GCG GGG GCG GGG CCC GAT CGG AAA GTA AGT TAC TTA CGt CT-3′) and St-PheR (5′-GAT CGG AAA GTA AGT TAC TTA CGg CA-3′) for PCR-1, and St-LeuR (5′-GCG GGC AGG GCG GCG GGG GCG GGG CCC GAT CGG AAA GTA AGT TAC TTA CGA gTA-3′) and St-SerR (5′-CGA TCG GAA AGT AAG TTA CTT ACG AtT G-3′) for PCR-2. The expected amplicon size formed by the allele-specific primers St-L/SR and St-PheR (with primer kdrF) in PCR-1 are 218 and 191 bp receptively. The expected size of amplicons in PCR-2 with allele specific primers St-LeuR and St-SerR are 218 and 192 bp, respectively.
Data management, transformation and analyses
The reading of the activity/content of the enzymes were done in a UV/visible microtitre plate reader (BioTek, USA) run under KC junior software and the data were directly extracted to Microsoft Excel for further analysis. Mean values of activity or contents of each enzyme of all populations were compared by ANOVA in conjunction with the Tukey’s statistical test using SPSS version 19 software. Enzyme ratios (ER) were calculated by dividing the mean activities or content of the enzymes in the field populations with those of the Beech susceptible strain.
The molecular analyses data were calculated based on the frequency of kdr mutations in An. stephensi populations from different provinces in dead or alive specimens following bioassays. Hardy Weinberg equilibrium software was used for analysis between different kdr genotypes .
Approximately 2250, 1800 and 2250 3–4 days old adult mosquitoes were used for bioassays in 2014, 2015 and 2017, respectively. Two hundred 3–4 days old adult mosquito specimens were used for biochemical assays. Dead and alive adult mosquitoes following bioassays with deltamethrin from each population were used for molecular analysis.
Considering the susceptibility threshold of a mortality rate above 98%, then only one population of An. stephensi from Kunar was susceptible to bendiocarb and the Laghman and Nangarhar populations were resistant to all insecticides tested. When bioassay mortalities in 2014 were compared to those of 2017, there was a reduction of 27% in susceptibility to permethrin in Kunar population, while the susceptibility to permethrin in the other two field populations remained more or less the same. Susceptibility to deltamethrin was the same in 2014 and 2017 in Kunar population, whereas an increase of about 30% in susceptibility to deltamethrin was observed in the other two field populations. Regarding susceptibility to malathions, reduction of 55%, 46% and 41% were monitored in the Kunar, Laghman and Nangarhar populations, respectively. The susceptibility to bendiocarb increased 3% in the Kunar population, changing the status of this population from being resistant to susceptible. The susceptibility to bendiocarb in the other two field populations decreased 43.5% and 16%. Surprisingly, the susceptibility to DDT increased 55% in the Kunar population, however, in Laghman population it reduced about 44% while it remained the same in the Nangarhar population (Table 2).
The esterase activity ratios of An. stephensi from Kunar, Laghman and Nangarhar Provinces in East Afghanistan compared to the susceptible Beech strain are given in Table 3. The mean activity of alpha- and beta-naphthyl acetate were 0.000983 and 0.000961 µM/min/mg protein in the Kunar population, 0.001 and 0.001 µM/min/mg protein in the Laghman population and 0.000911 and 0.000904 µM/min/mg protein in the Nangarhar population and 0.000517 and 0.000567 µM/min/mg protein in the susceptible Beech strain. The esterase activities of mosquitoes from the Kunar, Laghman and Nangarhar populations were statistically significantly higher than those of the susceptible Beech strain at the 5% level. The esterase activity in the Laghman population was marginally higher than those of the Kunar and Nangarhar populations (Table 3).
The activity of the GSTs in the Laghman and Nangarhar populations were significantly higher than that of the susceptible Beech strain, however, the differences between the GST activity of Kunar population was not significantly different from that of the susceptible strain (Table 3).
Cytochrome P450s contents
The ratio of cytochrome P450s in the Kunar, Laghman and Nangarhar populations were 2.68, 2.37 and 2.45 when compared with that of the susceptible Beech strain (Table 3). However, the differences of the cytochrome P450 contents between the Kunar, Laghman and Nangarhar populations were not statistically significant (Table 3).
The AChE inhibition rates with the carbamate propoxur were 68.63% in the susceptible Beech strain, 51.98% in the Kunar population, 48.82% in the Laghman population and 53.14% in the Nangarhar population (Table 3). The inhibition levels in all field populations were lower than the threshold of 60% set for considering the AChE insensitive to propoxur. There were significant differences between the three populations in AChE inhibition when compared with that of the susceptible Beech strain (p > 0.001). However, the differences between the inhibition rates of AChE in the field populations were not statistically significant from each other (p = 0.09). The frequency of individuals with iAChE (inhibition less than 60% with propoxur) were 37%, 64%, 78% and 75% in susceptible, Kunar, Laghman and Nangarhar populations of An. stephensi, respectively.
In summary, the differences between the activities/contents of all the enzymes measured in this study in the Kunar, Laghman and Nangarhar populations were statistically significant compared with those of the susceptible Beech strain. The Laghman population had a marginally higher enzymes activities/content compared with the other Afghan field populations.
Primers designed to amplify a flanking region encompassing the SII6 region containing the kdr codon of the vgsc of An. stephensi gave a 493 bp product, the results of sequencing this region are shown in Fig. 2. The cloned products of the kdr alleles were used as positive control in genotyping reactions.
Using the forward primer designed in this study and four reverse primers, a total of 180 mosquitoes were successfully genotyped for all three different alleles (homozygote and heterozygote) of wild type, kdr west and kdr east (Fig. 3). The results are summarized in Table 4.
The total frequency of kdr west in all populations was 4.44% (5 heterozygotes and 3 homozygotes). The frequency of kdr in the three populations of Kunar, Laghman and Nangarhar were not significantly different. The frequency of the kdr alleles in pyrethroid bioassay survivors were, however, significantly different from the dead mosquitoes in all populations tested. Seven out of 8 kdr west mutated individuals were among bioassay survivors. Individuals homozygote for kdr west were detected only in the bioassay survivors (one in each population).
The frequency of kdr east was significantly higher than kdr west. The overall frequency of kdr east was 29%. The frequency of kdr east in the alive and dead specimens of the studied populations were 38% and 26.7%, respectively. The frequency of kdr east in Nangarhar, Laghman and Kunar populations were 35%, 32.7% and 22.9% respectively. Three of the four homozygote kdr east mutants were among the bioassay survivors.
Regardless of the type of the mutated alleles, the frequency of kdr in all populations were 33.33%. The frequency of kdr east and west combined in the populations were in the order of Nangarhar 39.6%, Laghman 37.9% and Kunar 27.9%. The allelic frequency of Leucin, Serine and Phenylalanine were 0.814, 0.156 and 0.03 respectively, all within the HWE (p = 0.14).
A general reduction in the susceptibility of An. stephensi from the eastern provinces of Afghanistan in recent years is probably in part due to the development of agriculture [51, 52] and increased use of different classes of insecticides that exerted a selection pressure on mosquitoes, especially where the larval breeding sites are within or close to the rice fields . The increase in resistance was most obvious to organophosphates and carbamates. Interestingly, the susceptibility to deltamethrin increased in the Laghman and Nangarhar populations, whereas to bendiocarb, it decreased in these populations. It might be due to the recent use of bendiocarb in rotation with deltamethrin in the region. Although resistance to bendiocarb has been reported in different mosquito species [53,54,55,56], development of bendiocarb resistance in An. stephensi populations in Afghanistan is now a cause for concern for the malaria control programme, who use this in IRS as an alternative to deltamethrin, in combination with LNs to manage insecticide resistance . Malathion and bendiocarb share the same mode of action and have the same target site resistance mechanism [57, 58]. The recent selection of bendiocarb resistance combined with the longstanding resistance to malathion [4, 9] suggests multiple resistance mechanisms in the Laghman population of An. stephensi.
Differential patterns of selection of insecticide resistance in populations of An. stephensi in different provinces of Kunar, Laghman and Nangarhar were observed . For example, the Kunar population remained susceptible to bendiocarb whereas the populations in Laghman and Nangarhar developed resistance to this insecticide between 2014 and 2017. Resistance levels to the pyrethroid insecticides permethrin and deltamethrin remained low, suggesting that resistance is conferred only by the relatively weak kdr mechanism.
Biochemical assays suggested that multiple metabolic mechanisms of resistance may have been selected in An. stephensi from Afghanistan. GST-based metabolism is often associated with DDT resistance [59,60,61]. GSTs activities were the highest in Laghman population compared with the other two field populations, and DDT resistance levels in Laghman population were the highest among all three field populations. Besides involving in DDT resistance, GSTs may be secondarily involved in pyrethroid insecticide resistance, so care must be taken in interpretation of the dynamics of this enzyme group as its metabolic role is not restricted to DDT .
Cytochrome P450s content in the Kunar, Nangarhar and Laghman populations of An. stephensi is increasing. Implications of this information for looking at the potential for moving to the new generations of PBO LNs should be considered . There are evidence that PBO LNs improved control of malaria transmission compared with standard long-lasting insecticidal nets where pyrethroid resistance is prevalent .
The mean inhibition rates and frequency of iAChE individuals in the field populations of An. stephensi from Afghanistan are significantly higher than that of the susceptible Beech strain. The iAChE should confer resistance to malathion and bendiocarb in An. stephensi [65, 66].
The importance of different enzyme groups in conferring insecticide resistance in different insects, especially mosquitoes, is overwhelming [12, 21, 23, 24, 27, 29, 55]. General esterases and cytochrome P450s are involved in pyrethroid insecticide resistance in An. stephensi from Dubai, Iran and India [12, 20, 21, 23, 67]. Esterases also confer resistance to OPs and to a lesser extent to pyrethroid insecticides [19, 68]. Involvement of iAChE in OPs and carbamate insecticides resistance is evident in many insect groups including An. stephensi from Iran, Afghanistan and India [12, 23, 27, 28]. A similar pattern of AChE insensitivity was seen in Anopheles albimanus in Mexico , in Turkish populations of the Anopheles maculipennis , and An. stephensi from Iran . The reduction in susceptibility of An. stephensi from Afghanistan to bendiocarb in recent years  is of concern, as it is the insecticide of choice for IRS and insecticide resistance management in the East of Afghanistan is reliant on this in combination with LN distribution . Close monitoring is now required to ensure that the situation does not deteriorate and impact on the ability to control malaria.
The current study suggests that both metabolic and target site resistance have increased between 2014 and 2017 in Afghanistan. Knock down resistance mechanism of pyrethroid resistance in An. stephensi was first determined in 2003 by Enayati et al.  in the DUB-S strain. At that time, only L1014F (later known as kdr west) mutation was observed in the species. After the discovery of L1014S mutation (known as kdr east) in An. gambiae in Kenya , reports of the presence of this trait in different mosquito species emerged [37, 38, 40, 71]. The presence of kdr east in An. stephensi from India was reported for the first time by Singh et al. . The development of kdr mechanism in An. stephensi from Afghanistan was first reported in 2016 , a report that triggered this bigger study in populations of An. stephensi from Kunar, Laghman and Nangarhar.
The overall frequency of kdr east in Afghan An. stephensi has increased by 8% between 2016 and 2017 . The kdr mutants were significantly more frequent in bioassay survivors compared to dead (38% vs 27%), as expected from this relatively weak pyrethroid resistance mechanism. The allelic frequency of kdr east was much lower than that reported in an Indian population of An. stephensi , however, the allelic frequency of kdr west was higher. Kdr west and east are prevalent in different Anopheles species worldwide [72, 73]. In An. gambiae in many West African countries, the frequency of kdr west has reached fixation (total homozygocity) [54, 74,75,76,77]. In other parts of Africa, co-occurrence of kdr west and kdr east with relatively high frequency is reported [34, 37, 71, 73, 78]. It has been argued that even with high kdr mutation rates, LNs may still be effective in protecting people against malaria [79,80,81,82,83,84,85]. However, a systematic review with meta-analysis revealed that the efficacy of the LNs diminishes with high intensity insecticide resistance [86, 87].
The development of different kdr traits may reflect the history of phenotypic resistance to different classes of insecticides. It was shown that kdr west confers stronger resistance to pyrethroid insecticides compared with kdr east, as the latter confers higher resistance to DDT and lesser to pyrethroids . The susceptibility of Laghman population of An. stephensi to DDT (45%) is lower compared to Kunar (the highest 70%) and Nangarhar populations (60%).
The proportion of arable land to the total area of the country is in the order of Nangarhar (12%) > Kunar (7%) > Laghman (4%) . The intensive use of organophosphorus and carbamate insecticides in agriculture, coupled with the historical and current use of these insecticides in public health has selected increasing levels of resistance to malathion and bendiocarb in Afghan An. stephensi between 2014 and 2017 . Larviciding with temephos and IRS using bendiocarb may exacerbate the levels of resistance to these classes of insecticides in An. stephensi.
The results of biochemical assays support the development of resistance to all classes of insecticides in An. stephensi in these provinces. A detailed resistance management plan should now be prepared to address potential issues going forward before resistance problem start to reduce the effectiveness of malaria control.
Anopheles stephensi populations from Kunar, Laghman and Nangarhar developed a range of resistance to different insecticides. Resistance to DDT, malathion, bendiocarb and pyrethroid insecticides is evident in different populations of the mosquito. Different enzyme groups are involved in the resistance to insecticides in An. stephensi from Kunar, Laghman and Nangarhar Provinces of Afghanistan. The contents of cytochrome P450s and the levels of activities of general esterases, glutathione S-transferases and the insensitivity of acetylcholinesterase are increasing in these populations. Based on the results, it can be concluded that the strength of metabolic resistance in An. stephensi from Laghman is slightly higher to multiple insecticides than the other two field populations. Knockdown resistance (both L1014F and L1014S) gene frequency is also on the rise. Moreover, homozygote individuals for either kdr traits have been detected for the first time in the field populations. These observations have huge technical as well as practical implications to malaria control programmes in Afghanistan. Based on these observations, the following general recommendations can be made: (i) vector control programmes need to be evidence-based and to be guided by routine monitoring and evaluation of vector control interventions including susceptibility to insecticides; (ii) malaria vector control is implemented based on careful stratification which includes, among other things, vector susceptibility to insecticides and their underlying resistance mechanisms; (iii) close collaboration of public health and agriculture sectors is required for effective management of insecticide resistance. In other words, for sound insecticide resistance management, implementation of strategies such as IVM in public health and IPM in agriculture is recommended. On the other hand, in planning and implementation of control measures for vector-borne diseases especially malaria, considerations should be made not only to the growing insecticide resistance status in the vectors, but also to its mechanisms in An. stephensi in the East of Afghanistan. Given the differential development of insecticide resistance and its underlying mechanisms in An. stephensi populations from different provinces of Afghanistan, it is recommended that the malaria control planning be province specific e.g. applying bendiocarb IRS in Kunar where An. stephensi is susceptible to this insecticide, while deploying PBO nets in Laghman and Nangarhar Provinces where cytochrome P450 is involved in pyrethroid insecticides resistance. LNs with alternative insecticides e.g. chlorfenapyr, and primiphos methyl long lasting IRS may also be answers to some of the challenges of insecticide resistance in An. stephensi in Afghanistan.
Availability of data and materials
All data generated or analysed during this study are included in this published article.
global plan for insecticide resistance management
indoor residual spraying
long-lasting insecticidal nets
Ministry of Public Health
National Malaria and Leishmaniasis Control Programme
polymerase chain reaction
sodium dodecyl sulfate
voltage gated sodium channel
World Health Organization
World Health Organization Pesticide Evaluation Scheme
WHO. World malaria report. Geneva: World Health Organization; 2018.
Afghanistan Ministry of Public Health. National Malaria Strategic Plan 2013–2017. Islamic Republic of Afghanistan. Kabul: Afghanistan Ministry of Public Health; 2013.
Rowland M, Mohammed N, Rehman H, Hewitt S, Mendis C, Ahmad M, et al. Anopheline vectors and malaria transmission in eastern Afghanistan. Trans R Soc Trop Med Hyg. 2002;96:620–6.
Safi N, Hameed H, Sediqi W, Himmat E. National Malaria and Leishmaniasis Control Programme (NMLCP) Annual Report, 2008. Afghan Annu Malar J. 2009;1:8–14.
Singh OP, Dykes CL, Lather M, Agrawal OP, Adak T. Knockdown resistance (kdr)-like mutations in the voltage-gated sodium channel of a malaria vector Anopheles stephensi and PCR assays for their detection. Malar J. 2011;10:59.
Safi NHZ. Status of insecticide susceptibility in Laghman, Kunar and Nangarhar Provinces, Afghanistan in 2017. Kabul: Ministry of Public Health, Afghanistan; 2017.
Ministry of Public Health. National malaria strategic plan 2013–1017. National Malaria and Leishmaniasis Control Programme; 2013. http://moph.gov.af/Content/Media/Documents/NationalStrategicPlanforMalariaControlandElimination184201510534818553325325.doc. Accessed 26 Dec 2018.
Ahmad M. Evaluating insecticide resistance status in wild-caught Anopheles mosquitoes in Afghanistan. Kabul: HealthNet-TPO; 2015.
Barwa C. Status of Insecticide sSusceptibility in Afghanistan. (General Directorate of Preventive Medicine and Primary Health Care CDD, National Malaria and Leishmaniasis Control Programme, Afghanistan: Ministry of Public Health, Islamic Republic of Afghanistan; 2011. https://moph.gov.af/Content/Media/Documents/NationalStrategicPlanforMalariaControlandElimination184201510534818553325325.doc. Accessed 26 Dec 2018.
Singh OP, Dykes CL, Adak T. Kdr-like mutations in the voltage gated sodium channel of a malaria vector Anopheles stephensi and development of PCR-based assays for their detection. Malar J. 2010;9(Suppl 2):P67.
Soltani A, Vatandoost H, Oshadhi MA, Enayati AA, Raeisi A, Eshraghian MR, et al. Baseline susceptibility of different geographical strains of Anopheles stephensi (Diptera: Culicidae) to temephos in malarious areas of Iran. J Arthropod-Borne Dis. 2013;7:56–65.
Soltani A, Vatandoost H, Oshaghi MA, Ravasan NM, Enayati AA, Asgarian F. Resistance mechanisms of Anopheles stephensi (Diptera: Culicidae) to temephos. J Arthropod-Borne Dis. 2014;9:71–83.
Soltani A, Vatandoost H, Oshaghi MA, Ravasan NM, Enayati AA, Asgarian F. Resistance mechanisms of Anopheles stephensi (Diptera: Culicidae) to temephos. J Arthropod-borne Dis. 2015;9:71–83.
Dykes CL, Das MK, Eapen A, Batra CP, Ghosh SK, Vijayan VA, et al. Knockdown resistance (kdr) mutations in Indian Anopheles stephensi (Diptera: Culicidae) populations. J Med Ent. 2016;53:315–20.
Gorouhi MA, Vatandoost H, Oshaghi MA, Raeisi A, Enayati AA, Mirhendi H, et al. Current susceptibility status of Anopheles stephensi (Diptera: Culicidae) to different imagicides in a malarious area, Southeastern of Iran. J Arthropod-borne Dis. 2016;10:493–500.
Mogaddam MY, Haghi FM, Fazeli-Dinan M, Hoseini N, Enayati AA. A review of insecticide resistance in malaria vectors of Iran. J Mazandaran Univ Med Sci. 2016;25:394–411.
Enayati AA, Ladonni H. Biochemical assays baseline data of permethrin resistance in Anopheles stephensi (Diptera, Culicidae) from Iran. Pak J Biol Sci. 2006;9:1265–70.
Enayati AA, Ranson H, Hemingway J. Insect glutathione transferases and insecticide resistance. Insect Mol Biol. 2005;14:3–8.
Enayati AA, Vatandoost H, Ladonni H, Townson H, Hemingway J. Molecular evidence for a kdr-like pyrethroid resistance mechanism in the malaria vector mosquito Anopheles stephensi. Med Vet Entomol. 2003;17:138–44.
Ganesh K, Vijayan V, Urmila J, Gopalan N, Prakash S. Role of esterases and monooxygenase in the deltamethrin resistance in Anopheles stephensi Giles (1908), at Mysore. Indian J Exp Biol. 2002;40:583–8.
Ganesh KN, Urmila J, Vijayan VA. Pyrethroid susceptibility & enzyme activity in two malaria vectors, Anopheles stephensi (Liston) & Anopheles culicifacies (Giles) from Mysore, India. Indian J Med Res. 2003;117:30–8.
Gayathri V, Balakrishna Murthy P. Reduced susceptibility to deltamethrin and kdr mutation in Anopheles stephensi Liston, a malaria vector in India. J Am Mosq Control Assoc. 2006;22:678–88.
Hariprasad TP, Shetty NJ. Biochemical basis of alphamethrin resistance in different life stages of Anopheles stephensi strains of Bangalore, India. Pest Manag Sci. 2016;72:1689–701.
Hemingway J. The molecular basis of two contrasting metabolic mechanisms of insecticide resistance. Insect Biochem Mol Biol. 2000;30:1009–15.
Hemingway J, Ranson H. Insecticide resistance in insect vectors of human disease. Annu Rev Entomol. 2000;45:369–89.
Hodjati MH, Curtis CF. Genetics of resistance to permethrin in Anopheles stephensi. J Vector Borne Dis. 2003;40:60–3.
Prasad KM, Raghavendra K, Verma V, Velamuri PS, Pande V. Esterases are responsible for malathion resistance in Anopheles stephensi: A proof using biochemical and insecticide inhibition studies. J Vector Borne Dis. 2017;54:226–32.
Safi NHZ, Ahmadi AA, Nahzat S, Ziapour SP, Nikookar SH, Fazeli-Dinan M, et al. Evidence of metabolic mechanisms playing a role in multiple insecticides resistance in Anopheles stephensi populations from Afghanistan. Malar J. 2017;16:100.
Sanil D, Shetty V, Shetty N. Differential expression of glutathione s-transferase enzyme in different life stages of various insecticide-resistant strains of Anopheles stephensi: a malaria vector. J Vector Borne Dis. 2014;51:97–105.
Djègbè 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.
Aïzoun N, Aïkpon R, Akogbéto M. Evidence of increasing L1014F kdr mutation frequency in Anopheles gambiae s.l. pyrethroid resistant following a nationwide distribution of LLINs by the Beninese National Malaria Control Programme. Asian Pac J Trop Biomed. 2014;4:239–43.
Ibrahim SS, Manu YA, Tukur Z, Irving H, Wondji CS. High frequency of kdr L1014F is associated with pyrethroid resistance in Anopheles coluzzii in Sudan savannah of northern Nigeria. BMC Infect Dis. 2014;14:441.
Camara S, Koffi AA, Ahoua Alou LP, Koffi K, Kabran JPK, Koné A, et al. Mapping insecticide resistance in Anopheles gambiae (s.l.) from Côte d’Ivoire. Parasit Vectors. 2018;11:19.
Thiaw O, Doucouré S, Sougoufara S, Bouganali C, Konaté L, Diagne N, et al. Investigating insecticide resistance and knock-down resistance (kdr) mutation in Dielmo, Senegal, an area under long lasting insecticidal-treated nets universal coverage for 10 years. Malar J. 2018;17:123.
Zoh DD, Ahoua Alou LP, Toure M, Pennetier C, Camara S, Traore DF, et al. The current insecticide resistance status of Anopheles gambiae (s.l.) (Culicidae) in rural and urban areas of Bouaké. Côte d’Ivoire. Parasit Vectors. 2018;11:118.
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.
Pinto J, Lynd A, Elissa N, Donnelly MJ, Costa C, Gentile G, et al. 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 Ent. 2006;20:27–32.
Verhaeghen K, Van Bortel W, Roelants P, Backeljau T, Coosemans M. Detection of the East and West 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.
Namountougou M, Diabaté A, Etang J, Bass C, Sawadogo SP, Gnankinié O, et al. First report of the L1014S kdr mutation in wild populations of Anopheles gambiae M and S molecular forms in Burkina Faso (West Africa). Acta Trop. 2013;125:123–7.
Chouaïbou M, Kouadio FB, Tia E, Djogbenou L. First report of the East African kdr mutation in an Anopheles gambiae mosquito in Côte d’Ivoire. Wellcome Open Res. 2017;2:8.
Singh OP, Bali P, Hemingway J, Subbarao SK, Dash AP, Adak T. PCR-based methods for the detection of L1014 kdr mutation in Anopheles culicifacies sensu lato. Malar J. 2009;8:154.
Singh OP, Dykes CL, Das MK, Pradhan S, Bhatt RM, Agrawal OP, et al. Presence of two alternative kdr-like mutations, L1014F and L1014S, and a novel mutation, V1010L, in the voltage gated Na+ channel of Anopheles culicifacies from Orissa, India. Malar J. 2010;9:146.
Ahmad M, Buhler C, Pignatelli P, Ranson H, Nahzat SM, Naseem M, et al. Status of insecticide resistance in high-risk malaria provinces in Afghanistan. Malar J. 2016;15:98.
WHO. Malaria surveillance, monitoring & evaluation: a reference manual. Geneva: World Health Organization; 2018.
WHO. Global malaria programme. Test procedures for insecticide resistance monitoring in malaria vector mosquitoes. 2nd ed. Geneva: World Health Organization; 2016.
WHO. Global plan for insecticide resistance management in malaria vectors (GPIRM). Geneva: World Health Organization; 2012.
Glick J. Illustrated key to the female Anopheles of Southwestern Asia and Egypt (Diptera: Culicidae). Mosq Syst. 1992;4:125–53.
Hemingway J. Field and laboratory manual for the mechanistic detection of insecticide resistance in insects. WHO-CTD/MAL/986. Geneva: World Health Organization; 1998. p. 35.
Livak KJ. Organization and mapping of a sequence on the Drosophila melanogaster X and Y chromosomes that is transcribed during spermatogenesis. Genetics. 1984;107:611–34.
Rodriguez S, Gaunt TR, Day IN. Hardy-Weinberg equilibrium testing of biological ascertainment for mendelian randomization studies. Am J Epidemiol. 2009;6:359.
Harris W. Agricultural developement for afghansitan, agricultural review of Kunar Province. Afghanistan; 2011. http://www.fresnostate.edu/jcast/ifa/documents/pdf/adapt/KUNAR.pdf. Accessed 26 Dec 2018.
Afghanistan Research and Evaluation Unit. Provincial profiles: Regional Rural Economic Regeneration Assessment and Strategies (RRERS) study/National Area Based Development Programme, Ministry of Rural Rehabilitation and Development (MRRD), GRM International; 2006.
Aïkpon R, Agossa F, Ossè R, Oussou O, Aïzoun N, Oké-Agbo F, et al. Bendiocarb resistance in Anopheles gambiae s.l. populations from Atacora department in Benin, West Africa: a threat for malaria vector control. Parasit Vectors. 2013;6:192.
Protopopoff N, Matowo J, Malima R, Kavishe R, Kaaya R, Wright A, et al. High level of resistance in the mosquito Anopheles gambiae to pyrethroid insecticides and reduced susceptibility to bendiocarb in north-western Tanzania. Malar J. 2013;12:149.
Edi CV, Djogbenou L, Jenkins AM, Regna K, Muskavitch MA, Poupardin R, et al. CYP6 P450 enzymes and ACE-1 duplication produce extreme and multiple insecticide resistance in the malaria mosquito Anopheles gambiae. PLoS Genet. 2014;10:e1004236.
Raghavendra K, Barik TK, Sharma P, Bhatt RM, Srivastava HC, Sreehari U, et al. Chlorfenapyr: a new insecticide with novel mode of action can control pyrethroid resistant malaria vectors. Malar J. 2011;10:16.
Bigoga JD, Ndangoh DN, Awono-Ambene PH, Patchoke S, Fondjo E, Leke RG. Pyrethroid resistance in Anopheles gambiae from the rubber cultivated area of Niete, South Region of Cameroon. Acta Trop. 2012;124:210–4.
Fishel FM. IRAC’s Insecticide mode of Action Classification. https://edis.ifas.ufl.edu/pdffiles/PI/PI12100.pdf. Accessed 26 Dec 2018.
Djègbè I, Agossa FR, Jones CM, Poupardin R, Cornelie S, Akogbéto M, et al. Molecular characterization of DDT resistance in Anopheles gambiae from Benin. Parasit Vectors. 2014;7:409.
Etang J, Manga L, Toto JC, Guillet P, Fondjo E, Chandre F. Spectrum of metabolic-based to DDT and pyrethroids in Anopheles gambiae s.l. populations from Cameroon. J Vector Ecol. 2007;32:123–33.
Fossog Tene B, Poupardin R, Costantini C, Awono-Ambene P, Wondji CS, Ranson H, et al. Resistance to DDT in an urban setting: common mechanisms implicated in both M and S forms of Anopheles gambiae in the city of Yaoundé Cameroon. PLoS ONE. 2013;8:e61408.
Vontas JG, Small GJ, Hemingway J. Glutathione S-transferases as antioxidant defence agents confer pyrethroid resistance in Nilaparvata lugens. Biochem J. 2001;357:65–72.
Ketoh GK, Ahadji-Dabla KM, Chabi J, Amoudji AD, Apetogbo GY, Awokou F, et al. Efficacy of two PBO long lasting insecticidal nets against natural populations of Anopheles gambiae s.l. in experimental huts, Kolokopé, Togo. PLoS ONE. 2018;13:e0192492.
Protopopoff N, Mosha JF, Lukole E, Charlwood JD, Wright A, Mwalimu CD, et al. Effectiveness of a long-lasting piperonyl butoxide-treated insecticidal net and indoor residual spray interventions, separately and together, against malaria transmitted by pyrethroid-resistant mosquitoes: a cluster, randomised controlled, two-by-two factorial design trial. Lancet. 2018;391:1577–88.
Cisse MBM, Keita C, Dicko A, Dengela D, Coleman J, Lucas B, et al. Characterizing the insecticide resistance of Anopheles gambiae in Mali. Malar J. 2015;14:327.
Assogba BS, Djogbénou LS, Saizonou J, Milesi P, Djossou L, Djegbe I, et al. Phenotypic effects of concomitant insensitive acetylcholinesterase (ace-1R) and knockdown resistance (kdrR) in Anopheles gambiae: a hindrance for insecticide resistance management for malaria vector control. Parasit Vectors. 2014;7:548.
Nikpour F, Vatandoost H, Hanafi-Bojd AA, Raeisi A, Ranjbar M, Enayati AA, et al. Evaluation of deltamethrin in combination of piperonyl butoxide (PBO) against pyrethroid resistant malaria vector, Anopheles stephensi in IRS implementation: an experimental semi-filed trial in Iran. J Arthropod-borne Dis. 2017;11:469–81.
Montella IR, Schama R, Valle D. The classification of esterases: an important gene family involved in insecticide resistance-A review. Mem Inst Oswaldo Cruz. 2012;107:437–49.
Penilla RP, Rodrigues AD, Hemingway J, Torres JL, Arredondo-Jimenez JI, Rodriguez MH. Resistance management strategies in malaria vector mosquito control. Baseline data for a large-scale field trial against Anopheles albimanus in Mexico. Med Vet Ent. 1998;12:217–33.
Akiner MM. Malathion and propoxur resistance in Turkish populations of the Anopheles maculipennis Meigen (Diptera: Culicidae) and relation to the insensitive acetylcholinesterase. Turkiye Parazitol Derg. 2014;38:111.
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 resistance in Anopheles gambiae sensu lato population of Tanzania. Trop Med Int Health. 2014;19:331–41.
Protopopoff N, Verhaeghen K, Van Bortel W, Roelants P, Marcotty T, Baza D, et al. A significant increase in kdr in Anopheles gambiae is associated with an intensive vector control intervention in Burundi highlands. Trop Med Int Health. 2008;13:1479–87.
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:364.
Nardini L, Hunt RH, Dahan-Moss YL, Christie N, Christian RN, Coetzee M, et al. Malaria vectors in the Democratic Republic of the Congo: the mechanisms that confer insecticide resistance in Anopheles gambiae and Anopheles funestus. Malar J. 2017;16:448.
Chandre F, Darriet F, Duchon S, Finot L, Manguin S, Carnevale P, et al. Modifications of pyrethroid effects associated with kdr mutation in Anophelos gambiae. Med Vet Ent. 2000;14:81–8.
Mathias DK, Ochomo E, Atieli F, Ombok M, Nabie Bayoh M, 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.
Mawejje HD, Wilding CS, Rippon EJ, Hughes A, Weetman D, Donnelly MJ. Insecticide resistance monitoring of field-collected Anopheles gambiae s.l. populations from Jinja, eastern Uganda, identifies high levels of pyrethroid resistance. Med Vet Ent. 2013;27:276–83.
Mint Mohamed Lemine A, Ould Lemrabott MA, Niang EHA, Basco LK, Bogreau H, Faye O, et al. Pyrethroid resistance in the major malaria vector Anopheles arabiensis in Nouakchott. Mauritania. Parasit Vectors. 2018;11:344.
Darriet F, Guillet P, Guessan RNN, Doannio JMC, Koffi AA, Konan LY, et al. The impact of permethrin and deltamethrin resistance in Anopheles gambiae s.s. on the efficacy of insecticide-treated mosquito nets. Geneva: World Health Organization; 1999; https://apps.who.int/iris/bitstream/handle/10665/66463/WHO_MAL_99.1088.pdf?sequence=1&isAllowed=y. Accessed 26 Dec 2018.
Darriet F, Guillet P, N’Guessan R, Doannio JM, Koffi A, Konan LY, et al. Impact of resistance of Anopheles gambiae ss to permethrin and deltamethrin on the efficacy of impregnated mosquito nets. Med Trop (Mars). 1998;58:349–54 (in French).
Darriet F, N’Guessan R, Koffi AA, Konan L, Doannio JM, Chandre F, et al. Impact of pyrethrin resistance on the efficacity of impregnated mosquito nets in the prevention of malaria: results of tests in experimental cases with deltamethrin SC. Bull Soc Pathol Exot. 2000;93:131–4 (in French).
Darriet F, Robert V, Thovein N, Carnevalle P. Evaluation of the efficacy of permethrin-impregnated intact and perforated mosquito nets against vectors of malaria. Geneva: World Health Orgnization, WHO/VBC/84:899.
Henry MC, Doannio JM, Darriet F, Nzeyimana I, Carnevale P. Efficacy of permethrin-impregnated Olyset Net mosquito nets in a zone with pyrethroid resistance vectors. II. Parasitic and clinical evaluation. Med Trop (Mars). 1999;59:355–7 (in French).
N’Guessan R, Corbel V, Akogbéto M, Rowland M. Reduced efficacy of insecticide-treated nets and indoor residual spraying for malaria control in pyrethroid resistance area, Benin. Emerg Infect Dis. 2007;13:199–206.
N’Guessan R, Darriet F, Doannio JM, Chandre F, Carnevale P. Olyset Net efficacy against pyrethroid-resistant Anopheles gambiae and Culex quinquefasciatus after 3 years’ field use in C te d’Ivoire. Med Vet Entomol. 2001;15:97–104.
Strode C, Donegan S, Garner P, Enayati AA, Hemingway J. The Impact of pyrethroid resistance on the efficacy of insecticide-treated bed nets against African anopheline mosquitoes: systematic review and meta-analysis. PLoS Med. 2014;11:e1001619.
Enayati AA, Hemingway J. Pyrethroid insecticide resistance and treated bednets efficacy in malaria control. Pest Biochem Physiol. 2006;84:116–26.
Reimer L, Fondjo E, Patchoké S, Diallo B, Lee Y, Ng A, et al. Relationship between kdr mutation and resistance to pyrethroid and DDT insecticides in natural populations of Anopheles gambiae. J Med Entomol. 2008;45:260–6.
This project was supported financially by the WHO Regional Office for the Eastern Mediterranean and WHO country office Afghanistan. Thanks also to the personnel of the biochemistry of pesticide lab and Cellular and Molecular Biology Research Center at School of Public Health and Faculty of Medicine, Mazandaran University of Medical Sciences, Iran for their assistance in biochemical and molecular analyses of the specimens.
This project was supported financially by WHO Regional Office for the Eastern Mediterranean and technically by WHO country office Afghanistan and Iran.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.