Evolution of an Arsenal

An ancient, conserved gene regulatory network led to the rise of oral venom systems

Although oral venom systems are ecologically important characters, how they originated is still unclear. In this study, we show that oral venom systems likely originated from a gene regulatory network conserved across amniotes. This network, which we term the “metavenom network,” comprises over 3,000 housekeeping genes coexpressed with venom and play a role in protein folding and modification. Comparative transcriptomics revealed that the network is conserved between venom glands of snakes and salivary glands of mammals. This suggests that while these tissues have evolved different functions, they share a common regulatory core, that persisted since their common ancestor. We propose several evolutionary mechanisms that can utilize this common regulatory core to give rise to venomous animals from their nonvenomous ancestors.

Oral venom systems evolved multiple times in numerous vertebrates enabling the exploitation of unique predatory niches. Yet how and when they evolved remains poorly understood. Up to now, most research on venom evolution has focused strictly on the toxins. However, using toxins present in modern day animals to trace the origin of the venom system is difficult, since they tend to evolve rapidly, show complex patterns of expression, and were incorporated into the venom arsenal relatively recently. Here we focus on gene regulatory networks associated with the production of toxins in snakes, rather than the toxins themselves. We found that overall venom gland gene expression was surprisingly well conserved when compared to salivary glands of other amniotes. We characterized the “metavenom network,” a network of ∼3,000 nonsecreted housekeeping genes that are strongly coexpressed with the toxins, and are primarily involved in protein folding and modification. Conserved across amniotes, this network was coopted for venom evolution by exaptation of existing members and the recruitment of new toxin genes. For instance, starting from this common molecular foundation, Heloderma lizards, shrews, and solenodon, evolved venoms in parallel by overexpression of kallikreins, which were common in ancestral saliva and induce vasodilation when injected, causing circulatory shock. Derived venoms, such as those of snakes, incorporated novel toxins, though still rely on hypotension for prey immobilization. These similarities suggest repeated cooption of shared molecular machinery for the evolution of oral venom in mammals and reptiles, blurring the line between truly venomous animals and their ancestors.

Proteomic analysis of three medically important Nigerian Naja (Naja haje, Naja katiensis and Naja nigricollis) snake venoms

Proteomics technologies enable a comprehensive study of complex proteins and their functions. The venom proteomes of three medically important Nigerian Elapidae snakes Naja haje, Naja katiensis and Naja nigricollis was studied using HILIC coupled with LC-MS/MS analysis. Results revealed a total of 57, 55, and 46 proteins in the venoms of N. haje, N. katiensis, and N. nigricollis, respectively, with molecular mass ranging between 5 and 185  kDa. These snakes have 38 common proteins in addition to 3 uncommon proteins: actiflagelin, cathelicidin, and cystatin identified in their venoms. The identified proteins belonged to 14 protein families in N. haje and N. katiensis, and 12 protein families in N. nigricollis. Of the total venom proteins, 3FTx was the most abundant protein family, constituting 52% in N. haje and N. katiensis, and 41% in N. nigricollis, followed by PLA₂, constituting 37% in N. nigricollis, 26% in N. haje, and 24% in N. katiensis. Other protein families, including LAAO, CRISPs, VEGF, PLB, CVF, SVMP, SVH, AMP, PI, Globin, Actin, and C-type lectins, were also detected, although, at very low abundances. Quantification of the relative abundance of each protein revealed that alpha and beta fibrinogenase and PLA_2, which constituted 18-26% of the total proteome, were the most abundant. The 3 uncommon proteins have no known function in snake venom. However, actiflagelin activates sperm motility; cystatin inhibits angiogenesis, while cathelicidin exerts antimicrobial effects. The three Nigerian Naja genus proteomes displayed 70% similarity in composition, which suggests the possibility of formulating antivenom that may cross-neutralise the venoms of cobra species found in Nigeria. These data provide insights into clinically relevant peptides/proteins present in the venoms of these snakes. Data are available via ProteomeXchange with identifier PXD024627.

Comparative venomics and preclinical efficacy evaluation of a monospecific Hemachatus antivenom towards sub-Saharan Africa cobra venoms

Cobras are the most medically important elapid snakes in Africa. The African genera Naja and Hemachatus include snakes with neurotoxic and cytotoxic venoms, with shared biochemical, toxinological and antigenic characteristics. We have studied the antigenic cross-reactivity of four sub-Saharan Africa cobra venoms against an experimental monospecific Hemachatus haemachatus antivenom through comparative proteomics, preclinical assessment of neutralization, and third generation antivenomics. The venoms of H. haemachatus, N. annulifera, N. mossambica and N. nigricollis share an overall qualitative family toxin composition but depart in their proportions of three-finger toxin (3FTxs) classes, phospholipases A₂ (PLA₂s), snake venom metalloproteinases (SVMPs), and cysteine-rich secretory proteins (CRISPs). A monospecific anti-Hemachatus antivenom produced by Costa Rican Instituto Clodomiro Picado neutralized the lethal activity of the homologous and heterologous neuro/cytotoxic (H. haemachatus) and cyto/cardiotoxic (N. mossambica and N. nigricollis) venoms of the three spitting cobras sampled, while it was ineffective against the lethal and toxic activities of the neurotoxic venom of the non-spitting snouted cobra N. annulifera. The ability of the anti-Hemachatus-ICP antivenom to neutralize toxic (dermonecrotic and anticoagulant) and enzymatic (PLA₂) activities of spitting cobra venoms suggested a closer kinship of H. haemachatus and Naja subgenus Afrocobra spitting cobras than to Naja subgenus Uraeus neurotoxic taxa. These results were confirmed by third generation antivenomics. BIOLOGICAL SIGNIFICANCE: African Naja species represent the most widespread medically important elapid snakes across Africa. To gain deeper insight into the spectrum of medically relevant toxins, we compared the proteome of three spitting cobras (Hemachatus haemachatus, Naja mossambica and N. nigricollis) and one non-spitting cobra (N. annulifera). Three finger toxins and phospholipases A₂ are the two major protein families among the venoms analyzed. The development of antivenoms of broad species coverage is an urgent need in sub-Saharan Africa. An equine antivenom raised against H. haemachatus venom showed cross-reactivity with the venoms of H. haemachatus, N. mossambica and N. nigricollis, while having poor recognition of the venom of N. annulifera. This immunological information provides clues for the design of optimum venom mixtures for the preparation of broad spectrum antivenoms.

