Snake venoms are integrated systems, but abundant venom proteins evolve more rapidly

A Russian Doll of Resistance: Nested Gains and Losses of Venom Immunity in Varanid Lizards

The interplay between predator and prey has catalyzed the evolution of venom systems, with predators honing their venoms in response to the evolving resistance of prey. A previous study showed that the African varanid species Varanus exanthematicus has heightened resistance to snake venoms compared to the Australian species V. giganteus, V. komodoensis, and V. mertensi, likely due to increased predation by sympatric venomous snakes on V. exanthematicus. To understand venom resistance among varanid lizards, we analyzed the receptor site targeted by venoms in 27 varanid lizards, including 25 Australian varanids. The results indicate an active evolutionary arms race between Australian varanid lizards and sympatric neurotoxic elapid snakes. Large species preying on venomous snakes exhibit inherited neurotoxin resistance, a trait potentially linked to their predatory habits. Consistent with the 'use it or lose it' aspect of venom resistance, this trait was secondarily reduced in two lineages that had convergently evolved gigantism (V. giganteus and the V. komodoensis/V. varius clade), suggestive of increased predatory success accompanying extreme size and also increased mechanical protection against envenomation due to larger scale osteoderms. Resistance was completely lost in the mangrove monitor V. indicus, consistent with venomous snakes not being common in their arboreal and aquatic niche. Conversely, dwarf varanids demonstrate a secondary loss at the base of the clade, with resistance subsequently re-evolving in the burrowing V. acanthurus/V. storri clade, suggesting an ongoing battle with neurotoxic predators. Intriguingly, within the V. acanthurus/V. storri clade, resistance was lost again in V. kingorum, which is morphologically and ecologically distinct from other members of this clade. Resistance was also re-evolved in V. glebopalma which is terrestrial in contrast to the arboreal/cliff dwelling niches occupied by the other members of its clade (V. glebopalma, V. mitchelli, V. scalaris, V. tristis). This 'Russian doll' pattern of venom resistance underscores the dynamic interaction between dwarf varanids and Australian neurotoxic elapid snakes. Our research, which included testing Acanthophis (death adder) venoms against varanid receptors as models for alpha-neurotoxic interactions, uncovered a fascinating instance of the Red Queen Hypothesis: some death adders have developed more potent toxins specifically targeting resistant varanids, a clear sign of the relentless predator-prey arms race. These results offer new insight into the complex dynamics of venom resistance and highlight the intricate ecological interactions that shape the natural world.

Varying Modes of Selection Among Toxin Families in the Venoms of the Giant Desert Hairy Scorpions (Hadrurus)

Venoms are primarily believed to evolve under strong diversifying selection resulting from persistent coevolution between predator and prey. Recent research has challenged this hypothesis, proposing that venoms from younger venomous lineages (e.g., snakes and cone snails) are governed predominantly by diversifying selection, while venoms from older venomous lineages (e.g., centipedes, scorpions, and spiders) are under stronger purifying selection. However, most research in older lineages has tested selection at more diverse phylogenetic scales. Although these tests are important for evaluating broad macroevolutionary trends underlying venom evolution, they are less equipped to detect species-level evolutionary trends, which likely have large impacts on venom variation seen at more diverse phylogenetic scales. To test for selection among closely related species from an older venomous lineage, we generated high-throughput venom-gland transcriptomes and venom proteomes for four populations of Giant Desert Hairy Scorpions (Hadrurus), including three Hadrurus arizonensis populations and one Hadrurus spadix population. We detected significant episodic and pervasive diversifying selection across a highly abundant toxin family that likely has a major role in venom function ([Formula: see text]KTxs), providing a contrast to the stronger purifying selection identified from other studies on scorpion venoms. Conversely, we detected weak episodic diversifying and/or stronger purifying selection in four toxin families (non-disulfide bridged peptides, phospholipase A2s, scorpine-like antimicrobial peptides, and serine proteases), most of which were less abundant and likely have ancillary functional roles. Finally, although we detected several major toxin families at disproportionate transcriptomic and/or proteomic abundances, we did not identify significant sex-based variation in Hadrurus venoms.

Genome assembly and annotation of the Brown-Spotted Pit viper Protobothrops mucrosquamatus

The Brown-Spotted Pit viper (Protobothrops mucrosquamatus), also known as the Chinese habu, is a widespread and highly venomous snake distributed from Northeastern India to Eastern China. Genomics research can contribute to our understanding of venom components and natural selection in vipers. Here, we collected, sequenced and assembled the genome of a male P. mucrosquamatus individual from China. We generated a highly continuous reference genome, with a length of 1.53 Gb and 41.18% of repeat elements content. Using this genome, we identified 24,799 genes, 97.97% of which could be annotated. We verified the validity of our genome assembly and annotation process by generating a phylogenetic tree based on the nuclear genome single-copy genes of six other reptile species. The results of our research will contribute to future studies on Protobothrops biology and the genetic basis of snake venom.

The best of both worlds? Rattlesnake hybrid zones generate complex combinations of divergent venom phenotypes that retain high toxicity

Studying the consequences of hybridization between closely related species with divergent traits can reveal patterns of evolution that shape and maintain extreme trophic adaptations. Snake venoms are an excellent model system for examining the evolutionary and ecological patterns that underlie highly selected polymorphic traits. Here we investigate hybrid venom phenotypes that result from natural introgression between two rattlesnake species that express highly divergent venom phenotypes: Crotalus o. concolor and C. v. viridis. Though not yet documented, interbreeding between these species may lead to novel venom phenotypes with unique activities that break the typical trends of venom composition in rattlesnakes. The characteristics of these unusual phenotypes could unveil the roles of introgression in maintaining patterns of venom composition and variation, including the near ubiquitous dichotomy between neurotoxic or degradative venoms observed across rattlesnakes. We use RADseq data to infer patterns of gene flow and hybrid ancestry between these diverged lineages and link these genetic data with analyses of venom composition, biological activity, and whole animal model toxicity tests to understand the impacts of introgression on venom composition. We find that introgressed populations express admixed venom phenotypes that do not sacrifice biological activity (lethal toxicity) or overall abundance of dominant toxins compared to parental venoms. These hybridized venoms therefore do not represent a trade-off in functionality between the typical phenotypic extremes but instead represent a unique combination of characters whose expression appears limited to the hybrid zone.

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.

Phylogeny-Related Variations in Venomics: A Test in a Subset of Habu Snakes (Protobothrops)

We conducted a comparative analysis to unveil the divergence among venoms from a subset of Old World habu snakes (Protobothrops) in terms of venomic profiles and toxicological and enzymatic activities. A total of 14 protein families were identified in the venoms from these habu snakes, and 11 of them were shared among these venoms. The venoms of five adult habu snakes were overwhelmingly dominated by SVMP (32.56 ± 13.94%), PLA₂ (22.93 ± 9.26%), and SVSP (16.27 ± 4.79%), with a total abundance of over 65%, while the subadult P. mangshanensis had an extremely low abundance of PLA₂ (1.23%) but a high abundance of CTL (51.47%), followed by SVMP (22.06%) and SVSP (10.90%). Apparent interspecific variations in lethality and enzymatic activities were also explored in habu snake venoms, but no variations in myotoxicity were found. Except for SVSP, the resemblance of the relatives within Protobothrops in other venom traits was estimated to deviate from Brownian motion evolution based on phylogenetic signals. A comparative analysis further validated that the degree of covariation between phylogeny and venom variation is evolutionarily labile and varies among clades of closely related snakes. Our findings indicate a high level of interspecific variation in the venom proteomes of habu snakes, both in the presence or absence and the relative abundance of venom protein families, and that these venoms might have evolved under a combination of adaptive and neutral mechanisms.

Diversity Begets Diversity When Diet Drives Snake Venom Evolution, but Evenness Rather Than Richness Is What Counts

