Evidence of the Involvement of a Plus-C Odorant-Binding Protein HparOBP14 in Host Plant Selection and Oviposition of the Scarab Beetle Holotrichia parallela

Simple Summary The scarab beetle Holotrichia parallela is a serious underground pest and causes serious damages in China to a variety of crops. To reduce the use of pesticides, insect olfactory proteins attract more and more attention in the development of pollution-free control agents in plant protection. In this study, we evaluate the molecular mechanism in the scarab beetle to detect oviposition cues. We clone a leg biased gene HparOBP14 which encodes for an odorant-binding protein of the scarab beetle and demonstrate its involvement in binding, electrophysiological, and behavioral responses to the oviposition chemicals by the knockdown of HparOBP14 expression using RNA interference technique. Our study provides a strong theoretical basis for the development of environmentally acceptable strategies for H. parallela control. Abstract Holotrichia parallela is one of the agriculturally important scarab beetle pests in China. In this study, HparOBP14 was cloned, which is the most abundantly expressed among the OBP genes in the legs of female H. parallela adults. Sequence comparison and phylogenetic analysis showed that HparOBP14 has a Plus-C structure motif. The expression profile analysis revealed that HparOBP14 expression was the highest in the female antennae and then in the legs. The fluorescence competitive binding experiment of the recombinant HparOBP14 protein showed that HparOBP14 had an affinity with 6-methyl-5-heptene-2-one (plant volatile), 3-methylindole, p-cymene, methanol, formaldehyde, α-pinene, and geraniol (organic fertilizer volatile). Knockdown HparOBP14 expression decreased significantly the EAG response of the injected female adults to p-cymene, methanol, formaldehyde, α-pinene, and geraniol. Similarly, the injected female adults were significantly less attracted to geraniol and methanol. Therefore, HparOBP14 might bind organic matter volatiles during oviposition. These results are not only helpful to analyze the olfactory recognition mechanism of female adult H. parallela when choosing suitable oviposition sites, but also to provide target genes for green prevention and control of H. parallela in the future.


Introduction
The scarab beetle Holotrichia parallela Motschulsk (Coleoptera: Scarabaeida e) is a seriously harmful pest in agriculture and horticulture in China. The beetles principally damage the leaves of plants such as elm, willow, and apple tree. The larvae eat the roots of

Tissue-Specific RNA Extraction and cDNA Synthesis
The tissues (antennae, heads without antennae, thoraxes, abdomens, legs, and wings) of 9-11-day-old female beetles were collected separately into 1.5 mL centrifuge tubes, frozen quickly in liquid nitrogen, and then stored at −80 • C for further study. The total RNA was extracted using a total RNA extraction kit (Tianmo, Beijing, China) according to the manufacturer's instructions. The yield and quality of RNA were determined (A260/A280 = 1.8~2.0) with NanoDrop 2000c spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) and 1% agarose gel electrophoresis. First-strand cDNA was synthesized with FastKing gDNA Dispelling RT SuperMix (TIANGEN, Beijing, China). The reverse transcription reaction contains: 4 µL 5 × FastKing-RT SuperMix, 1 µL total RNA (100 ng/µL), 15 µL RNase-Free ddH 2 O. The condition of the reverse transcription reactions was as follows: 42 • C for 15 min (reverse transcription reaction), followed by 95 • C for 3 min (inactive enzymes). All cDNA samples were stored at −20 • C.

Phylogenetic Analysis and Expression Levels of OBP Genes
Twenty OBP genes were identified from results of the transcriptome of female H. parallela legs [35]. The gene-specific primers of HparOBP genes were designed by Primer Premier 5 software (PREMIER Biosoft, San Francisco, CA, USA) and are listed in Table  S1. The specificity of each primer set was validated by melting curve analysis, and the efficiency was calculated by analyzing the standard curves with a fivefold cDNA dilution series. The expression level of all 20 HparOBPs in the legs of H. parallela female adults was analyzed, and the expression of HparOBP14 in each tissue (antennae, heads without antennae, thoraxes, abdomens, legs, and wings) was analyzed using RT-qPCR with the gene-specific primers using the 7500 Real-Time PCR Detection System (ABI, Singapore). The first strand cDNAs synthesized by the reverse transcription in 2.2 described above were used as templates. The GAPDH (glyceraldehyde-3-phosphate dehydrogenase) was used as an internal reference gene for normalizing and calculating relative HparOBPs expression. According to the operation of GoTaq ® qPCR Master Mix Kit (Promega, Wisconsin, US), the RT-qPCR reaction was prepared on ice with 2 µL of cDNA template, 10 µL of the Master Mix, 0.5 µL of each of the primer pairs (10 µmol/L), and 7 µL of ddH2O. The RT-qPCR reaction was carried out by 95 • C desaturation for 2 min, followed by 40 cycles of 95 • C for 30 s, and 60 • C for 1 min. Each reaction was performed with three biological samples as biological replicates and repeated three times as technical replicates. The data of the obtained Ct values was analyzed by using 2 −∆∆CT method [36].

