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

Molecular regulation of reversible cardiac remodeling: lessons from species with extreme physiological adaptations

Some vertebrates evolved to have a remarkable capacity for anatomical and physiological plasticity in response to environmental challenges. One example of such plasticity can be found in the ambush-hunting snakes of the genus Python, which exhibit reversible cardiac growth with feeding. The predation strategy employed by pythons is associated with months-long fasts that are arrested by ingestion of large prey. Consequently, digestion compels a dramatic increase in metabolic rate and hypertrophy of multiple organs, including the heart. In this Review, we summarize the post-prandial cardiac adaptations in pythons at the whole-heart, cellular and molecular scales. We highlight circulating factors and cellular signaling pathways that are altered during digestion to affect cardiac form and function and propose possible mechanisms that may drive the post-digestion regression of cardiac mass. Adaptive physiological cardiac hypertrophy has also been observed in other vertebrates, including in fish acclimated to cold water, birds flying at high altitudes and exercising mammals. To reveal potential evolutionarily conserved features, we summarize the molecular signatures of reversible cardiac remodeling identified in these species and compare them with those of pythons. Finally, we offer a perspective on the potential of biomimetics targeting the natural biology of pythons as therapeutics for human heart disease.

Postprandial cardiac hypertrophy is sustained by mechanics, epigenetic, and metabolic reprogramming in pythons

Constricting pythons, known for their ability to consume infrequent, massive meals, exhibit rapid and reversible cardiac hypertrophy following feeding. Our primary goal was to investigate how python hearts achieve this adaptive response after feeding. Isolated myofibrils increased force after feeding without changes in sarcomere ultrastructure and without increasing energy cost. Ca2+ transients were prolonged after feeding with no changes in myofibril Ca2+ sensitivity. Feeding reduced titin-based tension, resulting in decreased cardiac tissue stiffness. Feeding also reduced the activity of sirtuins, a metabolically linked class of histone deacetylases, and increased chromatin accessibility. Transcription factor enrichment analysis on transposase-accessible chromatin with sequencing revealed the prominent role of transcription factors Yin Yang1 and NRF1 in postfeeding cardiac adaptation. Gene expression also changed with the enrichment of translation and metabolism. Finally, metabolomics analysis and adenosine triphosphate production demonstrated that cardiac adaptation after feeding not only increased energy demand but also energy production. These findings have broad implications for our understanding of cardiac adaptation across species and hold promise for the development of innovative approaches to address cardiovascular diseases.

Computational strategies for cross-species knowledge transfer and translational biomedicine

Research organisms provide invaluable insights into human biology and diseases, serving as essential tools for functional experiments, disease modeling, and drug testing. However, evolutionary divergence between humans and research organisms hinders effective knowledge transfer across species. Here, we review state-of-the-art methods for computationally transferring knowledge across species, primarily focusing on methods that utilize transcriptome data and/or molecular networks. We introduce the term "agnology" to describe the functional equivalence of molecular components regardless of evolutionary origin, as this concept is becoming pervasive in integrative data-driven models where the role of evolutionary origin can become unclear. Our review addresses four key areas of information and knowledge transfer across species: (1) transferring disease and gene annotation knowledge, (2) identifying agnologous molecular components, (3) inferring equivalent perturbed genes or gene sets, and (4) identifying agnologous cell types. We conclude with an outlook on future directions and several key challenges that remain in cross-species knowledge transfer.

Dynamic evolutionary interplay between ontogenetic skull patterning and whole-head integration

The arrangement and morphology of the vertebrate skull reflect functional and ecological demands, making it a highly adaptable structure. However, the fundamental developmental and macroevolutionary mechanisms leading to different vertebrate skull phenotypes remain unclear. Here we exploit the morphological diversity of squamate reptiles to assess the developmental and evolutionary patterns of skull variation and covariation in the whole head. Our geometric morphometric analysis of a complex squamate ontogenetic dataset (209 specimens, 169 embryos, 44 species), covering stages from craniofacial primordia to fully ossified bones, reveals that morphological differences between snake and lizard skulls arose gradually through changes in spatial relationships (heterotopy) followed by alterations in developmental timing or rate (heterochrony). Along with dynamic spatiotemporal changes in the integration pattern of skull bone shape and topology with surrounding brain tissues and sensory organs, we identify a relatively higher phenotypic integration of the developing snake head compared with lizards. The eye, nasal cavity and Jacobson's organ are pivotal in skull morphogenesis, highlighting the importance of sensory rearrangements in snake evolution. Furthermore, our findings demonstrate the importance of early embryonic, ontogenetic and tissue interactions in shaping craniofacial evolution and ecological diversification in squamates, with implications for the nature of cranio-cerebral relations across vertebrates.

The macroevolutionary singularity of snakes

Snakes and lizards (Squamata) represent a third of terrestrial vertebrates and exhibit spectacular innovations in locomotion, feeding, and sensory processing. However, the evolutionary drivers of this radiation remain poorly known. We infer potential causes and ultimate consequences of squamate macroevolution by combining individual-based natural history observations (>60,000 animals) with a comprehensive time-calibrated phylogeny that we anchored with genomic data (5400 loci) from 1018 species. Due to shifts in the dynamics of speciation and phenotypic evolution, snakes have transformed the trophic structure of animal communities through the recurrent origin and diversification of specialized predatory strategies. Squamate biodiversity reflects a legacy of singular events that occurred during the early history of snakes and reveals the impact of historical contingency on vertebrate biodiversity.

Convergent genomic signatures associated with vertebrate viviparity

BACKGROUND

Viviparity-live birth-is a complex and innovative mode of reproduction that has evolved repeatedly across the vertebrate Tree of Life. Viviparous species exhibit remarkable levels of reproductive diversity, both in the amount of care provided by the parent during gestation, and the ways in which that care is delivered. The genetic basis of viviparity has garnered increasing interest over recent years; however, such studies are often undertaken on small evolutionary timelines, and thus are not able to address changes occurring on a broader scale. Using whole genome data, we investigated the molecular basis of this innovation across the diversity of vertebrates to answer a long held question in evolutionary biology: is the evolution of convergent traits driven by convergent genomic changes?

RESULTS

We reveal convergent changes in protein family sizes, protein-coding regions, introns, and untranslated regions (UTRs) in a number of distantly related viviparous lineages. Specifically, we identify 15 protein families showing evidence of contraction or expansion associated with viviparity. We additionally identify elevated substitution rates in both coding and noncoding sequences in several viviparous lineages. However, we did not find any convergent changes-be it at the nucleotide or protein level-common to all viviparous lineages.

CONCLUSIONS

Our results highlight the value of macroevolutionary comparative genomics in determining the genomic basis of complex evolutionary transitions. While we identify a number of convergent genomic changes that may be associated with the evolution of viviparity in vertebrates, there does not appear to be a convergent molecular signature shared by all viviparous vertebrates. Ultimately, our findings indicate that a complex trait such as viviparity likely evolves with changes occurring in multiple different pathways.

