Combined venomics, antivenomics and venom gland transcriptome analysis of the monocoled cobra ( Naja kaouthia ) from China

Analysis of Naja kaouthia snake venom composition and in-vitro enzymatic activities of 29 specimens in captivity: Highlighting the importance of individual variation in venom pool production

Naja kaouthia is a medically important snake, widely distributed throughout Southeast Asia, with a diverse venom composition. N. kaouthia venom is subject to significant intraspecific variation, caused by several factors, such as the wide geographic distribution of the species, sexual and ontogenetic factors. However, individual variation is a factor that has only been studied with small sample size groups and/or with pooled samples. With this in mind, this study evaluates the composition and in-vitro enzymatic activities of 29 individual venom samples from specimens born in captivity, with a similar genetic background caused by inbreeding, using SDS-PAGE under reducing conditions, RP-HPLC profiles and enzymatic activities of PLA2, LAAO and proteolytic activity over azocasein. Even in this scenario, we were able to observe significant variations in abundance and activity of PLA2. Individual variations in LAAO activity, as well as a sexual dimorphism in which males present a significantly higher LAAO activity than females were observed. Phosphodiasterase and CRiSP abundance were also found and considered to have multiple effects in the clinical manifestations of envenomation by presenting synergistic effects with other proteins from N. kaouthia venom. The RP-HPLC profiles were better at detecting compositional differences than SDS-PAGE profiles and better correlated with enzymatic activities, being a better technique to screen variation profiles and reinforcing the importance of individual venom analysis prior to pooling.

Exploring the Venom Gland Transcriptome of Bothrops asper and Bothrops jararaca: De Novo Assembly and Analysis of Novel Toxic Proteins

Previous proteomic studies of viperid venom revealed that it is mainly composed of metalloproteinases (SVMPs), serine proteinases (SVSPs), phospholipase A2 (PLA2), and C-type lectins (CTLs). However, other proteins appear in minor amounts that affect prey and need to be identified. This study aimed to identify novel toxic proteins in the venom gland transcriptome of Bothrops asper and Bothrops jararaca, using data from NCBI. Bioinformatics tools were used to assemble, identify, and compare potentially novel proteins in both species, and we performed functional annotation with BLASTX against the NR database. While previous assemblies have been performed for B. jararaca, this is the first assembly of the B. asper venom gland transcriptome. Proteins with potentially novel functions were identified, including arylsulfatase and dihydroorotate dehydrogenase, among others, that could have implications for venom toxicity. These results suggest that the identified proteins may contribute to venom toxic variation and provide new opportunities for antivenom research. The study improves the understanding of the protein composition of Bothrops venom and suggests new possibilities for the development of treatments and antivenoms.

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.

Composition, pharmacology, and pathophysiology of the venom of monocled cobra (Naja kaouthia)- a medically crucial venomous snake of southeast Asia: An updated review

The Monocled Cobra (Naja kaouthia), a category one medically significant snake from the Elapidae family, inflicts severe envenomation in South and Southeast Asian countries. N. kaouthia is distributed throughout the eastern and northeastern parts of India, Nepal, Bangladesh, Myanmar, Thailand, Vietnam, Malaysia, and southwestern China. Envenomation by N. kaouthia is a medical emergency, and the primary clinical symptoms are neurotoxicity and localized tissue destruction. Unfortunately, data on the actual magnitude of N. kaouthia envenomation is scarce due to poor record keeping, lack of diagnostic kits, and region-wise well-coordinated epidemiological surveys. The present review highlights the diversity in the composition of N. Kaouthia venom (NKV) across various geographical regions, as revealed through biochemical and proteomic analyses. The qualitative and quantitative differences in the toxin isoforms result in differences in lethality and pathophysiological manifestation that may limit the effectiveness of antivenom therapy. Studies on commercial polyvalent antivenom (PAV) effectiveness against distinct NKV samples have revealed varying toxicity and enzymatic activity neutralization. Additionally, the identification of snake venom's poorly immunogenic toxins by mass spectrometry, quantification of venom-specific antibodies, and implications for antivenom therapy against snakebites are highlighted. Future directions involve clinical studies on NK envenomation where the snake is frequently encountered and the correlation of this data with NKV composition in that region. For more efficient and superior hospital management of NK envenomation, research should enhance the current immunization procedure to boost the development of antibodies against less immunogenic venom components of this snake.

Three finger toxins of elapids: structure, function, clinical applications and its inhibitors

The WHO lists snakebite as a "neglected tropical disease". In tropical and subtropical areas, envenoming is an important public health issue. This review article describes the structure, function, chemical composition, natural inhibitors, and clinical applications of Elapids' Three Finger Toxins (3FTX) using scientific research data. The primary venomous substance belonging to Elapidae is 3FTX, that targets nAChR. Three parallel β-sheets combine to create 3FTX, which has four or five disulfide bonds. The three primary types of 3FTX are short-chain, long-chain, and nonconventional 3FTX. The functions of 3FTX depend on the specific toxin subtype and the target receptor or ion channel. The well-known effect of 3FTX is probably neurotoxicity because of the severe consequences of muscular paralysis and respiratory failure in snakebite victims. 3FTX have also been studied for their potential clinical applications. α-bungarotoxin has been used as a molecular probe to study the structure and function of nAChRs (Nicotinic Acetylcholine Receptors). Acid-sensing ion channel (ASIC) isoforms 1a and 1b are inhibited by Mambalgins, derived from Black mamba venom, which hinders their function and provide an analgesic effect. α- Cobra toxin is a neurotoxin purified from Chinese cobra (Naja atra) binds to nAChR at the neuronal junction and causes an analgesic effect for moderate to severe pain. Some of the plants and their compounds have been shown to inhibit the activity of 3FTX, and their mechanisms of action are discussed.

Proteomic Analysis of Heloderma horridum horridum Venom: Assessment to Its Transcriptome and Newfound Proteins

Heloderma horridum horridum, a venomous reptile native to America, has a venom with potential applications in treating type II diabetes. In this work, H. h. horridum venom was extracted, lyophilized, and characterized using enzymatic assays for hyaluronidase, phospholipase, and protease. Proteomic analysis of the venom was conducted employing bottom-up/shotgun approaches, SDS-PAGE, high-pH reversed-phase chromatography, and fractionation of tryptic peptides using nano-LC-MS/MS. The proteins found in H. h. horridum venom were reviewed according to the classification of the transcriptome previously reported. The proteomic approach identified 101 enzymes, 36 other proteins, 15 protein inhibitors, 11 host defense proteins, and 1 toxin, including novel venom components such as calcium-binding proteins, phospholipase A2 inhibitors, serpins, cathepsin, subtilases, carboxypeptidase-like, aminopeptidases, glycoside hydrolases, thioredoxin transferases, acid ceramidase-like, enolase, multicopper oxidases, phosphoglucose isomerase (PGI), fructose-1,6-bisphosphatase class 1, pentraxin-related, peptidylglycine α-hydroxylating monooxygenase/peptidyl-hydroxyglycine α-amidating lyase, carbonic anhydrase, acetylcholinesterase, dipeptidylpeptidase, and lysozymes. These findings contribute to understanding the venomous nature of H. h. horridum and highlight its potential as a source of bioactive compounds. Data are available via PRoteomeXchange with the identifier PXD052417.

