Functional Validation of Endogenous Redox Partner Cytochrome P450 Reductase Reveals the Key P450s CYP6P9a/-b as Broad Substrate Metabolizers Conferring Cross-Resistance to Different Insecticide Classes in Anopheles funestus.

The versatility of cytochrome P450 reductase (CPR) in transferring electrons to P450s from other closely related species has been extensively exploited, e.g., by using An. gambiae CPR (AgCPR), as a homologous surrogate, to validate the role of An. funestus P450s in insecticide resistance. However, genomic variation between the AgCPR and An. funestus CPR (AfCPR) suggests that the full metabolism spectrum of An. funestus P450s might be missed when using AgCPR. To test this hypothesis, we expressed AgCPR and AfCPR side-by-side with CYP6P9a and CYP6P9b and functionally validated their role in the detoxification of insecticides from five different classes. Major variations were observed within the FAD- and NADP-binding domains of AgCPR and AfCPR, e.g., the coordinates of the second FAD stacking residue AfCPR-Y456 differ from that of AgCPR-His456. While no significant differences were observed in the cytochrome c reductase activities, when co-expressed with their endogenous AfCPR, the P450s significantly metabolized higher amounts of permethrin and deltamethrin, with CYP6P9b-AfCPR membrane metabolizing α-cypermethrin as well. Only the CYP6P9a-AfCPR membrane significantly metabolized DDT (producing dicofol), bendiocarb, clothianidin, and chlorfenapyr (bioactivation into tralopyril). This demonstrates the broad substrate specificity of An. funestus CYP6P9a/-b, capturing their role in conferring cross-resistance towards unrelated insecticide classes, which can complicate resistance management.

Gene Conversion Explains Elevated Diversity in the Immunity Modulating APL1 Gene of the Malaria Vector Anopheles funestus.

Leucine-rich repeat proteins and antimicrobial peptides are the key components of the innate immune response to Plasmodium and other microbial pathogens in Anopheles mosquitoes. The APL1 gene of the malaria vector Anopheles funestus has exceptional levels of non-synonymous polymorphism across the range of An. funestus, with an average πn of 0.027 versus a genome-wide average of 0.002, and πn is consistently high in populations across Africa. Elevated APL1 diversity was consistent between the independent pooled-template and target-enrichment datasets, however no link between APL1 diversity and insecticide resistance was observed. Although lacking the diversity of APL1, two further mosquito innate-immunity genes of the gambicin anti-microbial peptide family had πn/πs ratios greater than one, possibly driven by either positive or balancing selection. The cecropin antimicrobial peptides were expressed much more highly than other anti-microbial peptide genes, a result discordant with current models of anti-microbial peptide activity. The observed APL1 diversity likely results from gene conversion between paralogues, as evidenced by shared polymorphisms, overlapping read mappings, and recombination events among paralogues. In conclusion, we hypothesize that higher gene expression of APL1 than its paralogues is correlated with a more open chromatin formation, which enhances gene conversion and elevated diversity at this locus.

Multi-omics analysis identifies a CYP9K1 haplotype conferring pyrethroid resistance in the malaria vector Anopheles funestus in East Africa.

Metabolic resistance to pyrethroids is a menace to the continued effectiveness of malaria vector controls. Its molecular basis is complex and varies geographically across Africa. Here, we used a multi-omics approach, followed-up with functional validation to show that a directionally selected haplotype of a cytochrome P450, CYP9K1 is a major driver of resistance in Anopheles funestus. A PoolSeq GWAS using mosquitoes alive and dead after permethrin exposure, from Malawi and Cameroon, detected candidate genomic regions, but lacked consistency across replicates. Targeted sequencing of candidate resistance genes detected several SNPs associated with known pyrethroid resistance QTLs. The most significant SNPs were in the cytochrome P450 CYP304B1 (Cameroon), CYP315A1 (Uganda) and the ABC transporter gene ABCG4 (Malawi). However, when comparing field resistant mosquitoes to laboratory susceptible, the pyrethroid resistance locus rp1 and SNPs around the ABC transporter ABCG4 were consistently significant, except for Uganda where SNPs in the P450 CYP9K1 was markedly significant. In vitro heterologous metabolism assays with recombinant CYP9K1 revealed that it metabolises type II pyrethroid (deltamethrin; 64% depletion) but not type I (permethrin; 0%), while moderately metabolising DDT (17%). CYP9K1 exhibited reduced genetic diversity in Uganda underlying an extensive selective sweep. Furthermore, a glycine to alanine (G454A) amino acid change in CYP9K1 was fixed in Ugandan mosquitoes but not in other An. funestus populations. This study sheds further light on the evolution of metabolic resistance in a major malaria vector by implicating more genes and variants that can be used to design field-applicable markers to better track resistance Africa-wide.

