From genome to "venome": molecular origin and evolution of the snake venom proteome inferred from phylogenetic analysis of toxin sequences and related body proteins

The Clot Thickens: Differential Coagulotoxic and Cardiotoxic Activities of Anguimorpha Lizard Venoms

Despite their evolutionary novelty, lizard venoms are much less studied in comparison to the intense research on snake venoms. While the venoms of helodermatid lizards have long been assumed to be for defensive purposes, there is increasing evidence of toxic activities more useful for predation than defence (such as paralytic neurotoxicity). This study aimed to ascertain the effects of Heloderma, Lanthanotus, and Varanus lizard venoms on the coagulation and cardiovascular systems. Anticoagulant toxicity was demonstrated for the Varanus species studied, with the venoms prolonging clotting times in human and bird plasma due to the destructive cleavage of fibrinogen. In contrast, thromboelastographic analyses on human and bird plasmas in this study demonstrated a procoagulant bioactivity for Heloderma venoms. A previous study on Heloderma venom using factor-depleted plasmas as a proxy model suggested a procoagulant factor was present that activated either Factor XI or Factor XII, but could not ascertain the precise target. Our activation studies using purified zymogens confirmed FXII activation. Comparisons of neonate and adult H. exasperatum, revealed the neonates to be more potent in the ability to activate FXII, being more similar to the venom of the smaller species H. suspectum than the adult H. exasperatum. This suggests potent FXII activation a basal trait in the genus, present in the small bodied last common ancestor. This also indicates an ontogenetic difference in prey preferences in the larger Heloderma species paralleing the change in venom biochemistry. In addition, as birds lack Factor XII, the ability to clot avian plasma suggested an additional procoagulant site of action, which was revealed to be the activation of Factor VII, with H. horridum being the most potent. This study also examined the effects upon the cardiovascular system, including the liberation of kinins from kininogen, which contributes to hypotension induction. This form of toxicity was previously described for Heloderma venoms, and was revealed in this study was to also be a pathophysiological effect of Lanthanotus and Varanus venoms. This suggests that this toxic activity was present in the venom of the last common ancestor of the anguimorph lizards, which is consistent with kallikrein enzymes being a shared toxin trait. This study therefore uncovered novel actions of anguimorph lizard venoms, not only contributing to the evolutionary biology body of knowledge but also revealing novel activities to mine for drug design lead compounds.

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.

The evolutionary novelty of insect defensins: from bacterial killing to toxin neutralization

Insect host defense comprises two complementary dimensions, microbial killing-mediated resistance and microbial toxin neutralization-mediated resilience, both jointly providing protection against pathogen infections. Insect defensins are a class of effectors of innate immunity primarily responsible for resistance to Gram-positive bacteria. Here, we report a newly originated gene from an ancestral defensin via genetic deletion following gene duplication in Drosophila virilis, which confers an enhanced resilience to Gram-positive bacterial infection. This gene encodes an 18-mer arginine-rich peptide (termed DvirARP) with differences from its parent gene in its pattern of expression, structure and function. DvirARP specifically expresses in D. virilis female adults with a constitutive manner. It adopts a novel fold with a 310 helix and a two CXC motif-containing loop stabilized by two disulfide bridges. DvirARP exhibits no activity on the majority of microorganisms tested and only a weak activity against two Gram-positive bacteria. DvirARP knockout flies are viable and have no obvious defect in reproductivity but they are more susceptible to the DvirARP-resistant Staphylococcus aureus infection than the wild type files, which can be attributable to its ability in neutralization of the S. aureus secreted toxins. Phylogenetic distribution analysis reveals that DvirARP is restrictedly present in the Drosophila subgenus, but independent deletion variations also occur in defensins from the Sophophora subgenus, in support of the evolvability of this class of immune effectors. Our work illustrates for the first time how a duplicate resistance-mediated gene evolves an ability to increase the resilience of a subset of Drosophila species against bacterial infection.

Indian green pit vipers: A lesser-known snake group of north-east India

Green pit vipers are one of the most widely distributed group of venomous snakes in south-east Asia. In Indian, green pit vipers are found in the Northern and North-eastern states spreading across eastern and central India and one of the lesser studied venoms. High morphological similarity among them has been a long-established challenge for species identification, however, a total of six species of Indian green pit viper belonging to genus Trimeresurus, Popeia and Viridovipera has been reported from North-east India. Biochemical and biological studies have revealed that venom exhibits substantial variation in protein expression level along with functional variability. The symptoms of envenomation are painful swelling at bite site, bleeding, necrosis along with systemic toxicity such as prolonged coagulopathy. Clinical data of green pit viper envenomated patients from Demow community health centre, Assam advocated against the use of Indian polyvalent antivenom pressing the need for a suitable antivenom for the treatment of green pit viper envenomation. To design effective and specific antivenom for green pit vipers, unveiling the proteome profile of these snakes is needed. In this study, a comparative venomic of green pit vipers of Northern and North-eastern India, their clinical manifestation as well as treatment protocol has been reviewed.

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.

Venomous Noodles: The Evolution of Toxins in Nemertea through Positive Selection and Gene Duplication

Some, probably most and perhaps all, members of the phylum Nemertea are poisonous, documented so far from marine and benthic specimens. Although the toxicity of these animals has been long known, systematic studies on the characterization of toxins, mechanisms of toxicity, and toxin evolution for this group are scarce. Here, we present the first investigation of the molecular evolution of toxins in Nemertea. Using a proteo-transcriptomic approach, we described toxins in the body and poisonous mucus of the pilidiophoran Lineus sanguineus and the hoplonemertean Nemertopsis pamelaroeae. Using these new and publicly available transcriptomes, we investigated the molecular evolution of six selected toxin gene families. In addition, we also characterized in silico the toxin genes found in the interstitial hoplonemertean, Ototyphlonemertes erneba, a meiofaunal taxa. We successfully identified over 200 toxin transcripts in each of these species. Evidence of positive selection and gene duplication was observed in all investigated toxin genes. We hypothesized that the increased rates of gene duplications observed for Pilidiophora could be involved with the expansion of toxin genes. Studies concerning the natural history of Nemertea are still needed to understand the evolution of their toxins. Nevertheless, our results show evolutionary mechanisms similar to other venomous groups.

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

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

