Transcriptomic meta-signatures identified in Anopheles gambiae populations reveal previously undetected insecticide resistance mechanisms

High pyrethroid resistance is associated with high frequencies of 1014F and 1014S kdr mutations in Anopheles arabiensis (Diptera: Culicidae) from Ouagadougou, Burkina Faso

Malaria remains a major public health threat in Burkina Faso, as in most sub-Saharan Africa countries. Malaria control relies mainly on long-lasting insecticide-treated nets (LLINs) and indoor residual spraying. In Burkina Faso, an escalating of insecticide resistance has been observed over the last decades. This study aimed to investigate insecticide resistance and the underlying mechanisms in Anopheles gambiae complex in Ouagadougou. Anopheles gambiae s.l. larvae were collected from gutters and ponds, in Zogona, Tampouy and Tanghin, 3 localities in Ouagadougou from July to October 2018. The larvae were reared in the laboratory to adults stage and susceptibility profile to pyrethroid, carbamate, and organophosphate insecticides was assessed using WHO tube assays. Mosquito species and mutations linked with insecticide resistance, were identified through PCR. More than 95% of the collected An. gambiae s.l. were An. arabiensis. An. arabiensis displayed high resistance to permethrin and deltamethrin, with mortalities below 30%, but was fully susceptible to bendiocarb, fenitrothion, and malathion. A high-frequency of the pyrethroid resistance-associated kdr mutation 1014F (0.81) was recorded, while the frequency of 1014S mutation (0.18) was lower. However, the carbamate and organophosphate-associated Ace-1 119S mutation was not detected. Localities and breeding site type appear to influence pyrethroid resistance in the An. arabiensis population of Ouagadougou. The high resistance to pyrethroids in An. arabiensis of urban Ouagadougou is underpinned, at least in part by high-frequency kdr mutations. This result supports the switch to next-generation LLINs, in well-established pyrethroid resistance zones of Burkina Faso including Ouagadougout.

Omics Approaches in Understanding Insecticide Resistance in Mosquito Vectors

In recent years, the emergence of insecticide resistance has been a major challenge to global public health. Understanding the molecular mechanisms of this phenomenon in mosquito vectors is paramount for the formulation of effective vector control strategies. This review explores the current knowledge of insecticide resistance mechanisms through omics approaches. Genomic, transcriptomic, proteomic, and metabolomics approaches have proven crucial to understand these resilient vectors. Genomic studies have identified multiple genes associated with insecticide resistance, while transcriptomics has revealed dynamic gene expression patterns in response to insecticide exposure and other environmental stimuli. Proteomics and metabolomics offer insights into protein expression and metabolic pathways involved in detoxification and resistance. Integrating omics data holds immense potential to expand our knowledge on the molecular basis of insecticide resistance in mosquitoes via information obtained from different omics platforms to understand regulatory mechanisms and differential expression of genes and protein, and to identify the transcription factors and novel molecules involved in the detoxification of insecticides. Eventually, these data will help construct predictive models, identify novel strategies, and develop targeted interventions to control vector-borne diseases.

Chemosensory Proteins Protect Nilaparvata lugens from Imidacloprid by Sequestering the Insecticide and Facilitating Metabolic Detoxification

The involvement of chemosensory proteins (CSPs) in binding to insecticides has been implicated. However, our understanding of CSP-mediated insecticide resistance remains limited. Herein, 15 CSP genes were identified and characterized from Nilaparvata lugens. Expression analysis identified six CSPs with overexpression in the imidacloprid-resistant strain, whose involvement in imidacloprid resistance was validated by RNA interference. Among them, four CSPs were successfully expressed using a prokaryotic expression system, and their binding affinities to imidacloprid were confirmed through fluorescence competitive binding assays. Knockdown of them impaired the capacity of N. lugens to metabolize imidacloprid and inhibited the activity of metabolic detoxification pathways, while their overexpression in Escherichia coli enhanced bacterial metabolic efficiency toward imidacloprid. Furthermore, the transcriptional regulation of CSP2 and CSP15 was found to be mediated by AhR/ARNT and CncC/MafK. These findings suggest that the overexpression of CSPs in N. lugens promotes imidacloprid resistance by sequestering the insecticide and enhancing metabolic detoxification.

Targeted genomic surveillance of insecticide resistance in African malaria vectors

The emergence of insecticide resistance is threatening the efforts of malaria control programmes, which rely heavily on a limited arsenal of insecticidal tools, such as insecticide-treated bed nets. Importantly, genomic surveillance of malaria vectors can provide critical, policy-relevant insights into the presence and evolution of insecticide resistance, allowing us to maintain and extend the shelf life of these interventions. Yet the complex genetic architecture of resistance, combined with resource constraints in malaria-endemic settings, have thus far precluded the widespread use of genomics in routine surveillance. Meanwhile, stakeholders in sub-Saharan Africa are moving towards locally driven, decentralised generation of genomic data, underscoring the need for standardised and robust genomics workflows. To address this need, we demonstrate an approach to targeted genomic surveillance in Anopheles gambiae s.l with Illumina sequencing. We target 90 genomic loci in the Anopheles gambiae s.l genome, including 55 resistance-associated mutations and 35 ancestry informative markers. This protocol is coupled with advanced, automated software for accurate and reproducible variant analysis. We are able to elucidate population structure and ancestry in our cohorts and accurately identify most species in the An. gambiae species complex. We report frequencies of variants at insecticide-resistance loci and explore the continued evolution of the pyrethroid target site, the Voltage-gated sodium channel. Applying the platform to a recently established colony of field-caught resistant mosquitoes (Siaya, Kenya), we identified seven independent resistance-associated variants contributing to reduced efficacy of insecticide-treated nets in East Africa. Additionally, we leverage a machine learning algorithm (XGBoost) to demonstrate the possibility of predicting bioassay mortality using genotypes alone. This achieved very high accuracy (73%), demonstrating the potential of targeted genomics to predictively monitor insecticide resistance. Together these tools provide a practical, scalable solution for resistance monitoring while advancing the goal of building local genomic surveillance capacity in sub-Saharan Africa.

Insecticide resistant Anopheles from Ethiopia but not Burkina Faso show a microbiota composition shift upon insecticide exposure

BACKGROUND

Malaria remains a key contributor to mortality and morbidity across Africa, with the highest burden in children under 5. Insecticide-based vector control tools, which target the adult Anopheles mosquitoes, are the most efficacious tool in disease prevention. Due to the widespread use of these interventions, insecticide resistance to the most used classes of insecticides is now pervasive across Africa. Understanding the underlying mechanisms contributing to this phenotype is necessary to both track the spread of resistance and to design new tools to overcome it.

METHODS

Here, we compare the microbiota composition of insecticide-resistant populations of Anopheles gambiae, An. coluzzii and An. arabiensis from Burkina Faso, and in the latter case additionally from Ethiopia, to insecticide-susceptible populations.

RESULTS

We show that the microbiota composition between insecticide-resistant and -susceptible populations does not differ in Burkina Faso. This result is supported by data from laboratory colonies originating in Burkina Faso across two countries. In contrast, An. arabiensis from Ethiopia demonstrates clear differences in microbiota composition in those dying from and those surviving insecticide exposure. To further understand resistance in this An. arabiensis population, we performed RNAseq and saw differential expression of detoxification genes associated with insecticide resistance and changes in respiration, metabolism and synapse-related ion channels.

CONCLUSIONS

Our results indicate that, in addition to changes in the transcriptome, microbiota can contribute to insecticide resistance in certain settings.

Mechanistic exploration of odorant binding protein-mediated chlorpyrifos resistance in Nilaparvata lugens: Insights from insecticide sequestration and transcriptional regulation

The effectiveness and sustainable application of insecticides are severely threatened by the rapid evolution of resistance in agricultural pests. Recent research indicates that odorant binding proteins (OBPs) may be involved in facilitating insecticide resistance, while the specific mechanisms remain poorly understood. Herein, 11 OBPs were identified from Nilaparvata lugens. Among them, OBP5 exhibited high and specific expression in the head, and showed constitutive overexpression in the chlorpyrifos-resistant strain. Knockdown of OBP5 notably restored susceptibility to chlorpyrifos in N. lugens, while overexpression of OBP5 in Escherichia coli significantly enhanced bacterial tolerance to chlorpyrifos. Fluorescence competitive binding assay confirmed the strong binding affinities of OBP5 to chlorpyrifos and its active metabolite chlorpyrifos-oxon. Molecular docking studies proposed a critical interacting amino acid (Lys147) in the binding site, which was further validated by comparative binding studies between wildtype OBP5 and the mutated protein OBP5K147A. Furthermore, Lim1β that also presented overexpression pattern in the resistant strain, was found to regulate expression of OBP5 through a dual-luciferase reporter assay. Our findings demonstrate that the overexpression of OBP5 contributes to chlorpyrifos resistance by binding and sequestering the insecticides, shedding light on the sequestration resistance mechanism conferred by OBPs and offering potential targets for resistance management.

