Molecular evolution of color vision of zebra finch

We have isolated and sequenced the RH1(Tg), RH2(Tg), SWS2(Tg), and LWS(Tg) opsin cDNAs from zebra finch retinas. Upon binding to 11-cis-retinal, these opsins regenerate the corresponding photosensitive molecules, visual pigments. The absorption spectra of visual pigments have a broad bell shape, with the peak being called lambda(max). Previously, SWS1(Tg) opsin cDNA was isolated from zebra finch retinal RNA, expressed in cultured COS1 cells, reconstituted with 11-cis-retinal, and the lambda(max) of the resulting visual pigment was shown to be 359nm. Here, the lambda(max) values of the RH1(Tg), RH2(Tg), SWS2(Tg), and LWS(Tg) pigments are determined to be 501, 505, 440, and 560nm, respectively. Molecular evolutionary analyses suggest that specific amino acid replacements in the SWS1 and SWS2 pigments, resulting from accelerated evolution, must have been responsible for their functional divergences among the avian pigments.     

Gene duplication and spectral diversification of cone visual pigments of zebrafish

Zebrafish is becoming a powerful animal model for the study of vision but the genomic organization and variation of its visual opsins have not been fully characterized. We show here that zebrafish has two red (LWS-1 and LWS-2), four green (RH2-1, RH2-2, RH2-3, and RH2-4), and single blue (SWS2) and ultraviolet (SWS1) opsin genes in the genome, among which LWS-2, RH2-2, and RH2-3 are novel. SWS2, LWS-1, and LWS-2 are located in tandem and RH2-1, RH2-2, RH2-3, and RH2-4 form another tandem gene cluster. The peak absorption spectra (lambdamax) of the reconstituted photopigments from the opsin cDNAs differed markedly among them: 558 nm (LWS-1), 548 nm (LWS-2), 467 nm (RH2-1), 476 nm (RH2-2), 488 nm (RH2-3), 505 nm (RH2-4), 355 nm (SWS1), 416 nm (SWS2), and 501 nm (RH1, rod opsin). The quantitative RT-PCR revealed a considerable difference among the opsin genes in the expression level in the retina. The expression of the two red opsin genes and of three green opsin genes, RH2-1, RH2-3, and RH2-4, is significantly lower than that of RH2-2, SWS1, and SWS2. These findings must contribute to our comprehensive understanding of visual capabilities of zebrafish and the evolution of the fish visual system and should become a basis of further studies on expression and developmental regulation of the opsin genes.     

Genetic basis of spectral tuning in the violet-sensitive visual pigment of African clawed frog, Xenopus laevis

Ultraviolet (UV) and violet vision in vertebrates is mediated by UV and violet visual pigments that absorb light maximally (lambdamax) at approximately 360 and 390-440 nm, respectively. So far, a total of 11 amino acid sites only in transmembrane (TM) helices I-III are known to be involved in the functional differentiation of these short wavelength-sensitive type 1 (SWS1) pigments. Here, we have constructed chimeric pigments between the violet pigment of African clawed frog (Xenopus laevis) and its ancestral UV pigment. The results show that not only are the absorption spectra of these pigments modulated strongly by amino acids in TM I-VII, but also, for unknown reasons, the overall effect of amino acid changes in TM IV-VII on the lambdamax-shift is abolished. The spectral tuning of the contemporary frog pigment is explained by amino acid replacements F86M, V91I, T93P, V109A, E113D, L116V, and S118T, in which V91I and V109A are previously unknown, increasing the total number of critical amino acid sites that are involved in the spectral tuning of SWS1 pigments in vertebrates to 13.     

Importance of the N-terminal sequence of the extrinsic 23 kDa polypeptide in Photosystem II in ion retention in oxygen evolution

