Human Interventions: Driving Forces of Mosquito Evolution

Behavioral resistance to insecticides: current understanding, challenges, and future directions

Identifying and understanding behavioral resistance to insecticides is vital for maintaining global food security, public health, and ecological balance. Behavioral resistance has been documented to occur in a multitude of insect taxa dating back to the 1940s, but has not received significant research attention due primarily to the complexities of studying insect behavior and a lack of any clear definition of behavioral resistance. In recent years, a systematic effort to investigate the mechanism(s) of behavioral resistance in pest taxa (e.g. the German cockroach and the house fly) has been undertaken. Here, we practically define behavioral resistance, describe the efforts taken by research groups to elucidate resistance mechanisms, and provide insight on designing appropriate bioassays for investigating behavioral resistance mechanisms in the future.

Lethal and sublethal impacts of membrane-fed ivermectin are concentration dependent in Anopheles coluzzii

BACKGROUND

Ivermectin is a well-tolerated anthelminthic drug with wide clinical and veterinary applications. It also has lethal and sublethal effects on mosquitoes. Mass drug administration with ivermectin has therefore been suggested as an innovative vector control tool in efforts to curb emerging insecticide resistance and reduce residual malaria transition. To support assessments of the feasibility and efficacy of current and future formulations of ivermectin for vector control, we sought to establish the relationship between ivermectin concentration and its lethal and sublethal impacts in a primary malaria vector.

METHODS

The in vitro effects of ivermectin on daily mortality and fecundity, measured by egg production, were assessed up to 14 days post-blood feed in a laboratory colony of Anopheles coluzzii. Mosquitoes were fed ivermectin in blood meals delivered by membrane feeding at one of six concentrations: 0 ng/ml (control), 10 ng/ml, 15 ng/ml, 25 ng/ml, 50 ng/ml, 75 ng/ml, and 100 ng/ml.

RESULTS

Ivermectin had a significant effect on mosquito survival in a concentration-dependent manner. The LC50 at 7 days was 19.7 ng/ml. The time to median mortality at ≥ 50 ng/ml was ≤ 4 days, compared to 9.6 days for control, and 6.3-7.6 days for ivermectin concentrations between 10 and 25 ng/ml. Fecundity was also affected; no oviposition was observed in surviving females from the two highest concentration treatment groups. While females exposed to 10 to 50 ng/ml of ivermectin did oviposit, significantly fewer did so in the 50 ng/ml treatment group compared to the control, and they also produced significantly fewer eggs.

CONCLUSIONS

Our results showed ivermectin reduced mosquito survival in a concentration-dependent manner and at ≥ 50 ng/ml significantly reduced fecundity in An. coluzzii. Results indicate that levels of ivermectin found in human blood following ingestion of a single 150-200 μg/kg dose would be sufficient to achieve 50% mortality across 7 days; however, fecundity in survivors is unlikely to be affected. At higher doses, a substantial impact on both survival and fecundity is likely. Treating human populations with ivermectin could be used as a supplementary malaria vector control method to kill mosquito populations and supress their reproduction; however strategies to safely maintain mosquitocidal blood levels of ivermectin against all Anopheles species require development.

Genetic and behavioral differences between above and below ground Culex pipiens bioforms

Efficiency of mosquito-borne disease transmission is dependent upon both the preference and fidelity of mosquitoes as they seek the blood of vertebrate hosts. While mosquitoes select their blood hosts through multi-modal integration of sensory cues, host-seeking is primarily an odor-guided behavior. Differences in mosquito responses to hosts and their odors have been demonstrated to have a genetic component, but the underlying genomic architecture of these responses has yet to be fully resolved. Here, we provide the first characterization of the genomic architecture of host preference in the polymorphic mosquito species, Culex pipiens. The species exists as two morphologically identical bioforms, each with distinct avian and mammalian host preferences. Cx. pipiens females with empirically measured host responses were prepared into reduced representation DNA libraries and sequenced to identify genomic regions associated with host preference. Multiple genomic regions associated with host preference were identified on all 3 Culex chromosomes, and these genomic regions contained clusters of chemosensory genes, as expected based on work in Anopheles gambiae complex mosquitoes and in Aedes aegypti. One odorant receptor and one odorant binding protein gene showed one-to-one orthologous relationships to differentially expressed genes in A. gambiae complex members with divergent host preferences. Overall, our work identifies a distinct set of odorant receptors and odorant binding proteins that may enable Cx. pipiens females to distinguish between their vertebrate blood host species, and opens avenues for future functional studies that could measure the unique contributions of each gene to host preference phenotypes.

