Article | Open | Published:

# BdorOBP2 plays an indispensable role in the perception of methyl eugenol by mature males of Bactrocera dorsalis (Hendel)

## Abstract

Bactrocera dorsalis (Hendel) is a fruit-eating pest that causes substantial economic damage to the fresh produce industry in tropical and sub-tropical countries. Methyl eugenol (ME) is a powerful attractant for mature males of B. dorsalis, and has been widely used for detecting, luring and eradicating B. dorsalis populations worldwide. However, the molecular mechanism underlying the olfactory perception of ME remains largely unknown. Here, we analyzed the differential proteomics profiling of the antennae between ME-responsive and ME-non-responsive males by using isobaric tags for relative and absolute quantitation (iTRAQ). In total, 4622 proteins were identified, of which 277 proteins were significant differentially expressed, with 192 up-regulated and 85 down-regulated in responsive male antennae. Quantitative real-time PCR (qRT-PCR) analysis confirmed the authenticity and accuracy of the proteomic analysis. Based on the iTRAQ and qRT-PCR results, we found that the odorant-binding protein 2 (BdorOBP2) was abundantly expressed in responsive male antennae. Moreover, BdorOBP2 was significantly up-regulated by ME in male antennae. Mature males showed significantly greater taxis toward ME than did mature females. Silencing BdorOBP2 reduced mature males’ responsiveness to ME. These results indicate that BdorOBP2 may play an essential role in the molecular mechanism underlying B. dorsalis olfactory perception of ME.

## Introduction

The Oriental fruit fly, Bactrocera dorsalis (Hendel), is one of the most destructive fruit/vegetable-eating agricultural pests in the world, particularly in Asiatic countries such as China1,2. B. dorsalis is a typical polyphagous pest with larval stage feeding, multiple mating, long life span and great fecundity of adults, capable of causing severe damage to more than 250 commercially-valuable tropical and subtropical crops, especially some staple fruits, including mango, guava, orange, carambola, jujube, loquat, peach, etc.3,4,5. Considering the mining habits of the larval stage, pest control strategies focus mainly on the adults. Trapping technology based on flavor cues is an important strategy for controlling these fruit flies. The natural phenylpropanoid, methyl eugenol (ME), has been used widely and effectively for monitoring and eradicating males of B. dorsalis populations worldwide6,7,8. ME functions as a precursor for the synthesis of a sex pheromone in B. dorsalis males, which makes those males more attractive to females9. However, there are human health implications. ME is considered carcinogenic to humans by the National Toxicology Program (NTP), US Department of Health and Human Services, making it not suitable for long-term use10,11,12,13,14. In addition, ME can trap sexually mature male flies, but it fails to attract immature males or females. Although, several analogues of ME have been developed13,15,16,17, the trade-off with male attractancy makes them an unlikely improvement over ME for detecting the oriental fruit fly in the field. Therefore, clarifying the molecular mechanism underlying the olfactory perception of ME should provide much-needed information for the developing of efficient, simple, green and sustainable lures for monitoring and controlling B. dorsalis pest populations.

Insects have mastered the art of using semiochemicals, encompassing pheromones, kairomones, repellents and attractants, as communication signals and to process stimuli from the environment, including the whereabouts of food sources, reproductive partners, oviposition sites, hosts, or for detecting predators from a distance1,2,18. Insect olfaction is highly sensitive so as to reliably discriminate among semiochemicals to mediate these important behaviors. Not surprisingly, the olfactory system of insects is sophisticated and involves many kinds of proteins: namely, odorant-binding proteins (OBPs), chemosensory proteins (CSPs), odorant receptors (ORs), ionotropic receptors (IRs), odorant-degrading enzymes (ODEs) and sensory neuron membrane proteins (SNMPs)19,20. The initial steps in odor detection involve the binding of an odor to the ORs displayed on the dendrites of olfactory sensory neurons (OSNs)21,22. Insect OBPs that are located in the aqueous sensillum lymph of the antennae are thought to be the first proteins binding the odors to facilitate the transport of hydrophobic odorants through the aqueous surroundings to peripheral ORs, thereby starting the signal transduction cascade leading to behavioral outputs23,24,25,26,27. Recently, the use of OBPs as a molecular target for rapidly elucidating behaviorally active compounds and screening of novel species-specific repellents or attractants in insects is gaining attention28,29,30. It has been demonstrated that OBP1 plays an essential role in mediating indole recognition in the antennae of female Anopheles gambiae 31. Additionally, knock-down of OBP1 in the mosquito Culex quinquefasciatus reduced its antennal response to several oviposition attractants32. By using a molecular docking and molecular dynamics simulation, Kempraj et al.18 studied 25 semiochemicals’ binding potential to a GOBP of B. dorsalis for predicting and screening behaviorally active compounds. Silencing Orco and OBP83a-2 decreased the ability of B. dorsalis adults to find semiochemical lures in a timely manner14,33. In this regard, this unique research strategy may help chemical ecologists in their prospecting of active semiochemicals for eco-friendly use in pest management practices.

Isobaric tags for relative and absolute quantitation (iTRAQ), a new technique that has become popular in proteomic studies in recent years, can provide extremely reliable quantitative comparisons and identifications of molecular candidates among complex biological samples34,35. To date, this approach has been successfully used in many organisms, such as Bemisia tabaci 36, Apostichopus japonicas 35, Locusta migratoria 37, and Culex pipiens pallens 38, to name a few. Interestingly, experiments have revealed that the proportion of males with non-taxis to ME could be increased persistently via several generations of artificial selection under laboratory conditions11,39. These results suggest that the characteristics of defective olfactory sensory physiology to ME of non-responsive individuals may be passed on to offspring. Therefore, we speculated that the responsiveness to ME of male B. dorsalis flies is controlled by those genes associated with olfactory processing. In the present study, we utilized the iTRAQ approach combined with LC-MS/MS analyses to investigate the differential proteomics profiling of the antennae between responsive and non-responsive B. dorsalis mature males to ME. We identified an odorant-binding protein 2 (OBP2) highly expressed in the antennae of responsive male flies and determined its molecular characterization, expression profile, and function in ME perception.

## Results

### Non-responders to ME in the select and control lines

The frequency of non-responsive males to ME in the select and control lines for each generation are illustrated in Fig. 1. In the select lines, the mean proportions of non-responders in the first generation (F1) and second generation (F2) were respectively 14.0 ± 0.9% and 28.8 ± 1.5%, both significantly higher than those in the control lines, which were 6.0 ± 0.8% (t = 13.86; df = 2; P = 0.0052) and 5.3 ± 0.4% (t = 14.62; df = 2; P = 0.0046), respectively. These results indicated that the responsiveness of males to ME evidently decreased via the selection experiments. Further, the proportion of non-responders increased rapidly and remained consistently high (28.0% – 29.3%) for the select males between the second and fifth generations. By contrast, the proportion of non-responders in the control line were not statistically different among all generations (F = 0.24; df = 5; P = 0.9375).

### Identification of differentially expressed antennal proteins between the responsive and non-responsive male flies

