Field, genetic, and modeling approaches show strong positive selection acting upon an insecticide resistance mutation in Anopheles gambiae s.s

Vkorc1 polymorphisms of the Norway rats in China: Implications for rodent management and evolutionary origin of anticoagulant resistance mutations

The extensive utilization of second-generation anticoagulant rodenticides (SGARs) has raised concerns regarding non-target animal safety and environmental contamination. It is essential to assess the anticoagulant resistance level in rodent populations and prioritize the use of relative low toxic first-generation anticoagulant rodenticides (FGARs) in susceptible rodent populations. Mutations in the vitamin K epoxide reductase complex subunit 1 (Vkorc1) gene confer anticoagulant resistance in Norway rats. However, the Vkorc1 polymorphisms remain unclear in most Norway rat populations in China although anticoagulant rodenticides have been widely used in China since the 1980s. Analysis of the Vkorc1 polymorphisms in 489 rats across China, combined with in silico binding affinity analysis, revealed three potential resistance mutations A26T, C96Y, and A140T at three distinct locations. In the remaining locations, Vkorc1 resistance mutations were absent, indicating that the FGARs could be effective in these areas. Additional evolutionary analysis of different Vkorc1 mutations suggested that the three missense mutations identified in China might have evolved independently as de novo mutations, and the resistance mutations in Europe are unlikely to be pre-existing variations in China. Further analysis of Vkorc1 haplotypes among European resistant rat populations is essential for understanding the origin of these resistance mutations. These findings emphasize the importance of customizing rodent control strategies in China based on regional resistance levels and gaining insights into the origins of Vkorc1 mutations for more effective rodent management strategies.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Unraveling the secrets: Evolution of resistance mediated by membrane proteins

Membrane protein-mediated resistance is a multidisciplinary challenge that spans fields such as medicine, agriculture, and environmental science. Understanding its complexity and devising innovative strategies are crucial for treating diseases like cancer and managing resistant pests in agriculture. This paper explores the dual nature of resistance mechanisms across different organisms: On one hand, animals, bacteria, fungi, plants, and insects exhibit convergent evolution, leading to the development of similar resistance mechanisms. On the other hand, influenced by diverse environmental pressures and structural differences among organisms, they also demonstrate divergent resistance characteristics. Membrane protein-mediated resistance mechanisms are prevalent across animals, bacteria, fungi, plants, and insects, reflecting their shared survival strategies evolved through convergent evolution to address similar survival challenges. However, variations in ecological environments and biological characteristics result in differing responses to resistance. Therefore, examining these differences not only enhances our understanding of adaptive resistance mechanisms but also provides crucial theoretical support and insights for addressing drug resistance and advancing pharmaceutical development.

Signatures of selective sweeps in continuous-space populations

Selective sweeps describe the process by which an adaptive mutation arises and rapidly fixes in the population, thereby removing genetic variation in its genomic vicinity. The expected signatures of selective sweeps are relatively well understood in panmictic population models, yet natural populations often extend across larger geographic ranges where individuals are more likely to mate with those born nearby. To investigate how such spatial population structure can affect sweep dynamics and signatures, we simulated selective sweeps in populations inhabiting a two-dimensional continuous landscape. The maximum dispersal distance of offspring from their parents can be varied in our simulations from an essentially panmictic population to scenarios with increasingly limited dispersal. We find that in low-dispersal populations, adaptive mutations spread more slowly than in panmictic ones, while recombination becomes less effective at breaking up genetic linkage around the sweep locus. Together, these factors result in a trough of reduced genetic diversity around the sweep locus that looks very similar across dispersal rates. We also find that the site frequency spectrum around hard sweeps in low-dispersal populations becomes enriched for intermediate-frequency variants, making these sweeps appear softer than they are. Furthermore, haplotype heterozygosity at the sweep locus tends to be elevated in low-dispersal scenarios as compared to panmixia, contrary to what we observe in neutral scenarios without sweeps. The haplotype patterns generated by these hard sweeps in low-dispersal populations can resemble soft sweeps from standing genetic variation that arose from substantially older alleles. Our results highlight the need for better accounting for spatial population structure when making inferences about selective sweeps.

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.

LLIN Evaluation in Uganda Project (LLINEUP)-effects of a vector control trial on Plasmodium infection prevalence and genotypic markers of insecticide resistance in Anopheles vectors from 48 districts of Uganda

Pyrethroid bednets treated with the synergist piperonyl butoxide (PBO) offer the possibility of improved vector control in mosquito populations with metabolic resistance. In 2017-2019, we conducted a large-scale, cluster-randomised trial (LLINEUP) to evaluate long-lasting insecticidal nets (LLINs) treated with a pyrethroid insecticide plus PBO (PBO LLINs), as compared to conventional, pyrethroid-only LLINs across 104 health sub-districts (HSDs) in Uganda. In LLINEUP, and similar trials in Tanzania, PBO LLINs were found to provide greater protection against malaria than conventional LLINs, reducing parasitaemia and vector density. In the LLINEUP trial, we conducted cross-sectional household entomological surveys at baseline and then every 6 months for two years, which we use here to investigate longitudinal changes in mosquito infection rate and genetic markers of resistance. Overall, 5395 female Anopheles mosquitoes were collected from 5046 households. The proportion of mosquitoes infected (PCR-positive) with Plasmodium falciparum did not change significantly over time, while infection with non-falciparum malaria decreased in An. gambiae s.s., but not An. funestus. The frequency of genetic markers associated with pyrethroid resistance increased significantly over time, but the rate of change was not different between the two LLIN types. The knock-down resistance (kdr) mutation Vgsc-995S declined over time as Vgsc-995F, the alternative resistance mutation at this codon, increased. Vgsc-995F appears to be spreading into Uganda. Distribution of LLINs in Uganda was previously found to be associated with reductions in parasite prevalence and vector density, but here we show that the proportion of infective mosquitoes remained stable across both PBO and non-PBO LLINs, suggesting that the potential for transmission persisted. The increased frequency of markers of pyrethroid resistance indicates that LLIN distribution favoured the evolution of resistance within local vectors and highlights the potential benefits of resistance management strategies.Trial registration: This study is registered with ISRCTN, ISRCTN17516395. Registered 14 February 2017, http://www.isrctn.com/ISRCTN17516395 .

Insecticide resistance reduces the profitability of insect-resistant rice cultivars

INTRODUCTION

Preventing crop yield loss caused by pests is critical for global agricultural production. Agricultural pest control has largely relied on chemical pesticides. The interaction between insecticide resistance and the adaptation of herbivorous pests to host plants may represent an emerging threat to future food security.

OBJECTIVES

This study aims to unveil genetic evidence for the reduction in the profitability of resistant cultivars derived from insecticide resistance in target pest insects.

METHODS

An experimental evolution system encompassing resistant rice and its major monophagous pest, the brown planthopper Nilaparvata lugens, was constructed. Whole genome resequencing and selective sweep analysis were utilized to identify the candidate gene loci related to the adaptation. RNA interference and induced expression assay were conducted to validate the function of the candidate loci.

RESULTS

We found that the imidacloprid-resistant population of N. lugens rapidly adapted to resistant rice IR36. Gene loci related to imidacloprid resistance may contribute to this phenomenon. Multiple alleles in the nicotinic acetylcholine receptor (nAChR)-7-like and P450 CYP4C61 were significantly correlated with changes in virulence to IR36 rice and insecticide resistance of N. lugens. One avirulent/susceptible genotype and two virulent/resistant genotypes could be inferred from the corresponding alleles. Importantly, we found that the virulent/resistant genotypes already exist in the wild in China, exhibiting increasing frequencies along with insecticide usage. We validated the relevance of these genotypes and the virulence to three more resistant rice cultivars. Knockdown of the above two genes in N. lugens significantly decreased both the resistance to imidacloprid and the virulence towards resistant rice.

CONCLUSION

