Skip to main content

Advertisement

Sex- and tissue-specific transcriptome analyses and expression profiling of olfactory-related genes in Ceracris nigricornis Walker (Orthoptera: Acrididae)

Article metrics

Abstract

Background

The sophisticated insect olfactory system plays an important role in recognizing external odors and enabling insects to adapt to environment. Foraging, host seeking, mating, ovipositing and other forms of chemical communication are based on olfaction, which requires the participation of multiple olfactory genes. The exclusive evolutionary trend of the olfactory system in Orthoptera insects is an excellent model for studying olfactory evolution, but limited olfaction research is available for these species. The olfactory-related genes of Ceracris nigricornis Walker (Orthoptera: Acrididae), a severe pest of bamboos, have not yet been reported.

Results

We sequenced and analyzed the transcriptomes from different tissues of C. nigricornis and obtained 223.76 Gb clean data that were assembled into 43,603 unigenes with an N50 length of 2235 bp. Among the transcripts, 66.79% of unigenes were annotated. Based on annotation and tBLASTn results, 112 candidate olfactory-related genes were identified for the first time, including 20 odorant-binding proteins (OBPs), 10 chemosensory-binding proteins (CSPs), 71 odorant receptors (ORs), eight ionotropic receptors (IRs) and three sensory neuron membrane proteins (SNMPs). The fragments per kilobase per million mapped fragments (FPKM) values showed that most olfactory-related differentially expressed genes (DEGs) were enriched in the antennae, and these results were confirmed by detecting the expression of olfactory-related genes with quantitative real-time PCR (qRT-PCR). Among these antennae-enriched genes, some were sex-biased, indicating their different roles in the olfactory system of C. nigricornis.

Conclusions

This study provides the first comprehensive list and expression profiles of olfactory-related genes in C. nigricornis and a foundation for functional studies of these olfactory-related genes at the molecular level.

Background

Ceracris nigricornis Walker (Orthoptera: Acrididae) is a severe grasshopper pest of bamboos such as Phyllostachys heterocycla, Phyllostachys viridis and Phyllostachys glauca. C. nigricornis can also harm rice, corn, sorghum and other crops and can cause serious economic losses. Typically, the application of a substantial quantity of chemical insecticides, especially wide-spectrum insecticides, is the main method for controlling this pest. However, long-term application of pesticides may lead to pesticide resistance, pesticide residues, environmental pollution, and a decrease in the natural enemies of C. nigricornis [1,2,3,4]. In recent years, the use of eco-friendly nonhost plant volatiles to control phytophagous insects has increased; for example, plant volatiles from Trifolium repens L., Castanea mollissima Blume, Citrus reticulata Blanco, Kigelia africana (Lam.) and Myrica rubra (Lour.) have been used to interfere with the orientation and selection of plant volatiles of tea leaves in the olfactory system of Empoasca vitis, which reduces the level of E. vitis [5,6,7]. The ability of these nonhost plant volatiles to control the level of insects depends largely on the highly sensitive insect olfactory system [8]. Therefore, the elucidation of the molecular basis of the insect olfactory system is of great importance for new bio-pesticide development and pest control.

Olfaction is the primary sensory modality in insects and plays an important role in finding mating partners, food, oviposition sites and suitable habitats [9,10,11]. The insect olfactory system involves several different proteins, including binding proteins (odorant-binding proteins, OBPs; and chemosensory-binding proteins, CSPs), chemosensory membrane proteins (odorant receptors, ORs; ionotropic receptors, IRs; gustatory receptors, GRs; and sensory neuron membrane proteins, SNMPs), and odorant-degrading esterases (ODEs) [12, 13]. OBPs and CSPs are highly concentrated in the lymph of chemosensilla and are regarded as carriers of pheromones and odorants in insect chemoreception [13,14,15]. OBPs are small globular, water-soluble proteins that generally contain six highly conversed cysteine residues paired into three interlocking disulfide bridges [16, 17]. OBPs can bind and transport external odorant molecules to the olfactory receptors in the olfactory neuronal membrane, which is often considered the first step in olfactory recognition [18, 19]. CSPs are also small soluble proteins, also known as olfactory system of Drosophila melanogaster (OS-D)-like proteins or sensory appendage proteins, which contain only four conserved cysteine residues but have more conserved nucleotide sequences than OBPs across insect species [20,21,22]. CSPs are expressed in various chemosensory organs and have many functions. CSPs are also present in nonchemosensory organs and play a role in the transmission of pheromones, the solubility of nutrients, and the development of insecticide resistance [23,24,25].

The recognition and transmission of olfaction begins with the interaction between odorant molecules and ORs on the dendrites of olfactory receptor neurons (ORNs). Insect ORs were first identified in the D. melanogaster genome; ORs contain seven transmembrane domains (TMDs) and a membrane topology with an intracellular N-terminus and an extracellular C-terminus, and the membrane topology of insect ORs are reversed compared to that of vertebrate ORs [26]. ORs are nonselective cationic channels with high selectivity and specificity for odorant molecules. ORs can convert chemical signals of odorant molecules into electrical signals and play a role as a transit station in insect olfactory reactions. IR is a newly discovered gene family that was first studied in the olfactory system of D. melanogaster [27]. IRs evolved from the ionotropic glutamate receptor superfamily (iGluRs) and contains iGluRs conserved structural regions: three TMDs, a bipartite ligand-binding domain with two lobes and one ion channel pore [28]. IRs are expressed in coeloconic olfactory sensory neurons (OSNs) without ORs or coreceptors (ORcos) and mainly function in detecting acids, amines and other chemicals that cannot be recognized by ORs [29]. SNMPs are double transmembrane proteins that have a transmembrane domain at the C- and N-termini of the chain. SNMPs belong to the CD36 receptor family and are divided into two subfamilies, SNMP1 and SNMP2 [30]. The SNMP1 subfamily is coexpressed with pheromone receptors, and in situ hybridization indicated that it is associated with pheromone-sensitive neurons [31]. SNMP1 has been confirmed to participate in pheromone signal transduction. The SNMP2 subfamily was first identified from Manduca sexta and associates with pheromone-sensitive neurons, but it was expressed only in sensilla support cells [32].

Locusts and grasshoppers are major economic pests, but their genetic information is lacking, partly because their genomes are often very large. Currently, there are data for more than 100 genomes of Orthoptera species in the Genome Size Database (www.genomesize.com). The known variation in Orthoptera genome size ranges from 1.52 Gb for the cave cricket (Hadenoecus subterraneus) to 16.56 Gb for the mountain grasshopper (Podisma pedestris), with an 11-fold difference in size [33, 34]. Sequencing large genomes has higher requirements for sequencing technology and for human and material resources than sequencing small genomes, which explains why only the migratory locust Locusta migratoria (genome size is ~ 6.3 Gb) of Orthoptera has a complete genome sequence thus far [35]. Transcriptomic approaches offer an alternative to genomic approaches; transcriptomic approaches can generate almost all transcripts of a specific organ or tissue of a certain species in a comprehensive and rapid manner, and most molecular mechanisms of different biological processes are also elucidated in the transcriptome [36]. Transcriptomic sequencing results have less data and are more convenient for analysis than genomic sequencing results.

In recent years, there have been increasing reports on the transcriptome of the order Orthoptera, and such reports potentially provide resources for advancing the postgenomic research of Orthoptera insects; however, limited olfaction research is available for these species. Thus far, olfaction studies have been published regarding only L. migratoria [37,38,39,40,41], Oedaleus asiaticus [42], Oedaleus infernalis [43], Ceracris kiangsu [44] and Schistocerca gregaria [16, 45,46,47], and most investigations have focused on OBP genes. Here, we present a de novo transcriptome assembly for the bamboo grasshopper C. nigricornis and identified 112 putative olfactory-related genes comprising 20 OBPs, 10 CSPs, 71 ORs, 8 IRs and 3 SNMPs. Then, we evaluated the distribution of the expression patterns of these genes in different tissues of female and male adults by transcriptome analyses and quantitative real-time PCR (qRT-PCR). Our study provides the foundation for further studies of the molecular mechanism regulating the olfactory system in C. nigricornis.

Results

Sequencing and de novo assembly

