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

Genome assemblies and other genomic tools for understanding insect adaptation

Insects, the most diverse group of animals, exhibit remarkable adaptability, playing both crucial and problematic roles in ecosystems. Recent advancements in genomic technologies, such as high-throughput sequencing, have provided unprecedented insights into the genetic foundations of insect adaptation. This review explores key methodologies, including de novo and reference-guided genome assemblies, and highlights cutting-edge technologies like second- and third-generation sequencing and hybrid techniques. The article delves into the genetic mechanisms underlying insect adaptations, focusing on structural variants. Case studies, such as the Anopheles gambiae genome assembly and the genomic research on Drosophila melanogaster, demonstrate the practical applications of these technologies in understanding pesticide resistance, climate adaptation, and other evolutionary traits. This review underscores the transformative role of genomic tools in insect research, with significant implications for pest management, agriculture, and biodiversity conservation.

Evolutionary consequences of long-distance dispersal in mosquitoes

Long-distance dispersal (LDD) provides a means for mosquitoes to invade new regions and spread adaptive alleles, including those conferring insecticide resistance. Most LDD takes place on human transport vessels and will typically be rarer and more directionally constrained than active flight but can connect populations and regions that are otherwise mutually inaccessible. These features make LDD worthy of specific consideration in mosquito research. This paper reviews recent evolutionary research on LDD and its consequences for mosquito populations and mosquito control. LDD is the main source of mosquito range expansions, and genomic methods can now trace the origins of new invasions to specific towns or cities. Genomic methods can also give a rough indication of the number of invaders, which if very small may lead to the stochastic loss of advantageous alleles during invasion bottlenecks. Once invasions are established, LDD spreads adaptive alleles between populations. Emerging insights into insecticide resistance evolution indicate that LDD has repeatedly spread resistance mutations across global species ranges, but these broad patterns are convoluted by two other evolutionary processes: parallel adaptation at the same gene or gene cluster and polygenic adaptation at different genes in different populations. Together, these processes have produced patterns of similarity and dissimilarity at resistance genes that are decoupled from geographical distance. LDD within cities is less well studied but is important for planning and evaluating local control efforts. Urban investigations of LDD may help identify areas experiencing weaker selection pressures from insecticides and isolated areas to target for control.

Targeted genomic surveillance of insecticide resistance in African malaria vectors

The emergence of insecticide resistance is threatening the efforts of malaria control programmes, which rely heavily on a limited arsenal of insecticidal tools, such as insecticide-treated bed nets. Importantly, genomic surveillance of malaria vectors can provide critical, policy-relevant insights into the presence and evolution of insecticide resistance, allowing us to maintain and extend the shelf life of these interventions. Yet the complex genetic architecture of resistance, combined with resource constraints in malaria-endemic settings, have thus far precluded the widespread use of genomics in routine surveillance. Meanwhile, stakeholders in sub-Saharan Africa are moving towards locally driven, decentralised generation of genomic data, underscoring the need for standardised and robust genomics workflows. To address this need, we demonstrate an approach to targeted genomic surveillance in Anopheles gambiae s.l with Illumina sequencing. We target 90 genomic loci in the Anopheles gambiae s.l genome, including 55 resistance-associated mutations and 35 ancestry informative markers. This protocol is coupled with advanced, automated software for accurate and reproducible variant analysis. We are able to elucidate population structure and ancestry in our cohorts and accurately identify most species in the An. gambiae species complex. We report frequencies of variants at insecticide-resistance loci and explore the continued evolution of the pyrethroid target site, the Voltage-gated sodium channel. Applying the platform to a recently established colony of field-caught resistant mosquitoes (Siaya, Kenya), we identified seven independent resistance-associated variants contributing to reduced efficacy of insecticide-treated nets in East Africa. Additionally, we leverage a machine learning algorithm (XGBoost) to demonstrate the possibility of predicting bioassay mortality using genotypes alone. This achieved very high accuracy (73%), demonstrating the potential of targeted genomics to predictively monitor insecticide resistance. Together these tools provide a practical, scalable solution for resistance monitoring while advancing the goal of building local genomic surveillance capacity in sub-Saharan Africa.

