Profiling the venom gland transcriptomes of Costa Rican snakes by 454 pyrosequencing

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

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

Sequence Divergence in Venom Genes Within and Between Montane Pitviper (Viperidae: Crotalinae: Cerrophidion) Species is Driven by Mutation-Drift Equilibrium

Snake venom can vary both among and within species. While some groups of New World pitvipers-such as rattlesnakes-have been well studied, very little is known about the venom of montane pitvipers (Cerrophidion) found across the Mesoamerican highlands. Compared to most well-studied rattlesnakes, which are widely distributed, the isolated montane populations of Cerrophidion may facilitate unique evolutionary trajectories and venom differentiation. Here, we describe the venom gland transcriptomes for populations of C. petlalcalensis, C. tzotzilorum, and C. godmani from Mexico, and a single individual of C. sasai from Costa Rica. We explore gene expression variation in Cerrophidion and sequence evolution of toxins within C. godmani specifically. Cerrophidion venom gland transcriptomes are composed primarily of snake venom metalloproteinases, phospholipase A[Formula: see text]s (PLA[Formula: see text]s), and snake venom serine proteases. Cerrophidion petlalcalensis shows little intraspecific variation; however, C. godmani and C. tzotzilorum differ significantly between geographically isolated populations. Interestingly, intraspecific variation was mostly attributed to expression variation as we did not detect signals of selection within C. godmani toxins. Additionally, we found PLA[Formula: see text]-like myotoxins in all species except C. petlalcalensis, and crotoxin-like PLA[Formula: see text]s in the southern population of C. godmani. Our results demonstrate significant intraspecific venom variation within C. godmani and C. tzotzilorum. The toxins of C. godmani show little evidence of directional selection where variation in toxin sequence is consistent with evolution under a model of mutation-drift equilibrium. Cerrophidion godmani individuals from the southern population may exhibit neurotoxic venom activity given the presence of crotoxin-like PLA[Formula: see text]s; however, further research is required to confirm this hypothesis.

Snake Antivenoms-Toward Better Understanding of the Administration Route

Envenomations induced by animal bites and stings constitute a significant public health burden. Even though a standardized protocol does not exist, parenterally administered polyclonal antivenoms remain the mainstay in snakebite therapy. There is a prevailing opinion that their application by the i.m. route has poor efficacy and that i.v. administration should preferentially be chosen in order to achieve better accomplishment of the antivenom therapeutic activity. Recently, it has been demonstrated that neutralization not only in the systemic circulation but also in the lymphatic system might be of great importance for the clinical outcome since it represents another relevant body compartment through which the absorption of the venom components occurs. In this review, the present-day and summarized knowledge of the laboratory and clinical findings on the i.v. and i.m. routes of antivenom administration is provided, with a special emphasis on the contribution of the lymphatic system to the process of venom elimination. Until now, antivenom-mediated neutralization has not yet been discussed in the context of the synergistic action of both blood and lymph. A current viewpoint might help to improve the comprehension of the venom/antivenom pharmacokinetics and the optimal approach for drug application. There is a great need for additional dependable, practical, well-designed studies, as well as more practice-related experience reports. As a result, opportunities for resolving long-standing disputes over choosing one therapeutic principle over another might be created, improving the safety and effectiveness of snakebite management.

Commercial Antivenoms Exert Broad Paraspecific Immunological Binding and In Vitro Inhibition of Medically Important Bothrops Pit Viper Venoms

Snakebite envenoming is a life threatening neglected tropical disease that represents a considerable public health concern in the tropics. Viperid snakes of the genus Bothrops are among those of greatest medical importance in Latin America, and they frequently cause severe systemic haemotoxicity and local tissue destructive effects in human victims. Although snakebite antivenoms can be effective therapeutics, their efficacy is undermined by venom toxin variation among snake species. In this study we investigated the extent of paraspecific venom cross-reactivity exhibited by three distinct anti-Bothrops antivenoms (Soro antibotrópico-crotálico, BothroFav and PoliVal-ICP) against seven different Bothrops pit viper venoms from across Latin America. We applied a range of in vitro assays to assess the immunological binding and recognition of venom toxins by the antivenoms and their inhibitory activities against specific venom functionalities. Our findings demonstrated that, despite some variations, the monovalent antivenom BothroFav and the polyvalent antivenoms Soro antibotrópico-crotálico and PoliVap-ICP exhibited extensive immunological recognition of the distinct toxins found in the different Bothrops venoms, with Soro antibotrópico-crotálico generally outperformed by the other two products. In vitro functional assays revealed outcomes largely consistent with the immunological binding data, with PoliVap-ICP and BothroFav exhibiting the greatest inhibitory potencies against procoagulant and fibrinogen-depleting venom activities, though Soro antibotrópico-crotálico exhibited potent inhibition of venom metalloproteinase activities. Overall, our findings demonstrate broad levels of antivenom paraspecificity, with in vitro immunological binding and functional inhibition often highly comparable between venoms used to manufacture the antivenoms and those from related species, even in the case of the monovalent antivenom BothroFav. Our findings suggest that the current clinical utility of these antivenoms could possibly be expanded to other parts of Latin America that currently suffer from a lack of specific snakebite therapies.

Snake Venom-specific Phospholipase A2: A Diagnostic Marker for the Management of Snakebite Cases

Background

Snakebites are a common cause of morbidity and mortality, especially in tropical countries. Snakebites in any community are managed based on the clinical features and intravenous administration of antisnake venom (ASV). The administration of ASV is either deficient or given in excess based on clinical decisions and whole blood clotting test results. The present study is designed to analyze the level of snake venom component in the blood of snakebite in association with the clinical features.

Patients and methods