The function of the extrinsic 23 kDa polypeptide (OEC23) in Photosystem II (PS II) is to retain Ca(2+) and Cl(-) during the S-state turnover of the Mn cluster in photosynthetic oxygen evolution. Recombinant OEC23s from several plant species were produced in Escherichia coli to characterize the molecular mechanism of the OEC23 function then used in reconstitution experiments. One tobacco isoform, OEC23 (2AF), had much less oxygen-evolving activity than the spinach and cucumber OEC23s when PS II activities were reconstituted in salt-washed spinach PS II particles. The fact that the OEC23s had similar binding affinities for PS II particles suggests different ion-retention capacities for the individual OEC23s: The chimeric OEC23s produced between spinach and 2AF and those produced between cucumber and 2AF show that 58 N-terminal amino acid residues are important for PS II activity. Further dissection of the sequence and site-directed mutagenesis indicated the importance of 20 N-terminal amino acid residues for activity, in particular the asparagine at the 15th position. In spinach the N15D mutation decreased PS II activity, whereas in 2AF the D15N mutation increased it. This shows the importance of the N-terminal sequence of OEC23 in ion retention during the water-splitting process.     

Parallel and Convergent Evolution of the Dim-Light Vision Gene RH1 in Bats (Order: Chiroptera)

Rhodopsin, encoded by the gene Rhodopsin (RH1), is extremely sensitive to light, and is responsible for dim-light vision. Bats are nocturnal mammals that inhabit poor light environments. Megabats (Old-World fruit bats) generally have well-developed eyes, while microbats (insectivorous bats) have developed echolocation and in general their eyes were degraded, however, dramatic differences in the eyes, and their reliance on vision, exist in this group. In this study, we examined the rod opsin gene (RH1), and compared its evolution to that of two cone opsin genes (SWS1 and M/LWS). While phylogenetic reconstruction with the cone opsin genes SWS1 and M/LWS generated a species tree in accord with expectations, the RH1 gene tree united Pteropodidae (Old-World fruit bats) and Yangochiroptera, with very high bootstrap values, suggesting the possibility of convergent evolution. The hypothesis of convergent evolution was further supported when nonsynonymous sites or amino acid sequences were used to construct phylogenies. Reconstructed RH1 sequences at internal nodes of the bat species phylogeny showed that: (1) Old-World fruit bats share an amino acid change (S270G) with the tomb bat; (2) Miniopterus share two amino acid changes (V104I, M183L) with Rhinolophoidea; (3) the amino acid replacement I123V occurred independently on four branches, and the replacements L99M, L266V and I286V occurred each on two branches. The multiple parallel amino acid replacements that occurred in the evolution of bat RH1 suggest the possibility of multiple convergences of their ecological specialization (i.e., various photic environments) during adaptation for the nocturnal lifestyle, and suggest that further attention is needed on the study of the ecology and behavior of bats.     

Visual pigment evolution in Characiformes: The dynamic interplay of teleost whole‐genome duplication,surviving opsins and spectral tuning

Vision represents an excellent model for studying adaptation, given the genotype‐to‐phenotype map that has been characterized in a number of taxa. Fish possess a diverse range of visual sensitivities and adaptations to underwater light, making them an excellent group to study visual system evolution. In particular, some speciose but understudied lineages can provide a unique opportunity to better understand aspects of visual system evolution such as opsin gene duplication and neofunctionalization. In this study, we showcase the visual system evolution of neotropical Characiformes and the spectral tuning mechanisms they exhibit to modulate their visual sensitivities. Such mechanisms include gene duplications and losses, gene conversion, opsin amino acid sequence and expression variation, and A₁/A₂‐chromophore shifts. The Characiforms we studied utilize three cone opsin classes (SWS2, RH2, LWS) and a rod opsin (RH1). However, the characiform's entire opsin gene repertoire is a product of dynamic evolution by opsin gene loss (SWS1, RH2) and duplication (LWS, RH1). The LWS‐ and RH1‐duplicates originated from a teleost specific whole‐genome duplication as well as characiform‐specific duplication events. Both LWS‐opsins exhibit gene conversion and, through substitutions in key tuning sites, one of the LWS‐paralogues has acquired spectral sensitivity to green light. These sequence changes suggest reversion and parallel evolution of key tuning sites. Furthermore, characiforms' colour vision is based on the expression of both LWS‐paralogues and SWS2. Finally, we found interspecific and intraspecific variation in A₁/A₂‐chromophores proportions, correlating with the light environment. These multiple mechanisms may be a result of the diverse visual environments where Characiformes have evolved.     

