The Deep Origin and Recent Loss of Venom Toxin Genes in Rattlesnakes

VenomCap: An exon-capture probe set for the targeted sequencing of snake venom genes

Snake venoms are complex mixtures of toxic proteins that hold significant medical, pharmacological and evolutionary interest. To better understand the genetic diversity underlying snake venoms, we developed VenomCap, a novel exon-capture probe set targeting toxin-coding genes from a wide range of elapid snakes, with a particular focus on the ecologically diverse and medically important subfamily Hydrophiinae. We tested the capture success of VenomCap across 24 species, representing all major elapid lineages. We included snake phylogenomic probes in the VenomCap capture set, allowing us to compare capture performance between venom and phylogenomic loci and to infer elapid phylogenetic relationships. We demonstrated VenomCap's ability to recover exons from ~1500 target markers, representing a total of 24 known venom gene families, which includes the dominant gene families found in elapid venoms. We find that VenomCap's capture results are robust across all elapids sampled, and especially among hydrophiines, with respect to measures of target capture success (target loci matched, sensitivity, specificity and missing data). As a cost-effective and efficient alternative to full genome sequencing, VenomCap can dramatically accelerate the sequencing and analysis of venom gene families. Overall, our tool offers a model for genomic studies on snake venom gene diversity and evolution that can be expanded for comprehensive comparisons across the other families of venomous snakes.

Exploring venom diversity in Mixcoatlus browni and Mixcoatlus barbouri: A comparative analysis of two rare Mexican snake species with crotoxin-like presence

The genus Mixcoatlus is composed of three species: Mixcoatlus barbouri, M. browni, and M. melanurus, of which the venom composition of M. melanurus, the most common species of the three, has only recently been described. However, very little is known about the natural history of M. barbouri and M. browni, and the venom composition of these two species has remained thus far unexplored. In this study we characterize the proteomic profiles and the main biochemical and toxic activities of these two venoms. Proteomic data obtained by shotgun analysis of whole venom identified 12 protein families for M. barbouri, and 13 for M. browni. The latter venom was further characterized by using a quantitative 'venomics' protocol, which revealed that it is mainly composed of 51.1 % phospholipases A2 (PLA2), 25.5 % snake venom serine proteases (SVSP), 4.6 % l-amino oxidases (LAO), and 3.6 % snake venom metalloproteases (SVMP), with lower percentages other six protein families. Both venoms contained homologs of the basic and acidic subunits of crotoxin. However, due to limitations in M. barbouri venom availability, we could only characterize the crotoxin-like protein of M. browni venom, which we have named Mixcoatlutoxin. It exhibited a lethal potency in mice like that described for classical rattlesnake crotoxins. These findings expand knowledge on the distribution of crotoxin-like heterodimeric proteins in viper snake species. Further investigation of the bioactivities of the venom of M. barbouri, on the other hand, remains necessary.

Hyperglycosylation impairs the inhibitory activity of rCdtPLI2, the first recombinant beta-phospholipase A2 inhibitor

Crotoxin, a phospholipase A2 (PLA2) complex and the major Crotalus venom component, is responsible for the main symptoms described in crotalic snakebite envenomings and a key target for PLA2 inhibitors (PLIs). PLIs comprise the alpha, beta and gamma families, and, due to a lack of reports on beta-PLIs, this study aimed to heterologously express CdtPLI2 from Crotalus durissus terrificus venom gland to improve the knowledge of the neglected beta-PLI family. Thereby, recombinant CdtPLI2 (rCdtPLI2) was produced in the eukaryotic Pichia pastoris system to keep some native post-translational modifications. rCdtPLI2 (~41 kDa) presents both N- and O-linked glycans. Alpha-mannosidase digested-rCdtPLI2 (1 mol) strongly inhibited (73%) CB-Cdc catalytic activity (5 moles), demonstrating that glycosylations performed by P. pastoris affect rCdtPLI2 action. Digested-rCdtPLI2 also inhibited PLA2s from diverse Brazilian snake venoms. Furthermore, rCdtPLI2 (1 mol) abolished the catalytic activity of Lmr-PLA2 (5 moles) and reduced the CTx-Cdc (5 moles) enzyme activity by 65%, suppressing basic and acidic snake venom PLA2s. Additionally, crotalic antivenom did not recognize rCdtPLI2, suggesting a lack of neutralization by antivenom antibodies. These findings demonstrate that studying snake venom components may reveal interesting novel molecules to be studied in the snakebite treatment and help to understand these underexplored inhibitors.

Molecular Insights into the Low Complexity Secreted Venom of Calliactis polypus

Sea anemones are venomous animals that rely on their venom for prey capture, defense against predators, and intraspecific competition. Currently, comprehensive molecular and evolutionary analyses of the toxin repertoire for sea anemones are limited by a lack of proteomic data for most species. In this study, proteo-transcriptomic analysis was used to expand our knowledge of the proteinaceous components of sea anemone venom by determining the secreted venom proteome of Calliactis polypus. Electromechanical stimulation was used to obtain the secreted venom of C. polypus. We identified a low complexity proteome that was dominated by toxins with similarity to known neurotoxins, as well as six novel toxin candidates. The novel putative toxin candidates were found to be taxonomically restricted to species from the superfamily Metridioidea. Furthermore, the secreted venom of C. polypus had only three putative toxins in common with the venom of acontia from the same species and little similarity with the secreted venom of closely related species. Overall, this demonstrates that regionalized and lineage-specific variability in toxin abundance is common among sea anemone species. Moreover, the limited complexity of the toxin repertoire found in C. polypus supports the idea that peptide neurotoxins make up the dominant toxin arsenal found in the venom of sea anemones.

Diversely evolved xibalbin variants from remipede venom inhibit potassium channels and activate PKA-II and Erk1/2 signaling

BACKGROUND

The identification of novel toxins from overlooked and taxonomically exceptional species bears potential for various pharmacological applications. The remipede Xibalbanus tulumensis, an underwater cave-dwelling crustacean, is the only crustacean for which a venom system has been described. Its venom contains several xibalbin peptides that have an inhibitor cysteine knot (ICK) scaffold.

RESULTS

