A comparison of Drosophila melanogaster detoxification gene induction responses for six insecticides, caffeine and phenobarbital

Pb exposure causes non-linear accumulation of Pb in D. melanogaster controlled by metallothionein B and exerts ecological effects

Pb pollution can harm human health and the ecosystem. Therefore, it is worthwhile to study the metabolic processes of heavy metals in individual bodies and their influence on ecological systems. In this work, we analyzed the genetic responses and physiological changes of D. melanogaster which took diets exposed to different doses of Pb using transcriptomic analysis, ICP-MS, and various other physiological methods. We found that the Pb accumulated in D. melanogaster in a nonlinear pattern with the increase of Pb content in food. Metallothioneins (Mtns), especially the MtnB directly affects the accumulation and excretion of metal Pb in D. melanogaster, and causes the nonlinear accumulation. Metal regulatory transcription factor-1 (MTF-1) is involved in the regulation of Pb-induced high expressions of Mtns. Furthermore, an interaction between the metal metabolism pathway and xenobiotic response pathway leads to the cross-tolerances of Pb-exposed D. melanogaster to insecticides and other toxins. The oxidative stress induced by Pb toxicity may be the bridge between them. Our findings provide a physiological and molecular genetic basis for further study of the accumulation and metabolism of Pb in D. melanogaster.

Perspectives on the Drosophila melanogaster Model for Advances in Toxicological Science

The use of Drosophila melanogaster for studies of toxicology has grown considerably in the last decade. The Drosophila model has long been appreciated as a versatile and powerful model for developmental biology and genetics because of its ease of handling, short life cycle, low cost of maintenance, molecular genetic accessibility, and availability of a wide range of publicly available strains and data resources. These features, together with recent unique developments in genomics and metabolomics, make the fly model especially relevant and timely for the development of new approach methodologies and movements toward precision toxicology. Here, we offer a perspective on how flies can be leveraged to identify risk factors relevant to environmental exposures and human health. First, we review and discuss fundamental toxicologic principles for experimental design with Drosophila. Next, we describe quantitative and systems genetics approaches to resolve the genetic architecture and candidate pathways controlling susceptibility to toxicants. Finally, we summarize the current state and future promise of the emerging field of Drosophila metabolomics for elaborating toxic mechanisms. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC.

The transcriptomic signature of responses to larval crowding in Drosophila melanogaster

Intraspecific competition at the larval stage is an important ecological factor affecting life-history, adaptation and evolutionary trajectory in holometabolous insects. However, the molecular pathways underpinning these ecological processes are poorly characterized. We reared Drosophila melanogaster at three egg densities (5, 60, and 300 eggs/mL) and sequenced the transcriptomes of pooled third-instar larvae. We also examined emergence time, egg-to-adult viability, adult mass, and adult sex-ratio at each density. Medium crowding had minor detrimental effects on adult phenotypes compared to low density and yielded 24 differentially expressed genes (DEGs), including several chitinase enzymes. In contrast, high crowding had substantial detrimental effects on adult phenotypes and yielded 2107 DEGs. Among these, upregulated gene sets were enriched in sugar, steroid and amino acid metabolism as well as DNA replication pathways, whereas downregulated gene sets were enriched in ABC transporters, taurine, Toll/Imd signaling, and P450 xenobiotics metabolism pathways. Overall, our findings show that larval crowding has a large consistent effect on several molecular pathways (i.e., core responses) with few pathways displaying density-specific regulation (i.e., idiosyncratic responses). This provides important insights into how holometabolous insects respond to intraspecific competition during development.

The dynamic changes of genes revealed that persistently overexpressed genes drive the evolution of cyflumetofen resistance in Tetranychus cinnabarinus

Changes in gene expression are associated with the evolution of pesticide resistance in arthropods. In this study, transcriptome sequencing was performed in 3 different resistance levels (low, L; medium, M; and high, H) of cyflumetofen-resistant strain (YN-CyR). A total of 1 685 genes, including 97 detoxification enzyme genes, were upregulated in all 3 stages, of which 192 genes, including 11 detoxification enzyme genes, showed a continuous increase in expression level with resistance development (L to H). RNA interference experiments showed that overexpression of 7 genes (CYP392A1, TcGSTd05, CCE06, CYP389A1, TcGSTz01, CCE59, and CYP389C2) is involved in the development of cyflumetofen resistance in Tetranychus cinnabarinus. The recombinant CYP392A1 can effectively metabolize cyflumetofen, while CCE06 can bind and sequester cyflumetofen in vitro. We compared 2 methods for rapid screening of resistance molecular markers, including short-term induction and 1-time high-dose selection. Two detoxification enzyme genes were upregulated in the field susceptible strain (YN-S) by induction with 20% lethal concentration (LC₂₀ ) of cyflumetofen. However, 16 detoxification enzyme genes were upregulated by 1-time selection with LC₈₀ of cyflumetofen. Interestingly, the 16 genes were overexpressed in all 3 resistance stages. These results indicated that 1 685 genes that were upregulated at the L stage constituted the basis of cyflumetofen resistance, of which 192 genes in which upregulation continued to increase were the main driving force for the development of resistance. Moreover, the 1-time high-dose selection is an efficient way to rapidly obtain the resistance-related genes that can aid in the development of resistance markers and resistance management in mites.

Functional characterization of ABC transporters mediates multiple neonicotinoid resistance in a field population of Aphis gossypii Glover

The ATP-binding cassette (ABC) transporters C and G subfamilies have been reported to be involved in insecticide detoxification, with most studies showing increased gene transcript levels in response to insecticide exposure. Our previous studies have suggested that ABCC and G transporters participate in cyantraniliprole and thiamethoxam resistance of Aphis gossypii. In this study, we focused on the potential roles of the ABCC and G transporters of an A. gossypii field population (SDR) in neonicotinoid detoxification. The results of leaf dip bioassays showed 629.17- and 346.82-fold greater resistance to thiamethoxam and imidacloprid in the SDR strain, respectively, than in the susceptible strain (SS). Verapamil, an ABC inhibitor, was used for synergism bioassays, and the results showed synergistic effects with thiamethoxam, with synergistic ratios (SRs) of 2.07 and 6.68 in the SS and SDR strains, respectively. In addition to thiamethoxam, verapamil increased imidacloprid toxicity by 1.68- and 1.62-fold in the SS and SDR strains respectively. Then, the expression levels of several ABCC and G transporters were analyzed in different treatments. We found that the transcript levels of AgABCG4, AgABCG17, AgABCG26, AgMRP8 and AgMRP12 were higher in the SDR strain than in the SS strain. The mRNA expression of AgABCG4, AgABCG7, AgABCG13, AgABCG17, AgABCG26, AgMRP8 and AgMRP12 in the SDR strain was increased after thiamethoxam and imidacloprid exposure. The results of transgenic Drosophila melanogaster bioassays suggested that overexpression of AgABCG4, AgABCG7, AgABCG13, AgABCG17, AgABCG26, AgMRP8 and AgMRP12 in transgenic flies was sufficient to confer thiamethoxam and imidacloprid resistance, and AgABCG4, AgABCG7, AgABCG13, AgABCG26 and AgMRP12 may be related to α-cypermethrin cross-resistance with weak effects. In addition, the knockdown of AgABCG4, AgABCG13, AgABCG26, AgMRP8 and AgMRP12, and the knockdown of AgABCG7 and AgABCG26 increased thiamethoxam and imidacloprid mortality in the SDR strain, respectively. Our results suggest that changes in the expression levels of ABCC and G transporters may contribute to neonicotinoid detoxification in the SDR strain, and provide a foundation for clarify the potential roles of ABCC and G transporters in insecticide resistance.

Resistance to dicyclanil and imidacloprid in the sheep blowfly, Lucilia cuprina, in Australia

BACKGROUND

The sheep blowfly, Lucila cuprina, is a myiasis-causing parasite responsible for significant production losses and welfare issues for the Australian sheep industry. Control relies largely on the use of insecticides. The pyrimidine compound, dicyclanil, is the predominant control chemical, although other insecticides also are used, including imidacloprid, ivermectin, cyromazine and spinosad. We investigated in vitro resistance patterns and mechanisms in field-collected blowfly strains.

RESULTS

The Walgett 2019 strain showed significant levels of resistance to both dicyclanil and imidacloprid, with resistance factors at the IC₅₀ of 26- and 17-fold, respectively, in in vitro bioassays. Co-treatment with the cytochrome P450 inhibitor, aminobenzotriazole, resulted in significant levels of synergism for dicyclanil and imidacloprid (synergism ratios of 7.2- and 6.1-fold, respectively), implicating cytochrome P450 in resistance to both insecticides. Cyp12d1 transcription levels were increased up to 40-fold throughout the larval life stages in the resistant strain compared to a reference susceptible strain, whereas transcription levels of some other cyp genes (6g1, 4d1, 28d1) did not differ between the strains. Similar resistance levels also were observed in flies collected from the same property in two subsequent years.

