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.

Distinct regulatory networks control toxin gene expression in elapid and viperid snakes

BACKGROUND

Venom systems are ideal models to study genetic regulatory mechanisms that underpin evolutionary novelty. Snake venom glands are thought to share a common origin, but there are major distinctions between venom toxins from the medically significant snake families Elapidae and Viperidae, and toxin gene regulatory investigations in elapid snakes have been limited. Here, we used high-throughput RNA-sequencing to profile gene expression and microRNAs between active (milked) and resting (unmilked) venom glands in an elapid (Eastern Brown Snake, Pseudonaja textilis), in addition to comparative genomics, to identify cis- and trans-acting regulation of venom production in an elapid in comparison to viperids (Crotalus viridis and C. tigris).

RESULTS

Although there is conservation in high-level mechanistic pathways regulating venom production (unfolded protein response, Notch signaling and cholesterol homeostasis), there are differences in the regulation of histone methylation enzymes, transcription factors, and microRNAs in venom glands from these two snake families. Histone methyltransferases and transcription factor (TF) specificity protein 1 (Sp1) were highly upregulated in the milked elapid venom gland in comparison to the viperids, whereas nuclear factor I (NFI) TFs were upregulated after viperid venom milking. Sp1 and NFI cis-regulatory elements were common to toxin gene promoter regions, but many unique elements were also present between elapid and viperid toxins. The presence of Sp1 binding sites across multiple elapid toxin gene promoter regions that have been experimentally determined to regulate expression, in addition to upregulation of Sp1 after venom milking, suggests this transcription factor is involved in elapid toxin expression. microRNA profiles were distinctive between milked and unmilked venom glands for both snake families, and microRNAs were predicted to target a diversity of toxin transcripts in the elapid P. textilis venom gland, but only snake venom metalloproteinase transcripts in the viperid C. viridis venom gland. These results suggest differences in toxin gene posttranscriptional regulation between the elapid P. textilis and viperid C. viridis.

CONCLUSIONS

Our comparative transcriptomic and genomic analyses between toxin genes and isoforms in elapid and viperid snakes suggests independent toxin regulation between these two snake families, demonstrating multiple different regulatory mechanisms underpin a venomous phenotype.

Prevalent bee venom genes evolved before the aculeate stinger and eusociality

BACKGROUND

Venoms, which have evolved numerous times in animals, are ideal models of convergent trait evolution. However, detailed genomic studies of toxin-encoding genes exist for only a few animal groups. The hyper-diverse hymenopteran insects are the most speciose venomous clade, but investigation of the origin of their venom genes has been largely neglected.

RESULTS

Utilizing a combination of genomic and proteo-transcriptomic data, we investigated the origin of 11 toxin genes in 29 published and 3 new hymenopteran genomes and compiled an up-to-date list of prevalent bee venom proteins. Observed patterns indicate that bee venom genes predominantly originate through single gene co-option with gene duplication contributing to subsequent diversification.

CONCLUSIONS

Most Hymenoptera venom genes are shared by all members of the clade and only melittin and the new venom protein family anthophilin1 appear unique to the bee lineage. Most venom proteins thus predate the mega-radiation of hymenopterans and the evolution of the aculeate stinger.

Domain loss enabled evolution of novel functions in the snake three-finger toxin gene superfamily

Three-finger toxins (3FTXs) are a functionally diverse family of toxins, apparently unique to venoms of caenophidian snakes. Although the ancestral function of 3FTXs is antagonism of nicotinic acetylcholine receptors, redundancy conferred by the accumulation of duplicate genes has facilitated extensive neofunctionalization, such that derived members of the family interact with a range of targets. 3FTXs are members of the LY6/UPAR family, but their non-toxin ancestor remains unknown. Combining traditional phylogenetic approaches, manual synteny analysis, and machine learning techniques (including AlphaFold2 and ProtT5), we have reconstructed a detailed evolutionary history of 3FTXs. We identify their immediate ancestor as a non-secretory LY6, unique to squamate reptiles, and propose that changes in molecular ecology resulting from loss of a membrane-anchoring domain and changes in gene expression, paved the way for the evolution of one of the most important families of snake toxins.

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.

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.

Molecular mechanisms underlying the N-methyl-d-aspartate receptor antagonists: Highlighting their potential for transdiagnostic therapeutics

The so-called "psychedelic renaissance" has stimulated expanded interest in several classes of drugs that appear to possess transdiagnostic effects in the treatment of mental health disorders, specifically. N-methyl-d-aspartate receptor (NMDAR) antagonists are one such class with diverse therapeutic potential. NMDARs mediate excitatory postsynaptic signalling in the central nervous system (CNS) and are integral to normal neurobiological processes including neuronal development, synaptic transmission, and plasticity, and thus involved in learning and memory. However, NMDAR hyper-function is also implicated in acute CNS trauma, neuropsychiatric and neurodegenerative disorders, as well as chronic pain. The complex structure of NMDARs permits several locations for therapeutic inhibition, making these receptors a potential target for multiple drugs which modulate them in different ways. NMDAR antagonists, which may be competitive, non-competitive, or uncompetitive, either block glutamate from binding the receptor or modulate the response to glutamate binding. Despite longstanding concerns about side effects of NMDAR antagonists, recent research suggests that, when appropriately used, these agents have favourable safety profiles. Furthermore, their fast-acting mechanism of action, resulting in rapid effects compared to other therapeutic agents, makes them a promising class of drugs that may yield effective therapeutics for multiple CNS disorders.

