Next Article in Journal
Dominance of Fructose-Associated Fructobacillus in the Gut Microbiome of Bumblebees (Bombus terrestris) Inhabiting Natural Forest Meadows
Next Article in Special Issue
Comparison of the Ecological Traits and Boring Densities of Aromia bungii (Faldermann, 1835) (Coleoptera: Cerambycidae) in Two Host Tree Species
Previous Article in Journal
Evaluation of Optimal Reference Genes for qRT-PCR Analysis in Hyphantria cunea (Drury)
Previous Article in Special Issue
Characteristics of Trees Infested by the Invasive Primary Wood-Borer Aromia bungii (Coleoptera: Cerambycidae)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification and Expression Profile of Chemosensory Receptor Genes in Aromia bungii (Faldermann) Antennal Transcriptome

1
Anhui Provincial Key Laboratory of Microbial Control, Engineering Research Center of Fungal Biotechnology, Ministry of Education, School of Forestry & Landscape Architecture, Anhui Agricultural University, Hefei 230036, China
2
Structural Cellular Biology Unit, Okinawa Institute of Science and Technology Graduate University, 1919-1 Tancha, Onna-Son 904-0495, Japan
3
Sterling International, Inc., Spokane, WA 99216, USA
4
Graduate School of Bio-Applications and Systems Engineering, Tokyo University of Agriculture and Technology, Koganei 2-24-16, Tokyo 184-8588, Japan
*
Authors to whom correspondence should be addressed.
Submission received: 27 November 2021 / Revised: 9 January 2022 / Accepted: 11 January 2022 / Published: 14 January 2022
(This article belongs to the Collection Biology and Control of the Invasive Wood-Boring Beetles)

Abstract

:

Simple Summary

There are many chemosensory receptor genes involved in insect chemodetection, including odorant receptors (ORs), gustatory receptors (GRs) and ionotropic receptors (IRs). In contrast to the well-studied Lepidoptera chemosensory receptor genes, the molecular mechanisms of olfactory sensing in Coleoptera are much less understood. The olfactory system plays a crucial role in insect survival. Understanding the olfactory mechanism of insects in depth might provide theoretical guidance for the development of effective pest control measures. The red-necked longicorn beetle, Aromia bungii (Faldermann) (Coleoptera: Cerambycidae), is a wood-boring pest. In order to increase our understanding of the chemosensory receptor genes of the beetle, we first analyzed the transcriptome data of adult A. bungii antennae using bioinformatics, followed by the screening and identification of chemosensory receptor genes. Then, the expression of the chemosensory receptor genes of both male and female adults was examined using qRT-PCR. These findings will provide valuable information for the analysis of the role of chemosensory receptor genes in A. bungii.

Abstract

The red-necked longicorn beetle, Aromia bungii (Faldermann) (Coleoptera: Cerambycidae), is a major destructive, wood-boring pest, which is widespread throughout the world. The sex pheromone of A. bungii was reported earlier; however, the chemosensory mechanism of the beetle remains almost unknown. In this study, 45 AbunORs, 6 AbunGRs and 2 AbunIRs were identified among 42,197 unigenes derived from the antennal transcriptome bioinformatic analysis of A. bungii adults. The sequence of putative Orco (AbunOR25) found in this study is highly conserved with the known Orcos from other Coleoptera species, and these Orco genes might be potentially used as target genes for the future development of novel and effective control strategies. Tissue expression analysis showed that 29 AbunOR genes were highly expressed in antennae, especially in the antennae of females, which was consistent with the idea that females might express more pheromone receptors for sensing pheromones, especially the sex pheromones produced by males. AbunOR5, 29, 31 and 37 were clustered with the pheromone receptors of the cerambycid Megacyllene caryae, suggesting that they might be putative pheromone receptors of A. bungii. All six AbunGRs were highly expressed in the mouthparts, indicating that these GRs may be involved in the taste perception process. Both AbunIRs were shown to be female-mouthparts-biased, suggesting that they might also be related to the tasting processes. Our study provides some basic information towards a deeper understanding of the chemosensing mechanism of A. bungii at a molecular level.

1. Introduction

Insects have evolved a highly specialized and sensitive chemosensory system that can accurately identify some minor environmental changes and specific odorant and tastant materials that are complex in nature, and this ability is primarily reliant on a number of sensory (taste and smell) neurons distributed in their epidermis [1,2]. The chemosensory receptors that recognize chemical signals in insects are mainly distributed on antennae and taste-related organs [3]. Chemosensory recognition in insects is a complex process involving multiple chemosensory-related genes, including odorant-binding proteins (OBPs), odorant-degrading enzymes (ODEs), odorant receptors (ORs), ionotropic receptors (IRs), gustatory receptors (GRs) and sensory neuron membrane proteins (SNMPs) [4,5,6,7,8,9]. The general olfactory recognition process is that external odor molecules enter lymphatic fluid through pores in the cuticular surface of olfactory sensilla, combine with OBPs in the lymphatic fluid to form an OBP–odorant complex and arrive on the dendritic membrane of olfactory receptor neurons carried by OBPs [10]. Molecules then bind to the ORs or IRs on the dendritic membrane that convert the chemical signals into electrical signals. The electrical signals are transmitted to the central nervous system (CNS) through the axon on the other pole of the olfactory receptor neuron and then guide the insect to make the related physiological response. After the completion of the signal transmission, the excess odor molecules are degraded by the ODEs to avoid the damage of the olfactory receptor neurons due to continuous stimulation [3]. Gustatory sensilla have a similar structure, with only a single pore at the top of the sensory hair [11,12].
Olfactory-related receptor genes play a vital role in insects for the recognition of odorant changes in the external environment [3]. In 1999, three laboratories identified Drosophila melanogaster’s olfactory receptor genes almost simultaneously, which began a new chapter in the exploration of insect olfactory mechanisms [13]. This discovery also lays a foundation for the development of olfactory-based insect control technology via the identification of specific molecular targets and olfactory genes. ORs, GRs and IRs are three families of insect chemosensory receptors, all of which were first studied in D. melanogaster [14,15,16]. ORs function as a heterotetrameric receptor complex they form together with a conserved olfactory receptor coreceptor (Orco) [17] in the membrane of olfactory sensory neurones. Orco exists as a single and highly conserved orthologue in each species, and it is necessary for the function of the receptor complex [18]. To date, this OR–Orco system appears to only exist in insects, and the absence of Orco has been documented only from the earliest apterygotes [19]. Apart from Orco, ORs can be divided into two categories: odorant receptors (ORs) and pheromone receptors (PRs) [20]. Research on PR-related genes in Coleoptera has been limited.
Studies of olfactory mechanisms in insects are largely based on model insects (D. melanogaster, Anopheles gambiae, Bombyx mori, Tribolium castaneum, Schistocerca gregaria, Bemisia tabaci, Aedes aegypti, Apis mellifera and Hermetia illucens, etc.) [11,13,21,22,23,24,25,26,27]. Chemosensory genes can be targeted to develop environmentally friendly pest management strategies [28,29]. As the first coleopteran species with a sequenced genome, T. castaneum [30] helped further the study of chemosensory proteins in Coleoptera, such as Capnodis tenebrionis, Holotrichia parallela, Agrilus mali and Lissorhoptrus oryzophilus [31,32,33,34]. The maturity of the next-generation sequencing technology provides the possibility for the study of insect antennae transcriptomes. At present, chemosensory receptor genes from many coleopterans have been identified and analyzed by antennal transcriptomes. According to previous studies based on antennal transcriptomes, a total of 43 ORs, 6 GRs and 7 IRs were identified in the European Spruce bark Beetle, Ips typographus, 49 ORs, 2 GRs and 15 IRs in Dendroctonus ponderosae [35], 22 ORs, 4 GRs and 3 IRs from Dendroctonus valens [36], 43 ORs, 2 GRs and 5 IRs in Anomala corpulenta [37], 43 ORs, 10 GRs and 9 IRs in Colaphellus bowringi [38], 20 ORs and 6 IRs in Tenebrio molitor [39], and 63 ORs, 7 GRs and 28 IRs in Rhynchophorus palmarum [40]. Recently, several chemosensory genes from Cerambycidae were also identified, including Megacyllene caryae (57 ORs) [41], Anoplophora chinensis (53 ORs, 17 GRs and 4 IRs) [42], Apriona germari (42 ORs and 3 IRs) [43], Anoplophora glabripennis (37 ORs, 11 GRs and 7 IRs) [44] and Monochamus alternatus (9 ORs, 1 GRs and 7 IRs) [45]. The chemosensory perception mechanism in Aromia bungii is currently unknown.
The red-necked longicorn beetle, A. bungii, is an important wood-boring pest of peach, apricot, plum and other fruit trees [46]. A. bungii is difficult to control because their larvae live in a protected habitat beneath the bark of trees [47]. Larvae of A. bungii damage tree branches, stem phloem and xylem, cut off the transport tissues of trees and accumulate insect feces and sawdust, resulting in peach branch dry gum and weakening of the tree, which is destructive to fruit trees, forest production and urban landscaping [48]. The male-produced aggregation pheromone of A. bungii was identified as (E)-2-cis-6,7-epoxynonenal [46], while (R)-(+)-citronellal was reported as the female-produced sex pheromone component [49]. Both pheromones might have great potential as an attractant for monitoring and controlling this highly destructive cerambycid beetle.
In the current study, we conducted an antennal transcriptome analysis of A. bungii adults and identified a total of 53 putative chemosensory receptor genes, including 45 ORs, 6 GRs and 2 IRs. Moreover, the expression analysis of candidate chemosensory receptor genes in different tissues of both sexes, including antennae, mouthparts (maxillary palps and labial palps) and abdominal tips, was validated via quantitative Real-Time PCR (qRT-PCR) in order to discover the chemosensory genes that may play a key role in the life cycle of A. bungii. Our findings on these chemosensory receptors might lead to a new perspective for controlling this economically important longhorn beetle via the identification of specific molecular targets and chemosensory genes.

