Snake venomics of two poorly known Hydrophiinae: Comparative proteomics of the venoms of terrestrial Toxicocalamus longissimus and marine Hydrophis cyanocinctus

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.

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.

Polypharmacological Cell-Penetrating Peptides from Venomous Marine Animals Based on Immunomodulating, Antimicrobial, and Anticancer Properties

Complex pathological diseases, such as cancer, infection, and Alzheimer's, need to be targeted by multipronged curative. Various omics technologies, with a high rate of data generation, demand artificial intelligence to translate these data into druggable targets. In this study, 82 marine venomous animal species were retrieved, and 3505 cryptic cell-penetrating peptides (CPPs) were identified in their toxins. A total of 279 safe peptides were further analyzed for antimicrobial, anticancer, and immunomodulatory characteristics. Protease-resistant CPPs with endosomal-escape ability in Hydrophis hardwickii, nuclear-localizing peptides in Scorpaena plumieri, and mitochondrial-targeting peptides from Synanceia horrida were suitable for compartmental drug delivery. A broad-spectrum S. horrida-derived antimicrobial peptide with a high binding-affinity to bacterial membranes was an antigen-presenting cell (APC) stimulator that primes cytokine release and naïve T-cell maturation simultaneously. While antibiofilm and wound-healing peptides were detected in Synanceia verrucosa, APC epitopes as universal adjuvants for antiviral vaccination were in Pterois volitans and Conus monile. Conus pennaceus-derived anticancer peptides showed antiangiogenic and IL-2-inducing properties with moderate BBB-permeation and were defined to be a tumor-homing peptide (THP) with the ability to inhibit programmed death ligand-1 (PDL-1). Isoforms of RGD-containing peptides with innate antiangiogenic characteristics were in Conus tessulatus for tumor targeting. Inhibitors of neuropilin-1 in C. pennaceus are proposed for imaging probes or therapeutic delivery. A Conus betulinus cryptic peptide, with BBB-permeation, mitochondrial-targeting, and antioxidant capacity, was a stimulator of anti-inflammatory cytokines and non-inducer of proinflammation proposed for Alzheimer's. Conclusively, we have considered the dynamic interaction of cells, their microenvironment, and proportional-orchestrating-host- immune pathways by multi-target-directed CPPs resembling single-molecule polypharmacology. This strategy might fill the therapeutic gap in complex resistant disorders and increase the candidates' clinical-translation chance.

Bibliometric Analysis of Literature in Snake Venom-Related Research Worldwide (1933–2022)

Simple Summary

Around the world, snake envenomation poses a serious health risk. Proteins with pharmacological effects are present in snake venom. Recent studies elaborate snake venom and its potential application, including as a cancer drug and antibacterial substances. Our study aimed to analyze the global profile of the literature in snake venom research from documents indexed in the Scopus database between 1933 and 2022. In total, 2999 documents were published with Brazil showing the highest productivity. Antivenom, proteomics, and transcriptomics are emerging as hot topics on a global scale. The present study offers a distinctive overview of snake venom research conducted worldwide.

Abstract

Snake envenomation is a severe economic and health concern affecting countries worldwide. Snake venom carries a wide variety of small peptides and proteins with various immunological and pharmacological properties. A few key research areas related to snake venom, including its applications in treating cancer and eradicating antibiotic-resistant bacteria, have been gaining significant attention in recent years. The goal of the current study was to analyze the global profile of literature in snake venom research. This study presents a bibliometric review of snake venom-related research documents indexed in the Scopus database between 1933 and 2022. The overall number of documents published on a global scale was 2999, with an average annual production of 34 documents. Brazil produced the highest number of documents (n = 729), followed by the United States (n = 548), Australia (n = 240), and Costa Rica (n = 235). Since 1963, the number of publications has been steadily increasing globally. At a worldwide level, antivenom, proteomics, and transcriptomics are growing hot issues for research in this field. The current research provides a unique overview of snake venom research at global level from 1933 through 2022, and it may be beneficial in guiding future research.

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.

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.

Eye size in North American watersnakes (genus Nerodia) correlates with variation in feeding ecology

Visual acuity and sensitivity positively correlate to eye size in vertebrates, and eye size relates to the ecology of colubrid snakes. We investigated whether eye morphology of North American colubrids of the genus Nerodia correlates with ecology as well. Although all members of the genus utilize aquatic habits, they differ widely in the proportion of anurans they eat. We specifically tested whether eye size and placement is associated with the proportion of frogs in the diet to determine whether these two aspects of eye morphology relate to feeding ecology. Using phylogenetic comparative methods, we found a significantly positive association between eye size and the proportion of anurans eaten by Nerodia species. Although the evidence is equivocal, the anterior placement of relatively small eyes in one species may also enhance anurophagy. Although eye size may improve a snake’s ability to feed on frogs, eye size must compete with other selective forces on head shape in trade-offs that may also influence eye size.

