Search

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. The outdoor-biting malaria vector Anopheles arabiensis is of increasing concern for malaria transmission because it is apparently less susceptible to many indoor control interventions, yet knowledge of its mechanisms of resistance remains limited. 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 globally high resistance to pyrethroids. We performed a genome-wide association study to determine the molecular mechanisms of resistance to the pyrethroid deltamethrin (commonly used in bednets) and PM (widespread use for indoor spraying), in An. arabiensis from 2 regions in Tanzania. Genomic regions of positive selection in these populations were largely driven by copy number variants (CNVs) in gene families involved in metabolic resistance. We found evidence of a new gene cluster involved in resistance to PM, identifying a strong selective sweep tied to a CNV in the carboxylesterase genes Coeae2g - Coeae6g. Using complementary data from another malaria vector, 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 (Cyp6aap_Dup33). Against this background of metabolic resistance, resistance caused by mutations in the insecticide target sites was very rare or absent. Mutations in the pyrethroid target site Vgsc were at very low frequency in Tanzania, yet combining these samples with 3 An. arabiensis individuals from West Africa revealed a startling evolutionary diversity, 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., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Lucas et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

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 programs. The organophosphate (OP), pirimiphos-methyl, is a relatively new chemical in the vector control armory but is now widely used in indoor-residual spray campaigns. While 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 A. 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 OPs 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 A. gambiae, A. coluzzii and A. arabiensis. As in C. pipiens, copy number variants have arisen 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 A. 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 programs., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Resistance to insecticides and adaptation to a diverse range of environments present challenges to Anopheles gambiae s.l. mosquito control efforts in sub-Saharan Africa. Whole-genome-sequencing is often employed for identifying the genomic basis underlying adaptation in Anopheles, but remains expensive for large-scale surveys. Reduced coverage whole-genome-sequencing can identify regions of the genome involved in adaptation at a lower cost, but is currently untested in Anopheles mosquitoes. Here, we use reduced coverage WGS to investigate population genetic structure and identify signatures of local adaptation in Anopheles mosquitoes across southern Ghana. In contrast to previous analyses, we find no structuring by ecoregion, with Anopheles coluzzii and Anopheles gambiae populations largely displaying the hallmarks of large, unstructured populations. However, we find signatures of selection at insecticide resistance loci that appear ubiquitous across ecoregions in An. coluzzii, and strongest in forest ecoregions in An. gambiae. Our study highlights resistance candidate genes in this region, and validates reduced coverage WGS, potentially to very low coverage levels, for population genomics and exploratory surveys for adaptation in Anopheles taxa., (© 2024. The Author(s).)

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., (© 2023. Springer Nature Limited.)

Background: Resistance in malaria vectors to pyrethroids, the most widely used class of insecticides for malaria vector control, threatens the continued efficacy of vector control tools. Target-site resistance is an important genetic resistance mechanism caused by mutations in the voltage-gated sodium channel (Vgsc) gene that encodes the pyrethroid target-site. Understanding the geographic distribution of target-site resistance, and temporal trends across different vector species, can inform strategic deployment of vector control tools., Results: We develop a Bayesian statistical spatiotemporal model to interpret species-specific trends in the frequency of the most common resistance mutations, Vgsc-995S and Vgsc-995F, in three major malaria vector species Anopheles gambiae, An. coluzzii, and An. arabiensis over the period 2005-2017. The models are informed by 2418 observations of the frequency of each mutation in field sampled mosquitoes collected from 27 countries spanning western and eastern regions of Africa. For nine selected countries, we develop annual predictive maps which reveal geographically structured patterns of spread of each mutation at regional and continental scales. The results show associations, as well as stark differences, in spread dynamics of the two mutations across the three vector species. The coverage of ITNs was an influential predictor of Vgsc allele frequencies, with modelled relationships between ITN coverage and allele frequencies varying across species and geographic regions. We found that our mapped Vgsc allele frequencies are a significant partial predictor of phenotypic resistance to the pyrethroid deltamethrin in An. gambiae complex populations., Conclusions: Our predictive maps show how spatiotemporal trends in insecticide target-site resistance mechanisms in African An. gambiae vary across individual vector species and geographic regions. Molecular surveillance of resistance mechanisms will help to predict resistance phenotypes and track their spread., (© 2022. The Author(s).)

