Identification, genomic organization and expression pattern of glutathione S-transferase in the silkworm, Bombyx mori

Genome sequencing and assembly of Indian golden silkmoth, Antheraea assamensis Helfer (Saturniidae, Lepidoptera)

Muga silkworm (Antheraea assamensis), one of the economically important wild silkmoths, is unique among saturniid silkmoths. It is confined to the North-eastern part of India. Muga silk has the highest value among the other silks. Unlike other silkmoths, A. assamensis has a low chromosome number (n = 15), and ZZ/ZO sex chromosome system. Here, we report the first high-quality draft genome of A. assamensis, assembled by employing the Illumina and PacBio sequencing platforms. The assembled genome of A. assamensis is 501.18 Mb long, with 2697 scaffolds and an N50 of 683.23 Kb. The genome encompasses 18,385 protein-coding genes, 86.29% of which were functionally annotated. Phylogenetic analysis of A. assamensis revealed its divergence from other Antheraea species approximately 28.7 million years ago. Moreover, an investigation into detoxification-related gene families, CYP450, GST, and ABC-transporter, revealed a significant expansion in A. assamensis as compared to the Bombyx mori. This expansion is comparable to Spodoptera litura, suggesting adaptive responses linked to the polyphagous behavior observed in these insects. This study provides valuable insights into the molecular basis of evolutionary divergence and adaptations in muga silkmoth. The genome assembly reported in this study will significantly help in the functional genomics studies on A. assamensis and other Antheraea species along with comparative genomics analyses of Bombycoidea insects.

Comparative transcriptomics reveals the conservation and divergence of reproductive genes across three sympatric Tomicus bark beetles

Three tree-killing bark beetles belonging to the genus Tomicus, Tomicus yunnanensis, Tomicus brevipilosus and Tomicus minor (Coleoptera; Curculionidae, Scolytinae), are serious wood-borers with larvae feeding on the phloem tissues of Pinus yunnanensis. The three Tomicus beetles, in some cases, coexist in a same habitat, providing a best system for exploring the conservation and divergence of reproductive genes. Here, we applied comparative transcriptomics and molecular biology approaches to characterize reproductive-related genes in three sympatric Tomicus species. Illumina sequencing of female and male reproductive systems and residual bodies generated a large number of clean reads, representing 185,920,232 sequences in T. yunnanensis, 169,153,404 in T. brevipilosus and 178,493,176 in T. minor that were assembled into 32,802, 56,912 and 33,670 unigenes, respectively. The majority of the genes had detectable expression in reproductive tissues (FPKM >1), particularly those genes in T. brevipilosus accounting for 76.61 % of the total genes. From the transcriptomes, totally 838 genes encoding 463 detoxification enzymes, 339 chemosensory membrane proteins and 36 ionotropic glutamate receptors (iGluRs) were identified, including 622 reproductive tissue-expressed genes. Of these, members of carboxylesterases (COEs), ionotropic receptors (IRs), sensory neuron membrane proteins (SNMPs) and iGluRs were highly conserved in gene numbers and sequence identities across three Tomicus species. Further, expression profiling analyses revealed a number of genes expressed in reproductive tissues and the diverse expression characteristics in these beetles. The results provide evidence for the conservation and differences of reproductive genes among three sympatric closely related beetles, helping understand their different reproductive strategies and the maximization of the reproductive success.

Characterization of three glutathione S-transferases potentially associated with adaptation of the wheat blossom midge Sitodiplosis mosellana to host plant defense

BACKGROUND

Insect glutathione S-transferases (GSTs), a multifunctional protein family, play a crucial role in detoxification of plant defensive compounds. However, they have been rarely investigated in Sitodiplosis mosellana, a destructive pest of wheat worldwide. In this study, we characterized for the first time a delta (SmGSTd1) and two epsilon GST genes (SmGSTe1 and SmGSTe2) and analyzed their expression patterns and functions associated with adaptation to host plant defense in this species.

RESULTS

Expression of these SmGST genes greatly increased in S. mosellana larvae feeding on resistant wheat varieties Kenong1006, Shanmai139 and Jinmai47 which contain higher tannin and ferulic acid, the major defensive compounds of wheat against this pest, compared with those feeding on susceptible varieties Xinong822, Xinong88 and Xiaoyan22. Their expression was also tissue-specific, most predominant in larval midgut. Recombinant SmGSTs expressed in Escherichia coli could catalyze the conjugation of 1-chloro-2,4-dinitrobenzene, with activity peak at pH around 7.0 and temperature between 30 and 40 °C. Notably, they could metabolize tannin and ferulic acid, with the strongest metabolic ability by SmGSTe2 against two compounds, followed by SmGSTd1 on tannin, and SmGSTe1 on ferulic acid.

CONCLUSION

The results suggest that these SmGSTs are important in metabolizing wheat defensive chemicals during feeding, which may be related to host plant adaptation of S. mosellana. Our study has provided information for future investigation and development of strategies such as host-induced gene silencing of insect-detoxifying genes for managing pest adaptation. © 2023 Society of Chemical Industry.

Bmo-miR-3351 modulates glutathione content and inhibits BmNPV proliferation by targeting BmGSTe6 in Bombyx mori

MicroRNAs (miRNAs) are small non-coding RNAs that play pivotal roles in the host response to invading pathogens. Among these pathogens, Bombyx mori nucleopolyhedrovirus (BmNPV) is one of the main causes of substantial economic losses in sericulture, and there are relatively few studies on the specific functions of miRNAs in the B. mori-BmNPV interaction. Therefore, we conducted transcriptome sequencing to identify differentially expressed (DE) messenger RNAs (mRNAs) and miRNAs in the midgut of 2 B. mori strains (BmNPV-susceptible strain P50 and BmNPV-resistant strain A35) after BmNPV infection. Through correlation analysis of the miRNA and mRNA data, we identified a comprehensive set of 21 miRNAs and 37 predicted target mRNAs. Notably, miR-3351, which has high expression in A35, exhibited remarkable efficacy in suppressing BmNPV proliferation. Additionally, we confirmed that miR-3351 binds to the 3' untranslated region (3' UTR) of B. mori glutathione S-transferase epsilon 6 (BmGSTe6), resulting in its downregulation. Conversely, BmGSTe6 displayed an opposite expression pattern to miR-3351, effectively promoting BmNPV proliferation. Notably, BmGSTe6 levels were positively correlated with glutathione S-transferase activity, consequently influencing intracellular glutathione content in the infected samples. Furthermore, our investigation revealed the protective role of glutathione against BmNPV infection in BmN cells. In summary, miR-3351 modulates glutathione content by downregulating BmGSTe6 to inhibit BmNPV proliferation in B. mori. Our findings enriched the research on the role of B. mori miRNAs in the defense against BmNPV infection, and suggests that the antiviral molecule, glutathione, offers a novel perspective on preventing viral infection in sericulture.

Role of the epsilon glutathione S-transferases in xanthotoxin tolerance in Spodoptera litura

Spodoptera litura, a polyphagous lepidopteran pest, demonstrates a remarkable capacity to adapt to varying host plants by efficiently detoxifying phytochemicals. However, the underlying mechanism for this adaptation is not well understood. Herein, twenty eplison glutathione S-transferase genes (GSTes) were characterized and their roles in phytochemical tolerance were analyzed in S. litura. Most of the GSTe genes were mainly expressed in the larval midgut and fat body. Exposure to the phytochemicals, especially xanthotoxin, induced the expression of most GSTe genes. Molecular docking analysis revealed that xanthotoxin could form stable bonds with six xanthotoxin-responsive GSTes, with binding free energies ranging from -36.44 to -68.83 kcal mol-1. Knockdown of these six GSTe genes increased the larval susceptibility to xanthotoxin. Furthermore, xanthotoxin exposure significantly upregulated the expression of two transcription factor genes CncC and MafK. Silencing of either CncC or MafK reduced the expression of GSTe16, which exhibited the largest change in response to xanthotoxin. Additionally, analysis of the promoter sequence of GSTe16 revealed the presence of seven CncC/Maf binding sites. Luciferase reporter assays showed that CncC and MafK enhanced the expression of GSTe16, leading to the increased xanthotoxin tolerance in S. litura. These findings provide insight into the functions and transcriptional regulatory mechanisms of GSTes, thereby enhancing our understanding of the role of GSTs in the adaptation of lepidopteran pests to phytochemicals.

Exploring the proteins and metabolites associated with male antennae responses to female exposure of Antheraea pernyi (Lepidoptera: Saturniidae) moths

Detection of sex pheromones of insects relies on the antennae. The female pheromone signal transmission in the male antennae ultimately initiates the courtship and mating behaviors of males. To investigate the proteins and metabolites involved in this neural transduction, integrative proteomics and metabolomics analysis including tandem mass tag (TMT) proteomic quantification and liquid chromatography tandem mass spectrometry (LC/MS)-based metabolomics was adopted for comparing proteomic and metabolic changes between the antennae of male moths following stimulation by females and the non-stimulated males of Antheraea pernyi (Guérin-Méneville, Lepidoptera: Saturniidae) in this study. A total of 92 differentially expressed proteins (DEPs) containing 52 upregulated and 40 downregulated proteins and 545 differentially expressed metabolites (DEMs) including 218 upregulated and 327 downregulated metabolites were identified from the antennae of female-stimulated male moths based on the proteome and metabolome data, respectively. Bioinformatics analysis was performed for the 45 DEPs and 160 DEMs, including Gene Ontology (GO), Clusters of Orthologous Groups (COG), and Kyoto Encylopaedia of Genes and Genomes (KEGG) enrichment analysis and Human Metabolome Database (HMDB) annotation. A number of DEPs and DEMs related to neural transmission of female pheromone signals in the male antennae of A. pernyi were screened, including tyrosine hydroxylase, cryptochrome-1, tachykinin, arylalkylamine N-acetyltransferase, cadherin-23, glutathione S-transferase delta 3, tyramine, tryptamine, n-oleoyl dopamine, n-stearoyl dopamine, and n-stearoyl tyrosine. The altered expression levels of those proteins or metabolites were speculated involved in regulating the neuron activity for enhanced transmission of neural impulses and continuous perception, reception, and transduction of female pheromone signals. Our findings yielded novel insights into the potential mechanisms in the antennae of male A. pernyi responding to female attraction.