Advancing Energy Sustainability Through Solar-to-Fuel Technologies: From Materials to Devices and Systems

To achieve carbon neutrality and sustainable development, innovative solar-to-fuel systems have been designed through the integration of solar energy harvesting and electrochemical devices. Over the last decade, there have been notable advancements in enhancing the efficiency and durability of these solar-to-fuel systems. Despite the advancements, there remains significant potential for further improvements in the performance of systems. Enhancements can be achieved by optimizing electrochemical catalysts, advancing the manufacturing technologies of photovoltaics and electrochemical cells, and refining the overall design of these systems. In the realm of catalyst optimization, the effectiveness of materials can be significantly improved through active site engineering and strategic use of functional groups. Similarly, the performance of electrochemical devices can be enhanced by incorporating specific additives into electrolytes and optimizing gas diffusion electrodes. Improvements in solar harvesting devices are achievable through efficient passivant and self-assembled monolayers, which enhance the overall quality and efficiency of these systems. Additionally, optimizing the energy conversion efficiency involves the strategic use of DC converters, photoelectrodes, and redox media. This review aims to provide a comprehensive overview of the advancements in solar-powered electrochemical energy conversion systems, laying a solid foundation for future research and development in the field of energy sustainability.

Bio-efficacy of field aged novel class of long-lasting insecticidal nets, against pyrethroid-resistant malaria vectors in Tanzania: A series of experimental hut trials

New classes of long-lasting insecticidal nets (LLINs), have been recommended by the World Health Organization (WHO) to control malaria vectors resistant to pyrethroid insecticides. This study was nested in a large-scale cluster-randomized controlled trial conducted (cRCT) in Tanzania. A series of experimental hut trials (EHTs) aimed to evaluate the bio-efficacy of trial LLINs on mosquito indicators most pertinent to malaria transmission over 3 years of use in the community in order to better understand the outcomes of the cRCT. The following field-collected LLINs were assessed: 1) Olyset Plus (combining piperonyl butoxide synergist and permethrin), 2) Interceptor G2 (chlorfenapyr and alpha-cypermethrin), 3) Royal Guard (pyriproxyfen and alpha-cypermethrin), 4) Interceptor (alpha-cypermethrin only) conducted in parallel with 5) a new Interceptor, and 6) an untreated net. Thirty nets of each type were withdrawn from the community at 12, 24, and 36 months after distribution and used for the EHTs. Pre-specified outcomes were 72-hour mortality for Interceptor G2, 24-hour mortality for Olyset Plus, and fertility based on egg development stage for Royal Guard. Overall, Interceptor G2 LLINs induced higher 72-hour mortality compared to standard LLINs of the same age up to12 months (44% vs 21%, OR: 3.5, 95% CI: 1.9-6.6, p-value < 0.001), and 24-hour mortality was only significantly higher in Olyset Plus when new (OR: 13.6, 95%CI: 4.4-41.3, p-value < 0.001) compared to standard LLINs but not at 12 months (17% vs 13%; OR: 2.1, 95% CI: 1.0-4.3; p-value = 0.112). A small, non-significant effect of pyriproxyfen on Anopheles fertility was observed for Royal Guard up to 12 months (75% vs 98%, OR: 1.1, 95% CI: 0.0-24.9, p-value = 0.951). There was no evidence of a difference in the main outcomes for any of the new classes of LLINs at 24 and 36 months compared to standard LLINs. Interceptor G2 LLINs showed superior bio-efficacy compared to standard LLINs only up to 12 months, and the effect of Olyset Plus was observed when new for all species and 12 months for An. gambiae s.l. only. The pyriproxyfen component of Royal Guard had a short and limited effect on fertility. The decrease in effectiveness of Olyset Plus and Royal Guard LLINs in the EHTs aligns with findings from the cRCT, whereas efficacy of Interceptor G2 lasted for a longer period in the cRCT compared to the EHT. Further investigations are needed to understand the complete scope of chlorfenapyr mode of action. Additional EHT in various contexts will help confirm the residual efficacy of the dual active ingredient LLINs and support the development of longer-lasting nets.

A cub and sushi domain-containing protein with esterase-like activity confers insecticide resistance in the Indian malaria vector Anopheles stephensi

Chemical insecticides (organophosphates and pyrethroids) in the form of IRS (Indoor Residual Sprays) and LLINs (Long Lasting Insecticidal Nets) are the cornerstone for vector control, globally. However, their incessant use has resulted in widespread development of resistance in mosquito vectors, warranting continuous monitoring and investigation of the underlying mechanisms of resistance. Here, we identified a previously uncharacterized- Cub and Sushi Domain containing Insecticide Resistance (CSDIR) protein and generated evidence for its role in mediating insecticide resistance in the Anopheles stephensi. A strong binding affinity of the CSDIR protein towards different classes of insecticide molecules-malathion (KD 6.43 μM) and deltamethrin (KD 46.7 μM) were demonstrated using MD simulation studies and Surface Plasmon Resonance (SPR) experiments. Further, the recombinant CSDIR913-1190 protein exhibited potent esterase-like activity (α-naphthyl acetate (α-NA)- 1.356 ± 0.262 mM/min/mg and β-naphthyl acetate (β -NA)- 1.777 ± 0.220 mM/min/mg). Interestingly, dsRNA-mediated gene silencing of the CSDIR transcripts caused >60% mortality in resistant An. stephensi upon 1-h exposure to deltamethrin and malathion insecticides, compared to the control group. A significant reduction in the esterase-like activity was also observed against α-NA (p = 0.004) and β-NA (p = 0.025) in CSDIR silenced mosquitoes compared to the control group. Using computational analysis and experimental data, our results provided significant evidence of the involvement of the CSDIR protein in mediating insecticide resistance in Anopheles mosquitoes. Thereby making the CSDIR protein, a novel candidate for exploration of novel insecticide molecules. These data would also be helpful in further understanding the development of metabolic resistance by the Anopheles vector.

Mechanisms of transcriptional regulation in Anopheles gambiae revealed by allele-specific expression

Malaria control relies on insecticides targeting the mosquito vector, but this is increasingly compromised by insecticide resistance, which can be achieved by elevated expression of detoxifying enzymes that metabolize the insecticide. In diploid organisms, gene expression is regulated both in cis, by regulatory sequences on the same chromosome, and by trans acting factors, affecting both alleles equally. Differing levels of transcription can be caused by mutations in cis-regulatory modules (CRM), but few of these have been identified in mosquitoes. We crossed bendiocarb-resistant and susceptible Anopheles gambiae strains to identify cis-regulated genes that might be responsible for the resistant phenotype using RNAseq, and CRM sequences controlling gene expression in insecticide resistance relevant tissues were predicted using machine learning. We found 115 genes showing allele-specific expression (ASE) in hybrids of insecticide susceptible and resistant strains, suggesting cis-regulation is an important mechanism of gene expression regulation in A. gambiae. The genes showing ASE included a higher proportion of Anopheles-specific genes on average younger than genes with balanced allelic expression.

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

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

Indoor residual spraying of experimental huts in Cameroon highlights the potential of Fludora® Fusion to control wild pyrethroid-resistant malaria vectors

