Molecular diversity and accelerated evolution of C-type lectin-like proteins from snake venom

Snake venoms are integrated systems, but abundant venom proteins evolve more rapidly

Background

While many studies have shown that extracellular proteins evolve rapidly, how selection acts on them remains poorly understood. We used snake venoms to understand the interaction between ecology, expression level, and evolutionary rate in secreted protein systems. Venomous snakes employ well-integrated systems of proteins and organic constituents to immobilize prey. Venoms are generally optimized to subdue preferred prey more effectively than non-prey, and many venom protein families manifest positive selection and rapid gene family diversification. Although previous studies have illuminated how individual venom protein families evolve, how selection acts on venoms as integrated systems, is unknown.

Results

Using next-generation transcriptome sequencing and mass spectrometry, we examined microevolution in two pitvipers, allopatrically separated for at least 1.6 million years, and their hybrids. Transcriptomes of parental species had generally similar compositions in regard to protein families, but for a given protein family, the homologs present and concentrations thereof sometimes differed dramatically. For instance, a phospholipase A₂ transcript comprising 73.4 % of the Protobothrops elegans transcriptome, was barely present in the P. flavoviridis transcriptome (<0.05 %). Hybrids produced most proteins found in both parental venoms. Protein evolutionary rates were positively correlated with transcriptomic and proteomic abundances, and the most abundant proteins showed positive selection. This pattern holds with the addition of four other published crotaline transcriptomes, from two more genera, and also for the recently published king cobra genome, suggesting that rapid evolution of abundant proteins may be generally true for snake venoms. Looking more broadly at Protobothrops, we show that rapid evolution of the most abundant components is due to positive selection, suggesting an interplay between abundance and adaptation.

Conclusions

Given log-scale differences in toxin abundance, which are likely correlated with biosynthetic costs, we hypothesize that as a result of natural selection, snakes optimize return on energetic investment by producing more of venom proteins that increase their fitness. Natural selection then acts on the additive genetic variance of these components, in proportion to their contributions to overall fitness. Adaptive evolution of venoms may occur most rapidly through changes in expression levels that alter fitness contributions, and thus the strength of selection acting on specific secretome components.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1832-6) contains supplementary material, which is available to authorized users.

Testing the Toxicofera: comparative transcriptomics casts doubt on the single, early evolution of the reptile venom system

The identification of apparently conserved gene complements in the venom and salivary glands of a diverse set of reptiles led to the development of the Toxicofera hypothesis - the single, early evolution of the venom system in reptiles. However, this hypothesis is based largely on relatively small scale EST-based studies of only venom or salivary glands and toxic effects have been assigned to only some putative Toxicoferan toxins in some species. We set out to examine the distribution of these proposed venom toxin transcripts in order to investigate to what extent conservation of gene complements may reflect a bias in previous sampling efforts. Our quantitative transcriptomic analyses of venom and salivary glands and other body tissues in five species of reptile, together with the use of available RNA-Seq datasets for additional species, shows that the majority of genes used to support the establishment and expansion of the Toxicofera are in fact expressed in multiple body tissues and most likely represent general maintenance or "housekeeping" genes. The apparent conservation of gene complements across the Toxicofera therefore reflects an artefact of incomplete tissue sampling. We therefore conclude that venom has evolved multiple times in reptiles.

Proteomic investigation of Sri Lankan hump-nosed pit viper (Hypnale hypnale) venom

The hump-nosed pit viper, Hypanle hypnale, contributes to snakebite mortality and morbidity in Sri Lanka. Studies showed that the venom is hemotoxic and nephrotoxic, with some biochemical and antigenic properties similar to the venom of Calloselasma rhodostoma (Malayan pit viper). To further characterize the complexity composition of the venom, we investigated the proteome of a pooled venom sample from >10 Sri Lankan H. hypnale with reverse-phase high performance liquid chromatography (rp-HPLC), sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and peptide sequencing (tandem mass-spectrometry and/or N-terminal sequencing). The findings ascertained that two phospholipase A2 subtypes (E6-PLA2, W6-PLA2) dominate the toxin composition by 40.1%, followed by snake venom metalloproteases (36.9%), l-amino acid oxidase (11.9%), C-type lectins (5.5%), serine proteases (3.3%) and others (2.3%). The presence of the major toxins correlates with the venom's major pathogenic effects, indicating these to be the principal target toxins for antivenom neutralization. This study supports the previous finding of PLA2 dominance in the venom but diverges from the view that H. hypnale venom has low expression of large enzymatic toxins. The knowledge of the composition and abundance of toxins is essential to elucidate the pathophysiology of H. hypnale envenomation and to optimize antivenom formulation in the future.

A new C-type lectin (RVsnaclec) purified from venom of Daboia russelii russelii shows anticoagulant activity via inhibition of FXa and concentration-dependent differential response to platelets in a Ca²⁺-independent manner

