Evolution of an Arsenal

Sugar-coated survival: N-glycosylation as a unique bearded dragon venom resistance trait within Australian agamid lizards

In the ongoing evolutionary arms race between predators and prey, adaptive innovations often trigger a reciprocal response. For instance, the emergence of α-neurotoxins in snake venom has driven prey species targeted by these snakes to evolve sophisticated defense mechanisms. This study zeroes in on the particular motifs within the orthosteric sites of post-synaptic nicotinic acetylcholine receptors (nAChR) that confer resistance to α-neurotoxins, often through structural alterations of nAChR. This research examined Australian agamid lizards, a primary prey group for Australian elapid snakes, which are subject to predatory selection pressures. We previously showed that Pogona vitticeps (Central bearded dragon) was resistant to α-neurotoxic snake venoms through a steric hindrance form resistance evolving within the nAChR orthosteric, specifically through the 187-189NVT motif resulting in the presence of N-glycosylation, with the branching carbohydrate chains impeding the binding by the neurotoxins. This adaptive trait is thought to be a compensatory mechanism for the lizard's limited escape capabilities. Despite the significance of this novel adaptation, the prevalence and evolutionary roots of such venom resistance in Australian agamids have not been thoroughly investigated. To fill this knowledge gap, we undertook a comprehensive sequencing analysis of the nAChR ligand-binding domain across the full taxonomical diversity of Australian agamid species. Our findings reveal that the N-glycosylation resistance mechanism is a trait unique to the Pogona genus and absent in other Australian agamids. This aligns with Pogona's distinctive morphology, which likely increases vulnerability to neurotoxic elapid snakes, thereby increasing selective pressures for resistance. In contrast, biolayer interferometry experiments with death adder (Acanthophis species) venoms did not indicate any resistance-related binding patterns in other agamids, suggesting a lack of similar resistance adaptations, consistent with these lineages either being fast-moving, covered with large defensive spines, or being arboreal. This research not only uncovers a novel α-neurotoxin resistance mechanism in Australian agamids but also highlights the complex dynamics of the predator-prey chemical arms race. It provides a deeper understanding of how evolutionary pressures shape the interactions between venomous snakes and their prey.

Highly conserved and extremely variable: The paradoxical pattern of toxin expression revealed by comparative venom-gland transcriptomics of Phalotris (Serpentes: Dipsadidae)

Although non-front fanged snakes account for almost two-thirds of snake diversity, most studies on venom composition and evolution focus exclusively on front-fanged species, which comprise most of the clinically relevant accidents. Comprehensive reports on venom composition of non-front fanged snakes are still scarce for several groups. In this study, we address such shortage of knowledge by providing new insights about the venom composition among species of Phalotris, a poorly studied Neotropical dipsadid genus. Phalotris are known for their specialized venom delivery system and toxic venoms, which can cause life-threatening accidents in humans. We evaluate the venom-gland transcriptome of Phalotris, comparing the following three South American species: P. reticulatus for the Araucaria Pine forests, P. lemniscatus for the Pampa grasslands, and P. mertensi for the Brazilian Cerrado. Our results indicate similar venom profiles, in which they share a high expression level of Kunitz-type inhibitors (KUNZ). On the other hand, comparative analyses revealed substantial differences in the expression levels of C-type lectins (CTL) and snake venom metalloproteinases (SVMP). The diverse set of SVMP and CTL isoforms shows signals of positive selection, and we also identified truncated forms of type III SVMPs, which resemble type II and type I SVMPs of viperids. Additionally, we identified a CNP precursor hosting a proline-rich region containing a BPP motif resembling those commonly detected in viperid venoms with hypotensive activity. Altogether, our results suggest an evolutionary history favoring high expression levels of few KUNZ isoforms in Phalotris venoms, contrasting with a highly diverse set of SVMP and CTL isoforms. Such diversity can be comparable with the venom variability observed in some viperids. Our findings highlight the extreme phenotypic diversity of non-front fanged snakes and the importance to allocate greater effort to study neglected groups of Colubroidea.

Unveiling Novel Kunitz- and Waprin-Type Toxins in the Micrurus mipartitus Coral Snake Venom Gland: An In Silico Transcriptome Analysis

Kunitz-type peptide expression has been described in the venom of snakes of the Viperidae, Elapidae and Colubridae families. This work aimed to identify these peptides in the venom gland transcriptome of the coral snake Micrurus mipartitus. Transcriptomic analysis revealed a high diversity of venom-associated Kunitz serine protease inhibitor proteins (KSPIs). A total of eight copies of KSPIs were predicted and grouped into four distinctive types, including short KSPI, long KSPI, Kunitz-Waprin (Ku-WAP) proteins, and a multi-domain Kunitz-type protein. From these, one short KSPI showed high identity with Micrurus tener and Austrelaps superbus. The long KSPI group exhibited similarity within the Micrurus genus and showed homology with various elapid snakes and even with the colubrid Pantherophis guttatus. A third group suggested the presence of Kunitz domains in addition to a whey-acidic-protein-type four-disulfide core domain. Finally, the fourth group corresponded to a transcript copy with a putative 511 amino acid protein, formerly annotated as KSPI, which UniProt classified as SPINT1. In conclusion, this study showed the diversity of Kunitz-type proteins expressed in the venom gland transcriptome of M. mipartitus.

The role of venom proteomics and single-domain antibodies for antivenoms: Progress in snake envenoming treatment

