Molecular drivers of insecticide resistance in the Sahelo-Sudanian populations of a major malaria vector Anopheles coluzzii.

BACKGROUND: Information on common markers of metabolic resistance in malaria vectors from countries sharing similar eco-climatic characteristics can facilitate coordination of malaria control. Here, we characterized populations of the major malaria vector Anopheles coluzzii from Sahel region, spanning four sub-Saharan African countries: Nigeria, Niger, Chad and Cameroon. RESULTS: Genome-wide transcriptional analysis identified major genes previously implicated in pyrethroid and/or cross-resistance to other insecticides, overexpressed across the Sahel, including CYP450s, glutathione S-transferases, carboxylesterases and cuticular proteins. Several, well-known markers of insecticide resistance were found in high frequencies-including in the voltage-gated sodium channel (V402L, I940T, L995F, I1527T and N1570Y), the acetylcholinesterase-1 gene (G280S) and the CYP4J5-L43F (which is fixed). High frequencies of the epidemiologically important chromosomal inversion polymorphisms, 2La, 2Rb and 2Rc, were observed (~80% for 2Rb and 2Rc). The 2La alternative arrangement is fixed across the Sahel. Low frequencies of these inversions (<10%) were observed in the fully insecticide susceptible laboratory colony of An. coluzzii (Ngoussou). Several of the most commonly overexpressed metabolic resistance genes sit in these three inversions. Two commonly overexpressed genes, GSTe2 and CYP6Z2, were functionally validated. Transgenic Drosophila melanogaster flies expressing GSTe2 exhibited extremely high DDT and permethrin resistance (mortalities <10% in 24h). Serial deletion of the 5' intergenic region, to identify putative nucleotide(s) associated with GSTe2 overexpression, revealed that simultaneous insertion of adenine nucleotide and a transition (T->C), between Forkhead box L1 and c-EST putative binding sites, were responsible for the high overexpression of GSTe2 in the resistant mosquitoes. Transgenic flies expressing CYP6Z2 exhibited marginal resistance towards 3-phenoxybenzylalcohol (a primary product of pyrethroid hydrolysis by carboxylesterases) and a type II pyrethroid, α-cypermethrin. However, significantly higher mortalities were observed in CYP6Z2 transgenic flies compared with controls, on exposure to the neonicotinoid, clothianidin. This suggests a possible bioactivation of clothianidin into a toxic intermediate, which may make it an ideal insecticide against populations of An. coluzzii overexpressing this P450. CONCLUSIONS: These findings will facilitate regional collaborations within the Sahel region and refine implementation strategies through re-focusing interventions, improving evidence-based, cross-border policies towards local and regional malaria pre-elimination.

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.

Exploring the molecular mechanisms of increased intensity of pyrethroid resistance in Central African population of a major malaria vector Anopheles coluzzii.

Molecular mechanisms driving the escalation of pyrethroid resistance in the major malaria mosquitoes of Central Africa remain largely uncharacterized, hindering effective management strategies. Here, resistance intensity and the molecular mechanisms driving it were investigated in a population of Anopheles coluzzii from northern Cameroon. High levels of pyrethroid and organochloride resistance were observed in An. coluzzii population, with no mortality for 1× permethrin; only 11% and 33% mortalities for 5× and 10× permethrin diagnostic concentrations, and <2% mortalities for deltamethrin and DDT, respectively. Moderate bendiocarb resistance (88% mortality) and full susceptibility to malathion were observed. Synergist bioassays with piperonyl butoxide recovered permethrin susceptibility, with mortalities increasing to 53.39%, and 87.30% for 5× and 10× permethrin, respectively, implicating P450 monooxygenases. Synergist bioassays with diethyl maleate (DEM) recovered permethrin and DDT susceptibilities (mortalities increasing to 34.75% and 14.88%, respectively), implicating glutathione S-transferases. RNA-seq-based genome-wide transcriptional analyses supported by quantitative PCR identified glutathione S-transferase, GSTe2 (RNA-seqFC = 2.93 and qRT-PCRFC = 8.4, p < 0.0043) and CYP450, CYP6Z2 (RNA-seqFC = 2.39 and qRT-PCRFC = 11.7, p < 0.0177) as the most overexpressed detoxification genes in the pyrethroid-resistant mosquitoes, compared to mosquitoes of the susceptible Ngousso colony. Other overexpressed genes include P450s, CYP6M2 (FC = 1.68, p < 0.0114), CYP4G16 (FC = 2.02, p < 0.0005), and CYP4G17 (FC = 1.86, p < 0.0276). While high frequency of the 1014F kdr mutation (50%) and low frequencies of 1014S (6.61%) and 1575Y (10.29%) were observed, no ace-1 mutation was detected in bendiocarb-resistant populations, suggesting the preeminent role of metabolic mechanism. Overexpression of metabolic resistance genes (including GSTe2 and CYP6Z2 known to confer resistance to multiple insecticides) in An. coluzzii from the Sudan Savannah of Cameroon highlights the need for alternative management strategies to reduce malaria burden in northern Cameroon.