The total mass spectra number detected in the male antennae proteomes was 266460, representing 48059 peptide spectra and 26679 distinct peptides. The results also indicated that a total of 4622 proteins were successfully identified (Table S1). The global expression changes of all proteins with quantitative iTRAQ ratios are shown in Fig. 2B; the red circles are proteins that exhibited a differential expression pattern in responsive males compared with the non-responsive males. A heat map (Fig. 2C) was constructed according to the data of 277 differential proteins between responder and non-responder. Compared with the non-responsive males, 192 proteins showed an increased abundance, and 85 proteins showed a decreased abundance in the responsive males antennae (with a fold change of ≥1.20 or ≤ 0.83, P value < 0.05; FDR < 0.01) (Fig. 2A). Among the differentially expressed proteins, we identified five proteins associated with insect olfactory recognition, including four odorant binding proteins, and one odorant receptor (Table 1). The respective mass spectra and peptide sequences of these five olfactory proteins are shown in Figure S1. Notably, BdorOBP2 was up-regulated by 1.34-fold in the responsive males’ relative to the non-responsive males’ antennae. Further, proteins related to energy regulation, protein transportation and binding were also differentially expressed between the responders and non-responders (Supplementary Dataset S1).

### GO annotation of the differentially expressed proteins

For insights into the functional categories that were altered between the responsive and non-responsive male antennae, the Gene Ontology (GO) database and UniProt knowledge base were used to categorize the identified differential proteins. They covered a wide range of biological processes, molecular functions and cellular components, which could be classified into 14, 5, and 6 subcategory groups, respectively (Fig. 3). Specifically, the largest group within the biological process category was that of metabolic processes, followed by cellular processes and single-organism processes, whereas catalytic activity and binding were the most common categories for molecular function. The cellular component functions of these proteins were mainly related to the cell, organelles, membrane and extracellular region. Upon further investigation by the GO Enrichment analysis (P value < 0.05), all of the following were significantly affected: aminoglycan metabolic process, cuticle pigmentation, receptor binding, small molecule binding, neuropeptide hormone activity, protein heterodimerization activity, hydrolase activity acting on glycosyl bonds, growth factor receptor binding, and extracellular region (Figure S2).

### KEGG pathway analysis

To identify the active biological pathways in the ME perception process, we mapped the detected proteins to reference canonical pathways in the Kyoto Encyclopedia of Genes and Genomes (KEGG) by using KAAS (KEGG Automatic Annotation Server). We obtained 174 maps using the proteins, and the predictions for the most differentially expressed proteins suggest that they are involved in 19 pathways (Figure S3). The pathways with the greatest representation by unique proteins were those of purine metabolism, olfactory transporter and basal transcription factors, with six, six and five proteins respectively. Compared with non-responsive males, six proteins with either a greater (BdorOBP2, nucleotide phosphodiesterase) or lower (BdorOBP44a, BdorOBP69a, BdorOBA5, BdorOR94b) abundance in the responsive males’ antennae are related to olfactory transport (Table 1). The annotations provide a valuable resource for investigating the specific processes, functions and pathways involved in the ME perception process.

### Validation of differentially expressed proteins by qRT-PCR

To validate the results from the proteomic analysis, genes encoding differentially expressed levels of olfactory proteins in the antennae between the responsive and non-responsive males were detected by qRT-PCR (Fig. 4). Results showed that the expression of BdorOBP2 in the antennae of responsive males was remarkably higher (t = 97.87; df = 2; P = 0.0001). The transcript levels of BdorOBP44a, BdorOBA5, BdorOBP69a, BdorOR94b, BdorOB19A and BdorPBP4 were respectively 3.36, 2.76, 5.87, 7.17, 2.68 and 2.88 times higher in the non-responsive than in responsive male antennae. For the other genes tested, BdorOB99A, BdorOBP15, BdorOBP83a-1, BdorPBP2 and BdorIR84a were expressed at similar levels in the antennae of the non-responsive and responsive males. The mRNA results by qRT-PCR showed that 10 genes: BdorOBP2, BdorOBP44a, BdorOBA5, BdorOBP69a, BdorOR94b, BdorOB99A, BdorOBP15, BdorOBP83a-1, BdorPBP2 and BdorIR84a were consistent with the protein expression analyzed by iTRAQ. However, two proteins, BdorOB19A and BdorPBP4 expressed no significant difference at the proteomic level, which pointed to mRNA-protein expression inconsistencies.

### Phylogenetic analysis of BdorOBP2

The full-length BdorOBP2 complementary DNA (cDNA) consisted of 621 nucleotides (nt), with an open reading frame (ORF) of 417 bp, and an encoded polypeptide of 139 amino acids. To determine the phylogenetic relationship between BdorOBP2 and the 91 OBPs reported in Bactrocera tau, Bactrocera cucurbitae, Bactrocera latifrons and Ceratitis capitata, a neighbor-joining tree was constructed (Fig. 5). As expected, the BdorOBPs clustered together with orthologous OBPs from the Tephritidae species with the best BLASTP hit. The classic OBPs from B. dorsalis shared phylogenetic relationships with the OBP homologs from the Tephritidae species. Notably, BdorOBP2 is clustered in a branch with the orthologous gene OBP56.

### Sex and age affect the taxis of flies to ME and BdorOBP2 expression

Adult sex and age had a clear influence on the flies’ taxis to ME. As Fig. 6A shows, responsiveness increased with age, with newly emerged adults presenting the lowest taxis. The 15-day-old male responders numbered 180.3 ± 3.3, which was significantly higher than in 3-day-old males (16.7 ± 1.8) and 15-day-old females (28.3 ± 2.0) (F = 951.89; df = 2; P < 0.0001). Furthermore, in contrast to the low expression observed in the 3-day-old male flies, BdorOBP2 was highly expressed in 15-day-old males, being 8.33-fold greater than for 3-day-old males and 2.27-fold that of the 15-day-old females (F = 117.31; df = 2; P = 0.0003) (Fig. 6B).

### Confirmation of ME exposure for up-regulated BdorOBP2 expression

The ME treatment had a clear, positive influence on BdorOBP2 expression. As Fig. 7 shows, this was up-regulated in the mature male antennae to a level 2.72-fold greater than that under the mineral oil treatment after 1 h (t = 7.80; df = 2; P = 0.0161), a difference that reached c. 9.92-fold after 2 h (t = 5.93; df = 2; P = 0.0273).

### BdorOBP2 mediated the taxis of male flies to ME

The bioassay done on B. dorsalis showed that the average mortalities of dsBdorOBP2, dsGFP and the buffer treatment groups at 24 h were 13.6 ± 0.8%, 12.8 ± 1.0% and 15.2 ± 1.0% respectively (Fig. 8A) and the mortality increased to 17.6 ± 0.8%, 16.4 ± 1.5% and 18.0 ± 0.9% after 48 h continued treatment (Fig. 8B). However, the mortalities of the blank controls were comparatively much low, at 3.6 ± 0.8% for 24 h (F = 38.34; df = 3; P < 0.0001) and 4.4 ± 0.8% at 48 h (F = 33.91; df = 3; P < 0.0001). Notably, the mortality of the dsBdorOBP2-treated male flies was not significantly different from that of the dsGFP- and Buffer-treated male flies. The results indicate that the microinjection had a negative impact on B. dorsalis adults’ survival.