Blood samples were collected from the patients admitted to Karnataka Institute of Medical (KIMS) hospital with a history of snakebite considering the inclusion criteria. Serum was collected from the blood of snakebite patients before and after ASV and used to assess the level of venom-specific phospholipase A2 (PLA2) enzyme using the enzyme-linked immunosorbent assay (ELISA) method.

Results

Quantitative ELISA results revealed that the snake venom-specific PLA2 in the victim's blood was in the range of 0.3-1.27 mg/mL before the administration of ASV. However, the concentration of PLA2 after 24 hours of ASV administration was decreased in most of the patients. Further, it was observed that envenomation complications were directly proportional to the amount of snake venom-specific PLA2 found in the blood of the snakebite patient.

Conclusion

The study concludes that snake venom-specific PLA2 in the blood of snakebite patients could be used as a reliable venom marker, which helps in determination of appropriate ASV dosage in snakebite patients.

How to cite this article

Kaulgud RS, Hasan T, Vanti GL, Veeresh S, Uppar AP, Kurjogi MM. Snake Venom-specific Phospholipase A2: A Diagnostic Marker for the Management of Snakebite Cases. Indian J Crit Care Med 2022;26(12):1259-1266.

Bibliometric Analysis of Literature in Snake Venom-Related Research Worldwide (1933–2022)

Simple Summary

Around the world, snake envenomation poses a serious health risk. Proteins with pharmacological effects are present in snake venom. Recent studies elaborate snake venom and its potential application, including as a cancer drug and antibacterial substances. Our study aimed to analyze the global profile of the literature in snake venom research from documents indexed in the Scopus database between 1933 and 2022. In total, 2999 documents were published with Brazil showing the highest productivity. Antivenom, proteomics, and transcriptomics are emerging as hot topics on a global scale. The present study offers a distinctive overview of snake venom research conducted worldwide.

Abstract

Snake envenomation is a severe economic and health concern affecting countries worldwide. Snake venom carries a wide variety of small peptides and proteins with various immunological and pharmacological properties. A few key research areas related to snake venom, including its applications in treating cancer and eradicating antibiotic-resistant bacteria, have been gaining significant attention in recent years. The goal of the current study was to analyze the global profile of literature in snake venom research. This study presents a bibliometric review of snake venom-related research documents indexed in the Scopus database between 1933 and 2022. The overall number of documents published on a global scale was 2999, with an average annual production of 34 documents. Brazil produced the highest number of documents (n = 729), followed by the United States (n = 548), Australia (n = 240), and Costa Rica (n = 235). Since 1963, the number of publications has been steadily increasing globally. At a worldwide level, antivenom, proteomics, and transcriptomics are growing hot issues for research in this field. The current research provides a unique overview of snake venom research at global level from 1933 through 2022, and it may be beneficial in guiding future research.

Venom of the Annulated Sea Snake Hydrophis cyanocinctus: A Biochemically Simple but Genetically Complex Weapon

Given that the venom system in sea snakes has a role in enhancing their secondary adaption to the marine environment, it follows that elucidating the diversity and function of venom toxins will help to understand the adaptive radiation of sea snakes. We performed proteomic and de novo NGS analyses to explore the diversity of venom toxins in the annulated sea snake (Hydrophis cyanocinctus) and estimated the adaptive molecular evolution of the toxin-coding unigenes and the toxicity of the major components. We found three-finger toxins (3-FTxs), phospholipase A2 (PLA2) and cysteine-rich secretory protein (CRISP) in the venom proteome and 59 toxin-coding unigenes belonging to 24 protein families in the venom-gland transcriptome; 3-FTx and PLA2 were the most abundant families. Nearly half of the toxin-coding unigenes had undergone positive selection. The short- (i.p. 0.09 μg/g) and long-chain neurotoxin (i.p. 0.14 μg/g) presented fairly high toxicity, whereas both basic and acidic PLA2s expressed low toxicity. The toxicity of H. cyanocinctus venom was largely determined by the 3-FTxs. Our data show the venom is used by H. cyanocinctus as a biochemically simple but genetically complex weapon and venom evolution in H. cyanocinctus is presumably driven by natural selection to deal with fast-moving prey and enemies in the marine environment.

Venom-gland transcriptomic, venomic, and antivenomic profiles of the spine-bellied sea snake (Hydrophis curtus) from the South China Sea

Background

A comprehensive evaluation of the -omic profiles of venom is important for understanding the potential function and evolution of snake venom. Here, we conducted an integrated multi-omics-analysis to unveil the venom-transcriptomic and venomic profiles in a same group of spine-bellied sea snakes (Hydrophis curtus) from the South China Sea, where the snake is a widespread species and might generate regionally-specific venom potentially harmful to human activities. The capacity of two heterologous antivenoms to immunocapture the H. curtus venom was determined for an in-depth evaluation of their rationality in treatment of H. curtus envenomation. In addition, a phylogenetic analysis by maximum likelihood was used to detect the adaptive molecular evolution of full-length toxin-coding unigenes.

Results

A total of 90,909,384 pairs of clean reads were generated via Illumina sequencing from a pooled cDNA library of six specimens, and yielding 148,121 unigenes through de novo assembly. Sequence similarity searching harvested 63,845 valid annotations, including 63,789 non-toxin-coding and 56 toxin-coding unigenes belonging to 22 protein families. Three protein families, three-finger toxins (3-FTx), phospholipase A₂ (PLA₂), and cysteine-rich secretory protein, were detected in the venom proteome. 3-FTx (27.15% in the transcriptome/41.94% in the proteome) and PLA₂ (59.71%/49.36%) were identified as the most abundant families in the venom-gland transcriptome and venom proteome. In addition, 24 unigenes from 11 protein families were shown to have experienced positive selection in their evolutionary history, whereas four were relatively conserved throughout evolution. Commercial Naja atra antivenom exhibited a stronger capacity than Bungarus multicinctus antivenom to immunocapture H. curtus venom components, especially short neurotoxins, with the capacity of both antivenoms to immunocapture short neurotoxins being weaker than that for PLA₂s.

