Elucidation of phenotypic adaptations: Molecular analyses of dim-light vision proteins in vertebrates

A single amino acid substitution confers enhanced methylation activity of mammalian Dnmt3b on chromatin DNA

Dnmt3a and Dnmt3b are paralogous enzymes responsible for de novo DNA methylation but with distinguished biological functions. In mice, disruption of Dnmt3b but not Dnmt3a causes global DNA hypomethylation, especially in repetitive sequences, which comprise the large majority of methylated DNA in the genome. By measuring DNA methylation activity of Dnmt3a and Dnmt3b homologues from five species, we found that mammalian Dnmt3b possessed significantly higher methylation activity on chromatin DNA than Dnmt3a and non-mammalian Dnmt3b. Sequence comparison and mutagenesis experiments identified a single amino acid substitution (I662N) in mammalian Dnmt3b as being crucial for its high chromatin DNA methylation activity. Further mechanistic studies demonstrated this substitution markedly enhanced the binding of Dnmt3b to nucleosomes and hence increased the chromatin DNA methylation activity. Moreover, this substitution was crucial for Dnmt3b to efficiently methylate repetitive sequences, which increased dramatically in mammalian genomes. Consistent with our observation that Dnmt3b evolved more rapidly than Dnmt3a during the emergence of mammals, these results demonstrated that the I662N substitution in mammalian Dnmt3b conferred enhanced chromatin DNA methylation activity and contributed to functional adaptation in the epigenetic system.

The neutral theory of molecular evolution in the genomic era

The neutral theory of molecular evolution has been widely accepted and is the guiding principle for studying evolutionary genomics and the molecular basis of phenotypic evolution. Recent data on genomic evolution are generally consistent with the neutral theory. However, many recently published papers claim the detection of positive Darwinian selection via the use of new statistical methods. Examination of these methods has shown that their theoretical bases are not well established and often result in high rates of false-positive and false-negative results. When the deficiencies of these statistical methods are rectified, the results become largely consistent with the neutral theory. At present, genome-wide analyses of natural selection consist of collections of single-locus analyses. However, because phenotypic evolution is controlled by the interaction of many genes, the study of natural selection ought to take such interactions into account. Experimental studies of evolution will also be crucial.

CodonTest: modeling amino acid substitution preferences in coding sequences

Codon models of evolution have facilitated the interpretation of selective forces operating on genomes. These models, however, assume a single rate of non-synonymous substitution irrespective of the nature of amino acids being exchanged. Recent developments have shown that models which allow for amino acid pairs to have independent rates of substitution offer improved fit over single rate models. However, these approaches have been limited by the necessity for large alignments in their estimation. An alternative approach is to assume that substitution rates between amino acid pairs can be subdivided into rate classes, dependent on the information content of the alignment. However, given the combinatorially large number of such models, an efficient model search strategy is needed. Here we develop a Genetic Algorithm (GA) method for the estimation of such models. A GA is used to assign amino acid substitution pairs to a series of rate classes, where is estimated from the alignment. Other parameters of the phylogenetic Markov model, including substitution rates, character frequencies and branch lengths are estimated using standard maximum likelihood optimization procedures. We apply the GA to empirical alignments and show improved model fit over existing models of codon evolution. Our results suggest that current models are poor approximations of protein evolution and thus gene and organism specific multi-rate models that incorporate amino acid substitution biases are preferred. We further anticipate that the clustering of amino acid substitution rates into classes will be biologically informative, such that genes with similar functions exhibit similar clustering, and hence this clustering will be useful for the evolutionary fingerprinting of genes.

Evolution in protein-coding DNA sequences can be modeled at three levels: nucleotides, amino acids or codons that encode the amino acids. Codon models incorporate nucleotide and amino acid information, and allow the estimation of the rate at which amino acids are replaced () versus the rate at which they are preserved (). The ratio has been used in thousands of studies to detect molecular footprints of natural selection. A serious limitation of most codon models is the unrealistic assumption that all non-synonymous substitutions occur at the same rate. Indeed, amino acid models have consistently demonstrated that different residues are exchanged more or less frequently, depending on incompletely understood factors. We derive and validate a computational approach for inferring codon models which combine the power to investigate natural selection with data-driven amino acid substitution biases from alignments. The addition of amino acid properties can lead to more powerful and accurate methods for studying natural selection and the evolutionary history of protein-coding sequences. The pattern of amino acid substitutions specific to a given alignment can be used to compare and contrast the evolutionary properties of different genes, providing an evolutionary analog to protein family comparisons.