The transcriptomes of the antennae (A), head (antennae were cut off; H), abdomen-thorax (T), legs (L) and wings (W) of female and male C. nigricornis were separately sequenced using the Illumina HiSeq X Ten platform. After the low-quality reads were filtered, a total of 223.76 Gb clean data were obtained from all 30 tissue samples, and the clean data of each tissue sample reached 6.30 Gb with a Q30 percentage greater than 94% (Additional file 1: Table S1). After all of the samples were assembled, 43,603 unigenes were generated with an N50 length of 2235 bp, and among them, 20,914 unigenes (47.96%) had a length of over 1 kb (Additional file 1: Table S2). To assess the transcriptome assembly completeness, the benchmarking sets of universal single-copy orthologs (BUSCO) v3.0.2 completeness assessment tool was used together with the Insecta odb9 database with 1658 reference genes [48]. The result had a completeness score of 89.1%, a fragmented score of 2.5% and a missing BUSCO score of only 8.4% (Additional file 1: Table S3).

Functional annotation

A total of 24,832 (56.95%), 15,750 (36.12%), 13,150 (30.16%), 11,503 (26.38%), 18,100 (41.51%), 26,933 (61.77%), 22,295 (51.13%) and 12,102 (27.75%) unigenes were successfully annotated to National Center for Biotechnology Information (NCBI) nonredundant protein sequences (NR), Swiss-Prot (a manually annotated and reviewed protein sequence database), Kyoto Encyclopedia of Genes and Genomes (KEGG), euKaryotic Orthologous Groups of proteins (KOG), Clusters of Orthologous Groups of proteins (COG), EggNOG (A database of orthologous groups and functional annotation), Protein family (Pfam) and Gene Ontology (GO) databases, respectively, which covered a total of 29,122 (66.79%) unigenes (Additional file 1: Table S4).

A query of the NR database indicated that a high percentage of C. nigricornis sequences closely matched insect sequences (11,695, 74.94%). Among these sequences, the highest match sequence and percentage was identified with sequences of Cryptotermes secundus (2548, 21.79%), followed by sequences of Zootermopsis nevadensis (1765, 15.09%), Nilaparvata lugens (883, 7.55%), L. migratoria (729, 6.23%), Rhagoletis zephyria (223, 1.91%), Lasius niger (198, 1.69%), Bemisia tabaci (195, 1.67%), S. gregaria (144, 1.23%), and Tribolium castaneum (141, 1.21%) (Additional file 2: Figure S1).

Blast2GO was applied to classify the functional groups of all unigenes of C. nigricornis. As one unigene could align to multiple GO categories, the assigned GO term was apparently larger than the annotated loci. In total, 12,102 unigenes were classified into at least one of the three main GO categories: 8661 (71.57%) were assigned to biological process, 6828 (56.42%) were assigned to cellular component and 9766 (80.70%) were assigned to molecular function. For the biological process (including 22 subcategories) category, metabolic process (5787 unigenes), cellular process (5539 unigenes) and single-organism process (3361 unigenes) were the most highly enriched GO terms, whereas cell (4705 unigenes), cell part (4674 unigenes), and organelle (3263 unigenes) were the most predominant GO terms in the cellular component (including 17 subcategories) category. For molecular function (including 16 subcategories), the most represented GO terms were catalytic activity (5846 unigenes), binding (5140 unigenes) and structural molecule activity (1190 unigenes) (Additional file 2: Figure S2).

Olfactory-related gene identification

Based on functional annotation and tBLASTn results, a total of 20 candidate OBP genes (CnigOBP1–20) were identified in the transcriptome of C. nigricornis (Table 1). All of these candidate OBP genes had six conserved cysteine residues (Additional file 2: Figure S3). Among the 20 OBP genes, 17 had intact open reading frames (ORFs) with lengths ranging from 408 bp to 816 bp. Except for CnigOBP11 and CnigOBP16, all full-length OBPs had a predicted signal peptide (a signature of secretory proteins) at the N-terminal region (Table 1). The conserved domain prediction of these candidate OBP genes showed that all of them had the domain of a pheromone/general odorant-binding protein (PhBP or PBP_GOBP) (InterPro: IPR006170) (Additional file 1: Table S5).

Table 1 Summary of odorant binding proteins (OBPs) identified in C. nigricornis

A total of 10 candidate CSP genes (CnigCSP1–10) were identified from the transcriptomes of different tissues of C. nigricornis (Table 2). All of these candidate CSP genes had four conserved cysteine residues and a conserved OS-D domain (InterPro: IPR005055); however, for CnigCSP9, two OS-D domains were identified by the conserved domain prediction (Additional file 2: Figure S4 and Additional file 1: Table S6). Among the 10 CSP genes, six CSP genes had full-length ORFs, the remaining CSP genes were incomplete due to a lack of a 5′ or 3′ terminus. The SignalP tests showed that all full-length CSP genes had a predicted signal peptide (Table 2).

Table 2 Summary of chemosensory proteins (CSPs) identified in C. nigricornis

In the transcriptomes of C. nigricornis, 71 candidate OR genes were identified, including 70 conventional ORs (CnigOR1–70) and one ORco (CnigORco). Of these, only 15 candidate OR genes had complete ORFs with lengths longer than 394 amino acids and had 4–7 TMDs (Table 3). Eight candidate IR genes were identified (CnigIR1–5, CnigIR8a, CnigIR25a and CnigIR76b), and six contained a full-length ORF with lengths longer than 319 amino acids (Table 4). Three candidate SNMPs were identified and named CnigSNMP1, CnigSNMP2 and CnigSNMP2a. Only CnigSNMP1 had complete ORFs encoding 532 amino acids (Table 5).

Table 3 Summary of odorant receptors (ORs) identified in C. nigricornis
Table 4 Summary of ionotropic receptors (IRs) identified in C. nigricornis
Table 5 Summary of sensory neuron membrane proteins (SNMPs) identified in C. nigricornis

Homology relationship of olfactory-related genes

To reveal the homology relationships of all olfactory related genes of C. nigricornis with other insect gene sets, we conducted phylogenetic analyses based on the amino acid sequences of 121 OBPs from nine species (Additional file 2: Figure S5), 87 CSPs from seven species (Additional file 2: Figure S6), 293 ORs from three species (Additional file 2: Figure S7), 115 IRs from five species (Additional file 2: Figure S8), and 24 SNMPs from nine species (Additional file 2: Figure S9), respectively. All members of CSPs (Additional file 2: Figure S6), IRs (Additional file 2: Figure S8), and SNMPs (Additional file 2: Figure S9) show orthologous relationships with the counterparts from other orthopteran species. For OBPs, CnigOBP6 and CnigOBP8 are paralogous and may arouse by a recent gene duplication event; the remaining 18 CnigOBPs have orthologous relationships with the other orthopteran species (Additional file 2: Figure S5). For 71 CnigORs, 61 show orthologous relationships with orthopteran species, the other 10 CnigORs (CnigOR5/48, CnigOR17/18, CnigOR19/20, CnigOR35/36, CnigOR40/41) show 2:1 orthologous relationships with the orthopteran species, indicating in-paralogous or out-paralogous relationships among these CnigORs pairs.

Tissue-specific expression analyses by RNA-Seq

To fully understand the differential expression patterns of olfactory-related genes in different tissues, the Illumina reads of each RNA sample were mapped to the reference transcripts to determine the expression quantity. The average number of mapped reads was 75.26% (Additional file 1: Table S1). Fragments per kilobase per million mapped fragments (FPKM) values [49] were determined to measure the gene expression levels.

To detect olfactory-related differentially expressed genes (DEGs) of different tissues, we considered four combinations for comparison: antennae vs. head (A vs H), antennae vs. abdomen-thorax (A vs T), antennae vs. leg (A vs L) and antennae vs. wing (A vs W). A total of 21,484 transcripts were DEGs, and 6357 of them were assigned to GO terms. Among the 6357 DEGs, we found that cellular process, metabolic process and single-organism process represented a high percentage of the biological process category, catalytic activity and binding represented the majority of the molecular function category, cell and cell part represented the greatest proportion of the cellular component category (Fig. 1a). In the biological process category, significantly enriched GO terms were mainly associated with chemosensory perception, such as sensory perception of chemical stimulus, sensory perception of smell, detection of chemical stimulus involved in sensory perception of smell, detection of stimulus involved in sensory perception and sensory perception (Fig. 1b). This expression pattern suggests that chemosensory perception is differentially expressed in different tissues. To better understand these differences, we manually inspected the transcription of genes encoding binding proteins (OBPs and CSPs) and chemosensory membrane proteins (ORs, IRs and SNMPs) to find DEGs in different tissues. The hierarchical cluster analysis showed that all differentially expressed olfactory-related genes were clustered into three clusters (Fig. 1c). In Cluster I and Cluster III, a total of 105 genes had the highest expression in antennae, but the log2FPKM values of Cluster III were higher than those of Cluster I. In Cluster II, the expression detected in the head was higher than that detected in the antennae, legs, wings and abdomen-thorax.

Fig. 1
figure1

Differentially expressed genes (DEGs) analysis of C. nigricornis. a Gene ontology (GO) classifications of DEGs in C. nigricornis. b GO enrichment analysis in the top 25 DEGs of biological process category. c K-means clustering of 111 olfactory-related DEGs. A: antennae; H: head (antennae were cut off); L: leg; W: wing; T: abdomen-thorax. Log10 (FPKM+ 1) values were used, and FPKM values were the average values of each tissues, including female three biological repeats and male three biological repeats

To clearly describe the DEGs of each olfactory-related gene, we performed cluster analysis on their DEGs and represented in heatmap (Fig. 2). Based on the expression levels in different tissues, most of the differentially expressed OBP genes were highly expressed in the antennae and head tissues. Among them, 12 candidate OBPs (CnigOBP2/3/4/6/7/8/11/13/14/15/16/20) showed antennae-specific expression, and seven candidate OBPs (CnigOBP1/9/10/12/17/18/19) showed head-specific expression (Fig. 2b). The expression analysis showed that except for CnigCSP5 and CnigCSP8, all of the differentially expressed CSP genes were highly expressed in the antennae. Among those genes highly expressed in the antennae, CnigCSP6 and CnigCSP9 were also highly expressed in the legs, CnigCSP1 and CnigCSP10 were also highly expressed in the wings, and CnigCSP7 was also highly expressed in the abdomen-thorax. CnigCSP5 was relatively highly expressed in the legs and wings, and CnigCSP8 was highly expressed in the head (Fig. 2c). Except for CnigOR69 that showed relatively high expression levels in the legs, all differentially expressed ORs showed antennal-specific or antennal-biased expression (Fig. 2a). The IRs showed a similar expression profile to the ORs, but more IRs than ORs were detected in the legs and wings (Fig. 2d). All SNMPs were highly expressed in the antennae (Fig. 2e).

Fig. 2
figure2

Hierarchical clustering of differentially expressed olfactory-related genes with difference tissues in C. nigricornis. a odorant receptors (ORs); b odorant binding proteins (OBPs); c chemosensory proteins (CSPs); d ionotropic receptors (IRs); e sensory neuron membrane proteins (SNMPs). A: antennae; H: head (antennae were cut off); L: leg; W: wing; T: abdomen-thorax. Log10 (FPKM+ 1) values were used, and FPKM values were the average values of each tissues including female three biological repeats and male three biological repeats

Sex- and tissue-specific expression analyses by qRT-PCR

To explore the tissue distribution expression pattern of olfactory-related genes in female and male adults of C. nigricornis and to test the RNA-Seq results, we investigated the expression patterns of 12 OBPs, all CSPs and 12 ORs with qRT-PCR analyses. The results showed that the expression levels of the tested genes in different tissues were mostly consistent with the results of RNA-Seq. Based on the results of qRT-PCR assays, six OBP genes (CnigOBP3/4/6/14/16/20) had significantly increased expressed in the antennae, and CnigOBP3 and CnigOBP20 were predominantly expressed in the male antennae (P < 0.05) (Fig. 3). CnigOBP1, CnigOBP12 and CnigOBP19 were predominantly expressed in the head, and CnigOBP1 was significantly more highly expressed in the male head than in the female head among the three head-biased OBPs (Fig. 3). In addition, CnigOBP15 was highly expressed in the legs, CnigOBP11 was highly expressed in the antennae and legs, CnigOBP13 was highly expressed in three tissues (antennae, head and leg), and all 12 OBPs had no or little expression in the wings and abdomen-thorax (Fig. 3).

Fig. 3
figure3

Quantitative real-time PCR analysis of relative expression levels of odorant binding proteins (OBPs) from C. nigricornis. A: antennae; H: head (antennae were cut off); L: leg; W: wing; T: abdomen-thorax. The β-actin was used as the reference gene and CnigOBP3 in female antennae as a positive control. The error bars represent the standard error of three independent experiments. Different small letters above bars indicate significant differences among different tissues (P < 0.05). * indicates significant difference between both sexes in the same tissue (P < 0.05), and ns indicates no significant difference

The qRT-PCR results revealed that all 10 CSPs had significant differences in the expression levels among different tissues. Four CSP genes (CnigCSP1/2/3/4) were expressed more highly in the antennae than in other tissues (P < 0.05), and except for CnigCSP3, all antennae-biased CSPs were significantly more highly expressed in the male antennae than in the female antennae (Fig. 4). CnigCSP8 was predominantly expressed in the head, while CnigCSP5 was highly expressed in the legs. Meanwhile, there were significant differences in the expression levels of CnigCSP8 in the head and CnigCSP5 in the legs between males and females, and CnigCSP8 and CnigCSP5 were highly expressed in the female head and leg tissues, respectively (Fig. 4). The remaining four CSPs (CnigCSP6/7/9/10) were highly expressed in more than two tissues; among them, CnigCSP6, CnigCSP7 and CnigCSP9 were more highly expressed in the antennae and legs than in other tissues, and CnigCSP10 was more highly expressed in the antennae, head and legs than in the wings and abdomen-thorax (Fig. 4).

Fig. 4
figure4

Quantitative real-time PCR analysis of relative expression levels of chemosensory proteins (CSPs) from C. nigricornis. A: antennae; H: head (antennae were cut off); L: leg; W: wing; T: abdomen-thorax. The β-actin was used as the reference gene and CnigCSP1 in female antennae as a positive control. The error bars represent the standard error of three independent experiments. Different small letters above bars indicate significant differences among different tissues (P < 0.05). * indicates significant difference between both sexes in the same tissue (P < 0.05), and ns indicates no significant difference

All randomly selected OR genes were strongly expressed in the antennae, whereas they were not or were more faintly expressed in other tissues (Fig. 5). Among the 12 ORs, four CnigOR genes (CnigOR13/18/R37/70) were more highly expressed in the male antennae than in the female antennae, four CnigOR genes (CnigOR10/20/32/35) were more highly expressed in the female antennae than in the male antennae (P < 0.05), and the remaining four CnigOR genes (CnigOR6/55/60/65) were not significantly different between the male and female antennae (Fig. 5).

Fig. 5
figure5

Quantitative real-time PCR analysis of relative expression levels of odorant receptors (ORs) from C. nigricornis. A: antennae; H: head (antennae were cut off); L: leg; W: wing; T: abdomen-thorax. The β-actin was used as the reference gene and CnigOR37 in female antennae as a positive control. The error bars represent the standard error of three independent experiments. Different small letters above bars indicate significant differences among different tissues (P < 0.05). * indicates significant difference between both sexes in the same tissue (P < 0.05), and ns indicates no significant difference

Discussion

Orthoptera insects have an exclusive evolutionary trend in their olfactory system, which is an excellent model for studying the olfactory evolution of insects [50]. However, few investigations have focused on the molecular basis of olfaction in these species. In this study, we sequenced and analyzed the transcriptomes of C. nigricornis from the antennae, head (antennae were cut off), legs, wings and abdomen-thorax, and the results of this study further enriched the molecular biological foundation of C. nigricornis. Furthermore, we also identified 112 candidate olfactory-related genes in the transcriptomes of C. nigricornis for the first time, including 20 OBPs, 10 CSPs, 71 ORs, eight IRs and three SNMPs. The numbers of candidate olfactory-related genes identified in the transcriptomes of C. nigricornis were similar to the numbers of candidate olfactory-related genes identified in the transcriptomes of O. asiaticus (15 OBPs, 60 ORs, 6 IRs and 3 SNMPs) [42, 51], but less than that identified in the genome of L. migratoria (22 OBPs, 142 ORs and 32 IRs) [35, 52]. One possible reason for this finding is that the olfactory-related genes in C. nigricornis and O. asiaticus were identified from transcriptome data, and these genes are usually expressed at low levels in transcriptome studies. One limiting factor in transcriptome sequencing is the inaccurate detection of genes with low transcriptional abundance, which may lead to the deletion of genes with relatively low expression levels [53]. However, the sequencing method or depth used in the genomic data of L. migratoria may have allowed the detection of genes with lower expression levels.

It has been demonstrated that OBPs increase the sensitivity of odor for insects and mediate the recognition and discrimination of odor compounds [54]. In this study, we identified 20 OBPs in C. nigricornis, differing from 22, 18, 15 and 14 OBPs identified in L. migratoria [52], O. infernalis [43], O. asiaticus [42] and S. gregaria [55], respectively. This may reflect that different insects evolved different physiological behaviors in the process of adapting to various environments, which might lead to divergent evolutionary trajectories of olfactory genes of the same ancestry, resulting in different functional genes among species [18]. Orthopteran insects possess a significantly smaller expansion of the OBP family than Dipteran insects, such as D. melanogaster [56] containing 51 OBPs and Aedes aegypti [57] containing 66 OBPs. However, this reduced number of OBPs seems to be balanced by the expansion of the CSP family, of which a large number of CSPs (70) were reported in the oriental locust L. migratoria [58]. This finding indicates the specific physiological and evolutionary level of the Orthopteran olfactory system.

Based on the RNA-Seq and qRT-PCR results, most OBPs identified in C. nigricornis were higher expressed in antennae than in other tissues, indicating that these OBPs might play roles in the recognition of sex pheromones and host volatile compounds, as in other insect species (Figs. 2b and 3) [59, 60]. Moreover, the expression profiles of most OBPs in the antennae showed male-biased expression, suggesting their possible crucial roles in detecting female pheromone and mating behavior (Fig. 3) [19]. However, we also found that several OBPs were more highly expressed in the head and legs, suggesting that these OBPs might participate in other physiological functions [43]. Among them, CnigOBP1, CnigOBP12 and CnigOBP19 were much more highly expressed in head tissues than in other tissues. The head tissues used in this study still contained other chemosensory organs, such as labium, maxillary palp, labial palp and proboscis, even though their antennae were cut off. Previous studies have suggested that the OBPs expressed in mouthparts and proboscis are associated with gustatory responses [61, 62]. This indicated that CnigOBP1, CnigOBP12 and CnigOBP19 might play important roles in taste functions.

CSPs have similar functions to OBPs in chemical communication in insects, as they bind small molecules, such as pheromones and odorants, and transport them to chemoreceptors [21]. However, CSPs usually exhibit much broader expression profiles than OBPs both in olfactory tissues and nonolfactory tissues and play multiple important roles in biological processes such as ecdysis, leg regeneration and embryo maturation [15, 25, 63]. An interesting example of the diverse roles of CSPs was observed in the locust L. migratoria, in which several CSPs showed increased expression in gregarious locusts, but these CSPs displayed contrasting expression trends when the locusts underwent the physiological transformation from the ‘gregarious’ phase to the ‘solitary’ phase, which indicated that CSPs are one factor triggering locust phase shift [23]. Our results show that CnigCSP1, CnigCSP2, CnigCSP3 and CnigCSP4 were antennae-enriched, and these four CSPs might be involved in the chemosensory process (Fig. 4). In addition, in nonolfactory tissues, the expression profiles of C. nigricornis CSPs were much broader than those of C. nigricornis OBPs (Figs. 3 and 4). CnigCSP6/7/9 showed higher expression in both antennae and legs, CnigCSP5 was mainly expressed in the legs, CnigCSP8 was predominantly expressed in the head and CnigCSP10 was highly expressed in the antennae, head and legs. We speculate that the broad and diverse expression patterns of C. nigricornis CSPs might have other crucial physiological functions, which need further functional verification.

ORs are usually expressed in the dendrites of antennal sensilla and act as biotransducers to convert chemical signals of odorant molecules into electrical signals [64]. To transmit odor-evoked signals, two types of ORs are required: one is a highly conserved ORco, and the other is a specific OR, which varies according to ORN type [65, 66]. In this study, we identified 70 specific ORs and a conserved ORco, and the RNA-Seq results showed that all ORs were highly expressed in the antennae, except CnigOR69, which were relatively highly expressed in nonolfactory leg tissues (Fig. 2a). The antennae-biased ORs may play an important role in odorant reception in antennae or be involved in the olfactory sensing process [67], and the nonolfactory tissues highly expressed ORs, which suggests that they may participate in other physiological processes in addition to olfaction [53]. Moreover, we also investigated the distribution of OR expression patterns in female and male tissues with qRT-PCR analyses. The expression profiles showed that all 12 randomly selected ORs were strongly expressed in the antennae, which was consistent with the results of RNA-Seq. Among the 12 ORs, four CnigOR genes (CnigOR13/18/R37/70) showed male-biased expression, suggesting that they may play a role in female pheromone detection and mating behavior, whereas four CnigOR genes (CnigOR10/20/32/35) displayed female-biased expression predicted to function in oviposition-related odorant detection [68, 69] (Fig. 5).

IRs, a new subfamily of ORs, are involved in not only olfaction but also gustation, hygrosensation and thermosensation, which are widely distributed throughout the body, including the labellum, leg, pharynx, and wing sensilla [70,71,72,73,74]. We identified eight IRs in C. nigricornis, of which CnigIR2 were not differentially expressed in different tissues, suggesting that this IR may have multiple functions. Except for CnigIR5, all differentially expressed IRs were strongly expressed in antennae (Fig. 2d). Recent studies have shown that antennae-enriched IRs play important roles in odor and thermosensation [75], the antennae-enriched IRs of C. nigricornis that were identified in our study may have similar roles.

We also identified three SNMPs, and the FPKM-value showed that all SNMPs were strongly expressed in antennae, indicating that they may be involved in the process of olfaction. However, CnigSNMP1 and CnigSNMP2 were also broadly expressed in nonolfactory tissues, including head, leg, wing and abdomen-thorax tissues (Fig. 2e). Similar patterns were also observed in other studies of D. melanogaster, A. aegypti, Spodoptera litura and Cnaphalocrocis medinalis [76,77,78,79,80]. The broad expression patterns of SNMP1 and SNMP2 imply that they may not only function in odorant perception, but also have different functions specific to the various tissues [78]. The exact role for CnigSNMP1 and CnigSNMP2 in olfactory and nonolfactory tissues remains inconclusive and requires further investigation.

Conclusions

The 112 candidate olfactory-related genes that we identified from C. nigricornis compose the comprehensive list of olfactory-related genes in this bamboo pest of C. nigricornis. Sex- and tissue-specific expression profiling revealed that most of the candidate olfactory-related genes were antennae-enriched, but some were nonantennae-enriched, and some were sex-biased, indicating their different roles in the olfactory system of C. nigricornis. Our results provide a foundation to facilitate functional studies of these olfactory-related genes in C. nigricornis at the molecular level.

Methods

Species and tissues collection

All adult male and adult female specimens of C. nigricornis were collected in Xi’an, China, in August 2018. Different tissues, including antennae, head (antennae were cut off), abdomen-thorax, legs and wings, were dissected from female and male specimens, respectively. Each type of tissue sample was collected from nine individuals as an independent biological replicate for transcriptome sequencing. To minimize biological variance, three independent biological replicates were performed for each tissue sample. All tissue samples were immediately placed into enzyme-free centrifuge tubes, which were immersed in liquid nitrogen and stored at − 80 °C until further use.

RNA extraction, library construction and Illumina sequencing

Total RNA was extracted from each frozen tissue using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the protocol recommended by the manufacturer. RNA integrity and purity were determined with an Agilent Bioanalyzer 2100 (Agilent Technologies, CA, USA) and a Nano Drop 2000 (Thermo Scientific, Wilmington, DE, USA). The cDNA library construction and Illumina sequencing were carried out by the Biomarker Technology Company, Beijing, China. First, qualified RNA from each tissue was mixed in equal amounts. Then, the cDNA library was prepared using the NEBNext Ultra™ RNA Library Prep Kit for Illumina (NEB, USA), and the Agilent Bioanalyzer 2100 system was used to assess the library. Finally, the amplified mRNA library was sequenced on the Illumina HiSeq X Ten platform to generate 150 bp paired-end reads.

De novo assembly and functional annotation

Prior to assembly, low-quality reads and adaptor reads were removed from raw data through in-house Perl scripts to obtain clean reads. Then, the clean reads were de novo assembled using Trinity v2.4.0 [81] with default parameters except that ‘min_kmer_cov’ set to 2. Next, cd-hit-est v4.6.1 software [82] was used to remove duplicate and highly similar sequences to obtain nonredundant unigenes. Furthermore, the completeness of the nonredundant unigenes was assessed using BUSCO v3.0.2 software [48], and the insect lineage database (insecta_orthoDB9, created 13/02/2016) was used as a proxy for the minimum completeness assessment library.

BLAST v2.2.23 software [83] was used to search and annotate all assembled unigenes with the publicly available protein databases, including NR, Swiss-Prot, KEGG, KOG, COG and EggNOG, and an E-value cut-off of 1.0E-05. KOBAS2.0 was used to obtain KEGG Orthology results of unigenes in the KEGG annotation. Blast2GO [84] was used for GO annotation based on the protein annotation results of the NR database. For the GO classifications, the default parameters were used (E-value <1E-5, an annotation cut-off > 5 and a GO weight > 5). After the amino acid sequences of the unigenes were obtained, HMMER [85] software (E-value ≤1E-10) was used to BLAST search against the Pfam database to obtain the annotation information of unigenes.

Identification of olfactory-related genes

Olfactory-related genes from the transcriptomes of C. nigricornis were identified by the functional annotation results. To obtain more information, a tBLASTn program was also performed using available olfactory proteins from Orthoptera species as queries to identify candidate unigenes. All candidate olfactory genes were manually confirmed using the BLASTx program against the NR database at NCBI (https://blast.ncbi.nlm.nih.gov/Blast.cgi). The ORFs of candidate ORs, IRs, SNMPs, OBPs and CSPs were predicted using ORF Finder (https://www.ncbi.nlm.nih.gov/orffinder/). The N-terminal signal peptides and conserved domains of candidate OBPs and CSPs were predicted by the SignalP V 4.1 program (http://www.cbs.dtu.dk/services/SignalP/) [86] and SMART (simple modular architecture research tool, http://smart.embl.de/) [87], respectively. TMHMM Server Version 2.0 (http://www.cbs.dtu.dk/services/TMHMM/) was used to predict the TMDs of candidate ORs, IRs and SNMPs.

Phylogenetic analysis

Phylogenetic analysis was performed based on the amino acid sequences of C. nigricornis olfactory-related genes using MEGA7.0 software [88]. The maximum likelihood method with Jones-Taylor-Thornton (JTT) model was used to create the phylogenetic trees. Branch support was assessed by bootstrapping with 1000 replicates. All the amino acid sequences used to construct the phylogenetic tree are shown in Additional file 3.

Differential expression analysis

Using the assembled transcriptome as reference sequences, the clean data from various samples of C. nigricornis were mapped back onto the reference sequences using Bowtie2 v2.1.0 software [89]. The unigene expression levels among various samples were estimated by RSEM [90] according to the readcount values of the unigenes for each sample, which were obtained from the mapping results. FPKM values can eliminate the effects of the sequencing depth and transcript length on readcounts. To make the unigene expression levels of various samples comparable, all the readcounts were transformed into FPKM values.

To detect the differential expression of candidate olfactory-related genes among different tissues from transcriptomes of C. nigricornis, we used the antennae as the experimental tissue, and four other tissues (head, abdomen-thorax, leg, and wing) were used as control samples. Therefore, a total of four combinations were compared, including antennae vs. head, antennae vs. abdomen-thorax, antennae vs. leg and antennae vs. wing. The differential expression analysis of each comparison combination was analyzed by the DESeq R package (1.10.1) [91]. The resulting P values were adjusted using Benjamini and Hochberg’s approach for controlling the false discovery rate (FDR). In general, the fold change (log2 ratio) and adjusted P value (FDR) were used as the key indexes for screening DEGs. In this study, |log2 (Fold Change)| > 2 and adjusted P value < 0.05 were used as the thresholds for identifying significant DEGs. The log2-transformed FPKM values of the DEGs were used for hierarchical clustering using the pheatmap R package. GO enrichment analysis of DEGs was implemented by the topGO R package-based Kolmogorov-Smirnov test.

Expression profile analysis

To explore the tissue distribution expression pattern of olfactory-related genes in female and male adults and to verify the accuracy of transcriptome analysis, qRT-PCR analysis was performed to evaluate the expression profiles of putative olfactory-related genes in different tissues of male and female C. nigricornis adults. Total RNA was extracted using TRIzol reagent, and cDNA was synthesized by Goldenstar™ RT6 cDNA Synthesis Mix (TsingKe, Beijing, China) following the protocol recommended by the manufacturer. A total of 12 OBPs, all CSPs and 12 ORs were randomly selected for qRT-PCR analysis. To correct sample-to-sample variation, ß-actin was selected as the reference gene. The specific primers used in the qRT-PCR analysis were designed online (https://www.ncbi.nlm.nih.gov/tools/primer-blast/) (Additional file 1: Table S7). qRT-PCR was performed on the CFX96 Touch Real-Time PCR Detection System (Bio-Rad, USA) using 2 x T5 Fast qPCR Mix (TsingKe, Beijing, China) with the following conditions: initial denaturation at 95 °C for 1 min, followed by 40 cycles of 95 °C for 10 s and 60 °C for 15 s. Finally a melting curve of the PCR products was analyzed to confirm the specificity of the primers. Nontemplate negative control reactions (replacing cDNA with H2O) were included in each experiment. Three biological replicates with three technical replicates were conducted for all experiments.

The results were analyzed using the BIO-RAD CFX Manager software. The relative quantification among all tissues was calculated by the comparative 2-CT method [92]. A one-way analysis of variance (ANOVA) followed by Duncan’s test were performed for the comparative analyses of each gene among different tissues, and Student’s t-test was used to analyze each gene in the same tissue from females and males. Both analyses were performed by SPSS Statistics 18.0 software (SPSS Inc., Chicago, IL, USA).

Availability of data and materials

The authors declare that the data supporting the finding of this study are available in the article and its supplementary information files. The olfactory-related genes we identified from Ceracris nigricornis were all submitted to the National Center for Biotechnology Information (NCBI) (https://www.ncbi.nlm.nih.gov/) with the accession numbers were listed in Tables 1, 2, 3, 4 and 5. The Transcriptome Shotgun Assembly project has been deposited at DDBJ/EMBL/GenBank under the accession GHNZ00000000. The version described in this paper is the first version, GHNZ01000000. The raw data submitted to the NCBI with the accession numbers of SRR9089653, SRR9089654, SRR9089655, SRR9089656, SRR9089657, SRR9089658, SRR9089659, SRR9089660, SRR9089661, SRR9089662, SRR9089663, SRR9089664, SRR9089665, SRR9089666, SRR9089667, SRR9089668, SRR9089669, SRR9089670, SRR9089671, SRR9089672, SRR9089673, SRR9089674, SRR9089675, SRR9089676, SRR9089677, SRR9089678, SRR9089679, SRR9089680, SRR9089681, SRR9089682.

Abbreviations

CSPs:

Chemosensory binding proteins

DEGs:

Differentially expressed genes

FPKM:

Fragments per kilobase per million mapped fragments

GO:

Gene Ontology

IR:

Ionotropic receptor

OBP:

Odorant-binding protein

OR:

Odorant receptor

ORF:

Open reading frame

qRT-PCR:

Quantitative realtime polymerase chain reaction

SNMP:

Sensory neuron membrane protein

References

  1. 1.

    Qiu Z, Liu F, Lu H, Yuan H, Zhang Q, Huang Y. De novo assembly and characterization of the transcriptome of grasshopper Shirakiacris shirakii. Int J Mol Sci. 2016;17(7):1110.

  2. 2.

    Bian L, Sun X, Luo Z, Zhang Z, Chen Z. Design and selection of trap color for capture of the tea leafhopper, Empoasca vitis, by orthogonal optimization. Entomologia Experimentalis Et Applicata. 2014;151(3):247–58.

  3. 3.

    Cai X, Xu X, Bian L, Luo Z, Chen Z. Measurement of volatile plant compounds in field ambient air by thermal desorption-gas chromatography-mass spectrometry. Anal Bioanal Chem. 2015;407(30):9105–14.

  4. 4.

    Jin S, Chen ZM, Backus EA, Sun XL, Xiao B. Characterization of EPG waveforms for the tea green leafhopper, Empoasca vitis Göthe (Hemiptera: Cicadellidae), on tea plants and their correlation with stylet activities. J Insect Physiol. 2012;58(9):1235–44.

  5. 5.

    Song T, Wang K, Peng W, Wang J, Xiao R, Zeng F, Tang Y. Ecological effects of intercropping white clover on tea plantation in a subtropical hilly region. Acta Ecol Sin. 2006;26(11):3647–55.

  6. 6.

    Ji X, Wang M, Jiang L, Han B. Similarity and disparity of arthropod community compositions among ten tea-intercropping plantations with different species of fruits in the Dongtingshan Mountains, Taihu Lake. Chin J Appl Entomol. 2011;48(5):1471–8.

  7. 7.

    Ye H, Cui L, He X, Han B. Effect of intercropping tea with citrus, waxberry, or snake gourd on population density and spatial distribution of the tea green leafhopper and araneids. Acta Ecol Sinica. 2010;30:6019–26.

  8. 8.

    Vieira FG, Foreˆt S, He X, Rozas J, Field LM, Zhou J-J. Unique features of odorant-binding proteins of the parasitoid wasp Nasonia vitripennis revealed by genome annotation and comparative analyses. PLoS One. 2012;7(8):e43034.

  9. 9.

    Deisig N, Kropf J, Vitecek S, Pevergne D, Rouyar A, Sandoz J, Lucas P, Gadenne C, Anton S, Barrozo R. Differential interactions of sex pheromone and plant odour in the olfactory pathway of a male moth. PLoS One. 2012;7(3):e33159.

  10. 10.

    Binyameen M, Hussain A, Yousefi F, Birgersson G, Schlyter F. Modulation of reproductive behaviors by non-host volatiles in the polyphagous Egyptian cotton leafworm, Spodoptera littoralis. J Chem Ecol. 2013;39(10):1273–83.

  11. 11.

    Leal WS. Odorant reception in insects: roles of receptors, binding proteins, and degrading enzymes. Annu Rev Entomol. 2013;58:373–91.

  12. 12.

    Justice R, Dimitratos S, Walter M, Woods D, Biessmann H. Sexual dimorphic expression of putative antennal carrier protein genes in the malaria vector Anopheles gambiae. Insect Mol Biol. 2003;12(6):581–94.

  13. 13.

    Pelosi P, Iovinella I, Felicioli A, Dani FR. Soluble proteins of chemical communication: an overview across arthropods. Front Physiol. 2014;5:320.

  14. 14.

    Vieira FG, Rozas J. 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. 2011;3:476–90.

  15. 15.

    Pelosi P, Iovinella I, Zhu J, Wang G, Dani FR. Beyond chemoreception: diverse tasks of soluble olfactory proteins in insects. Biol Rev. 2018;93(1):184–200.

  16. 16.

    Angeli S, Ceron F, Scaloni A, Monti M, Monteforti G, Minnocci A, Petacchi R, Pelosi P. Purification, structural characterization, cloning and immunocytochemical localization of chemoreception proteins from Schistocerca gregaria. Eur J Biochem. 1999;262(3):745–54.

  17. 17.

    Zeng Y, Yang YT, Wu QJ, Wang SL, Xie W, Zhang YJ. Genome-wide analysis of odorant-binding proteins and chemosensory proteins in the sweet potato whitefly, Bemisia tabaci. Insect Sci. 2019;26(4):620–34.

  18. 18.

    Li X-M, Zhu X-Y, Wang Z-Q, Wang Y, He P, Chen G, Sun L, Deng D-G, Zhang Y-N. Candidate chemosensory genes identified in Colaphellus bowringi by antennal transcriptome analysis. BMC Genomics. 2015;16(1):1028.

  19. 19.

    Brito NF, Moreira MF, Melo AC. A look inside odorant-binding proteins in insect chemoreception. J Insect Physiol. 2016;95:51–65.

  20. 20.

    Wanner KW, Willis LG, Theilmann DA, Isman MB, Feng Q, Plettner E. Analysis of the insect os-d-like gene family. J Chem Ecol. 2004;30(5):889–911.

  21. 21.

    Pelosi P, Calvello M, Ban L. Diversity of odorant-binding proteins and chemosensory proteins in insects. Chem Senses. 2005;30(Suppl 1):i291–2.

  22. 22.

    Gong D-P, Zhang H-j, Zhao P, Lin Y, Xia Q-Y, Xiang Z-H. Identification and expression pattern of the chemosensory protein gene family in the silkworm, Bombyx mori. Insect Biochem Mol Biol. 2007;37(3):266–77.

  23. 23.

    Guo W, Wang X, Ma Z, Xue L, Han J, Yu D, Kang L. CSP and takeout genes modulate the switch between attraction and repulsion during behavioral phase change in the migratory locust. PLoS Genet. 2011;7(2):e1001291.

  24. 24.

    Zhang Y-N, Jin J-Y, Jin R, Xia Y-H, Zhou J-J, Deng J-Y, Dong S-L. Differential expression patterns in chemosensory and non-chemosensory tissues of putative chemosensory genes identified by transcriptome analysis of insect pest the purple stem borer Sesamia inferens (Walker). PLoS One. 2013;8(7):e69715.

  25. 25.

    Cheng D, Lu Y, Zeng L, Liang G, He X. Si-CSP9 regulates the integument and moulting process of larvae in the red imported fire ant, Solenopsis invicta. Sci Rep. 2015;5:9245.

  26. 26.

    Clyne PJ, Warr CG, Freeman MR, Lessing D, Kim J, Carlson JR. A novel family of divergent seven-transmembrane proteins: candidate odorant receptors in Drosophila. Neuron. 1999;22(2):327–38.

  27. 27.

    Benton R, Vannice KS, Gomez-Diaz C, Vosshall LB. Variant ionotropic glutamate receptors as chemosensory receptors in Drosophila. Cell. 2009;136(1):149–62.

  28. 28.

    Depan Cao YL, Wei J, Liao X, Walker WB, Li J, Wang G. Identification of candidate olfactory genes in Chilo suppressalis by antennal transcriptome analysis. Int J Biol Sci. 2014;10(8):846.

  29. 29.

    Zhao H, Du Y, Gao P, Wang S, Pan J, Jiang Y. Antennal transcriptome and differential expression analysis of five chemosensory gene families from the Asian honeybee Apis cerana cerana. PLoS One. 2016;11(10):e0165374.

  30. 30.

    Gomez-Diaz C, Bargeton B, Abuin L, Bukar N, Reina JH, Bartoi T, Graf M, Ong H, Ulbrich MH, Masson J-F. A CD36 ectodomain mediates insect pheromone detection via a putative tunnelling mechanism. Nat Commun. 2016;7:11866.

  31. 31.

    Forstner M, Gohl T, Gondesen I, Raming K, Breer H, Krieger J. Differential expression of SNMP-1 and SNMP-2 proteins in pheromone-sensitive hairs of moths. Chem Senses. 2008;33(3):291–9.

  32. 32.

    Rogers ME, Krieger J, Vogt RG. Antennal SNMPs (sensory neuron membrane proteins) of lepidoptera define a unique family of invertebrate CD36-like proteins. J Neurobiol. 2001;49(1):47–61.

  33. 33.

    Camacho JM, Ruiz-Ruano F, Martín-Blázquez R, López-León M, Cabrero J, Lorite P, Cabral-de-Mello D, Bakkali M. A step to the gigantic genome of the desert locust: chromosome sizes and repeated DNAs. Chromosoma. 2015;124(2):263–75.

  34. 34.

    Dufresne F, Jeffery N. A guided tour of large genome size in animals: what we know and where we are heading. Chromosom Res. 2011;19(7):925–38.

  35. 35.

    Wang X, Fang X, Yang P, Jiang X, Jiang F, Zhao D, Li B, Cui F, Wei J, Ma C. The locust genome provides insight into swarm formation and long-distance flight. Nat Commun. 2014;5:2957.

  36. 36.

    Hittinger CT, Johnston M, Tossberg JT, Rokas A. Leveraging skewed transcript abundance by RNA-Seq to increase the genomic depth of the tree of life. Proc Natl Acad Sci. 2010;107(4):1476–81.

  37. 37.

    Ban L, Scaloni A, D'ambrosio C, Zhang L, Yan Y, Pelosi P. Biochemical characterization and bacterial expression of an odorant-binding protein from Locusta migratoria. Cell Mol Life Sci. 2003;60(2):390–400.

  38. 38.

    Jin X, Brandazza A, Navarrini A, Ban L, Zhang S, Steinbrecht RA, Zhang L, Pelosi P. Expression and immunolocalisation of odorant-binding and chemosensory proteins in locusts. Cell Mol Life Sci. 2005;62(10):1156–66.

  39. 39.

    Yu Y, Cui X, Jiang Q, Jin X, Guo Z, Zhao X, Bi Y, Zhang L. New isoforms of odorant-binding proteins and potential semiochemicals of locusts. Arch Insect Biochem Physiol. 2007;65(1):39–49.

  40. 40.

    Yu Y, Zhang S, Zhang L, Zhao X. Developmental expression of odorant-binding proteins and chemosensory proteins in the embryos of Locusta migratoria. Arch Insect Biochem Physiol. 2009;71(2):105–15.

  41. 41.

    Li J, Zhang L, Wang X. An odorant-binding protein involved in perception of host plant odorants in locust Locusta migratoria. Arch Insect Biochem Physiol. 2016;91(4):221–9.

  42. 42.

    Zhang S, Pang B, Zhang L. Novel odorant-binding proteins and their expression patterns in grasshopper, Oedaleus asiaticus. Biochem Biophys Res Commun. 2015;460(2):274–80.

  43. 43.

    Zhang Y, Tan Y, Zhou X-R, Pang B-P. A whole-body transcriptome analysis and expression profiling of odorant binding protein genes in Oedaleus infernalis. Comp Biochem Physiol D Genomics Proteomics. 2018;28:134–41.

  44. 44.

    Li R, Jiang G-F, Dong S-Y. A head transcriptome provides insights into odorant binding proteins of the bamboo grasshopper. Genes Genomics. 2018;40(9):991–1000.

  45. 45.

    Tomaselli S, Crescenzi O, Sanfelice D, Wechselberger R, Angeli S, Scaloni A, Boelens R, Tancredi T, Pelosi P, Picone D. Solution structure of a chemosensory protein from the desert locust Schistocerca gregaria. Biochemistry. 2006;45(35):10606–13.

  46. 46.

    Yang Y, Krieger J, Zhang L, Breer H. The olfactory co-receptor Orco from the migratory locust (Locusta migratoria) and the desert locust (Schistocerca gregaria): identification and expression pattern. Int J Biol Sci. 2012;8(2):159.

  47. 47.

    Guo M, Krieger J, Große-Wilde E, Mißbach C, Zhang L, Breer H. Variant ionotropic receptors are expressed in olfactory sensory neurons of coeloconic sensilla on the antenna of the desert locust (Schistocerca gregaria). Int J Biol Sci. 2014;10(1):1.

  48. 48.

    Simão FA, Waterhouse RM, Ioannidis P, Kriventseva EV, Zdobnov EM. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics. 2015;31(19):3210–2.

  49. 49.

    Mortazavi A, Williams BA, Mccue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5(7):621.

  50. 50.

    Hansson BS, Stensmyr MC. Evolution of insect olfaction. Neuron. 2011;72(5):698–711.

  51. 51.

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

  52. 52.

    Wang Z, Yang P, Chen D, Jiang F, Li Y, Wang X, Kang L. Identification and functional analysis of olfactory receptor family reveal unusual characteristics of the olfactory system in the migratory locust. Cell Mol Life Sci. 2015;72(22):4429–43.

  53. 53.

    Li K, Wei H, Shu C, Zhang S, Cao Y, Luo C, Yin J. Identification and comparison of candidate odorant receptor genes in the olfactory and non-olfactory organs of Holotrichia oblita Faldermann by transcriptome analysis. Comp Biochem Physiol Part D Genomics Proteomics. 2017;24:1–11.

  54. 54.

    Laughlin JD, Ha TS, Jones DN, Smith DP. Activation of pheromone-sensitive neurons is mediated by conformational activation of pheromone-binding protein. Cell. 2008;133(7):1255–65.

  55. 55.

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

  56. 56.

    Hekmat-Scafe DS, Scafe CR, McKinney AJ, Tanouye MA. Genome-wide analysis of the odorant-binding protein gene family in Drosophila melanogaster. Genome Res. 2002;12(9):1357–69.

  57. 57.

    Zhou JJ, He XL, Pickett J, Field L. Identification of odorant-binding proteins of the yellow fever mosquito Aedes aegypti: genome annotation and comparative analyses. Insect Mol Biol. 2008;17(2):147–63.

  58. 58.

    Zhou X-H, Ban L-P, Iovinella I, Zhao L-J, Gao Q, Felicioli A, Sagona S, Pieraccini G, Pelosi P, Zhang L. Diversity, abundance, and sex-specific expression of chemosensory proteins in the reproductive organs of the locust Locusta migratoria manilensis. Biol Chem. 2013;394(1):43–54.

  59. 59.

    Swarup S, Williams TI, Anholt RR. Functional dissection of odorant binding protein genes in Drosophila melanogaster. Genes Brain Behav. 2011;10(6):648–57.

  60. 60.

    Guo W, Ren D, Zhao L, Jiang F, Song J, Wang X, Kang L. Identification of odorant-binding proteins (OBPs) and functional analysis of phase-related OBPs in the migratory locust. Frontiers in physiology. 2018;9:984.

  61. 61.

    DelCampo ML, Palmer S, Caillaud M. Characterization of a new odorant binding protein gene in gustatory organs of Manduca sexta larvae (Lepidoptera: Sphingidae). Ann Entomol Soc Am. 2011;104(2):319–25.

  62. 62.

    Shanbhag S, Park S-K, Pikielny C, Steinbrecht R. Gustatory organs of Drosophila melanogaster: fine structure and expression of the putative odorant-binding protein PBPRP2. Cell Tissue Res. 2001;304(3):423–37.

  63. 63.

    Maleszka J, Foret S, Saint R, Maleszka R. RNAi-induced phenotypes suggest a novel role for a chemosensory protein CSP5 in the development of embryonic integument in the honeybee (Apis mellifera). Dev Genes Evol. 2007;217(3):189–96.

  64. 64.

    Zhang D-D, Löfstedt C. Moth pheromone receptors: gene sequences, function, and evolution. Front Ecol Evol. 2015;3:105.

  65. 65.

    Bengtsson JM, Trona F, Montagné N, Anfora G, Ignell R, Witzgall P, Jacquin-Joly E. Putative chemosensory receptors of the codling moth, Cydia pomonella, identified by antennal transcriptome analysis. PLoS One. 2012;7(2):e31620.

  66. 66.

    Jones PL, Pask GM, Rinker DC, Zwiebel LJ. Functional agonism of insect odorant receptor ion channels. Proc Natl Acad Sci. 2011;108(21):8821–5.

  67. 67.

    Leal WS, Choo Y-M, Xu P, da Silva CS, Ueira-Vieira C. Differential expression of olfactory genes in the southern house mosquito and insights into unique odorant receptor gene isoforms. Proc Natl Acad Sci. 2013;110(46):18704–9.

  68. 68.

    Nakagawa T, Sakurai T, Nishioka T, Touhara K. Insect sex-pheromone signals mediated by specific combinations of olfactory receptors. Science. 2005;307(5715):1638–42.

  69. 69.

    Grosjean Y, Rytz R, Farine J-P, Abuin L, Cortot J, Jefferis GS, Benton R. An olfactory receptor for food-derived odours promotes male courtship in Drosophila. Nature. 2011;478(7368):236.

  70. 70.

    Depetris-Chauvin A, Galagovsky D, Grosjean Y. Chemicals and chemoreceptors: ecologically relevant signals driving behavior in Drosophila. Front Ecol Evol. 2015;3:41.

  71. 71.

    Ganguly A, Pang L, Duong V-K, Lee A, Schoniger H, Varady E, Dahanukar A. A molecular and cellular context-dependent role for Ir76b in detection of amino acid taste. Cell Rep. 2017;18(3):737–50.

  72. 72.

    Koh T-W, He Z, Gorur-Shandilya S, Menuz K, Larter NK, Stewart S, Carlson JR. The Drosophila IR20a clade of ionotropic receptors are candidate taste and pheromone receptors. Neuron. 2014;83(4):850–65.

  73. 73.

    Knecht ZA, Silbering AF, Cruz J, Yang L, Croset V, Benton R, Garrity PA. Ionotropic receptor-dependent moist and dry cells control hygrosensation in Drosophila. Elife. 2017;6:e26654.

  74. 74.

    Ni L, Klein M, Svec KV, Budelli G, Chang EC, Ferrer AJ, Benton R, Samuel AD, Garrity PA. The ionotropic receptors IR21a and IR25a mediate cool sensing in Drosophila. Elife. 2016;5:e13254.

  75. 75.

    Rimal S, Lee Y. The multidimensional ionotropic receptors of Drosophila melanogaster. Insect Mol Biol. 2018;27(1):1–7.

  76. 76.

    Benton R, Vannice KS, Vosshall LB. An essential role for a CD36-related receptor in pheromone detection in Drosophila. Nature. 2007;450(7167):289.

  77. 77.

    Jin X, Ha TS, Smith DP. SNMP is a signaling component required for pheromone sensitivity in Drosophila. Proc Natl Acad Sci. 2008;105(31):10996–1001.

  78. 78.

    Vogt RG, Miller NE, Litvack R, Fandino RA, Sparks J, Staples J, Friedman R, Dickens JC. The insect SNMP gene family. Insect Biochem Mol Biol. 2009;39(7):448–56.

  79. 79.

    Zhang J, Liu Y, Walker WB, Dong SL, Wang GR. Identification and localization of two sensory neuron membrane proteins from Spodoptera litura (Lepidoptera: Noctuidae). Insect Sci. 2015;22(3):399–408.

  80. 80.

    Liu S, Zhang YR, Zhou WW, Liang QM, Yuan X, Cheng J, Zhu ZR, Gong ZJ. Identification and characterization of two sensory neuron membrane proteins from Cnaphalocrocis medinalis (Lepidoptera: Pyralidae). Arch Insect Biochem Physiol. 2013;82(1):29–42.

  81. 81.

    Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644.

  82. 82.

    Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22(13):1658–9.

  83. 83.

    Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25(17):3389–402.

  84. 84.

    Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.

  85. 85.

    Nastou KC, Tsaousis GN, Papandreou NC, Hamodrakas SJ. MBPpred: proteome-wide detection of membrane lipid-binding proteins using profile hidden Markov models. Biochimica et Biophysica Acta (BBA)-proteins and. Proteomics. 2016;1864(7):747–54.

  86. 86.

    Thomas Nordahl P. S?Ren B, Gunnar VH, Henrik N. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8(10):785–6.

  87. 87.

    Letunic I, Bork P. 20 years of the SMART protein domain annotation resource. Nucleic Acids Res. 2017;46(D1):D493–6.

  88. 88.

    Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–4.

  89. 89.

    Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3):R25.

  90. 90.

    Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12(1):323.

  91. 91.

    Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):1–12.

  92. 92.

    Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods. 2001;25(4):402–8.

Download references

Acknowledgements

We thank Xue Zhang for assisting us in collecting specimens and Yingchun Lu for assisting us in plotting data.

Funding

This work was supported by the National Natural Science Foundation of China (Grant No. 31872217, 31372192) to Yuan Huang.

Author information

YH and HY conceived and designed the experiments. HY conducted the experiments and wrote the paper. HHC, LNZ, CY analyzed the data. YH, HHC and HY discussed the results. YH identified insects. All authors read and approved the final manuscript.

Correspondence to Yuan Huang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1. The summary of the Illumina sequencing data. Note: FA: female antennae; MA: male antennae; FH: female head (antennae were cut off); MH: male head (antennae were cut off); FL: female leg; ML: male leg; FW: female wing; MW: male wing; FT: female abdomen-thorax; MT: male abdomen-thorax. Representation of three biological repeats with Arabic numbers 1, 2 and 3. Table S2. Length distribution and quality metrics of C. nigricornis transcripts and unigenes. Table S3. BUSCO analyzed the assembly completeness of C. nigricornis. Table S4. The summary of functional annotation of C. nigricornis transcriptomes. Table S5. Conserved domains of odorant binding proteins (OBPs) in C. nigricornis. Table S6. Conserved domains of chemosensory proteins (CSPs) in C. nigricornis. Table S7. Primers used for qRT-PCR.

Additional file 2: Figure S1. Insect species distribution of C. nigricornis unigenes’ best-hit annotation term in NR database. Figure S2. Gene ontology (GO) classifications of C. nigricornis unigenes. Figure S3. Alignments of the C. nigricornis odorant binding proteins (OBPs). Boxes show the six conserved cysteine residues. Figure S4. Alignments of the C. nigricornis chemosensory proteins (CSPs). Boxes show the four conserved cysteine residues. Figure S5. Phylogenetic tree of odorant-binding proteins (OBPs) from C. nigricornis and other insects. C. kiangsu (Ckia), L. migratoria (Lmig), O. asiaticus (Oasi), O. infernalis (Oinf), S. gregaria (Sgre), A. glycines (Agly), D. ponderosae (Dpon) and H. armigera (Harm). The OBPs of C. nigricornis are represented by red font. Figure S6. Phylogenetic tree of chemosensory-binding proteins (CSPs) from C. nigricornis and other insects. Table S2. L. migratoria (Lmig), O. asiaticus (Oasi), O. infernalis (Oinf), A. gambiae (Agam), D. ponderosae (Dpon) and H. armigera (Harm). The CSPs of C. nigricornis are represented by red font. Figure S7. Phylogenetic tree of odorant receptors (ORs) from C. nigricornis and other insects. Table S3. L. migratoria (Lmig) and A. lineolatus (Alin). The ORs of C. nigricornis are represented by red font. Figure S8. Phylogenetic tree of ionotropic receptors (IRs) from C. nigricornis and other insects. L. migratoria (Lmig), O. asiaticus (Oasi), A. lineolatus (Alin) and D. melanogaster (Dmel). The IRs of C. nigricornis are represented by red font. Figure S9. Phylogenetic tree of sensory neuron membrane proteins (SNMPs) from C. nigricornis and other insects. A. lineolatus (Alin), A. aegypti (Aaeg), A. mellifera (Amel), B. mori (Bmor), D. melanogaster (Dmel), O. asiaticus (Oasi), S. gregaria (Sgre) and T. castaneum (Tcas). The SNMPs of C. nigricornis are represented by red font.

Additional file 3: Table S1. Amino acid sequences of 121 OBPs of Ceracris nigricornis and other insect species used to construct phylogenetic tree. Table S2. Amino acid sequences of 87 CSPs of Ceracris nigricornis and other insect species used to construct phylogenetic tree. Table S3. Amino acid sequences of 293 ORs of Ceracris nigricornis and other insect species used to construct phylogenetic tree. Table S4. Amino acid sequences of 115 IRs of Ceracris nigricornis and other insect species used to construct phylogenetic tree. Table S5. Amino acid sequences of 24 SNMPs of Ceracris nigricornis and other insect species used to construct phylogenetic tree.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Yuan, H., Chang, H., Zhao, L. et al. Sex- and tissue-specific transcriptome analyses and expression profiling of olfactory-related genes in Ceracris nigricornis Walker (Orthoptera: Acrididae). BMC Genomics 20, 808 (2019) doi:10.1186/s12864-019-6208-x

Download citation

Keywords

  • Ceracris nigricornis
  • Transcriptome
  • Expression profiles analysis
  • Odorant-binding protein
  • Chemosensory-binding protein
  • Odorant receptor
  • Ionotropic receptor
  • Sensory neuron membrane protein