Bothrops jararaca venom proteome rearrangement upon neonate to adult transition

Biochemical and toxicological profiles of venoms from an adult female South American bushmaster (Lachesis muta rhombeata) and her offspring

In this work, we compared the biochemical and toxicological profiles of venoms from an adult female specimen of Lachesis muta rhombeata (South American bushmaster) and her seven offspring born in captivity, based on SDS-PAGE, RP-HPLC, enzymatic, coagulant, and hemorrhagic assays. Although adult and juvenile venoms showed comparable SDS-PAGE profiles, juveniles lacked some chromatographic peaks compared with adult venom. Adult venom had higher proteolytic (caseinolytic) activity than juvenile venoms (p < 0.05), but there were no significant inter-venom variations in the esterase, PLA2, phosphodiesterase and L-amino acid oxidase (LAAO) activities, although the latter activity was highly variable among the venoms. Juveniles displayed higher coagulant activity on human plasma, with a minimum coagulant dose ∼42% lower than the adult venom (p < 0.05), but there were no age-related differences in thrombin-like activity. Adult venom was more fibrinogenolytic (based on the rate of fibrinogen chain degradation) and hemorrhagic than juvenile venoms (p < 0.05). The effective dose of Bothrops/Lachesis antivenom (produced by the Instituto Butantan) needed to neutralize the coagulant activity was ∼57% greater for juvenile venoms (p < 0.05), whereas antivenom did not attenuate the thrombin-like activity of juvenile and adult venoms. Antivenom significantly reduced the hemorrhagic activity of adult venom (400 μg/kg, i. d.), but not that of juvenile venoms. Overall, these data indicate a compositional and functional ontogenetic shift in L. m. rhombeata venom.

Proteomics and life-history variability of Endogenous Phospholipases A2 Inhibitors (PLIs) in Bothrops jararaca plasma

In Brazil, the genus Bothrops is responsible for most ophidian accidents. Snake venoms have a wide variety of proteins and peptides exhibiting a broad repertoire of pharmacological and toxic effects that elicit systemic injury and characteristic local effects. The snakes' natural resistance to envenomation caused by the presence of inhibitory compounds on their plasma have been extensively studied. However, the presence of these inhibitors in different developmental stages is yet to be further discussed. The aim of this study was to evaluate the ontogeny of Bothrops jararaca plasma inhibitor composition and, to this end, plasma samples of B. jararaca were obtained from different developmental stages (neonates, youngs, and adults) and sexes (female and male). SDS-PAGE, Western blotting, affinity chromatography, and mass spectrometry were performed to analyze the protein profile and interaction between B. jararaca plasma and venom proteins. In addition, the presence of γBjPLI, a PLA2 inhibitor previously identified and characterized in B. jararaca serum, was confirmed by Western blotting. According to our results, 9-17% of plasma proteins were capable of binding to venom proteins in the three developmental stages. The presence of different endogenous inhibitors and, more specifically, different PLA2 inhibitor (PLI) classes and antihemorrhagic factors were confirmed in specimens of B. jararaca from newborn by mass spectrometry. For the first time, the αPLI and βPLI were detected in B. jararaca plasma, although low or no ontogenetic and sexual correlation were found. The γPLI were more abundant in adult female, than in neonate and young female, but similar to neonate, young and adult male according to the results of mass spectrometry analysis. Our results suggest that there are proteins in the plasma of these animals that can help counteract the effects of self-envenomation from birth.

miRNAs derived from cobra venom exosomes contribute to the cobra envenomation

Currently, there is an increasing amount of evidence indicating that exosomes and the miRNAs they contain are crucial players in various biological processes. However, the role of exosomes and miRNAs in snake venom during the envenomation process remains largely unknown. In this study, fresh venom from Naja atra of different ages (2-month-old, 1-year-old, and 5-year-old) was collected, and exosomes were isolated through ultracentrifugation. The study found that exosomes with inactivated proteins and enzymes can still cause symptoms similar to cobra envenomation, indicating that substances other than proteins and enzymes in exosomes may also play an essential role in cobra envenomation. Furthermore, the expression profiles of isolated exosome miRNAs were analyzed. The study showed that a large number of miRNAs were co-expressed and abundant in cobra venom exosomes (CV-exosomes) of different ages, including miR-2904, which had high expression abundance and specific sequences. The specific miR-2094 derived from CV-exosomes (CV-exo-miR-2904) was overexpressed both in vitro and in vivo. As a result, CV-exo-miR-2904 induced symptoms similar to cobra envenomation in mice and caused liver damage, demonstrating that it plays a crucial role in cobra envenomation. These results reveal that CV-exosomes and the miRNAs they contain play a significant regulatory role in cobra envenomation. Our findings provide new insights for the treatment of cobra bites and the development of snake venom-based medicines.

The Venom Composition of the Snake Tribe Philodryadini: 'Omic' Techniques Reveal Intergeneric Variability among South American Racers

