Structures of genes encoding TATA box-binding proteins from Trimeresurus gramineus and T. flavoviridis snakes

Secreted Phospholipases A2: Drivers of Inflammation and Cancer

Secreted phospholipase 2 (sPLA2) is the largest family of phospholipase A2 (PLA2) enzymes with 11 mammalian isoforms. Each sPLA2 exhibits different localizations and specific properties, being involved in a very wide spectrum of biological processes. The enzymatic activity of sPLA2 has been well described; however, recent findings have shown that they could regulate different signaling pathways by acting directly as ligands. Arachidonic acid (AA) and its derivatives are produced by sPLA2 in collaboration with other molecules in the extracellular space, making important impacts on the cellular environment, being especially relevant in the contexts of immunity and cancer. For these reasons, this review focuses on sPLA2 functions in processes such as the promotion of EMT, angiogenesis, and immunomodulation in the context of tumor initiation and progression. Finally, we will also describe how this knowledge has been applied in the search for new sPLA2 inhibitory compounds that can be used for cancer treatment.

Molecular evolution of toxin genes in Elapidae snakes

The venom of the sea krait, Laticauda semifasciata, consists primarily of two toxic proteins, phospholipase A(2) (PLA(2)) and a three-finger-structure toxin. We have cloned both toxic protein genes, including the upstream region. PLA(2) genes contain three types of inserted sequences: an AG-rich region, a chicken repeat 1-like long interspersed nucleotide element sequence and an intron II 3' side repeat sequence. The molecular divergence of L. semifasciata PLA(2) genes was defined on the basis of the inserted sequences and their sequence homology. The length of intron I in the three-finger-structure toxin genes differs from species to species. The alignment analysis of intron I of the three-finger-structure toxin genes revealed that the intron I sequence of the ancestral gene comprised ten genetic regions. A hypothetical evolutionary process for the three-finger-structure toxin genes has also been developed.

Optimal Intron Analyses in the Trimeresurus Radiation of Asian Pitvipers

Nuclear introns are commonly used as phylogenetic markers, but a number of issues related to alignment strategies, indel treatments, and the incorporation of length-variant heterozygotes (LVHs) are not routinely addressed when generating phylogenetic hypotheses. Topological congruence in relation to an extensive mitochondrial DNA multigene phylogeny (derived from 2,423 bp of 12S, 16S, ND4, and CYTB genes) of the Asian pitviper Trimeresurus radiation was used to compare combinations of "by eye" and edited and unedited ClustalX 1.8 alignments of two nuclear introns. Indels were treated as missing data, fifth character states, and assigned simple and multistate codes. Upon recovery of the optimal alignment and indel treatment strategy, a total evidence approach was used to investigate the phylogenetic utility of the indels and test new generic arrangements within Trimeresurus. Approximately one third of the intron data partitions exhibited LVHs, suggesting that they are common in introns. Furthermore, a simple concatenation approach can facilitate the incorporation of LVHs into phylogenetic analyses to make use of all available data and investigate mechanisms of molecular evolution. Analyses of ClustalX 1.8-assisted alignments were generally more congruent than the "by eye" alignment and the analysis of a simple coded, edited ClustalX 1.8 (gap opening cost 5, gap extension cost 1) alignment revealed the most congruent tree. The total evidence approach supported the new arrangements within Trimeresurus, suggesting that the phylogeny should be considered as a working benchmark in Asian pitviper systematics. Finally, a critical appraisal of the diverse array of indels (56 to 57 per intron, ranging from 1 to 151 bp in length) suggested that they are a combination of Hennigian and homoplasious events unrelated to indel size or location within the intron. [Alignment; indels; intron analysis; length-variant heterozygotes; Trimeresurus.].

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.

Early vertebrate evolution of the TATA-binding protein, TBP

TBP functions in transcription initiation in all eukaryotes and in Archaebacteria. Although the 181-amino acid (aa) carboxyl (C-) terminal core of the protein is highly conserved, TBP proteins from different phyla exhibit diverse sequences in their amino (N-) terminal region. In mice, the TBP N-terminus plays a role in protecting the placenta from maternal rejection; however the presence of similar TBP N-termini in nontherian tetrapods suggests that this domain also has more primitive functions. To gain insights into the pretherian functions of the N-terminus, we investigated its phylogenetic distribution. TBP cDNAs were isolated from representative nontetrapod jawed vertebrates (zebrafish and shark), from more primitive jawless vertebrates (lamprey and hagfish), and from a prevertebrate cephalochordate (amphioxus). Results showed that the tetrapod N-terminus likely arose coincident with the earliest vertebrates. The primary structures of vertebrate N-termini indicates that, historically, this domain has undergone events involving intragenic duplication and modification of short oligopeptide-encoding DNA sequences, which might have provided a mechanism of de novo evolution of this polypeptide.

