The choreography of the chemical defensome response to insecticide stress: insights into the Anopheles stephensi transcriptome using RNA-Seq

Integrated de novo transcriptome of Culex pipiens mosquito larvae as a resource for genetic control strategies

We present a de novo transcriptome of the mosquito vector Culex pipiens, assembled by sequences of susceptible and insecticide resistant larvae. The high quality of the assembly was confirmed by TransRate and BUSCO. A mapping percentage until 94.8% was obtained by aligning contigs to Nr, SwissProt, and TrEMBL, with 27,281 sequences that simultaneously mapped on the three databases. A total of 14,966 ORFs were also functionally annotated by using the eggNOG database. Among them, we identified ORF sequences of the main gene families involved in insecticide resistance. Therefore, this resource stands as a valuable reference for further studies of differential gene expression as well as to identify genes of interest for genetic-based control tools.

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.

Time-series analysis of transcriptomic changes due to permethrin exposure reveals that Aedes aegypti undergoes detoxification metabolism over 24 h

Insecticide resistance is a multifaceted response and an issue across taxa. Aedes aegypti, the mosquito that vectors Zika, dengue, chikungunya, and yellow fever, demonstrates high levels of pyrethroid resistance across the globe, presenting a challenge to public health officials. To examine the transcriptomic shifts across time after exposure to permethrin, a 3'Tag-Seq analysis was employed on samples 6, 10, and 24 h after exposure along with controls. Differential expression analysis revealed significant shifts in detoxifying enzymes and various energy-producing metabolic processes. These findings indicate significant alterations in gene expression associated with key energy mobilization pathways within the system. These changes encompass a coordinated response involving lipolysis, beta-oxidation, and the citric acid cycle, required for the production of energetic molecules such as ATP, NADH, NADPH, and FADH. These findings highlight a complex interplay of metabolic processes that may have broader implications for understanding insect physiology and response to environmental stimuli. Among the upregulated detoxifying enzymes are cytochrome P450s, glutathione s-transferases and peroxidases, and ATP-binding cassette transporters. Additionally, eight heat shock genes or genes with heat shock domains exhibit the highest fold change across time. Twenty-four hours after exposure, samples indicate a global downregulation of these processes, though principal component analysis suggests lasting signatures of the response. Understanding the recovery response to insecticide exposure provides information on possible new genetic and synergist targets to explore.

ABCH2 transporter mediates deltamethrin uptake and toxicity in the malaria vector Anopheles coluzzii

Contact insecticides are primarily used for the control of Anopheles malaria vectors. These chemicals penetrate mosquito legs and other appendages; the first barriers to reaching their neuronal targets. An ATP-Binding Cassette transporter from the H family (ABCH2) is highly expressed in Anopheles coluzzii legs, and further induced upon insecticide exposure. RNAi-mediated silencing of the ABCH2 caused a significant increase in deltamethrin mortality compared to control mosquitoes, coincident with a corresponding increase in 14C-deltamethrin penetration. RT-qPCR analysis and immunolocalization revealed ABCH2 to be mainly localized in the legs and head appendages, and more specifically, the apical part of the epidermis, underneath the cuticle. To unravel the molecular mechanism underlying the role of ABCH2 in modulating pyrethroid toxicity, two hypotheses were investigated: An indirect role, based on the orthology with other insect ABCH transporters involved with lipid transport and deposition of CHC lipids in Anopheles legs which may increase cuticle thickness, slowing down the penetration rate of deltamethrin; or the direct pumping of deltamethrin out of the organism. Evaluation of the leg cuticular hydrocarbon (CHC) content showed no affect by ABCH2 silencing, indicating this protein is not associated with the transport of leg CHCs. Homology-based modeling suggested that the ABCH2 half-transporter adopts a physiological homodimeric state, in line with its ability to hydrolyze ATP in vitro when expressed on its own in insect cells. Docking analysis revealed a deltamethrin pocket in the homodimeric transporter. Furthermore, deltamethrin-induced ATP hydrolysis in ABCH2-expressing cell membranes, further supports that deltamethrin is indeed an ABCH2 substrate. Overall, our findings pinpoint ABCH2 participating in deltamethrin toxicity regulation.

Profiling of Dermanyssus gallinae genes involved in acaricide resistance

The poultry red mite (PRM), Dermanyssus gallinae, is a major threat for the poultry industry worldwide. Chemical compounds have been extensively used for PRM control, leading to selection of resistant mites. Molecular mechanisms of resistance have been investigated in arthropods, showing the role of target-site insensitivity and enhanced detoxification. Few studies are available about those mechanisms in D. gallinae, and none have yet focused on the expression levels of detoxification enzymes and other defense-related genes through RNA-seq. We tested PRM populations from Italy for their susceptibility to the acaricidal compounds phoxim and cypermethrin. Mutations in the voltage-gated sodium channel (vgsc) and in acetylcholinesterase (AChE) were investigated, detecting mutations known to be associated with acaricide/insecticide resistance in arthropods, including M827I and M918L/T in the vgsc and G119S in the AChE. RNA-seq analysis was performed to characterize metabolic resistance in fully susceptible PRM and in cypermethrin-resistant PRM exposed and unexposed to cypermethrin as well as phoxim resistant PRM exposed and unexposed to phoxim. Detoxification enzymes (including P450 monooxygenases and glutathione-S-transferases), ABC transporters and cuticular proteins were constitutively overexpressed in phoxim and cypermethrin resistant mites. In addition, heat shock proteins were found constitutively and inductively upregulated in phoxim resistant mites, while in cypermethrin resistant mites esterases and an aryl hydrocarbon receptor were constitutively highly expressed. The findings suggest that acaricide resistance in D. gallinae is due to both target-site insensitivity and overexpression of detoxification enzymes and other xenobiotic defense-related genes, which is mostly constitutive and not induced by treatment. Understanding the molecular basis of resistance could be useful to screen or test PRM populations in order to select targeted acaricides and to avoid the abuse/misuse of the few available compounds.

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

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

