# Glutathione-s-transferases in the olfactory organ of the noctuid moth spodoptera littoralis, diversity and conservation of chemosensory clades

Nicolas Durand, Marie-Anne Pottier, David Siaussat, Françoise F. Bozzolan, Martine Maïbèche, Thomas Chertemps

<figure><img src="image_1.png" /></figure>

## To cite this version:

Nicolas Durand, Marie-Anne Pottier, David Siaussat, Françoise F. Bozzolan, Martine Maïbèche, et al. Glutathione-s-transferases in the olfactory organ of the noctuid moth spodoptera littoralis, diversity and conservation of chemosensory clades. Frontiers in Physiology, 2018, 9, pp.1-11. 10.3389/fphys.2018.01283. hal-02621342

HAL Id: hal-02621342

https://hal.inrae.fr/hal-02621342

Submitted on 26 May 2020

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Check for
updates

# Glutathione-S-Transferases in the Olfactory Organ of the Noctuid Moth *Spodoptera littoralis*, Diversity and Conservation of Chemosensory Clades

Nicolas Durand, Marie-Anne Pottier, David Siaussat, Françoise Bozzolan, Martine Maïbèche and Thomas Chertemps*

Sorbonne Université, INRA, CNRS, UPEC, IRD, Univ. P7, Institute of Ecology and Environmental Sciences of Paris, Paris, France

Glutathione-S-transferases (GSTs) are conjugating enzymes involved in the detoxification of a wide range of xenobiotic compounds. The expression of GSTs as well as their activities have been also highlighted in the olfactory organs of several species, including insects, where they could play a role in the signal termination and in odorant clearance. Using a transcriptomic approach, we identified 33 putative GSTs expressed in the antennae of the cotton leafworm *Spodoptera littoralis*. We established their expression patterns and revealed four olfactory-enriched genes in adults. In order to investigate the evolution of antennal GST repertoires in moths, we re-annotated antennal transcripts corresponding to GSTs in two moth and one coleopteran species. We performed a large phylogenetic analysis that revealed an unsuspected structural—and potentially functional—diversity of GSTs within the olfactory organ of insects. This led us to identify a conserved clade containing most of the already identified antennal-specific and antennal-enriched GSTs from moths. In addition, for all the sequences from this clade, we were able to identify a signal peptide, which is an unusual structural feature for GSTs. Taken together, these data highlight the diversity and evolution of GSTs in the olfactory organ of a pest species and more generally in the olfactory system of moths, and also the conservation of putative extracellular members across multiple insect orders.

**Keywords:** GST (glutathione S transferase), *Spodoptera littoralis*, olfaction, detoxification, odorant degrading enzyme

## INTRODUCTION

Detoxification is a common process present in nearly all living organisms, from prokaryotes to eukaryotes, allowing the elimination of toxic substances of both exogenous or endogenous origin by sequential enzymatic reactions. The first step of detoxification consists in the introduction of functional groups into lipophilic xenobiotics, mainly by oxido-reduction and hydrolysis reactions performed by phase I enzymes, such as cytochromes P450s (CYPs) and carboxylesterases (CCEs). These phase I metabolites are then conjugated to small hydrophilic molecules by phase II enzymes, a large group of broad-specificity transferases, which in combination can metabolize almost any

**OPEN ACCESS**

**Edited by:**
Monique Gauthier,
Université Toulouse III Paul Sabatier,
France

**Reviewed by:**
Shannon Bryn Olsson,
National Centre for Biological
Sciences, India
Liang Sun,
Tea Research Institute (CAAS), China

*Correspondence:
Thomas Chertemps
thomas.chertemps@
sorbonne-universite.fr

**Specialty section:**
This article was submitted to
Invertebrate Physiology,
a section of the journal
Frontiers in Physiology

Received: 29 June 2018

Accepted: 27 August 2018

Published: 27 September 2018

**Citation:**
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
hydrophobic compound that contains nucleophilic or electrophilic groups (Bock, 2010). Two of the most important classes of this group are the Glutathione-S-transferases (GSTs, EC 2.5.1.18), and the Uridine diphosphate-glycosyltransferases (UGTs, EC 2.4.1.17) (Jakoby and Ziegler, 1990). These two steps of biotransformation led to more hydrophilic metabolites, facilitating their final excretion by various efflux transporters (phase III) such as multidrug resistance transporters (Dermauw and Van Leeuwen, 2014).

Among phase II enzymes, GSTs are highly diversified and play various functions in the detoxification of a wide range of xenobiotic compounds but also in normal cellular metabolism or in oxidative stress response (Li et al., 2007). They catalyze the conjugation of tri-peptide glutathione (GSH) to endogenous electrophilic molecules or to products from the phase I. GSTs are hetero- or homo-dimeric proteins of approximately 25 kDa in size. Each monomer has two domains joined by a variable linker region. The amino terminal domain is highly conserved and provides most of the GSH binding site (G-site) while the carboxyl terminal domain interacts with the hydrophobic substrate (H-site).

In insects, GSTs can be sorted based on their cellular localizations, i.e., mainly cytosolic or microsomal, on their substrate specificities and phylogenetic relationships (Enayati et al., 2005). The cytosolic GSTs are divided in six classes (Delta, Epsilon, Omega, Sigma, Theta, and Zeta) (Sheehan et al., 2001). Despite that they catalyze similar reactions than cytosolic GSTs, microsomal GSTs, now referred as membrane-associated proteins in eicosanoid and glutathione (MAPEG), are very different in origin and structure, as they are mainly trimeric transmembrane proteins (Toba and Aigaki, 2000). Insect GSTs are intensely studied for their role in insecticide resistance (Enayati et al., 2005), in particular, members of the delta and epsilon classes, which are specific to arthropods, have been implicated in resistance to various pesticides. MAPEGs also play a role in xenobiotic detoxification (Zhou et al., 2013) and seem in addition involved in aging process (Toba and Aigaki, 2000). Omega, theta, zeta and microsomal sub-groups appear to be involved in various cellular processes, including protection against oxidative stress (Tu and Akgül, 2005).

Insect GSTs are known to be specifically or preferentially expressed into major detoxification organs such as fat body, midgut but also in epidermis and in Malpighian tubules (Huang et al., 2011). However, the expression of GSTs as well as their activities have been also highlighted in the olfactory organs of several insect species. Antennal expressed GSTs have been indeed identified in various moth species, such as *Manduca sexta* (Rogers et al., 1999), *Helicoverpa armigera* (Wang et al., 2004), *Amyelois transitella* (Leal et al., 2009), *Bombyx mori* (Tan et al., 2014), *Chilo suppressalis* (Liu et al., 2015a), *Epiphyas postvittana* (Corcoran et al., 2015), *Cnaphalocris medinalis* (Liu et al., 2015b), and *Cydia pomonella* (Huang et al., 2017), but also in the fruit fly *Drosophila melanogaster* (Younus et al., 2014) or in the beetles *Agrilus planipennis* (Mamidala et al., 2013), *Dendroctonus valens* (Gu et al., 2015), and *Phyllotreta striolata* (Wu et al., 2016). This particular localization led to the hypothesis of a possible dual function of GSTs in antennae where,

