Metabolic and target-site mechanisms combine to confer strong DDT resistance in Anopheles gambiae

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.

Structural and Thermodynamic Insights into Dimerization Interfaces of Drosophila Glutathione Transferases

This study presents a comprehensive analysis of the dimerization interfaces of fly GSTs through sequence alignment. Our investigation revealed GSTE1 as a particularly intriguing target, providing valuable insights into the variations within Delta and Epsilon GST interfaces. The X-ray structure of GSTE1 was determined, unveiling remarkable thermal stability and a distinctive dimerization interface. Utilizing circular dichroism, we assessed the thermal stability of GSTE1 and other Drosophila GSTs with resolved X-ray structures. The subsequent examination of GST dimer stability correlated with the dimerization interface supported by findings from X-ray structural analysis and thermal stability measurements. Our discussion extends to the broader context of GST dimer interfaces, offering a generalized perspective on their stability. This research enhances our understanding of the structural and thermodynamic aspects of GST dimerization, contributing valuable insights to the field.

Uncovering the genetic diversity in Aedes aegypti insecticide resistance genes through global comparative genomics

Aedes aegypti is vector of many arboviruses including Zika, dengue, yellow fever, West Nile, and Chikungunya. Its control efforts are hampered by widespread insecticide resistance reported in the Americas and Asia, while data from Africa is more limited. Here we use publicly available 729 Ae. aegypti whole-genome sequencing samples from 15 countries, including nine in Africa, to investigate the genetic diversity in four insecticide resistance linked genes: ace-1, GSTe2, rdl and vgsc. Apart from vgsc, the other genes have been less investigated in Ae. aegypti, and almost no genetic diversity information is available. Among the four genes, we identified 1,829 genetic variants including 474 non-synonymous substitutions, some of which have been previously documented, as well as putative copy number variations in GSTe2 and vgsc. Global insecticide resistance phenotypic data demonstrated variable resistance in geographic areas with resistant genotypes. Overall, our work provides the first global catalogue and geographic distribution of known and new amino-acid mutations and duplications that can be used to guide the identification of resistance drivers in Ae. aegypti and thereby support monitoring efforts and strategies for vector control.

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.

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.

A Novel Approach Utilizing Domain Adversarial Neural Networks for the Detection and Classification of Selective Sweeps

The identification and classification of selective sweeps are of great significance for improving the understanding of biological evolution and exploring opportunities for precision medicine and genetic improvement. Here, a domain adaptation sweep detection and classification (DASDC) method is presented to balance the alignment of two domains and the classification performance through a domain-adversarial neural network and its adversarial learning modules. DASDC effectively addresses the issue of mismatch between training data and real genomic data in deep learning models, leading to a significant improvement in its generalization capability, prediction robustness, and accuracy. The DASDC method demonstrates improved identification performance compared to existing methods and excels in classification performance, particularly in scenarios where there is a mismatch between application data and training data. The successful implementation of DASDC in real data of three distinct species highlights its potential as a useful tool for identifying crucial functional genes and investigating adaptive evolutionary mechanisms, particularly with the increasing availability of genomic data.

Preventing multi-resistance: New insights for managing fungal adaptation

Sustainable crop protection is vital for food security, yet it is under threat due to the adaptation of a diverse and evolving pathogen population. Resistance can be managed by maximising the diversity of selection pressure through dose variation and the spatial and temporal combination of active ingredients. This study explores the interplay between operational drivers for maximising the sustainability of management strategies in relation to the resistance status of fungal populations. We applied an experimental evolution approach to three artificial populations of Zymoseptoria tritici, an economically significant wheat pathogen, each differing in initial resistance status. Our findings reveal that diversified selection pressure curtails the selection of resistance in naïve populations and those with low frequencies of single resistance. Increasing the number of modes of action most effectively delays resistance development, surpassing the increase in the number of fungicides, fungicide choice based on resistance risk, and temporal variation in fungicide exposure. However, this approach favours generalism in the evolved populations. The prior presence of multiple resistant isolates and their subsequent selection in populations override the effects of diversity in management strategies, thereby invalidating any universal ranking. Therefore, the initial resistance composition must be specifically considered in sustainable resistance management to address real-world field situations.

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.

Aedes aegypti CCEae3A carboxylase expression confers carbamate, organophosphate and limited pyrethroid resistance in a model transgenic mosquito

Insecticide resistance is a serious threat to our ability to control mosquito vectors which transmit pathogens including malaria parasites and arboviruses. Understanding the underlying mechanisms is an essential first step in tackling the challenges presented by resistance. This study aimed to functionally characterise the carboxylesterase, CCEae3A, the elevated expression of which has been implicated in temephos resistance in Aedes aegypti and Aedes albopictus larvae. Using our GAL4/UAS expression system, already established in insecticide-sensitive Anopheles gambiae mosquitoes, we produced transgenic An. gambiae mosquitoes that express an Ae. aegypti CCEae3A ubiquitously. This new transgenic line permits examination of CCEae3A expression in a background in which there is not a clear orthologue in Vectorbase and allows comparison with existing An. gambiae GAL4-UAS lines. Insecticide resistance profiling of these transgenic An. gambiae larvae indicated significant increases in resistance ratio for three organophosphate insecticides, temephos (6), chloropyriphos (6.6) and fenthion (3.2) when compared to the parental strain. Cross resistance to adulticides from three major insecticide classes: organophosphates (malathion, fenitrothion and pirimiphos methyl), carbamates (bendiocarb and propoxur) and pyrethroid (alpha-cypermethrin) was also detected. Resistance to certain organophosphates and carbamates validates conclusions drawn from previous expression and phenotypic data. However, detection of resistance to pirimiphos methyl and alphacypermethrin has not previously been formally associated with CCEae3A, despite occurring in Ae. aegypti strains where this gene was upregulated. Our findings highlight the importance of characterising individual resistance mechanisms, thereby ensuring accurate information is used to guide future vector control strategies.

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

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

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.

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.

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.

Recombinant expression and characterization of GSTd3 from a resistant population of Anopheles arabiensis and comparison of DDTase activity with GSTe2

The development of insecticide resistance in malaria vectors is a challenge for the global effort to control and eradicate malaria. Glutathione S-transferases (GSTs) are multifunctional enzymes involved in the detoxification of many classes of insecticides. For mosquitoes, it is known that overexpression of an epsilon GST, GSTe2, confers resistance towards DDT and pyrethroids. In addition to GSTe2, consistent overexpression of a delta class GST, GSTd3, has been observed in insecticide resistant populations of different malaria vector species. However, the functional role of GSTd3 towards DDT resistance has not yet been investigated. Here, we recombinantly expressed both GSTe2 and GSTd3 from Anopheles arabiensis and compared their metabolic activities against DDT. Both AaGSTd3 and AaGSTe2 exhibited CDNB-conjugating and glutathione peroxidase activity and DDT metabolism was observed for both GSTs. However, the DDT dehydrochlorinase activity exhibited by AaGSTe2 was much higher than for AaGSTd3, and AaGSTe2 was also able to eliminate DDE although the metabolite could not be identified. Molecular modeling revealed subtle differences in the binding pocket of both enzymes and a better fit of DDT within the H-site of AaGSTe2. The overexpression but much lower DDT metabolic activity of AaGSTd3, might suggest that AaGSTd3 sequesters DDT. These findings highlight the complexity of insecticide resistance in the major malaria vectors and the difficulties associated with control of the vectors using DDT, which is still used for indoor residual spraying.

