Snake venomics and toxicological profiling of the arboreal pitviper Bothriechis supraciliaris from Costa Rica

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.

Individual variations in the protein profiles and functional activities of the eyelash palm pit-viper (Bothriechis schlegelii) venom from the Colombian southwest region

Bothriechis schlegelii is a venomous snake found in Central and South America, mainly sighted in regions devoted to agriculture. However, in Colombia, little is known about its contribution to the total envenoming cases. Furthermore, there are no reports of the biochemical and functional activities of venoms from the southwest populations, and the differences respecting other populations are unknown. This study analyzed the protein profiles of venom samples obtained from three specimens originating from this region of Colombia using electrophoresis and chromatography. The lethality, edema-induction, hemorrhagic, defibrinating, coagulant, and indirect hemolytic activities were also evaluated. As a result, venoms were composed of proteins with a wide range of molecular weights, most of them below <37 kDa, with differences between male and female electrophoretic and chromatographic profiles. These variations were also observed in the evaluation of venom functional activities such as pro-coagulant, indirect hemolytic, and edema-inducing activities, whereas neither hemorrhagic nor defibrinating activities were detected. These results are also different considering reports with venom samples from other geographical locations, restating the existence of high intraspecific variability in B. schlegelii venoms, which could have relevant pathophysiological and therapeutic implications.

Boyer L, D. Possani L, López-Casillas F, Olvera A, Romero C, Zamudio F, Alagón A. Neotropical Rattlesnake (Crotalus simus) Venom Pharmacokinetics in Lymph and Blood Using an Ovine Model

The most abundant protein families in viper venoms are Snake Venom Metalloproteases (SVMPs), Snake Venom Serine Proteases (SVSPs) and Phospholipases (PLA₂s). These are primarily responsible for the pathophysiology caused by the bite of pit-vipers; however, there are few studies that analyze the pharmacokinetics (PK) of whole venom (WV) and its protein families. We studied the pathophysiology, PK profile and differential absorption of representative toxins from venom of Neotropical Rattlesnake (Crotalus simus) in a large animal model (ovine). Toxins studied included crotoxin (the main lethal component), which causes moderate to severe neurotoxicity; SVSPs, which deplete fibrinogen; and SVMPs, which cause local tissue damage and local and systemic hemorrhage. We found that Whole Venom (WV) was highly bioavailable (86%) 60 h following intramuscular (IM) injection, and extrapolation suggests that bioavailability may be as high as 92%. PK profiles of individual toxins were consistent with their physicochemical properties and expected clinical effects. Lymph cannulated animals absorbed 1.9% of WV through lymph during the first 12 h. Crotoxin was minimally detectable in serum after intravenous (IV) injection; however, following IM injection it was detected in lymph but not in blood. This suggests that crotoxin is quickly released from the blood toward its tissue targets.

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).

Investigation of the estuarine stonefish (Synanceia horrida) venom composition

The Estuarine stonefish (Synanceia horrida) is recognised as one of the most venomous fish species in the world but the overall venom composition has yet to be investigated using in-depth transcriptomic and proteomic methods. To date, known venom components are restricted to a hyaluronidase and a large, pore-forming toxin known as Stonustoxin (SNTX). Transcriptomic sequencing of the venom gland resulted in over 170,000 contigs with only 0.4% that were homologous to putative venom proteins. Integration of the transcriptomic data with proteomic data from the S. horrida venom confirmed the hyaluronidase and SNTX to be present, together with several other protein families including major contributions from C-type lectins. Other protein families observed included peroxiredoxin and several minor protein families such as Golgi-associated plant pathogenesis related proteins, tissue pathway factor inhibitors, and Kazal-type serine protease inhibitors that, although not putative venom proteins, may contribute to the venom's adverse effects. BIOLOGICAL SIGNIFICANCE: Proteomic analysis of milked Synanceia horrida venom, paired with transcriptomic analysis of the venom gland tissue revealed for the first time the composition of one of the world's most dangerous fish venoms. The results demonstrate that the venom is relatively less complex compared to other well-studied venomous animals with a number of unique proteins not previously found in animal venoms.

Comparative venom-gland transcriptomics and venom proteomics of four Sidewinder Rattlesnake (Crotalus cerastes) lineages reveal little differential expression despite individual variation

Changes in gene expression can rapidly influence adaptive traits in the early stages of lineage diversification. Venom is an adaptive trait comprised of numerous toxins used for prey capture and defense. Snake venoms can vary widely between conspecific populations, but the influence of lineage diversification on such compositional differences are unknown. To explore venom differentiation in the early stages of lineage diversification, we used RNA-seq and mass spectrometry to characterize Sidewinder Rattlesnake (Crotalus cerastes) venom. We generated the first venom-gland transcriptomes and complementary venom proteomes for eight individuals collected across the United States and tested for expression differences across life history traits and between subspecific, mitochondrial, and phylotranscriptomic hypotheses. Sidewinder venom was comprised primarily of hemorrhagic toxins, with few cases of differential expression attributable to life history or lineage hypotheses. However, phylotranscriptomic lineage comparisons more than doubled instances of significant expression differences compared to all other factors. Nevertheless, only 6.4% of toxins were differentially expressed overall, suggesting that shallow divergence has not led to major changes in Sidewinder venom composition. Our results demonstrate the need for consensus venom-gland transcriptomes based on multiple individuals and highlight the potential for discrepancies in differential expression between different phylogenetic hypotheses.

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.

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.

Isolation and characterization of four medium-size disintegrins from the venoms of Central American viperid snakes of the genera Atropoides, Bothrops, Cerrophidion and Crotalus