Single-domain antibodies (sdAbs) hold promise for developing new biopharmaceuticals to treat neglected tropical diseases (NTDs), including snakebites, which are severe and occur frequently. In addition, limitations of conventional snakebite treatments, especially in terms of local action, and the global antivenom crisis incentivize the use of this biotechnological tool to design next-generation snakebite antivenoms. Conventional antivenoms for snakebite treatment are usually composed of immunoglobulin G or F(ab')2 fragments derived from the plasma of immunized animals. sdAbs, the smallest antigen-binding fragments, are derived from the variable domains of camelid heavy-chain antibodies. sdAbs may have some advantages over conventional antivenoms for local toxicity, such as better penetration into tissues due to their small size, and high solubility and affinity for venom antigens due to their unique antigen-binding loops and ability to access cryptic epitopes. We present an overview of current antivenom therapy in the context of sdAb development for toxin neutralization. Furthermore, strategies are presented for identifying snake venom's major toxins as well as for developing antisnake toxin sdAbs by employing proteomic tools for toxin neutralization.

High-Voltage Toxin'Roll: Electrostatic Charge Repulsion as a Dynamic Venom Resistance Trait in Pythonid Snakes

The evolutionary interplay between predator and prey has significantly shaped the development of snake venom, a critical adaptation for subduing prey. This arms race has spurred the diversification of the components of venom and the corresponding emergence of resistance mechanisms in the prey and predators of venomous snakes. Our study investigates the molecular basis of venom resistance in pythons, focusing on electrostatic charge repulsion as a defense against α-neurotoxins binding to the alpha-1 subunit of the postsynaptic nicotinic acetylcholine receptor. Through phylogenetic and bioactivity analyses of orthosteric site sequences from various python species, we explore the prevalence and evolution of amino acid substitutions that confer resistance by electrostatic repulsion, which initially evolved in response to predatory pressure by Naja (cobra) species (which occurs across Africa and Asia). The small African species Python regius retains the two resistance-conferring lysines (positions 189 and 191) of the ancestral Python genus, conferring resistance to sympatric Naja venoms. This differed from the giant African species Python sebae, which has secondarily lost one of these lysines, potentially due to its rapid growth out of the prey size range of sympatric Naja species. In contrast, the two Asian species Python brongersmai (small) and Python bivittatus (giant) share an identical orthosteric site, which exhibits the highest degree of resistance, attributed to three lysine residues in the orthosteric sites. One of these lysines (at orthosteric position 195) evolved in the last common ancestor of these two species, which may reflect an adaptive response to increased predation pressures from the sympatric α-neurotoxic snake-eating genus Ophiophagus (King Cobras) in Asia. All these terrestrial Python species, however, were less neurotoxin-susceptible than pythons in other genera which have evolved under different predatory pressure as: the Asian species Malayopython reticulatus which is arboreal as neonates and juveniles before rapidly reaching sizes as terrestrial adults too large for sympatric Ophiophagus species to consider as prey; and the terrestrial Australian species Aspidites melanocephalus which occupies a niche, devoid of selection pressure from α-neurotoxic predatory snakes. Our findings underline the importance of positive selection in the evolution of venom resistance and suggest a complex evolutionary history involving both conserved traits and secondary evolution. This study enhances our understanding of the molecular adaptations that enable pythons to survive in environments laden with venomous threats and offers insights into the ongoing co-evolution between venomous snakes and their prey.

The macroevolutionary singularity of snakes

Snakes and lizards (Squamata) represent a third of terrestrial vertebrates and exhibit spectacular innovations in locomotion, feeding, and sensory processing. However, the evolutionary drivers of this radiation remain poorly known. We infer potential causes and ultimate consequences of squamate macroevolution by combining individual-based natural history observations (>60,000 animals) with a comprehensive time-calibrated phylogeny that we anchored with genomic data (5400 loci) from 1018 species. Due to shifts in the dynamics of speciation and phenotypic evolution, snakes have transformed the trophic structure of animal communities through the recurrent origin and diversification of specialized predatory strategies. Squamate biodiversity reflects a legacy of singular events that occurred during the early history of snakes and reveals the impact of historical contingency on vertebrate biodiversity.

Optimizing drug discovery for snakebite envenoming via a high-throughput phospholipase A2 screening platform

Snakebite envenoming is a neglected tropical disease that causes as many as 1.8 million envenomings and 140,000 deaths annually. To address treatment limitations that exist with current antivenoms, the search for small molecule drug-based inhibitors that can be administered as early interventions has recently gained traction. Snake venoms are complex mixtures of proteins, peptides and small molecules and their composition varies substantially between and within snake species. The phospholipases A2 (PLA2) are one of the main pathogenic toxin classes found in medically important viper and elapid snake venoms, yet varespladib, a drug originally developed for the treatment of acute coronary syndrome, remains the only PLA2 inhibitor shown to effectively neutralise venom toxicity in vitro and in vivo, resulting in an extremely limited drug portfolio. Here, we describe a high-throughput drug screen to identify novel PLA2 inhibitors for repurposing as snakebite treatments. We present method optimisation of a 384-well plate, colorimetric, high-throughput screening assay that allowed for a throughput of ∼2,800 drugs per day, and report on the screening of a ∼3,500 post-phase I repurposed drug library against the venom of the Russell's viper, Daboia russelii. We further explore the broad-spectrum inhibitory potential and efficacy of the resulting top hits against a range of medically important snake venoms and demonstrate the utility of our method in determining drug EC50s. Collectively, our findings support the future application of this method to fully explore the chemical space to discover novel PLA2-inhibiting drugs of value for preventing severe pathology caused by snakebite envenoming.

A Review of Rattlesnake Venoms

Venom components are invaluable in biomedical research owing to their specificity and potency. Many of these components exist in two genera of rattlesnakes, Crotalus and Sistrurus, with high toxicity and proteolytic activity variation. This review focuses on venom components within rattlesnakes, and offers a comparison and itemized list of factors dictating venom composition, as well as presenting their known characteristics, activities, and significant applications in biosciences. There are 64 families and subfamilies of proteins present in Crotalus and Sistrurus venom. Snake venom serine proteases (SVSP), snake venom metalloproteases (SVMP), and phospholipases A2 (PLA2) are the standard components in Crotalus and Sistrurus venom. Through this review, we highlight gaps in the knowledge of rattlesnake venom; there needs to be more information on the venom composition of three Crotalus species and one Sistrurus subspecies. We discuss the activity and importance of both major and minor components in biomedical research and drug development.

Resistance Is Not Futile: Widespread Convergent Evolution of Resistance to Alpha-Neurotoxic Snake Venoms in Caecilians (Amphibia: Gymnophiona)

Predatory innovations impose reciprocal selection pressures upon prey. The evolution of snake venom alpha-neurotoxins has triggered the corresponding evolution of resistance in the post-synaptic nicotinic acetylcholine receptors of prey in a complex chemical arms race. All other things being equal, animals like caecilians (an Order of legless amphibians) are quite vulnerable to predation by fossorial elapid snakes and their powerful alpha-neurotoxic venoms; thus, they are under strong selective pressure. Here, we sequenced the nicotinic acetylcholine receptor alpha-1 subunit of 37 caecilian species, representing all currently known families of caecilians from across the Americas, Africa, and Asia, including species endemic to the Seychelles. Three types of resistance were identified: (1) steric hindrance from N-glycosylated asparagines; (2) secondary structural changes due to the replacement of proline by another amino acid; and (3) electrostatic charge repulsion of the positively charged neurotoxins, through the introduction of a positively charged amino acid into the toxin-binding site. We demonstrated that resistance to alpha-neurotoxins convergently evolved at least fifteen times across the caecilian tree (three times in Africa, seven times in the Americas, and five times in Asia). Additionally, as several species were shown to possess multiple resistance modifications acting synergistically, caecilians must have undergone at least 20 separate events involving the origin of toxin resistance. On the other hand, resistance in non-caecilian amphibians was found to be limited to five origins. Together, the mutations underlying resistance in caecilians constitute a robust signature of positive selection which strongly correlates with elapid presence through both space (sympatry with caecilian-eating elapids) and time (Cenozoic radiation of elapids). Our study demonstrates the extent of convergent evolution that can be expected when a single widespread predatory adaptation triggers parallel evolutionary arms races at a global scale.

The Venom Composition of the Snake Tribe Philodryadini: 'Omic' Techniques Reveal Intergeneric Variability among South American Racers

Snakes of the Philodryadini tribe are included in the Dipsadidae family, which is a diverse group of rear-fanged snakes widespread in different ecological conditions, including habitats and diet. However, little is known about the composition and effects of their venoms despite their relevance for understanding the evolution of these snakes or even their impact on the occasional cases of human envenoming. In this study, we integrated venom gland transcriptomics, venom proteomics and functional assays to characterize the venoms from eight species of the Philodryadini tribe, which includes the genus Philodryas, Chlorosoma and Xenoxybelis. The most abundant components identified in the venoms were snake venom metalloproteinases (SVMPs), cysteine-rich secretory proteins (CRISPs), C-type lectins (CTLs), snake endogenous matrix metalloproteinases type 9 (seMMP-9) and snake venom serinoproteinases (SVSPs). These protein families showed a variable expression profile in each genus. SVMPs were the most abundant components in Philodryas, while seMMP-9 and CRISPs were the most expressed in Chlorosoma and Xenoxybelis, respectively. Lineage-specific differences in venom composition were also observed among Philodryas species, whereas P. olfersii presented the highest amount of SVSPs and P. agassizii was the only species to express significant amounts of 3FTx. The variability observed in venom composition was confirmed by the venom functional assays. Philodryas species presented the highest SVMP activity, whereas Chlorosoma species showed higher levels of gelatin activity, which may correlate to the seMMP-9 enzymes. The variability observed in the composition of these venoms may be related to the tribe phylogeny and influenced by their diets. In the presented study, we expanded the set of venomics studies of the Philodryadini tribe, which paves new roads for further studies on the evolution and ecology of Dipsadidae snakes.

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.

Rhomb-I, a P–I metalloproteinase from Lachesis muta rhombeata venom degrades vessel extra cellular matrix components and impairs platelet aggregation

Rhomb-I, a 23-kDa metalloproteinase was isolated from L. m. rhombeata venom. Its dimethylcasein proteolysis was abolished by metal chelators, and slightly enhanced by Ca²⁺ and Mg²⁺ ions, but inhibited by Co²⁺, Zn²⁺ and α2-macroglobulin. In aqueous solution, rhomb-I autoproteolyzed to a 20- and 11-kDa fragments at 37 °C. The amino acid sequence showed high homology with other snake venom metalloproteinases. Rhomb-I causes hemorrhage that may be ascribed to hydrolysis of essential basement membrane, extracellular matrix and plasma proteins. It preferentially cleaves the α-chains of fibrin (ogen). Rhomb-I inhibited convulxin- and von Willebrand factor (vWF)-induced aggregation on human platelets without significant effect on collagen-stimulated aggregation or other effectors. It digests vWF into a low-molecular-mass multimers of vWF and a rvWF-A1 domain to a 27-kDa fragment as revealed by western blotting with mouse anti-rvWF A1-domain IgG. Incubation of platelets with rhomb-I resulted in adhesion to and cleavage of platelet receptors glycoprotein (GP)Ibα and GPVI to release a 55-kDa soluble form. Both membrane glycoproteins GPIbα that binds vWF, together with GPVI which binds collagen, play a key role in mediating platelet adhesion/activation and can initiate (patho)physiological thrombus formation. Conclusions: rhomb-I is implicated in the pathophysiology of Lachesis envenoming by disrupting vasculature, hemostasis and platelet aggregation through impairing vWF-GPIb axis and blocking GPVI-collagen binding.

Metabolome-Based Classification of Snake Venoms by Bioinformatic Tools

Snakebite is considered a neglected tropical disease, and it is one of the most intricate ones. The variability found in snake venom is what makes it immensely complex to study. These variations are present both in the big and the small molecules found in snake venom. This study focused on examining the variability found in the venom's small molecules (i.e., mass range of 100-1000 Da) between two main families of venomous snakes-Elapidae and Viperidae-managing to create a model able to classify unknown samples by means of specific features, which can be extracted from their LC-MS data and output in a comprehensive list. The developed model also allowed further insight into the composition of snake venom by highlighting the most relevant metabolites of each group by clustering similarly composed venoms. The model was created by means of support vector machines and used 20 features, which were merged into 10 principal components. All samples from the first and second validation data subsets were correctly classified. Biological hypotheses relevant to the variation regarding the metabolites that were identified are also given.

Spatial Venomics─Cobra Venom System Reveals Spatial Differentiation of Snake Toxins by Mass Spectrometry Imaging

Among venomous animals, toxic secretions have evolved as biochemical weapons associated with various highly specialized delivery systems on many occasions. Despite extensive research, there is still limited knowledge of the functional biology of most animal toxins, including their venom production and storage, as well as the morphological structures within sophisticated venom producing tissues that might underpin venom modulation. Here, we report on the spatial exploration of a snake venom gland system by matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI), in combination with standard proteotranscriptomic approaches, to enable in situ toxin mapping in spatial intensity maps across a venom gland sourced from the Egyptian cobra (Naja haje). MALDI-MSI toxin visualization on the elapid venom gland reveals a high spatial heterogeneity of different toxin classes at the proteoform level, which may be the result of physiological constraints on venom production and/or storage that reflects the potential for venom modulation under diverse stimuli.

