Large-scale snake genome analyses provide insights into vertebrate development.

Snakes are a remarkable squamate lineage with unique morphological adaptations, especially those related to the evolution of vertebrate skeletons, organs, and sensory systems. To clarify the genetic underpinnings of snake phenotypes, we assembled and analyzed 14 de novo genomes from 12 snake families. We also investigated the genetic basis of the morphological characteristics of snakes using functional experiments. We identified genes, regulatory elements, and structural variations that have potentially contributed to the evolution of limb loss, an elongated body plan, asymmetrical lungs, sensory systems, and digestive adaptations in snakes. We identified some of the genes and regulatory elements that might have shaped the evolution of vision, the skeletal system and diet in blind snakes, and thermoreception in infrared-sensitive snakes. Our study provides insights into the evolution and development of snakes and vertebrates.

Snake Venom Gland Organoids.

Wnt dependency and Lgr5 expression define multiple mammalian epithelial stem cell types. Under defined growth factor conditions, such adult stem cells (ASCs) grow as 3D organoids that recapitulate essential features of the pertinent epithelium. Here, we establish long-term expanding venom gland organoids from several snake species. The newly assembled transcriptome of the Cape coral snake reveals that organoids express high levels of toxin transcripts. Single-cell RNA sequencing of both organoids and primary tissue identifies distinct venom-expressing cell types as well as proliferative cells expressing homologs of known mammalian stem cell markers. A hard-wired regional heterogeneity in the expression of individual venom components is maintained in organoid cultures. Harvested venom peptides reflect crude venom composition and display biological activity. This study extends organoid technology to reptilian tissues and describes an experimentally tractable model system representing the snake venom gland.

Evolution of Cortical Neurogenesis in Amniotes Controlled by Robo Signaling Levels.

Cerebral cortex size differs dramatically between reptiles, birds, and mammals, owing to developmental differences in neuron production. In mammals, signaling pathways regulating neurogenesis have been identified, but genetic differences behind their evolution across amniotes remain unknown. We show that direct neurogenesis from radial glia cells, with limited neuron production, dominates the avian, reptilian, and mammalian paleocortex, whereas in the evolutionarily recent mammalian neocortex, most neurogenesis is indirect via basal progenitors. Gain- and loss-of-function experiments in mouse, chick, and snake embryos and in human cerebral organoids demonstrate that high Slit/Robo and low Dll1 signaling, via Jag1 and Jag2, are necessary and sufficient to drive direct neurogenesis. Attenuating Robo signaling and enhancing Dll1 in snakes and birds recapitulates the formation of basal progenitors and promotes indirect neurogenesis. Our study identifies modulation in activity levels of conserved signaling pathways as a primary mechanism driving the expansion and increased complexity of the mammalian neocortex during amniote evolution.

Progressive Loss of Function in a Limb Enhancer during Snake Evolution.

The evolution of body shape is thought to be tightly coupled to changes in regulatory sequences, but specific molecular events associated with major morphological transitions in vertebrates have remained elusive. We identified snake-specific sequence changes within an otherwise highly conserved long-range limb enhancer of Sonic hedgehog (Shh). Transgenic mouse reporter assays revealed that the in vivo activity pattern of the enhancer is conserved across a wide range of vertebrates, including fish, but not in snakes. Genomic substitution of the mouse enhancer with its human or fish ortholog results in normal limb development. In contrast, replacement with snake orthologs caused severe limb reduction. Synthetic restoration of a single transcription factor binding site lost in the snake lineage reinstated full in vivo function to the snake enhancer. Our results demonstrate changes in a regulatory sequence associated with a major body plan transition and highlight the role of enhancers in morphological evolution. PAPERCLIP.

Synchrotron tomography of a stem lizard elucidates early squamate anatomy.

