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.

Pseudocerastes urarachnoides is a fascinating viper and as yet has been reported only in western Iran. An elaborated arachnid-like caudal structure is a unique feature of this viper, hence gives it the common name “Iranian spider-tailed viper”. During tail wagging, the structure is reminiscent of a moving spider. Tail movements are used for two different purposes in snakes: defense via tail vibration and hunting via both caudal luring and caudal distraction. Caudal luring in snakes is the wriggling or wagging of the posterior part of tail, in the presence of a potential prey, with conspicuous color pattern while the rest of body is cryptically colored. Previous studies have speculated on the role of caudal structure of P. urarachnoides in hunting. Our 2.5-year study has revealed that development of the structure of the caudal lure is commenced after birth and is linearly correlated to snout-vent length. The caudal lure attracts some species of birds. Caudal luring behavior is carried out both in the presence and absence of birds. The findings are reported for the first time and confirmed by direct observation of undisturbed individuals in the field.

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.