Venom Proteome of Spine-Bellied Sea Snake (Hydrophis curtus) from Penang, Malaysia: Toxicity Correlation, Immunoprofiling and Cross-Neutralization by Sea Snake Antivenom

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 and Immunoprofiling of Persian Horned Viper Venom, Pseudocerastes Persicus, from Central Part of Iran

Numerous species of venomous snakes of medical importance exist in Iran. Pseudocerastes persicus (P. persicus), one of the medically important snakes, also called the Persian horned viper, has a geographical spread that extends to the east, southwest, and central areas of Iran and is endemic across the wider region. As a result, this species is responsible for many snakebite occurrences. Venom from P. persicus found in the central province of Semnan contains phospholipase A2 and L-amino acid oxidase activities, and high toxic potency. The venom was fractionated by reverse-phase high-performance liquid chromatography (HPLC) and analyzed by Sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE), Western blotting and two-dimensional electrophoresis. Using liquid chromatography with tandem mass spectrometry (LC-MS/MS), a range of components were identified, consistent with the biochemical and toxicological properties of the venom. Proteins identified from 2D electrophoresis and shotgun methods included metallo- and serine proteases, phospholipases, oxidases, and Kunitz trypsin inhibitors, along with many other components at lower qualitative abundance. This study provides a more detailed understanding of the protein profile of Iranian P. persicus venom, which can be effective in the production of an effective antidote against it. The analysis of the resulting data shows that there is a wide range of proteins in the venom of the Persian horned viper. This information can provide a better understanding of how venom is neutralized by polyclonal antivenom. Considering the wide presence of this snake and its related species in Iran and surrounding countries, knowing the venom protein profile of this family can be of great support to antivenom producers such as Razi Vaccine & Serum Research Institute in the preparation of regional antivenoms.

Medically important snakes and snakebite envenoming in Iran

Snakebite is a common health condition in Iran with a diverse snake fauna, especially in tropical southern and mountainous western areas of the country with plethora of snake species. The list of medically important snakes, circumstances and effects of their bite, and necessary medical care require critical appraisal and should be updated regularly. This study aims to review and map the distributions of medically important snake species of Iran, re-evaluate their taxonomy, review their venomics, describe the clinical effects of envenoming, and discuss medical management and treatment, including the use of antivenom. Nearly 350 published articles and 26 textbooks with information on venomous and mildly venomous snake species and snakebites of Iran, were reviewed, many in Persian (Farsi) language, making them relatively inaccessible to an international readership. This has resulted in a revised updated list of Iran's medically important snake species, with taxonomic revisions of some, compilation of their morphological features, remapping of their geographical distributions, and description of species-specific clinical effects of envenoming. Moreover, the antivenom manufactured in Iran is discussed, together with treatment protocols that have been developed for the hospital management of envenomed patients.

De Novo Venom Gland Transcriptome Assembly and Characterization for Calloselasma rhodostoma (Kuhl, 1824), the Malayan Pit Viper from Malaysia: Unravelling Toxin Gene Diversity in a Medically Important Basal Crotaline

In Southeast Asia, the Malayan Pit Viper (Calloselasma rhodostoma) is a venomous snake species of medical importance and bioprospecting potential. To unveil the diversity of its toxin genes, this study de novo assembled and analyzed the venom gland transcriptome of C. rhodostoma from Malaysia. The expression of toxin genes dominates the gland transcriptome by 53.78% of total transcript abundance (based on overall FPKM, Fragments Per Kilobase Million), in which 92 non-redundant transcripts belonging to 16 toxin families were identified. Snake venom metalloproteinase (SVMP, PI > PII > PIII) is the most dominant family (37.84% of all toxin FPKM), followed by phospholipase A₂ (29.02%), bradykinin/angiotensin-converting enzyme inhibitor-C-type natriuretic peptide (16.30%), C-type lectin (CTL, 10.01%), snake venom serine protease (SVSP, 2.81%), L-amino acid oxidase (2.25%), and others (1.78%). The expressions of SVMP, CTL, and SVSP correlate with hemorrhagic, anti-platelet, and coagulopathic effects in envenoming. The SVMP metalloproteinase domains encode hemorrhagins (kistomin and rhodostoxin), while disintegrin (rhodostomin from P-II) acts by inhibiting platelet aggregation. CTL gene homologues uncovered include rhodocytin (platelet aggregators) and rhodocetin (platelet inhibitors), which contribute to thrombocytopenia and platelet dysfunction. The major SVSP is a thrombin-like enzyme (an ancrod homolog) responsible for defibrination in consumptive coagulopathy. The findings provide insight into the venom complexity of C. rhodostoma and the pathophysiology of envenoming.

