Venom down under: dynamic evolution of Australian elapid snake toxins

VenomCap: An exon-capture probe set for the targeted sequencing of snake venom genes

Snake venoms are complex mixtures of toxic proteins that hold significant medical, pharmacological and evolutionary interest. To better understand the genetic diversity underlying snake venoms, we developed VenomCap, a novel exon-capture probe set targeting toxin-coding genes from a wide range of elapid snakes, with a particular focus on the ecologically diverse and medically important subfamily Hydrophiinae. We tested the capture success of VenomCap across 24 species, representing all major elapid lineages. We included snake phylogenomic probes in the VenomCap capture set, allowing us to compare capture performance between venom and phylogenomic loci and to infer elapid phylogenetic relationships. We demonstrated VenomCap's ability to recover exons from ~1500 target markers, representing a total of 24 known venom gene families, which includes the dominant gene families found in elapid venoms. We find that VenomCap's capture results are robust across all elapids sampled, and especially among hydrophiines, with respect to measures of target capture success (target loci matched, sensitivity, specificity and missing data). As a cost-effective and efficient alternative to full genome sequencing, VenomCap can dramatically accelerate the sequencing and analysis of venom gene families. Overall, our tool offers a model for genomic studies on snake venom gene diversity and evolution that can be expanded for comprehensive comparisons across the other families of venomous snakes.

Unveiling Novel Kunitz- and Waprin-Type Toxins in the Micrurus mipartitus Coral Snake Venom Gland: An In Silico Transcriptome Analysis

Kunitz-type peptide expression has been described in the venom of snakes of the Viperidae, Elapidae and Colubridae families. This work aimed to identify these peptides in the venom gland transcriptome of the coral snake Micrurus mipartitus. Transcriptomic analysis revealed a high diversity of venom-associated Kunitz serine protease inhibitor proteins (KSPIs). A total of eight copies of KSPIs were predicted and grouped into four distinctive types, including short KSPI, long KSPI, Kunitz-Waprin (Ku-WAP) proteins, and a multi-domain Kunitz-type protein. From these, one short KSPI showed high identity with Micrurus tener and Austrelaps superbus. The long KSPI group exhibited similarity within the Micrurus genus and showed homology with various elapid snakes and even with the colubrid Pantherophis guttatus. A third group suggested the presence of Kunitz domains in addition to a whey-acidic-protein-type four-disulfide core domain. Finally, the fourth group corresponded to a transcript copy with a putative 511 amino acid protein, formerly annotated as KSPI, which UniProt classified as SPINT1. In conclusion, this study showed the diversity of Kunitz-type proteins expressed in the venom gland transcriptome of M. mipartitus.

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.

The Eastern Bandy Bandy Vermicella annulata, expresses high abundance of SVMP, CRiSP and Kunitz protein families in its venom proteome

The Australian elapid snake radiation (Hydrophiinae) has evolved in the absence of competition from other advanced snakes. This has resulted in ecological specialisation in Australian elapids and the potential for venom proteomes divergent to other elapids. We characterised the venom of the Australian elapid Vermicella annulata (eastern bandy bandy). The venom was analysed using a two-dimensional fractionation process consisting of reverse-phase high-performance liquid chromatography then sodium dodecyl sulphate polyacrylamide gel electrophoresis, followed by bottom-up proteomics. Resulting peptides were matched to a species-specific transcriptome and 87% of the venom was characterised. We identified 11 toxins in the venom from six families: snake venom metalloproteinases (SVMP; 24.2%; two toxins) that are class P-III SVMPs containing a disintegrin-like domain, three-finger toxins (3FTx; 21.6%; five toxins), kunitz peptides (KUN; 19.5%; one toxin), cysteine-rich secretory proteins (CRiSP; 18%; one toxin), and phospholipase A2 (PLA2; 4%; two toxins). The venom had low toxin diversity with five protein families having one or two toxins, except for 3FTx with five different toxins. V. annulata expresses an unusual venom proteome, with high abundances of CRiSP, KUN and SVMP, which are not normally highly expressed in elapid venoms. This unusual venom composition could be an adaptation to its specialised diet. BIOLOGICAL SIGNIFICANCE: Although the Australian elapid radiation represents the most extensive speciation event of elapids on any continent, with 100 terrestrial species, the venom composition of these snakes has rarely been investigated, with only five species currently characterised. Here we provide the venom proteome of a sixth species, Vermicella annulata. The venom of this species could be particularly informative from an evolutionary perspective, as it is an extreme dietary specialist, only preying on blind snakes (Typhlopidae). We show that V. annulata expresses a highly unusual venom for an elapid, due to the high abundance of the protein families SVMP, CRiSP, and KUN, which together make up 61% of the venom. When averaged across all species, a typical elapid venom is 82% PLA2 and 3FTx. This is the second recorded instance of an Australian elapid having evolved highly divergent venom expression.

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

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

Identification of waprin and its microbicidal activity: A novel protein component of honeybee (Apis mellifera) venom

Bee venom is a rich source of biologically and pharmacologically active proteins. Waprin is a protein component of venoms; however, waprin has yet to be identified in bee venom. Moreover, the biological functions of waprin in venoms remain poorly characterized. Thus, in this study, we have identified and characterized waprin: a novel protein component from the venom of honeybees (Apis mellifera). The waprin in A. mellifera venom (Amwaprin) was found to consist of an 80-amino acid mature peptide, in which the whey acidic protein domain contains four conserved disulfide bonds. We discovered the presence of the Amwaprin protein in secreted venom by using an antibody against recombinant Amwaprin produced in baculovirus-infected insect cells. Recombinant Amwaprin exhibited inhibitory activity against microbial serine proteases and elastases but not thrombin or plasmin. It recognized carbohydrates in the microbial cell wall molecules and bound to the live microbial surfaces. The binding action of Amwaprin produced its microbicidal activity by inducing structural damage to bacterial and fungal cell walls. In addition, recombinant Amwaprin is heat-stable and contains no hemolytic activity. These findings demonstrate that Amwaprin acts as a microbicidal and anti-elastolytic agent.

