Levels of insecticide resistance to deltamethrin, malathion, and temephos, and associated mechanisms in Aedes aegypti mosquitoes from the Guadeloupe and Saint Martin islands (French West Indies)
© The Author(s). 2017
Received: 14 July 2016
Accepted: 6 February 2017
Published: 10 February 2017
In the Guadeloupe and Saint Martin islands, Aedes aegypti mosquitoes are the only recognized vectors of dengue, chikungunya, and Zika viruses. For around 40 years, malathion was used as a mosquito adulticide and temephos as a larvicide. Since the European Union banned the use of these two insecticide molecules in the first decade of the 21st century, deltamethrin and Bacillus thuringiensis var. israelensis are the remaining adulticide and larvicide, respectively, used in Guadeloupe. In order to improve the management of vector control activities in Guadeloupe and Saint Martin, we investigated Ae. aegypti resistance to and mechanisms associated with deltamethrin, malathion, and temephos.
Ae. aegypti mosquitoes were collected from six different localities of Guadeloupe and Saint Martin. Larvae were used for malathion and temephos bioassays, and adult mosquitoes for deltamethrin bioassays, following World Health Organization recommendations. Knockdown resistance (Kdr) genotyping for V1016I and F1534C mutations, and expression levels of eight enzymes involved in detoxification mechanisms were examined in comparison with the susceptible reference Bora Bora strain.
Resistance ratios (RR50) calculated for Ae. aegypti larvae showed high resistance levels to temephos (from 8.9 to 33.1-fold) and low resistance levels to malathion (from 1.7 to 4.4-fold). Adult females displayed moderate resistance levels to deltamethrin regarding the time necessary to affect 50% of individuals, varying from 8.0 to 28.1-fold. Molecular investigations on adult mosquitoes showed high resistant allele frequencies for V1016I and F1534C (from 85 to 96% and from 90 to 98%, respectively), as well as an overexpression of the glutathione S-transferase gene, GSTe2, the carboxylesterase CCEae3a, and the cytochrome genes 014614, CYP6BB2, CYP6M11, and CYP9J23.
Ae. aegypti populations from Guadeloupe and Saint Martin exhibit multiple resistance to organophosphates (temephos and malathion), and pyrethroids (deltamethrin). The mechanisms associated with these resistance patterns show strong frequencies of F1534C and V1016I Kdr mutations, and an over-expression of CCEae3a, GSTe2, and four cytochrome P450 genes (014614, CYP9J23, CYP6M11, CYP6BB2). These results will form the baseline for a deeper understanding of the insecticide resistance levels and associated mechanisms of Ae. aegypti populations and will be used to improve vector control strategies in Guadeloupe and Saint Martin.
KeywordsAedes aegypti Mosquitoes Insecticide resistance Deltamethrin Malathion Temephos Guadeloupe Saint Martin
Please see Additional file 1 for translations of the abstract into the six official working languages of the United Nations.
Recently, from December 2013 to January 2015, the chikungunya virus (CHIKV) intensely hit the Guadeloupe islands of the French West Indies, which resulted in a severe outbreak in which an estimated 40% (160,000 people) of the population became infected [1, 2]. Most of the Latin American countries and Caribbean islands were still reporting CHIKV outbreaks when a new threat was reported from Brazil—the arrival of the Zika virus (ZIKV) . This arbovirus, new for the American region, was first considered as a mild disease until the discovery of Guillain-Barré and microcephaly syndromes associated with ZIKV infections . Since then, the number of severe cases has drastically increased in Brazil and in most of the Latin American countries and Caribbean islands where ZIKV has been reported, and the first severe neurological cases associated with the ZIKV in the French West Indies were reported from Martinique in early 2016 . In the French West Indies, both CHIKV and ZIKV are transmitted by the mosquito species Aedes aegypti, the vector of the dengue viruses that cause epidemics in the region every 2 to 3 years .
Since no specific treatment and/or vaccines are commonly available against CHIKV and ZIKV, and since the dengue vaccine has been deployed only very recently in a small number of experimental areas , the control of vectors and personal protection against mosquito bites remain the only available tools for preventing and controlling these emerging arboviral diseases [8, 9].
Historically, vector control activities in Guadeloupe have been carried out by a Vector Control Agency from the French Ministry of Health local delegation (Agence Régionale de Santé, ARS) and by ensuring larval reduction through the elimination of breeding sites and larvicide treatments, as well as the spraying of adulticide outdoors and indoors . Biological larval control has also been tested using larvivorous fish in bigger tanks (J. Gustave, personal communication). The insecticides commonly used in the past have been temephos as a larvicide, and malathion and deltamethrin as adulticides . However, since 2009, temephos has been withdrawn from the list of vector control insecticides, following European Union recommendations, and has been replaced by biological products using Bacillus thuringiensis var. israelensis (Bti) . In 2010, the adulticide malathion was also withdrawn, consequently, the vector control activities in Guadeloupe were reorganized and now comprise routine elimination of breeding sites, the use of Bti as a larvicide, and the use of deltamethrin as an indoor adulticide.
In the French overseas departments of Guadeloupe, French Guiana, and Martinique, the surveillance of arbovirus epidemics is organized around a specific system called the dengue epidemics surveillance, alert and management program (PSAGE). Surveillance and control measures increase if necessary, according to the epidemiological situation, which is monitored all year round through a network of physicians in each territory who report weekly numbers of clinical cases for dengue, CHIKV, and now ZIKV. The PSAGE system proposes five operational situations [http://opac.invs.sante.fr/doc_num.php?explnum_id=3517]: phase 1) inter-epidemic situation with sporadic transmission; phase 2) some clusters of cases, seasonal increase, outdoor spraying of deltamethrin in and around clusters; phase 3) pre-epidemic alert with a significant increase in the number of clinical cases, outdoor spraying of deltamethrin in selected areas, reinforcement of entomological surveillance; phase 4) outbreak alert, intensification of outdoor spraying of deltamethrin in all places, communicating messages adapted to the epidemic context (including personal protection and reinforcement of personal elimination of breeding sites) to populations; and the final phase 5) end of the epidemic and evaluation of the lessons learned.
The large deployment of insecticides to control mosquito vectors has allowed the development of resistance worldwide . For the Guadeloupe islands, very few studies have been done on the resistance levels of the mosquito Ae. aegypti , and resistance mechanisms have not been investigated. The resistance of mosquitoes to chemical insecticides can be attributed to two main mechanisms . The first mechanism is the modification of the molecular site targeted by the chemical, and the second is an increased metabolism of the chemical through mutations and/or over-expression of detoxifying enzymes.
First report of Ae. aegypti insecticide resistance in Guadeloupe and Saint Martin are almost 20 years old, and resistance levels have probably significantly evolved with the increasing insecticide dosages. Thus, an update of Ae. aegypti resistance levels is urgently needed to improve the management of vector control activities in the Guadeloupe and Saint-Martin islands. Therefore, in this study, we investigated insecticide resistance levels and their associated molecular mechanisms for six Ae. aegypti populations collected in Guadeloupe and Saint Martin islands between January 2014 and October 2015.
Mosquito collection and rearing
Ae. aegypti mosquito populations used in the study
Breeding sites: type (number)
Date of collection
14 Jan 2014
tires (multiple); casks (4); buckets (2)
23 Jan 2015
30 Jan 2015
casks (2); bucket (1)
30 Jan 2015
aquatic plants (2); plant cutting (1)
16 Jan 2014
plant cutting (2)
16 Jan 2014
buckets (9); casks (3); small wastes (4); pot dish (1)
24 Sep 2015
plant cutting (1); bucket (1)
22 Jan 2014
16°15'30.9"N 61° 35'14.6"W
watering can (1)
22 Jan 2014
plant cutting (1); bucket (1)
13 Feb 2015
13 Feb 2015
Saint Martin East French part
abandoned boat (1); abandoned jacuzzi (1); tires (multiple); casks (2); discarded small containers (1)
04–05 Feb 2014
Saint-Martin West French part
tires (2); casks (3); small waste (1); plant cutting (1); pot dish (1); buckets (3)
04–05 Feb 2014
abandoned boat (1); tires (multiple); abandoned jet ski (1); discarded small containers (1)
25–26 Nov 2014
pot dish (2); flower vase (1); casks (4); watering can (1); bucket (2); tires (5)
07 Oct 2015
Larval and Adult test
casks (2); bucket (1); flower pot (1); watering can (1); tires (multiple)
06 Oct 2015
Larval and Adult test
Larval bioassays were performed following WHO recommendations . The late third- and early fourth-instar larvae were used for each mosquito population. Four replicates per concentration and five concentrations were tested with 25 larvae per replicate and per concentration. The insecticide alcohol dilutions were provided by the WHO Malaysia manufacturer (Vector Control Research Unit, School of Biological Sciences, Universiti Sains Malaysia). The insecticide solutions received were temephos (156.25 μg/ml and 31.25 μg/ml) and malathion (156.25 μg/ml and 31.25 μg/ml). To test the mosquito field samples and the Bora- Bora strain, concentrations were adjusted to include different percentages of survival/mortality, varying between 0 and 100%. Larvae were not fed during the 24 h of insecticide exposure. The concentrations tested for field samples were 1.5, 0.3, 0.15, 0.05, 0.015 μg/ml, and 1, 0.4, 0.2, 0.1, 0.04 μg/ml, for temephos and malathion, respectively (from insecticide solutions at 156.25 μg/ml); and 0.009, 0.0075, 0.006, 0.004, 0.0025 μg/ml (temephos), and 0.1, 0.075, 0.05, 0.04, 0.025 μg/ml (malathion) for the reference Bora Bora strain (from insecticide solutions at 31.25 μg/ml). The results were analyzed with XLSTAT-Biomed software (XLSTAT 2015 France) to determine the lethal concentration for 50% (LC50) and 95% (LC95) of the populations. Resistance ratios (RR50 and RR95) were calculated using LC50 and LC95 rates from Ae. aegypti-field-sampled populations compared with the LC50 and LC95 rates of the susceptible Bora Bora strain. The resistance levels were ranked into three categories: low resistance (RR50 < 5), medium or moderate resistance (5 ≤ RR50 ≤ 10), and high resistance (RR50 > 10) [17, 18].
Adult bioassays were performed following WHO recommendations adapted for this study [19, 20]. Deltamethrin-impregnated papers were used at three different concentrations: 0.05, 0.06, and 0.08%. For each mosquito population and each concentration, four replicates and two negative controls of 25 females each were tested. The 0.05% deltamethrin-impregnated papers were ordered from the WHO Malaysia provider, and the 0.06 and 0.08% deltamethrin papers were impregnated in the laboratory using the deltamethrin, PESTANAL® (Sigma-Aldrich, Inc. Missouri, USA) in a solution of two-thirds acetone (Sigma-Aldrich, Inc.) and one-third silicone (VWR® Pennsylvanie, USA) with Whatman® Grade 1 Qualitative Filtration Paper (Sigma-Aldrich, Inc. Missouri, USA). The impregnation was done following WHO procedures  using 2 ml of insecticide solution to impregnate one 12 × 15 cm paper that was dried for one night and stored at 4 °C. Adult mosquitoes were exposed by tarsal contact for 1 h to the 0.06 and 0.08% deltamethrin-impregnated papers. The knockdown (KD) (after 1 h of exposure) and mortality rates (after 24 h) were determined. A kinetic of KD rate during 2 h (by accounting the number of KD mosquitoes every 5 min) was also conducted for each population using 0.05% impregnated papers to estimate the KDT50 (Time necessary to affect 50% of mosquitoes) using XLSTAT software (dose-effect option). The KRR50 is the ratio between KDT50 of mosquitoes from the field and KDT50 of mosquitoes from the Bora Bora strain. KRR50 was scaled as follows: KRR50<1 = susceptible, 1≤KRR50<10 = low resistance, 10≤KRR50≤30 = moderate resistance, 30<KRR50<100 = high resistance, and KRR50≥100 = very high resistance .
Knockdown resistance (Kdr) genotyping
Primer sequence (5' – 3')
The allele-specific real-time quantitative polymerase chain reaction (PCR) using SYBR® Green dye (Applied Biosystems®, Californie, USA) was done on an Applied Biosystems® 7500 thermal cycler (Californie, USA). The amplification consisted of a 95 °C 3 min holding stage, 40 cycles at 95 °C for 15 s, a 60 °C 31 s cycling stage, and a melt curve stage. For the PCR, a 15 μl solution comprising 7.5 μl of SYBR® Green PCR Master Mix, 0.4 μl of 10 μM reverse primer, 0.2 μl of 10 μM of both forward primers (C1534-f and F1534-f, see Table 2), 3.7 μl of H2O, and 3 μl of genomic DNA was made.
A Ile1016/Ile1016 homozygote (resistant) has a single peak at 76 °C; a Val1016/Val1016 homozygote (susceptible) has double peaks at 79 °C and 83 °C; and a Val1016/Ile1016 heterozygote has triple peaks at 76 °C, 79 °C, and 83 °C. A Cys1534/Cys1534 homozygote (resistant) has a single peak at 82 °C; a Phe1534/Phe1534 (susceptible) homozygote has a single peak at 78 °C; and a Phe1534/Cys1534 heterozygote has double peaks at 78 °C and 82 °C. Lastly, the distribution of the three genotypes for both mutations and allele frequencies were calculated.
Gene expression study
The 12 genes used in relative quantitative PCR 
Name in this study
Carboxyl/cholinesterase alpha esterase
Carboxyl/cholinesterase alpha esterase
60S ribosomal protein L8
40S ribosomal protein S7
Chloride channel protein 2
The same DNA that was used for Kdr PCR was also used for the gene amplification study. Relative quantitative PCR was performed for five of the eight genes, i.e. CCEae3a, 014614, CYP6M11, CYP9J23, and GSTe2. The PCR was done as previously described in the Gene expression study section above. Results show genes relative quantification for field mosquito populations compared to Bora Bora strain and using both 007808 and 005950 housekeeping genes (DataAssist™ v3.01 software).
Statistical analyses were performed with STATISTICA 8 software (Statsoft Inc., Oklahoma, USA) and XLSTAT 2015 software. STATISTICA was used to perform ANOVA Kruskal-Wallis tests (to compare adult KD/mortality rates regarding deltamethrin and amplification gene ratios); Kruskal-Wallis Z tests (to compare larval mortality rates regarding temephos and malathion, adult mortality rates regarding deltamethrin, and expression and amplification genes ratios); and Mann–Whitney U tests (to compare gene expression ratios of Bora Bora strain and each population). XLSTAT was used to perform log-probit logistical regression to investigate the dose-effect relation regarding insecticide larval tests data, as well as to estimate KDT values using adult kinetic data. XLSTAT was also used to perform Spearman’s rank correlation tests (between expression and amplification ratios) and the principal component analysis (PCA).
Resistance status of Ae. aegypti larvae from Guadeloupe and Saint Martin, compared to the reference Bora Bora strain, to malathion and temephos
LC50 95% CI
LC95 95% CI
LC50 95% CI
LC95 95% CI
The differences in the levels of resistance observed between Ae. aegypti populations for temephos were tested for statistical significance using the Kruskal-Wallis Z test, which highlighted significantly higher resistance levels among mosquitoes from Les Abymes and Baie-Mahault, compared to Saint-François and Saint Martin (P < 0.001). In addition, significantly higher resistance levels were also observed in the three mosquito populations from Anse-Bertrand, Saint-François, and Deshaies, compared to Saint Martin West (P = 0.000 for Anse-Bertrand, P = 0.027 for Saint-François, and P = 0.000 for Deshaies). Finally, the mosquitoes from Saint Martin West had a significantly higher mortality rate than all mosquito populations from Guadeloupe.
For malathion, the results for all field mosquito populations show some resistance in comparison with the susceptible Bora Bora strain (see Table 4), with the RR50 ranging from 1.7 for Anse-Bertrand and Deshaies, to 4.4 for Les Abymes. The RR95 varied from 2.0 for Deshaies to 3.9 for Les Abymes and Baie-Mahault.
The differences in the levels of resistance between Ae. aegypti populations were tested for statistical significance using the Kruskal-Wallis Z test. The results show higher mortality rates suggesting lower resistance for the mosquito populations of Anse-Bertrand and Deshaies compared to all other field populations (P < 0.01) except for Saint-François (P = 0.325 for Anse-Bertrand and P = 0.227 for Deshaies). The mosquito population of Saint-François showed a higher mortality rate than those from Les Abymes (P = 0.001) and Baie-Mahault (P = 0.002), which suggests that these two latter populations have a higher resistance to malathion.
Resistance status of Ae. aegypti females from Guadeloupe and Saint Martin to deltamethrin
Deltamethrin concentration/time of exposure
KRR50 (ci 95%)
The mortality rates after 24 h following a 1-h exposure to a 0.06% deltamethrin concentration were significantly different between mosquito populations from Anse-Bertrand (0.64) and Baie-Mahault (0.90) (Kruskal-Wallis Z test, P = 0.040). The mortality rates after 24 h following a 1-h exposure to a 0.08% deltamethrin concentration were significantly different, with mosquito populations from Saint-François and Les Abymes exhibiting the lowest mortality rates, thus suggesting higher resistance levels (ANOVA Kruskal-Wallis: P < 0.0001).
KRR50 values were variable, ranging from 8.0 (Les Abymes) to 28.1 (Anse-Bertrand), indicating a low to moderate resistance level for deltamethrin when compared to the susceptible Bora Bora strain. Two levels of resistance were detected in the six field mosquito populations. The mosquito population of Les Abymes had the lowest KRR50 (8.0), suggesting a low resistance level to deltamethrin. All others mosquito populations had a KRR50 of between 12 and 28, and are considered moderately resistant.
Allele frequencies for the V1016I and F1534C Kdr mutations for each Ae. aegypti population
Mutation V1016I ƒ[I]
Mutation V1016I ƒ[V]
Mutation F1534C ƒ[C]
Mutation F1534C ƒ[F]
Detoxification enzyme levels
The cytochrome P450 monooxygenase, 014614, was overexpressed in the six mosquito populations compared to the susceptible Bora Bora strain, and ranged from 4.3 for Anse-Bertrand to 21.8 for Saint-François, with significant differences in expression ratios between these two locations (Kruskal-Wallis Z test: P = 0.011). The cytochrome P450 monooxygenase, CYP6BB2, was found to be overexpressed only in the Saint-François and Deshaies populations with similar ratios of 3.0 and 2.5, respectively (Mann–Whitney U test: P = 0.512). The cytochrome P450 monooxygenase, CYP6M11, was over-transcribed for mosquito populations of Saint Martin, Baie-Mahault, and Saint-François, with ratios of 2.5, 3.4, and 4.0, respectively, and similar transcription levels (ANOVA Kruskal-Wallis: P = 0.201). The cytochrome P450 monooxygenase, CYP9J23, was significantly overexpressed for all mosquito populations with the exception of Saint Martin, with ratios ranging from 2.0 (Deshaies) to 11.5 (Saint-François). The transcription levels were significantly higher for Saint-François than for Deshaies (Kruskal-Wallis Z test: P = 0.026).
The carboxyl/cholinesterase 3A (CCEea3a) was over-transcribed in the six mosquito populations, with expression ratios ranging from 6.02 for Saint Martin to 46.5 for Saint-François. Mosquito populations of Saint-François (46.5-fold), Deshaies (32.7-fold), and Baie-Mahault (20.0-fold) exhibited higher levels of over-transcription than the rest of the populations (ANOVA Kruskal-Wallis: P = 0.011).
The glutathione S-transferase 2 (GSTe2) was over-transcribed only in the mosquito populations of Les Abymes and Saint-François with similar ratios of 2.3 and 2.4, respectively (Mann–Whitney U test: P = 0.218).
CCEae6a and the cytochrome P450 monooxygenase, CYP9J28, were not found to be over-transcribed in any of the populations when compared to the susceptible Bora Bora strain.
Gene amplification study
The cytochrome P450 monooxygenase 014614 gene was amplified in all Ae. aegypti populations, with ratios ranging from 20.7 for Saint-Martin to 48.9 for Baie-Mahault. Significant differences in 014614 gene amplification were found between Saint Martin and both Baie-Mahault (Kruskal-Wallis Z test: P = 0.001) and Saint-François populations (Kruskal-Wallis Z test: P = 0.005). The cytochrome P450 monooxygenase, CYP6M11, was significantly amplified in the mosquito populations of Les Abymes, Baie-Mahault, and Deshaies, with ratios of 2.8, 2.2, and 2.0 respectively (ANOVA Kruskal-Wallis test: P = 0.036). The cytochrome P450 monooxygenase, CYP9J23, was amplified similarly in all Ae. aegypti populations with ratios of 5.2 to 7.0 (ANOVA Kruskal-Wallis test: P = 0.329). CCEae3a was amplified with ratios of >19 in all Ae. aegypti populations. Significant differences in amplification levels were observed between mosquito populations of Anse-Bertrand with a ratio of 19.2 and both Baie-Mahault and Les Abymes populations with ratios of 60.4 (Kruskal-Wallis Z test: P = 0.003) and 40.0 (Kruskal-Wallis Z test: P = 0.013), respectively. GSTe2 was similarly amplified in mosquito populations from Baie-Mahault (4.5-fold), Deshaies (3.4-fold), Anse-Bertrand (2.3-fold), and Les Abymes (2.1-fold) (ANOVA Kruskal-Wallis test: P = 0.115).
Correlation between overexpression and gene amplification of detoxification genes
Correlations between mean values of the expression and amplification ratios were investigated using the Spearman’s rank correlation test in XLSTAT software for each gene and all mosquito populations studied (see Additional file 3), however, no significant correlation was found. Nevertheless, the highest levels of association between gene amplification and expression ratios were found for the cytochrome P450 monooxygenases CYP9J23 and 014614 (r = 0.543; P = 0.297), followed by CCEae3a (r = 0.200; P = 0.714). Finally, the cytochrome CYP6M11 and GSTe2 showed the weakest associations between gene amplification and expression levels (r = −0.143; P = 0.803).
Association analysis between phenotypes, Kdr mutations, and detoxification enzymes regarding insecticide resistance
The first three PCA axes allow to represent 76.77% of the starting information with 34.38, 26.54, and 19.85%, respectively. The expression of the 014614 gene was mainly represented on the first PCA axe (Cos2 = 0.811) and was found to be significantly correlated to those of CYP6BB2 and CCEae3a genes (r = 0.870 and 0.853, respectively). The amplification ratio of CCEae3a was mainly represented on the second PCA axe (Cos2 = 0.901) and principally correlated with the GSTe2 amplification ratio (r = 0.808) and RR50 for temephos (r = 0.769). The CYP9J23 amplification ratio was mainly represented on the third axe (Cos2 = 0.633), and was found to be positively correlated with the 014614 and GSTe2 amplification ratios (r = 0.712 and r = 0.680, respectively). Furthermore, the expression ratios of CYP9J23 and CYP6M11 genes were significantly correlated (r = 0.850). The allele frequencies of resistant Kdr 1534C and 1016I mutations were positively correlated (r = 0.611), but only the 1534C resistant mutant allele was correlated with deltamethrin resistance KRR50 (r = 0.611). Temephos resistance was significantly correlated with the amplification ratio of the CYP6M11 gene (r = 0.849). Finally, the mortality rate after 24 h following a 1-h exposure to 0.08% deltamethrin was significantly and negatively correlated with the expression ratio of the GSTe2 gene (r = −0.960).
The larval tests carried out on Ae. aegypti mosquito populations from Guadeloupe and Saint Martin showed that these mosquitoes have high resistance to temephos and low resistance to malathion. Tests done with adult mosquitoes showed moderate resistance levels to deltamethrin for KRR50s.
The chemical larvicide temephos was used for about 40 years, from the late 1960s to 2010, and resistance was reported as early as 1984 (Pasteur Institute archives). For adult vector control, malathion was used from the late 1960s to 2009, and deltamethrin from the early 1980s until the present day (Pasteur Institute archives). The mosquito populations of Guadeloupe and Saint Martin are reported to be resistant (weakly) to malathion for the first time in Guadeloupe (in this study), but are known to be resistant to deltamethrin since the late 1990s (Pasteur Institute and ARS archives). The long duration of use of temephos as a larvicide in Guadeloupe and Saint Martin, and the rapid increase of deltamethrin resistance may be explained by cross-resistance between these two molecules .
For temephos, the results of our study (i.e. RR50 between 8.9 and 33) are similar to those obtained for Ae. aegypti populations from the region, including in Martinique , Cuba, Venezuela, Costa Rica, Panama, Nicaragua, and Jamaica . The detoxification enzymes commonly associated with organophosphates resistance are carboxylesterases (COEs) [28–30]. The gene expression study performed on adult mosquitoes revealed an overexpression for CCEae3a. Furthermore, the PCA showed a positive correlation between RR50 for temephos and CCEae3a amplification (r = 0.769). These results agree with previous studies demonstrating the implication of CCEae3a upregulation on larval temephos resistance . In 1999, the higher resistance levels found for Ae. aegypti mosquito populations  were likely a consequence of the strong selective pressure exerted by the use of temephos as the main larvicide at the time especially due to its low price, efficacy and remanence. Since 2010, only Bti has been used in Guadeloupe as a larvicide, which explains the decrease in resistance levels and the lower resistance ratios found for temephos in this study.
For malathion, the resistance levels observed on larvae are low, with the RR50 ranging from 1.7 to 4.4 for all mosquito populations tested, when compared to the Bora Bora strain. Although malathion has been used for more than 40 years to control adult mosquitoes in Guadeloupe, the mosquito larvae did not develop a strong resistance to this molecule. This observation can be explained by the alternate use of malathion and deltamethrin since the late 1980s/early 1990s in Guadeloupe (Pasteur Institute archives), which may have delayed the emergence of higher resistance levels. Furthermore, the use of malathion in Guadeloupe has been intensified mostly during epidemic outbreaks. Another factor that may have mitigated the resistance level for malathion in mosquito populations is the possible dependency on a specific carboxylesterase mechanism with no cross-resistance with temephos. Our results on resistance levels to malathion are in agreement with recent studies conducted on Ae. aegypti larvae from Martinique , Cuba, Jamaica, Panama, Costa Rica , and adult mosquitoes from Venezuela .
Deltamethrin has been used since the late 1980s/early 1990s and is currently the only adulticide authorized in European countries, including the French overseas territories of Guadeloupe and Saint Martin. The observed knockdown deltamethrin resistance ratios (8≤KRR50≤28) suggest that the Ae. aegypti populations of Guadeloupe and Saint Martin are more resistant than those of Martinique . On the contrary, when compared to French Guiana , Ae. aegypti populations from Guadeloupe were found to be less resistant in terms of KD (93–100% for Guadeloupe versus 27–37% for French Guiana) and mortality rates (64–90% for Guadeloupe versus 14–30% for French Guiana) obtained after exposure to 0.06% deltamethrin. Although our results contrast with those obtained by Dusfour and colleagues in 2015  (in which lower KD and mortality rates, 37 and 42%, respectively, were obtained for the Ae. aegypti populations from Baie-Mahault), our study and the study of Dusfour and colleagues in 2015 show that the Ae. aegypti populations from Guadeloupe have lower resistance levels to deltamethrin than those from French Guiana.
It has been previously reported that COEs, cytochrome P450 monooxygenases, and GSTs are involved in pyrethroids resistance and metabolization [35–39]. The adult gene expression study performed on eight detoxification enzymes revealed a significant overexpression of the genes GSTe2, CCEae3a, and the cytochrome P450 monooxygenases 014614, CYP9J23, CYP6M11, and CYP6BB2. The PCA revealed a significant and negative correlation between the GSTe2 expression ratio and the mortality rates obtained 24 h after a 1-h exposure to 0.08% deltamethrin. The overexpression of GSTe2 in Ae. aegypti resistant to DichloroDiphenylTrichloroethane (DDT) and permethrin had previously been reported in Thailand .
Kdr mutations are known to be strongly associated with organochloride and pyrethroid resistance for many mosquito species including Ae. aegypti [41, 42]. Kdr genotyping of the six mosquito populations studied revealed a highly resistant allele frequency of V1016I mutation ranging from 85 to 96%. These results are in agreement with those found for mosquito populations from Martinique and French Guiana, where frequencies ranged from 87 to 97% and from 65 to 92%, respectively [26, 33]. The mutant allele frequency of F1534C was also extremely high in populations from Guadeloupe, ranging from 90 to 98%. Similar frequencies of these mutations have been found in Mexico and Grand Cayman islands [43, 44]. A V1016I mutation has often been associated with permethrin resistance , while a F1534C mutation has been associated with a permethrin and DDT resistance . The PCA revealed a positive correlation between the deltamethrin resistance ratio and the mutant allele frequency of F1534C (r = 0.611), but a low correlation with V1016I (r = 0.262). Similar observations have already been reported in 2013 for the mosquito populations from Martinique regarding the V1016I mutation .
In conclusion, our results highlight that deltamethrin resistance levels seem to be more linked to Kdr mutation 1534C and GSTe2 at high concentrations, rather than other detoxification enzyme mechanisms. In terms of temephos and malathion, resistance levels were found to be more related to detoxification mechanisms involving esterase (CCEae3a) and cytochrome (CYP6M11) enzymes rather than Kdr mutations.
The levels of resistance of Ae. aegypti populations of the Guadeloupe and Saint Martin islands to the different types of insecticides are likely a consequence of vector control and agricultural activities. In Guadeloupe, lindane (ɣHCH) and chlordecone, two organochlorine pesticides, have been used in banana cultures to eliminate weevils in 1965–1974 and 1972–1993, respectively . Furthermore, the principal active products of herbicides used in local sugar cane cultures are generally organophosphates and carbamates [http://e-phy.agriculture.gouv.fr/], while certain insecticides used in horticulture contain pyrethroids such as deltamethrin. Today, the most abundant pesticide in water in the French overseas departments is organochlorine [http://www.statistiques.developpement-durable.gouv.fr/lessentiel/ar/246/211/contamination-globale-eaux-souterraines-pesticides.html]. The long history of pesticide use in agriculture in Guadeloupe with the same classes of molecules applied for vector control could thus have contributed to the type of insecticide resistance in Ae. aegypti populations. This phenomenon has been commonly described for malaria vectors [48–50].
Finally, gene amplification is commonly associated with an overexpression of detoxification enzymes [13, 51, 52]. In our study, no significant correlation was found between expression and amplification ratios for any of the five genes tested, however, the strongest associations were recorded for CYP9J23 and 014614 (r = 0.543).
In the present study, the populations of Ae. aegypti mosquitoes—the main vector of dengue, CHIKV, and ZIKV—from the Guadeloupe and Saint Martin islands were found to have different resistance levels to several insecticides, with a high resistance to the former larvicide temephos, moderate resistance to the adulticide deltamethrin, and low resistance to the former adulticide malathion. The resistance levels are associated with strong frequencies of F1534C and V1016I Kdr mutations, as well as the over-transcription of CCEae3a, GSTe2, and four cytochromes P450 (014614, CYP9J23, CYP6M11, CYP6BB2) genes.
Mosquito resistance to deltamethrin is a serious challenge for vector control authorities since this product is the only adulticide authorized in the French overseas territories, in particular during epidemics. In addition, deltamethrin has also been used in agriculture for decades. In this context, there is an urgent need to carefully choose the molecules used for vector control and agricultural activities, as well as reduce their utilization in order to limit the chance of resistance. It is also urgent to identify alternative molecules that are effective and respectful of the environment that could be used in a vector control strategy based on a rotation and/or combination of insecticides (which leads to resistance development slowing down), as recommended by the WHO .
Our knowledge of the molecular basis of Ae. aegypti insecticide resistance has contributed to real progress in the past few years, but it is still not developed enough either to be used as a diagnostic tool for identifying resistance without the necessity to develop heavy testing procedures or understanding how the resistance mechanisms are developing. The present study elucidated that it could be interesting to: i) assess the activity of some families of detoxification enzymes using biochemical/enzymatic assays; ii) use synergists such as Triphenyl phosphate (COEs inhibitor), piperonyl butoxide (P450s inhibitor), or ethacrynic acid or diethyl maleate (GSTs inhibitor) in order to specify enzyme families and mechanisms involved in resistance ; iii) segregate tested mosquitoes and conduct quantitative PCR to understand induced gene expression and thus to specify the role of some genes individually; iv) add more detoxification genes in the analysis (i.e. GSTe7, CYP6Z) ; and v) investigate more Kdr mutations as I1011V which has already been associated to pyrethroid resistance in Cuba and Mexico , as well as I1011M and G923V in Brazil . A better understanding of the resistant genes and their roles is needed in order to improve vector control strategies, mitigate insecticide resistance levels, and prevent and control the emergence and spread of Ae. aegypti-borne diseases.
Agence Régionale de Santé
Bacillus thuringiensis var. israelensis
First laboratory generation mosquitoes
Second laboratory generation mosquitoes
Principal component analysis
Polymerase chain reaction
Programme de Surveillance, d’Alerte et de Gestion des epidémies de dengue
World Health Organization
We thank the regional health agency of Guadeloupe for its collaboration regarding mosquito sampling. We thank the vector control center of Martinique for graciously offering the Bora-Bora susceptible reference strain and for the training of DG on insecticide resistance tests. We thank the Pasteur Institute of French Guiana for graciously offering Kdr PCR controls and the International Division of Pasteur Institute for supporting DG’s training in Alpine ecology laboratory (LECA Grenoble, France).
This study was financially supported by the Fonds Européen de Développement Economique et Régional project FED-1-1.4 32932 and by the ARS through the collaboration convention no° 2014-140116. DG was supported by La Région Guadeloupe, Le Fond Social Européen, and the University of Antilles.
Availability of data and materials
DG, AVR, and FF conceived the study. DG, CD, AG and AVR collected biological materials, performed the experiments, and analyzed the data. DG, AVR and FF wrote the paper. CR and JG provided technical support for biological material collection. JPD, FFC and TG provided PCR primers, protocols and technical support for the screening of insecticide resistance mechanisms. All authors read and approved the final paper.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Informed consent was obtained from the owners of the houses that were investigated for mosquito collections. This consent was managed by the Vector Control Agency of Guadeloupe, through the Ethical Committee of the French Ministry of Health, because all mosquito collections were done in collaboration with their agents.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Leparc-Goffart I, Nougairede A, Cassadou S, Prat C, de Lamballerie X. Chikungunya in the Americas. Lancet. 2014;383:514.View ArticlePubMedGoogle Scholar
- Rollé A, Schepers K, Cassadou S, Curlier E, Madeux B, Hermann-storck C, et al. Severe sepsis and septic shock associated with Chikungunya virus infection, Guadeloupe, 2014. Emerg Infect Dis. 2016;22:891–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Zanluca C, de Melo VCA, Mosimann ALP, dos Santos GIV, dos Santos CND, Luz K. First report of autochthonous transmission of Zika virus in Brazil. Mem Inst Oswaldo Cruz. 2015;110:569–72.View ArticlePubMedPubMed CentralGoogle Scholar
- Petersen LR, Jamieson DJ, Powers AM, Honein MA. Zika virus. N Engl J Med. 2016;374(16):1552–63.View ArticlePubMedGoogle Scholar
- CIRE ANTILLES-GUYANE. Emergence du virus Zika aux Antilles Guyane. Situation épidémiologique. INVS, Point épidémiologique du 28 avril 2016: 16/2016.
- Larrieu S, Cassadou S, Rosine J, Chappert JL, Blateau A, Ledrans M, et al. Lessons raised by the major 2010 dengue epidemics in the French West Indies. Acta Trop Elsevier BV. 2014;131:37–40.View ArticleGoogle Scholar
- Gottschamel J, Lössl A, Ruf S, Wang Y, Skaugen M, Bock R, et al. Production of dengue virus envelope protein domain III-based antigens in tobacco chloroplasts using inducible and constitutive expression systems. Plant Mol Biol. 2016;91(4–5):497–512.View ArticlePubMedGoogle Scholar
- Weaver S, Costa F, Garcia-Blanco M, Ko A, Ribeiro G, Saade G, et al. Zika virus: history, emergence, biology, and prospects for control. Antivir Res. 2016;130:69–80.View ArticlePubMedGoogle Scholar
- Debboun M, Strickman D. Insect repellents and associated personal protection for a reduction in human disease. Med Vet Entomol. 2013;27:1–9.View ArticlePubMedGoogle Scholar
- Rosine J. Résistance d’Aedes aegypti et de Culex pipiens quinquefasciatus aux insecticides organophospohorés, biologiques et aux pyréthrinoïdes en Martinique et en Guadeloupe. Diplôme d’Etudes Approfondies. Univ. Pierre et Marie Curie (Paris VI) 51p. 1999.
- Settimi L, Orford R, Davanzo F, Hague C, Desel H, Pelclova D, et al. Development of a new categorization system for pesticides exposure to support harmonized reporting between EU Member States. Env Int. 2016;91:332–40.View ArticleGoogle Scholar
- ANSES. Recherche d’insecticides potentiellement utilisables en lutte anti-vectorielle. Avis de l’Anses. Rapp. d’expertise Collect. Ed. Sci. Saisine n°2009-SA-0338. ISBN 978-2-11-128697-9. Anses Ed. 2011;158p.
- Faucon F, Dusfour I, Gaude T, Navratil V, Boyer F, Chandre F, et al. Identifying genomic changes associated with insecticide resistance in the dengue mosquito Aedes aegypti by deep targeted sequencing. Genome Res. 2015;25:1347–59.View ArticlePubMedPubMed CentralGoogle Scholar
- Hemingway J, Hawkes NJ, McCarroll L, Ranson H. The molecular basis of insecticide resistance in mosquitoes. Insect Biochem Mol Biol. 2004;34:653–65.View ArticlePubMedGoogle Scholar
- Marcombe S, Poupardin R, Darriet F, Reynaud S, Bonnet J, Strode C, et al. Exploring the molecular basis of insecticide resistance in the dengue vector Aedes aegypti: a case study in Martinique Island (French West Indies). BMC Genomics. 2009;10:494.View ArticlePubMedPubMed CentralGoogle Scholar
- World Health Organization. Guidelines for laboratory and field testing of mosquito larvicides, Geneva. 2005 (WHO/CDS/WHOPES/GCDPP/2005.13).
- Bisset JA, Rodríguez MM, Ricardo Y, Ranson H, Pérez O, Moya M, et al. Temephos resistance and esterase activity in the mosquito Aedes aegypti in Havana, Cuba increased dramatically between 2006 and 2008. Med Vet Entomol. 2011;25:233–9.View ArticlePubMedGoogle Scholar
- Mazzarri M, Georghiou G. Characterization of resistance to organophosphate, carbamate, and pyrethroid insecticides in field populations of Aedes aegypti from Venezuela. J Am Mosq Control Assoc. 1995;11:315–22.PubMedGoogle Scholar
- World Health Organization. Guidelines for testing mosquito adulticides for indoor residual spraying and treatment of mosquito nets, Geneva. 2006 (WHO/CDS/NTD/WHOPES/GCDPP/2006.3).
- World Health Organization. Test procedures: for insecticide resistance monitoring in malaria vector mosquitoes, Geneva. 2013.
- Khan HAA, Akram W, Shehzad K, Shaalan EA. First report of field evolved resistance to agrochemicals in dengue mosquito, Aedes albopictus (Diptera: Culicidae), from Pakistan. Parasit Vectors BioMed Central Ltd. 2011;4:146.View ArticleGoogle Scholar
- Yanola J, Somboon P, Prapanthadara L. A novel point mutation in the Aedes aegypti voltage-gated sodium channel gene associated with permethrin resistance. 2nd Int. Conf. Dengue Dengue Haemorhagic Fever, Oct. 2008. p. 15–7
- Yanola J, Somboon P, Walton C, Nachaiwieng W, Prapanthadara L. A novel F1552/C1552 point mutation in the Aedes aegypti voltage-gated sodium channel gene associated with permethrin resistance. Pestic. Biochem. Physiol. Elsevier Inc.; 2010
- Yanola J, Somboon P, Walton C, Nachaiwieng W, Somwang P, Prapanthadara L. 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.View ArticlePubMedGoogle Scholar
- Rodriguez MM, Bisset J, Ruiz M, Soca A. Cross-resistance to pyrethroid and organophosphorus insecticides induced by selection with temephos in Aedes aegypti (Ditera:Culicidae) from Cuba. J Med Entomol. 2002;39:882–8.View ArticlePubMedGoogle Scholar
- Marcombe S, Mathieu RB, Pocquet N, Riaz MA, Poupardin R, Sélior S, et al. Insecticide resistance in the dengue vector Aedes aegypti from martinique: distribution, mechanisms and relations with environmental factors. PLoS One. 2012;7:e30989.View ArticlePubMedPubMed CentralGoogle Scholar
- Rodriguez MM, Bisset JA, Fernandez D. Levels of insecticide resistance and resistance mechanisms in Aedes aegypti from some Latin American countries. J Am Mosq Control Assoc. 2007;23:420–9.View ArticlePubMedGoogle Scholar
- Saavedra-Rodriguez K, Strode C, Flores AE, Garcia-Luna S, Reyes-Solis G, Ranson H, et al. Differential transcription profiles in Aedes aegypti detoxification genes following temephos selection. Insect Mol Biol. 2014;23:199–215.View ArticlePubMedGoogle Scholar
- Paiva MHS, Lovin DD, Mori A, Melo-santos MAV, Severson DW, Ayres FJ. Identification of a major Quantitative Trait Locus determining resistance to the organophosphate temephos in the dengue vector mosquito Aedes aegypti. Genomics Elsevier BV. 2015;107:40–8.View ArticleGoogle Scholar
- Poupardin R, Srisukontarat W, Yunta C, Ranson H. Identification of carboxylesterase genes implicated in temephos resistance in the dengue vector Aedes aegypti. PLoS Negl Trop Dis. 2014;8:e2743.View ArticlePubMedPubMed CentralGoogle Scholar
- Etienne M. Etude de la bioécologie d’Aedes aegypti à la Martinique en relation avec l’épidémiologie de la dengue. 2006
- Alvarez LC, Ponce G, Oviedo M, Lopez B, Flores AE. Resistance to malathion and deltamethrin in Aedes aegypti (Diptera : Culicidae) from western Venezuela. J Med Entomol. 2013;50:1031–9.View ArticlePubMedGoogle Scholar
- Dusfour I, Thalmensy V, Gaborit P, Issaly J, Carinci R, Girod R. Multiple insecticide resistance in Aedes aegypti (Diptera: Culicidae) populations compromises the effectiveness of dengue vector control in French Guiana. Mem Inst Oswaldo Cruz. 2011;106:346–52.View ArticlePubMedGoogle Scholar
- Dusfour I, Zorrilla P, Guidez A, Issaly J, Girod R. Deltamethrin resistance mechanisms in Aedes aegypti populations from three French overseas territories worldwide. PLoS Negl Trop Dis. 2015;9:1–17.View ArticleGoogle Scholar
- Somwang P, Yanola J, Suwan W, Walton C, Lumjuan N, Prapanthadara P, Somboon LA. Enzymes-based resistant mechanism in pyrethroid resistant and susceptible Aedes aegypti strains from northern Thailand. Parasitol Res. 2011;109:531–7.View ArticlePubMedGoogle Scholar
- Stevenson BJ, Pignatelli P, Nikou D, Paine MJI. Pinpointing P450s associated with pyrethroid metabolism in the dengue vector, Aedes aegypti: developing new tools to combat insecticide resistance. PLoS Negl Trop Dis. 2012;6:e1595.View ArticlePubMedPubMed CentralGoogle Scholar
- Chandor-Proust A, Bibby J, Régent-Kloeckner M, Roux J, Guittard-Crilat E, Poupardin R, et al. The central role of mosquito cytochrome P450 CYP6Zs in insecticide detoxification revealed by functional expression and structural modelling. Biochem J. 2013;455:75–85.View ArticlePubMedPubMed CentralGoogle Scholar
- David J-P, Ismail HM, Chandor-proust A, Paine MJI. Role of cytochrome P450s in insecticide resistance : impact on the control of mosquito-borne diseases and use of insecticides on Earth. Philos Trans R Soc Lond B Biol Sci. 2013;368:1–12.View ArticleGoogle Scholar
- Maia RT, Nadvorny D. Molecular docking of Anopheles gambiae and Aedes aegypti glutathione S-Transferases Epsilon 2 (GSTE2) against usnic acid: an evidence of glutathione conjugation. Brazilian Arch Biol Technol. 2014;57:689–94.View ArticleGoogle Scholar
- Lumjuan N, McCarroll L, Prapanthadara LA, Hemingway J, Ranson H. Elevated activity of an Epsilon class glutathione transferase confers DDT resistance in the dengue vector, Aedes aegypti. Insect Biochem Mol Biol. 2005;35:861–71.View ArticlePubMedGoogle Scholar
- Donnelly MJ, Corbel V, Weetman D, Wilding CS, Williamson MS, Black WC. Does kdr genotype predict insecticide-resistance phenotype in mosquitoes? Trends Parasitol. 2009;25:213–9.View ArticlePubMedGoogle Scholar
- Saavedra-Rodriguez K, Strode C, Flores Suarez A, Fernandez Salas I, Ranson H, Hemingway J, et al. Quantitative trait loci mapping of genome regions controlling permethrin resistance in the mosquito Aedes aegypti. Genetics. 2008;180:1137–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Vera-Maloof FZ, Saavedra-Rodriguez K, Elizondo-Quiroga AE, Lozano-Fuentes S, Black IV WC. 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Harris AF, Rajatileka S, Ranson H. Pyrethroid resistance in Aedes aegypti from Grand Cayman. Am J Trop Med Hyg. 2010;83:277–84.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Marcombe S, Paris M, Paupy C, Bringuier C, Yebakima A, Chandre F, et al. Insecticide-driven patterns of genetic variation in the dengue vector Aedes aegypti in Martinique Island. PLoS One. 2013;8:e77857.View ArticlePubMedPubMed CentralGoogle Scholar
- Cabidoche Y-M, Lesueur-Jannoyer M. Long term pollution of soils in French West Indies : how to manage chlordecone contamination? Innov Agron. 2011;16:117–33.Google Scholar
- Philbert A, Lyantagaye SL, Nkwengulila G. A review of agricultural pesticides use and the selection for resistance to insecticides in malaria vectors. Adv Entomol. 2014;2:120–8.View ArticleGoogle Scholar
- Diabaté A, Baldet T, Chandre F, Akogbéto M, Guiguemde TR, Darriet F, et al. The role of agricultural use of insecticides in resistance to pyrethroids in Anopheles gambiae s.l. in Burkina Faso. Am J Trop Med Hyg. 2002;67:617–22.PubMedGoogle Scholar
- Nkya T, Poupardin R, Laporte F, Akhouayri I, Mosha F, Magesa S, et al. Impact of agriculture on the selection of insecticide resistance in the malaria vector Anopheles gambiae : a multigenerational study in controlled conditions. Parasit Vectors. 2014;7:480.PubMedPubMed CentralGoogle Scholar
- Bariami V, Jones CM, Poupardin R, Vontas J, Ranson H. Gene amplification, ABC transporters and cytochrome P450s: unraveling the molecular basis of pyrethroid resistance in the dengue vector, Aedes aegypti. PLoS Negl Trop Dis. 2012;6:e1692.View ArticlePubMedPubMed CentralGoogle Scholar
- Li X, Schuler M, Berenbaum M. Molecular mechanisms of metabolic resistance to synthetic and natural xenobiotics. Annu Rev Entomol. 2007;52:231–53.View ArticlePubMedGoogle Scholar
- WHO. The technical basis for coordinated action against insecticide resistance: preserving the effectiveness of modern malaria vector control. Glob. Malar. Program. WHO Headquarters. Draft Meet. Rep. 2010.
- Bisset JA, Marin R, Rodriguez MM, Severson DW, Ricardo Y, French L, et al. Insecticide resistance in two Aedes aegypti (Diptera: Culicidae) strains from Costa Rica. J Med Entomol. 2013;50:352–61.View ArticlePubMedGoogle Scholar
- Brengues C, Hawkes NJ, Chandre F, McCarroll L, Duchon S, Guillet P, et al. Pyrethroid and DDT cross-resistance in Aedes aegypti is correlated with novel mutatuions in the voltage-gated sodium channel gene. Med Vet Entomol. 2003;17:87–94.View ArticlePubMedGoogle Scholar