Pulmonary involvement from animal toxins: the cellular mechanisms

Venomous animals and their venom have always been of human interest because, despite species differences, coevolution has made them capable of targeting key physiological components of our bodies. Respiratory failure from lung injury is one of the serious consequences of envenomation, and the underlying mechanisms are rarely discussed. This review aims to demonstrate how toxins affect the pulmonary system through various biological pathways. Herein, we propose the common underlying cellular mechanisms of toxin-induced lung injury: interference with normal cell function and integrity, disruption of normal vascular function, and provocation of excessive inflammation. Viperid snakebites are the leading cause of envenomation-induced lung injury, followed by other terrestrial venomous animals such as scorpions, spiders, and centipedes. Marine species, particularly jellyfish, can also inflict such injury. Common pulmonary manifestations include pulmonary edema, pulmonary hemorrhage, and exudative infiltration. Severe envenomation can result in acute respiratory distress syndrome. Pulmonary involvement suggests severe envenomation, thus recognizing these mechanisms and manifestations can aid physicians in providing appropriate treatment.

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.

A novel broad spectrum venom metalloproteinase autoinhibitor in the rattlesnake Crotalus atrox evolved via a shift in paralog function

The complexity of snake venom composition reflects adaptation to the diversity of prey and may be driven at times by a coevolutionary arms race between snakes and venom-resistant prey. However, many snakes are also resistant to their own venom due to serum-borne inhibitors of venom toxins, which raises the question of how snake autoinhibitors maintain their efficacy as venom proteins evolve. To investigate this potential three-way arms race among venom, prey, and autoinhibitors, we have identified and traced the evolutionary origin of serum inhibitors of snake venom metalloproteinases (SVMPs) in the Western Diamondback rattlesnake Crotalus atrox which possesses the largest known battery of SVMP genes among crotalids examined. We found that C. atrox expresses five members of a Fetuin A-related metalloproteinase inhibitor family but that one family member, FETUA-3, is the major SVMP inhibitor that binds to approximately 20 different C. atrox SVMPs and inhibits activities of all three SVMP classes. We show that the fetua-3 gene arose deep within crotalid evolution before the origin of New World species but, surprisingly, fetua-3 belongs to a different paralog group than previously identified SVMP inhibitors in Asian and South American crotalids. Conversely, the C. atrox FETUA-2 ortholog of previously characterized crotalid SVMP inhibitors shows limited activity against C. atrox SVMPs. These results reveal that there has been a functional evolutionary shift in the major SVMP inhibitor in the C. atrox lineage as the SVMP family expanded and diversified in the Crotalus lineage. This broad-spectrum inhibitor may be of potential therapeutic interest.

CAP superfamily proteins from venomous animals: Who we are and what to do?

Cysteine-rich secretory proteins (CRISPs), antigen 5 (Ag5), and pathogenesis-related (PR-1) superfamily proteins (CAP superfamily proteins) are found in diverse species across the bacterial, fungal, plant, mammalian, and venomous animal kingdoms. Notably, CAP proteins are found in a remarkable range of species across the venomous animal kingdom and are present almost ubiquitously in venoms, even when venoms are produced in very small quantities. Meanwhile, in comparison to mammals, venomous animals are underappreciated and easy to ignore. Overwhelming evidence suggests that CAP proteins derived from venomous animals exhibit diverse activities, including ion channel, inflammatory, proteolysis, and immune regulatory activities. To understand the potential biological functions of CAP proteins in venom more effectively, we need to examine the significance of the evolution of venomous animals in the animal kingdom, for their survival. In this article, we will review the current status of research on CAP proteins in venomous animals, including their isolation, characterization, known biological activities, and sequence alignments. We will also discuss the rapid evolution of CAP proteins with varied subtypes in venomous animals. A treasure trove of information can be obtained by studying the CAP proteins in venomous animals; hence, it is necessary to explore these proteins further.

The roles of balancing selection and recombination in the evolution of rattlesnake venom

The origin of snake venom involved duplication and recruitment of non-venom genes into venom systems. Several studies have predicted that directional positive selection has governed this process. Venom composition varies substantially across snake species and venom phenotypes are locally adapted to prey, leading to coevolutionary interactions between predator and prey. Venom origins and contemporary snake venom evolution may therefore be driven by fundamentally different selection regimes, yet investigations of population-level patterns of selection have been limited. Here, we use whole-genome data from 68 rattlesnakes to test hypotheses about the factors that drive genomic diversity and differentiation in major venom gene regions. We show that selection has resulted in long-term maintenance of genetic diversity within and between species in multiple venom gene families. Our findings are inconsistent with a dominant role of directional positive selection and instead support a role of long-term balancing selection in shaping venom evolution. We also detect rapid decay of linkage disequilibrium due to high recombination rates in venom regions, suggesting that venom genes have reduced selective interference with nearby loci, including other venom paralogues. Our results provide an example of long-term balancing selection that drives trans-species polymorphism and help to explain how snake venom keeps pace with prey resistance.

Venom Variation of Neonate and Adult Chinese Cobras in Captivity Concerning Their Foraging Strategies

The venom and transcriptome profile of the captive Chinese cobra (Naja atra) is not characterized until now. Here, LC-MS/MS and illumine technology were used to unveil the venom and trascriptome of neonates and adults N. atra specimens. In captive Chinese cobra, 98 co-existing transcripts for venom-related proteins was contained. A total of 127 proteins belong to 21 protein families were found in the profile of venom. The main components of snake venom were three finger toxins (3-FTx), snake venom metalloproteinase (SVMP), cysteine-rich secretory protein (CRISP), cobra venom factor (CVF), and phosphodiesterase (PDE). During the ontogenesis of captive Chinese cobra, the rearrangement of snake venom composition occurred and with obscure gender difference. CVF, 3-FTx, PDE, phospholipase A2 (PLA2) in adults were more abundant than neonates, while SVMP and CRISP in the neonates was richer than the adults. Ontogenetic changes in the proteome of Chinese cobra venom reveals different strategies for handling prey. The levels of different types of toxin families were dramatically altered in the wild and captive specimens. Therefore, we speculate that the captive process could reshape the snake venom composition vigorously. The clear comprehension of the composition of Chinese cobra venom facilitates the understanding of the mechanism of snakebite intoxication and guides the preparation and administration of traditional antivenom and next-generation drugs for snakebite.

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.

