Rapid high throughput SYBR green assay for identifying the malaria vectors Anopheles arabiensis, Anopheles coluzzii and Anopheles gambiae s.s. Giles

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Pyrethroid resistance and gene expression profile of a new resistant An. gambiae colony from Uganda reveals multiple resistance mechanisms and overexpression of Glutathione-S-Transferases linked to survival of PBO-pyrethroid combination

Background

The effectiveness of long-lasting insecticidal nets (LLINs) are being threatened by growing resistance to pyrethroids. To restore their efficacy, a synergist, piperonyl butoxide (PBO) which inhibits cytochrome P450s has been incorporated into pyrethroid treated nets. A trial of PBO-LLINs was conducted in Uganda from 2017 and we attempted to characterize mechanisms of resistance that could impact intervention efficacy.

Methods

We established an Anopheles gambiae s.s colony in 2018 using female mosquitoes collected from Busia district in eastern Uganda. We first assessed the phenotypic resistance profile of this colony using WHO tube and net assays using a deltamethrin dose-response approach. The Busia colony was screened for known resistance markers and RT-qPCR targeting 15 genes previously associated with insecticide resistance was performed.

Results

The Busia colony had very high resistance to deltamethrin, permethrin and DDT. In addition, the colony had moderate resistance to alpha-cypermethrin and lambda-cyhalothrin but were fully susceptible to bendiocarb and fenitrothion. Exposure to PBO in combination with permethrin and deltamethrin resulted in higher mortality rates in both net and tube assays, with a higher mortality observed in net assays than tube assays. The kdr marker, Vgsc-995S was at very high frequency (91.7-98.9%) whilst the metabolic markers Coeae1d and Cyp4j5-L43F were at very low (1.3% - 11.5%) and moderate (39.5% - 44.7%) frequencies respectively. Our analysis showed that gene expression pattern in mosquitoes exposed to deltamethrin, permethrin or DDT only were similar in comparison to the susceptible strain and there was significant overexpression of cytochrome P450s, glutathione-s-transferases (GSTs) and carboxyl esterases (COEs). However, mosquitoes exposed to both PBO and pyrethroid strikingly and significantly only overexpressed closely related GSTs compared to unexposed mosquitoes while major cytochrome P450s were underexpressed.

Conclusions

The high levels of pyrethroid resistance observed in Busia appears associated with a wide range of metabolic gene families.

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

To keep ahead of the evolution of resistance to insecticides in mosquitoes, national malaria control programmes must make use of a range of insecticides, both old and new, while monitoring resistance mechanisms. Knowledge of the mechanisms of resistance remains limited in Anopheles arabiensis, which in many parts of Africa is of increasing importance because it is apparently less susceptible to many indoor control interventions. Furthermore, comparatively little is known in general about resistance to non-pyrethroid insecticides such as pirimiphos-methyl (PM), which are crucial for effective control in the context of resistance to pyrethroids. We performed a genome-wide association study to determine the molecular mechanisms of resistance to deltamethrin (commonly used in bednets) and PM, in An. arabiensis from two regions in Tanzania. Genomic regions of positive selection in these populations were largely driven by copy number variants (CNVs) in gene families involved in resistance to these two insecticides. We found evidence of a new gene cluster involved in resistance to PM, identifying a strong selective sweep tied to a CNV in the Coeae2g-Coeae6g cluster of carboxylesterase genes. Using complementary data from An. coluzzii in Ghana, we show that copy number at this locus is significantly associated with PM resistance. Similarly, for deltamethrin, resistance was strongly associated with a novel CNV allele in the Cyp6aa / Cyp6p cluster. Against this background of metabolic resistance, target site resistance was very rare or absent for both insecticides. Mutations in the pyrethroid target site Vgsc were at very low frequency in Tanzania, yet combining these samples with three An. arabiensis individuals from West Africa revealed a startling diversity of evolutionary origins of target site resistance, with up to 5 independent origins of Vgsc-995 mutations found within just 8 haplotypes. Thus, despite having been first recorded over 10 years ago, Vgsc resistance mutations in Tanzanian An. arabiensis have remained at stable low frequencies. Overall, our results provide a new copy number marker for monitoring resistance to PM in malaria mosquitoes, and reveal the complex picture of resistance patterns in An. arabiensis.

