The genetic architecture of target-site resistance to pyrethroid insecticides in the African malaria vectors Anopheles gambiae and Anopheles coluzzii

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.

Variation in spatial population structure in the Anopheles gambiae species complex

Anopheles gambiae , Anopheles coluzzii , and Anopheles arabiensis are three of the most widespread vectors of malaria parasites, with geographical ranges stretching across wide swaths of Africa. Understanding the population structure of these closely related species, including the extent to which populations are connected by gene flow, is essential for understanding how vector control implemented in one location might indirectly affect vector populations in other locations. Here, we assessed the population structure of each species based on a combined data set of publicly available and newly processed whole-genome sequences. The data set included single nucleotide polymorphisms from whole genomes of 2,410 individual mosquitoes sampled from 128 locations across 19 African countries. We found that A. gambiae sampled from several countries in West and Central Africa showed low genetic differentiation from each other according to principal components analysis (PCA) and ADMIXTURE modeling. Using Estimated Effective Migration Surfaces (EEMS), we showed that this low genetic differentiation indicates high effective migration rates for A. gambiae across this region. Outside of this region, we found eight groups of sampling locations from Central, East, and Southern Africa for which A. gambiae showed higher genetic differentiation, and lower effective migration rates, between each other and the West/Central Africa group. These results indicate that the barriers to and corridors for migration between populations of A. gambiae differ across the geographical range of this malaria vector species. Using the same methods, we found higher genetic differentiation and lower migration rates between populations of A. coluzzii in West and Central Africa than for A. gambiae in the same region. In contrast, we found lower genetic differentiation and higher migration rates between populations of A. arabiensis in Tanzania, compared to A. gambiae in the same region. These differences between A. gambiae , A. coluzzii , and A. arabiensis indicate that migration barriers and corridors may vary, even between very closely related species. Overall, our results demonstrate that migration rates vary both within and between species of Anopheles mosquitoes, presumably based on species-specific responses to the ecological or environmental conditions that may impede or facilitate migration, and the geographical patterns of these conditions across the landscape. Together with previous findings, this study provides robust evidence that migration rates between populations of malaria vectors depend on the ecological context, which should be considered when planning surveillance of vector populations, monitoring for insecticide resistance, and evaluating interventions.

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.

Whole-genome sequencing of major malaria vectors reveals the evolution of new insecticide resistance variants in a longitudinal study in Burkina Faso

BACKGROUND

Intensive deployment of insecticide based malaria vector control tools resulted in the rapid evolution of phenotypes resistant to these chemicals. Understanding this process at the genomic level is important for the deployment of successful vector control interventions. Therefore, longitudinal sampling followed by whole genome sequencing (WGS) is necessary to understand how these evolutionary processes evolve over time. This study investigated the change in genetic structure and the evolution of the insecticide resistance variants in natural populations of Anopheles gambiae over time and space from 2012 to 2017 in Burkina Faso.

METHODS

New genomic data have been generated from An. gambiae mosquitoes collected from three villages in the western part of Burkina Faso between 2012 and 2017. The samples were whole-genome sequenced and the data used in the An. gambiae 1000 genomes (Ag1000G) project as part of the Vector Observatory. Genomic data were analysed using the analysis pipeline previously designed by the Ag1000G project.

RESULTS

The results showed similar and consistent nucleotide diversity and negative Tajima's D between An. gambiae sensu stricto (s.s.) and Anopheles coluzzii. Principal component analysis (PCA) and the fixation index (FST) showed a clear genetic structure in the An. gambiae sensu lato (s.l.) species. Genome-wide FST and H12 scans identified genomic regions under divergent selection that may have implications in the adaptation to ecological changes. Novel voltage-gated sodium channel pyrethroid resistance target-site alleles (V402L, I1527T) were identified at increasing frequencies alongside the established alleles (Vgsc-L995F, Vgsc-L995S and N1570Y) within the An. gambiae s.l.

POPULATIONS

Organophosphate metabolic resistance markers were also identified, at increasing frequencies, within the An. gambiae s.s. populations from 2012 to 2017, including the SNP Ace1-G280S and its associated duplication. Variants simultaneously identified in the same vector populations raise concerns about the long-term efficacy of new generation bed nets and the recently organophosphate pirimiphos-methyl indoor residual spraying in Burkina Faso.

CONCLUSION

These findings highlighted the benefit of genomic surveillance of malaria vectors for the detection of new insecticide resistance variants, the monitoring of the existing resistance variants, and also to get insights into the evolutionary processes driving insecticide resistance.

Population genomic evidence of a putative 'far-west' African cryptic taxon in the Anopheles gambiae complex

The two main Afrotropical malaria vectors - Anopheles coluzzii and An. gambiae - are genetically distinct and reproductively isolated across West Africa. However, populations at the western extreme of their range are assigned as "intermediate" between the two species by whole genome sequence (WGS) data, and as hybrid forms by conventional molecular diagnostics. By exploiting WGS data from 1190 specimens collected across west Africa via the Anopheles gambiae 1000 Genomes network, we identified a putative taxon in the far-west (provisionally named Bissau molecular form), which did not arise by admixture but rather may have originated at the same time as the split between An. coluzzii and An. gambiae. Intriguingly, this taxon lacks insecticide resistance mechanisms commonly observed in the two main species. These findings lead to a change of perspective on malaria vector species in the far-west region with potential for epidemiological implications, and a new challenge for genetic-based mosquito control approaches.

Comparative Genomics Uncovers the Evolutionary Dynamics of Detoxification and Insecticide Target Genes Across 11 Phlebotomine Sand Flies

Sand flies infect more than 1 million people annually with Leishmania parasites and other bacterial and viral pathogens. Progress in understanding sand fly adaptations to xenobiotics has been hampered by the limited availability of genomic resources. To address this gap, we sequenced, assembled, and annotated the transcriptomes of 11 phlebotomine sand fly species. Subsequently, we leveraged these genomic resources to generate novel evolutionary insights pertaining to their adaptations to xenobiotics, including those contributing to insecticide resistance. Specifically, we annotated over 2,700 sand fly detoxification genes and conducted large-scale phylogenetic comparisons to uncover the evolutionary dynamics of the five major detoxification gene families: cytochrome P450s (CYPs), glutathione-S-transferases (GSTs), UDP-glycosyltransferases (UGTs), carboxyl/cholinesterases (CCEs), and ATP-binding cassette (ABC) transporters. Using this comparative approach, we show that sand flies have evolved diverse CYP and GST gene repertoires, with notable lineage-specific expansions in gene groups evolutionarily related to known xenobiotic metabolizers. Furthermore, we show that sand flies have conserved orthologs of (i) CYP4G genes involved in cuticular hydrocarbon biosynthesis, (ii) ABCB genes involved in xenobiotic toxicity, and (iii) two primary insecticide targets, acetylcholinesterase-1 (Ace1) and voltage gated sodium channel (VGSC). The biological insights and genomic resources produced in this study provide a foundation for generating and testing hypotheses regarding the molecular mechanisms underlying sand fly adaptations to xenobiotics.

First detection of kdr L1014F allele in Anopheles ziemanni and Anopheles pharoensis in Cameroon and distribution of the allele in members of the Anopheles gambiae complex

BACKGROUND

Knockdown resistance (kdr) is one of the primary resistance mechanisms present in anopheline species. Although this mutation is largely spread across the Anopheles gambiae s.l. members, its prevalence in other species is still not well documented.

METHODS

