Guang-Wei Li1,2†
Xiu-Lin Chen
Jun-Xiang Wu- 1Shaanxi Province Key Laboratory of Jujube, Yan’an University, Yan’an, China
- 2College of Life Sciences, Yan’an University, Yan’an, China
- 3Key Laboratory of Plant Protection Resources and Pest Management (Northwest A&F University), Ministry of Education, Yangling, China
Chemosensory proteins (CSPs) belong to a family of small water-soluble proteins that can selectively bind and transport odorant molecules for olfactory communication in insects. To date, their definite physiological functions in olfaction remain controversial when compared with odorant binding proteins (OBPs). To investigate the functions of CSPs in the oriental fruit moth Grapholita molesta, we determined the tissue expression patterns and binding properties of the CSP, GmolCSP8. The key binding sites of GmolCSP8 with a representative ligand were evaluated using molecular flexible docking, site-directed mutagenesis and ligand-binding experiments. Multiple sequence alignment and phylogenetic analysis showed that GmolCSP8 possesses a typical conserved four cysteines motif and shares high sequence identity with some CSP members of other Lepidopteran insects. GmolCSP8 was predominantly expressed in the wings and antennae of both male and female adults and may be involve in contact chemoreception. Recombinant GmolCSP8 (rGmolCSP8) exhibited specific-binding affinities to small aliphatic alcohols (C4–12) and had the strongest binding affinity to 1-hexanol. The three-dimensional structure of GmolCSP8 was constructed using the structure of sgCSP4 as a template. Site-directed mutagenesis and ligand-binding experiments confirmed that Thr27 is the key binding site in GmolCSP8 for 1-hexanol binding, because this residue can form hydrogen bond with the oxygen atom of the hydroxyl group in 1-hexanol, and Leu30 may play an important role in binding to 1-hexanol. We found that pH significantly affected the binding affinities of rGmolCSP8 to ligand, revealing that ligand-binding and -release by this protein is related to a pH-dependent conformational transition. Based on these results, we infer that GmolCSP8 may participate in the recognition and transportation of 1-hexanol and other small aliphatic alcohols.
Introduction
The oriental fruit moth, Grapholita molesta Busck (Lepidoptera: Tortricidae), is a major pest species of fruit trees belonging to the family Rosacea, and has caused substantial losses in fruit yields throughout the world (Myers et al., 2007). The adult moths exhibit pronounced seasonal host-transfer behavior in adjacent areas of peach and pear, or peach and apple orchards (Rajapakse et al., 2006). The larvae mainly infest peach shoots in the early growing season and then shift to damage pear or apple fruits in the vicinity of peach orchards in the mid-to-late season (Piñero and Dorn, 2009; Lu et al., 2012; Najar-Rodriguez et al., 2013a). Selection of host location and oviposition behaviors of G. molesta are mediated by semiochemicals [e.g., (Z)-3-hexen-1-ol, 1-hexanol, (E)-2-hexenal, hexanal, heptanal, octanal, nonanal, benzaldehyde, (Z)-3-hexen-1-yl acetate, butyl hexanoate, (Z)-β-ocimene, and (E)-β-farnesene] emitted from host-plants (Natale et al., 2003; Il’ichev et al., 2009; Piñero et al., 2008). Changes in the components of host-plant volatiles are the primary cause of host transfer for this pest (Najar-Rodriguez et al., 2013b). G. molesta is predominantly monitored and controlled by sex pheromone trapping of male moths. However, this species is able to undergo multiple mating (up to fours), and the females can present a higher mating rate in orchards controlled by sex attractants (de Morais et al., 2012). In addition, the females, especially gravid females exhibit stronger flight performance than males, and have the capacity for inter-orchard flights (Steiner and Yetter, 1933). Female G. molesta are considered the main colonists and can cause a serious threat to pear and apple orchards in the vicinity of peach crops. Development of a control strategy for female G. molesta based on host-plant-volatiles is very promising. Thus, investigating the function of olfaction-related genes will help to better understand the olfactory system of G. molesta.
The insect olfactory system is highly specific and sensitive, and can detect and discriminate semiochemicals from complex environments and transfer signals to the nervous system to initiate specific behavioral responses, such as identifying conspecifics, finding mates, mating with partners, locating host plants, selecting oviposition sites, and avoiding toxins and predators (Sato and Touhara, 2009; Antony et al., 2016). Olfaction perception in insects is regulated through numerous olfactory proteins present in the sensilla lymph of the antennae, including chemosensory proteins (CSPs), odorant binding proteins (OBPs), odorant receptors (ORs), ionotropic receptors (IRs), sensory neuron membrane proteins (SNMPs), and odorant degrading enzymes (ODEs) (Bruyne and Baker, 2008). CSPs and OBPs are small-soluble carrier proteins, which function during the initial stages of chemoreception in insects. Members of these two groups of proteins can selectively bind and transport hydrophobic odorant molecules across the aqueous sensillum lymph to ORs located on the dendrite membranes of olfactory receptor neuron (ORN) cells, activating olfactory-related signal transduction (Laughlin et al., 2008).
Compared with the expression patterns of OBPs, which are specific or highly expressed in antennae, CSPs are expressed in chemosensory tissues, including antennae (Jacquin-Joly et al., 2001; Liu et al., 2010; Zhang et al., 2013, 2014), labial palps (Jin et al., 2005), maxillary palps (Maleszka and Stange, 1997), proboscis (Nagnan-Le et al., 2000), legs (Wanner et al., 2004), and wings (Zhou et al., 2008), and in non-chemosensory tissues, such as the ejaculatory bulb (Dyanov and Dzitoeva, 1995; Sirot et al., 2008), and pheromone gland (Jacquin-Joly et al., 2001) in numerous insect species, suggesting multiple roles, beyond chemosensation. Abundant CSPs in antennae are involved in the perception and recognition of sex pheromones or host-plant volatiles. For example, Monochamus alternatus MaltCSP5 exhibited strong binding affinities to pine volatile components myrcene, (+)-β-pinene, and (-)isolongifolene (Ali et al., 2019); Mythimna separata MsepCSP8 can bind to plant volatiles hexanal and α-terpinene (Younas et al., 2018); and Microplitis mediator MmedCSP3 showed high binding affinities to sex pheromone components cis-11-hexadecenyl aldehyde (Z11-16: Ald), cis-11-hexadecanol (Z11-16: OH), and trans-11-tetradecenyl acetate (E11-14: Ac) (Peng et al., 2017). RNAi and behavioral tests showed that Nilaparvata lugens NlugCSP8 is involved in the perception of plant volatiles nerolidol, hexanal, and trans-2-hexenal (Waris et al., 2018). The non-olfactory functions of insect CSPs have also been documented, and include embryonic integument formation in Apis mellifera (Maleszka et al., 2007), leg regeneration in Periplaneta americana (Kitabayashi et al., 1998), tissue-remodeling in Drosophila melanogaster (Sabatier et al., 2003), and regulation of circadian cycles in Frankliniella occidentalis (Claridge-Chang et al., 2001). The proposed olfactory function of insect CSPs is supported by genes expression in the lymph of chemosensilla and ligand-binding properties of recombinant proteins to odor molecules observed in vitro, however, owing to a lack of direct evidence on the molecular recognition of semiochemicals by CSPs, their exact physiological functions and molecular mechanisms still remains unclear.
In our previous study, we identified 17 CSP genes from the antennae of G. molesta using RNA-Sequencing (RNA-seq); five of these were selected with reads per kilobase per million mapped reads (RPM) value that was high in the antennae transcriptome, and successfully expressed in a prokaryotic expression system. Ligand-binding assays demonstrated that four recombinant GmolCSPs (rGmolCSPs) exhibited only weak binding affinities to host plant volatiles (Li G.W. et al., 2015). In this study, we found that a CSP, GmolCSP8, demonstrated binding specificity to small aliphatic alcohols. The main objective of this study was to identify the functions of GmolCSP8 during chemoreception. We performed real-time quantitative PCR (qRT-PCR) to monitor the expression of GmolCSP8 in different tissues of the adult males and females. Binding properties of GmolCSP8 were assessed using four sex pheromones and host-plant volatiles as putative ligands. Based on the results of ligand-binding assays, we combined homology modeling, molecular docking, site-directed mutagenesis, and ligand-binding assays to investigate the binding sites of GmolCSP8 to 1-hexanol and proposed a possible ligand-binding mechanism.
Materials and Methods
Insect Rearing and Tissue Collection
Grapholita molesta in all experiments were obtained from our laboratory colony, which has been maintained for more than 90 generations at the College of Plant Protection, Northwest A&F University, Yangling, China. The larvae were reared on an artificial diet and kept in an artificial climate chest under conditions of 25 ± 1°C, 70 ± 5% relative humidity (RH), and a 15:9 (light: dark) photoperiod (Du et al., 2010). Pupae were sexed, and male and female moths were kept in separate glass jar after emergence. The adults were fed a 5% honey solution and maintained under the conditions described above. To determine the tissue expression profiles of the target CSP gene in G. molesta, the antennae (200 moths per sample), heads (without antennae; 40 moths per sample), thoraxes (without wings and legs; 10 moths each sample), abdomens (five moths per sample), legs (100 moths per sample), and wings (100 moths per sample) from both sexes were excised (three replicates) from 3 to 4-day-old moths and immediately transferred to 1.5 mL centrifuge tubes immersed in liquid nitrogen and prepared for total RNA isolation.
RNA Extraction and cDNA Synthesis
Total RNA from each sample was isolated using RNAiso Plus extraction reagent (TaKaRa, Shiga, Japan). Total RNA integrity was evaluated using 1.5% agarose electrophoresis, and the concentration was quantified by a spectrophotometry (SimpliNano, GE, United States). Residual genomic DNA was removed using DNase I (Thermo Scientific, United States), and 1 μg of RNA was reverse-transcribed into first-strand cDNA with an oligo (dT)18 primer using Revertaid First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, Waltham, MA, United States). The products were stored at -80°C.
ORF Amplification and Phylogenetic Analysis
Based on the annotated Unigene from the antennal transcriptome of G. molesta (Li G.W. et al., 2015), full-length cDNA encoding GmolCSP8 was amplified with female antennae cDNA using ExTaq DNA polymerase (Thermo Scientific, United States) and a pair of specific-primers (Table 1). PCR conditions were as follows: initial denaturation at 95°C for 5 min, followed by 35 cycles of 95°C for 30 s, 58°C for 30 s, and 72°C for 30 s, and final extension at 72°C for 10 min. The PCR products were gel-purified with DNA Purification Kit (TianGen, Beijing, China), and ligated into a pMD®19-T vector (TaKaRa, Dalian, China), and then were transformed into DH5α Escherichia coli competent cells (TianGen, Beijing, China). These were then sequenced (Aoke Biotech Company, Xi’an, China) to obtain the correct sequence. Sequence similarities between GmolCSP8 and various homologs from different insects were calculated via NCBI BLASTX programs. Amino acid sequences were aligned by DNAMAN 6.0 software (Lynnon Biosoft, San Ramon, CA, United States). N-terminal signal peptides were predicted using the online SignalP 4.0 server1 (Petersen et al., 2011). The isoelectric point (pI) and molecular weight (Mw) of GmolCSP8 were calculated with the ExPASy server program “Compute pI/Mw”2 (Wilkins et al., 1999). Based on sequence alignment obtained via MAFFT 6.0 software (CBRC, Misawa, Miyata, Japan) (Kuraku et al., 2013), a phylogenetic tree was constructed by Seaview 4.0 software (Prabi, Paris, French) (Gouy et al., 2010) using the neighbor-joining method with a p-distance model and pairwise gap deletion. Bootstrapping was performed to estimate the reliability of branches using 1,000 replicates. GenBank accession numbers of 104 CSPs and corresponding amino acid sequences used in the phylogenetic analyses are listed in Supplementary Data Sheet S1.
Table 1. Gene-specific primers used for tissue expression analyses and site-directed mutagenesis key amino acid sites of GmolCSP8.
Expression Analysis Using qRT-PCR
The relative transcriptional levels of GmolCSP8 in different tissues were assessed by using qRT-PCR on a CFX96 Real-Time PCR Detection System (Bio-Rad, United States). All qRT-PCR assays were conducted according to the Minimum Information for Publication of Quantitative Real-Time PCR Experiments (MIQE) Guidelines (Bustin et al., 2009). The reference gene β-actin (GenBank No: KF022227.1) of G. molesta was used as the endogenous control to normalize target gene expression and correct for samples. The stability of β-actin gene expression had been verified in different temperatures and tissues (Chen et al., 2014, 2018). Specific-primers for target and reference genes were designed using the online program Primer3-blast3 (Ye et al., 2012) and are listed in Table 1. The amplification efficiencies of target and reference genes were validated by analyzing standard curves with a 5-fold cDNA serial dilution of a female antennae template (three replicates), with resultant efficiency of 90–105%. The specificity of each primer was assessed by melting-curve analysis, which confirmed that no non-specific product was produced. To check reproducibility, each sample was analyzed using three technical and three biological replicates. Each qPCR was performed in a 20 μL reaction volume contained 10 μL of 2 × SYBR® Premix Ex TaqTM II mixture (TaKaRa, Dalian, China), 0.8 μL of each primer (10 μM), 2.0 μL of sample cDNA, and 6.4 μL of sterilized ddH2O. Samples without template cDNA served as negative controls. The qPCR cycling conditions consisted of an initial denaturing cycle of 95°C for 30 s, followed by 40 cycles of 95°C for 10 s, 60°C for 30 s and 72°C for 30 s. The relative expression of GmolCSP8 was quantified using the comparative 2-ΔΔCt method (Livak and Schmittgen, 2001). To estimate the relative change in expression between different tissues, the abdomen sample was used as a calibrator for comparison between tissues. Significant differences in expression between tissues were assessed through one-way analysis of variance (ANOVA) (Tukey’s HSD tests) with a critical level of α = 0.05. Differences between male and female moths were measured by a paired t-test. All data were analyzed by using the SPSS 18.0 software (SPSS Inc., Chicago, IL, United States).
Production and Purification of rGmolCSP8
Specific primers with a Bam HI restriction enzyme site in the sense primer and a Hind III site in the antisense primer (Table 1) were used to amplify cDNA encoding the mature GmolCSP8 (without signal peptide). The PCR product was first ligated into the pMD®19-T vector (TaKaRa Co., Dalian, China) by T/A Pairing, and the resulting construct was transformed into DH5α competent cells. Positive clones were grown in LB liquid medium with ampicillin (10 mg/mL) and confirmed by sequencing. pMD®19-T plasmids containing the target sequences were digested with Bam HI and Hind III restriction enzymes for 4 h at 37°C. The digestion products were gel-purified with a DNA Purification Kit (TianGen, Beijing, China) on agarose gels, and then ligated into the pET28a (+) (Novagen, Madison, WI, United States) expression vector with T4 DNA Ligase (Thermo Fisher Scientific, Waltham, MA, United States), previously linearized with the same restriction endonucleases. The correct plasmids containing GmolCSP8 genes were transferred to BL21 (DE3) competent cells for overexpression. Single colonies were grown in LB liquid medium containing 50 μg/mL of kanamycin with shaking at 220 rpm for 14 h at 37°C. The culture was diluted 1:100 with fresh LB medium and incubated under the same conditions until the OD600 reached approximately 0.6–0.8. Expression was induced by adding 100 mM stock solution of isopropyl-β-D-thiogalactoside (IPTG) to a final concentration of 0.5 mM for an additional 10 h at 28°C and 160 rpm. The bacterial cells were harvested by centrifugation at 8,000 g for 10 min and resuspended in lysis buffer (20 mM Tris–HCl pH 7.4, 1 mM phenylmethanesulfonyl fluoride, and 250 mM NaCl). Bacterial cells were lysed with lysozyme (10 mg/mL) for 30 min at room temperature. The cultures was sonicated and centrifuged again at 13000 g for 30 min. SDS-PAGE analysis revealed the presence of soluble rGmolCSP8, and the supernatant was enriched by complete His-Tag Purification Resin (Roche, Mannheim, Germany). Purified rGmolCSP8 was desalted through extensive dialysis with 20 mM Tris–HCl buffer (pH 7.4) until the imidazole was eliminated. The fusion tag of the pET28a (+) expression vector, which is located on the N-terminus of GmolCSP8, was not removed owing to its small size and its lack of influence on the binding-ligand ability (Chen et al., 2018; Cui et al., 2018). Concentration of soluble protein was quantified as recommended in the BCA Protein Assay Kit (Beyotime, Shanghai, China). Purified rGmolCSP8 was identified an anti-His tag monoclonal antibody (Catalog # CW0286) (Cwbio, Beijing, China) via western blot analysis, as previously described (Chen et al., 2018; Li et al., 2018).
Fluorescence Competitive Binding Assays
To determine whether GmolCSP8 is involved in odorant binding, fluorescence competitive binding assays were performed (Zhong et al., 2012). The fluorescence spectra were recorded using a fluorescence spectrophotometer (F-4500, Hitachi, Tokyo, Japan) in a 1 cm light path fluorimeter quartz cuvette with 10 nm slits for excitation and 20 nm slits for emission. The selected chemicals are listed in Table 2, and included four sex pheromone components and 34 volatiles from different host plants (peach and pear trees) of G. molesta based on previous investigations (Natale et al., 2004; Lu et al., 2012; Najar-Rodriguez et al., 2013b). The fluorescent probe N-phenyl-1-naphthylamine (1-NPN) and all putative ligands used in binding assays were dissolved in methanol (HPLC purity grade) to obtain a 1 mM stock solution. The rGmolCSP8 protein was diluted to in 20 mM Tris–HCl to obtain a 2 μM solution, pH 7.4. To measure the binding affinity of 1-NPN for GmolCSP8, a 2 μM solution of protein was titrated with aliquots of 1 mM 1-NPN prepared to final concentrations ranging from 1 to 22 μM. The fluorescent probe was excited at 337 nm and emission spectra were recorded from 370 to 550 nm. The affinities of selected ligands were measured with the protein and 1-NPN at 2 μM, and the sex pheromone components were added to final concentrations of 0, 1, 2, 4, 6, 8, 10, and 12 μM, and of synthetic plant volatiles to final concentrations of 0, 2, 4, 6, 8, 12, 16, and 24 μM. Binding data were obtained from three replicates assays for each concentration of each tested ligand. Dissociation constants (Kd) of protein to 1-NPN were calculated by non-linear regression, using GraphPad Prism 4.00 software (GraphPad Software Inc.). The inhibition constants (Ki) of the competitors were determined using the corresponding IC50 values (concentration of competitor that reduced the initial fluorescence intensity by half) based on the following equation: Ki = [IC50]/(1+[1-NPN]/K1-NPN), where [1-NPN] is the free concentration of 1-NPN and K1-NPN is the dissociation constant of the complex GmolCSP8/1-NPN. In this study, if the IC50 value was less than 10 μM, the tested ligands were considered to exhibit strong binding affinities to rGmolCSP8; an IC50 value between 10 and 20 μM indicated that ligands demonstrated medium binding ability to protein, an IC50 exceeding 20 μM indicated that ligands have a weak binding capacity to protein. If the IC50 value exceeded 24 μM, no further calculation of the binding affinity was considered.
Table 2. Binding affinities of rGmolCSP8 to sex pheromone and host-plant volatile compounds were detected via fluorescence binding assays.
Effect of pH on the Binding Affinity of rGmolCSP8
Based on the fluorescent binding assays described above, we further assessed the effects of pH on the binding abilities of GmolCSP8 with strongly bound ligand. The binding affinities of rGmolCSP8 to the representative ligand 1- hexanol were measured at six different pH values (4.0, 5.0, 6.8, 7.4, 8.0, and 8.8) in 20 mM Tris–HCl buffer. pH was adjusted by titrating concentrated hydrochloric acid using a PHS-3E pH meter (Shanghai Precision Scientific Instrument Co., Ltd., Shanghai, China). The binding assays and data processing were performed as described previously.
Three-Dimensional (3D) Structure Modeling and Molecular Docking of Ligand
A protein blast search was performed for the GmolCSP8 amino acid sequence against the Protein Data Bank (PDB) to identify structural templates for homology modeling. Based on sequence similarity and coverage of the target and template proteins (Supplementary Table S1), we selected the crystal structure of CSPsg4 from the desert locust Schistocerca gregaria (PDB ID: 2gvs.1) as a template to generate the 3D structure of GmolCSP8 using Modeller 9.10 software (Gunda et al., 2016) with an automated approach. A 1-hexanol model was generated using ChemBioDraw12.0 (Cousins, 2011), and the three-dimensional structure was preserved after optimization via MM2 energy minimization. The molecular docking of GmolCSP8 binding 1-hexanol was performed using Autodock 4.0 software tool as the following methods: first, the structures of 1-hexanol and 3D models of GmolCSP8 were optimized by adding hydrogens using Kollman charges. Then, the model was prepared by optimizing torsion angles and saved in PDBQT format. Finally, flexible docking of GmolCSP8 with 1-hexanol was performed by Autodock software, without considering the effect of water and other solvents. Amino acids of 3Å around 1-hexanol were considered to be important residues of the binding ligand.
Site-Directed Mutagenesis and Expression of Mutants
To verify the four predicted key amino acid sites of GmolCSP8 binding for 1-hexanol, recombinant plasmids for four CSP8 mutants were obtained via site-directed mutagenesis experiments. Specific primers introducing mutation sites of the four key amino acid CSP8 plasmids (GmolCSP8-T27A, GmolCSP8-L30A, GmolCSP8-V40A, and GmolCSP8-L80A) were designed and are listed in Table 1. The four mutants were mutated by using the TaKaRa MutanBEST Kit (TaKaRa, Shiga, Japan) with the pMD®19-T/GmolCSP8 plasmid DNA as template. All plasmids containing mutations were validated by sequencing with five randomly selected single colonies. pMD®19-T plasmids including four CSP8 mutants were digested with Bam HI and Hind III restriction enzymes. Then, the target sequences were ligated into the expression vector pET28a (+) and transformed into BL21 (DE3) competent cells for expression. After sequencing verification, the corresponding rGmolCSP8 mutant proteins were induced and purified, and the binding affinities to 1-hexanol and its analogs were verified via ligand-binding assays according to the method described above.
Results
Cloning and Phylogenetic Analysis of GmolCSP8
The full-length cDNA of GmolCSP8 was cloned and then confirmed by sequencing. Three incorrect basic groups were identified in the Unigene of GmolCSP8, obtained by antennal transcriptome sequencing, and the verified sequence was deposited in GenBank under accession number KR003781. GmolCSP8, encoded by 119 amino acids with an open reading frame (ORF) of 360 bp, was found to contain a predicted signal peptide of 15 amino acid residues at the N-terminus (Figure 1A). The calculated molecular weight (MW) and isoelectric point (pI) of mature GmolCSP8 is 13.53 kDa and 8.22, respectively.
Figure 1. Cloning and phylogenetic analysis of GmolCSP8 from Grapholita molesta. (A) Nucleotide and deduced amino acid sequence of GmolCSP8 from G. molesta. The four conserved cysteines are marked by an amaranthine circle. The predicted signal peptide sequence is underlined. (B) Sequence alignment of GmolCSP8 with other insect species. The four conserved cysteines are marked with five-pointed red star. The names and the GenBank accession number of the 15 CSPs are as follows: GmolCSP8, Grapholita molesta, ALC79595.1; LbotCSP8.1, Lobesia botrana, AXF48705.1; OfurCSP6, Ostrinia furnacalis, BAV56810.1; SlitCSP12, Spodoptera litura, ALJ30223.1; CmedCSP31, Cnaphalocrocis medinalis, ALT31613.1; DkikCSP, Dendrolimus kikuchii, AII01029.1; DpunCSP10, Dendrolimus punctatus, ARO70314.1; CpunCSP8, Conogethes punctiferalis, AHX37222.1; EoblCSP20, Ectropis obliqua, ALS03845.1; DhouCSP, Dendrolimus houi, AII01021.1; AipsCSP4, Agrotis ipsilon, AGR39574.1; PxutCSP, Papilio xuthus, BAF91717.1; CsupCSP, Chilo suppressalis, AHC05674.1; CmedCSP, Cnaphalocrocis medinalis, AIX97832.1; BmorCSP13, Bombyx mori, NP_001037180.1.
Moreover, multiple sequence alignment of GmolCSP8 with homologous CSPs from other Lepidopteran species revealed the presence of a typical conserved motif containing four cysteines (C1-X6-8-C2-X16-21-C3-X2-C4, X represent any amino acids except cysteine) (Figure 1B). Furthermore, using the NCBI BLASTx program, we found that GmolCSP8 shared relatively high sequence similarity with other insect CSPs, with 82% identity with Lobesia botrana LbotCSP8 (GenBank No. AXF48705.1), and almost 60% identities with Ostrinia furnacalis OfurCSP6 (GenBank No. BAV56810.1), Spodoptera litura SlitCSP12 (GenBank No. ALJ30223.1), Cnaphalocrocis medinalis CmedCSP31, Dendrolimus punctatus DpunCSP10 (GenBank No. ARO70314.1), Conogethes punctiferalis CpunCSP8 (GenBank No. AHX37222.1), and Ectropis obliqua EoblCSP20 (GenBank No. ALS03845.1). GmolCSP8 shared relatively low sequence identity (less than 50%) with the other 14 CSPs from intraspecific species (Supplementary Table S2). Subsequently, 104 CSPs obtained from different order insects comprised of Graphilita molesta (15 sequences), Laodelphax striatella (12 sequences), Aedes aegypti (43 sequences), Bombyx mori (16 sequences) and Tribolium castaneum (18 sequences) were selected and used to construct phylogenetic tree. Phylogenetic analysis showed that GmolCSP8 was clustered into a small branch close to BmorCSP13 from B. mori and TcasCSP6 from T. castaneum (Figure 2). The majority of CSPs of A. aegypti and T. castaneum clustered into the same branch with intraspecific CSPs, while the CSPs of G. molesta, B. mori, and L. striatella were highly divergent and clustered into dozens of branches of the in Neighbor-joining tree (Figure 2).
Figure 2. Phylogenetic tree of CSPs from Grapholita molesta and other insects through the Neighbor-joining method with a bootstrap replication of 1,000. The GenBank accession number and the corresponding amino acid sequences used in the phylogenetic analyses were listed in Supplementary Data Sheet S1. The unrooted tree was constructed using the BioNJ algorithm in Seaview v.4.0 based on the sequence alignment produced via MAFFT version 6. Gmol, Graphilita molesta; Lstr, Laodelphax striatella; Cqui, Culex quinquefasciatus; Bmor, Bombyx mori; Tcas, Tribolium castaneum.
Tissue Expression of GmolCSP8
A pilot experiment showed that the amplification efficiencies of target and reference genes were between 90 and 105%, with no non-specific amplification (Supplementary Figure S1), confirming that the two primer pairs and selected amplification conditions could be used for subsequent quantification. The qRT-PCR results showed that GmolCSP8 was predominantly expressed in the wings and antennae of both males and females, but lowly expressed in abdomen and legs of both males and females, as well as in heads and thoraxes of males (Figure 3). In addition, the expression of GmolCSP8 in the wings, heads, thoraxes and abdomen differed significantly between males and females, with 5.60- and 23.48-fold higher expression in the heads and thoraxes of female than those of male. It is speculated that GmolCSP8 may participate in chemoreception and other physiological functions based on its expression characteristics.
Figure 3. Relative expression level of GmolCSP8 in different tissues of Grapholita molesta. An: antennae, He: heads, Th: thoraxes, Ab: abdomen, Le: legs, Wi: wings. The fold changes are relative to the expression levels in the abdomen. Different lowercase and capital letters above each bar indicate significantly differences (P < 0.05, ANOVA, Tukey) among different tissues in female and male, respectively. ∗Indicate significant differences in expression levels between two sexes in the same tissues (Independent t-test, α = 0.05).
Fluorescent Binding Assays
rGmolCSP8 was successfully expressed in soluble forms in a prokaryotic expression system with a high yield (about 0.8 mg/mL). Subsequently, the expression and purification of rGmolCSP8 was examined by 15% SDS-PAGE (Figure 4A) and western blot (Figure 4B), respectively. In addition, the binding affinities of GmolCSP8 to various ligands were determined by ligand-binding assays. Among those, binding of the fluorescent probe 1-NPN to GmolCSP8 was evaluated and a gradual saturation and linear Scatchard plot was obtained, yielding a Kd value of 1-NPN with rGmolCSP8 of 13.85 μM (Figure 5). In summary, 38 chemicals were selected as putative ligands to detect the binding characteristics of rGmolCSP8, which was comprised of four sex pheromones and 34 volatiles from peach and pear trees (Table 2). The results showed that rGmolCSP8 exhibited specific binding properties and could only bind to a few ligands, such as alcohols and terpenes. Of those, rGmolCSP8 showed medium binding affinity to the minor sex pheromone component 1-dodecanol (12:OH) with a Ki value of 12.90 μM. In addition, rGmolCSP8 displayed the strongest binding affinities to green leaf volatile 1-hexanol with a Ki value of 2.11 μM, as well as the other small aliphatic alcohols comprised of 1-butanol (6.65 μM), 1-pentanol (8.18 μM), 3-methyl-1-butanol (6.32 μM), cis-3-hexen-1-ol (8.51 μM), 1-heptanol (6.20 μM), and 1-decanol (6.60 μM), respectively (Figure 6). However, for α-ocimene, rGmolCSP8 showed a weak binding affinity with a Ki value of 17.58 μM. However, rGmolCSP8 did not bind to the tested aldehydes, esters, benzonitriles, or alkanes (Table 2).
Figure 4. SDS-PAGE (A) and Western blot (B) analysis of the expressed products and purified of rGmolCSP8. M1, M2: Protein molecular weight marker; 1. Non-induced pET28a(+)/GmolCSP8; 2. Induced pET28a(+)/GmolCSP8; 3. pET28a(+)/GmolCSP8 supernatant; 4. pET28a(+)/GmolCSP8 pellet; 5. Purified pET28a (+)/GmolCSP8; 6. Western Blot analysis of purified rGmolCSP8.
Figure 5. Binding curve and relative Scatchard plot of 1-NPN to rGmolCSP8. A 2 μM of rGmolCSP8 protein in 20 mM Tris–HCl buffer (pH 7.4) was titrated with aliquots of 1 mM 1-NPN stock solution to final concentrations of 0 to 22 μM, and the emission spectra were recorded between 370 and 550 nm. The calculated dissociation constant of rGmolCSP8 was 13.85 μM.
Figure 6. Binding curves of rGmolCSP8 with the small aliphatic alcohols. Mixtures of 2 μM rGmolCSP8 and 2 μM 1-NPN were titrated with 1 mM stock solution of each ligand to final concentrations of 0–24 μM. Relative fluorescence intensities are reported as percentages of the initial fluorescence values without competitors.
Effect of pH on the Binding Affinity of rGmolCSP8
To determine the effect of pH on the binding affinity of rGmolCSP8 to ligand, the binding abilities of GmolCSP8 in different pH buffered solutions (20 mM Tris–HCl) were evaluated with the representative ligand 1-hexanol in vitro. The resulting binding values indicated that pH significantly affected the binding affinities of rGmolCSP8 to 1-hexanol; the binding affinity was higher at pH 6.8–8.0 than at pH 4.0–5.0 and 8.8 (Figure 7), revealing that ligand binding and protein release were related to a pH-dependent conformational transition.
Figure 7. Effects of pH on the binding affinities of rGmolCSP8 with 1-hexanol. The mixture solution of 2 μM rGmolCSP8 protein and 2 μM 1-NPN was titrated with 1 mM 1-hexanol stock solution to final concentrations of 22 μM. All data represent a mean of three independent measurements, and the bars denote mean ± SEM.
Structural Modeling and Molecular Docking
Structural templates for GmolCSP8 from G. molesta were identified by sequence alignments. The results demonstrated that GmolCSP8 shared more than 30% sequence identities with four structurally determined CSPs, including Schistocerca gregaria CSPsg4 (2gvs.1), Mamestra brassicae MbraCSP2 (1k19.1.A), M. brassicae MbraCSPA6 (1kx9.1.A) and Bombyx mori BmorCSP1 (2jnt.1.A). Importantly, GmolCSP8 was found to share the highest similarity and coverage with CSPsg4, with values of 44.7 and 95.0%, respectively (Figure 8A). Finally, the 3D structure of GmolCSP8 was generated using the CSPsg4 structure as a template. The generated structure was validated by PROCHECK. The model presented 89.2% of amino acid residues in the core region, 8.6% of amino acid residues in the additionally allowed region, 1.1% of amino acid residues in the generously favored region, and 1.1% of amino acid residues in the disallowed region. A Ramachandran plot of the modeled protein is shown in Supplementary Figure S2.
Figure 8. Modeled three-dimensional structure and molecular docking of GmolCSP8. (A) Sequence alignment of GmolCSP8 and CSPsg4 (2gvs.1). α-helices are displayed as spiral coils. The same residues are highlighted with a red background. (B) Three-dimensional structure of GmolCSP8. Two disulfide bridges are colored in yellow. α-helices, N-terminal, and C-terminal are labeled. (C) Molecular docking of GmolCSP8 with 1-hexanol. 1-hexanol is displayed as a stick model with the hydroxyl oxygen in red. (D) The interaction diagram of amino acid residues to 1- hexanol. The top eight potential key residues (Tyr26, Thr27, Val40, Leu43, Lys44, Leu30, Val31, and Leu80) are shown with black font.
The simulated 3D structure of GmolCSP8 had an overall globular shape and consisted of six α-helices, represented as α1 (Leu13-Gly18), α2 (Glu20-Val31), α3 (Val40-Glu52), α4 (Asp60-Lys76), α5 (Pro78-Tyr88), and α6 (His95-Ala97). Of those, helices α1 and α2 (residues 13–31), together with helices α4 and α5 (residues 60–88) formed two paralleled V-shaped structures; and helix α3 was perpendicular to the planes determined by two V-shaped structures and positioned between the four ends of the two V-shaped structures (Figure 8B). In addition, helix α6 was located on the outer surface of the α4–α5 helices and did not participate in the formation of the binding pocket. Four conserved cysteines are linked, in turn, to form two pairs of disulfide bridges (connected Cys29 and Cys36, Cys55, and Cys58) to maintain the spatial structure of protein. The residue Cys29 is located on the helices, and the residues Cys36, Cys55, and Cys58 are located on loops.
Molecular flexible docking was conducted to predict the key binding sites of GmolCSP8 to 1-hexanol (Figure 8C). Interactions between functional groups of ligand and binding sites of GmolCSP8 are shown in Figure 8D. The docking simulation showed that Tyr26, Thr27, Val40, Leu43, Lys44, Leu30, Val31, and Leu80 may be responsible for the binding to 1-hexanol and Thr27 might be involved in specific recognition via H-bonding with hydroxyl oxygen of 1-hexanol. To determine whether each predicted residue had an important role in 1-hexanol binding, four residues (Thr27, Leu30, Val40, and Leu80) with lower interaction energy to 1-hexanol were mutated to alanine using site-directed mutagenesis. Then, the functions of putative key residues were validated by measuring the binding affinity of each mutant protein to 1-hexanol.
Site-Directed Mutagenesis of GmolCSP8 and Binding Affinities of Mutants
Using site-directed mutagenesis, four predicted key amino acid residues (Thr27, Leu30, Val40, and Leu80) of GmolCSP8 were substituted with alanine (Ala), yielding the four mutants T27A, L30A, V40A, and L80A (Figure 9). There was no significant difference in the calculated dissociation constant (Kd) of the four mutants to the fluorescence probe 1-NPN (Supplementary Figure S3). To evaluate the functions of these residues in the binding process, the binding affinities of mutant CSP8-T27A, CSP8-L30A, CSP8-V40A, and CSP8-L80A proteins to 1-hexanol and its analogs (Figure 10F) were measured using binding assays (Table 3). Compared with wild-type rGmolCSP8, the ability of CSP8-T27A to bind to 1-hexanol, 1-butanol, and 1-pentanol were significantly reduced and did not displace 50% 1-NPN from the CSP8-T27A/1-NPN complex solution at concentrations of ligands up to 24 μM. It also demonstrated a 1.5- to 2.0-fold decrease toward 1-heptanol and 1-decanol, indicating the Thr27 was one of the key amino acid residues of GmolCSP8 involved in the binding to 1-hexanol, Thr27 may engage in H-bond with this ligand (Figures 10A,B). Mutant CSP8-L30A lost its ability to bind to 1-pentanol, 1-heptanol, and 1-decanol and demonstrated a 6-fold decrease in its ability to bind 1-hexanol (Figure 10C), revealing that Leu30 is an important residue in the binding site of GmolCSP8 for 1-hexanol and its homologs. The mutants CSP8-V40A and CSP8-L80A showed only a slightly change in binding to 1-hexanol and its analogs (Figures 10D,E); a possible explanation is that Val40 and Leu80 do not participate in recognition of the small aliphatic alcohols.
Figure 9. SDS-PAGE analysis of rGmolCSP8 wild-type (WT) and mutants T27A, L30A, V40A, and L80A. M represents protein molecular weight marker.
Table 3. Binding affinities of 1-hexanol and its analogs to rGmolCSP8 Wild-Type and mutants.
Discussion
Analysis of the tissue expression profiles of CSP genes can help to understand their function in insects. GmolCSP8 was highly expressed in antennae, suggesting potential roles in olfactory chemoreception. Antennae are the most important olfactory organs in insects; they function as a switch that regulates the entry of volatile semiochemicals from the outside environment into the olfactory system (Gao and Chess, 1999; Leal, 2013). Several studies have shown that antenna-specific or -highly expressed CSPs can selectively perceive and recognize sex pheromones and general odorants (Calvello et al., 2005; Peng et al., 2017). In insects, antennae, wings, and legs belong to contact organs while the head, thorax, and abdomen are classified as non-contact organs (Gong et al., 2007;He et al., 2011; Hua et al., 2013). GmolCSP8 was also highly expressed in wings; the relative expression in female wings was significantly higher than that in male wings, indicating that it may be involved in contact chemoreception, especially for females. The biased or high expression of CSPs in female wings has also been found in other insects; for example, Apolygus lucorum AlucCSP2 and AlucCSP3 were specifically expressed in female wings (Hua et al., 2013), Nilaparvata lugens NlugCSP5 was biased and highly expressed in female wings (Yang et al., 2014). The transcript level of GmolCSP8 in the abdomen of female and male adults was negligible; combined with the recombinant protein not binding the major sex pheromones; we speculated that this gene may not participate in the transport of sex pheromones in the gonads (Li Z.Q. et al., 2015). The ubiquitously expression of GmolCSP8 in G. molesta indicates that this gene may be involved in functions other than chemosensing (Guo et al., 2011; Ali et al., 2019).
To elucidate the olfactory function of GmolCSP8, 38 compounds, including four sex pheromones and 34 plant volatiles emitted from peach and pear trees, were selected for ligand-binding assays using 20 mM Tris–HCl (at pH 7.4) as a buffer solution (Cardé et al., 1975, 1979; Natale et al., 2003; Najar-Rodriguez et al., 2013b; Lu et al., 2015). GmolCSP8 exhibited specific-binding affinities to green leaf volatile 1-hexanol, as well as its analogs 1-butanol, 3-methyl-1-butanol, 1-heptanol, 1-decanol, and 1-dodecanol. Conversely, others behaviorally important ligands such as (Z)-3-hexenyl acetate, (E)-2-hexenal, hexanal, heptanal, benzaldehyde, butyl hexanoate, α-pinene and benzonitrile, cannot bind to GmolCSP8. For example, a three-compound mixture of (Z)-3-hexenyl acetate, (Z)-3-hexen-1-ol and benzaldehyde in a 4:1:1 proportion elicited a similar attractant effect on female G. molesta as the natural blend from peach shoot volatiles (Natale et al., 2003); butyl hexanoate was clearly attractive for virgin and mated female G. molesta at doses of 0.32–1.35 mg (Natale et al., 2004). Five aldehydes, including hexanal, heptanal, octanal, and nonanal, and benzaldehyde, play essential roles in female attraction to host plants (Najar-Rodriguez et al., 2013b). The special binding property with small linear aliphatic alcohols suggested that GmolCSP8 may possess a suitable binding site for ligands with a hydroxyl functional group. The length of the carbon chain significantly affects the affinity of protein with ligands; GmolCSP8 exhibited the strongest affinity to 1-hexanol (C6) with a Ki value of 2.11 μM, and also displayed relative higher binding abilities to other aliphatic alcohols with 4–12 carbon atoms with a Ki value of 6.60–12.90 μM. It is unable to bind to 1-tetradecanol when the number of carbon atoms in the ligand increases to 14. 1-dodecanol is a minor sex pheromone component of G. molesta, and functions as a synergist of sex attractant, which increased the frequency of male landing and as a stimulant that induces the occurrence of mating behavior when adult moths of both sexes are in close proximity (Cardé et al., 1975, 1979). The results of ligand-binding assays showed that GmolCSP8 has medium binding abilities to 1-dodecanol with a Ki value of 12.90 μM in vitro. The binding affinities of this protein to 1-dodecanol were significantly lower than those of GmolOBPs to 1-dodecanol (e.g., the Ki values of GmolGOBP2, GmolOBP8, and GmolOBP15 to 1-dodecanol were 1.22, 4.39, and 3.24 μM, respectively) (Li et al., 2016a,b). Further validation combining gene knockout experiments and behavioral assays is required to confirm whether GmolCSP8 is involved in the perception and recognition of 1-dodecanol.
Figure 10. Binding curves of Wild-Type GmolCSP8 and select mutants. (A) Wild-Type GmolCSP8, (B) mutant T27A, (C) mutant L30A, (D) mutant V40A, (E) mutant L80A, and (F) ligand compounds and their structural formulas. The mixture solution of 2 μM each recombinant protein and 2 μM 1-NPN was titrated with 1 mM ligand stock solution to final concentrations of 2 μM. Relative fluorescence intensities are reported as percentages of the initial fluorescence values without competitors.
The binding values of GmolCSP8 to its preferred ligand 1-hexanol under different pH environments revealed that a pH-dependent conformational transition is involved in binding and releasing protein to ligands. GmolCSP8 displayed a strong binding affinity to 1-hexanol at pH 6.8–8.0, while it showed a very weak affinity at pH 4.0. The binding and release mechanism of GmolCSP8 was similar to that of BmorPBP1 and BmorGOBP2 of Bombyx mori, ApolPBP of Antheraea polyphemus, and AtraPBP1 of Amyelois transitella (Mohanty et al., 2004; Zhou et al., 2009; Xu et al., 2010). The longer C-terminal tail of BmorPBP1 can form an additional helix, α7, which is inserted between the first and second α-helix and occupies a bombykol-binding site at low pH resulting in ligand release, The helix α7 converted to an unfolded conformation under a neutral pH value and exhibited ligand-binding activity (Sandler et al., 2000). Another explanation is that the protonation of sidechains causes a general unraveling of the BmorPBP1 conformation as the pH is reduced, and the loss of secondary structure resulted in the opening of the binding pocket to release sex pheromones (Lartigue et al., 2002). Several CSPs have demonstrated that a pH-induced conformational change triggers ligand-binding and -release, including MaltCSP5 from Monochamus alternatus and MsepCSP8 from Mythimna separate (Younas et al., 2018; Ali et al., 2019). The α-helical content of CmedCSP33 was decreased at low pH values, suggesting the protein conformation loosened, which may facilitate odorant binding (Duan et al., 2019). However, the pH-dependent mechanism of CSPs in the binding and release of odorants remains unclear. Compared with the three-dimensional structure of OBPs, five helices (α1, α2, α4, α5, and α3) in CSPs form narrower and elongated channels, which conveniently bind aliphatic ligands. In contrast to the longer C-terminal of OBPs, CSPs possess a disordered and long N-terminal; whether the N-terminal peptide segment can form a regular helix at some pH values and trigger binding and release in a comparable way to the C-terminal helix of OBP remains to be explored.
Chemosensory proteins are important odorant recognition proteins during the initial process of olfactory perception, and the detailed interactions of these proteins with specific odor molecules have been elucidated via site-directed mutagenesis (Tan et al., 2018). Several predicted key binding site residues of CSPs form a hydrogen bond, including Tyr26 of CSPMbraA6 with 12-bromo-dodecanol (Lartigue et al., 2002), Trp83 of CSPsg4 with oleamide (Tomaselli et al., 2006), and Cys 60 of Ac-ASP3 with isoamyl acetate (Liu et al., 2013). Some insect CSPs use van der Waals forces to bind odor molecules, such as Phe29 of CmedCSP33 with nerolidol (Duan et al., 2019), and F44 of AmelCSP1 to ionone (Tan et al., 2018). Compared with the wild type rGmolCSP8 protein, the T27A mutation abolished the ability to bind 1-hexanol. The most credible interpretation is that the oxygen atom of the hydroxyl group in 1-hexanol forms a hydrogen bond with Thr27; this ligand cannot be recognized by the mutant protein due to the loss of hydrogen bonding. The L30A mutant showed an approximate 6-fold decreased in binding to 1-hexanol; we speculated that the Van der Waals force formed by Leu30 of GmolCSP8 and 1-hexanol, and possibly by Thr27 and Leu30, combined to recognize 1-hexanol. V40A and L80A mutants revealed no significant change in the binding affinities to 1-hexanol and its analogs, suggesting that Val40 and L80A are not involved in the recognition of these ligands. Overall, Thr27 is the key binding site of GmolCSP8, involved in the initial recognition of 1-hexanol, and Leu30 also plays an important role in binding to this ligand.
The selective binding and transport of external volatiles by insects is accomplished through the synergistic action of various OBPs and CSPs (Liu et al., 2012). 1-hexanol is the main green leaf volatile in peach shoots and pear fruits (Piñero and Dorn, 2007; Najar-Rodriguez et al., 2013b). Electrophysiological experiments have shown that 1-hexanol can elicit strong EAG responses in male and female antennae of G. molesta (Li et al., 2016b). 1-hexanol is also a key component of a mixture (1-hexanol, nonanal, ethyl butanoate, butyl acetate, 3-methylbutyl acetate, ethyl hexanoate, and hexyl acetate, ratio 3:1:100:8:4:12:4) with attractive activity to adult G. molesta (Lu et al., 2014). Preliminary evidence from fluorescence competitive binding assays confirmed that GmolCSP8 participates in the recognition and transport of 1-hexanol. GmolCSP8 may serve as a target gene for G. molesta host-locating behavior.
Author Contributions
G-WL, X-LC, W-QW, and J-XW conceived and designed the experiments. G-WL and X-LC performed the experiments. L-HC, G-WL, and W-QW analyzed and processed data. G-WL and J-XW wrote and edited the manuscript. All authors agreed to publish this manuscript.
Funding
This study was supported and funded by the National Natural Science Foundation of China (Grant Nos. 31272043 and 31860506), The Project of Industry-University-Research Collaboration of Yan’an University in 2017 (Grant No. 2017cxy04), and Ph.D. Research Startup Foundation of Yan’an University (Grant No. YDBK2017-37).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors appreciated Dr. Jian Gao (College of Advanced Agriculture and Bioengineering, Yangtze Normal University, Fuling, China) for his suggestion on our manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2019.00552/full#supplementary-material
FIGURE S1 | Melting curves and Standard curves of reference and target genes in qRT-PCR. (A) and (B) are the melting curves of Gmolβ-actin and GmolCSP8 genes, respectively. (C,D) represents the standard curves of Gmolβ-actin and GmolCSP8 genes, respectively.
FIGURE S2 | Ramachandran plot of predicted of GmolCSP8 (KR003781) of Grapholita molesta.
FIGURE S3 | Binding curves of 1-NPN to GmolCSP8 Wild-Type and its four mutants. rGmolCSP8. A 2 μM of each protein in 20 mM Tris–HCl buffer (pH 7.4) was titrated with aliquots of 1 mM 1-NPN stock solution to final concentrations of 0 μM to 24 μM, and the emission spectra were recorded between 370 and 550 nm. The calculated dissociation constants of rGmolCSP8 Wild-Type, T27A, L30A, V40A, and L80A were 13.85, 11.19, 14.38, 11.82, and 11.71 μM, respectively.
TABLE S1 | Similarity and coverage of GmolCSP8 with other known crystal structure insect CSPs.
TABLE S2 | Sequence identity analysis of GmolCSP8 and 14 other CSPs of Grapholita molesta.
DATA SHEET S1 | Amino acid sequences of CSPs were used in phylogenetic analyses.
Footnotes
- ^ http://www.cbs.dtu.dk/services/SignalP-4.0/
- ^ https://web.expasy.org/compute_pi/
- ^ https://www.ncbi.nlm.nih.gov/tools/primer-blast
References
Ali, S., Ahmed, M. Z., Li, N., Ali, S. A. I., and Wang, M. Q. (2019). Functional characteristics of chemosensory proteins in the sawyer beetle Monochamus alternatus hope. Bull. Entomol. Res. 109, 34–42. doi: 10.1017/S0007485318000123
Antony, B., Soffan, A., Jakše, J., Abdelazim, M. M., Aldosari, S. A., Aldawood, A. S., et al. (2016). Identification of the genes involved in odorant reception and detection in the palm weevil Rhynchophorus ferrugineus, an important quarantine pest, by antennal transcriptome analysis. BMC Genomics 17:69. doi: 10.1186/s12864-016-2362-6
Bruyne, M., and Baker, T. C. (2008). Odor detection in insects: volatile codes. J. Chem. Ecol. 34, 882–897. doi: 10.1007/s10886-008-9485-4
Bustin, S. A., Benes, V., Garson, J. A., Hellemans, J., Huggett, J., Kubista, M., et al. (2009). The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin. Chem. 55, 611–622. doi: 10.1373/clinchem.2008.112797
Calvello, M., Brandazza, A., Navarrini, A., Dani, F. R., Turillazzi, S., Felicioli, A., et al. (2005). Expression of odorant-binding proteins and chemosensory proteins in some hymenoptera. Insect Biochem. Mol. Biol. 35, 297–307. doi: 10.1016/j.ibmb.2005.01.002
Cardé, A. M., Baker, T. C., and Cardé, R. T. (1979). Identification of a four-component sex pheromone of the female oriental fruit moth, Grapholita molesta (Lepidoptera: Tortricidae). J. Chem. Ecol. 5, 423–427. doi: 10.1007/BF00987927
Cardé, R. T., Baker, T. C., and Roelofs, W. L. (1975). Behavioural role of individual components of a multichemical attractant system in the oriental fruit moth. Nature 253, 348–349. doi: 10.1038/253348a0
Chen, H., Xu, X. L., Li, Y. P., and Wu, J. X. (2014). Characterization of heat shock protein 90, 70 and their transcriptional expression patterns on high temperature in adult of Grapholita molesta (Busck). Insect Sci. 21, 439–448. doi: 10.1111/1744-7917.12057
Chen, X. L., Su, L., Li, B. L., Li, G. W., and Wu, J. X. (2018). Molecular and functional characterization of three odorant binding proteins from the oriental fruit moth Grapholita molesta (Busck) (Lepidoptera: Tortricide). Arch. Insect Biochem. Physiol. 98:e21456. doi: 10.1002/arch.21456
Claridge-Chang, A., Wijnen, H., Naef, F., Boothroyd, C., Rajewsky, N., and Young, M. W. (2001). Circadian regulation of gene expression systems in the Drosophila head. Neuron 32, 657–671. doi: 10.1016/S0896-6273(01)00515-3
Cousins, K. R. (2011). Computer review of chemdraw ultra 12.0. J. Am. Chem. Soc. 133:8388. doi: 10.1021/ja204075s
Cui, X. N., Liu, D. G., Sun, K. K., He, Y., and Shi, X. Q. (2018). Expression profiles and functional characterization of two odorant-binding proteins from the apple buprestid beetle Agrilus mali (Coleoptera: Buprestidae). J. Econ. Entomol. 111, 1420–1432. doi: 10.1093/jee/toy066
de Morais, R. M., Redaelli, L. R., and Sant’Ana, J. (2012). Age and multiple mating effects on reproductive success of Grapholita molesta (Busck) (Lepidoptera, Tortricidae). Rev. Bras. Entomol. 56, 319–324. doi: 10.1590/S0085-56262012000300009
Du, J., Wang, Y. R., and Wu, J. X. (2010). Effect of four different artificial diets on development and reproduction of Grapholita molesta (Lepidoptera: Tortricidae). J. Shanxi Agric. Univ. 30, 229–231. doi: 10.13842/j.cnki.issn1671-8151.2010.03.015
Duan, S. G., Li, D. Z., and Wang, M. Q. (2019). Chemosensory protein used as target for screening behaviourally active compounds in rice pest Cnaphalocrocis medinalis (Lepidoptera: Pyralidae). Insect Mol. Biol. 28, 123–135. doi: 10.1111/imb.12532
Dyanov, H., and Dzitoeva, S. (1995). Method for attachment of microscopic preparations on glass for in situ hybridization, PRINS and in situ PCR studies. Biotechniques 18, 822–826. doi: 10.1021/ic020451f
Gao, Q., and Chess, A. (1999). Identification of candidate Drosophila olfactory receptors from genomic DNA sequence. Genomics 60, 31–39. doi: 10.1006/geno.1999.5894
Gong, D. P., Zhang, H. J., Zhao, P., Lin, Y., Xia, Q. Y., and Xiang, Z. H. (2007). Identification and expression pattern of the chemosensory protein gene family in the silkworm. bombyx mori. Insect Biochem. Mol. Biol. 37, 266–277. doi: 10.1016/j.ibmb.2006.11.012
Gouy, M., Guindon, S., and Gascuel, O. (2010). SeaView version 4: A multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol. Biol. Evol. 27, 221–224. doi: 10.1093/molbev/msp259
Gunda, S. K., Chary Rudroju, S., Bandi, S., and Bodhankar, R. (2016). Homology modeling of human asparaginase like protein and its molecular interaction studies with natural flavonoids. Int. J. Pharm. Sci. Res. 7, 287–291. doi: 10.13040/IJPSR.0975-8232.7(1).287-91
Guo, W., Wang, X., Ma, Z., Xue, L., Han, J., Yu, D., et al. (2011). CSP and takeout genes modulate the switch between attraction and repulsion during behavioral phase change in the migratory locust. PLoS Genet. 7:e1001291. doi: 10.1371/journal.pgen
He, P., Zhang, J., Liu, N. Y., Zhang, Y. N., Yang, K., and Dong, S. L. (2011). Distinct expression profiles and different functions of odorant binding proteins in Nilaparvata lugens Stål. PLoS One 6:e28921. doi: 10.1371/journal.pone.0028921
Hua, J. F., Zhang, S., Cui, J. J., Wang, D. J., Wang, C. Y., Luo, J. Y., et al. (2013). Functional characterizations of one odorant binding protein and three chemosensory proteins from Apolygus lucorum (Meyer-Dur) (Hemiptera: Miridae) legs. J. Insect Physiol. 59, 690–696. doi: 10.1016/j.jinsphys.2013.04.013
Il’ichev, A. L., Kugimiya, S., Williams, D. G., and Takabayashi, J. (2009). Volatile compounds from young peach shoots attract males of oriental fruit moth in the field. J. Plant Interact. 4, 289–294. doi: 10.1080/17429140903267814
Jacquin-Joly, E., Vogt, R. G., Francois, M. C., and Nagnan-Le, M. P. (2001). Functional and expression pattern analysis of chemosensory proteins expressed in antennae and pheromonal gland of Mamestra brassicae. Chem. Senses 26, 833–844. doi: 10.1093/chemse/26.7.833
Jin, X., Brandazza, A., Navarrini, A., Ban, L., Zhang, S., Steinbrecht, R. A., et al. (2005). Expression and immunolocalisation of odorant-binding and chemosensory proteins in locusts. Cell. Mol. Life Sci. 62, 1156–1166. doi: 10.1007/s00018-005-5014-6
Kitabayashi, A. N., Arai, T., Kubo, T., and Natori, S. (1998). Molecular cloning of cDNA for p10, a novel protein that increases in the regenerating legs of Periplaneta americana (American cockroach). Insect Biochem. Mol. Biol. 28, 785–790. doi: 10.1016/S0965-1748(98)00058-7
Kuraku, S., Zmasek, C. M., Nishimura, O., and Katoh, K. (2013). aLeaves facilitates on-demand exploration of metazoan gene family trees on MAFFT sequence alignment server with enhanced interactivity. Nucleic Acids Res. 41, 22–28. doi: 10.1093/nar/gkt389
Lartigue, A., Campanacci, V., Roussel, A., Larsson, A. M., Jones, T. A., Tegoni, M., et al. (2002). X-ray structure and ligand binding study of a moth chemosensory protein. J. Biol. Chem. 277, 32094–32098. doi: 10.1074/jbc.M204371200
Laughlin, J. D., Ha, T. S., Jones, D. N. M., and Smith, D. P. (2008). Activation of pheromone-sensitive neurons is mediated by conformational activation of pheromone binding protein. Cell 133, 1255–1265. doi: 10.1016/j.cell.2008.04.046
Leal, W. S. (2013). Odorant reception in insects: roles of receptors, binding proteins, and degrading enzymes. Annu. Rev. Entomol. 58, 373–391. doi: 10.1146/annurev-ento-120811-153635
Li, G. W., Chen, X. L., Li, B. L., Zhang, G. H., Li, Y. P., and Wu, J. X. (2016a). Binding properties of general odorant binding proteins from the oriental fruit moth. Grapholita molesta (Busck) (Lepidoptera: Tortricidae). PLoS One 11:e0155096. doi: 10.1371/journal.pone.0155096
Li, G. W., Zhang, Y., Li, Y. P., Wu, J. X., and Xu, X. L. (2016b). Cloning, expression, and functional of three odorant binding proteins of the oriental fruit moth, Grapholita molesta (Busck) (Lepidoptera: Tortricidae). Arch. Inscet Biochem. Physiol. 91, 67–87. doi: 10.1002/arch.21309
Li, G. W., Chen, X. L., Xu, X. L., and Wu, J. X. (2018). Degradation of sex pheromone and plant volatile components by an antennal glutathione S-transferase in the oriental fruit moth, Grapholita molesta Busck (Lepidoptera: Tortricidae). Arch. Insect Biochem. Physiol. 99:e21512. doi: 10.1002/arch.21512
Li, G. W., Du, J., Li, Y. P., and Wu, J. X. (2015). Identification of putative olfactory genes from the oriental fruit moth Grapholita molesta via anantennal transcriptome analysis. PLoS One 10:e0142193. doi: 10.1371/journal.pone.0142193
Li, Z. Q., Zhang, S., Luo, J. Y., Zhu, J., Cui, J. J., and Dong, S. L. (2015). Expression analysis and binding assays in the chemosensory protein gene family indicate multiple roles in Helicoverpa armigera. J. Chem. Ecol. 41, 473–485. doi: 10.1007/s10886-015-0574-x
Liu, Q. J., Wang, H., Li, H. L., Zhang, J., Zhuang, S. L., Zhang, F. N., et al. (2013). Impedance sensing and molecular modeling of an olfactory biosensor based on chemosensory proteins of honeybee. Biosens. Bioelectron. 40, 174–179. doi: 10.1016/j.bios.2012.07.011
Liu, X. L., Luo, Q. A., Zhong, G. H., Rizwan-ul-Haq, M., and Hu, M. Y. (2010). Molecular characterization and expression pattern of four chemosensory proteins from diamondback moth, Plutella xylostella (Lepidoptera: Plutellidae). J. Biochem. 148, 189–200. doi: 10.1093/jb/mvq050
Liu, Y., Gu, S. H., Zhang, Y. J., Guo, Y. Y., and Wang, G. R. (2012). Candidate olfaction genes identified within the Helicoverpa armigera antennal transcriptome. PLoS One 7:e48260. doi: 10.1371/journal.pone.0048260
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262
Lu, P. F., Huang, L. Q., and Wang, C. Z. (2012). Identification and field evaluation of pear fruit volatiles attractive to the oriental fruit moth, Cydia molesta. J. Chem. Ecol. 38, 1003–1016. doi: 10.1007/s10886-012-0152-4
Lu, P. F., Qiao, H. L., Xu, Z. C., Cheng, J., Zong, S. X., and Luo, Y. Q. (2014). Comparative analysis of peach and pear fruit volatiles attractive to the oriental fruit moth, Cydia molesta. J. Plant Interact. 9, 388–395. doi: 10.1080/17429145.2013.843724
Lu, P. F., Wang, R., Wang, C. Z., Luo, Y. Q., and Qiao, H. L. (2015). Sexual differences in electrophysiological and behavioral responses of Cydia molesta to peach and pear volatiles. Entomol. Exp. Appl. 157, 279–290. doi: 10.1111/eea.12362
Maleszka, J., Forêt, S., Saint, R., and Maleszka, R. (2007). RNAi-induced phenotypes suggest a novel role for a chemosensory protein CSP5 in the development of embryonic integument in the honeybee (Apis mellifera). Dev. Genes Evol. 217, 189–196. doi: 10.1007/s00427-006-0127-y
Maleszka, R., and Stange, G. (1997). Molecular cloning, by a novel approach, of a cDNA encoding a putative olfactory protein in the labial palps of the moth Cactoblastis cactorum. Gene 202, 39–43. doi: 10.1016/S0378-1119(97)00448-4
Mohanty, S., Zubkov, S., and Gronenborn, A. M. (2004). The solution NMR structure of Antheraea polyphemus PBP provides new insight into pheromone recognition by pheromone-binding proteins. J. Mol. Biol. 337, 443–452. doi: 10.1016/j.jmb.2004.01.009
Myers, C. T., Hull, L. A., and Krawczyk, G. (2007). Effects of orchard host plants (apple and peach) on development of oriental fruit moth (Lepidoptera: Tortricidae). J. Econ. Entomol. 100, 421–430. doi: 10.1093/jee/100.2.421
Nagnan-Le, M. P., Cain, A. H., Jacquin-Joly, E., Francois, M. C., Ramachandran, S., Maida, R., et al. (2000). Chemosensory proteins from the proboscis of Mamestra brassicae. Chem. Senses 25, 541–553. doi: 10.1093/chemse/25.5.541
Najar-Rodriguez, A., Bellutti, N., and Dorn, S. (2013a). Larval performance of the oriental fruit moth across fruits from primary and secondary hosts. Physiol. Entomol. 38, 63–70. doi: 10.1111/phen.12003
Najar-Rodriguez, A., Orschel, B., and Dorn, S. (2013b). Season-long volatile emissions from peach and pear trees in situ, overlapping profiles, and olfactory attraction of an oligophagous fruit moth in the laboratory. J. Chem. Ecol. 39, 418–429. doi: 10.1007/s10886-013-0262-7
Natale, D., Mattiacci, L., Hern, A., Pasqualini, E., and Dorn, S. (2003). Response of female Cydia molesta (Lepidoptera: Tortricidae) to plant derived volatiles. Bull. Entomol. Res. 93, 335–342. doi: 10.1079/BER2003250
Natale, D., Mattiacci, L., Pasqualini, E., and Dorn, S. (2004). Apple and peach fruit volatiles and the apple constituent butyl hexanoate attract female oriental fruit moth, Cydia molesta, in the laboratory. J. Appl. Entomol. 128, 22–27. doi: 10.1046/j.1439-0418.2003.00802.x
Peng, Y., Wang, S. N., Li, K. M., Liu, J. T., Zheng, Y., Shan, S., et al. (2017). Identification of odorant binding proteins and chemosensory proteins in Microplitis mediator as well as functional characterization of chemosensory protein 3. PLoS One 12:e0180775. doi: 10.1371/journal.pone.0180775
Petersen, T. N., Brunak, S., von Heijne, G., and Nielsen, H. (2011). SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat. Methods 8, 785–786. doi: 10.1038/nmeth.1701
Piñero, J. C., and Dorn, S. (2007). Synergism between aromatic compounds and green leaf volatiles derived from the host plant underlies female attraction in the oriental fruit moth. Entomol. Exp. Appl. 125, 185–194. doi: 10.1111/j.1570-7458.2007.00614.x
Piñero, J. C., and Dorn, S. (2009). Response of female oriental fruit moth to volatiles from apple and peach trees at three phonological stages. Entomol. Exp. Appl. 131, 67–74. doi: 10.1111/j.1570-7458.2009.00832.x
Piñero, J. C., Giovanni, G. C., and Dorn, S. (2008). Synergistic behavioral responses of female oriental fruit moths (Lepidoptera:Tortricidae) to synthetic host plant-derived mixtures are mirrored by odor-evoked calcium activity in their antennal lobes. J. Insect Physiol. 54, 333–343. doi: 10.1016/j.jinsphys.2007.10.002
Rajapakse, C. N. K., Walter, G. H., Moore, C. J., Hull, C. D., and Cribb, B. W. (2006). Host recognition by a polyphagous lepidopteran (Helicoverpa armigera): primary host plants, host produced volatiles and neurosensory stimulation. Physiol. Entomol. 31, 270–277. doi: 10.1111/j.1365-3032.2006.00517.x
Sabatier, L., Jouanguy, E., Dostert, C., Zachary, D., Dimarcq, J. L., Bulet, P., et al. (2003). Pherokine-2 and -3 : two Drosophila molecules related to pheromone/odor-binding proteins induced by viral and bacterial infection. Eur. J. Biochem. 270, 3398–3407. doi: 10.1046/j.1432-1033.2003.03725.x
Sandler, B. H., Nikonova, L., Leal, W. S., and Clardy, J. (2000). Sexual attraction in the silkworm moth: structure of the pheromone-binding-protein-bombykol complex. Chem. Biol. 7, 143–151. doi: 10.1016/S1074-5521(00)00078-8
Sato, K., and Touhara, K. (2009). Insect olfaction: receptors, signal transduction, and behavior. Results Probl. Cell Differ. 47, 121–138. doi: 10.1007/400_2008_10
Sirot, L. K., Poulson, R. L., Mckenna, M. C., Girnary, H., Wolfner, M. F., and Harrington, L. C. (2008). Identity and transfer of male reproductive gland proteins of the dengue vector mosquito, Aedes aegypti: potential tools for control of female feeding and reproduction. Insect Biochem. Mol. Biol. 38, 176–189. doi: 10.1016/j.ibmb.2007.10.007
Steiner, L. F., and Yetter, W. P. (1933). Second report on the efficiency of bait traps for the oriental fruit moth as indicated by the release and capture of marked adults. J. Econ. Entomol. 26, 774–788. doi: 10.1093/jee/26.4.774
Tan, J., Song, X. M., Fu, X. B., Wu, F., Hu, F. L., and Li, H. L. (2018). Combinatorial multispectral, thermodynamics, docking and site directed mutagenesis reveal the cognitive characteristics of honey bee chemosensory protein to plant semiochemical. Spectrochim. Acta A Mol. Biomol. Spectrosc. 201, 346–353. doi: 10.1016/j.saa.2018.04.074
Tomaselli, S., Crescenzi, O., Sanfelice, D., Eiso, A. B., Wechselberger, R., Angeli, S., et al. (2006). Solution structure of a chemosensory protein from the desert locust Schistocerca gregaria. Biochemistry 45, 10606–10613. doi: 10.1021/bi060998w
Wanner, K. W., Willis, L. G., Theilmann, D. A., Isman, M. B., Feng, Q., and Plettner, E. (2004). Analysis of the insect OS-D-like gene family. J. Chem. Ecol. 30, 889–911. doi: 10.1023/B:JOEC.0000028457.51147.d4
Waris, M. I., Younas, A., Ul Qamar, M. T., Hao, L., Ameen, A., Ali, S., et al. (2018). Silencing of chemosensory protein gene NlugCSP8 by RNAi induces declining behavioral responses of Nilaparvata lugens. Front. Physiol. 9:379. doi: 10.3389/fphys.2018.00379
Wilkins, M. R., Gasteiger, E., Bairoch, A., Sanchez, J. C., and Williams, K. L. (1999). Protein identification and analysis tools in the ExPASy server. Methods Mol. Biol. 112, 531–532. doi: 10.1385/1-59259-890-0:571
Xu, X. Z., Xu, W., Rayo, J., Ishida, Y., Leal, W. S., and Ames, J. B. (2010). NMR structure of navel orangeworm moth pheromone-binding protein (AtraPBP1): implications for pH-sensitive pheromone detection. Biochemistry 49, 1469–1476. doi: 10.1021/bi9020132
Yang, K., He, P., and Dong, S. L. (2014). Different expression profiles suggest functional differentiation among chemosensory proteins in Nilaparvata lugens (Hemiptera: Delphacidae). J. Insect Sci. 14:270. doi: 10.1093/jisesa/ieu132
Ye, J., Coulouris, G., Zaretskaya, I., Cutcutache, I., Rozen, S., and Madden, T. (2012). Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. BMC Bioinformatics 13:134. doi: 10.1186/1471-2105-13-134
Younas, A., Waris, M. I., Tahir Ul Qamar, M., Shaaban, M., Prager, S. M., and Wang, M. Q. (2018). Functional analysis of the chemosensory protein MsepCSP8 from the oriental armyworm Mythimna separata. Front. Physiol. 9:872. doi: 10.3389/fphys.2018.00872
Zhang, Y. N., Jin, J. Y., Jin, R., Xia, Y. H., Zhou, J. J., Deng, J. Y., et al. (2013). Differential expression patterns in chemosensory and non-chemosensory tissues of putative chemosensory genes identified by transcriptome analysis of insect pest the purple stem borer Sesamia inferens (Walker). PLoS One 8:e69715. doi: 10.1371/journal.pone.0069715
Zhang, Y. N., Ye, Z. F., Yang, K., and Dong, S. L. (2014). Antenna-predominant and male-biased CSP19 of Sesamia inferens is able to bind the female sex pheromones and host plant volatiles. Gene 536, 279–286. doi: 10.1016/j.gene.2013.12.011
Zhong, T., Yin, J., Deng, S. S., Li, K. B., and Cao, Y. Z. (2012). Fluorescence competition assay for the assessment of green leaf volatiles and trans-beta-farnesene bound to three odorant-binding proteins in the wheat aphid Sitobion avenae (Fabricius). J. Insect Physiol. 58, 771–781. doi: 10.1016/j.jinsphys.2012.01.011
Zhou, J. J., Robertson, G., He, X., Dufour, S., Hooper, A. M., Pickett, J. A., et al. (2009). Characterization of Bombyx mori odorant-binding proteins reveals that a general odorant-binding protein discriminates between sex pheromone components. J. Mol. Biol. 389, 529–545. doi: 10.1016/j.jmb.2009.04.015
Keywords: Grapholita molesta, chemosensory protein, chemoreception, fluorescence ligand-binding assays, molecular docking, site-directed mutagenesis
Citation: Li G-W, Chen X-L, Chen L-H, Wang W-Q and Wu J-X (2019) Functional Analysis of the Chemosensory Protein GmolCSP8 From the Oriental Fruit Moth, Grapholita molesta (Busck) (Lepidoptera: Tortricidae). Front. Physiol. 10:552. doi: 10.3389/fphys.2019.00552
Received: 20 December 2018; Accepted: 18 April 2019;
Published: 07 May 2019.
Edited by:
Peng He, Guizhou University, ChinaReviewed by:
Guan-Heng Zhu, University of Kentucky, United StatesYa-Nan Zhang, Huaibei Normal University, China
Joe Hull, United States Department of Agriculture, United States
Copyright © 2019 Li, Chen, Chen, Wang and Wu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wen-Qiang Wang, wenqiangwang@yau.edu.cn Jun-Xiang Wu, Junxw@nwsuaf.edu.cn
†These authors have contributed equally to this work