- Open Access
Dynamics of multiple insecticide resistance in the malaria vector Anopheles gambiae in a rice growing area in South-Western Burkina Faso
Malaria Journalvolume 7, Article number: 188 (2008)
Insecticide resistance of the main malaria vector, Anopheles gambiae, has been reported in south-western Burkina Faso, West Africa. Cross-resistance to DDT and pyrethroids was conferred by alterations at site of action in the sodium channel, the Leu-Phe kdr mutation; resistance to organophosphates and carbamates resulted from a single point mutation in the oxyanion hole of the acetylcholinesterase enzyme designed as ace-1R.
An entomological survey was carried out during the rainy season of 2005 at Vallée du Kou, a rice growing area in south-western Burkina Faso. At the Vallée du Kou, both insecticide resistance mechanisms have been previously described in the M and S molecular forms of An. gambiae. This survey aimed i) to update the temporal dynamics and the circumsporozoite infection rate of the two molecular forms M and S of An. gambiae ii) to update the frequency of the Leu-Phe kdr mutation within these forms and finally iii) to investigate the occurrence of the ace-1Rmutation.
Mosquitoes collected by indoor residual collection and by human landing catches were counted and morphologically identified. Species and molecular forms of An. gambiae, ace-1Rand Leu-Phe kdr mutations were determined using PCR techniques. The presence of the circumsporozoite protein of Plasmodium falciparum was determined using ELISA.
Anopheles gambiae populations were dominated by the M form. However the S form occurred in relative important proportion towards the end of the rainy season with a maximum peak in October at 51%. Sporozoite rates were similar in both forms. The frequency of the Leu-Phe kdr mutation in the S form reached a fixation level while it is still spreading in the M form. Furthermore, the ace-1Rmutation prevailed predominately in the S form and has just started spreading in the M form. The two mutations occurred concomitantly both in M and S populations.
These results showed that the Vallée du Kou, a rice growing area formerly occupied mainly by M susceptible populations, is progressively colonized by S resistant populations living in sympatry with the former. As a result, the distribution pattern of insecticide resistance mutations shows the occurrence of both resistance mechanisms concomitantly in the same populations. The impact of multiple resistance mechanisms in M and S populations of An. gambiae on vector control measures against malaria transmission, such as insecticide-treated nets (ITNs) and indoor residual spraying (IRS), in this area is discussed.
Malaria transmission in tropical Africa is mainly dominated by Anopheles gambiae complex, including An. gambiae s.s., the most anthropophilic vector transmitting malaria in sub-Saharan Africa . Formerly considered as a single species, it began early to accumulate genetic heterogeneity. In the 1980s, cytogenetic studies based on chromosomal inversion arrangements found five incipient chromosomal forms. The suspected role of these inversions was to restrict gene flow among the forms and to provide adaptation to different ecological settings [2, 3]. In Burkina Faso, Mopti and Savanna chromosomal forms dominated An. gambiae population structure . These chromosomal forms appeared more or less genetically isolated in the field, presumably through prezygotic barriers since viable and fertile hybrids have been obtained in the laboratory [5–7]. However, cytogenetic analysis is not a precise way to evaluate the degree of hybridization between forms because of the presence of cryptic 'heterokaryotypes'. Recent studies based on molecular markers such as X-linked ribosomal DNA suggested the existence of only two entities within An. gambiae, referred to as M and S molecular forms . So far, in Burkina Faso and Mali, in savannah environments of West Africa, Mopti and Savanna chromosomal forms coincide respectively with M and S molecular forms . As a main malaria major vector with high level of polymorphism, An. gambiae has been a subject of many investigations in West Africa, such as bio-ecology and insecticide resistance studies [9–12].
In Burkina Faso, a study carried out in the mid-1980s in Vallée du Kou , showed that the Mopti chromosomal form was dominating. Molecular-based identification of An. gambiae s.s. populations conducted in this area in 1999 and 2000 confirmed that the M molecular form predominated throughout the year [11, 12, 14], although with some temporal variations.
The S molecular form occurred in low frequency until the end of the rainy season (October/November), when it peaked around 30% , as found previously by Robert et al . A similar pattern was found in the same environment in Mali [7, 9]. More recently, in 2004 the S form was observed in Vallée du Kou towards the end of the rainy season at a frequency of 50% . This study gave no results regarding the frequencies of the two forms throughout the malaria transmission season.
In West Africa, the main mechanism involved in pyrethroid-resistance in An. gambiae is caused by target site insensitivity through a knockdown resistance (kdr)-like mutation caused by a single point mutation (Leu-Phe) in the para-sodium channel gene . Preliminary surveys done in Vallée du Kou in 1999 indicated that the Leu-Phe kdr mutation has been found almost only in the S form at high allelic frequency (0.95) compared to just 0.006 in the M form . However, the spread of the mutation in the M population seemed an ongoing process in Vallée du Kou as it increased to a frequency of 0.02 in 2000. Nowadays, the population structure of An.gambiae and their pyrethroid resistance status are probably modified with the changing in agricultural practices needing intensive use of insecticides (cotton and vegetable cropping) and also the increasing of man made breeding sites as puddles throughout the village.
Furthermore, it has been recently noted in An. gambiae from the same area the occurrence of a single point mutation (glycine to serine at position 119) in the oxyanion hole of the acetylcholinesterase enzyme . This mutation named ace-1Rmutation was associated with insensitivity of An. gambiae to organophosphates and carbamates .
The objective of the present study was to gather recent information on the dynamics of the two molecular forms of An. gambiae throughout the malaria transmission season in this area with particular attention to resistance mechanisms. This information is crucial for a proper evaluation of new insecticides or vector control tools expected to be involved in malaria control and resistance management.
Materials and methods
Vallée du Kou (4°25' W, 11°24' N) is a rice-growing valley covering 1,200 ha and comprising seven villages, with a total of 4,470 habitants surrounded by humid savannah. The rainy season extends from June to October and the dry season from November to May. The Kou River is a permanent source of irrigation water and there are two rice crops per year from July to November and from January to May. Few insecticides are used on rice, but huge amounts of insecticides are used extensively in cotton fields' located on the outside of rice fields. During the last two years, some producers started to grow vegetables involving intensive use of insecticides. Mosquito collection had been carried out in the village numbered as seven (VK7), which is located at the end of the rice fields (Figure 1). VK7 has about 600 inhabitants, mainly farmers. Sheep, goats, pigs, and a few cows are also present. Cotton and maize fields surrounded this village.
The irrigation system and rice fields provide year-round mosquito breeding. Additional breeding sites are created by rains in the depressions and ponds. High densities of An. gambiae (up to 200 bites/person/night) are recorded during the rainy season . The two molecular forms M and S of An. gambiae occur in sympatry notably at the end of the rainy season [12, 15]. The kdr-based mechanism conferring resistance to pyrethroids and DDT and the ace-1Rmutation conferring resistance to organophosphates and carbamates predominate in the S form [11, 12, 17].
Anopheline mosquitoes were sampled during the rainy season from July to November 2005 mainly by indoor manual collection for resting mosquitoes and, secondarily by human landing catches.
Indoor resting mosquitoes were collected regularly in four houses in the village early in the morning by manual aspirators four times a month, on four consecutive days. These indoor collections carried out monthly from July to November were used to establish the temporal dynamics of the two molecular forms of An. gambiae and their resistance status by PCR.
Human landing catches were ensured only in August and in October to evaluate the sporozoite infection rate of each molecular form. These two periods correspond to the peak of the M form in August and of the S form in October. Catches were carried out between 20.00 and 06.00, during two consecutive nights for each period inside and just outside of four houses of the village. They were performed by informed volunteers who were provided free and rapid treatment when suspected clinical signs of malaria according to World Health Organization (WHO)-recommended regimen on the basis of fever and detectable Plasmodium falciparum parasitemia.
Laboratory processing of mosquitoes
Anophelines were sorted and assigned to species based on morphological characters using standard identification keys . Females referring to An. gambiae complex were processed by PCR concomitantly for identification of species and molecular forms of An. gambiae [20, 8]. Detection of Leu-Phe kdr and ace-1Rmutations were performed on indoor resting specimens by PCR from genomic DNA following Martinez-Torres et al  and Weill et al , respectively.
The head-thoraces of anopheline females issued from human landing catches were tested for the presence of the circumsporozoite protein (CSP) of P. falciparum, the major malarial parasite occurring in the study area, by enzyme-linked immunosorbent analysis (ELISA) according to Wirtz et al . Samples of August were used for ELISA process because this month corresponded to the peak of the M form whereas no or few individuals of S form were found in indoor collections.
The sporozoite rate was defined as the proportion of mosquitoes found positive for P. falciparum CS protein. The entomological inoculation (EIR) was calculated as the product of HBR and the sporozoite rate of mosquitoes caught on landing collections.
Dynamics of the M and S molecular forms
Overall, 330 mosquitoes were analysed for identification of species and molecular forms of An. gambiae (Figure 2). Anopheles arabiensis was absent. Both M and S molecular forms occurred. The overall frequency of the molecular M form was higher than that of the S form reaching respectively 69% versus 31%, but the relative prevalence of these forms throughout the collecting period showed some monthly varying frequencies. During the first part of the rainy season (July and August), the M form predominated, whereas the two forms are found in similar proportion during the second part of the rainy season (from September to November). Indeed in July all mosquitoes analysed (n = 32) were only of the M form. The S form started appearing toward the end of August at low frequency of 4% (3/84) reaching a maximum peak of 51% (38/75) in October equalling the frequency of the M form.
The sporozoite rate was determined in An. gambiae females issued from indoor human landing catches carried out in two sampling sets: August and October 2005 (Table 1). Because in August no specimens of the S form were found, sporozoite rate was determined only by the M form reaching 1.19%.
In October, the M and S forms were collected in equal numbers. Results yielded an average sporozoite rate of 3.08% with no significant difference between the two forms (Fisher's Exact test, P = 0.33).
Distribution of the kdr mutation
An average of 49 mosquitoes issued from indoor resting fauna were analysed monthly for the Leu-Phe kdr mutation (Table 2). The kdr mutation occurred in both M and S forms with varying allelic frequencies in the M form. Kdr frequency was higher in the S form irrespective of the month (χ2 = 6.14, df = 2, P < 0.02) reaching 0.93 in mean and ranged from 0.50 to 0.98 throughout its occurrence period. None homozygous sensible individual was identified from the S form. The mean frequency of the kdr mutation in the form M was 0.096 and ranged from 0 to 0.23 throughout the five months. Only 8 mosquitoes from the M form were identified in September and November as homozygous RR for this mutation vs. 15 hybrids RS overall among 161 specimens analysed. In the S form the homozygous resistant predominated irrespective of the month contributing thus to achieve the fixation process of this allele in the S population.
Distribution of the ace-1Rmutation
58 mosquitoes issued from indoor resting fauna were analysed for the ace-1Rmutation from September and October 2005, corresponding to 29 specimens per month (Table 3). The ace-1Rmutation was detected at low frequency in the M form with only one heterozygous RS per month corresponding respectively to a mean allelic frequency of 0.031. Conversely in the S form, this mutation occurred in a relative higher frequency comparing to the M form (χ2 = 6.75, df = 1, P < 0.01) with an equal mean of 0.37 per month.
One individual from the S form was detected as a homozygous resistant for this mutation in September. In the S form, all specimens heterozygous for the ace-1Rmutation (17 females) had also the kdr mutation in heterozygous status. The one individual homozygous ace-1R/ace-1Rwas also homozygous kdrr/kdrr. Inversely in the M form, the two individual detected as heterozygous ace-1r/ace-1swere homozygous kdrs/kdrs.
In Burkina Faso, current chromosomal and molecular forms of An. gambiae s.s. are tightly correlated [4, 10]. Formerly Robert et al.  studying the distribution of An. gambiae s.s. cytotypes in the rice field area of Vallée du Kou in 1984 showed a predominance of the Mopti chromosomal form. With the progress in molecular genetic, this distribution has been updated in 1999 and 2000 by Diabate et al [11, 12], pointing out the predominance of the M molecular form corresponding to the Mopti chromosomal form. However the occurrence of the S molecular form (corresponding to the Savanna chromosomal form) has been observed toward the end of the rainy season. The results obtained in 2005 followed the same pattern of distribution, but the overall proportion of the S form has increased further reaching a maximum of 51% vs. 24% at the same period in 2000 . Taking to account that both studies were performed in the same place (VK7) and during the same period (from July to November), it appears that the relative frequency of the S form has increased. Future studies to confirm this trend are encouraged.
The increase in the S form could probably be due to human activities creating new habitats for the S population in this area, such as house building using bricks of banco and muddles taken out from the soil. Such activities increase the number of temporary rain-filled breeding sites like pits, ponds and puddles that are favourable to S form development throughout the village.
In general, in West Africa, the S form is not well adapted to rice paddies, whilst the M form develops rather well [24, 25]. However, some of the rice paddies used to grow vegetables and some irrigation canals not well managed could constitute during the rainy season, suitable habitats for the S larvae.
The changing in the vector population structure may presumably increase the malaria transmission level. Even though the malaria vector population was dominated by the M form throughout the year, this form had a low sporozoite rate because of the low parity rate observed within its population . During the sympatry period (October), the sporozoite index in the S form was not statistically different to that of the M form, but these are preliminary results that would need to be confirmed on a bigger sample and during different periods of the year.
The main mechanism conferring resistance of An. gambiae to pyrethroids in West Africa, the Leu-Phe kdr mutation, did not vary in the S form (93%) compared to its frequency in 1999 and 2000 , as expected.
The spread of the kdr mutation is an ongoing process in the M form as its allelic frequency in 2005 was five fold higher than in 2000. Indeed six years ago, the kdr mutation was found occurring only in the S form and some investigations conducted in VK7 during the rainy season 1999 failed to identify this mutation in the M one .
This difference persisted because of the strong reproductive barriers between the two forms relevantly pointed out by some authors [4, 26, 27]. In November 1999, the kdr mutation was identified for the first time in the M form at a very low frequency (0.006) with only one heterozygous RS among 161 specimens tested . The following year, the frequency in the M form increased to 0.02 and all individuals kdr-positive were only heterozygous (RS). Now, the proportion of homozygous resistant in the M form is increasing which enhances the overall frequency of this mutation.
The high frequency of the kdr mutation in the S form is presumably due to the long-term and extensive use of insecticide for cotton crop protection, DDT in the 1960–1970s replaced by pyrethroids in the 1980s . Then, its spreading from the S to the M form through introgression  is a recent and ongoing process, limited in savannah environments of West Africa to the place where An. gambiae M/Mopti and S/Savanna forms are found in sympatry at high densities .
Cotton crops are located just on the outside of rice fields. Insecticides applied on cotton during the rainy season may drift or be washed to the rice field during heavy rains conferring a selection pressure even for the M form. Moreover, during the last decade some rice producers also started to grow vegetables in the paddies using the same cotton's insecticides.
These agricultural uses of insecticides, mainly pyrethroids, exert the main selective pressure on the mosquito populations because of the quantity applied (e.g. six rounds of treatment each two weeks on cotton crop during the rainy season) and its action on the larval stage on large population of both sex. The domestic use of pyrethroids, through coils and ITNs, is more selective (it acts only on the anthrophilic fraction of biting females) and probably plays a secondary role on resistance selection in this rural area.
The present results showed that the ace-1Rmutation is mostly present in the S form and less frequently in the M one. This finding suggests that the ace-1Rresistance allele is evolving along the same pathway like the kdr mutation in this area. As for kdr, it occurred probably prior in the S form and may acquire by the M form through introgressive hybridization . The selection of the ace-1Rmutation in the S form could be related to the increasing use of organophosphates in cotton treatment in mixture with pyrethroids since the end of the 1990s. This insecticide resistance management (IRM) strategy was implemented at a large scale in West Africa to manage the pyrethroid resistance of the main cotton pest, Helicoverpa armigera .
With the exposure of the M populations to insecticide pressure, the ace-1Rmutation began to spread within this form. Similarly to the kdr mutation, it will probably increases in frequency within the M populations in the coming years. Additional studies are crucial to determine precisely the origin of this gene among the M form and gene flow pattern between the two forms in natural populations. Alternatively the reproductive fitness associated with this mutation in An. gambiae both S and M population remains to be evaluated .
The reported change in malaria vectors population structure is mainly driven by human activities and will call for modified malaria control strategies. The increasing of the S form proportion and the emergence of the ace-1Rmutation concomitantly with the Leu-Phe kdr mutation among the same populations of An. gambiae s.s. is an atypical ecological pattern in an irrigated rice growing area. With the expansion of agricultural practices such as vegetable growing, the application of pesticides with different mechanisms of action is rising. This may favour the development of multiple resistance mechanisms in An. gambiae [33, 34], which is a dynamic process that needs to be carefully monitored at the molecular form level and through designed spatial and seasonal surveys.
Further studies are needed to determine: (i) the phenotypic effects, particularly when the two mutations occur concomitantly and taking into account if metabolic-based resistance is present and (ii) the operational impact of both mutations on the efficacy of pyrethroid or organophosphate/carbamate based vector control.
Until recently, several studies conducted in savannah environment of Ivory Coast (West Africa) showed that pyrethroids treated nets still achieved a good control of kdr resistant An. gambiae either in experimental huts  or in field trials . Indeed recent paper from N'Guessan et al. in southern Benin  indicated that the kdr target insensitivity present at high frequency in M/Forest population of An. gambiae is associated with the decreased efficacy of ITNs and pyrethroid based IRS. Some alternatives to pyrethroids on ITNs are therefore necessary. Preliminary studies using organophosphate and carbamate treated nets in experimental huts have already shown good results in areas of kdr resistance [38, 39].
The presence of multiple resistance mechanism in An. gambiae in south-west Burkina Faso may constitute an obstacle for the future success of malaria control programmes based on ITNs or IRS with pyrethroids or organophosphates/carbamates. The present study should provide useful information for small and large-scale field trials on insecticide efficacy in this study area.
Coetze M, Craig M, Le Sueur D: Distribution of African malaria mosquitoes belonging to the Anopheles gambiae complex. Parasitol Today. 2000, 16: 74-77. 10.1016/S0169-4758(99)01563-X.
Coluzzi M, Petrarca V, Di Deco MA: Chromosomal inversion intergradation and incipient speciation in Anopheles gambiae. Bollet Zool. 1985, 52: 45-63.
Coluzzi M, Sabatini A, della Torre A, Di Deco MA, Petrarca V: A polytene chromosome analysis of the Anopheles gambiae complex. Science. 2002, 298: 1415-1418. 10.1126/science.1077769.
della Torre A, Fanello C, Akogbeto M, Dossou-yovo J, Favia G, Petrarca V, Coluzzi M: Molecular evidence of incipient speciation within Anopheles gambiae s.s. in West Africa. Insect Mol Biol. 2001, 10: 9-18. 10.1046/j.1365-2583.2001.00235.x.
Di Deco MA, Petrarca V, Villani F, Coluzzi M: Polimorfismo cromosomico da inversioni paracentriche e eccesso degli eterocariotipi in ceppi di Anopheles allevati in laboratorio. Parassitologia. 1980, 22: 304-306.
Persiani A, Di Deco MA, Petrangeli G: Osserzioni di laboratorio su polimorfismi da inversione originati da incroci tra popolazioni diverse di Anopheles gambiae s.s. Ann Istit Sup San. 1986, 22: 221-224.
Touré YT, Petrarca V, Traoré SF, Coulibaly A, Maïga HM, Sankaré O, Sow M, Di Deco MA, Coluzzi M: The distribution and inversion polymorphism of chromosomally recognised taxa of the Anopheles gambiae complex in Mali, West Africa. Parassitologia. 1998, 40: 477-511.
Favia G, Lanfrancotti A, Spanos L, Sideén-Kiamos I, Louis C: Molecular characterization of ribosomal DNA polymorphisms discriminating among chromosomal forms of Anopheles gambiae s.s. Insect Mol Biol. 2001, 10: 5-3. 10.1046/j.1365-2583.2001.00236.x.
Touré YT, Petrarca V, Traoré SF, Coulibaly A, Maïga HM, Sankaré O, Sow M, Di Deco MA, Coluzzi M: Ecological genetic studies in the chromosomal form Mopti of Anopheles gambiae in Mali, West Africa. Genetica. 1994, 9: 213-223. 10.1007/BF01443435.
della Torre A, Tu Z, Petrarca V: On the distribution and genetic differentiation of Anopheles gambiae s.s. molecular forms. Insect Biochem Mol Biol. 2005, 35: 755-69. 10.1016/j.ibmb.2005.02.006.
Diabaté A, Baldet T, Chandre F, Akogbeto M, Darriet F, Brengues C, Guiguemdé TR, Guillet P, Hemingway J, Hougard JM: 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-622.
Diabaté A, Baldet T, Chandre F, Dabiré KR, Kengne P, Simard F, Guiguemdé TR, Guillet P, Hemingway J, Hougard JM: Kdr mutation, a genetic marker to assess events of introgression between the molecular M and S forms of Anopheles gambiae (Diptera: Culicidae) in the tropical savannah area of West Africa. J Med Entomol. 2003, 40: 195-198.
Robert V, Gazin P, Boudin C, Molez JF, Ouedraogo V, Carnevale P: La transmission du paludisme en zone de savane arborée et en zone rizicole dans les environs de Bobo-Dioulasso, Burkina Faso. Ann Soc Belge Med Trop. 1985, 65: 201-214.
Baldet T, Diabaté A, Guiguemdé TR: Etude de la transmission du paludisme en 1999 dans la zone rizicole de la Vallée du Kou (Bama) Burkina Faso. Cahiers Santé. 2003, 13: 56-60.
Diabaté A, Chandre F, Rowland M, N'Guessan R, Duchon S, Dabiré KR, Hougard JM: The indoor use of plastic sheeting pre-impregnated with insecticide for control of malaria vectors. Trop Med Int Health. 2006, 11: 597-603. 10.1111/j.1365-3156.2006.01605.x.
Chandre F, Darriet F, Manga L, Akogbeto M, Faye O, Mouchet J, Guillet P: Status of pyrethroid resistance in Anopheles gambiae sensu lato. Bull World Health Organ. 1999, 77: 230-234.
Djogbénou L, Dabiré R, Diabaté A, Kengne P, Akogbéto M, Hougard JM, Chandre F: Identification and geographic distribution of the ace-1Rmutation in the malaria vector Anopheles gambiae in south-western Burkina Faso, West Africa. Am J Trop Med Hyg. 2008, 78: 298-302.
Weill M, Lutfalla G, Mogensen K, Chandre F, Berthomieu A, Berticat C, Pasteur N, Philips A, Fort P, Raymond M: Comparative genomics: Insecticide resistance in mosquito vectors. Nature. 2003, 423: 136-137. 10.1038/423136b.
Gillies MT, De Meillon B: The Anophelinae of Africa South of the Sahara. South Africa Institute of Medical Research Johannesburg, South Africa. 1968
Scott JA, Brogdon WG, Collins FH: Identification of single specimens of Anopheles gambiae complex by polymerase chain reaction. Am J Trop Med Hyg. 1993, 49: 520-529.
Martinez-Torrez D, Chandre F, Williamson MS, Darriet F, Berge JB, Devonshire AL, Guillet P, Pasteur N, Pauron D: Molecular characterisation of pyrethroid knockdown resistance (kdr) in the major malaria vector Anopheles gambiae s.s. Insect Mol Biol. 1998, 7: 179-184. 10.1046/j.1365-2583.1998.72062.x.
Weill M, Malcolm C, Chandre F, Mogensen K, Berthomieu A, Marquine M, Raymond M: The unique mutation in ace-1 giving high insecticide resistance is easily detectable in mosquito vectors. Insect Mol Biol. 2004, 13: 1-7. 10.1111/j.1365-2583.2004.00452.x.
Wirtz RA, Zavala F, Charoenvit Y, Campbell GH, Burkot TR, Schneider I, Esser KM, Beaudoin RL, Andre RG: Comparative testing of monoclonal antibodies against Plasmodium falciparum sporozoites for ELISA development. Bull World Health Organ. 1987, 65: 39-45.
Touré YT, Petrarca V, Traoré SF, Coulibaly A, Maïga HM, Sankaré O, Sow M, Di Deco MA, Coluzzi M: Ecological genetic studies in the chromosomal form Mopti of Anopheles gambiae s.s. in Mali, West Africa. Genetica. 1994, 94: 213-223. 10.1007/BF01443435.
Touré YT, Petrarca V, Traoré SF, Coulibaly A, Maïga HM, Sankaré O, Sow M, Di Deco MA, Coluzzi M: The distribution and inversion polymorphism of chromosomally recognized taxa of the Anopheles gambiae complex in Mali, West Africa. Parassitologia. 1998, 40: 477-511.
Fanello CM, Akogbeto M, della Torre A: Distribution of the pyrethroid knock down resistance (kdr) in Anopheles gambiae s.l. from Benin. Trans R Soc Trop Med Hyg. 2000, 94: 112-
Gentile G, Slotman M, Ketmaeir V, Powell JR, Caccone A: Attempts to molecularly distinguish cryptic taxa in Anopheles gambiae s.s. Insect Mol Biol. 2001, 10: 25-32. 10.1046/j.1365-2583.2001.00237.x.
Weill M, Chandre F, Brengues C, Manguin S, Akogbeto M, Pasteur N, Guillet P, Raymond M: The kdr mutation occurs in the Mopti form of Anopheles gambiae s.s. through introgression. Insect Mol Biol. 2000, 9: 451-455. 10.1046/j.1365-2583.2000.00206.x.
Diabaté A, Brengues C, Baldet T, Dabiré KR, Hougard JM, Akogbeto M, Kengne P, Simard F, Guillet P, Hemingway J, Chandre F: The spread of the Leu-Phe kdr mutation through Anopheles gambiae complex in Burkina Faso: genetic introgression and de novo phenomena. Trop Med Int Health. 2004, 9: 1267-73. 10.1111/j.1365-3156.2004.01336.x.
Djogbenou L, Chandre F, Berthomieu A, Dabiré R, Koffi A, Alout H, Weill M: Evidence of Introgression of the ace-1RMutation and of the ace-1RDuplication in West African Anopheles gambiae s.s. PLoS One. 2008
Martin T, Ochoub GO, Djihinto A, Traore D, Togola M, Vassal JM, Vaissayre M, Fournier D: Controlling an insecticide-resistant bollworm in West Africa. Agr Eco Environment. 2005, 107: 409-411. 10.1016/j.agee.2004.11.006.
Okoye PN, Brooke BD, Hunt RH, Coetzee M: Relative developmental and reproductive fitness associated with pyrethroid resistance in the major southern African malaria vector, Anopheles funestus. Bull Entomol Res. 2007, 97: 599-605. 10.1017/S0007485307005317.
Chouaïbou M, Etang J, Brévault T, Nwane P, Hinzoumbe CK, Mimpfoundi R, Simard F: 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. 2008, 13 (4): 476-486.
Corbel V, N'Guessan R, Brengues C, Chandre F, Djogbenou L, Martin T, Akogbéto M, Hougard JM, Rowland M: Multiple insecticide resistance mechanisms in Anopheles gambiae and Culex quinquefasciatus from Benin, West Africa. Acta Tropica. 2007, 101: 207-216. 10.1016/j.actatropica.2007.01.005.
Darriet F, N'Guessan R, Koffi AA, Konan L, Doannio JM, Chandre F, Carnevale P: Impact of pyrethrin resistance on the efficacy of impregnated mosquito nets in the prevention of malaria: results of tests in experimental cases with deltamethrin SC. Bull Soc Path Exot. 2000, 93: 131-134.
Henry MC, Assi SB, Rogier C, Dossou-Yovo J, Chandre F, Guillet P, Carnevale P: Protective efficacy of lambdacyhalothrin treated nets in Anopheles gambiae pyrethroid resistance areas of Côte d'Ivoire. Am J Trop Med Hyg. 2005, 73: 859-864.
N'Guessan R, Corbel V, Akogbeto M, Rowland M: Reduced efficacy of insecticide treated nets and indoor residual spraying for malaria control in pyrethroid resistance area, Benin. Emerg Infect Dis. 2007, 13 (2): 199-206.
Kolaczinski JH, Fanello C, Hervé JP, Conway DJ, Carnevale P, Curtis CF: Experimental and molecular genetic analysis of the impact of pyrethroid and non-pyrethroid insecticide impregnated bednets for mosquito control in an area of pyrethroid resistance. Bull Entomol Res. 2000, 90: 125-132.
Asidi AN, N'Guessan R, Koffi AA, Curtis CF, Hougard JM, Chandre F, Corbel V, Darriet F, Zaim M, Rowland MW: Experimental hut evaluation of bednets treated with an organophosphate (chlorpyrifos-methyl) or a pyrethroid (lambdacyhalothrin) alone and in combination against insecticide-resistant Anopheles gambiae and Culex quinquefasciatus mosquitoes. Malar J. 2005, 4: 25-30. 10.1186/1475-2875-4-25.
This study was financially supported by a grant of TDR/RCS re-entry project n°A30821 and a partial contribution of Sumitomo via the "Anopheles Biology and Control (ABC) network. We are grateful to villagers of Vallée du Kou who made this study possible. We thank Professor Martin Akogbeto, Head of the Centre de Recherche Entomologique de Cotonou (CREC, Cotonou, Benin) for his support, Joseph Chabi and Sébastien Koudenoukpo for their help during the genotyping. Finally we wish to thank Dr O. Skovmand from Intelligent Insect Control (IIC, Montpellier, France) for his helpful suggestions and correction made to the manuscript.
The authors declare that they have no competing interests.
DKR participated to the study design, undertook the field study, analysed the data and wrote the paper. DA participated to the study design, the data analysis and the manuscript drafting. LD participated to the study design and the samples analysis in the laboratory. OA participated to the field study and the sample collection. NR participated to the study design. OJB is the administrative authority who facilitated the implementation of the study. JMH designed the study, participated to the data analysis and the drafting of the paper. FC participated to the study design and the data analysis. TB participated to the data analysis and interpretation, the drafting and the revision to the paper. All authors read and approved the final manuscript.