Dissecting an ancient stress resistance trait syndrome in the compost yeast Kluyveromyces marxianus

In the search to understand how evolution builds new traits, ancient events are often the hardest to dissect. Species-unique traits pose a particular challenge for geneticists-cases in which a character arose long ago and, in the modern day, is conserved within a species, distinguishing it from reproductively isolated relatives. In this work, we have developed the budding yeast genus Kluyveromyces as a model for mechanistic dissection of trait variation across species boundaries. Phenotypic profiling revealed robust heat and chemical-stress tolerance phenotypes that distinguished the compost yeast K. marxianus from the rest of the clade. We used culture-based, transcriptomic, and genetic approaches to characterize the metabolic requirements of the K. marxianus trait syndrome. We then generated a population-genomic resource for K. marxianus and harnessed it in molecular-evolution analyses, which found hundreds of housekeeping genes with evidence for adaptive protein variation unique to this species. Our data support a model in which, in the distant past, K. marxianus underwent a vastly complex remodeling of its proteome to achieve stress resistance. Such a polygenic architecture, involving nucleotide-level allelic variation on a massive scale, is consistent with theoretical models of the mechanisms of long-term adaptation, and suggests principles of broad relevance for interspecies trait genetics.

A lizard is never late: Squamate genomics as a recent catalyst for understanding sex chromosome and microchromosome evolution

In 2011, the first high-quality genome assembly of a squamate reptile (lizard or snake) was published for the green anole. Dozens of genome assemblies were subsequently published over the next decade, yet these assemblies were largely inadequate for answering fundamental questions regarding genome evolution in squamates due to their lack of contiguity or annotation. As the "genomics age" was beginning to hit its stride in many organismal study systems, progress in squamates was largely stagnant following the publication of the green anole genome. In fact, zero high-quality (chromosome-level) squamate genomes were published between the years 2012 and 2017. However, since 2018, an exponential increase in high-quality genome assemblies has materialized with 24 additional high-quality genomes published for species across the squamate tree of life. As the field of squamate genomics is rapidly evolving, we provide a systematic review from an evolutionary genomics perspective. We collated a near-complete list of publicly available squamate genome assemblies from more than half-a-dozen international and third-party repositories and systematically evaluated them with regard to their overall quality, phylogenetic breadth, and usefulness for continuing to provide accurate and efficient insights into genome evolution across squamate reptiles. This review both highlights and catalogs the currently available genomic resources in squamates and their ability to address broader questions in vertebrates, specifically sex chromosome and microchromosome evolution, while addressing why squamates may have received less historical focus and has caused their progress in genomics to lag behind peer taxa.

The Need for Next-Generation Antivenom for Snakebite Envenomation in India

The limitations posed by currently available antivenoms have emphasized the need for alternative treatments to counteract snakebite envenomation. Even though exact epidemiological data are lacking, reports have indicated that most global snakebite deaths are reported in India. Among the many problems associated with snakebite envenomation, issues related to the availability of safer and more efficient antivenoms are of primary concern. Since India has the highest number of global snakebite deaths, efforts should be made to reduce the burden associated with snakebite envenoming. Alternative methods, including aptamers, camel antivenoms, phage display techniques for generating high-affinity antibodies and antibody fragments, small-molecule inhibitors, and natural products, are currently being investigated for their effectiveness. These alternative methods have shown promise in vitro, but their in vivo effectiveness should also be evaluated. In this review, the issues associated with Indian polyvalent antivenoms in neutralizing venom components from geographically distant species are discussed in detail. In a nutshell, this review gives an overview of the current drawbacks of using animal-derived antivenoms and several alternative strategies that are currently being widely explored.

A frameshift variant in the melanophilin gene is associated with loss of pigment from shed skin in ball pythons ( Python regius )

Melanophilin is a myosin adaptor required for transporting the pigment melanin within cells. Loss of melanophilin in fish, birds, and mammals causes pigmentation defects, but little is known about the role of melanophilin in non-avian reptiles. Here we show that a frameshift in the melanophilin gene in ball python ( P. regius ) is associated with loss of pigment from shed skin. This variant is predicted to remove the myosin-binding domain of melanophilin and thereby impair transport of melanin-containing organelles. Our study represents the first description of a melanophilin variant in a non-avian reptile and confirms the role of melanophilin across vertebrates.

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.

Stripes and loss of color in ball pythons (Python regius) are associated with variants affecting endothelin signaling

Color patterns in nonavian reptiles are beautifully diverse, but little is known about the genetics and development of these patterns. Here, we investigated color patterning in pet ball pythons (Python regius), which have been bred to show color phenotypes that differ dramatically from the wildtype form. We report that several color phenotypes in pet animals are associated with putative loss-of-function variants in the gene encoding endothelin receptor EDNRB1: (1) frameshift variants in EDNRB1 are associated with conversion of the normal mottled color pattern to skin that is almost fully white, (2) missense variants affecting conserved sites of the EDNRB1 protein are associated with dorsal, longitudinal stripes, and (3) substitutions at EDNRB1 splice donors are associated with subtle changes in patterning compared to wildtype. We propose that these phenotypes are caused by loss of specialized color cells (chromatophores), with loss ranging from severe (fully white) to moderate (dorsal striping) to mild (subtle changes in patterning). Our study is the first to describe variants affecting endothelin signaling in a nonavian reptile and suggests that reductions in endothelin signaling in ball pythons can produce a variety of color phenotypes, depending on the degree of color cell loss.

The State of Squamate Genomics: Past, Present, and Future of Genome Research in the Most Speciose Terrestrial Vertebrate Order

Squamates include more than 11,000 extant species of lizards, snakes, and amphisbaenians, and display a dazzling diversity of phenotypes across their over 200-million-year evolutionary history on Earth. Here, we introduce and define squamates (Order Squamata) and review the history and promise of genomic investigations into the patterns and processes governing squamate evolution, given recent technological advances in DNA sequencing, genome assembly, and evolutionary analysis. We survey the most recently available whole genome assemblies for squamates, including the taxonomic distribution of available squamate genomes, and assess their quality metrics and usefulness for research. We then focus on disagreements in squamate phylogenetic inference, how methods of high-throughput phylogenomics affect these inferences, and demonstrate the promise of whole genomes to settle or sustain persistent phylogenetic arguments for squamates. We review the role transposable elements play in vertebrate evolution, methods of transposable element annotation and analysis, and further demonstrate that through the understanding of the diversity, abundance, and activity of transposable elements in squamate genomes, squamates can be an ideal model for the evolution of genome size and structure in vertebrates. We discuss how squamate genomes can contribute to other areas of biological research such as venom systems, studies of phenotypic evolution, and sex determination. Because they represent more than 30% of the living species of amniote, squamates deserve a genome consortium on par with recent efforts for other amniotes (i.e., mammals and birds) that aim to sequence most of the extant families in a clade.

The genome assembly and annotation of the many-banded krait, Bungarus multicinctus