Squamates (lizards and snakes) include more than 10,000 living species, descended from an ancestor that diverged more than 240 million years ago from that of their closest living relative, Sphenodon. However, a deficiency of fossil evidence1-7, combined with serious conflicts between molecular and morphological accounts of squamate phylogeny8-13 (but see ref. 14), has caused uncertainty about the origins and evolutionary assembly of squamate anatomy. Here we report the near-complete skeleton of a stem squamate, Bellairsia gracilis, from the Middle Jurassic epoch of Scotland, documented using high-resolution synchrotron phase-contrast tomography. Bellairsia shares numerous features of the crown group, including traits related to cranial kinesis (an important functional feature of many extant squamates) and those of the braincase and shoulder girdle. Alongside these derived traits, Bellairsia also retains inferred ancestral features including a pterygoid-vomer contact and the presence of both cervical and dorsal intercentra. Phylogenetic analyses return strong support for Bellairsia as a stem squamate, suggesting that several features that it shares with extant gekkotans are plesiomorphies, consistent with the molecular phylogenetic hypothesis that gekkotans are early-diverging squamates. We also provide confident support of stem squamate affinities for the enigmatic Oculudentavis. Our findings indicate that squamate-like functional features of the suspensorium, braincase and shoulder girdle preceded the origin of their palatal and vertebral traits and indicate the presence of advanced stem squamates as persistent components of terrestrial assemblages up to at least the middle of the Cretaceous period.

Recombinant snake antivenoms get closer to the clinic.

Snakebite envenomings kill ~100 000 victims each year and leave many more with permanent sequelae. Antivenoms have been available for more than 125 years but are in need of innovation. A new study by Khalek et al. highlights broadly neutralizing human monoclonal antibodies (mAbs) that might be used to develop recombinant antivenoms with superior therapeutic benefits.

Geometric phase predicts locomotion performance in undulating living systems across scales.

Self-propelling organisms locomote via generation of patterns of self-deformation. Despite the diversity of body plans, internal actuation schemes and environments in limbless vertebrates and invertebrates, such organisms often use similar traveling waves of axial body bending for movement. Delineating how self-deformation parameters lead to locomotor performance (e.g. speed, energy, turning capabilities) remains challenging. We show that a geometric framework, replacing laborious calculation with a diagrammatic scheme, is well-suited to discovery and comparison of effective patterns of wave dynamics in diverse living systems. We focus on a regime of undulatory locomotion, that of highly damped environments, which is applicable not only to small organisms in viscous fluids, but also larger animals in frictional fluids (sand) and on frictional ground. We find that the traveling wave dynamics used by mm-scale nematode worms and cm-scale desert dwelling snakes and lizards can be described by time series of weights associated with two principal modes. The approximately circular closed path trajectories of mode weights in a self-deformation space enclose near-maximal surface integral (geometric phase) for organisms spanning two decades in body length. We hypothesize that such trajectories are targets of control (which we refer to as "serpenoid templates"). Further, the geometric approach reveals how seemingly complex behaviors such as turning in worms and sidewinding snakes can be described as modulations of templates. Thus, the use of differential geometry in the locomotion of living systems generates a common description of locomotion across taxa and provides hypotheses for neuromechanical control schemes at lower levels of organization.

A bifurcation integrates information from many noisy ion channels and allows for milli-Kelvin thermal sensitivity in the snake pit organ.

In various biological systems, information from many noisy molecular receptors must be integrated into a collective response. A striking example is the thermal imaging organ of pit vipers. Single nerve fibers in the organ reliably respond to milli-Kelvin (mK) temperature increases, a thousand times more sensitive than their molecular sensors, thermo-transient receptor potential (TRP) ion channels. Here, we propose a mechanism for the integration of this molecular information. In our model, amplification arises due to proximity to a dynamical bifurcation, separating a regime with frequent and regular action potentials (APs), from a regime where APs are irregular and infrequent. Near the transition, AP frequency can have an extremely sharp dependence on temperature, naturally accounting for the thousand-fold amplification. Furthermore, close to the bifurcation, most of the information about temperature available in the TRP channels' kinetics can be read out from the times between consecutive APs even in the presence of readout noise. A key model prediction is that the coefficient of variation in the distribution of interspike times decreases with AP frequency, and quantitative comparison with experiments indeed suggests that nerve fibers of snakes are located very close to the bifurcation. While proximity to such bifurcation points typically requires fine-tuning of parameters, we propose that having feedback act from the order parameter (AP frequency) onto the control parameter robustly maintains the system in the vicinity of the bifurcation. This robustness suggests that similar feedback mechanisms might be found in other sensory systems which also need to detect tiny signals in a varying environment.

DeTox: a pipeline for the detection of toxins in venomous organisms.

