Diet and snake venom evolution

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

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

Generation of chicken-based IgY polyclonal antibodies against Dendroaspis polylepis and preclinical evaluation of envenomation-neutralizing efficacy vis-à-vis selected commercial antivenoms

The Black mamba, D. polylepis, is one of the many venomous snakes found in Kenya, and known to account for some snakebite incidents. The Kenyan Ministry of Health data reveals annual 15,000 snakebites occurrences. Also, 1 in 15 people in Kenya gets bitten by a snake, and tragically, 1 in 147 of these individuals die of snakebite yearly. Traditionally, antivenoms for treatment are produced from horse or sheep but have complicated and expensive production issues. Alternative production approaches, such as using IgY antibodies derived from chicken egg yolks, may overcome disadvantages with traditional antivenom manufacturing techniques. In this current study, D. polylepis specific IgY polyclonal antibodies were purified from the egg yolks of chickens immunized with D. polylepis venom. These antibodies were subsequently assessed for their in-vivo neutralizing capacity vis-à-vis commercial antivenoms, PANAF-Premium and VINS. The IgY antibodies were purified by ammonium sulfate precipitation and affinity-chromatography, with quality and specificity determined by SDS-PAGE and ELISA. The LD50 of D. polylepis was found to be 0.54 mg/kg in chicks, and 0.34 mg/kg in mice, respectively. Pool of extracted IgY yielded 2.8 mg/mL concentration. Purified IgY under non-reducing and reducing conditions on SDS-PAGE exhibited a single-protein band of about 183 kDa and two bands (67 kDa and 25 kDa), respectively. The minimum-edematogenic dose was 0.05 μg. Anti-D. polylepis IgY antibodies and two antivenoms demonstrated the capacity to neutralize the toxic activities of D. polylepis venom. This study confirms a successful IgY generation against Black mamba venom for the first time, and observed toxic effects of the venom as well as neutralizing capacity of antivenoms.

Assessing the effects of temperature, diet and threat conditions on defensive behaviour and venom regeneration in scorpion (Buthusatlantis)

Understanding animal's behaviour and adaptation in the face of threats and predators under different biotic and abiotic conditions is fundamental in ecology. In this study we examined defensive behaviour of Buthus atlantis scorpion in order to assess how various factors such as temperature, prey type, and threatening conditions influence stinging behaviour, venom usage and regeneration. Our study had revealed that stings frequency was significantly lower in cooler temperature compared to the medium and warm temperature. Threatening condition had no significant effect in medium and warmer temperature, the difference between the two conditions was only significant in the cooler temperature. Conversely, we had shown that venom expenditure in B. atlantis is regulated by both temperature and threatening conditions. Our results show that scorpions maintained in higher temperatures yielded the greatest amount of venom compared to those in lower temperatures. Analyses of proteins concentration according to temperature and diet variation had revealed that scorpions placed in intermediate (25 °C) and warmer temperature (40 °C) had a significantly higher venom proteins concentration when compared to the cooler temperature (10 °C). Results also showed that scorpions adjust their venom usage based on their perception of danger, which can be influenced by temperature.

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.

Host and venom evolution in parasitoid wasps: does independently adapting to the same host shape the evolution of the venom gland transcriptome?

BACKGROUND

Venoms have repeatedly evolved over 100 occasions throughout the animal tree of life, making them excellent systems for exploring convergent evolutionary novelty. Growing evidence supports that venom evolution is predominantly driven by prey or host-related selection pressures, and the expression patterns of venom glands reflect adaptive evolution. However, it remains elusive whether the evolution of expression patterns in venom glands is likewise a convergent evolution driven by their prey/host species.

RESULTS

We utilized parasitoid wasps that had independently adapted to Drosophila hosts as models to investigate the convergent evolution of venom gland transcriptomes in 19 hymenopteran species spanning ~ 200 million years of evolution. Comparative transcriptome analysis reveals that the global expression patterns among the venom glands of Drosophila parasitoid wasps do not achieve higher similarity compared to non-Drosophila parasitoid wasps. Further evolutionary analyses of expression patterns at the single gene, orthogroup, and Gene Ontology (GO) term levels indicate that some orthogroups/GO terms show correlation with the Drosophila parasitoid wasps. However, these groups rarely include genes highly expressed in venom glands or putative venom genes in the Drosophila parasitoid wasps.

CONCLUSIONS

Our study suggests that convergent evolution may not play a predominant force shaping gene expression levels in the venom gland of the Drosophila parasitoid wasps, offering novel insights into the co-evolution between venom and prey/host.

Prey Shifts Drive Venom Evolution in Cone Snails

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

Selection Across the Three-Dimensional Structure of Venom Proteins from North American Scolopendromorph Centipedes

Gene duplication followed by nucleotide differentiation is one of the simplest mechanisms to develop new functions for genes. However, the evolutionary processes underlying the divergence of multigene families remain controversial. We used multigene families found within the diversity of toxic proteins in centipede venom to test two hypotheses related to venom evolution: the two-speed mode of venom evolution and the rapid accumulation of variation in exposed residues (RAVER) model. The two-speed mode of venom evolution proposes that different types of selection impact ancient and younger venomous lineages with negative selection being the predominant form in ancient lineages and positive selection being the dominant form in younger lineages. The RAVER hypothesis proposes that, instead of different types of selection acting on different ages of venomous lineages, the different types of selection will selectively contribute to amino acid variation based on whether the residue is exposed to the solvent where it can potentially interact directly with toxin targets. This hypothesis parallels the longstanding understanding of protein evolution that suggests that residues found within the structural or active regions of the protein will be under negative or purifying selection, and residues that do not form part of these areas will be more prone to positive selection. To test these two hypotheses, we compared the venom of 26 centipedes from the order Scolopendromorpha from six currently recognized species from across North America using both transcriptomics and proteomics. We first estimated their phylogenetic relationships and uncovered paraphyly among the genus Scolopendra and evidence for cryptic diversity among currently recognized species. Using our phylogeny, we then characterized the diverse venom components from across the identified clades using a combination of transcriptomics and proteomics. We conducted selection-based analyses in the context of predicted three-dimensional properties of the venom proteins and found support for both hypotheses. Consistent with the two-speed hypothesis, we found a prevalence of negative selection across all proteins. Consistent with the RAVER hypothesis, we found evidence of positive selection on solvent-exposed residues, with structural and less-exposed residues showing stronger signal for negative selection. Through the use of phylogenetics, transcriptomics, proteomics, and selection-based analyses, we were able to describe the evolution of venom from an ancient venomous lineage and support principles of protein evolution that directly relate to multigene family evolution.

Venom variation among the three subspecies of the North African mountain viper Vipera monticola Saint Girons 1953

