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

Venom variation and ontogenetic changes in the Crotalus molossus complex: Insights into composition, activities, and antivenom neutralization

The Crotalus molossus complex consists of five to seven phylogenetically related lineages of black-tailed rattlesnakes widely distributed in Mexico. While previous studies have noted venom variation within specific lineages of the Crotalus molossus complex, a comprehensive characterization of interspecific and ontogenetic venom variations, their functional implications, and the neutralizing ability of the Mexican antivenom against these variants remains largely unexamined. Herein, using two proteomic approaches for five lineages (C. basiliscus, C. m. molossus, C. m. nigrescens, C. m. oaxacus, and C. ornatus) of the C. molossus complex we characterized the number of toxins and their relative abundance in the venom of individuals of varying sizes. All five lineages undergo ontogenetic venom composition shifts associated with snake length. However, the pattern of ontogenetic shifts varied among lineages. In some lineages, these shifts led to significant differences in proteolytic, phospholipase A2, and fibrinogenolytic activities. Venom in smaller C. basiliscus, C. m. nigrescens, and C. m. oaxacus individuals had lower LD50 values (more lethal) in mice. Whereas the venom lethality of C. m. nigrescens (both juvenile and adult) and C. m. oaxacus (adult) was several times higher in a mammalian (mouse) model than in a reptilian (iguana) model. Antivipmyn® showed different neutralizing potencies toward venom pools. Overall, our results indicated that even among closely related rattlesnake lineages, venom phenotypes may vary greatly, impacting their function and the efficacy of antivenom neutralization.

Phylogenomic Discordance is Driven by Wide-Spread Introgression and Incomplete Lineage Sorting During Rapid Species Diversification Within Rattlesnakes (Viperidae: Crotalus and Sistrurus)

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

Where the "ruber" Meets the Road: Using the Genome of the Red Diamond Rattlesnake to Unravel the Evolutionary Processes Driving Venom Evolution

Understanding the proximate and ultimate causes of phenotypic variation is fundamental in evolutionary research, as such variation provides the substrate for selection to act upon. Although trait variation can arise due to selection, the importance of neutral processes is sometimes understudied. We presented the first reference-quality genome of the Red Diamond Rattlesnake (Crotalus ruber) and used range-wide 'omic data to estimate the degree to which neutral and adaptive evolutionary processes shaped venom evolution. We characterized population structure and found substantial genetic differentiation across two populations, each with distinct demographic histories. We identified significant differentiation in venom expression across age classes with substantially reduced but discernible differentiation across populations. We then used conditional redundancy analysis to test whether venom expression variation was best predicted by neutral divergence patterns or geographically variable (a)biotic factors. Snake size was the most significant predictor of venom variation, with environment, prey availability, and neutral sequence variation also identified as significant factors, though to a lesser degree. By directly including neutrality in the model, our results confidently highlight the predominant, yet not singular, role of life history in shaping venom evolution.

Fab Antivenom Reversal of Neurotoxicity Caused by a Juvenile Crotalus horridus Lacking Canebrake Toxin

We present a case of neurotoxic effects in a pediatric patient after envenomation by a timber rattlesnake (Crotalus horridus) in the Appalachian upstate of South Carolina. Though some members of this species are capable of primarily neurotoxic envenomation, there is heterogeneity in venom composition, and neurotoxic timber rattlesnakes are not endemic to the Appalachian region. However, neurotoxic effects caused by C horridus species lacking typical neurotoxins have been suspected, though not previously confirmed in the medical literature. This case presents a patient who was envenomated by a genotypically confirmed non-neurotoxic C horridus but who nevertheless presented with symptoms consistent with primary neurotoxicity. Neurotoxic effects can be variable in their response to traditional antivenom, though this patient demonstrated rapid response to treatment, representing a novel case in the literature of neurotoxic effects from a snake lacking typical neurotoxins with documented improvement with traditional antivenom.

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

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

Neurologic and hematologic sequalae following a timber rattlesnake (Crotalus horridus) envenomation in a dachshund

A 2-year-old female Dachshund had a witnessed timber rattlesnake envenomation. Although rattlesnake envenomations are a common, potentially life-threatening event in companion animals, timber rattlesnake envenomations in the dog are rarely reported. This dog described in this case report had significant hematologic and neurologic clinical derangements consistent with Types A and B rattlesnake venom and a suspected hypersensitivity reaction to the venom. This patient was treated aggressively with antivenom and fully recovered without any persistent neurologic signs at follow-up.

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

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

A current perspective on snake venom composition and constituent protein families

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