Polymorphisms in the osteopontin promoter affect its transcriptional activity
Francesca Giacopelli1,
Renato Marciano1,
Angela Pistorio2,
Paolo Catarsi3,
Silvana Canini4,
Gerard Karsenty5 and
Roberto Ravazzolo1,6
1 Laboratory of Molecular Genetics G. Gaslini Institute, Genova, Italy
2 Epidemiology and Statistics Unit, G. Gaslini Institute, Genova, Italy
3 Laboratory of Nephrology, G. Gaslini Institute, Genova, Italy
4 Laboratory of Chemical-Clinical Analysis, G. Gaslini Institute, Genova, Italy
5 Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas
6 Department of Pediatrics and Center of Excellence for Biomedical Research, University of Genova, Genova, Italy
 |
ABSTRACT
|
---|
Understanding the molecular mechanisms that underlie regulation of transcription of the human osteopontin encoding gene (OPN) may help to clarify several processes, such as fibrotic evolution of organ damage, tumorigenesis and metastasis, and immune response, in which OPN overexpression is observed. With the aim to evaluate variants with functional effect on transcription, we have analyzed the promoter region and focused our investigation on three common variants present in the first 500 bp upstream of the transcription start site. Transfection of constructs carrying the four most frequent haplotypes relative to variants at 66, 156, and 443 fused to the luciferase reporter gene in a panel of different cell lines showed that one haplotype conferred a significantly reduced level of reporter gene expression in all tested cell lines. We describe that the 66 polymorphism modifies the binding affinity for the SP1/SP3 transcription factors, the 156 polymorphism is included in a yet uncharacterized RUNX2 binding site, and the 443 polymorphism causes differential binding of an unknown factor. The finding of differential effects of various combination of variants in haplotypes may contribute to explain data of association studies reported in several already published articles. Future association studies using haplotypes instead of single OPN variants will allow to achieve more accurate results referable to differential expression of OPN in several common diseases, in which OPN is considered a candidate susceptibility gene.
haplotype; gene regulation; Sp1; RUNX2
 |