A current perspective on snake venom composition and constituent protein families

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

An immunoinformatic approach to assessing the immunogenic capacity of alpha-neurotoxins in elapid snake venoms

Introduction: Most elapid snakes produce venoms that contain alpha-neurotoxins (α-NTXs), which are proteins that cause post-synaptic blockade and paralysis in snakebite envenoming. However, existing elapid antivenoms are known for their low potency in neutralizing the neurotoxic activity of α-NTXs, while the immunological basis has not been elucidated. Methods: In this study, a structure-based major histocompatibility complex II (MHCII) epitope predictor of horse (Equus caballus), complemented with DM-editing determinant screening algorithm was adopted to assess the immunogenicity of α-NTXs in the venoms of major Asiatic elapids (Naja kaouthia, Ophiophagus hannah, Laticauda colubrina, Hydrophis schistosus, Hydrophis curtus). Results: The scoring metric M₂R, representing the relative immunogenic performance of respective α-NTXs, showed all α-NTXs have an overall low M₂R of <0.3, and most of the predicted binders feature non-optimal P1 anchor residues. The M₂R scores correlate strongly (R² = 0.82) with the potency scores (p-score) generated based on the relative abundances of α-NTXs and the neutralization potency of commercial antivenoms. Discussion: The immunoinformatic analysis indicates that the inferior antigenicity of α-NTXs is not only due to their small molecular size but also the subpar immunogenicity affected by their amino acid composition. Structural modification with conjugation and synthetic epitope as immunogen may potentially enhance the immunogenicity for improved antivenom potency against α-NTXs of elapid snakes.

In Vitro neurotoxicity and myotoxicity of Malaysian Naja sumatrana and Naja kaouthia venoms: Neutralization by monovalent and Neuro Polyvalent Antivenoms from Thailand

Naja sumatrana and Naja kaouthia are medically important elapids species found in Southeast Asia. Snake bite envenoming caused by these species may lead to morbidity or mortality if not treated with the appropriate antivenom. In this study, the in vitro neurotoxic and myotoxic effects N. sumatrana and N. kaouthia venoms from Malaysian specimens were assessed and compared. In addition, the neutralizing capability of Cobra Antivenom (CAV), King Cobra Antivenom (KCAV) and Neuro Polyvalent Antivenom (NPAV) from Thailand were compared. Both venoms produced concentration-dependent neurotoxic and myotoxic effects in the chick biventer cervicis nerve-muscle preparation. Based on the time to cause 90% inhibition of twitches (i.e. t₉₀) N. kaouthia venom displayed more potent neurotoxic and myotoxic effects than N. sumatrana venom. All three of the antivenoms significantly attenuated venom-induced twitch reduction of indirectly stimulated tissues when added prior to venom. When added after N. sumatrana venom, at the t₉₀ time point, CAV and NPAV partially restored the twitch height but has no significant effect on the reduction in twitch height caused by N. kaouthia venom. The addition of KCAV, at the t₉₀ time point, did not reverse the attenuation of indirectly stimulated twitches caused by either venom. In addition, none of the antivenoms, when added prior to venom, prevented attenuation of directly stimulated twitches. Differences in the capability of antivenoms, especially NPAV and CAV, to reverse neurotoxicity and myotoxicity indicate that there is a need to isolate and characterize neurotoxins and myotoxins from Malaysian N. kaouthia and N. sumatrana venoms to improve neutralization capability of the antivenoms.

Equatorial Spitting Cobra (Naja sumatrana) from Malaysia (Negeri Sembilan and Penang), Southern Thailand, and Sumatra: Comparative Venom Proteomics, Immunoreactivity and Cross-Neutralization by Antivenom

