The tuatara genome reveals ancient features of amniote evolution.

The tuatara (Sphenodon punctatus)-the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana1,2-is an iconic species that is endemic to New Zealand2,3. A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes2,4. Here we analyse the genome of the tuatara, which-at approximately 5 Gb-is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing.

Publisher Correction: The tuatara genome reveals ancient features of amniote evolution.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Diversity and Molecular Evolution of Nonvisual Opsin Genes across Environmental, Developmental, and Morphological Adaptations in Frogs.

Nonvisual opsins are transmembrane proteins expressed in the eyes and other tissues of many animals. When paired with a light-sensitive chromophore, nonvisual opsins form photopigments involved in various nonvisual, light-detection functions including circadian rhythm regulation, light-seeking behaviors, and seasonal responses. Here, we investigate the molecular evolution of nonvisual opsin genes in anuran amphibians (frogs and toads). We test several evolutionary hypotheses including the predicted loss of nonvisual opsins due to nocturnal ancestry and potential functional differences in nonvisual opsins resulting from environmental light variation across diverse anuran ecologies. Using whole-eye transcriptomes of 81 species, combined with genomes, multitissue transcriptomes, and independently annotated genes from an additional 21 species, we identify which nonvisual opsins are present in anuran genomes and those that are also expressed in the eyes, compare selective constraint among genes, and test for potential adaptive evolution by comparing selection between discrete ecological classes. At the genomic level, we recovered all 18 ancestral vertebrate nonvisual opsins, indicating that anurans demonstrate the lowest documented amount of opsin gene loss among ancestrally nocturnal tetrapods. We consistently found expression of 14 nonvisual opsins in anuran eyes and detected positive selection in a subset of these genes. We also found shifts in selective constraint acting on nonvisual opsins in frogs with differing activity periods, habitats, distributions, life histories, and pupil shapes, which may reflect functional adaptation. Although many nonvisual opsins remain poorly understood, these findings provide insight into the diversity and evolution of these genes across anurans, filling an important gap in our understanding of vertebrate opsins and setting the stage for future research on their functional evolution across taxa.

Diversity and Evolution of Frog Visual Opsins: Spectral Tuning and Adaptation to Distinct Light Environments.

Visual systems adapt to different light environments through several avenues including optical changes to the eye and neurological changes in how light signals are processed and interpreted. Spectral sensitivity can evolve via changes to visual pigments housed in the retinal photoreceptors through gene duplication and loss, differential and coexpression, and sequence evolution. Frogs provide an excellent, yet understudied, system for visual evolution research due to their diversity of ecologies (including biphasic aquatic-terrestrial life cycles) that we hypothesize imposed different selective pressures leading to adaptive evolution of the visual system, notably the opsins that encode the protein component of the visual pigments responsible for the first step in visual perception. Here, we analyze the diversity and evolution of visual opsin genes from 93 new eye transcriptomes plus published data for a combined dataset spanning 122 frog species and 34 families. We find that most species express the four visual opsins previously identified in frogs but show evidence for gene loss in two lineages. Further, we present evidence of positive selection in three opsins and shifts in selective pressures associated with differences in habitat and life history, but not activity pattern. We identify substantial novel variation in the visual opsins and, using microspectrophotometry, find highly variable spectral sensitivities, expanding known ranges for all frog visual pigments. Mutations at spectral-tuning sites only partially account for this variation, suggesting that frogs have used tuning pathways that are unique among vertebrates. These results support the hypothesis of adaptive evolution in photoreceptor physiology across the frog tree of life in response to varying environmental and ecological factors and further our growing understanding of vertebrate visual evolution.

Scaling up Functional Analyses of the G Protein-Coupled Receptor Rhodopsin.

