Diet Breadth Mediates the Prey Specificity of Venom Potency in Snakes

Prey Shifts Drive Venom Evolution in Cone Snails

Venom systems are complex traits that have independently emerged multiple times in diverse plant and animal phyla. Within each venomous lineage there typically exists interspecific variation in venom composition where several factors have been proposed as drivers of variation, including phylogeny and diet. Understanding these factors is of broad biological interest and has implications for the development of antivenom therapies and venom-based drug discovery. Because of their high species richness and the presence of several major evolutionary prey shifts, venomous marine cone snails (genus Conus) provide an ideal system to investigate drivers of interspecific venom variation. Here, by analyzing the venom gland expression profiles of ∼3,000 toxin genes from 42 species of cone snail, we elucidate the role of prey-specific selection pressures in shaping venom variation. By analyzing overall venom composition and individual toxin structures, we demonstrate that the shifts from vermivory to piscivory in Conus are complemented by distinct changes in venom composition independent of phylogeny. In vivo injections of venom from piscivorous cone snails in fish further showed a higher potency compared with venom of nonpiscivores demonstrating a selective advantage. Together, our findings provide compelling evidence for the role of prey shifts in directing the venom composition of cone snails and expand our understanding of the mechanisms of venom variation and diversification.

Prey specificity of predatory venoms

Venom represents a key adaptation of many venomous predators, allowing them to immobilise prey quickly through chemical rather than physical warfare. Evolutionary arms races between prey and a predator are believed to be the main factor influencing the potency and composition of predatory venoms. Predators with narrowly restricted diets are expected to evolve specifically potent venom towards their focal prey, with lower efficacy on alternative prey. Here, we evaluate hypotheses on the evolution of prey-specific venom, focusing on the effect of restricted diet, prey defences, and prey resistance. Prey specificity as a potential evolutionary dead end is also discussed. We then provide an overview of the current knowledge on venom prey specificity, with emphasis on snakes, cone snails, and spiders. As the current evidence for venom prey specificity is still quite limited, we also overview the best approaches and methods for its investigation and provide a brief summary of potential model groups. Finally, possible applications of prey-specific toxins are discussed.

Venomous Noodles: The Evolution of Toxins in Nemertea through Positive Selection and Gene Duplication

Some, probably most and perhaps all, members of the phylum Nemertea are poisonous, documented so far from marine and benthic specimens. Although the toxicity of these animals has been long known, systematic studies on the characterization of toxins, mechanisms of toxicity, and toxin evolution for this group are scarce. Here, we present the first investigation of the molecular evolution of toxins in Nemertea. Using a proteo-transcriptomic approach, we described toxins in the body and poisonous mucus of the pilidiophoran Lineus sanguineus and the hoplonemertean Nemertopsis pamelaroeae. Using these new and publicly available transcriptomes, we investigated the molecular evolution of six selected toxin gene families. In addition, we also characterized in silico the toxin genes found in the interstitial hoplonemertean, Ototyphlonemertes erneba, a meiofaunal taxa. We successfully identified over 200 toxin transcripts in each of these species. Evidence of positive selection and gene duplication was observed in all investigated toxin genes. We hypothesized that the increased rates of gene duplications observed for Pilidiophora could be involved with the expansion of toxin genes. Studies concerning the natural history of Nemertea are still needed to understand the evolution of their toxins. Nevertheless, our results show evolutionary mechanisms similar to other venomous groups.

Inhibition of bacterial biofilms by the snake venom proteome

Snake venoms possess a range of pharmacological and toxicological activities. Here we evaluated the antibacterial and anti-biofilm activity against methicillin-susceptible and methicillin-resistant Staphylococcus aureus (MSSA and MRSA) of venoms from the Samar spitting cobra Naja samarensis and the Puff adder Bitis arietans. Both venoms prevented biofilm production by pathogenic S. aureus in a growth-independent manner, with the B. arietans venom being most potent. Fractionation showed the active molecule to be heat-labile and >10 kDa in size. Proteomic profiles of N. samarensis venom revealed neurotoxins and cytotoxins, as well as an abundance of serine proteases and three-finger toxins, while serine proteases, metalloproteinases and C-lectin types were abundant in B. arietans venom. These enzymes may have evolved to prevent bacteria colonising the snake venom gland. From a biomedical biotechnology perspective, they have valuable potential for anti-virulence therapy to fight antibiotic resistant microbes.

Intersexual Differences in the Gene Expression of Phoneutria depilata (Araneae, Ctenidae) Toxins Revealed by Venom Gland Transcriptome Analyses

The wandering spider, Phoneutria depilata, is one of Colombia's most active nocturnal arthropod predators of vertebrates and invertebrates. Its venom has been a relevant subject of study in the last two decades. However, the scarcity of transcriptomic data for the species limits our knowledge of the distinct components present in its venom for linking the mainly neurotoxic effects of the spider venom to a particular molecular target. The transcriptome of the P. depilata venom gland was analyzed to understand the effect of different diets or sex and the impact of these variables on the composition of the venom. We sequenced venom glands obtained from ten males and ten females from three diet treatments: (i) invertebrate: Tenebrio molitor, (ii) vertebrate: Hemidactylus frenatus, and (iii) mixed (T. molitor + H. frenatus). Of 17,354 assembled transcripts from all samples, 65 transcripts relating to venom production differed between males and females. Among them, 36 were classified as neurotoxins, 14 as serine endopeptidases, 11 as other proteins related to venom production, three as metalloprotease toxins, and one as a venom potentiator. There were no differences in transcripts across the analyzed diets, but when considering the effect of diets on differences between the sexes, 59 transcripts were differentially expressed. Our findings provide essential information on toxins differentially expressed that can be related to sex and the plasticity of the diet of P. depilata and thus can be used as a reference for venomics of other wandering spider species.

