Comparative venom gland transcriptomics of Naja kaouthia (monocled cobra) from Malaysia and Thailand: elucidating geographical venom variation and insights into sequence novelty

De Novo Assembly of Venom Gland Transcriptome of Tropidolaemus wagleri (Temple Pit Viper, Malaysia) and Insights into the Origin of Its Major Toxin, Waglerin

The venom proteome of Temple Pit Viper (Tropidolaemus wagleri) is unique among pit vipers, characterized by a high abundance of a neurotoxic peptide, waglerin. To further explore the genetic diversity of its toxins, the present study de novo assembled the venom gland transcriptome of T. wagleri from west Malaysia. Among the 15 toxin gene families discovered, gene annotation and expression analysis reveal the dominating trend of bradykinin-potentiating peptide/angiotensin-converting enzyme inhibitor-C-type natriuretic peptide (BPP/ACEI-CNP, 76.19% of all-toxin transcription) in the transcriptome, followed by P-III snake venom metalloproteases (13.91%) and other toxins. The transcript TwBNP01 of BPP/ACEI-CNP represents a large precursor gene (209 amino acid residues) containing the coding region for waglerin (24 residues). TwBNP01 shows substantial sequence variations from the corresponding genes of its sister species, Tropidolaemus subannulatus of northern Philippines, and other viperid species which diversely code for proline-rich small peptides such as bradykinin-potentiating peptides (BPPs). The waglerin/waglerin-like peptides, BPPs and azemiopsin are proline-rich, evolving de novo from multiple highly diverged propeptide regions within the orthologous BPP/ACEI-CNP genes. Neofunctionalization of the peptides results in phylogenetic constraints consistent with a phenotypic dichotomy, where Tropidolaemus spp. and Azemiops feae convergently evolve a neurotoxic trait while vasoactive BPPs evolve only in other species.

Action mechanism of snake venom l-amino acid oxidase and its double-edged sword effect on cancer treatment: Role of pannexin 1-mediated interleukin-6 expression

Snake venom l-amino acid oxidases (svLAAOs) have been recognized as promising candidates for anticancer therapeutics. However, multiple aspects of their catalytic mechanism and the overall responses of cancer cells to these redox enzymes remain ambiguous. Here, we present an analysis of the phylogenetic relationships and active site-related residues among svLAAOs and reveal that the previously proposed critical catalytic residue His 223 is highly conserved in the viperid but not the elapid svLAAO clade. To gain further insight into the action mechanism of the elapid svLAAOs, we purify and characterize the structural, biochemical, and anticancer therapeutic potentials of the Thailand elapid snake Naja kaouthia LAAO (NK-LAAO). We find that NK-LAAO, with Ser 223, exhibits high catalytic activity toward hydrophobic l-amino acid substrates. Moreover, NK-LAAO induces substantial oxidative stress-mediated cytotoxicity with the magnitude relying on both the levels of extracellular hydrogen peroxide (H2O2) and intracellular reactive oxygen species (ROS) generated during the enzymatic redox reactions, but not being influenced by the N-linked glycans on its surface. Unexpectedly, we discover a tolerant mechanism deployed by cancer cells to dampen the anticancer activities of NK-LAAO. NK-LAAO treatment amplifies interleukin (IL)-6 expression via the pannexin 1 (Panx1)-directed intracellular calcium (iCa2+) signaling pathway to confer adaptive and aggressive phenotypes on cancer cells. Accordingly, IL-6 silencing renders cancer cells vulnerable to NK-LAAO-induced oxidative stress together with abrogating NK-LAAO-stimulated metastatic acquisition. Collectively, our study urges caution when using svLAAOs in cancer treatment and identifies the Panx1/iCa2+/IL-6 axis as a therapeutic target for improving the effectiveness of svLAAOs-based anticancer therapies.

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.

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.

Identification of waprin and its microbicidal activity: A novel protein component of honeybee (Apis mellifera) venom

Bee venom is a rich source of biologically and pharmacologically active proteins. Waprin is a protein component of venoms; however, waprin has yet to be identified in bee venom. Moreover, the biological functions of waprin in venoms remain poorly characterized. Thus, in this study, we have identified and characterized waprin: a novel protein component from the venom of honeybees (Apis mellifera). The waprin in A. mellifera venom (Amwaprin) was found to consist of an 80-amino acid mature peptide, in which the whey acidic protein domain contains four conserved disulfide bonds. We discovered the presence of the Amwaprin protein in secreted venom by using an antibody against recombinant Amwaprin produced in baculovirus-infected insect cells. Recombinant Amwaprin exhibited inhibitory activity against microbial serine proteases and elastases but not thrombin or plasmin. It recognized carbohydrates in the microbial cell wall molecules and bound to the live microbial surfaces. The binding action of Amwaprin produced its microbicidal activity by inducing structural damage to bacterial and fungal cell walls. In addition, recombinant Amwaprin is heat-stable and contains no hemolytic activity. These findings demonstrate that Amwaprin acts as a microbicidal and anti-elastolytic agent.