Eukaryotic cells use G protein-coupled receptors (GPCRs) to convert external stimuli into internal signals to elicit cellular responses. However, how mutations in GPCR-coding genes affect GPCR activation and downstream signaling pathways remain poorly understood. Approaches such as deep mutational scanning show promise in investigations of GPCRs, but a high-throughput method to measure rhodopsin activation has yet to be achieved. Here, we scale up a fluorescent reporter assay in budding yeast that we engineered to study rhodopsin's light-activated signal transduction. Using this approach, we measured the mutational effects of over 1200 individual human rhodopsin mutants, generated by low-frequency random mutagenesis of the GPCR rhodopsin (RHO) gene. Analysis of the data in the context of rhodopsin's three-dimensional structure reveals that transmembrane helices are generally less tolerant to mutations compared to flanking helices that face the lipid bilayer, which suggest that mutational tolerance is contingent on both the local environment surrounding specific residues and the specific position of these residues in the protein structure. Comparison of functional scores from our screen to clinically identified rhodopsin disease variants found many pathogenic mutants to be loss of function. Lastly, functional scores from our assay were consistent with a complex counterion mechanism involved in ligand-binding and rhodopsin activation. Our results demonstrate that deep mutational scanning is possible for rhodopsin activation and can be an effective method for revealing properties of mutational tolerance that may be generalizable to other transmembrane proteins.

Transcriptomic evidence for visual adaptation during the aquatic to terrestrial metamorphosis in leopard frogs.

BACKGROUND: Differences in morphology, ecology, and behavior through ontogeny can result in opposing selective pressures at different life stages. Most animals, however, transition through two or more distinct phenotypic phases, which is hypothesized to allow each life stage to adapt more freely to its ecological niche. How this applies to sensory systems, and in particular how sensory systems adapt across life stages at the molecular level, is not well understood. Here, we used whole-eye transcriptomes to investigate differences in gene expression between tadpole and juvenile southern leopard frogs (Lithobates sphenocephalus), which rely on vision in aquatic and terrestrial light environments, respectively. Because visual physiology changes with light levels, we also tested the effect of light and dark exposure. RESULTS: We found 42% of genes were differentially expressed in the eyes of tadpoles versus juveniles and 5% for light/dark exposure. Analyses targeting a curated subset of visual genes revealed significant differential expression of genes that control aspects of visual function and development, including spectral sensitivity and lens composition. Finally, microspectrophotometry of photoreceptors confirmed shifts in spectral sensitivity predicted by the expression results, consistent with adaptation to distinct light environments. CONCLUSIONS: Overall, we identified extensive expression-level differences in the eyes of tadpoles and juveniles related to observed morphological and physiological changes through metamorphosis and corresponding adaptive shifts to improve vision in the distinct aquatic and terrestrial light environments these frogs inhabit during their life cycle. More broadly, these results suggest that decoupling of gene expression can mediate the opposing selection pressures experienced by organisms with complex life cycles that inhabit different environmental conditions throughout ontogeny.

Ocular lens morphology is influenced by ecology and metamorphosis in frogs and toads.

The shape and relative size of an ocular lens affect the focal length of the eye, with consequences for visual acuity and sensitivity. Lenses are typically spherical in aquatic animals with camera-type eyes and axially flattened in terrestrial species to facilitate vision in optical media with different refractive indices. Frogs and toads (Amphibia: Anura) are ecologically diverse, with many species shifting from aquatic to terrestrial ecologies during metamorphosis. We quantified lens shape and relative size using 179 micro X-ray computed tomography scans of 126 biphasic anuran species and tested for correlations with life stage, environmental transitions, adult habits and adult activity patterns. Across broad phylogenetic diversity, tadpole lenses are more spherical than those of adults. Biphasic species with aquatic larvae and terrestrial adults typically undergo ontogenetic changes in lens shape, whereas species that remain aquatic as adults tend to retain more spherical lenses after metamorphosis. Further, adult lens shape is influenced by adult habit; notably, fossorial adults tend to retain spherical lenses following metamorphosis. Finally, lens size relative to eye size is smaller in aquatic and semiaquatic species than other adult ecologies. Our study demonstrates how ecology shapes visual systems, and the power of non-invasive imaging of museum specimens for studying sensory evolution.

