Skip to main content

Impact of deltamethrin selection on kdr mutations and insecticide detoxifying enzymes in Aedes aegypti from Mexico

Abstract

Background

Insecticide resistance is a serious problem for vector control programmes worldwide. Resistance is commonly attributed to mutations at the insecticide’s target site or increased activity of detoxification enzymes.

Methods

We determined the knockdown concentration (KC50) and lethal concentration (LC50) of deltamethrin in six natural populations of adult Aedes aegypti from southeastern Mexico. These populations were then selected over five generations using the LC50 from the preceding generation that underwent selection, and the heritability of deltamethrin resistance was quantified. For each generation, we also determined the frequency of the kdr alleles L410, I1016 and C1534, and the levels of activity of three enzyme families (α- and β-esterases, mixed-function oxidases and glutathione S-transferases (GST)) associated with insecticide detoxification.

Results

There was an increase in KC50 and LC50 values in the subsequent generations of selection with deltamethrin (FS5vs FS0). According to the resistance ratios (RRs), we detected increases in LC50 ranging from 1.5 to 5.6 times the values of the parental generation and in KC50 ranging from 1.3–3.8 times the values of the parental generation. Triple homozygous mutant individuals (tri-locus, LL/II/CC) were present in the parental generations and increased in frequency after selection. The frequency of L410 increased from 1.18-fold to 2.63-fold after selection with deltamethrin (FS5vs FS0) in the populations analyzed; for I1016 an increase between 1.19-fold to 2.79-fold was observed, and C1534 was fixed in all populations after deltamethrin selection. Enzymatic activity varied significantly over the generations of selection. However, only α- esterase activity remained elevated in multiple populations after five generations of deltamethrin selection. We observed an increase in the mean activity levels of GSTs in two of the six populations analyzed.

Conclusions

The high levels of resistance and their association with high frequencies of kdr mutations (V410L, V1016I and F1534C) obtained through artificial selection, suggest an important role of these mutations in conferring resistance to deltamethrin. We highlight the need to implement strategies that involve the monitoring of kdr frequencies in insecticide resistance monitoring and management programmes.

Background

Aedes aegypti (L.) is an arbovirus vector of great public health importance, and is the principal vector transmitting dengue, chikungunya and Zika viruses in the Americas [1,2,3]. Currently, there are no widely licensed vaccines or specific treatments for any of these diseases, rendering vector control the principal strategy for preventing their transmission. The reduction or elimination of larval habitats and chemical control with the use of larvicides and adulticides are the main strategies to reduce vector-human contact and break the transmission cycles of these diseases [4]. The insecticides used for such purposes worldwide are mainly from the following chemical groups: organophosphates, organochlorines, carbamates and pyrethroids [5, 6]. Pyrethroids are often the first line of defense to break the transmission of these diseases, and they are the most commonly used worldwide due to their low toxicity to mammals including humans [7, 8].

In Mexico, the application of pyrethroid insecticides is regulated by the Ministry of Health, and their use dates back to 1999 when the use of DDT was banned in the health sector. At this time, permethrin became the insecticide of choice for space spraying and deltamethrin for indoor residual spraying [9]. In 2011, the use of permethrin was discontinued; however, deltamethrin remained on the list of insecticides approved for use as an adulticide in indoor residual spraying. In 2012 other pyrethroids were added to the list; however, deltamethrin remained, and in 2016 its use for space spraying applications was also approved, a practice that remains to date [10]. The use of pyrethroids in Mexico has been favored for more than 20 years, which led to widespread resistance in mosquito populations [11,12,13].

Insecticide resistance is one of the main obstacles to the success of vector control programmes [12, 14, 15]. Two of the principal mechanisms that cause resistance in mosquitoes are alterations at the target site of the insecticide, including knockdown resistance (kdr) mutations on the voltage-gated sodium channel gene (vgsc), and increased metabolic activity [16, 17]. Target-site resistance arising from kdr mutations can cause resistance to both DDT and pyrethroids, since both types of insecticides target sodium channels at the axon level [18].

Since the discovery of the L1014F mutation in Musca domestica and since the determination of its role in insecticide resistance [19], more than 50 mutations have been described for the vgsc gene in different pests and disease vectors [20]. Thirteen of these sequence changes are commonly found in Ae. aegypti in ten different positions covering five regions of the vgsc gene (IIS4-5, IIS5-6, IIS-6, IIIS6 and IVS5), identified at positions 410 (V → L), 419 (V → L), 923 (G → I), 982 (L → W), 989 (S → P), 1011 (I → M or V), 1016 (V → G or I), 1520 (T → I), 1534 (F → C or L), and 1763 (D → Y) [21,22,23,24,25,26,27,28,29]. It should be noted that only five of the mutations described for this species have been functionally associated with sensitivity to pyrethroids and DDT by heterologous expression and electrophysiology assays; only S989P, I1011M, V1016G, F1534C and V410L have been confirmed [20, 27, 30, 31]. Similarly, it is important to note that mutations are usually associated with a specific geographic sector, for example, mutation V1016G has been described in Southeast Asia but not in the Americas [32].

In Mexico, mutations F1534C, V1016I and V410L have been reported in populations of Ae. aegypti in different states, sites where their temporal distribution indicates the absence of mutant alleles in 2000, with a gradual increase in these mutations during the period of 2002–2008, up to fixation of the mutant alleles in 2016, showing an allelic frequency for the three simultaneous changes of 0.47 [30]. Of these mutations, the change from phenylalanine to cysteine at position 1534 was described initially by [25], who associated its occurrence with pyrethroid resistance; however, confirmation by heterologous expression and electrophysiology has demonstrated that this change reduces sensitivity of VGSC to type I but not type II pyrethroids [33]. The presence of this mutation has been linked to the development of low resistance levels and its association with other gene changes to the establishment of a more resistant phenotype. Such is the case of its association with V1016I, described by Saavedra et al. [22], where its participation in the reduction of sensitivity of VGSC has not been identified when appearing as a single mutation [31]. An analysis of the frequency of V1016I and F1534C in 24 populations of Ae. aegypti in Mexico showed that the most likely evolutionary process is the appearance of F1534C and then V1016I, this being corroborated by the low fitness of haplotype I1016/F1534 [34] and the increase in insecticide resistance expressed in Xenopus oocytes when inducing the two mutations [31]. Temporal distribution studies of this mutation in Mexico identified an increase in the frequency of the I1016 mutant allele during the period of 1996 to 2000 with a frequency of 0.04% up to an increase to 33.2% in 2007 to 2009 [35].

The V410L mutation has been recently described by Haddi et al. [27], who determined their association with resistance to type I and II pyrethroids. This mutation has been studied in Mexico in populations of Ae. aegypti collected between 2000 and 2016 where its appearance was reported in 2002 in heterozygous individuals collected in Coatzacoalcos, Veracruz, and where the increase in the frequency of L410 was also documented as being as high as 0.9 in populations collected in 2014 in Merida, Yucatan [30]. A similar pattern was found for 26 populations of Ae. aegypti collected in the eastern part of Mexico, showing a frequency of the L410 allele close to the one in populations from Minatitlan and Jose Cardel in Veracruz and Cancun in Quintana Roo (0.99, 0.97 and 0.93, respectively), but the average values, as well as the interval of frequencies observed in populations ranged from 0.3 to 0.99 [36].

Metabolic resistance to pyrethroids is mediated mainly by glutathione S-transferases (GSTs), esterases and mixed-function oxidases (MFOs) [37,38,39]. In Mexico, Flores et al. [14] previously reported high levels of α- and β-esterase activity related to permethrin selection in populations of Ae. aegypti from Baja California Norte and Sur. The increased activity of these enzymes has also been reported in populations from the state of Sonora in a study of permethrin resistance and in the state of Veracruz in a study of chlorpyrifos resistance [40, 41]. A further study reported high levels of GST activity in the state of Guerrero, which was found to be related to resistance to DDT [12].