Assessing the Binding of Venoms from Aquatic Elapids to the Nicotinic Acetylcholine Receptor Orthosteric Site of Different Prey Models

The evolution of an aquatic lifestyle from land dwelling venomous elapids is a radical ecological modification, bringing about many evolutionary changes from morphology to diet. Diet is an important ecological facet which can play a key role in regulating functional traits such as venom composition and prey-specific targeting of venom. In addition to predating upon novel prey (e.g., fish, fish eggs and invertebrates), the venoms of aquatic elapids also face the challenge of increased prey-escape potential in the aquatic environment. Thus, despite the independent radiation into an aquatic niche on four separate occasions, the venoms of aquatic elapids are evolving under convergent selection pressures. Utilising a biolayer interferometry binding assay, this study set out to elucidate whether crude venoms from representative aquatic elapids were target-specific to the orthosteric site of postsynaptic nicotinic acetylcholine receptor mimotopes of fish compared to other terrestrial prey types. Representatives of the four aquatic lineages were: aquatic coral snakes representative was Micrurus surinamensis;, sea kraits representative was Laticauda colubrina; sea snakes representatives were two Aipysurus spp. and eight Hydrophis spp; and water cobras representative was Naja annulata. No prey-specific differences in crude venom binding were observed from any species tested, except for Aipysurus laevis, which showed slight evidence of prey-potency differences. For Hydrophis caerulescens, H. peronii, H. schistosus and M. surinamensis, there was a lack of binding to the orthosteric site of any target lineage. Subsequent testing on the in vitro chick-biventer cervicis muscle preparation suggested that, while the venoms of these species bound postsynaptically, they bound to allosteric sites rather than orthosteric. Allosteric binding is potentially a weaker but faster-acting form of neurotoxicity and we hypothesise that the switch to allosteric binding is likely due to selection pressures related to prey-escape potential. This research has potentially opened up the possibility of a new functional class of toxins which have never been assessed previously while shedding light on the selection pressures shaping venom evolution.

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.

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.

Adaptive evolution of distinct prey-specific toxin genes in rear-fanged snake venom

Venom proteins evolve rapidly, and as a trophic adaptation are excellent models for predator-prey evolutionary studies. The key to a deeper understanding of venom evolution is an integrated approach, combining prey assays with analysis of venom gene expression and venom phenotype. Here, we use such an approach to study venom evolution in the Amazon puffing snake, Spilotes sulphureus, a generalist feeder. We identify two novel three-finger toxins: sulditoxin and sulmotoxin 1. These new toxins are not only two of the most abundant venom proteins, but are also functionally intriguing, displaying distinct prey-specific toxicities. Sulditoxin is highly toxic towards lizard prey, but is non-toxic towards mammalian prey, even at greater than 22-fold higher dosage. By contrast, sulmotoxin 1 exhibits the reverse trend. Furthermore, evolutionary analysis and structural modelling show highest sequence variability in the central loop of these proteins, probably driving taxon-specific toxicity. This is, to our knowledge, the first case in which a bimodal and contrasting pattern of toxicity has been shown for proteins in the venom of a single snake in relation to diet. Our study is an example of how toxin gene neofunctionalization can result in a venom system dominated by one protein superfamily and still exhibit flexibility in prey capture efficacy.

Spine-bellied sea snake (Hydrophis curtus) venom shows greater skeletal myotoxicity compared with cardiac myotoxicity

For the first time the impedance-based xCELLigence real-time cell analysis system was used to measure the myotoxicity of sea snake venom. With a focus on the spine-bellied sea snake (Hydrophis curtus), the venom of four sea snake species and three terrestrial snake species were compared for myotoxicity against a human skeletal muscle cell line (HSkMC). Hydrophis curtus venom was also tested on a human cardiac muscle cell line (HCM). Surprisingly, all four sea snake venoms tested on HSkMC produced an initial 100-280% rise in xCELLigence cell index that peaked within the first two hours before falling. The cell index rise of H. curtus venom was correlated with the WST-1 cell proliferation assay, which demonstrated an increase in mitochondrial metabolism. The myotoxicity of H. curtus was 4.7-8.2 fold less potent than the other sea snakes tested, the Australian beaked sea snake (Hydrophis zweifeli), the elegant sea snake (Hydrophis elegans) and the olive sea snake (Aipysurus laevis). If our cell-based results translate to H. curtus envenomations, this implies that H. curtus would be less myotoxic than the other three. Yet the myotoxicity of H. curtus venom to cardiac muscle cells was nine times weaker than for skeletal muscle cells, providing evidence that the venom has a selective effect on skeletal muscle cells. This evidence, combined with the slow-acting nature of the venom, supports a digestive role for sea snake myotoxins.

Synergistic strategies of predominant toxins in snake venoms