Elevated resistance to pyrethroids in major malaria vectors has led to the introduction of novel insecticides including neonicotinoids. There is a fear that efficacy of these new insecticides could be impacted by cross-resistance mechanisms from metabolic resistance to pyrethroids. In this study, after evaluating the resistance to deltamethrin, clothianidin and mixture of clothianidin + deltamethrin in the lab using CDC bottle assays, the efficacy of the new IRS formulation Fludora® Fusion was tested in comparison to clothianidin and deltamethrin applied alone using experimental hut trials against wild free-flying pyrethroid-resistant Anopheles funestus from Elende and field An. gambiae collected from Nkolondom reared in the lab and released in the huts. Additionally, cone tests on the treated walls were performed each month for a period of twelve months to evaluate the residual efficacy of the sprayed products. Furthermore, the L1014F-kdr target-site mutation and the L119F-GSTe2 mediated metabolic resistance to pyrethroids were genotyped on a subset of mosquitoes from the EHT to assess the potential cross-resistance. All Anopheles species tested were fully susceptible to clothianidin and clothianidin + deltamethrin mixture in CDC bottle assay while resistance was noted to deltamethrin. Accordingly, Fludora® Fusion (62.83% vs 42.42%) and clothianidin (64.42% vs 42.42%) induced significantly higher mortality rates in EHT than deltamethrin (42.42%) against free flying An. funestus from Elende in month 1 (M1) and no significant difference in mortality was observed between the first (M1) and sixth (M6) months of the evaluation (P > 0.05). However, lower mortality rates were recorded against An. gambiae s.s from Nkolondom (mortality rates 50%, 45.56% and 26.68%). In-situ cone test on the wall showed a high residual efficacy of Fludora® Fusion and clothianidin on the susceptible strain KISUMU (> 12 months) and moderately on the highly pyrethroid-resistant An. gambiae strain from Nkolondom (6 months). Interestingly, no association was observed between the L119F-GSTe2 mutation and the ability of mosquitoes to survive exposure to Fludora® Fusion, whereas a trend was observed with the L1014F-kdr mutation. This study highlights that Fludora® Fusion, through its clothianidin component, has good potential of controlling pyrethroid-resistant mosquitoes with prolonged residual efficacy. This could be therefore an appropriate tool for vector control in several malaria endemic regions.

Parallel Evolution in Mosquito Vectors-A Duplicated Esterase Locus is Associated With Resistance to Pirimiphos-methyl in Anopheles gambiae

The primary control methods for the African malaria mosquito, Anopheles gambiae, are based on insecticidal interventions. Emerging resistance to these compounds is therefore of major concern to malaria control programs. The organophosphate (OP), pirimiphos-methyl, is a relatively new chemical in the vector control armory but is now widely used in indoor-residual spray campaigns. While generally effective, phenotypic resistance has developed in some areas in malaria vectors. Here, we used a population genomic approach to identify novel mechanisms of resistance to pirimiphos-methyl in A. gambiae s.l mosquitoes. In multiple populations, we found large and repeated signals of selection at a locus containing a cluster of detoxification enzymes, some of whose orthologs are known to confer resistance to OPs in Culex pipiens. Close examination revealed a pair of alpha-esterases, Coeae1f and Coeae2f, and a complex and diverse pattern of haplotypes under selection in A. gambiae, A. coluzzii and A. arabiensis. As in C. pipiens, copy number variants have arisen at this locus. We used diplotype clustering to examine whether these signals arise from parallel evolution or adaptive introgression. Using whole-genome sequenced phenotyped samples, we found that in West Africa, a copy number variant in A. gambiae is associated with resistance to pirimiphos-methyl. Overall, we demonstrate a striking example of contemporary parallel evolution which has important implications for malaria control programs.

Genetic surveillance of insecticide resistance in African Anopheles populations to inform malaria vector control

Insecticide resistance in malaria vector populations poses a major threat to malaria control, which relies largely on insecticidal interventions. Contemporary vector-control strategies focus on combatting resistance using multiple insecticides with differing modes of action within the mosquito. However, diverse genetic resistance mechanisms are present in vector populations, and continue to evolve. Knowledge of the spatial distribution of these genetic mechanisms, and how they impact the efficacy of different insecticidal products, is critical to inform intervention deployment decisions. We developed a catalogue of genetic-resistance mechanisms in African malaria vectors that could guide molecular surveillance. We highlight situations where intervention deployment has led to resistance evolution and spread, and identify challenges in understanding and mitigating the epidemiological impacts of resistance.

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

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

Distribution and insecticide resistance profile of the major malaria vector Anopheles funestus group across the African continent

There has been significant progress in malaria control in the last 2 decades, with a decline in mortality and morbidity. However, these gains are jeopardised by insecticide resistance, which negatively impacts the core interventions, such as insecticide-treated nets (ITN) and indoor residual spraying (IRS). While most malaria control and research efforts are still focused on Anopheles gambiae complex mosquitoes, Anopheles funestus remains an important vector in many countries and, in some cases, contributes to most of the local transmission. As countries move towards malaria elimination, it is important to ensure that all dominant vector species, including An. funestus, an important vector in some countries, are targeted. The objective of this review is to compile and discuss information related to A. funestus populations' resistance to insecticides and the mechanisms involved across Africa, emphasising the sibling species and their resistance profiles in relation to malaria elimination goals. Data on insecticide resistance in An. funestus malaria vectors in Africa were extracted from published studies. Online bibliographic databases, including Google Scholar and PubMed, were used to search for relevant studies. Articles published between 2000 and May 2023 reporting resistance of An. funestus to insecticides and associated mechanisms were included. Those reporting only bionomics were excluded. Spatial variation in species distribution and resistance to insecticides was recorded from 174 articles that met the selection criteria. It was found that An. funestus was increasingly resistant to the four classes of insecticides recommended by the World Health Organisation for malaria vector control; however, this varied by country. Insecticide resistance appears to reduce the effectiveness of vector control methods, particularly IRS and ITN. Biochemical resistance due to detoxification enzymes (P450s and glutathione-S-transferases [GSTs]) in An. funestus was widely recorded. However, An. funestus in Africa remains susceptible to other insecticide classes, such as organophosphates and neonicotinoids. This review highlights the increasing insecticide resistance of An. funestus mosquitoes, which are important malaria vectors in Africa, posing a significant challenge to malaria control efforts. While An. funestus has shown resistance to the recommended insecticide classes, notably pyrethroids and, in some cases, organochlorides and carbamates, it remains susceptible to other classes of insecticides such as organophosphates and neonicotinoids, providing potential alternative options for vector control strategies. The study underscores the need for targeted interventions that consider the population structure and geographical distribution of An. funestus, including its sibling species and their insecticide resistance profiles, to effectively achieve malaria elimination goals.

New and emerging mechanisms of insecticide resistance

The continuous use of insecticides over the last eight decades has led to the development of resistance to these insecticides. Research in the last few decades showed that the mechanisms underlying resistance are diverse but can generally be classified under several modes of resistance such as target-site resistance, metabolic resistance, and penetration resistance. In this review, we highlight new discoveries in insecticide resistance research made over the past few years, including an emerging new mode of resistance, sequestration resistance, where the overexpression of olfactory proteins binds and sequesters insecticides in resistant strains, as well as recent research on how posttranscriptional regulation can impact resistance. Future research will determine the generality of these emerging mechanisms across insect species.

A perspective on the impacts of microplastics on mosquito biology and their vectorial capacity

Microplastics (plastic particles <5 mm) permeate aquatic and terrestrial ecosystems and constitute a hazard to animal life. Although much research has been conducted on the effects of microplastics on marine and benthic organisms, less consideration has been given to insects, especially those adapted to urban environments. Here, we provide a perspective on the potential consequences of exposure to microplastics within typical larval habitat on mosquito biology. Mosquitoes represent an ideal organism in which to explore the biological effects of microplastics on terrestrial insects, not least because of their importance as an infectious disease vector. Drawing on evidence from other organisms and knowledge of the mosquito life cycle, we summarise some of the more plausible impacts of microplastics including physiological, ecotoxicological and immunological responses. We conclude that although there remains little experimental evidence demonstrating any adverse effect on mosquito biology or pathogen transmission, significant knowledge gaps remain, and there is now a need to quantify the effects that microplastic pollution could have on such an important disease vector.

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

Background

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

Methods

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

Results

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

Conclusions

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

Transcriptomic analysis of Anopheles gambiae from Benin reveals overexpression of salivary and cuticular proteins associated with cross-resistance to pyrethroids and organophosphates

BACKGROUND

Insecticide resistance (IR) is one of the major threats to malaria vector control programs in endemic countries. However, the mechanisms underlying IR are poorly understood. Thus, investigating gene expression patterns related to IR can offer important insights into the molecular basis of IR in mosquitoes. In this study, RNA-Seq was used to characterize gene expression in Anopheles gambiae surviving exposure to pyrethroids (deltamethrin, alphacypermethrin) and an organophosphate (pirimiphos-methyl).

RESULTS

