The Origins of Novel Protein Interactions during Animal Opsin Evolution

Background Biologists are gaining an increased understanding of the genetic bases of phenotypic change during evolution. Nevertheless, the origins of phenotypes mediated by novel protein-protein interactions remain largely undocumented. Methodology/Principle Findings Here we analyze the evolution of opsin visual pigment proteins from the genomes of early branching animals, including a new class of opsins from Cnidaria. We combine these data with existing knowledge of the molecular basis of opsin function in a rigorous phylogenetic framework. We identify adaptive amino acid substitutions in duplicated opsin genes that correlate with a diversification of physiological pathways mediated by different protein-protein interactions. Conclusions/Significance This study documents how gene duplication events early in the history of animals followed by adaptive structural mutations increased organismal complexity by adding novel protein-protein interactions that underlie different physiological pathways. These pathways are central to vision and other photo-reactive phenotypes in most extant animals. Similar evolutionary processes may have been at work in generating other metazoan sensory systems and other physiological processes mediated by signal transduction.


INTRODUCTION
Documenting the specific genetic changes driving phenotypic evolution is a fundamental goal of current biology. Genetic changes are known to modify phenotype during evolution by altering the interactions between a protein and its ecological or biochemical environment [1][2][3][4], by modulating existing proteinprotein interactions [5], or by changing protein-DNA interactions through regulatory mutations [6][7][8][9]. However, the specific genetic changes that give rise to the evolutionary origins of novel proteinprotein interactions have rarely been documented in detail [10], but see [11].
Animal phototransduction pathways offer great opportunity for elucidating the genetic basis for evolutionary novelty for a number of reasons. First, a diversity of presumably ancient phototransduction pathways exists in animals [12]. Second, the composition of these cascades has been the subject of numerous functional biochemical studies. This is especially important because experimental demonstration that specific mutations cause phenotypic changes is often the most difficult aspect of a full documentation of the causal genetic changes driving phenotypic evolution [13]. Third, the proteins of animal phototransduction are amenable to phylogenetic study.
Fundamental to animal phototransduction pathways are the opsin visual pigment proteins, which bind to light reactive chromophores. As members of the G protein-coupled-receptor (GPCR) family [6], the various clades of animal opsins activate alternative G-proteins, resulting in three major phototransduction networks: ciliary, rhabdomeric, and G o -coupled [14][15][16]. Ciliary opsins initiate signaling through binding of a G i/t a subunit of the G-protein [15], rhabdomeric opsins utilize a G q a subunit in signaling [17] and class-specific G o a has been identified in the G ocoupled opsin signaling pathway [18][19][20]. Another class of opsins, including Retinal G-protein-coupled Receptor (RGR) and retinochromes, probably do not signal through any G-protein. Instead they are involved in the re-activation of the light reactive chromophore.
Here, we couple specific functional knowledge about opsin's role in signal transduction with new phylogenetic analyses. These analyses elucidate specific genetic changes that were likely involved in the origins of the different animal phototransduction networks, which mediate various light responses of animals. Our phylogenetically based analyses indicate a significant correlation between opsins' G-protein binding phenotypes and amino acid positions in the fourth cytoplasmic loop of opsin, especially positions homologous to 310 and 312 of bovine rhodopsin. Previous biochemical analyses demonstrate that these same amino acid positions are involved in G-protein binding function. Our additional analyses indicate that G-protein binding phenotypes likely diversified at the time of opsin gene duplication events before the origin of bilaterians, and the specific amino acid changes involved retain a pattern consistent with purifying selection.