Composition and toxicity of venom produced by araneophagous white-tailed spiders (Lamponidae: Lampona sp.)

Prey-specialised spiders are adapted to capture specific prey items, including dangerous prey. The venoms of specialists are often prey-specific and less complex than those of generalists, but their venom composition has not been studied in detail. Here, we investigated the venom of the prey-specialised white-tailed spiders (Lamponidae: Lampona), which utilise specialised morphological and behavioural adaptations to capture spider prey. We analysed the venom composition using proteo-transcriptomics and taxon-specific toxicity using venom bioassays. Our analysis identified 208 putative toxin sequences, comprising 103 peptides < 10 kDa and 105 proteins > 10 kDa. Most peptides belonged to one of two families characterised by scaffolds containing eight or ten cysteine residues. Toxin-like proteins showed similarity to galectins, leucine-rich repeat proteins, trypsins and neprilysins. The venom of Lampona was shown to be more potent against the preferred spider prey than against alternative cricket prey. In contrast, the venom of a related generalist was similarly potent against both prey types. These data provide insights into the molecular adaptations of venoms produced by prey-specialised spiders.

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.

Cytotoxicity of Venoms and Cytotoxins from Asiatic Cobras (Naja kaouthia, Naja sumatrana, Naja atra) and Neutralization by Antivenoms from Thailand, Vietnam, and Taiwan

Envenoming by cobras (Naja spp.) often results in extensive local tissue necrosis when optimal treatment with antivenom is not available. This study investigated the cytotoxicity of venoms and purified cytotoxins from the Monocled Cobra (Naja kaouthia), Taiwan Cobra (Naja atra), and Equatorial Spitting Cobra (Naja sumatrana) in a mouse fibroblast cell line, followed by neutralization of the cytotoxicity by three regional antivenoms: the Thai Naja kaouthia monovalent antivenom (NkMAV), Vietnamese snake antivenom (SAV) and Taiwanese Neuro bivalent antivenom (NBAV). The cytotoxins of N. atra (NA-CTX) and N. sumatrana (NS-CTX) were identified as P-type cytotoxins, whereas that of N. kaouthia (NK-CTX) is S-type. All venoms and purified cytotoxins demonstrated varying concentration-dependent cytotoxicity in the following trend: highest for N. atra, followed by N. sumatrana and N. kaouthia. The antivenoms moderately neutralized the cytotoxicity of N. kaouthia venom but were weak against N. atra and N. sumatrana venom cytotoxicity. The neutralization potencies of the antivenoms against the cytotoxins were varied and generally low across NA-CTX, NS-CTX, and NK-CTX, possibly attributed to limited antigenicity of CTXs and/or different formulation of antivenom products. The study underscores the need for antivenom improvement and/or new therapies in treating local tissue toxicity caused by cobra envenomings.

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.

Local Cytotoxic Effects in Cobra Envenoming: A Pilot Study

The cobra (genus Naja (N.)) is one of the most common venomous snakes. Due to its frequency and deadly complications of muscle paralysis, local necrosis, and chronic musculoskeletal disability, it should not be ignored. The pathology of devastating tissue destruction, even though specific antivenoms exist, is not fully clear. Here, we attempted to dig in envenomed tissues to study the clinical toxicology of cobra venom. Four cases of N. atra snake envenomation, in which the subjects developed advanced tissue injury, were involved in this study. We used enzyme-ligand sandwich immunoassay (ELISA) to assay the whole venom, cytotoxin A3 and short-chain neurotoxin (sNTX) in blood, bullae, wound discharge, and debrided tissue. We found that persistently high concentrations of venom and toxins, especially cytotoxin A3, were detected in bullae, wound discharge fluid and necrotic tissue of these patients even after large doses of specific antivenom treatment, and wide excision and advanced debridement could largely remove these toxins, lessen the size of necrosis, and promote wound healing. We also found that the point-of-care apparatus, ICT-Cobra kit, might be used to promptly monitor the wound condition and as one of the indicators of surgical intervention in cases of cobra envenomation in Taiwan.

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.