Synergism is a significant phenomenon present in snake venoms that may be an evolving strategy to potentiate toxicities. Synergism exists between different toxins or toxin complexes in various snake venoms, with phospholipaseA₂s (PLA₂s) (toxins or subunits) the main enablers. The predominant toxins, snake venom PLA₂s, metalloproteases (SVMPs), serine proteases (SVSPs) and three-finger toxins (3FTxs), play essential roles in synergistic processes. The hypothetical mechanisms of synergistic effect can be generalized under the effects of amplification and chaperoning. The Toxicity Score is among the few quantitative methods to assess synergism. Selection of toxins involved in synergistically enhanced toxicity as the targets are important for development of novel antivenoms or inhibitors.

The Venom of the Spine-Bellied Sea Snake (Hydrophis curtus): Proteome, Toxin Diversity and Intraspecific Variation

The spine-bellied sea snake (Hydrophis curtus) is known to cause human deaths, yet its venom composition has not yet been proteomically characterised. An in-depth proteomic analysis was performed on H. curtus venom from two different seasons, January and June, corresponding to adults and subadults, respectively. Venoms from adult and subadult H. curtus individuals were compared using reversed-phase high-performance liquid chromatography (RP-HPLC), matrix-assisted laser desorption ionisation-time of flight (MALDI-TOF) mass spectrometry and liquid chromatography electrospray ionisation mass spectrometry (LC-ESI-MS) to detect intraspecific variation, and the molecular weight data obtained with ESI-MS were used to assess toxin diversity. RP-HPLC and LC-ESI-MS/MS were used to characterise the venom proteome and estimate the relative abundances of protein families present. The most abundant protein family in January and June venoms is phospholipase A₂ (PLA₂: January 66.7%; June 54.5%), followed by three-finger toxins (3FTx: January 30.4%; June 40.4%) and a minor component of cysteine-rich secretory proteins (CRISP: January 2.5%; June 5%). Trace amounts of snake venom metalloproteinases (SVMP), C-type lectins and housekeeping and regulatory proteins were also found. Although the complexity of the venom is low by number of families present, each family contained a more diverse set of isoforms than previously reported, a finding that may have implications for the development of next-generation sea snake antivenoms. Intraspecific variability was shown to be minor with one obvious exception of a 14,157-Da protein that was present in some January (adult) venoms, but not at all in June (subadult) venoms. There is also a greater abundance of short-chain neurotoxins in June (subadult) venom compared with January (adult) venom. These differences potentially indicate the presence of seasonal, ontogenetic or sexual variation in H. curtus venom.

A Review and Database of Snake Venom Proteomes

Advances in the last decade combining transcriptomics with established proteomics methods have made possible rapid identification and quantification of protein families in snake venoms. Although over 100 studies have been published, the value of this information is increased when it is collated, allowing rapid assimilation and evaluation of evolutionary trends, geographical variation, and possible medical implications. This review brings together all compositional studies of snake venom proteomes published in the last decade. Compositional studies were identified for 132 snake species: 42 from 360 (12%) Elapidae (elapids), 20 from 101 (20%) Viperinae (true vipers), 65 from 239 (27%) Crotalinae (pit vipers), and five species of non-front-fanged snakes. Approximately 90% of their total venom composition consisted of eight protein families for elapids, 11 protein families for viperines and ten protein families for crotalines. There were four dominant protein families: phospholipase A₂s (the most common across all front-fanged snakes), metalloproteases, serine proteases and three-finger toxins. There were six secondary protein families: cysteine-rich secretory proteins, l-amino acid oxidases, kunitz peptides, C-type lectins/snaclecs, disintegrins and natriuretic peptides. Elapid venoms contained mostly three-finger toxins and phospholipase A₂s and viper venoms metalloproteases, phospholipase A₂s and serine proteases. Although 63 protein families were identified, more than half were present in <5% of snake species studied and always in low abundance. The importance of these minor component proteins remains unknown.

From Marine Venoms to Drugs: Efficiently Supported by a Combination of Transcriptomics and Proteomics

The potential of marine natural products to become new drugs is vast; however, research is still in its infancy. The chemical and biological diversity of marine toxins is immeasurable and as such an extraordinary resource for the discovery of new drugs. With the rapid development of next-generation sequencing (NGS) and liquid chromatography-tandem mass spectrometry (LC-MS/MS), it has been much easier and faster to identify more toxins and predict their functions with bioinformatics pipelines, which pave the way for novel drug developments. Here we provide an overview of related bioinformatics pipelines that have been supported by a combination of transcriptomics and proteomics for identification and function prediction of novel marine toxins.

How the Cobra Got Its Flesh-Eating Venom: Cytotoxicity as a Defensive Innovation and Its Co-Evolution with Hooding, Aposematic Marking, and Spitting