RESULTS
Multiple opsin genes are present in cnidarians but absent in a demosponge multiple unique opsins which all lack introns. Six of these cnidarian opsins were found in public Expressed Sequence Tag (EST) databases (Table S1). Consistent with opsin status, Hydra2 is expressed in the nerve net of Hydra based on in situ hybridization and includes sequence motif hallmarks of opsin (figure 1).
While multiple opsins were found in the cnidarian genomes, our screens for opsins in trace data from the poriferan Amphimedon queenslandica did not produce any putative opsins. These screens did produce several non-opsin, rhodopsin-class GPCR genes from Amphimedon (data not shown). We were also unable to obtain opsin sequences from the trace genome data of the placozoan Trichoplax or the choanoflagellate Monosiga. Animal (type II) opsins are also unknown from numerous fungal genomes [21].
A new class of cnidarian opsins helps resolve phylogenetic relationships In our discussion of unrooted phylogenies, we refrain from using the common terms ''sister group'' and ''clade'' in favor of the terms ''adjacent group'' and ''clan'' as the former terms imply a rooting hypothesis a priori [22].
Our phylogenetic analyses reveal a new clan of opsins known only in cnidarians, which we have named ''cnidops''. Together, metazoan opsins form two major clans in unrooted analyses. One clan unites rhabdomeric and RGR/G o -coupled opsins in an adjacent group and the second consists of the newly identified cnidops family plus the ciliary opsins, which we here confirm by phylogenetic analysis to include a non-bilaterian representative (Nematostella4; figure 2).
Without non-opsin outgroups, our analyses of metazoan opsin phylogeny yielded well supported topologies. In order to root our tree, we used a combination of likelihood comparisons, reconciled tree analyses and parametric bootstrapping (summarized in figure 3 and described in Methods). The majority of available evidence indicates that opsin phylogeny is best rooted between the ciliary and non-ciliary clades (figure 4A).

Timing of origin of major opsin clades
Using our rooted opsin phylogeny, reconciled tree analyses (RTA) provide new information on the timing of origin of the various opsin clades (figure 5). RTA indicates that the ciliary and cnidops lineages originated by gene duplication at or before the cnidarianbilaterian common ancestor (Eumetazoa), approximately 600 million years ago [23,24].
Of the two ancient metazoan opsin lineages, one is represented in extant taxa by the newly described cnidops opsins of cnidaria, RGR/G o and rhabdomeric opsins. Our results concur with previous minima for the origins of RGR/G o and rhabdomeric opsins at the protostome-deuterostome common ancestor (Urbilateria) [14,25] and provide a new maximum for the origin of these clades by a gene duplication event younger than the eumetazoan ancestor (figure 5). A second ancient opsin lineage survives as the ciliary opsins and includes one cnidarian opsin (Nematostella4) and those of the vertebrate rods, cones, pineal and parapineal, and of invertebrate extra-ocular cells [25,26]. Despite finding the Nematostella gene, we did not find any Hydra genes from this ciliary opsin clade.

The early evolution of G protein binding partners
Our phylogeny provides a framework for understanding key transitions in the evolution of animal phototransduction pathways. Opsins interact with their corresponding G proteins, in part through binding at conserved sequence motifs located in opsin's fourth cytoplasmic loop [27]. We obtained clear reconstructed ancestral states for most of the residues in a conserved tri-peptide motif for the ciliary, rhabdomeric and RGR/G o nodes. For the most part, the remainder of the residues in the fourth cytoplasmic loop can be unequivocally reconstructed to the level of Dayhoff classes (i.e., C/HRK/FYW/DENQ/LIVM/GATSP) [28] ( figure 6). In addition, many of the residues of this region significantly co-vary with G protein a subunit interaction on our phylogeny. Two of the highest scoring residues in the co-variation analysis, residues 310 and 312 corresponding to positions 1 and 3 of the G protein-binding tripeptide motif [27], possessed a signature of selection using algorithms in DIVERGE2 [29] (figure 7).
We also explored the character history of G protein binding interactions by reconstructing their ancestral states across our phylogeny ( figure 6). While we cannot conclusively resolve the G protein interaction in the ancestor of rhabdomeric+G o opsins or the ancestor of ciliary+rhabdomeric+G o clades, our reconstructed ancestral G protein interaction states for ciliary, rhabdomeric and G o are highly significant (P.0.95 for each) and are congruent with previous empirical biochemical studies [17,19,30] (figure 6).

DISCUSSION
Our findings provide several new insights into the evolution of animal phototransduction cascades. First, we identify a new class of opsins that appears to be unique to cnidarians. Second, we find no evidence of any opsins in genomes of a demosponge, a parazoan, a choanoflagellate or fungi. We conclude that opsinmediated phototransduction cascades originated early in eumetazoan evolution, after this lineage diverged from the demosponge lineage. The rhabdomeric and G o -coupled cascades appear to be bilaterian innovations. Third, the origins of novel phototransduction cascades, as defined by the evolution of alternative opsin-G protein interactions, co-occur on our phylogeny with the gene duplications that have given rise to the major opsin classes. Fourth, the specific mutational events that have lead to the functional diversification of the major opsin lineages are suggested by covariation analyses and tests for selection. These results demonstrate how gene duplications, together with structural mutations that lead to novel protein-protein interactions, can contribute to the evolution of novel physiological traits. Each of these points is discussed further below.

