Functional characterizations of venom phenotypes in the eastern diamondback rattlesnake (Crotalus adamanteus) and evidence for expression-driven divergence in toxic activities among populations

Platelet and Red Cell Responses to Three North American Pit Vipers

We investigated the hemotoxic effects of three North American pit vipers in healthy human donor blood. Using experiments focusing on platelet and red blood cell activity, we found differential effects of these venoms on these cellular components. Platelet aggregation was most induced by C. adamanteus. Platelet activation was highest with C. atrox. Red blood cells had calcium expression and erythrocyte formation most induced by C. adamanteus and A. piscivorus. These results demonstrate the complex interplay of individual cellular effects with clinical presentations seen in envenomings from these species.

Venom diversity in Naja mossambica: Insights from proteomic and immunochemical analyses reveal intraspecific differences

BACKGROUND

Intraspecific variations in snake venom composition have been extensively documented, contributing to the diverse clinical effects observed in envenomed patients. Understanding these variations is essential for developing effective snakebite management strategies and targeted antivenom therapies. We aimed to comprehensively investigate venoms from three distinct populations of N. mossambica from Eswatini, Limpopo, and KwaZulu-Natal regions in Africa in terms of their protein composition and reactivity with three commercial antivenoms (SAIMR polyvalent, EchiTAb+ICP, and Antivipmyn Africa).

METHODOLOGY/PRINCIPAL FINDINGS

Naja mossambica venoms from Eswatini region exhibited the highest content of neurotoxic proteins, constituting 20.70% of all venom proteins, compared to Limpopo (13.91%) and KwaZulu-Natal (12.80%), and was characterized by the highest diversity of neurotoxic proteins, including neurotoxic 3FTxs, Kunitz-type inhibitors, vespryns, and mamba intestinal toxin 1. KwaZulu-Natal population exhibited considerably lower cytotoxic 3FTx, higher PLA2 content, and significant diversity in low-abundant proteins. Conversely, Limpopo venoms demonstrated the least diversity as demonstrated by electrophoretic and mass spectrometry analyses. Immunochemical assessments unveiled differences in venom-antivenom reactivity, particularly concerning low-abundance proteins. EchiTAb+ICP antivenom demonstrated superior reactivity in serial dilution ELISA assays compared to SAIMR polyvalent.

CONCLUSIONS/SIGNIFICANCE

Our findings reveal a substantial presence of neurotoxic proteins in N. mossambica venoms, challenging previous understandings of their composition. Additionally, the detection of numerous peptides aligning to uncharacterized proteins or proteins with unknown functions underscores a critical issue with existing venom protein databases, emphasizing the substantial gaps in our knowledge of snake venom protein components. This underscores the need for enhanced research in this domain. Moreover, our in vitro immunological assays suggest EchiTAb+ICP's potential as an alternative to SAIMR antivenom, requiring confirmation through prospective in vivo neutralization studies.

Snakes on a plain: biotic and abiotic factors determine venom compositional variation in a wide-ranging generalist rattlesnake

BACKGROUND

Snake venoms are trophic adaptations that represent an ideal model to examine the evolutionary factors that shape polymorphic traits under strong natural selection. Venom compositional variation is substantial within and among venomous snake species. However, the forces shaping this phenotypic complexity, as well as the potential integrated roles of biotic and abiotic factors, have received little attention. Here, we investigate geographic variation in venom composition in a wide-ranging rattlesnake (Crotalus viridis viridis) and contextualize this variation by investigating dietary, phylogenetic, and environmental variables that covary with venom.

RESULTS

Using shotgun proteomics, venom biochemical profiling, and lethality assays, we identify 2 distinct divergent phenotypes that characterize major axes of venom variation in this species: a myotoxin-rich phenotype and a snake venom metalloprotease (SVMP)-rich phenotype. We find that dietary availability and temperature-related abiotic factors are correlated with geographic trends in venom composition.

CONCLUSIONS

Our findings highlight the potential for snake venoms to vary extensively within species, for this variation to be driven by biotic and abiotic factors, and for the importance of integrating biotic and abiotic variation for understanding complex trait evolution. Links between venom variation and variation in biotic and abiotic factors indicate that venom variation likely results from substantial geographic variation in selection regimes that determine the efficacy of venom phenotypes across populations and snake species. Our results highlight the cascading influence of abiotic factors on biotic factors that ultimately shape venom phenotype, providing evidence for a central role of local selection as a key driver of venom variation.