CONCLUSION

This study indicates that in vitro resistance to both dicyclanil and imidacloprid in this field-collected blowfly strain is likely mediated by cytochrome P450, with Cyp12d1 implicated as the enzyme responsible; however, it remains possible that another P450 also may be involved. A common resistance mechanism for the two drugs has important implications for drug rotation strategies designed to prolong the useful life of flystrike control chemicals. © 2022 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Expression Patterns of Drosophila Melanogaster Glutathione Transferases

Glutathione transferases (GSTs) are ubiquitous enzymes that catalyze the conjugation of glutathione to various molecules. Among the 42 GSTs identified in Drosophila melanogaster, Delta and Epsilon are the largest classes, with 25 members. The Delta and Epsilon classes are involved in different functions, such as insecticide resistance and ecdysone biosynthesis. The insect GST number variability is due mainly to these classes. Thus, they are generally considered supports during the evolution for the adaptability of the insect species. To explore the link between Delta and Epsilon GST and their evolution, we analyzed the sequences using bioinformatic tools. Subgroups appear within the Delta and Epsilon GSTs with different levels of diversification. The diversification also appears in the sequences showing differences in the active site. Additionally, amino acids essential for structural stability or dimerization appear conserved in all GSTs. Quantitative real-time polymerase chain reaction (qRT-PCR) analysis revealed that the transcripts corresponding to these two classes are heterogeneously expressed within D. melanogaster. Some GSTs, such as GSTD1, are highly expressed in all tissues, suggesting their general function in detoxification. Conversely, some others, such as GSTD11 or GSTE4, are specifically expressed at a high level specifically in antennae, suggesting a potential role in olfaction.

Powerful, efficient QTL mapping in Drosophila melanogaster using bulked phenotyping and pooled sequencing

Despite the value of recombinant inbred lines for the dissection of complex traits, large panels can be difficult to maintain, distribute, and phenotype. An attractive alternative to recombinant inbred lines for many traits leverages selecting phenotypically extreme individuals from a segregating population, and subjecting pools of selected and control individuals to sequencing. Under a bulked or extreme segregant analysis paradigm, genomic regions contributing to trait variation are revealed as frequency differences between pools. Here, we describe such an extreme quantitative trait locus, or extreme quantitative trait loci, mapping strategy that builds on an existing multiparental population, the Drosophila Synthetic Population Resource, and involves phenotyping and genotyping a population derived by mixing hundreds of Drosophila Synthetic Population Resource recombinant inbred lines. Simulations demonstrate that challenging, yet experimentally tractable extreme quantitative trait loci designs (≥4 replicates, ≥5,000 individuals/replicate, and selecting the 5-10% most extreme animals) yield at least the same power as traditional recombinant inbred line-based quantitative trait loci mapping and can localize variants with sub-centimorgan resolution. We empirically demonstrate the effectiveness of the approach using a 4-fold replicated extreme quantitative trait loci experiment that identifies 7 quantitative trait loci for caffeine resistance. Two mapped extreme quantitative trait loci factors replicate loci previously identified in recombinant inbred lines, 6/7 are associated with excellent candidate genes, and RNAi knock-downs support the involvement of 4 genes in the genetic control of trait variation. For many traits of interest to drosophilists, a bulked phenotyping/genotyping extreme quantitative trait loci design has considerable advantages.

Xenobiotic responses in insects

Insects have evolved a powerful detoxification system to protect themselves against environmental and anthropogenic xenobiotics including pesticides and nanoparticles. The resulting tolerance to insecticides is an immense problem in agriculture. In this study, we summarize advances in our understanding of insect xenobiotic responses: the detoxification strategies and the regulation mechanisms against xenobiotics including nanoparticles, the problem of response specificity and the potential usefulness of this study field for an elaborate pest management. In particular, we highlight that versatility of the detoxification system relies on the relatively unspecific recognition of a broad range of potential toxic substances that trigger either of various canonical xenobiotic responses signaling pathways, including CncC/Keap1, HR96, AHR/ARNT, GPCR, and MAPK/CREB. However, it has emerged that the actual response to an inducer may nevertheless be specific. There are two nonexclusive possibilities that may explain response specificity: (1) differential cross-talk between the known pathways and (2) additional, yet unidentified regulators and pathways of detoxification. Hence, a deeper and broader understanding of the regulation mechanisms of xenobiotic response in insects in the future might facilitate the development and application of highly efficient and environmentally friendly pest control methods, allowing us to face the challenge of the world population growth.

Evaluation of tolerance to λ-cyhalothrin and response of detoxification enzymes in silkworms reared on artificial diet

A representative silkworm rearing mode of Ⅰ-Ⅲ instars reared on artificial diet and Ⅳ-Ⅴ instars reared on fresh mulberry leaves has been recognized in some sericultural areas of China. Under this rearing mode, silkworms are prone to be poisoned by pesticide residues on mulberry leaves at the Ⅳ and Ⅴ instar stages. As one of the most widely applied insecticides, λ-cyhalothrin was used to study the insecticide tolerance of silkworm reared on artificial diet (referred as the AD group). Our results showed that the newly ecdysized Ⅳ instar larvae in the AD group were less tolerant to λ-cyhalothrin compared to the mulberry leaves reared group (referred as the ML group). After continuous exposure to trace λ-cyhalothrin, the weight gain and the survival rate of silkworms were significantly lower in the AD group than those in the ML group, even though compensatory growth was observed in the control of the AD group. Histopathology and ultrastructure of fat body showed that λ-cyhalothrin induced more severe cell injuries in the AD group, such as shrunken nucleus, dilatation of endoplasmic reticulum, and mitochondrial swelling. The transcription levels of detoxification related genes (CYP4M5, CYP6AB4, CarE2, CarE5, GSTe1 and GSTe3) and the enzyme activities of P450s, CarEs and GSTs were inducible by trace λ-cyhalothrin in a time-specific manner, and the data showed that the response of P450 enzyme activity was retarded in the AD group, indicating a potential reason for a higher sensitivity to λ-cyhalothrin. Our results provided a new clue for the study of the relationship between feed nutrition and detoxification ability, and also provided an important reference for the development of modern silkworm rearing mode.

Short exposure to nitenpyram pesticide induces effects on reproduction, development and metabolic gene expression profiles in Drosophila melanogaster (Diptera: Drosophilidae)

Although the toxicity of neonicotinoid insecticides has been demonstrated in several studies, the information on metabolism, behavior, and health risk remains limited and has raised concerns about its potential toxicity. Thus, in this study we assessed the effects of nitenpyram using different sublethal concentrations (one-third and one-tenth of the acute LC₅₀ values) on various developmental and metabolic parameters from gene expression regulation in Drosophila melanogaster (model system used worldwide in ecotoxicological studies). As a result, nitenpyram sublethal concentrations prolonged the developmental time for both pupation and eclosion. Additionally, nitenpyram sublethal concentrations significantly decreased the lifespan, pupation rate, eclosion rate, and production of eggs of D. melanogaster. Moreover, the mRNA expression of genes relevant for development and metabolism was significantly elevated after exposure. Mixed function oxidase enzymes (Cyp12d1), (Cyp9f2), and (Cyp4ae1), hemocyte proliferation (RyR), and immune response (IM4) genes were upregulated, whereas lifespan (Atg7), male mating behavior (Ple), female fertility (Ddc), and lipid metabolism (Sxe2) genes were downregulated. These findings support a solid basis for further research to determine the hazardous effects of nitenpyram on health and the environment.

Transcriptomic analysis of salivary gland and proteomic analysis of oral secretion in Helicoverpa armigera under cotton plant leaves, gossypol, and tannin stresses

Gossypol and tannin are involved in important chemical defense processes in cotton plants. In this study, we used transcriptomics and proteomics to explore the changes in salivary gland functional genes and oral secretion (OS) proteins after feeding with artificial diet (containing gossypols and tannins) and cotton plant leaves. We found that dietary cotton plant leaves, gossypols and tannins exerted adverse impacts on the genes that regulated the functions of peptidase, GTPase, glycosyl hydrolases in the salivary glands of the Helicoverpa armigera (H. armigera). However, GST, UGT, hydrolases, and lipase genes were up-regulated to participate in the detoxification and digestive of H. armigera. The oral secretory proteins of H. armigera were significantly inhibited under the stress of gossypol and tannin, such as enzyme activity, but some proteins (such as PZC71358.1) were up-regulated and involved in immune and digestive functions. The combined analysis of transcriptomics and metabolomics showed a weak correlation, and the genes and proteins involved were mainly in digestive enzyme activities. Our work clarifies the deleterious physiological impacts of gossypols and tannins on H. armigera and reveals the mechanism by which H. armigera effectively mitigate the phytotoxic effects through detoxification and immune systems.

Short-term metabolic resistance inductive effect of different agrochemical groups on Anopheles gambiae mosquitoes