INTRODUCTION
|
---|
OSTEOPONTIN (OPN) is a phosphorylated acidic glycoprotein that displays several functions in different physiological and pathological processes, including bone remodeling, cell-mediated immunity, maintenance or reconfiguration of tissue integrity during inflammatory processes, coronary restenosis, and tumor cell metastasis. The protein contains an Arg-Gly-Asp (RGD) motif, is found both as an immobilized extracellular matrix molecule in mineralized tissues and as a cytokine in body fluids, is able to interact with various receptors, including
v-ß3 and other integrins, and it may also be a ligand for certain variant forms of CD44. These interactions induce cellular signaling pathways, allowing OPN to mediate cell-matrix and cell-cell interactions (6, 40).
Overexpression of OPN has been described in several conditions in which basic inflammatory processes are activated, such as arthritis (24, 27, 43), myocardial remodeling after infarction (37), kidney interstitial fibrosis after obstructive uropathy and other renal insults (42), and wound healing (21, 29).
The OPN encoding gene [secreted phosphoprotein 1 (SPP1), MIM #166490] is mapped on human chromosome 4q21-q25, together with other members of the so-called SIBLING family of proteins, bone sialoprotein and dentin matrix protein-1, which share some structural characteristics (10). The OPN gene gives rise to different mRNA transcripts because of alternative spliced isoforms (30), and, moreover, numerous OPN isoforms are also due to posttranslational processing, such as phosphorylation and glycosylation and proteolytic cleavage by thrombin and matrix metalloproteinases (1, 33, 41).
OPN knockout mice have provided several pieces of information about OPN function. The first report of targeted inactivation of OPN described impaired wound healing (21). Several other articles followed, describing reduced macrophage infiltration and interstitial fibrosis in the kidney in mouse models of renal fibrosis (26, 29); reduced loss of bone mineral in ovariectomized null mice (45); resistance to progression in an experimental mouse model of multiple sclerosis (4, 18); impaired type-1 immunity to viral and bacterial infection and granulomatous response (3); abnormal tissue remodeling after myocardial infarction (38); and increased arterial calcification (34). On the other hand, in OPN-overexpressing transgenic mice, OPN plays an important role in the development of vascular medial thickening without injury and in neointima formation after arterial injury in vivo (16).
Genetic variations in the OPN gene have been described, and associations of these polymorphisms with lupus erythematosus, multiple sclerosis, urolithiasis, primary biliary cirrhosis, and autoimmune lymphoproliferative syndrome (ALPS) have been reported (5, 11, 19, 23, 44). In previous work we have described a polymorphic variant in the first intron of human OPN (12) that, although devoid of functional significance, might be a marker for another or other linked functional sites.
Here we present results about the functional characterization of other polymorphisms in a 474-bp region upstream of the OPN transcription start site, estimation of linkage disequilibrium (LD) between all polymorphic sites, and haplotype frequencies.
 |
MATERIALS AND METHODS
|
---|
Analysis of polymorphisms in OPN regulatory regions.
Genomic DNA was extracted by standard procedure from peripheral blood of unrelated individuals of the same geographical origin (Italian), collected after informed consent.
The entire regulatory region from 2236 to +1200 was amplified by PCR with paired primers and sequenced with both forward and reverse primers; oligonucleotides used for PCR amplification and sequencing and PCR annealing temperature are reported in Table 1. PCR was performed using 50 ng DNA as a template under the following conditions: 94°C for 7 min, then 30 cycles of 94°C for 55 s, annealing temperature for 55 s and 72°C for 90 s, final extension at 72°C for 10 min. After affinity membrane purification, the PCR products were subjected to cycle sequencing with the respective forward and reverse primers. Sequence analysis was performed using an automated ABI PRISM 3100 DNA sequencer.
The intronic polymorphism at +245 (TG/TGTG) was analyzed as previously described (12).
Haplotype reconstruction and statistical analysis.
Pairwise LD was measured by D' using the Arlequin software ver. 2.000 (32), by which the Hardy-Weinberg equilibrium was also tested. Haplotypes definition and frequency estimation were obtained by using the computer program PHASE 2.02 (35, 36). Estimates were obtained using 1,000 iterations, 10 thinning intervals, and 1,000 burn-in cycles.
Descriptive statistics were reported either in text or in figures in terms of medians and first and third quartiles for quantitative data. As the number of observations was quite low, all nonparametric tests were used to perform statistical analysis. Mann-Whitney U test was used to compare parameters in the two groups of data (for example: luciferase activity B haplotype transfected with Runx2 vs. control vector lacking the Runx2 cDNA). The one-way nonparametric analysis of variance (Kruskal-Wallis test) was used to compare parameters in more than two groups of data, and the Newman-Keuls test was applied for multiple comparisons to explore post-hoc differences between pairs of experiments; whenever the sample size of the experiments was not the same, the Dunn test was performed. All statistical tests were two-sided; a P value of less than 0.05 was considered statistically significant; the statistical package "Statistica for Windows" (release 6, StatSoft) was used to perform all the analyses.
Plasmid construction.
Fragments from 1206 to +42 were prepared by PCR amplification from genomic DNA (sense primer 1206 OPNF and antisense primer 2304 R as in Ref. 12) of homozygous individuals for the four most common haplotypes and cloned in pCR II Topo TA vector (Invitrogen). These plasmids were then digested with HindIII, which cuts at 471 in the SSP1 promoter and in multiple cloning site of pCR II Topo vector. Fragments obtained were cloned in the pGL3 Basic Vector (Promega), giving rise to A, B, C, and D plasmids which contain fragments of different alleles of the human OPN gene from 474 to +42. Site-specific mutagenesis was performed according to a PCR-based procedure (13) to obtain E and F plasmid constructs, using C and D plasmids as templates.
The pCMV-Osf2 plasmid, containing the full-length mouse Runx2 cDNA (9), and pCMV-SP1 plasmid, containing the SP1 cDNA (46) under the control of the cytomegalovirus (CMV) promoter, were kindly provided by Drs. P. Ducy and F. Ramirez, respectively.
All the plasmids were carefully verified by restriction analysis and complete DNA sequencing.
Cell culture and transfection.
All cell lines were maintained in Dulbeccos modified Eagles medium (DMEM) except for HeLa, which was maintained in minimum essential medium, supplemented with 10% (vol/vol) fetal bovine serum (GIBCO-BRL) , 2 mM L-glutamine, 100 U/ml penicillin, and 100 µg/ml streptomycin at 37°C in a humidified atmosphere with 5% CO2. Nonessential amino acids (Euroclone) were included for HEK293 and HeLa cells. Human proximal tubular epithelial cells (HK2) were cultured in DMEM/F-12 containing 5 µg/ml insulin, 5 µg/ml transferrin, 5 ng/ml selenium, 5 ng/ml hydrocortisone, 5 pg/ml 3,5,3'-triiodothyronine, 5 pg/ml prostaglandin E1, 10 ng/ml epidermal growth factor, 5% (vol/vol) fetal bovine serum, 2 mM L-glutamine, and 100 U/ml penicillin and streptomycin.
All cell lines were available in the laboratory and originally obtained from the American Type Culture Collection.
Transfections were performed using the polyethylenimine (PEI) cationic polymer, as previously described (25), except for HeLa cells which were transfected using Lipofectamine 2000 reagent (Invitrogen) according to manufacturers instructions.
Cells were plated into six-well dishes at 5060% confluence and transfected at an estimated 8090% confluence after 24 h. The amount of DNA used for transfection was 2.1 µg for each haplotype plasmid, 3.0 µg for each expression vector, and 0.1 µg for Renilla luciferase control construct pRL-CMV (Promega).
Cell lysates were prepared with the Luciferase Assay System kit (Promega), and luciferase activity was determined with a luminometer (Turner BioSystems).
The data obtained represent median values within first and third quartile of at least three independent experiments, each performed either in triplicate or in duplicate. Luciferase activity was calculated as normalized activity by comparison to a cotransfected Renilla luciferase control construct pRL-CMV (Promega). All data were also confirmed using at least two different plasmid preparations.
Preparation of nuclear extracts and electrophoretic mobility shift assays.
Nuclear extracts were prepared as described previously (7) from U2OS, COS7, and HeLa cell lines, and protein content was determined using a Bradford protein assay kit (Bio-Rad). The probes used for electrophoretic mobility shift assay (EMSA) are listed in Table 1 (upper strand, nucleotide variant in bold). The probes were labeled with [
-32P]ATP. Labeled oligonucleotide probe (15 fmol) were incubated with 510 µg of nuclear extract for 20 min at room temperature in the presence of 1x gel shift binding buffer containing 10 mM HEPES, pH 7.9, 0.1 mM EDTA, 50 mM KCl, 0.5 mM dithiothreitol, 10% glycerol, 1 µg of poly(dI-dC)/poly(dI-dC), and 1x protease inhibitor cocktail (Roche Molecular Biochemicals). Competition in EMSA was performed with either 50-, 150-, or 250-fold molar excess of unlabeled oligonucleotide. For the supershift assay, Sp1, Sp3 (Santa Cruz Biotech), and Runx2 antibodies (8) were included in the incubation mixture for 20 min in ice. Complexes were separated in a 5% nondenaturing polyacrylamide gel. After electrophoresis, the gel was dried on Whatman 3MM paper and subjected to autoradiography.
 |
RESULTS
|
---|
SNP identification and haplotype reconstruction.
In a previous article, we described a polymorphism located at +245 with respect to the transcription start site in the first intron of the OPN gene, for which we did not find a direct functional role (12). To verify whether other variants, in linkage disequilibrium with the above polymorphism, might play some role on expression or function of the OPN protein, we decided to investigate on polymorphisms located in the 5' flanking region in search for variants which might have a functional role in transcriptional regulation.
To identify naturally occurring sequence variations within OPN regulatory sequences, a fragment flanking the 5' end and including the first intron (2100 to +1015) of the gene was amplified and sequenced in six Italian individuals. A total of nine variants were found: the genomic structure and locations of variant sites are shown in Fig. 1A; Fig. 1B shows the sequence between 462 and +3. Three of these SNPs (+245, +954, +1020) were located in the first intron: the first one (+245) was previously described by us, and the two others were already reported in the Japanese population (17). Six novel SNPs were found in the sequence upstream of the transcription start site and also confirmed in a recent Japanese study (22). Five of the observed variants (443, 156, 66, +245, +2080) were further analyzed in 50 Italian individuals, and allele frequency was determined. Allele frequencies of all polymorphic variants are shown in Table 2.

View larger version (30K):
[in this window]
[in a new window]
|
Fig. 1. Location of single nucleotide polymorphisms (SNPs) in the human osteopontin (OPN) gene regulatory region from 2236 to +1200 and sequence of the OPN promoter. A: SNPs and their respective location are shown (not in scale) relative to the transcription start site (arrow, +1). B: sequence of 500 bp of the OPN promoter: location of polymorphism 443 C/T, G/GG 156, and T/G 66 is shown with nucleotide in bold; sequences of oligonucleotides used in gel shift assay are underlined.
|
|
Calculation of linkage disequilibrium between pairs of variants (D' value), obtained by the Arlequin software, resulted in highly significant values, as shown in Table 3. Haplotypes reconstruction, obtained by the computer program PHASE, resulted in finding 14 haplotypes as shown in Table 4 with their respective frequencies. Considering the 443 C/T, 156 G/GG, and 66 T/G polymorphisms in the immediate 5' flanking region, four haplotype combinations accounted for 96% of the total haplotype frequencies, the remaining 4% consisting of three haplotypes with frequency lower than 2%. The four most frequent haplotypes are indicated as A, B, C, and D, as shown in Fig. 2.

View larger version (19K):
[in this window]
[in a new window]
|
Fig. 2. Most frequent haplotypes of the human osteopontin gene. The four most frequent haplotypes and the derived E and F haplotypes are shown in the diagram, which describes the pGL3 constructs in which they were inserted. These constructs were used for transfection experiments into a panel of cell lines.
|
|
Transfection experiments.
The location of three polymorphic variants in the immediate 500 bp at the 5' end of the OPN gene prompted us to investigate on their effect in promoter activity; therefore, DNA fragments (from 474 to +42 relative to the transcription start site) corresponding to the four common haplotypes were inserted in pGL3 expression plasmids upstream of the luciferase reporter gene (Fig. 2). Transient transfection experiments were performed in a panel of cell lines, chosen because these are related to OPN expression and regulation: HeLa (epithelial, derived from human uterine cervical carcinoma), MCF7 (epithelial, derived from human breast adenocarcinoma), U2OS (osteoblast-like, derived from human osteogenic sarcoma), COS7 (fibroblast, derived from monkey kidney), HEK293 (fibroblast, derived from human embryonic kidney), and HK2 (human immortalized epithelial cells derived from normal proximal convoluted tubule).
Figure 3 shows that in all tested cell lines the A haplotype showed the lowest promoter activity, significantly different from the remaining three haplotypes, with a maximum of 4.87-fold with respect to B haplotype for the HEK293 cell line. Differences among B, C, and D haplotypes, more evident in some of the tested cell lines, did not reach statistical significance. Since the A haplotype differs from B, C, and D at the 66 position, to further confirm that this G/T variant has a major effect on promoter activity, we prepared two additional plasmids by mutagenesis of C and D constructs, carrying G in place of T at the 66 position (see Fig. 2). These plasmids, named E and F, were transfected in COS7 and U2OS cell lines. In both tested cell lines, luciferase activity of haplotypes A, E, and F did not show statistical difference, whereas E vs. C and F vs. D differed significantly (data not shown), consistent with a major role of the T/G variant at 66.

View larger version (17K):
[in this window]
[in a new window]
|
Fig. 3. Transfection of different cell lines with constructs containing haplotypes of the OPN promoter. The indicated cell lines were transiently transfected with constructs containing the 4 most frequent haplotypes fused to the luciferase reporter gene in the pGL3 vector. Values are expressed as luciferase units at various scales after normalization for cotransfected Renilla luciferase control construct pRL-CMV. Data shown represent median values (circles) and first and third quartile (error bars) of luciferase activity; each unit of observation was the result of at least three independent experiments, performed either in triplicate or in duplicate. *P < 0.05. **P < 0.01. ***P < 0.001. All comparisons were performed with respect to haplotype A.
|
|
Transcription factor binding sites and variant alleles.
An analysis of sequences corresponding to the three variant locations was carried out by using computer comparison with consensus binding sites for known transcription factors, both by the MatInspector (Genomatix) online software and by the use of a personal binding site data base inserted in the MacVector program. This analysis suggested putative binding site for the SP1 factor around the 66 SNP, for CBFA1/RUNX2 factor around the 156 small insertion and for the MYT1 zinc finger factor around the 443 SNP.
The nucleotide sequence of the human osteopontin promoter between 170 and 30, compared with the rat, mouse, and porcine sequence, is shown in Fig. 4A in which we observed that the positions of the SNPs 66 and 156 are located in a conserved region. Considering RUNX2 binding sites, which are very important for regulation of OPN expression in bone tissue, we observed that the 156 SNP generates a RUNX2 binding site very close to another RUNX2 binding site located 14 bp downstream, which is highly conserved between species, as shown in Fig. 4B.

View larger version (33K):
[in this window]
[in a new window]
|
Fig. 4. Alignment of nucleotide sequences of the promoter regions from 170 to 30 of the human, mouse, rat, and porcine osteopontin genes. A: the nucleotide numbers correspond to the human sequence. Asterisks have been inserted in sequences to maximize alignment. Dashes in mouse, rat, and porcine sequences indicate nucleotides conserved with respect to the human sequence. Polymorphic variants identified in our study are set in bold and underscore font. B: detail of the 163 to 129 sequence to highlight RUNX2 binding sites (marked with underscore).
|
|
We then designed double-stranded oligonucleotides reproducing the sequences corresponding to the above loci in both allelic forms (see Table 1) and performed EMSA using nuclear extracts from different cell lines.
Figure 5 shows that the 82 to 47 oligonucleotides, including the 66 T or G alleles formed retarded complexes with U2OS nuclear extracts, which were also observed in comparable way with HEK293 nuclear extracts (not shown). These complexes were sequence specific since they could be competed by the respective unlabeled oligonucleotides. Differences in intensity were appreciable for the two alternative probes and cross-competitions highlighted that the oligonucleotide with the T allele contained a higher affinity binding site than the G allele. Addition of specific anti-SP1 and anti-SP3 antibodies in the EMSA caused appearance of a specific supershifted complex for SP1 and disappearance of the specific complex for SP3. A nonrelated control antibody did not affected the EMSA pattern. EMSA with HeLa nuclear extract, known to contain abundant SP1/SP3 factors confirmed this finding (not shown).

View larger version (78K):
[in this window]
[in a new window]
|
Fig. 5. Gel retardation assay with a probe including the 66 SNP. Gel retardation assay with a double-stranded oligonucleotide that reproduces the OPN promoter from 88 to 48 containing the G/T polymorphism, and nuclear extracts from U2OS cells: lanes 7 and 8, in the absence, and lanes 16 and 914, in the presence of competitor oligonucleotides at 50-, 100-, or 250-fold molar excess in the two allelic forms. Presence of antibody against SP1 generates a supershifted complex (lane 17), whereas antibody against SP3 shows disappearance of a complex as indicated with asterisk (lane 18). No effect observed in the absence of antibody or in the presence of a nonspecific antibody (lanes 15 and 16).
|
|
Figure 6 shows the EMSA experiments using the 172 to 137 probes with 156 G or GG alleles, in which a retarded complex was observed. Sequence specificity was demonstrated by competition with the respective unlabeled oligonucleotide. This complex differs for the two alternative alleles, as can be clearly observed in cross-competition, more evident with the HEK293 nuclear extract. A double component is present in the complex when the probe contains the G allele. These two components can be completely competed only by the G oligonucleotide, whereas the GG oligonucleotide competes only for one of the two components, and, when used as a probe, the resulting complex has higher intensity (see Fig. 6C for detail). This was also confirmed with the U2OS nuclear extract, although at lower intensity. As anticipated above, computer analysis revealed that the GG allele could give rise to a putative RUNX2 binding site on the minus strand. Based on this indication, we used an anti-murine Runx2 antibody in the EMSA and obtained a supershifted complex, more abundant with the probe carrying the GG allele (Fig. 6D). This supershifted complex was sequence specific since it was competed by the unlabeled oligonucleotide and was not observed when a nonrelated antibody was used in the assay (not shown). We exclude that the 14 bp downstream RUNX2 binding site (at 140/134) could have any role in this assay since this site is interrupted in the probe oligonucleotide used.

View larger version (106K):
[in this window]
[in a new window]
|
Fig. 6. Gel retardation assay with a probe including the 156 SNP. A and B: a gel retardation assay with a double-stranded oligonucleotide that reproduces the 88 to 48 sequence containing the G/GG polymorphism, and nuclear extracts from HEK293 and U2OS cells, respectively. The two central lanes, in the absence, and flanking lanes, in the presence, of competitor oligonucleotides at 50-, 100-, or 250-fold molar excess in the two allelic forms, as indicated. C: cross-competition with allelic oligonucleotides, with HEK293 nuclear extract. D: supershift obtained with the specific anti-Runx2 antibody.
|
|
EMSA experiments were also performed using the 460 to 425 probes with 443 T or C alleles. This probe, predicted to contain a putative binding site for the MYT1 zinc finger factor (14), gave rise to retarded complexes that differed according to nuclear extracts used and, in part, according to T or C allele (not shown). Lack of specific antibodies prevented us from proving binding of a specific transcription factor.
Cotransfection experiments with SP1 and Runx2 expression vectors.
SP1 overexpression in U2OS cotransfected with haplotypes A and B, which differ at the 66 position only, and a CMV-SP1 expression plasmid resulted in around 10-fold activation of promoter activity with a 1.24-fold higher activation of the B plasmid compared with A (average of two experiments performed in triplicate: 1.211.27), consistent with higher binding affinity of SP1 to the T allele at 66, located in the B haplotype.
To test whether the 156 polymorphism could underlie a different response to the RUNX2 transcription factor, cotransfections with reporter plasmids containing two promoter haplotypes differing at 156 (haplotypes B and D) and an expression vector carrying the Runx2 murine cDNA were performed in U2OS, MCF7, and HEK293 cell lines. When we tested the values obtained by comparing the two haplotypes in the three different cell lines, we obtained two types of information, as shown in Fig. 7: 1) U2OS and MCF7 cells displayed statistically significant response to RUNX2, whereas HEK293 were virtually unresponsive; 2) the haplotype with the GG allele showed higher response than the one with the G allele, in both U2OS (P = 0.0078) and MCF7 cells (P = 0.0001).

