Evolution of an Arsenal

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.

Are Fireworms Venomous? Evidence for the Convergent Evolution of Toxin Homologs in Three Species of Fireworms (Annelida, Amphinomidae)

Amphinomids, more commonly known as fireworms, are a basal lineage of marine annelids characterized by the presence of defensive dorsal calcareous chaetae, which break off upon contact. It has long been hypothesized that amphinomids are venomous and use the chaetae to inject a toxic substance. However, studies investigating fireworm venom from a morphological or molecular perspective are scarce and no venom gland has been identified to date, nor any toxin characterized at the molecular level. To investigate this question, we analyzed the transcriptomes of three species of fireworms-Eurythoe complanata, Hermodice carunculata, and Paramphinome jeffreysii-following a venomics approach to identify putative venom compounds. Our venomics pipeline involved de novo transcriptome assembly, open reading frame, and signal sequence prediction, followed by three different homology search strategies: BLAST, HMMER sequence, and HMMER domain. Following this pipeline, we identified 34 clusters of orthologous genes, representing 13 known toxin classes that have been repeatedly recruited into animal venoms. Specifically, the three species share a similar toxin profile with C-type lectins, peptidases, metalloproteinases, spider toxins, and CAP proteins found among the most highly expressed toxin homologs. Despite their great diversity, the putative toxins identified are predominantly involved in three major biological processes: hemostasis, inflammatory response, and allergic reactions, all of which are commonly disrupted after fireworm stings. Although the putative fireworm toxins identified here need to be further validated, our results strongly suggest that fireworms are venomous animals that use a complex mixture of toxins for defense against predators.

Characterization of a double WAP domain-containing protein from the red swamp crayfish Procambarus clarkii

Crustaceans express multiple whey acidic protein (WAP) domain containing proteins which are components of host immunity. In the present study, a new double WAP domain containing protein was identified from red swamp crayfish Procambarus clarkii, designated Pc-DWD. The ORF is 387 bp, encoding 128 amino acids consisting of signal peptide of 18 residues, and two tandem WAP domains of 38 and 44 residues. Multiple alignment indicates the presence of conserved motifs in both WAP domains, and phylogenetic analysis shows that Pc-DWD is a new member of the type-IV crustin family. Pc-DWD transcripts were found most abundantly in hemocytes, gills, intestine and heart, and induced by Vibrio anguillarum, Staphylococcus aureus and white spot syndrome virus challenge. RNAi knockdown of Pc-DWD expression led to increased expression of white spot syndrome virus genes and increased crayfish mortality after virus infection. Recombinant Pc-DWD exhibited strong protease inhibitory activity towards commercial subtilicin A and protease K. Pc-DWD inhibited the crude proteases from V. anguillarum and S. aureus cultures and from the crayfish tissue extracts. We infer that Pc-DWD acts in crayfish bacterial and viral immunity.

Has snake fang evolution lost its bite? New insights from a structural mechanics viewpoint

Venomous snakes-the pinnacle of snake evolution-are characterized by their possession of venom-conducting fangs ranging from grooved phenotypes characterizing multiple lineages of rear-fanged taxa to tubular phenotypes present in elapids, viperids and atractaspidines. Despite extensive research, controversy still exists on the selective pressures involved in fang phenotype diversification. Here, we test the hypothesis that larger fangs and consequently a shift to an anterior position in the maxilla evolved to compensate for the costs of structural changes, i.e. higher stress upon impact in tubular fangs compared to grooved fangs. Direct voxel-based stress simulations conducted on high-resolution µCT scans, analysed within a phylogenetic framework, showed no differences in stress distribution between the three fang phenotypes, despite differences in (relative) fang length. These findings suggest that additional compensatory mechanisms are responsible for the biomechanical optimization and that fang length might instead be related to differential striking behaviour strategies.

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.

Enter the Dragon: The Dynamic and Multifunctional Evolution of Anguimorpha Lizard Venoms