Coagulopathy after snake envenomation in South Korea

OBJECTIVE

This retrospective study investigated the nature and severity of venom-induced consumption coagulopathy (VICC) and determined the clotting factors involved in VICC in patients after envenomation by South Korea's snakes. Additionally, we studied the effectiveness of antivenom for the treatment of VICC after envenomation.

METHODS

Included patients were divided into three groups according to the severity of VICC (no VICC, partial VICC, and complete VICC). Data, including changes in coagulation parameters during hospitalization and clotting factors at presentation, were collected and analyzed.

RESULTS

One hundred nineteen patients who presented at our emergency department within 3 h after snake envenomation were included. VICC developed in 34 patients (27 patients with partial VICC and 7 patients with complete VICC). Two of 34 patients with VICC required blood transfusions. Five patients with complete VICC had an undetectable fibrinogen concentration at presentation. Three patients with complete VICC had an unmeasurable INR and aPTT within 24 h. The median times of the most extreme values were 10 h for INR, 12 h for aPTT, and 16 h for fibrinogen after presentation in the VICC group. The D-dimer concentration peaked at a median of 63.5 h after presentation. The activities of factors II and X were significantly reduced in the complete VICC group (factor II: 88 (84-99.3)% in the non-VICC group vs. 69 (49.5-83.5)% in the complete VICC group; factor X:94 (83-102) in the non-VICC group vs. 70 (66.5-79.8)% in the complete VICC group), while there was no difference in factor V activity at presentation. The time from bite to first antivenom administration did not correlate with the time course and most extreme concentrations for fibrinogen and D-dimer within the VICC groups.

DISCUSSION AND CONCLUSION

VICC occurs in approximately one-quarter of snakebite patients in South Korea; however, VICC itself does not appear to lead to clinical deterioration. Fibrinogen is an early diagnostic maker for complete VICC. Clotting factors II and X are involved in VICC. Future investigations should explore the mechanism of VICC from Korean snakebites and the effect of antivenom on VICC.

Electric Blue: Molecular Evolution of Three-Finger Toxins in the Long-Glanded Coral Snake Species Calliophis bivirgatus

The genus Calliophis is the most basal branch of the family Elapidae and several species in it have developed highly elongated venom glands. Recent research has shown that C. bivirgatus has evolved a seemingly unique toxin (calliotoxin) that produces spastic paralysis in their prey by acting on the voltage-gated sodium (NaV) channels. We assembled a transcriptome from C. bivirgatus to investigate the molecular characteristics of these toxins and the venom as a whole. We find strong confirmation that this genus produces the classic elapid eight-cysteine three-finger toxins, that δδ-elapitoxins (toxins that resemble calliotoxin) are responsible for a substantial portion of the venom composition, and that these toxins form a distinct clade within a larger, more diverse clade of C. bivirgatus three-finger toxins. This broader clade of C. bivirgatus toxins also contains the previously named maticotoxins and is somewhat closely related to cytotoxins from other elapids. However, the toxins from this clade that have been characterized are not themselves cytotoxic. No other toxins show clear relationships to toxins of known function from other species.

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

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

A Transcriptomic Approach to the Recruitment of Venom Proteins in a Marine Annelid

The growing number of known venomous marine invertebrates indicates that chemical warfare plays an important role in adapting to diversified ecological niches, even though it remains unclear how toxins fit into the evolutionary history of these animals. Our case study, the Polychaeta Eulalia sp., is an intertidal predator that secretes toxins. Whole-transcriptome sequencing revealed proteinaceous toxins secreted by cells in the proboscis and delivered by mucus. Toxins and accompanying enzymes promote permeabilization, coagulation impairment and the blocking of the neuromuscular activity of prey upon which the worm feeds by sucking pieces of live flesh. The main neurotoxins ("phyllotoxins") were found to be cysteine-rich proteins, a class of substances ubiquitous among venomous animals. Some toxins were phylogenetically related to Polychaeta, Mollusca or more ancient groups, such as Cnidaria. Some toxins may have evolved from non-toxin homologs that were recruited without the reduction in molecular mass and increased specificity of other invertebrate toxins. By analyzing the phylogeny of toxin mixtures, we show that Polychaeta is uniquely positioned in the evolution of animal venoms. Indeed, the phylogenetic models of mixed or individual toxins do not follow the expected eumetazoan tree-of-life and highlight that the recruitment of gene products for a role in venom systems is complex.

Illegal online pet trade in venomous snakes and the occurrence of snakebites in Brazil

Greater access to the Internet has boosted the online pet trade and especially enabled the trade in niche-targeted groups, such as venomous species. Despite their fearsome reputation, an array of venomous animals, such as snakes, spiders and scorpions arise interest among pet hobbyists, which exposes owners, sellers, and others involved in their transport and maintenance to potentially serious accidents by envenomation. To assess the potential risk to human health posed by the trade and ownership of venomous pet snakes, we examined social media posts trading or portraying native and exotic venomous species as pets (Facebook™ and YouTube™) and official seizures in Brazil between 2015 and 2020. In addition, we recorded all venomous snakebite events occurring in the country during the same period from the Brazilian Ministry of Health database. We compiled 114,931 venomous snakebite events resulting in over 600 human deaths. Bothrops spp. and Crotalus spp. were responsible for 86% and 10% of these events. We recorded 241 individual venomous snakes kept as pets of at least 16 species (seven native and nine non-native). The taxa with higher numbers of snakebites were also those most often kept as pets. Venomous pet snakes were sold at low prices, US$ 71.70 ± 18.44 on average. Our findings highlight a dangerous market, given the similarity of taxa raised as pets and those causing snakebites, combined with low availability of anti-venom for exotic species in Brazil. Additionally, several of the venomous pet snakes recorded are invasive species in many countries. Trade in venomous snakes is prohibited in Brazil, but the law is not well-enforced. To curb this illegal market and discourage consumer demand, we suggest that tougher penalties for sellers and owners should be considered, along with the development of awareness campaigns on the consequences of the snakebite injuries and the lack of antivenoms.

