From molecules to macroevolution: Venom as a model system for evolutionary biology across levels of life

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.

Molecular diversity assessed by MALDI mass spectrometry of two scorpion species venom from two different locations in Morocco

Scorpion venom is a cocktail of molecules whose composition is remarkably plastic, controlled by several factors. The Moroccan scorpion fauna is characterized by its richness and high rate of endemism and the venom molecular variability of many species is not yet well characterized. The aim of the present study was to highlight the molecular variability of the venom composition of Androctonus amoreuxi and Buthacus stockmanni (endemic species), both belonging to the Buthidae family, collected from two Moroccan regions, Zagora and Tan-tan. Characterization of the molecular mass fingerprints (MFPs) of each specimen was performed by Matrix-assisted laser desorption ionization-mass spectrometry (MALDI-MS) using a sandwich (Sand) and a dried-droplet (DD) sample preparation and dilutions. Considering these two methods, a total of 828 ion signals were detected, and Sand method produced more adducts (56%) than DD (44%). We observed interspecific variations in the venom composition between these two species showing they share 235 ion signals, while 226 and 367 are specific for these two species, respectively. Moreover, B. stockmanni specimens showed a clear difference in their MFPs between the two geographical areas studied, suggesting intraspecific variations. Moreover, specimens from each population also show an intraspecific variability. In addition, for the same individual, a variation in the venom composition was also recorded depending on the milking frequency. Our results confirmed the presence of characteristic components in each extracted venom sample. In conclusion, MFPs assessed by MALDI-MS represent a fast, non-supervised, sensitive, reliable and cost-efficient approach for taxonomic identification and molecular variability characterization. This study undoubtedly represents a step forward for understanding the scorpion venom plasticity, intra/inter variations, and their temporal and geographical variability.

Comprehensive isoform-level analysis reveals the contribution of alternative isoforms to venom evolution and repertoire diversity

Animal venom systems have emerged as valuable models for investigating how novel polygenic phenotypes may arise from gene evolution by varying molecular mechanisms. However, a significant portion of venom genes produce alternative mRNA isoforms that have not been extensively characterized, hindering a comprehensive understanding of venom biology. In this study, we present a full-length isoform-level profiling workflow integrating multiple RNA sequencing technologies, allowing us to reconstruct a high-resolution transcriptome landscape of venom genes in the parasitoid wasp Pteromalus puparum Our findings demonstrate that more than half of the venom genes generate multiple isoforms within the venom gland. Through mass spectrometry analysis, we confirm that alternative splicing contributes to the diversity of venom proteins, acting as a mechanism for expanding the venom repertoire. Notably, we identified seven venom genes that exhibit distinct isoform usages between the venom gland and other tissues. Furthermore, evolutionary analyses of venom serpin3 and orcokinin further reveal that the co-option of an ancient isoform and a newly evolved isoform, respectively, contributes to venom recruitment, providing valuable insights into the genetic mechanisms driving venom evolution in parasitoid wasps. Overall, our study presents a comprehensive investigation of venom genes at the isoform level, significantly advancing our understanding of alternative isoforms in venom diversity and evolution and setting the stage for further in-depth research on venoms.

Diversity Begets Diversity When Diet Drives Snake Venom Evolution, but Evenness Rather Than Richness Is What Counts

Snake venoms are primarily used to subjugate prey, and consequently, their evolution has been shown to be predominantly driven by diet-related selection pressure. Venoms tend to be more lethal to prey than non-prey species (except in cases of toxin resistance), prey-specific toxins have been identified, and preliminary work has demonstrated an association between the diversity of diet classes and that of toxicological activities of whole venom. However, venoms are complex mixtures of many toxins, and it remains unclear how toxin diversity is driven by diet. Prey-specific toxins do not encompass the molecular diversity of venoms, and whole venom effects could be driven by one, few, or all components, so the link between diet and venom diversity remains minimally understood. Here, we collated a database of venom composition and diet records and used a combination of phylogenetic comparative methods and two quantitative diversity indices to investigate whether and how diet diversity relates to the toxin diversity of snake venoms. We reveal that venom diversity is negatively related to diet diversity using Shannon's index but positively related using Simpson's index. Since Shannon's index predominantly considers the number of prey/toxins, whereas Simpson's index more strongly reflects evenness, we provide insights into how the diet-venom diversity link is driven. Specifically, species with low diet diversity tend to have venoms dominated by a few abundant (possibly specialised) toxin families, whereas species with diverse diets tend to 'hedge their bets' by having venoms with a more even composition of different toxin classes.