Conclusions

Our study clarified the venom-gland transcriptomic and venomic profiles along with the within-group divergence of a H. curtus population from the South China Sea. Adaptive evolution of most venom components driven by natural selection appeared to occur rapidly during evolutionary history. Notably, the utility of commercial N. atra and B. multicinctus antivenoms against H. curtus toxins was not comprehensive; thus, the development of species-specific antivenom is urgently needed.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07824-7.

Differentiation of Central Slovenian and Moscow populations of Rana temporaria frogs using peptide biomarkers of temporins family

Skin secretion represents the only means of defense for the majority of frog species. That phenomenon is based on the fact that the main components of the secretion are peptides demonstrating greatly varying types of bioactivity. They fulfill regulatory functions, fight microorganisms and may be even helpful against predators. These peptides are considered to be rather promising pharmaceuticals of future generation as according to the present knowledge microorganisms are unlikely to develop resistance to them. Mass spectrometry sequencing of these peptides is the most efficient first step of their study providing reliably their primary structures, i.e., amino acids sequence and S-S bond motif. Besides discovering new bioactive peptides, mass spectrometry appears to be an efficient tool of taxonomy studies, allowing for distinguishing not only between closely related species, but also between populations of the same species. Application of several tandem mass spectrometry tools (CID, HCD, ETD, EThcD) available with Orbitrap mass analyzer allowed us to obtain full sequence of about 60 peptides in the secretion of Slovenian population of brown ranid frog Rana temporaria. The problem of sequence inside C-terminal cycle formed by two Cys and differentiation of isomeric Leu and Ile residues was done in top-down mode without any derivatization steps. Besides general biomarkers of Rana temporaria species, Central Slovenian population of Rana temporaria demonstrates six novel temporins and one brevinin 1, which may be treated as biomarkers of that population.

Kinetic analysis of effects of temperature and time on the regulation of venom expression in Bungarus multicinctus

Venom gland is a highly efficient venom production system to maintain their predatory arsenal. Venom toxins mRNA has been shown to increase abruptly in snake after venom expenditure, while the dynamics of venom accumulation during synthesis are poorly understood. Here, PacBio long-read sequencing, Illumina RNA sequencing (RNA-seq), and label-free proteome quantification were used to investigate the composition landscape and time- and temperature-dependent dynamics changes of the Bungarus multicinctus venom gland system. Transcriptome data (19.5223 Gb) from six adult B. multicinctus tissues were sequenced using PacBio RS II to generate a reference assembly, and average 7.28 Gb of clean RNA-seq data was obtained from venom glands by Illumina sequencing. Differentially expressed genes (DEGs) mainly were protein processing rather than venom toxins. Post-translational modifications provided the evidence of the significantly different proportions of toxins in the venom proteome with the changing of replenishment time and temperature, but constant of venom toxins mRNA in the venom gland transcriptome. Dynamic of toxins and genes involved in venom gland contraction suggesting the formation of the mature venom gland system would take at least 9 days. In addition, 59 toxin processing genes were identified, peptidylprolyl isomerase B of which underwent positive selection in Toxicofera. These results provide a reference for determining the extraction time of venom, production of polyclonal and monoclonal antibody for precise treatment plans of venomous snakebites, and construction of an in vitro synthesis system for snake venom protein.

A Comprehensive Multi-Omic Approach Reveals a Relatively Simple Venom in a Diet Generalist, the Northern Short-Tailed Shrew, Blarina brevicauda

Animals that use venom to feed on a wide diversity of prey may evolve a complex mixture of toxins to target a variety of physiological processes and prey-defense mechanisms. Blarina brevicauda, the northern short-tailed shrew, is one of few venomous mammals, and is also known to eat evolutionarily divergent prey. Despite their complex diet, earlier proteomic and transcriptomic studies of this shrew's venom have only identified two venom proteins. Here, we investigated with comprehensive molecular approaches whether B. brevicauda venom is more complex than previously understood. We generated de novo assemblies of a B. brevicauda genome and submaxillary-gland transcriptome, as well as sequenced the salivary proteome. Our findings show that B. brevicauda's venom composition is simple relative to their broad diet and is likely limited to seven proteins from six gene families. Additionally, we explored expression levels and rate of evolution of these venom genes and the origins of key duplications that led to toxin neofunctionalization. We also found three proteins that may be involved in endogenous self-defense. The possible synergism of the toxins suggests that vertebrate prey may be the main target of the venom. Further functional assays for all venom proteins on both vertebrate and invertebrate prey would provide further insight into the ecological relevance of venom in this species.

Bothrops atrox, the most important snake involved in human envenomings in the amazon: How venomics contributes to the knowledge of snake biology and clinical toxinology

Bothrops atrox snakes are mostly endemic of the Amazon rainforest and is certainly the South American pit viper responsible for most of the snakebites in the region. The composition of B. atrox venom is significantly known and has been used to trace the relevance of the venom phenotype for snake biology and for the impacts in the clinics of human patients involved in accidents by B. atrox. However, in spite of the wide distribution and the great medical relevance of B. atrox snakes, B. atrox taxonomy is not fully resolved and the impacts of the lack of taxonomic resolution on the studies focused on venom or envenoming are currently unknown. B. atrox venom presents different degrees of compositional variability and is generally coagulotoxic, inducing systemic hematological disturbances and local tissue damage in snakebite patients. Antivenoms are the effective therapy for attenuating the clinical signs. This review brings a comprehensive discussion of the literature concerning B. atrox snakes encompassing from snake taxonomy, diet and venom composition, towards clinical aspects of snakebite patients and efficacy of the antivenoms. This discussion is highly supported by the contributions that venomics and antivenomics added for the advancement of knowledge of B. atrox snakes, their venoms and the treatment of accidents they evoke.

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

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