besides their original implication in xenobiotic metabolism, they could play a role in the signal termination and in odorant clearance, as Odorant-Degrading Enzymes (ODEs, Vogt and Riddiford, 1981; Chertemps, 2017). Moreover, the biochemical characterization of an antennal-restricted GST in *M. sexta* (*msexGST-msolf*) confirmed the activity of such enzymes in odorant conjugation and highlighted the possible contribution of GSTs to the detoxification of compounds that might interfere with odorant detection (Rogers et al., 1999).

The cotton leaf worm *Spodoptera littoralis*, is a highly polyphagous crop pest in a broad area of distribution around the Mediterranean basin (Salama et al., 1971). In this pest species, transcriptomic approaches (Legeai et al., 2011; Jacquin-Joly et al., 2012) have already led to the identification of various olfactory gene repertoires, such as Odorant Receptors and Odorant-Binding Proteins, but also of several phase I and II enzymes, such as CCEs (Durand et al., 2010b), CYPs (Pottier et al., 2012), or UGTs (Bozzolan et al., 2014). Here, using more recent and complete transcriptomic data (Poivet et al., 2013), we identified 33 putative GSTs expressed in the antennae of this species. Moreover, their expression patterns were studied using both qualitative and quantitative PCR and revealed four olfactory-enriched genes in *S. littoralis* adults. In order to investigate the evolution and phylogenetic relationships of antennal GST repertoires in moths and the relative conservation of such diversity, we built a phylogenetic analysis based on available sequences from 18 insect species, including 13 lepidopteran species. We first re-annotated antennal transcripts corresponding to GSTs in two other moth species belonging to basal lepidopteran taxa, i.e., the peach fruit moth, *Carposina sasakiii* from the Carposinidae family and the purplish birch-miner moth *Eriocrania semipurpurella*, as a member of the Eriocraniidae. We then annotated GSTs from an antennal transcriptome of a coleopteran species, the european spruce bark beetle *Ips typographus*, to compare with non-lepidopteran species. This large phylogenetic analysis revealed an unsuspected structural—and potentially functional—diversity of GSTs within the olfactory organ. Surprisingly, inside the delta class, we identified a conserved clade containing most of the already identified antennal-specific and antennal-enriched GSTs from moths. In addition, for all the sequences from this clade, we were able to identify a signal peptide, which is an unusual structural feature for GSTs. Taken together, these data highlight the diversity and evolution of GSTs in the olfactory organ of a pest species and more generally in the olfactory system of moths, with in particular the finding of some conserved putative extracellular members across multiple insect orders.

# MATERIALS AND METHODS

## Insects and Tissue Collection

Insects were reared on semi-artificial diet at 23°C, 60-70% relative humidity, and under a 16:8 h light: dark (LD) photoperiod. Adults were kept under an inverted LD regime and provided with a 10% sucrose solution. Male and female antennae and various tissues (proboscis, brain, legs, thorax and abdomen) from 2 day-old males as well as from 7th instar larvae (heads, guts
and carcasses) were dissected and stored at -80°C until RNA extraction.

## RNA Isolation and cDNA Synthesis

Total RNAs were extracted with TRIzol® Reagent (Invitrogen, Carlsbad, CA, United States) and were quantified by spectrophotometry at 260 nm. Single-stranded cDNAs were synthesized from total RNAs (5 µg) from the various tissues using Superscript II reverse transcriptase (Gibco BRL, Invitrogen) and an oligo(dT)18 primer and they were treated with DNase I (Roche, Basel, Switzerland).

## RT-PCR and qRT-PCR

