Abstract: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.

Venom Gland Transcriptomes of Two Elapid Snakes (bungarus Multicinctus and Naja Atra) and Evolution of Toxin Genes

Molecular Diversification of Peptide Toxins from the Tarantula Haplopelma Hainanum (ornithoctonus Hainana) Venom Based on Transcriptomic, Peptidomic, and Genomic Analyses.

Genomic, transcriptomic, and epigenomic analysis of a medicinal snake, Bungarus multicinctus, to provides insights into the origin of Elapidae neurotoxins

Venom-Related Transcripts from Bothrops jararaca Tissues Provide Novel Molecular Insights into the Production and Evolution of Snake Venom

The structural and functional divergence of a neglected three-finger toxin subfamily in lethal elapids

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

Transcriptome Analysis of Deinagkistrodon Acutus Venomous Gland Focusing on Cellular Structure and Functional Aspects Using Expressed Sequence Tags

Distinct regulatory networks control toxin gene expression in elapid and viperid snakes

Venom-gland transcriptomic, venomic, and antivenomic profiles of the spine-bellied sea snake (Hydrophis curtus) from the South China Sea

Hierarchical chromatin features reveal the toxin production in Bungarus multicinctus

A Catalog for Transcripts in the Venom Gland of the Agkistrodon Acutus: Identification of the Toxins Potentially Involved in Coagulopathy

Proteomic characteristics of six snake venoms from the Viperidae and Elapidae families in China and their relation to local tissue necrosis

[Cloning and Characterization of Cdnas of Cardiotoxin-Like Protein in Bungarus Multicinctus Venom Gland].

Proteomic characterization of two snake venoms: Naja naja atra and Agkistrodon halys.

Characterization of Venoms of Deinagkistrodon Acutus and Bungarus Multicinctus Using Proteomics and Peptidomics

Two Non-Coding Transcripts of Toxins from the Scorpion But-hus Martensii Karsch: Evidence for the Existence of NMD in Scorpion Venom Gland Cells

Cdna Cloning and Sequence Analysis of Six Neurotoxin‐like Proteins from Chinese Continental Banded Krait

Evolution of Three-Finger Toxin Genes in Neotropical Colubrine Snakes (Colubridae)

Genomic Sequence Analysis and Organization of BmKαTx11 and BmKαTx15 from Buthus Martensii Karsch: Molecular Evolution of Α-Toxin Genes

Proteomic characterization and comparison of Malaysian Bungarus candidus and Bungarus fasciatus venoms

Exploring the Venom Gland Transcriptome of Bothrops asper and Bothrops jararaca: De Novo Assembly and Analysis of Novel Toxic Proteins