Detection of a reduced susceptibility to chlorfenapyr in the malaria vector Anopheles gambiae contrasts with full susceptibility in Anopheles funestus across Africa

New insecticides have recently been produced to help control pyrethroid-resistant malaria vectors including the pyrrole, chlorfenapyr. Monitoring the susceptibility of mosquito populations against this new product and potential cross-resistance with current insecticides is vital for better resistance management. In this study, we assessed the resistance status of the major malaria vectors Anopheles gambiae and Anopheles funestus to chlorfenapyr across Africa and explored potential cross-resistance with known pyrethroid resistance markers. Efficacy of chlorfenapyr 100 µg/ml against An. gambiae and An. funestus from five Cameroonian locations, the Democratic Republic of Congo, Ghana, Uganda, and Malawi was assessed using CDC bottle assays. Synergist assays were performed with PBO (4%), DEM (8%) and DEF (0.25%) and several pyrethroid-resistant markers were genotyped in both species to assess potential cross-resistance between pyrethroids and chlorfenapyr. Resistance to chlorfenapyr was detected in An. gambiae populations from DRC (Kinshasa) (mortality rate: 64.3 ± 7.1%) Ghana (Obuasi) (65.9 ± 7.4%), Cameroon (Mangoum; 75.2 ± 7.7% and Nkolondom; 86.1 ± 7.4). In contrast, all An. funestus populations were fully susceptible. A negative association was observed between the L1014F-kdr mutation and chlorfenapyr resistance with a greater frequency of homozygote resistant mosquitoes among the dead mosquitoes after exposure compared to alive (OR 0.5; P = 0.02) whereas no association was found between GSTe2 (I114T in An. gambiae; L119F in An. funestus) and resistance to chlorfenapyr. A significant increase of mortality to chlorfenapyr 10 µg/ml was observed in An. funestus after to PBO, DEM and DEF whereas a trend for a decreased mortality was observed in An. gambiae after PBO pre-exposure. This study reveals a greater risk of chlorfenapyr resistance in An. gambiae populations than in An. funestus. However, the higher susceptibility in kdr-resistant mosquitoes points to higher efficacy of chlorfenapyr against the widespread kdr-based pyrethroid resistance.

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.

Evidence of High Frequencies of Insecticide Resistance Mutations in Aedes aegypti (Culicidae) Mosquitoes in Urban Accra, Ghana: Implications for Insecticide-based Vector Control of Aedes-borne Arboviral Diseases

The most widespread arboviral diseases such as Dengue, Chikungunya, and Zika are transmitted mainly by Aedes mosquitoes. Due to the lack of effective therapeutics for most of these diseases, vector control remains the most effective preventative and control measure. This study investigated and compared the species composition, insecticide susceptibility, and resistance mechanisms in Aedes mosquito populations from a forest reserve converted to an eco-park and a peri-domestic sites in urban Accra, Ghana. Immature Aedes were sampled from the study sites, raised to adults, and exposed to deltamethrin, permethrin, DDT, fenitrothion, bendiocarb, permethrin + PBO, and deltamethrin + PBO using WHO tube assays. Melting curve analyses were performed for F1536C, V1016I, and V410L genetic mutations in surviving and dead mosquitoes following exposure to deltamethrin and permethrin. Microplate assay was used to access enzyme activity levels in adult mosquitoes from both populations. Aedes aegypti was found to be the dominant species from both study populations. The susceptibility test results revealed a high frequency of resistance to all the insecticides except fenitrothion. F1534C mutations were observed in 100% and 97% of mosquitoes from the peri-domestic and forest population, respectively but were associated with pyrethroid resistance only in the forest population (P < 0.0001). For the first time in Aedes mosquitoes in Ghana, we report the existence V410L mutations, mostly under selection only in the forest population (HWE P < 0.0001) and conclude that Aedes vectors in urban Accra have developed resistance to many commonly used insecticides. This information is important for the formulation of vector control strategies for Aedes control in Ghana.

Tandem duplication of a genomic region encoding glutathione S-transferase epsilon-2 and -4 genes in DDT-resistant Anopheles stephensi strain from India

The glutathione S-transferases (GST) genes are a multigene family of enzymes involved in the metabolism of endogenous and xenobiotic compounds by catalysing the conjugation of the reduced form of glutathione to the substrate. The epsilon class of GST (GSTe), unique to arthropods, is known to be involved in the detoxification process of several classes of insecticides, and GSTe2 in particular is known to have DDT dehydrochlorinase activity. This communication reports a tandem duplication of a genomic region encoding GSTe2 and GSTe4 genes in a laboratory-colonized DDT-resistant Anopheles stephensi. We identified duplication breakpoints and the organization of gene duplication through Sanger sequencing performed on long-PCR products. Manual annotation of sequences revealed a tandemly-arrayed duplication of a 3.62 kb segment of GST epsilon gene clusters comprised of five genes: a partial GSTe1, GSTe2, GSTe2-pseudogene, GSTe4 and partial GSTe5, interconnected by a conserved 2.42 kb DNA insert segment major part of which is homologous to a genomic region located on a different chromosome. The tandemly duplicated array contained a total of two GSTe2 and three GSTe4 functional paralog genes. Read-depth coverage and split-read analysis of Illumina-based whole-genome sequence reads confirmed the presence of duplication in the corresponding region of the genome. The increased gene dose in mosquitoes as a result of the GSTe gene-duplication may be an adaptive process to increase levels of detoxifying enzymes to counter insecticide pressure.

Evolution of the Ace-1 and Gste2 Mutations and Their Potential Impact on the Use of Carbamate and Organophosphates in IRS for Controlling Anopheles gambiae s.l., the Major Malaria Mosquito in Senegal

Widespread of insecticide resistance amongst the species of the Anopheles gambiae complex continues to threaten vector control in Senegal. In this study, we investigated the presence and evolution of the Ace-1 and Gste2 resistance genes in natural populations of Anopheles gambiae s.l., the main malaria vector in Senegal. Using historical samples collected from ten sentinel health districts, this study focused on three different years (2013, 2017, and 2018) marking the periods of shift between the main public health insecticides families (pyrethroids, carbamates, organophosphates) used in IRS to track back the evolutionary history of the resistance mutations on the Ace-1 and Gste2 loci. The results revealed the presence of four members of the Anopheles gambiae complex, with the predominance of An. arabiensis followed by An. gambiae, An. coluzzii, and An. gambiae–coluzzii hybrids. The Ace-1 mutation was only detected in An. gambiae and An. gambiae–coluzzii hybrids at low frequencies varying between 0.006 and 0.02, while the Gste2 mutation was found in all the species with a frequency ranging between 0.02 and 0.25. The Ace-1 and Gste2 genes were highly diversified with twenty-two and thirty-one different haplotypes, respectively. The neutrality tests on each gene indicated a negative Tajima’s D, suggesting the abundance of rare alleles. The presence and spread of the Ace-1 and Gste2 resistance mutations represent a serious threat to of the effectiveness and the sustainability of IRS-based interventions using carbamates or organophosphates to manage the widespread pyrethroids resistance in Senegal. These data are of the highest importance to support the NMCP for evidence-based vector control interventions selection and targeting.

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.