The cytotoxicity of the venom of 25 species of Old World elapid snake was tested and compared with the morphological and behavioural adaptations of hooding and spitting. We determined that, contrary to previous assumptions, the venoms of spitting species are not consistently more cytotoxic than those of closely related non-spitting species. While this correlation between spitting and non-spitting was found among African cobras, it was not present among Asian cobras. On the other hand, a consistent positive correlation was observed between cytotoxicity and utilisation of the defensive hooding display that cobras are famous for. Hooding and spitting are widely regarded as defensive adaptations, but it has hitherto been uncertain whether cytotoxicity serves a defensive purpose or is somehow useful in prey subjugation. The results of this study suggest that cytotoxicity evolved primarily as a defensive innovation and that it has co-evolved twice alongside hooding behavior: once in the Hemachatus + Naja and again independently in the king cobras (Ophiophagus). There was a significant increase of cytotoxicity in the Asian Naja linked to the evolution of bold aposematic hood markings, reinforcing the link between hooding and the evolution of defensive cytotoxic venoms. In parallel, lineages with increased cytotoxicity but lacking bold hood patterns evolved aposematic markers in the form of high contrast body banding. The results also indicate that, secondary to the evolution of venom rich in cytotoxins, spitting has evolved three times independently: once within the African Naja, once within the Asian Naja, and once in the Hemachatus genus. The evolution of cytotoxic venom thus appears to facilitate the evolution of defensive spitting behaviour. In contrast, a secondary loss of cytotoxicity and reduction of the hood occurred in the water cobra Naja annulata, which possesses streamlined neurotoxic venom similar to that of other aquatic elapid snakes (e.g., hydrophiine sea snakes). The results of this study make an important contribution to our growing understanding of the selection pressures shaping the evolution of snake venom and its constituent toxins. The data also aid in elucidating the relationship between these selection pressures and the medical impact of human snakebite in the developing world, as cytotoxic cobras cause considerable morbidity including loss-of-function injuries that result in economic and social burdens in the tropics of Asia and sub-Saharan Africa.

Venomics: integrative venom proteomics and beyond*

Venoms are integrated phenotypes that evolved independently in, and are used for predatory and defensive purposes by, a wide phylogenetic range of organisms. The same principles that contribute to the evolutionary success of venoms, contribute to making the study of venoms of great interest in such diverse fields as evolutionary ecology and biotechnology. Evolution is profoundly contingent, and nature also reinvents itself continuosly. Changes in a complex phenotypic trait, such as venom, reflect the influences of prior evolutionary history, chance events, and selection. Reconstructing the natural history of venoms, particularly those of snakes, which will be dealt with in more detail in this review, requires the integration of different levels of knowledge into a meaningful and comprehensive evolutionary framework for separating stochastic changes from adaptive evolution. The application of omics technologies and other disciplines have contributed to a qualitative and quantitative advance in the road map towards this goal. In this review we will make a foray into the world of animal venoms, discuss synergies and complementarities of the different approaches used in their study, and identify current bottlenecks that prevent inferring the evolutionary mechanisms and ecological constraints that molded snake venoms to their present-day variability landscape.

Sea snakes (Elapidae, Hydrophiinae) in their westernmost extent: an updated and illustrated checklist and key to the species in the Persian Gulf and Gulf of Oman

The Persian Gulf is known as the westernmost distribution limit for sea snakes, except for Hydrophis platurus (Linnaeus, 1766) that reaches southeastern Africa. Previous identification guides for sea snakes of the Persian Gulf and its adjacent waters in the Gulf of Oman were based on old data and confined mostly to written descriptions. Therefore, a series of field surveys were carried out in 2013 and 2014 through Iranian coastal waters of both gulfs to provide a comprehensive sampling of sea snakes in the area. This paper presents an illustrated and updated checklist and identification tool for sea snakes in the Persian Gulf and Gulf of Oman, which are based on new material and a review of the literature. This checklist includes ten species of marine hydrophiines, of which one, Microcephalophis cantoris (Günther, 1864), is a new record for the area. All specimens examined herein are deposited and available at the Zoological Museum of Shahid Bahonar University of Kerman, Kerman province, Iran.

Combined proteomics and transcriptomics identifies sting-related toxins of jellyfish Cyanea nozakii

Jellyfish sting has become a worldwide issue of critical concern to human health and safety in coastal areas in recent decades. Cyanea nozakii is one of the dominant blooming species and dangerous stingers in China. However, it remains unclear how many and what types of toxins are present in the venom. So, we used a combined transcriptomics and proteomics approach to investigate the venom composition of jellyfish C. nozakii. In total 4,608,524 Illumina valid reads were obtained to de novo assemble to 40,434 unigenes in the transcriptomics analysis. And, a total of 311,635 MS/MS spectra with 12,247 unique MS/MS spectra were generated to 1556 homologous proteins in the proteomics analysis. 174 potential toxin proteins were identified, with 27 proteins homology to the toxins from venomous animals, including phospholipase A2, zinc metalloproteinase-disintegrin agkistin, serine protease inhibitor, plancitoxin-1, alpha-latrocrustotoxin-Lt1a, etc. This study described the transcriptomics and venom proteomics of jellyfish C. nozakii for the first time. Our findings provide a comprehensive understanding of the venom composition of C. nozakii. Furthermore, the results may also be very helpful for the discovery of novel bioactive proteins, as well as the development of effective treatments for jellyfish sting in the future.