Montane Rattlesnakes in México: Venoms of Crotalus tancitarensis and Related Species within the Crotalus intermedius Group

The Crotalus intermedius group is a clade of rattlesnakes consisting of several species adapted to a high elevation habitat, primarily in México. Crotalus tancitarensis was previously classified as C. intermedius, until individuals occurring on Cerro Tancítaro in Michoacán, México, were reevaluated and classified as a new species (C. tancitarensis) based on scale pattern and geographic location. This study aimed to characterize the venom of C. tancitarensis and compare the venom profile to those of other species within the Crotalus intermedius group using gel electrophoresis, biochemical assays, reverse-phase high performance liquid chromatography, mass spectrometry, and lethal toxicity (LD50) assays. Results show that the venom profiles of species within the Crotalus intermedius group are similar, but with distinct differences in phospholipase A2 (PLA2), metalloproteinase PI (SVMP PI), and kallikrein-like serine proteinase (SVSP) activity and relative abundance. Proteomic analysis indicated that the highland forms produce venoms with 50-60 protein isoforms and a composition typical of type I rattlesnake venoms (abundant SVMPs, lack of presynaptic PLA2-based neurotoxins), as well as a diversity of typical Crotalus venom components such as serine proteinases, PLA2s, C-type lectins, and less abundant toxins (LAAOs, CRiSPs, etc.). The overall venom profile of C. tancitarensis appears most similar to C. transversus, which is consistent with a previous mitochondrial DNA analysis of the Crotalus intermedius group. These rattlesnakes of the Mexican highlands represent a radiation of high elevation specialists, and in spite of divergence of species in these Sky Island habitats, venom composition of species analyzed here has remained relatively conserved. The majority of protein family isoforms are conserved in all members of the clade, and as seen in other more broadly distributed rattlesnake species, differences in their venoms are largely due to relative concentrations of specific components.

The Chemosensory Repertoire of the Eastern Diamondback Rattlesnake (Crotalus adamanteus) Reveals Complementary Genetics of Olfactory and Vomeronasal-Type Receptors

Pitviper sensory perception incorporates diverse stimuli through the integration of trichromatic color vision, bifocal heat-sensing, and dual-system chemoperception. Chemoperception, or olfaction, is mediated by chemoreceptors in the olfactory bulb and the vomeronasal organ, but the true genomic complexity of the gene families and their relative contributions is unknown. A full genomic accounting of pitviper chemoperception directly complements our current understanding of their venoms by generating a more complete polyphenic representation of their predatory arsenal. To characterize the genetic repertoire of pitviper chemoperception, we analyzed a full-genome assembly for Crotalus adamanteus, the eastern diamondback rattlesnake. We identified hundreds of genes encoding both olfactory receptors (ORs; 362 full-length genes) and type-2 vomeronasal receptors (V2Rs; 430 full-length genes). Many chemoreceptor genes are organized into large tandem repeat arrays. Comparative analysis of V2R orthologs across squamates demonstrates how gene array expansion and contraction underlies the evolution of the chemoreceptor repertoire, which likely reflects shifts in life history traits. Chromosomal assignments of chemosensory genes identified sex chromosome specific chemoreceptor genes, providing gene candidates underlying observed sex-specific chemosensory-based behaviors. We detected widespread episodic evolution in the extracellular, ligand-binding domains of both ORs and V2Rs, suggesting the diversification of chemoreceptors is driven by transient periods of positive selection. We provide a robust genetic framework for studying pitviper chemosensory ecology and evolution.

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.

Effect of dilution of canine blood samples on the specificity of saline agglutination tests for immune-mediated hemolysis

BACKGROUND

Saline agglutination tests (SATs) are widely recommended for diagnosis of immune-mediated hemolytic anemia in dogs, but there are frequent false-positive results.

OBJECTIVES

Specificity of SATs will improve at higher saline-to-blood ratios.

ANIMALS

One hundred fifty dogs treated at a veterinary referral hospital with hematocrits ≤30%.

METHODS