Our findings provide direct genetic evidence to the eco-evolutionary consequence of insecticide resistance, and suggest an urgent need for the implementation of predictably sustainable pest management.

LLIN Evaluation in Uganda Project (LLINEUP) - Plasmodium infection prevalence and genotypic markers of insecticide resistance in Anopheles vectors from 48 districts of Uganda

Background

In 2017-2019, we conducted a large-scale, cluster-randomised trial (LLINEUP) to evaluate long-lasting insecticidal nets (LLINs) treated with a pyrethroid insecticide plus the synergist piperonyl butoxide (PBO LLINs), as compared to conventional, pyrethroid-only LLINs across 104 health sub-districts (HSDs) in Uganda. In LLINEUP, and similar trials in Tanzania, PBO LLINs were found to provide greater protection against malaria than conventional LLINs, reducing parasitaemia and vector density. In the LLINEUP trial, cross-sectional entomological surveys were carried out at baseline and then every 6 months for two years. In each survey, ten households per HSD were randomly selected for indoor household entomological collections.

Results

Overall, 5395 female Anopheles mosquitoes were collected from 5046 households. The proportion of mosquitoes infected with Plasmodium falciparum did not change significantly over time, while infection with non-falciparum malaria decreased in An. gambiae s.s, but not An. funestus. The frequency of genetic markers associated with pyrethroid resistance increased significantly over time, but the rate of change was not different between the two LLIN types. The knock-down resistance (kdr) mutation Vgsc-995S declined over time as Vgsc-995F, the alternative resistance mutation at this codon, increased. Vgsc-995F appears to be spreading into Uganda.

Conclusions

Distribution of LLINs in Uganda was associated with reductions in parasite prevalence and vector density, but the proportion of infective mosquitoes remained stable, suggesting that the potential for transmission persisted. The increased frequency of markers of pyrethroid resistance indicates that LLIN distribution favoured the evolution of resistance within local vectors and highlights the potential benefits of resistance management strategies.Trial registration:: This study is registered with ISRCTN, ISRCTN17516395. Registered 14 February 2017, http://www.isrctn.com/ISRCTN17516395.

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.

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.

Historical samples reveal a combined role of agriculture and public-health applications in vector resistance to insecticides

BACKGROUND

Insecticide resistance is the major threat to vector control and for the prevention of vector-borne diseases. Because almost all insecticides used against insect vectors are or have been used in agriculture, a connection between agricultural insecticide use and resistance in insect vectors has been hypothesized. However, it is challenging to find a causal link between past agricultural use of insecticides and current resistance in vector populations without historical data series. Here we investigated the relative contribution across time of agricultural and public-health insecticide applications in selecting for diflubenzuron (DFB) resistance in Culex pipiens populations. Using DNA sequencing, we looked for DFB resistant mutations in current and historical mosquito samples, dating back to the 1980s-1990s, when DFB was used in agriculture but not yet in mosquito control.

RESULTS

In the samples collected before the introduction of DFB in vector control, we found the resistant mutation I1043M in rural regions but not any of the neighboring urban and natural areas, indicating that the selection pressure was derived by agriculture. However, after the introduction of DFB for vector control, the resistant mutations were found across all study areas showing that the initial selection from agriculture was further boosted by the selection pressure imposed by the mosquito control applications in the 2000s.

CONCLUSIONS

Our findings support a combined role of agricultural and public-health use of insecticides in vector resistance across time and call for specific actions in integrated resistance management, including increased communication between agriculture and health practitioners. © 2022 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

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.

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

Resistance to pyrethroid insecticides is a major concern for malaria vector control. Pyrethroids target the voltage-gated sodium channel (VGSC), an essential component of the mosquito nervous system. Substitutions in the amino acid sequence can induce a resistance phenotype. We use whole-genome sequence data from phase 2 of the Anopheles gambiae 1000 Genomes Project (Ag1000G) to provide a comprehensive account of genetic variation in the Vgsc gene across 13 African countries. In addition to known resistance alleles, we describe 20 other non-synonymous nucleotide substitutions at appreciable population frequency and map these variants onto a protein model to investigate the likelihood of pyrethroid resistance phenotypes. Thirteen of these novel alleles were found to occur almost exclusively on haplotypes carrying the known L995F kdr (knock-down resistance) allele and may enhance or compensate for the L995F resistance genotype. A novel mutation I1527T, adjacent to a predicted pyrethroid-binding site, was found in tight linkage with V402L substitutions, similar to allele combinations associated with resistance in other insect species. We also analysed genetic backgrounds carrying resistance alleles, to determine which alleles have experienced recent positive selection, and describe ten distinct haplotype groups carrying known kdr alleles. Five of these groups are observed in more than one country, in one case separated by over 3000 km, providing new information about the potential for the geographical spread of resistance. Our results demonstrate that the molecular basis of target-site pyrethroid resistance in malaria vectors is more complex than previously appreciated, and provide a foundation for the development of new genetic tools for insecticide resistance management.

Pyrethroid Resistance in Anopheles gambiae Not Associated with Insecticide-Treated Mosquito Net Effectiveness Across Sub-Saharan Africa

Pyrethroid resistance is a major concern for malaria vector control programs that predominantly rely on insecticide-treated mosquito nets (ITNs). Contradictory results of the impact of resistance have been observed during field studies. We combined continent-wide estimates of pyrethroid resistance in Anopheles gambiae from 2006 to 2017, with continent-wide survey data to assess the effect of increasing pyrethroid resistance on the effectiveness of ITNs to prevent malaria infections in sub-Saharan Africa. We used a pooled-data approach and a meta-regression of survey regions to assess how pyrethroid resistance affects the association between ITN ownership and malaria outcomes for children 6 to 59 months of age. ITN ownership reduced the risk of malaria outcomes according to both the pooled and meta-regression approaches. According to the pooled analysis, there was no observed interaction between ITN ownership and estimated level of pyrethroid resistance (likelihood ratio [LR] test, 1.127 for malaria infection confirmed by the rapid diagnostic test, P = 0.2885; LR test = 0.161 for microscopy-confirmed malaria infection, P = 0.161; LR test = 0.646 for moderate or severe anemia, P = 0.4215). Using the meta-regression approach to determine the level of pyrethroid resistance did not explain any of the variance in subnational estimates of ITN effectiveness for any of the outcomes. ITNs decreased the risk of malaria independent of the levels of pyrethroid resistance in malaria vector populations.

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.

First Report of Target Site Insensitivity in Pyrethroid Resistant Anopheles gambiae from Southern Guinea Savanna, Northern-Nigeria

Background:

Malaria is a major public health problem and life threatening parasitic vector-borne disease. For the first time, we established and report the molecular mechanism responsible for Anopheles gambiae s.l. resistance to pyrethroids and DDT from Yamaltu Deba, Southern Guinea Savanna, Northern-Nigeria.

Methods:

The susceptibility profile of An. gambiae s.l. to four insecticides (DDT 4%, bendiocarb 0.1%, malathion 5% and deltamethrin 0.05%) using 2–3 days old females from larvae collected from study area between August and November, 2018 was first established. Genomic DNA was then extracted from 318 mosquitoes using Livak DNA extraction protocol for specie identification and kdr genotyping. The mosquitoes were identified to species level and then 96 genotyped for L1014F and L1014S kdr target site mutations.

Results:

The mosquitoes were all resistant to DDT, bendiocarb and deltamethrin but fully susceptible to malathion. An. coluzzii was found to be the dominant sibling species (97.8%) followed by An. arabiensis (1.9%) and An. gambiae s.s (0.3%). The frequency of the L1014F kdr mutation was relatively higher (83.3%) than the L1014S (39%) in the three species studied. The L1014F showed a genotypic frequency of 75% resistance (RR), 17% heterozygous (RS) and 8% susceptible (SS) with an allelic frequency of 87% RR and 13% SS while the L1014S showed a genotypic frequency of RR (16%), RS (38%) and SS (46%) with an allelic frequency of 40% RR and 60% SS, respectively.

