Skip to main content

A simple method for studying the molecular mechanisms of ultraviolet and violet reception in vertebrates



Many vertebrate species use ultraviolet (UV) reception for such basic behaviors as foraging and mating, but many others switched to violet reception and improved their visual resolution. The respective phenotypes are regulated by the short wavelength-sensitive (SWS1) pigments that absorb light maximally (λmax) at ~360 and 395–440 nm. Because of strong epistatic interactions, the biological significance of the extensive mutagenesis results on the molecular basis of spectral tuning in SWS1 pigments and the mechanisms of their phenotypic adaptations remains uncertain.


The magnitudes of the λmax-shifts caused by mutations in a present-day SWS1 pigment and by the corresponding forward mutations in its ancestral pigment are often dramatically different. To resolve these mutagenesis results, the A/B ratio, in which A and B are the areas formed by amino acids at sites 90, 113 and 118 and by those at sites 86, 90 and 118 and 295, respectively, becomes indispensable. Then, all critical mutations that generated the λmax of a SWS1 pigment can be identified by establishing that 1) the difference between the λmax of the ancestral pigment with these mutations and that of the present-day pigment is small (3 ~ 5 nm, depending on the entire λmax-shift) and 2) the difference between the corresponding A/B ratios is < 0.002.


Molecular adaptation has been studied mostly by using comparative sequence analyses. These statistical results provide biological hypotheses and need to be tested using experimental means. This is an opportune time to explore the currently available and new genetic systems and test these statistical hypotheses. Evaluating the λmaxs and A/B ratios of mutagenized present-day and their ancestral pigments, we now have a method to identify all critical mutations that are responsible for phenotypic adaptation of SWS1 pigments. The result also explains spectral tuning of the same pigments, a central unanswered question in phototransduction.


Early vertebrate ancestors used ultraviolet (UV) reception [1], which have been maintained in a wide range of modern species and play key roles in mate choice, foraging, predator avoidance, communication and migration [26]. In many other species, however, UV reception has been replaced by violet (or blue) reception during evolution [1]. Retinas of these animals do not receive UV light [7, 8] and are protected from UV damage and also their newly acquired violet reception have improved visual resolution and subtle contrast detection [9]. The avian ancestor also lost UV reception, but certain modern species regained it [1], which are again useful for various behaviors, including orientation based on the sun compass during migration [3, 6, 10].

UV and violet reception are regulated by the short wavelength-sensitive (SWS1) pigments that absorb light maximally (λmax) at ~360 and 395–440 nm, respectively [1]. To study the molecular basis of spectral tuning in a present-day pigment, vision scientists introduce mutations into orthologous as well as paralogous pigments and try to convert the λmaxs of these mutants to the λmax of interest [1113]. However, identical mutations in different pigments can generate different λmaxs and therefore variable sets of mutations can achieve similar λmaxs [14, 15]. Similarly, to recapitulate the evolution of a certain present-day pigment, molecular evolutionists infer the evolution in a reverse fashion by introducing mutations into the present-day pigment, but the magnitudes of λmax-shift caused by these mutations and the corresponding forward mutations in its ancestral pigment can differ drastically [1113, 16, 17]. Hence, using the traditional experimental approach, multiple sets of mutations can explain an observed λmax and the evolutionary mechanisms inferred can be misleading or even erroneous [14, 15, 18, 19]. This “multiple-solution” problem occurs because 1) pigment-specific non-additive (epistatic) interactions are ignored and 2) the current mutagenesis experiments are used to search for “any” mutations that can achieve targeted λmaxs. The first problem must be corrected by reconstructing appropriate ancestral pigments and manipulating them rather than present-day pigments [14, 15, 20, 21]. To address the second problem, we need to establish the one-to-one relationship between the SWS1 pigment and the dichotomous phenotypes; for that, a new approach is required.

Each visual pigment consists of an opsin and the retinal, either 11-cis-retinal or 11-cis-3, 4-dehydroretinal [22, 23]. At the chemical level, each visual pigment consists of a mixture of pigments with protonated Schiff base (SB) nitrogen-linked retinals (PSBR) and those with unprotonated SB nitrogen-linked retinals (SBR). When the SBR form is energetically more stable than PSBR, the SWS1 pigment is UV-sensitive; otherwise it is violet-sensitive [2426]. The protonation status of SB nitrogen is influenced strongly by the hydrogen-bond network (HBN) that is formed by amino acids at sites 86, 90, 113, 114, 118 and 295 and two water molecules [25]; throughout the paper, the amino acid site numbers follow those of bovine rhodopsin (GenBank accession no. M21606). Therefore, the HBN structure offers an opportunity to establish the one-to-one relationship between a SWS1 pigment and its dichotomous phenotype and identify all critical mutations unambiguously.


Ancestral pigments

Previously, we inferred the amino acid sequences of ancestral pigments by applying maximum likelihood-based Bayesian method (PAML [27]) to a phylogenetic tree of 21 SWS1 pigments and genetically engineered pigments of the ancestors of 1) jawed vertebrates (AncVertebrate), 2) Tetrapods (AncTetrapod), 3) Amphibians (AncAmphibian), 4) Amniotes (AncAmniote), 5) Sauropsids (AncSauropsid), 6) Aves (AncBird) and 7) Boreoeutheria (AncBoreotheria) [1]. Here we applied the PAML with JTT and WAG substitution models to a new phylogenetic tree of 33 representative SWS1 pigments (Fig. 1a) and inferred all ancestral pigments. The amino acid sequences inferred are highly consistent and have very similar posterior probabilities (PP). For AncVertebrate, for example, the two models predict the same amino acids at 244 (86 %) out of a total of 282 sites considered (see below) with PP ≥ 0.95 and at 23 of the remaining 38 sites with PP ≥ 0.70; different amino acids are predicted only at three sites (Additional file 1: Table S1).

Fig. 1

Construction of ancestral SWS1 pigments from 33 present-day sequences. (a) A composite phylogenetic tree with the λmaxs at 10 critical nodes. The numbers beside various branches indicate the numbers of mutations introduced into present-day and ancestral pigments. Arrows indicate the direction of the reconstruction of ancestral pigments. The numbers after each organism (and their SWS1 pigments) indicate their λmaxs. UV- and violet-sensitive pigments are distinguished by black and blue rectangles, respectively. (b) The neutral and adaptive mutations in ancestral pigments (black and blue circles, respectively). (c) Statistically inferred ancestral amino acids that have PP < 95 % (in black circles). The 2-D models are after Palczewski [97]

Since JTT and WAG models predict very similar ancestral sequences, we used the amino acid sequences predicted by JTT model (Additional file 2: Figure S1) and introduced the necessary amino acid changes into the internal segment (sites 31–312) in the pMT5 containing the N- and C-termini of the chameleon-359 (for the justification, see Methods) and engineered AncVertebrate, AncTetrapod, AncAmphibian, AncAmniote as well as those of Euteleosts (AncEuteleost), Mammals (AncMammal) and Eutherians (AncEutheria) (Additional file 2: Figure S1; Fig. 1a). The in vitro assays [28] show that the λmaxs of AncVertebrate, AncTetrapod, AncAmphibian and AncAmniote which are based on the two different phylogenetic trees have λmaxs of ~ 360 nm and, furthermore, with the exception of AncBird-393, the ancestral pigments (AncVertebrate-361, AncEuteleost-364, AncTetrapod-359, AncAmphibian-359, AncAmniote-359, AncSauropsid-360, AncMammal-359, AncEutheria-360 and AncBoreotheria-360 [29]) are all UV-sensitive, where the numbers after pigment names indicate their λmaxs (Fig. 1a, Additional file 3: Figure S2).

When the amino acids at a total of 282 sites of the 10 ancestral SWS1 pigment sequences are compared, we can find 85 polymorphic sites (Fig. 1b, Additional file 2: Figure S1), among which only amino acid changes F49V, F86S, L116V and S118A (or F49V/F86S/L116V/118S) cause the significant λmax-shift between AncSauropsid-360 and AncBird-393 [1]. Hence, amino acid changes at 81 out of the 85 sites (~95 %) do not cause any significant λmax-shift and may be considered as “selectively neutral” changes (Fig. 1b). Virtually identical proportions of “neutral” changes have been observed for paralogous RH1 (or rhodopsin) and middle and long wavelength-sensitive (M/LWS) pigments [30, 31]. It is also comforting to learn that, despite having a total of 70 amino acid sites with PP < 0.95, none of the λmaxs of the 10 ancestral pigments are affected significantly by the uncertainty involved in the statistical inference (Fig. 1c).

Mutagenesis-induced phenotypes; a small number of mutations

