Integration of transcriptomic and proteomic approaches for snake venom profiling

Dietary Breadth Predicts Toxin Expression Complexity in the Venoms of North American Gartersnakes

Selection on heritable phenotypic variation has played a prominent role in shaping the remarkable adaptations found across the Tree of Life. Complex ecological traits, such as snake venoms, are thought to be the products of selection because they directly link to fitness and survival. Snake venom increases the efficiency of prey capture and processing and is thus likely under intense selection. While many studies of snake venom have investigated the relationship between venom and diet, they have primarily focused on medically relevant front-fanged snakes. However, recent work has suggested that many non-front fanged snakes also rely on venom for subduing prey, despite having reduced toxicity toward humans. Here, we set out to uncover variation in toxin-producing genes, along with the ecological and evolutionary pressures impacting snake venom characteristics in the North American gartersnakes (Squamata: Natricidae: Thamnophis), a model group of non-front-fanged snakes. We annotated and analyzed Duvernoy's venom gland transcriptomes from 16 species representing all the major lineages within Thamnophis. We then generated measures of complexity of both toxins and dietary breadth. We found strong correlations between the complexity of toxin gene expression and phylogenetic diversity of diet, but no relationship between the complexity of the genetic makeup of the transcriptomes (allelic or sequence variation) and diet complexities. We also found phylogenetic signal associated with venom complexity, suggesting some influence of ancestry on venom characteristics. We suggest that, in non-front-fanged snakes, expression of toxins rather than sequence complexity is under strong selection by dietary diversity. These findings contradict similar studies from front-fanged snakes where increased transcriptomic complexity varies positively with dietary diversity, exposing a potential novel relationship between a complex phenotype-toxin expression-and its selective pressures-diet.

Progress and Challenges in the Field of Snakebite Envenoming Therapeutics

Snakebite envenoming kills and maims hundreds of thousands of people every year, especially in the rural settings of tropical regions. Envenomings are still treated with animal-derived antivenoms, which have prevented many lives from being lost but which are also medicines in need of innovation. Strides are being made to improve envenoming therapies, with promising efforts made toward optimizing manufacturing and quality aspects of existing antivenoms, accelerating research and development of recombinant antivenoms based on monoclonal antibodies, and repurposing of small-molecule inhibitors that block key toxins. Here, we review the most recent advances in these fields and discuss therapeutic opportunities and limitations for different snakebite treatment modalities. Finally, we discuss challenges related to preclinical and clinical evaluation, regulatory pathways, large-scale manufacture, and distribution and access that need to be addressed to fulfill the goals of the World Health Organization's global strategy to prevent and control snakebite envenoming.

V-ToCs (Venom Toxin Clustering): A tool for the investigation of sequence and structure similarities in snake venom toxins

Recently, there has been a major push toward the development of next-generation treatments against snakebite envenoming. However, unlike current antivenoms that rely on animal-derived polyclonal antibodies, most of these novel approaches are reliant on an in-depth understanding of the over 2000 known snake venom toxins. Indeed, by identifying similarities (i.e., conserved epitopes) across these different toxins, it is possible to design cross-reactive treatments, such as broadly-neutralising antibodies, that target these similarities. Therefore, in this project, we built an automated pipeline that generates sequence and structural distance matrices and homology trees across all available snake venom toxin sequences and structures. To facilitate analysis, we also developed a user-friendly and high-throughput visualisation tool, coined "Venom TOxin CluStering" (V-ToCs). This tool allows researchers to easily investigate sequence and structure patterns in snake venom toxins for a wide array of purposes, such as elucidating toxin evolution, and will also hopefully help guide the discovery and development of increasingly broadly-neutralising antivenoms in the near future.

Proteomic Analysis and Biochemical Characterization of the Nematocyst Extract of the Hydrozoan Velella velella

The venom contained within cnidarian nematocysts has a complex composition and holds significant potential for biotechnological applications. In this context, one of the most effective methods for studying nematocyst contents is the proteomic approach, which can detect even trace amounts of compounds while minimizing the need for large-scale animal collection, thus helping to preserve ecosystem integrity. This study aimed to provide a comprehensive proteomic and biochemical characterization of the crude nematocyst extract from the common hydrozoan Velella velella. Despite not being harmful to humans, the analysis of the crude venom extract from V. velella brought to the identification of 783 different proteins, categorized into structural components, enzymes, and potential toxins, revealing a qualitative composition of the venom similar to that of other more toxic cnidarians. Biochemical assays confirmed the presence of various active hydrolytic enzymes within the extract, including proteases, phospholipases, hyaluronidases, DNases, and chitinases. These findings pave the road for future studies involving the pharmacological applications of Velella velella venom components through recombinant production and functional testing.

Proteogenomics-guided functional venomics resolves the toxin arsenal and activity of Deinagkistrodon acutus venom