BIOLOGICAL SIGNIFICANCE

Jellyfish Cyanea nozakii is one of the most dangerous stingers in the coast of china. Hundreds of thousands of people would be stung every year and victims suffered a severe pain, itch, swelling, inflammation, wheal and even more serious consequence. However, it remains unclear how many and what types of toxins are present as well as the relationship between the clinical symptoms and toxins. Our combined transcriptomics and proteomics findings can provide a comprehensive understanding of the venom composition of C. nozakii and will also be helpful for the development of effective treatments for jellyfish sting in the future.

Full-Length Venom Protein cDNA Sequences from Venom-Derived mRNA: Exploring Compositional Variation and Adaptive Multigene Evolution

Envenomation of humans by snakes is a complex and continuously evolving medical emergency, and treatment is made that much more difficult by the diverse biochemical composition of many venoms. Venomous snakes and their venoms also provide models for the study of molecular evolutionary processes leading to adaptation and genotype-phenotype relationships. To compare venom complexity and protein sequences, venom gland transcriptomes are assembled, which usually requires the sacrifice of snakes for tissue. However, toxin transcripts are also present in venoms, offering the possibility of obtaining cDNA sequences directly from venom. This study provides evidence that unknown full-length venom protein transcripts can be obtained from the venoms of multiple species from all major venomous snake families. These unknown venom protein cDNAs are obtained by the use of primers designed from conserved signal peptide sequences within each venom protein superfamily. This technique was used to assemble a partial venom gland transcriptome for the Middle American Rattlesnake (Crotalus simus tzabcan) by amplifying sequences for phospholipases A₂, serine proteases, C-lectins, and metalloproteinases from within venom. Phospholipase A₂ sequences were also recovered from the venoms of several rattlesnakes and an elapid snake (Pseudechis porphyriacus), and three-finger toxin sequences were recovered from multiple rear-fanged snake species, demonstrating that the three major clades of advanced snakes (Elapidae, Viperidae, Colubridae) have stable mRNA present in their venoms. These cDNA sequences from venom were then used to explore potential activities derived from protein sequence similarities and evolutionary histories within these large multigene superfamilies. Venom-derived sequences can also be used to aid in characterizing venoms that lack proteomic profiles and identify sequence characteristics indicating specific envenomation profiles. This approach, requiring only venom, provides access to cDNA sequences in the absence of living specimens, even from commercial venom sources, to evaluate important regional differences in venom composition and to study snake venom protein evolution.

This work demonstrates that full-length venom protein messenger RNAs are present in secreted venoms and can be used to acquire full-length protein sequences of toxins from both front-fanged (Elapidae, Viperidae) and rear-fanged (Colubridae) snake venoms, eliminating the need to use venom glands. Full-length transcripts were obtained from venom samples that were fresh, newly lyophilized, old, field desiccated or commercially prepared, representing a significant advance over previous attempts which produced only partial sequence transcripts. Transcripts for all major venom protein families (metalloproteinases, serine proteases, C-type lectins, phospholipases A₂ and three-finger toxins) responsible for clinically significant snakebite symptoms were obtained from venoms. These sequences aid in the identification and characterization of venom proteome profiles, allowing for the identification of peptide sequences, specific isoforms, and novel venom proteins. The application of this technique will help to provide venom protein sequences for many snake species, including understudied rear-fanged snakes. Venom protein transcripts offer important insights into potential snakebite envenomation profiles and the molecular evolution of venom protein multigene families. By requiring only venom to obtain venom protein cDNAs, the approach detailed here will provide access to cDNA-based protein sequences from commercial and other venom sources, facilitating study of snake venom protein composition and evolution.

Venomic Analysis of the Poorly Studied Desert Coral Snake, Micrurus tschudii tschudii, Supports the 3FTx/PLA₂ Dichotomy across Micrurus Venoms

The venom proteome of the poorly studied desert coral snake Micrurus tschudii tschudii was unveiled using a venomic approach, which identified ≥38 proteins belonging to only four snake venom protein families. The three-finger toxins (3FTxs) constitute, both in number of isoforms (~30) and total abundance (93.6% of the venom proteome), the major protein family of the desert coral snake venom. Phospholipases A₂ (PLA₂s; seven isoforms, 4.1% of the venom proteome), 1–3 Kunitz-type proteins (1.6%), and 1–2 l-amino acid oxidases (LAO, 0.7%) complete the toxin arsenal of M. t. tschudii. Our results add to the growing evidence that the occurrence of two divergent venom phenotypes, i.e., 3FTx- and PLA₂-predominant venom proteomes, may constitute a general trend across the cladogenesis of Micrurus. The occurrence of a similar pattern of venom phenotypic variability among true sea snake (Hydrophiinae) venoms suggests that the 3FTx/PLA₂ dichotomy may be widely distributed among Elapidae venoms.