Extracellular Vesicles from Bothrops jararaca Venom Are Diverse in Structure and Protein Composition and Interact with Mammalian Cells

Snake venoms are complex cocktails of non-toxic and toxic molecules that work synergistically for the envenoming outcome. Alongside the immediate consequences, chronic manifestations and long-term sequelae can occur. Recently, extracellular vesicles (EVs) were found in snake venom. EVs mediate cellular communication through long distances, delivering proteins and nucleic acids that modulate the recipient cell's function. However, the biological roles of snake venom EVs, including possible cross-organism communication, are still unknown. This knowledge may expand the understanding of envenoming mechanisms. In the present study, we isolated and characterized the EVs from Bothrops jararaca venom (Bj-EVs), giving insights into their biological roles. Fresh venom was submitted to differential centrifugation, resulting in two EV populations with typical morphology and size range. Several conserved EV markers and a subset of venom related EV markers, represented mainly by processing enzymes, were identified by proteomic analysis. The most abundant protein family observed in Bj-EVs was 5'-nucleotidase, known to be immunosuppressive and a low abundant and ubiquitous toxin in snake venoms. Additionally, we demonstrated that mammalian cells efficiently internalize Bj-EVs. The commercial antibothropic antivenom partially recognizes Bj-EVs and inhibits cellular EV uptake. Based on the proteomic results and the in vitro interaction assays using macrophages and muscle cells, we propose that Bj-EVs may be involved not only in venom production and processing but also in host immune modulation and long-term effects of envenoming.

From Snake Venoms to Therapeutics: A Focus on Natriuretic Peptides

Snake venom is a cocktail of multifunctional biomolecules that has evolved with the purpose of capturing prey and for defense. These biomolecules are classified into different classes based on their functions. They include three-finger toxins, natriuretic peptides, phospholipases and metalloproteinases. The focus for this review is on the natriuretic peptide (NP), which is an active component that can be isolated from the venoms of vipers and mambas. In these venoms, NPs contribute to the lowering of blood pressure, causing a rapid loss of consciousness in the prey such that its mobility is reduced, paralyzing the prey, and often death follows. Over the past 30 years since the discovery of the first NP in the venom of the green mamba, venom NPs have shown potential in the development of drug therapy for heart failure. Venom NPs have long half-lives, different pharmacological profiles, and may also possess different functions in comparison to the mammalian NPs. Understanding their mechanisms of action provides the strategies needed to develop new NPs for treatment of heart failure. This review summarizes the venom NPs that have been identified over the years and how they can be useful in drug development.

Biochemistry and toxicology of proteins and peptides purified from the venom of Vipera berus berus

The isolation and characterization of individual snake venom components is important for a deeper understanding of the pathophysiology of envenomation and for improving the therapeutic procedures of patients. It also opens possibilities for the discovery of novel toxins that might be useful as tools for understanding cellular and molecular processes. The variable venom composition, toxicological and immunological properties of the common vipers (Vipera berus berus) have been reviewed. The combination of venom gland transcriptomics, bottom-up and top-down proteomics enabled comparison of common viper venom proteomes from multiple individuals. V. b. berus venom contains proteins and peptides belonging to 10–15 toxin families: snake venom metalloproteinase, phospholipases A₂ (PLA₂), snake venom serine proteinase, aspartic protease, L-amino acid oxidase (LAAO), hyaluronidase, 5′-nucleotidase, glutaminyl-peptide cyclotransferase, disintegrin, C-type lectin (snaclec), nerve growth factor, Kunitz type serine protease inhibitor, snake venom vascular endothelial growth factor, cysteine-rich secretory protein, bradykinin potentiating peptide, natriuretic peptides. PLA₂ and LAAO from V. b. berus venom produce more pronounced cytotoxic effects in cancer cells than normal cells, via induction of apoptosis, cell cycle arrest and suppression of proliferation. Proteomic data of V. b. berus venoms from different parts of Russia and Slovakian Republic have been compared with analogous data for Vipera nikolskii venom. Proteomic studies demonstrated quantitative differences in the composition of V. b. berus venom from different geographical regions. Differences in the venom composition of V. berus were mainly driven by the age, sex, habitat and diet of the snakes. The venom variability of V. berus results in a loss of antivenom efficacy against snakebites. The effectiveness of antibodies is discussed. This review presents an overview with a special focus on different toxins that have been isolated and characterized from the venoms of V. b. berus. Their main biochemical properties and toxic actions are described.

•

Vipera berus berus venom composition is variable among different populations.

•

Venom contains about 15 protein/peptide families.

•

It disturbs blood coagulation inducing pro- or anticoagulant effects.

•

Venom contains different types of blood factor X activators.

•

PLA₂ and L-amino acid oxidase produce cytotoxic effects in cancer cells.

Proteomic insight into the venom composition of the largest European rear-fanged snake, Malpolon monspessulanus monspessulanus