Taipan Natriuretic Peptides Are Potent and Selective Agonists for the Natriuretic Peptide Receptor A

Cardiovascular ailments are a major cause of mortality where over 1.3 billion people suffer from hypertension leading to heart-disease related deaths. Snake venoms possess a broad repertoire of natriuretic peptides with therapeutic potential for treating hypertension, congestive heart failure, and related cardiovascular disease. We now describe several taipan (Oxyuranus microlepidotus) natriuretic peptides TNPa-e which stimulated cGMP production through the natriuretic peptide receptor A (NPR-A) with higher potencies for the rat NPR-A (rNPR-A) over human NPR-A (hNPR-A). TNPc and TNPd were the most potent, demonstrating 100- and 560-fold selectivity for rNPR-A over hNPR-A. In vivo studies found that TNPc decreased diastolic and systolic blood pressure (BP) and increased heart rate (HR) in conscious normotensive rabbits, to a level that was similar to that of human atrial natriuretic peptide (hANP). TNPc also enhanced the bradycardia due to cardiac afferent stimulation (Bezold-Jarisch reflex). This indicated that TNPc possesses the ability to lower blood pressure and facilitate cardiac vagal afferent reflexes but unlike hANP does not produce tachycardia. The 3-dimensional structure of TNPc was well defined within the pharmacophoric disulfide ring, displaying two turn-like regions (RMSD = 1.15 Å). Further, its much greater biological stability together with its selectivity and potency will enhance its usefulness as a biological tool.

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.

The Target Selects the Toxin: Specific Amino Acids in Snake-Prey Nicotinic Acetylcholine Receptors That Are Selectively Bound by King Cobra Venoms

Snake venom is an adaptive ecological trait that has evolved primarily as a form of prey subjugation. Thus, the selection pressure for toxin diversification is exerted by the prey’s physiological targets, with this pressure being particularly acute for specialist feeders, such as the King Cobra species, all of which are snake-prey specialists. However, while extensive research has been undertaken to elucidate key amino acids that guide toxin structure–activity relationships, reciprocal investigations into the specific sites guiding prey-lineage selective effects have been lacking. This has largely been due to the lack of assay systems amenable to systematic amino acid replacements of targeted proteins in the prey’s physiological pathways. To fill this knowledge gap, we used a recently described approach based upon mimotope peptides corresponding to the orthosteric site of nicotinic acetylcholine receptor alpha-1 subunits, a major binding site for snake venom neurotoxins that cause flaccid paralysis. We investigated the venoms of four different types of King Cobra (Cambodian, Javan, Malaysian, and Thai). This approach allowed for the determination of the key amino acid positions in King Cobra snake prey that are selectively bound by the toxins, whereby replacing these amino acids in the snake-prey orthosteric site with those from lizards or rats resulted in a significantly lower level of binding by the venoms, while conversely replacing the lizard or rat amino acids with those from the snake at that position increased the binding. By doing such, we identified three negatively charged amino acids in the snake orthosteric site that are strongly bound by the positively charged neurotoxic three-finger toxins found in King Cobra venom. This study, thus, sheds light on the selection pressures exerted by a specialist prey item for the evolution of lineage-selective toxins.

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 A2 (PLA2s) in venom-induced effects was assessed with the sPLA2 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 sPLA2 inhibitor varespladib, we conclude that PLA2 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.

A retrospective evaluation of eastern coral snake envenomation and antivenom administration in cats: 30 cases (2012-2019)

Thirty cats were identified to be have been suspected to have a potential coral snake envenomation after searching medical records from 2012 to 2019 at a university teaching hospital. The records were reviewed and evaluated for signalment, date and time of the snake encounter, elapsed time between encounter and hospital examination, presenting complaint, initial physical examination findings, initial laboratory findings, antivenom dose and duration of administration, adverse reactions to antivenom, additional treatments administered, progression of clinical signs, length of hospitalization, and outcome. Thirteen cats presented with clinical signs consistent with envenomation while 17 cats were treated for possible asymptomatic envenomation, as defined by the owner discovering a live or dead coral snake in their home or on their property. Initial physical examination findings included tachypnea with short shallow breaths and use of accessory muscles; tetraparesis with normal or decreased to absent spinal reflexes; cranial nerve deficits including decreased to absent gag, slow pupillary light reflexes, and absent physiologic nystagmus; and normal or altered mentation. Laboratory findings included hypercapnia, hyperglycemia, hypokalemia, increased aspartate aminotransferase activity, increased alanine aminotransferase activity, echinocytosis, leukocytosis, azotemia, and hyperlactatemia. Twenty-eight cats received antivenom; two cats received two vials while twenty-six cats received one vial. Antivenom reaction was suspected in one cat that developed facial swelling during administration of the drug. Average length of hospitalization was 1 day for cats without clinical signs and 3 days for cats with clinical signs. Twenty-nine cats survived to discharge. Due to the inclusion criteria of the study, cats euthanized on presentation or discharged without receiving antivenom may have been unintentionally excluded from the study. Diagnosis of eastern coral snake envenomation should be suspected in the cat that has acute onset of lower motor neuron neuropathy. Prognosis with treatment is considered good with 97% of cats surviving to discharge. Antivenom reaction occurred in 3.5% of administrations with none being fatal. Monitoring of hypercapnia was critical in making the decision to mechanically ventilate patients. Supportive care that includes antivenom administration, recumbency care, and mechanical ventilation if needed are the mainstays of therapy.

Sweet as honey, bitter as bile: Mitochondriotoxic peptides and other therapeutic proteins isolated from animal tissues, for dealing with mitochondrial apoptosis