Characterization of the glutathione S-transferase genes in the sand flies Phlebotomus papatasi and Lutzomyia longipalpis shows expansion of the novel glutathione S-transferase xi (X) class

Leishmaniasis control often relies upon insecticidal control of phlebotomine sandfly vector populations. Such methods are vulnerable to the evolution of insecticide resistance via a range of molecular mechanisms. There is evidence that two major resistance mechanisms, target site insensitivity and metabolic resistance, have evolved in some sandfly populations and further genetic characterization of resistance would be useful to understand and combat it. To facilitate the study of the mechanisms of metabolic resistance, here we improved the annotation and characterized a major detoxification gene family, the glutathione-s-transferases (GST), in the genomes of two sand fly species: Phlebotomus papatasi and Lutzomyia longipalpis. The compositions of the GST gene family differ markedly from those of Aedes and Anopheles mosquitoes. Most strikingly, the xi (X) class of GSTs appears to have expanded in both sand fly genomes. Our results provide a basis for further studies of metabolic resistance mechanisms in these important disease vector species.

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.

Multi-omics analysis identifies a CYP9K1 haplotype conferring pyrethroid resistance in the malaria vector Anopheles funestus in East Africa

Metabolic resistance to pyrethroids is a menace to the continued effectiveness of malaria vector controls. Its molecular basis is complex and varies geographically across Africa. Here, we used a multi‐omics approach, followed‐up with functional validation to show that a directionally selected haplotype of a cytochrome P450, CYP9K1 is a major driver of resistance in Anopheles funestus. A PoolSeq GWAS using mosquitoes alive and dead after permethrin exposure, from Malawi and Cameroon, detected candidate genomic regions, but lacked consistency across replicates. Targeted sequencing of candidate resistance genes detected several SNPs associated with known pyrethroid resistance QTLs. The most significant SNPs were in the cytochrome P450 CYP304B1 (Cameroon), CYP315A1 (Uganda) and the ABC transporter gene ABCG4 (Malawi). However, when comparing field resistant mosquitoes to laboratory susceptible, the pyrethroid resistance locus rp1 and SNPs around the ABC transporter ABCG4 were consistently significant, except for Uganda where SNPs in the P450 CYP9K1 was markedly significant. In vitro heterologous metabolism assays with recombinant CYP9K1 revealed that it metabolises type II pyrethroid (deltamethrin; 64% depletion) but not type I (permethrin; 0%), while moderately metabolising DDT (17%). CYP9K1 exhibited reduced genetic diversity in Uganda underlying an extensive selective sweep. Furthermore, a glycine to alanine (G454A) amino acid change in CYP9K1 was fixed in Ugandan mosquitoes but not in other An. funestus populations. This study sheds further light on the evolution of metabolic resistance in a major malaria vector by implicating more genes and variants that can be used to design field‐applicable markers to better track resistance Africa‐wide.

Comparative study of the effect of solvents on the efficacy of neonicotinoid insecticides against malaria vector populations across Africa

Background

New insecticides with a novel mode of action such as neonicotinoids have recently been recommended for public health by WHO. Resistance monitoring of such novel insecticides requires a robust protocol to monitor the development of resistance in natural populations. In this study, we comparatively used three different solvents to assess the susceptibility of malaria vectors to neonicotinoids across Africa.

Methods

Mosquitoes were collected from May to July 2021 from three agricultural settings in Cameroon (Njombe-Penja, Nkolondom, and Mangoum), the Democratic Republic of Congo (Ndjili-Brasserie), Ghana (Obuasi), and Uganda (Mayuge). Using the CDC bottle test, we compared the effect of three different solvents (ethanol, acetone, MERO) on the efficacy of neonicotinoids against Anopheles gambiae s.l. In addition, TaqMan assays were used to genotype key pyrethroid-resistant markers in An. gambiae and odds ratio based on Fisher exact test were used to evaluate potential cross-resistance between pyrethroids and clothianidin.

Results

Lower mortality was observed when using absolute ethanol or acetone alone as solvent for clothianidin (11.4‒51.9% mortality in Nkolondom, 31.7‒48.2% in Mangoum, 34.6‒56.1% in Mayuge, 39.4‒45.6% in Obuasi, 83.7‒89.3% in Congo and 71.1‒95.9% in Njombe pendja) compared to acetone + MERO for which 100% mortality were observed for all the populations. Similar observations were done for imidacloprid and acetamiprid. Synergist assays (PBO, DEM and DEF) with clothianidin revealed a significant increase of mortality suggesting that metabolic resistance mechanisms are contributing to the reduced susceptibility. A negative association was observed between the L1014F-kdr mutation and clothianidin resistance with a greater frequency of homozygote resistant mosquitoes among the dead than among survivors (OR = 0.5; P = 0.02). However, the I114T-GSTe2 was in contrast significantly associated with a greater ability to survive clothianidin with a higher frequency of homozygote resistant among survivors than other genotypes (OR = 2.10; P = 0.013).

Conclusions

This study revealed a contrasted susceptibility pattern depending on the solvents with ethanol/acetone resulting to lower mortality, thus possibly overestimating resistance, whereas the MERO consistently showed a greater efficacy of neonicotinoids but it could prevent to detect early resistance development. Therefore, we recommend monitoring the susceptibility using both acetone alone and acetone + MERO (4 µg/ml for clothianidin) to capture the accurate resistance profile of the mosquito populations.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s40249-022-00962-4.

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.

Expression of pyrethroid metabolizing P450 enzymes characterizes highly resistant Anopheles vector species targeted by successful deployment of PBO-treated bednets in Tanzania

Long lasting insecticidal nets (LLINs) are a proven tool to reduce malaria transmission, but in Africa efficacy is being reduced by pyrethroid resistance in the major vectors. A previous study that was conducted in Muleba district, Tanzania indicated possible involvement of cytochrome P450 monooxygenases in a pyrethroid resistance in An. gambiae population where pre-exposure to piperonyl butoxide (PBO) followed by permethrin exposure in CDC bottle bioassays led to partial restoration of susceptibility. PBO is a synergist that can block pyrethroid-metabolizing enzymes in a mosquito. Insecticide resistance profiles and underlying mechanisms were investigated in Anopheles gambiae and An. funestus from Muleba during a cluster randomized trial. Diagnostic dose bioassays using permethrin, together with intensity assays, suggest pyrethroid resistance that is both strong and very common, but not extreme. Transcriptomic analysis found multiple P450 genes over expressed including CYP6M2, CYP6Z3, CYP6P3, CYP6P4, CYP6AA1 and CYP9K1 in An. gambiae and CYP6N1, CYP6M7, CYP6M1 and CYP6Z1 in An. funestus. Indeed, very similar suites of P450 enzymes commonly associated with resistant populations elsewhere in Africa were detected as over expressed suggesting a convergence of mechanisms across Sub-Saharan African malaria vectors. The findings give insight into factors that may correlate with pyrethroid PBO LLIN success, broadly supporting model predictions, but revision to guidelines previously issued by the World Health Organization is warranted.