Anopheles gambiae larvae's ability to grow and emerge in water containing lethal concentrations of clothianidin, acetamiprid, or imidacloprid is consistent with cross-resistance to neonicotinoids

BACKGROUND

For decades, various agrochemicals have been successfully repurposed for mosquito control. However, preexisting resistance caused in larval and adult populations by unintentional pesticide exposure or other cross-resistance mechanisms poses a challenge to the efficacy of this strategy. A better understanding of larval adaptation to the lethal and sublethal effects of residual pesticides in aquatic habitats would provide vital information for assessing the efficacy of repurposed agrochemicals against mosquitoes.

METHODS

We reared field-collected mosquito larvae in water containing a concentration of agrochemical causing 100% mortality in susceptible mosquitoes after 24 h (lethal concentration). Using this experimental setup, we tested the effect of lethal concentrations of a pyrrole (chlorfenapyr, 0.10 mg/l), a pyrethroid (deltamethrin, 1.5 mg/l), and three neonicotinoids including imidacloprid (0.075 mg/l), acetamiprid (0.15 mg/l), and clothianidin (0.035 mg/l) on mortality rates, growth, and survival in third-instar larvae of the two sibling species Anopheles gambiae and Anopheles coluzzii collected from Yaoundé, Cameroon.

RESULTS

We found that An. gambiae and An. coluzzii larvae were susceptible to chlorfenapyr and were killed within 24 h by a nominal concentration of 0.10 mg/l. Consistent with strong resistance, deltamethrin induced low mortality in both species. Lethal concentrations of acetamiprid, imidacloprid, and clothianidin strongly inhibited survival, growth, and emergence in An. coluzzii larvae. By contrast, depending on the active ingredient and the population tested, 5-60% of immature stages of An. gambiae were able to grow and emerge in water containing a lethal concentration of neonicotinoids, suggesting cross-resistance to this class of insecticides.

CONCLUSIONS

These findings corroborate susceptibility profiles observed in adults and suggest that unintentional pesticide exposure or other cross-resistance processes could contribute to the development of resistance to neonicotinoids in some Anopheles populations.

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

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

Evaluation of the inheritance and dominance of behavioral resistance to imidacloprid in the house fly (Musca domestica L.) (Diptera: Muscidae)

The house fly, Musca domestica, is a cosmopolitan species known for its pestiferous nature and potential to mechanically vector numerous human and animal pathogens. Control of adult house flies often relies on insecticides formulated into food baits. However, due to the overuse of these baits, insecticide resistance has developed to all insecticide classes currently registered for use in the United States. Field populations of house flies have developed resistance to imidacloprid, the most widely used neonicotinoid insecticide for fly control, through both physiological and behavioral resistance mechanisms. In the current study, we conducted a comprehensive analysis of the inheritance and dominance of behavioral resistance to imidacloprid in a lab-selected behaviorally resistant house fly strain. Additionally, we conducted feeding preference assays to assess the feeding responses of genetic cross progeny to imidacloprid. Our results confirmed that behavioral resistance to imidacloprid is inherited as a polygenic trait, though it is inherited differently between male and female flies. We also demonstrated that feeding preference assays can be instrumental in future genetic inheritance studies as they provide direct insight into the behavior of different strains under controlled conditions that reveal, interactions between the organism and the insecticide. The findings of this study carry significant implications for pest management and underscore the need for integrated pest control approaches that consider genetic and ecological factors contributing to resistance.

Identifying potential insecticide resistance markers through genomic-level comparison of Bemisia tabaci (Gennadius) lines