The timing and relationships of animal opsins
The inter-relationships of animal opsins reported here are based on thorough phylogenetic analyses that included standard phylogenetic tree estimation, comprehensive likelihood comparisons, RTA, and parametric bootstrapping (figure 3). Taken together, the majority of evidence favors the hypothesis depicted in figure 4A, although the trees in figures 4B and C remain viable alternative hypotheses. However, because these latter hypotheses require additional gene duplication and loss events to explain the rooted topology, and because parametric bootstrapping simulations indicate that long branch lengths can contribute to the erroneous resolution of specifically these topologies, we favor the former hypothesis (figure 4A). Given present sampling and the assumption of monophyly of major opsin clades, any other phylogenetic hypotheses are significantly rejected by the data. Even if figure 4A does not represent the true history of opsin, the two viable alternative hypotheses (figures 4B and C) do not drastically alter our interpretation of the timing of origins of animal opsin in RTA analysis. Irrespective of which hypothesis we consider, animal opsins originate at the same node (Eumetazoa). If we assume the hypothesis depicted in figure 4B our conclusion that RGR/G o -coupled and rhabdomeric opsins are bilaterian synapomorphies also remains intact, but requires an additional loss of cnidops prior to Bilateria. The latter is also true of the topology shown in figure 4C.
A maximum date for the origin of animal opsins (when a GPCR first bound a light reactive chromophore) is equivocal because we did not obtain any opsins in the genome sequence of an earlier branching animal, the demosponge Amphimedon queenslandica, or from any non-animal genomes. While the larvae of Amphimedon are capable of phototaxis, studies of spectral sensitivity have suggested the activities of a flavin-based photoreceptor in this behavior [31], as opposed to retinal as used by opsins. Further, we uncovered no animal (type II) opsins from bioinformatic screens of the placozoan Trichoplax or the choanoflagelate Monosiga genomes and no true homolog of type II opsins have been described outside of Metazoa [32]. While we cannot rule out an earlier opsin origin coupled with loss in lineages leading to the species with genome sequences, we strongly favor the hypothesis that animal opsins were present in the eumetazoan ancestor of cnidarians and bilaterians, but not at earlier nodes in eukaryote phylogeny. Since porifera, Monosiga, and fungi seem to lack opsins, yet possess rhodopsin-class GPCR proteins without the Lysine residue to bind a light reactive chromophore, one possibility is that a rhodopsin-class GPCR gained the ability to bind a chromophore by substitution to a Lysine residue during animal evolution.
Our results contrast with previous studies on opsin relationships that placed ciliary and G o -coupled opsins as sister groups [14], or adjacent groups [15], and those that placed ciliary and rhabdomeric opsins as sister groups [12]. However, in lacking data from nonbilaterian animals, these previous works were able to include only three of the four animal opsin lineages discussed here. The new clade of opsins from the Cnidarian, ''cnidops'', has aided in the resolution of the overall opsin phylogeny, especially because non-opsin outgroups seem to destabilize ingroup topology.
The presence of an opsin from the cnidarian Nematostella that in our analyses groups within the ciliary clade has been foreshadowed by previous morphological descriptions of ciliary photoreceptors [33,34], but no previous opsin data from the group have confirmed this hypothesis. Our ability to find only cnidops and not ciliary opsins in Hydra probably indicates that ciliary clade opsins were lost during this species' evolutionary history.