Our screenings revealed that all tested xibalbin variants particularly inhibit potassium channels. Xib1 and xib13 with their eight-cysteine domain similar to spider knottins also inhibit voltage-gated sodium channels. No activity was noted on calcium channels. Expanding the functional testing, we demonstrate that xib1 and xib13 increase PKA-II and Erk1/2 sensitization signaling in nociceptive neurons, which may initiate pain sensitization. Our phylogenetic analysis suggests that xib13 either originates from the common ancestor of pancrustaceans or earlier while xib1 is more restricted to remipedes. The ten-cysteine scaffolded xib2 emerged from xib1, a result that is supported by our phylogenetic and machine learning-based analyses.

CONCLUSIONS

Our functional characterization of synthesized variants of xib1, xib2, and xib13 elucidates their potential as inhibitors of potassium channels in mammalian systems. The specific interaction of xib2 with Kv1.6 channels, which are relevant to treating variants of epilepsy, shows potential for further studies. At higher concentrations, xib1 and xib13 activate the kinases PKA-II and ERK1/2 in mammalian sensory neurons, suggesting pain sensitization and potential applications related to pain research and therapy. While tested insect channels suggest that all probably act as neurotoxins, the biological function of xib1, xib2, and xib13 requires further elucidation. A novel finding on their evolutionary origin is the apparent emergence of X. tulumensis-specific xib2 from xib1. Our study is an important cornerstone for future studies to untangle the origin and function of these enigmatic proteins as important components of remipede but also other pancrustacean and arthropod venoms.

Diverse Gene Regulatory Mechanisms Alter Rattlesnake Venom Gene Expression at Fine Evolutionary Scales

Understanding and predicting the relationships between genotype and phenotype is often challenging, largely due to the complex nature of eukaryotic gene regulation. A step towards this goal is to map how phenotypic diversity evolves through genomic changes that modify gene regulatory interactions. Using the Prairie Rattlesnake (Crotalus viridis) and related species, we integrate mRNA-seq, proteomic, ATAC-seq and whole-genome resequencing data to understand how specific evolutionary modifications to gene regulatory network components produce differences in venom gene expression. Through comparisons within and between species, we find a remarkably high degree of gene expression and regulatory network variation across even a shallow level of evolutionary divergence. We use these data to test hypotheses about the roles of specific trans-factors and cis-regulatory elements, how these roles may vary across venom genes and gene families, and how variation in regulatory systems drive diversity in venom phenotypes. Our results illustrate that differences in chromatin and genotype at regulatory elements play major roles in modulating expression. However, we also find that enhancer deletions, differences in transcription factor expression, and variation in activity of the insulator protein CTCF also likely impact venom phenotypes. Our findings provide insight into the diversity and gene-specificity of gene regulatory features and highlight the value of comparative studies to link gene regulatory network variation to phenotypic variation.

Assessing Target Specificity of the Small Molecule Inhibitor MARIMASTAT to Snake Venom Toxins: A Novel Application of Thermal Proteome Profiling

New treatments that circumvent the pitfalls of traditional antivenom therapies are critical to address the problem of snakebite globally. Numerous snake venom toxin inhibitors have shown promising cross-species neutralization of medically significant venom toxins in vivo and in vitro. The development of high-throughput approaches for the screening of such inhibitors could accelerate their identification, testing, and implementation and thus holds exciting potential for improving the treatments and outcomes of snakebite envenomation worldwide. Energetics-based proteomic approaches, including thermal proteome profiling and proteome integral solubility alteration (PISA) assays, represent "deep proteomics" methods for high throughput, proteome-wide identification of drug targets and ligands. In the following study, we apply thermal proteome profiling and PISA methods to characterize the interactions between venom toxin proteoforms in Crotalus atrox (Western Diamondback Rattlesnake) and the snake venom metalloprotease (SVMP) inhibitor marimastat. We investigate its venom proteome-wide effects and characterize its interactions with specific SVMP proteoforms, as well as its potential targeting of non-SVMP venom toxin families. We also compare the performance of PISA thermal window and soluble supernatant with insoluble precipitate using two inhibitor concentrations, providing the first demonstration of the utility of a sensitive high-throughput PISA-based approach to assess the direct targets of small molecule inhibitors for snake venom.

Phylogenomic discordance is driven by wide-spread introgression and incomplete lineage sorting during rapid species diversification within rattlesnakes (Viperidae: Crotalus and Sistrurus)

Phylogenomics allows us to uncover the historical signal of evolutionary processes through time and estimate phylogenetic networks accounting for these signals. Insight from genome-wide data further allows us to pinpoint the contributions to phylogenetic signal from hybridization, introgression, and ancestral polymorphism across the genome. Here we focus on how these processes have contributed to phylogenetic discordance among rattlesnakes (genera Crotalus and Sistrurus), a group for which there are numerous conflicting phylogenetic hypotheses based on a diverse array of molecular datasets and analytical methods. We address the instability of the rattlesnake phylogeny using genomic data generated from transcriptomes sampled from nearly all known species. These genomic data, analyzed with coalescent and network-based approaches, reveal numerous instances of rapid speciation where individual gene trees conflict with the species tree. Moreover, the evolutionary history of rattlesnakes is dominated by incomplete speciation and frequent hybridization, both of which have likely influenced past interpretations of phylogeny. We present a new framework in which the evolutionary relationships of this group can only be understood in light of genome-wide data and network-based analytical methods. Our data suggest that network radiations, like seen within the rattlesnakes, can only be understood in a phylogenomic context, necessitating similar approaches in our attempts to understand evolutionary history in other rapidly radiating species.

Higher evolutionary dynamics of gene copy number for Drosophila glue genes located near short repeat sequences

BACKGROUND

During evolution, genes can experience duplications, losses, inversions and gene conversions. Why certain genes are more dynamic than others is poorly understood. Here we examine how several Sgs genes encoding glue proteins, which make up a bioadhesive that sticks the animal during metamorphosis, have evolved in Drosophila species.

RESULTS