The effects of RNAi of the BdorOBP2 precursor transcript on male responsiveness to ME were analyzed. The qRT-PCR demonstrated that the knockdown efficiencies were significant for the BdorOBP2 mRNA levels at 24 h and 48 h after the male flies had been injected with dsBdorOBP2 (Fig. 8C,D). After treatment for 24 h, BdorOBP2 expression decreased approximately 90% compared with those of the dsGFP, Buffer and blank control groups (F = 37.85; df = 3; P < 0.0001). Consequently, the dsBdorOBP2-treated males were much less trapped by ME than were the dsGFP- and Buffer-treated flies. This difference was strongest at 24 h, when the males trapped rate by ME was 12.0 ± 2.0% in dsBdorOBP2 treatment group, significantly less than dsGFP-treated group (88.1 ± 1.4%), Buffer-treated group (89.2 ± 1.2%) and control group (91.7 ± 1.9%) (F = 548.80; df = 3; P < 0.0001) (Fig. 8E). After treatment for 48 h, though the males trapped in the dsBdorOBP2 treatment group had noticeably increased (to 56.8 ± 4.1%), it was still lower, at c. two-thirds that of the other treatment groups and the blank control (F = 23.49; df = 3; P < 0.0001) (Fig. 8F). These results demonstrated that the male flies could not be significantly trapped by ME after BdorOBP2 was silenced.

## Discussion

The Oriental fruit fly is a biologically interesting and economically relevant insect herbivore, whose mature males are effectively trapped by ME. Over the past few decades, numerous studies have been focused on the factors influencing the ME-luring efficiency on males, and the metabolic pathways or metabolites of ME after feeding by B. dorsalis adults7,9,12,40. However, the molecular mechanisms underlying the olfactory perception of ME remain poorly understood. Through the iTRAQ technique, we identified some differentially expressed olfactory proteins of the antennae between responsive and non-responsive B. dorsalis males to ME. Among these olfactory proteins, BdorOBP2 was abundantly expressed in the responsive mature males’ antennae of B. dorsalis and its expression level was remarkably up-regulated by ME. Moreover, silencing BdorOBP2 dramatically decreased the number of attracted males to ME. Taken together, results obtained here indicate that BdorOBP2 likely has an essential role to play in mediating the taxis of ME to B. dorsalis male adults. Hence, our present findings further enrich and deepen our understanding of the molecular mechanism underpinning ME perception by B. dorsalis mature males.

Notably, the responsiveness to ME in B. dorsalis males was not uniform throughout the day: it peaked in the morning, declined in the afternoon, and it then dropped markedly at dusk40. This daily fluctuation in the male responsiveness to ME suggests an inverse relationship with the daily cycle of their mating behaviour. For this reason, we performed all of our behavioural assays in the morning. We found that the proportion of select males non-responsive to ME increases as the generation increases. However, the highest proportion would stabilize at about 28%. Yet more remarkably, the selection experiment did not result in the complete disappearance of the male flies’ tropism to ME. This agrees with previous findings by Shelly11 and Guo et al.39, whose showed that the responsiveness of male B. dorsalis to ME could be reduced via artificial selection (6–12 generations) under laboratory conditions. Nonetheless, the laboratory experiment by Zheng et al.14 demonstrated that a small proportion of males became attracted more slowly, with about 10% of males not even attracted to ME after 24 h. This may be attributed to the individual differences in epigenetics, such as in defective olfactory sensory physiology.

Although the gene and protein databases of B. dorsalis were incomplete, 277 proteins were successfully found as differentially expressed between the responsive and non-responsive males’ antennae by iTRAQ technology. The functions of these differentially expressed proteins are likely involved in olfactory transduction, amino acid metabolism, basal transcription factors, MAPK signaling pathway, PI3K-Akt signaling pathway, endocytosis, phosphatidylinositol signaling system, cGMP-PKG sinaling pathway, calcium signaling pathway and PPAR signaling pathway. The calcium/calmodulin-dependent 3′,5′-cyclic nucleotide phosphodiesterase 1 C, which plays a pivotal role in balancing intracellular Ca2+/CaM and cGMP signaling41, was found abundantly expressed in the antennae of responsive males. It is preliminarily conjectured that this protein may be involved in how B. dorsalis males perceive ME, but this needs further investigation. In addition, four OBPs (BdorOBP2, BdorOBP44a, BdorOBP69a, BdorOBA5) and one OR (BdorOR94b) had significant differential expression between responders and non-responders. But only BdorOBP2 was highly expressed in the antennae of responsive males. The consistency in the results with the mRNA expression by qRT-PCR analysis implies there is a real significance of BdorOBP2 in the ME perception by mature male flies. While a diverse pattern between mRNA and protein expression in two genes (BdorOB19A and BdorPBP4) may reflect the lack of a direct relationship between mRNA and protein expression, it could also be attributed to post-transitional effects and/or other regulatory mechanisms, such as a lack of synchronization42,43,44.

Immature male flies were obviously not attracted to ME, but the highest attraction was observed in 15-day-old males, a result consistent with Karunaratne and Karunaratne40. That BdorOBP2 was abundantly expressed in the mature male antennae agrees with prior findings by Zheng et al.1. Accordingly, very few of the 15-day-old females were attracted by ME. As a corollary, the expression level of BdorOBP2 in the antennae of mature male flies was up-regulated by ME in just 1 h, which strongly suggests that BdorOBP2 upregulation may be closely related to the ME attraction for male adults.

Intriguingly, prior researches have demonstrated that proteins in the GOBP2 class have a conserved sequence across different species, and that they possess high binding activity to major pheromones in several insect species45,46,47. Previous study on Lepidopteran insects has proved that MsexGOBP2 could be photoaffinity-labeled by the pheromone analogue, (6E, 11Z)-hexadecadienyl diazoacetate48. CsupGOBP2 has a high specificity for a major pheromone component, 11Z-hexadecenal49, and BmorGOBP2 can bind to the B. mori sex pheromone component, (10E,12Z)-hexadecadien-1-ol46. The long-chained chemical, trans-11- tetradecen-1-yl acetate, a sex pheromone component of Loxostege sticticalis Pyralidae, has a high binding affinity to LstiGOBP2 47. Silencing of AgOBP2 impaired the host recognition and detection of oviposition attractants in Aphis gossypii Glover. Simultaneously, OBP2 could potentially serve as a practicable target for RNAi-mediated gene silencing in the pest control of Hemipteran insect pests50. In the present study, silencing BdorOBP2 via the injection of dsRNA reduced the attractive effects of ME to mature males. In other words, our results demonstrate for the first time that BdorOBP2 participates in mediating the taxis to ME in male B. dorsalis flies, which sheds light on the mechanisms underlying the strong attraction of its males to ME. Meanwhile, the microinjected dsRNA decreased the survival of B. dorsalis male adults. A plausible explanation for this important phenomenon is that microinjection caused physical damage leading to the death of the tested flies.