While snake venoms have been the subject of intense study, comparatively little work has been done on lizard venoms. In this study, we have examined the structural and functional diversification of anguimorph lizard venoms and associated toxins, and related these results to dentition and predatory ecology. Venom composition was shown to be highly variable across the 20 species of Heloderma, Lanthanotus, and Varanus included in our study. While kallikrein enzymes were ubiquitous, they were also a particularly multifunctional toxin type, with differential activities on enzyme substrates and also ability to degrade alpha or beta chains of fibrinogen that reflects structural variability. Examination of other toxin types also revealed similar variability in their presence and activity levels. The high level of venom chemistry variation in varanid lizards compared to that of helodermatid lizards suggests that venom may be subject to different selection pressures in these two families. These results not only contribute to our understanding of venom evolution but also reveal anguimorph lizard venoms to be rich sources of novel bioactive molecules with potential as drug design and development lead compounds.

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.

Haemotoxic snake venoms: their functional activity, impact on snakebite victims and pharmaceutical promise

Snake venoms are mixtures of numerous proteinacious components that exert diverse functional activities on a variety of physiological targets. Because the toxic constituents found in venom vary from species to species, snakebite victims can present with a variety of life-threatening pathologies related to the neurotoxic, cytotoxic and haemotoxic effects of venom. Of the 1·8 million people envenomed by snakes every year, up to 125 000 die, while hundreds of thousands survive only to suffer with life-changing long-term morbidity. Consequently, snakebite is one of the world's most severe neglected tropical diseases. Many snake venoms exhibit strong haemotoxic properties by interfering with blood pressure, clotting factors and platelets, and by directly causing haemorrhage. In this review we provide an overview of the functional activities of haemotoxic venom proteins, the pathologies they cause in snakebite victims and how their exquisite selectivity and potency make them amenable for use as therapeutic and diagnostic tools relevant for human medicine.

Disseminated intravascular coagulation caused by moojenactivase, a procoagulant snake venom metalloprotease

Snake venom toxins that activate coagulation factors are key players in the process of venom-induced coagulopathy, and account for severe clinical manifestations. The present study applies a variety of biochemical, hematological, and histopathological approaches to broadly investigate the intravascular and systemic effects of moojenactivase (MooA), the first described PIIId subclass metalloprotease isolated from Bothrops sp. venom that activates coagulation factors. MooA induced consumption coagulopathy with high toxic potency, characterized by prolongation of prothrombin and activated partial thromboplastin time, consumption of fibrinogen and the plasma coagulation factors X and II, and thrombocytopenia. MooA promoted leukocytosis and expression of the proinflammatory cytokines interleukin-6 and tumor necrosis factor-α, accompanied by tissue factor-dependent procoagulant activity in peripheral blood mononuclear cells. This metalloprotease also caused intravascular hemolysis, elevated plasma levels of creatine kinase-MB, aspartate transaminase, and urea/creatinine, and induced morphopathological alterations in erythrocytes, heart, kidney, and lungs associated with thrombosis and hemorrhage. Diagnosis of MooA-induced disseminated intravascular coagulation represents an important approach to better understand the pathophysiology of Bothrops envenomation and develop novel therapeutic strategies targeting hemostatic disturbances.

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.

Selection To Increase Expression, Not Sequence Diversity, Precedes Gene Family Origin and Expansion in Rattlesnake Venom

Gene duplication is the primary mechanism leading to new genes and phenotypic novelty, but the proximate evolutionary processes underlying gene family origin, maintenance, and expansion are poorly understood. Although sub- and neofunctionalization provide clear long-term advantages, selection does not act with foresight, and unless a redundant gene copy provides an immediate fitness advantage, the copy will most likely be lost. Many models for the evolution of genes immediately following duplication have been proposed, but the robustness and applicability of these models is unclear because of the lack of data at the population level. We used qPCR, protein expression data, genome sequencing, and hybrid enrichment to test three competing models that differ in whether selection favoring the spread of duplicates acts primarily on expression level or sequence diversity for specific toxin-encoding loci in the eastern diamondback rattlesnake (Crotalus adamanteus). We sampled 178 individuals and identified significant inter- and intrapopulation variation in copy number, demonstrated that copy number was significantly and positively correlated with protein expression, and found little to no sequence variation across paralogs in all populations. Collectively, these results demonstrate that selection for increased expression, not sequence diversity, was the proximate evolutionary process underlying gene family origin and expansion, providing data needed to resolve the debate over which evolutionary processes govern the fates of gene copies immediately following duplication.