Snake venoms are primarily used to subjugate prey, and consequently, their evolution has been shown to be predominantly driven by diet-related selection pressure. Venoms tend to be more lethal to prey than non-prey species (except in cases of toxin resistance), prey-specific toxins have been identified, and preliminary work has demonstrated an association between the diversity of diet classes and that of toxicological activities of whole venom. However, venoms are complex mixtures of many toxins, and it remains unclear how toxin diversity is driven by diet. Prey-specific toxins do not encompass the molecular diversity of venoms, and whole venom effects could be driven by one, few, or all components, so the link between diet and venom diversity remains minimally understood. Here, we collated a database of venom composition and diet records and used a combination of phylogenetic comparative methods and two quantitative diversity indices to investigate whether and how diet diversity relates to the toxin diversity of snake venoms. We reveal that venom diversity is negatively related to diet diversity using Shannon's index but positively related using Simpson's index. Since Shannon's index predominantly considers the number of prey/toxins, whereas Simpson's index more strongly reflects evenness, we provide insights into how the diet-venom diversity link is driven. Specifically, species with low diet diversity tend to have venoms dominated by a few abundant (possibly specialised) toxin families, whereas species with diverse diets tend to 'hedge their bets' by having venoms with a more even composition of different toxin classes.

Full sequencing and comparison of five venom metalloproteases of Trimeresurus gracilis: The PI-enzyme is most similar to okinalysin but the PIII-enzyme is most similar to Crotalus venom enzymes

The cDNAs encoding the Zn⁺²-metalloproteases (SVMPs) of Trimeresurus gracilis (abbreviated as Tgc), a pitviper endemic to Taiwan, were cloned from venom glands and sequenced. The amino-acid sequences of five novel SVMPs, including one P-III, three P-II and one P-I class enzymes, were thus deduced and subjected to BLAST-analyses. The P-III enzyme (designated as Tgc-PIII) is structurally most similar to the PIII-SVMPs of New World pitvipers but not similar to the PIII-SVMP of Ovophis okinavensis. Sequence-similarity analysis of 22 homologous PIII-SVMPs reveal three major structural subtypes of the pitviper PIII-SVMPs, which possibly have different substrate specificities. In addition, Tgc-PIII and the PI-class SVMP (named Tgc-MP) were isolated from the venom and verified by mass spectrometry. All the three deduced sequences of PII-SVMPs (Tgc-PIIs) contain an abnormal Zn⁺²-binding-site in their catalytic-domain, and an identical "long-disintegrin" domain. The predicted 85-residues disintegrin, gracilisin, bears high similarities to some long-disintegrins of the New-World pitvipers and salmosin3. By BLAST search and comparison, Tgc-MP is 96% similar to okinalysin, the hemorrhagic PI-SVMP of O. okinavensis, rather than any other PI-SVMPs in the databanks. Our results confirm the fast evolution of Tgc-SVMPs as well as their structural similarities to different SVMP-classes of the New-World pitvipers and of O. okinavensis, respectively. The implications of our findings are discussed along with our previous sequence comparisons of venom phospholipases A₂ and ten venom serine proteases of Tgc.

The Biological Hierarchy, Time, and Temporal 'Omics in Evolutionary Biology: A Perspective

Sequencing data-genomics, transcriptomics, epigenomics, proteomics, and metabolomics-have revolutionized biological research, enabling a more detailed study of processes, ranging from subcellular to evolutionary, that drive biological organization. These processes, collectively, are responsible for generating patterns of phenotypic variation and can operate over dramatically different timescales (milliseconds to billions of years). While researchers often study phenotypic variation at specific levels of biological organization to isolate processes operating at that particular scale, the varying types of sequence data, or 'omics, can also provide complementary inferences to link molecular and phenotypic variation to produce an integrated view of evolutionary biology, ranging from molecular pathways to speciation. We briefly describe how 'omics has been used across biological levels and then demonstrate the utility of integrating different types of sequencing data across multiple biological levels within the same study to better understand biological phenomena. However, single-time-point studies cannot evaluate the temporal dynamics of these biological processes. Therefore, we put forward temporal 'omics as a framework that can better enable researchers to study the temporal dynamics of target processes. Temporal 'omics is not infallible, as the temporal sampling regime directly impacts inferential ability. Thus, we also discuss the role the temporal sampling regime plays in deriving inferences about the environmental conditions driving biological processes and provide examples that demonstrate the impact of the sampling regime on biological inference. Finally, we forecast the future of temporal 'omics by highlighting current methodological advancements that will enable temporal 'omics to be extended across species and timescales. We extend this discussion to using temporal multi-omics to integrate across the biological hierarchy to evaluate and link the temporal dynamics of processes that generate phenotypic variation.

Biological and Medical Aspects Related to South American Rattlesnake Crotalus durissus (Linnaeus, 1758): A View from Colombia

In Colombia, South America, there is a subspecies of the South American rattlesnake Crotalus durissus, C. d. cumanensis, a snake of the Viperidae family, whose presence has been reduced due to the destruction of its habitat. It is an enigmatic snake from the group of pit vipers, venomous, with large articulated front fangs, special designs on its body, and a characteristic rattle on its tail. Unlike in Brazil, the occurrence of human envenomation by C. durisus in Colombia is very rare and contributes to less than 1% of envenomation caused by snakes. Its venom is a complex cocktail of proteins with different biological effects, which evolved with the purpose of paralyzing the prey, killing it, and starting its digestive process, as well as having defense functions. When its venom is injected into humans as the result of a bite, the victim presents with both local tissue damage and with systemic involvement, including a diverse degree of neurotoxic, myotoxic, nephrotoxic, and coagulopathic effects, among others. Its biological effects are being studied for use in human health, including the possible development of analgesic, muscle relaxant, anti-inflammatory, immunosuppressive, anti-infection, and antineoplastic drugs. Several groups of researchers in Brazil are very active in their contributions in this regard. In this work, a review is made of the most relevant biological and medical aspects related to the South American rattlesnake and of what may be of importance for a better understanding of the snake C. d. cumanensis, present in Colombia and Venezuela.

The Fast and the Furriest: Investigating the Rate of Selection on Mammalian Toxins

The evolution of venom and the selection pressures that act on toxins have been increasingly researched within toxinology in the last two decades, in part due to the exceptionally high rates of diversifying selection observed in animal toxins. In 2015, Sungar and Moran proposed the 'two-speed' model of toxin evolution linking evolutionary age of a group to the rates of selection acting on toxins but due to a lack of data, mammals were not included as less than 30 species of venomous mammal have been recorded, represented by elusive species which produce small amounts of venom. Due to advances in genomics and transcriptomics, the availability of toxin sequences from venomous mammals has been increasing. Using branch- and site-specific selection models, we present the rates of both episodic and pervasive selection acting upon venomous mammal toxins as a group for the first time. We identified seven toxin groups present within venomous mammals, representing Chiroptera, Eulipotyphla and Monotremata: KLK1, Plasminogen Activator, Desmallipins, PACAP, CRiSP, Kunitz Domain One and Kunitz Domain Two. All but one group (KLK1) was identified by our results to be evolving under both episodic and pervasive diversifying selection with four toxin groups having sites that were implicated in the fitness of the animal by TreeSAAP (Selection on Amino Acid Properties). Our results suggest that venomous mammal ecology, behaviour or genomic evolution are the main drivers of selection, although evolutionary age may still be a factor. Our conclusion from these results indicates that mammalian toxins are following the two-speed model of selection, evolving predominately under diversifying selection, fitting in with other younger venomous taxa like snakes and cone snails-with high amounts of accumulating mutations, leading to more novel adaptions in their toxins.

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.

Hybridization between Crotalus aquilus and Crotalus polystictus Species: A Comparison of Their Venom Toxicity and Enzymatic Activities