A whole transcriptomic approach provides novel insights into the molecular basis of organophosphate and pyrethroid resistance in Anopheles arabiensis from Ethiopia

The development of insecticide resistance in malaria vectors is of increasing concern in Ethiopia because of its potential implications for vector control failure. To better elucidate the specificity of resistance mechanisms and to facilitate the design of control strategies that minimize the likelihood of selecting for cross-resistance, a whole transcriptomic approach was used to explore gene expression patterns in a multi-insecticide resistant population of Anopheles arabiensis from Oromia Region, Ethiopia. This field population was resistant to the diagnostic doses of malathion (average mortality of 71.9%) and permethrin (77.4%), with pools of survivors and unexposed individuals analyzed using Illumina RNA-sequencing, alongside insecticide susceptible reference strains. This population also demonstrated deltamethrin resistance but complete susceptibility to alpha-cypermethrin, bendiocarb and propoxur, providing a phenotypic basis for detecting insecticide-specific resistance mechanisms. Transcriptomic data revealed overexpression of genes including cytochrome P450s, glutathione-s-transferases and carboxylesterases (including CYP4C36, CYP6AA1, CYP6M2, CYP6M3, CYP6P4, CYP9K1, CYP9L1, GSTD3, GSTE2, GSTE3, GSTE4, GSTE5, GSTE7 and two carboxylesterases) that were shared between malathion and permethrin survivors. We also identified nineteen highly overexpressed cuticular-associated proteins (including CYP4G16, CYP4G17 and chitinase) and eighteen salivary gland proteins (including D7r4 short form salivary protein), which may be contributing to a non-specific resistance phenotype by either enhancing the cuticular barrier or promoting binding and sequestration of insecticides, respectively. These findings provide novel insights into the molecular basis of insecticide resistance in this lesser well-characterized major malaria vector species.

Evaluation of the interaction between insecticide resistance-associated genes and malaria transmission in Anopheles gambiae sensu lato in central Côte d'Ivoire

Background

There is evidence that the knockdown resistance gene (Kdr) L1014F and acetylcholinesterase-1 gene (Ace-1^R) G119S mutations involved in pyrethroid and carbamate resistance in Anopheles gambiae influence malaria transmission in sub-Saharan Africa. This is likely due to changes in the behaviour, life history and vector competence and capacity of An. gambiae. In the present study, performed as part of a two-arm cluster randomized controlled trial evaluating the impact of household screening plus a novel insecticide delivery system (In2Care Eave Tubes), we investigated the distribution of insecticide target site mutations and their association with infection status in wild An. gambiae sensu lato (s.l.) populations.

Methods

Mosquitoes were captured in 40 villages around Bouaké by human landing catch from May 2017 to April 2019. Randomly selected samples of An. gambiae s.l. that were infected or not infected with Plasmodium sp. were identified to species and then genotyped for Kdr L1014F and Ace-1^R G119S mutations using quantitative polymerase chain reaction assays. The frequencies of the two alleles were compared between Anopheles coluzzii and Anopheles gambiae and then between infected and uninfected groups for each species.

Results

The presence of An. gambiae (49%) and An. coluzzii (51%) was confirmed in Bouaké. Individuals of both species infected with Plasmodium parasites were found. Over the study period, the average frequency of the Kdr L1014F and Ace-1^R G119S mutations did not vary significantly between study arms. However, the frequencies of the Kdr L1014F and Ace-1^R G119S resistance alleles were significantly higher in An. gambiae than in An. coluzzii [odds ratio (95% confidence interval): 59.64 (30.81–131.63) for Kdr, and 2.79 (2.17–3.60) for Ace-1^R]. For both species, there were no significant differences in Kdr L1014F or Ace-1^R G119S genotypic and allelic frequency distributions between infected and uninfected specimens (P > 0.05).

Conclusions

Either alone or in combination, Kdr L1014F and Ace-1^R G119S showed no significant association with Plasmodium infection in wild An. gambiae and An. coluzzii, demonstrating the similar competence of these species for Plasmodium transmission in Bouaké. Additional factors including behavioural and environmental ones that influence vector competence in natural populations, and those other than allele measurements (metabolic resistance factors) that contribute to resistance, should be considered when establishing the existence of a link between insecticide resistance and vector competence.

Graphical Abstract

The need for practical insecticide-resistance guidelines to effectively inform mosquito-borne disease control programs

Monitoring local mosquito populations for insecticide resistance is critical for effective vector-borne disease control. However, widely used phenotypic assays, which are designed to monitor the emergence and spread of insecticide resistance (technical resistance), do not translate well to the efficacy of vector control products to suppress mosquito numbers in the field (practical resistance). This is because standard testing conditions such as environmental conditions, exposure dose, and type of substrate differ dramatically from those experienced by mosquitoes under field conditions. In addition, field mosquitoes have considerably different physiological characteristics such as age and blood-feeding status. Beyond this, indirect impacts of insecticide resistance and/or exposure on mosquito longevity, pathogen development, host-seeking behavior, and blood-feeding success impact disease transmission. Given the limited number of active ingredients currently available and the observed discordance between resistance and disease transmission, we conclude that additional testing guidelines are needed to determine practical resistance-the efficacy of vector control tools under relevant local conditions- in order to obtain programmatic impact.

CRISPR/Cas9 modified An. gambiae carrying kdr mutation L1014F functionally validate its contribution in insecticide resistance and combined effect with metabolic enzymes

Insecticide resistance in Anopheles mosquitoes is a major obstacle in maintaining the momentum in reducing the malaria burden; mitigating strategies require improved understanding of the underlying mechanisms. Mutations in the target site of insecticides (the voltage gated sodium channel for the most widely used pyrethroid class) and over-expression of detoxification enzymes are commonly reported, but their relative contribution to phenotypic resistance remain poorly understood. Here we present a genome editing pipeline to introduce single nucleotide polymorphisms in An. gambiae which we have used to study the effect of the classical kdr mutation L1014F (L995F based on An. gambiae numbering), one of the most widely distributed resistance alleles. Introduction of 1014F in an otherwise fully susceptible genetic background increased levels of resistance to all tested pyrethroids and DDT ranging from 9.9-fold for permethrin to >24-fold for DDT. The introduction of the 1014F allele was sufficient to reduce mortality of mosquitoes after exposure to deltamethrin treated bednets, even as the only resistance mechanism present. When 1014F was combined with over-expression of glutathione transferase Gste2, resistance to permethrin increased further demonstrating the critical combined effect between target site resistance and detoxification enzymes in vivo. We also show that mosquitoes carrying the 1014F allele in homozygosity showed fitness disadvantages including increased mortality at the larval stage and a reduction in fecundity and adult longevity, which can have consequences for the strength of selection that will apply to this allele in the field.