Snakes of the Philodryadini tribe are included in the Dipsadidae family, which is a diverse group of rear-fanged snakes widespread in different ecological conditions, including habitats and diet. However, little is known about the composition and effects of their venoms despite their relevance for understanding the evolution of these snakes or even their impact on the occasional cases of human envenoming. In this study, we integrated venom gland transcriptomics, venom proteomics and functional assays to characterize the venoms from eight species of the Philodryadini tribe, which includes the genus Philodryas, Chlorosoma and Xenoxybelis. The most abundant components identified in the venoms were snake venom metalloproteinases (SVMPs), cysteine-rich secretory proteins (CRISPs), C-type lectins (CTLs), snake endogenous matrix metalloproteinases type 9 (seMMP-9) and snake venom serinoproteinases (SVSPs). These protein families showed a variable expression profile in each genus. SVMPs were the most abundant components in Philodryas, while seMMP-9 and CRISPs were the most expressed in Chlorosoma and Xenoxybelis, respectively. Lineage-specific differences in venom composition were also observed among Philodryas species, whereas P. olfersii presented the highest amount of SVSPs and P. agassizii was the only species to express significant amounts of 3FTx. The variability observed in venom composition was confirmed by the venom functional assays. Philodryas species presented the highest SVMP activity, whereas Chlorosoma species showed higher levels of gelatin activity, which may correlate to the seMMP-9 enzymes. The variability observed in the composition of these venoms may be related to the tribe phylogeny and influenced by their diets. In the presented study, we expanded the set of venomics studies of the Philodryadini tribe, which paves new roads for further studies on the evolution and ecology of Dipsadidae snakes.

Sialic acid-containing glycans play a role in the activity of snake venom proteases

Structural variability is a feature of snake venom proteins, and glycosylation is a post-translational modification that contributes to the diversification of venom proteomes. Studies by our group have shown that Bothrops venoms are distinctly defined by their glycoprotein content, and that most hybrid/complex N-glycans identified in these venoms contain sialic acid. Considering that metalloproteases and serine proteases are abundant components of Bothrops venoms and essential in the envenomation process, and that these enzymes contain several glycosylation sites, the role of sialic acid in venom proteolytic activity was evaluated. Here we show that removal of sialic acid by treatment of nine Bothrops venoms with neuraminidase (i) altered the pattern of gelatinolysis in zymography of most venoms and reduced the gelatinolytic activity of all venoms, (ii) decreased the proteolytic activity of some venoms on fibrinogen and the clotting activity of human plasma of all venoms, and (iii) altered the proteolysis profile of plasma proteins by B. jararaca venom, suggesting that sialic acid may play a role in the interaction of proteases with their protein substrates. In contrast, the profile of venom amidolytic activity on Bz-Arg-pNA did not change after removal of sialic acid, indicating that this monosaccharide is not essential in N-glycans of serine proteases acting on small substrates. In summary, these results expand the knowledge about the variability of the subproteomes of Bothrops venom proteases, and for the first time point to the importance of carbohydrate chains containing sialic acid in the enzymatic activities of venom proteases relevant in human envenomation.

Extracellular Vesicles from Bothrops jararaca Venom Are Diverse in Structure and Protein Composition and Interact with Mammalian Cells

Snake venoms are complex cocktails of non-toxic and toxic molecules that work synergistically for the envenoming outcome. Alongside the immediate consequences, chronic manifestations and long-term sequelae can occur. Recently, extracellular vesicles (EVs) were found in snake venom. EVs mediate cellular communication through long distances, delivering proteins and nucleic acids that modulate the recipient cell's function. However, the biological roles of snake venom EVs, including possible cross-organism communication, are still unknown. This knowledge may expand the understanding of envenoming mechanisms. In the present study, we isolated and characterized the EVs from Bothrops jararaca venom (Bj-EVs), giving insights into their biological roles. Fresh venom was submitted to differential centrifugation, resulting in two EV populations with typical morphology and size range. Several conserved EV markers and a subset of venom related EV markers, represented mainly by processing enzymes, were identified by proteomic analysis. The most abundant protein family observed in Bj-EVs was 5'-nucleotidase, known to be immunosuppressive and a low abundant and ubiquitous toxin in snake venoms. Additionally, we demonstrated that mammalian cells efficiently internalize Bj-EVs. The commercial antibothropic antivenom partially recognizes Bj-EVs and inhibits cellular EV uptake. Based on the proteomic results and the in vitro interaction assays using macrophages and muscle cells, we propose that Bj-EVs may be involved not only in venom production and processing but also in host immune modulation and long-term effects of envenoming.

Profilings of subproteomes of lectin-binding proteins of nine Bothrops venoms reveal variability driven by different glycan types