This is the first report on the characterization of a snaclec (RVsnaclec) purified from Daboia russelii russelii venom. The RVsnaclec is a heterodimer of two subunits, α (15.1 kDa) and β (9 kDa). These subunits are covalently linked to form multimeric (αβ)₂ and (αβ)₄ structures. Peptide mass fingerprinting analysis of RVsnaclec via LC-MS/MS demonstrated its similarity to snaclecs purified from other viperid snake venoms. Two tryptic peptide sequences of RVsnaclec revealed the putative conserved domains of C-type lectin (CTL). RVsnaclec dose-dependently increased the Ca-clotting time and prothrombin time of platelet-poor plasma (PPP); however, it did not affect the partial thromboplastin time (APTT) or thrombin time of PPP. The in vitro and in vivo anticoagulant activity of RVsnaclec is correlated to its binding and subsequent uncompetitive inhibition of FXa (Ki = 0.52 μmole) in a Ca(2+)-independent manner; however, supplementation with 0.25 mM Ca(2+) enhanced the Xa binding potency of RVsnaclec. Monovalent or polyvalent antivenom failed to neutralize its anticoagulant potency, and RVsnaclec did not inhibit trypsin, chymotrypsin, thrombin or plasmin. RVsnaclec was devoid of hemolytic activity or cytotoxicity against several human cancer cell lines, demonstrated concentration-dependent aggregation and deaggregation of human platelets, and inhibited the ADP-induced aggregation of platelet. RVsnaclec (5.0 mg/kg body weight) was non-lethal to mice and showed no adverse pharmacological effects, suggesting that it has potential as a lead compound for future therapeutic applications in cardiovascular disorders.

Phosphodiesterase from Daboia russelli russelli venom: purification, partial characterization and inhibition of platelet aggregation

Phosphodiesterases (PDEs) belong to a super-family of enzymes that have multiple roles in the metabolism of extracellular nucleotides and regulation of nucleotide-based intercellular signalling. A PDE from Russell's viper (Daboia russelli russelli) venom (DR-PDE) was purified by gel filtration, ion exchange and affinity chromatographies. Homogeneity of the preparation was verified by SDS-PAGE, SE-HPLC and mass spectrometry. It was free from 5'-nucleotidase, alkaline phosphatase and protease activities. Identity of the enzyme was ensured from partial sequence homology with other PDEs. DR-PDE was inactivated by polyvalent anti-venom serum and metal chelators. The enzyme was partially inhibited by the root extracts of four medicinal plants but remained unaffected by inhibitors of intracellular PDEs. DR-PDE hydrolyses ADP and thus, strongly inhibits ADP-induced platelet aggregation in human platelet rich plasma. This study leads to better understanding of a component of Russell's viper venom that affects homoeostatic system of the victim.

Lebecin, a new C-type lectin like protein from Macrovipera lebetina venom with anti-tumor activity against the breast cancer cell line MDA-MB231

C-type lectins like proteins display various biological activities and are known to affect especially platelet aggregation. Few of them have been reported to have anti-tumor effects. In this study, we have identified and characterized a new C-type lectin like protein, named lebecin. Lebecin is a heterodimeric protein of 30 kDa. The N-terminal amino acid sequences of both subunits were determined by Edman degradation and the entire amino acid sequences were deduced from cDNAs. The precursors of both lebecin subunits contain a 23-amino acid residue signal peptide and the mature α and β subunits are composed of 129 and 131 amino acids, respectively. Lebecin is shown to be a potent inhibitor of MDA-MB231 human breast cancer cells proliferation. Furthermore, lebecin dose-dependently inhibited the integrin-mediated attachment of these cells to different adhesion substrata. This novel C-type lectin also completely blocked MDA-MB231 cells migration towards fibronectin and fibrinogen in haptotaxis assays.

Crystal structures of carbohydrate recognition domain of blood dendritic cell antigen-2 (BDCA2) reveal a common domain-swapped dimer

We report on crystal structures of a carbohydrate recognition domain (CRD) of human C-type lectin receptor blood dendritic cell antigen-2 (BDCA2). Three different crystal forms were obtained at 1.8-2.3 Å resolution. In all three, the CRD has a basic C-type lectin fold, but a long loop extends away from the core domain to form a domain-swapped dimer. The structures of the dimers from the three different crystal forms superimpose well, indicating that domain swapping and dimer formation are energetically stable. The structure of the dimer is compared with other domain-swapped proteins, and a possible regulation mechanism of BDCA2 is discussed.

Quantitative high-throughput profiling of snake venom gland transcriptomes and proteomes (Ovophis okinavensis and Protobothrops flavoviridis)

Background

Advances in DNA sequencing and proteomics have facilitated quantitative comparisons of snake venom composition. Most studies have employed one approach or the other. Here, both Illumina cDNA sequencing and LC/MS were used to compare the transcriptomes and proteomes of two pit vipers, Protobothrops flavoviridis and Ovophis okinavensis, which differ greatly in their biology.

Results

Sequencing of venom gland cDNA produced 104,830 transcripts. The Protobothrops transcriptome contained transcripts for 103 venom-related proteins, while the Ovophis transcriptome contained 95. In both, transcript abundances spanned six orders of magnitude. Mass spectrometry identified peptides from 100% of transcripts that occurred at higher than contaminant (e.g. human keratin) levels, including a number of proteins never before sequenced from snakes. These transcriptomes reveal fundamentally different envenomation strategies. Adult Protobothrops venom promotes hemorrhage, hypotension, incoagulable blood, and prey digestion, consistent with mammalian predation. Ovophis venom composition is less readily interpreted, owing to insufficient pharmacological data for venom serine and metalloproteases, which comprise more than 97.3% of Ovophis transcripts, but only 38.0% of Protobothrops transcripts. Ovophis venom apparently represents a hybrid strategy optimized for frogs and small mammals.

Conclusions

This study illustrates the power of cDNA sequencing combined with MS profiling. The former quantifies transcript composition, allowing detection of novel proteins, but cannot indicate which proteins are actually secreted, as does MS. We show, for the first time, that transcript and peptide abundances are correlated. This means that MS can be used for quantitative, non-invasive venom profiling, which will be beneficial for studies of endangered species.