Widespread and Differential Neurotoxicity in Venoms from the Bitis Genus of Viperid Snakes

Research into the neurotoxic activity of venoms from species within the snake family Viperidae is relatively neglected compared with snakes in the Elapidae family. Previous studies into venoms from the Bitis genus of vipers have identified the presence of presynaptic phospholipase A₂ neurotoxins in B. atropos and B. caudalis, as well as a postsynaptic phospholipase A₂ in B. arietans. Yet, no studies have investigated how widespread neurotoxicity is across the Bitis genus or if they exhibit prey selectivity of their neurotoxins. Utilising a biolayer interferometry assay, we were able to assess the binding of crude venom from 14 species of Bitis to the neuromuscular α-1 nAChR orthosteric site across a wide range of vertebrate taxa mimotopes. Postsynaptic binding was seen for venoms from B. arietans, B. armata, B. atropos, B. caudalis, B. cornuta, B. peringueyi and B. rubida. To further explore the types of neurotoxins present, venoms from the representatives B. armata, B. caudalis, B. cornuta and B. rubida were additionally tested in the chick biventer cervicis nerve muscle preparation, which showed presynaptic and postsynaptic activity for B. caudalis and only presynaptic neurotoxicity for B. cornuta and B. rubida, with myotoxicity also evident for some species. These results, combined with the biolayer interferometry results, indicate complex neurotoxicity exerted by Bitis species, which varies dramatically by lineage tested upon. Our data also further support the importance of sampling across geographical localities, as significant intraspecific variation of postsynaptic neurotoxicity was reported across the different localities.

Replacement and Parallel Simplification of Nonhomologous Proteinases Maintain Venom Phenotypes in Rear-Fanged Snakes

Novel phenotypes are commonly associated with gene duplications and neofunctionalization, less documented are the cases of phenotypic maintenance through the recruitment of novel genes. Proteolysis is the primary toxic character of many snake venoms, and ADAM metalloproteinases, named snake venom metalloproteinases (SVMPs), are largely recognized as the major effectors of this phenotype. However, by investigating original transcriptomes from 58 species of advanced snakes (Caenophidia) across their phylogeny, we discovered that a different enzyme, matrix metalloproteinase (MMP), is actually the dominant venom component in three tribes (Tachymenini, Xenodontini, and Conophiini) of rear-fanged snakes (Dipsadidae). Proteomic and functional analyses of these venoms further indicate that MMPs are likely playing an "SVMP-like" function in the proteolytic phenotype. A detailed look into the venom-specific sequences revealed a new highly expressed MMP subtype, named snake venom MMP (svMMP), which originated independently on at least three occasions from an endogenous MMP-9. We further show that by losing ancillary noncatalytic domains present in its ancestors, svMMPs followed an evolutionary path toward a simplified structure during their expansion in the genomes, thus paralleling what has been proposed for the evolution of their Viperidae counterparts, the SVMPs. Moreover, we inferred an inverse relationship between the expression of svMMPs and SVMPs along the evolutionary history of Xenodontinae, pointing out that one type of enzyme may be substituting for the other, whereas the general (metallo)proteolytic phenotype is maintained. These results provide rare evidence on how relevant phenotypic traits can be optimized via natural selection on nonhomologous genes, yielding alternate biochemical components.

Size Matters: An Evaluation of the Molecular Basis of Ontogenetic Modifications in the Composition of Bothrops jararacussu Snake Venom

Ontogenetic changes in venom composition have been described in Bothrops snakes, but only a few studies have attempted to identify the targeted paralogues or the molecular mechanisms involved in modifications of gene expression during ontogeny. In this study, we decoded B. jararacussu venom gland transcripts from six specimens of varying sizes and analyzed the variability in the composition of independent venom proteomes from 19 individuals. We identified 125 distinct putative toxin transcripts, and of these, 73 were detected in venom proteomes and only 10 were involved in the ontogenetic changes. Ontogenetic variability was linearly related to snake size and did not correspond to the maturation of the reproductive stage. Changes in the transcriptome were highly predictive of changes in the venom proteome. The basic myotoxic phospholipases A₂ (PLA₂s) were the most abundant components in larger snakes, while in venoms from smaller snakes, PIII-class SVMPs were the major components. The snake venom metalloproteinases (SVMPs) identified corresponded to novel sequences and conferred higher pro-coagulant and hemorrhagic functions to the venom of small snakes. The mechanisms modulating venom variability are predominantly related to transcriptional events and may consist of an advantage of higher hematotoxicity and more efficient predatory function in the venom from small snakes.

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

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

Widespread Evolution of Molecular Resistance to Snake Venom α-Neurotoxins in Vertebrates

Venomous snakes are important subjects of study in evolution, ecology, and biomedicine. Many venomous snakes have alpha-neurotoxins (α-neurotoxins) in their venom. These toxins bind the alpha-1 nicotinic acetylcholine receptor (nAChR) at the neuromuscular junction, causing paralysis and asphyxia. Several venomous snakes and their predators have evolved resistance to α-neurotoxins. The resistance is conferred by steric hindrance from N-glycosylated asparagines at amino acids 187 or 189, by an arginine at position 187 that has been hypothesized to either electrostatically repulse positively charged neurotoxins or sterically interfere with α-neurotoxin binding, or proline replacements at positions 194 or 197 of the nAChR ligand-binding domain to inhibit α-neurotoxin binding through structural changes in the receptor. Here, we analyzed this domain in 148 vertebrate species, and assessed its amino acid sequences for resistance-associated mutations. Of these sequences, 89 were sequenced de novo. We find widespread convergent evolution of the N-glycosylation form of resistance in several taxa including venomous snakes and their lizard prey, but not in the snake-eating birds studied. We also document new lineages with the arginine form of inhibition. Using an in vivo assay in four species, we provide further evidence that N-glycosylation mutations reduce the toxicity of cobra venom. The nAChR is of crucial importance for normal neuromuscular function and is highly conserved throughout the vertebrates as a result. Our research shows that the evolution of α-neurotoxins in snakes may well have prompted arms races and mutations to this ancient receptor across a wide range of sympatric vertebrates. These findings underscore the inter-connectedness of the biosphere and the ripple effects that one adaption can have across global ecosystems.