The genetic basis for opsin-G protein interactions
Multiple lines of evidence support the hypothesis that amino acid substitutions in the 4th cytoplasmic loop of duplicated opsins were involved in the origins of novel opsin-G protein interactions. First, specific amino acid states significantly co-vary with G-protein activation phenotypes during opsin phylogeny (figure 7). Because we used a Bayesian approach that estimated the posterior probability of all possible trees, the estimated correlation accounts for uncertainty in opsin phylogeny. The highest p-values in these correlation analyses correspond to residues formerly demonstrated by site-directed mutagenesis to mediate opsin-G protein interaction in ciliary opsins (i.e., residues 310 and 312 [27]).
These amino acid changes retain the signature of selection, as given by site-specific shifts in evolutionary rates, a statistical test that uncovers significant conservation within gene clades compared to diversification between gene clades [29] (figure 7).  Table 1. Orange X indicates that tree had significantly lower likelihood in opsin-only dataset. Red X indicates tree had significantly lower likelihood in opsin+outgroup dataset. Yellow X indicates that tree implies additional gene duplication/loss events compared to minimum (tree 1 is minimum with 2 duplications 0 losses implied). Grey X indicates tree with cnidops as earliest branching opsin group-a result inferred in parametric bootstrap analyses, which incorrectly grouped cnidops+outgroup opsins because of long branch-attraction (see supplementary methods). Cil = ciliary; Rh = rhabdomeric; RG = Go/RGR. doi:10.1371/journal.pone.0001054.g003 Although this pattern cannot be used to determine if the mutations were originally fixed by natural selection, conservation after an ancient diversification indicates a long period of stasis, consistent with purifying selection on adaptive, alternative amino acid states. The combination of correlation, functional tests, and a pattern consistent with purifying selection strongly suggests that substitutions in residues homologous to bovine opsin 310 and 312 contributed to the origin of novel phototransduction pathways with alternative G-protein binding partners.
The pattern of early divergence followed by conservation is also evident in our estimates of ancestral opsin-G protein binding phenotypes. Using available data, we were able to unequivocally resolve the ancestral G-protein usage phenotypes for the ciliary, rhabdomeric and RGR/G o -coupled opsin clades (figure 6). That we can resolve the G protein binding specificities at the base of each of these clades, but not at deeper nodes, suggests an early diversification of G protein binding specificities, coupled with relative stasis of protein interactions once the major opsin clades were established.
Because we find that opsins were present in extant eumetazoans but not in animals representing earlier branching lineages (i.e. demosponges), the evolution of novel phototransduction cascades in animals probably proceeded by evolving specialized proteinprotein interactions between new opsin paralogs and existing G protein signaling intermediaries, an evolutionary process previously described as molecular exploitation [4]. Each of the G a subunit paralogs known to interact with opsins in bilaterians has also been described from plants, fungi, protists, cnidarians and sponges indicating a long pre-animal history for the different G-proteins [35,36].
While our inclusion of new data from the Cnidaria has greatly enhanced our ability to reconstruct opsin phylogeny, functional or phenotypic data from the newly reported cnidops clade do not exist. Further studies on the biochemistry of signaling by cnidops are required if we are to better understand the history of opsin-G

Implications for early animal sensory systems
Our finding of an adaptive, early partitioning of opsin-mediated signaling into discrete functional classes in the course of their evolution provides insight into the organismal properties of the ancestor of cnidarian and bilaterian animals. The activities of ciliary, rhabdomeric and G o -coupled opsins are known to exist in close proximity [37], if not the same cell [30], in numerous taxa. The early, adaptive evolution of the different photosensory pathways suggested by our analyses may be an indication of the functional necessity to transduse these signaling cascades through non-overlapping channels in the eumetazoan ancestor, an animal that by our analysis possessed ancestral ciliary and ancestral cnidops+rhabdomeric+RGR/G o -clade opsins. Such early functional requirements could have included the need to delineate light information used for vision from that used to entrain the circadian clock cycle [38] or to differentiate between light of differing wavelengths, two core processes accomplished by specific cell types in bilaterian animals [12]. The proposed relationship between gene duplication and biological complexity has been invoked and re-invoked at various points throughout the past century of biological investigation [39], and was most famously articulated by Ohno [40]. Despite the great level of interest in models for the evolution of complexity based on gene duplication, data demonstrating how gene duplication can be implicated in the origins of novel protein-protein interactions are not common. Here we begin to fill this gap and show that a group of familiar signaling pathways, animal phototransduction cascades, provides a useful model for the study of the origins of novel phenotypes.