Parallel evolution in mosquito vectors - a duplicated esterase locus is associated with resistance to pirimiphos-methyl in An. gambiae.

The primary control methods for the African malaria mosquito, Anopheles gambiae, are based on insecticidal interventions. Emerging resistance to these compounds is therefore of major concern to malaria control programmes. The organophosphate, pirimiphos-methyl, is a relatively new chemical in the vector control armoury but is now widely used in indoor residual spray campaigns. Whilst generally effective, phenotypic resistance has developed in some areas in malaria vectors. Here, we used a population genomic approach to identify novel mechanisms of resistance to pirimiphos-methyl in Anopheles gambiae s.l mosquitoes. In multiple populations, we found large and repeated signals of selection at a locus containing a cluster of detoxification enzymes, some of whose orthologs are known to confer resistance to organophosphates in Culex pipiens. Close examination revealed a pair of alpha-esterases, Coeae1f and Coeae2f, and a complex and diverse pattern of haplotypes under selection in An. gambiae, An. coluzzii and An. arabiensis. As in Cx. pipiens, copy number variation seems to play a role in the evolution of insecticide resistance at this locus. We used diplotype clustering to examine whether these signals arise from parallel evolution or adaptive introgression. Using whole-genome sequenced phenotyped samples, we found that in West Africa, a copy number variant in Anopheles gambiae is associated with resistance to pirimiphos-methyl. Overall, we demonstrate a striking example of contemporary parallel evolution which has important implications for malaria control programmes.

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.

Discovery of knock-down resistance in the major African malaria vector Anopheles funestus.

A major mechanism of insecticide resistance in insect pests is knock-down resistance (kdr) caused by mutations in the voltage-gated sodium channel (Vgsc) gene. Despite being common in most malaria Anopheles vector species, kdr mutations have never been observed in Anopheles funestus, the principal malaria vector in Eastern and Southern Africa. While monitoring 10 populations of An. funestus in Tanzania, we unexpectedly found resistance to DDT, a banned insecticide, in one location. Through whole-genome sequencing of 333 An. funestus samples from these populations, we found 8 novel amino acid substitutions in the Vgsc gene, including the kdr variant, L976F (L1014F in An. gambiae), in tight linkage disequilibrium with another (P1842S). The mutants were found only at high frequency in one region, with a significant decline between 2017 and 2023. Notably, kdr L976F was strongly associated with survivorship to the exposure to DDT insecticide, while no clear association was noted with a pyrethroid insecticide (deltamethrin). Further study is necessary to identify the origin and spread of kdr in An. funestus, and the potential threat to current insecticide-based vector control in Africa.

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.

Mechanisms of transcriptional regulation in Anopheles gambiae revealed by allele specific expression.

Malaria control relies on insecticides targeting the mosquito vector, but this is increasingly compromised by insecticide resistance, which can be achieved by elevated expression of detoxifying enzymes that metabolize the insecticide. In diploid organisms, gene expression is regulated both in cis, by regulatory sequences on the same chromosome, and by trans acting factors, affecting both alleles equally. Differing levels of transcription can be caused by mutations in cis-regulatory modules (CRM), but few of these have been identified in mosquitoes. We crossed bendiocarb resistant and susceptible Anopheles gambiae strains to identify cis-regulated genes that might be responsible for the resistant phenotype using RNAseq, and cis-regulatory module sequences controlling gene expression in insecticide resistance relevant tissues were predicted using machine learning. We found 115 genes showing allele specific expression in hybrids of insecticide susceptible and resistant strains, suggesting cis regulation is an important mechanism of gene expression regulation in Anopheles gambiae. The genes showing allele specific expression included a higher proportion of Anopheles specific genes on average younger than genes those with balanced allelic expression.

Genome-wide association studies reveal novel loci associated with pyrethroid and organophosphate resistance in Anopheles gambiae s.l.