The North African mountain viper (Vipera monticola) is a medically relevant venomous snake distributed in Morocco, Algeria, and Tunisia. Three subspecies of V. monticola, exhibiting differences in morphotypes and dietary regimes, are currently recognised: V. m. monticola, V. m. atlantica, and V. m. saintgironsi. Through the application of snake venomics, we analysed the venoms of specimens of Moroccan origin belonging to each of the three subspecies. Snake venom metalloproteinases (svMP), snake venom serine proteases (svSP), C-type lectin and C-type lectin-related proteins (CTL), and phospholipases A2 (PLA2) were predominant, with PLA2 being the most abundant toxin family overall. Disintegrins (DI) and cysteine-rich secretory proteins (CRISP) were exclusive to V. m. monticola and V. m. atlantica, while l-amino-acid oxidases (LAAO) were only found in V. m. saintgironsi. The differences detected in the venom profiles, as well as in presence/absence and relative abundances of toxin families, indicate the occurrence of intraspecific venom variation within V. monticola. The identified patterns of venom similarity between subspecies seem to align more with their phylogenetic relationships than with the reported differences in their feeding habits.

Seasonality in Crotalus durissus venom

Rattlesnakes belonging to the genus Crotalus are widely distributed throughout the Americas. In Brazil, symptoms commonly associated with envenomation by Crotalus durissus collilineatus include myalgia, rhabdomyolysis, renal failure, neurotoxicity, and progressive paralysis, which are related to the protein composition of this venom. Snake venom composition exhibits compositional variability that may reflect geographic distribution, age, captivity, diet, sex, and even individual genetics. Although seasonality is also considered a possible source of variation, there are few reports of such variability in snake venom. In this work, venoms of the same eight C. durissus collilineatus were extracted every three months for two years, to analyze seasonal changes in composition and activities. To this end, venom composition was analyzed by protein quantification, SDS-PAGE, and HPLC, and the LAAO, PLA2 and coagulant activities were measured. Venoms of these C. d. collilineatus showed minor seasonal differences in venom activities and no composition differences were found. LAAO and coagulant activities displayed a pattern of seasonal change, while PLA2 activity seemed to have no seasonality tendency. Also, there are sexual differences, in which males seem to be more stable than females in regard to some activities. Individual variability occurs even in seasonal variation of activities, highlighting the importance of controlling circumstances of venom extraction before comparing results between groups of snakes.

A Comparative Analysis of the Cytotoxic and Vascular Activity Effects of Western Diamondback Rattlesnake (Crotalus atrox) and Eastern Diamondback Rattlesnake (Crotalus adamanteus) Venoms Using a Chick Embryo Model

Crotalus snakebites induce various toxicological effects, encompassing neurological, myotoxic, and cytotoxic symptoms, with potentially fatal outcomes. Investigating venom toxicity is essential for public health, and developing new tools allows for these effects to be studied more comprehensively. The research goals include the elucidation of the physiological consequences of venom exposure and the assessment of toxicity using animal models. Chicken embryos serve as valuable models for assessing venom toxicity through the chick embryotoxicity screening test (CHEST) and the chick chorioallantoic membrane (CAM) assay, particularly useful for evaluating vascular impacts. C. adamanteus venom application resulted in higher embryotoxicity and morphological abnormalities, such as Siamese twins. The CAM assay demonstrated the hemorrhagic effects of venom, varying with venom type and concentration. The irritant potential of both venom types was classified as slight or moderate depending on their concentration. Additionally, acetylcholinesterase (AChE) activity was performed to receive information about organ toxicity. The results show that both venoms induced changes in the whole embryo, heart, and liver weights, but the C. adamanteus venom was identified as more toxic. Specific venom concentrations affected AChE activity in embryonic tissues. These findings underscore the embryotoxic and vasoactive properties of Crotalus venoms, providing valuable insights into their mechanisms of toxicity and potential applications in biomedicine.

Indian green pit vipers: A lesser-known snake group of north-east India

Green pit vipers are one of the most widely distributed group of venomous snakes in south-east Asia. In Indian, green pit vipers are found in the Northern and North-eastern states spreading across eastern and central India and one of the lesser studied venoms. High morphological similarity among them has been a long-established challenge for species identification, however, a total of six species of Indian green pit viper belonging to genus Trimeresurus, Popeia and Viridovipera has been reported from North-east India. Biochemical and biological studies have revealed that venom exhibits substantial variation in protein expression level along with functional variability. The symptoms of envenomation are painful swelling at bite site, bleeding, necrosis along with systemic toxicity such as prolonged coagulopathy. Clinical data of green pit viper envenomated patients from Demow community health centre, Assam advocated against the use of Indian polyvalent antivenom pressing the need for a suitable antivenom for the treatment of green pit viper envenomation. To design effective and specific antivenom for green pit vipers, unveiling the proteome profile of these snakes is needed. In this study, a comparative venomic of green pit vipers of Northern and North-eastern India, their clinical manifestation as well as treatment protocol has been reviewed.

Variability in antivenom neutralization of Mexican viperid snake venoms

BACKGROUND