Genome-Wide Transcriptional Analysis and Functional Validation Linked a Cluster of Epsilon Glutathione S-Transferases with Insecticide Resistance in the Major Malaria Vector Anopheles funestus across Africa.

Resistance is threatening the effectiveness of insecticide-based interventions in use for malaria control. Pinpointing genes associated with resistance is crucial for evidence-based resistance management targeting the major malaria vectors. Here, a combination of RNA-seq based genome-wide transcriptional analysis and RNA-silencing in vivo functional validation were used to identify key insecticide resistance genes associated with DDT and DDT/permethrin cross-resistance across Africa. A cluster of glutathione-S-transferase from epsilon group were found to be overexpressed in resistant populations of Anopheles funestus across Africa including GSTe1 [Cameroon (fold change, FC: 2.54), Ghana (4.20), Malawi (2.51)], GSTe2 [Cameroon (4.47), Ghana (7.52), Malawi (2.13)], GSTe3 [Cameroon (2.49), Uganda (2.60)], GSTe4 in Ghana (3.47), GSTe5 [Ghana (2.94), Malawi (2.26)], GSTe6 [Cameroun (3.0), Ghana (3.11), Malawi (3.07), Uganda (3.78)] and GSTe7 (2.39) in Ghana. Validation of GSTe genes expression profiles by qPCR confirmed that the genes are differentially expressed across Africa with a greater overexpression in DDT-resistant mosquitoes. RNAi-based knock-down analyses supported that five GSTe genes are playing a major role in resistance to pyrethroids (permethrin and deltamethrin) and DDT in An. funestus, with a significant recovery of susceptibility observed when GSTe2, 3, 4, 5 and GSTe6 were silenced. These findings established that GSTe3, 4, 5 and 6 contribute to DDT resistance and should be further characterized to identify their specific genetic variants, to help design DNA-based diagnostic assays, as previously done for the 119F-GSTe2 mutation. This study highlights the role of GSTes in the development of resistance to insecticides in malaria vectors and calls for actions to mitigate this resistance.

Influence of a Major Mountainous Landscape Barrier (Mount Cameroon) on the Spread of Metabolic (GSTe2) and Target-Site (Rdl) Resistance Alleles in the African Malaria Vector Anopheles funestus.

Increased levels of insecticide resistance in major malaria vectors such as Anopheles funestus threaten the effectiveness of insecticide-based control programmes. Understanding the landscape features impacting the spread of resistance makers is necessary to design suitable resistance management strategies. Here, we examined the influence of the highest mountain in West Africa (Mount Cameroon; 4095 m elevation) on the spread of metabolic and target-site resistance alleles in An. funestus populations. Vector composition varied across the four localities surveyed along the altitudinal cline with major vectors exhibiting high parity rate (80.5%). Plasmodium infection rates ranged from 0.79% (An. melas) to 4.67% (An. funestus). High frequencies of GSTe2R (67-81%) and RdlR (49-90%) resistance alleles were observed in An. funestus throughout the study area, with GSTe2R frequency increasing with altitude, whereas the opposite is observed for RdlR. Patterns of genetic diversity and population structure analyses revealed high levels of polymorphisms with 12 and 16 haplotypes respectively for GSTe2 and Rdl. However, the reduced diversity patterns of resistance allele carriers revealed signatures of positive selection on the two genes across the study area irrespective of the altitude. Despite slight variations associated with the altitude, the spread of resistance alleles suggest that control strategies could be implemented against malaria vectors across mountainous landscapes.

Investigation of DDT resistance mechanisms in Anopheles funestus populations from northern and southern Benin reveals a key role of the GSTe2 gene.