Resistance to insecticides in Anopheles mosquitoes threatens the effectiveness of the most widespread tools currently used to control malaria. The genetic underpinnings of resistance are still only partially understood, with much of the variance in resistance phenotype left unexplained. We performed a multi-country large scale genome-wide association study of resistance to two insecticides widely used in malaria control: deltamethrin and pirimiphos-methyl. Using a bioassay methodology designed to maximise the phenotypic difference between resistant and susceptible samples, we sequenced 969 phenotyped female An. gambiae and An. coluzzii from ten locations across four countries in West Africa (Benin, Côte d'Ivoire, Ghana and Togo), identifying single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) segregating in the populations. The patterns of resistance association were highly multiallelic and variable between populations, with different genomic regions contributing to resistance, as well as different mutations within a given region. While the strongest and most consistent association with deltamethrin resistance came from the region around Cyp6aa1 , this resistance was based on a combination of several independent CNVs in An. coluzzii , and on a non-CNV bearing haplotype in An. gambiae . Further signals involved a range of cytochrome P450, mitochondrial, and immunity genes. Similarly, for pirimiphos-methyl, while the strongest signal came from the region of Ace1 , more widespread signals included cytochrome P450s, glutathione S-transferases, and a subunit of the nAChR target site of neonicotinoid insecticides. The regions around Cyp9k1 and the Tep family of immune genes were associated with resistance to both insecticide classes, suggesting possible cross-resistance mechanisms. These locally-varying, multigenic and multiallelic patterns highlight the challenges involved in genomic monitoring and surveillance of resistance, and form the basis for improvement of methods used to detect and predict resistance. Based on simulations of resistance variants, we recommend that yet larger scale studies, exceeding 500 phenotyped samples per population, are required to better identify associated genomic regions.

Genome-wide association studies reveal novel loci associated with pyrethroid and organophosphate resistance in Anopheles gambiae and Anopheles coluzzii.

Resistance to insecticides in Anopheles mosquitoes threatens the effectiveness of malaria control, but the genetics of resistance are only partially understood. We performed a large scale multi-country genome-wide association study of resistance to two widely used insecticides: deltamethrin and pirimiphos-methyl, using sequencing data from An. gambiae and An. coluzzii from ten locations in West Africa. Resistance was highly multi-genic, multi-allelic and variable between populations. While the strongest and most consistent association with deltamethrin resistance came from Cyp6aa1, this was based on several independent copy number variants (CNVs) in An. coluzzii, and on a non-CNV haplotype in An. gambiae. For pirimiphos-methyl, signals included Ace1, cytochrome P450s, glutathione S-transferases and the nAChR target site of neonicotinoid insecticides. The regions around Cyp9k1 and the Tep family of immune genes showed evidence of cross-resistance to both insecticides. These locally-varying, multi-allelic patterns highlight the challenges involved in genomic monitoring of resistance, and may form the basis for improved surveillance methods.

Capturing the transcription factor interactome in response to sub-lethal insecticide exposure.

The increasing levels of pesticide resistance in agricultural pests and disease vectors represents a threat to both food security and global health. As insecticide resistance intensity strengthens and spreads, the likelihood of a pest encountering a sub-lethal dose of pesticide dramatically increases. Here, we apply dynamic Bayesian networks to a transcriptome time-course generated using sub-lethal pyrethroid exposure on a highly resistant Anopheles coluzzii population. The model accounts for circadian rhythm and ageing effects allowing high confidence identification of transcription factors with key roles in pesticide response. The associations generated by this model show high concordance with lab-based validation and identifies 44 transcription factors putatively regulating insecticide-responsive transcripts. We identify six key regulators, with each displaying differing enrichment terms, demonstrating the complexity of pesticide response. The considerable overlap of resistance mechanisms in agricultural pests and disease vectors strongly suggests that these findings are relevant in a wide variety of pest species.

High concentrations of membrane-fed ivermectin are required for substantial lethal and sublethal impacts on Aedes aegypti.

BACKGROUND: With widespread insecticide resistance in mosquito vectors, there is a pressing need to evaluate alternatives with different modes of action. Blood containing the antihelminthic drug ivermectin has been shown to have lethal and sub-lethal effects on mosquitoes. Almost all work to date has been on Anopheles spp., but impacts on other anthropophagic vectors could provide new options for their control, or additional value to anti-malarial ivermectin programmes. METHODS: Using dose-response assays, we evaluated the effects of ivermectin delivered by membrane feeding on daily mortality (up to 14 days post-blood feed) and fecundity of an Indian strain of Aedes aegypti. RESULTS: The 7-day lethal concentration of ivermectin required to kill 50% of adult mosquitoes was calculated to be 178.6 ng/ml (95% confidence intervals 142.3-218.4) for Ae. aegypti, which is much higher than that recorded for Anopheles spp. in any previous study. In addition, significant effects on fecundity and egg hatch rates were only recorded at high ivermectin concentrations (≥ 250 ng/ul). CONCLUSION: Our results suggest that levels of ivermectin present in human blood at current dosing regimes in mass drug administration campaigns, or even those in a recent higher-dose anti-malaria trial, are unlikely to have a substantial impact on Ae. aegypti. Moreover, owing to the strong anthropophagy of Ae. aegypti, delivery of higher levels of ivermectin in livestock blood is also unlikely to be an effective option for its control. However, other potential toxic impacts of ivermectin metabolites, accumulation in tissues, sublethal effects on behaviour, or antiviral action might increase the efficacy of ivermectin against Ae. aegypti and the arboviral diseases it transmits, and require further investigation.