Escalation of pyrethroid resistance in Anopheles mosquitoes threatens to reduce the effectiveness of our most important tools in malaria control. Studying the mechanisms underlying insecticide resistance is critical to design mitigation strategies. Here, using genome modified mosquitoes, we functionally characterize the most prevalent mutation in resistant mosquitoes, showing that it confers substantial levels of resistance to all tested pyrethroids and undermines the performance of pyrethroid-treated nets. Furthermore, we show that combining this mutation with elevated levels of a detoxification enzyme further increases resistance. The pipeline we have developed provides a robust approach to quantifying the contribution of different combinations of resistance mechanisms to the overall phenotype, providing the missing link between resistance monitoring and predictions of resistance impact.

Influence of insecticide resistance on the biting and resting preferences of malaria vectors in the Gambia

Background

The scale-up of indoor residual spraying and long-lasting insecticidal nets, together with other interventions have considerably reduced the malaria burden in The Gambia. This study examined the biting and resting preferences of the local insecticide-resistant vector populations few years following scale-up of anti-vector interventions.

Method

Indoor and outdoor-resting Anopheles gambiae mosquitoes were collected between July and October 2019 from ten villages in five regions in The Gambia using pyrethrum spray collection (indoor) and prokopack aspirator from pit traps (outdoor). Polymerase chain reaction assays were performed to identify molecular species, insecticide resistance mutations, Plasmodium infection rate and host blood meal.

Results

A total of 844 mosquitoes were collected both indoors (421, 49.9%) and outdoors (423, 50.1%). Four main vector species were identified, including An. arabiensis (indoor: 15%, outdoor: 26%); An. coluzzii (indoor: 19%, outdoor: 6%), An. gambiae s.s. (indoor: 11%, outdoor: 16%), An. melas (indoor: 2%, outdoor: 0.1%) and hybrids of An. coluzzii-An. gambiae s.s (indoors: 3%, outdoors: 2%). A significant preference for outdoor resting was observed in An. arabiensis (Pearson X² = 22.7, df = 4, P<0.001) and for indoor resting in An. coluzzii (Pearson X² = 55.0, df = 4, P<0.001). Prevalence of the voltage-gated sodium channel (Vgsc)-1014S was significantly higher in the indoor-resting (allele freq. = 0.96, 95%CI: 0.78–1, P = 0.03) than outdoor-resting (allele freq. = 0.82, 95%CI: 0.76–0.87) An. arabiensis population. For An. coluzzii, the prevalence of most mutation markers was higher in the outdoor (allele freq. = 0.92, 95%CI: 0.81–0.98) than indoor-resting (allele freq. = 0.78, 95%CI: 0.56–0.86) mosquitoes. However, in An. gambiae s.s., the prevalence of Vgsc-1014F, Vgsc-1575Y and GSTe2-114T was high (allele freq. = 0.96–1), but did not vary by resting location. The overall sporozoite positivity rate was 1.3% (95% CI: 0.5–2%) in mosquito populations. Indoor-resting An. coluzzii had mainly fed on human blood while indoor-resting An. arabiensis fed on animal blood.

Conclusion

In this study, high levels of resistance mutations were observed that could be influencing the mosquito populations to rest indoors or outdoors. The prevalent animal-biting behaviour demonstrated in the mosquito populations suggest that larval source management could be an intervention to complement vector control in this setting.

Unravelling population structure heterogeneity within the genome of the malaria vector Anopheles gambiae

Background

Whole genome re-sequencing provides powerful data for population genomic studies, allowing robust inferences of population structure, gene flow and evolutionary history. For the major malaria vector in Africa, Anopheles gambiae, other genetic aspects such as selection and adaptation are also important. In the present study, we explore population genetic variation from genome-wide sequencing of 765 An. gambiae and An. coluzzii specimens collected from across Africa. We used t-SNE, a recently popularized dimensionality reduction method, to create a 2D-map of An. gambiae and An. coluzzii genes that reflect their population structure similarities.

Results

The map allows intuitive navigation among genes distributed throughout the so-called “mainland” and numerous surrounding “island-like” gene clusters. These gene clusters of various sizes correspond predominantly to low recombination genomic regions such as inversions and centromeres, and also to recent selective sweeps. Because this mosquito species complex has been studied extensively, we were able to support our interpretations with previously published findings. Several novel observations and hypotheses are also made, including selective sweeps and a multi-locus selection event in Guinea-Bissau, a known intense hybridization zone between An. gambiae and An. coluzzii.

Conclusions

Our results present a rich dataset that could be utilized in functional investigations aiming to shed light onto An. gambiae s.l genome evolution and eventual speciation. In addition, the methodology presented here can be used to further characterize other species not so well studied as An. gambiae, shortening the time required to progress from field sampling to the identification of genes and genomic regions under unique evolutionary processes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07722-y.

Mechanisms of amitraz resistance in a Rhipicephalus microplus strain from southern Brazil

Amitraz is one of the most used acaricides for the control of ticks of domestic animals, however, extensive use of this active ingredient has favored the development of resistant populations of Rhipicephalus microplus worldwide. The possible mechanisms of metabolic and/or target-site alterations mechanisms of amitraz resistance were investigated in a Brazilian field population of R. microplus (São Gabriel strain). Bioassays with the synergists piperonylbutoxide, triphenylphosphate and diethyl-maleate were used to evaluate the metabolic mechanisms involved. Target-site insensitivity was investigated by amplification and sequencing of a fragment of the octopamine/tyramine (OCT/TYR) receptor gene. Piperonylbutoxide synergism (synergism ratio = 2.8) indicated the participation of the P450 pathway in the detoxification of amitraz. Previously reported single nucleotide polymorphisms that confer amino acid changes in the OCT/TYR receptor, threonine to proline (T8P) and leucine to serine (L22S), were found in the amitraz-resistant strain but not in the susceptible reference strain. The results suggest that amitraz resistance in the studied strain is multi-factorial and may result from cytochrome P450 detoxification and mutations in octopamine receptors.

Genome-Wide Transcriptional Analysis and Functional Validation Linked a Cluster of Epsilon Glutathione S-Transferases with Insecticide Resistance in the Major Malaria Vector Anopheles funestus across Africa

Resistance is threatening the effectiveness of insecticide-based interventions in use for malaria control. Pinpointing genes associated with resistance is crucial for evidence-based resistance management targeting the major malaria vectors. Here, a combination of RNA-seq based genome-wide transcriptional analysis and RNA-silencing in vivo functional validation were used to identify key insecticide resistance genes associated with DDT and DDT/permethrin cross-resistance across Africa. A cluster of glutathione-S-transferase from epsilon group were found to be overexpressed in resistant populations of Anopheles funestus across Africa including GSTe1 [Cameroon (fold change, FC: 2.54), Ghana (4.20), Malawi (2.51)], GSTe2 [Cameroon (4.47), Ghana (7.52), Malawi (2.13)], GSTe3 [Cameroon (2.49), Uganda (2.60)], GSTe4 in Ghana (3.47), GSTe5 [Ghana (2.94), Malawi (2.26)], GSTe6 [Cameroun (3.0), Ghana (3.11), Malawi (3.07), Uganda (3.78)] and GSTe7 (2.39) in Ghana. Validation of GSTe genes expression profiles by qPCR confirmed that the genes are differentially expressed across Africa with a greater overexpression in DDT-resistant mosquitoes. RNAi-based knock-down analyses supported that five GSTe genes are playing a major role in resistance to pyrethroids (permethrin and deltamethrin) and DDT in An. funestus, with a significant recovery of susceptibility observed when GSTe2, 3, 4, 5 and GSTe6 were silenced. These findings established that GSTe3, 4, 5 and 6 contribute to DDT resistance and should be further characterized to identify their specific genetic variants, to help design DNA-based diagnostic assays, as previously done for the 119F-GSTe2 mutation. This study highlights the role of GSTes in the development of resistance to insecticides in malaria vectors and calls for actions to mitigate this resistance.

