Proteomic analysis of venom variability and ontogeny across the arboreal palm-pitvipers (genus Bothriechis)

Envenomation by Eyelash Viper Bothriechis schlegelii (Berthold, 1846) in Southwestern Colombia

INTRODUCTION

Bothriechis schlegelii is a Crotaline viperid species of Central America and Northern South America. The characteristics of its envenomation have not been well established. We present clinical characteristics of human cases evaluated and treated in a hospital in southwestern Colombia.

METHODS

We evaluated data from patients who suffered Bothriechis schlegelii envenomation and were seen at Fundación Valle del Lili Hospital, Cali, Colombia between 2011 and 2022.

RESULTS

Eight patients were included, with a median age of 24 years. Snakebites occurred in rural areas. Six (75%) patients were bitten on the upper extremities in relation to the arboreal habits of this animal. The most common symptoms were pain and edema (N = 8, 100%), ecchymoses (N = 2, 25%), and paresthesia (N = 2, 25%). The most common systemic findings were hypofibrinogenemia (N = 8, 100%) and prolonged prothrombin time in five patients (N = 5, 62.5%). All were treated with polyvalent antivenom for Colombian snakes, with a good response and outcome.

CONCLUSIONS

Most bite sites from B. schlegelii were on the upper limbs. All patients had both local manifestations, including edema, pain, and systemic effects with hypofibrinogenemia, but none had systemic bleeding. Every patient received antivenom and had favorable outcomes.

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.

Bibliometric Analysis of Literature in Snake Venom-Related Research Worldwide (1933-2022)

Snake envenomation is a severe economic and health concern affecting countries worldwide. Snake venom carries a wide variety of small peptides and proteins with various immunological and pharmacological properties. A few key research areas related to snake venom, including its applications in treating cancer and eradicating antibiotic-resistant bacteria, have been gaining significant attention in recent years. The goal of the current study was to analyze the global profile of literature in snake venom research. This study presents a bibliometric review of snake venom-related research documents indexed in the Scopus database between 1933 and 2022. The overall number of documents published on a global scale was 2999, with an average annual production of 34 documents. Brazil produced the highest number of documents (n = 729), followed by the United States (n = 548), Australia (n = 240), and Costa Rica (n = 235). Since 1963, the number of publications has been steadily increasing globally. At a worldwide level, antivenom, proteomics, and transcriptomics are growing hot issues for research in this field. The current research provides a unique overview of snake venom research at global level from 1933 through 2022, and it may be beneficial in guiding future research.

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.

