Expression Profiling and Functional Analysis of Candidate Odorant Receptors in Galeruca daurica

Simple Summary Odorant receptors (ORs) play an important role in the olfactory system in insects. However, there is no functional research on the odorant receptors of Galeruca daurica. In this study, 21 OR genes were identified from the transcriptome database of G. daurica adults. Most GdauORs were mainly expressed in antennae, and the expression levels of GdauORs in adults were affected by age. When GdauOR4, GdauOR15, and GdauORco were silenced by RNAi, the electrophysiological responses to host plant volatiles were significantly decreased. Abstract Galeruca daurica (Joannis) is an oligophagous pest in the grasslands of Inner Mongolia, China, which feed mainly on Allium spp. Odorant receptors (ORs) play an important role in the olfactory system in insects, and function together with olfactory co-receptor (ORco). In this study, 21 OR genes were identified from the transcriptome database of G. daurica adults, and named GdauOR1-20 and GdauORco. The expression profiles were examined by RT-qPCR and RNA interference (RNAi) and electroantennogram (EAG) experiments were conducted to further identify the olfactory functions of GdauOR4, GdauOR11, GdauOR15, and GdauORco. It was found that 15 GdauORs (OR1, OR3-6, OR8, OR11-13, OR15, OR17-20, and ORco) were mainly expressed in antennae, and the expression levels of GdauORs in adults were affected by age. When GdauOR4, GdauOR15, and GdauORco were silenced by RNAi, the electrophysiological responses to host plant volatiles were significantly decreased in G. daurica. This study lays a necessary foundation for clarifying the mechanism on finding host plants in G. daurica.


Introduction
A sophisticated olfactory system is a key physiological element for survival and reproduction in insects [1]. In the olfactory system, insects detect semiochemicals through interactions with various olfactory proteins, including odorant-binding proteins (OBPs) and chemosensory proteins (CSPs) which bind odors in lymph fluid, odorant receptors (ORs), and ionotropic receptors (IRs) that convert chemical signals into nerve electrical signals through ion channels, as well as odorant-degrading enzymes (ODEs) that are considered to decompose odorants [1,2]. Olfactory co-receptor (Orco) was initially identified as a member of the OR family in Drosophila [3]. During the initial olfaction process, odorant receptors (ORs) expressed on the dendrites of olfactory sensory neurons (OSNs) play a central role in converting the chemical signals into electrical signals via the formation of heteromeric complexes that operate as odorant-gated ion channels with an Orco [4]. Most species of insects express just one Orco and a distinct complement of OR, ranging from just four members in Calopteryx splendens to more than 300 in Tribolium castaneum [5][6][7]. Numerous studies have shown that insect ORs are mainly expressed in olfactory organs, for