Transcriptome and Metabolome Analysis Reveals the Importance of Amino-Acid Metabolism in Spodoptera Frugiperda Exposed to Spinetoram

Pests are inevitably exposed to sublethal and lethal doses in the agroecosystem following the application of pesticides indispensable to protect food sources. The effect of spinetoram on amino-acid metabolism of fall armyworm, Spodoptera frugiperda (J.E. Smith), was investigated, at the dose of LC10 and LC90, by transcriptome and LC-MS/MS analysis. Using statistics-based analysis of both POS and NEG mode, a total of 715,501 metabolites in S. frugiperda were significantly changed after spinetoram treatment. The enhancement of glucose metabolism provides energy support for detoxification in larvae. The decrease in valine and isoleucine is associated with an increase in leucine, without maintaining the conservation of citric acid in the larvae. The down-regulation of phenylalanine may retard the tricarboxylic acid cycle to produce GTP. The abundance of lysine was decreased in response to spinetoram exposure, which damages the nervous system of the larvae. The abundance of arginine increases and causes non-functional contraction of the insect's muscles, causing the larva to expend excess energy. Tryptophan provides an important substrate for eliminating ROS. The changes in glutamic acid, aspartic acid, and lysine cause damage to the nerve centers of the larvae. The results of transcriptome and LC-MS/MS analysis revealed the effects of pesticide exposure on amino-acid metabolism of S. frugiperda successfully and provide a new overview of the response of insect physio-biochemistry against pesticides.

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.

Knockdown or inhibition of arginine kinases enhances susceptibility of Tribolium castaneum to deltamethrin

Metabolism of insecticides is an energy-consuming process. As an important component of the intracellular energy buffering system, arginine kinase (AK) plays an important role in insect cellular energy homeostasis and environmental stress response, but the involvement of AKs in the response to chemical stressors (insecticides) remains largely unknown. In this study, using Tribolium castaneum as a model organism, we found that deltamethrin treatment significantly upregulated the expression of TcAK1 and TcAK2 and decreased the whole body ATP content. The knockdown of TcAK1 or TcAK2 significantly enhances the deltamethrin-induced ATP depletion and increase the susceptibility of T. castaneum to deltamethrin. In addition, pretreatment with two AK inhibitors, rutin and quercetin, significantly decreased the lifespan of beetles treated with deltamethrin. These results suggest that AKs might be involved in detoxification of insecticides by regulating cellular energy balance.

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.

Transcriptome Analysis of Aedes albopictus (Diptera: Culicidae) Larvae Exposed With a Sublethal Dose of Haedoxan A

Aedes albopictus is the vector of arbovirus diseases including yellow fever, dengue, Zika virus, and chikungunya fever, and it poses an enormous threat to human health worldwide. Previous studies have revealed that haedoxan A (HA), which is an insecticidal sesquilignan from Phryma leptostachya L., is a highly effective natural insecticide for managing mosquitoes and houseflies; however, the mechanisms underlying the response of Ae. albopictus after treatment with sublethal concentrations of HA is not clear. Here, high-throughput sequencing was used to analyze the gene expression changes in Ae. albopictus larvae after treatment with the LC30 of HA. In total, 416 differentially expressed genes (DEGs) were identified, including 328 upregulated genes and 88 downregulated genes. Identification and verification of related DEGs were performed by RT-qPCR. The results showed that two P450 unigenes (CYP4C21 and CYP304A1), one carboxylesterase, and one ABC transporter (ABCG1) were induced by HA, which indicated that these detoxifying enzyme genes might play a major role in the metabolic and detoxification processes of HA. Additionally, acetylcholine receptor subunit ɑ2 (AChRα2), AChRα5, AChRα9, and the glutamate receptor ionotropic kainate 2 (GRIK2) were found to be upregulated in HA-treated larvae, suggesting that HA affected the conduction of action potentials and synaptic transmission by disrupting the function of neural receptors. These results provide a foundation for further elucidating the target of HA and the mechanism of detoxification metabolism in Ae. albopictus.

Genomic and Transcriptomic Analysis Reveals Cuticular Protein Genes Responding to Different Insecticides in Fall Armyworm Spodoptera frugiperda