Four disintegrins were isolated from the venoms of the Central American viperid snakes Atropoides mexicanus (atropoimin), Bothrops asper (bothrasperin), Cerrophidion sasai (sasaimin), and Crotalus simus (simusmin). Purifications were performed by reverse-phase HPLC. The four disintegrins have biochemical characteristics, i.e. molecular mass and location of Cys, which allow their classification within the group of medium-size disintegrins. All of them present the canonical RGD sequence, which determines their interaction with integrins in cell membranes. The disintegrins inhibited ADP and collagen-induced human platelet aggregation, with similar IC50s in the nM range. In addition, disintegrins inhibited the adhesion of an endothelial cell line and a melanoma cell line to the extracellular matrix proteins type I collagen, laminin, fibronectin, and vitronectin, albeit showing variable ability to exert this activity. This study expands the inventory of this family of viperid venom proteins, and reports, for the first time, disintegrins from the venoms of species of the genera Atropoides and Cerrophidion.

Venomics and antivenomics of Bothrops erythromelas from five geographic populations within the Caatinga ecoregion of northeastern Brazil

The Caatinga lancehead, Bothrops erythromelas, is a medically relevant species, responsible for most of the snakebite accidents in most parts of its distribution range in northeastern Brazil. The spectrum and geographic variability of its venom toxins were investigated applying a venomics approach to venom pools from five geographic areas within the Caatinga ecoregion. Despite its wide habitat, populations of B. erythromelas from Ceará, Pernambuco, Juazeiro, Paraiba, and Ilha de Itaparica exhibit highly conserved venom proteomes. Mirroring their compositional conservation, the five geographic venom pools also showed qualitatively and quantitatively overlapping antivenomic profiles against antivenoms generated in Vital Brazil (BR) and Clodomiro Picado (CR) Institutes, using different venoms in the immunization mixtures. The paraspecificity exhibited by the Brazilian SAB and the Costa Rican BCL antivenoms against venom toxins from B. erythromelas indicates large immunoreactive epitope conservation across genus Bothrops during the last ~14 million years, thus offering promise for the possibility of generating a broad-spectrum bothropic antivenom. Biological Significance Accidental snakebite envenomings represent an important public health hazard in Brazil. Ninety per cent of the yearly estimated 20-30,000 snakebite accidents are caused by species of the Bothrops genus. Bothrops erythromelas, a small, moderately stocky terrestrial venomous snake, is responsible for most of the snakebite accidents in its broad distribution range in the Caatinga, a large ecoregion in northeastern Brazil. To gain a deeper insight into the spectrum of medically important toxins present in the venom of the Caatinga lancehead, we applied a venomics approach to define the proteome and geographic variability of adult B. erythromelas venoms from five geographic regions. Although intraspecific compositional variation between venoms among specimens from different geographic regions has long been appreciated by herpetologists and toxinologists as a general feature of highly adaptable and widely distributed snake species, the five B. erythromelas populations investigated exhibit highly conserved venom proteomes. The overall toxin profile of the Caatinga lancehead's venom explains the local and systemic effects observed in envenomations by B. erythromelas. The five geographic venom pools sampled also showed qualitatively and quantitatively overlapping antivenomic profiles against antivenoms generated using different bothropic venoms in the immunization mixtures. The large immunoreactive epitope conservation across genus Bothrops offers promise for the generation of a broad-spectrum bothropic antivenom.

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.

Proteomic and functional profiling of the venom of Bothrops ayerbei from Cauca, Colombia, reveals striking interspecific variation with Bothrops asper venom

Bothrops ayerbei, a pitviper inhabiting the Patía River's basin (Valle Alto del Río Patía) in the Southwestern Department of Cauca, Colombia, was considered as a variant form of Bothrops asper prior to being proposed as a new species in 2010, on the basis of subtle morphological differences. This study reports the proteomic and functional profiling of B. ayerbei venom. Its most striking feature is an almost complete absence (0.7%) of phospholipases A2 (PLA2), which is in contrast to the high proportion of these enzymes (25.3%) in the venom of B. asper from Cauca, as well as in other species of Bothrops. The predominant proteins in B. ayerbei venom are metalloproteinases (53.7%), in agreement with its higher hemorrhagic and lethal activities compared to B. asper venom. Moreover, the negligible content of PLA2s in B. ayerbei venom correlates with its weaker myotoxic effect, in contrast to B. asper venom, here shown to contain abundant Asp49- and Lys49-type PLA2s responsible for its strong myotoxic activity. Other components identified in B. ayerbei venom include bradykinin-potentiating-like peptides and proteins belonging to the C-type lectin/lectin-like, serine proteinase, l-amino acid oxidase, disintegrin, cysteine-rich secretory protein, nerve growth factor, and phosphodiesterase families. The venom composition of B. ayerbei resembles that of neonate specimens of B. asper, which shows a predominance of metalloproteinases, with only low amounts of PLA2s. Therefore, the present findings suggest that the expression of venom proteins in B. ayerbei, in contrast to B. asper, might retain a marked 'paedomorphic' condition. Altogether, the proteomic and toxicological characterization of the venom of B. ayerbei here reported argues in favor of its taxonomical separation from B. asper in Cauca, Colombia.

BIOLOGICAL SIGNIFICANCE

B. ayerbei, a pitviper found in Cauca, Colombia, had been considered as a variant form of B. asper, but was recently described as a new species on the basis of subtle morphological differences. Our study provides the first detailed proteomic and functional analysis of the venom of B. ayerbei, revealing striking interspecific variation from B. asper, thus arguing in favor of their taxonomical separation. In addition, the observed composition of the venom of B. ayerbei correlates well with its functional and toxicological properties, helping to predict the main clinical manifestations in envenomings by this species, which inflicts a considerable number of snakebites in the Southwestern regions of Colombia.