Diversification in vipers: Phylogenetic relationships, time of divergence and shifts in speciation rates

Paraspecificity of Mexican antivipmyn TRI antivenom in envenomation by Chinese Protobothrops mangshanensis (Mangshan pit viper) in France: A case report and experimental neutralization of venom procoagulant effect

This case report presents an exotic envenomation by a Chinese snake, Protobothrops mangshanensis. Its venom exhibited potent activity against plasma and fibrinogen, among other enzymatic activities. The patient initially presented with edema of the right upper limb, without tissue necrosis. There were no signs of bleeding; however, severe hypofibrinogenemia was observed (nadir value at 0.4 g/L), with a marked increase in fibrinogen degradation products and D-dimers, without any other coagulation disturbances. In the absence of a specific antivenom available against Asian Crotalinae venoms, the patient was treated at the 29th hour after bite with six vials of Antivipmyn™ TRI (Instituto Bioclon, Mexico, Mexico), a Mexican antivenom initially intended for American Crotalinae venoms, i.e., Bothrops asper, Lachesis muta and Crotalus durissus. Fibrinogen began to rise 6 hours after the antivenom infusion and was within the normal range 38 hours later. The report also underscores the utility of ClotPro® (Haemonetics ®USA), a viscoelastic test, for real-time monitoring of the snakebite-related coagulopathy. The clotting time was extended to 188 seconds on the EX-test while the MCF was decreased to 31 mm on the EX-test and the AP-test and was not measurable on the FIB-test, confirming severe hypofibrinogenemia. In order to confirm the paraspecificity of antivenom on the venom of P. mangshanensis, we studied the experimental neutralization of the venom procoagulant effect by Antivipmyn TRI and Green Pit Viper antivenom, which has been used in previous published clinical cases of P. mangshanensis envenomation. Both Antivipmyn™ TRI and Green Pit Viper antivenom corrected the procoagulant effect induced by P. mangshanensis venom. These findings suggest that Antivipmyn™ TRI cross-reacts with Protobothrops mangshanensis venom. In the absence of antivenom covering Asian Crotalinae, Antivipmyn TRI should be considered to treat an envenomation by Protobothrops spp.

Phylogenomic insights into the diversity and evolution of Palearctic vipers

Despite decades of molecular research, phylogenetic relationships in Palearctic vipers (genus Vipera) still essentially rely on a few loci, such as mitochondrial barcoding genes. Here we examined the diversity and evolution of Vipera with ddRAD-seq data from 33 representative species and subspecies. Phylogenomic analyses of ∼ 1.1 Mb recovered nine major clades corresponding to known species/species complexes which are generally consistent with the mitochondrial phylogeny, albeit with a few deep discrepancies that highlight past hybridization events. The most spectacular case is the Italian-endemic V. walser, which is grouped with the alpine genetic diversity of V. berus in the nuclear tree despite carrying a divergent mitogenome related to the Caucasian V. kaznakovi complex. Clustering analyses of SNPs suggest potential admixture between diverged Iberian taxa (V. aspis zinnikeri and V. seoanei), and confirm that the Anatolian V. pontica corresponds to occasional hybrids between V. (ammodytes) meridionalis and V. kaznakovi. Finally, all analyzed lineages of the V. berus complex (including V. walser and V. barani) form vast areas of admixture and may be delimited as subspecies. Our study sets grounds for future taxonomic and phylogeographic surveys on Palearctic vipers, a group of prime interest for toxinological, ecological, biogeographic and conservation research.

Novel phylogenomic inference and 'Out of Asia' biogeography of cobras, coral snakes and their allies

Estimation of evolutionary relationships among lineages that rapidly diversified can be challenging, and, in such instances, inaccurate or unresolved phylogenetic estimates can lead to erroneous conclusions regarding historical geographical ranges of lineages. One example underscoring this issue has been the historical challenge posed by untangling the biogeographic origin of elapoid snakes, which includes numerous dangerously venomous species as well as species not known to be dangerous to humans. The worldwide distribution of this lineage makes it an ideal group for testing hypotheses related to historical faunal exchanges among the many continents and other landmasses occupied by contemporary elapoid species. We developed a novel suite of genomic resources, included worldwide sampling, and inferred a robust estimate of evolutionary relationships, which we leveraged to quantitatively estimate geographical range evolution through the deep-time history of this remarkable radiation. Our phylogenetic and biogeographical estimates of historical ranges definitively reject a lingering former 'Out of Africa' hypothesis and support an 'Out of Asia' scenario involving multiple faunal exchanges between Asia, Africa, Australasia, the Americas and Europe.

Successful use of anti-venom cross-neutralization effects in the clinical management of Shore Pit Viper envenomation

As the landscape becomes more urbanized, snakebites have increasingly become uncommon presentations to the emergency departments in Singapore, while snakebites causing significant envenomation are even rarer. In this case report, we discuss a 55-year-old man who had significant envenomation from a Shore Pit Viper (Trimeresurus Purpureomaculatus) and who was successfully treated with haemato-toxic polyvalent antivenom (HPAV). He initially presented with pain, swelling and bleeding over his wound. Due to a deterioration in his coagulation profile, he was given two doses of HPAV which is typically reserved for viperid snakes instead. Following administration of the anti-venom, the patient's coagulation profile improved, and the local soft tissue effects of the venom resolved. He did not manifest any adverse effects and was discharged uneventfully about 72 h after the snakebite. The cross-neutralization potential of HPAV for Shore Pit Viper (Trimeresurus Purpureomaculatus) venom in this case study suggests that there may be a possible common underlying chemical structure and pathophysiology among the venom proteins of various snake species. Given that Trimeresurus-specific antivenom is unavailable in most countries, this cross-neutralization strategy deserves further consideration and evaluation in similar circumstances.

A Comparative Analysis of the Cytotoxic and Vascular Activity Effects of Western Diamondback Rattlesnake (Crotalus atrox) and Eastern Diamondback Rattlesnake (Crotalus adamanteus) Venoms Using a Chick Embryo Model