Isolation and Characterization of A2-EPTX-Nsm1a, a Secretory Phospholipase A2 from Malaysian Spitting Cobra (Naja sumatrana) Venom

Phospholipase A₂ (PLA₂) toxins are one of the main toxin families found in snake venom. PLA₂ toxins are associated with various detrimental effects, including neurotoxicity, myotoxicity, hemostatic disturbances, nephrotoxicity, edema, and inflammation. Although Naja sumatrana venom contains substantial quantities of PLA₂ components_, there is limited information on the function and activities of PLA₂ toxins from the venom. In this study, a secretory PLA₂ from the venom of Malaysian N. sumatrana, subsequently named A2-EPTX-Nsm1a, was isolated, purified, and characterized. A2-EPTX-Nsm1a was purified using a mass spectrometry-guided approach and multiple chromatography steps. Based on LC-MSMS, A2-EPTX-Nsm1a was found to show high sequence similarity with PLA₂ from venoms of other Naja species. The PLA₂ activity of A2-EPTX-Nsm1 was inhibited by 4-BPB and EDTA. A2-EPTX-Nsm1a was significantly less cytotoxic in a neuroblastoma cell line (SH-SY5Y) compared to crude venom and did not show a concentration-dependent cytotoxic activity. To our knowledge, this is the first study that characterizes and investigates the cytotoxicity of an Asp49 PLA₂ isolated from Malaysian N. sumatrana venom in a human neuroblastoma cell line.

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.

Proteomics and neutralization of Bungarus multicinctus (Many-banded Krait) venom: Intra-specific comparisons between specimens from China and Taiwan

The Many-banded Krait (Bungarus multicinctus) is a medically important venomous snake in East Asia. This study investigated the venom proteomes of B. multicinctus from Guangdong, southern China (BM-China) and insular Taiwan (BM-Taiwan), and the neutralization activities of two antivenom products (produced separately in China and Taiwan) against the lethal effect of the venoms. The venom proteomes of both specimens contained similar toxin families, notwithstanding small variations in the subtypes and abundances of minor components. More than 90% of the total venom proteins belong to three-finger toxins (3FTx, including alpha-neurotoxins) and phospholipases A₂ (PLA₂, including beta-bungarotoxins), supporting their key involvement in the pathophysiology of krait envenomation which manifests as pre- and post-synaptic neurotoxicity. The venoms exhibited potent neurotoxic and lethal effects with extremely low i.v. LD₅₀ of 0.027 μg/g (Bm-China) and 0.087 μg/g (Bm-Taiwan), respectively, in mice. Bungarus multicinctus monovalent antivenom (BMMAV) produced in China and Neuro bivalent antivenom (NBAV) produced in Taiwan were immunoreactive toward both venoms and their toxin fractions. The antivenoms neutralized the venom lethality variably, with BMMAV being more efficacious than NBAV by approximately two-fold. Findings suggest that the monovalent antivenom has a higher potency presumably due to its species-specificity toward the krait venom.

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.

Genome of the pincer wasp Gonatopus flavifemur reveals unique venom evolution and a dual adaptation to parasitism and predation

Background

Hymenoptera comprise extremely diverse insect species with extensive variation in their life histories. The Dryinidae, a family of solitary wasps of Hymenoptera, have evolved innovations that allow them to hunt using venom and a pair of chelae developed from the fore legs that can grasp prey. Dryinidae larvae are also parasitoids of Auchenorrhyncha, a group including common pests such as planthoppers and leafhoppers. Both of these traits make them effective and valuable for pest control, but little is yet known about the genetic basis of its dual adaptation to parasitism and predation.

Results

We sequenced and assembled a high-quality genome of the dryinid wasp Gonatopus flavifemur, which at 636.5 Mb is larger than most hymenopterans. The expansion of transposable elements, especially DNA transposons, is a major contributor to the genome size enlargement. Our genome-wide screens reveal a number of positively selected genes and rapidly evolving proteins involved in energy production and motor activity, which may contribute to the predatory adaptation of dryinid wasp. We further show that three female-biased, reproductive-associated yellow genes, in response to the prey feeding behavior, are significantly elevated in adult females, which may facilitate the egg production. Venom is a powerful weapon for dryinid wasp during parasitism and predation. We therefore analyze the transcriptomes of venom glands and describe specific expansions in venom Idgf-like genes and neprilysin-like genes. Furthermore, we find the LWS2-opsin gene is exclusively expressed in male G. flavifemur, which may contribute to partner searching and mating.