Venomous organisms have independently evolved the ability to produce toxins 101 times during their evolutionary history, resulting in over 200 000 venomous species. Collectively, these species produce millions of toxins, making them a valuable resource for bioprospecting and understanding the evolutionary mechanisms underlying genetic diversification. RNA-seq is the preferred method for characterizing toxin repertoires, but the analysis of the resulting data remains challenging. While early approaches relied on similarity-based mapping to known toxin databases, recent studies have highlighted the importance of structural features for toxin detection. The few existing pipelines lack an integration between these complementary approaches, and tend to be difficult to run for non-experienced users. To address these issues, we developed DeTox, a comprehensive and user-friendly tool for toxin research. It combines fast execution, parallelization and customization of parameters. DeTox was tested on published transcriptomes from gastropod mollusks, cnidarians and snakes, retrieving most putative toxins from the original articles and identifying additional peptides as potential toxins to be confirmed through manual annotation and eventually proteomic analysis. By integrating a structure-based search with similarity-based approaches, DeTox allows the comprehensive characterization of toxin repertoire in poorly-known taxa. The effect of the taxonomic bias in existing databases is minimized in DeTox, as mirrored in the detection of unique and divergent toxins that would have been overlooked by similarity-based methods. DeTox streamlines toxin annotation, providing a valuable tool for efficient identification of venom components that will enhance venom research in neglected taxa.

Comparative analysis of human, rodent and snake deltavirus replication.

The recent discovery of Hepatitis D (HDV)-like viruses across a wide range of taxa led to the establishment of the Kolmioviridae family. Recent studies suggest that kolmiovirids can be satellites of viruses other than Hepatitis B virus (HBV), challenging the strict HBV/HDV-association dogma. Studying whether kolmiovirids are able to replicate in any animal cell they enter is essential to assess their zoonotic potential. Here, we compared replication of three kolmiovirids: HDV, rodent (RDeV) and snake (SDeV) deltavirus in vitro and in vivo. We show that SDeV has the narrowest and RDeV the broadest host cell range. High resolution imaging of cells persistently replicating these viruses revealed nuclear viral hubs with a peculiar RNA-protein organization. Finally, in vivo hydrodynamic delivery of viral replicons showed that both HDV and RDeV, but not SDeV, efficiently replicate in mouse liver, forming massive nuclear viral hubs. Our comparative analysis lays the foundation for the discovery of specific host factors controlling Kolmioviridae host-shifting.

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.

The population genetics of the causative agent of snake fungal disease indicate recent introductions to the USA.

Snake fungal disease (SFD; ophidiomycosis), caused by the pathogen Ophidiomyces ophiodiicola (Oo), has been documented in wild snakes in North America and Eurasia, and is considered an emerging disease in the eastern United States of America. However, a lack of historical disease data has made it challenging to determine whether Oo is a recent arrival to the USA or whether SFD emergence is due to other factors. Here, we examined the genomes of 82 Oo strains to determine the pathogen's history in the eastern USA. Oo strains from the USA formed a clade (Clade II) distinct from European strains (Clade I), and molecular dating indicated that these clades diverged too recently (approximately 2,000 years ago) for transcontinental dispersal of Oo to have occurred via natural snake movements across Beringia. A lack of nonrecombinant intermediates between clonal lineages in Clade II indicates that Oo has actually been introduced multiple times to North America from an unsampled source population, and molecular dating indicates that several of these introductions occurred within the last few hundred years. Molecular dating also indicated that the most common Clade II clonal lineages have expanded recently in the USA, with time of most recent common ancestor mean estimates ranging from 1985 to 2007 CE. The presence of Clade II in captive snakes worldwide demonstrates a potential mechanism of introduction and highlights that additional incursions are likely unless action is taken to reduce the risk of pathogen translocation and spillover into wild snake populations.

The vomeronasal organ mediates interspecies defensive behaviors through detection of protein pheromone homologs.

Potential predators emit uncharacterized chemosignals that warn receiving species of danger. Neurons that sense these stimuli remain unknown. Here we show that detection and processing of fear-evoking odors emitted from cat, rat, and snake require the function of sensory neurons in the vomeronasal organ. To investigate the molecular nature of the sensory cues emitted by predators, we isolated the salient ligands from two species using a combination of innate behavioral assays in naive receiving animals, calcium imaging, and c-Fos induction. Surprisingly, the defensive behavior-promoting activity released by other animals is encoded by species-specific ligands belonging to the major urinary protein (Mup) family, homologs of aggression-promoting mouse pheromones. We show that recombinant Mup proteins are sufficient to activate sensory neurons and initiate defensive behavior similarly to native odors. This co-option of existing sensory mechanisms provides a molecular solution to the difficult problem of evolving a variety of species-specific molecular detectors.

