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.

Pharmacological Characterisation of Pseudocerastes and Eristicophis Viper Venoms Reveal Anticancer (Melanoma) Properties and a Potentially Novel Mode of Fibrinogenolysis.

Venoms are a rich source of potential lead compounds for drug discovery, and descriptive studies of venom form the first phase of the biodiscovery process. In this study, we investigated the pharmacological potential of crude Pseudocerastes and Eristicophis snake venoms in haematological disorders and cancer treatment. We assessed their antithrombotic potential using fibrinogen thromboelastography, fibrinogen gels with and without protease inhibitors, and colourimetric fibrinolysis assays. These assays indicated that the anticoagulant properties of the venoms are likely induced by the hydrolysis of phospholipids and by selective fibrinogenolysis. Furthermore, while most fibrinogenolysis occurred by the direct activity of snake venom metalloproteases and serine proteases, modest evidence indicated that fibrinogenolytic activity may also be mediated by selective venom phospholipases and an inhibitory venom-derived serine protease. We also found that the Pseudocerastes venoms significantly reduced the viability of human melanoma (MM96L) cells by more than 80%, while it had almost no effect on the healthy neonatal foreskin fibroblasts (NFF) as determined by viability assays. The bioactive properties of these venoms suggest that they contain a number of toxins suitable for downstream pharmacological development as candidates for antithrombotic or anticancer agents.

Extensive Variation in the Activities of Pseudocerastes and Eristicophis Viper Venoms Suggests Divergent Envenoming Strategies Are Used for Prey Capture.

Snakes of the genera Pseudocerastes and Eristicophis (Viperidae: Viperinae) are known as the desert vipers due to their association with the arid environments of the Middle East. These species have received limited research attention and little is known about their venom or ecology. In this study, a comprehensive analysis of desert viper venoms was conducted by visualising the venom proteomes via gel electrophoresis and assessing the crude venoms for their cytotoxic, haemotoxic, and neurotoxic properties. Plasmas sourced from human, toad, and chicken were used as models to assess possible prey-linked venom activity. The venoms demonstrated substantial divergence in composition and bioactivity across all experiments. Pseudocerastes urarachnoides venom activated human coagulation factors X and prothrombin and demonstrated potent procoagulant activity in human, toad, and chicken plasmas, in stark contrast to the potent neurotoxic venom of P. fieldi. The venom of E. macmahonii also induced coagulation, though this did not appear to be via the activation of factor X or prothrombin. The coagulant properties of P. fieldi and P. persicus venoms varied among plasmas, demonstrating strong anticoagulant activity in the amphibian and human plasmas but no significant effect in that of bird. This is conjectured to reflect prey-specific toxin activity, though further ecological studies are required to confirm any dietary associations. This study reinforces the notion that phylogenetic relatedness of snakes cannot readily predict venom protein composition or function. The significant venom variation between these species raises serious concerns regarding antivenom paraspecificity. Future assessment of antivenom is crucial.

A new species of the genus Eirenis Jan, 1863 (Squamata: Colubridae) from Kerman Province in South-central Iran.

A new species of Dwarf Snakes, Eirenis rafsanjanicus sp. n., is described from south-central Iran. This species is well distinguished from other congeners in both molecular and morphological characters, and its description increases the number of described Eirenis species to 26. Eirenis rafsanjanicus sp. n. is genetically closest (6% in Cytb) to the recently described E. yassujicus, which inhabits high habitats in southern Iran. Morphologically, E. rafsanjanicus sp. n. is more similar to E. modestus than to other congeners. A revised key to the genus is presented.

Extreme venom variation in Middle Eastern vipers: a proteomics comparison of Eristicophis macmahonii, Pseudocerastes fieldi and Pseudocerastes persicus.

Venoms of the viperid sister genera Eristicophis and Pseudocerastes are poorly studied despite their anecdotal reputation for producing severe or even lethal envenomations. This is due in part to the remote and politically unstable regions that they occupy. All species contained are sit and wait ambush feeders. Thus, this study examined their venoms through proteomics techniques in order to establish if this feeding ecology, and putatively low levels of gene flow, have resulted in significant variations in venom profile. The techniques indeed revealed extreme venom variation. This has immediate implications as only one antivenom is made (using the venom of Pseudocerastes persicus) yet the proteomic variation suggests that it would be of only limited use for the other species, even the sister species Pseudocerastes fieldi. The high degree of variation however also points toward these species being rich resources for novel compounds which may have use as lead molecules in drug design and development. BIOLOGICAL SIGNIFICANCE: These results show extreme venom variation between these closely related snakes. These results have direct implications for the treatment of the envenomed patient.

Molecular phylogeny and historical biogeography of the Anatolian lizard Apathya (Squamata, Lacertidae).

Apathya is a lacertid genus occurring mainly in south-east Turkey and its adjacent regions (part of Iran and Iraq). So far two morphological species have been attributed to the genus; A. cappadocica (with five subspecies, A. c.cappadocica, A. c.muhtari, A. c.schmidtlerorum, A. c. urmiana and A. c.wolteri) and A.yassujica. The first species occupies most of the genus' distribution range, while A. yassujica is endemic of the Zagros Mountains. Here, we explored Apathya's taxonomy and investigated the evolutionary history of the species by employing phylogenetic and phylogeographic approaches and using both mitochondrial (mtDNA) and nuclear markers. The phylogenetic relationships and the genetic distances retrieved, revealed that Apathya is a highly variable genus, which parallels its high morphological variation. Such levels of morphological and genetic differentiation often exceed those between species of other Lacertini genera that are already treated as full species, suggesting the necessity for a taxonomic revision of Apathya. The phylogeographical scenario emerging from the genetic data suggests that the present distribution of the genus was determined by a combination of dispersal and vicariance events between Anatolia and Southwest Asia dating back to the Miocene and continuing up to the Pleistocene. Key geological events for the understanding of the phylogeography of the genus are the movement of the Arabian plate that led to the configuration of Middle East (orogenesis of the mountain ranges of Turkey and Iran) and the formation of Anatolian Diagonal.