Next Article in Journal
Evaluation of Rearing Factors Affecting Clanis bilineata tsingtauica Mell Larvae Fed by Susceptible Soybean Variety NN89-29 in Spring and Autumn Sowing
Next Article in Special Issue
Alpha-Cypermethrin Resistance in Musca domestica: Resistance Instability, Realized Heritability, Risk Assessment, and Insecticide Cross-Resistance
Previous Article in Journal
Labial Sensory Organs of Two Leptoglossus Species (Hemiptera: Coreidae): Their Morphology and Supposed Function
Previous Article in Special Issue
Analysis of the Feeding Behavior and Life Table of Nilaparvata lugens and Sogatella furcifera (Hemiptera: Delphacidae) under Sublethal Concentrations of Imidacloprid and Sulfoxaflor
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial Distribution of Pyrethroid Resistance and kdr Mutations in Aedes aegypti from La Guajira, Colombia

by
Ronald Maestre-Serrano
1,
Zulibeth Flórez-Rivadeneira
1,2,
Juan M. Castro-Camacho
1,
Eva Soto-Arenilla
1,
Doris Gómez-Camargo
3,
Paula Pareja-Loaiza
4,
Gustavo Ponce-Garcia
5,†,
Alan E. Juache-Villagrana
5 and
Adriana E. Flores
5,*
1
Facultad de Ciencias de la Salud, Universidad Libre Seccional Barranquilla, Km 7 Antigua Via Puerto Colombia, Barranquilla 080001, Colombia
2
Secretaria de Salud Departamental, Gobernacion de La Guajira, Calle 12 # 8-19, Riohacha 440001, Colombia
3
Facultad de Medicina—Sede Zaragocilla, Universidad de Cartagena, Calle 30 N° 48-152, Cartagena de Indias 130001, Colombia
4
Facultad de Ciencias de la Salud, Universidad Simon Bolivar, Carrera 59 No. 59-92, Barranquilla 080002, Colombia
5
Facultad de Ciencias Biologicas, Universidad Autonoma de Nuevo Leon, Av. Universidad s/n Cd. Universitaria, San Nicolas de los Garza, NL 66455, Mexico
*
Author to whom correspondence should be addressed.
Deceased.
Insects 2023, 14(1), 31; https://doi.org/10.3390/insects14010031
Submission received: 14 December 2022 / Revised: 24 December 2022 / Accepted: 26 December 2022 / Published: 29 December 2022

Abstract

:

Simple Summary

Current strategies to suppress arbovirus outbreaks include insecticide use against larvae and adult mosquitoes. The control of Aedes aegypti by insecticides is challenging due to a rapid increase in resistance. In Colombia, pyrethroids have been one of the most widely used insecticides to control adult forms of Ae. aegypti because of its low impact on the environment, low toxicity to mammals, and greater effectiveness. We detected the frequency and intensity of resistance to permethrin, deltamethrin, lambda-cyhalothrin, and associated kdr mutations in Ae. aegypti from La Guajira, Colombia. Modeling the spatial distribution of pyrethroid resistance and the mechanisms (e.g., kdr mutations) provoking it could enhance current national and departmental vector control programs by classifying areas according to the insecticide resistance status of Ae. aegypti populations and orient strategies such as rotations in endangered areas. Thus, benefits such as susceptibility recovery might be obtained.

Abstract

Dengue, chikungunya, and Zika are of great concern to the public health of Colombia. One of the main control strategies for these diseases is the application of insecticides directed at the Aedes aegypti vector. However, insecticide resistance has been increasingly recorded in the country, making control measures difficult. Here, we evaluated the resistance profiles for pyrethroids in populations of Ae. aegypti from La Guajira, Colombia. The frequency (diagnostic dose, DD) and intensity (2×, 5×, and 10× DD) of resistance to permethrin, deltamethrin, and lambda-cyhalothrin were determined in 15 populations of Ae. aegypti from La Guajira, Colombia, using the bottle bioassay. The kdr mutations V1016I, F1534C, and V410L, were identified, and their allele and genotype frequencies were calculated. Finally, the mortality values for the analyzed pyrethroids were interpolated following the IDW method for predicting pyrethroid resistance. The populations of Ae. aegypti showed a high frequency of resistance to permethrin with a low to moderate intensity, which was associated with the triple-resistant haplotype LL410/II1016/CC1534. They remain susceptible to deltamethrin and, in some populations, expressed the risk of developing resistance to lambda-cyhalothrin.

Graphical Abstract

1. Introduction

Aedes aegypti (Diptera: Culicidae) is a species of interest to public health as it is a competent vector of arboviral diseases such as dengue, Zika, yellow fever, and chikungunya, causing morbidity and mortality impacts worldwide [1]. In the Americas between the year 2020 and the epidemiological week 22 of 2022, 4,131,433 cases of arbovirosis were recorded, of which 3,883,586 corresponded to dengue cases, 221,735 cases of chikungunya and 26,112 cases of Zika [2,3,4,5]. Colombia is one of the countries that contributes the most to the incidence of these arboviral diseases in the Americas; between 1995 and 2019, 1,371,722 cases of dengue were recorded in the country, of which 252,623 cases (18.4%) were in the Caribbean region [6,7,8,9,10], with 4910 cases (1.9%) being reported in the department of La Guajira, which in the last decade has shown an upward trend in the incidence of this disease [11].
On the other hand, the first case of chikungunya reported in Colombia was identified in week 37 of 2014 [12]. From then until 2019, 489,528 cases of this disease have been reported in the country, with 119,104 cases (24.3%) occurring in the Caribbean region and 11,364 (9.5%) in La Guajira. The first outbreak of Zika in the country was reported in October 2015. From then until the epidemiological close of 2019, 109,999 cases were recorded in the country, of which 19,662 cases (17.8%) corresponded to the Caribbean region and 706 (3.5%) of these were for the department of La Guajira [13].
In the absence of a vaccine available for these arboviruses, control measures have been aimed at Ae. aegypti, such as physical control to eliminate domestic breeding sites such as pools, tires, tanks, and pots, among others; biological control through fish, and bacteria, such as the species Bacillus thuringiensis var israelensis; growth regulators and chemical control with the application of insecticides [14].
In Colombia, organophosphate, organochlorine, carbamate, and pyrethroid insecticides have been used to control Ae. aegypti for more than five decades as the main tool for the national and departmental programs for the prevention and control of vector-transmitted diseases (VTD) [15]. However, resistance to these insecticides has been increasingly recorded, which has made control actions difficult within VTD programs in different regions of the country [16]. In Colombia, pyrethroids have been one of the most widely used insecticides to control adult forms of Ae. aegypti, because of its low impact on the environment, low toxicity to mammals, and greater effectiveness. In La Guajira, through the Vector Control Program, malathion and temephos have been used to control dengue, chikungunya, and Zika outbreaks; however, as of 2018, alpha-cypermethrin, pirimiphos-methyl, and lambda-cyhalothrin have been used [17].
The target site of action of pyrethroids is the voltage-gated sodium channel (VGSC): a transmembrane protein made up of four homologous domains (I–IV), each with six hydrophobic segments (S1–S6) [18]. Conformational changes in the structure of this channel have been identified by mutations that confer resistance to Ae. aegypti to pyrethroids, and these are called kdr (knockdown resistance) mutations causing knockdown resistance [19]. In the Caribbean region of Colombia, resistance to the pyrethroids lambda-cyhalothrin, deltamethrin, cyfluthrin, and permethrin has been found in Ae. aegypti populations with the identification of the kdr mutations V1016I, F1534C, and V410L [20,21,22,23,24,25]. In La Guajira, few works have evaluated the susceptibility status of vector populations. However, resistance to permethrin, deltamethrin, cyfluthrin, and lambda-cyhalothrin has been identified in the municipality of San Juan del Cesar, whereas resistance to lambda-cyhalothrin has been identified in the district of Riohacha [21,23]. These findings have allowed the interdisciplinary team of the vector program to develop new strategies that consist of rotating insecticide molecules to control these arboviral diseases. Despite these results, the susceptibility of Ae. aegypti to pyrethroids and their underlying resistance mechanisms are unknown for most of the municipalities in this department. All the above information is essential, observing the endemic behavior of dengue and the introduction of chikungunya and Zika viruses in recent years in La Guajira. Likewise, it is important to consider that the susceptibility of insects to pyrethroids varies in space and time between populations according to the selection pressure exerted by these insecticides between the different municipalities.
The main objective of this study was to evaluate the resistance status of Ae. aegypti to the pyrethroids permethrin, deltamethrin, and lambda-cyhalothrin and to assess the presence and contribution of the kdr mutations V410L, V1016I, and F1534C to pyrethroid resistance in the 15 municipalities of the department of La Guajira, Colombia.