Coordinating tiny limbs and long bodies: Geometric mechanics of lizard terrestrial swimming.

Although typically possessing four limbs and short bodies, lizards have evolved diverse morphologies, including elongate trunks with tiny limbs. Such forms are hypothesized to aid locomotion in cluttered/fossorial environments but propulsion mechanisms (e.g., the use of body and/or limbs to interact with substrates) and potential body/limb coordination remain unstudied. Here, we use biological experiments, a geometric theory of locomotion, and robophysical models to investigate body-limb coordination in diverse lizards. Locomotor field studies in short-limbed, elongate lizards (Brachymeles and Lerista) and laboratory studies of fully limbed lizards (Uma scoparia and Sceloporus olivaceus) and a snake (Chionactis occipitalis) reveal that body-wave dynamics can be described by a combination of standing and traveling waves; the ratio of the amplitudes of these components is inversely related to the degree of limb reduction and body elongation. The geometric theory (which replaces laborious calculation with diagrams) helps explain our observations, predicting that the advantage of traveling-wave body undulations (compared with a standing wave) emerges when the dominant thrust-generation mechanism arises from the body rather than the limbs and reveals that such soil-dwelling lizards propel via "terrestrial swimming" like sand-swimming lizards and snakes. We test our hypothesis by inducing the use of traveling waves in stereotyped lizards via modulating the ground-penetration resistance. Study of a limbed/undulatory robophysical model demonstrates that a traveling wave is beneficial when propulsion is generated by body-environment interaction. Our models could be valuable in understanding functional constraints on the evolutionary processes of elongation and limb reduction as well as advancing robot designs.

No link between population isolation and speciation rate in squamate reptiles.

Rates of species formation vary widely across the tree of life and contribute to massive disparities in species richness among clades. This variation can emerge from differences in metapopulation-level processes that affect the rates at which lineages diverge, persist, and evolve reproductive barriers and ecological differentiation. For example, populations that evolve reproductive barriers quickly should form new species at faster rates than populations that acquire reproductive barriers more slowly. This expectation implicitly links microevolutionary processes (the evolution of populations) and macroevolutionary patterns (the profound disparity in speciation rate across taxa). Here, leveraging extensive field sampling from the Neotropical Cerrado biome in a biogeographically controlled natural experiment, we test the role of an important microevolutionary process-the propensity for population isolation-as a control on speciation rate in lizards and snakes. By quantifying population genomic structure across a set of codistributed taxa with extensive and phylogenetically independent variation in speciation rate, we show that broad-scale patterns of species formation are decoupled from demographic and genetic processes that promote the formation of population isolates. Population isolation is likely a critical stage of speciation for many taxa, but our results suggest that interspecific variability in the propensity for isolation has little influence on speciation rates. These results suggest that other stages of speciation-including the rate at which reproductive barriers evolve and the extent to which newly formed populations persist-are likely to play a larger role than population isolation in controlling speciation rate variation in squamates.

Caecilian Genomes Reveal the Molecular Basis of Adaptation and Convergent Evolution of Limblessness in Snakes and Caecilians.

We present genome sequences for the caecilians Geotrypetes seraphini (3.8 Gb) and Microcaecilia unicolor (4.7 Gb), representatives of a limbless, mostly soil-dwelling amphibian clade with reduced eyes, and unique putatively chemosensory tentacles. More than 69% of both genomes are composed of repeats, with retrotransposons being the most abundant. We identify 1,150 orthogroups that are unique to caecilians and enriched for functions in olfaction and detection of chemical signals. There are 379 orthogroups with signatures of positive selection on caecilian lineages with roles in organ development and morphogenesis, sensory perception, and immunity amongst others. We discover that caecilian genomes are missing the zone of polarizing activity regulatorysequence (ZRS) enhancer of Sonic Hedgehog which is also mutated in snakes. In vivo deletions have shown ZRS is required for limb development in mice, thus, revealing a shared molecular target implicated in the independent evolution of limblessness in snakes and caecilians.