Crotalus snakebites induce various toxicological effects, encompassing neurological, myotoxic, and cytotoxic symptoms, with potentially fatal outcomes. Investigating venom toxicity is essential for public health, and developing new tools allows for these effects to be studied more comprehensively. The research goals include the elucidation of the physiological consequences of venom exposure and the assessment of toxicity using animal models. Chicken embryos serve as valuable models for assessing venom toxicity through the chick embryotoxicity screening test (CHEST) and the chick chorioallantoic membrane (CAM) assay, particularly useful for evaluating vascular impacts. C. adamanteus venom application resulted in higher embryotoxicity and morphological abnormalities, such as Siamese twins. The CAM assay demonstrated the hemorrhagic effects of venom, varying with venom type and concentration. The irritant potential of both venom types was classified as slight or moderate depending on their concentration. Additionally, acetylcholinesterase (AChE) activity was performed to receive information about organ toxicity. The results show that both venoms induced changes in the whole embryo, heart, and liver weights, but the C. adamanteus venom was identified as more toxic. Specific venom concentrations affected AChE activity in embryonic tissues. These findings underscore the embryotoxic and vasoactive properties of Crotalus venoms, providing valuable insights into their mechanisms of toxicity and potential applications in biomedicine.

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

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

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.

Toxinological profile and histopathological alterations induced by Bothrocophias campbelli venom from Colombia

Snakebite envenomings most frequently reported in Colombia are caused by snakes of the genera Bothrops, Bothriechis, Bothrocophias, and Porthidium. Their venoms induce local and systemic pathophysiological effects, sometimes leading to permanent sequelae such as reduced mobility of the limbs, amputations, besides the risk of death. The genus Bothrocophias includes nine species, among which B. campbelli has a distribution restricted to the department of Nariño in Colombia. In this work we determined the toxinological profile its venom, by performing assays for the lethal, hemorrhagic, edematogenic, and myotoxic activities in mouse models, as well as for in vitro coagulant activity on human plasma. The lethal toxicity of the venom was 142.7 µg venom/mouse (111.4-179.8 µg/mouse; 6.6-10.6 µg/g body weight) by intraperitoneal route. Its hemorrhagic activity (minimum hemorrhagic dose: 12.7 ± 2.3 µg) is generally weaker compared to other South American vipers, but edematogenic (minimum edematogenic dose 1.0 ± 0.3 µg), and myotoxic (minimum myotoxic dose 3.9 ± 2.5 µg) activities are very potent. Histopathological examination of the injected mouse gastrocnemius muscle showed prominent disorganization of the myofibrils, myonecrosis, and an intense inflammatory leukocyte infiltrate. In vitro, the minimal coagulant dose was 12.3 ± 0.5 µg. Overall, this toxinological profile would predict that the clinical picture of envenomings by B. campbelli might be characterized by moderate disturbances in the coagulation cascade, mild local hemorrhage, and, conversely, severe myonecrosis and edema, which could potentially lead to compartment syndrome and gangrene.

Hidden in the sand: Phylogenomics unravel an unexpected evolutionary history for the desert-adapted vipers of the genus Cerastes

The desert vipers of the genus Cerastes are a small clade of medically important venomous snakes within the family Viperidae. According to published morphological and molecular studies, the group is comprised by four species: two morphologically similar and phylogenetically sister taxa, the African horned viper (Cerastes cerastes) and the Arabian horned viper (Cerastes gasperettii); a more distantly related species, the Saharan sand viper (Cerastes vipera), and the enigmatic Böhme's sand viper (Cerastes boehmei), only known from a single specimen in captivity allegedly captured in Central Tunisia. In this study, we sequenced one mitochondrial marker (COI) as well as genome-wide data (ddRAD sequencing) from 28 and 41 samples, respectively, covering the entire distribution range of the genus to explore the population genomics, phylogenomic relationships and introgression patterns within the genus Cerastes. Additionally, and to provide insights into the mode of diversification of the group, we carried out niche overlap analyses considering climatic and habitat variables. Both nuclear phylogenomic reconstructions and population structure analyses have unveiled an unexpected evolutionary history for the genus Cerastes, which sharply contradicts the morphological similarities and previously published mitochondrial approaches. Cerastes cerastes and C. vipera are recovered as sister taxa whilst C. gasperettii is a sister taxon to the clade formed by these two species. We found a relatively high niche overlap (OI > 0.7) in both climatic and habitat variables between C. cerastes and C. vipera, contradicting a potential scenario of sympatric speciation. These results are in line with the introgression found between the northwestern African populations of C. cerastes and C. vipera. Finally, our genomic data confirms the existence of a lineage of C. cerastes in Arabia. All these results highlight the importance of genome-wide data over few genetic markers to study the evolutionary history of species.

The Brazilian Atlantic Bushmaster Lachesis (Linnaeus, 1766) Mitogenome With Insights On Snake Evolution And Divergence (Serpentes: Viperidae: Crotalinae)

This study presents the first complete mitogenome of the Brazilian Atlantic bushmaster Lachesis with insights into snake evolution. The total length was 17,177 bp, consisting of 13 PCGs, 22 tRNAs, two rRNAs and a duplicate control region (CRs). Almost all genes were encoded by the heavy-strand, except for the ND6 gene and eight tRNAs (tRNA-Gln, Ala, Asn, Cys, Tyr, Ser[TGA anticodon], Glu, Pro). Only ATG, ATA, and ATC were starting codons for protein-coding sequences. Stop codons mainly were TAA, AGA, AGG, and TAG; whereas ND1, ND3, and CYTB terminated with incomplete stop codons. Phylogeny retrieved Lachesis within the Crotalinae as the sister group of Agkistrodon; and the Lachesis+Agkistrodon clade as the sister group of (Sistrurus+Crotalus)+Bothrops. The tree supports Crotalinae, Viperinae, and Azemiopinae in the Viperidae family, being sister taxa of Colubridae+(Elapidae+Psammophiidae). The mean genetic distance across 15 snake families and 57 nucleotide sequences was 0.37. The overall mean value of genetic distance across the Crotalinae was 0.23, with Lachesis muta exhibiting the shortest distance of 0.2 with Agkistrodon piscivorus, Protobothrops dabieshanensis and P. flavoviridis and the greatest 0.25 with Gloydius blomhoffii, Trimeresurus albolabris, S. miliarius, and Deinagkistrodon acutus. The complete Atlantic L. muta mitogenome presented herein is only the third annotated mitogenome from more than 430 described Brazilian snake species.

Comparative Analysis of Alpha-1 Orthosteric-Site Binding by a Clade of Central American Pit Vipers (Genera Atropoides, Cerrophidion, Metlapilcoatlus, and Porthidium)