Tissue distribution of *S. littoralis* GSTs was first investigated by RT-PCR. The ubiquitous ribosomal *SlitRpl13* gene, which presents a constant expression in all tissues tested (Durand et al., 2010a), was used as control gene. Primer pairs and PCR conditions are indicated in Table S1. PCR products were loaded on 1% agarose gels and visualized using Gel Red (VWR, Radnor, PE, United States) (Figure S2). Amplification by qPCR of 8 *S. littoralis* GSTs, namely *SlitGSTd2*, *SlitGSTe6*, *SlitGSTe8*, *SlitGSTe9*, *SlitGSTe15*, *SlitGSTs6*, *SlitGSTo3*, *SlitMGST1-3*, and the reference gene *SlitRpl13* was performed as described in detail in Durand et al. (2010a) using the LightCycler® 480 real-time PCR system (Roche). Data were analyzed with LightCycler 480® software (Roche). The crossing point values (C~p~-values) were first determined for the reference genes with a run formed by the 5-fold dilution series, the measuring points, and three negative controls. The normalized *S. littoralis* GST expressions were calculated with Q-Gene software (Simon, 2003) using *SlitRpl13* as reference. This gene has been already demonstrated as the best reference gene in these conditions (Durand et al., 2010b). Each reaction was run in triplicate (technical replicate) with three independent biological replicates. Statistical analyses have been performed with GraphPad Prism® 5 software (ANOVA, post-test Tukey's multiple comparison).

## Identification of Antennal GSTs

Putative partial GST cDNAs were identified from a *de novo* transcriptome of *S. littoralis* (Poivet et al., 2013) by tBLASTn against a dataset of 47 *Spodoptera litura* GST sequences (Zhang et al., 2016) and using known GST genes from insect non-redundant sequence databases (National Center for Biotechnology Information, NCBI). Sequences were completed with the *de novo* transcriptome of *S. littoralis* composed of tissues from various origin, including larvae (Poivet et al., 2013). Antennal enriched sequences where then deduced from antennal specific transcriptomes (Legeai et al., 2011), and confirmed with PCR methods (see above). We named all *S. littoralis* sequences according to the corresponding *S. litura* and *S. frugiperda* sequences (Zhang et al., 2016; Gouin et al., 2017). Subsequently the 7 *S. littoralis* GSTs that were previously published in (Lalouette et al., 2016), namely *SIGSTe1*, *e2*, *e3*, *e4*, *d1*, *d2*, and *d3* were respectively renamed *SlitGSTe15*, *e8*, *e14*, *e12*, *d2*, *d3*, and *SlitMGST1-3*. All sequences have been deposited in Genbank with reference accession number from MH177577 to 177613. GST sequences from other antennal transcriptomes were retrieve

from high quality datasets (regarding length and quality of transcripts and/or number of overall sequences) and selected in various clades of holometabolous insects for comparison. *Ips typographus*, *E. semipurpurella*, and *C. sasakii* were manually annotated from published antennal transcriptomes [PRJEB3262, PRJNA377940, (Yuvaraj et al., 2017), and PRJNA383289 respectively] using the same protocol as *S. littoralis* sequences and named according our phylogenetic analysis (Figure S3).

## Multiple Sequence Alignments and Phylogenetic Analysis

Amino acid sequences were aligned using MAFFT (using L-INS-i option) (Katoh and Standley, 2013) implemented in the Geneious software (http://www.geneious.com, Kearse et al., 2012). Phylogenetic trees were constructed using PhyML (Guindon et al., 2010) based on the LG+G+I+F substitution model as determined by the SMS server (Lefort et al., 2017), using Nearest Neighbor Interchange (NNI). Branch supports were estimated by a Bayesian-like transformation of aLRT (aBayes) (Anisimova et al., 2011). A dendrogram was created and colored using FigTree software (http://tree.bio.ed.ac.uk/software/figtree/). Our final dataset included 415 sequences, including 271 sequences from 13 lepidopteran species for the cytosolic GST sequences, and 15 sequences from 6 lepidopteran species for the MAPEGs.

## Identification of Predicted Signal Peptides in GST Sequences

SignalP4.1 software (Petersen et al., 2011) was used with default D-cutoff value to predict the presence and location of signal peptide (SP) cleavage sites in the GST amino acid sequences. We performed a tBlastn analysis on the nucleotide collection and transcriptome shotgun assembly databases available on the NCBI website. To provide a clear representation of the diversity of GSTs with SP, we aligned a selection of sequences corresponding to the different organism families using MUSCLE (Edgar, 2004) (Figure S4).

## RESULTS

### Identification of *S. littoralis* Antennal GSTs

A total of 37 full-length sequences encoding putative GST proteins from *S. littoralis* (*SlitGSTs*) were identified in our transcriptomic analysis, including 30 cytosolic GSTs and 3 MAPEGs expressed in antennae (Table 1 and Figure 1). Their molecular characteristics such as peptide lengths, estimated molecular masses as well as isoelectric points are indicated in Table S2. As a comparison, we analyzed several available transcriptomes from various holometabolous insects. As shown by the comparison with antennal sequences from other species listed in Table 1, this is the highest number of putative GSTs identified in an insect antenna. An overlook of GSTs in other holometabolous insects, such as Diptera and Coleoptera, confirmed this high number with only 19 sequences in *P. striolata* (Wu et al., 2016), 17 in *I. typographus* (this study) and 31 in *D. melanogaster* (Younus et al., 2014). The high number of antennal GSTs in *S. littoralis* is mainly due to an expansion of
the epsilon clade, with 15 GSTe sequences for only one to six in the other moth species (Table 1), a phenomenon also observed in *D. melanogaster* with 12 GSTe (Younus et al., 2014).

## Phylogeny of Lepidopteran Antennal GSTs

As shown by phylogenetic analysis (Figure 1), all the identified sequences from *S. littoralis* were assigned to the 6 already known insect cytosolic GST clades, except 2 sequences assigned to two distinct unclassified clades (Un.1 and Un.2; Figure 1 and Table S2). In addition, we identified 3 microsomal sequences, that represent the highest number of MAPEGs in antennal transcriptome described so far. In *E. semipurpurella* and *C. sasaki* we also identified respectively 1 and 2 antennal MAPEGs (Table 1). Absence of MAPEG report in the other moth species could probably be due to previous incomplete annotations restricted to cytosolic GSTs.

According to our phylogenetic tree, epsilon, sigma and delta classes represent the major part of the identified sequences in *S. littoralis*, distributed into 5 well supported sub-clades for epsilon GSTs and 2 sub-clades for sigma and delta GSTs. As shown by the short branch length in the phylogenetic tree, antennal GSTz1 sequences are much conserved in moths, even between non-dytrisian and ditrysian species (90.45% of identity between *E. semipurpurella* and *S. littoralis* corresponding sequences). The ratios of non-synonymous to synonymous substitutions were estimated for 10 GSTz1 (Table S3). Their values far <1.0 indicate that these genes are under strong purifying selection pressure, suggesting a functional conservation. The selection pressure for the GSTz2 genes seems not as important with only 67.1% amino acid identity between BmorGSTz2 and SlitGSTz2 (Figure 1).

Amongst SlitGSTs from the delta class, SlitGSTd2 is of particular interest as it falls into a basal clade containing one or two sequences from every species. Moreover, all the delta GSTs from moths known to have an antennal specific expression, such as MsexGST-msolf from *M. sexta* (Robertson et al., 1999; Rogers et al., 1999), BmGSTd4 from *B. mori* (Tan et al., 2014), GST-haolf (Wang et al., 2004) from *H. armigera* and CpomGSTd2 (Huang et al., 2017) from *C. pomonella* fall within this clade.

## Tissue-Related Expression of *S. littoralis* GSTs

Seven adult and larval tissues were tested in order to precise the expression patterns of *SlitGSTs* by RT-PCR and qPCR. Most of them are expressed in all the tissues examined (Figure 2A). *SlitGSTe15* is preferentially expressed in adults whereas *SlitGSTs2* and *SlitGSTs6* are preferentially expressed in larval tissues (Figures 2A,B). *SlitGSTd2*, *e9*, *e6*, and *SlitMGST1-3* seemed preferentially expressed in adult antennae, as confirmed by qPCR on the same tissues but including male legs (Figure 2B). These genes are also expressed at low level in other chemosensory or nervous adult tissues, such as in the proboscis and in the legs for *SlitGSTd2*, in the legs for *SlitGSTe9* and in the brain for *SlitMGST1-3*. All three genes seem more expressed in male antennae than in female antennae but they are not male enriched.

<table><caption>TABLE 1 | Comparison of antennal GSTs repertoire in insect species.</caption><thead><tr><th>Family</th><th>Species</th><th>Delta</th><th>Epsilon</th><th>Zeta</th><th>Theta</th><th>Sigma</th><th>Omega</th><th>Unclassified</th><th>Microsomal</th><th>Total</th><th>Reference</th></tr></thead><tbody><tr><td rowspan="2">Coleoptera</td><td>Phyllostreta striolata</td><td>3(1)</td><td>6</td><td>1</td><td>1</td><td>4</td><td>1</td><td>1</td><td>2(2)</td><td>19</td><td>Wu et al., 2016*</td></tr><tr><td>ips typographus</td><td>1</td><td>5</td><td>-</td><td>1</td><td>5</td><td>2</td><td>1</td><td>2</td><td>17</td><td>This study</td></tr><tr><td>Drosophilidae</td><td>Drosophila melanogaster</td><td>8(2)</td><td>12(2)</td><td>2</td><td>4</td><td>1</td><td>3</td><td>1</td><td>-</td><td>31</td><td>Younus et al., 2014</td></tr><tr><td rowspan="2">Eriocraniidae</td><td>Eriocrania semipurpurella</td><td>4</td><td>4</td><td>2</td><td>1</td><td>2</td><td>3</td><td>2</td><td>1</td><td>19</td><td>This study</td></tr><tr><td>Carposina sasaki</td><td>3</td><td>3</td><td>2</td><td>1</td><td>4</td><td>4</td><td>1</td><td>2</td><td>20</td><td>This study</td></tr><tr><td>Tortricidae</td><td>Cydia pomonella</td><td>1(1)</td><td>1</td><td>1</td><td>1</td><td>1</td><td>2</td><td>1</td><td>2</td><td>10</td><td>Huang et al., 2017</td></tr><tr><td rowspan="2">Crambidae</td><td>Epiphyas postvittana</td><td>3</td><td>5</td><td>-</td><td>1</td><td>6</td><td>3</td><td>2</td><td>-</td><td>20</td><td>Corcoran et al., 2015</td></tr><tr><td>Chilo suppressalis</td><td>4</td><td>3</td><td>1</td><td>1</td><td>2</td><td>4</td><td>1</td><td>-</td><td>16</td><td>Liu et al., 2015a</td></tr><tr><td rowspan="2">Noctuidae</td><td>Cnaphatoecrosis medinalis</td><td>4</td><td>5</td><td>2</td><td>-</td><td>5</td><td>3</td><td>2</td><td>-</td><td>21</td><td>Liu et al., 2015b</td></tr><tr><td>Spodoptera littoralis</td><td>3(1)</td><td>15(2)</td><td>2</td><td>1</td><td>4</td><td>3</td><td>2</td><td>3(1)</td><td>33</td><td>This study</td></tr></tbody></table>

Numbers of adult antennal enriched GSTs are indicated between parentheses. *Corrected number. Original number was 27 but 9 sequences are not insect GSTs and were removed from the analysis (Table S5).
<figure><img src="image_4.png" /></figure>

## Identification of Predicted Signal Peptide in Insect GST Sequences

As seen previously, *SlitGStd2* clusters in a conserved subgroup which includes *BmGSTD4* and *CpomGSTD2* for which predicted signal peptide (SP) sequences have been identified (Tan et al., 2014; Huang et al., 2017). *SlitGStd2* possesses also a predicted 22 amino acid SP at the N-terminus of the protein sequence, and its presence was confirmed by using RT-PCR with specific primers (Figure S1). Screening of the 36 remaining S. littoralis GST sequences confirmed that *SlitGStd2* was the only sequence presenting this feature.

To test whether this characteristic was shared with the other members of the subclade, we first searched for SP in the sequence of *AtraGST* (Leal et al., 2009), *MsexGST-msolf* (Robertson et al., 1999; Rogers et al., 1999), *CsupGSTD1* (Liu et al., 2015a) and *CmedGSTD1* (Liu et al., 2015b) and we indeed identified putative SPs for all them. We then completed by bioinformatics the N-terminus sequences for the 5 other genes from this clade (*SfruGSTD2*, *HarmGSTD02*, *BmGSTD1*, *EposGST11*, and *PxylGSTD1*) and we also identified predicted SP for each of them. All the members of the *SlitGStd2* cluster thus contain a predicted SP.
<figure><img src="image_5.png" /><figcaption>**FIGURE 2** | Analysis of *Spodoptera littoralis* GST expression throughout the body. (A) RT-PCR products visualized on agarose gel. (B) Expression of *SlitGSTd2*, *SlitGSTe6*, *SlitGSTe9*, *SlitGSTe15*, *SlitGSTs6*, and *SlitMGST1-3* by qPCR. Data were obtained from triplicate experiments and are given as the Mean ± SD. Templates were from adult female antennae, adult male antennae, proboscis, brains, abdomen and larval heads, guts and carcasses. *SlitRpl13* was used as the control gene.</figcaption></figure>

We then searched the public databases and the literature for insect GSTs possessing this unusual feature and identified predicted SP in 67 full-length GST sequences (Table S4). Putative SP is present in the sequences of GSTs from the delta clade in 38 lepidopteran species, including 6 non-dytrisian species, in 6 hemipteran and 9 phasmopteran species. SP was also isolated in a sigma GST from a Coleoptera, the mountain pine beetle *Dendroctonus ponderosae*. Finally, we found 12 sigma and omega GSTs with SP in nematodes, another group of Ecdysozoa. Alignment of SP sequences from some insect and nematode GSTs (Figure 3) revealed that these signal peptides are 15–26 amino acid long and generally end with one to three alanine residues as amino acids with hydrophobic side chain. In Lepidoptera, the amino acids around the cleavage site are well conserved with a N-A-A/X-A followed by a RSK motif.

DISCUSSION

Diversity of GSTs in S. littoralis Antenna

In the cotton leafworm *S. littoralis*, an unsuspected number
of CCEs, UGTs, and CYPs had been previously described in

antenna (Durand et al., 2010b; Pottier et al., 2012; Bozzolan
et al., 2014), leading to the hypothesis that antennal enzymes
could participate in signal inactivation and odorant clearance
as Odorant-Degrading Enzymes, but also in detoxification
processes. Indeed, various airborne compounds, such as toxic
allelochemicals emitted by plants or anthropic xenobiotics could
enter the olfactory sensilla and reach the olfactory receptor
neurons (ORNs) and potentially harm them (Siaussat et al.,
2014). In this study, we demonstrate an unsuspected diversity of
GSTs in *S. littoralis* antenna. With 33 described GSTs, the number
of antennal-expressed genes is more abundant than any other
Lepidopteran described so far.

Without any genome-scale information, the total number of
GSTs from *S. littoralis* cannot be predicted. However, with the
recent completion of the genome of *S. litura*, a sister species
(47 GSTs annotated, Cheng et al., 2017), and from a closely
related species *S. frugiperda* (46 GSTs, Gouin et al., 2017), we can
speculate that at least two third of total GSTs could be present in
*S. littoralis* antennae. This massive expansion, driven by multiple
gene duplications and polymorphism, is suspected to be part of
*S. littoralis* ability to detect and detoxify many plant compounds,
<figure><img src="image_6.png" /><figcaption>**FIGURE 3** | Partial alignment of 35 GSTs from insects and nematodes showing several sequences from delta, sigma and omega classes with predicted signal peptides. Each sequence is representative of one family. Slit, *Spodoptera littoralis*; Bsup, *Biston suppressaria*; Msex, *Manduca sexta*; Atra, *Amyelois transitella*; Bany, *Bicyclus anynana*; Ehip, *Eogystia hippophaecolus*; Ztau, *Zygaena fausta*; Cpom, *Cydia pomonella*; Pxyl, *Plutella xylostella*; Tbis, *Tineola bisselliella*; Tque, *Tischeria quercitella*; Alam, *Andesiana lamellata*; Psp, *Ptyssoptera* sp.; Esem, *Eriocrania semipurpurea*; Llin, *Lygus lineolaris*; Meup, *Macrosiphum euphorbiae*; Ssip, *Sipyloidea* sipylus; Pphi, *Phyllium philippinicum*; Aasp, *Aretaon asperrimus*; Psch, *Peruphasma schultei*; Dpon, *Dendroctonus ponderosae*; Ovol, *Onchocerca volvulus*; Tcan, *Toxocara canis*; Agal, *Ascaridia galli*; Bxyl, *Bursaphelenchus xylophilus*. BxylGSTs full name is BUX.s000647.112 (Espada et al., 2016a). Identifiers in bold indicate sequences linked to representative publications. Pful, *Phragmatobia fulliginosa*; Bm, *Bombay mori*; Sric, *Samia ricini*; Lsti, *Loxostege sticticalis*; Ppol, *Papilio polytes*; Acon, *Argyresthia conjugella*; Etia, *Extatosoma tiaratum*.</figcaption></figure>

enabling this species to cope with numerous host plants and may contribute to insecticide resistance (Elghar et al., 2005).

Phylogenetic analysis based on available antennal transcriptomes and genomes from representative lepidopterans, including non-ditrysian basal orders showed that all SlitGST sequences fall into specific clades, with notable expansions in Noctuid family. This great number of orthologous GST groups in moths suggests recent radiation/expansion events in these species.

## Microsomal Antennal GSTs

In insects, the MAPEG clade contains in general few gene duplicates but seems involved in various metabolic pathways, including aging and pesticide detoxification. Indeed, a mutation in *D. melanogaster* MGST reduced lifespan, and in *Nilaparvata lugens*, GSTm2 expression level is induced by several pesticides (Toba and Aigaki, 2000; Zhou et al., 2013). If *S. littoralis* MAPEGs are very conserved, their expression pattern is remarkably

different. *SlitMGST1-1* and *SlitMGST1-2* were indeed broadly expressed in all tissues and development stages tested whereas *SlitMGST1-3* was restricted to the antenna suggesting that these three genes could have evolved different functions after a noctuid-specific duplication event. A previous study showed that *SlitMGST1-3* expression is induced after deltamethrin exposure (Lalouette et al., 2016), thus suggesting a possible specialization of *SlitMGST1-3* toward toxic molecules. In the beetle *P. striolata* two MGSTs are also preferentially expressed in the antennae of adult insects (Wu et al., 2016) meaning that members of this protein family could have conserved important function linked to olfaction.

## Sigma GSTs

We identified 4 sigma GSTs in *S. littoralis* antenna, which correspond to the second major expansion described in our phylogenetic analysis. Their expression pattern is similar in any tested tissue, they are highly expressed both in larval and
adult tissues with the exception of *SlitGSTs6* whose expression is almost restricted in the larval midgut and carcass. Moreover, *SlitGSTs6* clustered in a conserved group with other larval specific enzymes, like *PxylGSTs1* (You et al., 2015). These results are in agreement with the already described expression patterns of other sigma GSTs whose function is often associated with xenobiotic detoxification in insect larval midgut (Huang et al., 2011; Qin et al., 2014). In particular, their expression level is affected after pesticide exposure suggesting a role in response toward toxic compounds and subsequent generated oxidative stress.

Overall, our results suggest that sigma GSTs could be separate in two clades: one with a broad expression, including antenna, associated with basal conjugation functions and one larval-specific, which activity could be linked with oxidative stress response.

## Omega GSTs

Omega GSTs define a particular clade as their catalytic properties differ from other GSTs. Indeed, GSTo have a specific dehydroascorbate reductase and thiol transferase activities conferring an oxidative stress protection (Yamamoto et al., 2009a). This conserved function is associated with an apparent 1:1 orthology relationship of the three GSTo groups identified in our phylogenetic tree, with exception of *B. mori* (4 GSTo) and *P. xylostella* (5 GSTo). Moreover, expression analysis confirms a consistent pattern throughout any conditions. A GSTo in the silkworm is highly expressed in an insecticide-resistant strain and shows high affinity with organophosphate insecticides, indicating that it may contribute to insecticide resistance and oxidative stress responses, a potential conserved role across lepidopteran GSTs (Yamamoto et al., 2011).

## Zeta and Theta GSTs

The role of the zeta class GSTs had been first linked with phenylalanine and tyrosine catabolism, as maleylacetoacetate isomerases, suggesting a constitutive expression during every life stages (Board et al., 1997). However, *BmGSTz* had been associated with permethrin resistance with a predominant expression in fat body (Yamamoto et al., 2009b). This two opposite results illustrate the potential role of the two GSTz identified in Lepidoptera, which defined two highly supported clades with a strict 1:1 orthology between species. Thus, we can speculate about their function, with extremely conserved zeta1 GSTs as maleylacetoacetate isomerases (under strong purifying selection pressure) and zeta2 GSTs involved in insecticide resistance. We observed a ubiquitous expression of *SlitGSTz1* and *z2*, in agreement with their proposed functions.

As expected, we also identified a single theta GST in *S. littoralis*, widely distributed in various tested tissues. In *B. mori*, GSTt1 has been shown to possess a role in defense mechanisms against oxidative stress and in the metabolism of lipid peroxidation products (Yamamoto et al., 2005).

## Unclassified GSTs

According to our phylogeny, unclassified GSTs segregate in two paralogs groups, each composed of 1:1 orthologs from each

lepidopteran species. Functional information regarding those clades is scarce. *L. migratoria* *GSTu1* is expressed in Malpighian tubules and its downregulation using RNAi leads to a higher sensitivity to carbaryl and chlorpyrifos insecticides (Qin et al., 2013). In silkworm, *BmGSTu2* is induced in a resistant strain and is able to conjugate glutathione to the organophosphate insecticide diazinion (Yamamoto and Yamada, 2016). *SlitGSTu1* and *u2* are found to be expressed both at larval and adult stages, in all tissues tested, with a predominant expression of *SlitGSTu1*. It is likely that those genes could share similar functions than the one observed in other insects.

## Epsilon and Delta GSTs

Epsilon and delta clades are the most common GSTs in insects. They are widely recognized to have specific detoxification functions related to resistance to various insecticides (Enayati et al., 2005). In *S. littoralis*, epsilon clade accounts for almost half of the described sequences, with various expression patterns ranging from ubiquitous to antennal enriched genes. SlitGSTe clustered with lepidopteran GSTe functionally involved in insecticide conjugation and protection against oxidative stress in *B. mori*, *H. armigera*, *S. litura*, and *S. exigua* (Huang et al., 2011; Yamamoto et al., 2013; Liu et al., 2015b; Xu et al., 2015; You et al., 2015; Zhou et al., 2015; Wan et al., 2016; Zhang et al., 2016; Hirowatari et al., 2017; Labade et al., 2018). Moreover, three *P. xylostella* GSTe are also preferentially expressed in the antennae (He et al., 2017), suggesting a potential role in olfaction. *SlitGSTe9* and *SlitGSTe6* were predominantly expressed in the adult antenna, legs and larval midgut; in addition, *SlitGSTe15* was restricted to adult life stages and restricted to chemosensory tissues. This expression pattern suggests a possible role in olfaction and gustation aside from classical detoxification processes encountered in larval midgut.

Of particular interest is a clade containing single sequences from each lepidopteran species where *SlitGSTe16* appeared as a putative ortholog of *B. mori* *noperra-bo* (*BmorGSTe7*), a GST with cholesterol transporter activity involved in ecdysteroid biosynthesis (Enya et al., 2015). As *SlitGSTe16* is expressed in all tissues, a potential role in endocrine plasticity could be mediated by this enzyme, including in antennae, as ecdysteroids have been shown to modulate *S. littoralis* olfactory response (Bigot et al., 2012).

Delta GSTs sit in two conserved clades; the first one with *SlitGSTd1/d3* includes enzymes capable to metabolize pesticides in *B. mori* and *H. armigera* (Yamamoto et al., 2012; Labade et al., 2018) and exhibiting ubiquitous expression pattern. However, as far as expression information is available, the second one with *SlitGSTd2* is only composed of antennal-specific enzymes, such as *B. mori*, *M. sexta*, and *A. transitella* and *C. pomonella* GSTs (Rogers et al., 1999; Leal et al., 2009; Tan et al., 2014; Huang et al., 2017). MsexGST-msolf has been shown to degrade odorants *in vitro* and may have a role of ODE, especially toward aldehyde odorants, whereas CpomGSTd2 is active toward insecticides. SlitGSTd2 and the other GSTs from this clade may have evolved a function in odorant degradation and/or in protection of the ORNs toward toxic molecules in moth antennae. According to our analysis, *SlitGSTd2* is overexpressed in antenna and has the
highest expression level, compared to any other genes tested here. Such an expression pattern in consistent with an ODE function as olfaction associated proteins, like OBPs, use to be highly expressed in antenna and to have a rapid turnover in this tissue (Leal, 2013).

In addition, our bioinformatics analysis revealed that all the delta-2-like GST sequences from this clade possess a signal peptide signature, suggesting that they may be secreted proteins. More globally, our extensive analysis of signal peptide presence in GSTs from various insect orders revealed that this structural feature is more widely spread than suspected and not restricted to delta GSTs, as we found also a SP in an antennal GST sigma from a coleopteran species. Extracellular GSTs have been previously characterized in several Nematodes species (Sommer et al., 2001; Liebau et al., 2008; Espada et al., 2016b), and in the pine wood nematode *Bursaphelenchus xylophilus*; one of these enzyme metabolize various terpenoid compounds (Espada et al., 2016a), also known as common odorant molecules for insects. These antennal GSTs could thus be secreted in the sensillar lymph surrounding the sensory neurons, where they would directly interact with their relative substrate. Several mode of action are likely to occur in this aqueous environment: GSTs can conjugate their substrate with glutathione, as demonstrated in vertebrates where both GSTs and GSH are found in the olfactory mucus (Krishna et al., 1992; Debat et al., 2007); alternatively the binding properties of GSTs (as ligandins, Gonzalez et al., 2018) could act as a scavenger of odorants and harmful compounds.

We have revealed in this present work the occurrence of a high diversity of GST genes expressed in the olfactory organ of a pest moth. Phylogenetic analysis showed that these genes were distributed amongst the well-defined insect GSTs clades, in agreement with different cellular localization. The SlitGST structural diversity together with their different relative

spatial and developmental expression probably reflects their functional divergence and substrate specificities. Amongst this large repertoire, antennal GSTs could play a dual function in this tissue; first as detoxifying enzymes, where they could protect this delicate organ toward harmful compounds, but they could also play a role in the dynamic of olfactory signal: the conjugation of odorants (or their relative metabolites) could induce the termination of olfactory signaling. Moreover, the conjugation of such molecules may play a crucial role in odorant clearance, with the removal of any olfactory-active compounds in the sensillar lymph. Of particular interest is SlitGStd2, as it is antennal specific, probably secreted and likely to be involved in dual mechanisms in olfactory together with detoxification functions. Further studies using in vitro biochemical assays will reveal SliGStd2 function and substrate specificity, and will decipher if this enzyme is more related to toxic compounds, odorants or both. Overall future characterization, in particular biochemically but also physiologically, will allow to understand the precise function of all this enzymatic diversity in such a specialized organ, and to unravel their precise role in insect's biology as xenobiotic metabolizing enzymes and/or odorant degrading enzymes.

## AUTHOR CONTRIBUTIONS

ND, FB, and TC performed analysis. M-AP and DS provide material. ND, MM, and TC wrote the manuscript.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys. 2018.01283/full#supplementary-material

## REFERENCES

Anisimova, M., Gil, M., Dufayard, J. F., Dessimoz, C., and Gascuel, O. (2011). Survey of branch support methods demonstrates accuracy, power, and robustness of fast likelihood-based approximation schemes. *Syst. Biol.* 60, 685–699. doi: 10.1093/sysbio/syr041

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, 22–31. doi: 10.1016/j.ibmb.2011.10.003

Board, P. G., Baker, R. T., Chelvanayagam, G., and Jermiin, L. S. (1997). Zeta, a novel class of glutathione transferases in a range of species from plants to humans. *Biochem. J.* 328(Pt 3), 929–935.

Bock, K. W. (2010). Functions and transcriptional regulation of adult human hepatic UDP-glucuronosyl-transferases (UGTs): mechanisms responsible for interindividual variation of UGT levels. *Biochem. Pharmacol.* 80, 771–777. doi: 10.1016/j.bcp.2010.04.034

Bozzolan, F., Siaussat, D., Maria, A., Durand, N., Pottier, M. A., Chertemps, T., et al. (2014). Antennal uridine diphosphate (UDP)-glycosyltransferases in a pest insect: diversity and putative function in odorant and xenobiotics clearance. *Insect Mol. Biol.* 23, 539–549. doi: 10.1111/imb.12100

Cheng, T., Wu, J., Wu, Y., Chilukuri, R. V., Huang, L., Yamamoto, K., et al. (2017). Genomic adaptation to polyphagy and insecticides in a major East Asian noctuid pest. *Nat. Ecol. Evol.* 1, 1747–1756. doi: 10.1038/s41559-017-0314-4

Chertemps, T. (2017). “Molecular basis of pheromone detection in insects?,” in *Reference Module in Life Sciences* (Elsevier). Availabl online at: http://www.sciencedirect.com/science/article/pii/B9780128096338040383

Corcoran, J. A., Jordan, M. D., Thrimawithana, A. H., Crowhurst, R. N., and Newcomb, R. D. (2015). The peripheral olfactory repertoire of the lightbrown apple moth, *Epiphysas postvittana*. *PLoS ONE* 10:e0128596. doi: 10.1371/journal.pone.0128596

Debat, H., Eloit, C., Blon, F., Sarazin, B., Henry, C., Huet, J. C., et al. (2007). Identification of human olfactory cleft mucus proteins using proteomic analysis. *J. Proteome Res.* 6, 1985–1996. doi: 10.1021/pr0606575

Dermauw, W., and Van Leeuwen, T. (2014). The ABC gene family in arthropods: comparative genomics and role in insecticide transport and resistance. *Insect. Biochem. Mol. Biol.* 45, 89–110. doi: 10.1016/j.ibmb.2013.11.001

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:e15026. doi: 10.1371/journal.pone.0015026

Durand, N., Carot-Sans, G., Chertemps, T., Montagne, 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, 87–97. doi: 10.1111/j.1365-2583.2009.00939.x

Edgar, R. C. (2004). MUSCLE: multiple sequence alignment with high accuracy and high throughput. *Nucleic Acids Res.* 32, 1792–1797. doi: 10.1093/nar/gkh340
Elghar, G. E. A., Elbermawy, Z. A., Yousef, A. G., and Elhady, H. K. A. (2005). Monitoring and characterization of insecticide resistance in the cotton leafworm, *Spodoptera littoralis* (Boisd.) (*Lepidoptera: Noctuidae*). *J. Asia Pac. Entomol.* 8, 397-410. doi: 10.1016/S1226-8615(08)60263-0

Enayati, A. A., Ranson, H., and Hemingway, J. (2005). Insect glutathione transferases and insecticide resistance. *Insect Mol. Biol.* 14, 3–8. doi: 10.1111/j.1365-2583.2004.00529.x

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

Espada, M., Jones, J. T., and Mota, M. (2016a). Characterization of glutathione S-transferases from the pine wood nematode, *Bursaphelenchus xylophilus*. *Nematology* 18, 697. doi: 10.1163/15685411-00002985

Espada, M., Silva, A. C., van den Akker, S. E., Cock, P. J. A., Mota, M., and Jones, J. T. (2016b). Identification and characterization of parasitism genes from the pinewood nematode *Bursaphelenchus xylophilus* reveals a multilayered detoxification strategy. *Mol. Plant Pathol.* 17, 286–295. doi: 10.1111/mpp. 12280

Gonzalez, D., Fraichard, S., Grassein, P., Delarue, P., Senet, P., Nicolai, A., et al. (2018). Characterization of a *Drosophila* glutathione transferase involved in isothiocyanate detoxification. *Insect Biochem. Mol. Biol.* 95, 33–43. doi: 10.1016/j.ibmb.2018.03.004

Gouin, A., Bretaudieu, A., Nam, K., Gimenez, S., Aury, J. M., Duvic, B., et al. (2017). Two genomes of highly polyphagous lepidopteran pests (*Spodoptera frugiperda*, *Noctuidae*) with different host-plant ranges. *Sci. Rep.* 7, 11816. doi: 10.1038/s41598-017-10461-4

Gu, X. C., Zhang, Y. N., Kang, K., Dong, S. L., and Zhang, L. W. (2015). Antennal transcriptome analysis of odorant reception genes in the Red Turpentine Beetle (RTB), *Dendroctonus valens*. PLoS ONE 10:e0125159. doi: 10.1371/journal.pone.0125159

Guindon, S., Dufayard, J. F., Lefort, V., Anisimova, M., Hordijk, W., and Gascuel, O. (2010). New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. *Syst. Biol.* 59, 307–321. doi: 10.1093/sysbio/syq010

He, P., Zhang, Y. F., Hong, D. Y., Wang, J., Wang, X. L., Zuo, L. H., et al. (2017). A reference gene set for sex pheromone biosynthesis and degradation genes from the diamondback moth, *Plutella xylostella*, based on genome and transcriptome digital gene expression analyses. *BMC Genomics* 18:219. doi: 10.1186/s12864-017-3592-y

Hirowatari, A., Chen, Z., Mita, K., and Yamamoto, K. (2017). Enzymatic characterization of two epsilon-class glutathione S-transferases of *Spodoptera litura*. *Arch. Insect Biochem. Physiol.* 97:e21443. doi: 10.1002/arch.21443

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, 409–416. doi: 10.1093/aesa/sax041

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, 1033–1044. doi: 10.1016/j.jinsphys.2011.05.001

Jacquin-Joly, E., Legeai, F., Montagne, N., Monsempes, C., Francois, M. C., Poulain, J., et al. (2012). Candidate chemosensory Genes in Female Antennae of the Noctuid Moth *Spodoptera littoralis*. *Int. J. Biol. Sci.* 8, 1036–1050. doi: 10.7150/ijbs.4469

Jakoby, W. B., and Ziegler, D. M. (1990). The enzymes of detoxication. *J. Biol. Chem.* 265, 20715–20718.

Katoh, K., and Standley, D. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. *Mol. Biol. Evol.* 30, 772–780. doi: 10.1093/molbev/mst010

Kearse, M., Moir, R., Wilson, A., Stones-Havas, S., Cheung, M., Sturrock, S., et al. (2012). Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. *Bioinformatics* 28, 1647–1649. doi: 10.1093/bioinformatics/bts199

Krishna, N. S., Getchell, M. L., Tate, S. S., Margolis, F. L., and Getchell, T. V. (1992). Glutathione and gamma-glutamyl transpeptidase are differentially distributed in the olfactory mucosa of rats. *Cell Tissue Res.* 270, 475–484.

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., Wycké, 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, 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

Leal, W. S., Ishida, Y., Pelletier, J., Xu, W., Rayo, J., Xu, X., et al. (2009). Olfactory proteins mediating chemical communication in the navel orangeworm moth, *Amyelois transitella*. PLoS ONE 4:e7235. doi: 10.1371/journal.pone.0007235

Lefort, V., Longueville, J. E., and Gascuel, O. (2017). SMS: Smart Model Selection in PhyML. *Mol. Biol. Evol.* 34, 2422–2424. doi: 10.1093/molbev/msx149

Legeai, F., Malpel, S., Montagne, N., Monsempes, C., Cousserans, F., Merlin, C., et al. (2011). An expressed sequence Tag collection from the male antennae of the Noctuid moth *Spodoptera littoralis*: a resource for olfactory and pheromone detection research. *BMC Genomics* 12:86. doi: 10.1186/1471-2164-12-86

Li, X. C., Schuler, M. A., and Berenbaum, M. R. (2007). Molecular mechanisms of metabolic resistance to synthetic and natural xenobiotics. *Annu. Rev. Entomol.* 52, 231–253. doi: 10.1146/annurev.ento.51.110104.151104

Liebau, E., Hoppner, J., Muhlmeister, M., Burmeister, C., Luersen, K., Perbandt, M., et al. (2008). The secretory omega-class glutathione transferase OvgST3 from the human pathogenic parasite *Onchocerca volvulus*. *Fesb J.* 275, 3438–3453. doi: 10.1111/j.1742-4658.2008.06494.x

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: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 leaf folder, *Cnaphalocrocis medinalis* (*Lepidoptera: Pyralidae*): identification and expression profiles. *Arch. Insect Biochem. Physiol.* 90, 1–13. doi: 10.1002/arch.21240

Mamidala, P., Wijeratne, A. J., Wijeratne, S., Poland, T., Qazi, S. S., Doucet, D., et al. (2013). Identification of odor-processing genes in the emerald ash borer, *Agriplus planipennis*. PLoS ONE 8:e56555. doi: 10.1371/journal.pone.0056555

Petersen, T. N., Brunak, S., von Heijne, G., and Nielsen, H. (2011). SignalP 4.0: discriminating signal peptides from transmembrane regions. *Nat. Methods* 8, 785–786. doi: 10.1038/nmeth.1701

Poivet, E., Gallot, A., Montagne, N., Glaser, N., Legeai, F., and Jacquin-Joly, E. (2013). A comparison of the olfactory gene repertoires of adults and larvae in the noctuid moth *Spodoptera littoralis*. PLoS ONE 8:e60263. doi: 10.1371/journal.pone.0060263

Pottier, M. A., Bozzolan, F., Chertemps, T., Jacquin-Joly, E., Lalouette, L., Siaussat, D., et al. (2012). Cytochrome P450s and cytochrome P450 reductase in the olfactory organ of the cotton leafworm *Spodoptera littoralis*. *Insect Mol. Biol.* 21, 568–580. doi: 10.1111/j.1365-2583.2012.01160.x

Qin, G., Jia, M., Liu, T., Zhang, X., Guo, Y., Zhu, K. Y., et al. (2013). Characterization and functional analysis of four glutathione S-transferases from the migratory locust, *Locusta migratoria*. PLoS ONE 8:e58410. doi: 10.1371/journal.pone.0058410

Qin, G., Liu, T., Guo, Y., Zhang, X., Ma, E., and Zhang, J. (2014). Effects of chlorpyrifos on glutathione S-transferase in migratory locust, *Locusta migratoria*. *Pestic. Biochem. Physiol.* 109, 1–5. doi: 10.1016/j.pestbp.2013. 12.008

Robertson, H. M., Martos, R., Sears, C. R., Todres, E. Z., Walden, K. K., and Nardi, J. B. (1999). Diversity of odourant binding proteins revealed by an expressed sequence tag project on male *Manduca sexta* moth antennae. *Insect Mol. Biol.* 8, 501–518.

Rogers, M. E., Jani, M. K., and Vogt, R. G. (1999). An olfactory-specific glutathione-S-transferase in the sphinx moth *Manduca sexta*. *J. Exp. Biol.* 202(Pt 12), 1625–1637.

Salama, H. S., Dimetry, N. Z. and Salem, S. A. (1971). On the host preference and biology of the cotton leaf worm *Spodoptera littoralis* Bois. *Z Angew Entomol.* 67, 261–266. doi: 10.1111/j.1439-0418.1971.tb02122.x
Sheehan, D., Meade, G., Foley, V. M., and Dowd, C. A. (2001). Structure, function and evolution of glutathione transferases: implications for classification of non-mammalian members of an ancient enzyme superfamily. *Biochem. J.* 360(Pt 1), 1–16. doi: 10.1042/bj360001

Siaussat, D., Chertemps, T., and Maïbèche, M. (2014). "Detoxification, stress and immune responses in insect antenna: new insights from transcriptomics," in *Short Views on Insect Molecular Biology*, ed R. Chandrasekar (Tiruchirappalli: Bharathidasan University) 75–98.

Simon, P. (2003). Q-Gene: processing quantitative real-time RT-PCR data. *Bioinformatics* 19, 1439–1440. doi: 10.1093/bioinformatics/btg157

Sommer, A., Nimtz, M., Conradt, H. S., Brattig, N., Boettcher, K., Fischer, P., et al. (2001). Structural analysis and antibody response to the extracellular glutathione S-transferases from *Onchocerca volvulus*. *Infect. Immun.* 69, 7718–7728. doi: 10.1128/IAI.69.12.7718-7728.2001

Tan, X., Hu, X. M., Zhong, X. W., Chen, Q. M., Xia, Q. Y., and Zhao, P. (2014). Antenna-specific glutathione S-transferase in male silkworm *Bombyx mori*. *Int. J. Mol. Sci.* 15, 7429–7443. doi: 10.3390/ijms150 57429

Toba, G., and Aigaki, T. (2000). Disruption of the microsomal glutathione S-transferase-like gene reduces life span of *Drosophila melanogaster*. *Gene* 253, 179–187. doi: 10.1016/S0378-1119(00)00246-8

Tu, C. P., and Akgül, B. (2005). *Drosophila* glutathione S-transferases. *Meth. Enzymol.* 401, 204–226. doi: 10.1016/S0076-6879(05) 01013-X

Vogt, R. G., and Riddiford, L. M. (1981). Pheromone binding and inactivation by moth antennae. *Nature* 293, 161. doi: 10.1038/293161a0

Wan, H., Zhan, S., Xia, X., Xu, P., You, H., Jin, B. R., et al. (2016). Identification and functional characterization of an epsilon glutathione S-transferase from the beet armyworm (*Spodoptera exigua*). *Pestic. Biochem. Physiol.* 132, 81–88. doi: 10.1016/j.pestbp.2015.09.009

Wang, G. Y., Guo, Y.-Y., and Kong-Ming, W. (2004). Cloning of a cDNA fragment of an antenna-specific gene in *Helicoverpa armigera*. *Chin. J. Agric. Biotechnol.* 1, 37–43. doi: 10.1079/CJAB20048

Wu, Z. Z., Bin, S. Y., He, H. L., Wang, Z. B., Li, M., and Lin, J. T. (2016). Differential expression analysis of chemoreception genes in the striped flea beetle *Phyllotreta striolata* using a transcriptomic approach. *PLoS ONE* 11:e0153067. doi: 10.1371/journal.pone.0153067

Xu, Z. B., Zou, X. P., Zhang, N., Feng, Q. L., and Zheng, S. C. (2015). Detoxification of insecticides, allechemicals and heavy metals by glutathione S-transferase SIGSTE1 in the gut of *Spodoptera litura*. *Insect Sci.* 22, 503–511. doi: 10.1111/1744-7917.12142

Yamamoto, K., Aso, Y., and Yamada, N. (2013). Catalytic function of an epsilon-class glutathione S-transferase of the silkworm. *Insect Mol. Biol.* 22, 523–531. doi: 10.1111/imb.12041

Yamamoto, K., Nagaoka, S., Banno, Y., and Aso, Y. (2009a). Biochemical properties of an omega-class glutathione S-transferase of the silkworm, *Bombyx mori*. *Comp. Biochem. Physiol. C. Toxicol. Pharmacol.* 149, 461–467. doi: 10.1016/j.cbpc.2008.10.108

Yamamoto, K., Shigeoka, Y., Aso, Y., Banno, Y., Kimura, M., and Nakashima, T. (2009b). Molecular and biochemical characterization of a Zeta-class glutathione S-transferase of the silkworm. *Pestic. Biochem. Physiol.* 94, 30–35. doi: 10.1016/j.pestbp.2009.02.008

Yamamoto, K., Teshiba, S., Shigeoka, Y., Aso, Y., Banno, Y., Fujiki, T., et al. (2011). Characterization of an omega-class glutathione S-transferase in the stress response of the silkworm. *Insect Mol. Biol.* 20, 379–386. doi: 10.1111/j.1365-2583.2011.01073.x

Yamamoto, K., Usuda, K., Kakuta, Y., Kimura, M., Higashiura, A., Nakagawa, A., et al. (2012). Structural basis for catalytic activity of a silkworm Delta-class glutathione transferase. *Biochim. Biophys. Acta* 1820, 1469–1474. doi: 10.1016/j.bbagen.2012.04.022

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

Yamamoto, K., Zhang, P., Miake, F., Kashige, N., Aso, Y., Banno, Y., et al. (2005). Cloning, expression and characterization of theta-class glutathione S-transferase from the silkworm, *Bombyx mori*. *Comp. Biochem. Physiol. B Biochem. Mol. Biol.* 141, 340–346. doi: 10.1016/j.cbpc.2005.04.012

You, Y., Xie, M., Ren, N., Cheng, X., Li, J., Ma, X., et al. (2015). Characterization and expression profiling of glutathione S-transferases in the diamondback moth, *Plutella xylostella* (L.). *BMC Genomics* 16: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*. *Insect Biochem. Mol. Biol.* 53, 30–43. doi: 10.1016/j.ibmb.2014.07.003

Yuvaraj, J. K., Corcoran, J. A., Andersson, M. N., Newcomb, R. D., Anderbrant, O., and Lofstedt, C. (2017). Characterization of odorant receptors from a Non-ditrysian moth, *Eriocrania semipurpurella* sheds light on the origin of sex pheromone receptors in lepidoptera. *Mol. Biol. Evol.* 34, 2733–2746. doi: 10.1093/molbev/msx215

Zhang, N., Liu, J., Chen, S. N., Huang, L. H., Feng, Q. L., and Zheng, S. C. (2016). Expression profiles of glutathione S-transferase superfamily in *Spodoptera litura* tolerated to sublethal doses of chlorpyrifos. *Insect Sci.* 23, 675–687. doi: 10.1111/1744-7917.12202

Zhou, L., Fang, S. M., Huang, K., Yu, Q. Y., and Zhang, Z. (2015). Characterization of an epsilon-class glutathione S-transferase involved in tolerance in the silkworm larvae after long term exposure to insecticides. *Ecotoxicol. Environ. Saf.* 120, 20–26. doi: 10.1016/j.ecoenv.2015.05.027

Zhou, W. W., Liang, Q. M., Xu, Y., Gurr, G. M., Bao, Y. Y., Zhou, X. P., et al. (2013). Genomic Insights into the Glutathione S-transferase gene family of two rice planthoppers, *Nilaparvata lugens* (Stal) and *Sogatella furcifera* (Horvath) (Hemiptera: Delphacidae). *PLoS ONE* 8:e56604. doi: 10.1371/journal.pone.0056604

**Conflict of Interest Statement:** 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.

Copyright © 2018 Durand, Pottier, Siaussat, Bozzolan, Maïbèche and Chertemps.
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.