The structural and functional divergence of a neglected three-finger toxin subfamily in lethal elapids

Bungarus multicinctus is a widely distributed and medically important elapid snake that produces lethal neurotoxic venom. To study and enhance existing antivenom, we explore the complete repertoire of its toxin genes based on de novo chromosome-level assembly and multi-tissue transcriptome data. Comparative genomic analyses suggest that the three-finger toxin family (3FTX) may evolve through the neofunctionalization of flanking LY6E. A long-neglected 3FTX subfamily (i.e., MKA-3FTX) is also investigated. Only one MKA-3FTX gene, which evolves a different protein conformation, is under positive selection and actively transcribed in the venom gland, functioning as a major toxin effector together with MKT-3FTX subfamily homologs. Furthermore, this lethal snake may acquire self-resistance to its β-bungarotoxin via amino acid replacements on fast-evolving KCNA2. This study provides valuable resources for further evolutionary and structure-function studies of snake toxins, which are fundamental for the development of effective antivenoms and drug candidates.

Origins, genomic structure and copy number variation of snake venom myotoxins

Crotamine, myotoxin a and homologs are short peptides that often comprise major fractions of rattlesnake venoms and have been extensively studied for their bioactive properties. These toxins are thought to be important for rapidly immobilizing mammalian prey and are implicated in serious, and sometimes fatal, responses to envenomation in humans. While high quality reference genomes for multiple venomous snakes are available, the loci that encode myotoxins have not been successfully assembled in any existing genome assembly. Here, we integrate new and existing genomic and transcriptomic data from the Prairie Rattlesnake (Crotalus viridis viridis) to reconstruct, characterize, and infer the chromosomal locations of myotoxin-encoding loci. We integrate long-read transcriptomics (Pacific Bioscience's Iso-Seq) and short-read RNA-seq to infer gene sequence diversity and characterize patterns of myotoxin and paralogous β-defensin expression across multiple tissues. We also identify two long non-coding RNA sequences which both encode functional myotoxins, demonstrating a newly discovered source of venom coding sequence diversity. We also integrate long-range mate-pair chromatin contact data and linked-read sequencing to infer the structure and chromosomal locations of the three myotoxin-like loci. Further, we conclude that the venom-associated myotoxin is located on chromosome 1 and is adjacent to non-venom paralogs. Consistent with this locus contributing to venom composition, we find evidence that the promoter of this gene is selectively open in venom gland tissue and contains transcription factor binding sites implicated in broad trans-regulatory pathways that regulate snake venoms. This study provides the best genomic reconstruction of myotoxin loci to date and raises questions about the physiological roles and interplay between myotoxin and related genes, as well as the genomic origins of snake venom variation.

The chemistry of snake venom and its medicinal potential

The fascination and fear of snakes dates back to time immemorial, with the first scientific treatise on snakebite envenoming, the Brooklyn Medical Papyrus, dating from ancient Egypt. Owing to their lethality, snakes have often been associated with images of perfidy, treachery and death. However, snakes did not always have such negative connotations. The curative capacity of venom has been known since antiquity, also making the snake a symbol of pharmacy and medicine. Today, there is renewed interest in pursuing snake-venom-based therapies. This Review focuses on the chemistry of snake venom and the potential for venom to be exploited for medicinal purposes in the development of drugs. The mixture of toxins that constitute snake venom is examined, focusing on the molecular structure, chemical reactivity and target recognition of the most bioactive toxins, from which bioactive drugs might be developed. The design and working mechanisms of snake-venom-derived drugs are illustrated, and the strategies by which toxins are transformed into therapeutics are analysed. Finally, the challenges in realizing the immense curative potential of snake venom are discussed, and chemical strategies by which a plethora of new drugs could be derived from snake venom are proposed.

Genetic characterization of potential venom resistance proteins in California ground squirrels ( Otospermophilus beecheyi ) using transcriptome analyses

Understanding the molecular basis of adaptations in coevolving species requires identifying the genes that underlie reciprocally selected phenotypes, such as those involved in venom in snakes and resistance to the venom in their prey. In this regard, California ground squirrels (CGS; Otospermophilus beecheyi) are eaten by northern Pacific rattlesnakes (Crotalus oreganus oreganus), but individual squirrels may still show substantial resistance to venom and survive bites. A recent study using proteomics identified venom interactive proteins (VIPs) in the blood serum of CGS. These VIPs represent possible resistance proteins, but the sequences of genes encoding them are unknown despite the value of such data to molecular studies of coevolution. To address this issue, we analyzed a de novo assembled transcriptome from CGS liver tissue-where many plasma proteins are synthesized-and other tissues from this species. We then examined VIP sequences in terms of three characteristics that identify them as possible resistance proteins: evidence for positive selection, high liver expression, and nonsynonymous variation across CGS populations. Based on these characteristics, we identified five VIPs (i.e., α-2-macroglobulin, α-1-antitrypsin-like protein GS55-LT, apolipoprotein A-II, hibernation-associated plasma protein HP-20, and hibernation-associated plasma protein HP-27) as the most likely candidates for resistance proteins among VIPs identified to date. Four of these proteins have been previously implicated in conferring resistance to the venom in mammals, validating our approach. When combined with the detailed information available for rattlesnake venom proteins, these results set the stage for future work focused on understanding coevolutionary interactions at the molecular level between these species.

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.

Genomic Confirmation of the P-IIIe Subclass of Snake Venom Metalloproteinases and Characterisation of Its First Member, a Disintegrin-Like/Cysteine-Rich Protein

Disintegrin-like/cysteine-rich (DC) proteins have long been regarded just as products of proteolysis of P-III snake venom metalloproteinases (SVMPs). However, here we demonstrate that a DC protein from the venom of Vipera ammodytes (Vaa; nose-horned viper), VaaMPIII-3, is encoded per se by a P-III SVMP-like gene that has a deletion in the region of the catalytic metalloproteinase domain and in part of the non-catalytic disintegrin-like domain. In this way, we justify the proposal of the introduction of a new subclass P-IIIe of SVMP-derived DC proteins. We purified VaaMPIII-3 from the venom of Vaa in a series of chromatographic steps. A covalent chromatography step based on thiol-disulphide exchange revealed that VaaMPIII-3 contains an unpaired Cys residue. This was demonstrated to be Cys6 in about 90% and Cys19 in about 10% of the VaaMPIII-3 molecules. We further constructed a three-dimensional homology model of VaaMPIII-3. From this model, it is evident that both Cys6 and Cys19 can pair with Cys26, which suggests that the intramolecular thiol-disulphide exchange has a regulatory function. VaaMPIII-3 is an acidic 21-kDa monomeric glycoprotein that exists in at least six N-glycoforms, with isoelectric points ranging from pH 4.5 to 5.1. Consistent with the presence of an integrin-binding motif in its sequence, SECD, VaaMPIII-3 inhibited collagen-induced platelet aggregation. It also inhibited ADP- and arachidonic-acid-induced platelet aggregation, but not ristocetin-induced platelet agglutination and the blood coagulation cascade.

Recurrent Co-Option and Recombination of Cytokine and Three Finger Proteins in Multiple Reproductive Tissues Throughout Salamander Evolution

Reproductive proteins evolve at unparalleled rates, resulting in tremendous diversity of both molecular composition and biochemical function between gametes of different taxonomic clades. To date, the proteomic composition of amphibian gametes is largely a molecular mystery, particularly for Urodeles (salamanders and newts) for which few genomic-scale resources exist. In this study, we provide the first detailed molecular characterization of gametes from two salamander species (Plethodon shermani and Desmognathus ocoee) that are models of reproductive behavior. Long-read PacBio transcriptome sequencing of testis and ovary of both species revealed sex-specific expression of many genes common to vertebrate gametes, including a similar expression profile to the egg coat genes of Xenopus oocytes. In contrast to broad conservation of oocyte genes, major testis transcripts included paralogs of salamander-specific courtship pheromones (PRF, PMF, and SPF) that were confirmed as major sperm proteins by mass spectrometry proteomics. Sperm-specific paralogs of PMF and SPF are likely the most abundant secreted proteins in P. shermani and D. ocoee, respectively. In contrast, sperm PRF lacks a signal peptide and may be expressed in cytoplasm. PRF pheromone genes evolved independently multiple times by repeated gene duplication of sperm PRF genes with signal peptides recovered through recombination with PMF genes. Phylogenetic analysis of courtship pheromones and their sperm paralogs support that each protein family evolved for these two reproductive contexts at distinct evolutionary time points between 17 and 360 million years ago. Our combined phylogenetic, transcriptomic and proteomic analyses of plethodontid reproductive tissues support that the recurrent co-option and recombination of TFPs and cytokine-like proteins have been a novel driving force throughout salamander evolution and reproduction.

Taxon-selective venom variation in adult and neonate Daboia russelii (Russell's Viper), and antivenom efficacy

Major variations in venom composition can occur between juvenile and adult venomous snakes. However, due to logistical constraints, antivenoms are produced using adult venoms in immunising mixtures, possibly resulting in limited neutralisation of juvenile snake venoms. Daboia russelii is one of the leading causes of snakebite death across South Asia. Its venom is potently procoagulant, causing stroke in prey animals but causing in humans consumptive coagulopathy-a net anticoagulant state-and sometimes death resulting from hemorrhage. In this in vitro study, we compared the venom activity of-and antivenom efficacy against-six 2-week-old D. russelii relative to that of their parents. Using a coagulation analyser, we quantified the relative coagulotoxicity of these venoms in human, avian, and amphibian plasma. The overall potency on human plasma was similar across all adult and neonate venoms, and SII (Serum Institute of India) antivenom was equipotent in neutralising these coagulotoxic effects. In addition, all venoms were also similar in their action upon avian plasma. In contrast, the neonate venoms were more potent on amphibian plasma, suggesting amphibians make up a larger proportion of neonate diet than adult diet. A similar venom potency in human and avian plasmas but varying selectivity for amphibian plasma suggests ontogenetic differences in toxin isoforms within the factor X or factor V activating classes, thereby providing a testable hypothesis for future transcriptomics work. By providing insights into the functional venom differences between adult and neonate D. russelii venoms, we hope to inform clinical treatment of patients envenomated by this deadly species and to shed new light on the natural history of these extremely medically important snakes.

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.

Co-option of the same ancestral gene family gave rise to mammalian and reptilian toxins

Background

Evolution can occur with surprising predictability when organisms face similar ecological challenges. For most traits, it is difficult to ascertain whether this occurs due to constraints imposed by the number of possible phenotypic solutions or because of parallel responses by shared genetic and regulatory architecture. Exceptionally, oral venoms are a tractable model of trait evolution, being largely composed of proteinaceous toxins that have evolved in many tetrapods, ranging from reptiles to mammals. Given the diversity of venomous lineages, they are believed to have evolved convergently, even though biochemically similar toxins occur in all taxa.

Results

Here, we investigate whether ancestral genes harbouring similar biochemical activity may have primed venom evolution, focusing on the origins of kallikrein-like serine proteases that form the core of most vertebrate oral venoms. Using syntenic relationships between genes flanking known toxins, we traced the origin of kallikreins to a single locus containing one or more nearby paralogous kallikrein-like clusters. Additionally, phylogenetic analysis of vertebrate serine proteases revealed that kallikrein-like toxins in mammals and reptiles are genetically distinct from non-toxin ones.

Conclusions

Given the shared regulatory and genetic machinery, these findings suggest that tetrapod venoms evolved by co-option of proteins that were likely already present in saliva. We term such genes ‘toxipotent’—in the case of salivary kallikreins they already had potent vasodilatory activity that was weaponized by venomous lineages. Furthermore, the ubiquitous distribution of kallikreins across vertebrates suggests that the evolution of envenomation may be more common than previously recognized, blurring the line between venomous and non-venomous animals.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01191-1.

Snake Venoms in Diagnostic Hemostasis and Thrombosis

Snake venoms have evolved primarily to immobilize and kill prey, and consequently, they contain some of the most potent natural toxins. Part of that armory is a range of hemotoxic components that affect every area of hemostasis, which we have harnessed to great effect in the study and diagnosis of hemostatic disorders. The most widely used are those that affect coagulation, such as thrombin-like enzymes unaffected by heparin and direct thrombin inhibitors, which can help confirm or dispute their presence in plasma. The liquid gold of coagulation activators is Russell's viper venom, since it contains activators of factor X and factor V. It is used in a range of clotting-based assays, such as assessment of factor X and factor V deficiencies, protein C and protein S deficiencies, activated protein C resistance, and probably the most important test for lupus anticoagulants, the dilute Russell's viper venom time. Activators of prothrombin, such as oscutarin C from Coastal Taipan venom and ecarin from saw-scaled viper venom, are employed in prothrombin activity assays and lupus anticoagulant detection, and ecarin has a valuable role in quantitative assays of direct thrombin inhibitors. Snake venoms affecting primary hemostasis include botrocetin from the jararaca, which can be used to assay von Willebrand factor activity, and convulxin from the cascavel, which can be used to detect deficiency of the platelet collagen receptor, glycoprotein VI. This article takes the reader to every area of the diagnostic hemostasis laboratory to appreciate the myriad applications of snake venoms available in diagnostic practice.

From Poison to Promise: The Evolution of Tetrodotoxin and Its Potential as a Therapeutic

Tetrodotoxin (TTX) is a potent neurotoxin that was first identified in pufferfish but has since been isolated from an array of taxa that host TTX-producing bacteria. However, determining its origin, ecosystem roles, and biomedical applications has challenged researchers for decades. Recognized as a poison and for its lethal effects on humans when ingested, TTX is primarily a powerful sodium channel inhibitor that targets voltage-gated sodium channels, including six of the nine mammalian isoforms. Although lethal doses for humans range from 1.5-2.0 mg TTX (blood level 9 ng/mL), when it is administered at levels far below LD₅₀, TTX exhibits therapeutic properties, especially to treat cancer-related pain, neuropathic pain, and visceral pain. Furthermore, TTX can potentially treat a variety of medical ailments, including heroin and cocaine withdrawal symptoms, spinal cord injuries, brain trauma, and some kinds of tumors. Here, we (i) describe the perplexing evolution and ecology of tetrodotoxin, (ii) review its mechanisms and modes of action, and (iii) offer an overview of the numerous ways it may be applied as a therapeutic. There is much to be explored in these three areas, and we offer ideas for future research that combine evolutionary biology with therapeutics. The TTX system holds great promise as a therapeutic and understanding the origin and chemical ecology of TTX as a poison will only improve its general benefit to humanity.

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.

Experimental Manipulation of Corticosterone Does Not Affect Venom Composition or Functional Activity in Free-Ranging Rattlesnakes

AbstractVenom is an integral feeding trait in many animal species. Although venom often varies ontogenetically, little is known about the proximate physiological mediators of venom variation within individuals. The glucocorticoid hormone corticosterone (CORT) can alter the transcription and activation of proteins, including homologues of snake venom components such as snake venom metalloproteinases (SVMPs) and phospholipase A₂ (PLA₂). CORT is endogenously produced by snakes, varies seasonally and also in response to stress, and is a candidate endogenous mediator of changes in venom composition and functional activity. Here, we tested the hypothesis that CORT induces changes in snake venom by sampling the venom of wild adult rattlesnakes before and after they were treated with either empty (control) or CORT-filled (treatment) Silastic implants. We measured longitudinal changes in whole-venom composition, whole-venom total protein content, and enzymatic activity of SVMP and PLA₂ components of venom. We also assessed the within-individual repeatability of venom components. Despite successfully elevating plasma CORT in the treatment group, we found no effect of CORT treatment or average plasma CORT level on any venom variables measured. Except for total protein content, venom components were highly repeatable within individuals ([Formula: see text]). Our results indicate that the effects of CORT, a hormone commonly associated with stress and metabolic functions, in adult rattlesnake venom are negligible. Our findings bode well for venom researchers and biomedical applications that rely on the consistency of venoms produced from potentially stressed individuals and provide an experimental framework for future studies of proximate mediators of venom variation across an individual's life span.

Origin and Evolution of the Gene Family of Proteinaceous Pheromones, the Exocrine Gland-Secreting Peptides, in Rodents

The exocrine-gland secreting peptide (ESP)gene family encodes proteinaceous pheromones that are recognized by the vomeronasal organ in mice. For example, ESP1 is a male pheromone secreted in tear fluid that regulates socio-sexual behavior, and ESP22 is a juvenile pheromone that suppresses adult sexual behavior. The family consists of multiple genes and has been identified only in mouse and rat genomes. The coding region of a mouse ESP gene is separated into two exons, each encoding signal and mature sequences. Here, we report the origin and evolution of the ESP gene family. ESP genes were found only in the Muridea and Cricetidae families of rodents, suggesting a recent origin of ESP genes in the common ancestor of murids and cricetids. ESP genes show a great diversity in number, length, and sequence among different species as well as mouse strains. Some ESPs in rats and golden hamsters are expressed in the lacrimal gland and the salivary gland. We also found that a mature sequence of an ESP gene showed overall sequence similarity to the α-globin gene. The ancestral ESP gene seems to be generated by recombination of a retrotransposed α-globin gene with the signal-encoding exon of the CRISP2 gene located adjacent to the ESP gene cluster. This study provides an intriguing example of molecular tinkering in rapidly evolving species-specific proteinaceous pheromone genes.

The evolutionary dynamics of venom toxins made by insects and other animals

Animal venoms are recognised as unique biological systems in which to study molecular evolution. Venom use has evolved numerous times among the insects, and insects today use venom to capture prey, defend themselves from predators, or to subdue and modulate host responses during parasitism. However, little is known about most insect venom toxins or the mode and tempo by which they evolve. Here, I review the evolutionary dynamics of insect venom toxins, and argue that insects offer many opportunities to examine novel aspects of toxin evolution. The key questions addressed are: How do venomous animals evolve from non-venomous animals, and how does this path effect the composition and pharmacology of the venom? What genetic processes (gene duplication, co-option, neofunctionalisation) are most important in toxin evolution? What kinds of selection pressures are acting on toxin-encoding genes and their cognate targets in envenomated animals? The emerging evidence highlights that venom composition and pharmacology adapts quickly in response to changing selection pressures resulting from new ecological interactions, and that such evolution occurs through a stunning variety of genetic mechanisms. Insects offer many opportunities to investigate the evolutionary dynamics of venom toxins due to their evolutionary history rich in venom-related adaptations, and their quick generation time and suitability for culture in the laboratory.

Tracking the recruitment and evolution of snake toxins using the evolutionary context provided by the Bothrops jararaca genome

Venom is a key adaptive innovation in snakes, and how nonvenom genes were co-opted to become part of the toxin arsenal is a significant evolutionary question. While this process has been investigated through the phylogenetic reconstruction of toxin sequences, evidence provided by the genomic context of toxin genes remains less explored. To investigate the process of toxin recruitment, we sequenced the genome of Bothrops jararaca, a clinically relevant pitviper. In addition to producing a road map with canonical structures of genes encoding 12 toxin families, we inferred most of the ancestral genes for their loci. We found evidence that 1) snake venom metalloproteinases (SVMPs) and phospholipases A2 (PLA2) have expanded in genomic proximity to their nonvenomous ancestors; 2) serine proteinases arose by co-opting a local gene that also gave rise to lizard gilatoxins and then expanded; 3) the bradykinin-potentiating peptides originated from a C-type natriuretic peptide gene backbone; and 4) VEGF-F was co-opted from a PGF-like gene and not from VEGF-A. We evaluated two scenarios for the original recruitment of nontoxin genes for snake venom: 1) in locus ancestral gene duplication and 2) in locus ancestral gene direct co-option. The first explains the origins of two important toxins (SVMP and PLA2), while the second explains the emergence of a greater number of venom components. Overall, our results support the idea of a locally assembled venom arsenal in which the most clinically relevant toxin families expanded through posterior gene duplications, regardless of whether they originated by duplication or gene co-option.

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.