Selection Across the Three-Dimensional Structure of Venom Proteins from North American Scolopendromorph Centipedes

Gene duplication followed by nucleotide differentiation is one of the simplest mechanisms to develop new functions for genes. However, the evolutionary processes underlying the divergence of multigene families remain controversial. We used multigene families found within the diversity of toxic proteins in centipede venom to test two hypotheses related to venom evolution: the two-speed mode of venom evolution and the rapid accumulation of variation in exposed residues (RAVER) model. The two-speed mode of venom evolution proposes that different types of selection impact ancient and younger venomous lineages with negative selection being the predominant form in ancient lineages and positive selection being the dominant form in younger lineages. The RAVER hypothesis proposes that, instead of different types of selection acting on different ages of venomous lineages, the different types of selection will selectively contribute to amino acid variation based on whether the residue is exposed to the solvent where it can potentially interact directly with toxin targets. This hypothesis parallels the longstanding understanding of protein evolution that suggests that residues found within the structural or active regions of the protein will be under negative or purifying selection, and residues that do not form part of these areas will be more prone to positive selection. To test these two hypotheses, we compared the venom of 26 centipedes from the order Scolopendromorpha from six currently recognized species from across North America using both transcriptomics and proteomics. We first estimated their phylogenetic relationships and uncovered paraphyly among the genus Scolopendra and evidence for cryptic diversity among currently recognized species. Using our phylogeny, we then characterized the diverse venom components from across the identified clades using a combination of transcriptomics and proteomics. We conducted selection-based analyses in the context of predicted three-dimensional properties of the venom proteins and found support for both hypotheses. Consistent with the two-speed hypothesis, we found a prevalence of negative selection across all proteins. Consistent with the RAVER hypothesis, we found evidence of positive selection on solvent-exposed residues, with structural and less-exposed residues showing stronger signal for negative selection. Through the use of phylogenetics, transcriptomics, proteomics, and selection-based analyses, we were able to describe the evolution of venom from an ancient venomous lineage and support principles of protein evolution that directly relate to multigene family evolution.

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.

Proteome of monocled cobra (Naja kaouthia) venom and potent anti breast cancer peptide from trypsin hydrolyzate of the venom protein

Anticancer peptide is one of the target in the development of new anticancer drug. Bioactive peptide can be originated from isolated free peptide or produced by hydrolysis of protein. Protein is the main component of Naja kaouthia venom, and due to the toxicity of the venom, it can be assessed as the source of anticancer peptides. This study aims to characterize the venom protein and to identify peptides from the snake venom of N. kaouthia as anticancer. Proteome analysis was employed trypsin hydrolysis of N. kaouthia venom protein completed with HRMS analysis protein database query. Preparative tryptic hydrolysis of the protein followed by reverse-phased fractionation and anti breast cancer activity testing were performed to identify the potent anticancer from the hydrolysate. Proteomic analysis by high-resolution mass spectrometry revealed that there are 20 enzymatic and non-enzymatic proteins in N. kaouthia venom. The 25% methanol peptide fraction had the most active anticancer activity against MCF-7 breast cancer cells and showed promising selectivity (selectivity index = 12.87). Amino acid sequences of eight peptides were identified as potentially providing anticancer compounds. Molecular docking analysis showed that WWSDHR and IWDTIEK peptides gave specific interactions and better binding affinity energy with values of -9.3 kcal/mol and -8.4 kcal/mol, respectively. This study revealed peptides from the snake venom of N. kaouthia became a potent source of new anticancer agents.

Multilevel Comparison of Indian Naja Venoms and Their Cross-Reactivity with Indian Polyvalent Antivenoms

Snake envenoming is caused by many biological species, rather than a single infectious agent, each with a multiplicity of toxins in their venom. Hence, developing effective treatments is challenging, especially in biodiverse and biogeographically complex countries such as India. The present study represents the first genus-wide proteomics analysis of venom composition across Naja species (N. naja, N. oxiana, and N. kaouthia) found in mainland India. Venom proteomes were consistent between individuals from the same localities in terms of the toxin families present, but not in the relative abundance of those in the venom. There appears to be more compositional variation among N. naja from different locations than among N. kaouthia. Immunoblotting and in vitro neutralization assays indicated cross-reactivity with Indian polyvalent antivenom, in which antibodies raised against N. naja are present. However, we observed ineffective neutralization of PLA2 activities of N. naja venoms from locations distant from the source of immunizing venoms. Antivenom immunoprofiling by antivenomics revealed differential antigenicity of venoms from N. kaouthia and N. oxiana, and poor reactivity towards 3FTxs and PLA2s. Moreover, there was considerable variation between antivenoms from different manufacturers. These data indicate that improvements to antivenom manufacturing in India are highly desirable.

Recurrent neurotoxity in Naja kaouthia envenomation: A case report from Assam, India

A 35 year old, male patient, bitten by Naja kaouthia with mild pain was admitted in Demow Government Community Health Centre. After 90 min post bite he developed neurotoxic symptoms. As per standard protocol, the patient was treated with 25 vials of antivenom and two doses of glycopyrrolate and neostigmine. Subsequently, he was seemingly devoid of any neurotoxic symptoms and showed signs of recovery. However, after 70 h, the neurotoxic symptoms recurred, and the patient was again treated with an additional 10 vials of ASV along with one dose of glycopyrrolate and neostigmine. Subsequently, the patient recovered completely from all the other symptoms of envenomation. This is the first report of recurrence of neurotoxic symptoms in a patient envenomed by Naja kaouthia in Assam, India and supports the need for greater attention and careful documentation of management of snakebite in the region.

Current Insights in the Mechanisms of Cobra Venom Cytotoxins and Their Complexes in Inducing Toxicity: Implications in Antivenom Therapy

Cytotoxins (CTXs), an essential class of the non-enzymatic three-finger toxin family, are ubiquitously present in cobra venoms. These low-molecular-mass toxins, contributing to about 40 to 60% of the cobra venom proteome, play a significant role in cobra venom-induced toxicity, more prominently in dermonecrosis. Structurally, CTXs contain the conserved three-finger hydrophobic loops; however, they also exhibit a certain degree of structural diversity that dictates their biological activities. In their mechanism, CTXs mediate toxicity by affecting cell membrane structures and membrane-bound proteins and activating apoptotic and necrotic cell death pathways. Notably, some CTXs are also responsible for depolarizing neurons and heart muscle membranes, thereby contributing to the cardiac failure frequently observed in cobra-envenomed victims. Consequently, they are also known as cardiotoxins (CdTx). Studies have shown that cobra venom CTXs form cognate complexes with other components that potentiate the toxic effects of the venom's individual component. This review focuses on the pharmacological mechanism of cobra venom CTXs and their complexes, highlighting their significance in cobra venom-induced pathophysiology and toxicity. Furthermore, the potency of commercial antivenoms in reversing the adverse effects of cobra venom CTXs and their complexes in envenomed victims has also been discussed.