Taxon-selective venom variation in adult and neonate Daboia russelii (Russell's Viper), and antivenom efficacy

Major variations in venom composition can occur between juvenile and adult venomous snakes. However, due to logistical constraints, antivenoms are produced using adult venoms in immunising mixtures, possibly resulting in limited neutralisation of juvenile snake venoms. Daboia russelii is one of the leading causes of snakebite death across South Asia. Its venom is potently procoagulant, causing stroke in prey animals but causing in humans consumptive coagulopathy-a net anticoagulant state-and sometimes death resulting from hemorrhage. In this in vitro study, we compared the venom activity of-and antivenom efficacy against-six 2-week-old D. russelii relative to that of their parents. Using a coagulation analyser, we quantified the relative coagulotoxicity of these venoms in human, avian, and amphibian plasma. The overall potency on human plasma was similar across all adult and neonate venoms, and SII (Serum Institute of India) antivenom was equipotent in neutralising these coagulotoxic effects. In addition, all venoms were also similar in their action upon avian plasma. In contrast, the neonate venoms were more potent on amphibian plasma, suggesting amphibians make up a larger proportion of neonate diet than adult diet. A similar venom potency in human and avian plasmas but varying selectivity for amphibian plasma suggests ontogenetic differences in toxin isoforms within the factor X or factor V activating classes, thereby providing a testable hypothesis for future transcriptomics work. By providing insights into the functional venom differences between adult and neonate D. russelii venoms, we hope to inform clinical treatment of patients envenomated by this deadly species and to shed new light on the natural history of these extremely medically important snakes.

Individual variations in the protein profiles and functional activities of the eyelash palm pit-viper (Bothriechis schlegelii) venom from the Colombian southwest region

Bothriechis schlegelii is a venomous snake found in Central and South America, mainly sighted in regions devoted to agriculture. However, in Colombia, little is known about its contribution to the total envenoming cases. Furthermore, there are no reports of the biochemical and functional activities of venoms from the southwest populations, and the differences respecting other populations are unknown. This study analyzed the protein profiles of venom samples obtained from three specimens originating from this region of Colombia using electrophoresis and chromatography. The lethality, edema-induction, hemorrhagic, defibrinating, coagulant, and indirect hemolytic activities were also evaluated. As a result, venoms were composed of proteins with a wide range of molecular weights, most of them below <37 kDa, with differences between male and female electrophoretic and chromatographic profiles. These variations were also observed in the evaluation of venom functional activities such as pro-coagulant, indirect hemolytic, and edema-inducing activities, whereas neither hemorrhagic nor defibrinating activities were detected. These results are also different considering reports with venom samples from other geographical locations, restating the existence of high intraspecific variability in B. schlegelii venoms, which could have relevant pathophysiological and therapeutic implications.

Mutual enlightenment: A toolbox of concepts and methods for integrating evolutionary and clinical toxinology via snake venomics and the contextual stance

Snakebite envenoming is a neglected tropical disease that may claim over 100,000 human lives annually worldwide. Snakebite occurs as the result of an interaction between a human and a snake that elicits either a defensive response from the snake or, more rarely, a feeding response as the result of mistaken identity. Snakebite envenoming is therefore a biological and, more specifically, an ecological problem. Snake venom itself is often described as a "cocktail", as it is a heterogenous mixture of molecules including the toxins (which are typically proteinaceous) responsible for the pathophysiological consequences of envenoming. The primary function of venom in snake ecology is pre-subjugation, with defensive deployment of the secretion typically considered a secondary function. The particular composition of any given venom cocktail is shaped by evolutionary forces that include phylogenetic constraints associated with the snake's lineage and adaptive responses to the snake's ecological context, including the taxa it preys upon and by which it is predated upon. In the present article, we describe how conceptual frameworks from ecology and evolutionary biology can enter into a mutually enlightening relationship with clinical toxinology by enabling the consideration of snakebite envenoming from an "ecological stance". We detail the insights that may emerge from such a perspective and highlight the ways in which the high-fidelity descriptive knowledge emerging from applications of -omics era technologies - "venomics" and "antivenomics" - can combine with evolutionary explanations to deliver a detailed understanding of this multifactorial health crisis.

Clinical syndromes associated with Viperidae family snake envenomation in southwestern Colombia

BACKGROUND

In southwestern Colombia there is a notable variety of snakes that belong to the Viperidae family (vipers). The particular clinical manifestation related to species is poorly reported.

METHODS

Based on a prospective study about envenomation caused by vipers from 2011 to 2019 at the Fundación Valle del Lili Hospital, Cali, in southwest Colombia, we selected cases of admitted patients in which the snakes responsible were fully identified. They were cataloged by clinical syndrome according to prevalent signs (edema-inducing, necrotizing, blister-inducing, procoagulant, anticoagulant or myotoxic) and were related to the species that caused the envenomation.

RESULTS

From a cohort of 53 patients, 21 patients (16 males [72.7%]) with an average age of 35 (3-69) y were included. The syndromes associated with envenomation were anticoagulant and necrotizing effects of Bothrops asper (five patients [22.7%]), blister-inducing and anticoagulant effects of Bothrops rhombeatus (five [22.7%]), anticoagulant effects of Bothrops punctatus (three patients [13.6%]), edema-inducing and anticoagulant effects of Bothriechis schlegelii (five [22.7%]), edema-inducing and myotoxic effects of Bothrocophias colombianus (one [4.5%]), edema-inducing and myotoxic effects of Bothrocophias myersi (one [4.5%]) and edema-inducing effects of Porthidium nasutum (one [4.5%]).

CONCLUSION

In southwestern Colombia there is notable variety in species of snakes belonging to the family Viperidae (vipers) whose envenomation causes various clinical syndromes.

From birth to adulthood: An analysis of the Brazilian lancehead (Bothrops moojeni) venom at different life stages

The Brazilian lancehead (Bothrops moojeni) has a wide distribution in Brazil and represents a serious public health hazard. Previous works reported that the symptoms of snakebites caused by B. moojeni juveniles’ bites were mainly related to coagulation, while those caused by adults’ bites had a more prominent local damage. In this work, we analyzed the venoms of B. moojeni at different life stages to better understand the ontogeny shift in this species. Snakes were grouped by age and sex, and venom pools were formed accordingly. Compositional analyses by one-dimensional electrophoresis (1-DE), chromatography, and mass spectrometry revealed that ontogenetic changes might be mostly related to phospholipase A₂ (PLA₂) and metalloproteases. Regarding the venoms functional aspect, proteolytic, L-amino acid oxidase, PLA₂, and coagulant in vitro activities were assayed, but only the first and the last ones showed age-related changes, with the venom of snakes up to 1 year-old displaying lower proteolytic and higher coagulant activities, while those from 2 years-old onward presented the opposite relation. The venoms of 3 years-old snakes were exceptions to the compositional and functional pattern of adults as both venoms presented profiles similar to neonates. Sex-related differences were observed in specific groups and were not age-related. In vivo experiments (median lethal dose and hemorrhagic activity) were statistically similar between neonates and adults, however we verified that the adult venom killed mice faster comparing to the neonates. All venoms were mostly recognized by the antibothropic serum and displayed similar profiles to 1-DE in western blotting. In conclusion, the Brazilian lancehead venom showed ontogenetic shift in its composition and activities. Furthermore, this change occurred in snakes from 1 to 2 years-old, and interestingly the venom pools from 3 years-old snakes had particular characteristics, which highlights the importance of comprehensive studies to better understand venom variability.

Clinical implications of ontogenetic differences in the coagulotoxic activity of Bothrops jararacussu venoms

Is snake venom activity influenced by size? This is a long-standing question that can have important consequences for the treatment of snake envenomation. Ontogenetic shifts in venom composition are a well-documented characteristic of numerous snake species. Although snake venoms can cause a range of pathophysiological disturbances, establishing the coagulotoxic profiles related to such shifts is a justified approach because coagulotoxicity can be deadly, and its neutralisation is a challenge for current antivenom therapy. Thus, we aimed to assess the coagulotoxicity patterns on plasma and fibrinogen produced by B othrops jararacussu venoms from individuals of different sizes and sex, and the neutralisation potential of SAB (anti bothropic serum produced by Butantan Institute). The use of a metalloproteinase inhibitor (Prinomastat) and a serine proteinase inhibitor (AEBSF) enabled us to determine the toxin class responsible for the observed coagulopathy: activity on plasma was found to be metalloprotease driven, while the activity on fibrinogen is serine protease driven. To further explore differences in venom activity, the activation of Factor X and prothrombin as a function of snake size was also evaluated. All the venoms exhibited a potent procoagulant effect upon plasma and were less potent in their pseudo-procoagulant clotting effect upon fibrinogen. On human plasma, the venoms from smaller snakes produced more rapid clotting than the larger ones. In contrast, the venom activity on fibrinogen had no relation with size or sex. The difference in procoagulant potency was correlated with the bigger snakes being proportionally better neutralized by antivenom due to the lower levels of procoagulant toxins, than the smaller. Thus, while the antivenom ultimately neutralized the venoms, proportionally more would be needed for an equal mass of venom from a small snake than a large one. Similarly, the neutralisation by SAB of the pseudo-procoagulant clotting effects was also correlated with relative potency, with the smaller and bigger snakes being neutralized proportional to potency, but with no correlation to size. Thromboelastography (TEG) tests on human and toad plasma revealed that small snakes' venoms acted quicker than large snakes' venom on both plasmas, with the action upon amphibian plasma consistent with smaller snakes taking a larger proportion of anuran prey than adults. Altogether, the ontogenetic differences regarding coagulotoxic potency and corresponding impact upon relative antivenom neutralisation of snakes with different sizes were shown, underscoring the medical importance of investigating ontogenetic changes in order to provide data crucial for evidence-based design of clinical management strategies.

Structural and bioinformatic analyses of Azemiops venom serine proteases reveal close phylogeographic relationships to pitvipers from eastern China and the New World

The semi-fossil and pit-less Azemiops feae is possibly the most primitive crotalid species. Here, we have cloned and sequenced cDNAs encoding four serine proteases (vSPs) from the venom glands of Chinese A. feae. Full amino-acid sequences of the major vSP (designated as AzKNa) and three minor vSPs (designated as AzKNb, AzKNc and Az-PA) were deduced. Using Protein-BLAST search, the ten most-similar vSPs for each Azemiops vSP have been selected for multiple sequence alignment, and all the homologs are crotalid vSPs. The results suggest that the A. feae vSPs are structurally most like those of eastern-Chinese Gloydius, Viridovipera, Protobothrops and North American pitvipers, and quite different from more-specialized vSPs such as Agkistrodon venom Protein-C activators. The vSPs from Chinese A. feae and those from Vietnamese A. feae show significant sequence variations. AzKNa is acidic and contains six potential N-glycosylation sites and its surface-charge distribution differs greatly from that of AzKNb, as revealed by 3D-modeling. AzKNb and AzKNc do not contain N-glycosylation sites although most of their close homologs contain one or two. Az-PA belongs to the plasminogen-activator subtype with a conserved N²⁰-glycosylation site. The evolution of this subtype of vSPs in Azemiops and related pitvipers has been traced by phylogenetic analysis.

Snake venom color and L-amino acid oxidase: An evidence of long-term captive Crotalus durissus terrificus venom plasticity

The venom color variation of Crotalus durissus terrificus (Cdt) is attributed to the presence of the toxin L-amino acid oxidase (LAAO). During the venom milking routine of Instituto Butantan, we have noticed that most venoms of captive Cdt specimens show a yellowish color, while most venoms of wild specimens are white. Here we describe a comparative analysis of long-term captive (LTC) and recently wild-caught (RWC) Cdt, focusing on LAAO variation. For the identification of LAAO in individual venoms, four different approaches were employed: evaluation of the enzymatic activity, SDS-PAGE, Western blotting, and ELISA. In addition, mass spectrometry analysis was performed using pooled samples. Although some variation among these methodologies was observed, it was possible to notice that the presence of LAAO was significantly higher in the venom of LTC individuals. LAAO was identified in 60-80% LTC specimens and in only 10-12% of RWC specimens. Furthermore, this enzyme accounts for 5.6% of total venom proteins of LTC Cdt pooled venom, while it corresponds to only 0.7% of RWC Cdt pooled venom. These findings strongly suggest that captive maintenance increases the expression of LAAO in Cdt venom.

Gradual and Discrete Ontogenetic Shifts in Rattlesnake Venom Composition and Assessment of Hormonal and Ecological Correlates

Ontogenetic shifts in venom occur in many snakes but establishing their nature as gradual or discrete processes required additional study. We profiled shifts in venom expression from the neonate to adult sizes of two rattlesnake species, the eastern diamondback and the timber rattlesnake. We used serial sampling and venom chromatographic profiling to test if ontogenetic change occurs gradually or discretely. We found evidence for gradual shifts in overall venom composition in six of eight snakes, which sometimes spanned more than two years. Most chromatographic peaks shift gradually, but one quarter shift in a discrete fashion. Analysis of published diet data showed gradual shifts in overall diet composition across the range of body sizes attained by our eight study animals, while the shifts in abundance of different prey classes varied in form from gradual to discrete. Testosterone concentrations were correlated with the change in venom protein composition, but the relationship is not strong enough to suggest causation. Venom research employing simple juvenile versus adult size thresholds may be failing to account for continuous variation in venom composition lifespan. Our results imply that venom shifts represent adaptive matches to dietary shifts and highlight venom for studies of alternative gene regulatory mechanisms.

A brief review on the natural history, venomics and the medical importance of bushmaster (Lachesis) pit viper snakes

Snakes of the genus Lachesis, commonly known as bushmasters, are the largest venomous snakes in the Americas. Because these snakes have their habitats in areas of remote forests they are difficult to find, and consequently there are few studies of Lachesis taxa in their natural ecosystems. Bushmasters are distributed in tropical forest areas of South and Central America. In Brazil they can be found in the Amazon Rainforest and the Atlantic Forest. Despite the low incidence of cases, laquetic envenoming causes severe permanent sequelae due to the high amount of inoculated venom. These accidents are characterized by local pain, hemorrhage and myonecrosis that can be confused with bothropic envenomings. However, victims of Lachesis bites develop symptoms characteristic of Lachesis envenoming, known as vagal syndrome. An important message of this bibliographic synthesis exercise is that, despite having the proteomic profiles of all the taxa of the genus available, very few structure-function correlation studies have been carried out. Therefore the motivation for this review was to fill a gap in the literature on the genus Lachesis, about which there is no recent review. Here we discuss data scattered in a number of original articles published in specialized journals, spanning the evolutionary history and extant phylogeographic distribution of the bushmasters, their venom composition and diet, as well as the pathophysiology of their bites to humans and the biological activities and possible biotechnological applicability of their venom toxins.

Comparative characterization of Viperidae snake venoms from Perú reveals two compositional patterns of phospholipase A2 expression

Snake species within the Bothrops complex (sensu lato) are of medical relevance in Latin America, but knowledge on their venom characteristics is limited, or even unavailable, for some taxa. Perú harbors 17 species of pit vipers, within the genera Bothrops, Bothriechis, Bothrocophias, Porthidium, Crotalus, and Lachesis. This study compared the venoms of twelve species, through chromatographic and electrophoretic profiles, as well as proteolytic and phospholipase A₂ (PLA₂) activities. Also, proteomic profiles were analyzed for nine of the venoms using a shotgun approach. Results unveiled conspicuous differences in the expression of venom PLA₂s among species, six of them presenting scarce levels as judged by RP-HPLC profiles. Since most species within the bothropoid lineage possess venoms with high to intermediate abundances of this protein family, our findings suggest the existence of a phenotypic duality in the expression of venom PLA₂s within the Bothrops (sensu lato) complex. Bothrops barnetti and Bothrocophias andianus venoms, very scarce in PLA₂s, were shown to lack significant myotoxic activity, highlighting that the observed variability in PLA₂ expression bears toxicological correlations with effects attributed to these proteins. Finally, an attempt to identify phylogenetic relationships of bothropoid species from Perú presenting low- or high-PLA₂ venom phenotypes showed an interspersed pattern, thus precluding a simple phylogenetic interpretation of this venom compositional dichotomy.

•

Venoms from 12 viperids of Perú were compared.

•

Conspicuous differences in the expression of PLA₂ were found.

•

Venoms presenting scarce levels of PLA₂ lack myotoxicity.

•

A new phenotypic dichotomy in venom PLA₂ expression is described within Bothrops (sensu lato).

Repurposing Cancer Drugs Batimastat and Marimastat to Inhibit the Activity of a Group I Metalloprotease from the Venom of the Western Diamondback Rattlesnake, Crotalus atrox

Snakebite envenomation causes over 140,000 deaths every year, predominantly in developing countries. As a result, it is one of the most lethal neglected tropical diseases. It is associated with incredibly complex pathophysiology due to the vast number of unique toxins/proteins present in the venoms of diverse snake species found worldwide. Here, we report the purification and functional characteristics of a Group I (PI) metalloprotease (CAMP-2) from the venom of the western diamondback rattlesnake, Crotalus atrox. Its sensitivity to matrix metalloprotease inhibitors (batimastat and marimastat) was established using specific in vitro experiments and in silico molecular docking analysis. CAMP-2 shows high sequence homology to atroxase from the venom of Crotalus atrox and exhibits collagenolytic, fibrinogenolytic and mild haemolytic activities. It exerts a mild inhibitory effect on agonist-induced platelet aggregation in the absence of plasma proteins. Its collagenolytic activity is completely inhibited by batimastat and marimastat. Zinc chloride also inhibits the collagenolytic activity of CAMP-2 by around 75% at 50 μM, while it is partially potentiated by calcium chloride. Molecular docking studies have demonstrated that batimastat and marimastat are able to bind strongly to the active site residues of CAMP-2. This study demonstrates the impact of matrix metalloprotease inhibitors in the modulation of a purified, Group I metalloprotease activities in comparison to the whole venom. By improving our understanding of snake venom metalloproteases and their sensitivity to small molecule inhibitors, we can begin to develop novel and improved treatment strategies for snakebites.

Predicting antibacterial activity from snake venom proteomes

The continued evolution of antibiotic resistance has increased the urgency for new antibiotic development, leading to exploration of non-traditional sources. In particular, snake venom has garnered attention for its potent antibacterial properties. Numerous studies describing snake venom proteomic composition as well as antibiotic efficacy have created an opportunity to synthesize relationships between venom proteomes and their antibacterial properties. Using literature reported values from peer-reviewed studies, our study generated models to predict efficacy given venom protein family composition, snake taxonomic family, bacterial Gram stain, bacterial morphology, and bacterial respiration strategy. We then applied our predictive models to untested snake species with known venom proteomic compositions. Overall, our results provide potential protein families that serve as accurate predictors of efficacy as well as promising organisms in terms of antibacterial properties of venom. The results from this study suggest potential future research trajectories for antibacterial properties in snake venom by offering hypotheses for a variety of taxa.

Comparative compositional and functional analyses of Bothrops moojeni specimens reveal several individual variations

Snake venoms are complex protein mixtures with different biological activities that can act in both their preys and human victims. Many of these proteins play a role in prey capture and in the digestive process of these animals. It is known that some snakes are resistant to the toxicity of their own venom by mechanisms not yet fully elucidated. However, it was observed in the Laboratory of Herpetology of Instituto Butantan that some Bothrops moojeni individuals injured by the same snake species showed mortalities caused by envenoming effects. This study analyzed the biochemical composition of 13 venom and plasma samples from Bothrops moojeni specimens to assess differences in their protein composition. Application of sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) showed distinct venom protein profiles, but very homogeneous plasma profiles. Western Blotting (WB) was performed with plasma samples, which were submitted to incubation with the respective venom. Some individuals showed an immunorecognized band zone around 25 kDa, indicating interaction between the same individual plasma and venom proteins. Crossed-WB assay using non-self-plasma and venom showed that this variability is due to venom protein composition instead of plasma composition. These venoms presented higher caseinolytic, collagenolytic and coagulant activities than the venoms without these regions recognized by WB. Mass spectrometry analyses performed on two individuals revealed that these individuals present, in addition to higher protein concentrations, other exclusive proteins in their composition. When these same two samples were tested in vivo, the results also showed higher lethality in these venoms, but lower hemorrhagic activity than in the venoms without these regions recognized by WB. In conclusion, some Bothrops moojeni specimens differ in venom composition, which may have implications in envenomation. Moreover, the high individual venom variability found in this species demonstrates the importance to work with individual analyses in studies involving intraspecific venom variability and venom evolution.