Prospective diagnostic accuracy study. Immune-mediated hemolysis (IMH) was considered present if a gel direct antiglobulin test (DAT) was positive and there was clinical evidence of hemolysis (n = 9), absent if another mechanism for anemia was identified and the DAT was negative or there was no hemolysis (n = 138), and if IMH status was unclear, dogs were excluded (n = 3). Saline agglutination tests were prepared at 1 : 1, 4 : 1, 9 : 1, and 49 : 1 saline-to-blood ratios, and microscopic agglutination was considered a positive result.

RESULTS

Specificity for IMH increased from 29% (95% confidence interval 20-38) at a 1 : 1 dilution to 97% (93-99) at a 49 : 1 dilution. Sensitivity was 88% (47-100) at 1 : 1 and 4 : 1 dilutions and 67% (30-93%) at 9 : 1 and 49 : 1 dilutions. Diagnostic accuracy increased from 33% (24-42) at 1 : 1 dilution to 95% (90-98) at 49 : 1 dilution.

CONCLUSIONS AND CLINICAL IMPORTANCE

If performed using a 49 : 1 saline-to-blood ratio, SATs achieve high specificity for IMH. Based on a gold standard of positive DAT and evidence of hemolysis, lower saline-to-blood ratio results should not be used because false-positive results are common.

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.

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.

Convergent recruitment of adamalysin-like metalloproteases in the venom of the red bark centipede (Scolopocryptops sexspinosus)

Many venom proteins have presumably been convergently recruited by taxa from diverse venomous lineages. These toxic proteins have characteristics that allow them to remain stable in solution and have a high propensity for toxic effects on prey and/or potential predators. Despite this well-established convergent toxin recruitment, some toxins seem to be lineage specific. To further investigate the toxic proteins found throughout venomous lineages, venom proteomics and venom-gland transcriptomics were performed on two individual red bark centipedes (Scolopocryptops sexspinosus). Combining the protein phenotype with the transcript genotype resulted in the first in-depth venom characterization of S. sexspinosus, including 72 venom components that were identified in both the transcriptome and proteome and 1468 nontoxin transcripts identified in the transcriptome. Ten different toxin families were represented in the venom and venom gland with the majority of the toxins belonging to metalloproteases, CAPS (cysteine-rich secretory protein, antigen 5, and pathogenesis-related 1 proteins), and β-pore-forming toxins. Nine of these toxin families shared a similar proteomic structure to venom proteins previously identified from other centipedes. However, the most highly expressed toxin family, the adamalysin-like metalloproteases, has until now only been observed in the venom of snakes. We confirmed adamalysin-like metalloprotease activity by means of in vivo functional assays. The recruitment of an adamalysin-like metalloprotease into centipede venom represents a striking case of convergent evolution.

A Single Mutation Unlocks Cascading Exaptations in the Origin of a Potent Pitviper Neurotoxin

Evolutionary innovations and complex phenotypes seemingly require an improbable amount of genetic change to evolve. Rattlesnakes display two dramatically different venom phenotypes. Type I venoms are hemorrhagic with low systemic toxicity and high expression of tissue-destroying snake venom metalloproteinases. Type II venoms are highly neurotoxic and lack snake venom metalloproteinase expression and associated hemorrhagic activity. This dichotomy hinges on Mojave toxin (MTx), a phospholipase A2 (PLA2) based β-neurotoxin expressed in Type II venoms. MTx is comprised of a nontoxic acidic subunit that undergoes extensive proteolytic processing and allosterically regulates activity of a neurotoxic basic subunit. Evolution of the acidic subunit presents an evolutionary challenge because the need for high expression of a nontoxic venom component and the proteolytic machinery required for processing suggests genetic changes of seemingly little immediate benefit to fitness. We showed that MTx evolved through a cascading series of exaptations unlocked by a single nucleotide change. The evolution of one new cleavage site in the acidic subunit unmasked buried cleavage sites already present in ancestral PLA2s, enabling proteolytic processing. Snake venom serine proteases, already present in the venom to disrupt prey hemostasis, possess the requisite specificities for MTx acidic subunit proteolysis. The dimerization interface between MTx subunits evolved by exploiting a latent, but masked, hydrophobic interaction between ancestral PLA2s. The evolution of MTx through exaptation of existing functional and structural features suggests complex phenotypes that depend on evolutionary innovations can arise from minimal genetic change enabled by prior evolution.

