Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Re-Visiting Insecticide Resistance Status in Anopheles gambiae from Côte d'Ivoire: A Nation-Wide Informative Survey

Abstract

Insecticide resistance constitutes a major threat that may undermine current gain in malaria control in most endemic countries. National Malaria Control Programmes (NMCPs) need as much information as possible on the resistance status of malaria vectors and underlying mechanisms in order to implement the most relevant and efficient control strategy. Bioassays, biochemical and molecular analysis were performed on An. gambiae collected in six sentinel sites in Côte d'Ivoire. The sites were selected on the basis of their bioclimatic status and agricultural practices. An. gambiae populations across sites showed high levels of resistance to organochloride, pyrethroid and carbamate insecticides. The kdr and ace-1R mutations were detected in almost all sentinel sites with mosquitoes on the coastal and cotton growing areas mostly affected by these mutations. At almost all sites, the levels of detoxifying enzymes (mixed-function oxidases (MFOs), non-specific esterases (NSE) and glutathione-S-transferases (GSTs)) in An. gambiae populations were significantly higher than the levels found in the susceptible strain Kisumu. Pre-exposure of mosquitoes to PBO, an inhibitor of MFOs and NSEs, significantly increased mortality rates to pyrethroids and carbamates in mosquitoes but resistance in most cases was not fully synergised by PBO, inferring a residual role of additional mechanisms, including kdr and ace-1 site insensitivity. The large distribution of resistance in Côte d'Ivoire raises an important question of whether to continue to deploy pyrethroid-based long-lasting insecticidal nets (LLINs) and insecticide residual spraying (IRS) towards which resistance continues to rise with no guarantee that the level of resistance would not compromise their efficacy. Innovative strategies that combine insecticide and synergists in LLINs or spatially LLIN and an effective non-pyrethroid insecticide for IRS could be in the short term the best practice for the NMCP to manage insecticide resistance in malaria vectors in Côte d'Ivoire and other endemic countries facing resistance.

Background

Malaria vector control programmes rely largely on long-lasting insecticidal nets (LLINs) and indoor residual spraying (IRS) in sub-Saharan Africa [1]. For LLINs thirteen brands impregnated with only three pyrethroid insecticides (permethrin, deltamethrin and α-cypermethrin) are recommended by the World Health Organization Pesticide Scheme (WHOPES) [2]. WHOPES also recommends 12 insecticide formulations for IRS, involving four insecticide classes: organochlorides, organophosphates, carbamates and pyrethroids [3].

Unfortunately insecticide resistance mechanisms are spreading faster than thought across Africa [4]. This threatens to undermine the efficacy of IRS and LLINs that contain the four insecticide classes cited above. Four types of resistance mechanisms against public health insecticides have been identified: metabolic resistance, target site resistance, penetration resistance and behavioural resistance [4][6]. Nevertheless only metabolic and target site resistance mechanisms have been extensively investigated at both the phenotypic and genetic levels [7].

So far, all insecticides used in malaria vector control are neurotoxic. The acetylcholinesterase which hydrolyses the neurotransmitter, acetylcholine in the synaptic gap is the target of organophosphate and carbamate insecticides. The voltage-gated sodium channel that triggers impulse along the neuron membrane via discharge of Na+ is the target site of pyrethroid and organochlorine insecticides. The c-aminobutyric acid (GABA) is a neurotransmitter of the inhibitory synapses in the insect nervous system [8]. Non-silent point mutations that render these targets less sensitive to insecticides are respectively: the ace-1R, the kdr and the Rdl mutation [9][11]. Three main groups of enzymes are involved in metabolic resistance mainly through an overproduction process: carboxylesterases (efficient against organophosphate and carbamate insecticides), glutathione-S-transferases or GSTs (efficient against organophosphate, organochlorine, and pyrethroid insecticides) and cytochrome P450-dependent monooxygenases (efficient against most insecticide classes, frequently in conjunction with other enzymes). The overproduction might be the result of gene amplification increasing the gene's copy number or gene expression modification.

These mechanisms might have strong negative impact on the operational efficacy of both ITNs and IRS [12][14]. The study of forces driving these mechanisms evolution, their detection and their monitoring are crucial in order to implement the most relevant insecticide management strategy.

National Malaria Control Programmes (NMCP) in every endemic country have to revise and update their strategic plan every 3 or 5 years. In Côte d'Ivoire, the national strategy for malaria vector control relies essentially on the scaling up of systematic distribution of LLINs. Between 2007 and 2011 the NMCP distributed 9,671,246 LLINs to the targeted population of pregnant women and children under five years. Since 2011 the LLINs distribution campaign with 7,429,470 LLINs aims to scale the universal coverage up. The NMCP also plans to implement a vast IRS campaign to supplement universal coverage of LLINs [15].

In Cote d'Ivoire, there have been a number of studies investigating insecticide resistance mechanisms but they were all far from being representative because just two bioclimatic zones (guinean and wet savannah) were surveyed [16][20]. In order to implement the most relevant strategy in all bioclimatic zones, the NMCP developed and supported the present study to generate wider baseline data on resistance status in malaria vector populations in the country. The resistance levels to the four classes of insecticide so far deployed for ITNs and potentially for IRS in the country were investigated. The underlying mechanisms driving resistance were characterized.

