Evolution and origin of vomeronasal-type odorant receptor gene repertoire in fishes
© Hashiguchi and Nishida; licensee BioMed Central Ltd. 2006
Received: 11 July 2006
Accepted: 03 October 2006
Published: 03 October 2006
In teleost fishes that lack a vomeronasal organ, both main odorant receptors (ORs) and vomeronasal receptors family 2 (V2Rs) are expressed in the olfactory epithelium, and used for perception of water-soluble chemicals. In zebrafish, it is known that both ORs and V2Rs formed multigene families of about a hundred copies. Whereas the contribution of V2Rs in zebrafish to olfaction has been found to be substantially large, the composition and structure of the V2R gene family in other fishes are poorly known, compared with the OR gene family.
To understand the evolutionary dynamics of V2R genes in fishes, V2R sequences in zebrafish, medaka, fugu, and spotted green pufferfish were identified from their draft genome sequences. There were remarkable differences in the number of intact V2R genes in different species. Most V2R genes in these fishes were tightly clustered in one or two specific chromosomal regions. Phylogenetic analysis revealed that the fish V2R family could be subdivided into 16 subfamilies that had diverged before the separation of the four fishes. Genes in two subfamilies in zebrafish and another subfamily in medaka increased in their number independently, suggesting species-specific evolution in olfaction. Interestingly, the arrangements of V2R genes in the gene clusters were highly conserved among species in the subfamily level. A genomic region of tetrapods corresponding to the region in fishes that contains the V2R cluster was found to have no V2R gene in any species.
Our results have indicated that the evolutionary dynamics of fish V2Rs are characterized by rapid gene turnover and lineage-specific phylogenetic clustering. In addition, the present phylogenetic and comparative genome analyses have shown that the fish V2Rs have expanded after the divergence between teleost and tetrapod lineages. The present identification of the entire V2R repertoire in fishes would provide useful foundation to the future functional and evolutionary studies of fish V2R gene family.
Olfaction is a sense for recognizing environmental chemicals. In many animals, olfaction plays crucial roles in various activities, such as foraging, migration, and mating. In vertebrates, odor chemicals are perceived by three evolutionary distinct groups of seven-transmembrane G protein-coupled receptors (GPCRs). Genes encoding the main odorant receptors (ORs) form the largest multigene family in vertebrates. For example, over 1,000 distinct OR copies have been identified in the mouse genome [1, 2]. In addition to ORs, vertebrates have two distinct families of GPCRs for chemical receptors, called vomeronasal receptors family 1 (V1Rs) and family 2 (V2Rs). In mammals, these receptors are mainly expressed in the vomeronasal organ [3–6], and are considered to be used for detecting pheromones. The V1R gene repertoire has been described in several mammalian species, the numbers of intact genes varying from a few to over 150 [7, 8] among species. The repertoire of V2R genes has been described in mice and rats . The number of intact V2R genes is counted 61 in mice and 57 in rats .
In teleost fishes that lack a vomeronasal organ, on the other hand, both main odorant and vomeronasal receptors are expressed in the olfactory epithelium [10–12]. Recent database studies have revealed that fish ORs and V2Rs form multigene families with a hundred copies, respectively [13–15]. In addition, one V1R homolog has also been found in several fishes, and its expression confirmed in the olfactory epithelium .
In the previous study, we identified 88 V2R genes and pseudogenes in the zebrafish genome . This number is not small compared with the number of OR sequences in this species (133 copies) . Thus, in zebrafish, the contribution of V2R to olfaction seems substantially large. The repertoire of V2R genes in each fish is considered to reflect the ability of olfaction in fish species, given that different V2Rs bind to different sets of odor chemicals. The repertoire of V2R genes in fishes, however, is almost unknown except in zebrafish  and fugu [16, 17].
It has been indicated that fish V2Rs recognize mainly amino acids [18–20]. Consequently, they are considered to be receptors for naturally occurring odors, not pheromones, because, in fishes, amino acids are common odorant substances found in natural waters . However, recent studies on mammals have indicated that some V2Rs recognize peptides released by individuals and are used for chemical communication. For example, a peptide pheromone secreted from the extraorbital lacrimal gland of male mice (the peptide pheromone was contained in the tears of male mice) was suggested to be recognized by V2Rs . V2Rs may also be used as receptors for small peptides that serve as ligands for major histocompatibility complex (MHC) molecules . MHC-based sexual selection is also known to involve olfactory mechanisms in fishes [24, 25]. For instance, female three-spined stickleback Gasterosteus aculeatus has been suggested to assess the degree of MHC diversity of their potential partners by sensing peptides for MHC ligands . V2Rs in fish are possibly also involved in chemical communication by small peptides like MHC ligands as in mammals. Therefore, understanding the evolutionary dynamics of fish V2Rs may provide some insights into the mechanisms of sexual selection and speciation in fishes.
Repertoires of V2R genes in zebrafish, medaka, fugu, and pufferfish
As a result of the database search and gene prediction, we identified 57, 39, 38, and 23 V2R sequences in zebrafish, medaka, fugu, and pufferfish genomes, respectively. In these V2R sequences, 46 zebrafish, 22 medaka, 30 fugu, and 12 pufferfish genes were considered to be putatively functional. The proportion of putatively functional members in the V2R sequences in zebrafish (82%) and fugu (78%) were higher than in medaka (56%) and pufferfish (52%). However, these percentages may not be reliable because some of the V2R partial sequences may be attributed to the redundancy of the draft genome sequences. Indeed, several putative V2R partial sequences identified in the previous version of the zebrafish genome assembly (Zv. 4)  turned out to be redundant in the latest version of the genome assembly (Zv. 5) examined in this study.
Phylogenetic relationships of fish V2R genes
Figure 3 also indicated that all fish V2Rs formed a monophyletic clade with high bootstrap probability (96%). We defined "subfamily" within the fish V2R family on the basis of phylogenetic grouping into clades that diverged before the divergence between the lineages of zebrafish and other fishes. Within the fish V2R family, 11 subfamilies could simply be recognized (2, 3, 4, 7, 8, 9, 10, 12, 13, 14, and 16; open circles in Figure 3). In addition, five monophyletic clades specific to zebrafish, or two or three of other fishes were found (1, 5, 6, 11, and 15). These clades also seem to have originated before the divergence of the lineages of zebrafish and other fishes, because the divergences of these clades predated those of adjacent subfamilies (Figure 3). Consequently, 16 virtually equivalent subfamilies could be identified within the fish V2R clade (Figure 3) and any intact fish V2Rs were included in one of these subfamilies. Almost all of these subfamilies could also be identified in the MP and ML trees [see Additional files 1 and 2].
In addition, to determine the gene subfamily to which each of "V2R pseudogenes" and "partial sequences" belongs, a BLASTP search was conducted against putatively functional fish V2R amino acid sequences by using the translated pseudogene sequences of fish V2Rs as a query. Query sequences were assigned to the subfamily to which the best hit of the query belonged with more than 80% amino acid identity. As a result, all V2R pseudogenes identified in this study could be assigned to one of the 16 subfamilies defined by phylogenetic analysis.
Numbers of V2R genes belonging to different subfamilies in the four fishes.
Genomic distributions of V2R genes in zebrafish, medaka, fugu, and pufferfish
Comparison of the gene clusters among four fishes revealed similarities in the arrangements of fish V2R genes in the gene subfamily level (Figure 4). First, the arrangement of V2R genes in the middle of the cluster begins with subfamilies 2, 14, 4, and end with subfamily 16, in all four fishes. Second, the position and orientation of one subfamily 12 gene also seems common to all species. Such concordance is likely to reflect the gene order of the V2R cluster in the common ancestor of the four fishes.
Comparison of the genomic region containing V2R gene cluster in fish with corresponding regions in tetrapods
Large variation of V2R gene repertoire among fish species
Numbers of odorant/pheromone receptors, taste receptors, and trace amine receptors identified in the four fishes.
GPCR gene family
Family size variation (ratio)
Trace amine-associated receptor
Recently, the trace amine-associated receptor (TAAR) family has been shown to be used as odorant receptors in mice . Homologs of TAAR genes are also found in fishes , and some of them are suggested to be expressed in the olfactory epithelium in zebrafish . Therefore, also in fishes, TAAR family is likely to be used for perception of odor chemicals. Interestingly, the size of the functional TAAR repertoire varies by 7.1-fold between zebrafish and fugu (Table 2). It is possible that the differences of V2R and TAAR gene repertoires mainly cause the difference of odor sensitivity among fish species.
Evolution of the fish V2R gene family
Phylogenetic analysis of V2Rs and several "V2R-related" GPCRs (family 3 GPCRs)  revealed that the four GPCR families, GPRC6As, CaSR, family C V2R, and V2Rs (family A+B), were separated before the divergence of teleost fish and tetrapods (Figure 3). In fishes, GPRC6As and CaSRs appear to be functionally distinct from V2Rs (family A+B) because their expression profiles seem different from the fish family A+B V2R genes. It is known that in fish, CaSR is expressed mainly in the kidneys and functions as a salinity sensor [37, 38]. In goldfish, 5.24 receptor (GPRC6A ortholog) is expressed in olfactory cells, responding to naturally occurring amino acids [28, 29]. However, goldfish 5.24 receptor is expressed not only in the olfactory epithelium, but also in the gills, tongue, and other tissues . In human, GPRC6A is reported to be widely expressed in the brain and peripheral tissues with highest levels in the kidney, skeletal muscle, testis, and leucocytes . Additionally, GPRC6A and CaSR were single copy genes in the four fishes. This also suggest that these receptors are not involved in the distinction of the diverse array of odor chemicals. Fish family C V2R is also a single copy gene in all four fishes (Table 1), and thus, may be functionally distinct from other fish V2Rs (family A+B); however, the expression of the fish family C V2R is not yet known and further studies are needed.
The present phylogenetic analysis also indicated that fish V2Rs were further subdivided into at least 16 subfamilies that diverged before the divergence between the lineages of zebrafish and three other fishes (Figure 3). The number of V2R genes in each subfamily was remarkably different among species as discussed above. Generally, the evolution of vertebrate odorant and vomeronasal receptors is characterized by rapid gene turnover and lineage-specific phylogenetic clustering [8, 9, 13, 39]. These characteristics are also observed in fish V2Rs. Indeed, the phylogenetic tree shows that fish V2R genes tend to form species-specific gene clades (e.g. zebrafish subfamilies 9 and 16, medaka subfamily 4; Figures 3 and 4). These V2R subfamilies unique to each species may be related to the adaptation to species-specific environmental odor and/or pheromonal signals. On the other hand, almost one-to-one orthologous relationships were found among the four fishes in four subfamilies, 2, 7, 10, and 12 (Figure 3). These V2Rs may be used to perceive common odor substances to these fishes, like amino and nucleic acids .
In fishes, most V2R genes were located on one specific genomic region as a gene cluster (Figure 4). Interestingly, the arrangements of V2R genes in the cluster appear to be well conserved among species at the gene subfamily level (Figure 4). This suggests that the ancestral fish V2R gene cluster was already present in the most recent common ancestor of the four fishes. In zebrafish, two gene clusters were found on chromosomes 17 and 18. Two indirect pieces of evidence strongly suggest that the two zebrafish V2R gene clusters originated from one chromosomal translocation occurred in the zebrafish lineage. First, genes in each of these clusters were not phylogenetically clustered (Figure 3). Second, a genomic region near one gene in subfamily 12 (12-1) that contains several genes in subfamilies 7, 9, and 10 were lacking only in the zebrafish chromosome 18 gene cluster (Figure 4), and these genes were located on the chromosome 17 gene cluster (Figure 4).
From the synteny of neprilysin, PLC-η2, and several genes between fishes and tetrapods, we could identify the genomic region of tetrapods that corresponds to the fish genomic region containing the V2R gene cluster (Figure 5). Physical maps of this region (Figure 5) show that the V2R gene cluster located between neprilysin and PLC-η2 is fish-specific and are not found in tetrapods. This indicates that the V2R gene cluster in fishes has been originated after the separation of teleost and tetrapod lineages. The maps also show that one or more family C V2R genes are located near SLC33A1 and GMPS genes commonly found in the zebrafish, frog, mouse, and human in common. This suggests that the family A+B of fish V2Rs has been originated by a tandem duplication of the ancestral V2R gene occurred before the divergence of teleosts and tetrapods (Figure 3), and the family C of V2Rs is the sister of the family A+B (see also Figure 3). This may also indicate that one of the ancestral V2R family A+B located in the genomic region between neprilysin and PLC-η2 was increased in their number by tandem gene duplications, only in the lineage of teleosts, but not of tetrapods.
In the last part of this section, we summarize the long-term evolutionary scenario of the fish V2R gene family. Before the divergence of teleost fish and tetrapod lineages, an ancestral V2R (family A+B) gene and family C V2R gene had diverged from one common ancestral V2R gene by tandem gene duplication. After teleost-tetrapod divergence (ca. 480 MYA; Figure 1) , in the common ancestor of the four fishes, prototypes of fish V2R gene subfamilies observed in the present fishes (Figure 3) had been generated by tandem gene duplications from the ancestral V2R (family A+B) gene. Then, these prototype V2R genes formed a gene cluster. During the evolution of the lineages of the four fishes, gene duplications and losses occurred in each of these lineages. In the zebrafish lineage that had diverged from other fish lineages at ca. 320 MYA, the V2R gene cluster separated into two, and genes in subfamilies 9 and 16 have increased their numbers within the gene clusters in chromosomes 17 and 18, respectively. In the medaka lineage (diverged from fugu and pufferfish lineage at ca. 190 MYA), subfamily 4 genes have increased their copies within the gene cluster. In the fugu and pufferfish lineage, it appears that V2R sequences in some subfamilies are pseudogenized only in the pufferfish (Table 1). For instance, in subfamily 4, fugu has five putatively functional genes, but in pufferfish, four out of five are pseudogenes. Similarly, in subfamily 16, fugu has seven genes and one pseudogene, whereas pufferfish has four pseudogenes and no functional gene. This implies that, in the pufferfish lineage, after the separation from the fugu lineage (ca. 85 MYA), some V2Rs might lose their function by some biological reasons, such as the adaptation to freshwater environment. In the future, to elucidate the more detailed evolutionary process of V2R (family A+B) gene family, data from nonteleost fishes, namely "ancient fish" , seem to be useful. In addition, to clarify the divergence period of V2R family to the other family 3 GPCR families, examination of Elasmobranchii fishes (i.e. sharks and rays) may be one of important steps of studies.
Functional significance of the genes located near the fish V2R gene cluster
Two protein-coding genes, neprilysin and PLC-η2, are located near the fish V2R gene cluster (Figure 5). Interestingly, it appears that the products of these genes can interact with V2Rs. Langenau and colleagues  reported that the transcription of neprilysin was upregulated in the ovary of yellow perch (Perca flavescens), incubating it with 17α, 20β-dihydroxy-4-pregnen-3-one (17, 20-P) that stimulated fish ovulation. Neprilysin is a membrane-bound neutral peptidase expressed in various tissues . Some peptides that are cleaved by neprilysin may be released into the water from individuals and spawned eggs, and perceived by some V2Rs. Such peptides might be used as a signal in reproduction and/or identity in fishes. Recent studies in mammals have indicated that V2Rs are receptors for peptide sex pheromones  and MHC ligands . These studies may be consistent with the hypothesis that fish V2Rs are also used to recognize peptides that serve as reproductive and/or individuality signals. It is also interesting that 17, 20-P is one of the components of the pre-ovulately pheromone which affects male goldfish .
PLC-η2 is one of the η-type phospholipase C, localized specifically in neurons. V2Rs are included in the family 3 GPCRs, which is known to be coupled with phospholipase C . Western blot analysis showed that, in mouse, PLC-η2 protein was not localized in the olfactory epithelium or the vomeronasal organ . This suggests that PLC-η2 is not coupled with V2Rs; however, it was also shown that PLC-η2 was localized in the main- and sub-olfactory bulbs . This implies some functional relationship between PLC-η2 and the odorant recognition system in vertebrates.
The close linkage between the genes encoding V2Rs and neprilysin/PLC-η2 in fishes seems to be an intriguing finding in light of the function of fish V2R genes. Functional studies of neprilysin and PLC-η2 could provide further information on the evolutionary mechanisms of fish V2Rs.
In this paper, we have presented the repertoire of V2R gene family in the four teleost fish species. The number of intact V2R genes in these fishes varied from 12 in pufferfish to 46 in zebrafish. Phylogenetic analysis has shown that the evolution of fish V2Rs are characterized by rapid gene turnover and lineage-specific clustering, as known in ORs and V1Rs [8, 9, 13, 39]. Such evolutionary patterns may suggest that fish V2Rs are involved in the lineage-specific adaptation to different odor environments for each fish species. In addition, we elucidated substantial parts of the long-term evolution of fish V2R gene family, by phylogenetic and comparative genome analyses. We believe that our results provide one of the foundation to the future functional and evolutionary studies of fish V2R gene family.
Identification of fish V2R genes from draft genome sequences
The 6.5× coverage of the zebrafish genome sequence (Zv.5), 8.7× coverage of the fugu genome sequence (Fugu version 4), and 8.3× coverage of the pufferfish genome sequence (Genoscope Tetraodon 7) are available at ENSEMBL . The 6.7× medaka genome sequence is available from the National Institute of Genetics (NIG) DNA Sequencing Center, Mishima, Japan . First, a TBLASTN search was conducted with E-value 10-10 against the genomic data by using the TM domain of several representative V2R amino acid sequences known in fugu, goldfish, and zebrafish as queries. Obtained sequences were confirmed by BLASTP searches against the NCBI non-redundant (nr) database. In each of these sequences, if the best hit in this search was a previously known fish V2R, it was considered a V2R sequence. Second, each region of BLAST similarity was extended in 5' and 3' directions to perform a detailed prediction of V2R coding sequences. Each region of BLAST similarity was extended 4 kb in the 5' and 1 kb to the 3' direction. Gene prediction, described as follows, was conducted in each of these extended regions. After the first round of gene prediction, if the total coding sequence of the V2R gene was not included in this region, the 5' sequence was further extended up to 10 kb.
For each of these regions, V2R coding sequences were estimated based on the profile HMM-based gene prediction with the program WISE2 . Seven full-length fugu  and goldfish  V2Rs were aligned using the program ClastalW  with the default settings. A profile HMM was constructed from the alignment by using the program HMMER software package , and used for gene prediction. In addition to V2R (families A+B and C) homologs, CaSRs and a homolog of goldfish 5.24 amino acid receptors (GPRC6A)  were also identified. These "V2R-related" GPCR family 3 receptors were included in phylogenetic analysis. The obtained putative protein sequences were examined by the TMHMM method  for the presence of seven transmembrane domains. Chromosomal positions of putative V2R genes and pseudogenes in zebrafish and pufferfish could be determined by mapping them onto chromosome contigs. The list and sequences of fish V2Rs are available as supplementary information [Additional files 3, 4, 5, 6, 7, 8, 9, 10].
Deduced amino acid sequences of 90 putatively functional V2R genes (including family C V2Rs) in zebrafish, medaka, fugu, and pufferfish and several related GPCRs (CaSRs, GPRC6As, and T1R taste receptors in some vertebrates) were aligned by using the program FFT-NS-i (MAFFT 5.731)  and slightly modified by eyes. Seventeen vertebrate T1Rs were used as outgroups. Phylogenetic tree was constructed using the neighbor-joining method  with JTT matrix distances  implemented in the program MEGA 3.1 . The reliability of each tree node was assessed by the bootstrap method with 1,000 replications. The MP and ML trees were reconstructed by using reduced dataset consisted of 45 representative fish V2R (including family C) genes, three tetrapod V2R genes, and nine related GPCR family 3 genes, because in these methods, it is virtually impossible to reconstruct trees from the full dataset (134 OTUs). Initially, an alignment of nucleotide sequences of these genes were constructed based on the alignment of corresponding amino acid sequences. For the ML methods, the optimal model of sequence evolution was determined to be the Tamura-Nei + I + gamma model by the hierarchical likelihood-ratio test implemented in MODELTEST 3.7 . The ML tree was searched by a heuristic-search algorism implemented in PAUP 4.10 b . Similar to the ML tree reconstruction, the MP tree was also obtained by a heuristic-search algorism implemented in PAUP 4.10 b. The reliability of MP tree nodes was assessed by the bootstrap method with 1,000 replications. In the MP and ML trees, the fugu T1R1 sequence was used as an outgroup.
Identification of tetrapod genomic regions corresponding to the region in fishes with the V2R gene cluster
Initially, several protein-coding genes were characterized around V2R gene clusters in zebrafish, medaka, and pufferfish, by using ENSEMBL genome browser . Fugu was not included in this analysis because fugu V2R clusters were separated into at least four unconnected scaffolds and thus we could not identify the gene order around V2R clusters. Next, putative orthologs of these genes were identified in the genomic sequences of the frog, chicken, mouse, and human by TBLASTN searches. The identities of these putative orthologs were confirmed by "reciprocal best hits" of TBLASTN searches between fishes and tetrapods.
The zebrafish, fugu, pufferfish, and medaka sequence data were produced by the Sanger Institute, International Fugu Genome Consortium, Genoscope and the Broad Institute at MIT, and the National Institute of Genetics (NIG) at Mishima, Japan, respectively. We thank Y. Yamanoue for permission to use the fish illustrations. We also thank R. Kawahara and H. Takeshima for helpful discussion, and two anonymous reviewers for variable comments on the earlier version of the manuscript. This work was partially supported by Grants-in-Aid from the Japan Society for the Promotion of Science to MN.
- 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.View ArticlePubMedGoogle Scholar
- Godfrey PA, Malnic B, Buck LB: The mouse olfactory receptor gene family. Proc Natl Acad Sci USA. 2004, 101: 2156-2161. 10.1073/pnas.0308051100.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Herrada G, Dulac C: A novel family of putative pheromone receptors in mammals with a topographically originated and sexually dimorphic distribution. Cell. 1997, 90: 763-773. 10.1016/S0092-8674(00)80536-X.View ArticlePubMedGoogle Scholar
- Ryba NJP, Tilindelli R: A new multigene family of putative pheromone receptor. Neuron. 1997, 19: 371-379. 10.1016/S0896-6273(00)80946-0.View ArticlePubMedGoogle Scholar
- Rodriguez I, Del Punta K, Rothman A, Ishii T, Mombaerts P: Multiple new and isolated families within the mouse superfamily of V1r vomeronasal receptors. Nat Neurosci. 2002, 5: 134-140. 10.1038/nn795.View ArticlePubMedGoogle Scholar
- Grus WE, Shi P, Zhang Y-p, Zhang J: Dramatic variation of the vomeronasal pheromone receptor gene repertoire among five orders of placental and marsupial mammals. Proc Natl Acad Sci USA. 2005, 102: 5767-5772. 10.1073/pnas.0501589102.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang H, Shi P, Zhang Y-p, 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.View ArticlePubMedGoogle Scholar
- Cao Y, Oh BC, Stryer L: Cloning and localization of two multigene receptor families in goldfish olfactory epithelium. Proc Natl Acad Sci USA. 1998, 95: 11987-11992. 10.1073/pnas.95.20.11987.PubMed CentralView ArticlePubMedGoogle Scholar
- Asano-Miyoshi M, Suda T, Yasuoka A, Osima S, Yamashita S, Abe K, Emori Y: Random expression of main and vomeronasal olfactory receptor genes in immature and mature olfactory epithelia of Fugu rubripes. J Biochem Tokyo. 2000, 127: 915-924.View ArticlePubMedGoogle Scholar
- Pfister P, Rodriguez I: Olfactory expression of a single and highly variable V1r pheromone receptor-like gene in fish species. Proc Natl Acad Sci USA. 2005, 102: 5489-5494. 10.1073/pnas.0402581102.PubMed CentralView ArticlePubMedGoogle Scholar
- Niimura Y, Nei M: Evolutionary dynamics of olfactory receptor genes in fishes and tetrapods. Proc Natl Acad Sci USA. 2005, 102: 6039-6044. 10.1073/pnas.0501922102.PubMed CentralView ArticlePubMedGoogle Scholar
- Alioto TS, Ngai J: The odorant receptor repertoire of teleost fish. BMC Genomics. 2005, 6: 173-10.1186/1471-2164-6-173.PubMed CentralView ArticlePubMedGoogle Scholar
- Hashiguchi Y, Nishida M: Evolution of vomeronasal-type odorant receptor genes in the zebrafish genome. Gene. 2005, 362: 19-28. 10.1016/j.gene.2005.07.044.View ArticlePubMedGoogle Scholar
- Naito T, Saito Y, Yamamoto J, Nozaki Y, Tomura K, Hanzawa M, Nakanishi S, Brenner S: Putative pheromone receptors related to the Ca2+-sensing receptor in Fugu. Proc Natl Acad Sci USA. 1998, 95: 5178-5181. 10.1073/pnas.95.9.5178.PubMed CentralView ArticlePubMedGoogle Scholar
- Bjarnadottir TK, Fredriksson R, Schioth HB: The gene repertoire and the common evolutionary history of glutamate, pheromone (V2R), taste(1) and other related G protein-coupled receptors. Gene. 2005, 362: 70-84. 10.1016/j.gene.2005.07.029.View ArticlePubMedGoogle Scholar
- Hansen A, Rolen SH, Anderson K, Morita Y, Caprio J, Finger TE: Correlation between olfactory receptor cell type and function in the channel catfish. J Neurosci. 2003, 23 (28): 9328-9339.PubMedGoogle Scholar
- Hansen A, Anderson KT, Finger TE: Differential distribution of olfactory receptor neurons in goldfish: structures and molecular correlates. J Comp Neurol. 2004, 477: 347-359. 10.1002/cne.20202.View ArticlePubMedGoogle Scholar
- Eisthen HL: The goldfish knows: Olfactory receptor cell morphology predicts receptor gene expression. J Comp Neurol. 2004, 477: 341-346. 10.1002/cne.20258.View ArticlePubMedGoogle Scholar
- Hara TJ: Olfaction and gustation in fish: an overview. Acta Physiol Scand. 1994, 152: 207-217.View ArticlePubMedGoogle Scholar
- Kimoto H, Haga S, Sato K, Touhara K: Sex-specific peptides from exocrine glands stimulate mouse vomeronasal sensory neurons. Nature. 2005, 437: 898-901. 10.1038/nature04033.View ArticlePubMedGoogle Scholar
- Leinders-Zufall T, Brennan P, Widmayer P, Chndramani SP, Maul-Pavicic A, Jager M, Li X-H, Breer H, Zufall F, Boehm T: MHC class I peptides as chemosensory signals in the vomeronasal organ. Science. 2004, 306: 1033-1037. 10.1126/science.1102818.View ArticlePubMedGoogle Scholar
- Reusch TBH, Haberll MA, Aeschlimann PB, Milinski M: Female sticklebacks count alleles in a strategy of sexual selection explaining MHC polymorphism. Nature. 2001, 414: 300-302. 10.1038/35104547.View ArticlePubMedGoogle Scholar
- Aeschlimann PB, Haberli MA, Reusch TBH, Boehm T, Milinski M: Female stickleback Gasterosteus aculeatus use self-reference to optimize MHC allele number during mate selection. Behav Ecol Sociobiol. 2003, 54: 119-126.Google Scholar
- Milinski M, Griffiths S, Wegner KM, Reusch TBH, Hass-Assembaum A, Boehm T: Mate choice decisions of stickleback females predictably modified by MHC peptide ligands. Proc Natl Acad Sci USA. 2005, 102: 4414-4418. 10.1073/pnas.0408264102.PubMed CentralView ArticlePubMedGoogle Scholar
- Yamanoue Y, Miya M, Inoue JG, Matsuura K, Nishida M: The mitochondrial genome of spotted green pufferfish Tetraodon nigroviridis (Teleostei: Tetraodontiformes) and divergence time estimation among model organisms in fishes. Genes Genet Syst. 2006, 81: 29-39. 10.1266/ggs.81.29.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Luu P, Acher F, Bertrand H-O, 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.View ArticlePubMedGoogle Scholar
- Wellendorph P, Brauner-Osborne H: Molecular cloning, expression, and sequence analysis of GPRC6A, a novel family C G-protein-coupled receptor. Gene. 2004, 335: 37-46. 10.1016/j.gene.2004.03.003.View ArticlePubMedGoogle Scholar
- Shi P, Zhang J: Contrasting modes of evolution between vertebrate sweet/umami receptor genes and bitter receptor genes. Mol Biol Evol. 2006, 23: 292-300. 10.1093/molbev/msj028.View ArticlePubMedGoogle Scholar
- Ishimaru Y, Okada S, Naito H, Nagai T, Yasuoka A, Matsumoto I, Abe K: Two families of candidate taste receptors in fishes. Mech Dev. 2005, 122: 1310-1321. 10.1016/j.mod.2005.07.005.View ArticlePubMedGoogle Scholar
- Gloriam DEI, Bjarnadottir TK, Yan Y-L, Postlethwait JH, Schioth HB, Fredriksson RF: The repertoire of trace amine G-protein-coupled receptors: large expansion in zebrafish. Mol Phylogenet Evol. 2005, 35: 470-482. 10.1016/j.ympev.2004.12.003.View ArticlePubMedGoogle Scholar
- Mundy NI: Genetic Basis of olfactory communication in primates. Am J Primatol. 2002, 11: 535-546.Google Scholar
- Liberles SD, Buck L: A second class of chemosensory receptors in the olfactory epithelium. Nature. 2006, 442: 645-650. 10.1038/nature05066.View ArticlePubMedGoogle Scholar
- Pin J-P, Galvez T, Prezeau L: Evolution, structure, and activation mechanism of family 3/C G-protein-coupled receptors. Pharmacol Ther. 2003, 98: 325-354. 10.1016/S0163-7258(03)00038-X.View ArticlePubMedGoogle Scholar
- Nearing J, Betka M, Quinn S, Hentschel H, Elger M, Baum M, Bai M, Chattopadyhay N, Brown EM, Hebert SC, Harris HW: Polyvalent cation receptor proteins (CaRs) are salinity sensors in fish. Proc Natl Acad Sci USA. 2002, 99: 9231-9236. 10.1073/pnas.152294399.PubMed CentralView ArticlePubMedGoogle Scholar
- Loretz CA, Pollina C, Hyodo S, Takei Y, Changl W, Shoback D: cDNA cloning and functional expression of a Ca2+-sensing receptor with truncated C-terminal tail from the Mozambique tilapia (Oreochromis mossambicus). J Biol Chem. 2004, 279: 53288-53297. 10.1074/jbc.M410098200.View ArticlePubMedGoogle Scholar
- Grus WE, Zhang J: Rapid turnover and species-specificity of vomeronasal pheromone receptor genes in mice and rats. Gene. 2004, 340: 303-312. 10.1016/j.gene.2004.07.037.View ArticlePubMedGoogle Scholar
- Inoue JG, Miya M, Tsukamoto K, Nishida M: Basal actinopterygian relationships: a mitogenomic perspective on the phylogeny of the "ancient fish". Mol Phylogenet Evol. 2003, 26: 110-120. 10.1016/S1055-7903(02)00331-7.View ArticlePubMedGoogle Scholar
- Langenau DM, Goetz FW, Roberts SB: The upregulation of messenger ribonucleic acids during 17α, 20β-dihydroxy-4-pregnen-3-one-induced ovulation in the perch ovary. J Mol Endocrinol. 1999, 23: 137-152. 10.1677/jme.0.0230137.View ArticlePubMedGoogle Scholar
- Sorensen PW, Stacey NE: Brief review of fish pheromones and discussion of their possible uses in the control of non-indigenous teleost fishes. NZ J Mar Freshwat Res. 2004, 38: 399-417.View ArticleGoogle Scholar
- Nakahara M, Shinozawa M, Nakamura Y, Irino Y, Morita M, Kudo Y, Fukami K: A novel phospholipase C, PLC-η2, is a neuron-specific isozyme. J Biol Chem. 2005, 280: 29128-29134. 10.1074/jbc.M503817200.View ArticlePubMedGoogle Scholar
- ENSEMBL Genome browser. [http://www.ensembl.org]
- National Institute of Genetics DNA Sequencing Center (Medaka site). [http://dolphin.lab.nig.ac.jp/medaka]
- Birney E, Clamp M, Durbin R: GeneWise and Genomewise. Genome Res. 2004, 14: 988-995. 10.1101/gr.1865504.PubMed CentralView ArticlePubMedGoogle Scholar
- Thompson JD, Higgins GD, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, positions-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680.PubMed CentralView ArticlePubMedGoogle Scholar
- HMMER Software Package. [http://hmmer.wustl.edu]
- Sonnhammer EL, von Heijne G, Krogh A: A hidden Markov model for predicting transmembrane helices in protein sequences. Proc Int Conf Intell Syst Mol Biol. 1998, 6: 175-182.PubMedGoogle Scholar
- Katoh K, Misawa K, Kuma K, Miyata T: MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002, 30: 3059-3066. 10.1093/nar/gkf436.PubMed CentralView ArticlePubMedGoogle Scholar
- Saitou N, Nei M: The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.PubMedGoogle Scholar
- Jones DT, Taylor WR, Thornton JM: The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci. 1992, 8: 275-282.PubMedGoogle Scholar
- Kumar S, Tamura K, Nei M: MEGA 3: Integrated software for molecular evolutionary genetics analysis and sequence alignment. Briefings in Bioinformatics. 2004, 5: 150-163. 10.1093/bib/5.2.150.View ArticlePubMedGoogle Scholar
- Posada D, Crandall KA: Modeltest: testing the model of DNA substitution. Bioinformatics. 1998, 14: 187-188. 10.1093/bioinformatics/14.9.817.View ArticleGoogle Scholar
- Swofford DL: PAUP*: Phylogenetic analysis using parsimony, v4.0b10. 2001, Smithsonian Institution and Sinauer Associates, SunderlandGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.