Structural and bioinformatic analyses of Azemiops venom serine proteases reveal close phylogeographic relationships to pitvipers from eastern China and the New World

The semi-fossil and pit-less Azemiops feae is possibly the most primitive crotalid species. Here, we have cloned and sequenced cDNAs encoding four serine proteases (vSPs) from the venom glands of Chinese A. feae. Full amino-acid sequences of the major vSP (designated as AzKNa) and three minor vSPs (designated as AzKNb, AzKNc and Az-PA) were deduced. Using Protein-BLAST search, the ten most-similar vSPs for each Azemiops vSP have been selected for multiple sequence alignment, and all the homologs are crotalid vSPs. The results suggest that the A. feae vSPs are structurally most like those of eastern-Chinese Gloydius, Viridovipera, Protobothrops and North American pitvipers, and quite different from more-specialized vSPs such as Agkistrodon venom Protein-C activators. The vSPs from Chinese A. feae and those from Vietnamese A. feae show significant sequence variations. AzKNa is acidic and contains six potential N-glycosylation sites and its surface-charge distribution differs greatly from that of AzKNb, as revealed by 3D-modeling. AzKNb and AzKNc do not contain N-glycosylation sites although most of their close homologs contain one or two. Az-PA belongs to the plasminogen-activator subtype with a conserved N²⁰-glycosylation site. The evolution of this subtype of vSPs in Azemiops and related pitvipers has been traced by phylogenetic analysis.

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.

Patagonin-CRISP: Antimicrobial Activity and Source of Antimicrobial Molecules in Duvernoy's Gland Secretion (Philodryas patagoniensis Snake)

Snake venom contains a variety of toxins with a range of biological activity, among these toxins cysteine-rich secreted proteins (CRISPs) can be found. The proteins of this family have masses of 20-30 kDa and display homologous amino acid sequences containing 16 cysteine residues, forming eight disulfide bonds. Some of these proteins have been explored, characterized, and described in terms of their activity; however, little is known about their range of activities. A search for new antimicrobial molecules is ongoing, as the number of microbial strains resistant to available antibiotics is increasing. We identified antimicrobial activity in the secretion of Duvernoy's gland of the rear-fanged Philodryas patagoniensis. Fractions of this venom were subjected to reverse-phase high performance liquid chromatography and analyzed to determine their antimicrobial activity with a liquid broth inhibition assay. One of the fractions presented activity against a Gram-negative bacterium and a filamentous fungus. This fraction was analyzed with LC-MS/MS, and a protein of 24,848.8 Da was identified. Database searches allowed us to identify it as a CRISP due to the presence of some unique fragments in the molecule. We called it patagonin-CRISP, as the same protein in the venom of P. patagoniensis had previously been characterized as having a different biological activity. Patagonin-CRISP presented activity at very low concentrations and showed no cytotoxic activity. This is the first time that antimicrobial activity has been identified for P. patagoniensis venom or for a CRISP family protein.