The distribution and relative potency of post-synaptic neurotoxic activity within Crotalinae venoms has been the subject of less investigation in comparison with Elapidae snake venoms. No previous studies have investigated post-synaptic neurotoxic activity within the Atropoides, Metlapilcoatlus, Cerrophidion, and Porthidium clade. Given the specificity of neurotoxins to relevant prey types, we aimed to uncover any activity present within this clade of snakes that may have been overlooked due to lower potency upon humans and thus not appearing as a clinical feature. Using biolayer interferometry, we assessed the relative binding of crude venoms to amphibian, lizard, bird, rodent and human α-1 nAChR orthosteric sites. We report potent alpha-1 orthosteric site binding in venoms from Atropoides picadoi, Metlapilcoatlus occiduus, M. olmec, M. mexicanus, M. nummifer. Lower levels of binding, but still notable, were evident for Cerrophidion godmani, C. tzotzilorum and C. wilsoni venoms. No activity was observed for Porthidium venoms, which is consistent with significant alpha-1 orthosteric site neurotoxicity being a trait that was amplified in the last common ancestor of Atropoides/Cerrophidion/Metlapilcoatlus subsequent to the split by Porthidium. We also observed potent taxon-selective activity, with strong selection for non-mammalian targets (amphibian, lizard, and bird). As these are poorly studied snakes, much of what is known about them is from clinical reports. The lack of affinity towards mammalian targets may explain the knowledge gap in neurotoxic activity within these species, since symptoms would not appear in bite reports. This study reports novel venom activity, which was previously unreported, indicating toxins that bind to post-synaptic receptors may be more widespread in pit vipers than previously considered. While these effects appear to not be clinically significant due to lineage-specific effects, they are of significant evolutionary novelty and of biodiscovery interest. This work sets the stage for future research directions, such as the use of in vitro and in vivo models to determine whether the alpha-1 orthosteric site binding observed within this study confers neurotoxic venom activity.

Speciation rates are unrelated to the formation of population structure in Malagasy gemsnakes

Speciation rates vary substantially across the tree of life. These rates should be linked to the rate at which population structure forms if a continuum between micro and macroevolutionary patterns exists. Previous studies examining the link between speciation rates and the degree of population formation in clades have been shown to be either correlated or uncorrelated depending on the group, but no study has yet examined the relationship between speciation rates and population structure in a young group that is constrained spatially to a single-island system. We examine this correlation in 109 gemsnakes (Pseudoxyrhophiidae) endemic to Madagascar and originating in the early Miocene, which helps control for extinction variation across time and space. We find no relationship between rates of speciation and the formation rates of population structure over space in 33 species of gemsnakes. Rates of speciation show low variation, yet population structure varies widely across species, indicating that speciation rates and population structure are disconnected. We suspect this is largely due to the persistence of some lineages not susceptible to extinction. Importantly, we discuss how delimiting populations versus species may contribute to problems understanding the continuum between shallow and deep evolutionary processes.

Crosstalk of Inflammation and Coagulation in Bothrops Snakebite Envenoming: Endogenous Signaling Pathways and Pathophysiology

Snakebite envenoming represents a major health problem in tropical and subtropical countries. Considering the elevated number of accidents and high morbidity and mortality rates, the World Health Organization reclassified this disease to category A of neglected diseases. In Latin America, Bothrops genus snakes are mainly responsible for snakebites in humans, whose pathophysiology is characterized by local and systemic inflammatory and degradative processes, triggering prothrombotic and hemorrhagic events, which lead to various complications, organ damage, tissue loss, amputations, and death. The activation of the multicellular blood system, hemostatic alterations, and activation of the inflammatory response are all well-documented in Bothrops envenomings. However, the interface between inflammation and coagulation is still a neglected issue in the toxinology field. Thromboinflammatory pathways can play a significant role in some of the major complications of snakebite envenoming, such as stroke, venous thromboembolism, and acute kidney injury. In addition to exacerbating inflammation and cell interactions that trigger vaso-occlusion, ischemia-reperfusion processes, and, eventually, organic damage and necrosis. In this review, we discuss the role of inflammatory pathways in modulating coagulation and inducing platelet and leukocyte activation, as well as the inflammatory production mediators and induction of innate immune responses, among other mechanisms that are altered by Bothrops venoms.

Venom Proteomics of Trimeresurus gracilis, a Taiwan-Endemic Pitviper, and Comparison of Its Venom Proteome and VEGF and CRISP Sequences with Those of the Most Related Species

Trimeresurus gracilis is an endemic alpine pitviper in Taiwan with controversial phylogeny, and its venom proteome remains unknown. In this study, we conducted a proteomic analysis of T. gracilis venom using high-performance liquid chromatography-tandem mass spectrometry and identified 155 toxin proteoforms that belong to 13 viperid venom toxin families. By searching the sequences of trypsin-digested peptides of the separated HPLC fractions against the NCBI database, T. gracilis venom was found to contain 40.3% metalloproteases (SVMPs), 15.3% serine proteases, 6.6% phospholipases A2, 5.0% L-amino acid oxidase, 4.6% Cys-rich secretory proteins (CRISPs), 3.2% disintegrins, 2.9% vascular endothelial growth factors (VEGFs), 1.9% C-type lectin-like proteins, and 20.2% of minor toxins, nontoxins, and unidentified peptides or compounds. Sixteen of these proteoforms matched the toxins whose full amino-acid sequences have been deduced from T. gracilis venom gland cDNA sequences. The hemorrhagic venom of T. gracilis appears to be especially rich in PI-class SVMPs and lacks basic phospholipase A2. We also cloned and sequenced the cDNAs encoding two CRISP and three VEGF variants from T. gracilis venom glands. Sequence alignments and comparison revealed that the PI-SVMP, kallikrein-like proteases, CRISPs, and VEGF-F of T. gracilis and Ovophis okinavensis are structurally most similar, consistent with their close phylogenetic relationship. However, the expression levels of some of their toxins were rather different, possibly due to their distinct ecological and prey conditions.

The Phylogenomic and Biogeographic History of the Gartersnakes, Watersnakes, and Allies (Natricidae: Thamnophiini)

North American Thamnophiini (gartersnakes, watersnakes, brownsnakes, and swampsnakes) are an ecologically and phenotypically diverse temperate clade of snakes representing 61 species across 10 genera. In this study, we estimate phylogenetic trees using ∼3,700 ultraconserved elements (UCEs) for 76 specimens representing 75% of all Thamnophiini species. We infer phylogenies using multispecies coalescent methods and time calibrate them using the fossil record. We also conducted ancestral area estimation to identify how major biogeographic boundaries in North America affect broadscale diversification in the group. While most nodes exhibited strong statistical support, analysis of concordant data across gene trees reveals substantial heterogeneity. Ancestral area estimation demonstrated that the genus Thamnophis was the only taxon in this subfamily to cross the Western Continental Divide, even as other taxa dispersed southward toward the tropics. Additionally, levels of gene tree discordance are overall higher in transition zones between bioregions, including the Rocky Mountains. Therefore, the Western Continental Divide may be a significant transition zone structuring the diversification of Thamnophiini during the Neogene and Pleistocene. Here we show that despite high levels of discordance across gene trees, we were able to infer a highly resolved and well-supported phylogeny for Thamnophiini, which allows us to understand broadscale patterns of diversity and biogeography.

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

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