To better understand the nature of the currently available mutagenesis results of SWS1 pigments [1315], we added our own 23 new data (Methods). Using these old and new mutagenesis results, we first establish that the λmax-shifts caused by a small number of mutations often deviate drastically from the additive mutational effects and depend strongly on their molecular backgrounds (Fig. 2, Additional file 4: Table S2). Probably the most unusual mutation is the deletion of the codon 86, which encoded F86 in AncVertebrate-361, in scabbardfish-423 of Lepidopus fitchi. When we consider present-day pigments, the violet-sensitivity of scabbardfish-423 and the UV-sensitivities of bfin killifish-355 (Lucania goodie) and lampfish-371 (Stenobrachius leucepsarus) can be interchanged largely by deleting F86 from the UV pigments and inserting F86 into scabbardfish-423, respectively, which suggest that F86 deletion has significantly contributed to scabbardfish-423 evolution. However, the λmax-shifts caused by the two deletion mutants differ by ~20 nm and, furthermore, when F86 is deleted from AncVertebrate-361, the λmax of the mutant is 43 nm shorter than that of scabbardfish-423 [26]. Much to our surprise, when F86 is deleted from AncEuteleost-364, which is more closely related to scabbardfish-423 than AncVertebrate-361 (Fig. 1a), the mutant pigment becomes structurally unstable and its λmax cannot be evaluated. All of these results show that additional mutations are involved in the scabbardfish-423 evolution.

Fig. 2

The λmax-shifts generated by various mutations. The lengths of arrows represent the λmax-shifts and filled circles indicate that no λmax-shift occurred. Red and black arrows show the mutational effects of ancestral and present-day UV pigments, respectively, whereas blue arrows indicate those of violet pigments. The different λmax-shifts of bovine mutants with Y86F have been evaluated by using dark spectra [84] and dark–light spectra [92] and only the former result is shown. AncBird*-393 is identical to AncSauropsid-360 with mutations F49V/F86S/L116V/S118A

One major characteristic of the mutagenesis results is that the magnitude of a λmax-shift (or Δλmax) caused by mutations tends to be much smaller in UV pigments than in violet pigments. We can see this in the 11 sets of comparisons (Fig. 2): 1) F86 deletion in AncVertebrate-361 and F86 insertion in scabbardfish-423 (Δλmax = 19 vs −60 nm, respectively); 2) F86L in AncBoreotheria-360 and L86F in human-414 (0 vs −32 nm); 3) F86M in AncAmphibian-359 and M86F in frog-423 (Xenopus laevis) (2 vs −25 nm); 4) F86S in AncSauropsid-360 and AncEutheria-360 and S86F in AncBird-393 and elephant-419 (Loxodonta africana) (17 vs −38 and 14 vs −52 nm, respectively); 5) F86Y in AncBoreotheria-360 and Y86F in bovine-438 (Bos taurus) and squirrel-440 (Sciurus carolinensis) (48 vs −71 ~ −76 nm); 6) C86S in zebra finch-359 (Taeniopygia guttata) and S86C in the functionally equivalent violet pigments (AncBird-393 and AncSauropsid-360 with F49V/F86S/L116V/118S, or AncBird*-393) (1 vs −7 or −27 nm); 7) T93I in AncBoreotheria-360 and I93T in bovine-438 (−4 vs −22 nm); 8) T93I/S97T in AncBoreotheria-360 and I93T/T97S in bovine-438 (1 vs −39 nm); 9) T93I/L116V in AncEutheria-360 and I93T/V116L in elephant-419 (0 vs −8 nm); 10) F86S/L116V in AncEutheria-360 and S86F/V116L in elephant-419 (38 vs −59 nm) and 11) F86Y/T93I/S97T in AncBoreotheria-360 and the reverse changes in bovine-438 (48 vs −79 nm).

Hence, epistatic interactions tend to operate more strongly in UV pigments than in violet pigments. Furthermore, identical mutations cause variable λmax-shifts among orthologous pigments. S90C contributed significantly to the evolution of UV-sensitivities of some modern avian pigments (e.g. zebra finch-359, canary-366 (Serinus canaria) and budgerigar-360 (Melopsittacus undulatus)) from the violet-sensitive AncBird-393 [1, 32, 33] (Fig. 1a). S90C in AncBird*-393, pigeon-393, chicken-415, frog-423 and bovine-438 decreases their λmaxs by 7–46 nm, whereas the identical mutation causes no λmax-shift in mouse-359 (Additional file 4: Table S2). Variable Δλmaxs with the same mutations can also be found in 1) deletion of F86 from AncVertebrate-361, lampfish-371 and bfin killifish-354 (Δλmax = 19–58 nm), 2) S86F in AncBird-393 and elephant-419 (Δλmax = −38 and −52 nm, respectively), 3) F86Y in AncBoreotheria-360, mouse-359 and goldfish-360 (Δλmax = 48–66 nm), 4) S86C in phenotypically identical AncBird-393 and AncBird*-393 (Δλmax = −7 and −27 nm, respectively) and 5) I93T in AncEutheria-360, AncBoreotheria-360, mouse-359, elephant-419 and bovine-438 (Δλmax = 0–22 nm).

In summary, 1) identical mutations can cause, sometimes drastically, different λmax-shifts, 2) forward and reverse mutations can shift the λmax by different magnitudes to the opposite directions, or even in the same direction, and 3) UV pigments, particularly ancestral UV pigments, tend to be less responsive to mutations than violet pigments to the corresponding reverse changes. Two sets of forward and reverse mutations shift the λmax in the same direction: 1) T93I in AncBoreotheria-360 and I93T in elephant-419 and bovine-438 and 2) E113D in AncAmphibian-359 and D113E in frog-423 (Additional file 4: Table S2). The differential effects of forward and reverse mutations clearly show that the evolutionary mechanisms of UV and violet reception must be studied by using ancestral pigments rather than present-day pigments. One notable exception is Y86F in wallaby-420 (Macropus eugenii) and F86Y in AncMammal-359, which fully interchange the two original λmaxs (Fig. 2; Additional file 4: Table S2).

At the chemical level, each SWS1 pigment consists of a mixture of PSBR and SBR (see Background). The major λmax-shifts of SWS1 pigments are caused by changes in the relative ground-state energies of the pigments with the two retinal groups. The calculated relative ground-state energies of a SWS1 pigment with SBR subtracted from that with PSBR (ΔE) is positive (varying between 2.5 and 8.3 kcal/mol) for a UV pigment while it is negative for a violet pigment (varying between −2.4 and −5.6 kcal/mol) [25]. The wider ΔE range explains the functionally conservative nature of UV pigments.

Multiple mutations

As the number of critical mutations identified increases, the magnitudes of λmax-shifts caused by forward and reverse mutations tend to become similar. Since epistatic interactions are reflected better by multiple mutations than by single mutations, this observation may be expected. This trend can be seen in 1) F86S/T93I in AncEutheria-360 and S86F/I93T in elephant-419 (Δλmax = 55 vs −60 nm, respectively), 2) F86Y/T93I in mouse-359 and Y86F/I93T in bovine-438 (66 vs −68 nm) and 3) F86S/T93I/L116V in AncEutheria-360 and the reverse mutations in elephant-419 (50 vs −59 nm) (Fig. 2, Additional file 4: Table S2). We can find three examples of excellent symmetry between the λmax-shifts caused by forward mutations in an ancestral pigment and reverse mutations in a corresponding present-day pigment: 1) F49V/F86S/L116V/S118A in AncSauropsid-360 and the reverse mutations in AncBird-393 (33 vs −33 nm); 2) F86M/V91I/T93P/V109A/E113D/L116V/S118T in AncAmphibian-359 and the reverse mutations in frog-423 (62 vs −68 nm) and 3) F46T/F49L/T52F/F86L/T93P/A114G/S118T in AncBoreotheria-360 and the reverse mutations in human-414 (52 vs −54 nm) (Fig. 2).

The goal of all of these mutagenesis analyses is to find the molecular mechanisms of spectral tuning and evolution of a present-day pigment. A weakness of this traditional approach becomes apparent from the mutagenesis analyses of elephant-419 evolution. F86S/T93I in AncEutheria-360 and S86F/I93T in elephant-419 achieve Δλmaxs of 55 and −60 nm, respectively (Additional file 4: Table S2), which interchange the λmaxs of the two pigments reasonably well and elephant-419 seems to have evolved from AncEutheria-360 by F86S/T93I. However, elephant-419 has incorporated 14 additional mutations and AncEutheria-360 with F86S/T93I/L116V attains a Δλmax of 50 nm (Additional file 4: Table S2), which moves further away from the λmax of elephant-419, which show that neither F86S/T93I nor F86S/T93I/L116V explain elephant-419 evolution. Hence, to identify all critical mutations, it is necessary, but not sufficient, to manipulate and compare the λmaxs of present-day pigments and their ancestral pigments. To alleviate this type of problem, we may check whether mutations that attained the desired λmax-shift also achieve the key protein structural change.

Molecular modelling of Hydrogen-Bond Network (HBN): AMBER models

We divided the HBN region into two parts: one area formed by amino acids at sites 90, 113 and 118 (area A) and another area determined by those at sites 86, 90 and 118 and 295 (area B), each of which is further subdivided into three parts A1-A3 and B1-B3 (Fig. 3a). Since the amino acids at sites 113 and 114 are bonded to each other, their effects on the HBN structure and λmax-shift are correlated [25]. Therefore, one of these sites, i.e., 113, is included in area A. The distance between SB nitrogen that lies in the area A (Fig 3a) and the counterion E113 in AncVertebrate-361 is 1.99 Å, but the corresponding distances in AncBird-393 and human-414 are ~1.80 Å [25]. The wider distances contribute in making the SBN forms of UV pigments energetically more stable than PSBR forms, or vice versa for violet pigments. As the distance between E113 and the SB nitrogen decreases, the A/B ratio tends to decrease and the λmax of a pigment is expected to increase (Fig. 3a).