Snakes are a vital component of wildlife resources and are widely distributed across the globe. The many-banded krait Bungarus multicinctus is a highly venomous snake found across Southern Asia and central and southern China. Snakes are an ancient reptile group, and their genomes can provide important clues for understanding the evolutionary history of reptiles. Additionally, genomic resources play a crucial role in comprehending the evolution of all species. However, snake genomic resources are still scarce. Here, we present a highly contiguous genome of B. multicinctus with a size of 1.51 Gb. The genome contains a repeat content of 40.15%, with a total length exceeding 620 Mb. Additionally, we annotated a total of 24,869 functional genes. This research is of great significance for comprehending the evolution of B. multicinctus and provides genomic information on the genes involved in venom gland functions.

Piebaldism and chromatophore development in reptiles are linked to the tfec gene

Reptiles display great diversity in color and pattern, yet much of what we know about vertebrate coloration comes from classic model species such as the mouse and zebrafish.1,2,3,4 Captive-bred ball pythons (Python regius) exhibit a remarkable degree of color and pattern variation. Despite the wide range of Mendelian color phenotypes available in the pet trade, ball pythons remain an overlooked species in pigmentation research. Here, we investigate the genetic basis of the recessive piebald phenotype, a pattern defect characterized by patches of unpigmented skin (leucoderma). We performed whole-genome sequencing and used a case-control approach to discover a nonsense mutation in the gene encoding the transcription factor tfec, implicating this gene in the leucodermic patches in ball pythons. We functionally validated tfec in a lizard model (Anolis sagrei) using the gene editing CRISPR/Cas9 system and TEM imaging of skin. Our findings show that reading frame mutations in tfec affect coloration and lead to a loss of iridophores in Anolis, indicating that tfec is required for chromatophore development. This study highlights the value of captive-bred ball pythons as a model species for accelerating discoveries on the genetic basis of vertebrate coloration.

Genome Evolution and the Future of Phylogenomics of Non-Avian Reptiles

Non-avian reptiles comprise a large proportion of amniote vertebrate diversity, with squamate reptiles-lizards and snakes-recently overtaking birds as the most species-rich tetrapod radiation. Despite displaying an extraordinary diversity of phenotypic and genomic traits, genomic resources in non-avian reptiles have accumulated more slowly than they have in mammals and birds, the remaining amniotes. Here we review the remarkable natural history of non-avian reptiles, with a focus on the physical traits, genomic characteristics, and sequence compositional patterns that comprise key axes of variation across amniotes. We argue that the high evolutionary diversity of non-avian reptiles can fuel a new generation of whole-genome phylogenomic analyses. A survey of phylogenetic investigations in non-avian reptiles shows that sequence capture-based approaches are the most commonly used, with studies of markers known as ultraconserved elements (UCEs) especially well represented. However, many other types of markers exist and are increasingly being mined from genome assemblies in silico, including some with greater information potential than UCEs for certain investigations. We discuss the importance of high-quality genomic resources and methods for bioinformatically extracting a range of marker sets from genome assemblies. Finally, we encourage herpetologists working in genomics, genetics, evolutionary biology, and other fields to work collectively towards building genomic resources for non-avian reptiles, especially squamates, that rival those already in place for mammals and birds. Overall, the development of this cross-amniote phylogenomic tree of life will contribute to illuminate interesting dimensions of biodiversity across non-avian reptiles and broader amniotes.

Burmese pythons in Florida: A synthesis of biology, impacts, and management tools

Burmese pythons (Python molurus bivittatus) are native to southeastern Asia, however, there is an established invasive population inhabiting much of southern Florida throughout the Greater Everglades Ecosystem. Pythons have severely impacted native species and ecosystems in Florida and represent one of the most intractable invasive-species management issues across the globe. The difficulty stems from a unique combination of inaccessible habitat and the cryptic and resilient nature of pythons that thrive in the subtropical environment of southern Florida, rendering them extremely challenging to detect. Here we provide a comprehensive review and synthesis of the science relevant to managing invasive Burmese pythons. We describe existing control tools and review challenges to productive research, identifying key knowledge gaps that would improve future research and decision making for python control.

Reference genome of the rubber boa, Charina bottae (Serpentes: Boidae)

The rubber boa, Charina bottae is a semi-fossorial, cold-temperature adapted snake that ranges across the wetter and cooler ecoregions of the California Floristic Province. The rubber boa is 1 of 2 species in the family Boidae native to California and currently has 2 recognized subspecies, the Northern rubber boa C. bottae bottae and the Southern rubber boa C. bottae umbratica. Recent genomic work on C. bottae indicates that these 2 subspecies are collectively composed of 4 divergent lineages that separated during the late Miocene. Analysis of habitat suitability indicates that C. bottae umbratica montane sky-island populations from southern California will lose the majority of their habit over the next 70 yr, and is listed as Threatened under the California Endangered Species Act. Here, we report a new, chromosome-level assembly of C. bottae bottae as part of the California Conservation Genomics Project (CCGP). Consistent with the reference genome strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology to produce a de novo assembled genome. The assembly comprises 289 scaffolds covering 1,804,944,895 bp, has a contig N50 of 37.3 Mb, a scaffold N50 of 97 Mb, and BUSCO completeness score of 96.3%, and represents the first reference genome for the Boidae snake family. This genome will enable studies of genetic differentiation and connectivity among C. bottae bottae and C. bottae umbratica populations across California and help manage locally endemic lineages as they confront challenges from human-induced climate warming, droughts, and wildfires across California.

Structure and evolution of the squamate major histocompatibility complex as revealed by two Anolis lizard genomes

The major histocompatibility complex (MHC) is an important genomic region for adaptive immunity and has long been studied in ecological and evolutionary contexts, such as disease resistance and mate and kin selection. The MHC has been investigated extensively in mammals and birds but far less so in squamate reptiles, the third major radiation of amniotes. We localized the core MHC genomic region in two squamate species, the green anole (Anolis carolinensis) and brown anole (A. sagrei), and provide the first detailed characterization of the squamate MHC, including the presence and ordering of known MHC genes in these species and comparative assessments of genomic structure and composition in MHC regions. We find that the Anolis MHC, located on chromosome 2 in both species, contains homologs of many previously-identified mammalian MHC genes in a single core MHC region. The repetitive element composition in anole MHC regions was similar to those observed in mammals but had important distinctions, such as higher proportions of DNA transposons. Moreover, longer introns and intergenic regions result in a much larger squamate MHC region (11.7 Mb and 24.6 Mb in the green and brown anole, respectively). Evolutionary analyses of MHC homologs of anoles and other representative amniotes uncovered generally monophyletic relationships between species-specific homologs and a loss of the peptide-binding domain exon 2 in one of two mhc2β gene homologs of each anole species. Signals of diversifying selection in each anole species was evident across codons of mhc1, many of which appear functionally relevant given known structures of this protein from the green anole, chicken, and human. Altogether, our investigation fills a major gap in understanding of amniote MHC diversity and evolution and provides an important foundation for future squamate-specific or vertebrate-wide investigations of the MHC.

A community-science approach identifies genetic variants associated with three color morphs in ball pythons (Python regius)