Cone visual pigments in two marsupial species: the fat-tailed dunnart (Sminthopsis crassicaudata) and the honey possum (Tarsipes rostratus)

Uniquely for non-primate mammals, three classes of cone photoreceptors have been previously identified by microspectrophotometry in two marsupial species: the polyprotodont fat-tailed dunnart (Sminthopsis crassicaudata) and the diprotodont honey possum (Tarsipes rostratus). This report focuses on the genetic basis for these three pigments. Two cone pigments were amplified from retinal cDNA of both species and identified by phylogenetics as members of the short wavelength-sensitive 1 (SWS1) and long wavelength-sensitive (LWS) opsin classes. In vitro expression of the two sequences from the fat-tailed dunnart confirmed the peak absorbances at 363 nm in the UV for the SWS1 pigment and 533 nm for the LWS pigment. No additional expressed cone opsin sequences that could account for the middle wavelength cones could be amplified. However, amplification from the fat-tailed dunnart genomic DNA with RH1 (rod) opsin primer pairs identified two genes with identical coding regions but sequence differences in introns 2 and 3. Uniquely therefore for a mammal, the fat-tailed dunnart has two copies of an RH1 opsin gene. This raises the possibility that the middle wavelength cones express a rod rather than a cone pigment.     

Multiple origins of the green-sensitive opsin genes in fish

Vertebrate opsins are divided into four major groups: RH1 (rhodopsins), RH2 (rhodopsinlike with various absorption sensitivities), SWS (short-wavelength sensitive), and LWS/MWS (long and middle-wavelength sensitive) groups. The green opsin genes (g101 _Af and g101 _Af ) in a Mexican characin Astyanax fasciatus belong to the LWS/MWS group, whereas those in goldfish belong to the RH2 group (Yokoyama 1994, Mol Biol Evol 11:32–39). A newly isolated opsin gene (rh11 _Af ) from A. fasciatus contains five exons and four introns, spanning 4.2 kilobases from start to stop codons. This gene is most closely related to the two green opsin genes of goldfish and belongs to the RH2 group. In the LWS/MWS group, gene duplication of the ancestral red and green opsin genes predates the speciation between A. fasciatus and goldfish, suggesting that goldfish also has an additional gene which is orthologous to g101 _Af and g103 _Af .Correspondence to: S. Yokoyama     

Vertebrate opsins belonging to different classes vary in constitutively active properties resulting from salt-bridge mutations

Vertebrate opsins are classified into one of five classes on the basis of amino acid similarity. These classes are short wavelength sensitive 1 and 2 (SWS1, SWS2), medium/long wavelength sensitive (M/LWS), and rod opsin like 1 and 2 (RH1, RH2). In bovine rod opsin (RH1), two critical amino acids form a salt bridge in the apoprotein that maintains the opsin in an inactive state. These residues are K296, which functions as the chromophore binding site, and E113, which functions as the counterion to the protonated Schiff base. Corresponding residues in each of the other vertebrate opsin classes are believed to play similar roles. Previous reports have demonstrated that mutations in these critical residues result in constitutive activation of transducin by RH1 class opsins in the absence of chromophore. Additionally, recent reports have shown that an E113Q mutation in SWS1 opsin is constitutively active. Here we ask if the other classes of vertebrate opsins maintain activation characteristics similar to that of bovine RH1 opsin. We approach this question by making the corresponding substitutions which disrupt the K296/E113 salt bridge in opsins belonging to the other vertebrate opsin classes. The mutant opsins are tested for their ability to constitutively activate bovine transducin. We demonstrate that mutations disrupting this key salt bridge produce constitutive activation in all classes. However, the mutant opsins differ in their ability to be quenched in the dark state by the addition of chromophore as well as in their level of constitutive activation. The differences in constitutive activation profiles suggest that structural differences exist among the opsin classes that may translate into a difference in activation properties.     