Proteomic profiling, functional characterization, and immunoneutralization of the venom of Porthidium porrasi, a pitviper endemic to Costa Rica

The genus Porthidium includes nine pitviper species inhabiting Mexico, Central America, and northern South America. Porthidium porrasi is a species endemic to the Southwest of Costa Rica, for which no information on its venom was available. In this study, the proteomic composition and functional activities of P. porrasi venom are described. The most abundant venom proteins were identified as metalloproteinases (36.5%). In descending order of abundance, proteins belonging to the disintegrin, phospholipase A₂, serine proteinase, C-type lectin/lectin-like, vascular endothelial growth factor, Cysteine-rich secretory protein, L-amino acid oxidase, phospholipase B, and phosphodiesterase families were also identified. P. porrasi venom showed a weak lethal potency in mice (10 μg/g body weight by intraperitoneal route), induced marked hemorrhage and edema, and weak myotoxic effect. These in vivo activities, as well as those assayed in vitro (proteolytic and phospholipase A₂ activities) correlated with compositional data. A comparison of P. porrasi venom with those of three other Porthidium species studied to date reveals a generally conserved compositional and functional pattern in this pitviper genus. Importantly, the lethal effect of P. porrasi venom in mice was adequately cross-neutralized by a heterospecific polyvalent antivenom, supporting its use in the treatment of eventual envenomings by this species.