The Equatorial Spitting Cobra (Naja sumatrana) is a medically important venomous snake species in Southeast Asia. Its wide geographical distribution implies potential intra-specific venom variation, while there is no species-specific antivenom available to treat its envenoming. Applying a protein-decomplexing proteomic approach, the study showed that three-finger toxins (3FTX), followed by phospholipases A₂ (PLA₂), were the major proteins well-conserved across N. sumatrana venoms of different locales. Variations were noted in the subtypes and relative abundances of venom proteins. Of note, alpha-neurotoxins (belonging to 3FTX) are the least in the Penang specimen (Ns-PG, 5.41% of total venom proteins), compared with geographical specimens from Negeri Sembilan (Ns-NS, 14.84%), southern Thailand (Ns-TH, 16.05%) and Sumatra (Ns-SU, 10.81%). The alpha-neurotoxin abundance, in general, correlates with the venom’s lethal potency. The Thai Naja kaouthia Monovalent Antivenom (NkMAV) was found to be immunoreactive toward the N. sumatrana venoms and is capable of cross-neutralizing N. sumatrana venom lethality to varying degrees (potency = 0.49–0.92 mg/mL, interpreted as the amount of venom completely neutralized per milliliter of antivenom). The potency was lowest against NS-SU venom, implying variable antigenicity of its lethal alpha-neurotoxins. Together, the findings suggest the para-specific and geographical utility of NkMAV as treatment for N. sumatrana envenoming in Southeast Asia.

An immunodetection assay developed using cobra cytotoxin-specific antibodies: Potential diagnostics for cobra envenoming

Cobra (Naja spp.) envenoming is a life-threatening medical emergency, and a correct diagnosis is crucial to initiating timely and appropriate antivenom treatment. However, snakebite diagnostics remain unavailable in Southeast Asia. This study, therefore, developed an immunodetection assay with a potential diagnostic application for cobra envenoming. The cytotoxin of Naja kaouthia (Thai Monocled Cobra) (Nk-CTX) was purified from its venom to produce CTX-specific antibodies in rabbits and chickens. A double-antibody sandwich enzyme-linked immunosorbent assay was developed using the purified anti-Nk-CTX antibodies (immunoglobulin G and immunoglobulin Y), and its selectivity, specificity, and sensitivity for the venoms of five major cobra species in Southeast Asia (N. kaouthia, Naja sumatrana, Naja sputatrix, and Naja siamensis, Naja philippinensis) were studied. The results showed the immunoassay discriminates cobra venoms from other species commonly implicated in snakebites in Southeast Asia, i.e., the Malayan Krait, Many-banded Krait, King Cobra, Eastern Russell's Viper, Malayan Pit Viper and White-lipped Pit Viper. The immunoassay has a high sensitivity for the five cobra venoms, with detection limits (LoD) ranging from 0.6 to 2.6 ng/ml. Together, the findings suggest the potential diagnostic application of the cytotoxin immunoassay for cobra envenoming. The immunoassay was found to exhibit high immunoreactivity toward ten Asiatic cobra venoms (absorbance>1.5), in contrast to African cobra venoms with low immunoreactivity (absorbance<0.9). Considering the varying CTX antigenicity between Asiatic and African cobras, the immunoassay for African cobras should utilize antibodies produced specifically from the cytotoxins of African cobra venoms.

Separation and Analytical Techniques Used in Snake Venomics: A Review Article

The deleterious consequences of snake envenomation are due to the extreme protein complexity of snake venoms. Therefore, the identification of their components is crucial for understanding the clinical manifestations of envenomation pathophysiology and for the development of effective antivenoms. In addition, snake venoms are considered as libraries of bioactive molecules that can be used to develop innovative drugs. Numerous separation and analytical techniques are combined to study snake venom composition including chromatographic techniques such as size exclusion and RP-HPLC and electrophoretic techniques. Herein, we present in detail these existing techniques and their applications in snake venom research. In the first part, we discuss the different possible technical combinations that could be used to isolate and purify SV proteins using what is known as bioassay-guided fractionation. In the second part, we describe four different proteomic strategies that could be applied for venomics studies to evaluate whole venom composition, including the mostly used technique: RP-HPLC. Eventually, we show that to date, there is no standard technique used for the separation of all snake venoms. Thus, different combinations might be developed, taking into consideration the main objective of the study, the available resources, and the properties of the target molecules to be isolated.

On characterizing the Red-headed Krait (Bungarus flaviceps) venom: Decomplexation proteomics, immunoreactivity and toxicity cross-neutralization by hetero-specific antivenoms