T. Lalremsanga EO, Bhosale H, Gowande G, Patel H, Idiatullina SS, Poyarkov NA. Systematics of Trimeresurus popeiorum Smith, 1937 with a revised molecular phylogeny of Asian pitvipers of the genus Trimeresurus Lacépède, 1804 sensu lato

The pit viper snake genus Trimeresurus Lacépède, 1804 sensu lato, is a diverse group of nocturnal serpents comprising over 61 species. The genus is morphologically heterogeneous and has been divided into several subgenera. We present an updated phylogeny of Asian pitvipers and propose a revised classification. Additionally, we revise the taxonomy of T. popeiorum Smith, 1937 and propose taxonomic changes with support from molecular and morphological data. We restrict T. popeiorumsensu stricto to northeastern India, Bangladesh, southern China, and northern Myanmar; populations beyond these areas require further assessment. We also synonymize T. yingjiangensisChen et al., 2019 with T. popeiorum based on overlapping morphological characters, molecular data, and distribution. The findings shed new light on the taxonomy of T. popeiorum, warranting the need for assessing the population of T. popeiorum from southeast Asia.

Cytogenetic Analysis of the Bimodal Karyotype of the Common European Adder, Vipera berus (Viperidae)

Vipera berus is the species with the largest range of snakes on Earth and one of the largest among reptiles in general. It is also the only snake species found in the Arctic Circle. Vipera berus is the most involved species of the genus Vipera in the process of interspecific hybridization in nature. The taxonomy of the genus Vipera is based on molecular markers and morphology and requires clarification using SC-karyotyping. This work is a detailed comparative study of the somatic and meiotic karyotypes of V. berus, with special attention to DNA and protein markers associated with synaptonemal complexes. The karyotype of V. berus is a remarkable example of a bimodal karyotype containing both 16 large macrochromosomes and 20 microchromosomes. We traced the stages of the asynchronous assembly of both types of bivalents. The number of crossing-over sites per pachytene nucleus, the localization of the nucleolar organizer, and the unique heterochromatin block on the autosomal bivalent 6-an important marker-were determined. Our results show that the average number of crossing-over sites per pachytene nucleus is 49.5, and the number of MLH1 sites per bivalent 1 reached 11, which is comparable to several species of agamas.

Extreme Procoagulant Potency in Human Plasma of Venoms from the African Viperid Genera Atheris, Cerastes, and Proatheris and the Relative Efficacy of Antivenoms and Synthetic Enzyme-Inhibitors

The African viperid snake genera Atheris, Cerastes, and Proatheris are closely related, similar in size, but occupy extremely divergent ecological niches (arboreal in tropical rainforests, fossorial in deserts, and swamp-dwelling, respectively). Their venoms have not previously been subjected to comparative analyses for their action upon the coagulation of blood, most notably with significant data deficiencies from Atheris and Proatheris. In contrast, the closely related genus Echis is well-documented as capable of producing potent procoagulant effects. In light of this, we set out to compare the coagulotoxic actions of Atheris ceratophora, A. chlorechis, A. desaixi, A. nitschei, A. squamigera, C. cerastes, C. cerastes gasperettii, C. vipera, and Proatheris superciliaris and explore potential pharmacological interventions to reestablish normal blood coagulation. All venoms displayed extremely potent procoagulant effects, over twice as fast as the most potent Echis reported to date. Although Cerastes is used in the immunising mixture of two different regionally available antivenoms (Inoserp-MENA with C. cerastes, C. cerastes gasperettii, C. vipera and Saudi Arabian polyvalent with C. cerastes), none of the other species in this study are included in the immunising mixture of any antivenom. Notably, all the Cerastes species were only neutralised by the Inoserp-MENA antivenom. C. cerastes venom was not neutralised well by the Saudi Arabian antivenom, with the low levels of recognition for any of the Cerastes venoms suggesting a strong regional variation in the venom of this species, as the C. cerastes venom tested was of African (Tunisian) origin versus Saudi locality used in that antivenom's production. The other antivenoms (Micropharm EchiTAbG, ICP EchiTAb-Plus-ICP, Inosan Inoserp Pan-Africa, Premium Serums PANAF Sub-Sahara Africa, South African Vaccine Producers Echis, South African Vaccine Producers Polyvalent) all displayed trivial-to-no ability to neutralise the procoagulant toxicity of any of the Atheris, Cerastes, or Proatheris venoms. Comparative testing of the enzyme inhibitors DMPS, marimastat, and prinomastat, revealed a very potent neutralising capacity of marimastat, with prinomastat showing lower but still significant potency at the same molar concentration, while a 5× molar concentration of DMPS had no apparent effect on procoagulant venom effects normalized by the other inhibitors. These results and methods contribute to the body of knowledge of potential clinical effects and data necessary for evidence-based advancement of clinical management strategies.

Evolutionary allometry and ecological correlates of fang length evolution in vipers

Traits for prey acquisition form the phenotypic interface of predator-prey interactions. In venomous predators, morphological variation in venom delivery apparatus like fangs and stingers may be optimized for dispatching prey. Here, we determine how a single dimension of venom injection systems evolves in response to variation in the size, climatic conditions and dietary ecology of viperid snakes. We measured fang length in more than 1900 museum specimens representing 199 viper species (55% of recognized species). We find both phylogenetic signal and within-clade variation in relative fang length across vipers suggesting both general taxonomic trends and potential adaptive divergence in fang length. We recover positive evolutionary allometry and little static allometry in fang length. Proportionally longer fangs have evolved in larger species, which may facilitate venom injection in more voluminous prey. Finally, we leverage climatic and diet data to assess the global correlates of fang length. We find that models of fang length evolution are improved through the inclusion of both temperature and diet, particularly the extent to which diets are mammal-heavy diets. These findings demonstrate how adaptive variation can emerge among components of complex prey capture systems.

Differential Antivenom and Small-Molecule Inhibition of Novel Coagulotoxic Variations in Atropoides, Cerrophidion, Metlapilcoatlus, and Porthidium American Viperid Snake Venoms