The spread of resistance to insecticides in disease-carrying mosquitoes poses a threat to the effectiveness of control programmes, which rely largely on insecticide-based interventions. Monitoring mosquito populations is essential, but obtaining phenotypic measurements of resistance is laborious and error-prone. High-throughput genotyping offers the prospect of quick and repeatable estimates of resistance, while also allowing resistance markers to be tracked and studied. To demonstrate the potential of highly-mulitplexed genotypic screening for measuring resistance-association of mutations and tracking their spread, we developed a panel of 28 known or putative resistance markers in the major malaria vector Anopheles gambiae, which we used to screen mosquitoes from a wide swathe of Sub-Saharan Africa (Burkina Faso, Ghana, Democratic Republic of Congo (DRC) and Kenya). We found resistance association in four markers, including a novel mutation in the detoxification gene Gste2 (Gste2-119V). We also identified a duplication in Gste2 combining a resistance-associated mutation with its wild-type counterpart, potentially alleviating the costs of resistance. Finally, we describe the distribution of the multiple origins of kdr resistance, finding unprecedented diversity in the DRC. This panel represents the first step towards a quantitative genotypic model of insecticide resistance that can be used to predict resistance status in An. gambiae.

Background: Identification of variation in Ace-1 copy number and G119S mutation genotype from samples of Anopheles gambiae and Anopheles coluzzii across West Africa are important diagnostics of carbamate and organophosphate resistance at population and individual levels. The most widespread and economical method, PCR-RFLP, suffers from an inability to discriminate true heterozygotes from heterozygotes with duplication., Methods: In addition to PCR-RFLP, in this study three different molecular techniques were applied on the same mosquito specimens: TaqMan qPCR, qRTPCR and ddPCR. To group heterozygous individuals recorded from the PCR-RFLP analysis into different assumptive genotypes K-means clustering was applied on the Z-scores of data obtained from both the TaqMan and ddPCR methods. The qRTPCR analysis was used for absolute quantification of copy number variation., Results: The results indicate that most heterozygotes are duplicated and that G119S mutation must now be regarded as a complex genotype ranging from primarily single-copy susceptible Glycine homozygotes to balanced and imbalanced heterozygotes, and multiply-amplified resistant Serine allele homozygotes. Whilst qRTPCR-based gene copy analysis suffers from some imprecision, it clearly illustrates differences in copy number among genotype groups identified by TaqMan or ddPCR. Based on TaqMan method properties, and by coupling TaqMan and ddPCR methods simultaneously on the same type of mosquito specimens, it demonstrated that the TaqMan genotype assays associated with the K-means clustering algorithm could provide a useful semi-quantitative estimate method to investigate the level of allele-specific duplication in mosquito populations., Conclusions: Ace-1 gene duplication is evidently far more complex in An. gambiae and An. coluzzii than the better-studied mosquito Culex quinquefasciatus, which consequently can no longer be considered an appropriate model for prediction of phenotypic consequences. These require urgent further evaluation in Anopheles. To maintain the sustained effectiveness carbamates and organophosphates as alternative products to pyrethroids for malaria vector control, monitoring of duplicated resistant alleles in natural populations is essential to guide the rational use of these insecticides.