Contextual Constraints: Dynamic Evolution of Snake Venom Phospholipase A2

Venom is a dynamic trait that has contributed to the success of numerous organismal lineages. Predominantly composed of proteins, these complex cocktails are deployed for predation and/or self-defence. Many non-toxic physiological proteins have been convergently and recurrently recruited by venomous animals into their toxin arsenal. Phospholipase A₂ (PLA₂) is one such protein and features in the venoms of many organisms across the animal kingdom, including snakes of the families Elapidae and Viperidae. Understanding the evolutionary history of this superfamily would therefore provide insight into the origin and diversification of venom toxins and the evolution of novelty more broadly. The literature is replete with studies that have identified diversifying selection as the sole influence on PLA₂ evolution. However, these studies have largely neglected the structural/functional constraints on PLA₂s, and the ecology and evolutionary histories of the diverse snake lineages that produce them. By considering these crucial factors and employing evolutionary analyses integrated with a schema for the classification of PLA₂s, we uncovered lineage-specific differences in selection regimes. Thus, our work provides novel insights into the evolution of this major snake venom toxin superfamily and underscores the importance of considering the influence of evolutionary and ecological contexts on molecular evolution.

Role of Phospholipases A2 in Vascular Relaxation and Sympatholytic Effects of Five Australian Brown Snake, Pseudonaja spp., Venoms in Rat Isolated Tissues

Human envenoming by Australian brown snakes (Pseudonaja spp.) may result in potentially life-threatening hypotension and subsequent cardiovascular collapse. There have been relatively few studies of the cardiovascular and sympathetic effects of Pseudonaja spp. venoms. In this study, we have examined the effects of venom from five brown snake species-P. affinis, aspidorhyncha, inframacula, nuchalis, and textilis-on cardiac inotropic and chronotropic responses, vascular tone, and sympathetic nerve-induced vascular contractions in rat isolated tissues. The role of phospholipases A₂ (PLA₂s) in venom-induced effects was assessed with the sPLA₂ inhibitor varespladib. In rat isolated left and right atria, there were no physiologically relevant effects of Pseudonaja venoms (0.1-30 µg/ml) on left atrial force of contraction (inotropy) or right atrial rate (chronotropy). In contrast, in isolated small mesenteric arteries precontracted with a thromboxane mimetic, each of the five brown snake venoms (at 30 µg/ml) caused marked vasorelaxation (-60 to -90% of contractile tone). Pretreatment with varespladib (1 µM) significantly inhibited the vasorelaxation caused by P. aspidorhyncha, P. nuchalis, and P. textilis venoms. Electrically induced sympathetic nerve-mediated contractions of mesenteric arteries were significantly attenuated by only P. textilis, and P. affinis venoms (30 µg/ml) and these sympatholytic effects were inhibited by varespladib (1 µM). Based on their inhibition with the sPLA₂ inhibitor varespladib, we conclude that PLA₂ toxins in P. aspidorhyncha, P. nuchalis, and P. textilis venoms are involved in brown snake venom-induced vasorelaxation and the sympatholytic effects of P. affinis, and P. textilis venoms. Our study supports the promising potential role of varespladib as an initial (pre-referral) and/or adjunct (in combination with antivenom) therapeutic agent for brown snake envenoming.

How the Toxin got its Toxicity

Venom systems are functional and ecological traits, typically used by one organism to subdue or deter another. A predominant subset of their constituent molecules—“toxins”—share this ecological function and are therefore molecules that mediate interactions between organisms. Such molecules have been referred to as “exochemicals.” There has been debate within the field of toxinology concerning the evolutionary pathways leading to the “recruitment” of a gene product for a toxic role within venom. We review these discussions and the evidence interpreted in support of alternate pathways, along with many of the most popular models describing the origin of novel molecular functions in general. We note that such functions may arise with or without gene duplication occurring and are often the consequence of a gene product encountering a novel “environment,” i.e., a range of novel partners for molecular interaction. After stressing the distinction between “activity” and “function,” we describe in detail the results of a recent study which reconstructed the evolutionary history of a multigene family that has been recruited as a toxin and argue that these results indicate that a pluralistic approach to understanding the origin of novel functions is advantageous. This leads us to recommend that an expansive approach be taken to the definition of “neofunctionalization”—simply the origins of a novel molecular function by any process—and “recruitment”—the “weaponization” of a molecule via the acquisition of a toxic function in venom, by any process. Recruitment does not occur at the molecular level or even at the level of gene expression, but only when a confluence of factors results in the ecological deployment of a physiologically active molecule as a toxin. Subsequent to recruitment, the evolutionary regime of a gene family may shift into a more dynamic form of “birth-and-death.” Thus, recruitment leads to a form of “downwards causation,” in which a change at the ecological level at which whole organisms interact leads to a change in patterns of evolution at the genomic level.

