Pyrethroid Resistance in the Major Malaria Vector Anopheles funestus is Exacerbated by Overexpression and Overactivity of the P450 CYP6AA1 Across Africa

A cub and sushi domain-containing protein with esterase-like activity confers insecticide resistance in the Indian malaria vector Anopheles stephensi

Chemical insecticides (organophosphates and pyrethroids) in the form of IRS (Indoor Residual Sprays) and LLINs (Long Lasting Insecticidal Nets) are the cornerstone for vector control, globally. However, their incessant use has resulted in widespread development of resistance in mosquito vectors, warranting continuous monitoring and investigation of the underlying mechanisms of resistance. Here, we identified a previously uncharacterized- Cub and Sushi Domain containing Insecticide Resistance (CSDIR) protein and generated evidence for its role in mediating insecticide resistance in the Anopheles stephensi. A strong binding affinity of the CSDIR protein towards different classes of insecticide molecules-malathion (KD 6.43 μM) and deltamethrin (KD 46.7 μM) were demonstrated using MD simulation studies and Surface Plasmon Resonance (SPR) experiments. Further, the recombinant CSDIR913-1190 protein exhibited potent esterase-like activity (α-naphthyl acetate (α-NA)- 1.356 ± 0.262 mM/min/mg and β-naphthyl acetate (β -NA)- 1.777 ± 0.220 mM/min/mg). Interestingly, dsRNA-mediated gene silencing of the CSDIR transcripts caused >60% mortality in resistant An. stephensi upon 1-h exposure to deltamethrin and malathion insecticides, compared to the control group. A significant reduction in the esterase-like activity was also observed against α-NA (p = 0.004) and β-NA (p = 0.025) in CSDIR silenced mosquitoes compared to the control group. Using computational analysis and experimental data, our results provided significant evidence of the involvement of the CSDIR protein in mediating insecticide resistance in Anopheles mosquitoes. Thereby making the CSDIR protein, a novel candidate for exploration of novel insecticide molecules. These data would also be helpful in further understanding the development of metabolic resistance by the Anopheles vector.

Phenotypic evidence of deltamethrin resistance and identification of selective sweeps in Anopheles mosquitoes on the Bijagós Archipelago, Guinea-Bissau

Vector control in the Bijagós Archipelago of Guinea-Bissau currently relies on pyrethroid insecticide-treated nets. However, data on insecticide resistance in Guinea-Bissau is limited. This study identified deltamethrin resistance in the Anopheles gambiae sensu lato complex on Bubaque island using WHO tube tests in November 2022. Whole genome sequencing of An. gambiae sensu stricto mosquitoes identified six single nucleotide polymorphisms (SNPs) previously associated with, or putatively associated with, insecticide resistance: T791M, L995F, N1570Y, A1746S and P1874L in the vgsc gene, and L119V in the gste2 gene. Twenty additional non-synonymous SNPs were identified in insecticide-resistance associated genes. Four of these SNPs were present at frequencies over 5% in the population: T154S, I126F and G26S in the vgsc gene and A65S in ace1. Genome wide selection scans using Garud's H12 statistic identified two selective sweeps: one in chromosome X and one in chromosome 2R. Both selective sweeps overlap with metabolic genes previously associated with insecticide resistance, including cyp9k1 and the cyp6aa/cyp6p gene cluster. This study presents the first phenotypic testing for deltamethrin resistance and the first whole genome sequence data for Anophelesgambiae mosquitoes from the Bijagós, contributing data of significance for vector control policy in this region.

Cryptic population structure and insecticide resistance in Anopheles gambiae from the southern Democratic Republic of Congo

The Democratic Republic of Congo (DRC) suffers from one of the highest malaria burdens worldwide, but information on its Anopheles vector populations is relatively limited. Preventative malaria control in DRC is reliant on pyrethroid-treated nets, raising concerns over the potential impacts of insecticide resistance. We sampled Anopheles gambiae from three geographically distinct populations (Kimpese, Kapolowe and Mikalayi) in southern DRC, collecting from three sub-sites per population and characterising mosquito collections from each for resistance to pyrethroids using WHO tube bioassays. Resistance to each of three different pyrethroids was generally high in An. gambiae with < 92% mortality in all tests, but varied between collections, with mosquitoes from Kimpese being the most resistant. Whole genome sequencing of 165 An. gambiae revealed evidence for genetic differentiation between Kimpese and Kapolowe/Mikalayi, but not between the latter two sample sites despite separation of approximately 800 km. Surprisingly, there was evidence of population structure at a small spatial scale between collection subsites in Kimpese, despite separation of just tens of kilometres. Intra-population (H12) and inter-population (FST) genome scans identified multiple peaks corresponding to genes associated with insecticide resistance such as the voltage gated sodium channel (Vgsc) target site on chromosome 2L, a Cyp6 cytochrome P450 cluster on chromosome arm 2R, and the Cyp9k1 P450 gene on chromosome X. In addition, in the Kimpese subsites, the P450 redox partner gene Cpr showed evidence for contemporary selection (H12) and population differentiation (FST) meriting further exploration as a potential resistance associated marker.

Whole genome sequence analysis of population structure and insecticide resistance markers in Anopheles melas from the Bijagós Archipelago, Guinea-Bissau

BACKGROUND