Data mining
Publicly available trace genome sequence data from Hydra magnipapillata, Nematostella vectensis and Amphimedon queenslandica (http://www.ncbi.nlm.nih.gov/Traces) were subjected to tblastn [41] searches using a wide range of homologous opsin proteins from bilaterian animals (Table S1). We implemented an in-house assembly pipeline, where BLAST hits were extended to their largest possible contigs under stringent parameters, using only clean trace data as determined using PHRAP [42] and an 85% identity cut off for joining trace fragments. Genes were predicted from genomic contigs using GenomeScan [43]. A combination of BLAST searches, phylogenetic analyses and sequence motif analyses were used to establish the authenticity of opsin sequences recovered from genome analyses. In particular, the residue homologous to the bovine lysine 296, which binds a light reactive chromophore [44], has been used to validate opsin identity (figure 1). We corroborated the existence of our predicted proteins by using EST databases at either NCBI (http://www.ncbi.nlm. nih.gov/), for Hydra, or the Joint Genome Institute, for Nematostella (http://www.jgi.doe.gov/). Homologous proteins and outgroup data from bilaterian animals used in phylogenetic analyses were also obtained from NCBI (Table S1).

Animal culture and in situ hybridization
Hydra magnipapillata (UC Irvine strain) were reared in the laboratory under standard conditions at 18uC. The expression of one opsin, HM2, was investigated in Hydra magnipapillata using in situ hybridization following Grens et al. [45]. The following modifications were made: A pCRII vector (Invitrogen) containing HM2 opsin cDNA (,750 bp) was used as template for probe synthesis. DIG-labeled probes were synthesized using T3 and T7 polymerase (Roche). Both sense and anti-sense probes (0.3 mg/ml) were hybridized at 55uC for 48 to 60 hours. Sense probes produced no detectable signal.

Phylogenetic approaches
For phylogenetic analyses, only the 7-transmembrane region including intervening inter-and extra-cellular domains was included (330 amino acids), as it was difficult to ascertain homology of N-and C-termini due to sequence length variation and lack of conservation across genes. Taxon/gene selection was done in a manner that enhanced both the taxonomic representation across the Metazoa and the inclusion of well-studied subfamilies (i.e., arthropod Rh opsins and vertebrate visual opsins), but also allowed for thorough computational analyses. Protein sequences were aligned using T-COFFEE under default parameters [46] and alignment manipulations were done using Seaview [47]. Phylogenetic analyses were conducted using unweighted Maximum Parsimony (MP) implemented in PAUP* 4.0b10 [48], Maximum Likelihood (ML) implemented in PHYML v2.4.4 [49], and Bayesian Markov Chain Monte Carlo (BMCMC) implemented in MrBayes 3.1 [50]. Support for internal nodes was assessed with 1000 bootstrap replicates for MP and ML, and posterior probability for BMCMC analyses. ML and BMCMC approaches assumed best-fit models of protein evolution, as determined using ProtTest [51]. For ML tree calculations, both the proportion of invariant sites and the a parameter of the gamma distribution were estimated and both tree topology and branch lengths were optimized in PHYML. Bayesian analyses were conducted for 20 million generations using default priors and heating parameters for each of four chains. For the monitoring of progress in BMCMC runs, we assumed convergence of the Markov Chains when standard deviation of split frequencies (SDSF) fell below 0.01, indicating that the two independent runs resulted in similar phylogenies [50]. Following the BMCMC runs, burnin was assessed using Tracer [52] and by the SDSF of the two independent BMCMC runs. For unrooted analyses, the initial 11 million generations with an SDSF of 0.0015 or greater were discarded. For rooted analyses the first 12 million generations with an SDSF of 0.0024 or greater were removed prior to consensus tree calculation.
The evolutionary relationships between members of the larger GPCR class of proteins remain poorly understood. Because of this, our outgroup selection in rooted analyses was empirical. We selected the shortest-branched sequences from a pool of rhodopsin class GPCRs [32] from a variety of taxa (Table S1) after nonbootstrapped ML and shorter (i.e., 3-5 million steps) BMCMC analyses.

Tests of phylogenetic hypotheses
In order to address the lack of resolution among the major opsin clades encountered in phylogenetic analyses that included outgroups (figure S2), we tested the significance of a wide range of competing phylogenetic hypotheses for metazoan opsin. First, we assumed monophyly of each of the four major opsin clades (cnidarian opsin, ciliary, RGR/G o and rhabdomeric). This assumption is supported by the results of unrooted MP, ML and BMCMC analyses (figure S2) and by previous studies on opsin phylogeny [12,15]. Our topology contains four opsin clades that can be rearranged in 15 possible rooted binary trees. Second, we calculated the likelihood for each possible tree using the ''resolve multifurcations'' function in TREEFINDER [53]. Finally, we used CODEML, as included in PAML [54], to assess the significance of each of these trees. These results were further analyzed using CONSEL [55]. A similar approach was used to assess alternative topologies in unrooted trees where having four major clades involves three possible unrooted topologies (figure 4).

Parametric bootstrapping
Given the ambiguous results obtained when including non-opsin outgroups, we suspected that long branch-attraction (LBA) artifacts [56], where the cnidops sequences were pulled to the base of the tree by long branch outgroups, were confounding our analyses. We tested this hypothesis by using Huelsenbeck's [57] method of parametric bootstrapping. Here, a tree including a nonopsin outgroup was constrained to reflect the best topology consistent with Tree 1B in figure 3. Assuming this tree and the best-fit model of opsin molecular evolution from prot-test analysis [51] , we simulated 100 data sets using SeqGen [58]. Each of these simulated data sets was subjected to ML and MP analyses. Since the data were simulated, we expected to recover the assumed phylogeny unless the topology/model combination is sensitive to artifacts like long branch-attraction [57]. For ML analyses we found that the simulated cnidarian genes are incorrectly attracted to the base of the tree in 45% of replicates, whereas the correct, simulated topology was recovered in only 35% of the replicates. In MP analyses, the cnidarian opsin genes were incorrectly attracted Figure 7. Summary of co-variation analyses and tests for selection. Previously described G protein binding states, together with amino acid sequences of the 4 th cytoplasmic loop regions from our opsin dataset, were used to test the hypothesis that the two character sets co-vary across our opsin phylogeny. Co-variation was assessed using mutual information content (MIC). Shown here is a representation of these results. Predictive P-values (Bottom axis) based on MIC for co-variation between each residue in the 4 th cytoplasimic loop and their respective G protein interactions (Right axis) are shown. These analyses were conducted using Bayesian mutational simulation mapping in SIMMAP [64]. We also tested possibility that these residues had evolved under a selective regime using the criterion of site specific rate heterogeneity as implemented in DIVERGE2 [29]. The two highest scoring residues from co-variation analyses (310 and 312) also retain the signature of selection (asterisk). See table S1 for citations. P = predictive P-value, M = the M statistic given by MIC. doi:10.1371/journal.pone.0001054.g007 to the outgroup in 49% of the replicates, while the simulated phylogeny was again recovered in only 35% of cases. These simulations show that analyses including opsin outgroups are sensitive to artifacts that induce erroneous topology estimates. Our results are similar to those reported for analyses of strepsipteran relationships, a well-known case of LBA [57].

Reconciled tree analysis
The number of gene duplication and loss events that are implied by a given rooting hypothesis can be used to assess the chances that a given hypothesis is tenable [59]. We compared the number of duplications and losses implied by each possible root position of our unrooted phylogeny using NOTUNG [60,61]. For these analyses we assumed a conservative species-level phylogeny for the major taxa included in our analyses: Hydrozoa, Anthozoa, Cephalochordata, Urochordata, Vertebrata, Annelida, Mollusca, and Arthropoda [62] .

Character mapping and tests of co-variance
We scored G-protein interaction phenotypes for the opsins in our phylogenetic analysis as a discrete state character. Although these phenotypes are unknown for many opsins, especially those from non-model organisms, our analysis includes interaction data for 48% of the opsin sequences represented on our phylogeny (figures 6, S1).
Ancestral state reconstructions and tests for correlated character evolution were conducted using Bayesian mutational mapping [63] as implemented in SIMMAP [64]. Reconstructions were integrated over 18,000 trees that were sampled from our unrooted BMCMC analyses. Ancestral state reconstructions were conducted with 10 realizations sampled from the prior and 10 realizations sampled for each tree. We set an equal prior on the bias parameter but did not abstract a prior for the rate parameter, instead using branch lengths from BMCMC trees as relative rates. As there is no a priori reason to assume different prior parameters, we chose this model, as it is the least parameterized model for morphological evolution possible in SIMMAP [64]. We also tested a range of alternative prior settings and obtained similar P-values in these analyses (results not shown). Tests for selection were conducted in DIVERGE2 [29] using our amino acid alignment (see supplemental materials), rooted phylogeny and the structure of bovine opsin previously calculated by x-ray crystallography [65]. Our estimate for site specific rate shift (functional divergence) between ciliary and cnidops+RGR/G o +Rhabdomeric was significantly larger than 0: h 1 = 5.8260.05. Site-specific posterior probabilities for residues 310 and 312 were 0.45 and 0.37 respectively.