Odorant receptor co-receptor (Orco) is an atypical common odorant receptor. It is highly conserved among different insect species, and is co-expressed with a specific odorant receptor in the endomembrane system in olfactory sensory neurons. The Orco/OR complex confers odor-recognition specificity and promotes the functional reconstitution of the odor-evoked signaling transduction pathway in the sensory neurons51,52,53. In this complex, Orco is suggested to be responsible for the stabilization, localization, and for the correct protein folding of OR, thus making it obligatory for the detection of odorants54. In addition, Orco also works as a selective ion channel in the response to common odors, such as aldehydes, ketones, esters, and aromatics54,55. BdOrco was reportedly up-regulated by the sexual attractant ME and played a critical role in mediating the taxis of B. dorsalis males to ME14. Furthermore, gene silencing via the microinjection of dsRNA recently confirmed that BdorOBP83a-2 actively participates in the process of ME detection by B. dorsalis 33. The phylogenetic tree was generated to show how BdorOBP2 and BdorOBP83a-2 belong to different groups of orthologous proteins (Fig. 5). Olfaction in insects is extremely complex, consisting of various interactions of many classes of proteins and effectors to reliably translate any external semiochemicals in the environment, so as to ultimately produce a behavioral response in the insect56. Considering our findings alongside other research to date, we preliminary posit that ME odor molecules bind to either BdorOBP2 or BdorOBP83a-2, and then are transferred to OSNs, where BdOrco coordinates with a specific OR to bind to the ME odor molecules. Only then is the olfactory signal transduction pathway finally activated.

It is the conventional ORs which are activated by odorants that determine the specificity of the response in the olfactory neurons51,56. However, no OR protein that underwent significant up-regulation in the responsive males was identified by iTRAQ in our present study. We speculate that because there is no direct stimulation by ME on males, this did not cause the differential expression of OR proteins. Thus, these OR proteins, which are assumed to operate in the olfactory perception of ME in B. dorsalis, await further investigation. However, because we lack B. dorsalis genomic information, many proteins cannot yet be annotated or their functions and interactions remain poorly understood, especially for those olfactory-related proteins. Nevertheless, with the gradual completion of genome sequencing of B. dorsalis, the molecular mechanism and transport pathway underpinning how ME is able to lure mature male flies will become clearer.

## Materials and Methods

### Ethics statement

No specific permits were required for the described field studies, and no specific permissions were required for these locations/activities. We confirm that these locations are not privately owned or protected in any way and that the field studies did not involve endangered or protected species.

### Insect

The B. dorsalis genetic sexing strain (GSS) had been maintained in a laboratory for approximately thirty generations in South China Agricultural University. The female pupae are white and the male pupae are brown. The flies were reared at 27 ± 1 °C, 75 ± 1% relative humidity, and a photoperiod cycle of 14 h L/10 h D. Hatched larvae were maintained on an artificial diet includes sugar (8.99%), yeast (15.06%), nipagen (0.15%), sodium benzoate (0.15%), wheat germ oil (0.15%), citric acid (1.70%), and water (73.81%)57. Before pupation, larvae were transferred into small plastic boxes with sand. Pupae were kept at 27 ± 1 °C until adults emerged. Adult flies were reared in wooden cages (35 cm by 35 cm by 35 cm) and fed another artificial diet consisting of yeast extract: dry sugar at 1:1 (w/w).

### Artificial selection of non-responsive males to ME

We chose the male flies without taxis to ME according to the methods of Shelly11 and Guo et al.39. The select and control lines were established as follows. At the start of the study, the colony from the B. dorsalis genetic sexing strain had been maintained in our laboratory for approximately 30 generations. Male and female new pupae of the same generation were separated and reared in different screen cages (1.0 m × 1.0 m × 1.0 m). When the adult male flies were 15 days old, a fly trap containing 1.0 mL of pure ME (Energy Chemical Company, Shanghai, China) was put into the screen cage. After treatment for 2 h, we removed the trap and discarded the trapped male flies. The remaining male files were kept and reared with unlimited food and water. Three days after the first trapping trial, a second trial was conducted using the same procedure. The rest of male files in the screen cage were considered non-responsive to ME. For the select line, males that failed to get trapped by ME in the double-exposure experiment were used as sires for the next generation; females at the same emergence age were taken randomly from the same generation. For the control line, males and females were chosen haphazardly from among untested individuals in the control stock. Following the method described above, the flies were selected and bred for five generations.

For the select and control lines, the responsiveness of male progeny in each generation to ME was tested in the laboratory. Firstly, 200 male flies (15 days old) of the select and control lines were randomly chosen as the test subjects. Then the double-exposure experiment method described above was performed to detect the non-responsiveness ratio to ME. Finally, the non-responsiveness to ME ratio of male flies were calculated using the following equation (1). Each experiment consisted of three replicates and was conducted from 9 a.m. to 12 a.m. During experiments, temperature and relative humidity in the room were maintained at 27 ± 1 °C and 75 ± 1% RH, respectively.

$${\rm{N}}{\rm{o}}{\rm{n}}{\rm{ \mbox{-} }}{\rm{r}}{\rm{e}}{\rm{s}}{\rm{p}}{\rm{o}}{\rm{n}}{\rm{s}}{\rm{i}}{\rm{v}}{\rm{e}}{\rm{n}}{\rm{e}}{\rm{s}}{\rm{s}}\,{\rm{t}}{\rm{o}}\,{\rm{M}}{\rm{E}}\,{\rm{r}}{\rm{a}}{\rm{t}}{\rm{e}}\,( \% )=\frac{{\rm{N}}{\rm{o}}.{\rm{o}}{\rm{f}}\,{\rm{m}}{\rm{a}}{\rm{l}}{\rm{e}}\,{\rm{f}}{\rm{l}}{\rm{i}}{\rm{e}}{\rm{s}}\,{\rm{n}}{\rm{o}}{\rm{t}}\,{\rm{t}}{\rm{r}}{\rm{a}}{\rm{p}}{\rm{p}}{\rm{e}}{\rm{d}}\,{\rm{b}}{\rm{y}}\,{\rm{M}}{\rm{E}}\,}{{\rm{N}}{\rm{o}}.{\rm{o}}{\rm{f}}\,{\rm{t}}{\rm{e}}{\rm{s}}{\rm{t}}\,{\rm{m}}{\rm{a}}{\rm{l}}{\rm{e}}\,{\rm{f}}{\rm{i}}{\rm{l}}{\rm{e}}{\rm{s}}}{\rm{\times }}{\rm{100}}$$
(1)

### Identification of antennal differentially expressed proteins between the responsive and non-responsive males by iTRAQ

#### Protein extraction, digestion and iTRAQ labeling

The antennae were dissected separately from the responsive and non-responsive male adults (each was 18 days old; approximately 4500 individuals); these came from the sixth generation colony of the oriental fruit fly in the select line. Their antennae were frozen immediately in liquid nitrogen, and stored at −80 °C until extraction. Three independent biological replicates of the frozen antennae were prepared for the iTRAQ analysis. The antennae were dissected in a SDT lysis buffer (containing 4% sodium dodecyl sulfate [SDS] and 150 mM Tris-HCl [pH 8.0]), and transferred to 2 mL tubes with an amount of quartz sand (another 1/4 inch ceramic bead MP 6540-424 for the tissue samples). The lysate was homogenized by an MP homogenizer (24 × 2, 6.0 M/S, 60 s; done twice). The homogenate was sonicated and then boiled for 15 min. After being centrifuged at 14000 g for 40 min, the supernatant was filtered with 0.22 µm filters. The filtrate was quantified with a BCA Protein Assay Kit (Bio-Rad, USA). For a quality check, 20 µg of the total proteins per sample were mixed with 5 × Loading Buffer and boiled for 5 min. The proteins were separated on a 12.5% SDS-PAGE gel. Protein bands were visualized through Coomassie Blue R-250 staining.