The invasive whitefly (Bemisia tabaci) MED is one of the most economically damaging plant pests. The extensive use of insecticide over decades has led to that the invasive B. tabaci MED has developed resistance to a wide range of insecticide classes, but little is known about the genetic background associated with resistance. To this end, we conducted a comparative genome-wide analysis of single-base nucleotide polymorphisms between MED whitefly lines collected from fields that were recently infested and an insecticide-susceptible MED whitefly line collected in 1976. First, low-coverage genome sequencings were conducted on DNA isolated from individual whiteflies. The sequencing results were evaluated using an available B. tabaci MED genome as a reference. Significant genetic differences were discovered between MED whitefly lines collected from fields that were recently infested and an insecticide-susceptible MED whitefly line based on the principal component analyses. Top GO categories and KEGG pathways that might be involved in insecticide resistance development were identified, and several of them have not been previously associated with resistance. Additionally, we identified several genetic loci with novel variations including Cytochrome P450 monooxygenases (P450s), UDP-glucuronosyltransferases (UGTs), Glutathione S-transferases (GSTs), esterase, carboxyl-esterases (COE), ABC transporters, fatty acyl-CoA reductase, voltage-gated sodium channels, GABA receptor, and cuticle proteins (CPs) that were previously reported to have close associations with pesticide resistance in well-studied insect groups that provide an essential resource for the design of insecticide resistance-linked loci arrays insecticide. Our results was obtained solely on resequencing genome data sets, more pesticide bio-assays combined with omics datasets should be further used to verify the markers identified here.

Population Genomics for Insect Conservation

Insects constitute vital components of ecosystems. There is alarming evidence for global declines in insect species diversity, abundance, and biomass caused by anthropogenic drivers such as habitat degradation or loss, agricultural practices, climate change, and environmental pollution. This raises important concerns about human food security and ecosystem functionality and calls for more research to assess insect population trends and identify threatened species and the causes of declines to inform conservation strategies. Analysis of genetic diversity is a powerful tool to address these goals, but so far animal conservation genetics research has focused strongly on endangered vertebrates, devoting less attention to invertebrates, such as insects, that constitute most biodiversity. Insects' shorter generation times and larger population sizes likely necessitate different analytical methods and management strategies. The availability of high-quality reference genome assemblies enables population genomics to address several key issues. These include precise inference of past demographic fluctuations and recent declines, measurement of genetic load levels, delineation of evolutionarily significant units and cryptic species, and analysis of genetic adaptation to stressors. This enables identification of populations that are particularly vulnerable to future threats, considering their potential to adapt and evolve. We review the application of population genomics to insect conservation and the outlook for averting insect declines.

Concentration Dependent Feeding on Imidacloprid by Behaviorally Resistant House Flies, Musca domestica L. (Diptera: Muscidae)

The house fly (Musca domestica L.) is a cosmopolitan and synanthropic pest fly commonly associated with confined animal facilities, known to mechanically vector numerous disease-causing pathogens. Control of adult house flies often relies on insecticides formulated into insecticidal baits, though many baits have failed due to insecticide resistance. House fly resistance to imidacloprid, the most widely used neonicotinoid insecticide available for fly control, has evolved through physiological and behavioral mechanisms in field populations. Behavioral resistance to imidacloprid was documented in field populations of flies from southern California dairies. Lab colonies of these flies were established and behavioral resistance to imidacloprid was selected over several generations. The current study examined the ability of these lab-selected flies to feed on varying concentrations of imidacloprid formulated in sucrose, and if these flies would demonstrate a feeding preference for different concentrations of imidacloprid when exposed in bioassays. Behaviorally resistant flies preferred to feed on untreated sucrose as opposed to treated sucrose at concentrations greater than 25 µg/g imidacloprid when provided sucrose treated with and without imidacloprid. When provisioned with only sucrose treated with a low and high imidacloprid concentration, flies fed on the low concentrations (≤100 µg/g) imidacloprid but reduced feeding on either treatment when concentrations were >100 µg/g imidacloprid. The current study extends the body of knowledge on house fly behavioral resistance to imidacloprid, which could provide insights into future failures of granular fly baits.

Signs of Urban Evolution? Morpho-Functional Traits Co-variation Along a Nature-Urban Gradient in a Chagas Disease Vector