Color morphs in ball pythons (Python regius) provide a unique and largely untapped resource for understanding the genetics of coloration in reptiles. Here we use a community-science approach to investigate the genetics of three color morphs affecting production of the pigment melanin. These morphs-Albino, Lavender Albino, and Ultramel-show a loss of melanin in the skin and eyes, ranging from severe (Albino) to moderate (Lavender Albino) to mild (Ultramel). To identify genetic variants causing each morph, we recruited shed skins of pet ball pythons via social media, extracted DNA from the skins, and searched for putative loss-of-function variants in homologs of genes controlling melanin production in other vertebrates. We report that the Albino morph is associated with missense and non-coding variants in the gene TYR. The Lavender Albino morph is associated with a deletion in the gene OCA2. The Ultramel morph is associated with a missense variant and a putative deletion in the gene TYRP1. Our study is one of the first to identify genetic variants associated with color morphs in ball pythons and shows that pet samples recruited from the community can provide a resource for genetic studies in this species.

Genome of the endangered Guatemalan Beaded Lizard, Heloderma charlesbogerti, reveals evolutionary relationships of squamates and declines in effective population sizes

Many lizard species face extinction due to worldwide climate change. The Guatemalan Beaded Lizard, Heloderma charlesbogerti, is a member of the Family Helodermatidae that may be particularly imperiled; fewer than 600 mature individuals are believed to persist in the wild. In addition, H. charlesbogerti lizards are phenotypically remarkable. They are large in size, charismatically patterned, and possess a venomous bite. Here, we report the draft genome of the Guatemalan Beaded Lizard using DNA from a wild-caught individual. The assembled genome totals 2.31 Gb in length, similar in size to the genomes of related species. Single-copy orthologs were used to produce a novel molecular phylogeny, revealing that the Guatemalan Beaded Lizard falls into a clade with the Asian Glass Lizard (Anguidae) and in close association with the Komodo Dragon (Varanidae) and the Chinese Crocodile Lizard (Shinisauridae). In addition, we identified 31,411 protein-coding genes within the genome. Of the genes identified, we found 504 that evolved with a differential constraint on the branch leading to the Guatemalan Beaded Lizard. Lastly, we identified a decline in the effective population size of the Guatemalan Beaded Lizard approximately 400,000 years ago, followed by a stabilization before starting to dwindle again 60,000 years ago. The results presented here provide important information regarding a highly endangered, venomous reptile that can be used in future conservation, functional genetic, and phylogenetic analyses.

Comparative mitogenomic analyses provide evolutionary insights into the retrolateral tibial apophysis clade (Araneae: Entelegynae)

The retrolateral tibial apophysis (RTA) clade is the largest spider lineage within Araneae. To better understand the diversity and evolution, we newly determined mitogenomes of ten RTA species from six families and performed a comparative mitogenomics analysis by combining them with 40 sequenced RTA mitogenomes available on GenBank. The ten mitogenomes encoded 37 typical mitochondrial genes and included a large non-coding region (putative control region). Nucleotide composition and codon usage were well conserved within the RTA clade, whereas diversity in sequence length and structural features was observed in control region. A reversal of strand asymmetry in nucleotide composition, i.e., negative AT-skews and positive GC-skews, was observed in each RTA species, likely resulting from mitochondrial gene rearrangements. All protein-coding genes were evolving under purifying selection, except for atp8 whose Ka/Ks was larger than 1, possibly due to positive selection or selection relaxation. Both mutation pressure and natural selection might contribute to codon usage bias of 13 protein-coding genes in the RTA lineage. Phylogenetic analyses based on mitogenomic data recovered a family-level phylogeny within the RTA; {[(Oval calamistrum clade, Dionycha), Marronoid clade], Sparassidae}. This study characterized RTA mitogenomes and provided some new insights into the phylogeny and evolution of the RTA clade.

The Rattlesnake W Chromosome: A GC-Rich Retroelement Refugium with Retained Gene Function Across Ancient Evolutionary Strata

Sex chromosomes diverge after the establishment of recombination suppression, resulting in differential sex-linkage of genes involved in genetic sex determination and dimorphic traits. This process produces systems of male or female heterogamety wherein the Y and W chromosomes are only present in one sex and are often highly degenerated. Sex-limited Y and W chromosomes contain valuable information about the evolutionary transition from autosomes to sex chromosomes, yet detailed characterizations of the structure, composition, and gene content of sex-limited chromosomes are lacking for many species. In this study, we characterize the female-specific W chromosome of the prairie rattlesnake (Crotalus viridis) and evaluate how recombination suppression and other processes have shaped sex chromosome evolution in ZW snakes. Our analyses indicate that the rattlesnake W chromosome is over 80% repetitive and that an abundance of GC-rich mdg4 elements has driven an overall high degree of GC-richness despite a lack of recombination. The W chromosome is also highly enriched for repeat sequences derived from endogenous retroviruses and likely acts as a "refugium" for these and other retroelements. We annotated 219 putatively functional W-linked genes across at least two evolutionary strata identified based on estimates of sequence divergence between Z and W gametologs. The youngest of these strata is relatively gene-rich, however gene expression across strata suggests retained gene function amidst a greater degree of degeneration following ancient recombination suppression. Functional annotation of W-linked genes indicates a specialization of the W chromosome for reproductive and developmental function since recombination suppression from the Z chromosome.

Identification of Daboia siamensis venome using integrated multi-omics data

Snakebite, classified by World Health Organization as a neglected tropical disease, causes more than 100,000 deaths and 2 million injuries per year. Currently, available antivenoms do not bind with strong specificity to target toxins, which means that severe complications can still occur despite treatment. Moreover, the cost of antivenom is expensive. Knowledge of venom compositions is fundamental for producing a specific antivenom that has high effectiveness, low side effects, and ease of manufacture. With advances in mass spectrometry techniques, venom proteomes can now be analyzed in great depth at high efficiency. However, these techniques require genomic and transcriptomic data for interpreting mass spectrometry data. This study aims to establish and incorporate genomics, transcriptomics, and proteomics data to study venomics of a venomous snake, Daboia siamensis. Multiple proteins that have not been reported as venom components of this snake such as hyaluronidase-1, phospholipase B, and waprin were discovered. Thus, multi-omics data are advantageous for venomics studies. These findings will be valuable not only for antivenom production but also for the development of novel therapeutics.

The draft genome sequence of the Brahminy blindsnake Indotyphlops braminus

Blindsnakes of infraoder Scolecophidia (order Squamata) are the most basal group of extant snakes, comprising of more than 450 species with ecological and morphological features highly specialized to underground living. The Brahminy blindsnake, Indotyphlops braminus, is the only known obligate parthenogenetic species of snakes. Although the origin of I. braminus is thought to be South Asia, this snake has attracted worldwide attention as an alien species, as it has been introduced to all continents except Antarctica. In this study, we present the first draft genome assembly and annotation of I. braminus. We generated approximately 480 Gbp of sequencing data and produced a draft genome with a total length of 1.86 Gbp and N50 scaffold size of 1.25 Mbp containing 89.3% of orthologs conserved in Sauropsida. We also identified 0.98 Gbp (52.82%) of repetitive genome sequences and a total of 23,560 protein-coding genes. The first draft genome of I. braminus will facilitate further study of snake evolution as well as help to understand the emergence mechanism of parthenogenetic vertebrates.