Glycan evolution in response to collaboration, conflict, and constraint

Glycans, oligo- and polysaccharides secreted or attached to proteins and lipids, cover the surfaces of all cells and have a regulatory capacity and structural diversity beyond any other class of biological molecule. Glycans may have evolved these properties because they mediate cellular interactions and often face pressure to evolve new functions rapidly. We approach this idea two ways. First, we discuss evolutionary innovation. Glycan synthesis, regulation, and mode of chemical interaction influence the spectrum of new forms presented to evolution. Second, we describe the evolutionary conflicts that arise when alleles and individuals interact. Glycan regulation and diversity are integral to these biological negotiations. Glycans are tasked with such an amazing diversity of functions that no study of cellular interaction can begin without considering them. We propose that glycans predominate the cell surface because their physical and chemical properties allow the rapid innovation required of molecules on the frontlines of evolutionary conflict.

Identification of α2β1 integrin inhibitor VP-i with anti-platelet properties in the venom of Vipera palaestinae

Integrins are receptors of the extracellular matrix (ECM), playing a vital role in pathophysiological processes. They bind to ECM ligands like collagens and can mediate wound healing as well as tumor metastasis and thrombosis, thus being a part of cell adhesion and migration as well as platelet aggregation. For this reason, identifying α2β1 integrin-specific antagonists can assist in the development of drugs to treat tumor progression, angiogenesis, and cardiovascular diseases. Snake venoms have been shown to contain antagonists which target collagen-binding integrins. EMS16, rhodocetin, and VP12 are three toxins belonging to the C-type lectin-related protein family and have been proven to inhibit the α2β1 integrin, specifically the α2 integrin A domain. To specifically isolate antagonists targeting the α2β1 integrin A domain, we developed a protocol based on affinity chromatography. Using this novel approach, the toxin VP-i was isolated from Vipera palaestinae venom. We show that VP-i binds to the α2 integrin A domain and that it successfully inhibits adhesion of various cells to type I collagen as well as cell migration. Moreover, our results indicate that VP-i differs structurally from the previously purified VP12, although not functionally, and therefore is a further venom compound which can be utilized for drug development.

Identification and structural characterization of a new three-finger toxin hemachatoxin from Hemachatus haemachatus venom

Snake venoms are rich sources of biologically active proteins and polypeptides. Three-finger toxins are non-enzymatic proteins present in elapid (cobras, kraits, mambas and sea snakes) and colubrid venoms. These proteins contain four conserved disulfide bonds in the core to maintain the three-finger folds. Although all three-finger toxins have similar fold, their biological activities are different. A new three-finger toxin (hemachatoxin) was isolated from Hemachatus haemachatus (Ringhals cobra) venom. Its amino acid sequence was elucidated, and crystal structure was determined at 2.43 Å resolution. The overall fold is similar to other three-finger toxins. The structure and sequence analysis revealed that the fold is maintained by four highly conserved disulfide bonds. It exhibited highest similarity to particularly P-type cardiotoxins that are known to associate and perturb the membrane surface with their lipid binding sites. Also, the increased B value of hemachotoxin loop II suggests that loop II is flexible and may remain flexible until its interaction with membrane phospholipids. Based on the analysis, we predict hemachatoxin to be cardiotoxic/cytotoxic and our future experiments will be directed to characterize the activity of hemachatoxin.

Cloning of cDNAs and molecular characterisation of C-type lectin-like proteins from snake venoms

C-type lectin-like proteins (CTLPs) isolated from snake venoms are the largest and most complex non-mammalian vertebrate C-type lectin-like domain family. In the present study, we simultaneously amplified four cDNAs encoding different types of CTLP subunits from the venoms of two different species of snakes by RT-PCR with a single sense primer and a nested universal primer - two CTLP subunit-encoding cDNAs were cloned from Deinagkistrodon acutus venom and two from Agkistrodon halys Pallas venom. All four cloned CTLP subunits exhibited typical motifs in their corresponding domain regions but with relatively-low sequence similarities to each other. Compared with previously-published CTLPs, the four cloned CTLPs subunits showed slight variations in the calcium-binding sites and the disulphide bonding patterns. To our knowledge, these data constitute the first example of co-expression of CTLP platelet glycoprotein Ib-binding subunits and coagulation factors in Agkistrodon halys Pallas venom.

The continuing saga of snake venom disintegrins

Disintegrins, a family of polypeptides released in the venoms of viperid snakes (vipers and rattlesnakes) by the proteolytic processing of multidomain metalloproteinases, selectively block the function of β(1) and β(3) integrin receptors. Few of the proteins isolated and characterized from snake venoms have proven to be more structural and functional versatile than the disintegrins. Not surprisingly, 25 years after their discovery, our knowledge on the evolutionary history and the molecular determinants modulating the integrin inhibitory activity of disintegrins still remain fragmentary. This paper highlights some seminal contributions, including personal accounts of pioneer authors, related to basic and applied research on disintegrins. Investigators have evaluated disintegrin applications in therapies for a number of pathologies in which integrin receptors play relevant roles, particularly myocardial infarction and inappropriate tumor angiogenesis. Completing the continuing story of the disintegrin family by applying novel research approaches may hold the key to learn how to use deadly toxins as therapeutic agents.

A novel method for the purification of low soluble recombinant C-type lectin proteins