Sequences and structural organization of phospholipase A2 genes from Vipera aspis aspis, V. aspis zinnikeri and Vipera berus berus venom. Identification of the origin of a new viper population based on ammodytin I1 heterogeneity

We used a PCR-based method to determine the genomic DNA sequences encoding phospholipases A2 (PLA2s) from the venoms of Vipera aspis aspis (V. a. aspis), Vipera aspis zinnikeri (V. a. zinnikeri), Vipera berus berus (V. b. berus) and a neurotoxic V. a. aspis snake (neurotoxic V. a. aspis) from a population responsible for unusual neurotoxic envenomations in south-east France. We sequenced five groups of genes, each corresponding to a different PLA2. The genes encoding the A and B chains of vaspin from the neurotoxic V. a. aspis, PLA2-I from V. a. zinnikeri, and the anticoagulant PLA2 from V. b. berus are described here. Single nucleotide differences leading to amino-acid substitutions were observed both between genes encoding the same PLA2 and between genes encoding different PLA2s. These differences were clustered in exons 3 and 5, potentially altering the biological activities of PLA2. The distribution and characteristics of the PLA2 genes differed according to the species or subspecies. We characterized for the first time genes encoding neurotoxins from the V. a. aspis and V. b. berus snakes of central France. Genes encoding ammodytins I1 and I2, described previously in Vipera ammodytes ammodytes (V. am. ammodytes), were also present in V. a. aspis and V. b. berus. Three different ammodytin I1 gene sequences were characterized: one from V. b. berus, the second from V. a. aspis, V. a. zinnikeri and the neurotoxic V. a. aspis, and the third from the neurotoxic V. a. aspis. This third sequence was identical with the reported sequence of the V. am. ammodytes ammodytin I1 gene. Genes encoding monomeric neurotoxins of V. am. ammodytes venom, ammodytoxins A, B and C, and the Bov-B LINE retroposon, a phylogenetic marker found in V. am. ammodytes genome, were identified in the genome of the neurotoxic V. a. aspis. These results suggest that the population of neurotoxic V. a. aspis snakes from south-east France may have resulted from interbreeding between V. a. aspis and V. am. ammodytes.

Fundamental cellular processes do not require vertebrate-specific sequences within the TATA-binding protein

The 180-amino acid core of the TATA-binding protein (TBPcore) is conserved from Archae bacteria to man. Vertebrate TBPs contain, in addition, a large and highly conserved N-terminal region that is not found in other phyla. We have generated a line of mice in which the tbp allele is replaced with a version, tbp(Delta N), which lacks 111 of 135 N-terminal amino acid residues. Most tbp(Delta N/Delta N) fetuses die in midgestation. To test whether a disruption of general cellular processes contributed to this fetal loss, primary fibroblast cultures were established from +/+, Delta N/+, and Delta N/Delta N fetuses. The cultures exhibited no genotype-dependent differences in proliferation or in expression of the proliferative markers dihydrofolate reductase (DHFR) mRNA (S phase-specific) and cdc25B mRNA (G(2)-specific). The mutation had no effect on transcription initiation site fidelity by either RNA polymerase II (pol II) or pol III. Moreover, the mutation did not cause differences in levels of U6 RNA, a pol III-dependent component of the splicing machinery, in mRNA splicing efficiency, in expression of housekeeping genes from either TATA-containing or TATA-less promoters, or in global gene expression. Our results indicated that general eukaryotic cell functions are unaffected by deletion of these vertebrate-specific sequences from TBP. Thus, all activities of this polypeptide domain must either be compensated for by redundant activities or be restricted to situations that are not represented by primary fibroblasts.

A comparative analysis of invaded sequences from group IA phospholipase A(2) genes provides evidence about the divergence period of genes groups and snake families