Conclusion:

This study reveals that both kdr mutations present simultaneously in Northern-Nigeria, however contribution of L1014F which is common in West Africa was more than twice of L1014S mutation found in East Africa.

Detection of kdr and ace-1 mutations in wild populations of Anopheles arabiensis and An. melas in a residual malaria transmission area of Senegal

In the central western Senegal, malaria transmission has been reduced low due to the combination of several effective control interventions. However, despite this encouraging achievement, residual malaria transmission still occurring in few areas, mainly ensured by An. arabiensis and An. melas. The resurgence or the persistence of the disease may have originated from the increase and the spread of insecticide resistance genes among natural malaria vectors populations. Therefore, assessing the status and mechanisms of insecticides resistance among targeted malaria vectors is of highest importance to better characterize factors underlying the residual transmission where it occurs. Malaria vectors were collected from three selected villages using nocturnal human landing catches (HLC) and pyrethrum spray collections (PSC) methods. An. gambiae s.l. specimens were identified at the species level then genotyped for the presence of kdr-west (L1014F), kdr-east (L1014S) and ace-1^R mutations by qPCR. An. arabiensis (69.36%) and An. melas (27.99%) were the most common species of the Gambiae complex in the study area. Among An. arabiensis population, the allelic frequency of the kdr-east (22.66%) was relatively higher than for kdr-west mutation (9.96%). While for An. melas populations, the overall frequencies of both mutations were very low, being respectively 1.12% and 0.40% for the L1014S and L1014F mutations. With a global frequency of 2%, only the heterozygous form of the G119S mutation was found only in An. arabiensis and in all the study sites. The widespread occurrence of the kdr mutation in both An. arabiensis and An. melas natural populations, respectively the main and focal vectors in the central-western Senegal, may have contributed to maintaining malaria transmission in the area. Thus, compromising the effectiveness of pyrethroids-based vector control measures and the National Elimination Goal. Therefore, monitoring and managing properly insecticide resistance became a key programmatic intervention to achieve the elimination goal where feasible, as aimed by Senegal. Noteworthy, this is the first report of the ace-1 mutation in natural populations of An. arabiensis from Senegal, which need to be closely monitored to preserve one of the essential insecticide classes used in IRS to control the pyrethroids-resistant populations.

Evolution of the Insecticide Target Rdl in African Anopheles Is Driven by Interspecific and Interkaryotypic Introgression

The evolution of insecticide resistance mechanisms in natural populations of Anopheles malaria vectors is a major public health concern across Africa. Using genome sequence data, we study the evolution of resistance mutations in the resistance to dieldrin locus (Rdl), a GABA receptor targeted by several insecticides, but most notably by the long-discontinued cyclodiene, dieldrin. The two Rdl resistance mutations (296G and 296S) spread across West and Central African Anopheles via two independent hard selective sweeps that included likely compensatory nearby mutations, and were followed by a rare combination of introgression across species (from A. gambiae and A. arabiensis to A. coluzzii) and across nonconcordant karyotypes of the 2La chromosomal inversion. Rdl resistance evolved in the 1950s as the first known adaptation to a large-scale insecticide-based intervention, but the evolutionary lessons from this system highlight contemporary and future dangers for management strategies designed to combat development of resistance in malaria vectors.

The need for new vector control approaches targeting outdoor biting Anopheline malaria vector communities

Since the implementation of Roll Back Malaria, the widespread use of insecticide-treated nets (ITNs) and indoor residual spraying (IRS) is thought to have played a major part in the decrease in mortality and morbidity achieved in malaria-endemic regions. In the past decade, resistance to major classes of insecticides recommended for public health has spread across many malaria vector populations. Increasingly, malaria vectors are also showing changes in vector behaviour in response to current indoor chemical vector control interventions. Changes in the time of biting and proportion of indoor biting of major vectors, as well as changes in the species composition of mosquito communities threaten the progress made to control malaria transmission. Outdoor biting mosquito populations contribute to malaria transmission in many parts of sub-Saharan Africa and pose new challenges as they cannot be reliably monitored or controlled using conventional tools. Here, we review existing and novel approaches that may be used to target outdoor communities of malaria vectors. We conclude that scalable tools designed specifically for the control and monitoring of outdoor biting and resting malaria vectors with increasingly complex and dynamic responses to intensifying malaria control interventions are urgently needed. These are crucial for integrated vector management programmes designed to challenge current and future vector populations.

Resting behaviour of malaria vectors in highland and lowland sites of western Kenya: Implication on malaria vector control measures

BACKGROUND

Understanding the interactions between increased insecticide resistance and resting behaviour patterns of malaria mosquitoes is important for planning of adequate vector control. This study was designed to investigate the resting behavior, host preference and rates of Plasmodium falciparum infection in relation to insecticide resistance of malaria vectors in different ecologies of western Kenya.

METHODS

Anopheles mosquito collections were carried out during the dry and rainy seasons in Kisian (lowland site) and Bungoma (highland site), both in western Kenya using pyrethrum spray catches (PSC), mechanical aspiration (Prokopack) for indoor collections, clay pots, pit shelter and Prokopack for outdoor collections. WHO tube bioassay was used to determine levels of phenotypic resistance of indoor and outdoor collected mosquitoes to deltamethrin. PCR-based molecular diagnostics were used for mosquito speciation, genotype for knockdown resistance mutations (1014S and 1014F) and to determine specific host blood meal origins. Enzyme-linked Immunosorbent Assay (ELISA) was used to determine mosquito sporozoite infections.

RESULTS

Anopheles gambiae s.l. was the most predominant species (75%, n = 2706) followed by An. funestus s.l. (25%, n = 860). An. gambiae s.s hereafter (An. gambiae) accounted for 91% (95% CI: 89-93) and An. arabiensis 8% (95% CI: 6-9) in Bungoma, while in Kisian, An. arabiensis composition was 60% (95% CI: 55-66) and An. gambiae 39% (95% CI: 34-44). The resting densities of An. gambiae s.l and An. funestus were higher indoors than outdoor in both sites (An. gambiae s.l; F1, 655 = 41.928, p < 0.0001, An. funestus; F1, 655 = 36.555, p < 0.0001). The mortality rate for indoor and outdoor resting An. gambiae s.l F1 progeny was 37% (95% CI: 34-39) vs 67% (95% CI: 62-69) respectively in Bungoma. In Kisian, the mortality rate was 67% (95% CI: 61-73) vs 76% (95% CI: 71-80) respectively. The mortality rate for F1 progeny of An. funestus resting indoors in Bungoma was 32% (95% CI: 28-35). The 1014S mutation was only detected in indoor resitng An. arabiensis. Similarly, the 1014F mutation was present only in indoor resting An. gambiae. The sporozoite rates were highest in An. funestus followed by An. gambiae, and An. arabiensis resting indoors at 11% (34/311), 8% (47/618) and 4% (1/27) respectively in Bungoma. Overall, in Bungoma, the sporozoite rate for indoor resting mosquitoes was 9% (82/956) and 4% (8/190) for outdoors. In Kisian, the sporozoite rate was 1% (1/112) for indoor resting An. gambiae. None of the outdoor collected mosquitoes in Kisian tested positive for sporozoite infections (n = 73).

CONCLUSION

The study reports high indoor resting densities of An. gambiae and An. funestus, insecticide resistance, and persistence of malaria transmission indoors regardless of the use of long-lasting insecticidal nets (LLINs). These findings underline the difficulties of controlling malaria vectors resting and biting indoors using the current interventions. Supplemental vector control tools and implementation of sustainable insecticide resistance management strategies are needed in western Kenya.

Status of Insecticide Resistance and Its Mechanisms in Anopheles gambiae and Anopheles coluzzii Populations from Forest Settings in South Cameroon