Quantitative Real-Time PCR (RT-qPCR) Measurement
Gene expression profiles were analyzed by RT-qPCR. All primers for RT-qPCR were designed by Primer Premier 5.0 (http://www.premierbiosoft.com/primerdesign/index. html, accessed on 10 December 2020) (Table S1). RT-qPCR were performed using the FTC-3000P Real-Time Quantitative Thermal Cycler (Funglyn Biotech, Markham, ON, Canada) with BRYT Green ® dye (GoTaq ® qRT-PCR Master Mix, Promega, Madison, WI, USA) as the fluorescence reporter. The succinate dehydrogenase complex gene (SDHA) was used as a reference gene [26]. A 10-fold dilution series of antennal cDNA was employed to construct a standard curve to determine the PCR efficiency. All amplification efficiencies in the RT-qPCR analysis ranged from 85.8% to 103%. Experiments were performed in a 10 µL reaction mixture, including 1 µL cDNA (10 ng/µL), 0.2 µL forward primer (10 µmol/L), 0.2 µL reverse primer (10 µmol/L), 5 µL BRYT Green dye, and 3.6 µL RNase free water, and repeated three times for each sample. The RT-qPCR reaction was carried out according to the following procedure: initial denaturation at 95 • C of 10 min, followed by 45 cycles of 95 • C for 15 s, and annealing at 60 • C for 1 min, and analysis of the melting curve at the last. There was a specific single peak in each reaction in the dissociation stage. The relative quantification of each gene expression level was analyzed using the 2 −∆∆Ct method [27]. The expression level of GdauOR1 in male antennae was used as the control in the expression profiles of different tissues. The expression level of each gene in 1-day-old male antennae was used as the control in the expression profiles of different days after eclosion.

RNA Interference of GdauOR4/GdauOR11/GdauOR15/GdauORco
RNA interference (RNAi) primers containing T7 promoter were designed to the coding sequence of GdauOR4, GdauOR11, GdauOR15, and GdauORco (Table S2). PCR was performed in a 25 µL reaction (1 µL cDNA, 1 µL forward primer, 1 µL reverse primer, 9.5 µL Premix Taq TM , and 12.5 µL RNase free water) with a thermocycler (T100 TM Thermal Cycler, Bio-Rad, Hercules, CA, USA), conditions consisting of: 94 • C for 3 min followed by 35 cycles of 94 • C for 30 s, 65 • C for 30 s (each gene used a primer-specific temperature), 72 • C for 1 min, and a final extension at 72 • C for 10 min. Amplification products were purified and connected to the pGEM-T Easy vector, and then transformed into the competent-cell Escherichia coli DH5α (TaKaRa, Dalian, China). Plasmids were sequence verified and used as templates to synthesize double-stranded RNAs (dsRNA) through the kit of T7 RiboMAX TM Express RNAi System (Promega, Madison, WI, USA). The quality of dsRNA was detected by 1.5% agarose gel electrophoresis. The concentration was detected by NanoDrop™ 2000 spectrophotometer. Finally, the dsRNA was diluted to 1000 ng/µL in enzyme-free water and stored at −80 • C until use.
The 3-day-old female and male adults were selected as test insects. Two µL dsRNA was injected into the intersegmental membrane between the fourth and fifth abdominal segments of G. daurica using a microinjector (Shimadzu, Kyoto City, Japan). The experiments were divided into 5 treatment groups (dsOR4-injected, dsOR11-injected, dsOR15-injected, dsORco-injected and dsGFP-injected) with 3 biological replicates per group and 20 individuals per replicate. The treated insects were reared under the above conditions. Antennae of GdauORco group were collected separately for interference efficiency measurement by RT-qPCR at 24,48,72, and 96 h after injection. Antennae of GdauOR4, GdauOR11, and GdauOR15 groups were collected for interference efficiency measurement at 48 h after injection. The group of dsGFP-injected was treated as control.

EAG Recordings
Electroantennogram (EAG) was used to record the antennal responses to 13 volatiles of A. mongolium (Table 1) [20]. These compounds were dissolved in mineral oil (Bio-Rad, USA) at 10 µg/µL. The mineral oil was also set up as a blank control. An antenna was cut off from the base, and cut 0.2 mm off the tip, then connected between the electrodes by two silver wires and two capillary glass tubes filled with physiological saline (0.9% sodium chloride injection). The antennal tip and base were inserted into the capillary tubes with 0.5 mm each. Filter paper strips (5 × 15 mm) were loaded with 10 µL of each solution, and transferred to the Pasteur pipette. The top of the pipette was inserted into the small hole in the wall of the tube, which was connected to an air stimulus controller (CS-55; SynTech, Kirchzarten, Germany). The signals were detected by a high-impedance amplifier (IDAC-2; SynTech, Kirchzarten, Germany) and analyzed using SynTech software (GC-EAD 2014 v1.2.5). The stimulation was started with a pulse duration of 0.2 s until the baseline stabilized. The flow rate was 4 mL/s. The EAG responses were recorded during odor stimulation, and six female or male antennae of 3-day-old adults at 48 h after injection were used for each volatile odorant recording. experiments were divided into 5 treatment groups (dsOR4-injected, dsOR11-injected, dsOR15-injected, dsORco-injected and dsGFP-injected) with 3 biological replicates per group and 20 individuals per replicate. The treated insects were reared under the above conditions. Antennae of GdauORco group were collected separately for interference efficiency measurement by RT-qPCR at 24, 48, 72, and 96 h after injection. Antennae of GdauOR4, GdauOR11, and GdauOR15 groups were collected for interference efficiency measurement at 48 h after injection. The group of dsGFP-injected was treated as control.

EAG Recordings
Electroantennogram (EAG) was used to record the antennal responses to 13 volatiles of A. mongolium (Table 1) [20]. These compounds were dissolved in mineral oil (Bio-Rad, USA) at 10 μg/μL. The mineral oil was also set up as a blank control. An antenna was cut off from the base, and cut 0.2 mm off the tip, then connected between the electrodes by two silver wires and two capillary glass tubes filled with physiological saline (0.9% sodium chloride injection). The antennal tip and base were inserted into the capillary tubes with 0.5 mm each. Filter paper strips (5 × 15 mm) were loaded with 10 μL of each solution, and transferred to the Pasteur pipette. The top of the pipette was inserted into the small hole in the wall of the tube, which was connected to an air stimulus controller (CS-55; SynTech, Kirchzarten, Germany). The signals were detected by a high-impedance amplifier (IDAC-2; SynTech, Kirchzarten, Germany) and analyzed using SynTech software (GC-EAD 2014 v1.2.5). The stimulation was started with a pulse duration of 0.2 s until the baseline stabilized. The flow rate was 4 mL/s. The EAG responses were recorded during odor stimulation, and six female or male antennae of 3-day-old adults at 48 h after injection were used for each volatile odorant recording. experiments were divided into 5 treatment groups (dsOR4-injected, dsOR11-injected, dsOR15-injected, dsORco-injected and dsGFP-injected) with 3 biological replicates per group and 20 individuals per replicate. The treated insects were reared under the above conditions. Antennae of GdauORco group were collected separately for interference efficiency measurement by RT-qPCR at 24, 48, 72, and 96 h after injection. Antennae of GdauOR4, GdauOR11, and GdauOR15 groups were collected for interference efficiency measurement at 48 h after injection. The group of dsGFP-injected was treated as control.

EAG Recordings
Electroantennogram (EAG) was used to record the antennal responses to 13 volatiles of A. mongolium (Table 1) [20]. These compounds were dissolved in mineral oil (Bio-Rad, USA) at 10 μg/μL. The mineral oil was also set up as a blank control. An antenna was cut off from the base, and cut 0.2 mm off the tip, then connected between the electrodes by two silver wires and two capillary glass tubes filled with physiological saline (0.9% sodium chloride injection). The antennal tip and base were inserted into the capillary tubes with 0.5 mm each. Filter paper strips (5 × 15 mm) were loaded with 10 μL of each solution, and transferred to the Pasteur pipette. The top of the pipette was inserted into the small hole in the wall of the tube, which was connected to an air stimulus controller (CS-55; SynTech, Kirchzarten, Germany). The signals were detected by a high-impedance amplifier (IDAC-2; SynTech, Kirchzarten, Germany) and analyzed using SynTech software (GC-EAD 2014 v1.2.5). The stimulation was started with a pulse duration of 0.2 s until the baseline stabilized. The flow rate was 4 mL/s. The EAG responses were recorded during odor stimulation, and six female or male antennae of 3-day-old adults at 48 h after injection were used for each volatile odorant recording. experiments were divided into 5 treatment groups (dsOR4-injected, dsOR11-injected, dsOR15-injected, dsORco-injected and dsGFP-injected) with 3 biological replicates per group and 20 individuals per replicate. The treated insects were reared under the above conditions. Antennae of GdauORco group were collected separately for interference efficiency measurement by RT-qPCR at 24, 48, 72, and 96 h after injection. Antennae of GdauOR4, GdauOR11, and GdauOR15 groups were collected for interference efficiency measurement at 48 h after injection. The group of dsGFP-injected was treated as control.

EAG Recordings
Electroantennogram (EAG) was used to record the antennal responses to 13 volatiles of A. mongolium (Table 1) [20]. These compounds were dissolved in mineral oil (Bio-Rad, USA) at 10 μg/μL. The mineral oil was also set up as a blank control. An antenna was cut off from the base, and cut 0.2 mm off the tip, then connected between the electrodes by two silver wires and two capillary glass tubes filled with physiological saline (0.9% sodium chloride injection). The antennal tip and base were inserted into the capillary tubes with 0.5 mm each. Filter paper strips (5 × 15 mm) were loaded with 10 μL of each solution, and transferred to the Pasteur pipette. The top of the pipette was inserted into the small hole in the wall of the tube, which was connected to an air stimulus controller (CS-55; SynTech, Kirchzarten, Germany). The signals were detected by a high-impedance amplifier (IDAC-2; SynTech, Kirchzarten, Germany) and analyzed using SynTech software (GC-EAD 2014 v1.2.5). The stimulation was started with a pulse duration of 0.2 s until the baseline stabilized. The flow rate was 4 mL/s. The EAG responses were recorded during odor stimulation, and six female or male antennae of 3-day-old adults at 48 h after injection were used for each volatile odorant recording. experiments were divided into 5 treatment groups (dsOR4-injected, dsOR11-injected, dsOR15-injected, dsORco-injected and dsGFP-injected) with 3 biological replicates per group and 20 individuals per replicate. The treated insects were reared under the above conditions. Antennae of GdauORco group were collected separately for interference efficiency measurement by RT-qPCR at 24, 48, 72, and 96 h after injection. Antennae of GdauOR4, GdauOR11, and GdauOR15 groups were collected for interference efficiency measurement at 48 h after injection. The group of dsGFP-injected was treated as control.

EAG Recordings
Electroantennogram (EAG) was used to record the antennal responses to 13 volatiles of A. mongolium (Table 1) [20]. These compounds were dissolved in mineral oil (Bio-Rad, USA) at 10 μg/μL. The mineral oil was also set up as a blank control. An antenna was cut off from the base, and cut 0.2 mm off the tip, then connected between the electrodes by two silver wires and two capillary glass tubes filled with physiological saline (0.9% sodium chloride injection). The antennal tip and base were inserted into the capillary tubes with 0.5 mm each. Filter paper strips (5 × 15 mm) were loaded with 10 μL of each solution, and transferred to the Pasteur pipette. The top of the pipette was inserted into the small hole in the wall of the tube, which was connected to an air stimulus controller (CS-55; SynTech, Kirchzarten, Germany). The signals were detected by a high-impedance amplifier (IDAC-2; SynTech, Kirchzarten, Germany) and analyzed using SynTech software (GC-EAD 2014 v1.2.5). The stimulation was started with a pulse duration of 0.2 s until the baseline stabilized. The flow rate was 4 mL/s. The EAG responses were recorded during odor stimulation, and six female or male antennae of 3-day-old adults at 48 h after injection were used for each volatile odorant recording.  experiments were divided into 5 treatment groups (dsOR4-injected, dsOR11-injected, dsOR15-injected, dsORco-injected and dsGFP-injected) with 3 biological replicates per group and 20 individuals per replicate. The treated insects were reared under the above conditions. Antennae of GdauORco group were collected separately for interference efficiency measurement by RT-qPCR at 24, 48, 72, and 96 h after injection. Antennae of GdauOR4, GdauOR11, and GdauOR15 groups were collected for interference efficiency measurement at 48 h after injection. The group of dsGFP-injected was treated as control.

EAG Recordings
Electroantennogram (EAG) was used to record the antennal responses to 13 volatiles of A. mongolium (Table 1) [20]. These compounds were dissolved in mineral oil (Bio-Rad, USA) at 10 μg/μL. The mineral oil was also set up as a blank control. An antenna was cut off from the base, and cut 0.2 mm off the tip, then connected between the electrodes by two silver wires and two capillary glass tubes filled with physiological saline (0.9% sodium chloride injection). The antennal tip and base were inserted into the capillary tubes with 0.5 mm each. Filter paper strips (5 × 15 mm) were loaded with 10 μL of each solution, and transferred to the Pasteur pipette. The top of the pipette was inserted into the small hole in the wall of the tube, which was connected to an air stimulus controller (CS-55; SynTech, Kirchzarten, Germany). The signals were detected by a high-impedance amplifier (IDAC-2; SynTech, Kirchzarten, Germany) and analyzed using SynTech software (GC-EAD 2014 v1.2.5). The stimulation was started with a pulse duration of 0.2 s until the baseline stabilized. The flow rate was 4 mL/s. The EAG responses were recorded during odor stimulation, and six female or male antennae of 3-day-old adults at 48 h after injection were used for each volatile odorant recording.  experiments were divided into 5 treatment groups (dsOR4-injected, dsOR11-injected, dsOR15-injected, dsORco-injected and dsGFP-injected) with 3 biological replicates per group and 20 individuals per replicate. The treated insects were reared under the above conditions. Antennae of GdauORco group were collected separately for interference efficiency measurement by RT-qPCR at 24, 48, 72, and 96 h after injection. Antennae of GdauOR4, GdauOR11, and GdauOR15 groups were collected for interference efficiency measurement at 48 h after injection. The group of dsGFP-injected was treated as control.

EAG Recordings
Electroantennogram (EAG) was used to record the antennal responses to 13 volatiles of A. mongolium (Table 1) [20]. These compounds were dissolved in mineral oil (Bio-Rad, USA) at 10 μg/μL. The mineral oil was also set up as a blank control. An antenna was cut off from the base, and cut 0.2 mm off the tip, then connected between the electrodes by two silver wires and two capillary glass tubes filled with physiological saline (0.9% sodium chloride injection). The antennal tip and base were inserted into the capillary tubes with 0.5 mm each. Filter paper strips (5 × 15 mm) were loaded with 10 μL of each solution, and transferred to the Pasteur pipette. The top of the pipette was inserted into the small hole in the wall of the tube, which was connected to an air stimulus controller (CS-55; SynTech, Kirchzarten, Germany). The signals were detected by a high-impedance amplifier (IDAC-2; SynTech, Kirchzarten, Germany) and analyzed using SynTech software (GC-EAD 2014 v1.2.5). The stimulation was started with a pulse duration of 0.2 s until the baseline stabilized. The flow rate was 4 mL/s. The EAG responses were recorded during odor stimulation, and six female or male antennae of 3-day-old adults at 48 h after injection were used for each volatile odorant recording.  experiments were divided into 5 treatment groups (dsOR4-injected, dsOR11-injected, dsOR15-injected, dsORco-injected and dsGFP-injected) with 3 biological replicates per group and 20 individuals per replicate. The treated insects were reared under the above conditions. Antennae of GdauORco group were collected separately for interference efficiency measurement by RT-qPCR at 24, 48, 72, and 96 h after injection. Antennae of GdauOR4, GdauOR11, and GdauOR15 groups were collected for interference efficiency measurement at 48 h after injection. The group of dsGFP-injected was treated as control.

EAG Recordings
Electroantennogram (EAG) was used to record the antennal responses to 13 volatiles of A. mongolium (Table 1) [20]. These compounds were dissolved in mineral oil (Bio-Rad, USA) at 10 μg/μL. The mineral oil was also set up as a blank control. An antenna was cut off from the base, and cut 0.2 mm off the tip, then connected between the electrodes by two silver wires and two capillary glass tubes filled with physiological saline (0.9% sodium chloride injection). The antennal tip and base were inserted into the capillary tubes with 0.5 mm each. Filter paper strips (5 × 15 mm) were loaded with 10 μL of each solution, and transferred to the Pasteur pipette. The top of the pipette was inserted into the small hole in the wall of the tube, which was connected to an air stimulus controller (CS-55; SynTech, Kirchzarten, Germany). The signals were detected by a high-impedance amplifier (IDAC-2; SynTech, Kirchzarten, Germany) and analyzed using SynTech software (GC-EAD 2014 v1.2.5). The stimulation was started with a pulse duration of 0.2 s until the baseline stabilized. The flow rate was 4 mL/s. The EAG responses were recorded during odor stimulation, and six female or male antennae of 3-day-old adults at 48 h after injection were used for each volatile odorant recording.  experiments were divided into 5 treatment groups (dsOR4-injected, dsOR11-injected, dsOR15-injected, dsORco-injected and dsGFP-injected) with 3 biological replicates per group and 20 individuals per replicate. The treated insects were reared under the above conditions. Antennae of GdauORco group were collected separately for interference efficiency measurement by RT-qPCR at 24, 48, 72, and 96 h after injection. Antennae of GdauOR4, GdauOR11, and GdauOR15 groups were collected for interference efficiency measurement at 48 h after injection. The group of dsGFP-injected was treated as control.

EAG Recordings
Electroantennogram (EAG) was used to record the antennal responses to 13 volatiles of A. mongolium (Table 1) [20]. These compounds were dissolved in mineral oil (Bio-Rad, USA) at 10 μg/μL. The mineral oil was also set up as a blank control. An antenna was cut off from the base, and cut 0.2 mm off the tip, then connected between the electrodes by two silver wires and two capillary glass tubes filled with physiological saline (0.9% sodium chloride injection). The antennal tip and base were inserted into the capillary tubes with 0.5 mm each. Filter paper strips (5 × 15 mm) were loaded with 10 μL of each solution, and transferred to the Pasteur pipette. The top of the pipette was inserted into the small hole in the wall of the tube, which was connected to an air stimulus controller (CS-55; SynTech, Kirchzarten, Germany). The signals were detected by a high-impedance amplifier (IDAC-2; SynTech, Kirchzarten, Germany) and analyzed using SynTech software (GC-EAD 2014 v1.2.5). The stimulation was started with a pulse duration of 0.2 s until the baseline stabilized. The flow rate was 4 mL/s. The EAG responses were recorded during odor stimulation, and six female or male antennae of 3-day-old adults at 48 h after injection were used for each volatile odorant recording.

Data Analysis
All data were statistically analyzed using SPSS Statistics 18.0. In the expression profile analysis, the differences between males and females were compared using T-test. Duncan's test was performed for the differences between different treatments using One-Way ANOVA, and p < 0.05 was taken as the level of statistical significance. RNA interference and EAG analysis were analyzed using T-test.

Identification of Putative OR Genes in G. daurica
A total of 21 putative OR genes were identified from the transcriptome database of G. daurica adults, which were named as GdauOR1-20 and GdauORco. These sequences data have been submitted to the GenBank Datebase of NCBI under accession numbers MK691770-MK691790 (Table 2).

Data Analysis
All data were statistically analyzed using SPSS Statistics 18.0. In the expression profile analysis, the differences between males and females were compared using T-test. Duncan's test was performed for the differences between different treatments using One-Way ANOVA, and p < 0.05 was taken as the level of statistical significance. RNA interference and EAG analysis were analyzed using T-test.

Identification of Putative OR Genes in G. daurica
A total of 21 putative OR genes were identified from the transcriptome database of G. daurica adults, which were named as GdauOR1-20 and GdauORco. These sequences data have been submitted to the GenBank Datebase of NCBI under accession numbers MK691770-MK691790 (Table 2).

Data Analysis
All data were statistically analyzed using SPSS Statistics 18.0. In the expression profile analysis, the differences between males and females were compared using T-test. Duncan's test was performed for the differences between different treatments using One-Way ANOVA, and p < 0.05 was taken as the level of statistical significance. RNA interference and EAG analysis were analyzed using T-test.

Identification of Putative OR Genes in G. daurica
A total of 21 putative OR genes were identified from the transcriptome database of G. daurica adults, which were named as GdauOR1-20 and GdauORco. These sequences data have been submitted to the GenBank Datebase of NCBI under accession numbers MK691770-MK691790 (Table 2).

Data Analysis
All data were statistically analyzed using SPSS Statistics 18.0. In the expression profile analysis, the differences between males and females were compared using T-test. Duncan's test was performed for the differences between different treatments using One-Way ANOVA, and p < 0.05 was taken as the level of statistical significance. RNA interference and EAG analysis were analyzed using T-test.

Identification of Putative OR Genes in G. daurica
A total of 21 putative OR genes were identified from the transcriptome database of G. daurica adults, which were named as GdauOR1-20 and GdauORco. These sequences data have been submitted to the GenBank Datebase of NCBI under accession numbers MK691770-MK691790 (Table 2).

Data Analysis
All data were statistically analyzed using SPSS Statistics 18.0. In the expression profile analysis, the differences between males and females were compared using T-test. Duncan's test was performed for the differences between different treatments using One-Way ANOVA, and p < 0.05 was taken as the level of statistical significance. RNA interference and EAG analysis were analyzed using T-test.

Identification of Putative OR Genes in G. daurica
A total of 21 putative OR genes were identified from the transcriptome database of G. daurica adults, which were named as GdauOR1-20 and GdauORco. These sequences data have been submitted to the GenBank Datebase of NCBI under accession numbers MK691770-MK691790 (Table 2).

Data Analysis
All data were statistically analyzed using SPSS Statistics 18.0. In the expression profile analysis, the differences between males and females were compared using T-test. Duncan's test was performed for the differences between different treatments using One-Way ANOVA, and p < 0.05 was taken as the level of statistical significance. RNA interference and EAG analysis were analyzed using T-test.

Identification of Putative OR Genes in G. daurica
A total of 21 putative OR genes were identified from the transcriptome database of G. daurica adults, which were named as GdauOR1-20 and GdauORco. These sequences data have been submitted to the GenBank Datebase of NCBI under accession numbers MK691770-MK691790 (Table 2).

Expression Profile Analysis of G. daurica OR Genes
The expression profiles of GdauORs in different tissues (Figure 1) showed that the expression levels of fifteen genes in antennae were significantly higher than those in other tissues, including GdauOR1, GdauOR3-6, GdauOR8, GdauOR11-13, GdauOR15, GdauOR17-20, and GdauORco. Among them, the expression levels of seven genes (GdauOR1, GdauOR5-6, GdauOR17, GdauOR19-20, and GdauORco) in female antennae were significantly higher than those in male antennae, while four genes (GdauOR3, GdauOR8, GdauOR11, and GdauOR12) were opposite. Significantly, the genes with the highest expression levels in antennae were GdauOR4 and GdauORco. In addition, the expressions of GdauOR3, GdauOR7, GdauOR16, GdauOR18 and GdauOR20 were significantly enriched in wings. Peculiarly, GdauOR2 was abundantly expressed in abdomens, and the expression level in males was significantly higher than that in females. While GdauOR9, GdauOR10, and GdauOR14 were not shown in the graph, because their expression levels were too low.   The expression profiling results of GdauORs at different days after eclosion ( Figure 2) showed that GdauOR1 had a higher expression level in males at 15 d, while GdauOR2 had a higher expression level in females at 15 d. GdauOR3 had a higher expression level in females at 45 d. GdauOR4 and GdauORco had a higher expression level in females at 90 d. GdauOR5 had a higher expression level in females at 25 d. GdauOR6 had a higher expression level in females at 7 d and males in 25 d as well as 45 d. GdauOR7 and GdauOR13 had a higher expression level in males at 1 d. GdauOR8 and GdauOR11 had a higher expression level in males at 90 d. GdauOR12 had a higher expression level in males at 7 d. The expression level of GdauOR15 increased with the days after eclosion, and highest in females at 90 d. GdauOR16 had a higher expression level in males at 45 d. GdauOR17 had a higher expression level in males at 1 d and 7 d. GdauOR18 had a higher expression level in males at 25 d. GdauOR19 had almost the same expression levels for all stages of adults with few exceptions. GdauOR20 had a higher expression level in females at 7 d. Each gene showed significant differences between male and female antennae at different days of eclosion. It is noteworthy that there were significant differences in expression levels of 13 genes between male and female antennae at 90 d, when they were in the mating stage. Among them, the expression levels of GdauOR4, GdauOR7, GdauOR13, GdauOR15, GdauOR16, GdauOR17, and GdauORco were significantly higher in females than in males, whereas GdauOR1, GdauOR3, GdauOR5, GdauOR6, GdauOR8, and GdauOR11 were opposite.

RNA Interference of GdauOR4/GdauOR11/GdauOR15/GdauORco
To examine the RNAi efficiency, RT-qPCR was performed. Briefly, injection of dsRNA significantly decreased expression levels of target genes. Compared with the dsGFPinjected, when dsORco were injected, the mRNA levels of GdauORco in female antennae were decreased by 56%, 76%, 89%, and 87% at 24, 48, 72, and 96 h post-injection, respectively ( Figure 3). The mRNA levels of GdauORco in male antennae were decreased by 45%, 80%, 83%, and 88% at 24, 48, 72, and 96 h post-injection, respectively. This result indicated that a reduction of RNAi efficiency to about 25% at 48 h post-injection is well processed for GdauORs. Similarly, the expression levels of dsOR4-injected, dsOR11-injected, and dsOR15injected at 48 h post-injection in both female and male antennae were reduced to less than 27% compared with control groups (Figure 4). Among them, RNAi reduced the expression levels of GdauOR4 to 19.9% and 27% in females and males, GdauOR11 was reduced to 14.5% and 19.5% in females and males, and GdauOR15 was reduced to 20.6% and 25.8% in females and males, respectively. There was no significant influence on expression levels of GdauOR11 and GdauOR15 when dsOR4 was injected (Figure 4a), and similarly, injection of dsOR11 did not affect the expression of GdauOR4 and GdauOR15 (Figure 4b); injection of dsOR15 also did not affect the expression of GdauOR4 and GdauOR11 (Figure 4c).  45-day-old; 90 d: 90-day-old. Error bars represent the standard error of three independent experiments. Different capital and small letters above bars indicate significant difference among different developmental days of female and male, respectively (Duncan's test; p < 0.05). The asterisk above bars indicates significant difference between males and females (T-test; *: p < 0.05; **: p < 0.01; ns: No significant difference).

Electroantennogram Analysis
In order to elucidate the physiological function of GdauOR4, GdauOR11, GdauOR15 and GdauORco in the perception of host plant volatiles, the electrophysiological responses after RNAi were detected by EAG experiments ( Figure 5). Compared with the control group, EAG activity of females in response to six volatiles was significantly reduced when dsOR4 was injected, including diallyl sulfide, diallyl disulphide, diallyl trisulfide, 2hexenal, disulfide methyl 2-propenyl (p < 0.01), and dimethyl trisulfide (p < 0.05). Males showed significantly lower EAG responses to four odors when dsOR4 was injected, including diallyl sulfide, diallyl disulphide, diallyl trisulfide, and 2-hexenal (p < 0.01). In the dsOR11-injected group, both females and males had no significant influence on the electrophysiological responses to all tested odors. In dsOR15-injected group, EAG activity of females in response to three volatiles was significantly reduced, including dimethyl trisulfide, myrcene, and disulfide methyl 2-propenyl (p < 0.01). However, there was no significant influence on the electrophysiological responses when dsOR15 was injected into males. When ORco was silenced, the electrophysiological responses to eight volatiles were significantly decreased in females, including diallyl sulfide, dimethyl trisulfide, diallyl disulphide, 2-hexen-1-ol, 2-hexenal, 1,3,5-cycloheptatriene, disulfide methyl 2-propenyl 45-day-old; 90 d: 90-day-old. Error bars represent the standard error of three independent experiments. Different capital and small letters above bars indicate significant difference among different developmental days of female and male, respectively (Duncan's test; p < 0.05). The asterisk above bars indicates significant difference between males and females (T-test; *: p < 0.05; **: p < 0.01; ns: No significant difference).

Electroantennogram Analysis
In order to elucidate the physiological function of GdauOR4, GdauOR11, GdauOR15 and GdauORco in the perception of host plant volatiles, the electrophysiological responses after RNAi were detected by EAG experiments ( Figure 5). Compared with the control group, EAG activity of females in response to six volatiles was significantly reduced when dsOR4 was injected, including diallyl sulfide, diallyl disulphide, diallyl trisulfide, 2hexenal, disulfide methyl 2-propenyl (p < 0.01), and dimethyl trisulfide (p < 0.05). Males showed significantly lower EAG responses to four odors when dsOR4 was injected, including diallyl sulfide, diallyl disulphide, diallyl trisulfide, and 2-hexenal (p < 0.01). In the dsOR11-injected group, both females and males had no significant influence on the electrophysiological responses to all tested odors. In dsOR15-injected group, EAG activity of females in response to three volatiles was significantly reduced, including dimethyl trisulfide, myrcene, and disulfide methyl 2-propenyl (p < 0.01). However, there was no significant influence on the electrophysiological responses when dsOR15 was injected into males. When ORco was silenced, the electrophysiological responses to eight volatiles were significantly decreased in females, including diallyl sulfide, dimethyl trisulfide, diallyl disulphide, 2-hexen-1-ol, 2-hexenal, 1,3,5-cycloheptatriene, disulfide methyl 2-propenyl

Discussion
In this study, we identified 21 OR genes from the adult G. daurica transcriptome, less than identified for most other Coleoptera [28][29][30]. We speculate that our use of a highthroughput sequencing approach using whole insects rather than antennae may result in a significant number of olfactory genes being buried. Expression profiling of chemosensory genes are of primary importance for exploring the function and the olfactory recognition mechanism of insects [31][32][33]. The tissue expression profiles of 21 ORs showed that most OR genes (GdauOR1, GdauOR3-6, GdauOR8, GdauOR11-13, GdauOR15, GdauOR17-20, and GdauORco) were highly expressed in antennae, which were similar to the results in other species [9]. This result suggested that these ORs might be involved in the chemosensory process [34][35][36]. Notably, GdauORco and GdauOR4 were extremely highly expressed in antennae, suggesting that they may play an important role in the olfactory system. While, GdauOR9, GdauOR10, and GdauOR14 had extremely low expression, presumably expressed in the other developmental stages [37,38]. In addition, the expressions of GdauOR3, GdauOR7, GdauOR16, GdauOR18, and GdauOR20 were significantly enriched in wings, and GdauOR2 were abundantly expressed in abdomens, which suggested that they might be involved in other physiological functions of G. daurica. Different days after eclosion expression analysis indicated that the expression levels of ORs in antennae were affected by age. The similar results were observed in other olfactory proteins [39]. The expression profile analysis lays a necessary foundation for revealing the olfactory recognition mechanism of G. daurica. Four ORs, GdauOR4, GdauOR11, GdauOR15, and GdauORco were selected for further research in odor detection, because they were highly expressed in antennae.
In a previous study, RNAi technology was effectively used in G. daurica by dsRNA injection [23]. Thus, RNAi injection experiments against four ORs were conducted. Same as the previous study, the expression levels were significantly decreased 48 h after dsRNA injection, indicating that dsRNA injection was suitable for target gene interference in G. daurica. Furthermore, it has been reported that dsRNA injection can lead to off-target, which means that the expression levels of other non-target genes may be reduced by RNAi [40]. As far as our current work is concerned, there were no significant influence on expression levels of GdauOR11 and GdauOR15 when dsOR4 was injected, and similarly, injection of dsOR11 did not affected the expression of GdauOR4 and GdauOR15, injection of dsOR15 also did not affected the expression of GdauOR4 and GdauOR11. This RT-qPCR result is consistent with previous study that the mRNA levels of HarmPBP1, HarmPBP2 and HarmPBP3 were not affected by each specific dsRNA injection in Helicoverpa armigera [41].
Many studies have shown that the involvement of genes in olfactory functions can be ultimately impaired by silencing individual OR genes to influence odor preference [35,42]. The silencing of a single OR gene by RNAi resulted in different electrophysiological changes of G. daurica to host plant volatiles. The electrophysiological response of females to six volatiles was significantly reduced when dsOR4 was injected, including diallyl sulfide, diallyl disulphide, diallyl trisulfide, 2-hexenal, disulfide methyl 2-propenyl, and dimethyl trisulfide. Males showed significantly lower EAG responses to four odors when dsOR4 were injected, including diallyl sulfide, diallyl disulphide, diallyl trisulfide, and 2-hexenal. Among them, diallyl sulfide, diallyl disulphide, diallyl trisulfide, disulfide methyl 2propenyl, and dimethyl trisulfide are sulfocompounds with a strong pungent odor, which are the symbolic components of Allium plants in Liliaceae such as onion and garlic [43][44][45]. In addition, diallyl disulphide and disulfide methyl 2-propenyl account for 43.16% of the total volatile of A. mongolicum [20]. Combined with the expression profile analysis, the expression level of GdauOR4 was the highest among all ORs. Thus, we speculate that GdauOR4 may be the pivotal receptor for host location of G. daurica. For instance, GmolOR9 was highly expressed in antennae, and silencing GmolOR9 resulted in reduced sensitivity of Grapholita molesta to host volatiles [46].
There was no significant influence on the electrophysiological responses of both females and male to host volatiles after silencing GdauOR11. It implied that it may be involved in the recognition of other semiochemicals. In addition, the expression level of GdauOR11 in male antennae was significantly higher than in females during the whole developmental period of adults. Based on these findings, we speculate that GdauOR11 may be involved in recognizing pheromones released by females. In Mythimna separata, MsepOR3 was specifically and abundantly expressed in male antennae, and the Xenopus oocytes expressing MsepOR3/ORco evinced dose dependent responses to the sex pheromone Z11-16: Ald [47].
EAG activity of females in response to three host plant volatiles (dimethyl trisulfide, myrcene, and disulfide methyl 2-propenyl) was significantly reduced after silencing GdauOR15. However, there was no significant influence on the electrophysiological responses in males. Expression profile analysis showed that GdauOR15 were mainly expressed in antennae. Moreover, the expression level increased gradually with the development of adults, and it was significantly higher in females than in females at 45-days-old and 90-days-old. This implies that GdauOR15 may have different functions in odor perception between males and females in G. daurica, and it may play more important roles in females searching for host plants or suitable oviposition sites. For instance, HparOR27 was expressed mainly in female antennae, and it was broadly responsive to three host plant volatiles in Holotrichia parallela [48].
When GdauORco was silenced, the electrophysiological responses to eight host plant volatiles were significantly decreased in females. Meanwhile, the EAG responses of males to six volatiles were reduced significantly. These results suggest that GdauORco is necessary for odorant responses of G. daurica. ORs cannot assemble, transport, or function in the absence of ORco and the loss of this single receptor results in dramatically impaired olfactory behavior [4]. For instance, the HEK293/MdesOR115 cell line which was devoid of ORco did not respond to compounds in Hessian fly pheromone reception [11]. RNAi reduced the expression of SaveOrco to 34.11% in aphids, resulting in weaker EAG responses to plant volatiles and aphid alarm pheromone [49]. Taken together, interference with GdauORco resulted in EAG response decreased to eight volatiles in females and six volatiles in males, respectively, these substances include those causing a reduction in EAG response after silencing other three GdauOR genes, except myrcene. It indicates that GdauORco functionally synergizes with GdauOR4, GdauOR11, and GdauOR15 in the olfactory recognition process of G. daurica. The response to myrcene may be due to incomplete RNAi or to other unidentified targets OR [50,51].
Although the silencing of a single OR gene by RNAi resulted in electrophysiological changes of G. daurica to host plant volatiles, whether it has an impact on behavior needs to be further verified in future experiments.

Conclusions
Twenty-one OR genes were identified from the adult G. daurica transcriptome. Most ORs were abundantly expressed in antennae and the expression levels of ORs in adults were affected by age. Based on RNAi and EAG experiments, we speculate that GdauOR4 may be a pivotal receptor for host location in G. daurica. GdauOR15 may play more important roles in females, but this needs further study for confirmation.

Supplementary Materials:
The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/insects13070563/s1, Table S1: List of primers used for qRT-PCR; Table S2: List of primers used in the dsRNA synthesis.