Molecular Cloning of HparOBP14
Primer Premier 5.0 was utilized to design specific primers with BamHI and HindIII restriction sites (Table S1). The RNA was extracted from female antennae and used to synthesize cDNA for molecular cloning. The full-length sequence of HparOBP14 was amplified using 2 × Taq PCR Master Mix kit (TIANGEN Biotech, Beijing, China). The PCR reaction mixture contained 1 µL cDNA template, 25 µL 2 × PCR Master Mix, 1.5 µL each of forward and reverse primers (10 µmol/L), and 21 µL H2O. The PCR reaction was carried out at 95 • C for 3 min, followed by 35 cycles of 94 • C for 30 s, 60 • C for 30 s, 72 • C for 30 s, and final 72 • C extension for 10 min. The amplified products were detected by 1% agarose gel electrophoresis and purified according to the instructions of the TIANel Midi Purification Kit (TIANGEN Biotech, Beijing, China). After digestion with BamHI and HindIII, they were ligated into pGEM-T4 Easy Vector at 4 • C overnight to make pGEM-T4-OBP14 constructs. The ligation reaction mixture (5 µL) was transformed into Trans5α chemically competent cell (50 µL) and plated on agar plates (contains a final concentration of 50 µg/mL kanamycin). The positive clones (white monoclone) were randomly selected for sequencing. The sequence obtained was analyzed by BLAST search with the National Center for Biotechnology Information (NCBI) BLAST server in 10 May 2019 (http://www.ncbi.nlm.nih.gov/). The amino acid sequence was determined from the longest open reading frame found by using NCBI ORF finder in 10 May 2019 (https://www.ncbi.nlm.nih.gov/orffinder/) and compared by using DNAMAN. Signal peptide was predicted using Signalp-5.0 Server in 10 May 2019 (http://www.cbs.dtu.dk/ services/SignalP/index.php).

Expression and Purification of Recombinant Protein
The pGEM-T4-OBP14 construct and pET-30a(+) plasmid were digested with BamHI and HindIII at 37 • C for 4~5 h. The target fragments were gel purified and ligated with T4 DNA ligase (Takara Bio, Beijing, China) to obtain the expression construct pET30a-OBP14. The pET30a-OBP14 construct was transformed into BL21(DE3) cells (TIANGEN, Beijing, China), and the positive clone was picked and cultured in 5 mL LB medium containing kanamycin (final concentration, 50 µg/mL) overnight at 37 • C with vigorous shaking. The bacteria solution was sent out for sequencing (Shenzhen Huada Genomics Co., Ltd.) and an aliquot was also stored at −20 • C. The sequenced bacteria solution that was shown to contain the correct HparOBP14 sequence was inoculated into 500 mL LB liquid medium (supplemented with final concentration 50 µg/mL kanamycin) at a ratio of 1:100 and cultivated at 37 • C with 220 r/min vigorous shaking until the OD 600 value reached 0.8~1.0, then isopropyl-β-D-1-thiogalactopyranoside (IPTG) was added at a final concentration of 0.25 mmol/L. The expression of HparOBP14 was induced for 8 h at 28 • C with 180 r/min shaking. SDS-PAGE was used to monitor the expression of the recombinant HparOBP14 protein. The protein expressed in the supernatant was purified by Ni ion affinity chromatography (ÄKTA™ Avant Cytiva, GE, USA). The eluent containing HparOBP14 protein was dialyzed overnight in a Tris-HCL buffer (50 mmol/L, pH 7.4). The His-tag of the purified recombinant protein was removed by using recombinant enterokinase (Shanghai Korain Biotech Co Ltd., Shanghai, China), then a second purification was performed with the Ni affinity chromatography. The protein concentration and purity were determined with FlexStation3 microplate reader (Molecular Devices Inc.,CA, USA), and the pure protein samples were stored at −80 • C for later use.

Fluorescence Displacement Binding Assay
To determine the affinity of HparOBP14 to odorant volatiles, the fluorescent displacement binding assay was carried out [36]. A total of 27 compounds was selected, including 7 compounds from host plant (elm) volatiles [38,39], 11 compounds from organic fertilizer volatiles [17], and 9 compounds common in both of them. The odorant compounds (purity ≥95%) were purchased from Sigma-Aldrich. The fluorescent probe N-phenyl-1naphthylamine (1-NPN) was dissolved in methanol (chromatographic grade) to a final concentration of 1 mM and used as the stock solution. The stock solution of the purified protein was diluted with Tris-HCl (50 mmol/L, pH 7.4) to a final concentration of 2 µM and added into a colorimetric plate cuvette. The 1-NPN stock solution was then successively added to the protein solution to the final concentrations ranging from 2 µM to 20 µM. The fluorescence changes were recorded to get K 1-NPN , the dissociation constant of the complex HparOBP14/1-NPN.
The binding ability of HparOBP14 to the odorant compounds was determined by titrating the mixture of the protein and the 1-NPN (final concentration 2 µM of both) with each odorant compound at final concentrations ranging from 2 µM to 40 µM. The competitive dissociation constant was calculated from the corresponding IC 50  The fluorescence changes during the binding of either 1-NPN or odorant compounds to the HparOBP14 protein were measured and recorded with the excitation wavelength of 337 nm and within a wavelength range of 350~550 nm by a 970CRT fluorescence spectrophotometer (Yidian Analytical Instrument Company, Shanghai, China).

RNAi-Mediated Gene Silencing
The dsRNAs of HparOBP14 and GFP (GenBank number: ACY56286) (for primers, see Table S1) were synthesized using a T7 Ribo MAX™ Express RNAi System (Promega, Madison, WI, USA) according to the manufacturer's instructions. The GFP-dsRNA was used as the negative control. Nine-to-eleven-day-old female adult beetles were injected with 2 µg dsRNA into the adult head with a 5 µL injection needle (7000 series syringes, Hamilton, USA), then kept under the same condition as 2.1 described above. The antennae were dissected at 24 h, 48 h, 72 h, and 96 h after injection, and the samples were stored at −80 • C for further RT-qPCR analyses.

Electroantennogram (EAG) Recording
The EAG recording was conducted on the adult antennae with the EAG micromanipulator MP-15 and analyzed with EAG-Pro (Syntech, Kirchzarten, Germany). Briefly, two ends of the antennae were set on electrode holders with a conductive adhesive (Spectra 360 electrode gel). An ashless filter paper (Hangzhou Special Paper Co., Ltd., Hangzhou, China) was cut into strips (3 cm × 0.5 cm). A 20 µL compound solution (concentration 1 µg/µL in paraffin oil) was applied onto the paper strip which was then placed in a Pasteur tube (17 cm × 12 cm) and applied to the antennae by purified air at a flow rate of 1 L/min. The paraffin oil (20 µL) was used as a control. The stimulus recording time was 0.5 s and the interval between two stimulations was more than 1 min.

Y-Tube Olfactometer Behavioral Bioassays
The Y-tube olfactometer mainly consisted of a Y-shaped glass tube (main arm, 30 cm; sidearm, 20 cm; inner diameter, 2.5 cm; side arm angle, 60 degrees), an air sampler, a drying tower, a gas scrubber, and a flow meter. Two gas-washing cylinders were successively connected with the flow meter, drying tower, air pump, and Y-shaped glass tube through rubber pipes. The testing compound (concentration 1 µg/µL in paraffin oil) and paraffin oil (as control) were taken and dropped respectively into a rubber septum (Pherobio Technology Co., Ltd., Beijing, China) and placed respectively into the side arms of the Y-tube. The flow rate of the atmospheric air sampler was adjusted to 500 mL/min.
A female adult beetle was introduced into the main arm at the downwind end of the Y-olfactometer. Each beetle was given 5 min to respond to the treatment, and the first choice that the beetle made for each of the lateral arms was recorded. The response was regarded as valid only if the beetle crawled more than 3 cm into a lateral arm and stayed for more than 1 min. The following measurements were recorded for all individual beetles: the number of beetles that crawled into a lateral arm of the Y-tube (responsive beetles), and the number of beetles that did not make any choice and stayed in the main arm for 5 min (non-responsive beetles).
For each compound, more than 90 beetles were tested individually until 90 responsive beetles were obtained. Each individual beetle was tested only once, and the Y-olfactometer was replaced after every 10 beetles were tested to eliminate the odor effect. The number of individual beetles that selected the compound-treated arm of the Y-olfactometer was statistically compared with the number of beetles that selected the solvent-treated arm using a Chi-square goodness of fit test.

Statistical Analysis
The data of HparOBP14 expression in different tissues was analyzed by one-way analysis of variance (ANOVA) followed by Tukey's multiple comparison test for the comparison of the means with the SPSS Statistics version 23 (IBM Corp., Armonk, NY, USA). The data of behavioral, electrophysiological experiments, and RNAi treatment were analyzed by independent-sample t-test. All figures were generated with Prism 8 (GraphPad, La Jolla, CA, USA).

Phylogenetic and Expression Level Analysis of 20 HparOBP Genes
RT-qPCR method was used to verify the differential expression of 20 HparOBP genes in the legs of female H. parallela adults. Relative to the expression of HparOBP1, HparOBP14 was the most highly expressed OBP gene in the legs, and the expression of six other genes (HparOBP6, HparOBP12, HparOBP13, HparOBP16, and HparOBP29) was also relatively high ( Figure 1). Therefore, HparOBP14, the most highly expressed OBP gene in the legs, was selected for subsequent studies. The mature protein sequences of 20 HparOBP genes were compared with 46 OBPs of Coleoptera in the phylogenetic tree constructed using MEGA 10.0 software by maximum likelihood (ML) method with 1000 bootstrap replicates ( Figure 2). In the phylogenetic analysis, HparOBP14 with HparOBP6, HparOBP16, HparOBP26, HparOBP39, AcorOBP11, AcorOBP13, DponOBP2, DhelOBP10, AchiOBP3, and CbowOBP25 all belong to Plus-C OBPs (Figure 2, pink area). HparOBP14 is clustered with AcorOBP11and HparOBP39 in the same branch, indicating that they are closest to each other. HparOBP27, HparOBP29, HoblOBP23, HoblOBP24, HoblOBP26, AcorOBP5, and CbowOBP4 all belong to Minus-C OBPs (Figure 2, yellow area). The analysis shows that 66 OBPs have high diversity, and there are only a few obvious homologous groups in the same species.

Molecular Cloning HparOBP14
The full-length cDNA encoding HparOBP14 (GenBank: KR733560.1) [40] was cloned and verified by sequencing. The HparOBP14 sequence contains a single open reading frame (ORF) of 627 nucleotides and encodes a protein of 208 amino acids with a signal peptide of 17 amino acids at the N-terminal ( Figure S1). The predicted molecular weight is 23.8 kDa. HparOBP14 and Anomala corpulenta OBP11 (GenBank: AKC58532.1) has the highest similarity of 49.06%, while it has a high similarity with H. parallela OBP16 (GenBank: AKI84374.1) and Anoplophora chinensis OBP3 (GenBank: AUF7295 0.1), with a similarity of 32.29% and 29.44%, respectively ( Figure S2). HparOBP14 has the characteristics of the Plus-C OBP subgroup of the insect OBP family. In addition to the six highly conserved cysteine residues in typical Classical OBPs, HparOBP14 has two conserved cysteine residues in the region before C1. Between C4 and C5, and after C6, there are two additional conserved cysteine, C4a and C6a. A conserved proline is located downstream of C6 ( Figure S2). The conserved motif was of HparOBP14 was C1-X24-C2-X3-C3-X43-C4-X12-C4a-X9-C5-X8-C6-P-X9-C6a-X21 (X is any amino acid followed by the numbers of amino acids) [41].

Tissue Expression Profile Analysis of HparOBP14
The expression level of HparOBP14 in different tissues of H. parallela was analyzed by RT-qPCR. HparOBP14 had the highest expression in the antennae, followed by an intermediate expression in the legs and little expression in other tissues (Figure 3). The specific expression of HparOBP14 in the antennae and legs of female adult suggests that it may play a role in some physiological functions involving antennae and legs of H. parallela.

Transcription Level of HparOBP14 after dsRNA Injection
RNAi by dsRNA injection is an important means to examine gene functions. In this study, the expression of HparOBP14 mRNA was detected at 24 h, 48 h, 72 h, and 96 h after dsHparOBP14 injection. The relative expressions of HparOBP14 in the female antennae were decreased to 20% (24 h), 30% (48 h), 56% (72 h), and 70% (96 h) by the dsHparOBP14 injection, which are significantly different from the GFP gene expressions by the dsGFP injection (p < 0.05) ( Figure 6). There was no significant change in the gene expression levels of HparOBP14 in the antennae between uninjected and dsGFP-injected insects (p > 0.05) ( Figure 6). The results (mean ± SEM) were evaluated using a 2 −∆∆CT method, and relative to the expression level of HparOBP14 in the female abdomen. The comparison was analyzed by independent-sample t-test (n = 3). The significant difference was presented by the stars on the top of the error bars (* for p < 0.05, ** for p < 0.01, *** for p < 0.001 and ns for no significant difference).

Electroantennogram Responses of Female H. parallela after dsRNA Injection
The female H. parallela adults injected with dsHparOBP14 for 48 h were used to measure their responses to a total of seven compounds (3-methylindole, p-cymene, α-pinene, geraniol, methanol, formaldehyde, and 6-methyl-5-heptene-2-one) in electrophysiological studies (Figure 7). The EAG responsive amplitudes to p-cymene, α-pinene, geraniol, methanol, and formaldehyde decreased significantly in the dsHparOBP14-injected insects (p < 0.05) compared to those of uninjected control insects. There are clear responses to 3-methylindole and 6-methyl-5-heptene-2-one, but their amplitudes do not change significantly after dsRNA injection (p > 0.05) (Figure 7).  3). The data was analyzed by independent-sample t-test (n = 3). The comparison between uninjected (control) or dsHparOBP14-injected insects was analyzed by independent-sample t-test (n = 3). The significant difference was presented by the stars on the top of the error bars (** for p < 0.01, *** for p < 0.001 and ns for no significant difference).

