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

A genus-wide study on venom proteome variation and phospholipase A2 inhibition in Asian lance-headed pit vipers (genus: Trimeresurus)

High molecular weight proteins are present abundantly in viperid venoms. The amino acid sequence can be highly variable, contributing to the structure and function diversity of snake venom protein. However, this variability remains poorly understood in many species. The study investigated the venom protein variability in a distinct clade of Asian pit vipers (Trimeresurus species complex) through comparative proteomics, applying gel electrophoresis (SDS-PAGE), liquid chromatography-tandem mass spectrometry (LCMS/MS), and bioinformatic approaches. The proteomes revealed a number of conserved protein families, within each are variably expressed protein paralogs that are unrelated to the snake phylogeny and geographic origin. The expression levels of two major enzymes, i.e., snake venom serine proteinase and metalloproteinase, correlate weakly with procoagulant and hemorrhagic activities, implying co-expression of other functionally versatile toxins in the venom. The phospholipase A2 (PLA2) abundance correlates strongly with its enzymatic activity, and a unique phenotype was discovered in two species expressing extremely little PLA2. The commercial mono-specific antivenom effectively neutralized the venoms' procoagulant and hemorrhagic effects but failed to inhibit the PLA2 activities. Instead, the PLA2 activities of all venoms were effectively inhibited by the small molecule inhibitor varespladib, suggesting its potential to be repurposed as a highly potent adjuvant therapeutic in snakebite envenoming.

Deciphering toxico-proteomics of Asiatic medically significant venomous snake species: A systematic review and interactive data dashboard

Snakebite envenomation (SBE) is a neglected tropical disease (NTD) with an approximate 1.8 million cases annually. The tremendous figure is concerning, and the currently available treatment for snakebite envenomation is antivenom. However, the current antivenom has limited cross-neutralisation activity due to the variations in snake venom composition across species and geographical locations. The proteomics of medically important venomous species is essential as they study the venom compositions within and among different species. The advancement of sophisticated proteomic approaches allows intensive investigation of snake venoms. Nevertheless, there is a need to consolidate the venom proteomics profiles and distribution analysis to examine their variability patterns. This review systematically analysed the proteomics and toxicity profiles of medically important venomous species from Asia across different geographical locations. An interactive dashboard - Asiatic Proteomics Interactive Datasets was curated to consolidate the distribution patterns of the venom compositions, serve as a comprehensive directory for large-scale comparative meta-analyses. The population proteomics demonstrate higher diversities in the predominant venom toxins. Besides, inter-regional differences were also observed in Bungarus sp., Naja sp., Calliophis sp., and Ophiophagus hannah venoms. The elapid venoms are predominated with three-finger toxins (3FTXs) and phospholipase A2 (PLA2). Intra-regional variation is only significantly observed in Naja naja venoms. Proteomics diversity is more prominent in viper venoms, with widespread dominance observed in snake venom metalloproteinase (SVMP) and snake venom serine protease (SVSP). Correlations exist between the proteomics profiles and the toxicity (LD50) of the medically important venomous species. Additionally, the predominant toxins, alongside their pathophysiological effects, were highlighted and discussed as well. The insights of interactive toxico-proteomics datasets provide comprehensive frameworks of venom dynamics and contribute to developing antivenoms for snakebite envenomation. This could reduce misdiagnosis of SBE and accelerate the researchers' data mining process.

Dual therapy with phospholipase and metalloproteinase inhibitors from Sinonatrix annularis alleviated acute kidney and liver injury caused by multiple snake venoms

Snakebite envenomation often induces acute kidney injury (AKI) and acute liver injury (ALI), leading to augmented injuries and poor rehabilitation. Phospholipase A2 (PLA2) and metalloproteinase (SVMP) present in venom are responsible for the envenomation-associated events. In this study, mice envenomed with Deinagkistrodon acutus, Naja atra, or Agkistrodon halys pallas venom exhibited typical AKI and ALI symptoms, including significantly increased plasma levels of myoglobin, free hemoglobin, uric acid, aspartate aminotransferase, and alanine aminotransferase and upregulated expression of kidney NGAL and KIM-1. These effects were significantly inhibited when the mice were pretreated with natural inhibitors of PLA2 and SVMP isolated from Sinonatrix annularis (SaPLIγ and SaMPI). The inhibitors protected the physiological structural integrity of the renal tubules and glomeruli, alleviating inflammatory infiltration and diffuse hemorrhage in the liver. Furthermore, the dual therapy alleviated oxidative stress and apoptosis in the kidneys and liver by mitigating mitochondrial damage, thereby effectively reducing the lethal effect of snake venom in the inhibitor-treated mouse model. This study showed that dual therapy with inhibitors of metalloproteinase and phospholipase can effectively prevent ALI and AKI caused by snake bites. Our findings suggest that intrinsic inhibitors present in snakes are prospective therapeutic agents for multi-organ injuries caused by snake envenoming.