The Red-headed Krait (Bungarus flaviceps) is a medically important venomous snake species in Southeast Asia, while there is no specific antivenom available for its envenoming. This study investigated the venom composition through a decomplexation proteomic approach, and examined the immunoreactivity as well as cross-neutralization efficacy of two hetero-specific krait antivenoms, Bungarus candidus Monovalent Antivenom (BcMAV) and Bungarus fasciatus Monovalent Antivenom (BfMAV), against the venom of B. flaviceps from Peninsular Malaysia. A total of 43 non-redundant proteoforms belonging to 10 toxin families were identified in the venom proteome, which is dominated by phospholipases A₂ including beta-bungarotoxin lethal subunit (56.20 % of total venom proteins), Kunitz-type serine protease inhibitors (19.40 %), metalloproteinases (12.85 %) and three-finger toxins (7.73 %). The proteome varied in quantitative aspect from the earlier reported Indonesian (Sumatran) sample, suggesting geographical venom variation. BcMAV and BfMAV were immunoreactive toward the B. flaviceps venom, with BcMAV being more efficacious in immunological binding. Both antivenoms cross-neutralized the venom lethality with varying efficacy, where BcMAV was more potent than BfMAV by ~13 times (normalized potency: 38.04 mg/g vs. 2.73 mg/g, defined as the venom amount completely neutralized by one-gram antivenom protein), supporting the potential utility of BcMAV for para-specific neutralization against B. flaviceps venom.

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

Preparation and Evaluation of a Horse Antiserum against the Venom of Sea Snake Hydrophis curtus from Hainan, China

Sea snake venom is extremely toxic, and it can induce severe respiratory failure and cause high mortality. The most effective first aid treatment for sea snake bites is to inject antivenom as soon as possible. However, in China, there are only four types of terrestrial snake antivenoms, none of which are effective in the treatment of sea snake bites. In order to develop an antivenom for the dominant species of sea snakes in Chinese seas, Hydrophis curtus venom (HcuV) was chosen as the antigen to immunize horses. From immune plasma, a high-titer Hydrophis curtus antivenom (HcuAV) was prepared. In vitro assessment showed that HcuAV had a cross-neutralizing capacity against HcuV and Hydrophis cyanocinctus venom (HcyV). In vivo assessment indicated that HcuAV injection could significantly improve the survival rates of the HcuV and HcyV envenomated mice (0% to 100% and 87.5%, respectively) when it was injected at a sufficient amount within the shortest possible time. In addition, HcuAV could also effectively alleviate multiple organ injuries caused by HcuV. These results provide experimental support for the future clinical application of HcuAV.

Snake Venomics: Fundamentals, Recent Updates, and a Look to the Next Decade

Venomic research, powered by techniques adapted from proteomics, transcriptomics, and genomics, seeks to unravel the diversity and complexity of venom through which knowledge can be applied in the treatment of envenoming, biodiscovery, and conservation. Snake venom proteomics is most extensively studied, but the methods varied widely, creating a massive amount of information which complicates data comparison and interpretation. Advancement in mass spectrometry technology, accompanied by growing databases and sophisticated bioinformatic tools, has overcome earlier limitations of protein identification. The progress, however, remains challenged by limited accessibility to samples, non-standardized quantitative methods, and biased interpretation of -omic data. Next-generation sequencing (NGS) technologies enable high-throughput venom-gland transcriptomics and genomics, complementing venom proteomics by providing deeper insights into the structural diversity, differential expression, regulation and functional interaction of the toxin genes. Venomic tissue sampling is, however, difficult due to strict regulations on wildlife use and transfer of biological materials in some countries. Limited resources for techniques and funding are among other pertinent issues that impede the progress of venomics, particularly in less developed regions and for neglected species. Genuine collaboration between international researchers, due recognition of regional experts by global organizations (e.g., WHO), and improved distribution of research support, should be embraced.

Marine Origin Ligands of Nicotinic Receptors: Low Molecular Compounds, Peptides and Proteins for Fundamental Research and Practical Applications