Exploring the Diversity and Novelty of Toxin Genes in Naja sumatrana, the Equatorial Spitting Cobra from Malaysia through De Novo Venom-Gland Transcriptomics

The equatorial spitting cobra, Naja sumatrana, is a distinct species of medically important venomous snakes, listed as WHO Category 1 in Southeast Asia. The diversity of its venom genes has not been comprehensively examined, although a few toxin sequences annotated to Naja sputatrix were reported previously through cloning studies. To investigate this species venom genes’ diversity, de novo venom-gland transcriptomics of N. sumatrana from West Malaysia was conducted using next-generation sequencing technology. Genes encoding toxins represented only 60 of the 55,396 transcripts, but were highly expressed, contributing to 79.22% of total gene expression (by total FPKM) in the venom-glands. The toxin transcripts belong to 21 families, and 29 transcripts were further identified as full-length. Three-finger toxins (3FTx) composed of long, short, and non-conventional groups, constituted the majority of toxin transcripts (91.11% of total toxin FPKM), followed by phospholipase A₂ (PLA₂, 7.42%)—which are putatively pro-inflammatory and cytotoxic. The remaining transcripts in the 19 families were expressed at extremely low levels. Presumably, these toxins were associated with ancillary functions. Our findings unveil the diverse toxin genes unique to N. sumatrana, and provide insights into the pathophysiology of N. sumatrana envenoming.

Comparative venom-gland transcriptomics and venom proteomics of four Sidewinder Rattlesnake (Crotalus cerastes) lineages reveal little differential expression despite individual variation

Changes in gene expression can rapidly influence adaptive traits in the early stages of lineage diversification. Venom is an adaptive trait comprised of numerous toxins used for prey capture and defense. Snake venoms can vary widely between conspecific populations, but the influence of lineage diversification on such compositional differences are unknown. To explore venom differentiation in the early stages of lineage diversification, we used RNA-seq and mass spectrometry to characterize Sidewinder Rattlesnake (Crotalus cerastes) venom. We generated the first venom-gland transcriptomes and complementary venom proteomes for eight individuals collected across the United States and tested for expression differences across life history traits and between subspecific, mitochondrial, and phylotranscriptomic hypotheses. Sidewinder venom was comprised primarily of hemorrhagic toxins, with few cases of differential expression attributable to life history or lineage hypotheses. However, phylotranscriptomic lineage comparisons more than doubled instances of significant expression differences compared to all other factors. Nevertheless, only 6.4% of toxins were differentially expressed overall, suggesting that shallow divergence has not led to major changes in Sidewinder venom composition. Our results demonstrate the need for consensus venom-gland transcriptomes based on multiple individuals and highlight the potential for discrepancies in differential expression between different phylogenetic hypotheses.

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.

Quantity, Not Quality: Rapid Adaptation in a Polygenic Trait Proceeded Exclusively through Expression Differentiation

A trait's genomic architecture can affect the rate and mechanism of adaptation, and although many ecologically-important traits are polygenic, most studies connecting genotype, phenotype, and fitness in natural populations have focused on traits with relatively simple genetic bases. To understand the genetic basis of polygenic adaptation, we must integrate genomics, phenotypic data, ecology, and fitness effects for a genetically tractable, polygenic trait; snake venoms provide such a system for studying polygenic adaptation because of their genetic tractability and vital ecological role in feeding and defense. We used a venom transcriptome-proteome map, quantitative proteomics, genomics, and fitness assays in sympatric prey to construct a genotype-phenotype-fitness map for the venoms of an island-mainland pair of rattlesnake populations. Reciprocal fitness experiments demonstrated that each population was locally adapted to sympatric prey. We identified significant expression differentiation with little to no coding-sequence variation across populations, demonstrating that expression differentiation was exclusively the genetic basis of polygenic adaptation. Previous research on the genetics of adaptation, however, has largely been biased toward investigating protein-coding regions because of the complexity of gene regulation. Our results showed that biases at the molecular level can be in the opposite direction, highlighting the need for more systematic comparisons of different molecular mechanisms underlying rapid, adaptive evolution in polygenic traits.

Biological and Proteolytic Variation in the Venom of Crotalus scutulatus scutulatus from Mexico