Mitochondria are subcellular organelles involved in cell metabolism and cell life-cycle. Their role in apoptosis regulation makes them an interesting target of new drugs for dealing with cancer or rare diseases. Several peptides and proteins isolated from animal and plant sources are known for their therapeutic properties and have been tested on cancer cell-lines and xenograft murine models, highlighting their ability in inducing cell-death by triggering mitochondrial apoptosis. Some of those molecules have been even approved as drugs. Conversely, many other bioactive compounds are still under investigation for their proapoptotic properties. In this review we report about a group of peptides, isolated from animal venoms, with potential therapeutic properties related to their ability in triggering mitochondrial apoptosis. This class of compounds is known with different names, such as mitochondriotoxins or mitocans.

Snake three-finger α-neurotoxins and nicotinic acetylcholine receptors: molecules, mechanisms and medicine

Snake venom three-finger α-neurotoxins (α-3FNTx) act on postsynaptic nicotinic acetylcholine receptors (nAChRs) at the neuromuscular junction (NMJ) to produce skeletal muscle paralysis. The discovery of the archetypal α-bungarotoxin (α-BgTx), almost six decades ago, exponentially expanded our knowledge of membrane receptors and ion channels. This included the localisation, isolation and characterization of the first receptor (nAChR); and by extension, the pathophysiology and pharmacology of neuromuscular transmission and associated pathologies such as myasthenia gravis, as well as our understanding of the role of α-3FNTxs in snakebite envenomation leading to novel concepts of targeted treatment. Subsequent studies on a variety of animal venoms have yielded a plethora of novel toxins that have revolutionized molecular biomedicine and advanced drug discovery from bench to bedside. This review provides an overview of nAChRs and their subtypes, classification of α-3FNTxs and the challenges of typifying an increasing arsenal of structurally and functionally unique toxins, and the three-finger protein (3FP) fold in the context of the uPAR/Ly6/CD59/snake toxin superfamily. The pharmacology of snake α-3FNTxs including their mechanisms of neuromuscular blockade, variations in reversibility of nAChR interactions, specificity for nAChR subtypes or for distinct ligand-binding interfaces within a subtype and the role of α-3FNTxs in neurotoxic envenomation are also detailed. Lastly, a reconciliation of structure-function relationships between α-3FNTx and nAChRs, derived from historical mutational and biochemical studies and emerging atomic level structures of nAChR models in complex with α-3FNTxs is discussed.

Cloning and identification of a new multifunctional Ascaris-type peptide from the hemolymph of Buthus martensii Karsch

Only a few work have been done for peptides from non-venom gland tissues of venomous animals. Here, with the help of the whole body transcriptomic and the hemolymph proteomic data of the Chinese scorpion Buthus martensii Karsch, we identified the first Ascaris-type peptide BmHDP from scorpion hemolymph. The precursor of BmHDP has 80 residues, including a 16 residue signal peptide and a 64 residue mature peptide. The mature peptide has 10 conserved cysteines and adopts a conserved Ascaris-type fold. Using combined inclusion body refolding and biochemical identification strategies, recombinant BmHDP was obtained successfully. Protease inhibitory assays showed that BmHDP inhibited chymotrypsin apparently at a concentration of 8 nM. Patch-clamp experiments showed that BmHDP inhibited the Kv1.3 potassium channel apparently at a concentration of 1000 nM. Coagulation experiment assays showed that BmHDP inhibited intrinsic coagulation pathway apparently at a concentration of 500 nM. To the best of our knowledge, BmHDP is the first Ascaris-type peptide from scorpion hemolymph. Our work highlighted a functional link between scorpion non-venom gland peptides and venom gland toxin peptides, and suggested that scorpion hemolymph might be a new source of bioactive peptides.

Evolutionary Interpretations of Nicotinic Acetylcholine Receptor Targeting Venom Effects by a Clade of Asian Viperidae Snakes

Ecological variability among closely related species provides an opportunity for evolutionary comparative studies. Therefore, to investigate the origin and evolution of neurotoxicity in Asian viperid snakes, we tested the venoms of Azemiops feae, Calloselasma rhodostoma, Deinagkistrodon acutus, Tropidolaeums subannulatus, and T. wagleri for their relative specificity and potency upon the amphibian, lizard, bird, rodent, and human α-1 (neuromuscular) nicotinic acetylcholine receptors. We utilised a biolayer interferometry assay to test the binding affinity of these pit viper venoms to orthosteric mimotopes of nicotinic acetylcholine receptors binding region from a diversity of potential prey types. The Tropidolaemus venoms were much more potent than the other species tested, which is consistent with the greater prey escape potential in arboreal niches. Intriguingly, the venom of C. rhodostoma showed neurotoxic binding to the α-1 mimotopes, a feature not known previously for this species. The lack of prior knowledge of neurotoxicity in this species is consistent with our results due to the bias in rodent studies and human bite reports, whilst this venom had a greater binding affinity toward amphibian and diapsid α-1 targets. The other large terrestrial species, D. acutus, did not display any meaningful levels of neurotoxicity. These results demonstrate that whilst small peptide neurotoxins are a basal trait of these snakes, it has been independently amplified on two separate occasions, once in Azemiops and again in Tropidolaemus, and with Calloselasma representing a third possible amplification of this trait. These results also point to broader sources of novel neuroactive peptides with the potential for use as lead compounds in drug design and discovery.

An Appetite for Destruction: Detecting Prey-Selective Binding of α-Neurotoxins in the Venom of Afro-Asian Elapids