Within Neotropical pit-vipers, the Mexican/Central-American clade consisting of Atropoides, Cerrophidion, Metlapilcoatlus, and Porthidium is a wide-ranging, morphologically and ecologically diverse group of snakes. Despite their prevalence, little is known of the functional aspects of their venoms. This study aimed to fill the knowledge gap regarding coagulotoxic effects and to examine the potential of different therapeutic approaches. As a general trait, the venoms were shown to be anticoagulant but were underpinned by diverse biochemical actions. Pseudo-procoagulant activity (i.e., thrombin-like), characterized by the direct cleavage of fibrinogen to form weak fibrin clots, was evident for Atropoides picadoi, Cerrophidiontzotzilorum, Metlapilcoatlus mexicanus, M. nummifer, M. occiduus, M. olmec, and Porthidium porrasi. In contrast, other venoms cleaved fibrinogen in a destructive (non-clotting) manner, with C. godmani and C. wilsoni being the most potent. In addition to actions on fibrinogen, clotting enzymes were also inhibited. FXa was only weakly inhibited by most species, but Cerrophidion godmani and C. wilsoni were extremely strong in their inhibitory action. Other clotting enzymes were more widely inhibited by diverse species spanning the full taxonomical range, but in each case, there were species that had these traits notably amplified relatively to the others. C. godmani and C. wilsoni were the most potent amongst those that inhibited the formation of the prothrombinase complex and were also amongst the most potent inhibitors of Factor XIa. While most species displayed only low levels of thrombin inhibition, Porthidium dunni potently inhibited this clotting factor. The regional polyvalent antivenom produced by Instituto Picado Clodomiro was tested and was shown to be effective against the diverse anticoagulant pathophysiological effects. In contrast to the anticoagulant activities of the other species, Porthidium volcanicum was uniquely procoagulant through the activation of Factor VII and Factor XII. This viperid species is the first snake outside of the Oxyuranus/Pseudonaja elapid snake clade to be shown to activate FVII and the first snake venom of any kind to activate FXII. Interestingly, while small-molecule metalloprotease inhibitors prinomastat and marimastat demonstrated the ability to prevent the procoagulant toxicity of P. volcanicum, neither ICP antivenom nor inhibitor DMPS showed this effect. The extreme variation among the snakes here studied underscores how venom is a dynamic trait and how this can shape clinical outcomes and influence evolving treatment strategies.

Phylogeography supports lineage divergence for an endemic rattlesnake (Crotalus ravus) of the Neotropical montane forest in the Trans-Mexican Volcanic Belt

The formation of the Trans-Mexican Volcanic Belt (TMVB) and Pleistocene climatic fluctuations have been shown to influence the diversification of lineages and species distributed throughout central Mexico. In some taxa, however, evidence of lineage diversification is not easily recognized, as often is the case in reptiles. Here we present a phylogeographic study on a Mexican endemic rattlesnake species (Crotalus ravus), with the aim of understanding how distinct lineages are distributed across the TMVB. Genetic (mtDNA) and genomic (ddRADseq) data were generated from samples across the species’ range to evaluate phylogeographic structure, estimate phylogenetic relationships and divergence times, and perform environmental niche modeling (ENM). Both datasets recover strong phylogeographic structuring of two distinct lineages on an east-west axis, with an estimated Pleistocene divergence (~1.47 Myr). The ENM suggest that the distribution of the two lineages experienced expansion and reduction events throughout recent evolutionary time. We attribute the diversification of C. ravus lineages to geological events associated with the formation of the TMVB, as well as Quaternary climate changes, both of which have been previously recognized in co-distributed taxa in the TMVB. This work emphasizes the existence of cryptic diversification processes in a morphologically conserved species distributed in a region of complex climatic and orogenic heterogeneity.

A tale about vipers’ tails: phylogeography of black-tailed rattlesnakes

The phylogenetic relationships among black-tailed rattlesnakes remain poorly understood and some authors indicated that the diversity of this group has been underestimated and additional analyses are required to clarify the biogeographic patterns throughout its distribution in Mexico. Therefore, the aim of this study was to elucidate the phylogenetic relationships among black-tailed rattlesnakes across their range, identifying relative divergence times among the main clades and reconstructing the biogeographical history of the group. Three partial mitochondrial genes (ND4, cytb and ATPase6) and one nuclear gene (RAG1) were sequenced to infer the phylogenetic relationships, through the maximum likelihood and Bayesian inference-based methods; demographic history reconstruction was investigated through Bayesian Skyline plot analysis and the ancestral area reconstruction was carried out considering a Bayesian framework. We found strong evidence that the black-tailed rattlesnakes’ group is composed of six clades, which is in agreement with subspecies previously reported. Divergence time estimation indicated that the origin of the C. molossus group could be traced to the middle of the Miocene (~7.71 Mya). Ancestral area reconstruction indicated that early divergence events occurred in Central Mexico, probably related to the geological dynamics of the Trans-Mexican Volcanic Belt. The lineage C. m. oaxacus is the basal member of the C. molossus group. Furthermore, the combination of geological events and changes in Quaternary vegetation may have contributed to the divergence of C. molossus clades. Our results suggest several clades within C. molossus complex could be potentially recognized as separate species.

Maximizing Molecular Data From Low-Quality Fluid-Preserved Specimens in Natural History Collections

Over the past decade, museum genomics studies have focused on obtaining DNA of sufficient quality and quantity for sequencing from fluid-preserved natural history specimens, primarily to be used in systematic studies. While these studies have opened windows to evolutionary and biodiversity knowledge of many species worldwide, published works often focus on the success of these DNA sequencing efforts, which is undoubtedly less common than obtaining minimal or sometimes no DNA or unusable sequence data from specimens in natural history collections. Here, we attempt to obtain and sequence DNA extracts from 115 fresh and 41 degraded samples of homalopsid snakes, as well as from two degraded samples of a poorly known snake, Hydrablabes periops. Hydrablabes has been suggested to belong to at least two different families (Natricidae and Homalopsidae) and with no fresh tissues known to be available, intractable museum specimens currently provide the only opportunity to determine this snake’s taxonomic affinity. Although our aim was to generate a target-capture dataset for these samples, to be included in a broader phylogenetic study, results were less than ideal due to large amounts of missing data, especially using the same downstream methods as with standard, high-quality samples. However, rather than discount results entirely, we used mapping methods with references and pseudoreferences, along with phylogenetic analyses, to maximize any usable molecular data from our sequencing efforts, identify the taxonomic affinity of H. periops, and compare sequencing success between fresh and degraded tissue samples. This resulted in largely complete mitochondrial genomes for five specimens and hundreds to thousands of nuclear loci (ultra-conserved loci, anchored-hybrid enrichment loci, and a variety of loci frequently used in squamate phylogenetic studies) from fluid-preserved snakes, including a specimen of H. periops from the Field Museum of Natural History collection. We combined our H. periops data with previously published genomic and Sanger-sequenced datasets to confirm the familial designation of this taxon, reject previous taxonomic hypotheses, and make biogeographic inferences for Hydrablabes. A second H. periops specimen, despite being seemingly similar for initial raw sequencing results and after being put through the same protocols, resulted in little usable molecular data. We discuss the successes and failures of using different pipelines and methods to maximize the products from these data and provide expectations for others who are looking to use DNA sequencing efforts on specimens that likely have degraded DNA.