Efficacy of chlorfenapyr-pyrethroid and piperonyl butoxide-pyrethroid long-lasting insecticidal nets (LLINs) compared to pyrethroid-only LLINs for malaria control in Côte d'Ivoire: a three group, cluster randomised trial

BACKGROUND

The massive scale-up of long-lasting insecticidal nets (LLIN) has led to a major reduction in malaria burden in many sub-Saharan African (SSA) countries. The World Health Organization (WHO) has recently issued a strong recommendation for the use of chlorfenapyr-pyrethroid LLINs compared to standard pyrethroid-only LLINs in areas of high insecticide resistance intensity. However, there is still a lack of conclusive evidence on the efficacy of piperonyl butoxide-pyrethroid (PBO-py) LLINs, especially in West Africa, where vector composition and resistance mechanisms may be different from vectors in East Africa.

METHODS

This is a three-arm, superiority, triple-blinded, cluster randomised trial, with village as the unit of randomisation. This study conducted in Côte d'Ivoire will evaluate the efficacy on epidemiological and entomological outcomes of (1) the control arm: MAGNet® LN, which contains the pyrethroid, alpha-cypermethrin, (2) VEERALIN® LN, a net combining the synergist PBO and alpha-cypermethrin, and (3) Interceptor® G2 LN, which incorporates chlorfenapyr and alpha-cypermethrin, two adulticides with different mechanisms of action. A total of 33 villages with an average of 200 households per village will be identified, mapped, and randomised in a ratio of 1:1:1. Nets will be distributed at a central point following national guidelines with 1 net for every 2 people. The primary outcome of the trial will be incidence of malaria cases (confirmed by rapid diagnostic test (RDT)) in a cohort of 50 children aged 6 months to 10 years in each cluster, followed for 12 months (active case detection). Secondary outcomes are cross-sectional community prevalence of malaria infection (confirmed by RDT) in the study population at 6 and 12 months post-intervention (50 randomly selected persons per cluster), vector density, entomological inoculation rate (EIR), and phenotypic and genotypic insecticide resistance at baseline and 12 months post-intervention in 3 sentinel villages in each treatment arm.

DISCUSSION

In addition to generating further evidence for next-generation LLINs, this study will also provide the first evidence for pyrethroid-PBO nets in a West African setting. This could further inform WHO recommendations on the pragmatic use of pyrethroid-PBO nets.

TRIAL REGISTRATION

ClinicalTrials.gov NCT05796193. Registered on April 3, 2023.

Susceptibility of Anopheles gambiae to Natural Plasmodium falciparum Infection: A Comparison between the Well-Established Anopheles gambiae s.s Line and a Newly Established Ugandan Anopheles gambiae s.s. Line

Much of our understanding of malaria transmission comes from mosquito feeding assays using Anopheles mosquitoes from colonies that are well adapted to membrane feeding. This raises the question whether results from colony mosquitoes lead to overestimates of outcomes in wild Anopheles mosquitoes. We successfully established an Anopheles colony using progeny of wild Anopheles gambiae s.s. mosquitoes (Busia mosquitoes) and directly compared their susceptibility to infection with Plasmodium falciparum with the widely used An. gambiae s.s. mosquitoes (Kisumu mosquitoes) using gametocyte-infected Ugandan donor blood. The proportion of infectious feeds did not differ between Busia (71.8%, 23/32) and Kisumu (68.8%, 22/32, P = 1.00) mosquitoes. When correcting for random effects of donor blood, we observed a 23% higher proportion of infected Busia mosquitoes than infected Kisumu mosquitoes (RR, 1.23; 95% CI, 1.10-1.38, P < 0.001). This study suggests that feeding assays with Kisumu mosquitoes do not overestimate outcomes in wild An. gambiae s.s. mosquitoes, the mosquito species most relevant to malaria transmission in Uganda.