A New Family of Diverse Skin Peptides from the Microhylid Frog Genus Phrynomantis

A wide range of frogs produce skin poisons composed of bioactive peptides for defence against pathogens, parasites and predators. While several frog families have been thoroughly screened for skin-secreted peptides, others, like the Microhylidae, have remained mostly unexplored. Previous studies of microhylids found no evidence of peptide secretion, suggesting that this defence adaptation was evolutionarily lost. We conducted transcriptome analyses of the skins of Phrynomantis bifasciatus and Phrynomantis microps, two African microhylid species long suspected to be poisonous. Our analyses reveal 17 evolutionary related transcripts that diversified from to those of cytolytic peptides found in other frog families. The 19 peptides predicted to be processed from these transcripts, named phrynomantins, show a striking structural diversity that is distinct from any previously identified frog skin peptide. Functional analyses of five phrynomantins confirm the loss of a cytolytic function and the absence of insecticidal or proinflammatory activity, suggesting that they represent an evolutionary transition to a new, yet unknown function. Our study shows that peptides have been retained in the defence poison of at least one microhylid lineage and encourages research on similarly understudied taxa to further elucidate the diversity and evolution of skin defence molecules.

Trait differentiation and modular toxin expression in palm-pitvipers

Background

Modularity is the tendency for systems to organize into semi-independent units and can be a key to the evolution and diversification of complex biological systems. Snake venoms are highly variable modular systems that exhibit extreme diversification even across very short time scales. One well-studied venom phenotype dichotomy is a trade-off between neurotoxicity versus hemotoxicity that occurs through the high expression of a heterodimeric neurotoxic phospholipase A₂ (PLA₂) or snake venom metalloproteinases (SVMPs). We tested whether the variation in these venom phenotypes could occur via variation in regulatory sub-modules through comparative venom gland transcriptomics of representative Black-Speckled Palm-Pitvipers (Bothriechis nigroviridis) and Talamancan Palm-Pitvipers (B. nubestris).

Results

We assembled 1517 coding sequences, including 43 toxins for B. nigroviridis and 1787 coding sequences including 42 toxins for B. nubestris. The venom gland transcriptomes were extremely divergent between these two species with one B. nigroviridis exhibiting a primarily neurotoxic pattern of expression, both B. nubestris expressing primarily hemorrhagic toxins, and a second B. nigroviridis exhibiting a mixed expression phenotype. Weighted gene coexpression analyses identified six submodules of transcript expression variation, one of which was highly associated with SVMPs and a second which contained both subunits of the neurotoxic PLA₂ complex. The sub-module association of these toxins suggest common regulatory pathways underlie the variation in their expression and is consistent with known patterns of inheritance of similar haplotypes in other species. We also find evidence that module associated toxin families show fewer gene duplications and transcript losses between species, but module association did not appear to affect sequence diversification.

Conclusion

Sub-modular regulation of expression likely contributes to the diversification of venom phenotypes within and among species and underscores the role of modularity in facilitating rapid evolution of complex traits.

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.

Exploring the Diversity and Novelty of Toxin Genes in Naja sumatrana, the Equatorial Spitting Cobra from Malaysia through De Novo Venom-Gland Transcriptomics

The equatorial spitting cobra, Naja sumatrana, is a distinct species of medically important venomous snakes, listed as WHO Category 1 in Southeast Asia. The diversity of its venom genes has not been comprehensively examined, although a few toxin sequences annotated to Naja sputatrix were reported previously through cloning studies. To investigate this species venom genes’ diversity, de novo venom-gland transcriptomics of N. sumatrana from West Malaysia was conducted using next-generation sequencing technology. Genes encoding toxins represented only 60 of the 55,396 transcripts, but were highly expressed, contributing to 79.22% of total gene expression (by total FPKM) in the venom-glands. The toxin transcripts belong to 21 families, and 29 transcripts were further identified as full-length. Three-finger toxins (3FTx) composed of long, short, and non-conventional groups, constituted the majority of toxin transcripts (91.11% of total toxin FPKM), followed by phospholipase A₂ (PLA₂, 7.42%)—which are putatively pro-inflammatory and cytotoxic. The remaining transcripts in the 19 families were expressed at extremely low levels. Presumably, these toxins were associated with ancillary functions. Our findings unveil the diverse toxin genes unique to N. sumatrana, and provide insights into the pathophysiology of N. sumatrana envenoming.

Snake venom proteome and immuno-profiling of the hundred-pace viper, Deinagkistrodon acutus, in Taiwan

Deinagkistrodon acutus, also known as the hundred-pace viper or Chinese moccasin, is a clinically significant venomous snake in Taiwan. To address the lack of knowledge on the venom proteome of D. acutus, the venom composition was studied by a bottom-up proteomic approach combining reverse phase high-performance liquid chromatography, SDS-PAGE, and LC-MS/MS analysis. The immunoreactivity and cross-reactivity of Taiwanese freeze-dried D. acutus antivenom (DA-AV) and hemorrhagic antivenom (FH-AV) were investigated, as well. The proteomic analysis revealed the presence of 29 distinct proteins from D. acutus venom belonging to 8 snake venom protein families. Snake venom metalloproteinase (SVMP, 46.86%), C-type lectin (CLEC, 37.59%), phospholipase A₂ (PLA₂, 7.33%) and snake venom serine protease (SVSP, 6.62%) were the most abundant proteins. In addition to DA-AV, FH-AV also showed a profile of broad immunorecognition toward the venom of D. acutus. Remarkably, both antivenoms specifically reacted with the HPLC fractions containing SVMPs, and the titer was 5-10 times higher than fractions of other components. This information helps us to deeply understand the pathophysiology of D. acutus envenomation and guide us to development of more effective antivenom for clinical treatment.