Diversifying selection and functional analysis of interleukin-4 suggests antagonism-driven evolution at receptor-binding interfaces

Background

Interleukin-4 (IL4) is a secreted immunoregulatory cytokine critically involved in host protection from parasitic helminths [1]. Reasoning that helminths may have evolved mechanisms to antagonize IL4 to maximize their dispersal, we explored mammalian IL4 evolution.

Results

This analysis revealed evidence of diversifying selection at 15 residues, clustered in epitopes responsible for IL4 binding to its Type I and Type II receptors. Such a striking signature of selective pressure suggested either recurrent episodes of pathogen antagonism or ligand/receptor co-evolution. To test the latter possibility, we performed detailed functional analysis of IL4 allotypes expressed by Mus musculus musculus and Mus musculus castaneus, which happen to differ at 5 residues (including three at positively selected sites) in and adjacent to the site 1 epitope that binds the IL4Rα subunit shared by the Type I and Type II IL4 receptors. We show that this intra-species variation affects the ability of IL4 neither to bind IL4 receptor alpha (IL4Rα) nor to signal biological responses through its Type I receptor.

Conclusions

Our results -- reminiscent of clustered positively selected sites revealing functionally important residues at host-virus interaction interfaces -- are consistent with IL4 having evolved to avoid recurrent pathogen antagonism, while maintaining the capacity to bind and signal through its cognate receptor. This work exposes what may be a general feature of evolutionary conflicts fought by pathogen antagonists at host protein-protein interaction interfaces involved in immune signaling: the emergence of receptor-binding ligand epitopes capable of buffering amino acid variation.

Benchmarking multi-rate codon models

The single rate codon model of non-synonymous substitution is ubiquitous in phylogenetic modeling. Indeed, the use of a non-synonymous to synonymous substitution rate ratio parameter has facilitated the interpretation of selection pressure on genomes. Although the single rate model has achieved wide acceptance, we argue that the assumption of a single rate of non-synonymous substitution is biologically unreasonable, given observed differences in substitution rates evident from empirical amino acid models. Some have attempted to incorporate amino acid substitution biases into models of codon evolution and have shown improved model performance versus the single rate model. Here, we show that the single rate model of non-synonymous substitution is easily outperformed by a model with multiple non-synonymous rate classes, yet in which amino acid substitution pairs are assigned randomly to these classes. We argue that, since the single rate model is so easy to improve upon, new codon models should not be validated entirely on the basis of improved model fit over this model. Rather, we should strive to both improve on the single rate model and to approximate the general time-reversible model of codon substitution, with as few parameters as possible, so as to reduce model over-fitting. We hint at how this can be achieved with a Genetic Algorithm approach in which rate classes are assigned on the basis of sequence information content.

The p1D4-hrGFP II expression vector: a tool for expressing and purifying visual pigments and other G protein-coupled receptors

The heterologous expression of membrane proteins such as G protein-coupled receptors can be a notoriously difficult task. We have engineered an expression vector, p1D4-hrGFP II, in order to efficiently express visual pigments in mammalian cell culture. This expression vector is based on pIRES-hrGFP II (Stratagene), with the addition of a C-terminal 1D4 epitope tag for immunoblotting and immunoaffinity purification. This vector employs the CMV promoter and hrGFP II, a co-translated reporter gene. We measured the effectiveness of pIRES-hrGFP II in expressing bovine rhodopsin, and showed a 3.9- to 5.7-fold increase in expression as measured by absorbance spectroscopy as compared with the pMT vector, a common choice for visual pigment expression. We then expressed zebrafish RH2-1 using p1D4-hrGFP II in order to assess its utility in expressing cone opsins, known to be less stable and more difficult to express than bovine rhodopsin. We show a λ(280)/λ(MAX) value of 3.3, one third of that reported in previous studies, suggesting increased expression levels and decreased levels of misfolded, non-functional visual pigment. Finally, we monitored HEK293T cell growth following transfection with pIRES-hrGFP II using fluorescence microscopy to illustrate the benefits of having a co-translated reporter during heterologous expression studies.

Allelic variation in Malawi cichlid opsins: a tale of two genera