The structural and functional divergence of a neglected three-finger toxin subfamily in lethal elapids

Bungarus multicinctus is a widely distributed and medically important elapid snake that produces lethal neurotoxic venom. To study and enhance existing antivenom, we explore the complete repertoire of its toxin genes based on de novo chromosome-level assembly and multi-tissue transcriptome data. Comparative genomic analyses suggest that the three-finger toxin family (3FTX) may evolve through the neofunctionalization of flanking LY6E. A long-neglected 3FTX subfamily (i.e., MKA-3FTX) is also investigated. Only one MKA-3FTX gene, which evolves a different protein conformation, is under positive selection and actively transcribed in the venom gland, functioning as a major toxin effector together with MKT-3FTX subfamily homologs. Furthermore, this lethal snake may acquire self-resistance to its β-bungarotoxin via amino acid replacements on fast-evolving KCNA2. This study provides valuable resources for further evolutionary and structure-function studies of snake toxins, which are fundamental for the development of effective antivenoms and drug candidates.

The rise of genomics in snake venom research: recent advances and future perspectives

Snake venoms represent a danger to human health, but also a gold mine of bioactive proteins that can be harnessed for drug discovery purposes. The evolution of snakes and their venom has been studied for decades, particularly via traditional morphological and basic genetic methods alongside venom proteomics. However, while the field of genomics has matured rapidly over the past 2 decades, owing to the development of next-generation sequencing technologies, snake genomics remains in its infancy. Here, we provide an overview of the state of the art in snake genomics and discuss its potential implications for studying venom evolution and toxinology. On the basis of current knowledge, gene duplication and positive selection are key mechanisms in the neofunctionalization of snake venom proteins. This makes snake venoms important evolutionary drivers that explain the remarkable venom diversification and adaptive variation observed in these reptiles. Gene duplication and neofunctionalization have also generated a large number of repeat sequences in snake genomes that pose a significant challenge to DNA sequencing, resulting in the need for substantial computational resources and longer sequencing read length for high-quality genome assembly. Fortunately, owing to constantly improving sequencing technologies and computational tools, we are now able to explore the molecular mechanisms of snake venom evolution in unprecedented detail. Such novel insights have the potential to affect the design and development of antivenoms and possibly other drugs, as well as provide new fundamental knowledge on snake biology and evolution.

Snake Venomics: Fundamentals, Recent Updates, and a Look to the Next Decade

Venomic research, powered by techniques adapted from proteomics, transcriptomics, and genomics, seeks to unravel the diversity and complexity of venom through which knowledge can be applied in the treatment of envenoming, biodiscovery, and conservation. Snake venom proteomics is most extensively studied, but the methods varied widely, creating a massive amount of information which complicates data comparison and interpretation. Advancement in mass spectrometry technology, accompanied by growing databases and sophisticated bioinformatic tools, has overcome earlier limitations of protein identification. The progress, however, remains challenged by limited accessibility to samples, non-standardized quantitative methods, and biased interpretation of -omic data. Next-generation sequencing (NGS) technologies enable high-throughput venom-gland transcriptomics and genomics, complementing venom proteomics by providing deeper insights into the structural diversity, differential expression, regulation and functional interaction of the toxin genes. Venomic tissue sampling is, however, difficult due to strict regulations on wildlife use and transfer of biological materials in some countries. Limited resources for techniques and funding are among other pertinent issues that impede the progress of venomics, particularly in less developed regions and for neglected species. Genuine collaboration between international researchers, due recognition of regional experts by global organizations (e.g., WHO), and improved distribution of research support, should be embraced.

Distinguishing Evolutionary Conservation from Derivedness

While the concept of “evolutionary conservation” has enabled biologists to explain many ancestral features and traits, it has also frequently been misused to evaluate the degree of changes from a common ancestor, or “derivedness”. We propose that the distinction of these two concepts allows us to properly understand phenotypic and organismal evolution. From a methodological aspect, “conservation” mainly considers genes or traits which species have in common, while “derivedness” additionally covers those that are not commonly shared, such as novel or lost traits and genes to evaluate changes from the time of divergence from a common ancestor. Due to these differences, while conservation-oriented methods are effective in identifying ancestral features, they may be prone to underestimating the overall changes accumulated during the evolution of certain lineages. Herein, we describe our recently developed method, “transcriptomic derivedness index”, for estimating the phenotypic derivedness of embryos with a molecular approach using the whole-embryonic transcriptome as a phenotype. Although echinoderms are often considered as highly derived species, our analyses with this method showed that their embryos, at least at the transcriptomic level, may not be much more derived than those of chordates. We anticipate that the future development of derivedness-oriented methods could provide quantitative indicators for finding highly/lowly evolvable traits.

Catecholamines are key modulators of ventricular repolarization patterns in the ball python (Python regius)

Ectothermic vertebrates experience daily changes in body temperature, and anecdotal observations suggest these changes affect ventricular repolarization such that the T-wave in the ECG changes polarity. Mammals, in contrast, can maintain stable body temperatures, and their ventricular repolarization is strongly modulated by changes in heart rate and by sympathetic nervous system activity. The aim of this study was to assess the role of body temperature, heart rate, and circulating catecholamines on local repolarization gradients in the ectothermic ball python (Python regius). We recorded body-surface electrocardiograms and performed open-chest high-resolution epicardial mapping while increasing body temperature in five pythons, in all of which there was a change in T-wave polarity. However, the vector of repolarization differed between individuals, and only a subset of leads revealed T-wave polarity change. RNA sequencing revealed regional differences related to adrenergic signaling. In one denervated and Ringer's solution-perfused heart, heating and elevated heart rates did not induce change in T-wave polarity, whereas noradrenaline did. Accordingly, electrocardiograms in eight awake pythons receiving intra-arterial infusion of the β-adrenergic receptor agonists adrenaline and isoproterenol revealed T-wave inversion in most individuals. Conversely, blocking the β-adrenergic receptors using propranolol prevented T-wave change during heating. Our findings indicate that changes in ventricular repolarization in ball pythons are caused by increased tone of the sympathetic nervous system, not by changes in temperature. Therefore, ventricular repolarization in both pythons and mammals is modulated by evolutionary conserved mechanisms involving catecholaminergic stimulation.

Identification of an integrated stress and growth response signaling switch that directs vertebrate intestinal regeneration

BACKGROUND

Snakes exhibit extreme intestinal regeneration following months-long fasts that involves unparalleled increases in metabolism, function, and tissue growth, but the specific molecular control of this process is unknown. Understanding the mechanisms that coordinate these regenerative phenotypes provides valuable opportunities to understand critical pathways that may control vertebrate regeneration and novel perspectives on vertebrate regenerative capacities.

