Venom Composition in a Phenotypically Variable Pit Viper ( Trimeresurus insularis) across the Lesser Sunda Archipelago

Quantification of snake venom proteomes by mass spectrometry-considerations and perspectives

The advent of soft ionization mass spectrometry-based proteomics in the 1990s led to the development of a new dimension in biology that conceptually allows for the integral analysis of whole proteomes. This transition from a reductionist to a global-integrative approach is conditioned to the capability of proteomic platforms to generate and analyze complete qualitative and quantitative proteomics data. Paradoxically, the underlying analytical technique, molecular mass spectrometry, is inherently nonquantitative. The turn of the century witnessed the development of analytical strategies to endow proteomics with the ability to quantify proteomes of model organisms in the sense of "an organism for which comprehensive molecular (genomic and/or transcriptomic) resources are available." This essay presents an overview of the strategies and the lights and shadows of the most popular quantification methods highlighting the common misuse of label-free approaches developed for model species' when applied to quantify the individual components of proteomes of nonmodel species (In this essay we use the term "non-model" organisms for species lacking comprehensive molecular (genomic and/or transcriptomic) resources, a circumstance that, as we detail in this review-essay, conditions the quantification of their proteomes.). We also point out the opportunity of combining elemental and molecular mass spectrometry systems into a hybrid instrumental configuration for the parallel identification and absolute quantification of venom proteomes. The successful application of this novel mass spectrometry configuration in snake venomics represents a proof-of-concept for a broader and more routine application of hybrid elemental/molecular mass spectrometry setups in other areas of the proteomics field, such as phosphoproteomics, metallomics, and in general in any biological process where a heteroatom (i.e., any atom other than C, H, O, N) forms integral part of its mechanism.

Assessing Target Specificity of the Small Molecule Inhibitor MARIMASTAT to Snake Venom Toxins: A Novel Application of Thermal Proteome Profiling

New treatments that circumvent the pitfalls of traditional antivenom therapies are critical to address the problem of snakebite globally. Numerous snake venom toxin inhibitors have shown promising cross-species neutralization of medically significant venom toxins in vivo and in vitro. The development of high-throughput approaches for the screening of such inhibitors could accelerate their identification, testing, and implementation and thus holds exciting potential for improving the treatments and outcomes of snakebite envenomation worldwide. Energetics-based proteomic approaches, including thermal proteome profiling and proteome integral solubility alteration (PISA) assays, represent "deep proteomics" methods for high throughput, proteome-wide identification of drug targets and ligands. In the following study, we apply thermal proteome profiling and PISA methods to characterize the interactions between venom toxin proteoforms in Crotalus atrox (Western Diamondback Rattlesnake) and the snake venom metalloprotease (SVMP) inhibitor marimastat. We investigate its venom proteome-wide effects and characterize its interactions with specific SVMP proteoforms, as well as its potential targeting of non-SVMP venom toxin families. We also compare the performance of PISA thermal window and soluble supernatant with insoluble precipitate using two inhibitor concentrations, providing the first demonstration of the utility of a sensitive high-throughput PISA-based approach to assess the direct targets of small molecule inhibitors for snake venom.

Indian green pit vipers: A lesser-known snake group of north-east India

Green pit vipers are one of the most widely distributed group of venomous snakes in south-east Asia. In Indian, green pit vipers are found in the Northern and North-eastern states spreading across eastern and central India and one of the lesser studied venoms. High morphological similarity among them has been a long-established challenge for species identification, however, a total of six species of Indian green pit viper belonging to genus Trimeresurus, Popeia and Viridovipera has been reported from North-east India. Biochemical and biological studies have revealed that venom exhibits substantial variation in protein expression level along with functional variability. The symptoms of envenomation are painful swelling at bite site, bleeding, necrosis along with systemic toxicity such as prolonged coagulopathy. Clinical data of green pit viper envenomated patients from Demow community health centre, Assam advocated against the use of Indian polyvalent antivenom pressing the need for a suitable antivenom for the treatment of green pit viper envenomation. To design effective and specific antivenom for green pit vipers, unveiling the proteome profile of these snakes is needed. In this study, a comparative venomic of green pit vipers of Northern and North-eastern India, their clinical manifestation as well as treatment protocol has been reviewed.

Proteome Decomplexation of Trimeresurus erythrurus Venom from Mizoram, India

