Snake Venomics of Bothriechis nigroviridis Reveals Extreme Variability among Palm Pitviper Venoms: Different Evolutionary Solutions for the Same Trophic Purpose

Exploring venom diversity in Mixcoatlus browni and Mixcoatlus barbouri: A comparative analysis of two rare Mexican snake species with crotoxin-like presence

The genus Mixcoatlus is composed of three species: Mixcoatlus barbouri, M. browni, and M. melanurus, of which the venom composition of M. melanurus, the most common species of the three, has only recently been described. However, very little is known about the natural history of M. barbouri and M. browni, and the venom composition of these two species has remained thus far unexplored. In this study we characterize the proteomic profiles and the main biochemical and toxic activities of these two venoms. Proteomic data obtained by shotgun analysis of whole venom identified 12 protein families for M. barbouri, and 13 for M. browni. The latter venom was further characterized by using a quantitative 'venomics' protocol, which revealed that it is mainly composed of 51.1 % phospholipases A2 (PLA2), 25.5 % snake venom serine proteases (SVSP), 4.6 % l-amino oxidases (LAO), and 3.6 % snake venom metalloproteases (SVMP), with lower percentages other six protein families. Both venoms contained homologs of the basic and acidic subunits of crotoxin. However, due to limitations in M. barbouri venom availability, we could only characterize the crotoxin-like protein of M. browni venom, which we have named Mixcoatlutoxin. It exhibited a lethal potency in mice like that described for classical rattlesnake crotoxins. These findings expand knowledge on the distribution of crotoxin-like heterodimeric proteins in viper snake species. Further investigation of the bioactivities of the venom of M. barbouri, on the other hand, remains necessary.

Envenomation by Eyelash Viper Bothriechis schlegelii (Berthold, 1846) in Southwestern Colombia

INTRODUCTION

Bothriechis schlegelii is a Crotaline viperid species of Central America and Northern South America. The characteristics of its envenomation have not been well established. We present clinical characteristics of human cases evaluated and treated in a hospital in southwestern Colombia.

METHODS

We evaluated data from patients who suffered Bothriechis schlegelii envenomation and were seen at Fundación Valle del Lili Hospital, Cali, Colombia between 2011 and 2022.

RESULTS

Eight patients were included, with a median age of 24 years. Snakebites occurred in rural areas. Six (75%) patients were bitten on the upper extremities in relation to the arboreal habits of this animal. The most common symptoms were pain and edema (N = 8, 100%), ecchymoses (N = 2, 25%), and paresthesia (N = 2, 25%). The most common systemic findings were hypofibrinogenemia (N = 8, 100%) and prolonged prothrombin time in five patients (N = 5, 62.5%). All were treated with polyvalent antivenom for Colombian snakes, with a good response and outcome.

CONCLUSIONS

Most bite sites from B. schlegelii were on the upper limbs. All patients had both local manifestations, including edema, pain, and systemic effects with hypofibrinogenemia, but none had systemic bleeding. Every patient received antivenom and had favorable outcomes.

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

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.

Differential coagulotoxic and neurotoxic venom activity from species of the arboreal viperid snake genus Bothriechis (palm-pitvipers)

The viperid snake genus Bothriechis consists of eleven species distributed among Central and South America, living across low and high-altitude habitats. Despite Bothriechis envenomations being prominent across the Central and South American region, the functional effects of Bothriechis venoms are poorly understood. Thus, the aim of this study was to investigate the coagulotoxic and neurotoxic activities of Bothriechis venoms to fill this knowledge gap. Coagulotoxic investigations revealed Bothriechis nigroviridis and B. schlegelii to have pseudo-procoagulant venom activity, forming weak clots that rapidly break down, thereby depleting fibrinogen levels and thus contributing to a net anticoagulant state. While one sample of B. lateralis also showed weaker pseudo-procoagulant activity, directly clotting fibrinogen, two samples of B. lateralis venom were anticoagulant through the inhibition of thrombin and factor Xa activity. Differential efficacy of PoliVal-ICP antivenom was also observed, with the pseudo-procoagulant effect of B. nigroviridis venom poorly neutralised, despite this same activity in the venom of B. schlegelii being effectively neutralised. Significant specificity of these fibrinogen cleaving toxins was also observed, with no activity upon model amphibian, avian, lizard or rodent plasma observed. However, upon avian plasma the venom of B. nigroviridis exerted a complete anticoagulant effect, in contrast to the pseudo-procoagulant effect seen on human plasma. Neurotoxic investigations revealed B. schlegelii to be unique among the genus in having potent binding to the orthosteric site of the alpha-1 postsynaptic nicotinic acetylcholine receptor (with B. lateralis having a weaker but still discernible effect). This represents the first identification of postsynaptic nAChR neurotoxic activity for Bothriechis. In conclusion this study identifies notable differential activity within the coagulotoxic and postsynaptic neurotoxic activity of Bothriechis venoms, supporting previous research, and highlights the need for further studies with respect to antivenom efficacy as well as coagulotoxin specificity for Bothriechis venoms.

Advances in venomics: Modern separation techniques and mass spectrometry

Snake venoms are complex chemical mixtures of biologically active proteins and non-protein components. Toxins have a wide range of targets and effects to include ion channels and membrane receptors, and platelet aggregation and platelet plug formation. Toxins target these effectors and effects at high affinity and selectivity. From a pharmacological perspective, snake venom compounds are a valuable resource for drug discovery and development. However, a major challenge to drug discovery using snake venoms is isolating and analyzing the bioactive proteins and peptides in these complex mixtures. Getting molecular information from complex mixtures such as snake venoms requires proteomic analyses, generally combined with transcriptomic analyses of venom glands. The present review summarizes current knowledge and highlights important recent advances in venomics with special emphasis on contemporary separation techniques and bioinformatics that have begun to elaborate the complexity of snake venoms. Several analytical techniques such as two-dimensional gel electrophoresis, RP-HPLC, size exclusion chromatography, ion exchange chromatography, MALDI-TOF-MS, and LC-ESI-QTOF-MS have been employed in this regard. The improvement of separation approaches such as multidimensional-HPLC, 2D-electrophoresis coupled to soft-ionization (MALDI and ESI) mass spectrometry has been critical to obtain an accurate picture of the startling complexity of venoms. In the case of bioinformatics, a variety of software tools such as PEAKS also has been used successfully. Such information gleaned from venomics is important to both predicting and resolving the biological activity of the active components of venoms, which in turn is key for the development of new drugs based on these venom components.