Fig. 3

The HBN regions of SWS1 pigments. (a) Area A (A1-A3) contains SB nitrogen and the counter-ion E113 and area B (B1-B3) includes two water molecules. The example shows the six amino acids of AncVertebrate-361 and the 11-cis-retinal with SB nitrogen (N). Dotted lines indicate the shortest distances between two amino acids. (b) The λmaxs and A/B ratios of the ancestral and present-day pigments. UV pigments are shown by black circles (group 1), whereas the violet-sensitive AncBird-393 and pigeon-393 (group 2) and the others (group 3) are shown by blue circles. The widths of light blue rectangles show the 95 % confidence intervals of the A/B ratios for group 1 (0.54–0.55) and group 3 (0.43–0.51), whereas the corresponding intervals for group 2 is a narrow strip at around 0.49 (not shown) (Additional file 5: Table S3). Scabbardfish-423 lacks F86 and the most closely located amino acid to the HBN region is V91 and therefore its A/B ratio is evaluated by considering V91 in place of F86

For five representative pigments (AncVertebrate-359, AncBird-393, zebra finch-359, budgerigar-363 and human-414), initial structures were obtained from their amino acid sequences through SWISS MODEL homology modelling [34] using the structure of bovine rhodopsin (pdb code: 1U19) as a template. After adding hydrogen atoms and fully optimizing the modelled structures at pure AMBER96 force field level, we evaluated the sizes of areas A and B as well as A/B ratios.

The results show that UV pigments (AncVertebrate-359, zebra finch-359 and budgerigar-363) have A/B ratios of 0.507–0.527, followed by 0.421 of AncBird-393 and 0.406 of human-414, in that order, which show an inverse relationship between the A/B ratio and λmax of a visual pigment (Table 1). To relate such structural variations with the λmax values of the visual pigments, we performed multivariable regression analysis between the λmax values (dependent variable) and structural parameters like A/B and B (independent variables) and obtained a relationship, λmax = 390.33–357.284 x (A/B) + 2.345 x B. The λmaxs predicted from this formula are very close to the actual values determined by the in-vitro assay (R 2 = 0.969) and standard error (S = 6.065) and the mean absolute error (MAE = 3.085) are satisfactorily small (Table 1). Hence, the λmax of a pigment decreases as the A/B ratio increases and B value decreases.

Table 1 Comparison of A/B values of HBNs using SWISS and AMBER models

Simpler HBN models: SWISS models

The A/B ratios were also obtained directly from SWISS MODEL method considering only the amino acid sequences. We refer to the HBN structures before and after AMBER geometry optimizations as SWISS models and AMBER models, respectively (Table 1). The A/B ratios of SWISS models are larger than those of AMBER models and the overestimation is largest for human-414. Because of this, the inverse relationship between A/B ratio and λmax does not hold for the less accurate SWISS models of the five pigments. However, when we consider the SWISS models of the 10 ancestral and all currently known 33 present-day SWS1 pigments, they can be distinguished roughly into three groups: 1) the A/B ratios of the SWISS models of the UV pigments with λmaxs of 354–371 nm (0.547 ± 0.003, group 1) are larger than those of AncBird-393 and pigeon-393 (0.492 ± 0.0003, group 2), which tend to be larger than the A/B ratios of violet pigments (0.469 ± 0.021, group 3) (Fig. 3b, Additional file 5: Table S3). Like those of AMBER models, the smallest A/B ratios of the group 3 (or violet) pigments are caused by the compressed A region plus the expanded B region and the intermediate A/B ratios of the SWISS models of group 2 pigments come from an expanded B region (Additional file 5: Table S3).

Human-414, Squirrel-440, bovine-438 and wallaby-420 have much larger A/B ratios than the rest of the group 3 pigments; similarly, zebra finch-359 and bfin killifish-355 have much larger A/B ratios than the other group 1 pigments (Fig. 3b, Additional file 5: Table S3). During the evolution of human-414 from AncBoreotheria-360, three critical changes (F86L, A114G and S118T) have been incorporated in the HBN region. These changes make the compression of A region and expansion of B region in human-414 less effective in the SWISS models than in AMBER models and generate the higher A/B ratio of its SWISS model (Table 1). For the same reason, F86Y in squirrel-440, bovine-438 and wallaby-420 as well as F86C and S90C in zebra finch-359 and S118A in bfin killifish-355 have generated the large A/B ratios of their SWISS models. The smallest A/B ratio of scabbardfish-423 comes from its unique protein structure, in which V91 needs to be considered in place of F86.

The major advantage of using the less accurate SWISS models is that they are readily accessible to everyone and, importantly, the A/B ratios of the SWISS models of UV pigments (0.535–0.577) can still be distinguished from those of violet pigments (0.355–0.526) (Fig. 3b). In analysing SWS1 pigments, the variable λmaxs and A/B values within each of the three pigment groups are irrelevant because we are concerned mainly with the major λmax-shifts among UV pigments (group 1), AncBird-393 (group 2) and violet pigments (group 3): group 1 → group 2, group 1 → group 3, group 2 → group 1 and group 2 → group 3 (Fig. 1a). For each of these phenotypic adaptive processes [1], we can establish the one-to-one relationship between A/B ratios and dichotomous phenotypes of SWS1 pigments.

Criteria for acceptable mutagenesis results

To examine whether or not the mutagenesis result of a certain present-day pigment reflects the epistatic interactions correctly, we evaluate the λmax and A/B ratio of its ancestral pigment subtracted from those of a mutant pigment (denoted as d(λmax) and d(A/B), respectively). Similarly, the validity of the mutagenesis result of an ancestral pigment can be examined by evaluating its d(λmax) and d(A/B) values by considering the λmax and A/B ratio of the corresponding present-day pigments. Following the traditional interpretation of mutagenesis results, it seems reasonable to consider that present-day and ancestral mutant pigments fully explain the λmaxs of the target (ancestral and present-day) pigments when d(λmax) ≤ 3 ~ 5 nm, depending on the magnitudes of total λmax-shift considered. Following the mutagenesis results of wallaby-420, AncBird-393, frog-423 and human-414 (see below), the A/B ratio of the target pigment may be considered to be fully converted when d(A/B) ≤ 0.002.

Searching for the critical mutations in SWS1 pigments

Considering d(λmax) and d(A/B) together, mutagenesis results of SWS1 pigments can be distinguished into three classes: 1) amino acid changes satisfy d(λmax) ≤ 3 ~ 5 nm and d(A/B) ≤ 0.002 (class I); 2) those satisfy only d(λmax) ≤ 3 ~ 5 nm (class II) and 3) those satisfy neither conditions (class III) (Table 2, Additional file 5: Table S3). When only a small number of mutations are considered, class I includes F86 insertion in scabbardfish-423 and Y86F in wallaby-420, both achieving d(λmax) = 2 nm and d(A/B) = 0.001. However, the F86 deletion mutants of AncVertebrate-361, lampfish-371 and bfin killifish-355 all belong to class III, confirming that scabbard-423 did not evolve by F86 deletion alone. On the other hand, F86Y in AncMammal-359 belongs to class I, establishing that wallaby-420 indeed evolved from AncMammal-359 by F86Y alone. Compared with these two examples, Y86F in squirrel-440 and F86Y in AncBoreotheria-360 belong to classes II and III, respectively, showing that squirrel-440 evolution did not occur by F86Y alone.

Table 2 Comparisons of d(λmax) and d(A/B) for different sets of pigments

Class I also includes three sets of reverse mutations: 1) V49F/S86F/V116L/A118S in AncBird-393, 2) M86F/I91V/P93T/A109V/D113E/V116L/T118S in frog-423 and 3) T46F/L49F/F52T/L86F/P93T/G114A/T118S in human-414. The corresponding forward mutations in AncSauropsid-360, AncAmphibian-359 and AncBoreotheria-360 also belong to class I (Table 2). Hence, AncBird-393 evolved from AncSauropsid-360 by four mutations, while frog-423 and human-414 evolved from their ancestral pigments by a different set of seven mutations. On the other hand, despite their significant magnitudes of λmax-shifts, individual mutations L86F in human-414 (Δλmax = −32 nm and d(λmax) = 22 nm) and M86F in frog-423 (Δλmax = −25 nm and d(λmax) = 39 nm) belong to class III (Additional file 4: Table S2). Moreover, Y86F in bovine-438 decreases the λmax by ~70 nm, but this mutation (d(λmax) = 7 nm) still belongs to class III and furthermore class III status of F86Y in AncBoreotheria-360 shows that the evolutionary mechanism of bovine-438 is still unsolved (Table 2).

Among the three classes, class II is particularly disconcerting because even when the λmaxs of present-day pigments can be converted to those of their ancestral pigments, these mutations do not achieve the key protein structural changes. Class II includes Y86F of squirrel-440 as well as S86F/I93T and S86F/I93T/V116L of elephant-419 (Table 2). Hence, either additional mutations can be involved or they might not have played significant roles during evolution (see Discussion). As suspected, class III includes many single mutations, which are represented by such mutations as L86F in human-414, M86F in frog-423, Y86F in bovine-438 and S86F in elephant-419.