A triple WAP domain containing protein acts in antibacterial immunity of weather loach, Misgurnus anguillicaudatus

Whey acidic protein domain (WAPD) occurs in a variety of proteins in animals and many of WAPD-containing proteins are involved in immunity. In the present study, a novel protein containing three WAPDs was identified from the weather loach, Misgurnus anguillicaudatus, designated as MaTWD. MaTWD share high identity with TWDs from fish but low identity with TWDs from other animal phyla. MaTWD transcripts mainly distributed in gills and head kidney responded to bacterial challenge with significant upregulation. In vitro assay with recombinant MaTWD protein revealed that MaTWD had antiprotease activity against bacterial proteases. Moreover, MaTWD exhibited bacterial binding capacity and antimicrobial activity. Most importantly, exogenous MaTWD protected loach against bacterial infection by reducing loach mortality. We infer that MaTWD participates in the antibacterial immunity of loach via its antiprotease and antimicrobial activities.

A symmetry or asymmetry: Functional and compositional comparison of venom from the left and right glands of the Indochinese spitting cobra (Naja siamensis)

Contralaterally positioned maxillary (upper jaw) venom glands in snakes are mechanically independent, being able to discharge venom from either gland separately. This has led some studies to test venom function and composition of each contralaterally positioned venom gland to investigate any differences. However, the data on the subject to-date derives from limited sample sizes, appearing somewhat contradictory, and thus still remains inconclusive. Here, we tested samples obtained from the left and right venom glands of four N. siamensis specimens for their relative binding to the orthosteric site of amphibian, lizard, snake, bird, and rodent alpha-1 nicotinic acetylcholine receptors. We also show the relative proteomic patterns displayed by reversed phase liquid chromatography – mass spectrometry. Our results indicate that three of the venom gland sets showed no difference in both functional binding and composition, whilst one venom gland set showed a slight difference in functional binding (but not in specificity patterns between prey types) or venom composition. We hypothesise that these differences in functional binding may be due to one gland having previously ejected venom at some time prior to venom extraction, whilst its contralateral counterpart did not. This might cause the differential rate of toxin replenishment to be unequal between glands, thus instigating the difference in potency, likely due to uneven toxin proportions between glands at the time of venom extraction. These results demonstrate that the separate venom producing glands in snakes remain under the same genetic control elements and produce identical venom components.

•

No differences in venom profile by contralateral venom glands of Naja siamensis.

•

No differences in alpha-1 binding potency by contralateral N. siamensis venom glands.

•

A single outlier in the data showed a small difference in relative potency.

•

Potency is likely affected by differential replenishment rates.

What makes a fang? Phylogenetic and ecological controls on tooth evolution in rear-fanged snakes

BACKGROUND

Fangs are a putative key innovation that revolutionized prey capture and feeding in snakes, and - along with their associated venom phenotypes - have made snakes perhaps the most medically-significant vertebrate animals. Three snake clades are known for their forward-positioned fangs, and these clades (Elapidae, Viperidae, and Atractaspidinae) contain the majority of snakes that are traditionally considered venomous. However, many other snakes are "rear-fanged": they possess potentially venom-delivering teeth situated at the rear end of the upper jaw. Quantification of fang phenotypes - and especially those of rear-fanged species - has proved challenging or impossible owing to the small size and relative rarity of many such snakes. Consequently, it has been difficult to understand the evolutionary history of both venom and prey-capture strategies across extant snakes. We quantified variation in the dentition of 145 colubriform ("advanced") snake species using microCT scanning and compared dental characters with ecological data on species' diet and prey capture method(s) to understand broader patterns in snake fang evolution.

RESULTS

Dental traits such as maxilla length, tooth number, and fang size show strong phylogenetic signal across Colubriformes. We find extreme heterogeneity and evolutionary lability in the rear-fanged phenotype in colubrid (colubrine, dipsadine, and natricine lineages) and lamprophiid snakes, in contrast to relative uniformity in the front fanged phenotypes of other groups (vipers and, to a lesser extent, elapids). Fang size and position are correlated with venom-use in vipers, elapids, and colubrid snakes, with the latter group shifting fangs anteriorly by shortening the entire maxillary bone. We find that maxilla length and tooth number may also be correlated with the evolution of dietary specialization. Finally, an ancestral state reconstruction suggests that fang loss is a recurring phenomenon in colubrid snakes, likely accompanied by shifts in diet and prey capture mode.

CONCLUSIONS

Our study provides a framework for quantifying the complex morphologies associated with venom use in snakes. Our results suggest that fang phenotypes, and particularly the rear-fanged phenotype, in snakes are both diverse and labile, facilitating a wide range of ecological strategies and contributing to spectacular radiations of these organisms in tropical and subtropical biomes worldwide.

Toxinology provides multidirectional and multidimensional opportunities: A personal perspective

In nature, toxins have evolved as weapons to capture and subdue the prey or to counter predators or competitors. When they are inadvertently injected into humans, they cause symptoms ranging from mild discomfort to debilitation and death. Toxinology is the science of studying venoms and toxins that are produced by a wide variety of organisms. In the past, the structure, function and mechanisms of most abundant and/or most toxic components were characterized to understand and to develop strategies to neutralize their toxicity. With recent technical advances, we are able to evaluate and determine the toxin profiles using transcriptomes of venom glands and proteomes of tiny amounts of venom. Enormous amounts of data from these studies have opened tremendous opportunities in many directions of basic and applied research. The lower costs for profiling venoms will further fuel the expansion of toxin database, which in turn will provide greater exciting and bright opportunities in toxin research.

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

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

Macroevolutionary Analyses Suggest That Environmental Factors, Not Venom Apparatus, Play Key Role in Terebridae Marine Snail Diversification