The extracted proteins were digested according to the previously reported method35,58,59. Briefly, 200 μg of proteins for each sample were digested with 4 μg Trypsin Gold (Promega, Madison, WI, USA) in 40 μL DS buffer overnight at 37 °C, and the resulting peptides were collected as a filtrate. The peptides of each sample were desalted on C18 Cartridges (Empore™ SPE Cartridges C18, bed I.D. 7 mm, volume 3 mL, Sigma), concentrated by vacuum centrifugation and reconstituted in 40 µl of 0.1% (v/v) formic acid. The peptide content was estimated by absorbance at 280 nm using an extinction coefficient of 1.1 of 0.1% (g/l) solution that was calculated on the basis of the frequency of tryptophan and tyrosine in vertebrate proteins. After digestion, the peptides were reconstituted in 0.5 M TEAB and labeled according to the manufacturer’s instructions for the iTRAQ 8-plex reagents (Applied Biosystems, Foster City, CA, USA).

#### Liquid chromatography (LC) separation and mass spectra (MS) quantification of peptides

iTRAQ labeled peptides were fractionated by SCX chromatography (Phenomenex, Inc, USA) using a Shimadzu LC-20AB HPLC Pump system. The peptides were eluted with a gradient of buffer A (10 mM NaH2PO4 in 25% ACN, pH 3.0) and buffer B (500 mM KCl, 10 mM KH2PO4 in 25% of ACN, pH 3.0). The specific fractionating procedures were as follows: a flow rate of 1 ml/min with a gradient of 0–8% buffer B for 22 min, 8–52% buffer B for 30 min, 52–100% buffer B for 50 min, 100% buffer B for 55 min, and buffer B was reset to 0% after 58 min. The elution was monitored by absorbance at 214 nm, and fractions were collected every 1 min. Collected fractions were desalted on C18 Cartridges and concentrated by vacuum centrifugation. Nano LC-MS/MS analysis on each of these fractions was performed using a Q-Exactive mass spectrometer (Thermo Fisher Scientific Inc. Rockford, IL., USA) equipped with nanoelectrospray ionization.