Each year, 3,800 cases of snakebite envenomation are reported in Mexico, resulting in 35 fatalities. The only scientifically validated treatment for snakebites in Mexico is the use of antivenoms. Currently, two antivenoms are available in the market, with one in the developmental phase. These antivenoms, produced in horses, consist of F(ab')2 fragments generated using venoms from various species as immunogens. While previous studies primarily focused on neutralizing the venom of the Crotalus species, our study aims to assess the neutralization capacity of different antivenom batches against pit vipers from various genera in Mexico.

METHODOLOGY

We conducted various biological and biochemical tests to characterize the venoms. Additionally, we performed neutralization tests using all three antivenoms to evaluate their effectiveness against lethal activity and their ability to neutralize proteolytic and fibrinogenolytic activities.

RESULTS

Our results reveal significant differences in protein content and neutralizing capacity among different antivenoms and even between different batches of the same product. Notably, the venom of Crotalus atrox is poorly neutralized by all evaluated batches despite being the primary cause of envenomation in the country's northern region. Furthermore, even at the highest tested concentrations, no antivenom could neutralize the lethality of Metlapilcoatlus nummifer and Porthidium yucatanicum venoms. These findings highlight crucial areas for improving existing antivenoms and developing new products.

CONCLUSION

Our research reveals variations in protein content and neutralizing potency among antivenoms, emphasizing the need for consistency in venom characteristics as immunogens. While Birmex neutralizes more LD50 per vial, Antivipmyn excels in specific neutralization. The inability of antivenoms to neutralize certain venoms, especially M. nummifer and P. yucatanicum, highlights crucial improvement opportunities, given the medical significance of these species.

Studying Venom Toxin Variation Using Accurate Masses from Liquid Chromatography-Mass Spectrometry Coupled with Bioinformatic Tools

This study provides a new methodology for the rapid analysis of numerous venom samples in an automated fashion. Here, we use LC-MS (Liquid Chromatography-Mass Spectrometry) for venom separation and toxin analysis at the accurate mass level combined with new in-house written bioinformatic scripts to obtain high-throughput results. This analytical methodology was validated using 31 venoms from all members of a monophyletic clade of Australian elapids: brown snakes (Pseudonaja spp.) and taipans (Oxyuranus spp.). In a previous study, we revealed extensive venom variation within this clade, but the data was manually processed and MS peaks were integrated into a time-consuming and labour-intensive approach. By comparing the manual approach to our new automated approach, we now present a faster and more efficient pipeline for analysing venom variation. Pooled venom separations with post-column toxin fractionations were performed for subsequent high-throughput venomics to obtain toxin IDs correlating to accurate masses for all fractionated toxins. This workflow adds another dimension to the field of venom analysis by providing opportunities to rapidly perform in-depth studies on venom variation. Our pipeline opens new possibilities for studying animal venoms as evolutionary model systems and investigating venom variation to aid in the development of better antivenoms.

Foraging predicts the evolution of warning coloration and mimicry in snakes

Warning coloration and Batesian mimicry are classic examples of Darwinian evolution, but empirical evolutionary patterns are often paradoxical. We test whether foraging costs predict the evolution of striking coloration by integrating genetic and ecological data for aposematic and mimetic snakes (Elapidae and Dipsadidae). Our phylogenetic comparison on a total of 432 species demonstrated that dramatic changes in coloration were well predicted by foraging strategy. Multiple tests consistently indicated that warning coloration and conspicuous mimicry were more likely to evolve in species where foraging costs of conspicuous appearance were relaxed by poor vision of their prey, concealed habitat, or nocturnal activity. Reversion to crypsis was also well predicted by ecology for elapids but not for dipsadids. In contrast to a theoretical prediction and general trends, snakes' conspicuous coloration was correlated with secretive ecology, suggesting that a selection regime underlies evolutionary patterns. We also found evidence that mimicry of inconspicuous models (pitvipers) may have evolved in association with foraging demand for crypsis. These findings demonstrate that foraging is an important factor necessary to understand the evolution, persistence, and diversity of warning coloration and mimicry of snakes, highlighting the significance of additional selective factors in solving the warning coloration paradox.

Comparative analysis of Deinagkistrodon acutus venom from Taiwan and China utilizing chromatographic, electrophoretic, and bioinformatic approaches, along with ELISA employing a monospecific antivenom

Deinagkistrodon acutus is a medically important pitviper inhabiting mainly South China and Taiwan. The hemorrhagic effects of its envenoming are compatible to its venom, which is abundant in metalloproteases (svMPs) and C-type lectin-like proteins. In this study, we investigated geographic variations in the venom of D. acutus collected from Taiwan and four Mainland Chinese provinces: Fujian, Jiangxi, Anhui, and Hunan. The variations were assessed through high-performance liquid chromatography, non-metric multidimensional scaling analysis, gel electrophoresis, and enzyme-linked immunosorbent assay (ELISA) with a monospecific antivenom (DaMAV) generated against the Taiwanese D. acutus venom, and discussed based on venom-protein sequences in databases and literature related to D. acutus venom. Additionally, the cross-reactivity of DaMAV against Crotalus horridus and Calloselasma rhodostoma venoms was investigated. We noted differential abundances of D. acutus venom metalloproteases, C-type lectin-like proteins, and phospholipase A2, along with point mutations and selective expression of serine protease isoforms. The ELISA results revealed that the venom from Taiwan was more reactive toward Taiwanese DaMAV than the four Mainland Chinese venoms, consistent with chromatographic profile differences, whereas C. horridus venom presented moderate cross-reactivity with DaMAV. The observed immunoreactivities of these venom with DaMAV can be attributed to the high prevalence of their PIII-svMPs, which are the dominant antigens, and the conservation of PIII-svMP epitopes.

Evaluation of venom diversity and antivenom quality from the venom of long-term captive vs recently wild captured Pseudocerastes persicus snake: An In vitro and In vivo study

Snakebite envenomation is a life-threatening condition and antivenoms are used as the most effective treatment. Venom obtained from snakes in long-term captivity showed some variations in comparison to the venom of the wild snakes. The objective of this study is to compare the venom of the Pseudocerastes persicus under long-term captivity and wild conditions as well as the antivenom obtained from these venoms. We have analyzed venom samples and produced trivalent antivenoms using the venom of long-term captive (LTC) or recently wild-captured (RWC) Pseudocerastes persicus, and RWC Macrovipera lebetina, and Echis carinatus. The HPLC analysis revealed that the RWC snakes' venom had three peaks that were not present in the LTC snake's venom. Further analysis using MALDI-TOF and MS/MS showed that the fraction with a retention time (RT) of 14 min contained a toxin from the Kunitz-type serine protease inhibitor (KUT) class, while the fraction with RT 21 a peptide identified within the snake venom metalloproteinase (SVMP) class. The third peak was identified as a sphingolipid. Interestingly, the in vivo preclinical tests showed no significant differences in the effectiveness of the antivenoms. which could be due to the cross-immunogenicity or cross-reactivity between different toxins in the venom. According to our results, small variations in the venom composition of a species do not lead to a decrease in the efficacy of the polyvalent antivenom.

Venom Composition of Neglected Bothropoid Snakes from the Amazon Rainforest: Ecological and Toxinological Implications

Snake venoms have evolved in several families of Caenophidae, and their toxins have been assumed to be biochemical weapons with a role as a trophic adaptation. However, it remains unclear how venom contributes to the success of venomous species for adaptation to different environments. Here we compared the venoms from Bothrocophias hyoprora, Bothrops taeniatus, Bothrops bilineatus smaragdinus, Bothrops brazili, and Bothrops atrox collected in the Amazon Rainforest, aiming to understand the ecological and toxinological consequences of venom composition. Transcriptomic and proteomic analyses indicated that the venoms presented the same toxin groups characteristic from bothropoids, but with distinct isoforms with variable qualitative and quantitative abundances, contributing to distinct enzymatic and toxic effects. Despite the particularities of each venom, commercial Bothrops antivenom recognized the venom components and neutralized the lethality of all species. No clear features could be observed between venoms from arboreal and terrestrial habitats, nor in the dispersion of the species throughout the Amazon habitats, supporting the notion that venom composition may not shape the ecological or toxinological characteristics of these snake species and that other factors influence their foraging or dispersal in different ecological niches.

Proteome analysis of Daboia russelii venom, a medically important snake from the Indian sub-continent

Daboia russelii is a category-I medically important snake throughout the Indian sub-continent contributing to majority of snakebite incidences in this part of the world. As such, extensive studies on its venom composition and search of efficient and appropriate interventions for its treatment become crucial. In this study, the proteome of Daboia russelii venom from Tanore, Rajshahi, Bangladesh was profiled using a combination of chromatographic and mass spectrometric techniques. A total of 37 different proteins belonging to 11 different snake venom protein families were detected. Proteomics analysis revealed the presence of major phospholipase A2 toxins. Daboiatoxin (both A and B subunits), the main lethal PLA2 toxin in the venom of Daboia siamensis (Myanmar viper) which is neurotoxic, myotoxic and cytotoxic was detected. Presence of Daboxin P, which is a major protein in the venom of Indian Daboia russelii with strong anticoagulant activity, was also observed. Inconsistent distribution of such lethal toxins in the venom of same species calls for more investigations of snake venoms from lesser explored regions and formulation of better alternatives to the current antivenom therapy for efficient treatment.

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.

The best of both worlds? Rattlesnake hybrid zones generate complex combinations of divergent venom phenotypes that retain high toxicity

Studying the consequences of hybridization between closely related species with divergent traits can reveal patterns of evolution that shape and maintain extreme trophic adaptations. Snake venoms are an excellent model system for examining the evolutionary and ecological patterns that underlie highly selected polymorphic traits. Here we investigate hybrid venom phenotypes that result from natural introgression between two rattlesnake species that express highly divergent venom phenotypes: Crotalus o. concolor and C. v. viridis. Though not yet documented, interbreeding between these species may lead to novel venom phenotypes with unique activities that break the typical trends of venom composition in rattlesnakes. The characteristics of these unusual phenotypes could unveil the roles of introgression in maintaining patterns of venom composition and variation, including the near ubiquitous dichotomy between neurotoxic or degradative venoms observed across rattlesnakes. We use RADseq data to infer patterns of gene flow and hybrid ancestry between these diverged lineages and link these genetic data with analyses of venom composition, biological activity, and whole animal model toxicity tests to understand the impacts of introgression on venom composition. We find that introgressed populations express admixed venom phenotypes that do not sacrifice biological activity (lethal toxicity) or overall abundance of dominant toxins compared to parental venoms. These hybridized venoms therefore do not represent a trade-off in functionality between the typical phenotypic extremes but instead represent a unique combination of characters whose expression appears limited to the hybrid zone.

Mammalian neurotoxins, Blarina paralytic peptides, cause hyperpolarization of human T-type Ca channel hCav3.2 activation

Among the rare venomous mammals, the short-tailed shrew Blarina brevicauda has been suggested to produce potent neurotoxins in its saliva to effectively capture prey. Several kallikrein-like lethal proteases have been identified, but the active substances of B. brevicauda remained unclear. Here, we report Blarina paralytic peptides (BPPs) 1 and 2 isolated from its submaxillary glands. Synthetic BPP2 showed mealworm paralysis and a hyperpolarization shift (-11 mV) of a human T-type Ca2+ channel (hCav3.2) activation. The amino acid sequences of BPPs were similar to those of synenkephalins, which are precursors of brain opioid peptide hormones that are highly conserved among mammals. However, BPPs rather resembled centipede neurotoxic peptides SLPTXs in terms of disulfide bond connectivity and stereostructure. Our results suggested that the neurotoxin BPPs were the result of convergent evolution as homologs of nontoxic endogenous peptides that are widely conserved in mammals. This finding is of great interest from the viewpoint of the chemical evolution of vertebrate venoms.

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.

The Need for Next-Generation Antivenom for Snakebite Envenomation in India

The limitations posed by currently available antivenoms have emphasized the need for alternative treatments to counteract snakebite envenomation. Even though exact epidemiological data are lacking, reports have indicated that most global snakebite deaths are reported in India. Among the many problems associated with snakebite envenomation, issues related to the availability of safer and more efficient antivenoms are of primary concern. Since India has the highest number of global snakebite deaths, efforts should be made to reduce the burden associated with snakebite envenoming. Alternative methods, including aptamers, camel antivenoms, phage display techniques for generating high-affinity antibodies and antibody fragments, small-molecule inhibitors, and natural products, are currently being investigated for their effectiveness. These alternative methods have shown promise in vitro, but their in vivo effectiveness should also be evaluated. In this review, the issues associated with Indian polyvalent antivenoms in neutralizing venom components from geographically distant species are discussed in detail. In a nutshell, this review gives an overview of the current drawbacks of using animal-derived antivenoms and several alternative strategies that are currently being widely explored.

Comparative Analysis of Alpha-1 Orthosteric-Site Binding by a Clade of Central American Pit Vipers (Genera Atropoides, Cerrophidion, Metlapilcoatlus, and Porthidium)

The distribution and relative potency of post-synaptic neurotoxic activity within Crotalinae venoms has been the subject of less investigation in comparison with Elapidae snake venoms. No previous studies have investigated post-synaptic neurotoxic activity within the Atropoides, Metlapilcoatlus, Cerrophidion, and Porthidium clade. Given the specificity of neurotoxins to relevant prey types, we aimed to uncover any activity present within this clade of snakes that may have been overlooked due to lower potency upon humans and thus not appearing as a clinical feature. Using biolayer interferometry, we assessed the relative binding of crude venoms to amphibian, lizard, bird, rodent and human α-1 nAChR orthosteric sites. We report potent alpha-1 orthosteric site binding in venoms from Atropoides picadoi, Metlapilcoatlus occiduus, M. olmec, M. mexicanus, M. nummifer. Lower levels of binding, but still notable, were evident for Cerrophidion godmani, C. tzotzilorum and C. wilsoni venoms. No activity was observed for Porthidium venoms, which is consistent with significant alpha-1 orthosteric site neurotoxicity being a trait that was amplified in the last common ancestor of Atropoides/Cerrophidion/Metlapilcoatlus subsequent to the split by Porthidium. We also observed potent taxon-selective activity, with strong selection for non-mammalian targets (amphibian, lizard, and bird). As these are poorly studied snakes, much of what is known about them is from clinical reports. The lack of affinity towards mammalian targets may explain the knowledge gap in neurotoxic activity within these species, since symptoms would not appear in bite reports. This study reports novel venom activity, which was previously unreported, indicating toxins that bind to post-synaptic receptors may be more widespread in pit vipers than previously considered. While these effects appear to not be clinically significant due to lineage-specific effects, they are of significant evolutionary novelty and of biodiscovery interest. This work sets the stage for future research directions, such as the use of in vitro and in vivo models to determine whether the alpha-1 orthosteric site binding observed within this study confers neurotoxic venom activity.

The Venom Composition of the Snake Tribe Philodryadini: 'Omic' Techniques Reveal Intergeneric Variability among South American Racers

Snakes of the Philodryadini tribe are included in the Dipsadidae family, which is a diverse group of rear-fanged snakes widespread in different ecological conditions, including habitats and diet. However, little is known about the composition and effects of their venoms despite their relevance for understanding the evolution of these snakes or even their impact on the occasional cases of human envenoming. In this study, we integrated venom gland transcriptomics, venom proteomics and functional assays to characterize the venoms from eight species of the Philodryadini tribe, which includes the genus Philodryas, Chlorosoma and Xenoxybelis. The most abundant components identified in the venoms were snake venom metalloproteinases (SVMPs), cysteine-rich secretory proteins (CRISPs), C-type lectins (CTLs), snake endogenous matrix metalloproteinases type 9 (seMMP-9) and snake venom serinoproteinases (SVSPs). These protein families showed a variable expression profile in each genus. SVMPs were the most abundant components in Philodryas, while seMMP-9 and CRISPs were the most expressed in Chlorosoma and Xenoxybelis, respectively. Lineage-specific differences in venom composition were also observed among Philodryas species, whereas P. olfersii presented the highest amount of SVSPs and P. agassizii was the only species to express significant amounts of 3FTx. The variability observed in venom composition was confirmed by the venom functional assays. Philodryas species presented the highest SVMP activity, whereas Chlorosoma species showed higher levels of gelatin activity, which may correlate to the seMMP-9 enzymes. The variability observed in the composition of these venoms may be related to the tribe phylogeny and influenced by their diets. In the presented study, we expanded the set of venomics studies of the Philodryadini tribe, which paves new roads for further studies on the evolution and ecology of Dipsadidae snakes.

Snakes on a plain: biotic and abiotic factors determine venom compositional variation in a wide-ranging generalist rattlesnake

BACKGROUND

Snake venoms are trophic adaptations that represent an ideal model to examine the evolutionary factors that shape polymorphic traits under strong natural selection. Venom compositional variation is substantial within and among venomous snake species. However, the forces shaping this phenotypic complexity, as well as the potential integrated roles of biotic and abiotic factors, have received little attention. Here, we investigate geographic variation in venom composition in a wide-ranging rattlesnake (Crotalus viridis viridis) and contextualize this variation by investigating dietary, phylogenetic, and environmental variables that covary with venom.

RESULTS

Using shotgun proteomics, venom biochemical profiling, and lethality assays, we identify 2 distinct divergent phenotypes that characterize major axes of venom variation in this species: a myotoxin-rich phenotype and a snake venom metalloprotease (SVMP)-rich phenotype. We find that dietary availability and temperature-related abiotic factors are correlated with geographic trends in venom composition.

CONCLUSIONS

Our findings highlight the potential for snake venoms to vary extensively within species, for this variation to be driven by biotic and abiotic factors, and for the importance of integrating biotic and abiotic variation for understanding complex trait evolution. Links between venom variation and variation in biotic and abiotic factors indicate that venom variation likely results from substantial geographic variation in selection regimes that determine the efficacy of venom phenotypes across populations and snake species. Our results highlight the cascading influence of abiotic factors on biotic factors that ultimately shape venom phenotype, providing evidence for a central role of local selection as a key driver of venom variation.

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.

Phylogeny-Related Variations in Venomics: A Test in a Subset of Habu Snakes (Protobothrops)

We conducted a comparative analysis to unveil the divergence among venoms from a subset of Old World habu snakes (Protobothrops) in terms of venomic profiles and toxicological and enzymatic activities. A total of 14 protein families were identified in the venoms from these habu snakes, and 11 of them were shared among these venoms. The venoms of five adult habu snakes were overwhelmingly dominated by SVMP (32.56 ± 13.94%), PLA₂ (22.93 ± 9.26%), and SVSP (16.27 ± 4.79%), with a total abundance of over 65%, while the subadult P. mangshanensis had an extremely low abundance of PLA₂ (1.23%) but a high abundance of CTL (51.47%), followed by SVMP (22.06%) and SVSP (10.90%). Apparent interspecific variations in lethality and enzymatic activities were also explored in habu snake venoms, but no variations in myotoxicity were found. Except for SVSP, the resemblance of the relatives within Protobothrops in other venom traits was estimated to deviate from Brownian motion evolution based on phylogenetic signals. A comparative analysis further validated that the degree of covariation between phylogeny and venom variation is evolutionarily labile and varies among clades of closely related snakes. Our findings indicate a high level of interspecific variation in the venom proteomes of habu snakes, both in the presence or absence and the relative abundance of venom protein families, and that these venoms might have evolved under a combination of adaptive and neutral mechanisms.

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

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

Toxic and antigenic characterization of Peruvian Micrurus surinamensis coral snake venom

Micrurus surinamensis is a semi-aquatic coral snake found in primary forest region and can cause relevant human accidents. In this work we investigated the toxic and antigenic activities of the Peruvian Micrurus surinamensis venom (MsV). We found that MsV show hyaluronidase activity but lack LAAO and PLA₂ enzymatic activities. Interestingly, MsV induce edematogenic responses but cannot cause nociceptive effects. Furthermore, MsV can reduce in vitro cell viability in MGSO-3 cell line derived from human breast cancer tissue. To evaluate its antigenic potential, rabbits were immunized with MsV, which proved to be immunogenic. ELISA, immunobloting and in vivo neutralization assays demonstrated that the specific rabbit anti-MsV antivenom is more efficient than the therapeutic Brazilian antivenom in recognizing and neutralizing the lethal activity of MsV. MsV differs in protein profile and biological activities from M. frontalis venom (MfV), used as control, which impairs its recognition and neutralization by Brazilian therapeutic anti-elapidic antivenom. We performed a SPOT immunoassay for the identification of B-cell linear epitopes in the main toxins described for MsV targeted by the elicited neutralizing antibodies previously produced. A membrane containing 15-mer peptides representing the sequences of five 3TFxs and five PLA₂s was produced and probed with anti- MsV antibodies. Results revealed important regions in 3FTx toxins for venom neutralization. Identifying the main MsV components and its biological activities can be helpful in guiding the production of antivenoms and in the optimization of treatment for coral snake envenomation in Brazil.

Molecular phylogeny reveals distinct evolutionary lineages of the banded krait, Bungarus fasciatus (Squamata, Elapidae) in Asia

The banded krait, Bungarus fasciatus is a widespread elapid snake, likely to comprise several distinct species in different geographic regions of Asia. Therefore, based on molecular phylogenetics and comparative morphology data, we present an overview of the systematic composition of the species to delimit potential biogeographic boundaries. Our phylogenetic analyses, based on four mitochondrial genes, reveal the existence of at least three evolutionary lineages within B. fasciatus, corresponding to Indo-Myanmar, Sundaic and eastern Asian lineages. We are convinced that there are at least three taxonomic entities within the nomen B. fasciatus and restrict the distribution of B. fasciatus sensu stricto to the Indo-Myanmar region. We also provide additional natural history data of the taxon from eastern India. Finally, we advocate further studies to establish the degree of reproductive isolation among these diverging evolutionary lineages and to reassess the systematic status of this species complex especially the Sundaic and eastern Asian lineages.

Killer Knots: Molecular Evolution of Inhibitor Cystine Knot Toxins in Wandering Spiders (Araneae: Ctenidae)

Venom expressed by the nearly 50,000 species of spiders on Earth largely remains an untapped reservoir of a diverse array of biomolecules with potential for pharmacological and agricultural applications. A large fraction of the noxious components of spider venoms are a functionally diverse family of structurally related polypeptides with an inhibitor cystine knot (ICK) motif. The cysteine-rich nature of these toxins makes structural elucidation difficult, and most studies have focused on venom components from the small handful of medically relevant spider species such as the highly aggressive Brazilian wandering spider Phoneutria nigriventer. To alleviate difficulties associated with the study of ICK toxins in spiders, we devised a comprehensive approach to explore the evolutionary patterns that have shaped ICK functional diversification using venom gland transcriptomes and proteomes from phylogenetically distinct lineages of wandering spiders and their close relatives. We identified 626 unique ICK toxins belonging to seven topological elaborations. Phylogenetic tests of episodic diversification revealed distinct regions between cysteine residues that demonstrated differential evidence of positive or negative selection, which may have structural implications towards the specificity and efficacy of these toxins. Increased taxon sampling and whole genome sequencing will provide invaluable insights to further understand the evolutionary processes that have given rise to this diverse class of toxins.

Highly Evolvable: Investigating Interspecific and Intraspecific Venom Variation in Taipans (Oxyuranus spp.) and Brown Snakes (Pseudonaja spp.)

Snake venoms are complex mixtures of toxins that differ on interspecific (between species) and intraspecific (within species) levels. Whether venom variation within a group of closely related species is explained by the presence, absence and/or relative abundances of venom toxins remains largely unknown. Taipans (Oxyuranus spp.) and brown snakes (Pseudonaja spp.) represent medically relevant species of snakes across the Australasian region and provide an excellent model clade for studying interspecific and intraspecific venom variation. Using liquid chromatography with ultraviolet and mass spectrometry detection, we analyzed a total of 31 venoms covering all species of this monophyletic clade, including widespread localities. Our results reveal major interspecific and intraspecific venom variation in Oxyuranus and Pseudonaja species, partially corresponding with their geographical regions and phylogenetic relationships. This extensive venom variability is generated by a combination of the absence/presence and differential abundance of venom toxins. Our study highlights that venom systems can be highly dynamical on the interspecific and intraspecific levels and underscores that the rapid toxin evolvability potentially causes major impacts on neglected tropical snakebites.

Montane Rattlesnakes in México: Venoms of Crotalus tancitarensis and Related Species within the Crotalus intermedius Group

The Crotalus intermedius group is a clade of rattlesnakes consisting of several species adapted to a high elevation habitat, primarily in México. Crotalus tancitarensis was previously classified as C. intermedius, until individuals occurring on Cerro Tancítaro in Michoacán, México, were reevaluated and classified as a new species (C. tancitarensis) based on scale pattern and geographic location. This study aimed to characterize the venom of C. tancitarensis and compare the venom profile to those of other species within the Crotalus intermedius group using gel electrophoresis, biochemical assays, reverse-phase high performance liquid chromatography, mass spectrometry, and lethal toxicity (LD50) assays. Results show that the venom profiles of species within the Crotalus intermedius group are similar, but with distinct differences in phospholipase A2 (PLA2), metalloproteinase PI (SVMP PI), and kallikrein-like serine proteinase (SVSP) activity and relative abundance. Proteomic analysis indicated that the highland forms produce venoms with 50-60 protein isoforms and a composition typical of type I rattlesnake venoms (abundant SVMPs, lack of presynaptic PLA2-based neurotoxins), as well as a diversity of typical Crotalus venom components such as serine proteinases, PLA2s, C-type lectins, and less abundant toxins (LAAOs, CRiSPs, etc.). The overall venom profile of C. tancitarensis appears most similar to C. transversus, which is consistent with a previous mitochondrial DNA analysis of the Crotalus intermedius group. These rattlesnakes of the Mexican highlands represent a radiation of high elevation specialists, and in spite of divergence of species in these Sky Island habitats, venom composition of species analyzed here has remained relatively conserved. The majority of protein family isoforms are conserved in all members of the clade, and as seen in other more broadly distributed rattlesnake species, differences in their venoms are largely due to relative concentrations of specific components.

Spatial Venomics─Cobra Venom System Reveals Spatial Differentiation of Snake Toxins by Mass Spectrometry Imaging

Among venomous animals, toxic secretions have evolved as biochemical weapons associated with various highly specialized delivery systems on many occasions. Despite extensive research, there is still limited knowledge of the functional biology of most animal toxins, including their venom production and storage, as well as the morphological structures within sophisticated venom producing tissues that might underpin venom modulation. Here, we report on the spatial exploration of a snake venom gland system by matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI), in combination with standard proteotranscriptomic approaches, to enable in situ toxin mapping in spatial intensity maps across a venom gland sourced from the Egyptian cobra (Naja haje). MALDI-MSI toxin visualization on the elapid venom gland reveals a high spatial heterogeneity of different toxin classes at the proteoform level, which may be the result of physiological constraints on venom production and/or storage that reflects the potential for venom modulation under diverse stimuli.