Larvae of An. gambiae s.s. collected from Bassila and Djougou in Benin were reared to adulthood and phenotyped for IR using a modified CDC intensity bottle bioassay. The results showed that mosquitoes from Djougou were more resistant to pyrethroids (5X deltamethrin: 51.7% mortality; 2X alphacypermethrin: 47.4%) than Bassila (1X deltamethrin: 70.7%; 1X alphacypermethrin: 77.7%), while the latter were more resistant to pirimiphos-methyl (1.5X: 48.3% in Bassila and 1X: 21.5% in Djougou). RNA-seq was then conducted on resistant mosquitoes, non-exposed mosquitoes from the same locations and the laboratory-susceptible An. gambiae s.s. Kisumu strain. The results showed overexpression of detoxification genes, including cytochrome P450s (CYP12F2, CYP12F3, CYP4H15, CYP4H17, CYP6Z3, CYP9K1, CYP4G16, and CYP4D17), carboxylesterase genes (COEJHE5E, COE22933) and glutathione S-transferases (GSTE2 and GSTMS3) in all three resistant mosquito groups analyzed. Genes encoding cuticular proteins (CPR130, CPR10, CPR15, CPR16, CPR127, CPAP3-C, CPAP3-B, and CPR76) were also overexpressed in all the resistant groups, indicating their potential role in cross resistance in An. gambiae. Salivary gland protein genes related to 'salivary cysteine-rich peptide' and 'salivary secreted mucin 3' were also over-expressed and shared across all resistant groups.

CONCLUSION

Our results suggest that in addition to metabolic enzymes, cuticular and salivary gland proteins could play an important role in cross-resistance to multiple classes of insecticides in Benin. These genes warrant further investigation to validate their functional role in An. gambiae resistance to insecticides.

Genomic Profiling of Insecticide Resistance in Malaria Vectors: Insights into Molecular Mechanisms

Malaria control faces challenges from widespread insecticide resistance in major Anopheles species. This study, employing a cross-species approach, integrates RNA-Sequencing, whole-genome sequencing, and microarray data to elucidate drivers of insecticide resistance in Anopheles gambiae complex and An. funestus. Findings show an inverse relationship between genetic diversity and gene expression, with highly expressed genes experiencing stronger purifying selection. These genes cluster physically in the genome, revealing potential coordinated regulation. We identified known and novel candidate insecticide resistance genes, enriched in metabolic, cuticular, and behavioural functions. We also present AnoExpress, a Python package, and an online interface for user-friendly exploration of resistance candidate expression. Despite millions of years of speciation, convergent gene expression responses to insecticidal selection pressures are observed across Anopheles species, providing crucial insights for malaria vector control. This study culminates in a rich dataset that allows us to understand molecular mechanisms, better enabling us to combat insecticide resistance effectively.

Functional Characterization of CYP6QE1 and CYP6FV21 in Resistance to λ-Cyhalothrin and Imidacloprid in Bradysia odoriphaga

Cytochrome P450 monooxygenases (P450s) belong to a family of metabolic enzymes that are involved in the detoxification of insecticides. In this study, our bioassay results showed that a field-collected population of Bradysia odoriphaga displayed a moderate resistance to λ-cyhalothrin and imidacloprid. Compared to susceptible population, CYP6QE1 and CYP6FV21 were significantly overexpressed in the field population. The expression of CYP6QE1 and CYP6FV21 was more abundant in the third and fourth larval stages, and CYP6QE1 and CYP6FV21 were most highly expressed in the midgut and Malpighian tubules. Exposure to λ-cyhalothrin and imidacloprid significantly increased the expression levels of CYP6QE1 and CYP6FV21. Furthermore, the silencing of CYP6QE1 and CYP6FV21 significantly increased the susceptibility of B. odoriphaga larvae to λ-cyhalothrin and imidacloprid. The overexpression of CYP6QE1 and CYP6FV21 significantly enhanced the tolerance of transgenic Drosophila melanogaster lines to λ-cyhalothrin and imidacloprid. In addition, molecular docking revealed that these two P450 proteins have strong binding affinity toward λ-cyhalothrin and imidacloprid insecticides. Taken together, these results indicate that the overexpression of CYP6QE1 and CYP6FV21 is responsible for resistance to λ-cyhalothrin and imidacloprid in B. odoriphaga.

Parallel evolution in mosquito vectors - a duplicated esterase locus is associated with resistance to pirimiphos-methyl in An. gambiae

The primary control methods for the African malaria mosquito, Anopheles gambiae, are based on insecticidal interventions. Emerging resistance to these compounds is therefore of major concern to malaria control programmes. The organophosphate, pirimiphos-methyl, is a relatively new chemical in the vector control armoury but is now widely used in indoor residual spray campaigns. Whilst generally effective, phenotypic resistance has developed in some areas in malaria vectors. Here, we used a population genomic approach to identify novel mechanisms of resistance to pirimiphos-methyl in Anopheles gambiae s.l mosquitoes. In multiple populations, we found large and repeated signals of selection at a locus containing a cluster of detoxification enzymes, some of whose orthologs are known to confer resistance to organophosphates in Culex pipiens. Close examination revealed a pair of alpha-esterases, Coeae1f and Coeae2f, and a complex and diverse pattern of haplotypes under selection in An. gambiae, An. coluzzii and An. arabiensis. As in Cx. pipiens, copy number variation seems to play a role in the evolution of insecticide resistance at this locus. We used diplotype clustering to examine whether these signals arise from parallel evolution or adaptive introgression. Using whole-genome sequenced phenotyped samples, we found that in West Africa, a copy number variant in Anopheles gambiae is associated with resistance to pirimiphos-methyl. Overall, we demonstrate a striking example of contemporary parallel evolution which has important implications for malaria control programmes.

Heat shock proteins, thermotolerance, and insecticide resistance in mosquitoes

Mosquitoes transmit pathogens that pose a threat to millions of people globally. Unfortunately, widespread insecticide resistance makes it difficult to control these public health pests. General mechanisms of resistance, such as target site mutations or increased metabolic activity, are well established. However, many questions regarding the dynamics of these adaptations in the context of developmental and environmental conditions require additional exploration. One aspect of resistance that deserves further study is the role of heat shock proteins (HSPs) in insecticide tolerance. Studies show that mosquitoes experiencing heat stress before insecticide exposure demonstrate decreased mortality. This is similar to the observed reciprocal reduction in mortality in mosquitoes exposed to insecticide prior to heat stress. The environmental shifts associated with climate change will result in mosquitoes occupying environments with higher ambient temperatures, which could enhance existing insecticide resistance phenotypes. This physiological relationship adds a new dimension to the problem of insecticide resistance and further complicates the challenges that vector control and public health personnel face. This article reviews studies illustrating the relationship between insecticide resistance and HSPs or hsp genes as well as the intersection of thermotolerance and insecticide resistance. Further study of HSPs and insecticide resistance could lead to a deeper understanding of how environmental factors modulate the physiology of these important disease vectors to prepare for changing climatic conditions and the development of novel strategies to prevent vector-borne disease transmission.

Mechanisms of transcriptional regulation in Anopheles gambiae revealed by allele specific expression

Malaria control relies on insecticides targeting the mosquito vector, but this is increasingly compromised by insecticide resistance, which can be achieved by elevated expression of detoxifying enzymes that metabolize the insecticide. In diploid organisms, gene expression is regulated both in cis, by regulatory sequences on the same chromosome, and by trans acting factors, affecting both alleles equally. Differing levels of transcription can be caused by mutations in cis-regulatory modules (CRM), but few of these have been identified in mosquitoes. We crossed bendiocarb resistant and susceptible Anopheles gambiae strains to identify cis-regulated genes that might be responsible for the resistant phenotype using RNAseq, and cis-regulatory module sequences controlling gene expression in insecticide resistance relevant tissues were predicted using machine learning. We found 115 genes showing allele specific expression in hybrids of insecticide susceptible and resistant strains, suggesting cis regulation is an important mechanism of gene expression regulation in Anopheles gambiae. The genes showing allele specific expression included a higher proportion of Anopheles specific genes on average younger than genes those with balanced allelic expression.

Deltamethrin and transfluthrin select for distinct transcriptomic responses in the malaria vector Anopheles gambiae

BACKGROUND

The widespread use of pyrethroid insecticides in Africa has led to the development of strong resistance in Anopheles mosquitoes. Introducing new active ingredients can contribute to overcome this phenomenon and ensure the effectiveness of vector control strategies. Transfluthrin is a polyfluorinated pyrethroid whose structural conformation was thought to prevent its metabolism by cytochrome P450 monooxygenases in malaria vectors, thus representing a potential alternative for managing P450-mediated resistance occurring in the field. In this study, a controlled selection was used to compare the dynamics of resistance between transfluthrin and the widely used pyrethroid deltamethrin in the mosquito Anopheles gambiae. Then, the associated molecular mechanisms were investigated using target-site mutation genotyping and RNA-seq.

