Research | Open | Published:
The current distribution and characterization of the L1014F resistance allele of the kdr gene in three malaria vectors (Anopheles gambiae, Anopheles coluzzii, Anopheles arabiensis) in Benin (West Africa)
Malaria Journalvolume 18, Article number: 175 (2019)
The fight against malaria faces various biological obstacles, including the resistance of parasites to anti-malarial drugs and the resistance of mosquito vectors to insecticides. The resistance of Anopheles gambiae sensu lato (s.l.) to pyrethroids, the only class of insecticides used to impregnate mosquito nets, is known in Benin; the expansion of this resistance is influenced by the existence of gene flow between species, otherwise by the presence or absence of the kdr mutation in them. The objective of this study is to determine the spatial distribution of An. gambiae and the level of expression of the pyrethroid resistance kdr gene in seven agro-ecological zones of Benin.
The study was conducted in 18 localities belonging to seven agro-ecological zones where environmental parameters varied. The sites represent the main areas of eco-epidemiological malaria in Benin. Anopheles gambiae larvae were collected in natural breeding sites using ladles and dipping method and reared under standard conditions. These larvae were reared under standard conditions of temperature and humidity (26 to 30 °C and 60 to 90%) at the insectarium of the Centre de Recherche Entomologique de Cotonou (CREC). Adult female mosquitoes having emerged are morphologically and molecularly identified. Homozygous resistant (1014F/1014F), homozygous sensitive (1014L/1014L) and heterozygous (1014F/1014L) genotypes of the L1014F kdr gene mutation are determined by PCR.
A total of 677 An. gambiae was subjected at the PCR. The results revealed the presence of three vector species of the An. gambiae complex, of which 409 Anopheles coluzzii, 259 An. gambiae, 5 hybrids (An. coluzzii/An. gambiae) and 4 Anopheles arabiensis in the different agro-ecological zones. The four An. arabiensis were only found in Dassa, a locality in the cotton zone of central Benin. The frequency of distribution of the L1014F allele of the kdr gene varies from 84.48 to 100% in An. gambiae, from 80 to 100% in An. coluzzii and from 0 to 75% in An. arabiensis in the different agro-ecological zones. Moreover, a significant difference is generally observed in the distribution of the L1014F allele (P < 0.05). By comparing in pairs the distribution frequencies of this allele in the two species by agro-ecological zone, only a significant difference is noted in the central cotton and fishery zones (P = 0.0496).
In summary, even if the data are in small portions, the An. Arabiensis species was found only in central Benin and the L1014F allele of the kdr gene is widespread and seems to fix in all the species recorded in the different agro-ecological zones. This situation amplifies the problem of resistance, which could eventually be a significant obstacle for the malaria vectors control. Similarly, a study of their genetic structure via the L1014F allele is necessary in order to put in place strategies to manage this resistance. These strategies will take into account both the ecology and the genetic diversity of the organisms involved to preserve the effectiveness of pyrethroids, the only insecticides used for the impregnation of mosquito nets.
The malaria control faces various biological obstacles, including the resistance of parasites to anti-malarial drugs  and of mosquito vectors to insecticides . So far, the fight against this disease is essentially based on chemotherapy targeting the parasite in humans using anti-malarial drugs, as well as actions aiming at simultaneously reducing the human-vector contact, the density and the longevity of Anopheles vectors . As a result of years of efforts to control or even eradicate malaria, it still remains a global public health problem .
In the absence of a vaccine to immunize humans against harmful actions of the pathogen, vector control appears for the moment as the best preventive tool . Among the different methods that are commonly used are long-lasting insecticidal nets (LLINs) and indoor residual spraying (IRS) . With the increased use of insecticides from different classes to combat mosquitoes effectively, vector populations resistant to these chemicals have become increasingly important . The rationales are to be found in the maintenance of ancestral polymorphism and the vertiginous selection of resistance alleles that are, the most advantageous in the new environment created by the misuse of these insecticides. Indeed, several cases of emergence of vector resistance to insecticides have been described in Eastern  and Western [7,8,9,10] Africa.
Studies carried out in Burkina Faso , Benin  and Côte d'Ivoire  have shown that the frequency of the L1014F resistance allele of the kdr gene is higher in Anopheles gambiae sensu lato (s.l.) in agricultural areas usually treated with insecticides, as compared to rural areas where farmers grow only food or local consumption products. The hypothesis of the contribution of certain agricultural practices to the selection and extension of resistance in vectors cannot be excluded. In the current context where pyrethroids resistance is very widespread in Africa , studies on the diversity of this insecticide resistance gene in vector populations are important to improve the control tools used against An. gambiae, major vector of malaria.
Globally, the larvae of An. gambiae grow in small, shallow, relatively clean and sunny water reservoirs (puddles of water, stagnant water) . The increasing presence of An. gambiae s.l. in ecological sites, which were originally inappropriate and unsuitable for this species , leads to believe that this species has developed mechanisms that allow it to adapt to xenobiotics present in these new environments, the said mechanisms being characteristic of the normal biological evolution of these vectors. Clearly, the link between resistance and these new adaptations is still not formally established. Benin is a country of great geographical (climatic, relief) and ecological variability . It is, therefore, imperative to develop more molecular studies on the vectors that are encountered in order to propose more effective and sustainable control methods that take into account the genetic diversity within the An. gambiae complex according to the different ecological and geographical niches. In addition to the fact that the resistance of An. gambiae to pyrethroids, the only class of insecticides used for the impregnation of mosquito nets [7,8,9], is known and well documented in Benin [7,8,9], the expansion of this genetic and hereditary trait is strongly influenced by the existence of gene flows maintained through migration, preferential mating or inbreeding between different populations of the same species. In view of the dynamics that characterize populations of living organisms, it is important to focus on the biological evolution in time and space of An. gambiae sensu stricto (s.s.), An. coluzzii and An. arabiensis, three sibling species taken into account in the current study. Thus, we will to study the spatial distribution of these three species and of the resistance allele L1014F of the kdr gene within them in seven agro-ecological zones in order to better orient the control methods.
Characteristics of the study area
The current study was implemented in Benin, a West African country located in the intertropical zone between the equator and the Tropic of Cancer, more precisely between the parallels 6°30′and 12°30′ of latitude North on the one hand, and the meridian 1° and 3°40′ East longitude on the other hand. Benin has a subequatorial climate with four seasons (two rainy seasons interspersed by two dry seasons) in the south; a Sudano-Guinean transitional climate with four seasons similar to that of the South in the Centre of the country and, a Sudanese climate with two seasons (a rainy and a dry) in the north . The relief of Benin is not very rugged but has a low coastal and sandy region limited by lagoons. It also presents a plateau of ferruginous clay, a silico-clay plateau dotted with some undergrowth . In the north-west of the country is the massive Atacora with a height of 800 m and in the north-east, the silico–clay fertile plains of Niger . Benin has a large hydrographic network with two important basins namely: the Niger basin with several tributaries and the coastal basin whose rivers (Ouémé, Zou, Mono) reach the Atlantic Ocean. Each basin is dotted with small permanent or temporary water bodies that maintain agro-pastoral activities . Overall, eight agro-ecological zones have been defined in Benin , namely: the far north zone of Benin (Zone 1), the northern cotton zone of Benin (Zone 2), the food crop production zone of South-Borgou (Zone 3), Western Atacora Zone (Zone 4), Cotton Zone of Centre (Zone 5), bar land Zone (Zone 6), Depression Zone (Zone 7) and Fisheries Zone (Zone 8). Seven of the eight Zones are subject to this study. They are classified on the basis of homogeneity taking into account the consideration of climatic and agro-pedological parameters, cropping systems, population density and vegetation cover . In total, eighteen study sites were selected from seven different agro-ecological zones (Fig. 1). Thus, the districts of Dogbo, Houéyogbé, Toviklin, Bohicon and Porto-Novo have been selected in the bar land Zone. Those of Dassa, Savè, Ouèssè, Djidja and Kétou were selected in the cotton zone of central Benin. The districts of Malanville, Kandi and Copargo are, respectively, located in the zones of the far north, Northern cotton and Western Atacora. The districts of Kalalé and Péhunco are considered as being in the food crop production zone of South Borgou, while those of Cotonou, Comè and Sô-Ava are in the fishery zone.
Sampling and treatment of mosquitoes
Larvae collections and mosquito rearing
Larvae collections were conducted in all study sites from May 2015 to December 2017. They are made only during rainy seasons. Larvae of An. gambiae were collected from the positive breeding sites encountered using the dipping method. The latter consists in collecting mosquito larvae from the surface of the breeding sites through a ladle. The collected larvae were filtered and kept in labeled jars according to the surveyed sites and then transported to the insectary of the Center for Entomological Research of Cotonou (CREC). These larvae were raised under standard conditions of temperature (27 °C ± 2) and humidity (65% ± 10) until adult stage.
Molecular identification and species distribution
Adult female mosquitoes resulting from reared larvae were morphologically identified using morphological identification keys of Gillies and Coetzee  and of, Gillies and De Meillon . They are then subjected to the polymerase chain reaction (PCR). Thus, at least 30 mosquitoes were analysed by site. The SINE 200 PCR of Santolamazza et al.  and the PCR protocol of Scott et al.  were carried out for the identification of the species of the An. gambiae complex.
Determination of the L1014F mutation of the kdr gene in the three species
The resistant allele L1014F of the kdr gene was determined according to the protocol of Martinez-Torres et al. . This PCR diagnostic test consists in using four oligonucleotides or primers (Agd1, Agd2, Agd3, Agd4) and a Taq polymerase to search by amplification for resistant or susceptible alleles on a DNA fragment coding for the voltage-dependent sodium channel in each tested mosquito. The Agd1/Agd2 primers pair flanks the kdr gene by amplifying a 293 bp product as a control. The pair of Agd3/Agd1 primer pairs only with the resistance portion of the kdr gene to amplify a 195pb fragment. The Agd4 /Agd2 pair associates only with the portion of the sensitive gene by amplifying a 137 bp fragment. According to the protocol of Martinez-Torres et al. , the nucleotide sequences of these primers are: Agd1: 5′-ATAGATTCCCCGACCATG-3′; Agd2: 5′-ACAAGGATGATGAACC-3′; Agd3: 5′-AATTTGCATTACTTACGACA-3′; Agd4: 5′-CTGTAGTGATAGGAAATTTA-3'.
Analysis of the geographical distribution of the L1014F allele of the kdr gene in Anopheles gambiae, Anopheles coluzzii and Anopheles arabiensis
The frequencies of the L1014F allele of the kdr gene were determined by agro-ecological zone. They were calculated using the F-statistics formula  of the Genepop software version 4.2. The confidence intervals of the allelic frequencies were calculated using the binomial formula. The proportion test  of the R software, version 3.3.3 was used to compare the distribution of the different species and the of L1014F allele frequencies within these species of the various agro-ecological zones. It was therefore assessed whether the species and their resistance allele are distributed in the same way in the different zones.
Molecular identification and distribution of the different species
Table 1 provides results obtained from the molecular identification of species. After molecular species identification performed on An.Gambiae, the presence of three sibling vector species was revealed. They were distributed in variable proportions with a total of 409 An. coluzzii (60.41% ± 3.68), 259 An. gambiae s.s. (38.26% ± 3.66), 05 hybrids ((0.74% ± 0.64) An. coluzzii/An. gambiae) and 04 An. arabiensis (0.59% ± 0.58).
Overall, the frequency of distribution ranged from 0 to 96.67% in An. gambiae s.s., from 3.33 to 100% in An. coluzzii, and from 0 to 2.67% for An. arabiensis according to the agro-ecological zones.
Anopheles gambiae was found in the proportions of 96.67% in the cotton zone of north and food crop production zone of South Borgou and, 87.10%, 41.33%, 4.35% and 39.06% respectively in the Western Atacora, Central Cotton, bar land and fisheries zones. The distribution frequencies of An. coluzzii, was 3.33% in the Northern cotton and food crop production zones of South Borgou, 100%, 12.90%, 55.33%, 59.38% and 95.11% respectively in the Far North, Western Atacora, central cotton, bar land and fisheries zones. Anopheles arabiensis was absent in the different surveyed zones except for of the central cotton zone where it was 2.67%. A significant difference is noted in the distribution of sibling species within different agro-ecological zones (p < 0.001).
Anopheles gambiae and An. coluzzii live in sympatry in six zones except the northern extreme zone. Anopheles gambiae is widespread in northern Benin while, An. coluzzii was the major species from Central to South. It is observed that the distribution is more or less homogeneous with characteristic ecological preferences. Thus, An. gambiae is more adapted to the ecology of the North with high frequencies in the zones of cotton production of North (96.67% ± 6.42), of food crop production of South Borgou (96.67% ± 4.54) and Western Atacora (87.10% ± 11.80). Anopheles coluzzii prefers the southern zones with high frequencies in the central cotton (55.33% ± 7.96), bar land (59.38% ± 6.94) and fisheries (95.11% ± 3.12) zones even though we noted its exclusive presence in the far north. Anopheles arabiensis was found only in Dassa, the central cotton zone, but with a low number, which may be an ecological indicator.
Distribution of the frequency of the L1014F kdr mutation in An. coluzzii, An. gambiae and An. arabiensis
Table 2 presents the frequency distribution of the L1014F resistant allele of the kdr gene in An. coluzzii and An. gambiae s.s. in the different surveyed agro-ecological zones. Overall, the frequency of the L1014F resistant allele of the kdr gene ranged from 84.48 to 100% in An. gambiae and from 80 to 100% in An. coluzzii.
In An. gambiae, the frequencies were 0%, 84.48%, 91.38%, 96.30%, 93.55%, 99.33%, and 100% respectively in the Far North, Northern Cotton, food crop producing of South Borgou, West-Atacora, Central Cotton, laterite and Fisheries zones.
In An. coluzzii, the frequencies were 100% in the northern cotton and food crop production zones of South Borgou, and 80%, 87.5%, 85.54%, 87.72%, 88.86% and 87.29% respectively in the far north, western Atacora, central cotton, bar land and fisheries zones. On the other hand, a significant difference is generally observed in the distribution of the L1014F allele (p < 0.05). The frequencies of the L1014F allele of the kdr gene are similar between An. gambiae and An. coluzzii in all agro-ecological zones except for the central cotton (p = 0.0496) and of the fishery (p < 0.001) zones. The small number (n = 4) of An. arabiensis detected did not allow a valid comparison of the frequency of the L1014F allele of the kdr gene of this species with the one of An. coluzzii and An. gambiae.
Table 2 shows that the frequency of the L1014F allele of the kdr gene is at least 80% ± 10.12 in the different species of the different agro-ecological zones. In the same agro-ecological zone, there is no significant difference (p > 0.05) in the distribution of this allele in the two species except in the central cotton (p = 0.0496) and bar land zones (p < 0.001). Similarly, within the agro-ecological zones of the northern (far North, Western Atacora, North cotton and food crop production of South Borgou zones), no significant difference is noted in the distribution of this allele among the different species (p > 0.05). On the other hand, in the southern part of Benin (central cotton and bar land zones), there is a significant difference in the distribution of this allele among the different species (p < 0.05) except for the fishery zone (p = 0.3181).
The results of this study showed an heterogeneous distribution of An. gambiae and An. coluzzii in the different geographical areas of Benin and a strong sympatry between these two species in some localities. The presence of hybrids (An. gambiae /An. coluzzii) recorded in some study areas corroborates the results previously obtained by several authors [25,26,27] and supports the hypothesis of gene flow between the two species despite the well advanced speciation. The presence of hybrid suggests incomplete reproductive isolation between An. gambiae and An. coluzzii . On the other hand, the fertility and the reproductive capacity of these hybrids remain to be proven since their number is small.
Overall, the results show that An. coluzzii was the predominant species in the southern regions of the country while An. gambiae was the major species in the northern and central regions. This distribution of species could be due to variable climatic conditions (a subequatorial climate with 4-seasons, 2 rainy seasons interspersed with two dry seasons in the south; a Sudano-Guinean transition climate with four seasons similar to that of the South at the center and, a Sudanese climate with a wet and a dry season in the north ) as well as with the ecological characteristics (physico-chemical properties of the mosquito breeding sites) of the selected areas and the sampling periods depending on the dry or rainy season. According to Mbida et al. , An. coluzzii is associated to permanent breeding sites and those resulting from human activity, while An. gambiae is more frequent in rain-dependent temporary breeding sites. The same authors also concluded that An. coluzzii prefers urban water collections and adapts quickly to pollution. These results corroborate those obtained by Djogbénou et al. . The absence of An. gambiae in the agro-ecological zone of the far north could be due to the small number of surveyed localities in the only district of Malanville or to the physico-chemical properties of the breeding sites from which larvae were sampled .
In addition, only four individuals of An. arabiensis in Dassa, a savanna species previously found in northern Benin . According to Djogbénou et al.  and Gnanguènon et al. , An. arabiensis would have moved to the centre of the country because of drought and human activities (destruction of forests concomitantly with the use of large areas for agriculture, hunting, grazing and timber harvesting and urbanization). Thus, the fact of finding it in the centre of the country (Dassa) testifies to its adaptation to the new climatic conditions, an evolutionary capacity peculiar to living beings. Although the effects of climate change such as drought and rain are not yet clearly quantifiable, it could be one of the causes of the displacement of living organisms. Insects are one of the best indicators to better appreciate the phenomenon. Species migration, such as the hypothesis given in the case of An. arabiensis for the present study, may contribute to the spread of this vector in Central and Southern Benin.
Considering the different surveyed zones, no significant difference was noted between the frequencies of the resistant allele L1014F of the kdr gene in the populations of An. gambiae and An. coluzzii (p > 0.05) except for those the central cotton and bar land zones (p < 0.05). Benin is a country where the most practiced activity is agriculture with similar agricultural practices (the use of pesticides and chemical fertilizers) from one locality to another. This could contribute to the wide distribution of the L1014F locus. Also, the climatic conditions of the laterite and central cotton agro-ecological zones are close or even identical. On the other hand, these two zones do not share the same climatic conditions as the agro-ecological zones of the north (far north of Benin, food crop producing of South-Borgou and Western Atacora) which are themselves close or even identical. This could also be a justifying argument. According to Akogbéto et al.  and Yadouléton et al. [9, 12], agricultural practices constitute a pressure of selection which spread the resistance allele in Benin. National coverage of insecticide-treated mosquito nets in the country could also contribute to the selection of the L1014F resistant allele. The presence and attachment of the L1014F allele of the kdr gene in some populations of An. coluzzii, which was not the case a few years ago  are evolutionary realities in our study area. This could possibly be generalized very soon. This is proof of the increase of the receptive potential the L1014F allele of the kdr gene in An. gambiaes.l. And certainly by genetic self-stop of the expansion of resistance genes in Benin. With this spectrum of expression of the resistance gene, It is inevitably moving towards a lack of effectiveness of the insecticide-based tools currently used in the fight against malaria vectors. The presence of this allele in An. coluzzii differs from one author to another, and several have reported its presence only in An. gambiae [25, 32]. The existence of the L1014F allele of the kdr gene in An. coluzzii could also be due to an ancestral phylogenetic inheritance or introgression process within its populations [33, 34], or to the migration of An. coluzzii carrying the resistant allele kdr "L1014F" from one locality to another  given the proximity of some of our areas. These results corroborate those obtained in Benin, Burkina Faso  and Mali . The L1014F allele of the kdr gene is not directly distributed according to geography, but the ecology of the vectors in question is an indispensable factor in the selection of the resistant allele of the gene [9, 12, 13, 26]. A study of the genetic structure of these populations is necessary in order to explore their genetic dynamics in the different agro-ecological zones.
In summary, for the moment that even if the data are in small portions, the An. Arabiensis species was found only in central Benin and the L1014F allele of the kdr gene is widespread and seems to fix in all the species recorded in the different agro-ecological zones. This situation amplifies the problem of resistance, which could eventually be a significant obstacle for the malaria vectors control. Similarly, a study of their genetic structure via the L1014F allele is necessary in order to put in place strategies to manage this resistance. These strategies will take into account both the ecology and the genetic diversity of the organisms involved to preserve the effectiveness of pyrethroids, the only insecticides used for the impregnation of mosquito nets.
Availability of data and materials
The data supporting the conclusions of this article are included in the article. Raw data will be available upon request to the corresponding author.
Cotonou Entomological Research Center
Programme National de Lutte contre le Paludisme
Word Health Organization
Plowe CV. The evolution of drug-resistant malaria. Trans R Soc Trop Med Hyg. 2009;103:11–4.
Ranson H, N’Guessan R, Lines J, Moiroux N, Nkuni Z, Corbel V. Pyrethroid resistance in African anopheline mosquitoes: what are the implications for malaria control? Trends Parasitol. 2011;27:91–8.
WHO. World Malaria Report 2015. Geneva: World Health Organization; 2015.
Kelly-Hope L, Ranson H, Hemingway J. Lessons from the past: managing insecticide resistance in malaria control and eradication programmes. Lancet Infect Dis. 2008;8:387–9.
Hemingway J, Hawkes NJ, Mccarroll L, Ranson H. The molecular basis of insecticide resistance in mosquitoes. Insect Biochem Mol Biol. 2004;34:653–65.
Vulule JM, Beach RF, Atieli FK, Roberts JM, Mount DL, Mwangi RW. Reduced susceptibility of Anopheles gambiae to permethrin associated with the use of permethrin-impregnated bednets and curtains in Kenya. Med Vet Entomol. 1994;8:71–5.
Akogbéto M, Yakoubou S. Resistance of malaria vectors to pyrethroids used for impregnated bednets, Benin, West Africa (in French). Bull Soc Pathol Exot. 1999;92:123–30.
Corbel V, N’Guessan R, Brengues C, Chandre F, Djogbenou L, Martin T, et al. Multiple insecticide mechanisms in Anopheles gambiae and Culex quinquefasciatus from Benin, West Africa. Acta Trop. 2007;101:207–16.
Yadouleton AW, Asidi A, Djouaka RF, Braima J, Agossou CD, Akogbéto MC. Development of vegetable farming: a cause of the emergence of insecticide resistance in populations of Anopheles gambiae in urban areas of Benin. Malar J. 2009;8:103.
Chandre F, Darrier F, Manga L, Akogbeto M, Faye O, Mouchet J, et al. Status of pyrethroid resistance in Anopheles gambiae sensu lato. Bull World Health Org. 1999;77:230–4.
Diabaté A, Baldet T, Chandre F, Akogbéto M, Guiguemba TR, Darriet F, et al. The role of agricultural use of insecticides in resistance to pyrethroids in Anopheles gambiae sl in Burkina Faso. Am J Trop Med Hyg. 2002;67:617–22.
Yadouléton AW, N’Guessan R, Allagbé H, Asidi A, Boko M, Ossè R, et al. The impact of the expansion of urban vegetable farming on malaria transmission in major cities of Benin. Parasit Vectors. 2010;3:118.
Tia E, Chouaibou M, Gbalégba CNG, Boby AMO, Koné M, Kadjo AK. Distribution of species and kdr gene frequency among Anopheles gambiae ss and Anopheles coluzzii populations in five agricultural sites in Côte d’Ivoire (in French). Bull Soc Pathol Exot. 2017;110:130–4.
Carnevale P, Robert V. Les anophèles: biologie, transmission du Plasmodium et lutte antivectorielle. IRD Editions. 2009.
Kudom AA. Larval ecology of Anopheles coluzzii in Cape Coast, Ghana: water quality, nature of habitat and implication for larval control. Malar J. 2015;14:447.
Rapport national sur le développement humain. Agriculture, Sécurité alimentaire et Développement Humain au Bénin; 2015. p. 54–8.
Adam KS, Boko M. Le Bénin. Nouvelle édition mise à jour. Cotonou : Editions du Flamboyant, Paris : EDICEF, 1993. p. 95.
Gillies MT, Coetzee MA. Supplement to the Anophelinae of Africa South of the Sahara (Afrotropical Region). Publ South Afr Inst Med Res. 1987;1987:55.
Gillies MT, De Meillon D. The Anophelinae of Africa South of the Sahara (Ethiopian Zoogeographical Region). Publ South Afr Inst Med Res. 1968;54:1–343.
Santolamazza F, Mancini E, Simard F, Qi Y, Tu Z, Della Torre A, et al. Insertion polymorphisms of SINE200 retrotransposons within speciation islands of Anopheles gambiae molecular forms. Malar J. 2008;7:163.
Scott JA, Broodon GW, Collins HF. Identification of single specimens of the Anopheles gambiae complex by the polymerase chain reaction. Am J Trop Med Hyg. 1993;49:520–9.
Martinez-Torres D, Chandre F, Williamson MS, Darriet F, Berge JB, Devonshire AL, et al. Molecular characterization of pyrethroid knockdown resistance (kdr) in the major malaria vector Anopheles gambiae ss. Insect Mol Biol. 1998;7:179–84.
Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution. 1984;38:1358–70.
R Development Core Team: R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria 2011. ISBN 3-900051-07-0, https://www.R-project.org/.
Della Torre A, Fanello C, Akogbeto M, Dossou-yovo J, Favia G, Petrarca V, et al. Molecular evidence of incipient within Anopheles gambiae ss in West Africa. Insect Mol Biol. 2001;10:9–18.
Péka M. Evaluation de la sensibilité aux pyréthrinoïdes des populations d’Anopheles gambiaes.l. de la zone côtière de Côte d’Ivoire. Mémoire DEA: CEMV, Université de Bouaké; 2001. p. 57.
Djogbenou L, Pasteur N, Bio-Bangana S, Baldet T, Irish SR, Akogbéto M. Malaria vectors in the Republic of Benin: distribution of species and molecular forms of the Anopheles gambiae complex. Acta Trop. 2010;114:116–22.
Taylor C, Touré YT, Carnahan J, Norris DE, Dolo G, Traoré SF, et al. Gene flow among populations of the malaria vector, Anopheles gambiae, in Mali. West Africa. Genetics. 2001;157:743–50.
Mbida Mbida A, Etang J, Akono Ntonga P, Eboumbou Moukoko C, Awono-Ambene P, Tagne D, et al. New insight into Anopheles coluzzii Coetzee & Wilkerson, 2013 larval ecology in the Wouri estuary, Littoral-Cameroon (in French). Bull Soc de Pathol Exot. 2017;110:92–101.
Gnanguenon V, Govoetchan R, Agossa RF, Ossè R, Oke-Agbo F, Azondekon R, et al. Transmission patterns of Plasmodium falciparum by Anopheles gambiae in Benin. Malar J. 2014;13:444.
Akogbeto MC, Djouaka RF, Kinde–Gazard DA. Screening of pesticide residues in soil and water samples from agricultural settings. Malar J. 2006;5:22.
Euzéby J. Grand dictionnaire illustré de parasitology médicale et vétérinaire. Éditions Lavoisier. 2008; 832.
Diabate A, Baldet T, Chandre F, Dabire RK, Kengue P, Guiguemde TR, et al. Kdr mutation, a genetic marquer to asses events of introgression between the molecular M and S forms of Anopheles gambiae (Diptera: Culicidae) in the tropical savanna area of West Africa. J Med Entomol. 2003;40:195–8.
Touré YT, Petrarca V, Traoré SF, Coulibaly A, Maiga HM, Sankaré O, et al. The distribution and inversion polymorphism of chromosomally recognized taxa of the Anopheles gambiae complex in Mali, West Africa. Parassitologia. 1998;40:477–511.
Edi CA, Koudou BG, Bellai L, Adja AM, Chouaibou M, Bonfoh B, et al. Long-term trends in Anopheles gambiae insecticide resistance in Côte d’Ivoire. Parasit Vectors. 2014;7:500.
Weill M, Chandre F, Brengues C, Manguin S, Akogbeto M, Pasteur N, et al. The kdr mutation occurs in the Mopti form of Anopheles gambiae ss through introgression. Insect Mol Biol. 2000;9:451–5.
Cisse MB, Keita CC, Dicko A, Dengela D, Coleman J, Lucas B, et al. Characterizing the insecticide resistance of Anopheles gambiae in Mali. Malar J. 2015;14:327.
We thank the anonymous reviewers for their constructive comments on manuscript. We also thank the Cotonou Entomological Research Center (CREC), mainly the director for his support of the study.
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.