Equipped for Sexual Stings? Male-Specific Venom Peptides in Euscorpius italicus

In the animal kingdom, intraspecific variation occurs, for example, between populations, different life stages, and sexes. For venomous animals, this can involve differences in their venom composition. In cases where venom is utilized in the context of mating, the differences in composition might be driven by sexual selection. In this regard, the genus Euscorpius is a promising group for further research, as some of these scorpions exhibit a distinct sexual dimorphism and are known to perform a sexual sting during mating. However, the venom composition of this genus remains largely unexplored. Here, we demonstrate that Euscorpius italicus exhibits a male-specific venom composition, and we identify a large fraction of the substances involved. The sex specificity of venom peptides was first determined by analyzing the presence/absence patterns of ion signals in MALDI-TOF mass spectra of venom samples from both sexes and juveniles. Subsequently, a proteo-transcriptomic analysis provided sequence information on the relevant venom peptides and their corresponding precursors. As a result, we show that several potential toxin precursors are down-regulated in male venom glands, possibly to reduce toxic effects caused to females during the sexual sting. We have identified the precursor of one of the most prominent male-specific venom peptides, which may be an ideal candidate for activity tests in future studies. In addition to the description of male-specific features in the venom of E. italicus, this study also includes a general survey of venom precursors in this species.