Conclusions

Our results provide new insights into the genome evolution, predatory adaptation, venom evolution, and sex-biased genes in G. flavifemur, and present genomic resources for future in-depth comparative analyses of hymenopterans that may benefit pest control.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01081-6.

Proteo-Trancriptomic Analyses Reveal a Large Expansion of Metalloprotease-Like Proteins in Atypical Venom Vesicles of the Wasp Meteorus pulchricornis (Braconidae)

Meteorus pulchricornis (Ichneumonoidea, Braconidae) is an endoparasitoid wasp of lepidopteran caterpillars. Its parasitic success relies on vesicles (named M. pulchricornis Virus-Like Particles or MpVLPs) that are synthesized in the venom gland and injected into the parasitoid host along with the venom during oviposition. In order to define the content and understand the biogenesis of these atypical vesicles, we performed a transcriptome analysis of the venom gland and a proteomic analysis of the venom and purified MpVLPs. About half of the MpVLPs and soluble venom proteins identified were unknown and no similarity with any known viral sequence was found. However, MpVLPs contained a large number of proteins labelled as metalloproteinases while the most abundant protein family in the soluble venom was that of proteins containing the Domain of Unknown Function DUF-4803. The high number of these proteins identified suggests that a large expansion of these two protein families occurred in M. pulchricornis. Therefore, although the exact mechanism of MpVLPs formation remains to be elucidated, these vesicles appear to be “metalloproteinase bombs” that may have several physiological roles in the host including modifying the functions of its immune cells. The role of DUF4803 proteins, also present in the venom of other braconids, remains to be clarified.

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.

Proteomics, toxicity and antivenom neutralization of Sri Lankan and Indian Russell's viper (Daboia russelii) venoms

BACKGROUND

The western Russell's viper (Daboia russelii) is widely distributed in South Asia, and geographical venom variation is anticipated among distant populations. Antivenoms used for Russell's viper envenomation are, however, raised typically against snakes from Southern India. The present study investigated and compared the venom proteomes of D. russelii from Sri Lanka (DrSL) and India (DrI), the immunorecognition of Indian VINS Polyvalent Antivenom (VPAV) and its efficacy in neutralizing the venom toxicity.

METHODS

The venoms of DrSL and DrI were decomplexed with C18 high-performance liquid chromatography and SDS-polyacrylamide gel electrophoresis under reducing conditions. The proteins fractionated were identified through nano-ESI-liquid chromatography-tandem mass spectrometry (LCMS/MS). The immunological studies were conducted with enzyme-linked immunosorbent assay. The neutralization of the venom procoagulant effect was evaluated in citrated human plasma. The neutralization of the venom lethality was assessed in vivo in mice adopting the WHO protocol.

RESULTS

DrSL and DrI venom proteomes showed comparable major protein families, with phospholipases A2 (PLA2) being the most abundant (> 60% of total venom proteins) and diverse (six protein forms identified). Both venoms were highly procoagulant and lethal (intravenous median lethal dose in mice, LD50 = 0.24 and 0.32 µg/g, for DrSL and DrI, respectively), while lacking hemorrhagic and anticoagulant activities. VPAV was immunoreactive toward DrSL and DrI venoms, indicating conserved protein antigenicity in the venoms. The high molecular weight venom proteins were, however, more effectively immunorecognized than small ones. VPAV was able to neutralize the coagulopathic and lethal effects of the venoms moderately.

CONCLUSION

Considering that a large amount of venom can be injected by Russell's viper during envenomation, the potency of antivenom can be further improved for optimal neutralization and effective treatment. Region-specific venoms and key toxins may be incorporated into the immunization procedure during antivenom production.

Proteome Analysis of Toxic Fractions of Iranian Cobra (Naja naja Oxiana) Snake Venom Using Two-Dimensional Electrophoresis and Mass Spectrometry