Proteomic characteristics of six snake venoms from the Viperidae and Elapidae families in China and their relation to local tissue necrosis

Patients envenomed by snakes from the Viperidae and Elapidae families in China often have varying degrees of local tissue necrosis. Due to the relative clinical characteristics of local tissue necrosis and ulceration following envenoming, this study has analyzed the proteome of six snake venoms from the Viperidae and Elapidae family, and the toxin profiles of each snake were compared and correlated with the clinical manifestations that follow cytotoxic envenoming. Deinagkistrodon acutus and Naja atra envenomation induce severe ulceration, which is absent in Bungarus multicinctus envenomation and mild in the other three vipers. It is interesting to note that the proportion of c-type lectins (CTL) (20.63%) in Deinagkistrodon acutus venom was relatively high, which differs from the venom of other vipers. In addition, three-fingered toxin (3FTx) (2.15%) is present in the venom of Deinagkistrodon acutus, but has not been detected in the remaining three vipers. Snake venom metalloprotease (SVMP) (34.4%-44.7%), phospholipase A2 (PLA2) (9.81%-40.83%), and snake venom serine protease (SVSP) (9.44%-16.2%) represent the most abundant families of toxin in Viperidae venom. The Elapidae venom proteome was mainly composed of neurotoxins and cytotoxins, including 3FTx (39.28%-60.08%) and PLA2 (8.24%-58.95%) toxins, however, the proportion of CRISPS (26.36%) in Naja atra venom was relatively higher compared to Bungarus multicinctus venom. Significant differences in SVMP, SVSP, and 3FTx expression levels exist between the Viperidae and the Elapidae family. The main toxins responsible for the development of tissue necrosis and ulcerations following Viperidae envenoming are hematotoxins (SVSMP, SVSP) and myotoxins (PLA2). Deinagkistrodon acutus venom contains high levels of CTL and traces of 3FTx, leading to more severe local necrosis. However, Naja atra venom can also cause severe local necrosis through the effects of myotoxin (3FTx, CRISP, PLA2). Bungarus multicinctus venom does not contain myotoxins, resulting in pure systemic neurological manifestations no obvious necrosis of local tissue in patients.

Venom Proteomics of Trimeresurus gracilis, a Taiwan-Endemic Pitviper, and Comparison of Its Venom Proteome and VEGF and CRISP Sequences with Those of the Most Related Species

Trimeresurus gracilis is an endemic alpine pitviper in Taiwan with controversial phylogeny, and its venom proteome remains unknown. In this study, we conducted a proteomic analysis of T. gracilis venom using high-performance liquid chromatography-tandem mass spectrometry and identified 155 toxin proteoforms that belong to 13 viperid venom toxin families. By searching the sequences of trypsin-digested peptides of the separated HPLC fractions against the NCBI database, T. gracilis venom was found to contain 40.3% metalloproteases (SVMPs), 15.3% serine proteases, 6.6% phospholipases A2, 5.0% L-amino acid oxidase, 4.6% Cys-rich secretory proteins (CRISPs), 3.2% disintegrins, 2.9% vascular endothelial growth factors (VEGFs), 1.9% C-type lectin-like proteins, and 20.2% of minor toxins, nontoxins, and unidentified peptides or compounds. Sixteen of these proteoforms matched the toxins whose full amino-acid sequences have been deduced from T. gracilis venom gland cDNA sequences. The hemorrhagic venom of T. gracilis appears to be especially rich in PI-class SVMPs and lacks basic phospholipase A2. We also cloned and sequenced the cDNAs encoding two CRISP and three VEGF variants from T. gracilis venom glands. Sequence alignments and comparison revealed that the PI-SVMP, kallikrein-like proteases, CRISPs, and VEGF-F of T. gracilis and Ovophis okinavensis are structurally most similar, consistent with their close phylogenetic relationship. However, the expression levels of some of their toxins were rather different, possibly due to their distinct ecological and prey conditions.