View larger version (17K):
[in this window]
[in a new window]
|
Fig. 7. Cotransfection of vectors containing different promoter haplotypes with a Runx2 expression vector. U2OS, HEK293, and MCF7 cells were transiently transfected with constructs containing either haplotype D or B, differing for alleles at 156, in the presence of Runx2 or control vector pCMV5. All transfections were performed in triplicate; luciferase activities for each haplotype are presented as ratio between median of haplotype transfected with Runx2 and haplotype transfected with control vector pCMV5. **P < 0.01. ***P < 0.001. n.s., not significant.
|
|
 |
DISCUSSION
|
---|
Changes of expression of the OPN gene, mostly overexpression, have been described in a variety of conditions; therefore, investigating mechanisms of OPN transcriptional regulation can provide useful information on the molecular basis of expression variation.
Several polymorphisms in the human OPN encoding gene have been identified in different populations: some have been reported in the 5' flanking region in a Japanese population (22), others in exons and introns and in the 3' untranslated region (12, 17). Our work was aimed at verifying whether alleles at polymorphic sites in the 5' flanking region were able to influence the level of promoter activity and characterizing their ability to bind transcription factors.
Our analysis, focused on a 1-kb region upstream of the transcription start site, confirmed the presence of polymorphic variants, the same as those observed in the Japanese population. The reconstruction of haplotypes showed that four of them have frequencies above 10% in our population with a high degree of linkage disequilibrium inside this genomic region.
The functional characterization of three of the promoter polymorphic variants was carried out by preparation of expression plasmids in which the luciferase reporter gene was fused to the OPN promoter sequence in the four most frequent haplotypes found in our population. Transfection of these plasmids in a panel of different cell lines showed that the combination of variants within haplotypes can significantly affect promoter activity. Haplotype A showed consistently lower promoter activity compared with haplotypes B, C, and D in all tested cell lines. This was also confirmed with haplotypes E and F, carrying the G allele at 66, like the A haplotype, which also displayed low activity as the A haplotype, and significantly different from the respective C and D haplotypes, carrying the T allele at 66.
We therefore examined the effect of different alleles on the ability to bind nuclear factors at the level of the three single polymorphic variants.
The SNP at 66 was predicted to bind the SP1 transcription factor both by the knowledge that the sequence between 68 and 59 is an SP1 binding site in the human OPN promoter (39) and by promoter analysis with the MatInspector software. This was confirmed by our EMSA in which we showed that SP1 and SP3 factors recognize this site and that the allele with the T nucleotide inside the SP1 recognition sequence has a higher binding affinity. Consistent with this finding, overexpression of SP1 induced higher promoter activation on a plasmid carrying the T allele compared with the equivalent plasmid carrying the G allele.
The promoter sequence from 94 to 62 is essential for the human OPN promoter activity (12, 39) and, in the portion related to the SP1 binding site, is highly conserved between species with the T allele in it. Thus the finding that a SNP in this sequence is able to reduce promoter activity to at least 50% in different cell lines appears as particularly significant when considering the molecular basis of individual variation of gene expression. Moreover, this SP1 polymorphic site could have an important role in different conditions in which SP1 interacts with other factors such as AP1, which finds an immediately upstream recognition site, or RUNX2, whose binding sites are located less than 100 bp upstream. It is well known that OPN transcription can be activated by various stimuli both through AP1 (20, 28, 31) and RUNX2 (15) factors and that both these factors could interact with SP1 in a cooperative way.
For the SNP at 156, our sequence analysis suggested that a binding site for RUNT factors, although not exactly matching the reported consensus sequence (9), might be created by the presence of two G nucleotides inside a long stretch of Ts, which, on the complementary strand, gives rise to the AACCAAA sequence. We therefore conducted an analysis, by both EMSA and cotransfection experiments with a Runx2 expression plasmid, which showed that an antibody to Runx2 was able to supershift a retarded complex and that Runx2 expression was able to activate the OPN promoter in a differential way, according to the allele. Both the G and the GG containing putative RUNT site are able to bind the factor with an apparent better ability for the GG allele, as shown by higher intensity of the retarded complex detectable when this becomes the only component (Fig. 6C), and higher intensity of the supershifted band. This was verified with nuclear extracts from two cell lines, U2OS and HEK293. Cotransfection experiments showed that expression of Runx2 activated the OPN promoter in U2OS and MCF7 cell lines, whereas no activation was observed in the HEK293 cell line, suggesting that, although a RUNT factor is present in nuclear extracts of HEK293, promoter activation depends on the presence of some additional factor absent in this last cell line. It is interesting to note that besides U2OS, an osteoblast-like cell line derived from osteosarcoma, activation took place in MCF7 cells, derived from human breast cancer, which is a tumor highly prone to bone metastasis. Another interesting observation comes from recent reports about variants in binding sites for the RUNX1 gene, which belongs to the RUNT family of transcription factors as well as RUNX2, found associated with common autoimmune diseases (2). Since the RUNX factors share their DNA binding sites and can heterodimerize, our finding raises various hypotheses on possible regulation of OPN related to susceptibility to disease in which the immune response has a relevant role.
The 443 SNP, at our EMSA, showed differential binding of a yet unidentified nuclear factor. Our sequence analysis suggests that this may be a putative binding site for the MYT1 zinc finger transcription factor, which, up to now, has been reported as involved in neurogenesis and pancreas development (14), or a related factor belonging to the same zinc finger subfamily. Although our results show differential binding of a nuclear factor to the sequence according to the polymorphic alleles, further work is required to clarify this point.
Our results clearly suggest that different haplotypes related to common variants in the promoter of OPN gene may influence the individual degree of gene expression, although other levels of regulation have to be considered especially in the case of this multifunctional protein. Moreover, individual variability in response to various stimuli, such as those triggered by inflammation, necrosis, growth factors, and hormones, is even more complicated, although variability in sequences that are targets of different factors is likely to play a role.
Study of genetic factors involved in susceptibility to complex traits relies on analysis of DNA variants at the susceptibility locus and, subsequently, on correlation of alleles with functional properties. The finding that haplotypes in the OPN encoding gene affect promoter activity will be of great help in the study of diseases for which osteopontin is a good candidate as a susceptibility gene.
 |