Venom Variation of Neonate and Adult Chinese Cobras in Captivity Concerning Their Foraging Strategies

The venom and transcriptome profile of the captive Chinese cobra (Naja atra) is not characterized until now. Here, LC-MS/MS and illumine technology were used to unveil the venom and trascriptome of neonates and adults N. atra specimens. In captive Chinese cobra, 98 co-existing transcripts for venom-related proteins was contained. A total of 127 proteins belong to 21 protein families were found in the profile of venom. The main components of snake venom were three finger toxins (3-FTx), snake venom metalloproteinase (SVMP), cysteine-rich secretory protein (CRISP), cobra venom factor (CVF), and phosphodiesterase (PDE). During the ontogenesis of captive Chinese cobra, the rearrangement of snake venom composition occurred and with obscure gender difference. CVF, 3-FTx, PDE, phospholipase A2 (PLA2) in adults were more abundant than neonates, while SVMP and CRISP in the neonates was richer than the adults. Ontogenetic changes in the proteome of Chinese cobra venom reveals different strategies for handling prey. The levels of different types of toxin families were dramatically altered in the wild and captive specimens. Therefore, we speculate that the captive process could reshape the snake venom composition vigorously. The clear comprehension of the composition of Chinese cobra venom facilitates the understanding of the mechanism of snakebite intoxication and guides the preparation and administration of traditional antivenom and next-generation drugs for snakebite.

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.

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.