The three-finger toxin fold: a multifunctional structural scaffold able to modulate cholinergic functions

Three-finger fold toxins are miniproteins frequently found in Elapidae snake venoms. This fold is characterized by three distinct loops rich in β-strands and emerging from a dense, globular core reticulated by four highly conserved disulfide bridges. The number and diversity of receptors, channels, and enzymes identified as targets of three-finger fold toxins is increasing continuously. Such manifold diversity highlights the specific adaptability of this fold for generating pleiotropic functions. Although this toxin superfamily disturbs many biological functions by interacting with a large diversity of molecular targets, the most significant target is the cholinergic system. By blocking the activity of the nicotinic and muscarinic acetylcholine receptors or by inhibiting the enzyme acetylcholinesterase, three-finger fold toxins interfere most drastically with neuromuscular junction functioning. Several of these toxins have become powerful pharmacological tools for studying the function and structure of their molecular targets. Most importantly, since dysfunction of these receptors/enzyme is involved in many diseases, exploiting the three-finger scaffold to create novel, highly specific therapeutic agents may represent a major future endeavor. This is an article for the special issue XVth International Symposium on Cholinergic Mechanisms.

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.

Viper Venom Botox: The Molecular Origin and Evolution of the Waglerin Peptides Used in Anti-Wrinkle Skin Cream

The molecular origin of waglerin peptides has remained enigmatic despite their industrial application in skin cream products to paralyse facial muscles and thus reduce the incidence of wrinkles. Here we show that these neurotoxic peptides are the result of de novo evolution within the prepro region of the C-type natriuretic peptide gene in Tropidolaemus venoms, at a site distinct from the domain encoding for the natriuretic peptide. It is the same region that yielded the azemiopsin peptides from Azemiops feae, indicative of a close relationship of this toxin gene between these two genera. The precursor region for the molecular evolution is a biodiversity hotspot that has yielded other novel bioactive peptides with novel activities. We detail the diversity of components in this and other species in order to explore what characteristics enable it to be such a biodiscovery treasure trove. The unusual function of Tropidolaemus venoms may have been selected for due to evolutionary pressures brought about by a high likelihood of prey escape.

Small Packages, Big Returns: Uncovering the Venom Diversity of Small Invertebrate Conoidean Snails

Venomous organisms used in research were historically chosen based on size and availability. This opportunity-driven strategy created a species bias in which snakes, scorpions, and spiders became the primary subjects of venom research. Increasing technological advancements have enabled interdisciplinary studies using genomics, transcriptomics, and proteomics to expand venom investigation to animals that produce small amounts of venom or lack traditional venom producing organs. One group of non-traditional venomous organisms that have benefitted from the rise of -omic technologies is the Conoideans. The Conoidean superfamily of venomous marine snails includes, the Terebridae, Turridae (s.l), and Conidae. Conoidea venom is used for both predation and defense, and therefore under strong selection pressures. The need for conoidean venom peptides to be potent and specific to their molecular targets has made them important tools for investigating cellular physiology and bioactive compounds that are beneficial to improving human health. A convincing case for the potential of Conoidean venom is made with the first commercially available conoidean venom peptide drug Ziconotide (Prialt®), an analgesic derived from Conus magus venom that is used to treat chronic pain in HIV and cancer patients. Investigation of conoidean venom using -omics technology provides significant insights into predator-driven diversification in biodiversity and identifies novel compounds for manipulating cellular communication, especially as it pertains to disease and disorders.

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.

[Venoms and medical research]

Animal venoms are complex chemical cocktails, comprising a wide range of biologically active reticulated peptides that target with high selectivity and efficacy a variety of enzymes, membrane receptors, ion channels...Venoms can therefore be seen as large natural libraries of biologically active molecules that are continuously selected and highly refined by the evolution process, up to the point where every molecule is endowed with pharmacological properties that are highly valuable in the context of human use and drug development. Therefore, venom exploration constitutes a prerequisite to drug discovery. However, mass spectrometry and transcriptomics via NGS (Next Generation Sequencing) studies have shown the presence of up to 1000 peptides in the venom of single species of cone snails and spiders. Therefore the global animal venom resource can be seen as a collection of more than 50 to 100 000 000 peptides and proteins of which only ~5000 are known. That extraordinary "Eldorado" of bio-optimized compounds justifies the development of more global and cutting-edge strategies and technologies to explore this resource more efficiently than actually. De novo developed approaches and recently obtained results will be described.

The Deep Origin and Recent Loss of Venom Toxin Genes in Rattlesnakes

The genetic origin of novel traits is a central but challenging puzzle in evolutionary biology. Among snakes, phospholipase A2 (PLA2)-related toxins have evolved in different lineages to function as potent neurotoxins, myotoxins, or hemotoxins. Here, we traced the genomic origin and evolution of PLA2 toxins by examining PLA2 gene number, organization, and expression in both neurotoxic and non-neurotoxic rattlesnakes. We found that even though most North American rattlesnakes do not produce neurotoxins, the genes of a specialized heterodimeric neurotoxin predate the origin of rattlesnakes and were present in their last common ancestor (∼22 mya). The neurotoxin genes were then deleted independently in the lineages leading to the Western Diamondback (Crotalus atrox) and Eastern Diamondback (C. adamanteus) rattlesnakes (∼6 mya), while a PLA2 myotoxin gene retained in C. atrox was deleted from the neurotoxic Mojave rattlesnake (C. scutulatus; ∼4 mya). The rapid evolution of PLA2 gene number appears to be due to transposon invasion that provided a template for non-allelic homologous recombination.

Trends in the Evolution of Snake Toxins Underscored by an Integrative Omics Approach to Profile the Venom of the Colubrid Phalotris mertensi

Only few studies on snake venoms were dedicated to deeply characterize the toxin secretion of animals from the Colubridae family, despite the fact that they represent the majority of snake diversity. As a consequence, some evolutionary trends observed in venom proteins that underpinned the evolutionary histories of snake toxins were based on data from a minor parcel of the clade. Here, we investigated the proteins of the totally unknown venom from Phalotris mertensi (Dipsadinae subfamily), in order to obtain a detailed profile of its toxins and to appreciate evolutionary tendencies occurring in colubrid venoms. By means of integrated omics and functional approaches, including RNAseq, Sanger sequencing, high-resolution proteomics, recombinant protein production, and enzymatic tests, we verified an active toxic secretion containing up to 21 types of proteins. A high content of Kunitz-type proteins and C-type lectins were observed, although several enzymatic components such as metalloproteinases and an L-amino acid oxidase were also present in the venom. Interestingly, an arguable venom component of other species was demonstrated as a true venom protein and named svLIPA (snake venom acid lipase). This finding indicates the importance of checking the actual protein occurrence across species before rejecting genes suggested to code for toxins, which are relevant for the discussion about the early evolution of reptile venoms. Moreover, trends in the evolution of some toxin classes, such as simplification of metalloproteinases and rearrangements of Kunitz and Wap domains, parallel similar phenomena observed in other venomous snake families and provide a broader picture of toxin evolution.

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.

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.

Stabilising the Integrity of Snake Venom mRNA Stored under Tropical Field Conditions Expands Research Horizons

BACKGROUND

Snake venoms contain many proteinaceous toxins that can cause severe pathology and mortality in snakebite victims. Interestingly, mRNA encoding such toxins can be recovered directly from venom, although yields are low and quality is unknown. It also remains unclear whether such RNA contains information about toxin isoforms and whether it is representative of mRNA recovered from conventional sources, such as the venom gland. Answering these questions will address the feasibility of using venom-derived RNA for future research relevant to biomedical and antivenom applications.

METHODOLOGY/PRINCIPAL FINDINGS

Venom was extracted from several species of snake, including both members of the Viperidae and Elapidae, and either lyophilized or immediately added to TRIzol reagent. TRIzol-treated venom was incubated at a range of temperatures (4-37°C) for a range of durations (0-48 hours), followed by subsequent RNA isolation and assessments of RNA quantity and quality. Subsequently, full-length toxin transcripts were targeted for PCR amplification and Sanger sequencing. TRIzol-treated venom yielded total RNA of greater quantity and quality than lyophilized venom, and with quality comparable to venom gland-derived RNA. Full-length sequences from multiple Viperidae and Elapidae toxin families were successfully PCR amplified from TRIzol-treated venom RNA. We demonstrated that venom can be stored in TRIzol for 48 hours at 4-19°C, and 8 hours at 37°C, at minimal cost to RNA quality, and found that venom RNA encoded multiple toxin isoforms that seemed homologous (98-99% identity) to those found in the venom gland.

CONCLUSIONS/SIGNIFICANCE

The non-invasive experimental modifications we propose will facilitate the future investigation of venom composition by using venom as an alternative source to venom gland tissue for RNA-based studies, thus obviating the undesirable need to sacrifice snakes for such research purposes. In addition, they expand research horizons to rare, endangered or protected snake species and provide more flexibility to performing fieldwork on venomous snakes in tropical conditions.

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.

[Relationships between venomous function and innate immune function]

Venomous function is investigated in relation to innate immune function in two cases selected from scorpion venom and serpent venom. In the first case, structural analysis of scorpion toxins and defensins reveals a close interrelation between both functions (toxic and innate immune system function). In the second case, structural and functional studies of natural inhibitors of toxic snake venom phospholipases A2 reveal homology with components of the innate immune system, leading to a similar conclusion. Although there is a clear functional distinction between neurotoxins, which act by targeting membrane ion channels, and the circulating defensins which protect the organism from pathogens, the scorpion short toxins and defensins share a common protein folding scaffold with a conserved cysteine-stabilized alpha-beta motif of three disulfide bridges linking a short alpha helix and an antiparallel beta sheet. Genomic analysis suggests that these proteins share a common ancestor (long venom toxins were separated from an early gene family which gave rise to separate short toxin and defensin families). Furthermore, a scorpion toxin has been experimentally synthetized from an insect defensin, and an antibacterial scorpion peptide, androctonin (whose structure is similar to that of a cone snail venom toxin), was shown to have a similar high affinity for the postsynaptic acetylcholine receptor of Torpedo sp. Natural inhibitors of phospholipase A2 found in the blood of snakes are associated with the resistance of venomous snakes to their own highly neurotoxic venom proteins. Three classes of phospholipases A2 inhibitors (PLI-α, PLI-β, PLI-γ) have been identified. These inhibitors display diverse structural motifs related to innate immune proteins including carbohydrate recognition domains (CRD), leucine rich repeat domains (found in Toll-like receptors) and three finger domains, which clearly differentiate them from components of the adaptive immune system. Thus, in structure, function and phylogeny, venomous function in both vertebrates and invertebrates are clearly interrelated with innate immune function.

Expression Differentiation Is Constrained to Low-Expression Proteins over Ecological Timescales

Protein expression level is one of the strongest predictors of protein sequence evolutionary rate, with high-expression protein sequences evolving at slower rates than low-expression protein sequences largely because of constraints on protein folding and function. Expression evolutionary rates also have been shown to be negatively correlated with expression level across human and mouse orthologs over relatively long divergence times (i.e., ∼100 million years). Long-term evolutionary patterns, however, often cannot be extrapolated to microevolutionary processes (and vice versa), and whether this relationship holds for traits evolving under directional selection within a single species over ecological timescales (i.e., <5000 years) is unknown and not necessarily expected. Expression is a metabolically costly process, and the expression level of a particular protein is predicted to be a tradeoff between the benefit of its function and the costs of its expression. Selection should drive the expression level of all proteins close to values that maximize fitness, particularly for high-expression proteins because of the increased energetic cost of production. Therefore, stabilizing selection may reduce the amount of standing expression variation for high-expression proteins, and in combination with physiological constraints that may place an upper bound on the range of beneficial expression variation, these constraints could severely limit the availability of beneficial expression variants. To determine whether rapid-expression evolution was restricted to low-expression proteins owing to these constraints on highly expressed proteins over ecological timescales, we compared venom protein expression levels across mainland and island populations for three species of pit vipers. We detected significant differentiation in protein expression levels in two of the three species and found that rapid-expression differentiation was restricted to low-expression proteins. Our results suggest that various constraints on high-expression proteins reduce the availability of beneficial expression variants relative to low-expression proteins, enabling low-expression proteins to evolve and potentially lead to more rapid adaptation.

Identification and molecular characterization of five putative toxins from the venom gland of the snake Philodryas chamissonis (Serpentes: Dipsadidae)