Ecological venomics: How genomics, transcriptomics and proteomics can shed new light on the ecology and evolution of venom

Animal venom is a complex cocktail of bioactive chemicals that traditionally drew interest mostly from biochemists and pharmacologists. However, in recent years the evolutionary and ecological importance of venom is realized as this trait has direct and strong influence on interactions between species. Moreover, venom content can be modulated by environmental factors. Like many other fields of biology, venom research has been revolutionized in recent years by the introduction of systems biology approaches, i.e., genomics, transcriptomics and proteomics. The employment of these methods in venom research is known as 'venomics'. In this review we describe the history and recent advancements of venomics and discuss how they are employed in studying venom in general and in particular in the context of evolutionary ecology. We also discuss the pitfalls and challenges of venomics and what the future may hold for this emerging scientific field.

Danger in the reef: Proteome, toxicity, and neutralization of the venom of the olive sea snake, Aipysurus laevis

Four specimens of the olive sea snake, Aipysurus laevis, were collected off the coast of Western Australia, and the venom proteome was characterized and quantitatively estimated by RP-HPLC, SDS-PAGE, and MALDI-TOF-TOF analyses. A. laevis venom is remarkably simple and consists of phospholipases A2 (71.2%), three-finger toxins (3FTx; 25.3%), cysteine-rich secretory proteins (CRISP; 2.5%), and traces of a complement control module protein (CCM; 0.2%). Using a Toxicity Score, the most lethal components were determined to be short neurotoxins. Whole venom had an intravenous LD50 of 0.07 mg/kg in mice and showed a high phospholipase A2 activity, but no proteinase activity in vitro. Preclinical assessment of neutralization and ELISA immunoprofiling showed that BioCSL Sea Snake Antivenom was effective in cross-neutralizing A. laevis venom with an ED50 of 821 μg venom per mL antivenom, with a binding preference towards short neurotoxins, due to the high degree of conservation between short neurotoxins from A. laevis and Enhydrina schistosa venom. Our results point towards the possibility of developing recombinant antibodies or synthetic inhibitors against A. laevis venom due to its simplicity.

Bioinformatics-Aided Venomics

Venomics is a modern approach that combines transcriptomics and proteomics to explore the toxin content of venoms. This review will give an overview of computational approaches that have been created to classify and consolidate venomics data, as well as algorithms that have helped discovery and analysis of toxin nucleic acid and protein sequences, toxin three-dimensional structures and toxin functions. Bioinformatics is used to tackle specific challenges associated with the identification and annotations of toxins. Recognizing toxin transcript sequences among second generation sequencing data cannot rely only on basic sequence similarity because toxins are highly divergent. Mass spectrometry sequencing of mature toxins is challenging because toxins can display a large number of post-translational modifications. Identifying the mature toxin region in toxin precursor sequences requires the prediction of the cleavage sites of proprotein convertases, most of which are unknown or not well characterized. Tracing the evolutionary relationships between toxins should consider specific mechanisms of rapid evolution as well as interactions between predatory animals and prey. Rapidly determining the activity of toxins is the main bottleneck in venomics discovery, but some recent bioinformatics and molecular modeling approaches give hope that accurate predictions of toxin specificity could be made in the near future.

Venomics of the beaked sea snake, Hydrophis schistosus: A minimalist toxin arsenal and its cross-neutralization by heterologous antivenoms

The venom proteome of Hydrophis schistosus (syn: Enhydrina schistosa) captured in Malaysian waters was investigated using reverse-phase HPLC, SDS-PAGE and high-resolution liquid chromatography-tandem mass spectrometry. The findings revealed a minimalist profile with only 18 venom proteins. These proteins belong to 5 toxin families: three-finger toxin (3FTx), phospholipase A2 (PLA2), cysteine-rich secretory protein (CRISP), snake venom metalloprotease (SVMP) and L-amino acid oxidase (LAAO). The 3FTxs (3 short neurotoxins and 4 long neurotoxins) constitute 70.5% of total venom protein, 55.8% being short neurotoxins and 14.7% long neurotoxins. The PLA2 family consists of four basic (21.4%) and three acidic (6.1%) isoforms. The minor proteins include one CRISP (1.3%), two SVMPs (0.5%) and one LAAO (0.2%). This is the first report of the presence of long neurotoxins, CRISP and LAAO in H. schistosus venom. The neurotoxins and the basic PLA2 are highly lethal in mice with an intravenous median lethal dose of <0.2 μg/g. Cross-neutralization by heterologous elapid antivenoms (Naja kaouthia monovalent antivenom and Neuro polyvalent antivenom) was moderate against the long neurotoxin and basic PLA2, but weak against the short neurotoxin, indicating that the latter is the limiting factor to be overcome for improving the antivenom cross-neutralization efficacy.