Hybridization is defined as the interbreeding of individuals from two populations distinguishable by one or more heritable characteristics. Snake hybridization represents an interesting opportunity to analyze variability and how genetics affect the venom components between parents and hybrids. Snake venoms exhibit a high degree of variability related to biological and biogeographical factors. The aim of this work is to analyze the protein patterns and enzymatic activity of some of the main hemotoxic enzymes in snake venoms, such as serine proteases (trypsin-like, chymotrypsin-like, and elastase-like), metalloproteases, hyaluronidases, and phospholipase A₂. The lethal dose of 50 (LD₅₀) of venom from the Crotalus aquilus (Cabf) and Crotalus polystictus (Cpbm) parents and their hybrids in captivity was determined, and phenetic analysis is also conducted, which showed a high similarity between the hybrids and C. polystictus. The protein banding patterns and enzymatic activity analyze by zymography resulted in a combination of proteins from the parental venoms in the hybrids, with variability among them. In some cases, the enzymatic activity is higher in the hybrids with a lower LD₅₀ than in the parents, indicating higher toxicity. These data show the variability among snake venoms and suggest that hybridization is an important factor in changes in protein concentration, peptide variability, and enzymatic activity that affect toxicity and lethality.

Active Expression of Genes for Protein Modification Enzymes in Habu Venom Glands

Genes encoding snake venom toxins have been studied extensively. However, genes involved in the modification and functioning of venom proteins are little known. Protobothrops is a genus of pit vipers, which are venomous and inhabit the Nansei (Southwest) islands of Japan, Taiwan China, Vietnam, Thailand, Myanmar, Nepal, Bhutan, and India. Our previous study decoded the genome of Protobothrops flavoviridis, a species endemic to the Nansei Islands, Japan, and revealed unique evolutionary processes of some venom genes. In this study, we analyzed genes that are highly expressed in venom glands to survey genes for candidate enzymes or chaperone proteins involved in toxin folding and modification. We found that, in addition to genes that encode venom proteins and ribosomal proteins, genes that encode protein disulfide isomerase (PDI) family members (orthologs of human P4HB and PDIA3), Selenoprotein M (SELENOM), and Calreticulin (CALR) are highly expressed in venom glands. Since these enzymes or chaperones are involved in protein modification and potentially possess protein folding functions, we propose that P4HB, SELENOM, CALR, and PDIA3 encode candidate enzymes or chaperones to confer toxic functions upon the venom transcriptome.

Venom gene sequence diversity and expression jointly shape diet adaptation in pitvipers

Understanding the joint roles of protein sequence variation and differential expression during adaptive evolution is a fundamental, yet largely unrealized goal of evolutionary biology. Here, we use phylogenetic path analysis to analyze a comprehensive venom-gland transcriptome dataset spanning three genera of pitvipers to identify the functional genetic basis of a key adaptation (venom complexity) linked to diet breadth (DB). The analysis of gene-family-specific patterns reveals that, for genes encoding two of the most important venom proteins (snake venom metalloproteases and snake venom serine proteases), there are direct, positive relationships between sequence diversity (SD), expression diversity (ED), and increased DB. Further analysis of gene-family diversification for these proteins showed no constraint on how individual lineages achieved toxin gene SD in terms of the patterns of paralog diversification. In contrast, another major venom protein family (PLA₂s) showed no relationship between venom molecular diversity and DB. Additional analyses suggest that other molecular mechanisms—such as higher absolute levels of expression—are responsible for diet adaptation involving these venom proteins. Broadly, our findings argue that functional diversity generated through sequence and expression variations jointly determine adaptation in the key components of pitviper venoms, which mediate complex molecular interactions between the snakes and their prey.

Exploring the five-paced viper (Deinagkistrodon acutus) venom proteome by integrating a combinatorial peptide ligand library approach with shotgun LC-MS/MS

BACKGROUND

Snake venoms are complex mixtures of toxic proteins or peptides encoded by various gene families that function synergistically to incapacitate prey. In the present study, in order to unravel the proteomic repertoire of Deinagkistrodon acutus venom, some trace abundance components were analyzed.

METHODS

Shotgun proteomic approach combined with shotgun nano-LC-ESI-MS/MS were employed to characterize the medically important D. acutus venom, after collected samples were enriched with the combinatorial peptide ligand library (CPLL).

RESULTS

This avenue helped us find some trace components, undetected before, in D. acutus venom. The results indicated that D. acutus venom comprised 84 distinct proteins from 10 toxin families and 12 other proteins. These results are more than twice the number of venom components obtained from previous studies, which were only 29 distinct proteins obtained through RP-HPLC for the venom of the same species. The present results indicated that in D. acutus venom, the most abundant components (66.9%) included metalloproteinases, serine proteinases, and C-type lectin proteins; the medium abundant components (13%) comprised phospholipases A2 (PLA2) and 5'-nucleotidases and nucleases; whereas least abundant components (6%) were aminopeptidases, L-amino acid oxidases (LAAO), neurotoxins and disintegrins; and the trace components. The last were undetected before the use of conventional shotgun proteomics combined with shotgun nano-LC-ESI-MS/MS, such as cysteine-rich secretory proteins Da-CRPa, phospholipases B-like 1, phospholipases B (PLB), nerve growth factors (NGF), glutaminyl-peptide cyclortransferases (QC), and vascular non-inflammatory molecules 2 (VNN2).

CONCLUSION

These findings demonstrated that the CPLL enrichment method worked well in finding the trace toxin proteins in D. acutus venom, in contrast with the previous venomic characterization of D. acutus by conventional LC-MS/MS. In conclusion, this approach combined with the CPLL enrichment was effective for allowing us to explore the hidden D. acutus venomic profile and extended the list of potential venom toxins.

Proteo-Transcriptomic Characterization of Sirex nitobei (Hymenoptera: Siricidae) Venom

The wood-boring woodwasp Sirex nitobei is a native pest in Asia, infecting and weakening the host trees in numerous ecological and commercial coniferous forest plantations. In China, hosts of S. nitobei are diverse, so the pest has spread to several provinces of China, resulting in considerable economic and ecological damage. During female oviposition, S. nitobei venom along with arthrospores of the symbiotic fungus Amylostereum areolatum or A. chaetica is injected into host trees, and the combination of these two biological factors causes the death of xylem host trees. The presence of venom alone causes only the yellowing and wilting of needles. In this study, we constructed the venom gland transcriptome of S. nitobei for the first time and a total of 15,036 unigenes were acquired. From the unigenes, 11,560 ORFs were identified and 537 encoding protein sequences with signal peptides at the N-terminus. Then, we used the venomics approach to characterize the venom composition of female S. nitobei and predicted 1095 proteins by liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis. We focused on seven proteins that were both highly expressed in the venom gland transcriptome and predicted in the crude venom proteome. These seven proteins are laccase-2, laccase-3, a protein belonging to the Kazal family, chitooligosaccharidolytic β-N-acetylglucosaminidase, beta-galactosidase, icarapin-like protein, and waprin-Thr1-like protein. Using quantitative real-time PCR (qRT-PCR), we also proved that the genes related to these seven proteins are specifically expressed in the venom glands. Finally, we revealed the functional role of S. nitobei venom in the physiological response of host trees. It can not only promote the colonization of symbiotic fungus but contribute to the development of eggs and larvae. This study provides a deeper understanding of the molecular mechanism of the woodwasp–pine interaction.

Venom-gland transcriptomic, venomic, and antivenomic profiles of the spine-bellied sea snake (Hydrophis curtus) from the South China Sea

Background

A comprehensive evaluation of the -omic profiles of venom is important for understanding the potential function and evolution of snake venom. Here, we conducted an integrated multi-omics-analysis to unveil the venom-transcriptomic and venomic profiles in a same group of spine-bellied sea snakes (Hydrophis curtus) from the South China Sea, where the snake is a widespread species and might generate regionally-specific venom potentially harmful to human activities. The capacity of two heterologous antivenoms to immunocapture the H. curtus venom was determined for an in-depth evaluation of their rationality in treatment of H. curtus envenomation. In addition, a phylogenetic analysis by maximum likelihood was used to detect the adaptive molecular evolution of full-length toxin-coding unigenes.

Results