Phylogeny-Related Variations in Venomics: A Test in a Subset of Habu Snakes (Protobothrops)

We conducted a comparative analysis to unveil the divergence among venoms from a subset of Old World habu snakes (Protobothrops) in terms of venomic profiles and toxicological and enzymatic activities. A total of 14 protein families were identified in the venoms from these habu snakes, and 11 of them were shared among these venoms. The venoms of five adult habu snakes were overwhelmingly dominated by SVMP (32.56 ± 13.94%), PLA₂ (22.93 ± 9.26%), and SVSP (16.27 ± 4.79%), with a total abundance of over 65%, while the subadult P. mangshanensis had an extremely low abundance of PLA₂ (1.23%) but a high abundance of CTL (51.47%), followed by SVMP (22.06%) and SVSP (10.90%). Apparent interspecific variations in lethality and enzymatic activities were also explored in habu snake venoms, but no variations in myotoxicity were found. Except for SVSP, the resemblance of the relatives within Protobothrops in other venom traits was estimated to deviate from Brownian motion evolution based on phylogenetic signals. A comparative analysis further validated that the degree of covariation between phylogeny and venom variation is evolutionarily labile and varies among clades of closely related snakes. Our findings indicate a high level of interspecific variation in the venom proteomes of habu snakes, both in the presence or absence and the relative abundance of venom protein families, and that these venoms might have evolved under a combination of adaptive and neutral mechanisms.

Clinical features of snake envenomation in South Korea

INTRODUCTION

Three venomous snakes of the Gloydius genus belonging to the Viperidae family cause most snake envenomations in South Korea. Envenomation signs often include local swelling, coagulopathy, and rhabdomyolysis. The benefit of additional antivenom after the initial does is unclear.

METHODS

This retrospective study divided patients into four groups according to the presence of rhabdomyolysis (creatine kinase ≥1000 IU/L) and coagulopathy, which were defined using the Korean Society on Thrombosis and Hemostasis disseminated intravascular coagulation score (rhabdomyolysis, coagulopathy, combination, and local effects groups). We describe the clinical features of envenomation and the antivenom response.

RESULTS

Greater local swelling predicted more severe snakebite pain. Ninety of the 231 enrolled patients (38.9%) developed rhabdomyolysis. The patients with severe rhabdomyolysis in the combination group displayed higher peak creatine kinase activity than the rhabdomyolysis group. Seven patients with rhabdomyolysis, including two patients requiring kidney replacement therapy, developed acute kidney injury, but the incidence of acute kidney injury did not differ between the combination group and rhabdomyolysis group. Bleeding developed in 3.5% of the patients, but its incidence did not differ between the combination and coagulopathy groups. Approximately half of all patients needed repeated antivenom administration, mainly due to the local envenomation effect. Earlier administration of additional antivenom for progressive local swelling did not reduce the hospitalization duration.

CONCLUSION

Rhabdomyolysis is one of the major effects of Gloydius snake envenomation in South Korea, although it is not associated with the same risk of clinical deterioration as coagulopathy. Additionally, the ability of antivenom to ameliorate local swelling should be investigated to prevent unnecessary antivenom administration in South Korea.

Differences between Two Groups of Burmese Vipers (Viperidae: Azemiops) in the Proteomic Profiles, Immunoreactivity and Biochemical Functions of Their Venoms

