Figures
Abstract
The oriental fruit fly, Bactrocera dorsalis, is a devastating fruit fly pest in tropical and sub-tropical countries. Like other insects, this fly uses its chemosensory system to efficiently interact with its environment. However, our understanding of the molecular components comprising B. dorsalis chemosensory system is limited. Using next generation sequencing technologies, we sequenced the transcriptome of four B. dorsalis developmental stages: egg, larva, pupa and adult chemosensory tissues. A total of 31 candidate odorant binding proteins (OBPs), 4 candidate chemosensory proteins (CSPs), 23 candidate odorant receptors (ORs), 11 candidate ionotropic receptors (IRs), 6 candidate gustatory receptors (GRs) and 3 candidate sensory neuron membrane proteins (SNMPs) were identified. The tissue distributions of the OBP and CSP transcripts were determined by RT-PCR and a subset of nine genes were further characterized. The predicted proteins from these genes shared high sequence similarity to Drosophila melanogaster pheromone binding protein related proteins (PBPRPs). Interestingly, one OBP (BdorOBP19c) was exclusively expressed in the sex pheromone glands of mature females. RT-PCR was also used to compare the expression of the candidate genes in the antennae of male and female B. dorsalis adults. These antennae-enriched OBPs, CSPs, ORs, IRs and SNMPs could play a role in the detection of pheromones and general odorants and thus could be useful target genes for the integrated pest management of B. dorsalis and other agricultural pests.
Citation: Wu Z, Zhang H, Wang Z, Bin S, He H, Lin J (2015) Discovery of Chemosensory Genes in the Oriental Fruit Fly, Bactrocera dorsalis. PLoS ONE 10(6): e0129794. https://doi.org/10.1371/journal.pone.0129794
Academic Editor: Joseph Clifton Dickens, United States Department of Agriculture, Beltsville Agricultural Research Center, UNITED STATES
Received: February 15, 2015; Accepted: May 13, 2015; Published: June 12, 2015
Copyright: © 2015 Wu 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 relevant data are within the paper and its Supporting Information files.
Funding: This work was supported by the National Natural Science Foundation of China (31171852).
Competing interests: The authors have declared that no competing interests exist.
Introduction
Chemoreception plays a crucial role in insects such as agricultural pests, disease vectors and social insects. These insects use two sensations, olfaction and gustation, to evaluate and locate food sources, shelter, mates, and oviposition sites as well as to avoid predators and other dangers [1–5]. The major molecular components of insect olfaction include odorant-binding proteins (OBPs), odorant receptors (ORs), ionotropic receptors (IRs), sensory neuron membrane proteins (SNMPs) and odorant-degrading enzyme (ODEs) [6], and the major gustatory or contact chemosensation-related proteins are gustatory receptors (GRs) [7,8]. In addition, chemosensory proteins (CSPs) are also found in olfactory and gustatory organs of insects and are involved in the detection of chemicals [9–14].
Chemosensory proteins are widely used by tephritid fruit flies to locate host plants and thereby cause major losses in fruits and vegetables worldwide. Because of their devastating impact on agriculture they are often the target of intense insecticide applications in order to protect commercial production of agricultural crops. The oriental fruit fly, Bactrocera dorsalis, is the main fruit fly pest in tropical and sub-tropical countries, and is reported to feed on >117 host species, in 76 genera and 37 families [15]. Since sporadic outbreaks of the pest have been reported worldwide, this fly has been the target of global integrated pest management [16,17,18].
Previous studies have shown that B. dorsalis exhibits sexually dimorphic behavior, which is influenced by olfactory cues expressed during specific developmental stages [19,20,21]. When the male flies reach sexual maturity, they are strongly attracted to and compulsively feed on Methyl eugenol (ME) (non-host compounds), and this behavior is used to control B. dorsalis via male annihilation technique through mass trapping [22]. Behavioral assays with electro-physiologically active compounds from mango as a host plant, revealed that γ-octalactone induced oviposition by gravid B. dorsalis females [20]. In addition, the pest control strategy based on the behavioral manipulation of B. dorsalis still relies on ME based male lure [23]. However, the molecular mechanisms of behavior-based pest control are not clearly understood. Thus, there is a need to understand the molecular basis of chemoreception in tephritid fruit flies.
Currently, understanding of the molecular components in the B. dorsalis chemosensory system is limited with only 10 known OBPs [24], and reports that ME increases the gene expression level of OR co-receptor [25] and that BdorCSP2 is involved in the chemoreception of Rhodojaponin-III, an antifeedant [11]. Recently, using “computational reverse chemical ecology,” one B. dorsalis OBP protein (GenBank ID: ACB56577.1) expressed in the antennae of gravid females was shown to be an attractant to semiochemicals [26]. Therefore, systematic research on chemoreception may provide valuable information that could be used for the rapid screening of potential semiochemicals. In this study, we applied a transcriptomic approach to identify a large array of candidate chemosensory genes in B. dorsalis. We used next generation sequencing (454 Life Sciences) and evaluated the presence of chemosensory genes in all the developmental stages of B. dorsalis: eggs collected within 24 h of oviposition; larvae (first, second and third instars); pupae (1 d-old, 4 d-old and 7 d-old pupae) and newly-emerged adult chemosensory tissues, including antenna, leg and head (within six days of eclosion) in a 1:1 female: male ratio. Our results show the presence of a number of chemosensory gene transcripts in B. dorsalis and the presence of antennae-specific OBPs that could be used effectively towards the control of agricultural pests.
Materials and Methods
Ethics Statement
The oriental fruit fly, B. dorsalis is not included in the ‘‘List of Endangered and Protected Animals in China” because it is a major fruit fly pest in tropical and sub-tropical countries. All experiments were performed in compliance with the general ethical guidelines in order to minimize pain and discomfort to the insects.
Insect Rearing
The oriental fruit fly, B. dorsalis, was obtained from a laboratory-reared stock colony (Institute for Management of Invasive Alien Species, Zhongkai University of Agriculture and Engineering, Guangzhou, PR China) maintained at 28°C, 70% relative humidity, and a 14: 10 (L: D) photoperiod for the past 8 years. Adult flies were reared on an artificial diet mix described previously [27], and newly hatched larvae were reared on banana in the laboratory [28].
RNA Extraction, cDNA Library Preparation and Sequencing
Total RNA was isolated from the following developmental stages: eggs collected within 24h of oviposition; larvae (first, second and third instar larvae; ratio 1: 1: 1); pupae (1d-old, 4d-old and 7d-old pupae; ratio 1: 1: 1); and adult chemosensory tissues, including antenna, leg and head (within six days of eclosion; ratio 1: 1: 1) in a 1:1 female: male ratio. All samples were snap-frozen in liquid nitrogen and stored at –80°C until total RNA was extracted. Construction of normalized cDNA libraries from the four B. dorsalis samples, and 454 pyrosequencing were carried out as follows. First, total RNA were extracted from each sample using TRIzol reagent (Life Science Technologies-Invitrogen), and the quantity and quality of RNA were assessed by spectrophotometry and gel electrophoresis. Then, mRNA was isolated from 20 μg of each total RNA using the Oligotex mRNA Mini kit (Qiagen, CA). First strand cDNA was synthetized from 1μg mRNA with SuperScript III reverse transcriptase using dT15VN2 primer (Invitrogen) under the following conditions: 5 minutes at 65°C, 2 minutes at 4°C, 1 h at 42°C and 10 min at 70°C in a PCR machine (Bio-Rad). The second strand was synthesized from 1 μl of the first strand cDNA reaction mix using DNA Ligase, DNA polymerase I and RNaseH from E. coli according to the manufacturer’s instructions (Invitrogen). T4 DNA polymerase was added and incubated for 5 min at 16°C in a PCR machine. The synthesized double stranded cDNA were purified with the QIAquick PCR Purification kit (Qiagen, Valencia, CA, USA), and the yield was determined using the TBS 380 Fluorometer (Turner Biosystems). Subsequently, cDNA was fragmented by sonication and the cDNA samples ranging in size from 100 bp to 800 bp were purified on a 2% agarose gel. Then, DNA concentration in each cDNA sample was determined using the Bioanalyzer DNA1000 kit (Agilent, USA). Each purified cDNA sample was then used to synthesize single-strand template DNA (sstDNA) libraries using the GS20 DNA Library Preparation kit (Roche Applied Science) following the manufacturer's recommendations (1/4 run for each sample). Library quality was assessed on an Agilent Bioanalyzer High Sensitivity DNA chip. Finally, each library was normalized in equimolar concentrations and diluted to 1x106 molecules/μl. Emulsion based clonal amplification and sequencing were performed on the 454 Genome Sequencer FLX Titanium system according to the manufacturer’s instructions (454 Life Sciences, Branford, CT). The raw data from 454 reads are deposited in the NCBI Short Read Archive under the accession numbers SRX862648, SRX862768, SRX862771, SRX862773, respectively.
De novo Assembly
The raw 454 sequences in SFF files were extracted using the Python script sff_extract.py developed by COMAV (http://bioinf.comav.upv.es). All the raw sequences were then processed to remove low quality and adaptor sequences using programs SeqClean Lastest86_64 [29], Newbler 2.5.3 [30] and LUCY 1.20p [31]. The resulting sequences were then screened against the NCBI UniVec database and bacterial genome sequences to remove possible contaminants. Cleaned reads shorter than 50 bp were discarded. De novo assembly of the high quality 454 sequences from each B. dorsalis sample was performed by Newbler version 2.5.3 using default parameters under the cDNA option (Roche, Branford, CT, USA).
Gene Annotation
Amino acid sequences predicted from the assembled 454 sequences were compared to protein sequences in the NCBI non-redundant (nr) protein database on a local server using the BLASTALL program with the cutoff e value of 10−5 [32]. GO annotation was performed using Blast2GO. GO association was done by BLASTX comparison against the NCBI nr database [33,34]. To specifically annotate OBPs, CSPs, ORs, IRs, GRs and SNMPs in B. dorsalis, assembled sequences were analyzed using TBLASTN and TBLASTX programs against custom-made databases consisting of insect sequences processed using the BioEdit program [35]. Sequences whose best TBLASTN hits corresponded to OBPs, CSPs, ORs, IRs, GRs and SNMPs were then retained as candidate B. dorsalis chemosensory transcripts and their translation was manually verified and corrected if needed. Finally, families of all candidate B. dorsalis chemosensory protein sequences were analyzed on Pfam [36].
RACE-PCR, Cloning and Sequence Analysis
To obtain the full-length coding sequences of the candidate transcripts, the SMART RACE-PCR kit (Clontech) was used with gene-specific primers (S1 Table) designed using Primer Premier 6 (PREMIER Biosoft International, CA, USA) following the manufacturer’s instructions. The amplified products were separated on a 2% agarose gel prior to purifying the products using the Agarose Gel DNA Purification Kit (TAKARA, China). The amplified fragments were then cloned into the pMD20-T vector (TAKARA, China) and sequenced from both directions. Then, open reading frames (ORF) in the assembled full-length unigenes were identified using the ORF finder (http://www.ncbi.nlm.nih.gov/gorf/gorf.html). The signal peptide of OBPs and CSPs were predicted using SignalP 4.0 [37]. Transmembrane domains of candidate ORs, IRs, GRs and SNMPs were predicted using TMHMM 2.0 [38]. The deduced protein sequences were further confirmed by searching the Pfam database with default parameters and e-value 1.0 [39]. Based on these searches, putative chemosensory genes in the B. dorsalis transcriptome were named after their Drosophila homologues. Transcripts with the highest similarity to the same Drosophila genes were differentiated with a numerical postscript (S2 Table).
Comparative Analysis of Chemosensory Genes between B. dorsalis Developmental Stages
Following assembly, transcripts were assigned an RPKM [40] value based on the number of uniquely mapping reads aligning to each transcript using SOAP software (release 2.21). The RPKM of chemosensory gene transcripts from eggs, larvae, pupae and adults were compared the differential expression of chemosensory genes in the various developmental stages.
Phylogenetic Analyses
Phylogenetic analyses of the B. dorsalis chemosensory genes were reconstructed based on the amino sequences after removal of the signal peptides and the data set collected from NCBI. The OBP data set contained 52 sequences from D. melanogaster [41,42,43], 16 sequences from Ceratitis capitata [44], 15 OBPs from R. pomonella [45] and 9 OBPs from R. suavis [46]. The CSP data set contained 47 sequences from 12 Drosophila sp. [41] and 5 sequences from Glossina morsitans morsitans [14]. The OR data set contained 63 OR sequences from D. melanogaster [47,48] and 76 ORs from A. gambiae [49]. The iGluR and IR data sets contained 66 IR sequences from D. melanogaster and 55 IR sequences from A. gambiae [50]. The SNMP data set contained 26 SNMP sequences identified in Diptera and Lepidoptera [51].
For all proteins analyzed, their respective amino acid sequences were aligned using MAFFT v.6 (E-INS-I parameter set for OBPs and CSPs; FFT-NS-2 parameter set for ORs and IRs) [52]. For each data set, the best-fit model of protein evolution was selected by MEGA 6.0 using the Akaike information criterion (the LG+I+G model of OBP data set; the LG+G model of CSP data set; the LG+I+F model of OR data set; the LG+G model of IR data set; the LG+I+G model of SNMP data set). Dendrograms were calculated using maximum likelihood analysis with MEGA 6.0 [53] with both SPR (Subtree Pruning and Regrafting) and MP (Maximum Parsimony) methods for tree topology improvement. Robustness of the branches was assessed with 1000 bootstrap pseudo-replicates. Dendrograms were viewed and edited in FigTree (http://tree.bio.ed.ac.uk/software/figtree/).
Analysis of Chemosensory Gene Expression by RT-PCR
RT-PCR was employed to investigate and compare the expression of candidate chemosensory genes in different B. dorsalis tissues. Total RNA from different tissues was extracted as described above and treated with DNase I (TAKARA, China) to remove trace amounts of genomic DNA. Then, first strand cDNA synthesized using the First strand cDNA synthesis kit (TAKARA, China) was used as a template in PCR reactions with gene-specific primers designed using Primer3 (http://primer3.ut.ee/) [54] (S3 Table). The B. dorsalis α-tublin gene (GenBank Acc. GU269902) was used as the control [55]. Because OBPs and CSPs are not restricted to the olfactory tissues and are known to participate in other physiological functions [56,57,58,59,60,61], two schemes of RT- PCR analyses (RT-PCR and qRT-PCR) were used (S1 Text). Each RT-PCR was repeated three times using three independently isolated RNA samples. PCR amplification products were separated on a 1.5% agarose gel and verified by direct DNA sequencing (Invitrogen, China).
qRT-PCR Analysis
qRT-PCR was used to quantify expression levels of the OBP genes and CSP genes that were antennae-rich or antennae-specific. Total RNA isolated from 100 antennae pairs and two whole bodies of male and female flies was used to synthesize first strand cDNA as described above. Primers for the OBP genes were designed using Primer3 (S3 Table). The reactions were performed with 2 μl of the cDNA as template in a LightCycler 480 System (Roche Applied Science) using the SYBR Premix EX Taq (TAKARA, China). The B. dorsalis α-tublin (α-TUB) (GenBank Acc. GU269902) was used as an internal control for normalization. Negative controls without cDNA template or transcriptase were included in each experiment. To check reproducibility, each qRT-PCR reaction had three technical replicates and three biological replicates. Relative expression of the genes in the various tissues was estimated using the 2 -ΔΔCT method [62]. Statistical analyses of the relative expression data were performed using Prism 5.0 (GraphPad Software, CA). Statistical significance of the temporal expression was analyzed by ANOVA followed by a Tukey multiple comparison test. A value of P<0.05 was considered statistically significant.
Results
Sequencing, Assembly and Annotation
The transcriptome of B. dorsalis eggs, larvae, pupae and adult chemosensory tissues generated using the GS/FLX 454 technology yielded a total of 1,122,242 raw reads (Table 1). After assembly, 22,934 unigenes were generated in each sample (see Table 1). Among the unigenes, the majority of B. dorsalis transcripts were assigned to the “binding” and “catalytic activity” in the molecular function GO category (Fig 1) in all four developmental stages with each category having 1561 to 2280 reads (S4 Table). The category “Transporter Activity” had 117 to 201 reads while the rest including “Receptor Activity,” “Molecular Transducer Activity,” “Protein Binding Transcription Factor Activity,” “Nucleic Acid binding Transcription Factor Activity,” “Structural Molecule Activity,” “Enzyme Regulator Activity,” and “Electron Carrier Activity Antioxidant Activity” had less than 100 reads each (S4 Table).
Identification and Characterization of Chemosensory Genes
In order to identify chemosensory genes from B. dorsalis, we performed transcriptome sequencing of four developmental stages in B. dorsalis (eggs, larvae, pupae and adults chemosensory tissues). Using homologous searches, a total of 78 putative chemosensory genes were identified, including 31 OBPs, 4 CSPs 23 ORs, 11 IRs, 6 GRs and 3 SNMPs (S5 Table and S6 Table, fasta format file in S2 Text). The complete mRNA sequences of these genes were obtained by performing RACE-PCR.
Among these were 8 OBPs (GenBank accession no. KC559112.1, KC559113.1, KC559114.1, KC559117.1, KC559118.1, KC559119.1, KC559121.1) [24] and 1 CSP (KC897022) [11,24] previously reported in B. dorsalis. However, our transcriptomes did not contain three previously reported B. dorsalis OBPs (KC559115.1, KC559116.1 and KC559120.1). The remaining OBP and CSP sequences were considered as unique and their predicted protein sequences were named as BdorOBPs and BdorCSPs, respectively.
The length of the complete BdorOBPs ranged from 134 (BdorOBP56e) to 274 amino acids (BdorOBP83ef). Prediction of signal peptide revealed that these BdorOBPs were secretory proteins, except for three (BdorOBP56e, BdorOBP57c and BdorOBP69a), which were intracellular proteins. Based on the presence of conserved cysteine profiles, 20 BdorOBPs were classified as classic OBPs, with the six conserved cysteine residues characteristic to insect OBPs [43,63]; five were identified as minus-C OBPs, which encoded putative polypeptides with four or five conserved cysteine residues [43]; two were Plus-C OBPs, with two conserved cysteines plus one proline; two were dimer OBPs, with two six-cysteine signatures; and two were atypical OBPs, with 9–10 cysteines and a long C-terminus (S5 Table). Among the classic OBPs, BdorOBP84a-1 was slightly different with one extra cysteine residue located between C5 and C6, despite having a highly conserved OBP secondary structure. Similarly, BdorOBP50c and BdorOBP50e had the general characteristics of Plus-C but also had one additional cysteine before two successive conserved cysteine residues (S5 Table).
The length of the complete BdorCSPs ranged from 111 (BdorCSP4) to 156 amino acids (BdorCSP2). Four BdorCSPs were predicted as secretory proteins and all four had the characteristics of insect CSP gene families, with four high cysteine profiles.
Among the ORs, BdorORCO, BdorOR7a-1 and BdorOR63a-2a belonged to the highly conserved OR family with seven transmembrane domains, which is a characteristic of insect ORs, while the others belonged to a divergent member of the OR family with 4, 5 and 8 transmembrane domains (S6 Table). Among the IRs, BdorIR41a, BdorIR75d and BdorIR100a belonged to the highly conserved ionotropic glutamate receptors (iGluRs) family with three transmembrane domains (S6 Table), and others belonged to a divergent group of IRs. Bioinformatic analysis identified 6 candidate GRs. The insect GRs contained seven transmembrane domains. TMHMM2.0 predicted 1 candidate GR (BdorGr63a) with seven transmembrane domains (S6 Table). Not surprisingly, three SNMP sequences, BdorSNMP1-1, BdorSNMP1-2 and BdorSNMP2 that shared high amino acid identity (86 to 91%) with the conserved insect CD36 family were also identified.
Expression profile of Chemosensory Genes
Based on the RPKM value, it was evident that all these candidate odorant-binding proteins were expressed at a high level, while the candidate chemosensory membrane proteins were expressed at a low level in the different development stages (S7 Table). To provide functional clues, the tissue distribution of the OBP and CSP gene transcripts were examined using semi-quantitative RT-PCR (Fig 2). The results showed the robust expression of 9 OBP genes (BdorOBPlush, BdorOBP19a, BdorOBP56h, BdorOBP69a, BdorOBP83a-1, BdorOBP83a-2, BdorOBP84a-1 and BdorOBP84a-2) and BdorCSP3 exclusively in the male and female antennae. Although the PBRP homologs with D. melanogaster, BdorOBP19d-1, BdorOBP19d-2 and BdorOBP28, were present in the antennae, they were also abundant in other tissues (Fig 2). Based on such differential expression of these OBPs and CSPs in male and female antennae, we presume that these OBPs in males and females may be involved in detecting sex pheromones while in females they may play important roles in locating suitable host plants and oviposition sites. We further quantified the B. dorsalis OBP and CSP gene transcripts using qRT-PCR and compared their expression levels in different tissues between sexes. The results suggested that 3 OBP gene transcripts (BdorOBPlush, BdorOBP83a-2, BdorOBP84a-1) were expressed higher in the male antennae than in the female antennae, while BdorOBP19a was expressed higher in the female antennae than in the male antennae (p<0.01) (Fig 3). In addition, the results suggested that 4 OBP gene transcripts BdorOBP69a, BdorOBP19d-1 and BdorOBP83a-1 and BdorCSP3 were expressed equally in adult male and female antennae. However, two OBPs, BdorOBP56e and BdorOBP83cd, were not found in adult tissues but rather in other developmental stages i.e. eggs, larvae and pupae (S1 Fig). Interestingly, BdorOBP19c appears to be a pheromonal gland OBP present in the rectal glands of mature females, but not in immature female rectal glands (Fig 2 and S1 Fig).
MA: male antennae, MB: male body, FA: female antennae, FB: female body. Error bars represent standard error. Different letters (a, b) above each bar denote significant differences between samples (p<0.05).
Expression patterns of 23 ORs, 3 SNMPs, 11 IRs and 2 GRs in male antennae, female antennae and legs were analyzed by semi-quantitative RT-PCR. The results demonstrated that all 23 ORs, 2 SNMPs (BdorSNMP1-1 and BdorSNMP1-2), 6 IRs (BdorIR40a, BdorIR41a, BdorIR75d, BdorIR76b, BdorIR84s and BdorIR92a), and 2 GRs (GR63 and GR21) were specifically expressed in antennae (Fig 4) suggesting that these receptors may play important roles in the detection of odorants. RT-PCR analyses also showed that the receptors for all 33 above proteins (ORs, SNMPs, GRs and IRs) were expressed in both sexes with some differentially expressed in male or female antennae.
MA: male antennae, FA: female antennae, L: legs.
Phylogenetic Analysis of Chemosensory Genes
In order to assign putative functions to chemosensory genes, we determined the phylogenetic relationship between the 31 BdorOBPs identified in this study, and 52 OBPs previously reported in D. melanogaster and other tephritid species (the Mediterranean Fruit Fly, C. capitata; the Northern walnut husk fly, Rhagoletis suavis; and the apple maggot Rhagoletis pomonella) [41,42]. The results are presented as a Maximum Likelihood mid-point rooted tree in Fig 5. As expected, the BdorOBPs clustered together with orthologous OBPs from Drosophila and other tephritids with the best BLASTP hit. The classic OBPs from B. dorsalis shared phylogenetic relationships with homologs OBP from Drosophila and other tephritid species, which were previously classified as PBPRPs [44,64]. The nine B. dorsalis PBPRPs (BdorOBP69a, BdorOBP83a-1, BdorOBP83a-2, Bdor19d-1, Bdor19d-2, Bdor19d-3, BdorOBP28a, BdorOBP84a-1 and BdorOBP84a-2) were distributed in four well distinct clades together with homologous genes from the tephritid species. Bdor19d-1, Bdor19d-2, Bdor19d-3 and BdorOBP28a grouped with DmelOBP19d/PBPRP2 and DmelOBP28a/PBPRP5 because of the high degree of similarity between DmelOBP19d/PBPRP2 and DmelOBP28a/PBPRP5 [64]. In both cases, the four B. dorsalis Minus-C OBPs (BdorOBP99c-1, BdorOBP99c-2, BdorOBP99c-3 and BdorOBP99c-4) clustered with the D. melanogaster and other tephritid Minus-C ortholog clade. BdorOBPlush and BdorOBP19a, which had robust expression levels in the antennae, also clustered together with homologous OBPs.
Further, the phylogenetic relationship between four BdorCSPs identified in this study, and previously reported 5 GmmCSPs [14] and 47 D. melanogaster CSPs [41] are shown in the Maximum Likelihood mid-point rooted tree in Fig 6. The four BdorCSPs (BdorCSP1-4) were distributed separately into four well distinct clades together with tephritid homologs. Interestingly, the antenna-specific BdorCSP3 appeared to be closely related to GmmCSP2, which was previously reported to be involved in olfaction [14]. Moreover, the two antenna-specific CSPs BdorCSP3 and GmmCSP2 were orthologs of DmelCSP3 (DmelA10 or OS-D), which is expressed in the antennal segment 3 of D. melanogaster sensillum coeloconicum [65,66]. BdorCSP1, GmmCSP4 and 12 Drosophila CSPs diverged from the rest of the genes.
Sequence similarity analysis of the B. dorsalis ORs revealed that BdorORCO grouped into a conserved clade containing olfactory co-receptors from A. gambiae and D. melanogaster (Fig 7). The other BdorORs clustered together with the Drosophila ORs that produced the best BLASTX hits. In addition, BdorGR21 and BdorGR63 were found in a clade with two carbon dioxide receptors from A. gambiae [67,68] and D. melanogaster [69,70], respectively (Fig 7).
Phylogenetic analysis of all IRs generated an ML tree (Fig 8), which showed clustering of BdorIRs with ‘‘divergent IR” or “Non-NMDA iGluRs” in the clade. A limitation here is the possible identification of false positives due to low expression of some transcripts in the whole insect. Such low expression may have resulted in the lack of identification of two conserved IRs (IR25 and IR8a) typically found in all insects.
The BdorSNMPs grouped together with orthologs from Drosophila sp. (D. melanogaster and D. pseudoobscura) (Fig 9), BdorSNMP1-1 and BdorSNMP1-2 were found in a clade with Drosophila SNMP1, and BdorSNMP2 was found in a clade with SNMP2 from Drosophila sp.
Discussion
The B. dorsalis transcriptomes reported thus far have focused on genes related to development [71,72], digestion, detoxification[73,74,75], sexual dimorphism and reproduction [76]. However, the B. dorsalis chemosensory genes have not been characterized previously at the transcriptome level. To identify the chemosensory genes in B. dorsalis, we sequenced the transcriptome of all the developmental stages including adult chemosensory tissues. Our study identified 31 OBPs, 4 CSPs 23 ORs, 11 IRs, 6 GRs and 3 SNMPs. It is noteworthy to mention that in addition to the new OBPs and new CSPs in our study, we also identified chemosensory membrane proteins (23 ORs, 11 IRs, 6 GRs and 3 SNMPs) previously not reported. Interestingly, we did not identify any IR-coreceptors (IR8a and IR25a) likely due to the limited transcriptome coverage in our study (one total Run for RNA-seq) and/or the low abundance of IR-coreceptors in this species. To some degree, the RNA-seq data would only provide limited reference information.
Based on current research, insect OBPs and CSPs have been assigned two different functions: olfaction and non-olfaction [61]. Most olfaction-related OBPs and CSPs are abundant in the sensillum lymph of olfactory organs (antennae and maxillary palp), and play a critical role as solubilizers and/or carriers of odorants and pheromones[63,64,77–82]. However, non-olfaction-related OBPs and CSPs have been found in the pheromone gland secretions involved in the delivery of semiochemicals (example, pheromones)[83–89], and secretions of the reproductive organs involved in egg and embryo development[90,91]. The expression patterns of both olfaction and non-olfaction related B. dorsalis OBP and CSP genes determined in this study could provide insights to the functions of these proteins. Importantly, nine BdorOBPs and BdorCSP3 showed the highest expression in antenna, suggesting an olfactory role for these genes with antenna being the major olfactory organ. The PBRP homologs, BdorOBP19d-1, BdorOBP19d-2 and BdorOBP28, were also highly expressed in the antennae and in non-chemosensory tissues (Fig 2) indicating that these OBPs may also have other non-olfactory functions. qRT-PCR analyses of antennae-rich or antennae-specific genes in adult males and females revealed that BdorOBPlush, BdorOBP83a-2 and BdorOBP84a-1 were expressed higher in the male antennae than in female antennae (Fig 3). These genes are more likely to play a role in the odorant perception of sex pheromones or as male attractants.
Recently, it was shown that the GOBP protein (100% amino acid sequence similarity with BdorOBP84a-1) extracted from gravid B. dorsalis female antennae had a high affinity to the male attractant, ME, which is used widely as a lure along with insecticides to attract insects for pest control. However, ME is a powerful attractant for B. dorsalis males rather than females. In contrast, BdorOBP19a was expressed higher in the female antennae (Fig 3) suggesting that it could play a role in odorant perception of sex pheromone or oviposition behavior. This OBP can be a potential target to attract female B. dorsalis flies.
In general, flies express four CSPs and the D. melanogaster CSP3 homolog has been demonstrated to play a role in the detection of odorants in Glossina morsitans morsitans [14]. We found that BdorCSP3 had antennae-specific expression profile, which may be critical for the perception of some host volatiles as reported previously through binding assays and RNAi coupled with electrophysiological tests [10]. In addition, BdorOBP56e was expressed only in specific developmental stages (eggs to 1d-pupae) (S1 Fig) suggesting its involvement in egg and pupal development. Interestingly, BdorOBP19c was specifically expressed in the mature female pheromone gland rather than immature pheromone gland indicating a possible involvement of this OBP in the binding and transportation of female specific sex pheromones and their precursors [61].
Compared to OBP and CSP transcripts, ORs are highly restricted to the antennae and are expressed in low levels. Consistently, we found that all the B. dorsalis ORs were also specifically expressed in the antennae of both sexes. However, a new member of the chemosensory receptor family, Ionotropic receptor, did not follow this expression pattern but was expressed in other tissues. Among 11 BdorIRs we discovered, 6 IRs (BdorIR40a, BdorIR41a, BdorIR75d, BdorIR76b, BdorIR84s and BdorIR92a) were specifically expressed in the antennae of both sexes suggesting that they could be involved in odorant detection.
Moreover, phylogenetic analysis indicated clustering of BdorIR40a with A. gambiae and D. melanogaster IR40a, which could directly detect DEET and is a target of insect repellents [92]. Thus, BdorIR40a could be a good target for pest control. In addition, the extent of this shared “DEET repellency” could account for the low degree of species-specific diversity of IRs among the dipterans analyzed here. In addition, two GRs (BdorGR63 and BdorGR21) homologous to insect CO2 receptors [67–70] were specifically expressed in the antennae of both sexes suggesting their involvement in the detection of CO2. Interestingly, two SNMP transcripts, BdorSNMP1-1 and BdorSNMP1-2, displayed a high antennae biased expression profile while BdorSNMP2 was highly expressed in the legs. It is plausible that the two SNMP1s are important for chemosensory function, while SNMP2 may have functions in addition to chemosensation.
Conclusions
By sequencing the transcriptome from various B. dorsalis developmental stages, we identified a variety of genes potentially involved in olfactory signal detection and pheromone biosynthesis in this notorious fruit fly pest. Expression profile analysis revealed that 9 OBPs, 1 CSPs, 23 ORs, 2 SNMPs, 6 IRs and 2 GRs are specifically or mainly expressed in the male and female antennae. The antennae-enriched OBPs, CSPs, ORs, IRs and SNMPs could play a role in the detection of pheromones and general odorants. The identified OBP (BdorOBP19c) could play a role in the binding and transportation of female specific sex pheromones and their precursors. The chemosensory genes identified in our study will provide the basis for functional studies.
Supporting Information
S2 Table. BLASTX analyses of B. dorsalis chemosensory genes compared to D. melanogaster peptide database and suggested chemosensory gene names.
https://doi.org/10.1371/journal.pone.0129794.s002
(XLSX)
S3 Table. Primers used in RT-PCR and qRT-PCR.
https://doi.org/10.1371/journal.pone.0129794.s003
(XLSX)
S4 Table. GO analyses of the transcriptome data from different B. dorsali developmental stages.
https://doi.org/10.1371/journal.pone.0129794.s004
(XLSX)
S5 Table. List of B. dorsalis chemosensory genes putatively involved in odorant binding.
https://doi.org/10.1371/journal.pone.0129794.s005
(XLSX)
S6 Table. List of B. dorsalis chemosensory genes putatively involved in chemosensory reception.
https://doi.org/10.1371/journal.pone.0129794.s006
(XLSX)
S7 Table. The RPKM value of candidate chemosensory genes in different B. dorsalis development stages.
https://doi.org/10.1371/journal.pone.0129794.s007
(XLSX)
S1 Fig. Developmental stage-specific expression of BdorOBP19c,BdorOBP56e and BdorOBP83cd.
https://doi.org/10.1371/journal.pone.0129794.s008
(TIF)
S2 Text. Chemosensory sequences identified in B. dorsalis in Fasta format.
https://doi.org/10.1371/journal.pone.0129794.s010
(DOCX)
Acknowledgments
We thank Mengqiu Qu for insect rearing, Mei Li and Renzhao Xu for insect dissection, Wanyu Xiao for dendrograms drawing, Lei Chen for bioinformatics support.
Author Contributions
Conceived and designed the experiments: JL Z. Wu. Performed the experiments: Z. Wu HZ Z. Wang. Analyzed the data: Z. Wu. Contributed reagents/materials/analysis tools: SB HH. Wrote the paper: Z. Wu JL.
References
- 1. Matsuo T, Sugaya S, Yasukawa J, Aigaki T, Fuyama Y (2007) Odorant-binding proteins OBP57d and OBP57e affect taste perception and host-plant preference in Drosophila sechellia. PLoS Biol; 5: 985–996.
- 2. Asahina K, Pavlenkovich V, Vosshall LB (2008) The survival advantage of olfaction in a competitive environment. Curr Biol; 18: 1153–1155. pmid:18674910
- 3. Whiteman NK, Pierce NE (2008) Delicious poison: genetics of Drosophila host plant preference. Trends Ecol Evol; 23: 473–478. pmid:18657878
- 4. Montell C (2009) A taste of the Drosophila gustatory receptors. Curr Opin Neurobiol; 19: 345–353. pmid:19660932
- 5. Hallem EA, Dahanukar A, Carlson JR (2006) Insect odor and taste receptors. Annu Rev Entomol; 51: 113–135. pmid:16332206
- 6. Leal WS (2013) Odorant Reception in Insects: Roles of Receptors, Binding Proteins, and Degrading Enzymes. Annu Rev Entomol; 58: 373–391. pmid:23020622
- 7. Robertson HM, Warr CG, Carlson JR (2003) Molecular evolution of the insect chemoreceptor gene superfamily in Drosophila melanogaster. Proc Natl Acad Sci SA; 1002: 14537–14542.
- 8. Scott K, Brady RJ, Cravchik A, Morozov P, Rzhetsky A, Zuker C, et al. (2001) A chemosensory gene family encoding candidate gustatory and olfactory receptors in Drosophila. Cell; 104: 661–673. pmid:11257221
- 9. Zhang NY, Ye ZF, Yang K, Dong SL (2014) Antenna-predominant and male-biased CSP19 of Sesamia inferensis able to bind the female sex pheromones and host plant volatiles. Gene; 536: 279–286. pmid:24361960
- 10. Yi X, Wang P, Wang Z, Cai J, Hu M, Zhong G. (2014) Involvement of a specific chemosensory protein from Bactrocera dorsalis in perceiving host plant volatiles. J Chem Ecol; 40: 267–275. pmid:24627092
- 11. Yi X, Zhao H, Dong X, Wang P, Hu M, Zhong G. (2013) BdorCSP2 Is Important for Antifeed and Oviposition-Deterring Activities Induced by Rhodojaponin-III against Bactrocera dorsalis. PloS One; 8: e77295. pmid:24155937
- 12. Zhang T, Wang W, Zhang Z, Zhang Y, Guo Y (2013) Functional Characteristics of a Novel Chemosensory Protein in the Cotton Bollworm Helicoverpa armigera (Hubner). J Integr Agr; 12: 853–861.
- 13. Gu SH, Wang SY, Zhang XY, Ji P, Liu JT, Wang GR, et al. (2012) Functional Characterizations of Chemosensory Proteins of the Alfalfa Plant Bug Adelphocoris lineolatus Indicate Their Involvement in Host Recognition. PLoS One; 7: e42871. pmid:22900060
- 14. Liu R, He X, Lehane S, Lehane M, Hertz-Fowler C, Berriman M, et al. (2012) Expression of chemosensory proteins in the tsetse fly Glossina morsitans morsitans is related to female host-seeking behaviour. Insect Mol Biol; 21: 41–48. pmid:22074189
- 15. Allwood AJ, Chinajariyawong A, Kritsaneepaiboon S, Drew RAI, Hamacek EL, Hancock DL, et al. (1999) Host plant records for fruit flies (Diptera: Tephritidae) in south east Asia. Raffles Bulletin of Zoology Supplement; 7:1–92.
- 16.
Schutze MK, Aketarawong N, Amornsak W, Armstrong KF, Augustinos AA, Barr N, et al. (2014) Synonymization of key pest species within the Bactrocera dorsalis species complex (Diptera: Tephritidae): taxonomic changes based on a review of 20 years of integrative morphological, molecular, cytogenetic, behavioural and chemoecological data. Systematic Entomology; 1–16.
- 17. Wan X, Liu Y, Zhang B (2012) Invasion history of the oriental fruit fly, Bactrocera dorsalis, in the Pacific-Asia region: two main invasion routes. PLoS One; 7: e36176. pmid:22567138
- 18. Wu ZZ, Li HM, Bin SY, Ma J, He HL, Li XF, et al. (2014) Sequence analysis of mitochondrial ND1 gene can reveal the genetic structure and origin of Bactrocera dorsalis s.s. BMC Evol Biol; 14: 55. pmid:24655832
- 19. Kumaran N, Balagawi S, Schutze MK, Clarke AR (2013) Evolution of lure response in tephritid fruit flies: phytochemicals as drivers of sexual selection. Anim Behav; 85: 781–789.
- 20. Pagadala Damodaram KJ, Kempraj V, Aurade RM, Venkataramanappa RK, Nandagopal B, Verghese A, et al. (2014) Oviposition site-selection by Bactrocera dorsalis is mediated through an innate recognition template tuned to γ-octalactone. PloS One; 9: e85764. pmid:24465690
- 21. Kamala Jayanthi PD, Kempraj V, Aurade RM, Venkataramanappa RK, Nandagopal B, Verghese A, et al. (2014) Specific volatile compounds from mango elicit oviposition in gravid Bactrocera dorsalis females. J Chem Ecol; 40 (3): 259–66. pmid:24623046
- 22. Tan KH, Nishida R (2012) Methyl eugenol: Its occurrence, distribution, and role in nature, especially in relation to insect behavior and pollination. J Insect Sci; 12: 1–74. pmid:23465075
- 23.
Vargas RI, Shelly TE, Leblanc L, Pinero JC (2010) Recent Advances in Methyl Eugenol and Cue-Lure Technologies for Fruit Fly Detection, Monitoring, and Control in Hawaii. Vitamins and Hormones; Pheromones: 575–595.
- 24. Zheng W, Peng W, Zhu C, Zhang Q, Saccone G, Zhang H. (2013) Identification and Expression Profile Analysis of Odorant Binding Proteins in the Oriental Fruit Fly Bactrocera dorsalis. Int J Mol Sci; 14: 14936–14949. pmid:23867609
- 25. Zheng W, Zhu C, Peng T, Zhang H (2012) Odorant receptor co-receptor Orco is upregulated by methyl eugenol in male Bactrocera dorsalis (Diptera: Tephritidae). J Insect Physiol; 58: 1122–1127. pmid:22634470
- 26. Jayanthi KP, Kempraj V, Aurade RM, Roy TK, Shivashankara KS, Verghese A. (2014) Computational reverse chemical ecology: virtual screening and predicting behaviorally active semiochemicals for Bactrocera dorsalis. BMC Genomics; 15: 209. pmid:24640964
- 27. Li X, Zhang M, Zhang H (2011) RNA interference of four genes in adult Bactrocera dorsalis by feeding their dsRNAs. PLoS One; 6: e17788. pmid:21445257
- 28. Jayanthi PDK, Verghese A (2002) A simple and cost-effective mass rearing technique for the tephritid fruit fly, Bactrocera dorsalis (Hendel). Current Science; 82(3): 266–268.
- 29. Chen Y, Lin C, Wang C, Wu H, Hwang P (2007) An optimized procedure greatly improves EST vector contamination removal. BMC Genomics; 8: 416. pmid:17997864
- 30. Miller JR, Koren S, Sutton G (2010) Assembly algorithms for next-generation sequencing data. Genomics; 95: 315–327. pmid:20211242
- 31. Chou H, Holmes MH (2001) DNA sequence quality trimming and vector removal. Bioinformatics; 17: 1093–1104. pmid:11751217
- 32. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, et al. (1997) Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res; 25: 3389–3402. pmid:9254694
- 33. Götz S, García-Gómez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, et al. (2008) High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res; 36: 3420–3435. pmid:18445632
- 34. Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M, et al. (2005) Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics; 21: 3674–3676. pmid:16081474
- 35. Xu YL, He P, Zhang L, Fang SQ, Dong SL, Zhang YJ, et al. 2009. Large-scale identification of odorant-binding proteins and chemosensory proteins from expressed sequence tags in insects. BMC genomics,10: 632 pmid:20034407
- 36. Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, et al. (2014) Pfam: the protein families database. Nucleic Acids Res; 42: 222–230.
- 37. Petersen TN, Brunak S, von Heijne G, Nielsen H (2011) SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods; 8: 785–786. pmid:21959131
- 38. Krogh A, Larsson B, von Heijne G, Sonnhammer EL (2001) Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol; 305: 567–580. pmid:11152613
- 39. Finn RD, Mistry J, Tate J, Coggill P, Heger A, Pollington JE, et al. (2010) The Pfam protein families database. Nucleic Acids Res; 38: 211–222.
- 40. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B (2008) Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods; 5: 621–628. pmid:18516045
- 41. Vieira FG, Rozas J (2011) Comparative Genomics of the Odorant-Binding and Chemosensory Protein Gene Families across the Arthropoda: Origin and Evolutionary History of the Chemosensory System. Genome Biol Evol; 3: 476–490. pmid:21527792
- 42. Vieira FG, Sanchez-Gracia A, Rozas J (2007) Comparative genomic analysis of the odorant-binding protein family in 12 Drosophila genomes: purifying selection and birth-and-death evolution. Genome Biol; 8: R235. pmid:18039354
- 43. Hekmat-Scafe DS, Scafe CR, McKinney AJ, Tanouye MA (2002) Genome-wide analysis of the odorant-binding protein gene family in Drosophila melanogaster. Genome Res; 12: 1357–1369. pmid:12213773
- 44. Siciliano P, Scolari F, Gomulski LM, Falchetto M, Manni M, Gabrieli P, et al. (2014) Sniffing Out Chemosensory Genes from the Mediterranean Fruit Fly, Ceratitis capitata. PLoS One; 9: e85523. pmid:24416419
- 45. Schwarz D, Robertson HM, Feder JL, Varala K, Hudson ME, Ragland GJ, et al. (2009) Sympatric ecological speciation meets pyrosequencing: sampling the transcriptome of the apple maggot Rhagoletis pomonella. BMC Genomics; 10: 633. pmid:20035631
- 46. Ramsdell KM, Lyons-Sobaski SA, Robertson HM, Walden KK, Feder JL, Wanner K, et al. (2010) Expressed sequence tags from cephalic chemosensory organs of the northern walnut husk fly, Rhagoletis suavis, including a putative canonical odorant receptor. J Insect Sci; 10: 51. pmid:20569128
- 47. Robertson HM, Warr CG, Carlson JR (2003) Molecular evolution of the insect chemoreceptor gene superfamily in Drosophila melanogaster. Proc Natl Acad Sci U S A; 100: 14537–14542. pmid:14608037
- 48. Gao Q, Chess A (1999) Identification of candidate Drosophila olfactory receptors from genomic DNA sequence. Genomics; 60: 31–39. pmid:10458908
- 49. Hill CA, Fox AN, Pitts RJ, Kent LB, Tan PL, Chrystal MA, et al. (2002) G protein-coupled receptors in Anopheles gambiae. Science; 298: 176–178. pmid:12364795
- 50. Croset V, Cummins SF, Benton R (2010) Ancient Protostome Origin of Chemosensory Ionotropic Glutamate Receptors and the Evolution of Insect Taste and Olfaction. J Neurogenet; 24: 30–31.
- 51. Vogt RG, Miller NE, Litvack R, Fandino RA, Sparks J, Staples J, et al. (2009) The insect SNMP gene family. Insect Biochem Molec; 39: 448–456. pmid:19364529
- 52. Katoh K, Toh H (2010) Parallelization of the MAFFT multiple sequence alignment program. Bioinformatics; 26: 1899–1900. pmid:20427515
- 53. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S (2013) MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol; 30: 2725–2729. pmid:24132122
- 54. Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, et al. (2012) Primer3-new capabilities and interfaces. Nucleic Acids Res; 40: e115. pmid:22730293
- 55. Shen G, Jiang H, Wang X, Wang J (2010) Evaluation of endogenous references for gene expression profiling in different tissues of the oriental fruit fly Bactrocera dorsalis (Diptera: Tephritidae). BMC Mol Biol; 11: 76. pmid:20923571
- 56. Findlay GD, Yi X, Maccoss MJ, Swanson WJ (2008) Proteomics reveals novel Drosophila seminal fluid proteins transferred at mating. PLoS Biol; 6: e178. pmid:18666829
- 57. Pelosi P, Zhou JJ, Ban LP, Calvello M (2006) Soluble proteins in insect chemical communication. Cell Mol Life Sci; 63: 1658–1676. pmid:16786224
- 58. Pophof B (2004) Pheromone-binding proteins contribute to the activation of olfactory receptor neurons in the silkmoths Antheraea polyphemus and Bombyx mori. Chem Senses; 29: 117–125. pmid:14977808
- 59. Graham LA, Brewer D, Lajoie G, Davies PL (2003) Characterization of a subfamily of beetle odorant-binding proteins found in hemolymph. Mol Cell Proteomics; 2: 541–549. pmid:12883044
- 60. Kaissling KE (2001) Olfactory perireceptor and receptor events in moths: a kinetic model. Chem Senses; 26: 125–150. pmid:11238244
- 61. Pelosi P, Iovinella I, Felicioli A, Dani FR (2014) Soluble proteins of chemical communication: an overview across arthropods. Front Physiol; 5: 320. pmid:25221516
- 62. Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods; 25: 402–408. pmid:11846609
- 63. Pelosi P, Maida R (1995) Odorant-binding proteins in insects. Comp Biochem Phys B; 111: 503–514. pmid:7613772
- 64. Pikielny CW, Hasan G, Rouyer F, Rosbash M (1994) Members of a family of Drosophila putative odorant-binding proteins are expressed in different subsets of olfactory hairs. Neuron; 12: 35–49. pmid:7545907
- 65. McKenna MP, Hekmat-Scafe DS, Gaines P, Carlson JR (1994) Putative Drosophila pheromone-binding proteins expressed in a subregion of the olfactory system. J Biol Chem; 269: 16340–16347. pmid:8206941
- 66. Pikielny CW, Hasan G, Rouyer F, Rosbash M (1994) Members of a family of Drosophila putative odorant-binding proteins are expressed in different subsets of olfactory hairs. Neuron; 12: 35–49. pmid:7545907
- 67. Turner SL, Ray A (2009) Modification of CO2 avoidance behaviour in Drosophila by inhibitory odorants. Nature; 461: 159–277.
- 68. Tauxe GM, MacWilliam D, Boyle SM, Guda T, Ray A (2013) Targeting a dual detector of skin and CO2 to modify mosquito host seeking. Cell; 155: 1365–1379. pmid:24315103
- 69. Jones WD, Cayirlioglu P, Kadow IG, Vosshall LB (2007) Two chemosensory receptors together mediate carbon dioxide detection in Drosophila. Nature; 445: 86–90. pmid:17167414
- 70. Kwon JY, Dahanukar A, Weiss LA, Carlson JR (2007) The molecular basis of CO2 reception in Drosophila. Proc Natl Acad Sci U S A; 104: 3574–3578. pmid:17360684
- 71. Zheng W, Peng T, He W, Zhang H (2012) High-throughput sequencing to reveal genes involved in reproduction and development in Bactrocera dorsalis (Diptera: Tephritidae). PLoS One; 7: e36463. pmid:22570719
- 72. Shen GM, Dou W, Niu JZ, Jiang HB, Yang WJ, Jia FX, et al. (2011) Transcriptome analysis of the oriental fruit fly (Bactrocera dorsalis). PLoS One; 6: e29127. pmid:22195006
- 73. Yang WJ, Yuan GR, Cong L, Xie YF, Wang JJ (2014) De novo cloning and annotation of genes associated with immunity, detoxification and energy metabolism from the fat body of the oriental fruit fly, Bactrocera dorsalis. PLoS One; 9: e94470. pmid:24710118
- 74. Hsu JC, Chien TY, Hu CC, Chen MJ, Wu WJ, Feng HT, et al. (2012) Discovery of genes related to insecticide resistance in Bactrocera dorsalis by functional genomic analysis of a de novo assembled transcriptome. PLoS One; 7: e40950. pmid:22879883
- 75. Shen GM, Dou W, Huang Y, Jiang XZ, Smagghe G, Wang JJ. (2013) In silico cloning and annotation of genes involved in the digestion, detoxification and RNA interference mechanism in the midgut of Bactrocera dorsalis [Hendel (Diptera: Tephritidae)]. Insect Mol Biol; 22: 354–365. pmid:23577657
- 76. Geib SM, Calla B, Hall B, Hou S, Manoukis NC (2014) Characterizing the developmental transcriptome of the oriental fruit fly, Bactrocera dorsalis (Diptera: Tephritidae) through comparative genomic analysis with Drosophila melanogaster utilizing modENCODE datasets. BMC Genomics; 15: 942. pmid:25348373
- 77. Zhou JJ (2010) Odorant-Binding Proteins in Insects. Vitamins and Hormones; 83: 241–272. pmid:20831949
- 78. Grosse-Wilde E, Svatos A, Krieger J (2006) A pheromone-binding protein mediates the bombykol-induced activation of a pheromone receptor in vitro. Chem Senses; 31: 547–555. pmid:16679489
- 79. Pelosi P, Zhou JJ, Ban LP, Calvello M (2006) Soluble proteins in insect chemical communication. Cell Mol Life Sci; 63: 1658–1676. pmid:16786224
- 80. Xu PX, Atkinson R, Jones D, Smith DP (2005) Drosophila OBP LUSH is required for activity of pheromone-sensitive neurons. Neuron; 45: 193–200. pmid:15664171
- 81. Nakagawa T, Sakurai T, Nishioka T, Touhara K (2005) Insect sex-pheromone signals mediated by specific combinations of olfactory receptors. Science; 307: 1638–1642. pmid:15692016
- 82. Vogt RG, Riddiford LM (1981) Pheromone binding and inactivation by moth antennae. Nature; 293: 161–163. pmid:18074618
- 83. Zhou X, Ban L, Iovinella I, Zhao L, Gao Q, Felicioli A, et al. (2013) Diversity, abundance, and sex-specific expression of chemosensory proteins in the reproductive organs of the locust Locusta migratoria manilensis. Biol Chem; 394: 43–54. pmid:23096575
- 84. Gu SH, Wu KM, Guo YY, Pickett JA, Field LM, Zhou JJ, et al. (2013) Identification of genes expressed in the sex pheromone gland of the black cutworm Agrotis ipsilon with putative roles in sex pheromone biosynthesis and transport. BMC Genomics; 14: 636. pmid:24053512
- 85. Dani FR, Michelucci E, Francese S, Mastrobuoni G, Cappellozza S, Marca GL, et al. (2011) Odorant-binding proteins and chemosensory proteins in pheromone detection and release in the silkmoth Bombyx mori. Chem Senses; 36: 335–344. pmid:21220518
- 86. Jacquin-Joly E, Vogt RG, Francois MC, Nagnan-Le MP (2001) Functional and expression pattern analysis of chemosensory proteins expressed in antennae and pheromonal gland of Mamestra brassicae. Chem Senses; 26: 833–844. pmid:11555479
- 87. Sun YL, Huang LQ, Pelosi P, Wang CZ (2012) Expression in antennae and reproductive organs suggests a dual role of an odorant-binding protein in two sibling Helicoverpa species. PLoS One; 7: e30040. pmid:22291900
- 88. Iovinella I, Dani FR, Niccolini A, Sagona S, Michelucci E, Gazzano A, et al. (2011) Differential expression of odorant-binding proteins in the mandibular glands of the honey bee according to caste and age. J Proteome Res; 10: 3439–3449. pmid:21707107
- 89. Li S, Picimbon JF, Ji S, Kan Y, Chuanling Q, Zhou JJ, et al. (2008) Multiple functions of an odorant-binding protein in the mosquito Aedes aegypti. Biochem Biophys Res Commun; 372: 464–468. pmid:18502197
- 90. Marinotti O, Ngo T, Kojin BB, Chou SP, Nguyen B, Juhn J, et al. (2014) Integrated proteomic and transcriptomic analysis of the Aedes aegypti eggshell. BMC Dev Biol; 14: 15. pmid:24707823
- 91. Maleszka J, Foret S, Saint R, 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. pmid:17216269
- 92. Kain P, Boyle SM, Tharadra SK, Guda T, Pham C, Dahanukar A, et al. (2013) Odour receptors and neurons for DEET and new insect repellents. Nature; 502: 507–512. pmid:24089210