The role of sequence variation in the spectral tuning of color vision is well established in many systems. This includes the cichlids of Lake Victoria where sequence variation has been linked to environmental light gradients and speciation. The cichlids of Lake Malawi are a similar model for visual evolution, but the role of gene sequence variation in visual tuning between closely related species is unknown. This work describes such variation in multiple species of two rock-dwelling genera: Metriaclima and Labidochromis. Genomic DNA for seven cone opsin genes was sequenced and the structure of the opsin proteins was inferred. Retinal binding pocket polymorphisms were identified and compared to available data regarding spectral absorbance shifts. Sequence variation with known or potential effects on absorbance spectra were found in four genes: SWS1 (UV sensitive), SWS2B (violet sensitive), RH2Abeta (green sensitive), and LWS (red sensitive). Functional variation was distributed such that each genus had both a variable short-wavelength and long-wavelength sensitive opsin. This suggests spectral tuning is important at the margins of the cichlid visual spectrum. Further, there are two SWS1 opsin alleles that differ in sensitivity by 10 nm and are >2 MY divergent. One of these occurs in a haplotype block >1 kb. Potential haplotype blocks were found around the RH2 opsin loci. These data suggest that molecular diversification has resulted in functionally unique alleles and changes to the visual system. These data also suggest that opsin sequence variation tunes spectral sensitivities between closely related species and that the specific regions of spectral tuning are genus-specific.

Contrasting modes of evolution of the visual pigments in Heliconius butterflies

The adult compound eyes of passion-vine butterflies in the genus Heliconius contain one more UV opsin than other butterflies. Together with an 11-cis-3-hydroxyretinal chromophore, their four opsin genes UVRh1, UVRh2, BRh, and LWRh produce four rhodopsins that are UV-, blue-, or long wavelength absorbing. One of the Heliconius UV opsin genes, UVRh2, was found to have evolved under positive selection following recent gene duplication, using the branch-site test of selection. Using a more conservative test, the small-sample method, we confirm our prior finding of positive selection of UVRh2 and provide new statistical evidence of episodic evolution, that is, positive selection followed by purifying selection. We also newly note that one of the positively selected amino acid sites contains substitutions with known spectral tuning effects in avian ultraviolet- and violet-sensitive visual pigments. As this is one of a handful of described examples of positive selection of any specific gene in any butterfly where functional variation between copies has been characterized, we were interested in examining the molecular and physiological context of this adaptive event by examining the UV opsin genes in contrast to the other visual pigment genes. We cloned BRh and LWRh from 13 heliconiine species and UVRh1 and UVRh2 from Heliconius elevatus. In parallel, we performed in vivo epi-microspectrophotometric experiments to estimate the wavelength of peak absorbance, λ(max), of several rhodopsins in seven heliconiine species. In contrast to UVRh2, we found both physiological and statistical evidence consistent with purifying selection on UVRh1, BRh, and LWRh along the branch leading to the common ancestor of Heliconius. These results underscore the utility of combining molecular and physiological experiments in a comparative context for strengthening evidence for adaptive evolution at the molecular level.

Analyzing protein structure and function using ancestral gene reconstruction

Protein families with functionally diverse members can illuminate the structural determinants of protein function and the process by which protein structure and function evolve. To identify the key amino acid changes that differentiate one family member from another, most studies have taken a horizontal approach, swapping candidate residues between present-day family members. This approach has often been stymied, however, by the fact that shifts in function often require multiple interacting mutations; chimeric proteins are often nonfunctional, either because one lineage has amassed mutations that are incompatible with key residues that conferred a new function on other lineages, or because it lacks mutations required to support those key residues. These difficulties can be overcome by using a vertical strategy, which reconstructs ancestral genes and uses them as the appropriate background in which to study the effects of historical mutations on functional diversification. In this review, we discuss the advantages of the vertical strategy and highlight several exemplary studies that have used ancestral gene reconstruction to reveal the molecular underpinnings of protein structure, function, and evolution.

Mutation-selection models of coding sequence evolution with site-heterogeneous amino acid fitness profiles

Modeling the interplay between mutation and selection at the molecular level is key to evolutionary studies. To this end, codon-based evolutionary models have been proposed as pertinent means of studying long-range evolutionary patterns and are widely used. However, these approaches have not yet consolidated results from amino acid level phylogenetic studies showing that selection acting on proteins displays strong site-specific effects, which translate into heterogeneous amino acid propensities across the columns of alignments; related codon-level studies have instead focused on either modeling a single selective context for all codon columns, or a separate selective context for each codon column, with the former strategy deemed too simplistic and the latter deemed overparameterized. Here, we integrate recent developments in nonparametric statistical approaches to propose a probabilistic model that accounts for the heterogeneity of amino acid fitness profiles across the coding positions of a gene. We apply the model to a dozen real protein-coding gene alignments and find it to produce biologically plausible inferences, for instance, as pertaining to site-specific amino acid constraints, as well as distributions of scaled selection coefficients. In their account of mutational features as well as the heterogeneous regimes of selection at the amino acid level, the modeling approaches studied here can form a backdrop for several extensions, accounting for other selective features, for variable population size, or for subtleties of mutational features, all with parameterizations couched within population-genetic theory.