Prospective randomized clinical study of a new topical formulation for face wrinkle reduction and dermal regeneration

BACKGROUND

A new topical formulation (TF) based on 3 main lifting components has been developed to reduce superficial facial wrinkles.

OBJECTIVES

Determine the effectiveness of this new TF in reducing superficial face wrinkles and restructuring the dermal matrix.

METHODS

Women, aged 30-65 y.o. with moderate to severe crow's feet wrinkles were included.

EXCLUSION CRITERIA

men; younger than 30 or older than 65 years old; smokers. Patients received 15 IU of botulinum toxin on crow's feet and 2 creams. Fifty patients (Group 1) applied the TP (Product A) and 50 (Group 2) a placebo (Product B). Assessments were made by digital macro-photography's, Antera 3D, and a patient satisfaction questionnaire.

RESULTS

From April to June 2019, 100 women were enrolled in the study and were divided into two homogeneous groups. No major or minor side effects were reported. In group 1, wrinkles, texture, static and dynamic crow's feet wrinkles improved significantly at 3 and 6 months. Patients were very satisfied at 3 months and satisfied at 6 months. In group 2, wrinkles and texture improved significantly at 3 months but did not improve at 6 months. Static and dynamic crow's feet wrinkles improved significantly at 1 and 3 months but did not improve significantly at 6 months.