The fall armyworm (FAW), Spodoptera frugiperda, is a serious pest of crucial crops causing great threats to the food security of the world. It has evolved resistance to various insecticides, while the underlying molecular mechanisms remain largely unknown. Cuticular proteins (CPs), as primary components in cuticle, play an important role in insects' protection against environmental stresses. Few of them have been documented as participating in insecticide resistance in several insect species. In order to explore whether CP genes of the FAW exhibit a functional role in responding to insecticides stress, a total of 206 CPs, classified into eight families, were identified from the genome of the FAW through a homology-based approach coupled with manual efforts. The temporal expression profiles of all identified CP genes across developmental stages and their responses to 23 different insecticides were analyzed using the RNA-seq data. Expression profiling indicated that most of the CP genes displayed stage-specific expression patterns. It was found that the expression of 51 CP genes significantly changed after 48 h exposure to 17 different insecticides. The expression of eight CP genes responding to four insecticides were confirmed by RT-PCR analysis. The results showed that their overall expression profiles were consistent with RNA-seq analysis. The findings provide a basis for further functional investigation of CPs implied in insecticide stress in FAW.

Resolving the Rules of Robustness and Resilience in Biology Across Scales

Why do some biological systems and communities persist while others fail? Robustness, a system's stability, and resilience, the ability to return to a stable state, are key concepts that span multiple disciplines within and outside the biological sciences. Discovering and applying common rules that govern the robustness and resilience of biological systems is a critical step toward creating solutions for species survival in the face of climate change, as well as the for the ever-increasing need for food, health, and energy for human populations. We propose that network theory provides a framework for universal scalable mathematical models to describe robustness and resilience and the relationship between them, and hypothesize that resilience at lower organization levels contribute to robust systems. Insightful models of biological systems can be generated by quantifying the mechanisms of redundancy, diversity, and connectivity of networks, from biochemical processes to ecosystems. These models provide pathways towards understanding how evolvability can both contribute to and result from robustness and resilience under dynamic conditions. We now have an abundance of data from model and non-model systems and the technological and computational advances for studying complex systems. Several conceptual and policy advances will allow the research community to elucidate the rules of robustness and resilience. Conceptually, a common language and data structure that can be applied across levels of biological organization needs to be developed. Policy advances such as cross-disciplinary funding mechanisms, access to affordable computational capacity, and the integration of network theory and computer science within the standard biological science curriculum will provide the needed research environments. This new understanding of biological systems will allow us to derive ever more useful forecasts of biological behaviors and revolutionize the engineering of biological systems that can survive changing environments or disease, navigate the deepest oceans, or sustain life throughout the solar system.

High-resolution genetic mapping reveals cis-regulatory and copy number variation in loci associated with cytochrome P450-mediated detoxification in a generalist arthropod pest

Chemical control strategies are driving the evolution of pesticide resistance in pest populations. Understanding the genetic mechanisms of these evolutionary processes is of crucial importance to develop sustainable resistance management strategies. The acaricide pyflubumide is one of the most recently developed mitochondrial complex II inhibitors with a new mode of action that specifically targets spider mite pests. In this study, we characterize the molecular basis of pyflubumide resistance in a highly resistant population of the spider mite Tetranychus urticae. Classical genetic crosses indicated that pyflubumide resistance was incompletely recessive and controlled by more than one gene. To identify resistance loci, we crossed the resistant population to a highly susceptible T. urticae inbred strain and propagated resulting populations with and without pyflubumide exposure for multiple generations in an experimental evolution set-up. High-resolution genetic mapping by a bulked segregant analysis approach led to the identification of three quantitative trait loci (QTL) linked to pyflubumide resistance. Two QTLs were found on the first chromosome and centered on the cytochrome P450 CYP392A16 and a cluster of CYP392E6-8 genes. Comparative transcriptomics revealed a consistent overexpression of CYP392A16 and CYP392E8 in the experimental populations that were selected for pyflubumide resistance. We further corroborated the involvement of CYP392A16 in resistance by in vitro functional expression and metabolism studies. Collectively, these experiments uncovered that CYP392A16 N-demethylates the toxic carboxamide form of pyflubumide to a non-toxic compound. A third QTL coincided with cytochrome P450 reductase (CPR), a vital component of cytochrome P450 metabolism. We show here that the resistant population harbors three gene copies of CPR and that this copy number variation is associated with higher mRNA abundance. Together, we provide evidence for detoxification of pyflubumide by cytochrome P450s that is likely synergized by gene amplification of CPR.

Our understanding of the causal genetic variants that drive the evolution of quantitative traits, such as polygenic pesticide resistance, remains very limited. Here, we followed a high-resolution genetic mapping approach to localize the genetic variants that cause pyflubumide resistance in the two-spotted spider mite Tetranychus urticae. Three well-supported QTL were uncovered and pointed towards a major role for cytochrome P450-mediated detoxification. Cis-regulatory variation for cytochrome P450s was observed, and in vitro cytochrome P450 experiments showed that pyflubumide was metabolized into a non-toxic derivate. A third QTL centered on cytochrome P450 reductase (CPR), which is required for cytochrome P450 activity, and is amplified in pyflubumide resistant populations. Our results indicate that pyflubumide resistance is mediated by cytochrome P450 detoxification that is enhanced by gene amplification at the CPR locus.