Genetic analyses of visual pigments of the pigeon (Columba livia)

We isolated five classes of retinal opsin genes rh1(Cl), rh2(Cl), sws1(Cl), sws2(Cl), and lws(Cl) from the pigeon; these encode RH1(Cl), RH2(Cl), SWS1(Cl), SWS2(Cl), and LWS(Cl) opsins, respectively. Upon binding to 11-cis-retinal, these opsins regenerate the corresponding photosensitive molecules, visual pigments. The absorbance spectra of visual pigments have a broad bell shape with the peak, being called lambdamax. Previously, the SWS1(Cl) opsin cDNA was isolated from the pigeon retinal RNA, expressed in cultured COS1 cells, reconstituted with 11-cis-retinal, and the lambdamax of the resulting SWS1(Cl) pigment was shown to be 393 nm. In this article, using the same methods, the lambdamax values of RH1(Cl), RH2(Cl), SWS2(Cl), and LWS(Cl) pigments were determined to be 502, 503, 448, and 559 nm, respectively. The pigeon is also known for its UV vision, detecting light at 320-380 nm. Being the only pigments that absorb light below 400 nm, the SWS1(Cl) pigments must mediate its UV vision. We also determined that a nonretinal P(Cl) pigment in the pineal gland of the pigeon has a lambdamax value at 481 nm.     

Salmonid Opsin Sequences Undergo Positive Selection and Indicate an Alternate Evolutionary Relationship in Oncorhynchus

Positive selection can be demonstrated by statistical analysis when non-synonymous nucleotide substitutions occur more frequently than synonymous substitutions (dN>dS). This pattern of sequence evolution has been observed in the rhodopsin gene of cichlids. Mutations in opsin genes resulting in amino acid (AA) replacement appear to be associated with the evolution of specific color patterns and the evolution of courtship behaviors. Within fish, AA replacements in opsin proteins have improved vision at great depths and have occurred in deep-sea species. Salmonids experience diverse photic environments during their life history. Furthermore, sexual selection has resulted in species-specific male and female coloration during spawning. To look for evidence of positive selection in salmonid opsins, we sequenced the RH1, RH2, LWS, SWS1, and SWS2 genes from six Pacific salmon species as well as the Atlantic salmon. These salmonids include landlocked and migratory species and species that vary in their coloration during spawning. In each opsin gene comparison from all species sampled, traditional dN:dS analysis did not indicate positive selection. However, the more sensitive Creevey–McInerney statistical analysis indicates that RH1 and RH2 experienced positive selection early in the evolution and speciation of salmonids.     

The significant conservative and variable regions of the homologous protein sequences.

A method of identification of significant conservative and variable regions in homologous protein sequences is presented. A set of aligned homologous sequences is divided into two groups consisting of m and n most related sequences. Each pair of sequences from different group is compared using unitary similarity matrix. The superposition of pairwise comparisons scanned by a window of 10 amino acid residues gives intergroup local variability profile (VP). Area S of the figure between the VP and its mean value line is compared with averaged area S(r) of 1000 VPs of artificial homologous protein families. The difference (S-S(r)) given in standard deviation units sigma r is believed to be the amino acid substitution overall irregularity along the homologous protein sequences OI = (S-S(r))/sigma r. If OI greater than 2, the real VP extrema containing the surplus of area S-(S(r) + 2 sigma r) are cut off. The cut off stretches are likely to be significant conservative and variable regions. The significant conservative and variable regions of six homologous sequence families (phospholipases A2, cytochromes b, alpha-subunits of Na, K-ATPase, L- and M-subunits of photosynthetic bacteria photoreaction centre and human rhodopsins) were identified. It was shown that for artificial homologous protein sequences derived by k-fold lengthening of natural proteins the OI value rises as square root of k. To compare the degree of substitution irregularity in homologous protein sequence families of different length L the value of standard substitution overall irregularity for L = 250 is proposed.     