Protease Activity Profiling of Snake Venoms Using High-Throughput Peptide Screening

Snake venom metalloproteinases (SVMPs) and snake venom serine proteinases (SVSPs) are among the most abundant enzymes in many snake venoms, particularly among viperids. These proteinases are responsible for some of the clinical manifestations classically seen in viperid envenomings, including hemorrhage, necrosis, and coagulopathies. The objective of this study was to investigate the enzymatic activities of these proteins using a high-throughput peptide library to screen for the proteinase targets of the venoms of five viperid (Echis carinatus, Bothrops asper, Daboia russelii, Bitis arietans, Bitis gabonica) and one elapid (Naja nigricollis) species of high medical importance. The proteinase activities of these venoms were each tested against 360 peptide substrates, yielding 2160 activity profiles. A nonlinear regression model that accurately described the observed enzymatic activities was fitted to the experimental data, allowing for the comparison of cleavage rates across species. In this study, previously unknown protein targets of snake venom proteinases were identified, potentially implicating novel human and animal proteins that may be involved in the pathophysiology of viper envenomings. The functional relevance of these targets was further evaluated and discussed. These new findings may contribute to our understanding of the clinical manifestations and underlying biochemical mechanisms of snakebite envenoming by viperid species.

Venom of prey-specialized spiders is more toxic to their preferred prey: A result of prey-specific toxins