Snakebite primarily impacts rural communities of Africa, Asia, and Latin America. The sharp-nosed viper (Deinagkistrodon acutus) is among the snakes of highest medical importance in Asia. Despite various studies on its venom using modern venomics techniques, a comprehensive understanding of composition and function of this species' venom remains lacking. We combined proteogenomics with extensive bioactivity profiling to present the first genome-level catalogue of D. acutus venom proteins and their exochemistry. Our analysis identified an unusually simple venom containing 45 components from 20 distinct protein families. Relative toxin abundances indicate that C-type lectin and C-type lectin-related protein (CTL), snake venom metalloproteinase (svMP), snake venom serine protease (svSP), and phospholipase A2 (PLA2) constitute 90 % of the venom. Bioassays targeting key aspects of viperid envenomation showed considerable concentration-dependent cytotoxicity, particularly in kidney and lung cells, and potent protease and PLA2 activity. Factor Xa and thrombin activities were minor, and no plasmin activity was observed. Effects on haemolysis, intracellular calcium (Ca2+) release, and nitric oxide (NO) synthesis were negligible. Our analysis provides the first holistic genome-based overview of the toxin arsenal of D. acutus, predicting the molecular and functional basis of its life-threatening effects, and opens novel avenues for treating envenomation by this highly dangerous snake.

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

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

Acetylcholine-Binding Protein Affinity Profiling of Neurotoxins in Snake Venoms with Parallel Toxin Identification

Snakebite is considered a concerning issue and a neglected tropical disease. Three-finger toxins (3FTxs) in snake venoms primarily cause neurotoxic effects since they have high affinity for nicotinic acetylcholine receptors (nAChRs). Their small molecular size makes 3FTxs weakly immunogenic and therefore not appropriately targeted by current antivenoms. This study aims at presenting and applying an analytical method for investigating the therapeutic potential of the acetylcholine-binding protein (AChBP), an efficient nAChR mimic that can capture 3FTxs, for alternative treatment of elapid snakebites. In this analytical methodology, snake venom toxins were separated and characterised using high-performance liquid chromatography coupled with mass spectrometry (HPLC-MS) and high-throughput venomics. By subsequent nanofractionation analytics, binding profiling of toxins to the AChBP was achieved with a post-column plate reader-based fluorescence-enhancement ligand displacement bioassay. The integrated method was established and applied to profiling venoms of six elapid snakes (Naja mossambica, Ophiophagus hannah, Dendroaspis polylepis, Naja kaouthia, Naja haje and Bungarus multicinctus). The methodology demonstrated that the AChBP is able to effectively bind long-chain 3FTxs with relatively high affinity, but has low or no binding affinity towards short-chain 3FTxs, and as such provides an efficient analytical platform to investigate binding affinity of 3FTxs to the AChBP and mutants thereof and to rapidly identify bound toxins.

Comparing morphological and secretory aspects of cephalic glands among the New World coral snakes brings novel insights on their biological roles

Oral and other cephalic glands have been surveyed by several studies with distinct purposes. Despite the wide diversity and medical relevance of the New World coral snakes, studies focusing on understanding the biological roles of the glands within this group are still scarce. Specifically, the venom glands of some coral snakes were previously investigated but all other cephalic glands remain uncharacterized. In this sense, performing morphological and molecular analysis of these glands may help better understand their biological role. Here, we studied the morphology of the venom, infralabial, rictal, and harderian glands of thirteen species of Micrurus and Micruroides euryxanthus. We also performed a molecular characterization of these glands from selected species of Micrurus using transcriptomic and proteomic approaches. We described substantial morphological variation in the cephalic glands of New World coral snakes and structural evidence for protein-secreting cells in the inferior rictal glands. Our molecular analysis revealed that the venom glands, as expected, are majorly devoted to toxin production, however, the infralabial and inferior rictal glands also expressed some toxin genes at low to medium levels, despite the marked morphological differences. On the other hand, the harderian glands were dominated by the expression of lipocalins, but do not produce toxins. Our integrative analysis, including the prediction of biological processes and pathways, helped decipher some important traits of cephalic glands and better understand their biology.

miRNAs derived from cobra venom exosomes contribute to the cobra envenomation