Anopheles melas is an understudied malaria vector with a potential role in malaria transmission on the Bijagós Archipelago of Guinea-Bissau. This study presents the first whole-genome sequencing and population genetic analysis for this species from the Bijagós. To our knowledge, this also represents the largest population genetic analysis using WGS data from non-pooled An. melas mosquitoes.

METHODS

WGS was conducted for 30 individual An. melas collected during the peak malaria transmission season in 2019 from six different islands on the Bijagós Archipelago. Bioinformatics tools were used to investigate the population structure and prevalence of insecticide resistance markers in this mosquito population.

RESULTS

Insecticide resistance mutations associated with pyrethroid resistance in Anopheles gambiae s.s. from the Bijagós were absent in the An. melas population, and no signatures of selective sweeps were identified in insecticide resistance-associated genes. Analysis of structural variants identified a large duplication encompassing the cytochrome-P450 gene cyp9k1. Phylogenetic analysis using publicly available mitochondrial genomes indicated that An. melas from the Bijagós split into two phylogenetic groups because of differentiation on the mitochondrial genome attributed to the cytochrome C oxidase subunits COX I and COX II and the NADH dehydrogenase subunits 1, 4, 4L and 5.

CONCLUSIONS

This study identified an absence of insecticide-resistant SNPs common to An. gambiae in the An. melas population, but did identify structural variation over insecticide resistance-associated genes. Furthermore, this study presents novel insights into the population structure of this malaria vector using WGS analysis. Additional studies are required to further understand the role of this vector in malaria transmission.

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.

Distribution and insecticide resistance profile of the major malaria vector Anopheles funestus group across the African continent

There has been significant progress in malaria control in the last 2 decades, with a decline in mortality and morbidity. However, these gains are jeopardised by insecticide resistance, which negatively impacts the core interventions, such as insecticide-treated nets (ITN) and indoor residual spraying (IRS). While most malaria control and research efforts are still focused on Anopheles gambiae complex mosquitoes, Anopheles funestus remains an important vector in many countries and, in some cases, contributes to most of the local transmission. As countries move towards malaria elimination, it is important to ensure that all dominant vector species, including An. funestus, an important vector in some countries, are targeted. The objective of this review is to compile and discuss information related to A. funestus populations' resistance to insecticides and the mechanisms involved across Africa, emphasising the sibling species and their resistance profiles in relation to malaria elimination goals. Data on insecticide resistance in An. funestus malaria vectors in Africa were extracted from published studies. Online bibliographic databases, including Google Scholar and PubMed, were used to search for relevant studies. Articles published between 2000 and May 2023 reporting resistance of An. funestus to insecticides and associated mechanisms were included. Those reporting only bionomics were excluded. Spatial variation in species distribution and resistance to insecticides was recorded from 174 articles that met the selection criteria. It was found that An. funestus was increasingly resistant to the four classes of insecticides recommended by the World Health Organisation for malaria vector control; however, this varied by country. Insecticide resistance appears to reduce the effectiveness of vector control methods, particularly IRS and ITN. Biochemical resistance due to detoxification enzymes (P450s and glutathione-S-transferases [GSTs]) in An. funestus was widely recorded. However, An. funestus in Africa remains susceptible to other insecticide classes, such as organophosphates and neonicotinoids. This review highlights the increasing insecticide resistance of An. funestus mosquitoes, which are important malaria vectors in Africa, posing a significant challenge to malaria control efforts. While An. funestus has shown resistance to the recommended insecticide classes, notably pyrethroids and, in some cases, organochlorides and carbamates, it remains susceptible to other classes of insecticides such as organophosphates and neonicotinoids, providing potential alternative options for vector control strategies. The study underscores the need for targeted interventions that consider the population structure and geographical distribution of An. funestus, including its sibling species and their insecticide resistance profiles, to effectively achieve malaria elimination goals.

Copy number variants underlie the major selective sweeps in insecticide resistance genes in Anopheles arabiensis from Tanzania

To keep ahead of the evolution of resistance to insecticides in mosquitoes, national malaria control programmes must make use of a range of insecticides, both old and new, while monitoring resistance mechanisms. Knowledge of the mechanisms of resistance remains limited in Anopheles arabiensis, which in many parts of Africa is of increasing importance because it is apparently less susceptible to many indoor control interventions. Furthermore, comparatively little is known in general about resistance to non-pyrethroid insecticides such as pirimiphos-methyl (PM), which are crucial for effective control in the context of resistance to pyrethroids. We performed a genome-wide association study to determine the molecular mechanisms of resistance to deltamethrin (commonly used in bednets) and PM, in An. arabiensis from two regions in Tanzania. Genomic regions of positive selection in these populations were largely driven by copy number variants (CNVs) in gene families involved in resistance to these two insecticides. We found evidence of a new gene cluster involved in resistance to PM, identifying a strong selective sweep tied to a CNV in the Coeae2g-Coeae6g cluster of carboxylesterase genes. Using complementary data from An. coluzzii in Ghana, we show that copy number at this locus is significantly associated with PM resistance. Similarly, for deltamethrin, resistance was strongly associated with a novel CNV allele in the Cyp6aa / Cyp6p cluster. Against this background of metabolic resistance, target site resistance was very rare or absent for both insecticides. Mutations in the pyrethroid target site Vgsc were at very low frequency in Tanzania, yet combining these samples with three An. arabiensis individuals from West Africa revealed a startling diversity of evolutionary origins of target site resistance, with up to 5 independent origins of Vgsc-995 mutations found within just 8 haplotypes. Thus, despite having been first recorded over 10 years ago, Vgsc resistance mutations in Tanzanian An. arabiensis have remained at stable low frequencies. Overall, our results provide a new copy number marker for monitoring resistance to PM in malaria mosquitoes, and reveal the complex picture of resistance patterns in An. arabiensis.