In specialized predators, a variety of adaptations have evolved to such a level of specificity that they allow very effective exploitation of focal prey. Venom is an essential adaptive trait of predatory venomous species, such as spiders, yet our knowledge of spider venom is incomplete. In agreement with the prey preference hypothesis, we expected that the venom of spider specialists should be more toxic to focal than to alternative prey, because it is composed of prey-specific toxins. Here we used spiders with three types of trophic specializations: specialists that were ant-eating, termite-eating and spider-eating. We compared the efficacy of prey capture of preferred and alternative prey (measured as paralysis latency) with that of related generalists and profiled the venom of the studied species using proteomic methods. We used 22 spider species: six myrmecophagous, two termitophagous, three araneophagous and 11 euryphagous generalist species belonging to different families. We found that ten of the eleven specialist species induced significantly shorter paralysis latency in preferred prey than in alternative prey. Generalists exhibited either similar efficiency on both prey types or slightly higher efficiency on preferred prey. Multivariate analysis of proteomic profiles (peptides and proteins) revealed significant differences between trophic specializations, particularly in peptides. Specialists appear to have venom composed of unique specific compounds as revealed by the multivariate ordination and indicator analysis. These components are likely prey-specific toxins.

Proteomic Deep Mining the Venom of the Red-Headed Krait, Bungarus flaviceps