Snake venoms contain a complex mixture of many biological molecules including proteins. The purification of recombinant proteins is a key step in studying their function and structure with affinity chromatography as the common method used in their purification. In bacterial expression systems, hydrophobic recombinant proteins are usually precipitated into inclusion bodies, and contaminants are typically associated with tagged proteins after purification. The purpose of this study was to develop a procedure to purify hydrophobic recombinant proteins without an affinity tag. Snake venom mature C-type lectin-like proteins (CLPs) with a tag were cloned, expressed, and purified by repeated sonication and wash steps. The effects of the signal peptide on the expression and solubility of the recombinant protein were investigated. The CLPs in washed inclusion bodies were solubilized and refolded by dialysis. The CLPs without a tag were successfully purified with a yield 38 times higher than the traditional method, and inhibited blood platelet aggregation with an IC(50) of 100.57 μM in whole blood. This novel procedure is a rapid, and inexpensive method to purify functional recombinant hydrophobic CLPs from snake venoms useful in the development of drug therapies.

Recombinant expression of mutants of the Frankenstein disintegrin, RTS-ocellatusin. Evidence for the independent origin of RGD and KTS/RTS disintegrins

The requirements to transform a short disintegrin of the RGD clade into an RTS disintegrin, were investigated through the generation of recombinant mutants of ocellatusin in which the RGD tripeptide was substituted for RTS in different positions along the integrin-specificity loop. Any attempt to create an active integrin α(1)β(1) inhibitory motif within the specificity loop of ocellatusin was unsuccessful. Replacing the whole RGD-loop of ocellatusin by the RTS-loop of jerdostatin was neither sufficient for confering α(1)β(1) binding specificity to this ocellatusin-RTS Frankenstein(2) mutant. Factors other than the integrin-binding loop sequence per se are thus required to transform a disintegrin scaffold from the RGD clade into another scaffold from the RTS/KTS clade. Moreover, our results provide evidences, that the RTS/KTS short disintegrins have potentially been recruited into the venom gland of Eurasian vipers independently from the canonical neofunctionalization pathway of the RGD disintegrins. PCR-amplifications of jerdostatin-like sequences from a number of taxa across reptiles, including snakes (Crotalinae, Viperinae, and Elapidae taxa) and lizards (Lacertidae and Iguanidae) clearly showed that genes coding for RTS/KTS disintegrins existed long before the split of Lacertidae and Iguania, thus predating the recruitment of the SVMP precursors of disintegrins, providing strong support for the view of an independent evolutionary history of the RTS/KTS and the RGD clades of short disintegrins.

Cloning, characterization and mutagenesis of Russell's viper venom L-amino acid oxidase: Insights into its catalytic mechanism

To investigate the structure-function relationships and geographic variations of L-amino acid oxidase (LAAO) from Daboia venoms, a single LAAO (designated as DrLAO) was purified from eastern Indian Daboia russelii venom and characterized. The purified DrLAO showed subunit molecular mass of 60-64kDa; its N-terminal sequence (1-20) was identical to those of several true viper LAAOs. Its preferred substrates were hydrophobic l-amino acids and the kinetic specificities were ordered as follows: Phe, Tyr, Met, Leu, and Trp. Enzyme assay and Western blotting showed that the venom LAAO contents of D. russelii were higher than those of Daboia siamensis. DrLAO dose-dependently inhibited ADP- and collagen-induced platelet aggregation with IC(50) values of 0.27 and 0.82μM, respectively. Apparently, DrLAO may synergize with other venom components to prolong and enhance bleeding symptoms after Daboia envenoming. The full sequence of DrLAO was deduced from its cDNA sequence and then confirmed by peptide mass fingerprinting. Molecular phylogenetic analysis revealed that SV-LAAO family members could be differentiated not only by snake taxonomy but also by the variations at position 223, and they divided into H223, S223, N223, and D223 subclasses. We have further prepared recombinant DrLAO and mutants by the Pichia expression system. Mutagenic analyses of DrLAO His223 revealed that this residue bound substrates instead of serving as an essential base in the catalytic steps. Our results suggest a direct hydride transfer from substrate to FAD as the mechanism for SV-LAAOs.

SNPs of hemocyanin C-terminal fragment in shrimp Litopenaeus vannamei

In this study, we identified a variable region in the C-terminus of hemocyanin from the shrimp Litopenaeus vannamei (2288-2503bp, HcSC) by sequence alignments. A total of 13 SNPs were identified by PCR-SSCP and HcSC clone sequencing. The SSCP patterns of HcSC could be modulated in Vibro parahaemolyticus-treated shrimps. A novel SSCP band with four SNP sites was identified in V. parahaemolyticus-resistant shrimps. More importantly, three of these four SNPs introduced variations in amino acid sequence and possibly secondary structure of the HcSC polypeptide and resulted in a higher agglutinative activity against seven pathogenic bacteria. These results suggest that the C-terminus of shrimp L. vannamei hemocyanin possesses SNPs, which may be related to shrimp resistance to different pathogens.

Diversity of metalloproteinases in Bothrops neuwiedi snake venom transcripts: evidences for recombination between different classes of SVMPs

Background