Intraspecific Differences in the Venom of Crotalus durissus cumanensis from Colombia

Biochemical and biological differences in the venom of Crotalus durissus cumanensis from three ecoregions of Colombia were evaluated. Rattlesnakes were collected from the geographic areas of Magdalena Medio (MM), Caribe (CA) and Orinoquía (OR). All three regionally distributed venoms contain proteases, PLA₂s and the basic subunit of crotoxin. However, only crotamine was detected in the CA venom. The highest lethality, coagulant, phospholipase A₂ and hyaluronidase activities were found in the MM venom. Also, some differences, observed by western blot and immunoaffinity, were found in all three venoms when using commercial antivenoms. Furthermore, all three eco-regional venoms showed intraspecific variability, considering the differences in the abundance and intensity of their components, in addition to the activity and response to commercial antivenoms.

Venom composition of Trimeresurus albolabris, T. insularis, T. puniceus and T. purpureomaculatus from Indonesia

Background:

Several studies have been published on the characterization of Trimeresurus venoms. However, there is still limited information concerning the venom composition of Trimeresurus species distributed throughout Indonesia, which contributes to significant snakebite envenomation cases. The present study describes a comparative on the composition of T. albolabris, T. insularis, T. puniceus, and T. purpureomaculatus venoms originated from Indonesia.