Evaluation of diffuse reflectance spectroscopy for predicting age, species, and cuticular resistance of Anopheles gambiae s.l under laboratory conditions

Mid-infrared spectroscopy (MIRS) combined with machine learning analysis has shown potential for quick and efficient identification of mosquito species and age groups. However, current technology to collect spectra is destructive to the sample and does not allow targeting specific tissues of the mosquito, limiting the identification of other important biological traits such as insecticide resistance. Here, we assessed the use of a non-destructive approach of MIRS for vector surveillance, micro diffuse reflectance spectroscopy (µDRIFT) using mosquito legs to identify species, age and cuticular insecticide resistance within the Anopheles gambiae s.l. complex. These mosquitoes are the major vectors of malaria in Africa and the focus on surveillance in malaria control programs. Legs required significantly less scanning time and showed more spectral consistence compared to other mosquito tissues. Machine learning models were able to identify An. gambiae and An. coluzzii with an accuracy of 0.73, two ages groups (3 and 10 days old) with 0.77 accuracy and we obtained accuracy of 0.75 when identifying cuticular insecticide resistance. Our results highlight the potential of different mosquito tissues and µDRIFT as tools for biological trait identification on mosquitoes that transmit malaria. These results can guide new ways of identifying mosquito traits which can help the creation of innovative surveillance programs by adapting new technology into mosquito surveillance and control tools.

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.

RNA-Seq-Pop: Exploiting the sequence in RNA sequencing-A Snakemake workflow reveals patterns of insecticide resistance in the malaria vector Anopheles gambiae

We provide a reproducible and scalable Snakemake workflow, called RNA-Seq-Pop, which provides end-to-end analysis of RNA sequencing data sets. The workflow allows the user to perform quality control, perform differential expression analyses and call genomic variants. Additional options include the calculation of allele frequencies of variants of interest, summaries of genetic variation and population structure, and genome-wide selection scans, together with clear visualizations. RNA-Seq-Pop is applicable to any organism, and we demonstrate the utility of the workflow by investigating pyrethroid resistance in selected strains of the major malaria mosquito, Anopheles gambiae. The workflow provides additional modules specifically for An. gambiae, including estimating recent ancestry and determining the karyotype of common chromosomal inversions. The Busia laboratory colony used for selections was collected in Busia, Uganda, in November 2018. We performed a comparative analysis of three groups: a parental G24 Busia strain; its deltamethrin-selected G28 offspring; and the susceptible reference strain Kisumu. Measures of genetic diversity reveal patterns consistent with that of laboratory colonization and selection, with the parental Busia strain exhibiting the highest nucleotide diversity, followed by the selected Busia offspring, and finally, Kisumu. Differential expression and variant analyses reveal that the selected Busia colony exhibits a number of distinct mechanisms of pyrethroid resistance, including the Vgsc-995S target-site mutation, upregulation of SAP genes, P450s and a cluster of carboxylesterases. During deltamethrin selections, the 2La chromosomal inversion rose in frequency (from 33% to 86%), supporting a previous link with pyrethroid resistance. RNA-Seq-Pop is hosted at: github.com/sanjaynagi/rna-seq-pop. We anticipate that the workflow will provide a useful tool to facilitate reproducible, transcriptomic studies in An. gambiae and other taxa.

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.

Effects of Elevated Temperatures on the Growth and Development of Adult Anopheles gambiae (s.l.) (Diptera: Culicidae) Mosquitoes