Prey-selective venoms and toxins have been documented across only a few species of snakes. The lack of research in this area has been due to the absence of suitably flexible testing platforms. In order to test more species for prey specificity of their venom, we used an innovative taxonomically flexible, high-throughput biolayer interferometry approach to ascertain the relative binding of 29 α-neurotoxic venoms from African and Asian elapid representatives (26 Naja spp., Aspidelaps scutatus, Elapsoidea boulengeri, and four locales of Ophiophagus hannah) to the alpha-1 nicotinic acetylcholine receptor orthosteric (active) site for amphibian, lizard, snake, bird, and rodent targets. Our results detected prey-selective, intraspecific, and geographical differences of α-neurotoxic binding. The results also suggest that crude venom that shows prey selectivity is likely driven by the proportions of prey-specific α-neurotoxins with differential selectivity within the crude venom. Our results also suggest that since the α-neurotoxic prey targeting does not always account for the full dietary breadth of a species, other toxin classes with a different pathophysiological function likely play an equally important role in prey immobilisation of the crude venom depending on the prey type envenomated. The use of this innovative and taxonomically flexible diverse assay in functional venom testing can be key in attempting to understanding the evolution and ecology of α-neurotoxic snake venoms, as well as opening up biochemical and pharmacological avenues to explore other venom effects.

Fulditoxin, representing a new class of dimeric snake toxins, defines novel pharmacology at nicotinic ACh receptors

BACKGROUND AND PURPOSE

Animal toxins have contributed significantly to our understanding of the neurobiology of receptors and ion channels. We studied the venom of the coral snake Micrurus fulvius fulvius and identified and characterized the structure and pharmacology of a new homodimeric neurotoxin, fulditoxin, that exhibited novel pharmacology at nicotinic ACh receptors (nAChRs).

EXPERIMENTAL APPROACH

Fulditoxin was isolated by chromatography, chemically synthesized, its structure determined by X-ray crystallography, and its pharmacological actions on nAChRs characterized by organ bath assays and two-electrode voltage clamp electrophysiology.

KEY RESULTS

Fulditoxin's distinct 1.95-Å quaternary structure revealed two short-chain three-finger α-neurotoxins (α-3FNTxs) non-covalently bound by hydrophobic interactions and an ability to bind metal and form tetrameric complexes, not reported previously for three-finger proteins. Although fulditoxin lacked all conserved amino acids canonically important for inhibiting nAChRs, it produced postsynaptic neuromuscular blockade of chick muscle at nanomolar concentrations, comparable to the prototypical α-bungarotoxin. This neuromuscular blockade was completely reversible, which is unusual for snake α-3FNTxs. Fulditoxin, therefore, interacts with nAChRs by utilizing a different pharmacophore. Unlike short-chain α-3FNTxs that bind only to muscle nAChRs, fulditoxin utilizes dimerization to expand its pharmacological targets to include human neuronal α4β2, α7, and α3β2 nAChRs which it blocked with IC₅₀ values of 1.8, 7, and 12 μM respectively.

CONCLUSIONS AND IMPLICATIONS

Based on its distinct quaternary structure and unusual pharmacology, we named this new class of dimeric Micrurus neurotoxins represented by fulditoxin as Σ-neurotoxins, which offers greater insight into understanding the interactions between nAChRs and peptide antagonists.

Activity of two key toxin groups in Australian elapid venoms show a strong correlation to phylogeny but not to diet

Background

The relative influence of diet and phylogeny on snake venom activity is a poorly understood aspect of snake venom evolution. We measured the activity of two enzyme toxin groups – phospholipase A₂ (PLA₂), and L-amino acid oxidase (LAAO) – in the venom of 39 species of Australian elapids (40% of terrestrial species diversity) and used linear parsimony and BayesTraits to investigate any correlation between enzyme activity and phylogeny or diet.

Results

PLA₂ activity ranged from 0 to 481 nmol/min/mg of venom, and LAAO activity ranged from 0 to 351 nmol/min/mg. Phylogenetic comparative methods, implemented in BayesTraits showed that enzyme activity was strongly correlated with phylogeny, more so for LAAO activity. For example, LAAO activity was absent in both the Vermicella and Pseudonaja/Oxyuranus clade, supporting previously proposed relationships among these disparate taxa. There was no association between broad dietary categories and either enzyme activity. There was strong evidence for faster initial rates of change over evolutionary time for LAAO (delta parameter mean 0.2), but no such pattern in PLA₂ (delta parameter mean 0.64). There were some exceptions to the phylogenetic patterns of enzyme activity: different PLA₂ activity in the ecologically similar sister-species Denisonia devisi and D. maculata; large inter-specific differences in PLA₂ activity in Hoplocephalus and Austrelaps.

Conclusions

We have shown that phylogeny is a stronger influence on venom enzyme activity than diet for two of the four major enzyme families present in snake venoms. PLA₂ and LAAO activities had contrasting evolutionary dynamics with the higher delta value for PLA₂ Some species/individuals lacked activity in one protein family suggesting that the loss of single protein family may not incur a significant fitness cost.

Omics Technologies for Profiling Toxin Diversity and Evolution in Snake Venom: Impacts on the Discovery of Therapeutic and Diagnostic Agents

Snake venoms are primarily composed of proteins and peptides, and these toxins have developed high selectivity to their biological targets. This makes venoms interesting for exploration into protein evolution and structure-function relationships. A single venom protein superfamily can exhibit a variety of pharmacological effects; these variations in activity originate from differences in functional sites, domains, posttranslational modifications, and the formations of toxin complexes. In this review, we discuss examples of how the major venom protein superfamilies have diversified, as well as how newer technologies in the omics fields, such as genomics, transcriptomics, and proteomics, can be used to characterize both known and unknown toxins.Because toxins are bioactive molecules with a rich diversity of activities, they can be useful as therapeutic and diagnostic agents, and successful examples of toxin applications in these areas are also reviewed. With the current rapid pace of technology, snake venom research and its applications will only continue to expand.

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.

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.

Engineering varied serine protease inhibitors by converting P1 site of BF9, a weakly active Kunitz-type animal toxin