Life Science Identifier (Hydrablabes periops)

urn:lsid:zoobank.org:pub:F2AA44 E2-D2EF-4747-972A-652C34C2C09D.

Contextual Constraints: Dynamic Evolution of Snake Venom Phospholipase A2

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

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

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

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

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

Clinical and Evolutionary Implications of Dynamic Coagulotoxicity Divergences in Bothrops (Lancehead Pit Viper) Venoms

Despite coagulotoxicity being a primary weapon for prey capture by Bothrops species (lancehead pit vipers) and coagulopathy being a major lethal clinical effect, a genus-wide comparison has not been undertaken. To fill this knowledge gap, we used thromboelastography to compare 37 venoms, from across the full range of geography, taxonomy, and ecology, for their action upon whole plasma and isolated fibrinogen. Potent procoagulant toxicity was shown to be the main venom effect of most of the species tested. However, the most basal species (B. pictus) was strongly anticoagulant; this is consistent with procoagulant toxicity being a novel trait that evolved within Bothrops subsequent to their split from anticoagulant American pit vipers. Intriguingly, two of the arboreal species studied (B. bilineatus and B. taeniatus) lacked procoagulant venom, suggesting differential evolutionary selection pressures. Notably, some terrestrial species have secondarily lost the procoagulant venom trait: the Mogi Mirim, Brazil locality of B. alternatus; San Andres, Mexico locality of B. asper; B. diporus; and the São Roque of B. jararaca. Direct action on fibrinogen was extremely variable; this is consistent with previous hypotheses regarding it being evolutionary decoupled due to procoagulant toxicity being the primary prey-capture weapon. However, human patients live long enough for fibrinogen depletion to be clinically significant. The extreme variability may be reflective of antivenom variability, with these results thereby providing a foundation for such future work of clinical relevance. Similarly, the venom diversification trends relative to ecological niche will also be useful for integration with natural history data, to reconstruct the evolutionary pressures shaping the venoms of these fascinating snakes.

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.

Efficacy and Limitations of Chemically Diverse Small-Molecule Enzyme-Inhibitors against the Synergistic Coagulotoxic Activities of Bitis Viper Venoms

Snakebite remains a significant public health burden globally, disproportionately affecting low-income and impoverished regions of the world. Recently, researchers have begun to focus on the use of small-molecule inhibitors as potential candidates for the neutralisation of key snake venom toxins and as potential field therapies. Bitis vipers represent some of the most medically important as well as frequently encountered snake species in Africa, with a number of species possessing anticoagulant phospholipase A2 (PLA2) toxins that prevent the prothrombinase complex from inducing clot formation. Additionally, species within the genus are known to exert pseudo-procoagulant activity, whereby kallikrein enzymatic toxins cleave fibrinogen to form a weak fibrin clot that rapidly degrades, thereby depleting fibrinogen levels and contributing to the net anticoagulant state. Utilising well-validated coagulation assays measuring time until clot formation, this study addresses the in vitro efficacy of three small molecule enzyme inhibitors (marimastat, prinomastat and varespladib) in neutralising these aforementioned activities. The PLA2 inhibitor varespladib showed the greatest efficacy for the neutralisation of PLA2-driven anticoagulant venom activity, with the metalloproteinase inhibitors prinomastat and marimastat both showing low and highly variable degrees of cross-neutralisation with PLA2 anticoagulant toxicity. However, none of the inhibitors showed efficacy in neutralising the pseudo-procoagulant venom activity exerted by the venom of B. caudalis. Our results highlight the complex nature of snake venoms, for which single-compound treatments will not be universally effective, but combinations might prove highly effective. Despite the limitations of these inhibitors with regards to in vitro kallikrein enzyme pseudo-procoagulant venom activity, our results further support the growing body of literature indicating the potential use of small molecule inhibitors to enhance first-aid treatment of snakebite envenoming, particularly in cases where hospital and thus antivenom treatment is either unavailable or far away.

Isolated by dry lands: integrative analyses unveil the existence of a new species and a previously unknown evolutionary lineage of Brazilian Lanceheads (Serpentes: Viperidae: Bothrops) from a Caatinga moist-forest enclave

Here we describe a new species of Lancehead (Bothrops jabrensis sp. nov.) based on three individuals sampled from a previously unknown population from Pico do Jabre, an isolated and small Caatinga moist-forest enclave (CMFE) located in northeastern Brazil. Although this new species has an external morphology resembling those found in representatives of the Bothrops jararaca (Wied-Neuwied, 1824) species group, B. jabrensis can be diagnosed by a combination of meristic and color characters. Molecular phylogenetic analysis indicates the new species represents a unique and highly divergent lineage within Bothrops revealing the existence of a previously unknown phylogenetic lineage that has been evolving as an independent unit for more than 8 million years. Additionally, the estimated divergence time of this lineage conflicts with some proposed scenarios of historical processes associated with the evolution of CMFEs. Finally, the uniqueness of this species indicates its relevance for the maintenance of the phylogenetic diversity of Lanceheads in South America. Like other CMFEs, Pico do Jabre is consistently threatened by poaching, illegal fires, deforestation for agricultural purposes, and illegal logging. The restricted distribution of B. jabrensis, in a small and disturbed CMFE, strongly suggests that this species is critically endangered and is likely approaching extinction as a natural population.

Dynamic genetic differentiation drives the widespread structural and functional convergent evolution of snake venom proteinaceous toxins

Background

The explosive radiation and diversification of the advanced snakes (superfamily Colubroidea) was associated with changes in all aspects of the shared venom system. Morphological changes included the partitioning of the mixed ancestral glands into two discrete glands devoted for production of venom or mucous respectively, as well as changes in the location, size and structural elements of the venom-delivering teeth. Evidence also exists for homology among venom gland toxins expressed across the advanced snakes. However, despite the evolutionary novelty of snake venoms, in-depth toxin molecular evolutionary history reconstructions have been mostly limited to those types present in only two front-fanged snake families, Elapidae and Viperidae. To have a broader understanding of toxins shared among extant snakes, here we first sequenced the transcriptomes of eight taxonomically diverse rear-fanged species and four key viperid species and analysed major toxin types shared across the advanced snakes.

Results