BACKGROUND: Understanding the molecular basis of insecticide resistance in mosquito, such as Anopheles funestus, is an important step in developing strategies to mitigate the resistance problem. This study aims to assess the role of the GSTe2 gene in DDT resistance and determine the genetic diversity of this gene in An. funestus. METHODS: Gene expression analysis was performed using microarrays and PCR while the potential mutation associated with resistance was determined using sequencing. RESULTS: Low expression level of GSTe2 gene was recorded in Burkina-Faso samples with a fold change of 3.3 while high expression (FC 35.6) was recorded in southern Benin in Pahou (FC 35.6) and Kpome (FC 13.3). The sequencing of GSTe2 gene in six localities showed that L119F-GSTe2 mutation is almost getting fixed in highly DDT-resistant Benin (Pahou, Kpome, Doukonta) and Nigeria (Akaka Remo) mosquitoes with a low mutation rate observed in Tanongou (Benin) and Burkina-Faso mosquitoes. CONCLUSION: This study shows the key role of the GSTe2 gene in DDT resistant An. funestus in Benin. Polymorphism analysis of this gene across Benin revealed possible barriers to gene flow, which could impact the design and implementation of resistance management strategies in the country.

A 6.5-kb intergenic structural variation enhances P450-mediated resistance to pyrethroids in malaria vectors lowering bed net efficacy.

Elucidating the complex evolutionary armory that mosquitoes deploy against insecticides is crucial to maintain the effectiveness of insecticide-based interventions. Here, we deciphered the role of a 6.5-kb structural variation (SV) in driving cytochrome P450-mediated pyrethroid resistance in the malaria vector, Anopheles funestus. Whole-genome pooled sequencing detected an intergenic 6.5-kb SV between duplicated CYP6P9a/b P450s in pyrethroid-resistant mosquitoes through a translocation event. Promoter analysis revealed a 17.5-fold higher activity (p < .0001) for the SV- carrying fragment than the SV- free one. Quantitative real-time PCR expression profiling of CYP6P9a/b for each SV genotype supported its role as an enhancer because SV+/SV+ homozygote mosquitoes had a significantly greater expression for both genes than heterozygotes SV+/SV- (1.7- to 2-fold) and homozygotes SV-/SV- (4-to 5-fold). Designing a PCR assay revealed a strong association between this SV and pyrethroid resistance (SV+/SV+ vs. SV-/SV-; odds ratio [OR] = 2,079.4, p < .001). The 6.5-kb SV is present at high frequency in southern Africa (80%-100%) but absent in East/Central/West Africa. Experimental hut trials revealed that homozygote SV mosquitoes had a significantly greater chance to survive exposure to pyrethroid-treated nets (OR 27.7; p < .0001) and to blood feed than susceptible mosquitoes. Furthermore, mosquitoes homozygote-resistant at the three loci (SV+/CYP6P9a_R/CYP6P9b_R) exhibited a higher resistance level, leading to a far superior ability to survive exposure to nets than those homozygotes susceptible at the three loci, revealing a strong additive effect. This study highlights the important role of structural variations in the development of insecticide resistance in malaria vectors and their detrimental impact on the effectiveness of pyrethroid-based nets.

Investigating the molecular basis of multiple insecticide resistance in a major malaria vector Anopheles funestus (sensu stricto) from Akaka-Remo, Ogun State, Nigeria.