2. Materials and Methods

2.1. Study Area

The study was carried out during 2020 and 2021 in the 15 municipalities of the department of La Guajira (Figure 1). La Guajira is bordered on the north by the Caribbean Sea, on the east by the Caribbean Sea and Venezuela, on the south by the department of Cesar, and on the west by the department of Magdalena and the Caribbean Sea. It has an area of 20,848 km2, representing 1.8% of the Colombian territory.
The inclusion criteria for the selection of the neighborhoods where the parental strains of Ae. aegypti were collected in the selected municipalities were high indices of Ae. aegypti, the frequent application of insecticides for the control of this mosquito, and a high incidence of dengue in the last ten years prior to the start of the project.

2.2. Collection of Mosquitoes and Obtaining F1/F2

Entomological inspections were carried out for the collection of immature forms of Ae. aegypti in potential breeding sites such as swimming pools, plastic/metallic bins, tires, and pots. The collected entomological material was transported in 500-mL plastic bottles to the insectarium of the Barranquilla Campus of Universidad Libre, where the F1 and F2 generations were obtained under controlled conditions of temperatures (28 ± 2 °C), relative humidity (60 ± 10%), and 12 h:12 h (light:dark) photoperiod.

2.3. Bioassays

The pyrethroids permethrin, deltamethrin, and lambda-cyhalothrin were evaluated using the stock solutions of the technical grade insecticides (ChemService, West Chester, PA, USA) using acetone as a diluent. The diagnostic doses (DD) of permethrin (15 µg/bottle), deltamethrin (10 µg/bottle), and lambda-cyhalothrin (10 µg/bottle) were used, considering the diagnostic time (DT) of 30 min for the three cases [26]. The bottle bioassay [27] used a 250-mL Schott bottle containing 1 mL of an acetone solution of technical-grade insecticide at the DD of each insecticide. The bottle was capped and shaken to ensure uniform coverage and then allowed to stand and dry for 24 h at room temperature with protection from light.
Between 15 and 25 unfed 3-day-old adult females were placed in the bottles until the diagnostic time was completed. The number of dead mosquitoes at the time of diagnosis was recorded to determine the level of susceptibility in the study populations. This procedure was repeated four times. A bottle with acetone, without insecticide, was used as a control bottle. In all cases, the correction was applied by Abbott’s formula when mortality in the control group was between 5 and 20%, and the bioassay was invalid when mortality exceeded 20% [28].
All the procedures described above were performed in all field populations of Ae. aegypti and the susceptible Rockefeller strain as a reference.
The frequency of resistance was calculated, and the populations were designated susceptible when mortality was ≥98%, while between 90 and 97% suggested the possibility of resistance that requires confirmation and mortality with <90% indicating resistance [29].
In the populations where resistance was determined at the diagnostic dose and time, the intensity of the resistance was analyzed through new bioassays that consisted of evaluating two (2×), five (5×), and ten (10×) times the DD [29]. The results were interpreted as follows, ≥98% mortality with 5× DD exposure was considered low-intensity resistance, and <98% mortality was considered moderate to high-intensity resistance. For the case of exposure to 10x the DD, a mortality ≥98% was considered moderate intensity resistance, while <98% mortality was considered high-intensity resistance.
A Kruskal–Wallis test followed by Dunn’s multiple comparisons test (α = 0.05) was used to compare the mortality rates produced by the DD of permethrin, deltamethrin, and lambda-cyhalothrin.

2.4. Identification of kdr Mutations

Forty mosquitoes of the F0 generation (field) were selected from each population and evaluated. DNA extraction was performed using the Extracta Kit (Quanta Biosciences 95091-250, Gaithersburg, MD, USA) and was later quantified using a UV-Vis NanoDrop One Microvolume spectrophotometer (Thermo Scientific; Waltham, MA, USA).
PCRs were performed in a CFX90 real-time thermal cycler (Bio-Rad; Hercules, CA, USA). The genotype for loci 1016, 1534, and 410 was determined by melting curve analysis. The amplification of the V1016I mutation in the VGSC gene was performed according to the method described by Saavedra et al. [30]; the final volume of each reaction was 20 µL, which included: 10 μL of Perfecta Sybr® Green Supermix (Quanta 95054-500, Beverly, MA, USA), 1 μL of both primers V1016(r) (5′-CGGGCAGGGCGGCGGGGGCGGGGCCACAAATTGTTTCCCACCCGCACCGG-3′), I1016(f) (5′-GCGGGCACAATTGTTTCCCACCCGCACTGA-3′) and I1016(r) (5′-GGATGAACCGAAATTGGACAAAAGC-3′), 6 µL of water, and 1 µL of the DNA template. Amplification reactions were performed according to the following thermal profile: an initial denaturation cycle at 95 °C for 3 min, 40 cycles at 95 °C for 10 s, 60 °C for 10 s, and 72 °C for 30 s. A final extension was performed at 95 °C for 10 s. The melting curves were determined by a denaturation gradient from 65 °C to 95 °C with an increase of 0.2 °C every 10 s.
The F1534C mutation was identified according to the method described by Yanola et al. [31]; the final volume of each reaction was 20 μL, which included 9 μL of Perfecta Sybr® Green Supermix (Quanta 95054-500, Beverly, MA, USA), 0.66 μL of the primers F1534(f) (5′-GCGGGCTCTACTTTGTGTTCTTCATCATATT-3′) and F1534(r) (5′-TCTGCTCGTTGAAGTTGTCGAT-3′), 0.65 μL of primer C1534(f) (5′-GCGGGCAGGGCGGCG GGGGCGGGGCCTCTACTTTGTGTTCTTCATCATGTG-3′), 7.5 µL of water, and 2 µL of the DNA template. Amplification reactions were performed with the following thermal profile: an initial denaturation cycle at 95 °C for 3 min, 37 cycles at 95 °C for 10 s, 57 °C for 30 s and 72 °C for 30 s, and one final extension at 95 °C for 10 s. The melting curves were determined using a denaturation gradient from 65 to 95 °C with an increase of 0.5 °C every 5 s.
The V410L mutation was identified according to Hadii et al. [32] with a final volume of 21 μL containing 10 μL of Perfecta Sybr® Green Supermix (Quanta 95054-500, Beverly, MA, USA), 0.1 μL of each of the primers L410(f) (5′-GCGGGCATCTTCTTGGGTTCGTTCTACCATT-3′) and V410(f) (5′-GCGGGCAGGGCGGCGGGGGCGGGGCCATCTTCTTGGGTTCGTTCTACCGTG-3′), 0.2 μL of primer L410(r) (5′-TTCTTCCTCGGCGGCCTCTT-3′), 9.6 μL of water, and 1.0 μL of the DNA template. Amplification reactions were performed with the following thermal profile: an initial denaturation cycle at 95 °C for 3 min, followed by 39 cycles at 95 °C for 10 s, 95 °C for 10 s, and 60 °C for 10 s. The melting curves were determined by a denaturation gradient from 65 to 95 °C with an increase of 0.2 °C every 10 s.
Homozygous-recessive (mutant) mosquitoes for the V1016I, F1534C, and V410L mutations were used as positive controls.
The melting curve results for the V1016I mutation were interpreted as follows. A single peak at 83 °C meant a dominant homozygous mosquito (V/V), that is, a susceptible genotype. A single peak at 76 °C corresponded to a homozygous recessive (I/I) mosquito, that is, a resistant genotype. Finally, the presence of the two peaks at 76 and 83 °C corresponded to a heterozygous mosquito (V/I). The melting curve results for the F1534C mutation were interpreted as follows. A single peak at 80 °C corresponded to a homozygous dominant mosquito (F/F), that is, a susceptible genotype. A single peak seen at 85 °C corresponded to a homozygous recessive (C/C) mosquito, that is, a resistant genotype. Finally, the presence of the two peaks at 80 and 85 °C meant a heterozygous mosquito (F/C). The melting curve results for the V410L mutation were interpreted as follows. A peak at 86.5 °C corresponded to a dominant homozygote (V/V); a peak at 83 °C corresponded to a recessive homozygote (L/L), susceptible and resistant genotypes, respectively. A peak at 83 °C and another at 86.5 °C corresponded to a heterozygote (V/L).