The present study investigated the distribution and allelic frequencies of kdr in An. gambiae s.l., An. pharoensis, and An. ziemanni samples collected in 2022 and 2023 in nine sites spread across five ecogeographical settings in Cameroon. Members of the An. gambiae complex were identified molecularly by polymerase chain reaction (PCR). kdr L1014F and L1014S alleles were screened by PCR and confirmed by sequencing.

RESULTS

An. gambiae (49.9%), An. coluzzii (36.5%), and An. arabiensis (13%) were identified, and the frequency of the kdr L1014F was high in both An. gambiae and An. coluzzii in all sites. The kdr L1014F allele was detected for the first time in 8 out of 14 An. ziemanni samples examined and in 5 out of 22 An. pharoensis samples examined. The kdr L1014S allele was scarce and found only in the heterozygote "RS" state in An. arabiensis and An. gambiae in Yangah and Santchou.

CONCLUSIONS

The present study sheds light on the rapid expansion of the kdr L1014F allele in malaria vectors in Cameroon and stresses the need for surveillance activities also targeting secondary malaria vectors to improve the control of malaria transmission.

Uridine diphosphate (UDP)-glycosyltransferases (UGTs) are associated with insecticide resistance in the major malaria vectors Anopheles gambiae s.l. and Anopheles funestus

Malaria remains one of the highest causes of morbidity and mortality, with 249 million cases and over 608,000 deaths in 2022. Insecticides, which target the Anopheles mosquito vector, are the primary method to control malaria. The widespread nature of resistance to the most important insecticide class, the pyrethroids, threatens the control of this disease. To reverse the stall in malaria control there is urgent need for new vector control tools, which necessitates understanding the molecular basis of pyrethroid resistance. In this study we utilised multi-omics data to identify uridine-diphosphate (UDP)-glycosyltransferases (UGTs) potentially involved in resistance across multiple Anopheles species. Phylogenetic analysis identifies sequence similarities between Anopheline UGTs and those involved in agricultural pesticide resistance to pyrethroids, pyrroles and spinosyns. Expression of five UGTs was characterised in An. gambiae and An. coluzzii to determine constitutive over-expression, induction, and tissue specificity. Furthermore, a UGT inhibitor, sulfinpyrazone, restored susceptibility to pyrethroids and DDT in An. gambiae, An. coluzzii, An. arabiensis and An. funestus, the major African malaria vectors. Taken together, this study provides clear association of UGTs with pyrethroid resistance as well as highlighting the potential use of sulfinpyrazone as a novel synergist for vector control.

Longitudinal survey of insecticide resistance in a village of central region of Burkina Faso reveals co-occurrence of 1014F, 1014S and 402L mutations in Anopheles coluzzii and Anopheles arabiensis

BACKGROUND

Pyrethroid resistance is one of the major threats for effectiveness of insecticide-treated bed nets (ITNs) in malaria vector control. Genotyping of mutations in the voltage gated sodium channel (VGSC) gene is widely used to easily assess the evolution and spread of pyrethroid target-site resistance among malaria vectors. L1014F and L1014S substitutions are the most common and best characterized VGSC mutations in major African malaria vector species of the Anopheles gambiae complex. Recently, an additional substitution involved in pyrethroid resistance, i.e. V402L, has been detected in Anopheles coluzzii from West Africa lacking any other resistance alleles at locus 1014. The evolution of target-site resistance mutations L1014F/S and V402L was monitored in An. coluzzii and Anopheles arabiensis specimens from a Burkina Faso village over a 10-year range after the massive ITN scale-up started in 2010.

METHODS

Anopheles coluzzii (N = 300) and An. arabiensis (N = 362) specimens collected both indoors and outdoors by different methods (pyrethrum spray catch, sticky resting box and human landing collections) in 2011, 2015 and 2020 at Goden village were genotyped by TaqMan assays and sequencing for the three target site resistance mutations; allele frequencies were statistically investigated over the years.

RESULTS

A divergent trend in resistant allele frequencies was observed in the two species: 1014F decreased in An. coluzzii (from 0.76 to 0.52) but increased in An. arabiensis (from 0.18 to 0.70); 1014S occurred only in An. arabiensis and slightly decreased over time (from 0.33 to 0.23); 402L increased in An. coluzzii (from 0.15 to 0.48) and was found for the first time in one An. arabiensis specimen. In 2020 the co-occurrence of different resistance alleles reached 43% in An. coluzzii (alleles 410L and 1014F) and 32% in An. arabiensis (alleles 1014F and 1014S).

CONCLUSIONS

Overall, an increasing level of target-site resistance was observed among the populations with only 1% of the two malaria vector species being wild type at both loci, 1014 and 402, in 2020. This, together with the co-occurrence of different mutations in the same specimens, calls for future investigations on the possible synergism between resistance alleles and their phenotype to implement local tailored intervention strategies.

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.

Mutations of voltage-gated sodium channel contribute to pyrethroid resistance in Panonychus citri

Insecticide resistance in Panonychus citri is a major obstacle to mite control in citrus orchards. Pyrethroid insecticides are continually used to control mites in China, although resistance to pyrethroids has evolved in some populations. Here, the resistance to the pyrethroid fenpropathrin was investigated and 7 out of 8 field-collected populations of P. citri exhibited a high level of resistance, ranging from 171-fold to 15 391-fold higher than the susceptible (SS) comparison strain. Three voltage-gated sodium channel (VGSC) mutations were identified in the tested populations: L1031V, F1747L, and F1751I. Amplicon sequencing was used to evaluate the frequency of these mutations in the 19 field populations. L1031V and F1751I were present in all populations at frequencies of 11.6%-82.1% and 0.5%-31.8%, respectively, whereas the F1747L mutation was only present in 12 populations from Chongqing, Sichuan, Guangxi, and Yunnan provinces. Introduction of these mutations singly or in combination into transgenic flies significantly increased their resistance to fenpropathrin and these flies also exhibited reduced mortality after exposure to the pyrethroids permethrin and β-cypermethrin. Panonychus citri VGSC homology modeling and ligand docking indicate that F1747 and F1751 form direct binding contacts with pyrethroids, which are lost with mutation, whereas L1031 mutation may diminish pyrethroid effects through an allosteric mechanism. Overall, the results provide molecular markers for monitoring pest resistance to pyrethroids and offer new insights into the basis of pyrethroid actions on sodium channels.

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.

AnoPrimer: Primer Design in malaria vectors informed by range-wide genomic variation

The major malaria mosquitoes, Anopheles gambiae s.l and Anopheles funestus, are some of the most studied organisms in medical research and also some of the most genetically diverse. When designing polymerase chain reaction (PCR) or hybridisation-based molecular assays, reliable primer and probe design is crucial. However, single nucleotide polymorphisms (SNPs) in primer binding sites can prevent primer binding, leading to null alleles, or bind suboptimally, leading to preferential amplification of specific alleles. Given the extreme genetic diversity of Anopheles mosquitoes, researchers need to consider this genetic variation when designing primers and probes to avoid amplification problems. In this note, we present a Python package, AnoPrimer, which exploits the Ag1000G and Af1000 datasets and allows users to rapidly design primers in An. gambiae or An. funestus, whilst summarising genetic variation in the primer binding sites and visualising the position of primer pairs. AnoPrimer allows the design of both genomic DNA and cDNA primers and hybridisation probes. By coupling this Python package with Google Colaboratory, AnoPrimer is an open and accessible platform for primer and probe design, hosted in the cloud for free. AnoPrimer is available here https://github.com/sanjaynagi/AnoPrimer and we hope it will be a useful resource for the community to design probe and primer sets that can be reliably deployed across the An. gambiae and funestus species ranges.