Identification and functional study of detoxification-related genes in response to tolfenpyrad stress in Glyphodes pyloalis Walker (Lepidoptera: Pyralidae)

Glyphodes pyloalis Walker (G. pyloalis) is a common destructive mulberry pest. Due to the long-term and frequent use of insecticides, it has developed tolerance to commonly used insecticides. Tolfenpyrad (TFP) is a novel pyrazole heterocyclic insecticide. In order to understand the TFP detoxification mechanism of G. pyloalis larvae, we first estimated the LC30 dose of TFP for 3rd instar G. pyloalis larvae. Next, we identified genes that were differentially expressed in 3rd instar G. pyloalis larvae treated with TFP compared to the control group by transcriptome sequencing. In total, 86,949,569 and 67,442,028 clean reads were obtained from TFP-treated and control G. pyloalis larvae, respectively. A total of 5588 differentially expressed genes (DEGs) were identified in TFP-treated and control G. pyloalis larvae, of which 3084 genes were upregulated and 2504 genes were downregulated. We analyzed the expression of 43 candidate detoxification enzyme genes associated with insecticide tolerance using qPCR. According to the spatiotemporal expression pattern of DEGs, we found that CYP6ABE1, CYP333A36 and GST-epsilon8 were highly expressed in the midgut, while CarEs14 was strongly expressed in haemolymph. Furthermore, we successfully knocked down these genes by RNA interference. After silencing CYP6ABE1 and CYP333A36, bioassay showed that the mortality rate of TFP-treated G. pyloalis larvae was significantly higher compared to the control group. This study provides a theoretical foundation for understanding the sensitivity of G. pyloalis to TFP and establish the basis for the effective and green management of this pest.

Antennae-enriched expression of candidate odorant degrading enzyme genes in the turnip aphid, Lipaphis erysimi

Aphids heavily rely on their olfactory system for foraging behavior. Odorant-degrading enzymes (ODEs) are essential in preserving the olfactory acuity of aphids by removing redundant odorants in the antennae. Certain enzymes within this group stand out as being enriched and/or biased expressed in the antennae, such as carboxylesterases (CXEs), cytochrome P450 (CYPs), glutathione S-transferases (GSTs), and UDP-glycosyltransferases (UGTs). Here, we performed a comparative transcriptome analysis of antennae and body tissue to isolate the antennal ODE genes of turnip aphid Lipaphis erysimi. A dataset of one CXE, seven CYPs, two GSTs, and five UGTs enriched in the antennae was identified and subjected to sequence analysis. Furthermore, qRT-PCR analyses showed that 13 ODE genes (LeCXE6, LeCYP4c1, LeCYP6a2, LeCYP6a13, LeCYP6a14.2, LeCYP6k1, LeCYP18a1, LeGST1, LeUGT1-7, LeUGT2B7, LeUGT2B13, LeUGT2C1.1, and LeUGT2C1.2) were specifically or significantly elevated in antennal tissues. Among these antennae-enriched ODEs, LeCYP4c1, LeCYP6a2, LeCYP6a13, LeCYP6a14.2, LeCYP18a1, LeUGT2B7, and LeUGT2B13 were found to exhibit significantly higher expression levels in alate aphids compared to apterous and nymph aphids, suggesting their putative role in detecting new host plant location. The results presented in this study highlight the identification and expression of ODE genes in L. erysimi, paving the path to investigate their functional role in odorant degradation during the olfactory processes.

Identification and Characterization of Glutathione S-transferase Genes in Spodoptera frugiperda (Lepidoptera: Noctuidae) under Insecticides Stress

Insect glutathione S-transferases (GSTs) serve critical roles in insecticides and other forms of xenobiotic chemical detoxification. The fall armyworm, Spodoptera frugiperda (J. E. Smith), is a major agricultural pest in several countries, especially Egypt. This is the first study to identify and characterize GST genes in S. frugiperda under insecticidal stress. The present work evaluated the toxicity of emamectin benzoate (EBZ) and chlorantraniliprole (CHP) against the third-instar larvae of S. frugiperda using the leaf disk method. The LC₅₀ values of EBZ and CHP were 0.029 and 1.250 mg/L after 24 h of exposure. Moreover, we identified 31 GST genes, including 28 cytosolic and 3 microsomal SfGSTs from a transcriptome analysis and the genome data of S. frugiperda. Depending on the phylogenetic analysis, sfGSTs were divided into six classes (delta, epsilon, omega, sigma, theta, and microsomal). Furthermore, we investigated the mRNA levels of 28 GST genes using qRT-PCR under EBZ and CHP stress in the third-instar larvae of S. frugiperda. Interestingly, SfGSTe10 and SfGSTe13 stood out with the highest expression after the EBZ and CHP treatments. Finally, a molecular docking model was constructed between EBZ and CHP using the most upregulated genes (SfGSTe10 and SfGSTe13) and the least upregulated genes (SfGSTs1 and SfGSTe2) of S. frugiperda larvae. The molecular docking study showed EBZ and CHP have a high binding affinity with SfGSTe10, with docking energy values of -24.41 and -26.72 kcal/mol, respectively, and sfGSTe13, with docking energy values of -26.85 and -26.78 kcal/mol, respectively. Our findings are important for understanding the role of GSTs in S. frugiperda regarding detoxification processes for EBZ and CHP.

Adaptive Changes in Detoxification Metabolism and Transmembrane Transport of Bombyx mori Malpighian Tubules to Artificial Diet

The high adaptability of insects to food sources has contributed to their ranking among the most abundant and diverse species on Earth. However, the molecular mechanisms underlying the rapid adaptation of insects to different foods remain unclear. We explored the changes in gene expression and metabolic composition of the Malpighian tubules as an important metabolic excretion and detoxification organ in silkworms (Bombyx mori) fed mulberry leaf and artificial diets. A total of 2436 differentially expressed genes (DEGs) and 245 differential metabolites were identified between groups, with the majority of DEGs associated with metabolic detoxification, transmembrane transport, and mitochondrial function. Detoxification enzymes, such as cytochrome P450 (CYP), glutathione-S-transferase (GST), and UDP-glycosyltransferase, and ABC and SLC transporters of endogenous and exogenous solutes were more abundant in the artificial diet group. Enzyme activity assays confirmed increased CYP and GST activity in the Malpighian tubules of the artificial diet-fed group. Metabolome analysis showed increased contents of secondary metabolites, terpenoids, flavonoids, alkaloids, organic acids, lipids, and food additives in the artificial diet group. Our findings highlight the important role of the Malpighian tubules in adaptation to different foods and provide guidance for further optimization of artificial diets to improve silkworm breeding.

Genome-Wide Identification of Detoxification Genes in Wild Silkworm Antheraea pernyi and Transcriptional Response to Coumaphos

For a half-century, the commercial wild silkworm, Antheraea pernyi, has been protected by coumaphos, which is an internal organophosphorus insecticide used to kill the potential parasitic fly larvae inside. Knowledge about the detoxification genes of A. pernyi as well as the detoxification mechanism for this species remains severely limited. In this study, we identified 281 detoxification genes (32 GSTs, 48 ABCs, 104 CYPs, and 97 COEs) in the genome of this insect, which are unevenly distributed over 46 chromosomes. When compared to the domesticated silkworm, Bombyx mori, a lepidopteran model species, A. pernyi has a similar number of ABCs, but a greater number of GSTs, CYPs, and COEs. By transcriptome-based expression analysis, we found that coumaphos at a safe concentration level significantly changed the pathways related to ATPase complex function and the transporter complex in A. pernyi. KEGG functional enrichment analysis indicated that protein processing in the endoplasmic reticulum was the most affected pathway after coumaphos treatment. Finally, we identified four significantly up-regulated detoxification genes (ABCB1, ABCB3, ABCG11, and ae43) and one significantly down-regulated detoxification gene (CYP6AE9) in response to coumaphos treatment, suggesting that these five genes may contribute to detoxification of coumaphos in A. pernyi. Our study provides the first set of detoxification genes for wild silkworms from Saturniidae and highlights the importance of detoxification gene repertoire in insect pesticide tolerance.

Genome-wide analyses of Glutathione S-transferase gene family and expression profiling under deltamethrin exposure in non-biting midge Propsilocerus akamusi

Glutathione S-transferases (GSTs) are major enzymes in detoxification phase II, and have been functioned in resistance to various insecticides or oxidative stress. Herein, we selected the non-biting midge, Propsilocerus akamusi, widespread in Asian aquatic ecosystems, to uncover the gene location, structure, and phylogenetics relationship of GSTs at genome scale first time. Thirty-three cytosolic and four microsomal GST genes were identified and located on the four chromosomes. The cytosolic GSTs involved in the eight subclasses and five GST genes were unclassified. The expansion of GST genes in P. akamusi experienced duplication events on the delta, theta, xi, iota, and unclassified subclasses. The RNA-Seq analyses and RT-qPCR validation showed that the expression of PaGSTt2 gene is significantly elevated, with deltamethrin concentration increasing. The tertiary structure of PaGSTt2 enzyme was reconstructed, which was different from the other theta gene in the active site. In addition, the GST genes of six chironomids were first described based on the assembled genomes to explore the difference of those in the adaptation to kinds of environments. The GST frame for P. akmusi and its expression profiles provide valuable resources to understand their role in insecticide resistance of this species, as well as those of other biting midges.