2.5. Spatial Distribution of Resistance to Pyrethroids

The insecticide resistance distribution was analyzed using the inverse distance weighting (IDW) interpolation procedure. The distance between the sampled sites and the mortality for each pyrethroid was used to estimate the mortality in the unsampled regions. IDW interpolation relies on the assumption that closer points have more similar values than locations separated by greater distances [33]. This association between distance and feature value is regulated by the power parameter in the IDW formula. To select the power value that gives the best estimates of predicted mortalities, the root squared mean error (RSME) was calculated by comparing true data values to interpolated values. Four different values (2, 2.5, 3, 3.5, and 4) were tested, and the power value with the lower RMSE was selected as it showed the lowest difference between the interpolated and true values. All analyses were conducted in R version 4.1.0 and QGIS version 3.24.
Allele and genotype frequencies were calculated; we verified that the populations were in Hardy–Weinberg equilibrium using a χ2 test. Wright’s FIS inbreeding coefficient was estimated, along with Wald’s correction [34,35]. In addition, the frequencies of the tri-locus haplotypes in the study populations were determined.
To assess whether the resistant haplotypes in Ae. aegypti were related to the frequency and intensity of resistance, a Spearman’s rank correlation analysis (α = 0.05) was performed between its frequency and the mortality rate caused by pyrethroid insecticides.

3. Results

3.1. Frequency and Intensity of Pyrethroid Resistance

The frequency of mortality obtained after exposure for 30 min to permethrin DD was less than 90% in 12 of the 15 populations analyzed, with mortality ranging from 30 to 76% (Figure 2). Only the populations from Manaure and Uribia showed values of 92 and 97% mortality, respectively, and the population of Barrancas exhibited 98% mortality showing susceptibility for the three pyrethroids.
All populations showed susceptibility to deltamethrin with 99–100% mortality (Figure 3).
In the case of lambda-cyhalothrin, the mortality frequency for all the populations was greater than the 90% threshold, suggesting the beginning of the emergence of resistance for the populations of Albania, Fonseca, Maicao, Riohacha, San Juan del Cesar, and Villanueva since they showed mortalities between 93 and 97%. The rest of the populations showed susceptibility to this insecticide (Figure 4).
When analyzing the intensity of resistance to permethrin, 100% mortality was evidenced at 2× DD in the populations from Manaure, Riohacha, and Uribia; only the Albania population showed 94% mortality, and the rest of the populations displayed mortality fluctuating between 68 and 90%. The populations from Distracción, Hatonuevo, La Jagua del Pilar, and Urumita showed a moderate intensity of resistance with mortalities of 92–96% when exposed to 5× DD. In contrast, for this same insecticide, 100% mortality was observed at the concentration of 10× in these same populations (Figure 2).
Regarding the intensity of resistance to lambda-cyhalothrin, mortalities between 99 and 100% were observed at the 2× dose in the populations from Albania, Fonseca, Maicao, Riohacha, San Juan del Cesar, and Villanueva (Figure 4).
The populations showed a higher frequency of resistance to permethrin than deltamethrin and lambda-cyhalothrin (Figure 5). The trends in mortality were significantly different between permethrin and deltamethrin and between permethrin and lambda-cyhalothrin (p < 0.01). Power values of 2, 3.5, and 4 resulted in lower RSME values for lambda-cyhalothrin, permethrin, and deltamethrin, respectively. IDW interpolation (Figure 5) showed that deltamethrin and lambda-cyhalothrin shared the same pattern where predicted mortalities were located between 90 and 100%. However, for permethrin, northern locations displayed higher predicted mortalities than southern populations, where the predicted values could be below 50%.

3.2. Allele and Genotype Frequencies of the kdr Mutations V410L, V1016I, and F1534C