Phenotypic Insecticide Resistance in Anopheles gambiae (Diptera: Culicidae): Specific Characterization of Underlying Resistance Mechanisms Still Matters

An effective control of malaria vectors requires an extensive knowledge of mechanisms underlying the resistance-phenotypes developed by these vectors against insecticides. We investigated Anopheles gambiae mosquitoes from Benin and Togo for their intensity of insecticide resistance and we discussed the involvement of genotyped mechanisms in the resistance-phenotypes observed. Three- to five-day-old adult mosquitoes emerged from field and laboratory An. gambiae larvae were assayed using WHO tube intensity tests against various doses of deltamethrin: 1× (0.05%); 2× (0.1%); 5× (0.25%); 7.5× (0.375%) and those of pirimiphos-methyl: 0.5× (0.125%); 1× (0.25%). Members of An. gambiae complex were screened in field populations using polymerase chain reaction (PCR) assays. The presence of kdrR(1014F/1014S) and ace-1R(119S) mutations was also investigated using TaqMan and PCR-RFLP techniques, respectively. Anopheles gambiae from field were very resistant to deltamethrin, whereas KisKdr and AcerKdrKis strains displayed 100% mortality rates at 2× the diagnostic dose. In contrast, the field mosquitoes displayed a low resistance-intensity against 1× the diagnostic dose of pirimiphos-methyl, whereas AcerKis and AcerKdrKis strains showed susceptibility at 0.5× the diagnostic dose. Anopheles gambiae s.s., Anopheles coluzzii, and Anopheles arabiensis were identified. Allelic frequencies of kdrR (1014F) and ace-1R (119S) mutations in the field populations varied from 0.65 to 1 and 0 to 0.84, respectively. The field An. gambiae displayed high-resistance levels against deltamethrin and pirimiphos-methyl when compared with those of the laboratory An. gambiae-resistant strains. These results exhibit the complexity of underlying insecticide resistance mechanisms in these field malaria vectors.

High pyrethroid/DDT resistance in major malaria vector Anopheles coluzzii from Niger-Delta of Nigeria is probably driven by metabolic resistance mechanisms

Entomological surveillance of local malaria vector populations is an important component of vector control and resistance management. In this study, the resistance profile and its possible mechanisms was characterised in a field population of the major malaria vector Anopheles coluzzii from Port Harcourt, the capital of Rivers state, in the Niger-Delta Region of Nigeria. Larvae collected in Port-Harcourt, were reared to adulthood and used for WHO bioassays. The population exhibited high resistance to permethrin, deltamethrin and DDT with mortalities of 6.7% ± 2.4, 37.5% ± 3.2 and 6.3% ± 4.1, respectively, but were fully susceptible to bendiocarb and malathion. Synergist bioassays with piperonylbutoxide (PBO) partially recovered susceptibility, with mortalities increasing to 53% ± 4, indicating probable role of CYP450s in permethrin resistance (χ2 = 29.48, P < 0.0001). Transcriptional profiling revealed five major resistance-associated genes overexpressed in the field samples compared to the fully susceptible laboratory colony, Ngoussou. Highest fold change (FC) was observed with GSTe2 (FC = 3.3 in permethrin exposed and 6.2 in unexposed) and CYP6Z3 (FC = 1.4 in exposed and 4.6 in unexposed). TaqMan genotyping of 32 F0 females detected the 1014F and 1575Y knockdown resistance (kdr) mutations with frequencies of 0.84 and 0.1, respectively, while 1014S mutation was not detected. Sequencing of a fragment of the voltage-gated sodium channel, spanning exon 20 from 13 deltamethrin-resistant and 9 susceptible females revealed only 2 distinct haplotypes with a low haplotype diversity of 0.33. The findings of high pyrethroid resistance but with a significant degree of recovery after PBO synergist assay suggests the need to move to PBO-based nets. This could be complemented with carbamate- or organophosphate-based indoor residual spraying in this area.

Combined over-expression of two cytochrome P450 genes exacerbates the fitness cost of pyrethroid resistance in the major African malaria vector Anopheles funestus

Metabolic resistance driven by multiple P450 genes is worsening insecticide resistance in malaria vectors. However, it remains unclear whether such multiple over-expression imposes an additive fitness cost in the vectors. Here, we showed that two highly over-expressed P450 genes (CYP6P9a and CYP6P9b) combine to impose additive fitness costs in pyrethroid-resistant Anopheles funestus. Genotyping of the CYP6P9b resistance allele in hybrid mosquitoes from a pyrethroid-resistant FUMOZ-R and the susceptible FANG strains revealed that this gene imposes a fitness cost in resistant mosquitoes similar to CYP6P9a. Homozygote susceptible CYP6P9b_S (SS) significantly lay more eggs than the resistant (OR = 2.2, P = 0.04) and with greater hatching rate (p < 0.04). Homozygote resistant larvae CYP6P9b_R (RR) developed significantly slower than homozygote susceptible from L1-L4 (χ² = 7.2; P = 0.007) with a late pupation observed for RR compared to both heterozygotes and homozygotes susceptible (χ² = 11.17; P = 0.0008). No difference was observed between genotypes for adult longevity with no change in allele frequency and gene expression across the lifespan. Furthermore, we established that CYP6P9b combines with CYP6P9a to additively exacerbate the fitness cost of pyrethroid resistance with a greater reduction in fecundity/fertility and increased developmental time of double homozygote resistant mosquitoes. Moreover, an increased proportion of double homozygote susceptible individuals was noted over 10 generations in the insecticide-free environment (χ² = 6.3; P = 0.01) suggesting a reversal to susceptibility in the absence of selection. Such greater fitness cost imposed by multiple P450 genes shows that resistance management strategy based on rotation could help slow the spread of resistance.

Genome variation and population structure among 1142 mosquitoes of the African malaria vector species Anopheles gambiae and Anopheles coluzzii

Mosquito control remains a central pillar of efforts to reduce malaria burden in sub-Saharan Africa. However, insecticide resistance is entrenched in malaria vector populations, and countries with a high malaria burden face a daunting challenge to sustain malaria control with a limited set of surveillance and intervention tools. Here we report on the second phase of a project to build an open resource of high-quality data on genome variation among natural populations of the major African malaria vector species Anopheles gambiae and Anopheles coluzzii We analyzed whole genomes of 1142 individual mosquitoes sampled from the wild in 13 African countries, as well as a further 234 individuals comprising parents and progeny of 11 laboratory crosses. The data resource includes high-confidence single-nucleotide polymorphism (SNP) calls at 57 million variable sites, genome-wide copy number variation (CNV) calls, and haplotypes phased at biallelic SNPs. We use these data to analyze genetic population structure and characterize genetic diversity within and between populations. We illustrate the utility of these data by investigating species differences in isolation by distance, genetic variation within proposed gene drive target sequences, and patterns of resistance to pyrethroid insecticides. This data resource provides a foundation for developing new operational systems for molecular surveillance and for accelerating research and development of new vector control tools. It also provides a unique resource for the study of population genomics and evolutionary biology in eukaryotic species with high levels of genetic diversity under strong anthropogenic evolutionary pressures.