BACKGROUND: Understanding the mechanisms used by Anopheles mosquitoes to survive insecticide exposure is key to manage existing insecticide resistance and develop more suitable insecticide-based malaria vector control interventions as well as other alternative integrated tools. To this regard, the molecular basis of permethrin, DDT and dieldrin resistance in Anopheles funestus (sensu stricto) at Akaka-Remo was investigated. METHODS: Bioassays were conducted on 3-5-day-old adult An. funestus (s.s.) mosquitoes for permethrin, DDT and dieldrin susceptibility test. The molecular mechanisms of mosquito resistance to these insecticides were investigated using microarray and reverse transcriptase PCR techniques. The voltage-gated sodium channel region of mosquitoes was also screened for the presence of knockdown resistance mutations (kdr west and east) by sequencing method. RESULTS: Anopheles funestus (s.s.) population was resistant to permethrin (mortality rate of 68%), DDT (mortality rate of 10%) and dieldrin (mortality rate of 8%) insecticides. Microarray and RT-PCR analyses revealed the overexpression of glutathione S-transferase genes, cytochrome P450s, esterase, trypsin and cuticle proteins in resistant mosquitoes compared to control. The GSTe2 was the most upregulated detoxification gene in permethrin-resistant (FC = 44.89), DDT-resistant (FC = 57.39) and dieldrin-resistant (FC = 41.10) mosquitoes compared to control population (FC = 22.34). The cytochrome P450 gene, CYP6P9b was also upregulated in both permethrin- and DDT-resistant mosquitoes. The digestive enzyme, trypsin (hydrolytic processes) and the cuticle proteins (inducing cuticle thickening leading to reduced insecticides penetration) also showed high involvement in insecticide resistance, through their overexpression in resistant mosquitoes compared to control. The kdr east and west were absent in all mosquitoes analysed, suggesting their non-involvement in the observed mosquito resistance. CONCLUSIONS: The upregulation of metabolic genes, especially the GSTe2 and trypsin, as well as the cuticle proteins is driving insecticide resistance of An. funestus (s.s.) population. However, additional molecular analyses, including functional metabolic assays of these genes as well as screening for a possible higher cuticular hydrocarbon and lipid contents, and increased procuticle thickness in resistant mosquitoes are needed to further describe their distinct roles in mosquito resistance.

An Africa-wide genomic evolution of insecticide resistance in the malaria vector Anopheles funestus involves selective sweeps, copy number variations, gene conversion and transposons.

Insecticide resistance in malaria vectors threatens to reverse recent gains in malaria control. Deciphering patterns of gene flow and resistance evolution in malaria vectors is crucial to improving control strategies and preventing malaria resurgence. A genome-wide survey of Anopheles funestus genetic diversity Africa-wide revealed evidences of a major division between southern Africa and elsewhere, associated with different population histories. Three genomic regions exhibited strong signatures of selective sweeps, each spanning major resistance loci (CYP6P9a/b, GSTe2 and CYP9K1). However, a sharp regional contrast was observed between populations correlating with gene flow barriers. Signatures of complex molecular evolution of resistance were detected with evidence of copy number variation, transposon insertion and a gene conversion between CYP6P9a/b paralog genes. Temporal analyses of samples before and after bed net scale up suggest that these genomic changes are driven by this control intervention. Multiple independent selective sweeps at the same locus in different parts of Africa suggests that local evolution of resistance in malaria vectors may be a greater threat than trans-regional spread of resistance haplotypes.

Exploring the Mechanisms of Multiple Insecticide Resistance in a Highly Plasmodium-Infected Malaria Vector Anopheles funestus Sensu Stricto from Sahel of Northern Nigeria.

The Nigerian Government is scaling up the distribution of insecticide-treated bed nets for malaria control, but the lack of surveillance data, especially in the Sudan/Sahel region of the country, may hinder targeting priority populations. Here, the vectorial role and insecticide resistance profile of a population of a major malaria vector Anopheles funestus sensu stricto from Sahel of Nigeria was characterised. An. funestus s.s. was the only vector found, with a high human blood index (100%) and a biting rate of 5.3/person/night. High Plasmodium falciparum infection was discovered (sporozoite rate = 54.55%). The population is resistant to permethrin (mortality = 48.30%, LT50 = 65.76 min), deltamethrin, DDT (dichlorodiphenyltrichloroethane) and bendiocarb, with mortalities of 29.44%, 56.34% and 54.05%, respectively. Cone-bioassays established loss of efficacy of the pyrethroid-only long-lasting insecticidal nets (LLINs); but 100% recovery of susceptibility was obtained for piperonylbutoxide (PBO)-containing PermaNet®3.0. Synergist bioassays with PBO and diethyl maleate recovered susceptibility, implicating CYP450s (permethrin mortality = 78.73%, χ2 = 22.33, P < 0.0001) and GSTs (DDT mortality = 81.44%, χ2 = 19.12, P < 0.0001). A high frequency of 119F GSTe2 mutation (0.84) was observed (OR = 16, χ2 = 3.40, P = 0.05), suggesting the preeminent role of metabolic resistance. These findings highlight challenges associated with deployment of LLINs and indoor residual spraying (IRS) in Nigeria.