Causes and Consequences of Snake Venom Variation

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

Functional venomics of the Big-4 snakes of Pakistan

In South Asia, the "Big-4" venomous snakes Naja naja, Bungarus caeruleus, Daboia russelii, and Echis carinatus are so-called because they are the most medically important snakes in the region. Antivenom is the only effective treatment option for snakebite envenoming but antivenom is not produced domestically in Pakistan making the country reliant on polyvalent products imported from India and Saudi Arabia. The present study investigated the toxin composition and activity of the venoms of Pakistani specimens by means of proteomic and physio/pharmacological experiments. To evaluate the composition of venoms, 1D/2D-PAGE of crude venoms and RP-HPLC followed by SDS-PAGE were performed. Enzymatic, hemolytic, coagulant and platelet aggregating activities of crude venoms were assayed and were concordant with expectations based on the abundance of protein species in each. Neutralization assays were performed using Bharat polyvalent antivenom (BPAV), a product raised against venoms from Big-4 specimens from southern India. BPAV exhibited cross-reactivity against the Pakistani venoms, however, neutralization of clinically relevant activities was variable and rarely complete. Cumulatively, the presented data not only highlight geographical variations present in the venoms of the Big-4 snakes of South Asia, but also demonstrate the neutralization potential of Indian polyvalent against the venom of Pakistani specimens. Given the partial neutralization observed, it is clear that whilst BPAV is a life-saving product in Pakistan, in future it is hoped that a region-specific product might be manufactured domestically, using venoms of local snakes in the immunising mixture.

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.

Proteomic and functional variation within black snake venoms (Elapidae: Pseudechis)

Pseudechis (black snakes) is an Australasian elapid snake genus that inhabits much of mainland Australia, with two representatives confined to Papua New Guinea. The present study is the first to analyse the venom of all 9 described Pseudechis species (plus one undescribed species) to investigate the evolution of venom composition and functional activity. Proteomic results demonstrated that the typical Pseudechis venom profile is dominated by phospholipase A₂ toxins. Strong cytotoxicity was the dominant function for most species. P. porphyriacus, the most basal member of the genus, also exhibited the most divergent venom composition, being the only species with appreciable amounts of procoagulant toxins. The relatively high presence of factor Xa recovered in P. porphyriacus venom may be related to a predominantly amphibian diet. Results of this study provide important insights to guide future ecological and toxinological investigations.

Rattling the border wall: Pathophysiological implications of functional and proteomic venom variation between Mexican and US subspecies of the desert rattlesnake Crotalus scutulatus

While some US populations of the Mohave rattlesnake (Crotalus scutulatus scutulatus) are infamous for being potently neurotoxic, the Mexican subspecies C. s. salvini (Huamantlan rattlesnake) has been largely unstudied beyond crude lethality testing upon mice. In this study we show that at least some populations of this snake are as potently neurotoxic as its northern cousin. Testing of the Mexican antivenom Antivipmyn showed a complete lack of neutralisation for the neurotoxic effects of C. s. salvini venom, while the neurotoxic effects of the US subspecies C. s. scutulatus were time-delayed but ultimately not eliminated. These results document unrecognised potent neurological effects of a Mexican snake and highlight the medical importance of this subspecies, a finding augmented by the ineffectiveness of the Antivipmyn antivenom. These results also influence our understanding of the venom evolution of Crotalus scutulatus, suggesting that neurotoxicity is the ancestral feature of this species, with the US populations which lack neurotoxicity being derived states.

Catch a tiger snake by its tail: Differential toxicity, co-factor dependence and antivenom efficacy in a procoagulant clade of Australian venomous snakes

A paradigm of venom research is adaptive evolution of toxins as part of a predator-prey chemical arms race. This study examined differential co-factor dependence, variations relative to dietary preference, and the impact upon relative neutralisation by antivenom of the procoagulant toxins in the venoms of a clade of Australian snakes. All genera were characterised by venoms rich in factor Xa which act upon endogenous prothrombin. Examination of toxin sequences revealed an extraordinary level of conservation, which indicates that adaptive evolution is not a feature of this toxin type. Consistent with this, the venoms did not display differences on the plasma of different taxa. Examination of the prothrombin target revealed endogenous blood proteins are under extreme negative selection pressure for diversification, this in turn puts a strong negative selection pressure upon the toxins as sequence diversification could result in a drift away from the target. Thus this study reveals that adaptive evolution is not a consistent feature in toxin evolution in cases where the target is under negative selection pressure for diversification. Consistent with this high level of toxin conservation, the antivenom showed extremely high-levels of cross-reactivity. There was however a strong statistical correlation between relative degree of phospholipid-dependence and clotting time, with the least dependent venoms producing faster clotting times than the other venoms even in the presence of phospholipid. The results of this study are not only of interest to evolutionary and ecological disciplines, but also have implications for clinical toxinology.

Correlation between ontogenetic dietary shifts and venom variation in Australian brown snakes (Pseudonaja)