The role of P450 enzymes in malaria and other vector-borne infectious diseases

Vector-borne infectious diseases are still an important global health problem. Malaria is the most important among them, mainly pediatric, life-threatening disease. Malaria and other vector-borne disorders caused by parasites, bacteria, and viruses have a strong impact on public health and significant economic costs. Most vector-borne diseases could be prevented by vector control, with attention to the ecological and biodiversity conservation aspects. Chemical control with pesticides and insecticides is widely used as a measure of prevention although increasing resistance to insecticides is a serious issue in vector control. Metabolic resistance is the most common mechanism and poses a big challenge. Insect enzyme systems, including monooxygenase CYP P450 enzymes, are employed by vectors mainly to metabolize insecticides thus causing resistance. The discovery and application of natural specific inhibitors/blockers of vector P450 enzymes as synergists for commonly used pesticides will contribute to the "greening" of insecticides. Besides vector CYPs, host CYP enzymes could also be exploited to fight against vector-borne diseases: using mostly their detoxifying properties and involvement in the immune response. Here, we review published research data on P450 enzymes from all players in vector-borne infections, that is, pathogens, vectors, and hosts, regarding the potential role of CYPs in disease. We discuss strategies on how to exploit cytochromes P450 in vector-borne disease control.

Deltamethrin and transfluthrin select for distinct transcriptomic responses in the malaria vector Anopheles gambiae

BACKGROUND

The widespread use of pyrethroid insecticides in Africa has led to the development of strong resistance in Anopheles mosquitoes. Introducing new active ingredients can contribute to overcome this phenomenon and ensure the effectiveness of vector control strategies. Transfluthrin is a polyfluorinated pyrethroid whose structural conformation was thought to prevent its metabolism by cytochrome P450 monooxygenases in malaria vectors, thus representing a potential alternative for managing P450-mediated resistance occurring in the field. In this study, a controlled selection was used to compare the dynamics of resistance between transfluthrin and the widely used pyrethroid deltamethrin in the mosquito Anopheles gambiae. Then, the associated molecular mechanisms were investigated using target-site mutation genotyping and RNA-seq.

METHODS

A field-derived line of An. gambiae carrying resistance alleles at low frequencies was used as starting material for a controlled selection experiment. Adult females were selected across 33 generations with deltamethrin or transfluthrin, resulting in three distinct lines: the Delta-R line (selected with deltamethrin), the Transflu-R line (selected with transfluthrin) and the Tiassale-S line (maintained without selection). Deltamethrin and transfluthrin resistance levels were monitored in each selected line throughout the selection process, as well as the frequency of the L1014F kdr mutation. At generation 17, cross-resistance to other public health insecticides was investigated and transcriptomes were sequenced to compare gene transcription variations and polymorphisms associated with adaptation to each insecticide.

RESULTS

A rapid increase in resistance to deltamethrin and transfluthrin was observed throughout the selection process in each selected line in association with an increased frequency of the L1014F kdr mutation. Transcriptomic data support a broader response to transfluthrin selection as compared to deltamethrin selection. For instance, multiple detoxification enzymes and cuticle proteins were specifically over-transcribed in the Transflu-R line including the known pyrethroid metabolizers CYP6M2, CYP9K1 and CYP6AA1 together with other genes previously associated with resistance in An. gambiae.

CONCLUSION

This study confirms that recurrent exposure of adult mosquitoes to pyrethroids in a public health context can rapidly select for various resistance mechanisms. In particular, it indicates that in addition to target site mutations, the polyfluorinated pyrethroid transfluthrin can select for a broad metabolic response, which includes some P450s previously associated to resistance to classical pyrethroids. This unexpected finding highlights the need for an in-depth study on the adaptive response of mosquitoes to newly introduced active ingredients in order to effectively guide and support decision-making programmes in malaria control.

Resilience of transfluthrin to oxidative attack by duplicated CYP6P9 variants known to confer pyrethroid resistance in the major malaria mosquito Anopheles funestus