Evolution, inactivation and loss of short wavelength-sensitive opsin genes during the diversification of Neotropical cichlids.

Natural variation in the number, expression and function of sensory genes in an organism's genome is often tightly linked to different ecological and evolutionary forces. Opsin genes, which code for the first step in visual transduction, are ideal models for testing how ecological factors such as light environment may influence visual system adaptation. Neotropical cichlid fishes are a highly ecologically diverse group that evolved in a variety of aquatic habitats, including black (stained), white (opaque) and clear waters. We used cross-species exon capture to sequence Neotropical cichlid short wavelength-sensitive (SWS) opsins, which mediate ultraviolet (UV) to blue visual sensitivity. Neotropical cichlid SWS1 opsin (UV-sensitive) underwent a relaxation of selective constraint during the early phases of cichlid diversification in South America, leading to pseudogenization and loss. Conversely, SWS2a (blue-sensitive) experienced a burst of episodic positive selection at the base of the South American cichlid radiation. This burst coincides with SWS1 relaxation and loss, and is consistent with findings in ecomorphological studies characterizing a period of extensive ecological divergence in Neotropical cichlids. We use ancestral sequence reconstruction and protein modelling to investigate mutations along this ancestral branch that probably modified SWS2a function. Together, our results suggest that variable light environments played a prominent early role in shaping SWS opsin diversity during the Neotropical cichlid radiation. Our results also illustrate that long-term evolution under light-limited conditions in South America may have reduced visual system plasticity; specifically, early losses of UV sensitivity may have constrained the evolutionary trajectory of Neotropical cichlid vision.

Eye size and investment in frogs and toads correlate with adult habitat, activity pattern and breeding ecology.

Frogs and toads (Amphibia: Anura) display diverse ecologies and behaviours, which are often correlated with visual capacity in other vertebrates. Additionally, anurans exhibit a broad range of relative eye sizes, which have not previously been linked to ecological factors in this group. We measured relative investment in eye size and corneal size for 220 species of anurans representing all 55 currently recognized families and tested whether they were correlated with six natural history traits hypothesized to be associated with the evolution of eye size. Anuran eye size was significantly correlated with habitat, with notable decreases in eye investment among fossorial, subfossorial and aquatic species. Relative eye size was also associated with mating habitat and activity pattern. Compared to other vertebrates, anurans have relatively large eyes for their body size, indicating that vision is probably of high importance. Our study reveals the role that ecology and behaviour may have played in the evolution of anuran visual systems and highlights the usefulness of museum specimens, and importance of broad taxonomic sampling, for interpreting macroecological patterns.

Evolution of nonspectral rhodopsin function at high altitudes.

High-altitude environments present a range of biochemical and physiological challenges for organisms through decreases in oxygen, pressure, and temperature relative to lowland habitats. Protein-level adaptations to hypoxic high-altitude conditions have been identified in multiple terrestrial endotherms; however, comparable adaptations in aquatic ectotherms, such as fishes, have not been as extensively characterized. In enzyme proteins, cold adaptation is attained through functional trade-offs between stability and activity, often mediated by substitutions outside the active site. Little is known whether signaling proteins [e.g., G protein-coupled receptors (GPCRs)] exhibit natural variation in response to cold temperatures. Rhodopsin (RH1), the temperature-sensitive visual pigment mediating dim-light vision, offers an opportunity to enhance our understanding of thermal adaptation in a model GPCR. Here, we investigate the evolution of rhodopsin function in an Andean mountain catfish system spanning a range of elevations. Using molecular evolutionary analyses and site-directed mutagenesis experiments, we provide evidence for cold adaptation in RH1. We find that unique amino acid substitutions occur at sites under positive selection in high-altitude catfishes, located at opposite ends of the RH1 intramolecular hydrogen-bonding network. Natural high-altitude variants introduced into these sites via mutagenesis have limited effects on spectral tuning, yet decrease the stability of dark-state and light-activated rhodopsin, accelerating the decay of ligand-bound forms. As found in cold-adapted enzymes, this phenotype likely compensates for a cold-induced decrease in kinetic rates-properties of rhodopsin that mediate rod sensitivity and visual performance. Our results support a role for natural variation in enhancing the performance of GPCRs in response to cold temperatures.