Environmental change (i.e., urbanization) impacts species in contrasting ways, with some species experiencing benefits given their way of life (i.e., blood-sucking insects). How these species respond to such change is not well understood and for species involved in human diseases, this “how” question is particularly important. Most Triatominae bug species inhabit tropical and subtropical forests where their vertebrate hosts’ temporal abundance depends on climate seasonality. However, in human encroached landscapes, triatomines can benefit from resource stability which may lead to adaptive phenotypic change to track novel hosts. We tested for an association between different landscapes and morpho-functional traits linked to sensory, motion, and feeding functions in Triatoma dimidiata and compared fecundity (i.e., number of eggs) in each landscape as a proxy of fitness. Using geometric and traditional morphometric tools, we predicted a morphological simplification in bugs inhabiting urbanized areas. While wing morphology or proboscis were not influenced by landscape class, the opposite occurred for thorax morphology and number of sensilla. Wing and thorax morphology did not covary under modified landscape scenarios, yet we detected a morpho-functional convergence for thorax size and antennal phenotype in both sexes, with a simplification trend, from nature to urban settings. Given no fecundity differences across landscapes, there is no potential reproductive costs. Moreover, the convergence of thorax size and antennal phenotype suggests differences in flight/locomotion performance and host/environment perception, as a possible adaptive response to relaxed selective pressures of the bug’s native habitat. These results imply that T. dimidiata could be adapting to urbanized areas.

Transposable element variants and their potential adaptive impact in urban populations of the malaria vector Anopheles coluzzii

Anopheles coluzzii is one of the primary vectors of human malaria in sub-Saharan Africa. Recently, it has spread into the main cities of Central Africa threatening vector control programs. The adaptation of An. coluzzii to urban environments partly results from an increased tolerance to organic pollution and insecticides. Some of the molecular mechanisms for ecological adaptation are known, but the role of transposable elements (TEs) in the adaptive processes of this species has not been studied yet. As a first step toward assessing the role of TEs in rapid urban adaptation, we sequenced using long reads six An. coluzzii genomes from natural breeding sites in two major Central Africa cities. We de novo annotated TEs in these genomes and in an additional high-quality An. coluzzii genome, and we identified 64 new TE families. TEs were nonrandomly distributed throughout the genome with significant differences in the number of insertions of several superfamilies across the studied genomes. We identified seven putatively active families with insertions near genes with functions related to vectorial capacity, and several TEs that may provide promoter and transcription factor binding sites to insecticide resistance and immune-related genes. Overall, the analysis of multiple high-quality genomes allowed us to generate the most comprehensive TE annotation in this species to date and identify several TE insertions that could potentially impact both genome architecture and the regulation of functionally relevant genes. These results provide a basis for future studies of the impact of TEs on the biology of An. coluzzii.

Understanding Mosquito Surveillance Data for Analytic Efforts: A Case Study

Mosquito surveillance data can be used for predicting mosquito distribution and dynamics as they relate to human disease. Often these data are collected by independent agencies and aggregated to state and national level portals to characterize broad spatial and temporal dynamics. These larger repositories may also share the data for use in mosquito and/or disease prediction and forecasting models. Assumed, but not always confirmed, is consistency of data across agencies. Subtle differences in reporting may be important for development and the eventual interpretation of predictive models. Using mosquito vector surveillance data from Arizona as a case study, we found differences among agencies in how trapping practices were reported. Inconsistencies in reporting may interfere with quantitative comparisons if the user has only cursory familiarity with mosquito surveillance data. Some inconsistencies can be overcome if they are explicit in the metadata while others may yield biased estimates if they are not changed in how data are recorded. Sharing of metadata and collaboration between modelers and vector control agencies is necessary for improving the quality of the estimations. Efforts to improve sharing, displaying, and comparing vector data from multiple agencies are underway, but existing data must be used with caution.

Selection, Reversion, and Characterization of House Fly (Diptera: Muscidae) Behavioral Resistance to the Insecticide Imidacloprid

