Accelerated evolution of snake venom phospholipase A2 isozymes for acquisition of diverse physiological functions

The nucleotide sequences of two cDNAs and four genes encoding Trimeresurus gramineus venom gland phospholipase A2 (PLA2) isozymes were determined and compared internally and externally with those encoding Trimeresurus flavoviridis venom gland PLA2 isozymes. It was revealed that the protein-coding regions are much more diversified than the 5' and 3' untranslated regions (UTRs) and the introns except for the signal peptide domain. The numbers of nucleotide substitutions per site (KN) for the UTRs and the introns were approximately one-quarter of the numbers of nucleotide substitutions per synonymous site (KS) for the protein-coding regions and were at the same level as the KN value of T. gramineus and T. flavoviridis TATA box-binding protein (TBP) genes, indicating that the protein-coding regions of PLA2 isozyme genes are unusually variable and that the UTRs including the introns of venom gland PLA2 isozyme genes have evolved at similar rate to those of non-venomous genes. The numbers of nucleotide substitutions per non-synonymous site (KA) values were close to or larger than the KS values for the protein-coding regions in venom gland PLA2 isozyme genes, indicating that the protein-coding regions of snake venom gland PLA2 isozyme genes have evolved via accelerated evolution. Furthermore, the evolutionary trees derived from the combined sequences of the 5' and 3' UTRs and the signal peptide domain of cDNAs were in accord with the consequences from taxonomy. In contrast, the evolutionary trees from the mature protein-coding region sequences of cDNAs and from the amino acid sequences showed random patterns. Estimations of nucleotide divergence of genes and the phylogenetic analysis reveal that snake venom group IJ PLA2 isozyme genes have been evolving under adaptive pressure to acquire new physiological activities.