Comparative characterization of Viperidae snake venoms from Perú reveals two compositional patterns of phospholipase A2 expression

Snake species within the Bothrops complex (sensu lato) are of medical relevance in Latin America, but knowledge on their venom characteristics is limited, or even unavailable, for some taxa. Perú harbors 17 species of pit vipers, within the genera Bothrops, Bothriechis, Bothrocophias, Porthidium, Crotalus, and Lachesis. This study compared the venoms of twelve species, through chromatographic and electrophoretic profiles, as well as proteolytic and phospholipase A₂ (PLA₂) activities. Also, proteomic profiles were analyzed for nine of the venoms using a shotgun approach. Results unveiled conspicuous differences in the expression of venom PLA₂s among species, six of them presenting scarce levels as judged by RP-HPLC profiles. Since most species within the bothropoid lineage possess venoms with high to intermediate abundances of this protein family, our findings suggest the existence of a phenotypic duality in the expression of venom PLA₂s within the Bothrops (sensu lato) complex. Bothrops barnetti and Bothrocophias andianus venoms, very scarce in PLA₂s, were shown to lack significant myotoxic activity, highlighting that the observed variability in PLA₂ expression bears toxicological correlations with effects attributed to these proteins. Finally, an attempt to identify phylogenetic relationships of bothropoid species from Perú presenting low- or high-PLA₂ venom phenotypes showed an interspersed pattern, thus precluding a simple phylogenetic interpretation of this venom compositional dichotomy.

•

Venoms from 12 viperids of Perú were compared.

•

Conspicuous differences in the expression of PLA₂ were found.

•

Venoms presenting scarce levels of PLA₂ lack myotoxicity.

•

A new phenotypic dichotomy in venom PLA₂ expression is described within Bothrops (sensu lato).

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.

Proteomic profiling, functional characterization, and immunoneutralization of the venom of Porthidium porrasi, a pitviper endemic to Costa Rica

The genus Porthidium includes nine pitviper species inhabiting Mexico, Central America, and northern South America. Porthidium porrasi is a species endemic to the Southwest of Costa Rica, for which no information on its venom was available. In this study, the proteomic composition and functional activities of P. porrasi venom are described. The most abundant venom proteins were identified as metalloproteinases (36.5%). In descending order of abundance, proteins belonging to the disintegrin, phospholipase A₂, serine proteinase, C-type lectin/lectin-like, vascular endothelial growth factor, Cysteine-rich secretory protein, L-amino acid oxidase, phospholipase B, and phosphodiesterase families were also identified. P. porrasi venom showed a weak lethal potency in mice (10 μg/g body weight by intraperitoneal route), induced marked hemorrhage and edema, and weak myotoxic effect. These in vivo activities, as well as those assayed in vitro (proteolytic and phospholipase A₂ activities) correlated with compositional data. A comparison of P. porrasi venom with those of three other Porthidium species studied to date reveals a generally conserved compositional and functional pattern in this pitviper genus. Importantly, the lethal effect of P. porrasi venom in mice was adequately cross-neutralized by a heterospecific polyvalent antivenom, supporting its use in the treatment of eventual envenomings by this species.

Venom of prey-specialized spiders is more toxic to their preferred prey: A result of prey-specific toxins

In specialized predators, a variety of adaptations have evolved to such a level of specificity that they allow very effective exploitation of focal prey. Venom is an essential adaptive trait of predatory venomous species, such as spiders, yet our knowledge of spider venom is incomplete. In agreement with the prey preference hypothesis, we expected that the venom of spider specialists should be more toxic to focal than to alternative prey, because it is composed of prey-specific toxins. Here we used spiders with three types of trophic specializations: specialists that were ant-eating, termite-eating and spider-eating. We compared the efficacy of prey capture of preferred and alternative prey (measured as paralysis latency) with that of related generalists and profiled the venom of the studied species using proteomic methods. We used 22 spider species: six myrmecophagous, two termitophagous, three araneophagous and 11 euryphagous generalist species belonging to different families. We found that ten of the eleven specialist species induced significantly shorter paralysis latency in preferred prey than in alternative prey. Generalists exhibited either similar efficiency on both prey types or slightly higher efficiency on preferred prey. Multivariate analysis of proteomic profiles (peptides and proteins) revealed significant differences between trophic specializations, particularly in peptides. Specialists appear to have venom composed of unique specific compounds as revealed by the multivariate ordination and indicator analysis. These components are likely prey-specific toxins.

Venom Ontogeny in the Mexican Lance-Headed Rattlesnake (Crotalus polystictus)

As trophic adaptations, rattlesnake venoms can vary in composition depending on several intrinsic and extrinsic factors. Ontogenetic changes in venom composition have been documented for numerous species, but little is known of the potential age-related changes in many rattlesnake species found in México. In the current study, venom samples collected from adult and neonate Crotalus polystictus from Estado de México were subjected to enzymatic and electrophoretic analyses, toxicity assays (LD₅₀), and MALDI-TOF mass spectrometry, and a pooled sample of adult venom was analyzed by shotgun proteomics. Electrophoretic profiles of adult males and females were quite similar, and only minor sex-based variation was noted. However, distinct differences were observed between venoms from adult females and their neonate offspring. Several prominent bands, including P-I and P-III snake venom metalloproteinases (SVMPs) and disintegrins (confirmed by MS/MS) were present in adult venoms and absent/greatly reduced in neonate venoms. Age-dependent differences in SVMP, kallikrein-like, phospholipase A₂ (PLA₂), and L-amino acid oxidase (LAAO) activity levels were confirmed by enzymatic activity assays, and like many other rattlesnake species, venoms from adult snakes have higher SVMP activity than neonate venoms. Conversely, PLA₂ activity was approximately 2.5 × greater in venoms from neonates, likely contributing to the increased toxicity (neonate venom LD₅₀ = 4.5 μg/g) towards non-Swiss albino mice when compared to adult venoms (LD₅₀ = 5.5 μg/g). Thrombin-like (TLE) and phosphodiesterase activities did not vary significantly with age. A significant effect of sex (between adult male and adult female venoms) was also observed for SVMP, TLE, and LAAO activities. Analysis of pooled adult venom by LC-MS/MS identified 14 toxin protein families, dominated by bradykinin-inhibitory peptides, SVMPs (P-I, P-II and P-III), disintegrins, PLA₂s, C-type-lectins, CRiSPs, serine proteinases, and LAAOs (96% of total venom proteins). Neonate and adult C. polystictus in this population consume almost exclusively mammals, suggesting that age-based differences in composition are related to physical differences in prey (e.g., surface-to-volume ratio differences) rather than taxonomic differences between prey. Venoms from adult C. polystictus fit a Type I pattern (high SVMP activity, lower toxicity), which is characteristic of many larger-bodied rattlesnakes of North America.