Snake envenomations constitute a worldwide neglected tropical disease, with the vast majority of lethal bites inflicted by front-fanged snakes from the viperid and elapid groups. Rear-fanged snakes (colubrids) were often considered harmless and as a result, are much less studied, but several documented deaths have suggested potent venom in this group as well. The largest European snake (Malpolon monspessulanus monspessulanus), known as the “Montpellier snake”, is such a rear-fanged snake that belongs to the Lamprophiidae family. Its venom remains largely unknown but cases of envenomation with neurological symptoms have been reported. Here, we provide the first insights into the composition of its venom using mass spectrometry methods. First, liquid chromatography coupled mass spectrometry analysis of the manually collected venom samples reveals a complex profile, with the majority of masses encompassing the range 500–3000 Da, 4000–8000 Da, and 10 000–30 000 Da. Next, shotgun proteomics allowed the identification of a total of 42 different known families of proteins, including snake venom metalloproteinases, peptidase M1, and cysteine-rich secretory proteins, as the most prominent. Interestingly, three-finger toxins were not detected, suggesting that neurotoxicity may occur via other, yet to be determined, toxin types. Overall, our results provide the basis for a better understanding of the effects of a peculiar snake venom on human symptomatology, but also on the main prey consumed by this species.

•

We investigate the venom composition of the largest European venomous snake.

•

LC-MS analysis of the crude manually collected venom revealed a complex profile.

•

Shotgun proteomic analysis identified 42 different protein families.

•

Major components include SVMPs, consistent with the clinical features.

Differences between Two Groups of Burmese Vipers (Viperidae: Azemiops) in the Proteomic Profiles, Immunoreactivity and Biochemical Functions of Their Venoms

Two recently revised Azemiops snakes with apparent differences in their external appearances and skeletal morphologies but unclear genetic boundaries have been proposed. Some researchers have refrained from using the newly proposed taxonomy because these two “species” might be two clades corresponding to different geographical populations of Azemiops feae. To improve the understanding of the kinship of these two Burmese viper groups, more of their characteristics should be explored in depth. We performed a comparative analysis of the proteomic profiles and biochemical activities of snake venoms from these two groups (Sichuan A. feae and Zhejiang A. feae) and evaluated the immunorecognition capacity of commercial antivenoms toward them. Eight protein families were identified in venoms from these two groups, while phospholipase B was only detected in venom from Sichuan A. feae. These protein families displayed varying degrees of differences in relative abundance between venoms, and phospholipase A₂ (Sichuan A. feae: 57.15%; Zhejiang A. feae: 65.94%) was the predominated component. Gloydius brevicaudus antivenom exhibited the strongest capacity to immunologically recognize these two venoms, but this was mainly limited to components with high molecular masses, some of which differed between venoms. Additionally, Zhejiang A. feae venom was more toxic than Sichuan A. feae venom, and the venoms expressed remarkable differences in enzymatic activities, probably resulting from the variation in the relative abundance of specific protein families. Our findings unveil differences between the two Burmese viper groups in terms of proteomic profiles, immunoreactivity, and the biochemical functions of their venoms. This information will facilitate the management of snakebites caused by these snakes.

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.

The Target Selects the Toxin: Specific Amino Acids in Snake-Prey Nicotinic Acetylcholine Receptors That Are Selectively Bound by King Cobra Venoms

Snake venom is an adaptive ecological trait that has evolved primarily as a form of prey subjugation. Thus, the selection pressure for toxin diversification is exerted by the prey’s physiological targets, with this pressure being particularly acute for specialist feeders, such as the King Cobra species, all of which are snake-prey specialists. However, while extensive research has been undertaken to elucidate key amino acids that guide toxin structure–activity relationships, reciprocal investigations into the specific sites guiding prey-lineage selective effects have been lacking. This has largely been due to the lack of assay systems amenable to systematic amino acid replacements of targeted proteins in the prey’s physiological pathways. To fill this knowledge gap, we used a recently described approach based upon mimotope peptides corresponding to the orthosteric site of nicotinic acetylcholine receptor alpha-1 subunits, a major binding site for snake venom neurotoxins that cause flaccid paralysis. We investigated the venoms of four different types of King Cobra (Cambodian, Javan, Malaysian, and Thai). This approach allowed for the determination of the key amino acid positions in King Cobra snake prey that are selectively bound by the toxins, whereby replacing these amino acids in the snake-prey orthosteric site with those from lizards or rats resulted in a significantly lower level of binding by the venoms, while conversely replacing the lizard or rat amino acids with those from the snake at that position increased the binding. By doing such, we identified three negatively charged amino acids in the snake orthosteric site that are strongly bound by the positively charged neurotoxic three-finger toxins found in King Cobra venom. This study, thus, sheds light on the selection pressures exerted by a specialist prey item for the evolution of lineage-selective toxins.

Profilings of subproteomes of lectin-binding proteins of nine Bothrops venoms reveal variability driven by different glycan types

Snake venom proteomes have long been investigated to explore a multitude of biologically active components that are used for prey capture and defense, and are involved in the pathological effects observed upon mammalian envenomation. Glycosylation is a major protein post-translational modification in venoms and contributes to the diversification of proteomes. We have shown that Bothrops venoms are markedly defined by their content of glycoproteins, and that most N-glycan structures of eight Bothrops venoms contain sialic acid, while bisected N-acetylglucosamine was identified in Bothrops cotiara venom. To further investigate the mechanisms involved in the generation of different venoms by related snakes, here the glycoproteomes of nine Bothrops venoms (Bothrops atrox, B. cotiara, Bothrops erythromelas, Bothrops fonsecai, B. insularis, Bothrops jararaca, Bothrops jararacussu, Bothrops moojeni and Bothrops neuwiedi) were comparatively analyzed by enrichment with three lectins of different specificities, recognizing bisecting N-acetylglucosamine- and sialic acid-containing glycoproteins, and mass spectrometry. The lectin capture strategy generated venom fractions enriched with several glycoproteins, including metalloprotease, serine protease, and L- amino acid oxidase, in addition to various types of low abundant enzymes. The different contents of lectin-enriched proteins underscore novel aspects of the variability of the glycoprotein subproteomes of Bothrops venoms and point to the role of distinct types of glycan chains in generating different venoms by closely related snake species.

Snake venom gene expression is coordinated by novel regulatory architecture and the integration of multiple co-opted vertebrate pathways

Understanding how regulatory mechanisms evolve is critical for understanding the processes that give rise to novel phenotypes. Snake venom systems represent a valuable and tractable model for testing hypotheses related to the evolution of novel regulatory networks, yet the regulatory mechanisms underlying venom production remain poorly understood. Here, we use functional genomics approaches to investigate venom regulatory architecture in the prairie rattlesnake and identify cis-regulatory sequences (enhancers and promoters), trans-regulatory transcription factors, and integrated signaling cascades involved in the regulation of snake venom genes. We find evidence that two conserved vertebrate pathways, the extracellular signal-regulated kinase and unfolded protein response pathways, were co-opted to regulate snake venom. In one large venom gene family (snake venom serine proteases), this co-option was likely facilitated by the activity of transposable elements. Patterns of snake venom gene enhancer conservation, in some cases spanning 50 million yr of lineage divergence, highlight early origins and subsequent lineage-specific adaptations that have accompanied the evolution of venom regulatory architecture. We also identify features of chromatin structure involved in venom regulation, including topologically associated domains and CTCF loops that underscore the potential importance of novel chromatin structure to coevolve when duplicated genes evolve new regulatory control. Our findings provide a model for understanding how novel regulatory systems may evolve through a combination of genomic processes, including tandem duplication of genes and regulatory sequences, cis-regulatory sequence seeding by transposable elements, and diverse transcriptional regulatory proteins controlled by a co-opted regulatory cascade.

Transcriptomic Insights into the Diversity and Evolution of Myxozoa (Cnidaria, Endocnidozoa) Toxin-like Proteins