The Anopheles coluzzii range extends into Kenya: detection, insecticide resistance profiles and population genetic structure in relation to conspecific populations in West and Central Africa

BACKGROUND

Anopheles coluzzii is a primary vector of malaria found in West and Central Africa, but its presence has hitherto never been documented in Kenya. A thorough understanding of vector bionomics is important as it enables the implementation of targeted and effective vector control interventions. Malaria vector surveillance efforts in the country have tended to focus on historically known primary vectors. The current study sought to determine the taxonomic status of samples collected from five different malaria epidemiological zones in Kenya as well as describe the population genetic structure and insecticide resistance profiles in relation to other An. coluzzii populations.

METHODS

Mosquitoes were sampled as larvae from Busia, Kwale, Turkana, Kirinyaga and Kiambu counties, representing the range of malaria endemicities in Kenya, in 2019 and 2021 and emergent adults analysed using Whole Genome Sequencing (WGS) data processed in accordance with the Anopheles gambiae 1000 Genomes Project phase 3. Where available, historical samples from the same sites were included for WGS. Comparisons were made with An. coluzzii cohorts from West and Central Africa.

RESULTS

This study reports the detection of An. coluzzii for the first time in Kenya. The species was detected in Turkana County across all three time points from which samples were analyzed and its presence confirmed through taxonomic analysis. Additionally, there was a lack of strong population genetic differentiation between An. coluzzii from Kenya and those from the more northerly regions of West and Central Africa, suggesting they represent a connected extension to the known species range. Mutations associated with target-site resistance to DDT and pyrethroids and metabolic resistance to DDT were found at high frequencies up to 64%. The profile and frequencies of the variants observed were similar to An. coluzzii from West and Central Africa but the ace-1 mutation linked to organophosphate and carbamate resistance present in An. coluzzii from coastal West Africa was absent in Kenya.

CONCLUSIONS

These findings emphasize the need for the incorporation of genomics in comprehensive and routine vector surveillance to inform on the range of malaria vector species, and their insecticide resistance status to inform the choice of effective vector control approaches.

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.

Whole transcriptomic analysis reveals overexpression of salivary gland and cuticular proteins genes in insecticide-resistant Anopheles arabiensis from Western Kenya

BACKGROUND

Effective vector control is key to malaria prevention. However, this is now compromised by increased insecticide resistance due to continued reliance on insecticide-based control interventions. In Kenya, we have observed heterogenous resistance to pyrethroids and organophosphates in Anopheles arabiensis which is one of the most widespread malaria vectors in the country. We investigated the gene expression profiles of insecticide resistant An. arabiensis populations from Migori and Siaya counties in Western Kenya using RNA-Sequencing. Centers for Disease Control and Prevention (CDC) bottle assays were conducted using deltamethrin (DELTA), alphacypermethrin (ACYP) and pirimiphos-methyl (PMM) to determine the resistance status in both sites.

RESULTS

Mosquitoes from Migori had average mortalities of 91%, 92% and 58% while those from Siaya had 85%, 86%, and 30% when exposed to DELTA, ACYP and PMM, respectively. RNA-Seq analysis was done on pools of mosquitoes which survived exposure ('resistant'), mosquitoes that were not exposed, and the insecticide-susceptible An. arabiensis Dongola strain. Gene expression profiles of resistant mosquitoes from both Migori and Siaya showed an overexpression mainly of salivary gland proteins belonging to both the short and long form D7 genes, and cuticular proteins (including CPR9, CPR10, CPR15, CPR16). Additionally, the overexpression of detoxification genes including cytochrome P450s (CYP9M1, CYP325H1, CYP4C27, CYP9L1 and CYP307A1), 2 carboxylesterases and a glutathione-S-transferase (GSTE4) were also shared between DELTA, ACYP, and PMM survivors, pointing to potential contribution to cross resistance to both pyrethroid and organophosphate insecticides.

CONCLUSION

This study provides novel insights into the molecular basis of insecticide resistance in An. arabiensis in Western Kenya and suggests that salivary gland proteins and cuticular proteins are associated with resistance to multiple classes of insecticides.

Speciation within the Anopheles gambiae complex: high-throughput whole genome sequencing reveals evidence of a putative new cryptic taxon in 'far-west' Africa

The two main Afrotropical malaria vectors - Anopheles coluzzii and An. gambiae - are genetically distinct and reproductively isolated across West Africa. However, populations at the western extreme of their range are assigned as "intermediate" between the two species by whole genome sequence (WGS) data, and as hybrid forms by conventional molecular diagnostics. By exploiting WGS data from 1,190 specimens collected across west Africa via the Anopheles gambiae 1000 Genomes network, we identify a novel putative taxon in the far-west (provisionally named Bissau molecular form), which did not arise by admixture but rather originated at the same time as the split between An. coluzzii and An. gambiae. Intriguingly, these populations lack insecticide resistance mechanisms commonly observed in the two main species. These findings lead to a change of perspective on malaria vector species in the far-west region with potential for epidemiological implications, and a new challenge for genetic-based mosquito control approaches.

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.

The Anopheles coluzzii range extends into Kenya: Detection, insecticide resistance profiles and population genetic structure in relation to conspecific populations in West and Central Africa

Background

Anopheles coluzzii is a primary vector of malaria found in West and Central Africa, but its presence has hitherto never been documented in Kenya. A thorough understanding of vector bionomics is important as it enables the implementation of targeted and effective vector control interventions. Malaria vector surveillance efforts in the country have tended to focus on historically known primary vectors. In the current study, we sought to determine the taxonomic status of samples collected from five different malaria epidemiological zones in Kenya as well asdescribe the population genetic structure and insecticide resistance profiles in relation to other An. coluzzi populations.

Methods

Mosquitoes were sampled as larvae from Busia, Kwale, Turkana, Kirinyaga and Kiambu counties, representing the range of malaria endemicities in Kenya, in 2019 and 2021 and emergent adults analysed using Whole Genome Sequencing data processed in accordance with the Anopheles gambiae 1000 Genomes Project phase 3. Where available, historical samples from the same sites were included for WGS.

Results

This study reports the detection of Anopheles coluzzii for the first time in Kenya. The species was detected in Turkana County across all three time points sampled and its presence confirmed through taxonomic analysis. Additionally, we found a lack of strong population genetic differentiation between An. coluzzii from Kenya and those from the more northerly regions of West and Central Africa, suggesting they represent a connected extension to the known species range. Mutations associated with target-site resistance to DDT and pyrethroids and metabolic resistance to DDT were found at high frequencies of ~60%. The profile and frequencies of the variants observed were similar to An. coluzzii from West and Central Africa but the ace-1 mutation linked to organophosphate and carbamate resistance present in An. coluzzii from coastal West Africa was absent in Kenya.

Conclusions

These findings emphasise the need for the incorporation of genomics in comprehensive and routine vector surveillance to inform on the range of malaria vector species, and their insecticide resistance status to inform the choice of effective vector control approaches.

Malaria transmission potential of Anopheles gambiae s.l. in indoor residual spraying areas with clothianidin 50 WG in northern Benin