How species diversification occurs remains an unanswered question in predatory marine invertebrates, such as sea snails of the family Terebridae. However, the anatomical disparity found throughput the Terebridae provides a unique perspective for investigating diversification patterns in venomous predators. In this study, a new dated molecular phylogeny of the Terebridae is used as a framework for investigating diversification of the family through time, and for testing the putative role of intrinsic and extrinsic traits, such as shell size, larval ecology, bathymetric distribution, and anatomical features of the venom apparatus, as drivers of terebrid species diversification. Macroevolutionary analysis revealed that when diversification rates do not vary across Terebridae clades, the whole family has been increasing its global diversification rate since 25 Ma. We recovered evidence for a concurrent increase in diversification of depth ranges, while shell size appeared to have undergone a fast divergence early in terebrid evolutionary history. Our data also confirm that planktotrophy is the ancestral larval ecology in terebrids, and evolutionary modeling highlighted that shell size is linked to larval ecology of the Terebridae, with species with long-living pelagic larvae tending to be larger and have a broader size range than lecithotrophic species. Although we recovered patterns of size and depth trait diversification through time and across clades, the presence or absence of a venom gland (VG) did not appear to have impacted Terebridae diversification. Terebrids have lost their venom apparatus several times and we confirm that the loss of a VG happened in phylogenetically clustered terminal taxa and that reversal is extremely unlikely. Our findings suggest that environmental factors, and not venom, have had more influence on terebrid evolution.

Development of a generic high-throughput screening assay for profiling snake venom protease activity after high-resolution chromatographic fractionation

Snakebites cause upwards of 1.8 million envenomings, 138,000 deaths and 500,000 cases of long term morbidity each year. Viper snake venoms (family Viperidae) generally contain a high proportion of proteases which can cause devastating effects such as hemorrhage, coagulopathy, edema, necrosis, and severe pain, in envenomed victims. In this study, analytical techniques were combined with enzymatic assays to develop a novel method for the detection of snake venom protease activity by using rhodamine-110-peptide substrate. In the so called at-line nanofractionation set up, crude venoms were first separated with reversed phase liquid chromatography, after which fractions were collected onto 384-well plates. Protease activity assays were then performed in the 384-well plates and bioassay chromatograms were constructed revealing protease activity. Parallel obtained UV absorbance, MS and proteomics data from a previous study facilitated toxin identification. The application of the rhodamine-110-peptide substrate assay showed significantly greater sensitivity compared to prior assays using casein-FITC as the substrate. Moreover, cross referencing UV and MS data and resulted in the detection of a number of tentative proteases suspected to exhibit protease activity, including snake venom serine proteases from Calloselasma rhodostoma and Daboia russelli venom and a snake venom metalloproteinase from the venom of Echis ocellatus. Our data demonstrate that his methodology can be a useful tool for selectively identifying snake venom proteases, and can be applied to provide a better understanding of protease-induced pathologies and the development of novel therapeutics for treating snakebite.

High throughput screening and identification of coagulopathic snake venom proteins and peptides using nanofractionation and proteomics approaches

Snakebite is a neglected tropical disease that results in a variety of systemic and local pathologies in envenomed victims and is responsible for around 138,000 deaths every year. Many snake venoms cause severe coagulopathy that makes victims vulnerable to suffering life-threating haemorrhage. The mechanisms of action of coagulopathic snake venom toxins are diverse and can result in both anticoagulant and procoagulant effects. However, because snake venoms consist of a mixture of numerous protein and peptide components, high throughput characterizations of specific target bioactives is challenging. In this study, we applied a combination of analytical and pharmacological methods to identify snake venom toxins from a wide diversity of snake species that perturb coagulation. To do so, we used a high-throughput screening approach consisting of a miniaturised plasma coagulation assay in combination with a venom nanofractionation approach. Twenty snake venoms were first separated using reversed-phase liquid chromatography, and a post-column split allowed a small fraction to be analyzed with mass spectrometry, while the larger fraction was collected and dispensed onto 384-well plates. After fraction collection, any solvent present in the wells was removed by means of freeze-drying, after which it was possible to perform a plasma coagulation assay in order to detect coagulopathic activity. Our results demonstrate that many snake venoms simultaneously contain both procoagulant and anticoagulant bioactives that contribute to coagulopathy. In-depth identification analysis from seven medically-important venoms, via mass spectrometry and nanoLC-MS/MS, revealed that phospholipase A2 toxins are frequently identified in anticoagulant venom fractions, while serine protease and metalloproteinase toxins are often associated with procoagulant bioactivities. The nanofractionation and proteomics approach applied herein seems likely to be a valuable tool for the rational development of next-generation snakebite treatments by facilitating the rapid identification and fractionation of coagulopathic toxins, thereby enabling specific targeting of these toxins by new therapeutics such as monoclonal antibodies and small molecule inhibitors.

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

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

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.

Alteration of Bumblebee Venom Composition toward Higher Elevation

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

Solenodon genome reveals convergent evolution of venom in eulipotyphlan mammals

Venom systems are key adaptations that have evolved throughout the tree of life and typically facilitate predation or defense. Despite venoms being model systems for studying a variety of evolutionary and physiological processes, many taxonomic groups remain understudied, including venomous mammals. Within the order Eulipotyphla, multiple shrew species and solenodons have oral venom systems. Despite morphological variation of their delivery systems, it remains unclear whether venom represents the ancestral state in this group or is the result of multiple independent origins. We investigated the origin and evolution of venom in eulipotyphlans by characterizing the venom system of the endangered Hispaniolan solenodon (Solenodon paradoxus). We constructed a genome to underpin proteomic identifications of solenodon venom toxins, before undertaking evolutionary analyses of those constituents, and functional assessments of the secreted venom. Our findings show that solenodon venom consists of multiple paralogous kallikrein 1 (KLK1) serine proteases, which cause hypotensive effects in vivo, and seem likely to have evolved to facilitate vertebrate prey capture. Comparative analyses provide convincing evidence that the oral venom systems of solenodons and shrews have evolved convergently, with the 4 independent origins of venom in eulipotyphlans outnumbering all other venom origins in mammals. We find that KLK1s have been independently coopted into the venom of shrews and solenodons following their divergence during the late Cretaceous, suggesting that evolutionary constraints may be acting on these genes. Consequently, our findings represent a striking example of convergent molecular evolution and demonstrate that distinct structural backgrounds can yield equivalent functions.

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

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

Clinical implications of coagulotoxic variations in Mamushi (Viperidae: Gloydius) snake venoms