An Experimental Hut Evaluation of PBO-Based and Pyrethroid-Only Nets against the Malaria Vector Anopheles funestus Reveals a Loss of Bed Nets Efficacy Associated with GSTe2 Metabolic Resistance.

Growing insecticide resistance in malaria vectors is threatening the effectiveness of insecticide-based interventions, including Long Lasting Insecticidal Nets (LLINs). However, the impact of metabolic resistance on the effectiveness of these tools remains poorly characterized. Using experimental hut trials and genotyping of a glutathione S-transferase resistance marker (L119F-GSTe2), we established that GST-mediated resistance is reducing the efficacy of LLINs against Anopheles funestus. Hut trials performed in Cameroon revealed that Piperonyl butoxide (PBO)-based nets induced a significantly higher mortality against pyrethroid resistant An. funestus than pyrethroid-only nets. Blood feeding rate and deterrence were significantly higher in all LLINs than control. Genotyping the L119F-GSTe2 mutation revealed that, for permethrin-based nets, 119F-GSTe2 resistant mosquitoes have a greater ability to blood feed than susceptible while the opposite effect is observed for deltamethrin-based nets. For Olyset Plus, a significant association with exophily was observed in resistant mosquitoes (OR = 11.7; p < 0.01). Furthermore, GSTe2-resistant mosquitoes (cone assays) significantly survived with PermaNet 2.0 (OR = 2.1; p < 0.01) and PermaNet 3.0 (side) (OR = 30.1; p < 0.001) but not for Olyset Plus. This study shows that the efficacy of PBO-based nets (e.g., blood feeding inhibition) against pyrethroid resistant malaria vectors could be impacted by other mechanisms including GST-mediated metabolic resistance not affected by the synergistic action of PBO. Mosaic LLINs incorporating a GST inhibitor (diethyl maleate) could help improve their efficacy in areas of GST-mediated resistance.

Cis-regulatory CYP6P9b P450 variants associated with loss of insecticide-treated bed net efficacy against Anopheles funestus.

Elucidating the genetic basis of metabolic resistance to insecticides in malaria vectors is crucial to prolonging the effectiveness of insecticide-based control tools including long lasting insecticidal nets (LLINs). Here, we show that cis-regulatory variants of the cytochrome P450 gene, CYP6P9b, are associated with pyrethroid resistance in the African malaria vector Anopheles funestus. A DNA-based assay is designed to track this resistance that occurs near fixation in southern Africa but not in West/Central Africa. Applying this assay we demonstrate, using semi-field experimental huts, that CYP6P9b-mediated resistance associates with reduced effectiveness of LLINs. Furthermore, we establish that CYP6P9b combines with another P450, CYP6P9a, to additively exacerbate the reduced efficacy of insecticide-treated nets. Double homozygote resistant mosquitoes (RR/RR) significantly survive exposure to insecticide-treated nets and successfully blood feed more than other genotypes. This study provides tools to track and assess the impact of multi-gene driven metabolic resistance to pyrethroids, helping improve resistance management.

A marker of glutathione S-transferase-mediated resistance to insecticides is associated with higher Plasmodium infection in the African malaria vector Anopheles funestus.