A key factor affecting malaria vector control efforts in Cameroon is the rapid expansion of insecticide resistance in Anopheles gambiae s.l (An. gambiae) populations; however, mechanisms involved in insecticide resistance in forest mosquito populations are still not well documented yet. The present study was conducted to screen molecular mechanisms conferring insecticide resistance in An. gambiae s.l. populations from the South Cameroon forest region. WHO bioassays were conducted with F0 An. gambiae females aged three to four days from forest (Sangmelima, Nyabessan, and Mbandjock) and urban sites (Yaoundé (Bastos and Nkolondom)), against pyrethroids (permethrin 0.75% and deltamethrin 0.05%) and carbamates (bendiocarb 0.1%). Members of the An. Gambiae s.l. species complex were identified using molecular diagnostic tools. TaqMan assays were used to screen for target site mutations. The expression profiles of eight genes implicated in insecticide resistance were assessed using RT-qPCR. Cuticle hydrocarbon lipids were measured to assess their potential implication in insecticide resistance. Both An. Gambiae and An. coluzzii were detected. An. gambiae was highly prevalent in Sangmelima, Nyabessan, Mbandjock, and Nkolondom. An. coluzzii was the only species found in the Yaoundé city center (Bastos). Low mortality rate to both pyrethroids and bendiocarb was recorded in all sites. High frequency of L1014F allele (75.32–95.82%) and low frequencies of L1014S (1.71–23.05%) and N1575Y (5.28–12.87%) were recorded. The G119S mutation (14.22–35.5%) was detected for the first time in An. gambiae populations from Cameroon. This mutation was rather absent from An. coluzzii populations. The detoxification genes Cyp6m2, Cyp9k1, Cyp6p4, Cyp6z1, as well as Cyp4g16 which catalyzes epicuticular hydrocarbon biosynthesis, were found to be overexpressed in at least one population. The total cuticular hydrocarvbon content, a proxy of cuticular resistance, did not show a pattern associated with pyrethroid resistance in these populations. The rapid emergence of multiple resistance mechanisms in An. Gambiae s.l. population from the South Cameroon forest region is of big concern and could deeply affect the sustainability of insecticide-based interventions strategies in this region.

A Re-Evaluation of Olive Fruit Fly Organophosphate-Resistant Ace Alleles in Iberia, and Field-Testing Population Effects after in-Practice Dimethoate Use

The management of the olive fruit fly (Bactrocera oleae) is traditionally based upon the use of organophosphate insecticides, mainly dimethoate. In this evolutionary arms race between man and pest, the flies have adapted a pesticide resistance, implying two point-mutations of the Ace gene -I214V and G488S- and a 9bp deletion -Δ3Q. We revisited 11 Iberian locations to evaluate this adaptation of organophosphate (OP)-resistant alleles through amplicon sequencing. Screening for populations where the wild type is prevalent allows an identification of hotspots for targeted mitigation measures; we have hence refined the scale to the region with the lowest OP-resistant alleles frequency 71 locations were sampled and individuals checked using a fast and low-cost allele-specific-primer polymerase chain reaction (ASP-PCR) method]. An increase in Ace gene point-mutations was observed, and the Δ3Q mutation remains undetected. The lowest frequencies of the OP-resistant alleles remain in the west, underlining the hypothesis of an introduction of resistance from eastern Mediterranean areas. A field test was performed by sampling the fly population before and after in-practice dimethoate application. A clear reduction in olive fruit fly numbers was observed, with no relevant changes in the genotypic frequencies of the resistance alleles. The findings are discussed in frame of the type and intensity of the selection pressure that has led to the adaptation to resistance and its consequences from the producer perspective.

The evolutionary origins of pesticide resistance

Durable crop protection is an essential component of current and future food security. However, the effectiveness of pesticides is threatened by the evolution of resistant pathogens, weeds and insect pests. Pesticides are mostly novel synthetic compounds, and yet target species are often able to evolve resistance soon after a new compound is introduced. Therefore, pesticide resistance provides an interesting case of rapid evolution under strong selective pressures, which can be used to address fundamental questions concerning the evolutionary origins of adaptations to novel conditions. We ask: (i) whether this adaptive potential originates mainly from de novo mutations or from standing variation; (ii) which pre-existing traits could form the basis of resistance adaptations; and (iii) whether recurrence of resistance mechanisms among species results from interbreeding and horizontal gene transfer or from independent parallel evolution. We compare and contrast the three major pesticide groups: insecticides, herbicides and fungicides. Whilst resistance to these three agrochemical classes is to some extent united by the common evolutionary forces at play, there are also important differences. Fungicide resistance appears to evolve, in most cases, by de novo point mutations in the target-site encoding genes; herbicide resistance often evolves through selection of polygenic metabolic resistance from standing variation; and insecticide resistance evolves through a combination of standing variation and de novo mutations in the target site or major metabolic resistance genes. This has practical implications for resistance risk assessment and management, and lessons learnt from pesticide resistance should be applied in the deployment of novel, non-chemical pest-control methods.

Adaptive Introgression: An Untapped Evolutionary Mechanism for Crop Adaptation

Global environmental changes strongly impact wild and domesticated species biology and their associated ecosystem services. For crops, global warming has led to significant changes in terms of phenology and/or yield. To respond to the agricultural challenges of this century, there is a strong need for harnessing the genetic variability of crops and adapting them to new conditions. Gene flow, from either the same species or a different species, may be an immediate primary source to widen genetic diversity and adaptions to various environments. When the incorporation of a foreign variant leads to an increase of the fitness of the recipient pool, it is referred to as "adaptive introgression". Crop species are excellent case studies of this phenomenon since their genetic variability has been considerably reduced over space and time but most of them continue exchanging genetic material with their wild relatives. In this paper, we review studies of adaptive introgression, presenting methodological approaches and challenges to detecting it. We pay particular attention to the potential of this evolutionary mechanism for the adaptation of crops. Furthermore, we discuss the importance of farmers' knowledge and practices in shaping wild-to-crop gene flow. Finally, we argue that screening the wild introgression already existing in the cultivated gene pool may be an effective strategy for uncovering wild diversity relevant for crop adaptation to current environmental changes and for informing new breeding directions.

Insecticide resistance in Anopheles gambiae from the northern Democratic Republic of Congo, with extreme knockdown resistance (kdr) mutation frequencies revealed by a new diagnostic assay

BACKGROUND

Mutations in the voltage-gated sodium channel at codon 1014 confer knock-down resistance (kdr) to pyrethroids in a wide range of insects. Anopheles gambiae exhibits two mutant alleles at codon 1014, serine and phenylalanine; and both are now widespread across Africa. Existing screening methods only allow for one resistant allele to be detected per assay. A new locked nucleic acid (LNA) qPCR assay was developed for the simultaneous detection of both mutant alleles and the wild type allele in a single assay. This tri-allelic detection assay was assessed as part of a study of the insecticide resistance in An. gambiae sensu stricto (s.s.) in the previously un-sampled area of Nord Ubangi, Democratic Republic of the Congo.

METHODS

Samples from three sites were tested for insecticide susceptibility using WHO bioassays, with and without the synergist PBO preceding pyrethroid exposures, and were subsequently analysed for frequency and resistance-association of the Vgsc-1014 and Vgsc-N1575Y mutations. Results from the LNA-kdr 1014 assay were compared to results from standard TaqMan-kdr assays.

RESULTS

Anopheles gambiae sensu lato (s.l.) was by far the predominant vector captured (84%), with only low frequencies of Anopheles funestus s.l. (9%) detected in Nord Ubangi. Molecular identification found An. gambiae s.s. to be the principal vector (99%) although Anopheles coluzzii was detected at very low frequency. Anopheles gambiae were susceptible to the carbamate insecticide bendiocarb, but resistant to DDT and to the pyrethroids permethrin and deltamethrin. Susceptibility to both pyrethroids was partially restored with prior exposure to PBO suggesting likely involvement of metabolic resistance. Anopheles gambiae s.s. was homozygous for kdr resistant alleles with both the L1014F and L1014S mutations present, and the N1575Y polymorphism was present at low frequency. The LNA-kdr assay simultaneously detected both resistant alleles and gave results entirely consistent with those from the two TaqMan-kdr assays.