Background: In order to assess the impact of the different groups of agricultural pesticides used in Côte d'Ivoire on the increase of mosquitoes resistance to insecticides, the expression profiles of 7 P450 cytochromes and one GSTE2 of Anopheles gambiae involved in mosquito resistance to insecticides were studied. The goal of this study was to determine the effect of short exposure of mosquito larvae to different groups of agricultural pesticides on mosquito resistance. Methods: Three groups of pesticides were selected: (i) agricultural insecticide solutions, (ii) non-insecticide pesticide solutions (a mixture of herbicides and fungicides), and (iii) a mixture of the first two. A fourth non-pesticide solution was used as a control. Four groups of each stage 2 larvae (strain Kisumu, male and female) were exposed to 20% concentrated solution for 24 hours. Susceptibility tests for dichlorodiphenyltrichloroethane (DDT) and Deltamethrin were carried out on adults aged 2-5 days. Quantitative reverse transcription polymerase chain reaction (RT-qPCR) was performed to quantify the expression of eight metabolic genes involved in mosquito resistance to insecticides.  Results: Susceptibility to DDT showed a similar increase in the time required to knock down 50% of mosquitoes (kdt50) in l colonies exposed to insecticides and non-insecticides compared to the control colony. As for deltamethrin, kdt50 was higher in the colonies exposed to insecticides and the pesticide mixture compared to the colony exposed to non-insecticides. Of all the genes studied in all colonies, except for CYP6P1 induced only in the colony consisting of the pesticide mixture, no genes were induced. Conclusions: This study confirmed that induction is influenced by the duration, the concentration of the solution and the type of xenobiotic used as an inducer. The overexpression of CYP6P1 confirmed the inductive effect that a short exposure of mosquito larvae to agricultural pesticides could have.

Identification and characterization of CYPs induced in the Drosophila antenna by exposure to a plant odorant

Members of the cytochrome p450 (CYP) enzyme family are abundantly expressed in insect olfactory tissues, where they are thought to act as Odorant Degrading Enzymes (ODEs). However, their contribution to olfactory signaling in vivo is poorly understood. This is due in part to the challenge of identifying which of the dozens of antennal-expressed CYPs might inactivate a given odorant. Here, we tested a high-throughput deorphanization strategy in Drosophila to identify CYPs that are transcriptionally induced by exposure to odorants. We discovered three CYPs selectively upregulated by geranyl acetate using transcriptional profiling. Although these CYPs are broadly expressed in the antenna in non-neuronal cells, electrophysiological recordings from CYP mutants did not reveal any changes in olfactory neuron responses to this odorant. Neurons were desensitized by pre-exposing flies to the odorant, but this effect was similar in CYP mutants. Together, our data suggest that the induction of a CYP gene by an odorant does not necessarily indicate a role for that CYP in neuronal responses to that odorant. We go on to show that some CYPs have highly restricted expression patterns in the antenna, and suggest that such CYPs may be useful candidates for further studies on olfactory CYP function.

Overexpression of ATP-binding cassette transporters associated with sulfoxaflor resistance in Aphis gossypii glover

BACKGROUND

Sulfoxaflor is a new insecticide for controlling against Aphis gossypii in the field. ATP-binding cassette (ABC) transporters belong to a large superfamily of proteins and play an important role in the detoxification process. However, the potential role of ABC transporters in sulfoxaflor resistance in A. gossypii is unknown.

RESULTS

In this study, an ABC transporter inhibitor, verapamil, dramatically increased the toxicity of sulfoxaflor in the resistant population with a synergistic ratio of 8.55. However, verapamil did not synergize sulfoxaflor toxicity in the susceptible population. The contents of ABC transporters were significantly increased in the Sul-R population. Based on RT-qPCR analysis, 10 of 23 ABC transcripts, ABCA1, ABCA2, ABCB1, ABCB5, ABCD1, ABCG7, ABCG16, ABCG26, ABCG27, and MRP7, were up-regulated in the Sul-R population compared to the Sus population. Meanwhile, inductive effects of ABCA1, ABCD1, ABCG7 and ABCG26 by sulfoxaflor were found in A. gossypii. Furthermore, knockdown of ABCA1 and ABCD1 using RNAi significantly increased the sulfoxaflor sensitivity in Sul-R aphids.

CONCLUSION

These results suggested that ABC transporters, especially the ABCA1 and ABCD1 genes, might be related with sulfoxaflor resistance in A. gossypii. This study will promote further work to validate the functional roles of these ABCs in sulfoxaflor resistance and might be helpful for the management of sulfoxaflor-resistant A. gossypii.

Pyrethroid metabolism by eleven Helicoverpa armigera P450s from the CYP6B and CYP9A subfamilies

Lepidopteran P450s of the CYP6B and CYP9A subfamilies are thought to play important roles in host plant adaptation and insecticide resistance. An increasing number of paralogs and orthologs with high levels of sequence identity have been found in these subfamilies by mining recent genome projects. However, the biochemical function of most of them remains unknown. A better understanding of the evolution of P450 genes and of the catalytic competence of the enzymes they encode is needed to facilitate studies of host plant use and insecticide resistance. Here, we focused on the full complement of CYP6B (4 genes) and CYP9A (7 genes) in the generalist herbivore, Helicoverpa armigera. These P450s were heterologously expressed in Sf9 cells and compared functionally. In vitro assays showed that all CYP6B and CYP9A P450s can metabolize esfenvalerate efficiently, except for the evolutionarily divergent CYP6B43. A new 2'-hydroxy-metabolite of esfenvalerate was identified and found to be the main metabolite produced by CYP9A12. All tested P450s showed only low induction responses to esfenvalerate. To put these results from H. armigera P450s in perspective, 158 complete CYP6B and 100 complete CYP9A genes from 34 ditrysian species were manually curated. The CYP9A subfamily was more widespread than the CYP6B subfamily and the latter showed dramatic gains and losses, with ten species lacking CYP6B genes. Two adjacent CYP6B loci were found on chromosome 21, with different fates during the evolution of Lepidoptera. The diversity and functional redundancy of CYP6B and CYP9A genes challenge resistance management and pest control strategies as many P450s are available to insects to cope with chemical stresses they encounter.

Honey bee queen health is unaffected by contact exposure to pesticides commonly found in beeswax

Honey bee queen health is crucial for colony health and productivity, and pesticides have been previously associated with queen loss and premature supersedure. Prior research has investigated the effects of indirect pesticide exposure on queens via workers, as well as direct effects on queens during development. However, as adults, queens are in constant contact with wax as they walk on comb and lay eggs; therefore, direct pesticide contact with adult queens is a relevant but seldom investigated exposure route. Here, we conducted laboratory and field experiments to investigate the impacts of topical pesticide exposure on adult queens. We tested six pesticides commonly found in wax: coumaphos, tau-fluvalinate, atrazine, 2,4-DMPF, chlorpyriphos, chlorothalonil, and a cocktail of all six, each administered at 1, 4, 8, 16, and 32 times the concentrations typically found in wax. We found no effect of any treatment on queen mass, sperm viability, or fat body protein expression. In a field trial testing queen topical exposure of a pesticide cocktail, we found no impact on egg-laying pattern, queen mass, emergence mass of daughter workers, and no proteins in the spermathecal fluid were differentially expressed. These experiments consistently show that pesticides commonly found in wax have no direct impact on queen performance, reproduction, or quality metrics at the doses tested. We suggest that previously reported associations between high levels of pesticide residues in wax and queen failure are most likely driven by indirect effects of worker exposure (either through wax or other hive products) on queen care or queen perception.

Short-term metabolic resistance inductive effect of different agrochemical groups on Anopheles gambiae mosquitoes

Background: In order to assess the impact of the different groups of agricultural pesticides used in Côte d'Ivoire on the increase of mosquitoes resistance to insecticides, the expression profiles of 7 P450 cytochromes and one GSTE2 of Anopheles gambiae involved in mosquito resistance to insecticides were studied. The goal of this study was to determine the effect of short exposure of mosquito larvae to different groups of agricultural pesticides on mosquito resistance. Methods: Three groups of pesticides were selected: (i) agricultural insecticide solutions, (ii) none-insecticide pesticide solutions (a mixture of herbicides and fungicides), and (iii) a mixture of the first two. A fourth non-pesticide solution was used as a control. Four groups of each stage 2 larvae (strain Kisumu, male and female) were exposed to 20% concentrated solution for 24 hours. Susceptibility tests for dichlorodiphenyltrichloroethane (DDT) and Deltamethrin were carried out on adults aged 2-5 days. Quantitative reverse transcription polymerase chain reaction (RT-qPCR) was performed to quantify the expression of eight metabolic genes involved in mosquito resistance to insecticides.  Results: Susceptibility to DDT showed a similar increase in the time required to knock down 50% of mosquitoes (kdt50) in l colonies exposed to insecticides and none-insecticides compared to the control colony. As for deltamethrin, kdt50 was higher in the colonies exposed to insecticides and the pesticide mixture compared to the colony exposed to none-insecticides. Of all the genes studied in all colonies, except for CYP6P1 induced only in the colony consisting of the pesticide mixture, no genes were induced. Conclusions: This study confirmed that induction is influenced by the duration, the concentration of the solution and the type of xenobiotic used as an inducer. The overexpression of CYP6P1 confirmed the inductive effect that a short exposure of mosquito larvae to agricultural pesticides could have.