Metabolic resistance to insecticides is threatening malaria control in Africa. However, the extent to which it impacts malaria transmission remains unclear. Here, we investigated the association between a marker of glutathione S-transferase mediated metabolic resistance and Plasmodium infection in field population of Anopheles funestus s.s. in comparison to the A296S-RDL target site mutation. The 119F-GSTe2 resistant allele was present in southern (Obout) (56%) and central (Mibellon) (25%) regions of Cameroon whereas the 296S-RDL resistant allele was detected at 98.5% and 15% respectively. The whole mosquito Plasmodium and sporozoite infection rates were 57% and 14.8% respectively in Obout (n = 508) and 19.7% and 5% in Mibellon (n = 360). No association was found between L119F-GSTe2 genotypes and whole mosquito infection status. However, when analyzing oocyst and sporozoite infection rates separately, the resistant homozygote 119F/F genotype was significantly more associated with Plasmodium infection in Obout than both heterozygote (OR = 2.5; P = 0.012) and homozygote susceptible (L/L119) genotypes (OR = 2.10; P = 0.013). In contrast, homozygote RDL susceptible mosquitoes (A/A296) were associated more frequently with Plasmodium infection than other genotypes (OR = 4; P = 0.03). No additive interaction was found between L119F and A296S. Sequencing of the GSTe2 gene showed no association between the polymorphism of this gene and Plasmodium infection. Glutathione S-transferase metabolic resistance is potentially increasing the vectorial capacity of resistant An. funestus mosquitoes. This could result in a possible exacerbation of malaria transmission in areas of high GSTe2-based metabolic resistance to insecticides.

Escalation of Pyrethroid Resistance in the Malaria Vector Anopheles funestus Induces a Loss of Efficacy of Piperonyl Butoxide-Based Insecticide-Treated Nets in Mozambique.

BACKGROUND: Insecticide resistance poses a serious threat to insecticide-based interventions in Africa. There is a fear that resistance escalation could jeopardize malaria control efforts. Monitoring of cases of aggravation of resistance intensity and its impact on the efficacy of control tools is crucial to predict consequences of resistance. METHODS: The resistance levels of an Anopheles funestus population from Palmeira, southern Mozambique, were characterized and their impact on the efficacy of various insecticide-treated nets established. RESULTS: A dramatic loss of efficacy of all long-lasting insecticidal nets (LLINs), including piperonyl butoxide (PBO)-based nets (Olyset Plus), was observed. This An. funestus population consistently (2016, 2017, and 2018) exhibited a high degree of pyrethroid resistance. Molecular analyses revealed that this resistance escalation was associated with a massive overexpression of the duplicated cytochrome P450 genes CYP6P9a and CYP6P9b, and also the fixation of the resistance CYP6P9a_R allele in this population in 2016 (100%) in contrast to 2002 (5%). However, the low recovery of susceptibility after PBO synergist assay suggests that other resistance mechanisms could be involved. CONCLUSIONS: The loss of efficacy of pyrethroid-based LLINs with and without PBO is a concern for the effectiveness of insecticide-based interventions, and action should be taken to prevent the spread of such super-resistance.

A cytochrome P450 allele confers pyrethroid resistance on a major African malaria vector, reducing insecticide-treated bednet efficacy.

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.

Exposure to the insecticide-treated bednet PermaNet 2.0 reduces the longevity of the wild African malaria vector Anopheles funestus but GSTe2-resistant mosquitoes live longer.

BACKGROUND: Despite the increased report of insecticide resistance in malaria vectors, its impact on mosquito's life-traits after exposure to insecticide-treated nets remains under investigated. Here, we assessed the effects of exposure to PermaNet 2.0 on several life traits of An. gambiae s.l. and An. funestus s.l. field mosquitoes in Cameroon. METHODOLOGY: Female Anopheles mosquitoes were collected indoor using electric aspirators in southern Cameroon (Obout) in 2016. After assessing the resistance status of F1 from the field collected-mosquitoes, progeny of the first generation (An. funestus s.l.) and seventh generation (An. gambiae s.l.) were used to assess the long-term effect of exposure to PermaNet 2.0 on several life-traits of these vectors (longevity, blood feeding ability, fecundity and fertility) in comparison to untreated net. In addition, the L119F-GSTe2 mutation associated with DDT/pyrethroids resistance in An. funestus was genotyped to assess its association with increased life-span post-exposure. PRINCIPAL FINDINGS: Both An. funestus and An. gambiae were resistant to pyrethroids and DDT with a greater level in the latter. Pyrethroid-only nets PermaNet 2.0 (17.5% mortality) and Olyset (0% mortality) exhibited a significantly reduced efficacy against An. funestus in contrast to a greater efficacy for PBO-based Nets Olyset Plus (65% mortality), PermaNet 3.0 top (100% mortality). In both species, mosquitoes that survived exposure to PermaNet 2.0 exhibited a significantly reduced longevity than those non-exposed (6.95 days vs 12.46 for An. funestus P<0.001; 8.87 vs 11.25 days for An. gambiae; P<0.001). However, no significant difference was observed for blood feeding and fecundity in both species. In addition, molecular analysis of the L119F-GSTe2 mutation revealed that this mutation is associated with an increase in the chance of surviving after exposure to this net in An. funestus. CONCLUSIONS: These results show that although the PermaNet 2.0 presents a reduced efficacy against resistant populations, it remains efficient after exposure by reducing the life expectancy of the vectors which could contribute in the reduction of malaria incidence.

