The taxonomic position and the unexpected divergence of the Habu viper, Protobothrops among Japanese subtropical islands.

There are four Habu species currently recognized in Japan: Protobothrops flavoviridis from the Amami Islands and the Okinawa Islands, P. tokarensis from the Tokara Islands, P. elegans from the Yaeyama Islands and Ovophis okinabvensis from the Amami Islands and the Okinawa Islands. To clarify their taxonomic positions, we determined the complete mitochondria genome sequence (approx. 17kb) from two specimens from two different islands each for P. flavoviridis, P. tokarensis and P. elegans as well as one specimen of O. okinavensis and reconstructed the molecular phylogeny of Protobothrops using the published sequences of related species. The maximum likelihood tree showed four major species groups within Protbothrops: Group I consisting of P. cornutus, P. dabieshanensis, P. jerdonii and P. xiangchengensis; Group II consisting of P. flavoviridis and P. tokarensis; Group III consisting of P. maolensis, P. mucrosquamatus and P. elegans; Group IV consisting of P. himalayanus and P. kaubacki. Since we observed an unexpected divergence and the paraphyly of the two samples of P. flavoviridis collected from different islands, Amami-Oshima and Okinawajima within the Group II, we expanded the analysis by increasing the number of P. flavoviridis and P. tokarensis collected from 10 islands: Amami-Oshima (5 specimens), Kakeromajima (4) and Tokunoshima (4) from the Amami Islands, Okinawajima (4), Iheyajima (4), Iejima (4), Tokashikijima (4) and Kumejima (4) from the Okinawa Islands, Kodakarajima (P. tokarensis) (4) and Takarajima (P. tokarensis) (4) from the Tokara Islands. The maximum likelihood tree of the 44 samples replicated the significant divergence of P. flavoviridis between the Amami Clade including Amami-Oshima, Kakeromajima and Tokunoshima and the Okinawa Clade including Okinawajima, Iheyajima, Iejima, Tokashikijima and Kumejima. The Amami Clade also include all specimens from the Tokara Islands currently known as an independent species, P. tokarensis, suggesting the paraphyly of the taxon, P. flavoviridis. In contrast, we observed a distinct lineage of the two specimens from the Yaeyama Islands, supporting the validity of the taxon, P. elegans as an independent species. By MCMC method, we estimated the divergence time between the Amami Clade and the Okinawa Clade to be 6.51MYA, suggesting that the vicariance of the two clades preceded the geological separation of the Amami Islands and the Okinawa Islands (∼1.5MYA). As expected from the limited mobility of terrestrial reptiles including snakes, we observed high genetic divergence in Habu mtDNA among Japanese subtropical island populations.

Structural characteristics and evolution of the Protobothrops elegans pancreatic phospholipase A2 gene in contrast with those of Protobothrops genus venom phospholipase A2 genes.

The nucleotide sequence of the gene encoding Protobothrops elegans (Crotalinae) pancreatic phospholipase A(2) (PLA(2)), abbreviated PePancPLA(2), was determined by means of inverted PCR techniques. Since its deduced amino acid sequence contains a pancreatic loop and shows high similarity to that of Laticauda semifasciata (Elapinae) group IB pancreatic PLA(2), PePancPLA(2) is classified into group IB PLA(2). The nucleotide sequences of the PePancPLA(2) gene, the L. semifasciata group IB pancreatic PLA(2) gene, and the L. semifasciata group IA venom PLA(2) gene are similar to one another but greatly dissimilar to those of Protobothrops genus (Crotalinae) group II venom PLA(2) genes, suggesting that the Elapinae group IB PLA(2) gene and the group IA PLA(2) gene appeared after Elapinae was established, and that the Crotalinae group II venom PLA(2) genes came into existence before Elapinae and Crotalinae diverged. A phylogenetic analysis of their amino acid sequences confirms this.

Structural characteristics and evolution of a novel venom phospholipase A2 gene from Protobothrops flavoviridis.

