- Open Access
Evidence of a multiple insecticide resistance in the malaria vector Anopheles funestus in South West Nigeria
Malaria Journalvolume 15, Article number: 565 (2016)
Knowing the extent and spread of insecticide resistance in malaria vectors is vital to successfully manage insecticide resistance in Africa. This information in the main malaria vector, Anopheles funestus sensu stricto, is completely lacking in the most populous country in Africa, Nigeria. This study reports the insecticide susceptibility status and the molecular basis of resistance of An. funestus as well as its involvement in malaria transmission in Akaka-Remo, a farm settlement village in southwest Nigeria.
Plasmodium infection analysis using TaqMan protocol coupled with a nested PCR revealed an infection rate of 8% in An. funestus s.s. from Akaka-Remo. WHO susceptibility tests showed this species has developed multiple resistance to insecticides in the study area. Anopheles funestus s.s. population in Akaka-Remo is highly resistant to organochlorines: dieldrin (8%) and DDT (10%). Resistance was also observed against pyrethroids: permethrin (68%) and deltamethrin (87%), and the carbamate bendiocarb (84%). Mortality rate with DDT slightly increased (from 10 to 30%, n = 45) after PBO pre-exposure indicating that cytochrome P450s play little role in DDT resistance while high mortalities were recorded after PBO pre-exposure with permethrin (from 68 to 100%, n = 70) and dieldrin (from 8 to 100%, n = 48) suggesting the implication of P450s in the observed permethrin and dieldrin resistance. High frequencies of resistant allele, 119F in F0 (77%) and F1 (80% in resistant and 72% in susceptible) populations with an odd ratio of 1.56 (P = 0.1859) show that L119F-GSTe2 mutation is almost fixed in the population. Genotyping of the A296S-RDL mutation in both F0 and F1 samples shows an association with dieldrin resistance with an odd ratio of 81 (P < 0.0001) (allelic frequency (R) = 76% for F0; for F1, 90 and 10% were observed in resistant and susceptible populations, respectively) as this mutation is not yet fixed in the population.
The study reports multiple insecticide resistance in An. funestus from Akaka Remo. It is, therefore, necessary to pay more attention to this major malaria vector for effective malaria control in Nigeria.
Malaria remains the most severe infectious disease and a major public health challenge in Nigeria [1, 2]. It is the main cause of morbidity and mortality in this most populous Africa country, with 97% of the national population at risk: it is responsible for an estimated 300,000 deaths annually in Nigeria; and it contributes to an estimated 11% maternal mortality as well as 25% of infant mortality [1, 3]. Malaria transmission in Nigeria has been attributed mainly to Anopheles gambiae sensu stricto (s.s.) and Anopheles funestus s.s. [4–6] with consistent Plasmodium infection rates of 1.0–2.7% (An. funestus) and 3.0–8.1% (An. gambiae) previously reported in case studies in Ogun, Oyo and Lagos states [4, 5]. Although, there was also a high sporozoite infection rate of 25% reported in Lagos state . In Nigeria, malaria control relies hugely on the use of indoor residual spraying (IRS) and insecticide-treated nets (ITNs) [2, 3]. However, resistance against the main insecticides used in public health (pyrethroids, carbamates and organochlorines) in malaria vectors is threatening the effectiveness of these control tools. Anopheles gambiae s.s resistance to insecticides, notably against pyrethroids , DDT [9, 10] and bendiocarb , has been documented in Nigeria, however, little is known so far concerning the insecticide susceptibility of the other major malaria vector An. funestus s.s. in the country. Pyrethroid insecticide is the class of insecticide mainly used in Nigeria for both ITNs and IRS . Two types of pyrethroids are mainly used in Nigeria for insecticide nets treatment: permethrin (Type 1) and deltamethrin (Type 2). In recent years, An. funestus s.s. populations have increasingly been reported to be resistant to these insecticides in other African countries, such as Uganda in East Africa ; Mozambique, Zambia, Zimbabwe and Malawi in Southern Africa [13–19], Cameroon in Central Africa [20, 21] and some West African countries, including Benin , Ghana  and Burkina-Faso . Resistance patterns against these insecticides vary significantly across Africa. For example, An. funestus was resistant to pyrethroids and carbamate but fully susceptible to DDT and dieldrin in southern Africa [20, 25]. However, a recent study in Malawi showed that this mosquito species has now began to develop resistance against organochlorines (dieldrin and DDT) . Anopheles funestus is resistant to pyrethroids and DDT, but remains susceptible to carbamate in Uganda and western Kenya . High resistance profiles were recorded with dieldrin in Cameroon . In the neighbouring country of Benin, resistance firstly reported in 2011  from the coast (Pahou) and was recently shown to have extended to the inland as the Kpome population was shown to be resistant to all insecticide classes apart from organophosphates . It remains to be established whether these resistances are also present in Nigeria and if yes information on the resistance pattern will be useful for the malaria control programs especially on the suitable insecticides to use for the control of this species.
Metabolic resistance mechanisms have so far been implicated in insecticide-resistant An. funestus across Africa [12, 22, 27] with cytochrome P450 genes conferring pyrethroid resistance and also cross-resistance to carbamates in southern African  as previously reported also for An. gambiae . DDT resistance mechanisms in An. funestus on the other hand have been associated with an up-regulation of glutathione S-transferases notably GSTe2 coupled with a point mutation L119F . No L1014F-kdr mutation has been implicated in pyrethroids and DDT resistance [12, 22], and no association exists between G119S and F455W mutations of the Ace-1 gene and carbamate resistance in this mosquito species [12, 22, 25]. However, the recent discovery of a new Ace-1 mutation (N485I) associated with carbamate resistance in southern African An. funestus populations  coupled with the presence of the A296S-RDL mutation in the GABA receptor of An. funestus  are evidence that target-site resistance mechanism also play a role in insecticide resistance profiles recorded in this malaria vector.
In order to help malaria control programmes, to design evidence-based strategies to control An. funestus in Nigeria, and to manage potential existing resistance, this study aims to establish the insecticides susceptibility profile and investigate the molecular basis of resistance of this species population in Akaka Remo: a farm settlement in southwest Nigeria.
No ethical permit was required for this study. However, there was a focus group discussion with the community and household heads where verbal consent was obtained for mosquito collections in the community after the study aims and objectives were explained.
Study site and mosquito collection
Study site description
Akaka-Remo (6°57′N, 3°43′E) is a rural locality in Remo-North local government of Ogun state in the Southwest of Nigeria, a region of about 71.4 km from Lagos and about 215 km from Pahou in Benin where resistance has previously been reported. This locality is surrounded with a permanent medium-size slow moving stream, called Erititi stream, that leads to the popular river Ona in Ibadan (Oyo state) with vegetation such as bananas, vegetables, maize, shrubs, trees and crops bordering the water bodies at almost all the locations, which serves as suitable breeding sites for An. funestus. The inhabited area of Akaka-Remo is about 0.25 square kilometres and its habitants are mainly the Yorubas and a small community of the Eguns. The main commercial activity is agriculture, which has attracted the use of pesticides (insecticides and herbicides) in this locality. Houses here are mainly made of mud and very few are made of cement, and are constructed at an average of 5 m away from one other. Most houses have either detached/destroyed or no ceilings. The selection criteria described above were mainly entomological as the main target for this research was to characterise populations of malaria vectors in this part of Nigeria.
Adult female Anopheles mosquitoes resting indoor were collected in thirty rooms, with the use of electric aspirators and torches between 06.00 a.m. and 10.00 a.m. from October, 2014 end of rainy season to April, 2015 beginning of rainy season except in January, 2015 due to intense harmattan (a short period of a very dry and dusty wind observed between the end of November and early March in West Africa). The 30 rooms were selected in a way to cover the various micro-ecologies found at Akaka Remo. The room was defined as a demarcated area in the house where inhabitants do sleep. Blood-fed mosquitoes collected were kept into cups until fully gravid before being subjected to the forced-egg laying technique  at a temperature of 25–28 °C with a relative humidity of 80% in the insectary. Hatched eggs were pooled and reared together in a mineral water, which was replaced every two days to reduce mortality as resulting larvae were daily fed with Tetramin™ baby fish. During these periods, a good number of eggs were also sent via courier to the Liverpool School of Tropical Medicine (LSTM) for rearing into F1 and for subsequent experiments.
Seasonal determination of mosquito densities per room
Mosquito density per room was estimated during four (4) annual climatic seasons: rainy season, transition from rainy to dry season, dry season and transition from dry to rainy season. The total number of An. funestus collected for each season were pooled and counted to estimate the seasonal number of An. funestus per room. The estimated density of An. funestus was now obtained by dividing the number of mosquitoes collected during each season by the number of rooms surveyed during that same season. This estimation of the density was done for the 4 surveyed seasons of the year.
A total of 96 mosquito females that were morphologically identified as An. funestus group  and had oviposited, were identified to species level using the PCR cocktail for An. funestus group described by Koekemoer et al.  after the genomic DNA was extracted .
Estimation of Plasmodium infection in wild caught (F0) Anopheles funestus
Ninety-three (93) F0 adult female An. funestus were analysed for Plasmodium infection using the TaqMan assay as previously described . Briefly, a plate was run at one cycle of 95 °C for 10 min in the first segment and the second segment was 40 cycles at 92 °C for 15 s and 60 °C for 1 min. Two fluorophore-labelled specific TaqMan probes (Applied Biosystems, California, USA) were used: FAM to detect Plasmodium falciparum, while HEX was used to detect the combination of Plasmodium ovale, Plasmodium vivax and Plasmodium malariae. A negative control (water) and positive controls (known FAM and OVM) were also used. A nested PCR  was subsequently performed for the samples to validate the TaqMan analysis.
Insecticide susceptibility tests
2–5 day old F1 adult female and male mosquitoes pooled from different F0 mosquitoes were used for this test according to the WHO . 20–25 mosquitoes per tube with at least 4 replicates were exposed to insecticide-impregnated or control papers for 1 h before transferring into clean holding tube with 10% sugar solution where mortality was determined after 24 h post insecticide exposure . Six insecticides belonging to the four classes of insecticides used for malaria vector control were tested: the pyrethroids permethrin (0.75%) and deltamethrin (0.05%), the organochlorines DDT (4%) and dieldrin (4%), the carbamate bendiocarb (0.1%) and the organophosphate malathion (5%).
PBO synergist tests
Due to the level of resistance observed against DDT, dieldrin and permethrin and because of previous studies showing strong involvement of P450 genes in pyrethroids resistance as well as its potential contribution to DDT resistance, 2–5 days old F1 adult mosquitoes were pre-exposed to 4% piperonyl butoxide (PBO) paper for 1 h and immediately exposed to 0.75% permethrin, 4% DDT and 4% dieldrin for 1 h. Although, there is no previous data linking P450 families to dieldrin resistance in this mosquito population but with the high resistance observed in this study, it became necessary to assess the potential effect of oxidase in diedrin resistance. Mortalities were later assessed after exposure; synergized group was compared to the un-synergized group after 24 h post-exposure. This comparison was used to evaluate the potential role of cytochrome P450 genes in the observed resistance. Two controls were used during this experiment: control 1 was constituted of mosquitoes exposed to papers neither with insecticides nor with PBO while control 2 was constituted of mosquitoes exposed to papers treated with PBO only.
Genotyping of resistance markers L119F-GSTe2 and A296S-RDL in females of An. funestus from Akaka-Remo
TaqMn assay  was used to genotype L119F-GSTe2 as a potential DDT resistance marker in Akaka-Remo, which was recently shown to confer DDT resistance in Benin  and also used to screen for A296S-RDL mutation known to confer dieldrin resistance . F0 and F1 (alive and dead) samples were used for both genotyping. Two fluorophore-labelled specific TaqMan probes were used: FAM to detect the homozygous resistant genotype, HEX to detect the homozygous susceptible genotype while both FAM and HEX were used to detect the heterozygous genotype. A negative control (water) and positive controls (known FAM, HEX and both) were also used in a 10 µl volume that also contains the SensiMix (Applied Biosystems, California, USA). The endpoint fluorescence was evaluated using the Agilent MXPro software and the relationship between the frequency of the resistant alleles and the insecticides (DDT and Dieldrin) resistance phenotypes were assessed.
Resistance status of mosquito classified as recommended by WHO  are as follows:
Susceptible mosquito population = Mortality >98%.
Suspected resistance in mosquito population = Mortality ranging from 90 to 98%.
Resistant mosquito population = Mortality <90%.
Chi square (using R software) was used to test for significant difference in percentage mortalities between female and male mosquito populations used for WHO susceptibility test and the distribution of the genotype frequencies (F1) between the resistant and susceptible mosquito samples. Had2know online statistical software  was used to test for significant difference between the observed genotypic frequencies (F0) and to confirm if the observed genotypic frequencies are according to Hardy–Weinberg equilibrium. Excel was used to compute percentage mortalities and standard errors while VassarStats online statistical software  was used to generate odd and risk ratios, and confidence levels of the frequency data.
Molecular (PCR) analysis of ninety-six (96) morphologically identified female An. funestus sensu lato collected from Akaka-Remo between October, 2014 (late rainy season) and April, 2015 (early rainy season) revealed that they all belong to An. funestus s.s. Anopheles funestus is the most abundant mosquito species (84%; n = 315 from a total of 376 mosquitoes collected) amongst other mosquito species and much more abundant (92%; An. funestus = 315 and 8%; An. gambiae = 26) than An. gambiae when compared within the Anopheles group. Figure 1 shows the seasonal variation of An. funestus s.s. per room at Akaka-Remo. The seasonal density of An. funestus per room (m/r) out of thirty (30) rooms aspirated are as follows: 0.03 m/r for rainy season, 1.8 m/r during transition from rainy to dry, 4 m/r during dry and 4.7 m/r during the transition from dry to rainy season. Other mosquito species also collected during these periods include An. gambiae spp. (7%; n = 26), Culex spp. (6%; n = 21), Mansonia spp. (2%; n = 9) and Aedes spp. (1%; n = 5).
Plasmodium infection rates
Seven, 7 (8%) out of ninety-three (93) wild-caught (F0) An. funestus s.s. analysed were positive for Plasmodium infection (Table 1). Six (7%) mosquitoes of which were infected with P. falciparum, while a mixed infection of P. ovale, P. vivax and P. malariae was found in only one mosquito (1%). The nested PCR analysis showed the presence of P. falciparum in 4 (4%) mosquitoes.
WHO susceptibility tests
A total of 96 F0 An. funestus oviposited out of the 315 samples collected on the field producing 1269 F1 adults (679 females and 590 males), which were all exposed to six different insecticides (Fig. 2). The highest level of resistance was recorded with organochlorines. Dieldrin exposure resulted into mortalities of 8% ± 3.24 (females) and 22% ± 1.73 (males). Likewise, DDT exposure produced mortalities of 10% ± 2.66 in females and 17% ± 2.45 in male populations. Resistance was also observed against both type I and II pyrethroids (without and with cyano group), with a mortality of 68% ± 5.64 in females (85% ± 3.15 for males) for permethrin (type I) and a mortality of 87% ± 10.96 (94% ± 3.98 for males) for deltamethrin (type II). In addition, bendiocarb (carbamate) resistance was also observed with mortalities of 84% ± 5.67 in females and 90% ± 2.36 for males. In contrast, a full susceptibility of 100% mortality was recorded in both females and males populations exposed to the organophosphate malathion. Overall, there was no significant difference (χ 2 = 7.73, df = 5, P = 0.172) in the percentage mortalities between the exposed females and males mosquitoes.
Synergist tests with PBO
There was a recovery of susceptibility to permethrin as mortality rose from 68 to 100% (n = 70) when permethrin was combined with the P450 inhibitor, PBO (Fig. 3). This suggests a likely significant role of cytochrome P450s in the pyrethroid resistance. Similarly, 100% mortality was recorded when PBO was combined with dieldrin. This unexpected recovery of susceptibility from 8 to 100% (n = 48) also implicates oxidases in dieldrin resistance. However, the combination of DDT with PBO only showed a slow increase in mortality from 10 to 30% (n = 45) suggesting a limited implication of P450 s in DDT resistance. No mortality was observed in the control mosquitoes exposed to control paper with no insecticide or only to PBO.
Genotyping and allelic distribution of L119F-GSTe2 mutation in the An. funestus population from Akaka-Remo
The L119F-GSTe2 mutation was detected in 94% of the F0 mosquitoes (n = 88) that were genotyped (Fig. 4a). Over half (52) of the total mosquitoes analysed were homozygous resistant RR, 31 were heterozygous RS while just 5 were homozygous susceptible, SS with allelic frequencies of R = 77% and S = 23%. Similarly, when the F1 generations (25 resistant and 25 susceptible after bioassays with DDT) were screened for L119F-GSTe2 mutation, a genotypic frequency of 64% RR, 32% RS, and 4% SS and 48% RR, 48% RS and 4% SS were produced in the resistant and susceptible populations respectively. These resulted into allelic frequencies (119F) of 80% in the resistant and 72% in the susceptible populations (Fig. 4b). The observed genotypic frequency was shown to be at Hardy–Weinberg equilibrium (P = 0.8935). However, there was no significant difference (χ 2 = 1.37, df = 2, P = 0.5037) in the frequency of L119F-GSTe2 mutation between the susceptible and resistant samples and consequently the correlation was also not significant (OR = 1.56; P = 0.1859).
Genotyping and allelic distribution of A296S-RDL mutation in An. funestus s.s. from Akaka-Remo
The A296S-RDL mutation was high in the F0 population (92 individuals genotyped) with homozygote RR claiming over half of the total mosquito analysed (50). Likewise, 40 samples were heterozygous RS and just 2 being homozygote susceptible SS (Fig. 5a). The high presence of A296S-RDL mutation in F0 mosquitoes correlates with the elevated phenotypic dieldrin resistance in this population (8 and 22% for females and males, respectively). But when F1 mosquitoes generated after bioassays with dieldrin (15 alive and 5 dead) were genotyped for A296S-RDL mutation, there was a high presence of the mutation in resistant population (genotypic frequency of 80% RR and 20% RS) with a relatively low presence in susceptible (genotypic frequency of 10% RS) sample. These produced allelic frequencies (296S) of 90% and 10% in the resistant and susceptible populations respectively (Fig. 5b). The observed genotypic frequency was also shown to be at Hardy–Weinberg equilibrium (P = 0.0617). There was a significant difference (χ 2 = 16, df = 2, P = 0.00034) in the frequency of A296S-RDL mutation in the resistant population compared to the susceptible sample and consequently correlation was also significant (OR = 81; P < 0.0001).
Role of An. funestus in malaria transmission at Akaka-Remo
Anopheles funestus is the most abundant mosquito species (84%) recorded at Akaka-Remo during the sampling period. Other mosquito species identified include An. gambiae s.l., Culex spp, Aedes spp. and Mansonia spp.; these represent 26% of sampled mosquitoes. Among the malaria vectors identified in this locality, An. funestus density was over 10 times higher that An. gambiae (92 and 8% respectively). This present report is in contrast with the report of Oyewole et al.  in 2005, where An. funestus collected, n = 85 was nowhere near that of An. gambiae, n = 500. This study supports previous reports in Ogun , Lagos , Oyo and Kwara states  that An. funestus s.s. is a major malaria vector in Nigeria with a confirmation of 8% Plasmodium infection in this mosquito species. It therefore emphasizes the importance of this vector and the threat it could pose to malaria transmission in the Southwest of Nigeria. The 8% infection rate observed in Akaka-Remo is similar to high levels of infection rates recorded previously for An. funestus across the continent such as the 20%  and 50%  observed in Burkina Faso, the 13.6%  and 18%  observed in Benin and 12.5% in Ghana . Although, some of the variations between these rates could be down to the differences in the methods used (TaqMan, ELISA and Nested-PCR) and the consistent high levels support a high vectorial capacity of An. funestus exhibits across the continent.
Another common member of the An. funestus group, Anopheles rivulorum that was previously identified both indoor and outdoor at Akaka-Remo  was absent in this study potentially due to a change in resting preference of this species although more entomological studies are needed to explain this change. Previous reports on the other known malaria vector, An. gambiae in Nigeria have shown that infection rates range mostly between 2 and 8.1% [4–6, 43]. This research has shown the in-houses abundant presence of An. funestus at Akaka-Remo compared to An. gambiae during all climatic seasons. It has also revealed the consistent level of mosquito infections with Plasmodium species (8% infected mosquitoes) in this locality hence, highlighting the fact that An. funestus plays a significant role in malaria transmission in this community in southwest Nigeria.
Multiple insecticide resistance at Akaka-Remo
This study reports that An. funestus s.s. in Akaka-Remo has developed resistance to most common public health insecticides. Results from this research highlight the presence of multiple insecticide resistance in this malaria vector. Most studies on insecticide resistance in Nigeria have focused on An. gambiae with less interest on An. funestus as this vector was essentially thought to be susceptible to pyrethroids, the main insecticide used for malaria vector control. This study has shown that this An. funestus population has not only become resistant to pyrethroids but to a wide range of public health insecticides from other families. The multi-insecticide resistance pattern observed in Akaka-Remo is similar to what was previously reported in Benin [22, 26].
The level of DDT resistance observed from Akaka Remo is as high as the case in Benin. This present report is higher than what was observed in Uganda and Kenya (40–42% mortality) , and in Malawi (69.9% mortality)  among others. The high DDT resistance recorded in Benin (Pahou) and now in Nigeria (Akaka Remo), both in West Africa compare to a relatively lower resistance profiles in the East and the Southern Africa might be as a result of different genetic make-up of this species between regions of Africa. Resistance recorded with dieldrin is not only the highest of all the six insecticides tested but also the highest recorded in Africa until now. Resistance level (8% mortality) is higher than Burkina-Faso (30% mortality)  and Benin (93.3% mortality) populations . In Southern Africa, dieldrin susceptibility has been frequently observed until a recent report of resistance in Malawi (83.9% mortality) . Organochlorines resistance recorded in this vector at Akaka-Remo coupled with reports of DDT resistance in An. gambiae  will obviously further disapproves the re-introduction of this insecticide family as an alternative to pyrethroids for mosquito control in Nigeria. DDT and dieldrin resistance recorded in this study could be associated with the residual effect of the long historical usage of this insecticide family (organochlorine) in agriculture when this sector was a key source of income in Nigeria . The oil boom in the 1970s latter shifted the national attention from agriculture to the oil and gas sectors. This economic sector (oil production) was more lucrative than agriculture and became the main asset for Nigeria . One could also argue that resistance may be due to the poor attitudes and/or ignorance of farmers towards observing good farming practices when using this insecticide family to control pests [46, 47]. During mosquito collections, it was observed that some of the villagers even make use of these agrochemicals to control insects at home. Such ignorance could add-on to the high DDT and dieldrin resistance in this study area. It is therefore important to determine the extent to which An. funestus has developed resistance to this insecticide family by investigating its spread across different geographical regions in Nigeria.
In Nigeria, agrochemicals use are approved by the National Agency for Food and Drug Administration and Control, NAFDAC . It is possible that the misuse and/or over use of these chemicals by farmers could be fundamental for the multi-resistance selection in this locality. The indiscriminate use of agro-chemicals by farmers could have also generated high chemical residues and other environmental pollutants that are washed into the water bodies (mosquito breeding sites) generating several xenobiotics that exercise a resistance selection in mosquitoes at larval stage [49–52]. Similarly, spilled petroleum products found in several mosquito breeding sites in the south–western Nigeria  might have also contributed to resistance selection in An. funestus from Akaka Remo through cross resistance mechanisms. Both environmental factors (generation of chemical pesticides and spillage of petroleum products) are common in Nigeria and can certainly contribute to the local selection of the observed insecticide resistance profiles. Further assessments are still needed to clearly map out the actual factor(s) contributing to the multi-insecticide resistance in An. funestus population from Akaka-Remo.
Pyrethroid resistance on the other hand is high but not as the previous two insecticides (DDT and dieldrin). Susceptibility test with pyrethroids is important because the country depends on this insecticide family for malaria vector control . Mortalities recorded with permethrin (68%) and deltamethrin (87%) are similar to what was reported in Pahou (permethrin = 66.7% and deltamethrin = 88.60%)  but higher than Kpome (permethrin = 13.03% and deltamethrin = 46.49) . These reports show that pyrethroids resistance is increasing in the West Africa population of An. funestus. The pattern of pyrethroids resistance in Western Africa is different from East  and Southern Africa , where resistance to deltamethrin is higher than permethrin. Pyrethroids resistance recorded in An. funestus from Akaka-Remo is a great concern for malaria control programs and there is a risk that this mosquito species would have developed resistance to pyrethroids in different regions of Nigeria due to the current heavy use of this class of insecticide both in agriculture and public health all over the country. If this happens, it will constitute more ordeals for future malaria vector control interventions through the use of pyrethroid-based insecticides for both ITNs and IRS. These findings therefore suggest further studies to determine the extent of pyrethroid resistance in An. funestus populations in Nigeria.
Bendiocarb resistance (84% mortality) in this species population is also a concern because carbamate-based insecticide interventions were recently introduced as an alternative to pyrethroid-based in West Africa. Bendiocarb resistance was also reported in Pahou (65% mortality) in 2011. Bendiocarb resistance is higher in Southern Africa: Zambia, Zimbabwe and Mozambique [16, 25, 54] but lower in East Africa  compared to West Africa. The current level of bendiocarb resistance raises an alarm as it might affect the success of the recently introduced bendiocarb-based IRS in West Africa. Hence, it will be important to have more information on the potential spread of resistance across different regions and some underlying factors that might be responsible. This information will guide the malaria control programs to improve its subsequent release of bendiocarb-based IRS. The organophosphate malathion is really proving to be the most reliable insecticide considering similar records of full susceptibility in An. funestus population all over Africa. This insecticide could, therefore, be used as an alternative insecticide to manage resistance.
Underlying mechanisms of the observed multiple insecticide resistance patterns at Akaka-Remo
The proven absence of the kdr mutation in An. funestus populations from Benin (a neighbouring country)  and other regions in Africa  coupled with the high mortality observed when permethrin was combined with PBO in this study support the fact that pyrethroid resistance is still driven by metabolic resistance mechanisms. The use of the synergist PBO revealed the role of oxidase, notably cytochrome P450s, in pyrethroid resistance of this mosquito population similar to what have been reported so far in Africa [12, 26]. However, the glutathione S-transferase gene, GSTe2 was shown to play a higher role in DDT resistance in Benin through over-transcription and also the selection of the resistant allele L119F . The high frequency of L119F in Akaka-Remo (77% in F0) associated with the high DDT resistance level support a significant role for the L119F-GSTe2 mutation in the DDT resistance in this An. funestus population. However, this high frequency of L119F in this location is probably the reason why a lack of correlation was observed when comparing resistant and susceptible samples as observed regularly for kdr mutations such as L1014F in An. gambiae in situation when the 1014F resistant allele is nearly fixed . The L119F-GSTe2 mutation has also been detected in other DDT resistant populations such as Ghana (44.2%) and Burkina-Faso (25%) in West Africa, and Cameroon (48.2%) in Central Africa [20, 23] as well as Uganda (20.4%) and Kenya (7.8%) in East Africa . In Southern Africa (Malawi) however, L119F allele is absent despite the recent detection of DDT resistance in this region suggesting a different DDT resistance mechanism in this An. funestus population .
Dieldrin resistance on the other hand showed a strong association with oxidase. This is unexpected because dieldrin resistance has always been linked with target-site insensitivity . Also, mutation detected on GABA receptor (A296S-RDL) gene in the parent and first filial generation of this An. funestus population indicates that this mosquito species adopts more than one mechanism for dieldrin resistance. More screening of dieldrin susceptibility should be done and further studies should be conducted to determine the geographical distribution of dieldrin resistance in An. funestus from Nigeria.
This study reports the presence of multiple insecticide resistance in An. funestus population from Akaka-Remo in the southwestern Nigeria. Molecular analysis conducted in the course of this research have revealed that this An. funestus population have developed multiple resistance mechanisms to withstand lethal doses of insecticides used in public health. The consistent implication of An. funestus in malaria transmission at Akaka-Remo was also established in this study. Nevertheless, further studies are needed to determine the spread of insecticide resistance and to conduct more investigations on underlying mechanisms of insecticides resistance for improved malaria control strategies in Nigeria.
National Agency for Food and Drug Administration and Control
mosquito per room
Liverpool School of Tropical Medicine
resistance to dieldrin
polymerase chain reaction
World Health Organization
CDC. Malaria Operational Plan FY. President’s Initiative in Nigeria. 2015.
World Health Organization. World malaria report. Geneva: World Health Organization; 2016.
CDC. Malaria Operational Plan FY. President’s Initiative in Nigeria. 2016.
Oyewole IO, Ibidapo CA, Oduola AO, Obansa JB, Awolola TS. Anthropophilic mosquitoes and malaria transmission in a tropical rain forest area of Nigeria. J Life Phys Sci. 2005;2:6–10.
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:525–31.
Oyewole IO, Awolola TS. Impact of urbanisation on bionomics and distribution of malaria vectors in Lagos, southwestern Nigeria. J Vector Borne Dis. 2006;43:173–8.
Okwa O, Rasheed A, Adeyemi A, Omoyeni M, Oni L, Fayemi A, et al. Anopheles species abundances, composition and vectoral competence in six areas of Lagos: Nigeria. J Cell Anim Biol. 2007;1:19–23.
Djouaka RF, Bakare AA, Coulibaly ON, Akogbeto MC, Ranson H, Hemingway J, et al. 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 Genom. 2008;9:538.
Oduola AO, Olojede JB, Ashiegbu CO, Olufemi A, Otubanjo OA, Awolola TS. High level of DDT resistance in the malaria mosquito: Anopheles gambiae s.l. from rural, semi urban and urban communities in Nigeria. J Rural Trop Public Health. 2010;9:114–20.
Oduola AO, Idowu ET, Oyebola MK, Adeogun AA, Olojede JB, Otubanjo OA, et al. Evidence of carbamate resistance in urban populations of Anopheles gambiae s.s. mosquitoes resistant to DDT and deltamethrin insecticides in Lagos. Int Pest Control. 2012;54:206.
Ibrahim KT, Popoola KO, Adewuyi OR. Susceptibility of Anopheles gambiae sensu lato (Diptera:Culicidae) to permethrin, deltamethrin and bendiocarb in Ibadan City, Southwest Nigeria. Curr Res J Biol Sci. 2013;5:42–8.
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:e110058.
Amenya DA, Naguran R, Lo T-CM, Ranson H, Spillings BL, Wood OR, et al. Over expression of a cytochrome P450 (CYP6P9) in a major African malaria vector, Anopheles funestus, resistant to pyrethroids. Insect Mol Biol. 2008;17:19–25.
Wondji CS, Coleman M, Kleinschmidt I, Mzilahowa T, Irving H, Ndula M, et al. Impact of pyrethroid resistance on operational malaria control in Malawi. Proc Natl Acad Sci USA. 2012;109:19063–70.
Vezenegho SB, Chiphwanya J, Hunt RH, Coetzee M, Bass C, Koekemoer LL. Characterization of the Anopheles funestus group, including Anopheles funestus-like, from Northern Malawi. Trans R Soc Trop Med Hyg. 2013;107:753–62.
Choi KS, Christian R, Nardini L, Wood OR, Agubuzo E, Muleba M, et al. Insecticide resistance and role in malaria transmission of Anopheles funestus populations from Zambia and Zimbabwe. Parasit Vectors. 2014;7:464.
Riveron JM, Chiumia M, Menze BD, Barnes KG, Irving H, Ibrahim SS, et al. Rise of multiple insecticide resistance in Anopheles funestus in Malawi: a major concern for malaria vector control. Malar J. 2015;14:344.
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:181–9.
Casimiro S, Coleman M, Mohloai P, Hemingway J, Sharp B. Insecticide resistance in Anopheles funestus (Diptera: Culicidae) from Mozambique. J Med Entomol. 2006;43:267–75.
Wondji CS, Dabire RK, Tukur Z, Irving H, Djouaka R, 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.
Menze BD, Riveron JM, Ibrahim SS, Irving H, Antonio-Nkondjio C, Awono-Ambene PH, et al. Multiple insecticide resistance in the malaria vector Anopheles funestus from Northern Cameroon is mediated by metabolic resistance alongside potential target site insensitivity mutations. PLoS ONE. 2016;11:e163261.
Djouaka R, Irving H, Tukur Z, Wondji CS. Exploring mechanisms of multiple insecticide resistance in a population of the malaria vector Anopheles funestus in Benin. PLoS ONE. 2011;6:e27760.
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:591–8.
Dabiré KR, Baldet T, Diabaté 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:990–7.
Cuamba N, Morgan JC, Irving H, Steven A, Wondji CS. High level of pyrethroid resistance in an Anopheles funestus population of the chokwe district in Mozambique. PLoS ONE. 2010;5:e11010.
Djouaka R, Riveron JM, Yessoufou A, Tchigossou G, Akoton R, Irving H, et al. Multiple insecticide resistance in an infected population of the malaria vector Anopheles funestus in Benin. Parasit Vectors. 2016;9:453.
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:27.
Ibrahim SS, Ndula M, Riveron JM, Irving H, Wondji CS. 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;25:3436–52.
Edi CV, Djogbénou L, Jenkins AM, Regna K, Muskavitch MAT, Poupardin R, et al. CYP6 P450 enzymes and ACE-1 duplication produce extreme and multiple insecticide resistance in the malaria mosquito Anopheles gambiae. PLoS Genet. 2014;10:e1004236.
Morgan JC, Irving H, Okedi LM, Steven A, Wondji CS. Pyrethroid resistance in an Anopheles funestus population from Uganda. PLoS ONE. 2010;5:e11872.
Gillies MT, Coetzee M. A Supplement to the Anophelinae of Africa South of the Sahara. Publ South African Inst Med Res. 1987;55:63.
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:804–11.
Livak KJ. Organization and mapping of a sequence on the Drosopila melanogaster X and Y chromosomes that is transcribed during spermatogenesis. Genetics. 1984;107:611–34.
Bass C, Nikou D, Blagborough AM, Vontas J, Sinden RE, Williamson MS, et al. PCR-based detection of Plasmodium in Anopheles mosquitoes: a comparison of a new high-throughput assay with existing methods. Malar J. 2008;7:177.
Snounou G, Viriyakosol S, Zhu XP, Jarra W, Pinheiro L, do Rosario VE, Thaithong SBK. High sensitivity of detection of human malaria parasites by the use of nested polymerase chain reaction. Mol Biochem Parasitol. 1993;61:315–20.
WHO. Informal consultation report. Geneva: World Health Organization; 1998.
WHO. World malaria report. Geneva: World Health Organization; 2013.
Had2know online statistical software. http://www.had2know.com/academics/hardy-weinberg-equilibrium-calculator-2-alleles.html.
VassarStats online statistical software. http://www.vassarstats.net/odds2x2.htm.
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:595–611.
Sandeu MM, Moussiliou A, Moiroux N, Padonou GG, Massougbodji A, Corbel V, et al. Optimized Pan-species and speciation duplex real-time PCR assays for Plasmodium parasites detection in malaria vectors. PLoS ONE. 2012;7:e52719.
Riveron JM, Osae M, Egyir-Yawson A, Irving H, Ibrahim SS, Wondji CS. Multiple insecticide resistance in the major malaria vector Anopheles funestus in southern Ghana: implications for malaria control. Parasit Vectors. 2016;9:504.
Oyewole IO, Ibidapo CA, Okwa OO, Oduola AO, Adeoye GO, Okoh HI, et al. Species composition and role of Anopheles mosquitoes in malaria transmission along Badagry Axis of Lagos Lagoon, Lagos, Nigeria. Int J Insect Sci. 2010;2:51–7.
Ogen O. The agricultural sector and Nigeria’s development: comparative perspectives from the Brazilian Agro-Industrial Economy, 1960–1995; 2007. p. 184–94.
Oil Boom Era (1971–1977). 2003. Accessed 18 Nov 2016.
Farm irrigation system and system Design fundamentals. p. 61–109.
Reid MC, McKenzie FE. The contribution of agricultural insecticide use to increasing insecticide resistance in African malaria vectors. Malar J. 2016;15:107.
NAFDAC. Guidelines for Registration of Imported Drug Products in Nigeria, National Agency for Food and Drug Administration and Control; 2004. p 1–7.
Akogbéto MC, Djouaka RF, Kindé-Gazard DA. Screening of pesticide residues in soil and water samples from agricultural settings. Malar J. 2006;5:22.
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.
Tene Fossog B, Kopya E, Ndo C, Menze-Djantio B, Costantini C, Njiokou F, et al. Water quality and Anopheles gambiae larval tolerance to pyrethroids in the cities of Douala and Yaounde (Cameroon). J Trop Med. 2012;2012.
Philbert A, Lyantagaye SL, Nkwengulila G. A review of agricultural pesticides use and the selection for resistance to insecticides in malaria vectors. Adv Entomol. 2014;120–8.
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.
Chanda E, Hemingway J, Kleinschmidt I, Rehman AM, Ramdeen V, Phiri FN, et al. Insecticide resistance and the future of malaria control in Zambia. PLoS ONE. 2011;6:e24336.
Okorie PN, Ademowo GO, Irving H, Kelly-Hope LA, Wondji CS. Insecticide susceptibility of Anopheles coluzzii and Anopheles gambiae mosquitoes in Ibadan South–West Nigeria. Med Vet Entomol. 2015;29:44–50.
RD and CSW designed the study; SMA and RA carried out mosquito collection and SMA, RA, GMT and JMR reared the mosquitoes; SMA, GMT, RA and JMR performed insecticide susceptibility and synergist tests; SMA, GMT, RD and HI carried out all laboratory experiments; SMA analysed all the data with assistance from JMR and CSW; MOK offered significant contributions to laboratory works and addressing reviewers comments; AAB gave advise on the study design and offered significant insight to finalize the manuscript; SMA, JMR, CSW and RD wrote the manuscript with contributions from all authors. All authors read, made inputs and approved the final manuscript. The authors wish to state that RD and SMA have equal contribution to the manuscript. All authors read and approved the final manuscript.
We appreciate Akaka-Remo community for their cooperation and Isaac Oyewole and innocent Djegbe for assistance during field work. We thank Eric Tossou, Gareth Weedall, Claude Gande and Murielle Soglo for their technical assistance and relevant advice.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated and analysed during this study are included in the published article.
Ethics approval and consent to participate
Not applicable. However, verbal consent was received from the community and household heads after the study aim and objectives were explained to them. During this research study, we neither used insecticides spraying for mosquito collections nor night collections of mosquitoes on human baits. All mosquitoes were sampled during daytime using electrical aspirators activated with batteries.
This work is supported by the Wellcome Trust grants Reference 099864/Z/12/Z awarded to RD and a Wellcome Trust Senior Research Fellowship in Biomedical Sciences to CSW (101893/Z/13/Z).
Rousseau J. Djouaka and Seun M. Atoyebi contributed equally to this work