Escalation of Pyrethroid Resistance in the Malaria Vector Anopheles funestus Induces a Loss of Efficacy of Piperonyl Butoxide­Based Insecticide-Treated Nets in Mozambique

Insecticide resistance poses a serious threat to insecticide-based interventions in Africa. There is a fear that resistance escalation could jeopardize malaria control efforts. Monitoring of cases of aggravation of resistance intensity and its impact on the efficacy of control tools is crucial to predict consequences of resistance. The resistance levels of an Anopheles funestus population from Palmeira, southern Mozambique, were characterized and their impact on the efficacy of various insecticide-treated nets established. A dramatic loss of efficacy of all long-lasting insecticidal nets (LLINs), including piperonyl butoxide (PBO)­based nets (Olyset Plus), was observed. This An. funestus population consistently (2016, 2017, and 2018) exhibited a high degree of pyrethroid resistance. Molecular analyses revealed that this resistance escalation was associated with a massive overexpression of the duplicated cytochrome P450 genes CYP6P9a and CYP6P9b, and also the fixation of the resistance CYP6P9a_R allele in this population in 2016 (100%) in contrast to 2002 (5%). However, the low recovery of susceptibility after PBO synergist assay suggests that other resistance mechanisms could be involved. The loss of efficacy of pyrethroid-based LLINs with and without PBO is a concern for the effectiveness of insecticide-based interventions, and action should be taken to prevent the spread of such super-resistance

Contrasting patterns of gene expression indicate differing pyrethroid resistance mechanisms across the range of the New World malaria vector Anopheles albimanus.

Decades of unmanaged insecticide use and routine exposure to agrochemicals have left many populations of malaria vectors in the Americas resistant to multiple classes of insecticides, including pyrethroids. The molecular basis of pyrethroid resistance is relatively uncharacterised in American malaria vectors, preventing the design of suitable resistance management strategies. Using whole transcriptome sequencing, we characterized the mechanisms of pyrethroid resistance in Anopheles albimanus from Peru and Guatemala. An. albimanus were phenotyped as either deltamethrin or alpha-cypermethrin resistant. RNA from 1) resistant, 2) unexposed, and 3) a susceptible laboratory strain of An. albimanus was sequenced and analyzed using RNA-Seq. Expression profiles of the three groups were compared based on the current annotation of the An. albimanus reference genome. Several candidate genes associated with pyrethroid resistance in other malaria vectors were found to be overexpressed in resistant An. albimanus. In addition, gene ontology terms related to serine-type endopeptidase activity, extracellular activity and chitin metabolic process were also commonly overexpressed in the field caught resistant and unexposed samples from both Peru and Guatemala when compared to the susceptible strain. The cytochrome P450 CYP9K1 was overexpressed 14x in deltamethrin and 8x in alpha-cypermethrin-resistant samples from Peru and 2x in deltamethrin-resistant samples from Guatemala, relative to the susceptible laboratory strain. CYP6P5 was overexpressed 68x in deltamethrin-resistant samples from Peru but not in deltamethrin-resistant samples from Guatemala. When comparing overexpressed genes between deltamethrin-resistant and alpha-cypermethrin-resistant samples from Peru, a single P450 gene, CYP4C26, was overexpressed 9.8x (p<0.05) in alpha-cypermethrin-resistant samples. In Peruvian deltamethrin-resistant samples, the knockdown resistance mutation (kdr) variant alleles at position 1014 were rare, with approximately 5% frequency, but in the alpha-cypermethrin-resistant samples, the frequency of these alleles was approximately 15-30%. Functional validation of the candidate genes and the kdr mutation as a resistance marker for alpha-cypermethrin will confirm the role of these mechanisms in conferring pyrethroid resistance.