CONCLUSION

This study provides rare data on insecticide resistance and mechanisms in Anopheles from the centre of Africa, with the first detection of N1575Y. Nord Ubangi populations of An. gambiae s.s. show insecticide resistance mediated by both metabolic mechanisms and Vgsc mutations. The LNA-kdr assay is particularly suitable for use in populations in which both 1014S and 1014F kdr alleles co-occur and provides robust results, with higher throughput and at a quarter of the cost of TaqMan assays.

Development of a rapid field-applicable molecular diagnostic for knockdown resistance (kdr) markers in An. gambiae

Background

The spread of insecticide resistance (IR) is a major threat to vector control programmes for mosquito-borne diseases. Early detection of IR using diagnostic markers could help inform these programmes, especially in remote locations where gathering reliable bioassay data is challenging. Most current molecular tests for genetic IR markers are only suitable for use in well-equipped laboratory settings. There is an unmet need for field-applicable diagnostics.

Methods

A single-cartridge test was designed to detect key IR mutations in the major African vector of malaria, Anopheles gambiae. Developed on the portable, rapid, point-of-care compatible PCR platform - Genedrive® (genedrive® plc), the test comprises two assays which target single nucleotide polymorphisms (SNPs) in the voltage gated sodium channel (VGSC) gene that exert interactive effects on knockdown resistance (kdr): L1014F, L1014S and N1575Y.

Results

Distinct melt peaks were observed for each allele at each locus. Preliminary validation of these assays using a test panel of 70 An. gambiae samples showed complete agreement of our assays with the widely-used TaqMan assays, achieving a sensitivity and specificity of 100%.

Conclusion

Here we show the development of an insecticide resistance detection assay for use on the Genedrive® platform that has the potential to be the first field-applicable diagnostic for kdr.

Insecticide resistance in Anopheles arabiensis populations from Dakar and its suburbs: role of target site and metabolic resistance mechanisms

Background

Urban malaria is an increasing concern in most of the sub-Saharan Africa countries. In Dakar, the capital city of Senegal, the malaria epidemiology has been complicated by recurrent flooding since 2005. The main vector control measure for malaria prevention in Dakar is the community use of long-lasting insecticide-treated nets. However, the increase of insecticide resistance reported in this area needs to be better understood for suitable resistance management. This study reports the situation of insecticide resistance and underlying mechanisms in Anopheles arabiensis populations from Dakar and its suburbs.

Results

All the populations tested showed resistance to almost all insecticides except organophosphates families, which remain the only lethal molecules. Piperonil butoxide (PBO) and ethacrinic acid (EA) the two synergists used, have respectively and significantly restored the susceptibility to DDT and permethrin of Anopheles population. Molecular identification of specimens revealed the presence of An. arabiensis only. Kdr genotyping showed the presence of the L1014F mutation (kdr-West) as well as L1014S (kdr-East). This L1014S mutation was found at very high frequencies (89.53%) in almost all districts surveyed, and in association with the L1014F (10.24%).

Conclusion

Results showed the contribution of both target-site and metabolic mechanisms in conferring pyrethroid resistance to An. arabiensis from the flooded areas of Dakar suburbs. These data, although preliminary, stress the need for close monitoring of the urban An. arabiensis populations for a suitable insecticide resistance management system to preserve core insecticide-based vector control tools in this flooded area.

Candidate-gene based GWAS identifies reproducible DNA markers for metabolic pyrethroid resistance from standing genetic variation in East African Anopheles gambiae

Metabolic resistance to pyrethroid insecticides is widespread in Anopheles mosquitoes and is a major threat to malaria control. DNA markers would aid predictive monitoring of resistance, but few mutations have been discovered outside of insecticide-targeted genes. Isofemale family pools from a wild Ugandan Anopheles gambiae population, from an area where operational pyrethroid failure is suspected, were genotyped using a candidate-gene enriched SNP array. Resistance-associated SNPs were detected in three genes from detoxification superfamilies, in addition to the insecticide target site (the Voltage Gated Sodium Channel gene, Vgsc). The putative associations were confirmed for two of the marker SNPs, in the P450 Cyp4j5 and the esterase Coeae1d by reproducible association with pyrethroid resistance in multiple field collections from Uganda and Kenya, and together with the Vgsc-1014S (kdr) mutation these SNPs explained around 20% of variation in resistance. Moreover, the >20 Mb 2La inversion also showed evidence of association with resistance as did environmental humidity. Sequencing of Cyp4j5 and Coeae1d detected no resistance-linked loss of diversity, suggesting selection from standing variation. Our study provides novel, regionally-validated DNA assays for resistance to the most important insecticide class, and establishes both 2La karyotype variation and humidity as common factors impacting the resistance phenotype.

Emerging Pyrethroid Resistance among Anopheles arabiensis in Kenya

Vector control programs, particularly in the form of insecticide-treated bed nets (ITNs), are essential for achieving malaria elimination goals. Recent reports of increasing knockdown resistance (kdr) mutation frequencies for Anopheles arabiensis in Western Kenya heightens the concern on the future effectiveness of ITNs in Kenya. We examined resistance in An. arabiensis populations across Kenya through kdr mutations and World Health Organization–recommended bioassays. We detected two kdr alleles, L1014F and L1014S. Kdr mutations were found in five of the 11 study sites, with mutation frequencies ranging from 3% to 63%. In two Western Kenya populations, the kdr L1014F allele frequency was as high as 10%. The L1014S frequency was highest at Chulaimbo at 55%. Notably, the kdr L1014F mutation was found to be associated with pyrethroid resistance at Port Victoria, but kdr mutations were not significantly associated with resistance at Chulaimbo, which had the highest kdr mutation frequency among all sites. This study demonstrated the emerging pyrethroid resistance in An. arabiensis and that pyrethroid resistance may be related to kdr mutations. Resistance monitoring and management are urgently needed for this species in Kenya where resistance is emerging and its abundance is becoming predominant. Kdr mutations may serve as a biomarker for pyrethroid resistance in An. arabiensis.

Human Interventions: Driving Forces of Mosquito Evolution

One of the most common strategies for controlling mosquito-borne diseases relies on the use of chemical pesticides to repel or kill the mosquito vector. Pesticide exposure interferes with several key biological functions in the mosquito and triggers a variety of adaptive responses whose underlying mechanisms are only partially elucidated. The availability of reference genome sequences opens up the possibility of tracking signatures of evolutionary changes, including the most recent, across the genomes of many vector species. In this review, we highlight the recent genomic changes, which contribute to the fascinating adaptation of malaria vectors of the sub-Saharan African region to intensive insecticide-based interventions. We emphasize the operational significance of detailed genomic knowledge for current monitoring and decision-making.

Malaria Vector Control Still Matters despite Insecticide Resistance

Mosquito vectors' resistance to insecticides is usually considered a major threat to the recent progresses in malaria control. However, studies measuring the impact of interventions and insecticide resistance reveal inconsistencies when using entomological versus epidemiological indices. First, evaluation tests that do not reflect the susceptibility of mosquitoes when they are infectious may underestimate insecticide efficacy. Moreover, interactions between insecticide resistance and vectorial capacity reveal nonintuitive outcomes of interventions. Therefore, considering ecological interactions between vector, parasite, and environment highlights that the impact of insecticide resistance on the malaria burden is not straightforward and we suggest that vector control still matters despite insecticide resistance.

Knockdown resistance mutations predict DDT resistance and pyrethroid tolerance in the visceral leishmaniasis vector Phlebotomus argentipes

BACKGROUND

Indoor residual spraying (IRS) with DDT has been the primary strategy for control of the visceral leishmaniasis (VL) vector Phlebotomus argentipes in India but efficacy may be compromised by resistance. Synthetic pyrethroids are now being introduced for IRS, but with a shared target site, the para voltage-gated sodium channel (VGSC), mutations affecting both insecticide classes could provide cross-resistance and represent a threat to sustainable IRS-based disease control.