METHODS

A field-derived line of An. gambiae carrying resistance alleles at low frequencies was used as starting material for a controlled selection experiment. Adult females were selected across 33 generations with deltamethrin or transfluthrin, resulting in three distinct lines: the Delta-R line (selected with deltamethrin), the Transflu-R line (selected with transfluthrin) and the Tiassale-S line (maintained without selection). Deltamethrin and transfluthrin resistance levels were monitored in each selected line throughout the selection process, as well as the frequency of the L1014F kdr mutation. At generation 17, cross-resistance to other public health insecticides was investigated and transcriptomes were sequenced to compare gene transcription variations and polymorphisms associated with adaptation to each insecticide.

RESULTS

A rapid increase in resistance to deltamethrin and transfluthrin was observed throughout the selection process in each selected line in association with an increased frequency of the L1014F kdr mutation. Transcriptomic data support a broader response to transfluthrin selection as compared to deltamethrin selection. For instance, multiple detoxification enzymes and cuticle proteins were specifically over-transcribed in the Transflu-R line including the known pyrethroid metabolizers CYP6M2, CYP9K1 and CYP6AA1 together with other genes previously associated with resistance in An. gambiae.

CONCLUSION

This study confirms that recurrent exposure of adult mosquitoes to pyrethroids in a public health context can rapidly select for various resistance mechanisms. In particular, it indicates that in addition to target site mutations, the polyfluorinated pyrethroid transfluthrin can select for a broad metabolic response, which includes some P450s previously associated to resistance to classical pyrethroids. This unexpected finding highlights the need for an in-depth study on the adaptive response of mosquitoes to newly introduced active ingredients in order to effectively guide and support decision-making programmes in malaria control.

Chlorpyrifos and chlorfenapyr resistance in Spodoptera frugiperda (Lepidoptera: Noctuidae) relies on UDP-glucuronosyltransferases

Fall armyworm, Spodoptera frugiperda (J. E. Smith), has become an important agricultural pest worldwide. S. frugiperda is mainly controlled by the chemical insecticides, whereas the frequent application of insecticides would result in the resistance development. Insect uridine diphosphate-glucuronosyltransferases (UGTs), as phase II metabolism enzymes, play vital roles in the breakdown of endobiotic and xenobiotics. In this study, 42 UGT genes were identified by RNA-seq, including 29 UGT genes were elevated compared to the susceptible population, and the transcript levels of 3 UGTs (UGT40F20, UGT40R18, and UGT40D17) were increased by more than 2.0-fold in the field populations. Expression pattern analysis revealed that S. frugiperda UGT40F20, UGT40R18, and UGT40D17 were increased by 6.34-, 4.26-, and 8.28-fold, compared the susceptible populations, respectively. The expression of UGT40D17, UGT40F20, and UGT40R18 was affected after exposure to phenobarbital, chlorpyrifos, chlorfenapyr, sulfinpyrazone, and 5-nitrouracil. The induced expression of UGT genes may have improved UGT enzymatic activity, while the inhibition of UGTs genes expression may decreased UGT enzymatic activity. Sulfinpyrazone, and 5-nitrouracil, significantly increased the toxicity of chlorpyrifos and chlorfenapyr, as well as phenobarbital significantly reduced the toxicity of chlorpyrifos and chlorfenapyr against the susceptible populations and field populations of S. frugiperda. The suppression of UGTs (UGT40D17, UGT40F20, and UGT40R18) significantly increased the insensitivity of the field populations to chlorpyrifos and chlorfenapyr. These findings strongly supported our viewpoint that UGTs may play a critical role in insecticide detoxification. This study provides a scientific basis for the management of S. frugiperda.

The chemosensory protein 1 contributes to indoxacarb resistance in Plutella xylostella (L.)

BACKGROUND

Insecticide resistance continuously poses a threat to agricultural production. Chemosensory protein-mediated resistance is a new mechanism of insecticide resistance discovered in recent years. In-depth research on resistance mediated by chemosensory proteins (CSPs) provides new insight into aid insecticide resistance management.

RESULTS

Chemosensory protein 1 in Plutella xylostella (PxCSP1) was overexpressed in the two indoxacarb-resistant field populations and PxCSP1 has a high affinity with indoxacarb. PxCSP1 was upregulated when exposed to indoxacarb and the knockdown of this gene elevated sensitivity to indoxacarb, which demonstrate that PxCSP1 is involved in the indoxacarb resistance. Considering that CSPs may confer resistance in insects via binding or sequestering, we explored the binding mechanism of indoxacarb in PxCSP1-mediated resistance. Using molecular dynamics simulations and site-directed mutation, we found that indoxacarb forms a solid complex with PxCSP1 mainly through van der Waals interactions and electrostatic interactions. Between these, the electrostatic interaction provided by the Lys100 side chain in PxCSP1, and especially the hydrogen bonding between the NZ atom and the O of the carbamoyl carbonyl group of indoxacarb, are the key factors for the high affinity of PxCSP1 to indoxacarb.

CONCLUSIONS

The overexpression of PxCPS1 and its high affinity to indoxacarb is partially responsible for indoxacarb resistance in P. xylostella. Modification of indoxacarb's carbamoyl group has the potential to alleviate indoxacarb resistance in P. xylostella. These findings will contribute to solving chemosensory protein-mediated indoxacarb resistance and provide a better understanding of the insecticide resistance mechanism. © 2023 Society of Chemical Industry.

Characterising the intensity of insecticide resistance: A novel framework for analysis of intensity bioassay data

Insecticide resistance is a growing problem that risks harming the progress made by vector control tools in reducing the malaria burden globally. New methods for quantifying the extent of resistance in wild populations are urgently needed to guide deployment of interventions to improve disease control. Intensity bioassays measure mosquito mortality at a range of insecticide doses and characterise phenotypic resistance in regions where resistance is already detected. These data are increasingly being collected but tend to exhibit high measurement error and there is a lack of formal guidelines on how they should be analysed or compared. This paper introduces a novel Bayesian framework for analysing intensity bioassay data, which uses a flexible statistical model able to capture a wide variety of relationships between mortality and insecticide dose. By accounting for background mortality of mosquitoes, our approach minimises the impact of this source of measurement noise resulting in more precise quantification of resistance. It outputs a range of metrics for describing the intensity and variability in resistance within the sample and quantifies the level of measurement error in the assay. The functionality is illustrated with data from laboratory-reared mosquitoes to show how the lethal dose varies within and between different strains. The framework can also be used to formally test hypotheses by explicitly considering the high heterogeneity seen in these types of data in field samples. Here we show that the intensity of resistance (as measured by the median lethal dose (LC50) of insecticide) increases over 7 years in mosquitoes from one village in Burkina Faso but remains constant in another. This work showcases the benefits of statistically rigorous analysis of insecticide bioassay data and highlights the additional information available from this and other dose-response data.

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

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

Chemosensory proteins are associated with thiamethoxam tolerance in bird cherry-oat aphid Rhopalosiphum padi

Rhopalosiphum padi (L.) is an important cosmopolitan pest of cereal crops. Thiamethoxam is widely used for control R. padi in some regions. Chemosensory proteins (CSPs) are a class of transporter proteins in arthropods which play a key role in various physiological processes including response to insecticide exposure. However, the role of R. padi CSPs (RpCSPs) in insecticide binding and susceptibility has not been well clarified. In this study, we found that the expression levels of RpCSP1, RpCSP4, RpCSP5, RpCSP7, RpCSP10 were dramatically upregulated after exposure to thiamethoxam. Suppression of RpCSP4 and RpCSP5 transcription by RNA interference significantly enhanced the susceptibility of R. padi to thiamethoxam. Molecular docking and fluorescence competitive binding showed that RpCSP4 and RpCSP5 had high binding affinity with thiamethoxam. The present results prove that RpCSP4 and RpCSP5 are related to insecticide resistance through high binding affinity to reduce the toxicity of insecticide.

Transcriptome profile of halofuginone resistant and sensitive strains of Eimeria tenella