Methods

Ethics statement

No specific permissions were required for the larvae sampling in urban areas because mosquito breeding sites were located in public areas. In contrast in rural areas, owners of the sampled cultivated areas gave their permission before the larvae collections. Field collections did not involve endangered or protected species.

Study sites

The study was conducted in six sentinel sites currently involved in malaria surveillance by the NMCP in Côte d'Ivoire. These sites were selected on the basis of bioclimatic parameters and agricultural practices in order to capture a broader range of ecological patterns found in Côte d'Ivoire (Figure 1 and Table 1). Where possible, mosquito larvae were sampled in breeding sites from both the urban and rural within sentinel sites:

thumbnail
Figure 1. Map of Côte d'Ivoire showing the localities in the different ecological zones where anopheline mosquitoes were collected.

https://doi.org/10.1371/journal.pone.0082387.g001

thumbnail
Table 1. Description of sampling sites, agricultural practices, rural/urban status, bioclimatic zones and the distribution of the molecular forms of An. gambiae s.s.

https://doi.org/10.1371/journal.pone.0082387.t001

  1. San Pedro and Abidjan, belonging to the Guinean bioclimatic zone were characterized by hyper-ombrophilous rainy forest. The average annual rainfall is between 1400 mm and 2000 mm, with two rainy seasons in a year (April to July and September to November) with an average temperature of 26±2°C according to records from the local wheather stations. The main agricultural practices include coffee, cocoa, rubber, palm and coconuts. Breeding sites investigated at San Pedro and three localities selected within the economic capital city, Abidjan, (i.e. Yopougon, Port-Bouët and Bingerville) were situated in urban or peri-urban areas.
  2. Man, is a city located in the Western zone of the country consisting of tropical rainy forest with hills and mountains. It is characterized by a single rainy season starting in February up to November. The average annual rainfall ranged from 1600 mm to 2500 mm and the annual temperature averaged 24.5°C. Coffee and cocoa are grown but at a small scale. There are valleys suitable for the production of rice. A nearby rural area, Zele, with intense farming was also considered for investigation.
  3. Abengourou, situated in the Eastern part of Côte d'Ivoire has a Guinean bioclimatic profile too. It has one rainy season with an average annual rainfall between 1200 mm and 1700 mm and an annual average temperature of 25.8°C. Coffee and cocoa are the main cash crops produced in the area. The breeding sites surveyed in Abengourou are all from urban areas.
  4. Yamoussoukro, located in the centre of the country is a transitional zone characterized by a wet savannah. This region has one rainy season and an average annual rainfall of 1200 mm and an annual average temperature of 25.8°C. Quantities of rice and vegetables are produced in the area for local consumption Breeding sites sampled there derived both from urban and sub-urban areas.
  5. Korhogo, in the North is a city located in a sub-Sudanian bioclimatic zone. The rainy season generally starts between May-June, peaks in August and ends in late October. The mean annual rainfall generally plateaus at 1200 mm with temperature varying over the seasons and the bioclimatic zone matching the savannah type. Cotton is the major commercial cash crop grown in this area. Some vegetables, rice and mangoes are also grown but at a small scale. Kaforo a nearby rural area where cotton is also intensively produced was included.

Mosquito collection

The sampling of An. gambiae larvae was done between June and July 2012 at the different sentinel sites (urban and rural). Larvae were collected from natural breeding sites such as ponds, puddles, footprints maintained by rainfall, and in rice and vegetable farms. They were then brought to the insectary at Institut Pierre Richet (IPR) and reared to adults.

Insecticide susceptibility tests

Bioassays on adult mosquitoes were conducted using WHO test kits [21]. Filter papers were impregnated with diagnostic concentrations of eight insecticides from 4 distinct classes as follows:

  1. Pyrethroids: permethrin (0.75%), deltamethrin (0.05%) and α-cypermethrin (0.05%);
  2. Pseudo-pyrethroid: etofenprox (0.05%);
  3. Organochlorides: Dichlorodiphenyltrichloroethane (DDT) (4%) and dieldrin (4%);
  4. Carbamate: carbosulfan (0.4%);
  5. Organophosphate: pyrimiphos-methyl (1%).

Bioassays were performed with batches of 25 unfed females of An. gambiae, 2–3 days old (four replicates per concentration). Mosquitoes were exposed to insecticide-impregnated papers for 60 min at 25±2°C and 80% relative humidity (RH). At the end of the exposure period, mosquitoes were transferred to observation tubes, supplied with 10% honey solution and held for 24 h before scoring delayed mortality. Batches of mosquitoes from each site exposed to untreated papers were used as negative controls.

In order to assess the involvement of detoxifying enzymes in the resistance phenotypes, complementary tests were performed with a 1 hour pre-exposure to piperonyl butoxide (PBO), an inhibitor of multiple function oxidases (MFO) and non-specific esterases (NSE). The wild mosquito populations were compared to the susceptible reference strain of An. gambiae s.s. Kisumu. All control survivor specimens (including the susceptible reference mosquito) were stored at −80°C for biochemical analysis. All samples exposed to the different insecticides were kept individually with silicagel at −20°C for molecular analysis.