Transcriptomes were constructed for the following families and species: Colubridae - Helicops leopardinus, Heterodon nasicus, Rhabdophis subminiatus; Homalopsidae – Homalopsis buccata; Lamprophiidae - Malpolon monspessulanus, Psammophis schokari, Psammophis subtaeniatus, Rhamphiophis oxyrhynchus; and Viperidae – Bitis atropos, Pseudocerastes urarachnoides, Tropidolaeumus subannulatus, Vipera transcaucasiana. These sequences were combined with those from available databases of other species in order to facilitate a robust reconstruction of the molecular evolutionary history of the key toxin classes present in the venom of the last common ancestor of the advanced snakes, and thus present across the full diversity of colubroid snake venoms. In addition to differential rates of evolution in toxin classes between the snake lineages, these analyses revealed multiple instances of previously unknown instances of structural and functional convergences. Structural convergences included: the evolution of new cysteines to form heteromeric complexes, such as within kunitz peptides (the beta-bungarotoxin trait evolving on at least two occasions) and within SVMP enzymes (the P-IIId trait evolving on at least three occasions); and the C-terminal tail evolving on two separate occasions within the C-type natriuretic peptides, to create structural and functional analogues of the ANP/BNP tailed condition. Also shown was that the de novo evolution of new post-translationally liberated toxin families within the natriuretic peptide gene propeptide region occurred on at least five occasions, with novel functions ranging from induction of hypotension to post-synaptic neurotoxicity. Functional convergences included the following: multiple occasions of SVMP neofunctionalised in procoagulant venoms into activators of the clotting factors prothrombin and Factor X; multiple instances in procoagulant venoms where kunitz peptides were neofunctionalised into inhibitors of the clot destroying enzyme plasmin, thereby prolonging the half-life of the clots formed by the clotting activating enzymatic toxins; and multiple occasions of kunitz peptides neofunctionalised into neurotoxins acting on presynaptic targets, including twice just within Bungarus venoms.

Conclusions

We found novel convergences in both structural and functional evolution of snake toxins. These results provide a detailed roadmap for future work to elucidate predator–prey evolutionary arms races, ascertain differential clinical pathologies, as well as documenting rich biodiscovery resources for lead compounds in the drug design and discovery pipeline.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01208-9.

Morphology and seasonality of the sexual segment of the kidney in genus Bothrops snakes

In Squamata, the sexual segment of the male kidney (SSK) produces secretory granules that are transmitted to the female tract during copulation. Here, we investigate the morphology of SSK in various species of the pitviper-genus Bothrops, a medically relevant clade of Neotropical snakes, using histology, histochemistry, as well as scanning and transmission electron microscopy. We also evaluated the relationship of SSK hypertrophy with spermatogenesis and storage of sperm. The SSK in Bothrops consists of tall columnar epithelial cells filled with electron-dense secretory granules. These granules are apparently released into the lumen by an apocrine secretory process. In general, the ultrastructural morphology of SSK in Bothrops is similar to those described for other pitvipers in the genus Agkistrodon. In most Bothrops males, the SSK produces carbohydrates and proteins, and occurs in all lobes of the kidney, occupying most of the kidney mass. In contrast, intersexual females (with male copulatory organs) of Bothrops insularis had a smaller SSK diameter than males and did not show hypertrophy. The seasonal development of SSK in Bothrops species seems to be driven by environmental factors rather than phylogeny, since sister species exhibited periods of hypertrophy at different months. In most Bothrops males, hypertrophy occurs at various seasons of the year (spring, summer, autumn), even in winter, and seasonal changes in SSK do not necessarily accompany spermatogenesis.

Sequence determination and bioinformatic comparison of ten venom serine proteases of Trimeresurus gracilis, a Taiwanese endemic pitviper with controversial taxonomy

Trimeresurus gracilis (Tgc) is endemic to Taiwan and shown to be closely related with Ovophis okinavensis by previous phylogenetic analyses, but their taxonomic status remain controversial. Here, we cloned and sequenced ten of its venom serine-proteases (designated as Tgc-vSPs). All the Tgc-vSPs conserve the catalytic triads, six appear to be kallikrein-like (KNs) and four are plasminogen-activator homologs (PAHs and PAs). They are studied under four structural categories: (1) highly similar Tgc-KN1, Tgc-KN2 and Tgc-KN3, with four predicted N-glycosylation sites; (2) Tgc-KN4, with a single N -glycosylation site; (3) Tgc-KN5 and Tgc-KN6, with two distinct N-glycosylation sites; (4) Tgc-PAH1/PAH2, TgcPA3, and Tgc-PA4, with two conserved N-glycosylation sites. Additionally, Tgc-KN1, Tgc-KN4 and Tgc-PAH1 were purified by reversed-phase HPLC and identified by peptide-mass-fingerprinting. Results of BLAST and sequence alignments reveal that Tgc-KN1∼3 and Tgc-KN6 are most like the vSPs of rattlesnakes, while the sequences of Tgc-KN4, KN5 and Tgc-PAH1/PAH2 match closely to the partial sequences of three O. okinavensis vSPs. Thus, our results reveal non-overlapping similarities of Tgc-vSPs to the O. okinavensis vSPs and vSPs of the New World pitvipers. In addition, molecular phylogenetic analyses of the plasminogen-activator like vSPs reveal separate evolution of two clusters of the enzymes with distinct functions.

Proteomics and immunocharacterization of Asian mountain pit viper (Ovophis monticola) venom

The venomic profile of Asian mountain pit viper Ovophis monticola is clarified in the present study. Using mass spectrometry-based proteomics, 247 different proteins were identified in crude venom of O. monticola found in Thailand. The most abundant proteins were snake venom metalloproteases (SVMP) (36.8%), snake venom serine proteases (SVSP) (31.1%), and phospholipases A2 (PLA2) (12.1%). Less abundant proteins included L-amino acid oxidase (LAAO) (5.7%), venom nerve growth factor (3.6%), nucleic acid degrading enzymes (3.2%), C-type lectins (CTL) (1.6%), cysteine-rich secretory proteins (CRISP) (1.2%) and disintegrin (1.2%). The immunoreactivity of this viper's venom to a monovalent antivenom against green pit viper Trimeresurus albolabris, or to a polyvalent antivenom against hemotoxic venom was investigated by indirect ELISA and two-dimensional (2D) immunoblotting. Polyvalent antivenom showed substantially greater reactivity levels than monovalent antivenom. A titer for the monovalent antivenom was over 1:1.28x107 dilution while that of polyvalent antivenom was 1:5.12x107. Of a total of 89 spots comprising 173 proteins, 40 spots of predominantly SVMP, SVSP and PLA2 were specific antigens for antivenoms. The 49 unrecognized spots containing 72 proteins were characterized as non-reactive proteins, and included certain types of CTLs and CRISPs. These neglected venom constituents could limit the effectiveness of antivenom-based therapy currently available for victims of pit viper envenomation.

