Skip to main content

ORIGINAL RESEARCH article

Front. Physiol., 26 September 2019
Sec. Invertebrate Physiology
This article is part of the Research Topic Invertebrate Neurobiology: Sensory Systems, Information Integration, Locomotor- and Behavioral Output View all 24 articles

Identification and Expression Analyses of Olfactory Gene Families in the Rice Grasshopper, Oxya chinensis, From Antennal Transcriptomes

  • Guang Zhou City Key Laboratory of Subtropical Fruit Tree Outbreak Control, Zhongkai University of Agriculture and Engineering, Guangzhou, China

The rice grasshopper Oxya chinensis is an important agricultural pest of rice and other gramineous plants. Chemosensory genes are crucial factors in direct interactions with odorants in the olfactory process. Here we identified genes encoding 18 odorant-binding proteins (OBPs), 13 chemosensory proteins (CSPs), 94 olfactory receptors (ORs), 12 ionotropic receptors (IRs), and two sensory neuron membrane proteins (SNMPs) from O. chinensis using an transcriptomic approach. Semi-quantitative RT-PCR assays revealed that six OBP-encoding genes (OchiOBP4, 5, 8, 9, 10, and 14), one CSP gene (OchiOBP10) and two IR genes (OchiIR28 and 29) were exclusively expressed in antennae, suggesting their roles in olfaction. Real-time quantitative PCR analyses revealed that genes expressed exclusively or predominantly in antennae also displayed significant differences in expression levels between males and females. Among the differentially expressed genes, 17 OR-encoding genes, one CSP- and one SNMP-gene showed female-biased expression, suggesting that they may be involved in some female-specific behaviors such as seeking oviposition site; whereas the three remaining OR-encoding genes showed male-biased expression, indicating their possible roles in sensing female sex pheromones. Our results laid a solid foundation for future studies to reveal olfactory mechanisms as well as designing strategies for controlling this rice pest.

Introduction

Olfaction plays an important role in regulating various physiological behaviors, such as foraging, mating, oviposition and avoiding predators (Leal, 2013). The process of olfactory perception is mediated by a series of peripheral olfactory proteins involving odorant-binding proteins (OBPs), chemosensory proteins (CSPs), sensory neuron membrane proteins (SNMPs), olfactory receptors (ORs), gustatory receptors (GRs), and ionotropic receptors (IRs) (Clyne et al., 1999, 2000; Vosshall et al., 1999; Galindo and Smith, 2001; Benton et al., 2007, 2009; Vieira and Rozas, 2011). OBPs and CSPs are suggested to be involve in binding and trafficking of hydrophobic odorant molecules across the sensillum lymph surrounding the olfactory sensory neurons (OSNs) on the sensilla of antennae and other chemosensory organs (Vogt et al., 1985; Sandler et al., 2000; Xu et al., 2005; Gomez-Diaz et al., 2013). ORs and GRs are the principle receptor proteins responsible for the detection of odorants and tastants, whereas IRs are involved in the detection of chemosensory, thermo- and hygro-sensory stimuli (Ai et al., 2010, 2013; Grosjean et al., 2011; Silbering et al., 2011; Kain et al., 2013; Su and Carlson, 2013; Koh et al., 2014; Chen et al., 2015; Stewart et al., 2015; Gorter et al., 2016; Hussain et al., 2016; Knecht et al., 2016, 2017; Ni et al., 2016). SNMPs play important roles in pheromone detection, which are expressed on the dendritic membranes of pheromone sensitive neurons (Rogers et al., 1997, 2001; Li et al., 2014; Jiang et al., 2016).

Gene families involved in olfactory perception have been extensively reported in many species in Lepidoptera, Coleoptera, Diptera, Hemiptera, and Hymenoptera during recent decades (Grosse-Wilde et al., 2011; Andersson et al., 2013, 2014; Brand et al., 2015; Paula et al., 2016). However, in the Orthoptera order, ORs have been identified from only three species: the genome of Locusta migratoria (142 ORs) (Wang et al., 2015), the antennal transcriptome of Schistocerca gregaria (119 ORs) (Pregitzer et al., 2017), and Oedaleus asiaticus (60 ORs) (Zhou et al., 2019). Members of the IR families have been identified in L. migratoria (32 IRs) (Wang et al., 2015) and O. asiaticus (6 IRs) (Zhou et al., 2019), the OBP families in S. gregaria (14 OBPs) (Jiang et al., 2017), Oedaleus infernalis (18 OBPs) (Zhang et al., 2018) and O. asiaticus (15 OBPs) (Zhou et al., 2019), the CSP families in L. migratoria (58 CSPs) and S. gregaria (42 CSPs) (Martin-Blazquez et al., 2017), and SNMP families in S. gregaria (2 SNMPs) (Jiang et al., 2016) and O. asiaticus (3 SNMPs) (Zhou et al., 2019). Consequently, additional Orthoptera species need to be investigated to reach a better understanding on olfactory receptive mechanisms for insect olfactory system.

Oxya chinensis is an oligophagous grasshopper pest and primarily feeds on graminaceous grasses. Electroantennogram (EAG) and behavioral bioassay have showed that O. chinensis adults have significantly higher olfactory sensitivity to geraniol compared to other host volatiles (Lu et al., 2008). Identification of target genes that are involved in olfactory perception may lead to environmentally-friendly approaches for controlling this pest (Venthur and Zhou, 2018). The objectives of this study were to generate antennal transcriptomes to identify major olfactory-related gene families (OBPs, CSPs, ORs, IRs, and SNMPs) from O. chinensis, examine expression profiles of the genes in various tissues to identify the antenna-specific genes, and compare their sex-specific expression patterns to identify female- and male-specific olfactory-related genes.

Materials and Methods

Insects, Tissue Collections, and RNA Extraction

The colony of O. chinensis used in this study was derived from insects originally collected from rice fields in the environs of Leizhou, Guangdong, China (110°05′E, 20°34′N). The colony was maintained in the laboratory on bristle grass in climate-controlled chamber under 30°C ± 1°C with 60% relative humidity and a photoperiod of 16 h light versus 8 h.

For transcriptomic analyses, 100 pairs of antennae were dissected separately from O. chinensis adult females and males. For semi-quantitative RT-PCR analysis, antennae, maxillary palps, foreleg tarsus, wings, and genitals were separately collected from adults (male: female = 1:1, n = 25 each). Two replicates were included for each tissue. For RT-qPCR analyses, 50 pairs of antennae were obtained from both females and males, separately along with other body parts including 20 pairs of maxillary palps, 20 pairs of foreleg tarsus, 20 pairs of wings, and 10 genitals from both males and females. Three replicates were included for each tissue sample. Tissues were homogenized to powder immediately in liquid nitrogen and stored at −80°C for further analyses.

Total RNA was extracted using Trizol reagent (Invitrogen, Carlsbad, CA, United States) and treated with RNase-free DNase I (Takara, Dalian, China) to remove potential genomic DNA contamination. RNA concentration, quality and quantity were analyzed on a NanoDrop ND-2000 Spectrophotometer (Nanodrop Technologies, United States) and a Qubit®2.0 Fluorometer (Invitrogen, Life Technologies, United States). RNA integrity was assessed with an Agilent 2100 Bioanalyzer (Agilent Technologies, United States).

cDNA Library Construction, Illumina Sequencing and de novo Assembly

