Molecular and Functional Characterization of Odorant Binding Protein 7 From the Oriental Fruit Moth Grapholita molesta (Busck) (Lepidoptera: Tortricidae)

Odorant-binding proteins (OBPs) are widely and abundantly distributed in the insect sensillar lymph and are essential for insect olfactory processes. The OBPs can capture and transfer odor molecules across the sensillum lymph to odorant receptors and trigger the signal transduction pathway. In this study, a putative OBP gene, GmolOBP7, was cloned using specific-primers, based on the annotated unigene which forms the antennal transcriptome of Grapholita molesta. Real-time PCR (qRT-PCR) analysis revealed that GmolOBP7 was highly expressed in the wings of males and the antennae of both male and female adult moths, while low levels were expressed in other tissues. The recombinant GmolOBP7 (rGmolOBP7) was successfully expressed and purified via Ni-ion affinity chromatography. The results of binding assays revealed that rGmolOBP7 exhibited a high binding affinity to the minor sex pheromone 1-dodecanol containing Ki of 7.48 μM and had high binding capacities to the host-plant volatiles, such as pear ester, lauraldehyde and α-ocimene. RNA-interference experiments were performed to further assess the function of GmolOBP7. qRT-PCR showed that the levels of mRNA transcripts significantly declined in 1 and 2 day old male and female moths, treated with GmolOBP7 dsRNA, compared with non-injection controls. The EAG responses of dsRNA-injected males and females to pear ester, as well as the EAG responses of dsRNA-injected males to 1-dodecanol, were significantly reduced compared to the GFP-dsRNA-injected and non-injected controls. We therefore infer that GmolOBP7 has a dual function in the perception and recognition of the host-plant volatiles and sex pheromones.