RESULTS

Here, we integrate a comprehensive set of phenotypic, transcriptomic, proteomic, and phosphoproteomic data from boa constrictors to identify the mechanisms that orchestrate shifts in metabolism, nutrient uptake, and cellular stress to direct phases of the regenerative response. We identify specific temporal patterns of metabolic, stress response, and growth pathway activation that direct regeneration and provide evidence for multiple key central regulatory molecules kinases that integrate these signals, including major conserved pathways like mTOR signaling and the unfolded protein response.

CONCLUSION

Collectively, our results identify a novel switch-like role of stress responses in intestinal regeneration that forms a primary regulatory hub facilitating organ regeneration and could point to potential pathways to understand regenerative capacity in vertebrates.

Eye-Transcriptome and Genome-Wide Sequencing for Scolecophidia: Implications for Inferring the Visual System of the Ancestral Snake

Molecular genetic data have recently been incorporated in attempts to reconstruct the ecology of the ancestral snake, though this has been limited by a paucity of data for one of the two main extant snake taxa, the highly fossorial Scolecophidia. Here we present and analyze vision genes from the first eye-transcriptomic and genome-wide data for Scolecophidia, for Anilios bicolor, and A. bituberculatus, respectively. We also present immunohistochemistry data for retinal anatomy and visual opsin-gene expression in Anilios. Analyzed in the context of 19 lepidosaurian genomes and 12 eye transcriptomes, the new genome-wide and transcriptomic data provide evidence for a much more reduced visual system in Anilios than in non-scolecophidian (=alethinophidian) snakes and in lizards. In Anilios, there is no evidence of the presence of 7 of the 12 genes associated with alethinophidian photopic (cone) phototransduction. This indicates extensive gene loss and many of these candidate gene losses occur also in highly fossorial mammals with reduced vision. Although recent phylogenetic studies have found evidence for scolecophidian paraphyly, the loss in Anilios of visual genes that are present in alethinophidians implies that the ancestral snake had a better-developed visual system than is known for any extant scolecophidian.

A non-lethal method for studying scorpion venom gland transcriptomes, with a review of potentially suitable taxa to which it can be applied

Scorpion venoms are mixtures of proteins, peptides and small molecular compounds with high specificity for ion channels and are therefore considered to be promising candidates in the venoms-to-drugs pipeline. Transcriptomes are important tools for studying the composition and expression of scorpion venom. Unfortunately, studying the venom gland transcriptome traditionally requires sacrificing the animal and therefore is always a single snapshot in time. This paper describes a new way of generating a scorpion venom gland transcriptome without sacrificing the animal, thereby allowing the study of the transcriptome at various time points within a single individual. By comparing these venom-derived transcriptomes to the traditional whole-telson transcriptomes we show that the relative expression levels of the major toxin classes are similar. We further performed a multi-day extraction using our proposed method to show the possibility of doing a multiple time point transcriptome analysis. This allows for the study of patterns of toxin gene activation over time a single individual, and allows assessment of the effects of diet, season and other factors that are known or likely to influence intraindividual venom composition. We discuss the gland characteristics that may allow this method to be successful in scorpions and provide a review of other venomous taxa to which this method may potentially be successfully applied.

Clutch may predict growth of hatchling Burmese pythons better than food availability or sex

Identifying which environmental and genetic factors affect growth pattern phenotypes can help biologists predict how organisms distribute finite energy resources in response to varying environmental conditions and physiological states. This information may be useful for monitoring and managing populations of cryptic, endangered, and invasive species. Consequently, we assessed the effects of food availability, clutch, and sex on the growth of invasive Burmese pythons (Python bivittatus Kuhl) from the Greater Everglades Ecosystem in Florida, USA. Though little is known from the wild, Burmese pythons have been physiological model organisms for decades, with most experimental research sourcing individuals from the pet trade. Here, we used 60 hatchlings collected as eggs from the nests of two wild pythons, assigned them to High or Low feeding treatments, and monitored growth and meal consumption for 12 weeks, a period when pythons are thought to grow very rapidly. None of the 30 hatchlings that were offered food prior to their fourth week post-hatching consumed it, presumably because they were relying on internal yolk stores. Although only two clutches were used in the experiment, we found that nearly all phenotypic variation was explained by clutch rather than feeding treatment or sex. Hatchlings from clutch 1 (C1) grew faster and were longer, heavier, in better body condition, ate more frequently, and were bolder than hatchlings from clutch 2 (C2), regardless of food availability. On average, C1 and C2 hatchling snout-vent length (SVL) and weight grew 0.15 cm d⁻¹ and 0.10 cm d⁻¹, and 0.20 g d⁻¹ and 0.03 g d⁻¹, respectively. Additional research may be warranted to determine whether these effects remain with larger clutch sample sizes and to identify the underlying mechanisms and fitness implications of this variation to help inform risk assessments and management.

This article has an associated First Person interview with the first author of the paper.

Summary: Hatchling pythons from an invasive population displayed substantial phenotypic variation in morphometrics, growth rates, and behaviors. This information may be useful for managing populations of cryptic, endangered, and invasive species.

A chromosome-level genome assembly for the eastern fence lizard (Sceloporus undulatus), a reptile model for physiological and evolutionary ecology

Background

High-quality genomic resources facilitate investigations into behavioral ecology, morphological and physiological adaptations, and the evolution of genomic architecture. Lizards in the genus Sceloporus have a long history as important ecological, evolutionary, and physiological models, making them a valuable target for the development of genomic resources.

Findings

We present a high-quality chromosome-level reference genome assembly, SceUnd1.0 (using 10X Genomics Chromium, HiC, and Pacific Biosciences data), and tissue/developmental stage transcriptomes for the eastern fence lizard, Sceloporus undulatus. We performed synteny analysis with other snake and lizard assemblies to identify broad patterns of chromosome evolution including the fusion of micro- and macrochromosomes. We also used this new assembly to provide improved reference-based genome assemblies for 34 additional Sceloporus species. Finally, we used RNAseq and whole-genome resequencing data to compare 3 assemblies, each representing an increased level of cost and effort: Supernova Assembly with data from 10X Genomics Chromium, HiRise Assembly that added data from HiC, and PBJelly Assembly that added data from Pacific Biosciences sequencing. We found that the Supernova Assembly contained the full genome and was a suitable reference for RNAseq and single-nucleotide polymorphism calling, but the chromosome-level scaffolds provided by the addition of HiC data allowed synteny and whole-genome association mapping analyses. The subsequent addition of PacBio data doubled the contig N50 but provided negligible gains in scaffold length.

Conclusions

These new genomic resources provide valuable tools for advanced molecular analysis of an organism that has become a model in physiology and evolutionary ecology.

The Genomic Processes of Biological Invasions: From Invasive Species to Cancer Metastases and Back Again