Resistance to common pyrethroids, such as deltamethrin and permethrin is widespread in the malaria mosquito Anopheles funestus and mainly conferred by upregulated cytochrome P450 monooxygenases (P450s). In the pyrethroid resistant laboratory strain An. funestus FUMOZ-R the duplicated genes CYP6P9a and CYP6P9b are highly upregulated and have been shown to metabolize various pyrethroids, including deltamethrin and permethrin. Here, we recombinantly expressed CYP6P9a and CYP6P9b from An. funestus using a baculovirus expression system and evaluated the interaction of the multifluorinated benzyl pyrethroid transfluthrin with these enzymes by different approaches. First, by Michaelis-Menten kinetics in a fluorescent probe assay with the model substrate 7-benzyloxymethoxy-4-trifluoromethylcoumarin (BOMFC), we showed the inhibition of BOMFC metabolism by increasing concentrations of transfluthrin. Second, we tested the metabolic capacity of recombinantly expressed CYP6P9 variants to degrade transfluthrin utilizing UPLC-MS/MS analysis and detected low depletion rates, explaining the virtual lack of resistance of strain FUMOZ-R to transfluthrin observed in previous studies. However, as both approaches suggested an interaction of CYP6P9 variants with transfluthrin, we analyzed the oxidative metabolic fate and failed to detect hydroxylated transfluthrin, but low amounts of an M-2 transfluthrin metabolite. Based on the detected metabolite we hypothesize oxidative attack of the gem-dimethyl substituted cyclopropyl moiety, resulting in the formation of an allyl cation upon ring opening. In conclusion, these findings support the resilience of transfluthrin to P450-mediated pyrethroid resistance, and thus, reinforces its employment as an important resistance-breaking pyrethroid in resistance management strategies to control the major malaria vector An. funestus.

The duplicated P450s CYP6P9a/b drive carbamates and pyrethroids cross-resistance in the major African malaria vector Anopheles funestus

Cross-resistance to insecticides in multiple resistant malaria vectors is hampering resistance management. Understanding its underlying molecular basis is critical to implementation of suitable insecticide-based interventions. Here, we established that the tandemly duplicated cytochrome P450s, CYP6P9a/b are driving carbamate and pyrethroid cross-resistance in Southern African populations of the major malaria vector Anopheles funestus. Transcriptome sequencing revealed that cytochrome P450s are the most over-expressed genes in bendiocarb and permethrin-resistant An. funestus. The CYP6P9a and CYP6P9b genes are overexpressed in resistant An. funestus from Southern Africa (Malawi) versus susceptible An. funestus (Fold change (FC) is 53.4 and 17 respectively), while the CYP6P4a and CYP6P4b genes are overexpressed in resistant An. funestus in Ghana, West Africa, (FC is 41.1 and 17.2 respectively). Other up-regulated genes in resistant An. funestus include several additional cytochrome P450s (e.g. CYP9J5, CYP6P2, CYP6P5), glutathione-S transferases, ATP-binding cassette transporters, digestive enzymes, microRNA and transcription factors (FC<7). Targeted enrichment sequencing strongly linked a known major pyrethroid resistance locus (rp1) to carbamate resistance centering around CYP6P9a/b. In bendiocarb resistant An. funestus, this locus exhibits a reduced nucleotide diversity, significant p-values when comparing allele frequencies, and the most non-synonymous substitutions. Recombinant enzyme metabolism assays showed that both CYP6P9a/b metabolize carbamates. Transgenic expression of CYP6P9a/b in Drosophila melanogaster revealed that flies expressing both genes were significantly more resistant to carbamates than controls. Furthermore, a strong correlation was observed between carbamate resistance and CYP6P9a genotypes with homozygote resistant An. funestus (CYP6P9a and the 6.5kb enhancer structural variant) exhibiting a greater ability to withstand bendiocarb/propoxur exposure than homozygote CYP6P9a_susceptible (e.g Odds ratio = 20.8, P<0.0001 for bendiocarb) and heterozygotes (OR = 9.7, P<0.0001). Double homozygote resistant genotype (RR/RR) were even more able to survive than any other genotype combination showing an additive effect. This study highlights the risk that pyrethroid resistance escalation poses to the efficacy of other classes of insecticides. Available metabolic resistance DNA-based diagnostic assays should be used by control programs to monitor cross-resistance between insecticides before implementing new interventions.

In Vitro and In Vivo Validation of CYP6A14 and CYP6N6 Participation in Deltamethrin Metabolic Resistance in Aedes albopictus

The extensive use of chemical insecticides for public health and agricultural purposes has increased the occurrence and development of insecticide resistance. This study used transcriptome sequencing to screen 10 upregulated metabolic detoxification enzyme genes from Aedes albopictus resistant strains. Of these, CYP6A14 and CYP6N6 were found to be substantially overexpressed in the deltamethrin-induced expression test, indicating their role in deltamethrin resistance in Ae. albopictus. Furthermore, the corresponding 60-kDa recombinant proteins, CYP6A14 and CYP6N6, were successfully expressed using the Escherichia coli expression system. Enzyme activity studies revealed that CYP6A14 (5.84 U/L) and CYP6N6 (6.3 U/L) have cytochrome P450 (CYP450) enzyme activity. In vitro, the metabolic analysis revealed that the recombinant proteins degraded deltamethrin into 1-oleoyl-sn-glycero-3-phosphoethanolamine and 2',2'-dibromo-2'-deoxyguanosine. Subsequently, the CYP450 genes in larvae of Ae. albopictus were silenced by RNA interference technology to study deltamethrin resistance in vivo. The silencing of CYP6A14 and CYP6N6 increased the mortality rate of mosquitoes without affecting their survival time, spawning quantity, hatching rate, and other normal life activities. Altogether, CYP6A14 and CYP6N6 belong to the CYP6 family and mutually increase deltamethrin resistance in Ae. albopictus.

Detection of a reduced susceptibility to chlorfenapyr in the malaria vector Anopheles gambiae contrasts with full susceptibility in Anopheles funestus across Africa