Although there were a lot of weakly active animal toxins in the venoms, their values and applications are still mysterious, such as BF9, which is a Kunitz-type toxin isolated from the venom of the elapid snake Bungarus fasciatus. Here, we used BF9 to be a molecular scaffold, and engineered eight BF9-derived peptides by changing P1 site Asn17 of BF9, such as BF9-N17Y and BF9-N17T designed from the polar subfamily, BF9-N17L and BF9-N17G designed from the Non-polar subfamily, BF9-N17D designed from acidic subfamily, and BF9-N17H, BF9-N17K and BF9-N17R designed from basic subfamily. Through enzyme inhibitor experiment assays, we found a potent and selective chymotrypsin inhibitor BF9-N17Y, a potent and selective coagulation factor XIa inhibitor BF9-N17H, and two highly potent coagulation factor XIa inhibitors BF9-N17K and BF9-N17. APTT and PT assays further showed that BF9-N17H, BF9-N17K and BF9-N17R were three novel anticoagulants with selectively intrinsic coagulation pathway inhibitory activity. Considering that natural weakly active animal toxins are also a huge peptide resource, our present work might open a new window about pharmacological applications of weakly active animal toxins, which might be good templates for potent and selective molecular probe and lead drug designs.

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.

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.

The medical threat of mamba envenoming in sub-Saharan Africa revealed by genus-wide analysis of venom composition, toxicity and antivenomics profiling of available antivenoms

Mambas (genus Dendroaspis) are among the most feared of the medically important elapid snakes found in sub-Saharan Africa, but many facets of their biology, including the diversity of venom composition, remain relatively understudied. Here, we present a reconstruction of mamba phylogeny, alongside genus-wide venom gland transcriptomic and high-resolution top-down venomic analyses. Whereas the green mambas, D. viridis, D. angusticeps, D. j. jamesoni and D. j. kaimosae, express 3FTx-predominant venoms, black mamba (D. polylepis) venom is dominated by dendrotoxins I and K. The divergent terrestrial ecology of D. polylepis compared to the arboreal niche occupied by all other mambas makes it plausible that this major difference in venom composition is due to dietary variation. The pattern of intrageneric venom variability across Dendroaspis represented a valuable opportunity to investigate, in a genus-wide context, the variant toxicity of the venom, and the degree of paraspecific cross-reactivity between antivenoms and mamba venoms. To this end, the immunological profiles of the five mamba venoms were assessed against a panel of commercial antivenoms generated for the sub-Saharan Africa market. This study provides a genus-wide overview of which available antivenoms may be more efficacious in neutralising human envenomings caused by mambas, irrespective of the species responsible. The information gathered in this study lays the foundations for rationalising the notably different potency and pharmacological profiles of Dendroaspis venoms at locus resolution. This understanding will allow selection and design of toxin immunogens with a view to generating a safer and more efficacious pan-specific antivenom against any mamba envenomation.

BIOLOGICAL SIGNIFICANCE

The mambas (genus Dendroaspis) comprise five especially notorious medically important venomous snakes endemic to sub-Saharan Africa. Their highly potent venoms comprise a high diversity of pharmacologically active peptides, including extremely rapid-acting neurotoxins. Previous studies on mamba venoms have focused on the biochemical and pharmacological characterisation of their most relevant toxins to rationalize the common neurological and neuromuscular symptoms of envenomings caused by these species, but there has been little work on overall venom composition or comparisons between them. Only very recently an overview of the composition of the venom of two Dendroaspis species, D. angusticeps and D. polylepis, has been unveiled through venomics approaches. Here we present the first genus-wide transcriptomic-proteomic analysis of mamba venom composition. The transcriptomic analyses described in this paper have contributed 29 (D. polylepis), 23 (D. angusticeps), 40 (D. viridis), 25 (D. j. jamesoni) and 21 (D. j. kaimosae), novel full-length toxin sequences to the non-redundant Dendroaspis sequence database. The mamba genus-wide venomic analysis demonstrated that major D. polylepis venom components are Kunitz-fold family toxins. This feature is unique in relation to the relatively conserved three-finger toxin (3FTx)-dominated venom compositions of the green mambas. Venom variation was interpreted in the context of dietary variation due to the divergent terrestrial ecology of D. polylepis compared to the arboreal niche occupied by all other mambas. Additionally, the degree of cross-reactivity conservation of mamba venoms was assessed by antivenomics against a panel of commercial antivenoms generated for the sub-Saharan Africa market. This study provides a genus-wide overview to infer which available antivenoms may be capable of neutralising human envenomings caused by mambas, irrespective of the species responsible. The information gathered in this study lays the foundations for rationalising the pharmacological profiles of mamba venoms at locus resolution. This understanding will contribute to the generation of a safer and more efficacious pan-Dendroaspis therapeutic antivenom against any mamba envenomation.

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.

Absolute venomics: Absolute quantification of intact venom proteins through elemental mass spectrometry

We report the application of a hybrid element and molecular MS configuration for the parallel absolute quantification of μHPLC-separated intact sulfur-containing venom proteins, via ICP triple quadrupole MS and ³²S/³⁴S isotope dilution analysis, and identification by ESI-QToF-MS of the toxins of the medically important African black-necked spitting cobra, Naja nigricollis (Tanzania); New Guinea small-eyed snake, Micropechis ikaheka; and Papuan black snake, Pseudechis papuanus. The main advantage of this approach is that only one generic sulfur-containing standard is required to quantify each and all intact Cys- and/or Met-containing toxins of the venom proteome. The results of absolute quantification are in reasonably good agreement with previously reported relative quantification of the most abundant protein families. However, both datasets depart in the quantification of the minor ones, showing a tendency for this set of proteins to be underestimated in standard peptide-centric venomics approaches. The molecular identity, specific toxic activity, and concentration in the venom, are the pillars on which the toxicovenomics-aimed discovery of the most medically-relevant venom toxins, e.g. those that need to be neutralized by an effective therapeutic antivenom, should be based. The pioneering venom proteome-wide absolute quantification shown in this paper represents thus a significant advance towards this goal. The potential of ICP triple quadrupole MS in proteomics in general, and venomics in particular, is critically discussed.