Six hundred mosquitoes from all populations (40/population) were genotyped for V410L, V1016I, and F1534C mutations. For the V410L mutations, three genotypes, VV 410, VL 410, and LL 410, were detected in each field population, except for the populations from El Molino and Villanueva, where the LL 410 genotype was not found. The L 410 mutant allele was more frequent in the population from Uribia, with 0.61, and the least frequent in the population from El Molino 0.10. All populations were found to be in a Hardy–Weinberg equilibrium except for the Manaure population. An excess of heterozygotes was evidenced in the populations from Barrancas, Dibulla, El Molino, La Jagua del Pilar, Manaure, and San Juan del Cesar (Table 1).
For the V1016I mutation, the three genotypes VV 1016, VI 1016, and II 1016 were detected in each population, except for the populations from El Molino and Urumita, where genotype II 1016 was not observed. The I1016 mutant allele was the most frequent in the population from Uribia, with 0.61, and the least frequent in the population from El Molino 0.09. All populations were found to be in the Hardy–Weinberg equilibrium, except for the Hatonuevo and Urumita populations, and an excess of heterozygotes was evidenced in the populations from Barrancas, El Molino, La Jagua del Pilar, Maicao, Manaure, San Juan del Cesar, Urumita, and Villanueva (Table 2).
Regarding the F1534C mutation, the CC 1534 genotype was found in all the populations examined. The allelic frequency for the populations from Manaure, Distracción, Riohacha, and Albania ranged between 0.91 and 0.95, and for the population from Villanueva, an allele frequency of 0.69 was observed. The mutation was found to be fixed in the population from Uribia. Most of the populations were in the Hardy–Weinberg equilibrium, except those from Dibulla, Distraction, and Uribia. On the other hand, the populations from Albania, Dibulla, Fonseca, Hatonuevo, Maicao, Manaure, Urumita, and Villanueva showed an excess of heterozygotes (Table 3).
A total of 13 tri-locus combinations were detected among the populations (Figure 6). The homozygous triple-resistant haplotype (II1016/CC1534/LL410) occurred in 11 populations with very low frequencies. The total frequency was 6.6% of the total number of mosquitoes analyzed. Regarding the co-occurrence of kdr mutations, a total of triple-homozygous mutants was found in 6.6% of the 600 mosquitoes analyzed, where this haplotype was found in the populations from Albania (7.5%), Barrancas (5%), Distracción (5%), Fonseca (2.5%), Hatonuevo (5%), La Jagua del Pilar (2.5%), Maicao (5%), Manaure (10%), Riohacha (10%), San Juan del Cesar (7.5%), and Uribia (40%).
The most frequent haplotypes were VV410/VV1016/CC1534, VL410/VI1016/CC1534, and VV410/VV1016/FC1534, with 32, 31, and 16% of the total mosquitoes, respectively (Figure 6).
When analyzing the association of the mortality rate when exposing the mosquitoes to 1× DD with the frequencies of the tri-locus haplotypes for the three insecticides, a significant positive association was only found for permethrin with the VL410/VV1016/CC1534 haplotype (r = 0.217, p < 0.01). The same was seen in the case of deltamethrin, with the maximum association value for the VL410/VI1016/CC1534 haplotype (r = 0.25, p < 0.01), and in the case of lambda-cyhalothrin, with the maximum association value for the VL/VV/CC haplotype (r = 0.24, p < 0.01). Frequencies of the triple-resistant haplotype (LL410/II1016/CC1534) were significantly associated with mortality rates when the mosquitoes were exposed to 2X DD permethrin (r = 0.66, p < 0.05).

4. Discussion

Due to the selection pressure exerted by lambda-cyhalothrin and deltamethrin in Colombia, resistance to these and other pyrethroid-type molecules has been reported in different populations of Ae. aegypti [15]. Specifically, in the department of La Guajira, there has been a history of resistance to lambda-cyhalothrin, deltamethrin, and permethrin in the populations from San Juan del Cesar and for lambda-cyhalothrin in Riohacha [21,23]. The results found in this work are consistent with these previous studies, except for the susceptibility result found for deltamethrin in the population from San Juan del Cesar. This may be due to the fact that low resistance to deltamethrin had previously been reported in this population through the determination of the resistance factor by a lethal concentration 50 (LC50). By contrast, in this work, susceptibility was determined through DD [21].
On the other hand, the results of the present work constitute the first record that expands the knowledge of the susceptibility status of the thirteen remaining populations in the department of La Guajira, generally finding resistance to permethrin, the populations in danger of developing resistance to lambda-cyhalothrin, and susceptibility to deltamethrin, except in the population of Barrancas, which exhibited susceptibility to the three insecticides. Resistance to lambda-cyhalothrin has also been widely reported in Colombia in the departments of Chocó, Antioquia, Putumayo, Atlántico, Cesar, La Guajira, Sucre, Córdoba, Bolívar, Cundinamarca, Santander, Caquetá, Guaviare, Meta, and Casanare [21,23,36,37,38,39]. For deltamethrin in Colombia, most Ae. aegypti populations have been reported as susceptible; there is also a record of resistance in populations from the departments of Bolívar, Cesar, Córdoba, Atlántico, Cundinamarca, Caquetá, and Casanare [21,36,39]. This variability in the status of susceptibility to lambda-cyhalothrin and deltamethrin may be due to the difference in the selection pressure exerted with these insecticides, considering the differences in the incidence of dengue in Colombia, if deemed that there is a heterogeneous pattern in the transmission of this disease in the country when finding hyperendemic, endemic, and hypoendemic areas [40]. In addition, it is important to note that the use of insecticides in agriculture is low for the department of La Guajira, considering that most of the rural land use is for livestock production (82%) and only approximately 9% is for agricultural production [41]. However, from 1960 to the late 1980s, the Cesar River Valley, south of the department of La Guajira, was the main cotton district in Colombia; DDT, chlorpyrifos, and methomyl were the main insecticides used for pest control [42].
The resistance to permethrin found in the populations evaluated has also been recorded in other populations of the vector in Colombia [25,37,39]. It may be due to cross-resistance with other pyrethroid-type molecules such as lambda-cyhalothrin and/or organochlorine DDT, considering that in Colombia, this insecticide was suspended at the beginning of the 1990s by Resolution 10255 of the Health Ministry [21]. Cross-resistance has also been reported with other insecticides, such as alpha-cypermethrin, in populations resistant to pyrethroids and DDT in Colombia [24]. All this evidence demonstrates the importance of evaluating over time both the insecticides that have been applied to control this species, as well as other insecticides not used by government agencies, to contribute to the prevention of resistance and provide other alternatives for vector control.
Regarding the spatial distribution of insecticide resistance, we showed the resistance profile for three pyrethroids in unsampled regions within La Guajira using IDW interpolation. It is noteworthy that permethrin was the only insecticide with mortalities below 90% for southern areas. However, all results must be interpreted with caution, considering the limitations of the interpolation technique. Although highly implemented, the IDW method is a deterministic procedure that does not allow the calculation of variances of predicted, unmeasured sites. Additionally, this procedure relies on the distance between the sampled locations and the non-variable power value to predict the values of unmeasured sites [43]. It is also important to mention that other factors besides distance (e.g., insecticide application and mosquito population dynamics) influence the distribution of insecticide resistance.
Pyrethroid resistance is mainly associated with the presence of kdr mutations in the para gene of VGSC in insects such as Ae. aegypti. In this species, the kdr mutations G923V, L982W, I1011M, V1016G, I1011V, V1016I, F1534C, S989P, D1763Y, T1520I, and V410L have been identified worldwide with differences in allele and genotype frequencies between Ae. aegypti populations reported [30,31,44,45,46,47,48].
In the Americas, the V1016I, F1534C, and V410L mutations have been mainly associated with resistance to pyrethroids in populations in Mexico [30,49,50,51], Brazil [32,52,53], Colombia [21,23,24,25], Perú [54], Costa Rica [55], Puerto Rico [56], and Venezuela [57].
These mutations have also been reported and associated with resistance in populations of this vector in Colombia and in the departments of Antioquia, Valle del Cauca, Atlántico, Cesar, Bolívar, Magdalena, Sucre, Córdoba, Meta, Santander, and Quindío, with allele frequencies for V1016I between 0.02 and 0.72, for F1534C between 0.44 and 1.0, and for V410L between 0.05 and 0.72 [21,22,23,24,25,58,59]. Specifically, in the department of La Guajira, the kdr mutations identified were V1016I and F1534C in the San Juan del Cesar populations and V1016I, F1534C, and V410L in Riohacha [21,23,24]. In the case of the San Juan del Cesar population, allele frequencies of 0.28 and 0.77 have been reported for I1016 and C1534, respectively. By contrast, in the present study, these frequencies increased to 0.33 and 0.78, respectively. In the Riohacha population, allele frequencies were 0.25 for I016 and 0.71 for C1534, while in the present study, there was an increase between 0.31 and 0.91, respectively.
Regarding the allele frequency for L410 in a population of this same municipality, a frequency of 0.36 was determined, and we found in the present study a decrease to 0.31 [21,23,24]. The other populations had no history of identifying kdr mutations, so the results of this work constitute the first record for these populations, in which allele frequencies for L410 were between 0.10 and 0.61, for I1016 between 0.15 and 0.60, and for C1534 between 0.68 and 1.0. It is important to highlight that the highest frequency for the mutations was identified in the Uribia population, which is an isolated population from a geographical point of view in the upper Guajira that corresponds to a desert area and has been subjected in recent years to high selection pressure due to the increase in dengue incidence and the introduction of chikungunya and Zika viruses.
For the department of La Guajira, there are only two previous studies on the co-occurrence of kdr mutations in populations of Ae. aegypti; when comparing the results of these works with those obtained in the present study, an increase in the frequency of the double-homozygous mutant (II/CC) haplotype was observed in the municipality of San Juan del Cesar [24] and triple-homozygous mutant (II/CC/LL) in the Riohacha district [23]. Considering the results obtained, the evident resistance to permethrin, and the risk of developing resistance to lambda-cyhalothrin, the vector control program in La Guajira recommends malathion, pirimiphos-methyl, and deltamethrin for adult control [17].