Retracing evolution of red fluorescence in GFP-like proteins from Faviina corals

Proteins of the green fluorescent protein family represent a convenient experimental model to study evolution of novelty at the molecular level. Here, we focus on the origin of Kaede-like red fluorescent proteins characteristic of the corals of the Faviina suborder. We demonstrate, using an original approach involving resurrection and analysis of the library of possible evolutionary intermediates, that it takes on the order of 12 mutations, some of which strongly interact epistatically, to fully recapitulate the evolution of a red fluorescent phenotype from the ancestral green. Five of the identified mutations would not have been found without the help of ancestral reconstruction, because the corresponding site states are shared between extant red and green proteins due to their recent descent from a dual-function common ancestor. Seven of the 12 mutations affect residues that are not in close contact with the chromophore and thus must exert their effect indirectly through adjustments of the overall protein fold; the relevance of these mutations could not have been anticipated from the purely theoretical analysis of the protein's structure. Our results introduce a powerful experimental approach for comparative analysis of functional specificity in protein families even in the cases of pronounced epistasis, provide foundation for the detailed studies of evolutionary trajectories leading to novelty and complexity, and will help rational modification of existing fluorescent labels.

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.

Evolutionary dynamics of rhodopsin type 2 opsins in vertebrates

Among the five groups of visual pigments in vertebrates, the rhodopsin type 2 (RH2) group shows the largest number of gene duplication events. We have isolated three intact and one nonfunctional RH2 opsin genes each from Northern lampfish (Stenobrachius leucopsarus) and scabbardfish (Lepidopus fitchi). Using the deduced amino acid sequences of these and other representative RH2 opsin genes in vertebrates, we have estimated the divergence times and evolutionary rates of amino acid substitution at various stages of RH2 opsin evolution. The results show that the duplications of the lampfish and scabbardfish RH2 opsins have occurred approximately 60 and approximately 30 million years ago (Ma), respectively. The evolutionary rates of RH2 opsins in the early vertebrate ancestors were approximately 0.25 x 10(-9)/site/year, which increased to approximately 1 x 10(-9) to 3 x 10(-9)/site/year in euteleost lineages and to approximately 0.3 x 10(-9) to 0.5 x 10(-9)/site/year in coelacanth and tetrapods.

Molecular mechanism of long-range synergetic color tuning between multiple amino acid residues in conger rhodopsin

The synergetic effects of multiple rhodopsin mutations on color tuning need to be completely elucidated. Systematic genetic studies and spectroscopy have demonstrated an interesting example of synergetic color tuning between two amino acid residues in conger rhodopsin's ancestral pigment (p501): -a double mutation at one nearby and one distant residue led to a significant λ(max) blue shift of 13 nm, whereas neither of the single mutations at these two sites led to meaningful shifts.To analyze the molecular mechanisms of this synergetic color tuning, we performed homology modeling, molecular simulations, and electronic state calculations. For the double mutant, N195A/A292S, in silico mutation analysis demonstrated conspicuous structural changes in the retinal chromophore, whereas that of the single mutant, A292S, was almost unchanged. Using statistical ensembles of QM/MM optimized structures, the excitation energy of retinal chromophore was evaluated for the three visual pigments. As a result, the λ(max) shift of double mutant (DM) from p501 was -8 nm, while that of single mutant (SM) from p501 was +1 nm. Molecular dynamics simulation for DM demonstrated frequent isomerization between 6-s-cis and 6-s-trans conformers. Unexpectedly, however, the two conformers exhibited almost identical excitation energy, whereas principal component analysis (PCA) identified the retinal-counterion cooperative change of BLA (bond length alternation) and retinal-counterion interaction lead to the shift.

Rhodopsin molecular evolution in mammals inhabiting low light environments