Functional analysis of ABCG and ABCH transporters from the red flour beetle, Tribolium castaneum

BACKGROUND

ATP-binding cassette transporter (ABC transporter) subfamilies ABCA-C and ABCG-H have been implicated in insecticide detoxification, mostly based on findings of elevated gene expression in response to insecticide treatment. We previously characterized TcABCA-C genes from the model beetle and pest Tribolium castaneum and demonstrated that TcABCA and TcABCC genes are involved in the elimination of diflubenzuron, because RNA interference (RNAi)-mediated gene silencing increased susceptibility. In this study, we focused on the potential functions of TcABCG and TcABCH genes in insecticide detoxification.

RESULTS

When we silenced the expression of TcABCG-H genes using RNAi, we noticed a previously unreported developmental RNAi phenotype for TcABCG-4F, which is characterized by 50% mortality and ecdysial arrest during adult moult. When we knocked down the Drosophila brown orthologue TcABCG-XC, we did not obtain apparent eye colour phenotypes but did observe a loss of riboflavin uptake by Malpighian tubules. Next, we determined the expression profiles of all TcABCG-H genes in different tissues and developmental stages and analysed transcript levels in response to treatment with four chemically unrelated insecticides. We found that some genes were specifically upregulated after insecticide treatment. However, when we determined insecticide-induced mortalities in larvae that were treated by double-stranded RNA injection to silence those TcABCG-H genes that were upregulated, we did not observe a significant increase in susceptibility to insecticides.

CONCLUSION

Our findings suggest that the observed insecticide-dependent induction of TcABCG-H gene expression reflects an unspecific stress response, and hence underlines the significance of functional studies on insecticide detoxification. © 2021 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Overexpression of Multiple UDP-Glycosyltransferase Genes Involved in Sulfoxaflor Resistance in Aphis gossypii Glover

UDP-glycosyltransferases (UGTs) are major phase II enzymes involved in the metabolic detoxification of xenobiotics. In this study, two UGT-inhibitors, 5-nitrouracil and sulfinpyrazone, significantly increased sulfoxaflor toxicity against sulfoxaflor-resistant (Sul-R) Aphis gossypii, whereas there were no synergistic effects in susceptible (Sus) A. gossypii. The activity of UGTs in the Sul-R strain was significantly higher (1.35-fold) than that in the Sus strain. Further, gene expression determination demonstrated that 11 of 23 UGT genes were significantly upregulated (1.40- to 5.46-fold) in the Sul-R strain, among which the expression levels of UGT350A2, UGT351A4, UGT350B2, UGT342C2, and UGT343C2 could be induced by sulfoxaflor. Additionally, knockdown of UGT350A2, UGT351A4, UGT350B2, and UGT343C2 using RNA interference (RNAi) significantly increased sensitivity (1.57- to 1.76-fold) to sulfoxaflor in the Sul-R strain. These results suggested that UGTs might be involved in sulfoxaflor resistance in A. gossypii. These findings will facilitate further work to validate the functional roles of these UGT genes in sulfoxaflor resistance.

De Novo Transcriptomic Analyses Revealed Some Detoxification Genes and Related Pathways Responsive to Noposion Yihaogong® 5% EC (Lambda-Cyhalothrin 5%) Exposure in Spodoptera frugiperda Third-Instar Larvae

Simple Summary

Insect pest resistance to synthetic insecticides is a major problem that limits efficient management and thus decreases productivity for farmers and increases the use of harmful materials that pollute the environment and endanger humans and beneficial organisms. A major approach for resistance management is understanding how insect pest field populations develop resistance at molecular levels. To provide a comprehensive insight into the resistance mechanisms of Spodoptera frugiperda larvae to lambda-cyhalothrin 5%, we investigated the molecular basis of resistance mechanism in field collected population of fall armyworm (Spodoptera frugiperda) to lambda-cyhalothrin 5% insecticide, a pyrethroid insecticide by using de novo transcriptomics analysis. We found that resistance to lambda-cyhalothrin 5% can be metabolic by increasing the levels of detoxifying enzymes such as P450, GST and UGT and related genes to insecticide resistance in the field population. The obtained transcriptome information provides large gene resources available for further studying the resistance development of Spodoptera frugiperda to pesticides. The DGE data provide comprehensive insights into the gene expression profiles of fall armyworm (Spodoptera frugiperda) to lambda-cyhalothrin 5% and will facilitate the study of the role of each gene in lambda-cyhalothrin resistance development.

Abstract

The fall armyworm, Spodoptera frugiperda (J.E. Smith) (Lepidoptera: Noctuidae), is a polyphagous, invasive insect pest which causes significant losses in important crops wherever it has spread. The use of pesticides in agriculture is a key tool in the management of many important crop pests, including S. frugiperda, but continued use of insecticides has selected for various types of resistance, including enzyme systems that provide enhanced mechanisms of detoxification. In the present study, we analyzed the de novo transcriptome of S. frugiperda larvae exposed to Noposion Yihaogong^® 5% emulsifiable concentrate (EC) insecticide focusing on detoxification genes and related pathways. Results showed that a total of 1819 differentially expressed genes (DEGs) were identified in larvae after being treated with Noposion Yihaogong^® 5% EC insecticide, of which 863 were up- and 956 down-regulated. Majority of these differentially expressed genes were identified in numerous Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, including metabolism of xenobiotics and drug metabolism. Furthermore, many of S. frugiperda genes involved in detoxification pathways influenced by lambda-cyhalothrin stress support their predicted role by further co-expression network analysis. Our RT-qPCR results were consistent with the DEG’s data of transcriptome analysis. The comprehensive transcriptome sequence resource attained through this study enriches the genomic platform of S. frugiperda, and the identified DEGs may enable greater molecular underpinnings behind the insecticide-resistance mechanism caused by lambda-cyhalothrin.

Genomic and transcriptomic analyses in Drosophila suggest that the ecdysteroid kinase-like (EcKL) gene family encodes the 'detoxification-by-phosphorylation' enzymes of insects

Phosphorylation is a phase II detoxification reaction that, among animals, occurs near exclusively in insects, but the enzymes responsible have never been cloned or otherwise identified. We propose the hypothesis that members of the arthropod-specific ecdysteroid kinase-like (EcKL) gene family encode detoxicative kinases. To test this hypothesis, we annotated the EcKL gene family in 12 species of Drosophila and explored their evolution within the genus. Many ancestral EcKL clades are evolutionarily unstable and have experienced repeated gene gain and loss events, while others are conserved as single-copy orthologs. Leveraging multiple published gene expression datasets from D. melanogaster, and using the cytochrome P450s-a classical detoxification family-as a test case, we demonstrate relationships between xenobiotic induction, detoxification tissue-enriched expression and evolutionary instability in the EcKLs and the P450s. We devised a systematic method for identifying candidate detoxification genes in large gene families that is concordant with experimentally determined functions of P450 genes in D. melanogaster. Applying this method to the EcKLs suggested a significant proportion of these genes play roles in detoxification, and that the EcKLs may constitute a detoxification gene family in insects. Additionally, we estimate that between 11 and 16 uncharacterised D. melanogaster P450s are strong detoxification candidates. Lastly, we also found previously unreported genomic and transcriptomic variation in a number of EcKLs and P450s associated with toxic stress phenotypes using a targeted phenome-wide association study (PheWAS) approach in D. melanogaster, presenting multiple future avenues of research for detoxification genetics in this species.

Transcriptomic identification and characterization of genes responding to sublethal doses of three different insecticides in the western flower thrips, Frankliniella occidentalis

Pretreatment with sublethal concentrations (LC₁₀) of three insecticides (chlorfenapyr, dinotefuran, and spinosad) enhanced tolerance to a lethal dose of the respective insecticide in the Western flower thrips, Frankliniella occidentalis. To identify genes responding to sublethal treatment with insecticides, transcriptome analysis was conducted for thrips treated with LC₁₀ of the three insecticides. When based on a fold change >1.5 or < -1.5 as a selection criterion, 199 transcripts were commonly up-regulated, whereas 31 transcripts were commonly down-regulated following all three insecticide treatments. The differential expression levels of representative genes were validated by quantitative PCR. Most over-transcribed transcripts could be categorized as basic biological processes, such as proteolysis and lipid metabolism. Detoxification genes, such as one glutathione S transferase S1, three UDP-glucuronosyltransferases, four CYP450s, and one ABC transporter G family member 20, were commonly overexpressed in all three insecticide-treated groups. Knockdown of the five representative commonly overexpressed genes via ingestion RNA interference increased mortalities to all the three test insecticides, supporting their common role in tolerance induction. In contrast, three C₂H₂-type zinc finger-containing proteins were significantly down-regulated in all insecticide-treated thrip groups. Since the tested insecticides have distinct structures and modes of action, the roles of commonly expressed genes in tolerance were discussed.