INTRODUCTION
The sophisticated olfactory system plays an essential role in an insect's survival and reproduction. Adult insects greatly depend on olfactory cues to locate mates and optimal host plants and avoiding predators (Takken and Knols, 1999;Leal, 2013;Suh et al., 2015). In the early events of olfactory processing, airborne chemical signals must pass through the aqueous barrier of the sensillum lymph, surrounding the dendrites of the olfactory receptor neuron (ORNs) cells (Li et al., 2015).
Odorant-binding proteins (OBPs), a kind of transport protein, can selectively bind and carry hydrophobic odorants across the sensillum lymph to odorant receptors (ORs) and trigger the signal transduction pathway (Pelosi et al., 2005). The converted electrophysiological signals are then sequentially processed in the antennal lobes, mushroom bodies and central nervous area, to induce a behavioral response in specific semiochemicals of insects (Feng and Prestwich, 1997;Helfrich-Förster, 2000;Hallem et al., 2006;Pelosi et al., 2006;Leal, 2013;Yi et al., 2014). OBPs are responsible for the connection between the external environment and ORNs in the odorant-molecule recognition process. OBPs also mediates the first stage of the physiological process involved in the sensing of the external environment, by insects (Willett and Harrison, 1999;Laughlin et al., 2008;Pelosi et al., 2014;Leal and Leal, 2015).
OBPs belong to a class of small water-soluble proteins that are impregnated in the sensillum lymph at extremely high concentrations (up to 10 mM; Vogt and Riddiford, 1981;Klein, 1987;Steinbrecht et al., 1992). The first insect OBP was identified in the antennae of male Antheraea ployphemus. By using a radiolabeled photo-affinity analog, this protein was designated as a pheromone binding protein (PBP) as it specifically bound to the female sex pheromone E6, Z11-hexadecadienyl acetate (Vogt and Riddiford, 1981). Since then, OBPs have been discovered in various insect orders (Hansson and Stensmyr, 2011;Antony et al., 2018;Fleischer et al., 2018). Lepidopteran OBPs are usually subdivided into three subfamilies including PBPs, general OBPs, and antennal binding proteins (ABPX), on the basis of aminoacid sequence homologies (Hekmat-Scafe et al., 2002). The PBPs are located in the sensilla trichodea and exhibit specific binding to female sex pheromones (Bette et al., 2002;Lautenschlager et al., 2007). GOBPs (further classified as GOBP1 and GOBP2) are primarily distributed in the sensilla basiconica and their function is mainly involved in the detection of general odorants (e.g., host plant volatiles; Vogt et al., 2002;Nardi et al., 2003;Maida et al., 2005;Liu et al., 2015). In some Lepidopteran species, GOBP2 also showed high-binding affinities to sex pheromones in addition to general odorants (Liu et al., 2010(Liu et al., , 2012Li et al., 2016a). ABPX are more divergent among insects and its functions may play a similar role than PBPs or GBOPs in the discrimination and transportation of semiochemicals (Tian et al., 2018).
The binding affinities of insect OBPs to odorant molecules have been measured via fluorescence competitive binding assays with N-phenyl-1-naphthylamine (1-NPN) as a probe (Pelosi et al., 2006;Zhou, 2010). For example, Helicoverpa armigera HarmOBP17 and HarmOBP18 have strong binding capacities to β-ionone (Li et al., 2013). Locusta migratoria LmigOBP1 exhibited specific-binding affinities to pentadecanol and 2pentadecanone, where Asn74 formed the key binding site in these two ligands (Jiang et al., 2009). Grapholita molesta GmolGOBP2 had specific binding ability to the minor sex pheromone component 1-dodecanol (Li et al., 2016a). The fluorescence competitive binding assay is only in vitro and the binding functions of OBPs still need to be verified further by experiments in vivo. RNAi experiments demonstrated that OBPs are indispensable in the olfactory communication of insects. For example, the electroantennogram (EAG) responses of female Adelpocoris lineolatus to tridecanal and 1-hexanol were drastically reduced after the double-stranded RNA (dsRNA) of AlinOBP4 was injected into both female and male adult insects . The EAG values of AgosOBP2-dsRNAtreated Aphis gossypii to cotton-derived volatiles were remarkably lower than those of non-injected controls (Rebijith et al., 2016). By silencing the RferOBP1768 gene of Rhynchophorus ferrugineus, the adults apparently lose the ability to recognize the aggregation pheromone compounds 4-methyl-5-nonanol and 4methyl-5-nonanone (Antony et al., 2018).
The oriental fruit moth Grapholita molesta, is a destructive fruit pest species that causes considerable economic losses in fruit yields on a global scale (Rothschild and Vickers, 1991). The first three moth generations mainly infest peach shoots in the early growing season, whereas the third generation begins to shift and attack pear and apple orchards in the late growing season. The migration of the adults is predominantly guided by the change in volatile components emitted by these host plants (Myers et al., 2007). At present, monitoring the G. molesta mainly depends on the pheromone trapping of male moths. However, the females have multiple mating abilities and their flight capabilities are three to six times greater than that of males and the females also have higher mating rates in the pheromone trapping orchards (Hughes and Dorn, 2002;Il'ichev et al., 2007;Zhang et al., 2012). Therefore, a strategy to monitor both female and male moths, based on olfactory cues emitted from host plants, is desirable. For example, a three-compound mixture of (Z)-3-hexen-1-ol, (Z)-3-hexen-1-yl acetate, and benzaldehyde in proportion 1:4:1 can attract female G.molesta just as well as the natural blend from peach shoots can (Natale et al., 2003).
In this study, GmolOBP7 was cloned using specific-primers based on the annotated unigene from the antennal transcriptome of G. molesta. qRT-PCR was performed to determine the expression patterns of GmolOBP7 in different tissues, genders, and developmental stages of the G. molesta. The binding affinities of the rGmolOBP7 with sex pheromone components and the host plants' volatiles, were measured via fluorescence binding assays. Furthermore, the ligand-binding functions of GmolOBP7 were further verified in vivo by knocking down the GmolOBP7 gene. The olfactory mechanism of the oriental fruit moth was further explicated to provide a theoretical basis for the design and implementation of control strategies against this fruit pest.

Insect Samples
G. molesta individuals were obtained from the College of Plant Protection, Northwest A&F University, Yangling, Shaanxi, China. The laboratory colony has been maintained for more than 90 generations. The larvae were reared on an artificial diet at 25 ± 1 • C, 70% ± 5% RH under a day/night cycle of 15:9, until pupation (Du et al., 2010). After pupation, male, and female pupae were placed in separate glass tubes and maintained under the same conditions described above. The adults were fed 5% honey solution daily. To detect the tissue distribution of GmolOBP7 in adult moths, various tissues (including antennae, heads without antennae, thoraces, abdomens, legs, and wings) were collected from 3-day-old males and females and immediately transferred to 1.5 mL Eppendorf tubes immersed in liquid nitrogen. All samples were stored at −80 • C prior to use. In order to determine the transcript level of GmolOBP7 in different developmental stages of the G. molesta, samples of eggs, larvae (including 1st, 2nd, 3rd, 4th, and 5th instars), pupae (including prepupae and later-pupae) and adults (including 1-d-old, 3-d-old, and 5-d-old adults) were collected and stored at −80 • C prior to use.

RNA Extraction, OBP Cloning, and Sequencing
Total RNA of all samples was extracted using a RNAiso Plus reagent (TaKaRa, Daian, China) according to the manufactures' instructions. The residual genomic DNA in the total RNA was removed using DNase I (Thermo Scientific, USA), and the first-strand cDNA was synthesized in accordance with the recommended protocols of the RevertAid First Strand cDNA Synthesis Kit (Thermo Scientific, USA). The products were stored at −80 • C.
The predicted coding region of GmolOBP7 was cloned using specific-primers based on the annotated unigene from the antennal transcriptome of G. molesta ( Table 1). The predicted results showed that GmolOBP7 had no signal peptide at the N-terminus of the amino acid sequence. In order to confirm whether we acquired the complete coding sequence of GmolOBP7, gene-specific primers were synthesized and used for 5' RACE (rapid amplification of cDNA ends; Table 1) Refers to a procedure in a previous study (Luo et al., 2011). A first 41-cycle touchdown PCR was performed using 5' RACE outer primers (named outer5F and outer5R; Table 1). A 25 µL PCR reaction system contained 12.5 µL of 2 × Super Pfx MasterMix (CWBIO, Beijing, China), 0.8 µL of each primer (10 µM), 1 µL of sample cDNA, and 9.9 µL of nuclease free water. The thermocycling program included denaturation at 95 • C for 5 min, followed by 16 cycles of 30 s at 95 • C, 1 min at 65 • C, and 2 min at 72 • C, and the annealing temperature was decreased 1 • C each four cycles. The remaining 25 cycles consisted of 30 s at 95 • C, 1 min at 61 • C, and 2 min at 72 • C, and a final extension step of 72 • C for 10 min. The PCR products were diluted 80 times with sterilized ddH 2 O. The second 41-cycle touchdown PCR was conducted using 5' RACE inner primers (named inner5F and inter5R; Table 1) and the template with diluted PCR products. The reaction system and procedure was the same as the first round of PCR. The amplified product was purified with an Universal DAN Purification Kit (TianGen, Beijing, China), and cloned into the pMD R 19-T cloning vector (TaKaRa, Dalian, China) and then transformed into DH5α Escherichia coli competent cells (TianGen, Beijing, China). Five positive clones were randomly selected for sequencing at the Aoke Biotech Company (Aoke, Xi'an, China).

Sequence and Phylogenetic Analyses
The online programs of ORF Finder (http://www.ncbi.nlm. nih.gov/gorf/gorf.html), SignalP 4.0 (http://www.cbs.dtu. dk/services/SignalP/), and ExPASy server (https://web. expasy.org/compute_pi/) were used to predict the ORFs, signal peptides and the molecular weight and isoelectric point of mature protein of GmolOBP7, respectively. The amino-acid sequences were aligned using ClustalX 1.83 software. A phylogenetic tree was established by the MEGA6.0 software using the neighbor-joining method (NJ) with 1,000 bootstrap replications, and the tree was drawn using Adobe Photoshop CS5.

Expression Analysis Using qRT-PCR
The expression levels of GmolOBP7, in different tissues and developmental stages of G. molesta, were measured via qRT-PCR. All qRT-PCR experiments were performed according to the MIQE Guidelines (Bustin et al., 2009). Specific primers were designed using the program Primer3-blast (https://www. ncbi.nlm.nih.gov/tools/primer-blast/), available online ( Table 1). The elongation factor 1-alpha gene (EF1-α) (GenBank No: KT363835.1) and the β-actin gene (GenBank No: KF022227.1) were used as reference genes. The reactions were performed on a CFX96 Real-Time PCR Detection System (Bio-Rad, USA). Each amplification reaction was conducted using a 20 µL reaction system containing 10 µL of 2×SYBR R Premix Ex Taq TM II mixture (TaKaRa, Dalian, China), 0.8 µL of each primer (10 µM), 1 µL of sample cDNA, and 7.4 µL of nuclease-free H 2 O. Samples without a template cDNA served as negative controls.
To check reproducibility, test samples and negative controls were performed in triplicates. qRT-PCR was performed via initial denaturation at 95 • C for 30 s, followed by 40 cycles of 95 • C for 5 s, 60 • C for 30 s and 72 • C for 30 s. The melting curves were used to examine primer specificity, and the standard curves were used to determinate the amplification efficiencies of target and reference genes. The expression levels of GmolOBP7 in different adult tissues and development stages were performed based on previous methods (Livak and Schmittgen, 2001;Liu et al., 2016). The expression level of all samples of the GmolOBP7 was calculated using the value of the amplification efficiency (E) and the value of the cycle threshold (Ct) (Equation 1). The normalized expression level of the tested samples was calculated by the geometric means of the expression level of the reference genes (β-actin and EF-1α) (Equation 2) (Vandesompele et al., 2002). The significant differences in different tissues and developmental stages were analyzed by the Tukey's HSD tests with a critical level of α = 0.05. The paired t-test was used to measure the impacts of the expression of GmolOBP7 between male and female moths. All the data were analyzed using SPSS 18.0 software (SPSS Inc., Chicago, IL, USA).
Normalized expression level of target gene (2)

Expression Vector Construction
Specific primers with restriction enzyme sites were designed to clone the coding region of the GmolOBP7 (Table 1), and the OBP7i-reverse TAATACGACTCACTATAGGGACCAAACTGAGCAGCGTTTT GFP-forward TAATACGACTCACTATAGGGGTGTTCAATGCTTTTCCCGT 315

GFP-reverse TAATACGACTCACTATAGGGCAATGTTGTGGCGAATTTTG
The restriction endonucleases are in parentheses after each primer, and the restriction sites are underlined. "-" means the amplified DNA fragment is an unknown in size.
PCR product was then cloned into the pMD R 19-T cloning vector (TaKaRa, Dalian, China) and sequenced. The recombinant plasmid pMD R 19-T/GmolOBP7 and the expression vector pET32a(+) (Novagen, Madison, WI, USA) were digested with the same restriction endonucleases, and the released DNA fragment was cloned into pET-32a(+) and then transformed into BL21 E. coli competent cells (Tiangen, Beijing, China). A positive clone containing pET32a(+)/GmolOBP7 was further confirmed by sequencing.

Protein Expression and Purification
The overnight bacterial solution was diluted with 750 mL of LB medium (with 100 mg/mL ampicillin) and cultured at 37 • C until its cell density reached a value of OD 60 = 0.6. The cultures were induced by adding isopropyl-β-Dthiogalactoside (IPTG) at a final concentration of 0.5 mM for an additional 5 h at 37 • C, 220 rpm. The bacterial cells were harvested by centrifugation (10 min at 8,000 rpm, 4 • C), and the pellets were then sonicated in a lysis buffer (1 mM phenylmethanesulfonyl fluoride, 250 mM NaCl, and 20 mM Tris-HCl pH 7.4) and centrifuged again (13,000 g, 30 min, 4 • C). A sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) analysis revealed that the rGmolOBP7 was mainly present in supernatants. The supernatant of rGmolOBP7 was enriched by a Ni-NTA His·Bind Resin column (7 sea Pharmatech Co., Shanghai, China) in accordance with the manufacturer's instructions. To avoid the effects of His-tag on subsequent experiments, the His-tag was cleaved by a recombinant enterokinase (NEB, Beijing, China) and removed using the column mentioned above. The target protein was purified using affinity chromatography, and the concentration was determined by the BCA protein kit (Beyotine, Shanghai, China).

Fluorescence Binding Assays
Fluorescence intensity was detected on a spectrophotofluorometer (F-4500, Hitachi, Japan) at room temperature using a quartz cuvette with a 1 cm light path. The silt width of excitation and emissions were all 10 nm. The fluorescence probe 1-NPN was excited at 337 nm and the emission spectra were recorded between 370 and 550 nm. Four sex pheromone components and 31 volatiles, derived from peach shoots and pear fruits, were selected for binding assays ( Table 2). The probe 1-NPN and tested ligands were all dissolved in spectrophotometric-grade methanol to obtain a 1 mM stock solution. The binding affinity of rGmolOBP7 with 1-NPN was measured by adding aliquots of 1-NPN to a 2 µM protein solution (diluted with 20 mM Tris-HCl pH 7.4) to final concentrations of 0 to 18 µM.
To test the binding affinities of GmolOBP7 to various ligands, 2 µM solutions of rGmolOBP7 and 1-NPN were titrated with the 1 mM solution of each ligand to a final concentration of 0-14 µM for sex pheromones and 0-35 µM for host volatiles. The corresponding florescence intensity values were collected as three Frontiers in Physiology | www.frontiersin.org

dsGmolOBP7 and dsGFP Synthesis
The specific-primers, including T7 RNA polymerase promoter, were designed to clone DNA fragments of GmolOBP7 for 317 bp and green fluorescent protein (GFP) for 315 bp ( Table 1).
The purified PCR products were used as a template for dsRNA (dsGmolOBP7 and dsGFP) synthesis using the T7 RiboMAXTM Express RNAi System kit (Promega, USA) according to the manufacturer's instructions. The purified dsRNA was quantified via spectrophotometry (SimpliNano, GE, USA), and the dsRNA integrity was monitored by electrophoresis on 1.5% agarose gel. Each dsRNA sample was dissolved in nuclease-free water to a final concentration of 3,500 ng/µL.

dsRNA Microinjection
Based on the expression patterns observed at different insect stages, 5-day-old G. molesta pupae (later-pupae) were selected to receive a dsRNA microinjection. The conjunctivum between the prothorax and mesothorax, the conjunctivum between the mesothorax and metathorax and the conjunctivum between the thorax and abdomen were initially selected as putative injection sites. Then, 39, 69, and 138 nL of RNase-free H 2 O were injected into different conjunctiva, respectively. A total of 69 nL (approximately 241.5 ng dsRNA) of dsGmolOBP7 or dsGFP was injected into the appropriate injection sites of each 5-dayold pupa, by using a PL1-100 Pico-Injector (Harvard Apparatus, Holliston, MA, USA) operated by an MP-255 Micromanipulator (Sutter, Novato, CA, USA). Each type of dsRNA was injected into 600 male and 600 female moths. The heads (with antenna) of 1-, 2-, 3-, and 4 day-old adults were dissected and immediately stored at −80 • C prior to use. The total RNA and first-strand cDNA were obtained in accordance with previous methods. The specimens were used for qRT-PCR analysis to test the reduction in GmolOBP7 transcription. Experiments were performed with three biological replicates and three technical replicates.

EAG Assays
EAG responses of dsRNA-injected (including dsGmolOBP7 and dsGFP) moths and non-injected controls, to sex pheromones and host plant volatiles were detected using Electroantennography. All stimulants were diluted with liquid paraffin to the final concentration of 10 mg/mL. Liquid paraffin and cis-3-hexenyl acetate were used as the blank and reference control, respectively. Both ends of adult antennae were cut and blocked with a drop of Spectra R 360 Electrode Gel (Parker Laboratories, Fairfield, USA). The basal section was connected to the reference electrode while the distal end was linked to the recording electrode. A filtered humidified air stream was delivered by a Syntech stimulus controller (CS55 model, Syntech, Germany) at a constant flow rate of 50 cm/s, and the time of stimuli flow was 0.5 s. Filter paper strips (0.6 cm × 4.5 cm) were dripped with 15 µL of chemical solution as a stimulus source and inserted into a 1.5 mL micropipet tip. A set of stimulants consisted of four sex pheromones and five host-plant volatiles which can be strongly bound with rGmolOBP7. Each antenna measured the group of randomly arranged stimulants described above. The antennae were stimulated one time with liquid paraffin and cis-3-hexenyl acetate and dissolved in solvent, before and after each group stimulation, in order to ensure the tested antenna were activated and the connecting pipe was not contaminated by stimulants. Recordings per stimulant were taken and the antennal responses were recorded. Eight male and female antennae were measured for each stimulant. The paired t-test was used to determine whether the difference in EAG values were significant between the dsRNA injected moths (dsGmolOBP7 and dsGFP) and the non-injected control. All the data were analyzed using SPSS 18.0 software (SPSS Inc., Chicago, IL, USA).

Identification of GmolOBP7 in G. molesta
We obtained the ORF of GmolOBP7 (GeneBank No. MF066359) using ordinary PCR and 5' RACE PCR based on the annotated unigene from the antennal transcriptome of G. molesta. The ORF of GmolOBP7 is 504 bp in length and encodes 167 amino acids ( Figure S1). GmolOBP7 possesses a common characteristic of known classical-OBPs with six-conserved cysteine motifs (Figure 1). The mature protein has a predicted molecular weight of 21.47 kDa and a theoretical pI of 6.85. The SignalP 4.1 server prediction indicated that GmolOBP7 did not have signal peptides at the N-terminus of the amino acid sequence. GmolOBP7 shares the highest identities with Ectropis obliqua EoblOBP1 and Spodoptera exigua SexiOBP11, with an identity of 67 and 64%, respectively. Phylogenetic analysis showed that GmolOBP7 were clustered into a small branch close to ABPXs from Ectropis obliqua and Cnaphalocrocis medinalis (Figure 2).

Expression Profiles of GmolOBP7 in G. molesta
The amplification efficiency and melting curve of target and reference genes showed that the two specific-primers were similar in amplification speed and no nonspecific products were produced ( Figure S2), therefore, the primers can be used for relative quantification. GmolOBP7 was detected in all tested tissues of both female and male moths ( Figure 3A), but the expression quantity was higher in the wings of males and antennae of both sexes, than that in other tested tissues. The expressed quantity of GmolOBP7 was different between males and females, with the quantity in male wings and legs being significantly higher than in female wings and legs (with about 13.31 and 2.99 times differences, respectively), and was higher in female thoraces than in male thoraces (about 2.43-fold higher). The expression levels of GmolOBP7 in different developmental stages were performed by qRT-PCR (Figure 3B), the highest expression quantity was found in adults, followed by the eggs, 1st instar larvae, and later pupae (5-day-old pupae), the expression levels of GmolOBP7 were extremely low in second-to fourthinstar larvae.

Expression and Purification of GmolOBP7
Recombinant GmolOBP7 (rGmolOBP7) was successfully expressed in E. coli as a soluble protein ( Figure 4A). After being purified, about 38 kDa recombinant protein with His-tag was obtained. To avoid the effect by the His-tag on subsequent binding assays, the His-tag of pET32(a+)/GmolOBP7 was cleaved by enterokinase and then removed by Ni-NTA His·Bind Resin and SDS-PAGE analysis showed that rGmolOBP7 had a higher purity after being purified for a second time via affinity chromatography ( Figure 4B). The purified recombinant protein was then tested for the binding affinity with various ligands.

Effect of RNAi Treatment on the Expression Level of GmolOBP7
Pilot experiments showed low eclosion rates when the pupae were injected with 138 nL of water at all putative conjunctiva (<30%). The emergence rates of the pupae injected with 39 and 69 nL of water at different candidate conjunctiva ranged from 78.6 to 84.2%. Thus, the conjunctivum between the prothorax and mesothorax was selected as the appropriate injection site, and 69 nL was selected as the appropriate dosage.
qRT-PCR analysis revealed that the transcription levels of dsGFP-injected moths had no significant differences compared to the non-injected male and female moths. The transcript levels of GmolOBP7 decreased to 52.57% (with 1-d eclosion) and 67.68% (with 2-d eclosion) in GmolOBP7-dsRNA-injected males compared to that in the dsGFP-treated and non-treated controls (Figure 7A), and decreased to 59.50% (with 1-d eclosion) and 77.17% (with 2-d eclosion) in GmolOBP7-dsRNA-injected females compared to that in the controls (Figure 7B). The transcript levels of GmolOBP7 in dsRNA-treated moths were increased to normal values after 3-days of eclosion. Thus, 1-dayold adult moths were selected for subsequent EAG assays.

Electrophysiological Experiments
The EAG response values of dsGFP-treated moths to nine tested stimulants had no significant differences to the non-injected male and female moths (Figure 8). The t-tests showed that the responses of both female and male moths, to pear ester were significantly reduced (P < 0.05) after injection with GmolOBP7-dsRNA, and the response value of male moths to 12:OH was also significantly decreased. However, the response to (Z)-8-dodecenyl acetate, (E)-8-dodecenyl acetate, (Z)-8-dodecenyl alcohol, (Z)-3-hexenyl acetate, lauraldehyde, α-pinene, and αocimene was not significantly different between dsRNA-treated moths and the non-injected control.

DISCUSSION
The OBP family genes are composed of many highly differentiated subfamily genes, the olfactory function of those OBPs highly-expressed in the antennae, such as GOBPs and PBPs, has been studied extensively (Zhou et al., 2009;Yin et al., 2012;Khuhro et al., 2017). However, the ligand-binding capacities of OBPs which had no antenna-specific expression or were lowly-expressed in antennae, remains poorly understood. We identified 26 OBPs from the antennal transcriptome of G. molesta (Li et al., 2015), the antenna-highly-expressed OBPs (PBP1-3, GOBP1-2, OBP8, OBP11, and OBP15) all with their preferred odorant ligands, such as GmolPBP2, which can bind specifically to the major sex pheromone components Z8-12:Ac and E8-12:Ac (Song et al., 2014), GmolGOBP1 have strong binding affinities to the major sex pheromone component Z8-12:OH and plant volatile decane (Li et al., 2016a), while hexanal was the preferred ligand of GmolOBP15 (Li et al., 2016b). We speculated that the relatively low-expression antennal OBPs may play a role in capturing and transporting the specific compounds of the host plant volatiles. The RPKM (reads per kilobase per million mapped reads) value of GmolOBP7 ranked 18th in 28 OBPs (Li et al., 2015), and can be expressed in soluble forms in a prokaryotic system. Therefore, we selected GmolOBP7 to evaluate its role in perceiving and recognizing the trace components emitted from peach shoots and pear fruits.
The expression profiles of olfactory-related genes in different tissues and sexes can provide clues to understand their physiological function (Ju et al., 2014). Numerous experiments have revealed that the antennae-enriched OBPs play an important role in detecting sex pheromones and host plant compounds (Sun et al., 2014;Yang et al., 2016;Khuhro et al., 2017). GmolOBP7 was expressed at relatively high levels in the antennae compared to other tissues, and might have potential functions in olfactory chemoreception. G. molesta reached the peak of mating after emergence (2-to 3-days) after which the flourishing period of oviposition of 3-to 5-day-old female adults occurred. The transcript levels of GmolOBP7 were slightly higher in 3-day-old male adults than in females of the same age, and the expressed levels were enhanced slightly in 3-dayold female adults. These expression characteristics implied the GmolOBP7 may be involved in the detection of sex pheromones and host-plant volatiles. In addition to antennae, GmolOBP7 was also abundantly expressed in the male wings of G. molesta, while similar expression profiles were found in BodoOBP17 from Bradysia odoriphaga (Zhao et al., 2018), MsepOBP19 from FIGURE 3 | Expression profiles of GmolOBP7 in different tissues (A) and developmental stages (B) of male and female moths. An, antennae; He, heads; Th, thoraces; Ab, abdomens; Le, Legs; Wi, Wings; 1st, first-instar larvae; 2nd, second-instar larvae; 3rd, third-instar larvae; 4th, fourth-instar larvae; 5th, fifth-instar larvae; Pup, prepupae; Later Pup, 5-d-old pupae; 1♂, 1-d-old adult males; 1♀, 1-d-old adult females; 3♂,3-d-old adult males; 3♀, 3-d-old adult females; 5♂, 5-d-old adult males; 5♀, 5-d-old adult females. Different lowercase and capital letters indicate significantly different expression levels among different tissues of female and male, respectively (Tukey's test, α = 0.05). Asterisks indicate significant different expression levels of GmolOBP7 between two sexes in the same tissue (Independent t-test, α = 0.05). Mythimna separata (Chang et al., 2017), AmalOBP8 from Agrilus mali (Cui et al., 2018), and AlucOBP6 from Apolygus lucorum (Hua et al., 2012). The chemoreception sensilla have been found on the wings of A. mali, as well as the taste organ and taste bristles and were also located on the wings of Drosophila melanogaster (Galindo and Smith, 2014). We speculated that the GmolOBP7 may play an important role in olfactory or gustatory perception, and further studies with non-volatile secondary metabolites of host plants are needed to verify this.
G. molesta thrives mainly on plants of the rosaceae family, and the peach and pear are considered the optimal host plants (Rice et al., 1972;Rajapakse et al., 2006). Plant volatiles serve as olfactory cues for G. molesta orientation, and guide the adults to switch from peach orchards to pear orchards during the growing season (Zhao et al., 1989;Najar-Rodriguez et al., 2013). We selected four sex pheromone components and 31 potential host-plant volatiles or its analogs, to determine the binding characteristics of rGmolOBP7. The sex pheromone components have been identified and widely used in the sexual trapping of male G. molesta (Cardé et al., 1975(Cardé et al., , 1979Reinke et al., 2014). The tested volatiles are known to be emitted from peach shoots and pear fruits. EAG studies on G. molesta have shown that some saturated and unsaturated volatile components of aldehydes, alcohols, acetate esters, terpenes and benzonitriles can effectively elicit responses from the antennal lobes of adult moths (Natale et al., 2003;Piñero and Dorn, 2009;Lu et al., 2015). Behavior response assays also indicated that the individual volatile component or mixture of several volatile compounds caused obvious attraction in adult moths (Natale et al., 2004;Piñero et al., 2008;Il'ichev et al., 2009;Yu et al., 2015).
The binding assays showed that GmolOBP7 has broad binding activities to various ligands including aldehydes, alcohols, esters, terpenoids and nitriles compounds. Pear ester, lauraldehyde, and dodecanol were the first three strongest ligands that bound to rGmolOBP7. Previous reports confirmed that the minor sex pheromone component 12:OH only elicited a weak EAG response to male antennae of G.moletsa. Its main function is a synergist attractant, that increases the frequency of male landing and is a stimulus that induces mating behavior when the male and female are close to each other, or when the male is close to pheromone lures (Cardé et al., 1975(Cardé et al., , 1979. EAG responses of GmolOBP7-dsRNA-treated males to 12:OH were significantly reduced compared with GFP-dsRNA-injected and non-injected controls. The simplest explanation is that GmolOBP7 may be involved in the perception of the sex pheromone 12:OH, and a behavioral response test of GmolOBP7-daRNA-treated to 12:OH is required to confirm this in future studies. Pear ester (Ethyl (E,Z)-2,4-decadienoate) belongs to a volatile derived from pear fruits and is widely applied in trapping female codling moth, Cydia pomonella, which is a closely related species of G. molesta (Vanessa et al., 2008). Pear ester exhibited the strongest binding affinity with GmolOBP7, and the EAG response values of dsRNA-treated males and females to pear ester were significantly decreased. GmolOBP7 may play the same role in perception of pear ester in male and female moths. rGmolOBP7 showed strong binding ability to lauraldehyde, but the EAG responses of dsRNA-treated male and female moths, to this compound, were not significantly different compared to non-injected controls. OBPs have a binding pocket formed by a six-α-helix fold, and usually have similar binding affinities to the ligand with the same structure and size. For example, Locusta migratoria LmigOBP1 binds to pentadecanol (C15), 2-pentadecanone (C15) and ethyl tridecanoate (C15) (Jiang et al., 2009), Bombyx mori BmorGOBP2 binds to (10E,12Z)hexadecadien-1-ol (bombykol) and (10E,12Z)-hexadecadienal (bombykal) (Zhou et al., 2009), Loxostege sticticalis LstiGOBP2 binds to 1-hexanol and 1-hexanal (Yin et al., 2012). We speculated that GmolOBP7 bound to lauraldehyde because of its size. Similar to pear ester and 12:OH, the lauraldehyde is also a derivative of a linear aliphatic hydrocarbon with 12 carbon atoms in the main chain. The molecular size of these three compounds are similar.
The binding assays were performed as recombinant OBPs expressed in vitro and the binding of OBPs with the ligands are affected by the shape and amino acid residues of the binding pocket of proteins, as well as the carbon-chain lengths, functional groups, isomers, and C = C bonds of ligands (Sandler et al., 2000;Mohanty et al., 2004;Wogulis et al., 2006;Li et al., 2008;Christina et al., 2017). OBPs may bind to many tested ligands  Table 2. FIGURE 7 | Expression level of GmolOBP7 in dsRNA-treated and no-injected moths at different ages of adult moths. (A) male moths; (B) female moths. About 240 ng (69 nL) of GmolOBP7 dsRNA and dsGFP were injected into 5-day-old pupae, respectively. Asterisk Different letters indicate significantly different expression levels of GmolOBP7 between dsRNA-treated moths and non-injected moths (independent t-test, α =0.05).
with similar structures or sizes. Whether the odorants with strong binding activity to OBPs play a role in chemoreception such as mating and host selecting in insects, still needs to be verified by electrophysiological and behavioral assays. The methods of the RNAi combined with an EAG assay is an effective way to verify whether the binding-active odorants can be recognized by insects . We found that GmolOBP7 exhibited binding activities in 21 of 35 tested ligands. The EAG assays preliminary revealed that GmolOBP7 may be involved in the detection of 12:OH and pear ester, however, FIGURE 8 | Electrophysiological response of Grapholita molesta to nine stimulants after RNAi knockdown. (A) male and (B) female moths. dsGmolOBP7 and dsGFP indicates treated moths injected with GmolOBP7 dsRNA and GFP dsRNA, respectively. Each treatment included 8 months, the antennae were stimulated with 15 µL volatile compound dissolved in liquid paraffin (ck), and ck and (Z)-3-hexenyl acetate were used to stimulate the antennae before and after a group of volatiles stimulation. Different letters indicate significant differences between dsRNA-treated moths and non-injected moths (independent t-test, α = 0.05).
whether GmolOBP7 participates in the perception of other remaining binding-active odorants, requires further functional verification.

AUTHOR CONTRIBUTIONS
X-LC, G-WL, X-LX, and J-XW conceived and designed the experimental plan. X-LC and G-WL performed the experiments. X-LC, J-XW, and X-LX analyzed and processed the data. X-LC wrote the paper. All authors read and agreed to publish this paper.