The Harderian gland transcriptomes of Caraiba andreae, Cubophis cantherigerus and Tretanorhinus variabilis, three colubroid snakes from Cuba

The Harderian gland is a cephalic structure, widely distributed among vertebrates. In snakes, the Harderian gland is anatomically connected to the vomeronasal organ via the nasolacrimal duct, and in some species can be larger than the eyes. The function of the Harderian gland remains elusive, but it has been proposed to play a role in the production of saliva, pheromones, thermoregulatory lipids and growth factors, among others. Here, we have profiled the transcriptomes of the Harderian glands of three non-front-fanged colubroid snakes from Cuba: Caraiba andreae (Cuban Lesser Racer); Cubophis cantherigerus (Cuban Racer); and Tretanorhinus variabilis (Caribbean Water Snake), using Illumina HiSeq2000 100 bp paired-end. In addition to ribosomal and non-characterized proteins, the most abundant transcripts encode putative transport/binding, lipocalin/lipocalin-like, and bactericidal/permeability-increasing-like proteins. Transcripts coding for putative canonical toxins described in venomous snakes were also identified. This transcriptional profile suggests a more complex function than previously recognized for this enigmatic organ.

Integrated Venomics and Venom Gland Transcriptome Analysis of Juvenile and Adult Mexican Rattlesnakes Crotalus simus, C. tzabcan, and C. culminatus Revealed miRNA-modulated Ontogenetic Shifts

Adult rattlesnakes within genus Crotalus express one of two distinct venom phenotypes, type I (hemorrhagic) and type II (neurotoxic). In Costa Rican Central American rattlesnake, ontogenetic changes in the concentration of miRNAs modulate venom type II to type I transition. Venomics and venom gland transcriptome analyses showed that adult C. simus and C. tzabcan expressed intermediate patterns between type II and type I venoms, whereas C. culminatus had a canonical type I venom. Neonate/juvenile and adult Mexican rattlesnakes showed notable inter- and intraspecific variability in the number, type, abundance and ontogenetic shifts of the transcriptional and translational venom gland activities. These results support a role for miRNAs in the ontogenetic venom compositional changes in the three congeneric Mexican rattlesnakes. It is worth noting the finding of dual-action miRNAs, which silence the translation of neurotoxic heterodimeric PLA₂ crotoxin and acidic PLA₂ mRNAs while simultaneously up-regulating SVMP-targeting mRNAs. Dual transcriptional regulation potentially explains the existence of mutually exclusive crotoxin-rich (type-II) and SVMP-rich (type-I) venom phenotypic dichotomy among rattlesnakes. Our results support the hypothesis that alterations of the distribution of miRNAs, modulating the translational activity of venom gland toxin-encoding mRNAs in response to an external cue, may contribute to the mechanism generating adaptive venom variability.

Comparative venom gland transcriptomics of Naja kaouthia (monocled cobra) from Malaysia and Thailand: elucidating geographical venom variation and insights into sequence novelty

Background

The monocled cobra (Naja kaouthia) is a medically important venomous snake in Southeast Asia. Its venom has been shown to vary geographically in relation to venom composition and neurotoxic activity, indicating vast diversity of the toxin genes within the species. To investigate the polygenic trait of the venom and its locale-specific variation, we profiled and compared the venom gland transcriptomes of N. kaouthia from Malaysia (NK-M) and Thailand (NK-T) applying next-generation sequencing (NGS) technology.

Methods

The transcriptomes were sequenced on the Illumina HiSeq platform, assembled and followed by transcript clustering and annotations for gene expression and function. Pairwise or multiple sequence alignments were conducted on the toxin genes expressed. Substitution rates were studied for the major toxins co-expressed in NK-M and NK-T.

Results and discussion

The toxin transcripts showed high redundancy (41–82% of the total mRNA expression) and comprised 23 gene families expressed in NK-M and NK-T, respectively (22 gene families were co-expressed). Among the venom genes, three-finger toxins (3FTxs) predominated in the expression, with multiple sequences noted. Comparative analysis and selection study revealed that 3FTxs are genetically conserved between the geographical specimens whilst demonstrating distinct differential expression patterns, implying gene up-regulation for selected principal toxins, or alternatively, enhanced transcript degradation or lack of transcription of certain traits. One of the striking features that elucidates the inter-geographical venom variation is the up-regulation of α-neurotoxins (constitutes ∼80.0% of toxin’s fragments per kilobase of exon model per million mapped reads (FPKM)), particularly the long-chain α-elapitoxin-Nk2a (48.3%) in NK-T but only 1.7% was noted in NK-M. Instead, short neurotoxin isoforms were up-regulated in NK-M (46.4%). Another distinct transcriptional pattern observed is the exclusively and abundantly expressed cytotoxin CTX-3 in NK-T. The findings suggested correlation with the geographical variation in proteome and toxicity of the venom, and support the call for optimising antivenom production and use in the region. Besides, the current study uncovered full and partial sequences of numerous toxin genes from N. kaouthia which have not been reported hitherto; these include N. kaouthia-specific l-amino acid oxidase (LAAO), snake venom serine protease (SVSP), cystatin, acetylcholinesterase (AChE), hyaluronidase (HYA), waprin, phospholipase B (PLB), aminopeptidase (AP), neprilysin, etc. Taken together, the findings further enrich the snake toxin database and provide deeper insights into the genetic diversity of cobra venom toxins.

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.

Proteomic analysis of venom variability and ontogeny across the arboreal palm-pitvipers (genus Bothriechis)

