Relationship between insecticide resistance profiles in Anopheles gambiae sensu lato and agricultural practices in Côte d’Ivoire
Parasites & Vectors volume 16, Article number: 270 (2023)
Insecticide-based malaria vector control is increasingly undermined due to the development of insecticide resistance in mosquitoes. Insecticide resistance may partially be related to the use of pesticides in agriculture, while the level and mechanisms of resistance might differ between agricultural practices. The current study aimed to assess whether phenotypic insecticide resistance and associated molecular resistance mechanisms in Anopheles gambiae sensu lato differ between agricultural practices.
We collected An. gambiae s.l. larvae in six sites with three different agricultural practices, including rice, vegetable and cocoa cultivation. We then exposed the emerging adult females to discriminating concentrations of bendiocarb (0.1%), deltamethrin (0.05%), DDT (4%) and malathion (5%) using the standard World Health Organization insecticide susceptibility test. To investigate underlying molecular mechanisms of resistance, we used multiplex TaqMan qPCR assays. We determined the frequency of target-site mutations, including Vgsc-L995F/S and Vgsc-N1570Y, and Ace1-G280S. In addition, we measured the expression levels of genes previously associated with insecticide resistance in An. gambiae s.l., including the cytochrome P450-dependent monooxygenases CYP4G16, CYP6M2, CYP6P1, CYP6P3, CYP6P4, CYP6Z1 and CYP9K1, and the glutathione S-transferase GSTe2.
The An. gambiae s.l. populations from all six agricultural sites were resistant to bendiocarb, deltamethrin and DDT, while the populations from the two vegetable cultivation sites were additionally resistant to malathion. Most tested mosquitoes carried at least one mutant Vgsc-L995F allele that is associated with pyrethroid and DDT resistance. In the cocoa cultivation sites, we observed the highest 995F frequencies (80–87%), including a majority of homozygous mutants and several in co-occurrence with the Vgsc-N1570Y mutation. We detected the Ace1 mutation most frequently in vegetable-growing sites (51–60%), at a moderate frequency in rice (20–22%) and rarely in cocoa-growing sites (3–4%). In contrast, CYP6M2, CYP6P3, CYP6P4, CYP6Z1 and CYP9K1, previously associated with metabolic insecticide resistance, showed the highest expression levels in the populations from rice-growing sites compared to the susceptible Kisumu reference strain.
In our study, we observed intriguing associations between the type of agricultural practices and certain insecticide resistance profiles in the malaria vector An. gambiae s.l. which might arise from the use of pesticides deployed for protecting crops.
Vector control is the main strategy for controlling malaria and has shown success in Africa over many years. Unfortunately, vectors are becoming increasingly resistant to the insecticides used in public health. This threatens the efficacy of long-lasting insecticidal nets (LLINs) and indoor residual spraying (IRS), which are the main tools for controlling malaria mosquitoes . The classical active compounds used for these tools are pyrethroids, organochlorines, organophosphates and carbamates. However, other active ingredients have recently been repurposed, including neonicotinoids and pyrroles . Resistance to commonly used insecticides has been reported in malaria vectors in several sub-Saharan African countries [3,4,5,6,7,8,9].
In addition to behavioural changes, different physiological mechanisms are involved in insecticide resistance. The most important ones described in African malaria vectors are target-site resistance leading to alterations of the insecticide target sites, preventing the binding of the insecticide ; metabolic resistance, characterised by changes in insect enzyme systems leading to rapid detoxification or sequestration of insecticides ; and cuticular resistance that reduces the amount of insecticide penetrating the insect .
Target-site mutations in the voltage-gated sodium channel (Vgsc) are associated to pyrethroid and dichlorodiphenyltrichloroethane (DDT) resistance and are also known as knockdown resistance (kdr) mutations [6, 13, 14]. In Anopheles gambiae sensu lato, kdr is predominantly conferred by two mutations at the same codon position Vgsc-L995F/S [15, 16]. An additional mutation, called ‘super kdr ’ (Vgsc-N1570Y), co-occurs with the L995F allele and increases resistance further . In West Africa and particularly in Côte d’Ivoire, kdr is among the most commonly reported resistance mechanisms, with the Vgsc-995F allele predominating [6, 18,19,20]. Another point mutation in the acetylcholinesterase gene that causes a glycine to serine substitution (Ace1-G280S) and confers resistance to carbamate and organophosphate  has also been found in Culex pipiens quinquefasciatus  and An. gambiae [23, 24] mosquitoes from Côte d’Ivoire. While the kdr and Ace1 target-site mutations are involved in insecticide resistance, they alone do not explain the highly resistant phenotypes observed in mosquitoes .
The other major physiological mechanism of insecticide resistance in mosquitoes is the overexpression of metabolic enzymes that detoxify or sequester insecticides . Metabolic resistance includes, among others, enzymes from three major families of enzymes: the cytochrome P450-dependent monooxygenases (P450s or CYPs), carboxylesterase and glutathione S-transferases (GSTs). Overexpression of several detoxification genes related to insecticide resistance has been detected in field mosquito populations in Côte d’Ivoire, mostly from rice-cultivating areas [14, 26,27,28]. In 2014, Edi et al.  showed that some genes belonging to the CYP6 P450 family metabolise pyrethroids and other insecticides in field mosquitoes from Tiassalé in southern Côte d’Ivoire. The involvement of P450s in the resistance of malaria vectors from the same locality has been repeatedly demonstrated in subsequent studies using synergist bioassays , quantitative polymerase chain reaction (qPCR) and RNA sequencing .
Although insecticide resistance in malaria vectors is attributed to the use of insecticides in public health, an increasing number of studies suggest that the use of pesticides in agriculture also contributes to the selection of insecticide resistance in mosquitoes [1, 18, 29, 30]. Growing human populations, particularly in Africa, exert increasing pressure on agricultural productivity, leading to intensified use of pesticides . Indeed, importations of pesticides between 2005 and 2015 nearly tripled in the West African region, particularly in the three largest pesticide markets of Côte d'Ivoire, Ghana and Nigeria . Furthermore, the uncontrolled and improper use of agrochemicals can lead to the development of insecticide resistance in non-target insects, including malaria vectors breeding in agricultural areas, since the compounds used to control crop pests often have the same active ingredients and molecular targets as those used in public health. . Therefore, farming irrigation greatly increases the risk of malaria for nearby communities as well as the development of insecticide resistance by pesticides intended to control crop pests [33,34,35]. Cross-resistance to public health insecticides in mosquitoes is now a real obstacle for the current vector control methods adopted by the World Health Organization (WHO) Global Plan for Insecticide Resistance Management (GPIRM) . These methods are primarily based on IRS and LLINs, and require a better understanding of the association between agricultural practices and insecticide resistance.
Previous studies have reported widespread use of pesticides in different types of agriculture and have shown that vector resistance to insecticides varies from one site to another in Côte d’Ivoire . Although both target-site resistance and upregulation of P450 genes have been described in An. gambiae s.l. in the country, the association between different agricultural practices and insecticide resistance has not been studied. Therefore, assessing the current status of malaria vector resistance from different agricultural practices becomes a necessity to facilitate more effective planning of control strategies based on the type of agriculture. Here, we assessed the phenotypic resistance to discriminating concentrations of insecticides using the standard WHO insecticide susceptibility test in An. gambiae s.l. collected from rice, vegetable and cocoa cultivation areas in western Côte d'Ivoire. We further measured the frequency of target-site resistance alleles and characterised the expression levels of metabolic genes previously associated with insecticide resistance in the country.
Larval collections were carried out in six study sites in Côte d’Ivoire with three different agricultural practices, including cocoa, rice and vegetable cultivation. Each type of cultivation included two sites in the study. The sites with irrigated rice cultivation were in Agboville (latitude: 5.935496°, longitude: −4.223084°) and Tiassalé (latitude: 5.904263°, longitude: −4.826142°) located in the forest zones of Côte d’Ivoire. The vegetable sites were Azaguié (latitude: 5.633333°, longitude: −4.083333°) and Dabou (latitude: 5.316667°, longitude −4.383333°), also located in the forest zone, while the cocoa sites were Issia (latitude: 6.487614°, longitude: −6.583677°) and Soubré (latitude: 5.786623°, longitude: −6.589017°), which are characterised by evergreen forests (Fig. 1).
The climate in the study areas is humid and equatorial , characterised by four seasons: (i) a long rainy season that brings heavy rains from May to June; (ii) a short rainy season with rains from August to September in the southern part and from August to October in the more western part; (iii) a short dry season from October to November; and (iv) a main dry season from December to April. The average temperature varies between 21 °C and 33 °C.
To control dipteran and lepidopteran insects damaging the crops, mainly pyrethroids and neonicotinoids are deployed in the rice fields of Agboville and Tiassalé and the vegetable fields of Azaguié and Dabou . Both agricultural practices are carried out throughout the year in Côte d’Ivoire. Issia and Soubré are located in the largest cocoa production region of Côte d’Ivoire. Cocoa production is the main economic activity in the western part of the country and one of the major pillars of Côte d’Ivoire’s economy. In cocoa-growing areas, neonicotinoids combined with pyrethroids are the insecticides predominantly used to protect cocoa fields against mirid infestations that are considered to be the most important pest problem in cocoa cultivation [18, 38]. In our study sites, the LLINs distributed in the irrigated rice and vegetable areas contain deltamethrin, while those distributed in the cocoa area contain alpha-cypermethrin .
Larval collections and insecticide susceptibility bioassay
We collected An. gambiae s.l. larvae in southern Côte d’Ivoire (i.e. Agboville, Azaguié, Dabou and Tiassalé) from June to July in 2018 and in the western Côte d’Ivoire (i.e. Issia, Soubré) in June 2019. In Agboville and Tiassalé, we collected larvae in irrigated rice fields, while in Dabou and Azaguié, they were collected in water pits in the vegetable fields. In Soubré and Issia, we made collections in water puddles in cocoa fields (Fig. 1). Once back to Abidjan, we reared the larvae to the adult stage in the insectary of the Centre Suisse de Recherches Scientifiques en Côte d’Ivoire (CSRS) under standard conditions of 26 ± 2 °C, 75 ± 10% relative humidity and a 12 h/12 h light/dark photoperiod. Larvae from each site were reared separately and fed each morning with 0.075 g of TetraMin® fish food powder (Tetra, Melle, Germany). Emerging adults had access to a 10% honey solution.
To determine the phenotypic insecticide susceptibility to the four classic insecticide classes commonly used in public health, we performed the standard WHO insecticide susceptibility test  with 2- to 5-day-old, non-blood-fed adult female mosquitoes that had emerged from the field-collected larvae. We tested their susceptibility against the WHO discriminating concentrations of bendiocarb (0.1%), deltamethrin (0.05%), DDT (4%) and malathion (5%) on treated filter papers sourced from WHO. To determine the knockdown ratio, we conducted bioassays with adults taken from an insecticide susceptible An. gambiae sensu stricto. Kisumu colony and exposed them to deltamethrin (0.05%) and DDT (4%). For each combination of insecticide and field site, we exposed six batches of 20–25 females including four batches exposed to insecticide-impregnated filter papers and two batches serving as negative controls and exposed to control papers containing only the insecticide carrier oil. During the 1 h exposure time, we recorded how many mosquitoes were knocked down at 5-min intervals. After an hour of exposure, we transferred the mosquitoes back to the holding tubes and allowed them to feed on 10% honey solution ad libitum, while the delayed mortality was recorded 24 h post-exposure.
Preparation of samples and extraction of nucleic acids for molecular analysis
We killed the mosquitoes that were still alive 24 h post-exposure in absolute ethanol and then blotted away any excess ethanol with a paper towel. We then gently transferred the mosquitoes by batches of 10 to 100 individuals into 1.5 ml microcentrifuge tubes that contained 0.7–1.4 ml of RNAlater® (Ambion, Inc., Austin, TX, USA) depending on the number of mosquitoes. We kept RNAlater tubes with mosquitoes overnight at 4 °C to allow for thorough penetration of the tissue. The following day, the excess RNAlater® was removed and the tubes were stored at −20 °C until further processing for DNA and RNA extraction.
For the extraction of total RNA and DNA, we randomly picked 50 RNAlater®-preserved individuals from the controls that had not been exposed to insecticides in the bioassays and processed them using the MagnaMedics magnetic bead-based kit (MagnaMedics GmbH, Aachen, Germany). We ground the mosquitoes individually in 1.5-ml microcentrifuge tubes by adding 200 µl TE buffer (10 mM tris(hydroxymethyl) aminomethane hydrochloride [Tris–HCl], 1 mM ethylenediaminetetraacetic acid [EDTA], pH 8.0) per tube and using a pestle, driven by a hand-held tissue grinder. Then, we added 150 µl lysis buffer, mixed everything by vortexing for 15 s and incubated the mixture for 10 min at room temperature while vortexing the tubes every 2 min for 15 s. Following the incubation, we spun down the non-lysed mosquito debris by centrifugation at 16,000×g for 2 min. Next, we transferred the clear supernatant into a new 1.5-ml tube, added 20 µl magnetic beads and 440 µl binding buffer, and vortexed the mixture for 15 s. We incubated the mixture again for 10 min at room temperature while vortexing the tubes in between as described above for the lysis step. The tubes were placed on a magnetic rack for 2 min, allowing for the magnetic beads to sediment, and then discarded the supernatant. We washed the remaining beads twice by adding 200 µl wash buffer, vortexed the tubes and let the mixture incubate for 1 min at room temperature before placing it for 2 min back on the magnetic rack and discarding the supernatant. Finally, we extracted the nucleic acids by adding 180 µl elution buffer, vortexed the tubes and incubated the mixture for 10 min in a water bath set at 50 °C and vortexed in between as above. After the tubes were removed from the water bath, we vortexed them again, spun them down and placed the tubes on the magnetic rack for 2 min. Finally, we collected the supernatant that now contained the purified DNA and RNA, transferred the solution to new 1.5-ml tubes and stored them at −80 °C.
To identify the sibling species of the An. gambiae s.l. complex, we performed two TaqMan® multiplex qPCR assays [41, 42] with modifications to the original protocols as described in Wipf et al. . In brief, the first assay differentiates Anopheles coluzzii and An. gambiae s.s. as a group (Ag+) from Anopheles bwambae, Anopheles melas, Anopheles merus and Anopheles quadriannulatus (Aq+) and from Anopheles arabiensis (Aa+). The second assay distinguishes between An. coluzzii (former molecular M-form) and An. gambiae s.s. (former molecular S-form) based on the SINE 200 X6.1 locus that is fixed in An. coluzzii and absent in An. gambiae s.s. We used the common primers designed by Santolamazza et al.  and the probes described in Wipf et al. . We ran the second assay on the DNA from the samples that turned out to be Ag+ in the first assay.
In addition to the species identification assays, we applied diagnostic TaqMan® qPCR assays to detect the target-site kdr mutations of the voltage-gated sodium channel (i.e. Vgsc-L995F/S and Vgsc-N1570Y) and the acetylcholinesterase mutation Ace1-G280S, following the protocol of Bass et al.  with the adaptations described in Mavridis et al. .
In addition to the diagnostic qPCRs for species identification and the detection of target-site mutations, we measured the expression levels of genes that have previously been associated with insecticide resistance in An. gambiae s.l., including the cytochrome P450-dependent monooxygenases CYP4G16, CYP6M2, CYP6P1, CYP6P3, CYP6P4, CYP6Z1 and CYP9K1, and the glutathione S-transferase GSTe2. As a reference for the overall gene expression, we measured additionally the expression levels of the housekeeping gene encoding for the ribosomal protein S7 (RPS7). Again, we used a series of messenger RNA (mRNA)-specific TaqMan® reverse transcription qPCRs (RT-qPCRs) that were developed by Mavridis et al. . Instead of pooling RNA or DNA, we measured the gene expression levels separately for each individual with 50 individuals per field site and 50 individuals from the insecticide-susceptible An. gambiae s.s. Kisumu colony. This susceptible strain from Kenya is kept in the insectary as a control for insecticide susceptibility bioassays.
We ran the qPCR reactions in volumes of 10 µl, containing 1 µl of template nucleic acid extract and 9 µl of master mix comprising primers and probes at final concentrations as published previously . The master mix reagents were supplied by Fast-Track Diagnostics (FTD, Esch-sur-Alzette, Luxembourg). All reactions were performed on a C1000 Touch/CFX96™ Real-Time PCR System (Bio-Rad Laboratories, Hercules, CA, USA) in 96-well plates (Sarstedt, Nümbrecht, Germany; catalogue number: 72.1980.202). The thermal cycle parameters were 15 min for the reverse transcription at 50 °C, RTase inactivation and initial denaturation at 95 °C for 3 min, followed by 40 cycles of denaturation at 95 °C for 3 s and annealing/extension steps at 60 °C for 30 s.
For the insecticide susceptibility assay, we scored both the immediate knockdown at 5-min intervals up to 1 h and the delayed mortality rates at 24 h post-exposure. Survival analysis was performed to visualise pyrethroid and DDT knockdown time, comparing the knockdown effect on the field mosquito populations to the one on the insecticide susceptible Kisumu colony. For the interpretation of the mortality rates in terms of phenotypic resistance status, we followed the WHO criteria : a mortality rate below 90% indicates resistance; a mortality rate equal to or above 98% indicates susceptibility; and a mortality rate between 90 and 97% is suggestive of possible resistance that needs to be confirmed.
The expression levels of the genes of interest relative to the internal control (RPS7) were calculated using the comparative threshold cycle (CT) method, also referred to as the 2−ΔCT method, where ΔCT = CT gene of interest—CT internal control . To assess whether gene expression levels differ between mosquito populations, we ran a linear regression model for each gene with 2−ΔCT as the dependent variable. We used the susceptible Kisumu strain as the reference. We set the level of significance at α = 0.05 and adjusted the P-values for multiple testing using the Bonferroni adjustment method .
We ran all data analysis in R version 4.0.3  using RStudio version 1.3.1093 . We used the ‘tidyverse’ R package for data tidying and visualisation . The R packages ‘survival’ and ‘survminer’  were used to plot the survival curves and to calculate the Kaplan–Meier estimates .
Insecticide susceptibility assays
The knockdown rates for DDT and deltamethrin varied significantly between the sites (DDT: P-value < 0.001; deltamethrin: P-value < 0.0001) but were substantially lower in all field populations when compared to the insecticide-susceptible Kisumu colony (P-value < 0.0001) (Fig. 2) confirming the insecticide resistance observed in the delayed 24 h mortality rates (Fig. 3). Deltamethrin exposure knocked down mosquitoes from all tested field and lab populations more rapidly than DDT (Fig. 2). Among the field populations, the population from the Dabou vegetable site showed the fastest and highest knockdown rate with both deltamethrin and DDT although endpoint knockdown rates 1 h post-exposure did not exceed 25% and 12%, respectively. The lowest knockdown rate with deltamethrin was observed in the Issia cocoa field population. In contrast, the majority of the Kisumu mosquitoes were already knocked down after 10 min exposure to deltamethrin and all of them after 40 min (Fig. 2). The knockdown probability with DDT was the lowest (0%) in Soubré cocoa field populations and highest (12%) in Dabou vegetable field populations (P-value < 0.0001) (Fig. 2). In comparison, half of the susceptible Kisumu mosquitoes were knocked down by DDT after 25 min and all after 50 min exposure time (Fig. 2).
All mosquito populations from the three agricultural practices were resistant to deltamethrin, DDT and bendiocarb. However, the populations from the rice and cocoa sites were still susceptible to malathion (Fig. 3). The Kisumu reference colony was fully susceptible to DDT and deltamethrin. The 24 h mortality rates in the non-exposed control groups were below 5% in each bioassay and thus we did not correct for control mortalities.
Across all sites, the qPCR assays confirmed the presence of both An. gambiae s.s. and An. coluzzii with a 96% predominance of An. coluzzii (n = 295). All mosquitoes from the rice fields of Agboville and Tiassalé as well as those from the vegetable farms of Azaguié and Dabou were identified as An. coluzzii, except for one An. coluzzii/An. gambiae s.s. hybrid from Azaguié. While An. coluzzii was the predominant species in the cocoa plantations of Issia and Soubré, we identified 12% and 10% of the specimens as An. gambiae s.s., respectively.
We identified several target-site mutations, including the kdr loci Vgsc-L995F and Vgsc-N1570Y and the acetylcholinesterase mutation Ace1-G280S, while the allelic frequency varied by site and agricultural practice (Fig. 4). In contrast, we did not detect the mutant Vgsc-L995S allele. The mutation Vgsc-L995F that is associated with pyrethroid and DDT resistance was the most frequent allele across all sites (Fig. 4). We detected high Vgsc-L995F allelic frequencies in the cocoa-growing sites of Issia and Soubré, whereas we found moderate frequencies in the vegetable and rice sites. The Vgsc-N1570Y mutation was present at low frequencies in Azaguié, Issia, Soubré and Tiassalé and undetected in Agboville and Dabou (Fig. 4). The highest frequencies were found in both cocoa-growing sites.
We found the Ace1-G280S mutation associated with carbamate and organophosphate resistance at the highest allelic frequencies in both the Dabou (60%) and Azaguié (51%) vegetable-growing sites but mostly in Dabou, while it was under 25% in the rice field populations of Agboville and Tiassalé, and under 5% in the cocoa fields of Issia and Soubré (Fig. 4).
The cytochrome P450-dependent monooxygenases CYP4G16, CYP6M2, CYP6P1, CYP6P3, CYP6P4, CYP6Z1 and CYP9K1 were significantly overexpressed compared to Kisumu with variations across the different agricultural sites (Figs. 5 and 6). Cocoa-growing sites showed overall low expression levels, while rice-growing sites had increased expression levels and vegetable-growing sites showed a mixed picture (Fig. 6). The three P450s CYP9K1, CYP6Z1 and CYP6M2 were significantly overexpressed in all six field sites while being most notably overexpressed in the rice-growing sites of Agboville and Tiassalé followed by the vegetable site of Dabou (Fig. 6). The highest fold changes were detected for CYP9K1 in Agboville (5.3-fold up) and Dabou (5.2-fold up). Interestingly, while CYP4G16 was upregulated in Dabou, the same gene was not significantly differentially expressed in the other vegetable-growing site Azaguié and the cocoa-growing site Issia.
Regardless of the prevailing agricultural practice, An. gambiae s.l. populations from all six Ivorian sites were resistant to deltamethrin, DDT and bendiocarb, whereas only populations from vegetable-growing areas were additionally resistant to malathion. While kdr mutations altering the pyrethroid and organochlorine target site appeared as the driving insecticide resistance mechanism in cocoa sites, we identified high expression levels of insecticide-metabolising P450 enzymes in the irrigated rice sites. In the two vegetable sites, we found both kdr and Ace1 target-site mutations associated with phenotypic resistance but the sites showed contrasting expression patterns of metabolic genes.
In Côte d’Ivoire—as in other African countries—it has been suggested that the use of pesticides in agriculture is one of the causes indirectly related to insecticide resistance in malaria vectors, although few studies have focused on the relationship between agriculture and insecticide resistance [29, 30, 53, 54]. Most of the previous studies conducted in Côte d’Ivoire have shown a very high trend of resistance to pyrethroids and DDT and also resistance to carbamates and organophosphates is progressively increasing [27, 28, 55]. It was found in a monitoring study that pyrethroids are predominant among the insecticides used in vegetable, rice and cocoa fields while vectors from these agricultural practices showed resistance to DDT, deltamethrin and bendiocarb . In 2016, Chouaïbou et al.  detected malathion resistance in malaria vectors collected from vegetable fields in Dabou and found that organophosphates accounted for 9% of insecticides used in this site compared to only 2% in the rice fields of Tiassalé, where mosquitoes were susceptible to malathion, like in our study. Thus, agricultural pesticides in mosquito breeding habitats may favour the selection of insecticide resistance in An. gambiae s.l. mosquitoes against compounds that were not used in vector control interventions.
In the search for the mechanisms causing the observed phenotypic resistance, we found a very high allelic frequency of the Vgsc-L995F mutation associated to pyrethroid and DDT resistance in cocoa areas. Indeed, Côte d’Ivoire being the first world cocoa producer, its western part is in the process of becoming the new zone of high cocoa production of the country . Thus, pesticides such as pyrethroids are intensively applied in cocoa farming to protect them against crop pests . These activities and the deployment of alpha-cypermethrin-treated nets during interventions of the national malaria control programme in this region  likely exerted enormous selection pressure, leading to the development of pyrethroid resistance in the vectors. Moreover, we detected at moderate frequency the presence of the Vgsc-1570Y allele that amplifies the effects of the Vgsc-L995F mutation. To our knowledge, this is the first time that the Vgsc-N1570Y mutation has been detected in Issia and Soubré cocoa fields. The findings of this study show that the phenotypic resistance observed in the cocoa region is probably associated with the high frequency of the Vgsc-L995F and Vgsc-N1570Y mutation, likely confirming the crucial role of the kdr mutation in conferring pyrethroid and DDT resistance.
We observed resistance to carbamate and organophosphate insecticides in An. gambiae s.l. populations from vegetable fields of Dabou and Azaguié. This may be indicative of intensive use of carbamate and organophosphate, which can pollute mosquito breeding sites in the cultivated environment. This hypothesis is supported by the study of Chouaïbou et al.  that found carbamate and organophosphate residues in mosquito breeding sites and in the soil in Dabou vegetable fields as well as phenotypic resistance in An. gambiae s.l. to bendiocarb and malathion. Reassuringly, we found the highest frequencies of the insensitive acetylcholinesterase Ace1-G280S gene in mosquitoes from vegetable areas, confirming the cross-resistance observed between carbamate and organophosphate insecticides. In 2016, Chouaïbou et al.  detected malathion resistance in malaria vectors collected from vegetable fields in Dabou and found that organophosphates accounted for 9% of insecticides used in this site compared to only 2% in the rice fields of Tiassalé, where mosquitoes were susceptible to malathion as in our study. Therefore, the presence of Ace1 mutation conferring cross-resistance to these insecticides represents an important threat for carbamate and organophosphate-based vector control strategies. The use of newly developed insecticides with different modes of action and their combination with older classes in large-scale interventions is urgently needed to manage insecticide resistance.
Finally, gene expression analysis revealed that several P450s were overexpressed in rice fields relative to Kisumu. We found that P450s CYP9K1, CYP6M2, CYP6Z1, CYP6P4 and CYP6P3 were the most upregulated among the genes detected in rice-growing sites and vegetable-growing sites. In fact, frequent exposure of larvae to agricultural pollutants and several xenobiotics in the water of rice and vegetable fields could induce metabolic stress resulting in insecticide detoxification in the insect. Recent studies have shown that recurrent exposure of An. gambiae s.l. and Aedes aegypti larvae to agrochemicals induced significant tolerance to insecticides due to a stimulation of multiple genes responsible for target-site resistance and metabolic resistance, including P450 genes [56, 57]. Thus, the practice of rice growing and vegetable growing several times a year and the intensive use of insecticides in these crops, as well as the xenobiotics present in the crop water pits favouring metabolic gene expression, may be why P450 mechanisms are more strongly selected in rice-growing sites and in the vegetable-growing site of Dabou than the cocoa-growing sites. The overexpression of several P450s revealed in both the rice-growing sites and vegetable-growing site of Dabou, alongside the target-site mutations, may explain the low DDT, deltamethrin and bendiocarb mortalities observed. We found that CYP6M2 was upregulated in all sites, but the gene was mostly overexpressed in Tiassalé compared to Kisumu. In fact, Tiassalé field mosquitoes were found resistant to deltamethrin, whereas they were still susceptible to malathion despite the presence of Ace1 mutations, which might be an indication for negative cross-resistance as discussed in Wipf et al. . The evidence that P450 enzymes (CYP6M2) can confer negative cross-resistance has also been provided directly in an in vivo study with African malaria vectors by Adolfi et al. . Such increased susceptibility to malathion may have a positive impact on insecticide resistance management—especially for the improvement of malaria vector control tools .
Although we observed some associations between insecticide resistance mechanisms and agricultural practice, we acknowledge that several other factors may have influenced the studied outcome. A major caveat is that the majority of our field sites with the same agricultural type were geographically closer to each other than to those of another type. Additionally, the larval collections in the cocoa-growing areas were done a year later than in the rice- and vegetable-growing areas. Both limitations above may have led to similar environmental factors—in addition to those stemming from agriculture—influencing resistance selection in malaria vectors. For example, the LLINs distributed in the different geographical areas were treated with different pyrethroid insecticides . Future studies would greatly benefit from measuring the strength of resistance to different insecticides using intensity or dose–response bioassays and including the novel insecticides prequalified by WHO: broflanilide, chlorfenapyr and clothianidin .
Although recent studies have shown overexpression of detoxification genes in some of our study sites but mostly in pooled mosquito samples, this study went a step further by measuring the gene expression levels of each individual separately, strengthening our results and providing entomological baseline data for rice-, vegetable- and cocoa-growing areas.
The insecticide resistance observed in all cultivation areas would constitute an obstruction to the various control strategies that are essentially based on LLINs and IRS if nothing is done to reverse the current situation. The involvement of other environmental factors such as particular pollutants and waste from industries established around crop areas needs to be thoroughly investigated to determine whether these factors contribute to the selection pressure of the metabolic and target site resistance mechanism in malaria vectors around farms. In conclusion, the present study revealed intriguing associations between agricultural practice and the type of resistance mechanisms in malaria vectors that merit further exploration. Therefore, it is suggested that national malaria control programmes collaborate more closely with the agricultural sector to jointly develop integrated risk management and vector control strategies involving farmers.
Availability of data and materials
All relevant data are included within the paper. Raw data are available from the corresponding author upon reasonable request.
- Ace1 :
- An. gambiae :
- An. gambiae s.l.:
Anopheles gambiae sensu lato
- An. gambiae s.s.:
Anopheles gambiae sensu stricto
Centre Suisse de Recherches Scientifiques en Côte d’Ivoire
- CT :
Global Plan for Insecticide Resistance Management
Indoor residual spraying
- Kdr :
Long-lasting insecticidal nets
Cytochrome P450-dependent monooxygenases
Reverse-transcription quantitative (or real-time) polymerase chain reaction
Tris(hydroxymethyl) aminomethane hydrochloride
- Vgsc :
Voltage-gated sodium channel
World Health Organization
Hemingway J. The role of vector control in stopping the transmission of malaria: threats and opportunities. Philos Trans R Soc Lond B Biol Sci. 2014;369:20130431.
WHO. WHO guidelines for malaria. Geneva: World Health Organization; 2023.
Antonio-Nkondjio C, Poupardin R, Tene BF, Kopya E, Costantini C, Awono-Ambene P, et al. Investigation of mechanisms of bendiocarb resistance in Anopheles gambiae populations from the city of Yaounde, Cameroon. Malar J. 2016;15:424.
Chouaibou MS. Increase in susceptibility to insecticides with aging of wild Anopheles gambiae mosquitoes from Côte d’Ivoire. BMC Infect Dis. 2012. https://doi.org/10.1186/1471-2334-12-214.
Diabate A, Baldet T, Chandre F, Akoobeto 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.
Edi CA, Koudou BG, Bellai L, Adja AM, Chouaibou M, Bonfoh B, et al. Long-term trends in Anopheles gambiae insecticide resistance in Cote d’Ivoire. Parasit Vectors. 2014;7:500.
Elanga-Ndille E, Nouage L, Ndo C, Binyang A, Assatse T, Nguiffo-Nguete D, et al. The G119S Acetylcholinesterase (Ace-1) target site mutation confers carbamate resistance in the major malaria vector Anopheles gambiae from Cameroon: a challenge for the coming IRS implementation. Genes. 2019;10:790.
Ngufor C, Fongnikin A, Rowland M, N’Guessan R. Indoor residual spraying with a mixture of clothianidin (a neonicotinoid insecticide) and deltamethrin provides improved control and long residual activity against pyrethroid resistant Anopheles gambiae sl in Southern Benin. PLoS ONE. 2017;12:e0189575.
Oxborough RM, N’Guessan R, Jones R, Kitau J, Ngufor C, Malone D, et al. The activity of the pyrrole insecticide chlorfenapyr in mosquito bioassay: towards a more rational testing and screening of non-neurotoxic insecticides for malaria vector control. Malar J. 2015;14:124.
Liu N. Insecticide resistance in mosquitoes: impact, mechanisms, and research directions. Annu Rev Entomol. 2015;60:537–59.
Poupardin R, Reynaud S, Strode C, Ranson H, Vontas J, David JP. Cross-induction of detoxification genes by environmental xenobiotics and insecticides in the mosquito Aedes aegypti: impact on larval tolerance to chemical insecticides. Insect Biochem Mol Biol. 2008;38:540–51.
Balabanidou V, Kampouraki A, MacLean M, Blomquist GJ, Tittiger C, Juarez MP, et al. Cytochrome P450 associated with insecticide resistance catalyzes cuticular hydrocarbon production in Anopheles gambiae. Proc Natl Acad Sci USA. 2016;113:9268–73.
Chandre F, Darriet F, Manguin S, Brengues C, Carnevale P, Guillet P. Pyrethroid cross resistance spectrum among populations of Anopheles gambiae s.s. from Côte d’Ivoire. J Am Mosq Control Assoc. 1999;15:53–9.
Wipf NC, Duchemin W, Kouadio FA, Fodjo BK, Sadia CG, Mouhamadou CS, et al. Multi-insecticide resistant malaria vectors in the field remain susceptible to malathion, despite the presence of Ace1 point mutations. PLoS Genet. 2022;18:e1009963.
Martinez-Torres D, Chandre F, Williamson MS, Darriet F, Bergé JB, Devonshire AL, et al. Molecular characterization of pyrethroid knockdown resistance (kdr) in the major malaria vector Anopheles gambiae s.s. Insect Mol Biol. 1998;7:179–84.
Ranson H, Jensen B, Vulule JM, Wang X, Hemingway J, Collins FH. Identification of a point mutation in the voltage-gated sodium channel gene of Kenyan Anopheles gambiae associated with resistance to DDT and pyrethroids. Insect Mol Biol. 2000;9:491–7.
Jones CM, Liyanapathirana M, Agossa FR, Weetman D, Ranson H, Donnelly MJ, et al. Footprints of positive selection associated with a mutation (N1575Y) in the voltage-gated sodium channel of Anopheles gambiae. Proc Natl Acad Sci USA. 2012;109:6614–9.
Fodjo BK, Koudou BG, Tia E, Saric J, N’Dri PB, Zoh MG, et al. Insecticides Resistance Status of An gambiae in Areas of Varying Agrochemical Use in Cote D’Ivoire. Biomed Res Int. 2018. https://doi.org/10.1155/2018/2874160.
Edi AVC, N’Dri BP, Chouaibou M, Kouadio FB, Pignatelli P, Raso G, et al. First detection of N1575Y mutation in pyrethroid resistant Anopheles gambiae in Southern Cote d’Ivoire. Wellcome Open Res. 2017;2:71.
Koffi AA, Alou LP, Adja MA, Koné M, Chandre F, N’Guessan R. Update on resistance status of Anopheles gambiae s.s. to conventional insecticides at a previous WHOPES field site, “Yaokoffikro”, 6 years after the political crisis in Côte d’Ivoire. Parasit Vectors. 2012. https://doi.org/10.1186/1756-3305-5-68.
Weill M, Lutfalla G, Mogensen K, Chandre F, Berthomieu A, Berticat C, et al. Comparative genomics: insecticide resistance in mosquito vectors. Nature. 2003;423:136–7.
Chandre F, Darriet F, Doannio JM, Rivière F, Pasteur N, Guillet P. Distribution of organophosphate and carbamate resistance in Culex pipiens quinquefasciatus (Diptera: Culicidae) in West Africa. J Med Entomol. 1997;34:664–71.
Edi CV, Koudou BG, Jones CM, Weetman D, Ranson H. Multiple-insecticide resistance in Anopheles gambiae mosquitoes Southern Cote d’Ivoire. Emerg Infect Dis. 2012;18:1508–11.
Ngufor C, Chouaïbou M, Tchicaya E, Loukou B, Kesse N, N’Guessan R, et al. Combining organophosphate-treated wall linings and long-lasting insecticidal nets fails to provide additional control over long-lasting insecticidal nets alone against multiple insecticide-resistant Anopheles gambiae in Côte d’Ivoire: an experimental hut trial. Malar J. 2014;13:396.
Machani MG, Ochomo E, Zhong D, Zhou G, Wang X, Githeko AK, et al. Phenotypic, genotypic and biochemical changes during pyrethroid resistance selection in Anopheles gambiae mosquitoes. Sci Rep. 2020;10:19063.
Chouaibou M, Zivanovic GB, Knox TB, Jamet HP, Bonfoh B. Synergist bioassays: a simple method for initial metabolic resistance investigation of field Anopheles gambiae s.l. populations. Acta Trop. 2014;130:108–11.
Edi CV, Djogbenou L, Jenkins AM, Regna K, Muskavitch MA, Poupardin R, et al. CYP6 P450 enzymes and ACE-1 duplication produce extreme and multiple insecticide resistance in the malaria mosquito Anopheles gambiae. PLoS Genet. 2014;10:e1004236.
Oumbouke WA, Pignatelli P, Barreaux AMG, Tia IZ, Koffi AA, Ahoua Alou LP, et al. Fine scale spatial investigation of multiple insecticide resistance and underlying target-site and metabolic mechanisms in Anopheles gambiae in central Cote d’Ivoire. Sci Rep. 2020;10:15066.
Chouaibou MS, Fodjo BK, Fokou G, Allassane OF, Koudou BG, David JP, et al. Influence of the agrochemicals used for rice and vegetable cultivation on insecticide resistance in malaria vectors in southern Cote d’Ivoire. Malar J. 2016;15:426.
Nkya TE, 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.
IFAD. Rural development report 2016. Rome: IFAD; 2016.
Haggblade S, Diarra A, Traoré A. Regulating agricultural intensification: lessons from West Africa’s rapidly growing pesticide markets. Develop Policy Rev. 2022;40:e12545.
Chouaibou MS, Chabi J, Bingham GV, Knox TB, N’Dri L, Kesse NB, et al. Increase in susceptibility to insecticides with aging of wild Anopheles gambiae mosquitoes from Côte d’Ivoire. BMC Infect Dis. 2012;12:214.
Koudou BG, Adja AM, Matthys B, Doumbia M, Cissé G, Koné M, et al. Agricultural activities and malaria transmission in two eco-epidemiological settings in central Côte d’Ivoire. Bull Soc Pathol Exot. 2007;100:124–6.
Chan K, Tusting LS, Bottomley C, Saito K, Djouaka R, Lines J. Malaria transmission and prevalence in rice-growing versus non-rice-growing villages in Africa: a systematic review and meta-analysis. Lancet Planet Health. 2022;6:e257–69.
WHO. Global plan for insecticide resistance management. Geneva: World Health Organization; 2012.
Yao Sadaiou Sabas B, Gislain Danmo K, Akoua Tamia Madeleine K, Jan B. Cocoa production and forest dynamics in Ivory Coast from 1985 to 2019. Land. 2020;9:524.
Wouter V, Luc ANGT, Alexandre K. Evaluation of cocoa mirid (Distantiella theobroma Dist. and Sahlbergella singularis Hagl.) control practices in Côte d’Ivoire. African J Agric Res. 2019;14:2006–16.
NMCP. Le Paludisme en Côte d'Ivoire 2020 http://www.pnlpcotedivoire.org/paludisme-en-cote-d-ivoire/.
WHO. WHO insecticide resistance test procedures. Geneva: World Health Organization; 2016.
Bass C, Williamson MS, Field LM. Development of a multiplex real-time PCR assay for identification of members of the Anopheles gambiae species complex. Acta Trop. 2008;107:50–3.
Santolamazza F, Mancini E, Simard F, Qi Y, Tu Z, Della Torre A. Insertion polymorphisms of SINE200 retrotransposons within speciation islands of Anopheles gambiae molecular forms. Malar J. 2008;7:163.
Bass C, Nikou D, Donnelly MJ, Williamson MS, Ranson H, Ball A, et al. Detection of knockdown resistance (kdr) mutations in Anopheles gambiae: a comparison of two new high-throughput assays with existing methods. Malar J. 2007;6:111.
Mavridis K, Wipf N, Muller P, Traore MM, Muller G, Vontas J. Detection and monitoring of insecticide resistance mutations in Anopheles gambiae: individual vs pooled specimens. Genes. 2018;9:479.
Mavridis K, Wipf N, Medves S, Erquiaga I, Muller P, Vontas J. Rapid multiplex gene expression assays for monitoring metabolic resistance in the major malaria vector Anopheles gambiae. Parasit Vectors. 2019;12:9.
Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008;3:1101–8.
Jafari M, Ansari-Pour N. Why, when and how to adjust your P values? Cell J. 2019;20:604–7.
R Core Team. R: A language and environment for statistical computing. Austria: R Foundation for Statistical Computing; 2020.
RStudio Team. RStudio: Integrated Development Environment for R. PBC, Boston, MA: RStudio; 2020.
Wickham H, Averick M, Bryan J, Chang W, McGowan L, François R, et al. Welcome to the Tidyverse. J Open Source Software. 2019;4:1686.
Alboukadel Kassambara [aut c, Marcin Kosinski [aut], Przemyslaw Biecek [aut], Scheipl Fabian [ctb]. Drawing Survival Curves using 'ggplot2'. 2021.
Stalpers LJA, Kaplan EL, Edward L. Kaplan and the Kaplan-Meier Survival Curve. BSHM Bulletin: J British Soc Hist Mathemat. 2018;33:109–35.
Antonio-Nkondjio C, Fossog BT, Ndo C, Djantio BM, Togouet SZ, Awono-Ambene P, et al. Anopheles gambiae distribution and insecticide resistance in the cities of Douala and Yaoundé (Cameroon): influence of urban agriculture and pollution. Malar J. 2011;10:154.
Nwane P, Etang J, Chouaїbou M, Toto JC, Koffi A, Mimpfoundi R, et al. Multiple insecticide resistance mechanisms in Anopheles gambiae s.l. populations from Cameroon, Central Africa. Parasites Vectors. 2013;6:41.
Sadia CG, Kouadio F-PA, Fodjo BK, Oyou SK, Beatrice A-GA, Koudou BG, et al. Short-term metabolic resistance inductive effect of different agrochemical groups on Anopheles gambiae mosquitoes. Wellcome Open Res. 2021;6:147.
Nkya TE, Akhouayri I, Poupardin R, Batengana B, Mosha F, Magesa S, et al. Insecticide resistance mechanisms associated with different environments in the malaria vector Anopheles gambiae: a case study in Tanzania. Malar J. 2014;13:28.
Riaz MA, Poupardin R, Reynaud S, Strode C, Ranson H, David JP. Impact of glyphosate and benzo[a]pyrene on the tolerance of mosquito larvae to chemical insecticides. Role of detoxification genes in response to xenobiotics. Aquat Toxicol. 2009;93:61–9.
Adolfi A, Poulton B, Anthousi A, Macilwee S, Ranson H, Lycett GJ. Functional genetic validation of key genes conferring insecticide resistance in the major African malaria vector, Anopheles gambiae. Proc Natl Acad Sci USA. 2019;116:25764–72.
Wanjala CL, Mbugi JP, Ototo E, Gesuge M, Afrane YA, Atieli HE, et al. Pyrethroid and DDT resistance and organophosphate susceptibility among Anopheles spp. Mosquitoes. Western Kenya Emerg Infect Dis. 2015;21:2178–81.
NMCP. Le paludisme en Côte d'Ivoire: National Malaria Control Programme; 2020 https://www.pnlpcotedivoire.org/paludisme-en-cote-d-ivoire/.
WHO. Prequalified vector control products: World Health Organization; 2023 https://extranet.who.int/pqweb/vector-control-products/prequalified-product-list.
We are grateful to all the farmers of Agboville, Tiassalé, Azaguié, Dabou, Issia and Soubré who gave us permission to collect larvae in their fields. We thank Laura Vavassori from Swiss TPH and Zoh Marius from CSRS for providing relevant criticism on the first version of the manuscript. Our acknowledgement goes to Sébastien Oyou Kere and Williams Adienin for helping with the larval collections. We are grateful to Dobri Laurent Didier, Akoupo Alain Sylvain and the late Nestor Kesse for their technical assistance at the CSRS insectary.
This project was funded by the European Union Horizon 2020 Framework Program (Grant Agreement No. 688207-DMC-MALVEC) and the Wellcome Trust .
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Kouadio, FP.A., Wipf, N.C., Nygble, A.S. et al. Relationship between insecticide resistance profiles in Anopheles gambiae sensu lato and agricultural practices in Côte d’Ivoire. Parasites Vectors 16, 270 (2023). https://doi.org/10.1186/s13071-023-05876-0