Snake venom metalloproteinases (SVMPs) are widely distributed in snake venoms and are versatile toxins, targeting many important elements involved in hemostasis, such as basement membrane proteins, clotting proteins, platelets, endothelial and inflammatory cells. The functional diversity of SVMPs is in part due to the structural organization of different combinations of catalytic, disintegrin, disintegrin-like and cysteine-rich domains, which categorizes SVMPs in 3 classes of precursor molecules (PI, PII and PIII) further divided in 11 subclasses, 6 of them belonging to PII group. This heterogeneity is currently correlated to genetic accelerated evolution and post-translational modifications.

Results

Thirty-one SVMP cDNAs were full length cloned from a single specimen of Bothrops neuwiedi snake, sequenced and grouped in eleven distinct sequences and further analyzed by cladistic analysis. Class P-I and class P-III sequences presented the expected tree topology for fibrinolytic and hemorrhagic SVMPs, respectively. In opposition, three distinct segregations were observed for class P-II sequences. P-IIb showed the typical segregation of class P-II SVMPs. However, P-IIa grouped with class P-I cDNAs presenting a 100% identity in the 365 bp at their 5' ends, suggesting post-transcription events for interclass recombination. In addition, catalytic domain of P-IIx sequences segregated with non-hemorrhagic class P-III SVMPs while their disintegrin domain grouped with other class P-II disintegrin domains suggesting independent evolution of catalytic and disintegrin domains. Complementary regions within cDNA sequences were noted and may participate in recombination either at DNA or RNA levels. Proteins predicted by these cDNAs show the main features of the correspondent classes of SVMP, but P-IIb and P-IIx included two additional cysteines cysteines at the C-termini of the disintegrin domains in positions not yet described.

Conclusions

In B. neuwiedi venom gland, class P-II SVMPs were represented by three different types of transcripts that may have arisen by interclass recombination with P-I and P-III sequences after the divergence of the different classes of SVMPs. Our observations indicate that exon shuffling or post-transcriptional mechanisms may be driving these recombinations generating new functional possibilities for this complex group of snake toxins.

A limited role for gene duplications in the evolution of platypus venom

Gene duplication followed by adaptive selection is believed to be the primary driver of venom evolution. However, to date, no studies have evaluated the importance of gene duplications for venom evolution using a genomic approach. The availability of a sequenced genome and a venom gland transcriptome for the enigmatic platypus provides a unique opportunity to explore the role that gene duplication plays in venom evolution. Here, we identify gene duplication events and correlate them with expressed transcripts in an in-season venom gland. Gene duplicates (1,508) were identified. These duplicated pairs (421), including genes that have undergone multiple rounds of gene duplications, were expressed in the venom gland. The majority of these genes are involved in metabolism and protein synthesis not toxin functions. Twelve secretory genes including serine proteases, metalloproteinases, and protease inhibitors likely to produce symptoms of envenomation such as vasodilation and pain were detected. Only 16 of 107 platypus genes with high similarity to known toxins evolved through gene duplication. Platypus venom C-type natriuretic peptides and nerve growth factor do not possess lineage-specific gene duplicates. Extensive duplications, believed to increase the potency of toxic content and promote toxin diversification, were not found. This is the first study to take a genome-wide approach in order to examine the impact of gene duplication on venom evolution. Our findings support the idea that adaptive selection acts on gene duplicates to drive the independent evolution and functional diversification of similar venom genes in venomous species. However, gene duplications alone do not explain the "venome" of the platypus. Other mechanisms, such as alternative splicing and mutation, may be important in venom innovation.

Adaptive evolution of the venom-targeted vWF protein in opossums that eat pitvipers

The rapid evolution of venom toxin genes is often explained as the result of a biochemical arms race between venomous animals and their prey. However, it is not clear that an arms race analogy is appropriate in this context because there is no published evidence for rapid evolution in genes that might confer toxin resistance among routinely envenomed species. Here we report such evidence from an unusual predator-prey relationship between opossums (Marsupialia: Didelphidae) and pitvipers (Serpentes: Crotalinae). In particular, we found high ratios of replacement to silent substitutions in the gene encoding von Willebrand Factor (vWF), a venom-targeted hemostatic blood protein, in a clade of opossums known to eat pitvipers and to be resistant to their hemorrhagic venom. Observed amino-acid substitutions in venom-resistant opossums include changes in net charge and hydrophobicity that are hypothesized to weaken the bond between vWF and one of its toxic snake-venom ligands, the C-type lectin-like protein botrocetin. Our results provide the first example of rapid adaptive evolution in any venom-targeted molecule, and they support the notion that an evolutionary arms race might be driving the rapid evolution of snake venoms. However, in the arms race implied by our results, venomous snakes are prey, and their venom has a correspondingly defensive function in addition to its usual trophic role.

Cloning and molecular modeling of Litopenaeus vannamei (Penaeidae) C-type lectin homologs with mutated mannose binding domain-2

C-type lectins are animal proteins that contain at least one carbohydrate recognition domain (CRD) capable of mediating sugar and calcium binding. Carbohydrate recognition is directly required for some biological functions, including the innate immune response. We cloned two novel C-type lectin (CTL) precursors from the commercial marine shrimp Litopenaeus vannamei. The cloned cDNAs encompass ORFs of 1044 nucleotides and encode highly similar two-domain polypeptides of 347 residues. The predicted proteins, LvCTL-br1 and -br2, contain the consensus triad that recognizes galactose (-GlnProAsp-) in CRD1 but also contain a mutated mannose-binding site (-GluProAsn-) in the second domain (CRD2). Phylogenetic analysis of LvCTL-br1 and -br2 and hundreds of CTL-like domain-containing proteins have allowed grouping of penaeid shrimp CTLs into three functional clusters. Reverse transcription coupled to PCR indicated that LvCTL-br1 expression is induced in shrimp gills upon IHHNV infection. Computational molecular modeling of LvCTL-br1 and -br2 revealed that three amino acid substitutions in CRD1 occur near the sugar binding site. Also, the 3-D models show a long loop of LvCTL-br1 CRD2 that might accommodate complex sugars. The structural data, evolutionary history and functional analysis support the hypothesis that gene duplication and accelerated evolution have caused functional diversification of penaeid shrimp C-type lectins.