Snake venoms are mostly composed of various proteins and peptides with toxicity and pharmacological effects depending on their geographical sources. Naja naja oxiana is one of the most medically important venomous snakes in Iran and Central Asia. The bite of this type of snake can cause severe pain and swelling, as well as neurotoxicity. Without medical treatment, symptoms quickly worsen and death can occur soon. A detailed understanding of venom components can provide new insight into the production of antivenom against toxic agents instead of crude venom. Specific antibodies against toxic fractions are of utmost importance in neutralizing crude venom. Therefore, the proteome profile of these fractions of Naja naja oxidana venom was analyzed using fractionation by gel filtration, two-dimensional electrophoresis, mass spectrometry, and data mining. Base on the results, in total, 32 spots were detected and categorized into three protein families, namely three-finger toxin (3FTx), phospholipase, and Cysteine-rich secretory proteins (CRISP). These proteins consist of more than 70% crude venom all with a molecular weight below 25 kDa. The 3FTx as a highly diverse constituent in the venom of Naja species was in large quantity in this district. Short-chain neurotoxins, including short neurotoxin, cytotoxin, and muscarinic toxin-like protein, were in abundance, respectively. In conclusion, the recognition of toxic fractions of Naja naja oxiana in this region could be of great help in the production of an effective antivenom against similar compositions. It can also help the medical care department to find out the clinical sign of cobra venom. To the best of our knowledge, this was the first study to report the proteomic of toxic fractions of Naja naja oxiana in Iran.

Identifying and revealing the geographical variation in Nemopilema nomurai venom metalloprotease and phospholipase A2 activities

Venom geographical variation is common among venomous animals. This phenomenon presents problems in the development of clinical treatments and medicines against envenomation. The venomous giant jellyfish Nemopilema nomurai, Scyphozoan, is a blooming jellyfish species in the Yellow Sea and the East China Sea that causes numerous jellyfish sting cases every year. Metalloprotease and phospholipase A₂ (PLA₂) are the main components in Nemopilema nomurai venom and may activate many toxicities, such as hemolysis, inflammation and lethality. Geographical variation in the content and activity of these enzymes may cause different symptoms and therapeutic problems. For the first time, we verified metalloprotease and PLA₂ geographical variation in Nemopilema nomurai venom by performing a comparative analysis of 31 venom samples by SDS-PAGE, analyzing protease zymography, enzymatic activity, and drawing contour maps. Band locations and intensities of SDS-PAGE and protease zymograms showed geographical differences. The enzymatic activities of both metalloprotease and PLA₂ showed a trend of geographic regularity. The distribution patterns of these activities are directly shown in contour maps. Metalloproteinase activity was lower near the coast. PLA₂-like activity was lower in the Southern Yellow Sea. We surmised that metalloproteinase and PLA₂-like activities might be related to venom ontogeny and species abundance respectively, and influenced by similar environmental factors. This study provides a theoretical basis for further ecological and medical studies of Nemopilema nomurai jellyfish venom.

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.

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.

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.

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.

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

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

Interspecific and intraspecific venom enzymatic variation among cobras (Naja sp. and Ophiophagus hannah)

The genera Ophiophagus and Naja comprise part of a clade of snakes referred to as cobras, dangerously venomous front-fanged snakes in the family Elapidae responsible for significant human mortality and morbidity throughout Asia and Africa. We evaluated venom enzyme variation for eleven cobra species and three N. kaouthia populations using SDS-PAGE venom fingerprinting and numerous enzyme assays. Acetylcholinesterase and PLA₂ activities were the most variable between species, and PLA₂ activity was significantly different between Malaysian and Thailand N. kaouthia populations. Venom metalloproteinase activity was low and significantly different among most species, but levels were identical for N. kaouthia populations; minor variation in venom L-amino acid oxidase and phosphodiesterase activities were seen between cobra species. Naja siamensis venom lacked the α-fibrinogenolytic activity common to other cobra venoms. In addition, venom from N. siamensis had no detectable metalloproteinase activity and exhibited an SDS-PAGE profile with reduced abundance of higher mass proteins. Venom profiles from spitting cobras (N. siamensis, N. pallida, and N. mossambica) exhibited similar reductions in higher mass proteins, suggesting the evolution of venoms of reduced complexity and decreased enzymatic activity among spitting cobras. Generally, the venom proteomes of cobras show highly abundant three-finger toxin diversity, followed by large quantities of PLA₂s. However, PLA₂ bands and activity were very reduced for N. haje, N. annulifera and N. nivea. Venom compositionalenzy analysis provides insight into the evolution, diversification and distribution of different venom phenotypes that complements venomic data, and this information is critical for the development of effective antivenoms and snakebite treatment.

Quantitative proteomics of Naja annulifera (sub-Saharan snouted cobra) venom and neutralization activities of two antivenoms in Africa

