The nucleotide sequences of four genes encoding Trimeresurus gramineus (green habu snake, crotalinae) venom gland phospholipase A2 (PLA2; phosphatidylcholine 2-acylhydrolase, EC 3.1.1.4) isozymes were compared internally and externally with those of six genes encoding Trimeresurus flavoviridis (habu snake, crotalinae) venom gland PLA2 isozymes. The numbers of nucleotide substitutions per site (KN) for the noncoding regions including introns were one-third to one-eighth of the numbers of nucleotide substitutions per synonymous site (KS) for the protein-coding regions of exons, indicating that the noncoding regions are much more conserved than the protein-coding regions. The KN values for the introns were found to be nearly equivalent to those of introns of T. gramineus and T. flavoviridis TATA box-binding protein genes, which are assumed to be a general (nonvenomous) gene. Thus, it is evident that the introns of venom gland PLA2 isozyme genes have evolved at a similar rate to those of nonvenomous genes. The numbers of nucleotide substitutions per nonsynonymous site (KA) were close to or larger than the KS values for the protein-coding regions in venom gland PLA2 isozyme genes. All of the data combined reveal that Darwinian-type accelerated evolution has universally occurred only in the protein-coding regions of crotalinae snake venom PLA2 isozyme genes.
Six Trimeresurus flavoviridis (Habu snake) venom gland phospholipase A2 (PLA2) isozyme genes were found to consist of four exons and three introns and to encode proteins of 138 amino acid residues, including the signal sequence of 16 amino acid residues. Comparison of the nucleotide sequences showed that the introns are much more homologous than the protein-coding regions of exons except for the signal peptide-coding region of the first exon. The numbers of nucleotide substitutions per site (KN) for introns are approximately one-fourth of the numbers of nucleotide substitutions per synonymous site (KS) for the protein-coding regions, indicating that the introns are unusually conserved. The absence of an apparent functional role for the introns suggests that the protein-coding regions, except for the signal peptide-coding domains, have evolved at greater substitution rates than introns. The fact that the numbers of nucleotide substitutions per nonsynonymous site (KA) are close to or larger than KS values for relevant pairs of genes revealed that Darwinian-type accelerated substitutions have occurred in the protein-coding regions or exons. This is compatible with the presence of PLA2 species with diverse physiological activities in the venom.