BIOLOGICAL SIGNIFICANCE

Animal venoms provide excellent model systems for investigating interactions between predators and prey, and the molecular mechanisms that contribute to adaptive protein evolution. On the other hand, numerous cases of snake bites occur yearly by encounters of humans and snakes in their shared natural environment. Snakebite envenoming is a serious global public health issue that affects the most impoverished and geopolitically disadvantaged rural communities in many tropical and subtropical countries. Unveiling the temporal and spatial patterns of venom variability is of fundamental importance to understand the molecular basis of envenoming, a prerequisite for developing therapeutic strategies against snakebite envenoming. Research on venoms has been continuously enhanced by advances in technology. The combined application of next-generation transcriptomic and venomic workflows has demonstrated unparalleled capabilities for venom characterization in unprecedented detail. However, mass spectrometry is not inherently quantitative, and this analytical limitation has sparked the development of methods to determine absolute abundance of proteins in biological samples. Here we show the potential of a hybrid element and molecular MS configuration for the parallel ESI-QToF-MS and ICP-QQQ detection and absolute quantification of intact sulfur-containing venom proteins via ³²S/³⁴S isotope dilution analysis. This configuration has been applied to quantify the toxins of the medically important African snake Naja nigricollis (Tanzania), and the Papuan species Micropechis ikaheka and Pseudechis papuanus.

Coagulating Colubrids: Evolutionary, Pathophysiological and Biodiscovery Implications of Venom Variations between Boomslang (Dispholidus typus) and Twig Snake (Thelotornis mossambicanus)

Venoms can deleteriously affect any physiological system reachable by the bloodstream, including directly interfering with the coagulation cascade. Such coagulopathic toxins may be anticoagulants or procoagulants. Snake venoms are unique in their use of procoagulant toxins for predatory purposes. The boomslang (Dispholidus typus) and the twig snakes (Thelotornis species) are iconic African snakes belonging to the family Colubridae. Both species produce strikingly similar lethal procoagulant pathologies. Despite these similarities, antivenom is only produced for treating bites by D. typus, and the mechanisms of action of both venoms have been understudied. In this study, we investigated the venom of D. typus and T. mossambicanus utilising a range of proteomic and bioactivity approaches, including determining the procoagulant properties of both venoms in relation to the human coagulation pathways. In doing so, we developed a novel procoagulant assay, utilising a Stago STA-R Max analyser, to accurately detect real time clotting in plasma at varying concentrations of venom. This approach was used to assess the clotting capabilities of the two venoms both with and without calcium and phospholipid co-factors. We found that T. mossambicanus produced a significantly stronger coagulation response compared to D. typus. Functional enzyme assays showed that T. mossambicanus also exhibited a higher metalloprotease and phospholipase activity but had a much lower serine protease activity relative to D. typus venom. The neutralising capability of the available boomslang antivenom was also investigated on both species, with it being 11.3 times more effective upon D. typus venom than T. mossambicanus. In addition to being a faster clotting venom, T. mossambicanus was revealed to be a much more complex venom composition than D. typus. This is consistent with patterns seen for other snakes with venom complexity linked to dietary complexity. Consistent with the external morphological differences in head shape between the two species, CT and MRI analyses revealed significant internal structural differences in skull architecture and venom gland anatomy. This study increases our understanding of not only the biodiscovery potential of these medically important species but also increases our knowledge of the pathological relationship between venom and the human coagulation cascade.

What killed Karl Patterson Schmidt? Combined venom gland transcriptomic, venomic and antivenomic analysis of the South African green tree snake (the boomslang), Dispholidus typus

Background

Non-front-fanged colubroid snakes comprise about two-thirds of extant ophidian species. The medical significance of the majority of these snakes is unknown, but at least five species have caused life-threatening or fatal human envenomings. However, the venoms of only a small number of species have been explored.

Methods

A combined venomic and venom gland transcriptomic approach was employed to characterise of venom of Dispholidus typus (boomslang), the snake that caused the tragic death of Professor Karl Patterson Schmidt. The ability of CroFab™ antivenom to immunocapture boomslang venom proteins was investigated using antivenomics.

Results

Transcriptomic-assisted proteomic analysis identified venom proteins belonging to seven protein families: three-finger toxin (3FTx); phospholipase A₂ (PLA₂); cysteine-rich secretory proteins (CRISP); snake venom (SV) serine proteinase (SP); C-type lectin-like (CTL); SV metalloproteinases (SVMPs); and disintegrin-like/cysteine-rich (DC) proteolytic fragments. CroFab™ antivenom efficiently immunodepleted some boomslang SVMPs.

Conclusions

The present work is the first to address the overall proteomic profile of D. typus venom. This study allowed us to correlate the toxin composition with the toxic activities of the venom. The antivenomic analysis suggested that the antivenom available at the time of the unfortunate accident could have exhibited at least some immunoreactivity against the boomslang SVMPs responsible for the disseminated intravascular coagulation syndrome that caused K.P. Schmidt's fatal outcome.

General significance

This study may stimulate further research on other non-front-fanged colubroid snake venoms capable of causing life-threatening envenomings to humans, which in turn should contribute to prevent fatal human accidents, such as that unfortunately suffered by K.P. Schmidt.

•

The venom proteome of Dispholidus typus (boomslang) is reported.

•

Transcriptomic-assisted proteomic analysis identified venom proteins belonging to seven protein families.

•

Boomslang venom proteome is dominated (75%) by snake venom PIII-metalloproteinases (PIII-SVMPs).

•

CroFab™ antivenom efficiently immunodepleted some boomslang PIII-SVMPs.

Rapid Radiations and the Race to Redundancy: An Investigation of the Evolution of Australian Elapid Snake Venoms