Bothriechis is a genus of eleven currently recognized slender and arboreal venomous snakes, commonly called palm-pitvipers that range from southern Mexico to northern South America. Despite dietary studies suggesting that palm-pitvipers are generalists with an ontogenetic shift toward endothermic prey, venom proteomic analyses have revealed remarkable divergence between the venoms of the Costa Rican species, B. lateralis, B. schlegelii, B. supraciliaris, and B. nigroviridis. To achieve a more complete picture of the venomic landscape across Bothriechis, the venom proteomes of biodiversity of the northern Middle American highland palm-pitvipers, B. thalassinus, B. aurifer, and B. bicolor from Guatemala, B. marchi from Honduras, and neonate Costa Rican B. lateralis and B. schlegelii, were investigated. B. thalassinus and B. aurifer venoms are comprised by similar toxin arsenals dominated by SVMPs (33-39% of the venom proteome), CTLs (11-16%), BPP-like molecules (10-13%), and CRISPs (5-10%), and are characterized by the absence of PLA₂ proteins. Conversely, the predominant (35%) components of B. bicolor are D49-PLA₂ molecules. The venom proteome of B. marchi is similar to B. aurifer and B. thalassinus in that it is rich in SVMPs and BPPs, but also contains appreciable amounts (14.3%) of PLA₂s. The major toxin family found in the venoms of both neonate B. lateralis and B. schlegelii, is serine proteinase (SVSP), comprising about 20% of their toxin arsenals. The venom of neonate B. schlegelii is the only palm-pitviper venom where relative high amounts of Kunitz-type (6.3%) and γPLA₂ (5.2%) inhibitors have been identified. Despite notable differences between their proteomes, neonate venoms are more similar to each other than to adults of their respective species. However, the ontogenetic changes taking place in the venom of B. lateralis strongly differ from those that occur in the venom of B. schlegelii. Thus, the ontogenetic change in B. lateralis produces a SVMP-rich venom, whereas in B. schlegelii the age-dependent compositional shift generates a PLA₂-rich venom. Overall, genus-wide venomics illustrate the high evolvability of palm-pitviper venoms. The integration of the pattern of venom variation across Bothriechis into a phylogenetic and biogeographic framework may lay the foundation for assessing, in future studies, the evolutionary path that led to the present-day variability of the venoms of palm-pitvipers.

SIGNIFICANCE

Bothriechis represents a monophyletic basal genus of eleven arboreal palm-pitvipers that range from southern Mexico to northern South America. Despite palm-pitvipers' putative status as diet generalists, previous proteomic analyses have revealed remarkable divergence between the venoms of Costa Rican species, B. lateralis, B. schlegelii, B. supraciliaris, and B. nigroviridis. Our current proteomic study of Guatemalan species, B. thalassinus, B. aurifer, and B. bicolor, Honduran B. marchi, and neonate B. lateralis and B. schlegelii from Costa Rica was undertaken to deepen our understanding of the evolutionary pattern of venom proteome diversity across Bothriechis. Ancestral characters are often, but not always, preserved in an organism's development. Venoms of neonate B. lateralis and B. schlegelii are more similar to each other than to adults of their respective species, suggesting that the high evolvability of palm-pitviper venoms may represent an inherent feature of Bothriechis common ancestor. Our genus-wide data identified four nodes of venom phenotype differentiation across the phylogeny of Bothriechis. Integrated into a phylogenetic and biogeographic framework, the pattern of venom variation across Bothriechis may lay the groundwork to establish whether divergence was driven by selection for efficient resource exploitation in arboreal 'islands', thereby contributing to the ecological speciation of the genus.

RNA sequence analyses of r-Moj-DM treated cells: TXNIP is required to induce apoptosis of SK-Mel-28

RNA sequencing of untreated and r-Moj-DM treated SK-Mel-28 cells was performed after 6 h, to begin unraveling the apoptotic pathway induced by r-Moj-DM. Bioinformatic analyses of RNA sequencing data yielded 40 genes that were differentially expressed. Nine genes were upregulated and 31 were downregulated. qRT-PCR was used to validate differential expression of 13 genes with known survival or apoptotic-inducing activities. Expression of BNiP3, IGFBP3, PTPSF, Prune 2, TGF-ß, and TXNIP were compared from cells treated with r-Moj-DN (a strong apoptotic inducer) or r-Moj-DA (a non-apoptotic inducer) for 1 h, 2 h, 4 h, and 6 h after treatment. Our results demonstrate that significant differences in expression are only detected after 4 h of treatment. In addition, expression of TXNIP (an apoptotic inducer) remains elevated at 4 h and 6 h only in r-Moj-DN treated cells. Based on the consistency of elevated TXNIP expression, we further studied TXNIP as a novel target of disintegrin activation. Confocal microscopy of anti-TXNIP stained SK-Mel-28 cells suggests nuclear localization of TXNIP after r-Moj-DM treatment. A stable TXNIP knockdown SK-Mel-28 cell line was produced to test TXNIP' role in the apoptotic induction by r-Moj-DM. High cell viability (74.3% ±9.1) was obtained after r-Moj-DM treatment of TXNIP knocked down SK-Mel-28 cells, compared to 34% ±0.187 for untransduced cells. These results suggest that TXNIP is required early in the apoptotic-inducing pathway resulting from r-Moj-DM binding to the αv integrin subunit.

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.

Full-Length Venom Protein cDNA Sequences from Venom-Derived mRNA: Exploring Compositional Variation and Adaptive Multigene Evolution