Snake venom proteomes have long been investigated to explore a multitude of biologically active components that are used for prey capture and defense, and are involved in the pathological effects observed upon mammalian envenomation. Glycosylation is a major protein post-translational modification in venoms and contributes to the diversification of proteomes. We have shown that Bothrops venoms are markedly defined by their content of glycoproteins, and that most N-glycan structures of eight Bothrops venoms contain sialic acid, while bisected N-acetylglucosamine was identified in Bothrops cotiara venom. To further investigate the mechanisms involved in the generation of different venoms by related snakes, here the glycoproteomes of nine Bothrops venoms (Bothrops atrox, B. cotiara, Bothrops erythromelas, Bothrops fonsecai, B. insularis, Bothrops jararaca, Bothrops jararacussu, Bothrops moojeni and Bothrops neuwiedi) were comparatively analyzed by enrichment with three lectins of different specificities, recognizing bisecting N-acetylglucosamine- and sialic acid-containing glycoproteins, and mass spectrometry. The lectin capture strategy generated venom fractions enriched with several glycoproteins, including metalloprotease, serine protease, and L- amino acid oxidase, in addition to various types of low abundant enzymes. The different contents of lectin-enriched proteins underscore novel aspects of the variability of the glycoprotein subproteomes of Bothrops venoms and point to the role of distinct types of glycan chains in generating different venoms by closely related snake species.

Bothrops Jararaca Snake Venom Modulates Key Cancer-Related Proteins in Breast Tumor Cell Lines

Cancer is characterized by the development of abnormal cells that divide in an uncontrolled way and may spread into other tissues where they may infiltrate and destroy normal body tissue. Several previous reports have described biochemical anti-tumorigenic properties of crude snake venom or its components, including their capability of inhibiting cell proliferation and promoting cell death. However, to the best of our knowledge, there is no work describing cancer cell proteomic changes following treatment with snake venoms. In this work we describe the quantitative changes in proteomics of MCF7 and MDA-MB-231 breast tumor cell lines following treatment with Bothrops jararaca snake venom, as well as the functional implications of the proteomic changes. Cell lines were treated with sub-toxic doses at either 0.63 μg/mL (low) or 2.5 μg/mL (high) of B. jararaca venom for 24 h, conditions that cause no cell death per se. Proteomics analysis was conducted on a nano-scale liquid chromatography coupled on-line with mass spectrometry (nLC-MS/MS). More than 1000 proteins were identified and evaluated from each cell line treated with either the low or high dose of the snake venom. Protein profiling upon venom treatment showed differential expression of several proteins related to cancer cell metabolism, immune response, and inflammation. Among the identified proteins we highlight histone H3, SNX3, HEL-S-156an, MTCH2, RPS, MCC2, IGF2BP1, and GSTM3. These data suggest that sub-toxic doses of B. jararaca venom have potential to modulate cancer-development related protein targets in cancer cells. This work illustrates a novel biochemical strategy to identify therapeutic targets against cancer cell growth and survival.

Old World Vipers-A Review about Snake Venom Proteomics of Viperinae and Their Variations

Fine-tuned by millions of years of evolution, snake venoms have frightened but also fascinated humanity and nowadays they constitute potential resources for drug development, therapeutics and antivenoms. The continuous progress of mass spectrometry techniques and latest advances in proteomics workflows enabled toxinologists to decipher venoms by modern omics technologies, so-called ‘venomics’. A tremendous upsurge reporting on snake venom proteomes could be observed. Within this review we focus on the highly venomous and widely distributed subfamily of Viperinae (Serpentes: Viperidae). A detailed public literature database search was performed (2003–2020) and we extensively reviewed all compositional venom studies of the so-called Old-World Vipers. In total, 54 studies resulted in 89 venom proteomes. The Viperinae venoms are dominated by four major, four secondary, six minor and several rare toxin families and peptides, respectively. The multitude of different venomics approaches complicates the comparison of venom composition datasets and therefore we differentiated between non-quantitative and three groups of quantitative workflows. The resulting direct comparisons within these groups show remarkable differences on the intra- and interspecies level across genera with a focus on regional differences. In summary, the present compilation is the first comprehensive up-to-date database on Viperinae venom proteomes and differentiating between analytical methods and workflows.

The distinct N-terminomes of Bothrops jararaca newborn and adult venoms

Using approaches of transcriptomics and proteomics we have shown that the phenotype of Bothrops jararaca venom undergoes a significant rearrangement upon neonate to adult transition. Most regulatory processes in biology are intrinsically related to modifications of protein structure, function, and abundance. However, it is unclear to which extent intrinsic proteolysis affects toxins and snake venom phenotypes upon ontogenesis. Here we assessed the natural N-terminome of Bothrops jararaca newborn and adult venoms and explored the degree of N-terminal protein truncation in ontogenetic-based proteome variation. To this end we applied the Terminal Amine Isotopic Labeling of Substrates (TAILS) technology to characterize venom collected in the presence of proteinase inhibitors. We identified natural N-terminal sequences in the newborn (71) and adult (84) venoms, from which only 37 were common to both. However, truncated toxins were found in higher number in the newborn (212) than in the adult (140) venom. Moreover, sequences N-terminally blocked by pyroglutamic acid were identified in the newborn (55) and adult (49) venoms. Most toxin classes identified by their natural N-terminal sequences showed a similar number of unique peptides in the newborn and adult venoms, however, those of serine proteinases and C-type lectins were more abundant in the adult venom. Truncated sequences from at least ten toxin classes were detected, however the catalytic and cysteine-rich domains of metalloproteinases were the most prone to proteolysis, mainly in the newborn venom. Our results underscore the pervasiveness of truncations in most toxin classes and highlight variable post-translational events in newborn and adult venoms.

