Insecticide resistance status of Aedes aegypti in Bangladesh
Parasites & Vectors volume 13, Article number: 622 (2020)
Arboviral diseases, including dengue and chikungunya, are major public health concerns in Bangladesh where there have been unprecedented levels of transmission reported in recent years. The primary approach to control these diseases is to control the vector Aedes aegypti using pyrethroid insecticides. Although chemical control has long been practiced, no comprehensive analysis of Ae. aegypti susceptibility to insecticides has been conducted to date. The aim of this study was to determine the insecticide resistance status of Ae. aegypti in Bangladesh and investigate the role of detoxification enzymes and altered target site sensitivity as resistance mechanisms.
Eggs of Aedes mosquitoes were collected using ovitraps from five districts across Bangladesh and in eight neighborhoods of the capital city Dhaka, from August to November 2017. CDC bottle bioassays were conducted for permethrin, deltamethrin, malathion, and bendiocarb using 3- to 5-day-old F0–F2 non-blood-fed female mosquitoes. Biochemical assays were conducted to detect metabolic resistance mechanisms, and real-time PCR was performed to determine the frequencies of the knockdown resistance (kdr) mutations Gly1016, Cys1534, and Leu410.
High levels of resistance to permethrin were detected in all Ae. aegypti populations, with mortality ranging from 0 to 14.8% at the diagnostic dose. Substantial resistance continued to be detected against higher (2×) doses of permethrin (5.1–44.4% mortality). Susceptibility to deltamethrin and malathion varied between populations while complete susceptibility to bendiocarb was observed in all populations. Significantly higher levels of esterase and oxidase activity were detected in most of the test populations as compared to the susceptible reference Rockefeller strain. A significant association was detected between permethrin resistance and the presence of Gly1016 and Cys1534 homozygotes. The frequency of kdr (knockdown resistance) alleles varied across the Dhaka Aedes populations. Leu410 was not detected in any of the tested populations.
The detection of widespread pyrethroid resistance and multiple resistance mechanisms highlights the urgency for implementing alternate Ae. aegypti control strategies. In addition, implementing routine monitoring of insecticide resistance in Ae. aegypti in Bangladesh will lead to a greater understanding of susceptibility trends over space and time, thereby enabling the development of improved control strategies.
Aedes (Stegomyia) aegypti (Linnaeus, 1762) is an important vector of arboviral diseases, principally dengue, chikungunya, and Zika. These increasingly common arboviral infections cause severe febrile illness and short- to long-term physical or cognitive impairments and even death. Dengue is the most prevalent and rapidly spreading arboviral disease worldwide, with an estimated 390 million annual infections and 3.9 billion people at risk . Chikungunya is also increasingly prevalent, and the prolonged pain and rheumatism resulting from infection can result in long-term physical problems and impaired daily life [2, 3]. Zika recently caused a major global pandemic in 2015–2016, leading to congenital malformations, Guillain–Barre syndrome, and other severe neurological complications .
The burden of arboviral diseases in Bangladesh is not well documented. The first major outbreak of dengue took place during the monsoon of 2000 and caused 5521 officially reported cases with 93 deaths . Since then, thousands of infections are reported each year although these numbers represent a fraction of the actual burden since only patients admitted to some selected hospitals are officially reported . Recent estimates suggest that 40 million people have been infected nationally, with an average of 2.4 million infections annually. Cases are mostly concentrated in the capital city Dhaka, where the seropositivity ranges from 36 to 85% . In 2019, Bangladesh experienced its largest outbreak, with 101,354 confirmed cases and 164 deaths . Since 2008, sporadic infections with chikungunya virus have been reported across Bangladesh, with the largest outbreak occurring in 2017 during which hundreds of thousands of inhabitants of Dhaka were infected . Zika virus transmission has not been widely reported in Bangladesh, with only a single confirmed case in 2016 in a 67-year old man from Chittagong who had not traveled outside of Bangladesh. Although a few additional Zika virus infections have been detected by antibody tests, there is no further evidence of Zika in Bangladesh [10, 11].
Aedes aegypti is highly abundant throughout Bangladesh, especially in Dhaka . In 2018, the Breteau index (BI; the number of Aedes-positive containers per 100 houses inspected) was > 100 in some parts of Dhaka . Recent studies in Dhaka have confirmed that plastic containers (plastic drums, buckets, plastic bags, bottles, and disposable cups) and discarded vehicle and construction materials (tires, battery shells, and cement mixers) are key containers for Aedes production. These are typical of the domestic and industrial detritus that facilitate the proliferation of Ae. aegypti across the globe. High Aedes abundance in Dhaka is also strongly associated with favorable climatic factors, including rainfall, temperature, and humidity .
In the absence of effective therapeutic drugs and vaccines, Ae. aegypti control is presently the only approach for preventing and controlling the transmission of most Aedes-borne arboviruses. Aedes aegypti control strategies rely heavily on the application of a limited number of chemical insecticides approved for public health use, principally pyrethroids, organochlorines, organophosphates, and carbamates . Of these, pyrethroid insecticides, such as deltamethrin, cypermethrin, and permethrin, are commonly used because of their low toxicity to mammals and their high efficacy against vectors. However, resistance to many insecticides has emerged in Ae. aegypti across the globe, representing a serious threat to control programs [15,16,17,18,19].
Resistance to insecticides is a dynamic evolutionary process that is driven by insecticide selection pressures . Resistance can be caused by physiological changes, including (i) changes to the mosquito cuticle so insecticides cannot penetrate, (ii) increased activity of insecticide detoxification enzymes, and/or (iii) structural modifications at the target site of the insecticide; it can also result from behavioral adaptations, such as insecticide avoidance .
Alterations in the target sites as a result of resistance to pyrethroids and dichlorodiphenyltrichloroethane (DDT) are often caused by mutations in the voltage-gated sodium channel (VGSC) transmembrane protein and are broadly referred to as ‘knockdown resistance’ (kdr) mutations. There are several point mutations on the VGSC gene known to confer kdr-type insecticide resistance in Ae. aegypti, most notably at positions 410, 989, 1016, and 1534 [18, 22, 23]. Increased enzyme activity resulting in metabolic resistance typically involves any of the three main groups of detoxification enzymes: carboxylesterases, mixed-function oxidases (MFOs), and glutathione S-transferases (GSTs) . Understanding the mechanisms of resistance and their specificity among insecticides is important to devising strategies to mitigate and manage insecticide resistance when it is detected.
Although there has been a recognized increase in Aedes-borne arboviruses in Bangladesh over the last 20 years, little or no organized use of insecticides against Ae. aegypti has occurred. Regular control activities are mostly carried out only in Dhaka, targeting the nuisance biting mosquitoes Culex quinquefasciatus and Aedes by thermal fogging with a combination of pyrethroid insecticides, including permethrin, prallethrin, and tetramethrin/bioallethrin. The increasing number of reported cases of Aedes-borne viral diseases indicate that the insecticides being used are having little impact. Development of resistance against commonly used insecticides in local Aedes populations may contribute to the failure of the vector control strategy. Occasional source reduction is also carried out by community engagement by both government and private initiatives. However, gaining access to all premises and achieving sufficient coverage of myriad oviposition sites in densely populated cities like Dhaka is a huge challenge . There are also structural challenges to control activities related to management, evaluation, and budget [26, 27].
The insecticide resistance status of Ae. aegypti has not previously been comprehensively assessed in Bangladesh. The purpose of this study was to assess the insecticide resistance status and resistance mechanisms of key Ae. aegypti populations in Bangladesh to better inform future insecticide choices for vector control.
Mosquitoes were collected from five districts throughout Bangladesh. Of these, the capital city, Dhaka, and the port city, Chittagong are high-transmission settings and the city of Rajshahi and district of Chapai Nawabganj are low-transmission settings [9, 28,29,30,31,32] (Fig. 1). The other district, Bandarban, was selected as it is a malaria-endemic region, and as preventative measures deltamethrin-impregnated long-lasting insecticidal nets (LLINs) are regularly distributed among the population and there is seasonal sporadic indoor residual spraying (IRS) [33, 34]. Since the majority of Aedes-borne arboviral infections are reported from Dhaka, eight areas within Dhaka City were selected for sampling .
Collection of Aedes eggs
Eggs were collected using ovitraps baited with a grass infusion. Ovitraps consisted of black 2-l containers made of plastic and an oviposition substrate of brown seed germination paper. The ovitraps were filled with 50 ml of 2- to 3-day-old grass infusion and 1200 ml of tap or rainwater. After verbal consent was obtained from household owners, the ovitraps were placed primarily indoors, including the main living area (under beds), behind refrigerators, under stairways, in garages, and on balconies. When these sites were not suitable, ovitraps were set in the yards under sheds close to the house. Within Dhaka City, the number of ovitraps varied from 50 to 70 per location, whereas in the areas outside of Dhaka (non-Dhaka), approximately 100 ovitraps were set in each location. For all non-Dhaka districts except Chittagong, eggs were collected from one urban and one rural location. All eggs were collected in 2017 during the traditional peak dengue transmission months of August to November.
Ovitraps were collected after 6 days in situ. Upon collection, the germination papers were dried and sent to the insectary at the Animal Research Facilities, International Centre for Diarrhoeal Disease Research, Bangladesh (icddr,b) in Dhaka. Due to unexpectedly long time required to prepare the rearing facility, mosquito rearing was delayed until December 2017. Given this delay, hatching rates were low for several locations, and eggs from adjacent locations were sometimes merged into a single population (Table 1). Mosquitoes were reared at a constant temperature (26–28 °C) and humidity (70–80%). When possible, mosquitoes were reared to the F2 generation to obtain sufficient numbers for a wide range of susceptibility tests. Artificial blood-feeding was provided using the methods described by Costa-da-Silva et al. . Adult mosquitoes were provided with 10% sucrose solution. In addition to the field populations, the ‘Rockefeller’ (ROCK) Ae. aegypti reference strain, which is insecticide susceptible, was obtained from the U.S. Centers for Disease Control and Prevention (CDC, Atlanta, GA, USA) and reared as a susceptible control in the bioassays.
Insecticide susceptibility testing
Susceptibility tests were conducted following the CDC bottle bioassay protocol  using 3- to 5-day-old, non-blood-fed female mosquitoes. Four insecticides belonging to three major chemical classes were tested for each population when sufficient mosquitoes were available: (i) the pyrethroids permethrin and deltamethrin; (ii) the organophosphate malathion; (iii) the carbamate bendiocarb. Mosquitoes were exposed to diagnostic dose of each insecticide, which was 15 µg/bottle for permethrin, 10 µg/bottle for deltamethrin, 50 µg/bottle for malathion, and 12.5 μg/bottle for bendiocarb. When resistance was detected and sufficient mosquitoes were available, resistance intensity assays were also conducted by exposing mosquitoes to insecticide levels that were two- and fivefold (2× and 5×, respectively) higher than the respective diagnostic dose. All bioassays comprised > 100 mosquitoes per insecticide per population across four test bottles and 15–25 mosquitoes in an untreated control bottle. Susceptibility status was recorded after 0, 15, and 30 min of insecticide exposure; at each time point, mosquitoes unable to stand were considered to be dead . Mortality data were interpreted according to World Health Organization (WHO) recommendations, with < 90% mortality in a population corresponding to resistance .
Biochemical assays were performed to detect potential metabolic mechanisms of resistance through the altered activity of detoxifying enzymes. From each population, thirty 1- to 2-day-old female mosquitoes were tested for non-specific β-esterases (β-EST), MFOs, acetylcholine esterase (AChE), and insensitive acetylcholine esterase (IAChE), with a protein assay conducted for each mosquito to control for differences in body size. All mosquitoes were freeze killed and kept at − 20 °C until analysis. Briefly, mosquitoes were individually homogenized in 100 µl of potassium phosphate buffer followed by dilution to 2 ml with additional buffer. For all tests, mosquito homogenates were run in triplicate on 96-well round-bottom microplates (Corning Inc., Corning, NY, USA). Homogenates of the ROCK insecticide-susceptible Ae. aegypti reference strain were used as a comparator.
For the β-EST assay, 100 µl of mosquito homogenate was added to each well followed by 100 µl β-naphthyl acetate. The plate was then incubated at room temperature for 20 min, following which 100 µl Fast Blue was added in each well and the plate was further incubated at room temperature for 4 min at which time absorbance was read on a spectrophotometer (BioTek Instruments Inc., Winooski, VT, USA) using a 540-nm filter.
For the MFO assay, 100 µl of mosquito homogenate was added to each well followed by 200 µl of 3,3,5,5-tetramethylbenzidine (TMBZ) and 25 µl 3% hydrogen peroxide. The plate was incubated for 10 min and absorbance was read on a spectrophotometer using a 620-nm filter.
For the AChE assay, 100 µl of mosquito homogenate was added to each well followed by 100 µl of acetylthiocholine iodide (ATCH) and 100 µl dithio-bis-2-nitrobenzoic acid (DTNB). Absorbance was read immediately (T0) using a 414-nm filter, and a second reading was taken after 10 min (T10). The absorbance at T0 was subtracted from T10 and used as the value for data analysis.
The IAChE assay was similar to the AChE assay, with the addition of propoxur to the ATCH to quantify the extent to which propoxur inhibited the reaction.
The total protein content of each mosquito was measured by adding 20 µl of the homogenate to a well together with 80 µl of potassium phosphate and 200 µl of protein dye. Absorbance was read immediately using a 620-nm filter. Due to unforeseen circumstances, we were unable to run a standard curve with the same reagents employed in the total protein assay of the mosquitoes, precluding precise quantification of protein content.
DNA extraction was carried out using the REDExtract-N-Amp™ tissue kit (Merck KGaA, Darmstadt, Germany) according to the manufacturer’s protocol. Briefly, individual female mosquitoes were placed in 1.5-ml microcentrifuge tubes and mixed with 100 µl of extraction solution and 25 µl of tissue preparation solution. The tubes were then incubated at room temperature for 10 min followed by further incubation for 3 min at 95 °C. Then, 100 µl of neutralization solution B was added to the sample and the sample mixed by vortexing.
Detection of kdr alleles (Gly1016, Cys1534, and Leu410)
To understand the correlation between phenotypic resistance and the presence of the kdr alleles Gly1016, Cys1534C, and Leu410, phenotyped mosquitoes exposed to permethrin and deltamethrin in the bioassays were assayed by real-time PCR. An additional 30 non-phenotyped mosquitoes from each of the six Dhaka populations were analyzed to estimate allele frequencies at the population level.
The Gly1016 PCR was performed following the protocol described by Saavedra-Rodriguez et al. . Each reaction contained 4.5 µl of iQ-SYBR Green Supermix (Bio-Rad Laboratories Inc., Hercules, CA, USA), 0.45 µl of each primer, one common Gly forward primer (5ʹ-ACC GAC AAA TTG TTT CCC-3ʹ), one reverse primer for either Val (5ʹ-GCG GGC AGC AAG GCT AAG AAA AGG TTA ATT A-3ʹ) or Gly (5ʹ-GCG GGC AGG GCG GGG GCG GGG CCA GCA AGG CTA AGA AAA GGT TAA CTC-3ʹ), 1 µl of template DNA, and ddH2O for a final reaction volume of 9 µl. Thermal cycling conditions were: 95 °C for 3 min; 40 cycles of 95 °C for 10 s, 58 °C for 10 s, 72 °C for 30 s; 95 °C for 10 s; and a ramp from 65 °C to 95 °C at a rate of 0.2 °C/10 s for melting curve analysis.
The Cys1534 PCR was based on the protocol described by Yanola et al. . Each reaction contained 4.5 µl of iQ-SYBR Green Supermix (Bio-Rad Laboratories Inc.), 0.45 µl Cys forward primer (5ʹ-GCG GGC AGG GCG GCG GGG GCG GGG CCT CTA CTT TGT GTT CTT CAT CAT GTG-3ʹ), and 0.45 µl each of Phe forward primer (5′-GCG GGC TCT ACT TTG TGT TCT TCA TCA TAT T-3′) and a common reverse primer (5′-TCT GCT CGT TGA AGT TGT CGA T-3′), 1 µl of template DNA and ddH2O for a final reaction volume of 9 µl. Thermal cycling conditions were: 95 °C for 3 min; 40 cycles of 95 °C for 10 s, 57 °C for 10 s, 72 °C for 30 s; 95 °C for 10 s; and a ramp from 65 °C to 95 °C at a rate of 0.5 °C/5 s for melting curve analysis.
The Leu410 PCR was performed based on the protocol described by Saavedra-Rodriguez et al. . Each reaction contained 4.5 µl of iQ-SYBR Green Supermix (Bio-Rad Laboratories Inc.), 0.45 µl of each primer (Val forward primer [5ʹGCG GGC AGG GCG GCG GGG GCG GGG CCA TCT TCT TGG GTT CGT TCT ACC GTG-3ʹ], Leu forward primer [5′-GCG GGC ATC TTC TTG GGT TCG TTC TAC CAT T-3′], and a common reverse primer [5′-TTC TTC CTC GGC GGC CTC TT-3′]), 1 µl of template DNA and ddH2O for a final reaction volume of 9 µl. Thermal cycling conditions were: 95 °C for 3 min; 40 cycles of 95 °C for 10 s, 60 °C for 10 s, 72 °C for 30 s; 95 °C for 10 s; and a ramp from 65 °C to 95 °C at a rate of 0.2 °C/10 s for melting curve analysis.
Percentage mortality at the diagnostic time of 30 min was used to describe the susceptibility status of the mosquito populations tested. Populations were classified as resistant and susceptible based on WHO and CDC guidelines [37, 38]: when mortality was < 90%, the population was considered to be resistant; when mortality was ≥ 98%, the population was considered to be susceptible; mortality ranging from 90 to 97% suggested that the population was developing resistance. Percentage mortality was calculated with the 95% CI from the bioassays and for allele frequencies.
Interquartile ranges of the mean of the optical density (OD) values from the biochemical assays were compared between study populations and the susceptible reference strain. Regression analyses were performed to measure the statistical significance of differences between the mean OD values between populations.
Pearson chi-square tests were performed to understand the associations between Gly1016 and Cys1534 genotypes and phenotypes of bioassayed mosquitoes. The population-level allele frequencies were calculated using the following equation :
The linkage disequilibrium, departures from the Hardy–Weinberg equilibrium (HWE), and the P value for Gly1016 and Cys1534 in each population were assessed using Fisher’s exact test in the population genetics software package GENEPOP version 4.2 (https://genepop.curtin.edu.au/) . Statistical analyses were conducted in Microsoft Excel 2016 (Microsoft Inc., Redmond, WA, USA) and Stata 15 (StataCorp LLC, College Station, TX, USA).
In the mosquito populations collected from Dhaka City, Ae. aegypti mortality ranged between 0% in Malibagh to 6.7% in Gulshan & Karail at the diagnostic dose of permethrin. A higher dose of permethrin (2× the diagnostic dose) was tested with the populations of Dhanmondi & Mohammadpur and Gulshan & Karail but still resulted in < 50% mortality at the diagnostic time point. In contrast, Ae. aegypti mortality at the diagnostic dose of deltamethrin varied between the sampling areas of Dhaka City, ranging from 49.0% (95% CI ± 7.3) in Gulshan & Karail to 100% (95% CI ± 1.6) in Uttara. Susceptibility to malathion was tested in three populations from Dhaka City: the Gulshan & Karail population was resistant (62.9% mortality, 95% CI ± 2.7), and the Dhanmondi & Mohammadpur (98.1% mortality, 95% CI ± 2.7) and Uttara (100% mortality, 95% CI ± 2.1) populations were susceptible. All Dhaka City populations tested against bendiocarb were susceptible (100% mortality in all populations) (Fig. 1).
The Ae. aegypti populations sampled from the non-Dhaka locations were also highly resistant to permethrin, with mortality ranging from 0% in the Chapai Nawabganj population to 14.8% (95% CI ± 2.0) in the Rajshahi population. When the concentration of permethrin was increased to 2× in Chittagong, mortality was still < 50%. However, when the permethrin concentration was increased to 5× in Bandarban, the population was fully susceptible (100% mortality). While the Chapai Nawabganj (100% mortality, CI ± 0.79) and Chittagong (99.0% mortality, 95% CI ± 0.79) populations were susceptible to deltamethrin, the Bandarban population was resistant to deltamethrin at the diagnostic dose (67% mortality, 95% CI ± 6.1) but susceptible when the concentration was increased to 2× (99.1% mortality, 95% CI ± 0.79). The Bandarban population was also resistant to the diagnostic dose of malathion (75.7% mortality, 95% CI ± 4.6) but susceptible to 2× higher malathion (100% mortality). The ROCK reference strain was confirmed to be fully susceptible to the diagnostic doses of the four insecticides. A summary of bioassay data is presented in Fig. 2.
All Ae. aegypti populations tested from field collections had significantly higher (P < 0.0001) MFO levels compared to the insecticide-susceptible ROCK reference strain. β-EST activity levels of Ae. aegypti populations from Azimpur, Uttara, Dhanmondi & Mohammadpur, Gulshan & Karail, Malibagh, Mirpur, and Bandarban were significantly (P < 0.0001) higher than those of the ROCK reference strain, but β-EST levels in Ae. aegypti populations from the non-Dhaka sites of Chapai Nawabganj, Chittagong, and Rajshahi were significantly lower than those of the ROCK strain (p < 0.0001). In the case of AChE activity, populations from Azimpur (P < 0.042), Chittagong (P < 0.019), and Gulshan & Karail (P < 0.0001) had significantly higher activity levels than the ROCK reference strain, and those from Dhanmondi & Mohammadpur (P < 0.001) and Malibagh (P < 0.001) had significantly lower activity levels than the ROCK strain. The estimated levels of IAChE were significantly higher (P < 0.0001) in the Gulshan & Karail population compared to the ROCK strain, which suggests that AChE insensitivity may exist in the former population. In comparison, IAChE levels were low across the remaining populations, suggesting that the target site remains sensitive. However, it is noteworthy that levels were significantly lower than those in the ROCK strain in Bandarban, Chapai Nawabganj, Mirpur, and Uttara (P < 0.001). When total protein content was compared between mosquito populations, with the exception of Azimpur and Mirpur, total protein content in all populations was significantly (P < 0.026) lower than in the ROCK strain, suggesting that body size was generally smaller for most of the field populations (Fig. 3).
Knockdown resistance (kdr) genotyping
A total of 142 phenotyped mosquitoes from bioassays at diagnostic dose (1×) permethrin and 59 phenotyped mosquitoes from deltamethrin 1X bioassays were analyzed for the Gly1016 mutation. From the Dhaka mosquito populations exposed to permethrin, 37.8% (28/74) of the survivors (assessed as being alive at the time points) were mutant homozygotes (GG) and 29.7% (22/74) were wild-type homozygotes (VV). The correlations between genotype and phenotype of permethrin-exposed Dhaka mosquitoes were statistically significant (P < 0.0001). Most of the dead mosquitoes were wild-type homozygotes (12/14, 85.7%). Among the mosquitoes from sites outside of Dhaka, more than half of the permethrin survivors were heterozygotes (23/44, 52.3%,) and there was an equal number (5/10, 50.0%) of wild-type homozygotes and heterozygotes among the dead mosquitoes. For the deltamethrin bioassays, only dead mosquitoes were genotyped due to limitations at the time of the bioassay. Similar genotype frequencies were seen in the mosquitoes from Dhaka that were dead after exposure to deltamethrin. However, the mosquitoes from outside of Dhaka did not include any mutant homozygotes (Table 2).
Of the 170 mosquitoes screened for the Cys1534 mutation, 110 mosquitoes were from permethrin bioassays and the remaining were from deltamethrin bioassays from both Dhaka and non-Dhaka populations. From the permethrin phenotyped Dhaka mosquitoes, 54.1% (33/61) of the resistant (surviving) mosquitoes were 1534 mutant homozygotes (CC) and 41.0% (25/61) were wild-type homozygotes (FF). In the case of permethrin-susceptible mosquitoes, 90.0% (9/10) were FF and the remaining individual was CC. From the non-Dhaka populations, 37.9% (11/29) of permethrin-resistant mosquitoes were CC and 27.6% (8/29) were heterozygotes (FC). Interestingly, none of the permethrin-susceptible mosquitoes from the non-Dhaka sites was FF, and eight of ten were CC. The correlations between genotype and phenotypes of permethrin exposed mosquitoes in both populations were statistically significant (P < 0.016 for Dhaka and P < 0.043 for non-Dhaka). A total of 60 dead mosquitoes from the deltamethrin bioassays were analyzed for Cys1534. Interestingly, most of the mosquitoes were wild-type homozygotes (FF; 19/30, 65.5%) in Dhaka, whereas the opposite was seen for non-Dhaka populations (Table 3).
All mosquitoes (n = 264) from the permethrin and deltamethrin bioassays (1× and 2×) genotyped for Leu410 were found to be wild-type homozygotes.
Of the 177 non-phenotyped mosquitoes from the Dhaka populations, more than half were V1016G heterozygotes (90/177, 51%). The highest Gly1016 homozygote (GG) frequency was observed in the Gulshan & Karail (23/30, 77%) population followed by the Mirpur (14/30, 47%) and Malibagh (11/29, 38%)populations (Fig. 4). In the case of Cys1534, the largest group of the mosquitoes were homozygous wild type (FF) (77/177, 43.5%,). The highest mutant homozygote (CC) frequency was recorded in the Dhanmondi & Mohammadpur population (2/29, 41.4%) (Fig. 5).
The overall allele frequency of Gly1016 and Cys1534 was 57.1% (95% CI ± 8.41) and 38.4% (95% CI ± 5.66), respectively. The Gulshan & Karail population had highest frequency of Gly1016, 85.0% (95% CI ± 30.42), and the Dhanmondi & Mohammadpur population had the highest frequency of Cys1534, 63.8% (95% CI ± 23.22) (Table 4). The HWE test revealed that three populations had significant departures from HWE: the Azimpur population for Gly1016, the Uttara population for Cys1534, and the Mirpur population for both (Table 4).
When all three kdr alleles were considered together across the phenotyped mosquitoes, the most common tri-locus genotype was VV/CC/VV (n = 53/158) (1016/1534/410) followed by GG/FF/VV (n = 45/158). When overall genotype was linked to phenotype, permethrin survivors were most commonly GG/FF/VV (n = 33/87). For mosquitoes that were phenotypically susceptible, the most common genotype for those killed by both deltamethrin and permethrin was VV/CC/VV (n = 24/73) (Fig. 6).
The application of chemical insecticides either in the form of space sprays, thermal fogging, or LLINs has been carried out for many years in Bangladesh. However, documented evidence of mosquito susceptibility to insecticides is scanty. Some information can be obtained from the ‘Malaria Threat Map’ website on insecticide resistance in some Anopheles species . A recent article reported permethrin and deltamethrin resistance in Anopheles vagus in the highest malarious region in Bangladesh . However, these reports are limited to phenotypic characteristics, and no clear understanding of resistance mechanisms for any mosquito species is available.
Despite the increasing prevalence of Aedes-borne diseases in Bangladesh, the insecticide resistance status of Ae. aegypti has previously not been assessed. The results reported here provide a comprehensive overview of insecticide resistance across Dhaka City and in several other sites of high epidemiological importance. We report a high frequency and intensity of permethrin resistance in all mosquito populations that were studied. However, despite this high level of permethrin resistance, susceptibility to deltamethrin was still present in several of the populations. This difference in susceptibility to different pesticides suggests that the underlying mechanisms causing resistant phenotypes in these populations may not be shared across the pyrethroid class of chemicals.
The increased activity of enzymes, including β-ESTs, and increasing levels of MFOs in the populations suggests an important role of metabolic mechanisms in conferring resistance. All of the Dhaka mosquito populations had elevated levels of oxidases and increased esterase activity and they were resistant to permethrin. Outside of Dhaka, esterase activity was notably lower in Chapai Nawabganj and Rajshahi, and while both populations were resistant to permethrin, the latter population remained susceptible to deltamethrin (the former was not tested). Increased activity of esterases and increased levels of oxidases may also be associated with the malathion resistance that was detected previously in the Gulshan & Karail and Bandarban populations [45,46,47,48]. In addition, AChE activity was elevated in the Gulshan & Karail population, possibly also contributing to the malathion resistance that was detected there. An important limitation of the biochemical assay data presented in this study is the lack of rigorous normalization of the data based on total protein content. As such, all comparisons were conducted relative to the insecticide-susceptible ROCK control strain. However, the fact that almost all of the field populations had significantly lower total protein content than the susceptible ROCK strain against which they were compared suggests that the increased enzyme levels and activities that were detected could not simply be attributed to greater protein content in the field populations. Additional information on alpha esterases and GSTs could be a focus of future studies. A growing body of evidence suggests that these are important mechanisms in pyrethroid resistance in Ae. aegypti , with GSTe2 associated with resistance to both permethrin and deltamethrin, and GSTe7 associated with resistance to deltamethrin [50, 51].
The kdr mutations Gly1016 and Cys1534 were found at varying frequencies in the mosquito populations collected across Dhaka City. This fine-scale spatial heterogeneity suggests that selection pressures for insecticide resistance are variable across small spatial scales within Dhaka City, reflecting trends that have been reported elsewhere [42, 52]. Historically, Aedes control in Dhaka and other major cities in Bangladesh solely depends on thermal fogging using a combination of pyrethroid insecticides. Pyrethroids are also commonly used in households via commercially available coils and aerosols. Both operational and domestic insecticide use may contribute to insecticide resistance selection pressures in Ae. aegypti .
The Gly1016 and Cys1534 kdr mutations have been widely reported in mosquito populations in Asia [40, 54, 55]. An additional mutation, Leu410, has also been reported in association with pyrethroid resistance, but its prevalence in Asia has not yet been well studied . Expression of insect sodium channels in Xenopus oocytes coupled with electro-neurophysiological measurements has demonstrated that Gly1016, Cys1534, and Leu410 reduce the sensitivity of the VGSC to permethrin and deltamethrin [23, 56]. However, Leu410 was not detected in any of the mosquito populations in the current study. This was an unexpected result, as previous research has suggested the parallel evolution of this mutation together with the polymorphisms at positions 1016 and 1534 , both of which were detected at moderate to high frequencies in our study. The co-occurrence of Pro989 with Gly1016 conferring high pyrethroid resistance in Ae. aegypti has been reported previously . However, our study did not any analysis of the S989P kdr mutation.
In Dhaka City, 1016G and 1534C homozygous mutants were mostly associated with survival in the permethrin bioassays. It is also worth noting that the population from the Dhaka City neighborhood of Mirpur was resistant to permethrin yet susceptible to deltamethrin and was also the population with the highest frequency of Val1016 and Phe1534 wild-type homozygotes. These findings suggest that while kdr alleles may be contributing to the insecticide resistance that was detected, they are not the only mechanism and such relationship is not rare [57, 58].
From an operational perspective, the data presented here will be important in guiding the choice of vector control tools. Given the widespread and intense permethrin resistance that was detected, vector control products containing alternative compounds should be used. Although some populations have remained susceptible to deltamethrin, given the high degree of permethrin resistance, it would be prudent to search for alternatives outside of the pyrethroid class. Particularly notable was the detection of deltamethrin resistance in mosquitoes collected in the Bandarban region, where deltamethrin-treated bed nets are routinely used for malaria control . Bandarban was also the only non-Dhaka site to show significantly elevated esterase activity, suggesting that the population was experiencing comparatively greater selective pressure across multiple mechanisms as compared to the other non-Dhaka sites. Vector control activities have focused largely on malaria vectors and have not routinely targeted Aedes in this part of Bangladesh. The finding that the Aedes population was resistant to the insecticide relied upon for malaria control highlights the importance of implementing strategies based on integrated vector management in Bandarban.
The only insecticide to which every population tested was susceptible was bendiocarb. However, there is no product registered in Bangladesh that could be employed for Aedes control that contains bendiocarb as an active ingredient. The next best candidate was malathion, with public health agencies desperately seeking alternatives to pyrethroids. Nevertheless, malathion resistance was detected in several of the populations studied, both inside and outside of Dhaka city. Also, malathion has been used in agriculture for many years in Bangladesh, so selection pressure outside of vector control already exists to a certain degree [59, 60]. In such a scenario as we detected in Bangladesh with a patchwork of insecticide-resistant phenotypes, it will be challenging to find a ‘one size fits all’ solution for Aedes control.
This current study provides evidence of insecticide resistance in Ae. aegypti and data on resistance mechanisms, including detoxification enzymes and kdr mutations, in Bangladesh. High pyrethroid resistance may be compromising the existing Aedes control strategies, and the presence of multiple resistance mechanisms poses further challenges regarding alternatives. Continuous surveillance of insecticide resistance will enable trends in susceptibility to be monitored over space and time and will provide a more robust evidence base upon which to select the most effective vector control tools and strategies. In cities like Dhaka, where operational control faces challenges posed by insecticide resistance, in addition to the rational use of chemicals, sustainable and alternative tools like biocontrol approaches should be considered.
The preliminary results were disseminated among different stakeholders and mosquito control authorities immediately after the data were analyzed. Followed by the outbreak of dengue during the monsoon season of 2019, these research findings and recommendations were reinvestigated by the policymakers. As a result, permethrin has been replaced by malathion for the control of adult mosquitoes in Dhaka city [61, 62].
Indoor residual spraying
Insensitive acetylcholine esterase
International Centre for Diarrhoeal Disease Research, Bangladesh
Long-lasting insecticidal nets
U.S. Centers for Disease Control and Prevention
Voltage-gated sodium channel
Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, et al. The global distribution and burden of dengue. Nature. 2013;496(7446):504–7.
Sissoko D, Malvy D, Ezzedine K, Renault P, Moscetti F, Ledrans M, et al. Post-epidemic chikungunya disease on Reunion Island: course of rheumatic manifestations and associated factors over a 15-month period. PLoS Negl Trop Dis. 2009;3(3):e389.
Couturier E, Guillemin F, Mura M, Leon L, Virion JM, Letort MJ, et al. Impaired quality of life after chikungunya virus infection: a 2-year follow-up study. Rheumatology. 2012;51(7):1315–22.
Paixao ES, Barreto F, Teixeira Mda G, Costa Mda C, Rodrigues LC. History, epidemiology, and clinical manifestations of Zika: a systematic review. Am J Public Health. 2016;106(4):606–12.
Banu S, Hu W, Hurst C, Guo Y, Islam MZ, Tong S. Space-time clusters of dengue fever in Bangladesh. Trop Med Int Health. 2012;17(9):1086–91.
Sharmin S, Viennet E, Glass K, Harley D. The emergence of dengue in Bangladesh: epidemiology, challenges and future disease risk. Trans R Soc Trop Med Hyg. 2015;109(10):619–27.
Salje H, Paul KK, Paul R, Rodriguez-Barraquer I, Rahman Z, Alam MS, et al. Nationally-representative serostudy of dengue in Bangladesh allows generalizable disease burden estimates. Elife. 2019;8:e42869.
Institute of Epidemiology, Disease Control and Research. Dengue current situation. https://www.iedcr.gov.bd/index.php/dengue/219. Accessed 26 Apr 2020.
Hossain MS, Hasan MM, Islam MS, Islam S, Mozaffor M, Khan MAS, et al. Chikungunya outbreak (2017) in Bangladesh: clinical profile, economic impact and quality of life during the acute phase of the disease. PLoS Negl Trop Dis. 2018;12(6):e0006561.
Muraduzzaman AKM, Sultana S, Shirin T, Khatun S, Islam M, Rahman M. Introduction of Zika virus in Bangladesh: an impending public health threat. Asian Pac J Trop Med. 2017;10(9):925–8.
GeurtsvanKessel CH, Islam Z, Islam MB, Kamga S, Papri N, van de Vijver D, et al. Zika virus and Guillain-Barre syndrome in Bangladesh. Ann Clin Transl Neurol. 2018;5(5):606–15.
Aktaruzzaman MM. Dissemination workshop on monsoon Aedes survey in Dhaka North & South City Corporations. Directorate General of Health Services (DGHS); 2018.
Paul KK, Dhar-Chowdhury P, Haque CE, Al-Amin HM, Goswami DR, Kafi MAH, et al. Risk factors for the presence of dengue vector mosquitoes, and determinants of their prevalence and larval site selection in Dhaka, Bangladesh. PLoS ONE. 2018;13(6):e0199457.
World Health Organization. Dengue: guidelines for diagnosis, treatment, prevention and control. Geneva: World Health Organization; 2009.
Marcombe S, Darriet F, Tolosa M, Agnew P, Duchon S, Etienne M, et al. Pyrethroid resistance reduces the efficacy of space sprays for dengue control on the island of Martinique (Caribbean). PLoS Negl Trop Dis. 2011;5(6):e1202.
Maciel-de-Freitas R, Avendanho FC, Santos R, Sylvestre G, Araujo SC, Lima JB, et al. Undesirable consequences of insecticide resistance following Aedes aegypti control activities due to a dengue outbreak. PLoS ONE. 2014;9(3):e92424.
Plernsub S, Saingamsook J, Yanola J, Lumjuan N, Tippawangkosol P, Walton C, et al. Temporal frequency of knockdown resistance mutations, F1534C and V1016G, in Aedes aegypti in Chiang Mai city, Thailand and the impact of the mutations on the efficiency of thermal fogging spray with pyrethroids. Acta Trop. 2016;162:125–32.
Moyes CL, Vontas J, Martins AJ, Ng LC, Koou SY, Dusfour I, et al. Contemporary status of insecticide resistance in the major Aedes vectors of arboviruses infecting humans. PLoS Negl Trop Dis. 2017;11(7):e0005625.
Dusfour I, Vontas J, David JP, Weetman D, Fonseca DM, Corbel V, et al. Management of insecticide resistance in the major Aedes vectors of arboviruses: advances and challenges. PLoS Negl Trop Dis. 2019;13(10):e0007615.
Devine GJ, Denholm I. Insecticide and acaricide resistance. In: Resh V, Carde R, editors. Encyclopedia of insects. San Diego: Academic Press; 2009. p. 505–11.
World Health Organization. Techniques to detect insecticide resistance mechanisms (field and laboratory manual). Geneva: World Health Organization; 1998.
Dong K, Du Y, Rinkevich F, Nomura Y, Xu P, Wang L, et al. Molecular biology of insect sodium channels and pyrethroid resistance. Insect Biochem Mol Biol. 2014;50:1–17.
Haddi K, Tome HVV, Du Y, Valbon WR, Nomura Y, Martins GF, et al. Detection of a new pyrethroid resistance mutation (V410L) in the sodium channel of Aedes aegypti: a potential challenge for mosquito control. Sci Rep. 2017;7:46549.
Hemingway J, Ranson H. Insecticide resistance in insect vectors of human disease. Annu Rev Entomol. 2000;45:371–91.
Usuga AF, Zuluaga-Idarraga LM, Alvarez N, Rojo R, Henao E, Rua-Uribe GL. Barriers that limit the implementation of thermal fogging for the control of dengue in Colombia: a study of mixed methods. BMC Public Health. 2019;19(1):669.
Staff Correspondence. Bangladesh set an example of failure in dengue prevention. Dhaka: The Daily Star; 2019.
Alam H. No fund for mosquito control activities. Dhaka: The Daily Star; 2019.
Anonymous. First identified outbreak of chikungunya in Bangladesh, 2008, vol. 7. Dhaka: icddr,b; 2009. p. 1–6.
Khatun S, Chakraborty A, Rahman M, Nasreen Banu N, Rahman MM, Hasan SM, et al. An outbreak of chikungunya in rural Bangladesh, 2011. PLoS Negl Trop Dis. 2015;9(7):e0003907.
Salje H, Lessler J, Paul KK, Azman AS, Rahman MW, Rahman M, et al. How social structures, space, and behaviors shape the spread of infectious diseases using chikungunya as a case study. Proc Natl Acad Sci USA. 2016;113(47):13420–5.
Haque F, Rahman M, Banu NN, Sharif AR, Jubayer S, Shamsuzzaman A, et al. An epidemic of chikungunya in northwestern Bangladesh in 2011. PLoS ONE. 2019;14(3):e0212218.
Hsan K, Hossain MM, Sarwar MS, Wilder-Smith A, Gozal D. Unprecedented rise in dengue outbreaks in Bangladesh. Lancet Infect Dis. 2019;19(12):1287.
Chowdhury R, Chowdhury V, Faria S, Islam S, Maheswary NP, Akhter S, et al. Indoor residual spraying for kala-azar vector control in Bangladesh: a continuing challenge. PLoS Negl Trop Dis. 2018;12(10):e0006846.
Alam MS, Al-Amin HM, Khan WA, Haque R, Nahlen BL, Lobo NF. Preliminary report of pyrethroid resistance in Anopheles vagus, an important malaria vector in Bangladesh. Am J Trop Med Hyg. 2020;103(2):810–11. https://doi.org/10.4269/ajtmh.20-0146.
Anonymous. Chikungunya: 23 areas at higher risk in Dhaka. In: bdnews24com. Dhaka: bdnews24.com; 2017.
Costa-da-Silva AL, Navarrete FR, Salvador FS, Karina-Costa M, Ioshino RS, Azevedo DS, et al. Glytube: a conical tube and parafilm M-based method as a simplified device to artificially blood-feed the dengue vector mosquito, Aedes aegypti. PLoS ONE. 2013;8(1):e53816.
Brogdon W, Chan A. Guidelines for evaluating insecticide resistance in vectors using the CDC bottle bioassay/methods in Anopheles research. CDC Technical Report. U.S. Centers for Disease Control and Prevention; 2010.
World Health Organization. Test procedures for insecticide resistance monitoring in malaria vector mosquitoes, 2 edn. Geneva: World Health Organization; 2016.
Saavedra-Rodriguez K, Urdaneta-Marquez L, Rajatileka S, Moulton M, Flores AE, Fernandez-Salas I, et al. A mutation in the voltage-gated sodium channel gene associated with pyrethroid resistance in Latin American Aedes aegypti. Insect Mol Biol. 2007;16(6):785–98.
Yanola J, Somboon P, Walton C, Nachaiwieng W, Somwang P, Prapanthadara LA. High-throughput assays for detection of the F1534C mutation in the voltage-gated sodium channel gene in permethrin-resistant Aedes aegypti and the distribution of this mutation throughout Thailand. Trop Med Int Health. 2011;16(4):501–9.
Saavedra-Rodriguez K, Maloof FV, Campbell CL, Garcia-Rejon J, Lenhart A, Penilla P, et al. Parallel evolution of vgsc mutations at domains IS6, IIS6 and IIIS6 in pyrethroid resistant Aedes aegypti from Mexico. Sci Rep. 2018;8(1):6747.
Deming R, Manrique-Saide P, Medina Barreiro A, Cardena EU, Che-Mendoza A, Jones B, et al. Spatial variation of insecticide resistance in the dengue vector Aedes aegypti presents unique vector control challenges. Parasites Vectors. 2016;9:67.
Rousset F. genepop’007: a complete re-implementation of the genepop software for Windows and Linux. Mol Ecol Resour. 2008;8(1):103–6.
World Health Organization. Vector insecticide resistance. https://apps.who.int/malaria/maps/threats/?theme=prevention&mapType=prevention%3A0&bounds=%5B%5B87.66488246615313%2C20.717484946543323%5D%2C%5B98.56829160856773%2C26.61260734428494%5D%5D&insecticideClass=PYRETHROIDS&insecticideTypes=&assayTypes=MOLECULAR_ASSAY%2CBIOCHEMICAL_ASSAY%2CSYNERGIST-INSECTICIDE_BIOASSAY&synergistTypes=&species=&vectorSpecies=&surveyTypes=&deletionType=HRP2_PROPORTION_DELETION&plasmodiumSpecies=P._FALCIPARUM&drug=DRUG_AL&mmType=1&endemicity=false&countryMode=false&storyMode=false&storyModeStep=0&filterOpen=true&filtersMode=filters&years=2010%2C2020. Accessed 03 July 2020.
Paeporn P, Supaphathom K, Srisawat R, Komalamisra N, Deesin V, Ya-umphan P, et al. Biochemical detection of pyrethroid resistance mechanism in Aedes aegypti in Ratchaburi province, Thailand. Trop Biomed. 2004;21(2):145–51.
Rodriguez MM, Bisset JA, Fernandez D. Levels of insecticide resistance and resistance mechanisms in Aedes aegypti from some Latin American countries. J Am Mosq Control Assoc. 2007;23(4):420–9.
Alvarez LC, Ponce G, Oviedo M, Lopez B, Flores AE. Resistance to malathion and deltamethrin in Aedes aegypti (Diptera: Culicidae) from western Venezuela. J Med Entomol. 2013;50(5):1031–9.
Ngoagouni C, Kamgang B, Brengues C, Yahouedo G, Paupy C, Nakoune E, et al. Susceptibility profile and metabolic mechanisms involved in Aedes aegypti and Aedes albopictus resistant to DDT and deltamethrin in the Central African Republic. Parasites Vectors. 2016;9(1):599.
Aponte A, Penilla RP, Rodriguez AD, Ocampo CB. Mechanisms of pyrethroid resistance in Aedes (Stegomyia) aegypti from Colombia. Acta Trop. 2019;191:146–54.
Lumjuan N, McCarroll L, Prapanthadara L-A, Hemingway J, Ranson H. Elevated activity of an Epsilon class glutathione transferase confers DDT resistance in the dengue vector, Aedes aegypti. Insect Biochem Mol Biol. 2005;35:861–71.
Lumjuan N, Rajatileka S, Changsom D, Wicheer J, Leelapat P, Prapanthadara LA, et al. The role of the Aedes aegypti Epsilon glutathione transferases in conferring resistance to DDT and pyrethroid insecticides. Insect Biochem Mol Biol. 2011;41(3):203–9.
Brito LP, Carrara L, de Freitas RM, Lima JBP, Martins AJ. Levels of resistance to pyrethroid among distinct kdr alleles in Aedes aegypti laboratory lines and frequency of kdr alleles in 27 natural populations from Rio de Janeiro, Brazil. Biomed Res Int. 2018;2018:2410819.
Gray L, Florez SD, Barreiro AM, Vadillo-Sanchez J, Gonzalez-Olvera G, Lenhart A, et al. Experimental evaluation of the impact of household aerosolized insecticides on pyrethroid resistant Aedes aegypti. Sci Rep. 2018;8(1):12535.
Stenhouse SA, Plernsub S, Yanola J, Lumjuan N, Dantrakool A, Choochote W, et al. Detection of the V1016G mutation in the voltage-gated sodium channel gene of Aedes aegypti (Diptera: Culicidae) by allele-specific PCR assay, and its distribution and effect on deltamethrin resistance in Thailand. Parasites Vectors. 2013;6(1):253.
Ishak IH, Jaal Z, Ranson H, Wondji CS. Contrasting patterns of insecticide resistance and knockdown resistance (kdr) in the dengue vectors Aedes aegypti and Aedes albopictus from Malaysia. Parasites Vectors. 2015;8:181.
Hirata K, Komagata O, Itokawa K, Yamamoto A, Tomita T, Kasai S. A single crossing-over event in voltage-sensitive Na+ channel genes may cause critical failure of dengue mosquito control by insecticides. PLoS Negl Trop Dis. 2014;8(8):e3085.
Kasai S, Komagata O, Itokawa K, Shono T, Ng LC, Kobayashi M, et al. Mechanisms of pyrethroid resistance in the dengue mosquito vector, Aedes aegypti: target site insensitivity, penetration, and metabolism. PLoS Negl Trop Dis. 2014;8(6):e2948.
Martins AJ, Lins RM, Linss JG, Peixoto AA, Valle D. Voltage-gated sodium channel polymorphism and metabolic resistance in pyrethroid-resistant Aedes aegypti from Brazil. Am J Trop Med Hyg. 2009;81(1):108–15.
Hasanuzzaman M, Rahman MA, Salam MA. Identification and quantification of pesticide residues in water samples of Dhamrai Upazila, Bangladesh. Appl Water Sci. 2017;7(6):2681–8.
Numan A, Khan M, Uddin R, Rahman M, Bhuiyan M, Akter N. Risk assessment of commonly used major pesticides for tomato (Solanum lycopersicum L.) cultivation in Bangladesh. Adv Nutr Fd Sci: ANAFS-109. 2018.
Alam H. Insecticides ineffective on mosquitoes: icddr,b. Dhaka: The Daily Star; 2019.
Star Online Report. DNCC to spray new insecticide from tomorrow to kill mosquito. Dhaka: The Daily Star; 2019.
icddr,b acknowledges with gratitude the commitment of CDC to its research efforts. icddr,b is also grateful to the Governments of Bangladesh, Canada, Sweden, and the UK for providing core/unrestricted support.
The findings and conclusions in this paper are those of the authors and do not necessarily represent the official position of the CDC.
This study was funded by U.S. Centers for Disease Control and Prevention (CDC, Atlanta, GA, USA), Grant Number 3U01GH001207-03S1.
Ethics approval and consent to participate
The study was approved by the icddr,b Ethical Review Committee (ERC), PR-17050. Informed oral consent was obtained from the household owners before ovitrap installation.
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.
About this article
Cite this article
Al-Amin, H.M., Johora, F.T., Irish, S.R. et al. Insecticide resistance status of Aedes aegypti in Bangladesh. Parasites Vectors 13, 622 (2020). https://doi.org/10.1186/s13071-020-04503-6
- Aedes aegypti
- Insecticide resistance