The genetic regulatory architecture and epigenomic basis for age-related changes in rattlesnake venom

The gene regulatory mechanisms shaping the heterogeneity of venom production in the Cape coral snake

BACKGROUND

Venoms and their associated glands and delivery structures have evolved numerous times among animals. Within these venom systems, the molecular, cellular, and morphological components interact and co-evolve to generate distinct, venom phenotypes that are increasingly recognized as models for studying adaptive evolution. However, toxins are often unevenly distributed across venom-producing tissues in patterns that are not necessarily adaptive but instead likely result from constraints associated with protein secretion.

RESULTS

We generate a high-quality draft genome of the Cape coral snake (Aspidelaps lubricus) and combine analyses of venom gland single-cell RNA-seq data with spatial venom gland in situ toxin distributions. Our results reveal that while different toxin families are produced by distinct populations of cells, toxin expression is fine-tuned by regulatory modules that result in further specialization of toxin production within each cell population. We also find that the evolution of regulatory elements closely mirrors the evolution of their associated toxin genes, resulting in spatial association of closely related and functionally similar toxins in the venom gland. While this compartmentalization is non-adaptive, the modularity of the underlying regulatory network likely facilitated the repeated evolution of defensive venom in spitting cobras.

CONCLUSIONS

Our results provide new insight into the variability of toxin regulation across snakes, reveal the molecular mechanisms underlying the heterogeneous toxin production in snake venom glands, and provide an example of how constraints can result in non-adaptive character states that appear to be adaptive, which may nevertheless facilitate evolutionary innovation and novelty.

Unique physiological and regulatory activity drives divergent toxin and non-toxin gene expression in rattlesnake accessory venom glands

Understanding the mechanisms by which organs and tissues evolve new physiological functions is central to understanding the evolution of novelty. This is particularly interesting in the context of related tissues that evolve specialized, yet complementary, functions. Snake venom glands are an attractive system to test hypotheses related to the evolution and specialization of novel physiological function, as these modified salivary glands have evolved over ∼60 MY to synthesize and store venom. Front-fanged venomous snakes (elapids and viperids) possess two types of venom glands: the main and accessory glands. The larger main gland produces greater quantities of venom toxins and has been studied extensively, while the smaller accessory gland has received less attention. Here, we explore gene expression differences between main and accessory venom glands across three rattlesnake species (Crotalus cerberus, C. oreganus concolor and C. viridis). Our findings indicate that accessory glands express most venom genes at significantly lower levels than the main gland, with a few exceptions that may represent biologically relevant contributions of accessory glands to venom. The two glands also exhibit distinct trans-regulatory environments that we link to key differences in their underlying physiology and secretory roles. Our results further suggest that two signaling pathways that regulate venom, the unfolded protein response (UPR) and extracellular signal-regulated kinase (ERK), show significantly lower activation in the accessory gland. These findings provide insight into the physiological and functional diversification of snake venom systems, highlighting how distinct glandular systems have evolved contrasting and complementary roles driven by distinct physiological and regulatory mechanisms.

Life history and chromosome organization determine chemoreceptor gene expression in rattlesnakes

Predatory species who hunt for their prey rely on a suite of integrated characters, including sensory traits that are also used for nonpredatory behaviors. Linking the evolution of sensory traits to specific selection pressures therefore requires a deep understanding of the underlying genetics and molecular mechanisms producing these complex phenotypes. However, this relationship remains poorly understood for complex sensory systems that consist of proteins encoded by large gene families. The chemosensory repertoire of rattlesnakes includes hundreds of type-2 vomeronasal receptors and olfactory receptors, representing the two largest gene families found in the genome. To investigate the biological importance of this chemoreceptor diversity, we assessed gene expression in the eastern diamondback rattlesnake (Crotalus adamanteus) and identified sex- and age-biased genes. We found type-2 vomeronasal receptor expression in the vomeronasal epithelium was limited to juvenile snakes, suggesting the sensory programming of this tissue may be correlated with early life development. In the olfactory epithelium, we found subtle expression biases that were more indicative of life history rather than development. We also found transcriptional evidence for dosage compensation of sex-linked genes and trait integration in the expression of transcription factors. We overlay our molecular characterizations in Crotalus adamanteus onto updated olfactory receptor and type-2 vomeronasal receptor phylogenies, providing a genetic road map for future research on these receptors. Finally, we investigated the deeper macroevolutionary context of the most highly expressed type-2 vomeronasal receptor gene spanning the rise of tetrapods and estimated the strength of positive selection for individual amino acid residues in the predicted protein structure. We hypothesize that this gene may have evolved as a conserved signaling subunit to ensure consistent G-protein coupled receptor functionality, potentially relaxing signaling constraints on other type-2 vomeronasal receptor paralogs and promoting ligand binding specificity.

A Segregating Structural Variant Defines Novel Venom Phenotypes in the Eastern Diamondback Rattlesnake

Of all mutational mechanisms contributing to phenotypic variation, structural variants are both among the most capable of causing major effects as well as the most technically challenging to identify. Intraspecific variation in snake venoms is widely reported, and one of the most dramatic patterns described is the parallel evolution of streamlined neurotoxic rattlesnake venoms from hemorrhagic ancestors by means of deletion of snake venom metalloproteinase (SVMP) toxins and recruitment of neurotoxic dimeric phospholipase A2 (PLA2) toxins. While generating a haplotype-resolved, chromosome-level genome assembly for the eastern diamondback rattlesnake (Crotalus adamanteus), we discovered that our genome animal was heterozygous for a ∼225 Kb deletion containing six SVMP genes, paralleling one of the two steps involved in the origin of neurotoxic rattlesnake venoms. Range-wide population-genomic analysis revealed that, although this deletion is rare overall, it is the dominant homozygous genotype near the northwestern periphery of the species' range, where this species is vulnerable to extirpation. Although major SVMP deletions have been described in at least five other rattlesnake species, C. adamanteus is unique in not additionally gaining neurotoxic PLA2s. Previous work established a superficially complementary north-south gradient in myotoxin (MYO) expression based on copy number variation with high expression in the north and low in the south, yet we found that the SVMP and MYO genotypes vary independently, giving rise to an array of diverse, novel venom phenotypes across the range. Structural variation, therefore, forms the basis for the major axes of geographic venom variation for C. adamanteus.

Dietary Breadth Predicts Toxin Expression Complexity in the Venoms of North American Gartersnakes

Selection on heritable phenotypic variation has played a prominent role in shaping the remarkable adaptations found across the Tree of Life. Complex ecological traits, such as snake venoms, are thought to be the products of selection because they directly link to fitness and survival. Snake venom increases the efficiency of prey capture and processing and is thus likely under intense selection. While many studies of snake venom have investigated the relationship between venom and diet, they have primarily focused on medically relevant front-fanged snakes. However, recent work has suggested that many non-front fanged snakes also rely on venom for subduing prey, despite having reduced toxicity toward humans. Here, we set out to uncover variation in toxin-producing genes, along with the ecological and evolutionary pressures impacting snake venom characteristics in the North American gartersnakes (Squamata: Natricidae: Thamnophis), a model group of non-front-fanged snakes. We annotated and analyzed Duvernoy's venom gland transcriptomes from 16 species representing all the major lineages within Thamnophis. We then generated measures of complexity of both toxins and dietary breadth. We found strong correlations between the complexity of toxin gene expression and phylogenetic diversity of diet, but no relationship between the complexity of the genetic makeup of the transcriptomes (allelic or sequence variation) and diet complexities. We also found phylogenetic signal associated with venom complexity, suggesting some influence of ancestry on venom characteristics. We suggest that, in non-front-fanged snakes, expression of toxins rather than sequence complexity is under strong selection by dietary diversity. These findings contradict similar studies from front-fanged snakes where increased transcriptomic complexity varies positively with dietary diversity, exposing a potential novel relationship between a complex phenotype-toxin expression-and its selective pressures-diet.