Knockdown of AMP-activated protein kinase increases the insecticidal efficiency of pymetrozine to Nilaparvata lugens

Insecticides are the main tools used to control Nilaparvata lugens (Stål), a serious pest of rice in Asia. However, repeated application of insecticides has caused many negative effects. Reducing the amount of insecticide used, while maintaining good pest population control, would be valuable. AMP-activated protein kinase (AMPK), a sensor of cellular energy status, helps to maintain insect energy balance at the cellular and whole-body level. The role of AMPK in insect response to insecticide stimulation is unknown. We studied the functions of AMPK catalytic subunit alpha (NlAMPKα) in the development of N. lugens and in response to pymetrozine, an insecticide used to control insect pests with piercing-sucking mouthparts. A phylogenetic analysis of protein sequences from 12 species in six orders showed that insects have only the AMPKα 2 subtype. RNA interference against NlAMPKα demonstrated that blocking the AMPK pathway led to a decrease in the systemic ATP level and an increase in N. lugens mortality. NlAMPKα responded to the energy stress caused by pymetrozine treatment, which activated downstream energy metabolic pathways to compensate for the energy imbalance. However, the ATP level in pymetrozine- treated nymphs was not increased, suggesting that ATP is consumed more than synthesized. When NlAMPKα expression was reduced in pymetrozine-treated nymphs by RNAi, the ATP level was decreased and the mortality was significantly increased. At day eight post 0.5 g/3 L of pymetrozine and dsNlAMPKα treatment, nymph survival was 29.33%, which was similar to the 27.33% survival of 1 g/3 L pymetrozine-treated nymphs. Addition of dsNlAMPKα can reduce the concentration of pymetrozine used by 50% while providing comparable efficacy. These results indicate that AMPK helps maintain the energy metabolism of N. lugens in response to pymetrozine treatment. Knockdown of NlAMPKα increases the insecticidal efficiency of pymetrozine to N. lugens.

The chemical defensome of five model teleost fish

How an organism copes with chemicals is largely determined by the genes and proteins that collectively function to defend against, detoxify and eliminate chemical stressors. This integrative network includes receptors and transcription factors, biotransformation enzymes, transporters, antioxidants, and metal- and heat-responsive genes, and is collectively known as the chemical defensome. Teleost fish is the largest group of vertebrate species and can provide valuable insights into the evolution and functional diversity of defensome genes. We have previously shown that the xenosensing pregnane x receptor (pxr, nr1i2) is lost in many teleost species, including Atlantic cod (Gadus morhua) and three-spined stickleback (Gasterosteus aculeatus), but it is not known if compensatory mechanisms or signaling pathways have evolved in its absence. In this study, we compared the genes comprising the chemical defensome of five fish species that span the teleosteii evolutionary branch often used as model species in toxicological studies and environmental monitoring programs: zebrafish (Danio rerio), medaka (Oryzias latipes), Atlantic killifish (Fundulus heteroclitus), Atlantic cod, and three-spined stickleback. Genome mining revealed evolved differences in the number and composition of defensome genes that can have implication for how these species sense and respond to environmental pollutants, but we did not observe any candidates of compensatory mechanisms or pathways in cod and stickleback in the absence of pxr. The results indicate that knowledge regarding the diversity and function of the defensome will be important for toxicological testing and risk assessment studies.

Transcriptional regulation of xenobiotic detoxification genes in insects - An overview

Arthropods have well adapted to the vast array of chemicals they encounter in their environment. Whether these xenobiotics are plant allelochemicals or anthropogenic insecticides one of the strategies they have developed to defend themselves is the induction of detoxification enzymes. Although upregulation of detoxification enzymes and efflux transporters in response to specific inducers has been well described, in insects, yet, little is known on the transcriptional regulation of these genes. Over the past twenty years, an increasing number of studies with insects have used advanced genetic tools such as RNAi, CRISPR/Cas9 and reporter gene assays to dissect the genomic grounds of their xenobiotic response and hence contributed substantially in improving our knowledge on the players involved. Xenobiotics are partly recognized by various "xenobiotic sensors" such as membrane-bound or nuclear receptors. This initiates a molecular reaction cascade ultimately leading to the translocation of a transcription factor to the nucleus that recognizes and binds to short sequences located upstream their target genes to activate transcription. To date, a number of signaling pathways were shown to mediate the upregulation of detoxification enzymes in arthropods and to play a role in either metabolic resistance to insecticides or host-plant adaptation. These include nuclear receptors AhR/ARNT and HR96, GPCRs, CncC and MAPK/CREB. Recent work reveals that upregulation and activation of some components of these pathways as well as polymorphism in the binding motifs of transcription factors are linked to insects' adaptive processes. The aim of this mini-review is to summarize and describe recent work that shed some light on the main regulatory routes of detoxification gene expression in insects.

Screening of Helicoverpa armigera Mobilome Revealed Transposable Element Insertions in Insecticide Resistance Genes

Simple Summary