We examined high-quality genome assemblies of 24 Drosophila species to study the evolutionary dynamics of four glue genes that are present in D. melanogaster and are part of the same gene family - Sgs1, Sgs3, Sgs7 and Sgs8 - across approximately 30 millions of years. We annotated a total of 102 Sgs genes and grouped them into 4 subfamilies. We present here a new nomenclature for these Sgs genes based on protein sequence conservation, genomic location and presence/absence of internal repeats. Two types of glue genes were uncovered. The first category (Sgs1, Sgs3x, Sgs3e) showed a few gene losses but no duplication, no local inversion and no gene conversion. The second group (Sgs3b, Sgs7, Sgs8) exhibited multiple events of gene losses, gene duplications, local inversions and gene conversions. Our data suggest that the presence of short "new glue" genes near the genes of the latter group may have accelerated their dynamics.

CONCLUSIONS

Our comparative analysis suggests that the evolutionary dynamics of glue genes is influenced by genomic context. Our molecular, phylogenetic and comparative analysis of the four glue genes Sgs1, Sgs3, Sgs7 and Sgs8 provides the foundation for investigating the role of the various glue genes during Drosophila life.

ToxCodAn-Genome: an automated pipeline for toxin-gene annotation in genome assembly of venomous lineages

BACKGROUND

The rapid development of sequencing technologies resulted in a wide expansion of genomics studies using venomous lineages. This facilitated research focusing on understanding the evolution of adaptive traits and the search for novel compounds that can be applied in agriculture and medicine. However, the toxin annotation of genomes is a laborious and time-consuming task, and no consensus pipeline is currently available. No computational tool currently exists to address the challenges specific to toxin annotation and to ensure the reproducibility of the process.

RESULTS

Here, we present ToxCodAn-Genome, the first software designed to perform automated toxin annotation in genomes of venomous lineages. This pipeline was designed to retrieve the full-length coding sequences of toxins and to allow the detection of novel truncated paralogs and pseudogenes. We tested ToxCodAn-Genome using 12 genomes of venomous lineages and achieved high performance on recovering their current toxin annotations. This tool can be easily customized to allow improvements in the final toxin annotation set and can be expanded to virtually any venomous lineage. ToxCodAn-Genome is fast, allowing it to run on any personal computer, but it can also be executed in multicore mode, taking advantage of large high-performance servers. In addition, we provide a guide to direct future research in the venomics field to ensure a confident toxin annotation in the genome being studied. As a case study, we sequenced and annotated the toxin repertoire of Bothrops alternatus, which may facilitate future evolutionary and biomedical studies using vipers as models.

CONCLUSIONS

ToxCodAn-Genome is suitable to perform toxin annotation in the genome of venomous species and may help to improve the reproducibility of further studies. ToxCodAn-Genome and the guide are freely available at https://github.com/pedronachtigall/ToxCodAn-Genome.

Web of venom: exploration of big data resources in animal toxin research