Shifts in Selective Pressures on Snake Phototransduction Genes Associated with Photoreceptor Transmutation and Dim-Light Ancestry.

The visual systems of snakes are heavily modified relative to other squamates, a condition often thought to reflect their fossorial origins. Further modifications are seen in caenophidian snakes, where evolutionary transitions between rod and cone photoreceptors, termed photoreceptor transmutations, have occurred in many lineages. Little previous work, however, has focused on the molecular evolutionary underpinnings of these morphological changes. To address this, we sequenced seven snake eye transcriptomes and utilized new whole-genome and targeted capture sequencing data. We used these data to analyze gene loss and shifts in selection pressures in phototransduction genes that may be associated with snake evolutionary origins and photoreceptor transmutation. We identified the surprising loss of rhodopsin kinase (GRK1), despite a low degree of gene loss overall and a lack of relaxed selection early during snake evolution. These results provide some of the first evolutionary genomic corroboration for a dim-light ancestor that lacks strong fossorial adaptations. Our results also indicate that snakes with photoreceptor transmutation experienced significantly different selection pressures from other reptiles. Significant positive selection was found primarily in cone-specific genes, but not rod-specific genes, contrary to our expectations. These results reveal potential molecular adaptations associated with photoreceptor transmutation and also highlight unappreciated functional differences between rod- and cone-specific phototransduction proteins. This intriguing example of snake visual system evolution illustrates how the underlying molecular components of a complex system can be reshaped in response to changing selection pressures.

Functional Shifts in Bat Dim-Light Visual Pigment Are Associated with Differing Echolocation Abilities and Reveal Molecular Adaptation to Photic-Limited Environments.

Bats are excellent models for studying the molecular basis of sensory adaptation. In Chiroptera, a sensory trade-off has been proposed between the visual and auditory systems, though the extent of this association has yet to be fully examined. To investigate whether variation in visual performance is associated with echolocation, we experimentally assayed the dim-light visual pigment rhodopsin from bat species with differing echolocation abilities. While spectral tuning properties were similar among bats, we found that the rate of decay of their light-activated state was significantly slower in a nonecholocating bat relative to species that use distinct echolocation strategies, consistent with a sensory trade-off hypothesis. We also found that these rates of decay were remarkably slower compared with those of other mammals, likely indicating an adaptation to dim light. To examine whether functional changes in rhodopsin are associated with shifts in selection intensity upon bat Rh1 sequences, we implemented selection analyses using codon-based likelihood clade models. While no shifts in selection were identified in response to diverse echolocation abilities of bats, we detected a significant increase in the intensity of evolutionary constraint accompanying the diversification of Chiroptera. Taken together, this suggests that substitutions that modulate the stability of the light-activated rhodopsin state were likely maintained through intensified constraint after bats diversified, being finely tuned in response to novel sensory specializations. Our study demonstrates the power of combining experimental and computational approaches for investigating functional mechanisms underlying the evolution of complex sensory adaptations.

Evolutionary transformation of rod photoreceptors in the all-cone retina of a diurnal garter snake.