Transposable elements (TEs) are mobile DNA sequences that can copy themselves within a host genome. TE-mediated changes in regulation can lead to massive and rapid changes in expression, responses that are potentially highly adaptive when an organism is faced with a mortality agent in the environment, such as an insecticide. Helicoverpa armigera shows a hight number of reported cases of insecticide resistance worldwide, having evolved resistance against pyrethroids, organophosphates, carbamates, organochlorines, and recently to macrocyclic lactone spinosad and several Bacillus thuringiensis toxins. In the present study, we conducted a TE annotation using combined approaches, and the results revealed a total of 8521 TEs, representing 236,132 copies, covering 12.86% of the H. armigera genome. In addition, we underlined TE insertions in defensome genes and we successfully identified nine TE insertions belonging to the RTE, R2, CACTA, Mariner and hAT superfamilies.

Abstract

The cotton bollworm Helicoverpa armigera Hübner (Lepidoptera: Noctuidae) is an important pest of many crops that has developed resistance to almost all groups of insecticides used for its management. Insecticide resistance was often related to Transposable Element (TE) insertions near specific genes. In the present study, we deeply retrieve and annotate TEs in the H. armigera genome using the Pipeline to Retrieve and Annotate Transposable Elements, PiRATE. The results have shown that the TE library consists of 8521 sequences representing 236,132 TE copies, including 3133 Full-Length Copies (FLC), covering 12.86% of the H. armigera genome. These TEs were classified as 46.71% Class I and 53.29% Class II elements. Among Class I elements, Short and Long Interspersed Nuclear Elements (SINEs and LINEs) are the main families, representing 21.13% and 19.49% of the total TEs, respectively. Long Terminal Repeat (LTR) and Dictyostelium transposable element (DIRS) are less represented, with 5.55% and 0.53%, respectively. Class II elements are mainly Miniature Inverted Transposable Elements (MITEs) (49.11%), then Terminal Inverted Repeats (TIRs) (4.09%). Superfamilies of Class II elements, i.e., Transib, P elements, CACTA, Mutator, PIF-harbinger, Helitron, Maverick, Crypton and Merlin, were less represented, accounting for only 1.96% of total TEs. In addition, we highlighted TE insertions in insecticide resistance genes and we successfully identified nine TE insertions belonging to RTE, R2, CACTA, Mariner and hAT superfamilies. These insertions are hosted in genes encoding cytochrome P450 (CyP450), glutathione S-transferase (GST), and ATP-binding cassette (ABC) transporter belonging to the G and C1 family members. These insertions could therefore be involved in insecticide resistance observed in this pest.

Reduced stress defence responses contribute to the higher toxicity of a pesticide under warming

There is a pressing need to identify the molecular mechanisms underlying the, often magnifying, interactive effects between contaminants and natural stressors. Here we test our hypothesis that lower general stress defence responses contribute to synergistic interactions between stressors. We focus on the widespread pattern that many contaminants are more toxic at higher temperatures. Specifically, we tested the effects of an environmentally realistic low-effect and high-effect concentration of the pesticide chlorpyrifos under warming at the gene expression level in the northern house mosquito Culex pipiens molestus (Forskal, 1775). By applying the independent action model for combined stressors on RNA-sequencing data, we identified interactive gene expression patterns under combined exposure to chlorpyrifos and warming for general stress defence responses: protection of macromolecules, antioxidant processes, detoxification and energy metabolism/allocation. Most of these general stress defence response genes showed upregulated antagonistic interactions (i.e., were less upregulated than expected under the independent action model). This indicates that when pesticide exposure was combined with warming, the general stress defence responses were no longer buffering increased stress levels, which may contribute to a higher sensitivity to toxicants under warming. These upregulated antagonistic interactions were stronger for the high-effect chlorpyrifos concentration, indicating that exposure to this concentration under warming was most stressful. Our results highlight that quantitative analysis of the frequency and strength of the interaction types of general stress defence response genes, specifically focusing on antagonistic upregulations and synergistic downregulations, may advance our understanding of how natural stressors modify the toxicity of contaminants.

MAPK-directed activation of the whitefly transcription factor CREB leads to P450-mediated imidacloprid resistance

The evolution of insect resistance to pesticides poses a continuing threat to agriculture and human health. While much is known about the proximate molecular and biochemical mechanisms that confer resistance, far less is known about the regulation of the specific genes/gene families involved, particularly by trans-acting factors such as signal-regulated transcription factors. Here we resolve in fine detail the trans-regulation of CYP6CM1, a cytochrome P450 that confers resistance to neonicotinoid insecticides in the whitefly Bemisia tabaci, by the mitogen-activated protein kinase (MAPK)-directed activation of the transcription factor cAMP-response element binding protein (CREB). Reporter gene assays were used to identify the putative promoter of CYP6CM1, but no consistent polymorphisms were observed in the promoter of a resistant strain of B. tabaci (imidacloprid-resistant, IMR), which overexpresses this gene, compared to a susceptible strain (imidacloprid-susceptible, IMS). Investigation of potential trans-acting factors using in vitro and in vivo assays demonstrated that the bZIP transcription factor CREB directly regulates CYP6CM1 expression by binding to a cAMP-response element (CRE)-like site in the promoter of this gene. CREB is overexpressed in the IMR strain, and inhibitor, luciferase, and RNA interference assays revealed that a signaling pathway of MAPKs mediates the activation of CREB, and thus the increased expression of CYP6CM1, by phosphorylation-mediated signal transduction. Collectively, these results provide mechanistic insights into the regulation of xenobiotic responses in insects and implicate both the MAPK-signaling pathway and a transcription factor in the development of pesticide resistance.