Venom gland transcriptomes of two elapid snakes (Bungarus multicinctus and Naja atra) and evolution of toxin genes

BACKGROUND

Kraits (genus Bungarus) and cobras (genus Naja) are two representative toxic genera of elapids in the old world. Although they are closely related genera and both of their venoms are very toxic, the compositions of their venoms are very different. To unveil their detailed venoms and their evolutionary patterns, we constructed venom gland cDNA libraries and genomic bacterial artificial chromosome (BAC) libraries for Bungarus multicinctus and Naja atra, respectively. We sequenced about 1500 cDNA clones for each of the venom cDNA libraries and screened BAC libraries of the two snakes by blot analysis using four kinds of toxin probes; i.e., three-finger toxin (3FTx), phospholipase A2 (PLA2), kunitz-type protease inhibitor (Kunitz), and natriuretic peptide (NP).

RESULTS

In total, 1092 valid expressed sequences tags (ESTs) for B. multicinctus and 1166 ESTs for N. atra were generated. About 70% of these ESTs can be annotated as snake toxin transcripts. 3FTx (64.5%) and β bungarotoxin (25.1%) comprise the main toxin classes in B. multicinctus, while 3FTx (95.8%) is the dominant toxin in N. atra. We also observed several less abundant venom families in B. multicinctus and N. atra, such as PLA2, C-type lectins, and Kunitz. Peculiarly a cluster of NP precursors with tandem NPs was detected in B. multicinctus. A total of 71 positive toxin BAC clones in B. multicinctus and N. atra were identified using four kinds of toxin probes (3FTx, PLA2, Kunitz, and NP), among which 39 3FTx-positive BACs were sequenced to reveal gene structures of 3FTx toxin genes.

CONCLUSIONS

Based on the toxin ESTs and 3FTx gene sequences, the major components of B. multicinctus venom transcriptome are neurotoxins, including long chain alpha neurotoxins (α-ntx) and the recently originated β bungarotoxin, whereas the N. atra venom transcriptome mainly contains 3FTxs with cytotoxicity and neurotoxicity (short chain α-ntx). The data also revealed that tandem duplications contributed the most to the expansion of toxin multigene families. Analysis of nonsynonymous to synonymous nucleotide substitution rate ratios (dN/dS) indicates that not only multigene toxin families but also other less abundant toxins might have been under rapid diversifying evolution.

Characterisation of a mannose-binding C-type lectin from Oxyuranus scutellatus snake venom

C-type lectins are calcium-dependent sugar binding proteins and are distributed ubiquitously amongst vertebrate organisms. As part of a wider study on Australian snake venom components, we have identified and characterised a C-type lectin from the venom of Oxyuranus scutellatus (Australian coastal taipan) with mannose-binding activity. This protein exhibited a subunit molecular mass of 15 kDa and was found to bind mannose and also bind to and agglutinate erythrocytes in a Ca(2+)-dependent manner. cDNA transcripts coding for C-lectin proteins were cloned and sequenced from six Australian elapid snake species and an antibody generated against the O. scutellatus mannose-binding C-lectin identified C-lectin proteins in the venom of 13 Australian elapid snakes by immunoblotting. Experimental evidence and molecular modelling also suggest that this protein exhibits a unique dimeric structure. This is the first confirmed example of a snake venom C-lectin with mannose-binding activity.

Comparative venom gland transcriptome surveys of the saw-scaled vipers (Viperidae: Echis) reveal substantial intra-family gene diversity and novel venom transcripts

Background

Venom variation occurs at all taxonomical levels and can impact significantly upon the clinical manifestations and efficacy of antivenom therapy following snakebite. Variation in snake venom composition is thought to be subject to strong natural selection as a result of adaptation towards specific diets. Members of the medically important genus Echis exhibit considerable variation in venom composition, which has been demonstrated to co-evolve with evolutionary shifts in diet. We adopt a venom gland transcriptome approach in order to investigate the diversity of toxins in the genus and elucidate the mechanisms which result in prey-specific adaptations of venom composition.

Results

Venom gland transcriptomes were created for E. pyramidum leakeyi, E. coloratus and E. carinatus sochureki by sequencing ~1000 expressed sequence tags from venom gland cDNA libraries. A standardised methodology allowed a comprehensive intra-genus comparison of the venom gland profiles to be undertaken, including the previously described E. ocellatus transcriptome. Blast annotation revealed the presence of snake venom metalloproteinases, C-type lectins, group II phopholipases A₂, serine proteases, L-amino oxidases and growth factors in all transcriptomes throughout the genus. Transcripts encoding disintegrins, cysteine-rich secretory proteins and hyaluronidases were obtained from at least one, but not all, species. A representative group of novel venom transcripts exhibiting similarity to lysosomal acid lipase were identified from the E. coloratus transcriptome, whilst novel metallopeptidases exhibiting similarity to neprilysin and dipeptidyl peptidase III were identified from E. p. leakeyi and E. coloratus respectively.

