Background: Metabolic resistance to insecticides is the biggest threat to the continued effectiveness of malaria vector control. However, its underlying molecular basis, crucial for successful resistance management, remains poorly characterized. Results: Here, we demonstrate that the single amino acid change L119F in an upregulated glutathione S-transferase gene, GSTe2, confers high levels of metabolic resistance to DDT in the malaria vector Anopheles funestus. Genome-wide transcription analysis revealed that GSTe2 was the most over-expressed detoxification gene in DDT and permethrin-resistant mosquitoes from Benin. Transgenic expression of GSTe2 in Drosophila melanogaster demonstrated that over-transcription of this gene alone confers DDT resistance and cross-resistance to pyrethroids. Analysis of GSTe2 polymorphism established that the point mutation is tightly associated with metabolic resistance to DDT and its geographical distribution strongly correlates with DDT resistance patterns across Africa. Functional characterization of recombinant GSTe2 further supports the role of the L119F mutation, with the resistant allele being more efficient at metabolizing DDT than the susceptible one. Importantly, we also show that GSTe2 directly metabolizes the pyrethroid permethrin. Structural analysis reveals that the mutation confers resistance by enlarging the GSTe2 DDT-binding cavity, leading to increased DDT access and metabolism. Furthermore, we show that GSTe2 is under strong directional selection in resistant populations, and a restriction of gene flow is observed between African regions, enabling the prediction of the future spread of this resistance. Conclusions: This first DNA-based metabolic resistance marker in mosquitoes provides an essential tool to track the evolution of resistance and to design suitable resistance management strategies.

Pyrethroid insecticides are critical for malaria control in Africa. However, resistance to this insecticide class in the malaria vector Anopheles funestus is spreading rapidly across Africa, threatening the success of ongoing and future malaria control programs. The underlying resistance mechanisms driving the spread of this resistance in wild populations remain largely unknown. Here, we show that increased expression of two tandemly duplicated P450 genes, CYP6P9a and CYP6P9b , is the main mechanism driving pyrethroid resistance in Malawi and Mozambique, two southern African countries where this insecticide class forms the mainstay of malaria control. Genome-wide transcription analysis using microarray and quantitative RT-PCR consistently revealed that CYP6P9a and CYP6P9b are the two genes most highly overexpressed (>50-fold; q < 0.01) in permethrin-resistant mosquitoes. Transgenic expression of CYP6P9a and CYP6P9b in Drosophila melanogaster demonstrated that elevated expression of either of these genes confers resistance to both type I (permethrin) and type II (deltamethrin) pyrethroids. Functional characterization of recombinant CYP6P9b confirmed that this protein metabolized both type I (permethrin and bifenthrin) and type II (deltamethrin and Lambda-cyhalothrin) pyrethroids but not DDT. Variability analysis identified that a single allele of each of these genes is predominantly associated with pyrethroid resistance in field populations from both countries, which is suggestive of a single origin of this resistance that has since spread across the region. Urgent resistance management strategies should be implemented in this region to limit a further spread of this resistance and minimize its impact on the success of ongoing malaria control programs.

Metabolic resistance to insecticides such as pyrethroids in mosquito vectors threatens control of malaria in Africa. Unless it is managed, recent gains in reducing malaria transmission could be lost. To improve monitoring and assess the impact of insecticide resistance on malaria control interventions, we elucidated the molecular basis of pyrethroid resistance in the major African malaria vector, Anopheles funestus. We showed that a single cytochrome P450 allele (CYP6P9a_R) in A. funestus reduced the efficacy of insecticide-treated bednets for preventing transmission of malaria in southern Africa. Expression of key insecticide resistance genes was detected in populations of this mosquito vector throughout Africa but varied according to the region. Signatures of selection and adaptive evolutionary traits including structural polymorphisms and cis-regulatory transcription factor binding sites were detected with evidence of selection due to the scale-up of insecticide-treated bednet use. A cis-regulatory polymorphism driving the overexpression of the major resistance gene CYP6P9a allowed us to design a DNA-based assay for cytochrome P450–mediated resistance to pyrethroid insecticides. Using this assay, we tracked the spread of pyrethroid resistance and found that it was almost fixed in mosquitoes from southern Africa but was absent from mosquitoes collected elsewhere in Africa. Furthermore, a field study in experimental huts in Cameroon demonstrated that mosquitoes carrying the resistance CYP6P9a_R allele survived and succeeded in blood feeding more often than did mosquitoes that lacked this allele. Our findings highlight the need to introduce a new generation of insecticide-treated bednets for malaria control that do not rely on pyrethroid insecticides.