Thamnophis sirtalis and their toxic relationship: Testing for intraspecific venom variation in Common Garter Snakes

Intraspecific phenotypic variation can be used as a window into the ecological differences among individuals of a species and lead to a better understanding of adaptive evolution. Adaptive traits, such as venom, that play an important ecological role for a species are useful models for understanding the sources of intraspecific variation. Intraspecific studies on front-fanged venomous snakes have offered deeper insights into the diverse mechanisms and adaptations that support the effectiveness of venom across species. Despite the extensive research on front-fanged venomous snakes, rear-fanged snakes, representing two-thirds of all snake species, have been largely overlooked. To test for sex and age-based intraspecific venom variation, we sequenced the messenger RNA from the Duvernoy's gland of 9 male and 10 female Common Garter Snakes, Thamnophis sirtalis, of different sizes from a single location. Our data represent the most venom gland transcriptomes of any venomous snake species from a single location and represent the first Duvernoy's venom gland transcriptomes for Thamnophis sirtalis. We found four toxin families dominate the Thamnophis sirtalis transcriptome: Snake Venom Metalloproteinases (SVMPs), Three-finger toxins (3FTxs), Cysteine-Rich Secretory Proteins (CRISPs), and C-type lectins (CTLs). Thamnophis sirtalis exhibits a unique balance in toxin expression, with approximately 30% each of neurotoxic (3FTx-dominated) and enzymatic (SVMP-dominated) components. No other published RFS Duvernoy's gland transcriptome displays this ratio, rather they are dominated by one or the other. Additionally, venom expression varies with sex and size, with differences in toxin gene expression between males and females as they grow. Our study provides new insights on venom composition in a RFS species and highlights the amount of intraspecific variation possible among individuals from a single population.

VenomCap: An exon-capture probe set for the targeted sequencing of snake venom genes

Snake venoms are complex mixtures of toxic proteins that hold significant medical, pharmacological and evolutionary interest. To better understand the genetic diversity underlying snake venoms, we developed VenomCap, a novel exon-capture probe set targeting toxin-coding genes from a wide range of elapid snakes, with a particular focus on the ecologically diverse and medically important subfamily Hydrophiinae. We tested the capture success of VenomCap across 24 species, representing all major elapid lineages. We included snake phylogenomic probes in the VenomCap capture set, allowing us to compare capture performance between venom and phylogenomic loci and to infer elapid phylogenetic relationships. We demonstrated VenomCap's ability to recover exons from ~1500 target markers, representing a total of 24 known venom gene families, which includes the dominant gene families found in elapid venoms. We find that VenomCap's capture results are robust across all elapids sampled, and especially among hydrophiines, with respect to measures of target capture success (target loci matched, sensitivity, specificity and missing data). As a cost-effective and efficient alternative to full genome sequencing, VenomCap can dramatically accelerate the sequencing and analysis of venom gene families. Overall, our tool offers a model for genomic studies on snake venom gene diversity and evolution that can be expanded for comprehensive comparisons across the other families of venomous snakes.

Current Technologies in Snake Venom Analysis and Applications

This comprehensive review explores the cutting-edge advancements in snake venom research, focusing on the integration of proteomics, genomics, transcriptomics, and bioinformatics. Highlighting the transformative impact of these technologies, the review delves into the genetic and ecological factors driving venom evolution, the complex molecular composition of venoms, and the regulatory mechanisms underlying toxin production. The application of synthetic biology and multi-omics approaches, collectively known as venomics, has revolutionized the field, providing deeper insights into venom function and its therapeutic potential. Despite significant progress, challenges such as the functional characterization of toxins and the development of cost-effective antivenoms remain. This review also discusses the future directions of venom research, emphasizing the need for interdisciplinary collaborations and new technologies (mRNAs, cryo-electron microscopy for structural determinations of toxin complexes, synthetic biology, and other technologies) to fully harness the biomedical potential of venoms and toxins from snakes and other animals.

