Search

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., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

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., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

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., (© 2024. The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature.)

Animals with enhanced dim-light sensitivity are at higher risk of light-induced retinal degeneration when exposed to bright light conditions. 1 , 2 , 3 , 4 This trade-off is mediated by the rod photoreceptor sensory protein, rhodopsin (RHO), and its toxic vitamin A chromophore by-product, all-trans retinal. 5 , 6 , 7 , 8 Rod arrestin (Arr-1) binds to RHO and promotes sequestration of excess all-trans retinal, 9 , 10 which has recently been suggested as a protective mechanism against photoreceptor cell death. 2 , 11 We investigated Arr-1 evolution in animals at high risk of retinal damage due to periodic bright-light exposure of rod-dominated retinas. Here, we find the convergent evolution of enhanced Arr-1/RHO all-trans-retinal sequestration in owls and deep-diving whales. Statistical analyses reveal a parallel acceleration of Arr-1 evolutionary rates in these lineages, which is associated with the introduction of a rare Arr-1 mutation (Q69R) into the RHO-Arr-1 binding interface. Using in vitro assays, we find that this single mutation significantly enhances RHO-all-trans-retinal sequestration by ∼30%. This functional convergence across 300 million years of evolutionary divergence suggests that Arr-1 and RHO may play an underappreciated role in the photoprotection of the eye, with potentially vast clinical significance., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 Elsevier Inc. All rights reserved.)

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.

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., (© 2022. The Author(s).)