Identification of transcriptome and fluralaner responsive genes in the common cutworm Spodoptera litura Fabricius, based on RNA-seq

Background

Fluralaner is a novel isoxazoline insecticide with a unique action site on the γ-aminobutyric acid receptor (GABAR), shows excellent activity on agricultural pests including the common cutworm Spodoptera litura, and significantly influences the development and fecundity of S. litura at either lethal or sublethal doses. Herein, Illumina HiSeq Xten (IHX) platform was used to explore the transcriptome of S. litura and to identify genes responding to fluralaner exposure.

Results

A total of 16,572 genes, including 451 newly identified genes, were observed in the S. litura transcriptome and annotated according to the COG, GO, KEGG and NR databases. These genes included 156 detoxification enzyme genes [107 cytochrome P450 enzymes (P450s), 30 glutathione S-transferases (GSTs) and 19 carboxylesterases (CarEs)] and 24 insecticide-targeted genes [5 ionotropic GABARs, 1 glutamate-gated chloride channel (GluCl), 2 voltage-gated sodium channels (VGSCs), 13 nicotinic acetylcholine receptors (nAChRs), 2 acetylcholinesterases (AChEs) and 1 ryanodine receptor (RyR)]. There were 3275 and 2491 differentially expressed genes (DEGs) in S. litura treated with LC₃₀ or LC₅₀ concentrations of fluralaner, respectively. Among the DEGs, 20 related to detoxification [16 P450s, 1 GST and 3 CarEs] and 5 were growth-related genes (1 chitin and 4 juvenile hormone synthesis genes). For 26 randomly selected DEGs, real-time quantitative PCR (RT-qPCR) results showed that the relative expression levels of genes encoding several P450s, GSTs, heat shock protein (HSP) 68, vacuolar protein sorting-associated protein 13 (VPSAP13), sodium-coupled monocarboxylate transporter 1 (SCMT1), pupal cuticle protein (PCP), protein takeout (PT) and low density lipoprotein receptor adapter protein 1-B (LDLRAP1-B) were significantly up-regulated. Conversely, genes encoding esterase, sulfotransferase 1C4, proton-coupled folate transporter, chitinase 10, gelsolin-related protein of 125 kDa (GRP), fibroin heavy chain (FHC), fatty acid synthase and some P450s were significantly down-regulated in response to fluralaner.

Conclusions

The transcriptome in this study provides more effective resources for the further study of S. litura whilst the DEGs identified sheds further light on the molecular response to fluralaner.

Age- and sex-related ABC transporter expression in pyrethroid-susceptible and -resistant Aedes aegypti

Resistance mechanisms to synthetic insecticides often include point mutations and increased expression of genes encoding detoxification enzymes. Since pyrethroids are the main adulticides used against Aedes aegypti, which vectors pathogens such as Zika virus, understanding resistance to this insecticide class is of significant relevance. We focused on adenosine triphosphate (ATP)-binding cassette (ABC) transporters in the pyrethroid-resistant Puerto Rico (PR) strain of Ae. aegypti. We investigated the expression patterns of six ABC transporters previously characterized as differentially expressed in insecticide-challenged mosquitoes, or increased mRNA expression in pyrethroid-resistant Ae. aegypti, by comparing PR to the Rockefeller (Rock) susceptible strain. No constitutive differential expression between strains was detected, but expression differences for these genes was influenced by sex and age, suggesting that their role is independent from resistance in PR. Instead, ABC transporters may be induced after insecticide exposure. Challenging mosquitoes with deltamethrin, with or without ABC transporter modulators, showed that Rock and PR responded differently, but a contribution of ABC transporters to deltamethrin toxicity is suspected. Moreover, the effect of dexamethasone, which enhanced the inhibition of nerve firing by deltamethrin, was observed using a Drosophila central nervous system preparation, showing synergy of these two compounds through the potential inhibition of ABC transporters.

Gene silencing through RNAi and antisense Vivo-Morpholino increases the efficacy of pyrethroids on larvae of Anopheles stephensi

Background

Insecticides are still at the core of insect pest and vector control programmes. Several lines of evidence indicate that ABC transporters are involved in detoxification processes against insecticides, including permethrin and other pyrethroids. In particular, the ABCG4 gene, a member of the G subfamily, has consistently been shown to be up-regulated in response to insecticide treatments in the mosquito malaria vector Anopheles stephensi (both adults and larvae).

Methods

To verify the actual involvement of this transmembrane protein in the detoxification process of permethrin, bioassays on larvae of An. stephensi, combining the insecticide with a siRNA, specifically designed for the inhibition of ABCG4 gene expression were performed. Administration to larvae of the same siRNA, labeled with a fluorescent molecule, was effected to investigate the systemic distribution of the inhibitory RNA into the larval bodies. Based on siRNA results, similar experiments using antisense Vivo-Morpholinos (Vivo-MOs) were effected. These molecules, compared to siRNA, are expected to guarantee a higher stability in environmental conditions and in the insect gut, and present thus a higher potential for future in-field applications.