Philodryas chamissonis is a rear-fanged snake endemic to Chile. Its bite produces mild to moderate symptoms with proteolytic and anti-coagulant effects. Presently, the composition of the venom, as well as, the biochemical and structural characteristics of its toxins, remains unknown. In this study, we cloned and reported the first full-length sequences of five toxin-encoding genes from the venom gland of this species: Type III snake venom metalloprotease (SVMP), snake venom serine protease (SVSP), Cysteine-rich secretory protein (CRISP), α and β subunits of C-type lectin-like protein (CLP) and C-type natriuretic peptide (NP). These genes are highly expressed in the venom gland and their sequences exhibited a putative signal peptide, suggesting that these are components of the venom. These putative toxins had different evolutionary relationships with those reported for some front-fanged snakes, being SVMP, SVSP and CRISP of P. chamissonis closely related to the toxins present in Elapidae species, while NP was more related to those of Viperidae species. In addition, analyses suggest that the α and β subunits of CLP of P. chamissonis might have a α-subunit scaffold in common with Viperidae species, whose highly variable C-terminal region might have allowed the diversification in α and β subunits. Our results provide the first molecular description of the toxins possibly implicated in the envenomation of prey and humans by the bite of P. chamissonis.

Discovery of toxin-encoding genes from the false viper Macropisthodon rudis, a rear-fanged snake, by transcriptome analysis of venom gland

Although rear-fanged snakes are often considered as non-threatening to humans, some species are lethal or medically hazardous. The toxin components and bioactivities of front-fanged snakes have been extensively studied; however, only limited research has explored the venoms of rear-fanged snakes. The false viper, Macropisthodon rudis, is widespread in southern China, but little is known about the toxins that this snake produces. Here, we analyzed the transcriptome of the venom gland of M. rudis using high-throughput sequencing with an illumina HiSeq 2000. The raw data were assembled and annotated using public databases. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and gene ontology (GO) were analyzed. Using sequence comparisons, snake venom metalloproteinases (SVMPs) and a phosphodiesterase (PDE) were discovered in the venom gland of M. rudis.

Cobra venom proteome and glycome determined from individual snakes of Naja atra reveal medically important dynamic range and systematic geographic variation

Recent progress in snake venomics has shed much light on the intra-species variation among the toxins from different geographical regions and has provided important information for better snakebite management. Most previous reports on snake venomics were based on venoms pooled from different snakes. In this study, we present the proteomic and glycomic profiles of venoms from individual Naja atra snakes. The results reveal wide dynamic range of three-finger toxins. Systematic classification based on cardiotoxin (CTX-) profiles of A2/A4 and A6, respectively, allowed the identification of two putative subspecies of Taiwan cobra from the eastern and western regions. We also identified four major N-glycan moieties on cobra snake venom metalloproteinase on the bi-antennary glycan core. ELISA showed that these glycoproteins (<3%) could elicit much higher antibody response in antiserum when compared to other high-abundance cobra venom toxins such as small molecular weight CTXs (~60%). By removing these high-molecular weight glycoproteins from the immunogen, we demonstrated better protection than that achieved with conventional crude venom immunization in mice challenged by crude venom. We conclude that both intra-species and inter-individual variations of proteomic and glycomic profiles of snake venomics should be considered to provide better antivenomic approach for snakebite management.

BIOLOGICAL SIGNIFICANCE

Based on the proteomic and glycomic profiles of venoms obtained from individual snakes, we demonstrated a surprisingly wide dynamic range and geographical variation of three-finger toxins in cobra venomics. This provides a reasonable explanation for the variable neutralization effects of antivenom treatment on victims suffering from cobra snakebite and suggests a simple and economic method to produce potent antivenom with better efficacy. Since two major venomic profiles with distinct dynamic ranges were observed for Taiwan cobra venoms isolated from the eastern and western regions, the current venomic profile should be used as a quality control for future production of antivenom in clinical applications.

Cabinet of Curiosities: Venom Systems and Their Ecological Function in Mammals, with a Focus on Primates