New insecticides have recently been produced to help control pyrethroid-resistant malaria vectors including the pyrrole, chlorfenapyr. Monitoring the susceptibility of mosquito populations against this new product and potential cross-resistance with current insecticides is vital for better resistance management. In this study, we assessed the resistance status of the major malaria vectors Anopheles gambiae and Anopheles funestus to chlorfenapyr across Africa and explored potential cross-resistance with known pyrethroid resistance markers. Efficacy of chlorfenapyr 100 µg/ml against An. gambiae and An. funestus from five Cameroonian locations, the Democratic Republic of Congo, Ghana, Uganda, and Malawi was assessed using CDC bottle assays. Synergist assays were performed with PBO (4%), DEM (8%) and DEF (0.25%) and several pyrethroid-resistant markers were genotyped in both species to assess potential cross-resistance between pyrethroids and chlorfenapyr. Resistance to chlorfenapyr was detected in An. gambiae populations from DRC (Kinshasa) (mortality rate: 64.3 ± 7.1%) Ghana (Obuasi) (65.9 ± 7.4%), Cameroon (Mangoum; 75.2 ± 7.7% and Nkolondom; 86.1 ± 7.4). In contrast, all An. funestus populations were fully susceptible. A negative association was observed between the L1014F-kdr mutation and chlorfenapyr resistance with a greater frequency of homozygote resistant mosquitoes among the dead mosquitoes after exposure compared to alive (OR 0.5; P = 0.02) whereas no association was found between GSTe2 (I114T in An. gambiae; L119F in An. funestus) and resistance to chlorfenapyr. A significant increase of mortality to chlorfenapyr 10 µg/ml was observed in An. funestus after to PBO, DEM and DEF whereas a trend for a decreased mortality was observed in An. gambiae after PBO pre-exposure. This study reveals a greater risk of chlorfenapyr resistance in An. gambiae populations than in An. funestus. However, the higher susceptibility in kdr-resistant mosquitoes points to higher efficacy of chlorfenapyr against the widespread kdr-based pyrethroid resistance.

Sequential phase I metabolism of pyrethroids by duplicated CYP6P9 variants results in the loss of the terminal benzene moiety and determines resistance in the malaria mosquito Anopheles funestus

Pyrethroid resistance in Anopheles funestus is threatening the eradication of malaria. One of the major drivers of pyrethroid resistance in An. funestus are cytochrome P450 monooxygenases CYP6P9a and CYP6P9b, which are found upregulated in resistant An. funestus populations from Sub-Saharan Africa and are known to metabolise pyrethroids. Here, we have functionally expressed CYP6P9a and CYP6P9b variants and investigated their interactions with azole-fungicides and pyrethroids. Some azole fungicides such as prochloraz inhibited CYP6P9a and CYP6P9b at nanomolar concentrations, whereas pyrethroids were weak inhibitors (>100 μM). Amino acid sequence comparisons suggested that a valine to isoleucine substitution at position 310 in the active site cavity of CYP6P9a and CYP6P9b, respectively, might affect substrate binding and metabolism. We therefore swapped the residues by site directed mutagenesis to produce CYP6P9a^I310V and CYP6P9b^V310I. CYP6P9b^V310I produced stronger metabolic activity towards coumarin substrates and pyrethroids, particularly permethrin. The V310I mutation was previously also detected in a pyrethroid resistant field population of An. funestus in Benin. Additionally, we found the first metabolite of permethrin and deltamethrin after hydroxylation, 4'OH permethrin and 4'OH deltamethrin, were also suitable substrates for CYP6P9-variants, and were depleted by both enzymes to a higher extent than as their respective parent compounds (approximately 20% more active). Further, we found that both metabolites were toxic against An. funestus FANG (pyrethroid susceptible) but not towards FUMOZ-R (pyrethroid resistant) mosquitoes, the latter suggesting detoxification by overexpressed CYP6P9a and CYP6P9b. We confirmed by mass-spectrometric analysis that CYP6P9a and CYP6P9b are capable of cleaving phenoxybenzyl-ethers in type I pyrethroid permethrin and type II pyrethroid deltamethrin and that both enzymes preferentially metabolise trans-permethrin. This provides new insight into the metabolism of pyrethroids and a greater understanding of the molecular mechanisms of pyrethroid resistance in An. funestus.

Identification of a rapidly-spreading triple mutant for high-level metabolic insecticide resistance in Anopheles gambiae provides a real-time molecular diagnostic for antimalarial intervention deployment