Myxozoa is a speciose group of endoparasitic cnidarians that can cause severe ecological and economic effects. Their cnidarian affinity is affirmed by genetic relatedness and the presence of nematocysts, historically called “polar capsules”. Previous studies have revealed the presence of toxin-like proteins in myxozoans; however, the diversity and evolution of venom in Myxozoa are not fully understood. Here, we performed a comparative analysis using the newly sequenced transcriptomes of five Myxobolidae species as well as some public datasets. Toxin mining revealed that myxozoans have lost most of their toxin families, while most species retained Kunitz, M12B, and CRISP, which may play a role in endoparasitism. The venom composition of Endocnidozoa (Myxozoa + Polypodium) differs from that of free-living cnidarians and may be influenced by ecological and environmental factors. Phylogenetic analyses showed that toxin families of myxozoans and free-living cnidarians were clustered into different clades. Selection analyses showed that purifying selection was the dominant evolutionary pressure in toxins, while they were still influenced by episodic adaptive selection. This suggests that the potency or specificity of a particular toxin or species might increase. Overall, our findings provide a more comprehensive framework for understanding the diversity and evolution of Myxozoa venoms.

Nanobiotechnology with Therapeutically Relevant Macromolecules from Animal Venoms: Venoms, Toxins, and Antimicrobial Peptides

Some diseases of uncontrolled proliferation such as cancer, as well as infectious diseases, are the main cause of death in the world, and their causative agents have rapidly developed resistance to the various existing treatments, making them even more dangerous. Thereby, the discovery of new therapeutic agents is a challenge promoted by the World Health Organization (WHO). Biomacromolecules, isolated or synthesized from a natural template, have therapeutic properties which have not yet been fully studied, and represent an unexplored potential in the search for new drugs. These substances, starting from conglomerates of proteins and other substances such as animal venoms, or from minor substances such as bioactive peptides, help fight diseases or counteract harmful effects. The high effectiveness of these biomacromolecules makes them promising substances for obtaining new drugs; however, their low bioavailability or stability in biological systems is a challenge to be overcome in the coming years with the help of nanotechnology. The objective of this review article is to describe the relationship between the structure and function of biomacromolecules of animal origin that have applications already described using nanotechnology and targeted delivery.

The rise of genomics in snake venom research: recent advances and future perspectives

Snake venoms represent a danger to human health, but also a gold mine of bioactive proteins that can be harnessed for drug discovery purposes. The evolution of snakes and their venom has been studied for decades, particularly via traditional morphological and basic genetic methods alongside venom proteomics. However, while the field of genomics has matured rapidly over the past 2 decades, owing to the development of next-generation sequencing technologies, snake genomics remains in its infancy. Here, we provide an overview of the state of the art in snake genomics and discuss its potential implications for studying venom evolution and toxinology. On the basis of current knowledge, gene duplication and positive selection are key mechanisms in the neofunctionalization of snake venom proteins. This makes snake venoms important evolutionary drivers that explain the remarkable venom diversification and adaptive variation observed in these reptiles. Gene duplication and neofunctionalization have also generated a large number of repeat sequences in snake genomes that pose a significant challenge to DNA sequencing, resulting in the need for substantial computational resources and longer sequencing read length for high-quality genome assembly. Fortunately, owing to constantly improving sequencing technologies and computational tools, we are now able to explore the molecular mechanisms of snake venom evolution in unprecedented detail. Such novel insights have the potential to affect the design and development of antivenoms and possibly other drugs, as well as provide new fundamental knowledge on snake biology and evolution.

Venom production and secretion in reptiles

The venom glands of reptiles, particularly those of front-fanged advanced snakes, must satisfy conflicting biological demands: rapid synthesis of potentially labile and highly toxic proteins, storage in the gland lumen for long periods, stabilization of the stored secretions, immediate activation of toxins upon deployment and protection of the animal from the toxic effects of its own venom. This dynamic system could serve as a model for the study of a variety of different phenomena involving exocrine gland activation, protein synthesis, stabilization of protein products and secretory mechanisms. However, these studies have been hampered by a lack of a long-term model that can be propagated in the lab (as opposed to whole-animal studies). Numerous attempts have been made to extend the lifetime of venom gland secretory cells, but only recently has an organoid model been shown to have the requisite qualities of recapitulation of the native system, self-propagation and long-term viability (>1 year). A tractable model is now available for myriad cell- and molecular-level studies of venom glands, protein synthesis and secretion. However, venom glands of reptiles are not identical, and many differ very extensively in overall architecture, microanatomy and protein products produced. This Review summarizes the similarities among and differences between venom glands of helodermatid lizards and of rear-fanged and front-fanged snakes, highlighting those areas that are well understood and identifying areas where future studies can fill in significant gaps in knowledge of these ancient, yet fascinating systems.

Venomics Reveals a Non-Compartmentalised Venom Gland in the Early Diverged Vermivorous Conus distans

The defensive use of cone snail venom is hypothesised to have first arisen in ancestral worm-hunting snails and later repurposed in a compartmentalised venom duct to facilitate the dietary shift to molluscivory and piscivory. Consistent with its placement in a basal lineage, we demonstrate that the C. distans venom gland lacked distinct compartmentalisation. Transcriptomics revealed C. distans expressed a wide range of structural classes, with inhibitory cysteine knot (ICK)-containing peptides dominating. To better understand the evolution of the venom gland compartmentalisation, we compared C. distans to C. planorbis, the earliest diverging species from which a defence-evoked venom has been obtained, and fish-hunting C. geographus from the Gastridium subgenus that injects distinct defensive and predatory venoms. These comparisons support the hypothesis that venom gland compartmentalisation arose in worm-hunting species and enabled repurposing of venom peptides to facilitate the dietary shift from vermivory to molluscivory and piscivory in more recently diverged cone snail lineages.

Venom of Viperidae: A Perspective of its Antibacterial and Antitumor Potential

:

The emergence of multi-drug resistant bacteria and limitations on cancer treatment represent two important challenges in modern medicine. Biological compounds have been explored with a particular focus on venoms. Although they can be lethal or cause considerable damage to humans, venom is also a source rich in components with high therapeutic potential.

:

Viperidae family is one of the most emblematic venomous snake families and several studies highlighted the antibacterial and antitumor potential of viper toxins. According to the literature, these activities are mainly associated to five protein families - svLAAO, Disintegrins, PLA2, SVMPs and C-type lectins- that act through different mechanisms leading to the inhibition of the growth of bacteria, as well as, cytotoxic effects and inhibition of metastasis process. In this review, we provide an overview of the venom toxins produced by species belonging to the Viperidae family, exploring their roles during the envenoming and their pharmacological properties, in order to demonstrate its antibacterial and antitumor potential.

Correlations between rotational thromboelastometry (ROTEM) and standard coagulation tests following viper snakebites

BACKGROUND

A high prevalence of venom-induced consumption coagulopathy has been reported in individuals with viper snakebites. Rotational thromboelastometry (ROTEM) is a rapid technique that could be advantageous in assessing and monitoring coagulation disorders.

PURPOSE

To explore correlations between ROTEM and standard coagulation tests.

PATIENTS AND METHODS

This prospective observational study was performed among 41 patients with viper envenomation admitted to the Vietnam Poison Control Center from April 2016 to October 2017. Standard coagulation measurements [platelet count, prothrombin time (PT), activated partial thromboplastin time (aPTT), and fibrinogen level] and ROTEM indicators [clotting time (CT), amplitude (at set time: 5 and 10 minutes), clot information time (CFT) and maximum clot firmness (MCF) for extrinsic (EXTEM), intrinsic (INTEM), and fibrin based (FIBTEM) ROTEM] were obtained.

RESULTS

For INTEM, EXTEM, the FIBTEM, proportions of patients with prolonged CT were 34.1%, 63.4%, and 61.0% respectively and the proportions of patients with decreased MCF were 62.2%, 62.2%, and 35.5%, respectively. Moderate correlations were observed between PT and EXTEM CT (r = 0.627), aPTT and INTEM CT (r = 0.626), fibrinogen and FIBTEM MCF (r = 0.723), and platelet count and EXTEM MCF (0.60).