Green pit vipers are the largest group of venomous vipers in tropical and subtropical Asia, which are responsible for most of the bite cases across this region. Among the green pit vipers of the Indian subcontinent, Trimeresurus erythrurus is the most prevalent; however, limited knowledge is available about its venomics. Proteome decomplexation of T. erythrurus venom using mass spectrometry revealed a blend of 53 different proteins/peptides belonging to 10 snake venom protein families. Phospholipase A₂ and snake venom serine proteases were found to be the major enzymatic families, and Snaclec was the major nonenzymatic family in this venom. These protein families might be responsible for consumptive coagulopathy in victims. Along with these, snake venom metalloproteases, l-amino acid oxidases, disintegrins, and cysteine-rich secretory proteins were also found, which might be responsible for inducing painful edema, tissue necrosis, blistering, and defibrination in patients. Protein belonging to C-type lectins, C-type natriuretic peptides, and glutaminyl-peptide cyclotransfreases were also observed as trace proteins. The crude venom shows platelet aggregation in the absence of any agonist, suggesting their role in alterations in platelet functions. This study is the first proteomic analysis of T. erythrurus venom, contributing an overview of different snake venom proteins/peptides responsible for various pathophysiological disorders obtained in patients. Data are available via ProteomeXchange with the identifier PXD038311.

A current perspective on snake venom composition and constituent protein families

Snake venoms are heterogeneous mixtures of proteins and peptides used for prey subjugation. With modern proteomics there has been a rapid expansion in our knowledge of snake venom composition, resulting in the venom proteomes of 30% of vipers and 17% of elapids being characterised. From the reasonably complete proteomic coverage of front-fanged snake venom composition (179 species-68 species of elapids and 111 species of vipers), the venoms of vipers and elapids contained 42 different protein families, although 18 were only reported in < 5% of snake species. Based on the mean abundance and occurrence of the 42 protein families, they can be classified into 4 dominant, 6 secondary, 14 minor, and 18 rare protein families. The dominant, secondary and minor categories account for 96% on average of a snake's venom composition. The four dominant protein families are: phospholipase A₂ (PLA₂), snake venom metalloprotease (SVMP), three-finger toxins (3FTx), and snake venom serine protease (SVSP). The six secondary protein families are: L-amino acid oxidase (LAAO), cysteine-rich secretory protein (CRiSP), C-type lectins (CTL), disintegrins (DIS), kunitz peptides (KUN), and natriuretic peptides (NP). Venom variation occurs at all taxonomic levels, including within populations. The reasons for venom variation are complex, as variation is not always associated with geographical variation in diet. The four dominant protein families appear to be the most important toxin families in human envenomation, being responsible for coagulopathy, neurotoxicity, myotoxicity and cytotoxicity. Proteomic techniques can be used to investigate the toxicological profile of a snake venom and hence identify key protein families for antivenom immunorecognition.

Diverge and Conquer: Phylogenomics of southern Wallacean forest skinks (Genus: Sphenomorphus) and their colonization of the Lesser Sunda Archipelago

The archipelagos of Wallacea extend between the Sunda and Sahul Shelves, serving as a semipermeable two-way filter influencing faunal exchange between Asia and Australo-Papua. Forest skinks (Genus Sphenomorphus) are widespread throughout southern Wallacea and exhibit complex clinal, ontogenetic, sexual, and seasonal morphological variation, rendering species delimitation difficult. We screened a mitochondrial marker for 245 Sphenomorphus specimens from this area to inform the selection of 104 samples from which we used targeted sequence capture to generate a dataset of 1154 nuclear genes (∼1.8 Mb) plus complete mitochondrial genomes. Phylogenomic analyses recovered many deeply divergent lineages, three pairs of which are now sympatric, that began to diversify in the late Miocene shortly after the oldest islands are thought to have become emergent. We infer a complex and nonstepping-stone pattern of island colonization, with the group having originated in the Sunda Arc islands before using Sumba as a springboard for colonization of the Banda Arcs. Estimates of population structure and gene flow across the region suggest total isolation except between two Pleistocene Aggregate Island Complexes that become episodically land-bridged during glacial maxima. These historical processes have resulted in at least 11 Sphenomorphus species in the region, nine of which require formal description. This fine-scale geographic partitioning of undescribed species highlights the importance of utilizing comprehensive genomic studies for defining biodiversity hotspots to be considered for conservation protection.

Venom composition of Trimeresurus albolabris, T. insularis, T. puniceus and T. purpureomaculatus from Indonesia

Background:

Several studies have been published on the characterization of Trimeresurus venoms. However, there is still limited information concerning the venom composition of Trimeresurus species distributed throughout Indonesia, which contributes to significant snakebite envenomation cases. The present study describes a comparative on the composition of T. albolabris, T. insularis, T. puniceus, and T. purpureomaculatus venoms originated from Indonesia.

Methods:

Protein content in the venom of four Trimeresurus species was determined using Bradford assay, and the venom proteome was elucidated using one-dimension SDS PAGE nano-ESI- LCMS/MS shotgun proteomics.

Results:

The venom of T. albolabris contained the highest protein content of 11.1 mg/mL, followed by T. puniceus, T. insularis and T. purpureomaculatus venom with 10.7 mg/mL, 8.9 mg/mL and 5.54 mg/mL protein, respectively. In total, our venomic analysis identified 65 proteins belonging to 16 protein families in T. purpureomaculatus; 64 proteins belonging to 18 protein families in T. albolabris; 58 different proteins belonging to 14 protein families in T. puniceus; and 48 different proteins belonging to 14 protein familiesin T. insularis. Four major proteins identified in all venoms belonged to snake venom metalloproteinase, C-type lectin, snake venom serine protease, and phospholipase A2. There were 11 common proteins in all venoms, and T. puniceus venom has the highest number of unique proteins compared to the other three venoms. Cluster analysis of the proteins and venoms showed that T. puniceus venom has the most distinct venom composition.

Conclusions:

Overall, the results highlighted venom compositional variation of four Trimeresurus spp. from Indonesia. The venoms appear to be highly similar, comprising at least four protein families that correlate with venom’s toxin properties and function. This study adds more information on venom variability among Trimeresurus species within the close geographic origin and may contribute to the development of optimum heterologous antivenom.

Potential para-specific and geographical utility of Thai Green Pit Viper (Trimeresurus albolabris) Monovalent Antivenom: Neutralization of procoagulant and hemorrhagic activities of diverse Trimeresurus pit viper venoms

The Trimeresurus complex consists of diverse medically important venomous pit vipers that cause snakebite envenomation. Antivenoms, however, are in limited supply, and are specific to only two out of the many species across Asia. This study thus investigated the immunoreactivities of regional pit viper antivenoms toward selected Trimeresurus pit viper venoms, and examined the neutralization of their hemotoxic activities. Trimeresurus albolabris Monovalent Antivenom (TaMAV, Thailand) exhibited a higher immunoreactivity than Hemato Bivalent Antivenom (HBAV, raised against Trimeresurus stejnegeri and Protobothrops mucrosquamatus, Taiwan) and Gloydius brevicaudus Monovalent Antivenom (GbMAV, China), attributed to its monovalent nature and conserved antigens in the Trimeresurus pit viper venoms. The venoms showed moderate-to-strong in vitro procoagulant and in vivo hemorrhagic effects consistent with hemotoxic envenomation, except for the Sri Lankan Trimeresurus trigonocephalus venom which lacked hemorrhagic activity. TaMAV was able to differentially neutralize both in vitro and in vivo hemotoxic effects of the venoms, with the lowest efficacy shown against the procoagulant effect of T. trigonocephalus venom. The findings suggest that TaMAV is a potentially useful treatment for envenomation caused by hetero-specific Trimeresurus pit vipers, in particular those in Southeast Asia and East Asia. Clinical study is warranted to establish its spectrum of para-specific effectiveness, and dosages need be tailored to the different species in respective regions.

Snake venom proteomics and antivenomics of two Sundaic lance-headed pit vipers: Trimeresurus wiroti (Malaysia) and Trimeresurus puniceus (Indonesia)