The antiparasitic drug halofuginone is important for controlling apicomplexan parasites. However, the occurrence of halofuginone resistance is a major obstacle for it to the treatment of apicomplexan parasites. Current studies have identified the molecular marker and drug resistance mechanisms of halofuginone in Plasmodium falciparum. In this study, we tried to use transcriptomic data to explore resistance mechanisms of halofuginone in apicomplexan parasites of the genus Eimeria (Apicomplexa: Eimeriidae). After halofuginone treatment of E. tenella parasites, transcriptome analysis was performed using samples derived from both resistant and sensitive strains. In the sensitive group, DEGs associated with enzymes were significantly downregulated, whereas the DNA damaging process was upregulated after halofuginone treatment, revealing the mechanism of halofuginone-induced parasite death. In addition, 1,325 differentially expressed genes (DEGs) were detected between halofuginone resistant and sensitive strains, and the DEGs related to translation were significantly downregulated after halofuginone induction. Overall, our results provide a gene expression profile for further studies on the mechanism of halofuginone resistance in E. tenella.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Comparative analysis of detoxification-related gene superfamilies across five hemipteran species

Background

Hemiptera is one of the most speciose orders of insects, and the most speciose considering Hemimetabola. Through their evolutive history, hemipterans with different feeding habits have adapted to deal with different chemical challenges. Three major gene families are involved in xenobiotic detoxification in insects: the cytochromes P450 (CYPs), carboxyl/cholinesterases (CCEs), and glutathione transferases (GSTs). Here we perform a comparative analysis on the complement of these gene superfamilies across five hemipteran species; four heteropterans (the pentatomid plant feeders Nezara viridula and Halyomorpha halys; the hematophagous Cimex lectularius, Cimicidae, and Rhodnius prolixus, Reduviidae), and one Auchenorrhyncha plant feeder (Nilaparvata lugens).

Results

Our results point to an expansion of several enzyme families associated with xenobiotic detoxification in heteropterans with respect to other species and the existence of a dynamic evolution pattern including CYP3 clan, hormone and pheromone processing class in the CCE superfamily, and sigma class in GST superfamily. Other detoxification-related families are reduced in the hemipteran species analyzed here: reduction or even absence of epsilon class and reduced delta class in GST superfamily; absence of mitochondrial CYP12 family; absence of CYP9 family in CYP3 clan; and reduction or even absence of some dietary/detoxification groups of CCEs. Interestingly, the most polyphagous species analyzed here (H. halys) is also the one that presents the largest repertoire of detoxification enzymes. Gene cluster analysis suggests that this could be due to gene duplication events.

Conclusions

The evolutionary analysis performed here reveals characteristics that are both common and particular for heteropterans. The composition and organization of detoxification-related gene families could shed light on evolutionary forces that shaped their divergence. These families are important for both the detoxification of diet products and for conferring tolerance or resistance to synthetic insecticides. Furthermore, we present the first comprehensive analysis of detoxification gene superfamilies in N. viridula, an understudied species in spite of its economic relevance as a crop pest. The information obtained is of interest for basic insect science as well as for the control of harmful species and the management of insecticide resistance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08974-y.

Functional Inquiry into ATP-Binding Cassette Transporter Genes Contributing to Spirotetramat Resistance in Aphis gossypii Glover

ATP-binding cassette (ABC) transporters regulate the efflux of a broad spectrum of substrates to extracellular transporting, which play an important role in the detoxification process in arthropods. Here, we described a comprehensive approach to explore the involvement of ABC transporters in spirotetramat resistance in cotton aphids. In this study, synergism bioassays showed 17.05% and 35.42% increases in the toxicity to spirotetramat with the ABC inhibitor verapamil in adult and 3rd instar nymph aphids of the SR strain, respectively. In a competitive assay based on the microinjection of a fluorescent ABC transporter substrate, verapamil (a general ABC inhibitor) and spirotetramat significantly inhibited the elimination of Texas Red. Based on transcriptome data of midguts of spirotetramat-susceptible (SS) and -resistant (SR) strains, the expression levels of ABCB4, ABCB5, ABCF2, MRP11, and MRP12 were significantly upregulated in the SR strain midgut compared to that of the SS strain. Gene functional analysis based on ectopic expression and RNA interference (RNAi) proved that ABCB4, ABCB5, ABCF2, MRP11, and MRP12 were involved in the tolerance of cotton aphids to spirotetramat. Moreover, the upregulated ABCF2, ABCB4, and ABCB5 in the midgut of the SR strain contributed more to the resistance of spirotetramat in in vitro functional analysis. In summary, these results demonstrate that candidate ABC transporter genes in the midgut tissue were involved in spirotetramat resistance, which will help reveal the relationship between ABC transporters and the development of spirotetramat resistance in field populations.

Exposure of Anopheles gambiae larvae to a sub-lethal dose of an agrochemical mixture induces tolerance to adulticides used in vector control management

The heavy use of pesticides in agricultural areas often leads to the contamination of nearby mosquito larvae breeding sites. Exposure to complex mixtures of agrochemicals can affect the insecticide sensitivity of mosquito larvae. Our study objective was to determine whether agrochemical residues in Anopheline larval breeding sites can affect the tolerance of adults to commonly used adulticides. We focussed on Fludora® Fusion, a vector control insecticide formulation combining two insecticides (deltamethrin and clothianidin) with different modes of action. An. gambiae larvae were exposed to a sub-lethal dose of a mixture of agrochemical pesticides used in a highly active agricultural area on the Ivory Coast. Comparative bioassays with Fludora Fusion mixture and its two insecticide components (deltamethrin and clothianidin) were carried out between adult mosquitoes exposed or not to the agrochemicals at the larval stage. A transcriptomic analysis using RNA sequencing was then performed on larvae and adults to study the molecular mechanisms underlying the phenotypic changes observed. Bioassays revealed a significantly increased tolerance of adult females to clothianidin (2.5-fold) and Fludora Fusion mixture (2.2-fold) following larval exposure to agrochemicals. Significantly increased tolerance to deltamethrin was not observed suggesting that insecticide exposure affects the adult efficacy of the Fludora Fusion mixture mainly through mechanisms acting on clothianidin. Transcriptomic analysis revealed the potential of agrochemicals to induce various resistance mechanisms including cuticle proteins, detoxification action and altered insecticide sequestration. These results suggest that although the Fludora Fusion mixture is effective for adult vector control, its efficacy may be locally affected by the ecological context. The present study also suggests that, although the complex interactions between the use of agrochemicals and vector control insecticides are difficult to decipher in the field, they still must be considered in the context of insecticide resistance management programmes.

Transcriptomic modulation in response to an intoxication with deltamethrin in a population of Triatoma infestans with low resistance to pyrethroids

Background

Triatoma infestans is the main vector of Chagas disease in the Southern Cone. The resistance to pyrethroid insecticides developed by populations of this species impairs the effectiveness of vector control campaigns in wide regions of Argentina. The study of the global transcriptomic response to pyrethroid insecticides is important to deepen the knowledge about detoxification in triatomines.

Methodology and findings

We used RNA-Seq to explore the early transcriptomic response after intoxication with deltamethrin in a population of T. infestans which presents low resistance to pyrethroids. We were able to assemble a complete transcriptome of this vector and found evidence of differentially expressed genes belonging to diverse families such as chemosensory and odorant-binding proteins, ABC transporters and heat-shock proteins. Moreover, genes related to transcription and translation, energetic metabolism and cuticle rearrangements were also modulated. Finally, we characterized the repertoire of previously uncharacterized detoxification-related gene families in T. infestans and Rhodnius prolixus.

Conclusions and significance

Our work contributes to the understanding of the detoxification response in vectors of Chagas disease. Given the absence of an annotated genome from T. infestans, the analysis presented here constitutes a resource for molecular and physiological studies in this species. The results increase the knowledge on detoxification processes in vectors of Chagas disease, and provide relevant information to explore undescribed potential insecticide resistance mechanisms in populations of these insects.

Chemosensory proteins confer adaptation to the ryanoid anthranilic diamide insecticide cyantraniliprole in Aphis gossypii glover