Venom phenotype conservation suggests integrated specialization in a lizard-eating snake

Biological specialization reduces the size of niche space while increasing efficiency in the use of available resources. Specialization often leads to phenotypic changes via natural selection aligning with niche space constraints. Commonly observed changes are in size, shape, behavior, and traits associated with feeding. One often selected trait for dietary specialization is venom, which, in snakes, often shows variation dependent on diet across and within species. The Neotropical Blunt-headed Treesnake (Imantodes cenchoa) is a highly specialized, rear-fanged, arboreal, lizard hunter that displays a long thin body, enlarged eyes, and a large Duvernoy's gland. However, toxin characterization of I. cenchoa has never been completed. Here, we use RNA-seq and mass spectrometry to assemble, annotate, and analyze the venom gland transcriptomes of four I. cenchoa from across their range. We find a lack of significant venom variation at the sequence and expression levels, suggesting venom conservation across the species. We propose this conservation provides evidence of a specialized venom repertoire, adapted to maximize efficiency of capturing and processing lizards. Importantly, this study provides the most complete venom gland transcriptomes of I. cenchoa and evidence of venom specialization in a rear-fanged snake, giving insight into selective pressures of venom across all snake species.

Pangenomics of the death cap mushroom Amanita phalloides, and of Agaricales, reveals dynamic evolution of toxin genes in an invasive range

The poisonous European mushroom Amanita phalloides (the "death cap") is invading California. Whether the death caps' toxic secondary metabolites are evolving as it invades is unknown. We developed a bioinformatic pipeline to identify the MSDIN genes underpinning toxicity and probed 88 death cap genomes from an invasive Californian population and from the European range, discovering a previously unsuspected diversity of MSDINs made up of both core and accessory elements. Each death cap individual possesses a unique suite of MSDINs, and toxin genes are significantly differentiated between Californian and European samples. MSDIN genes are maintained by strong natural selection, and chemical profiling confirms MSDIN genes are expressed and result in distinct phenotypes; our chemical profiling also identified a new MSDIN peptide. Toxin genes are physically clustered within genomes. We contextualize our discoveries by probing for MSDINs in genomes from across the order Agaricales, revealing MSDIN diversity originated in independent gene family expansions among genera. We also report the discovery of an MSDIN in an Amanita outside the "lethal Amanitas" clade. Finally, the identification of an MSDIN gene and its associated processing gene (POPB) in Clavaria fumosa suggest the origin of MSDINs is older than previously suspected. The dynamic evolution of MSDINs underscores their potential to mediate ecological interactions, implicating MSDINs in the ongoing invasion. Our data change the understanding of the evolutionary history of poisonous mushrooms, emphasizing striking parallels to convergently evolved animal toxins. Our pipeline provides a roadmap for exploring secondary metabolites in other basidiomycetes and will enable drug prospecting.

Diet of a threatened rattlesnake (eastern massasauga) revealed by DNA metabarcoding

Characterizing the diet of imperiled species using minimally invasive methods is crucial to understanding their ecology and conservation requirements. Here, we apply a DNA metabarcoding approach to study the diet of the eastern massasauga rattlesnake (Sistrurus catenatus), a Federally Threatened snake found throughout the Great Lakes region. Eighty-three fecal samples collected across 10 different massasauga populations located in Michigan, USA, were sequenced, with 70 samples containing prey DNA. We used universal metazoan primers and developed a host-specific oligonucleotide blocker to characterize their diet. We identified at least 12 different prey species, with eastern massasaugas exhibiting opportunistic feeding and a strong preference towards small mammals. Meadow voles (Microtus pennsylvanicus) were the most common prey item (70% of diet) followed by the northern short-tailed shrew (Blarina brevicauda) and masked shrew (Sorex cinereus; 15.7% of diet each), along with occasional bird and snake prey. Adult individuals exhibited a more generalized diet, consuming a larger number of prey taxa on average. Younger snakes consumed a smaller variety of prey items and tended to consume smaller-sized mammals such as masked shrews (Sorex cinereus) and northern short-tailed shrews (Blarina brevicauda). We conclude that small mammals are a crucial part of eastern massasauga rattlesnake diet and recommend this be taken into consideration when conservation strategies are developed. The methods developed in this study can be applied to other reptile species, providing an accurate, minimally invasive, and thorough diet assessment for at-risk reptile species.

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.

Composition and toxicity of venom produced by araneophagous white-tailed spiders (Lamponidae: Lampona sp.)

