|
Copyright Pfister et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Divergent Evolution among Teleost V1r Receptor Genes Suzannah Rutherford, Academic Editor Fred Hutchinson Cancer Research Center, United States of America * To whom correspondence should be addressed. E-mail: ivan.rodriguez/at/zoo.unige.ch Conceived and designed the experiments: IR JM PP. Performed the experiments: JM PP JR. Analyzed the data: IR JM PP JR. Wrote the paper: IR JM PP. Received February 9, 2007; Accepted March 20, 2007. | |||||||||||||||
Abstract The survival of vertebrate species is dependent on the ability of individuals to adequately interact with each other, a function often mediated by the olfactory system. Diverse olfactory receptor repertoires are used by this system to recognize chemicals. Among these receptors, the V1rs, encoded by a very large gene family in most mammals, are able to detect pheromones. Teleosts, which also express V1r receptors, possess a very limited V1r repertoire. Here, taking advantage of the possibility to unequivocally identify V1r orthologs in teleosts, we analyzed the olfactory expression and evolutionary constraints of a pair of clustered fish V1r receptor genes, V1r1 and V1r2. Orthologs of the two genes were found in zebrafish, medaka, and threespine stickleback, but a single representative was observed in tetraodontidae species. Analysis of V1r1 and V1r2 sequences from 12 different euteleost species indicate different evolutionary rates between the two paralogous genes, leading to a highly conserved V1r2 gene and a V1r1 gene under more relaxed selective constraint. Moreover, positively-selected sites were detected in specific branches of the V1r1 clade. Our results suggest a conserved agonist specificity of the V1R2 receptor between euteleost species, its loss in the tetraodontidae lineage, and the acquisition of different chemosensory characteristics for the V1R1 receptor. | |||||||||||||||
In vertebrates, interindividual interactions related to reproduction are largely dependent on pheromonal communication. Chemosensory structures located in mammals in the nasal cavity and in teleosts in the olfactory rosette represent the major tools used for such exchanges. These structures contain thousands of olfactory sensory neurons, organized in pseudostratified neuroepithelia. Each olfactory sensory neuron extends a single dendrite towards the outside world and an axon which directly connects to the olfactory bulb, in the brain. A single or a few members of a remarkably large olfactory receptor gene repertoire are expressed by each olfactory sensory neuron. Thus, chemical information from the outside world is transported by multiple parallel lines to the olfactory bulb, where it is processed, directed towards various brain areas, and translated into a conscious or unconscious perception. Four different classes of olfactory receptors have been described in teleosts: odorant, trace amine-associated, vomeronasal type 2 (V2r) and vomeronasal type 1 (V1r) receptors. Odorant receptors number 143 in zebrafish [1], and represent, as in mammals, one of the largest gene families in the genome. Trace amine-associated receptors [2], are also present in fish species (57 genes in zebrafish) [3]. V2rs react to amino acids in teleosts, and are relatively numerous (at least 24 potentially functional V2rs in zebrafish) [4]–[7]. V1rs mediate pheromone detection in mice [8], [9] and are thought to play a similar role in other species. In mammals, V1r genes are highly divergent, often species-specific, and in rodents number over 100 [10]–[14]. The role played by the size and variability of these V1r repertoires is not understood yet. Thus, basic questions, such as the degree of potential redundancy between paralogous V1r genes, remain unanswered. However, positive Darwinian pressure has been shown to act on some of these receptor genes, between paralogous or orthologous sequences [12], [15]–[19], indicating, at least that for some novel genes, novel functions were acquired. The variability and profusion of V1r genes in some species result from a high rate of gene birth and death[12], [13], [19]. This unusual characteristic renders the identification of orthologous V1r pairs, if existent at the functional level, difficult to perform with certainty, even in species as closely related as mice and rats. Despite the extensive genomic duplications which affected teleost genomes [20]–[24], fish V1r genes appear not to have gone through the flourishing diversification and amplification characterizing mammalian pheromone receptor genes. This peculiar situation allows the study of orthologous V1rs between euteleost species (Table 1) whose radiation dates back over 110 mio. years (Figure 1) [25]–[28]. We here report such an analysis.
| |||||||||||||||
A Linked Pair of Teleost V1r Receptor Genes A single V1r gene, based on the zebrafish (Danio rerio, Dr) (Zv4, 2005), medaka (Oryzias latipes) (Medaka 1, 2005), tetraodon (Tetraodon nigroviridis, Tn) (Tetraodon 7, 2004) and fugu (Takifugu rubripes, Tr) (Fugu 2.0, 2004) assemblies, was previously identified in teleosts [16]. Taking advantage of the newly released zebrafish (Zv6, 2006), medaka (Medaka 1, 2006) and threespine stickleback (Gasterosteus aculeatus, Ga) (BROAD S1, 2006) assemblies (see Materials and Methods), and including the fugu (Fugu 4.0, 2005) and tetraodon databases, we performed a search for potential novel V1r sequences. Databases were mined with the TBLASTN and BLASTN algorithms using as queries all known fish V1r, mouse V1rb2, V1rf3 and V1re4 sequences. Our criteria for inclusion were an uninterrupted coding sequence of at least 850 nucleotides contained in a single exon (all known vertebrate V1r genes exhibit such characteristics), a potential translated product containing 7 transmembrane segments, and an expected value inferior to 10−2.We identified two different V1r sequences fulfilling these criteria in Dr, and found their corresponding orthologs in Ga and Ol, one of them corresponding to a previously described zebrafish V1r [16]. A single V1r gene was found in Tn and Tr. We named the V1r common to all tested species V1r1 and the one missing in the tetraodontidae lineage V1r2. We hypothesized the existence of a potential pseudogenized or distant V1r2 sequence in tetraodontidae species, and tested this hypothesis by searching for sequences with a potential to encode 10 amino acid triplets conserved in Ol, Dr, and Ga V1r2: HLA (49), LTR (59), PQT (64), GCK (80), NMA (140), APR (153), GFC (167), TRA (228), FGI (245) and PVV (263), in sequences surrounding tetraodontidae V1r1 genes (100kb upstream and downstream of tetraodon V1r1, and 11kb upstream and 10kb downstream of fugu V1r1). We estimated that for a sequence to be considered potentially related to a V1r2 gene, at least 4 of the 10 conserved triplets were to be present, irrespective of the reading frame or the transcription orientation, and with a distance between each triplet not exceeding+/−12 amino acids relative to the one in the known V1r2 sequences. No sequence fulfilling these criteria was identified, suggesting the absence of a V1r2 gene in the vicinity of the V1r1 tetraodontidae gene (but without excluding the potential presence of a highly degenerate V1r2 pseudogene). We also performed PCRs at low stringency and/or used degenerate primers (see Materials and Methods) on genomic Tr and Tn DNA, using V1r2 primers located at positions conserved in all other identified fish V1r2 genes. No amplicon corresponding to a potential V1r2 gene or pseudogene was obtained, again supporting the potential absence of this gene in Tetraodontidae. However, our own experience makes us wary of such negative observations [16]. Mammalian V1Rs exhibit 10 highly conserved amino acid residues (some of which are exclusively observed in V1Rs), and a conserved glycosylation site in extracellular loop 2 [14]. 7 of the conserved residues and the glycosylation site were present in all teleost V1r sequences (Figure 2).
To dispose of a larger and therefore more reliable dataset, we then expanded our analysis to species pertaining to a single genus, Danio. Using PCR primers specific for Dr V1r1 and V1r2 sequences, we investigated the presence of the two genes in the following species: D. albolineatus (Da), D. choprae (Dc), D. dangila (Dd), D. kerri (Dk), D. kyathit (Dky), D. malabaricus (Dm), D. nigrofasciatus (Dn), D. pantheri (Dp) and D. yoma (Dy) (Table 1).V1r1 and V1r2 sequences were found in all tested Danio species. Analysis of the genomic location of the two V1r genes indicated that they were clustered on chromosomes 22 in zebrafish, 5 in medaka, and 17 in stickleback. Their transcription directions were inverted in all cases observed, in a bidirectional orientation, with an unusual short distance between coding sequences [14], ranging from 1.1kb (Dm) to 3.4 kb (Ol) (Figure 3).
A Conserved Sequence in the 5′UTR of V1r1 V1r receptor genes in mice exhibit a characteristic not observed in odorant receptor genes: an unusual conservation of non-coding sequences between members of a given family, which includes transcribed and non-transcribed regions [18]. We evaluated the potential homology of non-coding sequences between orthologous teleost V1r1 genes (Dr, Ol, Ga, Tr, Dm, and Tn) by pairwise comparisons of the sequences using the Pipmaker analysis tool. We identified a 75 nucleotide segment located 5′ to the V1r1 coding sequence containing 31 invariant nucleotides (Figure 4). The position of this conserved sequence relative to the start site ranged from -81 (Dr) to -185 (Ol) (Figure 3). 5′ RACE was performed and one transription start site was identified 404 bp upstream of the V1r1 first ATG codon (Figure 3), indicating that the conserved sequence was part of the V1r1 5′UTR.
The conserved segment was not found duplicated 5′ toV1r2 sequence, and was not identified associated with any V1r mouse gene. No significant conservation of other non-coding segments was identified neighboring V1r2 sequences. Expression of V1r1 and V1r2 in the Zebrafish Olfactory Rosette We investigated the potential expression of V1r1 and V1r2 in the olfactory rosette using two approaches. We first performed RT-PCRs on olfactory rosette extracts from male and female adult zebrafish using primers specific for V1r1 or V1r2. V1r1 and V1r2 transcripts were consistently observed in this tissue (Figure 5a). No transcripts were found in other organs, including barbels, lips, heart, gills, muscle and brain (data not shown). Second, and in order to obtain a more precise picture of the cells expressing the V1r receptor genes, we performed in situ hybridizations on olfactory rosette sections. Single neurons located in the sensory neuroepithelium did express V1r1 transcripts, similarly to what was observed for cells expressing V1r2 (Figure 5b–d). No evidence of sexual dimorphism in the expression of the V1r receptor genes was observed (data not shown).
Unequal Evolutionary Rates Between V1r1 and V1r2 Two phylogenetic trees, including 24 V1r1 and V1r2 teleost genes or their corresponding proteins, were generated using Maximum Likelihood methods (Figure 6). We evaluated potential evolutionary differences between the V1r1 and V1r2 clades by performing a Relative Rate Test, and found that the V1r1 clade evolved significantly faster than the V1r2 clade (p![]() ![]()
To determine the origin of the acceleration detected in the V1r1 clade, we performed a series of tests of increasing sensitivity, looking at the whole dataset or at specific clades, at full-length sequences or codon by codon, and aimed at measuring the selective pressures acting on these genes. Two potential explanations for the observed differential evolutionary rates, which are not mutually exclusive, involve relaxation of the purifying selection, and positive Darwinian selection. The dN/dS ratio (ω) for the overall dataset was first calculated using the one-ratio model (M0) as implemented in the PAML software (see Materials and Methods). An ω of 0.099 was estimated, suggesting that a strong purifying selection is globally acting on V1r receptors. Using the two-ratio model, the ω estimates for the V1r1 and V1r2 clades were 0.147 and 0.050 respectively (Figure 6). The Likelihood Ratio Test (LRT) indicated that the two-ratio model fit the dataset significantly better than the one-ratio model (Table 2), confirming that the V1r1 clade has a 3×times higher rate of non-synonymous mutations. However, no evidence for positive selection (ω>1) was found.
To further investigate potential changes in selective pressure on specific residues, we applied site-specific models to our dataset (see Materials and Methods). The comparison of different site-specific models suggested heterogeneous evolutionary rates within V1r sequences. Thus, the M3 model (discrete with K We then analyzed our dataset with a site-specific approach, but looking at specific branches. We used the branch-site model D (see Materials and Methods) to identify potential selective pressure differences for a category of sites between the two paralogous clades. The LRTs involving model D (with K Thus, models with additional categories of ω, applied either to the whole dataset or to specific groups within the tree, were systematically preferred to models of lower complexities (Table 2). The potential existence of rare positively selected sites in a subset of branches was then assessed. Using again a branch-site model, model A (see Materials and Methods), we analyzed the 21 longest internal and terminal branches of the V1r1/V1r2 clades. We found 5 branches exhibiting positively selected sites in the V1r1 clade and a single branch in the V1r2 clade. Among these branches, we were able to identify 9 significant positively selected sites, 5 of which were located in transmembrane regions (Figure 2). All significant sites pertained to the 5 branches of the V1r1 clade (Table 3, Figure 6).
In accordance with our previous estimates (results of model M8), the percentage of these residues was limited, ranging from 0.4% to 7.0% (Table 3). 1 to 3 sites per branch displayed an ω significantly above 1 (BEB calculations) (Table 3). Our approach, involving the use of models of increasing finesse, thus led to the confirmation of differential evolutionary rates between paralogous V1r clades, and to the detection of positive selection between orthologous V1r genes. | |||||||||||||||
We here describe two linked teleost V1r genes, V1r1 and V1r2, and identify their orthologs in multiple fish species, taking advantage of their physical proximity and relative conservation. Apparently, not all teleosts share the same V1r repertoire; we were indeed unable to identify a V1r2 receptor gene in the two tested tetraodontidae species. This negative observation may naturally reflect a failure to identify an existing sequence, due to a still incomplete coverage of the tetraodon and fugu genomes. This seems however unlikely since these genomes are close to being fully sequenced. Moreover, since V1r1 and V1r2 genes are found in close association in all teleost species tested (less than 3.5 kb apart), and since relatively large genomic sequences surrounding the tetraodon and fugu V1r1 were available and searched, we propose that our findings reflect a gene loss, which specifically affected the tetraodontidae lineage. This potential genus-specific gene loss is consistent with the remarkably compact genome of members of this lineage, genome that can be as small as 390 Mb in fugu and 340 Mb in tetraodon (while zebrafish and medaka genome sizes are 1.7 Gb and 800 Mb respectively). This limited genomic material translates into a contraction of gene number, shrinking known to affect many gene families, including olfactory receptor gene repertoires. Thus, a comparison between fugu and zebrafish odorant receptor genes indicates a drastic reduction in the fugu repertoire size (143 vs. 44 and 42 odorant receptor genes in zebrafish, fugu and tetraodon respectively) [1]. We identified a short sequence in the 5′UTR of V1r1 transcripts, conserved between multiple teleost species, including some whose common ancestor dates back over 100 mio. years. This observation is reminiscent of the remarkable 5′UTR conservation present in mouse V1ra and V1rb transcripts [18]. It may point to a transcriptional or translational regulatory element acting on fish V1r1 genes, and possibly also on V1r2. We also report a peculiar genomic organization of the fish V1r1 and V1r2 genes. As little as 1.1 kb separates the coding sequences of the two genes, leaving limited space for two 5′UTRs and two promoters since their transcriptional direction is inverted. The zebrafish V1r intergenic sequence, whose very limited size renders it easily amenable to genetic manipulation, will surely prove to be of interest in the study of the mechanisms regulating V1r expression, which at present remain completely unknown. We observed that unequal selective pressures affected the evolution of the V1r1 and V1r2 clades. The V1r1 clade exhibits a relatively relaxed negative pressure, with on given branches, specific residues under positive Darwinian selection, a situation observed for some mammalian V1rs. This finding suggests a possible variable role played by orthologous V1R1 receptors, or at least non-homogenous agonist response profiles between the latter (due for example to species-specific coevolution of agonist-receptor pairs). The situation is clearly different for the V1r2 clade: a strong purifying pressure is at work on V1R2 receptors, suggesting a common role and/or agonist profile between orthologous V1R2 proteins. Such theories will be put to test as soon as we identify natural agonists for fish and mammalian V1rs. But the road is apparently long before we dispose of a general view of agonist-V1r receptor pairs in vertebrates. Indeed, despite an identification of the first V1r genes over a decade ago [29], only a single such pair has been identified [8]. This situation, partly resulting from the difficulties faced to express these receptors in vitro, also reflects a more problematic situation, which is the very limited number of known potential agonists. | |||||||||||||||
Identification of V1r genes The following ENSEMBL databases were searched using the BLASTN and TBLASTN algorithms: Dr Zv6 (August 2006 release, 7×coverage), Ga BROAD S1 (June 2006 release, 11×coverage), Ol MEDAKA1 (May 2006 release, 6.7×coverage), Tn Tetraodon 7 (September 2004 release, 8.3×coverage), and Tr Fugu 4.0 (May 2005 release, 90% coverage). Database mined sequences from Dr, Ol and Ga were PCR amplified from genomic DNA and sequenced for confirmation. Da, Dc, Dd, Dk, Dky, Dm, Dn, Dp and Dy V1r1 and V1r2 were PCR amplified from genomic DNA using the following primers: TTC CCG GTA ACA CCG CTG TCA TCT G, GAA GAA GAG CCG GAG GTC AGT GAT CAG, TAT GGA CCT GTG TGT CAC, CAG CCT TTA TGA AAC ACA TTC AC, TGG ACC TGT GTG TCA CCA TCA AAA GG, TCA GTT CTT GCC GCT GGA GTT CTT GCC, and sequenced. The Dm intergenic sequence was amplified from genomic DNA and sequenced. Conditions for degenerate PCRs performed on fugu and tetraodon genomic DNA were: 4 min. at 95°C followed by 48 cycles of 45 sec. at 95°C, 4 min. at 50°C and 3 min at 72°C, followed by a final extension time of 10 min at 72°C, using the following primers: TCT GSM TCA CCT KCA TKC TGA GYG YST WCC A and GGG CTG AGV GMK GCR TAS MAY GAG GAR AAA AA. Low stringency amplifications were performed with an annealing temperature of 52°C. Sequences were deposited in GenBank with the following accession numbers: DQ887609, DQ887610, DQ887610, DQ887612, DQ887613, AY279523, Q887614, DQ887615, DQ887616, 880989, 884592, 884596, 884602, 884604, 884606, 881043, 881043, 881049, 881053, 884622, and 884628.RT-PCR mRNA was extracted (RNeasy extraction kit, Qiagen) from 40 male and female adult olfactory rosettes. cDNA was synthesized using a random hexamer primer. Primers used for Dr V1r2, β-actin and OMP genes were previously described [16]. Primers for Dr V1r1 were the following: TGT TTC TGT CCC TGG TGC TGG T and AGC CGG AGG TCA GTG ATC AG. PCR conditions were previously described [16].5′ RACE 5′ RACE (Rapid Amplification of cDNA Ends) was performed following the Roche Applied Science RACE kit protocol. First strand cDNA synthesis, and subsequent PCR cycles were performed using the following Dr V1r1 primers: AGC ACG ATG GCA TCA GCA GGC GAG A, TCG CGG CGC ACC AGC AAG ATG AAG, and CGC AGA TGA CAG CGG TGT TAC C. Amplification conditions were: 3 min. at 95°C, followed by 34 cycles of 15 sec. at 95°C, 30 sec. at 55°C (or 60°C depending on the primer), and 1 min. at 72°C, with a final extension of 10 min. at 72°C.5′UTR Sequence Logo The 75 base pair sequences corresponding to the V1r1 5′UTR homology region of Dr, Dm, Ol, Ga, Tr, and Tn, initially identified using the Pipmaker analysis tool [30] available at http://pipmaker.bx.psu.edu/pipmaker, were aligned. A logo was generated using the sequence logo software [31] available at http://weblogo.berkeley.edu/logo.cgi.In situ hybridizations Dr V1r1 and V1r2 probes corresponded to genomic amplicons obtained with the following primers: ATG TGG TGC CGA TGG TGC TGA TGG, AGC CGT GTG TTG GAG ACG AGG ATC AG, TAT GGA CCT GTG TGT CAC, and TCA TGG AAG TCC ACA TGG CAG AAG. The conserved 5′UTR V1r1 box probe was flanked by the primers AGT GGA AAT GCA GTG TGC GC and CAA TTA CCT GAA TGA TGT GC. RNA digoxygenin labeled probes were synthesized according to the DIG RNA labeling kit supplier protocol (Roche Molecular Biochemicals). 10 µm cryosections were fixed for 20 minutes with 4% paraformaldehyde at 4°C. Hybridizations were performed in hybridization buffer (50% deionized formamide, 10% dextran sulphate, 1 µg/µl tRNA, 1×Denhardt's, 1×salt solution (2M NaCl, 100mM TRIS, 50mM NaH2PO4, 50 mM Na2HPO4, 50mM EDTA, pH 7.5)) overnight at 65°C. Primary anti-DIG antibody (1![]() Phylogenetic trees Nucleotide sequences of 24 teleost V1r sequences and mouse V1rb2, V1re4 and V1rf3 were aligned using ClustalX [32] and manually rearranged using the Bioedit alignment software (v. 7.0.5.3).The MODELTEST v3.7 program [33] was used to determine the model of DNA sequence evolution that fit our data best. The best fit model was the General Time Reversible (GTR) with a gamma shape distribution of evolutionary rates (α Evolutionary rates and selective pressures Differences in the evolutionary rates of the two paralogous V1r copies were statistically tested using the RRtree program [35]. Selective pressures acting on the V1r receptor genes were assessed using different models as implemented in the PAML software v.3.15 [36]. Couples of nested models were compared using Likelihood Ratio Tests (LRTs) to statistically determine the best model. Twice the difference of the likelihood values of the tree obtained under each model approximately follows a chi-square distribution and, together with the number of degrees of freedom (df) separating the two models, allows the calculation of the associated p-value.To evaluate a possible unequal ratio of non-synonymous vs. synonymous substitutions (ω) in the two paralogous V1r copies, we performed an LRT comparing the two-ratio model (one ω per clade) with the one-ratio model (M0) where ω is constant. Varying selective pressure along the genes was assessed by a LRT test comparing the site-specific model M3 (discrete) with 2 or 3 classes of sites displaying different ω, against the M0 model (one-ratio). The branch-site model D [37], an extension of the model M3 which allows selective pressure at one class of sites to be different in two parts of the phylogeny, was used to test for divergent selective pressure between the two paralogous V1r genes, by applying an LRT on models D vs. M3. Positive selection Positive selection was first assessed using two LRTs involving site-specific models of increasing complexity. Model M2a (which allows one group of sites to have an ω>1), was compared to model M1a (where ω cannot exceed 1). Discrete model M8 (in which the ω value can vary according to a beta distribution (between 0 and 1) and which includes an extra category of sites with ω>1), was compared to model M7 (in which ω can only vary according to the beta distribution). Branch-site models were used to identify positive selection acting on sites along specific branches of the phylogenetic tree. Model A [38] was used to identify positive selection acting on sites along specific branches of the tree. This model, which allows one category of sites to evolve under positive selection (ω2>1), was compared by LRT to that same model with an ω2 fixed at 1 (Test 2). This LRT test is the most reliable test according to Zhang et al. [38]. Significance of positively selected sites was evaluated using Bayes Empirical Bayes (BEB) calculations [39].Animals Animals were housed and handled accordingly to the guidelines and regulations of the institution and of the state of Geneva. | |||||||||||||||
Acknowledgments We thank J. Wittbrodt for medaka genomic DNA, and C. Bauer for help on the confocal imaging platform. | |||||||||||||||
Footnotes Competing Interests: The authors have declared that no competing interests exist. Funding: This study was funded by the Swiss National Science Foundation. | |||||||||||||||
References 1. Alioto, TS; Ngai, J. The odorant receptor repertoire of teleost fish. BMC Genomics. 2005;6:173. [PubMed] 2. Liberles, SD; Buck, LB. A second class of chemosensory receptors in the olfactory epithelium. Nature. 2006;442:645–650. [PubMed] 3. Gloriam, DE; Bjarnadottir, TK; Yan, YL; Postlethwait, JH; Schioth, HB, et al. The repertoire of trace amine G-protein-coupled receptors: large expansion in zebrafish. Mol Phylogenet Evol. 2005;35:470–482. [PubMed] 4. Hashiguchi, Y; Nishida, M. Evolution of vomeronasal-type odorant receptor genes in the zebrafish genome. Gene. 2005;362:19–28. [PubMed] 5. 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. [PubMed] 6. Alioto, TS; Ngai, J. The repertoire of olfactory C family G protein-coupled receptors in zebrafish: candidate chemosensory receptors for amino acids. BMC Genomics. 2006;7:309. [PubMed] 7. 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. [PubMed] 8. Boschat, C; Pelofi, C; Randin, O; Roppolo, D; Luscher, C, et al. Pheromone detection mediated by a V1r vomeronasal receptor. Nat Neurosci. 2002;5:1261–1262. [PubMed] 9. Del Punta, K; Leinders-Zufall, T; Rodriguez, I; Jukam, D; Wysocki, CJ, et al. Deficient pheromone responses in mice lacking a cluster of vomeronasal receptor genes. Nature. 2002;419:70–74. [PubMed] 10. Young, JM; Kambere, M; Trask, BJ; Lane, RP. Divergent V1R repertoires in five species: Amplification in rodents, decimation in primates, and a surprisingly small repertoire in dogs. Genome Res. 2005;15:231–240. [PubMed] 11. Lane, RP; Young, J; Newman, T; Trask, BJ. Species specificity in rodent pheromone receptor repertoires. Genome Res. 2004;14:603–608. [PubMed] 12. Grus, W; Shi, P; Zhang, Y; Zhang, J. Dramatic variation of the vomeronasal pheromone gene repertoire among five orders of placental and marsupial mammals. Proc Natl Acad Sci U S A. 2005;102:5767–5772. [PubMed] 13. Grus, WE; Zhang, J. Rapid turnover and species-specificity of vomeronasal pheromone receptor genes in mice and rats. Gene. 2004;340:303–312. [PubMed] 14. 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. [PubMed] 15. Zhang, X; Rodriguez, I; Mombaerts, P; Firestein, S. Odorant and vomeronasal receptor genes in two mouse genome assemblies. Genomics. 2004;83:802–811. [PubMed] 16. 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. [PubMed] 17. Mundy, NI; Cook, S. Positive Selection During the Diversification of Class I Vomeronasal Receptor-like (V1RL) Genes, Putative Pheromone Receptor Genes, in Human and Primate Evolution. Mol Biol Evol. 2003;20:1805–1810. [PubMed] 18. Lane, RP; Cutforth, T; Axel, R; Hood, L; Trask, BJ. Sequence analysis of mouse vomeronasal receptor gene clusters reveals common promoter motifs and a history of recent expansion. Proc Natl Acad Sci U S A. 2002;99:291–296. [PubMed] 19. Shi, P; Bielawski, JP; Yang, H; Zhang, YP. Adaptive diversification of vomeronasal receptor 1 genes in rodents. J Mol Evol. 2005;60:566–576. [PubMed] 20. Meyer, A; Schartl, M. Gene and genome duplications in vertebrates: the one-to-four (-to-eight in fish) rule and the evolution of novel gene functions. Curr Opin Cell Biol. 1999;11:699–704. [PubMed] 21. Sidow, A. Gen(om)e duplications in the evolution of early vertebrates. Curr Opin Genet Dev. 1996;6:715–722. [PubMed] 22. Suga, H; Koyanagi, M; Hoshiyama, D; Ono, K; Iwabe, N, et al. Extensive gene duplication in the early evolution of animals before the parazoan-eumetazoan split demonstrated by G proteins and protein tyrosine kinases from sponge and hydra. J Mol Evol. 1999;48:646–653. [PubMed] 23. Vandepoele, K; De Vos, W; Taylor, JS; Meyer, A; Van de Peer, Y. Major events in the genome evolution of vertebrates: paranome age and size differ considerably between ray-finned fishes and land vertebrates. Proc Natl Acad Sci U S A. 2004;101:1638–1643. [PubMed] 24. Christoffels, A; Koh, EG; Chia, JM; Brenner, S; Aparicio, S, et al. Fugu genome analysis provides evidence for a whole-genome duplication early during the evolution of ray-finned fishes. Mol Biol Evol. 2004;21:1146–1151. [PubMed] 25. Wittbrodt, J; Shima, A; Schartl, M. Medaka–a model organism from the far East. Nat Rev Genet. 2002;3:53–64. [PubMed] 26. Crnogorac-Jurcevic, T; Brown, JR; Lehrach, H; Schalkwyk, LC. Tetraodon fluviatilis, a new puffer fish model for genome studies. Genomics. 1997;41:177–184. [PubMed] 27. Chen, WJ; Orti, G; Meyer, A. Novel evolutionary relationship among four fish model systems. Trends Genet. 2004;20:424–431. [PubMed] 28. Furutani-Seiki, M; Wittbrodt, J. Medaka and zebrafish, an evolutionary twin study. Mech Dev. 2004;121:629–637. [PubMed] 29. Dulac, C; Axel, R. A novel family of genes encoding putative pheromone receptors in mammals. Cell. 1995;83:195–206. [PubMed] 30. Schwartz, S; Zha ng, Z; Frazer, KA; Smit, A; Riemer, C, et al. PipMaker--a web server for aligning two genomic DNA sequences. Genome Res. 2000;10:577–586. [PubMed] 31. Crooks, GE; Hon, G; Chandonia, JM; Brenner, SE. WebLogo: a sequence logo generator. Genome Res. 2004;14:1188–1190. [PubMed] 32. Thompson, JD; Higgins, DG; Gibson, TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22:4673–4680. [PubMed] 33. Posada, D; Crandall, KA. MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998;14:817–818. [PubMed] 34. Guindon, S; Gascuel, O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003;52:696–704. [PubMed] 35. Robinson-Rechavi, M; Huchon, D. RRTree: relative-rate tests between groups of sequences on a phylogenetic tree. Bioinformatics. 2000;16:296–297. [PubMed] 36. Yang, Z. PAML: a program package for phylogenetic analysis by maximum likelihood. Comput Appl Biosci. 1997;13:555–556. [PubMed] 37. Bielawski, JP; Yang, Z. A maximum likelihood method for detecting functional divergence at individual codon sites, with application to gene family evolution. J Mol Evol. 2004;59:121–132. [PubMed] 38. Zhang, J; Nielsen, R; Yang, Z. Evaluation of an improved branch-site likelihood method for detecting positive selection at the molecular level. Mol Biol Evol. 2005;22:2472–2479. [PubMed] 39. Yang, Z; Wong, WS; Nielsen, R. Bayes empirical bayes inference of amino acid sites under positive selection. Mol Biol Evol. 2005;22:1107–1118. [PubMed] | |||||||||||||||