Venom delivery systems (VDS) are common in the animal kingdom, but rare amongst mammals. New definitions of venom allow us to reconsider its diversity amongst mammals by reviewing the VDS of Chiroptera, Eulipotyphla, Monotremata, and Primates. All orders use modified anterior dentition as the venom delivery apparatus, except Monotremata, which possesses a crural system. The venom gland in most taxa is a modified submaxillary salivary gland. In Primates, the saliva is activated when combined with brachial gland exudate. In Monotremata, the crural spur contains the venom duct. Venom functions include feeding, intraspecific competition, anti-predator defense and parasite defense. Including mammals in discussion of venom evolution could prove vital in our understanding protein functioning in mammals and provide a new avenue for biomedical and therapeutic applications and drug discovery.

Transcriptomics and venomics: implications for medicinal chemistry

Over the last three decades, transcriptomic studies of venom gland cells have continuously evolved, opening up new possibilities for exploring the molecular diversity of animal venoms, a prerequisite for the discovery of new drug candidates and molecular phylogenetics. The molecular complexity of animal venoms is much greater than initially thought. In this review, we describe the different technologies available for transcriptomic studies of venom, from the original individual cloning approaches to the more recent global Next Generation Sequencing strategies. Our understanding of animal venoms is evolving, with the discovery of complex and diverse bio-optimized cocktails of compounds, including mostly peptides and proteins, which are now beginning to be studied by academic and industrial researchers.

Synthesis and characterization of an antibacterial and non-toxic dimeric peptide derived from the C-terminal region of Bothropstoxin-I

Infectious diseases are among the leading global causes of death, increasing the search for novel antibacterial agents. Among these, biologically active peptides are an excellent research tool. Using solid-phase peptide synthesis (SPPS), this work aimed to synthesize the peptide derived from the C-terminal region of Bothropstoxin-I (BthTX-I) (p-BthTX-I, sequence: KKYRYHLKPFCKK), and its disulfide-linked dimeric form, obtained via air oxidation (p-BthTX-I)2. Two other peptides were synthesized to evaluate the dimerization effect on antimicrobial activity. In both sequences, the cysteine (Cys) residue was replaced by the serine (Ser) residue, differing, however, in their C-terminus position. The antimicrobial activity of the peptides against gram-negative (Escherichia (E.) coli) and gram-positive (Staphylococcus (S.) aureus) bacteria and yeast (Candida (C.) albicans) was evaluated. Interestingly, only peptides containing the Cys residue showed antimicrobial activity, suggesting the importance of Cys residue and its dimerization for the observed activity. Apparently, p-BthTX-I and (p-BthTX-I)2 did not promote lysis or form pores and were not able to interact with membranes. Furthermore, they neither showed antifungal activity against C. albicans nor toxicity against erythrocytes, epithelial cells, or macrophages, indicating a potential specificity against prokaryotic cells.

Weaponization of a Hormone: Convergent Recruitment of Hyperglycemic Hormone into the Venom of Arthropod Predators

Arthropod venoms consist primarily of peptide toxins that are injected into their prey with devastating consequences. Venom proteins are thought to be recruited from endogenous body proteins and mutated to yield neofunctionalized toxins with remarkable affinity for specific subtypes of ion channels and receptors. However, the evolutionary history of venom peptides remains poorly understood. Here we show that a neuropeptide hormone has been convergently recruited into the venom of spiders and centipedes and evolved into a highly stable toxin through divergent modification of the ancestral gene. High-resolution structures of representative hormone-derived toxins revealed they possess a unique structure and disulfide framework and that the key structural adaptation in weaponization of the ancestral hormone was loss of a C-terminal α helix, an adaptation that occurred independently in spiders and centipedes. Our results raise a new paradigm for toxin evolution and highlight the value of structural information in providing insight into protein evolution.

Privileged frameworks from snake venom

Venom as a form of chemical prey capture is a key innovation that has underpinned the explosive radiation of the advanced snakes (Caenophidia). Small venom proteins are often rich in disulfide bonds thus facilitating stable molecular scaffolds that present key functional residues on the protein surface. New toxin types are initially developed through the venom gland over-expression of normal body proteins, their subsequent gene duplication and diversification that leads to neofunctionalisation as random mutations modify their structure and function. This process has led to preferentially selected (privileged) cysteine-rich scaffolds that enable the snake to build arrays of toxins many of which may lead to therapeutic products and research tools. This review focuses on cysteine-rich small proteins and peptides found in snake venoms spanning natriuretic peptides to phospholipase enzymes, while highlighting their three-dimensional structures and biological functions as well as their potential as therapeutic agents or research tools.

