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.

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.

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 allcone 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. rhodopsin evolution | visual evolution | reptile vision | snake photoreceptors | visual pigment H ow complex structures can arise has long fascinated evolutionary biologists, and the evolution of the eye, as noted by Charles Darwin (1), is perhaps the most famous example. Within the vertebrate eye, the light-sensing photoreceptors are complex, highly specialized cellular structures that can be divided into two general types based on their distinct morphologies and functions: cones, which are active during the day and contain cone opsin pigments, and rods, which mediate dim-light vision and contain rhodopsin (RH1) (2-4). The visual pigments contained in cone photoreceptors are classified into four different subtypes that mediate vision across the visible spectrum from the UV to the red (2). Although most vertebrate retinas are duplex, containing both cones and rods, squamate reptiles (lizards and snakes) are unusual, not only in having highly variable photoreceptor morphologies, but also for several instances of the absence of an entire class of photoreceptors, resulting in simplex retinas composed of only cones or rods (4).In a seminal book published in 1942, Walls (4) hypothesized that, during evolution, vertebrate photoreceptors could transform from one type to another, a process that he termed photoreceptor "transm...

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.

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.