Conclusion

The comparison of Echis venom gland transcriptomes revealed substantial intrageneric venom variation in representations and cluster numbers of the most abundant venom toxin families. The expression profiles of established toxin groups exhibit little obvious association with venom-related adaptations to diet described from this genus. We suggest therefore that alterations in isoform diversity or transcript expression levels within the major venom protein families are likely to be responsible for prey specificity, rather than differences in the representation of entire toxin families or the recruitment of novel toxin families, although the recruitment of lysosomal acid lipase as a response to vertebrate feeding cannot be excluded. Evidence of marked intrageneric venom variation within the medically important genus Echis strongly advocates further investigations into the medical significance of venom variation in this genus and its impact upon antivenom therapy.

Exploring the venom proteome of the western diamondback rattlesnake, Crotalus atrox, via snake venomics and combinatorial peptide ligand library approaches

We report the proteomic characterization of the venom of the medically important North American western diamondback rattlesnake, Crotalus atrox, using two complementary approaches: snake venomics (to gain an insight of the overall venom proteome), and two solid-phase combinatorial peptide ligand libraries (CPLL), followed by 2D electrophoresis and mass spectrometric characterization of in-gel digested protein bands (to capture and "amplify" low-abundance proteins). The venomics approach revealed approximately 24 distinct proteins belonging to 2 major protein families (snake venom metalloproteinases, SVMP, and serine proteinases), which represent 69.5% of the total venom proteins, 4 medium abundance families (medium-size disintegrin, PLA(2), cysteine-rich secretory protein, and l-amino acid oxidase) amounting to 25.8% of the venom proteins, and 3 minor protein families (vasoactive peptides, endogenous inhibitor of SVMP, and C-type lectin-like). This toxin profile potentially explains the cytotoxic, myotoxic, hemotoxic, and hemorrhagic effects evoked by C. atrox envenomation. Further, our results showing that C. atrox exhibits a similar level of venom variation as Sistrurus miliarius points to a "diversity gain" scenario in the lineage leading to the Sistrurus catenatus taxa. On the other hand, the two combinatorial hexapeptide libraries captured distinct sets of proteins. Although the CPLL-treated samples did not retain a representative venom proteome, protein spots barely, or not at all, detectable in the whole venom were enriched in the two CPLL-treated samples. The amplified low copy number C. atrox venom proteins comprised a C-type lectin-like protein, several PLA(2) molecules, PIII-SVMP isoforms, glutaminyl cyclase isoforms, and a 2-cys peroxiredoxin highly conserved across the animal kingdom. Peroxiredoxin and glutaminyl cyclase may participate, respectively, in redox processes leading to the structural/functional diversification of toxins, and in the N-terminal pyrrolidone carboxylic acid formation required in the maturation of bioactive peptides such as bradykinin-potentiating peptides and endogenous inhibitors of metalloproteases. Our findings underscore the usefulness of combinatorial peptide libraries as powerful tools for mining below the tip of the iceberg, complementing thereby the data gained using the snake venomics protocol toward a complete visualization of the venom proteome.

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.

Role of accelerated segment switch in exons to alter targeting (ASSET) in the molecular evolution of snake venom proteins

Background

Snake venom toxins evolve more rapidly than other proteins through accelerated changes in the protein coding regions. Previously we have shown that accelerated segment switch in exons to alter targeting (ASSET) might play an important role in its functional evolution of viperid three-finger toxins. In this phenomenon, short sequences in exons are radically changed to unrelated sequences and hence affect the folding and functional properties of the toxins.

Results

Here we analyzed other snake venom protein families to elucidate the role of ASSET in their functional evolution. ASSET appears to be involved in the functional evolution of three-finger toxins to a greater extent than in several other venom protein families. ASSET leads to replacement of some of the critical amino acid residues that affect the biological function in three-finger toxins as well as change the conformation of the loop that is involved in binding to specific target sites.

Conclusion

ASSET could lead to novel functions in snake venom proteins. Among snake venom serine proteases, ASSET contributes to changes in three surface segments. One of these segments near the substrate binding region is known to affect substrate specificity, and its exchange may have significant implications for differences in isoform catalytic activity on specific target protein substrates. ASSET therefore plays an important role in functional diversification of snake venom proteins, in addition to accelerated point mutations in the protein coding regions. Accelerated point mutations lead to fine-tuning of target specificity, whereas ASSET leads to large-scale replacement of multiple functionally important residues, resulting in change or gain of functions.

Transcriptomic basis for an antiserum against Micrurus corallinus (coral snake) venom

BACKGROUND

Micrurus corallinus (coral snake) is a tropical forest snake belonging to the family Elapidae. Its venom shows a high neurotoxicity associated with pre- and post-synaptic toxins, causing diaphragm paralysis, which may result in death. In spite of a relatively small incidence of accidents, serum therapy is crucial for those bitten. However, the adequate production of antiserum is hampered by the difficulty in obtaining sufficient amounts of venom from a small snake with demanding breeding conditions. In order to elucidate the molecular basis of this venom and to uncover possible immunogens for an antiserum, we generated expressed sequences tags (ESTs) from its venom glands and analyzed the transcriptomic profile. In addition, their immunogenicity was tested using DNA immunization.

RESULTS

