*Department of Biochemistry, Faculty of Medicine and Bioscience Center, Faculty of Science, National University of Singapore, Singapore;
and
Department of Biochemistry, Faculty of Medicine, University Malaya, Malaysia
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
The venom PLA2s, on the contrary, while possessing a digestive function, exhibit a wide variety of pharmacological properties such as antiplatelet, anticoagulant, hemolytic, neurotoxic (presynaptic), myotoxic, edema-inducing, hemorrhagic, cytolytic, cardiotoxic, and muscarinic inhibitor activities (Rosenberg 1990
; Harris 1991
; Hawgood and Bon 1991
; Yang 1994
). The PLA2s that have been extensively characterized at protein level both structurally and functionally are the extracellular forms, such as the pancreatic PLA2 from mammals (human [Verheij et al. 1983
], bovine [Fleer, Verheij, and De Haas 1978
], rat [O'Hara et al. 1986
], and canine [Kerfelec et al. 1986
]) and the venom phospholipases from a bee (Apis mellifera; Kuchler et al. 1989
) and snakes (viperids [Ritonja and Gubensek 1985
; Wang et al. 1992
], crotalids [Fukagawa et al. 1993
; Faure et al. 1991
], elapids [Tan and Armugam 1989
; Rowan, Harvey, and Menez 1991
]).
Traditionally, the venom secretory PLA2s are classified into three main groups (groups IIII). They are cysteine-rich (five to seven disulphide bonds) proteins with low molecular mass (1318 kDa) which require Ca2+ for catalysis (Heinrikson, Kruger, and Keim 1977
). Many (over 100) such PLA2s have been sequenced and found to be more than 40% identical to each other. The amino acid residues involved in catalysis have always been found to be conserved.
The authors of a phylogenetic analysis using the amino acid sequences of PLA2s (Davidson and Dennis 1990
) proposed that at least two PLA2 genes existed at the time the ancestral lines of reptiles and mammals diverged, long before snakes came into existence. The two major groups of poisonous snakes, the Proteroglypha (Elapidae, including cobras and kraits) and the Solenoglypha (Viperidae, containing rattlesnakes and vipers), each have been considered to have acquired one of the two ancestral PLA2 genes for the production of their venoms. Heinrikson, Kruger, and Keim (1977
) also classified the PLA2s into two groups, based on their Cys residues, as group I PLA2s, having a disulphide bridge between the half cysteines at positions 11 and 69 (found in the venom of Elapidae snakes and in the pancreatic juices of mammals), and group II PLA2s, characterized by a C-terminal extension containing a half cysteine linked to a half cysteine at position 50 (found in Viperidae snake venoms and in mammalian platelets, liver, and spleen).
While there have been a number of studies on the cloning of cDNAs encoding PLA2s (Smith 1990
; Pan, Chang, and Chiou 1994
; Moura-da Silva et al. 1995
; Armugam et al. 1997
) there are only a few reports available on the cloning and characterization of PLA2 genes (Seilhamer et al. 1986
; Kordis and Gubensek 1996
; Ohno et al. 1998
). The structural organizations of the mammalian group I and group II PLA2 genes have been found to be different from each other. The former consists of four exons and three introns (Seilhamer et al. 1986
; Kerfelec et al. 1990
) and the latter of five exons and four introns (Seilhamer et al. 1989
; Komada, Kudo, and Inoue 1990
). The genes of Crotalinae venom PLA2 (Ohno et al. 1998
), however, showed a structure similar to that of group I PLA2 genes, although the protein structure conformed with group II enzymes. The genes of Viperidae snake venom PLA2s (Kordis and Gubensek 1996
) have been found to posses the group II PLA2 gene structure. Kordis and Gubensek later explained that this irregular structure of the Crotalinae group II PLA2 gene is due to a 40-bp deletion in exon 1 of the gene which resulted in abolishment of the splicing of the first intron. This apparent structural irregularity and the absence of information on the group I venom PLA2 genes from elapids prompted us to examine the PLA2 genes in the spitting cobra Naja sputatrix (synonyms: 1827 Naja sputatrix Boie, 1907 Naja naja sputatrix Stejneger, and 1989 Naja sputatrix Wuster and Thorpe; Wuster and Thorpe 1991
). In this report, we present the structure, organization, and promoter analysis of genes encoding two isoforms (acidic and neutral) of PLA2s in the venom of N. sputatrix (Armugam et al. 1997
). Both of the genes contain four exons and three introns and constitute the missing link in the biology of PLA2s and form the first report on the nonpancreatic group I PLA2 from a venomous snake. Besides identifying some of the regulatory elements and proteins that are involved in the control of expression of the elapid PLA2 gene, we also analyzed for the first time the phylogenetic relationship of both group I and group II PLA2s based on the information of whole genes. Group I and group II PLA2s appear to have evolved in separate ways from a common ancestor, whereas the human and elapid group I PLA2s showed a common path of evolution.
![]() |
Materials and Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Amplification of Genomic DNA
Custom-designed oligonucleotide primers, based on the sequence of the highly conserved regions in the PLA2 cDNA of N. sputatrix, were supplied by Oswell DNA Service (United Kingdom). PCR was performed using the Advantage Genomic PCR kit (Clontech). Amplification of PLA2 genes was carried out on 10 ng of genomic DNA of the snake liver in a 50-µl reaction mixture. The thermal profile involved a hot start at 94°C for 1 min, followed by 30 cycles of 30 s at 94°C, 30 s at 55°C, and 3 min at 68°C, followed by a final extension at 68°C for 10 min using a Perkin Elmer Cetus thermal cycler (model 480). The following primers were used in the amplification of the PLA2 gene.
Sense primers were X224F (5'-ATgAATCCTgCTCACCTTCTg-3'), X347F (5'-CTTggTggCATTTTgCggACTACggTTgCTAC-3'), and X223F (5'-TTCAAgACCTATTCATACgAgTgT-3'), and antisense primers were X230R (5'-CgTCCgCAgTAgCAACCgTAgT-3'), X346R (5'-gAATAggTCTTgAAgTAgggCCAgCATCTggA-3'), and X225R (5'-gCCTTgAggTCgATATTgTAgTTg-3'). Primers were designed based on cDNA sequences of N. sputatrix PLA2 (Armugam et al. 1997
) using the PrimerSelect program of DNASTAR Inc. (Madison, Wis.). A fourth set of primers, PGF (sense) (5'-ggCACTgAggATgggATTg-3') and PGR (antisense) (5'-AAggTCCCTgTTgggTCCCTggTgC-3'), were used to amplify the complete PLA2 gene(s).
Amplification of 5'- and 3'-flanking Regions
The 5' and 3' ends of the PLA2 gene were amplified using the Universal GenomeWalker kit (Clontech). Briefly, the procedure involved the construction of the adaptor-ligated libraries (GenomeWalker uncloned libraries) made by separate restriction digestion of genomic DNA with DraI, EcoRV, PvuII, ScaI, and StuI, followed by ligation to a special adaptor provided in the kit. Two sets of primers were used for the amplification of the 5' region: the adaptor primer 1 (AP1) (5'gTAATACgACTCACTATAgggC-3'), provided in the kit, and the gene-specific reverse primer GSP1 (X267) (5'-AATCTCCACTTACCTgCTgCCAggAT-3') for primary PCR and the nested adaptor primer 2 (AP2) (5'-ACTATAgggCACgCgTggT-3') and the nested gene-specific reverse primer GSP2 (X268) (5'-gCTGCCAggATCAgAAggTgAgCAgg-3') for secondary PCR.
Primers used for mapping the 3' end of PLA2 gene(s) were again AP1 and the gene-specific forward primers GSP3 (X276) (5'-gACCgCTTggCAgCCATCTgCTTCg-3') for primary PCR, followed by AP2 and the nested gene-specific forward primer GSP4 (X277) (5'-ATCTgCTTCgCCggAgCCCCTTACAA-3') for secondary PCR. "Touchdown" PCR (Roux 1995
) was conducted for the amplification of the 3' and 5' regions.
Subcloning and DNA Sequencing
The PCR products were analyzed on agarose gel. The appropriate bands were purified (Qiaquick gel extraction kit, Qiagen) and subcloned into pT-Adv vector (Clontech Advantage cloning kit). The ligated products were then transformed into Escherichia coli strain TOP10F', and the recombinants were selected on LB-Amp (50 µg/ml) plates supplemented with IPTG and X-Gal. Putative recombinant clones were then subjected to Sanger dideoxy DNA sequencing (Sanger, Nicklen, and Coulson 1977
) using M13/pUC universal primers as well as gene-specific primers where appropriate. All double-stranded sequencing was performed with the ABI Prism Dye Terminator Cycle Sequencing Ready Reaction kit (Perkin-Elmer) using an Automated DNA sequencer (Applied Biosystems, model 373).
Sequence Analysis
Nucleotide sequence homology searches of nonredundant databases in GenBank (National Centre for Biotechnology Information) were performed using the BLAST program. DNA sequence alignments were carried out using the DNASIS software package from Hitachi Software Engineering. The molecular evolutionary relationships among the phospholipase A2 genes were examined using the MegAlign program from DNASTAR software. Acidic and neutral PLA2s reported in this paper have been assigned GenBank accession numbers AF101235 and AF101236, respectively.
Primer Extension Analysis
Primer extension analysis was performed to identify the transcription initiation site(s) of the PLA2 gene(s). Total cellular RNA was extracted from the N. sputatrix venom glands using the guanidine isothiocyanate method (Chomczynski and Sacchi 1987
), and the integrity of total RNA was analyzed by denaturing formaldehyde agarose electrophoresis (Sambrook, Fritsch, and Maniatis 1989
). The primer extension was carried out according to Lachumanan et al. (1998)
. Ten picomoles of 5'-end 33P-labeled reverse primer GSP1 (5'-AATCTCCACTTACCTgCTgCCAggAT-3') was annealed to RNA (10 µg) at 65°C, and a primer extension reaction was carried out at 42°C for 1 h in the presence of 20 U of MuMLV reverse transcriptase. The reverse-transcribed products were electrophoresed on 6% denaturing polyacrylamide DNA sequencing gel. A template harboring the 5'-flanking region of the PLA2 gene was used as a size marker.
Transfection of Mammalian Cells and Promoter Analysis
Human hepatoma (HepG2) and Chinese hamster ovary (CHO) cells were grown in alpha-MEM media and transfected using the calcium phosphate precipitation method, enhanced by glycerol shock (Promega Corp., Madison, Wis.). Ten micrograms of pMAMneo CAT containing the 5' PLA2 promoter region its deletion constructs were cotransfected separately with 6 µg pSV-ß-Gal for checking the transfection efficiency. Cells were harvested 48 h after transfection and freeze-thawed in 250 mM Tris-HCl (pH 8.0). The CAT enzyme reaction was performed in the presence of 14C-chloramphenicol (0.05 mCi/ml), n-Butyryl CoA in 250 mM Tris-Cl (pH 8.0). The reaction product was separated on silica-coated plates by thin-layer chromatography (TLC), using chloroform/methanol (97:3) as the mobile phase. The TLC plate was then autoradiographed to identify the product. The CAT activity was calculated based on the counts per minute of monoacetylated species counted in the COBRA AutoGamma counter (Packard Instruments Co. Inc.).
DNase I Footprinting
The DNase 1 footprinting analysis was done using the Sure Track footprinting kit (Pharmacia Biotech, Sweden) as described by Garabedian et al. (1993)
. Briefly, probes were prepared by digesting the deletion and the full promoter constructs in pT7 Blue(R) vector (Novagen, Madison, Wis.) with HindIII and SacI restriction enzymes and labeling the HindIII end using Klenow DNA polymerase and [
-32P]dATP and [
-32P]dGTP. The radiolabeled probes were incubated separately with CHO, snake venom gland, snake liver, and HepG2 nuclear extracts for 30 min on ice in binding buffer containing 10 mM Tris-Cl (pH 7.5), 50 mM NaCl, 2.5 mM MgCl2, 0.5 mM DTT, and 4% glycerol. The protein-bound DNA fragment was treated with 0.3 U of DNase 1 for 1 min at room temperature and extracted with phenol/chloroform. The protected DNA was electrophoresed on an 8% polyacrylamide gel under denaturing conditions. The gel was autoradiographed on Kodak X-Omat AR5 film at -80°C.
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Primers X346 and X347, which overlap the reverse primer X230, amplified the regions of exon 2, intron II, and a part of exon 3 of the gene (nucleotides 16322740). The remainder of exon 3, intron III, and exon 4 (~800-bp fragment) were obtained using primers X233 and X225 (nucleotides 27273396). From the nucleotide sequences of the above fragments, carried out on both strands, a complete sequence of the group I PLA2 gene from N. sputatrix was obtained. The gene comprised four exons, exon 1 (34 bp), exon 2 (172 bp), exon 3 (113 bp), and exon 4 (122 bp), and three introns, intron I (1.45 kb), intron II (973 bp), and intron III (508 bp), for a total length of 3.379 kb from the ATG initiation codon to the TAG termination codon (nucleotides 303409).
The 5'-flanking region of the gene was amplified by PCR using an adaptor primer (AP1) and a reverse primer of the signal peptide. Similarly, the 3' end of the gene was obtained by PCR using a forward primer selected from exon 4 and the adaptor primer. PCR products of 395 and 2.08 kb were obtained for the 5' region and the 3' region, respectively, from the EcoRV uncloned library. Amplification of the complete gene was achieved using a primer set at the 5' and 3' regions of the gene. The 4.1-kb fragment obtained was sequenced and confirmed to be the PLA2 gene. Analysis of the gene sequences showed the presence of two PLA2 genes, one of them encoding for an acidic PLA2 (APLA) and the other for a neutral PLA2 (NPLA). Both genes showed identical introns, as well as 5' and 3' regions (GenBank accession numbers AF101235 and AF101236). The variant amino acids (and the exons encoding them) Asp20His20 and Asp39
Gln39 (exon 2); Val46
Ileu46 and Gly52
Asn52 (exon 3), and Asp83
Asn83 (exon 4) for acidic and neutral PLA2s, respectively (Armugam et al. 1997
), can be observed on the deduced amino acid sequences. Except for these regions, both APLA and NPLA genes exhibited identical sequences on all of the introns and exons, as well as the 5' and 3' regions of the gene.
Phylogenetic Analysis of the PLA2 Genes
Nucleotide sequences of the Trimeresurus flavoviridis, Trimeresurus gramineus, Vipera ammodytes, Crotalus scutulatus scutulatus, human, and N. sputatrix PLA2 genes were aligned using the MegAlign DNASTAR software package. The noncoding regions appeared to be highly conserved among the families, while the coding regions showed significant variations. Our phylogenetic analysis (fig. 1 ) using the whole gene sequences showed that the two groups of snakes and mammals originated from a common ancestor and that the Elapidae and Viperidae diverged later from each other.
|
|
|
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
The codons for the conserved residues, important for function, are located in exon 2 of the N. sputatrix gene. The highly conserved substrate recognition amino terminal segment of lipophilic residues Leu2, Phe5, and Ile9 and the overlapping Ca2+-binding loop, containing the glycine-rich sequence Tyr24-Gly25-Cys26-Tyr27-Cys28-Gly29-Arg30-Gly31-Gly32-Ser33-Gly34 (Kramer et al. 1989
), are found in exon 2. Exon 3 contains the coding regions for the characteristic cobra loop (Glu53, Ala54, Glu55) found in cobra PLA2 (Heinrikson, Kruger, and Keim 1977
), as well as the active-site residues His48, Asp49, Tyr52, and Tyr64, important for the formation of the catalytic network at its proximal end. The anticoagulant region, corresponding to residues 5477 (Kini and Evans 1989
), is located at the end of exon 3. The codons for the presynaptic cluster residues 5565 and 8089 (Arriagada and Cid 1989
) are separated by intron III. The codons for both neurotoxicity (residues 80100; Kini and Iwanaga 1986a
) and myotoxicity (residues 8088; Kini and Iwanaga 1986b
) and the region contributing to lethality (Phe100Tyr105; Gubensek et al. 1994
) are located in exon 4. Hence, exons 3 and 4 contain the residues which are important for the pathophysiological activities described for snake venom PLA2.
The nucleotide substitution in the gene was found to be restricted to the protein-coding region (exons 2, 3, and 4). Exon 1 and the introns, as well as the 5' and 3' ends of the gene, have been found to be highly conserved among the two PLA2 genes of N. sputatrix. Exons 2, 3, and 4 exhibited only nonsynonymous nucleotide substitutions leading to amino acid changes in the PLA2 protein. This observation has also been reported for PLA2 genes from Viperidae (Ogawa et al. 1996
). Thus, the elapid group I PLA2 genes could have evolved via positive Darwinian selection in a manner similar to that of Viperidae PLA2 genes.
The total intervening sequences (introns) in T. flavoviridis, V. ammodytes, C. s. scutulatus, N. sputatrix, human, and canine pancreatic PLA2 genes have been found to be 1.28, 1.34, 1.84, 2.93, 5.00, and 8.19 kb, respectively (fig. 4
). The disparity in sizes of the genes among crotalids, elapids, and mammals were found to be mainly due to the differences in sizes of introns 1 and 3. The introns of viperid PLA2 gene are comparatively smaller (163693 bp) than the others. Intron 1 of the N. sputatrix (1.45 kb) PLA2 gene is of approximately the same size as its mammalian counterpart (1.6 kb) and is six times as large as that of crotalids. On the other hand, introns 3 of N. sputatrix (508 bp) and C. s. scutulatus (884 bp) PLA2 genes are significantly smaller (five- and threefold, respectively) than that of the human pancreatic PLA2 gene (fig. 4
). The smaller intron 3 may have provided an added advantage in adaptive evolution of the PLA2 genes in snakes, as illustrated by the additional pharmacological properties possibly encoded by exons 3 and 4 in the venom PLA2 proteins. We speculate that the PLA2 gene acquired longer introns in its sequence during the evolution of mammals from reptiles in order to minimize the rate of mutation and to protect its fundamental role as a catalytic protein. The evolution of multiple copies of cardiotoxin (Lachumanan et al. 1998
) and neurotoxin (Afifiyan et al. 1999
) genes from N. sputatrix containing shorter introns, as well as ion channel toxins encoded by intronless genes found in lower forms of animals such as sea anemone (Gendeh, Chung, and Jeyaseelan 1997
) and scorpions (Delabre et al. 1995
; Becerril, Marangoni, and Possani 1997
), provides support for our hypothesis.
Based on the BLAST analysis, the similarities between human pancreatic and N. sputatrix PLA2 intron I and II were 29% and 31% respectively in the opposite (plus and minus) strands of the DNA, while the Viperidae PLA2 showed 24%37% homology (in the plus strand). ART2 retrotransposons (Kordis and Gubensek 1995
) and CRI-like long terminal repeats (Nobuhisa et al. 1998
) reported to be present in intron 4 and the 3' region of Crotalinae PLA2 genes and postulated to have been acquired by homologous recombination during evolution could not be found in the PLA2 gene of N. sputatrix. Nevertheless, the first 200 bp (minus strand) of intron II of the N. sputatrix PLA2 gene showed 90%100% similarity to intron II of the PLA2 inhibitor gene (AB003473), the PLA2 gene (D31779), and the TATA boxbinding protein gene (D31782, D31777) of Trimeresurus sp. The 3' end of intron II showed 80% homology to the 3' UTR of the serine protease type toxins from the habu snake (D67080), KN-BJ2 from Bothrops jararaca (AB004067), plasminogen activator from T. stejnegeri (U21903), calobin from Agkistrodon ussuriensis (U32937), and salmobin from Gloydius halys (AF056033). This observation shows that while the exons code for the mature protein of PLA2, the intron, particularly intron II, shows some features of "cross-talk" among the introns of related and unrelated genes. The significance of this has yet to be determined. The intron II also shows an 11-bp (GTCTTCTAGTC) direct repeat sequence flanking a 123-bp nucleotide sequence (nucleotides 17671889). Insertion elements such as
-actin-processed pseudogene and endogenous retroviral-like elements (Ting et al. 1992
) have been found to be involved in the tissue-specific expression of the salivary amylase gene. Similarly, insertion of HERV-H family of retroviral elements (Kowalski et al. 1997
) in an intron of the PLA2-like gene (PLA2L) has been shown to alter the expression of the gene in a tissue-specific manner. Although the 123-bp insertion element in intron II of the N. sputatrix PLA2 gene could not be identified, it is possible that it may have a functional role in the expression of the PLA2 gene in the venom gland of this snake.
Ohno et al. (1998)
, using the combined sequences of the 5' and 3' UTRs and the signal peptide region of the Viperidae (group II PLA2), showed that Viperidae which evolved from a common ancestor later gave rise to the Crotalinae and Viperinae groups of PLA2s. Our phylogenetic analysis using the secretory PLA2 whole-gene sequences from humans, Proteroglypha (Elapidae; N. sputatrix), and Solenoglypha (Viperidae; V. ammodytes, T. flavoviridis, and C. s. scutulatus), while confirming the origin of Viperidae genes (Ohno et al. 1998
), provided concrete evidence that the same ancestral gene was responsible for the evolution of both the group I and the group II PLA2s, as postulated by Davidson and Dennis (1990)
. However, Davidson and Dennis (1990)
proposed that at least two PLA2 genes existed at the time the ancestral lines of reptiles and mammals diverged. In our analysis, based on gene sequences, the mammalian and elapid PLA2 genes appear to have diverged from the same ancestral gene. Further evidence for the possible evolution of mammalian and elapid PLA2 can be seen by examining the distribution of TGT/C codons encoding cysteine amino acid residues that are important for the maintenance of the tertiary structure of PLA2. The TGT/C nucleotides are distributed as 3, 6, and 5 codons among exons 2, 3, and 4, respectively, in both N. sputatrix and mammalian PLA2 genes, unlike the 2, 6, and 6 distribution of these codons in the corresponding exons of viperids and crotalids.
The presence of an extra loop is characteristic of the pancreatic PLA2 enzymes in mammals. This loop, referred to as the "pancreatic loop," has been thought to be present only in mammalian PLA2. In recent years, there have been reports on the presence of a pancreatic loop in venom PLA2s of the elapids such as the king cobra (Huang et al. 1996
), Brazilian coral snakes (Francis et al. 1997
), and Australian snakes (Jeyaseelan et al. 1998
). The pancreatic loop in mammalian group I PLA2s is encoded by exon 3 of the gene. However, this loop is absent in the N. sputatrix PLA2s. The loss of the pancreatic loop was suggested as an added adaptive advantage for the development of the toxic properties among venom PLA2s (Davidson and Dennis 1990
). Hence, it is possible that snake PLA2, especially that of N. sputatrix, via a positive Darwinian type of accelerated evolution (Ogawa et al. 1996
), could have lost its pancreatic loop from its exon 3 in order to acquire additional toxic properties, while the mammalian PLA2s which have been considered to have evolved under neutrality retained the pancreatic loop to minimize the rate of mutation of the gene. It is interesting to note that the PLA2s of Australian snakes and the king cobra, when compared with those of the modern elapids N. sputatrix and N. atra (Jeyaseelan et al. 1998
), appear at the lower end of the evolutionary tree, indicating that the presence of the pancreatic loop in PLA2 is an ancestral property. Thus, it appears that exons 3 and 4, which are separated by a short intron III in N. sputatrix, form the hot spots for mutations leading to diverse pathophysiological properties while retaining the other fundamental features of the PLA2 enzymes, such as the catalytic property of the exon 2 of the gene.
Promoter activity studies of the CHO and HepG2 cell lines showed that the expression of the PLA2 gene is tissue-specific, as the construct 5'PLA-233 exhibited the highest promoter activity in CHO cells and the lowest activity in HepG2 cells. This deletion construct contains, apart from the TATA and the TIS, four Sp1, two -IRE, one NF-IL6, and one AP-2 cis-acting element. The promoter region between nucleotides -116 and -233 appears to contain crucial cis-acting elements, which are involved in the up/down regulation of the gene expression. This region appears to recognize and differentiate between a normal cell and a hepatoma cell. DNase 1 footprinting analysis using the nuclear extracts from snake venom gland, snake liver, CHO, and HepG2 cells showed that an 87-bp region (nucleotides -146 to -233) containing one proximal Sp1 (positions -153 and -171) and two
-IRE sites (positions -165 and -187) has been protected (fig. 3A
). Footprinting using CHO cell nuclear proteins shows protection at nucleotides -220 to -233, containing CF1 elements, as well as nucleotides -147 to -206 (two
-IRE sites and an Sp1 site) and an 8-bp region upstream of the CF1 site. However, only the -147 to -206 region is protected by the HepG2 nuclear protein. It appears that both the snake venom gland and snake liver nuclear proteins show interactions similar to that of CHO cell nuclear extract. Whether the proteins binding to these regions are identical or not remains to be elucidated. The binding of HepG2 nuclear proteins to this region (Sp1 and flanking
-IREs) appears to repress the promoter activity, as demonstrated by the CAT assay of 5'PLA-233 (fig. 2C
). The
-IRE elements which are known to activate the human cytoplasmic PLA2 promoter activity (Wu et al. 1994
) were shown to inhibit the transcription activity of the venom secretory PLA2. Thus, the Sp1 site and its two flanking
-IRE elements at the upstream region of AP2 (position -113) behave as silencer elements in HepG2 cells. While only a slight increase has been observed by the removal of IRE elements (5'PLA-266 to 5'PLA-233), the removal of the 2Sp1 sites (5'PLA-266 to 5'PLA-116) showed a significant increase in CAT activity in HepG2 cells. A similar observation has been reported by Hayashi et al. (1998)
for
-IRE and Sp1 elements on PCI gene expression. Our promoter analysis using CHO cells demonstrates that these
-IRE and Sp1 elements, in combination with the CF1 element and an unknown factor (nucleotides -220 to -233) have an effect (enhancement) opposite that observed in HepG2 cells (fig. 2C
), thus confirming the tissue specificity of the PLA2 promoter. Hayashi et al. (1998)
observed that the AP-2 element in the 5'-flanking region is capable of exhibiting an enhancer/repressor activity depending on its position with respect to Sp1 sites. In our study, an AP2 element situated at position -113 on the promoter plays a suppressive role in HepG2 cells (fig. 2C
). A basal promoter activity can be observed with TATA and the TIS (5'PLA-51 clone). Inclusion of two Sp1 sites at positions -65 and -68 in 5'PLA-90 yielded an increase in promoter activity in HepG2 cells and a decrease in activity in CHO cells. Similarly, Park, Plummer, and Krystal (1998)
reported that selective binding of Sp1 at the promoter region is critical for maximal activity of the human c-kit promoter. Addition of an AP2 element to the 5'PLA-90 construct (5'PLA-116) lowered the level of activity by 30% in a tissue-specific manner in HepG2 cells. Thus, the AP2 site, the two SpI sites (at positions -153 and -171) and the two
-IRE elements (at positions -165 and -187) function as suppressors in HepG2 cells and possibly as enhancers in CHO cells. The downstream Sp1 sites (at positions -65 and -68) have a suppressive role in HepG2 and an enhancing role in CHO cells. Thus, it appears that the Sp1 elements could function in combination with the nearby enhancer/repressor-binding proteins as promoter-specific factors to regulate gene expression in a tissue-specific manner, as in the case of the snake PLA2 gene. The functions of the third
-IRE element (at position -246) and the NF-IL6 elements (at positions -246 and -134) are still unknown.
Alternating purine-pyrimidine dinucleotide repeats are known to form Z-DNA conformations in vitro; hence, it has been speculated that their existence in the 5'-flanking regions of genes may confer a regulatory effect on gene transcription (Brahmachari et al. 1995
). The CA-repeats in the 5'-flanking region have been reported to cause an inhibitory effect on the gene expression of the rat prolactin gene (Naylor and Clark 1990
) and the human cytoplasmic PLA2 gene (Hayashi et al. 1998
). Singer and Gottschling (1994)
also reported that telomeres (TG repeats) binding to telomerase while maintaining genomic integrity could function as silencer elements. Seeding of telomeres has been shown to require a telomeric sequence motif, GGGTTA, in vertebrates (Konig, Fairall, and Rhodes 1998
). The telomeric repeatbinding proteins are also known to contain a DNA-binding motif similar to that of the Myb-like domain. Interestingly, a (TG12) repeat, a cMyb-binding site, and a GGGTTA site can be observed in the 5' region of the PLA2 gene of N. sputatrix (fig. 2A
). Removal of this (TG)12 repeat from the promoter region brought about an increase in promoter activity in both CHO and HepG2 cells (fig. 2C
), thus confirming its role as a silencer in gene regulation. A second (TG)12 repeat can be found at the 3' end of the N. sputatrix PLA2 gene. Its function is still unknown. It is possible that this (TG)12 repeat could also be involved in the regulation of the gene activity.
The 3' region of N. sputatrix PLA2 contains two functional poly-A signal sites. Studies of 3' rapid amplification of cDNA ends (results not shown) showed that the N. sputatrix PLA2 mRNA could be terminated at either one of these two sites, which are about 50 bp apart from each other. An AT-rich region known to destabilize mRNA transcripts (Higgins 1991
) is located in the 3' region of the gene. This AT-rich region, together with the (TG)12 repeats found in both the 5' and the 3' ends of the gene, could possibly explain the lower turnover of PLA2 message in the venom gland cells in comparison with the cardiotoxin messages that we observed previously in a separate study (Lachumanan et al. 1999
). Thus, the regulation of PLA2 gene expression (in a tissue-specific manner) appears to be controlled by concerted activities of the 5' regulatory elements, as well as the 3' AT- and TG-rich regions.
![]() |
Conclusions |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Footnotes |
---|
1 Keywords: phospholipase A2 gene
group I PLA
Naja sputatrix,
spitting cobra
phylogeny
adaptive evolution
2 Address for correspondence and reprints: Kandiah Jeyaseelan, Department of Biochemistry, Faculty of Medicine, National University of Singapore, 10 Medical Drive, Singapore 119260. E-mail: bchjeya{at}nus.edu.sg
![]() |
literature cited |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Afifiyan, F., A. Armugam, C. H. Tan, P. Gopalakrishnakone, and K. Jeyaseelan. 1999. Postsynaptic alpha-neurotoxin gene of the spitting cobra, Naja naja sputatrix: structure, organization, and phylogenetic analysis. Genome Res. 9:259266.
Armugam, A., L. Earnest, M. C. M. Chung, P. Gopalakrishnakone, N. H. Tan, C. H. Tan, and K. Jeyaseelan. 1997. Cloning and characterization of cDNAs encoding three isoforms of phospholipase A2 in Malayan spitting cobra (Naja naja sputatrix) venom. Toxicon 35:2737.
Arriagada, E., and H. Cid. 1989. Search for a toxic-site in snake venom phospholipase A2. Arch. Biol. Med. Exp. 22:97105.[ISI][Medline]
Becerril, B., S. Marangoni, and L. D. Possani. 1997. Toxins and genes isolated from scorpions of the genus Tityus. Toxicon 35:821835.
Brahmachari, S. K., G. Meera, P. S. Sarkar, P. Balagurumoorthy, J. Tripathi, S. Raghavan, U. Shaligram, and S. Pataskar. 1995. Simple repetitive sequences in the genome: structure and functional significance. Electrophoresis 16:17051714.
Breathnach, R., and P. Chambon. 1981. Organization and expression of eukaryotic split genes coding for proteins. Annu. Rev. Biochem. 50:349383.[ISI][Medline]
Chomczynski, P., and N. Sacchi. 1987. Single step method of RNA isolation by acid-guanidinium-thiocyanate-phenol-chloroform extraction. Ann. Biochem. 162:156159.
Davidson, F. F., and E. A. Dennis. 1990. Evolutionary relationships and implications for the regulation of phospholipase A2 from snake venom to human secreted forms. J. Mol. Evol. 31:228238.[ISI][Medline]
Delabre, M. L., P. Pasero, M. Marilley, and P. E. Bougis. 1995. Promoter structure and intron-exon organization of a scorpion alpha-toxin gene. Biochemistry 34:67266736.
Dennis, E. A. 1983. Phospholipases. Pp. 307353 in P. D. Boyer, ed. The enzymes. Vol. 16, 3rd edition. Academic Press, New York.
. 1994. Diversity of group types, regulation and function of phospholipase A2. J. Biol. Chem. 269:1305713060.
Dennis, E. A., R. A. Deems, and L. Yu. 1992. Extracellular phospholipase A2. Adv. Exp. Med. Biol. 318:3539.[Medline]
Dennis, E. A., S. G. Rhee, M. M. Billah, and Y. A. Hannum. 1991. Role of phospholipase in generating second messengers in signal transduction. FASEB J. 5:20682077.
Faure, G., J. L. Guillame, L. Camoin, B. Saliou, and C. Bon. 1991. Multiplicity of acidic subunit isoforms of crotoxin, the phospholipase A2 neurotoxin from Crotalus durissus terrificus venom, results from post-translational modifications. Biochemistry 30:8074.
Fleer, E. A. M., H. M. Verheij, and G. H. De Haas. 1978. The primary structure of bovine pancreatic phospholipase A2. Eur. J. Biochem. 82:261269.[Abstract]
Francis, B. R., N. J. Da Silva Junior, C. Seebart, L. L. Casais De Silva, J. J. Schmidt, and I. I. Kaiser. 1997. Toxins isolated from the venom of Brazilian coral snake (Micrurus frontalis frontalis) include hemorrhagic type phospholipase A2 and postsynaptic neurotoxins. Toxicon 35:193203.
Fukagawa, T., T. Nose, Y. Shimohigashi, T. Ogawa, N. Oda, K. I. Nakashima, C. C. Chang, and M. Ohno. 1993. Purification, sequencing and characterization of single amino acid-substituted phospholipase A2 isozymes from Trimeresurus gramineus (green habu snake) venom. Toxicon 31:957.
Garabedian, M. J., J. La Baer, W. H. Liu, and J. R. Thomas. 1993. Analysis of protein-DNA interactions. Pp. 243293 in B. D. Hames and S. J. Higgins, eds. Gene transcription: a practical approach. IRL, Oxford, England.
Gendeh, G. S., M. C. M. Chung, and K. Jeyaseelan. 1997. Genomic structure of a potassium channel toxin from Heteractis magnifica. FEBS Lett. 418:183188.
Gubensek, F., N. S. Liang, J. Pungercar, B. Strukelj, V. Curin-Serbec, and I. Krizaj. 1994. Presynaptically acting phospholipase A2 from Vipera ammodytes venom. Ann. N.Y. Acad. Sci. 710:120125.[ISI][Medline]
Harris, J. B. 1991. Phospholipases in snake venom and their effects on nerve and muscle. Pp. 91129 in A. L. Harvey, ed. Snake toxins. Pergamon Press, London.
Hawgood, B. J., and C. Bon. 1991. Snake venom presynaptic toxins. P. 3 in A. T. Tu, ed. Handbook of natural toxins, Vol. 5. Reptile venoms and toxins. Marcell Dekker, New York.
Hayashi, T., M. Usui, J. Nishioka, Z. X. Zhang, and K. Suzuki. 1998. Regulation of the human protein C inhibitor gene expression in HepG2 cells. Role of Sp1 and AP2. Biochem. J. 332:573582.
Heinrikson, R. L., E. Kruger, and P. S. Keim. 1977. Amino acid sequence of phospholipase A2- from the venom of Crotalus adamenteus. A new classification of phospholipase A2 based upon structural determinants. J. Biol. Chem. 252:4913.
Higgins, C. F. 1991. Stability and degradation of mRNA. Curr. Opin. Cell Biol. 3:10131018.[Medline]
Huang, M. Z., P. Gopalakrishnakone, M. C. M. Chung, and R. M. Kini. 1996. Complete amino acid sequence of an acidic, cardiotoxic phospholipase A2 from the venom of Ophiophagus hannah (King cobra): a novel cobra venom enzyme with pancreatic loop. Arch. Biochem. Biophys. 332:150156.
Jeyaseelan, K., A. Armugam, R. P. Nair et al. (14 co-authors). 1998. Snake toxin genes: cloning, structure-function analysis and phylogeny. From venoms to drugs, Inaugral Conference, Aug 1998, Queensland, Australia. Toxicon 37:1484 [Abstract].
John, T. R., J. J. Smith, and I. I. Kaiser. 1994. Genomic sequences encoding the acidic and basic subunits of Mojave toxin: unusually high sequence identity of non-coding region. Gene 139:229234.
. 1996. A phospholipase A2-like pseudogene retaining the highly conserved introns of Mojave toxin and other snake venom group II PLA2s, but having different exons. DNA Cell Biol. 15:661668.[ISI][Medline]
Kerfelec, R., K. S. La Forge, A. Puigserver, and G. Scheele. 1986. Primary structures of canine pancreatic lipase and phospholipase A2 messenger RNAs. Pancreas 1:430437.
Kerfelec, B., K. S. Laforge, P. Vasiloudes, A. Puigserver, and G. A. Scheele. 1990. Isolation and sequence of the canine pancreatic phospholipase A2 gene. Eur. J. Biochem. 190:299304.[Abstract]
Kini, R. M., and H. J. Evans. 1989. Structure-function relationships of phospholipase A2. The anticoagulant region of phospholipase A2. J. Biol. Chem. 262:1440214407.
Kini, R. M., and S. Iwanaga. 1986a. Structure-function relationships of phospholipases A2I: prediction of presynaptic neurotoxicity. Toxicon 24:527541.
. 1986b. Structure-function relationships of phospholipases A2II: charge density distribution and the myotoxicity of presynaptically neurotoxic phospholipase A2. Toxicon 24:895905.
Komada, M., I. Kudo, and K. Inoue. 1990. Structure of gene coding for rat group II phospholipase A2. Biochem. Biophys. Res. Commun. 168:10591065.[ISI][Medline]
Konig, P., L. Fairall, and D. Rhodes. 1998. Sequence specific DNA recognition by the Myb-like domain of the human telomere binding protein TRF-1: a model for the protein DNA complex. Nucleic Acids Res. 26:17311740.
Kordis, D., and F. Gubensek. 1995. Horizontal SINE transfer between vertebrate classes. Nat. Genet. 10:131132.[ISI][Medline]
. 1996. Ammodytoxin C gene helps to elucidate the irregular structure of Crotalinae group II phospholipase A2 genes. Eur. J. Biochem. 240:8390.[Abstract]
Kowalski, P., J. D. Freeman, D. T. Nelson, and D. L. Mager. 1997. Genomic structure and evolution of a novel gene (PLA2L) with duplicated PLA2-like domains. Genomics 39:3846.
Kramer, M., C. Hession, B. Johansen, G. Hayes, P. Mcgray, E. P. Chow, R. Tizard, and B. Pepinsky. 1989. Structure and properties of a human non-pancreatic phospholipase A2. J. Biol. Chem. 264:57685775.
Kuchler, K., M. Gmachl, M. J. Sippl, and G. Kreil. 1989. Analysis of the cDNA for phospholipase A2 from honeybee venom glands. Eur. J. Biochem. 184:249254.[Abstract]
Lachumanan, R., A. Armugam, D. Ponraj, P. Gopalakrishnakone, C. H. Tan, and K. Jeyaseelan. 1999. In situ hybridization and immunohistochemical analysis of the expression of cardiotoxin and neurotoxin genes in Naja naja sputatrix. J. Histochem. Cytochem. 47:551560.
Lachumanan, R., A. Armugam, C. H. Tan, and K. Jeyaseelan. 1998. Structure and organization of the cardiotoxin genes in Naja naja sputatrix. FEBS Lett. 433:119124.
Moura-da Silva, A. M., M. J. I. Paine, M. R. V. Diniz, R. D. G. Theakston, and J. M. Crampton. 1995. The molecular cloning of a phospholipase A2 from Bothrops jararacussu venom: evolution of venom group II phospholipase A2 may imply gene duplications. J. Mol. Evol. 41:174179.[ISI][Medline]
Mukherjee, A. B., L. Miele, and N. Pattabiraman. 1994. Phospholipases A2 enzymes: regulation and physiological role. Biochem. Pharmacol. 48:110.[ISI][Medline]
Nakashima, K. I., I. Nobushisa, M. Deshimaru et al. (11 co- authors). 1995. Accelerated evolution in the protein coding regions is universal in crotalinae snake venom gland phospholipase A2 isozyme genes. Proc. Natl. Acad. Sci. USA 92:56055609.
Nakashima, K. I., T. Ogawa, N. Oda, M. Hattori, Y. Sasaki, H. Kihara, and M. Ohno. 1993. Accelerated evolution of Trimeresurus flavoviridis venom gland phospholipase A2 isozymes. Proc. Natl. Acad. Sci. USA 90:59645968.
Naylor, L. H., and E. M. Clark. 1990. d(TG)n.d(CA)n sequences upstream of the rat prolactin gene from Z-DNA and inhibit gene transcription. Nucleic Acids Res. 18:15951601.[Abstract]
Nobuhisa, I., T. Ogawa, M. Deshimaru, T. Chijiwa, K. I. Nakashima, Y. Chuman, Y. Shimohigashi, Y. Fukumaki, S. Hattori, and M. Ohno. 1998. Retrotransposable CR1-like elements in Crotalinae snake genomes. Toxicon 36:915920.
Ogawa, T., K. Nakashima, I. Nobuhisa, M. Deshimaru, Y. Shimohigashi, Y. Fukumaki, Y. Sakaki, S. Hattori, and M. Ohno. 1996. Accelerated evolution of snake venom phospholipase A2 isozymes for acquisition of diverse physiological functions. Toxicon 34:12291236.
O'Hara, O., M. Tamiki, E. Nakamura, Y. Tsuruta, Y. Fujii, M. Shin, H. Teraoka, and M. Okamoto. 1986. Dog and rat pancreatic phospholipase A2: complete amino acid sequences deduced from complementary DNAs. J. Biochem. 99:733739.[Abstract]
Ohno, M., R. Menez, T. Ogawa et al. (12 co-authors). 1998. Molecular evolution of snake toxins: is the functional diversity of snake toxins associated with a mechanism of accelerated evolution? Prog. Nucleic Acids Res. Mol. Biol. 59:307364.
Pan, F. M., W. C. Chang, and S. H. Chiou. 1994. cDNA and protein sequences coding for the precursor of phospholipase A2 from Taiwan cobra Naja naja atra. Biochem. Mol. Biol. Int. 33:187194.
Park, G. H., H. K. Plummer, and G. W. Krystal. 1998. Selective binding is critical for maximal activity of the Human c-kit promoter. Blood 92:41384149.
Ritonja, A., and F. Gubensek. 1985. Ammodytoxin A, a highly lethal phospholipase A2 from Vipera ammodytes ammodytes venom. Biochim. Biophys. Acta 828:306.
Rosenberg, P. 1990. Phospholipases. P. 67 in W. T. Shier and D. Mebs, eds. Handbook of toxinology. Marcel Dekker, New York.
Roux, K. H. 1995. Optimization and troubleshooting in PCR. PCR Methods Appl. 4:51855194.
Rowan, E. G., A. L. Harvey, and A. Menez. 1991. Neuromuscular effects of nigexine, a basic phospholipase A2 from Naja nigricollis venom. Toxicon 29:371374.
Sambrook, J., E. F. Fritsch, and T. Maniatis. 1989. Molecular cloning: a laboratory manual. 2nd edition. Cold Spring Harbor Laboratory Press, New York.
Sanger, F., S. Nicklen, and A. R. Coulson. 1977. DNA sequencing with chain terminating inhibitors. Proc. Natl. Acad. Sci. USA 74:54635467.
Scheele, G., and R. Jacoby. 1983. Proteolytic processing of presecretory proteins is required for development of biological activities in pancreatic exocrine proteins. J. Biol. Chem. 258:20052009.
Seilhamer, J., W. Pruzanski, P. Vadas, S. Plant, J. A. Miller, J. Kloss, and L. K. Johnson. 1989. Cloning and recombinant expression of phospholipase A2 present in rheumatoid arthritic synovial fluid. J. Biol. Chem. 264:53355338.
Seilhamer, J. J., T. L. Randall, M. Yamanaka, and L. K. Johnson. 1986. Pancreatic phospholipase A2: isolation of the human gene and cDNAs from porcine pancreas and human lung. DNA 5:519527.
Singer, M. S., and D. E. Gottschling. 1994. TLC1, the template RNA component of the Saccharomyces cerevisiae telomerase. Science 266:404409.
Smith, L. A. 1990. Cloning, characterization and expression of animal toxin genes for vaccine development. J. Toxicol. Toxins Rev. 9:243283.
Tan, N. H., and A. Armugam. 1989. The anticoagulant activity of Malayan cobra (Naja naja sputatrix) venom and venom phospholipase A2 enzyme. Biochem. Int. 19:803810.[ISI][Medline]
Ting, C. N., M. P. Rosenberg, C. M. Snow, L. L. Samuelson, and M. H. Meisler. 1992. Endogenous retroviral sequences are required for tissue specific expression of a human salivary amylase gene. Genes Dev. 6:14571465.[Abstract]
Verheij, H. M., J. Westerman, B. Sternby, and G. De Haas. 1983. The complete primary sequence of phospholipase A2 from human pancreas. Biochim. Biophys. Acta 747:9399.
Wang, Y. M., P. J. Lu, C. L. Ho, and I. H. Tsai. 1992. Characterization and molecular cloning of neurotoxic phospholipase A2 from Taiwan viper (Vipera russelli formonsensis). Eur. J. Biochem. 209:635641.[Abstract]
Wu, T., T. Ikuzono, C. W. Angus, and J. H. Sheilhamer. 1994. Characterization of the promoter for the human 85kDa cytosolic phospholipase A2 gene. Nucleic Acids Res. 22:50935098.[Abstract]
Wuster, W., and R. S. Thorpe. 1991. Asiatic cobras: systematics and snakebite. Experientia 47:205209.
Yang, C. C. 1994. Structure-function relationship of phospholipase A2 from snake venoms. J. Toxicol. Toxin Rev. 13:125177.[ISI]