Rattlesnake venoms may be classified according to the presence/absence and relative abundance of the neurotoxic phospholipases A2s (PLA2s), such as Mojave toxin, and snake venom metalloproteinases (SVMPs). In Mexico, studies to determine venom variation in Mojave Rattlesnakes (Crotalus scutulatus scutulatus) are limited and little is known about the biological and proteolytic activities in this species. Tissue (34) and venom (29) samples were obtained from C. s. scutulatus from different locations within their distribution in Mexico. Mojave toxin detection was carried out at the genomic (by PCR) and protein (by ELISA) levels for all tissue and venom samples. Biological activity was tested on representative venoms by measuring LD50 and hemorrhagic activity. To determine the approximate amount of SVMPs, 15 venoms were separated by RP-HPLC and variation in protein profile and proteolytic activity was evaluated by SDS-PAGE (n = 28) and Hide Powder Azure proteolytic analysis (n = 27). Three types of venom were identified in Mexico which is comparable to the intraspecific venom diversity observed in the Sonoran Desert of Arizona, USA: Venom Type A (∼Type II), with Mojave toxin, highly toxic, lacking hemorrhagic activity, and with scarce proteolytic activity; Type B (∼Type I), without Mojave toxin, less toxic than Type A, highly hemorrhagic and proteolytic; and Type A + B, containing Mojave toxin, as toxic as venom Type A, variable in hemorrhagic activity and with intermediate proteolytic activity. We also detected a positive correlation between SVMP abundance and hemorrhagic and proteolytic activities. Although more sampling is necessary, our results suggest that venoms containing Mojave toxin and venom lacking this toxin are distributed in the northwest and southeast portions of the distribution in Mexico, respectively, while an intergradation in the middle of both zones is present.

Catch a tiger snake by its tail: Differential toxicity, co-factor dependence and antivenom efficacy in a procoagulant clade of Australian venomous snakes

A paradigm of venom research is adaptive evolution of toxins as part of a predator-prey chemical arms race. This study examined differential co-factor dependence, variations relative to dietary preference, and the impact upon relative neutralisation by antivenom of the procoagulant toxins in the venoms of a clade of Australian snakes. All genera were characterised by venoms rich in factor Xa which act upon endogenous prothrombin. Examination of toxin sequences revealed an extraordinary level of conservation, which indicates that adaptive evolution is not a feature of this toxin type. Consistent with this, the venoms did not display differences on the plasma of different taxa. Examination of the prothrombin target revealed endogenous blood proteins are under extreme negative selection pressure for diversification, this in turn puts a strong negative selection pressure upon the toxins as sequence diversification could result in a drift away from the target. Thus this study reveals that adaptive evolution is not a consistent feature in toxin evolution in cases where the target is under negative selection pressure for diversification. Consistent with this high level of toxin conservation, the antivenom showed extremely high-levels of cross-reactivity. There was however a strong statistical correlation between relative degree of phospholipid-dependence and clotting time, with the least dependent venoms producing faster clotting times than the other venoms even in the presence of phospholipid. The results of this study are not only of interest to evolutionary and ecological disciplines, but also have implications for clinical toxinology.

Integrated Venomics and Venom Gland Transcriptome Analysis of Juvenile and Adult Mexican Rattlesnakes Crotalus simus, C. tzabcan, and C. culminatus Revealed miRNA-modulated Ontogenetic Shifts

Adult rattlesnakes within genus Crotalus express one of two distinct venom phenotypes, type I (hemorrhagic) and type II (neurotoxic). In Costa Rican Central American rattlesnake, ontogenetic changes in the concentration of miRNAs modulate venom type II to type I transition. Venomics and venom gland transcriptome analyses showed that adult C. simus and C. tzabcan expressed intermediate patterns between type II and type I venoms, whereas C. culminatus had a canonical type I venom. Neonate/juvenile and adult Mexican rattlesnakes showed notable inter- and intraspecific variability in the number, type, abundance and ontogenetic shifts of the transcriptional and translational venom gland activities. These results support a role for miRNAs in the ontogenetic venom compositional changes in the three congeneric Mexican rattlesnakes. It is worth noting the finding of dual-action miRNAs, which silence the translation of neurotoxic heterodimeric PLA₂ crotoxin and acidic PLA₂ mRNAs while simultaneously up-regulating SVMP-targeting mRNAs. Dual transcriptional regulation potentially explains the existence of mutually exclusive crotoxin-rich (type-II) and SVMP-rich (type-I) venom phenotypic dichotomy among rattlesnakes. Our results support the hypothesis that alterations of the distribution of miRNAs, modulating the translational activity of venom gland toxin-encoding mRNAs in response to an external cue, may contribute to the mechanism generating adaptive venom variability.