5. Conclusions

The populations from La Guajira, Colombia, were resistant to permethrin with moderate intensity, associated with the frequency of the triple-mutated haplotype (LL410/II1016/CC1534), and were susceptible to deltamethrin and most to lambda-cyhalothrin, with six populations at risk of developing resistance.

Author Contributions

Conceptualization, R.M.-S., P.P.-L., G.P.-G. and A.E.F.; Methodology, R.M.-S., Z.F.-R., J.M.C.-C. and E.S.-A.; Software, A.E.F. and A.E.J.-V.; Validation, R.M.-S., P.P.-L., G.P.-G. and D.G.-C.; Formal analysis, A.E.F. and A.E.J.-V.; Investigation, R.M.-S., P.P.-L., G.P.-G., D.G.-C. and A.E.F.; Resources, R.M.-S.; Data curation, R.M.-S., Z.F.-R., J.M.C.-C., E.S.-A., D.G.-C., P.P.-L., G.P.-G., A.E.J.-V. and A.E.F.; Writing—original draft preparation, R.M.-S., P.P.-L. and A.E.F.; Writing—review, and editing, R.M.-S., Z.F.-R., J.M.C.-C., E.S.-A., D.G.-C., P.P.-L., A.E.J.-V. and A.E.F.; Visualization, A.E.F. and R.M.-S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by Universidad Libre seccional Barranquilla, Colombia, through internal funding.

Data Availability Statement

All the required data relevant to the presented study are included in the manuscript.

Acknowledgments