Studies of insecticide resistance provide insights into the capacity of populations to show rapid evolutionary responses to contemporary selection. Malaria control remains heavily dependent on pyrethroid insecticides, primarily in long lasting insecticidal nets (LLINs). Resistance in the major malaria vectors has increased in concert with the expansion of LLIN distributions. Identifying genetic mechanisms underlying high-level resistance is crucial for the development and deployment of resistance-breaking tools. Using the Anopheles gambiae 1000 genomes (Ag1000g) data we identified a very recent selective sweep in mosquitoes from Uganda which localized to a cluster of cytochrome P450 genes. Further interrogation revealed a haplotype involving a trio of mutations, a nonsynonymous point mutation in Cyp6p4 (I236M), an upstream insertion of a partial Zanzibar-like transposable element (TE) and a duplication of the Cyp6aa1 gene. The mutations appear to have originated recently in An. gambiae from the Kenya-Uganda border, with stepwise replacement of the double-mutant (Zanzibar-like TE and Cyp6p4-236 M) with the triple-mutant haplotype (including Cyp6aa1 duplication), which has spread into the Democratic Republic of Congo and Tanzania. The triple-mutant haplotype is strongly associated with increased expression of genes able to metabolize pyrethroids and is strongly predictive of resistance to pyrethroids most notably deltamethrin. Importantly, there was increased mortality in mosquitoes carrying the triple-mutation when exposed to nets cotreated with the synergist piperonyl butoxide (PBO). Frequencies of the triple-mutant haplotype remain spatially variable within countries, suggesting an effective marker system to guide deployment decisions for limited supplies of PBO-pyrethroid cotreated LLINs across African countries.

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.

Pyrethroid resistance in the New World malaria vector Anopheles albimanus is mediated by cytochrome P450 CYP6P5

Pyrethroid resistance in the malaria vector Anopheles albimanus presents an obstacle to malaria elimination in the Americas. Here, An. albimanus CYP6P5 (the most overexpressed P450 in a Peruvian population) was functionally characterized. Recombinant CYP6P5 metabolized the type II pyrethroids, deltamethrin and α-cypermethrin with comparable affinities (K_M of 3.3 μM ± 0.4 and 3.6 μM ± 0.5, respectively), but exhibited a 2.7-fold higher catalytic rate for α-cypermethrin (k_cat of 6.02 min^-1 ± 0.2) versus deltamethrin (2.68 min^-1 ± 0.09). Time-course assays revealed progressive depletion of the above pyrethroids with production of four HPLC-detectable metabolites. Low depletion was obtained with type I pyrethroid, permethrin. Transgenic expression in Drosophila melanogaster demonstrated that overexpression of CYP6P5 alone conferred type II pyrethroid resistance, with only 16% and 55.3% mortalities in flies exposed to 0.25% α-cypermethrin and 0.15% deltamethrin, respectively. Synergist bioassays using P450 inhibitor piperonylbutoxide significantly recovered susceptibility (mortality = 73.6%, p < 0.001) in synergized flies exposed to 4% piperonylbutoxide, plus 0.25% α-cypermethrin, compared to non-synergized flies (mortality = 4.9%). Moderate resistance was also observed towards 4% DDT. These findings established the preeminent role of CYP6P5 in metabolic resistance in An. albimanus, highlighting challenges associated with deployment of insecticide-based control tools in the Americas.

Biochemical profiling of functionally expressed CYP6P9 variants of the malaria vector Anopheles funestus with special reference to cytochrome b5 and its role in pyrethroid and coumarin substrate metabolism

Cytochrome P450 monooxygenases (P450s) are well studied enzymes catalyzing the oxidative metabolism of xenobiotics in insects including mosquitoes. Their duplication and upregulation in agricultural and public health pests such as anopheline mosquitoes often leads to an enhanced metabolism of insecticides which confers resistance. In the laboratory strain Anopheles funestus FUMOZ-R the duplicated P450s CYP6P9a and CYP6P9b are highly upregulated and proven to confer pyrethroid resistance. Microsomal P450 activity is regulated by NADPH cytochrome P450 oxidoreductase (CPR) required for electron transfer, whereas the modulatory role of cytochrome b₅ (CYB5) on insect P450 activity is less clear. In previous studies CYP6P9a and CYP6P9b were recombinantly expressed in tandem with An. gambiae CPR using E. coli-expression systems and CYB5 added to the reaction mix to enhance activity. However, the precise role of CYB5 on substrate turn-over when combined with CYP6P9a and CYP6P9b remains poorly investigated, thus one objective of our study was to address this knowledge gap. In contrast to the CYP6P9 variants, the expression levels of both CYB5 and CPR were not upregulated in the pyrethroid resistant FUMOZ-R strain when compared to the susceptible FANG strain, suggesting no immediate regulatory role of these genes in pyrethroid resistance in FUMOZ-R. Here, for the first time we recombinantly expressed CYP6P9a and CYP6P9b from An. funestus in a baculovirus expression system using High-5 insect cells. Co-expression of each enzyme with CPR from either An. gambiae or An. funestus did not reveal noteworthy differences in catalytic capacity. Whereas the co-expression of An. funestus CYB5 - tested at different multiplicity of infection (MOI) ratios - resulted in a significantly higher metabolization of coumarin substrates as measured by fluorescence assays. This was confirmed by Michaelis-Menten kinetics using the most active substrate, 7-benzyloxymethoxy-4-trifluoromethylcoumarin (BOMFC). We observed a similar increase in coumarin substrate turnover by adding human CYB5 to the reaction mix. Finally, we compared by UPLC-MS/MS analysis the depletion rate of deltamethrin and the formation of 4'OH-deltamethrin by recombinantly expressed CYP6P9a and CYP6P9b with and without CYB5 and detected no difference in the extent of deltamethrin metabolism. Our results suggest that co-expression (or addition) of CYB5 with CYP6P9 variants, recombinantly expressed in insect cells, can significantly enhance their metabolic capacity to oxidize coumarins, but not deltamethrin.