Currently, there is an increasing amount of evidence indicating that exosomes and the miRNAs they contain are crucial players in various biological processes. However, the role of exosomes and miRNAs in snake venom during the envenomation process remains largely unknown. In this study, fresh venom from Naja atra of different ages (2-month-old, 1-year-old, and 5-year-old) was collected, and exosomes were isolated through ultracentrifugation. The study found that exosomes with inactivated proteins and enzymes can still cause symptoms similar to cobra envenomation, indicating that substances other than proteins and enzymes in exosomes may also play an essential role in cobra envenomation. Furthermore, the expression profiles of isolated exosome miRNAs were analyzed. The study showed that a large number of miRNAs were co-expressed and abundant in cobra venom exosomes (CV-exosomes) of different ages, including miR-2904, which had high expression abundance and specific sequences. The specific miR-2094 derived from CV-exosomes (CV-exo-miR-2904) was overexpressed both in vitro and in vivo. As a result, CV-exo-miR-2904 induced symptoms similar to cobra envenomation in mice and caused liver damage, demonstrating that it plays a crucial role in cobra envenomation. These results reveal that CV-exosomes and the miRNAs they contain play a significant regulatory role in cobra envenomation. Our findings provide new insights for the treatment of cobra bites and the development of snake venom-based medicines.

Evolution of Three-Finger Toxin Genes in Neotropical Colubrine Snakes (Colubridae)

Snake venom research has historically focused on front-fanged species (Viperidae and Elapidae), limiting our knowledge of venom evolution in rear-fanged snakes across their ecologically diverse phylogeny. Three-finger toxins (3FTxs) are a known neurotoxic component in the venoms of some rear-fanged snakes (Colubridae: Colubrinae), but it is unclear how prevalent 3FTxs are both in expression within venom glands and more broadly among colubrine species. Here, we used a transcriptomic approach to characterize the venom expression profiles of four species of colubrine snakes from the Neotropics that were dominated by 3FTx expression (in the genera Chironius, Oxybelis, Rhinobothryum, and Spilotes). By reconstructing the gene trees of 3FTxs, we found evidence of putative novel heterodimers in the sequences of Chironius multiventris and Oxybelis aeneus, revealing an instance of parallel evolution of this structural change in 3FTxs among rear-fanged colubrine snakes. We also found positive selection at sites within structural loops or "fingers" of 3FTxs, indicating these areas may be key binding sites that interact with prey target molecules. Overall, our results highlight the importance of exploring the venoms of understudied species in reconstructing the full evolutionary history of toxins across the tree of life.

A Review on Genotoxic and Genoprotective Effects of Biologically Active Compounds of Animal Origin

Envenomation by animal venoms remains a serious medical and social problem, especially in tropical countries. On the other hand, animal venoms are widely used as a source of biologically active compounds for the development of novel drugs. Numerous derivatives of animal venoms are already used in clinical practice. When analysing the mechanisms of action of animal venoms, attention is usually focused on the main target of the venom's enzymes and peptides such as neurotoxic, cytotoxic or haemorrhagic effects. In the present review, we would like to draw attention to the "hidden" effects of animal venoms and their derivatives in regard to DNA damage and/or protection against DNA damage. Alkaloids and terpenoids isolated from sponges such as avarol, ingenamine G or variolin B manifest the capability to bind DNA in vitro and produce DNA breaks. Trabectidin, isolated from a sea squirt, also binds and damages DNA. A similar action is possible for peptides isolated from bee and wasp venoms such as mastoparan, melectin and melittin. However, DNA lesions produced by the crude venoms of jellyfish, scorpions, spiders and snakes arise as a consequence of cell membrane damage and the subsequent oxidative stress, whereas certain animal venoms or their components produce a genoprotective effect. Current research data point to the possibility of using animal venoms and their components in the development of various potential therapeutic agents; however, before their possible clinical use the route of injection, molecular target, mechanism of action, exact dosage, possible side effects and other fundamental parameters should be further investigated.

A secretory system for extracellular production of spider neurotoxin huwentoxin-I in Escherichia coli

Due to their advantages in structural stability and versatility, cysteine-rich peptides, which are secreted from the venom glands of venomous animals, constitute a naturally occurring pharmaceutical arsenal. However, the correct folding of disulfide bonds is a challenging task in the prokaryotic expression system like Escherichia coli due to the reducing environment. Here, a secretory expression plasmid pSE-G1M5-SUMO-HWTX-I for the spider neurotoxin huwentoxin-I (HWTX-I) with three disulfides as a model of cysteine-rich peptides was constructed. By utilizing the signal peptide G1M5, the fusion protein 6 × His-SUMO-HWTX-I was successfully secreted into extracellular medium of BL21(DE3). After enrichment using cation-exchange chromatography and purification utilizing the Ni-NTA column, 6 × His-SUMO-HWTX-I was digested via Ulp1 kinase to release recombinant HWTX-I (rHWTX-I), which was further purified utilizing RP-HPLC. Finally, both impurities with low and high molecular weights were completely removed. The molecular mass of rHWTX-I was identified as being 3750.8 Da, which was identical to natural HWTX-I with three disulfide bridges. Furthermore, by utilizing whole-cell patch clamp, the sodium currents of hNav1.7 could be inhibited by rHWTX-I and the IC50 value was 419 nmol/L.