Snake bite is currently one of the most neglected tropical diseases affecting much of the developing world. Asian pit vipers are responsible for a considerable amount of envenomations annually and bites can cause a multitude of clinical complications resulting from coagulopathic and neuropathic effects. While intense research has been undertaken for some species of Asian pit viper, functional coagulopathic effects have been neglected for others. We investigated their effects upon the human clotting cascade using venoms of four species of Gloydius and Ovophis okinavensis, a species closely to Gloydius. All species of included within this investigation displayed varying fibrinogenolytic effects, resulting in a net anticoagulant outcome. Gloydius saxatilis and Gloydius ussuriensis displayed the most variable effects from differing localities, sampled from Russia and Korea. As this Gloydius investigation includes some geographical variation, notable results indicate key variations of these species that point to possible limitations in antivenom cross-reactivities, which may have implications for the clinical care of victims envenomed by these snakes.

Snake Venoms in Drug Discovery: Valuable Therapeutic Tools for Life Saving

Animal venoms are used as defense mechanisms or to immobilize and digest prey. In fact, venoms are complex mixtures of enzymatic and non-enzymatic components with specific pathophysiological functions. Peptide toxins isolated from animal venoms target mainly ion channels, membrane receptors and components of the hemostatic system with high selectivity and affinity. The present review shows an up-to-date survey on the pharmacology of snake-venom bioactive components and evaluates their therapeutic perspectives against a wide range of pathophysiological conditions. Snake venoms have also been used as medical tools for thousands of years especially in tradition Chinese medicine. Consequently, snake venoms can be considered as mini-drug libraries in which each drug is pharmacologically active. However, less than 0.01% of these toxins have been identified and characterized. For instance, Captopril^® (Enalapril), Integrilin^® (Eptifibatide) and Aggrastat^® (Tirofiban) are drugs based on snake venoms, which have been approved by the FDA. In addition to these approved drugs, many other snake venom components are now involved in preclinical or clinical trials for a variety of therapeutic applications. These examples show that snake venoms can be a valuable source of new principle components in drug discovery.

Preliminary Biochemical and Venomic Characterization of the Venom of Phalotris lemniscatus (Serpentes, Colubridae)

BACKGROUND

For many decades, research on snake venom toxinology focused mainly on the venoms of Viperidae and Elapidae species, which were traditionally the only ones considered as venomous. However, much less interest has been given to the venom produced by opisthoglyphous colubrid snakes, since they were typically considered of no clinical relevance.

OBJECTIVE

The aim of this work is to perform a preliminary biochemical and venomic characterization of the venom of the colubrid snake Phalotris lemniscatus, a species that has been responsible for two relevant cases of envenomation in Uruguay.

METHODS

We extracted venom from collected specimens and performed different biochemical and proteomic assays to understand its toxin composition.

RESULTS

We found that the venom of P. lemniscatus is composed of protein families typically present in snake venoms, such as metallo and serine preoteases, L-amino acid oxidases, phospholipases A2s, Ctype lectines-like, Kunitz-type proteins and three-finger toxins. Activity assays demonstrated a highly active gelatinolytic component as well as a potent capability to induce blood coagulation.

CONCLUSION

The results indicate that the venom of P. lemniscatus contains hemotoxic activities and components that resemble those found in Viperidae (Bothrops) snakes and that can induce a clinically relevant accident. Further studies are needed to better understand the venom composition of this colubrid snake and its most active compounds.

Venom composition of adult Western Diamondback Rattlesnakes (Crotalus atrox) maintained under controlled diet and environmental conditions shows only minor changes

Many species of snakes produce venom as a chemical means of procuring potentially fractious prey. Studies have increasingly focused on venom compositional variation between and within individual snakes of the same species/subspecies, with significant differences often being observed. This variation in composition has been attributed to differences in age, season, diet, and environment, suggesting that these factors could help explain the inter- and intra-specific variation found in some snake venoms, perhaps via some type of feedback mechanism(s). To address several of these possible sources of variation, this study utilized wild-caught Western Diamondback Rattlesnakes (Crotalus atrox) from Cochise Co., AZ. Sixteen adult C. atrox were maintained in the lab on a diet of NSA mice for eight months to determine whether venom composition changed in captivity under a static diet in a stable environment. Reducing 1-D SDS-PAGE, fibrinogen degradation assays, reversed-phase HPLC, and MALDI-TOF mass spectrometry revealed only minor differences over time within individuals. Venom L-amino acid oxidase (LAAO) and phosphodiesterase activities significantly increased over the course of captivity, with no changes occurring in azocasein metalloproteinase, kallikrein-like serine proteinase (KLSP), or thrombin-like serine proteinase (TLSP) activities. Snake total length was positively correlated with TLSP activity and negatively correlated with LAAO and KLSP activity. There was typically a much higher degree of variation between individuals than within individuals for all analyses performed and measurements collected. Because the overall "fingerprint" of each snake's venom remained more/less constant, it is concluded that biologically significant changes in venom composition did not occur within individual C. atrox as a function of captivity/diet. However, this study does indicate that differences in activity levels do occur in minor venom enzyme components, but the differences observed are likely to be of minimal significance to the production of antivenom or to subsequent treatment of human envenomations.

Serum-based inhibition of pitviper venom by eastern indigo snakes (Drymarchon couperi)

When organisms possess chemical defenses, their predators may eventually evolve resistance to their toxins. Eastern indigo snakes (Drymarchon couperi; EIS) prey on pitvipers and are suspected to possess physiological resistance to their venom. In this study, we formally investigated this hypothesis using microassays that measured the ability of EIS blood sera to inhibit (A) hemolytic and (B) snake venom metalloproteinase (SVMP) activity of copperhead (Agkistrodon contortrix) venom. To serve as controls, we also tested the inhibitory ability of sera from house mice (Mus musculus) and checkered gartersnakes (Thamnophis marcianus), a snake that does not feed on pitvipers. Sera from both EIS and gartersnakes inhibited over 60% of SVMP activity, while only EIS sera also inhibited venom hemolytic activity (78%). Our results demonstrate that EIS serum is indeed capable of inhibiting two of the primary classes of toxins found in copperhead venom, providing the first empirical evidence suggesting that EIS possess physiological resistance to venom upon injection. Because we documented resistance to hemolytic components of pitviper venom within EIS but not gartersnakes, we speculate this resistance may be driven by selection from feeding on pitvipers while resistance to SVMP may be relatively widespread among snakes.