Venom Variation of Neonate and Adult Chinese Cobras in Captivity Concerning Their Foraging Strategies

The venom and transcriptome profile of the captive Chinese cobra (Naja atra) is not characterized until now. Here, LC-MS/MS and illumine technology were used to unveil the venom and trascriptome of neonates and adults N. atra specimens. In captive Chinese cobra, 98 co-existing transcripts for venom-related proteins was contained. A total of 127 proteins belong to 21 protein families were found in the profile of venom. The main components of snake venom were three finger toxins (3-FTx), snake venom metalloproteinase (SVMP), cysteine-rich secretory protein (CRISP), cobra venom factor (CVF), and phosphodiesterase (PDE). During the ontogenesis of captive Chinese cobra, the rearrangement of snake venom composition occurred and with obscure gender difference. CVF, 3-FTx, PDE, phospholipase A2 (PLA2) in adults were more abundant than neonates, while SVMP and CRISP in the neonates was richer than the adults. Ontogenetic changes in the proteome of Chinese cobra venom reveals different strategies for handling prey. The levels of different types of toxin families were dramatically altered in the wild and captive specimens. Therefore, we speculate that the captive process could reshape the snake venom composition vigorously. The clear comprehension of the composition of Chinese cobra venom facilitates the understanding of the mechanism of snakebite intoxication and guides the preparation and administration of traditional antivenom and next-generation drugs for snakebite.

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.

The myth of cobra venom cytotoxin: More than just direct cytolytic actions

Cobra venom cytotoxin (CTX) is a non-enzymatic three-finger toxin that constitutes 40-60% of cobra venom. Thus, it plays an important role in the pathophysiology of cobra envenomation, especially in local dermonecrosis. The three-finger hydrophobic loops of CTX determine the cytotoxicity. Nevertheless, the actual mechanisms of cytotoxicity are not fully elucidated as they involve not only cytolytic actions but also intracellular signalling-mediated cell death pathways. Furthermore, the possible transition cell death pattern remains to be explored. The actual molecular mechanisms require further studies to unveil the relationship between different CTXs from different cobra species and cell types which may result in differential cell death patterns. Here, we discuss the biophysical interaction of CTX with the cell membrane involving four binding modes: electrostatic interaction, hydrophobic partitioning, isotropic phase, and oligomerisation. Oligomerisation of CTX causes pore formation in the membrane lipid bilayer. Additionally, the CTX-induced apoptotic pathway can be executed via death receptor-mediated extrinsic pathways and mitochondrial-mediated intrinsic pathways. We also discuss lysosomal-mediated necrosis and the occurrence of necroptosis following CTX action. Collectively, we provided an insight into concentration-dependent transition of cell death pattern which involves different mechanistic actions. This contributes a new direction for further investigation of cytotoxic pathways activated by the CTXs for future development of biotherapeutics targeting pathological effects caused by CTX.

Preclinical efficacy testing of three antivenoms against Naja ashei venom-induced lethality

This study aimed to determine the efficacy of Inoserp, Vins bioproducts, and South African Institute of Medical Research (SAIMR) polyvalent antivenoms in neutralizing Naja ashei venom-induced lethality in mice. The neutralization efficacy of the antivenoms were expressed as effective dose, median effective ratio, potency, normalized potency, volume, and the number of vials of antivenom required to neutralize 100 mg of Naja ashei venom (NAV).

Snake Venom Proteomics of Samar Cobra (Naja samarensis) from the Southern Philippines: Short Alpha-Neurotoxins as the Dominant Lethal Component Weakly Cross-Neutralized by the Philippine Cobra Antivenom

The Samar Cobra, Naja samarensis, is endemic to the southern Philippines and is a WHO-listed Category 1 venomous snake species of medical importance. Envenomation caused by N. samarensis results in neurotoxicity, while there is no species-specific antivenom available for its treatment. The composition and neutralization of N. samarensis venom remain largely unknown to date. This study thus aimed to investigate the venom proteome of N. samarensis for a comprehensive profiling of the venom composition, and to examine the immunorecognition as well as neutralization of its toxins by a hetero-specific antivenom. Applying C₁₈ reverse-phase high-performance liquid chromatography (RP-HPLC) and tandem mass spectrometry (LC-MS/MS), three-finger toxins (3FTx) were shown to dominate the venom proteome by 90.48% of total venom proteins. Other proteins in the venom comprised snake venom metalloproteinases, phospholipases A_2, cysteine-rich secretory proteins, venom nerve growth factors, L-amino acid oxidases and vespryn, which were present at much lower abundances. Among all, short-chain alpha-neurotoxins (SαNTX) were the most highly expressed toxin within 3FTx family, constituting 65.87% of the total venom proteins. The SαNTX is the sole neurotoxic component of the venom and has an intravenous median lethal dose (LD₅₀) of 0.18 μg/g in mice. The high abundance and low LD₅₀ support the potent lethal activity of N. samarensis venom. The hetero-specific antivenom, Philippine Cobra Antivenom (PCAV, raised against Naja philippinensis) were immunoreactive toward the venom and its protein fractions, including the principal SαNTX. In efficacy study, PCAV was able to cross-neutralize the lethality of SαNTX albeit the effect was weak with a low potency of 0.20 mg/ml (defined as the amount of toxin completely neutralized per milliliter of the antivenom). With a volume of 5 ml, each vial of PCAV may cross-neutralize approximately 1 mg of the toxin in vivo. The findings support the potential para-specific use of PCAV in treating envenomation caused by N. samarensis while underscoring the need to improve the potency of its neutralization activity, especially against the highly lethal alpha-neurotoxins.

Compositional and toxicological investigation of pooled venom from farm-raised Naja atra

Background:

Naja atra is a venomous snake species medically relevant in China. In the current study, we evaluated the composition and toxicological profile of venom collected from farm-raised N. atra.

Methods:

Venom was collected from third-generation captive bred N. atra on a snake farm in Hunan Province, China. The venom was analyzed using sodium dodecyl sulfate polyacrylamide gel electrophoresis and nano-liquid chromatography with electrospray ionization tandem mass spectrometry. In addition, hemolytic activity, median lethal dose, serum biochemical and histopathological parameters were accessed.

Results:

N. atra venom proteome was dominated by phospholipase A₂ (46.5%) and three-finger toxins (41.4 %), and a set of common low relative abundance proteins, including cysteine-rich secretory proteins (4.7%), NGF-beta (2.4%), snake venom metalloproteinase (1.5%), glutathione peroxidase (0.6%), vespryn (0.3%), and 5ʹ-nucleotidases (0.2%) were also found. Furthermore, the venom exhibited direct hemolytic activity, neurotoxicity, myotoxicity, and high lethal potency in mice, with a subcutaneous median lethal dose of 1.02 mg/kg. Histopathological analysis and serum biochemical tests revealed that venom caused acute hepatic, pulmonary and renal injury in mice.

Conclusion:

This study revealed the composition and toxicity of venom collected from farm-raised N. atra, thereby providing a reference for the analysis of venom samples collected from captive-born venomous snakes in the future.

A strategy for the enrichment and characterization of disulfide bond-contained proteins from Chinese cobra (Naja atra) venom

A new strategy combined gold-coated magnetic nanocomposites assisted enrichment with mass spectrometry was developed for the characterization of disulfide bond-contained proteins from Chinese cobra (Naja atra) venom. In this work, core-shell nanocomposites were synthesized by the seed-mediated growth method and used for the enrichment of snake venom proteins containing disulfide bonds. A total of 3545 tryptic digested peptides derived from 96 venom proteins in Naja atra venom were identified. The venom proteins comprised 14 toxin families including three-finger toxins, phospholipase A₂ , snake venom metalloproteinase, cobra venom factor, and so forth. Extra 16 venom proteins were detected exclusively in the nanocomposites set, among which 11 venom proteins were from the three-finger toxins family. In the present study, the proposed simple and efficient protocol replaced the tedious and laborious technologies commonly used for pre-separating crude snake venom, suggesting widely implementation in low-abundance or trace disulfide bond-contained proteins or peptides characterization.

Elucidating the Venom Diversity in Sri Lankan Spectacled Cobra (Naja naja) through De Novo Venom Gland Transcriptomics, Venom Proteomics and Toxicity Neutralization

Inadequate effectiveness of Indian antivenoms in treating envenomation caused by the Spectacled Cobra/Indian Cobra (Naja naja) in Sri Lanka has been attributed to geographical variations in the venom composition. This study investigated the de novo venom-gland transcriptomics and venom proteomics of the Sri Lankan N. naja (NN-SL) to elucidate its toxin gene diversity and venom variability. The neutralization efficacy of a commonly used Indian antivenom product in Sri Lanka was examined against the lethality induced by NN-SL venom in mice. The transcriptomic study revealed high expression of 22 toxin genes families in NN-SL, constituting 46.55% of total transcript abundance. Three-finger toxins (3FTX) were the most diversely and abundantly expressed (87.54% of toxin gene expression), consistent with the dominance of 3FTX in the venom proteome (72.19% of total venom proteins). The 3FTX were predominantly S-type cytotoxins/cardiotoxins (CTX) and α-neurotoxins of long-chain or short-chain subtypes (α-NTX). CTX and α-NTX are implicated in local tissue necrosis and fatal neuromuscular paralysis, respectively, in envenomation caused by NN-SL. Intra-species variations in the toxin gene sequences and expression levels were apparent between NN-SL and other geographical specimens of N. naja, suggesting potential antigenic diversity that impacts antivenom effectiveness. This was demonstrated by limited potency (0.74 mg venom/ml antivenom) of the Indian polyvalent antivenom (VPAV) in neutralizing the NN-SL venom. A pan-regional antivenom with improved efficacy to treat N. naja envenomation is needed.

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.

Utilisation of compounds from venoms in drug discovery

Difficult drug targets are becoming the normal course of business in drug discovery, sometimes due to large interacting surfaces or only small differences in selectivity regions. For these, a different approach is merited: compounds lying somewhere between the small molecule and the large antibody in terms of many properties including stability, biodistribution and pharmacokinetics. Venoms have evolved over millions of years to be complex mixtures of stable molecules derived from other somatic molecules, the stability comes from the pressure to be ready for delivery at a moment's notice. Snakes, spiders, scorpions, jellyfish, wasps, fish and even mammals have evolved independent venom systems with complex mixtures in their chemical arsenal. These venom-derived molecules have been proven to be useful tools, such as for the development of antihypotensive angiotensin converting enzyme (ACE) inhibitors and have also made successful drugs such as Byetta® (Exenatide), Integrilin® (Eptifibatide) and Echistatin. Only a small percentage of the available chemical space from venoms has been investigated so far and this is growing. In a new era of biological therapeutics, venom peptides present opportunities for larger target engagement surface with greater stability than antibodies or human peptides. There are challenges for oral absorption and target engagement, but there are venom structures that overcome these and thus provide substrate for engineering novel molecules that combine all desired properties. Venom researchers are characterising new venoms, species, and functions all the time, these provide great substrate for solving the challenges presented by today's difficult targets.

Engineering of single-domain antibodies for next-generation snakebite antivenoms

Given the magnitude of the global snakebite crisis, strategies to ensure the quality of antivenom, as well as the availability and sustainability of its supply are under development by several research groups. Recombinant DNA technology has allowed the engineering of monoclonal antibodies and recombinant fragments as alternatives to conventional antivenoms. Besides having higher therapeutic efficacy, with broad neutralization capacity against local and systemic toxicity, novel antivenoms need to be safe and cost-effective. Due to the biological and physical chemical properties of camelid single-domain antibodies, with high volume of distribution to distal tissue, their modular format, and their versatility, their biotechnological application has grown considerably in recent decades. This article presents the most up-to-date developments concerning camelid single-domain-based antibodies against major toxins from snake venoms, the main venomous animals responsible for reported envenoming cases and related human deaths. A brief discussion on the composition, challenges, and perspectives of antivenoms is presented, as well as the road ahead for next-generation antivenoms based on single-domain antibodies.

Snake venom vascular endothelial growth factors (svVEGFs): Unravelling their molecular structure, functions, and research potential

Vascular endothelial growth factor (VEGF) is a key regulator of angiogenesis, a physiological process characterized by the formation of new vessels from a preexisting endothelium. VEGF has also been implicated in pathologic states, such as neoplasias, intraocular neovascular disorders, among other conditions. VEGFs are distributed in seven different families: VEGF-A, B, C, D, and PIGF (placental growth factor), which are identified in mammals; VEGF-E, which are encountered in viruses; and VEGF-F or svVEGF (snake venom VEGF) described in snake venoms. This is the pioneer review of svVEGF family, exploring its distribution among the snake venoms, molecular structure, main functions, and potential applications.

Remarkable intrapopulation venom variability in the monocellate cobra (Naja kaouthia) unveils neglected aspects of India's snakebite problem