The study objective was to assess the frequency of the kdr-L995F and ace-1 G280S genetic mutations in Anopheles gambiae s.l. mosquitoes and examine their ability to transmit Plasmodium falciparum in areas where indoor residual spraying (IRS) was implemented with Clothianidin 50 WG. The study was conducted in six communes in the Alibori and Donga departments of which four were IRS-treated and two were untreated and served as control. Post-IRS monthly samples of adult mosquitoes were collected in study communes using human landing catches (HLC). An. gambiae s.l. specimens were processed to detect kdr-L995F and ace-1 G280S mutations via PCR as well as Plasmodium falciparum infectivity through CSP ELISA. Our data revealed a high and similar allelic frequency for the kdr-L995F mutation in both treated and control communes (79% vs. 77%, p = 0.14) whilst allelic frequency of the ace-1 G280S mutation was lower across the study area (2-3%, p = 0.58). The sporozoite rate was 2.6% and 2.4% respectively in treated and untreated communes (p = 0.751). No association was found between Plasmodium falciparum infection in Anopheles gambiae s.l. vectors and carriage of kdr-L995F and ace-1 G280S mutations regardless of genotypes. The study findings underline the need for an integrated approach to malaria control, combining different control methods to effectively target transmission. Regular monitoring of insecticide resistance and genetic mutations is essential to guide control strategies.

Genomic surveillance of Anopheles mosquitoes on the Bijagós Archipelago using custom targeted amplicon sequencing identifies mutations associated with insecticide resistance

BACKGROUND

Insecticide resistance is reducing the efficacy of vector control interventions, consequently threatening efforts to control vector-borne diseases, including malaria. Investigating the prevalence of molecular markers of resistance is a useful tool for monitoring the spread of insecticide resistance in disease vectors. The Bijagós Archipelago (Bijagós) in Guinea-Bissau is a region of stable malaria transmission where insecticide-treated nets are the mainstay for malaria control. However, the prevalence of molecular markers of insecticide resistance in malaria vectors is not well understood.

METHODS

A total of 214 Anopheles mosquitoes were analysed from 13 islands across the Bijagós. These mosquitoes were collected using CDC light traps in November 2019, during the peak malaria transmission season. High-throughput multiplex amplicon sequencing was used to investigate the prevalence of 17 different molecular markers associated with insecticide resistance in four genes: vgsc, rdl, ace1 and gste2.

RESULTS

Of the 17 screened mutations, four were identified in mosquitoes from the Bijagós: vgsc L995F (12.2%), N1570Y (6.2%) and A1746S (0.7%) and rdl A269G (1.1%). This study is the first to report the L995F knock-down resistance (kdr)-west allele in Anopheles melas on the Archipelago. An additional eight non-synonymous single-nucleotide polymorphisms were identified across the four genes which have not been described previously. The prevalences of the vgsc L995F and N1570Y mutations were higher on Bubaque Island than on the other islands in this study; Bubaque is the most populous island in the archipelago, with the greatest population mobility and connection to continental Guinea-Bissau.

CONCLUSIONS

This study provides the first surveillance data for genetic markers present in malaria vectors from islands across the Bijagós Archipelago. Overall prevalence of insecticide resistance mutations was found to be low. However, the identification of the vgsc L995F and N1570Y mutations associated with pyrethroid resistance warrants further monitoring. This is particularly important as the mainstay of malaria control on the islands is the use of pyrethroid insecticide-treated nets.

The role of P450 enzymes in malaria and other vector-borne infectious diseases

Vector-borne infectious diseases are still an important global health problem. Malaria is the most important among them, mainly pediatric, life-threatening disease. Malaria and other vector-borne disorders caused by parasites, bacteria, and viruses have a strong impact on public health and significant economic costs. Most vector-borne diseases could be prevented by vector control, with attention to the ecological and biodiversity conservation aspects. Chemical control with pesticides and insecticides is widely used as a measure of prevention although increasing resistance to insecticides is a serious issue in vector control. Metabolic resistance is the most common mechanism and poses a big challenge. Insect enzyme systems, including monooxygenase CYP P450 enzymes, are employed by vectors mainly to metabolize insecticides thus causing resistance. The discovery and application of natural specific inhibitors/blockers of vector P450 enzymes as synergists for commonly used pesticides will contribute to the "greening" of insecticides. Besides vector CYPs, host CYP enzymes could also be exploited to fight against vector-borne diseases: using mostly their detoxifying properties and involvement in the immune response. Here, we review published research data on P450 enzymes from all players in vector-borne infections, that is, pathogens, vectors, and hosts, regarding the potential role of CYPs in disease. We discuss strategies on how to exploit cytochromes P450 in vector-borne disease control.

Molecular surveillance of insecticide resistance in Phlebotomus argentipes targeted by indoor residual spraying for visceral leishmaniasis elimination in India

Molecular surveillance of resistance is an increasingly important part of vector borne disease control programmes that utilise insecticides. The visceral leishmaniasis (VL) elimination programme in India uses indoor residual spraying (IRS) with the pyrethroid, alpha-cypermethrin to control Phlebotomus argentipes the vector of Leishmania donovani, the causative agent of VL. Prior long-term use of DDT may have selected for knockdown resistance (kdr) mutants (1014F and S) at the shared DDT and pyrethroid target site, which are common in India and can also cause pyrethroid cross-resistance. We monitored the frequency of these marker mutations over five years from 2017-2021 in sentinel sites in eight districts of north-eastern India covered by IRS. Frequencies varied markedly among the districts, though finer scale variation, among villages within districts, was limited. A pronounced and highly significant increase in resistance-associated genotypes occurred between 2017 and 2018, but with relative stability thereafter, and some reversion toward more susceptible genotypes in 2021. Analyses linked IRS with mutant frequencies suggesting an advantage to more resistant genotypes, especially when pyrethroid was under-sprayed in IRS. However, this advantage did not translate into sustained allele frequency changes over the study period, potentially because of a relatively greater net advantage under field conditions for a wild-type/mutant genotype than projected from laboratory studies and/or high costs of the most resistant genotype. Further work is required to improve calibration of each 1014 genotype with resistance, preferably using operationally relevant measures. The lack of change in resistance mechanism over the span of the study period, coupled with available bioassay data suggesting susceptibility, suggests that resistance has yet to emerge despite intensive IRS. Nevertheless, the advantage of resistance-associated genotypes with IRS and under spraying, suggest that measures to continue monitoring and improvement of spray quality are vital, and consideration of future alternatives to pyrethroids for IRS would be advisable.

Detection of insecticide resistance markers in Anopheles funestus from the Democratic Republic of the Congo using a targeted amplicon sequencing panel

Vector control strategies have been successful in reducing the number of malaria cases and deaths globally, but the spread of insecticide resistance represents a significant threat to disease control. Insecticide resistance has been reported across Anopheles (An.) vector populations, including species within the An. funestus group. These mosquitoes are responsible for intense malaria transmission across sub-Saharan Africa, including in the Democratic Republic of the Congo (DRC), a country contributing > 12% of global malaria infections and mortality events. To support the continuous efficacy of vector control strategies, it is essential to monitor insecticide resistance using molecular surveillance tools. In this study, we developed an amplicon sequencing ("Amp-seq") approach targeting An. funestus, and using multiplex PCR, dual index barcoding, and next-generation sequencing for high throughput and low-cost applications. Using our Amp-seq approach, we screened 80 An. funestus field isolates from the DRC across a panel of nine genes with mutations linked to insecticide resistance (ace-1, CYP6P4, CYP6P9a, GSTe2, vgsc, and rdl) and mosquito speciation (cox-1, mtND5, and ITS2). Amongst the 18 non-synonymous mutations detected, was N485I, in the ace-1 gene associated with carbamate resistance. Overall, our panel represents an extendable and much-needed method for the molecular surveillance of insecticide resistance in An. funestus populations.