Interrogating the Venom of the Viperid Snake Sistrurus catenatus edwardsii by a Combined Approach of Electrospray and MALDI Mass Spectrometry

The complete sequence characterization of snake venom proteins by mass spectrometry is rather challenging due to the presence of multiple isoforms from different protein families. In the present study, we investigated the tryptic digest of the venom of the viperid snake Sistrurus catenatus edwardsii by a combined approach of liquid chromatography coupled to either electrospray (online) or MALDI (offline) mass spectrometry. These different ionization techniques proved to be complementary allowing the identification a great variety of isoforms of diverse snake venom protein families, as evidenced by the detection of the corresponding unique peptides. For example, ten out of eleven predicted isoforms of serine proteinases of the venom of S. c. edwardsii were distinguished using this approach. Moreover, snake venom protein families not encountered in a previous transcriptome study of the venom gland of this snake were identified. In essence, our results support the notion that complementary ionization techniques of mass spectrometry allow for the detection of even subtle sequence differences of snake venom proteins, which is fundamental for future structure-function relationship and possible drug design studies.

Jellyfish venomics and venom gland transcriptomics analysis of Stomolophus meleagris to reveal the toxins associated with sting

Jellyfish Stomolophus meleagris is a very dangerous animal because of its strong toxicity. However, the composition of the venom is still unclear. Both proteomics and transcriptomics approaches were applied in present study to investigate the major components and their possible relationships to the sting. The proteomics of the venom from S. meleagris was conducted by tryptic digestion of the crude venom followed by RP-HPLC separation and MS/MS analysis of the tryptic peptides. The venom gland transcriptome was analyzed using a high-throughput Illumina sequencing platform HiSeq 2000 with de novo assembly. A total of 218 toxins were identified including C-type lectin, phospholipase A₂ (PLA₂), potassium channel inhibitor, protease inhibitor, metalloprotease, hemolysin and other toxins, most of which should be responsible for the sting. Among them, serine protease inhibitor, PLA₂, potassium channel inhibitor and metalloprotease are predominant, representing 28.44%, 21.56%, 16.06% and 15.14% of the identified venom proteins, respectively. Overall, our combined proteomics and transcriptomics approach provides a systematic overview of the toxins in the venom of jellyfish S. meleagris and it will be significant to understand the mechanism of the sting.

BIOLOGICAL SIGNIFICANCE

Jellyfish Stomolophus meleagris is a very dangerous animal because of its strong toxicity. It often bloomed in the coast of China in recent years and caused thousands of people stung and even deaths every year. However, the components which caused sting are still unknown yet. In addition, no study about the venomics of jellyfish S. meleagris has been reported. In the present study, both proteomics and transcriptomics approaches were applied to investigate the major components related to the sting. The result showed that major component included C-type lectin, phospholipase A₂, potassium channel inhibitor, protease inhibitor, metalloprotease, hemolysin and other toxins, which should be responsible for the effect of sting. This is the first research about the venomics of jellyfish S. meleagris. It will be significant to understand the mechanism of the biological effects and helpful to develop ways to deal with the sting.

Two color morphs of the pelagic yellow-bellied sea snake, Pelamis platura, from different locations of Costa Rica: snake venomics, toxicity, and neutralization by antivenom

The yellow-bellied sea snake, Pelamis platura, is the most broadly distributed snake species. Despite being endowed with a highly lethal venom, a proteomic analysis of its toxin composition was unavailable. The venoms of specimens collected in Golfo de Papagayo and Golfo Dulce (Costa Rica), where two distinctive color morphs occur, were chromatographically compared. The latter inhabits a fjord-like gulf where the transit of oceanic sea snakes into and from the basin is restricted, thus possibly affecting gene flow. RP-HPLC evidenced a conserved venom protein profile in both populations, despite their divergent color phenotypes. Following a trend observed in other sea snakes, P. platura venom is relatively simple, being composed of proteins of the three-finger toxin (3FTx), phospholipase A2 (PLA2), cysteine-rich secretory protein (CRISP), 5'-nucleotidase, and metalloproteinase families. The first three groups represent 49.9%, 32.9%, and 9.1% of total venom protein, respectively. The most abundant component (~26%) is pelamitoxin (P62388), a short-chain 3FTx, followed by a major basic PLA2 (~20%) and a group of three isoforms of CRISPs (~9%). Whereas isolated pelamitoxin was highly lethal to mice, neither the PLA2 nor the CRISP fraction caused death. However, the PLA2 rapidly increased plasma creatine kinase activity after intramuscular injection, indicating its myotoxic action. Differing from myotoxic PLA2s of viperids, this PLA2 was not cytolytic to murine myogenic cells in vitro, suggesting possible differences in its mechanism of action. The median lethal dose (LD50) estimates for P. platura crude venom in mice and in three species of fishes did not differ significantly. The sea snake antivenom manufactured by CSL Ltd. (Australia), which uses Enhydrina schistosa as immunogen, cross-recognized the three major components of P. platura venom and, accordingly, neutralized the lethal activity of crude venom and pelamitoxin, therefore being of potential usefulness in the treatment of envenomations by this species.