Higher temperatures expected in a future warmer climate could adversely affect the growth and development of mosquitoes. This study investigated the effects of elevated temperatures on longevity, gonotrophic cycle length, biting rate, fecundity, and body size of Anopheles gambiae (s.l.) (Diptera: Culicidae) mosquitoes. Anopheles gambiae (s.l.) eggs obtained from laboratory established colonies were reared under eight temperature regimes (25, 28, 30, 32, 34, 36, 38, and 40°C), and 80 ± 10% RH. All adults were allowed to feed on a 10% sugar solution soaked in cotton wool; however, some mosquitoes were provided blood meal using guinea pig. Longevity was estimated for both blood-fed and non-blood-fed mosquitoes and analyzed using the Kaplan-Meier survival analysis. One-way ANOVA was used to test the effect of temperature on gonotrophic cycle length, biting rate, and fecundity. Adult measurement data were log-transformed and analyzed using ordinary least square regression with robust standard errors. Increasing temperature significantly decreased the longevity of both blood-fed (Log-rank test; X2(4) = 904.15, P < 0.001) and non-blood-fed (Log-rank test; X2(4) = 1163.60, P < 0.001) mosquitoes. In addition, the fecundity of mosquitoes decreased significantly (ANOVA; F(2,57) = 3.46, P = 0.038) with an increase in temperature. Body size (β = 0.14, 95% CI, 0.16, 0.12, P < 0.001) and proboscis length (β = 0.13, 95% CI, 0.17, 0.09, P < 0.001) significantly decreased with increasing temperature from 25 to 34°C. Increased temperatures expected in a future warmer climate could cause some unexpected effects on mosquitoes by directly influencing population dynamics and malaria transmission.

Relationship between temperature and Anopheles gambiae sensu lato mosquitoes' susceptibility to pyrethroids and expression of metabolic enzymes

Background

Malaria remains one of the most devastating diseases globally, and the control of mosquitoes as the vector is mainly dependent on chemical insecticides. Elevated temperatures associated with future warmer climates could affect mosquitoes' metabolic enzyme expression and increase insecticide resistance, making vector control difficult. Understanding how mosquito rearing temperatures influence their susceptibility to insecticide and expression of metabolic enzymes could aid in the development of novel tools and strategies to control mosquitoes in a future warmer climate. This study evaluated the effects of temperature on the susceptibility of Anopheles gambiae sensu lato (s.l.) mosquitoes to pyrethroids and their expression of metabolic enzymes.

Methods

Anopheles gambiae s.l. eggs obtained from laboratory-established colonies were reared under eight temperature regimes (25, 28, 30, 32, 34, 36, 38, and 40 °C). Upon adult emergence, 3- to 5-day-old female non-blood-fed mosquitoes were used for susceptibility tests following the World Health Organization (WHO) bioassay protocol. Batches of 20–25 mosquitoes from each temperature regime (25–34 °C) were exposed to two pyrethroid insecticides (0.75% permethrin and 0.05% deltamethrin). In addition, the levels of four metabolic enzymes (α-esterase, β-esterase, glutathione S-transferase [GST], and mixed-function oxidase [MFO]) were examined in mosquitoes that were not exposed and those that were exposed to pyrethroids.

Results

Mortality in An. gambiae s.l. mosquitoes exposed to deltamethrin and permethrin decreased at temperatures above 28 °C. In addition, mosquitoes reared at higher temperatures were more resistant and had more elevated enzyme levels than those raised at low temperatures. Overall, mosquitoes that survived after being exposed to pyrethroids had higher levels of metabolic enzymes than those that were not exposed to pyrethroids.

Conclusions

This study provides evidence that elevated temperatures decreased An. gambiae s.l. mosquitoes' susceptibility to pyrethroids and increased the expression of metabolic enzymes. This evidence suggests that elevated temperatures projected in a future warmer climate could increase mosquitoes' resistance to insecticides and complicate malaria vector control measures. This study therefore provides vital information, and suggests useful areas of future research, on the effects of temperature variability on mosquitoes that could guide vector control measures in a future warmer climate.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05273-z.

Effects of elevated temperatures on the development of immature stages of Anopheles gambiae (s.l.) mosquitoes

OBJECTIVE

This study investigated the effects of temperature on the development of the immature stages of An. gambiae (s.l.) mosquitoes.

METHODS