CONCLUSION

ROTEM indicated a hypocoagulation state in patients with viper snakebite and was moderately correlated with standard coagulation parameters.

Dynamic genetic differentiation drives the widespread structural and functional convergent evolution of snake venom proteinaceous toxins

Background

The explosive radiation and diversification of the advanced snakes (superfamily Colubroidea) was associated with changes in all aspects of the shared venom system. Morphological changes included the partitioning of the mixed ancestral glands into two discrete glands devoted for production of venom or mucous respectively, as well as changes in the location, size and structural elements of the venom-delivering teeth. Evidence also exists for homology among venom gland toxins expressed across the advanced snakes. However, despite the evolutionary novelty of snake venoms, in-depth toxin molecular evolutionary history reconstructions have been mostly limited to those types present in only two front-fanged snake families, Elapidae and Viperidae. To have a broader understanding of toxins shared among extant snakes, here we first sequenced the transcriptomes of eight taxonomically diverse rear-fanged species and four key viperid species and analysed major toxin types shared across the advanced snakes.

Results

Transcriptomes were constructed for the following families and species: Colubridae - Helicops leopardinus, Heterodon nasicus, Rhabdophis subminiatus; Homalopsidae – Homalopsis buccata; Lamprophiidae - Malpolon monspessulanus, Psammophis schokari, Psammophis subtaeniatus, Rhamphiophis oxyrhynchus; and Viperidae – Bitis atropos, Pseudocerastes urarachnoides, Tropidolaeumus subannulatus, Vipera transcaucasiana. These sequences were combined with those from available databases of other species in order to facilitate a robust reconstruction of the molecular evolutionary history of the key toxin classes present in the venom of the last common ancestor of the advanced snakes, and thus present across the full diversity of colubroid snake venoms. In addition to differential rates of evolution in toxin classes between the snake lineages, these analyses revealed multiple instances of previously unknown instances of structural and functional convergences. Structural convergences included: the evolution of new cysteines to form heteromeric complexes, such as within kunitz peptides (the beta-bungarotoxin trait evolving on at least two occasions) and within SVMP enzymes (the P-IIId trait evolving on at least three occasions); and the C-terminal tail evolving on two separate occasions within the C-type natriuretic peptides, to create structural and functional analogues of the ANP/BNP tailed condition. Also shown was that the de novo evolution of new post-translationally liberated toxin families within the natriuretic peptide gene propeptide region occurred on at least five occasions, with novel functions ranging from induction of hypotension to post-synaptic neurotoxicity. Functional convergences included the following: multiple occasions of SVMP neofunctionalised in procoagulant venoms into activators of the clotting factors prothrombin and Factor X; multiple instances in procoagulant venoms where kunitz peptides were neofunctionalised into inhibitors of the clot destroying enzyme plasmin, thereby prolonging the half-life of the clots formed by the clotting activating enzymatic toxins; and multiple occasions of kunitz peptides neofunctionalised into neurotoxins acting on presynaptic targets, including twice just within Bungarus venoms.

Conclusions

We found novel convergences in both structural and functional evolution of snake toxins. These results provide a detailed roadmap for future work to elucidate predator–prey evolutionary arms races, ascertain differential clinical pathologies, as well as documenting rich biodiscovery resources for lead compounds in the drug design and discovery pipeline.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01208-9.

Revisiting the potential of South American rattlesnake Crotalus durissus terrificus toxins as therapeutic, theranostic and/or biotechnological agents

The potential biotechnological and biomedical applications of the animal venom components are widely recognized. Indeed, many components have been used either as drugs or as templates/prototypes for the development of innovative pharmaceutical drugs, among which many are still used for the treatment of human diseases. A specific South American rattlesnake, named Crotalus durissus terrificus, shows a venom composition relatively simpler compared to any viper or other snake species belonging to the Crotalus genus, although presenting a set of toxins with high potential for the treatment of several still unmet human therapeutic needs, as reviewed in this work. In addition to the main toxin named crotoxin, which is under clinical trials studies for antitumoral therapy and which has also anti-inflammatory and immunosuppressive activities, other toxins from the C. d. terrificus venom are also being studied, aiming for a wide variety of therapeutic applications, including as antinociceptive, anti-inflammatory, antimicrobial, antifungal, antitumoral or antiparasitic agent, or as modulator of animal metabolism, fibrin sealant (fibrin glue), gene carrier or theranostic agent. Among these rattlesnake toxins, the most relevant, considering the potential clinical applications, are crotamine, crotalphine and gyroxin. In this narrative revision, we propose to organize and present briefly the updates in the accumulated knowledge on potential therapeutic applications of toxins collectively found exclusively in the venom of this specific South American rattlesnake, with the objective of contributing to increase the chances of success in the discovery of drugs based on toxins.

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.

A non-lethal method for studying scorpion venom gland transcriptomes, with a review of potentially suitable taxa to which it can be applied

Scorpion venoms are mixtures of proteins, peptides and small molecular compounds with high specificity for ion channels and are therefore considered to be promising candidates in the venoms-to-drugs pipeline. Transcriptomes are important tools for studying the composition and expression of scorpion venom. Unfortunately, studying the venom gland transcriptome traditionally requires sacrificing the animal and therefore is always a single snapshot in time. This paper describes a new way of generating a scorpion venom gland transcriptome without sacrificing the animal, thereby allowing the study of the transcriptome at various time points within a single individual. By comparing these venom-derived transcriptomes to the traditional whole-telson transcriptomes we show that the relative expression levels of the major toxin classes are similar. We further performed a multi-day extraction using our proposed method to show the possibility of doing a multiple time point transcriptome analysis. This allows for the study of patterns of toxin gene activation over time a single individual, and allows assessment of the effects of diet, season and other factors that are known or likely to influence intraindividual venom composition. We discuss the gland characteristics that may allow this method to be successful in scorpions and provide a review of other venomous taxa to which this method may potentially be successfully applied.

Venom proteomic analysis of medically important Nigerian viper Echis ocellatus and Bitis arietans snake species

Snakebite envenoming remains a neglected tropical disease which poses severe health hazard, especially for the rural inhabitants in Africa. In Nigeria, vipers are responsible for the highest number of deaths. Hydrophilic interaction liquid chromatography coupled with LC-MS/MS was used to analyze the crude venoms of Echis ocellatus (Carpet viper) and Bitis arietans (Puff adder) in order to understand their venom proteomic identities. Results obtained revealed that gel-free proteomic analysis of the crude venoms led to the identification of 85 and 79 proteins, respectively. Seventy-eight (78) proteins were common between the two snake species with a 91.8% similarity score. The identified proteins belong to 18 protein families in E. ocellatus and 14 protein families in B. arietans. Serine proteases (22.31%) and metalloproteinases (21.06%) were the dominant proteins in the venom of B. arietans; while metalloproteinases (34.84%), phospholipase A₂s (21.19%) and serine proteases (15.50%) represent the major toxins in the E. ocellatus venom. Other protein families such as three-finger toxins and cysteine-rich venom proteins were detected in low proportions. This study provides an insight into the venom proteomic analysis of the two Nigerian viper species, which could be useful in identifying the toxin families to be neutralized in case of envenomation.

•

Venom proteomic of Nigeria's most medically important snakes is presented.

•

SVMP, SVSP and PLA₂ were the major toxin families in E. ocellatus and B. arietans.

•

The venom proteomes of these vipers displayed 91.8% similarity in composition.

Rapid increase in snake dietary diversity and complexity following the end-Cretaceous mass extinction