Independent Recruitment of Different Types of Phospholipases A2 to the Venoms of Caenophidian Snakes: The Rise of PLA2-IIE within Pseudoboini (Dipsadidae).

Snake venoms harbor a wide and diverse array of enzymatic and nonenzymatic toxic components, allowing them to exert myriad effects on their prey. However, they appear to trend toward a few optimal compositional scaffolds, dominated by four major toxin classes: SVMPs, SVSPs, 3FTxs, and PLA2s. Nevertheless, the latter appears to be restricted to vipers and elapids, as it has never been reported as a major venom component in rear-fanged species. Here, by investigating the original transcriptomes from 19 species distributed in eight genera from the Pseudoboini tribe (Dipsadidae: Xenodontinae) and screening among seven additional tribes of Dipsadidae and three additional families of advanced snakes, we discovered that a novel type of venom PLA2, resembling a PLA2-IIE, has been recruited to the venom of some species of the Pseudoboini tribe, where it is a major component. Proteomic and functional analyses of these venoms further indicate that these PLA2s play a relevant role in the venoms from this tribe. Moreover, we reconstructed the phylogeny of PLA2s across different snake groups and show that different types of these toxins have been recruited in at least five independent events in caenophidian snakes. Additionally, we present the first compositional profiling of Pseudoboini venoms. Our results demonstrate how relevant phenotypic traits are convergently recruited by different means and from homologous and nonhomologous genes in phylogenetically and ecologically divergent snake groups, possibly optimizing venom composition to overcome diverse adaptative landscapes.

Regionalization of the vertebral column and its correlation with heart position in snakes: Implications for evolutionary pathways and morphological diversification.

Spinal regionalization has important implications for the evolution of vertebrate body plans. We determined the variation in the number and morphology of vertebrae across the vertebral column (i.e., vertebral formula) for 63 snake species representing 13 families using intracolumnar variation in vertebral shape. Vertebral counts were used to determine the position of the heart, pylorus, and left kidney for each species. Across all species we observed a conspicuous midthoracic transition in vertebral shape, indicating four developmental domains of the precloacal vertebral column (cervical, anterior thoracic, posterior thoracic, and lumbar). Using phylogenetic analyses, the boundary between the anterior and posterior thoracic vertebrae was correlated with heart position. No associations were found between shifts in morphology of the vertebral column and either the pylorus or left kidney. We observed that among taxa, the number of preapex and postapex vertebrae could change independently from one another and from changes in the total number of precloacal vertebrae. Ancestral state reconstruction of the preapex and postapex vertebrae illustrated several evolutionary pathways by which diversity in the vertebral column and heart position have been attained. In addition, no conspicuous pattern was observed among the heart, pylorus, or kidney indicating that their relative positions to each other evolve independently. We conclude that snakes exhibit four morphologically distinct regions of the vertebral column. We discuss the implications of the forebody and hindbody vertebral formula on the morphological diversification of snakes.

Comparative morphology of oral glands in snakes of the family Homalopsidae reveals substantial variation and additional independent origins of salt glands within Serpentes.

Using diffusible iodine-based contrast-enhanced computed tomography (diceCT), we examined the morphology of the oral glands of 12 species of the family Homalopsidae. Snakes of this family exhibit substantial interspecific morphological variation in their oral glands. Particular variables are the venom glands, ranging from large (e.g., Subsessor bocourti) to small (e.g., Erpeton tentaculatum). The supra- and infralabial glands are more uniform in morphology, being the second most developed in almost all the sampled species. Premaxillary glands distinct from the supralabial glands were observed in five species (Myron richardsonii, Bitia hydroides, Cantoria violacea, Fordonia leucobalia, and Gerarda prevostiana), in addition to Cerberus rynchops, the only species in which this condition was previously documented associated with the excretion of salt. In the three species of the saltwater group of homalopsids (C. violacea, F. leucobalia, and G. prevostiana), the premaxillary glands also extend posteriorly, occupying a large area above the supralabial gland, a condition not observed in any other species of snake studied thus far. Character evolution analyses indicate that premaxillary glands differentiated from the supralabial gland and evolved independently three or four times in the family, always in lineages that invaded marine habitats. Our results suggest that the differentiated premaxillary glands are likely salt glands, as is the case in C. rynchops. If corroborated, this increases to six or seven the number of independent evolutionary origins of salt glands in snakes that have undergone an evolutionary transition to marine life.