Fangmei Zhang1,2
Yijun Chen1,3
Xiaocen Zhao1Shibao Guo1Feng Hong1Yanan Zhi1Li Zhang1Zhou Zhou1Yunhui Zhang2
Xuguo Zhou4
Xiangrui Li2*- 1College of Agriculture, Xinyang Agriculture and Forestry University, Xinyang, China
- 2State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing, China
- 3College of Agriculture, Xinjiang Agricultural University, Urumqi, China
- 4Department of Entomology, University of Kentucky, Lexington, KY, United states
Introduction: Carboxylesterases (CXEs) and glutathione S-transferases (GSTs) can terminate olfactory signals during chemosensation by rapid degradation of odorants in the vicinity of receptors. The tea grey geometrid, Ectropis grisescens (Lepidoptera, Geometridae), one of the most devastating insect herbivores of tea plants in China, relies heavily on plant volatiles to locate the host plants as well as the oviposition sites. However, CXEs and GSTs involved in signal termination and odorant clearance in E. grisescens remains unknown.
Methods: In this study, identification and spatial expression profiles of CXEs and GSTs in this major tea pest were investigated by transcriptomics and qRT-PCR, respectively.
Results: As a result, we identified 28 CXEs and 16 GSTs from female and male antennal transcriptomes. Phylogenetic analyses clustered these candidates into several clades, among which antennal CXEs, mitochondrial and cytosolic CXEs, and delta group GSTs contained genes commonly associated with odorants degradation. Spatial expression profiles showed that most CXEs (26) were expressed in antennae. In comparison, putative GSTs exhibited a diverse expression pattern across different tissues, with one GST expressed specifically in the male antennae.
Disscussion: These combined results suggest that 12 CXEs (EgriCXE1, 2, 4, 6, 8, 18, 20-22, 24, 26, and 29) and 5 GSTs (EgriGST1 and EgriGST delta group) provide a major source of candidate genes for odorants degradation in E. grisescens.
1 Introduction
Insects have a highly specific and sensitive chemosensory system, which is extremely critical for sensing various chemical signals and regulating a series of behaviors (Elgar et al., 2018; Fleischer et al., 2018; Kang et al., 2021). Olfactory perception involves various proteins, including odorant binding proteins (OBPs), chemosensory proteins (CSPs), odorant receptors (ORs), gustatory receptors (GRs), ionotropic receptors (IRs), sensory neuron membrane proteins (SNMPs), and odorant degrading enzymes (ODEs) (Fleischer et al., 2018; Sun et al., 2020; Liu et al., 2021). Briefly, odorant molecules are bound and transported by OBPs onto ORs, and then ORs are activated, leading to signal transduction. Soon after, odor molecules are rapidly degraded by various ODEs to terminate the stimulation, ensuring that insect olfactory sensing systems keep stability and sensitivity for odor identity (Leal, 2013; Suh et al., 2014; Li et al., 2018; He et al., 2019).
ODEs, as multiple enzyme families expressed in the sensillar lymph, include carboxylesterases (CXEs) (Liu et al., 2019; Yi et al., 2021), glutathione S-transferases (GSTs) (Liu et al., 2021; Xia et al., 2022), cytochrome P450 monooxygenases (P450s) (Baldwin et al., 2021; Blomquist et al., 2021), UDP-glucuronosyltransferases (UGTs) (Zhang et al., 2017a), alcohol dehydrogenases (ADs) (Huang et al., 2016), and aldehyde oxidases (AOXs) (Zhang et al., 2017b; Wang et al., 2017). The first identified ODE, ApolPDE, a kind of CXE, could effectively degrade sex pheromone E6Z11-160Ac in Antheraea polyphemus (Vogt and Riddiford, 1981). Antennal-specific GSTs could quickly remove or degrade the odorants from ORs to maintain chemoreceptor sensitivity (Younus et al., 2014; Durand et al., 2018). The antennal-enriched P450s could degrade plant volatiles, insecticides, and pheromones (Cano-Ramirez et al., 2013; Keeling et al., 2013). Antennal ADs have roles in olfaction, which further necessitated investigating its odorant degradation function (Huang et al., 2016). In vitro functional studies clarified the odorant inactivation role of antennal AOXs, such as degrading sex pheromones and plant volatile aldehydes in Amyelois transitella (Choo et al., 2013).
Of these ODEs, CXEs and GSTs are the most well-studied and involved in degrading pheromone/odorant degradation and harmful volatile xenobiotics to maintain the sensitivity of the olfactory receptor neurons (ORNs). So far, many insect CXEs and GSTs have been identified, including Chilo suppressalis (Liu et al., 2015a), Drosophila melanogaster (Chertemps et al., 2015), Ectropis obliqua (Sun et al., 2017), Spodoptera littoralis (Durand et al., 2018), Spodoptera exigua (He et al., 2019), Plodia interpunctella (Liu et al., 2019; Liu et al., 2021), Holotrichia parallela (Yi et al., 2021), and Sitophilus zeamais (Xia et al., 2022), and their functions in insect olfactory perception have been characterized.
CXEs as multifunctional enzymes widely exist in insects, microbes, and plants (Guo and Wong, 2020). In insects, most CXEs are involved in detoxification of exogenous chemicals and are responsible for insecticide resistance and metabolic resistance (Mao et al., 2021). CXEs commonly share conserved active residues, such as the pentapeptide “G-X-S-X-G” (X represents any amino acid), oxyanion hole, glutamate (E), and histidine (H) residues (Godoy et al., 2021). Insect antennal CXEs are characterized as ODEs because they occur in the sensilla and can inactivate odor. They could be divided into ten major clades: mitochondrial and cytosolic esterases, dipteran microsomal α-esterases, cuticular and antennal esterases, β-esterases and pheromone esterases, Lepidopteran juvenile hormone esterases (JHEs), non-Lepidopteran JHEs, moth antennal esterases, neuroligins, neuroreceptors, and gliotactins (Durand et al., 2010a; Oakeshott et al., 2015). To date, many insect antennae CXEs have been identified and functionally characterized for their involvement in the degradation of pheromones or/and plant volatiles. For example, in the genus Spodoptera, two CXEs from S. littoralis (SlCXE7 and 10) (Durand et al., 2010b; Durand et al., 2011) and three from S. exigua (SexiCXE4, 10, and 14) (He et al., 2014a; He et al., 2014b; He et al., 2015) were functionally characterized and degraded both sex pheromones and plant volatiles. A similar phenomenon was also observed in D. melanogaster (EST6) (Chertemps et al., 2012; Chertemps et al., 2015) and Plutella xylostella (PxylCCE16c) (Wang et al., 2021). Additionally, a previous study indicated that CXEs modulated insect mating and foraging behaviors through inactivation of sex pheromones and host volatiles, which made them a novel target for pest behavioral inhibition. For example, the knockdown of GmolCXE1 and 5 had an impact on the mating behaviors of male moths of Grapholita molesta (Wei et al., 2020). Another CXE gene jhedup (duplication of the Juvenile hormone esterase gene) regulated the electrophysiological response and food-seeking ability by hydrolyzing various ester odorants of jhedup mutant D. melanogaster (Steiner et al., 2017).
Similar to CXEs, GSTs also are a diverse family of multifunctional enzymes with conserved amino-terminal domain and carboxyl-terminal domain, which have conserved GSH-binding site (G-site) and hydrophobic substrate (H-site), respectively (Ketterman et al., 2011; Durand et al., 2018). Insect GSTs are divided into six subclasses: delta, epsilon, omega, sigma, theta, and zeta, with some GSTs remaining unclassified, of which the delta- and epsilon-class GSTs are insect-specific (Shi et al., 2012; Labade et al., 2018). Insect GSTs are implicated in the detoxification of xenobiotic compounds (Huang et al., 2011). However, GST expression and their activities demonstrated that they cause signal termination within the olfactory organs (Durand et al., 2018; Tang et al., 2019). For example, in G. molesta, the biochemical characterization of GmolGSTD1 had confirmed its high degradation preference for pheromone component (Z)-8-dodecenyl alcohol and the host plant volatile butyl hexanoate, which showed that it could inactivate odorant molecules and maintain sensitivity to olfactory communication of G. molesta (Li et al., 2018). The antenna-highly expressed CpomGSTd2 in the codling moth, Cydia pomonella, interfered with odorant detection by depredating the odorant (Huang et al., 2017). The antenna-specific SzeaGSTd1 in S. zeamais also played a crucial role in host location by degrading the host volatile, capryl alcohol (Xia et al., 2022).
Ectropis grisescens (Lepidoptera, Geometridae), one of the most destructive pests in tea plantations, causing serious economic losses, is more widely distributed in major tea plantations in China than its counterparts (E. obliqua Prout) (Zhang et al., 2014; Zhang et al., 2016; Li et al., 2019). At present, control of this pest mainly depends on chemical insecticides, leading to environmental pollution and pest resistance (Pan et al., 2021). Moreover, pesticide residues affect the safety of drinking tea (Cao et al., 2018), so using insecticides is forbidden on organic tea plantations (Ma et al., 2016). Therefore, development of novel, effective, and environmentally friendly strategies to control this pest is urgently needed. ODEs thus play an important role in the termination of odorant signals and allow restoration of sensitivity of the olfactory system (Li et al., 2018; He et al., 2019). E. grisescens mainly depends on plant volatiles to search for host plants and locate oviposition sites. Thus, analyses of its mechanism of signal termination and odorant clearance in these important behaviors are necessary. Previous studies considered antennae CXEs or GSTs as potential molecular targets for developing novel pest management strategies based on the manipulation of chemoreception-driven behaviors (Wei et al., 2020; Xia et al., 2022). However, there is no relative report about CXEs or GSTs in E. grisescens yet.
In the present study, we sequenced and analyzed the antennal transcriptome of E. grisescens using Illumina sequencing. Then, CXE and GST gene families and subfamilies were identified and cloned; sequence architecture and phylogenetic analysis were carried out; and finally, quantitative real-time PCR (qRT-PCR) was used to profile their expression patterns in various tissues from both sexes. In addition, the potential roles of the identified CXEs and GSTs in signal termination of olfaction or other physiological processes were discussed. Our study on antennae-specific CXEs and GSTs is particularly important for understanding the molecular mechanism underlying odorant inactivation and subsequent development of ODE-based pest control strategies in E. grisescens.
2 Materials and methods
2.1 Insect rearing and sample collection
Ectropis grisescens larvae were collected from Mount Zhenlei (32°37′N, 114°42′E), Xinyang, Henan, China, and cultivated on fresh tea leaves in the laboratory under the constant conditions of 24°C ± 1°C, 65% relative humidity, and a 16:8 h L: D photoperiod. Emerged adults were fed with 10% honey solution. Antenna, head (without antennae), thorax, abdomen, wing, and leg tissues from 2-day-old unmated male and female insects were dissected, frozen in liquid nitrogen immediately, and stored at −80°C in a refrigerator.
2.2 RNA extraction, cDNA library construction, and sequencing
Total RNA was extracted from male and female antennae (n = 50, three replicates, respectively) using TRIzol reagent (Life Technologies, Carlsbad, CA, United States). The concentration of RNA samples was determined using a NanoDrop ND-1000 spectrophotometer (Thermo Scientific, Wilmington, DE, United States). RNA integrity was assessed using the RNA Nano 6000 Assay Kit of the Bioanalyzer 2100 system (Agilent Technologies, CA, United States). The cDNA libraries were constructed from RNA samples for Illumina sequencing following the Illumina protocol. Sequencing was carried out on the Illumina Novaseq platform (Novogene Co., Ltd., Beijing), and 150-bp paired-end reads were generated.
2.3 Sequence assembly and functional annotation
Raw data (raw reads) of FASTQ format were first modified into clean data (clean reads) through in-house Perl scripts, and clean data (clean reads) were obtained by removing reads containing adapter, reads containing ploy-N, and low-quality reads from raw data. At the same time, Q20, Q30, and GC content of the clean data was calculated. All the downstream analyses were based on the clean data with high quality.
Transcriptome assembly was performed using Trinity with min_kmer_cov set to 2 (Grabherr et al., 2011) by default and all other parameters set to default. Unigene functions were annotated based on NCBI NR, NT, KO, Swiss-Prot, Pfam, GO, and KEGG using Blastx and Blastn searches (E-value < 10−5), retrieving proteins with the highest sequence similarity for each transcript and their protein functional annotations. KEGG Automatic Annotation Server (KASS) was used to search KEGG with E-value = 10−10 (Götz et al., 2008), and Blast2GO v2.5 was used for GO annotation with E-value = 10−6 (Moriya et al., 2007).
Coding sequences (CDSs) were predicted through aligning transcriptome sequences to the NR and Swiss-Prot databases or using ESTScan 3.0.3 (Iseli et al., 1999). The read count for each gene was obtained by mapping clean reads to the assembled transcriptome using RSEM (Bowtie2 parameters: mismatch 0). The read count was calculated as Fragments Per Kilobase of transcript per Million mapped reads (FPKM) (Mortazavi et al., 2008).
2.4 Identification and bioinformatics analyses of candidate CXEs and GSTs
Candidate EgriCXEs and EgriGSTs were identified from the transcriptome data. Furthermore, all candidate EgriCXEs and EgriGSTs were manually checked by the BLASTx program at the NCBI. The complete coding regions were predicted by ORF finder (https://www.ncbi.nlm.nih.gov/orffinder/). Putative signal peptides were predicted with SignalP 5.0 (http://www.cbs.dtu.dk/services/SignalP). The isoelectric point (pI) and molecular weight (Mw) were computed by the ExPASy tool (https://web.expasy.org/compute_pi/). Multiple sequence alignment of the EgriCXE protein was performed using the Clustal program in the Jalview (v2.11.20) software with default parameters (Waterhouse et al., 2009). Phylogenetic trees were constructed in MEGA 11.0 software using the neighbor-joining method with 1000-fold bootstrap resampling (Kumar et al., 2018). The phylogenetic tree image was created by EvolView online software (https://www.evolgenius.info/).
2.5 Tissue expression analysis by qRT-PCR
Total RNA was extracted separately from antennae (n = 50), heads (without antennae; n = 10), thoraxes (n = 5), abdomens (n = 5), wings (n = 5), and legs (n = 10) from both sexes using TRIzol reagent (Life Technologies, Carlsbad, CA, United States). cDNA was synthesized by using TransScript One-Step gDNA Removal and cDNA Synthesis SuperMix (Tiangen, Beijing, China) following the manufacturer’s instructions. Primers specific for EgriCXEs and EgriGSTs were designed using Primer Premier 5.0 software (Premier Biosoft International, Palo Alto, CA, United States; Supplementary Table S1) and synthesized by Sangon Biotech (Shanghai, China).
The reaction volume of 20 µL was prepared using TB Green Premix Ex Taq (Tli RNase H Plus) (TaKaRa, Beijing) by following instructions from the manual. qRT-PCR was conducted using an Applied Biosystems 7500 Fast Real-Time PCR System (Applied Biosystems, Carlsbad, CA) under the following conditions: 95°C for 30 s, then 40 cycles of 95°C for 5 s and 60°C for 34 s, last 95°C for 15 s, 60°C for 1 min, 95°C for 15 s followed by the melting curve analysis. Three biological replicates and three technical replications were carried out. The housekeeping gene glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used as an internal control to normalize the data.
2.6 Statistical analysis
The relative quantification was calculated by the comparative 2−ΔΔCT method (Livak and Schmittgen, 2001). The significance of each candidate EgriCXEs and EgriGSTs among various tissues was determined using a one-way analysis of variance (ANOVA). The significances of EgriCXEs and EgriGSTs from different tissues between female and male adult insects were assessed using a two-sample t-test in SAS statistical software 9.2 (SAS Institute Inc., Cary, NC, United States), with thresholds set at a p < 0.05.
3 Results
3.1 Transcriptome analysis and assembly
After filtering the low-quality and adapter reads, a total of 21,733,434 (97.27%), 23,300,358 (97.54%), 22,180,396 (96.87%), 22,713,765 (96.91%), 22,797,995 (98.15%), and 22,378,252 (97.27%) clean reads were generated from three replicates of female and male antennal libraries of E. grisescens (Supplementary Table S2). The total bases of sequence data were approximately 6.52–6.99 Gb from male and female samples. The average error rates of the sequences were 0.03%. The Q20 and Q30 values for each library exceeded 97% and 92%, respectively, with a GC content of 42.68%–45.79% (Supplementary Table S2). After merging and clustering, the 124,287 transcripts and 52,856 unigenes with a mean length of 1,518 and 1,233 bp, with N50 length of 2,552 and 2,170 bp, were identified, respectively (Supplementary Table S3). The length of transcriptomes and unigenes ranged from 301 to 30,857 bp, with an average length of 1,518 and 1,233 bp, respectively (Figure 1A, Supplementary Table S3).
FIGURE 1. Characterization and transcriptome analysis of RNA sequences from antennae of E. grisescens. (A) Length distributions of the transcriptome and unigenes of E. grisescens. (B) Species distribution from the Blastx results of E. grisescens unigenes in the NR database. (C) GO classifications of E. grisescens unigenes. (D) KEGG classifications of E. grisescens unigenes; a: cellular processes, b: environmental information processing, c: genetic information processing, d: metabolism, e: organismal systems.
3.2 Functional annotation
A total of 15,500 (29.32%), 7,321 (13.85%), 11,512 (21.77%), 14,325 (27.10%), 14,323 (27.09%), and 6,119 (11.57%) had BLASTn hits in the NT, KO, Swiss-Prot, PFAM, GO, and KOG databases, respectively (Supplementary Table S4). BLASTx results showed 19,118 (36.16%) unigenes had the best hits in the non-redundant protein (NR) database. Moreover, most of the annotated unigenes closely matched to Lepidoptera insect sequences (16,780), including Helicoverpa armigera (1,841), Spodoptera litura (1,681), Trichoplusia ni (1,618), Heliothis virescens (1,591), C. suppressalis (1,486), Ostrinia furnacalis (1,293), H. virescens (1,141), Amyelois transitella (817), Bombyx mori (808), and Hyposmocoma kahamanoa (665), as shown in Figure 1B.
Based on the GO annotations, 60,045 unigenes could be annotated into three functional categories: biological processes (50.37%), molecular function (26.98%), and cellular components (22.65%) (Figure 1C). A total of 42 GO terms were identified based on GO level 2, including the odorant recognition process, e.g., binding, catalytic activity, and transporter activity in the molecular function ontology, and localization, signaling, and response to stimulus in the biological process ontology. In the KEGG annotation, 7,936 unigenes were divided into five metabolic pathways: cellular processes, environmental information processing, genetic information processing, metabolism, and organismal systems. Most unigenes were assigned to signal transduction (11.98%), signaling molecules and interaction (1.51%), and environment adaptation (1.11%) involved in recognizing olfaction in insects (Figure 1D).
3.3 Identification of CXEs in E. grisescens
A total of 28 candidate EgriCXEs were identified from the antennal transcriptome of E. grisescens. The sequences were named EgriCXE1–6, 8–9, 11, 13, 15, 17–22, 24–27, 29–30, 32–34, and 36–37 according to their presumptive orthologs of E. obliqua (Sun et al., 2017) and were deposited in the GenBank database under accession numbers OQ296948 to OQ296975 (Table 1). All putative EgriCXEs shared relatively high identities (>60.93%) with their respective orthologs from other species, particularly its counterpart E. obliqua. The amino acid identity between these EgriCXEs ranged from 7.84% to 60.11% (Supplementary Table S5). EgriCXE sequences encoded 401 to 620 amino acid residues with molecular weight ranging from 45.25 to 69.63 kDa and the pI from 5.12 to 9.60. Furthermore, 11 EgriCXEs (EgriCXE1–2, 4–6, 13, 19, 21–22, 29, and 33) were predicted to have putative N-terminal signal peptides (Table 1). Multiple sequence alignment analyses showed that a conserved pentapeptide Gly–X–Ser–X–Gly motif (“X” represents any residue) and more variable oxyanion hole residues (Gly–Gly–Ala), conserved serine (S) residues, and glutamate (E) and histidine (H) residues of the catalytic triad were found (Supplementary Table S6).
TABLE 1. CXE identification and bioinformatics analysis of E. grisescens antennal transcriptomes.
3.4 Identification of GSTs in E. grisescens
A total of 16 candidate EgriGSTs were identified and named EgriGSTe1 to EgriGSTu1, and they were deposited in the GenBank database under the accession numbers OQ296976 to OQ296991 (Table 2). EgriGSTs encoded 179 to 298 amino acid residues, with molecular weight ranging from 20.06 to 33.67 kDa and the pI from 4.95 to 9.46. Blastx search of the best hits showed that all EgriGST sequences shared relatively high sequence identities (53.14%–100.00%) with their respective orthologs from other insects (Table 2). The sequence identities of these EgriGSTs range from 5.50% to 86.76% (Supplementary Table S7). Multiple sequence alignment analyses of the EgriGSTs showed that a conserved G-site can be found in the N-terminal domain and a more variable H-site can be observed with a low sequence identity in the C-terminal domain (Supplementary Figure S1).
TABLE 2. GST identification and bioinformatics analysis of E. grisescens antennal transcriptomes.
3.5 Phylogenetic analysis
The phylogenetic tree of EgriCXEs was constructed with 127 CXE sequences from six kinds of Lepidoptera insects: E. obliqua, S. littoralis, P. interpunctella, S. exigua, S. litura, and Sesamia inferens (Figure 2). The results showed that EgriCXEs were distributed in eight different clades: 1) EgriCXE1, 2, 4, 6, 8, 17–18, and 20–22 were clustered with the moth antennal esterase group; 2) EgriCXE13 and 33 were assigned into β-esterase and pheromone esterase group; 3) EgriCXE5 and 4) EgriCXE32 were distributed into cuticular and antennal esterases and neuroligins, respectively; 5) two EgriCXEs (EgriCXE9 and 37) and 6) four EgriCXEs (EgriCXE 9, 11, 34, and 36) were assigned into neuroreceptor and microsomal α-esterases, respectively; 7) EgriCXE3, 24–27, 29, and 30 were clustered with mitochondrial and cytosolic esterases; 8) EgriCXE15 was assigned into Lepidopteran juvenile hormone esterases.
FIGURE 2. Phylogenetic analysis of candidate EgriCXEs with other insect CXEs. MCEs: mitochondrial and cytosolic esterases; JHEs: Lepidopteran juvenile hormone esterases; NR: neuroreceptor; NL: neuroligin; C/AEs: cuticular and antennal esterases; β-EPEs: β-esterases and pheromone esterases; MAEs: moth antennal esterases; M-α-Es: microsomal α-esterases. Egri, Ectropis grisescens (N = 28); Eobl, Ectropis obliqua (N = 34); Slitt, Spodoptera littoralis (N = 19); Pint, Plodia interpunctella (N = 19); Sexi, Spodoptera exigua (N = 8); Slitu, Spodoptera litura (N = 3); Sinf, Sesamia inferens (N = 16). Candidate EgriCXEs are indicated by red. The GenBank accession numbers of the 127 CXEs protein used in the phylogenetic analysis are listed in Supplementary Table S8.
The phylogenetic tree of EgriGSTs was constructed with 159 GST sequences from eleven insect species: P. interpunctella, P. xylostella, C. pomonella, B. mori, C. suppressalis, Acyrthosiphon pisum, D. melanogaster, Anopheles gambiae, Tribolium castaneum, G. molesta, and S. zeamais (Figure 3). The results showed that the EgriGSTs were divided into six different GST groups: 1) five epsilon EgriGSTs (EgriGSTe1–5); 2) four delta EgriGSTs (EgriGSTd1–4); 3) four sigma EgriGSTs (EgriGSTs1–4); 4) a theta EgriGSTt1; 5) an omega EgriGSTo1, and 6) an unclassified EgriGSTu1.
FIGURE 3. Phylogenetic analysis of candidate EgriGSTs with other insect GSTs. Egri, Ectropis grisescens (N = 16); Pint, Plodia interpunctella (N = 12); Pxyl, Plutella xylostella (N = 19); Cpom, Cydia pomonella (N = 11); Bmor, Bombyx mori (N = 18); Csup, Chilo suppressalis (N = 15); Apis, Acyrthosiphon pisum (N = 14); Dmel, Drosophila melanogaster (N = 20); Agam, Anopheles gambiae (N = 16); Tcas, Tribolium castaneum (N = 12); Gmol, Grapholita molesta (N = 1); Szea, Sitophilus zeamais (N = 1). Candidate EgriGSTs are indicated by red. The GenBank accession numbers of the 159 GSTs protein used in the phylogenetic analysis are listed in Supplementary Table S9.
3.6 Expression profiles of EgriCXEs and EgriGSTs
qRT-PCR results of EgriCXEs (Figure 4) in antennae, heads, thoraxes, abdomens, legs, and wings of both sexes showed that 26 EgriCXEs (EgriCXE1–6, 8–9, 11, 13, 17–22, 24, 26–27, 29–30, 32–34, and 36–37) displayed significant antennal bias expression pattern, except for EgriCXE15 and EgriCXE25. Of them, EgriCXE1–2, 4, 8, 11, 13, 18, 20–22, 24, 26, 29, 32, 34, and 36–37 were significantly highly expressed in female antennae compared to the male antennae (tEgriCXE1 = 17.21, p < 0.0001; tEgriCXE2 = 38.90, p < 0.0001; tEgriCXE4 = 12.52, p = 0.0002; tEgriCXE8 = 7.20, p = 0.0020; tEgriCXE11 = 9.48, p = 0.0007; tEgriCXE13 = 12.77, p = 0.0002; tEgriCXE18 = 4.67, p = 0.0095; tEgriCXE20 = 3.71, p = 0.0207; tEgriCXE21 = 5.06, p = 0.0072; tEgriCXE22 = 3.28, p = 0.0304; tEgriCXE24 = 5.66, p = 0.0048; tEgriCXE26 = 4.67, p = 0.0095; tEgriCXE29 = 19.63, p < 0.0001; tEgriCXE32 = 10.92, p = 0.0004; tEgriCXE34 = 8.03, p = 0.0013; tEgriCXE36 = 7.10, p = 0.0021; and tEgriCXE37 = 8.49, p = 0.0011). Only EgriCXE6 were highly expressed in male antennae, showing male-specific expression (tEgriCXE6 = −8.02, p = 0.0013). Other EgriCXEs (EgriCXE3, 5, 9, 17, 19, 27, 30, and 33) showed no significant difference between the two sexes (tEgriCXE3 = −2.08, p = 0.1061; tEgriCXE5 = 1.98, p = 0.1857; tEgriCXE9 = 0.25, p = 0.8130; tEgriCXE17 = 1.49, p = 0.2115; tEgriCXE19 = −2.19, p = 0.0934; tEgriCXE27 = −0.45, p = 0.6873; tEgriCXE30 = 2.16, p = 0.0966; and tEgriCXE33 = 2.14, p = 0.0989), whereas, EgriCXE15 and EgriCXE25 were highly expressed in a non-chemosensory organ, heads and wings of both sexes, respectively. Furthermore, all EgriCXE genes were detected in other non-chemosensory organs with less expression levels.
FIGURE 4. Tissue- and sex-specific expression pattern of the EgriCXEs. A total of 28 EgriCXEs were clustered using relative expression values from each of the tissues. Expression values are relative to female head tissue (one-fold), with red and blue representing the highest and the lowest values, respectively. FAn: female antennae; FH: female head; FT: female thorax; FAb: female abdomen; FL: female leg; FW: female wing; MAn: male antennae; MH: male head; MT: male thorax; MAb: male abdomen; ML: male leg; MW: male wing.
The putative EgriGSTs showed a wide range of expression patterns (Figure 5): EgriGSTe1, s4, and u1 were highly expressed in male thoraxes and female wings, whereas EgriGSTe2 and d1 were strongly expressed in female thoraxes and male wings. EgriGSTe3 was highly expressed in male thoraxes and female abdomens. EgriGSTe4, s2, t1, and o1 were highly expressed in the wings of both sexes. EgriGSTd2, d3, and d4 were strongly expressed in the heads of both sexes. EgriGSTe5 and s3 were strongly expressed in the thoraxes of both sexes. EgriGSTs1 was specifically expressed in male antennae, with a significant difference between the two sexes (tEgriGSTs1 = −26.43, p = 0.0014). Some EgriGSTs were widely distributed and detected in various tissues, e.g., EgriGSTs2 and EgriGSTs3.
FIGURE 5. Tissue- and sex-specific expression pattern of the EgriGSTs. A total of 16 EgriGSTs were clustered using relative expression values from each of the tissues.
4 Discussion
In the present research, we identified 52,856 unigenes with a mean length of 1,233 bp from the male and female E. grisescens antennal transcriptome, indicating the high quality and great depth of sequencing at the transcriptome level. The results of Blastx homology search in the NCBI database revealed that E. grisescens unigenes shared relatively high homology with sequences from other Lepidoptera species. Ultimately, 28 candidate CXE genes and 16 candidate GST genes were identified in the E. grisescens antennal transcriptome.
In GO annotation of the transcriptome, several annotations were associated with olfaction in insects such as binding, catalytic activity, and transporter activity in the molecular function ontology, localization, signaling, and response to stimulus in the biological process ontology, which are vital steps of odorant recognition in insects (Schmidt and Benton, 2020; Zhang et al., 2020; Rihani et al., 2021; Tian et al., 2021). KEGG pathway analysis also had similar function annotations about recognizing olfaction, such as signal transduction and environment adaptation. The aforementioned results of function annotation were similar to the finding in Bemisia tabaci MED (Wang et al., 2017) and Athetis dissimilis (Song et al., 2021), which further showed that the identified EgriCXEs and EgriGSTs might participate in various chemical communications of E. grisescens.
The number of CXE genes identified in E. grisescens was the same as that in Cnaphalocrocis medinalis (Zhang et al., 2017b) and Athetis lepigone (Zhang et al., 2017c). However, this number was significantly greater than that of CXE genes identified in other species reported, such as C. pomonella (12) (Huang et al., 2016), C. suppressalis (19) (Liu et al., 2015a), H. parallela (20) (Yi et al., 2021), and A. lepigone (20) (Zhang et al., 2017c). The difference in gene numbers among different species might depend on the evolution of divergent behaviors in the long term, which resulted in gene duplication and loss (Zhang et al., 2017b; Huang et al., 2021). In addition, the sample preparation and sequencing method/depth might also be a reason. Multiple sequence alignment analyses showed that most EgriCXEs had the oxyanion hole residues Gly–Gly–Ala, the catalytic triad Ser–Glu–His, and the conserved pentapeptide Gly–X–Ser-–X–Gly. These characteristics indicated that most of the identified EgriCXEs might encode functional enzymes and play a vital role in the catalytic activity of CXEs (Oakeshott et al., 2010).
The number of GST genes identified in E. grisescens was different from that of S. zeamais (13) (Xia et al., 2022), P. interpunctella (17) (Liu et al., 2021), S. littoralis (33) (Durand et al., 2018), and C. medinalis (23) (Liu et al., 2015b). This massive increase in the number of of GSTs in insects enables them to detect plant compounds and resist the damage caused by insecticides, as described in some reports (Durand et al., 2018; Liu et al., 2021). The results of the EgriGST sequence analysis showed a conserved G-site can be found in the N-terminal domain, indicating function as GSH-binding. However, a more variable H-site could be observed with a low sequence identity in the C-terminal domain, which enabled GSTs to accommodate various substrate selectivities (Lerksuthirat and Ketterman, 2008). Furthermore, we also found that some full-length EgriCXEs or EgriGSTs had low identities of amino acid sequences, which suggested that these genes evolved rapidly during long-term adaptation to various environmental factors (Campanini and De Brito, 2016).
The phylogenetic tree showed that the EgriCXEs could be divided into eight groups using classifications from the previous studies (Durand et al., 2010a; Oakeshott et al., 2015). The majority of the antennal EgriCXEs were assigned to the clade that contained members of the esterase family. Among them, ten EgriCXEs (EgriCXE1, 2, 4, 6, 8, 17–18, and 20–22) constituted the biggest groups, moth antennal esterase branch, together with the SexiCXE14 of S. exigua (He et al., 2014c), and SlittCXE7 of S. littoralis (Durand et al., 2011), which caused the degradation of plant volatiles and pheromone compounds. Of them, nine EgriCXEs (EgriCXE1, 2, 4, 6, 8, 18, and 20–22) were significantly expressed in female or male antenna, suggesting these CXEs might have a similar function in odorant degradation. Seven EgriCXEs (EgriCXE3, 24–27, 29, and 30) were clustered into mitochondrial and cytosolic esterase clades with two well-characterized CXEs, SlittCXE10 and SexiCXE10, which were specifically active to plant volatiles pentyl acetate and Z3-6: Ac, respectively (Durand et al., 2010b; He et al., 2015). Of them, three EgriCXEs (EgriCXE 24, 26, and 29) showed high antennal bias, indicating that these EgriCXEs were potentially involved in odorant degradation.
Additionally, two EgriCXEs (EgriCXE13 and 33), along with annotated pheromone and plant volatiles degrading enzymes (like SexiCXE13 for plant volatiles pentyl acetate and the acetate sex pheromone Z9E12-14: Ac, (He et al., 2014a)), were assigned to the β-esterase and pheromone esterase clades. However, the expression pattern of EgriCXE13 and 33 showed no antennal bias in our qRT-PCR. This phenomenon was also observed in the microsomal α-esterase clade, containing four EgriCXEs (EgriCXE9, 11, 34, and 36) without antennal bias and SexiCXE11 which had a high degradation activity against sex pheromones Z9-14: Ac and plant volatile esters pentyl acetate (He et al., 2019). Therefore, functional confirmation of these CXEs would require further functional analyses using in vitro and in vivo methods.
EgriCXE32 constituted Lepidopteran JHEs, which could control juvenile hormone (JH) titer and regulate larval to adult transition in insects (Gu et al., 2015). Two EgriCXEs (EgriCXE19 and 37) and EgriCXE5 were assigned into neuroreceptor and neuroligins clades, respectively, which were mainly involved in neurological and sensory developmental function (Durand et al., 2017). The cuticular and antennal esterase clade only contained one EgriCXE (EgriCXE15) without antennal-biased expression, and this gene lacked CXE conserved characteristic, the oxyanion hole-forming residues (Supplementary Table S6), indicating that it might play other roles, such as detoxification of insecticides.
The phylogenetic tree showed that the EgriGSTs could be divided into six groups. The epsilon class group was common GSTs in insects, which was widely recognized to have a detoxification function (Lalouette et al., 2016; Labade et al., 2018), mediating endocrine plasticity (Bigot et al., 2012) and cholesterol transporter activity (Enya et al., 2015). For instance, SlitGSTe1 and SlitGSTe2 in S. littoralis antennae were induced by sublethal doses of deltamethrin and were involved in protecting against insecticides (Lalouette et al., 2016). The expression pattern suggested EgriGSTe1–5 was highly expressed in thoraxes, abdomens, and wings, which further verified their function as degraders of non-volatile substances.
Additionally, GSTs of the insect delta group in antennae were commonly associated with odorant degradation, which had been studied and verified in some moths. For example, PintGSTd1 of P. interpunctella could efficiently degrade the sex pheromone component Z9-12: Ac and host volatile hexanal (Liu et al., 2021). GmolGSTD1 of G. molesta exhibited high degradation activity to the sex pheromone component (Z)-8-dodecenyl alcohol and the host plant volatile butyl hexanoate (Li et al., 2018). SzeaGSTd1 of S. zeamais could degrade the volatile of the host capryl alcohol (Xia et al., 2022). Unexpectedly, we found that none of the delta genes in E. grisescens was restricted to antenna-based qRT-PCR results, which were mainly expressed in the male heads or ubiquitously expressed in tissues tested. However, phylogenetic tree analysis showed that the delta genes in E. grisescens were distributed among the well-defined insect GST clades, together with GmolGSTD1, PintGSTd1, and SzeaGSTd1. Considering this point, we hypothesized that the delta group mediated the degradation of odorants in E. grisescens, which needed to be further investigated.
The sigma class genes had diverse functions, e.g., recognizing invasive pathogens (Huang et al., 2011) and detoxifying insecticides (Qin et al., 2012; Qin et al., 2013). The expression pattern showed that the sigma group genes in E. grisescens exhibited ubiquitous expression patterns, with the exception of EgriGSTs1, whose expression was almost restricted in the male antennae, suggesting that EgriGSTs1 might play a crucial role in inactivating the chemical signals or pheromone-degrading enzymes (PDEs). Other classes (theta, omega, and unclassified group) were involved in the detoxification of xenobiotics (Qin et al., 2013; Yamamoto and Yamada, 2016; Durand et al., 2018).
In conclusion, we characterized 28 carboxylesterases and 16 glutathione S-transferases encoding odorant-degrading enzymes (ODEs) from E. grisescens antennal transcriptome. Furthermore, the expression patterns of carboxylesterases and glutathione S-transferases in different tissues were investigated to identify antennal enriched genes. Finally, 12 EgriCXEs (EgriCXE1, 2, 4, 6, 8, 18, 20–22, 24, 26, and 29) and five GSTs (EgriGST1 and EgriGST delta group) were identified as candidate target genes involved in odorant degradation of E. grisescens. The findings of this study revealed the potential involvement of ODEs in the olfactory system of E. grisescens.
Data availability statement
The data presented in the study are deposited in the Sequence Read Archive repository, accession number PRJNA784387. The data link: https://www.ncbi.nlm.nih.gov/sra/?term=PRJNA784387.
Author contributions
Conceptualization: FZ and XL; methodology: FH, FZ, and XL; software: YC and FZ; validation: XZ and YZ; investigation: SG; data curation: YC; writing—original draft preparation: YC and FZ; writing—review and editing: XL and XZ; supervision: YZ; project administration: ZZ and LZ; funding acquisition: SG and FZ. All authors have read and agreed to the published version of the manuscript.
Funding
This study was supported by the Natural Science Foundation of Henan Province (No. 212300410229), the Key Project for University Excellent Young Talents of Henan Province (No. 2020GGJS260), the Project of Science and Technology Innovation Team (Nos. XNKJTD-007 and KJCXTD-202001), special funds for Henan Province Scientific and Technological Development guided by the Central Government (Z20221341063), and Youth Foundation (Nos. 2019LG004 and 20200106), Xinyang Agriculture and Forestry University, P.R. China.
Acknowledgments
The authors would like to thank Y Tang and T Zhang (Institute of Plant Protection, Hebei Academy of Agriculture and Forestry Sciences, Baoding, Hebei, China) for their help during data analysis.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer WX declared a shared parent affiliation with the authors YZ and XL to the handling editor at the time of review.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2023.1183610/full#supplementary-material
References
Baldwin, S. R., Mohapatra, P., Nagalla, M., Sindvani, R., Amaya, D., Dickson, H. A., et al. (2021). Identification and characterization of CYPs induced in the Drosophila antenna by exposure to a plant odorant. Sci. Rep. 11 (1), 20530. doi:10.1038/s41598-021-99910-9
Bigot, L., Shaik, H. A., Bozzolan, F., Party, V., Lucas, P., Debernard, S., et al. (2012). Peripheral regulation by ecdysteroids of olfactory responsiveness in male Egyptian cotton leaf worms, Spodoptera littoralis. Insect biochem. Mol. Biol. 42 (1), 22–31. doi:10.1016/j.ibmb.2011.10.003
Blomquist, G. J., Tittiger, C., MacLean, M., and Keeling, C. I. (2021). Cytochromes P450: Terpene detoxification and pheromone production in bark beetles. Curr. Opin. Insect Sci. 43, 97–102. doi:10.1016/j.cois.2020.11.010
Campanini, E. B., and De Brito, R. A. (2016). Molecular evolution of odorant-binding proteins gene family in two closely related Anastrepha fruit flies. BMC Evol. Biol. 16 (1), 198. doi:10.1186/s12862-016-0775-0
Cano-Ramirez, C., Lopez, M. F., Cesar-Ayala, A. K., Pineda-Martinez, V., Sullivan, B. T., and Zúñiga, G. (2013). Isolation and expression of cytochrome P450 genes in the antennae and gut of pine beetle Dendroctonus rhizophagus (Curculionidae: Scolytinae) following exposure to host monoterpenes. Gene 520 (1), 47–63. doi:10.1016/j.gene.2012.11.059
Cao, P., Yang, D. J., Zhu, J. H., Liu, Z. P., Jiang, D. G., and Xu, H. B. (2018). Estimated assessment of cumulative dietary exposure to organophosphorus residues from tea infusion in China. Environ. Health Prev. 23 (1), 7. doi:10.1186/s12199-018-0696-1
Chertemps, T., François, A., Durand, N., Rosell, G., Dekker, T., Lucas, P., et al. (2012). A carboxylesterase, Esterase-6, modulates sensory physiological and behavioral response dynamics to pheromone in Drosophila. BMC Biol. 10, 56. doi:10.1186/1741-7007-10-56
Chertemps, T., Younus, F., Steiner, C., Durand, N., Coppin, C. W., Pandey, G., et al. (2015). An antennal carboxylesterase from Drosophila melanogaster, esterase 6, is a candidate odorant-degrading enzyme toward food odorants. Front. Physiol. 6, 315. doi:10.3389/fphys.2015.00315
Choo, Y. M., Pelletier, J., Atungulu, E., and Leal, W. S. (2013). Identification and characterization of an antennae-specific aldehyde oxidase from the navel orangeworm. PLoS One 8 (6), e67794. doi:10.1371/journal.pone.0067794
Durand, N., Carot-Sans, G., Bozzolan, F., Rosell, G., Siaussat, D., Debernard, S., et al. (2011). Degradation of pheromone and plant volatile components by a same odorant-degrading enzyme in the cotton leafworm, Spodoptera littoralis. PLoS One 6 (12), e29147. doi:10.1371/journal.pone.0029147
Durand, N., Carot-Sans, G., Chertemps, T., Bozzolan, F., Party, V., Renou, M., et al. (2010a). Characterization of an antennal carboxylesterase from the pest moth Spodoptera littoralis degrading a host plant odorant. PLoS One 5 (11), e15026. doi:10.1371/journal.pone.0015026
Durand, N., Carot-Sans, G., Chertemps, T., Montagné, N., Jacquin-Joly, E., Debernard, S., et al. (2010b). A diversity of putative carboxylesterases are expressed in the antennae of the noctuid moth Spodoptera littoralis. Insect Mol. Biol. 19 (1), 87–97. doi:10.1111/j.1365-2583.2009.00939.x
Durand, N., Chertemps, T., Bozzolan, F., and Maïbèche, M. (2017). Expression and modulation of neuroligin and neurexin in the olfactory organ of the cotton leaf worm Spodoptera littoralis. Insect Sci. 24 (2), 210–221. doi:10.1111/1744-7917.12312
Durand, N., Pottier, M. A., Siaussat, D., Bozzolan, F., Maïbèche, M., and Chertemps, T. (2018). Glutathione-S-Transferases in the olfactory organ of the Noctuid moth Spodoptera littoralis, diversity and conservation of chemosensory clades. Front. Physiol. 9, 1283. doi:10.3389/fphys.2018.01283
Elgar, M. A., Zhang, D., Wang, Q., Wittwer, B., Thi Pham, H., Johnson, T. L., et al. (2018). Insect antennal morphology: The evolution of diverse solutions to odorant perception. Yale J. Biol. Med. 91 (14), 457–469.
Enya, S., Daimon, T., Igarashi, F., Kataoka, H., Uchibori, M., Sezutsu, H., et al. (2015). The silkworm glutathione S-transferase gene noppera-bo is required for ecdysteroid biosynthesis and larval development. Insect biochem. Mol. Biol. 61, 1–7. doi:10.1016/j.ibmb.2015.04.001
Fleischer, J., Pregitzer, P., Breer, H., and Krieger, J. (2018). Access to the odor world: Olfactory receptors and their role for signal transduction in insects. Cell Mol. Life Sci. 75 (3), 485–508. doi:10.1007/s00018-017-2627-5
Godoy, R., Machuaca, J., Venthur, H., Quiroz, A., and Mutis, A. (2021). An overview of antennal esterases in Lepidoptera. Front. Physiol. 12, 643281. doi:10.3389/fphys.2021.643281
Götz, S., García-Gómez, J. M., Terol, J., Williams, T. D., Nagaraj, S. H., Nueda, M. J., et al. (2008). High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 36 (10), 3420–3435. doi:10.1093/nar/gkn176
Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29 (7), 644–652. doi:10.1038/nbt.1883
Gu, X., Kumar, S., Kim, E., and Kim, Y. (2015). A whole genome screening and RNA interference identify a juvenile hormone esterase-like gene of the diamondback moth, Plutella xylostella. Plutella Xylostella. J. Insect Physiol. 80, 81–87. doi:10.1016/j.jinsphys.2015.02.001
Guo, S., and Wong, S. M. (2020). A conserved carboxylesterase inhibits tobacco mosaic virus (TMV) accumulation in Nicotiana benthamiana Plants. Viruses 12 (2), 195. doi:10.3390/v12020195
He, P., Li, Z. Q., Liu, C. C., Liu, S. J., and Dong, S. L. (2014a). Two esterases from the genus Spodoptera degrade sex pheromones and plant volatiles. Genome 57 (4), 201–208. doi:10.1139/gen-2014-0041
He, P., Mang, D. Z., Wang, H., Wang, M. M., Ma, Y. F., Wang, J., et al. (2019). Molecular characterization and functional analysis of a novel candidate of cuticle carboxylesterase in Spodoptera exigua degradating sex pheromones and plant volatile esters. Pestic. Biochem. Physiol. 163, 227–234. doi:10.1016/j.pestbp.2019.11.022
He, P., Zhang, J., Li, Z. Q., Zhang, Y. N., Yang, K., Dong, S. L., et al. (2014b). Functional characterization of an antennal esterase from the noctuid moth, Spodoptera exigua. Arch. Insect Biochem. Physiol. 86 (2), 85–99. doi:10.1002/arch.21164
He, P., Zhang, Y. N., Li, Z. Q., Yang, K., Zhu, J. Y., Liu, S. J., et al. (2014c). An antennae-enriched carboxylesterase from Spodoptera exigua displays degradation activity in both plant volatiles and female sex pheromones. Insect Mol. Biol. 23 (4), 475–486. doi:10.1111/imb.12095
He, P., Zhang, Y. N., Yang, K., Li, Z. Q., and Dong, S. L. (2015). An antenna-biased carboxylesterase is specifically active to plant volatiles in Spodoptera exigua. Pestic. Biochem. Physiol. 123, 93–100. doi:10.1016/j.pestbp.2015.03.009
Huang, C., Zhang, X., He, D. F., Wu, Q., Tang, R., Xing, L. S., et al. (2021). Comparative fenomics provide insights into function and evolution of odorant binding proteins in Cydia pomonella. Front. Physiol. 12, 690185. doi:10.3389/fphys.2021.690185
Huang, X. L., Fan, D. S., Liu, L., and Feng, J. N. (2017). Identification and characterization of glutathione S-transferase genes in the antennae of codling moth (Lepidoptera: Tortricidae). Ann. Entomol. Soc. Am. 110 (4), 409–416. doi:10.1093/aesa/sax041
Huang, X. L., Liu, L., Su, X. J., and Feng, J. N. (2016). Identification of biotransformation enzymes in the antennae of codling moth Cydia pomonella. Gene 580 (1), 73–79. doi:10.1016/j.gene.2016.01.008
Huang, Y. F., Xu, Z. B., Lin, X. Y., Feng, Q. L., and Zheng, S. C. (2011). Structure and expression of glutathione S-transferase genes from the midgut of the common cutworm, Spodoptera litura (Noctuidae) and their response to xenobiotic compounds and bacteria. J. Insect Physiol. 57 (7), 1033–1044. doi:10.1016/j.jinsphys.2011.05.001
Iseli, C., Jongeneel, C. V., and Bucher, P. (1999). ESTScan: A program for detecting, evaluating, and reconstructing potential coding regions in EST sequences. Proc. Int. Conf. Intell. Syst. Mol. Biol. 7, 138–148.
Kang, Z. W., Liu, F. H., Xu, Y. Y., Cheng, J. H., Lin, X. L., Jing, X. F., et al. (2021). Identification of candidate odorant-degrading enzyme genes in the antennal transcriptome of Aphidius gifuensis. Entomol. Res. 51 (1), 36–54. doi:10.1111/1748-5967.12489
Keeling, C. I., Henderson, H., Li, M., Dullat, H. K., Ohnishi, T., and Bohlmann, J. (2013). CYP345E2, an antenna-specific cytochrome P450 from the mountain pine beetle, Dendroctonus ponderosae Hopkins, catalyses the oxidation of pine host monoterpene volatiles. Insect biochem. Mol. Biol. 43 (12), 1142–1151. doi:10.1016/j.ibmb.2013.10.001
Ketterman, A. J., Saisawang, C., and Wongsantichon, J. (2011). Insect glutathione transferases. Drug Metab. Rev. 43 (2), 253–265. doi:10.3109/03602532.2011.552911
Kumar, S., Stecher, G., Li, M., Knyaz, C., and Tamura, K. (2018). Mega X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35 (6), 1547–1549. doi:10.1093/molbev/msy096
Labade, C. P., Jadhav, A. R., Ahire, M., Zinjarde, S. S., and Tamhane, V. A. (2018). Role of induced glutathione-S-transferase from Helicoverpa armigera (Lepidoptera: Noctuidae) HaGST-8 in detoxification of pesticides. Ecotoxicol. Environ. Saf. 147, 612–621. doi:10.1016/j.ecoenv.2017.09.028
Lalouette, L., Pottier, M. A., Wycke, M. A., Boitard, C., Bozzolan, F., Maria, A., et al. (2016). Unexpected effects of sublethal doses of insecticide on the peripheral olfactory response and sexual behavior in a pest insect. Environ. Sci. Pollut. Res. Int. 23 (4), 3073–3085. doi:10.1007/s11356-015-5923-3
Leal, W. S. (2013). Odorant reception in insects: Roles of receptors, binding proteins, and degrading enzymes. Annu. Rev. Entomol. 58, 373–391. doi:10.1146/annurev-ento-120811-153635
Lerksuthirat, T., and Ketterman, A. J. (2008). Characterization of putative hydrophobic substrate binding site residues of a Delta class glutathione transferase from Anopheles dirus. Arch. Biochem. Biophys. 479 (1), 97–103. doi:10.1016/j.abb.2008.08.006
Li, G. W., Chen, X. L., Xu, X. L., and Wu, J. X. (2018). Degradation of sex pheromone and plant volatile components by an antennal glutathione S- transferase in the oriental fruit moth, Grapholita molesta Busck (Lepidoptera: Tortricidae). Arch. Insect Biochem. Physiol. 99 (4), e21512. doi:10.1002/arch.21512
Li, Z. Q., Cai, X. M., Luo, Z. X., Bian, L., Xin, Z. J., Liu, Y., et al. (2019). Geographical distribution of Ectropis grisescens (Lepidoptera: Geometridae) and Ectropis obliqua in China and description of an efficient identification method. J. Econ. Entomol. 112 (1), 277–283. doi:10.1093/jee/toy358
Liu, H. M., Tang, Y., Wang, Q. Y., Shi, H. Z., Yin, J., and Li, C. J. (2021). Identification and characterization of an antennae-specific Glutathione S-Transferase from the Indian meal moth. Front. Physiol. 12, 727619. doi:10.3389/fphys.2021.727619
Liu, H. M., Lei, X. P., Du, L. X., Yin, J. Y., Shi, H. Z., Zhang, T., et al. (2019). Antennae-specific carboxylesterase genes from Indian meal moth: Identification, tissue distribution and the response to semiochemicals. J. Stored Prod. Res. 84, 101528. doi:10.1016/j.jspr.2019.101528
Liu, S., Gong, Z. J., Rao, X. J., Li, M. Y., and Li, S. G. (2015a). Identification of putative carboxylesterase and glutathione S-transferase genes from the antennae of the Chilo suppressalis (Lepidoptera: Pyralidae). J. Insect Sci. 15 (1), 103. doi:10.1093/jisesa/iev082
Liu, S., Rao, X. J., Li, M. Y., Feng, M. F., He, M. Z., and Li, S. G. (2015b). Glutathione S-transferase genes in the rice leaffolder, Cnaphalocrocis medinalis (Lepidoptera: Pyralidae): Identification and expression profiles. Arch. Insect Biochem. Physiol. 90 (1), 1–13. doi:10.1002/arch.21240
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 25 (4), 402–408. doi:10.1006/meth.2001.1262
Ma, T., Xiao, Q., Yu, Y. G., Wang, C., Zhu, C. Q., Sun, Z. H., et al. (2016). Analysis of tea Geometrid (Ectropis grisescens) pheromone gland extracts using GC-EAD and GC×GC/TOFMS. J. Agric. Food Chem. 64 (16), 3161–3166. doi:10.1021/acs.jafc.6b00339
Mao, K. K., Ren, Z. J., Li, W. H., Cai, T. W., Qin, X. Y., Wan, H., et al. (2021). Carboxylesterase genes in nitenpyram-resistant Brown planthoppers, Nilaparvata lugens. Insect Sci. 28 (4), 1049–1060. doi:10.1111/1744-7917.12829
Moriya, Y., Itoh, M., Okudam, S., Yoshizawam, A. C., and Kanehisa, M. (2007). Kaas: An automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 35, W182–W185. doi:10.1093/nar/gkm321
Mortazavi, A., Williams, B. A., McCue, K., Schaefer, L., and Wold, B. (2008). Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat. Methods 5 (7), 621–628. doi:10.1038/nmeth.1226
Oakeshott, J. G., Claudianos, C., Campbell, P., Newcomb, R., and Russell, R. J. (2015). “Biochemical genetics and genomics of insect esterases,” in Comprehensive molecular insect science. Editors I. G. Lawrence, I. Kostas, and S. G. Sarjeet (Amsterdam: Elsevier), 309–381.
Oakeshott, J. G., Johnson, R. M., Berenbaum, M. R., Ranson, H., Cristino, A. S., and Claudianos, C. (2010). Metabolic enzymes associated with xenobiotic and chemosensory responses in Nasonia vitripennis. Insect Mol. Biol. 19, 147–163. doi:10.1111/j.1365-2583.2009.00961.x
Pan, Y. J., Fang, G. Q., Wang, Z. B., Cao, Y. H., Liu, Y. J., Li, G. Y., et al. (2021). Chromosome-level genome reference and genome editing of the tea geometrid. Mol. Ecol. Res. 21 (6), 2034–2049. doi:10.1111/1755-0998.13385
Qin, G. H., Jia, M., Liu, T., Zhang, X. Y., Guo, Y. P., Zhu, K. Y., et al. (2013). Characterization and functional analysis of four glutathione S-transferases from the migratory locust, Locusta migratoria. PLoS One 8 (3), e58410. doi:10.1371/journal.pone.0058410
Qin, G. H., Jia, M., Liu, T., Zhang, X. Y., Guo, Y. P., Zhu, K. Y., et al. (2012). Heterologous expression and characterization of a sigma glutathione S-transferase involved in carbaryl detoxification from oriental migratory locust, Locusta migratoria manilensis (Meyen). J. Insect Physiol. 58 (2), 220–227. doi:10.1016/j.jinsphys.2011.10.011
Rihani, K., Ferveur, J. F., and Briand, L. (2021). The 40-year mystery of insect odorant-binding proteins. Biomolecules 11 (4), 509. doi:10.3390/biom11040509
Schmidt, H. R., and Benton, R. (2020). Molecular mechanisms of olfactory detection in insects: Beyond receptors. Open Biol. 10, 200252. doi:10.1098/rsob.200252
Shi, H. X., Pei, L. H., Gu, S. S., Zhu, S. C., Wang, Y. Y., Zhang, Y., et al. (2012). Glutathione S-transferase (GST) genes in the red flour beetle, Tribolium castaneum, and comparative analysis with five additional insects. Genomics 100 (5), 327–335. doi:10.1016/j.ygeno.2012.07.010
Song, Y. Q., Song, Z. Y., Dong, J. F., Lv, Q. H., Chen, Q. X., and Sun, H. Z. (2021). Identification and comparative expression analysis of odorant-binding proteins in the reproductive system and antennae of Athetis dissimilis. Sci. Rep. 11 (1), 13941. doi:10.1038/s41598-021-93423-1
Steiner, C., Bozzolan, F., Montagne, N., Maibeche, M., and Chertemps, T. (2017). Neofunctionalization of "juvenile hormone esterase duplication" in Drosophila as an odorant-degrading enzyme towards food odorants. Sci. Rep. 7 (1), 12629. doi:10.1038/s41598-017-13015-w
Suh, E., Bohbot, J., and Zwiebel, L. J. (2014). Peripheral olfactory signaling in insects. Curr. Opin. Insect Sci. 6, 86–92. doi:10.1016/j.cois.2014.10.006
Sun, D. D., Huang, Y., Qin, Z. J., Zhan, H. X., Zhang, J. P., Liu, Y., et al. (2020). Identification of candidate olfactory genes in the antennal transcriptome of the stink bug Halyomorpha halys. Front. Physiol. 11, 876. doi:10.3389/fphys.2020.00876
Sun, L., Wang, Q., Wang, Q., Zhang, Y. X., Tang, M. J., Guo, H. W., et al. (2017). Identification and expression patterns of putative diversified carboxylesterases in the tea geometrid Ectropis obliqua prout. Front. Physiol. 8, 1085. doi:10.3389/fphys.2017.01085
Tang, Q. F., Shen, C., Zhang, Y., Yang, Z. P., Han, R. R., and Wang, J. (2019). Antennal transcriptome analysis of the maize weevil Sitophilus zeamais: Identification and tissue expression profiling of candidate odorant-binding protein genes. Arch. Insect Biochem. Physiol. 101, e21542. doi:10.1002/arch.21542
Tian, J. H., Zhan, H. X., Dewer, Y., Zhang, B. Y., Qu, C., Luo, C., et al. (2021). Whitefly network analysis reveals gene modules involved in host plant selection, development and evolution. Front. Physiol. 12, 656649. doi:10.3389/fphys.2021.656649
Vogt, R. G., and Riddiford, L. M. (1981). Pheromone binding and inactivation by moth antennae. Nature 293 (5828), 161–163. doi:10.1038/293161a0
Wang, M. M., Long, G. J., Guo, H., Liu, X. Z., Wang, H., Dewer, Y., et al. (2021). Two carboxylesterase genes in Plutella xylostella associated with sex pheromones and plant volatiles degradation. Pest Manag. Sci. 77 (6), 2737–2746. doi:10.1002/ps.6302
Wang, R., Li, F. Q., Zhang, W., Zhang, X. M., Qu, C., Tetreau, G., et al. (2017). Identification and expression profile analysis of odorant binding protein and chemosensory protein genes in Bemisia tabaci MED by head transcriptome. PLoS One 12 (2), e0171739. doi:10.1371/journal.pone.0171739
Waterhouse, A. M., Procter, J. B., Martin, D. M. A., Clamp, M., and Barton, G. J. (2009). Jalview Version 2-a multiple sequence alignment editor and analysis workbench. Bioinformatics 25 (9), 1189–1191. doi:10.1093/bioinformatics/btp033
Wei, H. S., Tan, S. Q., Li, Z., Li, J. C., Moural, T. W., Zhu, F., et al. (2020). Odorant degrading carboxylesterases modulate foraging and mating behaviors of Grapholita molesta. Chemosphere 270, 128647. doi:10.1016/j.chemosphere.2020.128647
Xia, D. S., Zheng, R. W., Huang, J. H., Lu, S. H., and Tang, Q. F. (2022). Identification and functional analysis of Glutathione S-Transferases from Sitophilus zeamais in olfactory organ. Insects 13 (3), 259. doi:10.3390/insects13030259
Yamamoto, K., and Yamada, N. (2016). Identification of a diazinon-metabolizing glutathione S-transferase in the silkworm, Bombyx mori. Sci. Rep. 6, 30073. doi:10.1038/srep30073
Yi, J. K., Wang, S., Wang, Z., Wang, X., Li, G. F., Zhang, X. X., et al. (2021). Identification of candidate carboxylesterases associated with odorant degradation in Holotrichia parallela antennae based on transcriptome analysis. Front. Physiol. 12, 674023. doi:10.3389/fphys.2021.674023
You, Y. C., Xie, M., Ren, N. N., Cheng, X. M., Li, J. Y., Ma, X. L., et al. (2015). Characterization and expression profiling of glutathione S-transferases in the diamondback moth, Plutella xylostella (L.). BMC Genomics 16 (1), 152. doi:10.1186/s12864-015-1343-5
Younus, F., Chertemps, T., Pearce, S. L., Pandey, G., Bozzolan, F., Coppin, C. W., et al. (2014). Identification of candidate odorant degrading gene/enzyme systems in the antennal transcriptome of Drosophila melanogaster. Drosophila melanogaster Insect Biochem. Mol. Biol. 53, 30–43. doi:10.1016/j.ibmb.2014.07.003
Zhang, F. M., Merchant, A., Zhao, Z. B., Zhang, Y. H., Zhang, J., Zhang, Q. W., et al. (2020). Characterization of MaltOBP1, a Minus-C odorant-binding protein, from the Japanese pine sawyer beetle, Monochamus alternatus Hope (Coleoptera: Cerambycidae). Front. Physiol. 11, 212. doi:10.3389/fphys.2020.00212
Zhang, G. H., Yuan, Z. J., Yin, K. S., Fu, J. Y., Tang, M. J., and Xiao, Q. (2016). Asymmetrical reproduce interference between two sibling species of the looper: Ectropis grisescens and Ectropis obliqua. Bull. Enotomol. Re. 11, 1–8. doi:10.1017/S0007485316000602
Zhang, G. H., Yuan, Z. J., Zhang, C. X., Yin, K. S., Tang, M. J., Guo, H. W., et al. (2014). Detecting deep divergence in seventeen populations of tea geometrid (Ectropis obliqua Prout) in China by COI mtDNA and cross-breeding. PloS One 9 (6), e99373. doi:10.1371/journal.pone.0099373
Zhang, Y. N., Li, Z. Q., Zhu, X. Y., Qian, J. L., Dong, Z. P., Xu, L., et al. (2017c). Identification and tissue distribution of carboxylesterase (CXE) genes in Athetis lepigone (Lepidoptera: Noctuidae) by RNA-seq. J. Asia Pac. Entomol. 20 (4), 1150–1155. doi:10.1016/j.aspen.2017.08.016
Zhang, Y. N., Ma, J. F., Xu, L., Dong, Z. P., Xu, J. W., Li, M. Y., et al. (2017a). Identification and expression patterns of UDP-glycosyltransferase (UGT) genes from insect pest Athetis lepigone (Lepidoptera: Noctuidae). J. Asia Pac. Entomol. 20 (1), 253–259. doi:10.1016/j.aspen.2017.01.008
Zhang, Y. X., Wang, W. L., Li, M. Y., Li, S. G., and Liu, S. (2017b). Identification of putative carboxylesterase and aldehyde oxidase genes from the antennae of the rice leaffolder, Cnaphalocrocis medinalis (Lepidoptera: Pyralidae). J. Asia Pac. Entomol. 20 (3), 907–913. doi:10.1016/j.aspen.2017.06.001
Keywords: Ectropis grisescens, antennal transcriptome, odorant-degrading enzyme, carboxylesterases, glutathione S-transferases
Citation: Zhang F, Chen Y, Zhao X, Guo S, Hong F, Zhi Y, Zhang L, Zhou Z, Zhang Y, Zhou X and Li X (2023) Antennal transcriptomic analysis of carboxylesterases and glutathione S-transferases associated with odorant degradation in the tea gray geometrid, Ectropis grisescens (Lepidoptera, Geometridae). Front. Physiol. 14:1183610. doi: 10.3389/fphys.2023.1183610
Received: 10 March 2023; Accepted: 21 March 2023;
Published: 04 April 2023.
Edited by:
Fengqi Li, Guizhou University, ChinaReviewed by:
Wen Xie, Insititute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, ChinaMan Zhao, Henan Agricultural University, China
Copyright © 2023 Zhang, Chen, Zhao, Guo, Hong, Zhi, Zhang, Zhou, Zhang, Zhou and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xiangrui Li, xrli@ippcaas.cn