- Research article
- Open Access
Molecular evolution of Odorant-binding proteins gene family in two closely related Anastrepha fruit flies
© The Author(s). 2016
- Received: 7 April 2016
- Accepted: 28 September 2016
- Published: 7 October 2016
Odorant-binding proteins (OBPs) are of great importance for survival and reproduction since they participate in initial steps of the olfactory signal transduction cascade, solubilizing and transporting chemical signals to the olfactory receptors. A comparative analysis of OBPs between closely related species may help explain how these genes evolve and are maintained under natural selection and how differences in these proteins can affect olfactory responses. We studied OBP genes in the closely related species Anastrepha fraterculus and A. obliqua, which have different host preferences, using data from RNA-seq cDNA libraries of head and reproductive tissues from male and female adults, aiming to understand the speciation process occurred between them.
We identified 23 different OBP sequences from Anastrepha fraterculus and 24 from A. obliqua, which correspond to 20 Drosophila melanogaster OBP genes. Phylogenetic analysis separated Anastrepha OBPs sequences in four branches that represent four subfamilies: classic, minus-C, plus-C and dimer. Both species showed five plus-C members, which is the biggest number found in tephritids until now. We found evidence of positive selection in four genes and at least one duplication event that preceded the speciation of these two species. Inferences on tertiary structures of putative proteins from these genes revealed that at least one positively selected change involves the binding cavity (the odorant binding region) in the plus-C OBP50a.
A. fraterculus and A. obliqua have a bigger OBP repertoire than the other tephritids studied, though the total number of Anastrepha OBPs may be larger, since we studied only a limited number of tissues. The contrast of these closely related species reveals that there are several amino acid changes between the homologous genes, which might be related to their host preferences. The plus-C OBP that has one amino acid under positive selection located in the binding cavity may be under a selection pressure to recognize and bind a new odorant. The other positively selected sites found may be involved in important structural and functional changes, especially ones in which site-specific changes would radically change amino acid properties.
- Sister species
- Recent speciation event
- Phylogenetic analyses
- OBP subfamilies
- OBP putative tertiary structure
The study of the genetics of species differences generally requires the identification of fixed genes between lineages. Most of these genes are simply “ordinary” traits that diverged between lineages with no direct role on isolation , while some could be involved with new changes or adaptations driven by selection. Up until recently, studies on speciation would then have to investigate these genes by interspecific crosses or other means to assess reproductive isolation  but new technical advances have enabled their identification by looking at selection at the gene level. Genes involved in reproduction and mate choice tend to evolve rapidly and display signatures of positive selection [3, 4], as well as genes facilitating chemoreception [5–8], because they are needed to interpret information about the environment, such as the presence of food or predators. Olfactory information, specifically, also controls social and sexual interactions between individuals of the same species, such as the detection of odors and pheromones essential for survival and reproduction . In insects, the solubilization and transport of chemical signals through the aqueous lymph of insect’s sensilla to the olfactory receptors is the initial step of the transduction cascade of olfactory signals, and is mediated by the Odorant-binding proteins – OBPs .
OBPs are small soluble globular proteins (10–30 kDa), with six highly conserved cysteine residues (referred to as C1 to C6) in characteristic positions of all known insect OBPs. These cysteines form three disulphide bonds that stabilize their tertiary structure  and help define a hydrophobic binding cavity [12–14]. Different OBPs have specific affinities to odors and are present in distinctive portions of the sensilla, suggesting that they participate in odor detection by restricting the available spectrum of odors into adjacent receptors . OBPs’ conserved cysteine residues, associated with functional information and phylogenetic relationships, have been used to classify OBPs into subfamilies: classic OBPs (six cysteines), plus-C OBPs (more than six cysteines) and minus-C OBPs (some members with less than six cysteines) [15–17]. Additionally, OBP transcripts encoding two complete OBP domains were identified in Drosophila melanogaster, and classified as a subfamily named dimer . Other subfamilies were also defined in Drosophila: PBP, ABPI, ABPII, CRLBP and D7 that are unevenly distributed among dipterans [15, 18, 19]. It is expected that more generalist species would have a larger OBP repertoire to recognize various plant chemicals, in contrast to more specialist species, in which the power of natural selection to maintain a large OBP repertoire may have been weaker .
OBP genes are part of a gene family with low sequence similarity among its members, but high conservation at the structural level . OBP genes are located in large gene clusters, suggesting that they arose by tandem duplication  and evolve mainly via the birth-and-death model , in which newly duplicated members progressively diverge in sequence and function, or may be lost to deletion or pseudogenization. The high disparity of gene sequences among OBPs implies a rapid rate of evolution in this gene family, suggesting that these genes might have evolved under the influence of positive selection [5, 8, 24], whereby even the conserved cysteine residues may be lost .
Here, we studied OBP genes in two closely related species of Tephritidae, Anastrepha fraterculus and A. obliqua, which are important fruit pests in South America. These species have diverged recently and exhibit a limited number of morphological  and genetic distinguishing characters . Though A. fraterculus has been associated with a wide number of hosts, it prefers several Myrtaceae fruits , being considered one of the main economic pests in South America. A. obliqua, on the other hand, though an important pest species as well, has been associated to a smaller number of hosts, several of those Anacardiaceae . Because OBPs are important targets for natural and sexual selection, their role in host and mate choice has previously been established in several species [29, 30]. The investigation of OBP evolution in these closely related species may provide clues about this group’s differentiation and host preference. We identified A. fraterculus and A. obliqua OBP members using RNA-seq data from reproductive and cephalic tissues of several different reproductive stages in male and female adults: before and after copulation and females after oviposition, and analyzed the patterns of molecular evolution between these two species. Several OBP genes were found to be evolving under positive selection and we speculate how these amino acid changes would affect protein structure and their consequences for adaptation.
Transcriptome libraries, assembly and annotation
We used populations of A. fraterculus from Southeast (22° 01′ 03″ S, 47 o 53′ 27″ W) and of A. obliqua from Midwest (16° 41′ 58″ S, 49 o 16′ 35″ W) regions of Brazil, that have been maintained in a controlled environment room at 25 °C ± 5 °C (60–90 % humidity) and natural photoperiod. Transcriptome libraries were generated for each species separately, using two reproductive stages for both sexes (virgin and post-mating), and a third one for females (post-oviposition). All profiles were made with biological replication, totaling 10 reproductive profiles per species. For each profile, we extracted the total RNA from head and reproductive tissues, totaling 20 libraries for each species. We used a pool of 10 individuals per library, amounting to 200 individuals in total.
More details on sample preparation, molecular procedures, assembly and annotation of the RNA-seq data are described elsewhere . Briefly, total RNA was extracted using the TRizol/chloroform protocol . RNA-seq libraries were constructed from four μg of total RNA using the TruSeq Stranded Total RNA Sample Prep Kit (Illumina) protocol, according to the manufacturer’s instructions. Pools of 12 libraries were run on an Illumina HiSeq2000 on a lane with runs of 2 × 100 bp paired-end reads. All reads were trimmed for quality and length with SeqyClean , keeping only reads with a minimum sequence length of 50 pb, a minimum of 0.01 for the parameter ‘max-avg-error’ and 0.05 for ‘max-error-at-ends’, and an average Phred quality score ≥ 20. Processed reads were assembled in two independent transcriptomes, one for each species, using the Trinity short read assembler (release 2013-02-25) , using default parameters.
Odorant-binding proteins identified in A. fraterculus and A. obliqua transcriptomes
A. fraterculus OBPs
ORF length (aa)
Signal Peptide (aa)
A. obliqua OBPs
ORF length (aa)
Signal Peptide (aa)
Alignments and phylogenetic analysis
We initially performed multiple alignment of the inferred OBP amino acid sequences for each species with MAFFT , using default settings. Average amino acid identities were obtained from the Percent Identity Matrix in the MAFFT alignment results. The multiple alignment of nucleotide sequences was manually adjusted in Bioedit v.184.108.40.206 , using the amino acid alignment as a guide. Three different methods, MaxChi , GENECONV  and RDP , were implemented in the RDP3 program  to investigate for recombination events.
Because OBPs are very divergent, we inferred the phylogenetic relationships in two steps. First, we aligned all A. fraterculus and A. obliqua OBPs here identified with OBP sequences of Ceratitis capitata and of D. melanogaster (obtained from GenBank - Additional file 1), using MAFFT as previously described. These species were chosen because D. melanogaster has the best curated genome and C. capitata is the closest species to Anastrepha with OBP sequences available. We used jModelTest  to estimate the best-fitting nucleotide model of substitution, that was used to infer Maximum likelihood phylogenetic relationships among OBPs with PhyML ver.3.0 . This first phylogenetic tree was reconstructed in order to corroborate the annotation and, consequently, the allocation in subfamilies.
Based on the confirmation of the subfamilies’ division in our first phylogenetic tree, we performed a new alignment combining all A. fraterculus and A. obliqua OBPs and included OBP sequences from the tephritids C. capitata, Bactrocera dorsalis , Rhagoletis pomonella , and Rhagoletis suavis . This alignment was performed with MAFFT separating and aligning OBP sequences by subfamily, in such a way that the alignments by subfamilies were combined in one general alignment using the Profile Alignment Mode in ClustalX2 , which produced a better alignment and a second phylogenetic tree, which was inferred as previously described. MAFFT also provided us with an identity matrix, which we used to estimate the average amino acid identity by subfamily. Due to the high divergence in the N-terminal region, for the phylogenetic and evolutionary analysis we removed the region before the first cysteine residue in all sequences prior to the alignments, as it was done in other studies [19, 50–52].
Since OBPs from different subfamilies show great divergence, we chose to perform the evolutionary analysis by subfamilies. We used separate alignments and phylogenetic trees for each one of the four OBP subfamilies using sequences of A. fraterculus, A. obliqua, C. capitata and D. melanogaster (Additional file 2). Before phylogenetic reconstruction, we tested the alignments for sequence saturation using DAMBE 5 . We investigated patterns of molecular evolution and positive selection with the strict branch-site test comparing models A vs. A-null , using the software CODEML, implemented in PAML v.4 . The comparison between models was tested using likelihood-ratio tests (LRTs) for hierarchical models , and we used a Bonferroni correction for the number of branches tested. A significantly higher likelihood for the alternative model than that of the null model would indicate evidence of positive selection. We also estimated pairwise Ka, Ks rates and its ratio Ka/Ks, using the KaKs_Calculator , for the putative orthologous genes using the MS model . Genes showing Ka/Ks rates higher than 0.5 were considered as potentially evolving under positive selection .
The Bayes Empirical Bayes (BEB) method implemented in PAML was used to estimate the posterior probability that a given site is evolving under positive selection. Furthermore, we analyzed whether amino acid replacements in positively selected sites alter physicochemical properties in the proteins using Conant-Stadler amino acid property set  with the PRIME method available in the Datamonkey web server . To make inferences about the position of the BEB positively selected sites on the tertiary structure of the proteins, we predicted the tertiary structure of Anastrepha OBPs under positive selection with PHYRE2 . The program used Ae. aegypti OBP1 (PDB: 3K1E, chain A ) as reference to infer the tertiary structure of Anastrepha OBP56h-1 and OBP56h-2, and A. gambiae OBP4 (PDB: 3Q8I, chain A ) and A. gambiae OBP48 (PDB: 4KYN, chain B ) as reference to OBP57c and OBP50a, respectively. Since gene sequences in both species were similar and because AfraOBP50a-1 and AfraOBP50a-2 were incomplete, we used A. obliqua sequences to make the predictions.
OBP genes in Anastrepha fraterculus and Anastrepha obliqua
We identified a similar number of sequences associated with OBPs in the two species: 23 different sequences in A. fraterculus and 24 in A. obliqua, which corresponded to 20 different D. melanogaster OBP genes (Table 1). In comparison, free-living insect OBP families described to date range from 11 genes in the beet armyworm Spodoptera exigua  to more than a hundred in some species of mosquitoes , although the parasitic body louse Pediculus humanus (which has a reduced genome), has only five OBPs . The length of complete open reading frames (Table 1) matched what has been described for several arthropod species . The longest OBPs are AfraOBP59a and AoblOBP59a. These sequences are unusual, because they have a specific region containing 134 amino acids that is situated between the first (C1) and the second (C2) cysteine, whereas others in the same subfamily have ~30 amino acids in this region. Even though an unusual length for OBP59a had already been demonstrated in other arthropods , in Anastrepha these OBPs are 106 amino acids longer than in D. melanogaster. Signal peptides, one of characteristic hallmarks of the OBP gene family, were predicted at the hydrophobic N-terminal for the almost all OBPs, except for AfraOBP19c and AoblOBP19c. Though not all OBPs described have signal peptides , we believe that the PrediSi program failed to find a match in the SwissProt database for both Anastrepha OBP19c signal peptide due to sequence divergence, not because these OBPs lack a signal peptide.
A comparison of putatively homologous OBPs in A. fraterculus and A. obliqua reveals great amino acid similarity, so much so that homologs of OBP19a and OBP99a had identical amino acid sequences in both species. In spite of that, the majority of orthologous OBPs differ by at least a few amino acids, similar to what was described for related species of aphids [52, 68]. However, even a single amino acid change may impact the overall 3D structure and/or the binding affinities of OBPs [69, 70]. For instance, a polymorphism in Obp57e was shown to be responsible for differences in host plant preferences between D. sechellia and D. melanogaster , and some polymorphisms were associated with natural variation in olfactory behavior in response to benzaldehyde in Obp99a, Obp99c, and Obp99d in D. melanogaster . Likewise, the few changes observed between A. fraterculus and A. obliqua OBPs may result in significant differences in olfactory responses, which is yet to be determined.
We did an alignment with all A. fraterculus OBP sequences, and another one with all A. obliqua OBP sequences. Two independent recombination tests were performed using these alignments, because recombination interferes with phylogenetic inferences and may generate higher rates of false positives in positive selection tests . Both tests performed failed to find significant evidence for recombination. These alignments revealed that even though putatively orthologous OBP copies in different species were very similar, there was great divergence among OBP genes, so much so that only the cysteine residues were conserved across all OBPs of the same species. Sequence divergence was higher than what has been described for other insects. The average amino acid identity among all OBPs was 16.65 % in A. fraterculus (ranging from 8.1 to 97.1 %) and 16.22 % in A. obliqua (varying from 8.0 to 98.7 %), compared with 20.4 % of identity between all OBPs in D. melanogaster , and 20 % between all OBPs found in Solenopsis invicta .
The higher intraspecific divergence among OBPs in A. fraterculus and A. obliqua, when compared with other insects, may suggest a greater selective pressure for OBP differentiation, or relaxed selection, which also led us to investigate for patterns of evolution in these genes. On the other hand, the high amino acid similarity found between the homologs in A. fraterculus and A. obliqua probably indicates that the recent time of divergence between these species may not have been enough to drastically alter their sequences, which could be a reflection of retention of similar physiological functions across species or simply of the recent time since their divergence. Since these species have diverged recently and have accumulated only a few changes, the current methodologies to investigate for positive selection between these species are not very efficient, but we may use these data in contrast to other Tephritidae and Drosophilidae to investigate for patterns of selection in OBPs as a whole. We chose to investigate Anastrepha OBPs because they may be involved with the group’s differentiation since they may be directly involved with host preference and mate choice. Therefore, patterns of selection in these genes may help us understand evolution and speciation not only in the genus Anastrepha, but also across this important family, which may make it useful in pest control programs.
Phylogenetic relationships among OBP genes
Attributes of Anastrepha Odorant-binding proteins subfamilies
Number of members
ORF length (aa)
Protein core region
Amino acid identity (%)
C1 - X24–33 - C2 - X3 - C3 - X35–47 - C4 - X8–18 - C5 - X8 - C6
C1 - X26 - C2 - X3 - C3 - X39–40 - C4 - X10 - C5 - X8 - C6 and C1 - X23–30 - C3 - X38 - C4 - X18–19 - C6
C1a - X11–13 - C1b - C1c - C1 - X21–48 - C2 - X3 - C3 - X43 - C4 - X20–33 - C5 - X9 - C6 - X8 - C6a - X10 - C6b - X9 - C6c
C1 - X28–29 - C2 - X3 - C3 - X31 - C4 - X10–11 - C5 - X8 - C6 - X17–28 - C1′ - X24–25 - C2′ - X3 - C3′ - X35 - C4′ - X9–12 - C5′ - X8 - C6′
C1 - X24–37 - C2 - X3 - C3 - X35–42 - C4 - X8–12 - C5 - X8 - C6
C1 - X26 - C2 - X3 - C3 - X39–40 - C4 - X10 - C5 - X8 - C6 and C1 - X23–30 - C3 - X38 - C4 - X18 - C6
C1a - X11–13 - C1b - C1c - C1 - X21–48 - C2 - X3 - C3 - X43 - C4 - X20–33 - C5 - X9 - C6 - X8 - C6a - X10 - C6b - X9 - C6c
C1 - X28–29 - C2 - X3 - C3 - X31 - C4 - X10–11 - C5 - X8 - C6 - X17–28 - C1′ - X24–25 - C2′ - X3 - C3′ - X35 - C4′ - X9–12 - C5′ - X8 - C6′
The average amino acid identity was higher in the dimer subfamily for both species, but the number of members in each subfamily may have inflated these values, since larger families show lower average identity (Table 2). Classic OBPs have the expected standard pattern of six conserved cysteine residues, whereas dimer OBPs have twelve cysteine residues, a pattern formed by the junction of two consecutive minus-C OBP domains . The dimer subfamily was identified in A. fraterculus, A. obliqua and D. melanogaster, though it has not been described for C. capitata yet. As in D. melanogaster, A. fraterculus and A. obliqua minus-C OBPs were divided in two lineages, one in which the members retain six conserved cysteine residues (a clade that includes OBP8a, OBP99c and OBP99d), and another one (that includes OBP83g, OBP99a and OBP99b) in which its members show only four conserved cysteine residues (Fig. 1). The reduced number of residues in some members is caused by the loss of the conserved pair of cysteines C2 and C5, which forms a disulphide bond . That happened after the minus subfamily diverged from the classic subfamily and might have a functional relevance, perhaps generating a more flexible structure [15, 19].
A. fraterculus and A. obliqua plus-C OBPs showed three additional conserved cysteine residues before C1 (referred to as C1a, C1b and C1c) and three others after C6 (C6a, C6b and C6c), making for a total of 12 conserved cysteine residues. Three sequences associated with OBP50a did not have the cysteine C6c. Similar to what was described for Drosophila, we found a conserved hydrophobic proline after the cysteine C6a in Anastrepha’s plus-C members, as well as nine residues between C5 and C6, instead of eight in the other subfamilies . Recently, Siciliano et al.  reported a plus-C OBP in C. capitata, related to DmelOBP49a, but in their phylogenetic tree, this OBP was grouped with classic OBPs. OBPs in the plus-C subfamily, such as the ones described here, generally show a substantial increase in length size compared with other OBPs, mostly because of the extended C-termini [5, 16]. The plus-C subfamily was also described in the pea aphid Acyrthosiphon pisum , and the high support in the tree indicates that their origin probably precedes the divergence of Diptera.
Even though the overall number of OBP genes seems to be equivalent among tephritids, we point out that species of Anastrepha have a larger OBPs repertoire than the other tephritids studied: 23 and 24 OBP genes, respectively. In comparison, 17 different OBPs were described for C. capitata, the closest species of A. fraterculus and A. obliqua studied, 13 of which corresponded to the classic OBP subfamily, three to the minus-C subfamily and one to the plus-C subfamily . We also found a larger number of plus-C OBPs when compared to other tephritids. We detected five apiece for A. fraterculus and A. obliqua, whereas two were described in R. pomonella and only one in R. suavis, B. dorsalis and C. capitata. The lack of genomic sequences for Anastrepha makes the analyses here not complete, since we may only rely on levels of divergence and phylogenetic expectations that could be verified solely upon investigation of the evolutionary patterns found in the species’ genome. It is possible that upon the availability of genome sequences, we may detect a wider range and number of OBP genes.
We observed some instances in which more than one sequence in the A. fraterculus and A. obliqua OBP repertoire was associated with the same D. melanogaster OBP. In A. fraterculus, two different sequences were associated with the genes OBP50a, OBP56d and OBP56h, and in A. obliqua, with the genes OBP49a, OBP56d, OBP56h and OBP99d. These findings could be due to gene duplication events, though we cannot rule out intraspecific variation, because we used a pool of flies from a single population to make the cDNA libraries in each species. Analyzing the values of pairwise identity between the two copies of the same OBP of a species and comparing copies of the same OBP between species, different Anastrepha sequences associated with the same OBP seem to be due to intraspecific variation, with the exception of two cases: OBP56h and OBP49a.
We believe that the copies homologous to OBP56h are consequence of a duplication event that preceded the divergence of A. fraterculus and A. obliqua. AfraOBP56h-1 and AoblOBP56h-1 have a pairwise identity of 93.6 %, and AfraOBP56h-2 and AoblOBP56h-2 an identity of 96.3 %, whereas the pairwise identity between OBP56h sequences of the same species is 29.3 % in A. fraterculus and 33.3 % in A. obliqua, clearly indicating them to be paralogous. Moreover, the phylogenetic analysis of plus-C members reveals that AoblOBP49a-2 groups with the other OBP49a, including DmelOBP49a (Additional file 2). However, AoblOBP49a-2 is 37 amino acids shorter than AoblOBP49a-1 and AfraOBP49a and quite different. While the pairwise amino acid identity between AoblOBP49a-1 and AfraOBP49a is 95 %, it is only 31 % between the two copies of A. obliqua (AoblOBP49a-1 and AoblOBP49a-2), suggesting as well that AoblOBP49a-2 is not an ortholog of DmelOBP49a, such as AoblOBP49a-1 and AfraOBP49a but, rather, a paralog. The fact that we failed to find this copy in A. fraterculus may indicate that it diverged after the separation of the species, which would be extremely interesting, but considering the high levels of divergence found among these paralogs, it is more likely that they diverged before the separation of Anastrepha species, and did not show significant levels of expression in A. fraterculus, which is also relevant.
Analysis of positive selection in OBPs
Similar to all insect OBPs studied so far , Anastrepha OBPs share low sequence similarity, rendering evolutionary analyses difficult. Therefore, we did not use the alignment that included all subfamilies because we identified that there was saturation at synonymous positions in more divergent comparisons (data not shown), which might lead to a higher number of false positives in the branch-site test of positive selection . We performed independent evolutionary analyses on each of the four OBP subfamily identified (Table 2; Additional file 2), which did not show evidence of sequence saturation. This, unfortunately, prevented us from investigating evidence of positive selection amongst the four subfamilies, but allowed us to investigate for positive selection within any of such subfamily.
Positive selection detected in Anastrepha OBP genes
Positively selected sitesb
40, 67, 68, 145, 152
5, 34, 61, 65, 76, 94, 146, 150
4, 56, 59, 146, 152, 156, 164, 177, 181
We investigated whether amino acid changes in OBPs detected by the BEB analysis were associated with site-specific changes related to five amino acid properties: chemical composition, polarity, volume, isoeletric point and hydropathy (p < 0.05), because not all non-synonymous substitutions are alike. If non-synonymous substitutions result in similar properties, the maintenance of the structure and/or chemical function of the region bearing such changes would be considered as conserved. On the other hand, radical non-synonymous changes drastically change important physicochemical attributes. When these changes involve sites under positive selection, it is more likely that they would be those that promote structural and functional changes in a protein. Among the five non-synonymous changes in the duplicated copy of Anastrepha OBP56h-2 we found two amino acid positions (68 and 145) that are associated with radical changes in the physicochemical properties, one associated to changes in isoelectric point, and another to changes in chemical composition. The high number of amino acid substitutions, particularly leading to radical changes, may be an indication of the effect of positive selection driving this recently duplicated copy to a different odor affinity.
Anastrepha OBP57c and Anastrepha OBP50a also have one amino acid position each associated with radical changes in the physicochemical properties. A change in the amino acid position 150 of Anastrepha OBP57c alters its isoeletric point, whereas a change in amino acid 59 of Anastrepha OBP50a affects its volume. Changes in chemical composition and volume of OBPs are important because they may be related to conformational changes in these globular proteins. Altering the structure of α-helices may modify their flexibility and, therefore, modify the interacting motifs of the protein . Likewise, substitutions that affect the isoeletric point are important since they may be involved with changes in the proteins’ solubility, which, in the aqueous lymph of sensillas, may be related with a greater efficiency in transport activity.
Estimates for Ka, Ks and Ka/Ks ratio between A. fraterculus and A. obliqua putatively orthologous OBP genes
OBP genes compared
AfraOBP8a - AoblOBP8a
AfraOBP19a - AoblOBP19a
AfraOBP19b - AoblOBP19b
AfraOBP19c - AoblOBP19c
AfraOBP19d - AoblOBP19d
AfraOBP47b - AoblOBP47b
AfraOBP49a - AoblOBP49a-1
AfraOBP49a - AoblOBP49a-2
AfraOBP50a-1 - AoblOBP50a
AfraOBP50a-2 - AoblOBP50a
AfraOBP50e - AoblOBP50e
AfraOBP56d-1 - AoblOBP56d-1
AfraOBP56d-2 - AoblOBP56d-2
AfraOBP56d-1 - AoblOBP56d-2
AfraOBP56d-2 - AoblOBP56d-1
AfraOBP56h-1 - AoblOBP56h-1
AfraOBP56h-2 - AoblOBP56h-2
AfraOBP56h-1 - AoblOBP56h-2
AfraOBP56h-2 - AoblOBP56h-1
AfraOBP57c - AoblOBP57c
AfraOBP59a - AoblOBP59a
AfraOBP83cd - AoblOBP83cd
AfraOBP83ef - AoblOBP83ef
AfraOBP83g - AoblOBP83g
AfraOBP99a - AoblOBP99a
AfraOBP99b - AoblOBP99b
AfraOBP99c - AoblOBP99c
AfraOBP99d - AoblOBP99d-1
AfraOBP99d - AoblOBP99d-2
Positively selected sites in OBP 3-D structures
Positively selected sites experience a faster rate of amino acid replacement that is a consequence of mutations being fixed at a higher rate than expected by chance. Therefore, knowing the location of these amino acid sites in the protein’s tertiary structure may help understand how these changes would affect their function. Structurally, OBPs share a common fold with α-helices connected by loops and interlinked by disulphide bonds, but despite their structural homology, they are predicted to bear binding cavities of different shapes . Positive selection in OBPs may promote functional divergence on the binding specificities, which is important because such structural diversity may enable OBPs to recognize and bind to a wider range of organic molecules and naturally occurring odorants .
The tertiary structure of A. gambiae OBP4 (AgamOBP4)  was selected by Phyre2 as a reference for AoblOBP57c and AfraOBP57c. AgamOBP4 is a protein with 124 amino acid residues, in which three disulphide bridges link α1 to α3, the start of α3 to α6, and the middle of α6 to α5. In Anastrepha (Fig. 4c), the positively selected sites were found in α2, α3, α4 and α6, and loops between α1 and α2, and between α3 and α4. Notably the positively selected amino acids located in α3 and especially in α6 may be structurally very important because they are located near disulphide bonds, specially the last one that causes changes in the chemical composition.
A. gambiae plus-C OBP48, selected by Phyre2 as reference for Anastrepha OBP50a, has 172 residues and eight α-helices that are stabilized by six disulphide bridges . AgamOBP48 is structurally divided into three domains: the “core”, the “NC-term” and the “cap”. The central “core” domain shares a common architecture with classic OBPs and consists of four bridged α-helices (α1, α3, α5 and α7), bearing six conserved cysteines that form three disulphide bridges connecting α1 to α3, α3 to the start of α7, and the middle of α7 to α5. This domain also contains two non-bridged α-helices α2 and α4, and a loop of eight amino acids (50s loop). Sites under positive selection in Anastrepha OBP50a sequences are found in the “core” region, in α3 and α7, in the loop between α2 and α3, and in the loop between α3 and β-sheet1 (Fig. 4d). A seventh cysteine in AgamOBP48 is located at the top of α3, but it is not involved in a disulphide bond, since not all cysteine residues are necessarily involved in disulphide bonds . For instance, A. fraterculus and A. obliqua plus-C OBP50a lost the twelfth conserved cysteine residue in a positively selected change that may have relevant structural consequences, but this loss does not seem to affect their function.
The “NC-term” domain comprises the N-terminal residues 1–25, the C-terminal 150 s loop (148–155) and the following C-terminal α8 helix (156–172), forming three disulphide bonds. The two subunits orient themselves so that the “NC-term” domain of one monomer inserts into the center of the neighboring monomer, to form a compact homodimer. The N-terminal component of this domain is characterized by two additional cysteines located in adjacent positions before C1 , pattern also observed in all Anastrepha plus-C OBPs. Anastrepha OBP50a has sites under positive selection in this domain, two ones in α8 and another one in the loop (Fig. 4d). Finally, the “Cap” domain, in which we failed to find sites under positive selection in Anastrepha, encompasses helix α6 and the loop near amino acid 120, and this interaction between the helix and the loop is characterized by multiple hydrogen-bonding that apparently serves to stabilize protein structure .
AgamOBP48 shows a single binding site, although it is argued that given the symmetry of the dimer, ligand binding to two equivalent binding sites, named “NC-term” pockets, cannot be excluded . We focused only on the single binding site to compare with Anastrepha OBP50a. The binding cavity is formed by 23 amino acid residues in the helices α5, α6 and α7, and in the loops between helices α2 and α3, and between α7 and α8. We found one amino acid under positive selection in the position exactly correspondent to the amino acid that forms the binding cavity in α7 of AgamOBP48. The other amino acid in α7, as well as the amino acids in the loop between α2 and α3, and in the loop between α7 and α8 did not correspond to the same position to the binding cavity, but they are juxtaposed to the amino acids that form it.
Because we are only using computational inferences, we cannot be certain about the 3D structures inferred for Anastrepha OBP56h, OBP57c and OBP50a, nor if the Anastrepha OBP binding cavity is formed in the same site positions that in the 3D structures of Ae. aegypti and A. gambiae OBPs used as reference. However, the OBPs’ high conservation at the structural level, observed in all insects studied, provide us some basis for this extrapolation [21, 22, 78]. The exception to this structure conservation is a unique feature of mosquito OBPs, their C-terminal loop covers the binding cavity, forming a “lid” for the release of ligands , therefore we expect that the C-terminal in Anastrepha may be quite different, potentially altering our 3D inference for the C-terminal part. Only Anastrepha OBP50a showed one amino acid under positive selection that possibly correspond to amino acids that form the binding cavity for odors. However, even if no replacement occurs in amino acids that are directly involved with the binding function, radical non-synonymous substitutions placed both in the α-helix or even in the loops are important because they might alter the size and shape of the binding cavity. For instance, by modifying the position of the disulphide bonds, as demonstrated elsewhere [21, 80].
Although the interaction mechanisms between OBPs and the OBP/ligand complex with ORs are still not well understood, studies showed that after binding the ligand, some OBPs are induced to a conformational change. These pH-dependent conformational changes were associated with changes in binding affinity, and it was reported to be common to some OBPs in distinct insects [12, 79, 81–84]. The OBP/ligand complex releases the ligand after they reach specific odorant receptors . Therefore, amino acid changes that lead to conformational changes, such as those found in A. fraterculus and A. obliqua OBPs may also interfere with the interactions between OBPs and ORs, even when they are not directly involved with the binding sites.
OBPs’ molecular evolution in A. fraterculus and A. obliqua may reflect specific adaptation
Despite the similar number of OBPs and their sequences, there is significant difference between OBPs of A. obliqua and A. fraterculus that may reflect specific adaptation. A wind tunnel test revealed that A. fraterculus antennae were responsive to Myrtaceae extracts, which also affected its oviposition rate . On the other hand, adults of A. obliqua were attracted to Anacardiaceae ripe fruits [86, 87], so much so that at least nine volatile compounds from the Anacardiaceae Spondias mombin elicited antennal response from both sexes of this species . Therefore, despite being considered generalists , both species have their preferences and show some host specificity for oviposition and feeding. Species of Anastrepha in general show lekking behavior , which has been described for A. obliqua [89, 90], as well as A. fraterculus . In these leks, a number of males compete for space and display to have access to females. Several factors have been associated with success in the leks, chief among them their pheromones . A comparison revealed that even though A. fraterculus and A. obliqua shared common chemical compounds on their pheromones, they also showed several different compounds emitted by calling males . These ecological and reproductive differences may have been the driving force behind the rapid rates of evolution we identified amongst their OBP sequences, and suggest that the evolution of OBP genes may have had a significant impact in the evolution of species differences in this group.
In this study, we used transcriptome data to identify over 23 different OBP genes in A. fraterculus and A. obliqua, which is the largest and the most diverse number of OBPs yet reported for a Tephritidae. We found great similarity in amino acid and DNA sequences among orthologous OBPs in A. fraterculus and A. obliqua, which may be a reflection of their recent divergence or evolutionary conservatism. However, OBPs from A. fraterculus and A. obliqua showed a faster rate of evolution when comparing to other genes among these species, a higher Ka/Ks ratio and evidence of positive selection on at least four Anastrepha OBP genes: OBP56h-1, OBP56h-2, OBP57c and OBP50a. We also found four positively selected sites in which site-specific changes would radically change amino acid properties, and likely promote structural and functional changes. One amino acid under positive selection in OBP50a is located in the binding cavity according the putative 3D-structure inference, which is important because such change may promote functional divergence of the binding specificities, and enable this protein to recognize and bind a new odorant. The other changes that are not directly involved with the binding function may also be important because they may alter the size and shape of the binding cavity or the solubility of the whole molecule. Considering that, as was shown in other insects, few amino acid changes in OBPs may result in significant differences in olfactory responses, our results stress out the importance of OBPs for the evolution and divergence of A. fraterculus and A. obliqua.
The authors would like to thank Federico G. Hoffmann and Michael W. Vandewege and our colleagues at the Population Genetics and Evolution Lab for helpful suggestions, and the anonymous reviewers from Axios for their invaluable comments and suggestions to improve the quality of the paper.
This study was supported by Fundação de Amparo à Pesquisa do Estado de São Paulo (grant numbers: 2012/17160-8 to E.B.C. and 2010/20455-4 to R.A.B.).
Availability of data and materials
EBC and RAB conceived and designed the research project, performed the analysis and wrote the manuscript. Both authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Orr HA, Irving S. Complex epistasis and the genetic basis of hybrid sterility in the Drosophila pseudoobscura Bogota-USA hybridization. Genetics. 2001;158(3):1089–100.PubMedPubMed CentralGoogle Scholar
- Noor MAF, Feder JL. Speciation genetics: evolving approaches. Nat Rev Genet. 2006;7(11):851–61.View ArticlePubMedGoogle Scholar
- Rice WR. Sexually antagonistic genes: experimental evidence. Science. 1992;256(5062):1436–9.View ArticlePubMedGoogle Scholar
- Swanson WJ, Vacquier VD. The rapid evolution of reproductive proteins. Nat Rev Genet. 2002;3(2):137–44.View ArticlePubMedGoogle Scholar
- Forêt S, Maleszka R. Function and evolution of a gene family encoding odorant binding-like proteins in a social insect, the honey bee (Apis mellifera). Genome Res. 2006;16(11):1404–13.View ArticlePubMedPubMed CentralGoogle Scholar
- Gardiner A, Barker D, Butlin RK, Jordan WC, Ritchie MG. Drosophila chemoreceptor gene evolution: selection, specialization and genome size. Mol Ecol. 2008;17(7):1648–57.View ArticlePubMedGoogle Scholar
- Kosiol C, Vinař T, Fonseca RR, Hubisz MJ, Bustamante CD, Nielsen R, et al. Patterns of positive selection in six mammalian genomes. PLoS Genet. 2008;4(8):e1000144.View ArticlePubMedPubMed CentralGoogle Scholar
- Sánchez-Gracia A, Rozas J. Divergent evolution and molecular adaptation in the Drosophila odorant-binding protein family: inferences from sequence variation at the OS-E and OS-F genes. BMC Evol Biol. 2008;8:323.View ArticlePubMedPubMed CentralGoogle Scholar
- Benton R. On the ORigin of smell: odorant receptors in insects. Cell Mol Life Sci. 2006;63(14):1579–85.View ArticlePubMedGoogle Scholar
- Vogt RG, Riddiford LM. Pheromone binding and inactivation by moth antennae. Nature. 1981;293:161–3.View ArticlePubMedGoogle Scholar
- Pelosi P, Maida R. Odorant-binding proteins in insects. Comp Biochem Physiol B Biochem Mol Biol. 1995;111(3):503–14.View ArticlePubMedGoogle Scholar
- Leite NR, Krogh R, Xu W, Ishida Y, Iulek J, Leal WS, et al. Structure of an odorant-binding protein from the mosquito Aedes aegypti suggests a binding pocket covered by a pH-sensitive ‘lid’. PLoS One. 2009;4(11):e8006.View ArticlePubMedPubMed CentralGoogle Scholar
- Lagarde A, Spinelli S, Tegoni M, He XL, Field L, Zhou J-J, et al. The crystal structure of odorant binding protein 7 from Anopheles gambiae exhibits an outstanding adaptability of its binding site. J Mol Biol. 2011;414(3):401–12.View ArticlePubMedGoogle Scholar
- Tsitsanou KE, Drakou CE, Thireou T, Gruber AV, Kythreoti G, Azem A, et al. Crystal and solution studies of the ‘Plus-C’ odorant-binding protein 48 from Anopheles gambiae: control of binding specificity through three-dimensional domain swapping. J Biol Chem. 2013;288(46):33427–38.View ArticlePubMedPubMed CentralGoogle Scholar
- Hekmat-Scafe DS, Scafe CR, McKinney AJ, Tanouye MA. Genome-wide analysis of the odorant-binding protein gene family in Drosophila melanogaster. Genome Res. 2002;12(9):1357–69.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhou J-J, Zhang GA, Huang W, Birkett MA, Field LM, Pickett JA, et al. Revisiting the odorant-binding protein LUSH of Drosophila melanogaster: evidence for odour recognition and discrimination. FEBS Lett. 2004;558(1–3):23–6.View ArticlePubMedGoogle Scholar
- Vieira FG, Forêt S, He XL, Rozas J, Field LM, et al. Unique features of odorant-binding proteins of the parasitoid wasp Nasonia vitripennis revealed by genome annotation and comparative analyses. PLoS One. 2012;7(8):e43034.View ArticlePubMedPubMed CentralGoogle Scholar
- Valenzuela JG, Charlab R, Gonzalez EC, Miranda-Santos IKF, Marinotti O, Francischetti IMB, et al. The D7 family of salivary proteins in blood sucking Diptera. Insect Mol Biol. 2002;11(2):149–55.View ArticlePubMedGoogle Scholar
- Vieira FG, Sánchez-Gracia A, Rozas J. Comparative genomic analysis of the odorant-binding protein family in 12 Drosophila genomes: purifying selection and birth-and-death evolution. Genome Biol. 2007;8(11):R235.View ArticlePubMedPubMed CentralGoogle Scholar
- Ohta S, Seto Y, Tamura K, Ishikawa Y, Matsuo T. Identification of odorant-binding protein genes expressed in the antennae and the legs of the onion fly, Delia antiqua (Diptera: Anthomyiidae). Appl Entomol Zool. 2014;49(1):89–95.View ArticleGoogle Scholar
- Tegoni M, Campanacci V, Cambillau C. Structural aspects of sexual attraction and chemical communication in insects. Trends Biochem Sci. 2004;29(5):257–64.View ArticlePubMedGoogle Scholar
- Galindo K, Smith DP. A large family of divergent Drosophila odorant-binding proteins expressed in gustatory and olfactory sensilla. Genetics. 2001;159(3):1059–72.PubMedPubMed CentralGoogle Scholar
- Nei M, Rooney AP. Concerted and birth-and-death evolution of multigene families. Annu Rev Genet. 2005;39(1):121–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Sánchez-Gracia A, Aguadé M, Rozas J. Patterns of nucleotide polymorphism and divergence in the odorant-binding protein genes OS-E and OS-F: analysis in the melanogaster species subgroup of Drosophila. Genetics. 2003;165(3):1279–88.PubMedPubMed CentralGoogle Scholar
- Perre P, Jorge LR, Lewinsohn TM, Zucchi RA. Morphometric differentiation of fruit fly pest species of the Anastrepha fraterculus group (Diptera: Tephritidae). Ann Entomol Soc Am. 2014;107(2):490–5.View ArticleGoogle Scholar
- Selivon D, Perondini ALP, Morgante JS. A genetic–morphological characterization of two cryptic species of the Anastrepha fraterculus complex (Diptera: Tephritidae). Ann Entomol Soc Am. 2005;98(3):367–81.View ArticleGoogle Scholar
- Pereira-Rêgo DRG, Jahnke SM, Redaelli LR, Schaffer N. Morfometria de Anastrepha fraterculus (Wied.) (Diptera: Tephritidae) relacionada a hospedeiros nativos, Myrtaceae. Arq Inst Biol. 2011;78:3743.Google Scholar
- Camargo CA, Odell E, Jirón LF. Interspecific interactions and host preference of Anastrepha obliqua and Ceratitis capitata (Diptera: Tephritidae), two pests of mango in Central America. Fla Entomol. 1996;79(2):266–8.View ArticleGoogle Scholar
- Matsuo T, Sugaya S, Yasukawa J, Aigaki T, Fuyama Y. Odorant-binding proteins OBP57d and OBP57e affect taste perception and host-plant preference in Drosophila sechellia. PLoS Biol. 2007;5(5):e118.View ArticlePubMedPubMed CentralGoogle Scholar
- Laughlin JD, Ha TS, Jones DNM, Smith DP. Activation of pheromone-sensitive neurons is mediated by conformational activation of pheromone-binding protein. Cell. 2008;133(7):1255–65.View ArticlePubMedPubMed CentralGoogle Scholar
- Rezende VB, Congrains CC, Lima ALA, Campanini EB, Nakamura AM, Oliveira JL, et al. Head transcriptomes of two closely related species of fruit flies of the Anastrepha fraterculus group reveals divergent genes in species with extensive gene flow. G3 Genes Genom Genet. 2016. doi:10.1534/g3.116.030486.Google Scholar
- Chomczynski P, Mackey K. Short technical reports. Modification of the TRI reagent procedure for isolation of RNA from polysaccharide- and proteoglycan-rich sources. BioTechniques. 1995;19(6):942–5.PubMedGoogle Scholar
- SeqyClean. https://github.com/ibest/seqyclean. Accessed 14 Jan 2014
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-Length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Evol. 1990;215:403–10.Google Scholar
- ORF finder. https://www.ncbi.nlm.nih.gov/orffinder/. Accessed 10 Feb 2014.
- Hiller K, Grote A, Scheer M, Münch R, Jahn D. PrediSi: prediction of signal peptides and their cleavage positions. Nucleic Acid Res. 2004;32:375–9.View ArticleGoogle Scholar
- Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30(4):772–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Hall TA. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999;41:95–8.Google Scholar
- Smith JM. Analyzing the mosaic structure of genes. J Mol Evol. 1992;34:126–9.PubMedGoogle Scholar
- Sawyer S. Statistical tests for detecting gene conversion. Mol Biol Evol. 1989;6(5):526–38.PubMedGoogle Scholar
- Martin D, Rybicki E. RDP: detection of recombination amongst aligned sequences. Bioinformatics. 2000;16(6):562–3.View ArticlePubMedGoogle Scholar
- Martin DP, Lemey P, Lott M, Moulton V, Posada D, Lefeuvre P. RDP3: a flexible and fast computer program for analyzing recombination. Bioinformatics. 2010;26(19):2462–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Posada D. jModelTest: phylogenetic model averaging. Mol Biol Evol. 2008;25(7):1253–6.View ArticlePubMedGoogle Scholar
- Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59(3):307–21.View ArticlePubMedGoogle Scholar
- Zheng W, Peng W, Zhu C, Zhang Q, Saccone G, Zhang H. Identification and expression profile analysis of odorant binding proteins in the oriental fruit fly Bactrocera dorsalis. Int J Mol Sci. 2013;14(7):14936–49.View ArticlePubMedPubMed CentralGoogle Scholar
- Schwarz D, Robertson HM, Feder JL, Varala K, Hudson ME, Ragland GJ, et al. Sympatric ecological speciation meets pyrosequencing: sampling the transcriptome of the apple maggot Rhagoletis pomonella. BMC Genomics. 2009;10(1):633.View ArticlePubMedPubMed CentralGoogle Scholar
- Ramsdell KMM, Lyons-Sobaski SA, Robertson HM, Walden KKO, Feder JL, Wanner K, et al. Expressed sequence tags from cephalic chemosensory organs of the northern walnut husk fly, Rhagoletis suavis, including a putative canonical odorant receptor. J Insect Sci. 2010;10:51.View ArticlePubMedPubMed CentralGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23(21):2947–8.View ArticlePubMedGoogle Scholar
- Zhou J-J, He XL, Pickett JA, Field LM. Identification of odorant-binding proteins of the yellow fever mosquito Aedes aegypti: genome annotation and comparative analyses. Insect Mol Biol. 2008;17(2):147–63.View ArticlePubMedGoogle Scholar
- Gong DP, Zhang HJ, Zhao P, Xia QY, Xiang ZH. The odorant binding protein gene family from the genome of silkworm, Bombyx mori. BMC Genomics. 2009;10(1):332.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhou J-J, Vieira FG, He XL, Smadja C, Liu R, Rozas J, et al. Genome annotation and comparative analyses of the odorant-binding proteins and chemosensory proteins in the pea aphid Acyrthosiphon pisum. Insect Mol Biol. 2010;19(Suppl2):113–22.View ArticlePubMedGoogle Scholar
- Xia X. DAMBE5: a comprehensive software package for data analysis in molecular biology and evolution. Mol Biol Evol. 2013;30(7):1720–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang J. Evaluation of an improved branch-site likelihood method for detecting positive selection at the molecular level. Mol Biol Evol. 2005;22(12):2472–9.View ArticlePubMedGoogle Scholar
- Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007;24(8):1586–91.View ArticlePubMedGoogle Scholar
- Anisimova M, Bielawsky JP, Yang Z. Accuracy and power of the likelihood ratio test in detecting adaptive molecular evolution. Mol Biol Evol. 2001;18:1585–92.View ArticlePubMedGoogle Scholar
- Zhang Z, Li J, Zhao X-Q, Wang J, Wong GK-S, Yu J. KaKs_Calculator: calculating Ka and Ks through model selection and model averaging. Genomics Proteomics Bioinformatics. 2006;4(4):259–63.View ArticlePubMedGoogle Scholar
- Posada D. Using MODELTEST and PAUP* to select a model of nucleotide substitution. In: Baxevanis AD, Davison DB, Page RDM, Petsko GA, Stein LD, Stormo GD, editors. Current protocols in bioinformatics. USA: John Wiley & Sons, Inc; 2003. p. 6.5.1–6.5.14.Google Scholar
- Swanson WJ, Wong A, Wolfner MF, Aquadro CF. Evolutionary expressed sequence tag analysis of Drosophila female reproductive tracts identifies genes subjected to positive selection. Genetics. 2004;168(3):457–65.View ArticleGoogle Scholar
- Conant GC, Wagner GP, Stadler PF. Modeling amino acid substitution patterns in orthologous and paralogous genes. Mol Phylogenet Evol. 2007;42(2):298–307.View ArticlePubMedGoogle Scholar
- Delport W, Poon AF, Frost SDW, Kosakovsky-Pond SL. Datamonkey 2010: a suite of phylogenetic analysis tools for evolutionary biology. Bioinformatics. 2010;26(19):2455–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Kelley LA, Sternberg MJE. Protein structure prediction on the web: a case study using the phyre server. Nat Protoc. 2009;4(3):363–71.View ArticlePubMedGoogle Scholar
- Davrazou F, Dong E, Murphy EJ, Jones DNM. Conformational ordering plays a key role in regulating heterodimeric interactions between odorant binding proteins from Anopheles gambiae. http://www.rcsb.org/pdb/explore.do?structureId=3q8i. Accessed 17 Sep 2014
- Zhu J-Y, Zhang L-F, Ze S-Z, Wang D-W, Yang B. Identification and tissue distribution of odorant binding protein genes in the beet armyworm, Spodoptera exigua. J Insect Physiol. 2013;59(7):722–8.View ArticlePubMedGoogle Scholar
- Manoharan M, Chong MNF, Vaïtinadapoulé A, Frumence E, Sowdhamini R, Offmann B. Comparative genomics of odorant binding proteins in Anopheles gambiae, Aedes aegypti, and Culex quinquefasciatus. Genome Biol Evol. 2013;5(1):163–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Kirkness EF, Haas BJ, Sun W, Braig HR, Perotti MA, Clark JM, et al. Genome sequences of the human body louse and its primary endosymbiont provide insights into the permanent parasitic lifestyle. Proc Natl Acad Sci U S A. 2010;107(27):12168–73.View ArticlePubMedPubMed CentralGoogle Scholar
- Vieira FG, Rozas J. 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. 2011;3:476–90.View ArticlePubMedPubMed CentralGoogle Scholar
- Qiao H, Tuccori E, He XL, Gazzano A, Field L, Zhou JJ, et al. Discrimination of alarm pheromone (E)-beta-farnesene by Aphid odorant-binding proteins. Insect Biochem Molec. 2009;39(5–6):414–9.View ArticleGoogle Scholar
- Wang S-Y, Gu S-H, Han L, Guo Y-Y, Zhou J-J, Zhang Y-J. Specific involvement of two amino acid residues in cis-nerolidol binding to odorant-binding protein 5 AlinOBP5 in the alfalfa plant bug, Adelphocoris lineolatus (Goeze). Insect Mol Biol. 2013;22(2):172–82.View ArticlePubMedGoogle Scholar
- Ahmed T, Zhang T-T, Wang Z-Y, He K-L, Bai S-X. Three amino acid residues bind corn odorants to McinOBP1 in the polyembryonic endoparasitoid of Macrocentrus cingulum brischke. PLoS One. 2014;9(4):e93501.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang P, Lyman RF, Shabalina SA, Mackay TFC, Anholt RRH. Association of polymorphisms in odorant-binding protein genes with variation in olfactory response to benzaldehyde in Drosophila. Genetics. 2007;177(3):1655–65.View ArticlePubMedPubMed CentralGoogle Scholar
- Anisimova M, Nielsen R, Yang Z. Effect of recombination on the accuracy of the likelihood method for detecting positive selection at amino acid sites. Genetics. 2003;164:1229–36.PubMedPubMed CentralGoogle Scholar
- Gotzek D, Robertson HM, 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. 2011;6(1):e16289.View ArticlePubMedPubMed CentralGoogle Scholar
- Siciliano P, Scolari F, Gomulski LM, Falchetto M, Manni M, Gabrieli P, Field LM, Zhou J-J, Gasperi G, Malacrida AR. Sniffing out chemosensory genes from the Mediterranean fruit fly, Ceratitis capitata. PLoS One. 2014;9(1):e85523.View ArticlePubMedPubMed CentralGoogle Scholar
- Fan J, Francis F, Liu Y, Chen JL, Cheng DF. An overview of odorant-binding protein functions in insect peripheral olfactory reception. Genet Mol Res. 2011;10(4):3056–69.View ArticlePubMedGoogle Scholar
- Fletcher W, Yang Z. The effect of insertions, deletions, and alignment errors on the branch-site test of positive selection. Mol Biol Evol. 2010;27(10):2257–67.View ArticlePubMedGoogle Scholar
- Kryazhimskiy S, Plotkin JB. The population genetics of dN/dS. Edited by Takashi Gojobori. PLoS Genet. 2008;4(12):e1000304.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhou J-J, Robertson G, He X, Dufour S, Hooper AM, Pickett JA, Keep NH, Field LM. Characterisation of Bombyx mori odorant-binding proteins reveals that a general odorant-binding protein discriminates between sex pheromone components. J Mol Biol. 2009;389:529–45.View ArticlePubMedGoogle Scholar
- Wogulis M, Morgan T, Ishida Y, Leal WS, Wilson DK. The crystal structure of an odorant binding protein from Anopheles gambiae: evidence for a common ligand release mechanism. Biochem Bioph Res Co. 2006;339(1):157–64.View ArticleGoogle Scholar
- Lartigue A. The crystal structure of a cockroach pheromone-binding protein suggests a new ligand binding and release mechanism. J Biol Chem. 2003;278(32):30213–8.View ArticlePubMedGoogle Scholar
- Sandler BH, Nikonova L, Leal WS, Clardy J. Sexual attraction in the silkworm moth: structure of the pheromone-binding-protein-bombykol complex. Chem Biol. 2000;7(2):143–51.View ArticlePubMedGoogle Scholar
- Xu W, Leal WS. Molecular switches for pheromone release from a moth pheromone-binding protein. Biochem Biophys Res Commun. 2008;372(4):559–64.View ArticlePubMedGoogle Scholar
- Mao Y, Xu X, Xu W, Ishida Y, Leal WS, Ames JB, Clardy J. Crystal and solution structures of an Odorant-binding protein from the southern house mosquito complexed with an oviposition pheromone. Proc Natl Acad Sci U S A. 2010;107(44):19102–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Xu X, Xu W, Rayo J, Ishida Y, Leal WS, Ames JB. NMR structure of navel orangeworm moth pheromone-binding protein (AtraPBP1): implications for pH-sensitive pheromone detection. Biochemistry. 2010;49(7):1469–76.View ArticlePubMedPubMed CentralGoogle Scholar
- Gregorio PLF, Sant’Ana J, Redaelli LR. Percepção química e visual de Anastrepha fraterculus (Diptera, Tephritidae) em laboratório. Iheringia Sér Zool. 2010;100(2):128–32.View ArticleGoogle Scholar
- Cruz-López L, Malo EA, Toledo J, Virgen A, Mazo A, Rojasm JC. A new potential attractant for Anastrepha obliqua from Spondias mombin fruits. J Chem Ecol. 2006;32(2):351–65.View ArticlePubMedGoogle Scholar
- López-Guillén G, López LC, Malo EA, Rojas JC. Olfactory responses of Anastrepha obliqua (Diptera: Tephritidae) to volatiles emitted by calling males. Fla Entomol. 2011;94(4):874–81.View ArticleGoogle Scholar
- Morgante JS, Selivon D, Solferini VN, Matioli SR. Evolutionary patterns in specialist and generalist species of Anastrepha. In: Aluja M, Liedo P, editors. Fruit flies: biology and management. New York: Springer; 1993. p. 15–20.View ArticleGoogle Scholar
- Aluja M, Piñeiro J, Jácome I, Díaz-Fleischer F, Sivinski J. Behavior of flies in the genus Anastrepha (Trypetinae: Toxotripanini). In: Aluja M, Norrbom AL, editors. Fruit flies (Tephritidae): phylogeny and evolution of behavior. Washington: CRC Press; 2000. p. 375–408.Google Scholar
- Aluja M, Cabrera M, Guillen J, Celedonio H, Ayora F. Behaviour of Anastrepha ludens, A. obliqua and A. serpentina (Diptera:Tephritidae) on a wild mango tree (Mangifera indica) harbouring three McPhail traps. Int J Trop Insect Sci. 1989;10(03):309–18.View ArticleGoogle Scholar
- Segura DF, Petit-Marty N, Sciurano R, Vera T, Calcagno G, Allinghi A, Cendra PG, Cladera JL, Vilardi J. Lekking behavior of Anastrepha fraterculus (Diptera: Tephritidae). Fla Entomol. 2007;90(1):154–62.View ArticleGoogle Scholar