METHODOLOGY/PRINCIPAL FINDINGS

A region of the Vgsc gene was sequenced in P. argentipes from the VL hotspot of Bihar, India. Two knockdown resistance (kdr) mutations were detected at codon 1014 (L1014F and L1014S), each common in mosquitoes, but previously unknown in phlebotomines. Both kdr mutations appear largely recessive, but as homozygotes (especially 1014F/F) or as 1014F/S heterozygotes exert a strong effect on DDT resistance, and significantly predict survivorship to class II pyrethroids in short-duration bioassays. The mutations are present at high frequency in wild P. argentipes populations from Bihar, with 1014F significantly more common in higher VL areas.

CONCLUSIONS/SIGNIFICANCE

The Vgsc mutations detected appear to be a primary mechanism underlying DDT resistance in P. argentipes and a contributory factor in reduced pyrethroid susceptibility, suggesting a potential impact if P. argentipes are subjected to suboptimal levels of pyrethroid exposure, or additional resistance mechanisms evolve. The assays to detect kdr frequency changes provide a sensitive, high-throughput monitoring tool to detecting spatial and temporal variation in resistance in P. argentipes.

A significant association between deltamethrin resistance, Plasmodium falciparum infection and the Vgsc-1014S resistance mutation in Anopheles gambiae highlights the epidemiological importance of resistance markers

Background

The success of malaria vector control is threatened by widespread pyrethroid insecticide resistance. However, the extent to which insecticide resistance impacts transmission is unclear. The objective of this study was to examine the association between the DDT/pyrethroid knockdown resistance mutation Vgsc-1014S, commonly termed kdr, and infection with Plasmodium falciparum sporozoites in Anopheles gambiae.

Methods

WHO standard methods were used to characterize susceptibility of wild female mosquitoes to 0.05 % deltamethrin. PCR-based molecular diagnostics were used to identify mosquitoes to species and to genotype at the Vgsc-L1014S locus. ELISAs were used to detect the presence of P.falciparum sporozoites and for blood meal identification.

Results

Anopheles mosquitoes were resistant to deltamethrin with mortality rates of 77.7 % [95 % CI 74.9–80.3 %]. Of 545 mosquitoes genotyped 96.5 % were A. gambiaes.s. and 3.5 % were Anopheles arabiensis. The Vgsc-1014S mutation was detected in both species. Both species were predominantly anthropophagic. In A. gambiaes.s., Vgsc-L1014S genotype was significantly associated with deltamethrin resistance (χ2 = 11.2; p < 0.001). The P. falciparum sporozoite infection rate was 4.2 %. There was a significant association between the presence of sporozoites and Vgsc-L1014S genotype in A. gambiaes.s. (χ2 = 4.94; p = 0.026).

Conclusions

One marker, Vgsc-1014S, was associated with insecticide resistance and P. falciparum infection in wild-caught mixed aged populations of A. gambiaes.s. thereby showing how resistance may directly impact transmission.

Malaria vector populations across ecological zones in Guinea Conakry and Mali, West Africa

Background

Malaria remains a pervasive public health problem in sub-Saharan West Africa. Here mosquito vector populations were explored across four sites in Mali and the Republic of Guinea (Guinea Conakry). The study samples the major ecological zones of malaria-endemic regions in West Africa within a relatively small distance.

Methods

Mosquito vectors were sampled from larval pools, adult indoor resting sites, and indoor and outdoor human-host seeking adults. Mosquitoes were collected at sites spanning 350 km that represented arid savannah, humid savannah, semi-forest and deep forest ecological zones, in areas where little was previously known about malaria vector populations. 1425 mosquito samples were analysed by molecular assays to determine species, genetic attributes, blood meal sources and Plasmodium infection status.

Results

Anopheles gambiae and Anopheles coluzzii were the major anophelines represented in all collections across the ecological zones, with A. coluzzii predominant in the arid savannah and A. gambiae in the more humid sites. The use of multiple collection methodologies across the sampling sites allows assessment of potential collection bias of the different methods. The L1014F kdr insecticide resistance mutation (kdr-w) is found at high frequency across all study sites. This mutation appears to have swept almost to fixation, from low frequencies 6 years earlier, despite the absence of widespread insecticide use for vector control. Rates of human feeding are very high across ecological zones, with only small fractions of animal derived blood meals in the arid and humid savannah. About 30 % of freshly blood-fed mosquitoes were positive for Plasmodium falciparum presence, while the rate of mosquitoes with established infections was an order of magnitude lower.

Conclusions

The study represents detailed vector characterization from an understudied area in West Africa with endemic malaria transmission. The deep forest study site includes the epicenter of the 2014 Ebola virus epidemic. With new malaria control interventions planned in Guinea, these data provide a baseline measure and an opportunity to assess the outcome of future interventions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-016-1242-5) contains supplementary material, which is available to authorized users.

Identification, Validation, and Application of Molecular Diagnostics for Insecticide Resistance in Malaria Vectors

Insecticide resistance is a major obstacle to control of Anopheles malaria mosquitoes in sub-Saharan Africa and requires an improved understanding of the underlying mechanisms. Efforts to discover resistance genes and DNA markers have been dominated by candidate gene and quantitative trait locus studies of laboratory strains, but with greater availability of genome sequences a shift toward field-based agnostic discovery is anticipated. Mechanisms evolve continually to produce elevated resistance yielding multiplicative diagnostic markers, co-screening of which can give high predictive value. With a shift toward prospective analyses, identification and screening of resistance marker panels will boost monitoring and programmatic decision making.

Diversity, resistance and vector competence of endophilic anophelines from southern Ghana

Background

As part of efforts to monitor the impact of vector control strategies so that they can be improved and more targeted, we collected baseline data on aspects of the bionomics of endophilic anophelines in southern Ghana.

Materials and Methods

Indoor resting anophelines were collected using mouth aspirators and pyrethroid spray catch. Anopheles females were identified to species level using morphological characteristics and sibling species were distinguished by PCR. The presence of the L1014F mutation, conferring resistance to insecticides, was determined in An. gambiae s.s. and An. coluzzii samples using TaqMan real-time PCR. Host blood meal sources were determined by PCR, and the presence of Plasmodium falciparum circumsporozoite proteins determined by ELISA.

Results

A total of 892 female Anopheles (31% An. gambiae, 41% An. coluzzii and 28% An. funestus) were collected from six villages. The L1014F mutation was almost fixed in all populations studied (allele frequencies: 0.87-1.00). Both An. gambiae s.l. and An. funestus fed mainly on humans, with a human blood index of 1, although some animal feeding was recorded in An. gambiae. P. falciparum was detected in all ecological zones and in all three major vector species, being 4.9% in An. funestus, 3.8% in An. gambiae s.s. and 1.1% in An. coluzzii.

Conclusions

These findings suggest that the three major vectors of malaria are present in all ecological zones of southern Ghana and contribute to disease transmission. The near fixation of the L1014F mutation in southern Ghana poses a great threat to vector control, thus highlighting the urgent need to implement measures to maintain the efficacy of current control tools and to develop novel control strategies.

The kdr-bearing haplotype and susceptibility to Plasmodium falciparum in Anopheles gambiae: genetic correlation and functional testing

Background

Members of the Anophelesgambiae species complex are primary vectors of human malaria in Africa. It is known that a large haplotype shared between An. gambiae and Anophelescoluzzii by introgression carries point mutations of the voltage-gated sodium channel gene para, including the L1014F kdr mutation associated with insensitivity to pyrethroid insecticides. Carriage of L1014F kdr is also correlated with higher susceptibility to infection with Plasmodium falciparum. However, the genetic mechanism and causative gene(s) underlying the parasite susceptibility phenotype are not known.

Methods