Functionally constrained genes are ideal insecticide targets because disruption is often fatal, and resistance mutations are typically costly. Synaptic acetylcholinesterase (AChE) is an essential neurotransmission enzyme targeted by insecticides used increasingly in malaria control. In Anopheles and Culex mosquitoes, a glycine-serine substitution at codon 119 of the Ace-1 gene confers both resistance and fitness costs, especially for 119S/S homozygotes. G119S in Anopheles gambiae from Accra (Ghana) is strongly associated with resistance, and, despite expectations of cost, resistant 119S alleles are increasing significantly in frequency. Sequencing of Accra females detected only a single Ace-1 119S haplotype, whereas 119G diversity was high overall but very low at non-synonymous sites, evidence of strong purifying selection driven by functional constraint. Flanking microsatellites showed reduced diversity, elevated linkage disequilibrium and high differentiation of 119S, relative to 119G homozygotes across up to two megabases of the genome. Yet these signals of selection were inconsistent and sometimes weak tens of kilobases from Ace-1. This unexpected finding is attributable to apparently ubiquitous amplification of 119S alleles as part of a large copy number variant (CNV) far exceeding the size of the Ace-1 gene, whereas 119G alleles were unduplicated. Ace-1 CNV was detectable in archived samples collected when the 119S allele was rare in Ghana. Multicopy amplification of resistant alleles has not been observed previously and is likely to underpin the recent increase in 119S frequency. The large CNV compromised localization of the strong selective sweep around Ace-1, emphasizing the need to integrate CNV analysis into genome scans for selection., (© 2015 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Adaptive introgression can provide novel genetic variation to fuel rapid evolutionary responses, though it may be counterbalanced by potential for detrimental disruption of the recipient genomic background. We examine the extent and impact of recent introgression of a strongly selected insecticide-resistance mutation (Vgsc-1014F) located within one of two exceptionally large genomic islands of divergence separating the Anopheles gambiae species pair. Here we show that transfer of the Vgsc mutation results in homogenization of the entire genomic island region (~1.5% of the genome) between species. Despite this massive disruption, introgression is clearly adaptive with a dramatic rise in frequency of Vgsc-1014F and no discernable impact on subsequent reproductive isolation between species. Our results show (1) how resilience of genomes to massive introgression can permit rapid adaptive response to anthropogenic selection and (2) that even extreme prominence of genomic islands of divergence can be an unreliable indicator of importance in speciation.

Malaria control relies heavily on pyrethroid insecticides, to which susceptibility is declining in Anopheles mosquitoes. To combat pyrethroid resistance, application of alternative insecticides is advocated for indoor residual spraying (IRS), and carbamates are increasingly important. Emergence of a very strong carbamate resistance phenotype in Anopheles gambiae from Tiassalé, Côte d'Ivoire, West Africa, is therefore a potentially major operational challenge, particularly because these malaria vectors now exhibit resistance to multiple insecticide classes. We investigated the genetic basis of resistance to the most commonly-applied carbamate, bendiocarb, in An. gambiae from Tiassalé. Geographically-replicated whole genome microarray experiments identified elevated P450 enzyme expression as associated with bendiocarb resistance, most notably genes from the CYP6 subfamily. P450s were further implicated in resistance phenotypes by induction of significantly elevated mortality to bendiocarb by the synergist piperonyl butoxide (PBO), which also enhanced the action of pyrethroids and an organophosphate. CYP6P3 and especially CYP6M2 produced bendiocarb resistance via transgenic expression in Drosophila in addition to pyrethroid resistance for both genes, and DDT resistance for CYP6M2 expression. CYP6M2 can thus cause resistance to three distinct classes of insecticide although the biochemical mechanism for carbamates is unclear because, in contrast to CYP6P3, recombinant CYP6M2 did not metabolise bendiocarb in vitro. Strongly bendiocarb resistant mosquitoes also displayed elevated expression of the acetylcholinesterase ACE-1 gene, arising at least in part from gene duplication, which confers a survival advantage to carriers of additional copies of resistant ACE-1 G119S alleles. Our results are alarming for vector-based malaria control. Extreme carbamate resistance in Tiassalé An. gambiae results from coupling of over-expressed target site allelic variants with heightened CYP6 P450 expression, which also provides resistance across contrasting insecticides. Mosquito populations displaying such a diverse basis of extreme and cross-resistance are likely to be unresponsive to standard insecticide resistance management practices.

Background: With high DDT resistance present throughout much of West Africa, carbamates and organophosphates are increasingly important alternatives to pyrethroids for indoor residual spraying (IRS). Though less widespread, resistance to both of these alternative insecticide classes has also been documented within the Anopheles gambiae species pair (formerly the M and S molecular forms) in West Africa. To manage insecticide efficacy, it is important to predict how and where resistance is likely to occur and spread, which could be aided by using molecular diagnostics with high predictive value., Methods: Anopheles coluzzii and An. gambiae s.s. were collected from 18 sites throughout southern Ghana and bioassayed with bendiocarb, the most commonly applied carbamate, and an organophosphate, fenitrothion. The Ace-1 target site substitution G119S was genotyped by qPCR., Results: Fenitrothion induced higher mortality than bendiocarb, though phenotypes correlated strongly across populations. Ace-1 119S was found at much higher frequency in An. gambiae s.s than An. coluzzii, exceeding 90% in a population from Greater Accra, the highest frequency reported to date. Ace-1 G119S was very strongly associated with resistance to both insecticides, providing high predictive power for diagnosis, though with some evidence for a differential effect between molecular forms for bendiocarb. Sequencing of the gene revealed a lack of variation in resistant alleles precluding determination of origin, but Ace-1 copy number variation was detected for the first time in Ghana., Conclusions: The results validate G119S as a useful diagnostic of organophosphate and carbamate resistance within and among populations, whilst highlighting the potential for an aggregate nature of Ace-1 genotypes, which may comprise both single-copy and duplicated genes. Further work is now required to determine the distribution and resistance-association of Ace-1 duplication.