Behavioral and Transcriptional Response to Selection for Olfactory Behavior in Drosophila

The detection, discrimination, and behavioral responses to chemical cues in the environment can have marked effects on organismal survival and reproduction, eliciting attractive or aversive behavior. To gain insight into mechanisms mediating this hedonic valence, we applied thirty generations of divergent artificial selection for Drosophila melanogaster olfactory behavior. We independently selected for positive and negative behavioral responses to two ecologically relevant chemical compounds: 2,3-butanedione and cyclohexanone. We also tested the correlated responses to selection by testing behavioral responses to other odorants and life history traits. Measurements of behavioral responses of the selected lines and unselected controls to additional odorants showed that the mechanisms underlying responses to these odorants are, in some cases, differentially affected by selection regime and generalization of the response to other odorants was only detected in the 2,3-butanedione selection lines. Food consumption and lifespan varied with selection regime and, at times, sex. An analysis of gene expression of both selection regimes identified multiple differentially expressed genes. New genes and genes previously identified in mediating olfactory behavior were identified. In particular, we found functional enrichment of several gene ontology terms, including cell-cell adhesion and sulfur compound metabolic process, the latter including genes belonging to the glutathione S-transferase family. These findings highlight a potential role for glutathione S-transferases in the evolution of hedonic valence to ecologically relevant volatile compounds and set the stage for a detailed investigation into mechanisms by which these genes mediate attraction and aversion.

Resource profitability, but not caffeine, affects individual and collective foraging in the stingless bee Plebeia droryana

Plants and pollinators form beneficial relationships, with plants offering resources in return for pollination services. Some plants, however, add compounds to nectar to manipulate pollinators. Caffeine is a secondary plant metabolite found in some nectars that affects foraging in pollinators. In honeybees, caffeine increases foraging and recruitment to mediocre food sources, which might benefit the plant, but potentially harms the colonies. For the largest group of social bees, the stingless bees, the effect of caffeine on foraging behaviour has not been tested yet, despite their importance for tropical ecosystems. More generally, recruitment and foraging dynamics are not well understood in most species. We examined whether caffeine affects the foraging behaviour of the stingless bee Plebeia droryana, which frequently visits plants that produce caffeinated nectar and pollen. We trained bees to food sources containing field-realistic concentrations of sugar and caffeine. Caffeine did not cause P. droryana to increase foraging frequency and persistence. We observed P. droryana recruiting to food sources; however, this behaviour was also not affected by caffeine. Instead we found that higher sugar concentrations caused bees to increase foraging effort. Thus, unlike in other pollinators, foraging behaviour in this stingless bee is not affected by caffeine. As the Brazilian P. droryana population that we tested has been exposed to coffee over evolutionary time periods, our results raise the possibility that it may have evolved a tolerance towards this central nervous system stimulant. Alternatively, stingless bees may show physiological responses to caffeine that differ from those of other bee groups.

Cis- and trans-acting variants contribute to survivorship in a naïve Drosophila melanogaster population exposed to ryanoid insecticides

Around the world insecticides are being deregistered and banned, as their environmental costs are deemed too great or their efficacy against pest insects is reduced through the evolution of insecticide resistance. With the introduction of replacement insecticides comes the responsibility to assess the way new insecticides perturb various levels of biological systems, from insect physiology to ecosystems. We used a systems genetics approach to identify genetic variants affecting survivorship of Drosophila melanogaster exposed to chlorantraniliprole. The study population was completely naïve to this insecticide chemistry and yet we find associations with variants in neuromuscular genes and coregulated detoxification genes. We predict that these variants will increase in populations of this “sentinel species” as these insecticides are applied in the environment.

Insecticide resistance is a paradigm of microevolution, and insecticides are responsible for the strongest cases of recent selection in the genome of Drosophila melanogaster. Here we use a naïve population and a novel insecticide class to examine the ab initio genetic architecture of a potential selective response. Genome-wide association studies (GWAS) of chlorantraniliprole susceptibility reveal variation in a gene of major effect, Stretchin Myosin light chain kinase (Strn-Mlck), which we validate with linkage mapping and transgenic manipulation of gene expression. We propose that allelic variation in Strn-Mlck alters sensitivity to the calcium depletion attributable to chlorantraniliprole’s mode of action. GWAS also reveal a network of genes involved in neuromuscular biology. In contrast, phenotype to transcriptome associations identify differences in constitutive levels of multiple transcripts regulated by cnc, the homolog of mammalian Nrf2. This suggests that genetic variation acts in trans to regulate multiple metabolic enzymes in this pathway. The most outstanding association is with the transcription level of Cyp12d1 which is also affected in cis by copy number variation. Transgenic overexpression of Cyp12d1 reduces susceptibility to both chlorantraniliprole and the closely related insecticide cyantraniliprole. This systems genetics study reveals multiple allelic variants segregating at intermediate frequency in a population that is completely naïve to this new insecticide chemistry and it foreshadows a selective response among natural populations to these chemicals.

Effects of cadmium on fecundity and defence ability of Drosophila melanogaster

Cadmium (chemical symbol, Cd) is an extremely common pollutant that poses a toxicity threat to organisms. Therefore, we tested Drosophila melanogaster fecundity, Cd accumulation, and activity of two enzymes following Cd stress and used quantitative real-time polymerase chain reaction (qPCR) to quantify the mRNA expression levels of several genes involved in fecundity and defence. D. melanogaster was placed in a medium containing different concentrations of Cd (13, 26, and 52 mg L^-1), following which, inductively coupled plasma atomic emission spectroscopy showed that Cd accumulation in Drosophila increased with the increase in its dietary intake. We also observed that Cd at these concentrations significantly prolonged the mating latency in females and reduced the number of eggs laid. However, the same Cd concentrations did not affect male fecundity. Acetylcholinesterase activity was only detected at 52 mg L^-1 Cd in both sexes, whereas glutathione S-transferase activity was inhibited at 26 and 52 mg L^-1 Cd in females. The results of qPCR indicated that exposure to 13-52 mg L^-1 Cd affected the expression of reproduction-related genes, including downregulation of enok and upregulation of dally and dpp. The same level of exposure also induced transcriptional responses from three defence-related genes (hsp70, gstd2, and gstd6). Taken together, the results revealed that Cd exposure might negatively affect the expression of genes associated with D. melanogaster reproduction and trigger the transcription of defence-related genes. We suggest that further analyses of fecundity and defence responses may help develop indicators of Cd toxicity and improve our understanding of antitoxin defences.

The expression of Spodoptera exigua P450 and UGT genes: tissue specificity and response to insecticides

Cytochrome P450 and UDP-glucosyltransferase (UGT) as phase I and phase II metabolism enzymes, respectively, play vital roles in the breakdown of endobiotics and xenobiotics. Insects can increase the expression of detoxification enzymes to cope with the stress from xenobiotics including insecticides. However, the molecular mechanisms for insecticide detoxification in Spodoptera exigua remain elusive, and the genes conferring insecticide metabolisms in this species are less well reported. In this study, 68 P450 and 32 UGT genes were identified. Phylogenetic analysis showed gene expansions in CYP3 and CYP4 clans of P450 genes and UGT33 family of this pest. P450 and UGT genes exhibited specific tissue expression patterns. Insecticide treatments in fat body cells of S. exigua revealed that the expression levels of P450 and UGT genes were significantly influenced by challenges of abamectin, lambda-cyhalothrin, chlorantraniliprole, metaflumizone and indoxacarb. Multiple genes for detoxification were affected in expression levels after insecticide exposures. The results demonstrated that lambda-cyhalothrin, chlorantraniliprole, metaflumizone and indoxacarb induced similar responses in the expression of P450 and UGT genes in fat body cells; eight P450 genes and four UGT genes were co-up-regulated significantly, and no or only a few CYP/UGT genes were down-regulated significantly by these four insecticides. However, abamectin triggered a distinct response for P450 and UGT gene expression; more P450 and UGT genes were down-regulated by abamectin than by the other four compounds. In conclusion, P450 and UGT genes from S. exigua were identified, and different responses to abamectin suggest a different mechanism for insecticide detoxification.

CYP4CJ1-mediated gossypol and tannic acid tolerance in Aphis gossypii Glover