Prey-specialised spiders are adapted to capture specific prey items, including dangerous prey. The venoms of specialists are often prey-specific and less complex than those of generalists, but their venom composition has not been studied in detail. Here, we investigated the venom of the prey-specialised white-tailed spiders (Lamponidae: Lampona), which utilise specialised morphological and behavioural adaptations to capture spider prey. We analysed the venom composition using proteo-transcriptomics and taxon-specific toxicity using venom bioassays. Our analysis identified 208 putative toxin sequences, comprising 103 peptides < 10 kDa and 105 proteins > 10 kDa. Most peptides belonged to one of two families characterised by scaffolds containing eight or ten cysteine residues. Toxin-like proteins showed similarity to galectins, leucine-rich repeat proteins, trypsins and neprilysins. The venom of Lampona was shown to be more potent against the preferred spider prey than against alternative cricket prey. In contrast, the venom of a related generalist was similarly potent against both prey types. These data provide insights into the molecular adaptations of venoms produced by prey-specialised spiders.

Divergent Specialization of Simple Venom Gene Profiles among Rear-Fanged Snake Genera (Helicops and Leptodeira, Dipsadinae, Colubridae)

Many venomous animals express toxins that show extraordinary levels of variation both within and among species. In snakes, most studies of venom variation focus on front-fanged species in the families Viperidae and Elapidae, even though rear-fanged snakes in other families vary along the same ecological axes important to venom evolution. Here we characterized venom gland transcriptomes from 19 snakes across two dipsadine rear-fanged genera (Leptodeira and Helicops, Colubridae) and two front-fanged genera (Bothrops, Viperidae; Micrurus, Elapidae). We compared patterns of composition, variation, and diversity in venom transcripts within and among all four genera. Venom gland transcriptomes of rear-fanged Helicops and Leptodeira and front-fanged Micrurus are each dominated by expression of single toxin families (C-type lectins, snake venom metalloproteinase, and phospholipase A2, respectively), unlike highly diverse front-fanged Bothrops venoms. In addition, expression patterns of congeners are much more similar to each other than they are to species from other genera. These results illustrate the repeatability of simple venom profiles in rear-fanged snakes and the potential for relatively constrained venom composition within genera.

Untangling interactions between Bitis vipers and their prey using coagulotoxicity against diverse vertebrate plasmas

Venom is a key evolutionary innovation which plays a primary role in prey subjugation by venomous snakes. However, while there is a growing body of literature indicating the composition and activity of snake venoms is under strong natural selection driven by differences in prey physiology, the majority of studies have historically focussed on the activity of snake venoms with regards only towards human or mammalian physiologies. This study aimed to use clotting assays measuring both time and strength of clotting to characterise the coagulotoxic activity of venoms from a taxonomically, morphologically, and ecologically diverse range of Bitis spp. of viperid snakes upon the plasma of model species: amphibian (Cane Toad, Rhinella marina); lizard (Blue-tongue Skink, Tiliqua scincoides); avian (Domestic Chicken, Gallus gallus); and rodent (Brown Rat, Rattus norvegicus). Significant variation in coagulotoxic activity across the different plasmas was observed between species and compared to the known affects upon human plasma. Bitis caudalis was notable in being active on all four plasmas, but in extremely divergent manners: accelerating clotting times and producing strong, stable clots upon amphibian plasma (consistent with true procoagulation); accelerating clotting time but producing weak, unstable clots upon lizard plasma (consistent with pseudo-procoagulation); delaying avian clotting time beyond machine maximum reading time (strong anticoagulation consistent with either inhibition of clotting enzymes or total destruction of fibrinogen, or both); and delaying clotting of rodent plasma (consistent with inhibition of clotting enzymes) and with only weak clots formed (consistent with destruction of fibrinogen). In contrast, the sister species B. peringueyi and B. schneideri displayed activity only upon the lizard plasma, slightly accelerating the clotting times to produce weak, unstable clots (consistent with pseudo-procoagulation). The other dwarf species, B. cornuta, displayed strong anticoagulation upon avian and rodent plasmas, delaying clotting beyond the machine maximum reading time (strong anticoagulation consistent with either inhibition of clotting enzymes or total destruction of fibrinogen, or both). In contrast, the giant species studied (B. gabonica) showed only a very weak pseudo-procoagulant activity upon lizard plasma. The wide range of variation seen within this study highlights the importance of studying venom activity on relevant models when making conclusions about the ecological role of venoms and the extreme limitation in extrapolating animal results to predict potential human clinical effects.

Differential coagulotoxic and neurotoxic venom activity from species of the arboreal viperid snake genus Bothriechis (palm-pitvipers)

The viperid snake genus Bothriechis consists of eleven species distributed among Central and South America, living across low and high-altitude habitats. Despite Bothriechis envenomations being prominent across the Central and South American region, the functional effects of Bothriechis venoms are poorly understood. Thus, the aim of this study was to investigate the coagulotoxic and neurotoxic activities of Bothriechis venoms to fill this knowledge gap. Coagulotoxic investigations revealed Bothriechis nigroviridis and B. schlegelii to have pseudo-procoagulant venom activity, forming weak clots that rapidly break down, thereby depleting fibrinogen levels and thus contributing to a net anticoagulant state. While one sample of B. lateralis also showed weaker pseudo-procoagulant activity, directly clotting fibrinogen, two samples of B. lateralis venom were anticoagulant through the inhibition of thrombin and factor Xa activity. Differential efficacy of PoliVal-ICP antivenom was also observed, with the pseudo-procoagulant effect of B. nigroviridis venom poorly neutralised, despite this same activity in the venom of B. schlegelii being effectively neutralised. Significant specificity of these fibrinogen cleaving toxins was also observed, with no activity upon model amphibian, avian, lizard or rodent plasma observed. However, upon avian plasma the venom of B. nigroviridis exerted a complete anticoagulant effect, in contrast to the pseudo-procoagulant effect seen on human plasma. Neurotoxic investigations revealed B. schlegelii to be unique among the genus in having potent binding to the orthosteric site of the alpha-1 postsynaptic nicotinic acetylcholine receptor (with B. lateralis having a weaker but still discernible effect). This represents the first identification of postsynaptic nAChR neurotoxic activity for Bothriechis. In conclusion this study identifies notable differential activity within the coagulotoxic and postsynaptic neurotoxic activity of Bothriechis venoms, supporting previous research, and highlights the need for further studies with respect to antivenom efficacy as well as coagulotoxin specificity for Bothriechis venoms.

Interpopulational variation and ontogenetic shift in the venom composition of Lataste's viper (Vipera latastei, Boscá 1878) from northern Portugal

Lataste's viper (Vipera latastei) is a venomous European viper endemic to the Iberian Peninsula, recognised as medically important by the World Health Organization. To date, no comprehensive characterisation of this species' venom has been reported. Here, we analysed the venoms of juvenile and adult specimens of V. latastei from two environmentally different populations from northern Portugal. Using bottom-up venomics, we produced six venom proteomes (three per population) from vipers belonging to both age classes (i.e., two juveniles and four adults), and RP-HPLC profiles of 54 venoms collected from wild specimens. Venoms from juveniles and adults differed in their chromatographic profiles and relative abundances of their toxins, suggesting the occurrence of ontogenetic changes in venom composition. Specifically, snake venom metalloproteinase (SVMP) was the most abundant toxin family in juvenile venoms, while snake venom serine proteinases (SVSPs), phospholipases A₂ (PLA₂s), and C-type lectin-like (CTLs) proteins were the main toxins comprising adult venoms. The RP-HPLC venom profiles were found to vary significantly between the two sampled localities, indicating geographic variability. Furthermore, the presence/absence of certain peaks in the venom chromatographic profiles appeared to be significantly correlated also to factors like body size and sex of the vipers. Our findings show that V. latastei venom is a variable phenotype. The intraspecific differences we detected in its composition likely mirror changes in the feeding ecology of this species, taking place during different life stages and under different environmental pressures. SIGNIFICANCE: Lataste's viper (Vipera latastei) is a medically important viper endemic to the Iberian Peninsula, inhabiting different habitats and undergoing a marked ontogenetic dietary shift. In the current study, we report the first proteomic analysis of V. latastei venom from two environmentally different localities in northern Portugal. Our bottom-up venomic analyses show that snake venom serine proteinases (SVSPs), phospholipases A2 (PLA2s), and C-type lectin-like (CTLs) proteins are the major components of adult V. latastei venom. The comparative analysis of young and adult venoms suggests the occurrence of ontogenetic shift in toxin abundances, with snake venom metalloproteinases (SVMPs) being the predominant toxins in juvenile venoms. Moreover, geographic venom variation between the two studied populations is also detected, with our statistical analyses suggesting that factors like body size and sex of the vipers are possibly at play in its determination. Our work represents the first assessment of the composition of V. latastei venom, and the first step towards a better understanding of the drivers behind its variability.

Scorpion Species with Smaller Body Sizes and Narrower Chelae Have the Highest Venom Potency

Scorpionism is a global health concern, with an estimation of over one million annual envenomation cases. Despite this, little is known regarding the drivers of scorpion venom potency. One widely held view is that smaller scorpions with less-developed chelae possess the most potent venoms. While this perception is often used as a guide for medical intervention, it has yet to be tested in a formal comparative framework. Here, we use a phylogenetic comparative analysis of 36 scorpion species to test whether scorpion venom potency, as measured using LD50, is related to scorpion body size and morphology. We found a positive relationship between LD50 and scorpion total length, supporting the perception that smaller scorpions possess more potent venoms. We also found that, independent of body size, scorpion species with long narrow chelae have higher venom potencies compared to species with more robust chelae. These results not only support the general perception of scorpion morphology and potency, but also the presence of an ecology trade-off with scorpions either selected for well-developed chelae or more potent venoms. Testing the patterns of venom variations in scorpions aids both our ecological understanding and our ability to address the global health burden of scorpionism.

Linking toxicity and predation in a venomous arthropod: the case of Tityus fuhrmanni (Scorpiones: Buthidae), a generalist predator scorpion

Background:

Scorpions are arachnids that have a generalist diet, which use venom to subdue their prey. The study of their trophic ecology and capture behavior is still limited compared to other organisms, and aspects such as trophic specialization in this group have been little explored.

Methods:

In order to determine the relationship between feeding behavior and venom toxicity in the scorpion species Tityus fuhrmanni, 33 specimens were offered prey with different morphologies and defense mechanisms: spiders, cockroaches and crickets. In each of the experiments we recorded the following aspects: acceptance rate, immobilization time and the number of capture attempts. The median lethal dose of T. fuhrmanni venom against the three different types of prey was also evaluated.

Results:

We found that this species does not have a marked difference in acceptance for any of the evaluated prey, but the number of capture attempts of spiders is higher when compared to the other types of prey. The immobilization time is shorter in spiders compared to other prey and the LD₅₀ was higher for cockroaches.

Conclusions:

These results indicate that T. fuhrmanni is a scorpion with a generalist diet, has a venom with a different potency among prey and is capable of discriminating between prey types and employing distinct strategies to subdue them.

Promoting co-existence between humans and venomous snakes through increasing the herpetological knowledge base

Snakebite incidence at least partly depends on the biology of the snakes involved. However, studies of snake biology have been largely neglected in favour of anthropic factors, with the exception of taxonomy, which has been recognised for some decades to affect the design of antivenoms. Despite this, within-species venom variation and the unpredictability of the correlation with antivenom cross-reactivity has continued to be problematic. Meanwhile, other aspects of snake biology, including behaviour, spatial ecology and activity patterns, distribution, and population demography, which can contribute to snakebite mitigation and prevention, remain underfunded and understudied. Here, we review the literature relevant to these aspects of snakebite and illustrate how demographic, spatial, and behavioural studies can improve our understanding of why snakebites occur and provide evidence for prevention strategies. We identify the large gaps that remain to be filled and urge that, in the future, data and relevant metadata be shared openly via public data repositories so that studies can be properly replicated and data used in future meta-analyses.

Armed stem to stinger: a review of the ecological roles of scorpion weapons

Scorpions possess two systems of weapons: the pincers (chelae) and the stinger (telson). These are placed on anatomically and developmentally well separated parts of the body, that is, the oral appendages and at the end of the body axis. The otherwise conserved body plan of scorpions varies most in the shape and relative dimensions of these two weapon systems, both across species and in some cases between the sexes. We review the literature on the ecological function of these two weapon systems in each of three contexts of usage: (i) predation, (ii) defense and (iii) sexual contests. In the latter context, we will also discuss their usage in mating. We first provide a comparative background for each of these contexts of usage by giving examples of other weapon systems from across the animal kingdom. Then, we discuss the pertinent aspects of the anatomy of the weapon systems, particularly those aspects relevant to their functioning in their ecological roles. The literature on the functioning and ecological role of both the chelae and the telson is discussed in detail, again organized by context of usage. Particular emphasis is given on the differences in morphology or usage between species or higher taxonomic groups, or between genders, as such cases are most insightful to understand the roles of each of the two distinct weapon systems of the scorpions and their evolutionary interactions. We aimed to synthesize the literature while minimizing conjecture, but also to point out gaps in the literature and potential future research opportunities.

Evidence for Resistance to Coagulotoxic Effects of Australian Elapid Snake Venoms by Sympatric Prey (Blue Tongue Skinks) but Not by Predators (Monitor Lizards)

Some Australian elapids possess potently procoagulant coagulotoxic venoms which activate the zymogen prothrombin into the functional enzyme thrombin. Although the activity of Australian elapid prothrombin-activators has been heavily investigated with respect to the mammalian, and in particular, human clotting cascades, very few studies have investigated the activity of their venom upon reptile plasmas. This is despite lizards representing both the primary diet of most Australian elapids and also representing natural predators. This study investigated the procoagulant actions of a diverse range of Australian elapid species upon plasma from known prey species within the genera Tiliqua (blue tongue skinks) as well as known predator species within the genera Varanus (monitor lizards). In addition to identifying significant variation in the natural responses of the coagulation cascade between species from the genera Tiliqua and Varanus relative to each other, as well as other vertebrate lineages, notable differences in venom activity were also observed. Within the genus Tiliqua, both T. rugosa and T. scincoides plasma displayed significant resistance to the procoagulant activity of Pseudechis porphyriacus venom, despite being susceptible to all other procoagulant elapid venoms. These results indicate that T. rugosa and T. scincoides have evolved resistance within their plasma to the coagulotoxic venom activity of the sympatric species P. porphyriacus. Other venoms were able to activate Tiliqua prothrombin, which suggests that the lessened activity of P. porphyriacus venom is not due to modifications of the prothrombin and may instead be due to a serum factor that specifically binds to P. porphyriacus toxins, as has been previously seen for squirrels resistant to rattlesnake venom. In contrast, none of the predatory lizards studied (Varanus giganteus, V. mertensi and V. varius) demonstrated resistance to the venom. This suggests that the mechanical protection afforded by thick osteodermic scales, and prey handling behaviour, removes a selection pressure for the evolution of resistance in these large predatory lizards. These results therefore reveal differential interactions between venoms of snakes with sympatric lizards that are on opposite sides of the predator-prey arms race.

A symphony of destruction: Dynamic differential fibrinogenolytic toxicity by rattlesnake (Crotalus and Sistrurus) venoms

What factors influence the evolution of a heavily selected functional trait in a diverse clade? This study adopts rattlesnakes as a model group to investigate the evolutionary history of venom coagulotoxicity in the wider context of phylogenetics, natural history, and biology. Venom-induced clotting of human plasma and fibrinogen was determined and mapped onto the rattlesnake phylogenetic tree to reconstruct the evolution of coagulotoxicity across the group. Our results indicate that venom phenotype is often independent of phylogenetic relationships in rattlesnakes, suggesting the importance of diet and/or other environmental variables in driving venom evolution. Moreover, the striking inter- and intraspecific variability in venom activity on human blood highlights the considerable variability faced by physicians treating envenomation. This study is the most comprehensive effort to date to describe and characterize the evolutionary and biological aspects of coagulotoxins in rattlesnake venom. Further research at finer taxonomic levels is recommended to elucidate patterns of variation within species and lineages.

Tentacle Morphological Variation Coincides with Differential Expression of Toxins in Sea Anemones

Phylum Cnidaria is an ancient venomous group defined by the presence of cnidae, specialised organelles that serve as venom delivery systems. The distribution of cnidae across the body plan is linked to regionalisation of venom production, with tissue-specific venom composition observed in multiple actiniarian species. In this study, we assess whether morphological variants of tentacles are associated with distinct toxin expression profiles and investigate the functional significance of specialised tentacular structures. Using five sea anemone species, we analysed differential expression of toxin-like transcripts and found that expression levels differ significantly across tentacular structures when substantial morphological variation is present. Therefore, the differential expression of toxin genes is associated with morphological variation of tentacular structures in a tissue-specific manner. Furthermore, the unique toxin profile of spherical tentacular structures in families Aliciidae and Thalassianthidae indicate that vesicles and nematospheres may function to protect branched structures that host a large number of photosynthetic symbionts. Thus, hosting zooxanthellae may account for the tentacle-specific toxin expression profiles observed in the current study. Overall, specialised tentacular structures serve unique ecological roles and, in order to fulfil their functions, they possess distinct venom cocktails.

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.

Duvernoy's Gland Transcriptomics of the Plains Black-Headed Snake, Tantilla nigriceps (Squamata, Colubridae): Unearthing the Venom of Small Rear-Fanged Snakes

The venoms of small rear-fanged snakes (RFS) remain largely unexplored, despite increased recognition of their importance in understanding venom evolution more broadly. Sequencing the transcriptome of venom-producing glands has greatly increased the ability of researchers to examine and characterize the toxin repertoire of small taxa with low venom yields. Here, we use RNA-seq to characterize the Duvernoy's gland transcriptome of the Plains Black-headed Snake, Tantilla nigriceps, a small, semi-fossorial colubrid that feeds on a variety of potentially dangerous arthropods including centipedes and spiders. We generated transcriptomes of six individuals from three localities in order to both characterize the toxin expression of this species for the first time, and to look for initial evidence of venom variation in the species. Three toxin families-three-finger neurotoxins (3FTxs), cysteine-rich secretory proteins (CRISPs), and snake venom metalloproteinases (SVMPIIIs)-dominated the transcriptome of T. nigriceps; 3FTx themselves were the dominant toxin family in most individuals, accounting for as much as 86.4% of an individual's toxin expression. Variation in toxin expression between individuals was also noted, with two specimens exhibiting higher relative expression of c-type lectins than any other sample (8.7-11.9% compared to <1%), and another expressed CRISPs higher than any other toxin. This study provides the first Duvernoy's gland transcriptomes of any species of Tantilla, and one of the few transcriptomic studies of RFS not predicated on a single individual. This initial characterization demonstrates the need for further study of toxin expression variation in this species, as well as the need for further exploration of small RFS venoms.

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 [...].

Potency Testing of Venoms and Antivenoms in Embryonated Eggs: An Ethical Alternative to Animal Testing

Venoms are complex mixtures of biologically active molecules that impact multiple physiological systems. Manufacture of antivenoms (AVs) therefore requires potency testing using in vivo models to ensure AV efficacy. As part of ongoing research to replace small animals as the standard model for AV potency testing, we developed an alternate in vivo method using the embryonated egg model (EEM). In this model, the survival of chicken embryos envenomated in ovo is determined prior to 50% gestation, when they are recognized as animals by animal welfare legislation. Embryos were found to be susceptible to a range of snake, spider, and marine venoms. This included funnel-web spider venom for which the only other vertebrate, non-primate animal model is newborn mice. Neutralization of venom with standard AV allowed correlation of AV potency results from the EEM to results from animal assays. Our findings indicate that the EEM provides an alternative, insensate in vivo model for the assessment of AV potency. The EEM may enable reduction or replacement of the use of small animals, as longer-term research that enables the elimination of animal use in potency testing continues.

Proteomic and Transcriptomic Techniques to Decipher the Molecular Evolution of Venoms

Nature's library of venoms is a vast and untapped resource that has the potential of becoming the source of a wide variety of new drugs and therapeutics. The discovery of these valuable molecules, hidden in diverse collections of different venoms, requires highly specific genetic and proteomic sequencing techniques. These have been used to sequence a variety of venom glands from species ranging from snakes to scorpions, and some marine species. In addition to identifying toxin sequences, these techniques have paved the way for identifying various novel evolutionary links between species that were previously thought to be unrelated. Furthermore, proteomics-based techniques have allowed researchers to discover how specific toxins have evolved within related species, and in the context of environmental pressures. These techniques allow groups to discover novel proteins, identify mutations of interest, and discover new ways to modify toxins for biomimetic purposes and for the development of new therapeutics.

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.

Utilising venom activity to infer dietary composition of the Kenyan horned viper (Bitis worthingtoni)