In summary, the goal of studying molecular basis of spectral tuning in a present-day pigment is to identify mutations that generated its λmax, while the mechanism of phenotypic adaptation of the same pigment is to find specific mutations that generated the λmax during evolution. These questions address the same phenomenon and can be solved simultaneously; for the latter problem, however, it would also be necessary to establish the relationship between the phenotypic changes and the changes in the organisms’ new environments (see the next section). Hence, among all mechanisms of spectral tuning and adaptive evolution of SWS1 pigments proposed to date, only those for AncBird-393, frog-423, human-414 and wallaby-420 can be supported.


Mutations in different molecular backgrounds can differ significantly in their contribution to phenotypic adaptation [14, 15]. Here we have seen that 1) mutagenesis results of present-day SWS1 pigments are highly pigment-specific and 2) the one-to-one relationship holds between A/B ratios of HBN region and dichotomous phenotypes (UV- and violet-sensitivities) of SWS1 pigments. We then developed a method for identifying all critical mutations that generated the λmaxs of present-day pigments by interchanging the λmaxs and A/B ratios of the present-day and their ancestral pigments. In applying this approach, the following three features are critical to keep in mind.

First, it is imperative to reconstruct proper ancestral pigments and manipulate them. To see the necessity of manipulating ancestral molecules and their phenotypes, we revisit the evolution of elephant-419 from AncEutheria-360. At present, this process is explained best by F86S/T93I/L116V, which are responsible for about 85 % of the entire λmax-shift and the triple mutant in AncEutheria-360 achieve neither the λmax nor A/B ratio of elephant-419 (class III), but the reverse mutant of elelphant-419 attains the λmax of AncEutheria-360, but not the A/B ratio (class II mutations) (Table 2). Introducing all possible combinations of the three mutations into elephant-419 and applying a linear model to their λmaxs and that of elephant-419 (λelephant), the individual and epistatic effects of these mutations on the λmax-shift (θ) were evaluated (Table 3). The results show that the major contributor is S86F (θS86F = −52 nm) and the effects of forward mutations in AncEutheria-360 may be inferred by reversing the sign of this θ value. Alternatively, S86F, I93T, V116L, S86F/I93T, S86F/V116L, T93I/V116L and S86F/I93T/V116L in elephant-419 may be regarded as functionally equivalent to T93I/L116V, F86S/L116V, F86S/T93I, L116V, T93I, F86S in AncEutheria-360 and the ancestral pigment, respectively. Then again, θF86S (51 nm) has a major impact in the elephant-419 evolution.

Table 3 Effects of mutations on the λmax-shift

Consequently, elephant-419 seems to have evolved mostly by F86S. However, this conclusion is incorrect. That is, when we introduce the corresponding forward mutations into AncEutheria-360, epistatic interactions (θF86SxT93I = 41 nm, θF86SxL113V = 25 nm and θF86SxT93IxL113V = −30 nm) have major impacts and the F86S-effect (14 nm) becomes less significant (Table 3), again showing that epistatic interactions are much stronger in the UV pigment than in elephant-419. This example demonstrates that the spectral tuning and evolutionary mechanism of a present-day pigment must be studied by manipulating its ancestral pigment.

Second, the A/B ratio can be useful for checking whether certain mutations that cause significant λmax-shifts were actually used for phenotypic (or functional) changes. For example, S86C and S90C in AncBird*-393 decrease the λmax to 366 and 360, respectively. In theory, both mutations explain the reversion from violet reception to UV reception in certain modern avian species. When the respective mutants are compared to budgerigar-359, d(A/B) values are 0.485–0.535 (= −0.050) and 0.528–0.535 (= −0.007); similarly, when they are compared to zebra finch-359, the d(A/B) values are 0.485–0.577 (= −0.092) and 0.528–0.577 (= −0.049), respectively (Additional file 4: Table S2 and Additional file 5: Table S3). In reality, therefore, the S86C and S90C mutants belong to classes III and II, respectively; moreover, the smaller d(λmax) and d(A/B) values suggest that S90C, not S86C, has contributed to the actual evolution of the avian UV pigments. Indeed, phylogenetic analyses strongly suggest that S90C preceded S86C and the effect of the latter mutation seems to have been insignificant during evolution [1].

Third, as suggested by S86C and S90C, d(λmax) and d(A/B) of a pigment are affected strongly by the order of mutation accumulations. For example, the seven critical mutations in AncAmphibian-359 and those in AncBoreotheria-360 shift the λmax individually only slightly, if any [19, 35], but as they start to accumulate, their contributions to the λmax-shift become much more significant (Additional file 5: Table S3). More specifically, when the transmembrane III of AncAmphibian-359 is replaced by that of frog-423, the mutant pigment increases the λmax to 410 nm (Δλmax = 51 nm). This change is caused by V109A, E113D, L116V and S118T, but E113D in AncAmphibian-359 decreases its λmax by 4 nm and V109A/L116V/S118T cause no λmax-shift at all [19]. This is not the whole story. Phylogenetic analyses strongly suggest that E113D occurred toward the end of frog-423 evolution; in that case, the four mutations contributed only 15 nm rather than 51 nm [35].

In addition to the SWS1 pigments, epistatic interactions have been found in RH1, RH1-like (RH2), SWS type 2 (SWS2) and M/LWS pigments as well [30, 3639]. The most extensively studied mutations are forward (A292S) and reverse (S292A) changes in various pigments [17]. The λmax-shifts caused by A292S in RH1, SWS2 and M/LWS pigments range between −10 and 1 nm (e.g. [30, 40]), between −8 and −4 nm (e.g. [37, 38]) and −33 nm [39], respectively, whereas S292A in RH1, RH2, SWS1 and M/LWS pigments increase the λmax by 4–12 nm [30], 7 nm [36], 0 nm [41] and 18–28 nm [40, 42], respectively. Hence, A292S and S292A mutants shift the λmax by −33 ~ 1 and 0 ~ 28 nm, respectively, reflecting the varying levels of epistatic interactions between site 292 and its molecular background.

In studying possible molecular adaptation, it is common to “claim evidence of adaptive evolution based on computational analyses alone [43].” These sequence analyses are based basically on the assumption that “to detect positive Darwinian selection, it is necessary to show that the number of nonsynonymous substitutions per nonsynonymous site (dn) is significantly greater than that of synonymous substitutions per synonymous site (ds) [44, 45].” Indeed, such statistical methods have been used widely particularly in comparative genomics ([46, 47]; for SWS1 genes see [48]). However, given high rates of false-positives and false-negatives [30, 49], their reliabilities are questionable [47, 5052]. Certainly, adaptive mutations identified using statistical methods generate biological hypotheses [52, 53], which must be tested eventually using experimental means [46, 47, 54].

We have seen that about 95 % of amino acid changes can be considered as “selectively neutral” changes. Currently, mutations at a total of 17 sites (positions 46, 49, 52, 86, 90, 91, 93, 97, 109, 113, 114, 116, 118, 179, 207, 256 and 277) can shift the λmaxs of various SWS1 pigments, which seem to have contributed for adaptive evolution of various vertebrate species [13, 17, 35]. The average rates of nucleotide substitution at all positions of 300 codons of violet and UV pigments are both ~0.70 x 10−9/site/year; hence, as suspected, the effects of the small proportion of adaptive sites on the evolutionary rates are buried among those of the neutral changes (Table 4). However, when we consider the 16 codon sites that are involved in the λmax-shift, the rate of nucleotide substitution is significantly higher for violet pigments (1.73 x 10−9/site/year) than for the functionally unchanged UV pigments (0.70 x 10−9/site/year) (Table 4). This suggests that adaptive sites tend to be more susceptible for mutation accumulation [35].

Table 4 Rates of nucleotide substitutions for pairs of violet and UV opsin genes that are evaluated by comparing to those of distantly related third opsin genes

The molecular analyses of phenotypic changes can also be found in such vertebrate systems as digestive enzymes secreted by the pancreas [55, 56], haemoglobins [5760], steroid receptors [6164] and olfactory receptors [6567]. If we are serious about testing the statistical hypotheses of molecular adaptation, then this is an opportune time not only to explore these and other new genetic systems but also to improve the procedures to study the molecular mechanisms of phenotypic adaptation. Since phenotypic changes do not necessarily mean that they are adaptive, their adaptive nature must be established by linking phenotypic differences to differences in animals’ ecological and physiological environments [20, 21, 68]. Characterizing visual pigments considering HBN regions and establishing the one-to-one relationship between SWS1 pigment and UV/violet reception is a significant development in that direction.


At present, molecular adaptations in vertebrates are studied almost always using comparative sequence analyses. These statistical results, however, provide only biological hypotheses and must be tested using experimental means. UV and violet reception, mediated by the SWS1 pigments with λmaxs at ~360 and 395–440 nm, respectively, offer such an opportunity. The mechanisms of spectral tuning (or λmax-shift) and phenotypic adaptation of a present-day SWS1 pigment can be elucidated not only by engineering its ancestral pigment but also by using a new signature of protein structure, the A/B ratio. We can identify critical mutations that generated a present-day pigment by matching the λmax and A/B ratio of an ancestral pigment with candidate mutations to that of the present-day pigment descended from it. Establishing the one-to-one relationship between SWS1 pigment and UV/violet reception is a significant development in actually testing statistical hypotheses of positive selection and also in understanding the mechanism of spectral tuning of the same pigment.