Among major vertebrate groups, anurans (frogs and toads) are understudied with regard to their visual systems, and little is known about variation among species that differ in ecology. We sampled North American anurans representing diverse evolutionary and life histories that likely possess visual systems adapted to meet different ecological needs. Using standard molecular techniques, visual opsin genes, which encode the protein component of visual pigments, were obtained from anuran retinas. Additionally, we extracted the visual opsins from publicly available genome and transcriptome assemblies, further increasing the phylogenetic and ecological diversity of our dataset to 33 species in total. We found that anurans consistently express four visual opsin genes ( RH1 , LWS , SWS1 , and SWS2 , but not RH2 ) even though reported photoreceptor complements vary widely among species. The proteins encoded by these genes showed considerable sequence variation among species, including at sites known to shift the spectral sensitivity of visual pigments in other vertebrates and had conserved substitutions that may be related to dim-light adaptation. Using molecular evolutionary analyses of selection (d N /d S ) we found significant evidence for positive selection at a subset of sites in the dim-light rod opsin gene RH1 and the long wavelength sensitive cone opsin LWS . The function of sites inferred to be under positive selection are largely unknown, but a few are likely to affect spectral sensitivity and other visual pigment functions based on proximity to previously identified sites in other vertebrates. We also found the first evidence of visual opsin duplication in an amphibian with the duplication of the LWS gene in the African bullfrog, which had distinct LWS copies on the sex chromosomes suggesting the possibility of sex-specific visual adaptation. Taken together, our results indicate that ecological factors, such as habitat and life history, as well as behavior, may be driving changes to anuran visual systems., Competing Interests: The authors declare no conflict of interest., (© 2022 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

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., (© 2021 John Wiley & Sons Ltd.)

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.

The tuatara (Sphenodon punctatus)-the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana 1,2 -is an iconic species that is endemic to New Zealand 2,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 amniotes 2,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.

Most vertebrates use a combination of rod and cone photoreceptors to enable vision in conditions ranging from starlight to direct sunlight. Nocturnal geckos, however, have simplex retinas that contain only rods in terms of morphology and physiology, but these rods are thought to be derived from cones through an evolutionary process known as photoreceptor transmutation. To investigate this, we generated eye transcriptomes and analyzed patterns of phototransduction gene evolution in geckos in comparison to other reptiles. We confirm that geckos have lost several major components of the rod phototransduction pathway, including rod opsin (RH1), which we identified as a pseudogene in multiple genomes. We also identified a partial rod transducin transcript, but found no evidence of the protein in retinal sections. However, we find that geckos express several complete rod phototransduction transcripts in the eye, which may contribute to the rod-like physiology of nocturnal gecko photoreceptors. Finally, we found surprising evidence that even though photoreceptor transmutation evolved independently in geckos and snakes, they have experienced parallel shifts in selective constraint on phototransduction genes. These results implicate adaptive change in the underlying molecular machinery of visual transduction, in addition to the convergent changes in cellular morphology, during photoreceptor transmutation., (© 2019 The Author(s). Evolution © 2019 The Society for the Study of Evolution.)

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.

Colubridae represents the most phenotypically diverse and speciose family of snakes, yet no well-assembled and annotated genome exists for this lineage. Here, we report and analyze the genome of the garter snake, Thamnophis sirtalis, a colubrid snake that is an important model species for research in evolutionary biology, physiology, genomics, behavior, and the evolution of toxin resistance. Using the garter snake genome, we show how snakes have evolved numerous adaptations for sensing and securing prey, and identify features of snake genome structure that provide insight into the evolution of amniote genomes. Analyses of the garter snake and other squamate reptile genomes highlight shifts in repeat element abundance and expansion within snakes, uncover evidence of genes under positive selection, and provide revised neutral substitution rate estimates for squamates. Our identification of Z and W sex chromosome-specific scaffolds provides evidence for multiple origins of sex chromosome systems in snakes and demonstrates the value of this genome for studying sex chromosome evolution. Analysis of gene duplication and loss in visual and olfactory gene families supports a dim-light ancestral condition in snakes and indicates that olfactory receptor repertoires underwent an expansion early in snake evolution. Additionally, we provide some of the first links between secreted venom proteins, the genes that encode them, and their evolutionary origins in a rear-fanged colubrid snake, together with new genomic insight into the coevolutionary arms race between garter snakes and highly toxic newt prey that led to toxin resistance in garter snakes.

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.

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., (© 2018 The Author(s).)

Convergent evolution in response to similar selective pressures is a well-known phenomenon in evolutionary biology. Less well understood is how selection drives convergence in protein function, and the underlying mechanisms by which this can be achieved. Here, we investigate functional convergence in the visual system of two distantly related lineages of high-altitude adapted Andean and Himalayan catfishes. Statistical analyses revealed in the two high-altitude lineages, a parallel acceleration of evolutionary rates in rhodopsin, the dim-light visual pigment. However, the elevated rates were found to be accompanied by substitutions at different sites in the protein. Experiments substituting Andean- or Himalayan-specific residues significantly accelerated the kinetic rates of rhodopsin, destabilizing the ligand-bound forms. As found in cold-adapted enzymes, this phenotype likely compensates for a cold-induced decrease in kinetic rates, properties of rhodopsin mediating rod sensitivity and visual performance. Our study suggests that molecular convergence in protein function can be driven by parallel shifts in evolutionary rates but via nonparallel molecular mechanisms. Signatures of natural selection may therefore be a powerful guide for identifying complex instances of functional convergence across a wider range of protein systems., (© 2017 The Author(s). Evolution © 2017 The Society for the Study of Evolution.)

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 Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

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., Competing Interests: The authors declare no conflict of interest.

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., Competing Interests: Competing interestsThe authors declare no competing or financial interests., (© 2017. Published by The Company of Biologists Ltd.)

Despite continued advances in sequencing technologies, there is a need for methods that can efficiently sequence large numbers of genes from diverse species. One approach to accomplish this is targeted capture (hybrid enrichment). While these methods are well established for genome resequencing projects, cross-species capture strategies are still being developed and generally focus on the capture of conserved regions, rather than complete coding regions from specific genes of interest. The resulting data is thus useful for phylogenetic studies, but the wealth of comparative data that could be used for evolutionary and functional studies is lost. Here, we design and implement a targeted capture method that enables recovery of complete coding regions across broad taxonomic scales. Capture probes were designed from multiple reference species and extensively tiled in order to facilitate cross-species capture. Using novel bioinformatics pipelines we were able to recover nearly all of the targeted genes with high completeness from species that were up to 200 myr divergent. Increased probe diversity and tiling for a subset of genes had a large positive effect on both recovery and completeness. The resulting data produced an accurate species tree, but importantly this same data can also be applied to studies of molecular evolution and function that will allow researchers to ask larger questions in broader phylogenetic contexts. Our method demonstrates the utility of cross-species approaches for the capture of full length coding sequences, and will substantially improve the ability for researchers to conduct large-scale comparative studies of molecular evolution and function., (© The Author(s) 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

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., (© 2017. Published by The Company of Biologists Ltd.)

Three complete mitochondrial genomes of South American electric fishes (Gymnotiformes), derived from high-throughput RNA sequencing (RNA-Seq), are reported herein. We report the complete mitochondrial genome of the bluntnose knifefish Brachyhypopomus n.sp. VERD, determined from newly sequenced data. We also provide the complete mitochondrial genomes for Sternopygus arenatus and the electric eel Electrophorus electricus , assembled from previously published transcriptome data. The mitochondrial genomes of Brachyhypopomus n.sp. VERD , Sternopygus arenatus and Electrophorus electricus have 13 protein-coding genes, 1 D-loop, 2 ribosomal RNAs and 22 transfer RNAs, and are 16,547, 16,667 and 16,906 bp in length, respectively. Phylogenetic analysis of the eight available mitochondrial genomes of gymnotiform fishes shows Apteronotus to be the sister lineage of other gymnotiformes, contradicting the "Sinusoidea" hypothesis that Apteronotidae and Sternopygidae are sister taxa., Competing Interests: The authors do not report any conflict of interests, and are responsible for the content of this publication., (© 2016 The Author(s). Published by Informa UK Limited, trading as Taylor & Francis Group.)

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.

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 nocturnal origin of mammals is a longstanding hypothesis that is considered instrumental for the evolution of endothermy, a potential key innovation in this successful clade. This hypothesis is primarily based on indirect anatomical inference from fossils. Here, we reconstruct the evolutionary history of rhodopsin--the vertebrate visual pigment mediating the first step in phototransduction at low-light levels--via codon-based model tests for selection, combined with gene resurrection methods that allow for the study of ancient proteins. Rhodopsin coding sequences were reconstructed for three key nodes: Amniota, Mammalia, and Theria. When expressed in vitro, all sequences generated stable visual pigments with λMAX values similar to the well-studied bovine rhodopsin. Retinal release rates of mammalian and therian ancestral rhodopsins, measured via fluorescence spectroscopy, were significantly slower than those of the amniote ancestor, indicating altered molecular function possibly related to nocturnality. Positive selection along the therian branch suggests adaptive evolution in rhodopsin concurrent with therian ecological diversification events during the Mesozoic that allowed for an exploration of the environment at varying light levels., (© 2015 The Author(s). Evolution © 2015 The Society for the Study of Evolution.)

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.

Protein N-glycosylation is found in all domains of life and has a conserved role in glycoprotein folding and stability. In animals, glycoproteins transit through the Golgi where the N-glycans are trimmed and rebuilt with sequences that bind lectins, an innovation that greatly increases structural diversity and redundancy of glycoprotein-lectin interaction at the cell surface. Here we ask whether the natural tension between increasing diversity (glycan-protein interactions) and site multiplicity (backup and status quo) might be revealed by a phylogenic examination of glycoproteins and NXS/T(X ≠ P) N-glycosylation sites. Site loss is more likely by mutation at Asn encoded by two adenosine (A)-rich codons, while site gain is more probable by generating Ser or Thr downstream of an existing Asn. Thus mutations produce sites at novel positions more frequently than the reversal of recently lost sites, and therefore more paths though sequence space are made available to natural selection. An intra-species comparison of secretory and cytosolic proteins revealed a departure from equilibrium in sequences one-mutation-away from NXS/T and in (A) content, indicating strong selective pressures and exploration of N-glycosylation positions during vertebrate evolution. Furthermore, secretory proteins have evolved at rates proportional to N-glycosylation site number, indicating adaptive interactions between the N-glycans and underlying protein. Given the topology of the genetic code, mutation of (A) is more often nonsynonomous, and Lys, another target of many PTMs, is also encoded by two (A)-rich codons. An examination of acetyl-Lys sites in proteins indicated similar evolutionary dynamics, consistent with asymmetry of the target and recognition portions of modified sites. Our results suggest that encoding asymmetry is an ancient mechanism of evolvability that increases diversity and experimentation with PTM site positions. Strong selective pressures on PTMs may have contributed to the A+T → G+C shift in genome-wide nucleotide composition during metazoan radiation.

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.

Taphonomic biases dictate how organisms are represented in the fossil record, but their effect on studies of vertebrate diversity dynamics is poorly studied. In contrast to the high diversity and abundance of small-bodied animals in extant ecosystems, small-bodied dinosaurs are less common than their large-bodied counterparts, but it is unclear whether this reflects unique properties of dinosaurian ecosystems or relates to taphonomic biases. A new, fully domed pachycephalosaurid dinosaur, Acrotholus audeti, from the Santonian of Alberta predates incompletely domed taxa, and provides important new information on pachycephalosaur evolution and the completeness of the ornithischian fossil record. Here we provide the first empirical evidence that the diversity of small-bodied ornithischian dinosaurs is strongly underestimated based on ghost lineages and the high proportion of robust and diagnostic frontoparietal domes compared with other pachycephalosaur fossils. This suggests preservational biases have a confounding role in attempts to decipher vertebrate palaeoecology and diversity dynamics through the Mesozoic.

Historically, studies of pachycephalosaurs have recognized plesiomorphically flat-headed taxa and apomorphically domed taxa. More recently, it has been suggested that the expression of the frontoparietal dome is ontogenetic and derived from a flat-headed juvenile morphology. However, strong evidence to support this hypothesis has been lacking. Here we test this hypothesis in a large, stratigraphically constrained sample of specimens assigned to Stegoceras validum, the best known pachycephalosaur, using multiple independent lines of evidence including conserved morphology of ornamentation, landmark-based allometric analyses of frontoparietal shape, and cranial bone histology. New specimens show that the diagnostic ornamentation of the parietosquamosal bar is conserved throughout the size range of the sample, which links flat-headed specimens to domed S. validum. High-resolution CT scans of three frontoparietals reveal that vascularity decreases with size and document a pattern that is consistent with previously proposed histological changes during growth. Furthermore, aspects of dome shape and size are strongly correlated and indicative of ontogenetic growth. These results are complementary and strongly support the hypothesis that the sample represents a growth series of a single taxon. Cranial dome growth is positively allometric, proceeds from a flat-headed to a domed state, and confirms the synonymy of Ornatotholus browni as a juvenile Stegoceras. This dataset serves as the first detailed model of growth and variation in a pachycephalosaur. Flat-headed juveniles possess three characters (externally open cranial sutures, tuberculate dorsal surface texture, and open supratemporal fenestrae) that are reduced or eliminated during ontogeny. These characters also occur in putative flat-headed taxa, suggesting that they may also represent juveniles of domed taxa. However, open cranial sutures and supratemporal fenestrae are plesiomorphic within Ornithischia, and thus should be expected in the adult stage of a primitive pachycephalosaur. Additional lines of evidence will be needed to resolve the taxonomic validity of flat-headed pachycephalosaur taxa.