Interpopulation venom variation has been widely documented in snakes across large geographical distances. This variability is known to markedly influence the effectiveness of snakebite therapy, as antivenoms manufactured against one population may not be effective against others. In contrast, the extent of intrapopulation venom variability, especially at finer geographical scales, remains largely uninvestigated. Moreover, given the historical focus on the 'big four' Indian snakes, our understanding of venom variation in medically important yet neglected snakes, such as the monocellate cobra (Naja kaouthia), remains unclear. To address this shortcoming, we investigated N. kaouthia venoms sampled across a small spatial scale (<50 km) in Eastern India. An interdisciplinary approach employed in this study unveiled considerable intrapopulation differences in the venom proteomic composition, pharmacological and biochemical activities, and toxicity profiles. Documentation of stark differences in venoms at such a finer geographical scale, despite the influence of similar ecological and environmental conditions, is intriguing. Furthermore, evaluation of in vitro and in vivo venom recognition and neutralisation potential of Indian polyvalent 'big four' antivenoms and Thai monovalent N. kaouthia antivenom revealed concerning deficiencies. These results highlight the negative impact of phylogenetic divergence and intrapopulation snake venom variation on the effectiveness of conventional antivenom therapy. SIGNIFICANCE: In contrast to our understanding of snake venom variation across large distances, which is theorised to be shaped by disparities in ecology and environment, intrapopulation variation at finer geographic scales remains scarcely investigated. Assessment of intrapopulation venom variability in Naja kaouthia at a small spatial scale (<50 km) in Eastern India unravelled considerable differences in venom compositions, activities and potencies. While the influence of subtle differences in prey preference and local adaptations cannot be ruled out, these findings, perhaps, also emphasise the role of accelerated molecular evolutionary regimes that rapidly introduce variations in evolutionarily younger lineages, such as advanced snakes. The inability of 'big four' Indian antivenoms and Thai N. kaouthia monovalent antivenom in countering these variations highlights the importance of phylogenetic considerations for the development of efficacious snakebite therapy. Thus, we provide valuable insights into the venoms of one of the most medically important yet neglected Indian snakes.

Convergent evolution of pain-inducing defensive venom components in spitting cobras

Convergent evolution provides insights into the selective drivers underlying evolutionary change. Snake venoms, with a direct genetic basis and clearly defined functional phenotype, provide a model system for exploring the repeated evolution of adaptations. While snakes use venom primarily for predation, and venom composition often reflects diet specificity, three lineages of cobras have independently evolved the ability to spit venom at adversaries. Using gene, protein, and functional analyses, we show that the three spitting lineages possess venoms characterized by an up-regulation of phospholipase A2 (PLA2) toxins, which potentiate the action of preexisting venom cytotoxins to activate mammalian sensory neurons and cause enhanced pain. These repeated independent changes provide a fascinating example of convergent evolution across multiple phenotypic levels driven by selection for defense.

Cytotoxicity of Snake Venoms and Cytotoxins From Two Southeast Asian Cobras (Naja sumatrana, Naja kaouthia): Exploration of Anticancer Potential, Selectivity, and Cell Death Mechanism

Venoms of cobras (Naja spp.) contain high abundances of cytotoxins, which contribute to tissue necrosis in cobra envenomation. The tissue-necrotizing activity of cobra cytotoxins, nevertheless, indicates anticancer potentials. This study set to explore the anticancer properties of the venoms and cytotoxins from Naja sumatrana (equatorial spitting cobra) and Naja kaouthia (monocled cobra), two highly venomous species in Southeast Asia. The cytotoxicity, selectivity, and cell death mechanisms of their venoms and cytotoxins (NS-CTX from N. sumatrana: NS-CTX; N. kaouthia: NK-CTX) were elucidated in human lung (A549), prostate (PC-3), and breast (MCF-7) cancer cell lines. Cytotoxins were purified through a sequential fractionation approach using cation-exchange chromatography, followed by C₁₈ reverse-phase high-performance liquid chromatography (HPLC) to homogeneity validated with sodium dodecyl sulfate-polyacrylamide gel electrophoresis, and identified by liquid chromatography-tandem mass spectrometry (LCMS/MS). The cobra venoms and their respective cytotoxins exhibited concentration-dependent growth inhibitory effects in all cell lines tested, with the cytotoxins being more potent compared to the corresponding whole venoms. NS-CTX and NK-CTX are, respectively, P-type and S-type isoforms of cytotoxin, based on the amino acid sequences as per LCMS/MS analysis. Both cytotoxins exhibited differential cytotoxic effects in the cell lines tested, with NS-CTX (P-type cytotoxin) being significantly more potent in inhibiting the growth of the cancer cells. Both cytotoxins demonstrated promising selectivity only for the A549 lung cancer cell line (selectivity index = 2.17 and 2.26, respectively) but not in prostate (PC-3) and breast (MCF-7) cancer cell lines (selectivity index < 1). Flow cytometry revealed that the A549 lung cancer cells treated with NS-CTX and NK-CTX underwent necrosis predominantly. Meanwhile, the cytotoxins induced mainly caspase-independent late apoptosis in the prostate (PC-3) and breast (MCF-7) cancer cells lines but lacked selectivity. The findings revealed the limitations and challenges that could be faced during the development of new cancer therapy from cobra cytotoxins, notwithstanding their potent anticancer effects. Further studies should aim to overcome these impediments to unleash the anticancer potentials of the cytotoxins.

A comparative analysis of the proteomes and biological activities of the venoms from two sea snakes, Hydrophis curtus and Hydrophis cyanocinctus, from Hainan, China

We characterized and compared the venom protein profiles of Hydrophis curtus (synonyms: Lapemis hardwickii, Lapemis curtus and Hydrophis hardwickii) and Hydrophis cyanocinctus, the two representatives of medically important venomous sea snakes in Chinese waters using proteomic approaches. A total of 47 and 38 putative toxins were identified in H. curtus venom (HcuV) and H. cyanocinctus venom (HcyV), respectively, and these toxins could be grouped into 15 functional categories, mainly proteinases, phospholipases, three-finger toxins (3FTxs), lectins, protease inhibitors, ion channel inhibitors, cysteine-rich venom proteins (CRVPs) and snake venom metalloproteases (SVMPs). The constituent ratio of each toxin category varied between HcuV and HcyV with 3FTx (54% in HcuV/69% in HcyV) and PLA2 (38% in HcuV/22% in HcyV) unanimously ranked as the top two most abundant families. Both HcuV and HcyV exhibited relatively high lethality (LD₅₀ values in mice of 0.34 μg/g and 0.24 μg/g, respectively), specific PLA2 activity and hemolytic activity. On the basis of several previous reports of HcuV and HcyV collected from other areas, these findings greatly expand our understanding of geographical variation and interspecies diversity of the two sea snake venoms and can provide a scientific basis for the development of specific sea snake antivenom in the future.

Preparation and detection of sea snake antisera raised in rabbits

