RNAseq-based gene expression profiling of the Anopheles funestus pyrethroid-resistant strain FUMOZ highlights the predominant role of the duplicated CYP6P9a/b cytochrome P450s

Metabolic resistance to pyrethroids with possible involvement of non-coding ribonucleic acids in Anopheles funestus, the major malaria vector in western Kenya

BACKGROUND

The resurgence of Anopheles funestus, a dominant vector of human malaria in western Kenya was partly attributed to insecticide resistance. However, evidence on the molecular basis of pyrethroid resistance in western Kenya is limited. Here, we reported metabolic resistance mechanisms and demonstrated that multiple non-coding Ribonucleic Acids (ncRNAs) could play a potential role in An. funestus resistance to pyrethroid in western Kenya. Anopheles funestus mosquitoes were sampled using aspiration methods in Bungoma, Teso, Siaya, Port Victoria and Kombewa in western Kenya. The F1 progenies were exposed to deltamethrin (0.05%), permethrin (0.75%), DDT (4%) and pirimiphos-methyl (0.25%) following WHO test guidelines. A synergist assay using piperonyl butoxide (PBO) (4%) was conducted to determine cytochrome P450s' role in pyrethroid resistance. RNA-seq was conducted on a combined pool of specimens that were resistant and unexposed, and the results were compared with those of the FANG susceptible reference strain. This approach aimed to uncover the molecular mechanisms underlying the observed phenotypic pyrethroid resistance.

RESULTS