Venom Profiling of the Insular Species Bothrops alcatraz: Characterization of Proteome, Glycoproteome, and N-Terminome Using Terminal Amine Isotopic Labeling of Substrates

Bothrops alcatraz, a species endemic to Alcatrazes Islands, is regarded as critically endangered due to its small area of occurrence and the declining quality of its habitat. We recently reported the identification of N-glycans attached to toxins of Bothrops species, showing similar compositions in venoms of the B. jararaca complex (B. jararaca, B. insularis, and B. alcatraz). Here, we characterized B. alcatraz venom using electrophoretic, proteomic, and glycoproteomic approaches. Electrophoresis showed that B. alcatraz venom differs from B. jararaca and B. insularis; however, N-glycan removal revealed similarities between them, indicating that the occupation of N-glycosylation sites contributes to interspecies variability in the B. jararaca complex. Metalloproteinase was the major toxin class identified in the B. alcatraz venom proteome followed by serine proteinase and C-type lectin, and overall, the adult B. alcatraz venom resembles that of B. jararaca juvenile specimens. The comparative glycoproteomic analysis of B. alcatraz venom with B. jararaca and B. insularis indicated that there may be differences in the utilization of N-glycosylation motifs among their different toxin classes. Furthermore, we prospected for the first time the N-terminome of a snake venom using the terminal amine isotopic labeling of substrates (TAILS) approach and report the presence of ∼30% of N-termini corresponding to truncated toxin forms and ∼37% N-terminal sequences blocked by pyroglutamic acid in B. alcatraz venom. These findings underscore a low correlation between venom gland transcriptomes and proteomes and support the view that post-translational processes play a major role in shaping venom phenotypes.

Deadly Proteomes: A Practical Guide to Proteotranscriptomics of Animal Venoms

Animal venoms are renowned for their toxicity, biochemical complexity, and as a source of compounds with potential applications in medicine, agriculture, and industry. Polypeptides underlie much of the pharmacology of animal venoms, and elucidating these arsenals of polypeptide toxins-known as the venom proteome or venome-is an important step in venom research. Proteomics is used for the identification of venom toxins, determination of their primary structure including post-translational modifications, as well as investigations into the physiology underlying their production and delivery. Advances in proteomics and adjacent technologies has led to a recent upsurge in publications reporting venom proteomes. Improved mass spectrometers, better proteomic workflows, and the integration of next-generation sequencing of venom-gland transcriptomes and venomous animal genomes allow quicker and more accurate profiling of venom proteomes with greatly reduced starting material. Technologies such as imaging mass spectrometry are revealing additional insights into the mechanism, location, and kinetics of venom toxin production. However, these numerous new developments may be overwhelming for researchers designing venom proteome studies. Here, the field of venom proteomics is reviewed and some practical solutions for simplifying mass spectrometry workflows to study animal venoms are offered.

A brief review on the natural history, venomics and the medical importance of bushmaster (Lachesis) pit viper snakes

Snakes of the genus Lachesis, commonly known as bushmasters, are the largest venomous snakes in the Americas. Because these snakes have their habitats in areas of remote forests they are difficult to find, and consequently there are few studies of Lachesis taxa in their natural ecosystems. Bushmasters are distributed in tropical forest areas of South and Central America. In Brazil they can be found in the Amazon Rainforest and the Atlantic Forest. Despite the low incidence of cases, laquetic envenoming causes severe permanent sequelae due to the high amount of inoculated venom. These accidents are characterized by local pain, hemorrhage and myonecrosis that can be confused with bothropic envenomings. However, victims of Lachesis bites develop symptoms characteristic of Lachesis envenoming, known as vagal syndrome. An important message of this bibliographic synthesis exercise is that, despite having the proteomic profiles of all the taxa of the genus available, very few structure-function correlation studies have been carried out. Therefore the motivation for this review was to fill a gap in the literature on the genus Lachesis, about which there is no recent review. Here we discuss data scattered in a number of original articles published in specialized journals, spanning the evolutionary history and extant phylogeographic distribution of the bushmasters, their venom composition and diet, as well as the pathophysiology of their bites to humans and the biological activities and possible biotechnological applicability of their venom toxins.

Comparison of Methods for Measuring Protein Concentration in Venom Samples

Snake venom is an extremely interesting natural mixture of proteins and peptides, characterized by both high diversity and high pharmacological potential. Much attention has been paid to the study of venom composition of different species and also detailed analysis of the properties of individual components. Since proteins and peptides are the active ingredients in venom, rapidly developing proteomic techniques are used to analyze them. During such analyses, one of the routine operations is to measure the protein concentration in the sample. The aim of this study was to compare five methods used to measure protein content in venoms of two snake species: the Viperids representative, Agkistrodon contortrix, and the Elapids representative, Naja ashei. The study showed that for A. contortrix venom, the concentration of venom protein measured by four methods is very similar and only the NanoDrop method clearly stands out from the rest. However, in the case of N. ashei venom, each technique yields significantly different results. We hope that this report will help to draw attention to the problem of measuring protein concentration, especially in such a complex mixture as animal venoms.

Evidence for Snake Venom Plasticity in a Long-Term Study with Individual Captive Bothrops atrox