A total of 90,909,384 pairs of clean reads were generated via Illumina sequencing from a pooled cDNA library of six specimens, and yielding 148,121 unigenes through de novo assembly. Sequence similarity searching harvested 63,845 valid annotations, including 63,789 non-toxin-coding and 56 toxin-coding unigenes belonging to 22 protein families. Three protein families, three-finger toxins (3-FTx), phospholipase A₂ (PLA₂), and cysteine-rich secretory protein, were detected in the venom proteome. 3-FTx (27.15% in the transcriptome/41.94% in the proteome) and PLA₂ (59.71%/49.36%) were identified as the most abundant families in the venom-gland transcriptome and venom proteome. In addition, 24 unigenes from 11 protein families were shown to have experienced positive selection in their evolutionary history, whereas four were relatively conserved throughout evolution. Commercial Naja atra antivenom exhibited a stronger capacity than Bungarus multicinctus antivenom to immunocapture H. curtus venom components, especially short neurotoxins, with the capacity of both antivenoms to immunocapture short neurotoxins being weaker than that for PLA₂s.

Conclusions

Our study clarified the venom-gland transcriptomic and venomic profiles along with the within-group divergence of a H. curtus population from the South China Sea. Adaptive evolution of most venom components driven by natural selection appeared to occur rapidly during evolutionary history. Notably, the utility of commercial N. atra and B. multicinctus antivenoms against H. curtus toxins was not comprehensive; thus, the development of species-specific antivenom is urgently needed.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07824-7.

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.

Venom Profiling of the Insular Species Bothrops alcatraz: Characterization of Proteome, Glycoproteome, and N-Terminome Using Terminal Amine Isotopic Labeling of Substrates

Bothrops alcatraz, a species endemic to Alcatrazes Islands, is regarded as critically endangered due to its small area of occurrence and the declining quality of its habitat. We recently reported the identification of N-glycans attached to toxins of Bothrops species, showing similar compositions in venoms of the B. jararaca complex (B. jararaca, B. insularis, and B. alcatraz). Here, we characterized B. alcatraz venom using electrophoretic, proteomic, and glycoproteomic approaches. Electrophoresis showed that B. alcatraz venom differs from B. jararaca and B. insularis; however, N-glycan removal revealed similarities between them, indicating that the occupation of N-glycosylation sites contributes to interspecies variability in the B. jararaca complex. Metalloproteinase was the major toxin class identified in the B. alcatraz venom proteome followed by serine proteinase and C-type lectin, and overall, the adult B. alcatraz venom resembles that of B. jararaca juvenile specimens. The comparative glycoproteomic analysis of B. alcatraz venom with B. jararaca and B. insularis indicated that there may be differences in the utilization of N-glycosylation motifs among their different toxin classes. Furthermore, we prospected for the first time the N-terminome of a snake venom using the terminal amine isotopic labeling of substrates (TAILS) approach and report the presence of ∼30% of N-termini corresponding to truncated toxin forms and ∼37% N-terminal sequences blocked by pyroglutamic acid in B. alcatraz venom. These findings underscore a low correlation between venom gland transcriptomes and proteomes and support the view that post-translational processes play a major role in shaping venom phenotypes.

Asymmetrical expression of toxins between the left and right venom glands of an individual prairie rattlesnake (Crotalus viridis viridis)

It is assumed that toxin expression is equivalent between left and right glands of a single snake. In the current study, we report venoms that differ in enzyme functionality and overall composition between the left and right gland of a single snake. The right gland produced venom of comparable composition to venom previously extracted from the same individual; however, the left gland produced venom with overall lower protein content and considerably less enzyme activity.

Toxicity of venom from the mamushi, Gloydius blomhoffii, (Squamata, Crotalinae) to centipedes

Vipers include approximately 300 species and usually feed on vertebrates, but over 30 species of them occasionally eat centipedes. Centipedes have been also known to occur in stomach contents of a Japanese pit viper, mamushi, Gloydius blomhoffii. Toxicity of the venom of mamushi to small mammals has been well studied, but there is no information concerning its toxicity to arthropods. Here, we studied the toxicity of the raw venom to the red-headed centipede, Scolopendra subspinipes mutilans, by comparing with the toxicity to two other common prey animals, a house mouse, Mus musclus, and a pond frog, Pelophylax nigromaculatus. The lethal doses for mice weighing around 21.5 g and frogs weighing around 3.78 g were less than 5 μl (equivalent to ca. 0.23 and 1.32 μl/g, respectively), which presumably corresponds to an approximate dose of mamushi's one envenomation. On the other hand, centipedes weighing around 1.86 g needed 10-36 μl of venom to die (16.0 μl/g on average). This result suggests that the centipedes are much more resistant to the venom than other prey animals, and it is difficult for mamushi to kill or incapacitate centipedes by the venom of a single envenomation.

Drift, selection and adaptive variation in small populations of a threatened rattlesnake

An important goal of conservation genetics is to determine if the viability of small populations is reduced by a loss of adaptive variation due to genetic drift. Here, we assessed the impact of drift and selection on direct measures of adaptive variation (toxin loci encoding venom proteins) in the eastern massasauga rattlesnake (Sistrurus catenatus), a threatened reptile that exists in small isolated populations. We estimated levels of individual polymorphism in 46 toxin loci and 1,467 control loci across 12 populations of this species, and compared the results with patterns of selection on the same loci following speciation of S. catenatus and its closest relative, the western massasauga (S. tergeminus). Multiple lines of evidence suggest that both drift and selection have had observable impacts on standing adaptive variation. In support of drift effects, we found little evidence for selection on toxin variation within populations and a significant positive relationship between current levels of adaptive variation and long- and short-term estimates of effective population size. However, we also observed levels of directional selection on toxin loci among populations that are broadly similar to patterns predicted from interspecific selection analyses that pre-date the effects of recent drift, and that functional variation in these loci persists despite small short-term effective sizes. This suggests that much of the adaptive variation present in populations may represent an example of "drift debt," a nonequilibrium state where present-day levels of variation overestimate the amount of functional genetic diversity present in future populations.

Interspecific and intraspecific venom enzymatic variation among cobras (Naja sp. and Ophiophagus hannah)

The genera Ophiophagus and Naja comprise part of a clade of snakes referred to as cobras, dangerously venomous front-fanged snakes in the family Elapidae responsible for significant human mortality and morbidity throughout Asia and Africa. We evaluated venom enzyme variation for eleven cobra species and three N. kaouthia populations using SDS-PAGE venom fingerprinting and numerous enzyme assays. Acetylcholinesterase and PLA₂ activities were the most variable between species, and PLA₂ activity was significantly different between Malaysian and Thailand N. kaouthia populations. Venom metalloproteinase activity was low and significantly different among most species, but levels were identical for N. kaouthia populations; minor variation in venom L-amino acid oxidase and phosphodiesterase activities were seen between cobra species. Naja siamensis venom lacked the α-fibrinogenolytic activity common to other cobra venoms. In addition, venom from N. siamensis had no detectable metalloproteinase activity and exhibited an SDS-PAGE profile with reduced abundance of higher mass proteins. Venom profiles from spitting cobras (N. siamensis, N. pallida, and N. mossambica) exhibited similar reductions in higher mass proteins, suggesting the evolution of venoms of reduced complexity and decreased enzymatic activity among spitting cobras. Generally, the venom proteomes of cobras show highly abundant three-finger toxin diversity, followed by large quantities of PLA₂s. However, PLA₂ bands and activity were very reduced for N. haje, N. annulifera and N. nivea. Venom compositionalenzy analysis provides insight into the evolution, diversification and distribution of different venom phenotypes that complements venomic data, and this information is critical for the development of effective antivenoms and snakebite treatment.

Toxin expression in snake venom evolves rapidly with constant shifts in evolutionary rates

