Cloning, expression and characterization of two C-type lectins from the venom gland of Bungarus multicinctus

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.

Protein complexes in snake venom

Snake venom contains mixture of bioactive proteins and polypeptides. Most of these proteins and polypeptides exist as monomers, but some of them form complexes in the venom. These complexes exhibit much higher levels of pharmacological activity compared to individual components and play an important role in pathophysiological effects during envenomation. They are formed through covalent and/or non-covalent interactions. The subunits of the complexes are either identical (homodimers) or dissimilar (heterodimers; in some cases subunits belong to different families of proteins). The formation of complexes, at times, eliminates the non-specific binding and enhances the binding to the target molecule. On several occasions, it also leads to recognition of new targets as protein-protein interaction in complexes exposes the critical amino acid residues buried in the monomers. Here, we describe the structure and function of various protein complexes of snake venoms and their role in snake venom toxicity.

Unusual accelerated rate of deletions and insertions in toxin genes in the venom glands of the pygmy copperhead (Austrelaps labialis) from Kangaroo island

Background

Toxin profiling helps in cataloguing the toxin present in the venom as well as in searching for novel toxins. The former helps in understanding potential pharmacological profile of the venom and evolution of toxins, while the latter contributes to understanding of novel mechanisms of toxicity and provide new research tools or prototypes of therapeutic agents.

Results

The pygmy copperhead (Austrelaps labialis) is one of the less studied species. In this present study, an attempt has been made to describe the toxin profile of A. labialis from Kangaroo Island using the cDNA library of its venom glands. We sequenced 658 clones which represent the common families of toxin genes present in snake venom. They include (a) putative long-chain and short-chain neurotoxins, (b) phospholipase A₂, (c) Kunitz-type protease inhibitor, (d) CRISPs, (e) C-type lectins and (f) Metalloproteases. In addition, we have also identified a novel protein with two Kunitz-type domains in tandem similar to bikunin.

Conclusion

Interestingly, the cDNA library reveals that most of the toxin families (17 out of 43 toxin genes; ~40%) have truncated transcripts due to insertion or deletion of nucleotides. These truncated products might not be functionally active proteins. However, cellular trancripts from the same venom glands are not affected. This unusual higher rate of deletion and insertion of nucleotide in toxin genes may be responsible for the lower toxicity of A. labialis venom of Kangroo Island and have significant effect on evolution of toxin genes.