Monkeying around with venom: an increased resistance to α-neurotoxins supports an evolutionary arms race between Afro-Asian primates and sympatric cobras

Background

Snakes and primates have a multi-layered coevolutionary history as predators, prey, and competitors with each other. Previous work has explored the Snake Detection Theory (SDT), which focuses on the role of snakes as predators of primates and argues that snakes have exerted a selection pressure for the origin of primates’ visual systems, a trait that sets primates apart from other mammals. However, primates also attack and kill snakes and so snakes must simultaneously avoid primates. This factor has been recently highlighted in regard to the movement of hominins into new geographic ranges potentially exerting a selection pressure leading to the evolution of spitting in cobras on three independent occasions.

Results

Here, we provide further evidence of coevolution between primates and snakes, whereby through frequent encounters and reciprocal antagonism with large, diurnally active neurotoxic elapid snakes, Afro-Asian primates have evolved an increased resistance to α-neurotoxins, which are toxins that target the nicotinic acetylcholine receptors. In contrast, such resistance is not found in Lemuriformes in Madagascar, where venomous snakes are absent, or in Platyrrhini in the Americas, where encounters with neurotoxic elapids are unlikely since they are relatively small, fossorial, and nocturnal. Within the Afro-Asian primates, the increased resistance toward the neurotoxins was significantly amplified in the last common ancestor of chimpanzees, gorillas, and humans (clade Homininae). Comparative testing of venoms from Afro-Asian and American elapid snakes revealed an increase in α-neurotoxin resistance across Afro-Asian primates, which was likely selected against cobra venoms. Through structure-activity studies using native and mutant mimotopes of the α-1 nAChR receptor orthosteric site (loop C), we identified the specific amino acids responsible for conferring this increased level of resistance in hominine primates to the α-neurotoxins in cobra venom.

Conclusion

We have discovered a pattern of primate susceptibility toward α-neurotoxins that supports the theory of a reciprocal coevolutionary arms-race between venomous snakes and primates.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01195-x.

A Strategy for Efficient Preparation of Genus-Specific Diagnostic Antibodies for Snakebites

As said by former United Nations Secretary-General Kofi Annan, "Snakebite is the most important tropical disease you've never heard of." Listed as a priority neglected tropical disease by the World Health Organization, snakebite envenoming (SBE) kills in excess of 125,000 people per year. However, due to the complexity and overlap of snake venom compositions, few reliable venom diagnostic methods for genus-/species-specific identification, which is crucial for successful SBE therapy, are available. Here, we develop a strategy to select and prepare genus-specific snake venom antibodies, which allows rapid and efficient clinical diagnosis of snakebite. Multi-omics approaches are used to choose candidate antigens from snake venoms and identify genus-specific antigenic epitope peptide fragments (GSAEPs) with ideal immunogenicity, specificity, and spatial accessibility. Double-antibody sandwich ELISA kit was established by matching a polyclonal antibody against a natural antigen and a monoclonal antibody that was prepared by natural protein as antigen and can specifically target the GSAEPs. The kit shows the ability to accurately identify venoms from similar genera of Trimeresurus and Protobothrops with a detection limit of 6.25 ng/ml on the snake venoms and a little cross-reaction, thus proving high feasibility and applicability.

Venomics of the poorly studied hognosed pitvipers Porthidium arcosae and Porthidium volcanicum

We report the first proteomics analyses of the venoms of two poorly studied snakes, the Manabi hognosed pitviper Porthidium arcosae endemic to the western coastal province of Manabí (Ecuador), and the Costa Rican hognosed pitviper P. volcanicum with distribution restricted to South Pacific Costa Rica and western Panamá. These venom proteomes share a conserved compositional pattern reported in four other congeneric species within the clade of South American Porthidium species, P. nasutum, P. lansbergii, P. ophryomegas, and P. porrasi. The paraspecific immunorecognition profile of antivenoms produced in Costa Rica (ICP polyvalent), Perú (Instituto Nacional de Salud) and Brazil (soro antibotrópico pentavalente, SAB, from Instituto Butantan) against the venom of P. arcosae was investigated through a third-generation antivenomics approach. The maximal venom-binding capacities of the investigated antivenoms were 97.1 mg, 21.8 mg, and 25.7 mg of P. arcosae venom proteins per gram of SAB, ICP, and INS-PERU antibody molecules, respectively, which translate into 28.4 mg, 13.1 mg, and 15.2 mg of total venom proteins bound per vial of SAB, ICP, and INS-PERU AV. The antivenomics results suggest that 21.8%, 7.8% and 6.1% of the SAB, ICP, and INS-PERU antibody molecules recognized P. arcosae venom toxins. The SAB antivenom neutralized P. arcosae venom's lethality in mice with an ED₅₀ of 31.3 mgV/g SAB AV. This preclinical neutralization paraspecificity points to Brazilian SAB as a promising candidate for the treatment of envenomings by Ecuadorian P. arcosae. BIOLOGICAL SIGNIFICANCE: Assessing the preclinical efficacy profile of antivenoms against homologous and heterologous medically relevant snake venoms represents an important goal towards defining the biogeographic range of their clinical utility. This is particularly relevant in regions, such as Mesoamerica, where a small number of pharmaceutical companies produce antivenoms against the venoms of a small number of species of maximum medical relevance among the local rich herpetofauna, leaving a wide range of snakes of secondary medical relevance, but also causing life-threatening human envenomings without nominal clinical coverage. This work is part of a larger project aiming at mapping the immunological characteristics of antivenoms generated in Latin American countries towards venoms of such poorly studied snakes of the local and neighboring countries' herpetofauna. Here we report the proteomics characterization of the Manabi hognosed pitviper Porthidium arcosae endemic to the western coastal province of Manabí (Ecuador), and the Costa Rican hognosed pitviper P. volcanicum with distribution restricted to southwestern Costa Rica, the antivenomics assessment of three bothropoid commercial antivenoms produced in Costa Rica, Perú, and Brazil against the venom components of P. arcosae, and the in vivo capacity of the Brazilian soro antibotrópico pentavalente (SAB) from Instituto Butantan to neutralize the murine lethality of P. arcosae venom. The preclinical paraspecific ED₅₀ of 31.3 mg of P. arcosae venom per gram of antivenom points to Brazilian SAB as a promising candidate for the treatment of envenomings by the Manabi hognosed pitviper P. arcosae.