Given the background levels of pyrethroid and DDT resistance in many parts of Mexico, the extent to which ongoing insecticide pressure selects for specific mechanisms of resistance is not well understood. To better understand how insecticide exposure selects for resistance we evaluate the effect of laboratory selection with deltamethrin on the underlying resistance mechanisms in populations of Ae. aegypti from southeastern Mexico, over five generations of selection.

Methods

Biological material

Eggs and larvae of Ae. aegypti were collected during 2014 from six sites: Merida, Progreso, Hunucma, and Hoctun in the State of Yucatan; and Agua Dulce and Jose Cardel, in the State of Veracruz (Table 1, Fig. 1). The criteria for site collection were: high incidence of dengue and frequent application of insecticides for its control; high infestation index of homes with Ae. aegypti; and previous registration of insecticide resistance. The collection of the entomological material was carried out in peridomestic sites. In each, entomological inspections were performed for the collection of stages of Ae. aegypti in at least 20 potential breeding places per site, such as containers, tires and flowerpots.

Table 1 Collection sites, geographical coordinates
Fig. 1
figure1

Collection sites of Aedes aegypti from Veracruz and Yucatan states from Mexico

The biological material obtained from the field was reared under insectary conditions at 26 ± 2 °C and 70–80% relative humidity with a 12h:12h (light:dark) photoperiod. The adults obtained from field-collected material were allowed to intermate, and the eggs resulting were designated FS0 (without previous selection). The New Orleans strain (NO) was used as a susceptible reference in the study, this strain was originally obtained from the CDC (Atlanta, GA, USA) and has been maintained since 2002.

Bioassays

The bioassays consisted of exposing 20–25 non-blood-fed, 2–3-day-old FS0 female mosquitoes to different concentrations of the pyrethroid deltamethrin (> 98% purity; Chem Service, West Chester, PA, USA) based on the bottle bioassay methodology described by the CDC [42]. Each bioassay consisted of bottles containing different concentrations of deltamethrin resulting in mortalities between 9–90%, with at least three replicates for each concentration, and an untreated control bottle. The numbers of knocked-down mosquitoes were recorded at 1 h. After 1 h of exposure, all mosquitoes were gently transferred to a recovery container without insecticide and were offered a cotton ball soaked in a sugar solution. Mortality was recorded at 24 h. Both the bottles and recovery containers were kept at 24 ± 2 °C and 70% RH.