Good vibrations: Assessing the stability of snake venom composition after researcher-induced disturbance in the laboratory

Phenotypic plasticity contributes to intraspecific variation in traits of many animal species. Venom is an integral trait to the success and survival of many snake species, and potential plasticity in venom composition is important to account for in the context of basic research as well as in human medicine for treating the various symptoms of snakebite and producing effective anti-venoms. Researchers may unknowingly induce changes in venom variation by subjecting snakes to novel disturbances and potential stressors. We explored phenotypic plasticity in snake venom composition over time in captive Pacific rattlesnakes (Crotalus oreganus) exposed to vibration treatment, compared to an undisturbed control group. Venom composition did not change significantly in response to vibration, nor was there a detectable effect of overall time in captivity, even though snakes re-synthesized venom stores while subjected to novel disturbance in the laboratory. This result indicates that venom composition is a highly repeatable phenotype over short time spans and that the composition of venom within adult individuals may be resistant to or unaffected by researcher-induced disturbance. On the other hand, the change in venom composition, measured as movement along the first principle component of venom phenotype space, was associated with baseline corticosterone (CORT) levels in the snakes. While differential forms of researcher-induced disturbance may not affect venom composition, significant changes in baseline CORT, or chronic stress, may affect the venom phenotype, and further investigations will be necessary to assess the nature of the relationship between CORT and venom protein expression.

Selection To Increase Expression, Not Sequence Diversity, Precedes Gene Family Origin and Expansion in Rattlesnake Venom

Gene duplication is the primary mechanism leading to new genes and phenotypic novelty, but the proximate evolutionary processes underlying gene family origin, maintenance, and expansion are poorly understood. Although sub- and neofunctionalization provide clear long-term advantages, selection does not act with foresight, and unless a redundant gene copy provides an immediate fitness advantage, the copy will most likely be lost. Many models for the evolution of genes immediately following duplication have been proposed, but the robustness and applicability of these models is unclear because of the lack of data at the population level. We used qPCR, protein expression data, genome sequencing, and hybrid enrichment to test three competing models that differ in whether selection favoring the spread of duplicates acts primarily on expression level or sequence diversity for specific toxin-encoding loci in the eastern diamondback rattlesnake (Crotalus adamanteus). We sampled 178 individuals and identified significant inter- and intrapopulation variation in copy number, demonstrated that copy number was significantly and positively correlated with protein expression, and found little to no sequence variation across paralogs in all populations. Collectively, these results demonstrate that selection for increased expression, not sequence diversity, was the proximate evolutionary process underlying gene family origin and expansion, providing data needed to resolve the debate over which evolutionary processes govern the fates of gene copies immediately following duplication.

The genetics of venom ontogeny in the eastern diamondback rattlesnake (Crotalus adamanteus)

The same selective forces that give rise to rapid inter- and intraspecific divergence in snake venoms can also favor differences in venoms across life-history stages. Ontogenetic changes in venom composition are well known and widespread in snakes but have not been investigated to the level of unambiguously identifying the specific loci involved. The eastern diamondback rattlesnake was previously shown to undergo an ontogenetic shift in venom composition at sexual maturity, and this shift accounted for more venom variation than geography. To characterize the genetics underlying the ontogenetic venom compositional change in C. adamanteus, we sequenced adult/juvenile pairs of venom-gland transcriptomes from five populations previously shown to have different adult venom compositions. We identified a total of 59 putative toxin transcripts for C. adamanteus, and 12 of these were involved in the ontogenetic change. Three toxins were downregulated, and nine were upregulated in adults relative to juveniles. Adults and juveniles expressed similar total levels of snake-venom metalloproteinases but differed substantially in their featured paralogs, and adults expressed higher levels of Bradykinin-potentiating and C-type natriuretic peptides, nerve growth factor, and specific paralogs of phospholipases A₂ and snake venom serine proteinases. Juvenile venom was more toxic to mice, indicating that the expression differences resulted in a phenotypically, and therefore potentially ecologically, significant difference in venom function. We also showed that adult and juvenile venom-gland transcriptomes for a species with known ontogenetic venom variation were equally effective at individually providing a full characterization of the venom genes of a species but that any particular individual was likely to lack several toxins in their transcriptome. A full characterization of a species’ venom-gene complement therefore requires sequencing more than one individual, although the ages of the individuals are unimportant.