Transcriptomes were generated through a commercial contract with Novogene Bioinformatics Technology, Co., Ltd. (Beijing, China). Libraries for sequencing were generated from 1.5 μg purified RNA per sample using a TruSeq RNA Sample Preparation Kit v2 (Illumina, San Diego, CA, United States) according to Illumina instructions and sequenced on an Illumina HiSeq 2500 platform (San Diego, CA, United States). Approximately 150 bp paired-end reads were generated.

Raw reads were firstly processed through in-house perl scripts. High-quality clean reads were obtained by removing reads containing adapter and unknown (poly-N) and low-quality reads. Clean reads were combined from both female and male antennae and were assembled into unigenes using Trinity (Version: r2013-11-10) with min_kmer_cov set to 2 under default settings (Grabherr et al., 2011). The clean reads from the O. chinensis female antennae and male antennae were deposited in the NCBI Sequence Read Archive (Female antennae: SAMN11484273; Male antennae: SAMN11484274).

Functional Annotation

Unigenes were searched against databases including NCBI non-redundant protein (nr), NCBI non-redundant nucleotide (nt), Swiss-Prot1, the Kyoto Encyclopedia of Genes and Genomes (KEGG)2, using BLASTx with a cut-of E-value of 10–5. Gene orthology (GO) and cluster of orthologous groups of proteins (COG)3 using Blast2GO program (Conesa et al., 2005; Gotz et al., 2008).

Gene Identification

To identify genes coding for OBPs, CSPs, ORs, IRs, and SNMPs, known protein sequences (OBPs, CSPs, ORs, IRs, and SNMPs) from other Orthopteran species were selected as queries to search the O. chinensis antennal transcriptomes. L. migratoria, S. gregaria, O. asiaticus, O. infernalis, and Ceracris kiangsu were used for OBPs; L. migratoria, S. gregaria, and O. asiaticus for CSPs; L. migratoria, S. gregaria, and O. asiaticus for ORs; L. migratoria and O. asiaticus for IRs; S. gregaria and O. asiaticus for SNMPs (Supplementary Table S1). tBLASTn was used to search and identify candidate chemosensory genes against the O. chinensis antennal transcriptomes, with a cut-of E-value of 10–5. Putative O. chinensis chemosensory genes were in turn used as queries to identify additional genes (tBLASTx and BLASTp). Repetitions were completed until no new candidates were identified. Candidate unigenes from the initial search were further manually examined using BLASTx against Genbank. Open reading frames (ORFs) were predicted with the ORF Finder in NCBI4. Amino acid sequences were aligned with MAFFT (version 7.308) (E-INS-I parameter set) (Katoh and Standley, 2013) and visualized with Geneious (version 9.1.3) (Kearse et al., 2012). Transmembrane domains, signal peptides, conserved cysteine locations in candidates were analyzed using the InterProScan tool plug-in in Geneious (Version: 9.1.3.) (Quevillon et al., 2005). Candidate genes coding for OBPs, CSPs, ORs, IRs, and SNMPs were listed in Supplementary Table S2, together with genetic characteristics, best matches in NCBI-nr database, protein domains and estimated expression levels.

Phylogenetic Analyses

Multiple amino acid sequence alignments were made with MAFFT (E-INS-I parameter) (Katoh and Standley, 2013). Phylogenetic trees were constructed using maximum likelihood analyses as implemented in FastTree2 [Jones-Taylor-Thornton (JTT) amino acid substitution model, 1,000 bootstrap replications] (Price et al., 2009, 2010). Dendrograms were created and colored in FigTree5. Sequences used for phylogenetic analysis are listed in Supplementary Table S3.

Analyses of Expression Levels and Differentially Expressed Genes

Gene expression levels were estimated using RSEM (Version: 1.2.15) (Li and Dewey, 2011) for each sample. The clean reads were mapped back to the assembled transcriptome. The read counts normalized using the Trimmed Mean of M-value normalization method (Robinson and Oshlack, 2010) for each mapped gene was used to calculate gene expression levels following the FPKM (fragments per kilobase per million read) method (Cock et al., 2010; Trapnell et al., 2010), and were used to identify differentially expressed genes between females and males using DEGseq edgeR package (Version: 3.4.2) (Robinson et al., 2010). Criteria for estimating significantly differentially expressed genes were set as q-value < 0.005 and the absolute value of log2 Fold Change > 1.

Semi-Quantitative RT-PCR and Real-Time Quantitative PCR

Semi-quantitative RT-PCR was carried out to compare expression levels of genes coding for OBPs, CSPs, and IRs in various tissues including antennae, maxillary palps, foreleg tarsus, wings and genitals. cDNA template was synthesized from total RNA using a PrimeScript RT reagent Kit (Takara, China). Two housekeeping genes, β-actin (Ochiβ-actin) and ribosomal protein 49 (Ochirp49) were used as controls. PCR reactions were conducted in a thermal cycler from Bio-Rad (CA, United States). PCR conditions were as follows: 95°C for 2 min, followed by 35 cycles of 95°C for 30 s, 56°C for 30 s, 72°C for 1 min, and a final extension for 10 min at 72°C. PCR products were separated on 1.5% agarose gels. Individual PCR reactions were repeated twice with independently isolated RNA samples.

Expression profiles of the antenna-predominant candidate genes for OBPs, CSPs, and IRs as well as all ORs and SNMPs were analyzed on a LightCycler 480 system (Roche Applied Science). RT-qPCR cycling parameters were set at 95°C for 5 min, followed by 45 cycles of 95°C for 10 s, and 60°C for 20 s. A melting curve analysis was then performed at 95°C for 20 s, 60°C for 30 s, and 95°C for 30 s in order to determine the specificity of primers. Negative controls without cDNA template (nuclease-free water) were included for each reaction. Three independent biological replications were performed with each biological replication measured in three technique replications. The results were analyzed using the LightCycler 480 Gene Scanning Software. Relative quantification was calculated following the 2–ΔΔCT method (Livak and Schmittgen, 2001), and normalized against the two reference genes (Ochiβ-actin and Ochirp49). Gene-specific primers were designed using Primer36 and listed in Supplementary Table S4.

Statistical Analysis

RT-qPCR data were analyzed and plotted using Prism 6.0 (GraphPad Software, CA, United States). The comparative analyses of each target gene among various tissues were determined using a one-way nested analysis of variance (ANOVA) followed by Duncan’s new multiple range test (α = 0.05) using SPSS 22.0 (SPSS Inc., Chicago, IL, United States). Values are presented as mean ± SE.

Results

Illumina Sequencing and de novo Assembly

A total of 53.4 and 44.8 million raw reads were obtained from male and female antennae cDNA libraries of O. chinensis, respectively. After removing adaptor sequences, low quality reads and contaminant sequences, 52.2 and 44.7 million of clean reads were generated. The reads were assembled into 120,803 unigenes with a mean length of 1,213 bp, and an N50 of 2,458 bp. The number of unigenes longer than 500 bp was 33,079, which accounted for 35.21% of all unigenes (Supplementary Table S5).

Functional Annotation