Identification of sibling species and M and S molecular of An. gambiae s.s.

Ribosomal DNA was extracted from individual mosquitoes following Collins et al [22] and used for polymerase chain reaction (PCR) analysis to determine the species from the An. gambiae complex following Scott et al [23] and M or S molecular form of An. gambiae s.s. according to Faviaet al [24].

PCR detection of the L1014F and L1014S kdr and G119S ace-1 mutations

The presence of L1014F and L1014S kdr alleles was assessed using hot oligonucleotide ligation assay (HOLA) technique according to the protocol of Lynd et al [25]. The PCR-RFLP diagnostic test was used to detect the presence of G119S mutation (ace-1 gene) as described by Weill et al [26].

Biochemical analysis

Biochemical assays were performed to compare the amount of MFO and the activity levels of NSE for both ß and α-naphtyl acetate. The production of glutathione S-transferases (GST) was also investigated in all field samples relative to the susceptible Kisumu strain [27]. Mosquitoes used for the biochemical analysis have not been exposed to any insecticides prior to the assay.

Data analysis

WHO criteria [28] were adopted to define the resistance status of the mosquito populations. When less than 80% mortality was observed the population was considered ‘resistant’, between 80% and 97% mortality the population was considered ‘tolerant’ (or ‘suspected of resistance’) and when the mortality was above 98% the population was considered ‘susceptible’. Biochemical assay data (amounts of MFO and enzymatic activity per mg protein of NSE and GST) of Kisumu and the An. gambiae populations from the sentinel sites were compared using Kruskall-Wallis non-parametric test (Statistica software). Conformity of L1014F and L1014S kdr and ace-1 mutations frequencies with Hardy-Weinberg equilibrium was tested for An. gambiae population from the sentinel sites using the exact probability test [29]. The frequency of the duplicated allele ace-1D was estimated as described in Lenormand et al [30]. Statistical significance was set at the 5% level.

Results

Species and molecular forms of An. gambiae in Côte d'Ivoire

Table 1 illustrates the distribution of M and S An. gambiae s.s. molecular forms at the different sentinel sites. In 318 mosquito samples analysed for sibling species identification across the eco-geographical sentinel sites, only An. gambiae s.s. were found. The S molecular form was predominant in the sub-Sudanian area in the north both in urban and rural areas (71.9% in Korhogo and 65.6% in Kaforo) and in the rural area of Zele (78.1%) in the Western part of Côte d'Ivoire. The distribution in the neighbouring urban area of Man was different, the S form representing 43.8% only of the An. gambiae population. By contrast, in the pre-forested (Yamoussoukro) and the rain forest areas (San Pedro, Bingerville, Yopougon, Port-Bouët), the M form was predominant (65.6% in Bingerville and 100% in San-Pedro, Yopougon and Port-Bouët).

Resistance status

Figure 2 shows the insecticide resistance status of An. gambiae populations from the sentinel sites compared with the susceptible reference strain Kisumu. All insecticides tested killed 99–100% of susceptible mosquitoes indicating the accuracy of the active ingredient deposits on the filter papers used for the bioassays.

thumbnail
Figure 2. Insecticidal effects of diagnostic concentrations of insecticides against Anopheles gambiae mosquitoes from different sentinel sites (60 min contact in WHO tube tests).

https://doi.org/10.1371/journal.pone.0082387.g002

An. gambiae populations from all sites showed strong resistance to permethrin and DDT with some levels of tolerance to permethrin observed in the urban areas of Man, San-Pedro and Yopougon. Dieldrin resistance was reported in the few representative sites tested (Man, Yopougon, Yamoussoukro and Kaforo) (Fig. 2). Populations from all 10 sites were more vulnerable to deltamethrin than any other insecticides tested; mortality ranged from 74–100% between locations, with Korhogo and Kaforo (cotton areas) showing susceptibility (98–100% mortality) despite they displayed strong resistance to permethrin and DDT.

There was resistance to alpha-cypermethrin at all sites, except populations from Yopougon. A high frequency of resistance to etofenprox was recorded in all populations of An. gambiae s.s. (Fig. 2).

An. gambiae s.s. was strongly resistant to carbosulfan (14–52% mortality) at all site but showed reduced susceptibility to pyrimiphos-methyl in Yamoussoukro, Korhogo, Yopougon and Port-Bouët (66–87%) and full susceptibility in Abengourou and Man (99–100%).

Detection of resistance genes by PCR

The frequency of the kdr mutations was investigated in every sample from all sentinel sites. The distribution of the L1014F kdr mutation in the molecular forms of An. gambiae s.s. populations is shown in Table 2. The L1014F kdr was detected in both M and S molecular forms. Its average frequency ranged from 0.47 to 0.86 with the highest values (>0.70) reported in the coastal localities (Bingerville, Port-Bouët, San-Pedro) and the areas with cotton in the North (Korhogo, Kaforo). No L1014S kdr type was detected in all samples analysed (N = 317).

thumbnail
Table 2. Distributionof genotypes and allelic frequency of the L1014F kdr mutation in An.gambiaes.s.populations from sentinel sites.

https://doi.org/10.1371/journal.pone.0082387.t002