Venom is a key evolutionary trait, as evidenced by its widespread convergent evolution across the animal kingdom. In an escalating prey-predator arms race, venoms evolve rapidly to guarantee predatory or defensive success. Variation in venom composition is ubiquitous among snakes. Here, we tested variation in venom activity on substrates relevant to blood coagulation among Pseudonaja (brown snake) species, Australian elapids responsible for the majority of medically important human envenomations in Australia. A functional approach was employed to elucidate interspecific variation in venom activity in all nine currently recognised species of Pseudonaja. Fluorometric enzymatic activity assays were performed to test variation in whole venom procoagulant activity among species. Analyses confirmed the previously documented ontogenetic shift from non-coagulopathic venom in juveniles to coagulopathic venom as adults, except for the case of P. modesta, which retains non-coagulopathic venom as an adult. These shifts in venom activity correlate with documented ontogenetic shifts in diet among brown snakes from specialisation on reptilian prey as juveniles (and throughout the life cycle of P. modesta), to a more generalised diet in adults that includes mammals. The results of this study bring to light findings relevant to both clinical and evolutionary toxinology.

How the Cobra Got Its Flesh-Eating Venom: Cytotoxicity as a Defensive Innovation and Its Co-Evolution with Hooding, Aposematic Marking, and Spitting

The cytotoxicity of the venom of 25 species of Old World elapid snake was tested and compared with the morphological and behavioural adaptations of hooding and spitting. We determined that, contrary to previous assumptions, the venoms of spitting species are not consistently more cytotoxic than those of closely related non-spitting species. While this correlation between spitting and non-spitting was found among African cobras, it was not present among Asian cobras. On the other hand, a consistent positive correlation was observed between cytotoxicity and utilisation of the defensive hooding display that cobras are famous for. Hooding and spitting are widely regarded as defensive adaptations, but it has hitherto been uncertain whether cytotoxicity serves a defensive purpose or is somehow useful in prey subjugation. The results of this study suggest that cytotoxicity evolved primarily as a defensive innovation and that it has co-evolved twice alongside hooding behavior: once in the Hemachatus + Naja and again independently in the king cobras (Ophiophagus). There was a significant increase of cytotoxicity in the Asian Naja linked to the evolution of bold aposematic hood markings, reinforcing the link between hooding and the evolution of defensive cytotoxic venoms. In parallel, lineages with increased cytotoxicity but lacking bold hood patterns evolved aposematic markers in the form of high contrast body banding. The results also indicate that, secondary to the evolution of venom rich in cytotoxins, spitting has evolved three times independently: once within the African Naja, once within the Asian Naja, and once in the Hemachatus genus. The evolution of cytotoxic venom thus appears to facilitate the evolution of defensive spitting behaviour. In contrast, a secondary loss of cytotoxicity and reduction of the hood occurred in the water cobra Naja annulata, which possesses streamlined neurotoxic venom similar to that of other aquatic elapid snakes (e.g., hydrophiine sea snakes). The results of this study make an important contribution to our growing understanding of the selection pressures shaping the evolution of snake venom and its constituent toxins. The data also aid in elucidating the relationship between these selection pressures and the medical impact of human snakebite in the developing world, as cytotoxic cobras cause considerable morbidity including loss-of-function injuries that result in economic and social burdens in the tropics of Asia and sub-Saharan Africa.

The Snake with the Scorpion's Sting: Novel Three-Finger Toxin Sodium Channel Activators from the Venom of the Long-Glanded Blue Coral Snake (Calliophis bivirgatus)

Millions of years of evolution have fine-tuned the ability of venom peptides to rapidly incapacitate both prey and potential predators. Toxicofera reptiles are characterized by serous-secreting mandibular or maxillary glands with heightened levels of protein expression. These glands are the core anatomical components of the toxicoferan venom system, which exists in myriad points along an evolutionary continuum. Neofunctionalisation of toxins is facilitated by positive selection at functional hotspots on the ancestral protein and venom proteins have undergone dynamic diversification in helodermatid and varanid lizards as well as advanced snakes. A spectacular point on the venom system continuum is the long-glanded blue coral snake (Calliophis bivirgatus), a specialist feeder that preys on fast moving, venomous snakes which have both a high likelihood of prey escape but also represent significant danger to the predator itself. The maxillary venom glands of C. bivirgatus extend one quarter of the snake's body length and nestle within the rib cavity. Despite the snake's notoriety its venom has remained largely unstudied. Here we show that the venom uniquely produces spastic paralysis, in contrast to the flaccid paralysis typically produced by neurotoxic snake venoms. The toxin responsible, which we have called calliotoxin (δ-elapitoxin-Cb1a), is a three-finger toxin (3FTx). Calliotoxin shifts the voltage-dependence of NaV1.4 activation to more hyperpolarised potentials, inhibits inactivation, and produces large ramp currents, consistent with its profound effects on contractile force in an isolated skeletal muscle preparation. Voltage-gated sodium channels (NaV) are a particularly attractive pharmacological target as they are involved in almost all physiological processes including action potential generation and conduction. Accordingly, venom peptides that interfere with NaV function provide a key defensive and predatory advantage to a range of invertebrate venomous species including cone snails, scorpions, spiders, and anemones. Enhanced activation or delayed inactivation of sodium channels by toxins is associated with the extremely rapid onset of tetanic/excitatory paralysis in envenomed prey animals. A strong selection pressure exists for the evolution of such toxins where there is a high chance of prey escape. However, despite their prevalence in other venomous species, toxins causing delay of sodium channel inhibition have never previously been described in vertebrate venoms. Here we show that NaV modulators, convergent with those of invertebrates, have evolved in the venom of the long-glanded coral snake. Calliotoxin represents a functionally novel class of 3FTx and a structurally novel class of NaV toxins that will provide significant insights into the pharmacology and physiology of NaV. The toxin represents a remarkable case of functional convergence between invertebrate and vertebrate venom systems in response to similar selection pressures. These results underscore the dynamic evolution of the Toxicofera reptile system and reinforces the value of using evolution as a roadmap for biodiscovery.