Antivenoms are currently the most effective medication used in the treatment of snakebites. However, there were relatively few studies on preparation of antivenoms targeting sea snakes, especially common sea snakes in China. In this study, we sought to prepare and detect mono- and bispecific antisera raised in rabbits against venoms of two sea snakes, Hydrophis cyanocinctus and H. curtus. The results of enzyme-linked immunosorbent assay showed that the rabbit antisera generally showed clearly detectable immunological cross-reactions after the third immunization and indicated that the strength of cross-reactions increased with an increase in the immunizing dose. Proteins within the H. cyanocinctus and H. curtus venoms showed similar profiles and were mainly concentrated in the low-molecular-weight region (8-25 kDa). Western blotting results revealed that the bands of these low-molecular weight proteins were dense and showed strong immunogenicity. Although we detected comparatively few bands of the high-molecular-weight proteins, these also showed strong immunogenicity. Our results indicate that both mono- and bispecific antisera both can neutralize H. cyanocinctus and H. curtus venoms, and in this regard, the monospecific H. curtus and bispecific antiserum were found to be superior to the H. cyanocinctus antiserum. Given the increasing frequency of snakebites worldwide, we believe that the findings of this study will have high practical applicability.

Naja kaouthia venom protein, Nk-CRISP, upregulates inflammatory gene expression in human macrophages

Cysteine-Rich Secretory Proteins (CRISP) are widespread in snake venoms and known to target ion channels. More recently, CRISPs have been shown to mediate inflammatory responses. Involvement of potential receptor in CRISP-induced inflammatory reactions, however, remains unknown. A CRISP protein named as Nk-CRISP, was isolated from the venom of Naja kaouthia. The molecular mass of the purified protein was found to be ~25 kDa and the primary sequence was determined by MALDI TOF-TOF. The involvement of this protein in proinflammatory effects were evaluated in THP-1 macrophage-like cells. Nk-CRISP treated cells induced up-regulation of several inflammatory marker genes in dose dependent manner. Toll like receptor 4 (TLR4)-myeloid differentiation factor 2 (MD2) complex are known to play crucial role in recognition of damage/pathogen-associated molecular patterns and activation of innate immune response. Therefore, we hypothesized that snake venom CRISP could also modulate the innate immune response via TLR4-MD2 complex. In-silico molecular docking study of cobra CRISP with TLR4-MD2 receptor complex reveals CRISP engages its cysteine-rich domain (CRD) to interact with complex. Inhibition of TLR4 signalling pathway using CLI-095 confirmed the role of TLR4 in Nk-CRISP induced inflammatory responses. Collectively, these findings imply that TLR4 initiates proinflammatory signalling following recognition of cobra CRISP and alteration of TLR4 receptor might improve or control CRISP induced inflammation.

The origin and diversification of a novel protein family in venomous snakes

The genetic origins of novelty are a central interest of evolutionary biology. Most new proteins evolve from preexisting proteins but the evolutionary path from ancestral gene to novel protein is challenging to trace, and therefore the requirements for and order of coding sequence changes, expression changes, or gene duplication are not clear. Snake venoms are important novel traits that are comprised of toxins derived from several distinct protein families, but the genomic and evolutionary origins of most venom components are not understood. Here, we have traced the origin and diversification of one prominent family, the snake venom metalloproteinases (SVMPs) that play key roles in subduing prey in many vipers. Genomic analyses of several rattlesnake (Crotalus) species revealed the SVMP family massively expanded from a single, deeply conserved adam28 disintegrin and metalloproteinase gene, to as many as 31 tandem genes in the Western Diamondback rattlesnake (Crotalus atrox) through a number of single gene and multigene duplication events. Furthermore, we identified a series of stepwise intragenic deletions that occurred at different times in the course of gene family expansion and gave rise to the three major classes of secreted SVMP toxins by sequential removal of a membrane-tethering domain, the cysteine-rich domain, and a disintegrin domain, respectively. Finally, we show that gene deletion has further shaped the SVMP complex within rattlesnakes, creating both fusion genes and substantially reduced gene complexes. These results indicate that gene duplication and intragenic deletion played essential roles in the origin and diversification of these novel biochemical weapons.

Anticoagulant activity of krait, coral snake, and cobra neurotoxic venoms with diverse proteomes are inhibited by carbon monoxide

BACKGROUND

A phenomena of interest is the in vitro anticoagulant effects of neurotoxins found in elapid venoms that kill by paralysis. These enzymes include phospholipase A2 (PLA2), and it has recently been demonstrated that carbon monoxide inhibits the PLA2-dependent neurotoxin contained in Mojave rattlesnake type A venom. The purpose of this investigation was to assess if the anticoagulant activity of elapid venoms containing PLA2 and/or three finger toxins could be inhibited by carbon monoxide.

METHODS

Venoms collected from Bungarus multicinctus, Micrurus fulvius, and five Naja species were exposed to carbon monoxide via carbon monoxide releasing molecule-2 prior to placement into human plasma. Coagulation kinetics were assessed via thrombelastography.

RESULTS

Compared with plasma without venom addition, all venoms had significant anticoagulant effects, with a 160-fold range of concentrations having similar anticoagulant effects in a species-specific manner. Carbon monoxide significantly inhibited the anticoagulant effect of all venoms tested, but inhibition was not complete in all cases.

CONCLUSION

Given that individual neurotoxin activity often depends on intact activity that includes anticoagulant action, it may be possible that carbon monoxide inhibits neurotoxicity. Future investigation is justified to assess such carbon monoxide mediated inhibition with purified neurotoxins in vitro and in vivo.

Omics Technologies for Profiling Toxin Diversity and Evolution in Snake Venom: Impacts on the Discovery of Therapeutic and Diagnostic Agents

Snake venoms are primarily composed of proteins and peptides, and these toxins have developed high selectivity to their biological targets. This makes venoms interesting for exploration into protein evolution and structure-function relationships. A single venom protein superfamily can exhibit a variety of pharmacological effects; these variations in activity originate from differences in functional sites, domains, posttranslational modifications, and the formations of toxin complexes. In this review, we discuss examples of how the major venom protein superfamilies have diversified, as well as how newer technologies in the omics fields, such as genomics, transcriptomics, and proteomics, can be used to characterize both known and unknown toxins.Because toxins are bioactive molecules with a rich diversity of activities, they can be useful as therapeutic and diagnostic agents, and successful examples of toxin applications in these areas are also reviewed. With the current rapid pace of technology, snake venom research and its applications will only continue to expand.

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.

Proteomic insights into short neurotoxin-driven, highly neurotoxic venom of Philippine cobra (Naja philippinensis) and toxicity correlation of cobra envenomation in Asia