Key innovations provide ecological opportunity by enabling access to new resources, colonization of new environments, and are associated with adaptive radiation. The most well-known pattern associated with adaptive radiation is an early burst of phenotypic diversification. Venoms facilitate prey capture and are widely believed to be key innovations leading to adaptive radiation. However, few studies have estimated their evolutionary rate dynamics. Here, we test for patterns of adaptive evolution in venom gene expression data from 52 venomous snake species. By identifying shifts in tempo and mode of evolution along with models of phenotypic evolution, we show that snake venom exhibits the macroevolutionary dynamics expected of key innovations. Namely, all toxin families undergo shifts in their rates of evolution, likely in response to changes in adaptive optima. Furthermore, we show that rapid-pulsed evolution modelled as a Lévy process better fits snake venom evolution than conventional early burst or Ornstein-Uhlenbeck models. While our results support the idea of snake venom being a key innovation, the innovation of venom chemistry lacks clear mechanisms that would lead to reproductive isolation and thus adaptive radiation. Therefore, the extent to which venom directly influences the diversification process is still a matter of contention.

Biological Activities and Proteomic Profile of the Venom of Vipera ursinii ssp., a very Rare Karst Viper from Croatia

The karst viper (Vipera ursinii ssp.) favours high-mountain dry grasslands in southern and south-eastern Croatia. It is medically less important than other Vipera species, because of its remote habitat and the very small amount of venom that it injects by its relatively short fangs. The scientific literature on Vipera ursinii deals mostly with the morphology, ecology and distribution range of this snake, due to the species’ conservation issues, while the toxinological aspects of its venom have not so far been investigated. Here we report on the composition and biological activity of the Vipera ursinii ssp. venom. Using a proteomics approach, we have identified 25 proteins in the venom that belong to seven protein families: snake venom metalloproteinase, serine protease, secreted phospholipase A₂, cysteine-rich secretory protein, snake C-type lectin-like protein, serine protease inhibitor and nerve growth factor. The Vipera ursinii ssp. venom was found to be distinctively insecticidal. Its lethal toxicity towards crickets was more than five times greater than that of Vipera ammodytes ammodytes venom, while the opposite held in mice. Interestingly, the mode of dying after injecting a mouse with Vipera ursinii ssp. venom may suggest the presence of a neurotoxic component. Neurotoxic effects of European vipers have so far been ascribed exclusively to ammodytoxins and ammodytoxin-like basic secreted phospholipases A₂. Structural and immunological analyses of the Vipera ursinii ssp. venom, however, confirmed that ammodytoxin-like proteins are not present in this venom.

Trait differentiation and modular toxin expression in palm-pitvipers

Background

Modularity is the tendency for systems to organize into semi-independent units and can be a key to the evolution and diversification of complex biological systems. Snake venoms are highly variable modular systems that exhibit extreme diversification even across very short time scales. One well-studied venom phenotype dichotomy is a trade-off between neurotoxicity versus hemotoxicity that occurs through the high expression of a heterodimeric neurotoxic phospholipase A₂ (PLA₂) or snake venom metalloproteinases (SVMPs). We tested whether the variation in these venom phenotypes could occur via variation in regulatory sub-modules through comparative venom gland transcriptomics of representative Black-Speckled Palm-Pitvipers (Bothriechis nigroviridis) and Talamancan Palm-Pitvipers (B. nubestris).

Results

We assembled 1517 coding sequences, including 43 toxins for B. nigroviridis and 1787 coding sequences including 42 toxins for B. nubestris. The venom gland transcriptomes were extremely divergent between these two species with one B. nigroviridis exhibiting a primarily neurotoxic pattern of expression, both B. nubestris expressing primarily hemorrhagic toxins, and a second B. nigroviridis exhibiting a mixed expression phenotype. Weighted gene coexpression analyses identified six submodules of transcript expression variation, one of which was highly associated with SVMPs and a second which contained both subunits of the neurotoxic PLA₂ complex. The sub-module association of these toxins suggest common regulatory pathways underlie the variation in their expression and is consistent with known patterns of inheritance of similar haplotypes in other species. We also find evidence that module associated toxin families show fewer gene duplications and transcript losses between species, but module association did not appear to affect sequence diversification.

Conclusion

Sub-modular regulation of expression likely contributes to the diversification of venom phenotypes within and among species and underscores the role of modularity in facilitating rapid evolution of complex traits.

Venomics of the ectoparasitoid wasp Bracon nigricans

BACKGROUND

Venom is one of the most important sources of regulation factors used by parasitic Hymenoptera to redirect host physiology in favour of the developing offspring. This has stimulated a number of studies, both at functional and "omics" level, which, however, are still quite limited for ectophagous parasitoids that permanently paralyze and suppress their victims (i.e., idiobiont parasitoids).

RESULTS

Here we present a combined transcriptomic and proteomic study of the venom of the generalist idiobiont wasp Bracon nigricans, an ectophagous larval parasitoid of different lepidopteran species, for which we recently described the host regulation strategy and the functional role of the venom in the induction of physiological changes in parasitized hosts. The experimental approach used led to the identification of the main components of B. nigricans venom involved in host regulation. Enzymes degrading lipids, proteins and carbohydrates are likely involved in the mobilization of storage nutrients from the fat body and may concurrently be responsible for the release of neurotoxic fatty acids inducing paralysis, and for the modulation of host immune responses.

CONCLUSION

The present work contributes to fill the gap of knowledge on venom composition in ectoparasitoid wasps, and, along with our previous physiological study on this species, provides the foundation on which to develop a functional model of host regulation, based both on physiological and molecular data. This paves the way towards a better understanding of parasitism evolution in the basal lineages of Hymenoptera and to the possible exploitation of venom as source of bioinsecticidal molecules.

The Indian cobra reference genome and transcriptome enables comprehensive identification of venom toxins

Snakebite envenoming is a serious and neglected tropical disease that kills ~100,000 people annually. High-quality, genome-enabled comprehensive characterization of toxin genes will facilitate development of effective humanized recombinant antivenom. We report a de novo near-chromosomal genome assembly of Naja naja, the Indian cobra, a highly venomous, medically important snake. Our assembly has a scaffold N50 of 223.35 Mb, with 19 scaffolds containing 95% of the genome. Of the 23,248 predicted protein-coding genes, 12,346 venom-gland-expressed genes constitute the ‘venom-ome’ and this included 139 genes from 33 toxin families. Among the 139 toxin genes were 19 ‘venom-ome-specific toxins’ (VSTs) that showed venom-gland-specific expression, and these probably encode the minimal core venom effector proteins. Synthetic venom reconstituted through recombinant VST expression will aid in the rapid development of safe and effective synthetic antivenom. Additionally, our genome could serve as a reference for snake genomes, support evolutionary studies and enable venom-driven drug discovery.

Analysis of a near-chromosomal genome assembly and transcriptome profiling of the Indian cobra identifies genes expressed in the venom glands. These data should help develop a new antivenom.

A serine protease-associated lectin in the cytolytic system of blowfly (Chrysomya megacephala) larvae: Evidence and characterization

Cytolytic activity against invading microorganisms is one of the innate forms of immunity in invertebrates. A serine protease-associated sialic acid-specific cytolytic lectin was purified using glutaraldehyde-fixed ox erythrocytes from the larval extract of blowfly (Chrysomya megacephala). The purified lectin lysed vertebrate erythrocytes with effective haemolysis of ox red blood cells (RBCs) in an isotonic medium. The degree of haemolytic (HL) activity of the purified cytolytic lectin depended on its concentration, pH, temperature, and calcium ions. It was sensitive to ethylenediaminetetraacetic acid. The native molecular mass of the C-type lectin was 260 ± 26 kDa, comprising four different polypeptide subunits of 75 kDa (pI ~8), 69 kDa (pI ~7.0), 61 kDa (pI ~5.3), and 55 kDa (pI ~4.6). The association between the C-type lectin and serine protease was confirmed by MALDI-TOF-MS analysis that revealed its homology in the same spectral peak as well as the proteases and phenylmethylsulphonyl fluoride inhibition of HL activity. Haemolysis inhibition by N-acetylneuraminic acid and other sugars revealed the properties of the lectin. The purified lectin distorted the integrity of ox RBCs and Paenalcaligenes hermetiae. This in vitro study documents the presence of a cytolytic system in blowfly (C. megacephala) larvae for the clearance of invading microbial pathogens in their feeding niche.