Envenomation of humans by snakes is a complex and continuously evolving medical emergency, and treatment is made that much more difficult by the diverse biochemical composition of many venoms. Venomous snakes and their venoms also provide models for the study of molecular evolutionary processes leading to adaptation and genotype-phenotype relationships. To compare venom complexity and protein sequences, venom gland transcriptomes are assembled, which usually requires the sacrifice of snakes for tissue. However, toxin transcripts are also present in venoms, offering the possibility of obtaining cDNA sequences directly from venom. This study provides evidence that unknown full-length venom protein transcripts can be obtained from the venoms of multiple species from all major venomous snake families. These unknown venom protein cDNAs are obtained by the use of primers designed from conserved signal peptide sequences within each venom protein superfamily. This technique was used to assemble a partial venom gland transcriptome for the Middle American Rattlesnake (Crotalus simus tzabcan) by amplifying sequences for phospholipases A₂, serine proteases, C-lectins, and metalloproteinases from within venom. Phospholipase A₂ sequences were also recovered from the venoms of several rattlesnakes and an elapid snake (Pseudechis porphyriacus), and three-finger toxin sequences were recovered from multiple rear-fanged snake species, demonstrating that the three major clades of advanced snakes (Elapidae, Viperidae, Colubridae) have stable mRNA present in their venoms. These cDNA sequences from venom were then used to explore potential activities derived from protein sequence similarities and evolutionary histories within these large multigene superfamilies. Venom-derived sequences can also be used to aid in characterizing venoms that lack proteomic profiles and identify sequence characteristics indicating specific envenomation profiles. This approach, requiring only venom, provides access to cDNA sequences in the absence of living specimens, even from commercial venom sources, to evaluate important regional differences in venom composition and to study snake venom protein evolution.

This work demonstrates that full-length venom protein messenger RNAs are present in secreted venoms and can be used to acquire full-length protein sequences of toxins from both front-fanged (Elapidae, Viperidae) and rear-fanged (Colubridae) snake venoms, eliminating the need to use venom glands. Full-length transcripts were obtained from venom samples that were fresh, newly lyophilized, old, field desiccated or commercially prepared, representing a significant advance over previous attempts which produced only partial sequence transcripts. Transcripts for all major venom protein families (metalloproteinases, serine proteases, C-type lectins, phospholipases A₂ and three-finger toxins) responsible for clinically significant snakebite symptoms were obtained from venoms. These sequences aid in the identification and characterization of venom proteome profiles, allowing for the identification of peptide sequences, specific isoforms, and novel venom proteins. The application of this technique will help to provide venom protein sequences for many snake species, including understudied rear-fanged snakes. Venom protein transcripts offer important insights into potential snakebite envenomation profiles and the molecular evolution of venom protein multigene families. By requiring only venom to obtain venom protein cDNAs, the approach detailed here will provide access to cDNA-based protein sequences from commercial and other venom sources, facilitating study of snake venom protein composition and evolution.

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.

Proteomic characterization and comparison of venoms from two elapid snakes (Bungarus multicinctus and Naja atra) from China

Bungarus multicinctus (many-banded krait) and Naja atra (Chinese cobra) are widely distributed and medically important venomous snakes in China; however, their venom proteomic profiles have not been fully compared. Here, we fractionated crude venoms and analyzed them using a combination of proteomic techniques. Three-finger toxins (3-FTx) and phospholipase A2 (PLA2) were most abundant in both species, respectively accounting for 32.6% and 66.4% of total B. multicinctus venom, and 84.3% and 12.2% of total N. atra venom. Venoms from these two species contained one common protein family and six less abundant species-specific protein families. The proteomic profiles of B. multicinctus and N. atra venoms and analysis of toxicological activity in mice suggested that 3-FTx and PLA2 are the major contributors to clinical symptoms caused by envenomation. The venoms differed in enzymatic activity, likely the result of inter-specific variation in the amount of related venom components. Antivenomics assessment revealed that a small number of venom components (3-FTxs and PLA2s in B. multicinctus, and 3-FTxs in N. atra) could not be immunocaptured completely, suggesting that we should pay attention to enhancing the immune response of these components in designing commercial antivenoms for B. multicinctus and N. atra.

BIOLOGICAL SIGNIFICANCE

The proteomic profiles of venoms from two medically important snake species - B. multicinctus and N. atra - have been explored. Quantitative and qualitative differences are evident in both venoms when proteomic profiles and transcriptomic results are compared; this is a reminder that combined approaches are needed to explore the precise composition of snake venom. Two protein families (3-FTx and PLA2) of high abundance in these snake venoms are major players in the biochemical and pharmacological effects of envenomation. Elucidation of the proteomic profiles of these snake venoms is helpful in understanding composition-function relationships and will facilitate the clinical application of antivenoms.

Ecological venomics: How genomics, transcriptomics and proteomics can shed new light on the ecology and evolution of venom

Animal venom is a complex cocktail of bioactive chemicals that traditionally drew interest mostly from biochemists and pharmacologists. However, in recent years the evolutionary and ecological importance of venom is realized as this trait has direct and strong influence on interactions between species. Moreover, venom content can be modulated by environmental factors. Like many other fields of biology, venom research has been revolutionized in recent years by the introduction of systems biology approaches, i.e., genomics, transcriptomics and proteomics. The employment of these methods in venom research is known as 'venomics'. In this review we describe the history and recent advancements of venomics and discuss how they are employed in studying venom in general and in particular in the context of evolutionary ecology. We also discuss the pitfalls and challenges of venomics and what the future may hold for this emerging scientific field.

Venom-gland transcriptome and venom proteome of the Malaysian king cobra (Ophiophagus hannah)

Background

The king cobra (Ophiophagus hannah) is widely distributed throughout many parts of Asia. This study aims to investigate the complexity of Malaysian Ophiophagus hannah (MOh) venom for a better understanding of king cobra venom variation and its envenoming pathophysiology. The venom gland transcriptome was investigated using the Illumina HiSeq™ platform, while the venom proteome was profiled by 1D-SDS-PAGE-nano-ESI-LCMS/MS.

Results

Transcriptomic results reveal high redundancy of toxin transcripts (3357.36 FPKM/transcript) despite small cluster numbers, implying gene duplication and diversification within restricted protein families. Among the 23 toxin families identified, three-finger toxins (3FTxs) and snake-venom metalloproteases (SVMPs) have the most diverse isoforms. These 2 toxin families are also the most abundantly transcribed, followed in descending order by phospholipases A₂ (PLA₂s), cysteine-rich secretory proteins (CRISPs), Kunitz-type inhibitors (KUNs), and L-amino acid oxidases (LAAOs). Seventeen toxin families exhibited low mRNA expression, including hyaluronidase, DPP-IV and 5’-nucleotidase that were not previously reported in the venom-gland transcriptome of a Balinese O. hannah. On the other hand, the MOh proteome includes 3FTxs, the most abundantly expressed proteins in the venom (43 % toxin sbundance). Within this toxin family, there are 6 long-chain, 5 short-chain and 2 non-conventional 3FTx. Neurotoxins comprise the major 3FTxs in the MOh venom, consistent with rapid neuromuscular paralysis reported in systemic envenoming. The presence of toxic enzymes such as LAAOs, SVMPs and PLA₂ would explain tissue inflammation and necrotising destruction in local envenoming. Dissimilarities in the subtypes and sequences between the neurotoxins of MOh and Naja kaouthia (monocled cobra) are in agreement with the poor cross-neutralization activity of N. kaouthia antivenom used against MOh venom. Besides, the presence of cobra venom factor, nerve growth factors, phosphodiesterase, 5’-nucleotidase, and DPP-IV in the venom proteome suggests its probable hypotensive action in subduing prey.

Conclusion

This study reports the diversity and abundance of toxins in the venom of the Malaysian king cobra (MOh). The results correlate with the pathophysiological actions of MOh venom, and dispute the use of Naja cobra antivenoms to treat MOh envenomation. The findings also provide a deeper insight into venom variations due to geography, which is crucial for the development of a useful pan-regional antivenom.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1828-2) contains supplementary material, which is available to authorized users.

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.

Whole Transcriptome of the Venom Gland from Urodacus yaschenkoi Scorpion

Australian scorpion venoms have been poorly studied, probably because they do not pose an evident threat to humans. In addition, the continent has other medically important venomous animals capable of causing serious health problems. Urodacus yaschenkoi belongs to the most widely distributed family of Australian scorpions (Urodacidae) and it is found all over the continent, making it a useful model system for studying venom composition and evolution. This communication reports the whole set of mRNA transcripts produced by the venom gland. U. yaschenkoi venom is as complex as its overseas counterparts. These transcripts certainly code for several components similar to known scorpion venom components, such as: alpha-KTxs, beta-KTxs, calcins, protease inhibitors, antimicrobial peptides, sodium-channel toxins, toxin-like peptides, allergens, La1-like, hyaluronidases, ribosomal proteins, proteasome components and proteins related to cellular processes. A comparison with the venom gland transcriptome of Centruroides noxius (Buthidae) showed that these two scorpions have similar components related to biological processes, although important differences occur among the venom toxins. In contrast, a comparison with sequences reported for Urodacus manicatus revealed that these two Urodacidae species possess the same subfamily of scorpion toxins. A comparison with sequences of an U. yaschenkoi cDNA library previously reported by our group showed that both techniques are reliable for the description of the venom components, but the whole transcriptome generated with Next Generation Sequencing platform provides sequences of all transcripts expressed. Several of which were identified in the proteome, but many more transcripts were identified including uncommon transcripts. The information reported here constitutes a reference for non-Buthidae scorpion venoms, providing a comprehensive view of genes that are involved in venom production. Further, this work identifies new putative bioactive compounds that could be used to seed research into new pharmacological compounds and increase our understanding of the function of different ion channels.

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.

Deciphering the venomic transcriptome of killer-wasp Vespa velutina

Wasp stings have been arising to be a severe public health problem in China in recent years. However, molecular information about lethal or toxic factors in wasp venom is extremely lacking. In this study, we used two pyrosequencing platforms to analyze the transcriptome of Vespa velutina, the most common wasp species native in China. Besides the substantial amount of transcripts encoding for allergens usually regarded as the major lethal factor of wasp sting, a greater abundance of hemostasis-impairing toxins and neurotoxins in the venom of V. velutina were identified, implying that toxic reactions and allergic effects are envenoming strategy for the dangerous outcomes. The pattern of differentially expressed genes before and after venom extraction clearly indicates that the manifestation of V. velutina stings depends on subtle regulations in the metabolic pathway required for toxin recruitment. This comparative analysis offers timely clues for developing clinical treatments for wasp envenoming in China and around the world.

Structural and Functional Diversity of Peptide Toxins from Tarantula Haplopelma hainanum (Ornithoctonus hainana) Venom Revealed by Transcriptomic, Peptidomic, and Patch Clamp Approaches

Spider venom is a complex mixture of bioactive peptides to subdue their prey. Early estimates suggested that over 400 venom peptides are produced per species. In order to investigate the mechanisms responsible for this impressive diversity, transcriptomics based on second generation high throughput sequencing was combined with peptidomic assays to characterize the venom of the tarantula Haplopelma hainanum. The genes expressed in the venom glands were identified, and the bioactivity of their protein products was analyzed using the patch clamp technique. A total of 1,136 potential toxin precursors were identified that clustered into 90 toxin groups, of which 72 were novel. The toxin peptides clustered into 20 cysteine scaffolds that included between 4 and 12 cysteines, and 14 of these groups were newly identified in this spider. Highly abundant toxin peptide transcripts were present and resulted from hypermutation and/or fragment insertion/deletion. In combination with variable post-translational modifications, this genetic variability explained how a limited set of genes can generate hundreds of toxin peptides in venom glands. Furthermore, the intraspecies venom variability illustrated the dynamic nature of spider venom and revealed how complex components work together to generate diverse bioactivities that facilitate adaptation to changing environments, types of prey, and milking regimes in captivity.