Summary: We confirm eastern indigo snakes possess a physiological resistance to pitviper venom and speculate that venom resistance in snakes may be a product of antagonistic interactions and phylogenetically conserved traits.

Coagulotoxic effects by brown snake (Pseudonaja) and taipan (Oxyuranus) venoms, and the efficacy of a new antivenom

Snakebite is a neglected tropical disease that disproportionately affects the poor. Antivenom is the only specific and effective treatment for snakebite, but its distribution is severely limited by several factors, including the prohibitive cost of some products. Papua New Guinea (PNG) is a snakebite hotspot but the high costs of Australian antivenoms (thousands of dollars per treatment) makes it unaffordable in PNG. A more economical taipan antivenom has recently been developed at the Instituto Clodomiro Picado (ICP) in Costa Rica for PNG and is currently undergoing clinical trials for the treatment of envenomations by coastal taipans (Oxyuranus scutellatus). In addition to potentially having the capacity to neutralise the effects of envenomations of non-PNG taipans, this antivenom may have the capacity to neutralise coagulotoxins in venom from closely related brown snakes (Pseudonaja spp.) also found in PNG. Consequently, we investigated the cross-reactivity of taipan antivenom across the venoms of all Oxyuranus and Pseudonaja species. In addition, to ascertain differences in venom biochemistry that influence variation in antivenom efficacy, we tested for relative cofactor dependence. We found that the new ICP taipan antivenom exhibited high selectivity for Oxyuranus venoms and only low to moderate cross-reactivity with any Pseudonaja venoms. Consistent with this genus level distinction in antivenom efficacy were fundamental differences in the venom biochemistry. Not only were the Pseudonaja venoms significantly more procoagulant, but they were also much less dependent upon the cofactors calcium and phospholipid. There was a strong correlation between antivenom efficacy, clotting time and cofactor dependence. This study sheds light on the structure-function relationships of the procoagulant toxins within these venoms and may have important clinical implications including for the design of next-generation antivenoms.

A Single Mutation Unlocks Cascading Exaptations in the Origin of a Potent Pitviper Neurotoxin

Evolutionary innovations and complex phenotypes seemingly require an improbable amount of genetic change to evolve. Rattlesnakes display two dramatically different venom phenotypes. Type I venoms are hemorrhagic with low systemic toxicity and high expression of tissue-destroying snake venom metalloproteinases. Type II venoms are highly neurotoxic and lack snake venom metalloproteinase expression and associated hemorrhagic activity. This dichotomy hinges on Mojave toxin (MTx), a phospholipase A2 (PLA2) based β-neurotoxin expressed in Type II venoms. MTx is comprised of a nontoxic acidic subunit that undergoes extensive proteolytic processing and allosterically regulates activity of a neurotoxic basic subunit. Evolution of the acidic subunit presents an evolutionary challenge because the need for high expression of a nontoxic venom component and the proteolytic machinery required for processing suggests genetic changes of seemingly little immediate benefit to fitness. We showed that MTx evolved through a cascading series of exaptations unlocked by a single nucleotide change. The evolution of one new cleavage site in the acidic subunit unmasked buried cleavage sites already present in ancestral PLA2s, enabling proteolytic processing. Snake venom serine proteases, already present in the venom to disrupt prey hemostasis, possess the requisite specificities for MTx acidic subunit proteolysis. The dimerization interface between MTx subunits evolved by exploiting a latent, but masked, hydrophobic interaction between ancestral PLA2s. The evolution of MTx through exaptation of existing functional and structural features suggests complex phenotypes that depend on evolutionary innovations can arise from minimal genetic change enabled by prior evolution.

Habu coagulotoxicity: Clinical implications of the functional diversification of Protobothrops snake venoms upon blood clotting factors

Venom can affect any part of the body reachable via the bloodstream. Toxins which specifically act upon the coagulation cascade do so either by anticoagulant or procoagulant mechanisms. Here we investigated the coagulotoxic effects of six species within the medically important pit viper genus Protobothrops (Habu) from the Chinese mainland and Japanese islands, a genus known to produce hemorrhagic shock in envenomed patients. Differential coagulotoxicity was revealed: P. jerdonii and P. mangshanensis produced an overall net anticoagulant effect through the pseudo-procoagulant clotting of fibrinogen; P. flavoviridis and P. tokarensis exhibit a strong anticoagulant activity through the destructive cleavage of fibrinogen; and while P. elegans and P. mucrosquamatus both cleaved the A-alpha and B-beta chains of fibrinogen they did not exhibit strong anticoagulant activity. These variations in coagulant properties were congruent with phylogeny, with the closest relatives exhibiting similar venom effects in their action upon fibrinogen. Ancestral state reconstruction indicated that anticoagulation mediated by pseudo-procoagulant cleavage of fibrinogen is the basal state, while anticoagulation produced by destructive cleavage of fibrinogen is the derived state within this genus. This is the first in depth study of its kind highlighting extreme enzymatic variability, functional diversification and clotting diversification within one genus surrounding one target site, governed by variability in co-factor dependency. The documentation that the same net overall function, anticoagulation, is mediated by differential underlying mechanics suggests limited antivenom cross-reactivity, although this must be tested in future work. These results add to the body of knowledge necessary to inform clinical management of the envenomed patient.