The purpose of our review is to briefly show what different compounds of marine origin, from low molecular weight ones to peptides and proteins, offer for understanding the structure and mechanism of action of nicotinic acetylcholine receptors (nAChRs) and for finding novel drugs to combat the diseases where nAChRs may be involved. The importance of the mentioned classes of ligands has changed with time; a protein from the marine snake venom was the first excellent tool to characterize the muscle-type nAChRs from the electric ray, while at present, muscle and α7 receptors are labeled with the radioactive or fluorescent derivatives prepared from α-bungarotoxin isolated from the many-banded krait. The most sophisticated instruments to distinguish muscle from neuronal nAChRs, and especially distinct subtypes within the latter, are α-conotoxins. Such information is crucial for fundamental studies on the nAChR revealing the properties of their orthosteric and allosteric binding sites and mechanisms of the channel opening and closure. Similar data are provided by low-molecular weight compounds of marine origin, but here the main purpose is drug design. In our review we tried to show what has been obtained in the last decade when the listed classes of compounds were used in the nAChR research, applying computer modeling, synthetic analogues and receptor mutants, X-ray and electron-microscopy analyses of complexes with the nAChRs, and their models which are acetylcholine-binding proteins and heterologously-expressed ligand-binding domains.

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.

De novo venom gland transcriptomics of Calliophis bivirgata flaviceps: uncovering the complexity of toxins from the Malayan blue coral snake

Background:

The Malayan blue coral snake, Calliophis bivirgata flaviceps, is a medically important venomous snake in Southeast Asia. However, the complexity and diversity of its venom genes remain little explored.

Methods:

To address this, we applied high-throughput next-generation sequencing to profile the venom gland cDNA libraries of C. bivirgata flaviceps. The transcriptome was de novo assembled, followed by gene annotation, multiple sequence alignment and analyses of the transcripts.

Results:

A total of 74 non-redundant toxin-encoding genes from 16 protein families were identified, with 31 full-length toxin transcripts. Three-finger toxins (3FTx), primarily delta-neurotoxins and cardiotoxin-like/cytotoxin-like proteins, were the most diverse and abundantly expressed. The major 3FTx (Cb_FTX01 and Cb_FTX02) are highly similar to calliotoxin, a delta-neurotoxin previously reported in the venom of C. bivirgata. This study also revealed a conserved tyrosine residue at position 4 of the cardiotoxin-like/cytotoxin-like protein genes in the species. These variants, proposed as Y-type CTX-like proteins, are similar to the H-type CTX from cobras. The substitution is conservative though, preserving a less toxic form of elapid CTX-like protein, as indicated by the lack of venom cytotoxicity in previous laboratory and clinical findings. The ecological role of these toxins, however, remains unclear. The study also uncovered unique transcripts that belong to phospholipase A₂ of Groups IA and IB, and snake venom metalloproteinases of PIII subclass, which show sequence variations from those of Asiatic elapids.

Conclusion:

The venom gland transcriptome of C. bivirgata flaviceps from Malaysia was de novo assembled and annotated. The diversity and expression profile of toxin genes provide insights into the biological and medical importance of the species.

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

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.

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

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

Immunoreactivity and neutralization capacity of Philippine cobra antivenom against Naja philippinensis and Naja samarensis venoms

BACKGROUND

The Philippine cobra (Naja philippinensis) and Samar cobra (Naja samarensis) are two WHO Category 1 medically important venomous snakes in the Philippines. Philippine cobra antivenom (PCAV) is the only antivenom available in the country, but its neutralization capacity against the venoms of N. philippinensis and hetero-specific N. samarensis has not been reported. This knowledge gap greatly hinders the optimization of antivenom use in the region.

METHODS

This study examined the immunological binding and neutralization capacity of PCAV against the two cobra venoms using WHO-recommended protocols.

RESULTS

In mice, both venoms were highly neurotoxic and lethal with a median lethal dose of 0.18 and 0.20 µg/g, respectively. PCAV exhibited strong and comparable immunoreactivity toward the venoms, indicating conserved venom antigenicity between the two allopatric species. In in vivo assay, PCAV was only moderately effective in neutralizing the toxicity of both venoms. Its potency was even lower against the hetero-specific N. samarensis venom by approximately two-fold compared with its potency against N. philippinensis venom.

CONCLUSION

The results indicated that PCAV could be used to treat N. samarensis envenomation but at a higher dose, which might increase the risk of hypersensitivity and worsen the shortage of antivenom supply in the field. Antivenom manufacturing should be improved by developing a low-dose, high-efficacy product against cobra envenomation.