Chemosensory proteins (CSPs) are a class of small transporter proteins expressed only in arthropods with various functions beyond chemoreception. Previous studies have been reported that CSPs are involved in the insecticide resistance. In this study, we found that AgoCSP1, AgoCSP4, and AgoCSP5 were constitutively overexpressed in an insecticide-resistant strain of Aphis gossypii and showed higher expression in broad body tissue (including fat bodies) than in the midgut but without tissue specificity. However, the function of these three upregulated AgoCSPs remains unknown. Here, we investigated the function of AgoCSPs in resistance to the diamide insecticide cyantraniliprole. Suppression of AgoCSP1, AgoCSP4 and AgoCSP5 transcription by RNAi significantly increased the sensitivity of resistant aphids to cyantraniliprole. Molecular docking and competitive binding assays indicated that these AgoCSPs bind moderate with cyantraniliprole. Transgenic Drosophila melanogaster expressing these AgoCSPs in the broad body or midgut showed higher tolerance to cyantraniliprole than control flies with the same genetic background; AgoCSP4 was more effective in broad body tissue, and AgoCSP1 and AgoCSP5 were more effective in the midgut, indicating that broad body and midgut tissues may be involved in the insecticide resistance mediated by the AgoCSPs examined. The present results strongly indicate that AgoCSPs participate in xenobiotic detoxification by sequestering and masking toxic insecticide molecules, providing insights into new factors involved in resistance development in A. gossypii.

Functional investigation of lncRNAs and target cytochrome P450 genes related to spirotetramat resistance in Aphis gossypii Glover

BACKGROUND

Spirotetramat is a tetramic acid derivative insecticide with novel modes of action for controlling Aphis gossypii Glover in the field. Previous studies have shown that long noncoding RNAs (lncRNAs) and cytochrome P450 monooxygenases (P450s) are involved in the detoxification process. However, the functions of lncRNAs in regulating P450 gene expression in spirotetramat resistance in A. gossypii are unknown.

RESULTS

In this study, we found CYP4CJ1, CYP6CY7 and CYP6CY21 expression levels to be significantly upregulated in a spirotetramat-resistant (SR) strain compared with a susceptible (SS) strain. Furthermore, knockdown of CYP4CJ1, CYP6CY7 and CYP6CY21 increased nymph and adult mortality in the SR strain following exposure to spirotetramat. Drosophila ectopically expressing CYP380C6, CYP4CJ1, CYP6DA2, CYP6CY7 and CYP6CY21 showed significantly decreased mortality after spirotetramat exposure, and CYP380C6, CYP4CJ1 and CYP6CY21 are putative targets of six lncRNAs. Silencing of lncRNAs MSTRG.36649.2/5 and MSTRG.71880.1 changed CYP6CY21 and CYP380C6 expression, altering the sensitivity of the SR strain to spirotetramat. Moreover, MSTRG.36649.2/5 did not compete for microRNA (miRNA) binding to regulate CYP6CY21 expression.

CONCLUSION

Our results confirm that CYP380C6, CYP4CJ1, CYP6DA2, CYP6CY7 and CYP6CY21 are potentially involved in the development of spirotetramat resistance in A. gossypii, and MSTRG.36649.2/5 and MSTRG.71880.1 probably regulate CYP6CY21 and CYP380C6 expression other than through the "sponge effect" of competing for miRNA binding. Our results provide a favorable molecular basis for studying cotton aphid P450 genes and lncRNA functions in spirotetramat resistance development.

V101I and R81T mutations in the nicotinic acetylcholine receptor β1 subunit are associated with neonicotinoid resistance in Myzus persicae

BACKGROUND

The peach-potato aphid, Myzus persicae (Sulzer) (Hemiptera: Aphididae), is a major pest worldwide. The intensive use of insecticides has led to the development of resistance against neonicotinoid insecticides. The R81T mutation in the nicotinic acetylcholine receptor (nAChR) beta1 subunit is considered a crucial mechanism adaptation to neonicotinoid resistance in M. persicae and Aphis gossypii.

RESULTS

Resistance-related mutations (R81T and V101I) were detected in the imidacloprid-resistant M. persicae AH19 population. The V101I mutation is reported for the first time. The V101I and R81T mutations existed separately, indicating that the two mutations evolved independently. Imidacloprid resistance in the AH19 population was stable without insecticide exposure. Four mutant strains were selected from the population with stable resistance. The resistance of the AH19-T, AH19-I, and AH19-TI strains to imidacloprid, thiamethoxam, and dinotefuran was significantly increased compared with the AH19-W strain. Synergism bioassays showed that the inhibition of three detoxification enzymes did not affect imidacloprid resistance in the AH19-T and AH19-I strains. Expression of nAChR β1 subunits in the AH19-W, AH19-T, and AH19-I strains remained unchanged.

CONCLUSION

The V101I mutation is associated with neonicotinoid resistance in M. persicae. The resistance of the AH19-T and AH19-I strains to neonicotinoids appears to be mainly due to the R81T and V101I mutations, whereas these mutations, together with changes in the cytochrome P450 monooxygenases and nAChR expression may be responsible for the development of neonicotinoid resistance in the AH19-TI strain. © 2021 Society of Chemical Industry.

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

Simple Summary

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

Abstract

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

Chemosensory Proteins Are Associated with Thiamethoxam and Spirotetramat Tolerance in Aphis gossypii Glover

Chemosensory proteins (CSPs) are a class of transporters in arthropods. Deeper research on CSPs showed that CSPs may be involved in some physiological processes beyond chemoreception, such as insect resistance to pesticides. We identified two upregulated CSPs in two resistant strains of Aphis gossypii Glover. To understand their role in the resistance of aphids to pesticides, we performed the functional verification of CSP1 and CSP4 in vivo and in vitro. Results showed that the sensitivity of the thiamethoxam-resistant strain to thiamethoxam increased significantly with the silencing of CSP1 and CSP4 by RNAi (RNA interference), and the sensitivity of the spirotetramat-resistant strain to spirotetramat increased significantly with the silencing of CSP4. Transgenic Drosophila melanogaster expressing CSPs exhibited stronger resistance to thiamethoxam, spirotetramat, and alpha-cypermethrin than the control did. In the bioassay of transgenic Drosophila, CSPs showed different tolerance mechanisms for different pesticides, and the overexpressed CSPs may play a role in processes other than resistance to pesticides. In brief, the present results prove that CSPs are related to the resistance of cotton aphids to insecticides.

Susceptibility of the Cowpea Aphid (Hemiptera: Aphididae) to Widely Used Insecticides in Australia

Globally, 27 aphid species have evolved resistance to almost 100 insecticide active ingredients. A proactive approach to resistance management in pest aphids is needed; this should include risk analysis, followed by regular baseline susceptibility assays for species deemed at high risk of evolving resistance. The cowpea aphid (Aphis craccivora Koch) has evolved insecticide resistance to multiple insecticides outside Australia and was recently identified as a high-risk species in Australia. In this study, we generated toxicity data against four insecticides (representing four unique chemical Mode of Action groups) for populations of A. craccivora collected across Australia. Alpha-cypermethrin was the most toxic chemical to A. craccivora in leaf-dip laboratory bioassays with an average LC50 value across nine populations of 0.008 mg a.i./L, which was significantly lower than dimethoate (1.17 mg a.i./L) and pirimicarb (0.89 mg a.i./L). Small, but significant, differences in sensitivity were detected in some populations against pirimicarb and dimethoate, whereas responses to alpha-cypermethrin and imidacloprid were not significantly different across all aphid populations examined in this study. For all insecticides, the field rate controlled 100% of individuals tested. The data generated will be important for future monitoring of insecticide responses of A. craccivora. Proactive management, including increased reliance on non-chemical pest management approaches and routine insecticide baseline sensitivity studies, is recommended for A. craccivora.

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

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

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

Genome resequencing reveals rapid, repeated evolution in the Colorado potato beetle

Insecticide resistance and rapid pest evolution threatens food security and the development of sustainable agricultural practices, yet the evolutionary mechanisms that allow pests to rapidly adapt to control tactics remains unclear. Here we examine how a global super-pest, the Colorado potato beetle (CPB), Leptinotarsa decemlineata, rapidly evolves resistance to insecticides. Using whole genome resequencing and transcriptomic data focused on its ancestral and pest range in North America, we assess evidence for three, non-mutually exclusive models of rapid evolution: pervasive selection on novel mutations, rapid regulatory evolution, and repeated selection on standing genetic variation. Population genomic analysis demonstrates that CPB is geographically structured, even among recently established pest populations. Pest populations exhibit similar levels of nucleotide diversity, relative to non-pest populations, and show evidence of recent expansion. Genome scans provide clear signatures of repeated adaptation across CPB populations, with especially strong evidence of selection on insecticide resistance genes in different populations. Analyses of gene expression show that constitutive upregulation of candidate insecticide resistance genes drives distinctive population patterns. CPB evolves insecticide resistance repeatedly across agricultural regions, leveraging similar genetic pathways but different genes, demonstrating a polygenic trait architecture for insecticide resistance that can evolve from standing genetic variation. Despite expectations, we do not find support for strong selection on novel mutations, or rapid evolution from selection on regulatory genes. These results suggest that integrated pest management practices must mitigate the evolution of polygenic resistance phenotypes among local pest populations, in order to maintain the efficacy and sustainability of novel control techniques.