Insecticide resistance in pest populations is an increasing problem in both urban and rural settings caused by over-application of insecticides and lack of rotation among chemical classes. The house fly (Musca domestica L.) is a cosmopolitan fly species implicated in the transmission of numerous pathogens, and which can be extremely pestiferous when present in high numbers. The evolution of insecticide resistance has long been documented in house flies, with resistance reported to all major insecticide classes. House fly resistance to imidacloprid, the most widely used neonicotinoid insecticide available for fly control, has been selected for in field populations through both physiological and behavioral resistance mechanisms. In the current study, house flies collected from a southern California dairy were selectively bred for behavioral resistance to imidacloprid, without increasing the physiological resistance profile of the selected flies. Flies were also successfully selected for behavioral susceptibility to imidacloprid. The rapid selection for either behavioral resistance or behavioral susceptibility suggests that inheritable alleles conferring behavioral resistance were already present in the wild-type fly population collected from the dairy site. The methods used for the specific selection of behavioral resistance (or susceptibility) in the fly population will be useful for further studies on the specific mechanisms conferring this resistance. House fly behavioral resistance was further investigated using behavioral observation and feeding preference assays, with resistance determined to be both contact-dependent and specific to the insecticide (imidacloprid) rather than to a non-insecticidal component of a bait matrix as previously documented.

Adaptation to agricultural pesticides may allow mosquitoes to avoid predators and colonize novel ecosystems

Human activities such as the application of agrochemicals may detrimentally disturb natural ecosystems, generating novel selection pressures. Here we examine how pesticides may influence community composition using the aquatic communities within bromeliad phytotelmata, and how adaptive responses to pesticides may influence community-level patterns. We first quantified the composition of macroinvertebrate communities from pesticide-free and pesticide-exposed locations. Complementary manipulative experiments where bromeliads were transplanted between pesticide-free and pesticide-exposed sites were then performed. Finally, pesticide bioassays on the most common predators (Mecistogaster modesta damselflies) and prey (Wyeomyia abebela mosquitoes) assessed a potential evolutionary mechanism that may influence community compositional differences. Our field survey revealed differences in W. abebela and M. modesta abundances between pesticide-free and pesticide-exposed areas. Our transplant experiment suggested compositional differences were not due to physical differences between bromeliads from different locations. Pesticide bioassays revealed that M. modesta from pesticide-free locations had higher innate pesticide tolerances than W. abebela from pesticide-free areas, but M. modesta larvae showed no evidence of adapted resistance as none were found where pesticides were used. Conversely, W. abebela larvae from pesticide-exposed locations had higher pesticide tolerances than individuals from pesticide-free sites, suggesting an adaptive response. This evolved resistance to pesticides may, therefore, allow W. abebela to colonize habitats free of the dominant predator in the system, explaining the higher W. abebela abundances in pesticide-exposed areas than in pesticide-free locations. We suggest that the total effect of novel stressors is driven by interactions between ecological and evolutionary processes.

Modelling the potential of genetic control of malaria mosquitoes at national scale

BACKGROUND

The persistence of malaria in large parts of sub-Saharan Africa has motivated the development of novel tools to complement existing control programmes, including gene-drive technologies to modify mosquito vector populations. Here, we use a stochastic simulation model to explore the potential of using a driving-Y chromosome to suppress vector populations in a 10⁶ km² area of West Africa including all of Burkina Faso.

RESULTS

The consequence of driving-Y introductions is predicted to vary across the landscape, causing elimination of the target species in some regions and suppression in others. We explore how this variation is determined by environmental conditions, mosquito behaviour, and the properties of the gene-drive. Seasonality is particularly important, and we find population elimination is more likely in regions with mild dry seasons whereas suppression is more likely in regions with strong seasonality.

CONCLUSIONS

Despite the spatial heterogeneity, we suggest that repeated introductions of modified mosquitoes over a few years into a small fraction of human settlements may be sufficient to substantially reduce the overall number of mosquitoes across the entire geographic area.

Integrated Mosquito Management: Is Precision Control a Luxury or Necessity?

The versatility of mosquito species that spread emerging arthropod-borne viruses such as Zika has highlighted the urgent need to re-evaluate mosquito-control standards. The prospect of using precise knowledge of the geographic distribution and vector status of local populations to guide targeted interventions has gained renewed attention, but the feasibility and utility of such an approach remain to be investigated. Using the example of mosquito management in the USA, we present ideas for designing, monitoring, and assessing precision vector control tailored to different environmental and epidemiological settings. We emphasize the technical adjustments that could be implemented in mosquito-control districts to enable targeted control while strengthening traditional management.