Role of Insect and Mammal Glutathione Transferases in Chemoperception

Glutathione transferases (GSTs) are ubiquitous key enzymes with different activities as transferases or isomerases. As key detoxifying enzymes, GSTs are expressed in the chemosensory organs. They fulfill an essential protective role because the chemosensory organs are located in the main entry paths of exogenous compounds within the body. In addition to this protective function, they modulate the perception process by metabolizing exogenous molecules, including tastants and odorants. Chemosensory detection involves the interaction of chemosensory molecules with receptors. GST contributes to signal termination by metabolizing these molecules. By reducing the concentration of chemosensory molecules before receptor binding, GST modulates receptor activation and, therefore, the perception of these molecules. The balance of chemoperception by GSTs has been shown in insects as well as in mammals, although their chemosensory systems are not evolutionarily connected. This review will provide knowledge supporting the involvement of GSTs in chemoperception, describing their localization in these systems as well as their enzymatic capacity toward odorants, sapid molecules, and pheromones in insects and mammals. Their different roles in chemosensory organs will be discussed in light of the evolutionary advantage of the coupling of the detoxification system and chemosensory system through GSTs.

Genome-wide screening and expression of glutathione S-transferase genes reveal that GSTe4 contributes to sensitivity against β-cypermethrin in Zeugodacus cucurbitae

Glutathione S-transferases (GSTs) are an essential multifunctional protein family with common detoxifying enzymes. In this study, 34 GST genes were identified from the melon fly, Zeugodacus cucurbitae, one of the most destructive pests worldwide. These GSTs include 32 cytosolic genes and two microsomal genes. Furthermore, these cytosolic GSTs were classified into six classes: 11 delta, 13 epsilon, three theta, one sigma, two zeta, and two omega. Most of these showed dynamic expression during the developmental stage, some of which showed stage-specific expression. The expression in various adult tissues showed that most of them were expressed in anti-stress-related tissues. The transcriptional response of the delta and epsilon families was determined when Z. cucurbitae was exposed to three insecticides, abamectin, dinotefuran, and β-cypermethrin. Seven genes were significantly up-regulated by abamectin exposure. Moreover, five and four genes were significantly up-regulated with dinotefuran and β-cypermethrin exposure, respectively, demonstrating their involvement in the detoxification of these such toxic substances in Z. cucurbitae. One example of these genes, ZcGSTe4, was randomly selected to explore its function in response to β-cypermethrin exposure. Over-expressed ZcGSTe4 in E. coli showed significant tolerance to β-cypermethrin, and RNAi-mediated suppression of ZcGSTe4 also increased the sensitivity of melon fly to this agent. This study provides a foundation for further studies on the mechanism of detoxification metabolism in the melon fly.

Toxicological and transcriptomic effects in Mythimna separata (Lepidoptera: Noctuidae) exposed to chlorantraniliprole and functional characterization of glutathione S-transferases

BACKGROUND

Chlorantraniliprole (CAP) is an efficient anthranilic diamide insecticide against economically important pests such as the oriental armyworm, Mythimna separata (Lepidoptera: Noctuidae). Resistance to CAP may develop due to enhanced enzymatic detoxification. The glutathione S-transferase (GST) superfamily in M. separata has not been systematically characterized. The aim of this study was therefore to explore the effects of lethal and sublethal doses of CAP on M. separata larvae, screen differentially expressed genes (DEGs) responding to CAP exposure, identify and characterize the GST superfamily, and analyze the metabolism of CAP by recombinant GSTs.

RESULTS

The toxicity bioassay showed that CAP was active against M. separata third-instar larvae. LC50 was 17.615, 3.127, and 1.336 mg/L after 24, 48, and 72 h, respectively. Poisoned larvae showed contracted somites and disrupted midgut. Total GST activity in larvae was significantly elevated 24 h after CAP exposure. RNA-sequencing generated 43 055 unigenes with an average length of 1010 bp, and 567 up-regulated and 692 down-regulated DEGs responding to CAP treatment were screened. Thirty-five GST genes were identified from unigenes, including 31 cytosolic, three microsomal, and one unclassified. The expression profile of GST genes was analyzed using samples from different developmental stages, adult tissues, and CAP treatments. Metabolic assays indicated that CAP was depleted by recombinant MseGSTe2 and MseGSTs6.

CONCLUSIONS

This study provides insight into the toxicological and transcriptomic effects in M. separata larvae exposed to CAP. The identification and functional characterization of the GST superfamily will improve our understanding of CAP detoxification by GSTs. © 2022 Society of Chemical Industry.

Glutathione S-transferase (GST) genes and their function associated with pyrethroid resistance in the malaria vector Anopheles sinensis

BACKGROUND

Glutathione S-transferases (GSTs), a multifunctional protein family, are involved in insecticide resistance. However, a systematic analysis of GSTs in Anopheles sinensis, an important vector for malaria transmission, is lacking. In this study, we investigated the diversity and characteristics of GST genes, and analyzed their expression patterns and functions associated with insecticide resistance in this species.

RESULTS

We identified 32 putative cytosolic and three putative microsomal GST genes in the An. sinensis genome. Transcriptome analysis showed that GSTs were highly expressed in larvae, and mainly expressed in the antennae, midgut and Malpighian tubules of adults. In addition, we found that GSTd2 and GSTe2 were significantly upregulated in four An. sinensis pyrethroid-resistant field populations. Furthermore, silencing of GSTd2 and GSTe2 significantly increased the susceptibility of An. sinensis to deltamethrin, and recombinant GSTd2 and GSTe2 exhibited high enzymatic activity in the metabolism of 1-chloro-2,4-dinitrobenzene and dichlorodiphenyltrichloroethane (DDT).

CONCLUSION

These results showed that GSTs are involved in the development of insecticide resistance in An. sinensis through transcriptional overexpression and enzymatic metabolization, facilitating our understanding of insecticide resistance in insects. © 2022 Society of Chemical Industry.

The Genomics and Population Genomics of the Light Brown Apple Moth, Epiphyas postvittana, an Invasive Tortricid Pest of Horticulture

Simple Summary

In this study, we produced a genomic resource for the light brown apple moth, Epiphyas postvittana, to understand the biological basis of adaptation to a high number of hosts (polyphagy) and the invasive nature of this and other lepidopteran pests. The light brown apple moth is an invasive pest of horticultural plants, with over 500 recorded plant hosts. With origins in Australia, the pest has subsequently spread to New Zealand, Hawaii, California and Europe, causing significant economic losses for fruit producers. Comparative genomic analyses with other lepidopteran genomes indicate that a high proportion of the genome is made up of repetitive sequences, with the majority of the known elements being DNA transposable elements and retrotransposons. Twenty gene families show significant expansions, including some likely to have a role in its pest status. Finally, population genomics, investigated by a RAD-tag approach, indicated likely patterns of invasion and admixture, with Californian moths most probably being derived from Australia.

Abstract

The light brown apple moth, Epiphyas postvittana is an invasive, polyphagous pest of horticultural systems around the world. With origins in Australia, the pest has subsequently spread to New Zealand, Hawaii, California and Europe, where it has been found on over 500 plants, including many horticultural crops. We have produced a genomic resource, to understand the biological basis of the polyphagous and invasive nature of this and other lepidopteran pests. The assembled genome sequence encompassed 598 Mb and has an N50 of 301.17 kb, with a BUSCO completion rate of 97.9%. Epiphyas postvittana has 34% of its assembled genome represented as repetitive sequences, with the majority of the known elements made up of longer DNA transposable elements (14.07 Mb) and retrotransposons (LINE 17.83 Mb). Of the 31,389 predicted genes, 28,714 (91.5%) were assigned to 11,438 orthogroups across the Lepidoptera, of which 945 were specific to E. postvittana. Twenty gene families showed significant expansions in E. postvittana, including some likely to have a role in its pest status, such as cytochrome p450s, glutathione-S-transferases and UDP-glucuronosyltransferases. Finally, using a RAD-tag approach, we investigated the population genomics of this pest, looking at its likely patterns of invasion.

Glutathione S-Transferase Genes are Involved in Lambda-Cyhalothrin Resistance in Cydia pomonella via Sequestration

Pest management is mostly accomplished by the use of insecticides. However, the overuse of insecticides has led to the development of resistance. Glutathione S-transferases (GSTs) are vital detoxification enzymes involved in insecticide resistance in insects. In this study, we report the involvement of GSTs in insecticide resistance to lambda-cyhalothrin in Cydia pomonella, a globally quarantined fruit pest. A total of 25 GST, including 22 cytosolic genes and 3 microsomal genes, are identified from the genome database of C. pomonella. These cytosolic genes are further classified into six classes, including four in delta, eight in epsilon, three in omega, three in sigma, one in theta, and one in zeta class, as well as two unclassified genes. The real-time quantitative polymerase chain reaction (RT-qPCR) shows that the majority of these genes are mainly expressed throughout the larval stage and in the midgut of the fourth-instar larvae. Exposure to an LD₁₀ dose of lambda-cyhalothrin resulted in the upregulation of 17 GST genes. Moreover, mRNA levels of most GST genes, with the exception of CpGSTe6, CpGSTd2, CpGSTd4, and CpGSTz1, are considerably higher in a lambda-cyhalothrin-resistant population (ZW_R) than those of susceptible strains. Recombinant CpGSTd1, CpGSTd3, CpGSTe3, and CpGSTs2 can bind and metabolize lambda-cyhalothrin, with the highest metabolic rate observed for CpGSTd3 but no metabolite(s) was detected, supporting the role of GSTs in sequestration of lambda-cyhalothrin. Molecular dynamics simulation analysis indicates that key residues of hydrophobic pocket-derived lipophilic energy S(lipo) interactions with a hydrophobic pharmacophore of lambda-cyhalothrin are crucial for metabolism by CpGSTd3 and further lead to resistance. Our study is the first to experimentally confirm the involvement of GSTs in lambda-cyhalothrin resistance via sequestration and provides new insights into resistance management in C. pomonella.

