- Research article
- Open Access
The odorant receptor repertoire of teleost fish
BMC Genomics volume 6, Article number: 173 (2005)
Vertebrate odorant receptors comprise three types of G protein-coupled receptors: the OR, V1R and V2R receptors. The OR superfamily contains over 1,000 genes in some mammalian species, representing the largest gene superfamily in the mammalian genome.
To facilitate an informed analysis of OR gene phylogeny, we identified the complete set of 143 OR genes in the zebrafish genome, as well as the OR repertoires in two pufferfish species, fugu (44 genes) and tetraodon (42 genes). Although the genomes analyzed here contain fewer genes than in mammalian species, the teleost OR genes can be grouped into a larger number of major clades, representing greater overall OR diversity in the fish.
Based on the phylogeny of fish and mammalian repertoires, we propose a model for OR gene evolution in which different ancestral OR genes or gene families were selectively lost or expanded in different vertebrate lineages. In addition, our calculations of the ratios of non-synonymous to synonymous codon substitutions among more recently expanding OR subgroups in zebrafish implicate residues that may be involved in odorant binding.
The perception and discrimination of thousands of different odorants by the vertebrate olfactory system results from the activation of specific odorant receptors expressed by olfactory neurons in the nose. The first odorant receptors were identified in the rat  and belong to what is now referred to as the "OR" superfamily of odorant receptors . ORs exhibit a predicted seven transmembrane topology and sequence motifs characteristic of the A family (rhodopsin-like or Class I) of G protein-coupled receptors. Subsequent to the initial discovery of the OR superfamily of odorant receptors, two unrelated types of G protein-coupled receptors (GPCRs) were identified in the mammalian vomeronasal organ, the V1R receptors  and the V2R receptors [4–7]. The vomeronasal V1R and V2R receptors are thought to subserve signaling to pheromonal compounds .
The OR gene superfamily is the largest multigene superfamily described in mammalian genomes. The completion of both the Celera and public consortium versions of the mouse genome confirmed the existence of about 1068 potential intact OR genes (comprising at least 228 subfamilies) and 334 pseudogenes [8, 9]. In humans, there are ~340 intact OR genes and ~300 pseudogenes [10–12]. By way of contrast, molecular cloning and genomic DNA blot hybridizations in fish species suggest an OR repertoire size approximately five- to ten-fold smaller than that of mammalian species [13, 14].
An understanding of how vertebrate olfactory receptor repertoires evolved can be gained from comparing the properties and organization of genes from divergent vertebrate species. In this regard, the zebrafish, Danio rerio, provides a useful model for comparative genomics studies. Recent studies have demonstrated that the zebrafish genome encodes only 1 V1R-like receptor  (T.S.A and J.N., unpublished results) and ~60 olfactory C-family (Class III) GPCRs related to the mammalian V2R family (T.S.A., P. Luu, E. VanName, and J.N., manuscript in preparation); one fish olfactory C-family receptor has been shown to be activated by amino acids [16, 17], which are potent odorants for fish.
In the OR superfamily, approximately 28 genes and 5 pseudogenes were identified previously in zebrafish using PCR and homology-based techniques [14, 18–22]. Although a number of phylogenetic reconstructions have been made [8, 9, 18, 23–28], a more accurate view of the OR superfamily's evolutionary history would be facilitated by comparisons between genomic datasets that include a more complete representation of member genes from each species (see also ).
For the present study, we carried out genome data mining on the zebrafish genome sequence provided by the Sanger Institute Danio rerio Sequencing Project and found 143 potentially intact genes belonging to the zebrafish OR superfamily. We find that despite the limited size of the zebrafish OR repertoire, it comprises eight diverse OR families, with family members sharing on average ~40% amino acid identity. In addition, OR genes from two pufferfish species – fugu and tetraodon – can be grouped into six families which overlap with the zebrafish gene families. Analysis of the ratio of possible non-synonymous to synonymous codon substitutions suggests that OR genes in general are under negative or purifying selection; only a small number of residues within the transmembrane domains – the likely sites of odorant binding – appear to have undergone positive selection. Based on these findings we propose a model for the evolution of the vertebrate OR repertoire.
Results and discussion
Prediction of zebrafish OR genes
The third (Zv3) and fourth (Zv4) draft zebrafish genome assemblies ftp://ftp.ensembl.org/pub/assembly/zebrafish/ of whole genome shotgun sequence (5.7 × coverage) were searched for OR gene sequences using a modification of the method described for identifying OR sequences from the mouse genome .
The protein coding sequences of the vast majority of known OR genes characterized to date are uninterrupted by introns, which obviates the need for splice site prediction in the identification of most OR genes. Our gene prediction strategy was to combine a low-threshold BLAST search with profile Hidden Markov Model- (HMM) based gene prediction with the program Genewise ftp://ftp.ebi.ac.uk/pub/software/unix/wise2/. The Genewise results were post-processed using custom Perl scripts to generate complete ORFs. This process was repeated in an iterative fashion, as follows. The zebrafish genome assembly was subjected to TBLASTN search with a representative set of known zebrafish OR s (<50% percent identity among members of this set). The gene prediction program Genewise was then run on the genomic sequences surrounding each unique BLAST hit using a profile hidden Markov model (HMM) of the OR superfamily (see Methods for details). In each round, the newly predicted OR genes were added as queries for the next BLAST search. They were also aligned to previous members, and a new profile HMM was constructed for use in the next round of gene prediction.
Initial query sequences included members of each of the OR9, OR2, OR5, OR4, and OR13 subfamilies described previously . Predicted genes were subjected to a set of criteria for inclusion in the final set of OR genes. First, all coding sequences were required to be longer than 700 base pairs in length and show highest sequence similarity to previously characterized OR sequences. To be considered full-length, the deduced amino acid sequence was required to be greater than 275 amino acids, contain seven predicted transmembrane domains and exhibit the presence of a conserved N-linked glycosylation site with the pattern N-X- [TS]-X (where X is any amino acid residue except for proline) at the N-terminus. Sequences failing to meet these criteria and those lacking start or stop codons due to assembly gaps were considered to be partial genes. FASTY comparison to known OR peptide sequences was used to generate conceptual translations of potential pseudogenes and to determine the number of disruptions (frame shifts, early stop codons and in-frame deletions). To compensate for disruptions in OR coding sequences due to possible sequencing and/or assembly errors, we adopted a classification system previously proposed for mining of the mouse genome . In the present study, a gene was considered intact if it had a complete coding sequence with up to one disruption, or a pseudogene if it was either a partial sequence with one or more disruptions or a full-length sequence with two or more disruptions.
Based on these criteria, our search identified 143 intact OR genes (136 with no disruptions), 7 partial genes, 10 pseudogenes greater than 700 bp in length, and 15 gene fragments shorter than 700 bp (these shorter fragments were excluded from further analysis) [see Additional file 2 and Additional file 10]. Thus, we believe our total OR gene count is a conservative estimate of the true size of the OR repertoire, with between 78% (136/175) and 86% ([143+7]/175) of identifiable OR sequences consisting of potentially functional OR genes.
How complete is the predicted OR gene repertoire? To address this question, we extracted 65 zebrafish OR sequences (28 non-redundant published OR genes or cDNAs, 4 unpublished full-length cDNAs, and 33 ESTs) from Genbank and determined whether or not they were represented in the set of predicted genes by aligning them to OR genomic sequence. The 33 ESTs were first consolidated into 23 clusters based on overlapping sequences. Eight of these corresponded to known genes, while 15 clusters were novel, bringing the non-redundant set of Genbank genes to 47. Forty-five out of 47 of these sequences were identified in our search [see Additional file 10]. The two that we were unable to identify, OR115-15 and OR128-14 (an unpublished cDNA that was identified in the more recent Zv5 genome assembly and an EST, respectively) belong to two of the largest subfamilies in the repertoire. Thus, we are confident that the repertoire of odorant receptors described here includes nearly all of the OR genes encoded in the zebrafish genome. However, because some gaps remain in the assembly near several small receptor gene clusters, we expect that as the genome sequence is refined, a few additional receptor genes will be found.
OR nomenclature and classification
OR families are typically defined as monophyletic groups with members that share greater than 40% amino acid identity, whereas subfamily members share greater than 60% amino acid identity . Using these operational definitions, we classified the zebrafish OR genes into families and subfamilies by reconstructing their phylogeny by neighbor-joining with 1000 bootstrap replicates. Clades of OR genes with less than ~40% and less than ~60% inter-branch amino acid identity were used to group genes into distinct families and subfamilies, respectively. The average percent identity between families is approximately 25% while the maximum observed percent identity between any two ORs of different families is 39%.
To unify the naming for zebrafish OR genes, we propose a revised nomenclature based on the following rationale. Both newly predicted and previously described OR genes were named (or re-named) according to subfamily membership. Subfamilies were numbered sequentially starting at the number 101 (to avoid confusion with previous zebrafish OR nomenclature) in a depth-first traversal of the phylogenetic tree. Within subfamilies, OR s were numbered sequentially according to genomic position, if known. The new nomenclature showing subfamily membership and correspondence to previously identified zebrafish OR genes are shown in Table S1 [see Additional file 10].
Genomic distribution of zebrafish OR genes
Previous studies have demonstrated that OR genes are clustered in vertebrate genomes [8, 9, 18, 31, 32]. In mammalian genomes, OR genes are distributed widely, residing on 18 chromosomes in the mouse  and 21 chromosomes in humans [10, 11]. From the zebrafish Zv3 and Zv4 assemblies, we found that 119 of the identified zebrafish OR genes are distributed in five major clusters containing between 14 and 31 genes each. There are two clusters on chromosome 15, two on chromosome 21, one on chromosome 10, several small clusters on chromosomes 8, 14 and 17, and in a few cases, genes exist as singletons (Figure 1) [see Additional file 10]. Subfamilies are largely contiguous (see below) and subfamily members usually share the same transcriptional orientation, suggesting tandem duplication as a mechanism of expansion within a subfamily [18, 33]. We were able to assign genomic locations for ~80% of the OR genes we identified (29 remain unassigned).
Phylogeny of zebrafish OR genes
Using a neighbor joining algorithm (see Methods), we constructed a phylogenetic tree of the 143 intact OR genes and 4 full-length pseudogenes identified in the two zebrafish genome assemblies, using the zebrafish melanocortin receptors as an outgroup (Figure 2a). Based on this analysis and the criteria set forth above, zebrafish ORs could be classified into 8 families (≥40% intra-family sequence identity) and 40 subfamilies (≥60% intra-subfamily sequence identity). An additional thirteen sequences comprising 7 partial genes and 6 pseudogene fragments were subsequently assigned to subfamilies based on their sequence similarities and additional phylogenetic analyses (data not shown). Most of the gene families contain between 12 and 40 genes each; the two smallest families, Family A and Family B, contain 6 and 1 genes each, respectively. The intra-subfamily identity threshold was lowered for three subfamilies, OR102, OR115 and OR125, to generate monophyletic clades [see Additional file 11]. High bootstrap support (Figure 2a) justify these classifications, with all subfamilies exhibiting bootstrap scores of 100%.
The topology of the phylogenetic tree shown in Figure 2a is supported by three additional lines of evidence. First, we calculated all possible pairwise identities both within and between different groups of OR sequences. With the three exceptions noted above for subfamilies OR102, OR115, and OR125, the minimum percent identity within each subfamily is ≥ 62% [see Additional file 11]. Importantly, the maximum inter-subfamily identity is 44% [see Additional file 12] and the maximum interfamily identity is 38% [see Additional file 13]; both of these values are well below the ≥ 62% identity typically observed between members within a given subfamily. Thus, the sorting of ORs by neighbor joining analysis into distinct families and subfamilies is supported by an analysis based on all possible pairwise identities. Second, highly related OR genes are tightly clustered in the zebrafish genome, with the members of a given subfamily residing adjacent to one another, uninterrupted by more distantly related genes [18, 33]. In the present analysis, we found that the assignment of OR subfamilies by neighbor joining analysis indeed is consistent with this genomic organization; out of 23 multigene subfamilies, members from only five (OR111, OR113, OR126, OR128, and OR133) are found in genomic clusters interrupted by genes from other subfamilies (Figure 1) [see Additional file 10]. Third, we constructed a phylogenetic tree using a maximum likelihood algorithm; at both the family and subfamily levels, maximum likelihood analysis yields tree topologies comparable to those derived by neighbor joining [see Additional file 5].
Comparison of fish and mammalian OR repertoires
To gain additional insight into how the OR gene superfamily evolved in vertebrates, the zebrafish ORs were aligned to additional sets of vertebrate OR sequences. OR genes were predicted from the genome sequences of two pufferfish species, fugu (Takifugu rubripes) http://genome.jgi-psf.org/fugu6/ and tetraodon (Tetraodon nigroviridis) http://www.genoscope.cns.fr/externe/tetranew/, using methods identical to those used for finding zebrafish OR s. Forty-four (3 with one disruption) and 42 (6 with one disruption) intact genes were found in fugu [27, 29] and tetraodon, respectively [see Additional file 3 and Additional file 4]. As the genome sequence data for these two species represent ~95% and ~92% coverage, respectively, we expect that the OR genes identified here comprise the majority of each species' OR repertoire. Thus, the OR repertoires of both pufferfish species appear to be only ~one-third the size of the zebrafish repertoire. A summary of the OR genes identified in zebrafish, fugu, and tetraodon genomes is provided in Table 2.
Nine-hundred-thirty-five mouse OR sequences [8, 9] were either downloaded from Genbank (864 genes) or extracted from MGSCv3 using published coordinates (71 genes) . Phylogenetic trees were computed for ORs from zebrafish and mouse (Figure 2b) [see Additional file 6] and zebrafish, fugu and tetraodon (Figure 2c) [see Additional file 7]. The location of the melanocortin receptor branch represents the root of each tree.
Mouse ORs can be classified into two groups, Class I and Class II, each showing on average greater than 40% intra-group sequence identity . Based on their greater similarity to the limited number of fish OR genes identified prior to the present study, Class I genes from amphibians and mammals have been referred to as "fish-like" [8, 10, 24]. However, our analysis of the complete set of zebrafish OR genes indicates that this view cannot be generalized to the entire fish OR repertoire. Mammalian Class I and Class II genes can in fact be grouped more closely with only two out of eight ~equidistantly-related zebrafish families; Class I genes show close similarity to only a small subset of zebrafish OR genes (OR112-1, OR113-1, OR113-2 and OR114-1, which together comprise Family A), and one zebrafish gene (OR101-1, comprising the single member Family B) clusters together with mammalian Class II genes (Figure 2b). We base these conclusions on phylogenetic reconstructions as determined by neighbor joining (Figure 2b) and maximum likelihood [see Additional file 6], as well as on a separate calculation of average pairwise identities of genes between families (Table 1) [see Additional file 13]. In all cases, the alignment of mouse and zebrafish genes was gap-minimized and trimmed to remove N- and C-terminal tails [see Additional file 2]. Overall, mouse Class I exhibits similar average pairwise identity to the zebrafish families (27.3 ± 4.8% identity [mean ± standard deviation]; range: 17 – 32%) as mouse Class II (27.7 ± 5.5%; range: 18 – 38%); the difference in mean values is not significant in a two-tailed t-test (p = 0.89). Calculations comparing consensus sequences representing each family yielded similar results (data not shown).
A comparison of teleost OR genes further reveals that six of the eight zebrafish OR families overlap with pufferfish families; families B and G do not appear to be present in pufferfish (Figure 2c) [see Additional file 7]. In our phylogeny comparing zebrafish and pufferfish genes by neighbor joining, we find low bootstrap support for Family E (score = 47), likely reflecting this family's closer proximity to Family F in the multi-species tree as compared to the tree generated with zebrafish OR genes alone (however, Family E has high bootstrap support by maximum likelihood analysis [see Additional file 7]). Interestingly, in zebrafish the most divergent family (Family H) shows only 17–19% identity to other families (versus 25–34% interfamily identity amongst the other families; Table 1). The location of the outgroup melanocortin receptor between Family H and the other families supports the conclusion that this family is the result of a very ancient gene duplication event. Based on the degree of divergence from other OR gene families, it is possible that the genes comprising Family H may not in fact encode bona fide odorant receptors. However, the predicted zebrafish Family H receptors retain one of the highly conserved OR signature motifs (see below), and one member of this family (OR137-7) was previously identified as an EST from a zebrafish olfactory epithelium cDNA library [see Additional file 10]. In addition, when zebrafish Family H sequences were used in BLAST searches of both the non-redundant protein sequence database and the mouse genome sequence, previously identified OR sequences were identified as the closest hits (data not shown). Family H also forms a cluster distinct from non-OR GPCRs in a phylogenetic tree comprising mouse and zebrafish ORs together with a set of 199 non-OR Type I (rhodopsin class) mouse GPCRs [see Additional file 8]. Thus, for the present purposes we consider the Family H sequences operationally as OR genes. More generally, this phylogenetic reconstruction based on OR and non-OR GPCRs reveals that the ORs as a group are distinct from the other Type I GPCRS.
A similar phylogeny for vertebrate OR genes was recently described . This study placed zebrafish, fugu, Xenopus and chicken OR genes into groups roughly comparable to those described here in Figure 2, with Family A corresponding to these authors' Group β (which clusters closely with human Class I genes); the single zebrafish gene comprising Family B falling within Group γ /human Class II; Family C corresponding to Group ε; Families D and G corresponding to Group ζ (which is not a monophyletic clade); Families E and F corresponding to Group δ; and Family H corresponding to Group η. Two highly divergent groups (not identified or retained in our search) – termed κ and θ – were also described, although their identities as OR genes are unclear .
Conserved motifs in predicted OR protein sequences
Previous studies of vertebrate ORs have identified a number of conserved sequence motifs characteristic of these receptors [8, 12, 35]. These include the following: an N-linked glycosylation site NX [TS]X in the N-terminal domain; the motif MA [FY] [DE]RYVAIC located at the third transmembrane domain (TM3)/second intracellular loop (IC2) junction which is thought to interact with G-proteins (specifically Golf); three conserved cysteine residues in the second extracellular loop (EC2) thought to partake in disulfide bonding; and the motif KAFSTCXSH in IC3 containing an intracellular cysteine conserved in GPCRs and potential phosphorylation sites. We found that these motifs are conserved in all the zebrafish OR families, with the exception of Family H, in which only the MAYDRYVAIC motif is conserved. This sequence conservation is illustrated by a sequence logo generated from the alignment of predicted full-length zebrafish OR coding sequences (Figure 3a). In this representation, the relative frequency with which an amino acid appears at a given position is reflected by the height of its one-letter amino acid code in the logo, with the total height at a given position proportional to the level of sequence conservation. Interestingly, when compared to the sequence logo representing the alignment of mouse Class I and Class II ORs (Figure 3b), the zebrafish OR logo shows lower conservation amongst the predicted zebrafish receptor sequences (reflected by the more numerous and shorter letters at individual positions in the logo), revealing the greater diversity within the zebrafish vs. mouse OR superfamily (Table 1).
Adaptive evolution of OR genes
What evolutionary processes might explain the diversity of OR gene sequences? This diversity could be the result of genetic drift, with polymorphisms in the population being fixed at a rate consistent with the absence of selective pressure. Alternatively, the sequence diversity could reflect true functional divergence, perhaps in the ligand specificities of the encoded receptors. Considering the diversity of OR proteins encoded in vertebrate genomes and the even greater diversity of compounds detected by these receptors, ligand binding sites within ORs may be expected to be under positive selective pressure as organisms evolve new receptor proteins to recognize odorant compounds. As for other Type 1 (rhodopsin class) GPCR ligands, odorants are thought to bind to OR proteins in the plane of the membrane, in contact with residues in the transmembrane domains . Consistent with this notion, previous studies have demonstrated that odorant receptor genes have been subject to positive selective pressure, especially in the transmembrane domains thought to coordinate odorant binding [13, 36] (however, see ). We therefore attempted to pinpoint the precise codon sites – and thus the amino acid residues – that may have been subjected to positive selection during the evolution of the zebrafish OR superfamily. To this end, we used the relative frequency of non-synonymous vs. synonymous codon substitutions to assess the selective processes acting on these receptor genes . Where there is no positive or negative selection on a sequence, the number of non-synonymous changes relative to the number of possible non-synonymous changes (dN) would be equal to the number of synonymous changes relative to the number of possible synonymous changes (dS) – i.e., dN/dS = 1. Significant deviations of dN/dS from unity reflect selection on the sequence; a dN/dS ratio > 1 indicates that a region has undergone positive selection, whereas a dN/dS ratio < 1 indicates negative or "purifying" selection . For our analysis, we aligned 136 full-length zebrafish OR coding sequences containing no disruptions and calculated dN/dS ratios based on a gap-minimized alignment (Table 3 and Figure 4a). We found that when the OR coding sequence is partitioned broadly into transmembrane domains (TMs) 1–7 and non-transmembrane domains (excluding the N- and C-terminal tails), none of these regions exhibits positive selection. Rather, with average dN/dS ratios <1, these protein regions all appear to be under negative or purifying selection. Interestingly, TMs 1, 3, 4, 5, and 6 display significantly higher average dN/dS ratios than the combined intracellular and extracellular loops (p < 1 × 103; see Figure 4a). The observation that these transmembrane domains were in general under less purifying selection than other regions of the protein is consistent with the possibility that they may have adapted to bind different odorants. In contrast, TM2 and TM7 display significantly lower average dN/dS ratios compared to the complete coding sequence (p < 0.05). The apparently stronger negative selection on TM2 and TM7 (as compared to the other transmembrane regions) suggests that these transmembrane domains subserve a common – perhaps structural – role in these receptors.
We hypothesized that specific codon sites corresponding to odorant-binding residues might have been positively selected as coding sequences diverged after gene duplication events. Accordingly, identification of these sites would suggest possible ligand binding sites. We therefore performed a site-by-site analysis of dN/dS ratios based on the alignment of the set of 136 full-length intact coding sequences used above. To avoid a potentially high rate of false positives common with pooled site methods (see ), we used the more conservative likelihood individual site (IS) method  based on the original proposed IS approach . The phylogenetic relationships between sequences were determined and a substitution model was estimated from the data. The ancestral codon sequences at each node were then reconstructed and the dN and dS values were calculated for each codon site. In Figure 4b, the probability of being under positive or negative selection (dN/dS values different than dN/dS = 1.0) for each codon site is indicated on a snake plot of a representative OR amino acid sequence, OR124-3 (see also Table 3). By these criteria, only two sites within the transmembrane domains (one in TM3 and one in TM4) appear to have been subjected to positive selection, consistent with the notion that they may play a role in contacting ligands. Interestingly, two adjacent sites in the short third extracellular loop (very close to the top of TM6) also exhibit dN/dS ratios > 1. Overall, our characterization of dN/dS ratios reveals a striking paucity of sites exhibiting signs of positive selection, possibly reflecting the dominating influence of negative selection throughout the receptor coding region. Alternatively, since non-synonymous substitutions may occur only sporadically over evolutionary time, the signatures of less recent substitutions may no longer be detected by this analysis of the entire zebrafish OR family.
Evolution of the vertebrate OR gene repertoire
The characterization of the complete OR repertoires from both fish and mammalian species allows an informed analysis of OR gene evolution in the vertebrate lineage. One noteworthy feature of our phylogenetic reconstruction is the presence of a group of zebrafish and pufferfish subfamilies which together form a putative OR family (Family H) more divergent than the other families are to each other (Figure 2c and Table 1); based on a BLAST search of the mouse genome using representative zebrafish Family H sequences, this family is absent from the mouse. We hypothesize that the node between this branch of the tree and the others is the root representing the most ancient gene duplication event observable in the teleost lineage. This is supported by the placement at this node of the melanocortin receptor (outgroup) branch. In addition, when we aligned five OR sequences from lamprey [25, 42] to the teleost ORs, they formed two additional families on either side of the melanocortin receptor branch, one which is clearly an OR family (more similar to teleost OR families A-G than to H) and one which appears as an outgroup (equidistant from all teleost families A-H and more dubious as an OR family) [see Additional file 9]. Since the lamprey diverged before the teleost/tetrapodon split, these observations provide further support for this node as the root of the tree. It should be noted, however, that until the lamprey genome has been fully sequenced, we will only have a partial picture of the ancestral OR repertoire. We expect that characterization of the entire lamprey OR repertoire will shed light on more ancient evolutionary events.
From our analysis of mammalian, teleost and lamprey OR sequences, we propose the following model for OR gene evolution in vertebrates. OR genes in present-day vertebrates likely descended from eight ancestral OR genes (or gene families) that existed at the time of the split between ray-finned and lobe-finned fish (the ancestors of teleosts and tetrapods, respectively) approximately 450 million years ago (mya) . A phylogenetic reconstruction based on mouse and zebrafish ORs and 199 mouse non-OR GPCRs [see Additional file 8] indicates that the ORs form a group distinct from all other Type I GPCRs, possibly reflecting a very ancient duplication event(s) and/or rapid divergence of the ORs in the evolution of Type I GPCRs. Our estimate of ancestral OR gene number is based on the identification of 8 OR gene families in teleosts, two of which show somewhat higher similarity to the 2 OR gene families in mammals. The grouping of zebrafish and pufferfish OR genes into common families indicates that the gene duplication events that gave rise to the major OR families probably occurred prior to the speciation of teleosts. In addition, the greater similarities between zebrafish Family A and mouse Class I, and between zebrafish Family B and mouse Class II infer that the ancestral genes for these families existed before the tetrapodon/teleost split. Our model therefore suggests a history during which ancestral genes or gene families were selectively lost during the evolution of the different vertebrate lineages. Of the ancestral families, zebrafish retained 8 families, fugu and tetraodon retained 6 families, and mammals retained 2 families. It should be noted that the low bootstrap score (47) for Family E in the comparison of zebrafish and pufferfish OR genes (Figure 2c) raises the possibility that Families E and F (which are adjacent to each other in the teleost phylogenetic tree) may have arisen from a more recent duplication in the teleost lineage. Alternatively, the genes in these groups may have been subjected to gene conversion events, with the effect of homogenizing the sequences between these two families.
It is also possible that the 4–6 gene families unique to teleosts descended from Family A/Class I and/or Family B/Class II ancestral genes, after the tetrapodon/teleost split. Such a scenario seems unlikely, however, considering the roughly equivalent degree of divergence exhibited between 7 out of the 8 teleost gene families (including Families A and B). Moreover, amphibian and avian OR genes can be grouped into 6 out of the 8 identified OR families (Families A, B/Class II, C, E, F and H), further implicating the presence of common ancestral genes for these families prior to the tetrapodon/teleost split .
Mechanisms of gene or family loss in a particular vertebrate lineage may have involved a number of processes, for example, gene conversion, pseudogenization of all genes in a family, unequal crossover recombination events during meiosis, or larger chromosomal rearrangements. From the available data we cannot infer the precise order and rate of OR gene family expansion and contraction, or speciation events. Nonetheless, six of the retained OR gene families were subject to a substantial net expansion and diversification in zebrafish (and to a lesser degree in the pufferfish), while the other two ancestors gave rise to the present-day mammalian Class I and Class II ORs as well as a small number of zebrafish genes. We hypothesize that relaxed selective pressure on a subset of the ancestral tetrapodon OR repertoire led to the loss of major OR gene families in the mammalian lineage. The expansion within the two remaining gene families was likely driven by the adaptation to the terrestrial odorous environment. Thus, different selective pressures found in the aquatic and terrestrial environments led to different sizes and shapes of the OR repertoires of fish and mammals.
It is generally thought that the diversity of OR sequences – as represented in the number of receptor families – underlies the diversity of chemical structures or "odor space" that can be detected by an organism's olfactory system. Thus, with ~6–8 OR gene families retained over evolutionary time (vs. 2 in mammals), fish may be capable of detecting a larger diversity of chemical structures than mammals. However, the larger total number of OR sequences in mammals (~1,000 vs. ~100 in fish) presumably allows a finer discrimination amongst the compounds that are detected by the mammalian olfactory system.
Iterative data mining
Genome-wide searches of the third (Zv3) and fourth (Zv4) draft zebrafish genome assemblies ftp://ftp.ensembl.org/pub/assembly/zebrafish/ made available by the Sanger Center on Nov 27, 2003, and July 12, 2004, respectively, were performed several times using the predicted ORs from each previous round to increase our querying power. This iterative data-mining approach has been published for finding OR genes in the mouse genome . A detailed description of our protocol is provided in the Supplement [see Additional file 1].
Alignment and tree construction
For multiple alignments of OR genes, ClustalX 1.81  was used with default parameters and gaps were inspected manually and edited in xced http://www.biophys.kyoto-u.ac.jp/~katoh/programs/align/xced/ to ensure integrity of transmembrane domains and proper alignment of anchoring OR motifs. N- and C-terminal tails were trimmed for all alignments. The neighbor-joining algorithm as implemented by PFAAT http://pfaat.sourceforge.net/ was used to generate unrooted phylogenetic trees from these alignments using the BLOSUM 50 similarity matrix; positions with greater than 40% gaps were excluded. One thousand bootstraps were performed to assess the support at each tree node. Trees were visualized with unrooted http://pbil.univ-lyon1.fr/software/njplot.html. Maximum likelihood analysis was carried out using PHYML http://atgc.lirmm.fr/phyml/ on the same processed amino acid alignments described above. Bootstrap analysis with 100 replicates was carried out using the JTT model of amino acid substitution. The consensus tree including bootstrap support for each node was plotted for each dataset using either ATV http://www.genetics.wustl.edu/eddy/atv/ or unrooted.
The sequences used for comparison to the zebrafish OR genes were obtained from Genbank and included the set of intact MORs [Genbank: AY072961] – [Genbank: AY074256]  plus 71 newly identified OR genes extracted from MGSCv3 using coordinates from the online supplement to , 5 full-length lamprey OR receptors [Genbank: AAC82383, Genbank: AAC82384, Genbank: AAC82385, Genbank: CAA10135, Genbank: CAA10136]. , the zebrafish melanocortin 1, 2, 3, 4, 5a and 5b receptors [Genbank: NP_851301.1, Genbank: NP_851302.1, Genbank: NP_851303.1, Genbank: NP_775385.1, Genbank: NP_775386.1, Genbank: NP_775387.1], and 199 mouse non-OR Class A GPCRs extracted from the GPCRDB http://www.gpcr.org/. Fugu and tetraodon OR sequences were predicted from the current genome assemblies [27, 34] using the methods described above for zebrafish.
For the calculation of percent identities, mouse and zebrafish amino acid sequences were multiply aligned and trimmed of their N- and C-terminal tails as described above. Calculations of average, minimum and maximum intra-family, inter-family, intra-subfamily and inter-subfamily percent identities were based on percent identities calculated for all pairs of amino acid sequences in this multiple alignment.
The dN/dS ratios for multi-codon regions (i.e. individual transmembrane domains or loop regions) of the odorant receptor coding sequence were determined using previously published methods . To make inferences about selective pressure (positive and negative selection) on individual codons (sites) within the coding sequence of the zebrafish OR genes, the Single Likelihood Ancestor Counting (SLAC) package http://www.datamonkey.org, which implements the Suzuki-Gojobori method , was used. Details regarding both of these methods are provided in the Supplement [see Additional file 1].
Genbank accession numbers
All sequences described in this study have been deposited in Genbank under accession numbers [Genbank: DQ305986] – [Genbank: DQ306145] (zebrafish), [Genbank:DQ306146] – [Genbank: DQ306202] (tetraodon), and [Genbank:DQ306203] – [Genbank: DQ306259] (fugu).
Buck L, Axel R: A novel multigene family may encode odorant receptors: a molecular basis for odor recognition. Cell. 1991, 65: 175-187. 10.1016/0092-8674(91)90418-X.
Mombaerts P: Genes and ligands for odorant, vomeronasal and taste receptors. Nat Rev Neurosci. 2004, 5: 263-278. 10.1038/nrn1365.
Dulac C, Axel R: A novel family of genes encoding putative pheromone receptors in mammals. Cell. 1995, 83: 195-206. 10.1016/0092-8674(95)90161-2.
Matsunami H, Buck LB: A multigene family encoding a diverse array of putative pheromone receptors in mammals. Cell. 1997, 90: 775-784. 10.1016/S0092-8674(00)80537-1.
Ryba NJ, Tirindelli R: A new multigene family of putative pheromone receptors. Neuron. 1997, 19: 371-379. 10.1016/S0896-6273(00)80946-0.
Herrada G, Dulac C: A novel family of putative pheromone receptors in mammals with a topographically organized and sexually dimorphic distribution. Cell. 1997, 90: 763-773. 10.1016/S0092-8674(00)80536-X.
Yang H, Shi P, Zhang YP, Zhang J: Composition and evolution of the V2r vomeronasal receptor gene repertoire in mice and rats. Genomics. 2005, 86: 306-315. 10.1016/j.ygeno.2005.05.012.
Zhang X, Firestein S: The olfactory receptor gene superfamily of the mouse. Nat Neurosci. 2002, 5: 124-133.
Zhang X, Rodriguez I, Mombaerts P, Firestein S: Odorant and vomeronasal receptor genes in two mouse genome assemblies. Genomics. 2004, 83: 802-811. 10.1016/j.ygeno.2003.10.009.
Glusman G, Yanai I, Rubin I, Lancet D: The complete human olfactory subgenome. Genome Res. 2001, 11: 685-702. 10.1101/gr.171001.
Malnic B, Godfrey PA, Buck LB: The human olfactory receptor gene family. Proc Natl Acad Sci U S A. 2004, 101: 2584-2589. 10.1073/pnas.0307882100.
Zozulya S, Echeverri F, Nguyen T: The human olfactory receptor repertoire. Genome Biol. 2001, 2: RESEARCH0018-10.1186/gb-2001-2-6-research0018.
Ngai J, Dowling MM, Buck L, Axel R, Chess A: The family of genes encoding odorant receptors in the channel catfish. Cell. 1993, 72: 657-666. 10.1016/0092-8674(93)90395-7.
Barth AL, Dugas JC, Ngai J: Noncoordinate expression of odorant receptor genes tightly linked in the zebrafish genome. Neuron. 1997, 19: 359-369. 10.1016/S0896-6273(00)80945-9.
Pfister P, Rodriguez I: Olfactory expression of a single and highly variable V1r pheromone receptor-like gene in fish species. Proc Natl Acad Sci U S A. 2005, 102: 5489-5494. 10.1073/pnas.0402581102.
Luu P, Acher F, Bertrand HO, Fan J, Ngai J: Molecular determinants of ligand selectivity in a vertebrate odorant receptor. J Neurosci. 2004, 24: 10128-10137. 10.1523/JNEUROSCI.3117-04.2004.
Speca DJ, Lin DM, Sorensen PW, Isacoff EY, Ngai J, Dittman AH: Functional identification of a goldfish odorant receptor. Neuron. 1999, 23: 487-498. 10.1016/S0896-6273(00)80802-8.
Dugas JC, Ngai J: Analysis and characterization of an odorant receptor gene cluster in the zebrafish genome. Genomics. 2001, 71: 53-65. 10.1006/geno.2000.6415.
Weth F, Nadler W, Korsching S: Nested expression domains for odorant receptors in zebrafish olfactory epithelium. Proc Natl Acad Sci USA. 1996, 93: 13,321-13,326. 10.1073/pnas.93.23.13321.
Vogt RG, Lindsay SM, Byrd CA, Sun M: Spatial patterns of olfactory neurons expressing specific odor receptor genes in 48-hour-old embryos of zebrafish Danio rerio. J Exp Biol. 1997, 200 ( Pt 3): 433-443.
Barth AL, Justice NJ, Ngai J: Asynchronous onset of odorant receptor expression in the developing zebrafish olfactory system. Neuron. 1996, 16: 23-34. 10.1016/S0896-6273(00)80020-3.
Byrd CA, Jones JT, Quattro JM, Rogers ME, Brunjes PC, Vogt RG: Ontogeny of odorant receptor gene expression in zebrafish, Danio rerio. J Neurobiol. 1996, 29: 445-458. 10.1002/(SICI)1097-4695(199604)29:4<445::AID-NEU3>3.0.CO;2-8.
Irie-Kushiyama S, Asano-Miyoshi M, Suda T, Abe K, Emori Y: Identification of 24 genes and two pseudogenes coding for olfactory receptors in Japanese loach, classified into four subfamilies: a putative evolutionary process for fish olfactory receptor genes by comprehensive phylogenetic analysis. Gene. 2004, 325: 123-135. 10.1016/j.gene.2003.10.011.
Freitag J, Krieger J, Strotmann J, Breer H: Two classes of olfactory receptors in Xenopus laevis. Neuron. 1995, 15: 1383-1392. 10.1016/0896-6273(95)90016-0.
Freitag J, Beck A, Ludwig G, von Buchholtz L, Breer H: On the origin of the olfactory receptor family: receptor genes of the jawless fish (Lampetra fluviatilis). Gene. 1999, 226: 165-174. 10.1016/S0378-1119(98)00575-7.
Freitag J, Ludwig G, Andreini I, Rossler P, Breer H: Olfactory receptors in aquatic and terrestrial vertebrates. J Comp Physiol [A]. 1998, 183: 635-650. 10.1007/s003590050287.
Aparicio S, Chapman J, Stupka E, Putnam N, Chia JM, Dehal P, Christoffels A, Rash S, Hoon S, Smit A, Gelpke MD, Roach J, Oh T, Ho IY, Wong M, Detter C, Verhoef F, Predki P, Tay A, Lucas S, Richardson P, Smith SF, Clark MS, Edwards YJ, Doggett N, Zharkikh A, Tavtigian SV, Pruss D, Barnstead M, Evans C, Baden H, Powell J, Glusman G, Rowen L, Hood L, Tan YH, Elgar G, Hawkins T, Venkatesh B, Rokhsar D, Brenner S: Whole-genome shotgun assembly and analysis of the genome of Fugu rubripes. Science. 2002, 297: 1301-1310. 10.1126/science.1072104.
Glusman G, Sosinsky A, Ben-Asher E, Avidan N, Sonkin D, Bahar A, Rosenthal A, Clifton S, Roe B, Ferraz C, Demaille J, Lancet D: Sequence, structure, and evolution of a complete human olfactory receptor gene cluster. Genomics. 2000, 63: 227-245. 10.1006/geno.1999.6030.
Niimura Y, Nei M: Evolutionary dynamics of olfactory receptor genes in fishes and tetrapods. Proc Natl Acad Sci U S A. 2005, 102: 6039-6044. 10.1073/pnas.0501922102.
Lancet D, Ben-Arie N: Olfactory Receptors. Current Biology. 1993, 3: 668-674. 10.1016/0960-9822(93)90064-U.
Ben-Arie N, Lancet D, Taylor C, Khen M, Walker N, Ledbetter DH, Carrozzo R, Patel K, Sheer D, Lehrach H, North MA: Olfactory receptor gene cluster on human chromosome 17: possible duplication of an ancestral receptor repertoire. Hum Mol Gen. 1994, 3: 229-235.
Young JM, Friedman C, Williams EM, Ross JA, Tonnes-Priddy L, Trask BJ: Different evolutionary processes shaped the mouse and human olfactory receptor gene families. Hum Mol Genet. 2002, 11: 535-546. 10.1093/hmg/11.5.535.
Kratz E, Dugas JC, Ngai J: Odorant receptor gene regulation: implications from genomic organization. Trends Genet. 2002, 18: 29-34. 10.1016/S0168-9525(01)02579-3.
Jaillon O, Aury JM, Brunet F, Petit JL, Stange-Thomann N, Mauceli E, Bouneau L, Fischer C, Ozouf-Costaz C, Bernot A, Nicaud S, Jaffe D, Fisher S, Lutfalla G, Dossat C, Segurens B, Dasilva C, Salanoubat M, Levy M, Boudet N, Castellano S, Anthouard V, Jubin C, Castelli V, Katinka M, Vacherie B, Biemont C, Skalli Z, Cattolico L, Poulain J, De Berardinis V, Cruaud C, Duprat S, Brottier P, Coutanceau JP, Gouzy J, Parra G, Lardier G, Chapple C, McKernan KJ, McEwan P, Bosak S, Kellis M, Volff JN, Guigo R, Zody MC, Mesirov J, Lindblad-Toh K, Birren B, Nusbaum C, Kahn D, Robinson-Rechavi M, Laudet V, Schachter V, Quetier F, Saurin W, Scarpelli C, Wincker P, Lander ES, Weissenbach J, Roest Crollius H: Genome duplication in the teleost fish Tetraodon nigroviridis reveals the early vertebrate proto-karyotype. Nature. 2004, 431: 946-957. 10.1038/nature03025.
Pilpel Y, Lancet D: The variable and conserved interfaces of modeled olfactory receptor proteins. Protein Sci. 1999, 8: 969-977.
Gilad Y, Bustamante CD, Lancet D, Paabo S: Natural selection on the olfactory receptor gene family in humans and chimpanzees. Am J Hum Genet. 2003, 73: 489-501. 10.1086/378132.
Gimelbrant AA, Skaletsky H, Chess A: Selective pressures on the olfactory receptor repertoire since the human-chimpanzee divergence. Proc Natl Acad Sci U S A. 2004, 101: 9019-9022. 10.1073/pnas.0401566101.
Nei M, Gojobori T: Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Mol Biol Evol. 1986, 3: 418-426.
Nei M: Selectionism and Neutralism in Molecular Evolution. Mol Biol Evol. 2005, Aug 24; [Epub ahead of print]:
Pond SL, Frost SD: A simple hierarchical approach to modeling distributions of substitution rates. Mol Biol Evol. 2005, 22: 223-234.
Suzuki Y, Gojobori T: A method for detecting positive selection at single amino acid sites. Mol Biol Evol. 1999, 16: 1315-1328.
Berghard A, Dryer L: A novel family of ancient vertebrate odorant receptors. J Neurobiol. 1998, 37: 383-392. 10.1002/(SICI)1097-4695(19981115)37:3<383::AID-NEU4>3.0.CO;2-D.
Hedges SB: The origin and evolution of model organisms. Nat Rev Genet. 2002, 3: 838-849. 10.1038/nrg929.
Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25: 4876-4882. 10.1093/nar/25.24.4876.
Perriere G, Gouy M: WWW-query: an on-line retrieval system for biological sequence banks. Biochimie. 1996, 78: 364-369. 10.1016/0300-9084(96)84768-7.
Guindon S, Lethiec F, Duroux P, Gascuel O: PHYML Online--a web server for fast maximum likelihood-based phylogenetic inference. Nucleic Acids Res. 2005, 33: W557-W559. 10.1093/nar/gki352.
Zmasek CM, Eddy SR: ATV: display and manipulation of annotated phylogenetic trees. Bioinformatics. 2001, 17: 383-384. 10.1093/bioinformatics/17.4.383.
Konvicka K, Campagne F, Weinstein H: Interactive construction of residue-based diagrams of proteins: the RbDe web service. Protein Eng. 2000, 13: 395-396. 10.1093/protein/13.6.395.
This work was supported by a grant from the National Institute on Deafness and Other Communications Disorders, National Institutes of Health (J.N.) and a genomics training grant from the National Institutes of Health (T.S.A.). We thank member of our lab for helpful discussions and K. Scott for comments on the manuscript.
TSA carried out the analysis. Both authors participated in the design of the study and writing of the manuscript.
Electronic supplementary material
About this article
Cite this article
Alioto, T.S., Ngai, J. The odorant receptor repertoire of teleost fish. BMC Genomics 6, 173 (2005). https://doi.org/10.1186/1471-2164-6-173
- Odorant Receptor
- Zebrafish Genome
- Vertebrate Lineage
- Codon Site