Multiple insecticide resistance in an infected population of the malaria vector Anopheles funestus in Benin
© The Author(s). 2016
Received: 16 December 2015
Accepted: 25 July 2016
Published: 17 August 2016
Knowledge on the spread and distribution of insecticide resistance in major malaria vectors such as Anopheles funestus is key to implement successful resistance management strategies across Africa. Here, by assessing the susceptibility status of an inland population of An. funestus Giles (Kpome) and investigating molecular basis of resistance, we show that multiple resistance and consistent plasmodium infection rate are present in Anopheles funestus populations from Kpome.
The insecticide susceptibility level of collected Anopheles funestus was assessed. Synergist (PBO) was used to screen resistance mechanisms. The TaqMan technique was used for genotyping of insecticide resistant alleles and detecting plasmodium infection levels. The nested PCR was used to further assess the plasmodium infection rate.
The TaqMan analysis of plasmodial infections revealed an infection rate (18.2 %) of An. funestus in this locality. The WHO bioassays revealed a multiple phenotypic resistance profile for An. funestus in Kpome. This population is highly resistant to pyrethroids (permethrin and deltamethrin), organochlorines (DDT), and carbamates (bendiocarb). A reduced susceptibility was observed with dieldrin. Mortalities did not vary after pre-exposure to PBO for DDT indicating that cytochrome P450s play little role in DDT resistance in Kpome. In contrast, we noticed, a significant increase in mortalities when PBO was combined to permethrin suggesting the direct involvement of P450s in pyrethroid resistance. A high frequency of the L119F-GSTe2 DDT resistance marker was observed in the wild DDT resistant population (9 %RS and 91 %RR) whereas the A296S mutation was detected at a low frequency (1 %RS and 99 %SS).
The presence of multiple resistance in An. funestus populations in the inland locality of Kpome is established in this study as recently documented in the costal locality of Pahou. Data from both localities suggest that resistance could be widespread in Benin and this highlights the need for further studies to assess the geographical distribution of insecticide resistance across Benin and neighboring countries as well as a more comprehensive analysis of the resistance mechanisms involved.
Malaria remains an important health issue in Benin where it is the main cause of morbidity and mortality with children under five and pregnant women being the more vulnerable groups . In this country as across Africa, long lasting insecticide nets (LLINs) and bendiocarb based indoor residual spraying (IRS) are the key tools used for malaria control . In addition to pyrethroids and DDT, organophosphates (Ops) and carbamates (like bendiocarb) are also used to some extent for IRS. However, the success of these control methods is jeopardized by the development by Anopheles species of resistance to insecticides such as pyrethroids and DDT which are used for mosquito control. Indeed, resistance against the main insecticides used in public health, such as pyrethroids, carbamates or DDT is increasingly reported in An. funestus Giles populations [2–5] from southern Africa (Mozambique  and Malawi [6, 7]); East Africa (Uganda) [8, 9]; Central Africa (Cameroon)  and West Africa (Ghana, Benin [4, 10]). In Benin, most insecticide susceptibility studies have been conducted on An. gambiae Giles with resistance to pyrethroids and DDT recorded in An. gambiae (sensu stricto) (s.s.) and An. arabiensis [11–14]. Anopheles gambiae (sensu lato) (s.l.) is a complex of sibling species including two main malaria vectors in sub-Saharan Africa (SSA): An. gambiae (s.s.) and An. arabiensis. As for An. gambiae (s.s.), two reproductive units previously known as “M” and “S” molecular forms were recently elevated to the species rank: the An. gambiae “S” is now designated as An. gambiae (s.s.) and the An. gambiae “M” form is now officially known as An. coluzzii . On the other hand, An. funestus, the other major malaria vector, has received very little attention. Anopheles funestus complex is a group of nine species . Of this complex, An. funestus (s.s.) is one of the main malaria vectors in the sub-Saharan Africa. Recent studies conducted in the coastal locality of Pahou in Benin revealed the presence of resistance of An. funestus to pyrethroids and DDT . Furthermore, Riveron et al.  demonstrated that DDT detoxification in An. funestus is mainly associated with a single mutation in the GSTe2 gene. Very little is known about the extent and spread of this resistance across Benin notably in the inland localities. Such information can neither be extrapolated from the recorded resistance patterns in coastal regions nor from the neighboring countries as the resistance profiles in An. funestus populations vary significantly across Africa. For example the resistance profiles observed in North Cameroon in 2007 (DDT and dieldrin resistance)  are different to those observed in southern African countries like Mozambique, Malawi and South Africa for pyrethroid, DDT or carbamate resistance [2, 7, 17, 18]. Pyrethroid resistance has been documented in South Africa . No cases of carbamate resistance have been found in this country. However, carbamate resistance has been well documented in the southern region of Mozambique around the border with South Africa . This resistance pattern in southern Africa is also different to that of East Africa (pyrethroid and DDT resistance but full susceptibility to carbamates)  or to that of Ghana (West Africa) (DDT resistance and pyrethroid resistance) . Mechanisms of resistance of An. funestus to pyrethroids, carbamates or DDT identified so far are mainly caused by the elevated expression of detoxification genes (metabolic resistance) as neither the L1014F kdr mutation nor the G119S Ace-1 mutation have been detected in this species [4, 8, 19]. However, the detection of A296S Rdl r mutation in the GABA receptor (involved in dieldrin resistance) of An. funestus indicates that target site resistance mechanism is also present in this species . Indeed, P450 genes have been found to be associated with pyrethroid resistance [8, 20, 21] and carbamates resistance as well [22, 23] while mechanisms for the DDT resistance detected in East Africa (Uganda) and West Africa (Benin) are mainly associated to the L119F and GSTe2 mutation . Benin is currently scaling up its malaria control program through Long Lasting Impregnated Nets (LLINs) and IRS . For a successful scaling up of activity, there is a great need to document the susceptibility profiles of malaria vectors to insecticides used in public health and the underlying mechanisms. Generated data will help guiding control programs on the selection of the most suitable insecticides to use and will facilitate the design of country tailored management strategies. However, complementary studies are required to assess the operational impact of the resistance on the efficacy/efficiency of vector control tools against malaria transmission and mortality/morbidity in local populations before giving any clear recommendation to malaria control programs.
In this study, we report the assessment of susceptibility levels to several insecticides used in public health of An. funestus population from the locality of Kpome, an inland locality in the southern Benin. We analyzed the infection rate with Plasmodium species in this malaria vector population and, investigated the underlying resistance mechanisms. This research will provide current information on the plasmodium infection rate of An. funestus and, the insecticide resistance patterns in this species These information will help to improve the control of malaria vectors in Benin.
Area of study and mosquito collection
Description of study site
Electric aspirators and torches were used to collect blood-fed An. funestus adult females resting indoors. Collections were conducted in houses between 06:00 and 10:00 am in Kpome. The collection was carried out 2 weeks per month and 30 to 40 rooms were aspirated each sampling day for the period between December 2013 and February 2014. The target sampling period corresponds to the dry season in the southern Benin when An. funestus densities are likely to be higher. Blood-fed and gravid females of An. funetus (F0) were kept in small cups until fully gravid. Eggs were obtained from F0 generation using the forced egg laying method . Eggs were sent via courier to the Liverpool School of Tropical Medicine (LSTM) insectary where they were allowed to hatch in small plastic cups with approximately 100 ml of water. L1 larvae emerging from eggs were later transferred into bowls (with 800–1000 ml of mineral water) for rearing and the production of the F1 generation. Larvae rearing was conducted under standard insectary conditions (26 ± 2 °C with a relative humidity of 80 %). The egg batch and emerging larvae from the same female mosquito were reared together and later pooled with larvae from other females if these females were found belonging to the same molecular species (F0 females underwent PCR species identification for facilitating larvae pooling). Larvae were fed daily with Tetramin™ baby fish food and the water of each larvae bowl was changed every two days to reduce the mortality rate. The F1 adults generated were randomly mixed in cages for subsequent experiments.
PCR-based species identification
All females used for individual oviposition were morphologically identified as belonging to the An. funestus group . Of these females, we subjected a total of 104 (about 75 % of females that oviposited) to DNA extractions using DNeasy extraction kits (Qiagen DNeasy kit) followed by PCR species identification. A PCR was carried out using the protocol described by Koekemoer et al.  for molecular identification of collected females. Anopheles gambiae collected in the same locality during same period were identified to species and molecular form using the Livak method  and the SINE-PCR protocol 
Infection of An. funestus with plasmodial species
The plasmodium infection rate was determined using the TaqMan assay described by Bass et al. . The reaction was performed in 10 μl final volume containing 1 × SensiMix (Bioline), 800 nM of each primer described by Bass et al.  and 200 nM of each probe labelled with fluorophores, FAM to detect Plasmodium falciparum, and HEX to detect P. ovale, P. vivax and P. malariae (P. ovm). Two P. falciparum samples and a mix of P. ovale, P. vivax and P. malariae were used as positive controls. The real-time PCR MX 3005 (Agilent, Santa) system was used for amplification with the following cycling conditions: 95 °C for 10 min for denaturation followed by 40 cycles of 15 s at 92 °C and 1 min at 60 °C. The same DNA material extracted from the indoor-caught An. funestus used for PCR-based species identification was used for the TaqMan Plasmodium infection analysis. Additionally, the infection rate of An. coluzzii was also estimated for comparison. A nested PCR was performed on all Taqman-Plasmodium positive samples to validate results from the TaqMan assay described by .
Insecticide susceptibility assays
Protocols and standard insecticide treated papers supplied by WHO  were used to test the susceptibility of An. funestus from Kpome to various insecticides. We assessed the susceptibility pattern to six insecticides belonging to the four major public health classes of insecticide: the pyrethroids type I permethrin (0.75 %) and type II deltamethrin (0.05 %); the organochlorines DDT (4 %) and dieldrin (4 %); the organophosphate malathion (5 %) and the carbamate bendiocarb (0.1 %). Knockdown was recorded after 1 h and a 10 % sugar solution was made available to survivors. Final mortality was scored 24 h post-exposure. Because of the absence of susceptible strains of An. funestus, the susceptible strain An. gambiae kisumu and the wild population of An. funestus were exposed to insecticide-treated and non-treated papers as controls . Prior to the experiment, we also confirmed the effectiveness of insecticide-treated papers by exposing the susceptible strain An. gambiae kisumu to insecticide-impregnated papers. A mortality in the range 98–100 % indicates susceptibility. A mortality of less than 98 % is suggestive of the existence of resistance and further investigation is needed. If the observed mortality (corrected if necessary) is between 90 and 97 %, the presence of resistant genes in the vector population must be confirmed through additional bioassay tests and/or by conducting molecular assays. If mortality is less than 90 %, confirmation of the existence of resistant genes in the test population with additional bioassays may not be necessary, as long as a minimum of 100 mosquitoes was tested .
To further assess the extent of the susceptibility levels to pyrethroids, the F1 population from Kpome was further exposed to permethrin for 90 min (30 min additional time). Due to the limited number of females, it was not possible to extend the exposure period for other insecticides. In addition to females, the level of insecticide susceptibility was also assessed in males of An. funestus. Each test included control mosquitoes exposed to untreated papers.
Following the high level of resistance recorded against 4 % DDT and 0.75 % permethrin, we tested the effect of the synergist PBO in combination with these insecticides to assess the potential role of monooxygenase enzymes. We assessed the synergist PBO because of previous reports of P450s involvement in DDT resistance in An. gambiae  and pyrethroid resistance in An. funestus . Fifty female and 50 male mosquitoes were pre-exposed to 4 % PBO paper for 1 h and immediately exposed to 4 % DDT and 0.75 % permethrin for 1 h. Mortalities following these combined exposures were recorded after 24 h and compared with the results obtained without PBO. Two controls were used for this experiment, one control with mosquitoes exposed to untreated papers and the second control with mosquitoes exposed to the synergist PBO only. Due to the limited number of emerging mosquitoes (F1), we could not test other synergists (e.g. esterase and GST) and other insecticides apart from DDT and permethrin (e.g. deltamethrin, carbamate etc.).
Genotyping of resistance markers L119F-GSTe2 and A296S-RDL
To assess the role of the L119F-GSTe2 in DDT resistance  and the A296S-RDL mutation conferring dieldrin resistance [3, 7], a TaqMan genotyping assay for these markers was used. DNA extracts from adult female An. funestus collected indoors (F0) were separately genotyped for the GSTe2 and the Rdl r allelic variants. In addition, female mosquitoes with known DDT susceptibility phenotypes (F1, dead post-exposure), as defined by the standard WHO protocol  were also screened for GSTe2. Taqman SNP genotyping assays were performed in 10 μl volume containing 1× Sensimix (Bioline), 80× primer/probe mix and 1 μl template DNA. Probes were labelled with two specific fluorophores FAM and HEX, FAM to detect the homozygous resistant genotype, HEX to detect the homozygous susceptible genotype and both FAM and HEX to detect the heterozygous genotype. The assay was performed on an Agilent MX3005 real-time PCR machine with cycling conditions of 95 °C for 10 min, followed by 40 cycles at 95 °C for 15 s and 60 °C for 1 min. FAM and HEX fluorescence was captured at the end of each cycle and genotypes called from endpoint fluorescence using the Agilent MXPro software.
Chi-square was used to test for significant differences in the percentage mortalities between female and male mosquitoes exposed to insecticides and Statistica package was used for TaqMan analysis. Genotype frequencies were confirmed to be according to Hardy-Weinberg equilibrium using the Had2know online statistical software (http://www.had2know.com/academics/hardy-weinberg-equilibriumcalculator-2-alleles.html).
Species identification and distribution
Results from PCR-based species identification of 104 morphologically identified females of An. funestus (s.l.) collected from Kpome during the dry season (from December 2013 to February 2014) revealed that they all belong to An. funestus (s.s.).
Plasmodium infection rates
Infectivity of An. funestus ss and An. coluzzii from Kpome to Plasmodium species (falciparum, ovale, vivax and malariae)
Infectivity using the Taqman protocol
Infectivity of An. funestus using the nested PCR protocol
+ve P. fal
+ve P. ovm
+ve P. fal & P. ovm
Total no. of infected Anopheles with Taqman (% positivity)
Total no. of infected An. funestus using nested PCR (12 out the 17 to Taqman)
+ve P. fal
+ve P. ov
+ve P. viv
+ve P. mal
An. funestus (s.s.)
17 ± 7.85 (18)
12 ± 6.83 (13)
Susceptibility to insecticides
Synergist assay with PBO
Genotyping of the GSTe2 alleles in the An. funestus population of Kpome
Following the observed high resistance profile of An. funestus Kpome to DDT, we screened for the GSTe2 mutation which was recently shown to confer DDT resistance in the coastal An. funestus population from Pahou in the Benin . Results showed a high presence of the L119F-GSTe2 mutation in the wild An. funestus population from Kpome with 91 % individuals homozygote for the 119 F resistant allele (84 out of 92 individuals tested) whereas 9 % were RS. The recorded genotypic frequency was found at Hardy-Weinberg (HW) equilibrium with a P-value of 0.6628 and the expected HW frequencies of 84.17 RR, 7.65 RS and 0.17 SS. No SS individual was recorded in Kpome showing that the GSTe2 is nearing fixation in the An. funestus population of Kpome as previously observed in Pahou  (Additional file 1: Figure S1).
Association of the L119F-GSTe2 mutation with DDT resistance
Genotyping of the A296S Rdl alleles in Kpome
Despite the relative susceptibility of the Kpome population to dieldrin, we genotyped the A296S mutation to confirm whether its frequency correlates with observed susceptibility in Kpome. Results showed a very low presence of Rdl r in the genotyped F0 population of Kpome; out of 92 individuals genotyped, 91 were SS (99 %) and only one was RS (1 %). No individual was recorded with RR mutation showing that the selection mechanism for this mutation in wild An. funestus populations of Kpome is still very low (Additional file 2: Figure S2). This low allelic frequency of Rdl r resistant allele in the population further correlates with the high susceptiblity of the An. funestus (s.s.) in Kpome to dieldrin (Fig. 2). The genotypic frequency of Rdl was found at Hardy-Weinberg equilibrum with a P-value of 0.9582 and the expected frequencies of 0 RR, 0.99 RS and 91 SS.
The susceptibility profile analysis of an inland population of An. funestus has revealed for the first time in Benin that multiple resistance patterns extent beyond the coast; the Kpome population having demonstrated resistance to several insecticide classes in our study.
Contribution of An. funestus to malaria transmission
This study further confirms the presence of An. funestus in the southern Benin as recently published by Djouaka et al.  and highlights its consistent infection level with Plasmodium species. The Kpome An. funestus population was found with consistent Plasmodium infections during the dry season (infection rate of 18.2 %). It is worth mentioning that the analysis of Plasmodium spp. in An. funestus also revealed that some mosquitoes could be co-infected with more than two species, namely P. falciparum and any of the P. ovale, P. malaria or P. vivax. When the infection rate of An. funestus was compared with that of An. coluzzii mosquitoes collected in Kpome during a similar period (dry season), we noticed a relatively high infection rate of indoor-collected An. funestus compared to indoor-collected An. coluzzii. This comparative analysis shows the high competitiveness of An. funestus in the transmission of malaria in southern Benin. The infection rates (infection rates > 15 %) recorded in both species (An. funestus and An. coluzzii) during the same period predicts possibilities of high malaria transmission and incidence in this locality during dry seasons. When similar samples were analyzed with the nested PCR, a relatively low number of infected mosquitoes was recorded; this certainly reveals a higher sensitivity of the Taqman technique compared to the nested PCR technique . The infection rate recorded in this An. funestus population is as high as previously recorded in other localities across Africa. For example, parasite rates of 22 %  and 27 %  have been recorded for this species in the past in South Africa and more recently, in West Africa, mean rates of infectivity between 3 % and 15 % have been observed in Burkina Faso [35, 36]. During the same studies, Costantini et al.  recorded up to 50 % An. funestus in one village in November 1991 (n = 56) positive for P. falciparum circumsporozoite and Dabire et al.  recorded up to 20 % of infected An. funestus in Lena in August 2000 (n = 40) using the same Elisa method. In southern Benin, a recent study using multiplex real-time PCR assays revealed P. falciparum infection and mixed infections with P. malariae and/or P. ovale in 13.6 % of An. funestus collected . Overall, high infection rates of An. funestus observed in some villages in West Africa in particular at the end of the rainy season and its related high anthropophily, contribute to the main role of this species in malaria transmission, surpassing in some cases other main vectors like An. gambiae (s.s.). The implication of An. funestus in malaria transmission during the dry season was also documented in Ghana , Nigeria , Burkina Faso [35, 36] and more recently, in southern Benin . No Plasmodium vivax was found in analysed mosquitoes from Benin. The infection rate recorded with the Taqman technique was higher than that obtained with the nested PCR; this may be due to the higher sensitivity of the Taqman assay compared to the nested PCR technique .
Multiple resistance could be widespread in An. funestus in Benin
The multiple resistance pattern of the An. funestus population in Kpome is similar to that of the coastal location of Pahou  suggesting that multiple resistance could be widespread in An. funestus populations in Benin at least in the southern region. This multiple resistance concerns several insecticide classes including the organochlorine DDT, pyrethroids and carbamates. The An. funestus population is highly resistant to both type I and II pyrethroids. The level of this resistance is consistently higher than that reported back in 2011 for the coastal population of Pahou. For example for females, the mortality rate after permethrin exposure was 66.7 % in Pahou in 2011 but only 13 % in Kpome in 2014. Similarly, for the type II pyrethroid, deltamethrin, a higher mortality rate of 88.6 % was observed in Pahou compared to only 46.5 % in Kpome. This difference could suggest an increase in the overall level of pyrethroid resistance in Benin in three years. This increase could be related to the increased ITNs coverage across Benin and additional selection factors such as agricultural use of pesticides [41, 42], and spilled petroleum products . High resistance was recently reported for this species in Malawi where a significant rise of levels of pyrethroid resistance was detected when comparing a 2014 population of southern Malawi to the 2009 population from the same location . Same increase was also reported in Uganda for An. funestus  suggesting that pyrethroid resistance is being further selected in field populations and this is a concern for the continued effectiveness of pyrethroid-based interventions such as LLINs. This concern is further compounded by the increase in resistance also observed in An. gambiae (s.l.) populations from Benin [14, 44], Burkina Faso and Cameroon [45–47].
The resistance recorded in Kpome is higher against permethrin (type I pyrethroid) than to deltamethrin (type II pyrethroid). This resistance pattern is different to that observed in southern or eastern Africa where resistance to type II pyrethroid is higher than to type I [6, 9, 19, 48]. This difference may underline the existence of some differences in the resistance mechanism against pyrethroids in Benin compared to the southern and Eastern Africa regions. This could be explained by the significant association established between the L119F-GSTe2 mutation and permethrin resistance but not for deltamethrin . Moderate permethrin resistance was also detected in other West African An. funestus populations from Obuasi in Ghana  but no pyrethroid resistance was reported in An. funestus from Soumossou in Burkina Faso in 2007, although it may have changed by now  indicating that this resistance may not yet be widely distributed across West Africa. This resistance to pyrethroids is of great concern for malaria control programs and, there is a risk that if such resistance is not managed properly, it can be further selected by ongoing control interventions such as the pyrethroid impregnated LLINs to a level that will seriously impact the success of future control programs. The current trend of pyrethroid resistance recorded with the main malaria vectors (An. gambiae and An. funestus) in Benin [4, 13, 14, 49] is a concern for the success of malaria control program focussed on pyrethroid-based LLINs. Beyond the evaluation of pyrethroid resistance profile of malaria vectors, a pilot Phase II entomological study examining the effectiveness of using ITNs at two sites in Benin has already given clear evidence of pyrethroid-based ITNs failing to control a pyrethroid-resistant An. gambiae population . Nevertheless, at the epidemiological level, a recent study performed in the South Benin revealed that the observed pyrethroid resistance of malaria vectors seems to date not to have affected the impact of LLINs; therefore, the use of LLINs remains associated with reduced malaria prevalence irrespective of resistance . The DDT resistance level observed in this population is one of the highest in Africa after the one documented at Pahou where no mortality was recorded after 1 h exposure . This confirms that DDT resistance level is high in Benin in contrast to southern Africa where such resistance has just been reported in Zambia and Malawi [7, 22]. The level of DDT resistance in Benin is also higher than the level observed in East Africa where mortality rates are above 40 % [9, 48]. The presence of DDT resistance within Benin is certainly going to reduce the options available for insecticide resistance management of An. funestus populations particularly with the observed resistance to pyrethroids and carbamates. This multi-resistance of An. funestus species to several groups of insecticides brings more complexities in the control of malaria as this high DDT resistance is also prevalent in several populations of An. gambiae in Benin [4, 49]. The presence of insecticide-resistant populations of An. funestus in costal urbanized and polluted localities (e.g. Pahou) and in the inlands (e.g. Kpome) with less urbanization and less pollution could be associated with gene flow among An. funestus populations in southern Benin. In addition, exogenous factors such as agriculture, urbanization (pollution including petroleum products) and ITNs/IRS (the use of public health insecticides) might have also contributed to the observed high resistance profile of mosquitoes [4, 11, 43, 52–54]. In the coastal localities of Benin all exogenous factors (ITNs/IRS, urban agriculture with high pesticide utilization, rapid urbanization/urban pollution and spilled petroleum products from vehicles, repair shops and road-side re-sellers) simultaneously act on mosquito population and intensely select for insecticide resistance, whereas in the inland localities the main dominant factor for selection of DDT resistance is the agricultural use of pesticides. An. funestus population from Kpome is also resistant to the carbamate bendiocarb. This observed bendiocarb resistance at Kpome confirms the presence of bendiocarb resistance in An. funestus in the southern Benin as earlier reported by Djouaka et al.  in Pahou. Also, the bendiocarb resistance was recently reported in An. gambiae (s.s.) in Atacora, in the northern part of Benin . This bendiocarb resistance recorded in this study raises a special concern for National Malaria Control programs which, because of high resistance to pyrethroids and DDT, are currently introducing bendiocarb based IRS for malaria vector control in West African countries . Such resistance was observed in An. gambiae population in Benin three years after the large implementation of bendiocarb based IRS . The rapid selection of bendiocarb resistance in An. gambiae was also observed in other West African countries, e.g. Burkina Faso following the introduction of IRS by the US PMI program . The same process probably affected An. funestus considering its high endophilic behavior. The full susceptibility to the organophosphate malathion is similar to results throughout Africa showing that An. funestus populations remain susceptible to this insecticide class which can therefore be used as alternative insecticide for resistance management strategies. However, the toxicity of malathion might limit its adoption by malaria vector control programs. Recently, other less toxic OP candidates like pirimiphos methyl  or products belonging to other classes of insecticide like the pyrrole insecticide chlorfenapyr  and the indoxacarb oxadiazine  revealed promising in Phase II studies against pyrethroid-resistant malaria vectors in Benin and elsewhere in sub-Saharan Africa.
The nearly full susceptibility to dieldrin correlated very well to the extremely low frequency of the Rdl r mutation in is similar to previous observation in Pahou  but contrasts to the high frequency observed in neighboring Burkina Faso (Vallée du Kou)  and Cameroon . In Burkina Faso, insecticide susceptibility tests performed in 1967 in the Soumousso village already reported resistance to dieldrin in An. funestus,  highlighting that insecticide resistance selected in the past could persist nowadays even in the absence of selection pressure (i.e. contemporary use of the same or related insecticide) in particular for resistance mechanisms that have no fitness cost for the insect.
Analysis of underlying mechanisms of the observed multiple resistance profiles
The underlying mechanisms of the multiple resistance pattern observed in this population was explored through various means. The synergist assay with PBO, an inhibitor of cytochrome P450 monooxygenases, indicated that this enzyme family does not play a major role in the high DDT resistance observed in Kpome. However, when mosquitoes where pre-exposed to PBO before being exposed to permethrin (PBO + permethrin), mortalities increased about 9-fold almost reaching 100 % for both males and females, showing that P450 monooxygenases are playing a major role in pyrethroid resistance at Kpome. This also further supports that pyrethroid resistance is mainly driven by metabolic resistance in field populations of An. funestus as observed throughout Africa with little evidence of knockdown resistance [19, 21]. In this study, because of the low number of mosquitoes collected, we did not carry out the biochemical analysis to determine levels/activities of known detoxifying enzyme systems including monooxygenases, non-specific esterases and GSTs as well as to detect the presence of an altered acetylcholinesterase associated with carbamate resistance. However, biochemical analysis performed on wild An. funestus from Ghana suggested that DDT and pyrethroid resistance may be metabolically mediated, although there were no clear correlations between enzyme levels/activities and insecticide resistance across families. Furthermore, an altered acetylcholinesterase conferring carbamate resistance was evidenced . The recovery of susceptibility to pyrethroids after PBO exposure is in line with the absence of any kdr mutation in Anopheles funestus as previously reported in the Pahou population . The very high frequency of the 119 F-GSTe2 resistant allele in Kpome, is similar to that observed in Pahou in 2011. Both data from Pahou and Kpome correlate the observed high level of DDT resistance in both localities and suggest that this mutation might play an important role in DDT resistance in West Africa. The results observed in Kpome are also in line with the detection of high frequency of this resistance allele in other DDT-resistant populations in West Africa notably in Cameroon (48.2 %), Ghana (44.2 %) and Burkina Faso (25 %), in accordance with the previously reported prevalence of DDT resistance in these countries . The resistant 119 F allele has also been detected in the eastern African countries of Uganda (20.4 %) and Kenya (7.8 %) but with lower frequencies, reflecting the moderate level of DDT resistance that was previously reported [8, 9]. However the above observations are in contrast with the complete absence of the 119 F in southern African populations despite the recent detection of DDT resistance in that region suggesting the presence of different mechanisms of DDT resistance in An. funestus [2, 7]. The observed high frequency of 119 F in both DDT susceptible and resistant phenotypes shows that this mutation is near fixation at Kpome; therefore, bioassays with DDT cannot be used for discriminating mutant individuals in this population. Such observations have also been documented with the kdr mutation in Anopheles gambiae in areas of high pyrethroid resistance and where this mutation (kdr) is almost fixed. We noticed in such areas that some dead mosquitoes post-pyrethroid bio-assays were carrying the kdr mutation [63, 64]. The 119 F resistant allele in Kpome was probably selected during large use of DDT for agriculture (i.e. against pests) [41, 42], and/or public health (against vectors) . Another possibility is that the 119 F allele could also be selected not by DDT but by pyrethroid-based interventions as it was shown that GSTe2 can metabolize permethrin .
This study has provided relevant information on the Plasmodium infection rate in An. funestus (s.s.) in the southern Benin. The study also highlighted the multiple resistance profile of An. funestus (s.s.) to insecticides used in public health such as permethrin, deltamethrin, DDT and bendiocarb. This multiple resistance profile previously observed in the Pahou population and now in the Kpome population highlights the need for further studies to assess the extent and geographical distribution of these resistances in An. funestus populations throughout Benin and West Africa. At the same time, these results call for the development of a more comprehensive analysis of the resistance mechanisms involved in order to improve the implementation and management of future control programs against this species and other Anopheles mosquitoes locally involved in malaria transmission.
DDT, dichlorodiphenyltrichloroethane; EST, esterase; WHO, World Health Organization
We thank Kpome community for their cooperation during field work. We also appreciate Claude Gande and Seun M. Atoyebi for their technical assistance and collaboration.
This work is supported by the Wellcome Trust (grant 099864/Z/12/Z) awarded to Dr. Rousseau Djouaka.
Availability of data and materials
All data generated and analysed during this study are included in the published article and its Additional files.
RD and CW carried out the study design, coordinated laboratory work, analyzed data and wrote the paper. RA, ID and RA conducted field sampling and contributed in laboratory analysis of sampled mosquitoes. JMR, GT and HI conducted laboratory analysis including advanced molecular work and contributed in the preparation of the manuscript. AY, KM, VM and MT advised on the study design and helped finalizing the manuscript. All authors read and approved the final version of the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
The request for ethical approval was not applicable. However, we obtained verbal consents of household heads before using aspirators for collection of mosquitoes in each room.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- NMCP. Rapport annuel d'activité. Benin: NMCP; 2014.Google Scholar
- Cuamba N, Helen I, Andrew S, Wondji CS. High level of pyrethroid resistance in an Anopheles funestus population of the Chokwe District in Mozambique. PLoS One. 2010;5(6):e11010.View ArticlePubMedPubMed CentralGoogle Scholar
- Wondji SS, Zainab T, Helen I, Rousseau D, Morgan JC. Identification and distribution of a GABA receptor mutation conferring dieldrin resistance in the malaria vector Anopheles funestus in Africa. Insect Biochem Mol Biol. 2011;41:484–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Djouaka R, Irving H, Tukur Z, CS W. Exploring mechanisms of multiple insecticide resistance in a population of the malaria vector Anopheles funestus in Benin. PLoS One. 2011;6(11):e27760.View ArticlePubMedPubMed CentralGoogle Scholar
- Riveron JM, Yunta C, Ibrahim SS, Djouaka R, Irving H, Menze BD, et al. A single mutation in the GSTe2 gene allows tracking of metabolically based insecticide resistance in a major malaria vector. Genome Biol. 2014;15(2):R27.View ArticlePubMedPubMed CentralGoogle Scholar
- Wondji CS, Kleinschmidt I, Mzilahowa T, Irving H, Ndula M, Rehman A, et al. Impact of pyrethroid resistance on operational malaria control in Malawi. Proc Natl Acad Sci USA. 2012;109(47):19063–70.View ArticlePubMedPubMed CentralGoogle Scholar
- Riveron JM, Menze BD, Barnes KG, Irving H, Ibrahim SS, Weedall GD, et al. Rise of multiple insecticide resistance in Anopheles funestus in Malawi: a major concern for malaria vector control. Mar J. 2015;14:334.Google Scholar
- Morgan JC, Irving H, Okedi LM, Steven A, Wondji CS. Pyrethroid resistance in an Anopheles funestus population from Uganda. PLoS One. 2010;5(7):e11872.View ArticlePubMedPubMed CentralGoogle Scholar
- Mulamba C, Riveron JM, Ibrahim SS, Irving H, Barnes KG, Mukwaya LG, et al. Widespread pyrethroid and DDT resistance in the major malaria vector Anopheles funestus in East Africa is driven by metabolic resistance mechanisms. Plos One. 2014;9(10):e110058.View ArticlePubMedPubMed CentralGoogle Scholar
- Okoye PN, Brooke BD, Koekemoer LL, Hunt RH, Coetzee M. Characterisation of DDT, pyrethroid and carbamate resistance in Anopheles funestus from Obuasi, Ghana. Trans R Soc Trop Med Hyg. 2008;102(6):591–8.View ArticlePubMedGoogle Scholar
- Djouaka RF, Bakare AA, Coulibaly ON, Akogbeto MC, Ranson H, Hemingway J, Strode C. Expression of the cytochrome P450s, CYP6P3 and CYP6M2 are significantly elevated in multiple pyrethroid resistant populations of Anopheles gambiae s.s. from Southern Benin and Nigeria. BMC Genomics. 2008;9:538.View ArticlePubMedPubMed CentralGoogle Scholar
- Yadouleton AW, Padonou G, Asidi A, Moiroux N, Banganna S, Corbel V, et al. Insecticide resistance status in Anopheles gambiae in southern Benin. Malar J. 2010;9(1):83.View ArticlePubMedPubMed CentralGoogle Scholar
- Djegbe I, Boussari O, Sidick A, Martin T, Ranson H, Chandre F, et al. 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. 2011;10:261.View ArticlePubMedPubMed CentralGoogle Scholar
- Virgile G, Kefilath B, Renaud G, Rodrigue A, Fredéric O-A, Roseric A, et al. Malaria vectors resistance to insecticides in Benin: current trends and mechanisms involved. Parasit Vectors. 2015;8:223.View ArticleGoogle Scholar
- Coetzee M, Hunt RH, Wilkerson R, Della Torre A, Coulibaly MB, Besansky NJ. Anopheles coluzzii and Anopheles amharicus, new members of the Anopheles gambiae complex. Zootaxa. 2013;3619(2):246.PubMedGoogle Scholar
- Coetzee M, Fontenille D. Advances in the study of Anopheles funestus, a major vector of malaria in Africa. Insect Biochem Mol Biol. 2004;34(7):599–605.View ArticlePubMedGoogle Scholar
- Hunt R, Edwardes M, Coetzee M. Pyrethroid resistance in southern African Anopheles funestus extends to Likoma Island in Lake Malawi. Parasit Vectors. 2010;3:122.View ArticlePubMedPubMed CentralGoogle Scholar
- Hargreaves K, Koekemoer LL, Brooke BD, Hunt RH, Mthembu J, Coetzee M. Anopheles funestus resistant to pyrethroid insecticides in South Africa. Med Vet Entomol. 2000;14(2):181–9.View ArticlePubMedGoogle Scholar
- Jacob M, Riveron HI, Miranda N, Barnes KG, Ibrahim SS, Paine MJI, Charles SW. Directionally selected cytochrome P450 alleles are driving the spread of pyrethroid resistance in the major malaria vector Anopheles funestus. Proc Natl Acad Sci USA. 2013;110(1):252–7.View ArticleGoogle Scholar
- Wondji CS, Irving H, Morgan J, Lobo NF, Collins FH, Hunt RH, et al. Two duplicated P450 genes are associated with pyrethroid resistance in Anopheles funestus, a major malaria vector. Genome Res. 2009;19(3):452–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Riveron JM, Ibrahim SS, Chanda E, Mzilahowa T, Cuamba N, Irving H, et al. The highly polymorphic CYP6M7 cytochrome P450 gene partners with the directionally selected CYP6P9a and CYP6P9b genes to expand the pyrethroid resistance front in the malaria vector Anopheles funestus in Africa. BMC Genomics. 2014;15:817.View ArticlePubMedPubMed CentralGoogle Scholar
- Edi CV, Djogbénou L, Jenkins AM, Regna K, Muskavitch MA, Poupardin R, et al. CYP6 P450 enzymes and ACE-1 duplication produce extreme and multiple insecticide resistance in the malaria mosquito Anopheles gambiae. Plos Genetics. 2014;10(3):e1004236.View ArticlePubMedPubMed CentralGoogle Scholar
- Sulaiman SI, Miranda N, Jacob MR, Helen I, Charles SW. The P450 CYP6Z1 confers carbamate/pyrethroid cross-resistance in a major African malaria vector beside a novel carbamate-insensitive N485I acetylcholinesterase-1 mutation. Mol Ecol. 2016.Google Scholar
- MT G. A supplement to the Anophelinae of Africa south of the Sahara (Afrotropical region). Johannesburg: South African Institute for Medical Research; 1987. p. 143.Google Scholar
- Koekemoer LL, Kamau L, Hunt RH, Coetzee M. A cocktail polymerase chain reaction assay to identify members of the Anopheles funestus (Diptera: Culicidae) group. Am J Trop Med Hyg. 2002;66(6):804–11.PubMedGoogle Scholar
- Livak KJ. Organization and mapping of a sequence on the Drosophila melanogaster X and Y chromosomes that is transcribed during spermatogenesis. Genetics. 1984;107(4):611–34.PubMedPubMed CentralGoogle Scholar
- Santolamazza F, Calzetta M, Etang J, Barrese E, Dia I, Caccone A, et al. Distribution of knock-down resistance mutations in Anopheles gambiae molecular forms in west and west-central Africa. Malar J. 2008;7:74.View ArticlePubMedPubMed CentralGoogle Scholar
- Bass CDN, Blagborough AM, Vontas J, Sinden RE, Williamson MS, Field LM. PCR-based detection of Plasmodium in Anopheles mosquitoes: a comparison of a new high-throughput assay with existing methods. Malar J. 2008;7:177.View ArticlePubMedPubMed CentralGoogle Scholar
- Snounou G, Viriyakosol S, Zhu XP, Jarra W, Pinheiro L, do Rosario VE, et al. High sensitivity of detection of human malaria parasites by the use of nested polymerase chain reaction. Mol Biochem Parasitol. 1993;61(2):315–20.PubMedGoogle Scholar
- WHO. Test procedures for insecticide resistance monitoring in malaria vector, bio-efficacy and persistence of insecticides on treated surfaces. Geneva: WHO/CDS/CPC/MAL/; 2013.Google Scholar
- Mitchell SN, Rigden DJ, Dowd AJ, Lu F, Wilding C, Weetman D, et al. Metabolic and target-site mechanisms combine to confer strong DDT resistance in Anopheles gambiae. PLoS One. 2014;9(3):e92662.View ArticlePubMedPubMed CentralGoogle Scholar
- Brooke BD, Kloke G, Hunt RH, Koekemoer LL, Temu EA, Taylor ME, et al. Bioassay and biochemical analyses of insecticide resistance in southern African Anopheles funestus (Diptera: Culicidae). Bull Entomol Res. 2001;91:265–72.View ArticlePubMedGoogle Scholar
- Meillon D. On Anopheles funestus and its allies in the Transvaa. Ann Trop Med Parasitol. 1933;27(1):83–97.View ArticleGoogle Scholar
- Swellengrebel NH, Annecke S, De Meillon B. Malaria investigations in some parts of the Transvaal and Zululand. South African Institute for Medical Research. 1931;245–274.Google Scholar
- Costantini C, Sagnon N, Ilboudo-Sanogo E, Coluzzi M, Boccolini D. Chromosomal and bionomic heterogeneities suggest incipient speciation in Anopheles funestus from Burkina Faso. Parassitologia. 1999;41(4):595–611.PubMedGoogle Scholar
- Dabire KR, Baldet T, Diabate A, Dia I, Costantini C, Cohuet A, et al. Anopheles funestus (Diptera: Culicidae) in a humid savannah area of western Burkina Faso: bionomics, insecticide resistance status, and role in malaria transmission. J Med Entomol. 2007;44(6):990–7.View ArticlePubMedGoogle Scholar
- Sandeu MM, Moussiliou A, Moiroux N, Padonou GG, Massougbodji A, Corbel V, Ndam NT. Optimized pan-species and speciation duplex real-time PCR assays for Plasmodium parasites detection in malaria vectors. PLoS One. 2012.Google Scholar
- Dadzie SK, Brenyah R, Appawu MA. Role of species composition in malaria transmission by the Anopheles funestus group (Diptera: Culicidae) in Ghana. J Vector Ecol. 2012;38(1):105–10.View ArticleGoogle Scholar
- Awolola TS, Oyewole IO, Koekemoer LL, Coetzee M. Identification of three members of the Anopheles funestus (Diptera: Culicidae) group and their role in malaria transmission in two ecological zones in Nigeria. Trans R Soc Trop Med Hyg. 2005;99(7):525–31.View ArticlePubMedGoogle Scholar
- Moiroux N, Boussari O, Djenontin A, Damien G, Cottrell G, Henry MC, et al. Dry season determinants of malaria disease and net use in Benin, West Africa. PLoS One. 2012;7(1), e30558.View ArticlePubMedPubMed CentralGoogle Scholar
- Akogbeto MC, Djouaka R, Noukpo H. Use of agricultural insecticides in Benin. Bull Soc Pathol Exot. 2005;98(5):400–5.PubMedGoogle Scholar
- Yadouleton A, Martin T, Padonou G, Chandre F, Asidi A, Djogbenou L, et al. Cotton pest management practices and the selection of pyrethroid resistance in Anopheles gambiae population in northern Benin. Parasites & vectors. 2011;4:60.View ArticleGoogle Scholar
- Djouaka RF, Bakare AA, Bankole HS, Doannio JM, Coulibaly ON, Kossou H, et al. Does the spillage of petroleum products in Anopheles breeding sites have an impact on the pyrethroid resistance? Malar J. 2007;6:159.View ArticlePubMedPubMed CentralGoogle Scholar
- Djogbenou L, Pasteur N, Akogbeto M, Weill M, Chandre F. Insecticide resistance in the Anopheles gambiae complex in Benin: a nationwide survey. Med Vet Entomol. 2011;25(3):256–67.View ArticlePubMedGoogle Scholar
- Antonio-Nkondjio C, Tene Fossog B, Kopya E, Poumachu Y, Menze Djantio B, Ndo C, et al. Rapid evolution of pyrethroid resistance prevalence in Anopheles gambiae populations from the cities of Douala and Yaoundé (Cameroon). Malar J. 2015;14:155.View ArticlePubMedPubMed CentralGoogle Scholar
- Kwiatkowska RMPN, Poupardin R, Irving H, Dabire RK, Mitchell S, Jones CM, et al. Dissecting the mechanisms responsible for the multiple insecticide resistance phenotype in Anopheles gambiae s.s, M form, from Vallée du Kou, Burkina Faso. Gene. 2013;519(1):98–106.View ArticlePubMedPubMed CentralGoogle Scholar
- Dabiré KR, Diabaté A, Namountougou M, Djogbenou L, Wondji C, Chandre F, et al. Trends in insecticide resistance in natural populations of malaria vectors in Burkina Faso, West Africa: 10 Years’ Surveys. In: Perveen F, editor. Insecticides - Pest Engineering. Rijeka: InTech; 2012. p. 479–502. ISBNGoogle Scholar
- Mulamba C, Irving H, Riveron JM, Mukwaya LG, Birungi J, Wondji CS. Contrasting Plasmodium infection rates and insecticide susceptibility profiles between the sympatric sibling species Anopheles parensis and Anopheles funestus s.s: a potential challenge for malaria vector control in Uganda. Parasit Vectors. 2014;7:71.View ArticlePubMedPubMed CentralGoogle Scholar
- Djègbè I, Agossa RF, Jones MC, Rodolphe P, Cornelie S, Akogbéto M, et al. Molecular characterization of DDT resistance in Anopheles gambiae from Benin. Parasites & Vectors. 2014;7:409.View ArticleGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Tokponnon FT, Ogouyémi AH, Sissinto Y, Sovi A, Gnanguenon V, Cornélie S, et al. Impact of long-lasting, insecticidal nets on anaemia and prevalence of Plasmodium falciparum among children under five years in areas with highly resistant malaria vectors. Malar J. 2014;13:76.View ArticlePubMedPubMed CentralGoogle Scholar
- Nkya TPRLF, Akhouayri I, Mosha F, Magesa S, Kisinza W, David JP. Impact of agriculture on the selection of insecticide resistance in the malaria vector Anopheles gambiae: a multigenerational study in controlled conditions. Parasit Vectors. 2014;7(1):480.PubMedPubMed CentralGoogle Scholar
- Antonio-Nkondjio C, Fossog BT, Ndo C, Djantio BM, Togouet SZ, Awono-Ambene P, et al. Anopheles gambiae distribution and insecticide resistance in the cities of Douala and Yaoundé (Cameroon): influence of urban agriculture and pollution. Malar J. 2011;10:154.View ArticlePubMedPubMed CentralGoogle Scholar
- Czeher C, Labbo R, Arzika I, Duchemin JB. Evidence of increasing Leu-Phe knockdown resistance mutation in Anopheles gambiae from Niger following a nationwide long-lasting insecticide-treated nets implementation. Malar J. 2008;7:189.View ArticlePubMedPubMed CentralGoogle Scholar
- Rock Aïkpon FA, Razaki O, Olivier O, Nazaire A, Frédéric O-A, Martin A. Bendiocarb resistance in Anopheles gambiae s.l. populations from Atacora department in Benin, West Africa: a threat for malaria vector control. Parasit Vectors. 2013;6:192.View ArticlePubMedPubMed CentralGoogle Scholar
- Akogbéto MC, Padonou GG, Gbénou D, Irish S, Yadouleton A. Bendiocarb, a potential alternative against pyrethroid resistant Anopheles gambiae in Benin, West Africa. 2010.Google Scholar
- Aïkpon R, Sèzonlin M, Ossè R, Akogbéto M. Evidence of multiple mechanisms providing carbamate and organophosphate resistance in field An. gambiae population from Atacora in Benin. Parasit Vectors. 2014;20147:568.View ArticleGoogle Scholar
- Badolo A, Traoré A, Jones CM, Sanou A, Flood L, Guelbéogo WM, et al. Three years of insecticide resistance monitoring in Anopheles gambiae in Burkina Faso: resistance on the rise? Malar J. 2012;11(1):232.View ArticlePubMedPubMed CentralGoogle Scholar
- Rowland M, Boko P, Odjo A, Asidi A, Akogbeto M, N’Guessan R. A new long-lasting indoor residual formulation of the organophosphate insecticide pirimiphos methyl for prolonged control of pyrethroid-resistant mosquitoes: An experimental hut trial in Benin. Plos One. 2013.Google Scholar
- Oxborough RM, N’Guessan R, Jones R, Kitau J, Ngufor C, Malone D, et al. The activity of the pyrrole insecticide chlorfenapyr in mosquito bioassay: towards a more rational testing and screening of non-neurotoxic insecticides for malaria vector control. Malar J. 2015;14:124.View ArticlePubMedPubMed CentralGoogle Scholar
- Oxborough RM, N’Guessan R, Kitau J, Tungu PK, Malone D, Mosha FW, Rowland WM. A new class of insecticide for malaria vector control: evaluation of mosquito nets treated singly with indoxacarb (oxadiazine) or with a pyrethroid mixture against Anopheles gambiae and Culex quinquefasciatus. Malar J. 2015;14:353.View ArticlePubMedPubMed CentralGoogle Scholar
- Hamon J, Sales S, Venard P, Coz J, Brengues J. The presence in southwest Upper Volta of populations of Anopheles funestus Giles resistant to dieldrin. Med Trop. 1968;28(2):221–6.Google Scholar
- Dabire KR, Diabate A, Agostinho F, Alves F, Manga L, Faye O, Baldet T. Distribution of the members of Anopheles gambiae and pyrethroid knock-down resistance gene (kdr) in Guinea-Bissau, West Africa. Bull Soc Pathol Exot. 2008;101(2):119–23.PubMedGoogle Scholar
- Natacha P, Johnson M, Robert M, Reginald K, Robert K, Alexandra W, et al. High level of resistance in the mosquito Anopheles gambiae to pyrethroid insecticides and reduced susceptibility to bendiocarb in north-western Tanzania. Malar J. 2013;12:149.View ArticleGoogle Scholar