Figures
Abstract
Odorant-binding proteins (OBPs) and chemosensory proteins (CSPs) of aphids are thought to be responsible for the initial molecular interactions during olfaction that mediate detection of chemical signals. Analysis of the diversity of proteins involved comprises critical basic research work that will facilitate the development of sustainable pest control strategies. To help us better understand differences in the olfactory system between winged and wingless grain aphids, we constructed an antennal transcriptome from winged and wingless Sitobion avenae (Fabricius), one of the most serious pests of cereal fields worldwide. Among the 133,331 unigenes in the antennal assembly, 13 OBP and 5 CSP putative transcripts were identified with 6 OBP and 3 CSP sequences representing new S. avenae annotations. We used qPCR to examine the expression profile of these genes sets across S. avenae development and in various tissues. We found 7 SaveOBPs and 1 SaveCSP were specifically or significantly elevated in antennae compared with other tissues, and that some transcripts (SaveOBP8, SaveCSP2 and SaveCSP5) were abundantly expressed in the legs of winged or wingless aphids. The expression levels of the SaveOBPs and SaveCSPs varied depending on the developmental stage. Possible physiological functions of these genes are discussed. Further molecular and functional studies of these olfactory related genes will explore their potential as novel targets for controlling S. avenae.
Citation: Xue W, Fan J, Zhang Y, Xu Q, Han Z, Sun J, et al. (2016) Identification and Expression Analysis of Candidate Odorant-Binding Protein and Chemosensory Protein Genes by Antennal Transcriptome of Sitobion avenae. PLoS ONE 11(8): e0161839. https://doi.org/10.1371/journal.pone.0161839
Editor: J Joe Hull, USDA Agricultural Research Service, UNITED STATES
Received: January 26, 2016; Accepted: August 12, 2016; Published: August 25, 2016
Copyright: © 2016 Xue et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All the sequences files of odorant binding proteins and chemosensory proteins identified from antennal transcriptomes of S. avenae are available from the Genbank database (accession numbers KU140605-22).
Funding: This work was supported by: 1. National Natural Science Foundation of China (http://www.nsfc.gov.cn/) 31371946. JLCH conceived and designed the experiments and made decision to publish, 2. National Natural Science Foundation of China (http://www.nsfc.gov.cn/) 31401740. JF conceived and designed the experiments and contributed to the writing of the manuscript, 3. The Agricultural Science and Technology Innovation Program (ASTIP). JLCH conceived and designed the experiments and made decision to publish, 4. The National Science and Technology Major Project (http://www.nmp.gov.cn/) 2014ZX0800201B. JF conceived and designed the experiments and contributed to the writing of the manuscript, and 5. Cooperation Project between China and Belgium from Ministry of Science and Technology (MOST) (http://www.istcp.org.cn/) 2014DFG32270. JLCH conceived and designed the experiments and decision to publish.
Competing interests: The authors have declared that no competing interests exist.
Introduction
The grain aphid, Sitobion avenae (Fabricius) is one of the most important pests of gramineous crops [1]. Damage from their sucking plant sap slows plant growth and reduces the number of tillers, severely diminishing the yield and quality of wheat [2]. In recent years, global climate warming, farming system changes and other factors have contributed to significant enhancement in the reproductive capacity and adaptability of this pest, resulting in great crop damage [3]. Therefore, research on environmentally safe prevention and control strategies for widespread use is extremely important [4, 5]. Detailed analysis of the S. avenae chemical sensing system can provide insights into aphid olfactory physiology and the molecular mechanisms used to detect semiochemicals. These insights into the chemoreception mechanism of aphids could aid in developing novel olfactory-based control strategies, such as repellents or attractants developing, and provide additional candidate genes for targeted disruption, which could interfere with both plant-aphid interactions and aphid responses of aphid to the external environment.
For aphids, olfaction plays an important role in distinguishing host plant volatiles from other environmental volatiles [6]. Antennae, which are largely the physiological basis for insect chemical ecology, are one of the principle organs that aphids use to recognize chemical information in the environment. Many olfactory-related proteins are responsible for discerning chemical information and regulating and controlling aphid behaviors such as selecting hosts and avoiding natural enemies [7–11]. A variety of odor-related proteins, such as odorant-binding proteins (OBPs), chemosensory proteins (CSPs), odorant receptors (ORs), sensory neuron membrane proteins (SNMPs) and odorant degrading enzymes (ODEs), contribute to the initiation of olfactory perception in insects. OBPs and CSPs, which are water-soluble, globular proteins that are concentrated (as high as 10 mM) in the sensillum lymph of insect antennae [12–17], are thought to provide the initial molecular interactions for chemical signals (semiochemicals) and to ferry the semiochemical molecules through the antennal sensillum lymph to the olfactory receptors (ORs).
Since the first insect OBP was identified from the antennae of male Antheraea polyphemus in 1981 [12], many new OBPs have been identified, the discovery of which has been greatly aided by new generation genome and transcriptome analyses. Genome sequencing facilitated the identification of 15 OBP genes in the hemipteran Acyrthosiphon pisum [18]. In contrast, mRNA-based transcriptome sequencing has led to the discovery of numerous OBPs from diverse hemipteran pests including Lygus lineolaris [19], Adelphocoris lineolatus [20], Aphis gossypii [21], Nilaparvata lugens [22, 23], Sogatella furcifera [24], and Apolygus lucorum [25]. Interestingly, the number of OBPs reported has differed depending on the species (eg. A. gossypii has 9 OBPs, L. lineolaris has 33 OBPs). Based on qPCR analyses, it has become clear that many of the identified OBP genes are highly expressed in antennae [18–25]. However, OBPs are also highly expressed in other tissues, such as legs and heads, which suggests that these OBPs might be associated with taste perception or participate in other physiological functions [26, 27].
OBPs can vary in the amino acid spacing between helices, in the lengths and positions of the loops that connect the helices and in the lengths of the C- and N-termini [28]. The nature and shape of the binding cavities can also differ, presumably representing specificities for different odorant molecules. For assessing the role of OBPs in olfaction in vitro, fluorescence competitive binding assays have been used to examine the binding capacity of recombinantly expressed OBPs for specific odorant ligands [15]. Although aphids are a major hemipteran group with numerous species throughout the world, a limited number of ligands have been assayed. The binding properties of three classical pea aphid OBPs (ApisOBP1, ApisOBP3 and ApisOBP8) have been investigated using 12 chemical compounds [29]. Research has confirmed that OBPs are keys in the perception of (E)-β-farnesene (EβF), an aphid alarm pheromone widely used in olfaction-based approaches to control aphid populations [30]. The chemical can interfere with an insect’s ability to find suitable hosts or olfactory cues regarding natural enemies, and has been used successfully in the “push–pull” pest control strategy [31, 32]. Different OBPs in various aphid species have been reported to play critical roles in discerning EβF in each respective species. Ligand binding assays suggested that ApisOBP3 and SaveOBP7 each have high binding affinities for EβF [29, 33].
Like OBPs, CSPs are also secreted proteins but they have a lower molecular weight (10–16 kDa). The CSP family contributes to the recognition of sex pheromones [34] and general odorants [35] and to other functions such as leg regeneration [36], development [37] and feeding [38]. Transcriptome sequence data for hemipterans are currently limited to 13 CSPs from A. pisum [18], 11 CSPs from N. lugens [39] and 9 CSPs from A. gossypii [21], and their exact functions are largely unknown. Examination of gene expression profiles, particularly in various parts of the body, and phylogenetic analyses could potentially provide important information concerning the function of CSP genes. Two CSPs identified in Megoura viciae, MvicOS-D1 and MvicOS-D2 did not bind any of 28 compounds known to elicit an electrophysiological response in electroantennograms or in single olfactory neuron preparations [40].
Our understanding of the molecular components comprising the S. avenae olfactory system is incomplete both in sequences and expression data. To understand the physiological mechanism underlying volatile detection in wheat aphids better, more research on olfactory related proteins has been needed. In our present work, we used an antennal transcriptome to identify OBPs and CSPs in winged and wingless S. avenae, and used quantitative real-time PCR to assess their expression in different tissues and ages. In addition, we discuss potential roles of the identified OBPs/CSPs transcripts in olfactory or other physiological processes.
Materials and Methods
Ethics statement
Grain aphids (S. avenae) were from a parthenogenetic colony initially collected from wheat at Langfang Farm, (the experimental farm of the Institute of Plant Protection, Chinese Academy of Agricultural Sciences), Hebei Province, China, which is not privately owned or protected. The grain aphid is neither endangered nor protected so no specific permission was required for its collection.
Aphid samples
S. avenae adults were collected from wheat fields at Langfang (39°30′42′′N, 116°36′7′′E), Hebei Province, China, in 2015, and a single wingless aphid was chosen to be reared as a parthenogenetic colony in the laboratory at 22 ± 1°C, with 75% relative humidity and 16h light/8 h dark. After 10 generations, antennae from 1000 winged and 1000 wingless adult aphids of S. avenae were collected, immediately frozen in liquid nitrogen, and stored at −80°C for transcriptome sequencing. The winged and wingless antennae for transcriptome analysis were collected five times, with ~200 aphids from each sample collected per time. The antennae were pulled off with forceps. Tissues from winged and wingless adults (1000 winged antennae, 1000 wingless antennae, 20 winged heads without antennae, 20 wingless heads without antennae, 20 winged thoraxes, 20 wingless thoraxes, 20 winged abdomens, 20 wingless abdomens, 1000 winged legs and 1000 wingless legs) and whole bodies from each development stage of aphid (15 individuals of 1st instar nymphs or 2nd instar nymphs, 10 of winged or wingless 3rd instar nymphs, 5 of winged or wingless 4th instar nymphs, 5 of winged adults, respectively) were collected in tubes as respective bulked samples with three replicates and immediately frozen in liquid nitrogen and stored at −80°C until needed.
Transcriptome sequencing, assembly and functional annotation
Antennal total RNA from winged or wingless antennae of S. avenae was extracted using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions. After assessing RNA quality on a spectrophotometer NanoDrop 2000c (Thermo Fisher Scientific, Waltham, MA, USA), a 3 μg RNA sample with standard quality ratios (1.8 < OD260/280 < 2.1) was purified using poly-T oligo-attached magnetic beads. Divalent cations under elevated temperature in NEBNext First Strand Synthesis Reaction Buffer (5×) was used for fragmentation. Single-stranded (ss) cDNA was synthesized using a random hexamer primer, M-MuLV Reverse Transcriptase and DNA Polymerase I and RNase H (NEB, USA). The 3' ends of the DNA fragments were adenylated and the NEBNext Adaptor was ligated to the fragments for hybridization. The library fragments were purified with the AMPure XP system (Beckman Coulter, Beverly, MA, USA) to size select cDNA fragments ~ 150–200 bp length. Then 3 μL USER Enzyme (NEB, USA) was used with size-selected, adaptor-ligated cDNA at 37°C for 15 min followed by 5 min at 95°C prior to PCR. PCR was performed with Phusion High-Fidelity DNA polymerase, Universal PCR primers and Index (X) Primer. The products were purified (AMPure XP system), and library quality was assessed using the Agilent Bioanalyzer 2100 system (Agilent Technologies, CA, USA). Clustering of the index-coded samples was performed on a cBot Cluster Generation System using TruSeq PE Cluster Kit v3-cBot-HS (Illumina, China) according to the manufacturer’s instructions. The library preparations were sequenced on an Illumina Hiseq 2500 platform and paired-end reads (the sequencing strategy was PE125) were generated after cluster generation. After sequencing, the raw reads were processed to remove low quality and adaptor sequences by ng_qc, and then assembled into unigenes using Trinity r20140413p1 min_kmer_cov:2 and other default parameters [41]. Then the unigenes were annotated using seven databases, including the non-redundant protein sequence (Nr, e-value = 1e-5), non-redundant nucleotide (Nt, e-value = 1e-5), Pfam (e-value = 0.01), Clusters of Orthologous Groups (KOG/COG, e-value = 1e-3), Swiss-Prot (e-value = 1e-5), Kyoto Encyclopedia of Genes and Genomes (KEGG, e-value = 1e-10) and Gene Ontology (GO, e-value = 1e-6) databases.
Identification and verification of transcripts encoding putative OBPs and CSPs
We used a motif search program consisting of C1-X15-39-C2-X3-C3-X21-44-C4-X7-12-C5-X8-C6 for OBPs [42] and C1-X6-8-C2-X16-21-C3-X2-C4 for CSPs [43] and the BLASTx program at the National Center for Biotechnology Information (NCBI, http://blast.ncbi.nlm.nih.gov/Blast.cgi) to confirm putative OBP and CSP genes. The candidate OBPs and CSPs were cloned and sequence validated. TransScript First-Strand cDNA Synthesis SuperMix (Transgen, Beijing, China) was used to synthesize template cDNA. PCR reactions were carried out with 2× Taq DNA polymerase (Transgen, Beijing, China) with an initial denaturation at 95°C for 5min; 35 cycles of 94°C for 45 s, 55°C for 30 s, 72°C for 1 min; and a final extension at 72°C for 10 min. The PCR products were gel-purified and subcloned into the pEASY-T1 Sample Cloning Vector (Transgen, Beijing, China) and sequenced using standard M13 primers. Gene-specific primers to clone ORF sequences of each OBP and CSP gene were designed using the program Primer3 (http://primer3.ut.ee/) (S1 Table).
Comparative analysis of transcripts for putative OBPs and CSPs
For comparing the differential expression of putative genes in the winged and wingless antennal transcriptomes of S. avenae, the read number for the OBPs and CSPs between different morph antennae of OBPs and CSPs was converted to RPKM (reads per kilobase per million mapped reads), using the formula: RPKM (A) = (1,000,000 × C × 1,000) / (N × L), where RPKM (A) is the expression of gene A, C is the number of reads uniquely aligned to gene A, N is the total number of reads uniquely aligned to all unigenes, and L is the number of bases in gene A. The RPKM method eliminates the influence of gene length and sequencing depth on the calculation of gene expression [44].
Sequence analysis and phylogenetic tree construction
The putative N-terminal signal peptides and most likely cleavage site were predicted by the SignalP 4.1 Server (http://www.cbs.dtu.dk/services/SignalP/). The amino acid sequence was deduced by the WebLab showorf program (http://weblab.cbi.pku.edu.cn/) and the candidate OBPs and CSPs were aligned using CLUSTAL Omega with default parameters (http://www.ebi.ac.uk/Tools/msa/clustalo/) and then arranged by BOXSHADE 3.21 (http://www.ch.embnet.org/software/BOX_form.html). The candidate OBPs and CSPs from S. avenae and related sequences from other aphid species were chosen for phylogenetic analysis. After sequences were aligned using ClustalX 2.1 with default gap penalty parameters of gap opening 10 and extension 0.2, all the phylogenetic trees were constructed using the neighbour-joining method implemented in MEGA 5.0 [45] with default settings and 1000 bootstrap replications.
A total of 167 OBP protein sequences from 21 hemipteran species were used for the phylogenetic analysis, sequences used are described in S2 Table and include 13 OBPs from S. avenae identified in the present study, 1 from Aphis craccivora, 2 from Aphis fabae, 10 from Aphis glycines, 13 from A. pisum, 9 from A. gossypi, 1 from Brevicoryne brassicae, 1 from Drepanosiphum platanoidis, 7 from Metopolophium dirhodum, 5 from Megoura viciae, 6 from Myzus persicae, 5 from Nasonovia ribisnigri, 5 from Pterocomma salicis, 2 from Lipaphis erysimi, 5 from Rhopalosiphum padi, 1 from Tuberolachnus salignus, 12 from A. lucorum, 14 from A. lineolatus, 33 from L. lineolaris, 10 from N. lugens and 12 from S. furcifera. In addition, 51 CSPs from seven hemipteran species were used for the phylogenetic analysis and include 5 CSPs from S. avenae identified in the present study, 7 from A. gossypi, 3 from M. persicae, 8 from A. lucorum, 8 from A. lineolatus, 9 from S. furcifera and 11 from N. lugens. The SaveOBP and SaveCSP accessions are listed in Table 1. The accession numbers for the other genes are in S2 Table.
Expression analysis of OBPs and CSPs in different-aged aphids and tissues
Total RNA from different tissues of winged and wingless aphids (antennae, heads without antennae, thoraxes, abdomens and legs) and whole bodies from each stage were extracted using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions. The ratio of OD260/280 was measured on a spectrophotometer NanoDrop 1000 (Thermo Fisher Scientific, Pittsburgh, PA, USA). Single-stranded cDNA templates using 1 μg RNA from various samples were synthetized using TransScript One-Step gDNA Removal and cDNA Synthesis SuperMix (Transgen, Beijing, China) according to the manufacturer’s instructions.
The expression level of each SaveOBP and SaveCSP transcript across the developmental stages stages and tissues in winged and wingless aphids was assessed using quantitative real-time PCR (qPCR). Specific primer pairs for qPCR were designed with Primer 3 (S3 Table), and qPCR was performed on an ABI 7500 Real-Time PCR System (Applied Biosystems, Carlsbad, CA). Two reference genes, β-actin and NADH dehydrogenase were used for normalizing target gene expression and correcting for sample-to-sample variation [46]. The qPCR reactions were performed in 20 μL reactions containing 10 μL SYBR Premix Ex Taq (TaKaRa, Beijing, China), 0.5 μL of each primer (10 μM), 0.4 μL Rox Reference Dye, 2 μL sample cDNA, and 6.6 μL sterilized H2O. The qPCR cycling parameters were 95°C for 30 s, followed by 40 cycles of 95°C for 15 s and 60°C for 30 s. Next the PCR products were heated to 95°C for 15 s, cooled to 60°C for 1 min and 95°C for 15 s to measure the dissociation curves. Negative controls without a template were included in each experiment to check reproducibility; each qPCR reaction for each sample was done in three technical replicates and three biological replicates for each transcript. Standard curves for reference genes and candidate genes were generated by gradient dilution to identify proper primers with 90–110% amplification efficiency and without nonspecific amplification. Relative quantities were calculated using the Vandesompele Method [47]. Differences in transcript expression in various tissues and ages were statistically analyzed with a one-way ANOVA using SAS 9.1 (SAS Institute, Cary, NC, USA) followed by the least-significant difference (LSD) method.
Results
Overview of transcriptomes
A total of 2.22 and 2.30 million raw reads were obtained from S. avenae antennae libraries from winged and wingless aphids, respectively. After removal of low-quality, adaptor, and contaminating sequences, 2.14 and 2.23 million clean reads were retained and assembled into 147,665 distinct transcripts (mean length = 652 bp) and 133,331 unigenes (mean length = 594 bp). The length distribution can be seen in S1 Fig.
In total, 89,452 (67.06% of all 133,331 unigenes), 42,393 (31.79%), 54,572(40.92%), 39,254 (29.44%), 55,639 (41.72%), 32,300 (24.22%) and 60,371 (45.27%) transcripts from S. avenae were annotated using the Nr, Nt, Pfam, KOG/COG, Swiss-Prot, KEGG and GO databases respectively (S4 Table).
In the GO annotation, biological process, cellular process and metabolic process were the most abundant GO terms. The cluster for biological process was the next largest group. Most transcripts that corresponded to molecular function were related to binding and catalytic activity (S2 Fig). In the KOG classification, unigenes clustered into 26 categories (S3 Fig). Among these categories, general function prediction was the dominant category, followed by signal transduction and post-translational modification, protein turnover and chaperon. All the unigenes annotated in the KO database were assigned to the 5 biological pathways described in the KEGG database: cellular processes, environmental information processing, genetic information processing, metabolism, and organismal systems (S4 Fig). The most common pathway was metabolism followed by genetic information processing, organismal systems and cellular processes. In the environmental information processing group, most genes (2776) were involved in signal transduction.
Identification and analysis of OBP genes in S. avenae
Thirteen putative OBPs were identified (Table 1) using a motif search and the NCBI BLASTx program. We named the OBP genes SaveOBP1 to 10 and SaveOBP13 to 15, following the nomenclature established for A. pisum [18, 21]. All OBP transcripts were confirmed by molecular cloning and sequencing. The sequencing results showed no differences with the trancriptomic data. Apart from SaveOBP13, the other 12 OBPs had complete open reading frames (ORFs) consisting of 400–750 bp nucleotides. Among the 13 OBPs, all have the characteristic insect OBP sequence motif [25] (Fig 1A). 9 SaveOBPs (SaveOBP2–3, 7–10 and 13–15) had the classic hemipteran OBP Cys motif (C1-X22-32-C2-X3-C3-X36-46-C4-X8-14-C5-X8-C6) [48]. SaveOBP4 had 49 amino acids between the first and second conserved cysteines and 21 amino acids between the fourth and fifth conserved cysteines. SaveOBP1 had 47 amino acids between the third and fourth conserved cysteines. 2 SaveOBPs (SaveOBP5 and SaveOBP6) belong to the ‘Plus-C’ OBP family and have the similar motif to the Cys spacing pattern C1-X20-41-C2-X3-C3-X41-46-C4-X19-29-C4a-X9-C5-X8-C6-P-X9-10-C6a-X9-10 [49] (Fig 1A). The signal peptide predictions are shown in Table 1. The 13 OBPs shared 6.21–30.24% amino acid identities with each other (S5 Table).
Sequences were aligned by Clustal Omega and edited using BoxShade. Black boxes show conserved cysteines. The conserved Cys residues are indicated. Shading represents sequence identity > 70%.
The phylogenetic tree for the 165 OBPs from bugs, aphids and planthoppers, revealed diverging relationships. Aphid OBPs clustered into 10 major groups, each containing several homologous OBPs from different aphid species [26]; amino acid identities for each group (OBP1–10) were 79.17%, 90.09%, 86.78%, 88.45%, 84.10%, 78.68%, 73.84%, 89.30%, 86.75% and 91.94%, respectively, indicating high conservation among the different aphid species. Most of the orthologous sequences in the S. avenae OBPs were largely limited to aphid OBPs with an average bootstrap value of 94%, and 4 had a high degree of similarity with A. pisum (Fig 2), which suggested that these sequences might have vertically descended from the same ancestors and were conserved for common functions in aphids. SaveOBP14 and SaveOBP15 were clustered with OBPs of bugs and planthoppers.
The tree was constructed using MEGA 5.0 with bootstrap support based on 1000 iterations. Aphid sequences are in red. Bug sequences are in blue. Planthopper sequences are in gray. Major clades for aphid OBPs are marked in a different color. Abbreviation for these aphid species: Save, S. avenae; Acra, A. craccivora; Psal, P. salicis; Brebr, B. brassicae; Apis, A. pisum; Dpla, D. platanoidis; Mper, M. persicae; Nrib, N. ribisnigri; Rpad, R. padi; Mvic, M. viciae; Tsal, T. salignus; Afab, A. fabae; Agos, A. gossypii; Mdir, M. dirhodum; Lery, L. erysimi; Agly, A. glycines; Aluc, A. lucorum; Alin, A. lineolatus; Llin, L. lineolaris; Nlug, N. lugens; Sfur, S. furcifera.
Identification of CSP genes in S. avenae
In total, we identified five transcripts belonging to the CSP family (Table 1); SaveCSP1–5. By BLASTx comparative analysis, all these genes were most similar to CSP sequences from aphids. All of the CSP transcripts were confirmed by molecular cloning and sequencing, showing that there was no difference between the transcriptome sequencing result and the sequence clones. The CSPs had structural features typical of insect CSPs with four conserved cysteines fitting the hemipteran Cys spacing motif C1-X5-6-C2-X18-19-C3-X2-C4 (Fig 1B) [48]. Only SaveCSP1 did not have a complete ORF, whereas the other four CSPs varied in length from 400–700 bp. The five CSPs shared 11.89–30.82% amino acid identities (S6 Table). The predicted results of the signal peptides are shown in Table 1.
To assign functions to each of the SaveCSPs, we constructed a phylogenetic tree using 51 identified CSPs from seven hemipteran species. SaveCSP1, SaveCSP3 and SaveCSP5 clustered with aphid CSPs with an amino acid identity of 65.67%, 95.68% and 59.51%, respectively (Fig 3). At present, there are few studies on aphid CSP genes in hemipteran species.
The tree was constructed using MEGA 5.0 with bootstrap support based on 1000 iterations. SaveOBP sequences are in bold. Abbreviation of these aphid spaces are as follows: Save, S. avenae; Mper, M. persicae; Agos, A. gossypii; Aluc, A. lucorum; Alin, A. lineolatus; Nlug, N. lugens; Sfur, S. furcifera.
Expression profile analysis of SaveOBPs
The results of qPCR showed that the expression levels of SaveOBP1, SaveOBP2, SaveOBP3, SaveOBP5 and SaveOBP7 genes were higher than those of other OBPs in all instars (Fig 4). Expression of SaveOBP8 and SaveOBP10 both increased when aphids became winged adults. At the same time, expression of SaveOBP6 increased when aphids became wingless adults (Fig 5).
Fold changes for 1st or 2nd instar nymphs are relative to transcript levels of SaveOBP15. Fold-changes for other stages are relative to transcript levels of wingless SaveOBP15 in wingless aphids of the same age. Differences in mean transcript levels were compared using one-way ANOVA, followed by the least-significant difference (LSD) method. Bars with different letters indicate significant differences (p < 0.05).
Fold-changes are relative to transcript levels in 1st instar nymphs. Differences in mean transcript levels were compared using one-way ANOVA, followed by the least-significant difference (LSD) method. Bars with different letters indicate significant differences (p < 0.05). 1 In: 1st instar nymph; 2 In: 2nd instar nymph; 3W In: 3rd winged instar nymph; 3WL In: 3rd wingless instar nymph; 4W In: 4th winged instar nymph; 4WL In: 4th wingless instar nymph; W: winged adult; WL: wingless adult.
Seven SaveOBP transcripts (SaveOBP6, SaveOBP7, SaveOBP9, SaveOBP10, SaveOBP13, SaveOBP14 and SaveOBP15) were specifically or highly expressed in antennae of both winged and wingless aphids. Based on RPKM, these seven genes were more abundant than the other SaveOBPs in winged and wingless antennae. Among 13 SaveOBP genes, SaveOBP10 had the highest RPKM value, followed by SaveOBP9, SaveOBP7, SaveOBP14, SaveOBP13, SaveOBP6 and SaveOBP15 (Table 1). Expression of SaveOBP10 in the antennae of winged aphids was higher than in those of wingless aphids (Fig 6).
Fold-changes are relative to transcript levels in abdomens of winged adult aphids. Differences in mean transcript levels were compared using one-way ANOVA, followed by the least-significant difference (LSD) method. Bars with different letters indicate significant differences (p < 0.05). An: antennae, H: heads, T: thoraxes, Ab: abdomens, L: legs.
We also found that expression of SaveOBP1 was significantly higher in the abdomens of both winged and wingless aphids than in other tissues. SaveOBP2 was primarily expressed in heads. SaveOBP3 was significantly elevated in wingless thorax/abdomen. The expression of SaveOBP8 in the legs of winged aphids was remarkably high (>40 fold higher than in other tissues) (Fig 6).
Expression profile analysis of SaveCSPs
In the developmental qPCR analysis of the five SaveCSPs across each stage (Fig 7), SaveCSP1 and SaveCSP2 had the highest expression level in 1st and 2nd instar nymphs. Expressions of SaveCSP2 and SaveCSP3 were higher in winged 3rd instar nymphs than in wingless 3rd instar nymphs. For 4th instar nymphs, the expression of SaveCSP2 and SaveCSP5 in winged nymphs was higher than in wingless nymphs. SaveCSP2 was significantly higher in winged and wingless adults. The only difference between winged and wingless adults was SaveCSP1, which had higher expression in the antennae of winged adult aphids (>200 fold higher than in the antenna of wingless aphids). SaveCSP2 was highest in legs of wingless adults, whereas SaveCSP5 was highest in the legs of winged aphids. SaveCSP4 was significantly higher in abdomens of winged and wingless adults (Fig 8).
Fold-changes for 1st and 2nd instar nymphs are relative to the transcript levels of SaveCSP1. Fold-changes for other stages are relative to transcript levels of SaveCSP1 in wingless aphids of the same stage. Differences in mean transcript levels were compared using one-way ANOVA, followed by the least-significant difference (LSD) method. Bars with different letters indicate significant differences (p < 0.05).
Fold-changes are relative to transcript levels in abdomens of winged adult aphids. Differences in mean transcript levels were compared using one-way ANOVA, followed by the least-significant difference (LSD) method. Bars with different letters indicate significant differences (p < 0.05). An: antennae, H: heads, T: thoraxes, Ab: abdomens, L: legs.
Discussion
In this study, we report on the sequencing, assembly and partial annotation of anntenal transcriptomes for winged and wingless grain aphids (S. avenae) and we identified 13 OBPs and five CSPs. Among them, six OBPs and three CSPs are newly reported transcripts. The number of OBPs is comparable to the 15 for A. pisum [18], 9 for A. gossypii [21], 10 for N. lugens [22] and 12 for S. furcifera [24]. However, nearly twice as many OBPs have been reported for A. lucorum (38 OBPs) and L. lineolaris (33 OBPs) [19, 25]. For CSPs, S. avenae is fewer than the other species in hemipteran, for example, 13 in A. pisum [18] and 9 in A. gossypii [21]. Surprisingly, both OBP11 and OBP12 which were found in A. pisum were not found in our S. avenae antennal transcriptome. This similar phenomenon also appeared in the antennal transcriptome analysis of A. gossypii in which OBP1 and CSP3 were not found [21]. The comprehensive factors of a complicated and heterogeneous host environment are likely responsible for these differences [18] or obvious seasonal alterations in host ranges might lead to the different number of genes for additional and more complex functions [50]. In addition, the cDNA libraries were constructed from antennae total RNA, genes with other functions expressed in other tissues might not appear, and some OBPs and CSPs may not yet have been identified in S. avenae as a result of the limited ESTs that were sequenced [21].
In the phylogenetic tree, most identified OBPs were clustered in highly conserved groups from different aphid species; the same family of OBPs segregated into different central clusters, distributed equally throughout the phylogenetic tree [51]. The conservation of the CSPs among different aphids might be low, only one of five SaveCSPs, SaveCSP3, clustered with MperCSP5 with over 70% identity in an amino acid sequence. The distribution of candidate orthologs in other aphid species suggests that these genes originating from a common ancestor may have similar functions or that they acquired novel functions via subfunctionalization [25]. These results suggest that OBP and CSP proteins in hemipteran insects undergo extensive gene duplication and divergence by natural selection, strongly indicating that they have diverse functions [24].
Antennae-restricted expression combined with age-restricted expression is a useful criterion to identify genes involved in specific chemoreception functions. Our study revealed that consistently high expression of SaveOBP1, SaveOBP2, SaveOBP3, SaveOBP5 and SaveOBP7 in all aphid instars suggests that the respective proteins have a basic conserved role, such as in feeding or identifying general volatiles. High expression of SaveOBP6 in wingless adults might suggest a functional role in the discrimination of egg-laying substrates and the perception of molecules related to new host-plant location [22]. SaveOBP8 and SaveOBP10 with high expression in winged adults might be involved in the perception of molecules related to new host-plant location [52].
Hemipterans typically have a high percentage of antennal expressed OBPs, for example, 5 of 9 OBPs in A. gossypii [21], 12 of 14 OBPs in A. lineolatus [20], at least 6 of 10 OBPs in N. lugens [22], and 21 of 33 OBPs in L. lineolaris [19]. Seven of 13 SaveOBPs (SaveOBP6, SaveOBP7, SaveOBP9, SaveOBP10, SaveOBP13, SaveOBP14 and SaveOBP15) were uniquely or primarily expressed in antennae compared with other tissues, indicating a vital olfactory role for these genes.
SaveOBP2 was dominantly expressed in heads of both winged and wingless aphids; it could be involved in gustatory function in insects [26, 27]. The elevated expression of SaveOBP1 in the abdomens of winged and wingless aphids could indicate a role in the storage and release of chemical compounds from specialized glands localized in the abdomen [25]. OBPs expressed in taste sensilla on legs have been reported to have a role in the behavioral adaptation of Drosophila sechellia [53, 54]. SaveOBP8 was abundantly expressed in legs of winged aphids, which could be related to the adaptation of S. avenae during migration.
According to previous reports on aphid olfaction, some SaveOBP genes showed the had expression profiles that were similar to orthologous genes in A. gossypii and A. pisum. SaveOBP6, AgosOBP6 and ApisOBP6 clustered in a branch with sequences highly expressed in antennae, which strongly suggests that these genes have the same or similar function [21, 55]. The same phenomena can also be found in the branches of SaveOBP7 and ApisOBP7, SaveOBP9 and AgosOBP9, SaveOBP10 and AgosOBP10 [21].
Synthesis of such a high concentration (up to 10 mM) of OBP proteins in insect antennae will consume a lot of energy, so the OBPs should also have some important physiological functions in addition to participating in initial recognition of olfactory signals. On the other hand, OBP expression in non-olfactory tissues has already been ascertained and suggests that they also may function as carriers of chemicals during developmental and physiological processes [25, 56–61]. Overall, OBPs have a very complex expression profile, both relative to body tissues and developmental stages, presumably in connection with their different roles in aphid behaviors [25].
These broad and diverse expression patterns also suggest that different CSPs serve varied functions, including chemosensation [62] and development [63], as well as other processes [37]. SaveCSP1, which was highly expressed in the antennae, could be involved in insect chemoreception. SaveCSP2 and SaveCSP5, both of which are higher expressed in legs, might participate in the process of taste or volatile reception or be indicative of olfactory sensilla on the legs [53, 54]. SaveCSP4 was widely expressed in chemosensory and non-chemosensory tissues. According to the phylogenetic tree, SaveCSP4 grouped in a clade with NlugCSP7, which suggests these two genes may have similar functions in physiological processes other than olfaction [39]. The expression levels of CSPs were affected by age, although without a clear pattern.
Considering all the quantitative qPCR results, it is clear that the olfactory system in winged adults of S. avenae differs across development stages. For example, SaveOBP10 and SaveCSP1 have relatively higher expression levels in winged adult antennae than in other tissues. In addition, expression of SaveOBP10 was highest in winged adults compared to the other developmental stages. SaveOBP8 and SaveCSP5 were also expressed abundantly in legs of winged aphids, about 40- and 20- fold more than in other tissues of winged or wingless aphids.
In summary, using next generation sequencing data for S. avenae, we identified 13 SaveOBP and 5 SaveCSP transcripts. The comprehensive comparison of expression patterns forms a basis for functional studies, especially in revealing major olfactory organ expression of OBP and CSP genes. On the basis of these data, biochemical analyses and behavioral studies will be done to better understand the significant diversity in the functional roles of one or more olfactory genes in the perception of a specific odor. A better understanding of the insect olfactory system and possible targets for insect pest control [47] should inform searches for eco-friendly pest control alternatives to conventional pesticides.
Supporting Information
S1 Fig. Length distribution of transcripts and unigenes.
https://doi.org/10.1371/journal.pone.0161839.s001
(TIF)
S2 Fig. Gene function classification (GO) of unigenes.
https://doi.org/10.1371/journal.pone.0161839.s002
(TIF)
S4 Fig. KEGG pathway classification of unigenes.
A: Cellular processes; B: Environmental information processing; C: Genetic information processing; D: Metabolism; E: Organismal systems.
https://doi.org/10.1371/journal.pone.0161839.s004
(TIF)
S5 Fig. OBP and CSP bands after qPCR in S. avenae.
https://doi.org/10.1371/journal.pone.0161839.s005
(TIF)
S2 Table. GenBank accession numbers of amino acid sequences used to construct the phylogenetic tree.
https://doi.org/10.1371/journal.pone.0161839.s007
(DOCX)
S5 Table. Sequence identity between SaveOBPs.
Calculations are based on amino acid sequence alignment by DNAMAN. The percentage identity for each pair is shown.
https://doi.org/10.1371/journal.pone.0161839.s010
(DOCX)
S6 Table. Sequence identity between SaveCSPs.
Calculations are based on amino acid sequence alignment by DNAMAN. The percentage identity of each pair is shown.
https://doi.org/10.1371/journal.pone.0161839.s011
(DOCX)
Acknowledgments
Our study was financially supported by the National Natural Science Foundation of China (31371946, 31401740), The Agricultural Science and Technology Innovation Program (ASTIP), The National S& T Major Project (2014ZX0800201B), and the Cooperation Project between China and Belgium from Ministry of Science and Technology (MOST) (2014DFG32270).
Author Contributions
- Conceptualization: WXX JLC.
- Data curation: JLC.
- Formal analysis: WXX YZ.
- Investigation: WXX.
- Methodology: WXX JF.
- Project administration: JLC JF.
- Resources: WXX JRS.
- Supervision: JLC.
- Validation: WXX QXX.
- Writing – original draft: WXX JF ZLH.
- Writing – review & editing: JLC JF.
References
- 1.
Blackman RL, Eastop VF. Aphids on the world herbaceous plants and shrubs: An identification and information guide. Chichester: Wiley; 2006.
- 2.
Wu JX. Agricultural Entomology. Beijing: China Agriculture Press; 2002.
- 3. Chi BJ, Zhu YF, Axel V, Chen JL, Liu Y. Demographic and quantitative food web analysis of Sitobion avenae and its natural enemies. Chinese Journal of Applied Entomology. 2014; 51(6): 1496–1503.
- 4.
Chen JL. Wheat aphids and control. Golden Shield Press, 2013.
- 5. Fan J, Liu Y, Zeng JG, Guo M, Sun JR, Cheng P, et al. Advancement of new prevent and control technologies for aphids in wheat and vegetable. Chinese Journal of Applied Entomology. 2014; 51(6): 1413−1434.
- 6. Visser JH, Piron PGM, Hardie J. The aphids’ peripheral perception of plant volatiles. Entomologia Experimentalis et Applicata. 1996; 80: 35–38.
- 7. Bruyne M, Baker TC. Odor detection in insects: volatile codes. J Chem Ecol. 2008; 34: 882–897. pmid:18535862.
- 8. Dahanukar A, Hallem EA, Carlson JR. Insect chemoreception. Curr Opin Neurobiol. 2005; 15: 423–430. pmid:16006118.
- 9. Tegoni M, Campanacci V, Cambillau C. Structural aspects of sexual attraction and chemical communication in insects. Trends Biochem Sci. 2004; 29:257–264. pmid:15130562.
- 10. Visser JH. Host odor perception in phytophagous insects. Annu Rev Entomol. 1986; 31: 121–144.
- 11.
Yan FM, Chen JL, Zhang YJ, Tang QB, Zhou HB. Advances in Insect Chemical Ecology. In “Report on Advances in Plant Protection during 2012–2013”, Edited by Guo Yuyuan. Chinese Science and Technology Press.; 2013.
- 12. Vogt RG, Riddiford LM. Pheromone binding and inactivation by moth antennae. Nature. 1981; 293: 161–163. pmid:18074618.
- 13. Klein U. Sensillum-lymph proteins from antennal olfactory hairs of the moth Antheraeapoly phemus (Saturniidae). Insect Biochem. 1987; 17: 1193–1204.
- 14. Pelosi P, Zhou JJ, Ban LP, Calvello M. Soluble proteins in insect chemical communication. Cell Mol Life Sci. 2006; 63: 1658–1676.
- 15. Zhou JJ. Odorant-binding proteins in insects. Vitam Horm. 2010; 83: 241–272. pmid:20831949.
- 16. Honson NS, Gong Y, Plettner E. Structure and function of insect odorant and pheromone-binding proteins (OBPs and PBPs) and chemosensory-specific proteins (CSPs). Recent Adv Phytochem. 2005; 39: 227–268.
- 17. Sachse S, Krieger J. Olfaction in insects. The primary processes of odor recognition and coding. e-Neuroforum. 2011; 2: 49–60.
- 18. Zhou JJ, Vieira FG, He XL, Smadja C, Liu R, Rozas J, et al. Genome annotation and comparative analyses of the odorant-binding proteins and chemosensory proteins in the pea aphid Acyrthosiphon pisum. Insect Mol Biol. 2010; 19 (Suppl. 2): 113–122. pmid:20482644.
- 19. Hull JJ, Perera OP, Snodgrass GL. Cloning and expression profiling of odorant-binding proteins in the tarnished plant bug, Lygus lineolaris. Insect Mol Biol. 2014; 23: 78–97. pmid:24224606
- 20. Gu SH, Wang SP, Zhang XY, Wu KM, Guo YY, Zhou JJ, et al. Identification and tissue distribution of odorant binding protein genes in the lucerne plant bug Adelphocoris lineolatus (Goeze). Insect Biochem Mol Biol. 2011; 41: 254–263. pmid:21232599.
- 21. Gu SH, Wu KM, Guo YY, Field LM, Pickett JA, Zhang YJ, et al. Identification and expression profiling of odorant binding proteins and chemosensory proteins between two wingless morphs and a winged morph of the cotton aphid Aphis gossypii Glover. PLoS ONE. 2013; 8: e73524. pmid:24073197.
- 22. Zhou SS, Sun Z, Ma W, Chen W, Wang MQ. De novo analysis of the Nilaparvata lugens (Stål) antenna transcriptome and expression patterns of olfactory genes. Comp Biochem Phys D. 2014; 9: 31–39. pmid:24440828.
- 23. He P, Zhang J, Liu NY, Zhang YN, Yang K, Dong SL. Distinct expression profiles and different functions of odorant-binding proteins in Nilaparvata lugens (Stål). PLoS ONE. 2011; 6(12): e28921. pmid:22174925
- 24. He M, He P. Molecular characterization, expression profiling, and binding properties of odorant binding protein genes in the white backed planthopper, Sogatella furcifera. Comp Biochem Phys B. 2014; 174: 1–8. pmid:24815350.
- 25. Yuan HB, Ding YX, Gu SH, Sun L, Zhu XQ, Liu HW, et al. Molecular characterization and expression profiling of odorant-binding proteins in Apolygus lucorum. PLoS ONE. 2015; 10(10): e0140562. pmid:26466366
- 26. Shanbhag SR, Park SK, Pikielny CW, Steinbrecht RA. Gustatory organs of Drosophila melanogaster: fine structure and expression of the putative odorant-binding protein PBPRP2. Cell Tissue Res. 2001; 304: 423–437. pmid:000169308900012.
- 27. Jeong YT, Shim J, Oh SR, Yoon HI, Kim CH, Moon SJ, et al. An odorant-binding protein required for suppression of sweet taste by bitter chemicals. Neuron. 2013; 79: 725–737.
- 28. Wogulis M, Morgan T, Ishida Y, Leal WS, Wilson DK. The crystal structure of an odorant binding protein from Anopheles gambiae: evidence for a common ligand release mechanism. Biochemical and Biophysical Research Communications. 2006; 339: 157–164. pmid:16300742.
- 29. Qiao HL, Tuccori E, He XL, Gazzano A, Field L. Discrimination of alarm pheromone (E)-β-farnesene by aphid odorant-binding proteins. Insect Biochem Mol Biol. 2009; 39: 414–419. pmid:19328854.
- 30.
Verheggen FJ, Haubruge E, Mescher MC. Chapter nine: Alarm pheromones-chemical signaling in response to danger. Vitam Horm, Academic Press. 2010; 83: 215–239. https://doi.org/10.1016/S0083-6729(10)83009-2 pmid:20831948.
- 31. Vandermoten S, Mescher MC, Francis F, Haubruge E, Verheggen FJ. Aphid alarm pheromone: An overview of current knowledge on biosynthesis and functions. Insect Biochem Mol Biol. 2012; 42: 155–163. pmid:22178597.
- 32. Pickett JA, Woodcock CM, Midega CAO, Khan ZR. Push-pull farming systems. Curr Opin Biotechnol. 2014; 26: 125–132. pmid:24445079.
- 33. Zhong T, Yin J, Deng SS, Li KB, Cao YZ. Fluorescence competition assay for the assessment of green leaf volatiles and trans-β-farnesene bound to three odorant-binding proteins in the wheat aphid Sitobion avenae (Fabricius). J Insect Physiol. 2012; 58: 771–781. pmid:22306433.
- 34. Jacquin-Joly E, Vogt RG, Francois MC, Nagnan-Le Meillour P. Functional and expression pattern analysis of chemosensory proteins expressed in antennae and pheromonal gland of Mamestra brassicae. Chem Senses. 2001; 26: 833–44. pmid:11555479.
- 35. Liu R, He X, Lehane S, Lehane M, Hertz-Fowler C, Berriman M, et al. Expression of chemosensory proteins in the tsetse fly Glossina morsitans morsitans is related to female host-seeking behaviour. Insect Mol Biol. 2012; 21: 41–8. pmid:22074189
- 36. Kitabayashi AN, Arai T, Kubo T, Natori S. () Molecular cloning of cDNA for p10, a novel protein that increases in the regenerating legs of Periplaneta Americana (American cockroach). Insect Biochem Mol Biol. 1998; 28: 785–90. pmid:9807224
- 37. Maleszka J, Foret S, Saint R, Maleszka R. () 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. 2007; 217: 189–196. pmid:17216269.
- 38. Liu YL, Guo H, Huang LQ, Pelosi P, Wang CZ. Unique function of a chemosensory protein in the proboscis of two Helicoverpa species. J Exp Biol. 2014; 217: 1821–1826. pmid:24625642
- 39. Yang K, He P, Dong SL. Different expression profiles suggest functional differentiation among chemosensory proteins in Nilaparvata lugens (Hemiptera: Delphacidae). Journal of Insect Science. 2014; 14: 270. pmid:25527582.
- 40. Jacobs SP, Liggins AP, Zhou JJ, Pickett JA, Jin X, Field LM. OS-D-like genes and their expression in aphids (Hemiptera: Aphididae). Insect Mol Biol. 2005; 14: 423–432. pmid:16033435.
- 41. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nature Biotechnology. 2011; 29: 644–652. pmid:21572440.
- 42. Zhou JJ, He XL, Pickett JA, Field LM. Identification of odorant-binding proteins of the yellow fever mosquito Aedesa egypti: genome annotation and comparative analyses. Insect Mol Biol. 2008; 17: 147–163. pmid:18353104.
- 43. Zhou JJ, Kan YC, Antoniw J, Pickett JA, Field LM. Genome and EST analyses and expression of a gene family with putative functions in insect chemoreception. Chem Senses. 2006; 31: 453–465. pmid:16581978.
- 44. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nature Methods. 2008; 5: 621–628. pmid:18516045
- 45. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011; 28: 2731–2739. pmid:21546353.
- 46. Yang C, Pan H, Liu Y, Zhou X. Selection of reference genes for expression analysis using quantitative real-time PCR in the pea aphid, Acyrthosiphon pisum (Harris) (Hemiptera, Aphidiae). PLoS ONE. 2014; 9(11): e110454. pmid:25423476
- 47. Vandesompele J, Katleen DP, Bruce P, Nadine VR, Anne DP, Rank S. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biology. 2002; 3(7): pmid:12184808.
- 48. Xu YL, He P, Zhang L, Fang SQ, Dong SL, Zhang YJ, et al. Large-scale identification of odorant-binding proteins and chemosensory proteins from expressed sequence tags in insects. BMC Genomics. 2009; 10: 632–645. pmid:20034407.
- 49. Zhou JJ, Huang WS, Zhang GA, Pickett JA, Field LM. “Plus-C” odorant binding protein genes in two Drosophila species and the malaria mosquito Anopheles gambiae. Gene. 2004; 327: 117–129. pmid:14960367.
- 50. Pan HS, Liu B, Lu YH, Wyckhuys KAG. Seasonal alterations in host range and fidelity in the polyphagous mirid bug, Apolygus lucorum (Heteroptera: Miridae). PLoS ONE. 2015; 10: e0117153. pmid:25692969.
- 51. Vogt RG, Callahan FE, Rogers ME, Dickens JC. Odorant binding protein diversity and distribution among the insect orders, as indicated by LAP, an OBP-related protein of the true bug Lygus lineolaris (Hemiptera, Heteroptera). Chem Senses. 1999; 24: 481–495. pmid:000083202400003.
- 52. Pickett JA, Wadhams LJ, Woodcock CM, Hardie J. The chemical ecology of aphids. Annu Rev Entomol. 1992; 37: 67–90.
- 53. Yasukawa J, Tomioka S, Aigaki T, Matsuo T. Matsuo Evolution of expression patterns of two odorant-binding protein genes, Obp57d and Obp57e, in Drosophila. Gene. 2010; 467(1–2): 25–34. pmid:20637846.
- 54. Harada E, Nakagawa J, Asano T, Taoka M, Sorimachi H, Ito Y, et al. Functional evolution of duplicated odorant-binding protein genes, Obp57d and Obp57e, in Drosophila. PLoS ONE 2012; 7: p. e29710. pmid:22238638
- 55. Filomena DB, Lea R, Daniele B, Annalisa G, Terenzio C, Sun YF, et al. Expression pattern analysis of odorant-binding proteins in the pea aphid Acyrthosiphon pisum. Insect Science. 2015; 22: 220–234. pmid:24591440.
- 56. Zhou SH, Zhang J, Zhang SG, Zhang L. Expression of chemosensory proteins in hairs on wings of Locusta migratoria (Orthoptera: Acrididae). J Appl Entomol. 2008; 132: 439–450.
- 57. Dani FR, Michelucci E, Francese S, Mastrobuoni G, Cappellozza S, La Marca G, et al. Odorant-binding proteins and chemosensory proteins in pheromone detection and release in the silkmoth Bombyx mori. Chem Senses. 2001; 36: 335–344. pmid:21220518.
- 58. Calvello M, Guerra N, Brandazza A, D’Ambrosio C, Scaloni A, Dani FR, et al. Soluble proteins of chemical communication in the social wasp Polistes dominulus. Cellular and Molecular Life Sciences. 2003; 60: 1933–1943. pmid:14523553.
- 59. Li S, Picimbon JF, Ji SD, Kan YC, Qiao CL, Zhou JJ, et al. Multiple functions of an odorant binding protein in the mosquito Aedes aegypti. Biochemical and Biophysical Research Communications. 2008; 372: 464–468. pmid:18502197.
- 60. Vogel H, Heidel AJ, Heckel DG, Groot AT. Transcriptome analysis of the sex pheromone gland of the noctuid moth Heliothis virescens. BMC Genomics. 2010; 11: 29. pmid:20074338.
- 61. Iovinella I, Dani FR, Niccolini A, Sagona S, Michelucci E, Gazzano A, et al. Differential expression of odorant-binding proteins in the mandibular glands of the honey bee according to caste and age. Journal of Proteomic Research. 2011; 10: 3439–3449. pmid:21707107.
- 62. Sun YL, Huang LQ, Pelosi P, Wang CZ. Expression in antennae and reproductive organs suggests a dual role of an odorant-binding protein in two sibling Helicoverpa species. PLoS ONE. 2012; 7: e30040. pmid:22291900.
- 63. González D, Zhao Q, McMahan C, Velasquez D, Haskins WE, Sponsel V, et al. The major antennal chemosensory protein of red imported fire ant workers. Insect Mol Biol. 2009; 18: 395–404. pmid:19523071.