Methods:

Protein content in the venom of four Trimeresurus species was determined using Bradford assay, and the venom proteome was elucidated using one-dimension SDS PAGE nano-ESI- LCMS/MS shotgun proteomics.

Results:

The venom of T. albolabris contained the highest protein content of 11.1 mg/mL, followed by T. puniceus, T. insularis and T. purpureomaculatus venom with 10.7 mg/mL, 8.9 mg/mL and 5.54 mg/mL protein, respectively. In total, our venomic analysis identified 65 proteins belonging to 16 protein families in T. purpureomaculatus; 64 proteins belonging to 18 protein families in T. albolabris; 58 different proteins belonging to 14 protein families in T. puniceus; and 48 different proteins belonging to 14 protein familiesin T. insularis. Four major proteins identified in all venoms belonged to snake venom metalloproteinase, C-type lectin, snake venom serine protease, and phospholipase A2. There were 11 common proteins in all venoms, and T. puniceus venom has the highest number of unique proteins compared to the other three venoms. Cluster analysis of the proteins and venoms showed that T. puniceus venom has the most distinct venom composition.

Conclusions:

Overall, the results highlighted venom compositional variation of four Trimeresurus spp. from Indonesia. The venoms appear to be highly similar, comprising at least four protein families that correlate with venom’s toxin properties and function. This study adds more information on venom variability among Trimeresurus species within the close geographic origin and may contribute to the development of optimum heterologous antivenom.

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

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

Pharmacological re-assessment of traditional medicinal plants-derived inhibitors as antidotes against snakebite envenoming: A critical review

ETHNOPHARMACOLOGICAL RELEVANCE

Traditional healers have used medicinal plants to treat snakebite envenomation worldwide; however, mostly without scientific validation. There have been many studies on the therapeutic potential of the natural products against snake envenomation.

AIM OF THE STUDY

This review has highlighted snake venom inhibitory activity of bioactive compounds and peptides from plants that have found a traditional use in treating snakebite envenomation. We have systematically reviewed the scenario of different phases of natural snake venom inhibitors characterization covering a period from 1994 until the present and critically analysed the lacuna of the studies if any, and further scope for their translation from bench to bedside.

MATERIALS AND METHODS

The medicinal plant-derived compounds used against snakebite therapy were reviewed from the available literature in public databases (Scopus, MEDLINE) from 1994 till 2020. The search words used were 'natural inhibitors against snakebite,' 'natural products as therapeutics against snakebite,' 'natural products as antidote against snake envenomation,' ' snake venom toxin natural inhibitors,' 'snake venom herbal inhibitors'. However, the scope of this review does not include computational (in silico) predictions without any wet laboratory validation and snake venom inhibitory activity of the crude plant extracts. In addition, we have also predicted the ADMET properties of the identified snake venom inhibitors to highlight their valuable pharmacokinetics for future clinical studies.

RESULTS

The therapeutic application of plant-derived natural inhibitors to treat snakebite envenomation as an auxiliary to antivenom therapy has been gaining significant momentum. Pharmacological reassessment of the natural compounds derived from traditional medicinal plants has demonstrated inhibition of the principal toxic enzymes of snake venoms at various extents to curb the lethal and/or deleterious effects of venomous snakebite. Nevertheless, such molecules are yet to be commercialized for clinical application in the treatment of snakebite. There are many obstacles in the marketability of the plant-derived natural products as snake envenomation antidote and strategies must be explored for the translation of these compounds from drug candidates to their clinical application.

CONCLUSION