Vertebrate retinas are generally composed of rod (dim-light) and cone (bright-light) photoreceptors with distinct morphologies that evolved as adaptations to nocturnal/crepuscular and diurnal light environments. Over 70 years ago, the "transmutation" theory was proposed to explain some of the rare exceptions in which a photoreceptor type is missing, suggesting that photoreceptors could evolutionarily transition between cell types. Although studies have shown support for this theory in nocturnal geckos, the origins of all-cone retinas, such as those found in diurnal colubrid snakes, remain a mystery. Here we investigate the evolutionary fate of the rods in a diurnal garter snake and test two competing hypotheses: (i) that the rods, and their corresponding molecular machinery, were lost or (ii) that the rods were evolutionarily modified to resemble, and function, as cones. Using multiple approaches, we find evidence for a functional and unusually blue-shifted rhodopsin that is expressed in small single "cones." Moreover, these cones express rod transducin and have rod ultrastructural features, providing strong support for the hypothesis that they are not true cones, as previously thought, but rather are modified rods. Several intriguing features of garter snake rhodopsin are suggestive of a more cone-like function. We propose that these cone-like rods may have evolved to regain spectral sensitivity and chromatic discrimination as a result of ancestral losses of middle-wavelength cone opsins in early snake evolution. This study illustrates how sensory evolution can be shaped not only by environmental constraints but also by historical contingency in forming new cell types with convergent functionality.

Accelerated Evolution and Functional Divergence of the Dim Light Visual Pigment Accompanies Cichlid Colonization of Central America.

Cichlids encompass one of the most diverse groups of fishes in South and Central America, and show extensive variation in life history, morphology, and colouration. While studies of visual system evolution in cichlids have focussed largely on the African rift lake species flocks, Neotropical cichlids offer a unique opportunity to investigate visual system evolution at broader temporal and geographic scales. South American cichlid colonization of Central America has likely promoted accelerated rates of morphological evolution in Central American lineages as they encountered reduced competition, renewed ecological opportunity, and novel aquatic habitats. To investigate whether such transitions have influenced molecular evolution of vision in Central American cichlids, we sequenced the dim-light rhodopsin gene in 101 Neotropical cichlid species, spanning the diversity of the clade. We find strong evidence for increased rates of evolution in Central American cichlid rhodopsin relative to South American lineages, and identify several sites under positive selection in rhodopsin that likely contribute to adaptation to different photic environments. We expressed a Neotropical cichlid rhodopsin protein invitro for the first time, and found that while its spectral tuning properties were characteristic of typical vertebrate rhodopsin pigments, the rate of decay of its active signalling form was much slower, consistent with dim light adaptation in other vertebrate rhodopsins. Using site-directed mutagenesis combined with spectroscopic assays, we found that a key amino acid substitution present in some Central American cichlids accelerates the rate of decay of active rhodopsin, which may mediate adaptation to clear water habitats.

The role of ecological factors in shaping bat cone opsin evolution.

Bats represent one of the largest and most striking nocturnal mammalian radiations, exhibiting many visual system specializations for performance in light-limited environments. Despite representing the greatest ecological diversity and species richness in Chiroptera, Neotropical lineages have been undersampled in molecular studies, limiting the potential for identifying signatures of selection on visual genes associated with differences in bat ecology. Here, we investigated how diverse ecological pressures mediate long-term shifts in selection upon long-wavelength (Lws) and short-wavelength (Sws1) opsins, photosensitive cone pigments that form the basis of colour vision in most mammals, including bats. We used codon-based likelihood clade models to test whether ecological variables associated with reliance on visual information (e.g. echolocation ability and diet) or exposure to varying light environments (e.g. roosting behaviour and foraging habitat) mediated shifts in evolutionary rates in bat cone opsin genes. Using additional cone opsin sequences from newly sequenced eye transcriptomes of six Neotropical bat species, we found significant evidence for different ecological pressures influencing the evolution of the cone opsins. While Lws is evolving under significantly lower constraint in highly specialized high-duty cycle echolocating lineages, which have enhanced sonar ability to detect and track targets, variation in Sws1 constraint was significantly associated with foraging habitat, exhibiting elevated rates of evolution in species that forage among vegetation. This suggests that increased reliance on echolocation as well as the spectral environment experienced by foraging bats may differentially influence the evolution of different cone opsins. Our study demonstrates that different ecological variables may underlie contrasting evolutionary patterns in bat visual opsins, and highlights the suitability of clade models for testing ecological hypotheses of visual evolution.