Mosquito eggs were obtained from laboratory established colonies and reared under eight temperature regimes (25, 28, 30, 32, 34, 36, 38 and 40 °C), and 80 ± 10% relative humidity. Larvae were checked daily for development to the next stage and for mortality. Pupation success, number of adults produced, and sex ratio of the newly emerged adults were recorded. Larval survival was monitored every 24 hours, and data were analyzed using Kaplan Meier survival analysis. Analysis of variance was used where data followed normal distribution, and a Kruskal-Wallis test where data were not normally distributed. Larval and pupal measurements were log-transformed and analyzed using ordinary least square regression with robust standard errors.

RESULTS

Increasing the temperature from 25 to 36 °C decreased the development time by 10.57 days. Larval survival (X² (6) = 5353.12, P < 0.001) and the number of adults produced (X² (5) = 28.16, P < 0.001) decreased with increasing temperature. Increasing temperatures also resulted in significantly smaller larvae and pupae (P < 0.001). At higher temperatures disproportionately more male than female mosquitoes were produced.

CONCLUSIONS

Increased temperature affected different developmental stages in the life cycle of An. gambiae (s.l.) mosquitoes, from larval to adult emergence. This study contributes to the knowledge on the relationship between temperature and Anopheles mosquitoes and provides useful information for modelling vector population dynamics in the light of climate change.

Resistance to pirimiphos-methyl in West African Anopheles is spreading via duplication and introgression of the Ace1 locus

Vector population control using insecticides is a key element of current strategies to prevent malaria transmission in Africa. The introduction of effective insecticides, such as the organophosphate pirimiphos-methyl, is essential to overcome the recurrent emergence of resistance driven by the highly diverse Anopheles genomes. Here, we use a population genomic approach to investigate the basis of pirimiphos-methyl resistance in the major malaria vectors Anopheles gambiae and A. coluzzii. A combination of copy number variation and a single non-synonymous substitution in the acetylcholinesterase gene, Ace1, provides the key resistance diagnostic in an A. coluzzii population from Côte d'Ivoire that we used for sequence-based association mapping, with replication in other West African populations. The Ace1 substitution and duplications occur on a unique resistance haplotype that evolved in A. gambiae and introgressed into A. coluzzii, and is now common in West Africa primarily due to selection imposed by other organophosphate or carbamate insecticides. Our findings highlight the predictive value of this complex resistance haplotype for phenotypic resistance and clarify its evolutionary history, providing tools to for molecular surveillance of the current and future effectiveness of pirimiphos-methyl based interventions.

High-Resolution Melting (HRM) Curve Assay for the Identification of Eight Fusarium Species Causing Ear Rot in Maize

Maize plants are often infected with fungal pathogens of the genus Fusarium. Taxonomic characterization of these species by microscopic examination of pure cultures or assignment to mating populations is time-consuming and requires specific expertise. Reliable taxonomic assignment may be strengthened by the analysis of DNA sequences. Species-specific PCR assays are available for most Fusarium pathogens, but the number of species that infect maize increases the labor and costs required for analysis. In this work, a diagnostic assay for major Fusarium pathogens of maize based on the analysis of melting curves of PCR amplicons was established. Short segments of genes RPB2 and TEF-1α, which have been widely used in molecular taxonomy of Fusarium, were amplified with universal primers in a real-time thermocycler and high-resolution melting (HRM) curves of the products were recorded. Among major Fusarium pathogens of maize ears, F. cerealis, F. culmorum, F. graminearum, F. equiseti, F. poae, F. temperatum, F. tricinctum, and F. verticillioides, all species except for the pair F. culmorum/F. graminearum could be distinguished by HRM analysis of a 304 bp segment of the RPB2 gene. The latter two species could be differentiated by HRM analysis of a 247 bp segment of the TEF-1α gene. The assay was validated with DNA extracted from pure cultures of fungal strains, successfully applied to total DNA extracted from infected maize ears and also to fungal mycelium that was added directly to the PCR master mix ("colony PCR"). HRM analysis thus offers a cost-efficient method suitable for the diagnosis of multiple fungal pathogens.

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.