Tissue localization and extracellular matrix degradation by PI, PII and PIII snake venom metalloproteinases: clues on the mechanisms of venom-induced hemorrhage

Snake venom hemorrhagic metalloproteinases (SVMPs) of the PI, PII and PIII classes were compared in terms of tissue localization and their ability to hydrolyze basement membrane components in vivo, as well as by a proteomics analysis of exudates collected in tissue injected with these enzymes. Immunohistochemical analyses of co-localization of these SVMPs with type IV collagen revealed that PII and PIII enzymes co-localized with type IV collagen in capillaries, arterioles and post-capillary venules to a higher extent than PI SVMP, which showed a more widespread distribution in the tissue. The patterns of hydrolysis by these three SVMPs of laminin, type VI collagen and nidogen in vivo greatly differ, whereas the three enzymes showed a similar pattern of degradation of type IV collagen, supporting the concept that hydrolysis of this component is critical for the destabilization of microvessel structure leading to hemorrhage. Proteomic analysis of wound exudate revealed similarities and differences between the action of the three SVMPs. Higher extent of proteolysis was observed for the PI enzyme regarding several extracellular matrix components and fibrinogen, whereas exudates from mice injected with PII and PIII SVMPs had higher amounts of some intracellular proteins. Our results provide novel clues for understanding the mechanisms by which SVMPs induce damage to the microvasculature and generate hemorrhage.

Local and systemic hemorrhage are typical manifestations of envenomings by viperid snakes. Hemorrhagic activity is due to the action of snake venom metalloproteinases (SVMPs) in the microvasculature, especially in capillaries. There are large differences in the hemorrhagic potential of SVMPs, depending on their domain composition. The present study compared PI, PII and PIII hemorrhagic SVMPs for their tissue distribution and their ability to cleave proteins of the extracellular matrix (ECM), especially those of the basement membrane (BM) that provides mechanical stability to microvessels. Observations indicate that PII and PIII SVMPs, which exert a high hemorrhagic activity, are preferentially located in microvessels, whereas PI SVMP is distributed in a more widespread fashion in the tissue. In addition, when these toxins are injected at doses that induce a similar hemorrhagic effect, they cleave type IV collagen to a similar extent, showing differences in the cleavage patterns of other ECM components, such as laminin, nidogen and type VI collagen. The analysis of the exudates resultant from the action of these SVMPs in the tissue revealed many similarities and some differences in the exudate proteomes. Overall our results indicate that hydrolysis of type IV collagen is a key event in the onset of microvessel damage, and that the ability of SVMPs to bind to microvessels greatly determines their hemorrhagic potential.

Restriction and recruitment-gene duplication and the origin and evolution of snake venom toxins

Snake venom has been hypothesized to have originated and diversified through a process that involves duplication of genes encoding body proteins with subsequent recruitment of the copy to the venom gland, where natural selection acts to develop or increase toxicity. However, gene duplication is known to be a rare event in vertebrate genomes, and the recruitment of duplicated genes to a novel expression domain (neofunctionalization) is an even rarer process that requires the evolution of novel combinations of transcription factor binding sites in upstream regulatory regions. Therefore, although this hypothesis concerning the evolution of snake venom is very unlikely and should be regarded with caution, it is nonetheless often assumed to be established fact, hindering research into the true origins of snake venom toxins. To critically evaluate this hypothesis, we have generated transcriptomic data for body tissues and salivary and venom glands from five species of venomous and nonvenomous reptiles. Our comparative transcriptomic analysis of these data reveals that snake venom does not evolve through the hypothesized process of duplication and recruitment of genes encoding body proteins. Indeed, our results show that many proposed venom toxins are in fact expressed in a wide variety of body tissues, including the salivary gland of nonvenomous reptiles and that these genes have therefore been restricted to the venom gland following duplication, not recruited. Thus, snake venom evolves through the duplication and subfunctionalization of genes encoding existing salivary proteins. These results highlight the danger of the elegant and intuitive “just-so story” in evolutionary biology.