Envenomation by Naja annulifera (snouted cobra), a non-spitting African cobra, can result in local tissue damage and fatal paralysis but a species-specific antivenom treatment is currently lacking. In this study, we investigated the quantitative proteome of N. annulifera venom, incorporating HPLC and LC-MS/MS to elucidate the venom toxicity. The immunoreactivities and in vivo neutralization activities of two hetero-specific antivenom products (Premium Serums Pan Africa polyvalent antivenom, PANAF and VINS African polyvalent antivenom, VAPAV) against the venom were subsequently examined. N. annulifera venom comprises 10 toxin families, with three-finger toxin (3FTx) being the most abundantly expressed (~78%). Within 3FTx, cytotoxin is the most dominant form and made up three-quarter of the venom bulk (~74%), whereas alpha-neurotoxins constitute <4% of the total venom proteins. Phospholipase A2 was undetected in the venom proteome, consistent with the unusual absence of PLA2 from the venoms of cobras in the Uraeus subgenus. In ELISA, PANAF and VAPAV showed comparable immunoreactivity toward the protein antigens of N. annulifera venom. These antivenoms, despite being raised against hetero-specific venoms, were capable of cross-neutralizing the lethal effect of N. annulifera venom in mice, with PANAF being marginally more potent.

Proteomics and antivenom immunoprofiling of Russell's viper (Daboia siamensis) venoms from Thailand and Indonesia

The Eastern Russell’s viper, Daboia siamensis, is a WHO Category 1 medically important venomous snake. It has a wide but disjunct distribution in Southeast Asia. The specific antivenom, D. siamensis Monovalent Antivenom (DsMAV-Thailand) is produced in Thailand but not available in Indonesia, where a heterologous trivalent antivenom, Serum Anti Bisa Ular (SABU), is used instead. This study aimed to investigate the geographical venom variation of D. siamensis from Thailand (Ds-Thailand) and Indonesia (Ds-Indonesia), and the immunorecognition of the venom proteins by antivenoms.

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.

Intraspecific sequence and gene expression variation contribute little to venom diversity in sidewinder rattlesnakes ( Crotalus cerastes)

Traits can evolve rapidly through changes in gene expression or protein-coding sequences. However, these forms of genetic variation can be correlated and changes to one can influence the other. As a result, we might expect traits lacking differential expression to preferentially evolve through changes in protein sequences or morphological adaptation. Given the lack of differential expression across the distribution of sidewinder rattlesnakes ( Crotalus cerastes), we tested this hypothesis by comparing the coding regions of genes expressed in the venom gland transcriptomes and fang morphology. We calculated Tajima's D and F_ST across four populations comparing toxin and nontoxin loci. Overall, we found little evidence of directional selection or differentiation between populations, suggesting that changes to protein sequences do not underlie the evolution of sidewinder venom or that toxins are under extremely variant selection pressures. Although low-expression toxins do not have higher sequence divergence between populations, they do have more standing variation on which selection can act. Additionally, we found significant differences in fang length among populations. The lack of differential expression and sequence divergence suggests sidewinders-given their generalist diet, moderate gene flow and environmental variation-are under stabilizing selection which functions to maintain a generalist phenotype. Overall, we demonstrate the importance of examining the relationship between gene expression and protein-coding changes to understand the evolution of complex traits.

Not so Dangerous After All? Venom Composition and Potency of the Pholcid (Daddy Long-Leg) Spider Physocyclus mexicanus

Pholcid spiders (Araneae: Pholcidae), officially "cellar spiders" but popularly known as "daddy long-legs," are renown for the potential of deadly toxic venom, even though venom composition and potency has never formally been studied. Here we detail the venom composition of male Physocyclus mexicanus using proteomic analyses and venom-gland transcriptomes ("venomics"). We also analyze the venom's potency on insects, and assemble available evidence regarding mammalian toxicity. The majority of the venom (51% of tryptic polypeptides and 62% of unique tryptic peptides) consists of proteins homologous to known venom toxins including enzymes (astacin metalloproteases, serine proteases and metalloendopeptidases, particularly neprilysins) and venom peptide neurotoxins. We identify 17 new groups of peptides (U1-17-PHTX) most of which are homologs of known venom peptides and are predicted to have an inhibitor cysteine knot fold; of these, 13 are confirmed in the proteome. Neprilysins (M13 peptidases), and astacins (M12 peptidases) are the most abundant venom proteins, respectively representing 15 and 11% of the individual proteins and 32 and 20% of the tryptic peptides detected in crude venom. Comparative evidence suggests that the neprilysin gene family is expressed in venoms across a range of spider taxa, but has undergone an expansion in the venoms of pholcids and may play a central functional role in these spiders. Bioassays of crude venoms on crickets resulted in an effective paralytic dose of 3.9 µg/g, which is comparable to that of crude venoms of Plectreurys tristis and other Synspermiata taxa. However, crickets exhibit flaccid paralysis and regions of darkening that are not observed after P. tristis envenomation. Documented bites on humans make clear that while these spiders can bite, the typical result is a mild sting with no long-lasting effects. Together, the evidence we present indicates pholcid venoms are a source of interesting new peptides and proteins, and effects of bites on humans and other mammals are inconsequential.

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.