Mosquitoes from the wild Burkina Faso population were challenged by feeding on natural P. falciparum gametocytes. Oocyst infection phenotypes were determined and were tested for association with SNP genotypes. Candidate genes in the detected locus were prioritized and RNAi-mediated gene silencing was used to functionally test for gene effects on P. falciparum susceptibility.

Results

A genetic locus, Pfin6, was identified that influences infection levels of P. falciparum in mosquitoes. The locus segregates as a ~3 Mb haplotype carrying 65 predicted genes including the para gene. The haplotype carrying the kdr allele of para is linked to increased parasite infection prevalence, but many single nucleotide polymorphisms on the haplotype are also equally linked to the infection phenotype. Candidate genes in the haplotype were prioritized and functionally tested. Silencing of para did not influence P. falciparum infection, while silencing of a predicted immune gene, serine protease ClipC9, allowed development of significantly increased parasite numbers.

Conclusions

Genetic variation influencing Plasmodium infection in wild Anopheles is linked to a natural ~3 megabase haplotype on chromosome 2L that carries the kdr allele of the para gene. Evidence suggests that para gene function does not directly influence parasite susceptibility, and the association of kdr with infection may be due to tight linkage of kdr with other gene(s) on the haplotype. Further work will be required to determine if ClipC9 influences the outcome of P. falciparum infection in nature, as well as to confirm the absence of a direct influence by para.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-015-0924-8) contains supplementary material, which is available to authorized users.

A genetic model of the effects of insecticide-treated bed nets on the evolution of insecticide-resistance

BACKGROUND AND OBJECTIVES

The evolution of insecticide-resistance in malaria vectors is emerging as a serious challenge for the control of malaria. Modelling the spread of insecticide-resistance is an essential tool to understand the evolutionary pressures and dynamics caused by the application of insecticides.

METHODOLOGY

We developed a population-genetic model of the spread of insecticide-resistance in a population of Anopheles vectors in response to insecticides used either as adulticides (focussing on insecticide-treated bed nets (ITNs)) or as larvicides (either for the control of malaria or, as an inadvertent side-product, in agriculture).

RESULTS

We show that indoor use of insecticides leads to considerably less selection pressure than their use as larvicides, supporting the idea that most resistance of malaria vectors is due to the agricultural use of the insecticides that are also used for malaria control. The reasons for the relatively low selection pressure posed by adulticides are (i) that males are not affected by the ITNs and, in particular, (ii) that the insecticides are also repellents, keeping mosquitoes at bay from contacting the insecticide but also driving them to bite either people who do not use the insecticide or alternative hosts.

CONCLUSION

We conclude by discussing the opposing public health benefits of high repellency at an epidemiological and an evolutionary timescale: whereas repellency is beneficial to delay the evolution of resistance, other models have shown that it decreases the population-level protection of the insecticide.

Incongruent nuclear and mitochondrial genetic structure of new world screwworm fly populations due to positive selection of mutations associated with dimethyl- and diethyl-organophosphates resistance

Livestock production is an important economic activity in Brazil, which has been suffering significant losses due to the impact of parasites. The New World screwworm (NWS) fly, Cochliomyia hominivorax, is an ectoparasite and one of the most important myiasis-causing flies endemic to the Americas. The geographic distribution of NWS has been reduced after the implementation of the Sterile Insect Technique (SIT), being eradicated in North America and part of Central America. In South America, C. hominivorax is controlled by chemical insecticides, although indiscriminate use can cause selection of resistant individuals. Previous studies have associated the Gly137Asp and Trp251Leu mutations in the active site of carboxylesterase E3 to resistance of diethyl and dimethyl-organophosphates insecticides, respectively. Here, we have sequenced a fragment of the carboxylesterase E3 gene (ChαE7), comprising part of intron iII, exon eIII, intron iIII and part of exon eIV, and three mitochondrial gene sequences (CR, COI and COII), of NWS flies from 21 locations in South America. These markers were used for population structure analyses and the ChαE7 gene was also investigated to gain insight into the selective pressures that have shaped its evolution. Analysis of molecular variance (AMOVA) and pairwise FST analysis indicated an increased genetic structure between locations in the ChαE7 compared to the concatenated mitochondrial genes. Discriminant analysis of principal components (DAPC) and spatial analysis of molecular variance (SAMOVA) indicated different degrees of genetic structure for all markers, in agreement with the AMOVA results, but with low correlation to geographic data. The NWS fly is considered a panmitic species based on mitochondrial data, while it is structured into three groups considering the ChαE7 gene. A negative association between the two mutations related to organophosphate resistance and Fay & Wu’s H significant negative values for the exons, suggest that these mutations evolved under positive selection.

Contemporary evolution of resistance at the major insecticide target site gene Ace-1 by mutation and copy number variation in the malaria mosquito Anopheles gambiae

Functionally constrained genes are ideal insecticide targets because disruption is often fatal, and resistance mutations are typically costly. Synaptic acetylcholinesterase (AChE) is an essential neurotransmission enzyme targeted by insecticides used increasingly in malaria control. In Anopheles and Culex mosquitoes, a glycine–serine substitution at codon 119 of the Ace-1 gene confers both resistance and fitness costs, especially for 119S/S homozygotes. G119S in Anopheles gambiae from Accra (Ghana) is strongly associated with resistance, and, despite expectations of cost, resistant 119S alleles are increasing significantly in frequency. Sequencing of Accra females detected only a single Ace-1 119S haplotype, whereas 119G diversity was high overall but very low at non-synonymous sites, evidence of strong purifying selection driven by functional constraint. Flanking microsatellites showed reduced diversity, elevated linkage disequilibrium and high differentiation of 119S, relative to 119G homozygotes across up to two megabases of the genome. Yet these signals of selection were inconsistent and sometimes weak tens of kilobases from Ace-1. This unexpected finding is attributable to apparently ubiquitous amplification of 119S alleles as part of a large copy number variant (CNV) far exceeding the size of the Ace-1 gene, whereas 119G alleles were unduplicated. Ace-1 CNV was detectable in archived samples collected when the 119S allele was rare in Ghana. Multicopy amplification of resistant alleles has not been observed previously and is likely to underpin the recent increase in 119S frequency. The large CNV compromised localization of the strong selective sweep around Ace-1, emphasizing the need to integrate CNV analysis into genome scans for selection.

Behavioral cost & overdominance in Anopheles gambiae

In response to the widespread use of control strategies such as Insecticide Treated Nets (ITN), Anopheles mosquitoes have evolved various resistance mechanisms. Kdr is a mutation that provides physiological resistance to the pyrethroid insecticides family (PYR). In the present study, we investigated the effect of the Kdr mutation on the ability of female An. gambiae to locate and penetrate a 1cm-diameter hole in a piece of netting, either treated with insecticide or untreated, to reach a bait in a wind tunnel. Kdr homozygous, PYR-resistant mosquitoes were the least efficient at penetrating an untreated damaged net, with about 51% [39-63] success rate compared to 80% [70-90] and 78% [65-91] for homozygous susceptible and heterozygous respectively. This reduced efficiency, likely due to reduced host-seeking activity, as revealed by mosquito video-tracking, is evidence of a recessive behavioral cost of the mutation. Kdr heterozygous mosquitoes were the most efficient at penetrating nets treated with PYR insecticide, thus providing evidence for overdominance, the rarely-described case of heterozygote advantage conveyed by a single locus. The study also highlights the remarkable capacity of female mosquitoes, whether PYR-resistant or not, to locate holes in bed-nets.

Evolution of insecticide resistance diagnostics in malaria vectors

Malaria control is reliant upon effective, programmatic-scale, anti-vector interventions. The widespread distribution of pyrethroid-treated bednets in sub-Saharan Africa has been a driver of morbidity and mortality reductions over the last decade. Unfortunately resistance to insecticides, and to pyrethroids in particular, is increasingly common in Anopheles malaria vectors, and is a major threat to continued control and future elimination. Here we argue that current methods to diagnose resistance often have limited utility and should be augmented or even partially replaced by wider application of DNA markers.