BIOLOGICAL SIGNIFICANCE

Integrative analyses of animal venoms that combine the power of proteomics (venomics) with the characterization of their functional and immunological properties are significantly expanding knowledge on these remarkable bioweapons, both from a basic and a medical perspective. Costa Rica harbors a unique population of the yellow-bellied sea snake, Pelamis platura, that is restricted to a fjord-like gulf (Golfo Dulce). This population differs markedly from oceanic populations found elsewhere along the Pacific coast of this country, by presenting a patternless bright yellow coloration, instead of the typical bicolored or tricolored pattern of this species. It has been suggested that the dominance of this yellow-morph in Golfo Dulce might reflect gene flow restrictions, caused by the oceanographic conditions at this location. The present study demonstrates that the remarkable phenotypic variation between the two color morphs inhabiting Golfo Dulce and Golfo de Papagayo, respectively, is not associated with differences in the expression of venom components, as shown by their conserved RP-HPLC profiles. Proteomic analysis revealed the relatively simple toxin composition of P. platura venom, which contains three predominant types of proteins: three-finger toxins (protein abundance: 49.9%), phospholipases A2 (32.9%), and cysteine-rich secretory proteins (9.1%), together with few minor components. Further, the involvement of these most abundant proteins in the toxic effects of the venom, and their cross-recognition and neutralization by a sea snake antivenom produced against the venom of Enhydrina schistosa, were analyzed.

Next-generation snake venomics: protein-locus resolution through venom proteome decomplexation

Venom research has been continuously enhanced by technological advances. High-throughput technologies are changing the classical paradigm of hypothesis-driven research to technology-driven approaches. However, the thesis advocated in this paper is that full proteome coverage at locus-specific resolution requires integrating the best of both worlds into a protocol that includes decomplexation of the venom proteome prior to liquid chromatography-tandem mass spectrometry matching against a species-specific transcriptome. This approach offers the possibility of proof-checking the species-specific contig database using proteomics data. Immunoaffinity chromatography constitutes the basis of an antivenomics workflow designed to quantify the extent of cross-reactivity of antivenoms against homologous and heterologous venom toxins. In the author's view, snake venomics and antivenomics form part of a biology-driven conceptual framework to unveil the genesis and natural history of venoms, and their within- and between-species toxicological and immunological divergences and similarities. Understanding evolutionary trends across venoms represents the Rosetta Stone for generating broad-ranging polyspecific antivenoms.

Elapid snake venom analyses show the specificity of the peptide composition at the level of genera Naja and Notechis

Elapid snake venom is a highly valuable, but till now mainly unexplored, source of pharmacologically important peptides. We analyzed the peptide fractions with molecular masses up to 10 kDa of two elapid snake venoms—that of the African cobra, N. m. mossambica (genus Naja), and the Peninsula tiger snake, N. scutatus, from Kangaroo Island (genus Notechis). A combination of chromatographic methods was used to isolate the peptides, which were characterized by combining complimentary mass spectrometric techniques. Comparative analysis of the peptide compositions of two venoms showed specificity at the genus level. Three-finger (3-F) cytotoxins, bradykinin-potentiating peptides (BPPs) and a bradykinin inhibitor were isolated from the Naja venom. 3-F neurotoxins, Kunitz/basic pancreatic trypsin inhibitor (BPTI)-type inhibitors and a natriuretic peptide were identified in the N. venom. The inhibiting activity of the peptides was confirmed in vitro with a selected array of proteases. Cytotoxin 1 (P01467) from the Naja venom might be involved in the disturbance of cellular processes by inhibiting the cell 20S-proteasome. A high degree of similarity between BPPs from elapid and viperid snake venoms was observed, suggesting that these molecules play a key role in snake venoms and also indicating that these peptides were recruited into the snake venom prior to the evolutionary divergence of the snakes.

Sequencing-grade de novo analysis of MS/MS triplets (CID/HCD/ETD) from overlapping peptides

Full-length de novo sequencing of unknown proteins remains a challenging open problem. Traditional methods that sequence spectra individually are limited by short peptide length, incomplete peptide fragmentation, and ambiguous de novo interpretations. We address these issues by determining consensus sequences for assembled tandem mass (MS/MS) spectra from overlapping peptides (e.g., by using multiple enzymatic digests). We have combined electron-transfer dissociation (ETD) with collision-induced dissociation (CID) and higher-energy collision-induced dissociation (HCD) fragmentation methods to boost interpretation of long, highly charged peptides and take advantage of corroborating b/y/c/z ions in CID/HCD/ETD. Using these strategies, we show that triplet CID/HCD/ETD MS/MS spectra from overlapping peptides yield de novo sequences of average length 70 AA and as long as 200 AA at up to 99% sequencing accuracy.