Rates of kdr and 24 h recovery were analyzed using a logistic regression model, QCal (https://sourceforge.net/projects/irmaproj/files/Qcal/) [43]. We determined the KC50 (concentration causing 50% knockdown) after 1 h of exposure. The LC50 (concentration causing 50% mortality) was estimated from mortality data 24 h after recovery. The confidence intervals were calculated using a significance level of α = 0.05. The mortalities were corrected according to Abbott’s formula [44] when mortality was observed in the control bottles. The resistance ratios (RR) were calculated by dividing the KC50 or LC50 by the KC50 or LC50 of the susceptible New Orleans (NO) Ae. aegypti reference strain. The magnitude of resistance was classified as high (RR > 10-fold), moderate (RR between 5–10-fold) or low (RR < 5-fold) according to the criteria proposed by Mazarri & Georghiou [45].

Selection with deltamethrin

Selection cohorts were generated according to the methodology described by Saavedra et al. [46]. The selection was carried out exposing 650–1000 FS1 males and females (1:4) for 1 h to the LC50 obtained for the previous generation (FS0) by the method described above. The survivors at 24 h were transferred to mosquito cages for breeding the subsequent generation (FS2). The dose response parameters were calculated again for the FS1 cohort and were used to select the next generation (FS2). This procedure was repeated for each generation of selection until FS5. The New Orleans strain was tested at the same time as the FS0 and was not tested in each generation of selection.

KC50 and LC50 values were compared between generations to monitor changes in resistance. The KC50 and LC50 values were considered significantly different if their 95% confidence intervals did not overlap.

Realized heritability

Resistance risk assessment was made by calculating realized heritability values (h2) in all selection cohorts for the knockdown (KC50) and lethal (LC50) parameters as described by Tabashnik [47]. The heritability index (h2= R/S) is calculated as the ratio of the response to selection (R) to the selection differential (S) according to the artificial selection technique of Falconer [48] and is related to the additive genetic variance for a trait. A low h2 predicts no additive genetic variance for a trait and a poor or very slow response to artificial selection, while a high h2 predicts a large additive genetic variance at one or a few loci that govern a trait and predicts a rapid response to artificial selection.

The selection differential is expressed as the product of selection intensity (i) and phenotypic standard deviation (S= iσ). The response to selection was determined for each generation as the difference in population means between subsequent generations with the probit analysis: R = log (Final KC50 or LC50) − log (Initial KC50 or LC50)/n, where final KC50 or LC50 area the values of the offspring after n generations of selection and initial KC50 or LC50 are the values of the parental generation before n generations of selection. The selection intensity was calculated from the proportion of surviving individuals in the examined population. The difference between KC50 or LC50 was calculated on a logarithmic scale because the logarithm of tolerance was assumed to be normally distributed. The phenotypic standard deviation at t-th generation (σt) was obtained as the inverse of the regression slope: σt= 1/bt. The parameters (R, i and σ) were determined at every generation and the realized heritability was estimated as the regression coefficient of cumulative responses on cumulative selection differentials [48].

Molecular assays

DNA was extracted from individual mosquitoes by the technique described by Coen et al. [49], and the DNA pellet was resuspended in 30 μl of ultrapure molecular grade water (Corning CellgroTM, Manassas, VA, USA). The quantity and quality of DNA were determined using a NanoDrop spectrophotometer 2000 (Thermo Fisher Scientific, Woonsocket, RI, USA).

The kdr alleles 1534C, 1016I and 410L were detected by endpoint PCR using a Bio-Rad (Hercules, CA, USA) T100TM thermocycler using the primers for each of the mutations described by Saavedra et al. [22, 30] and Yanola et al. [25].

The primers used to detect the V1016I mutation were: V1016fw (5′-GCG GGC AGG GCG GCG GGG GCG GGG CCA CAA ATT GTT TCC CAC CCG CAC CGG-3′); I1016fw (5′-GCG GGC ACA AAT TGT TTC CCA CCC GCA CTG A-3′); and I1016R (5′-TGA TGA ACC SGA ATT GGA CAA AAG C-3′). The PCR was carried out in a reaction mixture of 12.5 µl containing 1.25 µl of 10× buffer (Invitrogen, Carlsbad, CA, USA), 1.5 mM MgCl2, 0.2 mM dNTPs, 5 µM of each primer, 100 ng genomic DNA, and 2 U Taq polymerase (Invitrogen). The PCR reaction conditions were: 95 °C for 5 min; 29 cycles of 95 °C for 1 min, 60 °C for 1 min, 72 °C for 1 min 15 s; and a final extension step at 72 °C for 10 min [22].

The primers used to detect the F1534C mutation were: C1534fw (5′-GCG GGC AGG GCG GCG GGG GCG GGG CCT CTA CTT TGT GTT CTT CAT CAT GTG-3′); F1534fw (5′-GCG GGC TCT ACT TTG TGT TCT TCA TCA TAT T-3′); and F1534R (5′-TCT GCT CGT TGA AGT TGT CGA T-3′). The PCR was carried out in a reaction mixture of 12.5 µl containing 1.25 µl of 10× buffer (Invitrogen), 1.5 mM MgCl2, 0.1 mM dNTPs, 0.5 µM of each primer, 100 ng genomic DNA, and 1 U Taq polymerase (Invitrogen). The PCR reaction conditions were: 95 °C for 4 min; 35 cycles of 95 °C for 1 min, 57 °C for 1 min, 72 °C for 1 min; and a final extension step at 72 °C for 4 min [50].

The primers used to detect the V410L mutation were: V410fw (5′-GCG GGC AGG GCG GCG GGG GCG GGG CCA TCT TCT TGG GTT CGT TCT ACC GTG-3′), L410fw (5′-GCG GGC ATC TTC TTG GGT TCG TTC TAC CAT T-3′); and 410R (5′-TTC TTC CTC GGC GGC CTC TT-3′). Amplification was performed following the methodology described by Villanueva-Segura et al. [36]. In a 1.5 ml tube, a master mix was prepared as follows: 12.50 μl of GoTaq (Promega, Madison, WI, USA) was mixed with 11.35 μl of nuclease-free H2O (NFW, Promega), and 1 μM of each of the primers (V410fw, L410fw and 410rev) was then added. The master mix (24 μl) was then placed in 0.2 ml stoppered tubes followed by 1 μl of mosquito genomic DNA (~25 ng). The tubes were centrifuged for 1 min at 3300× rpm, then placed in a T100 thermocycler (Bio-Rad) with the following temperature program: 3 min at 95 °C; 30 cycles of 1 min at 95 °C, 20 s at 56 °C and 20 s at 72 °C; followed by 5 min at 72 °C.

DNA from the New Orleans susceptible strain was used as a negative control for the kdr assays, and previously genotyped individuals were used as positive controls.

The PCR products for the V410L, V1016I and F1534C assays were visualized on 2.5%, 3% and 4% agarose gels, respectively, using an UVITEC (Cambridge, UK) imaging system.

The frequencies of the alleles were determined for each population at each generation. We verified that the populations at the parental generation (FS0) were in Hardy-Weinberg equilibrium by means of Chi-square analysis. In addition, Wright’s FIS inbreeding coefficient was estimated, along with Wald’s correction [51, 52].

Biochemical assays

Thirty-two 2–4-day-old unfed female mosquitoes were individually homogenized in 2 ml of 0.01M phosphate buffer (pH 7.2). Aliquots of 100 μl of the homogenate were distributed in triplicate in flat-bottom microplates (Corning, Tewksbury, MA, USA) for testing the activity levels of the following families of enzymes: α- and β-esterases, MFOs, and GSTs, based on mosquito-specific biochemical assay protocols [53,54,55,56].

α- and β-esterases

To measure the activity levels of α- and β-esterases, 100 μl of either α- or β-naphthyl acetate (Sigma-Aldrich, St. Louis, MO,USA) dissolved in acetone (CTR Scientific, Monterrey, N.L., Mexico), and phosphate buffer (KPO4, pH 7.2) was added to each well and incubated for 20 min at room temperature. Subsequently, 100 μl of fast blue (tetrazotized O-dianisidine (Sigma-Aldrich) dissolved in distilled water) was added to each well and the microplates incubated for another 4 min, before absorbance was read at 540 nm using a spectrophotometer (ASYS Hitech GmbH, Eugendorf, Austria).

MFOs

200 μl of TMBZ (3,3,5,5-tetramethyl-benzidine dihydrochloride (Sigma-Aldrich) previously dissolved in methanol (Jalmek, Monterrey, N.L. Mexico) and 0.25 M sodium acetate buffer pH 5.0) were added to each well. Subsequently 25 μl of 3% hydrogen peroxide (H2O2) (Jalmek) was added to each well. After 10 min of incubation at room temperature the microplate was read at a wavelength of 620 nm.

GSTs

100 μl of reduced glutathione (GSH; Sigma-Aldrich) and 100 μl 1-chloro-2,4 dinitrobenzene (CDNB; Sigma-Aldrich) previously diluted in acetone and KPO4 buffer were added to each well. Absorbance was read at 340 nm immediately (T0) and after 10 min of incubation (T10). The values obtained after subtracting the initial reading (T0) from the reading at 10 min (T10) were used for the statistical analyses.

Total protein content was determined by the Bradford method [56], which was used to correct the activity values for all of the enzymes evaluated [57].

Positive and negative controls were included for MFOs and esterases. The same volume of homogenate which was used in the respective assays was used in running the controls. For the α- and β-esterases, α- and β -naphthyl acetate solution was used as positive controls, respectively. Cytochrome-C (Merck, Darmstadt, Germany) solution was the positive control for the MFO assay. KPO4 buffer was used as a negative control for each biochemical assay.

Calibration curves

The calibration curves were generated using control solutions for each enzyme family. α- or β-naphthol (Sigma-Aldrich) was used for α-esterases and β-esterases, with a concentration range of 0.3 to 5 μg/μl and 0.5 to 4 μg/μl, respectively. Cytochrome C concentrations from 0.002 to 0.65 μg/μl were used for MFOs, and 0.1 to 6.5 μg/μl BSA solutions were used for total protein determination.

The absorbance values obtained for each enzyme were used to calculate the mean absorbance values per mosquito, which was then converted to enzyme activity after accounting for the homogenization volume, total protein content and the activity unit considered for each enzyme. Each value was multiplied by a conversion factor that was obtained from the calibration curves [57, 58].

The mean activity values underwent an analysis of variance (ANOVA, P < 0.05) and Tukey’s multiple comparison of means (P < 0.05). Normality of the variance was verified by the Bartlett test. Statistical analyses were carried out with GraphPad Prism v.7 (GraphPad Software, Inc, Version 6.01, La Jolla, CA, USA; https://www.graphpad.com). Activity levels for all enzymes were determined in the parental population and in each deltamethrin-selected generation (FS0–FS5).

The activity values corresponding to the 99th percentile of the New Orleans reference strain were calculated for each enzyme family. The average enzymatic activity was classified as non-altered (NA) when < 15% of individuals did not exceed the 99th percentile of the reference strain, incipiently altered (IA) if 15–50% of the individuals exceeded the 99th percentile, and altered (A) when this percentage exceeded 50%, according to the criteria established by Montella et al. [59]. Additionally, the LC50 values and mean enzyme activity values of the parental generation and each selected generation (FS1–FS5) underwent linear regression analysis. Correlation coefficients were calculated to determine the degree of association between the two variables.

Results

Bioassays

Deltamethrin KC50 and LC50 were determined for the parental generation and in each selected generation (FS1-FS5) as well as for the insecticide-susceptible New Orleans (NO) reference strain. Although some variation was observed, KC50 and LC50 generally increased with respect to the parental generation (FS5vs FS0) (Tables 2, 3, Fig. 2).

Table 2 Deltamethrin knockdown concentrations (KC50 and KC90) and resistance ratios (RR) for Aedes aegypti females in the parental generation (FS0) and the deltamethrin-selected generations (FS1–FS5)

When comparing the LC50 values obtained from the parental generations to the NO strain, five of the populations (Merida, Progreso, Hoctun, Agua Dulce and Jose Cardel) were highly resistant to deltamethrin with an RRLC50 of 18–134-fold, and only one population (Hunucma) showed moderate resistance (RRLC50 of 6-fold) (Table 3).

Table 3 Deltamethrin knockdown concentrations (LC50 and LC90) and resistance ratios (RR) for Aedes aegypti females in the parental generation (FS0) and the deltamethrin-selected generations (FS1–FS5)
Fig. 2
figure2

a Knockdown concentrations values (KC50) as response to the selection with deltamethrin. b Lethal concentrations values (LC50) as response to the selection with deltamethrin in six populations of Aedes aegypti through five generations

With respect to knockdown in the parental generation, RRKC50 was high in two populations, Merida with RRKC50 of 20-fold and Jose Cardel with RRKC50 of 14-fold. Three populations (Hoctun, Progreso and Agua Dulce) showed moderate resistance with RRKC50 ranging between 5–9-fold, and one population (Hunucma) had a low RRKC50 of 3-fold (Table 2).

When the populations were selected with deltamethrin, the increase in LC50 was not significant for FS1 and FS2 for the Merida population (based on overlapping 95% CIs). However, there was a significant increase in LC50 from FS1 to FS5 in all other populations of Ae. aegypti (Table 3). For KC50, the increases were significant between all generations of selection (FS1-FS5) for the Progreso, Agua Dulce and Jose Cardel populations but were not significant for the Merida and Hunucma populations for FS1-FS2 and the Hoctun population for FS1-FS3 (Table 2).

When comparing the LC50 and KC50 values for the last deltamethrin-selected generation (FS5) vs the parental generation (FS0), we observed an increase of ~1.5–5.5 times for LC50 and ~1.4–2.7 times for KC50 in all populations (Fig. 2).

The results revealed that deltamethrin resistance increased in response to selection pressure (Tables 2, 3, Fig. 2). The knockdown resistance ratio (RRKC50) value increased from 20-fold to 38-fold from FS0-FS5 in the Merida population, from 6-fold to 23-fold for Progreso, from 3-fold to 8-fold for Hunucma, from 5-fold to 10-fold in Hoctun, and slight increases were observed in the Agua Dulce and Jose Cardel populations with 9–12-fold and 14–22-fold differences, respectively. Lethal concentrations also increased between FS0 and FS5, with the RRLC50 increasing in the Hunucma populations from 6-fold to 34-fold, doubling in the Agua Dulce and Progreso populations (25–49-fold and 18–36-fold, respectively) and with lower increases detected in the Merida, Hoctun and Jose Cardel populations (134–197-fold, 32–52-fold and 41–71-fold, respectively). When analyzing the relative increases in RRKC50 and RRLC50 during the selection, the Hunucma population showed the highest single-generation increase (> 0.3-fold) from FS1vs FS2 for KC50; Hunucma showed also the highest increase in RRLC50 (2.83) in the first selected generation relative to the parental generation (FS1vs FS0).

Realized heritability

The heritability of resistance (h2) estimated over the generations of deltamethrin selection was highest in the Hoctun and Merida populations, with values of 0.948 and 0.900, respectively, for LC50. In the case of KC50, the highest values were obtained for the Progreso population (0.779) and Merida populations (0.519). The response to the selection (R) was also high for Merida and Progreso, for both KC50 and LC50. The selection differential (S) was high in the Progreso population for LC50, but in the case of KC50, the highest values were obtained for Hoctun and Hunucma. The number of generations required for a 10-fold increase in LC50, (G), is the reciprocal of the response to selection (R) [47]. Thus, for the Progreso, Merida and Hoctun populations, an estimated mean of ~ 2 generations was necessary for a 10-fold increase in LC50; however, for Hunucma, Agua Dulce and Jose Cardel, a mean of > 4 generations was necessary. Regarding KC50, an estimated mean of ~ 2 generations was necessary for a 10-fold increase for the Merida and Progreso populations, but for Hunucma, Hoctun, Jose Cardel and Agua Dulce, means of 3, ~ 6, 7 and ~ 25 generations, respectively, were necessary for a 10-fold increase (Tables 4, 5).

Table 4 Estimation of heritability (h2) of resistance in Aedes aegypti after deltamethrin selection for knockdown concentration (KC50)
Table 5 Estimation of heritability (h2) of resistance in Aedes aegypti after deltamethrin selection for lethal concentration (LC50)

Molecular assays

Tables 6, 7, 8 summarize the frequency of L410, I1016 and C1534 alleles across all populations and generations (FS0–FS5). The three kdr mutations were present in all the parental populations, with frequencies ranging between 0.36–078 for the L410 allele and 0.34–0.77 for I1016. The highest frequencies were for C1534, ranging from 0.59–1. All populations were in Hardy-Weinberg equilibrium at the parental generation (FS0), with the exception of Jose Cardel for the mutation F1534C.

Table 6 V410L genotypes and allele frequencies in Aedes aegypti females in the parental and all selected generations. χ2 Hardy-Weinberg, inbreeding coefficients (FIS) and significance testing was calculated for the parental generation
Table 7 V1016I genotypes and allele frequencies in Aedes aegypti females in the parental and all selected generations. χ2 Hardy-Weinberg, inbreeding coefficients (FIS) and significance testing was calculated for the parental generation
Table 8 F1534C genotypes and allele frequencies in Aedes aegypti females in the parental and all selected generations. χ2 Hardy-Weinberg, inbreeding coefficients (FIS) and significance testing was calculated for the parental generation

A total of 20 combinations of tri-locus genotypes were detected among the populations in the parental generation (FS0). Figure 3 shows the frequency of each of the 20 tri-locus genotype combinations. Triple homozygote resistant genotype (LL410/II1016/CC1534) occurred in all parental (FS0) populations of Ae. aegypti with the highest frequency for the Merida population (0.63) and the lowest for Hoctun (0.06), and their frequency increased as generations were subsequently selected (Fig. 4).

Fig. 3
figure3

Frequencies of the tri-locus genotypes in the parental generation (FS0) of each population of Aedes aegypti. Genotypes order: 410/1016/1534

Fig. 4
figure4

Frequencies of the resistant trilocus genotype LL410/II1,016/CC1,534 across selection with deltamethrin (FS0-FS1) from six populations of Aedes aegypti from Mexico

Interestingly, the triple heterozygous genotype (VL410/VI1016/FC1534) was only found in the populations of Hunucma, Hoctu and Jose Cardel, in which the triple wild-type genotype (VV410/VI1016/FC1534) was also found.

In general, the frequencies of the three kdr alleles increased over the selected generations. The frequencies of L410 from the parental generation to the last selected generation increased from 0.78 to 0.92 in the Merida population, 0.65 to 0.92 in Progreso, 0.36 to 0.95 in Hunucma, 0.48 to 0.87 in Hoctun, 0.55 to 0.98 in Agua Dulce and 0.66 to 0.80 in Jose Cardel.

For the I1016 allele, increases of 0.77–0.92 were observed in the Merida population, 0.69–0.92 in Progreso, 0.45–0.95 in Hunucma, 0.34–0.88 in Hoctun, 0.56–0.96 in Agua Dulce and 0.63–0.92 in Jose Cardel.

The frequencies of the C1534 allele increased from 0.59 to 1 in the Hoctun population, 0.64 to 1 in Hunucma, 0.89 to 1 in Jose Cardel and 0.97 to 1 in Progreso; the C1534 allele was fixed at 1.0 in the Merida and Agua Dulce populations in all generations. This allele reached fixation after the first selection with deltamethrin (FS1) in the population from Progreso and Jose Cardel, in the third selected generation for Hoctun and in the fourth selected generation for Hunucma.

The L410 and I1016 alleles were selected in almost the same proportion from FS0 to the FS5 in all populations analyzed (Tables 6, 7, Fig. 5). Figure 6 shows that the increase in the frequencies of the L410 alleles in response to the selection with deltamethrin correlates significantly (P < 0.05) with the increase in the frequencies of the I106 allele by 98% in the population of Merida, 96% in Hunucma, 95% in Hoctun and 95% in Agua Dulce. Additionally, in the populations of Hoctun and Hunucma we found that the allelic frequency of C1534 increased significantly (P < 0.05) with the increase in allele frequency of I1016 and L410 throughout the selection in the populations of Hunucma and Hoctun (Fig. 7).

Fig. 5
figure5

Allele frequencies of L410 (a), I1016 (b) and C1534 (c) and their maximum and minimum frequencies in six Aedes aegypti populations across selection with deltamethrin (FS0-FS5)

Fig. 6
figure6

Regression analysis of allele frequencies of two para mutations: I1016 vs L410 across the selection with deltamethrin (FS0-FS5) in Aedes aegypti from Merida (a), Hunucma (b), Hoctun (c), Agua dulce (d)

Fig. 7
figure7

Regression analysis of allele frequencies of C1534 vs I016 (a, b) and C1534 vs L410 (b, c) across the selection with deltamethrin (FS0-FS5) in Aedes aegypti from Hunucma and Hoctun

The frequency of the mutant alleles for each mutation was correlated with the KC50 and LC50 values across the generations of selection (Table 9). KC50 values were significantly correlated with L410 frequencies in the Merida and Hoctun populations (P < 0.05); and with the frequencies of I1016 in the Merida, Hoctun and Jose Cardel (P < 0.05). Meanwhile, the frequencies of L410 and I1016 were significantly correlated with LC50 values in the Merida and Hunucma populations as well as C1534 correlated with LC50 values in Hunucma and Jose Cardel (P < 0.05). The presence and frequency of triple homozygous mutant individuals showed a significant correlation with KC50 only in the Merida population and LC50 in the Hunucma population (P < 0.05).

Table 9 Pearsonʼs correlation coefficients (r) (P-values in parenthesis) between the resistant allelic frequencies for L410, I1016, C1534 and with the frequencies of tri-locus genotypes vs KC50 and LC50 across the generations of selection with deltamethrin in five populations of Ae. aegypti from Mexico

Biochemical assays

Table 10 shows the mean enzyme activity levels detected in the biochemical assays in the parental generation and all deltamethrin-selected generations, where at least 30 individuals per population were analyzed. When comparing values from the parental generations with the 99th percentile of the New Orleans reference strain, five of the six populations analyzed (Merida, Progreso, Hunucma, Jose Cardel and Agua Dulce) showed altered enzyme activity for α-esterases (˃ 50% of the individuals exceeded the threshold of resistance), and Hoctun showed incipient altered enzyme activity (39% of the individuals exceeded the threshold). Merida, Agua Dulce and Jose Cardel showed altered enzyme activity for β-esterases (˃ 50%), and Progreso, Hunucma and Hoctun showed incipient enzyme alteration with percentages ranging between 16–43%. Altered activity of MFOs was detected in Progreso, Hoctun, Agua Dulce and José Cardel in 69–95% of the specimens analyzed. Only the Agua Dulce population showed altered enzyme activity for GSTs (71%) in the parental generation.

Table 10 Quantification of enzymatic activity in Aedes aegypti populations in the parental generation and deltamethrin-selected generations with respect to the New Orleans reference strain

The results show variations in enzyme activity across populations over the generations post-selection. However, when comparing FS5vs FS0, α-esterases remained altered in the Progreso, Hunucma and Agua Dulce; β-esterases in Progreso and Agua Dulce, and an increase in the mean activity levels of GSTs was observed only in the Hoctun and Agua Dulce populations (P < 0.0001) (Table 10). GSTs were the only enzyme family for which a significant correlation was detected between level of activity and LC50 values, but only in the FS4 generation (r = 0.83, P < 0.05).

Discussion

Target-site insensitivity and increased metabolic activity are key mechanisms of pyrethroid resistance in mosquito populations worldwide. The control of adult mosquitoes in Mexico has been carried out mainly with pyrethroids since 1999, and deltamethrin has been in use since the beginning of 2000 [9]. In the present study, five of the populations analyzed were found to be highly resistant to the pyrethroid deltamethrin in the parental generation according to the LC50 (RRLC50 > 10-fold), and only one population showed moderate resistance (RRCL50 5–10-fold). Likewise, high RRKC50 values were detected in two populations, three showed moderate resistance, and one showed low resistance. These results suggest that deltamethrin resistance is well-established in these populations.

In response to further deltamethrin selection, both LC50 and KC50 increased over the generations, with the concentrations capable of knocking down or killing 50% of the population at least doubling by FS5 compared to FS0. The highest heritability values (h2 ≥ 0.90) for LC50 were detected in the populations from Merida and Hoctun, followed by the Progreso population (0.61), suggesting high values of additive genetic variance. The highest value of h2 for KC50 was detected in the Progreso population (h2 > 0.70), followed by the Merida population (0.52). A high h2 value predicts a rapid response to selection, while a low value suggests a slow response to artificial selection [46].

The frequencies of the L410, I1016 and C1534 alleles increased with deltamethrin selection in all populations. When analyzing the frequencies in the last selected generation among the different populations, the frequencies of the L410 allele ranged between 0.80–0.98; this allele is able to significantly decrease the sensitivity of the sodium channel for both permethrin and deltamethrin [27]. The I1016 allele frequencies were 0.88–0.96 in FS5 and the C1534 allele reached fixation by FS5 in all populations analyzed. These results suggest that deltamethrin can select for C1534 more rapidly than I1016, corroborating the findings of Alvarez et al. [60].

The co-occurrence of I1016 and C1534 and high levels of pyrethroid resistance has been reported in other countries [61,62,63,64,65,66]. In Mexico, the co-occurrence of these two mutations and pyrethroid resistance has been also documented [34, 50]. Vera-Maloof et al. [34] originally suggested that in order to present resistance to pyrethroids, the sequential evolution of both mutations was necessary. They considered it unlikely that I1016 had evolved independently due to the low fitness exhibited by haplotype I1016/F1534; instead, they hypothesized that C1534 first evolved by conferring a low level of resistance individually and that I1016 arose from haplotype V1016/C1534 and selected quickly due to the high level of resistance conferred by the double mutant haplotype. The sequential selection of F1534C and V1016I was later confirm by Chen et al. [31].

The V410L mutation was reported for the first time in the sodium channels of mosquitoes by Haddi et al. [27]. They demonstrated that this mutation reduced the sensitivity of the vgsc gene to both type I (i.e. permethrin) and type II pyrethroids (i.e. deltamethrin). The presence of this mutation has been demonstrated in Ae. aegypti from Mexico for more than 16 years and it has been shown that L410 is in greater linkage disequilibrium with I1016 than with C1534. Our results show that L410 and I1016 responded to the selection with deltamethrin in a similar way in populations where C1534 was fixed; or in Hoctun and Hunucma where the increase in C1534 was associated with the increase in L410 and I1016. This is consistent with the proposed sequential model where both V410L and V1016I might have occurred independently on a C1534 haplotype and then became cis to C1534 by recombination. Or alternatively and considering that the F1534C mutation was fixed in practically four of the six basal populations (FS0), the three mutations arose independently at very low frequencies, and then by two recombination events, came to occur in a cis arrangement [30].

Increased activity of mixed-function oxidases (MFOs) and esterase activity have been associated with pyrethroid resistance in Ae. aegypti [38, 67]. Our results detected incipiently altered levels of MFO activity in four populations. Altered levels of α-esterase activity were detected in five populations, and at least four populations showed altered levels of β-esterase activity. However, there was no association between increased activities of these enzymes and phenotypic resistance, suggesting that these enzymes are not strongly associated with the metabolism of deltamethrin. However, contrary to the findings of Son et al. [68], who reported increased levels of MFOs in response to deltamethrin selection, we found that MFO activity decreased with successive selection. Given that the frequencies of the kdr mutations increased significantly with deltamethrin selection, the rapid selection of these mutations could have conditioned the activity of the enzymes in relation to deltamethrin resistance.

One criticism in our methodology is that we compare the enzymatic activity of each generation of deltamethrin in each of the populations with respect to the enzymatic activity of the susceptible New Orleans strain, when the most appropriate comparison should be with respect to the same population, in the same generation, without selection with deltamethrin.

In the metabolic resistance, three families of detoxifying enzymes, oxidases, esterases and glutathione transferases are mainly involved, all associated to confer resistance through overexpression, increased insecticide metabolism or a greater affinity for the chemical [6, 17]. An example of this contribution is reported by Lumjuan et al. [69] who determined an overexpression of GSTe2, GSTe5, GSTe7 and GSTE5-5 of epsilon glutathione transferase in individuals with resistance to DDT and pyrethroids. This fact is directly associated with the reduction in the susceptibility of these individuals to chemical control. Similarly, using quantitative trait loci (QTL) it was shown that temephos resistance is related to a QTL on chromosome two where a carboxylesterase cluster occurs; subsequent studies showed that there is an increase in the expression of these enzymes and therefore a possible intervention as a resistance mechanism [70]. Similar results were reported by Saavedra et al. [71], where the significant participation of an esterase marker discovered by QTL, CCEunk7o, and its relationship to resistance development in conjunction with mutations in the voltage-gated sodium channel was demonstrated.

Alongside gene expression studies, several investigations focus on studying the levels of detoxifying enzymatic activity. Proof of this is the increase in activity suffered by multiple-function oxidases, which includes cytochrome P450, in individuals selected with deltamethrin for 15 generations compared to the unselected susceptible group [60]; however, as for this mechanism to be in itself the one that confers insecticide resistance is questionable since other mechanisms were described for these populations.

To verify that the enzymatic activity has a direct effect on the decrease in susceptibility, synergists are used, chemical compounds that inhibit detoxifying enzymes. By inhibiting the enzymatic activity of multiple function oxidases, esterases and glutathione transferases directly involved as a resistance mechanism, proof of this is reported by Bharati & Saha [72] where an increase in enzymatic activities of multiple-function oxidases (CYP450) and carboxylesterases was identified as a mechanism of resistance, by inhibiting such enzymes, the susceptibility to the pyrethroid deltamethrin and the carbamate propoxur was recovered in populations of Ae. aegypti from Bangladesh. A similar pattern was demonstrated in populations of Ae. aegypti resistant to DDT and pyrethroids from Selangor, Malaysia, where the inhibition of multiple function oxidases and glutathione-S transferases through an assay with synergists correlated with an increase in the activity of such enzymes as a possible mechanism of resistance [73].

Identifying the direct and exclusive participation of metabolic resistance to insecticides is a complex process due to the simultaneous presence of mechanisms such as target site mutations and elevated levels of detoxifying enzymes. Despite this, the occurrence of Ae. albopictus populations free of insecticide-resistant target site mutations has been reported, demonstrating the direct and exclusive participation of cytochrome P450 genes such as CYP6P12 and the overexpression of cuticular genes as the main resistance mechanisms [74].

Conclusions

The high levels of resistance associated with high frequencies of kdr alleles L410, I1016 and C1534 obtained through artificial selection suggest the important role of these mutations in resistance to deltamethrin. The role of metabolic resistance was less clear, as the activity levels of key enzyme groups appeared to increase or decrease without clearly relating to the selection with deltamethrin. These findings highlight the importance of monitoring the susceptibility status of mosquito populations before choosing insecticide products for vector control, especially in areas where resistance may already be present and further selection could quickly result in fixation of kdr mutations.

Availability of data and materials

All data generated or analyzed during this study are included in the published article.

Abbreviations

kdr :

knockdown resistance

KC50 :

50% knockdown concentration

LC50 :

50% lethal concentration

RR:

resistance ratio

vgsc :

voltage gated sodium channel

CDC:

Centers for Disease Control and Prevention

WHO:

World Health Organization

DDT:

dichlorodiphenyltrichloroethane

PCR:

polymerase chain reaction

MFOs:

mixed-function oxidases

GSTs:

glutathione S-transferases

References

  1. 1.

    Smith L, Kasai S, Scott JG. Pyrethroid resistance in Aedes aegypti and Aedes albopictus: important mosquito vectors of human diseases. Pestic Biochem Physiol. 2016;133:1–12.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  2. 2.

    Huber JH, Childs ML, Caldwell JM, Mordecai EA. Seasonal temperature variation influences climate suitability for dengue, chikungunya, and Zika transmission. Plos Neglect Trop Dis. 2018;12:e0006451.

    Article  Google Scholar 

  3. 3.

    Kandel Y, Vulcan J, Rodriguez SD, Moore E, Chung HN, Mitra S, et al. Widespread insecticide resistance in Aedes aegypti L. from New Mexico, USA. PLoS One. 2019;14:e0212693.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  4. 4.

    WHO. WHO recommended insecticides for indoor residual spraying against malaria vectors. Geneva: World Health Organization; 2018. https://www.who.int/neglected_diseases/vector_ecology/vector-control/Insecticides_IRS_22_September_2018.pdf?ua=1. Accessed 19 Feb 2020.

  5. 5.

    WHO. Global strategy for dengue prevention and control 2012–2020. Geneva: World Health Organization; 2012. https://www.afro.who.int/publications/global-strategy-dengue-prevention-and-control-2012-2020. Accessed 19 Feb 2020.

  6. 6.

    Moyes C, Vontas J, Martins A, Ng L, Koou SY, Dusfour I. Contemporary status of insecticide resistance in the major Aedes vectors of arboviruses infecting humans. PLoS Neglect Trop Dis. 2017;11:e0005625.

    Article  CAS  Google Scholar 

  7. 7.

    Amelia-Yap ZH, Chen CD, Sofian-Azirun MM, Low VL. Pyrethroid resistance in the dengue vector Aedes aegypti in Southeast Asia: present situation and prospects for management. Parasit Vectors. 2018;11:332.

    PubMed  PubMed Central  Article  Google Scholar 

  8. 8.

    Ray DE, Fry JR. A reassessment of the neurotoxicity of pyrethroid insecticides. Pharmacol Therapeut. 2006;111:174–93.

    CAS  Article  Google Scholar 

  9. 9.

    DOF (Diario Oficial de la Federacion). NOM-032-SSA-2-2002 para la vigilancia epidemiológica, prevención y control de enfermedades transmitidas por vectores. Mexico. 2003. http://www.salud.gob.mx/unidades/cdi/nom/032ssa202.html. Accessed 19 Feb 2020.

  10. 10.

    CENAPRECE. Lista de Insumos Recomendados por el CENAPRECE para el combate de insectos vectores 2019. Mexico. 2019. https://www.gob.mx/cms/uploads/attachment/file/469289/Lista_de_Insumos_Recomendados_por_el_CENAPRECE.pdf. Accessed 19 Feb 2020.

  11. 11.

    Siller Q, Ponce G, Lozano S, Flores AE. Update on the frequency of Ile1016 mutation in voltage gated sodium channel gene of Aedes aegypti in Mexico. J Am Mosq Control Assoc. 2011;27:357–62.

    PubMed  Article  PubMed Central  Google Scholar 

  12. 12.

    Aponte HA, Penilla-Navarro P, Dzul-Manzanilla F, Che-Mendoza A, López AD, Solís F, et al. The pyrethroid resistance status and mechanisms in Aedes aegypti from Guerrero state, Mexico. Pest Biochem Phys. 2013;107:226–34.

    CAS  Article  Google Scholar 

  13. 13.

    Flores AE, Ponce-García G, Lopez-Monroy B, Villanueva-Segura OK, Rodriguez-Sánchez IP, Arredondo-Jímenez JI, Manrique-Saide P. Current status of the insecticide resistance in Aedes aegypti (Diptera: Culicidae) from Mexico. In: Stanislav T, editor. insecticide resistance. London: IntechOpen; 2016. p. 99–109.

    Google Scholar 

  14. 14.

    Flores AE, Grajales JS, Fernandez I, Ponce G, Loaiza H, Lozano S, et al. Mechanisms of insecticide resistance in field populations of Aedes aegypti (L.) from Quintana Roo, Southern Mexico. J Am Mosq Control. 2006;22:672–7.

    CAS  Article  Google Scholar 

  15. 15.

    Stevenson BJ, Pignatelli P, Nikou D, Paine MJ. Pinpointing P450s Associated with pyrethroid metabolism in the dengue vector, Aedes aegypti: developing new tools to combat insecticide resistance. PLoS Neglect Trop Dis. 2012;6:e1595.

    CAS  Article  Google Scholar 

  16. 16.

    Brogdon WG, McAllister JC. Insecticide resistance and vector control. Emerg Infect Dis. 1998;4:605–13.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  17. 17.

    Liu N. Insecticide resistance in mosquitoes: impact, mechanisms, and research directions. Ann Rev Entomol. 2015;60:537–59.

    CAS  Article  Google Scholar 

  18. 18.

    Narashashi T. Neuronal ion channels as the target sites of insecticides. Pharmacol Toxicol. 1996;79:1–14.

    Article  Google Scholar 

  19. 19.

    Rinkevich FD, Hedtke SM, Leichter CA, Harris SA, Su C, Brady SG, et al. Multiple origins of kdr-type resistance in the house fly, Musca domestica. PLoS One. 2012;7:e52761.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  20. 20.

    Du Y, Nomura Y, Satar G, Hu Z, Nauen R, He SY, et al. Molecular evidence for dual pyrethroid-receptor sites on a mosquito sodium channel. Proc Natl Acad Sci USA. 2013;110:11785–90.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  21. 21.

    Brengues C, Hawkes NJ, Chandre F, McCarroll L, Duchon S, et al. Pyrethroid and DDT cross-resistance in Aedes aegypti is correlated with novel mutations in the voltage-gated sodium channel gene. Med Vet Entomol. 2003;17:87–94.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  22. 22.

    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:785–98.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  23. 23.

    Chang C, Shen WK, Wang TT, Lin YH, Hsu EL, et al. A novel amino acid substitution in a voltage-gated sodium channel is associated with knockdown resistance to permethrin in Aedes aegypti. Insect Mol Biol. 2009;39:272–8.

    CAS  Article  Google Scholar 

  24. 24.

    Srisawat R, Komalamisra N, Eshita Y, Zheng M, Ono K, et al. Point mutations in domain II of the voltage-gated sodium channel gene in deltamethrin-resistant Aedes aegypti (Diptera: Culicidae). Appl Entomol Zool. 2010;45:275–82.

    CAS  Article  Google Scholar 

  25. 25.

    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:501–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  26. 26.

    Kushwah RBS, Dykes CL, Kapoor N, Adak T, Singh OP. Pyrethroid-resistance and presence of two knockdown resistance (kdr) mutations, F1534C and a novel mutation T1520I, in Indian Aedes aegypti. PLoS Neglect Trop Dis. 2015;9:e3332.

    Article  CAS  Google Scholar 

  27. 27.

    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.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  28. 28.

    Granada Y, Mejia-Jaramillo A, Strode C, Triana-Chavez O. A point mutation V419L in the sodium channel gene from natural populations of Aedes aegypti is involved in resistance to λ-cyhalothrin in Colombia. Insects. 2018;9:E23.

    PubMed  Article  PubMed Central  Google Scholar 

  29. 29.

    Kushwah RBS, Kaur T, Dykes CL, Kumar RH, Kapoor N, Singh OP. A new knockdown resistance (kdr) mutation F1534L in Aedes aegypti associated with insecticide resistance. BioRxiv. 2019;740829.

  30. 30.

    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:6747.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  31. 31.

    Chen M, Du Y, Wu S, Nomura Y, Zhu G, Zhorov BS, Dong K. Molecular evidence of sequential evolution of DDT- and pyrethroid-resistant sodium channel in Aedes aegypti. PLoS Neglect Trop Dis. 2019;13:e0007432.

    Article  CAS  Google Scholar 

  32. 32.

    Fernando SD, Hapugoda M, Perera R, Saavedra-Rodriguez K, Black WC, De Silva NK. First report of V1016G and S989P knockdown resistant (kdr) mutations in pyrethroid-resistant Sri Lankan Aedes aegypti mosquitoes. Parasit Vectors. 2018;11:526.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  33. 33.

    Hu Z, Du Y, Nomura Y, Dong K. A sodium channel mutation identified in Aedes aegypti selectively reduces cockroach sodium channel sensitivity to type I, but not type II pyrethroids. Insect Biochem Mol Biol. 2011;41:9–13.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  34. 34.

    Vera-Maloof FZ, Saavedra-Rodriguez K, Elizondo-Quiroga AE, Lozano-Fuentes S, Black WC IV. Coevolution of the Ile 1,016 and Cys1,534 mutations in the voltage gated sodium channel gene of Aedes aegypti in Mexico. PLoS Neglect Trop D. 2015;9:e0004263.

    Article  CAS  Google Scholar 

  35. 35.

    Ponce-Garcia G, Flores AE, Fernandez-Salas I, Saavedra-Rodriguez K, Reyes-Solis G, Lozano-Fuentes S. Recent rapid rise of a permethrin knock down resistance allele in Aedes aegypti in Mexico. PLoS Neglect Trop Dis. 2009;3:e531.

    Article  CAS  Google Scholar 

  36. 36.

    Villanueva-Segura OK, Ontiveros-Zapata K, Lopez-Monroy B, Ponce-Garcia G, Gutierrez-Rodriguez SM, Davila-Barboza JA, et al. Distribution and frequency of the kdr mutation V410L in natural populations of Aedes aegypti (L.) (Diptera: Culicidae) from eastern and southern Mexico. J Med Entomol. 2020;57:218–23.

    PubMed  Article  PubMed Central  Google Scholar 

  37. 37.

    Hemingway J, Ranson H. Insecticide resistance in insect vectors of human disease. Annu Rev Entomol. 2000;45:371–91.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  38. 38.

    Flores AE, Abeldaño W, Fernandez I, Badii MH, Loaiza H, Ponce G, et al. Elevated α-esterase levels associated with permethrin tolerance in Aedes aegypti (L.) from Baja California, México. Pestic Biochem Phys. 2005;82:66–78.

    CAS  Article  Google Scholar 

  39. 39.

    Bisset JA, Marin R, Rodriguez MM, Severson DW, Ricardo Y, French L, Díaz M, Perez O. Insecticide resistance in two Aedes aegypti (Diptera: Culicidae) strains from Costa Rica. J Med Entomol. 2013;50:352–61.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  40. 40.

    Flores AE, Reyes G, Fernández I, Sánchez FJ, Ponce G. Resistance to permethrin in Aedes aegypti (L.) in northern Mexico. Southwest Entomol. 2009;34:167–77.

    Article  Google Scholar 

  41. 41.

    Lopez B, Ponce G, Gonzalez JA, Gutierrez SM, Villanueva OK, Gonzalez G, et al. Susceptibility to chlorpyrifos in pyrethroid-resistant populations of Aedes aegypti (Diptera: Culicidae) from Mexico. J Med Entomol. 2014;51:644–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  42. 42.

    Brogdon W, Chan A. Guidelines for evaluating insecticide resistance in vectors using the CDC bottle bioassay. CDC technical report. Methods in Anopheles research, 2nd ed. Atlanta: Centers for Disease Control and Prevention; 2010.

  43. 43.

    Lozano FS, Saavedra RK, Black W IV, Eisen L. QCAL: a software application for the calculation of dose-response curves in insecticide resistance bioassays. J Am Mosq Control. 2012;28:59–91.

    Article  Google Scholar 

  44. 44.

    Abbott WS. A method of computing the effectiveness of an insecticide. Econ Entomol. 1925;18:265–7.

    CAS  Article  Google Scholar 

  45. 45.

    Mazzarri M, Georghiou G. Characterization of resistance to organophosphate, carbamate, and pyrethroid insecticides in field population of Aedes aegypti from Venezuela. J Am Mosq Control. 1995;11:315–22.

    CAS  Google Scholar 

  46. 46.

    Saavedra- Rodriguez K, Flores AE, Fernandez I, Strode C, Ranson H, Hemingway J, Black WC IV. Transcription of detoxification genes after permethrin selection in the mosquito Aedes aegypti. Insect Mol Biol. 2012;21:61–77.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  47. 47.

    Tabashnik BE. Resistance risk assessment: realized heritability of resistance to Bacillus thuringiensis in diamondback moth (Lepidoptera: Plutellidae), tobacco budworm (Lepidoptera: Noctuidae), and Colorado potato beetle (Coleoptera: Chrysomelidae). J Econ Entomol. 1992;85:1551–9.

    Article  Google Scholar 

  48. 48.

    Falconer DS. Introduction to quantitative genetics. 3rd ed. Harlow, Essex, UK/New York: Longmans Green/John Wiley & Sons; 1989.

    Google Scholar 

  49. 49.

    Coen E, Strachan T, Dover G. Dynamics of concerted evolution of ribosomal DNA and histone gene families in the melanogaster species subgroup of Drosophila. J Mol Biol. 1982;158:17–35.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  50. 50.

    Lopez-Monroy B, Gutierrez-Rodriguez SM, Villanueva-Segura OK, Ponce-Garcia G, Morales-Forcada F, Alvarez LC, Flores AE. Frequency and intensity of pyrethroid resistance through the CDC bottle bioassay and their association with the frequency of kdr mutations in Aedes aegypti (Diptera: Culicidae) from Mexico. Pest Manag Sci. 2018;74:2176–84.

    CAS  Article  Google Scholar 

  51. 51.

    Wright S. Systems of mating. II. The effects of inbreeding on the genetic composition of a population. Genetics. 1921;6:124–43.

    CAS  PubMed  PubMed Central  Google Scholar 

  52. 52.

    Agresti A, Coull BA. Approximate is better than “exact” for interval estimation of binominal proportions. Am Stat. 1998;52:119–26.

    Google Scholar 

  53. 53.

    Brogdon WG. Biochemical resistance detection: an alternative to bioassay. Parasitol Today. 1989;5:56–60.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  54. 54.

    Brogdon WG, Barber AM. Microplate assay of glutathione S-transferase activity for resistance detection in single mosquito triturates. Comp Biochem Physiol. 1990;96:339–42.

    CAS  Google Scholar 

  55. 55.

    Brogdon WG, McAllister JC, Vulule J. Heme peroxidase activity measured in single mosquitoes identifies individuals expressing an elevated oxidase for insecticide resistance. J Am Mosq Control Assoc. 1997;13:233–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  56. 56.

    Bradford M. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem. 1976;72:248–54.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  57. 57.

    Valle D, Montella IR, Ribeiro RA, Viana-Medeiros PF, Martins-Júnior AJ, Lima JBP. Metodologia para quantificação de atividade de enzimas relacionadas com a resistência a inseticidas em Aedes aegypti/Quantification methodology for enzyme activity related to insecticide resistance in Aedes aegypti. Rio de Janeiro, Brasília: Fundação Oswaldo Cruz, Secretaria de Vigilância em Saúde/Ministério da Saúde; 2006.

    Google Scholar 

  58. 58.

    Hemingway J. Techniques to detect insecticide resistance mechanisms (field and laboratory manual) WHO/CDC/CPC/MAL/98.6. Geneva: World Health Organization; 1998.

  59. 59.

    Montella I, Martins A, Fernandez V, Pereira J, Braga I, Valle D. Insecticide resistance mechanisms of Brazilian Aedes aegypti populations from 2001 to 2004. Am J Trop Med Hyg. 2007;77:467–77.

    PubMed  Article  PubMed Central  Google Scholar 

  60. 60.

    Alvarez LC, Ponce G, Saavedra-Rodriguez K, Lopez B, Flores AE. Frequency of V1016I and F1534C mutations in the voltage-gated sodium channel gene in Aedes aegypti in Venezuela. Pest Manag Sci. 2015;71:863–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  61. 61.

    Kawada H, Oo SZ, Thaung S, Kawashima E, Maung YN, Thu HM, et al. Co-occurrence of point mutations in the voltage-gated sodium channel of pyrethroid-resistant Aedes aegypti populations in Myanmar. PLoS Neglect Trop Dis. 2014;8:e3032.

    Article  CAS  Google Scholar 

  62. 62.

    Chapadense FG, Fernandes EK, Lima JB, Martins AJ, Silva LC, Rocha WT, et al. Phenotypic and genotypic profile of pyrethroid resistance in populations of the mosquito Aedes aegypti from Goiânia, Central West Brazil. Rev Soc Bras Med Trop. 2015;48:607–9.

    PubMed  Article  PubMed Central  Google Scholar 

  63. 63.

    Dusfour I, Zorrilla P, Guidez A, Issaly J, Girod R. Deltamethrin resistance mechanisms in Aedes aegypti populations from three French overseas territories worldwide. PLoS Neglect Trop Dis. 2015;9:e0004226.

    Article  CAS  Google Scholar 

  64. 64.

    Maestre-Serrano R, Pareja-Loaiza P, Gomez-Camargo D, Ponce-García G, Flores AE. Co-occurrence of V1016I and F1534C mutations in the voltage-gated sodium channel and resistance to pyrethroids in Aedes aegypti (L.) from the Colombian Caribbean region. Pest Manag Sci. 2019;75:1681–8.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  65. 65.

    Yaméogo F, Sanon A, McCall PJ, Shirozu T, Weetman D, et al. High frequencies of F1534C and V1016I kdr mutations and association with pyrethroid resistance in Aedes aegypti from Somgandé (Ouagadougou), Burkina Faso. Trop Med Health. 2019;47:4–11.

    Article  Google Scholar 

  66. 66.

    Deming R, Manrique-Saide P, Barreiro AM, Cardeña EU, Che-Mendoza A, Jones B, et al. Spatial variation of insecticide resistance in the dengue vector Aedes aegypti presents unique vector control challenges. Parasit Vectors. 2016;9:67.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  67. 67.

    Polson KA, Brogdon WG, Rawlins SC, Chadee DD. Characterization of insecticide resistance in Trinidadian strains of Aedes aegypti mosquitoes. Acta Trop. 2011;117:31–8.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  68. 68.

    Son-Un P, Choovattanapakorn N, Saingamsook J, Yanola J, Lumjuan N, et al. Effect of relaxation of deltamethrin pressure on metabolic resistance in a pyrethroid-resistant Aedes aegypti (Diptera: Culicidae) strain harboring fixed P989P and G1016G kdr alleles. J Med Entomol. 2018;55:975–81.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  69. 69.

    Lumjuan N, Rajatileka S, Changsom D, Wicheer J, Leelapat P, 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:203–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  70. 70.

    Paiva MHS, Lovin DD, Mori A, Melo-Santos MAV, Severson DW, et al. Identification of a major quantitative trait locus determining resistance to the organophosphate temephos in the dengue vector mosquito Aedes aegypti. Genomics. 2016;107:40–8.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  71. 71.

    Saavedra-Rodriguez K, Strode C, Flores Suarez A, Fernandez Salas I, Ranson H, et al. Quantitative trait loci mapping of genome regions controlling permethrin resistance in the mosquito Aedes aegypti. Genetics. 2008;180:1137–52.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  72. 72.

    Bharati M, Saha D. Multiple insecticide resistance mechanisms in primary dengue vector, Aedes aegypti (Linn.) from dengue endemic districts of sub-Himalayan West Bengal, India. PLoS One. 2018;13:e0203207.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  73. 73.

    Leong CS, Vythilingam I, Liew JWK, Wong ML, Wan-Yusof WS, et al. Enzymatic and molecular characterization of insecticide resistance mechanisms in field populations of Aedes aegypti from Selangor, Malaysia. Parasit Vectors. 2019;12:236.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  74. 74.

    Ishak IH, Riveron JM, Ibrahim SS, Stott R, Longbottom J, et al. The cytochrome P450 gene CYP6P12 confers pyrethroid resistance in kdr-free Malaysian populations of the dengue vector Aedes albopictus. Sci Rep. 2016;6:24707.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

Download references

Acknowledgements

We would like to thank Kevin Ontiveros-Zapata (Universidad Autonoma de Nuevo Leon, Mexico) for processing mosquitoes in the laboratory.

Disclaimer

The findings and conclusions in this paper are those of the authors and do not necessarily represent the official position of the Centers for Disease Control and Prevention.

Funding

This study was financially supported by CONACYT Mexico through the fund Problemas Nacionales, Project no. PN2016-2134.

Author information

Affiliations

Authors

Contributions

AEF conceived the study and obtained funding. YCP, PMS and GPG carried out collections of the populations included in the study. YCP performed the laboratory work. AEF, YCP, AL, GPG, KVS, BLM and IPRS analyzed the data. AEF, YCP and AL wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Adriana E. Flores.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Contreras-Perera, Y., Ponce-Garcia, G., Villanueva-Segura, K. et al. Impact of deltamethrin selection on kdr mutations and insecticide detoxifying enzymes in Aedes aegypti from Mexico. Parasites Vectors 13, 224 (2020). https://doi.org/10.1186/s13071-020-04093-3

Download citation

Keywords

  • Aedes aegypti
  • Heritability
  • kdr mutations
  • V410L
  • V1016I
  • F1534C
  • Metabolic detoxification