Two phospholipase A(2) (PLA(2)) genes classified into group IA were cloned from the genomic library of the sea snake Laticauda semifasciata. Eight clones were obtained by PCR cloning procedure from genomic DNA of Laticauda laticaudata (four clones) and Laticauda colubrina (four clones). The genes were 3.6-4.4kbp in length. Intron and exon organization of the group IA PLA(2) genes was the same as that of Naja sputatrix group IA PLA(2) genes (four exons and three introns). There were two kinds of repetitive sequences in the first and second introns of all sequenced PLA(2) genes. The differences in the length of these genes were derived from the length of their repetitive sequences. The chicken repeat-1 (CR1)-like long interspersed repeated DNA (LINE) sequences, different from the above repetitive sequences, were also found in all sequenced Laticauda PLA(2) genes. A comparative analysis of groups IA, IA' and IIA PLA(2)s genes suggests a period of CR1-like LINE integration during molecular and family evolution. The integration of CR1-like LINE into PLA(2) genes occurred after the divergence of groups I and II PLA(2)s but before the divergence of groups, IA and IA' PLA(2)s. These integration events occurred before the family divergence of Naja and Laticauda. The presence of CR1-like LINE and a comparison of intron and exon organization showed that the divergence of Naja and Bungarus occurred before the divergence of Laticauda and Naja.

The binding interaction of HMG-1 with the TATA-binding protein/TATA complex

High mobility protein-1 (HMG-1) has been shown to regulate transcription by RNA polymerase II. In the context that it acts as a transcriptional repressor, it binds to the TATA-binding protein (TBP) to form the HMG-1/TBP/TATA complex, which is proposed to inhibit the assembly of the preinitiation complex. By using electrophoretic mobility shift assays, we show that the acidic C-terminal domain of HMG-1 and the N terminus of human TBP are the domains that are essential for the formation of a stable HMG-1/TBP/TATA complex. HMG-1 binding increases the affinity of TBP for the TATA element by 20-fold, which is reflected in a significant stimulation of the rate of TBP binding, with little effect on the dissociation rate constant. In support of the binding target of HMG-1 being the N terminus of hTBP, the N-terminal polypeptide of human TBP competes with and inhibits HMG-1/TBP/TATA complex formation. Deletion of segments of the N terminus of human TBP was used to map the region(s) where HMG-1 binds. These findings indicate that interaction of HMG-1 with the Q-tract (amino acids 55-95) in hTBP is primarily responsible for stable complex formation. In addition, HMG-1 and the monoclonal antibody, 1C2, specific to the Q-tract, compete for the same site. Furthermore, calf thymus HMG-1 forms a stable complex with the TBP/TATA complex that contains TBP from either human or Drosophila but not yeast. This is again consistent with the importance of the Q-tract for this stable interaction and shows that the interaction extends over many species but does not include yeast TBP.

Analysis of the chicken TBP-like protein(tlp) gene: evidence for a striking conservation of vertebrate TLPs and for a close relationship between vertebrate tbp and tlp genes

TLP (TBP-like protein), which is a new protein dis-covered by us, has a structure similar to that of the C-terminal conserved domain (CCD) of TBP, although its function has not yet been elucidated. We isolated cDNA and genomic DNA that encode chicken TLP (cTLP) and determined their structures. The predicted amino acid sequence of cTLP was 98 and 91% identical to that of its mammalian and Xenopus counterparts, respectively, and its translation product was ubiquitously observed in chicken tissues. FISH detection showed that chicken tlp and tbp genes were mapped at 3q2.6-2.8 and 3q2.4-2.6 of the same chromosome, respectively. Genome analysis revealed that the chicken tlp gene was spliced with five introns. Interestingly, the vertebrate tbp genes were also found to be split by five introns when we focused on the CCDs, and their splicing points were similar to those of tlp. On the contrary, another TBP-resembling gene of Drosophila, trf1, is split by only one intron, as is the Drosophila 's tbp gene. These results support our earlier assumption that vertebrate TLPs did not directly descend from Drosophila TRF1. On the basis of these results together with phylogenetical exam-ination, we speculate that tlp diverged from an ancestral tbp gene through a process of gene duplication and point mutations.

Isolation and characterization of the gene coding for Artemia franciscana TATA-binding protein: expression in cryptobiotic and developing embryos

Genomic and cDNA clones coding for the Artemia franciscana homolog of the TATA box-binding protein (TBP) were isolated. The C-terminal region of the predicted protein displays up to 92% sequence identity with the conserved C-terminal regions of TBPs from other species. The gene is divided in seven exons that expand over a region of 33 kb. The position of the four introns located in the conserved C-terminal region has been compared with those of other species. Two of these introns have been generally conserved during evolution, another is an arthropod specific intron, present in Drosophila melanogaster and A. franciscana, and the other is only conserved between vertebrates and A. franciscana. Primer extension experiments detected several transcription initiation sites. Northern blot analyses showed the presence of four mRNAs of estimated sizes of 6.8, 2.6, 1.6 and 1.1 kb. Except for the low expression of the 6.8 and 2. 6 kb RNAs in encysted embryos, steady-state levels showed little variation during the activation of the encysted embryo and the first steps of embryonic and larval development. The amount of TBP protein expressed in encysted embryos and developing larvae has been analyzed by Western blot. Cryptobiotic embryos contain significant amounts of TBP although the level of expression increased almost twice during the first 20 h of development. The presence of TBP protein in cryptobiotic embryos suggests that TBP does not play, by itself, a critical role in the arrest of transcription characteristic of these resistance forms.