2. Materials and Methods

2.1. Insects and Tissue Collections

Newly emerging adults of A. bungii were collected from the campus of Anhui Agricultural University in Hefei, Anhui Province, in June 2019. The collected adults were maintained in the laboratory on sugarcane stems at 28–30 °C with a photoperiod of 14 h:10 h (light:dark) [50]. For RNA isolation and transcriptome analyses, antennae, mouthparts (maxillary palps and lower labial palps) and abdominal tips of healthy male and female adults were collected and placed in liquid nitrogen for quick freezing and then stored at −80 °C for subsequent experiments.

2.2. RNA Isolation and cDNA Library Construction

Total RNAs were separately extracted from adult antennae (both sexes) using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer’s instructions [43] and then treated with RNase-free DNase I (TaKaRa, Dalian, Liaoning, China). A NanoDrop ND-2000 Spectrophotometer (Nanodrop Technologies, Wilmington, DE, USA) was used to determine RNA concentration, and UV absorption values were recorded at 230/260 nm and 260/280 nm to test the purity of RNA products. RNA integrity was monitored on 1% agarose gel electrophoresis. The qualified total RNA of the antennae was frozen in liquid nitrogen and stored at −80 °C before being processed or placed in dry ice for transport to Novogene Co., Ltd. (Beijing, China) for transcriptome sequencing [50].

2.3. Transcriptome Sequencing and Data Assembly

Purified RNAs were prepared for the cDNA library using the TruSeq RNA Sample Preparation Kit (Illumina, San Diego, CA, USA) following the manufacturer’s instructions. After sequencing, paired-end reads were generated and firstly processed through Casava software (v1.7). Clean reads were obtained after the quality control of raw data and then spliced using Trinity software (v2.4.0) to obtain reference sequences for subsequent analysis. The longest transcript sequence (unigene) was annotated in NCBI non-redundant protein sequences (NCBI-nr), NCBI nucleotide sequences (NCBI-nt), Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG-Ontology, or KO), Protein family database (Pfam), EuKaryotic Ortholog Groups/Clusters of Orthologous Groups (KOG/COG) and Swiss-Prot databases by BLAST alignment with a cut-off E-value of 10−5. Preliminary candidate chemosensory-related gene nucleic acid sequences were obtained by screening the transcriptome data of adult antennae (males and females). Then, a local nucleic acid database was established in BioEdit software (v7.0.9.0), and amino acid sequences of chemosensory receptors related to Coleoptera were downloaded from the NCBI database as query sequences. The TBLASTN program was used to perform a local BLAST search (E-value of 10−5) on transcriptome data of A. bungii to receive the predicted chemosensory gene sequences [36]. BlastX (NCBI database) was used to compare these predicted sequences to screen out the candidate chemosensory genes, and relevant parameters were recorded, including the length of amino acid sequences, name of the species with the highest homology, name of the gene, entry number, E-values and gene similarity. We used the online prediction open reading frame (ORF Finder) website (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/orffinder/, accessed on 29 July 2021) to obtain the ORFs. The transmembrane prediction of receptor proteins was based on the TMHMM online website (https://services.healthtech.dtu.dk/service.php?TMHMM-2.0/, accessed on 12 August 2021) [51].

2.4. Phylogenetic Analysis

The evolutionary tree of candidate ORs was constructed from the protein sequences of A. germari [43], A. chinensis [42], A. glabripennis [44], Agrilus planipennis [52], T. castaneum [53], D. ponderosae [35], Ips typographus [35], A. corpulenta [37], M. caryae [41], M. alternatus [45], Phyllotreta striolata [54], D. valens [36], C. bowringi [38] and T. molitor [39]. The evolutionary tree of candidate GRs was built with the aligned protein sequences from A. chinensis, Ips typographus, A. glabripennis, D. melanogaster [5], D. ponderosae, B. mori, T. castaneum and P. striolata. The evolutionary tree of candidate IRs was established with the protein sequences from A. chinensis, A. glabripennis, D. ponderosae, D. melanogaster, P. striolata, A. corpulenta, Monochamus alternatus, M. caryae and T. molitor. ClustalX2.0 was used for a complete comparison of all sequences. The neighbor-joining (N-J) method was used to input 1000 replicates, and PHB results from ClustalX were placed in MEGA software (v5.0) to construct the evolutionary tree [55]. Finally, the phylogenetic tree was modified on the online website iTOL (https://itol.embl.de/, accessed on 14 September 2021) and annotated with PS technology.

2.5. Quantitative Real-Time PCR Analysis of Gene Expression

Total RNA was isolated from antennae, mouthparts (maxillary palp and labial palp) and abdominal tips (both sexes) of A. bungii adults using TRIzol reagent [42]. Quantitative real-time PCR (qRT-PCR) was used to evaluate the expression of candidate ORs, IRs and GRs in different tissues of each sex. Reverse RNA into cDNA was synthesized using a Prime Script RT Reagent Kit (Takara-bio, Shiga, Japan) with a gDNA Eraser (Perfect Real Time) according to the manufacturer’s instructions [43]. The cDNA was used as the template, and β-actin was used as the reference gene for qRT-PCR verification with different primers (Table S1). Three biological and three technical replicates were performed on each sample.

2.6. Statistical Analysis

The relative expression level of mRNA of each gene was normalized to those of the actin gene and calculated using the Q-gene method [56]. The significant difference among the experiment groups was analyzed using one-way ANOVA, followed by Duncan’s new multiple-range test (p = 0.05) [42], and values were presented as the mean ± SD. Graphics plot mapping was carried out using Graphpad Prism software v5.0 (GraphPad Inc., San Diego, CA, USA).

3. Results

3.1. Transcriptome Sequence and Homologous Assembly

The transcriptome information of the longicorn beetle, A. bungii, was characterized by constructing a cDNA library prepared from purified mRNA isolated from the adults’ antennae. Using Illumina sequencing, we obtained 45,642,924 raw reads and 43,302,906 (6.5 GB) clean reads after Trinity assembly. The clean reads were assembled into 79,280 transcripts. The longest cluster sequence was obtained via Corset hierarchical clustering as unigenes, and 42,197 unigenes were obtained (Table 1). The mean lengths of transcripts and unigenes were 783 bp and 1215 bp, respectively. Among the 42,197 unigenes, 31,489 were larger than 500 bp in length, accounting for 74.62%. The length distribution of transcripts and unigenes showed overall declines in the number of transcripts with increasing length (Supplementary Figure S1). However, the number of unigenes increased at first and then declined in the range between 1001 and 2000 bp. More than one-fifth of short reads (<301 bp) were assembled into a transcript, while reads over 500 bp were assembled into as many unigenes as there were transcripts, implying that longer reads (>500 bp) are more likely to be assembled into unigenes.

3.2. Functional Annotation and Expression Level

A total of 42,197 unigenes were functionally annotated in NCBI-nr, NCBI-nt, GO, KO, Pfam, KOG, Swiss-Prot and other databases (Table S2). Among them, 25,407 (60.21%) unigenes had homologous sequences in NCBI-Nr, 7292 (17.28%) unigenes in NCBI-Nt totaled 19,959 (47.29%) in GO, 12,154 (28.8%) in KO, 19,426 (46.03%) in Pfam, 11,787 (27.93%) in KOG and 19,226 (45.56%) in Swiss-Prot. Among the total of unigenes identified in the antennal transcriptome, 2696 genes were expressed at a very low level (FPKM = 0), and 8761 unigenes were highly expressed (FPKM > 10) (Table S3).
Among 25,407 (60.21%) unigenes that were compared to proteins in the NCBI Nr protein database, 84.5% of the annotated genes had more than 60% similarity with known proteins (Figure 1a). According to the E-value distribution diagram (Figure 1b), 59.6% of the annotated genes showed strong homology (E-value < 1 × 10−45), while 11.6% showed very low homology (1 × 10−15 < E-value < 1 × 10−5). The species distribution map showed that 60.4% of the transcripts had the highest homology with T. castaneum, followed by Dendroctonus ponderossae (17.7%), Acyrthosiphon pisum (1.2%), Leptinotarsa decemlineata (1.2%) and Camponotus floridanus (1.2%) (Figure 1c).
GO analysis was used to classify unigenes into different functional categories. There were 19,959 (47.29%) unigenes that were successfully classified in the biological processes, cellular compartment and molecular function. The most represented classes were assigned to different biological processes (54,453 unigenes). The rest were classified to the cellular components (36,327 unigenes) and the molecular functions (25,077 unigenes). These classifications were categorized into different subclasses. Among the categories of biological processes, the subclass of cellular processes (11,532 unigenes), metabolic processes (10,592 unigenes) and single-organism processes (9304 unigenes) were the most annotated. Among the categories of cell components, the most annotated were cell (7002 unigenes), cell part (7002 unigenes) and organelle (5002 unigenes). Among the categories of molecular functions, the most annotated ones were binding (11,906 unigenes) and catalytic activity (8608 unigenes) (Figure 2).
In the antennal transcriptome, 12,154 unigenes were KO annotated, and could be divided into five categories according to the KEGG metabolic pathway: cellular processes (Figure 3a), environmental information processing (Figure 3b), genetic information processing (Figure 3c), metabolism–metabolism (Figure 3d) and organizational systems (Figure 3e). Signal Transduction (1527 unigenes) was the most annotated gene pathway in the metabolic branch of environmental information processing, and it was also the one with the most annotated genes among all the participating metabolic pathways.
After comparing with KOG database, 11,787 unigenes (27.93%) were annotated. The abscess axis A–Z represents 26 types: the largest proportion of genes was involved in general function prediction alone (1929 unigenes; 16.37%); the number of unigenes in the KOG database was the largest, followed by signal transduction mechanisms (1424 unigenes; 12.08%) and posttranslational modification/protein conversion/chaperones (1265; 10.73%) (Figure 4).

3.3. Identification of Putative Odorant Receptors

A total of 45 putative AbunORs were identified based on the comparative analysis of the antennal transcriptome of A. bungii using BlastX database, of which, six AbunORs (AbunOR5, AbunOR9, AbunOR17, AbunOR25, AbunOR31 and AbunOR38) had full open reading frames (ORFs). These six AbunORs displayed seven transmembrane domains, a classical feature of ORs from insects (Table 2).
Based on the OR phylogenetic tree analysis that was constructed to evaluate the relationships of AbunORs with other reported coleopteran ORs, most AbunORs were clustered together with OR genes of A. chinensis [42], A. glabripennis [44], M. alternatus [45] and M. caryae [41] with a high sequence similarity. AbunOR25 was distributed in the Orco gene family and had a high sequence similarity with AchiOR1, AgerOR25, AcorOrco [37], MaltOR1 [45], TmolOrco [39], TcasOR1 [30] and McarOR1 [41], which was consistent with the high conservation of Orco gene among insects. As reported in the literature, OR genes in Coleoptera were divided into seven functional subgroups (clade 1–7) [53]. A total of 44 putative AbunORs except AbunOR25 were classified into five subgroups (clade 1–3, 5 and 7), among which, eight AbunORs (AbunOR6, 13, 20, 27, 28, 35, 37 and 38) were clustered in clade 1, nine AbunORs (AbunOR1, 2, 5, 8, 29, 31, 41, 43 and 44) were assigned to clade 2, twelve AbunORs (AbunOR4, 9, 11, 12, 14, 15, 17, 19, 21, 23, 32 and 45) were placed in clade 3, three AbunORs (AbunOR26, 33 and 39) were categorized into clade 5, and the remaining twelve AbunORs (AbunOR3, 7, 10, 16, 18, 22, 24, 30, 34, 36, 40 and 42) were included in clade 7. In addition, three AbunORa (AbunOR5, 29, 31 and 37) were clustered with the reported pheromone receptors (PRs) (McarOR3, 5 and 20) [41], suggesting that these three AbunORs might be the pheromone receptors for A. bungii (labeled as PR in Figure 5).

3.4. Tissue- and Sex-Specific Expression Analysis of Putative Odorant Receptors

The expression of OR genes in different tissues was studied via RT-qPCR. Different tissues of both sexes including antennae, mouthparts (maxillary palps, labial palps) and terminal (abdominal) tips were studied and analyzed (Figure 6). In antennae, 29 AbunORs (AbunOR1–2, 4–10, 12, 14, 16–17, 19, 21–23, 25, 29, 33–36, 39 and 41–45) showed female-biased expression. Nine AbunORs (AbunOR11, 13, 15, 18, 20, 26–28 and 31) were expressed similarly in the antennae of both sexes. All 45 AbunORs genes were expressed in the mouthparts, and 13 AbunORs (AbunOR4, 6, 12, 16, 18, 23–24, 26, 28, 31, 33, 37 and 40) were expressed significantly higher in the mouthparts of females than in those of males, while only two AbunORs (AbunOR5 and 38) were expressed significantly higher in the mouthparts of males than in those of females. All 45 AbunOR genes were expressed in the abdominal tips at very low levels and much lower than those in the antennae and mouthparts (maxillary palps, labial palps) (Figure 6 and Figure S2).

3.5. Identification of Putative Gustatory Receptors

Six AbunGRs (AbunGR1-6) were identified by bioinformatics analysis of the antennaal transcriptome. Five AbunGRs (AbunGR1–4, AbunGR6) had a complete open reading frame (ORF). Transmembrane domain prediction results showed that three AbunGRs (AbunGR4–6) had the predicted transmembrane helix structure (Table 3). In the GR evolutionary tree, the six AbunGRs were divided into four taste receptor gene families, namely sugar, fructose, bitter and carbon dioxide receptor families. Only AbunGR6 was clustered with the reported carbon dioxide receptors DmelGR21a and DmelGR63a [57] (Figure 7).

3.6. Tissue- and Sex-Specific Expression Analysis of Putative Gustatory Receptors

All six AbunGRs (AbunGR1–6) had a significantly high expression in the beetle mouthparts (Figure 8). Four of the six GR genes (AbunGR1–2, 4–5) were significantly female-biased in mouthparts. AbunGR2′s expression in the antennae of females was higher than in those of males. AbunGR6 showed similar expression levels between females and males in all the analyzed tissues. All six GRs genes (AbunGR1–6) were generally relatively female-biased in the antennae and mouthparts, while relative mRNA expression in abdominal tips was significantly low (Supplementary Figure S3).

3.7. Identification and Expression Analysis of Putative Ionotropic Receptors

Two AbunIRs were identified by analyzing the antennal transcriptome data (Table 4), and only AbunIR2 had a complete ORF with the most conserved three transmembrane domains. After BlastX homology comparison, AbunIR1 was compared to the IR of A. germari (52%), while AbunIR2 was compared to the IR gene of P. striolata (44%). According to the phylogenetic analysis of IRs from nine species (Figure 9a), the two AbunIRs were both mainly clustered with the IRs of coleopteran species, and the sequence similarity between AbunIR2, PstrIR47 and AcorIR75q was greater than 90% (Supplementary Figure S4). RT-qPCR results showed that AbunIR1 and AbunIR2 were specifically expressed in the mouthparts (maxillary palps and labial palps) of females (Figure 9b). Both AbunIRs (AbunIR1–2) were significantly overexpressed in the mouthparts of females compared to other tissues of both female and male A. bungii adults (Supplementary Figure S5).

4. Discussion

As the largest order of Insecta, Coleoptera is also one of the most important pest orders for forestry and agriculture. Among the groups of Coleoptera, many cerambycids (both native or invasive species) have become serious forest (or tree) pest insects, causing significant ecological and economical losses throughout the world [58]. Over the past decade, significant research progresses have been made regarding the chemical ecology of longhorn beetles, especially the identification of aggregation-sex pheromones (attractants) for some economically important species. However, compared to the well-studied lepidopterans, research on olfactory mechanisms at the molecular/gene levels in Cerambycidae still remains limited [59,60]. In the current study, we identified chemosensory receptor genes and studied their expression profiles based on the antennal transcriptome data of a highly damaging pest of Prunus fruit trees, the red-necked longhorn beetle, A. bungii.
A total of 42,197 unigenes were identified from the antennal transcriptome data, of which 74.62% were more than 500 bp in length. The NCBI-nr database had the largest number of unigenes (with homologous sequences) with a total of 25,407 unigene. A. bungii showed the highest homology with T. castaneum (60.4%), followed by D. ponderossae (17.7%). The number of unigenes successfully annotated in the GO database was lower than that in the NCBI-nr database, with a total of 19,959 unigenes. Among them, binding and catalytic activities were the most annotated, which are similar to the functions of olfactory-related genes. In addition, 12,154 unigenes were annotated in the KO database, and KEGG metabolic pathway classification showed that the subcategory of “signal transduction” is the most annotated gene pathway in the metabolic category of “environmental information processing” (Figure 2). These results can be related to the function of olfactory-related genes binding to odorant molecules to complete signal transduction [11,12,61]. Based on the annotated KOG database, the “general function prediction only” category was the largest, followed by “signal transduction mechanisms” (Figure 4). This result strongly suggested that these unigenes might be paralogs, or evolved from the same species of genes, and had molecular functions similar to those of the olfactory-related genes.
A total of 45 odorant receptor (ORs) genes, 6 gustatory receptors (GRs) genes and 2 ionotropic receptors (IRs) were obtained from transcriptome analysis of A. bungii antennae. Compared to the previous study on other coleopterans, the number of odorant receptor (OR) genes identified from A. bungii (45) was similar to those from Apriona germari (42) [43], from Holotrichia oblita (44) [62], from Anomala corpulenta (43) [37] and from I. typographus (43) [35]. As can be seen from the phylogenetic tree, AbunOR25 is distributed in the Orco gene family and has high sequence similarities with AchiOR1 [42], AgerOR25 [43], AcorOrco [37], MaltOR1 [45], TmolOrco [39], TcasOR1 [30], McarOR1 [41] and other genes from cerambycidae, which is consistent with the feature of high Orco conservatism among insects (Supplementary Figures S6 and S7). Orco exists as a single and highly conserved orthologue in each species, and it is necessary for the function of the receptor complex [63]. Since the sequence of Orco genes in Cerambycidae is highly conserved, the Orco genes of these species can be used as potential interference target genes for integrated pest management. So far, the odorant receptors in Coleoptera were divided into seven functional subgroups (clade from 1–7). Forty-four of the forty-five AbunORs (except AbunOR25) were distributed in five subgroups (clade 1–3, 5 and 7a). Previously, three pheromone receptors of M. caryae were identified: McarOR3 being sensitive to sex pheromone component (S)-2-methyl-1-butanol, McarOR5 being sensitive to 2-phenylethanol and McarOR20 being sensitive to (R)-3-hydroxyhexan-2-one and (2S, 3R)-2, 3-hexanediol [41,62]. Interestingly, AbunOR5 was clustered with McarOR3 in the phylogenetic tree, suggesting that AbunOR5 might be also sensitive to (S)-2-methyl-1-butanol. AbunOR29/31 that clustered with McarOR5 might respond to the sex pheromone component 2-phenylethanol, whereas AbunOR37 might respond to (R)-3-hydroxyhexan-2-one and (2S,3R)-2,3-hexanediol since it was clustered with McarOR20. In other words, these four AbunORs (5, 29, 31 and 37) were clustered with the three pheromone receptors of M. caryae and showed relatively similar amino acid sequences with them (Supplementary Figures S8–S10). These results suggested that AbunORs (5, 29, 31, 37) might be the putative pheromone receptors. These potential pheromone receptors were not all clustered in the same clade, strengthening the idea that pheromone receptors (PRs) from different beetle species do not cluster in specific clades, unlike PRs in Lepidoptera, where the majority of the characterized PRs are found in the so-called classical PR clade [64]. Instead, beetle PRs are scattered in the OR phylogeny, and OR clades include receptors detecting compounds from various ecological sources [65]. The divergence of this evolutionary branch of olfactory genes within species contributes the receptors to detect compounds derived from different ecological sources.
It is commonly known that odorant receptors are mainly expressed in the antennae. Since the pheromone of A. bungii is mainly produced by males, it is speculated and shown that the expression levels of OR genes in the females’ antennae are higher than that in males’ antennae. In this study, among the 45 putative AbunORs, 29 were highly expressed in the female antennae. AbunOR5, AbunOR29, AbunOR31 and AbunOR37, which were clustered together with three possible pheromone receptors (PRs) of M. caryae in the phylogenetic tree, were also significantly overexpressed in females. This is consistent with the idea that females may express more pheromone receptors for the sex pheromones released by males. In addition, all 45 AbunORs were expressed in the mouthparts, and most of them were expressed significantly more in the female mouthparts than in the male mouthparts, suggesting that they may be involved in the host selection and oviposition by females. The results of qRT-PCR analysis of T. castaneum showed that, with one exception, all TcasORs were expressed in the antennae but not in any taste organs [53]. Based on the tissue- and sex-specific expressions analysis of the A. germari, 27 of the 40 AgerORs expressed in antennae were female-biased, and only 2 AgerORs were expressed at the same level in both females and males. Three AgerORs were expressed as female-biased in maxillary palps. In addition, among the 17 AgerORs expressed in labial palps, 3 were highly expressed in the labial palps of females. In addition, 3 AgerORs were highly expressed in the abdominal tips of females, and 13 AgerORs showed high expression in the abdominal tips of males [43]. A total of 45 ORs were identified from A. chinensis antennal transcriptome. Forty-one putative OR genes were significantly expressed in the beetle antennae, of which eight AchiORs were significantly female-biased, while twenty-three ORs were significantly male-biased, and the remaining ten were expressed at similar levels in the antennae of both females and males. Moreover, only one AchiOR was highly expressed in the labial palps of both sexes [42]. These results indicated that the olfactory receptors of A. bungii were mainly overexpressed in female antennae, similar to most of the previous reports on Coleoptera-related olfactory receptor genes [32,34,43,62,66,67].
We identified six AbunGR genes, which were fewer than those identified from A. glabripennis (11) [44] and A. chinensis (17) [42], and more than those from I. typographus (6) and D. ponderosae (2) [35]. AbunGR6 is a potential homologue of carbon dioxide receptors of D. melanogaster, DmelGR21a and DmelGR63a, and might be involved in the recognition and detection of CO2 [68]. The remaining AbunGRs were not distributed in the other four known GR families (sugar, amino acids, salts and bitter compounds), and the similarity of GR gene sequences among different species was low in the phylogenetic tree (Figure 7). The divergence of the GR genes was remarkable, as the similarity between most receptor pairs was only 20% or less (at the amino acid sequence level) [57]. The expression of the T. castaneum chemoreceptor genes was investigated using qRT-PCR. All the predicted TcasGRs were expressed in the mouth parts and mostly in the prolegs of the adults; only seven TcasGRs were not expressed in the prolegs. The results of the qRT-PCR analysis of A. chinensis showed that most AchiGRs were prominently expressed in antennae [42]. In the current study, all six AbunGRs were highly expressed in gustatory organs based on the tissue- and sex-specific expressions analysis, suggesting that these AbunGRs were likely involved in the detection of soluble stimulants and feeding behaviors. Additionally, AbunGR2 was highly expressed in adult antennae, presumably because all olfactory and gustatory genes were derived from antennal transcriptome data rather than the complete genome.
Only two IR genes were identified from the antennal transcriptomic analysis in this study, which were fewer than those from P. striata (49) [54], D. ponderosae (15) [35] and B. longissima (19) [69] but were similar to those of A. germari (3) [43], A. corpulenta (5) [37], A. chinensis (4) [42] and A. glabripennis (4) [44]. IRs are divided into three subfamilies in the Drosophila genome: olfactory IRs, divergent IRs and co-receptor IRs [70]. Most olfactory IRs are specifically expressed in antennae [61], but not in other tissues or at low levels, so they are also called antennal IRs. Divergent IRs are to a large extent species-specific and, unlike olfactory IRs, are almost never expressed in antennae. Some divergent IRs are expressed in taste organs, suggesting a possible involvement in the taste perception [3,71]. As reported, AgerIRs were highly expressed in the abdominal tips and labial palps, while the expression of AchiIR in antennae was relatively high. In this study, the two putative AbunIRs (AbunIR1 and AbunIR2) were highly expressed in the mouthparts (maxillary palps and labial palps) of the females; thus, they might belong to divergent IRs and participate in the gustatory process in the taste organs.
In conclusion, 45 AbunORs, 6 AbunGRs and 2 AbunIRs were identified via antennal transcriptome bioinformatic analysis of A. bungii adults. The new Orco gene (AbunOR25) found in this study and other Orco genes of Cerambycidae species were all highly conservative; thus, we speculate that Orco genes might be used as potential interfering targets for further research and exploration of the development of a novel, viable and effective pest control strategy. Most olfactory receptors in A. bungii were significantly expressed in antennae, especially in the female’s antennae, suggesting that these receptor genes might be heavily involved in female-specific behaviors. Further analysis on their potential functionality showed that AbunOR5, AbunOR29, AbunOR31 and AbunOR37 were clustered with pheromone receptors of M. caryae, suggesting that they might be the pheromone receptors of A. bungii. All six AbunGRs were highly expressed in gustatory organs; thus, they are likely involved in taste perception. The two AbunIRs were both highly expressed in the female’s mouthparts, suggesting that they might also participate in the tasting processes.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/insects13010096/s1, Figure S1: Distribution of Unigenes and Transcripts length interval in the A. bungii transcriptome assembly, Figure S2: Analysis of the relative expression of 45 AbunORs in female and male tissues, Figure S3: Analysis of the relative expression of 6 AbunGRs in female and male tissues, Figure S4: Multiple amino acid sequence alignment of AbunIR2, PstrIR47and AcorIR75q, Figure S5: Analysis of the relative expression of 2 AbunIRs in female and male tissues, Figure S6: Amino acid sequence alignment of four known Orcos AgerOR25, AchiOR1, AglaOrco from Cerambycidae and AbunOR25 from A. bungii, Figure S7: Amino acid sequence alignment of nine known Orcos and a new Orco (AbunOR25) of A. bungii identified in this study, Figure S8: Amino acid sequence alignment of AbunOR3 and McarOR3 (Mcar, M. caryae), Figure S9: Amino acid sequence alignment of McarOR3 (Mcar, M. caryae), AbunOR29 and AbunOR31, Figure S10: Amino acid sequence alignment of McarOR3 (Mcar, M. caryae) and AbunOR37, Table S1: Primers of A. bungii chemosensory receptor genes used for RT-qPCR, Table S2: Summary of the gene annotation success ratio, Table S3: Summary of FPKM value.

Author Contributions

L.Z., D.M. and Z.W.: conceived and designed the experiments; Z.W., J.Q., J.Y., L.Z. and D.M.: searched the literature; Z.W., J.Q. and J.Y.: performed the experiments; Z.W., L.Z., D.M., Q.Z. and E.R.P.: analyzed the data; L.Z., D.M., Q.Z. and E.R.P.: interpreted the data; Z.W., L.Z., D.M., Q.Z. and E.R.P.: drafted the manuscript; Z.W., J.Q., J.Y., E.R.P., Q.Z., L.Z. and D.M.: revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (31870640, 31170616), National Key Research and Development Program (2017YFD0600101), Forest-Scientific Research in the Public Interest (No.201204506), Research Innovation Program for College Graduates of Anhui Agricultural University (Grant No. 2019YSJ-16).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data can be provided on request from the lead author.

Acknowledgments

We would like to thank Peng He, Tianzi Gu and Yanan Zhang for their helpful suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Montagné, N.; de Fouchier, A.; Newcomb, R.D.; Jacquin-Joly, E. Advances in the identification and characterization of olfactory receptors in insects. Prog. Mol. Biol. Transl. Sci. 2015, 130, 55–80. [Google Scholar] [CrossRef]
  2. Schmidt, H.R.; Benton, R. Molecular mechanisms of olfactory detection in insects: Beyond receptors. Open Biol. 2020, 10, 200252. [Google Scholar] [CrossRef]
  3. Wicher, D.; Miazzi, F. Functional properties of insect olfactory receptors: Ionotropic receptors and odorant receptors. Cell Tissue Res. 2021, 383, 7–19. [Google Scholar] [CrossRef]
  4. Vogt, R.G.; Riddiford, L.M. Pheromone binding and inactivation by moth antennae. Nature 1981, 293, 161–163. [Google Scholar] [CrossRef]
  5. Benton, R.; Vannice, K.S.; Gomez-Diaz, C.; Vosshall, L.B. Variant ionotropic glutamate receptors as chemosensory receptors in Drosophila. Cell 2009, 136, 149–162. [Google Scholar] [CrossRef] [Green Version]
  6. Younus, F.; Chertemps, T.; Pearce, S.L.; Pandey, G.; Bozzolan, F.; Coppin, C.W.; Russell, R.J.; Maïbèche-Coisne, M.; Oakeshott, J.G. Identification of candidate odorant degrading gene/enzyme systems in the antennal transcriptome of Drosophila melanogaster. Insect Biochem. Mol. Biol. 2014, 53, 30–43. [Google Scholar] [CrossRef] [Green Version]
  7. Liu, G.; Ma, H.; Xie, H.; Xuan, N.; Guo, X.; Fan, Z.; Rajashekar, B.; Arnaud, P.; Offmann, B.; Picimbon, J.F. Biotype Characterization, Developmental Profiling, Insecticide Response and Binding Property of Bemisia tabaci Chemosensory Proteins: Role of CSP in Insect Defense. PLoS ONE 2016, 11, e0154706. [Google Scholar] [CrossRef] [Green Version]
  8. Chen, Y.D.; Dahanukar, A. Recent advances in the genetic basis of taste detection in Drosophila. CMLS 2020, 77, 1087–1101. [Google Scholar] [CrossRef]
  9. Zhang, H.J.; Xu, W.; Chen, Q.M.; Sun, L.N.; Anderson, A.; Xia, Q.Y.; Papanicolaou, A. A phylogenomics approach to characterizing sensory neuron membrane proteins (SNMPs) in Lepidoptera. Insect Biochem. Mol. Biol. 2020, 118, 103313. [Google Scholar] [CrossRef]
  10. Ache, B.W.; Young, J.M. Olfaction: Diverse species, conserved principles. Neuron 2005, 48, 417–430. [Google Scholar] [CrossRef] [Green Version]
  11. Sánchez-Gracia, A.; Vieira, F.G.; Rozas, J. Molecular evolution of the major chemosensory gene families in insects. Heredity 2009, 103, 208–216. [Google Scholar] [CrossRef]
  12. Leal, W.S. Odorant reception in insects: Roles of receptors, binding proteins, and degrading enzymes. Ann. Rev. Entomol. 2013, 58, 373–391. [Google Scholar] [CrossRef]
  13. Pennisi, E. Fruit fly odor receptors found. Science 1999, 283, 1239. [Google Scholar] [CrossRef]
  14. Gao, Q.; Chess, A. Identification of candidate Drosophila olfactory receptors from genomic DNA sequence. Genomics 1999, 60, 31–39. [Google Scholar] [CrossRef]
  15. Clyne, P.J.; Warr, C.G.; Carlson, J.R. Candidate taste receptors in Drosophila. Science 2000, 287, 1830–1834. [Google Scholar] [CrossRef]
  16. Fleischer, J.; Pregitzer, P.; Breer, H.; Krieger, J. Access to the odor world: Olfactory receptors and their role for signal transduction in insects. Cell. Mol. Life Sci. 2018, 75, 485–508. [Google Scholar] [CrossRef]
  17. Zufall, F.; Domingos, A.I. The structure of Orco and its impact on our understanding of olfaction. J. Gen. Physiol. 2018, 150, 1602–1605. [Google Scholar] [CrossRef] [Green Version]
  18. Martin, F.; Alcorta, E. Regulation of olfactory transduction in the orco channel. Front. Cell. Neurosci. 2011, 5, 21. [Google Scholar] [CrossRef] [Green Version]
  19. Brand, P.; Robertson, H.M.; Lin, W.; Pothula, R.; Klingeman, W.E.; Jurat-Fuentes, J.L.; Johnson, B.R. The origin of the odorant receptor gene family in insects. eLife 2018, 7, e38340. [Google Scholar] [CrossRef]
  20. Zhang, X.; Liu, P.; Qin, Q.; Li, M.; Meng, R.; Zhang, T. Characterizing the Role of Orco Gene in Detecting Aggregation Pheromone and Food Resources in Protaetia brevitarsis Leiws (Coleoptera: Scarabaeidae). Front. Physiol. 2021, 12, 649590. [Google Scholar] [CrossRef]
  21. 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, 159–170. [Google Scholar] [CrossRef] [Green Version]
  22. Zeng, Y.; Yang, Y.T.; Wu, Q.J.; Wang, S.L.; Xie, W.; Zhang, Y.J. Genome-wide analysis of odorant-binding proteins and chemosensory proteins in the sweet potato whitefly, Bemisia tabaci. Insect Sci. 2019, 26, 620–634. [Google Scholar] [CrossRef] [Green Version]
  23. Bohbot, J.; Pitts, R.J.; Kwon, H.W.; Rützler, M.; Robertson, H.M.; Zwiebel, L.J. Molecular characterization of the Aedes aegypti odorant receptor gene family. Insect Mol. Biol. 2007, 16, 525–537. [Google Scholar] [CrossRef] [Green Version]
  24. Robertson, H.M.; Wanner, K.W. The chemoreceptor superfamily in the honey bee, Apis mellifera: Expansion of the odorant, but not gustatory, receptor family. Genome Res. 2006, 16, 1395–1403. [Google Scholar] [CrossRef] [Green Version]
  25. Scieuzo, C.; Nardiello, M.; Farina, D.; Scala, A.; Cammack, J.A.; Tomberlin, J.K.; Vogel, H.; Salvia, R.; Persaud, K.; Falabella, P. Hermetia illucens (L.) (Diptera: Stratiomyidae) Odorant Binding Proteins and Their Interactions with Selected Volatile Organic Compounds: An In Silico Approach. Insects 2021, 12, 814. [Google Scholar] [CrossRef]
  26. Carey, A.F.; Carlson, J.R. Insect olfaction from model systems to disease control. Proc. Natl. Acad. Sci. USA 2011, 108, 12987–12995. [Google Scholar] [CrossRef] [Green Version]
  27. Hansson, B.S.; Stensmyr, M.C. Evolution of insect olfaction. Neuron 2011, 72, 698–711. [Google Scholar] [CrossRef] [Green Version]
  28. Wang, G.; Carey, A.F.; Carlson, J.R.; Zwiebel, L.J. Molecular basis of odor coding in the malaria vector mosquito Anopheles gambiae. Proc. Natl. Acad. Sci. USA 2010, 107, 4418–4423. [Google Scholar] [CrossRef] [Green Version]
  29. Zhou, L.Y.; Li, W.; Liu, H.Y.; Xiang, F.; Kang, Y.K.; Yin, X.; Huang, A.P.; Wang, Y.J. Systemic identification and analyses of genes potentially involved in chemosensory in the devastating tea pest Basilepta melanopus. Comp. Biochem. Physiol. Part D Genom. Proteom. 2019, 31, 100586. [Google Scholar] [CrossRef] [PubMed]
  30. Richards, S.; Gibbs, R.A.; Weinstock, G.M.; Brown, S.J.; Denell, R.; Beeman, R.W.; Gibbs, R.; Beeman, R.W.; Brown, S.J.; Bucher, G.; et al. The genome of the model beetle and pest Tribolium castaneum. Nature 2008, 452, 949–955. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Bari, G.; Scala, A.; Garzone, V.; Salvia, R.; Yalcin, C.; Vernile, P.; Aresta, A.M.; Facini, O.; Baraldi, R.; Bufo, S.A.; et al. Chemical Ecology of Capnodis tenebrionis (L.) (Coleoptera: Buprestidae): Behavioral and Biochemical Strategies for Intraspecific and Host Interactions. Front. Physiol. 2019, 10, 604. [Google Scholar] [CrossRef]
  32. Yi, J.K.; Yang, S.; Wang, S.; Wang, J.; Zhang, X.X.; Liu, Y.; Xi, J.H. Identification of candidate chemosensory receptors in the antennal transcriptome of the large black chafer Holotrichia parallela Motschulsky (Coleoptera: Scarabaeidae). Comp. Biochem. Physiol. Part D Genom. Proteom. 2018, 28, 63–71. [Google Scholar] [CrossRef]
  33. Cui, X.; Liu, D.; Sun, K.; He, Y.; Shi, X. Expression Profiles and Functional Characterization of Two Odorant-Binding Proteins From the Apple Buprestid Beetle Agrilus mali (Coleoptera: Buprestidae). J. Econ. Entomol. 2018, 111, 1420–1432. [Google Scholar] [CrossRef] [PubMed]
  34. Zhang, X.; Yang, S.; Zhang, J.; Wang, X.; Wang, S.; Liu, M.; Xi, J. Identification and expression analysis of candidate chemosensory receptors based on the antennal transcriptome of Lissorhoptrus oryzophilus. Comp. Biochem. Physiol. Part D Genom. Proteom. 2019, 30, 133–142. [Google Scholar] [CrossRef]
  35. Andersson, M.N.; Grosse-Wilde, E.; Keeling, C.I.; Bengtsson, J.M.; Yuen, M.M.; Li, M.; Hillbur, Y.; Bohlmann, J.; Hansson, B.S.; Schlyter, F. Antennal transcriptome analysis of the chemosensory gene families in the tree killing bark beetles, Ips typographus and Dendroctonus ponderosae (Coleoptera: Curculionidae: Scolytinae). BMC Genom. 2013, 14, 198. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Gu, X.C.; Zhang, Y.N.; Kang, K.; Dong, S.L.; Zhang, L.W. Antennal Transcriptome Analysis of Odorant Reception Genes in the Red Turpentine Beetle (RTB), Dendroctonus valens. PLoS ONE 2015, 10, e0125159. [Google Scholar] [CrossRef] [Green Version]
  37. Li, X.; Ju, Q.; Jie, W.; Li, F.; Jiang, X.; Hu, J.; Qu, M. Chemosensory gene families in adult antennae of Anomala corpulenta Motschulsky (Coleoptera: Scarabaeidae: Rutelinae). PLoS ONE 2015, 10, e0121504. [Google Scholar] [CrossRef]
  38. Annemarie, H.; Hanghui, K.; Irina, B.; Michael, V.T.; Till, T.; Armin, G.T.; Elisabeth, W.; Wittko, F.; Ulrich, M.; Stefan, D. Deceptive Ceropegia dolichophylla fools its kleptoparasitic fly pollinators with exceptional floral scent. Front. Ecol. Environ. 2015, 3, 66. [Google Scholar] [CrossRef] [Green Version]
  39. Liu, S.; Rao, X.J.; Li, M.Y.; Feng, M.F.; He, M.Z.; Li, S.G. Identification of candidate chemosensory genes in the antennal transcriptome of Tenebrio molitor (Coleoptera: Tenebrionidae). Comp. Biochem. Physiol. Part D Genom. Proteom. 2015, 13, 44–51. [Google Scholar] [CrossRef]
  40. Gonzalez, F.; Johny, J.; Walker, W.B.; Guan, Q.; Mfarrej, S.; Jakše, J.; Montagné, N.; Jacquin-Joly, E.; Alqarni, A.S.; Al-Saleh, M.A.; et al. Antennal transcriptome sequencing and identification of candidate chemoreceptor proteins from an invasive pest, the American palm weevil, Rhynchophorus palmarum. Sci. Rep. 2021, 11, 8334. [Google Scholar] [CrossRef]
  41. Mitchell, R.F.; Hughes, D.T.; Luetje, C.W.; Millar, J.G.; Soriano-Agatón, F.; Hanks, L.M.; Robertson, H.M. Sequencing and characterizing odorant receptors of the cerambycid beetle Megacyllene caryae. Insect Biochem. Mol. Biol. 2012, 42, 499–505. [Google Scholar] [CrossRef] [Green Version]
  42. Sun, L.; Zhang, Y.N.; Qian, J.L.; Kang, K.; Zhang, X.Q.; Deng, J.D.; Tang, Y.P.; Chen, C.; Hansen, L.; Xu, T.; et al. Identification and Expression Patterns of Anoplophora chinensis (Forster) Chemosensory Receptor Genes from the Antennal Transcriptome. Front. Physiol. 2018, 9, 90. [Google Scholar] [CrossRef] [Green Version]
  43. Qian, J.L.; Mang, D.Z.; Lv, G.C.; Ye, J.; Li, Z.Q.; Chu, B.; Sun, L.; Liu, Y.J.; Zhang, L.W. Identification and Expression Profile of Olfactory Receptor Genes Based on Apriona germari (Hope) Antennal Transcriptome. Front. Physiol. 2020, 11, 807. [Google Scholar] [CrossRef]
  44. Hu, P.; Wang, J.; Cui, M.; Tao, J.; Luo, Y. Antennal transcriptome analysis of the Asian longhorned beetle Anoplophora glabripennis. Sci. Rep. 2016, 6, 26652. [Google Scholar] [CrossRef] [Green Version]
  45. Wang, J.; Li, D.Z.; Min, S.F.; Mi, F.; Zhou, S.S.; Wang, M.Q. Analysis of chemosensory gene families in the beetle Monochamus alternatus and its parasitoid Dastarcus helophoroides. Comp. Biochem. Physiol. Part D Genom. Proteom. 2014, 11, 1–8. [Google Scholar] [CrossRef]
  46. Xu, T.; Yasui, H.; Teale, S.A.; Fujiwara-Tsujii, N.; Wickham, J.D.; Fukaya, M.; Hansen, L.; Kiriyama, S.; Hao, D.; Nakano, A.; et al. Identification of a male-produced sex-aggregation pheromone for a highly invasive cerambycid beetle, Aromia bungii. Sci. Rep. 2017, 7, 7330. [Google Scholar] [CrossRef]
  47. Pea, E.D.L.; Schrader, G.; Vos, S. Pest survey card on Aromia bungii. EFSA Support. Publ. 2019, 16, EN-1731. [Google Scholar] [CrossRef] [Green Version]
  48. Yasui, H.; Fujiwara-Tsujii, N.; Yasuda, T.; Fukaya, M.; Kiriyama, S.; Nakano, A.; Watanabe, T.; Mori, K. Electroantennographic responses and field attraction of an emerging invader, the red-necked longicorn beetle Aromia bungii (Coleoptera: Cerambycidae), to the chiral and racemic forms of its male-produced aggregation-sex pheromone. Appl. Entomol. Zool. 2018, 54, 109–114. [Google Scholar] [CrossRef]
  49. Wang, W.C.; Cao, D.D.; Men, J. (R)-(+)-citronellal identified as a female-produced sex pheromone of Aromia bungii Faldermann (Coleoptera: Cerambycidae). Egypt. J. Biol. Pest Control 2018, 28, 77. [Google Scholar] [CrossRef]
  50. Antony, B.; Johny, J.; Aldosari, S.A. Silencing the Odorant Binding Protein RferOBP1768 Reduces the Strong Preference of Palm Weevil for the Major Aggregation Pheromone Compound Ferrugineol. Front. Physiol. 2018, 9, 252. [Google Scholar] [CrossRef]
  51. Krogh, A.; Larsson, B.; von Heijne, G.; Sonnhammer, E.L. Predicting transmembrane protein topology with a hidden Markov model: Application to complete genomes. J. Mol. Biol. 2001, 305, 567–580. [Google Scholar] [CrossRef] [Green Version]
  52. Mamidala, P.; Wijeratne, A.J.; Wijeratne, S.; Poland, T.; Qazi, S.S.; Doucet, D.; Cusson, M.; Beliveau, C.; Mittapalli, O. Identification of odor-processing genes in the emerald ash borer, Agrilus planipennis. PLoS ONE 2013, 8, e56555. [Google Scholar] [CrossRef] [Green Version]
  53. Engsontia, P.; Sanderson, A.P.; Cobb, M.; Walden, K.K.; Robertson, H.M.; Brown, S. The red flour beetle’s large nose: An expanded odorant receptor gene family in Tribolium castaneum. Insect Biochem. Mol. Biol. 2008, 38, 387–397. [Google Scholar] [CrossRef]
  54. Wu, Z.; Bin, S.; He, H.; Wang, Z.; Li, M.; Lin, J. Differential Expression Analysis of Chemoreception Genes in the Striped Flea Beetle Phyllotreta striolata Using a Transcriptomic Approach. PLoS ONE 2016, 11, e0153067. [Google Scholar] [CrossRef] [Green Version]
  55. Tamura, K.; Peterson, D.; Peterson, N.; Stecher, G.; Nei, M.; Kumar, S. MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 2011, 28, 2731–2739. [Google Scholar] [CrossRef] [Green Version]
  56. Simon, P. Q-Gene: Processing quantitative real-time RT-PCR data. Bioinformatics 2003, 19, 1439–1440. [Google Scholar] [CrossRef] [Green Version]
  57. Ebbs, M.L.; Amrein, H. Taste and pheromone perception in the fruit fly Drosophila melanogaster. Pflug. Arch. 2007, 454, 735–747. [Google Scholar] [CrossRef] [Green Version]
  58. Plewa, R.; Jaworski, T.; Hilszczański, J.; Horák, J. Investigating the biodiversity of the forest strata: The importance of vertical stratification to the activity and development of saproxylic beetles in managed temperate deciduous forests. For. Ecol. Manag. 2017, 402, 186–193. [Google Scholar] [CrossRef]
  59. Sweeney, J.D.; Silk, P.J.; Gutowski, J.M.; Wu, J.; Lemay, M.A.; Mayo, P.D.; Magee, D.I. Effect of chirality, release rate, and host volatiles on response of Tetropium fuscum (F.), Tetropium cinnamopterum Kirby, and Tetropium castaneum (L.) to the aggregation pheromone, Fuscumol. J. Chem. Ecol. 2010, 36, 1309–1321. [Google Scholar] [CrossRef]
  60. Ando, T.; Inomata, S.; Yamamoto, M. Lepidopteran sex pheromones. Top. Curr. Chem. 2004, 239, 51–96. [Google Scholar] [CrossRef]
  61. Abuin, L.; Bargeton, B.; Ulbrich, M.H.; Isacoff, E.Y.; Kellenberger, S.; Benton, R. Functional architecture of olfactory ionotropic glutamate receptors. Neuron 2011, 69, 44–60. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. 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 Genom. Proteom. 2017, 24, 1–11. [Google Scholar] [CrossRef]
  63. Mitchell, R.F.; Schneider, T.M.; Schwartz, A.M.; Andersson, M.N.; McKenna, D.D. The diversity and evolution of odorant receptors in beetles (Coleoptera). Insect Mol. Biol. 2020, 29, 77–91. [Google Scholar] [CrossRef]
  64. Zhao, T.; Ganji, S.; Schiebe, C.; Bohman, B.; Weinstein, P.; Krokene, P.; Borg-Karlson, A.K.; Unelius, C.R. Convergent evolution of semiochemicals across Kingdoms: Bark beetles and their fungal symbionts. ISME J. 2019, 13, 1535–1545. [Google Scholar] [CrossRef]
  65. Hou, X.Q.; Yuvaraj, J.K.; Roberts, R.E.; Zhang, D.D.; Unelius, C.R.; Löfstedt, C.; Andersson, M.N. Functional Evolution of a Bark Beetle Odorant Receptor Clade Detecting Monoterpenoids of Different Ecological Origins. Mol. Biol. Evol. 2021, 38, 4934–4947. [Google Scholar] [CrossRef]
  66. Abdel-Latief, M. A family of chemoreceptors in Tribolium castaneum (Tenebrionidae: Coleoptera). PLoS ONE 2007, 2, e1319. [Google Scholar] [CrossRef]
  67. Tang, Q.; Zhang, Y.; Shen, C.; Xia, D. Identification and Expression Profiling of Odorant Receptor Protein Genes in Sitophilus zeamais (Coleoptera: Curculionoidea) Using RT-qPCR. Neotrop. Entomol. 2019, 48, 538–551. [Google Scholar] [CrossRef]
  68. Xu, W.; Anderson, A. Carbon dioxide receptor genes in cotton bollworm Helicoverpa armigera. Naturwissenschaften 2015, 102, 11. [Google Scholar] [CrossRef]
  69. Bin, S.Y.; Qu, M.Q.; Li, K.M.; Peng, Z.Q.; Wu, Z.Z.; Lin, J.T. Antennal and abdominal transcriptomes reveal chemosensory gene families in the coconut hispine beetle, Brontispa longissima. Sci. Rep. 2017, 7, 2809. [Google Scholar] [CrossRef] [Green Version]
  70. Kozma, M.T.; Schmidt, M.; Ngo-Vu, H.; Sparks, S.D.; Senatore, A.; Derby, C.D. Chemoreceptor proteins in the Caribbean spiny lobster, Panulirus argus: Expression of Ionotropic Receptors, Gustatory Receptors, and TRP channels in two chemosensory organs and brain. PLoS ONE 2018, 13, e0203935. [Google Scholar] [CrossRef]
  71. Croset, V.; Rytz, R.; Cummins, S.F.; Budd, A.; Brawand, D.; Kaessmann, H.; Gibson, T.J.; Benton, R. Ancient protostome origin of chemosensory ionotropic glutamate receptors and the evolution of insect taste and olfaction. PLoS Genet. 2010, 6, e1001064. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Homology analysis of A. bungii unigenes. (a): Similarity distribution; (b): E-value distribution; (c): species classification. All unigenes that had BLASTX annotations within the NCBI nr database with a cutoff E-value of 10−5 were analyzed. The first hit of each sequence was used for analysis.
Figure 1. Homology analysis of A. bungii unigenes. (a): Similarity distribution; (b): E-value distribution; (c): species classification. All unigenes that had BLASTX annotations within the NCBI nr database with a cutoff E-value of 10−5 were analyzed. The first hit of each sequence was used for analysis.
Insects 13 00096 g001
Figure 2. Gene ontology (GO) classification of A. bungii unigenes. Gene ontology (GO) assignment of Aromia bungii unigenes. The GO classification map was made by uploading the GO ID numbers of genes for their involvement in biological processes, cellular components and molecular functions.
Figure 2. Gene ontology (GO) classification of A. bungii unigenes. Gene ontology (GO) assignment of Aromia bungii unigenes. The GO classification map was made by uploading the GO ID numbers of genes for their involvement in biological processes, cellular components and molecular functions.
Insects 13 00096 g002
Figure 3. KEGG classification of A. bungii unigenes. Kyoto Encyclopedia of Genes and Genomes (KEGG) classification of Aromia bungii unigenes. The x-axis indicates the percentage of annotated genes, and the y-axis indicates the KEGG categories. The capital letters against the colored bars indicate five main categories: (a) cellular processes, (b) environmental information processing, (c) genetic information processing, (d) metabolism and (e) organism systems.
Figure 3. KEGG classification of A. bungii unigenes. Kyoto Encyclopedia of Genes and Genomes (KEGG) classification of Aromia bungii unigenes. The x-axis indicates the percentage of annotated genes, and the y-axis indicates the KEGG categories. The capital letters against the colored bars indicate five main categories: (a) cellular processes, (b) environmental information processing, (c) genetic information processing, (d) metabolism and (e) organism systems.
Insects 13 00096 g003
Figure 4. KOG classification of A. bungii unigenes. Eukaryotic Ortholog Groups (KOG) were divided into 26 groups, and the genes with successful KOG annotations were classified according to the KOG group. The horizontal axis is the names of the 26 groups of KOG, and the vertical axis is the proportion of the number of genes annotated to the group to the total number of genes annotated.
Figure 4. KOG classification of A. bungii unigenes. Eukaryotic Ortholog Groups (KOG) were divided into 26 groups, and the genes with successful KOG annotations were classified according to the KOG group. The horizontal axis is the names of the 26 groups of KOG, and the vertical axis is the proportion of the number of genes annotated to the group to the total number of genes annotated.
Insects 13 00096 g004
Figure 5. Phylogenetic tree of A. bungii OR genes. Molecular phylogeny comparing AbunORs with odorant receptors (ORs) from 14 other insect species: a total of 45 ORs (AbunOR1–45) from Aromia bungii (Abun) and ORs from Apriona germari (Ager), Anoplophora chinensis (Achi), Anoplophora glabripennis (Agla), Agrilus planipennis (Apla), Tribolium castaneum (Tcas), Dendroctonus ponderosae (Dpon), Ips typographus (Ityp) Anomala corpulenta (Acor), Megacyllene caryae (Mcar), Monchamus alternatus (Malt), Phyllotreta striolata (Pstr), Dendroctonus valens (Dval), Colaphellus bowringi (Cbow) and Tenebrio molitor (Tmol).
Figure 5. Phylogenetic tree of A. bungii OR genes. Molecular phylogeny comparing AbunORs with odorant receptors (ORs) from 14 other insect species: a total of 45 ORs (AbunOR1–45) from Aromia bungii (Abun) and ORs from Apriona germari (Ager), Anoplophora chinensis (Achi), Anoplophora glabripennis (Agla), Agrilus planipennis (Apla), Tribolium castaneum (Tcas), Dendroctonus ponderosae (Dpon), Ips typographus (Ityp) Anomala corpulenta (Acor), Megacyllene caryae (Mcar), Monchamus alternatus (Malt), Phyllotreta striolata (Pstr), Dendroctonus valens (Dval), Colaphellus bowringi (Cbow) and Tenebrio molitor (Tmol).
Insects 13 00096 g005
Figure 6. Expression levels of A. bungii ORs in different tissues of female and male adults. Relative mRNA expression of AbunORs in Aromia bungii tissues. The relative mRNA levels were normalized to those of the actin gene and analyzed using the Q-gene method. All values are shown as mean ± SEM normalized. The data were analyzed via least significant difference test after one-way analysis of variance. FA: female antennae; MA: male antennae; FM: female mouthpart (maxillary palps and labial palps); MM: male mouthpart (maxillary palps and labial palps); FAb: female abdominal end; MAb: male abdominal end. Different letters (a–d) indicate significant differences between means (p < 0.05).
Figure 6. Expression levels of A. bungii ORs in different tissues of female and male adults. Relative mRNA expression of AbunORs in Aromia bungii tissues. The relative mRNA levels were normalized to those of the actin gene and analyzed using the Q-gene method. All values are shown as mean ± SEM normalized. The data were analyzed via least significant difference test after one-way analysis of variance. FA: female antennae; MA: male antennae; FM: female mouthpart (maxillary palps and labial palps); MM: male mouthpart (maxillary palps and labial palps); FAb: female abdominal end; MAb: male abdominal end. Different letters (a–d) indicate significant differences between means (p < 0.05).
Insects 13 00096 g006
Figure 7. Phylogenetic tree of A. bungii GR genes. Molecular phylogeny comparing AbunGRs with gustatory receptors (GRs) from 8 other insect species: A total of 6 GRs (AbunGR1–6) from Aromia bungii (Abun) and ORs from Anoplophora chinensis (Achi), Ips typographus (Ityp), Anoplophora glabripennis (Agla), Drosophila melanogaster (Dmel), Dendroctonus ponderosae (Dpon), Bombyx mori (Bmor), Tribolium castaneum (Tcas) and Phyllotreta striolata (Pstr).
Figure 7. Phylogenetic tree of A. bungii GR genes. Molecular phylogeny comparing AbunGRs with gustatory receptors (GRs) from 8 other insect species: A total of 6 GRs (AbunGR1–6) from Aromia bungii (Abun) and ORs from Anoplophora chinensis (Achi), Ips typographus (Ityp), Anoplophora glabripennis (Agla), Drosophila melanogaster (Dmel), Dendroctonus ponderosae (Dpon), Bombyx mori (Bmor), Tribolium castaneum (Tcas) and Phyllotreta striolata (Pstr).
Insects 13 00096 g007
Figure 8. Expression levels of A. bungii GRs in different tissues of female and male adults. Relative mRNA expression of AbunGRs in Aromia bungii tissues. The relative mRNA levels were normalized to those of the actin gene and analyzed using the Q-gene method. All values are shown as mean ± SEM normalized. The data were analyzed via least significant difference test after one-way analysis of variance. FA: female antennae; MA: male antennae; FM: female mouthpart (maxillary palps and labial palps); MM: male mouthpart (maxillary palps and labial palps); FAb: female abdominal end; MAb: male abdominal end. Different letters (a, b) indicate significant differences between means (p < 0.05).
Figure 8. Expression levels of A. bungii GRs in different tissues of female and male adults. Relative mRNA expression of AbunGRs in Aromia bungii tissues. The relative mRNA levels were normalized to those of the actin gene and analyzed using the Q-gene method. All values are shown as mean ± SEM normalized. The data were analyzed via least significant difference test after one-way analysis of variance. FA: female antennae; MA: male antennae; FM: female mouthpart (maxillary palps and labial palps); MM: male mouthpart (maxillary palps and labial palps); FAb: female abdominal end; MAb: male abdominal end. Different letters (a, b) indicate significant differences between means (p < 0.05).
Insects 13 00096 g008
Figure 9. Phylogenetic tree and relative mRNA expression of AbunIRs. (a) A total of 2 IRs (AbunIR1–2) from Aromia bungii (Abun) and Ors from Anoplophora chinensis (Achi), Anoplophora glabripennis (Agla), Dendroctonus ponderosae (Dpon), Drosophila melanogaster (Dmel), Phyllotreta striolata (Pstr), Anomala corpulenta (Acor), Monochamus alternatus (Malt), Megacyllene caryae (Mcar) and Tenebrio molitor (Tmol) were used to construct the phylogenetic tree. (b) Expression levels of A. bungii Irs in different tissues of female and male adults. Relative mRNA expression of AbunIRs in Aromia bungii tissues. The relative mRNA levels were normalized to those of the actin gene and analyzed using the Q-gene method. All values are shown as mean ± SEM normalized. The data were analyzed via least significant difference test after one-way analysis of variance. FA: female antennae; MA: male antennae; FM: female mouthpart (maxillary palps and labial palps); MM: male mouthpart (maxillary palps and labial palps); Fab: female abdominal end; Mab: male abdominal end. Different letters (a, b) indicate significant differences between means (p < 0.05).
Figure 9. Phylogenetic tree and relative mRNA expression of AbunIRs. (a) A total of 2 IRs (AbunIR1–2) from Aromia bungii (Abun) and Ors from Anoplophora chinensis (Achi), Anoplophora glabripennis (Agla), Dendroctonus ponderosae (Dpon), Drosophila melanogaster (Dmel), Phyllotreta striolata (Pstr), Anomala corpulenta (Acor), Monochamus alternatus (Malt), Megacyllene caryae (Mcar) and Tenebrio molitor (Tmol) were used to construct the phylogenetic tree. (b) Expression levels of A. bungii Irs in different tissues of female and male adults. Relative mRNA expression of AbunIRs in Aromia bungii tissues. The relative mRNA levels were normalized to those of the actin gene and analyzed using the Q-gene method. All values are shown as mean ± SEM normalized. The data were analyzed via least significant difference test after one-way analysis of variance. FA: female antennae; MA: male antennae; FM: female mouthpart (maxillary palps and labial palps); MM: male mouthpart (maxillary palps and labial palps); Fab: female abdominal end; Mab: male abdominal end. Different letters (a, b) indicate significant differences between means (p < 0.05).
Insects 13 00096 g009
Table 1. Summary of A. bungii antennae transcriptome.
Table 1. Summary of A. bungii antennae transcriptome.
Statistics ProjectNumber
Total Raw Reads45,642,924
Total Clean Reads43,302,906
Clean bases6.5G
Q20 percentage97.74%
Q30 percentage94.10%
GC percentage43.80%
Transcripts79,280
Mean length of transcripts783
N50 of transcripts1435
Unigenes42,197
Mean length of Unigenes1215
N50 of Unigenes1744
Table 2. The Blastx match of A. bungii OR genes.
Table 2. The Blastx match of A. bungii OR genes.
Gene NameORF
Length (bp)
Complete
ORF
Transmembrane HelixFPKM Best Blastx Match
NameAcc.numberSpeciesE-ValueIdentity (%)
AbunOR1330Yes11.64olfactory receptor OR1AJO62220.1Tenebrio molitor4 × 10−2243%
AbunOR21182Yes412.81odorant receptor Or1-likeXP_018566530.1Anoplophora glabripennis3 × 10−12148%
AbunOR3330Yes01.64odorant receptor Or2-likeXP_018567067.1Anoplophora glabripennis1 × 10−872%
AbunOR4402Yes21.17odorant receptor 4-likeXP_023027241.1Leptinotarsa decemlineata6 × 10−3749%
AbunOR51308Yes71.71odorant receptorAUF73018.1Anoplophora chinensis9 × 10−16255%
AbunOR61179Yes60.69odorant receptor Or2-likeXP_018579015.1Anoplophora glabripennis1 × 10−8139%
AbunOR7594Yes33.89odorant receptor Or2-likeXP_018576526.1Anoplophora glabripennis1 × 10−9569%
AbunOR8600Yes22.4odorant receptor Or1-likeXP_023310030.1Anoplophora glabripennis2 × 10−9078%
AbunOR91137Yes70.33odorant receptor 22b-likeXP_023311541.1Anoplophora glabripennis6 × 10−4533%
AbunOR10276Yes00.76odorant receptor 63a-likeXP_023016125.1Leptinotarsa decemlineata4 × 10−1746%
AbunOR11150Yes11.81odorant receptor 22c-likeXP_018579026.1Anoplophora glabripennis3 × 10−1053%
AbunOR12996Yes43.14odorant receptor 63a-likeXP_023016125.1Leptinotarsa decemlineata8 × 10−5035%
AbunOR13630Yes03.28odorant receptorAUF73043.1Anoplophora chinensis2 × 10−4437%
AbunOR1490Yes052.82odorant receptor 49b-likeXP_018577261.1Anoplophora glabripennis4 × 10−1374%
AbunOR15897NO43.27odorant receptor 49b-likeXP_022917715.1Onthophagus taurus8 × 10−2829%
AbunOR16507NO21.58odorant receptor 1APC94305.1Pyrrhalta aenescens3 × 10−3339%
AbunOR171086Yes73.44odorant receptor 94a-likeXP_018560823.1Anoplophora glabripennis5 × 10−3230%
AbunOR18486NO01.04odorant receptor 4-likeXP_018577142.1Anoplophora glabripennis2 × 10−3445%
AbunOR19411Yes32.27odorant receptorAUF73037.1Anoplophora chinensis4 × 10−3848%
AbunOR20321Yes13.51odorant receptor OR32ALR72575.1Colaphellus bowringi6 × 10−1735%
AbunOR21654Yes31.9odorant receptor 4-likeXP_023027241.1Leptinotarsa decemlineata3 × 10−3733%
AbunOR22366Yes13.25odorant receptor Or2-likeXP_018579015.1Anoplophora glabripennis6 × 10−4850%
AbunOR23576Yes23.7odorant receptor 94a-likeXP_018560823.1Anoplophora glabripennis4 × 10−4736%
AbunOR24297Yes02.21odorant receptor 63a-likeXP_023016125.1Leptinotarsa decemlineata6 × 10−1538%
AbunOR251440Yes749.99odorant receptor coreceptorXP_018568191.1Anoplophora glabripennis092%
AbunOR26936Yes46.25odorant receptor 18APC94230.1Pyrrhalta maculicollis9 × 10−9041%
AbunOR271158Yes60.66odorant receptor Or2-likeXP_018579015.2Anoplophora glabripennis5 × 10−9541%
AbunOR28894Yes42.79odorant receptorAUF73043.1Anoplophora chinensis2 × 10−6441%
AbunOR29735Yes52.76odorant receptor OR36ALR72579.1Colaphellus bowringi4 × 10−11361%
AbunOR30318Yes21.58odorant receptor Or2-likeXP_018576526.1Anoplophora glabripennis1 × 10−3155%
AbunOR311134Yes72.69odorant receptor 4-likeXP_018577142.1Anoplophora glabripennis2 × 10−8841%
AbunOR321143Yes61.34odorant receptor 19, partialAVN97831.1Anoplophora chinensis3 × 10−9943%
AbunOR33762Yes13.09odorant receptor 18APC94230.1Pyrrhalta maculicollis6 × 10−10554%
AbunOR34381Yes11.83odorant receptor 43a-likeXP_018573343.1Anoplophora glabripennis7 × 10−3246%
AbunOR351170Yes63.22odorant receptorAUF73043.1Anoplophora chinensis2 × 10−4528%
AbunOR36171Yes22.63odorant receptor 85b-likeXP_018564120.1Anoplophora glabripennis7 × 10−1357%
AbunOR37726Yes45.48odorant receptor OR24ALR72568.1Colaphellus bowringi3 × 10−9247%
AbunOR381104Yes73.11odorant receptor 49b-likeXP_018570955.1Anoplophora glabripennis1 × 10−17165%
AbunOR39627Yes20.96odorant receptor 286EFA01418.1Tribolium castaneum6 × 10−1726%
AbunOR40387Yes00.69odorant receptor OR20ALR72565.1Colaphellus bowringi4 × 10−2639%
AbunOR41840Yes41.06odorant receptor OR9ALR72554.1Colaphellus bowringi3 × 10−8743%
AbunOR42372NO21.61odorant receptor 19, partialAVN97831.1Anoplophora chinensis4 × 10−935%
AbunOR43333Yes21.87odorant receptor 26QNH68050.1Apriona germari4 × 10−7869%
AbunOR44537NO21.13odorant receptor Or2-likeXP_023311850.1Anoplophora glabripennis4 × 10−6856%
AbunOR45543Yes22.00odorant receptor 3, partialAVN97815.1Anoplophora chinensis5 × 10−3435%
Table 3. The Blastx match of A. bungii GR genes.
Table 3. The Blastx match of A. bungii GR genes.
Gene NameORF
Length
(bp)
Complete ORFTransmembrane HelixFPKM ValueBest Blastx Match
NameAcc.numberSpeciesE-ValueIdentity (%)
AbunGR1243NO02.59putative gustatory receptor GR9ALR72586.1Colaphellus bowringi1 × 10−1744%
AbunGR2198NO08.34gustatory receptor Gr83NP_001138948.1Tribolium castaneum1 × 10−1866%
AbunGR3297NO01.77putative gustatory receptor 2aXP_015840061.1Tribolium castaneum1 × 10−938%
AbunGR4597NO 33.22gustatory receptor 68a-likeXP_018567270.1Anoplophora glabripennis6 × 10−3739%
AbunGR5522NO41.53gustatory receptor 1EFA07594.2Tribolium castaneum2 × 10−11581%
AbunGR6510NO24.1gustatory receptor 3AKC58580.2Anomala corpulenta9 × 10−10383%
Table 4. The Blastx match of A. bungii IR genes.
Table 4. The Blastx match of A. bungii IR genes.
Gene NameORF
Length (bp)
Complete
ORF
Transmembrane HelixFPKM ValueBest Blastx Match
NameAcc.numberSpeciesE-ValueIdentity (%)
AbunIR1174NO02.59ionotropic receptor 1QNH68025.1Apriona germari1 × 10−1352%
AbunIR21257Yes36.03ionotropic receptor 3ANQ46495.1Phyllotreta striolata1 × 10−11044%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wu, Z.; Ye, J.; Qian, J.; Purba, E.R.; Zhang, Q.; Zhang, L.; Mang, D. Identification and Expression Profile of Chemosensory Receptor Genes in Aromia bungii (Faldermann) Antennal Transcriptome. Insects 2022, 13, 96. https://0-doi-org.brum.beds.ac.uk/10.3390/insects13010096

AMA Style

Wu Z, Ye J, Qian J, Purba ER, Zhang Q, Zhang L, Mang D. Identification and Expression Profile of Chemosensory Receptor Genes in Aromia bungii (Faldermann) Antennal Transcriptome. Insects. 2022; 13(1):96. https://0-doi-org.brum.beds.ac.uk/10.3390/insects13010096

Chicago/Turabian Style

Wu, Zhenchen, Jia Ye, Jiali Qian, Endang Rinawati Purba, Qinghe Zhang, Longwa Zhang, and Dingze Mang. 2022. "Identification and Expression Profile of Chemosensory Receptor Genes in Aromia bungii (Faldermann) Antennal Transcriptome" Insects 13, no. 1: 96. https://0-doi-org.brum.beds.ac.uk/10.3390/insects13010096

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop