When one phenotype is not enough: divergent evolutionary trajectories govern venom variation in a widespread rattlesnake species

Proteogenomics-guided functional venomics resolves the toxin arsenal and activity of Deinagkistrodon acutus venom

Snakebite primarily impacts rural communities of Africa, Asia, and Latin America. The sharp-nosed viper (Deinagkistrodon acutus) is among the snakes of highest medical importance in Asia. Despite various studies on its venom using modern venomics techniques, a comprehensive understanding of composition and function of this species' venom remains lacking. We combined proteogenomics with extensive bioactivity profiling to present the first genome-level catalogue of D. acutus venom proteins and their exochemistry. Our analysis identified an unusually simple venom containing 45 components from 20 distinct protein families. Relative toxin abundances indicate that C-type lectin and C-type lectin-related protein (CTL), snake venom metalloproteinase (svMP), snake venom serine protease (svSP), and phospholipase A2 (PLA2) constitute 90 % of the venom. Bioassays targeting key aspects of viperid envenomation showed considerable concentration-dependent cytotoxicity, particularly in kidney and lung cells, and potent protease and PLA2 activity. Factor Xa and thrombin activities were minor, and no plasmin activity was observed. Effects on haemolysis, intracellular calcium (Ca2+) release, and nitric oxide (NO) synthesis were negligible. Our analysis provides the first holistic genome-based overview of the toxin arsenal of D. acutus, predicting the molecular and functional basis of its life-threatening effects, and opens novel avenues for treating envenomation by this highly dangerous snake.

Comparative venom analysis between melanistic and normally coloured phenotypes of the common adder (Vipera berus)

Snake venom is an ecologically relevant functional trait directly linked with a snake's fitness and survival, facilitating predation and defence. Snake venom variation occurs at all taxonomic levels, but the study at the intraspecific level is still in its early stages. The common adder (Vipera berus) exhibits considerable variation in colour phenotypes across its distribution range. Melanistic (fully black) individuals are the subject of myths and fairytales, and in German folklore such 'hell adders' are considered more toxic than their normally coloured conspecifics despite any formal investigation. Here, we provide the first comparative analysis of venoms from melanistic and normally coloured common adders. Specifically, we compared the venom profiles by sodium dodecylsulfate polyacrylamide gel electrophoresis and reversed-phase high-performance liquid chromatography and tested the venoms' protease, phospholipase A2 and cytotoxic activities. Phospholipase A2 activity was similar in both phenotypes, whereas general protease activity was higher in the melanistic venom, which was also more cytotoxic at two concentrations (6.25 and 12.5 µg ml-1). These minor differences between the venoms of melanistic and normally coloured adders are unlikely to be of clinical relevance in the context of human envenomation. In light of our results, the claim that melanistic adders produce more toxic venom than their normally coloured conspecifics appears rooted entirely in folklore.

Cold kiss still hot: limited temperature effects on envenomation performance in predatory strikes of a Japanese pit viper (Gloydius blomhoffii)

Understanding how environmental factors affect the performance of predators can provide profound insights into predator-prey interactions from evolutionary and ecological perspectives and the global distributional patterns of each taxon. Almost all venomous predators are ectotherms, with muscle contraction properties depending on temperature. For predators having venom transportation systems driven by muscle contraction, temperature may have quite large effects on envenomation performance for prey subjugation. Here, we used videography and enzyme-linked immunosorbent assay to examine thermal effects on envenomation kinematics and venom expenditure in predatory strikes of a venomous snake, the Mamushi Gloydius blomhoffii, to its main rodent prey at various body temperatures under both field and laboratory experimental conditions. Unexpectedly, we found that the thermal effects on envenomation performance are limited over nearly the entire ecologically relevant range of temperature (from 13.2°C to 26.2°C). Although temperature statistically significantly affected the mass of venom injected under field conditions, temperature explained only a minor proportion of the variation in venom expenditure. These findings suggest that the Mamushi is able to maintain prey subjugation performance across a wide range of temperatures, which is highly advantageous for ectothermic predators. Further studies should examine the underlying mechanisms of the limited thermal effects and their ubiquity across venomous predators.

Venom diversity in Naja mossambica: Insights from proteomic and immunochemical analyses reveal intraspecific differences

BACKGROUND

Intraspecific variations in snake venom composition have been extensively documented, contributing to the diverse clinical effects observed in envenomed patients. Understanding these variations is essential for developing effective snakebite management strategies and targeted antivenom therapies. We aimed to comprehensively investigate venoms from three distinct populations of N. mossambica from Eswatini, Limpopo, and KwaZulu-Natal regions in Africa in terms of their protein composition and reactivity with three commercial antivenoms (SAIMR polyvalent, EchiTAb+ICP, and Antivipmyn Africa).

METHODOLOGY/PRINCIPAL FINDINGS

Naja mossambica venoms from Eswatini region exhibited the highest content of neurotoxic proteins, constituting 20.70% of all venom proteins, compared to Limpopo (13.91%) and KwaZulu-Natal (12.80%), and was characterized by the highest diversity of neurotoxic proteins, including neurotoxic 3FTxs, Kunitz-type inhibitors, vespryns, and mamba intestinal toxin 1. KwaZulu-Natal population exhibited considerably lower cytotoxic 3FTx, higher PLA2 content, and significant diversity in low-abundant proteins. Conversely, Limpopo venoms demonstrated the least diversity as demonstrated by electrophoretic and mass spectrometry analyses. Immunochemical assessments unveiled differences in venom-antivenom reactivity, particularly concerning low-abundance proteins. EchiTAb+ICP antivenom demonstrated superior reactivity in serial dilution ELISA assays compared to SAIMR polyvalent.

CONCLUSIONS/SIGNIFICANCE

Our findings reveal a substantial presence of neurotoxic proteins in N. mossambica venoms, challenging previous understandings of their composition. Additionally, the detection of numerous peptides aligning to uncharacterized proteins or proteins with unknown functions underscores a critical issue with existing venom protein databases, emphasizing the substantial gaps in our knowledge of snake venom protein components. This underscores the need for enhanced research in this domain. Moreover, our in vitro immunological assays suggest EchiTAb+ICP's potential as an alternative to SAIMR antivenom, requiring confirmation through prospective in vivo neutralization studies.

A Review of Rattlesnake Venoms

Venom components are invaluable in biomedical research owing to their specificity and potency. Many of these components exist in two genera of rattlesnakes, Crotalus and Sistrurus, with high toxicity and proteolytic activity variation. This review focuses on venom components within rattlesnakes, and offers a comparison and itemized list of factors dictating venom composition, as well as presenting their known characteristics, activities, and significant applications in biosciences. There are 64 families and subfamilies of proteins present in Crotalus and Sistrurus venom. Snake venom serine proteases (SVSP), snake venom metalloproteases (SVMP), and phospholipases A2 (PLA2) are the standard components in Crotalus and Sistrurus venom. Through this review, we highlight gaps in the knowledge of rattlesnake venom; there needs to be more information on the venom composition of three Crotalus species and one Sistrurus subspecies. We discuss the activity and importance of both major and minor components in biomedical research and drug development.

Hunting behavior and feeding ecology of Mojave rattlesnakes (Crotalus scutulatus), prairie rattlesnakes (Crotalus viridis), and their hybrids in southwestern New Mexico

Predators must contend with numerous challenges to successfully find and subjugate prey. Complex traits related to hunting are partially controlled by a large number of co-evolved genes, which may be disrupted in hybrids. Accordingly, research on the feeding ecology of animals in hybrid zones has shown that hybrids sometimes exhibit transgressive or novel behaviors, yet for many taxa, empirical studies of predation and diet across hybrid zones are lacking. We undertook the first such field study for a hybrid zone between two snake species, the Mojave rattlesnake (Crotalus scutulatus) and the prairie rattlesnake (Crotalus viridis). Specifically, we leveraged established field methods to quantify the hunting behaviors of animals, their prey communities, and the diet of individuals across the hybrid zone in southwestern New Mexico, USA. We found that, even though hybrids had significantly lower body condition indices than snakes from either parental group, hybrids were generally similar to non-hybrids in hunting behavior, prey encounter rates, and predatory attack and success. We also found that, compared to C. scutulatus, C. viridis was significantly more active while hunting at night and abandoned ambush sites earlier in the morning, and hybrids tended to be more viridis-like in this respect. Prey availability was similar across the study sites, including within the hybrid zone, with kangaroo rats (Dipodomys spp.) as the most common small mammal, both in habitat surveys and the frequency of encounters with hunting rattlesnakes. Analysis of prey remains in stomachs and feces also showed broad similarity in diets, with all snakes preying primarily on small mammals and secondarily on lizards. Taken together, our results suggest that the significantly lower body condition of hybrids does not appear to be driven by differences in their hunting behavior or diet and may instead relate to metabolic efficiency or other physiological traits we have not yet identified.

On the importance of types and the perils of en passant taxonomy: a brief history of the typification of Coluber naja Linnaeus, 1758 (Serpentes: Elapidae) and its implications, with the designation of a lectotype

In response to the recent in passing (en passant) taxonomic decision to split Naja naja (Linnaeus) and recognise the Sri Lankan populations as a separate species, N. polyocellata Deraniyagala, we analyse the evidence underlying the proposal and its nomenclatural implications. The proposed split is weakly supported by the available evidence, so that retaining N. naja as a single species seems appropriate until further analysis. Moreover, the proposal raises several issues concerning types, type locality and nomenclature. Linnaeus description of Coluber naja was based on a single preserved specimen seen by him (now lost) and several illustrations in Sebas Thesaurus. The specimens that were the basis of these illustrations constitute part of the type series. Two of the latter specimens, ZMB 2795 and 2796, have been rediscovered in the collections of the Museum fr Naturkunde, Berlin. Here, we describe them, and determine that both are of Sri Lankan origin. To settle the question of the type and type locality of this iconic taxon, we designate ZMB 2796 as lectotype for the species, thereby implicitly restricting the type locality to Sri Lanka. The name polyocellata thus becomes a subjective junior synonym of Coluber naja, and the name Naja brasiliensis Laurenti, 1768 an objective junior synonym thereof. Any taxonomic recognition of additional diversity within N. naja would thus require the renaming of Indian, not Sri Lankan spectacled cobras, but should await a significant body of convincing evidence. We caution against taxonomic decisions taken in passing, based on limited evidence and without in-depth assessment of their nomenclatural implications.

Malaysian and Chinese King Cobra Venom Cytotoxicity in Melanoma and Neonatal Foreskin Fibroblasts Is Mediated by Age and Geography

Snake venoms constitute a complex, rapidly evolving trait, whose composition varies between and within populations depending on geographical location, age and preys (diets). These factors have determined the adaptive evolution for predatory success and link venom heterogeneity with prey specificity. Moreover, understanding the evolutionary drivers of animal venoms has streamlined the biodiscovery of venom-derived compounds as drug candidates in biomedicine and biotechnology. The king cobra (Ophiophagus hannah; Cantor, 1836) is distributed in diverse habitats, forming independent populations, which confer differing scale markings, including between hatchlings and adults. Furthermore, king cobra venoms possess unique cytotoxic properties that are used as a defensive trait, but their toxins may also have utility as promising anticancer-agent candidates. However, the impact of geographical distribution and age on these potential venom applications has been typically neglected. In this study, we hypothesised that ontogenetic venom variation accompanies the morphological distinction between hatchlings and adults. We used non-transformed neonatal foreskin (NFF) fibroblasts to examine and compare the variability of venom cytotoxicity between adult captive breeding pairs from Malaysian and Chinese lineages, along with that of their progeny upon hatching. In parallel, we assessed the anticancer potential of these venoms in human-melanoma-patient-derived cells (MM96L). We found that in a geographical distribution and gender-independent manner, venoms from hatchlings were significantly less cytotoxic than those from adults (NFF; ~Log EC50: 0.5-0.6 vs. 0.2-0.35 mg/mL). This is consistent with neonates occupying a semifossorial habitat, while adults inhabit more above-ground habitats and are therefore more conspicuous to potential predators. We also observed that Malaysian venoms exhibited a slightly higher cytotoxicity than those from the Chinese cobra cohorts (NFF; Log EC50: 0.1-0.3 vs. 0.3-0.4 mg/mL), which is consistent with Malaysian king cobras being more strongly aposematically marked. These variations are therefore suggestive of differential anti-predator strategies associated with the occupation of distinct niches. However, all cobra venoms were similarly cytotoxic in both melanoma cells and fibroblasts, limiting their potential medical applications in their native forms.

Venomics of the milos viper (Macrovipera schweizeri) unveils patterns of venom composition and exochemistry across blunt-nosed viper venoms

Introduction: Snakebite is a neglected tropical disease and a globally important driver of death and morbidity. Vipers of the genus Macrovipera (Viperidae: Viperinae) are among the snakes of higher medical importance in the Old World. Despite the medical relevance of Macrovipera venoms, the knowledge regarding them is heterogeneously distributed with virtually all works conducted so far focusing on subspecies of Macrovipera lebetinus, while other species within the genus are largely overlooked. Here we present the first proteomic evaluation of the venom from the Greek endemic Milos viper (Macrovipera schweizeri). In line with clinical symptoms typically elicited by Macrovipera envenomations, Milos viper venom primarily comprises coagulotoxic and cytotoxic protein families, such as metalloproteinases (svMP) and serine proteases (svSP). Methods: We conducted comparative bioactivity assays on venoms from M. schweizeri and the M. lebetinus subspecies M. lebetinus cernovi, M. lebetinus obtusa, and M. lebetinus turanica, and showed that they all exhibit similarities in levels of cytotoxicity proteolytic activity, and inhibition of prokaryotic growth. Lastly, we compared Macrovipera venom profiles by 1D-SDS-PAGE and RP-HPLC, as well as our proteomic data with previously published Macrovipera venom proteomes. Results and discussion: The analyzes performed to reveal that a general venom profile seems to be conserved across blunt-nosed vipers, and that, M. schweizeri envenomations, similarly to those caused by other blunt-nosed vipers, are able to cause significant tissue damage. The present work represents an important starting point for the development of comparative studies across the full taxonomic range of the genus Macrovipera and can potentially help optimize the treatment of envenomations caused by M. schweizeri.

Identifying the geographic distribution pattern of venomous snakes and regions of high snakebite risk in Iran

We describe species richness patterns of venomous snakes in Iran in order to produce snakebite risk prediction maps and identify gaps in regional health care centers capable of managing snakebites. We digitized distribution maps from the literature, Global Biodiversity Information Facility (GBIF), and the results of our own field studies of 24 terrestrial venomous snake species (including 4 endemic to Iran). Species richness patterns were associated with eight environmental factors. The variables have been extracted from the WorldClim dataset (bio12 = annual precipitation, bio15 = precipitation seasonality, bio17 = precipitation of the driest quarter, bio2 = mean diurnal range, bio3 = isothermality (bio2/bio7), bio4 = temperature seasonality, bio9 = mean temperature of the driest quarter and slope). Based on spatial analyses, species richness in Iran is highly affected by three environmental variables (bio12, 15, and 17) associated with precipitation. The relationship patterns among these predictors and species richness were strong and linear. The hotspot regions for venomous snakes species are concentrated on the western to southwestern and north to northeastern regions of Iran, which is partially consistent with the known Irano-Anatolian biodiversity hotspot. Because of the high number of endemic species and climatic conditions on the Iranian Plateau, the venoms of snakes distributed in those areas may contain novel properties and components.

Sequence Divergence in Venom Genes Within and Between Montane Pitviper (Viperidae: Crotalinae: Cerrophidion) Species is Driven by Mutation-Drift Equilibrium

Snake venom can vary both among and within species. While some groups of New World pitvipers-such as rattlesnakes-have been well studied, very little is known about the venom of montane pitvipers (Cerrophidion) found across the Mesoamerican highlands. Compared to most well-studied rattlesnakes, which are widely distributed, the isolated montane populations of Cerrophidion may facilitate unique evolutionary trajectories and venom differentiation. Here, we describe the venom gland transcriptomes for populations of C. petlalcalensis, C. tzotzilorum, and C. godmani from Mexico, and a single individual of C. sasai from Costa Rica. We explore gene expression variation in Cerrophidion and sequence evolution of toxins within C. godmani specifically. Cerrophidion venom gland transcriptomes are composed primarily of snake venom metalloproteinases, phospholipase A[Formula: see text]s (PLA[Formula: see text]s), and snake venom serine proteases. Cerrophidion petlalcalensis shows little intraspecific variation; however, C. godmani and C. tzotzilorum differ significantly between geographically isolated populations. Interestingly, intraspecific variation was mostly attributed to expression variation as we did not detect signals of selection within C. godmani toxins. Additionally, we found PLA[Formula: see text]-like myotoxins in all species except C. petlalcalensis, and crotoxin-like PLA[Formula: see text]s in the southern population of C. godmani. Our results demonstrate significant intraspecific venom variation within C. godmani and C. tzotzilorum. The toxins of C. godmani show little evidence of directional selection where variation in toxin sequence is consistent with evolution under a model of mutation-drift equilibrium. Cerrophidion godmani individuals from the southern population may exhibit neurotoxic venom activity given the presence of crotoxin-like PLA[Formula: see text]s; however, further research is required to confirm this hypothesis.

One Size Fits All-Venomics of the Iberian Adder (Vipera seoanei, Lataste 1878) Reveals Low Levels of Venom Variation across Its Distributional Range

European vipers (genus Vipera) are medically important snakes displaying considerable venom variation, occurring at different levels in this group. The presence of intraspecific venom variation, however, remains understudied in several Vipera species. Vipera seoanei is a venomous snake endemic to the northern Iberian Peninsula and south-western France, presenting notable phenotypic variation and inhabiting several diverse habitats across its range. We analysed the venoms of 49 adult specimens of V. seoanei from 20 localities across the species' Iberian distribution. We used a pool of all individual venoms to generate a V. seoanei venom reference proteome, produced SDS-PAGE profiles of all venom samples, and visualised patterns of variation using NMDS. By applying linear regression, we then assessed presence and nature of venom variation between localities, and investigated the effect of 14 predictors (biological, eco-geographic, genetic) on its occurrence. The venom comprised at least 12 different toxin families, of which five (i.e., PLA2, svSP, DI, snaclec, svMP) accounted for about 75% of the whole proteome. The comparative analyses of the SDS-PAGE venom profiles showed them to be remarkably similar across the sampled localities, suggesting low geographic variability. The regression analyses suggested significant effects of biological and habitat predictors on the little variation we detected across the analysed V. seoanei venoms. Other factors were also significantly associated with the presence/absence of individual bands in the SDS-PAGE profiles. The low levels of venom variability we detected within V. seoanei might be the result of a recent population expansion, or of processes other than directional positive selection.

Interactions between physiology and behaviour provide insights into the ecological role of venom in Australian funnel-web spiders: Interspecies comparison

Australian funnel-web spiders are iconic species, characterized as being the most venomous spiders in the world. They are also valued for the therapeutics and natural bioinsecticides potentially hidden in their venom molecules. Although numerous biochemical and molecular structural approaches have tried to determine the factors driving venom complexity, these approaches have not considered behaviour, physiology and environmental conditions collectively, which can play a role in the evolution, complexity, and function of venom components in funnel-webs. This study used a novel interdisciplinary approach to understand the relationships between different behaviours (assessed in different ecological contexts) and morphophysiological variables (body condition, heart rate) that may affect venom composition in four species of Australian funnel-web spiders. We tested defensiveness, huddling behaviour, frequency of climbing, and activity for all species in three ecological contexts: i) predation using both indirect (puff of air) and direct (prodding) stimuli; ii) conspecific tolerance; and iii) exploration of a new territory. We also assessed morphophysiological variables and venom composition of all species. For Hadronyche valida, the expression of some venom components was associated with heart rate and defensiveness during the predation context. However, we did not find any associations between behavioural traits and morphophysiological variables in the other species, suggesting that particular associations may be species-specific. When we assessed differences between species, we found that the species separated out based on the venom profiles, while activity and heart rate are likely more affected by individual responses and microhabitat conditions. This study demonstrates how behavioural and morphophysiological traits are correlated with venom composition and contributes to a broader understanding of the function and evolution of venoms in funnel-web spiders.

A current perspective on snake venom composition and constituent protein families

Snake venoms are heterogeneous mixtures of proteins and peptides used for prey subjugation. With modern proteomics there has been a rapid expansion in our knowledge of snake venom composition, resulting in the venom proteomes of 30% of vipers and 17% of elapids being characterised. From the reasonably complete proteomic coverage of front-fanged snake venom composition (179 species-68 species of elapids and 111 species of vipers), the venoms of vipers and elapids contained 42 different protein families, although 18 were only reported in < 5% of snake species. Based on the mean abundance and occurrence of the 42 protein families, they can be classified into 4 dominant, 6 secondary, 14 minor, and 18 rare protein families. The dominant, secondary and minor categories account for 96% on average of a snake's venom composition. The four dominant protein families are: phospholipase A₂ (PLA₂), snake venom metalloprotease (SVMP), three-finger toxins (3FTx), and snake venom serine protease (SVSP). The six secondary protein families are: L-amino acid oxidase (LAAO), cysteine-rich secretory protein (CRiSP), C-type lectins (CTL), disintegrins (DIS), kunitz peptides (KUN), and natriuretic peptides (NP). Venom variation occurs at all taxonomic levels, including within populations. The reasons for venom variation are complex, as variation is not always associated with geographical variation in diet. The four dominant protein families appear to be the most important toxin families in human envenomation, being responsible for coagulopathy, neurotoxicity, myotoxicity and cytotoxicity. Proteomic techniques can be used to investigate the toxicological profile of a snake venom and hence identify key protein families for antivenom immunorecognition.

Diversity of Phospholipases A2 from Bothrops atrox Snake Venom: Adaptive Advantages for Snakes Compromising Treatments for Snakebite Patients

The evolution of snake venoms resulted in multigene toxin families that code for structurally similar isoforms eventually harboring distinct functions. PLA₂s are dominant toxins in viper venoms, and little is known about the impact of their diversity on human envenomings and neutralization by antivenoms. Here, we show the isolation of three distinct PLA₂s from B. atrox venom. FA1 is a Lys-49 homologue, and FA3 and FA4 are catalytic Asp-49 PLA₂s. FA1 and FA3 are basic myotoxic proteins, while FA4 is an acid non-myotoxic PLA₂. FA3 was the most potent toxin, inducing higher levels of edema, inflammatory nociception, indirect hemolysis, and anticoagulant activity on human, rat, and chicken plasmas. FA4 presented lower anticoagulant activity, and FA1 had only a slight effect on human and rat plasmas. PLA₂s presented differential reactivities with antivenoms, with an emphasis on FA3, which was not recognized or neutralized by the antivenoms used in this study. Our findings reveal the functional and antigenic diversity among PLA₂s from B. atrox venom, highlighting the importance of assessing venom variability for understanding human envenomations and treatment with antivenoms, particularly evident here as the antivenom fails to recognize FA3, the most active multifunctional toxin described.

Bacterial Adaptation to Venom in Snakes and Arachnida

Animal venoms are considered sterile sources of antimicrobial compounds with strong membrane-disrupting activity against multidrug-resistant bacteria. However, venomous bite wound infections are common in developing nations. Investigating the envenomation organ and venom microbiota of five snake and two spider species, we observed venom community structures that depend on the host venomous animal species and evidenced recovery of viable microorganisms from black-necked spitting cobra (Naja nigricollis) and Indian ornamental tarantula (Poecilotheria regalis) venoms. Among the bacterial isolates recovered from N. nigricollis, we identified two venom-resistant, novel sequence types of Enterococcus faecalis whose genomes feature 16 virulence genes, indicating infectious potential, and 45 additional genes, nearly half of which improve bacterial membrane integrity. Our findings challenge the dogma of venom sterility and indicate an increased primary infection risk in the clinical management of venomous animal bite wounds. IMPORTANCE Notwithstanding their 3 to 5% mortality, the 2.7 million envenomation-related injuries occurring annually-predominantly across Africa, Asia, and Latin America-are also major causes of morbidity. Venom toxin-damaged tissue will develop infections in some 75% of envenomation victims, with E. faecalis being a common culprit of disease; however, such infections are generally considered to be independent of envenomation. Here, we provide evidence on venom microbiota across snakes and arachnida and report on the convergent evolution mechanisms that can facilitate adaptation to black-necked cobra venom in two independent E. faecalis strains, easily misidentified by biochemical diagnostics. Therefore, since inoculation with viable and virulence gene-harboring bacteria can occur during envenomation, acute infection risk management following envenomation is warranted, particularly for immunocompromised and malnourished victims in resource-limited settings. These results shed light on how bacteria evolve for survival in one of the most extreme environments on Earth and how venomous bites must be also treated for infections.

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.

Convergent evolution of venom gland transcriptomes across Metazoa

Animals have repeatedly evolved specialized organs and anatomical structures to produce and deliver a mixture of potent bioactive molecules to subdue prey or predators-venom. This makes it one of the most widespread, convergent functions in the animal kingdom. Whether animals have adopted the same genetic toolkit to evolved venom systems is a fascinating question that still eludes us. Here, we performed a comparative analysis of venom gland transcriptomes from 20 venomous species spanning the main Metazoan lineages to test whether different animals have independently adopted similar molecular mechanisms to perform the same function. We found a strong convergence in gene expression profiles, with venom glands being more similar to each other than to any other tissue from the same species, and their differences closely mirroring the species phylogeny. Although venom glands secrete some of the fastest evolving molecules (toxins), their gene expression does not evolve faster than evolutionarily older tissues. We found 15 venom gland-specific gene modules enriched in endoplasmic reticulum stress and unfolded protein response pathways, indicating that animals have independently adopted stress response mechanisms to cope with mass production of toxins. This, in turn, activates regulatory networks for epithelial development, cell turnover, and maintenance, which seem composed of both convergent and lineage-specific factors, possibly reflecting the different developmental origins of venom glands. This study represents a first step toward an understanding of the molecular mechanisms underlying the repeated evolution of one of the most successful adaptive traits in the animal kingdom.

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.

Climate change and the increase of human population will threaten conservation of Asian cobras

Asian cobras (genus Naja) are venomous snakes distributed from the Middle East to Southeast Asia. Because cobras often live near humans settlements, they are responsible for a large part of snakebite incidents and as such pose a challenge for public health systems. In the light of growing human populations, correctly mapping the present and future ranges of Asian cobras is therefore important for both biological conservation and public health management. Here, we mapped the potential climatic niches of ten Asian cobra species for both the present and the future, with the aim to quantify changes in climate and human population densities relative to their current and future ranges. Our analyses reveal that cobras that are adapted to dry climates and inhabit islands have narrow climatic niches, while those of mainland species with larger geographic ranges are much wider. We also found a higher degree of fragmentation of future cobra distributions; within the next 50 years, Asian cobras will lose an average of around 60% of their current suitable climatic range. In the near future, Naja mandalayensis, N. sputatrix, N. samarensis, and N. philippinensis are likely to have no accessible suitable climate space left. Besides, a further increase of human populations in this region may also exponentially accelerate the effects of anthropogenic impacts. Solutions for conservation may involve awareness and appropriate use of law to overcome the rate of habitat degradation and the increase of animal trade of Asian cobras, while promoting investment on health systems to avoid snakebite fatalities.

Addressing the global snakebite crisis with geo-spatial analyses - Recent advances and future direction

Venomous snakebite is a neglected tropical disease that annually leads to hundreds of thousands of deaths or long-term physical and mental ailments across the developing world. Insufficient data on spatial variation in snakebite risk, incidence, human vulnerability, and accessibility of medical treatment contribute substantially to ineffective on-ground management. There is an urgent need to collect data, fill knowledge gaps and address on-ground management problems. The use of novel, and transdisciplinary approaches that take advantage of recent advances in spatio-temporal models, 'big data', high performance computing, and fine-scale spatial information can add value to snakebite management by strategically improving our understanding and mitigation capacity of snakebite. We review the background and recent advances on the topic of snakebite related geospatial analyses and suggest avenues for priority research that will have practical on-ground applications for snakebite management and mitigation. These include streamlined, targeted data collection on snake distributions, snakebites, envenomings, venom composition, health infrastructure, and antivenom accessibility along with fine-scale models of spatio-temporal variation in snakebite risk and incidence, intraspecific venom variation, and environmental change modifying human exposure. These measures could improve and 'future-proof' antivenom production methods, antivenom distribution and stockpiling systems, and human-wildlife conflict management practices, while simultaneously feeding into research on venom evolution, snake taxonomy, ecology, biogeography, and conservation.

Elucidating the Venom Diversity in Sri Lankan Spectacled Cobra (Naja naja) through De Novo Venom Gland Transcriptomics, Venom Proteomics and Toxicity Neutralization

Inadequate effectiveness of Indian antivenoms in treating envenomation caused by the Spectacled Cobra/Indian Cobra (Naja naja) in Sri Lanka has been attributed to geographical variations in the venom composition. This study investigated the de novo venom-gland transcriptomics and venom proteomics of the Sri Lankan N. naja (NN-SL) to elucidate its toxin gene diversity and venom variability. The neutralization efficacy of a commonly used Indian antivenom product in Sri Lanka was examined against the lethality induced by NN-SL venom in mice. The transcriptomic study revealed high expression of 22 toxin genes families in NN-SL, constituting 46.55% of total transcript abundance. Three-finger toxins (3FTX) were the most diversely and abundantly expressed (87.54% of toxin gene expression), consistent with the dominance of 3FTX in the venom proteome (72.19% of total venom proteins). The 3FTX were predominantly S-type cytotoxins/cardiotoxins (CTX) and α-neurotoxins of long-chain or short-chain subtypes (α-NTX). CTX and α-NTX are implicated in local tissue necrosis and fatal neuromuscular paralysis, respectively, in envenomation caused by NN-SL. Intra-species variations in the toxin gene sequences and expression levels were apparent between NN-SL and other geographical specimens of N. naja, suggesting potential antigenic diversity that impacts antivenom effectiveness. This was demonstrated by limited potency (0.74 mg venom/ml antivenom) of the Indian polyvalent antivenom (VPAV) in neutralizing the NN-SL venom. A pan-regional antivenom with improved efficacy to treat N. naja envenomation is needed.

Phenotypic and functional variation in venom and venom resistance of two sympatric rattlesnakes and their prey