The Cenozoic marked a period of dramatic ecological opportunity in Earth history due to the extinction of non-avian dinosaurs as well as to long-term physiographic changes that created new biogeographic theaters and new habitats. Snakes underwent massive ecological diversification during this period, repeatedly evolving novel dietary adaptations and prey preferences. The evolutionary tempo and mode of these trophic ecological changes remain virtually unknown, especially compared with co-radiating lineages of birds and mammals that are simultaneously predators and prey of snakes. Here, we assemble a dataset on snake diets (34,060 observations on the diets of 882 species) to investigate the history and dynamics of the multidimensional trophic niche during the global radiation of snakes. Our results show that per-lineage dietary niche breadths remained remarkably constant even as snakes diversified to occupy disparate outposts of dietary ecospace. Rapid increases in dietary diversity and complexity occurred in the early Cenozoic, and the overall rate of ecospace expansion has slowed through time, suggesting a potential response to ecological opportunity in the wake of the end-Cretaceous mass extinction. Explosive bursts of trophic innovation followed colonization of the Nearctic and Neotropical realms by a group of snakes that today comprises a majority of living snake diversity. Our results indicate that repeated transformational shifts in dietary ecology are important drivers of adaptive radiation in snakes and provide a framework for analyzing and visualizing the evolution of complex ecological phenotypes on phylogenetic trees.

The Cenozoic marked a period of dramatic ecological opportunity in Earth history due to the extinction of non-avian dinosaurs and long-term physiographic changes. This phylogenetic natural history study offers new insights into the evolution of snake ecological diversity after the end-Cretaceous mass extinction, as they took advantage of these new opportunities.

Plicidentine and the repeated origins of snake venom fangs

Snake fangs are an iconic exemplar of a complex adaptation, but despite striking developmental and morphological similarities, they probably evolved independently in several lineages of venomous snakes. How snakes could, uniquely among vertebrates, repeatedly evolve their complex venom delivery apparatus is an intriguing question. Here we shed light on the repeated evolution of snake venom fangs using histology, high-resolution computed tomography (microCT) and biomechanical modelling. Our examination of venomous and non-venomous species reveals that most snakes have dentine infoldings at the bases of their teeth, known as plicidentine, and that in venomous species, one of these infoldings was repurposed to form a longitudinal groove for venom delivery. Like plicidentine, venom grooves originate from infoldings of the developing dental epithelium prior to the formation of the tooth hard tissues. Derivation of the venom groove from a large plicidentine fold that develops early in tooth ontogeny reveals how snake venom fangs could originate repeatedly through the co-option of a pre-existing dental feature even without close association to a venom duct. We also show that, contrary to previous assumptions, dentine infoldings do not improve compression or bending resistance of snake teeth during biting; plicidentine may instead have a role in tooth attachment.

Impact of Repetitive DNA Elements on Snake Genome Biology and Evolution

The distinctive biology and unique evolutionary features of snakes make them fascinating model systems to elucidate how genomes evolve and how variation at the genomic level is interlinked with phenotypic-level evolution. Similar to other eukaryotic genomes, large proportions of snake genomes contain repetitive DNA, including transposable elements (TEs) and satellite repeats. The importance of repetitive DNA and its structural and functional role in the snake genome, remain unclear. This review highlights the major types of repeats and their proportions in snake genomes, reflecting the high diversity and composition of snake repeats. We present snakes as an emerging and important model system for the study of repetitive DNA under the impact of sex and microchromosome evolution. We assemble evidence to show that certain repetitive elements in snakes are transcriptionally active and demonstrate highly dynamic lineage-specific patterns as repeat sequences. We hypothesize that particular TEs can trigger different genomic mechanisms that might contribute to driving adaptive evolution in snakes. Finally, we review emerging approaches that may be used to study the expression of repetitive elements in complex genomes, such as snakes. The specific aspects presented here will stimulate further discussion on the role of genomic repeats in shaping snake evolution.

A symphony of destruction: Dynamic differential fibrinogenolytic toxicity by rattlesnake (Crotalus and Sistrurus) venoms

What factors influence the evolution of a heavily selected functional trait in a diverse clade? This study adopts rattlesnakes as a model group to investigate the evolutionary history of venom coagulotoxicity in the wider context of phylogenetics, natural history, and biology. Venom-induced clotting of human plasma and fibrinogen was determined and mapped onto the rattlesnake phylogenetic tree to reconstruct the evolution of coagulotoxicity across the group. Our results indicate that venom phenotype is often independent of phylogenetic relationships in rattlesnakes, suggesting the importance of diet and/or other environmental variables in driving venom evolution. Moreover, the striking inter- and intraspecific variability in venom activity on human blood highlights the considerable variability faced by physicians treating envenomation. This study is the most comprehensive effort to date to describe and characterize the evolutionary and biological aspects of coagulotoxins in rattlesnake venom. Further research at finer taxonomic levels is recommended to elucidate patterns of variation within species and lineages.

Development of the venom delivery system in Elapidae snake species: Naja siamensis and Oxyuranus microlepidotus

Many advanced snakes possess a unique venom delivery system which they utilise to subdue prey and for defence. Despite extensive efforts, the evolutionary differences in this key system between advanced snake families remains enigmatic. The current study has investigated the development of the venom delivery system using two oviparous Elapidae models, Naja siamensis and Oxyuranus microlepidotus. The development stages of the embryos in both models were detailed using previously standardised characterisation. Variations in the days post-oviposition between these stages was observed, despite a continuous development trajectory. These differences also translated to the development of the venom delivery system.

Fang evolution in venomous snakes: Adaptation of 3D tooth shape to the biomechanical properties of their prey

Venomous snakes are among the world's most specialized predators. During feeding, they use fangs to penetrate the body tissues of their prey, but the success of this penetration depends on the shape of these highly specialized teeth. Here, we examined the evolution of fang shape in a wide range of snakes using 3D geometric morphometrics (3DGM) and cross-sectional tooth sharpness measurements. We investigated the relationship of these variables with six diet categories based on the prey's biomechanical properties, and tested for evolutionary convergence using two methods. Our results show that slender elongate fangs with sharp tips are used by snakes that target soft-skinned prey (e.g., mammals), whereas fangs become more robust and blunter as the target's skin becomes scaly (e.g., fish and reptiles) and eventually hard-shelled (e.g., crustaceans), both with and without correction for evolutionary allometry. Convergence in fang shape is present, indicating that fangs of snakes with the same diet are more similar than those of closely related species with different diets. Establishing the relationship between fang morphology and diet helps to explain how snakes became adapted to different lifestyles, while also providing a proxy to infer diet in lesser known species or extinct snakes from the fossil record.

Phylogenetically diverse diets favor more complex venoms in North American pitvipers

The role of natural selection in the evolution of trait complexity can be characterized by testing hypothesized links between complex forms and their functions across species. Predatory venoms are composed of multiple proteins that collectively function to incapacitate prey. Venom complexity fluctuates over evolutionary timescales, with apparent increases and decreases in complexity, and yet the causes of this variation are unclear. We tested alternative hypotheses linking venom complexity and ecological sources of selection from diet in the largest clade of front-fanged venomous snakes in North America: the rattlesnakes, copperheads, cantils, and cottonmouths. We generated independent transcriptomic and proteomic measures of venom complexity and collated several natural history studies to quantify dietary variation. We then constructed genome-scale phylogenies for these snakes for comparative analyses. Strikingly, prey phylogenetic diversity was more strongly correlated to venom complexity than was overall prey species diversity, specifically implicating prey species' divergence, rather than the number of lineages alone, in the evolution of complexity. Prey phylogenetic diversity further predicted transcriptomic complexity of three of the four largest gene families in viper venom, showing that complexity evolution is a concerted response among many independent gene families. We suggest that the phylogenetic diversity of prey measures functionally relevant divergence in the targets of venom, a claim supported by sequence diversity in the coagulation cascade targets of venom. Our results support the general concept that the diversity of species in an ecological community is more important than their overall number in determining evolutionary patterns in predator trait complexity.

An ancient, conserved gene regulatory network led to the rise of oral venom systems