Peptides were identified by searching against the Uniprot_Bactrocera _dorsalis_19817_20160321.fasta (19817 sequences, download Mar 21, 2016) using a MS/MS data interpretation algorithm within MASCOT search engine (Matrix Science, London, UK; v.2.2) embedded into Proteome Discoverer 1.4 (Thermo, Pittsburgh, USA). When the MASCOT software was used to search the database, 4622 proteins were identified with a false discovery rate (FDR) of less than 1%. Differential expression ratios for proteins were obtained from MASCOT software (http://www.matrixscience.com), which calculates protein ratios using only ratios from the spectra that are distinct for each protein and excluding the shared peptides of protein isoforms. To calculate the differential expression ratios, all identified spectra from a protein were used to obtain an average protein ratio that was relative to the control label. Student t-test was used to analyze the differentially expressed proteins in antennae between the responsive and non-responsive B. dorsalis male flies. We only identified proteins with P-values <0.05 and fold changes ≥1.20 or ≤0.83 as being differentially expressed60.

#### Proteins functional and pathway enrichment analyses

Functional annotation of the proteins identified in B. dorsalis was carried out using Blast2GO, an integrated GO annotation and data mining tool that assigns GO through BLAST searches against online protein databases (http://www.geneontology.org)61,62. GO enrichment analysis was performed to provide all the GO terms that were significantly enriched in the differentially expressed proteins. Specifically, gene ontology enrichment analyses were performed using a hypergeometric test to map all differentially expressed proteins to GO terms in the database (http://www.geneontology.org/). The test uses the following equation (2):

$$P=1-\sum _{i=0}^{{\rm{m}}-1}\frac{(\begin{array}{c}M\\ i\end{array})(\begin{array}{c}N-M\\ n-i\end{array})}{(\begin{array}{c}N\\ n\end{array})}$$
(2)

where N is the number of all proteins with GO annotation; n is the number of differentially expressed proteins in N; M is the number of all proteins annotated to certain GO terms; and m is the number of differentially expressed proteins in M. The calculated P-value was first subjected to a Bonferroni correction, taking a corrected P-value of 0.05 as a threshold for statistical significance. GO terms fulfilling this condition were defined as significantly enriched GO terms in differentially expressed proteins. All identified proteins were mapped to pathway in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database using KEGG Automatic Annotation Server software (http://www.genome.ad.jp/kegg/). Significantly enriched metabolic pathways or signal transduction pathways in differentially expression proteins were identified using the same calculating formula as in GO analysis. Here N represents the number of all proteins with KEGG annotation, n is the number of differentially expressed proteins in N, M is the number of all proteins annotated to specific pathways, and m is number of differentially expressed proteins in M.

### Relative expression analysis of the candidate olfactory proteins at transcript levels by quantitative real-time PCR

The expression of protein transcripts was studied by qRT-PCR to validate the proteomic results. Specific olfaction genes that encode the olfactory proteins were detected following methods described previously2,63. An RNA extraction kit (Takara Biotechnology Co., Ltd., Japan) was used to extract the total RNA from the non-responsive and responsive male files’ antennae following the manufacturer’s instructions, and a gDNA eliminator spin column removed the genomic DNA. RNA was quantified by measuring the absorbance at 260 nm in a spectrophotometer (Thermo Nano DropTM 2000c; Santa Clara, USA). The purity of all RNA samples was assessed at an absorbance ratio of OD260/280, while the integrity of RNA was verified by 1% agarose gel electrophoresis. The cDNA was synthesized by using a PrimeScriptTM RT reagent Kit (Takara Biotechnology Co., Ltd., Japan), subsequently serving as a template for qRT-PCR.

The gene-specific primers were designed in Primer 5.0 (Premier, Canada) (Table S2). These were used to conduct quantitative real-time PCR to detect the differential expression olfaction proteins of responders compared with non-responders at the transcription level. By using a SYBR Premix ExTaq kit (Tiangen, Guangzhou, China) following the manufacturer’s instructions, with a Stratagene Mx3000 P thermal cycler (Agilent Technologies, Wilmington, Germany). The PCR efficiency of the genes was validated before gene expression analysis. The PCR master mix (20 μL) contained 10 μL of SYBR Green Supermix, 1 μL of cDNA templates, 1 μL of each forward and reverse primer (1 μmol/L), and 7 μL of double-distilled water. The following thermal program was executed: 95 °C for 15 min, followed by 40 cycles of 95 °C for 10 s, 55 °C for 20 s, 72 °C for 20 s, and a final melting cycle (from 55 °C to 95 °C). Three biological and three technical replicates were used and performed for each experiment, respectively. The α-tubulin gene of B. dorsalis was used as an internal control (GenBank accession number: XM_011212814). The relative gene expression levels were calculated by using the 2−ΔΔCT method as described previously64.

### Sequence analysis of BdorOBP2

The sequences information of BdorOBP2 was obtained from the National Center for Biotechnology Information (NCBI) database and submitted by Zheng et al.1 (GenBank accession number: KC559113). The amino acid sequences of BdorOBP2 identified in this study, along with the 91 known OBPs from other Tephritidae insects, were used to construct a phylogenetic tree with the MEGA 7.0 software (Molecular Evolutionary Genetics Analysis, Version 4.0, Sudhir Kumar, USA) and Interactive Tree Of Life (iTOL) web tool (http://itol.embl.de). The highly divergent signal peptides in the N terminal were removed. Branch support was assessed by the Neighbor Joining technique with 1000 bootstrap replications.

### Responsiveness of immature and mature flies to ME, and BdorOBP2 expression analysis by quantitative real-time PCR

Each sample of 200 3-day- and 15-day-old males, as well as the 15-day-old females of the control line, were randomly selected as the test subjects and released into a screen cage (1.0 m × 1.0 m × 1.0 m) lacking a trap. Approximately 30 min later, a fly trap containing 1.0 mL of pure ME was placed inside the screen cage. For the controls, traps were also placed but without ME. After trapping for 2 h, we removed the traps and counted the number of attracted flies. Three biological replicates were performed. Additionally, from the untested individuals in the control line, the total RNA of antennae from the 3-day- and 15-day-old males, and the 15-day-old females were extracted to determine their expression of BdorOBP2 (by using qRT-PCR, as described above).

### Effect of ME regulation on BdorOBP2 expression in male antennae

The protocol followed published procedures14. Each sample of 200 mature males of the control line were randomly selected and placed in a cage (35 cm × 35 cm × 35 cm) without a trap. After being starved for 12 h, 500 μL of ME (1:1 dilution with mineral oil), was applied on a filter paper spotted in a petri dish (diam. = 3.5 cm) to feed the hungry male flies, while the control groups were treated with an equal volume of just the mineral oil. Total RNA was extracted from the antennae after treatment for 1 h and 2 h, then reverse transcribed into single-chain cDNAs. Next, the qRT-PCR quantitative technique was applied to analyze the expression of BdorOBP2 after the ME treatment. Each treatment was replicated three times.

### RNA interference knock-down of BdorOBP2

An RNA interference experiment was performed to demonstrate the role of BdorOBP2 in ME perception. Double-stranded RNA (dsRNA) was synthesized and purified according to the MEGAscript® RNAi Kit (Thermo Fisher Scientific, USA). Primers used to synthesize dsRNA are listed in Table S3. The GFP gene was used as a control dsRNA (dsGFP) (GenBank accession number: AHE38523). The dsRNA expression and injection procedure for the male flies were performed as reported in previous studies65,66,67. Both dsBdorOBP2 and dsGFP were then purified and re-suspended in 1 × Injection Buffer (10 mM Tris-HCl pH 7, 1 mM EDTA). Finally, the quality of the dsRNA was determined by 2% TBE gel, and its concentration determined by a Nanodrop 1000 (Thermo, USA).

Each sample of 50 mature males (15 days-old) were randomly selected and placed into a cage (35 cm × 35 cm × 35 cm). The treatment males were then injected with 400 nL of dsBdorOBP2 (2000 ng/μL). Negative control groups were injected with an equal volume of dsGFP and 1 × Injection Buffer. The male flies were normally fed to serve as the blank control groups. All flies were reared in the cage and supplied with identical food and water. After treatment for 24 h and 48 h, the respective numbers of dead adults were counted. For the bioassay, males from the dsBdorOBP2 treatment groups, the negative control groups and the blank control groups were placed inside a screen cage (1.0 m × 1.0 m × 1.0 m) equipped with ME as a trap. The lured males were counted after 2 h, and the silencing efficiency of BdorOBP2 was detected by qRT-PCR. Each experiment was replicated five times, and 50 flies were tested in each repeat made.

### Statistical analysis

All the results from experimental replicates were expressed as the mean ± SE and analyzed by SAS 9.20 software (SAS Institute Inc. Cary. NC). All data were analyzed using Shapiro-Wilk and Levene’s tests for normal distribution and homogeneity of variances, respectively. If data were normally distributed and had similar variances, then means of measured variables were compared by one-way analysis of variance (ANOVA). Significant ANOVA results multiple comparison over three groups were assessed by Duncan’s multiple range test (DMRT, P = 0.05), and two-sample analysis was performed using Student’s t-test (P = 0.05). Non-normally distributed data were analyzed using a nonparametric Kruskal-Wallis test to compare medians; differences significant at the 0.05 significance level were subjected to a Mann-Whitney test for pairwise comparisons. P < 0.05 was considered to be statistically significant. Results were plotted with Origin 9.0.

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## References

1. 1.

Zheng, W. W. et al. Identification and Expression Profile Analysis of Odorant Binding Proteins in the Oriental Fruit Fly Bactrocera dorsalis. Int. J. Mol. Sci. 14, 14936–14949 (2013).

2. 2.

Liu, Z., Smagghe, G., Lei, Z. & Wang, J. J. Identification of male- and female-specific olfaction genes in antennae of the oriental fruit fly (Bactrocera dorsalis). PLoS One. 11, e0147783 (2016).

3. 3.

Clarke, A. R. et al. Invasive phytophagous pests arising through a recent tropical evolutionary radiation: The Bactrocera dorsalis complex of fruit flies. Annu. Rev. Entomol. 50, 293–319 (2005).

4. 4.

Stephens, A., Kriticos, D. & Leriche, A. The current and future potential geographical distribution of the oriental fruit fly, Bactrocera dorsalis (Diptera: Tephritidae). Bull. Entomol. Res. 97, 369–378 (2007).

5. 5.

Shen, G. M., Wang, X. N., Dou, W. & Wang, J. J. Biochemical and molecular characterization of acetylcholinesterase in four field populations of Bactrocera dorsalis (Hendel) (Diptera:Tephritidae). Pest. Manag. Sci. 68, 1553–1563 (2012).

6. 6.

Vargas, R. I. & Prokopy, R. Attraction and feeding responses of melon flies and oriental fruit flies (Diptera: Tephritidae) to various protein baits with and without toxicants. Hawaii. Ent. Soc. 38, 49–60 (2006).

7. 7.

Pagadala, D. K. J., Christine, M. W., John, C., Michael, A. B. & Toby, J. A. B. Isolation and identification of host cues from Mango, Mangifera indica, that attract gravid female oriental fruit fly, Bactrocera dorsalis. J. Chem. Ecol. 38, 361–369 (2012).

8. 8.

Shelly, T. E. Zingerone and the mating success and field attraction of male melon flies (Diptera: Tephritidae). J. Asia-Pac. Entomol. 20, 175–178 (2017).

9. 9.

Shelly, T. E., Edu, J. & McInnis, D. Prerelease consumption of methyl eugenol increases the mating competitiveness of sterile males of the oriental fruit fly, Bactrocera dorsalis, in large field enclosures. J. Insect. Sci. 10, 1–16 (2010).

10. 10.

Miller, E. C. et al. Structure-activity studies of the carcinogenicities in the mouse and rat of some naturally occurring and synthetic alkenyl benzene derivatives related to safrole and estragole. Cancer. Res. 43, 1124–1134 (1983).

11. 11.

Shelly, T. E. Selection for non-responsiveness to methyl eugenol in male oriental fruit flies (Diptera: Tephritidae). Fla. Entomol. 80, 248–253 (1997).

12. 12.

Smith, R. L. et al. Safety assessment of allylalkoxybenzene derivatives used as flavouring substances-methyl eugenol and estragole. Food. Chem. Toxicol. 40, 851–870 (2002).

13. 13.

Khrimian, A. et al. Ring-Fluorinated analog of methyl eugenol: attractiveness to and metabolism in the oriental Fruit Fly, Bactrocera dorsalis (Hendel). J. Chem. Ecol. 35, 209–218 (2009).

14. 14.

Zheng, Q. Q. et al. Odorant receptor co-receptor Orco is upregulated by methyl eugenol in male Bactrocera dorsalis (Diptera: Tephritidae). J. Insect. Physiol. 58, 1122–1127 (2012).

15. 15.

Khrimian, A. P., Demilo, A. B., Waters, R. M., Liqyido, N. J. & Nicholson, J. M. Monofluoroanalogs of eugenol methyl ether as novel attractants for the oriental fruit fly. J. Org. Chem. 59, 8034–8039 (1994).

16. 16.

Khrimian, A., Jang, E. B., Nagata, J. & Carvalho, L. Consumption and metabolism of 1, 2-Dimethoxy-4- (3-Fluoro-2-Propenyl) benzene, a fluorine analog of methyl eugenol, in the oriental fruit fly, Bactrocera dorsalis (Hendel). J. Chem. Ecol. 32, 1513–1526 (2006).

17. 17.

Jang, E. B., Khrimian, A. & Siderhurst, M. S. Di- and Tri-fluorinated analogs of methyl eugenol: attraction to and metabolism in the oriental fruit fly, Bactrocera dorsalis (Hendel). J. Chem. Ecol. 37, 553–564 (2011).

18. 18.

Kempraj, V. et al. Computational reverse chemical ecology: Virtual screening and predicting behaviorally active semiochemicals for Bactrocera dorsalis. BMC Genomics. 15, 209 (2014).

19. 19.

Justice, R. W., Biessmann, H., Walter, M. F., Dimitratos, S. D. & Woods, D. F. Genomics spawns novel approaches to mosquito control. Bioessays. 25, 1011–1020 (2003).

20. 20.

Leal, W. S. Odorant reception in insects: roles of receptors, binding proteins, and degrading enzymes. Annu. Rev. Entomol. 58, 373–391 (2013).

21. 21.

Vogt, R. G., Callahan, F. E., Rogers, M. E. & Dickens, J. C. 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. 24, 481–495 (1999).

22. 22.

Hallem, E. A., Dahanukar, A. & Carlson, J. R. Insect odor and taste receptors. Annu. Rev. Entomol. 51, 113–135 (2006).

23. 23.

Pelosi, P., Calvello, M. & Ban, L. Diversity of odorant-binding proteins and chemosensory proteins in insects. Chem. Senses. 30, i291–i292 (2005).

24. 24.

Sato, K., Pellegrino, M., Nakagawa, T., Vosshall, L. B. & Touhara, K. Insect olfactory receptors are heteromeric ligand-gated ion channels. Nature. 452, 1002–1006 (2008).

25. 25.

Taylor, A. J., Cook, D. J. & Scott, D. J. Role of odorant binding proteins: Comparing hypothetical mechanisms with experimental data. Chemosens. Percept. 1, 153–162 (2008).

26. 26.

Silbering, A. F. et al. Complementary function and integrated wiring of the evolutionarily distinct Drosophila Olfactory subsystems. J. Neurosci. 31, 13357–13375 (2011).

27. 27.

Siciliano, P. et al. Identification of pheromone components and their binding affinity to the odorant binding protein CcapOBP83a-2 of the Mediterranean fruit fly, Ceratitis capitata. Insect. Biochem. Molec. 48, 51–62 (2014).

28. 28.

Gotzek, D., Robertson, H. M., Wurm, Y. & Shoemaker, D. Odorant binding proteins of the red imported fire ant, Solenopsis invicta: an example of the problems facing the analysis of widely divergent proteins. PLoS One. 6, e16289 (2011).

29. 29.

Vandermoten, S., Francis, F., Haubruge, E. & Leal, W. S. Conserved odorant-binding proteins from aphids and eaves dropping predators. PLoS One. 6, e23608 (2011).

30. 30.

Tsitsanou, K. et al. Anopheles gambiae odorant binding protein crystal complex with the synthetic repellent DEET: implications for structure-based design of novel mosquito repellents. Cell. Mol. Life. Sci. 69, 283–297 (2012).

31. 31.

Biessmann, H. et al. The Anopheles gambiae odorant binding protein 1 (AgamOBP1) mediates indole recognition in the antennae of female mosquitoes. PLoS One. 5, e9471 (2010).

32. 32.

Pelletier, J., Guidolin, A., Syed, Z., Cornel, A. J. & Leal, W. S. Knockdown of a mosquito odorant-binding protein involved in the sensitive detection of oviposition attractants. J. Chem. Ecol. 36, 245–248 (2010).

33. 33.

Wu, Z. Z., Lin, J. T., Zhang, H. & Zeng, X. N. BdorOBP83a-2 mediates responses of the oriental fruit fly to semiochemicals. Front. Physiol. 7, 452 (2016).

34. 34.

Karp, N. A. et al. Addressing accuracy and precision issues in iTRAQ quantitation. Mol. Cell. Proteom. 9, 1885–1897 (2010).

35. 35.

Xu, D. X., Sun, L., Liu, S. L., Zhang, L. B. & Yang, H. S. Understanding the heat shock response in the sea cucumber Apostichopus japonicus, using iTRAQ-based proteomics. Int. J. Mol. Sci. 17, 150 (2016).

36. 36.

Yang, N. et al. Transcriptomic and proteomic responses of sweetpotato whitefly, Bemisia tabaci, to thiamethoxam. PLoS One. 8, e61820 (2013).

37. 37.

Tu, X. et al. Transcriptomic and proteomic analysis of pre-diapause and nondiapause eggs of migratory locust, Locusta migratoria L. (Orthoptera: Acridoidea). Sci. Rep. 5, 11402 (2015).

38. 38.

Wang, W. J. et al. Identification of proteins associated with pyrethroid resistance by iTRAQ-based quantitative proteomic analysis in Culex pipiens pallens. Parasit. Vectors. 8, 95 (2015).

39. 39.

Guo, Q. L. et al. Choosing male Bactrocera dorsalis (Hendel) without taxis to methyl eugenol. Chin. J. Trop. Crop. 31, 845–848 (2010).

40. 40.

Karunaratne, M. M. S. C. & Karunaratne, U. K. P. R. Factors influencing the responsiveness of male oriental fruit fly, Bactrocera dorsalis, to methyl eugenol (3,4-dimethoxyalyl benzene). Trop. Agric. Res. Exten. 15, 92–97 (2012).

41. 41.

Miller, C. L. et al. Role of Ca2+/Calmodulin-stimulated cyclic nucleotide phosphodiesterase 1 in mediating cardiomyocyte hypertrophy. Circ. Res. 105, 956–964 (2009).

42. 42.

Feng, M. et al. Antennal proteome comparison of sexually mature drone and forager honeybees. J. Proteome. Res. 10, 3246–3260 (2011).

43. 43.

Vogel, C. & Marcotte, E. M. Insights into the regulation of protein abundance from proteomic and transcriptomic analyses. Nat. Rev. Genet. 13, 227–232 (2012).

44. 44.

Zhang, L. Z. et al. Differential protein expression analysis following olfactory learning in Apis cerana. J. Comp. Physiol. A. 201, 1053–1061 (2015).

45. 45.

Zhou, J. J. et al. Characterisation of Bombyx mori Odorant-binding proteins reveals that a general odorant-binding protein discriminates between sex pheromone components. J. Mol. Biol. 389, 529–545 (2009).

46. 46.

He, X. L. et al. Binding of the general odorant binding protein of Bombyx mori BmorGOBP2 to the moth sex pheromone components. J. Chem. Ecol. 36, 1293–1305 (2010).

47. 47.

Yin, J. et al. Functional analysis of general odorant binding protein 2 from the meadow moth, Loxostege sticticalis L. (Lepidoptera: Pyralidae). PLoS One. 7, e33589 (2012).

48. 48.

Feng, L. & Prestwich, G. D. Expression and characterization of a lepidopteran general odorantbinding protein. Insect. Biochem. Molec. 27, 405–412 (1997).

49. 49.

Gong, Z. J. et al. Cloning, expression and functional analysis of a general odorant-binding protein 2 gene of the rice striped stem borer, Chilo suppressalis (Walker) (Lepidoptera: Pyralidae). J. Mol. Biol. 18, 405–417 (2009).

50. 50.

Rebijith, K. B. et al. RNA Interference of Odorant-binding protein 2 (OBP2) of the cotton aphid, Aphis gossypii (Glover), resulted in altered electrophysiological responses. Appl. Biochem. Biotechnol. 178, 251–266 (2016).

51. 51.

Benton, R., Sachse, S., Michnick, S. W. & Vosshall, L. B. A typical membrane topology and heteromeric function of Drosophila odorant receptors in vivo. PLoS Biol. 4, e20 (2006).

52. 52.

Patch, H. M., Velarde, R. A., Walden, K. K. & Robertson, H. M. A candidate pheromone receptor and two odorant receptors of the hawkmoth Manduca sexta. Chem. Senses. 34, 305–316 (2009).

53. 53.

Vosshall, L. B. & Hansson, B. S. A unified nomenclature system for the insect ofactory coreceptor. Chem. Senses. 36, 497–498 (2011).

54. 54.

Stengl, M. & Funk, N. W. The role of the coreceptor Orco in insect olfactory transduction. J. Comp. Physiol. A Neuroethol. Sens. Neural Behav. Physiol. 199, 897–909 (2013).

55. 55.

Larsson, M. C. et al. Or83b encodes a broadly expressed odorant receptor essential for Drosophila olfaction. Neuron. 43, 703–714 (2004).

56. 56.

Hallem, E. A. & Carlson, J. R. Coding of odors by a receptor repertoire. Cell. 125, 143–160 (2006).

57. 57.

Chang, C. L., Vargas, R. I., Caceres, C., Jang, E. & Cho, I. K. Development and assessment of a liquid larval diet for Bactrocera dorsalis (Diptera: Tephritidae). Ann. Entomol. Soc. Am. 99, 1191–1198 (2006).

58. 58.

Su, Y. et al. Alteration of intracellular protein expressions as a key mechanism of the deterioration of bacterial denitrification caused by copper oxide nanoparticles. Sci. Rep. 5, 15824 (2015).

59. 59.

Han, J. Z. et al. iTRAQ-based proteomic analysis of LI-F type peptides produced by Paenibacillus polymyxa JSa-9 mode of action against Bacillus cereus. J. Proteomics. 150, 130–140 (2017).

60. 60.

Briolant, S. et al. Research Plasmodium falciparum proteome changes in response to doxycycline treatment. Malar. J. 9, 141–154 (2010).

61. 61.

Conesa, A. et al. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 21, 3674–3676 (2005).

62. 62.

Götz, S. et al. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic. Acids. Res. 36, 3420–3435 (2008).

63. 63.

Nakamura, A. M. et al. Reference genes for accessing differential expression among developmental stages and analysis of differential expression of OBP genes in Anastrepha obliqua. Sci. Rep. 6, 17480 (2016).

64. 64.

LivaK, K. J. & Schmittgen, T. D. Analysis of relative gene expression data using real-time quantitive PCR and the 2−ΔΔCT method. Methods. 25, 402–408 (2001).

65. 65.

Li, X., Zhang, M. & Zhang, H. RNA interference of four genes in adult Bactrocera dorsalis by feeding their dsRNAs. PLoS One. 6, e17788 (2011).

66. 66.

Liu, G., Wu, Q., Li, J., Zhang, G. & Wan, F. RNAi-mediated knock-down of transformer and transformer 2 to generate male-only progeny in the oriental fruit fly, Bactrocera dorsalis (Hendel). PLoS One. 10, e0128892 (2015).

67. 67.

Dong, Y. C. et al. Bactrocera dorsalis male sterilization by targeted RNA interference of spermatogenesis: empowering sterile insect technique programs. Sci. Rep. 6, 35750 (2016).

## Acknowledgements

We sincerely thank the Charlesworth Group for correction of the manuscript language and the Shanghai Applied Protein Technology Co. Ltd for technology support and helpful advice. We are grateful to the members of our laboratory for their cooperation in oriental fruit fly rear and treatment. This research was supported by grants from Guangdong Province College High-Quality Professional Foundation (No. 246) and National key research and development project of China (No. 2016YFC1201200). The funders have no role in study design, data collection and analysis, decision to publish or preparation of the manuscript.

## Author information

### Affiliations

1. #### Department of Entomology, South China Agricultural University, Guangzhou, 510642, China

• Huan Liu
• , Xiao-Feng Zhao
• , Yi-Ye Han
• , Jin Chen
•  & Yong-Yue Lu

• Lang Fu

### Contributions

H.L. and Y.Y.L. conceived and designed the experiments. H.L., X.F.Z., L.F., Y.Y.H. and J.C. performed the experiments. H.L. analyzed the data. Y.Y.H. and J.C. provided material support. H.L. and Y.Y.L. wrote the manuscript. All authors discussed the results and reviewed the manuscript.

### Competing Interests

The authors declare that they have no competing interests.

### Corresponding author

Correspondence to Yong-Yue Lu.

## Electronic supplementary material

### DOI

https://doi.org/10.1038/s41598-017-15893-6
{"extraction_info": {"found_math": true, "script_math_tex": 0, "script_math_asciimath": 0, "math_annotations": 0, "math_alttext": 0, "mathml": 0, "mathjax_tag": 0, "mathjax_inline_tex": 0, "mathjax_display_tex": 1, "mathjax_asciimath": 0, "img_math": 0, "codecogs_latex": 0, "wp_latex": 0, "mimetex.cgi": 0, "/images/math/codecogs": 0, "mathtex.cgi": 0, "katex": 0, "math-container": 0, "wp-katex-eq": 0, "align": 0, "equation": 0, "x-ck12": 0, "texerror": 0, "math_score": 0.6535671353340149, "perplexity": 12831.219376377672}, "config": {"markdown_headings": true, "markdown_code": true, "boilerplate_config": {"ratio_threshold": 0.18, "absolute_threshold": 10, "end_threshold": 15, "enable": true}, "remove_buttons": true, "remove_image_figures": true, "remove_link_clusters": true, "table_config": {"min_rows": 2, "min_cols": 3, "format": "plain"}, "remove_chinese": true, "remove_edit_buttons": true, "extract_latex": true}, "warc_path": "s3://commoncrawl/crawl-data/CC-MAIN-2019-09/segments/1550247495367.60/warc/CC-MAIN-20190220170405-20190220192405-00175.warc.gz"}