There were 26,923 (45.5%), 15,092 (25.5%), 20,546 (34.7%), 5,662 (9.57%), 10,911 (18.4%), and 20,542 (34.7%) unigenes that had homologous sequences in NCBI-nr, NCBI-nt, Swiss-Prot, GO, COG, and KEGG databases, respectively. Based on the homologous sequences, a total of 29,067 (49.1%) unigenes were annotated and the remaining unigenes were unmappable at present (Supplementary Figure S1A). GO analyses categorized annotated genes into three functional categories, including “biological process,” “cellular component,” and “molecular function” (Supplementary Figure S1B). In “biological process,” major subcategories included “cellular,” “metabolic” and “single-organism.” In “cellular component,” the subcategories “cell” and “cell part” and “membrane” were major subcategories. In “molecular function,” the subcategories “binding” and “catalytic activity” contained the largest numbers of unigenes.

Candidate Genes Coding for Odorant Binding Proteins

Eighteen OBP unigenes were identified from the O. chinensis antennal transcriptomes (Supplementary Table S2: Sheet 1). All identified unigenes except OchiOBP18 had a full-length ORF with sizes ranged from 128 to 266 amino acid residues. Except three OBP genes (OchiOBP12, 14, and 17), all predicted proteins had a predicted signal peptide. Sequence identities of predicted OBPs with those from other Orthopteran insects in the NCBI-nr database ranged from 45.9 to 88.6%, with an average of 73.8%. Consistent with previous Orthopteran OBP classification by Jiang et al. (2017), twelve OBPs (OchiOBP1, 3-8, 10, 12-14, 16) were classic OBPs, with the typical six conserved C-residues (Supplementary Figure S2). Two (OchiOBP9 and 17) were atypical OBPs, with extraordinary long stretches between conserved C1 and C2 (Supplementary Figure S3). Four OBPs were plus-C OBPs, which were further divided into plus-C OBP type-A (OchiOBP2, 11, and 15) (with the extra C-residue both in C4-C5 and C5-C6) (Supplementary Figure S4) and plus-C OBP type-B (OchiOBP18) (with the extra C-residue in front of C1 and behind C6) (Supplementary Figure S5). Phylogenetic analyses of OBPs from O. chinensis and other five locust species (O. infernalis, O. asiaticus, L. migratoria, S. gregaria, and C. kiangsu) showed that all the locust OBPs formed distinct clades based on the number and position of cysteine residues, and segregated into the classic OBP, atypical OBP, and plus-C OBP sub-families as well as other locust OBPs (Figure 1).

FIGURE 1
www.frontiersin.org

Figure 1. Phylogenetic tree of putative odorant-binding proteins (OBPs). Branch support (circles at the branch nodes) was estimated using an approximate likelihood ratio test based on the scale indicated on the top left. Bars indicate branch lengths in proportion to amino acid substitutions per site. Classic OBPs are in red; plus-C OBPs type-A in green; plus-C OBPs type-B in blue; and atypical OBPs in purple.

Candidate Genes Coding for Chemosensory Proteins

Thirteen unigenes encoding CSPs were identified from the O. chinensis antennal transcriptomes (Supplementary Table S2: Sheet 2). All but two (OchiCSP11 and 13) had full-length ORFs encoding 120–297 amino acid residues. All predicted OchiCSPs contained four highly conserved four-cysteine profiles (Supplementary Figure S6). Except OchiCSP1, all predicted proteins possessed a signal peptide. A phylogenetic tree was built using all predicted OchiCSPs together with those from other locust species. Most OchiCSPs had orthologs with CSPs from other locust species, and no O. chinensis-specific CSP lineage was evident (Figure 2).

FIGURE 2
www.frontiersin.org

Figure 2. Phylogenetic tree of putative chemosensory proteins (CSPs). Branch support (circles at the branch nodes) was estimated using an approximate likelihood ratio test based on the scale indicated at the top left. Bars indicate branch lengths in proportion to amino acid substitutions per site.

Candidate Genes Coding for Odorant Receptors

Ninety-four OR-encoding unigenes were identified, including one Orco (OchiOR1) and 93 conventional OR genes (OchiOR2-94), which were classified to the 7-transmembrane receptor superfamily (Supplementary Table S2: Sheet 1). Among these OR unigenes, 54 had full-length ORFs encoding proteins with 367 to 487 amino acid residues. The highly conserved co-receptor OchiOR1 shared 94.25% identity to a co-receptor from L. migratoria (ALD51504.1), while other OchiORs shared 33.56 to 91.13% identity with average 68.7% with the respective L. migratoria ORs. A phylogenetic tree was generated using our identified ORs along with a data set containing representative ORs from L. migratoria and S. gregaria (Figure 3). The locust ORs formed three distinct clades, with the O. chinensis OR co-receptor (OchiOR1) formed a clade with the protein Orco from two other related species. The vast majority of OchiORs were clustered with orthologs from other locust species. Only a few species-specific clades and sister pairs were observed. Three OchiORs (OchiOR27, 53 and 56) were segregated into unique clades, while OchiOR36/60 and OchiOR28/76 formed the sister pairs.

FIGURE 3
www.frontiersin.org

Figure 3. Phylogenetic tree of putative odorant receptors (ORs). The distance tree was rooted by the conservative ORco gene orthologs (red). The species- specific clades and sister pairs are labeled with red dots. Branch support (circles at the branch nodes) was estimated using an approximate likelihood ratio test based on the scale indicated at the top left. Bars indicate branch lengths in proportion to amino acid substitutions per site. The clade I is in orange; clade II in green; and clade III in blue.

Candidate Genes Coding for Ionotropic Receptors (iGluRs/IRs)

Twelve putative iGluR/IR unigenes were identified, which were predicted to encode ligand-binding domain (S1 and S2) with three transmembrane domains (M1, M2, and M3) or portions of domains (Supplementary Table S2: Sheet 2). Of these iGluR/IR unigenes, eight had complete ORFs encoding at least 535 amino acid residues. Distinct clades were observed in a phylogenetic tree generated with our identified sequences and paralogs from other species including Drosophila melanogaster and L. migratoria (Figure 4). All identified iGluRs/IRs from O. chinensis were assigned to two phylogenetic groups, including N-Methyl-D-aspartic acid (NMDA) iGluRs (OchiGluR8, 9, and 12) and antennal IRs (OchiIR25a, 8a, 76b, 22, 6, 17, 20, 28, and 29). No candidates were assigned to non-NMDA iGluRs and divergent IRs (Figure 4). A set of “antennal IR” conserved among other species were absent from O. chinensis and L. migratoria, and only three “IR co-receptor” orthologs (OchiIR25a, 8a, and 76b) were clustered with D. melanogaster orthologs.

FIGURE 4
www.frontiersin.org

Figure 4. Phylogenetic tree of putative ionotropic receptors (IRs). Branch support (circles at the branch nodes) was estimated using an approximate likelihood ratio test based on the scale indicated at the top left. Bars indicate branch lengths in proportion to amino acid substitutions per site.

Candidate Genes Coding for Sensory Neuron Membrane Proteins

Two SNMP unigenes were identified, which matched the CD36 family. The genes contained complete ORFs that encoded proteins with two transmembrane domains (Supplementary Table S2: Sheet 3). Phylogenetic analyses revealed that all SNMPs were classified into two distinct subgroups, SNMP1 and SNMP2 (Figure 5). In the SNMP1 clade, OchiSNMP1 and OchiSNMP2 were clustered into two subclades along with those from other locust species.

FIGURE 5
www.frontiersin.org

Figure 5. Phylogenetic tree of putative sensory neuron membrane proteins (SNMPs). Branch support (circles at the branch nodes) was estimated using an approximate likelihood ratio test based on the scale indicated at the top left. Bars indicate branch lengths in proportion to amino acid substitutions per site.