Research on animal venoms and their components spans multiple disciplines, including biology, biochemistry, bioinformatics, pharmacology, medicine, and more. Manipulating and analyzing the diverse array of data required for venom research can be challenging, and relevant tools and resources are often dispersed across different online platforms, making them less accessible to nonexperts. In this article, we address the multifaceted needs of the scientific community involved in venom and toxin-related research by identifying and discussing web resources, databases, and tools commonly used in this field. We have compiled these resources into a comprehensive table available on the VenomZone website (https://venomzone.expasy.org/10897). Furthermore, we highlight the challenges currently faced by researchers in accessing and using these resources and emphasize the importance of community-driven interdisciplinary approaches. We conclude by underscoring the significance of enhancing standards, promoting interoperability, and encouraging data and method sharing within the venom research community.

Evolution towards simplicity in bacterial small heat shock protein system

Evolution can tinker with multi-protein machines and replace them with simpler single-protein systems performing equivalent functions in an equally efficient manner. It is unclear how, on a molecular level, such simplification can arise. With ancestral reconstruction and biochemical analysis, we have traced the evolution of bacterial small heat shock proteins (sHsp), which help to refold proteins from aggregates using either two proteins with different functions (IbpA and IbpB) or a secondarily single sHsp that performs both functions in an equally efficient way. Secondarily single sHsp evolved from IbpA, an ancestor specialized in strong substrate binding. Evolution of an intermolecular binding site drove the alteration of substrate binding properties, as well as the formation of higher-order oligomers. Upon two mutations in the α-crystallin domain, secondarily single sHsp interacts with aggregated substrates less tightly. Paradoxically, less efficient binding positively influences the ability of sHsp to stimulate substrate refolding, since the dissociation of sHps from aggregates is required to initiate Hsp70-Hsp100-dependent substrate refolding. After the loss of a partner, IbpA took over its role in facilitating the sHsp dissociation from an aggregate by weakening the interaction with the substrate, which became beneficial for the refolding process. We show that the same two amino acids introduced in modern-day systems define whether the IbpA acts as a single sHsp or obligatorily cooperates with an IbpB partner. Our discoveries illuminate how one sequence has evolved to encode functions previously performed by two distinct proteins.

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.

Prevalent bee venom genes evolved before the aculeate stinger and eusociality

BACKGROUND

Venoms, which have evolved numerous times in animals, are ideal models of convergent trait evolution. However, detailed genomic studies of toxin-encoding genes exist for only a few animal groups. The hyper-diverse hymenopteran insects are the most speciose venomous clade, but investigation of the origin of their venom genes has been largely neglected.

RESULTS

Utilizing a combination of genomic and proteo-transcriptomic data, we investigated the origin of 11 toxin genes in 29 published and 3 new hymenopteran genomes and compiled an up-to-date list of prevalent bee venom proteins. Observed patterns indicate that bee venom genes predominantly originate through single gene co-option with gene duplication contributing to subsequent diversification.

CONCLUSIONS

Most Hymenoptera venom genes are shared by all members of the clade and only melittin and the new venom protein family anthophilin1 appear unique to the bee lineage. Most venom proteins thus predate the mega-radiation of hymenopterans and the evolution of the aculeate stinger.

The Need for Next-Generation Antivenom for Snakebite Envenomation in India

The limitations posed by currently available antivenoms have emphasized the need for alternative treatments to counteract snakebite envenomation. Even though exact epidemiological data are lacking, reports have indicated that most global snakebite deaths are reported in India. Among the many problems associated with snakebite envenomation, issues related to the availability of safer and more efficient antivenoms are of primary concern. Since India has the highest number of global snakebite deaths, efforts should be made to reduce the burden associated with snakebite envenoming. Alternative methods, including aptamers, camel antivenoms, phage display techniques for generating high-affinity antibodies and antibody fragments, small-molecule inhibitors, and natural products, are currently being investigated for their effectiveness. These alternative methods have shown promise in vitro, but their in vivo effectiveness should also be evaluated. In this review, the issues associated with Indian polyvalent antivenoms in neutralizing venom components from geographically distant species are discussed in detail. In a nutshell, this review gives an overview of the current drawbacks of using animal-derived antivenoms and several alternative strategies that are currently being widely explored.

Sequence Divergence in Venom Genes Within and Between Montane Pitviper (Viperidae: Crotalinae: Cerrophidion) Species is Driven by Mutation-Drift Equilibrium

Snake venom can vary both among and within species. While some groups of New World pitvipers-such as rattlesnakes-have been well studied, very little is known about the venom of montane pitvipers (Cerrophidion) found across the Mesoamerican highlands. Compared to most well-studied rattlesnakes, which are widely distributed, the isolated montane populations of Cerrophidion may facilitate unique evolutionary trajectories and venom differentiation. Here, we describe the venom gland transcriptomes for populations of C. petlalcalensis, C. tzotzilorum, and C. godmani from Mexico, and a single individual of C. sasai from Costa Rica. We explore gene expression variation in Cerrophidion and sequence evolution of toxins within C. godmani specifically. Cerrophidion venom gland transcriptomes are composed primarily of snake venom metalloproteinases, phospholipase A[Formula: see text]s (PLA[Formula: see text]s), and snake venom serine proteases. Cerrophidion petlalcalensis shows little intraspecific variation; however, C. godmani and C. tzotzilorum differ significantly between geographically isolated populations. Interestingly, intraspecific variation was mostly attributed to expression variation as we did not detect signals of selection within C. godmani toxins. Additionally, we found PLA[Formula: see text]-like myotoxins in all species except C. petlalcalensis, and crotoxin-like PLA[Formula: see text]s in the southern population of C. godmani. Our results demonstrate significant intraspecific venom variation within C. godmani and C. tzotzilorum. The toxins of C. godmani show little evidence of directional selection where variation in toxin sequence is consistent with evolution under a model of mutation-drift equilibrium. Cerrophidion godmani individuals from the southern population may exhibit neurotoxic venom activity given the presence of crotoxin-like PLA[Formula: see text]s; however, further research is required to confirm this hypothesis.

The State of Squamate Genomics: Past, Present, and Future of Genome Research in the Most Speciose Terrestrial Vertebrate Order

Squamates include more than 11,000 extant species of lizards, snakes, and amphisbaenians, and display a dazzling diversity of phenotypes across their over 200-million-year evolutionary history on Earth. Here, we introduce and define squamates (Order Squamata) and review the history and promise of genomic investigations into the patterns and processes governing squamate evolution, given recent technological advances in DNA sequencing, genome assembly, and evolutionary analysis. We survey the most recently available whole genome assemblies for squamates, including the taxonomic distribution of available squamate genomes, and assess their quality metrics and usefulness for research. We then focus on disagreements in squamate phylogenetic inference, how methods of high-throughput phylogenomics affect these inferences, and demonstrate the promise of whole genomes to settle or sustain persistent phylogenetic arguments for squamates. We review the role transposable elements play in vertebrate evolution, methods of transposable element annotation and analysis, and further demonstrate that through the understanding of the diversity, abundance, and activity of transposable elements in squamate genomes, squamates can be an ideal model for the evolution of genome size and structure in vertebrates. We discuss how squamate genomes can contribute to other areas of biological research such as venom systems, studies of phenotypic evolution, and sex determination. Because they represent more than 30% of the living species of amniote, squamates deserve a genome consortium on par with recent efforts for other amniotes (i.e., mammals and birds) that aim to sequence most of the extant families in a clade.

Genomic, functional and structural analyses elucidate evolutionary innovation within the sea anemone 8 toxin family

BACKGROUND

The ShK toxin from Stichodactyla helianthus has established the therapeutic potential of sea anemone venom peptides, but many lineage-specific toxin families in Actiniarians remain uncharacterised. One such peptide family, sea anemone 8 (SA8), is present in all five sea anemone superfamilies. We explored the genomic arrangement and evolution of the SA8 gene family in Actinia tenebrosa and Telmatactis stephensoni, characterised the expression patterns of SA8 sequences, and examined the structure and function of SA8 from the venom of T. stephensoni.

RESULTS

We identified ten SA8-family genes in two clusters and six SA8-family genes in five clusters for T. stephensoni and A. tenebrosa, respectively. Nine SA8 T. stephensoni genes were found in a single cluster, and an SA8 peptide encoded by an inverted SA8 gene from this cluster was recruited to venom. We show that SA8 genes in both species are expressed in a tissue-specific manner and the inverted SA8 gene has a unique tissue distribution. While the functional activity of the SA8 putative toxin encoded by the inverted gene was inconclusive, its tissue localisation is similar to toxins used for predator deterrence. We demonstrate that, although mature SA8 putative toxins have similar cysteine spacing to ShK, SA8 peptides are distinct from ShK peptides based on structure and disulfide connectivity.

CONCLUSIONS

Our results provide the first demonstration that SA8 is a unique gene family in Actiniarians, evolving through a variety of structural changes including tandem and proximal gene duplication and an inversion event that together allowed SA8 to be recruited into the venom of T. stephensoni.

Venom comparisons of endemic and micro-endemic speckled rattlesnakes Crotalus mitchellii, C. polisi and C. thalassoporus from Baja California Peninsula

A high diversity of rattlesnake species can be found in the Baja California peninsula and the island of the Gulf of California, nevertheless, their venom has been poorly evaluated. The aim of this work was to present the first characterization of endemic Crotalus mitchellii, micro endemic C. polisi and C. thalassoporus venoms. All samples provoke human plasma coagulation showing doses in the rank of 2.3-41.0 μg and also produce rapid hydrolysis of the alpha chain of bovine fibrinogen while the beta chain is attacked at larger incubation periods by C. polisi and especially by C. thalassoporus. Phospholipase activity ranging from 23.2 to 173.8 U/mg. The venoms of C. thalassoporus and C. polisi show very high hemorrhagic activity (from 0.03 to 0.31 μg). A total of 130 toxin-related proteins were identified and classified into ten families. Crotalus mitchellii venom was characterized by high abundance of crotoxin-like and other phospholipase proteins (34.5%) and serine proteinases (29.8%). Crotalus polisi showed a similar proportion of metalloproteinases (34%) and serine proteinases (22.8%) components with important contribution of C-type lectins (14.3%) and CRiSP (14.0%) proteins. Venom of C. thalassoporus is dominated by metalloproteases that amount to more than 66% of total toxin proteins. These results provide a foundation for comprehending the biological, ecological and evolutionary significance of venom composition of speckled rattlesnake from the Baja California peninsula.

Micro and macroevolution of sea anemone venom phenotype

Venom is a complex trait with substantial inter- and intraspecific variability resulting from strong selective pressures acting on the expression of many toxic proteins. However, understanding the processes underlying toxin expression dynamics that determine the venom phenotype remains unresolved. By interspecific comparisons we reveal that toxin expression in sea anemones evolves rapidly and that in each species different toxin family dictates the venom phenotype by massive gene duplication events. In-depth analysis of the sea anemone, Nematostella vectensis, revealed striking variation of the dominant toxin (Nv1) diploid copy number across populations (1-24 copies) resulting from independent expansion/contraction events, which generate distinct haplotypes. Nv1 copy number correlates with expression at both the transcript and protein levels with one population having a near-complete loss of Nv1 production. Finally, we establish the dominant toxin hypothesis which incorporates observations in other venomous lineages that animals have convergently evolved a similar strategy in shaping their venom.

Highly Evolvable: Investigating Interspecific and Intraspecific Venom Variation in Taipans (Oxyuranus spp.) and Brown Snakes (Pseudonaja spp.)

Snake venoms are complex mixtures of toxins that differ on interspecific (between species) and intraspecific (within species) levels. Whether venom variation within a group of closely related species is explained by the presence, absence and/or relative abundances of venom toxins remains largely unknown. Taipans (Oxyuranus spp.) and brown snakes (Pseudonaja spp.) represent medically relevant species of snakes across the Australasian region and provide an excellent model clade for studying interspecific and intraspecific venom variation. Using liquid chromatography with ultraviolet and mass spectrometry detection, we analyzed a total of 31 venoms covering all species of this monophyletic clade, including widespread localities. Our results reveal major interspecific and intraspecific venom variation in Oxyuranus and Pseudonaja species, partially corresponding with their geographical regions and phylogenetic relationships. This extensive venom variability is generated by a combination of the absence/presence and differential abundance of venom toxins. Our study highlights that venom systems can be highly dynamical on the interspecific and intraspecific levels and underscores that the rapid toxin evolvability potentially causes major impacts on neglected tropical snakebites.

A current perspective on snake venom composition and constituent protein families

Snake venoms are heterogeneous mixtures of proteins and peptides used for prey subjugation. With modern proteomics there has been a rapid expansion in our knowledge of snake venom composition, resulting in the venom proteomes of 30% of vipers and 17% of elapids being characterised. From the reasonably complete proteomic coverage of front-fanged snake venom composition (179 species-68 species of elapids and 111 species of vipers), the venoms of vipers and elapids contained 42 different protein families, although 18 were only reported in < 5% of snake species. Based on the mean abundance and occurrence of the 42 protein families, they can be classified into 4 dominant, 6 secondary, 14 minor, and 18 rare protein families. The dominant, secondary and minor categories account for 96% on average of a snake's venom composition. The four dominant protein families are: phospholipase A₂ (PLA₂), snake venom metalloprotease (SVMP), three-finger toxins (3FTx), and snake venom serine protease (SVSP). The six secondary protein families are: L-amino acid oxidase (LAAO), cysteine-rich secretory protein (CRiSP), C-type lectins (CTL), disintegrins (DIS), kunitz peptides (KUN), and natriuretic peptides (NP). Venom variation occurs at all taxonomic levels, including within populations. The reasons for venom variation are complex, as variation is not always associated with geographical variation in diet. The four dominant protein families appear to be the most important toxin families in human envenomation, being responsible for coagulopathy, neurotoxicity, myotoxicity and cytotoxicity. Proteomic techniques can be used to investigate the toxicological profile of a snake venom and hence identify key protein families for antivenom immunorecognition.

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.

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.

Identification of Daboia siamensis venome using integrated multi-omics data

Snakebite, classified by World Health Organization as a neglected tropical disease, causes more than 100,000 deaths and 2 million injuries per year. Currently, available antivenoms do not bind with strong specificity to target toxins, which means that severe complications can still occur despite treatment. Moreover, the cost of antivenom is expensive. Knowledge of venom compositions is fundamental for producing a specific antivenom that has high effectiveness, low side effects, and ease of manufacture. With advances in mass spectrometry techniques, venom proteomes can now be analyzed in great depth at high efficiency. However, these techniques require genomic and transcriptomic data for interpreting mass spectrometry data. This study aims to establish and incorporate genomics, transcriptomics, and proteomics data to study venomics of a venomous snake, Daboia siamensis. Multiple proteins that have not been reported as venom components of this snake such as hyaluronidase-1, phospholipase B, and waprin were discovered. Thus, multi-omics data are advantageous for venomics studies. These findings will be valuable not only for antivenom production but also for the development of novel therapeutics.

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.

Contextual Constraints: Dynamic Evolution of Snake Venom Phospholipase A2

Venom is a dynamic trait that has contributed to the success of numerous organismal lineages. Predominantly composed of proteins, these complex cocktails are deployed for predation and/or self-defence. Many non-toxic physiological proteins have been convergently and recurrently recruited by venomous animals into their toxin arsenal. Phospholipase A₂ (PLA₂) is one such protein and features in the venoms of many organisms across the animal kingdom, including snakes of the families Elapidae and Viperidae. Understanding the evolutionary history of this superfamily would therefore provide insight into the origin and diversification of venom toxins and the evolution of novelty more broadly. The literature is replete with studies that have identified diversifying selection as the sole influence on PLA₂ evolution. However, these studies have largely neglected the structural/functional constraints on PLA₂s, and the ecology and evolutionary histories of the diverse snake lineages that produce them. By considering these crucial factors and employing evolutionary analyses integrated with a schema for the classification of PLA₂s, we uncovered lineage-specific differences in selection regimes. Thus, our work provides novel insights into the evolution of this major snake venom toxin superfamily and underscores the importance of considering the influence of evolutionary and ecological contexts on molecular evolution.

De Novo Genome Assembly Highlights the Role of Lineage-Specific Duplications in the Evolution of Venom in Fea’s Viper

Despite the medical significance to humans and important ecological roles filled by vipers, few high-quality genomic resources exist for these snakes outside of a few genera of pitvipers. Here we sequence, assemble, and annotate the genome of Fea’s Viper (Azemiops feae). This taxon is distributed in East Asia and belongs to a monotypic subfamily, sister to the pitvipers. The newly sequenced genome resulted in a 1.56 Gb assembly, a contig N50 of 1.59 Mb, with 97.6% of the genome assembly in contigs >50 Kb, and a BUSCO completeness of 92.4%. We found that A. feae venom is primarily composed of phospholipase A₂ (PLA₂) proteins expressed by genes that likely arose from lineage-specific PLA₂ gene duplications. Additionally, we show that renin, an enzyme associated with blood pressure regulation in mammals and known from the venoms of two viper species including A. feae, is expressed in the venom gland at comparative levels to known toxins and is present in the venom proteome. The cooption of this gene as a toxin may be more widespread in viperids than currently known. To investigate the historical population demographics of A. feae, we performed coalescent-based analyses and determined that the effective population size has remained stable over the last 100 kyr. This suggests Quaternary glacial cycles likely had minimal influence on the demographic history of A. feae. This newly assembled genome will be an important resource for studying the genomic basis of phenotypic evolution and understanding the diversification of venom toxin gene families.

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.

Modern venomics-Current insights, novel methods, and future perspectives in biological and applied animal venom research

Venoms have evolved >100 times in all major animal groups, and their components, known as toxins, have been fine-tuned over millions of years into highly effective biochemical weapons. There are many outstanding questions on the evolution of toxin arsenals, such as how venom genes originate, how venom contributes to the fitness of venomous species, and which modifications at the genomic, transcriptomic, and protein level drive their evolution. These questions have received particularly little attention outside of snakes, cone snails, spiders, and scorpions. Venom compounds have further become a source of inspiration for translational research using their diverse bioactivities for various applications. We highlight here recent advances and new strategies in modern venomics and discuss how recent technological innovations and multi-omic methods dramatically improve research on venomous animals. The study of genomes and their modifications through CRISPR and knockdown technologies will increase our understanding of how toxins evolve and which functions they have in the different ontogenetic stages during the development of venomous animals. Mass spectrometry imaging combined with spatial transcriptomics, in situ hybridization techniques, and modern computer tomography gives us further insights into the spatial distribution of toxins in the venom system and the function of the venom apparatus. All these evolutionary and biological insights contribute to more efficiently identify venom compounds, which can then be synthesized or produced in adapted expression systems to test their bioactivity. Finally, we critically discuss recent agrochemical, pharmaceutical, therapeutic, and diagnostic (so-called translational) aspects of venoms from which humans benefit.

Evolution, expression patterns and distribution of novel ribbon worm predatory and defensive toxins

Ribbon worms are active predators that use an eversible proboscis to inject venom into their prey and defend themselves with toxic epidermal secretions. Previous work on nemertean venom has largely focused on just a few species and has not investigated the different predatory and defensive secretions in detail. Consequently, our understanding of the composition and evolution of ribbon worm venoms is still very limited. Here, we present a comparative study of nemertean venom combining RNA-seq differential gene expression analyses of venom-producing tissues, tandem mass spectrometry-based proteomics of toxic secretions, and mass spectrometry imaging of proboscis sections, to shed light onto the composition and evolution of predatory and defensive toxic secretions in Antarctonemertes valida. Our analyses reveal a wide diversity of putative defensive and predatory toxins with tissue-specific gene expression patterns and restricted distributions to the mucus and proboscis proteomes respectively, suggesting that ribbon worms produce distinct toxin cocktails for predation and defense. Our results also highlight the presence of numerous lineage-specific toxins, indicating that venom evolution is highly divergent across nemerteans, producing toxin cocktails that might be finely tuned to subdue different prey. Our data also suggest that the hoplonemertean proboscis is a highly specialized predatory organ that seems to be involved in a variety of biological functions besides predation, including secretion and sensory perception. Overall, our results advance our knowledge into the diversity and evolution of nemertean venoms and highlight the importance of combining different types of data to characterize toxin composition in understudied venomous organisms.

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.

Snake Venomics: Fundamentals, Recent Updates, and a Look to the Next Decade

Venomic research, powered by techniques adapted from proteomics, transcriptomics, and genomics, seeks to unravel the diversity and complexity of venom through which knowledge can be applied in the treatment of envenoming, biodiscovery, and conservation. Snake venom proteomics is most extensively studied, but the methods varied widely, creating a massive amount of information which complicates data comparison and interpretation. Advancement in mass spectrometry technology, accompanied by growing databases and sophisticated bioinformatic tools, has overcome earlier limitations of protein identification. The progress, however, remains challenged by limited accessibility to samples, non-standardized quantitative methods, and biased interpretation of -omic data. Next-generation sequencing (NGS) technologies enable high-throughput venom-gland transcriptomics and genomics, complementing venom proteomics by providing deeper insights into the structural diversity, differential expression, regulation and functional interaction of the toxin genes. Venomic tissue sampling is, however, difficult due to strict regulations on wildlife use and transfer of biological materials in some countries. Limited resources for techniques and funding are among other pertinent issues that impede the progress of venomics, particularly in less developed regions and for neglected species. Genuine collaboration between international researchers, due recognition of regional experts by global organizations (e.g., WHO), and improved distribution of research support, should be embraced.

Differences in PLA2 Constitution Distinguish the Venom of Two Endemic Brazilian Mountain Lanceheads, Bothrops cotiara and Bothrops fonsecai

Interspecific differences in snake venom compositions can result from distinct regulatory mechanisms acting in each species. However, comparative analyses focusing on identifying regulatory elements and patterns that led to distinct venom composition are still scarce. Among venomous snakes, Bothrops cotiara and Bothrops fonsecai represent ideal models to complement our understanding of the regulatory mechanisms of venom production. These recently diverged species share a similar specialized diet, habitat, and natural history, but each presents a distinct venom phenotype. Here, we integrated data from the venom gland transcriptome and miRNome and the venom proteome of B. fonsecai and B. cotiara to better understand the regulatory mechanisms that may be acting to produce differing venom compositions. We detected not only the presence of similar toxin isoforms in both species but also distinct expression profiles of phospholipases A2 (PLA2) and some snake venom metalloproteinases (SVMPs) and snake venom serine proteinases (SVSPs) isoforms. We found evidence of modular expression regulation of several toxin isoforms implicated in venom divergence and observed correlated expression of several transcription factors. We did not find strong evidence for miRNAs shaping interspecific divergence of the venom phenotypes, but we identified a subset of toxin isoforms whose final expression may be fine-tuned by specific miRNAs. Sequence analysis on orthologous toxins showed a high rate of substitutions between PLA2s, which indicates that these toxins may be under strong positive selection or represent paralogous toxins in these species. Our results support other recent studies in suggesting that gene regulation is a principal mode of venom evolution across recent timescales, especially among species with conserved ecotypes.

Characterization of the venom and external morphology of a natural hybrid between Crotalus atrox and Crotalus mictlantecuhtli

Here we report, for the first time, a natural hybrid between Crotalus atrox and C. mictlantecuhtli based on intermediate characteristics of the external morphology and venom. Morphologically, the individual had characteristics of both parent species. The hybrid's venom exhibited an intermediate composition including the presence of crotoxin which has never been documented in C. atrox but is well documented in C. mictlantecuhtli. The hybrid's venom was highly toxic and showed an intermediate proteolytic activity between the parental species. The two Mexican antivenoms were able to neutralize the hybrid's venom's lethality.

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.

Varying Intensities of Introgression Obscure Incipient Venom-Associated Speciation in the Timber Rattlesnake (Crotalus horridus)

Ecologically divergent selection can lead to the evolution of reproductive isolation through the process of ecological speciation, but the balance of responsible evolutionary forces is often obscured by an inadequate assessment of demographic history and the genetics of traits under selection. Snake venoms have emerged as a system for studying the genetic basis of adaptation because of their genetic tractability and contributions to fitness, and speciation in venomous snakes can be associated with ecological diversification such as dietary shifts and corresponding venom changes. Here, we explored the neurotoxic (type A)–hemotoxic (type B) venom dichotomy and the potential for ecological speciation among Timber Rattlesnake (Crotalus horridus) populations. Previous work identified the genetic basis of this phenotypic difference, enabling us to characterize the roles geography, history, ecology, selection, and chance play in determining when and why new species emerge or are absorbed. We identified significant genetic, proteomic, morphological, and ecological/environmental differences at smaller spatial scales, suggestive of incipient ecological speciation between type A and type B C. horridus. Range-wide analyses, however, rejected the reciprocal monophyly of venom type, indicative of varying intensities of introgression and a lack of reproductive isolation across the range. Given that we have now established the phenotypic distributions and ecological niche models of type A and B populations, genome-wide data are needed and capable of determining whether type A and type B C. horridus represent distinct, reproductively isolated lineages due to incipient ecological speciation or differentiated populations within a single species.

Integration of transcriptomic and proteomic approaches for snake venom profiling

INTRODUCTION

Snake venoms contain many protein and peptide isoforms with high levels of sequence variation, even within a single species.

AREAS COVERED

In this review, we highlight several examples, from both published and unpublished work in our lab, demonstrating how a combined venom gland transcriptome and proteome methodology allows for comprehensive characterization of venoms, including those from understudied rear-fanged snake species, and we provide recommendations for using these approaches.

EXPERT OPINION

When characterizing venoms, peptide mass fingerprinting using databases built predominately from protein sequences originating from model organisms can be disadvantageous, especially when the intention is to document protein diversity. Therefore, the use of species-specific venom gland transcriptomes corrects for the absence of these unique peptide sequences in databases. The integration of transcriptomics and proteomics improves the accuracy of either approach alone for venom profiling.

Intraspecific venom variation of Mexican West Coast Rattlesnakes (Crotalus basiliscus) and its implications for antivenom production

Intraspecific variation in snake venoms has been widely documented worldwide. However, there are few studies on this subject in Mexico. Venom characterization studies provide important data used to predict clinical syndromes, to evaluate the efficacy of antivenoms and, in some cases, to improve immunogenic mixtures in the production of antivenoms. In the present work, we evaluated the intraspecific venom variation of Crotalus basiliscus, a rattlesnake of medical importance and whose venom is used in the immunization of horses to produce one of the Mexican antivenoms. Our results demonstrate that there is variation in biological and biochemical activities among adult venoms and that there is an ontogenetic change from juvenile to adult venoms. Juvenile venoms were more lethal and had higher percentages of crotamine and crotoxin, while adult venoms had higher percentages of snake venom metalloproteases (SVMPs). Additionally, we documented crotoxin-like PLA₂ variation in which specimens from Zacatecas, Sinaloa and Michoacán (except 1) lacked the neurotoxin, while the rest of the venoms had it. Finally, we evaluated the efficacy of three lots of Birmex antivenom and all three were able to neutralize the lethality of four representative venoms but were not able to neutralize crotamine. We also observed significant differences in the LD₅₀ values neutralized per vial among the different lots. Based on these results, we recommend including venoms containing crotamine in the production of antivenom for a better immunogenic mixture and to improve the homogeneity of lots.

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.

Tentacle Morphological Variation Coincides with Differential Expression of Toxins in Sea Anemones

Phylum Cnidaria is an ancient venomous group defined by the presence of cnidae, specialised organelles that serve as venom delivery systems. The distribution of cnidae across the body plan is linked to regionalisation of venom production, with tissue-specific venom composition observed in multiple actiniarian species. In this study, we assess whether morphological variants of tentacles are associated with distinct toxin expression profiles and investigate the functional significance of specialised tentacular structures. Using five sea anemone species, we analysed differential expression of toxin-like transcripts and found that expression levels differ significantly across tentacular structures when substantial morphological variation is present. Therefore, the differential expression of toxin genes is associated with morphological variation of tentacular structures in a tissue-specific manner. Furthermore, the unique toxin profile of spherical tentacular structures in families Aliciidae and Thalassianthidae indicate that vesicles and nematospheres may function to protect branched structures that host a large number of photosynthetic symbionts. Thus, hosting zooxanthellae may account for the tentacle-specific toxin expression profiles observed in the current study. Overall, specialised tentacular structures serve unique ecological roles and, in order to fulfil their functions, they possess distinct venom cocktails.

Gallic acid anti-myotoxic activity and mechanism of action, a snake venom phospholipase A2 toxin inhibitor, isolated from the medicinal plant Anacardium humile

Snakebite envenoming is the cause of an ongoing health crisis in several regions of the world, particularly in tropical and neotropical countries. This scenario creates an urgent necessity for new practical solutions to address the limitations of current therapies. The current study investigated the isolation, phytochemical characterization, and myotoxicity inhibition mechanism of gallic acid (GA), a myotoxin inhibitor obtained from Anacardium humile. The identification and isolation of GA was achieved by employing analytical chromatographic separation, which exhibited a compound with retention time and nuclear magnetic resonance spectra compatible with GA's commercial standard and data from the literature. GA alone was able to inhibit the myotoxic activity induced by the crude venom of Bothrops jararacussu and its two main myotoxins, BthTX-I and BthTX-II. Circular dichroism (CD), fluorescence spectroscopy (FS), dynamic light scattering (DLS), and interaction studies by molecular docking suggested that GA forms a complex with BthTX-I and II. Surface plasmon resonance (SPR) kinetics assays showed that GA has a high affinity for BthTX-I with a K_D of 9.146 × 10^-7 M. Taken together, the two-state reaction mode of GA binding to BthTX-I, and CD, FS and DLS assays, suggest that GA is able to induce oligomerization and secondary structure changes for BthTX-I and -II. GA and other tannins have been shown to be effective inhibitors of snake venoms' toxic effects, and herein we demonstrated GA's ability to bind to and inhibit a snake venom PLA₂, thus proposing a new mechanism of PLA₂ inhibition, and presenting more evidence of GA's potential as an antivenom compound.

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.

Identification of Novel Toxin Genes from the Stinging Nettle Caterpillar Parasa lepida (Cramer, 1799): Insights into the Evolution of Lepidoptera Toxins

Simple Summary

Many caterpillar species can produce toxins that cause harmful reactions to humans, varying from mild irritation to death. Currently, there is very limited knowledge about caterpillar toxin diversity, because only a few species have been investigated. We used the transcriptome technique to identify candidate toxin genes from the nettle caterpillar Parasa lepida (Cramer, 1799). It is a common pest of oil palm, coconut, and mango in South and South-East Asia, which can cause severe pain and allergic responses to those in contact with them. We reported 168 candidate toxin genes. Most of them are members of the toxin genes families commonly recruited in animal venoms such as serine protease and serine protease inhibitors. However, we identified 21 novel genes encoding knottin-like peptides expressed at a high level in the transcriptome. Their predicted 3D structures are similar to neurotoxins in scorpion and tarantula. Our study suggests that P. lepida venom contains diverse toxin proteins that potentially cause allergic reactions and pain. This study sheds light on the hidden diversity of toxin proteins in caterpillar lineage, which could be future fruitful new drug sources.

Abstract

Many animal species can produce venom for defense, predation, and competition. The venom usually contains diverse peptide and protein toxins, including neurotoxins, proteolytic enzymes, protease inhibitors, and allergens. Some drugs for cancer, neurological disorders, and analgesics were developed based on animal toxin structures and functions. Several caterpillar species possess venoms that cause varying effects on humans both locally and systemically. However, toxins from only a few species have been investigated, limiting the full understanding of the Lepidoptera toxin diversity and evolution. We used the RNA-seq technique to identify toxin genes from the stinging nettle caterpillar, Parasa lepida (Cramer, 1799). We constructed a transcriptome from caterpillar urticating hairs and reported 34,968 unique transcripts. Using our toxin gene annotation pipeline, we identified 168 candidate toxin genes, including protease inhibitors, proteolytic enzymes, and allergens. The 21 P. lepida novel Knottin-like peptides, which do not show sequence similarity to any known peptide, have predicted 3D structures similar to tarantula, scorpion, and cone snail neurotoxins. We highlighted the importance of convergent evolution in the Lepidoptera toxin evolution and the possible mechanisms. This study opens a new path to understanding the hidden diversity of Lepidoptera toxins, which could be a fruitful source for developing new drugs.

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.

The absence of thrombin-like activity in Bothrops erythromelas venom is due to the deletion of the snake venom thrombin-like enzyme gene

Snake venom thrombin-like enzymes (SVTLEs) are serine proteinases that clot fibrinogen. SVTLEs are distributed mainly in venoms from snakes of the Viperidae family, comprising venomous pit viper snakes. Bothrops snakes are distributed throughout Central and South American and are responsible for most venomous snakebites. Most Bothrops snakes display thrombin-like activity in their venoms, but it has been shown that some species do not present it. In this work, to understand SVTLE polymorphism in Bothrops snake venoms, we studied individual samples from two species of medical importance in Brazil: Bothrops jararaca, distributed in Southeastern Brazil, which displays coagulant activity on plasma and fibrinogen, and Bothrops erythromelas, found in Northeastern Brazil, which lacks direct fibrinogen coagulant activity but shows plasma coagulant activity. We tested the coagulant activity of venoms and the presence of SVTLE genes by a PCR approach. The SVTLE gene structure in B. jararaca is similar to the Bothrops atrox snake, comprising five exons. We could not amplify SVTLE sequences from B. erythromelas DNA, except for a partial pseudogene. These genes underwent a positive selection in some sites, leading to an amino acid sequence diversification, mostly in exon 2. The phylogenetic tree constructed using SVTLE coding sequences confirms that they are related to the chymotrypsin/kallikrein family. Interestingly, we found a B. jararaca specimen whose venom lacked thrombin-like activity, and its gene sequence was a pseudogene with SVTLE structure, presenting nonsense and frameshift mutations. Our results indicate an association of the lack of thrombin-like activity in B. jararaca and B. erythromelas venoms with mutations and deletions of snake venom thrombin-like enzyme genes.