Rapid light-induced shifts in opsin expression: finding new opsins, discerning mechanisms of change, and implications for visual sensitivity

Light-induced shifts in cone frequency and opsin expression occur in many aquatic species. Yet little is known about how quickly animals can alter opsin expression and, thereby, track their visual environments. Similarly, little is known about whether adult animals can alter opsin expression or whether shifts in opsin expression are limited to critical developmental windows. We took adult wild-caught bluefin killifish (Lucania goodei) from three different lighting environments (spring, swamp and variable), placed them under two different lighting treatments (clear vs. tea-stained water) and monitored opsin expression over 4 weeks. We measured opsin expression for five previously described opsins (SWS1, SWS2B, SWS2A, RH2-1 and LWS) as well as RH2-2 which we discovered via 454 sequencing. We used two different metrics of opsin expression. We measured expression of each opsin relative to a housekeeping gene and the proportional expression of each opsin relative to the total pool of opsins. Population and lighting environment had large effects on opsin expression which were present at the earliest time points indicating rapid shifts in expression. The two measures of expression produced radically different patterns. Proportional measures indicated large effects of light on SWS1 expression, whereas relative measures indicated no such effect. Instead, light had large effects on the relative expression of SWS2B, RH2-2, RH2-1 and LWS. We suggest that proportional measures of opsin expression are best for making inferences about colour vision, but that measures relative to a housekeeping gene are better for making conclusions about which opsins are differentially regulated.     

Divergent selection for opsin gene variation in guppy (Poecilia reticulata) populations of Trinidad and Tobago

The guppy is known to exhibit remarkable interindividual variations in spectral sensitivity of middle to long wavelength-sensitive (M/LWS) cone photoreceptor cells. The guppy has four M/LWS-type opsin genes (LWS-1, LWS-2, LWS-3 and LWS-4) that are considered to be responsible for this sensory variation. However, the allelic variation of the opsin genes, particularly in terms of their absorption spectrum, has not been explored in wild populations. Thus, we examined nucleotide variations in the four M/LWS opsin genes as well as blue-sensitive SWS2-B and ultraviolet-sensitive SWS1 opsin genes for comparison and seven non-opsin nuclear loci as reference genes in 10 guppy populations from various light environments in Trinidad and Tobago. For the first time, we discovered a potential spectral variation (180 Ser/Ala) in LWS-1 that differed at an amino acid site known to affect the absorption spectra of opsins. Based on a coalescent simulation of the nucleotide variation of the reference genes, we showed that the interpopulation genetic differentiation of two opsin genes was significantly larger than the neutral expectation. Furthermore, this genetic differentiation was significantly related to differences in dissolved oxygen (DO) level, and it was not explained by the spatial distance between populations. The DO levels are correlated with eutrophication that possibly affects the color of aquatic environments. These results suggest that the population diversity of opsin genes is significantly driven by natural selection and that the guppy could adapt to various light environments through color vision changes.     

Amino acid replacements and wavelength absorption of visual pigments in vertebrates

An important unanswered question in phototransduction is how visual pigments (VPs) regulate their wavelength of maximal absorption (lambda max). By constructing the evolutionary tree for 28 opsins with known lambda max values, we can identify the times and directions of lambda max shift of different VPs. A total of 55 amino acid changes are shown to correlate with the directions of lambda max shift and might have been important in determining lambda max of a VP. Among these, three amino acid changes are already proven to be responsible in modifying the green-sensitive VP to the red-sensitive VP. The present evolutionary analysis opens a new direction in understanding the mechanism for the regulation of wavelength absorption by a VP and, more generally, in studying molecular mechanism involved in adaptive evolution.      