Toxic Habits: An Analysis of General Trends and Biases in Snake Venom Research

Biases in snake venom research have been partially identified but seldomly quantified. Using the Google Scholar web search engine, we collected a total of 267 articles published between 1964 and 2021, and reviewed them to assess the main trends in this field of study. We developed a 4-category classification of the harmful potential of each of the 298 snake species retrieved from the analysed publications, and tested whether taxonomy, realm of origin, and/or assigned hazard category could affect how often each of them appeared in the articles considered. Overall, viperids were significantly more represented than any other snake taxon retrieved. The Neotropics were the most represented biogeographic realm for number of studied species, whereas information about the country of origin of the analysed specimens was often incomplete. The vast majority of the publications focused on snake venom characterisation, whereas more ecology-related topics were rarely considered. Hazard category and biogeographic realm of origin of each species had a significant effect on the number of articles dedicated to it, suggesting that a snake's harmful potential and place of origin influence its popularity in venom studies. Our analysis showed an overall positive trend in the number of snake venom studies published yearly, but also underlined severe neglect of snake families of supposedly minor medical relevance (e.g., Atractaspididae), underrepresentation of some of the areas most impacted by snakebite (i.e., Indomalayan and Afrotropic realms), and limited interest in the ecological and functional context of snake venom.

The roles of balancing selection and recombination in the evolution of rattlesnake venom

The origin of snake venom involved duplication and recruitment of non-venom genes into venom systems. Several studies have predicted that directional positive selection has governed this process. Venom composition varies substantially across snake species and venom phenotypes are locally adapted to prey, leading to coevolutionary interactions between predator and prey. Venom origins and contemporary snake venom evolution may therefore be driven by fundamentally different selection regimes, yet investigations of population-level patterns of selection have been limited. Here, we use whole-genome data from 68 rattlesnakes to test hypotheses about the factors that drive genomic diversity and differentiation in major venom gene regions. We show that selection has resulted in long-term maintenance of genetic diversity within and between species in multiple venom gene families. Our findings are inconsistent with a dominant role of directional positive selection and instead support a role of long-term balancing selection in shaping venom evolution. We also detect rapid decay of linkage disequilibrium due to high recombination rates in venom regions, suggesting that venom genes have reduced selective interference with nearby loci, including other venom paralogues. Our results provide an example of long-term balancing selection that drives trans-species polymorphism and help to explain how snake venom keeps pace with prey resistance.

Modern venomics-Current insights, novel methods, and future perspectives in biological and applied animal venom research

Venoms have evolved >100 times in all major animal groups, and their components, known as toxins, have been fine-tuned over millions of years into highly effective biochemical weapons. There are many outstanding questions on the evolution of toxin arsenals, such as how venom genes originate, how venom contributes to the fitness of venomous species, and which modifications at the genomic, transcriptomic, and protein level drive their evolution. These questions have received particularly little attention outside of snakes, cone snails, spiders, and scorpions. Venom compounds have further become a source of inspiration for translational research using their diverse bioactivities for various applications. We highlight here recent advances and new strategies in modern venomics and discuss how recent technological innovations and multi-omic methods dramatically improve research on venomous animals. The study of genomes and their modifications through CRISPR and knockdown technologies will increase our understanding of how toxins evolve and which functions they have in the different ontogenetic stages during the development of venomous animals. Mass spectrometry imaging combined with spatial transcriptomics, in situ hybridization techniques, and modern computer tomography gives us further insights into the spatial distribution of toxins in the venom system and the function of the venom apparatus. All these evolutionary and biological insights contribute to more efficiently identify venom compounds, which can then be synthesized or produced in adapted expression systems to test their bioactivity. Finally, we critically discuss recent agrochemical, pharmaceutical, therapeutic, and diagnostic (so-called translational) aspects of venoms from which humans benefit.