Transcript Abundance Based on FPKM

In this study, genes with FPKM values ≥1,000 were defined as highly expressed genes, those with values 300–1,000 were defined as moderately expressed genes, and those with values ≤300 were defined as weakly expressed genes. Based on these criteria, OchiOBP1-5 and OchiCSP1-4 were highly expressed in antennae (Supplementary Table S2: Sheets 1,2); OchiOR1 (Orco) was moderately expressed (Male FPKM: 318.52, Female FPKM: 391.73); other conventional ORs (OchiOR2-94), all iGluRs/IRs, and the two SNMPs were weakly expressed (RPKM ≤ 30) (Supplementary Table S2: Sheet 3). DEG analyses showed that OchiCSP10 and OchiOR2, 6, 8, 9, 11, 15, 33, 38, 39, 41, 43, 44, 48, 50, 53, 56, and 62 were female-predominant, whereas OchiCSP9, OchiOR34, 46 and 86, and OchiGluR8 and 12 were male-biased (Supplementary Table S2: Sheets 1–5).

Tissue- and Sex-Specific Expression

Based on semi-quantitative RT-PCR analyses, OchiOBP4, 5, 8, 9, 10, and 14 were almost exclusively expressed in antennae, while OchiOBP3 and 17 were abundant in antennae and maxillary palps (Figure 6A). The remaining OBP-encoding genes were abundant in multiple body parts. For CSP-encoding genes, OchiCSP10 was exclusively expressed in antennae, while other OchiCSPs were present multiple body parts (Figure 6B). For iGluR/IR-encoding genes, OchiIR28 and 29 were expressed predominately in antennae (Figure 6C).

FIGURE 6
www.frontiersin.org

Figure 6. Tissue-specific expression of OchiOBPs, OchiCSPs, and OchiIRs. (A) OBPs. (B) CSPs. (C) IRs. Genes exclusively or predominantly expressed in antennae are labeled in red dots. The β-actin (Ochiβ-actin) and ribosomal protein 49 (Ochirp49) were used as internal references to test the integrity of each cDNA template. An, antennae; M, maxillary palps; T, foreleg tarsus; W, wings; G, genitals.

All antenna-predominant genes were further analyzed using RT-qPCR (Figure 7). OchiOR2, 6, 8, 9, 11, 15, 33, 38, 39, 41, 43, 44, 48, 50, 53, 56, and 62 were significantly expressed at high levels in female antennae and OchiOR34, 46 and 86 were significantly expressed at higher levels in male antennae. The remaining OR-encoding unigenes were equally expressed in the antennae of both males and females (Figure 7A). For CSP-encoding genes, OchiCSP10 was expressed at significantly higher levels in female antennae. Additionally, the SNMP-encoding gene OchiSNMP1 were significantly upregulated in female antennae (Figure 7B). Other OBP- and IR-encoding genes were expressed largely equally in both males and females (Figure 7B).

FIGURE 7
www.frontiersin.org

Figure 7. Relative expression levels of the olfactory-related genes in the female antennae, male antennae, and other body parts from Oxya chinensis. (A) ORs. (B) SNMPs and the antenna-predominant candidates (OBPs, CSPs, and IRs). FA, female antennae; MA, male antennae; Bo, other body parts. The expression levels were estimated using the 2–ΔΔCT method. Bars represent standard error of three independent biological replicates with three technical duplicates for each replicate. Different small letters indicate statistically significant difference between tissues (p < 0.05, ANOVA, HSD).

Discussion

Identification and characterization of chemosensory genes are important steps toward understanding the evolution and primary functions of the insect olfactory system. In this study, we identified a set of candidate chemosensory genes from O. chinensis via analyzing antennal transcriptomes. Genetic and phylogenetic analyses of candidate chemosensory genes in O. chinensis was carried out to examine the similarities and differences of molecular components in chemosensory pathways. We further analyzed the expression profiles of chemosensory genes to identify olfaction-specific genes for future functional studies.

Odorant-binding proteins play vital roles in carrying odorants through the hemolymph to OR neurons, and transducing the resultant signals to downstream effectors in the olfactory system (Hallem et al., 2006). The number of OBPs identified in our antennal transcriptome was comparable to its related Locust species, namely 18 from O. infernalis, 16 from L. migratoria, 15 from O. asiaticus, and 14 from S. gregaria) (Ban et al., 2003; Xu et al., 2009; Yu et al., 2009; Zhang et al., 2015, 2018; Jiang et al., 2017). This may reflect physiological and evolutionary consistency between O. chinensis and other Locust species. OBPs from O. chinensis can be classified into four types: classic, atypical, plus-C (type-A and B), as well as Locust-specific OBPs (Jiang et al., 2017). Expression analysis revealed that five classic OBP (OchiOBP4, 5, 8, 10, and 14) and one atypical OBP (OchiOBP9) were expressed exclusively in antennae, suggesting that these genes are likely to have a role in antennal chemical-recognition processes. CSPs and OBPs can act as carriers for odorant molecules (Liu et al., 2012; Yi et al., 2013; Zhang T. et al., 2013; Zhang et al., 2014). However, in our RT-PCR analysis, only OchiCSP10 showed significantly antennae-biased expression, suggesting that this CSP could be involved in olfactory function.

Olfactory receptors, which connect binding proteins and OSNs to transduce olfactory signals, are the best known group of insect chemoreceptors. Like other locust species, O. chinensis possesses an expanded OR family (94) compared with Dipterans, Lepidopterans, and Hemipterans (Vieira and Rozas, 2011). A large number of OR-encoding genes may be associated with its ability for diverse host-odor perception to feed on a wider range of host plants. Great variation in the number of OR genes was found in different Locust species. Specifically, 94 OR genes were identified in O. chinensis, less than 120 from S. gregaria (Pregitzer et al., 2017), and 141 from L. migratoria (Wang et al., 2015), but significantly more than 60 from O. asiaticus (Zhou et al., 2019). The variation in gene numbers may reflect that O. chinensis inhabits a different ecological niche than other locust species. Generally, sexually dimorphic expression of ORs in antennae indicates possible pheromone receptors contributing sexual behaviors. Typically, Lepidopteran sex pheromones are released via females to attract males. Several moth sex pheromone for ORs have been functionally characterized, and most of them are expressed at higher levels in the male antennae (Krieger et al., 2004; Wanner et al., 2010; Zhang and Löfstedt, 2013). ORs expressed predominantly in female antennae are predicted to function in the detection of egg-laying-related odorant (Pelletier et al., 2010) or male released pheromones (or signals) (Anderson et al., 2009); ORs expressed evenly in the male and female antennae are predicted to function in general odorant perception (Yan et al., 2015). Therefore, we hypothesize that some or all of female-predominant OchiORs, including OchiOR2, 6, 8, 9, 11, 15, 33, 38, 39, 41, 43, 44, 48, 50, 53, 56, and 62), are likely involved in female specific behaviors such as finding plant hosts for oviposition. Male-predominant ORs (OchiOR34, 46, and 86) may be associated with detecting female sex pheromones. The remaining ORs with roughly equal expression levels in female and male antennae might be dedicated to general odorant detection.