Reconstruction of ancestral SWS1 pigments

Previously, applying PAML [27] to a composite phylogenetic tree of 21 SWS1 pigments, the ancestral SWS1 pigments of jawed vertebrates (AncVertebrate), tetrapods (AncTetrapod), amphibians (AncAmphibian), amniotes (AncAmniote), Sauropsids (AncSauropsid), birds (AncBird) and Boreoeutheria (AncBoreotheria) were inferred and reconstructed [1].

Here, we considered the composite evolutionary tree of 33 representative SWS1 pigments based on molecular (e.g. and paleontological data (Fig. 1a): lamprey-358 (Lamptera marinus, U67123; λmax = 358 [69]), goldfish-359 (Carassius auratus auratus, D85863; [70]), zebrafish-355 (Danio rerio, AB087810; [71]), scabbardfish-423 (Lepidopus fitchi, FJ443126; [26]), tilapia-360 (Oreochromis niloticus, AF191221; [72]), cichlid-Mzeb-368 (Maylandia zebra, AF191219; [73]), bfin killifish-355 (Lucania goodie, AY296735; [37]), medaka-356 (Oryzias latipes, AB223058; [74]), Pacific saury (Cololabis saira; KP099197), toothfish-369 (Dissostichus mawsoni, AY927651; [75]), lampfish-371 (Stenobrachius leucepsarus, FJ443127; [26]), frog-423 (Xenopus laevis, U23463; [19]), salamander-356 (Ambystoma tigrinum, AF038948; [76]), zebra finch-359 (Taeniopygia guttata, AF222331; [32]), canary-366 (Serinus canaria, AJ277922; [77]), budgerigar-363 (Melopsittacus undulates, Y11787; [33]), pigeon-393 (Columba livia, AF149234; [70]), chicken (Gallus gallus, M92039; [32]), chameleon-359 (Anolis carolinensis, AF134192; [78]), gecko-364 (Gekko gecko, AY024356; [79]), sunbeam-361 (Xenopeltis unicolor, FJ497234; [80]), human-414 (Homo sapiens, M13295; [18]), macaque-415 (Macaca fascicularis, AF158977; [81]), squirrel monkey-430 (Saimiri sciureus, U53875; [82]), marmoset-423 (Callithrix jacchus, L76201; [83]), bovine-438 (Bos taurus, U92557; [84]), mouse-359 (Mus musculus, U49720; [70]), rat-358 (Rattus norvegicus, U63972; [70]), squirrel-440 (Sciurus carolinensis, DQ302163; [85]), guinea pig-420 (Cavia porcellus, AY552608; [86]), elephant-419 (Loxodonta Africana, AY686753; [87]), wallaby-420 (Macropus eugenii, AY286017; [88]) and dunnart-363 (Sminthopsis crassicaudata, AY442173; [89]).

Using these sequences and those of RH1 pigment of bovine (Bos taurus, M21606), RH2 pigment of goldfish (Carassius auratus, L11865) and SWS2 pigment of zebrafish (Danio rerio, AB087809) as the outgroup, we inferred the amino acid sequences of SWS1 pigments at various nodes of the phylogenetic tree using PAML [27] using JTT and WAG substitution models.

When the amino acid sequences in the N and C termini of SWS1 pigments isolated from a wide range of vertebrate species are compared, the structures of the two segments differ significantly by including insertions and deletions. To evaluate the effects of these variable N and C termini on the λmax-shift, the N terminus (amino acids between sites 1 and 30) and the C terminus (those between sites 313 and 348) of the SWS1 pigment of mouse-359 were replaced by those of goldfish (Carassius auratus), chameleon (Anolis carolinensis) and human. The respective chimeric pigments had λmax values of 360, 359 and 360 nm, showing that the highly polymorphic amino acids at the two termini do not modify the λmax values of SWS1 pigments [1]. Hence, the previously reconstructed ancestral pigments contained the identical N and C termini of the chameleon SWS1 pigment plus the ancestral pigment-specific internal segment between sites 31 and 312 [1]. Following the same strategy, we reconstructed several ancestral pigments by introducing the necessary amino acid changes into the internal segment (sites 31 and 312) in the pMT5 containing the N- and C-termini of the chameleon-359.

All mutant opsins were generated by using QuickChange site-directed mutagenesis kits (Stratagene, La Jolla, CA). To rule out spurious mutations, the DNA fragment was sequenced by cycle sequencing reactions using the Sequitherm Excel II long-read kits (Epicentre Technologies, Madison, WI) with dye-labeled M13 forward and reverse primers. Reactions were run on a LI-COR (Lincoln, NE) 4300LD automated DNA sequencer.

Mutagenesis genesis data of ancestral and present-day pigments

Previously, the λmaxs of ancestral and present-day pigments as well as mutant pigments (mutations in subscripts) have been obtained by various authors: AncSauropsid-360, AncBird-393, AncBird*-393, AncSauropsid-360F86S, AncBird*-393S86C, AncBird*-393S90C, AncSauropsid-393A114G, AncSauropsid-360F49V/F86S/L116V/S118A, AncBird-393V49F/S86F/V116L/A118S [1], AncBoreotheria-360, AncBoreotheria-360F86L [35], bfin killifish-354 [37], goldfish-360, Pigeon-393, mouse-359 [70], lampfish-371, scabbardfish-423, AncVertebrate-361del (F86), lampfish-371del (F86), bfin killifish-355del (F86), scabbardfish-423ins (F86) [26], chameleon-359 [78], zebra finch-359, zebra finch-359C90S, chicken-415S90C, pigeon-393S90C [32, 90], budgerigar-363, budgerigar-363C90S [33], dunnart-363 [89], frog-423, AncAmphibian-359F86M, AncAmphibian-359E113D, AncAmphibian-359F86M/V91I/T93P/V109A/E113D/L116V/S118T, frog-423M86F/I91V/P93T/A109V/D113E/V116L/T118S [19], chicken-413 [91], human-414, mouse -359S90C, mouse-359A114G, human-414T46F/L49F/F52T/L86F/P93T/G114A/T118S [18], bovine-438, mouse -359F86Y, bovine-438Y86F, bovine-438S90C, mouse-359T93I, mouse-359F86Y/T93I, bovine-438Y86F/I93T, bovine-438Y86F/I93T/T97S [84], squirrel-440, squirrel-440Y86F [85], elephant-419, elephant-419S86F, elephant-419I93T, elephant-419S86F/I93T, elephant-419I93T/V116L, elephant-419S86F/V116L [87], wallaby-420 [88], bovine-438Y86F, goldfish-360F86Y [92], frog-423S90C [93] and frog-423D113E [94]. Note that two different Δλmax values have been obtained from a dark spectrum [84] or from a dark–light spectrum [92]. In our comparison, we considered the former results simply because it is a more direct measurement.

We also added the λmaxs of AncVertebrate-361, AncEuteleost-364, AncTetrapod-359, AncAmniote-359, AncAmphibian-359, AncMammal-359, AncEutheria-360, AncEuteleost-364del (F86), human-414L86F, frog-423M86F, AncEutheria-360F86S, AncBird-393S86F, AncMammal-359F86Y, wallaby-420Y86F, AncBoreotheria-360F86Y, AncBird-393S86C, zebra finch-359C86S, AncEutheria-360T93I, AncEutheria-360T93I, bovine-438I93T, bovine-438Y86F/I93T, AncEutheria-360F86S/T93I, AncBoreotheia-360T93I/S97T, bovine-438I93T/T97S, AncEutheria-360T93I/L116V, AncEutheria-360F86S/L116V, AncEutheria-360F86S/T93I/L116V, elephant-419S86F/I93T/V116L, AncBoreotheia-360F86Y/T93I/S97T and AncBoreotheria-360F46T/F49L/T52F/F86L/T93P/A114G/S118T.

The in vitro assay

Ancestral and other mutant opsins were expressed in COS1 cells by transient transfection [28]. The COS1 cell is a simian cell that is permissive for SV 40 replication and supports the replication of a recombinant DNA molecule containing a SV40 origin of replication that requires no ethical approvals or permissions. The pigments were regenerated by incubating the opsins with 11-cis-retinal (a gift from Dr. Rosalie K. Crouch at Storm Eye Institute, Medical University of South Carolina and National Eye Institutes) and were purified using immobilized 1D4 (The Culture Center, Minneapolis, MN) in buffer W1 (50 mM N-(2-hydroxyethyl) piperazine-N’-2-ethanesulfonic acid (HEPES) (pH 6.6), 140 mM NaCl, 3 mM MgCl2, 20 % (w/v) glycerol and 0.1 % dodecyl maltoside). UV visible spectra were recorded at 20 °C using a Hitachi U-3000 dual beam spectrophotometer. Visual pigments were bleached for 3 min using a 60 W standard light bulb equipped with a Kodak Wratten #3 filter at a distance of 20 cm. Data were analyzed using Sigmaplot software (Jandel Scientific, San Rafael, CA).

Linear models

To estimate the effects of individual and epistatic interactions of amino acid changes (θs) on the λmax-shift, we used linear statistical models, where the λmax of each visual pigment is expressed as a function of the λmax of AncEutheria-360 (or elephant-419) [87] and the effects of the appropriate single and multiple amino acid changes on the λmax-shift [95].

Modelling of protein structures

The initial structures of contemporary, ancestral and mutant pigments were obtained from their amino acid sequences with homology modelling (SWISS MODEL [34]) by taking the structure of bovine rhodopsin (pdb code: 1U19) as a template, which are called SWISS models. SWISS model structures of five representative pigments were further refined by 1) adding the missing hydrogen atoms and 2) full AMBER96 ( [24, 25] geometry optimizations, which are called AMBER models. In this protein modelling, we excluded the 11-cis-retinal not only because such protein structures can reveal whether or not the retinal can actually fit into the retinal-binding pocket but also because widely accessible protein modelling methods exclude the 11-cis-retinal. We also performed multivariable linear regression analysis between λmax values (dependent variable) and the structural parameters A, B, and A/B derived from the AMBER models (independent variable) [95].

Sequence analyses

The proportion (p) of different nucleotides was computed for each pairwise comparison and the number (d) of nucleotide substitutions per site was estimated from d = − (3/4) ln [1 – (4/3)p] [96]. For a set of violet opsin gene (A), UV opsin gene (B) and evolutionarily more distantly related third gene (C), the numbers of nucleotide substitutions per site between genes A and B (dAB), between genes A and C (dAC) and between genes B and C (dBC) were evaluated, and then the branch lengths for genes A and B after their divergence are determined by (1/2) (dAB + dAC – dBC) and (1/2) (dAB - dAC + dBC), respectively. The rates of nucleotide substitution for opsin genes A and B were evaluated by dividing their branch lengths by appropriate divergence times that were estimated from the timetree of life ( The divergence times between scabbardfish-423 and tilapia-360, between frog-423 and salamander-356, between chicken-415 and zebra finch-359, between human-414 and mouse-359, between bovine-438 and mouse-359, between squirrel-440 and mouse-359, between elephant-419 and mouse-359 and between wallaby-420 and dunnart-363 are taken as 114, 298, 102, 98, 98. 98, 105 and 62 MY ago, respectively ( Since F86 is missing from scabbardfish-423, 16 critical codon sites have been considered. Note that zebra finch-359 is listed under violet pigment because it went through the transition of UV pigment → violet pigment → UV pigment during evolution [1]. Standard errors for the evolutionary rates were estimated from [9p(1-p)/{(3–4p)2n}]1/2, where n is the number of nucleotide sites considered.

Ethics statement

All DNA sequence data were taken from publically available sources and no animal experiments were conducted in this study

Availability of supporting data

The authors confirm that all data underlying the findings are fully available without restriction. All relevant data are within the Methods and in the Additional files section.



Angstrom or 10−10 m


Assisted model building with energy refinement software package


The ancestral pigment of Amniotes


The ancestral pigment of Amphibians


The ancestral pigment of Aves


AncSauropsid-360 with mutations F49V/F86S/L116V/S118A


The ancestral pigment of Boreoeutherians


The ancestral pigment of Euteleosts, AncEutheria, the ancestral pigment of Eutherians


The ancestral pigment of the common ancestor of Placentals and Marsupials


The ancestral pigment of Sauropsids


The ancestral pigment of Tetrapods


The ancestral pigment of jawed vertebrates




The A/B ratio of an ancestral pigment subtracted from that of a mutant pigment


the λmax of an ancestral pigment subtracted from that of a mutant pigment

Δλmax :



Hydrogen-bond network



λelephant :

the λmax of elephant-419

λmax :

The wavelength of maximal absorption


Middle and long wavelength-sensitive




Phylogenetic analysis by maximum likelihood


Posterior probability


Pigments with protonated Schiff base nitrogen-linked retinals

RH1 pigment:



Pigments with unprotonated Schiff base nitrogen-linked retinals


Schiff base


A fully automated protein structure homology-modelling server


short wavelength-sensitive type 1


SWS type 2



θi :

the epistatic effect of mutation i on the λmax-shift


Whelan and Goldman


  1. 1.

    Shi Y, Yokoyama S. Molecular analysis of the evolutionary significance of ultraviolet vision in vertebrates. Proc Natl Acad Sci U S A. 2003;100(14):8308–13.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  2. 2.

    Fleishman LJ, Loew ER, Leal M. Ultraviolet vision in Lizard. Nature. 1993;365:397.

    Article  Google Scholar 

  3. 3.

    Viitala J, Korplmaki E, Palokangas P, Koivula M. Attraction of kestrels to vole scent marks visible in ultraviolet light. Nature. 1995;373:425–7.

    CAS  Article  Google Scholar 

  4. 4.

    Bennett ATD, Cuthill IC, Partridge JC, Maier EJ. Ultraviolet vision and mate choice in zebra finches. Nature. 1996;380:433–5.

    CAS  Article  Google Scholar 

  5. 5.

    Church SC, Bennett AT, Cuthill IC, Partridge JC. Ultraviolet cues affect the foraging behaviour of blue tits. Proc R Soc Lond B. 1998;265:1509–14.

    Article  Google Scholar 

  6. 6.

    Odeen A, Hastad O. The phylogenetic distribution of ultraviolet sensitivity in birds. MBMC Evol Biol. 2013;13:36.

    Article  Google Scholar 

  7. 7.

    Boettner EA, Wolter JR. Transmission of the ocular media. Invest Ophthalmol Vis Sci. 1962;1(6):776–83.

    Google Scholar 

  8. 8.

    Douglas RH, Jeffery G. The spectral transmission of ocular media suggests ultraviolet sensitivity is widespread among mammals. Proc R Soc Lond B Biol Sci. 2014;281:20132995.

    CAS  Article  Google Scholar 

  9. 9.

    Jacobs GH. Ultraviolet vision in vertebrates. Am Zool. 1992;32:544–54.

    Article  Google Scholar 

  10. 10.

    Bennett AT, Cuthill IC. Ultraviolet vision in birds: what is its function? Vision Res. 1994;34(11):1471–8.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Ebrey T, Koutalos Y. Vertebrate photoreceptors. Prog Retin Eye Res. 2001;20(1):49–94.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Ebrey TG, Takahashi Y. Photobiology of retinal proteins. In: Coohill TP, Valenzeno DP, editors. Photobiology for the 21st Century. Overland Park, KA: Valdenmar; 2002. p. 101–33.

    Google Scholar 

  13. 13.

    Hauser FE, Van Hazel I, Chang BSW. Spectral tuning in vertebrate short wavelength-sensitive 1 (SWS1) visual pigments: Can wavelength sensitivity be inferred from sequence data? J Exp Zool. 2014;322(7):529–39.

    CAS  Article  Google Scholar 

  14. 14.

    Yokoyama S. Synthesis of experimental molecular biology and evolutionary biology: An example from the world of vision. Bioscience. 2012;62(1):939–48.

    PubMed  PubMed Central  Google Scholar 

  15. 15.

    Yokoyama S. Synthetic biology of phenotypic adaptation in vertebrates: the next frontier. Mol Biol Evol. 2013;30(7):1495–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Yokoyama S. Molecular evolution of vertebrate visual pigments. Prog Retin Eye Res. 2000;19(4):385–419.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Yokoyama S. Evolution of dim-light and color visual pigments. Annu Rev Genom Human Genet. 2008;9:259–82.

    CAS  Article  Google Scholar 

  18. 18.

    Shi Y, Radlwimmer FB, Yokoyama S. Molecular genetics and the evolution of ultraviolet vision in vertebrates. Proc Natl Acad Sci U S A. 2001;98(20):11731–6.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Takahashi Y, Yokoyama S. Genetic basis of spectral tuning in the violet-sensitive visual pigment of African clawed frog, Xenopus laevis. Genetics. 2005;171(3):1153–60.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Yokoyama S. Molecular genetic basis of adaptive selection: examples from color vision in vertebrates. Annu Rev Genet. 1997;31:315–36.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Yokoyama S, Yokoyama R. Adaptive evolution of photoreceptors and visual pigments in vertebrates. Annu Rev Ecol Syst. 1996;27:543–67.

    Article  Google Scholar 

  22. 22.

    Wald G. Carotenoids and the visual cycle. J Gen Physiol. 1935;19:351–71.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Wald G. The photoreceptor process in vision. Am J Ophthalmol. 1955;40(5 Part 2):18–41.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Altun A, Yokoyama S, Morokuma K. Color tuning in short wavelength-sensitive human and mouse visual pigments: Ab initio Quantum Mechanics/Molecular Mechnics studies. J Phys Chem A. 2009;113(43):11685–92.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Altun A, Morokuma K, Yokoyama S. H-bond network around retinal regulates the evolution of ultraviolet and violet vision. ACS Chem Biol. 2011;6(8):775–80.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Tada T, Altun A, Yokoyama S. Evolutionary replacement of UV vision by violet vision in fish. Proc Natl Acad Sci U S A. 2009;106(41):17457–62.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007;24(8):1586–91.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Yokoyama S. Phylogenetic analysis and experimental approaches to study color vision in vertebrates. In: Palczewski K, editor. Methods Enzymol. 315th ed. San Diego: Academic; 2000. p. 312–25.

    Google Scholar 

  29. 29.

    Yokoyama S, Xing J, Liu Y, Faggionato D, Altun A, Starmer WT. Epistatic adaptive evolution of human color vision. PLoS Genet. 2014;10:e1004884.

    Article  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Yokoyama S, Tada T, Zhang H, Britt L. Elucidation of phenotypic adaptations: Molecular analyses of dim-light vision proteins in vertebrates. Proc Natl Acad Sci U S A. 2008;105(36):13480–5.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Yokoyama S, Radlwimmer FB. The molecular genetics and evolution of red and green color vision in vertebrates. Genetics. 2001;158(4):1697–710.

    CAS  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Yokoyama S, Radlwimmer FB, Blow NS. Ultraviolet pigments in birds evolved from violet pigments by a single amino acid change. Proc Natl Acad Sci U S A. 2000;97(13):7366–71.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Wilkie SE, Robinson PR, Cronin TW, Poopalasundaram S, Bowmaker JK, Hunt DM. Spectral tuning of avian violet- and ultraviolet-sensitive visual pigments. Biochemistry. 2000;39(27):7895–901.

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    Biasini M, Bienert S, Waterhouse A, Arnold K, Studer G, Schmidt T, Kiefer F, Gallo Cassarino T, Bertoni M, Bordoli L, et al. SWISS-MODEL: modelling protein tertiary and quaternary structure using evolutionary information. Nucleic Acids Res. 2014;42(W1):W252–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Yokoyama S, Altun A, Jia H, Yang H, Koyama T, Faggionato D, Liu Y, Starmer WT. Adaptive evolutionary paths from UV reception to sensing violet light by epistatic interactions. Sci Adv. 2015;1(8):e1500162.

    Article  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Takenaka N, Yokoyama S. Mechanisms of spectral tuning in the RH2 pigments of Tokay gecko and American chameleon. Gene. 2007;399(1):26–32.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Yokoyama S, Takenaka N, Blow N. A novel spectral tuning in the short wavelength-sensitive (SWS1 and SWS2) pigments of bluefin killifish (Lucania goodei). Gene. 2007;396(1):196–202.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  38. 38.

    Takahashi Y, Ebrey TG. Molecular basis of spectral tuning in the newt short wavelength sensitive visual pigment. Biochemistry. 2003;42(20):6025–34.

    CAS  Article  PubMed  Google Scholar 

  39. 39.

    Yokoyama S, Yang H, Starmer WT. Molecular basis of spectral tuning in the red- and green-sensitive (M/LWS) pigments in vertebrates. Genetics. 2008;179(4):2037–43.

    Article  PubMed  PubMed Central  Google Scholar 

  40. 40.

    Sun H, Macke JP, Nathans J. Mechanisms of spectral tuning in the mouse green cone pigment. Proc Natl Acad Sci U S A. 1997;94(16):8860–5.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Fasick JI, Lee N, Oprian DD. Spectral tuning in the human blue cone pigment. Biochemistry. 1999;38(36):11593–6.

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Fasick JI, Robsinson PR. Mechanism of spectral tuning in the dolphin visual pigments. Biochemistry. 1998;37(2):433–8.

    CAS  Article  PubMed  Google Scholar 

  43. 43.

    Hughes AL. The origin of adaptive phenotypes. Proc Natl Acad Sci U S A. 2008;105(36):13193–4.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  44. 44.

    Hughes AL, Nei M. Pattern of nucleotide substitution at major histocompatibility complex class I loci reveals overdominant selection. Nature. 1988;335(6186):167–70.

    CAS  Article  PubMed  Google Scholar 

  45. 45.

    Nei M, Kumar S. Molecular evolution and phylogenetics. Oxford: Oxford University Press; 2000.

    Google Scholar 

  46. 46.

    Nei M. Selectionism and neutralism in molecular evolution. Mol Biol Evol. 2005;22(12):2318–42.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  47. 47.

    Nei M. Mutation-driven evolution. Oxford: Oxford University Press; 2013.

    Google Scholar 

  48. 48.

    Zhang G, Li C, Li Q, Li B, Larkin DM, Lee C, Storz JF, Antunes A, Greenwold MJ, Meredith RW, et al. Comparative genomics reveals insights into avian genome evolution and adaptation. Science. 2014;346(6215):1311–20.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  49. 49.

    Nozawa M, Suzuki Y, Nei M. Reliabilities of identifying positive selection by the branch-site and the site-prediction methods. Proc Natl Acad Sci U S A. 2009;106(16):6700–5.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  50. 50.

    Nozawa M, Suzuki Y, Nei M: Response to Yang et al.: Problems with Bayesian methods of detecting positive selection at the DNA sequence level. Proc Natl Acad Sci U S A 2009, 106(36):10.1073/pnas.0906089106

  51. 51.

    Nozawa M, Suzuki Y, Nei M. Is positive selection responsible for the evolution of a duplicate UV-sensitive opsin gene in Heliconius butterflies? Proc Natl Acad Sci U S A. 2010;107(23):E96. author reply E97.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  52. 52.

    Yang Z, dos Reis M. Statistical properties of the branch-site test of positive selection. Mol Biol Evol. 2011;28(3):1217–28.

    CAS  Article  PubMed  Google Scholar 

  53. 53.

    Yang Z, Nielsen R, Goldman N: In defense of statistical methods for detecting positive selection. Proc Natl Acad Sci U S A 2009, 106(36):10.1073/pnas.0904550106

  54. 54.

    Nei M, Suzuki Y, Nozawa M. The neutral theory of molecular evolution in the genomic era. Annu Rev Genomics Hum Genet. 2010;11:265–89.

    CAS  Article  PubMed  Google Scholar 

  55. 55.

    Jermann TM, Opitz JG, Stackhouse J, Benner SA. Reconstructing the evolutionary history of the artiodactyl ribonuclease superfamily. Nature. 1995;374(6517):57–9.

    CAS  Article  PubMed  Google Scholar 

  56. 56.

    Zhang J, Zhang YP, Rosenberg HF. Adaptive evolution of a duplicated pancreatic ribonuclease gene in a leaf-eating monkey. Nat Genet. 2002;30(4):411–5.

    CAS  Article  PubMed  Google Scholar 

  57. 57.

    Perutz MF. Species adaptation in a protein molecule. Mol Biol Evol. 1983;1(1):1–28.

    CAS  PubMed  Google Scholar 

  58. 58.

    Jessen TH, Weber RE, Fermi G, Tame J, Braunitzer G. Adaptation of bird hemoglobins to high altitudes: demonstration of molecular mechanism by protein engineering. Proc Natl Acad Sci U S A. 1991;88(15):6519–22.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  59. 59.

    Weber RE, Jessen TH, Malte H, Tame J. Mutant hemoglobins (alpha 119-Ala and beta 55-Ser): functions related to high-altitude respiration in geese. J Appl Physiol. 1993;75(6):2646–55.

    CAS  PubMed  Google Scholar 

  60. 60.

    Natarajan C, Inoguchi N, Weber RE, Fago A, Moriyama H, Storz JF. Epistasis among adaptive mutations in deer mouse hemoglobin. Science. 2013;340:1324–7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  61. 61.

    Eick GN, Colucci JK, Harms MJ, Ortlund EA, Thornton JW. Evolution of minimal specificity and promiscuity in steroid hormone receptors. PLoS Genet. 2012;8(11):e1003072.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  62. 62.

    Harms MJ, Thornton JW. Analyzing protein structure and function using ancestral gene reconstruction. Curr Opin Struct Biol. 2010;20(3):360–6.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  63. 63.

    Harms MJ, Thornton JW. Evolutionary biochemistry: revealing the historical and physical causes of protein properties. Nat Rev Genet. 2013;14:559–71.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  64. 64.

    Harms MJ, Eick GN, Goswami D, Colucci JK, Griffin PR, Ortlund EA, Thorntona JW. Biophysical mechanisms for large-effect mutations in the evolution of steroid hormone receptors. Proc Natl Acad Sci USA. 2013;110(28):11475–80.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  65. 65.

    Zhuang H, Chien MS, Matsunami H. Dynamic functional evolution of an odorant receptor for sex-steroid-derived odors in primates. Proc Natl Acad Sci U S A. 2009;106(50):21247–51.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  66. 66.

    Adipietro KA, Mainland JD, Matsunami H. Functional evolution of mammalian odorant receptors. PLoS Genet. 2012;8(7):e1002821.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  67. 67.

    Mainland JD, Keller A, Li Y, Zhou T, Trimmer C, Snyder LL, Moberly AH, Adipietro KA, Liu WLL, Zhuang H, et al. The missense of smell: functional variability in the human odorant receptor repertoire. Nat Neuroscience. 2014;17:114–20.

    CAS  Article  PubMed  Google Scholar 

  68. 68.

    Carroll SB. The making of the fittest. New York: WW Norton & Co; 2006.

    Google Scholar 

  69. 69.

    Hunt DM, Carvalho LS, Cowing JA, Parry JW, Wilkie SE, Davies WL, Bowmaker JK. Spectral tuning of shortwave-sensitive visual pigments in vertebrates. Photochem Photobiol. 2007;83(2):303–10.

    CAS  Article  PubMed  Google Scholar 

  70. 70.

    Yokoyama S, Radlwimmer FB, Kawamura S. Regeneration of ultraviolet pigments of vertebrates. FEBS Lett. 1998;423(2):155–8.

    CAS  Article  PubMed  Google Scholar 

  71. 71.

    Chinen A, Hamaoka T, Yamada Y, Kawamura S. Gene duplication and spectral diversification of cone visual pigments of zebrafish. Genetics. 2003;163(2):663–75.

    CAS  PubMed  PubMed Central  Google Scholar 

  72. 72.

    Spady TC, Parry JW, Robinson PR, Hunt DM, Bowmaker JK, Carleton KL. Evolution of the cichlid visual palette through ontogenetic subfunctionalization of the opsin gene arrays. Mol Biol Evol. 2006;23(8):1538–47.

    CAS  Article  PubMed  Google Scholar 

  73. 73.

    Parry JW, Carleton KL, Spady T, Carboo A, Hunt DM, Bowmaker JK. Mix and match color vision: tuning spectral sensitivity by differential opsin gene expression in Lake Malawi cichlids. Curr Biol. 2005;15(19):1734–9.

    CAS  Article  PubMed  Google Scholar 

  74. 74.

    Matsumoto Y, Fukamachi S, Mitani H, Kawamura S. Functional characterization of visual opsin repertoire in Medaka (Oryzias latipes). Gene. 2006;371(2):268–78.

    CAS  Article  PubMed  Google Scholar 

  75. 75.

    Pointer MA, Cheng CH, Bowmaker JK, Parry JW, Soto N, Jeffery G, Cowing JA, Hunt DM. Adaptations to an extreme environment: retinal organisation and spectral properties of photoreceptors in Antarctic notothenioid fish. J Exp Biol. 2005;208(Pt 12):2363–76.

    CAS  Article  PubMed  Google Scholar 

  76. 76.

    Ma JX, Kono M, Xu L, Das J, Ryan JC, Hazard 3rd ES, Oprian DD, Crouch RK. Salamander UV cone pigment: sequence, expression, and spectral properties. Vis Neurosci. 2001;18(3):393–9.

    CAS  Article  PubMed  Google Scholar 

  77. 77.

    Das D, Wilkie SE, Hunt DM, Bowmaker JK. Visual pigments and oil droplets in the retina of a passerine bird, the canary Serinus canaria: microspectrophotometry and opsin sequences. Vision Res. 1999;39(17):2801–15.

    CAS  Article  PubMed  Google Scholar 

  78. 78.

    Kawamura S, Yokoyama S. Functional characterization of visual and nonvisual pigments of American chameleon (Anolis carolinensis). Vision Res. 1998;38(1):37–44.

    CAS  Article  PubMed  Google Scholar 

  79. 79.

    Loew ER. A third, ultraviolet-sensitive, visual pigment in the Tokay gecko (Gekko gekko). Vision Res. 1994;34(11):1427–31.

    CAS  Article  PubMed  Google Scholar 

  80. 80.

    Davies WL, Cowing JA, Bowmaker JK, Carvalho LS, Gower DJ, Hunt DM. Shedding light on serpent sight: the visual pigments of henophidian snakes. J Neurosci. 2009;29(23):7519–25.

    CAS  Article  PubMed  Google Scholar 

  81. 81.

    Bowmaker JK, Dartnall HJ, Mollon JD. Microspectrophotometric demonstration of four classes of photoreceptor in an old world primate, Macaca fascicularis. J Physiol. 1980;298:131–43.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  82. 82.

    Bowmaker JK, Jacobs GH, Mollon JD. Polymorphism of photopigments in the squirrel monkey: a sixth phenotype. Proc R Soc Lond B Biol Sci. 1987;231(1264):383–90.

    CAS  Article  PubMed  Google Scholar 

  83. 83.

    Travis DS, Bowmaker JK, Mollon JD. Polymorphism of visual pigments in a callitrichid monkey. Vision Res. 1988;28(4):481–90.

    CAS  Article  PubMed  Google Scholar 

  84. 84.

    Fasick JI, Applebury ML, Oprian DD. Spectral tuning in the mammalian short-wavelength sensitive cone pigments. Biochemistry. 2002;41(21):6860–5.

    CAS  Article  PubMed  Google Scholar 

  85. 85.

    Carvalho LS, Cowing JA, Wilkie SE, Bowmaker JK, Hunt DM. Shortwave visual sensitivity in tree and flying squirrels reflects changes in lifestyle. Curr Biol. 2006;16(3):R81–3.

    CAS  Article  Google Scholar 

  86. 86.

    Parry JW, Poopalasundaram S, Bowmaker JK, Hunt DM. A novel amino acid substitution is responsible for spectral tuning in a rodent violet-sensitive visual pigment. Biochemistry. 2004;43(25):8014–20.

    CAS  Article  PubMed  Google Scholar 

  87. 87.

    Yokoyama S, Takenaka N, Agnew DW, Shoshani J. Elephants and human color-blind deuteranopes have identical sets of visual pigments. Genetics. 2005;170(1):335–44.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  88. 88.

    Deeb SS, Wakefield MJ, Tada T, Marotte L, Yokoyama S, Marshall Graves JA. The cone visual pigments of an Australian marsupial, the tammar wallaby (Macropus eugenii): sequence, spectral tuning, and evolution. Mol Biol Evol. 2003;20(10):1642–9.

    CAS  Article  PubMed  Google Scholar 

  89. 89.

    Cowing JA, Arrese CA, Davies WL, Beazley LD, Hunt DM. Cone visual pigments in two marsupial species: the fat-tailed dunnart (Sminthopsis crassicaudata) and the honey possum (Tarsipes rostratus). Proc Biol Sci. 2008;275(1642):1491–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  90. 90.

    Yokoyama S, Blow NS, Radlwimmer FB. Molecular evolution of color vision of zebra finch. Gene. 2000;259(1–2):17–24.

    CAS  Article  PubMed  Google Scholar 

  91. 91.

    Okano T, Kojima D, Fukada Y, Shichida Y, Yoshizawa T. Primary structures of chicken cone visual pigments: vertebrate rhodopsins have evolved out of cone visual pigments. Proc Natl Acad Sci U S A. 1992;89(13):5932–6.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  92. 92.

    Cowing JA, Poopalasundaram S, Wilkie SE, Robinson PR, Bowmaker JK, Hunt DM. The molecular mechanism for the spectral shifts between vertebrate ultraviolet- and violet-sensitive cone visual pigments. Biochem J. 2002;367(Pt 1):129–35.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  93. 93.

    Dukkipati A, Vought BW, Singh D, Birge RR, Knox BE. Serine 85 in transmembrane helix 2 of short-wavelength visual pigments interacts with the retinylidene Schiff base counterion. Biochemistry. 2001;40(50):15098–108.

    CAS  Article  PubMed  Google Scholar 

  94. 94.

    Babu KR, Dukkipati A, Birge RR, Knox BE. Regulation of phototransduction in short-wavelength cone visual pigments via the retinylidene Schiff base counterion. Biochemistry. 2001;40(46):13760–6.

    CAS  Article  PubMed  Google Scholar 

  95. 95.

    Graybill F. An introduction to linear statistical models, vol. 1. New York: McGraw-Hill; 1961.

    Google Scholar 

  96. 96.

    Jukes TH, Cantor CR. Evolution of protein molecules. In: Munro HN, editor. Mammalian protein metabolism. New York: Academic; 1969. p. 21–132.

    Google Scholar 

  97. 97.

    Palczewski K. G protein-coupled receptor rhodopsin. Annu Rev Biochem. 2006;75:743–67.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

Download references


We thank Dr. Rosalie K. Crouch (Storm Eye Institute, Medical University of South Carolina) and National Eye Institute for supplying 11-cis-retinal and comments by Dr. Ruth Yokoyama and anonymous reviewers are greatly appreciated. S.Y. is supported by National Institutes of Health (EY016400) and Emory University.

Author information



Corresponding author

Correspondence to Shozo Yokoyama.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

SY designed, organized, performed SWISS MODEL and molecular evolutionary analyses; AA performed AMBER geometry optimization and SWISS model analyses; TT, YL and DF performed experiments; SY and AA wrote the paper. All authors have read and approved the final version of the manuscript.

Additional files

Additional file 1: Table S1.

Amino acids of AncVertebrate-361 with PP < 0.95 (in parentheses) inferred using PAML with JTT and WAG models. (DOCX 36 kb)

Additional file 2: Figure S1.

The amino acid sequences of 10 ancestral pigments. Ancestral pigments are reconstructed by introducing all different amino acids between sites 31 and 312 (see Methods). Amino acids with PP < 0.95 are indicated by bold italics. Four critical amino acids that caused the evolution of AncBird-393 from AncSauropsid-360 are indicated by black boxes. (PDF 33 kb)

Additional file 3: Figure S2.

The absorption spectra of seven ancestral pigments. (PDF 451 kb)

Additional file 4: Table S2.

The λmaxs and A/B ratios of HBNs of ancestral and present-day pigments. (DOCX 43 kb)

Additional file 5: Table S3.

Variable λmax-shifts and A/B ratios caused by various mutations. (DOCX 37 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Yokoyama, S., Tada, T., Liu, Y. et al. A simple method for studying the molecular mechanisms of ultraviolet and violet reception in vertebrates. BMC Evol Biol 16, 64 (2016).

Download citation


  • Visual pigments
  • UV and violet reception
  • Spectral tuning
  • Mutagenesis analyses
  • Hydrogen-bond network