Functional characterization of visual opsin repertoire in Medaka (Oryzias latipes)

A variety of visual pigment repertoires present in fish species is believed due to the great variation under the water of light environment. A complete set of visual opsin genes has been isolated and characterized for absorption spectra and expression in the retina only in zebrafish. Medaka (Oryzias latipes) is a fish species phylogenetically distant from zebrafish and has served as an important vertebrate model system in molecular and developmental genetics. We previously isolated a medaka rod opsin gene (RH1). In the present study we isolated all the cone opsin genes of medaka by genome screening of a lambda-phage and bacterial artificial chromosome (BAC) libraries. The medaka genome contains two red, LWS-A and LWS-B, three green, RH2-A, RH2-B and RH2-C, and two blue, SWS2-A and SWS2-B, subtype opsin genes as well as a single-copy of the ultraviolet, SWS1, opsin gene. Previously only one gene was believed present for each opsin type as reported in a cDNA-based study. These subtype opsin genes are closely linked and must be the products of local gene duplications but not of a genome-wide duplication. Peak absorption spectra (lambda(max)) of the reconstituted photopigments with 11-cis retinal varied greatly among the three green opsins, 452 nm for RH2-A, 516 nm for RH2-B and 492 nm for RH2-C, and between the two blue opsins, 439 nm for SWS2-A and 405 nm for SWS2-B. Zebrafish also has multiple opsin subtypes, but phylogenetic analysis revealed that medaka and zebrafish gained the subtype opsins independently. The lambda and BAC DNA clones isolated in this study could be useful for investigating the regulatory mechanisms and evolutionary diversity of fish opsin genes.     

A method of detecting conserved and variable segments in the amino acid sequences of homologous proteins

A set of aligned homologous protein sequences is divided into two groups consisting of m and n sequences. Each group contains sequences from the most related organisms. Value of the position dissimilarity of proteins from different groups of m and n sequences is defined as a number of mismatches in comparison of all possible m X n pairs of amino acid residues in the position (each from different group) divided by m X n. Ten position average of dissimilarity values is plotted vs. the first position number. Area of the figure between the profile of dissimilarity values and its mean value line characterizes the overall irregularity of amino acid substitutions along the protein sequences. If the area is greater than the average area for 1000 random profiles by more than two standard deviation units, the profile extrema containing the "surplus" of area are cut off. The cut-off stretches are likely to be variable and constant regions. If necessary, each of stretches may be separately tested and statistically estimated using a standard size sample of artificial protein families. Intergroup comparison of protein sequences reveals high overall irregularity of amino acid substitutions and identifies variable and conservative regions for all considered families of proteins: phospholipases A2, aspartate aminotransferases, alpha-subunits of Na+, K(+)-ATPase, L- and M-subunits of photosynthetic bacteria photoreaction centre, human rhodopsins.     

The molecular basis of adaptive evolution of squirrelfish rhodopsins

The wavelengths of maximal absorption (lambdamax) of the rhodopsins of nine squirrelfishes (N. sammara, N. argenteus, S. punctatissimum, S. microstoma, S. diadema, S. xantherythrum, S. spiniferum, N. aurolineatus, and S. tiere) and two soldierfishes (M. violacea and M. berndti) vary between 481 and 502 nm. Phylogenetic and mutagenesis analyses suggest that the common ancestor of these pigments had a lambdamax value of approximately 493 nm, and the contemporary lambdamax values were generated mostly by amino acid replacements E122M, F261Y, and A292S. The probability of observing all these amino acid replacements at specific branches of the phylogenetic tree is only 2.5 x 10(-9); it is highly unlikely that these changes have occurred by neutral evolution. Because of a close association between the lambdamax values of these pigments and the wavelengths of light available to the corresponding species, the excess number of amino acid changes at specific branches in the phylogenetic tree strongly suggests that the rhodopsins have undergone adaptive changes at various stages of the holocentrid evolution.