A Tricky Trait: Applying the Fruits of the "Function Debate" in the Philosophy of Biology to the "Venom Debate" in the Science of Toxinology

The "function debate" in the philosophy of biology and the "venom debate" in the science of toxinology are conceptually related. Venom systems are complex multifunctional traits that have evolved independently numerous times throughout the animal kingdom. No single concept of function, amongst those popularly defended, appears adequate to describe these systems in all their evolutionary contexts and extant variations. As such, a pluralistic view of function, previously defended by some philosophers of biology, is most appropriate. Venom systems, like many other functional traits, exist in nature as points on a continuum and the boundaries between "venomous" and "non-venomous" species may not always be clearly defined. This paper includes a brief overview of the concept of function, followed by in-depth discussion of its application to venom systems. A sound understanding of function may aid in moving the venom debate forward. Similarly, consideration of a complex functional trait such as venom may be of interest to philosophers of biology.

Canopy Venom: Proteomic Comparison among New World Arboreal Pit-Viper Venoms

Central and South American pitvipers, belonging to the genera Bothrops and Bothriechis, have independently evolved arboreal tendencies. Little is known regarding the composition and activity of their venoms. In order to close this knowledge gap, venom proteomics and toxin activity of species of Bothriechis, and Bothrops (including Bothriopsis) were investigated through established analytical methods. A combination of proteomics and bioactivity techniques was used to demonstrate a similar diversification of venom composition between large and small species within Bothriechis and Bothriopsis. Increasing our understanding of the evolution of complex venom cocktails may facilitate future biodiscoveries.

Firing the sting: chemically induced discharge of cnidae reveals novel proteins and peptides from box jellyfish (Chironex fleckeri) venom

Cnidarian venom research has lagged behind other toxinological fields due to technical difficulties in recovery of the complex venom from the microscopic nematocysts. Here we report a newly developed rapid, repeatable and cost effective technique of venom preparation, using ethanol to induce nematocyst discharge and to recover venom contents in one step. Our model species was the Australian box jellyfish (Chironex fleckeri), which has a notable impact on public health. By utilizing scanning electron microscopy and light microscopy, we examined nematocyst external morphology before and after ethanol treatment and verified nematocyst discharge. Further, to investigate nematocyst content or "venom" recovery, we utilized both top-down and bottom-up transcriptomics-proteomics approaches and compared the proteome profile of this new ethanol recovery based method to a previously reported high activity and recovery protocol, based upon density purified intact cnidae and pressure induced disruption. In addition to recovering previously characterized box jellyfish toxins, including CfTX-A/B and CfTX-1, we recovered putative metalloproteases and novel expression of a small serine protease inhibitor. This study not only reveals a much more complex toxin profile of Australian box jellyfish venom but also suggests that ethanol extraction method could augment future cnidarian venom proteomics research efforts.

Fossilized venom: the unusually conserved venom profiles of Heloderma species (beaded lizards and gila monsters)

Research into snake venoms has revealed extensive variation at all taxonomic levels. Lizard venoms, however, have received scant research attention in general, and no studies of intraclade variation in lizard venom composition have been attempted to date. Despite their iconic status and proven usefulness in drug design and discovery, highly venomous helodermatid lizards (gila monsters and beaded lizards) have remained neglected by toxinological research. Proteomic comparisons of venoms of three helodermatid lizards in this study has unravelled an unusual similarity in venom-composition, despite the long evolutionary time (~30 million years) separating H. suspectum from the other two species included in this study (H. exasperatum and H. horridum). Moreover, several genes encoding the major helodermatid toxins appeared to be extremely well-conserved under the influence of negative selection (but with these results regarded as preliminary due to the scarcity of available sequences). While the feeding ecologies of all species of helodermatid lizard are broadly similar, there are significant morphological differences between species, which impact upon relative niche occupation.

Extreme venom variation in Middle Eastern vipers: a proteomics comparison of Eristicophis macmahonii, Pseudocerastes fieldi and Pseudocerastes persicus