The Dual Prey-Inactivation Strategy of Spiders-In-Depth Venomic Analysis of Cupiennius salei

Most knowledge of spider venom concerns neurotoxins acting on ion channels, whereas proteins and their significance for the envenomation process are neglected. The here presented comprehensive analysis of the venom gland transcriptome and proteome of Cupiennius salei focusses on proteins and cysteine-containing peptides and offers new insight into the structure and function of spider venom, here described as the dual prey-inactivation strategy. After venom injection, many enzymes and proteins, dominated by α-amylase, angiotensin-converting enzyme, and cysteine-rich secretory proteins, interact with main metabolic pathways, leading to a major disturbance of the cellular homeostasis. Hyaluronidase and cytolytic peptides destroy tissue and membranes, thus supporting the spread of other venom compounds. We detected 81 transcripts of neurotoxins from 13 peptide families, whereof two families comprise 93.7% of all cysteine-containing peptides. This raises the question of the importance of the other low-expressed peptide families. The identification of a venom gland-specific defensin-like peptide and an aga-toxin-like peptide in the hemocytes offers an important clue on the recruitment and neofunctionalization of body proteins and peptides as the origin of toxins.

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.

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

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

A Protein Decomplexation Strategy in Snake Venom Proteomics

Snake venoms are complex mixtures of proteins and peptides that play vital roles in the survival of venomous snakes. As with their diverse pharmacological activities, snake venoms can be highly variable, hence the importance of understanding the compositional details of different snake venoms. However, profiling venom protein mixtures is challenging, in particular when dealing with the diversity of protein subtypes and their abundances. Here we described an optimized strategy combining a protein decomplexation method with in-solution trypsin digestion and mass spectrometry of snake venom proteins. The approach involves the integrated use of C₁₈ reverse-phase high-performance liquid chromatography (RP-HPLC), sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), and nano-electrospray ionization tandem mass spectrometry (nano-ESI-LC-MS/MS).

Snake venom proteome and immuno-profiling of the hundred-pace viper, Deinagkistrodon acutus, in Taiwan

Deinagkistrodon acutus, also known as the hundred-pace viper or Chinese moccasin, is a clinically significant venomous snake in Taiwan. To address the lack of knowledge on the venom proteome of D. acutus, the venom composition was studied by a bottom-up proteomic approach combining reverse phase high-performance liquid chromatography, SDS-PAGE, and LC-MS/MS analysis. The immunoreactivity and cross-reactivity of Taiwanese freeze-dried D. acutus antivenom (DA-AV) and hemorrhagic antivenom (FH-AV) were investigated, as well. The proteomic analysis revealed the presence of 29 distinct proteins from D. acutus venom belonging to 8 snake venom protein families. Snake venom metalloproteinase (SVMP, 46.86%), C-type lectin (CLEC, 37.59%), phospholipase A₂ (PLA₂, 7.33%) and snake venom serine protease (SVSP, 6.62%) were the most abundant proteins. In addition to DA-AV, FH-AV also showed a profile of broad immunorecognition toward the venom of D. acutus. Remarkably, both antivenoms specifically reacted with the HPLC fractions containing SVMPs, and the titer was 5-10 times higher than fractions of other components. This information helps us to deeply understand the pathophysiology of D. acutus envenomation and guide us to development of more effective antivenom for clinical treatment.

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