Fine scale spatial investigation of multiple insecticide resistance and underlying target-site and metabolic mechanisms in Anopheles gambiae in central Côte d'Ivoire

Routine monitoring of occurrence, levels and mechanisms of insecticide resistance informs effective management strategies, and should be used to assess the effect of new tools on resistance. As part of a cluster randomised controlled trial evaluating a novel insecticide-based intervention in central Côte d’Ivoire, we assessed resistance and its underlying mechanisms in Anopheles gambiae populations from a subset of trial villages. Resistance to multiple insecticides in An. gambiae s.s. and An. coluzzii was detected across villages, with dose–response assays demonstrating extremely high resistance intensity to the pyrethroid deltamethrin (> 1,500-fold), and mortality following exposure to pyrethroid-treated bednets was low (< 30% mortality in cone bioassays). The 1014F kdr mutation was almost fixed (≥ 90%) in all villages but the 1575Y kdr-amplifying mutation was relatively rare (< 15%). The carbamate and organophosphate resistance-associated Ace-1 G119S mutation was also detected at moderate frequencies (22–43%). Transcriptome analysis identified overexpression of P450 genes known to confer pyrethroid resistance (Cyp9K1, Cyp6P3, and Cyp6M2), and also a carboxylesterase (COEAE1F) as major candidates. Cyp6P3 expression was high but variable (up to 33-fold) and correlated positively with deltamethrin resistance intensity across villages (r² = 0.78, P = 0.02). Tools and strategies to mitigate the extreme and multiple resistance provided by these mechanisms are required in this area to avoid future control failures.

Anopheles metabolic proteins in malaria transmission, prevention and control: a review

The increasing resistance to currently available insecticides in the malaria vector, Anopheles mosquitoes, hampers their use as an effective vector control strategy for the prevention of malaria transmission. Therefore, there is need for new insecticides and/or alternative vector control strategies, the development of which relies on the identification of possible targets in Anopheles. Some known and promising targets for the prevention or control of malaria transmission exist among Anopheles metabolic proteins. This review aims to elucidate the current and potential contribution of Anopheles metabolic proteins to malaria transmission and control. Highlighted are the roles of metabolic proteins as insecticide targets, in blood digestion and immune response as well as their contribution to insecticide resistance and Plasmodium parasite development. Furthermore, strategies by which these metabolic proteins can be utilized for vector control are described. Inhibitors of Anopheles metabolic proteins that are designed based on target specificity can yield insecticides with no significant toxicity to non-target species. These metabolic modulators combined with each other or with synergists, sterilants, and transmission-blocking agents in a single product, can yield potent malaria intervention strategies. These combinations can provide multiple means of controlling the vector. Also, they can help to slow down the development of insecticide resistance. Moreover, some metabolic proteins can be modulated for mosquito population replacement or suppression strategies, which will significantly help to curb malaria transmission.

Insecticide resistance in indoor and outdoor-resting Anopheles gambiae in Northern Ghana

BACKGROUND

Selection pressure from continued exposure to insecticides drives development of insecticide resistance and changes in resting behaviour of malaria vectors. There is need to understand how resistance drives changes in resting behaviour within vector species. The association between insecticide resistance and resting behaviour of Anopheles gambiae sensu lato (s.l.) in Northern Ghana was examined.

METHODS

F1 progenies from adult mosquitoes collected indoors and outdoors were exposed to DDT, deltamethrin, malathion and bendiocarb using WHO insecticide susceptibility tests. Insecticide resistance markers including voltage-gated sodium channel (Vgsc)-1014F, Vgsc-1014S, Vgsc-1575Y, glutathione-S-transferase epsilon 2 (GSTe2)-114T and acetylcholinesterase (Ace1)-119S, as well as blood meal sources were investigated using PCR methods. Activities of metabolic enzymes, acetylcholine esterase (AChE), non-specific β-esterases, glutathione-S-transferase (GST) and monooxygenases were measured from unexposed F1 progenies using microplate assays.

RESULTS

Susceptibility of Anopheles coluzzii to deltamethrin 24 h post-exposure was significantly higher in indoor (mortality = 5%) than outdoor (mortality = 2.5%) populations (P = 0.02). Mosquitoes were fully susceptible to malathion (mortality: indoor = 98%, outdoor = 100%). Susceptibility to DDT was significantly higher in outdoor (mortality = 9%) than indoor (mortality = 0%) mosquitoes (P = 0.006). Mosquitoes were also found with suspected resistance to bendiocarb but mortality was not statistically different (mortality: indoor = 90%, outdoor = 95%. P = 0.30). Frequencies of all resistance alleles were higher in F1 outdoor (0.11-0.85) than indoor (0.04-0.65) mosquito populations, while Vgsc-1014F in F0 An. gambiae sensu stricto (s.s) was significantly associated with outdoor-resting behaviour (P = 0.01). Activities of non-specific β-esterase enzymes were significantly higher in outdoor than indoor mosquitoes (Mean enzyme activity: Outdoor = : 1.70/mg protein; Indoor = 1.35/mg protein. P < 0.0001). AChE activity was also more elevated in outdoor (0.62/mg protein) than indoor (0.57/mg protein) mosquitoes but this was not significant (P = 0.08). Human blood index (HBI) was predominantly detected in indoor (18%) than outdoor mosquito populations (3%).

CONCLUSIONS

The overall results did not establish that there was a significant preference of resistant malaria vectors to solely rest indoors or outdoors, but varied depending on the resistant alleles present. Phenotypic resistance was higher in indoor than outdoor-resting mosquitoes, but genotypic and metabolic resistance levels were higher in outdoor than the indoor populations. Continued monitoring of changes in resting behaviour within An. gambiae s.l. populations is recommended.

Genomic Variant Analyses in Pyrethroid Resistant and Susceptible Malaria Vector, Anopheles sinensis

Anopheles sinensis is a major malaria vector in Southeast Asia. Resistance to pyrethroid insecticides in this species has impeded malaria control in the region. Previous studies found that An. sinensis populations from Yunnan Province, China were highly resistant to deltamethrin and did not carry mutations in the voltage-gated sodium channel gene that cause knockdown resistance. In this study, we tested the hypothesis that other genomic variants are associated with the resistance phenotype. Using paired-end whole genome sequencing (DNA-seq), we generated 108 Gb of DNA sequence from deltamethrin -resistant and -susceptible mosquito pools with an average coverage of 83.3× depth. Using a stringent filtering method, we identified a total of 916,926 single nucleotide variants (SNVs), including 32,240 non-synonymous mutations. A total of 958 SNVs differed significantly in allele frequency between deltamethrin -resistant and -susceptible mosquitoes. Of these, 43 SNVs were present within 37 genes that code for immunity, detoxification, cuticular, and odorant proteins. A subset of 12 SNVs were randomly selected for genotyping of individual mosquitoes by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) and showed consistent allele frequencies with the pooled DNA-seq derived allele frequencies. In addition, copy number variations (CNVs) were detected in 56 genes, including 33 that contained amplification alleles and 23 that contained deletion alleles in resistant mosquitoes compared to susceptible mosquitoes. The genomic variants described here provide a useful resource for future studies on the genetic mechanism of insecticide resistance in this important malaria vector species.