A total of 1438 ESTs were generated and grouped into 611 clusters. Toxin transcripts represented 46% of the total ESTs. The two main toxin classes consisted of three-finger toxins (3FTx) (24%) and phospholipases A(2) (PLA(2)s) (15%). However, 8 other classes of toxins were present, including C-type lectins, natriuretic peptide precursors and even high-molecular mass components such as metalloproteases and L-amino acid oxidases. Each class included an assortment of isoforms, some showing evidence of alternative splicing and domain deletions. Five antigenic candidates were selected (four 3FTx and one PLA(2)) and used for a preliminary study of DNA immunization. The immunological response showed that the sera from the immunized animals were able to recognize the recombinant antigens.

CONCLUSION

Besides an improvement in our knowledge of the composition of coral snake venoms, which are very poorly known when compared to Old World elapids, the expression profile suggests abundant and diversified components that may be used in future antiserum formulation. As recombinant production of venom antigens frequently fails due to complex disulfide arrangements, DNA immunization may be a viable alternative. In fact, the selected candidates provided an initial evidence of the feasibility of this approach, which is less costly and not dependent on the availability of the venom.

Snake venom Vascular Endothelial Growth Factors (VEGF-Fs) exclusively vary their structures and functions among species

Vascular endothelial growth factor (VEGF-A) and its family proteins are crucial regulators of blood vessel formation and vascular permeability. Snake venom has recently been shown to be an exogenous source of unique VEGF (known as VEGF-F), and now, two types of VEGF-F with distinct biochemical properties have been reported. Here, we show that VEGF-Fs (venom-type VEGFs) are highly variable in structure and function among species, in contrast to endogenous tissue-type VEGFs (VEGF-As) of snakes. Although the structures of tissue-type VEGFs are highly conserved among venomous snake species and even among all vertebrates, including humans, those of venom-type VEGFs are extensively variegated, especially in the regions around receptor-binding loops and C-terminal putative coreceptor-binding regions, indicating that highly frequent variations are located around functionally key regions of the proteins. Genetic analyses suggest that venom-type VEGF gene may have developed from a tissue-type gene and that the unique sequence of its C-terminal region was generated by an alteration in the translation frame in the corresponding exons. We further verified that a novel venom-type VEGF from Bitis arietans displays unique properties distinct from already known VEGFs. Our results may provide evidence of a novel mechanism causing the generation of multiple snake toxins and also of a new model of molecular evolution.

Structural requirements of KTS-disintegrins for inhibition of alpha(1)beta(1) integrin

Obtustatin and viperistatin represent the shortest known snake venom monomeric disintegrins. In the present study, we have produced recombinant full-length wild-type and site-directed mutants of obtustatin to assess the role of the K(21)TS(23) tripeptide and C-terminal residues for specific inhibition of the alpha(1)beta(1) integrin. Thr(22) appeared to be the most critical residue for disintegrin activity, whereas substitution of the flanking lysine or serine residues for alanine resulted in a less pronounced decrease in the anti-alpha(1)beta(1) integrin activity of the disintegrin. The triple mutant A(21)AA(23) was devoid of blocking activity towards alpha(1)beta(1) integrin-mediated cell adhesion. The potency of recombinant KTS-disintegrins also depended on the residue C-terminally adjacent to the active motif. Substitution of Leu(24) of wild-type obtustatin for an alanine residue slightly decreased the inhibitory activity of the mutant, whereas an arginine residue in this position enhanced the potency of the mutant over wild-type obtustatin by 6-fold. In addition, the replacements L38V and P40Q may account for a further 25-fold increase in alpha(1)beta(1) inhibitory potency of viperistatin over KTSR-obtustatin.

Evolution of soldier-specific venomous protease in social aphids

In social aphids of the genus Tuberaphis a cysteine protease gene of the family cathepsin B exhibits soldier-specific expression and intestinal protease production. The product is orally excreted and injected by soldier nymphs into natural enemies, thereby exerting an insecticidal activity. In an attempt to gain insights into when and how the novel venomous protease for the altruistic caste has evolved, we investigated the soldier-specific type (S-type) and nonspecific type (N-type) cathepsin B genes from social and nonsocial aphids. All the social aphids examined, representing the genera Tuberaphis, Astegopteryx, and Cerataphis, possessed both the S-type and N-type genes. Phylogenetically distant nonsocial aphids also possessed cathepsin B genes allied to the S-type and the N-type, indicating the evolutionary origin of these genes in the common ancestor of extant aphids. In Tuberaphis species the S-type genes exhibited significant soldier-specific expression and accelerated molecular evolution whereas the N-type genes did not. In Astegopteryx and Cerataphis species, meanwhile, both the S-type and N-type genes exhibited neither remarkable soldier-specific expression nor accelerated molecular evolution. These results suggest that the S-type gene acquired the soldier-specific expression and the venom function after divergence of the genus Tuberaphis. On the structural model of the S-type protease of Tuberaphis styraci the accelerated molecular evolution was associated with the molecular surface rather than the catalytic cleft, suggesting that the venom activity was probably acquired by relatively minor modifications on the molecular surface rather than by generation of a novel active site. In Cerataphis jamuritsu the S-type gene was, although containing a stop codon, structurally almost intact and still transcribed, suggesting recent pseudogenization of the gene copy and possible relevance to relaxed functional constraint in the highly multiplied protease gene family. On the basis of these results we suggest that the massive amplification in aphid cathepsin B genes might have predisposed the evolution of venomous protease in the social aphid lineage and argue that gene duplication, accelerated molecular evolution, and acquisition of novel gene function must have played considerable roles in the evolution of complex biological systems including insect sociality.