Bitis are well known for being some of the most commonly encountered and medically important snake species in all of Africa. While the majority of species possess potently anticoagulant venom, only B. worthingtoni is known to possess procoagulant venom. Although known to be the basal species within the genus, B. worthingtoni is an almost completely unstudied species with even basic dietary information lacking. This study investigated various aspects of the unique procoagulant effects of B. worthingtoni venom. Coagulation assays determined the primary procoagulant effect to be driven by Factor X activating snake venom metalloprotease toxins. In addition to acting upon the mammalian blood clotting cascade, B. worthingtoni venom was also shown to clot amphibian plasma. As previous studies have shown differences in clotting factors between amphibian and mammalian plasmas, individual enzymes in snake venoms acting on plasma clotting factors can be taxon-selective. As venoms evolve under purifying selection pressures, this suggests that the procoagulant snake venom metalloprotease toxins present in B. worthingtoni have likely been retained from a recent common ancestor shared with the related amphibian-feeding Proatheris superciliaris, and that both amphibians and mammals represent a substantial proportion of B. worthingtoni current diet. Thus, taxon-specific actions of venoms may have utility in inferring dietary composition for rare or difficult to study species. An important caveat is that to validate this hypothesis field studies investigating the dietary ecology of B. worthingtoni must be conducted, as well as further investigations of its venom composition to reconstruct the molecular evolutionary history of the toxins present.

Assessing the Binding of Venoms from Aquatic Elapids to the Nicotinic Acetylcholine Receptor Orthosteric Site of Different Prey Models

The evolution of an aquatic lifestyle from land dwelling venomous elapids is a radical ecological modification, bringing about many evolutionary changes from morphology to diet. Diet is an important ecological facet which can play a key role in regulating functional traits such as venom composition and prey-specific targeting of venom. In addition to predating upon novel prey (e.g., fish, fish eggs and invertebrates), the venoms of aquatic elapids also face the challenge of increased prey-escape potential in the aquatic environment. Thus, despite the independent radiation into an aquatic niche on four separate occasions, the venoms of aquatic elapids are evolving under convergent selection pressures. Utilising a biolayer interferometry binding assay, this study set out to elucidate whether crude venoms from representative aquatic elapids were target-specific to the orthosteric site of postsynaptic nicotinic acetylcholine receptor mimotopes of fish compared to other terrestrial prey types. Representatives of the four aquatic lineages were: aquatic coral snakes representative was Micrurus surinamensis;, sea kraits representative was Laticauda colubrina; sea snakes representatives were two Aipysurus spp. and eight Hydrophis spp; and water cobras representative was Naja annulata. No prey-specific differences in crude venom binding were observed from any species tested, except for Aipysurus laevis, which showed slight evidence of prey-potency differences. For Hydrophis caerulescens, H. peronii, H. schistosus and M. surinamensis, there was a lack of binding to the orthosteric site of any target lineage. Subsequent testing on the in vitro chick-biventer cervicis muscle preparation suggested that, while the venoms of these species bound postsynaptically, they bound to allosteric sites rather than orthosteric. Allosteric binding is potentially a weaker but faster-acting form of neurotoxicity and we hypothesise that the switch to allosteric binding is likely due to selection pressures related to prey-escape potential. This research has potentially opened up the possibility of a new functional class of toxins which have never been assessed previously while shedding light on the selection pressures shaping venom evolution.

An Economic Dilemma Between Molecular Weapon Systems May Explain an Arachno-atypical Venom in Wasp Spiders (Argiope bruennichi)

Spiders use venom to subdue their prey, but little is known about the diversity of venoms in different spider families. Given the limited data available for orb-weaver spiders (Araneidae), we selected the wasp spider Argiope bruennichi for detailed analysis. Our strategy combined a transcriptomics pipeline based on multiple assemblies with a dual proteomics workflow involving parallel mass spectrometry techniques and electrophoretic profiling. We found that the remarkably simple venom of A. bruennichi has an atypical composition compared to other spider venoms, prominently featuring members of the cysteine-rich secretory protein, antigen 5 and pathogenesis-related protein 1 (CAP) superfamily and other, mostly high-molecular-weight proteins. We also detected a subset of potentially novel toxins similar to neuropeptides. We discuss the potential function of these proteins in the context of the unique hunting behavior of wasp spiders, which rely mostly on silk to trap their prey. We propose that the simplicity of the venom evolved to solve an economic dilemma between two competing yet metabolically expensive weapon systems. This study emphasizes the importance of cutting-edge methods to encompass the lineages of smaller venomous species that have yet to be characterized in detail, allowing us to understand the biology of their venom systems and to mine this prolific resource for translational research.

Venomics Approach Reveals a High Proportion of Lactrodectus-Like Toxins in the Venom of the Noble False Widow Spider Steatoda nobilis