Ionotropic receptors are conserved proteins and play key roles in the synaptic ligand-gated ion channels involved in olfaction, gustation, thermosensation and hygrosensation (Min et al., 2013; Zhang Y.V. et al., 2013; Koh et al., 2014; Chen et al., 2015; Stewart et al., 2015; Croset et al., 2016; Gorter et al., 2016; Hussain et al., 2016; Knecht et al., 2016, 2017; Ni et al., 2016; Prieto-Godino et al., 2016; Ganguly et al., 2017; Tauber et al., 2017). We identified nine antennal IR genes (OchiIR25a, 8a, 76b, 22, 6, 17, 20, 28, and 29) and three iGluRs (OchiGluR8, 9, and 12) from O. chinensis. In general, most antennal IRs were conserved across orders Diptera, Lepidoptera and Coleoptera, but only three IR co-receptor orthologs (OchiIR25a, 8a, and 76b) were found in O. chinensis, and there were no antennal IRs orthologs with D. melanogaster (Benton et al., 2009). In D. melanogaster, antennal IRs function as chemoreceptors (Benton et al., 2009) and are expressed in the peripheral olfactory neurons associated with the detection of amines, acids, general odors and sex pheromones (Ai et al., 2010; Silbering et al., 2011; Kain et al., 2013; Koh et al., 2014; Hussain et al., 2016). Notably, we found that two antennal IR genes (OchiIR28 and 29) were exclusively expressed in antennae, suggesting they potentially involved in odorant reception.

We also identified two SNMP unigenes. SNMPs are conserved throughout holometabolous insects and the SNMP1 subfamily are usually expressed in pheromone-sensitive OSNs, and mediate responses to lipid pheromones (Jin et al., 2008; Nichols and Vogt, 2008; Vogt et al., 2009; Gomez-Diaz et al., 2016). All SNMPs from locust species were clustered into the SNMP1 subclade. Both SNMPs were predominantly expressed in antennae in O. chinensis, supporting that these two SNMPs may perform functions in pheromone perception. Interestingly, OchiSNMP1 was predominantly expressed in female antennae (Figure 7B), suggesting that OchiSNMP1 may participate in sex pheromone reception.

Data Availability Statement

The datasets generated for this study can be found in the clean reads from the O. chinensis female antennae and male antennae were deposited in the NCBI Sequence Read Archive (Female antennae: SAMN11484273; Male antennae: SAMN11484274).

Author Contributions

YC and CK performed the experiments. ZW analyzed the data. ZW and JL wrote and revised the manuscript.

Funding

This work was supported by the Guang Zhou City Key Laboratory of Subtropical Fruit Tree Outbreak Control (Grant No. 201805010008) and the Serious Pest and Disease Mechanism and Sustainable Control Innovation Team of Guangdong Province (Grant No. 2017KCXTD018).

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.

Acknowledgments

We thank Dr. Mingshun Chen (Kansas State University, United States) for comments and editorial assistance on the manuscript.

Supplementary Material

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

FIGURE S1 | Transcriptome overview of O. chinensis antennae. (A) Unigenes annotated through different databases. (B) Gene ontology (GO) classification of the O. chinensis unigenes.

FIGURE S2 | The amino acid alignment of the predicted classic OBPs.

FIGURE S3 | The amino acid alignment of the predicted atypical OBPs.

FIGURE S4 | The amino acid alignment of the predicted plus-C OBPs type-A.

FIGURE S5 | The amino alignment of the predicted plus-C OBPs type-B.

FIGURE S6 | The amino alignment of the predicted CSPs.

TABLE S1 | List of OBPs, CSPs, ORs, IRs, and SNMPs in Orthopteran species.

TABLE S2 | A list of the identified OBPs, CSPs, ORs, IRs and SNMPs in O. chinensis, and detailed genetic characteristics, best matches in NCBI-nr database, protein domains and estimated expression levels found in each petal module.

TABLE S3 | The amino acid sequences used for phylogenetic analyses.

TABLE S4 | Gene-specific primers used for semi-quantitative RT-PCR and RT-qPCR.

TABLE S5 | Overview of the sequencing and assembly process.

Footnotes

  1. ^ http://www.ebi.ac.uk/uniprot
  2. ^ http://www.genome.jp/kegg/
  3. ^ http://clovr.org/docs/clusters-of-orthologous-groups-cogs/
  4. ^ https://www.ncbi.nlm.nih.gov/orffinder/
  5. ^ http://tree.bio.ed.ac.uk/software/figtree
  6. ^ http://primer3.ut.ee/

References

Ai, M., Blais, S., Park, J. Y., Min, S., Neubert, T. A., and Suh, G. S. (2013). Ionotropic glutamate receptors IR64a and IR8a form a functional odorant receptor complex in vivo in Drosophila. J. Neurosci. 33, 10741–10749. doi: 10.1523/JNEUROSCI.5419-12.2013

PubMed Abstract | CrossRef Full Text | Google Scholar

Ai, M., Min, S., Grosjean, Y., Leblanc, C., Bell, R., Benton, R., et al. (2010). Acid sensing by the Drosophila olfactory system. Nature 468, 691–695. doi: 10.1038/nature09537

PubMed Abstract | CrossRef Full Text | Google Scholar

Anderson, A. R., Wanner, K. W., Trowell, S. C., Warr, C. G., Jaquin-Joly, E., Zagatti, P., et al. (2009). Molecular basis of female-specific odorant responses in Bombyx mori. Insect Biochem. Mol. Biol. 39, 189–197. doi: 10.1016/j.ibmb.2008.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Andersson, M. N., Grosse-Wilde, E., Keeling, C. I., Bengtsson, J. M., Yuen, M. M., Li, M., et al. (2013). Antennal transcriptome analysis of the chemosensory gene families in the tree killing bark beetles, Ips typographus and Dendroctonus ponderosae (Coleoptera: Curculionidae: Scolytinae). BMC Genomics 14:198. doi: 10.1186/1471-2164-14-198

PubMed Abstract | CrossRef Full Text | Google Scholar

Andersson, M. N., Videvall, E., Walden, K. K., Harris, M. O., Robertson, H. M., and Löfstedt, C. (2014). Sex- and tissue-specific profiles of chemosensory gene expression in a herbivorous gall-inducing fly (Diptera: Cecidomyiidae). BMC Genomics 15:501. doi: 10.1186/1471-2164-15-501

PubMed Abstract | CrossRef Full Text | Google Scholar

Ban, L., Scaloni, A., D’Ambrosio, C., Zhang, L., Yahn, Y., and Pelosi, P. (2003). Biochemical characterization and bacterial expression of an odorant-binding protein from Locusta migratoria. Cell Mol. Life Sci. 60, 390–400. doi: 10.1007/s000180300032

PubMed Abstract | CrossRef Full Text | Google Scholar

Benton, R., Vannice, K. S., Gomez-Diaz, C., and Vosshall, L. B. (2009). Variant ionotropic glutamate receptors as chemosensory receptors in Drosophila. Cell 136, 149–162. doi: 10.1016/j.cell.2008.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Benton, R., Vannice, K. S., and Vosshall, L. B. (2007). An essential role for a CD36-related receptor in pheromone detection in Drosophila. Nature 450, 289–293. doi: 10.1038/nature06328

PubMed Abstract | CrossRef Full Text | Google Scholar