The Bov-B lines found in Vipera ammodytes toxic PLA2 genes are widespread in snake genomes

In the fourth intron of two toxic Vipera ammodytes PLA2 genes a Ruminantia specific 5'-truncated Bov-B LINE element was identified. Southern blot analysis of Bov-B LINE distribution in vertebrates shows that, apart from the Ruminantia, it is limited to Viperidae snakes (V. ammodytes, Vipera palaestinae, Echis coloratus, Bothrops alternatus, Trimeresurus flavoviridis and Trimeresurus gramineus). The copy number of the 3' end of Bov-B LINE in the V. ammodytes genome is between 62,000 and 75,000. At orthologous positions in other snake PLA2 genes the Bov-B LINE element is absent, indicating that its retrotransposition in the V. ammodytes PLA2 gene locus has occurred quite recently, about 5 Myr ago. The amplification of Bov-B LINEs in snakes may have occurred before the divergence of the Viperinae and Crotalinae subfamilies. Due to its wide distribution in Viperidae snakes it should be a valuable phylogenetic marker. The neighbour-joining phylogenetic tree shows two clusters of truncated Bov-B LINE, a Bovidae and a snake cluster, indicating an early horizontal transfer of this transposable element.

Retrotransposable CR1-like elements in crotalinae snake genomes

A part of the 3'-flanking region of BP-II gene, which is one of Trimeresurus flavoviridis venom gland phospholopase A2 (PLA2) isozyme genes, has a region homologous to avian chicken repeat 1 (CR1)-element. In the present study, ten CR1-like elements were further identified in T. gramineus venom gland PLA2 isozyme genes, T. flavoviridis PLA2 inhibitor (PLI) genes, and T. flavoviridis and T. gramineus TATA-box binding protein (TBP) genes. Southern blot analysis using a probe for CR1 showed that Crotalinae snake genomes contain a number of CR1-like elements.

Molecular evolution of snake toxins: is the functional diversity of snake toxins associated with a mechanism of accelerated evolution?

Recent studies revealed that animal toxins with unrelated biological functions often possess a similar architecture. To tentatively understand the evolutionary mechanisms that may govern this principle of functional prodigality associated with a structural economy, two complementary approaches were considered. One of them consisted of investigating the rates of mutations that occur in cDNAs and/or genes that encode a variety of toxins with the same fold. This approach was largely adopted with phospholipases A2 from Viperidae and to a lesser extent with three-fingered toxins from Elapidae and Hydrophiidae. Another approach consisted of investigating how a given fold can accommodate distinct functional topographies. Thus, a number of topologies by which three-fingered toxins exert distinct functions were investigated either by making chemical modifications and/or mutational analyses or by studying the three-dimensional structure of toxin-target complexes. This review shows that, although the two approaches are different, they commonly indicate that most if not all the surface of a snake toxin fold undergoes natural engineering, which may be associated with an accelerated rate of evolution. The biochemical process by which this phenomenon occurs remains unknown.

Bov-B long interspersed repeated DNA (LINE) sequences are present in Vipera ammodytes phospholipase A2 genes and in genomes of Viperidae snakes

Ammodytin L is a myotoxic Ser49 phospholipase A2 (PLA2) homologue, which is tissue-specifically expressed in the venom glands of Vipera ammodytes. The complete DNA sequence of the gene and its 5' and 3' flanking regions has been determined. The gene consists of five exons separated by four introns. Comparative analysis of the ammodytin L and ammodytoxin C genes shows that all intron and flanking sequences are considerably more conserved (93-97%) than the mature protein-coding exons. The pattern of nucleotide substitutions in protein-coding exons is not random but occurs preferentially on the first and the second positions of codons, which suggests positive Darwinian evolution for a new function. An Ruminantia specific ART-2 retroposon, recently recognised as a 5'-truncated Bov-B long interspersed repeated DNA (LINE) sequence, was identified in the fourth intron of both genes. This result suggests that ammodytin L and ammodytoxin C genes are derived by duplication of a common ancestral gene. The phylogenetic distribution of Bov-B LINE among vertebrate classes shows that, besides the Ruminantia, it is limited to Viperidae snakes (Vipera ammodytes, Vipera palaestinae, Echis coloratus, Bothrops alternatus, Trimeresurus flavoviridis and Trimeresurus gramineus). The copy number of the 3' end of Bov-B LINE in the Vipera ammodytes genome is between 62,000 and 75,000. The absence of Bov-B LINE at orthologous positions in other snake PLA2 genes indicates that its retrotransposition in the V. ammodytes PLA2 gene locus has occurred quite recently, about 5 My ago. The amplification of Bov-B LINEs in snakes may have occurred before the divergence of the Viperinae and Crotalinae subfamilies. Due to its wide distribution in Viperidae snakes it may be a valuable phylogenetic marker. The neighbor-joining phylogenetic tree shows two clusters of truncated Bov-B LINE, a Bovidae and a snake cluster, indicating an early horizontal transfer of this transposable element.