Omics Technologies for Profiling Toxin Diversity and Evolution in Snake Venom: Impacts on the Discovery of Therapeutic and Diagnostic Agents

Snake venoms are primarily composed of proteins and peptides, and these toxins have developed high selectivity to their biological targets. This makes venoms interesting for exploration into protein evolution and structure-function relationships. A single venom protein superfamily can exhibit a variety of pharmacological effects; these variations in activity originate from differences in functional sites, domains, posttranslational modifications, and the formations of toxin complexes. In this review, we discuss examples of how the major venom protein superfamilies have diversified, as well as how newer technologies in the omics fields, such as genomics, transcriptomics, and proteomics, can be used to characterize both known and unknown toxins.Because toxins are bioactive molecules with a rich diversity of activities, they can be useful as therapeutic and diagnostic agents, and successful examples of toxin applications in these areas are also reviewed. With the current rapid pace of technology, snake venom research and its applications will only continue to expand.

Alternative mRNA Splicing in Three Venom Families Underlying a Possible Production of Divergent Venom Proteins of the Habu Snake, Protobothrops flavoviridis

Snake venoms are complex mixtures of toxic proteins encoded by various gene families that function synergistically to incapacitate prey. A huge repertoire of snake venom genes and proteins have been reported, and alternative splicing is suggested to be involved in the production of divergent gene transcripts. However, a genome-wide survey of the transcript repertoire and the extent of alternative splicing still remains to be determined. In this study, the comprehensive analysis of transcriptomes in the venom gland was achieved by using PacBio sequencing. Extensive alternative splicing was observed in three venom protein gene families, metalloproteinase (MP), serine protease (SP), and vascular endothelial growth factors (VEGF). Eleven MP and SP genes and a VEGF gene are expressed as a total of 81, 61, and 8 transcript variants, respectively. In the MP gene family, individual genes are transcribed into different classes of MPs by alternative splicing. We also observed trans-splicing among the clustered SP genes. No other venom genes as well as non-venom counterpart genes exhibited alternative splicing. Our results thus indicate a potential contribution of mRNA alternative and trans-splicing in the production of highly variable transcripts of venom genes in the habu snake.

Venom Composition in a Phenotypically Variable Pit Viper ( Trimeresurus insularis) across the Lesser Sunda Archipelago

The genus Trimeresurus comprises a group of venomous pitvipers endemic to Southeast Asia and the Pacific Islands. Of these, Trimeresurus insularis, the White-lipped Island Pitviper, is a nocturnal, arboreal species that occurs on nearly every major island of the Lesser Sunda archipelago. In the current study, venom phenotypic characteristics of T. insularis sampled from eight Lesser Sunda Islands (Flores, Lembata, Lombok, Pantar, Sumba, Sumbawa, Timor, and Wetar) were evaluated via SDS-PAGE, enzymatic activity assays, fibrinogenolytic assays, gelatin zymography, and RP-HPLC, and the Sumbawa sample was characterized by venomic analysis. For additional comparative analyses, venoms were also examined from several species in the Trimeresurus complex, including T. borneensis, T. gramineus, T. puniceus, T. purpureomaculatus, T. stejnegeri, and Protobothrops flavoviridis. Despite the geographical isolation, T. insularis venoms from all eight islands demonstrated remarkable similarities in gel electrophoretic profiles and RP-HPLC patterns, and all populations had protein bands in the mass ranges of phosphodiesterases (PDE), l-amino acid oxidases (LAAO), P-III snake venom metalloproteinases (SVMP), serine proteases, cysteine-rich secretory proteins (CRISP), phospholipases A₂ (PLA₂), and C-type lectins. An exception was observed in the Lombok sample, which lacked protein bands in the mass range of serine protease and CRISP. Venomic analysis of the Sumbawa venom also identified these protein families, in addition to several proteins of lesser abundance (<1%), including glutaminyl cyclase, aminopeptidase, PLA₂ inhibitor, phospholipase B, cobra venom factor, 5'-nucleotidase, vascular endothelial growth factor, and hyaluronidase. All T. insularis venoms exhibited similarities in thrombin-like and PDE activities, while significant differences were observed for LAAO, SVMP, and kallikrein-like activities, though these differences were only observed for a few islands. Slight but noticeable differences were also observed with fibrinogen and gelatin digestion activities. Trimeresurus insularis venoms exhibited overall similarity to the other Trimeresurus complex species examined, with the exception of P. flavoviridis venom, which showed the greatest overall differentiation. Western blot analysis revealed that all major T. insularis venom proteins were recognized by Green Pitviper ( T. albolabris) antivenom, and reactivity was also seen with most venom proteins of the other Trimeresurus species, but incomplete antivenom-venom recognition was observed against P. flavoviridis venom proteins. These results demonstrate significant conservation in the venom composition of T. insularis across the Lesser Sunda archipelago relative to the other Trimeresurus species examined.

Snake venom NAD glycohydrolases: primary structures, genomic location, and gene structure

NAD glycohydrolase (EC 3.2.2.5) (NADase) sequences have been identified in 10 elapid and crotalid venom gland transcriptomes, eight of which are complete. These sequences show very high homology, but elapid and crotalid sequences also display consistent differences. As in Aplysia kurodai ADP-ribosyl cyclase and vertebrate CD38 genes, snake venom NADase genes comprise eight exons; however, in the Protobothrops mucrosquamatus genome, the sixth exon is sometimes not transcribed, yielding a shortened NADase mRNA that encodes all six disulfide bonds, but an active site that lacks the catalytic glutamate residue. The function of this shortened protein, if expressed, is unknown. While many vertebrate CD38s are multifunctional, liberating both ADP-ribose and small quantities of cyclic ADP-ribose (cADPR), snake venom CD38 homologs are dedicated NADases. They possess the invariant TLEDTL sequence (residues 144–149) that bounds the active site and the catalytic residue, Glu228. In addition, they possess a disulfide bond (Cys121–Cys202) that specifically prevents ADP-ribosyl cyclase activity in combination with Ile224, in lieu of phenylalanine, which is requisite for ADPR cyclases. In concert with venom phosphodiesterase and 5′-nucleotidase and their ecto-enzyme homologs in prey tissues, snake venom NADases comprise part of an envenomation strategy to liberate purine nucleosides, and particularly adenosine, in the prey, promoting prey immobilization via hypotension and paralysis.

Organic and Peptidyl Constituents of Snake Venoms: The Picture Is Vastly More Complex Than We Imagined

Small metabolites and peptides in 17 snake venoms (Elapidae, Viperinae, and Crotalinae), were quantified using liquid chromatography-mass spectrometry. Each venom contains >900 metabolites and peptides. Many small organic compounds are present at levels that are probably significant in prey envenomation, given that their known pharmacologies are consistent with snake envenomation strategies. Metabolites included purine nucleosides and their bases, neurotransmitters, neuromodulators, guanidino compounds, carboxylic acids, amines, mono- and disaccharides, and amino acids. Peptides of 2⁻15 amino acids are also present in significant quantities, particularly in crotaline and viperine venoms. Some constituents are specific to individual taxa, while others are broadly distributed. Some of the latter appear to support high anabolic activity in the gland, rather than having toxic functions. Overall, the most abundant organic metabolite was citric acid, owing to its predominance in viperine and crotaline venoms, where it chelates divalent cations to prevent venom degradation by venom metalloproteases and damage to glandular tissue by phospholipases. However, in terms of their concentrations in individual venoms, adenosine, adenine, were most abundant, owing to their high titers in Dendroaspis polylepis venom, although hypoxanthine, guanosine, inosine, and guanine all numbered among the 50 most abundant organic constituents. A purine not previously reported in venoms, ethyl adenosine carboxylate, was discovered in D. polylepis venom, where it probably contributes to the profound hypotension caused by this venom. Acetylcholine was present in significant quantities only in this highly excitotoxic venom, while 4-guanidinobutyric acid and 5-guanidino-2-oxopentanoic acid were present in all venoms.

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

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