The recent escalation in strength of pyrethroid resistance in Anopheles coluzzi in West Africa is linked to increased expression of multiple gene families

Background

Since 2011, the level of pyrethroid resistance in the major malaria mosquito, Anopheles coluzzi, has increased to such an extent in Burkina Faso that none of the long lasting insecticide treated nets (LLINs) currently in use throughout the country kill the local mosquito vectors. We investigated whether this observed increase was associated with transcriptional changes in field-caught Anopheles coluzzi using two independent whole-genome microarray studies, performed in 2011 and 2012.

Results

Mosquitoes were collected from south-west Burkina Faso in 2011 and 2012 and insecticide exposed or non-exposed insects were compared to laboratory susceptible colonies using whole-genome microarrays. Using a stringent filtering process we identified 136 genes, including the well-studied detoxification enzymes (p450 monoxygenases and esterases) and non-detoxification genes (e.g. cell transporters and cuticular components), associated with pyrethroid resistance, whose basal expression level increased during the timeframe of the study. A subset of these were validated by qPCR using samples from two study sites, collected over 3 years and marked increases in expression were observed each year. We hypothesise that these genes are contributing to this rapidly increasing resistance phenotype in An. coluzzi. A comprehensive analysis of the knockdown resistance (kdr) mutations (L1014S, L1014F and N1575Y) revealed that the majority of the resistance phenotype is not explained by target-site modifications.

Conclusions

Our data indicate that the recent and rapid increase in pyrethroid resistance observed in south-west Burkina Faso is associated with gene expression profiles described here. Over a third of these candidates are also overexpressed in multiple pyrethroid resistant populations of An. coluzzi from neighbouring Côte d’Ivoire. This suite of molecular markers can be used to track the spread of the extreme pyrethroid resistance phenotype that is sweeping through West Africa and to determine the functional basis of this trait.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1342-6) contains supplementary material, which is available to authorized users.

Adaptive introgression in an African malaria mosquito coincident with the increased usage of insecticide-treated bed nets

Animal species adapt to changes in their environment, including man-made changes such as the introduction of insecticides, through selection for advantageous genes already present in populations or newly arisen through mutation. A possible alternative mechanism is the acquisition of adaptive genes from related species via a process known as adaptive introgression. Differing levels of insecticide resistance between two African malaria vectors, Anopheles coluzzii and Anopheles gambiae, have been attributed to assortative mating between the two species. In a previous study, we reported two bouts of hybridization observed in the town of Selinkenyi, Mali in 2002 and 2006. These hybridization events did not appear to be directly associated with insecticide-resistance genes. We demonstrate that during a brief breakdown in assortative mating in 2006, A. coluzzii inherited the entire A. gambiae-associated 2L divergence island, which includes a suite of insecticide-resistance alleles. In this case, introgression was coincident with the start of a major insecticide-treated bed net distribution campaign in Mali. This suggests that insecticide exposure altered the fitness landscape, favoring the survival of A. coluzzii/A. gambiae hybrids, and provided selection pressure that swept the 2L divergence island through A. coluzzii populations in Mali. We propose that the work described herein presents a unique description of the temporal dynamics of adaptive introgression in an animal species and represents a mechanism for the rapid evolution of insecticide resistance in this important vector of human malaria in Africa.

Long-term trends in Anopheles gambiae insecticide resistance in Côte d'Ivoire

BACKGROUND

Malaria control is heavily dependent on the use of insecticides that target adult mosquito vectors via insecticide treated nets (ITNs) or indoor residual spraying (IRS). Four classes of insecticide are approved for IRS but only pyrethroids are available for ITNs. The rapid rise in insecticide resistance in African malaria vectors has raised alarms about the sustainability of existing malaria control activities. This problem might be particularly acute in Côte d'Ivoire where resistance to all four insecticide classes has recently been recorded. Here we investigate temporal trends in insecticide resistance across the ecological zones of Côte d'Ivoire to determine whether apparent pan-African patterns of increasing resistance are detectable and consistent across insecticides and areas.

METHODS

We combined data on insecticide resistance from a literature review, and bioassays conducted on field-caught Anopheles gambiae mosquitoes for the four WHO-approved insecticide classes for ITN/IRS. The data were then mapped using Geographical Information Systems (GIS) and the IR mapper tool to provide spatial and temporal distribution data on insecticide resistance in An. gambiae sensu lato from Côte d'Ivoire between 1993 and 2014.

RESULTS

Bioassay mortality decreased over time for all insecticide classes, though with significant spatiotemporal variation, such that stronger declines were observed in the southern ecological zone for DDT and pyrethroids than in the central zone, but with an apparently opposite effect for the carbamate and organophosphate. Variation in relative abundance of the molecular forms, coupled with dramatic increase in kdr 1014F frequency in M forms (An. coluzzii) seems likely to be a contributory factor to these patterns. Although records of resistance across insecticide classes have become more common, the number of classes tested in studies has also increased, precluding a conclusion that multiple resistance has also increased.

CONCLUSION

Our analyses attempted synthesis of 22 years of bioassay data from Côte d'Ivoire, and despite a number of caveats and potentially confounding variables, suggest significant but spatially-variable temporal trends in insecticide resistance. In the light of such spatio-temporal dynamics, regular, systematic and spatially-expanded monitoring is warranted to provide accurate information on insecticide resistance for control programme management.

Mosquito genomics. Highly evolvable malaria vectors: the genomes of 16 Anopheles mosquitoes

Variation in vectorial capacity for human malaria among Anopheles mosquito species is determined by many factors, including behavior, immunity, and life history. To investigate the genomic basis of vectorial capacity and explore new avenues for vector control, we sequenced the genomes of 16 anopheline mosquito species from diverse locations spanning ~100 million years of evolution. Comparative analyses show faster rates of gene gain and loss, elevated gene shuffling on the X chromosome, and more intron losses, relative to Drosophila. Some determinants of vectorial capacity, such as chemosensory genes, do not show elevated turnover but instead diversify through protein-sequence changes. This dynamism of anopheline genes and genomes may contribute to their flexible capacity to take advantage of new ecological niches, including adapting to humans as primary hosts.

Islands and stepping-stones: comparative population structure of Anopheles gambiae sensu stricto and Anopheles arabiensis in Tanzania and implications for the spread of insecticide resistance

Population genetic structures of the two major malaria vectors Anopheles gambiae s.s. and An. arabiensis, differ markedly across Sub-Saharan Africa, which could reflect differences in historical demographies or in contemporary gene flow. Elucidation of the degree and cause of population structure is important for predicting the spread of genetic traits such as insecticide resistance genes or artificially engineered genes. Here the population genetics of An. gambiae s.s. and An. arabiensis in the central, eastern and island regions of Tanzania were compared. Microsatellite markers were screened in 33 collections of female An. gambiae s.l., originating from 22 geographical locations, four of which were sampled in two or three years between 2008 and 2010. An. gambiae were sampled from six sites, An. arabiensis from 14 sites, and both species from two sites, with an additional colonised insectary sample of each species. Frequencies of the knock-down resistance (kdr) alleles 1014S and 1014F were also determined. An. gambiae exhibited relatively high genetic differentiation (average pairwise FST = 0.131), significant even between nearby samples, but without clear geographical patterning. In contrast, An. arabiensis exhibited limited differentiation (average FST = 0.015), but strong isolation-by-distance (Mantel test r = 0.46, p = 0.0008). Most time-series samples of An. arabiensis were homogeneous, suggesting general temporal stability of the genetic structure. An. gambiae populations from Dar es Salaam and Bagamoyo were found to have high frequencies of kdr 1014S (around 70%), with almost 50% homozygote but was at much lower frequency on Unguja Island, with no. An. gambiae population genetic differentiation was consistent with an island model of genetic structuring with highly restricted gene flow, contrary to An. arabiensis which was consistent with a stepping-stone model of extensive, but geographically-restricted gene flow.