Australia is the stronghold of the front-fanged venomous snake family Elapidae. The Australasian elapid snake radiation, which includes approximately 100 terrestrial species in Australia, as well as Melanesian species and all the world's sea snakes, is less than 12 million years old. The incredible phenotypic and ecological diversity of the clade is matched by considerable diversity in venom composition. The clade's evolutionary youth and dynamic evolution should make it of particular interest to toxinologists, however, the majority of species, which are small, typically inoffensive, and seldom encountered by non-herpetologists, have been almost completely neglected by researchers. The present study investigates the venom composition of 28 species proteomically, revealing several interesting trends in venom composition, and reports, for the first time in elapid snakes, the existence of an ontogenetic shift in the venom composition and activity of brown snakes (Pseudonaja sp.). Trends in venom composition are compared to the snakes' feeding ecology and the paper concludes with an extended discussion of the selection pressures shaping the evolution of snake venom.

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.

Colubrid Venom Composition: An -Omics Perspective

Snake venoms have been subjected to increasingly sensitive analyses for well over 100 years, but most research has been restricted to front-fanged snakes, which actually represent a relatively small proportion of extant species of advanced snakes. Because rear-fanged snakes are a diverse and distinct radiation of the advanced snakes, understanding venom composition among "colubrids" is critical to understanding the evolution of venom among snakes. Here we review the state of knowledge concerning rear-fanged snake venom composition, emphasizing those toxins for which protein or transcript sequences are available. We have also added new transcriptome-based data on venoms of three species of rear-fanged snakes. Based on this compilation, it is apparent that several components, including cysteine-rich secretory proteins (CRiSPs), C-type lectins (CTLs), CTLs-like proteins and snake venom metalloproteinases (SVMPs), are broadly distributed among "colubrid" venoms, while others, notably three-finger toxins (3FTxs), appear nearly restricted to the Colubridae (sensu stricto). Some putative new toxins, such as snake venom matrix metalloproteinases, are in fact present in several colubrid venoms, while others are only transcribed, at lower levels. This work provides insights into the evolution of these toxin classes, but because only a small number of species have been explored, generalizations are still rather limited. It is likely that new venom protein families await discovery, particularly among those species with highly specialized diets.

Venomic Analysis of the Poorly Studied Desert Coral Snake, Micrurus tschudii tschudii, Supports the 3FTx/PLA₂ Dichotomy across Micrurus Venoms

The venom proteome of the poorly studied desert coral snake Micrurus tschudii tschudii was unveiled using a venomic approach, which identified ≥38 proteins belonging to only four snake venom protein families. The three-finger toxins (3FTxs) constitute, both in number of isoforms (~30) and total abundance (93.6% of the venom proteome), the major protein family of the desert coral snake venom. Phospholipases A₂ (PLA₂s; seven isoforms, 4.1% of the venom proteome), 1–3 Kunitz-type proteins (1.6%), and 1–2 l-amino acid oxidases (LAO, 0.7%) complete the toxin arsenal of M. t. tschudii. Our results add to the growing evidence that the occurrence of two divergent venom phenotypes, i.e., 3FTx- and PLA₂-predominant venom proteomes, may constitute a general trend across the cladogenesis of Micrurus. The occurrence of a similar pattern of venom phenotypic variability among true sea snake (Hydrophiinae) venoms suggests that the 3FTx/PLA₂ dichotomy may be widely distributed among Elapidae venoms.

Proteomic comparisons of venoms of long-term captive and recently wild-caught Eastern brown snakes (Pseudonaja textilis) indicate venom does not change due to captivity

Snake venom is a highly variable phenotypic character, and its variation and rapid evolution are important because of human health implications. Because much snake antivenom is produced from captive animals, understanding the effects of captivity on venom composition is important. Here, we have evaluated toxin profiles from six long-term (LT) captive and six recently wild-caught (RC) eastern brown snakes, Pseudonaja textilis, utilizing gel electrophoresis, HPLC-MS, and shotgun proteomics. We identified proteins belonging to the three-finger toxins, group C prothrombin activators, Kunitz-type serine protease inhibitors, and phospholipases A2, among others. Although crude venom HPLC analysis showed LT snakes to be higher in some small molecular weight toxins, presence/absence patterns showed no correlation with time in captivity. Shotgun proteomics indicated the presence of similar toxin families among individuals but with variation in protein species. Although no venom sample contained all the phospholipase A2 subunits that form the textilotoxin, all did contain both prothrombin activator subunits. This study indicates that captivity has limited effects on venom composition, that venom variation is high, and that venom composition may be correlated to geographic distribution.

BIOLOGICAL SIGNIFICANCE

Through proteomic comparisons, we show that protein variation within LT and RC groups of snakes (Pseudonaja textilis) is high, thereby resulting in no discernible differences in venom composition between groups. We utilize complementary techniques to characterize the venom proteomes of 12 individual snakes from our study area, and indicate that individuals captured close to one another have more similar venom gel electrophoresis patterns than those captured at more distant locations. These data are important for understanding natural variation in and potential effects of captivity on venom composition.

Gene expression profiling of the venom gland from the Venezuelan mapanare (Bothrops colombiensis) using expressed sequence tags (ESTs)

BACKGROUND

Bothrops colombiensis is a highly dangerous pit viper and responsible for over 70% of snakebites in Venezuela. Although the composition in B. colombiensis venom has been identified using a proteome analysis, the venom gland transcriptome is currently lacking.

RESULTS