Comprehensive Snake Venomics of the Okinawa Habu Pit Viper, Protobothrops flavoviridis, by Complementary Mass Spectrometry-Guided Approaches

The Asian world is home to a multitude of venomous and dangerous snakes, which are used to induce various medical effects in the preparation of traditional snake tinctures and alcoholics, like the Japanese snake wine, named Habushu. The aim of this work was to perform the first quantitative proteomic analysis of the Protobothrops flavoviridis pit viper venom. Accordingly, the venom was analyzed by complimentary bottom-up and top-down mass spectrometry techniques. The mass spectrometry-based snake venomics approach revealed that more than half of the venom is composed of different phospholipases A2 (PLA₂). The combination of this approach and an intact mass profiling led to the identification of the three main Habu PLA₂s. Furthermore, nearly one-third of the total venom consists of snake venom metalloproteinases and disintegrins, and several minor represented toxin families were detected: C-type lectin-like proteins (CTL), cysteine-rich secretory proteins (CRISP), snake venom serine proteases (svSP), l-amino acid oxidases (LAAO), phosphodiesterase (PDE) and 5'-nucleotidase. Finally, the venom of P. flavoviridis contains certain bradykinin-potentiating peptides and related peptides, like the svMP inhibitors, pEKW, pEQW, pEEW and pENW. In preliminary MTT cytotoxicity assays, the highest cancerous-cytotoxicity of crude venom was measured against human neuroblastoma SH-SY5Y cells and shows disintegrin-like effects in some fractions.

The habu genome reveals accelerated evolution of venom protein genes

Evolution of novel traits is a challenging subject in biological research. Several snake lineages developed elaborate venom systems to deliver complex protein mixtures for prey capture. To understand mechanisms involved in snake venom evolution, we decoded here the ~1.4-Gb genome of a habu, Protobothrops flavoviridis. We identified 60 snake venom protein genes (SV) and 224 non-venom paralogs (NV), belonging to 18 gene families. Molecular phylogeny reveals early divergence of SV and NV genes, suggesting that one of the four copies generated through two rounds of whole-genome duplication was modified for use as a toxin. Among them, both SV and NV genes in four major components were extensively duplicated after their diversification, but accelerated evolution is evident exclusively in the SV genes. Both venom-related SV and NV genes are significantly enriched in microchromosomes. The present study thus provides a genetic background for evolution of snake venom composition.

In vitro effects of Crotalus atrox snake venom on chick and mouse neuromuscular preparations

The neuromuscular effect of venoms is not a major clinical manifestation shared between rattlesnakes native to the Americas, which showed two different venom phenotypes. Taking into account this dichotomy, nerve muscle preparations from mice and chicks were used to investigate the ability of Crotalus atrox venom to induce in vitro neurotoxicity and myotoxicity. Unlike crotalic venoms of South America, low concentrations of C. atrox venom did not result in significant effects on mouse neuromuscular preparations. The venom was more active on avian nerve-muscle, showing reduction of twitch heights after 120 min of incubation with 10, 30 and 100 μg/mL of venom with diminished responses to agonists and KCl. Histological analysis highlighted that C. atrox was myotoxic in both species of experimental animals; as evidenced by degenerative events, including edematous cells, delta lesions, hypercontracted fibers and muscle necrosis, which can lead to neurotoxic action. These results provide key insights into the myotoxicity and low neurotoxicity of C. atrox in two animal models, corroborating with previous genomic and proteomic findings and would be useful for a deeper understanding of venom evolution in snakes belonging to the genus Crotalus.

Population Genomic Analysis of a Pitviper Reveals Microevolutionary Forces Underlying Venom Chemistry

Venoms are among the most biologically active secretions known, and are commonly believed to evolve under extreme positive selection. Many venom gene families, however, have undergone duplication, and are often deployed in doses vastly exceeding the LD₅₀ for most prey species, which should reduce the strength of positive selection. Here, we contrast these selective regimes using snake venoms, which consist of rapidly evolving protein formulations. Though decades of extensive studies have found that snake venom proteins are subject to strong positive selection, the greater action of drift has been hypothesized, but never tested. Using a combination of de novo genome sequencing, population genomics, transcriptomics, and proteomics, we compare the two modes of evolution in the pitviper, Protobothrops mucrosquamatus. By partitioning selective constraints and adaptive evolution in a McDonald–Kreitman-type framework, we find support for both hypotheses: venom proteins indeed experience both stronger positive selection, and lower selective constraint than other genes in the genome. Furthermore, the strength of selection may be modulated by expression level, with more abundant proteins experiencing weaker selective constraint, leading to the accumulation of more deleterious mutations. These findings show that snake venoms evolve by a combination of adaptive and neutral mechanisms, both of which explain their extraordinarily high rates of molecular evolution. In addition to positive selection, which optimizes efficacy of the venom in the short term, relaxed selective constraints for deleterious mutations can lead to more rapid turnover of individual proteins, and potentially to exploration of a larger venom phenotypic space.

Structures of N-Glycans of Bothrops Venoms Revealed as Molecular Signatures that Contribute to Venom Phenotype in Viperid Snakes

The complexity of snake venoms has long been investigated to explore a myriad of biologically active proteins and peptides that are used for immobilizing or killing prey, and are responsible for the pathological effects observed on envenomation. Glycosylation is the main post-translational modification (PTM) of viperid venoms but currently there is little understanding of how protein glycosylation impacts the variation of venom proteomes. We have previously reported that Bothrops venom glycoproteomes contain a core of components that markedly define their composition and parallel their phylogenetic classification. Here we extend those observations to eight Bothrops species evaluating the N-glycomes by LC-MS as assigned cartoon structures and detailing those structures separately as methylated analogs using ion-trap mass spectrometry (MSⁿ). Following ion disassembly through multiple steps provided sequence and linkage isomeric details that characterized 52 unique compositions in Bothrops venoms. These occurred as 60 structures, of which 26 were identified in the venoms of the Jararaca Complex (B. alcatraz, B. insularis, and B. jararaca), 20 in B. erythromelas, B. jararacussu, B. moojeni and B. neuwiedi venoms, and 22 in B. cotiara venom. Further, quantitative analysis of these N-glycans showed variable relative abundances in the venoms. For the first time a comprehensive set of N-glycan structures present in snake venoms are defined. Despite the fact that glycosylation is not template-defined, the N-glycomes of these venoms mirror the phylogeny cladograms of South American bothropoid snakes reported in studies on morphological, molecular data and feeding habits, exhibiting distinct molecular signatures for each venom. Considering the complexity of N-glycan moieties generally found in glycoproteins, characterized by different degrees of branching, isomer structures, and variable abundances, our findings point to these factors as another level of complexity in Bothrops venoms, features that could dramatically contribute to their distinct biological activities.

Synergistic strategies of predominant toxins in snake venoms

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

Transcriptomics-guided bottom-up and top-down venomics of neonate and adult specimens of the arboreal rear-fanged Brown Treesnake, Boiga irregularis, from Guam