The voltage-gated sodium channel, para, limits Anopheles coluzzii vector competence in a microbiota dependent manner

The voltage-gated sodium channel, para, is a target of DDT and pyrethroid class insecticides. Single nucleotide mutations in para, called knockdown resistant or kdr, which contribute to resistance against DDT and pyrethroid insecticides, have been correlated with increased susceptibility of Anopheles to the human malaria parasite Plasmodium falciparum. However, a direct role of para activity on Plasmodium infection has not yet been established. Here, using RNA-mediated silencing, we provide in vivo direct evidence for the requirement of wild-type (wt) para function for insecticide activity of deltamethrin. Depletion of wt para, which is susceptible to insecticide, causes deltamethrin tolerance, indicating that insecticide-resistant kdr alleles are likely phenocopies of loss of para function. We then show that normal para activity in An. coluzzii limits Plasmodium infection prevalence for both P. falciparum and P. berghei. A transcriptomic analysis revealed that para activity does not modulate the expression of immune genes. However, loss of para function led to enteric dysbiosis with a significant increase in the total bacterial abundance, and we show that para function limiting Plasmodium infection is microbiota dependent. In the context of the bidirectional "enteric microbiota-brain" axis studied in mammals, these results pave the way for studying whether the activity of the nervous system could control Anopheles vector competence.

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.

Knockdown resistance allele L1014F introduced by CRISPR/Cas9 is not associated with altered vector competence of Anopheles gambiae for o'nyong nyong virus

Knockdown resistance (kdr) alleles conferring resistance to pyrethroid insecticides are widespread amongst vector populations. Previous research has suggested that these alleles are associated with changes in the vector competence of mosquitoes for arboviruses and Plasmodium, however non-target genetic differences between mosquito strains may have had a confounding effect. Here, to minimise genetic differences, the laboratory Anopheles gambiae Kisumu strain was compared to a CRISPR/Cas9 homozygous kdr L1014F mutant Kisumu-kdr line in order to examine associations with vector competence for o'nyong nyong virus (ONNV). Mosquitoes were infected using either blood feeds or intrathoracic microinjections. There were no significant differences in the prevalence of virus in mosquito body parts between kdr mutant and wildtype lines from either oral or intrathoracic injection routes. The ONNV titre was significantly higher in the legs of the wildtype strain at 7dpi following intrathoracic microinjection, but no other significant differences in viral titre were detected. ONNV was not detected in the saliva of mosquitoes from either strain. Our findings from per os infections suggest that the kdr L1014F allele is not associated with altered infection prevalence for ONNV, a key component of vector competence.

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.

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.

Strong Positive Selection in Aedes aegypti and the Rapid Evolution of Insecticide Resistance

Aedes aegypti vectors the pathogens that cause dengue, yellow fever, Zika virus, and chikungunya and is a serious threat to public health in tropical regions. Decades of work has illuminated many aspects of Ae. aegypti's biology and global population structure and has identified insecticide resistance genes; however, the size and repetitive nature of the Ae. aegypti genome have limited our ability to detect positive selection in this mosquito. Combining new whole genome sequences from Colombia with publicly available data from Africa and the Americas, we identify multiple strong candidate selective sweeps in Ae. aegypti, many of which overlap genes linked to or implicated in insecticide resistance. We examine the voltage-gated sodium channel gene in three American cohorts and find evidence for successive selective sweeps in Colombia. The most recent sweep encompasses an intermediate-frequency haplotype containing four candidate insecticide resistance mutations that are in near-perfect linkage disequilibrium with one another in the Colombian sample. We hypothesize that this haplotype may continue to rapidly increase in frequency and perhaps spread geographically in the coming years. These results extend our knowledge of how insecticide resistance has evolved in this species and add to a growing body of evidence suggesting that Ae. aegypti has an extensive genomic capacity to rapidly adapt to insecticide-based vector control.

The interplay between malaria vectors and human activity accounts for high residual malaria transmission in a Burkina Faso village with universal ITN coverage

BACKGROUND

Mosquito and human behaviour interaction is a key determinant of the maximum level of protection against malaria that can be provided by insecticide-treated nets (ITNs). Nevertheless, scant literature focuses on this interaction, overlooking a fundamental factor for efficient malaria control. This study aims to estimate malaria transmission risk in a Burkina Faso village by integrating vector biting rhythms with some key information about human habits.

METHODS

Indoor/outdoor human landing catches were conducted for 16 h (16:00-08:00) during 8 nights (September 2020) in Goden village. A survey about net usage and sleeping patterns was submitted to half the households (October-December 2020). A subsample of collected specimens of Anopheles gambiae sensu lato was molecularly processed for species identification, Plasmodium detection from heads-thoraxes and L1014F pyrethroid-resistance allele genotyping. Hourly mosquito abundance was statistically assessed by GLM/GAM, and the entomological inoculation rate (EIR) was corrected for the actual ITN usage retrieved from the questionnaire.

RESULTS

Malaria transmission was mainly driven by Anopheles coluzzii (68.7%) followed by A. arabiensis (26.2%). The overall sporozoite rate was 2% with L1014F estimated frequency of 0.68 (N = 1070 out of 15,201 A. gambiae s.l. collected). No major shift in mosquito biting rhythms in response to ITN or differences between indoor and outdoor catches were detected. Impressive high biting pressure (mean 30.3 mosquitoes/person/hour) was exerted from 20:00 to 06:00 with a peak at 4:00. Human survey revealed that nearly all inhabitants were awake before 20:00 and after 7:00 and at least 8.7% had no access to bednets. Adjusting for anthropological data, the EIR dropped from 6.7 to 1.2 infective bites/person/16 h. In a scenario of full net coverage and accounting only for the human sleeping patterns, the daily malaria transmission risk not targetable by ITNs was 0.69 infective bites.

CONCLUSIONS

The high mosquito densities and interplay between human/vector activities means that an estimated 10% of residual malaria transmission cannot be prevented by ITNs in the village. Locally tailored studies, like the current one, are essential to explore the heterogeneity of human exposure to infective bites and, consequently, to instruct the adoption of new vector control tools strengthening individual and community protection.

Widespread occurrence of copy number variants and fixation of pyrethroid target site resistance in Anopheles gambiae (s.l.) from southern Côte d'Ivoire