In order to minimize the adverse implications of snake envenomation, strategies must be developed for the smooth transition of these plant-derived small molecule inhibitors from bench to bedside. In this article we have presented an inclusive review and have critically analysed natural products for their therapeutic potential against snake envenomation, and have proposed a road map for use of natural products as antidote against snakebite.

A comparison of two different analytical workflows to determine the venom proteome composition of Naja kaouthia from North-East India and immunological profiling of venom against commercial antivenoms

The Indian monocled cobra (Naja kaouthia) is one of the most prevalent venomous snakes in northeast India (NEI) and is the cause of many fatalities. The composition of NEI N. kaouthia venom (NkV) was deciphered using two different proteomic approaches: (i) 1D SDS-PAGE coupled to label-free quantification of protein bands using stringent identification criteria and (ii) reversed-phase high-performance liquid chromatography (RP-HPLC) followed by quantification based on area under the RP-HPLC peaks. The proteomic data from both strategies were compared. Proteomic analyses from both workflows identified 32 proteins (toxins) distributed over 10-14 snake venom protein families in NEI NkV. The relative abundances of the venom proteins determined from the analytical workflows coincided with the densitometry band intensities of the NEI NkV. Phospholipase A₂ (13.1-16.0%) and three-finger toxins (58.5-64.2%) represented the most abundant enzymatic and non-enzymatic proteins in NEI NkV, respectively. Immuno-cross-reactivity studies by enzyme-linked immunoassay and immunoblot analyses pointed to the poor efficacy of commercial PAVs in recognizing the low molecular mass (<15 kDa) toxins of NEI NkV. Spectrofluorometric titration determined the presence of NEI NkV-specific antibodies in commercial PAV, at a level that was higher than that previously reported for eastern India NkV-specific antibodies in commercial antivenom.

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.

The rise of genomics in snake venom research: recent advances and future perspectives

Snake venoms represent a danger to human health, but also a gold mine of bioactive proteins that can be harnessed for drug discovery purposes. The evolution of snakes and their venom has been studied for decades, particularly via traditional morphological and basic genetic methods alongside venom proteomics. However, while the field of genomics has matured rapidly over the past 2 decades, owing to the development of next-generation sequencing technologies, snake genomics remains in its infancy. Here, we provide an overview of the state of the art in snake genomics and discuss its potential implications for studying venom evolution and toxinology. On the basis of current knowledge, gene duplication and positive selection are key mechanisms in the neofunctionalization of snake venom proteins. This makes snake venoms important evolutionary drivers that explain the remarkable venom diversification and adaptive variation observed in these reptiles. Gene duplication and neofunctionalization have also generated a large number of repeat sequences in snake genomes that pose a significant challenge to DNA sequencing, resulting in the need for substantial computational resources and longer sequencing read length for high-quality genome assembly. Fortunately, owing to constantly improving sequencing technologies and computational tools, we are now able to explore the molecular mechanisms of snake venom evolution in unprecedented detail. Such novel insights have the potential to affect the design and development of antivenoms and possibly other drugs, as well as provide new fundamental knowledge on snake biology and evolution.

Bitiscetin-3, a Novel C-Type Lectin-like Protein Cloned from the Venom Gland of the Viper Bitis arietans, Induces Platelet Agglutination and Inhibits Binding of Von Willebrand Factor to Collagen

Bitiscetin-1 (aka bitiscetin) and bitiscetin-2 are C-type lectin-like proteins purified from the venom of Bitis arietans (puff adder). They bind to von Willebrand factor (VWF) and—at least bitiscetin-1—induce platelet agglutination via enhancement of VWF binding to platelet glycoprotein Ib (GPIb). Bitiscetin-1 and -2 bind the VWF A1 and A3 domains, respectively. The A3 domain includes the major site of VWF for binding collagen, explaining why bitiscetin-2 blocks VWF-to-collagen binding. In the present study, sequences for a novel bitiscetin protein—bitiscetin-3—were identified in cDNA constructed from the B. arietans venom gland. The deduced amino acid sequences of bitiscetin-3 subunits α and β share 79 and 80% identity with those of bitiscetin-1, respectively. Expression vectors for bitiscetin-3α and -3β were co-transfected to 293T cells, producing the heterodimer protein recombinant bitiscetin-3 (rBit-3). Functionally, purified rBit-3 (1) induced platelet agglutination involving VWF and GPIb, (2) did not compete with bitiscetin-1 for binding to VWF, (3) blocked VWF-to-collagen binding, and (4) lost its platelet agglutination inducing ability in the presence of an anti-VWF monoclonal antibody that blocked VWF-to-collagen binding. These combined results suggest that bitiscetin-3 binds to the A3 domain, as does bitiscetin-2. Except for a small N-terminal fragment of a single subunit—which differs from that of both bitiscetin-3 subunits—the sequences of bitiscetin-2 have never been determined. Therefore, by identifying and analyzing bitiscetin-3, the present study is the first to present the full-length α- and β-subunit sequences and recombinant expression of a bitiscetin-family toxin that blocks the binding of VWF to collagen.

Immunorecognition and Neutralization of Crotalus durissus cumanensis Venom by a Commercial Antivenom Produced in Colombia

In Colombia, on average 2.9% of the nearly 5600 snakebite events that occur annually involve the rattlesnake Crotalus durissus cumanensis. The envenomation by this snake is mainly characterized by neurotoxicity and the main toxin is crotoxin (~64.7% of the total venom). The Instituto Nacional de Salud (INS) produces a polyvalent antivenom aimed at the treatment of bothropic, crotalid, and lachesic envenomations; nonetheless, its immune reactivity profile and neutralizing capacity over biological activities of the C. d. cumanensis venom has been poorly evaluated. In this sense, the study aims: (1) to describe an in-depth exploration of its immunoreactivity through second-generation antivenomics and HPLC fraction-specific ELISA immunoprofiles; and (2) to evaluate the neutralization pattern of the rattlesnake venom in vitro and in vivo biological activities. The results obtained showed a variable recognition of crotoxin subunits, in addition to a molecular mass-dependent immunoreactivity pattern in which the disintegrins were not recognized, and snake venom metalloproteinases and L-amino acid oxidases were the most recognized. Additionally, a high neutralization of proteolytic and coagulant activities was observed, but not over the PLA₂ activity. Further, the median effective dose against C. d. cumanensis venom lethality was 962 μL of antivenom per mg of venom. In conclusion, (1) the antivenom recognition over the crotoxin and the disintegrins of the C. d. cumanensis should be improved, thus aiming upcoming efforts for the exploration of new techniques and approaches in antivenom production in Colombia, and (2) the neutralization activity of the antivenom seems to follow the molecular mass-dependent recognition pattern, although other explanations should be explored.