The use of -omics technologies allows for the characterization of snake venom composition at a fast rate and at high levels of detail. In the present study, we investigated the protein content of Red-headed Krait (Bungarus flaviceps) venom. This analysis revealed a high diversity of snake venom protein families, as evidenced by high-throughput mass spectrometric analysis. We found all six venom protein families previously reported in a transcriptome study of the venom gland of B. flaviceps, including phospholipases A₂ (PLA₂s), Kunitz-type serine proteinase inhibitors (KSPIs), three-finger toxins (3FTxs), cysteine-rich secretory proteins (CRISPs), snaclecs, and natriuretic peptides. A combined approach of automated database searches and de novo sequencing of tandem mass spectra, followed by sequence similarity searches, revealed the presence of 12 additional toxin families. De novo sequencing alone was able to identify 58 additional peptides, and this approach contributed significantly to the comprehensive description of the venom. Abundant protein families comprise 3FTxs (22.3%), KSPIs (19%), acetylcholinesterases (12.6%), PLA₂s (11.9%), venom endothelial growth factors (VEGFs, 8.4%), nucleotidases (4.3%), and C-type lectin-like proteins (snaclecs, 3.3%); an additional 11 toxin families are present at significantly lower concentrations, including complement depleting factors, a family not previously detected in Bungarus venoms. The utility of a multifaceted approach toward unraveling the proteome of snake venoms, employed here, allowed detection of even minor venom components. This more in-depth knowledge of the composition of B. flaviceps venom facilitates a better understanding of snake venom molecular evolution, in turn contributing to more effective treatment of krait bites.

The habu genome reveals accelerated evolution of venom protein genes

Evolution of novel traits is a challenging subject in biological research. Several snake lineages developed elaborate venom systems to deliver complex protein mixtures for prey capture. To understand mechanisms involved in snake venom evolution, we decoded here the ~1.4-Gb genome of a habu, Protobothrops flavoviridis. We identified 60 snake venom protein genes (SV) and 224 non-venom paralogs (NV), belonging to 18 gene families. Molecular phylogeny reveals early divergence of SV and NV genes, suggesting that one of the four copies generated through two rounds of whole-genome duplication was modified for use as a toxin. Among them, both SV and NV genes in four major components were extensively duplicated after their diversification, but accelerated evolution is evident exclusively in the SV genes. Both venom-related SV and NV genes are significantly enriched in microchromosomes. The present study thus provides a genetic background for evolution of snake venom composition.

Venom Ontogeny in the Mexican Lance-Headed Rattlesnake (Crotalus polystictus)

As trophic adaptations, rattlesnake venoms can vary in composition depending on several intrinsic and extrinsic factors. Ontogenetic changes in venom composition have been documented for numerous species, but little is known of the potential age-related changes in many rattlesnake species found in México. In the current study, venom samples collected from adult and neonate Crotalus polystictus from Estado de México were subjected to enzymatic and electrophoretic analyses, toxicity assays (LD₅₀), and MALDI-TOF mass spectrometry, and a pooled sample of adult venom was analyzed by shotgun proteomics. Electrophoretic profiles of adult males and females were quite similar, and only minor sex-based variation was noted. However, distinct differences were observed between venoms from adult females and their neonate offspring. Several prominent bands, including P-I and P-III snake venom metalloproteinases (SVMPs) and disintegrins (confirmed by MS/MS) were present in adult venoms and absent/greatly reduced in neonate venoms. Age-dependent differences in SVMP, kallikrein-like, phospholipase A₂ (PLA₂), and L-amino acid oxidase (LAAO) activity levels were confirmed by enzymatic activity assays, and like many other rattlesnake species, venoms from adult snakes have higher SVMP activity than neonate venoms. Conversely, PLA₂ activity was approximately 2.5 × greater in venoms from neonates, likely contributing to the increased toxicity (neonate venom LD₅₀ = 4.5 μg/g) towards non-Swiss albino mice when compared to adult venoms (LD₅₀ = 5.5 μg/g). Thrombin-like (TLE) and phosphodiesterase activities did not vary significantly with age. A significant effect of sex (between adult male and adult female venoms) was also observed for SVMP, TLE, and LAAO activities. Analysis of pooled adult venom by LC-MS/MS identified 14 toxin protein families, dominated by bradykinin-inhibitory peptides, SVMPs (P-I, P-II and P-III), disintegrins, PLA₂s, C-type-lectins, CRiSPs, serine proteinases, and LAAOs (96% of total venom proteins). Neonate and adult C. polystictus in this population consume almost exclusively mammals, suggesting that age-based differences in composition are related to physical differences in prey (e.g., surface-to-volume ratio differences) rather than taxonomic differences between prey. Venoms from adult C. polystictus fit a Type I pattern (high SVMP activity, lower toxicity), which is characteristic of many larger-bodied rattlesnakes of North America.