The G119S mutation conferring resistance to organophosphates and carbamates was detected in 4 localities out of 6 in both molecular forms (Table 3). The M and S specimens carrying the mutation were heterozygotes, except for 5 specimens out of 26 (19.2%) found homozygous in Korhogo in the North and 1 out of 12 (8.3%°) in Man in Western Côte d'Ivoire. Based on the Hardy-Weinberg equilibrium, this excess of heterozygous supports the presence of a duplicated allele (ace-1D) in An. gambiae s.s. from the sentinel sites. Table 3 shows the predicted distribution of the ace-1D in the S and M forms of An. gambiae s.s.. The presence of the ace-1D is suspected in either M or S forms of An. gambiae s.s. from all sites except Abengourou and San-Pedro (Table 3). It is only significantly supported (after Bonferroni correction) in M samples from Port-Bouët and S from Bingerville, the model suggesting a frequency around 0.47 for those from Port-Bouët (p<0.01) and 0.57 from Bingerville (p = 0.01).

thumbnail
Table 3. ace-1 allele frequency distribution in An. gambiae s.s. molecular forms from Côte d'Ivoire.

https://doi.org/10.1371/journal.pone.0082387.t003

Synergist and biochemical analysis

Synergist.

Figure 3 shows the toxicity of permethrin and carbosulfan with and without PBO against An. gambiae s.s. from the sentinel sites. Pre-exposure of mosquitoes to PBO significantly increased the mortality rates to permethrin from all locations (p<0.05) except Yopougon (p = 0.539) (Figure 3a). Interestingly, PBO increased mortality to permethrin to full susceptibility level in populations from Korhogo (from 39% to 100%) and Abengourou (62% to 97%), suggesting that the resistance phenotype in these two areas is almost entirely mediated by the metabolic activities of MFOs and NSE.

thumbnail
Figure 3. Insecticidal effects of diagnostic concentrations of one pyrethroid insecticide (permethrin) and one carbamate insecticide (carbosulfan) (60 min contact in WHO tube tests) with or without a 60 min pre-exposition to PBO.

https://doi.org/10.1371/journal.pone.0082387.g003

Resistance to carbosulfan was largely pronounced across sites, mortality ranging from 14% to 63%. Pre-exposure of mosquitoes to PBO significantly increased mortality rates to carbosulfan at all sites but resistance was not fully synergised by this MFO and NSE inhibitor, inferring a residual role of additional mechanisms, including the ace-1 site insensitivity.

Biochemical assays.

Table 4 compares the mean amount of MFO and the mean activities of the NSE and GST found in An. gambiae s.s. from the sentinel sites relative to the levels in the susceptible An. gambiae Kisumu.

thumbnail
Table 4. Mean level of NSE. MFO and GST activity in An. gambiae s.s. populations from the sentinel sites relative to the susceptible reference strain Kisumu.

https://doi.org/10.1371/journal.pone.0082387.t004

All wild samples analyzed displayed significantly higher levels of esterase activity (using α-naphtyl acetate as a substrate) than that measured for Kisumu (p<0.05), except samples from Abengourou and Zele. The trend in activity for assays using ß-naphtyl acetate as a substrate was similar to the α-naphtyl assays. The highest α-and ß-esterase activities (>2.5-fold) were recorded on the coast (Yopougon, Port-Bouët) and the cotton growing area of Korhogo.

The production of MFO (cytochrome P450s) amount in An. gambiae from all sites but not Zele was significantly higher than the content measured in Kisumu (p<0.05), with broadly similar mean levels of activity expressed between locations.

The GST levels of activity in Yopougon, Port-Bouët, Korhogo, Yamoussoukro and Man populations were significantly higher than in Kisumu and the other field populations (p<0.05).

Overall a patchy distribution of overproduced quantities of NSE, MFO and GST was found in most An. gambiae s.s. populations analysed relative to the normal strain Kisumu; the most affected areas being on the coastal urban part of the country (Yopougon and Port-Bouët with vegetables and horticulture) and Korhogo in the North where cotton is produced.

Discussion

A nation-wide survey of insecticide resistance in An. gambiae s.s. in Côte d'Ivoire was conducted in order to have an overview of the resistance status and establish a baseline dataset that would guide the National Malaria Control Programme.

The bioassays data showed high resistance levels of An. gambiae s.s. to organochlorides, carbamates and pyrethroids but at a lesser extent toward deltamethrin. The resistance levels to the organophosphate (pyrimiphos-methyl) varied greatly from susceptibility to resistance across sites. The kdr and ace-1R mutations were highly expressed in all An. gambiae s.s. populations with the exception for the ace-1R in samples from Abengourou and San-Pedro. The ace-1D duplication was present in Bingerville and Port-Bouët respectively in S and M molecular forms confirming previous finding highlighting its presence in Côte d'Ivoire [31], [32]. Its spread seemed to be wider as it was suspected in the urban areas of Yopougon, Yamoussoukro and Man although the low sample size did not allow to confirm it.

The MFO quantities and, NSE and GST activities in almost all An. gambiae populations were significantly higher than in the susceptible strain Kisumu and the bioassays data with PBO indicated that at least the MFO and NSE were involved in the phenotypic expression of the resistance.