Establishing the extent, geographical distribution and mechanisms of insecticide resistance in malaria vectors is a prerequisite for resistance management. Here, we report a widespread distribution of insecticide resistance in the major malaria vector An. funestus across Uganda and western Kenya under the control of metabolic resistance mechanisms. Female An. funestus collected throughout Uganda and western Kenya exhibited a Plasmodium infection rate between 4.2 to 10.4%. Widespread resistance against both type I (permethrin) and II (deltamethrin) pyrethroids and DDT was observed across Uganda and western Kenya. All populations remain highly susceptible to carbamate, organophosphate and dieldrin insecticides. Knockdown resistance plays no role in the pyrethroid and DDT resistance as no kdr mutation associated with resistance was detected despite the presence of a F1021C replacement. Additionally, no signature of selection was observed on the sodium channel gene. Synergist assays and qRT-PCR indicated that metabolic resistance plays a major role notably through elevated expression of cytochrome P450s. DDT resistance mechanisms differ from West Africa as the L119F-GSTe2 mutation only explains a small proportion of the genetic variance to DDT resistance. The extensive distribution of pyrethroid and DDT resistance in East African An. funestus populations represents a challenge to the control of this vector. However, the observed carbamate and organophosphate susceptibility offers alternative solutions for resistance management.

Pyrethroid resistance in the major malaria vector Anopheles funestus is rapidly expanding across Southern Africa. It remains unknown whether this resistance has a unique origin with the same molecular basis or is multifactorial. Knowledge of the origin, mechanisms and evolution of resistance are crucial to designing successful resistance management strategies. Here, we established the resistance profile of a Zambian An. funestus population at the northern range of the resistance front. Similar to other Southern African populations, Zambian An. funestus mosquitoes are resistant to pyrethroids and carbamate, but in contrast to populations in Mozambique and Malawi, these insects are also DDT resistant. Genome-wide microarray-based transcriptional profiling and qRT-PCR revealed that the cytochrome P450 gene CYP6M7 is responsible for extending pyrethroid resistance northwards. Indeed, CYP6M7 is more over-expressed in Zambia [fold-change (FC) 37.7; 13.2 for qRT-PCR] than CYP6P9a (FC15.6; 8.9 for qRT-PCR) and CYP6P9b (FC11.9; 6.5 for qRT-PCR), whereas CYP6P9a and CYP6P9b are more highly over-expressed in Malawi and Mozambique. Transgenic expression of CYP6M7 in Drosophila melanogaster coupled with in vitro assays using recombinant enzymes and assessments of kinetic properties demonstrated that CYP6M7 is as efficient as CYP6P9a and CYP6P9b in conferring pyrethroid resistance. Polymorphism patterns demonstrate that these genes are under contrasting selection forces: the exceptionally diverse CYP6M7 likely evolves neutrally, whereas CYP6P9a and CYP6P9b are directionally selected. The higher variability of CYP6P9a and CYP6P9b observed in Zambia supports their lesser role in resistance in this country. Pyrethroid resistance in Southern Africa probably has multiple origins under different evolutionary forces, which may necessitate the design of different resistance management strategies.

Despite the recent progress in establishing the patterns of insecticide resistance in the major malaria vector Anopheles funestus, Central African populations of this species remain largely uncharacterised. To bridge this important gap and facilitate the implementation of suitable control strategies against this vector, we characterised the resistance patterns of An. funestus population from northern Cameroon. Collection of indoor-resting female mosquitoes in Gounougou (northern Cameroon) in 2012 and 2015 revealed a predominance of An. funestus during dry season. WHO bioassays performed using F₁ An. funestus revealed that the population was multiple resistant to several insecticide classes including pyrethroids (permethrin, deltamethrin, lambda-cyhalothrin and etofenprox), carbamates (bendiocarb) and organochlorines (DDT and dieldrin). However, a full susceptibility was observed against the organophosphate malathion. Bioassays performed with 2015 collection revealed that resistance against pyrethroids and DDT is increasing. PBO synergist assays revealed a significant recovery of susceptibility for all pyrethroids but less for DDT. Analysis of the polymorphism of a portion of the voltage-gated sodium channel gene (VGSC) revealed the absence of the L1014F/S kdr mutation but identified 3 novel amino acid changes I877L, V881L and A1007S. However, no association was established between VGSC polymorphism and pyrethroid/DDT resistance. The DDT resistant 119F-GSTe2 allele (52%) and the dieldrin resistant 296S-RDL allele (45%) were detected in Gounougou. Temporal analysis between 2006, 2012 and 2015 collections revealed that the 119F-GSTe2 allele was relatively stable whereas a significant decrease is observed for 296S-RDL allele. This multiple resistance coupled with the temporal increased in resistance intensity highlights the need to take urgent measures to prolong the efficacy of current insecticide-based interventions against An. funestus in this African region.