The ecological radiation of mammals to inhabit a variety of light environments is largely attributed to adaptive changes in their visual systems. Visual capabilities are conferred by anatomical features of the eyes as well as the combination and properties of their constituent light sensitive pigments. To test whether evolutionary switches to different niches characterized by dim-light conditions coincided with molecular adaptation of the rod pigment rhodopsin, we sequenced the rhodopsin gene in twenty-two mammals including several bats and subterranean mole-rats. We compared these to thirty-seven published mammal rhodopsin sequences, from species with divergent visual ecologies, including nocturnal, diurnal and aquatic groups. All taxa possessed an intact functional rhodopsin; however, phylogenetic tree reconstruction recovered a gene tree in which rodents were not monophyletic, and also in which echolocating bats formed a monophyletic group. These conflicts with the species tree appear to stem from accelerated evolution in these groups, both of which inhabit low light environments. Selection tests confirmed divergent selection pressures in the clades of subterranean rodents and bats, as well as in marine mammals that live in turbid conditions. We also found evidence of divergent selection pressures among groups of bats with different sensory modalities based on vision and echolocation. Sliding window analyses suggest most changes occur in transmembrane domains, particularly obvious within the pinnipeds; however, we found no obvious pattern between photopic niche and predicted spectral sensitivity based on known critical amino acids. This study indicates that the independent evolution of rhodopsin vision in ecologically specialised groups of mammals has involved molecular evolution at the sequence level, though such changes might not mediate spectral sensitivity directly.

Dynamic functional evolution of an odorant receptor for sex-steroid-derived odors in primates

Odorant receptors are among the fastest evolving genes in animals. However, little is known about the functional changes of individual odorant receptors during evolution. We have recently demonstrated a link between the in vitro function of a human odorant receptor, OR7D4, and in vivo olfactory perception of 2 steroidal ligands--androstenone and androstadienone--chemicals that are shown to affect physiological responses in humans. In this study, we analyzed the in vitro function of OR7D4 in primate evolution. Orthologs of OR7D4 were cloned from different primate species. Ancestral reconstruction allowed us to reconstitute additional putative OR7D4 orthologs in hypothetical ancestral species. Functional analysis of these orthologs showed an extremely diverse range of OR7D4 responses to the ligands in various primate species. Functional analysis of the nonsynonymous changes in the Old World Monkey and Great Ape lineages revealed a number of sites causing increases or decreases in sensitivity. We found that the majority of the functionally important residues in OR7D4 were not predicted by the maximum likelihood analysis detecting positive Darwinian selection.

The eyes have it: regulatory and structural changes both underlie cichlid visual pigment diversity

A major goal of evolutionary biology is to unravel the molecular genetic mechanisms that underlie functional diversification and adaptation. We investigated how changes in gene regulation and coding sequence contribute to sensory diversification in two replicate radiations of cichlid fishes. In the clear waters of Lake Malawi, differential opsin expression generates diverse visual systems, with sensitivities extending from the ultraviolet to the red regions of the spectrum. These sensitivities fall into three distinct clusters and are correlated with foraging habits. In the turbid waters of Lake Victoria, visual sensitivity is constrained to longer wavelengths, and opsin expression is correlated with ambient light. In addition to regulatory changes, we found that the opsins coding for the shortest- and longest-wavelength visual pigments have elevated numbers of potentially functional substitutions. Thus, we present a model of sensory evolution in which both molecular genetic mechanisms work in concert. Changes in gene expression generate large shifts in visual pigment sensitivity across the collective opsin spectral range, but changes in coding sequence appear to fine-tune visual pigment sensitivity at the short- and long-wavelength ends of this range, where differential opsin expression can no longer extend visual pigment sensitivity.

Rhodopsin population genetics and local adaptation: variable dim-light vision in sand gobies is illuminated

The visual pigments of fish are thought to be adapted to the variable spectral qualities of aquatic light environments. Most research on the role of natural selection on the evolution of rhodopsins and dim-light vision in fish has focused on variation among species and higher taxa. In this issue, Larmuseau et al. reveal substantial intraspecific sequence variation in RH1 (the rhodopsin gene) in sand gobies (Pomatoschistus minutus). Using population genetics and molecular evolution approaches, they detect positive selection on RH1 and find evidence for adaptation to local light conditions.

Evolutionary replacement of UV vision by violet vision in fish