This study provides a wider view of the spread of insecticide resistance in Côte d'Ivoire and adds a significant baseline knowledge to recent reports of high insecticide resistance levels in M'Bé [19], Yaokoffikro [20] near Bouaké in central Côte d'Ivoire, Tiassalé [18], [33] and Adzopé [17] in southern Côte d'Ivoire. All these reports including the present highlight the level of spread of resistance to all class of insecticides deployed for vector control, either as LLIN already in place or IRS under consideration for implementation in Côte d'Ivoire; an increasing trend commonly shared by several countries in West Africa[34][39].

The association between agricultural practices and the build up of insecticide resistance has been intensively investigated. It is worth noticing that An. gambiae s.s. population from areas with massive cotton production (Korhogo, Kaforo) displayed the highest resistance levels against almost all insecticides classes (organochlorides, pyrethroids and carbamates) as previously reported [34], [40][44]. Urban vegetable farming in Yopougon, Port-Bouët and Yamoussoukro areas was also associated with high levels of resistance, as previously detected [45][48]. The trend was not so clear within the rice-growing areas, generally associated with low application of insecticide [40], [44] but where moderate to high resistance level was found in An.gambiae [49].

There is a growing expansion of rubber production in Côte d'Ivoire and generally across part of Africa. So far reports on insecticide resistance in such context are scarce [50]. In-depth studies documenting the use of insecticides both socially and in the o agricultural sector is highly stressed in order to better appraise factors selecting and driving the evolution of insecticide resistance mechanisms [32].

The significance of the present study for the NMCP in Côte d'Ivoire raises an important question of whether to continue to deploy pyrethroid based LLINs and IRS towards which resistance continues to rise with no guarantee that the level of resistance seen in the country would not compromise their efficacy. No doubt it would be difficult to demonstrate the impact of resistance on the effectiveness of any of these interventions. There have been extensive randomized controlled trials (RCTs) (phase III) in part of Africa aiming at investigating the efficacy of ITNs for malaria prevention [14], [51] but very few have assessed how pyrethroid resistance might affect the effectiveness of such intervention. RCTs entail a set of communities randomly divided into groups, one that receives the novel form of vector control intervention, and comparison arms that often receive the old form of vector control tools or nothing. The key difficulty is that it is impossible to address the question to whether vector control would produce a smaller reduction in malaria if the vector mosquitoes are resistant than it would have done if they were susceptible, using RCT methods. This is simply because resistance is not an easy factor that can be allocated randomly to some communities and not to others. The distribution of resistance is patchy and its severity seems to differ from one area (locality) to another as seen in the present study. Moreover there may be more resistance or survival trend of mosquitoes in some villages than others because of variations in the quality of vector control operations [52] or in mosquito behavior [53]. This means there is no straight forward solution to address resistance impact. Some researchers are currently striving to bring up to the market novel insecticides or combination of these with pyrethroids to circumvent insecticide resistance but this may take several years before they are approved and made available for use. In the short term, NMCPs could consider adopting intervention strategies that e.g. combine a synergist such as PBO with pyrethroids in LLINs. The technology is intended to alleviate the load of resistance by removing the metabolic component of resistance due to MFOs and NSE [54]. The well known commercialized prototypes PermaNet® 3.0 [55][58] and more recently Olyset Plus® [59] of this kind have proved to be highly active against pyrethroid resistant mosquitoes and could be the product of choice for NMCPs. Ultimately, where LLINs are already in place and vectors survive them because of pyrethroid resistance [12], [13], NMCPs could consider spraying the homes with a non pyrethroid insecticide. In any event, insecticide combinations within homes for malaria control is an unavoidable reality. LLIN coverage is going universal and IRS with non-pyrethroid insecticides is being applied concurrently with LLINs as malaria control policy in many areas of high malaria transmission. Given the high resistance to carbamates detected across all sites but an appreciable level of susceptibility to the organophosphate, pyrimiphos methyl, it would be advisable to deploy this insecticide for IRS in combination with LLINs in Côte d'Ivoire. Recent evaluation of this insecticide applied as IRS in experimental huts showed extremely high level of control of pyrethroid resistant mosquitoes [60].

Conclusion

The insecticide resistance data presented in the sentinel survey are of great interest for the NMCP of Côte d'Ivoire and beyond, to countries also facing similar rise in the level of pyrethroid resistance in their local vectors. Innovative strategies that combine insecticide and synergists in LLINs or spatially LLIN and effective insecticide for IRS within homes could be in the short term the best practice for NMCPs to manage insecticide resistance in malaria vectors in endemic countries.

Acknowledgments

Sincere thanks go to Dr Koffi San, former Head of the National Malaria Control Programme of Côte d'Ivoire and the research unit of NMCP (Dr Serge B Assi., Dr Florence-Judith Kadjo). We are very grateful to all the staff at the Institut Pierre Richet, Bouaké, Côte d'Ivoire, particularly Aboubacar Koné, for his hard work during the field and laboratory experiments. We also thank to Aboubacar Sidick from the Institut de Recherche pour le Développement (IRD)/Centre de Recherche Entomologique de Cotonou (CREC), Benin, for his technical assistance in biochemical and molecular experiments. We are very grateful to Dr Pierrick Labbé for his crucial help on the ace-1D frequency calculations. Special thanks also go to Bamoro Coulibaly for mapping the study site.