Cytochrome P450 monooxygenases play a key role in herbivorous pest adaptation to host plants by the detoxification against plant allelochemicals. A new P450 gene (CYP4CJ1) was identified from Aphis gossypii, which displayed a positive response to plant allelochemicals. The transcript levels of CYP4CJ1 could be significantly induced by both gossypol and tannic acid. Knockdown of CYP4CJ1 increased the sensitivity of A. gossypii to these two plant allelochemicals. These results suggest that CYP4CJ1 could be involved in the tolerance of A. gossypii to some plant allelochemicals. Subsequently, we examined the regulatory mechanism of CYP4CJ1 based on the transcriptional and post-transcriptional level. A promoter region from -1422 to -1166 of CYP4CJ1 was identified, which was an essential plant allelochemical responsive region. In addition, miR-4133-3p was found to participate in the regulation of CYP4CJ1 post-transcriptionally. Our results suggest that the transcript abundance of CYP4CJ1, following the exposure of A. gossypii to gossypol and tannic acid can be attributed to both the transcriptional and post-transcriptional regulation mechanisms. These results are important for understanding the roles of P450s in the plant allelochemical tolerance of A. gossypii.

Impacts of Sub-lethal DDT Exposures on microRNA and Putative Target Transcript Expression in DDT Resistant and Susceptible Drosophila melanogaster Strains

Ten constitutively differentially expressed miRNAs were previously described between DDT-resistant 91-R and -susceptible control Drosophila melanogaster strains, and among their predicted target genes were those associated with metabolic DDT resistance mechanisms. The present study evaluated the inducibility of miRNA expression and putative downstream regulation of cytochrome P450s in response to DDT exposure in a time-dependent manner in 91-R and the susceptible Canton-S strain. Specifically, RT-qPCR analysis showed that DDT exposures led to the significant down-regulation (repression) of miR-310-3p, miR-311-3p, miR-312-3p, miR-313-3p, and miR-92a-3p levels in Canton-S. This is contrasted with the lack of significant changes in 91-R at most time-points following DDT exposure. The levels of expression among miRNAs exhibited opposite expression patterns compared to their corresponding putative target cytochrome P450s at the same time points after DDT exposure. Collectively, results from this study suggest that miR-310-3p, miR-311-3p, miR-312-3p, miR-313-3p, and miR-92a-3p might have a potential role in the control of DDT detoxification through the post-transcriptional regulation of target cytochrome P450s in Canton-S. Conversely, the lack of significant changes of these same miRNAs in 91-R following DDT-exposure suggests a possible adaptive mutation that removes repressive control mechanisms. These data are important for the understanding impact of adaptive changes in miRNA expression on post-transcriptional regulatory mechanism involved in the evolution of DDT resistance in 91-R.

Structural Variants and Selective Sweep Foci Contribute to Insecticide Resistance in the Drosophila Genetic Reference Panel

Patterns of nucleotide polymorphism within populations of Drosophila melanogaster suggest that insecticides have been the selective agents driving the strongest recent bouts of positive selection. However, there is a need to explicitly link selective sweeps to the particular insecticide phenotypes that could plausibly account for the drastic selective responses that are observed in these non-target insects. Here, we screen the Drosophila Genetic Reference Panel with two common insecticides; malathion (an organophosphate) and permethrin (a pyrethroid). Genome-wide association studies map survival on malathion to two of the largest sweeps in the D. melanogaster genome; Ace and Cyp6g1. Malathion survivorship also correlates with lines which have high levels of Cyp12d1, Jheh1 and Jheh2 transcript abundance. Permethrin phenotypes map to the largest cluster of P450 genes in the Drosophila genome, however in contrast to a selective sweep driven by insecticide use, the derived allele seems to be associated with susceptibility. These results underscore previous findings that highlight the importance of structural variation to insecticide phenotypes: Cyp6g1 exhibits copy number variation and transposable element insertions, Cyp12d1 is tandemly duplicated, the Jheh loci are associated with a Bari1 transposable element insertion, and a Cyp6a17 deletion is associated with susceptibility.

RNA interference validation of detoxification genes involved in ivermectin tolerance in Drosophila melanogaster

Previously, we observed increased transcription levels of specific cytochrome P450 monooxygenase (P450) and adenosine triphosphate binding cassette (ABC) transporter genes in human body lice, Pediculus humanus humanus, following exposure to ivermectin using the non-invasive induction assay, which resulted in tolerance. To confirm the roles of these genes in induction and tolerance, the robust genetic model insect Drosophila melanogaster was chosen. Orthologous genes corresponding to the body louse P450 (Cyp9f2, Cyp6g2 and Cyp9h1) and ABC transporter (Mrp1, GC1824 as an ABCB type and CG3327 as an ABCG type) genes were selected for in vivo bioassay. Following a brief treatment with a sublethal dose of ivermectin, the mortality response was significantly slower, indicating the presence of tolerance. Concurrently, the transcription levels of Cyp9f2 and Mrp1 at 3 h and those of Cyp6g2, Cyp9h1, Mrp1, CG1824 and CG3327 at 6 h post-treatment were upregulated, indicating gene induction. In behavioural bioassay using GAL4/UAS-RNA interference transgenic fly lines, increased susceptibility to ivermectin was observed following heat shock in the Cyp9f2 , Cyp6g2 , Cyp9h1 , Mrp1 or CG3327-knockdown flies. Considering that these five genes are orthologous to those which had the largest over-expression level following ivermectin-induced tolerance in the body louse, the current results suggest that they are also associated with ivermectin detoxification in D. melanogaster and that body lice and D. melanogaster are likely to share, in part, similar mechanisms of tolerance to ivermectin.

Characterization of a Drosophila glutathione transferase involved in isothiocyanate detoxification

Glutathione transferases (GSTs) are ubiquitous key enzymes that catalyse the conjugation of glutathione to xenobiotic compounds in the detoxification process. GSTs have been proposed to play a dual role in the signal termination of insect chemodetection by modifying odorant and tasting molecules and by protecting the chemosensory system. Among the 40 GSTs identified in Drosophila melanogaster, the Delta and Epsilon groups are insect-specific. GSTs Delta and Epsilon may have evolved to serve in detoxification, and have been associated with insecticide resistance. Here, we report the heterologous expression and purification of the D. melanogaster GST Delta 2 (GSTD2). We investigated the capacity of GSTD2 to bind tasting molecules. Among them, we found that isothiocyanates (ITC), insecticidal compounds naturally present in cruciferous plant and perceived as bitter, are good substrates for GSTD2. The X-ray structure of GSTD2 was solved, showing the absence of the classical Ser catalytic residue, conserved in the Delta and Epsilon GSTs. Using molecular dynamics, the interaction of ITC with the GSTD2 three-dimensional structure is analysed and discussed. These findings allow us to consider a biological role for GSTD2 in chemoperception, considering GSTD2 expression in the chemosensory organs and the potential consequences of insect exposure to ITC.

Association of detoxification enzymes with butene-fipronil in larvae and adults of Drosophila melanogaster

Insecticide resistance is a major challenge in successful insect pest control as the insects have the ability to develop resistance to various widely used insecticides. Butene-fipronil is a novel compound with high toxicity to insects and less toxicity to the non-target organisms. In the present study, the effect of butene-fipronil alone and in combination with three enzyme inhibitors, piperonyl butoxide (PBO), diethyl maleate (DEM), and triphenyl phosphate (TPP), was carried out on larvae and adults of Drosophilia melanogaster. Our results indicated that the co-toxicity indices of butene-fipronil + PBO, butene-fipronil + TPP, and butene-fipronil + DEM mixtures were 437.3, 335.0, and 210.3, respectively, in the second-instar larvae, while 186.6, 256.2, and 238.5, respectively, in the adults, indicating synergistic effects. Interestingly, butene-fipronil increased the expression of CYP28A5 in the larvae; CYP9F2, CYP304A1, CYP28A5, and CYP318A1 in the female adults; and CYP303A1 and CYP28A5 in the male adults. Furthermore, high-level expression of Est-7 was observed in the female adults compared to larvae and male adults. Our results suggest that there is no difference in butene-fipronil metabolism in larvae and male and female adults of D. melanogaster.

Transcriptome-based identification and characterization of genes commonly responding to five different insecticides in the diamondback moth, Plutella xylostella

When the 3rd instar larvae of the diamondback moth (DBM), Plutella xylostella, were pretreated with sublethal doses (LC₁₀) and then subsequently exposed to lethal doses (LC₅₀) of chlorantraniliprole, cypermethrin, dinotefuran, indoxacarb and spinosad via leaf dipping, their tolerance to insecticides was significantly enhanced. To identify genes that commonly respond to the treatment of different insecticides and are responsible for the tolerance enhancement, transcriptomic profiles of larvae treated with sublethal doses of the five insecticides were compared with that of untreated control. A total of 117,181 transcripts with a mean length of 662bp were generated by de novo assembly, of which 35,329 transcripts were annotated. Among them, 125, 143, 182, 215 and 149 transcripts were determined to be up-regulated whereas 67, 45, 60, 60 and 38 genes were down-regulated following treatments with chlorantraniliprole, cypermethrin, dinotefuran, indoxacarb and spinosad, respectively. Gene ontology (GO) analysis of differentially expressed genes (DEGs) revealed little differences in their GO profiles between treatments with different insecticides except for spinosad. Finally, the DEGs commonly responding to all insecticides were selected for further characterization, and some of their over-transcription levels were confirmed by quantitative PCR. The most notable examples of commonly responding over-transcribed genes were two cytochrome P450 genes (Cyp301a1 and Cyp9e2) and nine cuticular protein genes. In contrast, several genes composing the mitochondrial energy generation system were significantly down-regulated in all treated larvae. Considering the distinct structure and mode of action of the five insecticides tested, the differentially expressed genes identified in this study appear to be involved in general chemical defense at the initial stage of intoxication. Their possible roles in the tolerance/resistance development were discussed.