Transcriptomic analysis of resistance and short-term induction response to pyrethroids, in Anopheles coluzzii legs

BACKGROUND

Insecticide-treated bed nets and indoor residual spraying comprise the major control measures against Anopheles gambiae sl, the dominant vector in sub-Saharan Africa. The primary site of contact with insecticide is through the mosquitoes' legs, which represents the first barrier insecticides have to bypass to reach their neuronal targets. Proteomic changes and leg cuticle modifications have been associated with insecticide resistance that may reduce the rate of penetration of insecticides. Here, we performed a multiple transcriptomic analyses focusing on An. coluzzii legs.

RESULTS

Firstly, leg-specific enrichment analysis identified 359 genes including the pyrethroid-binder SAP2 and 2 other chemosensory proteins, along with 4 ABCG transporters previously shown to be leg enriched. Enrichment of gene families included those involved in detecting chemical stimuli, including gustatory and ionotropic receptors and genes implicated in hydrocarbon-synthesis. Subsequently, we compared transcript expression in the legs of a highly resistant strain (VK7-HR) to both a strain with very similar genetic background which has reverted to susceptibility after several generations without insecticide pressure (VK7-LR) and a lab susceptible population (NG). Two hundred thirty-two differentially expressed genes (73 up-regulated and 159 down-regulated) were identified in the resistant strain when compared to the two susceptible counterparts, indicating an over-expression of phase I detoxification enzymes and cuticular proteins, with decrease in hormone-related metabolic processes in legs from the insecticide resistant population. Finally, we analysed the short-term effect of pyrethroid exposure on An. coluzzii legs, comparing legs of 1 h-deltamethrin-exposed An. coluzzii (VK7-IN) to those of unexposed mosquitoes (VK7-HR) and identified 348 up-regulated genes including those encoding for GPCRs, ABC transporters, odorant-binding proteins and members of the divergent salivary gland protein family.

CONCLUSIONS

The data on An. coluzzii leg-specific transcriptome provides valuable insights into the first line of defense in pyrethroid resistant and short-term deltamethrin-exposed mosquitoes. Our results suggest that xenobiotic detoxification is likely occurring in legs, while the enrichment of sensory proteins, ABCG transporters and cuticular genes is also evident. Constitutive resistance is primarily associated with elevated levels of detoxification and cuticular genes, while short-term insecticide-induced tolerance is linked with overexpression of transporters, GPCRs and GPCR-related genes, sensory/binding and salivary gland proteins.

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

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

Overabundance of Asaia and Serratia Bacteria Is Associated with Deltamethrin Insecticide Susceptibility in Anopheles coluzzii from Agboville, Côte d'Ivoire

Insecticide resistance among mosquito species is now a pervasive phenomenon that threatens to jeopardize global malaria vector control efforts. Evidence of links between the mosquito microbiota and insecticide resistance is emerging, with significant enrichment of insecticide degrading bacteria and enzymes in resistant populations. Using 16S rRNA amplicon sequencing, we characterized and compared the microbiota of Anopheles coluzzii in relation to their deltamethrin resistance and exposure profiles. Comparisons between 2- and 3-day-old deltamethrin-resistant and -susceptible mosquitoes demonstrated significant differences in microbiota diversity. Ochrobactrum, Lysinibacillus, and Stenotrophomonas genera, each of which comprised insecticide-degrading species, were significantly enriched in resistant mosquitoes. Susceptible mosquitoes had a significant reduction in alpha diversity compared to resistant individuals, with Asaia and Serratia dominating microbial profiles. There was no significant difference in deltamethrin-exposed and -unexposed 5- to 6-day-old individuals, suggesting that insecticide exposure had minimal impact on microbial composition. Serratia and Asaia were also dominant in 5- to 6-day-old mosquitoes, which had reduced microbial diversity compared to 2- to 3-day-old mosquitoes. Our findings revealed significant alterations of Anopheles coluzzii microbiota associated with deltamethrin resistance, highlighting the potential for identification of novel microbial markers for insecticide resistance surveillance. qPCR detection of Serratia and Asaia was consistent with 16S rRNA sequencing, suggesting that population-level field screening of bacterial microbiota may be feasibly integrated into wider resistance monitoring, if reliable and reproducible markers associated with phenotype can be identified.

IMPORTANCE Control of insecticide-resistant vector populations remains a significant challenge to global malaria control and while substantial progress has been made elucidating key target site mutations, overexpressed detoxification enzymes and alternate gene families, the contribution of the mosquito microbiota to phenotypic insecticide resistance has been largely overlooked. We focused on determining the effects of deltamethrin resistance intensity on Anopheles coluzzii microbiota and identifying any microbial taxa associated with phenotype. We demonstrated a significant reduction in microbial diversity between deltamethrin-resistant and -susceptible mosquitoes. Insecticide degrading bacterial species belonging to Ochrobactrum, Lysinibacillus, and Stenotrophomonas genera were significantly enriched in resistant mosquitoes, while Asaia and Serratia dominated microbial profiles of susceptible individuals. Our results revealed significant alterations of Anopheles coluzzii microbiota associated with deltamethrin resistance, highlighting the potential for identification of novel microbial markers for surveillance and opportunities for designing innovative control techniques to prevent the further evolution and spread of insecticide resistance.

Development of interactive biological web applications with R/Shiny

Development of interactive web applications to deposit, visualize and analyze biological datasets is a major subject of bioinformatics. R is a programming language for data science, which is also one of the most popular languages used in biological data analysis and bioinformatics. However, building interactive web applications was a great challenge for R users before the Shiny package was developed by the RStudio company in 2012. By compiling R code into HTML, CSS and JavaScript code, Shiny has made it incredibly easy to build web applications for the large R community in bioinformatics and for even non-programmers. Over 470 biological web applications have been developed with R/Shiny up to now. To further promote the utilization of R/Shiny, we reviewed the development of biological web applications with R/Shiny, including eminent biological web applications built with R/Shiny, basic steps to build an R/Shiny application, commonly used R packages to build the interface and server of R/Shiny applications, deployment of R/Shiny applications in the cloud and online resources for R/Shiny.

Experimental evolution supports the potential of neonicotinoid-pyrethroid combination for managing insecticide resistance in malaria vectors

The introduction of neonicotinoids for managing insecticide resistance in mosquitoes is of high interest as they interact with a biochemical target not previously used in public health. In this concern, Bayer developed a combination of the neonicotinoid clothianidin and the pyrethroid deltamethrin (brand name Fludora Fusion) as a new vector control tool. Although this combination proved to be efficient against pyrethroid-resistant mosquitoes, its ability to prevent the selection of pyrethroid and neonicotinoid resistance alleles was not investigated. In this context, the objective of this work was to study the dynamics and the molecular mechanisms of resistance of An. gambiae to the separated or combined components of this combination. A field-derived An. gambiae line carrying resistance alleles to multiple insecticides at low frequencies was used as a starting for 33 successive generations of controlled selection. Resistance levels to each insecticide and target site mutation frequencies were monitored throughout the selection process. Cross resistance to other public health insecticides were also investigated. RNA-seq was used to compare gene transcription variations and polymorphisms across all lines. This study confirmed the potential of this insecticide combination to impair the selection of resistance as compared to its two separated components. Deltamethrin selection led to the rapid enrichment of the kdr L1014F target-site mutation. Clothianidin selection led to the over-transcription of multiple cytochrome P450s including some showing high homology with those conferring neonicotinoid resistance in other insects. A strong selection signature associated with clothianidin selection was also observed on a P450 gene cluster previously associated with resistance. Within this cluster, the gene CYP6M1 showed the highest selection signature together with a transcription profile supporting a role in clothianidin resistance. Modelling the impact of point mutations selected by clothianidin on CYP6M1 protein structure showed that selection retained a protein variant with a modified active site potentially enhancing clothianidin metabolism. In the context of the recent deployment of neonicotinoids for mosquito control and their frequent usage in agriculture, the present study highlights the benefit of combining them with other insecticides for preventing the selection of resistance and sustaining vector control activities.