Pyrethroid resistance was observed in all sites with an average mortality rate (MR) of 57.6%. Port Victoria had the highest level of pyrethroid resistance to permethrin (MR = 53%) and deltamethrin (MR = 11%. Teso had the lowest level of resistance to permethrin (MR = 70%) and deltamethrin (MR = 87%). Resistance to DDT was observed only in Kombewa (MR = 89%) and Port Victoria (MR = 85%). A full susceptibility to P-methyl (0.25%) was observed in all sites. PBO synergist assay revealed high susceptibility (> 98%) to pyrethroids in all the sites except for Port Victoria (MR = 96%). Whole transcriptomic analysis showed that most gene families associated with pyrethroid resistance comprised non-coding RNAs (67%), followed by immunity proteins (10%), cytochrome P450s (6%), cuticular proteins (5%), olfactory proteins (4%), glutathione S-transferases (3%), UDP-glycosyltransferases (2%), ATP-binding cassettes (2%) and carboxylesterases (1%).

CONCLUSION

This study unveils the molecular basis of insecticide resistance in An. funestus in western Kenya, highlighting for the first time the potential role of non-coding RNAs alongside metabolic detoxification in pyrethroid resistance. Targeting non-coding RNAs for intervention development could help in insecticide resistance management.

Temporal evolution of insecticide resistance and bionomics in Anopheles funestus, a key malaria vector in Uganda

Insecticide resistance escalation is decreasing the efficacy of vector control tools. Monitoring vector resistance is paramount in order to understand its evolution and devise effective counter-solutions. In this study, we monitored insecticide resistance patterns, vector population bionomics and genetic variants associated with resistance over 3 years from 2021 to 2023 in Uganda. Anopheles funestus s.s was the predominant species in Mayuge but with evidence of hybridization with other species of the An. funestus group. Sporozoite infection rates were relatively very high with a peak of 20.41% in March 2022. Intense pyrethroid resistance was seen against pyrethroids up to 10-times the diagnostic concentration but partial recovery of susceptibility in PBO synergistic assays. Among bednets, only PBO-based nets (PermaNet 3.0 Top and Olyset Plus) and chlorfenapyr-based net (Interceptor G2) had high mortality rates. Mosquitoes were fully susceptible to chlorfenapyr and organophosphates, moderately resistant to clothianidin and resistant to carbamates. The allele frequency of key P450, CYP9K1, resistance marker was constantly very high but that for CYP6P9A/b were very low. Interestingly, we report the first detection of resistance alleles for Ace1 gene (RS =  ~ 13%) and Rdl gene (RS =  ~ 21%, RR =  ~ 4%) in Uganda. The qRT-PCR revealed that Cytochrome P450s CYP9K1, CYP6P9A, CYP6P9b, CYP6P5 and CYP6M7 were consistently upregulated while a glutathione-S-transferase gene (GSTE2) showed low expression. Our study shows the complexity of insecticide resistance patterns and underlying mechanisms, hence constant and consistent spatial and temporal monitoring is crucial to rapidly detect changing resistance profiles which is key in informing deployment of counter interventions.

Two highly selected mutations in the tandemly duplicated CYP6P4a and CYP6P4b genes drive pyrethroid resistance in Anopheles funestus in West Africa

BACKGROUND

Gaining a comprehensive understanding of the genetic mechanisms underlying insecticide resistance in malaria vectors is crucial for optimising the effectiveness of insecticide-based vector control methods and developing diagnostic tools for resistance management. Considering the heterogeneity of metabolic resistance in major malaria vectors, the implementation of tailored resistance management strategies is essential for successful vector control. Here, we provide evidence demonstrating that two highly selected mutations in CYP6P4a and CYP6P4b are driving pyrethroid insecticide resistance in the major malaria vector Anopheles funestus, in West Africa.

RESULTS

Continent-wide polymorphism survey revealed escalated signatures of directional selection of both genes between 2014 and 2021. In vitro insecticide metabolism assays with recombinant enzymes from both genes showed that mutant alleles under selection exhibit higher metabolic efficiency than their wild-type counterparts. Using the GAL4-UAS expression system, transgenic Drosophila flies overexpressing mutant alleles exhibited increased resistance to pyrethroids. These findings were consistent with in silico predictions which highlighted changes in enzyme active site architecture that enhance the affinity of mutant alleles for type I and II pyrethroids. Furthermore, we designed two DNA-based assays for the detection of CYP6P4a-M220I and CYP6P4b-D284E mutations, showing their current confinement to West Africa. Genotype/phenotype correlation analyses revealed that these markers are strongly associated with resistance to types I and II pyrethroids and combine to drastically reduce killing effects of pyrethroid bed nets.

CONCLUSIONS

Overall, this study demonstrated that CYP6P4a and CYP6P4b contribute to pyrethroid resistance in An. funestus and provided two additional insecticide resistance molecular diagnostic tools that would contribute to monitoring and better management of resistance.

The Dual Resistance Mechanism of CYP325G4 and CYP6AA9 in Culex pipiens pallens Legs According to Transcriptome and Proteome Analysis

Mosquitoes within the Culex pipiens complex play a crucial role in human disease transmission. Insecticides, especially pyrethroids, are used to control these vectors. Mosquito legs are the main entry point and barrier for insecticides to gain their neuronal targets. However, the resistance mechanism in mosquito legs is unclear. Herein, we employed transcriptomic analyses and isobaric tags for relative and absolute quantitation techniques to investigate the resistance mechanism, focusing on Cx. pipiens legs. We discovered 2346 differentially expressed genes (DEGs) between deltamethrin-resistant (DR) and deltamethrin-sensitive (DS) mosquito legs, including 41 cytochrome P450 genes. In the same comparison, we identified 228 differentially expressed proteins (DEPs), including six cytochrome P450 proteins. Combined transcriptome and proteome analysis revealed only two upregulated P450 genes, CYP325G4 and CYP6AA9. The main clusters of DEGs and DEPs were associated with metabolic processes, such as cytochrome P450-mediated metabolism of drugs and xenobiotics. Transcription analysis revealed high CYP325G4 and CYP6AA9 expression in the DR strain at 72 h posteclosion compared with that in the DS strain, particularly in the legs. Mosquitoes knocked down for CYP325G4 were more sensitive to deltamethrin than the controls. CYP325G4 knockdown reduced the expression of several chlorinated hydrocarbon (CHC)-related genes, which altered the cuticle thickness and structure. Conversely, CYP6AA9 knockdown increased CHC gene expression without altering cuticle thickness and structure. P450 activity analysis demonstrated that CYP325G4 and CYP6AA9 contributed to metabolic resistance in the midgut and legs. This study identified CYP325G4 as a novel mosquito deltamethrin resistance factor, being involved in both metabolic and cuticular resistance mechanisms. The previously identified CYP6AA9 was investigated for its involvement in metabolic resistance and potential cuticular resistance in mosquito legs. These findings enhance our comprehension of resistance mechanisms, identifying P450s as promising targets for the future management of mosquito vector resistance, and laying a theoretical groundwork for mosquito 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.

Integrated de novo transcriptome of Culex pipiens mosquito larvae as a resource for genetic control strategies

We present a de novo transcriptome of the mosquito vector Culex pipiens, assembled by sequences of susceptible and insecticide resistant larvae. The high quality of the assembly was confirmed by TransRate and BUSCO. A mapping percentage until 94.8% was obtained by aligning contigs to Nr, SwissProt, and TrEMBL, with 27,281 sequences that simultaneously mapped on the three databases. A total of 14,966 ORFs were also functionally annotated by using the eggNOG database. Among them, we identified ORF sequences of the main gene families involved in insecticide resistance. Therefore, this resource stands as a valuable reference for further studies of differential gene expression as well as to identify genes of interest for genetic-based control tools.

Overexpression and nonsynonymous mutations of UDP-glycosyltransferases are potentially associated with pyrethroid resistance in Anopheles funestus

UDP-glycosyltransferases (UGTs) enzymes are pivotal in insecticide resistance by transforming hydrophobic substrates into more hydrophilic forms for efficient cell elimination. This study provides the first comprehensive investigation of Anopheles funestus UGT genes, their evolution, and their association with pyrethroid resistance. We employed a genome-wide association study using pooled sequencing (GWAS-PoolSeq) and transcriptomics on pyrethroid-resistant An. funestus, along with deep-targeted sequencing of UGTs in 80 mosquitoes Africa-wide. UGT310B2 was consistently overexpressed Africa-wide and significant gene-wise Fst differentiation was observed between resistant and susceptible populations: UGT301C2 and UGT302A3 in Malawi, and UGT306C2 in Uganda. Additionally, nonsynonymous mutations in UGT genes were identified. Gene-wise Tajima's D density curves provide insights into population structures within populations across these countries, supporting previous observations. These findings have important implications for current An. funestus control strategies facilitating the prediction of cross-resistance to other UGT-metabolised polar insecticides, thereby guiding more effective and targeted insecticide resistance management efforts.

Identification of putative promoter elements for epsilon glutathione s-transferases genes associated with resistance to DDT in the malaria vector mosquito anopheles arabiensis

The purpose of this study was to identify the putative regulatory elements in the promoter region of An. arabiensis strains which differed in susceptibility to DDT and compare with those identified in its sibling An. gambaie. Basal expression level of Epsilon class GSTs (Glutathione S - transferases) GSTe1 gene was 0.512 - 0.658 (95% CI) and GSTe2 0.672 - 1.204 (95% CI) in adults of DDT resistant KGB compared to 0.031 - 0.04 (95% CI) and 0.148 - 0.199 (95% CI) respectively in susceptible MAT strains of An. arabiensis. Induced mean expression of GSTe2 in larvae exposed to DDT for one hour was 0.901 - 1.172 (95% CI) in KGB and 0.475 - 0.724 (95% CI) in MAT strain. In present work, strain specific primers were used to amplify and sequenced the promoter regions of GSTe1 and GSTe2 in the KGB, MAT and field specimens. Computational analysis revealed presence of classical arthropod initiator sequence TCAGT and putative core promoter elements, GC, CAAT, TATA boxes. A typical TATA box was identified at 35 bp upstream Transcription Start Site (TSS) in GSTe1 but was absent in GSTe2. Several binding sites for regulatory elements downstream and multiple polymorphic sites were identified between strains. The role of these regulatory elements in transcription of these genes has not been determined. However, on comparison the 2 bp adenosine indel (insertion/deletion) which was essential in driving the promoter activity in An. gambiae was identified only DDT resistant KGB strain.

RNA-Seq-Pop: Exploiting the sequence in RNA sequencing-A Snakemake workflow reveals patterns of insecticide resistance in the malaria vector Anopheles gambiae

We provide a reproducible and scalable Snakemake workflow, called RNA-Seq-Pop, which provides end-to-end analysis of RNA sequencing data sets. The workflow allows the user to perform quality control, perform differential expression analyses and call genomic variants. Additional options include the calculation of allele frequencies of variants of interest, summaries of genetic variation and population structure, and genome-wide selection scans, together with clear visualizations. RNA-Seq-Pop is applicable to any organism, and we demonstrate the utility of the workflow by investigating pyrethroid resistance in selected strains of the major malaria mosquito, Anopheles gambiae. The workflow provides additional modules specifically for An. gambiae, including estimating recent ancestry and determining the karyotype of common chromosomal inversions. The Busia laboratory colony used for selections was collected in Busia, Uganda, in November 2018. We performed a comparative analysis of three groups: a parental G24 Busia strain; its deltamethrin-selected G28 offspring; and the susceptible reference strain Kisumu. Measures of genetic diversity reveal patterns consistent with that of laboratory colonization and selection, with the parental Busia strain exhibiting the highest nucleotide diversity, followed by the selected Busia offspring, and finally, Kisumu. Differential expression and variant analyses reveal that the selected Busia colony exhibits a number of distinct mechanisms of pyrethroid resistance, including the Vgsc-995S target-site mutation, upregulation of SAP genes, P450s and a cluster of carboxylesterases. During deltamethrin selections, the 2La chromosomal inversion rose in frequency (from 33% to 86%), supporting a previous link with pyrethroid resistance. RNA-Seq-Pop is hosted at: github.com/sanjaynagi/rna-seq-pop. We anticipate that the workflow will provide a useful tool to facilitate reproducible, transcriptomic studies in An. gambiae and other taxa.

Detect Cytochrome C Oxidase- and Glutathione-S-Transferase-Mediated Detoxification in a Permethrin-Resistant Population of Lygus lineolaris

Frequent sprays on cotton prompted resistance development in the tarnished plant bug (TPB). Knowledge of global gene regulation is highly desirable to better understand resistance mechanisms and develop molecular tools for monitoring and managing resistance. Novel microarray expressions of 6688 genes showed 3080 significantly up- or down-regulated genes in permethrin-treated TPBs. Among the 1543 up-regulated genes, 255 code for 39 different enzymes, and 15 of these participate in important pathways and metabolic detoxification. Oxidase is the most abundant and over-expressed enzyme. Others included dehydrogenases, synthases, reductases, and transferases. Pathway analysis revealed several oxidative phosphorylations associated with 37 oxidases and 23 reductases. One glutathione-S-transferase (GST LL_2285) participated in three pathways, including drug and xenobiotics metabolisms and pesticide detoxification. Therefore, a novel resistance mechanism of over-expressions of oxidases, along with a GST gene, was revealed in permethrin-treated TPB. Reductases, dehydrogenases, and others may also indirectly contribute to permethrin detoxification, while two common detoxification enzymes, P450 and esterase, played less role in the degradation of permethrin since none was associated with the detoxification pathway. Another potential novel finding from this study and our previous studies confirmed multiple/cross resistances in the same TPB population with a particular set of genes for different insecticide classes.

Tekoh T, S. Ibrahim S, Muhammad A, Kouamo M, Wondji MJ, Irving H, Hearn J, Wondji CS. 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.

Marked Effects of Larval Salt Exposure on the Life History and Gut Microbiota of the Malaria Vector Anopheles merus (Diptera: Culicidae)

Anopheles merus can breed in a range of saltwater concentrations. The consequences of this ability on the life history of adult An. merus are poorly understood. This study examined the effects of exposure to 0, 2.1875, 4.375, 8.75, and 17.5 g/L of sodium chloride on An. merus. The effects on larval development, adult longevity, fertility, and fecundity, as well as deltamethrin tolerance were examined. The effect of larval salt exposure on the expression of defensin-1 in adults was examined by quantitative Real-Time PCR. Finally, the effect of the larval salt concentration on microbial dynamics was assessed by 16S Next Generation Sequencing. High concentrations of saltwater increased larval development time and number of eggs laid, as well as deltamethrin tolerance. Larval exposure to salt also reduced the expression of defensin-1. The exposure also had a significant effect on microbial diversity in larvae and adults. The diversity of larvae decreased once adults emerged. Salt-tolerant bacterial genera predominated in larvae but were absent in adults. High salt concentrations resulted in greater abundance of Plasmodium-protective genera in adults. Although this study was conducted on a laboratory strain of An. merus, these data suggest that osmoregulation has a significant effect on the life history of the species with potential epidemiological consequences.

Experimental Hut Trials Reveal That CYP6P9a/b P450 Alleles Are Reducing the Efficacy of Pyrethroid-Only Olyset Net against the Malaria Vector Anopheles funestus but PBO-Based Olyset Plus Net Remains Effective

Malaria remains a major public health concern in Africa. Metabolic resistance in major malaria vectors such as An. funestus is jeopardizing the effectiveness of long-lasting insecticidal nets (LLINs) to control malaria. Here, we used experimental hut trials (EHTs) to investigate the impact of cytochrome P450-based resistance on the efficacy of PBO-based net (Olyset Plus) compared to a permethrin-only net (Olyset), revealing a greater loss of efficacy for the latter. EHT performed with progenies of F5 crossing between the An. funestus pyrethroid-resistant strain FUMOZ and the pyrethroid-susceptible strain FANG revealed that PBO-based nets (Olyset Plus) induced a significantly higher mortality rate (99.1%) than pyrethroid-only nets (Olyset) (56.7%) (p < 0.0001). The blood-feeding rate was higher in Olyset compared to Olyset Plus (11.6% vs. 5.6%; p = 0.013). Genotyping the CYP6P9a/b and the intergenic 6.5 kb structural variant (SV) resistance alleles showed that, for both nets, homozygote-resistant mosquitoes have a greater ability to blood-feed than the susceptible mosquitoes. Homozygote-resistant genotypes significantly survived more with Olyset after cone assays (e.g., CYP6P9a OR = 34.6; p < 0.0001) than homozygote-susceptible mosquitoes. A similar but lower correlation was seen with Olyset Plus (OR = 6.4; p < 0.001). Genotyping EHT samples confirmed that CYP6P9a/b and 6.5 kb_SV homozygote-resistant mosquitoes survive and blood-feed significantly better than homozygote-susceptible mosquitoes when exposed to Olyset. Our findings highlight the negative impact of P450-based resistance on pyrethroid-only nets, further supporting that PBO nets, such as Olyset Plus, are a better solution in areas of P450-mediated resistance to pyrethroids.