Ecological risks caused by neonicotinoid pesticides in sediments: A case study of freshwater basins in China

Neonicotinoid insecticides (NNIs) are extensively used in agricultural production in China due to their selective neurotoxicity towards target insects. In recent years, the rapid development of agriculture has increased the use and residue of NNIs. Consequently, the sediment environment, serving as the ultimate sink, is significantly impacted by NNIs. Upon release into the environment, NNIs can enter the human body through the food chain, posing potential ecological and human health risks. This study analyzed 79 sediment samples from two major river basins in North and South China, the Liaohe River basin in Liaoning Province and the Jianjiang River basin in Guangdong Province. The content, composition, distribution, and source of eight NNIs were analyzed, and assess the ecological and human health risks of the target compounds in these regions. The results indicated that the average concentration of NNIs in the sediments of the Jianjiang River basin (2.34 μg/kg) is slightly higher than that of the Liaohe River basin (2.32 μg/kg), and the sources of NNIs in the two areas were different, with differences in the sources of NNIs likely attributable to varying types of agricultural production. The risk assessment revealed that the ecotoxicological and public health risks were more pronounced in the Jianjiang River basin compared to the Liaohe River basin, underscoring the critical need for surveillance and management of hazardous substances like NNIs. The insights findings from this study can provide scientific guidance for the risk evaluation and environmental management of NNIs.

Genomic diversity of the African malaria vector Anopheles funestus

Anopheles funestus s.s. is a formidable human malaria vector across sub-Saharan Africa. To understand how the species is evolving, especially in response to malaria vector control, we sequenced 656 modern specimens (collected 2014-2018) and 45 historic specimens (collected 1927-1967) from 16 African countries. We find high levels of genetic variation with clear and stable continental patterns. Six segregating inversions might be involved in adaptation of local ecotypes. Strong recent signals of selection centred on canonical insecticide resistance genes are shared by multiple populations. A promising gene drive target in An. gambiae is highly conserved in An. funestus. This work represents a significant advance in our understanding of the genetic diversity and population structure of An. funestus and will enable smarter targeted malaria control.

Field-evolved resistance to neonicotinoids in the mosquito, Anopheles gambiae, is associated with mutations of nicotinic acetylcholine receptor subunits combined with cytochrome P450-mediated detoxification

New insecticides prequalified for malaria control interventions include modulators of nicotinic acetylcholine receptors that act selectively on different subunits leading to variable sensitivity among arthropods. This study aimed to investigate the molecular mechanisms underlying contrasting susceptibility to neonicotinoids observed in wild populations of two mosquito sibling species. Bioassays and a synergist test with piperonyl butoxide revealed that the sister taxa, Anopheles gambiae and An. coluzzii, from Yaounde, Cameroon, both have the potential to develop resistance to acetamiprid through cytochrome P450-mediated detoxification. However, contrary to An. coluzzii, An. gambiae populations are evolving cross-resistance to several active ingredients facilitated by mutations of nicotinic acetylcholine receptors (nAChRs). We sequenced coding regions on the β1 and α6 nAChR subunits where variants associated with resistance to neonicotinoids or to spinosyns have been found in agricultural pests and detected no mutation in An. coluzzii. By contrast, six nucleotide substitutions including an amino acid change in one of the loops that modulate ligand binding and affect sensitivity were present in the resistant species, An. gambiae. Allele frequency distributions were consistent with the spread of beneficial mutations that likely reduce the affinity of An. gambiae nAChRs for synthetic modulators. Our findings provide critical information for the application and resistance management of nAChR modulators in malaria prevention.

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

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.

Discovery of Knock-Down Resistance in the Major African Malaria Vector Anopheles funestus

A major insecticide resistance mechanism 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, with resistance mainly being conferred by detoxification enzymes. In a parallel study, we monitored 10 populations of An. funestus in Tanzania for insecticide resistance unexpectedly identified resistance to a banned insecticide, DDT, in the Morogoro region. Through whole-genome sequencing of 333 An. funestus samples from these populations, we found eight novel amino acid substitutions in the Vgsc gene, including the kdr variant, L976F (equivalent to L995F in An. gambiae), in tight linkage disequilibrium with another (P1842S). The mutants were found only at high frequency in one region and were accompanied by weak signatures of a selective sweep, with a significant decline between 2017 and 2023. Notably, kdr L976F was strongly associated with survivorship to exposure to DDT insecticide, while no clear association was noted with a pyrethroid insecticide (deltamethrin). The WHO prequalifies no DDT products for vector control, and the chemical is banned in Tanzania. Widespread DDT contamination and a legacy of extensive countrywide stockpiles may have selected for this mutation. Continued monitoring is necessary to understand the origin of kdr in An. funestus, and the threat posed to insecticide-based vector control in Africa.

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.

Mitochondrial Variation in Anopheles gambiae and Anopheles coluzzii: Phylogeographic Legacy and Mitonuclear Associations With Metabolic Resistance to Pathogens and Insecticides

Mitochondrial DNA has been a popular marker in phylogeography, phylogeny, and molecular ecology, but its complex evolution is increasingly recognized. Here, we investigated mitochondrial DNA variation in Anopheles gambiae and Anopheles coluzzii, in relation to other species in the Anopheles gambiae complex, by assembling the mitogenomes of 1,219 mosquitoes across Africa. The mitochondrial DNA phylogeny of the Anopheles gambiae complex was consistent with previously reported highly reticulated evolutionary history, revealing important discordances with the species tree. The three most widespread species (An. gambiae, An. coluzzii, and Anopheles arabiensis), known for extensive historical introgression, could not be discriminated based on mitogenomes. Furthermore, a monophyletic clustering of the three saltwater-tolerant species (Anopheles merus, Anopheles melas, and Anopheles bwambae) in the Anopheles gambiae complex also suggested that introgression and possibly selection shaped mitochondrial DNA evolution. Mitochondrial DNA variation in An. gambiae and An. coluzzii across Africa revealed significant partitioning among populations and species. A peculiar mitochondrial DNA lineage found predominantly in An. coluzzii and in the hybrid taxon of the African "far-west" exhibited divergence comparable to the interspecies divergence in the Anopheles gambiae complex, with a geographic distribution matching closely An. coluzzii's geographic range. This phylogeographic relict of the An. coluzzii and An. gambiae split was associated with population and species structure, but not with the rare Wolbachia occurrence. The lineage was significantly associated with single nucleotide polymorphisms in the nuclear genome, particularly in genes associated with pathogen and insecticide resistance. These findings underline potential mitonuclear coevolution history and the role played by mitochondria in shaping metabolic responses to pathogens and insecticides in Anopheles.

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 CRM sequences controlling gene expression in insecticide resistance relevant tissues were predicted using machine learning. We found 115 genes showing allele-specific expression (ASE) in hybrids of insecticide susceptible and resistant strains, suggesting cis-regulation is an important mechanism of gene expression regulation in A. gambiae. The genes showing ASE included a higher proportion of Anopheles-specific genes on average younger than genes with balanced allelic expression.

Field-evolved resistance to neonicotinoids in the mosquito, Anopheles gambiae, is associated with downregulation and mutations of nicotinic acetylcholine receptor subunits combined with cytochrome P450-mediated detoxification

Neonicotinoid insecticides act selectively on their nicotinic receptor targets leading to variable sensitivity among arthropods. This study aimed to investigate the molecular mechanisms underlying contrasting susceptibility to neonicotinoids observed in wild populations of two mosquito sibling species. Bioassays and a synergism test revealed that the sister taxa, Anopheles gambiae and An. coluzzii, from Yaounde, Cameroon, rely on cytochrome P450s to detoxify neonicotinoids and develop resistance. However, contrary to An. coluzzii, An. gambiae populations are evolving stronger resistance to several active ingredients facilitated by mutations and reduced expression of nicotinic acetylcholine receptors. Six mutations were detected in coding sequences of the β1 and α6 subunits, including two substitutions in one of the loops that modulate ligand binding and sensitivity. Allele frequencies were strongly correlated with a susceptibility gradient between An. coluzzii and An. gambiae suggesting that the mutations may play a key role in sensitivity. Messenger RNA expression levels of the β1, α3, and α7 subunits decreased dramatically, on average by 23.27, 17.50, 15.80-fold, respectively, in wild An. gambiae populations compared to a susceptible insectary colony. By contrast, only the β2 and α9-1 subunits were moderately downregulated (5.28 and 2.67-fold change, respectively) in field-collected An. coluzzii adults relative to susceptible colonized mosquitoes. Our findings provide critical information for the application and resistance management of neonicotinoids in malaria prevention.

Key gene modules and hub genes associated with pyrethroid and organophosphate resistance in Anopheles mosquitoes: a systems biology approach

Indoor residual spraying (IRS) and insecticide-treated nets (ITNs) are the main methods used to control mosquito populations for malaria prevention. The efficacy of these strategies is threatened by the spread of insecticide resistance (IR), limiting the success of malaria control. Studies of the genetic evolution leading to insecticide resistance could enable the identification of molecular markers that can be used for IR surveillance and an improved understanding of the molecular mechanisms associated with IR. This study used a weighted gene co-expression network analysis (WGCNA) algorithm, a systems biology approach, to identify genes with similar co-expression patterns (modules) and hub genes that are potential molecular markers for insecticide resistance surveillance in Kenya and Benin. A total of 20 and 26 gene co-expression modules were identified via average linkage hierarchical clustering from Anopheles arabiensis and An. gambiae, respectively, and hub genes (highly connected genes) were identified within each module. Three specific genes stood out: serine protease, E3 ubiquitin-protein ligase, and cuticular proteins, which were top hub genes in both species and could serve as potential markers and targets for monitoring IR in these malaria vectors. In addition to the identified markers, we explored molecular mechanisms using enrichment maps that revealed a complex process involving multiple steps, from odorant binding and neuronal signaling to cellular responses, immune modulation, cellular metabolism, and gene regulation. Incorporation of these dynamics into the development of new insecticides and the tracking of insecticide resistance could improve the sustainable and cost-effective deployment of interventions.

Parallel Evolution in Mosquito Vectors-A Duplicated Esterase Locus is Associated With Resistance to Pirimiphos-methyl in Anopheles 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 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.

Genetic surveillance of insecticide resistance in African Anopheles populations to inform malaria vector control

Insecticide resistance in malaria vector populations poses a major threat to malaria control, which relies largely on insecticidal interventions. Contemporary vector-control strategies focus on combatting resistance using multiple insecticides with differing modes of action within the mosquito. However, diverse genetic resistance mechanisms are present in vector populations, and continue to evolve. Knowledge of the spatial distribution of these genetic mechanisms, and how they impact the efficacy of different insecticidal products, is critical to inform intervention deployment decisions. We developed a catalogue of genetic-resistance mechanisms in African malaria vectors that could guide molecular surveillance. We highlight situations where intervention deployment has led to resistance evolution and spread, and identify challenges in understanding and mitigating the epidemiological impacts of resistance.

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.

Genetic markers associated with the widespread insecticide resistance in malaria vector Anopheles funestus populations across Tanzania

BACKGROUND

Anopheles funestus is a leading vector of malaria in most parts of East and Southern Africa, yet its ecology and responses to vector control remain poorly understood compared with other vectors such as Anopheles gambiae and Anopheles arabiensis. This study presents the first large-scale survey of the genetic and phenotypic expression of insecticide resistance in An. funestus populations in Tanzania.

METHODS

We performed insecticide susceptibility bioassays on An. funestus mosquitoes in nine regions with moderate-to-high malaria prevalence in Tanzania, followed by genotyping for resistance-associated mutations (CYP6P9a, CYP6P9b, L119F-GSTe2) and structural variants (SV4.3 kb, SV6.5 kb). Generalized linear models were used to assess relationships between genetic markers and phenotypic resistance. An interactive R Shiny tool was created to visualize the data and support evidence-based interventions.