Effect of dsRNA Injection on the Behavior of H. parallela
The role of HparOBP14 in H. parallela female adults was further verified through the Y-tube behavior experiment (Figure 8). The uninjected female beetles were attracted to the methanol, p-cymene, formaldehyde, 6-methyl-5-heptene-2-one (p < 0.05), α-pinene, geraniol, and 3-methylindole (p < 0.01). The attraction of the beetles to geraniol and methanol was markedly reduced by the injection of dsHparOBP14. The response of dsHparOBP14injected female adults to α-pinene and 6-methyl-5-heptene-2-one was significantly increased (p < 0.001), but the attraction to p-cymene, 3-methylindole, and formaldehyde had no obvious change (Figure 8). Figure 8. Behavioral response of female Holotrichia parallela adults to seven compounds. Seven compounds are α-pinene, geraniol, methanol, p-cymene, formaldehyde, 3-methylindole, and 6-methyl-5-heptene-2-one. The numbers in the green bars indicate the number of females that made the choice to the Y-tube olfactometer arm treated with odor compounds. The numbers in the orange bars indicate the number of females that made the choice to the Y-tube olfactometer arm treated with the solvent paraffin oil (control). The numbers in the white bars indicate the number of females that did not make any choice and stayed in the main arm for 5 min. The significant difference in the numbers between two lateral arms was indicated by * (p < 0.05), ** (p < 0.01), *** (p < 0.001), and analyzed with a Chi-square goodness of fit test.