The TBP gene from Aspergillus nidulans-structure and expression in Saccharomyces cerevisiae

The genomic and cDNA copy of the TATA-binding protein (TBP) gene from the filamentous fungus Aspergillus nidulans have been cloned. The gene is interrupted by four introns, one of which is in the long 5' untranslated region of 615 bp. The transcription initiation site was established and the levels of mRNA were analysed under diverse growth conditions and found to vary severalfold. The gene encodes a protein of 268 amino acids composed of an N-terminal domain of 88 amino acids with no significant homology to other TBPs and a C-terminal domain of 180 amino acids with about 95% homology to other fungal TBPs. A cDNA clone under the yeast ADH1 promoter was able to substitute for the yeast TBP gene in vivo; however, the transformants obtained grew poorly at 35 degrees C and on galactose and glycerol at 30 degrees C, though they could grow in the presence of copper ions or aminotriazole at this temperature. This phenotype may be the result of altered function of A. nidulans TBP in certain yeast transcription activation pathways.

Role for the amino-terminal region of human TBP in U6 snRNA transcription

Basal transcription from the human RNA polymerase III U6 promoter depends on a TATA box that recruits the TATA box-binding protein (TBP) and a proximal sequence element that recruits the small nuclear RNA (snRNA)-activating protein complex (SNAPc). TBP consists of a conserved carboxyl-terminal domain that performs all known functions of the protein and a nonconserved amino-terminal region of unknown function. Here, the amino-terminal region is shown to down-regulate binding of TBP to the U6 TATA box, mediate cooperative binding with SNAPc to the U6 promoter, and enhance U6 transcription.

Accelerated evolution of crotalinae snake venom gland serine proteases

Eight cDNAs encoding serine proteases isolated from Trimeresurus flavoviridis (habu snake) and T. gramineus (green habu snake) venom gland cDNA libraries showed that nonsynonymous nucleotide substitutions have accumulated in the mature protein-coding regions to cause amino acid changes. Southern blot analysis of T. flavoviridis genomic DNAs using two proper probes indicated that venom gland serine protease genes form a multigene family in the genome. These observations suggest that venom gland serine proteases have diversified their amino acid sequences in an accelerating manner. Since a similar feature has been previously discovered in crotalinae snake venom gland phospholipase A2 (PLA2) isozyme genes, accelerated evolution appears to be universal in plural isozyme families of crotalinae snake venom gland.

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.

Genomic structure of the maize TATA-box binding protein 1 (TBP-1): conserved exon/intron structure in eukaryotic TBP genes

The gene system of the TATA-box binding protein (TBP) is well suited for the study of the evolutionary conservation of essential components of eukaryotic transcription initiation. In this context we have isolated and sequenced the maize TBP gene for a comparison with TBP genes from other organisms. In particular, a molecular phylogenetic analysis of the exon/intron structure of these genes including the archaeal TBP homolog (Thermococcus celer) was performed, revealing that the intron insertion probably occurred after the early appearance of the characteristic tandem repeat within the highly conserved C-terminal domain of all known TBPs, but before separation of the eukaryotic progenitor into the different kingdoms.

Phosphorylation of the N-terminal domain of Xenopus TATA-box binding protein by DNA-dependent protein kinase depends on the C-terminal core domain

DNA-dependent protein kinase (DNA-PK) has been shown to phosphorylate several transcription factors in vitro, suggesting that this nuclear enzyme - in addition to its role in DNA repair and recombination - may be involved in transcriptional regulation. In the typical mechanism the DNA-bound kinase phosphorylates a substrate that is bound to the same DNA molecule. Here I report that the Xenopus TATA-box binding protein (xTBP) is hyperphosphorylated by DNA-PK in vitro. The phosphorylation is in the N-terminal domain of the protein but depends fully on the presence of the C-terminal core domain.

Accelerated evolution in the protein-coding regions is universal in crotalinae snake venom gland phospholipase A2 isozyme genes

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.