Ryan K. Schott

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

Significance The molecular basis of morphological and physiological adaptations in snakes is largely unknown. Here, we study these phenotypes using the genome of the Burmese python (Python molurus bivittatus), a model for extreme phenotypic plasticity and metabolic adaptation. We discovered massive rapid changes in gene expression that coordinate major changes in organ size and function after feeding. Many significantly responsive genes are associated with metabolism, development, and mammalian diseases. A striking number of genes experienced positive selection in ancestral snakes. Such genes were related to metabolism, development, lungs, eyes, heart, kidney, and skeletal structure—all highly modified features in snakes. Snake phenotypic novelty seems to be driven by the system-wide coordination of protein adaptation, gene expression, and changes in genome structure. 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.

The oldest North American pachycephalosaurid and the hidden diversity of small-bodied ornithischian dinosaurs

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.

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.

Cranial Ontogeny in Stegoceras validum (Dinosauria: Pachycephalosauria): A Quantitative Model of Pachycephalosaur Dome Growth and Variation

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.

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

Significance This study provides compelling evidence that the previously reported all-cone retina of a diurnal garter snake in fact contains a population of rod photoreceptors with the appearance, and presumably function, of cones. Our results suggest that the evolution of all-cone retinas occurred not through loss of rods but rather via the evolutionary transmutation of ancestral rods into more “cone-like” photoreceptors, to regain functionality that was lost during the early, possibly fossorial, origin of snakes. This study provides a better understanding of the process by which complex molecular/cellular structures and tissue types can evolve, and how, particularly for sensory systems, physiological constraints can be shaped by selective forces to produce evolutionary novelty. 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.

The anatomy and systematics of Colepiocephale lambei (Dinosauria: Pachycephalosauridae)

ABSTRACT Colepiocephale lambei from the Foremost Formation of Alberta is a problematic pachycephalosaurid that has recently been hypothesized as (1) synonymous with Stegoceras validum, (2) a distinct species within the Stegoceras clade, and (3) an unusual taxon more derived than Stegoceras. Conflict among these hypotheses arises from different interpretations of the parietosquamosal region of the dome. Re-evaluation of the morphology of C. lambei based on all known material from Alberta indicates that the posteromedial extension of the parietal is broken off at its base in every specimen, and what was interpreted as the posterior margin of the dome actually is the anterior margins of large supratemporal fenestrae. The squamosals are not preserved in any specimen, but the form of the parietal suggests that they were present in the typical pachycephalosaurid configuration. Colepiocephale lambei is considered to be a distinct taxon based on several autapomorphies and stratigraphic segregation from S. validum. Phylogenetic analyses that include new anatomical data confirm that C. lambei occupies a basal position within domed pachycephalosaurids, but its precise systematic position with respect to S. validum remains equivocal. This general interpretation is more consistent with its stratigraphic position as one of the oldest known members of Pachycephalosauridae. Additionally, this study results in the recognition of the first material referable to this taxon from the Judith River Formation of Montana, which results in a geographic range extension for this taxon and extends its stratigraphic range into the equivalent of the lower Oldman Formation in Alberta.

The molecular origin and evolution of dim-light vision in mammals.

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.

Encoding Asymmetry of the N-Glycosylation Motif Facilitates Glycoprotein Evolution

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.

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

ABSTRACT 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. Summary: An ancient duplication of the rhodopsin gene, rh1-2, is a functional visual pigment and is expressed in retinal photoreceptors.

Evolution of nonspectral rhodopsin function at high altitudes

Significance Protein evolution in response to different environments has long been of interest to both evolutionary biologists and biochemists. High-altitude specialist catfishes in the Andes mountains offer an opportunity to examine the molecular adaptations accompanying adaptation to cold environments. Rhodopsins and other visual pigments form the first step in vision and have long been a model system for studying the molecular basis of sensory adaptations; however, many of these studies have focused solely on spectral shifts. Recent studies suggest that other aspects of function are as important for visual performance. We demonstrate that high-altitude amino acid variants significantly accelerate RH1 kinetics. These results suggest that the activity–stability trade-off characterized in cold-adapted enzymes also affects adaptation of signaling proteins through similar molecular mechanisms. 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.