Cone-like rhodopsin expressed in the all-cone retina of the colubrid pine snake as a potential adaptation to diurnality.

Colubridae is the largest and most diverse family of snakes, with visual systems that reflect this diversity, encompassing a variety of retinal photoreceptor organizations. The transmutation theory proposed by Walls postulates that photoreceptors could evolutionarily transition between cell types in squamates, but few studies have tested this theory. Recently, evidence for transmutation and rod-like machinery in an all-cone retina has been identified in a diurnal garter snake (Thamnophis), and it appears that the rhodopsin gene at least may be widespread among colubrid snakes. However, functional evidence supporting transmutation beyond the existence of the rhodopsin gene remains rare. We examined the all-cone retina of another colubrid, Pituophis melanoleucus, thought to be more secretive/burrowing than Thamnophis We found that P. melanoleucus expresses two cone opsins (SWS1, LWS) and rhodopsin (RH1) within the eye. Immunohistochemistry localized rhodopsin to the outer segment of photoreceptors in the all-cone retina of the snake and all opsin genes produced functional visual pigments when expressed in vitro Consistent with other studies, we found that P. melanoleucus rhodopsin is extremely blue-shifted. Surprisingly, P. melanoleucus rhodopsin reacted with hydroxylamine, a typical cone opsin characteristic. These results support the idea that the rhodopsin-containing photoreceptors of P. melanoleucus are the products of evolutionary transmutation from rod ancestors, and suggest that this phenomenon may be widespread in colubrid snakes. We hypothesize that transmutation may be an adaptation for diurnal, brighter-light vision, which could result in increased spectral sensitivity and chromatic discrimination with the potential for colour vision.

A second visual rhodopsin gene, rh1-2, is expressed in zebrafish photoreceptors and found in other ray-finned fishes.

Rhodopsin (rh1) is the visual pigment expressed in rod photoreceptors of vertebrates that is responsible for initiating the critical first step of dim-light vision. Rhodopsin is usually a single copy gene; however, we previously discovered a novel rhodopsin-like gene expressed in the zebrafish retina, rh1-2, which we identified as a functional photosensitive pigment that binds 11-cis retinal and activates in response to light. Here, we localized expression of rh1-2 in the zebrafish retina to a subset of peripheral photoreceptor cells, which indicates a partially overlapping expression pattern with rh1 We also expressed, purified and characterized Rh1-2, including investigation of the stability of the biologically active intermediate. Using fluorescence spectroscopy, we found the half-life of the rate of retinal release of Rh1-2 following photoactivation to be more similar to that of the visual pigment rhodopsin than to the non-visual pigment exo-rhodopsin (exorh), which releases retinal around 5 times faster. Phylogenetic and molecular evolutionary analyses show that rh1-2 has ancient origins within teleost fishes, is under similar selective pressure to rh1, and likely experienced a burst of positive selection following its duplication and divergence from rh1 These findings indicate that rh1-2 is another functional visual rhodopsin gene, which contradicts the prevailing notion that visual rhodopsin is primarily found as a single copy gene within ray-finned fishes. The reasons for retention of this duplicate gene, as well as possible functional consequences for the visual system, are discussed.

Comparative sequence analyses of rhodopsin and RPE65 reveal patterns of selective constraint across hereditary retinal disease mutations.