Fitness Costs of the Glutathione S-Transferase Epsilon 2 (L119F-GSTe2) Mediated Metabolic Resistance to Insecticides in the Major African Malaria Vector Anopheles Funestus.

Metabolic resistance to insecticides threatens malaria control. However, little is known about its fitness cost in field populations of malaria vectors, thus limiting the design of suitable resistance management strategies. Here, we assessed the association between the glutathione S-transferase GSTe2-mediated metabolic resistance and life-traits of natural populations of Anopheles funestus. A total of 1200 indoor resting blood-fed female An. funestus (F0) were collected in Mibellon, Cameroon (2016/2017), and allowed to lay eggs individually. Genotyping of F1 mosquitoes for the L119F-GSTE2 mutation revealed that L/L119-homozygote susceptible (SS) mosquitoes significantly laid more eggs than heterozygotes L119F-RS (odds ratio (OR) = 2.06; p < 0.0001) and homozygote resistant 119F/F-RR (OR = 2.93; p < 0.0001). L/L119-SS susceptible mosquitoes also showed the higher ability for oviposition than 119F/F-RR resistant (OR = 2.68; p = 0.0002) indicating a reduced fecundity in resistant mosquitoes. Furthermore, L119F-RS larvae developed faster (nine days) than L119F-RR and L119F-SS (11 days) (X² = 11.052; degree of freedom (df) = 4; p = 0.02) suggesting a heterozygote advantage effect for larval development. Interestingly, L/L119-SS developed faster than 119F/F-RR (OR = 5.3; p < 0.0001) revealing an increased developmental time in resistant mosquitoes. However, genotyping and sequencing revealed that L119F-RR mosquitoes exhibited a higher adult longevity compared to RS (OR > 2.2; p < 0.05) and SS (OR > 2.1; p < 0.05) with an increased frequency of GSTe2-resistant haplotypes in mosquitoes of D30 after adult emergence. Additionally, comparison of the expression of GSTe2 revealed a significantly increased expression from D1-D30 after emergence of adults (Anova test (F) = 8; df= 3; p = 0.008). The negative association between GSTe2 and some life traits of An. funestus could facilitate new resistance management strategies. However, the increased longevity of GSTe2-resistant mosquitoes suggests that an increase in resistance could exacerbate malaria transmission.

Molecular basis of permethrin and DDT resistance in an Anopheles funestus population from Benin.

BACKGROUND: Insecticide resistance in Anopheles mosquitoes is threatening the success of malaria control programmes. In order to implement suitable insecticide resistance management strategies, it is necessary to understand the underlying mechanisms involved. To achieve this, the molecular basis of permethrin and DDT resistance in the principal malaria vector, Anopheles funestus from inland Benin (Kpome), was investigated. RESULTS: Here, using a microarray-based genome-wide transcription and qRT-PCR analysis, we showed that metabolic resistance mechanisms through over-expression of cytochrome P450 and glutathione S-transferase genes (GSTs) are a major contributor to DDT and permethrin resistance in Anopheles funestus from Kpome. The GSTe2 gene was the most upregulated detoxification gene in both DDT- [fold-change (FC: 16.0)] and permethrin-resistant (FC: 18.1) mosquitoes suggesting that upregulation of this gene could contribute to DDT resistance and cross-resistance to permethrin. CYP6P9a and CYP6P9b genes that have been previously associated with pyrethroid resistance were also significantly overexpressed with FC 5.4 and 4.8, respectively, in a permethrin resistant population. Noticeably, the GSTs, GSTd1-5 and GSTd3, were more upregulated in DDT-resistant than in permethrin-resistant Anopheles funestus suggesting these genes are more implicated in DDT resistance. The absence of the L1014F or L1014S kdr mutations in the voltage-gated sodium channel gene coupled with the lack of directional selection at the gene further supported that knockdown resistance plays little role in this resistance. CONCLUSIONS: The major role played by metabolic resistance to pyrethroids in this An. funestus population in Benin suggests that using novel control tools combining the P450 synergist piperonyl butoxide (PBO), such as PBO-based bednets, could help manage the growing pyrethroid resistance in this malaria vector in Benin.