Tiny but Mighty: Vipera ammodytes meridionalis (Eastern Long-Nosed Viper) Ontogenetic Venom Variations in Procoagulant Potency and the Impact on Antivenom Efficacies

The Eastern Long-Nosed Viper (Vipera ammodytes meridionalis) is considered one of the most venomous snakes in Europe. However, it is unknown whether ontogenetic variation in venom effects occurs in this subspecies and how this may impact antivenom efficacy. In this study, we compared the procoagulant activities of V. a. meridionalis venom on human plasma between neonate and adult venom phenotypes. We also examined the efficacy of three antivenoms-Viperfav, ViperaTAb, and Inoserp Europe-across our neonate and adult venom samples. While both neonate and adult V. a. meridionalis venoms produced procoagulant effects, the effects produced by neonate venom were more potent. Consistent with this, neonate venom was a stronger activator of blood-clotting zymogens, converting them into their active forms, with a rank order of Factor X >> Factor VII > Factor XII. Conversely, the less potent adult venom had a rank order of FXII marginally more activated than Factor VII, and both much more so than Factor X. This adds to the growing body of evidence that activation of factors besides FII (prothrombin) and FX are significant variables in reptile venom-induced coagulopathy. Although all three examined antivenoms displayed effective neutralization of both neonate and adult V. a. meridionalis venoms, they generally showed higher efficacy on adult venom than on neonate venom. The ranking of antivenom efficacy against neonate venom, from the most effective to the least effective, were Viperfav, Inoserp Europe, ViperaTAb; for adult venom, the ranking was Inoserp Europe, Viperfav, ViperaTAb. Our data reveal ontogenetic variation in V. a meridionalis, but this difference may not be of clinical concern as antivenom was effective at neutralizing both adult and neonate venom phenotypes. Regardless, our results highlight a previously undocumented ontogenetic shift, likely driven by the documented difference in prey preference observed for this species across age classes.

Where the "ruber" Meets the Road: Using the Genome of the Red Diamond Rattlesnake to Unravel the Evolutionary Processes Driving Venom Evolution

Understanding the proximate and ultimate causes of phenotypic variation is fundamental in evolutionary research, as such variation provides the substrate for selection to act upon. Although trait variation can arise due to selection, the importance of neutral processes is sometimes understudied. We presented the first reference-quality genome of the Red Diamond Rattlesnake (Crotalus ruber) and used range-wide 'omic data to estimate the degree to which neutral and adaptive evolutionary processes shaped venom evolution. We characterized population structure and found substantial genetic differentiation across two populations, each with distinct demographic histories. We identified significant differentiation in venom expression across age classes with substantially reduced but discernible differentiation across populations. We then used conditional redundancy analysis to test whether venom expression variation was best predicted by neutral divergence patterns or geographically variable (a)biotic factors. Snake size was the most significant predictor of venom variation, with environment, prey availability, and neutral sequence variation also identified as significant factors, though to a lesser degree. By directly including neutrality in the model, our results confidently highlight the predominant, yet not singular, role of life history in shaping venom evolution.

Web of venom: exploration of big data resources in animal toxin research

Research on animal venoms and their components spans multiple disciplines, including biology, biochemistry, bioinformatics, pharmacology, medicine, and more. Manipulating and analyzing the diverse array of data required for venom research can be challenging, and relevant tools and resources are often dispersed across different online platforms, making them less accessible to nonexperts. In this article, we address the multifaceted needs of the scientific community involved in venom and toxin-related research by identifying and discussing web resources, databases, and tools commonly used in this field. We have compiled these resources into a comprehensive table available on the VenomZone website (https://venomzone.expasy.org/10897). Furthermore, we highlight the challenges currently faced by researchers in accessing and using these resources and emphasize the importance of community-driven interdisciplinary approaches. We conclude by underscoring the significance of enhancing standards, promoting interoperability, and encouraging data and method sharing within the venom research community.