The venom proteome of Hydrophis curtus (synonym: Lapemis hardwickii) from Penang, Malaysia was investigated with nano-electrospray ionization-liquid chromatography tandem mass spectrometry (ESI-LCMS/MS) of the reverse-phase high-performance liquid chromatography (HPLC) venom fractions. Thirty distinct protein forms were identified as toxins from ten families. The three major protein families were phospholipase A₂ (PLA₂, 62.0% of total venom proteins), three-finger toxin (3FTX, 26.33%) and cysteine-rich secretory protein (CRiSP, 9.00%). PLA₂ comprises diverse homologues (11 forms), predominantly the acidic subtypes (48.26%). 3FTX composed of one short alpha-neurotoxin (SNTX, 22.89%) and four long alpha-neurotoxins (LNTX, 3.44%). Both SNTX and LNTX were lethal in mice (intravenous LD₅₀ = 0.10 and 0.24 μg/g, respectively) but the PLA₂ were non-lethal (LD₅₀ >1 μg/g). The more abundant and toxic SNTX appeared to be the main driver of venom lethality (holovenom LD₅₀ = 0.20 μg/g). The heterologous Sea Snake Antivenom (SSAV, Australia) effectively cross-neutralized the venom (normalized potency = 9.35 mg venom neutralized per g antivenom) and the two neurotoxins in vivo, with the LNTX being neutralized more effectively (normalized potency = 3.5 mg toxin/g antivenom) than SNTX (normalized potency = 1.57 mg/g). SSAV immunorecognition was strong toward PLA₂ but moderate-to-weak toward the alpha-neurotoxins, indicating that neutralization of the alpha-neurotoxins should be further improved.

Unlocking the secrets of banded coral snake (Calliophis intestinalis, Malaysia): A venom with proteome novelty, low toxicity and distinct antigenicity

The Asiatic coral snakes are basal in the phylogeny of coral snakes. Although envenoming by the Asiatic coral snakes is rarely fatal, little is known about their venom properties and variability from the American coral snakes. Integrating reverse-phase high performance liquid chromatography and nano-liquid chromatography-tandem mass spectrometry, we showed that the venom proteome of the Malaysian banded or striped coral snake (Calliophis intestinalis) was composed of mainly phospholipases A₂ (PLA₂, 43.4%) and three-finger toxins (3FTx, 20.1%). Within 3FTx, the cytotoxins or cardiotoxins (CTX) dominated while the neurotoxins' content was much lower. Its subproteomic details contrasted with the 3FTx profile of most Micrurus sp., illustrating a unique dichotomy of venom phenotype between the Old and the New World coral snakes. Calliophis intestinalis venom proteome was correlated with measured enzymatic activities, and in vivo it was myotoxic but non-lethal to mice, frogs and geckos at high doses (5-10 μg/g). The venom contains species-specific toxins with distinct sequences and antigenicity, and the antibodies raised against PLA₂ and CTX of other elapids showed poor binding toward its venom antigens. The unique venom proteome of C. intestinalis unveiled a repertoire of novel toxins, and the toxicity test supported the need for post-bite monitoring of myotoxic complication. SIGNIFICANCE: Malaysian banded or striped coral snake (Calliophis intestinalis) has a cytotoxin (CTX)-predominating venom proteome, a characteristic shared by its congener, the Malayan blue coral snake (Calliophis bivirgata). With little neurotoxins (NTX), it illustrates a CTX/NTX dichotomy of venom phenotype between the Old World and the New World coral snakes. The low toxicity of the venom imply that C. intestinalis bite envenoming can be managed via symptomatic relief of the mild to moderate pain with appropriate analgesia. Systemically, the serum creatine kinase level of patients should be monitored serially for potential complication of myotoxicity. The distinct antigenicity of the venom proteins implies that the empirical use of heterologous antivenom is mostly inappropriate and not recommended.

Transcriptomic and Proteomic Analyses Reveal the Diversity of Venom Components from the Vaejovid Scorpion Serradigitus gertschi

To understand the diversity of scorpion venom, RNA from venomous glands from a sawfinger scorpion, Serradigitus gertschi, of the family Vaejovidae, was extracted and used for transcriptomic analysis. A total of 84,835 transcripts were assembled after Illumina sequencing. From those, 119 transcripts were annotated and found to putatively code for peptides or proteins that share sequence similarities with the previously reported venom components of other species. In accordance with sequence similarity, the transcripts were classified as potentially coding for 37 ion channel toxins; 17 host defense peptides; 28 enzymes, including phospholipases, hyaluronidases, metalloproteases, and serine proteases; nine protease inhibitor-like peptides; 10 peptides of the cysteine-rich secretory proteins, antigen 5, and pathogenesis-related 1 protein superfamily; seven La1-like peptides; and 11 sequences classified as "other venom components". A mass fingerprint performed by mass spectrometry identified 204 components with molecular masses varying from 444.26 Da to 12,432.80 Da, plus several higher molecular weight proteins whose precise masses were not determined. The LC-MS/MS analysis of a tryptic digestion of the soluble venom resulted in the de novo determination of 16,840 peptide sequences, 24 of which matched sequences predicted from the translated transcriptome. The database presented here increases our general knowledge of the biodiversity of venom components from neglected non-buthid scorpions.