Carboxylesterase genes in pyrethroid resistant house flies, Musca domestica

Carboxylesterases are one of the major enzyme families involved in the detoxification of pyrethroids. Up-regulation of carboxylesterase genes is thought to be a major component of insecticide resistant mechanisms in insects. Based on the house fly transcriptome and genome database, a total of 39 carboxylesterase genes of different functional clades have been identified in house flies. In this study, eleven of these genes were found to be significantly overexpressed in the resistant ALHF house fly strain compared with susceptible aabys and wild-type CS strains. Eight up-regulated carboxylesterase genes with their expression levels were further induced to a higher level in response to permethrin treatments, indicating that constitutive and inductive overexpression of carboxylesterases are co-responsible for the enhanced detoxification of insecticides. Spatial expression studies revealed these up-regulated genes to be abundantly distributed in fat bodies and genetically mapped on autosome 2 or 3 of house flies, and their expression could be regulated by factors on autosome 1, 2 and 5. Taken together, these results demonstrate that multiple carboxylesterase genes are co-upregulated in resistant house flies, providing further evidence for their involvement in the detoxification of insecticides and development of insecticide resistance.

Detoxifying symbiosis: microbe-mediated detoxification of phytotoxins and pesticides in insects

Symbiotic microorganisms degrade natural and artificial toxic compounds, and confer toxin resistance on insect hosts.

De novo characterization of venom apparatus transcriptome of Pardosa pseudoannulata and analysis of its gene expression in response to Bt protein

Background

Pardosa pseudoannulata is a prevailing spider species, and has been regarded as an important bio-control agent of insect pests in farmland of China. However, the available genomic and transcriptomic databases of P. pseudoannulata and their venom are limited, which severely hampers functional genomic analysis of P. pseudoannulata. Recently high-throughput sequencing technology has been proved to be an efficient tool for profiling the transcriptome of relevant non-target organisms exposed to Bacillus thuringiensis (Bt) protein through food webs.

Results

In this study, the transcriptome of the venom apparatus was analyzed. A total of 113,358 non-redundant unigenes were yielded, among which 34,041 unigenes with complete or various length encoding regions were assigned biological function annotations and annotated with gene ontology and karyotic orthologous group terms. In addition, 3726 unigenes involved in response to stimulus and 720 unigenes associated with immune-response pathways were identified. Furthermore, we investigated transcriptomic changes in the venom apparatus using tag-based DGE technique. A total of 1724 differentially expressed genes (DEGs) were detected, while 75 and 372 DEGs were functionally annotated with KEGG pathways and GO terms, respectively. qPCR analyses were performed to verify the DEGs directly or indirectly related to immune and stress responses, including genes encoding heat shock protein, toll-like receptor, GST and NADH dehydrogenase.

Conclusion

This is the first study conducted to specifically investigate the venom apparatus of P. pseudoannulata in response to Bt protein exposure through tritrophic chain. A substantial fraction of transcript sequences was generated by high-throughput sequencing of the venom apparatus of P. pseudoannulata. Then a comparative transcriptome analysis showing a large number of candidate genes involved in immune response were identified by the tag-based DGE technology. This transcriptome dataset will provide a comprehensive sequence resource for furture molecular genetic research of the venom apparatus of P. pseudoannulata.

Electronic supplementary material

The online version of this article (10.1186/s12896-017-0392-z) contains supplementary material, which is available to authorized users.

The effect of insecticide synergist treatment on genome-wide gene expression in a polyphagous pest

Synergists can counteract metabolic insecticide resistance by inhibiting detoxification enzymes or transporters. They are used in commercial formulations of insecticides, but are also frequently used in the elucidation of resistance mechanisms. However, the effect of synergists on genome-wide transcription in arthropods is poorly understood. In this study we used Illumina RNA-sequencing to investigate genome-wide transcriptional responses in an acaricide resistant strain of the spider mite Tetranychus urticae upon exposure to synergists such as S,S,S-tributyl phosphorotrithioate (DEF), diethyl maleate (DEM), piperonyl butoxide (PBO) and cyclosporin A (CsA). Exposure to PBO and DEF resulted in a broad transcriptional response and about one third of the differentially expressed genes (DEGs), including cytochrome P450 monooxygenases and UDP-glycosyltransferases, was shared between both treatments, suggesting common transcriptional regulation. Moreover, both DEF and PBO induced genes that are strongly implicated in acaricide resistance in the respective strain. In contrast, CsA treatment mainly resulted in downregulation of Major Facilitator Superfamily (MFS) genes, while DEGs of the DEM treatment were not significantly enriched for any GO-terms.

Insights into DDT Resistance from the Drosophila melanogaster Genetic Reference Panel

Insecticide resistance is considered a classic model of microevolution, where a strong selective agent is applied to a large natural population, resulting in a change in frequency of alleles that confer resistance. While many insecticide resistance variants have been characterized at the gene level, they are typically single genes of large effect identified in highly resistant pest species. In contrast, multiple variants have been implicated in DDT resistance in Drosophila melanogaster; however, only the Cyp6g1 locus has previously been shown to be relevant to field populations. Here we use genome-wide association studies (GWAS) to identify DDT-associated polygenes and use selective sweep analyses to assess their adaptive significance. We identify and verify two candidate DDT resistance loci. A largely uncharacterized gene, CG10737, has a function in muscles that ameliorates the effects of DDT, while a putative detoxifying P450, Cyp6w1, shows compelling evidence of positive selection.

Interaction of Ferulic Acid with Glutathione S-Transferase and Carboxylesterase Genes in the Brown Planthopper, Nilaparvata lugens

Plant phenolics are crucial defense phytochemicals against herbivores and glutathione S-transferase (GST) and carboxylesterase (CarE) in herbivorous insects are well-known detoxification enzymes for such xenobiotics. To understand relationship between a plant phenolic and herbivore GST or CarE genes, we evaluated the relationship between a rice phenolic ferulic acid and resistance to brown planthopper (BPH, Nilaparvata lugens), and investigated the interaction of ferulic acid with GST or CarE genes in BPH. The results indicate that ferulic acid content in tested rice varieties was highly associated with resistance to BPH. Bioassays using artificial diets show that the phenolic acid toxicity to BPH was dose dependent and the LC₂₅ and LC₅₀ were 5.81 and 23.30 μg/ml at 72 hr, respectively. Activities of the enzymes BPH GST and CarE were increased at concentrations below the LC₅₀ of ferulic acid. Moreover, low ferulic acid concentrations (< LC₂₅) upregulated the transcriptional levels of NlGSTD1 and NlGSTE1 of the GST family and NlCE of the CarE family. By using dsRNA-induced gene silencing (DIGS) of GST or CarE, it was shown that suppressed expression levels of NlGSTD1, NlGSTE1 and NlCE were 14.6%-21.2%, 27.8%-34.2%, and 10.5%-19.8%, respectively. Combination of NlGSTD1, NlGSTE1 or NlCE knockdown with ferulic acid increased nymph mortality by 92.9%, 119.9%, or 124.6%, respectively. These results suggest that depletion of detoxification genes in herbivorous insects by plant-mediated RNAi technology might be a new potential resource for improving rice resistance to BPH.

Expression profiles of glutathione S-transferase superfamily in Spodoptera litura tolerated to sublethal doses of chlorpyrifos

Chlorpyrifos (CPF) is a broad-spectrum organophosphate insecticide. Glutathione S-transferases (GSTs) in insects are a family of detoxification enzymes and they play critical roles in CPF detoxification. Spodoptera litura is one of the most destructive agricultural pests in tropical and subtropical areas in the world. In this study, 37 Slgsts from 46 unique transcripts of gsts in S. litura transcriptome data, including eight previously reported GSTs, were identified and their expression patterns in susceptible and 12-generation-CPF-treated strains were analyzed to understand the roles of these Slgsts in sublethal doses of CPF tolerance. The results indicate that the members of the S. litura GST superfamily could be distinguished into three major groups: one group, including six cytosolic Slgsts (SlGSTe1, SlGSTe3, SlGSTe10, SlGSTe15, SlGSTo2 and SlGSTs5) and two microsomal Slgsts (SlMGST1-2 and SlMGST1-3), was directly responsible for CPF induction in both 12-generation-treated and susceptible strains; the second group, including three cytosolic Slgsts (SlGSTe13, SlGSTt1 and SlGSTz1) and one microsomal Slgst (SlMGST1-1), was induced only in the 12-generation-treated strain; the third group, including eight cytosolic Slgsts (two epsilon, three delta, one omega, one zeta and one unclassified Slgst), was expressed 1.52-5.15-fold higher in the 12-generation-treated strain than in the susceptible strain.