Results

Bioassays using two different concentrations of siRNA, associated with permethrin, led to an increase of larval mortality, compared with results with permethrin alone. These outcomes confirm that ABCG4 transporter plays a role in the detoxification process against the selected insecticide. Moreover, after fluorescent labelling, it was shown the systemic dissemination of siRNA in different body districts of An. stephensi larvae, which suggest a potential systemic effect of the molecule. At the same time, results of Vivo-MO experiments were congruent with those obtained using siRNA, thus confirming the potential of ABCG4 inhibition as a strategy to increase permethrin susceptibility in mosquitoes. For the first time, Vivo-MOs were administered in water to larvae, with evidence for a biological effect.

Conclusions

Targeting ABCG4 gene for silencing through both techniques resulted in an increased pyrethroid efficacy. These results open the way toward the possibility to exploit ABCG4 inhibition in the context of integrated programmes for the control An. stephensi mosquitoes and malaria transmission.

HSP superfamily of genes in the malaria vector Anopheles sinensis: diversity, phylogenetics and association with pyrethroid resistance

BACKGROUND

Heat shock proteins (HSPs) are molecular chaperones that are involved in many normal cellular processes and various kinds of environmental stress. There is still no report regarding the diversity and phylogenetics research of HSP superfamily of genes at whole genome level in insects, and the HSP gene association with pyrethroid resistance is also not well known. The present study investigated the diversity, classification, scaffold location, characteristics, and phylogenetics of the superfamily of genes in Anopheles sinensis genome, and the HSP genes associated with pyrethroid resistance.

METHODS

The present study identified the HSP genes in the An. sinensis genome, analysed their characteristics, and deduced phylogenetic relationships of all HSPs in An. sinensis, Anopheles gambiae, Culex quinquefasciatus and Aedes aegypti by bioinformatic methods. Importantly, the present study screened the HSPs associated with pyrethroid resistance using three field pyrethroid-resistant populations with RNA-seq and RT-qPCR, and looked over the HSP gene expression pattern for the first time in An. sinensis on the time-scale post insecticide treatment with RT-qPCR.

RESULTS

There are 72 HSP genes in An. sinensis genome, and they are classified into five families and 11 subfamilies based on their molecular weight, homology and phylogenetics. Both RNA-seq and qPCR analysis revealed that the expression of AsHSP90AB, AsHSP70-2 and AsHSP21.7 are significantly upregulated in at least one field pyrethroid-resistant population. Eleven genes are significantly upregulated in different period after pyrethroid exposure. The HSP90, sHSP and HSP70 families are proposed to be involved in pyrethroid stress response based in expression analyses of three field pyrethroid-resistant populations, and expression pattern on the time scale post insecticide treatment. The AsHSP90AB gene is proposed to be the essential HSP gene for pyrethroid stress response in An. sinensis.

CONCLUSIONS

This study provides the information frame for HSP superfamily of genes, and lays an important basis for the better understanding and further research of HSP function in insect adaptability to diverse environments.

Insecticide Exposure Triggers a Modulated Expression of ABC Transporter Genes in Larvae of Anopheles gambiae s.s

Insecticides remain a main tool for the control of arthropod vectors. The urgency to prevent the insurgence of insecticide resistance and the perspective to find new target sites, for the development of novel molecules, are fuelling the study of the molecular mechanisms involved in insect defence against xenobiotic compounds. In this study, we have investigated if ATP-binding cassette (ABC) transporters, a major component of the defensome machinery, are involved in defence against the insecticide permethrin, in susceptible larvae of the malaria vector Anopheles gambiae sensu stricto. Bioassays were performed with permethrin alone, or in combination with an ABC transporter inhibitor. Then we have investigated the expression profiles of five ABC transporter genes at different time points following permethrin exposure, to assess their expression patterns across time. The inhibition of ABC transporters increased the larval mortality by about 15-fold. Likewise, three genes were up-regulated after exposure to permethrin, showing different patterns of expression across the 48 h. Our results provide the first evidences of ABC transporters involvement in defence against a toxic in larvae of An. gambiae s.s. and show that the gene expression response is modulated across time, being continuous, but stronger at the earliest and latest times after exposure.

Transcriptional Responses in Defense-Related Genes of Sitobion avenae (Hemiptera: Aphididae) Feeding on Wheat and Barley

Many aphids can adapt on plants of differential resistance levels, but molecular interactions underlying host plant utilization and shift of aphids are still not well understood. Here, we sequenced the transcriptome and compared global gene expression profiles of Sitobion avenae (Fabricius) feeding on wheat (i.e., the susceptible plant) and barley (i.e., the resistant plant). In total, 65,200 high-quality unigenes were identified from the merged transcriptomes, and 861 of them were differentially expressed. Relative to their expression on barley, all differentially expressed P450 (cytochrome 450 monooxygenase) genes (e.g., SavCYP6A13, SavCYP4C1, and SavCYP4G15) in S. avenae on wheat were upregulated, suggesting the significance of P450s in xenobiotic metabolism for this aphid on wheat. For S. avenae on barley, some genes encoding for ATP-binding cassette (ABC) transporters (e.g., ABCG1, ABCG4, ABCB7, and ABCA5) and UDP-glucuronosyltransferases (UGTs; e.g., UGT2B1 and UGT2C1) showed a dramatic increase in expression, suggesting that ABC transporters and UGTs could be critical for detoxification metabolism in S. avenae on barley. In addition, the expression for nearly all differentially expressed proteases was upregulated for S. avenae on barley, reflecting xenobiotic challenges facing S. avenae on resistant barley. Thus, various detoxification and other defense-related genes (e.g., proteases and oxidases) showed highly inducible transcript regulation, suggesting great adaptive potential for S. avenae on different plants. All the abovementioned genes will be prime candidates for further studies of molecular interactions underlying host plant use and specialization in this aphid. Our results provide insights into aphids' defenses against variable phytochemicals, and the molecular basis underlying aphid-plant interactions as well.