ACKNOWLEDGMENTS
|
---|
We thank Drs. Renata Bocciardi and Marco Musso for helpful discussion and review of the manuscript and Loredana Velo for secretarial assistance.
This work was supported by a Fondo per gli Investimenti della Ricerca di Base (FIRB) grant of the Italian Ministry of University to R. Ravazzolo.
 |
FOOTNOTES
|
---|
Article published online before print. See web site for date of publication (http://physiolgenomics.physiology.org).
Address for reprint requests and other correspondence: R. Ravazzolo, G. Gaslini Institute, Laboratory of Molecular Genetics, Largo G. Gaslini 5, 16148 Genoa, Italy (E-mail: rravazzo{at}unige.it).
doi:10.1152/physiolgenomics.00138. 2004.
 |
REFERENCES
|
---|
- Agnihotri R, Crawford HC, Haro H, Matrisian LM, Havrda MC, and Liaw L. Osteopontin, a novel substrate for matrix metalloproteinase-3 (stromelysin-1) and matrix metalloproteinase-7 (matrilysin). J Biol Chem 276: 2826128267, 2001.[Abstract/Free Full Text]
- Alarcón-Riquelme ME. A RUNX trio with a taste for autoimmunity. Nat Genet 35: 299300, 2003.[CrossRef][ISI][Medline]
- Ashkar S, Weber GF, Panoutsakopoulou V, Sanchirico ME, Jansson M, Zawaideh S, Rittling SR, Denhardt DT, Glimcher MJ, and Cantor H. Eta-1 (osteopontin): an early component of type-1 (cell-mediated) immunity. Science 287: 860864, 2000.[Abstract/Free Full Text]
- Chabas D, Baranzini SE, Mitchell D, Bernard CC, Rittling SR, Denhardt DT, Sobel RA, Lock C, Karpuj M, Pedotti R, Heller R, Oksenberg JR, and Steinman L. The influence of the proinflammatory cytokine, osteopontin, on autoimmune demyelinating disease. Science 294: 17311735, 2001.[Abstract/Free Full Text]
- Chiocchetti A, Indelicato M, Bensi T, Mesturini R, Giordano M, Sametti S, Castelli L, Bottarel F, Mazzarino MC, Garbarini L, Giacopelli F, Valesini G, Santoro C, Dianzani I, Ramenghi U, and Dianzani U. High levels of osteopontin associated with polymorphisms in its gene are a risk factor for development of autoimmunity/lymphoproliferation. Blood 103: 13761382, 2004.[Abstract/Free Full Text]
- Denhardt DT, Noda M, ORegan AW, Pavlin D, and Berman JS. Osteopontin as a means to cope with environmental insults: regulation of inflammation, tissue remodeling, and cell survival. J Clin Invest 107: 10551061, 2001.[Free Full Text]
- Dignam JD, Lebovitz RM, and Roeder RG. Accurate transcription initiation by RNA polymerase II in a soluble extract from isolated mammalian nuclei. Nucleic Acids Res 11: 14751489, 1983.[Abstract]
- Ducy P, Starbuck M, Priemel M, Shen J, Pinero G, Geoffroy V, Amling M, and Karsenty G. A Cbfa1-dependent genetic pathway controls bone formation beyond embryonic development. Genes Dev 13: 10251036, 1999.[Abstract/Free Full Text]
- Ducy P, Zhang R, Geoffroy V, Ridall AL, and Karsenty G. Osf2/Cbfa1: a transcriptional activator of osteoblast differentiation. Cell 89: 747754, 1997.[ISI][Medline]
- Fisher LW and Fedarko NS. Six genes expressed in bones and teeth encode the current members of the SIBLING family of proteins. Connect Tissue Res 44, Suppl 1: 3340, 2003.
- Forton AC, Petri MA, Goldman D, and Sullivan KE. An osteopontin (SPP1) polymorphism is associated with systemic lupus erythematosus. Hum Mutat 19: 459, 2002.
- Giacopelli F, Rosatto N, Divizia MT, Cusano R, Caridi G, and Ravazzolo R. The first intron of the human osteopontin gene contains a C/EBP-beta-responsive enhancer. Gene Expr 11: 95104, 2003.[ISI][Medline]
- Higuchi R, Krummel B, and Saiki RK. A general method of in vitro preparation and specific mutagenesis of DNA fragments: study of protein and DNA interaction. Nucleic Acids Res 16: 73517367, 1988.[Abstract]
- Gu G, Wells JM, Dombkowski D, Preffer F, Aronow B, and Melton DA. Global expression analysis of gene regulatory pathways during endocrine pancreatic development. Development 131: 165179, 2004.[Abstract/Free Full Text]
- Inman CK and Shore P. The osteoblast transcription factor Runx2 is expressed in mammary epithelial cells and mediates osteopontin expression. J Biol Chem 278: 4868448689, 2003.[Abstract/Free Full Text]
- Isoda K, Nishikawa K, Kamezawa Y, Yoshida M, Kusuhara M, Moroi M, Tada N, and Ohsuzu F. Osteopontin plays an important role in the development of medial thickening and neointimal formation. Circ Res 91: 7782, 2002.[Abstract/Free Full Text]
- Iwasaki H, Shinohara Y, Ezura Y, Ishida R, Kodaira M, Kajita M, Nakajima T, Shiba T, and Emi M. Thirteen single-nucleotide polymorphisms in the human osteopontin gene identified by sequencing of the entire gene in Japanese individuals. J Hum Genet 46: 544546, 2001.[CrossRef][ISI][Medline]
- Jansson M, Panoutsakopoulou V, Baker J, Klein L and Cantor H. Attenuated experimental autoimmune encephalomyelitis in eta-1/osteopontin-deficient mice. J Immunol 168: 20962099, 2002.[Abstract/Free Full Text]
- Kikuchi K, Tanaka A, Miyakawa H, Kawashima Y, Kawaguchi N, Matsushita M, and Gershwin ME. Eta-1/osteopontin genetic polymorphism and primary biliary cirrhosis. Hepatol Res 26: 8790, 2003.[CrossRef][ISI][Medline]
- Kim HJ, Lee MH, Kim HJ, Shin HI, Choi JY, and Ryoo HM. Okadaic acid stimulates osteopontin expression through de novo induction of AP-1. J Cell Biochem 87: 93102, 2002.[CrossRef][ISI][Medline]
- Liaw L, Birk DE, Ballas CB, Whitsitt JS, Davidson JM, and Hogan BL. Altered wound healing in mice lacking a functional osteopontin gene (spp1). J Clin Invest 101: 14681478, 1998.[Abstract/Free Full Text]
- Mochida S, Hashimoto M, Matsui A, Naito M, Inao M, Nagoshi S, Nagano M, Egashira T, Mishiro S, and Fujiwara K. Genetic polymorphisms in promoter region of osteopontin gene may be a marker reflecting hepatitis activity in chronic hepatitis C patients. Biochem Biophys Res Commun 313: 10791085, 2004.[CrossRef][ISI][Medline]
- Niino M, Kikuchi S, Fukazawa T, Yabe I, and Tashiro K. Genetic polymorphisms of osteopontin in association with multiple sclerosis in Japanese patients. J Neuroimmunol 136: 125129, 2003.[CrossRef][ISI][Medline]
- Ohshima S, Kobayashi H, Yamaguchi N, Nishioka K, Umeshita-Sasai M, Mima T, Nomura S, Kon S, Inobe M, Uede T, and Saeki Y. Expression of osteopontin at sites of bone erosion in a murine experimental arthritis model of collagen-induced arthritis: possible involvement of osteopontin in bone destruction in arthritis. Arthritis Rheum 46: 10941101, 2002.[CrossRef][ISI][Medline]
- Patrone G, Puliti A, Bocciardi R, Ravazzolo R, and Romeo G. Sequence and characterisation of the RET proto-oncogene 5' flanking region: analysis of retinoic acid responsiveness at the transcriptional level. FEBS Lett 419: 7682, 1997.[CrossRef][ISI][Medline]
- Persy VP, Verhulst A, Ysebaert DK, De Greef KE, and De Broe ME. Reduced postischemic macrophage infiltration and interstitial fibrosis in osteopontin knockout mice. Kidney Int 63: 543553, 2003.[CrossRef][ISI][Medline]
- Petrow PK, Hummel KM, Schedel J, Franz JK, Klein CL, Muller-Ladner U, Kriegsmann J, Chang PL, Prince CW, Gay RE, and Gay S. Expression of osteopontin messenger RNA and protein in rheumatoid arthritis: effects of osteopontin on the release of collagenase 1 from articular chondrocytes and synovial fibroblasts. Arthritis Rheum 43: 15971605, 2000.[CrossRef][ISI][Medline]
- Renault MA, Jalvy S, Belloc I, Pasquet S, Sena S, Olive M, Desgranges C, and Gadeau AP. AP-1 is involved in UTP-induced osteopontin expression in arterial smooth muscle cells. Circ Res 93: 674681, 2003.[Abstract/Free Full Text]
- Rittling SR and Denhardt DT. Osteopontin function in pathology: lessons from osteopontin-deficient mice. Exp Nephrol 7: 103113, 1999.[CrossRef][ISI][Medline]
- Saitoh Y, Kuratsu J, Takeshima H, Yamamoto S, and Ushio Y. Expression of osteopontin in human glioma. Its correlation with the malignancy. Lab Invest 72: 5563, 1995.[ISI][Medline]
- Sakata R, Minami S, Sowa Y, Yoshida M, and Tamaki T. Trichostatin A activates the osteopontin gene promoter through AP1 site. Biochem Biophys Res Commun 315: 959963, 2004.[CrossRef][ISI][Medline]
- Schneider S, Roessli D, and Excoffier L. Arlequin ver. 2.000: A Software for Population Genetics Data Analysis. Genetics and Biometry Laboratory, University of Geneva, Switzerland, 2000.
- Sodek J, Ganss B, and McKee MD. Osteopontin. Crit Rev Oral Biol Med 11: 279303, 2000.[Abstract]
- Speer MY, McKee MD, Guldberg RE, Liaw L, Yang HY, Tung E, and Karsenty G, and Giachelli CM. Inactivation of the osteopontin gene enhances vascular calcification of matrix Gla protein-deficient mice: evidence for osteopontin as an inducible inhibitor of vascular calcification in vivo. J Exp Med 196: 10471055, 2002.[Abstract/Free Full Text]
- Stephens M, Smith NJ, and Donnelly P. A new statistical method for haplotype reconstruction from population data. Am J Hum Genet 68: 978989, 2001.[CrossRef][ISI][Medline]
- Stephens M and Donnelly P. A comparison of bayesian methods for haplotype reconstruction from population genotype data. Am J Hum Genet 73: 11621169, 2003.[CrossRef][ISI][Medline]
- Tamura A, Shingai M, Aso N, Hazuku T, and Nasu M. Osteopontin is released from the heart into the coronary circulation in patients with a previous anterior wall myocardial infarction. Circ J 67: 742744, 2003.[CrossRef][ISI][Medline]
- Trueblood NA, Xie Z, Communal C, Sam F, Ngoy S, Liaw L, Jenkins AW, Wang J, Sawyer DB, Bing OH, Apstein CS, Colucci WS, and Singh K. Exaggerated left ventricular dilation and reduced collagen deposition after myocardial infarction in mice lacking osteopontin. Circ Res 88: 10801087, 2001.[Abstract/Free Full Text]
- Wang D, Yamamoto S, Hijiya N, Benveniste EN, and Gladson CL. Transcriptional regulation of the human osteopontin promoter: functional analysis and DNA-protein interactions. Oncogene 19: 58015809, 2000.[CrossRef][ISI][Medline]
- Weber GF. The metastasis gene osteopontin: a candidate target for cancer therapy. Biochim Biophys Acta 1552: 6185, 2001.[CrossRef][ISI][Medline]
- Weber GF, Zawaideh S, Hikita S, Kumar VA, Cantor H, and Ashkar S. Phosphorylation-dependent interaction of osteopontin with its receptors regulates macrophage migration and activation. J Leukoc Biol 72: 752761, 2002.[Abstract/Free Full Text]
- Xie Y, Sakatsume M, Nishi S, Narita I, Arakawa M, and Gejyo F. Expression, roles, receptors, and regulation of osteopontin in the kidney. Kidney Int 60: 16451657, 2001.[CrossRef][ISI][Medline]
- Yamamoto N, Sakai F, Kon S, Morimoto J, Kimura C, Yamazaki H, Okazaki I, Seki N, Fujii T, and Uede T. Essential role of the cryptic epitope SLAYGLR within osteopontin in a murine model of rheumatoid arthritis. J Clin Invest 112: 181188, 2003.[Abstract/Free Full Text]
- Yamate T, Tsuji H, Amasaki N, Iguchi M, Kurita T, and Kohri K. Analysis of osteopontin DNA in patients with urolithiasis. Urol Res 28: 159166, 2000.[CrossRef][ISI][Medline]
- Yoshitake H, Rittling SR, Denhardt DT, and Noda M. Osteopontin-deficient mice are resistant to ovariectomy-induced bone resorption. Proc Natl Acad Sci USA 96: 81568160, 1999.[Abstract/Free Full Text]
- Zhang W, Ou J, Inagaki Y, Greenwel P, and Ramirez F. Synergistic cooperation between Sp1 and Smad3/Smad4 mediates transforming growth factor ß1 stimulation of
2(I)-collagen (COL1A2) transcription. J Biol Chem 275: 3923739245, 2000.[Abstract/Free Full Text]