Author Contributions

Conceived and designed the experiments: AAK LPAA. Performed the experiments: LPAA JPKK. Analyzed the data: AAK LPAA CP. Contributed reagents/materials/analysis tools: AAK CP. Wrote the paper: AAK LPAA RN CP.

References

  1. 1. WHO (2011) World Malaria Report. Geneva. 259 p.
  2. 2. WHO (2012) WHO recommended long-lastin insecticidal mosquito nets. Geneva. Available: http://www.who.int/whopes/Long_lasting_insecticidal_nets_Jul_2012.pdf p.
  3. 3. WHO (2009) WHO recommended insecticides for indoor residual spraying against malaria vectors. Geneva. Available: http://www.who.int/whopes/Insecticides_IRS_Malaria_09.pdf p.
  4. 4. Ranson H, N'Guessan R, Lines J, Moiroux N, Nkuni Z, et al.. (2010) Pyrethroid resistance in African anopheline mosquitoes: what are the implications for malaria control? Trends Parasitol.
  5. 5. Lapied B, Pennetier C, Apaire-Marchais V, Licznar P, Corbel V (2009) Innovative applications for insect viruses: towards insecticide sensitization. Trends Biotechnol 27: 190–198.
  6. 6. Rivero A, Vezilier J, Weill M, Read AF, Gandon S (2010) Insecticide control of vector-borne diseases: when is insecticide resistance a problem? PLoS Pathog 6..
  7. 7. Hemingway J, Ranson H (2000) Insecticide resistance in insect vectors of human disease. Annu Rev Entomol 45: 371–391.
  8. 8. Raymond-Delpech V, Matsuda K, Sattelle BM, Rauh JJ, Sattelle DB (2005) Ion channels: molecular targets of neuroactive insecticides. Invert Neurosci 5: 119–133.
  9. 9. Hemingway J, Hawkes NJ, McCarroll L, Ranson H (2004) The molecular basis of insecticide resistance in mosquitoes. Insect Biochem Mol Biol 34: 653–665.
  10. 10. Corbel V, N'Guessan R, Brengues C, Chandre F, Djogbenou L, et al. (2007) Multiple insecticide resistance mechanisms in Anopheles gambiae and Culex quinquefasciatus from Benin, West Africa. Acta Trop 101: 207–216.
  11. 11. VectorBase (2012) AnoBase: the Anopheles database. Insecticide resistance database search. Available: http://anobase.vectorbase.org/ir/. Accessed 2012 Oct 23.
  12. 12. N'Guessan R, Corbel V, Akogbeto M, Rowland M (2007) Reduced efficacy of insecticide-treated nets and indoor residual spraying for malaria control in pyrethroid resistance area, Benin. Emerging Infectious Diseases 13: 199–206.
  13. 13. Asidi A, N'Guessan R, Akogbeto M, Curtis C, Rowland M (2012) Loss of Household Protection from Use of Insecticide-Treated Nets against Pyrethroid-Resistant Mosquitoes, Benin. Emerg Infect Dis 18: 1101–1106.
  14. 14. Corbel V, Akogbeto M, Damien GB, Djenontin A, Chandre F, et al.. (2012) Combination of malaria vector control interventions in pyrethroid resistance area in Benin: a cluster randomised controlled trial. Lancet Infect Dis.
  15. 15. PNLP_Cote_d'Ivoire (2012) Plan stratégique 2012–2015 du programme national de lutte contre le paludisme de la Côte d'Ivoire. Abidjan. 110 p.
  16. 16. Edi CV, Koudou BG, Jones CM, Weetman D, Ranson H (2012) Multiple-insecticide resistance in Anopheles gambiae mosquitoes, Southern Cote d'Ivoire. Emerg Infect Dis 18: 1508–1511.
  17. 17. Ahoua Alou LP, Koffi AA, Adja MA, Assi SB, Kouassi PK, et al. (2012) Status of pyrethroid resistance in Anopheles gambiae s. s. M form prior to the scaling up of Long Lasting Insecticidal Nets (LLINs) in Adzope, Eastern Cote d'Ivoire. Parasit Vectors 5: 289.
  18. 18. Ahoua Alou LP, Koffi AA, Adja MA, Tia E, Kouassi PK, et al. (2010) Distribution of ace-1R and resistance to carbamates and organophosphates in Anopheles gambiae s.s. populations from Cote d'Ivoire. Malar J 9: 167.
  19. 19. Koffi AA, Ahoua Alou LP, Adja AM, Chandre F, Pennetier C (2013) Insecticide Resistance status of Anopheles gambiae s.s. population from M'Bé, a WHOPES labelled experimental hut station, 10 years after the political crisis in Côte d'Ivoire. Malaria Journal 12: 151.
  20. 20. Koffi AA, Alou LP, Adja MA, Kone M, Chandre F, et al. (2012) 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 Cote d'Ivoire. Parasit Vectors 5: 68.
  21. 21. WHO editor (2006) Guidelines for testing mosquito adulticides intended for Indoor Residual Spraying (IRS) and Insecticide Treated Nets (ITNs). WHO/CDS/NTD/WHOPES/GCDDP/2006.3.
  22. 22. Collins FH, Finnerty V, Petrarca V (1988) Ribosomal DNA-probes differentiate five cryptic species in the Anopheles gambiae complex. Parassitologia 30: 231–240.
  23. 23. Scott JA, Brogdon WG, Collins FH (1993) Identification of single specimens of the Anopheles gambiae complex by the polymerase chain reaction. Am J Trop Med Hyg 49: 520–529.
  24. 24. Favia G, Lanfrancotti A, Spanos L, Siden Kiamos I, Louis C (2001) Molecular characterization of ribosomal DNA polymorphisms discriminating among chromosomal forms of Anopheles gambiae s.s. Insect Molecular Biology 10: 19–23.
  25. 25. Lynd A, Ranson H, McCall PJ, Randle NP, Black WCt, et al. (2005) A simplified high-throughput method for pyrethroid knock-down resistance (kdr) detection in Anopheles gambiae. Malar J 4: 16.
  26. 26. Weill M, Malcolm C, Chandre F, Mogensen K, Berthomieu A, et al. (2004) The unique mutation in ace-1 giving high insecticide resistance is easily detectable in mosquito vectors. Insect Mol Biol 13: 1–7.
  27. 27. WHO (1998) Techniques to detect resistance mechanisms (Field and laboratory manual). Geneva. 42 p.
  28. 28. WHO (1998) Tests procedures for insecticide resistance monitoring in malaria vectors, bioefficacy and persistence of insecticides on treated surfaces. WHO/CDS/CPC/MAL/1998.12.
  29. 29. Raymond M, Rousset F (1995) GENEPOP Version 1.2 A population genetics software for exact tests and ecumenicism. Journal of Heredity: 248–249.
  30. 30. Lenormand T, Guillemaud T, Bourguet D, Raymond M (1998) Appearance and sweep of the gene duplication: adaptative response and potential for new functions in the mosquito Culex pipiens. Evolution 52: 1705–1712.
  31. 31. Djogbenou L, Labbe P, Chandre F, Pasteur N, Weill M (2009) Ace-1 duplication in Anopheles gambiae: a challenge for malaria control. Malar J 8: 70.
  32. 32. Koffi AA, Chandre F, Ahoua Alou LP, Labbe P, Diabate A, et al.. Agricultural practices and the distribution of ace-1 duplicated alleles in Anopheles gambiae s.s. from Côte-d'Ivoire, West Africa. BMC Infect Dis submitted.
  33. 33. Edi CV, Koudou BG, Jones CM, Weetman D, Ranson H (2012) Multiple-insecticide resistance in Anopheles gambiae mosquitoes, Southern Cote d'Ivoire. Emerg Infect Dis 18: 1508–1511.
  34. 34. Djogbenou L, Pasteur N, Akogbeto M, Weill M, Chandre F (2010) Insecticide resistance in the Anopheles gambiae complex in Benin: a nationwide survey. Med Vet Entomol.
  35. 35. Badolo A, Traore A, Jones CM, Sanou A, Flood L, et al. (2012) Three years of insecticide resistance monitoring in Anopheles gambiae in Burkina Faso: resistance on the rise? Malar J 11: 232.
  36. 36. Dabire RK, Namountougou M, Sawadogo SP, Yaro LB, Toe HK, et al. (2012) Population dynamics of Anopheles gambiae s.l. in Bobo-Dioulasso city: bionomics, infection rate and susceptibility to insecticides. Parasit Vectors 5: 127.
  37. 37. Namountougou M, Simard F, Baldet T, Diabate A, Ouedraogo JB, et al. (2012) Multiple Insecticide Resistance in Anopheles gambiae s.l. Populations from Burkina Faso, West Africa. PLoS ONE 7: e48412.
  38. 38. Namountougou M, Diabate A, Etang J, Bass C, Sawadogo SP, et al.. (2012) First report of the L1014S kdr mutation in wild populations of Anopheles gambiae M and S molecular forms in Burkina Faso (West Africa). Acta Trop.
  39. 39. Djegbe I, Boussari O, Sidick A, Martin T, Ranson H, et al. (2011) Dynamics of insecticide resistance in malaria vectors in Benin: first evidence of the presence of L1014S kdr mutation in Anopheles gambiae from West Africa. Malar J 10: 261.
  40. 40. Fane M, Cisse O, Traore CS, Sabatier P (2012) Anopheles gambiae resistance to pyrethroid-treated nets in cotton versus rice areas in Mali. Acta Trop 122: 1–6.
  41. 41. Yadouleton A, Martin T, Padonou G, Chandre F, Asidi A, et al. (2011) Cotton pest management practices and the selection of pyrethroid resistance in Anopheles gambiae population in northern Benin. Parasit Vectors 4: 60.
  42. 42. Chouaibou M, Etang J, Brevault T, Nwane P, Hinzoumbe CK, et al. (2008) Dynamics of insecticide resistance in the malaria vector Anopheles gambiae s.l. from an area of extensive cotton cultivation in Northern Cameroon. Trop Med Int Health 13: 476–486.
  43. 43. Muller P, Chouaibou M, Pignatelli P, Etang J, Walker ED, et al. (2008) Pyrethroid tolerance is associated with elevated expression of antioxidants and agricultural practice in Anopheles arabiensis sampled from an area of cotton fields in Northern Cameroon. Mol Ecol 17: 1145–1155.
  44. 44. Diabate A (2002) The Role of Agricultural Uses of Insecticides in Resistance to Prethroids in Anopheles gambiae S.L. in Burkina Faso. Am J Trop Med Hyg 67: 617–622.
  45. 45. Yadouleton A, N'Guessan R, Allagbe H, Asidi A, Boko M, et al. (2010) The impact of the expansion of urban vegetable farming on malaria transmission in major cities of Benin. Parasit Vectors 3: 118.
  46. 46. Tia E, Akogbeto M, Koffi A, Toure M, Adja AM, et al. (2006) [Pyrethroid and DDT resistance of Anopheles gambiae s.s. (Diptera: Culicidae) in five agricultural ecosystems from Cote-d'Ivoire]. Bull Soc Pathol Exot 99: 278–282.
  47. 47. Antonio-Nkondjio C, Fossog BT, Ndo C, Djantio BM, Togouet SZ, et al. (2011) Anopheles gambiae distribution and insecticide resistance in the cities of Douala and Yaounde (Cameroon): influence of urban agriculture and pollution. Malar J 10: 154.
  48. 48. Jones CM, Toe HK, Sanou A, Namountougou M, Hughes A, et al. (2012) Additional selection for insecticide resistance in urban malaria vectors: DDT resistance in Anopheles arabiensis from Bobo-Dioulasso, Burkina Faso. PLoS ONE 7: e45995.
  49. 49. Bigoga JD, Ndangoh DN, Awono-Ambene PH, Patchoke S, Fondjo E, et al. (2012) Pyrethroid resistance in Anopheles gambiae from the rubber cultivated area of Niete, South Region of Cameroon. Acta Trop 124: 210–214.
  50. 50. Bigoga JD, Nanfack FM, Awono-Ambene PH, Patchoke S, Atangana J, et al. (2012) Seasonal prevalence of malaria vectors and entomological inoculation rates in the rubber cultivated area of Niete, South Region of Cameroon. Parasit Vectors 5: 197.
  51. 51. Henry MC, Assi SB, Rogier C, Dossou-Yovo J, Chandre F, et al. (2005) Protective efficacy of lambda-cyhalothrin treated nets in Anopheles gambiae pyrethroid resistance areas of Cote d'Ivoire. Am J Trop Med Hyg 73: 859–864.
  52. 52. malERA , Aultman K, Burkot TR, Chandre F, Coetzee M, et al. (2011) A research agenda for malaria eradication: vector control. PLoS Med 8: e1000401.
  53. 53. Moiroux N, Gomez MB, Pennetier C, Elanga E, Djenontin A, et al. (2012) Changes in Anopheles funestus biting behavior following universal coverage of long-lasting insecticidal nets in Benin. J Infect Dis 206: 1622–1629.
  54. 54. Hemingway J (2000) The molecular basis of two contrasting metabolic mechanisms of insecticide resistance. Insect Biochem Mol Biol 30: 1009–1015.
  55. 55. Koudou BG, Koffi AA, Malone D, Hemingway J (2011) Efficacy of PermaNet(R) 2.0 and PermaNet(R) 3.0 against insecticide-resistant Anopheles gambiae in experimental huts in Cote d'Ivoire. Malar J 10: 172.
  56. 56. Corbel V, Chabi J, Dabire RD, Etang J, Nwane P, et al. (2010) Field efficacy of a new mosaic long-lasting mosquito net (PermaNet 3.0) against pyrethroid-resistant malaria vectors: a multi centre study in Western and Central Africa. Malar J 9: 113.
  57. 57. N'Guessan R, Asidi A, Boko P, Odjo A, Akogbeto M, et al. (2010) An experimental hut evaluation of PermaNet((R)) 3.0, a deltamethrin-piperonyl butoxide combination net, against pyrethroid-resistant Anopheles gambiae and Culex quinquefasciatus mosquitoes in southern Benin. Trans R Soc Trop Med Hyg 104: 758–765.
  58. 58. Van Bortel W, Chinh VD, Berkvens D, Speybroeck N, Trung HD, et al. (2009) Impact of insecticide-treated nets on wild pyrethroid resistant Anopheles epiroticus population from southern Vietnam tested in experimental huts. Malar J 8: 248.
  59. 59. Pennetier C, Bouraima A, Chandre F, Piameu M, Etang J, et al.. (2013) Efficacy of Olyset® Plus, a new long-lasting insecticidal net incorporating permethrin and piperonil-butoxide against multi-resistant malaria vectors. PLoS ONE submitted.
  60. 60. Rowland M, Boko P, Odjo A, Asidi A, Akogbeto M, et al. (2013) A new long-lasting indoor residual spray formulation of the organophosphate pirimiphos methyl for prolonged malaria vector control: an experimental hut trial against pyrethroid-resistant Anopheles gambiae and Culex quinquefasciatus in Southern Benin. PLoS ONE 8: e69516.