Retinitis pigmentosa (RP) comprises several heritable diseases that involve photoreceptor, and ultimately retinal, degeneration. Currently, mutations in over 50 genes have known links to RP. Despite advances in clinical characterization, molecular characterization of RP remains challenging due to the heterogeneous nature of causal genes, mutations, and clinical phenotypes. In this study, we compiled large datasets of two important visual genes associated with RP: rhodopsin, which initiates the phototransduction cascade, and the retinoid isomerase RPE65, which regenerates the visual cycle. We used a comparative evolutionary approach to investigate the relationship between interspecific sequence variation and pathogenic mutations that lead to degenerative retinal disease. Using codon-based likelihood methods, we estimated evolutionary rates (d N/d S) across both genes in a phylogenetic context to investigate differences between pathogenic and nonpathogenic amino acid sites. In both genes, disease-associated sites showed significantly lower evolutionary rates compared to nondisease sites, and were more likely to occur in functionally critical areas of the proteins. The nature of the dataset (e.g., vertebrate or mammalian sequences), as well as selection of pathogenic sites, affected the differences observed between pathogenic and nonpathogenic sites. Our results illustrate that these methods can serve as an intermediate step in understanding protein structure and function in a clinical context, particularly in predicting the relative pathogenicity (i.e., functional impact) of point mutations and their downstream phenotypic effects. Extensions of this approach may also contribute to current methods for predicting the deleterious effects of candidate mutations and to the identification of protein regions under strong constraint where we expect pathogenic mutations to occur.

The Burmese python genome reveals the molecular basis for extreme adaptation in snakes.

Snakes possess many extreme morphological and physiological adaptations. Identification of the molecular basis of these traits can provide novel understanding for vertebrate biology and medicine. Here, we study snake biology using the genome sequence of the Burmese python (Python molurus bivittatus), a model of extreme physiological and metabolic adaptation. We compare the python and king cobra genomes along with genomic samples from other snakes and perform transcriptome analysis to gain insights into the extreme phenotypes of the python. We discovered rapid and massive transcriptional responses in multiple organ systems that occur on feeding and coordinate major changes in organ size and function. Intriguingly, the homologs of these genes in humans are associated with metabolism, development, and pathology. We also found that many snake metabolic genes have undergone positive selection, which together with the rapid evolution of mitochondrial proteins, provides evidence for extensive adaptive redesign of snake metabolic pathways. Additional evidence for molecular adaptation and gene family expansions and contractions is associated with major physiological and phenotypic adaptations in snakes; genes involved are related to cell cycle, development, lungs, eyes, heart, intestine, and skeletal structure, including GRB2-associated binding protein 1, SSH, WNT16, and bone morphogenetic protein 7. Finally, changes in repetitive DNA content, guanine-cytosine isochore structure, and nucleotide substitution rates indicate major shifts in the structure and evolution of snake genomes compared with other amniotes. Phenotypic and physiological novelty in snakes seems to be driven by system-wide coordination of protein adaptation, gene expression, and changes in the structure of the genome.

Divergent positive selection in rhodopsin from lake and riverine cichlid fishes.

Studies of cichlid evolution have highlighted the importance of visual pigment genes in the spectacular radiation of the African rift lake cichlids. Recent work, however, has also provided strong evidence for adaptive diversification of riverine cichlids in the Neotropics, which inhabit environments of markedly different spectral properties from the African rift lakes. These ecological and/or biogeographic differences may have imposed divergent selective pressures on the evolution of the cichlid visual system. To test these hypotheses, we investigated the molecular evolution of the dim-light visual pigment, rhodopsin. We sequenced rhodopsin from Neotropical and African riverine cichlids and combined these data with published sequences from African cichlids. We found significant evidence for positive selection using random sites codon models in all cichlid groups, with the highest levels in African lake cichlids. Tests using branch-site and clade models that partitioned the data along ecological (lake, river) and/or biogeographic (African, Neotropical) boundaries found significant evidence of divergent selective pressures among cichlid groups. However, statistical comparisons among these models suggest that ecological, rather than biogeographic, factors may be responsible for divergent selective pressures that have shaped the evolution of the visual system in cichlids. We found that branch-site models did not perform as well as clade models for our data set, in which there was evidence for positive selection in the background. One of our most intriguing results is that the amino acid sites found to be under positive selection in Neotropical and African lake cichlids were largely nonoverlapping, despite falling into the same three functional categories: spectral tuning, retinal uptake/release, and rhodopsin dimerization. Taken together, these results would imply divergent selection across cichlid clades, but targeting similar functions. This study highlights the importance of molecular investigations of ecologically important groups and the flexibility of clade models in explicitly testing ecological hypotheses.