Resistance to pyrethroid and organophosphate insecticides in the malaria vector Anopheles gambiae (s.l.) is conferred by a variety of genetic mutations, including single nucleotide polymorphisms (SNPs) and copy number variants (CNVs). Knowledge of the distribution of these mutations in mosquito populations is a prerequisite for establishing better strategies for their management. In this study, a total of 755 Anopheles gambiae (s.l.) from southern Côte d'Ivoire were exposed to deltamethrin or pirimiphos-methyl insecticides and were screened to assess the distribution of SNPs and CNVs known or believed to confer resistance to one or other of the insecticide classes. Most individuals from the An. gambiae (s.l.) complex were identified by molecular tests as Anopheles coluzzii. Survival to deltamethrin (from 94% to 97%) was higher than to pirimiphos-methyl (from 10% to 49%). In An. gambiae (s.s.), the SNP in the Voltage Gated Sodium Channel (Vgsc) at the 995F locus (Vgsc-995F) was fixed, while other target site mutations were rare or absent (Vgsc-402L: 0%; Vgsc-1570Y: 0%, Acetylcholinesterase Acel-280S: 14%). In An. coluzzii, Vgsc-995F was the target site SNP found at highest frequency (65%) followed by other target site mutations (Vgsc-402L: 36%; Vgsc-1570Y: 0.33%; Acel-280S: 45%). The Vgsc-995S SNP was not present. The presence of the Ace1-280S SNP was found to be significantly linked to the presence of the Ace1-CNV, Ace1_AgDup. Significant association was found between the presence of the Ace1_AgDup and pirimiphos-methyl resistance in An. gambiae (s.s.) but not in An. coluzzii. The deletion Ace1_Del97 was found in one specimen of An. gambiae (s.s.). Four CNVs in the Cyp6aa/Cyp6p gene cluster, which contains genes of known importance for resistance, were detected in An. coluzzii, the most frequent being Dup 7 (42%) and Dup 14 (26%). While none of these individual CNV alleles were significantly associated with resistance, copy number in the Cyp6aa gene region in general was associated with increased resistance to deltamethrin. Elevated expression of Cyp6p3 was nearly associated with deltamethrin resistance, although there was no association of resistance with copy number. Use of alternative insecticides and control methods to arrest resistance spread in An. coluzzii populations is merited.

A novel V263I mutation in the glutamate-gated chloride channel of Plutella xylostella (L.) confers a high level of resistance to abamectin

The frequent and extensive use of insecticides leads to the evolution of insecticide resistance, which has become one of the constraints on global agricultural production. Avermectins are microbial-derived insecticides that target a wide number of insect pests, including the diamondback moth Plutella xylostella, an important global pest of brassicaceous vegetables. However, field populations of P. xylostella have evolved serious resistance to avermectins, including abamectin, thereby threatening the efficiency of these insecticides. In this study, a novel valine to isoleucine mutation (V263I) was identified in the glutamate-gated chloride channel (GluCl) of field P. xylostella populations, which showed different levels of resistance to abamectin. Electrophysiological analysis revealed that the V263I mutation significantly reduced the sensitivity of PxGluCl to abamectin by 6.9-fold. Genome-modified Drosophila melanogaster carrying the V263I mutation exhibited 27.1-fold resistance to abamectin. Then, a knockin strain (V263I-KI) of P. xylostella expressing the homozygous V263I mutation was successfully constructed using the CRISPR/Cas9. The V263I-KI had high resistance to abamectin (106.3-fold), but significantly reduced fecundity. In this study, the function of V263I mutation in PxGluCl was verified for the first time. These findings provide a more comprehensive understanding of abamectin resistance mechanisms and lay the foundation for providing a new molecular detection method for abamectin resistance monitoring.

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.

CRISPR-Mediated Cassette Exchange (CriMCE): A Method to Introduce and Isolate Precise Marker-Less Edits

The introduction of small unmarked edits to the genome of insects is essential to study the molecular underpinnings of important biological traits, such as resistance to insecticides and genetic control strategies. Advances in CRISPR genome engineering have made this possible, but prohibitively laborious for most laboratories due to low rates of editing and the lack of a selectable marker. To facilitate the generation and isolation of precise marker-less edits we have developed a two-step method based on CRISPR-mediated cassette exchange (CriMCE) of a marked placeholder for a variant of interest. This strategy can be used to introduce a wider range of potential edits compared with previous approaches while consolidating the workflow. We present proof-of-principle that CriMCE is a powerful tool by engineering three single nucleotide polymorphism variants into the genome of Anopheles gambiae, with 5-41 × higher rates of editing than homology-directed repair or prime editing.

High Levels of Admixture in Anopheles gambiae Populations from Côte d'Ivoire Revealed by Multilocus Genotyping

Anopheles coluzzii and An. gambiae-the two most recently radiated species of the An. gambiae complex and the major Afrotropical malaria vector species-are identified by markers in the X-centromeric IGS rDNA region. Putative IGS-hybrids are rarely found in the field, except in restricted areas where genomic studies have led to the hypothesis that the observed IGS-patterns are due to cryptic taxa rather than to hybridization between the two species. We investigated the genome-wide levels of admixture in two villages in Côte d'Ivoire where high levels of IGS-hybrids have been detected, confirming unparalleled high frequencies in the coastal village. Genotyping of 24 Ancestry Informative Markers (AIMs) along the three chromosomes produced discordant results between the IGS-marker and the multilocus genotype obtained for AIMs across the whole genome (29%) as well as AIMs on chromosome-X (considered to be fundamental for species reproductive isolation) only (21%). Results highlight a complicated pattern of admixture that deserves deeper genomic analyses to understand better possible underlying causes (from extensive processes of hybridization to the existence of different cryptic taxa), and stress the need of developing advanced diagnostics for An. coluzzii, An. gambiae and putative new taxa, instrumental for assessing taxon-specific epidemiological characters.

Analysis of the Genetic Variation of the Fruitless Gene within the Anopheles gambiae (Diptera: Culicidae) Complex Populations in Africa

Targeting genes involved in sexual determinism, for vector or pest control purposes, requires a better understanding of their polymorphism in natural populations in order to ensure a rapid spread of the construct. By using genomic data from An. gambiae s.l., we analyzed the genetic variation and the conservation score of the fru gene in 18 natural populations across Africa. A total of 34,339 SNPs were identified, including 3.11% non-synonymous segregating sites. Overall, the nucleotide diversity was low, and the Tajima’s D neutrality test was negative, indicating an excess of low frequency SNPs in the fru gene. The allelic frequencies of the non-synonymous SNPs were low (freq < 0.26), except for two SNPs identified at high frequencies (freq > 0.8) in the zinc-finger A and B protein domains. The conservation score was variable throughout the fru gene, with maximum values in the exonic regions compared to the intronic regions. These results showed a low genetic variation overall in the exonic regions, especially the male sex-specific exon and the BTB-exon 1 of the fru gene. These findings will facilitate the development of an effective gene drive construct targeting the fru gene that can rapidly spread without encountering resistance in wild populations.

Prospects for genomic surveillance for selection in schistosome parasites

Schistosomiasis is a neglected tropical disease caused by multiple parasitic Schistosoma species, and which impacts over 200 million people globally, mainly in low- and middle-income countries. Genomic surveillance to detect evidence for natural selection in schistosome populations represents an emerging and promising approach to identify and interpret schistosome responses to ongoing control efforts or other environmental factors. Here we review how genomic variation is used to detect selection, how these approaches have been applied to schistosomes, and how future studies to detect selection may be improved. We discuss the theory of genomic analyses to detect selection, identify experimental designs for such analyses, and review studies that have applied these approaches to schistosomes. We then consider the biological characteristics of schistosomes that are expected to respond to selection, particularly those that may be impacted by control programs. Examples include drug resistance, host specificity, and life history traits, and we review our current understanding of specific genes that underlie them in schistosomes. We also discuss how inherent features of schistosome reproduction and demography pose substantial challenges for effective identification of these traits and their genomic bases. We conclude by discussing how genomic surveillance for selection should be designed to improve understanding of schistosome biology, and how the parasite changes in response to selection.

Pyrethroid Resistance Situation across Different Eco-Epidemiological Settings in Cameroon