Polymorphisms in GSTE2 is associated with temephos resistance in Aedes aegypti

The glutathione S-transferases (GSTs) are enzymes involved in several distinct biological processes. In insects, the GSTs, especially delta and epsilon classes, play a key role in the metabolism of xenobiotics used to control insect populations. Here, we investigated its potential role in temephos resistance, examining the GSTE2 gene from susceptible (RecL) and resistant (RecR) strains of the mosquito Aedes aegypti, vector for several pathogenic arboviruses. Total GST enzymatic activity and the GSTE2 gene expression profile were evaluated, with the GSTE2 cDNA and genomic loci sequenced from both strains. Recombinant GSTE2 and mutants were produced in a heterologous expression system and assayed for enzyme kinetic parameters. These proteins also had their 3D structure predicted through molecular modeling. Our results showed that RecR has a profile of total GST enzymatic activity higher than RecL, with the expression of the GSTE2 gene in resistant larvae increasing six folds. Four exclusive RecR mutations were observed (L111S, I150V, E178A and A198E), which were absent in the laboratory susceptible strains. The enzymatic activity of the recombinant GSTE2 showed different kinetic parameters, with the GSTE2 RecR showing an enhanced ability to metabolize its substrate. The I150V mutation was shown to induce significant changes in catalytic parameters and a 3D modeling of GSTE2 mapped two of the RecR changes (L111S and I150V) near the enzyme's catalytic pocket, also implying an impact on its catalytic activity. Our results reinforce a potential role for GSTE2 in the metabolic resistance phenotype while contributing to the understanding of the molecular basis for the resistance mechanism.

Assessing connectivity despite high diversity in island populations of a malaria mosquito

Documenting isolation is notoriously difficult for species with vast polymorphic populations. High proportions of shared variation impede estimation of connectivity, even despite leveraging information from many genetic markers. We overcome these impediments by combining classical analysis of neutral variation with assays of the structure of selected variation, demonstrated using populations of the principal African malaria vector Anopheles gambiae. Accurate estimation of mosquito migration is crucial for efforts to combat malaria. Modeling and cage experiments suggest that mosquito gene drive systems will enable malaria eradication, but establishing safety and efficacy requires identification of isolated populations in which to conduct field testing. We assess Lake Victoria islands as candidate sites, finding one island 30 km offshore is as differentiated from mainland samples as populations from across the continent. Collectively, our results suggest sufficient contemporary isolation of these islands to warrant consideration as field-testing locations and illustrate shared adaptive variation as a useful proxy for connectivity in highly polymorphic species.

An Experimental Hut Evaluation of PBO-Based and Pyrethroid-Only Nets against the Malaria Vector Anopheles funestus Reveals a Loss of Bed Nets Efficacy Associated with GSTe2 Metabolic Resistance

Growing insecticide resistance in malaria vectors is threatening the effectiveness of insecticide-based interventions, including Long Lasting Insecticidal Nets (LLINs). However, the impact of metabolic resistance on the effectiveness of these tools remains poorly characterized. Using experimental hut trials and genotyping of a glutathione S-transferase resistance marker (L119F-GSTe2), we established that GST-mediated resistance is reducing the efficacy of LLINs against Anopheles funestus. Hut trials performed in Cameroon revealed that Piperonyl butoxide (PBO)-based nets induced a significantly higher mortality against pyrethroid resistant An. funestus than pyrethroid-only nets. Blood feeding rate and deterrence were significantly higher in all LLINs than control. Genotyping the L119F-GSTe2 mutation revealed that, for permethrin-based nets, 119F-GSTe2 resistant mosquitoes have a greater ability to blood feed than susceptible while the opposite effect is observed for deltamethrin-based nets. For Olyset Plus, a significant association with exophily was observed in resistant mosquitoes (OR = 11.7; p < 0.01). Furthermore, GSTe2-resistant mosquitoes (cone assays) significantly survived with PermaNet 2.0 (OR = 2.1; p < 0.01) and PermaNet 3.0 (side) (OR = 30.1; p < 0.001) but not for Olyset Plus. This study shows that the efficacy of PBO-based nets (e.g., blood feeding inhibition) against pyrethroid resistant malaria vectors could be impacted by other mechanisms including GST-mediated metabolic resistance not affected by the synergistic action of PBO. Mosaic LLINs incorporating a GST inhibitor (diethyl maleate) could help improve their efficacy in areas of GST-mediated resistance.

Functional genetic validation of key genes conferring insecticide resistance in the major African malaria vector, Anopheles gambiae

Insecticide resistance in Anopheles gambiae mosquitoes can derail malaria control programs, and to overcome it, we need to discover the underlying molecular basis. Here, we characterize 3 genes most often associated with insecticide resistance directly by their overproduction in genetically modified An. gambiae. We show that overexpression of each gene confers resistance to representatives of at least 1 insecticide class, and taken together, the 3 genes provide cross-resistance to all 4 major insecticide classes currently used in public health. These data validate the candidate genes as markers to monitor the spread of resistance in mosquito populations. The modified mosquitoes produced are also valuable tools to prescreen the efficacy of new insecticides against existing resistance mechanisms.

Resistance in Anopheles gambiae to members of all 4 major classes (pyrethroids, carbamates, organochlorines, and organophosphates) of public health insecticides limits effective control of malaria transmission in Africa. Increase in expression of detoxifying enzymes has been associated with insecticide resistance, but their direct functional validation in An. gambiae is still lacking. Here, we perform transgenic analysis using the GAL4/UAS system to examine insecticide resistance phenotypes conferred by increased expression of the 3 genes—Cyp6m2, Cyp6p3, and Gste2—most often found up-regulated in resistant An. gambiae. We report evidence in An. gambiae that organophosphate and organochlorine resistance is conferred by overexpression of GSTE2 in a broad tissue profile. Pyrethroid and carbamate resistance is bestowed by similar Cyp6p3 overexpression, and Cyp6m2 confers only pyrethroid resistance when overexpressed in the same tissues. Conversely, such Cyp6m2 overexpression increases susceptibility to the organophosphate malathion, presumably due to conversion to the more toxic metabolite, malaoxon. No resistant phenotypes are conferred when either Cyp6 gene overexpression is restricted to the midgut or oenocytes, indicating that neither tissue is involved in insecticide resistance mediated by the candidate P450s examined. Validation of genes conferring resistance provides markers to guide control strategies, and the observed negative cross-resistance due to Cyp6m2 gives credence to proposed dual-insecticide strategies to overcome pyrethroid resistance. These transgenic An. gambiae-resistant lines are being used to test the “resistance-breaking” efficacy of active compounds early in their development.

A high throughput multi-locus insecticide resistance marker panel for tracking resistance emergence and spread in Anopheles gambiae

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