Brand, P., Ramírez, S. R., Leese, F., Quezada-Euan, J. J., Tollrian, R., and Eltz, T. (2015). Rapid evolution of chemosensory receptor genes in a pair of sibling species of orchid bees (Apidae: Euglossini). BMC Evol. Biol. 15:176. doi: 10.1186/s12862-015-0451-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, C., Buhl, E., Xu, M., Croset, V., Rees, J. S., Lilley, K. S., et al. (2015). Drosophila ionotropic receptor 25a mediates circadian clock resetting by temperature. Nature 527, 516–520. doi: 10.1038/nature16148

PubMed Abstract | CrossRef Full Text | Google Scholar

Clyne, P. J., Warr, C. G., and Carlson, J. R. (2000). Candidate taste receptors in Drosophila. Science 287, 1830–1834. doi: 10.1126/science.287.5459.1830

PubMed Abstract | CrossRef Full Text | Google Scholar

Clyne, P. J., Warr, C. G., Freeman, M. R., Lessing, D., Kim, J., and Carlson, J. R. (1999). A novel family of divergent seven-transmembrane proteins: candidate odorant receptors in Drosophila. Neuron 22, 327–338. doi: 10.1016/s0896-6273(00)81093-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Cock, P. J., Fields, C. J., Goto, N., Heuer, M. L., and Rice, P. M. (2010). The sanger FASTQ file format for sequences with quality scores, and the Solexa/Illumina FASTQ variants. Nucleic Acids Res. 38, 1767–1771. doi: 10.1093/nar/gkp1137

PubMed Abstract | CrossRef Full Text | Google Scholar

Conesa, A., Gotz, S., Garcia-Gomez, J. M., Terol, J., Talon, M., and Robles, M. (2005). Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21, 3674–3676. doi: 10.1093/bioinformatics/bti610

PubMed Abstract | CrossRef Full Text | Google Scholar

Croset, V., Schleyer, M., Arguello, J. R., Gerber, B., and Benton, R. (2016). A molecular and neuronal basis for amino acid sensing in the Drosophila larva. Sci. Rep. 6:34871. doi: 10.1038/srep34871

PubMed Abstract | CrossRef Full Text | Google Scholar

Galindo, K., and Smith, D. P. (2001). A large family of divergent Drosophila odorant-binding proteins expressed in gustatory and olfactory sensilla. Genetics 159, 1059–1072.

PubMed Abstract | Google Scholar

Ganguly, A., Pang, L., Duong, V. K., Lee, A., Schoniger, H., Varady, E., et al. (2017). A molecular and cellular context-dependent role for Ir76b in detection of amino acid taste. Cell Rep. 18, 737–750. doi: 10.1016/j.celrep.2016.12.071

PubMed Abstract | CrossRef Full Text | Google Scholar

Gomez-Diaz, C., Bargeton, B., Abuin, L., Bukar, N., Reina, J. H., Bartoi, T., et al. (2016). A CD36 ectodomain mediates insect pheromone detection via a putative tunnelling mechanism. Nat. Commun. 7:11866. doi: 10.1038/ncomms11866

PubMed Abstract | CrossRef Full Text | Google Scholar

Gomez-Diaz, C., Reina, J. H., Cambillau, C., and Benton, R. (2013). Ligands for pheromone-sensing neurons are not conformationally activated odorant binding proteins. PLoS Biol. 11:e1001546. doi: 10.1371/journal.pbio.1001546

PubMed Abstract | CrossRef Full Text | Google Scholar

Gorter, J. A., Jagadeesh, S., Gahr, C., Boonekamp, J. J., Levine, J. D., and Billeter, J. C. (2016). The nutritional and hedonic value of food modulate sexual receptivity in Drosophila melanogaster females. Sci. Rep. 6:19441. doi: 10.1038/srep19441

PubMed Abstract | CrossRef Full Text | Google Scholar

Gotz, S., Garcia-Gomez, 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, 3420–3435. doi: 10.1093/nar/gkn176

PubMed Abstract | CrossRef Full Text | Google Scholar

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, 644–652. doi: 10.1038/nbt.1883

PubMed Abstract | CrossRef Full Text | Google Scholar

Grosjean, Y., Rytz, R., Farine, J. P., Abuin, L., Cortot, J., Jefferis, G. S., et al. (2011). An olfactory receptor for food-derived odours promotes male courtship in Drosophila. Nature 478, 236–240. doi: 10.1038/nature10428

PubMed Abstract | CrossRef Full Text | Google Scholar

Grosse-Wilde, E., Kuebler, L. S., Bucks, S., Vogel, H., Wicher, D., and Hansson, B. S. (2011). Antennal transcriptome of Manduca sexta. Proc. Natl. Acad. Sci. U.S.A. 108, 7449–7454. doi: 10.1073/pnas.1017963108

PubMed Abstract | CrossRef Full Text | Google Scholar

Hallem, E. A., Dahanukar, A., and Carlson, J. R. (2006). Insect odor and taste receptors. Annu. Rev. Entomol. 51, 113–135. doi: 10.1146/annurev.ento.51.051705.113646

PubMed Abstract | CrossRef Full Text | Google Scholar

Hussain, A., Zhang, M., Ucpunar, H. K., Svensson, T., Quillery, E., Gompel, N., et al. (2016). Ionotropic chemosensory receptors mediate the taste and smell of polyamines. PLoS Biol. 14:e1002454. doi: 10.1371/journal.pbio.1002454

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, X., Krieger, J., Breer, H., and Pregitzer, P. (2017). Distinct subfamilies of odorant binding proteins in locust (Orthoptera, Acrididae): molecular evolution, structural variation, and sensilla-specific expression. Front. Physiol. 8:734. doi: 10.3389/fphys.2017.00734

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, X., Pregitzer, P., Grosse-Wilde, E., Breer, H., and Krieger, J. (2016). Identification and characterization of two “Sensory Neuron Membrane Proteins” (SNMPs) of the desert locust, Schistocerca gregaria (Orthoptera: Acrididae). J. Insect Sci. 16:33. doi: 10.1093/jisesa/iew015

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, X., Ha, T. S., and Smith, D. P. (2008). SNMP is a signaling component required for pheromone sensitivity in Drosophila. Proc. Natl. Acad. Sci. U.S.A. 105, 10996–11001. doi: 10.1073/pnas.0803309105

PubMed Abstract | CrossRef Full Text | Google Scholar

Kain, P., Boyle, S. M., Tharadra, S. K., Guda, T., Christine, P., Dahanukar, A., et al. (2013). Odour receptors and neurons for DEET and new insect repellents. Nature 502, 507–512. doi: 10.1038/nature12594

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Knecht, Z. A., Silbering, A. F., Cruz, J., Yang, L., Croset, V., Benton, R., et al. (2017). Ionotropic receptor-dependent moist and dry cells control hygrosensation in Drosophila. eLife 6:e26654. doi: 10.7554/eLife.26654

PubMed Abstract | CrossRef Full Text | Google Scholar

Knecht, Z. A., Silbering, A. F., Ni, L., Klein, M., Budelli, G., Bell, R., et al. (2016). Distinct combinations of variant ionotropic glutamate receptors mediate thermosensation and hygrosensation in Drosophila. eLife 5:e17879. doi: 10.7554/eLife.17879

PubMed Abstract | CrossRef Full Text | Google Scholar

Koh, T. W., He, Z., Gorur-Shandilya, S., Menuz, K., Larter, N. K., Stewart, S., et al. (2014). The Drosophila IR20a clade of ionotropic receptors are candidate taste and pheromone receptors. Neuron 83, 850–865. doi: 10.1016/j.neuron.2014.07.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Krieger, J., Grosse-Wilde, E., Gohl, T., Dewer, Y. M., Raming, K., and Breer, H. (2004). Genes encoding candidate pheromone receptors in a moth (Heliothis virescens). Proc. Natl. Acad. Sci. U.S.A. 101, 11845–11850. doi: 10.1073/pnas.0403052101

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, B., and Dewey, C. N. (2011). RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics 12:323. doi: 10.1186/1471-2105-12-323

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Z., Ni, J. D., Huang, J., and Montell, C. (2014). Requirement for Drosophila SNMP1 for rapid activation and termination of pheromone-induced activity. PLoS Genet. 10:e1004600. doi: 10.1371/journal.pgen.1004600

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, R., He, X., Lehane, S., Lehane, M., Hertz-Fowler, C., Berriman, M., et al. (2012). Expression of chemosensory proteins in the tsetse fly Glossina morsitans morsitans is related to female host-seeking behaviour. Insect Mol. Biol. 21, 41–48. doi: 10.1111/j.1365-2583.2011.01114.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2-DELTADELTACT method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, F. P., Zhao, D. X., and Wang, A. P. (2008). Electroantennogram and behavioral responses of Oxya chinensis (Orthoptera: Acrididae) to plant volatiles. Chin. J. Trop. Crops 29, 225–230.

Google Scholar

Martin-Blazquez, R., Chen, B., Kang, L., and Bakkali, M. (2017). Evolution, expression and association of the chemosensory protein genes with the outbreak phase of the two main pest locusts. Sci. Rep. 7:6653. doi: 10.1038/s41598-017-07068-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Min, S., Ai, M., Shin, S. A., and Suh, G. S. (2013). Dedicated olfactory neurons mediating attraction behavior to ammonia and amines in Drosophila. Proc. Natl. Acad. Sci. U.S.A. 110, E1321–E1329. doi: 10.1073/pnas.1215680110

PubMed Abstract | CrossRef Full Text | Google Scholar

Ni, L., Klein, M., Svec, K. V., Budelli, G., Chang, E. C., Ferrer, A. J., et al. (2016). The ionotropic receptors IR21a and IR25a mediate cool sensing in Drosophila. eLife 5:e13254. doi: 10.7554/eLife.13254

PubMed Abstract | CrossRef Full Text | Google Scholar

Nichols, Z., and Vogt, R. G. (2008). The SNMP/CD36 gene family in Diptera, Hymenoptera and Coleoptera: Drosophila melanogaster, D-pseudoobscura, Anopheles gambiae, Aedes aegypti, Apis mellifera, and Tribolium castaneum. Insect Biochem. Mol. Biol. 38, 398–415. doi: 10.1016/j.ibmb.2007.11.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Paula, D. P., Togawa, R. C., Costa, M. M., Grynberg, P., Martins, N. F., and Andow, D. A. (2016). Identification and expression profile of odorant-binding proteins in Halyomorpha halys (Hemiptera: Pentatomidae). Insect Mol. Biol. 25, 580–594. doi: 10.1111/imb.12243

PubMed Abstract | CrossRef Full Text | Google Scholar

Pelletier, J., Hughes, D. T., Luetje, C. W., and Leal, W. S. (2010). An odorant receptor from the southern house mosquito Culex pipiens quinquefasciatus sensitive to oviposition attractants. PLoS One 5:e10090. doi: 10.1371/journal.pone.0010090

PubMed Abstract | CrossRef Full Text | Google Scholar

Pregitzer, P., Jiang, X., Grosse-Wilde, E., Breer, H., Krieger, J., and Fleischer, J. (2017). In search for pheromone receptors: certain members of the odorant receptor family in the desert locust Schistocerca gregaria (Orthoptera: Acrididae) Are Co-expressed with SNMP1. Int. J. Biol. Sci. 13, 911–922. doi: 10.7150/ijbs.18402

PubMed Abstract | CrossRef Full Text | Google Scholar

Price, M. N., Dehal, P. S., and Arkin, A. P. (2009). FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol. Biol. Evol. 26, 1641–1650. doi: 10.1093/molbev/msp077

PubMed Abstract | CrossRef Full Text | Google Scholar

Price, M. N., Dehal, P. S., and Arkin, A. P. (2010). FastTree 2–approximately maximum-likelihood trees for large alignments. PLoS One 5:e9490. doi: 10.1371/journal.pone.0009490

PubMed Abstract | CrossRef Full Text | Google Scholar

Prieto-Godino, L. L., Rytz, R., Bargeton, B., Abuin, L., Arguello, J. R., Peraro, M. D., et al. (2016). Olfactory receptor pseudo-pseudogenes. Nature 539, 93–97. doi: 10.1038/nature19824

PubMed Abstract | CrossRef Full Text | Google Scholar

Quevillon, E., Silventoinen, V., Pillai, S., Harte, N., Mulder, N., Apweiler, R., et al. (2005). InterProScan: protein domains identifier. Nucleic Acids Res. 33, W116–W120.

PubMed Abstract | Google Scholar

Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140. doi: 10.1093/bioinformatics/btp616

PubMed Abstract | CrossRef Full Text | Google Scholar

Robinson, M. D., and Oshlack, A. (2010). A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 11:R25. doi: 10.1186/gb-2010-11-3-r25

PubMed Abstract | CrossRef Full Text | Google Scholar

Rogers, M. E., Krieger, J., and Vogt, R. G. (2001). Antennal SNMPs (Sensory Neuron Membrane Proteins) of lepidoptera define a unique family of invertebrate CD36-like proteins. J. Neurobiol. 49, 47–61. doi: 10.1002/neu.1065

PubMed Abstract | CrossRef Full Text | Google Scholar

Rogers, M. E., Sun, M., Lerner, M. R., and Vogt, R. G. (1997). Snmp-1, a novel membrane protein of olfactory neurons of the silk moth Antheraea polyphemus with homology to the CD36 family of membrane proteins. J. Biol. Chem. 272, 14792–14799. doi: 10.1074/jbc.272.23.14792

PubMed Abstract | CrossRef Full Text | Google Scholar

Sandler, B. H., Nikonova, L., Leal, W. S., and Clardy, J. (2000). Sexual attraction in the silkworm moth: structure of the pheromone-binding-protein-bombykol complex. Chem. Biol. 7, 143–151. doi: 10.1016/S1074-5521(00)00078-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Silbering, A. F., Rytz, R., Grosjean, Y., Abuin, L., Ramdya, P., Jefferis, G. S., et al. (2011). Complementary function and integrated wiring of the evolutionarily distinct drosophila olfactory subsystems. J. Neurosci. 31, 13357–13375. doi: 10.1523/JNEUROSCI.2360-11.2011

PubMed Abstract | CrossRef Full Text | Google Scholar

Stewart, S., Koh, T. W., Ghosh, A. C., and Carlson, J. R. (2015). Candidate ionotropic taste receptors in the Drosophila larva. Proc. Natl. Acad. Sci. U.S.A. 112, 4195–4201. doi: 10.1073/pnas.1503292112

PubMed Abstract | CrossRef Full Text | Google Scholar

Su, C., and Carlson, J. R. (2013). Neuroscience. Circuit logic of avoidance and attraction. Science 340, 1295–1297. doi: 10.1126/science.1240139

PubMed Abstract | CrossRef Full Text | Google Scholar

Tauber, J. M., Brown, E. B., Li, Y., Yurgel, M. E., Masek, P., and Keene, A. C. (2017). A subset of sweet-sensing neurons identified by IR56d are necessary and sufficient for fatty acid taste. PLoS Genet. 13:e1007059. doi: 10.1371/journal.pgen.1007059

PubMed Abstract | CrossRef Full Text | Google Scholar

Trapnell, C., Williams, B. A., Pertea, G., Mortazavi, A., Kwan, G., van Baren, M. J., et al. (2010). Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 28, 511–515. doi: 10.1038/nbt.1621

PubMed Abstract | CrossRef Full Text | Google Scholar

Venthur, H., and Zhou, J. J. (2018). Odorant receptors and odorant-binding proteins as insect pest control targets: a comparative analysis. Front. Physiol. 9:1163. doi: 10.3389/fphys.2018.01163

CrossRef Full Text | Google Scholar

Vieira, F. G., and Rozas, J. (2011). Comparative genomics of the odorant-binding and chemosensory protein gene families across the arthropoda: origin and evolutionary history of the chemosensory system. Genome Biol. Evol. 3, 476–490. doi: 10.1093/gbe/evr033

PubMed Abstract | CrossRef Full Text | Google Scholar

Vogt, R. G., Miller, N. E., Litvack, R., Fandino, R. A., Sparks, J., Staples, J., et al. (2009). The insect SNMP gene family. Insect Biochem. Mol. Biol. 39, 448–456. doi: 10.1016/j.ibmb.2009.03.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Vogt, R. G., Riddiford, L. M., and Prestwich, G. D. (1985). Kinetic properties of a sex pheromone-degrading enzyme: the sensillar esterase of Antheraea polyphemus. Proc. Natl. Acad. Sci. U.S.A. 82, 8827–8831. doi: 10.1073/pnas.82.24.8827

PubMed Abstract | CrossRef Full Text | Google Scholar

Vosshall, L. B., Amrein, H., Morozov, P. S., Rzhetsky, A., and Axel, R. (1999). A spatial map of olfactory receptor expression in the Drosophila antenna. Cell 96, 725–736. doi: 10.1016/s0092-8674(00)80582-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Z., Yang, P., Chen, D., Jiang, F., Li, Y., Wang, X., et al. (2015). Identification and functional analysis of olfactory receptor family reveal unusual characteristics of the olfactory system in the migratory locust. Cell Mol. Life Sci. 72, 4429–4443. doi: 10.1007/s00018-015-2009-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Wanner, K. W., Nichols, A. S., Allen, J. E., Bunger, P. L., Garczynski, S. F. Jr., Linn, C. E., et al. (2010). Sex pheromone receptor specificity in the european corn borer moth, Ostrinia nubilalis. PLoS One 5:e8685. doi: 10.1371/journal.pone.0008685

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, P., Atkinson, R., Jones, D. N., and Smith, D. P. (2005). Drosophila OBP LUSH is required for activity of pheromone-sensitive neurons. Neuron 45, 193–200. doi: 10.1016/j.neuron.2004.12.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, Y. L., He, P., Zhang, L., Fang, S. Q., Dong, S. L., Zhang, Y. J., et al. (2009). Large-scale identification of odorant-binding proteins and chemosensory proteins from expressed sequence tags in insects. BMC Genomics 10:632. doi: 10.1186/1471-2164-10-632

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, S. W., Zhang, J., Liu, Y., Li, G. Q., and Wang, G. R. (2015). An olfactory receptor from Apolygus lucorum (Meyer-Dur) mainly tuned to volatiles from flowering host plants. J. Insect Physiol. 79, 36–41. doi: 10.1016/j.jinsphys.2015.06.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Yi, X., Zhao, H., Dong, X., Wang, P., Hu, M., and Zhong, G. (2013). BdorCSP2 is important for antifeed and oviposition-deterring activities induced by Rhodojaponin-III against Bactrocera dorsalis. PLoS One 8:e77295. doi: 10.1371/journal.pone.0077295

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, F., Zhang, S., Zhang, L., and Pelosi, P. (2009). Intriguing similarities between two novel odorant-binding proteins of locusts. Biochem. Biophys. Res. Commun. 385, 369–374. doi: 10.1016/j.bbrc.2009.05.074

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, D. D., and Löfstedt, C. (2013). Functional evolution of a multigene family: orthologous and paralogous pheromone receptor genes in the turnip moth, Agrotis segetum. PLoS One 8:e77345. doi: 10.1371/journal.pone.0077345

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, S., Pang, B., and Zhang, L. (2015). Novel odorant-binding proteins and their expression patterns in grasshopper, Oedaleus asiaticus. Biochem. Biophys. Res. Commun. 460, 274–280. doi: 10.1016/j.bbrc.2015.03.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, T., Wang, W., Zhang, Z., Zhang, Y., and Guo, Y. (2013). Functional characteristics of a novel chemosensory protein in the cotton bollworm Helicoverpa armigera (Hubner). J. Integr. Agric. 12, 853–861. doi: 10.1016/s2095-3119(13)60304-4

CrossRef Full Text | Google Scholar

Zhang, Y., Tan, Y., Zhou, X. R., and Pang, B. P. (2018). A whole-body transcriptome analysis and expression profiling of odorant binding protein genes in Oedaleus infernalis. Comp. Biochem. Physiol. Part D Genomics Proteomics 28, 134–141. doi: 10.1016/j.cbd.2018.08.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y. N., Ye, Z. F., Yang, K., and Dong, S. L. (2014). Antenna-predominant and male-biased CSP19 of Sesamia inferens is able to bind the female sex pheromones and host plant volatiles. Gene 536, 279–286. doi: 10.1016/j.gene.2013.12.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y. V., Ni, J., and Montell, C. (2013). The molecular basis for attractive salt-taste coding in Drosophila. Science 340, 1334–1338. doi: 10.1126/science.1234133

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Y. T., Li, L., Zhou, X. R., Tan, Y., and Pang, B. P. (2019). Identification and expression profiling of candidate chemosensory membrane proteins in the band-winged grasshopper, Oedaleus asiaticus. Comp. Biochem. Physiol. Part D Genomics Proteomics 30, 33–44. doi: 10.1016/j.cbd.2019.02.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: antennal transcriptome, olfactory gene, identification, expression analysis, Oxya chinensis

Citation: Cui Y, Kang C, Wu Z and Lin J (2019) Identification and Expression Analyses of Olfactory Gene Families in the Rice Grasshopper, Oxya chinensis, From Antennal Transcriptomes. Front. Physiol. 10:1223. doi: 10.3389/fphys.2019.01223

Received: 27 July 2019; Accepted: 09 September 2019;
Published: 26 September 2019.

Edited by:

Sylvia Anton, Institut National de la Recherche Agronomique (INRA), France

Reviewed by:

Pablo Pregitzer, University of Hohenheim, Germany
Dan-Dan Zhang, Lund University, Sweden

Copyright © 2019 Cui, Kang, Wu and Lin. 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: Zhongzhen Wu, emhvbmd6aGVuX3d1QDE2My5jb20=; Jintian Lin, bGluanRpYW5AMTYzLmNvbQ==

These authors have contributed equally to this work

Disclaimer: 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.