Venom gene sequence diversity and expression jointly shape diet adaptation in pitvipers

Understanding the joint roles of protein sequence variation and differential expression during adaptive evolution is a fundamental, yet largely unrealized goal of evolutionary biology. Here, we use phylogenetic path analysis to analyze a comprehensive venom-gland transcriptome dataset spanning three genera of pitvipers to identify the functional genetic basis of a key adaptation (venom complexity) linked to diet breadth (DB). The analysis of gene-family-specific patterns reveals that, for genes encoding two of the most important venom proteins (snake venom metalloproteases and snake venom serine proteases), there are direct, positive relationships between sequence diversity (SD), expression diversity (ED), and increased DB. Further analysis of gene-family diversification for these proteins showed no constraint on how individual lineages achieved toxin gene SD in terms of the patterns of paralog diversification. In contrast, another major venom protein family (PLA₂s) showed no relationship between venom molecular diversity and DB. Additional analyses suggest that other molecular mechanisms—such as higher absolute levels of expression—are responsible for diet adaptation involving these venom proteins. Broadly, our findings argue that functional diversity generated through sequence and expression variations jointly determine adaptation in the key components of pitviper venoms, which mediate complex molecular interactions between the snakes and their prey.

Old World Vipers-A Review about Snake Venom Proteomics of Viperinae and Their Variations

Fine-tuned by millions of years of evolution, snake venoms have frightened but also fascinated humanity and nowadays they constitute potential resources for drug development, therapeutics and antivenoms. The continuous progress of mass spectrometry techniques and latest advances in proteomics workflows enabled toxinologists to decipher venoms by modern omics technologies, so-called ‘venomics’. A tremendous upsurge reporting on snake venom proteomes could be observed. Within this review we focus on the highly venomous and widely distributed subfamily of Viperinae (Serpentes: Viperidae). A detailed public literature database search was performed (2003–2020) and we extensively reviewed all compositional venom studies of the so-called Old-World Vipers. In total, 54 studies resulted in 89 venom proteomes. The Viperinae venoms are dominated by four major, four secondary, six minor and several rare toxin families and peptides, respectively. The multitude of different venomics approaches complicates the comparison of venom composition datasets and therefore we differentiated between non-quantitative and three groups of quantitative workflows. The resulting direct comparisons within these groups show remarkable differences on the intra- and interspecies level across genera with a focus on regional differences. In summary, the present compilation is the first comprehensive up-to-date database on Viperinae venom proteomes and differentiating between analytical methods and workflows.

Phylogenetically diverse diets favor more complex venoms in North American pitvipers

The role of natural selection in the evolution of trait complexity can be characterized by testing hypothesized links between complex forms and their functions across species. Predatory venoms are composed of multiple proteins that collectively function to incapacitate prey. Venom complexity fluctuates over evolutionary timescales, with apparent increases and decreases in complexity, and yet the causes of this variation are unclear. We tested alternative hypotheses linking venom complexity and ecological sources of selection from diet in the largest clade of front-fanged venomous snakes in North America: the rattlesnakes, copperheads, cantils, and cottonmouths. We generated independent transcriptomic and proteomic measures of venom complexity and collated several natural history studies to quantify dietary variation. We then constructed genome-scale phylogenies for these snakes for comparative analyses. Strikingly, prey phylogenetic diversity was more strongly correlated to venom complexity than was overall prey species diversity, specifically implicating prey species' divergence, rather than the number of lineages alone, in the evolution of complexity. Prey phylogenetic diversity further predicted transcriptomic complexity of three of the four largest gene families in viper venom, showing that complexity evolution is a concerted response among many independent gene families. We suggest that the phylogenetic diversity of prey measures functionally relevant divergence in the targets of venom, a claim supported by sequence diversity in the coagulation cascade targets of venom. Our results support the general concept that the diversity of species in an ecological community is more important than their overall number in determining evolutionary patterns in predator trait complexity.

Special Issue: Evolutionary Ecology of Venom