We express our gratitude to the technicians of the Medical Entomology Laboratory of the Department of Health of La Guajira and to the students of the Bacteriology program from the Universidad Libre seccional Barranquilla, for their support in field and laboratory activities. The icon in the graphical abstract was made by Flowicon and Freepik from www.flaticon.com (accessed on 13 December 2022).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Arredondo-García, J.L.; Méndez-Herrera, A.; Medina-Cortina, H. Arbovirus en Latinoamérica. Acta Pediatr. Mex. 2016, 37, 111–131. Available online: https://www.medigraphic.com/pdfs/actpedmex/apm-2016/apm162h.pdf (accessed on 25 July 2022). [CrossRef] [Green Version]
  2. Instituto Nacional de Salud (INS). Gestión para la Vigilancia Entomológica y Control de la Transmisión de Dengue. 2010. Available online: https://www.paho.org/col/dmdocuments/Entomologia_DENGUE.pdf (accessed on 25 July 2022).
  3. Organización Panamericana de la Salud (OPS). Actualización Epidemiológica Dengue y otras Arbovirosis. 2020. Available online: https://www.paho.org/es/documentos/actualizacion-epidemiologica-dengue-otras-arbovirosis-10-junio-2020 (accessed on 25 July 2022).
  4. Organización Panamericana de la Salud (OPS). Actualización Epidemiológica. 2021. Available online: https://www.paho.org/es/documentos/actualizacion-epidemiologica-arbovirosis-contexto-covid-19-2-julio-2021 (accessed on 25 July 2022).
  5. Organización Panamericana de la Salud (OPS). Actualización Epidemiológica Anual para Dengue, Chikungunya y Zika. 2022. Available online: https://ais.paho.org/ha_viz/arbo/pdf/OPS%20Arbo%20Boletin%202022.pdf (accessed on 18 October 2022).
  6. Instituto Nacional de Salud (INS). Informe del Evento Dengue, Colombia, Hasta el Décimo Tercer Periodo Epidemiológico. 2015. Available online: https://www.ins.gov.co/buscador-eventos/Informesdeevento/Dengue%202015.pdf (accessed on 23 July 2022).
  7. Instituto Nacional de Salud (INS). Informe del Evento Dengue, Colombia, Hasta el Décimo Tercer Periodo Epidemiológico. 2016. Available online: https://www.ins.gov.co/buscador-eventos/Informesdeevento/Dengue%202016.pdf (accessed on 23 July 2022).
  8. Instituto Nacional de Salud (INS). Informe del Evento Dengue, Colombia. 2017. Available online: https://www.ins.gov.co/buscador-eventos/Informesdeevento/DENGUE%202017.pdf (accessed on 23 July 2022).
  9. Instituto Nacional de Salud (INS). Informe del Evento Dengue, Colombia. 2018. Available online: https://www.ins.gov.co/buscador-eventos/Informesdeevento/DENGUE_2018.pdf (accessed on 23 July 2022).
  10. Instituto Nacional de Salud (INS). Informe del Evento Dengue, Colombia. 2019. Available online: https://www.ins.gov.co/buscador-eventos/Informesdeevento/DENGUE_2019.pdf (accessed on 23 July 2022).
  11. Secretaria de Salud de la Guajira, A. (Secretaria de Salud de la Guajira, Riohacha, La Guajira, Colombia); Incidencia Dengue para el Departamento de La Guajira (Colombia). 2020; (Unpublished information). [Google Scholar]
  12. Informe Quincenal Epidemiológico Nacional (IQEN). Chikungunya en Colombia, el Inicio de la Transmisión Autóctona. 2014; Volume 19, pp. 261–279. Available online: https://www.ins.gov.co/buscador-eventos/IQEN/IQEN%20vol%2019%202014%20num%2018.pdf (accessed on 20 July 2022).
  13. Instituto Nacional de Salud. Informe de Evento, Enfermedad por virus Zika y Chikungunya. 2019. Available online: https://www.ins.gov.co/buscador-eventos/Informesdeevento/chikungunya%20y%20enfermedad%20por%20virus%20Zika_2019.pdf (accessed on 23 July 2022).
  14. World Health Organization (WHO). Pesticides and Their Application for the Control of Vectors and Pests of Public Health Importance. 2006. Available online: https://apps.who.int/iris/bitstream/handle/10665/69223/WHO_CDS_NTD_WHOPES_GCDPP_2006.1_eng.pdf?sequence=1&isAllowed=y (accessed on 22 July 2022).
  15. Maestre-Serrano, R. Susceptibility status of Aedes aegypti to insecticides in Colombia. In Insecticides-Pest Engineering; Perveeb, F., Ed.; InTech Publishers: Rijeka, Croatia, 2012; pp. 163–200. [Google Scholar] [CrossRef]
  16. Maestre, R.; Gómez, D. Dengue: Epidemiología, políticas públicas y resistencia de vectores a insecticidas. Rev. Cienc. Biomed. 2013, 4, 302–317. [Google Scholar]
  17. Secretaria de Salud de la Guajira, B. (Secretaria de Salud de la Guajira, Riohacha, La Guajira, Colombia); Insecticidas usados en el departamento de La Guajira para el control del dengue. 2022; (Unpublished information). [Google Scholar]
  18. Du, Y.; Nomura, Y.; Satar, G.; Hu, Z.; Nauen, R.; He, S.Y.; Zhorov, B.S.; Dong, K. Molecular evidence for dual pyrethroid-receptor sites on a mosquito sodium channel. Proc. Natl. Acad. Sci. USA 2013, 110, 11785–11790. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Catterall, W.A. Voltage-gated sodium channels at 60: Structure, function and pathophysiology. J. Physiol. 2012, 590, 2577–2589. [Google Scholar] [CrossRef] [PubMed]
  20. Maestre, R.; Rey, G.; De Las Salas, J.; Vergara, C.; Santacoloma, L.; Goenaga, S.; Carrasquilla, M.C. Susceptibility status of Aedes aegypti to insecticides in Atlántico (Colombia). Rev. Colomb. Entomol. 2010, 36, 242–248. [Google Scholar] [CrossRef]
  21. Maestre-Serrano, R.; Gomez-Camargo, D.; Ponce-Garcia, G.; Flores, A.E. Susceptibility to insecticides and resistance mechanisms in Aedes aegypti from the Colombian Caribbean Region. Pestic. Biochem. Physiol. 2014, 116, 63–73. [Google Scholar] [CrossRef]
  22. Atencia, M.C.; Pérez, M.D.J.; Jaramillo, M.C.; Caldera, S.M.; Cochero, S.; Bejarano, E.E. Primer reporte de la mutación F1534C asociada con resistencia cruzada a DDT y piretroides en Aedes aegypti en Colombia. Biomedica 2016, 36, 432–437. [Google Scholar] [CrossRef] [Green Version]
  23. Granada, Y.; Mejía-Jaramillo, A.M.; 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, 23. [Google Scholar] [CrossRef] [Green Version]
  24. Maestre-Serrano, R.; Pareja-Loaiza, P.; Gomez Camargo, D.; Ponce-García, G.; Flores, A.E. 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–1688. [Google Scholar] [CrossRef] [Green Version]
  25. Pareja-Loaiza, P.X.; Santacoloma Varon, L.; Rey Vega, G.; Gómez-Camargo, D.; Maestre-Serrano, R.; Lenhart, A. Mechanisms associated with pyrethroid resistance in populations of Aedes aegypti (Diptera: Culicidae) from the Caribbean coast of Colombia. PLoS ONE 2020, 15, e0228695. [Google Scholar] [CrossRef] [PubMed]
  26. Brogdon, W.; Chan, A. Guideline for Evaluating Insecticide Resistance in Vectors Using the CDC Bottle Bioassay (with Inserts 1 in 2012 and 2 in 2014). 2010. Available online: http://www.cdc.gov/malaria/features/bioassay.html (accessed on 15 January 2022).
  27. Brogdon, W.G.; McAllister, J.C. Simplification of adult mosquito bioassays through use of time-mortality determinations in glass bottles. J. Am. Mosq. Control Assoc. 1998, 14, 159–164. [Google Scholar] [PubMed]
  28. Abbott, W.S. A method of computing the effectiveness of an insecticide. J. Econ. Entomol. 1925, 18, 265–267. [Google Scholar] [CrossRef]
  29. World Health Organization (WHO). Test Procedures for Insecticide Resistance Monitoring in Malaria Vector Mosquitoes, 2nd ed.; World Health Organization, 2016; Available online: https://apps.who.int/iris/handle/10665/250677 (accessed on 12 January 2022).
  30. Saavedra-Rodriguez, K.; Urdaneta-Marquez, L.; Rajatileka, S.; Moulton, M.; Flores, A.E.; Fernandez-Salas, I.; Bisset, J.; Rodriguez, M.; McCall, P.J.; Donnelly, M.J.; 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–798. [Google Scholar] [CrossRef] [Green Version]
  31. Yanola, J.; Somboon, P.; Walton, C.; Nachaiwieng, W.; Somwang, P.; Prapanthadara, L.A. 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–509. [Google Scholar] [CrossRef]
  32. Haddi, K.; Tomé, H.V.V.; Du, Y.; Valbon, W.R.; Nomura, Y.; Martins, G.F.; Dong, K.; Oliveira, E.E. 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. [Google Scholar] [CrossRef] [Green Version]
  33. Wu, Y.H.; Hung, M.C. Comparison of Spatial Interpolation Techniques Using Visualization and Quantitative Assessment. In Applications of Spatial Statistics; IntechOpen: London, UK, 2016. [Google Scholar] [CrossRef] [Green Version]
  34. Wright, S. Systems of mating. II. The effects of inbreeding on the genetic composition of a population. Genetics 1921, 6, 124–143. [Google Scholar] [CrossRef]
  35. Agresti, A.; Coull, B.A. Approximate is better than “exact” for interval estimation of binominal proportions. Am. Stat. 1998, 52, 119–126. [Google Scholar]
  36. Santacoloma Varón, L.; Chaves Córdoba, B.; Brochero, H.L. Susceptibilidad de Aedes aegypti a DDT, deltametrina y lambdacialotrina en Colombia. Rev. Panam. Salud Pública 2010, 27, 66–73. [Google Scholar] [CrossRef] [PubMed]
  37. Fonseca-González, I.; Quiñones, M.L.; Lenhart, A.; Brogdon, W.G. Insecticide resistance status of Aedes aegypti (L.) from Colombia. Pest Manag. Sci. 2011, 67, 430–437. [Google Scholar] [CrossRef]
  38. Ocampo, C.B.; Salazar-Terreros, M.J.; Mina, N.J.; McAllister, J.; Brogdon, W. Insecticide resistance status of Aedes aegypti in 10 localities in Colombia. Acta Trop. 2011, 118, 37–44. [Google Scholar] [CrossRef] [PubMed]
  39. Ardila-Roldán, S.; Santacoloma, L.; Brochero, H. Estado de la sensibilidad a los insecticidas de uso en salud pública en poblaciones naturales de Aedes aegypti (Diptera: Culicidae) del departamento de Casanare, Colombia. Biomedica 2013, 33, 446–458. [Google Scholar] [CrossRef] [PubMed]
  40. Padilla, J.C.; Lizarazo, F.E.; Murillo, O.L.; Mendigaña, F.A.; Pachón, E.; Vera, M.J. Epidemiología de las principales enfermedades transmitidas por vectores en Colombia, 1990–2016. Biomedica 2017, 37, 27–40. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Guajira 360°. Available online: https://guajira360.org/informe-economico-reactivacion-del-sector-agropecuario-de-la-guajira/ (accessed on 22 December 2022).
  42. Wagner Medina, M. Las Huellas Ambientales del oro Blanco: La Expansión Algodonera en el Valle del río Cesar. (1950–1980). Master’s Dissertation/Tesis, Universidad de los Andes, Colombia, July 2011. Available online: https://repositorio.uniandes.edu.co/bitstream/handle/1992/11452/u470720.pdf?sequence=1&isAllowed=y (accessed on 21 December 2022).
  43. Lu, G.Y.; Wong, D.W. An adaptive inverse-distance weighting spatial interpolation technique. Comput. Geosci. 2008, 34, 1044–1055. [Google Scholar] [CrossRef]
  44. Chang, C.; Shen, W.K.; Wang, T.T.; Lin, Y.H.; Hsu, E.L.; Dai, S.M. A novel amino acid substitution in a voltage-gated sodium channel is associated with knockdown resistance to permethrin in Aedes aegypti. Insect Biochem. Mol. Biol. 2009, 39, 272–278. [Google Scholar] [CrossRef]
  45. Srisawat, R.; Komalamisra, N.; Apiwathnasorn, C.; Paeporn, P.; Roytrakul, S.; Rongsriyam, Y.; Eshita, Y. Field-collected permethrin-resistant Aedes aegypti from central Thailand contain point mutations in the domain IIS6 of the sodium channel gene (KDR). Southeast Asian J. Trop. Med. Public Health 2012, 43, 1380–1386. [Google Scholar]
  46. Kushwah, R.B.; Dykes, C.L.; Kapoor, N.; Adak, T.; Singh, O.P. Pyrethroid-resistance and presence of two knockdown resistance (kdr) mutations, F1534C and a novel mutation T1520I, in Indian Aedes aegypti. PLoS Negl. Trop. Dis. 2015, 9, e3332. [Google Scholar] [CrossRef] [Green Version]
  47. Brengues, C.; Hawkes, N.J.; Chandre, F.; McCarroll, L.; Duchon, S.; Guillet, P.; Manguin, S.; Morgan, J.C.; Hemingway, J. 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. [Google Scholar] [CrossRef] [Green Version]
  48. Collet, M.L.; Frizzo, C.; Orlandin, E.; Rona, L.D.P.; Nascimento, J.C.; Montano, M.A.E.; Müller, G.A.; Wagner, G. Frequency of the Val1016Ile mutation on the kdr gene in Aedes aegypti (Diptera: Culicidae) in south Brazil. GMR 2016, 15, gmr15048940. [Google Scholar] [CrossRef]
  49. Siller, Q.; Ponce, G.; Lozano, S.; Flores, A.E. 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–362. [Google Scholar] [CrossRef]
  50. Vera-Maloof, F.Z.; Saavedra-Rodriguez, K.; Elizondo-Quiroga, A.E.; Lozano-Fuentes, S.; Black, W.C., IV. Coevolution of the Ile1,016 and Cys1,534 Mutations in the Voltage Gated Sodium Channel Gene of Aedes aegypti in Mexico. PLoS Negl. Trop. Dis. 2015, 9, e0004263. [Google Scholar] [CrossRef] [PubMed]
  51. Villanueva-Segura, O.K.; Ontiveros-Zapata, K.A.; Lopez-Monroy, B.; Ponce-Garcia, G.; Gutierrez-Rodriguez, S.M.; Davila-Barboza, J.A.; Mora-Jasso, E.J.; Flores, A.E. 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–223. [Google Scholar] [CrossRef] [PubMed]
  52. Aguirre-Obando, O.A.; Pietrobon, A.J.; Dalla Bona, A.C.; Navarro-Silva, M.A. Contrasting patterns of insecticide resistance and knockdown resistance (kdr) in Aedes aegypti populations from Jacarezinho (Brazil) after a dengue outbreak. Rev. Bras. Entomol. 2016, 60, 94–100. [Google Scholar] [CrossRef] [Green Version]
  53. Melo Costa, M.; Campos, K.B.; Brito, L.P.; Roux, E.; Melo Rodovalho, C.; Bellinato, D.F.; Lima, J.B.P.; Martins, A.J. Kdr genotyping in Aedes aegypti from Brazil on a nation-wide scale from 2017 to 2018. Sci. Rep. 2020, 10, 13267. [Google Scholar] [CrossRef] [PubMed]
  54. Pinto, J.; Palomino, M.; Mendoza-Uribe, L.; Sinti, C.; Liebman, K.A.; Lenhart, A. Susceptibility to insecticides and resistance mechanisms in three populations of Aedes aegypti from Peru. Parasit. Vectors 2019, 12, 494. [Google Scholar] [CrossRef] [Green Version]
  55. Zardkoohi, A.; Castañeda, D.; Lol, J.C.; Castillo, C.; Lopez, F.; Marín Rodriguez, R.; Padilla, N. Co-occurrence of kdr mutations V1016I and F1534C and its association with phenotypic resistance to pyrethroids in Aedes aegypti (Diptera: Culicidae) populations from Costa Rica. J. Med. Entomol. 2020, 57, 830–836. [Google Scholar] [CrossRef]
  56. Ponce-García, G.; Del Río-Galvan, S.; Barrera, R.; Saavedra-Rodriguez, K.; Villanueva-Segura, K.; Felix, G.; Amador, M.; Flores, A.E. Knockdown Resistance Mutations in Aedes aegypti (Diptera: Culicidae) From Puerto Rico. J. Med. Entomol. 2016, 53, 1410–1414. [Google Scholar] [CrossRef]
  57. Alvarez, L.C.; Ponce-García, G.; Saavedra-Rodriguez, K.; Lopez, B.; Flores, A.E. Frequency of V1016I and F1534C mutations in the voltage-gated sodium channel gene in Aedes aegypti in Venezuela. Pest Manag. Sci. 2015, 71, 863–869. [Google Scholar] [CrossRef]
  58. Aguirre-Obando, O.A.; Dalla Bona, A.C.; Duque, L.J.E.; Navarro-Silva, M.A. Insecticide resistance and genetic variability in natural populations of Aedes (Stegomyia) aegypti (Diptera: Culicidae) from Colombia. Zoologia 2015, 32, 14–22. [Google Scholar] [CrossRef] [Green Version]
  59. Aponte, A.; Penilla, R.P.; Rodríguez, A.D.; Ocampo, C.B. Mechanisms of pyrethroid resistance in Aedes (Stegomyia) aegypti from Colombia. Acta Trop. 2019, 191, 146–154. [Google Scholar] [CrossRef]