The vertebrate ancestor possessed ultraviolet (UV) vision and many species have retained it during evolution. Many other species switched to violet vision and, then again, some avian species switched back to UV vision. These UV and violet vision are mediated by short wavelength-sensitive (SWS1) pigments that absorb light maximally (lambda(max)) at approximately 360 and 390-440 nm, respectively. It is not well understood why and how these functional changes have occurred. Here, we cloned the pigment of scabbardfish (Lepidopus fitchi) with a lambda(max) of 423 nm, an example of violet-sensitive SWS1 pigment in fish. Mutagenesis experiments and quantum mechanical/molecular mechanical (QM/MM) computations show that the violet-sensitivity was achieved by the deletion of Phe-86 that converted the unprotonated Schiff base-linked 11-cis-retinal to a protonated form. The finding of a violet-sensitive SWS1 pigment in scabbardfish suggests that many other fish also have orthologous violet pigments. The isolation and comparison of such violet and UV pigments in fish living in different ecological habitats will open an unprecedented opportunity to elucidate not only the molecular basis of phenotypic adaptations, but also the genetics of UV and violet vision.

An epistatic ratchet constrains the direction of glucocorticoid receptor evolution

The extent to which evolution is reversible has long fascinated biologists. Most previous work on the reversibility of morphological and life-history evolution has been indecisive, because of uncertainty and bias in the methods used to infer ancestral states for such characters. Further, despite theoretical work on the factors that could contribute to irreversibility, there is little empirical evidence on its causes, because sufficient understanding of the mechanistic basis for the evolution of new or ancestral phenotypes is seldom available. By studying the reversibility of evolutionary changes in protein structure and function, these limitations can be overcome. Here we show, using the evolution of hormone specificity in the vertebrate glucocorticoid receptor as a case-study, that the evolutionary path by which this protein acquired its new function soon became inaccessible to reverse exploration. Using ancestral gene reconstruction, protein engineering and X-ray crystallography, we demonstrate that five subsequent 'restrictive' mutations, which optimized the new specificity of the glucocorticoid receptor, also destabilized elements of the protein structure that were required to support the ancestral conformation. Unless these ratchet-like epistatic substitutions are restored to their ancestral states, reversing the key function-switching mutations yields a non-functional protein. Reversing the restrictive substitutions first, however, does nothing to enhance the ancestral function. Our findings indicate that even if selection for the ancestral function were imposed, direct reversal would be extremely unlikely, suggesting an important role for historical contingency in protein evolution.

Following evolutionary paths to protein-protein interactions with high affinity and selectivity

How do intricate multi-residue features such as protein-protein interfaces evolve? To address this question, we evolved a new colicin-immunity binding interaction. We started with Im9, which inhibits its cognate DNase ColE9 at 10(-14) M affinity, and evolved it toward ColE7, which it inhibits promiscuously (Kd > 10(-8) M). Iterative rounds of random mutagenesis and selection toward higher affinity for ColE7, and selectivity (against ColE9 inhibition), led to an approximately 10(5)-fold increase in affinity and a 10(8)-fold increase in selectivity. Analysis of intermediates along the evolved variants revealed that changes in the binding configuration of the Im protein uncovered a latent set of interactions, thus providing the key to the rapid divergence of new Im7 variants. Overall, protein-protein interfaces seem to share the evolvability features of enzymes, that is, the exploitation of promiscuous interactions and alternative binding configurations via 'generalist' intermediates, and the key role of compensatory stabilizing mutations in facilitating the divergence of new functions.

Characterization of photoreceptor cell types in the little brown bat Myotis lucifugus (Vespertilionidae)

We report the expression of three visual opsins in the retina of the little brown bat (Myotis lucifugus, Vespertilionidae). Gene sequences for a rod-specific opsin and two cone-specific opsins were cloned from cDNA derived from bat eyes. Comparative sequence analyses indicate that the two cone opsins correspond to an ultraviolet short-wavelength opsin (SWS1) and a long-wavelength opsin (LWS). Immunocytochemistry using antisera to visual opsins revealed that the little brown bat retina contains two types of cone photoreceptors within a rod-dominated background. However, unlike other mammalian photoreceptors, M. lucifugus cones and rods are morphologically indistinguishable by light microscopy. Both photoreceptor types have a thin, elongated outer segment. Using microspectrophotometry we classified the absorption spectrum for the ubiquitous rods. Similar to other mammals, bat rhodopsin has an absorption peak near 500 nm. Although we were unable to confirm a spectral range, cellular and molecular analyses indicate that M. lucifugus expresses two types of cone visual pigments located within the photoreceptor layer. This study provides important insights into the visual capacity of a nocturnal microchiropteran species.

The opsin repertoire of Jenynsia onca: a new perspective on gene duplication and divergence in livebearers

Background