Predator-prey interactions often lead to the coevolution of adaptations associated with avoiding predation and, for predators, overcoming those defences. Antagonistic coevolutionary relationships are often not simple interactions between a single predator and prey but rather a complex web of interactions between multiple coexisting species. Coevolution between venomous rattlesnakes and small mammals has led to physiological venom resistance in several mammalian taxa. In general, viperid venoms contain large quantities of snake venom metalloproteinase toxins (SVMPs), which are inactivated by SVMP inhibitors expressed in resistant mammals. We explored variation in venom chemistry, SVMP expression, and SVMP resistance across four co-distributed species (California Ground Squirrels, Bryant's Woodrats, Southern Pacific Rattlesnakes, and Red Diamond Rattlesnakes) collected from four different populations in Southern California. Our aim was to understand phenotypic and functional variation in venom and venom resistance in order to compare coevolutionary dynamics of a system involving two sympatric predator-prey pairs to past studies that have focused on single pairs. Proteomic analysis of venoms indicated that these rattlesnakes express different phenotypes when in sympatry, with Red Diamonds expressing more typical viperid venom (with a diversity of SVMPs) and Southern Pacifics expressing a more atypical venom with a broader range of non-enzymatic toxins. We also found that although blood sera from both mammals were generally able to inhibit SVMPs from both rattlesnake species, inhibition depended strongly on the snake population, with snakes from one geographic site expressing SVMPs to which few mammals were resistant. Additionally, we found that Red Diamond venom, rather than woodrat resistance, was locally adapted. Our findings highlight the complexity of coevolutionary relationships between multiple predators and prey that exhibit similar offensive and defensive strategies in sympatry.

A phylogeny-aware approach reveals unexpected venom components in divergent lineages of cone snails

Marine gastropods of the genus Conus are renowned for their remarkable diversity and deadly venoms. While Conus venoms are increasingly well studied for their biomedical applications, we know surprisingly little about venom composition in other lineages of Conidae. We performed comprehensive venom transcriptomic profiling for Conasprella coriolisi and Pygmaeconus traillii, first time for both respective genera. We complemented reference-based transcriptome annotation by a de novo toxin prediction guided by phylogeny, which involved transcriptomic data on two additional 'divergent' cone snail lineages, Profundiconus, and Californiconus. We identified toxin clusters (SSCs) shared among all or some of the four analysed genera based on the identity of the signal region-a molecular tag present in toxins. In total, 116 and 98 putative toxins represent 29 and 28 toxin gene superfamilies in Conasprella and Pygmaeconus, respectively; about quarter of these only found by semi-manual annotation of the SSCs. Two rare gene superfamilies, originally identified from fish-hunting cone snails, were detected outside Conus rather unexpectedly, so we further investigated their distribution across Conidae radiation. We demonstrate that both these, in fact, are ubiquitous in Conidae, sometimes with extremely high expression. Our findings demonstrate how a phylogeny-aware approach circumvents methodological caveats of similarity-based transcriptome annotation.

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.

Experimental Manipulation of Corticosterone Does Not Affect Venom Composition or Functional Activity in Free-Ranging Rattlesnakes

AbstractVenom is an integral feeding trait in many animal species. Although venom often varies ontogenetically, little is known about the proximate physiological mediators of venom variation within individuals. The glucocorticoid hormone corticosterone (CORT) can alter the transcription and activation of proteins, including homologues of snake venom components such as snake venom metalloproteinases (SVMPs) and phospholipase A₂ (PLA₂). CORT is endogenously produced by snakes, varies seasonally and also in response to stress, and is a candidate endogenous mediator of changes in venom composition and functional activity. Here, we tested the hypothesis that CORT induces changes in snake venom by sampling the venom of wild adult rattlesnakes before and after they were treated with either empty (control) or CORT-filled (treatment) Silastic implants. We measured longitudinal changes in whole-venom composition, whole-venom total protein content, and enzymatic activity of SVMP and PLA₂ components of venom. We also assessed the within-individual repeatability of venom components. Despite successfully elevating plasma CORT in the treatment group, we found no effect of CORT treatment or average plasma CORT level on any venom variables measured. Except for total protein content, venom components were highly repeatable within individuals ([Formula: see text]). Our results indicate that the effects of CORT, a hormone commonly associated with stress and metabolic functions, in adult rattlesnake venom are negligible. Our findings bode well for venom researchers and biomedical applications that rely on the consistency of venoms produced from potentially stressed individuals and provide an experimental framework for future studies of proximate mediators of venom variation across an individual's life span.

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.

Unexpected lack of specialisation in the flow properties of spitting cobra venom

Venom spitting is a defence mechanism based on airborne venom delivery used by a number of different African and Asian elapid snake species ('spitting cobras'; Naja spp. and Hemachatus spp.). Adaptations underpinning venom spitting have been studied extensively at both behavioural and morphological level in cobras, but the role of the physical properties of venom itself in its effective projection remains largely unstudied. We hereby provide the first comparative study of the physical properties of venom in spitting and non-spitting cobras. We measured the viscosity, protein concentration and pH of the venom of 13 cobra species of the genus Naja from Africa and Asia, alongside the spitting elapid Hemachatus haemachatus and the non-spitting viper Bitis arietans By using published microCT scans, we calculated the pressure required to eject venom through the fangs of a spitting and a non-spitting cobra. Despite the differences in the modes of venom delivery, we found no significant differences between spitters and non-spitters in the rheological and physical properties of the studied venoms. Furthermore, all analysed venoms showed a Newtonian flow behaviour, in contrast to previous reports. Although our results imply that the evolution of venom spitting did not significantly affect venom viscosity, our models of fang pressure suggests that the pressure requirements to eject venom are lower in spitting cobras than in non-spitting cobras.

Biogeographic venom variation in Russell's viper (Daboia russelii) and the preclinical inefficacy of antivenom therapy in snakebite hotspots

BACKGROUND

Snakebite in India results in over 58,000 fatalities and a vast number of morbidities annually. The majority of these clinically severe envenomings are attributed to Russell's viper (Daboia russelii), which has a near pan-India distribution. Unfortunately, despite its medical significance, the influence of biogeography on the composition and potency of venom from disparate D. russelii populations, and the repercussions of venom variation on the neutralisation efficacy of marketed Indian antivenoms, remain elusive.

METHODS

Here, we employ an integrative approach comprising proteomic characterisation, biochemical analyses, pharmacological assessment, and venom toxicity profiling to elucidate the influence of varying ecology and environment on the pan-Indian populations of D. russelii. We then conducted in vitro venom recognition experiments and in vivo neutralisation assays to evaluate the efficacy of the commercial Indian antivenoms against the geographically disparate D. russelii populations.

FINDINGS

We reveal significant intraspecific variation in the composition, biochemical and pharmacological activities and potencies of D. russelii venoms sourced from five distinct biogeographic zones across India. Contrary to our understanding of the consequences of venom variation on the effectiveness of snakebite therapy, commercial antivenom exhibited surprisingly similar neutralisation potencies against the majority of the investigated populations, with the exception of low preclinical efficacy against the semi-arid population from northern India. However, the ability of Indian antivenoms to counter the severe morbid effects of Daboia envenoming remains to be evaluated.

CONCLUSION

The concerning lack of antivenom efficacy against the north Indian population of D. russelii, as well as against two other 'big four' snake species in nearby locations, underscores the pressing need to develop pan-India effective antivenoms with improved efficacy in high snakebite burden locales.

Biogeographical venom variation in the Indian spectacled cobra (Naja naja) underscores the pressing need for pan-India efficacious snakebite therapy

Background

Snake venom composition is dictated by various ecological and environmental factors, and can exhibit dramatic variation across geographically disparate populations of the same species. This molecular diversity can undermine the efficacy of snakebite treatments, as antivenoms produced against venom from one population may fail to neutralise others. India is the world’s snakebite hotspot, with 58,000 fatalities and 140,000 morbidities occurring annually. Spectacled cobra (Naja naja) and Russell’s viper (Daboia russelii) are known to cause the majority of these envenomations, in part due to their near country-wide distributions. However, the impact of differing ecologies and environment on their venom compositions has not been comprehensively studied.

Methods

Here, we used a multi-disciplinary approach consisting of venom proteomics, biochemical and pharmacological analyses, and in vivo research to comparatively analyse N. naja venoms across a broad region (>6000 km; seven populations) covering India’s six distinct biogeographical zones.

Findings

By generating the most comprehensive pan-Indian proteomic and toxicity profiles to date, we unveil considerable differences in the composition, pharmacological effects and potencies of geographically-distinct venoms from this species and, through the use of immunological assays and preclinical experiments, demonstrate alarming repercussions on antivenom therapy. We find that commercially-available antivenom fails to effectively neutralise envenomations by the pan-Indian populations of N. naja, including a complete lack of neutralisation against the desert Naja population.

Conclusion

Our findings highlight the significant influence of ecology and environment on snake venom composition and potency, and stress the pressing need to innovate pan-India effective antivenoms to safeguard the lives, limbs and livelihoods of the country’s 200,000 annual snakebite victims.

Annually, India is burdened by the highest number of snake envenomations across the globe, with over 58,000 fatalities and three times the number of morbidities, predominantly affecting the rural agrarian communities. The spectacled cobra (Naja naja) and Russell’s viper (Daboia russelii) are responsible for the vast majority of envenomations in the country, in part, due to their near country-wide distributions. In this study, we unveil the astounding differences in venom composition of N. naja from six different biogeographical zones across the country (>6000 km). We provide a comprehensive account of their disparate venom proteomic profiles, biochemical and pharmacological effects, and the associated potencies. Our study uncovers alarming differences in the efficacy of the marketed polyvalent antivenoms in neutralising these venoms, thereby, emphasising the pressing need to develop dose-efficacious and pan-India effective antivenoms for the treatment of snakebites in the country. This study also highlights the significant influence of ecology and diverse environments on the venom variability, insinuating the necessity for innovating cost-effective and pan-India efficacious solutions to safeguard the lives, limbs and livelihoods of India’s two hundred thousand annual snakebite victims.

The Curious Case of the "Neurotoxic Skink": Scientific Literature Points to the Absence of Venom in Scincidae

In contrast to the clearly documented evolution of venom in many animal lineages, the origin of reptilian venom is highly debated. Historically, venom has been theorised to have evolved independently in snakes and lizards. However, some of the recent works have argued for the common origin of venom in “Toxicofera” reptiles, which include the order Serpentes (all snakes), and Anguimorpha and Iguania lizards. Nevertheless, in both these contrasting hypotheses, the lizards of the family Scincidae are considered to be harmless and devoid of toxic venoms. Interestingly, an unusual clinical case claiming neurotoxic envenoming by a scincid lizard was recently reported in Southern India. Considering its potentially significant medicolegal, conservation and evolutionary implications, we have summarised the scientific evidence that questions the validity of this clinical report. We argue that the symptoms documented in the patient are likely to have resulted from krait envenomation, which is far too frequent in these regions.

Extensive Variation in the Activities of Pseudocerastes and Eristicophis Viper Venoms Suggests Divergent Envenoming Strategies Are Used for Prey Capture

Snakes of the genera Pseudocerastes and Eristicophis (Viperidae: Viperinae) are known as the desert vipers due to their association with the arid environments of the Middle East. These species have received limited research attention and little is known about their venom or ecology. In this study, a comprehensive analysis of desert viper venoms was conducted by visualising the venom proteomes via gel electrophoresis and assessing the crude venoms for their cytotoxic, haemotoxic, and neurotoxic properties. Plasmas sourced from human, toad, and chicken were used as models to assess possible prey-linked venom activity. The venoms demonstrated substantial divergence in composition and bioactivity across all experiments. Pseudocerastes urarachnoides venom activated human coagulation factors X and prothrombin and demonstrated potent procoagulant activity in human, toad, and chicken plasmas, in stark contrast to the potent neurotoxic venom of P. fieldi. The venom of E. macmahonii also induced coagulation, though this did not appear to be via the activation of factor X or prothrombin. The coagulant properties of P. fieldi and P. persicus venoms varied among plasmas, demonstrating strong anticoagulant activity in the amphibian and human plasmas but no significant effect in that of bird. This is conjectured to reflect prey-specific toxin activity, though further ecological studies are required to confirm any dietary associations. This study reinforces the notion that phylogenetic relatedness of snakes cannot readily predict venom protein composition or function. The significant venom variation between these species raises serious concerns regarding antivenom paraspecificity. Future assessment of antivenom is crucial.

Variation in venom composition in the Australian funnel-web spiders Hadronyche valida

Mygalomorph venom properties and active components, which have importance in medicine, agronomy, venomics, ecology and evolution, have been widely studied, but only a small fraction have been characterised. Several studies have shown inter-individual variation in the composition of venom peptides based on ontogeny, sexual dimorphism, season and diet. However, intra-individual variation in venom composition, which could play a key role in the evolution, diversification and function of toxins, is poorly understood. In this study, we demonstrate significant intra- and inter-individual variation in venom composition in the Australian funnel-web spider Hadronyche valida, highlighting that individuals show different venom profiles over time. Fourteen (four juvenile and ten adult females) funnel-web spiders, maintained under the same environmental conditions and diet, were milked a total of four times, one month apart. We then used reversed-phase high performance liquid chromatography/electrospray ionisation mass spectrometry to generate venom fingerprints containing the retention time and molecular weights of the different toxin components in the venom. Across all individuals, we documented a combined total of 83 individual venom components. Only 20% of these components were shared between individuals. Individuals showed variation in the composition of venom peptides, with some components consistently present over time, while others were only present at specific times. When individuals were grouped using the Jaccard clustering index and Kernel Principal Component Analysis, spiders formed two distinct clusters, most likely due to their origin or time of collection. This study contributes to the understanding of variation in venom composition at different levels (intra-individual, and intra- and inter-specific) and considers some of the mechanisms of selection that may contribute to venom diversification within arachnids. In addition, inter-specific variation in venom composition can be highly useful as a chemotaxonomic marker to identify funnel-web species.

Integrating Top-Down and Bottom-Up Mass Spectrometric Strategies for Proteomic Profiling of Iranian Saw-Scaled Viper, Echis carinatus sochureki, Venom

Saw-scaled or carpet vipers (genus Echis) are considered to cause a higher global snakebite mortality than any other snake. Echis carinatus sochureki (ECS) is a widely distributed snake species, also found across the thirteen provinces of Iran, where it is assumed to be responsible for the most snakebite envenomings. Here, we collected the Iranian specimens of ECS from three different geographically distinct populations, investigated food habits, and performed toxicity assessment and venom proteome profiling to better understand saw-scaled viper life. Our results show that the prey items most commonly found in all populations were arthropods, with scorpions from the family Buthidae particularly well represented. LD₅₀ (median lethal dose) values of the crude venom demonstrate highly comparable venom toxicities in mammals. Consistent with this finding, venom characterization via top-down and bottom-up proteomics, applied to both crude venoms and size-exclusion chromatographic fractions, revealed highly comparable venom compositions among the different populations. By combining all proteomics data, we identified 22 protein families from 102 liquid chromatography and tandem mass spectrometry (LC-MS/MS) raw files, including the most abundant snake venom metalloproteinases (SVMPs, 29-34%); phospholipase A2 (PLA2s, 26-31%); snake venom serine proteinases (SVSPs, 11-12%); l-amino acid oxidases (LAOs, 8-11%), C-type lectins/lectin-like (CTLs, 7-9%) protein families, and many newly detected ones, e.g., renin-like aspartic proteases (RLAPs), fibroblast growth factors (FGFs), peptidyl-prolyl cis-trans isomerases (PPIs), and venom vasodilator peptides (VVPs). Furthermore, we identified and characterized methylated, acetylated, and oxidized proteoforms relating to the PLA2 and disintegrin toxin families and the site of their modifications. It thus seems that post-translational modifications (PTMs) of toxins, particularly target lysine residues, may play an essential role in the structural and functional properties of venom proteins and might be able to influence the therapeutic response of antivenoms, to be investigated in future studies.

Geographic variation of individual venom profile of Crotalus durissus snakes

Background:

South American rattlesnakes are represented in Brazil by a single species, Crotalus durissus, which has public health importance due to the severity of its envenomation and to its wide geographical distribution. The species is subdivided into several subspecies, but the current classification is controversial. In Brazil, the venoms of C. d. terrificus and C. d. collilineatus are used for hyperimmunization of horses for antivenom production, even though the distinction of these two subspecies are mostly by their geographical distribution. In this context, we described a comparative compositional and functional characterization of individual C. d. collilineatus and C. d. terrificus venoms from three Brazilian states.

Methods:

We compared the compositional patterns of C. d. terrificus and C. d. collilineatus individual venoms by 1-DE and RP-HPLC. For functional analyzes, the enzymatic activities of PLA₂, LAAO, and coagulant activity were evaluated. Finally, the immunorecognition of venom toxins by the crotalic antivenom produced at Butantan Institute was evaluated using Western blotting.

Results:

The protein profile of individual venoms from C. d. collilineatus and C. d. terrificus showed a comparable overall composition, despite some intraspecific variation, especially regarding crotamine and LAAO. Interestingly, HPLC analysis showed a geographic pattern concerning PLA₂. In addition, a remarkable intraspecific variation was also observed in PLA₂, LAAO and coagulant activities. The immunorecognition pattern of individual venoms from C. d. collilineatus and C. d. terrificus by crotalic antivenom produced at Butantan Institute was similar.

Conclusions:

The results highlighted the individual variability among the venoms of C. durissus ssp. specimens. Importantly, our data point to a geographical variation of C. durissus ssp. venom profile, regardless of the subspecies, as evidenced by PLA₂ isoforms complexity, which may explain the increase in venom neurotoxicity from Northeastern through Southern Brazil reported for the species.

Snake Venom Metalloproteinases (SVMPs): A structure-function update

Snake venom metalloproteinases (SVMPs) represent a diverse group of multi-domain proteins with several biological activities such as the ability to induce hemorrhage, proteolytic degradation of fibrinogen and fibrin, induction of apoptosis and inhibition of platelet aggregation. Due to these activities, SVMPs are responsible for many of the well-known pathological phenotypes in snake envenomations caused particularly by species from the Viperidae family and the Crotalinae subfamily. These proteins have been classified based on their size and domain structure into P–I, P-II and P-III classes. Comparatively, members of the P–I SVMPs possess the simplest structures, formed by the catalytic metalloproteinase domain only; the P-II SVMPs are moderately more complex, having the canonical disintegrin domain in addition to the metalloproteinase domain; members of the P-III class are more structurally varied, comprising the metalloproteinase, disintegrin-like, and cysteine-rich domains. Proteolytic cleavage, repeated domain loss and presence of other ancillary domains are responsible for structural diversities in the P-III class. However, studies continue to unveil the relationship between the structure and function of these proteins. In this review, we recovered evidences from literature on the structural peculiarities and functional classification of Snake Venom Metalloproteinases. In addition, we reflect on diversities that exist among each class while taking into account specific and up-to-date class-based activities.

Causes and Consequences of Snake Venom Variation

Snake venoms are mixtures of toxins that vary extensively between and within snake species. This variability has serious consequences for the management of the world’s 1.8 million annual snakebite victims. Advances in ‘omic’ technologies have empowered toxinologists to comprehensively characterize snake venom compositions, unravel the molecular mechanisms that underpin venom variation, and elucidate the ensuing functional consequences. In this review, we describe how such mechanistic processes have resulted in suites of toxin isoforms that cause diverse pathologies in human snakebite victims and we detail how variation in venom composition can result in treatment failure. Finally, we outline current therapeutic approaches designed to circumvent venom variation and deliver next-generation treatments for the world’s most lethal neglected tropical disease.

Interspecific and intraspecific venom enzymatic variation among cobras (Naja sp. and Ophiophagus hannah)

The genera Ophiophagus and Naja comprise part of a clade of snakes referred to as cobras, dangerously venomous front-fanged snakes in the family Elapidae responsible for significant human mortality and morbidity throughout Asia and Africa. We evaluated venom enzyme variation for eleven cobra species and three N. kaouthia populations using SDS-PAGE venom fingerprinting and numerous enzyme assays. Acetylcholinesterase and PLA₂ activities were the most variable between species, and PLA₂ activity was significantly different between Malaysian and Thailand N. kaouthia populations. Venom metalloproteinase activity was low and significantly different among most species, but levels were identical for N. kaouthia populations; minor variation in venom L-amino acid oxidase and phosphodiesterase activities were seen between cobra species. Naja siamensis venom lacked the α-fibrinogenolytic activity common to other cobra venoms. In addition, venom from N. siamensis had no detectable metalloproteinase activity and exhibited an SDS-PAGE profile with reduced abundance of higher mass proteins. Venom profiles from spitting cobras (N. siamensis, N. pallida, and N. mossambica) exhibited similar reductions in higher mass proteins, suggesting the evolution of venoms of reduced complexity and decreased enzymatic activity among spitting cobras. Generally, the venom proteomes of cobras show highly abundant three-finger toxin diversity, followed by large quantities of PLA₂s. However, PLA₂ bands and activity were very reduced for N. haje, N. annulifera and N. nivea. Venom compositionalenzy analysis provides insight into the evolution, diversification and distribution of different venom phenotypes that complements venomic data, and this information is critical for the development of effective antivenoms and snakebite treatment.

Toxin expression in snake venom evolves rapidly with constant shifts in evolutionary rates

Key innovations provide ecological opportunity by enabling access to new resources, colonization of new environments, and are associated with adaptive radiation. The most well-known pattern associated with adaptive radiation is an early burst of phenotypic diversification. Venoms facilitate prey capture and are widely believed to be key innovations leading to adaptive radiation. However, few studies have estimated their evolutionary rate dynamics. Here, we test for patterns of adaptive evolution in venom gene expression data from 52 venomous snake species. By identifying shifts in tempo and mode of evolution along with models of phenotypic evolution, we show that snake venom exhibits the macroevolutionary dynamics expected of key innovations. Namely, all toxin families undergo shifts in their rates of evolution, likely in response to changes in adaptive optima. Furthermore, we show that rapid-pulsed evolution modelled as a Lévy process better fits snake venom evolution than conventional early burst or Ornstein-Uhlenbeck models. While our results support the idea of snake venom being a key innovation, the innovation of venom chemistry lacks clear mechanisms that would lead to reproductive isolation and thus adaptive radiation. Therefore, the extent to which venom directly influences the diversification process is still a matter of contention.

Defensive Venoms: Is Pain Sufficient for Predator Deterrence?

Pain, though unpleasant, is adaptive in calling an animal’s attention to potential tissue damage. A long list of animals representing diverse taxa possess venom-mediated, pain-inducing bites or stings that work by co-opting the pain-sensing pathways of potential enemies. Typically, such venoms include toxins that cause tissue damage or disrupt neuronal activity, rendering painful stings honest indicators of harm. But could pain alone be sufficient for deterring a hungry predator? Some venomologists have argued “no”; predators, in the absence of injury, would “see through” the bluff of a painful but otherwise benign sting or bite. Because most algogenic venoms are also toxic (although not vice versa), it has been difficult to disentangle the relative contributions of each component to predator deterrence. Southern grasshopper mice (Onychomys torridus) are voracious predators of arthropods, feeding on a diversity of scorpion species whose stings vary in painfulness, including painful Arizona bark scorpions (Centruroides sculpturatus) and essentially painless stripe-tailed scorpions (Paravaejovis spinigerus). Moreover, southern grasshopper mice have evolved resistance to the lethal toxins in bark scorpion venom, rendering a sting from these scorpions painful but harmless. Results from a series of laboratory experiments demonstrate that painful stings matter. Grasshopper mice preferred to prey on stripe-tailed scorpions rather than bark scorpions when both species could sting; the preference disappeared when each species had their stingers blocked. A painful sting therefore appears necessary for a scorpion to deter a hungry grasshopper mouse, but it may not always be sufficient: after first attacking and consuming a painless stripe-tailed scorpion, many grasshopper mice went on to attack, kill, and eat a bark scorpion even when the scorpion was capable of stinging. Defensive venoms that result in tissue damage or neurological dysfunction may, thus, be required to condition greater aversion than venoms causing pain alone.

Fangs for the Memories? A Survey of Pain in Snakebite Patients Does Not Support a Strong Role for Defense in the Evolution of Snake Venom Composition

Animals use venoms for multiple purposes, most prominently for prey acquisition and self-defense. In snakes, venom composition often evolves as a result of selection for optimization for local diet. However, whether selection for a defensive function has also played a role in driving the evolution of venom composition has remained largely unstudied. Here, we use an online survey of snakebite victims to test a key prediction of a defensive function, that envenoming should result in the rapid onset of severe pain. From the analysis of 584 snakebite reports, involving 192 species of venomous snake, we find that the vast majority of bites do not result in severe early pain. Phylogenetic comparative analysis shows that where early pain after a bite evolves, it is often lost rapidly. Our results, therefore, do not support the hypothesis that natural selection for antipredator defense played an important role in the origin of venom or front-fanged delivery systems in general, although there may be intriguing exceptions to this rule.

When one phenotype is not enough: divergent evolutionary trajectories govern venom variation in a widespread rattlesnake species

Understanding the origin and maintenance of phenotypic variation, particularly across a continuous spatial distribution, represents a key challenge in evolutionary biology. For this, animal venoms represent ideal study systems: they are complex, variable, yet easily quantifiable molecular phenotypes with a clear function. Rattlesnakes display tremendous variation in their venom composition, mostly through strongly dichotomous venom strategies, which may even coexist within a single species. Here, through dense, widespread population-level sampling of the Mojave rattlesnake, Crotalus scutulatus, we show that genomic structural variation at multiple loci underlies extreme geographical variation in venom composition, which is maintained despite extensive gene flow. Unexpectedly, neither diet composition nor neutral population structure explain venom variation. Instead, venom divergence is strongly correlated with environmental conditions. Individual toxin genes correlate with distinct environmental factors, suggesting that different selective pressures can act on individual loci independently of their co-expression patterns or genomic proximity. Our results challenge common assumptions about diet composition as the key selective driver of snake venom evolution and emphasize how the interplay between genomic architecture and local-scale spatial heterogeneity in selective pressures may facilitate the retention of adaptive functional polymorphisms across a continuous space.

Reciprocal transplantation of the heterotrophic coral Tubastraea coccinea (Scleractinia: Dendrophylliidae) between distinct habitats did not alter its venom toxin composition

Tubastraea coccinea is an azooxanthellate coral species recorded in the Indian and Atlantic oceans and is presently widespread in the southwestern Atlantic with an alien status for Brazil. T. coccinea outcompete other native coral species by using a varied repertoire of biological traits. For example, T. coccinea has evolved potent venom capable of immobilizing and digesting zooplankton prey. Diversification and modification of venom toxins can provide potential adaptive benefits to individual fitness, yet acquired alteration of venom composition in cnidarians is poorly understood as the adaptive flexibility affecting toxin composition in these ancient lineages has been largely ignored. We used quantitative high-throughput proteomics to detect changes in toxin expression in clonal fragments of specimens collected and interchanged from two environmentally distinct and geographically separate study sites. Unexpectedly, despite global changes in protein expression, there were no changes in the composition and abundance of toxins from coral fragments recovered from either site, and following clonal transplantation between sites. There were also no apparent changes to the cnidome (cnidae) and gross skeletal or soft tissue morphologies of the specimens. These results suggest that the conserved toxin complexity of T. coccinea co-evolved with innovation of the venom delivery system, and its morphological development and phenotypic expression are not modulated by habitat pressures over short periods of time. The adaptive response of the venom trait to specific predatory regimes, however, necessitates further consideration.

Diet Breadth Mediates the Prey Specificity of Venom Potency in Snakes

Venoms are best known for their ability to incapacitate prey. In predatory groups, venom potency is predicted to reflect ecological and evolutionary drivers relating to diet. While venoms have been found to have preyspecific potencies, the role of diet breadth on venom potencies has yet to be tested at large macroecological scales. Here, using a comparative analysis of 100 snake species, we show that the evolution of prey-specific venom potencies is contingent on the breadth of a species' diet. We find that while snake venom is more potent when tested on species closely related to natural prey items, we only find this prey-specific pattern in species with taxonomically narrow diets. While we find that the taxonomic diversity of a snakes' diet mediates the prey specificity of its venom, the species richness of its diet was not found to affect these prey-specific potency patterns. This indicates that the physiological diversity of a species' diet is an important driver of the evolution of generalist venom potencies. These findings suggest that the venoms of species with taxonomically diverse diets may be better suited to incapacitating novel prey species and hence play an important role for species within changing environments.

Alteration of Bumblebee Venom Composition toward Higher Elevation

Venomous animals use venom, a complex biofluid composed of unique mixtures of proteins and peptides, for either predation or defense. Bumblebees, which occur in various habitats due to their unique thermoregulatory properties, mainly use venom for defense. Herein, we conducted an exploratory analysis of the venom composition of a bumblebee species (Bombus pascuorum) along an elevation gradient in the western Swiss Alps using shot-gun proteomic approaches to assess whether their defense mechanism varies along the gradient. The gradient was characterized by high temperatures and low humidity at low elevations and low temperatures and high humidity at high elevations. Venom composition is changing along the elevation gradient, with proteomic variation in the abundances of pain-inducing and allergenic proteins. In particular, the abundance of phospholipase A₂-like, the main component of bumblebee venom, gradually decreases toward higher elevation (lower temperature), suggesting venom alteration and thus a decrease in bumblebee defense towards harsher environments. Larger datasets may complement this study to validate the observed novel trends.

Coevolution of Snake Venom Toxic Activities and Diet: Evidence that Ecological Generalism Favours Toxicological Diversity

Snake venom evolution is typically considered to be predominantly driven by diet-related selection pressures. Most evidence for this is based on lethality to prey and non-prey species and on the identification of prey specific toxins. Since the broad toxicological activities (e.g., neurotoxicity, coagulotoxicity, etc.) sit at the interface between molecular toxinology and lethality, these classes of activity may act as a key mediator in coevolutionary interactions between snakes and their prey. Indeed, some recent work has suggested that variation in these functional activities may be related to diet as well, but previous studies have been limited in geographic and/or taxonomic scope. In this paper, we take a phylogenetic comparative approach to investigate relationships between diet and toxicological activity classes on a global scale across caenophidian snakes, using the clinically oriented database at toxinology.com. We generally find little support for specific prey types selecting for particular toxicological effects except that reptile-feeders are more likely to be neurotoxic. We find some support for endothermic prey (with higher metabolic rates) influencing toxic activities, but differently from previous suggestions in the literature. More broadly, we find strong support for a general effect of increased diversity of prey on the diversity of toxicological effects of snake venom. Hence, we provide evidence that selection pressures on the toxicological activities of snake venom has largely been driven by prey diversity rather than specific types of prey. These results complement and extend previous work to suggest that specific matching of venom characteristics to prey may occur at the molecular level and translate into venom lethality, but the functional link between those two is not constrained to a particular toxicological route.

Rapid and Differential Evolution of the Venom Composition of a Parasitoid Wasp Depending on the Host Strain

Parasitoid wasps rely primarily on venom to suppress the immune response and regulate the physiology of their host. Intraspecific variability of venom protein composition has been documented in some species, but its evolutionary potential is poorly understood. We performed an experimental evolution initiated with the crosses of two lines of Leptopilinaboulardi of different venom composition to generate variability and create new combinations of venom factors. The offspring were maintained for 10 generations on two strains of Drosophila melanogaster differing in resistance/susceptibility to the parental parasitoid lines. The venom composition of individuals was characterized by a semi-automatic analysis of 1D SDS-PAGE electrophoresis protein profiles whose accuracy was checked by Western blot analysis of well-characterized venom proteins. Results made evident a rapid and differential evolution of the venom composition on both hosts and showed that the proteins beneficial on one host can be costly on the other. Overall, we demonstrated the capacity of rapid evolution of the venom composition in parasitoid wasps, important regulators of arthropod populations, suggesting a potential for adaptation to new hosts. Our approach also proved relevant in identifying, among the diversity of venom proteins, those possibly involved in parasitism success and whose role deserves to be deepened.

Biochemical and immunochemical characterization of venoms from snakes of the genus Agkistrodon

In the present work, venoms from five species of the genus Agkistrodon were evaluated in terms of their enzymatic (Phospholipase A2 and caseinolytic) and biological (edema forming, hemorrhagic, procoagulant and lethal) effects. Horses were used to produce monovalent hyperimmune sera against each of three venoms (A. bilineatus, A. contortrix and A. piscivorus) and their neutralizing potency, expressed as Median Effective Dose (ED50), was determined against the venoms of all five species. In terms of PLA2 and caseinolytic activities, all venoms are extremely homogeneous. PLA2 activity is high, while caseinolytic activity is low when in contrast with that of the rattlesnake Crotalus simus. On the other hand, biological activities showed marked interspecific differences, particularly between the species from Mexico and those from the United States. Mexican species displayed higher edema-forming, hemorrhagic and lethal effects than US species, while none of the species studied presented procoagulant activity. All three monovalent hyperimmune sera showed good neutralizing potency against the analyzed venoms. Nonetheless, we observed relevant immunochemical differences among the venoms using ELISA and Western Blot assays. We conclude that the venoms of A. piscivorus (USA) and A. bilineatus would be ideal to use as immunogens for the production of a polyvalent antivenom with good neutralizing potency against the venoms of all the species of the genus.