Envenomation by two medically important Sundaic pit vipers, Trimeresurus wiroti (Malaysia) and Trimeresurus puniceus (Indonesia), causes hemotoxic syndrome with a potentially fatal outcome. Research on the compositions and antigenicity of these pit viper venoms is however lacking, limiting our understanding of the pathophysiology and treatment of envenomation. This study investigated the venom proteomes of both species through a protein decomplexation strategy, applying C₁₈ reverse-phase high-performance liquid chromatography followed by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), and protein identification through nano-electrospray ionization liquid chromatography-tandem mass spectrometry (nano-ESI-LCMS/MS) of trypsin-digested peptides. The venom antigenicity was profiled against the Thai Green Pit Viper Antivenom (GPVAV, a hetero-specific antivenom), using indirect enzyme-linked immunosorbent assay (ELISA). The venom proteomes of T. wiroti and T. puniceus consisted of 10 and 12 toxin families, respectively. The major proteins were of diverse snake venom serine proteases (19-30% of total venom proteins), snake venom metalloproteinases (17-26%), disintegrins (9-16%), phospholipases A₂ (8-28%) and C-type lectins (~8%). These were putative snake toxins implicated in hemorrhage and coagulopathy, consistent with clinical hemotoxicity. GPVAV showed strong immunorecognition toward high and medium molecular weight proteins (e.g., SVMP and PLA₂) in both venoms, while a lower binding activity was observed toward small proteins such as disintegrins. Conserved antigenicity in the major hemotoxins supported toxicity cross-neutralization by GPVAV and indicated that the immunorecognition of low molecular weight toxins may be optimized for improved binding efficacy. Taken together, the study provides insights into the pathophysiology and antivenom treatment of envenomation caused by T. wiroti and T. puniceus in the region.

Phylogenomic Analysis Reveals Dispersal-Driven Speciation and Divergence with Gene Flow in Lesser Sunda Flying Lizards (Genus Draco)

The Lesser Sunda Archipelago offers exceptional potential as a model system for studying the dynamics of dispersal-driven diversification. The geographic proximity of the islands suggests the possibility for successful dispersal, but this is countered by the permanence of the marine barriers and extreme intervening currents that are expected to hinder gene flow. Phylogenetic and species delimitation analyses of flying lizards (genus Draco) using single mitochondrial genes, complete mitochondrial genomes, and exome-capture data sets identified 9-11 deeply divergent lineages including single-island endemics, lineages that span multiple islands, and parapatrically-distributed non-sister lineages on the larger islands. Population clustering and PCA confirmed these genetic boundaries with isolation-by-distance playing a role in some islands or island sets. While gdi estimates place most candidate species comparisons in the ambiguous zone, migration estimates suggest 9 or 10 species exist with nuclear introgression detected across some intra-island contact zones. Initial entry of Draco into the archipelago occurred at 5.5-7.5 Ma, with most inter-island colonization events having occurred between 1-3 Ma. Biogeographical model testing favors scenarios integrating geographic distance and historical island connectivity, including an initial stepping-stone dispersal process from the Greater Sunda Shelf through the Sunda Arc as far eastward as Lembata Island. However, rather than reaching the adjacent island of Pantar by dispersing over the 15-km wide Alor Strait, Draco ultimately reached Pantar (and much of the rest of the archipelago) by way of a circuitous route involving at least five over-water dispersal events. These findings suggest that historical geological and oceanographic conditions heavily influenced dispersal pathways and gene flow, which in turn drove species formation and shaped species boundaries.

The Sequence and Three-Dimensional Structure Characterization of Snake Venom Phospholipases B

Snake venom phospholipases B (SVPLBs) are the least studied enzymes. They constitute about 1% of Bothrops crude venoms, however, in other snake venoms, it is present in less than 1%. These enzymes are considered the most potent hemolytic agent in the venom. Currently, no structural information is available about these enzymes from snake venom. To better understand its three-dimensional structure and mechanisms of envenomation, the current work describes the first model-based structure report of this enzyme from Bothrops moojeni venom named as B. moojeni phospholipase B (PLB_Bm). The structure model of PLB_Bm was generated using model building software like I-TESSER, MODELLER 9v19, and Swiss-Model. The build PLB_Bm model was validated using validation tools (PROCHECK, ERRAT, and Verif3D). The analysis of the PLB_Bm modeled structure indicates that it contains 491 amino acid residues that form a well-defined four-layer αββα sandwich core and has a typical fold of the N-terminal nucleophile aminohydrolase (Ntn-hydrolase). The overall structure of PLB_Bm contains 18 β-strands and 17 α-helices with many connecting loops. The structure divides into two chains (A and B) after maturation. The A chain is smaller and contains 207 amino acid residues, whereas the B chain is larger and contains 266 amino acid residues. The sequence and structural comparison among homologous snake venom, bacterial, and mammals PLBs indicate that differences in the length and sequence composition may confer variable substrate specificity to these enzymes. Moreover, the surface charge distribution, average volume, and depth of the active site cavity also vary in these enzymes. The present work will provide more information about the structure-function relationship and mechanism of action of these enzymes in snakebite envenomation.

Comparison of Methods for Measuring Protein Concentration in Venom Samples

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

Venomics and Cellular Toxicity of Thai Pit Vipers (Trimeresurus macrops and T. hageni)

The two venomous pit vipers, Trimeresurus macrops and T. hageni, are distributed throughout Thailand, although their abundance varies among different areas. No species-specific antivenom is available for their bite victims, and the only recorded treatment method is a horse antivenom raised against T. albolabris crude venom. To facilitate assessment of the cross-reactivity of heterologous antivenoms, protein profiles of T. macrops and T. hageni venoms were explored using mass-spectrometry-based proteomics. The results show that 185 and 216 proteins were identified from T. macrops and T. hageni venoms, respectively. Two major protein components in T. macrops and T. hageni venoms were snake venom serine protease and metalloproteinase. The toxicity of the venoms on human monocytes and skin fibroblasts was analyzed, and both showed a greater cytotoxic effect on fibroblasts than monocytic cells, with toxicity occurring in a dose-dependent rather than a time-dependent manner. Exploring the protein composition of snake venom leads to a better understanding of the envenoming of prey. Moreover, knowledge of pit viper venomics facilitates the selection of the optimum heterologous antivenoms for treating bite victims.

Clinical implications of coagulotoxic variations in Mamushi (Viperidae: Gloydius) snake venoms

Snake bite is currently one of the most neglected tropical diseases affecting much of the developing world. Asian pit vipers are responsible for a considerable amount of envenomations annually and bites can cause a multitude of clinical complications resulting from coagulopathic and neuropathic effects. While intense research has been undertaken for some species of Asian pit viper, functional coagulopathic effects have been neglected for others. We investigated their effects upon the human clotting cascade using venoms of four species of Gloydius and Ovophis okinavensis, a species closely to Gloydius. All species of included within this investigation displayed varying fibrinogenolytic effects, resulting in a net anticoagulant outcome. Gloydius saxatilis and Gloydius ussuriensis displayed the most variable effects from differing localities, sampled from Russia and Korea. As this Gloydius investigation includes some geographical variation, notable results indicate key variations of these species that point to possible limitations in antivenom cross-reactivities, which may have implications for the clinical care of victims envenomed by these snakes.

Clinical implications of differential antivenom efficacy in neutralising coagulotoxicity produced by venoms from species within the arboreal viperid snake genus Trimeresurus

Snake envenomation globally is attributed to an ever-increasing human population encroaching into snake territories. Responsible for many bites in Asia is the widespread genus Trimeresurus. While bites lead to haemorrhage, only a few species have had their venoms examined in detail. We found that Trimeresurus venom causes haemorrhaging by cleaving fibrinogen in a pseudo-procoagulation manner to produce weak, unstable, short-lived fibrin clots ultimately resulting in an overall anticoagulant effect due to fibrinogen depletion. The monovalent antivenom 'Thai Red Cross Green Pit Viper antivenin', varied in efficacy ranging from excellent neutralisation of T. albolabris venom through to T. gumprechti and T. mcgregori being poorly neutralised and T. hageni being unrecognised by the antivenom. While the results showing excellent neutralisation of some non-T. albolabris venoms (such as T. flavomaculaturs, T. fucatus, and T. macrops) needs to be confirmed with in vivo tests, conversely the antivenom failure T. hageni, and the very poor results against T. gumprechti and T. mcgregori, despite being conducted in the ideal scenario of preincubation of antivenom:venom, indicates that the likelihood of clinically relevant cross-reactivity for these species is low (T. gumprechti and T. mcgregori) to non-existent (T. hageni). These same latter three species were also not inhibited by the serine protease inhibitor AEBSF, suggesting that the toxins leading to a coagulotoxic effect in these species are non-serine proteases while in contrast T. albolabris coagulotoxicity was completely impeded by AEBSF, and thus driven by kallikrein-type serine proteases. There was a conspicuous lack of phylogenetic pattern in venom variation, with the most potent venoms (T. albolabris and T. hageni) being distant to each other on the organismal tree, and with the three most divergent and poorly neutralised venoms (T. gumprechti, T. hageni, and T. mcgregori) were also not each others closest relatives. This reinforces the paradigm that the fundamental dynamic evolution of venom results in organismal phylogeny being a poor predictor of venom potency or antivenom efficacy. This study provides a robust investigation on the differential venom effects from a wide range of Trimeresurus species on coagulation, highlighting differential fibrinogenolytic effects, while also investigating the relative antivenom neutralisation capabilities of the widely available Thai Red Cross Green Pit Viper antivenom. These results therefore have immediate, real-world implications for patients envenomed by Trimeresurus species.