Rapid emergence and spread of pyrethroid resistance in Anopheles gambiae populations is among the main factors affecting malaria vector control in Cameroon, but there is still not enough data on the exact pyrethroid resistance status across Cameroon. The present study assessed pyrethroid resistance profile in different eco-epidemiological settings in Cameroon. Susceptibility bioassay tests were performed with F0 An. gambiae females aged three to five days. Mosquito susceptibility to both permethrin and deltamethrin was assessed. Species of the An. gambiae s.l. complex were identified using molecular diagnostic tools. Target site mutations conferring resistance were detected using Taqman assays. Quantitative reverse transcription-real-time PCR (qRT-PCR) 3-plex TaqMan^® assays were used for the quantification of detoxification genes implicated in pyrethroid resistance. An. gambiae, An. coluzzii and An. arabiensis were identified in the different settings. An. gambiae was dominant in Santchou, Kékem, Bélabo, Bertoua and Njombé, while An. coluzzii was abundant in Tibati and Kaélé. High frequencies of the kdr L1014F allele ranging from 43% to 100% were recorded in almost all sites. The L1014S kdr allele was detected at low frequency (4.10–10%) only in mosquito populations from Njombé and Tibati. The N1575Y mutation was recorded in Kaélé, Santchou, Tibati and Bertoua with a frequency varying from 2.10% to 11.70%. Six Cytochrome P450 genes (Cyp6p3, Cyp6m2, Cyp9k1, Cyp6p4, Cyp6z1, and Cyp4g16) were found to be overexpressed in at least one population. Analysis of cuticular hydrocarbon lipids indicated a significant increase in CHC content in mosquito populations from Kaélé and Njombé compared to Kékem, Bélabo and Bertoua populations. The study indicated high pyrethroid resistance across different ecological settings in Cameroon with different profile of resistance across the country. The present situation calls for further actions in order to mitigate the impact of insecticide resistance on vector control measures.

Temporal variation of high-level pyrethroid resistance in the major malaria vector Anopheles gambiae s.l. in Yaoundé, Cameroon, is mediated by target-site and metabolic resistance

Constant assessment of insecticide resistance levels is mandatory to implement adequate malaria control tools, but little information is available on the annual dynamics of resistance. We, therefore, monitored variations in resistance in Anopheles gambiae s.l. over four seasons during 2 years in two localities of Yaoundé: urban Etoa-Meki and peri-urban Nkolondom. Mosquitoes were collected seasonally at larval stage and reared to adults for insecticide susceptibility tests and molecular analysis of resistance mechanisms. Anopheles coluzzii was found in Etoa-Meki and An. gambiae in Nkolondom. Low mortalities to pyrethroids were observed (permethrin <10%, deltamethrin <21%), and resistance extended to 5× and 10× diagnostic doses, revealing a marked increase compared to previous studies. A seasonal variation in resistance was observed with the highest levels within dry seasons in Etoa-Meki and rainy seasons in Nkolondom. The 1014F kdr allele shows a high frequency (0.9), associated with overexpression of metabolic genes (Cyp6M2, Cyp6P4, Cyp9K1, Cyp6Z1, and Cyp6Z2) varying significantly seasonally. This study reveals an escalation in resistance to pyrethroids in Yaoundé's malaria vectors with seasonal variations. An adequate choice of the implementation period of punctual vector control actions according to the resistance profile will help to potentiate the desired effect and thus improve its efficiency.

High-throughput barcoding method for the genetic surveillance of insecticide resistance and species identification in Anopheles gambiae complex malaria vectors

Surveillance of malaria vector species and the monitoring of insecticide resistance are essential to inform malaria control strategies and support the reduction of infections and disease. Genetic barcoding of mosquitoes is a useful tool to assist the high-throughput surveillance of insecticide resistance, discriminate between sibling species and to detect the presence of Plasmodium infections. In this study, we combined multiplex PCR, custom designed dual indexing, and Illumina next generation sequencing for high throughput single nucleotide polymorphism (SNP)-profiling of four species from the Anopheles (An.) gambiae complex (An. gambiae sensu stricto, An. coluzzii, An. arabiensis and An. melas). By amplifying and sequencing only 14 genetic fragments (500 bp each), we were able to simultaneously detect Plasmodium infection; insecticide resistance-conferring SNPs in ace1, gste2, vgsc and rdl genes; the partial sequences of nuclear ribosomal internal transcribed spacers (ITS1 and ITS2) and intergenic spacers (IGS), Short INterspersed Elements (SINE), as well as mitochondrial genes (cox1 and nd4) for species identification and genetic diversity. Using this amplicon sequencing approach with the four selected An. gambiae complex species, we identified a total of 15 non-synonymous mutations in the insecticide target genes, including previously described mutations associated with resistance and two new mutations (F1525L in vgsc and D148E in gste2). Overall, we present a reliable and cost-effective high-throughput panel for surveillance of An. gambiae complex mosquitoes in malaria endemic regions.

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.

Active genetics comes alive: Exploring the broad applications of CRISPR-based selfish genetic elements (or gene-drives): Exploring the broad applications of CRISPR-based selfish genetic elements (or gene-drives)

Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-based "active genetic" elements developed in 2015 bypassed the fundamental rules of traditional genetics. Inherited in a super-Mendelian fashion, such selfish genetic entities offered a variety of potential applications including: gene-drives to disseminate gene cassettes carrying desired traits throughout insect populations to control disease vectors or pest species, allelic drives biasing inheritance of preferred allelic variants, neutralizing genetic elements to delete and replace or to halt the spread of gene-drives, split-drives with the core constituent Cas9 endonuclease and guide RNA (gRNA) components inserted at separate genomic locations to accelerate assembly of complex arrays of genetic traits or to gain genetic entry into novel organisms (vertebrates, plants, bacteria), and interhomolog based copying systems in somatic cells to develop tools for treating inherited or infectious diseases. Here, we summarize the substantial advances that have been made on all of these fronts and look forward to the next phase of this rapidly expanding and impactful field.

The evolution of larger size in high-altitude Drosophila melanogaster has a variable genetic architecture

Important uncertainties persist regarding the genetic architecture of adaptive trait evolution in natural populations, including the number of genetic variants involved, whether they are drawn from standing genetic variation, and whether directional selection drives them to complete fixation. Here, we take advantage of a unique natural population of Drosophila melanogaster from the Ethiopian highlands, which has evolved larger body size than any other known population of this species. We apply a bulk segregant quantitative trait locus mapping approach to 4 unique crosses between highland Ethiopian and lowland Zambian populations for both thorax length and wing length. Results indicated a persistently variable genetic basis for these evolved traits (with largely distinct sets of quantitative trait loci for each cross), and at least a moderately polygenic architecture with relatively strong effects present. We complemented these mapping experiments with population genetic analyses of quantitative trait locus regions and gene ontology enrichment analysis, generating strong hypotheses for specific genes and functional processes that may have contributed to these adaptive trait changes. Finally, we find that the genetic architectures indicated by our quantitative trait locus mapping results for size traits mirror those from similar experiments on other recently evolved traits in this species. Collectively, these studies suggest a recurring pattern of polygenic adaptation in this species, in which causative variants do not approach fixation and moderately strong effect loci are present.

In vivo functional validation of the V402L voltage gated sodium channel mutation in the malaria vector An. gambiae

BACKGROUND

Pyrethroids are the most widely used insecticides for the control of malaria transmitting Anopheles gambiae mosquitoes and rapid increase in resistance to this insecticide class is of major concern. Pyrethroids target the Voltage Gated Sodium Channels (VGSCs), that have a key role in the normal function of the mosquitoes' nervous system. VGSC mutations L995F and L995S have long been associated with pyrethroid resistance and screening for their presence is routine in insecticide resistance management programs. Recently, a VGSC haplotype containing two amino acid substitutions associated with resistance in other species, V402L and I1527T, was identified. These two VGSC mutations are found in tight linkage and are mutually exclusive to the classical L995F/S mutations.