Variability in snake venom composition has been frequently reported and correlated to the adaptability of snakes to environmental conditions. Previous studies report plasticity for the venom phenotype. However, these observations are not conclusive, as the results were based on pooled venoms, which present high individual variability. Here we tested the hypothesis of plasticity by influence of confinement and single diet type in the venom composition of 13 adult specimens of Bothrops atrox snakes, maintained under captivity for more than three years. Individual variability in venom composition was observed in samples extracted just after the capture of the snakes. However, composition was conserved in venoms periodically extracted from nine specimens, which presented low variability restricted to the less abundant components. In a second group, composed of four snakes, drastic changes were observed in the venom samples extracted at different periods, mostly related to snake venom metalloproteinases (SVMPs), the core function toxins of B. atrox venom, which occurred approximately between 400 and 500 days in captivity. These data show plasticity in the venom phenotype during the lifetime of adult snakes maintained under captive conditions. Causes or functional consequences involved in the phenotype modification require further investigations.

Inflammation, angiogenesis and fibrogenesis are differentially modulated by distinct domains of the snake venom metalloproteinase jararhagin

Jararhagin, a metalloprotease from Bothrops jararaca snake venom, is a toxin containing the metalloproteinase, disintegrin-like and cysteine-rich domains; it causes acute inflammation and damage to vascular tissue. However, the actions of these domains on key components of chronic inflammation have not been determined. Our aim was to investigate the effects of jararhagin (Jar), jararhagin-C (Jar-C) and o-phenantrolin-treated jararhagin (Jar-Phe), on inflammatory response, blood vessel formation and extracellular matrix deposition in the murine sponge model. The polyether-polyurethane sponge matrix was implanted into Balb/c mice and injected daily with Jar (400 ng), Jar-Phe (400 ng), Jar-C (200 ng) or saline (control). Nine days after implantation, the sponge discs were removed and processed. In the Jar-treated implants, some of inflammatory markers (N-acetyl-β-d-glucosaminidase activity, CCL2 and TNF-α) and TGF-β1 levels were higher compared with the control group. In the Jar-C group, the inflammatory markers myeloperoxidase activity and CXCL1 were higher compared with the control. In this group, VEGF levels and collagen deposition were also higher. Jar-Phe treatment was able to inhibit the activity and/or production of MPO, CXCL1, CCL2 and TGF-β. The differential effects of these proteins in modulating the main components of fibrovascular tissue may be exploited in the management fibroproliferative diseases.

Proteopeptidomic, Functional and Immunoreactivity Characterization of Bothrops moojeni Snake Venom: Influence of Snake Gender on Venom Composition

Venom composition varies across snakes from all taxonomic levels and is influenced by the snakes’ age, habitat, diet, and sexual dimorphism. The present study reports the first in-depth investigation of venom composition in male and female Bothrops moojeni (B. moojeni) snakes (BmooM and BmooF, respectively) through three proteomics approaches associated with functional, cytotoxic, and immunoreactivity characterization. Compared with BmooM venom, BmooF venom exhibited weaker hyaluronidase, metalloproteinase, and phospholipase activity; stronger recognition by anti-bothropic serum; 1.4-fold stronger cytotoxicity; and greater number of peptides. The increased L-amino acid oxidase expression probably accounted for the stronger immunoreactivity and cytotoxicity of BmooF venom. BmooF and BmooM venom shared only 19% peptides. Some venom components were gender-specific, such as phospholipases B, phospholipase inhibitor, and hyaluronidases in BmooM, and cysteine-rich secretory proteins in BmooF. In conclusion, we describe herein the first proteomics study of B. moojeni snake venom and an in-depth characterization of gender-specific differences in venom composition. Altogether, our findings not only stress the importance of considering the snake’s gender during antivenom production, but also help to identify new potential drugs and biotechnological tools.

The Venom of the Spine-Bellied Sea Snake (Hydrophis curtus): Proteome, Toxin Diversity and Intraspecific Variation

The spine-bellied sea snake (Hydrophis curtus) is known to cause human deaths, yet its venom composition has not yet been proteomically characterised. An in-depth proteomic analysis was performed on H. curtus venom from two different seasons, January and June, corresponding to adults and subadults, respectively. Venoms from adult and subadult H. curtus individuals were compared using reversed-phase high-performance liquid chromatography (RP-HPLC), matrix-assisted laser desorption ionisation-time of flight (MALDI-TOF) mass spectrometry and liquid chromatography electrospray ionisation mass spectrometry (LC-ESI-MS) to detect intraspecific variation, and the molecular weight data obtained with ESI-MS were used to assess toxin diversity. RP-HPLC and LC-ESI-MS/MS were used to characterise the venom proteome and estimate the relative abundances of protein families present. The most abundant protein family in January and June venoms is phospholipase A₂ (PLA₂: January 66.7%; June 54.5%), followed by three-finger toxins (3FTx: January 30.4%; June 40.4%) and a minor component of cysteine-rich secretory proteins (CRISP: January 2.5%; June 5%). Trace amounts of snake venom metalloproteinases (SVMP), C-type lectins and housekeeping and regulatory proteins were also found. Although the complexity of the venom is low by number of families present, each family contained a more diverse set of isoforms than previously reported, a finding that may have implications for the development of next-generation sea snake antivenoms. Intraspecific variability was shown to be minor with one obvious exception of a 14,157-Da protein that was present in some January (adult) venoms, but not at all in June (subadult) venoms. There is also a greater abundance of short-chain neurotoxins in June (subadult) venom compared with January (adult) venom. These differences potentially indicate the presence of seasonal, ontogenetic or sexual variation in H. curtus venom.

Protein-species quantitative venomics: looking through a crystal ball

In this paper we discuss recent significant developments in the field of venom research, specifically the emergence of top-down proteomic applications that allow achieving compositional resolution at the level of the protein species present in the venom, and the absolute quantification of the venom proteins (the term “protein species” is used here to refer to all the different molecular forms in which a protein can be found. Please consult the special issue of Jornal of Proteomics “Towards deciphering proteomes via the proteoform, protein speciation, moonlighting and protein code concepts” published in 2016, vol. 134, pages 1-202). Challenges remain to be solved in order to achieve a compact and automated platform with which to routinely carry out comprehensive quantitative analysis of all toxins present in a venom. This short essay reflects the authors’ view of the immediate future in this direction for the proteomic analysis of venoms, particularly of snakes.

Strategies in 'snake venomics' aiming at an integrative view of compositional, functional, and immunological characteristics of venoms

This work offers a general overview on the evolving strategies for the proteomic analysis of snake venoms, and discusses how these may be combined through diverse experimental approaches with the goal of achieving a more comprehensive knowledge on the compositional, toxic, and immunological characteristics of venoms. Some recent developments in this field are summarized, highlighting how strategies have evolved from the mere cataloguing of venom components (proteomics/venomics), to a broader exploration of their immunological (antivenomics) and functional (toxicovenomics) characteristics. Altogether, the combination of these complementary strategies is helping to build a wider, more integrative view of the life-threatening protein cocktails produced by venomous snakes, responsible for thousands of deaths every year.

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

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

Prediction of Toxin Genes from Chinese Yellow Catfish Based on Transcriptomic and Proteomic Sequencing

Fish venom remains a virtually untapped resource. There are so few fish toxin sequences for reference, which increases the difficulty to study toxins from venomous fish and to develop efficient and fast methods to dig out toxin genes or proteins. Here, we utilized Chinese yellow catfish (Pelteobagrus fulvidraco) as our research object, since it is a representative species in Siluriformes with its venom glands embedded in the pectoral and dorsal fins. In this study, we set up an in-house toxin database and a novel toxin-discovering protocol to dig out precise toxin genes by combination of transcriptomic and proteomic sequencing. Finally, we obtained 15 putative toxin proteins distributed in five groups, namely Veficolin, Ink toxin, Adamalysin, Za2G and CRISP toxin. It seems that we have developed a novel bioinformatics method, through which we could identify toxin proteins with high confidence. Meanwhile, these toxins can also be useful for comparative studies in other fish and development of potential drugs.

Age-related association of venom gene expression and diet of predatory gastropods

Background

Venomous organisms serve as wonderful systems to study the evolution and expression of genes that are directly associated with prey capture. To evaluate the relationship between venom gene expression and prey utilization, we examined these features among individuals of different ages of the venomous, worm-eating marine snail Conus ebraeus. We determined expression levels of six genes that encode venom components, used a DNA-based approach to evaluate the identity of prey items, and compared patterns of venom gene expression and dietary specialization.

Results

C. ebraeus exhibits two major shifts in diet with age—an initial transition from a relatively broad dietary breadth to a narrower one and then a return to a broader diet. Venom gene expression patterns also change with growth. All six venom genes are up-regulated in small individuals, down-regulated in medium-sized individuals, and then either up-regulated or continued to be down-regulated in members of the largest size class. Venom gene expression is not significantly different among individuals consuming different types of prey, but instead is coupled and slightly delayed with shifts in prey diversity.

Conclusion

These results imply that changes in gene expression contribute to intraspecific variation of venom composition and that gene expression patterns respond to changes in the diversity of food resources during different growth stages.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0592-5) contains supplementary material, which is available to authorized users.

Ontogenetic Variation in Biological Activities of Venoms from Hybrids between Bothrops erythromelas and Bothrops neuwiedi Snakes

Lance-headed snakes are found in Central and South America, and they account for most snakebites in Brazil. The phylogeny of South American pitvipers has been reviewed, and the presence of natural and non-natural hybrids between different species of Bothrops snakes demonstrates that reproductive isolation of several species is still incomplete. The present study aimed to analyze the biological features, particularly the thrombin-like activity, of venoms from hybrids born in captivity, from the mating of a female Bothrops erythromelas and a male Bothrops neuwiedi, two species whose venoms are known to display ontogenetic variation. Proteolytic activity on azocoll and amidolytic activity on N-benzoyl-DL-arginine-p-nitroanilide hydrochloride (BAPNA) were lowest when hybrids were 3 months old, and increased over body growth, reaching values similar to those of the father when hybrids were 12 months old. The clotting activity on plasma diminished as hybrids grew; venoms from 3- and 6-months old hybrids showed low clotting activity on fibrinogen (i.e., thrombin-like activity), like the mother venom, and such activity was detected only when hybrids were older than 1 year of age. Altogether, these results point out that venom features in hybrid snakes are genetically controlled during the ontogenetic development. Despite the presence of the thrombin-like enzyme gene(s) in hybrid snakes, they are silenced during the first six months of life.

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.

Jararhagin-induced mechanical hyperalgesia depends on TNF-α, IL-1β and NFκB in mice

Jararhagin is a hemorrhagic metalloprotease from Bothrops jararaca snake venom. The hyperalgesic mechanisms of jararhagin were investigated focusing on the role of proinflammatory cytokines (TNF-α and IL-1β) and the transcription factor NFκB. Intraplantar administration of jararhagin (1, 10, 100 and 1000 ng/paw) induced mechanical hyperalgesia, and increased TNF-α levels at 1, 3 and 5 h, and IL-1β levels at 0.5, 1 and 3 h after its injection in the paw tissue. Pre-treatment with morphine (2, 6, 12 μg/paw) inhibited jararhagin-induced mechanical hyperagesia. The systemic or local pre-treatment with etanercept (10 mg/kg and 100 μg/paw) and IL-1ra (30 mg/kg and 100 pg/paw) inhibited jararhagin-induced mechanical hyperalgesia. Co-administration of jararhagin (0.1 ng/paw) and TNF-α (0.1 pg/paw) or jararhagin (0.1 ng/paw) and IL-1β (1 pg/paw) enhanced the mechanical hyperalgesia. The systemic or local pre-treatment with PDTC (NFκB inhibitor; 100 mg/kg and 100 μg/paw) inhibited jararhagin-induced mechanical hyperalgesia as well as PDTC decreased the jararhagin-induced production of TNF-α and IL-1β. Thus, these data demonstrate the involvement of pro-inflammatory cytokines TNF-α and IL-1β and nuclear transcription factor NFκB in jararhagin-induced mechanical hyperalgesia indicating that targeting these mechanisms might contribute to reduce the pain induced by B. jararaca snake venom.

Combined venomics, venom gland transcriptomics, bioactivities, and antivenomics of two Bothrops jararaca populations from geographic isolated regions within the Brazilian Atlantic rainforest

Bothrops jararaca is a slender and semi-arboreal medically relevant pit viper species endemic to tropical and subtropical forests in southern Brazil, Paraguay, and northern Argentina (Misiones). Within its geographic range, it is often abundant and is an important cause of snakebite. Although no subspecies are currently recognized, geographic analyses have revealed the existence of two well-supported B. jararaca clades that diverged during the Pliocene ~3.8Mya and currently display a southeastern (SE) and a southern (S) Atlantic rainforest (Mata Atlântica) distribution. The spectrum, geographic variability, and ontogenetic changes of the venom proteomes of snakes from these two B. jararaca phylogroups were investigated applying a combined venom gland transcriptomic and venomic analysis. Comparisons of the venom proteomes and transcriptomes of B. jararaca from the SE and S geographic regions revealed notable interpopulational variability that may be due to the different levels of population-specific transcriptional regulation, including, in the case of the southern population, a marked ontogenetic venom compositional change involving the upregulation of the myotoxic PLA2 homolog, bothropstoxin-I. This population-specific marker can be used to estimate the proportion of venom from the southern population present in the B. jararaca venom pool used for the Brazilian soro antibotrópico (SAB) antivenom production. On the other hand, the southeastern population-specific D49-PLA2 molecules, BinTX-I and BinTX-II, lend support to the notion that the mainland ancestor of Bothrops insularis was originated within the same population that gave rise to the current SE B. jararaca phylogroup, and that this insular species endemic to Queimada Grande Island (Brazil) expresses a pedomorphic venom phenotype. Mirroring their compositional divergence, the two geographic B. jararaca venom pools showed distinct bioactivity profiles. However, the SAB antivenom manufactured in Vital Brazil Institute neutralized the lethal effect of both venoms to a similar extent. In addition, immobilized SAB antivenom immunocaptured most of the venom components of the venoms of both B. jararaca populations, but did not show immunoreactivity against vasoactive peptides. The Costa Rican bothropic-crotalic-lachesic (BCL) antivenom showed the same lack of reactivity against vasoactive peptides but, in addition, was less efficient immunocapturing PI- and PIII-SVMPs from the SE venom, and bothropstoxin-I, a CRISP molecule, and a D49-PLA2 from the venom of the southern B. jararaca phylogroup. The remarkable paraspecificity exhibited by the Brazilian and the Costa Rican antivenoms indicates large immunoreactive epitope conservation across the natural history of Bothrops, a genus that has its roots in the middle Miocene. This article is part of a Special Issue entitled: Omics Evolutionary Ecolog.

Early significant ontogenetic changes in snake venoms

Snake venom plays a critical role in food acquisition, digestion, and defense. Venoms are known to change throughout the life of some snake species, but nothing is known about the venom composition of hatchling/neonate snakes prior to and just after their first shedding cycle, despite this being a critical time in the life of the snake. Using a cohort of Crotalus horridus and two cohorts of Crotalus adamanteus, we showed for the first time that snakes undergo significant changes in venom composition after the postnatal shedding event. The number of changes among cohorts ranged widely and there was wide variation in the direction of protein regulation, which appeared to be on a locus-specific level rather than protein-family level. These significant venom composition changes that take place in the first few weeks of life most likely play critical roles in venom economy and resource conservation and may partially explain the rare, post-birth maternal care found in some venomous species.

Omics meets biology: application to the design and preclinical assessment of antivenoms

Snakebite envenoming represents a neglected tropical disease that has a heavy public health impact worldwide, mostly affecting poor people involved in agricultural activities in Africa, Asia, Latin America and Oceania. A key issue that complicates the treatment of snakebite envenomings is the poor availability of the only validated treatment for this disease, antivenoms. Antivenoms can be an efficacious treatment for snakebite envenoming, provided they are safe, effective, affordable, accessible and administered appropriately. The shortage of antivenoms in various regions, particularly in Sub-Saharan Africa and some parts of Asia, can be significantly alleviated by optimizing the use of current antivenoms and by the generation of novel polyspecific antivenoms having a wide spectrum of efficacy. Complementing preclinical testing of antivenom efficacy using in vivo and in vitro functional neutralization assays, developments in venomics and antivenomics are likely to revolutionize the design and preclinical assessment of antivenoms by being able to test new antivenom preparations and to predict their paraspecific neutralization to the level of species-specific toxins.

Peptidomics combined with cDNA library unravel the diversity of centipede venom

Centipedes are one of the oldest venomous arthropods using toxin as their weapon to capture prey. But little attention was focused on them and only few centipede toxins were demonstrated with activity on ion channels. Therefore, more deep works are needed to understand the diversity of centipede venom. In the present study, we use peptidomics combined with cDNA library to uncover the diversity of centipede Scolopendra subspinipes mutilans L. Koch. 192 peptides were identified by LC-MS/MS and 79 precursors were deduced by cDNA library. Surprisingly, the signal peptides of centipede toxins were more complicated than any other animal toxins and even exhibited large differences in homologues. Meanwhile, a large number of variants generated by alternative cleavage sites were detected by mass spectra. Odd number of cystein (3, 5, 7) found in the mature peptides were seldom seen in peptide toxins. In additional, two novel cysteine frameworks (C-C-C-CCC, C-C-C-C-CC-CC) were identified from 16 different cysteine frameworks from centipede peptides. Only 29 precursors have clear targets, while others may provide a potential diversity function for centipede. These findings highlight the extensive diversity of centipede toxins and provide powerful tools to understand the capture and defense weapon of centipede.

BIOLOGICAL SIGNIFICANCE

Peptide toxins from venomous animal have attracted increasing attentions due to their extraordinary chemical and pharmacological diversity. Centipedes are one of the most used Chinese traditional medicines, but little was known about the active components. The venom of Scolopendra subspinipes mutilans L. Koch is first deeply analyzed in this work and most of peptides were never discovered before. Interestingly, the number and arrangement of cysteine showed a larger different to known peptide toxins such spider or scorpion toxins. Moreover, only 29 peptides from this centipede venom were identified with known function. It suggested that our work not only important to understand the composition of centipede venom, but also provide many valuable peptides for potential biological functions.

Proteomic and biochemical analyses of short-tailed pit viper (Gloydius brevicaudus) venom: age-related variation and composition-activity correlation

We conducted an in-depth analysis of the proteomic and biochemical profiles of the venom of neonate and adult short-tailed pit vipers (Gloydius brevicaudus). Identified proteins were assigned to a few main toxin families. Disintegrin, phospholipase A2 (PLA2), serine proteinase, cysteine-rich secretory protein, C-type lectin-like protein, l-amino acid oxidase and snake venom metalloproteinase (SVMP) were detected in both venoms, while 5'-nucleotidase was detected only in the adult venom. SVMP was the predominant protein family in both venoms (neonate: 65.7%; adult: 64.4%), followed by PLA2 (neonate: 13.4%; adult: 25.0%). Antivenomic analysis revealed that commercial G. brevicaudus antivenom almost neutralized the chromatographic peaks with medium and high molecular masses in both venoms, but did not completely recognize peaks with low molecular mass. Toxicological and enzymatic activities show remarkable age-related variation in G. brevicaudus venom, probably resulting from variation in venom composition. Our data demonstrate age-related variation across venomics, antivenomics and biochemical profiles of G. brevicaudus venom, and have implications for the management of G. brevicaudus bites, including improving antivenom preparation by combining both venoms.

BIOLOGICAL SIGNIFICANCE

This study investigates the composition and biochemical activity of neonate and adult Gloydius brevicaudus venoms. We found remarkable age-related variation in venom biological activity, likely the result of variation in venom composition. Antivenomics analysis was used to explore difference in neonate and adult G. brevicaudus venoms. Our findings have implications for the diagnosis and clinical management of G. brevicaudus bites, and the design of venom mixtures that will increase the efficacy of commercial antivenom. This article is part of a Special Issue entitled: Proteomics of non-model organisms.