The Philippine cobra, Naja philippinensis, is a WHO Category 1 venomous snake of medical importance responsible for fatal envenomation in the northern Philippines. To elucidate the venom proteome and pathophysiology of envenomation, N. philippinensis venom proteins were decomplexed with reverse-phase high-performance liquid chromatography, and protein fractions were subsequently digested with trypsin, followed by nano-liquid chromatography-tandem mass spectrometry analysis and data mining. Three-finger toxins (3FTX, 66.64% of total venom proteins) and phospholipases A₂ (PLA₂, 22.88%) constitute the main bulk of venom proteome. Other proteins are present at low abundances (<4% each); these include metalloproteinase, serine protease, cobra venom factor, cysteine-rich secretory protein, vespryn, phosphodiesterase, 5' nucleotidase and nerve growth factor. In the three-finger toxin family, the alpha-neurotoxins comprise solely short neurotoxins (SNTX, 44.55%), supporting that SNTX is the principal toxin responsible for neuromuscular paralysis and lethality reported in clinical envenomation. Cytotoxins (CTX) are the second most abundant 3FTX proteins in the venom (21.31%). The presence of CTX correlates with the venom cytotoxic effect, which is more prominent in murine cells than in human cells. From the practical standpoint, SNTX-driven neuromuscular paralysis is significant in N. philippinensis envenomation. Antivenom production and treatment should be tailored accordingly to ensure effective neutralization of SNTX. BIOLOGICAL SIGNIFICANCE: The venom proteome of Naja philippinensis, the Philippine cobra, is unravelled for the first time. Approximately half the protein bulk of the venom is made up of short neurotoxins (44.55% of the total venom proteins). As the only alpha-neurotoxins present in the venom, short neurotoxins are the causative toxins of the post-synaptic blockade and fast-onset neuromuscular paralysis in N. philippinensis envenomation. A substantial amount of cytotoxins (21.31%) was also detected in N. philippinensis venom, supporting that the venom can be cytotoxic although the effect is much weaker in human cells compared to murine cells. The finding is consistent with the low incidence of local tissue necrosis in N. philippinensis envenomation, although this does not negate the need for monitoring and care of bite wound in the patients.

Distinctive Distribution of Secretory Phospholipases A₂ in the Venoms of Afro-Asian Cobras (Subgenus: Naja, Afronaja, Boulengerina and Uraeus)

The protein abundances of phospholipases A₂ in cobra venom proteomes appear to vary among cobra species. To determine the unique distribution of snake venom phospholipases A₂ (svPLA₂) in the cobras, the svPLA₂ activities for 15 cobra species were examined with an acidimetric and a colorimetric assay, using egg yolk suspension and 4-nitro-3-octanoyloxy benzoic acid (NOBA) as the substrate. The colorimetric assay showed significant correlation between svPLA₂ enzymatic activities with the svPLA₂ protein abundances in venoms. High svPLA₂ activities were observed in the venoms of Asiatic spitting cobras (Naja sputatrix, Naja sumatrana) and moderate activities in Asiatic non-spitters (Naja naja, Naja atra, Naja kaouthia), African spitters (subgenus Afronaja), and forest cobra (subgenus Boulengerina). African non-spitting cobras of subgenus Uraeus (Naja haje, Naja annulifera, Naja nivea, Naja senegalensis) showed exceptionally low svPLA₂ enzymatic activities. The negligible PLA₂ activity in Uraeus cobra venoms implies that PLA₂ may not be ubiquitous in all snake venoms. The svPLA₂ in cobra envenoming varies depending on the cobra species. This may potentially influence the efficacy of cobra antivenom in specific use for venom neutralization.

Exploring the Diversity and Novelty of Toxin Genes in Naja sumatrana, the Equatorial Spitting Cobra from Malaysia through De Novo Venom-Gland Transcriptomics

The equatorial spitting cobra, Naja sumatrana, is a distinct species of medically important venomous snakes, listed as WHO Category 1 in Southeast Asia. The diversity of its venom genes has not been comprehensively examined, although a few toxin sequences annotated to Naja sputatrix were reported previously through cloning studies. To investigate this species venom genes’ diversity, de novo venom-gland transcriptomics of N. sumatrana from West Malaysia was conducted using next-generation sequencing technology. Genes encoding toxins represented only 60 of the 55,396 transcripts, but were highly expressed, contributing to 79.22% of total gene expression (by total FPKM) in the venom-glands. The toxin transcripts belong to 21 families, and 29 transcripts were further identified as full-length. Three-finger toxins (3FTx) composed of long, short, and non-conventional groups, constituted the majority of toxin transcripts (91.11% of total toxin FPKM), followed by phospholipase A₂ (PLA₂, 7.42%)—which are putatively pro-inflammatory and cytotoxic. The remaining transcripts in the 19 families were expressed at extremely low levels. Presumably, these toxins were associated with ancillary functions. Our findings unveil the diverse toxin genes unique to N. sumatrana, and provide insights into the pathophysiology of N. sumatrana envenoming.

Malaysian Cobra Venom: A Potential Source of Anti-Cancer Therapeutic Agents

Cancer is a deadly disease and there is an urgent need for the development of effective and safe therapeutic agents to treat it. Snake venom is a complex mixture of bioactive proteins that represents an attractive source of novel and naturally-derived anticancer agents. Malaysia is one of the world’s most biodiverse countries and is home to various venomous snake species, including cobras. Naja kaouthia, Naja sumatrana, and Ophiophagus hannah are three of the most common cobra species in Malaysia and are of medical importance. Over the past decades, snake venom has been identified as a potential source of therapeutic agents, including anti-cancer agents. This present review highlights the potential anticancer activity of the venom and purified venom protein of N. kaouthia, N. sumatrana, and O. hannah. In conclusion, this review highlights the important role of the venom from Malaysian cobras as an important resource that researchers can exploit to further investigate its potential in cancer treatment.

Antibody Cross-Reactivity in Antivenom Research

Antivenom cross-reactivity has been investigated for decades to determine which antivenoms can be used to treat snakebite envenomings from different snake species. Traditionally, the methods used for analyzing cross-reactivity have been immunodiffusion, immunoblotting, enzyme-linked immunosorbent assay (ELISA), enzymatic assays, and in vivo neutralization studies. In recent years, new methods for determination of cross-reactivity have emerged, including surface plasmon resonance, antivenomics, and high-density peptide microarray technology. Antivenomics involves a top-down assessment of the toxin-binding capacities of antivenoms, whereas high-density peptide microarray technology may be harnessed to provide in-depth knowledge on which toxin epitopes are recognized by antivenoms. This review provides an overview of both the classical and new methods used to investigate antivenom cross-reactivity, the advantages and disadvantages of each method, and examples of studies using the methods. A special focus is given to antivenomics and high-density peptide microarray technology as these high-throughput methods have recently been introduced in this field and may enable more detailed assessments of antivenom cross-reactivity.

An overview of Phoneutria nigriventer spider venom using combined transcriptomic and proteomic approaches