Discussion
The choice of the host plant by an insect female is a complex sequence that ends in oviposition, and it is a multi-modal behavior that involves, among other things, volatile and non-volatile chemical stimuli. The binding of these odor molecules to OBPs is an important step in the process of insect recognition of host plants and oviposition sites [37,42,43].
In this study, HparOBP14 was found highly expressed in the legs of H. parallela female adults among all 20 OBP genes. HparOBP14 is also highly expressed in the antennae of H. parallela. The two organs are well-known to be involved in different stages of reproductive behavior: antennae, main organs housing olfactory receptors for olfactory recognition of volatile plant compounds (VPCs), the legs carrying contact chemoreceptors allowing the recognition of non-volatile molecules once the female landed on the plant [44]. Our results show very clearly that HparOBP14 is involved in the responses of H. parallela female adults to the volatile compounds of the host plant or from the environment. The binding experiments indicate a low specificity of HparOBP14, the chemical structures of tested compounds being quite different. The olfactometer choice tests reveal a positive choice for the compounds tested upon female beetles. The EAG response of individual H. parallela female antennae to p-cymene, α-pinene, geraniol, methanol, and formaldehyde, but not to 6-methyl-5-hepten-2-one, was significantly reduced after HparOBP14 expression was knocked down. These results constitute an argument in favor of HparOBP14 s role in the recognition of the host plant volatiles.
It was reported that AlinOBP11 is the most abundant gene in the forelegs of Adelphocoris lineolatus adults, which bound to non-volatile plant secondary compounds and had an important function in taste [45,46]. However, in our EAG and Y-olfactometer experiments, the stimulus molecules were delivered to the insects by an air flow. Thus, these experiments can only assess the involvement of HparOBP14 in the perception of volatile molecules from host plants and oviposition sites. The data do not explain HparOBP14's function in the detection of non-volatile molecules by the legs.
Nevertheless, insect oviposition is mediated by attractive chemical cues emitted from oviposition sites in the environment [47]. Previous studies showed that the OBPs that are highly expressed in the legs of female insects play a role in the perception of plant volatiles, a behavioral step leading to oviposition [48]. OBP30 expressed in the legs of female Chilo suppressalis was speculated to be associated with oviposition-related compounds [49]. OBP57d and OBP57e expressed in Drosophila legs could recognize the oviposition-related volatile octanoic acid [50]. Research on the oviposition preferences found that the female beetles (Cyclocephala borealis, C. lurida, C. pasadenae, C. hirta, and C. parallela) prefer to lay eggs in soil that is rich in organic matters [51]. In this study, the leg-biased HparOBP14 had a relatively high affinity to methanol, geraniol, 3-methylindole, p-cymene, α-pinene, and formaldehyde, which are chemical cues for locating oviposition sites by H. parallela females [17]. Furthermore, knockdown of HparOBP14 expression significantly affected the behavior of H. parallela to decomposed organic fertilizer volatiles, in agreement with previous studies on the oviposition preferences of beetles [17,18,48]. The attractive response of H. parallela female adults was significantly decreased to methanol and geraniol by the HparOBP14 interference. Methanol and geraniol are volatile organic compounds emitted during the decomposition of plant litter and are also the potential volatiles to attract female H. parallela during oviposition [17,52]. Methanol is attractive not only to H. parallela, but also to the root-eating beetle Rhizophagus ferrugineus and H. oblita [18,53]. Geraniol, as a compound in bait, is attractive to both sexes of the Japanese beetle Popillia japonica [54,55] and to the leaf beetle Monolepta hieroglyphica as key plant volatiles [56,57]. These results support that HparOBP14 may play a role in oviposition behaviors of H. parallela female adults.
It is noted that there are differences between the results of the EAG experiments and those of the behavioral experiments. This may be due to different doses used in these two experiments. It is also noted that the in vitro binding results (Table 1) are different from those of the EAG and behavioral experiments. This could be due to the limitation of the fluoresce competitive binding assay used, namely 1-NPN displacement. It is not certain whether 1-NPN could bind only into the HparOBP14 binding pocket. The compounds tested could only displace some unspecific 1-NPN binding as the affinities obtained in this study are not as good as those (<10 µM) in other studies. Furthermore, the recombinant HparOBP14 also binds the host plant volatile 6-methyl-5-heptene-2-one, and the knockdown of HparOBP14 expression increased the attraction of the beetles to it but did not affect the EAG response of the female adults to this volatile compound. These responses could be mediated by several olfactory proteins; some of them could play a more important role than HparOBP14 in the response to 6-methyl-5-hepten-2-one. Therefore, HparOBP14 may be essential for female H. parallela to detect methanol and geraniol and may play a role in the process of finding suitable oviposition sites, as well as may participate in the process of identifying host plants. We cannot attribute a specific function (olfaction or taste) to HparOBP14 from the current study, but provide some evidence for the possible involvement of HparOBP14 in the perception of plant volatiles and oviposition.

Conclusions
In this study, we preliminarily identified the function of HparOBP14 which is highly expressed in the legs of female H. parallela. Based on the transcriptome data, we cloned the HparOBP14 gene and further described the tissue-specific expression and the binding spectrum as well as the changes of the electrophysiological and behavioral responses after the knockdown of HparOBP14 expression. These results provide the evidence that HparOBP14 may be involved in the oviposition behavior of H. parallela female adults, and also provide a theoretical basis for the research and development of environmentally friendly pest control strategy, which may ultimately help to reduce the reproduction of H. parallela.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/insects12050430/s1, Table S1. Primer pairs used in this study. Figure S1. Nucleotide sequence and deduced amino acid sequence of HparOBP14. Figure S2. Alignment of the deduced amino acid sequence of HparOBP14 with homologous proteins from other OBPs. Figure