RESULTS

Pyrethroid resistance was universal but reversible by piperonyl-butoxide (PBO). However, carbamate resistance was observed in only five of the nine districts, and dichloro-diphenyl-trichloroethane (DDT) resistance was found only in the Kilombero valley, south-eastern Tanzania. Conversely, there was universal susceptibility to the organophosphate pirimiphos-methyl in all sites. Genetic markers of resistance had distinct geographical patterns, with CYP6P9a-R and CYP6P9b-R alleles, and the SV6.5 kb structural variant absent or undetectable in the north-west but prevalent in all other sites, while SV4.3 kb was prevalent in the north-western and western regions but absent elsewhere. Emergent L119F-GSTe2, associated with deltamethrin resistance, was detected in heterozygous form in districts bordering Mozambique, Malawi and the Democratic Republic of Congo. The resistance landscape was most complex in western Tanzania, in Tanganyika district, where all five genetic markers were detected. There was a notable south-to-north spread of resistance genes, especially CYP6P9a-R, though this appears to be interrupted, possibly by the Rift Valley.

CONCLUSIONS

This study underscores the need to expand resistance monitoring to include An. funestus alongside other vector species, and to screen for both the genetic and phenotypic signatures of resistance. The findings can be visualized online via an interactive user interface and could inform data-driven decision-making for resistance management and vector control. Since this was the first large-scale survey of resistance in Tanzania's An. funestus, we recommend regular updates with greater geographical and temporal coverage.

Signatures of adaptation at key insecticide resistance loci in Anopheles gambiae in Southern Ghana revealed by reduced-coverage WGS

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.

Transcriptomic analysis of Anopheles gambiae from Benin reveals overexpression of salivary and cuticular proteins associated with cross-resistance to pyrethroids and organophosphates

BACKGROUND

Insecticide resistance (IR) is one of the major threats to malaria vector control programs in endemic countries. However, the mechanisms underlying IR are poorly understood. Thus, investigating gene expression patterns related to IR can offer important insights into the molecular basis of IR in mosquitoes. In this study, RNA-Seq was used to characterize gene expression in Anopheles gambiae surviving exposure to pyrethroids (deltamethrin, alphacypermethrin) and an organophosphate (pirimiphos-methyl).

RESULTS

Larvae of An. gambiae s.s. collected from Bassila and Djougou in Benin were reared to adulthood and phenotyped for IR using a modified CDC intensity bottle bioassay. The results showed that mosquitoes from Djougou were more resistant to pyrethroids (5X deltamethrin: 51.7% mortality; 2X alphacypermethrin: 47.4%) than Bassila (1X deltamethrin: 70.7%; 1X alphacypermethrin: 77.7%), while the latter were more resistant to pirimiphos-methyl (1.5X: 48.3% in Bassila and 1X: 21.5% in Djougou). RNA-seq was then conducted on resistant mosquitoes, non-exposed mosquitoes from the same locations and the laboratory-susceptible An. gambiae s.s. Kisumu strain. The results showed overexpression of detoxification genes, including cytochrome P450s (CYP12F2, CYP12F3, CYP4H15, CYP4H17, CYP6Z3, CYP9K1, CYP4G16, and CYP4D17), carboxylesterase genes (COEJHE5E, COE22933) and glutathione S-transferases (GSTE2 and GSTMS3) in all three resistant mosquito groups analyzed. Genes encoding cuticular proteins (CPR130, CPR10, CPR15, CPR16, CPR127, CPAP3-C, CPAP3-B, and CPR76) were also overexpressed in all the resistant groups, indicating their potential role in cross resistance in An. gambiae. Salivary gland protein genes related to 'salivary cysteine-rich peptide' and 'salivary secreted mucin 3' were also over-expressed and shared across all resistant groups.

CONCLUSION

Our results suggest that in addition to metabolic enzymes, cuticular and salivary gland proteins could play an important role in cross-resistance to multiple classes of insecticides in Benin. These genes warrant further investigation to validate their functional role in An. gambiae resistance to insecticides.

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.

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.