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

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.

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.

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.