This Special Issue of Toxins aims to increase the profile and understanding of how ecology shapes the evolution of venom systems, and also how venom influences the ecological attributes of and interactions among species [...].

Diet Diversity in Carnivorous Terebrid Snails Is Tied to the Presence and Absence of a Venom Gland

Predator-prey interactions are thought to play a driving role in animal evolution, especially for groups that have developed venom as their predatory strategy. However, how the diet of venomous animals influences the composition of venom arsenals remains uncertain. Two prevailing hypotheses to explain the relationship between diet and venom composition focus on prey preference and the types of compounds in venom, and a positive correlation between dietary breadth and the number of compounds in venom. Here, we examined venom complexity, phylogenetic relationship, collection depth, and biogeography of the Terebridae (auger snails) to determine if repeated innovations in terebrid foregut anatomy and venom composition correspond to diet variation. We performed the first molecular study of the diet of terebrid marine snails by metabarcoding the gut content of 71 terebrid specimens from 17 species. Our results suggest that the presence or absence of a venom gland is strongly correlated with dietary breadth. Specifically, terebrid species without a venom gland displayed greater diversity in their diet. Additionally, we propose a revision of the definition of venom complexity in conoidean snails to more accurately capture the breadth of ecological influences. These findings suggest that prey diet is an important factor in terebrid venom evolution and diversification and further investigations of other understudied organisms, like terebrids, are needed to develop robust hypotheses in this area.

Multi-level convergence of complex traits and the evolution of bioluminescence

Evolutionary convergence provides natural opportunities to investigate how, when, and why novel traits evolve. Many convergent traits are complex, highlighting the importance of explicitly considering convergence at different levels of biological organization, or 'multi-level convergent evolution'. To investigate multi-level convergent evolution, we propose a holistic and hierarchical framework that emphasizes breaking down traits into several functional modules. We begin by identifying long-standing questions on the origins of complexity and the diverse evolutionary processes underlying phenotypic convergence to discuss how they can be addressed by examining convergent systems. We argue that bioluminescence, a complex trait that evolved dozens of times through either novel mechanisms or conserved toolkits, is particularly well suited for these studies. We present an updated estimate of at least 94 independent origins of bioluminescence across the tree of life, which we calculated by reviewing and summarizing all estimates of independent origins. Then, we use our framework to review the biology, chemistry, and evolution of bioluminescence, and for each biological level identify questions that arise from our systematic review. We focus on luminous organisms that use the shared luciferin substrates coelenterazine or vargulin to produce light because these organisms convergently evolved bioluminescent proteins that use the same luciferins to produce bioluminescence. Evolutionary convergence does not necessarily extend across biological levels, as exemplified by cases of conservation and disparity in biological functions, organs, cells, and molecules associated with bioluminescence systems. Investigating differences across bioluminescent organisms will address fundamental questions on predictability and contingency in convergent evolution. Lastly, we highlight unexplored areas of bioluminescence research and advances in sequencing and chemical techniques useful for developing bioluminescence as a model system for studying multi-level convergent evolution.

Transcriptomic Analysis of Four Cerianthid (Cnidaria, Ceriantharia) Venoms

Tube anemones, or cerianthids, are a phylogenetically informative group of cnidarians with complex life histories, including a pelagic larval stage and tube-dwelling adult stage, both known to utilize venom in stinging-cell rich tentacles. Cnidarians are an entirely venomous group that utilize their proteinaceous-dominated toxins to capture prey and defend against predators, in addition to several other ecological functions, including intraspecific interactions. At present there are no studies describing the venom for any species within cerianthids. Given their unique development, ecology, and distinct phylogenetic-placement within Cnidaria, our objective is to evaluate the venom-like gene diversity of four species of cerianthids from newly collected transcriptomic data. We identified 525 venom-like genes between all four species. The venom-gene profile for each species was dominated by enzymatic protein and peptide families, which is consistent with previous findings in other cnidarian venoms. However, we found few toxins that are typical of sea anemones and corals, and furthermore, three of the four species express toxin-like genes closely related to potent pore-forming toxins in box jellyfish. Our study is the first to provide a survey of the putative venom composition of cerianthids and contributes to our general understanding of the diversity of cnidarian toxins.