Comparative venom gland transcriptomics of Naja kaouthia (monocled cobra) from Malaysia and Thailand: elucidating geographical venom variation and insights into sequence novelty

Background

The monocled cobra (Naja kaouthia) is a medically important venomous snake in Southeast Asia. Its venom has been shown to vary geographically in relation to venom composition and neurotoxic activity, indicating vast diversity of the toxin genes within the species. To investigate the polygenic trait of the venom and its locale-specific variation, we profiled and compared the venom gland transcriptomes of N. kaouthia from Malaysia (NK-M) and Thailand (NK-T) applying next-generation sequencing (NGS) technology.

Methods

The transcriptomes were sequenced on the Illumina HiSeq platform, assembled and followed by transcript clustering and annotations for gene expression and function. Pairwise or multiple sequence alignments were conducted on the toxin genes expressed. Substitution rates were studied for the major toxins co-expressed in NK-M and NK-T.

Results and discussion

The toxin transcripts showed high redundancy (41–82% of the total mRNA expression) and comprised 23 gene families expressed in NK-M and NK-T, respectively (22 gene families were co-expressed). Among the venom genes, three-finger toxins (3FTxs) predominated in the expression, with multiple sequences noted. Comparative analysis and selection study revealed that 3FTxs are genetically conserved between the geographical specimens whilst demonstrating distinct differential expression patterns, implying gene up-regulation for selected principal toxins, or alternatively, enhanced transcript degradation or lack of transcription of certain traits. One of the striking features that elucidates the inter-geographical venom variation is the up-regulation of α-neurotoxins (constitutes ∼80.0% of toxin’s fragments per kilobase of exon model per million mapped reads (FPKM)), particularly the long-chain α-elapitoxin-Nk2a (48.3%) in NK-T but only 1.7% was noted in NK-M. Instead, short neurotoxin isoforms were up-regulated in NK-M (46.4%). Another distinct transcriptional pattern observed is the exclusively and abundantly expressed cytotoxin CTX-3 in NK-T. The findings suggested correlation with the geographical variation in proteome and toxicity of the venom, and support the call for optimising antivenom production and use in the region. Besides, the current study uncovered full and partial sequences of numerous toxin genes from N. kaouthia which have not been reported hitherto; these include N. kaouthia-specific l-amino acid oxidase (LAAO), snake venom serine protease (SVSP), cystatin, acetylcholinesterase (AChE), hyaluronidase (HYA), waprin, phospholipase B (PLB), aminopeptidase (AP), neprilysin, etc. Taken together, the findings further enrich the snake toxin database and provide deeper insights into the genetic diversity of cobra venom toxins.

Venomics: integrative venom proteomics and beyond*

Venoms are integrated phenotypes that evolved independently in, and are used for predatory and defensive purposes by, a wide phylogenetic range of organisms. The same principles that contribute to the evolutionary success of venoms, contribute to making the study of venoms of great interest in such diverse fields as evolutionary ecology and biotechnology. Evolution is profoundly contingent, and nature also reinvents itself continuosly. Changes in a complex phenotypic trait, such as venom, reflect the influences of prior evolutionary history, chance events, and selection. Reconstructing the natural history of venoms, particularly those of snakes, which will be dealt with in more detail in this review, requires the integration of different levels of knowledge into a meaningful and comprehensive evolutionary framework for separating stochastic changes from adaptive evolution. The application of omics technologies and other disciplines have contributed to a qualitative and quantitative advance in the road map towards this goal. In this review we will make a foray into the world of animal venoms, discuss synergies and complementarities of the different approaches used in their study, and identify current bottlenecks that prevent inferring the evolutionary mechanisms and ecological constraints that molded snake venoms to their present-day variability landscape.