Although oral venom systems are ecologically important characters, how they originated is still unclear. In this study, we show that oral venom systems likely originated from a gene regulatory network conserved across amniotes. This network, which we term the “metavenom network,” comprises over 3,000 housekeeping genes coexpressed with venom and play a role in protein folding and modification. Comparative transcriptomics revealed that the network is conserved between venom glands of snakes and salivary glands of mammals. This suggests that while these tissues have evolved different functions, they share a common regulatory core, that persisted since their common ancestor. We propose several evolutionary mechanisms that can utilize this common regulatory core to give rise to venomous animals from their nonvenomous ancestors.

Oral venom systems evolved multiple times in numerous vertebrates enabling the exploitation of unique predatory niches. Yet how and when they evolved remains poorly understood. Up to now, most research on venom evolution has focused strictly on the toxins. However, using toxins present in modern day animals to trace the origin of the venom system is difficult, since they tend to evolve rapidly, show complex patterns of expression, and were incorporated into the venom arsenal relatively recently. Here we focus on gene regulatory networks associated with the production of toxins in snakes, rather than the toxins themselves. We found that overall venom gland gene expression was surprisingly well conserved when compared to salivary glands of other amniotes. We characterized the “metavenom network,” a network of ∼3,000 nonsecreted housekeeping genes that are strongly coexpressed with the toxins, and are primarily involved in protein folding and modification. Conserved across amniotes, this network was coopted for venom evolution by exaptation of existing members and the recruitment of new toxin genes. For instance, starting from this common molecular foundation, Heloderma lizards, shrews, and solenodon, evolved venoms in parallel by overexpression of kallikreins, which were common in ancestral saliva and induce vasodilation when injected, causing circulatory shock. Derived venoms, such as those of snakes, incorporated novel toxins, though still rely on hypotension for prey immobilization. These similarities suggest repeated cooption of shared molecular machinery for the evolution of oral venom in mammals and reptiles, blurring the line between truly venomous animals and their ancestors.

Proteomic analysis of three medically important Nigerian Naja (Naja haje, Naja katiensis and Naja nigricollis) snake venoms

Proteomics technologies enable a comprehensive study of complex proteins and their functions. The venom proteomes of three medically important Nigerian Elapidae snakes Naja haje, Naja katiensis and Naja nigricollis was studied using HILIC coupled with LC-MS/MS analysis. Results revealed a total of 57, 55, and 46 proteins in the venoms of N. haje, N. katiensis, and N. nigricollis, respectively, with molecular mass ranging between 5 and 185  kDa. These snakes have 38 common proteins in addition to 3 uncommon proteins: actiflagelin, cathelicidin, and cystatin identified in their venoms. The identified proteins belonged to 14 protein families in N. haje and N. katiensis, and 12 protein families in N. nigricollis. Of the total venom proteins, 3FTx was the most abundant protein family, constituting 52% in N. haje and N. katiensis, and 41% in N. nigricollis, followed by PLA₂, constituting 37% in N. nigricollis, 26% in N. haje, and 24% in N. katiensis. Other protein families, including LAAO, CRISPs, VEGF, PLB, CVF, SVMP, SVH, AMP, PI, Globin, Actin, and C-type lectins, were also detected, although, at very low abundances. Quantification of the relative abundance of each protein revealed that alpha and beta fibrinogenase and PLA_2, which constituted 18-26% of the total proteome, were the most abundant. The 3 uncommon proteins have no known function in snake venom. However, actiflagelin activates sperm motility; cystatin inhibits angiogenesis, while cathelicidin exerts antimicrobial effects. The three Nigerian Naja genus proteomes displayed 70% similarity in composition, which suggests the possibility of formulating antivenom that may cross-neutralise the venoms of cobra species found in Nigeria. These data provide insights into clinically relevant peptides/proteins present in the venoms of these snakes. Data are available via ProteomeXchange with identifier PXD024627.

Comparative venomics and preclinical efficacy evaluation of a monospecific Hemachatus antivenom towards sub-Saharan Africa cobra venoms

Cobras are the most medically important elapid snakes in Africa. The African genera Naja and Hemachatus include snakes with neurotoxic and cytotoxic venoms, with shared biochemical, toxinological and antigenic characteristics. We have studied the antigenic cross-reactivity of four sub-Saharan Africa cobra venoms against an experimental monospecific Hemachatus haemachatus antivenom through comparative proteomics, preclinical assessment of neutralization, and third generation antivenomics. The venoms of H. haemachatus, N. annulifera, N. mossambica and N. nigricollis share an overall qualitative family toxin composition but depart in their proportions of three-finger toxin (3FTxs) classes, phospholipases A₂ (PLA₂s), snake venom metalloproteinases (SVMPs), and cysteine-rich secretory proteins (CRISPs). A monospecific anti-Hemachatus antivenom produced by Costa Rican Instituto Clodomiro Picado neutralized the lethal activity of the homologous and heterologous neuro/cytotoxic (H. haemachatus) and cyto/cardiotoxic (N. mossambica and N. nigricollis) venoms of the three spitting cobras sampled, while it was ineffective against the lethal and toxic activities of the neurotoxic venom of the non-spitting snouted cobra N. annulifera. The ability of the anti-Hemachatus-ICP antivenom to neutralize toxic (dermonecrotic and anticoagulant) and enzymatic (PLA₂) activities of spitting cobra venoms suggested a closer kinship of H. haemachatus and Naja subgenus Afrocobra spitting cobras than to Naja subgenus Uraeus neurotoxic taxa. These results were confirmed by third generation antivenomics. BIOLOGICAL SIGNIFICANCE: African Naja species represent the most widespread medically important elapid snakes across Africa. To gain deeper insight into the spectrum of medically relevant toxins, we compared the proteome of three spitting cobras (Hemachatus haemachatus, Naja mossambica and N. nigricollis) and one non-spitting cobra (N. annulifera). Three finger toxins and phospholipases A₂ are the two major protein families among the venoms analyzed. The development of antivenoms of broad species coverage is an urgent need in sub-Saharan Africa. An equine antivenom raised against H. haemachatus venom showed cross-reactivity with the venoms of H. haemachatus, N. mossambica and N. nigricollis, while having poor recognition of the venom of N. annulifera. This immunological information provides clues for the design of optimum venom mixtures for the preparation of broad spectrum antivenoms.

Coagulopathy after snake envenomation in South Korea

OBJECTIVE

This retrospective study investigated the nature and severity of venom-induced consumption coagulopathy (VICC) and determined the clotting factors involved in VICC in patients after envenomation by South Korea's snakes. Additionally, we studied the effectiveness of antivenom for the treatment of VICC after envenomation.

METHODS

Included patients were divided into three groups according to the severity of VICC (no VICC, partial VICC, and complete VICC). Data, including changes in coagulation parameters during hospitalization and clotting factors at presentation, were collected and analyzed.

RESULTS

One hundred nineteen patients who presented at our emergency department within 3 h after snake envenomation were included. VICC developed in 34 patients (27 patients with partial VICC and 7 patients with complete VICC). Two of 34 patients with VICC required blood transfusions. Five patients with complete VICC had an undetectable fibrinogen concentration at presentation. Three patients with complete VICC had an unmeasurable INR and aPTT within 24 h. The median times of the most extreme values were 10 h for INR, 12 h for aPTT, and 16 h for fibrinogen after presentation in the VICC group. The D-dimer concentration peaked at a median of 63.5 h after presentation. The activities of factors II and X were significantly reduced in the complete VICC group (factor II: 88 (84-99.3)% in the non-VICC group vs. 69 (49.5-83.5)% in the complete VICC group; factor X:94 (83-102) in the non-VICC group vs. 70 (66.5-79.8)% in the complete VICC group), while there was no difference in factor V activity at presentation. The time from bite to first antivenom administration did not correlate with the time course and most extreme concentrations for fibrinogen and D-dimer within the VICC groups.

DISCUSSION AND CONCLUSION

VICC occurs in approximately one-quarter of snakebite patients in South Korea; however, VICC itself does not appear to lead to clinical deterioration. Fibrinogen is an early diagnostic maker for complete VICC. Clotting factors II and X are involved in VICC. Future investigations should explore the mechanism of VICC from Korean snakebites and the effect of antivenom on VICC.