Adenosine Triphosphate-Binding Cassette Transporters Are Not Involved In the Detoxification of Azadirachta indica Extracts In Anopheles stephensi Larvae

Detoxifying pathways of mosquitoes against the neem (Azadirachta indica) extracts are still unclear. The aim of the present study was to investigate the role of adenosine triphosphate-binding cassette (ABC) transporters in this process in Anopheles stephensi, one of the main malaria vectors in southern Asia. Third-stage larvae of An. stephensi were fed with fish food alone or in combination with neem extract at 0.5%, 1%, 5%, and 10%. Six ABC-transporter genes from 3 different subfamilies (B, C, and G) were analyzed to assess their relative expression compared with controls. A bioassay was also performed to assess larval mortality rate at different concentrations and in combination with verapamil, an ABC-transporter inhibitor. No significant variation in the expression levels of any transporter belonging to the B, C, and G subfamilies was detected. Furthermore, the use of verapamil did not induce an increase in mortality at any of the tested neem extract concentrations, indicating that ABC transporters are not involved in the detoxification of neem extracts in An. stephensi larvae.

Regulating resistance: CncC:Maf, antioxidant response elements and the overexpression of detoxification genes in insecticide resistance

Although genetic and genomic tools have greatly furthered our understanding of resistance-associated mutations in molecular target sites of insecticides, the genomic basis of transcriptional regulation of detoxification loci in insect pests and vectors remains relatively unexplored. Recent work using RNAi, reporter assays and comparative genomics are beginning to reveal the molecular architecture of this response, identifying critical transcription factors and their binding sites. Central to this is the insect ortholog of the mammalian transcription factor Nrf2, Cap 'n' Collar isoform-C (CncC) which as a heterodimer with Maf-S regulates the transcription of phase I, II and III detoxification loci in a range of insects, with CncC knockdown or upregulation directly affecting phenotypic resistance. CncC:Maf binds to specific antioxidant response element sequences upstream of detoxification genes to initiate transcription. Recent work is now identifying these binding sites for resistance-associated loci and, coupled with genome sequence data and reporter assays, enabling identification of polymorphisms in the CncC:Maf binding site which regulate the insecticide resistance phenotype.

Does host plant adaptation lead to pesticide resistance in generalist herbivores?

Most herbivorous arthropods feed on one or a few closely related plant species; however, certain insect and mite species have a greatly expanded host range. Several of these generalists also show a remarkable propensity to evolve resistance to chemical pesticides. In this review, we ask if the evolution of mechanisms to tolerate the diversity of plant secondary metabolites that generalist herbivores encounter, has pre-adapted them to resist synthetic pesticides. Critical examination of the evidence suggests that a generalist life-style per se is not a predictor of rapid resistance evolution to pesticides. Rather the prevalence of pesticide resistance in generalist herbivores probably reflects their economic importance as pests and thus the strong selection imposed by intensive pesticide use.

Gene expression modulation of ABC transporter genes in response to permethrin in adults of the mosquito malaria vector Anopheles stephensi

Living organisms have evolved an array of genes coding for detoxifying enzymes and efflux protein pumps, to cope with endogenous and xenobiotic toxic compounds. The study of the genes activated during toxic exposure is relevant to the area of arthropod vector control, since these genes are one of the targets upon which natural selection acts for the evolution of insecticide resistance. ATP-binding cassette (ABC) transporters participate to insecticide detoxification acting as efflux pumps, that reduce the intracellular concentration of toxic compounds, or of their metabolic derivatives. Here we analyzed the modulation of the expression of six genes coding for ABC transporters, after the exposure of adult females and males of the mosquito Anopheles stephensi, a major malaria vector in Asia, to permethrin. Male and female mosquitoes were exposed to insecticide for one hour, then the expression profiles of the ABC transporter genes AnstABCB2, AnstABCB3, AnstABCB4, AnstABCBmember6, AnstABCC11, and AnstABCG4 were analysed after one and 24h. Our results showed that three genes (AnstABCB2, AnstABCBmember6, AnstABCG4) were up-regulated in both sexes; two of these (AnstABCBmember6 and AnstABCG4) have previously been shown to be up-regulated also in larval stages of An. stephensi, supporting a role for these genes in permethrin defence in larvae as well as in adults. Finally, the same ABC transporter genes were activated both in females and males; however, the timing of gene induction was different, with a prompter induction in females than in males.