The concept of invasion is useful across a broad range of contexts, spanning from the fine scale landscape of cancer tumors up to the broader landscape of ecosystems. Invasion biology provides extraordinary opportunities for studying the mechanistic basis of contemporary evolution at the molecular level. Although the field of invasion genetics was established in ecology and evolution more than 50 years ago, there is still a limited understanding of how genomic level processes translate into invasive phenotypes across different taxa in response to complex environmental conditions. This is largely because the study of most invasive species is limited by information about complex genome level processes. We lack good reference genomes for most species. Rigorous studies to examine genomic processes are generally too costly. On the contrary, cancer studies are fortified with extensive resources for studying genome level dynamics and the interactions among genetic and non-genetic mechanisms. Extensive analysis of primary tumors and metastatic samples have revealed the importance of several genomic mechanisms including higher mutation rates, specific types of mutations, aneuploidy or whole genome doubling and non-genetic effects. Metastatic sites can be directly compared to primary tumor cell counterparts. At the same time, clonal dynamics shape the genomics and evolution of metastatic cancers. Clonal diversity varies by cancer type, and the tumors’ donor and recipient tissues. Still, the cancer research community has been unable to identify any common events that provide a universal predictor of “metastatic potential” which parallels findings in evolutionary ecology. Instead, invasion in cancer studies depends strongly on context, including order of events and clonal composition. The detailed studies of the behavior of a variety of human cancers promises to inform our understanding of genome level dynamics in the diversity of invasive species and provide novel insights for management.

Genomic bases underlying the adaptive radiation of core landbirds

BACKGROUND

Core landbirds undergo adaptive radiation with different ecological niches, but the genomic bases that underlie their ecological diversification remain unclear.

RESULTS

Here we used the genome-wide target enrichment sequencing of the genes related to vision, hearing, language, temperature sensation, beak shape, taste transduction, and carbohydrate, protein and fat digestion and absorption to examine the genomic bases underlying their ecological diversification. Our comparative molecular phyloecological analyses show that different core landbirds present adaptive enhancement in different aspects, and two general patterns emerge. First, all three raptorial birds (Accipitriformes, Strigiformes, and Falconiformes) show a convergent adaptive enhancement for fat digestion and absorption, while non-raptorial birds tend to exhibit a promoted capability for protein and carbohydrate digestion and absorption. Using this as a molecular marker, our results show relatively strong support for the raptorial lifestyle of the common ancestor of core landbirds, consequently suggesting a single origin of raptors, followed by two secondary losses of raptorial lifestyle within core landbirds. In addition to the dietary niche, we find at temporal niche that diurnal birds tend to exhibit an adaptive enhancement in bright-light vision, while nocturnal birds show an increased adaption in dim-light vision, in line with previous findings.

CONCLUSIONS

Our molecular phyloecological study reveals the genome-wide adaptive differentiations underlying the ecological diversification of core landbirds.

Why Do Some Vertebrates Have Microchromosomes?

With more than 70,000 living species, vertebrates have a huge impact on the field of biology and research, including karyotype evolution. One prominent aspect of many vertebrate karyotypes is the enigmatic occurrence of tiny and often cytogenetically indistinguishable microchromosomes, which possess distinctive features compared to macrochromosomes. Why certain vertebrate species carry these microchromosomes in some lineages while others do not, and how they evolve remain open questions. New studies have shown that microchromosomes exhibit certain unique characteristics of genome structure and organization, such as high gene densities, low heterochromatin levels, and high rates of recombination. Our review focuses on recent concepts to expand current knowledge on the dynamic nature of karyotype evolution in vertebrates, raising important questions regarding the evolutionary origins and ramifications of microchromosomes. We introduce the basic karyotypic features to clarify the size, shape, and morphology of macro- and microchromosomes and report their distribution across different lineages. Finally, we characterize the mechanisms of different evolutionary forces underlying the origin and evolution of microchromosomes.

Two reference-quality sea snake genomes reveal their divergent evolution of adaptive traits and venom systems

True sea snakes (Hydrophiini) are among the last and most successful clades of vertebrates that show secondary marine adaptation, exhibiting diverse phenotypic traits and lethal venom systems. To better understand their evolution, we generated the first chromosome-level genomes of two representative Hydrophiini snakes, Hydrophis cyanocinctus and H. curtus. Through comparative genomics we identified a great expansion of the underwater olfaction-related V2R gene family, consisting of more than 1,000 copies in both snakes. A series of chromosome rearrangements and genomic structural variations were recognized, including large inversions longer than 30 megabase (Mb) on sex chromosomes which potentially affect key functional genes associated with differentiated phenotypes between the two species. By integrating multiomics we found a significant loss of the major weapon for elapid predation, three-finger toxin genes, which displayed a dosage effect in H. curtus. These genetic changes may imply mechanisms that drove the divergent evolution of adaptive traits including prey preferences between the two closely related snakes. Our reference-quality sea snake genomes also enrich the repositories for addressing important issues on the evolution of marine tetrapods, and provide a resource for discovering marine-derived biological products.

Novel de Novo Genome of Cynopterus brachyotis Reveals Evolutionarily Abrupt Shifts in Gene Family Composition across Fruit Bats

Major novel physiological or phenotypic adaptations often require accompanying modifications at the genic level. Conversely, the detection of considerable contractions and/or expansions of gene families can be an indicator of fundamental but unrecognized physiological change. We sequenced a novel fruit bat genome (Cynopterus brachyotis) and adopted a comparative approach to reconstruct the evolution of fruit bats, mapping contractions and expansions of gene families along their evolutionary history. Despite a radical change in life history as compared with other bats (e.g., loss of echolocation, large size, and frugivory), fruit bats have undergone surprisingly limited change in their genic composition, perhaps apart from a potentially novel gene family expansion relating to telomere protection and longevity. In sharp contrast, within fruit bats, the new Cynopterus genome bears the signal of unusual gene loss and gene family contraction, despite its similar morphology and lifestyle to two other major fruit bat lineages. Most missing genes are regulatory, immune-related, and olfactory in nature, illustrating the diversity of genomic strategies employed by bats to contend with responses to viral infection and olfactory requirements. Our results underscore that significant fluctuations in gene family composition are not always associated with obvious examples of novel physiological and phenotypic adaptations but may often relate to less-obvious shifts in immune strategies.

Impact of Repetitive DNA Elements on Snake Genome Biology and Evolution

The distinctive biology and unique evolutionary features of snakes make them fascinating model systems to elucidate how genomes evolve and how variation at the genomic level is interlinked with phenotypic-level evolution. Similar to other eukaryotic genomes, large proportions of snake genomes contain repetitive DNA, including transposable elements (TEs) and satellite repeats. The importance of repetitive DNA and its structural and functional role in the snake genome, remain unclear. This review highlights the major types of repeats and their proportions in snake genomes, reflecting the high diversity and composition of snake repeats. We present snakes as an emerging and important model system for the study of repetitive DNA under the impact of sex and microchromosome evolution. We assemble evidence to show that certain repetitive elements in snakes are transcriptionally active and demonstrate highly dynamic lineage-specific patterns as repeat sequences. We hypothesize that particular TEs can trigger different genomic mechanisms that might contribute to driving adaptive evolution in snakes. Finally, we review emerging approaches that may be used to study the expression of repetitive elements in complex genomes, such as snakes. The specific aspects presented here will stimulate further discussion on the role of genomic repeats in shaping snake evolution.