The noble false widow spider Steatoda nobilis originates from the Macaronesian archipelago and has expanded its range globally. Outside of its natural range, it may have a negative impact on native wildlife, and in temperate regions it lives in synanthropic environments where it frequently encounters humans, subsequently leading to envenomations. S. nobilis is the only medically significant spider in Ireland and the UK, and envenomations have resulted in local and systemic neurotoxic symptoms similar to true black widows (genus Latrodectus). S. nobilis is a sister group to Latrodectus which possesses the highly potent neurotoxins called α-latrotoxins that can induce neuromuscular paralysis and is responsible for human fatalities. However, and despite this close relationship, the venom composition of S. nobilis has never been investigated. In this context, a combination of transcriptomic and proteomic cutting-edge approaches has been used to deeply characterise S. nobilis venom. Mining of transcriptome data for the peptides identified by proteomics revealed 240 annotated sequences, of which 118 are related to toxins, 37 as enzymes, 43 as proteins involved in various biological functions, and 42 proteins without any identified function to date. Among the toxins, the most represented in numbers are α-latrotoxins (61), δ-latroinsectotoxins (44) and latrodectins (6), all of which were first characterised from black widow venoms. Transcriptomics alone provided a similar representation to proteomics, thus demonstrating that our approach is highly sensitive and accurate. More precisely, a relative quantification approach revealed that latrodectins are the most concentrated toxin (28%), followed by α-latrotoxins (11%), δ-latroinsectotoxins (11%) and α-latrocrustotoxins (11%). Approximately two-thirds of the venom is composed of Latrodectus-like toxins. Such toxins are highly potent towards the nervous system of vertebrates and likely responsible for the array of symptoms occurring after envenomation by black widows and false widows. Thus, caution should be taken in dismissing S. nobilis as harmless. This work paves the way towards a better understanding of the competitiveness of S. nobilis and its potential medical importance.

Evolutionary Interpretations of Nicotinic Acetylcholine Receptor Targeting Venom Effects by a Clade of Asian Viperidae Snakes

Ecological variability among closely related species provides an opportunity for evolutionary comparative studies. Therefore, to investigate the origin and evolution of neurotoxicity in Asian viperid snakes, we tested the venoms of Azemiops feae, Calloselasma rhodostoma, Deinagkistrodon acutus, Tropidolaeums subannulatus, and T. wagleri for their relative specificity and potency upon the amphibian, lizard, bird, rodent, and human α-1 (neuromuscular) nicotinic acetylcholine receptors. We utilised a biolayer interferometry assay to test the binding affinity of these pit viper venoms to orthosteric mimotopes of nicotinic acetylcholine receptors binding region from a diversity of potential prey types. The Tropidolaemus venoms were much more potent than the other species tested, which is consistent with the greater prey escape potential in arboreal niches. Intriguingly, the venom of C. rhodostoma showed neurotoxic binding to the α-1 mimotopes, a feature not known previously for this species. The lack of prior knowledge of neurotoxicity in this species is consistent with our results due to the bias in rodent studies and human bite reports, whilst this venom had a greater binding affinity toward amphibian and diapsid α-1 targets. The other large terrestrial species, D. acutus, did not display any meaningful levels of neurotoxicity. These results demonstrate that whilst small peptide neurotoxins are a basal trait of these snakes, it has been independently amplified on two separate occasions, once in Azemiops and again in Tropidolaemus, and with Calloselasma representing a third possible amplification of this trait. These results also point to broader sources of novel neuroactive peptides with the potential for use as lead compounds in drug design and discovery.

Characterising Functional Venom Profiles of Anthozoans and Medusozoans within Their Ecological Context

This review examines the current state of knowledge regarding toxins from anthozoans (sea anemones, coral, zoanthids, corallimorphs, sea pens and tube anemones). We provide an overview of venom from phylum Cnidaria and review the diversity of venom composition between the two major clades (Medusozoa and Anthozoa). We highlight that the functional and ecological context of venom has implications for the temporal and spatial expression of protein and peptide toxins within class Anthozoa. Understanding the nuances in the regulation of venom arsenals has been made possible by recent advances in analytical technologies that allow characterisation of the spatial distributions of toxins. Furthermore, anthozoans are unique in that ecological roles can be assigned using tissue expression data, thereby circumventing some of the challenges related to pharmacological screening.

An Appetite for Destruction: Detecting Prey-Selective Binding of α-Neurotoxins in the Venom of Afro-Asian Elapids

Prey-selective venoms and toxins have been documented across only a few species of snakes. The lack of research in this area has been due to the absence of suitably flexible testing platforms. In order to test more species for prey specificity of their venom, we used an innovative taxonomically flexible, high-throughput biolayer interferometry approach to ascertain the relative binding of 29 α-neurotoxic venoms from African and Asian elapid representatives (26 Naja spp., Aspidelaps scutatus, Elapsoidea boulengeri, and four locales of Ophiophagus hannah) to the alpha-1 nicotinic acetylcholine receptor orthosteric (active) site for amphibian, lizard, snake, bird, and rodent targets. Our results detected prey-selective, intraspecific, and geographical differences of α-neurotoxic binding. The results also suggest that crude venom that shows prey selectivity is likely driven by the proportions of prey-specific α-neurotoxins with differential selectivity within the crude venom. Our results also suggest that since the α-neurotoxic prey targeting does not always account for the full dietary breadth of a species, other toxin classes with a different pathophysiological function likely play an equally important role in prey immobilisation of the crude venom depending on the prey type envenomated. The use of this innovative and taxonomically flexible diverse assay in functional venom testing can be key in attempting to understanding the evolution and ecology of α-neurotoxic snake venoms, as well as opening up biochemical and pharmacological avenues to explore other venom effects.