Search

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.

Hybrid zones provide valuable opportunities to understand the genomic mechanisms that promote speciation by providing insight into factors involved in intermediate stages of speciation. Here, we investigate introgression in a hybrid zone between two rattlesnake species (Crotalus viridis and Crotalus oreganus concolor) that have undergone historical allopatric divergence and recent range expansion and secondary contact. We use Bayesian genomic cline models to characterize genomic patterns of introgression between these lineages and identify loci potentially subject to selection in hybrids. We find evidence for a large number of genomic regions with biased ancestry that deviate from the genomic background in hybrids (i.e., excess ancestry loci), which tend to be associated with genomic regions with higher recombination rates. We also identify suites of excess ancestry loci that show highly correlated allele frequencies (including conspecific and heterospecific combinations) across physically unlinked genomic regions in hybrids. Our findings provide evidence for multiple multilocus evolutionary processes impacting hybrid fitness in this system.

Despite the extensive body of research on snake venom, many facets of snake venom systems, such as the physiology and regulation of the venom gland itself, remain virtually unstudied. Here, we use time series gene expression analyses of the rattlesnake venom gland in comparison with several non-venom tissues to characterize physiological and cellular processes associated with venom production and to highlight key distinctions of venom gland cellular and physiological function. We find consistent evidence for activation of stress response pathways in the venom gland, suggesting that mitigation of cellular stress is a crucial component of venom production. Additionally, we demonstrate evidence for an unappreciated degree of cellular and secretory activity in the steady state venom gland relative to other secretory tissues and identify vacuolar ATPases as the likely mechanisms driving acidification of the venom gland lumen during venom production and storage.

Organisms use several strategies to mitigate mitochondrial stress, including the activation of the mitochondrial unfolded protein response (UPR

Understanding how regulatory mechanisms evolve is critical for understanding the processes that give rise to novel phenotypes. Snake venom systems represent a valuable and tractable model for testing hypotheses related to the evolution of novel regulatory networks, yet the regulatory mechanisms underlying venom production remain poorly understood. Here, we use functional genomics approaches to investigate venom regulatory architecture in the prairie rattlesnake and identify cis-regulatory sequences (enhancers and promoters), trans-regulatory transcription factors, and integrated signaling cascades involved in the regulation of snake venom genes. We find evidence that two conserved vertebrate pathways, the extracellular signal-regulated kinase and unfolded protein response pathways, were co-opted to regulate snake venom. In one large venom gene family (snake venom serine proteases), this co-option was likely facilitated by the activity of transposable elements. Patterns of snake venom gene enhancer conservation, in some cases spanning 50 million yr of lineage divergence, highlight early origins and subsequent lineage-specific adaptations that have accompanied the evolution of venom regulatory architecture. We also identify features of chromatin structure involved in venom regulation, including topologically associated domains and CTCF loops that underscore the potential importance of novel chromatin structure to coevolve when duplicated genes evolve new regulatory control. Our findings provide a model for understanding how novel regulatory systems may evolve through a combination of genomic processes, including tandem duplication of genes and regulatory sequences, cis-regulatory sequence seeding by transposable elements, and diverse transcriptional regulatory proteins controlled by a co-opted regulatory cascade.

Over the past decade, the increasing availability and decreasing cost of genome sequencing technologies have led to significant advancements in our understanding of the evolution, variation and function of reptile venoms. This includes multiple annotated genomes for venomous reptile species as well as countless public data sets of venom gene sequences, gene expression data, and data availability from a growing diversity of species. These genomic resources have provided new insight into reptile venoms, and the continuing development of additional genomic resources and related data sets holds promise for continued advancement of our understanding of reptile venom systems. Here, we present an overview of currently available genomic resources for the study of reptile venoms, discuss genomic approaches used to study reptile venoms, and highlight recent advances in our understanding of how venom evolves and functions in reptiles based on genomic data. We also discuss ongoing progress and future work that will be important for closing major gaps in our current understanding of the genomic basis of venom variation, evolution and regulation.

High-quality genomic resources facilitate population-level and species-level comparisons to answer questions about behavioral ecology, morphological and physiological adaptations, as well as the evolution of genomic architecture. Squamate reptiles (lizards and snakes) are particularly diverse in characteristics that have intrigued evolutionary biologists, but high-quality genomic resources for squamates are relatively sparse. 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. We present a high-quality chromosome-level reference genome assembly, SceUnd1.0, (utilizing 10X Genomics Chromium, HiC, and PacBio data) and tissue/developmental stage transcriptomes for the Eastern Fence Lizard, Sceloporus undulatus. We performed synteny analysis with other available squamate chromosome-level assemblies to identify broad patterns of chromosome evolution including the fusion of micro- and macrochromosomes in S. undulatus. Using this new S. undulatus genome assembly we conducted reference-based assemblies for 34 other Sceloporus species to improve draft nuclear genomes assemblies from 1% coverage to 43% coverage on average. Across these species, typically >90% of reads mapped for species within 20 million years divergence from S. undulatus, this dropped to 75% reads mapped for species at 35 million years divergence. Finally we use RNAseq and whole genome resequencing data to compare the three assemblies as references, each representing an increased level of sequencing, cost and assembly efforts: Supernova Assembly with data from10X Genomics Chromium library; HiRise Assembly that added data from HiC library; and PBJelly Assembly that added data from PacBio sequencing. We found that the Supernova Assembly contained the full genome and was a suitable reference for RNAseq, but the chromosome-level scaffolds provided by the addition of the HiC data allowed the reference to be used for other whole genome analysis, including synteny and whole genome association mapping analyses. The addition of PacBio data provided negligible gains. Overall, these new genomic resources provide valuable tools for advanced molecular analysis of an organism that has become a model in physiology and evolutionary ecology.

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.

Invasive Burmese pythons (Python bivittatus Kuhl, 1820) have introduced a lung parasite, Raillietiella orientalis, (Hett, 1915) from the python's native range in Southeast Asia to its introduced range in Florida, where parasite spillover from pythons to two families and eight genera of native snakes has occurred. Because these novel host species present a diversity of ecological and morphological traits, and because these parasites attach to their hosts with hooks located on their cephalothorax, we predicted that R. orientalis would exhibit substantial, host-associated phenotypic plasticity in cephalothorax shape. Indeed, geometric morphometric analyses of 39 parasites from five host species revealed significant variation among host taxa in R. orientalis cephalothorax shape. We observed differences associated with host ecology, where parasites from semi-aquatic and aquatic snakes exhibited the greatest morphological similarity. Morphological analyses of R. orientalis recovered from invasive pythons, native pit vipers, and terrestrial snakes each revealed distinct shapes. Our results suggest R. orientalis can exhibit significant differences in morphology based upon host species infected, and this plasticity may facilitate infection with this non-native parasite in a wide array of novel squamate host species.