Jenynsia onca, commonly known as the one sided livebearer, is a member of the family Anablepidae. The opsin gene repertoires of J. onca's close relatives, the four-eyed fish (Anableps anableps) and the guppy (Poecilia reticulata), have been characterized and each found to include one unique LWS opsin. Currently, the relationship among LWS paralogs and orthologs in these species are unclear, making it difficult to test the hypotheses that link vision to morphology or life history traits. The phylogenetic signal appears to have been disrupted by gene conversion. Here we have sequenced the opsin genes of J. onca in order to resolve these relationships.

Findings

We identified nine visual opsins; LWS S180r, LWS S180, LWS P180, SWS1, SWS2A, SWS2B, RH1, RH2-1, and RH2-2. Key site analysis revealed only one unique haplotype, RH2-2, although this is unlikely to shift λ_max significantly. LWS P180 was found to be a product of a gene conversion event with LWS S180, followed by convergence to a proline residue at the 180 site.

Conclusion

Jenynsia onca has at least 9 visual opsins: three LWS, one RH1, two RH2, one SWS1 and two SWS2. The presence of LWS P180 moves the location of the LWS P180-S180 tandem duplication event back to the base of the Poeciliidae-Anablepidae clade, expanding the number of species possessing this unusual blue shifted LWS opsin. The presence of the LWS P180 gene also confirms that gene conversion events have homogenized opsin paralogs in fish, just as they have in humans.

A fish eye out of water: ten visual opsins in the four-eyed fish, Anableps anableps

The “four-eyed” fish Anableps anableps has numerous morphological adaptations that enable above and below-water vision. Here, as the first step in our efforts to identify molecular adaptations for aerial and aquatic vision in this species, we describe the A. anableps visual opsin repertoire. We used PCR, cloning, and sequencing to survey cDNA using unique primers designed to amplify eight sequences from five visual opsin gene subfamilies, SWS1, SWS2, RH1, RH2, and LWS. We also used Southern blotting to count opsin loci in genomic DNA digested with EcoR1 and BamH1. Phylogenetic analyses confirmed the identity of all opsin sequences and allowed us to map gene duplication and divergence events onto a tree of teleost fish. Each of the gene-specific primer sets produced an amplicon from cDNA, indicating that A. anableps possessed and expressed at least eight opsin genes. A second PCR-based survey of genomic and cDNA uncovered two additional LWS genes. Thus, A. anableps has at least ten visual opsins and all but one were expressed in the eyes of the single adult surveyed. Among these ten visual opsins, two have key site haplotypes not found in other fish. Of particular interest is the A. anableps-specific opsin in the LWS subfamily, S180γ, with a SHYAA five key site haplotype. Although A. anableps has a visual opsin gene repertoire similar to that found in other fishes in the suborder Cyprinodontoidei, the LWS opsin subfamily has two loci not found in close relatives, including one with a key site haplotype not found in any other fish species. A. anableps opsin sequence data will be used to design in situ probes allowing us to test the hypothesis that opsin gene expression differs in the distinct ventral and dorsal retinas found in this species.

Reliabilities of identifying positive selection by the branch-site and the site-prediction methods

Natural selection operating in protein-coding genes is often studied by examining the ratio (omega) of the rates of nonsynonymous to synonymous nucleotide substitution. The branch-site method (BSM) based on a likelihood ratio test is one of such tests to detect positive selection for a predetermined branch of a phylogenetic tree. However, because the number of nucleotide substitutions involved is often very small, we conducted a computer simulation to examine the reliability of BSM in comparison with the small-sample method (SSM) based on Fisher's exact test. The results indicate that BSM often generates false positives compared with SSM when the number of nucleotide substitutions is approximately 80 or smaller. Because the omega value is also used for predicting positively selected sites, we examined the reliabilities of the site-prediction methods, using nucleotide sequence data for the dim-light and color vision genes in vertebrates. The results showed that the site-prediction methods have a low probability of identifying functional changes of amino acids experimentally determined and often falsely identify other sites where amino acid substitutions are unlikely to be important. This low rate of predictability occurs because most of the current statistical methods are designed to identify codon sites with high omega values, which may not have anything to do with functional changes. The codon sites showing functional changes generally do not show a high omega value. To understand adaptive evolution, some form of experimental confirmation is necessary.

Directed enzyme evolution: climbing fitness peaks one amino acid at a time