Venoms of the viperid sister genera Eristicophis and Pseudocerastes are poorly studied despite their anecdotal reputation for producing severe or even lethal envenomations. This is due in part to the remote and politically unstable regions that they occupy. All species contained are sit and wait ambush feeders. Thus, this study examined their venoms through proteomics techniques in order to establish if this feeding ecology, and putatively low levels of gene flow, have resulted in significant variations in venom profile. The techniques indeed revealed extreme venom variation. This has immediate implications as only one antivenom is made (using the venom of Pseudocerastes persicus) yet the proteomic variation suggests that it would be of only limited use for the other species, even the sister species Pseudocerastes fieldi. The high degree of variation however also points toward these species being rich resources for novel compounds which may have use as lead molecules in drug design and development.

BIOLOGICAL SIGNIFICANCE

These results show extreme venom variation between these closely related snakes. These results have direct implications for the treatment of the envenomed patient.

A ray of venom: Combined proteomic and transcriptomic investigation of fish venom composition using barb tissue from the blue-spotted stingray (Neotrygon kuhlii)

Fish venoms remain almost completely unstudied despite the large number of species. In part this is due to the inherent nature of fish venoms, in that they are highly sensitive to heat, pH, lyophilisation, storage and repeated freeze-thawing. They are also heavily contaminated with mucus, which makes proteomic study difficult. Here we describe a novel protein-handling protocol to remove mucus contamination, utilising ammonium sulphate and acetone precipitation. We validated this approach using barb venom gland tissue protein extract from the blue-spotted stingray Neotrygon kuhlii. We analysed the protein extract using 1D and 2D gels with LC-MS/MS sequencing. Protein annotation was underpinned by a venom gland transcriptome. The composition of our N. kuhlii venom sample revealed a variety of protein types that are completely novel to animal venom systems. Notably, none of the detected proteins exhibited similarity to the few toxin components previously characterised from fish venoms, including those found in other stingrays. Putative venom toxins identified here included cystatin, peroxiredoxin and galectin. Our study represents the first combined survey of gene and protein composition from the venom apparatus of any fish and our novel protein handling method will aid the future characterisation of toxins from other unstudied venomous fish lineages.

BIOLOGICAL SIGNIFICANCE

These results show an efficient manner for removing mucus from fish venoms. These results are the first insights into the evolution of proteins present on stingrayvenom barbs.

Vintage venoms: proteomic and pharmacological stability of snake venoms stored for up to eight decades

For over a century, venom samples from wild snakes have been collected and stored around the world. However, the quality of storage conditions for "vintage" venoms has rarely been assessed. The goal of this study was to determine whether such historical venom samples are still biochemically and pharmacologically viable for research purposes, or if new sample efforts are needed. In total, 52 samples spanning 5 genera and 13 species with regional variants of some species (e.g., 14 different populations of Notechis scutatus) were analysed by a combined proteomic and pharmacological approach to determine protein structural stability and bioactivity. When venoms were not exposed to air during storage, the proteomic results were virtually indistinguishable from that of fresh venom and bioactivity was equivalent or only slightly reduced. By contrast, a sample of Acanthophis antarcticus venom that was exposed to air (due to a loss of integrity of the rubber stopper) suffered significant degradation as evidenced by the proteomics profile. Interestingly, the neurotoxicity of this sample was nearly the same as fresh venom, indicating that degradation may have occurred in the free N- or C-terminus chains of the proteins, rather than at the tips of loops where the functional residues are located. These results suggest that these and other vintage venom collections may be of continuing value in toxin research. This is particularly important as many snake species worldwide are declining due to habitat destruction or modification. For some venoms (such as N. scutatus from Babel Island, Flinders Island, King Island and St. Francis Island) these were the first analyses ever conducted and these vintage samples may represent the only venom ever collected from these unique island forms of tiger snakes. Such vintage venoms may therefore represent the last remaining stocks of some local populations and thus are precious resources. These venoms also have significant historical value as the Oxyuranus venoms analysed include samples from the first coastal taipan (Oxyuranus scutellatus) collected for antivenom production (the snake that killed the collector Kevin Budden), as well as samples from the first Oxyuranus microlepidotus specimen collected after the species' rediscovery in 1976. These results demonstrate that with proper storage techniques, venom samples can retain structural and pharmacological stability. This article is part of a Special Issue entitled: Proteomics of non-model organisms.

BIOLOGICAL SIGNIFICANCE

Molecular evolution of vertebrate neurotrophins: co-option of the highly conserved nerve growth factor gene into the advanced snake venom arsenalf

Neurotrophins are a diverse class of structurally related proteins, essential for neuronal development, survival, plasticity and regeneration. They are characterized by major family members, such as the nerve growth factors (NGF), brain-derived neurotrophic factors (BDNF) and neurotrophin-3 (NT-3), which have been demonstrated here to lack coding sequence variations and follow the regime of negative selection, highlighting their extremely important conserved role in vertebrate homeostasis. However, in stark contrast, venom NGF secreted as part of the chemical arsenal of the venomous advanced snake family Elapidae (and to a lesser extent Viperidae) have characteristics consistent with the typical accelerated molecular evolution of venom components. This includes a rapid rate of diversification under the significant influence of positive-selection, with the majority of positively-selected sites found in the secreted β-polypeptide chain (74%) and on the molecular surface of the protein (92%), while the core structural and functional residues remain highly constrained. Such focal mutagenesis generates active residues on the toxin molecular surface, which are capable of interacting with novel biological targets in prey to induce a myriad of pharmacological effects. We propose that caenophidian NGFs could participate in prey-envenoming by causing a massive release of chemical mediators from mast cells to mount inflammatory reactions and increase vascular permeability, thereby aiding the spread of other toxins and/or by acting as proapoptotic factors. Despite their presence in reptilian venom having been known for over 60 years, this is the first evidence that venom-secreted NGF follows the molecular evolutionary pattern of other venom components, and thus likely participates in prey-envenomation.

