Integrated "omics" profiling indicates that miRNAs are modulators of the ontogenetic venom composition shift in the Central American rattlesnake, Crotalus simus simus

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.

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.

Variability in antivenom neutralization of Mexican viperid snake venoms

BACKGROUND

Each year, 3,800 cases of snakebite envenomation are reported in Mexico, resulting in 35 fatalities. The only scientifically validated treatment for snakebites in Mexico is the use of antivenoms. Currently, two antivenoms are available in the market, with one in the developmental phase. These antivenoms, produced in horses, consist of F(ab')2 fragments generated using venoms from various species as immunogens. While previous studies primarily focused on neutralizing the venom of the Crotalus species, our study aims to assess the neutralization capacity of different antivenom batches against pit vipers from various genera in Mexico.

METHODOLOGY

We conducted various biological and biochemical tests to characterize the venoms. Additionally, we performed neutralization tests using all three antivenoms to evaluate their effectiveness against lethal activity and their ability to neutralize proteolytic and fibrinogenolytic activities.

RESULTS

Our results reveal significant differences in protein content and neutralizing capacity among different antivenoms and even between different batches of the same product. Notably, the venom of Crotalus atrox is poorly neutralized by all evaluated batches despite being the primary cause of envenomation in the country's northern region. Furthermore, even at the highest tested concentrations, no antivenom could neutralize the lethality of Metlapilcoatlus nummifer and Porthidium yucatanicum venoms. These findings highlight crucial areas for improving existing antivenoms and developing new products.

CONCLUSION

Our research reveals variations in protein content and neutralizing potency among antivenoms, emphasizing the need for consistency in venom characteristics as immunogens. While Birmex neutralizes more LD50 per vial, Antivipmyn excels in specific neutralization. The inability of antivenoms to neutralize certain venoms, especially M. nummifer and P. yucatanicum, highlights crucial improvement opportunities, given the medical significance of these species.

The genetic regulatory architecture and epigenomic basis for age-related changes in rattlesnake venom

Developmental phenotypic changes can evolve under selection imposed by age- and size-related ecological differences. Many of these changes occur through programmed alterations to gene expression patterns, but the molecular mechanisms and gene-regulatory networks underlying these adaptive changes remain poorly understood. Many venomous snakes, including the eastern diamondback rattlesnake (Crotalus adamanteus), undergo correlated changes in diet and venom expression as snakes grow larger with age, providing models for identifying mechanisms of timed expression changes that underlie adaptive life history traits. By combining a highly contiguous, chromosome-level genome assembly with measures of expression, chromatin accessibility, and histone modifications, we identified cis-regulatory elements and trans-regulatory factors controlling venom ontogeny in the venom glands of C. adamanteus. Ontogenetic expression changes were significantly correlated with epigenomic changes within genes, immediately adjacent to genes (e.g., promoters), and more distant from genes (e.g., enhancers). We identified 37 candidate transcription factors (TFs), with the vast majority being up-regulated in adults. The ontogenetic change is largely driven by an increase in the expression of TFs associated with growth signaling, transcriptional activation, and circadian rhythm/biological timing systems in adults with corresponding epigenomic changes near the differentially expressed venom genes. However, both expression activation and repression contributed to the composition of both adult and juvenile venoms, demonstrating the complexity and potential evolvability of gene regulation for this trait. Overall, given that age-based trait variation is common across the tree of life, we provide a framework for understanding gene-regulatory-network-driven life-history evolution more broadly.

Distinct regulatory networks control toxin gene expression in elapid and viperid snakes

BACKGROUND

Venom systems are ideal models to study genetic regulatory mechanisms that underpin evolutionary novelty. Snake venom glands are thought to share a common origin, but there are major distinctions between venom toxins from the medically significant snake families Elapidae and Viperidae, and toxin gene regulatory investigations in elapid snakes have been limited. Here, we used high-throughput RNA-sequencing to profile gene expression and microRNAs between active (milked) and resting (unmilked) venom glands in an elapid (Eastern Brown Snake, Pseudonaja textilis), in addition to comparative genomics, to identify cis- and trans-acting regulation of venom production in an elapid in comparison to viperids (Crotalus viridis and C. tigris).

RESULTS

Although there is conservation in high-level mechanistic pathways regulating venom production (unfolded protein response, Notch signaling and cholesterol homeostasis), there are differences in the regulation of histone methylation enzymes, transcription factors, and microRNAs in venom glands from these two snake families. Histone methyltransferases and transcription factor (TF) specificity protein 1 (Sp1) were highly upregulated in the milked elapid venom gland in comparison to the viperids, whereas nuclear factor I (NFI) TFs were upregulated after viperid venom milking. Sp1 and NFI cis-regulatory elements were common to toxin gene promoter regions, but many unique elements were also present between elapid and viperid toxins. The presence of Sp1 binding sites across multiple elapid toxin gene promoter regions that have been experimentally determined to regulate expression, in addition to upregulation of Sp1 after venom milking, suggests this transcription factor is involved in elapid toxin expression. microRNA profiles were distinctive between milked and unmilked venom glands for both snake families, and microRNAs were predicted to target a diversity of toxin transcripts in the elapid P. textilis venom gland, but only snake venom metalloproteinase transcripts in the viperid C. viridis venom gland. These results suggest differences in toxin gene posttranscriptional regulation between the elapid P. textilis and viperid C. viridis.

CONCLUSIONS

Our comparative transcriptomic and genomic analyses between toxin genes and isoforms in elapid and viperid snakes suggests independent toxin regulation between these two snake families, demonstrating multiple different regulatory mechanisms underpin a venomous phenotype.

Venom Composition of Neglected Bothropoid Snakes from the Amazon Rainforest: Ecological and Toxinological Implications

Snake venoms have evolved in several families of Caenophidae, and their toxins have been assumed to be biochemical weapons with a role as a trophic adaptation. However, it remains unclear how venom contributes to the success of venomous species for adaptation to different environments. Here we compared the venoms from Bothrocophias hyoprora, Bothrops taeniatus, Bothrops bilineatus smaragdinus, Bothrops brazili, and Bothrops atrox collected in the Amazon Rainforest, aiming to understand the ecological and toxinological consequences of venom composition. Transcriptomic and proteomic analyses indicated that the venoms presented the same toxin groups characteristic from bothropoids, but with distinct isoforms with variable qualitative and quantitative abundances, contributing to distinct enzymatic and toxic effects. Despite the particularities of each venom, commercial Bothrops antivenom recognized the venom components and neutralized the lethality of all species. No clear features could be observed between venoms from arboreal and terrestrial habitats, nor in the dispersion of the species throughout the Amazon habitats, supporting the notion that venom composition may not shape the ecological or toxinological characteristics of these snake species and that other factors influence their foraging or dispersal in different ecological niches.

Identifying the geographic distribution pattern of venomous snakes and regions of high snakebite risk in Iran

We describe species richness patterns of venomous snakes in Iran in order to produce snakebite risk prediction maps and identify gaps in regional health care centers capable of managing snakebites. We digitized distribution maps from the literature, Global Biodiversity Information Facility (GBIF), and the results of our own field studies of 24 terrestrial venomous snake species (including 4 endemic to Iran). Species richness patterns were associated with eight environmental factors. The variables have been extracted from the WorldClim dataset (bio12 = annual precipitation, bio15 = precipitation seasonality, bio17 = precipitation of the driest quarter, bio2 = mean diurnal range, bio3 = isothermality (bio2/bio7), bio4 = temperature seasonality, bio9 = mean temperature of the driest quarter and slope). Based on spatial analyses, species richness in Iran is highly affected by three environmental variables (bio12, 15, and 17) associated with precipitation. The relationship patterns among these predictors and species richness were strong and linear. The hotspot regions for venomous snakes species are concentrated on the western to southwestern and north to northeastern regions of Iran, which is partially consistent with the known Irano-Anatolian biodiversity hotspot. Because of the high number of endemic species and climatic conditions on the Iranian Plateau, the venoms of snakes distributed in those areas may contain novel properties and components.

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.

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

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

Venom biotechnology: casting light on nature's deadliest weapons using synthetic biology

Venoms are complex chemical arsenals that have evolved independently many times in the animal kingdom. Venoms have attracted the interest of researchers because they are an important innovation that has contributed greatly to the evolutionary success of many animals, and their medical relevance offers significant potential for drug discovery. During the last decade, venom research has been revolutionized by the application of systems biology, giving rise to a novel field known as venomics. More recently, biotechnology has also made an increasing impact in this field. Its methods provide the means to disentangle and study venom systems across all levels of biological organization and, given their tremendous impact on the life sciences, these pivotal tools greatly facilitate the coherent understanding of venom system organization, development, biochemistry, and therapeutic activity. Even so, we lack a comprehensive overview of major advances achieved by applying biotechnology to venom systems. This review therefore considers the methods, insights, and potential future developments of biotechnological applications in the field of venom research. We follow the levels of biological organization and structure, starting with the methods used to study the genomic blueprint and genetic machinery of venoms, followed gene products and their functional phenotypes. We argue that biotechnology can answer some of the most urgent questions in venom research, particularly when multiple approaches are combined together, and with other venomics technologies.

The Middle Eastern Cousin: Comparative Venomics of Daboia palaestinae and Daboia russelii

Among the medically most important snakes in the world, the species belonging to the genus Daboia have been attributed to the highest number of human envenomings, deaths and disabilities. Given their significant clinical relevance, the venoms of Russell's vipers (D. russelii and D. siamensis) have been the primary focus of research. In contrast, the composition, activity, ecology and evolution of venom of its congener, the Palestine viper (D. palaestinae), have remained largely understudied. Therefore, to unravel the factors responsible for the enhanced medical relevance of D. russelii in comparison to D. palaestinae, we comparatively evaluated their venom proteomes, biochemical activities, and mortality and morbidity inflicting potentials. Furthermore, the synthesis and regulation of venom in snakes have also remained underinvestigated, and the relative contribution of each venom gland remains unclear. We address this knowledge gap by sequencing the tissue transcriptomes of both venom glands of D. palaestinae, and comparatively evaluating their contribution to the secreted venom concoction. Our findings highlight the disparity in the venom composition, function and toxicities of the two Daboia species. We also show that toxin production is not partitioned between the two venom glands of D. palaestinae.

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

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

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

Simple Summary

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

Abstract

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

Intraspecific Differences in the Venom of Crotalus durissus cumanensis from Colombia

Biochemical and biological differences in the venom of Crotalus durissus cumanensis from three ecoregions of Colombia were evaluated. Rattlesnakes were collected from the geographic areas of Magdalena Medio (MM), Caribe (CA) and Orinoquía (OR). All three regionally distributed venoms contain proteases, PLA₂s and the basic subunit of crotoxin. However, only crotamine was detected in the CA venom. The highest lethality, coagulant, phospholipase A₂ and hyaluronidase activities were found in the MM venom. Also, some differences, observed by western blot and immunoaffinity, were found in all three venoms when using commercial antivenoms. Furthermore, all three eco-regional venoms showed intraspecific variability, considering the differences in the abundance and intensity of their components, in addition to the activity and response to commercial antivenoms.

Altered RNome expression in Murine Gastrocnemius Muscle following Exposure to Jararhagin, a Metalloproteinase from Bothrops jararaca Venom

Small RNAs (sRNAs) and microRNAs (miRNAs) are small endogenous noncoding single-stranded RNAs that regulate gene expression in eukaryotes. Experiments in mice and humans have revealed that a typical small RNA can affect the expression of a wide range of genes, implying that small RNAs function as global regulators. Here, we used small RNA deep sequencing to investigate how jararhagin, a metalloproteinase toxin produced from the venom of Bothrops jararaca, affected mmu-miRNAs expression in mice 2 hours (Jar 2hrs) and 24 hours (Jar 24hrs) after injection compared to PBS control. The findings revealed that seven mmu-miRNAs were substantially differentially expressed (p value (p (Corr) cut-off 0.05, fold change ≥ 2) at 2 hrs after jararhagin exposure and that the majority of them were upregulated when compared to PBS. In contrast to these findings, a comparison of Jar 24hrs vs. PBS 24hrs demonstrated that the majority of identified mmu-miRNAs were downregulated. Furthermore, the studies demonstrated that mmu-miRNAs can target the expression of several genes involved in the MAPK signaling pathway. The steady antithetical regulation of mmu-miRNAs may correlate with the expression of genes that trigger apoptosis via MAPK in the early stages, and this effect intensifies with time. The findings expand our understanding of the effects of jararhagin on local tissue lesions at the molecular level.

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

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

Biochemical and computational approaches to understand venom toxin-toxin interaction

Snake venoms are mainly composed of proteins and peptides (venom toxins). The venom transcriptomes and proteomes have been extensively investigated; however, venom toxin-toxin interactions remain poorly characterized. We detected the interaction of venom Asp49-PLA₂ and 3FTx using biochemical and computational approaches. A stable structure of Asp49-PLA₂-3FTx was identified, and the interface of Asp49-PLA₂ and 3FTx was analyzed. The approaches will shed light on understanding the venom complexity and deciphering the synergistic effects of venom toxins.

Interpopulational variation and ontogenetic shift in the venom composition of Lataste's viper (Vipera latastei, Boscá 1878) from northern Portugal

Lataste's viper (Vipera latastei) is a venomous European viper endemic to the Iberian Peninsula, recognised as medically important by the World Health Organization. To date, no comprehensive characterisation of this species' venom has been reported. Here, we analysed the venoms of juvenile and adult specimens of V. latastei from two environmentally different populations from northern Portugal. Using bottom-up venomics, we produced six venom proteomes (three per population) from vipers belonging to both age classes (i.e., two juveniles and four adults), and RP-HPLC profiles of 54 venoms collected from wild specimens. Venoms from juveniles and adults differed in their chromatographic profiles and relative abundances of their toxins, suggesting the occurrence of ontogenetic changes in venom composition. Specifically, snake venom metalloproteinase (SVMP) was the most abundant toxin family in juvenile venoms, while snake venom serine proteinases (SVSPs), phospholipases A₂ (PLA₂s), and C-type lectin-like (CTLs) proteins were the main toxins comprising adult venoms. The RP-HPLC venom profiles were found to vary significantly between the two sampled localities, indicating geographic variability. Furthermore, the presence/absence of certain peaks in the venom chromatographic profiles appeared to be significantly correlated also to factors like body size and sex of the vipers. Our findings show that V. latastei venom is a variable phenotype. The intraspecific differences we detected in its composition likely mirror changes in the feeding ecology of this species, taking place during different life stages and under different environmental pressures. SIGNIFICANCE: Lataste's viper (Vipera latastei) is a medically important viper endemic to the Iberian Peninsula, inhabiting different habitats and undergoing a marked ontogenetic dietary shift. In the current study, we report the first proteomic analysis of V. latastei venom from two environmentally different localities in northern Portugal. Our bottom-up venomic analyses show that snake venom serine proteinases (SVSPs), phospholipases A2 (PLA2s), and C-type lectin-like (CTLs) proteins are the major components of adult V. latastei venom. The comparative analysis of young and adult venoms suggests the occurrence of ontogenetic shift in toxin abundances, with snake venom metalloproteinases (SVMPs) being the predominant toxins in juvenile venoms. Moreover, geographic venom variation between the two studied populations is also detected, with our statistical analyses suggesting that factors like body size and sex of the vipers are possibly at play in its determination. Our work represents the first assessment of the composition of V. latastei venom, and the first step towards a better understanding of the drivers behind its variability.

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.

Promoting co-existence between humans and venomous snakes through increasing the herpetological knowledge base

Snakebite incidence at least partly depends on the biology of the snakes involved. However, studies of snake biology have been largely neglected in favour of anthropic factors, with the exception of taxonomy, which has been recognised for some decades to affect the design of antivenoms. Despite this, within-species venom variation and the unpredictability of the correlation with antivenom cross-reactivity has continued to be problematic. Meanwhile, other aspects of snake biology, including behaviour, spatial ecology and activity patterns, distribution, and population demography, which can contribute to snakebite mitigation and prevention, remain underfunded and understudied. Here, we review the literature relevant to these aspects of snakebite and illustrate how demographic, spatial, and behavioural studies can improve our understanding of why snakebites occur and provide evidence for prevention strategies. We identify the large gaps that remain to be filled and urge that, in the future, data and relevant metadata be shared openly via public data repositories so that studies can be properly replicated and data used in future meta-analyses.

Intraspecific variability of the Central American rattlesnake (Crotalus simus) venom and its usefulness to obtain a representative standard venom

Snake venoms are mixtures of proteins whose physicochemical features confer them toxicity and immunogenicity. Animals (e.g., horses or sheep) immunized with snake venoms produce antibodies towards the venom proteins. Since these antibodies can neutralize the venom toxicity, they have been used to formulate snake antivenoms. The efficacy of the antivenoms is widely accepted, and standard venoms are expected to be representative of the snake's population that inhabit in the region where the antivenom is intended to be used. The representativeness of a single venom collected from a Crotalus simus snake, and its usefulness as standard venom to produce an antivenom is evaluated. The use of an "average venom" might be as representative of the population intended to be used, as the standard venom composed by many venom samples. Variations in the relative abundance concentration of crotoxin in the C. simus leads to different clinical manifestations, as well as differences in the neutralization efficacy of the antivenoms. A monovalent anti-Cs antivenom was produced from a single venom C. simus specimen, and its efficacy in neutralizing the lethal activity of 30 C. simus snakes was tested. Despite the variations in the relative abundance content of crotoxin found in the proteomes, the monovalent anti-Cs antivenom was successful in neutralize the toxicity caused by the variations on the venom composition of three different snake population used. Interestingly, it seems that the sex is not a key factor in the lethality of the venoms tested. The concept of representative venom mixtures for immunization should be revised for the case of C. simus on the populations found in Costa Rica, since it might use as less as one representative individual whose venom covers the mainly toxic enzymes.

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.

Venom of the Annulated Sea Snake Hydrophis cyanocinctus: A Biochemically Simple but Genetically Complex Weapon

Given that the venom system in sea snakes has a role in enhancing their secondary adaption to the marine environment, it follows that elucidating the diversity and function of venom toxins will help to understand the adaptive radiation of sea snakes. We performed proteomic and de novo NGS analyses to explore the diversity of venom toxins in the annulated sea snake (Hydrophis cyanocinctus) and estimated the adaptive molecular evolution of the toxin-coding unigenes and the toxicity of the major components. We found three-finger toxins (3-FTxs), phospholipase A2 (PLA2) and cysteine-rich secretory protein (CRISP) in the venom proteome and 59 toxin-coding unigenes belonging to 24 protein families in the venom-gland transcriptome; 3-FTx and PLA2 were the most abundant families. Nearly half of the toxin-coding unigenes had undergone positive selection. The short- (i.p. 0.09 μg/g) and long-chain neurotoxin (i.p. 0.14 μg/g) presented fairly high toxicity, whereas both basic and acidic PLA2s expressed low toxicity. The toxicity of H. cyanocinctus venom was largely determined by the 3-FTxs. Our data show the venom is used by H. cyanocinctus as a biochemically simple but genetically complex weapon and venom evolution in H. cyanocinctus is presumably driven by natural selection to deal with fast-moving prey and enemies in the marine environment.

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

Background

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

Results

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

Conclusions

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

Supplementary Information

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

Mutual enlightenment: A toolbox of concepts and methods for integrating evolutionary and clinical toxinology via snake venomics and the contextual stance

Snakebite envenoming is a neglected tropical disease that may claim over 100,000 human lives annually worldwide. Snakebite occurs as the result of an interaction between a human and a snake that elicits either a defensive response from the snake or, more rarely, a feeding response as the result of mistaken identity. Snakebite envenoming is therefore a biological and, more specifically, an ecological problem. Snake venom itself is often described as a "cocktail", as it is a heterogenous mixture of molecules including the toxins (which are typically proteinaceous) responsible for the pathophysiological consequences of envenoming. The primary function of venom in snake ecology is pre-subjugation, with defensive deployment of the secretion typically considered a secondary function. The particular composition of any given venom cocktail is shaped by evolutionary forces that include phylogenetic constraints associated with the snake's lineage and adaptive responses to the snake's ecological context, including the taxa it preys upon and by which it is predated upon. In the present article, we describe how conceptual frameworks from ecology and evolutionary biology can enter into a mutually enlightening relationship with clinical toxinology by enabling the consideration of snakebite envenoming from an "ecological stance". We detail the insights that may emerge from such a perspective and highlight the ways in which the high-fidelity descriptive knowledge emerging from applications of -omics era technologies - "venomics" and "antivenomics" - can combine with evolutionary explanations to deliver a detailed understanding of this multifactorial health crisis.

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

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

Old World Vipers-A Review about Snake Venom Proteomics of Viperinae and Their Variations

Fine-tuned by millions of years of evolution, snake venoms have frightened but also fascinated humanity and nowadays they constitute potential resources for drug development, therapeutics and antivenoms. The continuous progress of mass spectrometry techniques and latest advances in proteomics workflows enabled toxinologists to decipher venoms by modern omics technologies, so-called ‘venomics’. A tremendous upsurge reporting on snake venom proteomes could be observed. Within this review we focus on the highly venomous and widely distributed subfamily of Viperinae (Serpentes: Viperidae). A detailed public literature database search was performed (2003–2020) and we extensively reviewed all compositional venom studies of the so-called Old-World Vipers. In total, 54 studies resulted in 89 venom proteomes. The Viperinae venoms are dominated by four major, four secondary, six minor and several rare toxin families and peptides, respectively. The multitude of different venomics approaches complicates the comparison of venom composition datasets and therefore we differentiated between non-quantitative and three groups of quantitative workflows. The resulting direct comparisons within these groups show remarkable differences on the intra- and interspecies level across genera with a focus on regional differences. In summary, the present compilation is the first comprehensive up-to-date database on Viperinae venom proteomes and differentiating between analytical methods and workflows.

From birth to adulthood: An analysis of the Brazilian lancehead (Bothrops moojeni) venom at different life stages

The Brazilian lancehead (Bothrops moojeni) has a wide distribution in Brazil and represents a serious public health hazard. Previous works reported that the symptoms of snakebites caused by B. moojeni juveniles’ bites were mainly related to coagulation, while those caused by adults’ bites had a more prominent local damage. In this work, we analyzed the venoms of B. moojeni at different life stages to better understand the ontogeny shift in this species. Snakes were grouped by age and sex, and venom pools were formed accordingly. Compositional analyses by one-dimensional electrophoresis (1-DE), chromatography, and mass spectrometry revealed that ontogenetic changes might be mostly related to phospholipase A₂ (PLA₂) and metalloproteases. Regarding the venoms functional aspect, proteolytic, L-amino acid oxidase, PLA₂, and coagulant in vitro activities were assayed, but only the first and the last ones showed age-related changes, with the venom of snakes up to 1 year-old displaying lower proteolytic and higher coagulant activities, while those from 2 years-old onward presented the opposite relation. The venoms of 3 years-old snakes were exceptions to the compositional and functional pattern of adults as both venoms presented profiles similar to neonates. Sex-related differences were observed in specific groups and were not age-related. In vivo experiments (median lethal dose and hemorrhagic activity) were statistically similar between neonates and adults, however we verified that the adult venom killed mice faster comparing to the neonates. All venoms were mostly recognized by the antibothropic serum and displayed similar profiles to 1-DE in western blotting. In conclusion, the Brazilian lancehead venom showed ontogenetic shift in its composition and activities. Furthermore, this change occurred in snakes from 1 to 2 years-old, and interestingly the venom pools from 3 years-old snakes had particular characteristics, which highlights the importance of comprehensive studies to better understand venom variability.

Snake venom color and L-amino acid oxidase: An evidence of long-term captive Crotalus durissus terrificus venom plasticity

The venom color variation of Crotalus durissus terrificus (Cdt) is attributed to the presence of the toxin L-amino acid oxidase (LAAO). During the venom milking routine of Instituto Butantan, we have noticed that most venoms of captive Cdt specimens show a yellowish color, while most venoms of wild specimens are white. Here we describe a comparative analysis of long-term captive (LTC) and recently wild-caught (RWC) Cdt, focusing on LAAO variation. For the identification of LAAO in individual venoms, four different approaches were employed: evaluation of the enzymatic activity, SDS-PAGE, Western blotting, and ELISA. In addition, mass spectrometry analysis was performed using pooled samples. Although some variation among these methodologies was observed, it was possible to notice that the presence of LAAO was significantly higher in the venom of LTC individuals. LAAO was identified in 60-80% LTC specimens and in only 10-12% of RWC specimens. Furthermore, this enzyme accounts for 5.6% of total venom proteins of LTC Cdt pooled venom, while it corresponds to only 0.7% of RWC Cdt pooled venom. These findings strongly suggest that captive maintenance increases the expression of LAAO in Cdt venom.

The Tiger Rattlesnake genome reveals a complex genotype underlying a simple venom phenotype

Variation in gene regulation is ubiquitous, yet identifying the mechanisms producing such variation, especially for complex traits, is challenging. Snake venoms provide a model system for studying the phenotypic impacts of regulatory variation in complex traits because of their genetic tractability. Here, we sequence the genome of the Tiger Rattlesnake, which possesses the simplest and most toxic venom of any rattlesnake species, to determine whether the simple venom phenotype is the result of a simple genotype through gene loss or a complex genotype mediated through regulatory mechanisms. We generate the most contiguous snake-genome assembly to date and use this genome to show that gene loss, chromatin accessibility, and methylation levels all contribute to the production of the simplest, most toxic rattlesnake venom. We provide the most complete characterization of the venom gene-regulatory network to date and identify key mechanisms mediating phenotypic variation across a polygenic regulatory network.

A Wolf in Another Wolf's Clothing: Post-Genomic Regulation Dictates Venom Profiles of Medically-Important Cryptic Kraits in India

The Common Krait (Bungarus caeruleus) shares a distribution range with many other ‘phenotypically-similar’ kraits across the Indian subcontinent. Despite several reports of fatal envenomings by other Bungarus species, commercial Indian antivenoms are only manufactured against B. caeruleus. It is, therefore, imperative to understand the distribution of genetically distinct lineages of kraits, the compositional differences in their venoms, and the consequent impact of venom variation on the (pre)clinical effectiveness of antivenom therapy. To address this knowledge gap, we conducted phylogenetic and comparative venomics investigations of kraits in Southern and Western India. Phylogenetic reconstructions using mitochondrial markers revealed a new species of krait, Romulus’ krait (Bungarus romulusi sp. nov.), in Southern India. Additionally, we found that kraits with 17 mid-body dorsal scale rows in Western India do not represent a subspecies of the Sind Krait (B. sindanus walli) as previously believed, but are genetically very similar to B. sindanus in Pakistan. Furthermore, venom proteomics and comparative transcriptomics revealed completely contrasting venom profiles. While the venom gland transcriptomes of all three species were highly similar, venom proteomes and toxicity profiles differed significantly, suggesting the prominent role of post-genomic regulatory mechanisms in shaping the venoms of these cryptic kraits. In vitro venom recognition and in vivo neutralisation experiments revealed a strong negative impact of venom variability on the preclinical performance of commercial antivenoms. While the venom of B. caeruleus was neutralised as per the manufacturer’s claim, performance against the venoms of B. sindanus and B. romulusi was poor, highlighting the need for regionally-effective antivenoms in India.

Pinniped Ontogeny as a Window into the Comparative Physiology and Genomics of Hypoxia Tolerance

Diving physiology has received considerable scientific attention as it is a central element of the extreme phenotype of marine mammals. Many scientific discoveries have illuminated physiological mechanisms supporting diving, such as massive, internally bound oxygen stores and dramatic cardiovascular regulation. However, the cellular and molecular mechanisms that support the diving phenotype remain mostly unexplored as logistic and legal restrictions limit the extent of scientific manipulation possible. With next-generation sequencing (NGS) tools becoming more widespread and cost-effective, there are new opportunities to explore the diving phenotype. Genomic investigations come with their own challenges, particularly those including cross-species comparisons. Studying the regulatory pathways that underlie diving mammal ontogeny could provide a window into the comparative physiology of hypoxia tolerance. Specifically, in pinnipeds, which shift from terrestrial pups to elite diving adults, there is potential to characterize the transcriptional, epigenetic, and posttranslational differences between contrasting phenotypes while leveraging a common genome. Here we review the current literature detailing the maturation of the diving phenotype in pinnipeds, which has primarily been explored via biomarkers of metabolic capability including antioxidants, muscle fiber typing, and key aerobic and anaerobic metabolic enzymes. We also discuss how NGS tools have been leveraged to study phenotypic shifts within species through ontogeny, and how this approach may be applied to investigate the biochemical and physiological mechanisms that develop as pups become elite diving adults. We conclude with a specific example of the Antarctic Weddell seal by overlapping protein biomarkers with gene regulatory microRNA datasets.

Size Matters: An Evaluation of the Molecular Basis of Ontogenetic Modifications in the Composition of Bothrops jararacussu Snake Venom

Ontogenetic changes in venom composition have been described in Bothrops snakes, but only a few studies have attempted to identify the targeted paralogues or the molecular mechanisms involved in modifications of gene expression during ontogeny. In this study, we decoded B. jararacussu venom gland transcripts from six specimens of varying sizes and analyzed the variability in the composition of independent venom proteomes from 19 individuals. We identified 125 distinct putative toxin transcripts, and of these, 73 were detected in venom proteomes and only 10 were involved in the ontogenetic changes. Ontogenetic variability was linearly related to snake size and did not correspond to the maturation of the reproductive stage. Changes in the transcriptome were highly predictive of changes in the venom proteome. The basic myotoxic phospholipases A₂ (PLA₂s) were the most abundant components in larger snakes, while in venoms from smaller snakes, PIII-class SVMPs were the major components. The snake venom metalloproteinases (SVMPs) identified corresponded to novel sequences and conferred higher pro-coagulant and hemorrhagic functions to the venom of small snakes. The mechanisms modulating venom variability are predominantly related to transcriptional events and may consist of an advantage of higher hematotoxicity and more efficient predatory function in the venom from small snakes.

Gradual and Discrete Ontogenetic Shifts in Rattlesnake Venom Composition and Assessment of Hormonal and Ecological Correlates

Ontogenetic shifts in venom occur in many snakes but establishing their nature as gradual or discrete processes required additional study. We profiled shifts in venom expression from the neonate to adult sizes of two rattlesnake species, the eastern diamondback and the timber rattlesnake. We used serial sampling and venom chromatographic profiling to test if ontogenetic change occurs gradually or discretely. We found evidence for gradual shifts in overall venom composition in six of eight snakes, which sometimes spanned more than two years. Most chromatographic peaks shift gradually, but one quarter shift in a discrete fashion. Analysis of published diet data showed gradual shifts in overall diet composition across the range of body sizes attained by our eight study animals, while the shifts in abundance of different prey classes varied in form from gradual to discrete. Testosterone concentrations were correlated with the change in venom protein composition, but the relationship is not strong enough to suggest causation. Venom research employing simple juvenile versus adult size thresholds may be failing to account for continuous variation in venom composition lifespan. Our results imply that venom shifts represent adaptive matches to dietary shifts and highlight venom for studies of alternative gene regulatory mechanisms.

A brief review on the natural history, venomics and the medical importance of bushmaster (Lachesis) pit viper snakes

Snakes of the genus Lachesis, commonly known as bushmasters, are the largest venomous snakes in the Americas. Because these snakes have their habitats in areas of remote forests they are difficult to find, and consequently there are few studies of Lachesis taxa in their natural ecosystems. Bushmasters are distributed in tropical forest areas of South and Central America. In Brazil they can be found in the Amazon Rainforest and the Atlantic Forest. Despite the low incidence of cases, laquetic envenoming causes severe permanent sequelae due to the high amount of inoculated venom. These accidents are characterized by local pain, hemorrhage and myonecrosis that can be confused with bothropic envenomings. However, victims of Lachesis bites develop symptoms characteristic of Lachesis envenoming, known as vagal syndrome. An important message of this bibliographic synthesis exercise is that, despite having the proteomic profiles of all the taxa of the genus available, very few structure-function correlation studies have been carried out. Therefore the motivation for this review was to fill a gap in the literature on the genus Lachesis, about which there is no recent review. Here we discuss data scattered in a number of original articles published in specialized journals, spanning the evolutionary history and extant phylogeographic distribution of the bushmasters, their venom composition and diet, as well as the pathophysiology of their bites to humans and the biological activities and possible biotechnological applicability of their venom toxins.

Neurotropin Inhibits Lipid Accumulation by Maintaining Mitochondrial Function in Hepatocytes via AMPK Activation

The accumulation of lipid droplets in the cytoplasm of hepatocytes, known as hepatic steatosis, is a hallmark of non-alcoholic fatty liver disease (NAFLD). Inhibiting hepatic steatosis is suggested to be a therapeutic strategy for NAFLD. The present study investigated the actions of Neurotropin (NTP), a drug used for chronic pain in Japan and China, on lipid accumulation in hepatocytes as a possible treatment for NAFLD. NTP inhibited lipid accumulation induced by palmitate and linoleate, the two major hepatotoxic free fatty acids found in NAFLD livers. An RNA sequencing analysis revealed that NTP altered the expression of mitochondrial genes. NTP ameliorated palmitate-and linoleate-induced mitochondrial dysfunction by reversing mitochondrial membrane potential, respiration, and β-oxidation, suppressing mitochondrial oxidative stress, and enhancing mitochondrial turnover. Moreover, NTP increased the phosphorylation of AMPK, a critical factor in the regulation of mitochondrial function, and induced PGC-1β expression. Inhibition of AMPK activity and PGC-1β expression diminished the anti-steatotic effect of NTP in hepatocytes. JNK inhibition could also be associated with NTP-mediated inhibition of lipid accumulation, but we did not find the association between AMPK and JNK. These results suggest that NTP inhibits lipid accumulation by maintaining mitochondrial function in hepatocytes via AMPK activation, or by inhibiting JNK.

Toxin transcripts in Crotalus atrox venom and in silico structures of toxins

The western diamondback rattlesnake (Crotalus atrox) is a common and widespread North American pit viper species, and its venom possesses medical applications. In this research, we identified 14 of the most common transcripts encoding 11 major venom toxins including transcripts for a three-finger toxin (3FTx) from the crude venom of C. atrox. In silico three-dimensional (3D) structures of 9 venom toxins were predicted by using deduced toxin amino acid sequences and a computer programme-MODELLER. The accuracy of all predicted toxin structures was evaluated by five stereochemical structure parameters including discrete optimised protein energy (DOPE) score, root mean square deviation (RMSD), Z-score, overall quality factor (ERRAT), and φ/ψ dihedral angle distribution of toxin backbone Cα residues, resulting that the overall predicted models are satisfied quality evaluation checks. Our present toxin transcripts and simulated individual toxin structures are important not only for revealing species-specific venom gene expression profiles, but also for predicting the toxin-toxin interactions and designing the structure-based toxin inhibitors for the treatment of snakebites.

Comparative characterization of Viperidae snake venoms from Perú reveals two compositional patterns of phospholipase A2 expression

Snake species within the Bothrops complex (sensu lato) are of medical relevance in Latin America, but knowledge on their venom characteristics is limited, or even unavailable, for some taxa. Perú harbors 17 species of pit vipers, within the genera Bothrops, Bothriechis, Bothrocophias, Porthidium, Crotalus, and Lachesis. This study compared the venoms of twelve species, through chromatographic and electrophoretic profiles, as well as proteolytic and phospholipase A₂ (PLA₂) activities. Also, proteomic profiles were analyzed for nine of the venoms using a shotgun approach. Results unveiled conspicuous differences in the expression of venom PLA₂s among species, six of them presenting scarce levels as judged by RP-HPLC profiles. Since most species within the bothropoid lineage possess venoms with high to intermediate abundances of this protein family, our findings suggest the existence of a phenotypic duality in the expression of venom PLA₂s within the Bothrops (sensu lato) complex. Bothrops barnetti and Bothrocophias andianus venoms, very scarce in PLA₂s, were shown to lack significant myotoxic activity, highlighting that the observed variability in PLA₂ expression bears toxicological correlations with effects attributed to these proteins. Finally, an attempt to identify phylogenetic relationships of bothropoid species from Perú presenting low- or high-PLA₂ venom phenotypes showed an interspersed pattern, thus precluding a simple phylogenetic interpretation of this venom compositional dichotomy.

•

Venoms from 12 viperids of Perú were compared.

•

Conspicuous differences in the expression of PLA₂ were found.

•

Venoms presenting scarce levels of PLA₂ lack myotoxicity.

•

A new phenotypic dichotomy in venom PLA₂ expression is described within Bothrops (sensu lato).

Venom Systems as Models for Studying the Origin and Regulation of Evolutionary Novelties

A central goal in biology is to determine the ways in which evolution repeats itself. One of the most remarkable examples in nature of convergent evolutionary novelty is animal venom. Across diverse animal phyla, various specialized organs and anatomical structures have evolved from disparate developmental tissues to perform the same function, that is, produce and deliver a cocktail of potent molecules to subdue prey or predators. Venomous organisms therefore offer unique opportunities to investigate the evolutionary processes of convergence of key adaptive traits, and the molecular mechanisms underlying the emergence of novel genes, cells, and tissues. Indeed, some venomous species have already proven to be highly amenable as models for developmental studies, and recent work with venom gland organoids provides manipulatable systems for directly testing important evolutionary questions. Here, we provide a synthesis of the current knowledge that could serve as a starting point for the establishment of venom systems as new models for evolutionary and molecular biology. In particular, we highlight the potential of various venomous species for the study of cell differentiation and cell identity, and the regulatory dynamics of rapidly evolving, highly expressed, tissue-specific, gene paralogs. We hope that this review will encourage researchers to look beyond traditional study organisms and consider venom systems as useful tools to explore evolutionary novelties.

Functional venomics of the Big-4 snakes of Pakistan

In South Asia, the "Big-4" venomous snakes Naja naja, Bungarus caeruleus, Daboia russelii, and Echis carinatus are so-called because they are the most medically important snakes in the region. Antivenom is the only effective treatment option for snakebite envenoming but antivenom is not produced domestically in Pakistan making the country reliant on polyvalent products imported from India and Saudi Arabia. The present study investigated the toxin composition and activity of the venoms of Pakistani specimens by means of proteomic and physio/pharmacological experiments. To evaluate the composition of venoms, 1D/2D-PAGE of crude venoms and RP-HPLC followed by SDS-PAGE were performed. Enzymatic, hemolytic, coagulant and platelet aggregating activities of crude venoms were assayed and were concordant with expectations based on the abundance of protein species in each. Neutralization assays were performed using Bharat polyvalent antivenom (BPAV), a product raised against venoms from Big-4 specimens from southern India. BPAV exhibited cross-reactivity against the Pakistani venoms, however, neutralization of clinically relevant activities was variable and rarely complete. Cumulatively, the presented data not only highlight geographical variations present in the venoms of the Big-4 snakes of South Asia, but also demonstrate the neutralization potential of Indian polyvalent against the venom of Pakistani specimens. Given the partial neutralization observed, it is clear that whilst BPAV is a life-saving product in Pakistan, in future it is hoped that a region-specific product might be manufactured domestically, using venoms of local snakes in the immunising mixture.

Fangs for the Memories? A Survey of Pain in Snakebite Patients Does Not Support a Strong Role for Defense in the Evolution of Snake Venom Composition

Animals use venoms for multiple purposes, most prominently for prey acquisition and self-defense. In snakes, venom composition often evolves as a result of selection for optimization for local diet. However, whether selection for a defensive function has also played a role in driving the evolution of venom composition has remained largely unstudied. Here, we use an online survey of snakebite victims to test a key prediction of a defensive function, that envenoming should result in the rapid onset of severe pain. From the analysis of 584 snakebite reports, involving 192 species of venomous snake, we find that the vast majority of bites do not result in severe early pain. Phylogenetic comparative analysis shows that where early pain after a bite evolves, it is often lost rapidly. Our results, therefore, do not support the hypothesis that natural selection for antipredator defense played an important role in the origin of venom or front-fanged delivery systems in general, although there may be intriguing exceptions to this rule.

Too Many False Targets for MicroRNAs: Challenges and Pitfalls in Prediction of miRNA Targets and Their Gene Ontology in Model and Non-model Organisms

Short ("seed") or extended base pairing between microRNAs (miRNAs) and their target RNAs enables post-transcriptional silencing in many organisms. These interactions allow the computational prediction of potential targets. In model organisms, predicted targets are frequently validated experimentally; hence meaningful miRNA-regulated processes are reported. However, in non-models, these reports mostly rely on computational prediction alone. Many times, further bioinformatic analyses such as Gene Ontology (GO) enrichment are based on these in silico projections. Here such approaches are reviewed, their caveats are highlighted and the ease of picking false targets from predicted lists is demonstrated. Discoveries that shed new light on how miRNAs evolved to regulate targets in various phyletic groups are discussed, in addition to the pitfalls of target identification in non-model organisms. The goal is to prevent the misuse of bioinformatic tools, as they cannot bypass the biological understanding of miRNA-target regulation.

Morphological alterations caused by manual venom extraction on the main venom gland of Bothrops asper and Crotalus simus snakes (Serpentes: Viperidae): Long-term implications for antivenom production

The only scientifically validated treatment for snakebite envenomation is the administration of antivenoms. For their production, small quantities of snake venom are injected in animals to elicit a specific antibody response. Snakes are kept in captivity, and their venom is regularly extracted to assure antivenom access. It has already been reported that the pressure exerted upon the venom gland during this extraction can cause tissue damage and fibrosis, leading to a decrease in the venom yield. We described the histopathology of venom glands for B. asper and C. simus snakes used for antivenom production. Based on these reported tissue abnormalities, we quantify the tissue injury by a generated damage-SCORE and fibrosis. A variety of histopathological damages were found such as fibrosis, edema, necrosis, hemorrhage, and formation of anomalous structures, especially in C. simus, which is more prone to suffer severe damage. The level and severity of the damage depend on the frequency and the number of venom extractions. Furthermore, we design an experimental intensive venom extraction scheme with which we could confirm the causality of these effects. In addition to the histopathological damages, the LD₅₀ and biochemical venom composition were also affected giving experimental evidence that the venom extraction not only causes tissue damage but also affects the composition stability and toxicity of the venom. In order to produce quality and effective antivenoms, an improvement of the management of snake collections could be established, such as rotation groups to assure the quality of the venom yielded.

Tipping the Scales: The Migration-Selection Balance Leans toward Selection in Snake Venoms

The migration-selection interaction is the strongest determinant of whether a beneficial allele increases in frequency within a population. Results of empirical studies examining the role of gene flow in an adaptive context, however, have largely been equivocal, with examples of opposing outcomes being repeatedly documented (e.g., local adaptation with high levels of gene flow vs. gene swamping). We compared neutral genomic and venom expression divergence for three sympatric pit vipers with differing ecologies to determine if and how migration-selection disequilibria result in local adaptation. We specifically tested whether neutral differentiation predicted phenotypic differentiation within an isolation-by-distance framework. The decoupling of neutral and phenotypic differentiation would indicate selection led to adaptive divergence irrespective of migration, whereas a significant relationship between neutral and venom expression differentiation would provide evidence in favor of the constraining force of gene flow. Neutral differentiation and geographic distance predicted phenotypic differentiation only in the generalist species, indicating that selection was the predominant mechanism in the migration-selection balance underlying adaptive venom evolution in both specialists. Dispersal is thought to be a stronger influence on genetic differentiation than specialization, but our results suggest the opposite. Greater specialization may lead to greater diversification rates, and extreme spatial and temporal variation in selective pressures can favor generalist phenotypes evolving under strong stabilizing selection. Our results are consistent with these expectations, suggesting that the equivocal findings of studies examining the role of gene flow in an adaptive context may be explained by ecological specialization theory.

Evidence for Snake Venom Plasticity in a Long-Term Study with Individual Captive Bothrops atrox

Variability in snake venom composition has been frequently reported and correlated to the adaptability of snakes to environmental conditions. Previous studies report plasticity for the venom phenotype. However, these observations are not conclusive, as the results were based on pooled venoms, which present high individual variability. Here we tested the hypothesis of plasticity by influence of confinement and single diet type in the venom composition of 13 adult specimens of Bothrops atrox snakes, maintained under captivity for more than three years. Individual variability in venom composition was observed in samples extracted just after the capture of the snakes. However, composition was conserved in venoms periodically extracted from nine specimens, which presented low variability restricted to the less abundant components. In a second group, composed of four snakes, drastic changes were observed in the venom samples extracted at different periods, mostly related to snake venom metalloproteinases (SVMPs), the core function toxins of B. atrox venom, which occurred approximately between 400 and 500 days in captivity. These data show plasticity in the venom phenotype during the lifetime of adult snakes maintained under captive conditions. Causes or functional consequences involved in the phenotype modification require further investigations.

Intact protein mass spectrometry reveals intraspecies variations in venom composition of a local population of Vipera kaznakovi in Northeastern Turkey

We report on the variable venom composition of a population of the Caucasus viper (Vipera kaznakovi) in Northeastern Turkey. We applied a combination of venom gland transcriptomics, de-complexing bottom-up and top-down venomics. In contrast to sole bottom-up venomics approaches and gel or chromatography based venom comparison, our combined approach enables a faster and more detailed comparison of venom proteomes from multiple individuals. In total, we identified peptides and proteins from 15 toxin families, including snake venom metalloproteinases (svMP; 37.8%), phospholipases A2 (PLA2; 19.0%), snake venom serine proteinases (svSP; 11.5%), C-type lectins (CTL; 6.9%) and cysteine-rich secretory proteins (CRISP; 5.0%), in addition to several low abundant toxin families. Furthermore, we identified intraspecies variations of the venom composition of V. kaznakovi, and find these were mainly driven by the age of the animals, with lower svSP abundance detected in juveniles. On the proteoform level, several small molecular weight toxins between 5 and 8 kDa in size, as well as PLA2s, drove the differences observed between juvenile and adult individuals. This study provides novel insights into the venom variability of V. kaznakovi and highlights the utility of intact mass profiling for fast and detailed comparison of snake venom. BIOLOGICAL SIGNIFICANCE: Population level and ontogenetic venom variation (e.g. diet, habitat, sex or age) can result in a loss of antivenom efficacy against snakebites from wide ranging snake populations. The current state of the art for the analysis of snake venoms are de-complexing bottom-up proteomics approaches. While useful, these have the significant drawback of being time-consuming and following costly protocols, and consequently are often applied to pooled venom samples. To overcome these shortcomings and to enable rapid and detailed profiling of large numbers of individual venom samples, we integrated an intact protein analysis workflow into a transcriptomics-guided bottom-up approach. The application of this workflow to snake individuals of a local population of V. kaznakovi revealed intraspecies variations in venom composition, which are primarily explained by the age of the animals, and highlighted svSP abundance to be one of the molecular drivers for the compositional differences observed.