Functionally characterized arthropod pest and pollinator cytochrome P450s associated with xenobiotic metabolism

The cytochrome P450 family (P450s) of arthropods includes diverse enzymes involved in endogenous essential physiological functions and in the oxidative metabolism of xenobiotics, insecticides and plant allelochemicals. P450s can also establish insecticide selectivity in bees and pollinators. Several arthropod P450s, distributed in different phylogenetic groups, have been associated with xenobiotic metabolism, and some of them have been functionally characterized, using different in vitro and in vivo systems. The purpose of this review is to summarize scientific publications on arthropod P450s from major insect and mite agricultural pests, pollinators and Papilio sp, which have been functionally characterized and shown to metabolize xenobiotics and/or their role (direct or indirect) in pesticide toxicity or resistance has been functionally validated. The phylogenetic relationships among these P450s, the functional systems employed for their characterization and their xenobiotic catalytic properties are presented, in a systematic approach, including critical aspects and limitations. The potential of the primary P450-based metabolic pathway of target and non-target organisms for the development of highly selective insecticides and resistance-breaking formulations may help to improve the efficiency and sustainability of pest control.

The cytochrome P450 CYP325A is a major driver of pyrethroid resistance in the major malaria vector Anopheles funestus in Central Africa

The overexpression and overactivity of key cytochrome P450s (CYP450) genes are major drivers of metabolic resistance to insecticides in African malaria vectors such as Anopheles funestus s.s. Previous RNAseq-based transcription analyses revealed elevated expression of CYP325A specific to Central African populations but its role in conferring resistance has not previously been demonstrated. In this study, RT-qPCR consistently confirmed that CYP325A is highly over-expressed in pyrethroid-resistant An. funestus from Cameroon, compared with a control strain and insecticide-unexposed mosquitoes. A synergist bioassay with PBO significantly recovered susceptibility for permethrin and deltamethrin indicating P450-based metabolic resistance. Analyses of the coding sequence of CYP325A Africa-wide detected high-levels of polymorphism, but with no predominant alleles selected by pyrethroid resistance. Geographical amino acid changes were detected notably in Cameroon. In silico homology modelling and molecular docking simulations predicted that CYP325A binds and metabolises type I and type II pyrethroids. Heterologous expression of recombinant CYP325A and metabolic assays confirmed that the most-common Cameroonian haplotype metabolises both type I and type II pyrethroids with depletion rate twice that the of the DR Congo haplotype. Analysis of the 1 kb putative promoter of CYP325A revealed reduced diversity in resistant mosquitoes compared to susceptible ones, suggesting a potential selective sweep in this region. The establishment of CYP325A as a pyrethroid resistance metabolising gene further explains pyrethroid resistance in Central African populations of An. funestus. Our work will facilitate future efforts to detect the causative resistance markers in the promoter region of CYP325A to design field applicable DNA-based diagnostic tools.

Assessing cross-resistance within the pyrethroids in terms of their interactions with key cytochrome P450 enzymes and resistance in vector populations

Background

It is important to understand whether the potential impact of pyrethroid resistance on malaria control can be mitigated by switching between different pyrethroids or whether cross-resistance within this insecticide class precludes this approach.

Methods

Here we assess the relationships among pyrethroids in terms of their binding affinity to, and depletion by, key cytochrome P450 enzymes (hereafter P450s) that are known to confer metabolic pyrethroid resistance in Anopheles gambiae (s.l.) and An. funestus, in order to identify which pyrethroids may diverge from the others in their vulnerability to resistance. We then investigate whether these same pyrethroids also diverge from the others in terms of resistance in vector populations.

Results

We found that the type I and II pyrethroids permethrin and deltamethrin, respectively, are closely related in terms of binding affinity to key P450s, depletion by P450s and resistance within vector populations. Bifenthrin, which lacks the common structural moiety of most pyrethroids, diverged from the other pyrethroids tested in terms of both binding affinity to key P450s and depletion by P450s, but resistance to bifenthrin has rarely been tested in vector populations and was not analysed here. Etofenprox, which also lacks the common structural moiety of most pyrethroids, diverged from the more commonly deployed pyrethroids in terms of binding affinity to key P450s and resistance in vector populations, but did not diverge from these pyrethroids in terms of depletion by the P450s. The analysis of depletion by the P450s indicated that etofenprox may be more vulnerable to metabolic resistance mechanisms in vector populations. In addition, greater resistance to etofenprox was found across Aedes aegypti populations, but greater resistance to this compound was not found in any of the malaria vector species analysed. The results for pyrethroid depletion by anopheline P450s in the laboratory were largely not repeated in the findings for resistance in malaria vector populations.

Conclusion

Importantly, the prevalence of resistance to the pyrethroids α-cypermethrin, cyfluthrin, deltamethrin, λ-cyhalothrin and permethrin was correlated across malaria vector populations, and switching between these compounds as a tool to mitigate against pyrethroid resistance is not advised without strong evidence supporting a true difference in resistance.

Cytochrome P450-based metabolic insecticide resistance in Anopheles and Aedes mosquito vectors: Muddying the waters