CONCLUSIONS

Our prospective and randomized study has shown that the new TF is safe and effective in reducing superficial face wrinkles and producing dermal regeneration. It, therefore, prolongs the duration of the botulinum toxin. Further controlled study would be necessary to compare the new TF to neurotoxin treatment, or its action alone.

The functions of CAP superfamily proteins in mammalian fertility and disease

BACKGROUND

Members of the cysteine-rich secretory proteins (CRISPS), antigen 5 (Ag5) and pathogenesis-related 1 (Pr-1) (CAP) superfamily of proteins are found across the bacterial, fungal, plant and animal kingdoms. Although many CAP superfamily proteins remain poorly characterized, over the past decade evidence has accumulated, which provides insights into the functional roles of these proteins in various processes, including fertilization, immune defence and subversion, pathogen virulence, venom toxicology and cancer biology.

OBJECTIVE AND RATIONALE

The aim of this article is to summarize the current state of knowledge on CAP superfamily proteins in mammalian fertility, organismal homeostasis and disease pathogenesis.

SEARCH METHODS

The scientific literature search was undertaken via PubMed database on all articles published prior to November 2019. Search terms were based on following keywords: 'CAP superfamily', 'CRISP', 'Cysteine-rich secretory proteins', 'Antigen 5', 'Pathogenesis-related 1', 'male fertility', 'CAP and CTL domain containing', 'CRISPLD1', 'CRISPLD2', 'bacterial SCP', 'ion channel regulator', 'CatSper', 'PI15', 'PI16', 'CLEC', 'PRY proteins', 'ASP proteins', 'spermatogenesis', 'epididymal maturation', 'capacitation' and 'snake CRISP'. In addition to that, reference lists of primary and review article were reviewed for additional relevant publications.

OUTCOMES

In this review, we discuss the breadth of knowledge on CAP superfamily proteins with regards to their protein structure, biological functions and emerging significance in reproduction, health and disease. We discuss the evolution of CAP superfamily proteins from their otherwise unembellished prokaryotic predecessors into the multi-domain and neofunctionalized members found in eukaryotic organisms today. At least in part because of the rapid evolution of these proteins, many inconsistencies in nomenclature exist within the literature. As such, and in part through the use of a maximum likelihood phylogenetic analysis of the vertebrate CRISP subfamily, we have attempted to clarify this confusion, thus allowing for a comparison of orthologous protein function between species. This framework also allows the prediction of functional relevance between species based on sequence and structural conservation.

WIDER IMPLICATIONS

This review generates a picture of critical roles for CAP proteins in ion channel regulation, sterol and lipid binding and protease inhibition, and as ligands involved in the induction of multiple cellular processes.

How the Toxin got its Toxicity

Venom systems are functional and ecological traits, typically used by one organism to subdue or deter another. A predominant subset of their constituent molecules—“toxins”—share this ecological function and are therefore molecules that mediate interactions between organisms. Such molecules have been referred to as “exochemicals.” There has been debate within the field of toxinology concerning the evolutionary pathways leading to the “recruitment” of a gene product for a toxic role within venom. We review these discussions and the evidence interpreted in support of alternate pathways, along with many of the most popular models describing the origin of novel molecular functions in general. We note that such functions may arise with or without gene duplication occurring and are often the consequence of a gene product encountering a novel “environment,” i.e., a range of novel partners for molecular interaction. After stressing the distinction between “activity” and “function,” we describe in detail the results of a recent study which reconstructed the evolutionary history of a multigene family that has been recruited as a toxin and argue that these results indicate that a pluralistic approach to understanding the origin of novel functions is advantageous. This leads us to recommend that an expansive approach be taken to the definition of “neofunctionalization”—simply the origins of a novel molecular function by any process—and “recruitment”—the “weaponization” of a molecule via the acquisition of a toxic function in venom, by any process. Recruitment does not occur at the molecular level or even at the level of gene expression, but only when a confluence of factors results in the ecological deployment of a physiologically active molecule as a toxin. Subsequent to recruitment, the evolutionary regime of a gene family may shift into a more dynamic form of “birth-and-death.” Thus, recruitment leads to a form of “downwards causation,” in which a change at the ecological level at which whole organisms interact leads to a change in patterns of evolution at the genomic level.

Pick Your Poison: Molecular Evolution of Venom Proteins in Asilidae (Insecta: Diptera)

Robber flies are an understudied family of venomous, predatory Diptera. With the recent characterization of venom from three asilid species, it is possible, for the first time, to study the molecular evolution of venom genes in this unique lineage. To accomplish this, a novel whole-body transcriptome of Eudioctria media was combined with 10 other publicly available asiloid thoracic or salivary gland transcriptomes to identify putative venom gene families and assess evidence of pervasive positive selection. A total of 348 gene families of sufficient size were analyzed, and 33 of these were predicted to contain venom genes. We recovered 151 families containing homologs to previously described venom proteins, and 40 of these were uniquely gained in Asilidae. Our gene family clustering suggests that many asilidin venom gene families are not natural groupings, as delimited by previous authors, but instead form multiple discrete gene families. Additionally, robber fly venoms have relatively few sites under positive selection, consistent with the hypothesis that the venoms of older lineages are dominated by negative selection acting to maintain toxic function.