A novel phospholipase A(2) (PLA(2)) gene, named PfPLA 6, was found in a 6,328-bp NIS-1(5')-a segment in the Protobothrops flavoviridis (Habu, Crotalinae) genome. A comparison of the aligned nucleotide sequences of Viperidae (Viperinae and Crotalinae) venom PLA(2) genes, including PfPLA 6, revealed the deletion of a 12-bp segment called S1EX 1 and a 55-bp segment called S2EX 1 in exon 1 and the interposition of a 219-bp segment called SINT 2 (SINE) in intron 2. A classification of Viperidae PLA(2) genes based on these structural modes indicated that the A-type genes (without SINE), including PfPLA 6, are evolutionarily ancestral to the B-type (Viperinae) and C-type (Crotalinae) PLA(2) genes (both with SINE). Since PfPLA 6 is a pseudogene, an active prototype of PfPLA 6 can be assumed to be the ancestral PLA(2) gene. Putative evolutionary processes from this A-type prototype PLA(2) gene to descendent PLA(2) genes are discussed.

A [Lys49]phospholipase A2 from Protobothrops flavoviridis venom induces caspase-independent apoptotic cell death accompanied by rapid plasma-membrane rupture in human leukemia cells.

Protobothrops flavoviridis venom contains plural phospholipase A(2) (PLA(2)) isozymes. A [Lys(49)]PLA(2) called BPII induced cell death in human leukemia cells. PLA2, an [Asp(49)]PLA(2) that has much stronger lipolytic activity than BPII, failed to induce cell death. BPII-treated cells showed morphological changes, DNA fragmentation, and nuclear condensation. This BPII-induced apoptotic cell death was neither inhibited by inhibitors of caspases 3 and 6 nor accompanied by activation of procaspase 3, indicating that BPII-induced cell death is caspase independent. Since inactive p-bromophenacylated BPII induced cell death, BPII-induced apoptotic cell death is independent of PLA(2) lipolytic activity. Rapid externalization of phosphatidylserine in BPII-treated cells was observed for fluorescein isothiocyanate (FITC)-labeled annexin V. In the cells treated with BPII, this spread over the cell membranes, implying that the cell toxicity of BPII is mediated via its cell-surface receptor.

Identification and evolution of venom phospholipase A2 inhibitors from Protobothrops elegans serum.

The cDNAs encoding venom phospholipase A(2) (PLA(2)) inhibitors (PLIs), named Protobothrops elegans (Pe)γPLI-A, PeγPLI-B, PeαPLI-A, and PeαPLI-B, were cloned from the P. elegans liver cDNA library. They were further divided into several constituents due to nucleotide substitutions in their open reading frames. For PeαPLI-A, two constituents, PeαPLI-A(a) and PeαPLI-A(b), were identified due to three nonsynonymous substitutions in exon 3. Far-western blot and mass-spectrometry analysis of the P. elegans serum proteins showed the presence of γPLIs, and αPLIs, which can bind venom PLA(2)s. In αPLIs from Protobothrops sera, A or B subtype-specific amino acid substitutions are concentrated only in exon 3. A comparison of γPLIs showed that γPLI-As are conserved and γPLI-Bs diversified. Mathematical analysis of the nucleotide sequences of Protobothrops γPLI-B cDNAs revealed that the particular loops in the three-finger motifs diversified by accelerated evolution. Such evolutionary features should have made serum PLIs acquire their respective inhibitory activities to adapt to venom PLA(2) isozymes.

Focused Proteomics Analysis of Habu Snake (Protobothrops flavoviridis) Venom Using Antivenom-Based Affinity Chromatography Reveals Novel Myonecrosis-Enhancing Activity of Thrombin-Like Serine Proteases.

Snakebites are one of the major causes of death and long-term disability in the developing countries due to the presence of various bioactive peptides and proteins in snake venom. In Japan, the venom of the habu snake (Protobothrops flavoviridis) causes severe permanent damage due to its myonecrotic toxins. Antivenom immunoglobulins are an effective therapy for snakebites, and antivenom was recently developed with effective suppressive activity against myonecrosis induced by snake venom. To compare the properties of an antivenom having anti-myonecrotic activity with those of conventional antivenom with no anti-myonecrotic activity, this study applied focused proteomics analysis of habu venom proteins using 2D gel electrophoresis. As a target protein for antivenom immunoglobulins with anti-myonecrotic activity, we identified a thrombin-like serine protease, TLSP2 (TLf2), which was an inactive proteolytic isoform due to the replacement of the active site, His43 with Arg. Additionally, we identified the unique properties and a novel synergistic function of pseudoenzyme TLf2 as a myonecrosis-enhancing factor. To our knowledge, this is the first report of a function of a catalytically inactive snake serine protease.

Discovery of the Gene Encoding a Novel Small Serum Protein (SSP) of Protobothrops flavoviridis and the Evolution of SSPs.

Small serum proteins (SSPs) are low-molecular-weight proteins in snake serum with affinities for various venom proteins. Five SSPs, PfSSP-1 through PfSSP-5, have been reported in Protobothrops flavoviridis ("habu", Pf) serum so far. Recently, we reported that the five genes encoding these PfSSPs are arranged in tandem on a single chromosome. However, the physiological functions and evolutionary origins of the five SSPs remain poorly understood. In a detailed analysis of the habu draft genome, we found a gene encoding a novel SSP, SSP-6. Structural analysis of the genes encoding SSPs and their genomic arrangement revealed the following: (1) SSP-6 forms a third SSP subgroup; (2) SSP-5 and SSP-6 were present in all snake genomes before the divergence of non-venomous and venomous snakes, while SSP-4 was acquired only by venomous snakes; (3) the composition of paralogous SSP genes in snake genomes seems to reflect snake habitat differences; and (4) the evolutionary emergence of SSP genes is probably related to the physiological functions of SSPs, with an initial snake repertoire of SSP-6 and SSP-5. SSP-4 and its derivative, SSP-3, as well as SSP-1 and SSP-2, appear to be venom-related and were acquired later.

Unique structure (construction and configuration) and evolution of the array of small serum protein genes of Protobothrops flavoviridis snake.

The nucleotide sequence of Protobothrops flavoviridis (Pf) 30534 bp genome segment which contains genes encoding small serum proteins (SSPs) was deciphered. The genome segment contained five SSP genes (PfSSPs), PfSSP-4, PfSSP-5, PfSSP-1, PfSSP-2, and PfSSP-3 in this order and had characteristic configuration and constructions of the particular nucleotide sequences inserted. Comparison between the configurations of the inserted chicken repeat-1 (CR1) fragments of P. flavoviridis and Ophiophagus hannah (Oh) showed that the nucleotide segment encompassing from PfSSP-1 to PfSSP-2 was inverted. The inactive form of PfSSP-1, named PfSSP-1δ(Ψ), found in the intergenic region (I-Reg) between PfSSP-5 and PfSSP-1 had also been destroyed by insertions of the plural long interspersed nuclear elements (LINEs) and DNA transposons. The L2 LINE inserted into the third intron or the particular repetitive sequences inserted into the second intron structurally divided five PfSSPs into two subgroups, the Long SSP subgroup of PfSSP-1, PfSSP-2 and PfSSP-5 or the Short SSP subgroup of PfSSP-3 and PfSSP-4 The mathematical analysis also showed that PfSSPs of the Long SSP subgroup evolved alternately in an accelerated and neutral manner, whereas those of the Short SSP subgroup evolved in an accelerated manner. Moreover, the ortholog analysis of SSPs of various snakes showed that the evolutionary emerging order of SSPs was as follows: SSP-5, SSP-4, SSP-2, SSP-1, and SSP-3 The unique interpretation about accelerated evolution and the novel idea that the transposable elements such as LINEs and DNA transposons are involved in maintaining the host genome besides its own transposition natures were proposed.

Alternative mRNA Splicing in Three Venom Families Underlying a Possible Production of Divergent Venom Proteins of the Habu Snake, Protobothrops flavoviridis.

Snake venoms are complex mixtures of toxic proteins encoded by various gene families that function synergistically to incapacitate prey. A huge repertoire of snake venom genes and proteins have been reported, and alternative splicing is suggested to be involved in the production of divergent gene transcripts. However, a genome-wide survey of the transcript repertoire and the extent of alternative splicing still remains to be determined. In this study, the comprehensive analysis of transcriptomes in the venom gland was achieved by using PacBio sequencing. Extensive alternative splicing was observed in three venom protein gene families, metalloproteinase (MP), serine protease (SP), and vascular endothelial growth factors (VEGF). Eleven MP and SP genes and a VEGF gene are expressed as a total of 81, 61, and 8 transcript variants, respectively. In the MP gene family, individual genes are transcribed into different classes of MPs by alternative splicing. We also observed trans-splicing among the clustered SP genes. No other venom genes as well as non-venom counterpart genes exhibited alternative splicing. Our results thus indicate a potential contribution of mRNA alternative and trans-splicing in the production of highly variable transcripts of venom genes in the habu snake.

A novel recombinant system for functional expression of myonecrotic snake phospholipase A(2) in Escherichia coli using a new fusion affinity tag.

A novel recombinant expression system in Escherichia coli was developed using conger eel galectin, namely, congerin II, as an affinity tag. This system was applied for the functional expression of myotoxic lysine-49-phospholipase A(2) ([Lys(49)]PLA(2)), termed BPII and obtained from Protobothrops flavoviridis (Pf) venom. Recombinant Pf BPII fused with a congerin tag has been successfully expressed as a soluble fraction and showed better quantitative yield when folded correctly. The solubility of the recombinant congerin II-tagged BPII increased up to >90% in E. coli strain JM109 when coexpressed with the molecular chaperones GroEL, GroES, and trigger factor (Tf). The tag protein was cleaved by digestion with restriction protease, such as alpha-thrombin or Microbacterium liquefaciens protease (MLP), to obtain completely active recombinant BPII. Thus, the congerin-tagged fusion systems containing the cleavage recognition site for alpha-thrombin or MLP were demonstrated to be highly efficient and useful for producing proteins of desired solubility and activity.

The habu genome reveals accelerated evolution of venom protein genes.

Evolution of novel traits is a challenging subject in biological research. Several snake lineages developed elaborate venom systems to deliver complex protein mixtures for prey capture. To understand mechanisms involved in snake venom evolution, we decoded here the ~1.4-Gb genome of a habu, Protobothrops flavoviridis. We identified 60 snake venom protein genes (SV) and 224 non-venom paralogs (NV), belonging to 18 gene families. Molecular phylogeny reveals early divergence of SV and NV genes, suggesting that one of the four copies generated through two rounds of whole-genome duplication was modified for use as a toxin. Among them, both SV and NV genes in four major components were extensively duplicated after their diversification, but accelerated evolution is evident exclusively in the SV genes. Both venom-related SV and NV genes are significantly enriched in microchromosomes. The present study thus provides a genetic background for evolution of snake venom composition.

Discovery of novel [Arg49]phospholipase A2 isozymes from Protobothrops elegans venom and regional evolution of Crotalinae snake venom phospholipase A2 isozymes in the southwestern islands of Japan and Taiwan.

Protobothrops (formerly Trimeresurus) elegans, a Crotalinae snake, inhabits Ishigaki and Iriomote islands of the Sakishima Islands of Japan which are located between Okinawa island of Japan and Taiwan. Two phospholipase A(2) (PLA(2)) isozymes were purified to homogeneity from P. elegans venom and sequenced. This led to a discovery of novel PLA(2) isozymes with Arg at position 49, that is, [Arg(49)]PLA(2) forms, named PeBP(R)-I and PeBP(R)-II. They are polymorphic at position 3, Val for PeBP(R)-I and Ile for PeBP(R)-II. The cDNAs encoding PeBP(R)-I and PeBP(R)-II were cloned. The cDNA encoding an [Asp(49)]PLA(2) named PePLA(2) was also obtained. In contrast to PLA(2) isozymes from Protobothrops genus with 122 amino acid residues, PeBP(R)-I and PeBP(R)-II are composed of 121 amino acid residues due to lack of Pro at position 90. They exhibited necrotic and edema-inducing activities but no hemorrhagic activity was detected. A phylogenetic tree constructed for venom PLA(2) isozymes of Protobothrops genus and of related genera in the southwestern islands of Japan and Taiwan revealed that PeBP(R)-I and PeBP(R)-II of P. elegans are evolutionarily much closer to PmK49PLA(2), a [Lys(49)]PLA(2), from P. mucrosquamatus (Taiwan) than BPI and BPII, both [Lys(49)]PLA(2) forms, from P. flavoviridis (Amami-Oshima and Tokunoshima islands of Japan). Such evolutionary relationships are also seen in neutral [Asp(49)]PLA(2) isozymes from the three Protobothrops species. Thus, P. elegans is the species much closer to P. mucrosquamatus than P. flavoviridis. Their evolutionary distances seem to be well related to geological history of the islands where they have lived. In addition, it was clearly noted that Ovophis okinavensis (Amami-Oshima), which had formerly belonged to the Trimeresurus genus, and Trimeresurus stejnegeri (Taiwan) are the species fairly distant from Protobothrops genus.

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

A number of C-type lectin-like proteins that affect thrombosis and hemostasis by inhibiting or activating specific platelet membrane receptors or blood coagulation factors have been isolated from the venom of various snake species and characterized and more than 80 have been sequenced. Recent data on the primary sequences and 3D structures of C-type lectins and C-type lectin-like proteins from snake venoms have enabled us to analyze their molecular evolution. Statistical analysis of their cDNA sequences shows that C-type lectin-like proteins, with some exceptions, have evolved in an accelerated manner to acquire their diverse functions. Phylogenetic analysis shows that the A and B chains of C-type lectin-like proteins are clearly separated from C-type lectins and that the A and B chains are further divided into a group of platelet receptor-binding proteins and a group of coagulation factor-binding proteins. Elucidation of the tertiary structures of several C-type lectin-like proteins led to the discovery of a unique domain-swapping interaction between heterodimeric subunits, which creates a concave surface for ligand binding.

Amino acid sequence of a basic aspartate-49-phospholipase A2 from Trimeresurus flavoviridis venom and phylogenetic analysis of Crotalinae venom phospholipases A2.

Trimeresurus flavoviridis snakes inhabit the southwestern islands of Japan: Amami-Oshima, Tokunoshima and Okinawa. A phospholipase A2 (PLA2) of basic nature (pI 8.5) was isolated from the venom of Amami-Oshima T. flavoviridis. Its amino acid sequence determined by the ordinary procedures was completely in accord with that predicted from the nucleotide sequence of the cDNA previously cloned from Amami-Oshima T. flavoviridis venom gland, which was named PLA-B'. It consists of 122 amino acid residues and has aspartate at position 49. It induced edema in a mouse footpad assay and caused necrosis in mouse skeletal muscles. PLA-B' is similar in sequence to PLA-B (Tokunoshima) and PL-Y (Okinawa), both basic [Asp49]PLA2s, with a few amino acid substitutions, indicating occurrence of interisland mutation. Although PLA2s of Crotalinae subfamily were phylogenetically classified into four types, PLA2 (acidic or neutral [Asp49]PLA2) type, basic [Asp49]PLA2 type, neurotoxic [Asp49]PLA2 type and [Lys49]PLA2 type, it was ascertained that PLA2s of PLA2 type and [Lys49]PLA2 type are most essential as toxic components for Crotalinae snake venoms and that basic [Asp49]PLA2-type PLA2s are uniquely contained only in the venoms of T. flavoviridis species. Prediction of physiological activities of some PLA2s was made based on their location in the phylogenetic tree. Relationship of divergence of PLA2s via accelerated evolution followed by less rapid mutation and physiological activities was discussed.

Properties and cDNA cloning of an antihemorrhagic factor (HSF) purified from the serum of Trimeresurus flavoviridis.

Habu serum factor (HSF) is a metalloproteinase inhibitor that is isolated from the serum of habu snake (Trimeresurus flavoviridis), and it can suppress snake venom-induced hemorrhage. In the present study, the inhibitory property and fundamental structure of HSF were analyzed in detail. HSF inhibited all the hemorrhagic and most of the non-hemorrhagic metalloproteinases tested from the venoms of T. flavoviridis and Gloydius halys brevicaudus. HSF was extremely stable in a broad range of temperature and pH, and the treatments with a temperature of 100 degrees C or pH ranging from 1 to 13 barely affects its reactivity against G. halys brevicaudus H6 protease. Gel filtration chromatography revealed that HSF binds to the H6 protease with a 1:1 molar ratio. A secondary structure profile of HSF that was monitored by circular dichroism spectrum remained unvaried up to 2 M urea. The activity of HSF was stoichiometrically abolished by chemical modification with 2,4,6-trinitrobenzene sulfonic acid and N-bromosuccinimide; this indicates that Lys and Trp residues in its sequence play a role in the inhibitory mechanism. In this study, the amino acid sequence of HSF that was obtained by cDNA cloning was identical to that reported previously, except for five substitutions. We concluded that these discrepancies reflect a difference in the places of capture of the snake specimens.

Purification, primary structures and evolution of coagulant proteases from Deinagkistrodon actus venom.

Deinagkistrodon (formerly Agkistrodon) actus (Taiwan) snake venom was found to contain at least seven closely related coagulant proteases. One of them, named actibin, was purified to homogeneity by means of four chromatographic steps. Actibin acted on fibrinogen to form fibrin clots with extremely high specific activity of 1,630 NIH units/mg and preferentially released fibrinopeptide A. Actibin was an acidic glycoprotein (pI 3.4) with molecular weight of 41,000, which was reduced to 28,800 after deglycosylation with N-glycanase. The k(cat)/K(m) values of actibin for hydrolysis of tosyl-l-arginine methyl ester and benzoyl-l-arginine p-nitroanilide were one-third to a half those for thrombin, reflecting a high potency of actibin in fibrinogen clotting. The amidase activities of actibin and its family proteases were inhibited by 3,4-dichloroisocoumarin, a serine protease inhibitor, indicating that actibin and its family proteases are serine proteases. Four cDNAs, named DaP1 and DaP7-DaP9, encoding D. actus coagulant proteases were cloned. All cDNAs contain an open reading frame of 780 bp coding for 260 amino acid residues, including a signal peptide of 24 amino acid residues. Their amino acid sequences predicted are highly homologous to one another with one to five amino acid substitutions. When four D. actus protease cDNAs were compared with the cDNAs coding for Trimeresurus flavoviridis and T. gramineus venom serine proteases, accelerated evolution was clearly observed. Similarity of the nucleotide sequences of four D. actus protease cDNAs with no synonymous and one to five nonsynonymous substitutions seems not to be in direct conformity with accelerated evolution. This possibly suggests that they have evolved to a similar direction to enhance their clotting activity rather than to produce other physiological activities.

Interisland variegation of venom [Lys(49)]phospholipase A2 isozyme genes in Protobothrops genus snakes in the southwestern islands of Japan.

Protobothrops tokarensis (Pt), a Crotalinae snake, inhabits only Takarajima and Kodakarajima islands of the Tokara Islands located in the immediate north of Amami-Oshima island of Japan. Kodakarajima P. tokarensis venom gland cDNA library gave four types of phospholipase A2 (PLA2) cDNAs encoding neutral [Asp(49)]PLA2, basic [Asp(49)]PLA2, highly basic [Asp(49)]PLA2, and [Lys(49)]PLA2. As the amino acid sequences encoded by their open reading frames (ORFs) were identical to those of PLA2, PLA-B, PLA-N, and BPI (a [Lys(49)]PLA2), respectively, from Amami-Oshima P. flavoviridis (Pf) venom, they were named PtPLA2, PtPLA-B, PtPLA-N, and PtBPI. Chromatography of P. tokarensis venom gave three PLA2 isozymes, PtPLA2, PtPLA-B, and PtBPI. However, BPII and BPIII ([Lys(49)]PLA2s) expressed in Amami-Oshima P. flavoviridis venom were not found in P. tokarensis venom. Genomic polymerase chain reaction (PCR) for P. tokarensis liver DNAs with the unique primers gave PtBPI gene. Notably it was found that LINE (long interspersed nuclear element)-1 fragment is inserted into second intron of PtBPI gene. The LINE-1 fragment may prevent duplication of PtBPI gene and thus formation of plural [Lys(49)]PLA2 genes in P. tokarensis genome. The interisland variegation of venom [Lys(49)]PLA2 isozyme genes in Protobothrops genus snakes in the southwestern islands of Japan is discussed.

Molecular evolution of myotoxic phospholipases A2 from snake venom.

After two decades of study, we draw the conclusion that venom-gland phospholipase A2 (PLA2) isozymes, including PLA2 myotoxins of Crotalinae snakes, have evolved in an accelerated manner to acquire their diverse physiological activities. In this review, we describe how accelerated evolution of venom PLA2 isozymes was discovered. This type of evolution is fundamental for other venom isozyme systems. Accelerated evolution of venom PLA2 isozyme genes is due to rapid change in exons, but not in introns and the flanking regions, being completely opposite to the case of the ordinary isozyme genes. The molecular mechanism by which proper base substitutions had occurred in the particular sites of venom isozyme genes is a puzzle to be solved in future studies. It should be noted that accelerated evolution occurred until the isozymes had acquired their particular function and, since then, they have evolved with less frequent mutation, possibly for functional conservation. We also found that interisland mutations occurred in venom PLA2 isozymes. The relationships between mutation and its driving force are speculative and the real mechanism remains a mystery.

Characterization, primary structure and molecular evolution of anticoagulant protein from Agkistrodon actus venom.

An anticoagulant protein named AaACP was isolated from Agkistrodon actus (hundred-pace snake of Taiwan, Viperidae) venom. AaACP inhibited the factor Xa-induced plasma coagulation in a concentration-dependent manner. Thus, AaACP seems to bind to factor Xa in prothrombinase complex. AaACP was composed of A and B chains linked by disulphide bond(s). The amino acid sequences of A and B chains of AaACP were analysed with a few residues unidentified which were complemented from the nucleotide sequences of their cDNAs. The A chain consisted of 129 amino acid residues and the B chain 123 amino acid residues. Their amino acid sequences were highly similar to those of A and B chains of a series of anticoagulant proteins which had been purified from the venoms of some Viperidae snakes. The A and B chains structurally belong to C-type lectin-like protein family of snake venom origin. Construction of phylogenetic tree of C-type lectins and C-type lectin-like proteins based on their amino acid sequences indicated that their A and B chains diverged before speciation of snake species. The comparison of the nucleotide sequences of the cDNAs encoding A and B chains of AaACP and of Trimeresurus flavoviridis (Viperidae) venom-gland factors IX/X-binding protein and factor IX-binding protein showed that the mature protein-coding region is much more variable than the signal peptide-coding domain and the 5'- and 3'-untranslated regions, being in contrast to the case of the ordinary isoprotein genes. The ratios of the numbers of nucleotide substitutions per nonsynonymous site (K(A)) and per synonymous site (K(S)) in the mature protein-coding region in the cDNA pairs were about three times greater than those for the ordinary isoprotein genes, suggesting that these genes have been evolving in an accelerated manner. Taking account of the functional diversities of venom-gland C-type lectins and C-type lectin-like proteins including factors IX and/or X-binding proteins, it can be said that their functional diversities have been acquired by accelerated evolution.

Discovery of a novel vascular endothelial growth factor (VEGF) with no affinity to heparin in Gloydius tsushimaensis venom.

Strong vascular permeability enhancing activity was found only in the venom of Gloydius tsushimaensis, in Tsushima island, Japan, when examined together with the venoms of G. blomhoffii snakes in several areas of Japan and of G. ussuriensis in South Korea. The active protein purified by using Superdex 75 and Mono Q columns showed no affinity to heparin, and migrated on SDS-PAGE with molecular weights of 26 and 13 kDa under nonreducing and reducing conditions, respectively, showing that it exists as homodimer. Its N-terminal amino acid sequence was highly homologous to those of snake venom vascular endothelial growth factors (VEGFs). The sequence of this protein, named GtVF, was inferred from the one base-substituted two cDNAs (438 bp) obtained via the 3' RACE. The phylogenetic analysis suggested the presence of a new type of snake venom VEGFs including GtVF with no affinity to heparin in addition to the known three types of snake venom VEGFs with high affinity to heparin. Since the vascular permeability enhancement by GtVF was inhibited by the antibody against kinase insert domain-containing receptor (KDR), the vascular permeability enhancing activity of GtVF arises through KDR but without heparin binding.