Comparative Venomics of the Vipera ammodytes transcaucasiana and Vipera ammodytes montandoni from Turkey Provides Insights into Kinship

The Nose-horned Viper (Vipera ammodytes) is one of the most widespread and venomous snakes in Europe, which causes high frequent snakebite accidents. The first comprehensive venom characterization of the regional endemic Transcaucasian Nose-horned Viper (Vipera ammodytes transcaucasiana) and the Transdanubian Sand Viper (Vipera ammodytes montandoni) is reported employing a combination of intact mass profiling and bottom-up proteomics. The bottom-up analysis of both subspecies identified the major snake protein families of viper venoms. Furthermore, intact mass profiling revealed the presence of two tripeptidic metalloprotease inhibitors and their precursors. While previous reports applied multivariate analysis techniques to clarify the taxonomic status of the subspecies, an accurate classification of Vipera ammodytestranscaucasiana is still part of the ongoing research. The comparative analysis of the viper venoms on the proteome level reveals a close relationship between the Vipera ammodytes subspecies, which could be considered to clarify the classification of the Transcaucasian Nose-horned Viper. However, the slightly different ratio of some venom components could be indicating interspecific variations of the two studied subspecies or intraspecies alternations based on small sample size. Additionally, we performed a bioactivity screening with the crude venoms against several human cancerous and non-cancerous cell lines, which showed interesting results against a human breast adenocarcinoma epithelial cell line. Several fractions of Vipera a. transcaucasiana demonstrated a strong cytotoxic effect on triple negative MDA MB 231 breast cancer cells.

A Review and Database of Snake Venom Proteomes

Advances in the last decade combining transcriptomics with established proteomics methods have made possible rapid identification and quantification of protein families in snake venoms. Although over 100 studies have been published, the value of this information is increased when it is collated, allowing rapid assimilation and evaluation of evolutionary trends, geographical variation, and possible medical implications. This review brings together all compositional studies of snake venom proteomes published in the last decade. Compositional studies were identified for 132 snake species: 42 from 360 (12%) Elapidae (elapids), 20 from 101 (20%) Viperinae (true vipers), 65 from 239 (27%) Crotalinae (pit vipers), and five species of non-front-fanged snakes. Approximately 90% of their total venom composition consisted of eight protein families for elapids, 11 protein families for viperines and ten protein families for crotalines. There were four dominant protein families: phospholipase A₂s (the most common across all front-fanged snakes), metalloproteases, serine proteases and three-finger toxins. There were six secondary protein families: cysteine-rich secretory proteins, l-amino acid oxidases, kunitz peptides, C-type lectins/snaclecs, disintegrins and natriuretic peptides. Elapid venoms contained mostly three-finger toxins and phospholipase A₂s and viper venoms metalloproteases, phospholipase A₂s and serine proteases. Although 63 protein families were identified, more than half were present in <5% of snake species studied and always in low abundance. The importance of these minor component proteins remains unknown.

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.

Biochemical and biological characterization of Bothriechis schlegelii snake venoms from Colombia and Costa Rica

Snakebites inflicted by the arboreal viperid snake Bothriechis schlegelii in humans are characterized by pain, edema, and ecchymosis at the site of the bite, rarely with blisters, local necrosis, or defibrination. Herein, a comparative study of Bothriechis schlegelii snake venoms from Colombia (BsCo) and Costa Rica (BsCR) was carried out in order to compare their main activities and to verify the efficacy of Bothrops antivenom produced in Brazil to neutralize them. Biochemical (SDS-PAGE and zymography) and biological parameters (edematogenic, lethal, hemorrhagic, nociceptive, and phospholipase A₂ activities) induced by BsCo and BsCR snake venoms were evaluated. The presence of antibodies in Bothrops antivenom that recognize BsCo and BsCR snake venoms by enzyme-linked immunosorbent assay and Western blotting, as well as the ability of this antivenom to neutralize the toxic activities were also verified. SDS-PAGE showed differences between venoms. Distinctive caseinolytic and hyaluronidase patterns were detected by zymography. BsCo and BsCR showed similar phospholipase A₂ activity. Strong cross-reactivity between BsCo and BsCR was detected using Bothrops antivenom with many components located between 150 and 35 kDa. BsCR was more edematogenic and almost fourfold more hemorrhagic than BsCo, and both venoms induced nociception. BsCR (LD₅₀ 5.60 mg/kg) was more lethal to mice than BsCo (LD₅₀ 9.24 mg/kg). Bothrops antivenom was effective in the neutralization of lethal and hemorrhagic activities of BsCo and BsCR and was partially effective in the neutralization of edematogenic and nociceptive activities. In conclusion, geographic distribution influences the composition and activities of Bothriechis schlegelii venoms. Bothrops antivenom cross-reacted with these venoms and was partially effective in neutralizing some toxic activities of BsCo and BsCR.

Venomics of Calloselasma rhodostoma, the Malayan pit viper: A complex toxin arsenal unraveled