SlGSTE9 participates in the stability of chlorpyrifos resistance in Spodoptera litura

BACKGROUND

Spodoptera litura is an important agricultural pest and has developed serious resistance to multiple insecticides. The resistance level to several insecticides is reported to be unstable, but the mechanism is less reported.

RESULTS

Chlorpyrifos and phoxim resistance level in a field-collected population of S. litura declined continuously from the first to the tenth generation and remained stable at the 11th and 12th generations without insecticide exposure. Synergist experiment showed that diethyl maleate and piperonyl butoxide significantly increased mortality to chlorpyrifos and phoxim in the first and sixth generations, but not in the 12th generation. The expression of 31 identified glutathione S-transferase (GST) genes in the third-instar larvae of S. litura in the first, sixth and 12th generations was determined, and eight genes were seen to decrease significantly in the sixth and 12th generations compared with the first generation. SlGSTe9 was selected for further functional study as it had higher abundance and significantly higher expression in the chlorpyrifos-resistant population than in the susceptible population. The recombinant protein of SlGSTE9 showed metabolism activity to chlorpyrifos in vitro and in Escherichia coli, but not to phoxim. Silencing of SlGSTe9 increased the cumulative mortality to chlorpyrifos significantly. SlGSTE9 also showed antioxidant activity to cumene hydroperoxide.

CONCLUSION

Our results suggest that SlGSTe9 is directly involved in chlorpyrifos resistance stability, but not in phoxim. SlGSTE9 may also participate in insecticides resistance by relieving the oxidase stress induced by insecticides. © 2021 Society of Chemical Industry.

Comparative Hessian Fly Larval Transcriptomics Provides Novel Insight into Host and Nonhost Resistance

The Hessian fly is a destructive pest of wheat. Employing additional molecular strategies can complement wheat's native insect resistance. However, this requires functional characterization of Hessian-fly-responsive genes, which is challenging because of wheat genome complexity. The diploid Brachypodium distachyon (Bd) exhibits nonhost resistance to Hessian fly and displays phenotypic/molecular responses intermediate between resistant and susceptible host wheat, offering a surrogate genome for gene characterization. Here, we compared the transcriptomes of Biotype L larvae residing on resistant/susceptible wheat, and nonhost Bd plants. Larvae from susceptible wheat and nonhost Bd plants revealed similar molecular responses that were distinct from avirulent larval responses on resistant wheat. Secreted salivary gland proteins were strongly up-regulated in all larvae. Genes from various biological pathways and molecular processes were up-regulated in larvae from both susceptible wheat and nonhost Bd plants. However, Bd larval expression levels were intermediate between larvae from susceptible and resistant wheat. Most genes were down-regulated or unchanged in avirulent larvae, correlating with their inability to establish feeding sites and dying within 4-5 days after egg-hatch. Decreased gene expression in Bd larvae, compared to ones on susceptible wheat, potentially led to developmentally delayed 2nd-instars, followed by eventually succumbing to nonhost resistance defense mechanisms.

Identification and Characterization of an Antennae-Specific Glutathione S-Transferase From the Indian Meal Moth

Insect glutathione-S-transferases (GSTs) play essential roles in metabolizing endogenous and exogenous compounds. GSTs that are uniquely expressed in antennae are assumed to function as scavengers of pheromones and host volatiles in the odorant detection system. Based on this assumption, antennae-specific GSTs have been identified and functionally characterized in increasing number of insect species. In the present study, 17 putative GSTs were identified from the antennal transcriptomic dataset of the Indian meal moth, Plodia interpunctella, a severe stored-grain pest worldwide. Among the GSTs, only PiGSTd1 is antennae-specific according to both Fragments Per Kilobase Million (FPKM) and quantitative real-time PCR (qRT-PCR) analysis. Sequence analysis revealed that PiGSTd1 has a similar identity as many delta GSTs from other moths. Enzyme kinetic assays using 1-chloro-2,4-dinitrobenzene (CDNB) as substrates showed that the recombinant PiGSTd1 gave a K_m of 0.2292 ± 0.01805 mM and a V_max of 14.02 ± 0.2545 μmol·mg⁻¹·min⁻¹ under the optimal catalytic conditions (35°C and pH = 7.5). Further analysis revealed that the recombinant PiGSTd1 could efficiently degrade the sex pheromone component Z9-12:Ac (75.63 ± 5.52%), as well as aldehyde volatiles, including hexanal (89.10 ± 2.21%), heptanal (63.19 ± 5.36%), (E)-2-octenal (73.58 ± 3.92%), (E)-2-nonenal (75.81 ± 1.90%), and (E)-2-decenal (61.13 ± 5.24%). Taken together, our findings suggest that PiGSTd1 may play essential roles in degrading and inactivating a variety of odorants, especially sex pheromones and host volatiles of P. interpunctella.

Characterization of glutathione S-transferase enzymes in Dictyostelium discoideum suggests a functional role for the GSTA2 isozyme in cell proliferation and development

In this report, we extend our previous characterization of Dictyostelium discoideum glutathione S-transferase (DdGST) enzymes that are expressed in the eukaryotic model organism. Transcript profiling of gstA1-gstA5 (alpha class) genes in vegetative, log phase cells identified gstA2 and gstA3 with highest expression (6-7.5-fold, respectively) when compared to other gstA transcripts. Marked reductions in all gstA transcripts occurred under starvation conditions, with gstA2 and gstA3 exhibiting the largest decreases (-96% and -86.6%, respectively). When compared to their pre-starvation levels, there was also a 60 percent reduction in total GST activity. Glutathione (GSH) pull-down assay and mass spectroscopy detected three isozymes (DdGSTA1, DdGSTA2 and DdGSTA3) that were predominantly expressed in vegetative cells. Biochemical and kinetic comparisons between rDdGSTA2 and rDdGSTA3 shows higher activity of rDdGSTA2 to the CDNB (1-chloro-2,4-dinitrobenzene) substrate. RNAi-mediated knockdown of endogenous DdGSTA2 caused a 60 percent reduction in proliferation, delayed development, and altered morphogenesis of fruiting bodies, whereas overexpression of rDdGSTA2 enzyme had no effect. These findings corroborate previous studies that implicate a role for phase II GST enzymes in cell proliferation, homeostasis, and development in eukaryotic cells.

Genome Mining and Expression Analysis of Carboxylesterase and Glutathione S-Transferase Genes Involved in Insecticide Resistance in Eggplant Shoot and Fruit Borer, Leucinodes orbonalis (Lepidoptera: Crambidae)

The shoot and fruit borer, Leucinodes orbonalis (Lepidoptera: Crambidae) is the major cause of low productivity in eggplant and insecticides being the mainstay of management of L. orbonalis. However, field control failures are widespread due to the evolution of insecticide resistance. Taking advantage of the whole genome sequence information, the present study investigated the level of insecticide resistance and the expression pattern of individual carboxylesterase (CE) and glutathione S-transferases (GSTs) genes in various field collected populations of L. orbonalis. Dose-mortality bioassays revealed a very high level of resistance development against fenvalerate (48.2–160-fold), phosalone (94-534.6-fold), emamectin benzoate (7.2–55-fold), thiodicarb (9.64–22.7-fold), flubendiamide (187.4–303.0-fold), and chlorantraniliprole (1.6–8.6-fold) in field populations as compared to laboratory-reared susceptible iso-female colony (Lo-S). Over-production of detoxification enzymes viz., CE and GST were evident upon enzyme assays. Mining of the draft genome of L. orbonalis yielded large number of genes potentially belonging to the CE and GST gene families with known history of insecticide resistance in other insects. Subsequent RT-qPCR studies on relative contribution of individual genes revealed over-expression of numerous GSTs and few CEs in field populations, indicating their possible involvement of metabolic enzymes in insecticide resistance. The genomic information will facilitate the development of novel resistance management strategies against this pest.

Construction and analysis of the protein-protein interaction network for the olfactory system of the silkworm Bombyx mori

Olfaction plays an essential role in feeding and information exchange in insects. Previous studies on the olfaction of silkworms have provided a wealth of information about genes and proteins, yet, most studies have only focused on a single gene or protein related to the insect's olfaction. The aim of the current study is to determine key proteins in the olfactory system of the silkworm, and further understand protein-protein interactions (PPIs) in the olfactory system of Lepidoptera. To achieve this goal, we integrated Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and network analyses. Furthermore, we selected 585 olfactory-related proteins and constructed a (PPI) network for the olfactory system of the silkworm. Network analysis led to the identification of several key proteins, including GSTz1, LOC733095, BGIBMGA002169-TA, BGIBMGA010939-TA, GSTs2, GSTd2, Or-2, and BGIBMGA013255-TA. A comprehensive evaluation of the proteins showed that glutathione S-transferases (GSTs) had the highest ranking. GSTs also had the highest enrichment levels in GO and KEGG. In conclusion, our analysis showed that key nodes in the biological network had a significant impact on the network, and the key proteins identified via network analysis could serve as new research targets to determine their functions in olfaction.