3D genomics across the tree of life reveals condensin II as a determinant of architecture type

We investigated genome folding across the eukaryotic tree of life. We find two types of three-dimensional (3D) genome architectures at the chromosome scale. Each type appears and disappears repeatedly during eukaryotic evolution. The type of genome architecture that an organism exhibits correlates with the absence of condensin II subunits. Moreover, condensin II depletion converts the architecture of the human genome to a state resembling that seen in organisms such as fungi or mosquitoes. In this state, centromeres cluster together at nucleoli, and heterochromatin domains merge. We propose a physical model in which lengthwise compaction of chromosomes by condensin II during mitosis determines chromosome-scale genome architecture, with effects that are retained during the subsequent interphase. This mechanism likely has been conserved since the last common ancestor of all eukaryotes.

Tracking the recruitment and evolution of snake toxins using the evolutionary context provided by the Bothrops jararaca genome

Venom is a key adaptive innovation in snakes, and how nonvenom genes were co-opted to become part of the toxin arsenal is a significant evolutionary question. While this process has been investigated through the phylogenetic reconstruction of toxin sequences, evidence provided by the genomic context of toxin genes remains less explored. To investigate the process of toxin recruitment, we sequenced the genome of Bothrops jararaca, a clinically relevant pitviper. In addition to producing a road map with canonical structures of genes encoding 12 toxin families, we inferred most of the ancestral genes for their loci. We found evidence that 1) snake venom metalloproteinases (SVMPs) and phospholipases A2 (PLA2) have expanded in genomic proximity to their nonvenomous ancestors; 2) serine proteinases arose by co-opting a local gene that also gave rise to lizard gilatoxins and then expanded; 3) the bradykinin-potentiating peptides originated from a C-type natriuretic peptide gene backbone; and 4) VEGF-F was co-opted from a PGF-like gene and not from VEGF-A. We evaluated two scenarios for the original recruitment of nontoxin genes for snake venom: 1) in locus ancestral gene duplication and 2) in locus ancestral gene direct co-option. The first explains the origins of two important toxins (SVMP and PLA2), while the second explains the emergence of a greater number of venom components. Overall, our results support the idea of a locally assembled venom arsenal in which the most clinically relevant toxin families expanded through posterior gene duplications, regardless of whether they originated by duplication or gene co-option.

Urates of colubroid snakes are different from those of boids and pythonids

Uricotelic species, such as squamate reptiles, birds and insects, effectively eliminate nitrogen as uric acid in a solid form commonly called urates. Observations made over a decade suggested that the voided urates produced by colubroids (modern snake species) exhibit remarkable differences from those of boids and pythons (ancient snake species). Here, we compare the urates generated by eight captive snake species fed the same diet. Although all fresh urates were wet at the time of excretion, those produced by modern snakes dried to a powdery solid, whereas those of ancient species dried to a rock-hard mass that was tightly adherent to surfaces. Powder X-ray diffraction and infrared spectroscopy analyses performed on voided urates produced by five modern and three ancient snakes confirmed their underlying chemical and structural differences. Urates excreted by ancient snakes were amorphous uric acid, whereas urates from modern snakes consisted primarily of ammonium acid urate, with some uric acid dihydrate. These compositional differences indicate that snakes have more than one mechanism to manage nitrogenous waste. Why different species use different nitrogen-handling pathways is not yet known, but the answer might be related to key differences in metabolism, physiology or, in the case of ancient snakes, the potential use of urates in social communication.

Venomous Snakes Reveal Ecological and Phylogenetic Factors Influencing Variation in Gut and Oral Microbiomes

The gastrointestinal tract (GIT) of vertebrates contains a series of organs beginning with the mouth and ending with the anus or cloacal opening. Each organ represents a unique environment for resident microorganisms. Due to their simple digestive anatomy, snakes are good models for studying microbiome variation along the GIT. Cloacal sampling captures the majority of the microbial diversity found in the GIT of snakes—yet little is known about the oral microbiota of snakes. Most research on the snake mouth and gut microbiota are limited to studies of a single species or captive-bred individuals. It therefore remains unclear how a host’s life history, diet, or evolutionary history correlate with differences in the microbial composition within the mouths and guts of wild snakes. We sampled the mouth and gut microbial communities from three species of Asian venomous snakes and utilized 16S rRNA microbial inventories to test if host phylogenetic and ecological differences correlate with distinct microbial compositions within the two body sites. These species occupy three disparate habitat types: marine, semi-arboreal, and arboreal, our results suggest that the diversity of snake mouth and gut microbial communities correlate with differences in both host ecology and phylogeny.

New Ther1-derived SINE Squam3 in scaled reptiles

BACKGROUND

SINEs comprise a significant part of animal genomes and are used to study the evolution of diverse taxa. Despite significant advances in SINE studies in vertebrates and higher eukaryotes in general, their own evolution is poorly understood.

RESULTS

We have discovered and described in detail a new Squam3 SINE specific for scaled reptiles (Squamata). The subfamilies of this SINE demonstrate different distribution in the genomes of squamates, which together with the data on similar SINEs in the tuatara allowed us to propose a scenario of their evolution in the context of reptilian evolution.

CONCLUSIONS

Ancestral SINEs preserved in small numbers in most genomes can give rise to taxa-specific SINE families. Analysis of this aspect of SINEs can shed light on the history and mechanisms of SINE variation in reptilian genomes.

Snake Recombination Landscapes Are Concentrated in Functional Regions despite PRDM9

Meiotic recombination in vertebrates is concentrated in hotspots throughout the genome. The location and stability of hotspots have been linked to the presence or absence of PRDM9, leading to two primary models for hotspot evolution derived from mammals and birds. Species with PRDM9-directed recombination have rapid turnover of hotspots concentrated in intergenic regions (i.e., mammals), whereas hotspots in species lacking PRDM9 are concentrated in functional regions and have greater stability over time (i.e., birds). Snakes possess PRDM9, yet virtually nothing is known about snake recombination. Here, we examine the recombination landscape and test hypotheses about the roles of PRDM9 in rattlesnakes. We find substantial variation in recombination rate within and among snake chromosomes, and positive correlations between recombination rate and gene density, GC content, and genetic diversity. Like mammals, snakes appear to have a functional and active PRDM9, but rather than being directed away from genes, snake hotspots are concentrated in promoters and functional regions-a pattern previously associated only with species that lack a functional PRDM9. Snakes therefore provide a unique example of recombination landscapes in which PRDM9 is functional, yet recombination hotspots are associated with functional genic regions-a combination of features that defy existing paradigms for recombination landscapes in vertebrates. Our findings also provide evidence that high recombination rates are a shared feature of vertebrate microchromosomes. Our results challenge previous assumptions about the adaptive role of PRDM9 and highlight the diversity of recombination landscape features among vertebrate lineages.