The venom of Malayan pit viper (Calloselasma rhodostoma) is highly toxic but also valuable in drug discovery. However, a comprehensive proteome of the venom that details its toxin composition and abundance is lacking. This study aimed to unravel the venom complexity through a multi-step venomic approach. At least 96 distinct proteins (29 basic, 67 acidic) in 11 families were identified from the venom. The venom consists of mainly snake venom metalloproteinases (SVMP, 41.17% of total venom proteins), within which the P-I (kistomin, 20.4%) and P-II (rhodostoxin, 19.8%) classes predominate. This is followed by C-type lectins (snaclec, 26.3%), snake venom serine protease (SVSP, 14.9%), L-amino acid oxidase (7.0%), phospholipase A2 (4.4%), cysteine-rich secretory protein (2.5%), and five minor toxins (nerve growth factor, neurotrophin, phospholipase B, 5' nucleotidase and phosphodiesterase, totaling 2.6%) not reported in the proteome hitherto. Importantly, all principal hemotoxins unveiled correlate with the syndrome: SVSP ancrod causes venom-induced consumptive coagulopathy, aggravated by thrombocytopenia caused by snaclec rhodocytin, a platelet aggregation inducer, while P-II rhodostoxin mediates hemorrhage, exacerbated by P-I kistomin and snaclec rhodocetin that inhibit platelet plug formation. These toxins exist in multiple isoforms and/or complex subunits, deserving further characterization for the development of an effective, polyspecific regional antivenom.

BIOLOGICAL SIGNIFICANCE

Advents in proteomics and bioinformatics have vigorously propelled the scientific discoveries of toxins from various lineages of venomous snakes. The Malayan pit viper, Calloselasma rhodostoma, is a medically important species in Southeast Asia as its bite can cause envenomation, while the venom is also a source of bioactive compounds for drug discovery. Detailed profiling of the venom, however, is inadequate possibly due to the complex nature of the venom and technical limitation in separating the constituents into details. Integrating a multi-step fractionation method, this study successfully revealed a comprehensive and quantitative profile of the composition of the venom of this medically important venomous snake. The relative abundance of the various venom proteins is determined in a global profile, providing useful information for understanding the pathogenic roles of the different toxins in C. rhodostoma envenomation. Notably, the principal hemotoxins were identified in great details, including the variety of toxin subunits and isoforms. The findings indicate that these toxins are the principal targets for effective antivenom neutralization, and should be addressed in the production of a pan-regional polyspecific antivenom. In addition, minor toxin components not reported previously in the venom were also detected in this study, enriching the current toxin database for the venomous snakes.

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.

Snake Venomics and Antivenomics of Bothrops diporus, a Medically Important Pitviper in Northeastern Argentina

Snake species within genus Bothrops are responsible for more than 80% of the snakebites occurring in South America. The species that cause most envenomings in Argentina, B. diporus, is widely distributed throughout the country, but principally found in the Northeast, the region with the highest rates of snakebites. The venom proteome of this medically relevant snake was unveiled using a venomic approach. It comprises toxins belonging to fourteen protein families, being dominated by PI- and PIII-SVMPs, PLA₂ molecules, BPP-like peptides, L-amino acid oxidase and serine proteinases. This toxin profile largely explains the characteristic pathophysiological effects of bothropic snakebites observed in patients envenomed by B. diporus. Antivenomic analysis of the SAB antivenom (Instituto Vital Brazil) against the venom of B. diporus showed that this pentabothropic antivenom efficiently recognized all the venom proteins and exhibited poor affinity towards the small peptide (BPPs and tripeptide inhibitors of PIII-SVMPs) components of the venom.

Heterodimeric l-amino acid oxidase enzymes from Egyptian Cerastes cerastes venom: Purification, biochemical characterization and partial amino acid sequencing

Two l-amino acid oxidase enzyme isoforms, Cc-LAAOI and Cc-LAAOII were purified to apparent homogeneity from Cerastes cerastes venom in a sequential two-step chromatographic protocol including; gel filtration and anion exchange chromatography. The native molecular weights of the isoforms were 115 kDa as determined by gel filtration on calibrated Sephacryl S-200 column, while the monomeric molecular weights of the enzymes were, 60, 56 kDa and 60, 53 kDa for LAAOI and LAAOII, respectively. The tryptic peptides of the two isoforms share high sequence homology with other snake venom l-amino acid oxidases. The optimal pH and temperature values of Cc-LAAOI and Cc-LAAOII were 7.8, 50 °C and 7, 60 °C, respectively. The two isoenzymes were thermally stable up to 70 °C. The K m and V max values were 0.67 mM, 0.135 μmol/min for LAAOI and 0.82 mM, 0.087 μmol/min for LAAOII. Both isoenzymes displayed high catalytic preference to long-chain, hydrophobic and aromatic amino acids. The Mn2 + ion markedly increased the LAAO activity for both purified isoforms, while Na+, K+, Ca2 +, Mg2 + and Ba2 + ions showed a non-significant increase in the enzymatic activity of both isoforms. Furthermore, Zn2 +, Ni2 +, Co2 +, Cu2 + and AL3 + ions markedly inhibited the LAAOI and LAAOII activities. l-Cysteine and reduced glutathione completely inhibited the LAAO activity of both isoenzymes, whereas, β-mercaptoethanol, O-phenanthroline and PMSF completely inhibited the enzymatic activity of LAAOII. Furthermore, iodoacitic acid inhibited the enzymatic activity of LAAOII by 46% and had no effect on the LAAOI activity.

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.

Proteomic investigation of Sri Lankan hump-nosed pit viper (Hypnale hypnale) venom

The hump-nosed pit viper, Hypanle hypnale, contributes to snakebite mortality and morbidity in Sri Lanka. Studies showed that the venom is hemotoxic and nephrotoxic, with some biochemical and antigenic properties similar to the venom of Calloselasma rhodostoma (Malayan pit viper). To further characterize the complexity composition of the venom, we investigated the proteome of a pooled venom sample from >10 Sri Lankan H. hypnale with reverse-phase high performance liquid chromatography (rp-HPLC), sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and peptide sequencing (tandem mass-spectrometry and/or N-terminal sequencing). The findings ascertained that two phospholipase A2 subtypes (E6-PLA2, W6-PLA2) dominate the toxin composition by 40.1%, followed by snake venom metalloproteases (36.9%), l-amino acid oxidase (11.9%), C-type lectins (5.5%), serine proteases (3.3%) and others (2.3%). The presence of the major toxins correlates with the venom's major pathogenic effects, indicating these to be the principal target toxins for antivenom neutralization. This study supports the previous finding of PLA2 dominance in the venom but diverges from the view that H. hypnale venom has low expression of large enzymatic toxins. The knowledge of the composition and abundance of toxins is essential to elucidate the pathophysiology of H. hypnale envenomation and to optimize antivenom formulation in the future.

Jellyfish venomics and venom gland transcriptomics analysis of Stomolophus meleagris to reveal the toxins associated with sting

Jellyfish Stomolophus meleagris is a very dangerous animal because of its strong toxicity. However, the composition of the venom is still unclear. Both proteomics and transcriptomics approaches were applied in present study to investigate the major components and their possible relationships to the sting. The proteomics of the venom from S. meleagris was conducted by tryptic digestion of the crude venom followed by RP-HPLC separation and MS/MS analysis of the tryptic peptides. The venom gland transcriptome was analyzed using a high-throughput Illumina sequencing platform HiSeq 2000 with de novo assembly. A total of 218 toxins were identified including C-type lectin, phospholipase A₂ (PLA₂), potassium channel inhibitor, protease inhibitor, metalloprotease, hemolysin and other toxins, most of which should be responsible for the sting. Among them, serine protease inhibitor, PLA₂, potassium channel inhibitor and metalloprotease are predominant, representing 28.44%, 21.56%, 16.06% and 15.14% of the identified venom proteins, respectively. Overall, our combined proteomics and transcriptomics approach provides a systematic overview of the toxins in the venom of jellyfish S. meleagris and it will be significant to understand the mechanism of the sting.

BIOLOGICAL SIGNIFICANCE

Jellyfish Stomolophus meleagris is a very dangerous animal because of its strong toxicity. It often bloomed in the coast of China in recent years and caused thousands of people stung and even deaths every year. However, the components which caused sting are still unknown yet. In addition, no study about the venomics of jellyfish S. meleagris has been reported. In the present study, both proteomics and transcriptomics approaches were applied to investigate the major components related to the sting. The result showed that major component included C-type lectin, phospholipase A₂, potassium channel inhibitor, protease inhibitor, metalloprotease, hemolysin and other toxins, which should be responsible for the effect of sting. This is the first research about the venomics of jellyfish S. meleagris. It will be significant to understand the mechanism of the biological effects and helpful to develop ways to deal with the sting.

Two color morphs of the pelagic yellow-bellied sea snake, Pelamis platura, from different locations of Costa Rica: snake venomics, toxicity, and neutralization by antivenom

The yellow-bellied sea snake, Pelamis platura, is the most broadly distributed snake species. Despite being endowed with a highly lethal venom, a proteomic analysis of its toxin composition was unavailable. The venoms of specimens collected in Golfo de Papagayo and Golfo Dulce (Costa Rica), where two distinctive color morphs occur, were chromatographically compared. The latter inhabits a fjord-like gulf where the transit of oceanic sea snakes into and from the basin is restricted, thus possibly affecting gene flow. RP-HPLC evidenced a conserved venom protein profile in both populations, despite their divergent color phenotypes. Following a trend observed in other sea snakes, P. platura venom is relatively simple, being composed of proteins of the three-finger toxin (3FTx), phospholipase A2 (PLA2), cysteine-rich secretory protein (CRISP), 5'-nucleotidase, and metalloproteinase families. The first three groups represent 49.9%, 32.9%, and 9.1% of total venom protein, respectively. The most abundant component (~26%) is pelamitoxin (P62388), a short-chain 3FTx, followed by a major basic PLA2 (~20%) and a group of three isoforms of CRISPs (~9%). Whereas isolated pelamitoxin was highly lethal to mice, neither the PLA2 nor the CRISP fraction caused death. However, the PLA2 rapidly increased plasma creatine kinase activity after intramuscular injection, indicating its myotoxic action. Differing from myotoxic PLA2s of viperids, this PLA2 was not cytolytic to murine myogenic cells in vitro, suggesting possible differences in its mechanism of action. The median lethal dose (LD50) estimates for P. platura crude venom in mice and in three species of fishes did not differ significantly. The sea snake antivenom manufactured by CSL Ltd. (Australia), which uses Enhydrina schistosa as immunogen, cross-recognized the three major components of P. platura venom and, accordingly, neutralized the lethal activity of crude venom and pelamitoxin, therefore being of potential usefulness in the treatment of envenomations by this species.

BIOLOGICAL SIGNIFICANCE

Integrative analyses of animal venoms that combine the power of proteomics (venomics) with the characterization of their functional and immunological properties are significantly expanding knowledge on these remarkable bioweapons, both from a basic and a medical perspective. Costa Rica harbors a unique population of the yellow-bellied sea snake, Pelamis platura, that is restricted to a fjord-like gulf (Golfo Dulce). This population differs markedly from oceanic populations found elsewhere along the Pacific coast of this country, by presenting a patternless bright yellow coloration, instead of the typical bicolored or tricolored pattern of this species. It has been suggested that the dominance of this yellow-morph in Golfo Dulce might reflect gene flow restrictions, caused by the oceanographic conditions at this location. The present study demonstrates that the remarkable phenotypic variation between the two color morphs inhabiting Golfo Dulce and Golfo de Papagayo, respectively, is not associated with differences in the expression of venom components, as shown by their conserved RP-HPLC profiles. Proteomic analysis revealed the relatively simple toxin composition of P. platura venom, which contains three predominant types of proteins: three-finger toxins (protein abundance: 49.9%), phospholipases A2 (32.9%), and cysteine-rich secretory proteins (9.1%), together with few minor components. Further, the involvement of these most abundant proteins in the toxic effects of the venom, and their cross-recognition and neutralization by a sea snake antivenom produced against the venom of Enhydrina schistosa, were analyzed.

Immunological profile of antivenoms: preclinical analysis of the efficacy of a polyspecific antivenom through antivenomics and neutralization assays

Parenteral administration of animal-derived antivenoms constitutes the mainstay in the treatment of snakebite envenomings. Despite the fact that this therapy has been available for over a century, the detailed understanding of the neutralizing and immunoreactivity profiles of the majority of antivenoms is pending. Currently, a combination of preclinical neutralization tests and 'antivenomics', i.e. a proteomic-based assessment of antivenom immunoreactivity, provides a powerful analytical platform to investigate the preclinical efficacy of antivenoms. In this review, the studies performed on the polyvalent antivenom manufactured by Instituto Clodomiro Picado, Costa Rica, are summarized. This antivenom is prepared by immunizing horses with a mixture of the venoms of Bothrops asper, Crotalus simus and Lachesis stenophrys, and is used in Central America for the treatment of envenomings by viperid species. Overall, the antivenom shows a widespread pattern of immunological reactivity against homologous and heterologous venoms, which correlates with its ability to neutralize lethal, hemorrhagic, myotoxic, coagulant, defibrinogenating, phospholipase A2 and proteinase activities of viperid venoms. At the same time, antivenomics detected several venom components against which the antivenom shows only partial or negligible immunorecognition, such as low molecular mass vasoactive peptides, disintegrins, and some phospholipases A2, P-I metalloproteinases and serine proteinases. Such information can be used to design strategies for enhancing the antibody response of horses against poorly immunogenic, toxicologically-relevant venom components in order to further improve the efficacy of this antivenom.

BIOLOGICAL SIGNIFICANCE

The timely parenteral administration of an appropriate antivenom remains, more than a century after the development of the first serum antivenimeux by Calmette and Phisalix and Bertrand, the only currently effective treatment for snakebite envenomings. A key technical issue in the generation of novel antivenoms is the design of optimized immunization venom mixtures that ensure that the resulting antidotes will be effective against the highest number of venoms from snakes of medical concern across the geographical range where they will be used. Antivenomics is a proteomics-based protocol developed to complement in vitro and in vivo standard preclinical tests in the qualitative and quantitative characterization of the immunological profile and the extent of cross-reactivity of antivenoms against homologous and heterologous venoms. Antivenomics is translational venomics. The combination of antivenomics and neutralization assays represents a powerful analytical platform to investigate the efficacy of antivenoms at the molecular and preclinical levels. This article is part of a Special Issue entitled: Proteomics of non-model organisms.

Venomous snakes of Costa Rica: biological and medical implications of their venom proteomic profiles analyzed through the strategy of snake venomics

In spite of its small territory of ~50,000km(2), Costa Rica harbors a remarkably rich biodiversity. Its herpetofauna includes 138 species of snakes, of which sixteen pit vipers (family Viperidae, subfamily Crotalinae), five coral snakes (family Elapidae, subfamily Elapinae), and one sea snake (Family Elapidae, subfamily Hydrophiinae) pose potential hazards to human and animal health. In recent years, knowledge on the composition of snake venoms has expanded dramatically thanks to the development of increasingly fast and sensitive analytical techniques in mass spectrometry and separation science applied to protein characterization. Among several analytical strategies to determine the overall protein/peptide composition of snake venoms, the methodology known as 'snake venomics' has proven particularly well suited and informative, by providing not only a catalog of protein types/families present in a venom, but also a semi-quantitative estimation of their relative abundances. Through a collaborative research initiative between Instituto de Biomedicina de Valencia (IBV) and Instituto Clodomiro Picado (ICP), this strategy has been applied to the study of venoms of Costa Rican snakes, aiming to obtain a deeper knowledge on their composition, geographic and ontogenic variations, relationships to taxonomy, correlation with toxic activities, and discovery of novel components. The proteomic profiles of venoms from sixteen out of the 22 species within the Viperidae and Elapidae families found in Costa Rica have been reported so far, and an integrative view of these studies is hereby presented. In line with other venomic projects by research groups focusing on a wide variety of snakes around the world, these studies contribute to a deeper understanding of the biochemical basis for the diverse toxic profiles evolved by venomous snakes. In addition, these studies provide opportunities to identify novel molecules of potential pharmacological interest. Furthermore, the establishment of venom proteomic profiles offers a fundamental platform to assess the detailed immunorecognition of individual proteins/peptides by therapeutic or experimental antivenoms, an evolving methodology for which the term 'antivenomics' was coined (as described in an accompanying paper in this special issue).

BIOLOGICAL SIGNIFICANCE

Venoms represent an adaptive trait and an example of both divergent and convergent evolution. A deep understanding of the composition of venoms and of the principles governing the evolution of venomous systems is of applied importance for exploring the enormous potential of venoms as sources of chemical and pharmacological novelty but also to fight the consequences of snakebite envenomings. Key to this is the identification of evolutionary and ecological trends at different taxonomical levels. However, the evolution of venomous species and their venoms do not always follow the same course, and the identification of structural and functional convergences and divergences among venoms is often unpredictable by a phylogenetic hypothesis. Snake venomics is a proteomic-centered strategy to deconstruct the complex molecular phenotypes the venom proteomes. The proteomic profiles of venoms from sixteen out of the 22 venomous species within the Viperidae and Elapidae families found in Costa Rica have been completed so far. An integrative view of their venom composition, including the identification of geographic and ontogenic variations, is hereby presented. Venom proteomic profiles offer a fundamental platform to assess the detailed immunorecognition of individual venom components by therapeutic or experimental antivenoms. This aspect is reviewed in the companion paper. This article is part of a Special Issue entitled: Proteomics of non-model organisms.

Comparison of phylogeny, venom composition and neutralization by antivenom in diverse species of bothrops complex

In Latin America, Bothrops snakes account for most snake bites in humans, and the recommended treatment is administration of multispecific Bothrops antivenom (SAB – soro antibotrópico). However, Bothrops snakes are very diverse with regard to their venom composition, which raises the issue of which venoms should be used as immunizing antigens for the production of pan-specific Bothrops antivenoms. In this study, we simultaneously compared the composition and reactivity with SAB of venoms collected from six species of snakes, distributed in pairs from three distinct phylogenetic clades: Bothrops, Bothropoides and Rhinocerophis. We also evaluated the neutralization of Bothrops atrox venom, which is the species responsible for most snake bites in the Amazon region, but not included in the immunization antigen mixture used to produce SAB. Using mass spectrometric and chromatographic approaches, we observed a lack of similarity in protein composition between the venoms from closely related snakes and a high similarity between the venoms of phylogenetically more distant snakes, suggesting little connection between taxonomic position and venom composition. P-III snake venom metalloproteinases (SVMPs) are the most antigenic toxins in the venoms of snakes from the Bothrops complex, whereas class P-I SVMPs, snake venom serine proteinases and phospholipases A₂ reacted with antibodies in lower levels. Low molecular size toxins, such as disintegrins and bradykinin-potentiating peptides, were poorly antigenic. Toxins from the same protein family showed antigenic cross-reactivity among venoms from different species; SAB was efficient in neutralizing the B. atrox venom major toxins. Thus, we suggest that it is possible to obtain pan-specific effective antivenoms for Bothrops envenomations through immunization with venoms from only a few species of snakes, if these venoms contain protein classes that are representative of all species to which the antivenom is targeted.

Snakebite envenomation is a serious health issue in Latin America, particularly in the Amazon, where antivenom administration may be delayed due to logistic constraints. Bothrops snakes are involved in most of the snakebite-related accidents in Brazil. This work reports a comparative study of the toxin composition and antigenicity of the Bothrops venoms used to prepare the commercial antivenom and its effectiveness against the venom from Bothrops atrox, a prevalent Amazon species that is not included in the pool. Our data show a lack of connection between Bothrops taxonomic identity and venom composition. We also show that different toxins display distinct reactivity with the tested antivenom. However, the antivenom reacted similarly with each class of toxin present in the venoms of the different snakes studied. Important evidence was the neutralization of the major toxic effects of B. atrox venom, not included in the mixture of antigens used to produce the antivenom. Based on the observed antigenicity of the distinct protein classes of toxins, we suggest that it is possible to obtain pan-specific and efficient Bothrops antivenoms via immunization with venoms from a few species of snakes that are representative of the protein composition of a large number of targeted species.

Snake venomics of two poorly known Hydrophiinae: Comparative proteomics of the venoms of terrestrial Toxicocalamus longissimus and marine Hydrophis cyanocinctus

The venom proteomes of Toxicocalamus longissimus and Hydrophis cyanocinctus, a fossorial and a marine species, respectively, of the Hydrophiinae genus of Elapidae, were investigated by Edman degradation of RP-HPLC isolated proteins, and de novo MS/MS sequencing of in-gel derived tryptic peptide ions. The toxin arsenal of T. longissimus is made up of 1-2 type-I PLA(2) molecules, which account for 6.5% of the venom proteins, a minor PIII-SVMP (1.4% of the venom toxins), and ~20 members of the 3FTx family comprising 92% of the venom proteome. Seventeen proteins (5 type-I PLA(2)s and 12 3FTxs) were found in the venom of H. cyanocinctus. Three-finger toxins and type-I PLA(2) proteins comprise, respectively, 81% and 19% of its venom proteome. The simplicity of the H. cyanocinctus venom proteome is highlighted by the fact that only 6 venom components (3 short-chain neurotoxins, two long-chain neurotoxins, and one PLA(2) molecule) exhibit relative abundances >5%. As expected from its high neurotoxin abundance, the LD(50) for mice of H. cyanocinctus venom was fairly low, 0.132μg/g (intravenous) and 0.172μg/g (intraperitoneal). Our data indicate that specialization towards a lethal cocktail of 3FTx and type-I PLA(2) molecules may represent a widely adopted trophic solution throughout the evolution of Elapidae. Our results also points to a minimization of the molecular diversity of the toxin arsenal of the marine snake Hydrophis cyanocinctus in comparison to the venom proteome of its terrestrial relatives, and highlight that the same evolutionary solution, economy of the toxin arsenal, has been convergently adopted by different taxa in response to opposite selective pressures, loss and gain of neurotoxicity.

Peptide fingerprinting of the neurotoxic fractions isolated from the secretions of sea anemones Stichodactyla helianthus and Bunodosoma granulifera. New members of the APETx-like family identified by a 454 pyrosequencing approach

Sea anemones are known to contain a wide diversity of biologically active peptides, mostly unexplored according to recent peptidomic and transcriptomic studies. In the present work, the neurotoxic fractions from the exudates of Stichodactyla helianthus and Bunodosoma granulifera were analyzed by reversed-phase chromatography and mass spectrometry. The first peptide fingerprints of these sea anemones were assessed, revealing the largest number of peptide components (156) so far found in sea anemone species, as well as the richer peptide diversity of B. granulifera in relation to S. helianthus. The transcriptomic analysis of B. granulifera, performed by massive cDNA sequencing with 454 pyrosequencing approach allowed the discovery of five new APETx-like peptides (U-AITX-Bg1a-e - including the full sequences of their precursors for four of them), which together with type 1 sea anemone sodium channel toxins constitute a very distinguishable feature of studied sea anemone species belonging to genus Bunodosoma. The molecular modeling of these new APETx-like peptides showed a distribution of positively charged and aromatic residues in putative contact surfaces as observed in other animal toxins. On the other hand, they also showed variable electrostatic potentials, thus suggesting a docking onto their targeted channels in different spatial orientations. Moreover several crab paralyzing toxins (other than U-AITX-Bg1a-e), which induce a variety of symptoms in crabs, were isolated. Some of them presumably belong to new classes of crab-paralyzing peptide toxins, especially those with molecular masses below 2kDa, which represent the smallest peptide toxins found in sea anemones.

Proteomics in evolutionary ecology: linking the genotype with the phenotype

The study of the proteome (proteomics), which includes the dynamics of protein expression, regulation, interactions and its function, has played a less prominent role in evolutionary and ecological investigations in comparison with the study of the genome and transcriptome. There are, however, a number of arguments suggesting that this situation should change. First, the proteome is closer to the phenotype than the genome or the transcriptome, and as such may be more directly responsive to natural selection, and thus closely linked to adaptation. Second, there is evidence of a low correlation between protein and transcript expression levels across genes in many different organisms. Finally, there have been some recent important technological improvements in proteomics methods that make them feasible, practical and useful to address a wide range of evolutionary questions even in nonmodel organisms. The different proteomic methods, their limitations and problems when interpreting empirical data are described and discussed. In addition, the proteomic literature pertaining to evolutionary ecology is reviewed with examples, and potential applications of proteomics in a variety of evolutionary contexts are outlined. New proteomic research trends such as the study of posttranslational modifications and protein-protein interactions, as well as the combined use of the different -omics approaches, are discussed in relation to the development of a more functional and integrated perspective, needed for achieving a more comprehensive knowledge of evolutionary change.

Venom peptide analysis of Vipera ammodytes meridionalis (Viperinae) and Bothrops jararacussu (Crotalinae) demonstrates subfamily-specificity of the peptidome in the family Viperidae

Snake venom peptidomes are valuable sources of pharmacologically active compounds. We analyzed the peptidic fractions (peptides with molecular masses < 10,000 Da) of venoms of Vipera ammodytes meridionalis (Viperinae), the most toxic snake in Europe, and Bothrops jararacussu (Crotalinae), an extremely poisonous snake of South America. Liquid chromatography/mass spectrometry (LC/MS), direct infusion electrospray mass spectrometry (ESI-MS) and matrix-assisted desorption/ionization time-of-flight mass spectrometry (MALDI-TOF-MS) were applied to characterize the peptides of both snake venoms. 32 bradykinin-potentiating peptides (BPPs) were identified in the Crotalinae venom and their sequences determined. 3 metalloproteinase inhibitors, 10 BPPs and a Kunitz-type inhibitor were observed in the Viperinae venom peptidome. Variability in the C-terminus of homologous BPPs was observed, which can influence the pharmacological effects. The data obtained so far show a subfamily specificity of the venom peptidome in the Viperidae family: BPPs are the major peptide component of the Crotalinae venom peptidome lacking Kunitz-type inhibitors (with one exception) while the Viperinae venom, in addition to BPPs, can contain peptides of the bovine pancreatic trypsin inhibitor family. We found indications for a post-translational phosphorylation of serine residues in Bothrops jararacussu venom BPP (S[combining low line]QGLPPGPPIP), which could be a regulatory mechanism in their interactions with ACE, and might influence the hypotensive effect. Homology between venom BPPs from Viperidae snakes and venom natriuretic peptide precursors from Elapidae snakes suggests a structural similarity between the respective peptides from the peptidomes of both snake families. The results demonstrate that the venoms of both snakes are rich sources of peptides influencing important physiological systems such as blood pressure regulation and hemostasis. The data can be used for pharmacological and medical applications.

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

BACKGROUND

A long term research goal of venomics, of applied importance for improving current antivenom therapy, but also for drug discovery, is to understand the pharmacological potential of venoms. Individually or combined, proteomic and transcriptomic studies have demonstrated their feasibility to explore in depth the molecular diversity of venoms. In the absence of genome sequence, transcriptomes represent also valuable searchable databases for proteomic projects.

RESULTS

The venom gland transcriptomes of 8 Costa Rican taxa from 5 genera (Crotalus, Bothrops, Atropoides, Cerrophidion, and Bothriechis) of pitvipers were investigated using high-throughput 454 pyrosequencing. 100,394 out of 330,010 masked reads produced significant hits in the available databases. 5.165,220 nucleotides (8.27%) were masked by RepeatMasker, the vast majority of which corresponding to class I (retroelements) and class II (DNA transposons) mobile elements. BLAST hits included 79,991 matches to entries of the taxonomic suborder Serpentes, of which 62,433 displayed similarity to documented venom proteins. Strong discrepancies between the transcriptome-computed and the proteome-gathered toxin compositions were obvious at first sight. Although the reasons underlaying this discrepancy are elusive, since no clear trend within or between species is apparent, the data indicate that individual mRNA species may be translationally controlled in a species-dependent manner. The minimum number of genes from each toxin family transcribed into the venom gland transcriptome of each species was calculated from multiple alignments of reads matched to a full-length reference sequence of each toxin family. Reads encoding ORF regions of Kazal-type inhibitor-like proteins were uniquely found in Bothriechis schlegelii and B. lateralis transcriptomes, suggesting a genus-specific recruitment event during the early-Middle Miocene. A transcriptome-based cladogram supports the large divergence between A. mexicanus and A. picadoi, and a closer kinship between A. mexicanus and C. godmani.

CONCLUSIONS

Our comparative next-generation sequencing (NGS) analysis reveals taxon-specific trends governing the formulation of the venom arsenal. Knowledge of the venom proteome provides hints on the translation efficiency of toxin-coding transcripts, contributing thereby to a more accurate interpretation of the transcriptome. The application of NGS to the analysis of snake venom transcriptomes, may represent the tool for opening the door to systems venomics.

Ending the drought: new strategies for improving the flow of affordable, effective antivenoms in Asia and Africa

The development of snake antivenoms more than a century ago should have heralded effective treatment of the scourge of snakebite envenoming in impoverished, mostly rural populations around the world. That snakebite still exists today, as a widely untreated illness that maims, kills and terrifies men, women and children in vulnerable communities, is a cruel anachronism. Antivenom can be an effective, safe and affordable treatment for snakebites, but apathy, inaction and the politicisation of public health have marginalised both the problem (making snakebite arguably the most neglected of all neglected tropical diseases) and its solution. For lack of any coordinated approach, provision of antivenoms has been pushed off the public health agenda, leading to an incongruous decline in demand for these crucial antidotes, excused and fed by new priorities, an absence of epidemiological data, and a poor regulatory framework. These factors facilitated the infiltration of poor quality products that degrade user confidence and undermine legitimate producers. The result is that tens of thousands are denied an essential life-saving medicine, allowing a toll of human suffering that is a summation of many individual catastrophes. No strategy has been developed to address this problem and to overcome the intransigence and inaction responsible for the global tragedy of snakebite. Attempts to engage with the broader public health community through the World Health Organisation (WHO), GAVI, and other agencies have failed. Consequently, the toxinology community has taken on a leadership role in a new approach, the Global Snakebite Initiative, which seeks to mobilise the resources, skills and experience of scientists and clinicians for whom venoms, toxins, antivenoms, snakes and snakebites are already fields of interest. Proteomics is one such discipline, which has embraced the potential of using venoms in bio-discovery and systems biology. The fields of venomics and antivenomics have recently evolved from this discipline, offering fresh hope for the victims of snakebites by providing an exciting insight into the complexities, nature, fundamental properties and significance of venom constituents. Such a rational approach brings with it the potential to design new immunising mixtures from which to raise potent antivenoms with wider therapeutic ranges. This addresses a major practical limitation in antivenom use recognised since the beginning of the 20th century: the restriction of therapeutic effectiveness to the specific venom immunogen used in production. Antivenomic techniques enable the interactions between venoms and antivenoms to be examined in detail, and if combined with functional assays of specific activity and followed up by clinical trials of effectiveness and safety, can be powerful tools with which to evaluate the suitability of current and new antivenoms for meeting urgent regional needs. We propose two mechanisms through which the Global Snakebite Initiative might seek to end the antivenom drought in Africa and Asia: first by establishing a multidisciplinary, multicentre, international collaboration to evaluate currently available antivenoms against the venoms of medically important snakes from specific nations in Africa and Asia using a combination of proteomic, antivenomic and WHO-endorsed preclinical assessment protocols, to provide a validated evidence base for either recommending or rejecting individual products; and secondly by bringing the power of proteomics to bear on the design of new immunising mixtures to raise Pan-African and Pan-Asian polyvalent antivenoms of improved potency and quality. These products will be subject to rigorous clinical assessment. We propose radically to change the basis upon which antivenoms are produced and supplied for the developing world. Donor funding and strategic public health alliances will be sought to make it possible not only to sustain the financial viability of antivenom production partnerships, but also to ensure that patients are relieved of the costs of antivenom so that poverty is no longer a barrier to the treatment of this important, but grossly neglected public health emergency.