Structures of N-Glycans of Bothrops Venoms Revealed as Molecular Signatures that Contribute to Venom Phenotype in Viperid Snakes

The complexity of snake venoms has long been investigated to explore a myriad of biologically active proteins and peptides that are used for immobilizing or killing prey, and are responsible for the pathological effects observed on envenomation. Glycosylation is the main post-translational modification (PTM) of viperid venoms but currently there is little understanding of how protein glycosylation impacts the variation of venom proteomes. We have previously reported that Bothrops venom glycoproteomes contain a core of components that markedly define their composition and parallel their phylogenetic classification. Here we extend those observations to eight Bothrops species evaluating the N-glycomes by LC-MS as assigned cartoon structures and detailing those structures separately as methylated analogs using ion-trap mass spectrometry (MSⁿ). Following ion disassembly through multiple steps provided sequence and linkage isomeric details that characterized 52 unique compositions in Bothrops venoms. These occurred as 60 structures, of which 26 were identified in the venoms of the Jararaca Complex (B. alcatraz, B. insularis, and B. jararaca), 20 in B. erythromelas, B. jararacussu, B. moojeni and B. neuwiedi venoms, and 22 in B. cotiara venom. Further, quantitative analysis of these N-glycans showed variable relative abundances in the venoms. For the first time a comprehensive set of N-glycan structures present in snake venoms are defined. Despite the fact that glycosylation is not template-defined, the N-glycomes of these venoms mirror the phylogeny cladograms of South American bothropoid snakes reported in studies on morphological, molecular data and feeding habits, exhibiting distinct molecular signatures for each venom. Considering the complexity of N-glycan moieties generally found in glycoproteins, characterized by different degrees of branching, isomer structures, and variable abundances, our findings point to these factors as another level of complexity in Bothrops venoms, features that could dramatically contribute to their distinct biological activities.

Comparative analyses of putative toxin gene homologs from an Old World viper, Daboia russelii

Availability of snake genome sequences has opened up exciting areas of research on comparative genomics and gene diversity. One of the challenges in studying snake genomes is the acquisition of biological material from live animals, especially from the venomous ones, making the process cumbersome and time-consuming. Here, we report comparative sequence analyses of putative toxin gene homologs from Russell’s viper (Daboia russelii) using whole-genome sequencing data obtained from shed skin. When compared with the major venom proteins in Russell’s viper studied previously, we found 45–100% sequence similarity between the venom proteins and their putative homologs in the skin. Additionally, comparative analyses of 20 putative toxin gene family homologs provided evidence of unique sequence motifs in nerve growth factor (NGF), platelet derived growth factor (PDGF), Kunitz/Bovine pancreatic trypsin inhibitor (Kunitz BPTI), cysteine-rich secretory proteins, antigen 5, andpathogenesis-related1 proteins (CAP) and cysteine-rich secretory protein (CRISP). In those derived proteins, we identified V11 and T35 in the NGF domain; F23 and A29 in the PDGF domain; N69, K2 and A5 in the CAP domain; and Q17 in the CRISP domain to be responsible for differences in the largest pockets across the protein domain structures in crotalines, viperines and elapids from the in silico structure-based analysis. Similarly, residues F10, Y11 and E20 appear to play an important role in the protein structures across the kunitz protein domain of viperids and elapids. Our study highlights the usefulness of shed skin in obtaining good quality high-molecular weight DNA for comparative genomic studies, and provides evidence towards the unique features and evolution of putative venom gene homologs in vipers.

A Review and Database of Snake Venom Proteomes

Advances in the last decade combining transcriptomics with established proteomics methods have made possible rapid identification and quantification of protein families in snake venoms. Although over 100 studies have been published, the value of this information is increased when it is collated, allowing rapid assimilation and evaluation of evolutionary trends, geographical variation, and possible medical implications. This review brings together all compositional studies of snake venom proteomes published in the last decade. Compositional studies were identified for 132 snake species: 42 from 360 (12%) Elapidae (elapids), 20 from 101 (20%) Viperinae (true vipers), 65 from 239 (27%) Crotalinae (pit vipers), and five species of non-front-fanged snakes. Approximately 90% of their total venom composition consisted of eight protein families for elapids, 11 protein families for viperines and ten protein families for crotalines. There were four dominant protein families: phospholipase A₂s (the most common across all front-fanged snakes), metalloproteases, serine proteases and three-finger toxins. There were six secondary protein families: cysteine-rich secretory proteins, l-amino acid oxidases, kunitz peptides, C-type lectins/snaclecs, disintegrins and natriuretic peptides. Elapid venoms contained mostly three-finger toxins and phospholipase A₂s and viper venoms metalloproteases, phospholipase A₂s and serine proteases. Although 63 protein families were identified, more than half were present in <5% of snake species studied and always in low abundance. The importance of these minor component proteins remains unknown.

Integrated Venomics and Venom Gland Transcriptome Analysis of Juvenile and Adult Mexican Rattlesnakes Crotalus simus, C. tzabcan, and C. culminatus Revealed miRNA-modulated Ontogenetic Shifts