RESULTS

We identify the presence of the V402L-I1527T haplotype in resistant An. coluzzii colonized strains and in field populations from Burkina Faso, at frequencies higher than previously reported; in some cases almost reaching fixation. Functional validation of V402L in insecticide resistance using a CRISPR/Cas9 genome modified line showed that it confers reduced mortality after exposure to all tested pyrethroids and DDT, but at lower levels compared to L995F. In contrast to L995F however, no fitness costs were identified for mosquitoes carrying V402L under laboratory conditions.

CONCLUSION

The V402L substitution confers pyrethroid resistance in An. gambiae in the absence of any other VGSC substitution and/or alternative resistance mechanisms. The lower fitness cost associated with this kdr mutation may provide a selective advantage over the classical kdr in some settings and genotyping at this locus should be added in the list of resistant alleles for routine screening.

Sympatric Populations of the Anopheles gambiae Complex in Southwest Burkina Faso Evolve Multiple Diverse Resistance Mechanisms in Response to Intense Selection Pressure with Pyrethroids

Simple Summary

Targeting mosquitoes with insecticides is one of the most effective methods to prevent malaria transmission. Although numbers of malaria cases have declined substantially this century, this pattern is not universal and Burkina Faso has one of the highest burdens of malaria; it is also a hotspot for the evolution of insecticide resistance in malaria vectors. We have established laboratory colonies from multiple species within the An. gambiae complex, the most efficient group of malaria vectors in the world, from larval collections in southwest Burkina Faso. Using bioassays with different insecticides widely used to control public health pests, we provide a profile of insecticide resistance in each of these colonies and, using molecular tools, reveal the genetic changes underpinning this resistance. We show that, whilst many resistance mechanisms are shared between species, there are some important differences which may affect resistance to current and future insecticide classes. The complexity, and diversity of resistance mechanisms highlights the importance of screening any potential new insecticide intended for use in malaria control against a wide range of populations. These stable laboratory colonies provide a valuable resource for insecticide discovery, and for further studies on the evolution and dispersal of insecticide resistance within and between species.

Abstract

Pyrethroid resistance in the Anopheles vectors of malaria is driving an urgent search for new insecticides that can be used in proven vector control tools such as insecticide treated nets (ITNs). Screening for potential new insecticides requires access to stable colonies of the predominant vector species that contain the major pyrethroid resistance mechanisms circulating in wild populations. Southwest Burkina Faso is an apparent hotspot for the emergence of pyrethroid resistance in species of the Anopheles gambiae complex. We established stable colonies from larval collections across this region and characterised the resistance phenotype and underpinning genetic mechanisms. Three additional colonies were successfully established (1 An. coluzzii, 1 An. gambiae and 1 An. arabiensis) to add to the 2 An. coluzzii colonies already established from this region; all 5 strains are highly resistant to pyrethroids. Synergism assays found that piperonyl butoxide (PBO) exposure was unable to fully restore susceptibility although exposure to a commercial ITN containing PBO resulted in 100% mortality. All colonies contained resistant alleles of the voltage gated sodium channel but with differing proportions of alternative resistant haplotypes. RNAseq data confirmed the role of P450s, with CYP6P3 and CYP6Z2 elevated in all 5 strains, and identified many other resistance mechanisms, some found across strains, others unique to a particular species. These strains represent an important resource for insecticide discovery and provide further insights into the complex genetic changes driving pyrethroid resistance.

Modelling spatiotemporal trends in the frequency of genetic mutations conferring insecticide target-site resistance in African mosquito malaria vector species

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.

Multi-insecticide resistant malaria vectors in the field remain susceptible to malathion, despite the presence of Ace1 point mutations

Insecticide resistance in Anopheles mosquitoes is seriously threatening the success of insecticide-based malaria vector control. Surveillance of insecticide resistance in mosquito populations and identifying the underlying mechanisms enables optimisation of vector control strategies. Here, we investigated the molecular mechanisms of insecticide resistance in three Anopheles coluzzii field populations from southern Côte d’Ivoire, including Agboville, Dabou and Tiassalé. All three populations were resistant to bendiocarb, deltamethrin and DDT, but not or only very weakly resistant to malathion. The absence of malathion resistance is an unexpected result because we found the acetylcholinesterase mutation Ace1-G280S at high frequencies, which would typically confer cross-resistance to carbamates and organophosphates, including malathion. Notably, Tiassalé was the most susceptible population to malathion while being the most resistant one to the pyrethroid deltamethrin. The resistance ratio to deltamethrin between Tiassalé and the laboratory reference colony was 1,800 fold. By sequencing the transcriptome of individual mosquitoes, we found numerous cytochrome P450-dependent monooxygenases – including CYP6M2, CYP6P2, CYP6P3, CYP6P4 and CYP6P5 – overexpressed in all three field populations. This could be an indication for negative cross-resistance caused by overexpression of pyrethroid-detoxifying cytochrome P450s that may activate pro-insecticides, thereby increasing malathion susceptibility. In addition to the P450s, we found several overexpressed carboxylesterases, glutathione S-transferases and other candidates putatively involved in insecticide resistance.

Insecticide-based mosquito control has saved millions of lives from malaria and other vector-borne diseases. However, the emergence and increase of insecticide resistant Anopheles populations seriously threaten to derail malaria control programmes. Surveillance of insecticide resistance and understanding the underlying molecular mechanisms are key for choosing effective vector control strategies. Here, we characterised the degree and mechanisms of resistance in three malaria vector populations from Côte d’Ivoire. Our key finding was that these multi-insecticide resistant malaria vectors largely remained susceptible to malathion, despite the presence of a mutation in the target enzyme of this organophosphate insecticide that would typically confer resistance. Intriguingly, we found overexpression of metabolic P450 enzymes that are known to detoxify insecticides and activate pro-insecticides such as malathion. It is highly probable that, here, we observed P450-mediated negative cross-resistance for the first time in Anopheles field populations. Negative cross-resistance merits further investigation as advantage could be taken of this phenomenon in the fight against multi-resistant malaria vectors.

Reversing insecticide resistance with allelic-drive in Drosophila melanogaster

A recurring target-site mutation identified in various pests and disease vectors alters the voltage gated sodium channel (vgsc) gene (often referred to as knockdown resistance or kdr) to confer resistance to commonly used insecticides, pyrethroids and DDT. The ubiquity of kdr mutations poses a major global threat to the continued use of insecticides as a means for vector control. In this study, we generate common kdr mutations in isogenic laboratory Drosophila strains using CRISPR/Cas9 editing. We identify differential sensitivities to permethrin and DDT versus deltamethrin among these mutants as well as contrasting physiological consequences of two different kdr mutations. Importantly, we apply a CRISPR-based allelic-drive to replace a resistant kdr mutation with a susceptible wild-type counterpart in population cages. This successful proof-of-principle opens-up numerous possibilities including targeted reversion of insecticide-resistant populations to a native susceptible state or replacement of malaria transmitting mosquitoes with those bearing naturally occurring parasite resistant alleles.

Insecticide resistance (IR) poses a major global health challenge. Here, the authors generate common IR mutations in laboratory Drosophila strains and use a CRISPR-based allelic-drive to replace an IR allele with a susceptible wild-type counterpart, providing a potent new tool for vector control.