Figure 1. Aedes aegypti collection sites in the department of La Guajira, Colombia.
Figure 1. Aedes aegypti collection sites in the department of La Guajira, Colombia.
Insects 14 00031 g001
Figure 2. Mortality (with standard deviations) of Ae. aegypti following exposure to 1× DD, 2× DD, 5× DD, and 10× DD of permethrin in bottle bioassays. The red lines indicate the resistance threshold (90% mortality).
Figure 2. Mortality (with standard deviations) of Ae. aegypti following exposure to 1× DD, 2× DD, 5× DD, and 10× DD of permethrin in bottle bioassays. The red lines indicate the resistance threshold (90% mortality).
Insects 14 00031 g002
Figure 3. Mortality (with standard deviation) of Ae. aegypti following exposure to 1× DD of deltamethrin in bottle bioassays. The red lines indicate the resistance threshold (90% mortality).
Figure 3. Mortality (with standard deviation) of Ae. aegypti following exposure to 1× DD of deltamethrin in bottle bioassays. The red lines indicate the resistance threshold (90% mortality).
Insects 14 00031 g003
Figure 4. Mortality (with standard deviation) of Ae. aegypti following exposure to 1× DD and 2× DD of lambda-cyhalothrin in bottle bioassays. The red lines indicate the resistance threshold (90% mortality).
Figure 4. Mortality (with standard deviation) of Ae. aegypti following exposure to 1× DD and 2× DD of lambda-cyhalothrin in bottle bioassays. The red lines indicate the resistance threshold (90% mortality).
Insects 14 00031 g004
Figure 5. Distribution of pyrethroid resistance (% of mortality to the DD) to (a) permethrin, (b) deltamethrin, and (c) lambda–cyhalothrin in Ae. aegypti populations from La Guajira, Colombia using IDW interpolation.
Figure 5. Distribution of pyrethroid resistance (% of mortality to the DD) to (a) permethrin, (b) deltamethrin, and (c) lambda–cyhalothrin in Ae. aegypti populations from La Guajira, Colombia using IDW interpolation.
Insects 14 00031 g005
Figure 6. Frequencies of tri-locus haplotypes in each population of Ae. aegypti from La Guajira, Colombia. Haplotype order: 410/1016/1534.
Figure 6. Frequencies of tri-locus haplotypes in each population of Ae. aegypti from La Guajira, Colombia. Haplotype order: 410/1016/1534.
Insects 14 00031 g006
Table 1. Allele and genotype frequencies of V410L in populations of Ae. aegypti from the department of La Guajira, Colombia.
Table 1. Allele and genotype frequencies of V410L in populations of Ae. aegypti from the department of La Guajira, Colombia.
PopulationnVVVLLLFrequency (IC 95%)FISx2 Hardy-Weinbergp Value
Manaure4092740.44 (0.30–0.59)−0.375.510.02
Uribia40717160.61 (0.46–0.75)0.100.430.50
Maicao40161860.38 (0.24–0.53)0.040.060.80
Riohacha40191740.31 (0.20–0.47)0.0100.94
Albania40241330.24 (0.13–0.39)0.100.420.51
Dibulla40201730.29 (0.17–0.44)−0.030.050.81
Hatonuevo4031720.14 (0.06–0.28)0.262.740.10
Barrancas40172120.31 (0.20–0.47)−0.221.960.16
Distracción40241330.24 (0.13–0.39)0.100.420.51
Fonseca40281020.18 (0.09–0.32)0.130.720.40
San Juan del Cesar40142150.39 (0.25–0.54)−0.100.440.50
El Molino4032800.10 (0.04–0.24)−0.110.490.50
Villanueva40211900.24 (0.13–0.39)−0.313.880.05
Urumita40162040.35 (0.22–0.50)−0.090.390.53
La Jagua del Pilar40241510.21 (0.11–0.36)−0.120.580.44
Abbreviations: n, sample size; VV, wild type; VL, heterozygotes; LL, homozygotes resistant.
Table 2. Allele and genotype frequencies of V1016I in populations of Ae. aegypti from La Guajira, Colombia.
Table 2. Allele and genotype frequencies of V1016I in populations of Ae. aegypti from La Guajira, Colombia.
PopulationnVVVIIIFrequency (IC 95%)FISx2 Hardy-Weinbergp Value
Manaure40162040.35 (0.22–0.50)−0.10.390.53
Uribia40717160.61 (0.46–0.75)0.10.440.51
Maicao40191920.29 (0.17–0.44)−0.161.020.31
Riohacha40191740.31 (0.19–0.47)0.0100.94
Albania40241330.24 (0.13–0.39)0.10.420.52
Dibulla40281110.16 (0.08–0.31)001
Hatonuevo4033520.11 (0.05–0.25)0.375.590.01
Barrancas40182020.30 (0.18–0.45)−0.191.450.23
Distracción40251320.21 (0.11–0.36)0.020.030.86
Fonseca40291010.15 (0.07–0.29)0.020.020.9
San Juan del Cesar40172030.33 (0.20–0.48)−0.140.780.38
El Molino4033700.09 (0.03–0.21)−0.0960.360.55
Villanueva40221710.24 (0.13–0.39)−0.171.20.27
Urumita40192100.26 (0.15–0.42)−0.355.060.02
La Jagua del Pilar40251410.20 (0.10–0.35)−0.090.350.55
Abbreviations: n, sample size; VV, wild type; VI, heterozygotes; II, homozygotes resistant.
Table 3. Allele and genotype frequencies of F1534C in populations of Ae. aegypti from La Guajira, Colombia.
Table 3. Allele and genotype frequencies of F1534C in populations of Ae. aegypti from La Guajira, Colombia.
PopulationnFFFCCCFrequency (IC95%)FISx2 Hardy-Weinbergp Value
Manaure4007330.91 (0.79–0.97)−0.10.370.54
Uribia4000401.00 (0.91–1.00)---
Maicao40011290.86 (0.72–0.94)−0.161.020.31
Riohacha4015340.91 (0.79–0.97)0.222.050.15
Albania4004360.95 (0.84–0.99)−0.050.110.74
Dibulla40122170.70 (0.55–0.82)−0.375.40.02
Hatonuevo40115240.79 (0.64–0.89)−0.120.580.44
Barrancas40110290.85 (0.71–0.93)0.020.020.9
Distracción4013360.94 (0.82–0.98)0.365.180.02
Fonseca40113260.81 (0.67–0.90)−0.070.180.67
San Juan del Cesar40312250.78 (0.63–0.88)0.140.780.38
El Molino40210280.83 (0.69–0.91)0.130.720.40
Villanueva40319180.69 (0.53–0.81)−0.110.440.5
Urumita40215230.76 (0.61–0.87)−0.030.050.82
La Jagua del Pilar40213250.79 (0.64–0.89)0.020.030.85
Abbreviations: n, sample size; FF, wild type; FC, heterozygotes; CC, homozygotes resistant.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Maestre-Serrano, R.; Flórez-Rivadeneira, Z.; Castro-Camacho, J.M.; Soto-Arenilla, E.; Gómez-Camargo, D.; Pareja-Loaiza, P.; Ponce-Garcia, G.; Juache-Villagrana, A.E.; Flores, A.E. Spatial Distribution of Pyrethroid Resistance and kdr Mutations in Aedes aegypti from La Guajira, Colombia. Insects 2023, 14, 31. https://doi.org/10.3390/insects14010031