De Novo Venom-Gland Transcriptomics of Spine-Bellied Sea Snake (Hydrophis curtus) from Penang, Malaysia-Next-Generation Sequencing, Functional Annotation and Toxinological Correlation

Envenomation resulted from sea snake bite is a highly lethal health hazard in Southeast Asia. Although commonly caused by sea snakes of Hydrophiinae, each species is evolutionarily distinct and thus, unveiling the toxin gene diversity within individual species is important. Applying next-generation sequencing, this study investigated the venom-gland transcriptome of Hydrophis curtus (spine-bellied sea snake) from Penang, West Malaysia. The transcriptome was de novo assembled, followed by gene annotation and sequence analyses. Transcripts with toxin annotation were only 96 in number but highly expressed, constituting 48.18% of total FPKM in the overall transcriptome. Of the 21 toxin families, three-finger toxins (3FTX) were the most abundantly expressed and functionally diverse, followed by phospholipases A₂. Lh_FTX001 (short neurotoxin) and Lh_FTX013 (long neurotoxin) were the most dominant 3FTXs expressed, consistent with the pathophysiology of envenomation. Lh_FTX001 and Lh_FTX013 were variable in amino acid compositions and predicted epitopes, while Lh_FTX001 showed high sequence similarity with the short neurotoxin from Hydrophis schistosus, supporting cross-neutralization effect of Sea Snake Antivenom. Other toxins of low gene expression, for example, snake venom metalloproteinases and L-amino acid oxidases not commonly studied in sea snake venom were also identified, enriching the knowledgebase of sea snake toxins for future study.

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.

A Neurotoxic Snake Venom without Phospholipase A2: Proteomics and Cross-Neutralization of the Venom from Senegalese Cobra, Naja senegalensis (Subgenus: Uraeus)

The Senegalese cobra, Naja senegalensis, is a non-spitting cobra species newly erected from the Naja haje complex. Naja senegalensis causes neurotoxic envenomation in Western Africa but its venom properties remain underexplored. Applying a protein decomplexation proteomic approach, this study unveiled the unique complexity of the venom composition. Three-finger toxins constituted the major component, accounting for 75.91% of total venom proteins. Of these, cardiotoxin/cytotoxin (~53%) and alpha-neurotoxins (~23%) predominated in the venom proteome. Phospholipase A₂, however, was not present in the venom, suggesting a unique snake venom phenotype found in this species. The venom, despite the absence of PLA₂, is highly lethal with an intravenous LD₅₀ of 0.39 µg/g in mice, consistent with the high abundance of alpha-neurotoxins (predominating long neurotoxins) in the venom. The hetero-specific VINS African Polyvalent Antivenom (VAPAV) was immunoreactive to the venom, implying conserved protein antigenicity in the venoms of N. senegalensis and N. haje. Furthermore, VAPAV was able to cross-neutralize the lethal effect of N. senegalensis venom but the potency was limited (0.59 mg venom completely neutralized per mL antivenom, or ~82 LD₅₀ per ml of antivenom). The efficacy of antivenom should be further improved to optimize the treatment of cobra bite envenomation in Africa.

Highlights of animal venom research on the geographical variations of toxin components, toxicities and envenomation therapy

Geographical variation of animal venom is common among venomous animals. This kind of intraspecific variation based on geographical location mainly concerned from envenomation cases and brought new problems in animal venom studies, including venom components regulatory mechanisms, differentiation of venom activities, and clinical treatment methods. At present, food is considered as the most related factor influencing venom development. Related research defined the variational venomous animal species by the comparison of venom components and activities in snakes, jellyfish, scorpions, cone snails, ants, parasitoid wasps, spiders and toads. In snake venom studies, researchers found that antivenom effectiveness was variated to different located venom samples. As described in some snake venom research, developing region-specific antivenom is the development trend. The difficulties of developing region-specific antivenom and theoretical solutions have been discussed. This review summarized biological studies of animal venom geographical variation by species, compared venom components and major biological activities of the vary venom from the same species, and listed the basic methods in comparing venom protein compositions and major toxicity differences to provide a comprehensive reference.

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.

Snake Venom Metalloproteinases (SVMPs): A structure-function update