We constructed a cDNA library from the venom gland of B. colombiensis, and a set of 729 high quality expressed sequence tags (ESTs) was identified. A total number of 344 ESTs (47.2% of total ESTs) was related to toxins. The most abundant toxin transcripts were metalloproteinases (37.5%), phospholipases A2s (PLA2, 29.7%), and serine proteinases (11.9%). Minor toxin transcripts were linked to waprins (5.5%), C-type lectins (4.1%), ATPases (2.9%), cysteine-rich secretory proteins (CRISP, 2.3%), snake venom vascular endothelium growth factors (svVEGF, 2.3%), L-amino acid oxidases (2%), and other putative toxins (1.7%). While 160 ESTs (22% of total ESTs) coded for translation proteins, regulatory proteins, ribosomal proteins, elongation factors, release factors, metabolic proteins, and immune response proteins. Other proteins detected in the transcriptome (87 ESTs, 11.9% of total ESTs) were undescribed proteins with unknown functions. The remaining 138 (18.9%) cDNAs had no match with known GenBank accessions.

CONCLUSION

This study represents the analysis of transcript expressions and provides a physical resource of unique genes for further study of gene function and the development of novel molecules for medical applications.

Deep venomics of the Pseudonaja genus reveals inter- and intra-specific variation

Australian elapid venom remains an under-investigated resource of novel bioactive peptides. In this study, the venom gland transcriptomes and proteomes of the Australian western brown snakes, Pseudonaja aspidorhyncha and Pseudonaja nuchalis, were compared to Pseudonaja textilis. A deep venomics strategy incorporating high throughput 454 pyrosequencing gave a total of 200,911 raw reads for the three venoms. Subsequent annotation identified 5716 transcripts from 20 different toxin families with inter-specific variation between species observed in eight of the less abundant families. Integration of each venom proteome with the corresponding annotated reads identified 65 isoforms from six toxin families; high sequence coverage highlighted subtle differences between sequences and intra and inter-specific variation between species. High quality MS/MS data identified unusual glycoforms with natriuretic peptides from P. aspidorhyncha and P. nuchaliscontaining O-linked trisaccharides with high homology to the glycosylated region of TNPc. Molecular evolutionary assessments indicated the accelerated evolution of all toxin families with the exception of both natriuretic peptides and P. aspidorhyncha PLA2s that were found to be evolutionarily constrained under purifying selection pressures. This study has revealed a wide range of novel peptide sequences from six bioactive peptide families and highlights the subtle differences between toxins in these closely related species.

BIOLOGICAL SIGNIFICANCE

Mining Australia's vastly untapped source of toxins from its venomous creatures has been significantly advanced by employing deep venomics methodology. Technological advances in transcriptome analysis using next generation sequencing platforms and proteome analysis by highly sensitive tandem mass spectrometry allowed a more comprehensive interrogation of three underinvestigated brown snake (Pseudonaja) venoms uncovering many novel peptide sequences that are unique to these closely related species. This generic strategy will provide invaluable information when applied to other venomous snakes for a deeper understanding of venom composition, envenomation, venom evolution, as well as identifying research tools and drug leads.

VenomKB, a new knowledge base for facilitating the validation of putative venom therapies

Animal venoms have been used for therapeutic purposes since the dawn of recorded history. Only a small fraction, however, have been tested for pharmaceutical utility. Modern computational methods enable the systematic exploration of novel therapeutic uses for venom compounds. Unfortunately, there is currently no comprehensive resource describing the clinical effects of venoms to support this computational analysis. We present VenomKB, a new publicly accessible knowledge base and website that aims to act as a repository for emerging and putative venom therapies. Presently, it consists of three database tables: (1) Manually curated records of putative venom therapies supported by scientific literature, (2) automatically parsed MEDLINE articles describing compounds that may be venom derived, and their effects on the human body, and (3) automatically retrieved records from the new Semantic Medline resource that describe the effects of venom compounds on mammalian anatomy. Data from VenomKB may be selectively retrieved in a variety of popular data formats, are open-source, and will be continually updated as venom therapies become better understood.

Phylogenetic and structural analysis of the phospholipase A2 gene family in vertebrates

The phospholipase A (PLA)2 family is the most complex gene family of phospholipases and plays a crucial role in a number of physiological activities. However, the phylogenetic background of the PLA2 gene family and the amino acid residues of the PLA2G7 gene following positive selection gene remain undetermined. In this study, we downloaded 49 genomic data sets of PLA from different species, including the human, house mouse, Norway rat, pig, dog, chicken, cattle, African clawed frog, Sumatran orangutan and the zebrafish species. Phylogenetic relationships were determined using the neighbor-joining (NJ), minimum evolution (ME) and maximum parsimony (MP) methods, as well as the Bayesian information criterion. The results were then presented as phylogenetic trees. Positive selection sites were detected using site, branch and branch‑site models. These methods led us to the following assumptions: i) closer lineages were observed between PLA2G16 and PLA2G6, PLA2G7 and PLA2G4, PLA2G3 and PLA2G12, as well as among PLA2G10, PLA2G5 and PLA2G15; ii) PLA2G5 appeared to be the origin of the PLA2 family, and PLA2G7 was one of the most evolutionarily distant PLA2 proteins; iii) 16 positive-selection sites were detected and were marked in the PLA2G7 protein sequence as 327D, 257Q, 276G, 34s, 66G, 67C, 319S, 28N, 50S, 54T, 58R, 75T, 88Q, 92R, 179H and 191K.

Venomous and poisonous Australian animals of veterinary importance: a rich source of novel therapeutics

Envenomation and poisoning by terrestrial animals (both vertebrate and invertebrate) are a significant economic problem and health risk for domestic animals in Australia. Australian snakes are some of the most venomous animals in the world and bees, wasps, ants, paralysis ticks, and cane toads are also present as part of the venomous and poisonous fauna. The diagnosis and treatment of envenomation or poisoning in animals is a challenge and can be a traumatic and expensive process for owners. Despite the potency of Australian venoms, there is potential for novel veterinary therapeutics to be modeled on venom toxins, as has been the case with human pharmaceuticals. A comprehensive overview of envenomation and poisoning signs in livestock and companion animals is provided and related to the potential for venom toxins to act as therapeutics.