Neuropeptide signalling systems - An underexplored target for venom drug discovery

Neuropeptides are signalling molecules mainly secreted from neurons that act as neurotransmitters or peptide hormones to affect physiological processes and modulate behaviours. In humans, neuropeptides are implicated in numerous diseases and understanding their role in physiological processes and pathologies is important for therapeutic development. Teasing apart the (patho)physiology of neuropeptides remains difficult due to ligand and receptor promiscuity and the complexity of the signalling pathways. The current approach relies on a pharmacological toolbox of agonists and antagonists displaying high selectivity for independent receptor subtypes, with the caveat that only few selective ligands have been discovered or developed. Animal venoms represent an underexplored source for novel receptor subtype-selective ligands that could aid in dissecting human neuropeptide signalling systems. Multiple endogenous-like neuropeptides as well as peptides acting on neuropeptide receptors are present in venoms. In this review, we summarise current knowledge on neuropeptides and discuss venoms as a source for ligands targeting neuropeptide signalling systems.

Proteomic Investigations of Two Pakistani Naja Snake Venoms Species Unravel the Venom Complexity, Posttranslational Modifications, and Presence of Extracellular Vesicles

Latest advancement of omics technologies allows in-depth characterization of venom compositions. In the present work we present a proteomic study of two snake venoms of the genus Naja i.e., Naja naja (black cobra) and Naja oxiana (brown cobra) of Pakistani origin. The present study has shown that these snake venoms consist of a highly diversified proteome. Furthermore, the data also revealed variation among closely related species. High throughput mass spectrometric analysis of the venom proteome allowed to identify for the N. naja venom 34 protein families and for the N. oxiana 24 protein families. The comparative evaluation of the two venoms showed that N. naja consists of a more complex venom proteome than N. oxiana venom. Analysis also showed N-terminal acetylation (N-ace) of a few proteins in both venoms. To the best of our knowledge, this is the first study revealing this posttranslational modification in snake venom. N-ace can shed light on the mechanism of regulation of venom proteins inside the venom gland. Furthermore, our data showed the presence of other body proteins, e.g., ankyrin repeats, leucine repeats, zinc finger, cobra serum albumin, transferrin, insulin, deoxyribonuclease-2-alpha, and other regulatory proteins in these venoms. Interestingly, our data identified Ras-GTpase type of proteins, which indicate the presence of extracellular vesicles in the venom. The data can support the production of distinct and specific anti-venoms and also allow a better understanding of the envenomation and mechanism of distribution of toxins. Data are available via ProteomeXchange with identifier PXD018726.

Snake three-finger α-neurotoxins and nicotinic acetylcholine receptors: molecules, mechanisms and medicine

Snake venom three-finger α-neurotoxins (α-3FNTx) act on postsynaptic nicotinic acetylcholine receptors (nAChRs) at the neuromuscular junction (NMJ) to produce skeletal muscle paralysis. The discovery of the archetypal α-bungarotoxin (α-BgTx), almost six decades ago, exponentially expanded our knowledge of membrane receptors and ion channels. This included the localisation, isolation and characterization of the first receptor (nAChR); and by extension, the pathophysiology and pharmacology of neuromuscular transmission and associated pathologies such as myasthenia gravis, as well as our understanding of the role of α-3FNTxs in snakebite envenomation leading to novel concepts of targeted treatment. Subsequent studies on a variety of animal venoms have yielded a plethora of novel toxins that have revolutionized molecular biomedicine and advanced drug discovery from bench to bedside. This review provides an overview of nAChRs and their subtypes, classification of α-3FNTxs and the challenges of typifying an increasing arsenal of structurally and functionally unique toxins, and the three-finger protein (3FP) fold in the context of the uPAR/Ly6/CD59/snake toxin superfamily. The pharmacology of snake α-3FNTxs including their mechanisms of neuromuscular blockade, variations in reversibility of nAChR interactions, specificity for nAChR subtypes or for distinct ligand-binding interfaces within a subtype and the role of α-3FNTxs in neurotoxic envenomation are also detailed. Lastly, a reconciliation of structure-function relationships between α-3FNTx and nAChRs, derived from historical mutational and biochemical studies and emerging atomic level structures of nAChR models in complex with α-3FNTxs is discussed.

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

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

The molecular basis of venom resistance in a rattlesnake-squirrel predator-prey system

Understanding how interspecific interactions mould the molecular basis of adaptations in coevolving species is a long-sought goal of evolutionary biology. Venom in predators and venom resistance proteins in prey are coevolving molecular phenotypes, and while venoms are highly complex mixtures it is unclear if prey respond with equally complex resistance traits. Here, we use a novel molecular methodology based on protein affinity columns to capture and identify candidate blood serum resistance proteins ("venom interactive proteins" [VIPs]) in California Ground Squirrels (Otospermophilus beecheyi) that interact with venom proteins from their main predator, Northern Pacific Rattlesnakes (Crotalus o. oreganus). This assay showed that serum-based resistance is both population- and species-specific, with serum proteins from ground squirrels showing higher binding affinities for venom proteins of local snakes compared to allopatric individuals. Venom protein specificity assays identified numerous and diverse candidate prey resistance VIPs but also potential targets of venom in prey tissues. Many specific VIPs bind to multiple snake venom proteins and, conversely, single venom proteins bind multiple VIPs, demonstrating that a portion of the squirrel blood serum "resistome" involves broad-based inhibition of nonself proteins and suggests that resistance involves a toxin scavenging mechanism. Analyses of rates of evolution of VIP protein homologues in related mammals show that most of these proteins evolve under purifying selection possibly due to molecular constraints that limit the evolutionary responses of prey to rapidly evolving snake venom proteins. Our method represents a general approach to identify specific proteins involved in co-evolutionary interactions between species at the molecular level.