Despite the substantial progress achieved in the characterization of cytochrome P450 (CYP) -based resistance mechanisms in mosquitoes, a number of questions remain unanswered. These include: (i) the regulation and physiology of resistance conferring CYPs; (ii) the actual contribution of CYPs in resistance alone or in combination with other detoxification partners or other resistance mechanisms; (iii) the association between overexpression levels and allelic variation, with the catalytic activity and the intensity of resistance and (iv) the true value of molecular diagnostics targeting CYP markers, for driving decision making in the frame of Insecticide Resistance Management applications. Furthermore, the translation of CYP - based insecticide resistance research in mosquitoes into practical applications, is being developed, but it is not fully exploited, as yet. Examples include the production of high throughput platforms for screening the liability (stability) or inhibition potential of novel insecticidal leads and synergists (add-ons), as well as the exploration of the negative cross resistance concept (i.e. detoxification of certain insecticides, but activation of others pro-insecticides). The goal of this review is to critically summarise the current knowledge and the gaps of the CYP-based metabolic insecticide resistance in Anopheles and Aedes mosquito vectors. The progress and limitations of the protein and the reverse/forward genetic approaches, the understanding and importance of molecular and physiological aspects, as well as the current and future exploitation routes of CYP research are discussed.

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.

Malaria control currently relies heavily on insecticide-based vector control interventions. Unfortunately, resistance to insecticides is threatening their continued effectiveness. Metabolic resistance has the greatest operational significance, yet it remains unclear how mosquito populations evolutionarily respond to the massive selection pressure from control interventions including insecticide-treated nets. Deciphering patterns of gene flow between populations of major malaria vectors such as Anopheles funestus and elucidating genomic signature of resistance evolution are crucial for designing resistance management strategies and preventing malaria resurgence. Here, we performed a genome-wide survey of An. funestus genetic diversity from across its continental range using reduced-genome representation (ddRADseq) and whole genome (PoolSeq) approaches revealing evidence of significant barriers to gene flow impacting the spread of insecticide resistance alleles. This study detected signatures of strong selective sweeps occurring in genomic regions controlling cytochrome P450-based and glutathione s-transferase metabolic resistance to insecticides in this species. Fine-scale analysis of the major pyrethroid resistance-associated genomic regions revealed complex molecular evolution with evidence of copy number variation, transposon insertion and gene conversion highlighting the risk that if this level of selection and spread of resistance continues unabated, our ability to control malaria with current interventions will be compromised.

Minimization of energy transduction confers resistance to phosphine in the rice weevil, Sitophilus oryzae

Infestation of phosphine (PH₃) resistant insects threatens global grain reserves. PH₃ fumigation controls rice weevil (Sitophilus oryzae) but not highly resistant insect pests. Here, we investigated naturally occurring strains of S. oryzae that were moderately resistant (MR), strongly resistant (SR), or susceptible (wild-type; WT) to PH₃ using global proteome analysis and mitochondrial DNA sequencing. Both PH₃ resistant (PH₃-R) strains exhibited higher susceptibility to ethyl formate-mediated inhibition of cytochrome c oxidase than the WT strain, whereas the disinfectant PH₃ concentration time of the SR strain was much longer than that of the MR strain. Unlike the MR strain, which showed altered expression levels of genes encoding metabolic enzymes involved in catabolic pathways that minimize metabolic burden, the SR strain showed changes in the mitochondrial respiratory chain. Our results suggest that the acquisition of strong PH₃ resistance necessitates the avoidance of oxidative phosphorylation through the accumulation of a few non-synonymous mutations in mitochondrial genes encoding complex I subunits as well as nuclear genes encoding dihydrolipoamide dehydrogenase, concomitant with metabolic reprogramming, a recognized hallmark of cancer metabolism. Taken together, our data suggest that reprogrammed metabolism represents a survival strategy of SR insect pests for the compensation of minimized energy transduction under anoxic conditions. Therefore, understanding the resistance mechanism of PH₃-R strains will support the development of new strategies to control insect pests.

Whole-genome sequencing reveals high complexity of copy number variation at insecticide resistance loci in malaria mosquitoes

Polymorphisms in genetic copy number can influence gene expression, coding sequence, and zygosity, making them powerful actors in the evolutionary process. Copy number variants (CNVs) are however understudied, being more difficult to detect than single-nucleotide polymorphisms. We take advantage of the intense selective pressures on the major malaria vector Anopheles gambiae, caused by the widespread use of insecticides for malaria control, to investigate the role of CNVs in the evolution of insecticide resistance. Using the whole-genome sequencing data from 1142 samples in the An. gambiae 1000 genomes project, we identified 250 gene-containing CNVs, encompassing a total of 267 genes of which 28 were in gene families linked to metabolic insecticide resistance, representing significant enrichment of these families. The five major gene clusters for metabolic resistance all contained CNVs, with 44 different CNVs being found across these clusters and multiple CNVs frequently covering the same genes. These 44 CNVs are widespread (45% of individuals carry at least one of them) and have been spreading through positive selection, indicated by their high local frequencies and extended haplotype homozygosity. Our results demonstrate the importance of CNVs in the response to selection, highlighting the urgent need to identify the contribution of each CNV to insecticide resistance and to track their spread as the use of insecticides in malaria endemic countries intensifies and as the operational deployment of next-generation bed nets targeting metabolic resistance gathers pace. Our detailed descriptions of CNVs found across the species range provide the tools to do so.

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.