Biochemical characterization of venom from Pseudoboa neuwiedii (Neuwied's false boa; Xenodontinae; Pseudoboini)

In this work, we examined the proteolytic and phospholipase A₂ (PLA₂) activities of venom from the opisthoglyphous colubrid Pseudoboa neuwiedii. Proteolytic activity (3 and 10 μg of venom) was comparable to that of Bothrops neuwiedii venom but less than Bothrops atrox. This activity was inhibited by EDTA and 1,10-phenanthroline but only slightly affected (≤30% inhibition) by PMSF and AEBSF, indicating it was mediated by snake venom metalloproteinases (SVMPs). The pH and temperature optima for proteolytic activity were 8.0 and 37 °C, respectively. The venom had no esterase activity, whereas PLA₂ activity was similar to B. atrox, greater than B. neuwiedii but less than B. jararacussu. SDS-PAGE revealed venom proteins >100 kDa, 45-70 kDa, 21-24 kDa and ~15 kDa, and mass spectrometry of protein bands revealed SVMPs, cysteine-rich secretory proteins (CRISPs) and PLA₂, but no serine proteinases. In gelatin zymography, the most active bands occurred at 65-68 kDa (seen with 0.05-0.25 μg of venom). Caseinolytic activity occurred at 50-66 kDa and was generally weaker than gelatinolytic activity. RP-HPLC of venom yielded 15 peaks, five of which showed gelatinolytic activity; peak 7 was the most active and apparently contained a P-III class SVMP. The venom showed α-fibrinogenase activity, without affecting the β and γ chains; this activity was inhibited by EDTA and 1,10-phenanthroline. The venom did not clot rat citrated plasma but reduced the rate and extent of coagulation after plasma recalcification. In conclusion, P. neuwiedii venom is highly proteolytic and could potentially affect coagulation in vivo by degrading fibrinogen via SVMPs.

Coagulotoxicity of Bothrops (Lancehead Pit-Vipers) Venoms from Brazil: Differential Biochemistry and Antivenom Efficacy Resulting from Prey-Driven Venom Variation

Lancehead pit-vipers (Bothrops genus) are an extremely diverse and medically important group responsible for the greatest number of snakebite envenomations and deaths in South America. Bothrops atrox (common lancehead), responsible for majority of snakebites and related deaths within the Brazilian Amazon, is a highly adaptable and widely distributed species, whose venom variability has been related to several factors, including geographical distribution and habitat type. This study examined venoms from four B. atrox populations (Belterra and Santarém, PA; Pres. Figueiredo, AM and São Bento, MA), and two additional Bothrops species (B. jararaca and B. neuwiedi) from Southeastern region for their coagulotoxic effects upon different plasmas (human, amphibian, and avian). The results revealed inter– and intraspecific variations in coagulotoxicity, including distinct activities between the three plasmas, with variations in the latter two linked to ecological niche occupied by the snakes. Also examined were the correlated biochemical mechanisms of venom action. Significant variation in the relative reliance upon the cofactors calcium and phospholipid were revealed, and the relative dependency did not significantly correlate with potency. Relative levels of Factor X or prothrombin activating toxins correlated with prey type and prey escape potential. The antivenom was shown to perform better in neutralising prothrombin activation activity than neutralising Factor X activation activity. Thus, the data reveal new information regarding the evolutionary selection pressures shaping snake venom evolution, while also having significant implications for the treatment of the envenomed patient. These results are, therefore, an intersection between evolutionary biology and clinical medicine.

Three-Finger Toxin Diversification in the Venoms of Cat-Eye Snakes (Colubridae: Boiga)

The Asian genus Boiga (Colubridae) is among the better studied non-front-fanged snake lineages, because their bites have minor, but noticeable, effects on humans. Furthermore, B. irregularis has gained worldwide notoriety for successfully invading Guam and other nearby islands with drastic impacts on the local bird populations. One of the factors thought to allow B. irregularis to become such a noxious pest is irditoxin, a dimeric neurotoxin composed of two three-finger toxins (3FTx) joined by a covalent bond between two newly evolved cysteines. Irditoxin is highly toxic to diapsid (birds and reptiles) prey, but roughly 1000 × less potent to synapsids (mammals). Venom plays an important role in the ecology of all species of Boiga, but it remains unknown if any species besides B. irregularis produce irditoxin-like dimeric toxins. In this study, we use transcriptomic analyses of venom glands from five species [B. cynodon, B. dendrophila dendrophila, B. d. gemmicincta, B. irregularis (Brisbane population), B. irregularis (Sulawesi population), B. nigriceps, B. trigonata] and proteomic analyses of B. d. dendrophila and a representative of the sister genus Toxicodryas blandingii to investigate the evolutionary history of 3FTx within Boiga and its close relative. We found that 92.5% of Boiga 3FTx belong to a single clade which we refer to as denmotoxin-like because of the close relation between these toxins and the monomeric denmotoxin according to phylogenetic, sequence clustering, and protein similarity network analyses. We show for the first time that species beyond B. irregularis secrete 3FTx with additional cysteines in the same position as both the A and B subunits of irditoxin. Transcripts with the characteristic mutations are found in B. d. dendrophila, B. d. gemmicincta, B. irregularis (Brisbane population), B. irregularis (Sulawesi population), and B. nigriceps. These results are confirmed by proteomic analyses that show direct evidence of dimerization within the venom of B. d. dendrophila, but not T. blandingii. Our results also suggest the possibility of novel dimeric toxins in other genera such as Telescopus and Trimorphodon. All together, this suggests that the origin of these peculiar 3FTx is far earlier than was appreciated and their evolutionary history has been complex.

Does size matter? Venom proteomic and functional comparison between night adder species (Viperidae: Causus) with short and long venom glands

Night adders (Causus species within the Viperidae family) are amphibian specialists and a common source of snakebite in Africa. Some species are unique in that they have the longest venom glands of any viper, extending approximately 10% of the body length. Despite their potential medical importance and evolutionary novelty, their venom has received almost no research attention. In this study, venoms from a short-glanded species (C. lichtensteinii) and from a long-glanded species (C. rhombeatus) were compared using a series of proteomic and bioactivity testing techniques to investigate and compare the toxin composition and functioning of the venoms of these two species. Both C. rhombeatus and C. lichtensteinii were similar in overall venom composition and inhibition of blood coagulation through non-clotting proteolytic cleavage of fibrinogen. While the 1D gel profiles were very similar to each other in the toxin types present, 2D gel analyses revealed isoformic differences within each toxin classes. This variation was congruent with differential efficacy of South African Institute for Medical Research snake polyvalent antivenom, with C. lichtensteinii unaffected at the dose tested while C. rhombeatus was moderately but significantly neutralized. Despite the variation within toxin classes, the similarity in overall venom biochemistry suggests that the selection pressure for the evolution of long glands served to increase venom yield in order to subjugate proportionally large anurans as a unique form of niche partitioning, and is not linked to significant changes in venom function. These results not only contribute to the body of venom evolution knowledge but also highlight the limited clinical management outcomes for Causus envenomations.