Transcriptome modulation during host shift is driven by secondary metabolites in desert Drosophila

High-throughput transcriptome studies are breaking new ground to investigate the responses that organisms deploy in alternative environments. Nevertheless, much remains to be understood about the genetic basis of host plant adaptation. Here, we investigate genome-wide expression in the fly Drosophila buzzatii raised in different conditions. This species uses decaying tissues of cactus of the genus Opuntia as primary rearing substrate and secondarily, the necrotic tissues of the columnar cactus Trichocereus terscheckii. The latter constitutes a harmful host, rich in mescaline and other related phenylethylamine alkaloids. We assessed the transcriptomic responses of larvae reared in Opuntia sulphurea and T. terscheckii, with and without the addition of alkaloids extracted from the latter. Whole-genome expression profiles were massively modulated by the rearing environment, mainly by the presence of T. terscheckii alkaloids. Differentially expressed genes were mainly related to detoxification, oxidation-reduction and stress response; however, we also found genes involved in development and neurobiological processes. In conclusion, our study contributes new data onto the role of transcriptional plasticity in response to alternative rearing environments.

CHARACTERIZATION AND EXPRESSION PROFILES OF FIVE POSSIBLE CYTOCHROME P450 GENES FROM Liposcelis entomophila (ENDERLEIN) (PSOCOPTERA: LIPOSCELIDIDAE)

In this study, the cDNAs of five cytochromes P450 genes (named CYP345P1, CYP358B1, CYP4FD2, CYP4CD2, and CYP6JN1) contained open reading frames from 1,500 to 1,554 nucleotides that encoded 499 to 517 amino acids were cloned from the psocid Liposcelis entomophila. They are characterized by predicted molecular weights from 57.67 to 59.64 kDa and theoretical isoelectric points of 5.57-9.07. Quantitative real-time PCR analysis showed these five genes were expressed at all tested developmental stages and higher expressions were observed in adults. CYP358B1 was expressed at higher levels in egg and adult compared to the larval stages. mRNA abundances of five genes were detected in both sexes and were relatively more abundant in adult females than in adult males. Synergism bioassay showed that the synergic ratio was 2.20 and 2.45 when insects were treated with the mixture of deltamethrin or malathion with the synergist piperonyl butoxide (PBO). Because PBO induces cytochrome P450s in some insects, this suggested to us that cytochromes P450 might participate in detoxification of these insecticides. The transcripts of the five cytochromes P450 genes in adult psocids could be induced to the highest level at 12 h after the exposure to malathion. After exposure to deltamethrin, CYP358B1 reached maximum expression at 24 h. The maximum expression of the other four genes occurred at 36 h. Treatments with the carbamate propoxur did not influence transcription of the cytochromes P450 gene. The induction profiles suggested that these five cytochrome P450 genes may be associated with deltamethrin and malathion metabolism in psocids.

Glutathione S-Transferase Regulation in Calanus finmarchicus Feeding on the Toxic Dinoflagellate Alexandrium fundyense

The effect of the dinoflagellate, Alexandrium fundyense, on relative expression of glutathione S-transferase (GST) transcripts was examined in the copepod Calanus finmarchicus. Adult females were fed for 5-days on one of three experimental diets: control (100% Rhodomonas spp.), low dose of A. fundyense (25% by volume, 75% Rhodomonas spp.), and high dose (100% A. fundyense). Relative expression of three GST genes was measured using RT-qPCR on days 0.5, 1, 2 and 5 in two independent experiments. Differential regulation was found for the Delta and the Sigma GSTs between 0.5 to 2 days, but not on day 5 in both experiments. The third GST, a microsomal, was not differentially expressed in either treatment or day. RT-qPCR results from the two experiments were similar, even though experimental females were collected from the Gulf of Maine on different dates and their reproductive output differed. In the second experiment, expression of 39 GSTs was determined on days 2 and 5 using RNA-Seq. Global gene expression analyses agreed with the RT-qPCR results. Furthermore, the RNA-Seq measurements indicated that only four GSTs were differentially expressed under the experimental conditions, and the response was small in amplitude. In summary, the A. fundyense diet led to a rapid and transient response in C. finmarchicus in three cytosolic GSTs, while a fourth GST (Omega I) was significantly up-regulated on day 5. Although there was some regulation of GSTs in response the toxic dinoflagellate, the tolerance to A. fundyense by C. finmarchicus is not dependent on the long-term up-regulation of specific GSTs.

Transcriptomic responses of the calanoid copepod Calanus finmarchicus to the saxitoxin producing dinoflagellate Alexandrium fundyense

In the Gulf of Maine, the copepod Calanus finmarchicus co-occurs with the neurotoxin-producing dinoflagellate, Alexandrium fundyense. The copepod is resistant to this toxic alga, but little is known about other effects. Gene expression profiles were used to investigate the physiological response of females feeding for two and five days on a control diet or a diet containing either a low or a high dose of A. fundyense. The physiological responses to the two experimental diets were similar, but changed between the time points. At 5-days the response was characterized by down-regulated genes involved in energy metabolism. Detoxification was not a major component of the response. Instead, genes involved in digestion were consistently regulated, suggesting that food assimilation may have been affected. Thus, predicted increases in the frequency of blooms of A. fundyense could affect C. finmarchicus populations by changing the individuals' energy budget and reducing their ability to build lipid reserves.

Cloning and Expression of Multiple Cytochrome P450 Genes: Induction by Fipronil in Workers of the Red Imported Fire Ant (Solenopsis invicta Buren)

Both exogenous and endogenous compounds can induce the expression of cytochrome P450 genes. The insect cytochrome P450 genes related to insecticide resistance are likely to be expressed as the “first line of defense” when challenged with insecticides. In this study, four cytochrome P450 genes, SinvCYP6B1, SinvCYP6A1, SinvCYP4C1, and SinvCYP4G15, were firstly isolated from workers of the red imported fire ant (Solenopsis invicta) through rapid amplification of cDNA ends (RACE) and sequenced. The fipronil induction profiles of the four cytochrome P450 genes and the two previously isolated CYP4AB1 and CYP4AB2 were characterized in workers. The results revealed that the expression of SinvCYP6B1, SinvCYP6A1, CYP4AB2, and SinvCYP4G15, increased 1.4-fold and 1.3-fold more than those of acetone control, respectively, after 24 h exposure to fipronil at concentrations of 0.25 μg mL⁻¹ (median lethal dose) and 0.56 μg mL⁻¹ (90% lethal dose), while no significant induction of the expression of CYP4AB1 and SinvCYP4C1 was detected. Among these genes, SinvCYP6B1 was the most significantly induced, and its maximum expression was 3.6-fold higher than that in acetone control. These results might suggest that multiple cytochrome P450 genes are co-up-regulated in workers of the fire ant through induction mechanism when challenged with fipronil. These findings indicated that cytochrome P450 genes play an important role in the detoxification of insecticides and provide a theoretical basis for the mechanisms of insecticide metabolism in the fire ant.

Identifying Loci Contributing to Natural Variation in Xenobiotic Resistance in Drosophila

Natural populations exhibit a great deal of interindividual genetic variation in the response to toxins, exemplified by the variable clinical efficacy of pharmaceutical drugs in humans, and the evolution of pesticide resistant insects. Such variation can result from several phenomena, including variable metabolic detoxification of the xenobiotic, and differential sensitivity of the molecular target of the toxin. Our goal is to genetically dissect variation in the response to xenobiotics, and characterize naturally-segregating polymorphisms that modulate toxicity. Here, we use the Drosophila Synthetic Population Resource (DSPR), a multiparent advanced intercross panel of recombinant inbred lines, to identify QTL (Quantitative Trait Loci) underlying xenobiotic resistance, and employ caffeine as a model toxic compound. Phenotyping over 1,700 genotypes led to the identification of ten QTL, each explaining 4.5-14.4% of the broad-sense heritability for caffeine resistance. Four QTL harbor members of the cytochrome P450 family of detoxification enzymes, which represent strong a priori candidate genes. The case is especially strong for Cyp12d1, with multiple lines of evidence indicating the gene causally impacts caffeine resistance. Cyp12d1 is implicated by QTL mapped in both panels of DSPR RILs, is significantly upregulated in the presence of caffeine, and RNAi knockdown robustly decreases caffeine tolerance. Furthermore, copy number variation at Cyp12d1 is strongly associated with phenotype in the DSPR, with a trend in the same direction observed in the DGRP (Drosophila Genetic Reference Panel). No additional plausible causative polymorphisms were observed in a full genomewide association study in the DGRP, or in analyses restricted to QTL regions mapped in the DSPR. Just as in human populations, replicating modest-effect, naturally-segregating causative variants in an association study framework in flies will likely require very large sample sizes.