Phoneutria nigriventer is one of the largest existing true spiders and one of the few considered medically relevant. Its venom contains several neurotoxic peptides that act on different ion channels and chemical receptors of vertebrates and invertebrates. Some of these venom toxins have been shown as promising models for pharmaceutical or biotechnological use. However, the large diversity and the predominance of low molecular weight toxins in this venom have hampered the identification and deep investigation of the less abundant toxins and the proteins with high molecular weight. Here, we combined conventional and next-generation cDNA sequencing with Multidimensional Protein Identification Technology (MudPIT), to obtain an in-depth panorama of the composition of P. nigriventer spider venom. The results from these three approaches showed that cysteine-rich peptide toxins are the most abundant components in this venom and most of them contain the Inhibitor Cysteine Knot (ICK) structural motif. Ninety-eight sequences corresponding to cysteine-rich peptide toxins were identified by the three methodologies and many of them were considered as putative novel toxins, due to the low similarity to previously described toxins. Furthermore, using next-generation sequencing we identified families of several other classes of toxins, including CAPs (Cysteine Rich Secretory Protein-CRiSP, antigen 5 and Pathogenesis-Related 1-PR-1), serine proteinases, TCTPs (translationally controlled tumor proteins), proteinase inhibitors, metalloproteinases and hyaluronidases, which have been poorly described for this venom. This study provides an overview of the molecular diversity of P. nigriventer venom, revealing several novel components and providing a better basis to understand its toxicity and pharmacological activities.

Studying Smaller and Neglected Organisms in Modern Evolutionary Venomics Implementing RNASeq (Transcriptomics)-A Critical Guide

Venoms are evolutionary key adaptations that species employ for defense, predation or competition. However, the processes and forces that drive the evolution of venoms and their toxin components remain in many aspects understudied. In particular, the venoms of many smaller, neglected (mostly invertebrate) organisms are not characterized in detail, especially with modern methods. For the majority of these taxa, even their biology is only vaguely known. Modern evolutionary venomics addresses the question of how venoms evolve by applying a plethora of -omics methods. These recently became so sensitive and enhanced that smaller, neglected organisms are now more easily accessible to comparatively study their venoms. More knowledge about these taxa is essential to better understand venom evolution in general. The methodological core pillars of integrative evolutionary venomics are genomics, transcriptomics and proteomics, which are complemented by functional morphology and the field of protein synthesis and activity tests. This manuscript focuses on transcriptomics (or RNASeq) as one toolbox to describe venom evolution in smaller, neglected taxa. It provides a hands-on guide that discusses a generalized RNASeq workflow, which can be adapted, accordingly, to respective projects. For neglected and small taxa, generalized recommendations are difficult to give and conclusions need to be made individually from case to case. In the context of evolutionary venomics, this overview highlights critical points, but also promises of RNASeq analyses. Methodologically, these concern the impact of read processing, possible improvements by perfoming multiple and merged assemblies, and adequate quantification of expressed transcripts. Readers are guided to reappraise their hypotheses on venom evolution in smaller organisms and how robustly these are testable with the current transcriptomics toolbox. The complementary approach that combines particular proteomics but also genomics with transcriptomics is discussed as well. As recently shown, comparative proteomics is, for example, most important in preventing false positive identifications of possible toxin transcripts. Finally, future directions in transcriptomics, such as applying 3rd generation sequencing strategies to overcome difficulties by short read assemblies, are briefly addressed.

Evaluating the Performance of De Novo Assembly Methods for Venom-Gland Transcriptomics

Venom-gland transcriptomics is a key tool in the study of the evolution, ecology, function, and pharmacology of animal venoms. In particular, gene-expression variation and coding sequences gained through transcriptomics provide key information for explaining functional venom variation over both ecological and evolutionary timescales. The accuracy and usefulness of inferences made through transcriptomics, however, is limited by the accuracy of the transcriptome assembly, which is a bioinformatic problem with several possible solutions. Several methods have been employed to assemble venom-gland transcriptomes, with the Trinity assembler being the most commonly applied among them. Although previous evidence of variation in performance among assembly software exists, particularly regarding recovery of difficult-to-assemble multigene families such as snake venom metalloproteinases, much work to date still employs a single assembly method. We evaluated the performance of several commonly used de novo assembly methods for the recovery of both nontoxin transcripts and complete, high-quality venom-gene transcripts across eleven snake and four scorpion transcriptomes. We varied k-mer sizes used by some assemblers to evaluate the impact of k-mer length on transcript recovery. We showed that the recovery of nontoxin transcripts and toxin transcripts is best accomplished through different assembly software, with SDT at smaller k-mer lengths and Trinity being best for nontoxin recovery and a combination of SeqMan NGen and a seed-and-extend approach implemented in Extender as the best means of recovering a complete set of toxin transcripts. In particular, Extender was the only means tested capable of assembling multiple isoforms of the diverse snake venom metalloproteinase family, while traditional approaches such as Trinity recovered at most one metalloproteinase transcript. Our work demonstrated that traditional metrics of assembly performance are not predictive of performance in the recovery of complete and high quality toxin genes. Instead, effective venom-gland transcriptomic studies should combine and quality-filter the results of several assemblers with varying algorithmic strategies.

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.

Computational Studies of Snake Venom Toxins

Most snake venom toxins are proteins, and participate to envenomation through a diverse array of bioactivities, such as bleeding, inflammation, and pain, cytotoxic, cardiotoxic or neurotoxic effects. The venom of a single snake species contains hundreds of toxins, and the venoms of the 725 species of venomous snakes represent a large pool of potentially bioactive proteins. Despite considerable discovery efforts, most of the snake venom toxins are still uncharacterized. Modern bioinformatics tools have been recently developed to mine snake venoms, helping focus experimental research on the most potentially interesting toxins. Some computational techniques predict toxin molecular targets, and the binding mode to these targets. This review gives an overview of current knowledge on the ~2200 sequences, and more than 400 three-dimensional structures of snake toxins deposited in public repositories, as well as of molecular modeling studies of the interaction between these toxins and their molecular targets. We also describe how modern bioinformatics have been used to study the snake venom protein phospholipase A2, the small basic myotoxin Crotamine, and the three-finger peptide Mambalgin.

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.

Combined venom profiling and cytotoxicity screening of the Radde's mountain viper (Montivipera raddei) and Mount Bulgar Viper (Montivipera bulgardaghica) with potent cytotoxicity against human A549 lung carcinoma cells

Here we report the first characterization of the endemic Mount Bulgar Viper (Montivipera bulgardaghica) and Radde's mountain viper (Montivipera raddei) venom by a combined approach using intact mass profiling and bottom-up proteomics. The cytotoxicity screening of crude venom as well as isolated serine proteases revealed a high activity against A549 human lung carcinoma cells. By means of intact mass profiling of native and reduced venom we observed basic and acidic phospholipases type A₂. Moreover, the analysis revealed snake venom metalloproteases, cysteine-rich secretory proteins, disintegrins, snake venom serine proteases, C-type lectins, a vascular endothelial growth factor and an L-amino acid oxidase.