Two recently revised Azemiops snakes with apparent differences in their external appearances and skeletal morphologies but unclear genetic boundaries have been proposed. Some researchers have refrained from using the newly proposed taxonomy because these two “species” might be two clades corresponding to different geographical populations of Azemiops feae. To improve the understanding of the kinship of these two Burmese viper groups, more of their characteristics should be explored in depth. We performed a comparative analysis of the proteomic profiles and biochemical activities of snake venoms from these two groups (Sichuan A. feae and Zhejiang A. feae) and evaluated the immunorecognition capacity of commercial antivenoms toward them. Eight protein families were identified in venoms from these two groups, while phospholipase B was only detected in venom from Sichuan A. feae. These protein families displayed varying degrees of differences in relative abundance between venoms, and phospholipase A₂ (Sichuan A. feae: 57.15%; Zhejiang A. feae: 65.94%) was the predominated component. Gloydius brevicaudus antivenom exhibited the strongest capacity to immunologically recognize these two venoms, but this was mainly limited to components with high molecular masses, some of which differed between venoms. Additionally, Zhejiang A. feae venom was more toxic than Sichuan A. feae venom, and the venoms expressed remarkable differences in enzymatic activities, probably resulting from the variation in the relative abundance of specific protein families. Our findings unveil differences between the two Burmese viper groups in terms of proteomic profiles, immunoreactivity, and the biochemical functions of their venoms. This information will facilitate the management of snakebites caused by these snakes.

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.

In Vitro Toxicity of Chinese Russell’s Viper (Daboia siamensis) Venom and Neutralisation by Antivenoms

Daboia siamensis (Russell’s viper) is a highly venomous and medically important snake in China, as well as much of Asia. There is minimal information on the pharmacological activity of the venom of the Chinese species, and currently no commercially available specific antivenom in China. This has led to the use of non-specific antivenoms to treat D. siamensis envenomation. In this study, the in vitro neurotoxicity and myotoxicity of D. siamensis venom was examined and the efficacy of four antivenoms was investigated, including the recently developed Chinese D. siamensis monovalent antivenom (C-DsMAV) and three commercially available antivenoms (Thai D. siamensis (Thai-DsMAV) monovalent antivenom, Deinagkistrodon acutus monovalent antivenom (DaAV), and Gloydius brevicaudus monovalent antivenom (GbAV). D. siamensis venom (10–30 µg/mL) caused the concentration-dependent inhibition of indirect twitches in the chick biventer cervicis nerve muscle preparation, without abolishing contractile responses to exogenous agonists ACh or CCh, indicating pre-synaptic neurotoxicity. Myotoxicity was also evident at these concentrations with inhibition of direct twitches, an increase in baseline tension, and the partial inhibition of ACh, CCh, and KCl responses. The prior addition of C-DsMAV or Thai-DsMAV prevented the neurotoxic and myotoxic activity of D. siamensis venom (10 µg/mL). The addition of non-specific antivenoms (GbAV and DaAV) partially prevented the neurotoxic activity of venom (10 µg/mL) but failed to neutralize the myotoxic effects. We have shown that D. siamensis venom exhibits in vitro weak presynaptic neurotoxicity and myotoxicity, which can be prevented by the pre-addition of the Chinese and Thai Russell’s viper antivenoms. Non-specific antivenoms were poorly efficacious. There should be further development of a monospecific antivenom against D. siamensis envenomation in China.

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

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

Snake Venom Proteomics, Immunoreactivity and Toxicity Neutralization Studies for the Asiatic Mountain Pit Vipers, Ovophis convictus, Ovophis tonkinensis, and Hime Habu, Ovophis okinavensis

Snakebite envenomation is a serious neglected tropical disease, and its management is often complicated by the diversity of snake venoms. In Asia, pit vipers of the Ovophis species complex are medically important venomous snakes whose venom properties have not been investigated in depth. This study characterized the venom proteomes of Ovophis convictus (West Malaysia), Ovophis tonkinensis (northern Vietnam, southern China), and Ovophis okinavensis (Okinawa, Japan) by applying liquid chromatography-tandem mass spectrometry, which detected a high abundance of snake venom serine proteases (SVSP, constituting 40–60% of total venom proteins), followed by phospholipases A₂, snake venom metalloproteinases of mainly P-III class, L-amino acid oxidases, and toxins from other protein families which were less abundant. The venoms exhibited different procoagulant activities in human plasma, with potency decreasing from O. tonkinensis > O. okinavensis > O. convictus. The procoagulant nature of venom confirms that consumptive coagulopathy underlies the pathophysiology of Ovophis pit viper envenomation. The hetero-specific antivenoms Gloydius brevicaudus monovalent antivenom (GbMAV) and Trimeresurus albolabris monovalent antivenom (TaMAV) were immunoreactive toward the venoms, and cross-neutralized their procoagulant activities, albeit at variably limited efficacy. In the absence of species-specific antivenom, these hetero-specific antivenoms may be useful in treating coagulotoxic envenomation caused by the different snakes in their respective regions.

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

Background

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

Results

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

Conclusions

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

Supplementary Information

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

Severe hypofibrinogenemia in patients bitten by Gloydius tsushimaensis in Tsushima Island, Nagasaki, Japan, and treatment strategy

Gloydius tsushimaensis is an endemic species inhabiting only Tsushima, a remote Japanese island, and is a distinct species from Gloydius blomhoffii widely distributed throughout mainland Japan and Gloydius brevicaudus and Gloydius ussuriensis which are geographically distributed in South Korea. This is the first multicenter retrospective study of G. tsushimaensis bites in Japan. A study of seventy-two patients who visited the former Izuhara Hospital, the former Naka Tsushima Hospital, Tsushima Hospital, and Kamitsushima Hospital during the fourteen years from January 1, 2005, to December 31, 2018, revealed the typical clinical characteristics of G. tsushimaensis bites. Five out of seventy-two cases (6.9%) showed severe hypofibrinogenemia, in which fibrinogen levels were below 100 mg/dl, which is an unreported clinical finding for G. blomhoffii bites. Generally, when fibrinogen levels are lower than 100 mg/dl, the bleeding risk increases, and it is perilous. Severe hypofibrinogenemia cases did not improve after G. blomhoffii antivenom administration. Additionally, all five cases had disseminated intravascular coagulation, and there were two cases of acute kidney injury and one death. five cases had a median maximum creatine kinase level of 5171 IU/l (Interquartile range: 4992-41,310). Although the mechanism is not precise, coagulation tests showed that the G. tsushimaensis venom contains a thrombin-like enzyme. Based on this research, we created an algorithm for the treatment of G. tsushimaensis bites and unified the treatment methods used on the island.

Development of an LC-MS multivariate nontargeted methodology for differential analysis of the peptide profile of Asian hornet venom (Vespa velutina nigrithorax): application to the investigation of the impact of collection period variation

Insect venom is a highly complex mixture of bioactive compounds, containing proteins, peptides, and small molecules. Environmental factors can alter the venom composition and lead to intraspecific variation in its bioactivity properties. The investigation of discriminating compounds caused by variation impacts can be a key to manage sampling and explore the bioactive compounds. The present study reports the development of a peptidomic methodology based on UHPLC-ESI-QTOF-HRMS analysis followed by a nontargeted multivariate analysis to reveal the profile variance of Vespa velutina venom collected in different conditions. The reliability of the approach was enhanced by optimizing certain XCMS data processing parameters and determining the sample peak threshold to eliminate the interfering features. This approach demonstrated a good repeatability and a criterion coefficient of variation (CV) > 30% was set for deleting nonrepeatable features from the matrix. The methodology was then applied to investigate the impact of collection period variation. PCA and PLS-DA models were used and validated by cross-validation and permutation tests. A slight discrimination was found between winter and summer hornet venom in two successive years with 10 common discriminating compounds. Graphical abstract.

Venomics of Trimeresurus (Popeia) nebularis, the Cameron Highlands Pit Viper from Malaysia: Insights into Venom Proteome, Toxicity and Neutralization of Antivenom

Trimeresurus nebularis is a montane pit viper that causes bites and envenomation to various communities in the central highland region of Malaysia, in particular Cameron’s Highlands. To unravel the venom composition of this species, the venom proteins were digested by trypsin and subjected to nano-liquid chromatography-tandem mass spectrometry (LC-MS/MS) for proteomic profiling. Snake venom metalloproteinases (SVMP) dominated the venom proteome by 48.42% of total venom proteins, with a characteristic distribution of P-III: P-II classes in a ratio of 2:1, while P-I class was undetected. Snaclecs constituted the second most venomous protein family (19.43%), followed by snake venom serine proteases (SVSP, 14.27%), phospholipases A₂ (5.40%), disintegrins (5.26%) and minor proteins including cysteine-rich secretory proteins, L-amino acid oxidases, phosphodiesterases, 5′-nucleotidases. The venomic profile correlates with local (painful progressive edema) and systemic (hemorrhage, coagulopathy, thrombocytopenia) manifestation of T. nebularis envenoming. As specific antivenom is unavailable for T. nebularis, the hetero-specific Thai Green Pit viper Monovalent Antivenom (GPVAV) was examined for immunological cross-reactivity. GPVAV exhibited good immunoreactivity to T. nebularis venom and the antivenom effectively cross-neutralized the hemotoxic and lethal effects of T. nebularis (lethality neutralizing potency = 1.6 mg venom per mL antivenom). The findings supported GPVAV use in treating T. nebularis envenoming.

Venom proteomics and antivenom neutralization for the Chinese eastern Russell's viper, Daboia siamensis from Guangxi and Taiwan

The eastern Russell's viper (Daboia siamensis) causes primarily hemotoxic envenomation. Applying shotgun proteomic approach, the present study unveiled the protein complexity and geographical variation of eastern D. siamensis venoms originated from Guangxi and Taiwan. The snake venoms from the two geographical locales shared comparable expression of major proteins notwithstanding variability in their toxin proteoforms. More than 90% of total venom proteins belong to the toxin families of Kunitz-type serine protease inhibitor, phospholipase A2, C-type lectin/lectin-like protein, serine protease and metalloproteinase. Daboia siamensis Monovalent Antivenom produced in Taiwan (DsMAV-Taiwan) was immunoreactive toward the Guangxi D. siamensis venom, and effectively neutralized the venom lethality at a potency of 1.41 mg venom per ml antivenom. This was corroborated by the antivenom effective neutralization against the venom procoagulant (ED = 0.044 ± 0.002 µl, 2.03 ± 0.12 mg/ml) and hemorrhagic (ED50 = 0.871 ± 0.159 µl, 7.85 ± 3.70 mg/ml) effects. The hetero-specific Chinese pit viper antivenoms i.e. Deinagkistrodon acutus Monovalent Antivenom and Gloydius brevicaudus Monovalent Antivenom showed negligible immunoreactivity and poor neutralization against the Guangxi D. siamensis venom. The findings suggest the need for improving treatment of D. siamensis envenomation in the region through the production and the use of appropriate antivenom.

Toxin-resolved antivenomics-guided assessment of the immunorecognition landscape of antivenoms

Snakebite envenoming represents a major issue in rural areas of tropical and subtropical regions across sub-Saharan Africa, South to Southeast Asia, Latin America and Oceania. Antivenoms constitute the only scientifically validated therapy for snakebite envenomings, provided they are safe, effective, affordable, accessible and administered appropriately. However, the lack of financial incentives in a technology that has remained relatively unchanged for more than a century, has contributed to some manufacturers leaving the market and others downscaling production or increasing the prices, leading to a decline in the availability and accessibility for these life-saving antidotes to millions of rural poor most at risk from snakebites in low income countries. The shortage of antivenoms can be significantly alleviated by optimizing the use of current antivenoms (through the assessment of their specific and paraspecific efficacy against the different medically relevant homologous and heterologous snake venoms) and by generating novel polyspecific antivenoms exhibiting broad clinical spectrum and wide geographic distribution range. Research on venoms has been continuously enhanced by advances in technology. Particularly, the last decade has witnessed the development of omics strategies for unravelling the toxin composition of venoms ("venomics") and to assess the immunorecognition profile of antivenoms ("antivenomics"). Here, we review recent developments and reflect on near future innovations that promise to revolutionize the mutually enlightening relationship between evolutionary and translational venomics.

Comparative analyses of putative toxin gene homologs from an Old World viper, Daboia russelii

Availability of snake genome sequences has opened up exciting areas of research on comparative genomics and gene diversity. One of the challenges in studying snake genomes is the acquisition of biological material from live animals, especially from the venomous ones, making the process cumbersome and time-consuming. Here, we report comparative sequence analyses of putative toxin gene homologs from Russell’s viper (Daboia russelii) using whole-genome sequencing data obtained from shed skin. When compared with the major venom proteins in Russell’s viper studied previously, we found 45–100% sequence similarity between the venom proteins and their putative homologs in the skin. Additionally, comparative analyses of 20 putative toxin gene family homologs provided evidence of unique sequence motifs in nerve growth factor (NGF), platelet derived growth factor (PDGF), Kunitz/Bovine pancreatic trypsin inhibitor (Kunitz BPTI), cysteine-rich secretory proteins, antigen 5, andpathogenesis-related1 proteins (CAP) and cysteine-rich secretory protein (CRISP). In those derived proteins, we identified V11 and T35 in the NGF domain; F23 and A29 in the PDGF domain; N69, K2 and A5 in the CAP domain; and Q17 in the CRISP domain to be responsible for differences in the largest pockets across the protein domain structures in crotalines, viperines and elapids from the in silico structure-based analysis. Similarly, residues F10, Y11 and E20 appear to play an important role in the protein structures across the kunitz protein domain of viperids and elapids. Our study highlights the usefulness of shed skin in obtaining good quality high-molecular weight DNA for comparative genomic studies, and provides evidence towards the unique features and evolution of putative venom gene homologs in vipers.

Tracing the evolution of venom phospholipases A2 in Gloydius strauchii and related pitvipers: A tale of two acidic isozymes

Two acidic Asp⁴⁹-PLA₂s with Glu⁶ substitution and a neutral Lys⁴⁹-PLA (designated Gst-K49) were cloned from G. strauchii venom glands, their full amino acid sequences were deduced. The predominant acidic PLA₂ (designated Gst-E6a) contains 124 residues and the M¹⁸W³⁰ substitutions, while the minor acidic PLA₂ (designated Gst-E6b) contains 122 residues and the V¹⁸A³⁰ substitutions. Their sequences are most similar to those of the respective orthologous PLA₂s of G. intermedius venom. Gst-E6a and Gst-E6b appear to be paralogs and possibly have different predatory targets or functions. The LC-MS/MS results indicate the presence of only three PLA₂ gene products in the crude venom, the relative expression levels were in the order of Gst-E6a ≫ Gst-E6b > Gst-K49, as confirmed by qPCR results. In contrast to other Gloydius, G. strauchii venom does not contain neurotoxic or basic anticoagulant Asp⁴⁹-PLA₂s, but Gst-K49 is the first Lys⁴⁹-PLA₂ identified in Gloydius venoms. However, its venom content is relatively low and its pI value 7.3 is much lower than those of other Lys⁴⁹-PLA₂s and. The Lys⁴⁹-PLA₂ genes appear to regress in the venom of most of Gloydius and related rattlesnake, and this evolutionary regression occurred before the dispersal of Asian pitvipers to the New World.

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.

Preparation of monoclonal antibodies against gamma-type phospholipase A2 inhibitors and immunodetection of these proteins in snake blood

Background

The gamma-type phospholipase A₂ inhibitor (PLIγ) is a natural protein commonly found in snake serum, which can neutralize pathophysiological effects of snake venom phospholipases A₂. Therefore, this protein is a potential candidate to the development of a novel antivenom. To the best of our knowledge, there is no antibody currently available for PLIγ identification and characterization.

Methods

Bioinformatics prediction of epitope using DNAStar software was performed based on the sequence of Sinonatrix annularis PLIγ (SaPLIγ). The best epitope ¹⁵¹CPVLRLSNRTHEANRNDLIKVA¹⁷² was chosen and synthesized, and then conjugated to keyhole limpet hemocyanin and bovine serum albumin for use as an immunogen and plate-coating antigen, respectively.

Results

Eighteen IgG anti-PLIγ mAb hybridoma cell strains were obtained, and all the mAbs had positive interaction with recombinant His6-PLIγ and natural SaPLIγ. Moreover, the mAb from 10E9 strain was also successfully used for the immunodetection of other snake serum PLIγs. cDNA sequence alignment of those PLIγs from different snake species showed that their epitope segments were highly homologous.

Conclusions

The successful preparation of anti-PLIγmAb is significant for further investigation on the relationship between the structure and function of PLIγs, as well as the interaction between PLIγs and PLA₂s.

Preclinical Evaluation of the Efficacy of Antivenoms for Snakebite Envenoming: State-of-the-Art and Challenges Ahead

Animal-derived antivenoms constitute the mainstay in the therapy of snakebite envenoming. The efficacy of antivenoms to neutralize toxicity of medically-relevant snake venoms has to be demonstrated through meticulous preclinical testing before their introduction into the clinical setting. The gold standard in the preclinical assessment and quality control of antivenoms is the neutralization of venom-induced lethality. In addition, depending on the pathophysiological profile of snake venoms, the neutralization of other toxic activities has to be evaluated, such as hemorrhagic, myotoxic, edema-forming, dermonecrotic, in vitro coagulant, and defibrinogenating effects. There is a need to develop laboratory assays to evaluate neutralization of other relevant venom activities. The concept of the 3Rs (Replacement, Reduction, and Refinement) in Toxinology is of utmost importance, and some advances have been performed in their implementation. A significant leap forward in the study of the immunological reactivity of antivenoms against venoms has been the development of “antivenomics”, which brings the analytical power of mass spectrometry to the evaluation of antivenoms. International partnerships are required to assess the preclinical efficacy of antivenoms against snake venoms in different regions of the world in order to have a detailed knowledge on the neutralizing profile of these immunotherapeutics.

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

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

Venomics of Tropidolaemus wagleri, the sexually dimorphic temple pit viper: Unveiling a deeply conserved atypical toxin arsenal

Tropidolaemus wagleri (temple pit viper) is a medically important snake in Southeast Asia. It displays distinct sexual dimorphism and prey specificity, however its venomics and inter-sex venom variation have not been thoroughly investigated. Applying reverse-phase HPLC, we demonstrated that the venom profiles were not significantly affected by sex and geographical locality (Peninsular Malaya, insular Penang, insular Sumatra) of the snakes. Essentially, venoms of both sexes share comparable intravenous median lethal dose (LD₅₀) (0.56-0.63 μg/g) and cause neurotoxic envenomation in mice. LCMS/MS identified six waglerin forms as the predominant lethal principles, comprising 38.2% of total venom proteins. Fourteen other toxin-protein families identified include phospholipase A₂, serine proteinase, snaclec and metalloproteinase. In mice, HPLC fractions containing these proteins showed insignificant contribution to the overall venom lethality. Besides, the unique elution pattern of approximately 34.5% of non-lethal, low molecular mass proteins (3-5 kDa) on HPLC could be potential biomarker for this primitive crotalid species. Together, the study unveiled the venom proteome of T. wagleri that is atypical among many pit vipers as it comprises abundant neurotoxic peptides (waglerins) but little hemotoxic proteinases. The findings also revealed that the venom is relatively well conserved intraspecifically despite the drastic morphological differences between sexes.

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.

Proteomic comparisons of venoms of long-term captive and recently wild-caught Eastern brown snakes (Pseudonaja textilis) indicate venom does not change due to captivity

Snake venom is a highly variable phenotypic character, and its variation and rapid evolution are important because of human health implications. Because much snake antivenom is produced from captive animals, understanding the effects of captivity on venom composition is important. Here, we have evaluated toxin profiles from six long-term (LT) captive and six recently wild-caught (RC) eastern brown snakes, Pseudonaja textilis, utilizing gel electrophoresis, HPLC-MS, and shotgun proteomics. We identified proteins belonging to the three-finger toxins, group C prothrombin activators, Kunitz-type serine protease inhibitors, and phospholipases A2, among others. Although crude venom HPLC analysis showed LT snakes to be higher in some small molecular weight toxins, presence/absence patterns showed no correlation with time in captivity. Shotgun proteomics indicated the presence of similar toxin families among individuals but with variation in protein species. Although no venom sample contained all the phospholipase A2 subunits that form the textilotoxin, all did contain both prothrombin activator subunits. This study indicates that captivity has limited effects on venom composition, that venom variation is high, and that venom composition may be correlated to geographic distribution.

BIOLOGICAL SIGNIFICANCE

Through proteomic comparisons, we show that protein variation within LT and RC groups of snakes (Pseudonaja textilis) is high, thereby resulting in no discernible differences in venom composition between groups. We utilize complementary techniques to characterize the venom proteomes of 12 individual snakes from our study area, and indicate that individuals captured close to one another have more similar venom gel electrophoresis patterns than those captured at more distant locations. These data are important for understanding natural variation in and potential effects of captivity on venom composition.

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

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

BIOLOGICAL SIGNIFICANCE

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

Advances in venomics

The term “venomics” was coined to describe the global study of venom and venom glands, targeting comprehensive characterization of the whole toxin profile of a venomous animal by means of proteomics, transcriptomics, genomics and bioinformatics studies.

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.