Directed evolution can generate a remarkable range of new enzyme properties. Alternate substrate specificities and reaction selectivities are readily accessible in enzymes from families that are naturally functionally diverse. Activities on new substrates can be obtained by improving variants with broadened specificities or by step-wise evolution through a sequence of more and more challenging substrates. Evolution of highly specific enzymes has been demonstrated, even with positive selection alone. It is apparent that many solutions exist for any given problem, and there are often many paths that lead uphill, one step at a time.

Positive Darwinian selection and the birth of an olfactory receptor clade in teleosts

Trace amine-associated receptors (TAARs) in mammals recently have been shown to function as olfactory receptors. We have delineated the taar gene family in jawless, cartilaginous, and bony fish (zero, 2, and >100 genes, respectively). We conclude that taar genes are evolutionary much younger than the related OR and ORA/V1R olfactory receptor families, which are present already in lamprey, a jawless vertebrate. The 2 cartilaginous fish genes appear to be ancestral for 2 taar classes, each with mammalian and bony fish (teleost) representatives. Unexpectedly, a whole new clade, class III, of taar genes originated even later, within the teleost lineage. Taar genes from all 3 classes are expressed in subsets of zebrafish olfactory receptor neurons, supporting their function as olfactory receptors. The highly conserved TAAR1 (shark, mammalian, and teleost orthologs) is not expressed in the olfactory epithelium and may constitute the sole remnant of a primordial, nonolfactory function of this family. Class III comprises three-fourths of all teleost taar genes and is characterized by the complete loss of the aminergic ligand-binding motif, stringently conserved in the other 2 classes. Two independent intron gains in class III taar genes represent extraordinary evolutionary dynamics, considering the virtual absence of intron gains during vertebrate evolution. The d(N)/d(S) analysis suggests both minimal global negative selection and an unparalleled degree of local positive selection as another hallmark of class III genes. The accelerated evolution of class III teleost taar genes conceivably might mark the birth of another olfactory receptor gene family.

The green-absorbing Drosophila Rh6 visual pigment contains a blue-shifting amino acid substitution that is conserved in vertebrates

The molecular mechanisms that regulate invertebrate visual pigment absorption are poorly understood. Through sequence analysis and functional investigation of vertebrate visual pigments, numerous amino acid substitutions important for this adaptive process have been identified. Here we describe a serine/alanine (S/A) substitution in long wavelength-absorbing Drosophila visual pigments that occurs at a site corresponding to Ala-292 in bovine rhodopsin. This S/A substitution accounts for a 10-17-nm absorption shift in visual pigments of this class. Additionally, we demonstrate that substitution of a cysteine at the same site, as occurs in the blue-absorbing Rh5 pigment, accounts for a 4-nm shift. Substitutions at this site are the first spectrally significant amino acid changes to be identified for invertebrate pigments sensitive to visible light and are the first evidence of a conserved tuning mechanism in vertebrate and invertebrate pigments of this class.

The evolution of animal chemosensory receptor gene repertoires: roles of chance and necessity

Chemosensory receptors are essential for the survival of organisms that range from bacteria to mammals. Recent studies have shown that the numbers of functional chemosensory receptor genes and pseudogenes vary enormously among the genomes of different animal species. Although much of the variation can be explained by the adaptation of organisms to different environments, it has become clear that a substantial portion is generated by genomic drift, a random process of gene duplication and deletion. Genomic drift also generates a substantial amount of copy-number variation in chemosensory receptor genes within species. It seems that mutation by gene duplication and inactivation has important roles in both the adaptive and non-adaptive evolution of chemosensation.

Detecting clusters of mutations

Positive selection for protein function can lead to multiple mutations within a small stretch of DNA, i.e., to a cluster of mutations. Recently, Wagner proposed a method to detect such mutation clusters. His method, however, did not take into account that residues with high solvent accessibility are inherently more variable than residues with low solvent accessibility. Here, we propose a new algorithm to detect clustered evolution. Our algorithm controls for different substitution probabilities at buried and exposed sites in the tertiary protein structure, and uses random permutations to calculate accurate P values for inferred clusters. We apply the algorithm to genomes of bacteria, fly, and mammals, and find several clusters of mutations in functionally important regions of proteins. Surprisingly, clustered evolution is a relatively rare phenomenon. Only between 2% and 10% of the genes we analyze contain a statistically significant mutation cluster. We also find that not controlling for solvent accessibility leads to an excess of clusters in terminal and solvent-exposed regions of proteins. Our algorithm provides a novel method to identify functionally relevant divergence between groups of species. Moreover, it could also be useful to detect artifacts in automatically assembled genomes.