AMA Style

Maestre-Serrano R, Flórez-Rivadeneira Z, Castro-Camacho JM, Soto-Arenilla E, Gómez-Camargo D, Pareja-Loaiza P, Ponce-Garcia G, Juache-Villagrana AE, Flores AE. Spatial Distribution of Pyrethroid Resistance and kdr Mutations in Aedes aegypti from La Guajira, Colombia. Insects. 2023; 14(1):31. https://doi.org/10.3390/insects14010031

Chicago/Turabian Style

Maestre-Serrano, Ronald, Zulibeth Flórez-Rivadeneira, Juan M. Castro-Camacho, Eva Soto-Arenilla, Doris Gómez-Camargo, Paula Pareja-Loaiza, Gustavo Ponce-Garcia, Alan E. Juache-Villagrana, and Adriana E. Flores. 2023. "Spatial Distribution of Pyrethroid Resistance and kdr Mutations in Aedes aegypti from La Guajira, Colombia" Insects 14, no. 1: 31. https://doi.org/10.3390/insects14010031

APA Style

Maestre-Serrano, R., Flórez-Rivadeneira, Z., Castro-Camacho, J. M., Soto-Arenilla, E., Gómez-Camargo, D., Pareja-Loaiza, P., Ponce-Garcia, G., Juache-Villagrana, A. E., & Flores, A. E. (2023). Spatial Distribution of Pyrethroid Resistance and kdr Mutations in Aedes aegypti from La Guajira, Colombia. Insects, 14(1), 31. https://doi.org/10.3390/insects14010031

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop