{'fraction_non_alphanumeric': 0.040793740670201264, 'fraction_numerical': 0.025454264683172904, 'mean_word_length': 5.2138173676635216, 'pattern_counts': {'":': 0, '<': 1, '<?xml version=': 0, '>': 2, 'https://': 3, 'lorem ipsum': 0, 'www.': 2, 'xml': 0}, 'pii_count': 1, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 8, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Identification and Expression Profile of Chemosensory Receptor Genes in Aromia bungii (Faldermann) Antennal Transcriptome

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.


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), odorantdegrading 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.
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,7epoxynonenal [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.

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.

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].

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 nonredundant 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 Insects 2022, 13, 96 4 of 20 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://www.ncbi.nlm.nih.gov/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].

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.

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).

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' anten- nae. 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.  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).
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).  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 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.

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.
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).

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). abscess axis A-Z represents 26 types: the largest proportion of genes was involved in ge eral function prediction alone (1929 unigenes; 16.37%); the number of unigenes in th KOG database was the largest, followed by signal transduction mechanisms (1424 un genes; 12.08%) and posttranslational modification/protein conversion/chaperones (126 10.73%) (Figure 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

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 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) (Figures 6 and S2). 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) (Figures 6 and S2).

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).

Identification of Putative Gustatory Receptors
Six AbunGRs (AbunGR1-6) were identified by bioinformatics analysis of the naal transcriptome. Five AbunGRs (AbunGR1-4, AbunGR6) had a complete open r frame (ORF). Transmembrane domain prediction results showed that three Ab (AbunGR4-6) had the predicted transmembrane helix structure (Table 3). In the G lutionary tree, the six AbunGRs were divided into four taste receptor gene fa namely sugar, fructose, bitter and carbon dioxide receptor families. Only AbunG clustered with the reported carbon dioxide receptors DmelGR21a and DmelGR6 (Figure 7).

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).

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).

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).

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).

Discussion
As the largest order of Insecta, Coleoptera is also one of the most important pe ders for forestry and agriculture. Among the groups of Coleoptera, many ceramb (both native or invasive species) have become serious forest (or tree) pest insects, cau significant ecological and economical losses throughout the world [58]. Over the past ade, significant research progresses have been made regarding the chemical ecolo longhorn beetles, especially the identification of aggregation-sex pheromones (attract for some economically important species. However, compared to the well-studied dopterans, research on olfactory mechanisms at the molecular/gene levels in Ceram dae still remains limited [59,60]. In the current study, we identified chemosensory rec genes and studied their expression profiles based on the antennal transcriptome dat highly damaging pest of Prunus fruit trees, the red-necked longhorn beetle, A. bungi A total of 42,197 unigenes were identified from the antennal transcriptome da which 74.62% were more than 500 bp in length. The NCBI-nr database had the la number of unigenes (with homologous sequences) with a total of 25,407 unigene. A. b showed the highest homology with T. castaneum (60.4%), followed by D. ponder (17.7%). The number of unigenes successfully annotated in the GO database was l than that in the NCBI-nr database, with a total of 19,959 unigenes. Among them, bin and catalytic activities were the most annotated, which are similar to the functions factory-related genes. In addition, 12,154 unigenes were annotated in the KO data and KEGG metabolic pathway classification showed that the subcategory of "signal t duction" is the most annotated gene pathway in the metabolic category of "environm information processing" (Figure 2). These results can be related to the functio

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 olfactoryrelated 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 sexspecific 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 CO 2 [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://www.mdpi.com/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 Table S1: Primers of A. bungii chemosensory receptor genes used for RT-qPCR, Table S2: Summary of the gene annotation success ratio,