Investigating the molecular basis of multiple insecticide resistance in a major malaria vector Anopheles funestus (sensu stricto) from Akaka-Remo, Ogun State, Nigeria

Background

Understanding the mechanisms used by Anopheles mosquitoes to survive insecticide exposure is key to manage existing insecticide resistance and develop more suitable insecticide-based malaria vector control interventions as well as other alternative integrated tools. To this regard, the molecular basis of permethrin, DDT and dieldrin resistance in Anopheles funestus (sensu stricto) at Akaka-Remo was investigated.

Methods

Bioassays were conducted on 3–5-day-old adult An. funestus (s.s.) mosquitoes for permethrin, DDT and dieldrin susceptibility test. The molecular mechanisms of mosquito resistance to these insecticides were investigated using microarray and reverse transcriptase PCR techniques. The voltage-gated sodium channel region of mosquitoes was also screened for the presence of knockdown resistance mutations (kdr west and east) by sequencing method.

Results

Anopheles funestus (s.s.) population was resistant to permethrin (mortality rate of 68%), DDT (mortality rate of 10%) and dieldrin (mortality rate of 8%) insecticides. Microarray and RT-PCR analyses revealed the overexpression of glutathione S-transferase genes, cytochrome P450s, esterase, trypsin and cuticle proteins in resistant mosquitoes compared to control. The GSTe2 was the most upregulated detoxification gene in permethrin-resistant (FC = 44.89), DDT-resistant (FC = 57.39) and dieldrin-resistant (FC = 41.10) mosquitoes compared to control population (FC = 22.34). The cytochrome P450 gene, CYP6P9b was also upregulated in both permethrin- and DDT-resistant mosquitoes. The digestive enzyme, trypsin (hydrolytic processes) and the cuticle proteins (inducing cuticle thickening leading to reduced insecticides penetration) also showed high involvement in insecticide resistance, through their overexpression in resistant mosquitoes compared to control. The kdr east and west were absent in all mosquitoes analysed, suggesting their non-involvement in the observed mosquito resistance.

Conclusions

The upregulation of metabolic genes, especially the GSTe2 and trypsin, as well as the cuticle proteins is driving insecticide resistance of An. funestus (s.s.) population. However, additional molecular analyses, including functional metabolic assays of these genes as well as screening for a possible higher cuticular hydrocarbon and lipid contents, and increased procuticle thickness in resistant mosquitoes are needed to further describe their distinct roles in mosquito resistance.

Molecular Cloning and Characterization of Five Glutathione S-Transferase Genes and Promoters from Micromelalopha troglodyta (Graeser) (Lepidoptera: Notodontidae) and Their Response to Tannic Acid Stress

Plants accumulate phenolic compounds such as tannic acid to resist insect herbivores. The survival of insects exposed to toxic secondary metabolites depends on the detoxification metabolism mediated by limited groups of glutathione S-transferases (GSTs). Micromelalopha troglodyta (Graeser) (Lepidoptera: Notodontidae) is an important foliar pest of poplar trees. GSTs play an important role in xenobiotic detoxification in M. troglodyta. Five GST genes were identified in M. troglodyta and were classified into five different cytosolic GST classes, delta, omega, sigma, theta, and zeta. Real-time fluorescent quantitative polymerase chain reaction (qPCR) was used to determine the mRNA expression of the five cloned GSTs in the midguts and fat bodies of M. troglodyta. The mRNA expression of the five GSTs was significantly induced when M. troglodyta was exposed to tannic acid. To further understand the tannic acid regulatory cascade, the 5′-flanking promoter sequences of the five MtGSTs were isolated by genome walking methods, and the promoters were very active and induced by tannic acid. In summary, the induction of GST mRNA expression was due to the response of five MtGST promoters to tannic acid. Therefore, MtGST promoters play an important role in the regulation of GST transcription.

An integrated approach to unravel a crucial structural property required for the function of the insect steroidogenic Halloween protein Noppera-bo

Ecdysteroids are the principal steroid hormones essential for insect development and physiology. In the last 18 years, several enzymes responsible for ecdysteroid biosynthesis encoded by Halloween genes were identified and genetically and biochemically characterized. However, the tertiary structures of these proteins have not yet been characterized. Here, we report the results of an integrated series of in silico, in vitro, and in vivo analyses of the Halloween GST protein Noppera-bo (Nobo). We determined crystal structures of Drosophila melanogaster Nobo (DmNobo) complexed with GSH and 17β-estradiol, a DmNobo inhibitor. 17β-Estradiol almost fully occupied the putative ligand-binding pocket and a prominent hydrogen bond formed between 17β-estradiol and Asp-113 of DmNobo. We found that Asp-113 is essential for 17β-estradiol–mediated inhibition of DmNobo enzymatic activity, as 17β-estradiol did not inhibit and physically interacted less with the D113A DmNobo variant. Asp-113 is highly conserved among Nobo proteins, but not among other GSTs, implying that this residue is important for endogenous Nobo function. Indeed, a homozygous nobo allele with the D113A substitution exhibited embryonic lethality and an undifferentiated cuticle structure, a phenocopy of complete loss-of-function nobo homozygotes. These results suggest that the nobo family of GST proteins has acquired a unique amino acid residue that appears to be essential for binding an endogenous sterol substrate to regulate ecdysteroid biosynthesis. To the best of our knowledge, ours is the first study describing the structural characteristics of insect steroidogenic Halloween proteins. Our findings provide insights relevant for applied entomology to develop insecticides that specifically inhibit ecdysteroid biosynthesis.

Identification and characterization of detoxification genes in two cerambycid beetles, Rhaphuma horsfieldi and Xylotrechus quadripes (Coleoptera: Cerambycidae: Clytini)

The longhorned beetles, Rhaphuma horsfieldi and Xylotrechus quadripes, are two polyphagous insects with larvae feeding on different host plants. In this study, we identified and characterized three gene superfamilies of cytochrome P450s (CYPs), carboxylesterases (COEs) and glutathione-S-transferases (GSTs) involved in the detoxification of endobiotics (e.g., hormones and steroids) and xenobiotics (e.g., insecticides, sex pheromones and plant allelochemicals) through a combination approach of bioinformatics, phylogenetics, expression profiles and genomics. Transcriptome analyses led to the identification of 281 transcripts encoding 135 P450s, 108 COEs and 38 GSTs from the two beetles, coupled with comparative studies of detoxification genes among coleopteran species, suggesting a correlation between host range and the sizes of P450 or COE gene repertoires. The P450s of two beetles were phylogenetically classified into four clades, representing the majority of genes in the CYP3 clan. The COEs from R. horsfieldi and X. quadripes were separately grouped into 11 and 10 clades, and the GST superfamily was assigned into six clades. Expression profiles revealed that the detoxification genes were broadly expressed in various tissues as an implication of functional diversities. Ultimately and more importantly, five alternative splicing events in the Epsilon GSTs, including RhorGSTe7.1/GSTe7.2 and XquaGSTe3.1/GST3.2, were acquired in Coleoptera, in which these genes and their orthologs shared highly conserved gene structure. Our current study has complemented the resources for the detoxification genes in the family Cerambycidae, and allows for functional experiments to identify candidate molecular targets involved in pest resistance to insecticides like organophosphates, organochlorines and pyrethroids.

Divergent molecular evolution in glutathione S-transferase conferring malathion resistance in the oriental fruit fly, Bactrocera dorsalis (Hendel)

Insect glutathione S-transferases (GSTs) are important in insecticide detoxification and Insect-specific GSTs, Epsilon and Delta, have largely expanded in insects. In this study, we functionally expressed and characterized an epsilon class GST gene (BdGSTe8), predominant in the adult Malpighian tubules of Bactrocera dorsalis. This gene may be associated with malathion resistance based on transcriptional studies of resistant and susceptible strains. RNA interference-mediated knockdown of this gene significantly recovered malathion susceptibility in the adults of a malathion-resistant strain, and overexpression of BdGSTe8 enhanced resistance in transgenic Drosophila. Analysis of BdGSTe8 polymorphism showed that several point mutations may be associated with metabolic resistance to malathion. A cytotoxicity assay in Escherichia coli indicated that both of the recombinant BdGSTe8 proteins may play a functional role in protecting cells from toxicity. The allele of BdGSTe8-B conferred higher levels of malathion detoxification capability. Liquid chromatography and ultra-performance liquid chromatography-tandem mass spectrometry analysis showed that the BdGSTe8-A allele did not metabolize malathion directly. However, the BdGSTe8-B allele was involved in the direct metabolism of malathion, which was caused by a mutation in V128A. Further analysis of the sequence suggests that BdGSTe8 evolved rapidly. It maybe play the role of a backup gene and could become a new gene in the future in order to retain the ability of detoxification of malathion, which was driven by positive selection. These results suggest that divergent molecular evolution in BdGSTe8 has played a role in metabolic resistance to malathion in B. dorsalis.

Functional characterization of a novel λ-cyhalothrin metabolizing glutathione S-transferase, CpGSTe3, from the codling moth Cydia pomonella

BACKGROUND

Recent work has shown that two codling moth (Cydia pomonella) glutathione S-transferase genes (GSTs), CpGSTd1 and CpGSTd3, can metabolize λ-cyhalothrin, one of the recommended insecticides for C. pomonella control worldwide. However, systematical characterization of delta and epsilon GSTs, especially their potential contributions in the metabolism of λ-cyhalothrin, is currently still lacking in C. pomonella.

RESULTS

In this study, a total of nine cDNA sequences were identified in C. pomonella, including four in the delta and five in the epsilon subclasses. RT-qPCR showed that seven GSTs were ubiquitously expressed at all developmental stages, and CpGSTe2, CpGSTe3, and CpGSTe4 were mainly expressed in larvae. The mRNA levels of CpGSTd2, CpGSTd4, and CpGSTe5 were significantly higher in male than in female adults. Tissue-specific expression analysis revealed that the CpGSTe2, CpGSTe3, and CpGSTe4 were highly expressed in the midgut while CpGSTd2 and CpGSTd4 were predominantly expressed in the Malpighian tubules. The transcripts of these GSTs (except CpGSTe1) were co-expressed following exposure to LD₁₀ of λ-cyhalothrin for 3 h. Recombinant CpGSTd4, CpGSTe2, and CpGSTe3 proteins expressed in Escherichia coli displayed glutathione-conjugating activity toward 1-chloro-2,4-dinitrobenzene. In addition, λ-cyhalothrin could inhibit the activity of recombinant CpGSTd4, CpGSTe2, and CpGSTe3 enzymes, but only recombinant CpGSTe3 showed λ-cyhalothrin metabolic capacity, with 21.88 ± 1.09% of parental compound being depleted within 1 h.

CONCLUSION

These data show that CpGSTe3 is a third GST gene, encoding an enzyme that metabolizes λ-cyhalothrin in C. pomonella. © 2019 Society of Chemical Industry.

Genomic and transcriptomic analyses of glutathione S-transferases in an endoparasitoid wasp, Pteromalus puparum

Pteromalus puparum is a gregarious pupal endoparasitoid with a wide host range. It deposits eggs into pierid and papilionid butterfly pupae. Glutathione S-transferases (GSTs) are a family of multifunctional detoxification enzymes that act in xenobiotic metabolism in insects. Insect genome projects have facilitated identification and characterization of GST family members. We identified 20 putative GSTs in the P. puparum genome, including 19 cytosolic and one microsomal. Phylogenetic analysis showed that P. puparum GSTs are clustered into Hymenoptera-specific branches. Transcriptomic data of embryos, larvae, female pupae, male pupae, female adults, male adults, venom glands, carcass, salivary glands, and ovaries revealed stage-, sex-, and tissue-specific expression patterns of GSTs in P. puparum. This is the most comprehensive study of genome-wide identification, characterization, and expression profiling of GST family in hymenopterans. Our results provide valuable information for understanding the metabolic adaptation of this wasp.

Comparative transcriptome profiling reveals candidate genes related to insecticide resistance of Glyphodes pyloalis

Glyphodes pyloalis Walker (Lepidoptera: Pyralididae) is a common pest in sericulture and has developed resistance to different insecticides. However, the mechanisms involved in insecticide resistance of G. pyloalis are poorly understood. Here, we present the first whole-transcriptome analysis of differential expression genes in insecticide-resistant and susceptible G. pyloalis. Clustering and enrichment analysis of DEGs revealed several biological pathways and enriched Gene Ontology terms were related to detoxification or insecticide resistance. Genes involved in insecticide metabolic processes, including cytochrome P450, glutathione S-transferases and carboxylesterase, were identified in the larval midgut of G. pyloalis. Among them, CYP324A19, CYP304F17, CYP6AW1, CYP6AB10, GSTs5, and AChE-like were significantly increased after propoxur treatment, while CYP324A19, CCE001c, and AChE-like were significantly induced by phoxim, suggesting that these genes were involved in insecticide metabolism. Furthermore, the sequence variation analysis identified 21 single nucleotide polymorphisms within CYP9A20, CYP6AB47, and CYP6AW1. Our findings reveal many candidate genes related to insecticide resistance of G. pyloalis. These results provide novel insights into insecticide resistance and facilitate the development of insecticides with greater specificity to G. pyloalis.

Effects of phenol on glutathione S-transferase expression and enzyme activity in Chironomus kiiensis larvae

Detoxifying enzyme mRNAs are potentially useful stress biomarkers. Glutathione S-transferase (GST) metabolises lipophilic organic contaminants and mitigates oxidative damage caused by environmental pollutants. Herein, 12 Chironomus kiiensis GSTs (CkGSTs1-6, CkGSTt1-2, CkGSTd1-2, CkGSTm1-2) were cloned and grouped into sigma, theta, delta and microsomal subclasses. Open reading frames (450-699 bp) encode 170-232 amino acid proteins with predicted molecular masses of 17.31-26.84 kDa and isoelectric points from 4.94 to 9.58. All 12 GSTs were expressed during all tested developmental stages, and 11 displayed higher expression in fourth-instar larvae than eggs. GST activity after 24 h of phenol exposure was used to estimate environmental phenol contamination. After exposure to sublethal concentrations of phenol for 48 h, expression and activity of CkGSTs were inhibited in C. kiiensis larvae. Expression of CkGSTd1-2 and CkGSTs1-2 varied with phenol concentration, indicating potential use as biomarkers for monitoring environmental phenol contamination.

High-quality genome assembly of the silkworm, Bombyx mori

In 2008, the genome assembly and gene models for the domestic silkworm, Bombyx mori, were published by a Japanese and Chinese collaboration group. However, the genome assembly contains a non-negligible number of misassembled and gap regions due to the presence of many repetitive sequences within the silkworm genome. The erroneous genome assembly occasionally causes incorrect gene prediction. Here we performed hybrid assembly based on 140 × deep sequencing of long (PacBio) and short (Illumina) reads. The remaining gaps in the initial genome assembly were closed using BAC and Fosmid sequences, giving a new total length of 460.3 Mb, with 30 gap regions and an N50 comprising 16.8 Mb in scaffolds and 12.2 Mb in contigs. More RNA-seq and piRNA-seq reads were mapped on the new genome assembly compared with the previous version, indicating that the new genome assembly covers more transcribed regions, including repetitive elements. We performed gene prediction based on the new genome assembly using available mRNA and protein sequence data. The number of gene models was 16,880 with an N50 of 2154 bp. The new gene models reflected more accurate coding sequences and gene sets than old ones. The proportion of repetitive elements was also reestimated using the new genome assembly, and was calculated to be 46.8% in the silkworm genome. The new genome assembly and gene models are provided in SilkBase (http://silkbase.ab.a.u-tokyo.ac.jp).

Genome-wide identification and expression of the entire 52 glutathione S-transferase (GST) subfamily genes in the Cu2+-exposed marine copepods Tigriopus japonicus and Paracyclopina nana

In this study, the entire glutathione S-transferases (GSTs), the major phase II detoxification enzyme, were identified in two marine copepod species Tigriopus japonicus and Paracyclopina nana. The genome-wide identification of GSTs in T. japonicus and P. nana resulted in 32 and 20 GSTs in total, respectively. Among the identified GSTs, two specific classes of GSTs, specifically sigma and delta/epsilon GSTs were the dominant form of cytosolic GSTs in T. japonicus, while delta/epsilon and mu classes were dominant cytosolic GSTs in P. nana. In addition, Membrane-Associated Proteins in Eicosanoid and Glutathione metabolism (MAPEG) family were found in relatively higher proportion compared to other classes. Moreover, sigma, delta/epsilon, and microsomal GSTs have shown to expand through tandem duplication. To validate the detoxification function of the identified GSTs, both copepods were exposed to copper (Cu²⁺) and the reactive oxygen species (ROS) level and GST activity were measured. With integration of phylogenetic analysis and xenobiotic-mediated GST mRNA expression patterns along with previous enzymatic activities, the functional divergence among species-specific GST genes was clearly observed. This study covers full identification of GST classes in two marine copepod species and their important role in marine environmental ecotoxicology.

CpGSTd3 is a lambda-Cyhalothrin Metabolizing Glutathione S-Transferase from Cydia pomonella (L.)

Little is known about the role of specific delta GST genes in the detoxification of lambda-cyhalothrin in the global quarantine fruit pest codling moth, Cydia pomonella (L.). Real-time quantitative PCR shows that CpGSTd3 was ubiquitously expressed at all developmental stages and is most abundant in the larval stage and lowest in the egg stage; the mRNA level of CpGSTd3 is higher in the midgut and Malpighian tubules of fourth-instar larvae and abdomens of adults than in other tissues. Exposure of fourth-instar larvae to an LD₁₀ dosage of lambda-cyhalothrin significantly induced the transcript of CpGSTd3 at 3 h, but the mRNA level was down-regulated after 12 h of treatment. Recombinant CpGSTd3 expressed in Escherichia coli was able to catalyze the conjugation of 1-chloro-2,4-dinitrobenzene (CDNB) and with an IC₅₀ value of 0.65 mM for lambda-cyhalothrin. Metabolism assays indicate that recombinant CpGSTd3 could metabolize lambda-cyhalothrin. These results suggest that CpGSTd3 is probably a lambda-cyhalothrin metabolizing GST in C. pomonella.

RNA-sequencing of a citrus bud-feeder, Podagricomela weisei (Coleoptera: Chrysomelidae), reveals xenobiotic metabolism/core RNAi machinery-associated genes and conserved miRNAs

The citrus leaf-mining beetle, Podagricomela weisei Heikertinger, is an important citrus pest that ingests the mesophyll and new shoots. The mechanism underlying the xenobiotic metabolism of P. weisei is not well understood, in part because of a lack of available genomic and transcriptomic data, which has hampered the development of novel pest management approaches [e.g., RNA interference (RNAi)]. In this study, we completed the deep sequencing of the P. weisei transcriptome to identify factors potentially involved in xenobiotic metabolism and the core RNAi machinery. The sequencing of the P. weisei transcriptome generated >27 million clean reads, ultimately yielding 90,410 unigenes with an N50 of 1065 bp. The unigenes were used as queries to search the Nr database, which revealed that 21,847 unigenes were homologous to known genes in various species. Transcripts encoding genes involved in xenobiotic metabolism were identified, including genes encoding cytochrome P450 monooxygenase (P450, 47 unigenes), glutathione S-transferase (GST, 12 unigenes), esterase (EST, 25 unigenes), and the ATP-binding cassette transporter (ABC transporter, 32 unigenes). A parallel sequencing of small RNAs detected 30 conserved miRNAs, with the most abundant being Pwe-miR-1-3p, with an expression level reaching 517,996 reads in the prepared library, followed by Pwe-miR-8-3p (149,402 reads). Genes encoding components of the miRNA, siRNA, and piRNA pathways were also identified, and the results indicated that P. weisei possesses only one of each gene in all three pathways. In summary, this is the first detailed analysis of the transcriptome and small RNAs of P. weisei. The datasets presented herein may form the basis for future molecular characterizations of P. weisei as well as the development of enhanced pest control strategies.

Genome-wide analysis of the Tetrahymena thermophila glutathione S-transferase gene superfamily

The ciliate Tetrahymena thermophila has a rapid response to detoxify xenobiotics, which presents opportunity to study the diversification of Glutathione S-Transferase superfamily. In-silico identification of putative GST genes were resulted with 70 GST genes; 49 TtGSTmu, 7 TtGSTomega, 5 TtGSTtheta, 2 TtGSTzeta, 4 TtMAPEG and 3 TtEF1G. TtGST superfamily has short intron carrying or intronless genes. The most expressed mRNAs of TtGST are limited to 4 members at all life stages. TtGST genes are widely distributed to all five micronuclear chromosomes with the highest diversified members from different classes in chromosome 4. The clustering and the orientation of some TtGSTs in the T. thermophila genome give clues about the recent gene duplication. Analysis of GSH affinity-purified GST proteins with Western blot and activity assay showed GST activity carrying purified TtGST populations. In conclusion, the enhanced genome capacity of TtGST superfamily may have evolved through improved GST enzymatic activity.

Evidence of peripheral olfactory impairment in the domestic silkworms: insight from the comparative transcriptome and population genetics

BACKGROUND

The insect olfactory system is a highly specific and sensitive chemical detector, which plays important roles in feeding, mating and finding an appropriate oviposition site. The ecological niche of Bombyx mori has changed greatly since domestication from B. mandarina, and its olfactory response to environmental odorants clearly decreased. However, the mechanisms that result in the olfactory impairment are largely unknown.

RESULTS

The antennal transcriptomes were compared between the domestic and wild silkworms. Comparison of the same sex between the domestic and wild silkworms revealed 1410 and 1173 differentially expressed genes (DEGs) in males and females, respectively. To understand the olfactory impairment, we mainly focused on the olfactory-related genes. In total, 30 olfactory genes and 19 odorant-degrading enzymes (ODEs) showed differential expression in the two comparisons, in which 19 and 14 were down-regulated in the domestic silkworm, respectively. Based on population genomic data, the down-regulated odorant receptors (ORs) showed a higher ratio of unique non-synonymous polymorphisms to synonymous polymorphisms (N/S ratio) in the domestic populations than that in the wild silkworms. Furthermore, one deleterious mutation was found in OR30 of the domestic population, which was located in transmembrane helix 6 (TM6).

CONCLUSIONS

Our results suggested that down-regulation of the olfactory-related genes and relaxed selection might be the major reasons for olfactory impairment of the domestic silkworm reared completely indoor environment. Reversely, wild silkworm may increase expression and remove deleterious polymorphisms of olfactory-related genes to retain sensitive olfaction.

The silk gland damage and the transcriptional response to detoxifying enzymes-related genes of Bombyx mori under phoxim exposure

Silk gland is a major organ of Bombyx mori for the synthesis and secretion of silk protein. Phoxim exposure can be toxic to B. mori and causes a decrease of fibroin synthesis, finally affecting the silk production in industry. To study the mechanism of metabolism and detoxification of silk gland under phoxim exposure, we measured the residual quantity of phoxim in silk gland and hemolymph after phoxim exposure, and the detoxifying enzymes-related genes and enzyme activity were also investigated. Results indicated that the residual amount of phoxim existed up to 24 h in silk gland compared with that in hemolymph, suggesting that phoxim can accumulate in the silk glands within a certain time course. The transcriptional levels of PI3K/Akt genes, including Akt, Tor1, p70s6k and 4e-bp, were up-regulated by 6.919, 1.358, 10.766 and 7.708-fold, respectively. The expression of two downstream genes (CncC and Keap1) was up-regulated by 1.939 and 3.373-fold, respectively. In addition, the transcriptional levels of detoxification-related genes including CYP6AB, CYP306A, CarE2, GST1 and GSTd1 were up-regulated by 1.731, 1.221, 1.366, 1.376 and 6.591-fold, respectively. The enzymatic activity of CYP450, CarE and GST were increased over time. These results provided possible insights into the injury of silk gland and the transcriptional response to detoxifying enzymes-related genes in silkworm after phoxim exposure.

Identification and expression profiles of twenty-six glutathione S-transferase genes from rice weevil, Sitophilus oryzae (Coleoptera: Curculionidae)

The rice weevil, Sitophilus oryzae, is one of the most destructive pests in stored cereals products. In this study, 26 cDNAs encoding glutathione S-transferases (GSTs) were sequenced and characterized in S. oryzae. Phylogenetic analysis displayed the categorization of 26 GSTs into six different cytosolic classes, including two in the delta, twelve in epsilon, three in omega, six in sigma, two in theta, and one in zeta class. RT-qPCR assay illustrated that the relative expression of ten GST genes was significantly higher in adult stages than in larval and pupal developmental stages. Tissue-specific expression analysis revealed that the SoGSTe5, SoGSTe7, SoGSTe12, and SoGSTz1 were up-regulated in the midgut, SoGSTe2, SoGSTe6, and SoGSTs2 were up-regulated in the fat body, and three GSTs (SoGSTd1, SoGSTd2 and SoGSTe4) were up-regulated in Malpighian tubules. RT-qPCR indicated that five GST genes were over expressed after exposure to phosphine at various times and concentrations. The increase in GST gene expressions after phosphine exposure in S. oryzae may lead to an improved tolerance for fumigations and xenobiotics.

Multiple Glutathione S-Transferase Genes in Heortia vitessoides (Lepidoptera: Crambidae): Identification and Expression Patterns

To elucidate the role of glutathione S-transferases (GSTs) in Heortia vitessoides Moore (Lepidoptera: Crambidae), one of the most destructive defoliating pests in Aquilaria sinensis (Lour.) Gilg (Thymelaeaceae) forests, 16 GST cDNAs were identified in the transcriptome of adult H. vitessoides. All cDNAs included a complete open reading frame and were designated HvGSTd1-HvGSTu2. A phylogenetic analysis showed that the 16 HvGSTs were classified into seven different cytosolic classes; three in delta, two in epsilon, three in omega, three in sigma, one in theta, two in zeta, and two in unclassified. The expression patterns of these HvGSTs in various larval and adult tissues, following exposure to half the lethal concentrations (LC50s) of chlorantraniliprole and beta-cypermethrin, were determined using real-time quantitative polymerase chain reaction (RT-qPCR). The expression levels of the 16 HvGSTs were found to differ among various larval and adult tissues. Furthermore, the RT-qPCR confirmed that the transcription levels of nine (HvGSTd1, HvGSTd3, HvGSTe2, HvGSTe3, HvGSTo3, HvGSTs1, HvGSTs3, HvGSTu1, and HvGSTu2) and six (HvGSTd1, HvGSTd3, HvGSTe2, HvGSTo2, HvGSTs1, and HvGSTu1) HvGST genes were significantly higher in the fourth-instar larvae following exposure to the insecticides chlorantraniliprole and beta-cypermethrin, respectively. These genes are potential candidates involved in the detoxification of these two insecticides. Further studies utilizing the RNA interference approach are required to enhance our understanding of the functions of these genes in this forest pest.

RNA interference of two glutathione S-transferase genes, Diaphorina citri DcGSTe2 and DcGSTd1, increases the susceptibility of Asian citrus psyllid (Hemiptera: Liviidae) to the pesticides fenpropathrin and thiamethoxam

BACKGROUND

The Asian citrus psyllid, Diaphorina citri Kuwayama, is an important agricultural pest of citrus globally. Foliar application of chemical insecticides is the most widely used option for reducing D. citri populations. Knockdown of glutathione S-transferase (GST) in several insect species leads to increased susceptibility to insecticides; however, information about the detoxifying role of GST genes in D. citri is unavailable.

RESULTS

Via a sequence homology search, we isolated and characterized three DcGST genes (DcGSTd1, DcGSTe1 and DcGSTe2) from D. citri. Phylogenetic analysis grouped DcGSTd1 into the delta class of GST genes, whereas DcGSTe1 and DcGSTe2 were clustered in the epsilon clade. Gene expression analysis revealed that chlorpyrifos treatment increased the mRNA levels of DcGSTe1 and fenpropathrin enhanced the expression level of DcGSTd1, while DcGSTe2 was significantly up-regulated after exposure to thiamethoxam at a dose of 30% lethal concentration (LC30). RNA interference (RNAi) of DcGSTe2 and DcGSTd1 followed by an insecticide bioassay increased the mortalities of thiamethoxam-treated psyllids by 23.0% and fenpropathrin-treated psyllids by 15.0%. In contrast, knockdown of DcGSTe1 did not significantly increase the susceptibility of D. citri to any of these three insecticides. Further, feeding with double-stranded RNA (dsDcGSTe2-d1) interfusion co-silenced DcGSTe2 and DcGSTd1 expression in D. citri, and led to an increase of susceptibility to both fenpropathrin and thiamethoxam.

CONCLUSION

The findings suggest that DcGSTe2 and DcGSTd1 play unique roles in detoxification of the pesticides thiamethoxam and fenpropathrin. In addition, co-silencing by creating a well-designed dsRNA interfusion against multiple genes was a good RNAi strategy in D. citri. © 2017 Society of Chemical Industry.

Identification of glutathione S-transferases in Bemisia tabaci (Hemiptera: Aleyrodidae) and evidence that GSTd7 helps explain the difference in insecticide susceptibility between B. tabaci Middle East-Minor Asia 1 and Mediterranean

The Bemisia tabaci (Gennadius) (Hemiptera:Aleyrodidae) species complex includes invasive and destructive pests of field crops, and the sibling species MEAM1 and MED are its two most damaging members. Previous research indicated that the replacement of Middle East-Minor Asia 1 (MEAM1) by Mediterranean (MED) as the dominant B. tabaci species in China can be mainly attributed to MED's greater tolerance to insecticides. Glutathione S-transferases (GSTs) play important roles in the detoxification of hydrophobic toxic compounds. To increase our understanding of differences in insecticide resistance between B. tabaci MEAM1 and MED, we searched the genomic and transcriptomic databases and identified 23 putative GSTs in both B. tabaci MEAM1 and MED. Through measuring mRNA levels of 18 of the GSTs after B. tabaci MEAM1 and MED adults were exposed to the insecticide imidacloprid, we found that the expression levels were increased more in B. tabaci MED than in MEAM1 (in particular, the expression level of GST-d7 was increased by 4.39-fold relative to the control). Knockdown of GST-d7 in B. tabaci MED but not in B. tabaci MEAM1 resulted in a substantial increase in the mortality of imidacloprid-treated adults. These results indicate that differences in GST-d7 may help explain why insecticide tolerance is greater in B. tabaci MED than in B. tabaci MEAM1.

The genome of the Hi5 germ cell line from Trichoplusia ni, an agricultural pest and novel model for small RNA biology

We report a draft assembly of the genome of Hi5 cells from the lepidopteran insect pest, Trichoplusia ni, assigning 90.6% of bases to one of 28 chromosomes and predicting 14,037 protein-coding genes. Chemoreception and detoxification gene families reveal T. ni-specific gene expansions that may explain its widespread distribution and rapid adaptation to insecticides. Transcriptome and small RNA data from thorax, ovary, testis, and the germline-derived Hi5 cell line show distinct expression profiles for 295 microRNA- and >393 piRNA-producing loci, as well as 39 genes encoding small RNA pathway proteins. Nearly all of the W chromosome is devoted to piRNA production, and T. ni siRNAs are not 2´-O-methylated. To enable use of Hi5 cells as a model system, we have established genome editing and single-cell cloning protocols. The T. ni genome provides insights into pest control and allows Hi5 cells to become a new tool for studying small RNAs ex vivo.

A common moth called the cabbage looper is becoming increasingly relevant to the scientific community. Its caterpillars are a serious threat to cabbage, broccoli and cauliflower crops, and they have started to resist the pesticides normally used to control them. Moreover, the insect’s germline cells – the ones that will produce sperm and eggs – are used in laboratories as ‘factories’ to artificially produce proteins of interest.

The germline cells also host a group of genetic mechanisms called RNA silencing. One of these processes is known as piRNA, and it protects the genome against ‘jumping genes’. These genetic elements can cause mutations by moving from place to place in the DNA: in germline cells, piRNA suppresses them before the genetic information is transmitted to the next generation. Not all germline cells grow equally well under experimental conditions, or are easy to use to examine piRNA mechanisms in a laboratory. The germline cells from the cabbage looper, on the other hand, have certain characteristics that would make them ideal to study piRNA in insects.

However, the genome of the moth had not yet been fully resolved. This hinders research on new ways of controlling the pest, on how to use the germline cells to produce more useful proteins, or on piRNA.

Decoding a genome requires several steps. First, the entire genetic information is broken in short sections that can then be deciphered. Next, these segments need to be ‘assembled’ – put together, and in the right order, to reconstitute the entire genome. Certain portions of the genome, which are formed of repeats of the same sections, can be difficult to assemble. Finally, the genome must be annotated: the different regions – such as the genes – need to be identified and labeled.

Here, Fu et al. assembled and annotated the genome of the cabbage looper, and in the process developed strategies that could be used for other species with a lot of repeated sequences in their genomes. Having access to the looper’s full genetic information makes it possible to use their germline cells to produce new types of proteins, for example for pharmaceutical purposes. Fu et al. went on to make working with these cells even easier by refining protocols so that modern research techniques, such as the gene-editing technology CRISPR-Cas9, can be used on the looper germline cells.

The mapping of the genome also revealed that the genes involved in removing toxins from the insects’ bodies are rapidly evolving, which may explain why the moths readily become resistant to insecticides. This knowledge could help finding new ways of controlling the pest.

Finally, the genes involved in RNA silencing were labeled: results show that an entire chromosome is the source of piRNAs. Combined with the new protocols developed by Fu et al., this could make cabbage looper germline cells the default option for any research into the piRNA mechanism. How piRNA works in the moth could inform work on human piRNA, as these processes are highly similar across the animal kingdom.

Identification and characterisation of seventeen glutathione S-transferase genes from the cabbage white butterfly Pieris rapae

Insect glutathione S-transferases (GSTs) play essential roles in the detoxification of insecticides and other xenobiotic compounds. The cabbage white butterfly, Pieris rapae, is an economically important agricultural pest. In this study, 17 cDNA sequences encoding putative GSTs were identified in P. rapae. All cDNAs include a complete open reading frame and were designated PrGSTd1-PrGSTz2. Based on phylogenetic analysis, PrGSTs were divided into six classes (delta, epsilon, omega, sigma, theta and zeta). The exon-intron organizations of these PrGSTs were also analysed. Recombinant proteins of eight PrGSTs (PrGSTD1, PrGSTD2, PrGSTE1, PrGSTE2, PrGSTO1, PrGSTS1, PrGSTT1 and PrGSTZ1) were heterologously expressed in Escherichia coli, and all of these proteins displayed glutathione-conjugating activity towards 1-chloro-2,4-dinitrobenzene (CDNB). Expression patterns in various larval tissues, at different life stages, and following exposure to sublethal doses of abamectin, chlorantraniliprole or lambda-cyhalothrin were determined by reverse transcription-quantitative PCR. The results showed that PrGSTe3, PrGSTs1, PrGSTs2, and PrGSTs4 were mainly transcribed in the fat body, while PrGSTe2 was expressed predominantly in the Malpighian tubules. Four genes (PrGSTe2, PrGSTo4, PrGSTs4 and PrGSTt1) were mainly expressed in fourth-instar larvae, while others were ubiquitously expressed in egg, larval, pupa and/or adult stages. Abamectin treatment significantly upregulated ten genes (PrGSTd1, PrGSTd3, PrGSTe1, PrGSTe2, PrGSTo1, PrGSTo3, PrGSTs1, PrGSTs3, PrGSTs4 and PrGSTt1). Chlorantraniliprole and lambda-cyhalothrin treatment significantly upregulated nine genes (PrGSTd1, PrGSTd2, PrGSTe1, PrGSTe2, PrGSTe3, PrGSTs1, PrGSTs3, PrGSTs4 and PrGSTz1) and ten genes (PrGSTd1, PrGSTd3, PrGSTe1, PrGSTe2, PrGSTo1, PrGSTo2, PrGSTs1, PrGSTs2, PrGSTs3 and PrGSTz2), respectively. These GSTs are potentially involved in the detoxification of insecticides.

Identification of glutathione S-transferase genes in Leptinotarsa decemlineata and their expression patterns under stress of three insecticides

Glutathione S-transferases (GSTs) is a family of multifunctional enzymes that are involved in detoxification of poisonous compounds. In the present paper, the Leptinotarsa decemlineata genome and transcriptome dataset were mined and 30 GST genes were identified. These GSTs belonged to cytosolic (29 genes) and microsomal (1 gene) classes. Among them 3 GSTs (LdGSTe2, LdGSTs4, and LdGSTo3) possessed splice variants. Of the 29 cytosolic LdGSTs, 3, 10, 5, 4, 4, and 1 members were classified as delta, epsilon, omega, sigma, theta, and zeta subclasses respectively, along with 2 unclassified genes. Phylogenetic analysis suggest that epsilon, omega and sigma subclasses appear to undergo species-specific bloom. Moreover, most epsilon, omega and sigma GSTs are tandemly arranged in three chromosome scaffolds. To find GST candidates involving in insecticide detoxification, we tested the mRNA levels of 20 GST transcripts under stress of cyhalothrin, fipronil or endosulfan. Out of them, LdGSTe2a, LdGSTe2b, LdGSTo5 and LdGSTt1 were significantly overexpressed after exposure to each of the three insecticides. Two other genes were respectively upregulated after cyhalothrin (LdGSTe10 and LdGSTu2) or endosulfan (LdGSTd1 and LdGSTu2) treatment. The diversified expression responses to insecticide exposure suggest that the LdGSTs may depend on a functionally complex system to detoxify different classes of insecticides. In addition, our findings provide a base for a better understanding of the evolution of insecticide resistance, and functional research on specific GST genes.