Functional characterization on invertebrate and vertebrate tissues of tachykinin peptides from octopus venoms

It has been previously shown that octopus venoms contain novel tachykinin peptides that despite being isolated from an invertebrate, contain the motifs characteristic of vertebrate tachykinin peptides rather than being more like conventional invertebrate tachykinin peptides. Therefore, in this study we examined the effect of three variants of octopus venom tachykinin peptides on invertebrate and vertebrate tissues. While there were differential potencies between the three peptides, their relative effects were uniquely consistent between invertebrate and vertebrae tissue assays. The most potent form (OCT-TK-III) was not only the most anionically charged but also was the most structurally stable. These results not only reveal that the interaction of tachykinin peptides is more complex than previous structure-function theories envisioned, but also reinforce the fundamental premise that animal venoms are rich resources of novel bioactive molecules, which are useful investigational ligands and some of which may be useful as lead compounds for drug design and development.

Squeezers and leaf-cutters: differential diversification and degeneration of the venom system in toxicoferan reptiles

Although it has been established that all toxicoferan squamates share a common venomous ancestor, it has remained unclear whether the maxillary and mandibular venom glands are evolving on separate gene expression trajectories or if they remain under shared genetic control. We show that identical transcripts are simultaneously expressed not only in the mandibular and maxillary glands, but also in the enigmatic snake rictal gland. Toxin molecular frameworks recovered in this study were three-finger toxin (3FTx), CRiSP, crotamine (beta-defensin), cobra venom factor, cystatin, epididymal secretory protein, kunitz, L-amino acid oxidase, lectin, renin aspartate protease, veficolin, and vespryn. We also discovered a novel low-molecular weight disulfide bridged peptide class in pythonid snake glands. In the iguanian lizards, the most highly expressed are potentially antimicrobial in nature (crotamine (beta-defensin) and cystatin), with crotamine (beta-defensin) also the most diverse. However, a number of proteins characterized from anguimorph lizards and caenophidian snakes with hemotoxic or neurotoxic activities were recruited in the common toxicoferan ancestor and remain expressed, albeit in low levels, even in the iguanian lizards. In contrast, the henophidian snakes express 3FTx and lectin toxins as the dominant transcripts. Even in the constricting pythonid and boid snakes, where the glands are predominantly mucous-secreting, low-levels of toxin transcripts can be detected. Venom thus appears to play little role in feeding behavior of most iguanian lizards or the powerful constricting snakes, and the low levels of expression argue against a defensive role. However, clearly the incipient or secondarily atrophied venom systems of these taxa may be a source of novel compounds useful in drug design and discovery.

Dracula's children: molecular evolution of vampire bat venom

While vampire bat oral secretions have been the subject of intense research, efforts have concentrated only on two components: DSPA (Desmodus rotundus salivary plasminogen activator) and Draculin. The molecular evolutionary history of DSPA has been elucidated, while conversely draculin has long been known from only a very small fragment and thus even the basic protein class was not even established. Despite the fact that vampire bat venom has a multitude of effects unaccounted by the documented bioactivities of DSPA and draculin, efforts have not been made to establish what other bioactive proteins are secreted by their submaxillary gland. In addition, it has remained unclear whether the anatomically distinct anterior and posterior lobes of the submaxillary gland are evolving on separate gene expression trajectories or if they remain under the shared genetic control. Using a combined proteomic and transcriptomic approach, we show that identical proteins are simultaneously expressed in both lobes. In addition to recovering the known structural classes of DSPA, we recovered a novel DSPA isoform as well as obtained a very large sequence stretch of draculin and thus established that it is a mutated version of the lactotransferrin scaffold. This study reveals a much more complex secretion profile than previously recognised. In addition to obtaining novel versions of scaffolds convergently recruited into other venoms (allergen-like, CRiSP, kallikrein, Kunitz, lysozyme), we also documented novel expression of small peptides related to calcitonin, PACAP, and statherin. Other overexpressed protein types included BPI-fold, lacritin, and secretoglobin. Further, we investigate the molecular evolution of various vampire bat venom-components and highlight the dominant role of positive selection in the evolution of these proteins. Conspicuously many of the proteins identified in the proteome were found to be homologous to proteins with known activities affecting vasodilation and platelet aggregation. We show that vampire bat venom proteins possibly evade host immune response by the mutation of the surface chemistry through focal mutagenesis under the guidance of positive Darwinian selection. These results not only contribute to the body of knowledge regarding haematophagous venoms but also provide a rich resource for novel lead compounds for use in drug design and development.

BIOLOGICAL SIGNIFICANCE

These results have direct implications in understanding the molecular evolutionary history of vampire bat venom. The unusual peptides discovered reinforce the value of studying such neglected taxon for biodiscovery.