The Brown Treesnake (Boiga irregularis) is an arboreal, nocturnal, rear-fanged venomous snake native to northern and eastern regions of Australia, Papua New Guinea and the Solomon Islands. It was inadvertently introduced onto the island of Guam during the late 1940's to early 1950's, and it has caused massive declines and extirpations of the native bird, lizard, and mammal populations. In the current study, we report the characterization of the venom proteome of an adult and a neonate B. irregularis specimens from Guam by a combination of venom gland transcriptomic and venomic analyses. Venom gland transcriptomic analysis of an adult individual identified toxins belonging to 18 protein families, with three-finger toxin isoforms being the most abundantly expressed transcripts, comprising 94% of all venom protein transcript reads. Transcripts for PIII-metalloproteinases, C-type lectins, cysteine-rich secretory proteins, acetylcholinesterases, natriuretic peptides, ficolins, phospholipase A₂ (PLA₂) inhibitors, PLA₂s, vascular endothelial growth factors, Kunitz-type protease inhibitors, cystatins, phospholipase Bs, cobra venom factors, waprins, SVMP inhibitors, matrix metalloproteinases, and hyaluronidases were also identified, albeit, at very low abundances ranging from 0.05% to 1.7% of the transcriptome. The venom proteomes of neonate and adult B. irregularis were also both overwhelmingly (78 and 84%, respectively) dominated by monomeric and dimeric 3FTxs, followed by moderately abundant (21% (N) and 13% (A)) CRISPs, low abundance (1% (N), 3% (A)) PIII-SVMPs, and very low abundance (<0.01%) PLA₂ and SVMP inhibitors. The differences in relative toxin abundances identified between neonate and adult snakes likely correlates to shifts in prey preference between the two age classes, from nearly-exclusively lizards to lizards, birds and small mammals. Immunoaffinity antivenomics with experimentally designed rabbit anti-Brown Treesnake (anti-BTS) venom IgGs against homologous venom from adult snakes demonstrated that CRISPs, PIII-SVMPs, and 60-70% of 3FTxs were effectively immunocaptured. Western blot analysis showed that all venom proteins were recognized by anti-BTS IgGs, and cross-reactivity with other rear-fanged snake venoms was also observed. Incubation of anti-BTS venom IgGs with crude B. irregularis venom resulted in a significant decrease in proteolytic (SVMP) activity against azocasein. These results provide the first comparative venomic and anti-venomic analysis of neonate and adult B. irregularis from Guam, further highlighting evolutionary trends in venom composition among rear-fanged venomous snakes.

SIGNIFICANCE PARAGRAPH

The Brown Treesnake (Boiga irregularis) has caused extensive ecological and economic damage to the island of Guam where it has become a classic example of the negative impacts of invasive species. In the current study, we report the first combined transcriptomic and proteomic analysis of B. irregularis venom of Guam origin. The transcriptome of an adult snake contained toxin sequences belonging to 18 protein families, with three-finger toxin (3FTx) isoforms being the most abundant and representing 94% of all venom protein transcript reads. Our bottom-up and top-down venomic analyses confirmed that 3FTxs are the major components of B. irregularis venom, and a comparative analysis of neonate and adult venoms demonstrate a clear ontogenetic shift in toxin abundance, likely driven by dietary variation between the two age classes. Second-generation antivenomics and Western blot analysis using purified anti-Brown Treesnake rabbit serum IgGs (anti-BTS IgGs) showed strong immunoreactivity toward B. irregularis venom. Interestingly, our anti-BTS IgGs did not cross-react with 3FTxs found in several other rear-fanged snake venoms, or against 3FTxs in the venom of the elapid Ophiophagus hannah, indicating that epitopes in these 3FTx molecules are quite distinct.

Computational Studies of Snake Venom Toxins

Most snake venom toxins are proteins, and participate to envenomation through a diverse array of bioactivities, such as bleeding, inflammation, and pain, cytotoxic, cardiotoxic or neurotoxic effects. The venom of a single snake species contains hundreds of toxins, and the venoms of the 725 species of venomous snakes represent a large pool of potentially bioactive proteins. Despite considerable discovery efforts, most of the snake venom toxins are still uncharacterized. Modern bioinformatics tools have been recently developed to mine snake venoms, helping focus experimental research on the most potentially interesting toxins. Some computational techniques predict toxin molecular targets, and the binding mode to these targets. This review gives an overview of current knowledge on the ~2200 sequences, and more than 400 three-dimensional structures of snake toxins deposited in public repositories, as well as of molecular modeling studies of the interaction between these toxins and their molecular targets. We also describe how modern bioinformatics have been used to study the snake venom protein phospholipase A2, the small basic myotoxin Crotamine, and the three-finger peptide Mambalgin.

Tracing the evolution of venom phospholipases A2 in Gloydius strauchii and related pitvipers: A tale of two acidic isozymes

Two acidic Asp⁴⁹-PLA₂s with Glu⁶ substitution and a neutral Lys⁴⁹-PLA (designated Gst-K49) were cloned from G. strauchii venom glands, their full amino acid sequences were deduced. The predominant acidic PLA₂ (designated Gst-E6a) contains 124 residues and the M¹⁸W³⁰ substitutions, while the minor acidic PLA₂ (designated Gst-E6b) contains 122 residues and the V¹⁸A³⁰ substitutions. Their sequences are most similar to those of the respective orthologous PLA₂s of G. intermedius venom. Gst-E6a and Gst-E6b appear to be paralogs and possibly have different predatory targets or functions. The LC-MS/MS results indicate the presence of only three PLA₂ gene products in the crude venom, the relative expression levels were in the order of Gst-E6a ≫ Gst-E6b > Gst-K49, as confirmed by qPCR results. In contrast to other Gloydius, G. strauchii venom does not contain neurotoxic or basic anticoagulant Asp⁴⁹-PLA₂s, but Gst-K49 is the first Lys⁴⁹-PLA₂ identified in Gloydius venoms. However, its venom content is relatively low and its pI value 7.3 is much lower than those of other Lys⁴⁹-PLA₂s and. The Lys⁴⁹-PLA₂ genes appear to regress in the venom of most of Gloydius and related rattlesnake, and this evolutionary regression occurred before the dispersal of Asian pitvipers to the New World.

A Review and Database of Snake Venom Proteomes

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

Coralsnake Venomics: Analyses of Venom Gland Transcriptomes and Proteomes of Six Brazilian Taxa

Venom gland transcriptomes and proteomes of six Micrurus taxa (M. corallinus, M. lemniscatus carvalhoi, M. lemniscatus lemniscatus, M. paraensis, M. spixii spixii, and M. surinamensis) were investigated, providing the most comprehensive, quantitative data on Micrurus venom composition to date, and more than tripling the number of Micrurus venom protein sequences previously available. The six venomes differ dramatically. All are dominated by 2-6 toxin classes that account for 91-99% of the toxin transcripts. The M. s. spixii venome is compositionally the simplest. In it, three-finger toxins (3FTxs) and phospholipases A₂ (PLA₂s) comprise >99% of the toxin transcripts, which include only four additional toxin families at levels ≥0.1%. Micrurus l. lemniscatus venom is the most complex, with at least 17 toxin families. However, in each venome, multiple structural subclasses of 3FTXs and PLA₂s are present. These almost certainly differ in pharmacology as well. All venoms also contain phospholipase B and vascular endothelial growth factors. Minor components (0.1-2.0%) are found in all venoms except that of M. s. spixii. Other toxin families are present in all six venoms at trace levels (<0.005%). Minor and trace venom components differ in each venom. Numerous novel toxin chemistries include 3FTxs with previously unknown 8- and 10-cysteine arrangements, resulting in new 3D structures and target specificities. 9-cysteine toxins raise the possibility of covalent, homodimeric 3FTxs or heterodimeric toxins with unknown pharmacologies. Probable muscarinic sequences may be reptile-specific homologs that promote hypotension via vascular mAChRs. The first complete sequences are presented for 3FTxs putatively responsible for liberating glutamate from rat brain synaptosomes. Micrurus C-type lectin-like proteins may have 6-9 cysteine residues and may be monomers, or homo- or heterodimers of unknown pharmacology. Novel KSPIs, 3× longer than any seen previously, appear to have arisen in three species by gene duplication and fusion. Four species have transcripts homologous to the nociceptive toxin, (MitTx) α-subunit, but all six species had homologs to the β-subunit. The first non-neurotoxic, non-catalytic elapid phospholipase A₂s are reported. All are probably myonecrotic. Phylogenetic analysis indicates that the six taxa diverged 15-35 million years ago and that they split from their last common ancestor with Old World elapines nearly 55 million years ago. Given their early diversification, many cryptic micrurine taxa are anticipated.