Adult rattlesnakes within genus Crotalus express one of two distinct venom phenotypes, type I (hemorrhagic) and type II (neurotoxic). In Costa Rican Central American rattlesnake, ontogenetic changes in the concentration of miRNAs modulate venom type II to type I transition. Venomics and venom gland transcriptome analyses showed that adult C. simus and C. tzabcan expressed intermediate patterns between type II and type I venoms, whereas C. culminatus had a canonical type I venom. Neonate/juvenile and adult Mexican rattlesnakes showed notable inter- and intraspecific variability in the number, type, abundance and ontogenetic shifts of the transcriptional and translational venom gland activities. These results support a role for miRNAs in the ontogenetic venom compositional changes in the three congeneric Mexican rattlesnakes. It is worth noting the finding of dual-action miRNAs, which silence the translation of neurotoxic heterodimeric PLA₂ crotoxin and acidic PLA₂ mRNAs while simultaneously up-regulating SVMP-targeting mRNAs. Dual transcriptional regulation potentially explains the existence of mutually exclusive crotoxin-rich (type-II) and SVMP-rich (type-I) venom phenotypic dichotomy among rattlesnakes. Our results support the hypothesis that alterations of the distribution of miRNAs, modulating the translational activity of venom gland toxin-encoding mRNAs in response to an external cue, may contribute to the mechanism generating adaptive venom variability.

Venom On-a-Chip: A Fast and Efficient Method for Comparative Venomics

Venom research has attracted an increasing interest in disparate fields, from drug development and pharmacology, to evolutionary biology and ecology, and rational antivenom production. Advances in “-omics” technologies have allowed the characterization of an increasing number of animal venoms, but the methodology currently available is suboptimal for large-scale comparisons of venom profiles. Here, we describe a fast, reproducible and semi-automated protocol for investigating snake venom variability, especially at the intraspecific level, using the Agilent Bioanalyzer on-chip technology. Our protocol generated a phenotype matrix which can be used for robust statistical analysis and correlations of venom variation with ecological correlates, or other extrinsic factors. We also demonstrate the ease and utility of combining on-chip technology with previously fractionated venoms for detection of specific individual toxin proteins. Our study describes a novel strategy for rapid venom discrimination and analysis of compositional variation at multiple taxonomic levels, allowing researchers to tackle evolutionary questions and unveiling the drivers of the incredible biodiversity of venoms.

Third Generation Antivenomics: Pushing the Limits of the In Vitro Preclinical Assessment of Antivenoms

Second generation antivenomics is a translational venomics approach designed to complement in vivo preclinical neutralization assays. It provides qualitative and quantitative information on the set of homologous and heterologous venom proteins presenting antivenom-recognized epitopes and those exhibiting impaired immunoreactivity. In a situation of worrying antivenom shortage in many tropical and sub-tropical regions with high snakebite mortality and morbidity rates, such knowledge has the potential to facilitate the optimal deployment of currently existing antivenoms and to aid in the rational design of novel broad specificity antidotes. The aim of the present work was to expand the analytical capability of the immunoaffinity second-generation antivenomics platform, endowing it with the ability to determine the maximal binding capacity of an antivenom toward the different toxins present in a venom, and to quantify the fraction of venom-specific antibodies present in a given antivenom. The application of this new platform, termed third generation (3G) antivenomics, in the preclinical evaluation of antivenoms is illustrated in this paper for the case of antivenom EchiTAb-Plus-ICP^® reactivity towards the toxins of homologous (B. arietans) and heterologous (N. melanoleuca) venoms.

Strategies in 'snake venomics' aiming at an integrative view of compositional, functional, and immunological characteristics of venoms

This work offers a general overview on the evolving strategies for the proteomic analysis of snake venoms, and discusses how these may be combined through diverse experimental approaches with the goal of achieving a more comprehensive knowledge on the compositional, toxic, and immunological characteristics of venoms. Some recent developments in this field are summarized, highlighting how strategies have evolved from the mere cataloguing of venom components (proteomics/venomics), to a broader exploration of their immunological (antivenomics) and functional (toxicovenomics) characteristics. Altogether, the combination of these complementary strategies is helping to build a wider, more integrative view of the life-threatening protein cocktails produced by venomous snakes, responsible for thousands of deaths every year.

Venomics: integrative venom proteomics and beyond*

Venoms are integrated phenotypes that evolved independently in, and are used for predatory and defensive purposes by, a wide phylogenetic range of organisms. The same principles that contribute to the evolutionary success of venoms, contribute to making the study of venoms of great interest in such diverse fields as evolutionary ecology and biotechnology. Evolution is profoundly contingent, and nature also reinvents itself continuosly. Changes in a complex phenotypic trait, such as venom, reflect the influences of prior evolutionary history, chance events, and selection. Reconstructing the natural history of venoms, particularly those of snakes, which will be dealt with in more detail in this review, requires the integration of different levels of knowledge into a meaningful and comprehensive evolutionary framework for separating stochastic changes from adaptive evolution. The application of omics technologies and other disciplines have contributed to a qualitative and quantitative advance in the road map towards this goal. In this review we will make a foray into the world of animal venoms, discuss synergies and complementarities of the different approaches used in their study, and identify current bottlenecks that prevent inferring the evolutionary mechanisms and ecological constraints that molded snake venoms to their present-day variability landscape.