Venom proteomic characterization and relative antivenom neutralization of two medically important Pakistani elapid snakes (Bungarus sindanus and Naja naja)

Intra- and interspecific variation in venom composition has been shown to have a major effect upon the efficacy of antivenoms. Due to the absence of domestically produced antivenoms, Pakistan is wholly reliant upon antivenoms produced in other countries, such as India. However, the efficacy of these antivenoms in neutralising the venoms of Pakistani snakes has not been ascertained. This is symptomatic of the general state of toxicological research in this country, which has a myriad of highly toxic and medically important venomous animals. Thus, there is a dire need for knowledge regarding the fundamental proteomics of these venoms and applied knowledge of the relative efficacy of foreign antivenoms. Here we present the results of our proteomic research on two medically important snakes of Pakistan: Bungarus sindanus and Naja naja. Indian Polyvalent Antivenom (Bharat Serums and Vaccines Ltd), which is currently marketed for use in Pakistan, was completely ineffective against either Pakistani species. In addition to the expected pre- and post-synaptic neurotoxic activity, the venom of the Pakistan population of N. naja was shown to be quite divergent from other populations of this species in being potently myotoxic. These results highlight the importance of studying divergent species and isolated populations, where the same data not only elucidates clinical problems in need of immediate attention, but also uncovers sources for novel toxins with potentially useful activities.

BIOLOGICAL SIGNIFICANCE

Pakistan Bungarus sindanus and Naja naja venoms are differentially complex. Naja naja is potently myotoxic. Neither venom is neutralized by Indian antivenom. These results have direct implications for the treatment of envenomed patients in Pakistan. The unusually myotoxic effects of Naja naja demonstrates the value of studying remote populations for biodiscovery.

Molecular phylogeny and evolution of the proteins encoded by coleoid (cuttlefish, octopus, and squid) posterior venom glands

In this study, we report for the first time a detailed evaluation of the phylogenetic history and molecular evolution of the major coleoid toxins: CAP, carboxypeptidase, chitinase, metalloprotease GON-domain, hyaluronidase, pacifastin, PLA2, SE-cephalotoxin and serine proteases, with the carboxypeptidase and GON-domain documented for the first time in the coleoid venom arsenal. We show that although a majority of sites in these coleoid venom-encoding genes have evolved under the regime of negative selection, a very small proportion of sites are influenced by the transient selection pressures. Moreover, nearly 70 % of these episodically adapted sites are confined to the molecular surface, highlighting the importance of variation of the toxin surface chemistry. Coleoid venoms were revealed to be as complex as other venoms that have traditionally been the recipient of the bulk of research efforts. The presence of multiple peptide/protein types in coleoids similar to those present in other animal venoms identifies a convergent strategy, revealing new information as to what characteristics make a peptide/protein type amenable for recruitment into chemical arsenals. Coleoid venoms have significant potential not only for understanding fundamental aspects of venom evolution but also as an untapped source of novel toxins for use in drug design and discovery.

Differential evolution and neofunctionalization of snake venom metalloprotease domains

Snake venom metalloproteases (SVMP) are composed of five domains: signal peptide, propeptide, metalloprotease, disintegrin, and cysteine-rich. Secreted toxins are typically combinatorial variations of the latter three domains. The SVMP-encoding genes of Psammophis mossambicus venom are unique in containing only the signal and propeptide domains. We show that the Psammophis SVMP propeptide evolves rapidly and is subject to a high degree of positive selection. Unlike Psammophis, some species of Echis express both the typical multidomain and the unusual monodomain (propeptide only) SVMP, with the result that a lower level of variation is exerted upon the latter. We showed that most mutations in the multidomain Echis SVMP occurred in the protease domain responsible for proteolytic and hemorrhagic activities. The cysteine-rich and disintegrin-like domains, which are putatively responsible for making the P-III SVMPs more potent than the P-I and P-II forms, accumulate the remaining variation. Thus, the binding sites on the molecule's surface are evolving rapidly whereas the core remains relatively conserved. Bioassays conducted on two post-translationally cleaved novel proline-rich peptides from the P. mossambicus propeptide domain showed them to have been neofunctionalized for specific inhibition of mammalian a7 neuronal nicotinic acetylcholine receptors. We show that the proline rich postsynaptic specific neurotoxic peptides from Azemiops feae are the result of convergent evolution within the precursor region of the C-type natriuretic peptide instead of the SVMP. The results of this study reinforce the value of studying obscure venoms for biodiscovery of novel investigational ligands.

The structural and functional diversification of the Toxicofera reptile venom system

The evolutionary origin and diversification of the reptilian venom system is described. The resolution of higher-order molecular phylogenetics has clearly established that a venom system is ancestral to snakes. The diversification of the venom system within lizards is discussed, as is the role of venom delivery in the behavioural ecology of these taxa (particularly Varanus komodoensis). The more extensive diversification of the venom system in snakes is summarised, including its loss in some clades. Finally, we discuss the contentious issue of a definition for "venom", supporting an evolutionary definition that recognises the homology of both the venom delivery systems and the toxins themselves.