Snake venom metalloproteinases (SVMPs) represent a diverse group of multi-domain proteins with several biological activities such as the ability to induce hemorrhage, proteolytic degradation of fibrinogen and fibrin, induction of apoptosis and inhibition of platelet aggregation. Due to these activities, SVMPs are responsible for many of the well-known pathological phenotypes in snake envenomations caused particularly by species from the Viperidae family and the Crotalinae subfamily. These proteins have been classified based on their size and domain structure into P–I, P-II and P-III classes. Comparatively, members of the P–I SVMPs possess the simplest structures, formed by the catalytic metalloproteinase domain only; the P-II SVMPs are moderately more complex, having the canonical disintegrin domain in addition to the metalloproteinase domain; members of the P-III class are more structurally varied, comprising the metalloproteinase, disintegrin-like, and cysteine-rich domains. Proteolytic cleavage, repeated domain loss and presence of other ancillary domains are responsible for structural diversities in the P-III class. However, studies continue to unveil the relationship between the structure and function of these proteins. In this review, we recovered evidences from literature on the structural peculiarities and functional classification of Snake Venom Metalloproteinases. In addition, we reflect on diversities that exist among each class while taking into account specific and up-to-date class-based activities.

A pan-specific antiserum produced by a novel immunization strategy shows a high spectrum of neutralization against neurotoxic snake venoms

Snakebite envenomation is a neglected tropical disease of high mortality and morbidity largely due to insufficient supply of effective and affordable antivenoms. Snake antivenoms are mostly effective against the venoms used in their production. It is thus crucial that effective and affordable antivenom(s) with wide para-specificity, capable of neutralizing the venoms of a large number of snakes, be produced. Here we studied the pan-specific antiserum prepared previously by a novel immunization strategy involving the exposure of horses to a ‘diverse toxin repertoire’ consisting of 12 neurotoxic Asian snake toxin fractions/ venoms from six species. This antiserum was previously shown to exhibit wide para-specificity by neutralizing 11 homologous and 16 heterologous venoms from Asia and Africa. We now show that the antiserum can neutralize 9 out of 10 additional neurotoxic venoms. Altogether, 36 snake venoms belonging to 10 genera from 4 continents were neutralized by the antiserum. Toxin profiles previously generated using proteomic techniques of these 36 venoms identified α-neurotoxins, β-neurotoxins, and cytotoxins as predominant toxins presumably neutralized by the antiserum. The bases for the wide para-specificity of the antiserum are discussed. These findings indicate that it is feasible to generate antivenoms of wide para-specificity against elapid neurotoxic venoms from different regions in the world and raises the possibility of a universal neurotoxic antivenom. This should reduce the mortality resulting from neurotoxic snakebite envenomation.

Comparative venom proteomics of banded krait (Bungarus fasciatus) from five geographical locales: Correlation of venom lethality, immunoreactivity and antivenom neutralization

The banded krait, Bungarus fasciatus is a medically important venomous snake in Asia. The wide distribution of this species in Southeast Asia and southern China indicates potential geographical variation of the venom which may impact the clinical management of snakebite envenomation. This study investigated the intraspecific venom variation of B. fasciatus from five geographical locales through a venom decomplexing proteomic approach, followed by toxinological and immunological studies. The venom proteomes composed of a total of 9 toxin families, comprising 22 to 31 proteoforms at varying abundances. The predominant proteins were phospholipase A₂ (including beta-bungarotoxin), Kunitz-type serine protease inhibitor (KSPI) and three-finger toxins (3FTx), which are toxins that cause neurotoxicity and lethality. The venom lethality varied with geographical origins of the snake, with intravenous median lethal doses (LD₅₀) ranging from 0.45-2.55 µg/g in mice. The Thai Bungarus fasciatus monovalent antivenom (BFMAV) demonstrated a dose-dependent increasing immunological binding activity toward all venoms; however, its in vivo neutralization efficacy varied vastly with normalized potency values ranging from 3 to 28 mg/g, presumably due to the compositional differences of dominant proteins in the different venoms. The findings support that antivenom use should be optimized in different geographical areas. The development of a pan-regional antivenom may be a more sustainable solution for the treatment of snakebite envenomation.

Geographical variations in king cobra (Ophiophagus hannah) venom from Thailand, Malaysia, Indonesia and China: On venom lethality, antivenom immunoreactivity and in vivo neutralization

The wide distribution of king cobra (Ophiophagus hannah), a medically important venomous snake in Asia could be associated with geographical variation in the toxicity and antigenicity of the venom. This study investigated the lethality of king cobra venoms (KCV) from four geographical locales (Malaysia, Thailand, Indonesia, China), and the immunological binding as well as in vivo neutralization activities of three antivenom products (Thai Ophiophagus hannah monovalent antivenom, OHMAV; Indonesian Serum Anti Bisa Ular, SABU; Chinese Naja atra monovalent antivenom, NAMAV) toward the venoms. The Indonesian and Chinese KCV were more lethal (median lethal dose, LD₅₀ ~0.5 μg/g) than those from Malaysia and Thailand (LD₅₀ ~1.0 μg/g). The antivenoms, composed of F(ab)'₂, were variably immunoreactive toward the KCV from all locales, with OHMAV exhibited the highest immunological binding activity. In mice, OHMAV neutralized the neurotoxic lethality of Thai KCV most effectively (normalized potency = 118 mg venom neutralized per g antivenom) followed by Malaysian, Indonesian and Chinese KCV. In comparison, the hetero-specific SABU was remarkably less potent by at least 6 to10 folds, whereas NAMAV appeared to be non-effective. The finding supports that a specific king cobra antivenom is needed for the effective treatment of king cobra envenomation in each region.

Diet Breadth Mediates the Prey Specificity of Venom Potency in Snakes

Venoms are best known for their ability to incapacitate prey. In predatory groups, venom potency is predicted to reflect ecological and evolutionary drivers relating to diet. While venoms have been found to have preyspecific potencies, the role of diet breadth on venom potencies has yet to be tested at large macroecological scales. Here, using a comparative analysis of 100 snake species, we show that the evolution of prey-specific venom potencies is contingent on the breadth of a species' diet. We find that while snake venom is more potent when tested on species closely related to natural prey items, we only find this prey-specific pattern in species with taxonomically narrow diets. While we find that the taxonomic diversity of a snakes' diet mediates the prey specificity of its venom, the species richness of its diet was not found to affect these prey-specific potency patterns. This indicates that the physiological diversity of a species' diet is an important driver of the evolution of generalist venom potencies. These findings suggest that the venoms of species with taxonomically diverse diets may be better suited to incapacitating novel prey species and hence play an important role for species within changing environments.

A Decoy-Receptor Approach Using Nicotinic Acetylcholine Receptor Mimics Reveals Their Potential as Novel Therapeutics Against Neurotoxic Snakebite

Snakebite is a neglected tropical disease that causes 138,000 deaths each year. Neurotoxic snake venoms contain small neurotoxins, including three-finger toxins (3FTxs), which can cause rapid paralysis in snakebite victims by blocking postsynaptic transmission via nicotinic acetylcholine receptors (nAChRs). These toxins are typically weakly immunogenic and thus are often not effectively targeted by current polyclonal antivenom therapies. We investigated whether nAChR mimics, also known as acetylcholine binding proteins (AChBPs), could effectively capture 3FTxs and therefore be developed as a novel class of snake-generic therapeutics for combatting neurotoxic envenoming. First, we identified the binding specificities of 3FTx from various medically important elapid snake venoms to nAChR using two recombinant nAChR mimics: the AChBP from Lymnaea stagnalis and a humanized neuronal α7 version (α7-AChBP). We next characterized these AChBP-bound and unbound fractions using SDS-PAGE and mass spectrometry. Interestingly, both mimics effectively captured long-chain 3FTxs from multiple snake species but largely failed to capture the highly related short-chain 3FTxs, suggesting a high level of binding specificity. We next investigated whether nAChR mimics could be used as snakebite therapeutics. We showed that while α7-AChBP alone did not protect against Naja haje (Egyptian cobra) venom lethality in vivo, it significantly prolonged survival times when coadministered with a nonprotective dose of antivenom. Thus, nAChR mimics are capable of neutralizing specific venom toxins and may be useful adjunct therapeutics for improving the safety and affordability of existing snakebite treatments by reducing therapeutic doses. Our findings justify exploring the future development of AChBPs as potential snakebite treatments.

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.

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.