|
|
||||||||
Molecular Diagnostics and Genetics |
1 Institut für Laboratoriums- und Transfusionsmedizin, Herz- und Diabeteszentrum Nordrhein-Westfalen, Universitätsklinik der Ruhr-Universität Bochum, 32545 Bad Oeynhausen, Germany.
2 Dermatologische Klinik, Krankenhaus Bethesda, 57258 Freudenberg, Germany.
3 Human genes: ABCC6, ATP-binding cassette subfamily C, member 6; CFTR, cystic fibrosis transmembrane conductance regulator (ATP-binding cassette subfamily C, member 7); SPP1, secreted phosphoprotein 1 (previously OPN, osteopontin); and XYLT1, XYLT2: xylosyltransferase I and II.
aAddress correspondence to this author at: Institut für Laboratoriums- und Transfusionsmedizin, Herz- und Diabeteszentrum Nordrhein-Westfalen, Georgstraße 11, 32545 Bad Oeynhausen, Germany. Fax +49-5731-97-2013; e-mail cgoetting{at}hdz-nrw.de.
| Abstract |
|---|
|
|
|---|
Methods: We screened an
2-kb region spanning the theoretical promoter of the SPP1 gene for sequence variations by denaturing HPLC and direct sequencing in 93 PXE patients. Sequence variations with a prevalence >5% were genotyped in 93 age- and sex-matched healthy controls. Statistical and haplotype association analyses were performed using Fisher exact test, PHASE v2.1.1, and Haploview 3.2.
Results: Mutational screening revealed 9 different sequence variations. Three SPP1 promoter polymorphisms (c.1748A>G, c.155_156insG, and c.244_245insTG) were significantly more frequent in PXE patients than in 93 age- and sex-matched healthy controls (Pcorrected < 0.05 each). The odds ratios (95% CI) for PXE among carriers of the 3 alleles were, respectively, 2.16 (1.343.48), 2.41 (1.513.82), and 1.97 (1.233.15). Haplotype analysis of 6 SPP1 promoter polymorphisms revealed 1 haplotype to be significantly reduced among PXE patients (Pcorrected = 0.035, odds ratio 1.80, 95% CI 1.192.71).
Conclusions: Polymorphisms in the SPP1 promoter are secondary genetic risk factors contributing to PXE susceptibility.
| Introduction |
|---|
|
|
|---|
The clinical course of PXE is highly variable. Although environmental influences may modify disease outcome, there might also exist secondary genetic variations (15). "Modifier genes" have been detected for cystic fibrosis, another heritable disorder with variable disease onset and progression (16). Cystic fibrosis is caused by mutations in CFTR (cystic fibrosis transmembrane conductance regulator), member 7 of the same ABC transporter subfamily as ABCC6. A candidate gene for PXE susceptibility is SPP1 (secreted phosphoprotein 1; previously OPN, osteopontin). SPP1 is a secreted, highly acidic phosphoprotein that is involved in immune cell activation, wound healing, and bone morphogenesis (17) and plays a major role in regulating mineralization processes in various tissues. Increased SPP1 expression is often associated with pathological calcification. Furthermore, SPP1 is a constitutive component of human skin and aorta, where it is localized to the elastic fiber and hypothesized to prevent calcification in the fibers (18). Skin biopsies from PXE patients show higher expression of SPP1 than do samples from unaffected regions or from healthy individuals (19). SPP1 expression is increased in mice suffering from dystrophic cardiac calcification (20), leading to the suggestion that high SPP1 expression is influenced by a trans-activator gene, the Dyscalc1 locus on chromosome 7 (20). The mouse abcc6 gene was identified as a potential candidate gene in this region (21).
SPP1 is a predominantly transcriptional regulated gene, and the SPP1 promoter is highly conserved among different species (22). Several polymorphisms in the SPP1 gene affect SPP1 expression and have been associated with various disorders, e.g., systemic lupus erythematosus and arteriosclerosis (23)(24)(25)(26).
We put forward the hypothesis that sequence variations in the SPP1 promoter region might account for the higher SPP1 expression observed in PXE patients and therefore promote disease outcome. We present data from a case-control association study on German PXE patients and an age- and sex-matched normal population with 6 sequence variations spanning the whole SPP1 promoter region.
| Materials and Methods |
|---|
|
|
|---|
|
dna extraction and mutational analysis in the promoter region of the spp1 gene
Genomic DNA was extracted from 200 µL EDTA-anticoagulated blood using the QIAamp blood reagent set (Qiagen) according to the manufacturers instructions. The entire regulatory region from 1864 to +317 of the SPP1 gene was initially amplified by PCR and then analyzed for sequence variations by partially denaturing HPLC (dHPLC). Nucleotide numbering refers to the SPP1 genomic DNA sequence (GenBank accession no. NT_016354), with the first nucleotide of exon 1 as the transcription initiation start site, referred to as nucleotide +1.
PCR was performed in a 50-µL reaction volume, containing
65 ng genomic DNA, 25 pmol of each primer (Biomers), 1.5 units HotStar Taq DNA polymerase (Qiagen) in 1 x PCR buffer supplied with the enzyme, and 0.25 mmol/L of each dNTP (Promega). The primer sequences, annealing temperatures, and sizes of the PCR products are summarized in Table 1 in the Data Supplement that accompanies the online version of this article at http://www.clinchem.org/content/vol53/issue5. The PCR conditions were as follows: initial denaturation at 95 °C for 15 min, 35 cycles of denaturation at 94 °C for 1 min, annealing for 1 min, and extension at 72 °C for 1 min, and a final extension at 72 °C for 15 min. dHPLC analysis was carried out on an automated HPLC device equipped with a DNA separation column (Wave System, Transgenomic) as described (28). Oven temperatures and the initial and final concentrations of buffer B for dHPLC analysis are given in Table 1 of the online Data Supplement. PCR products of samples showing aberrant peaks were purified via exonuclease I (1 unit) and shrimp alkaline phosphatase (1 unit) treatment per 5 µL amplicon at 37 °C for 30 min. Inactivation of the enzymes was carried out at 80 °C for 15 min. The purified PCR product served as a template for mutation identification by direct sequencing on both strands. DNA sequencing of the PCR products was performed on an ABI Prism 310 capillary sequencer using the Big Dye Terminator v1.1 cycle sequencing reagent set (Perkin-Elmer Applied Biosystems), 3 µL purified PCR product, and 2.5 pmol of the same PCR primers used for amplification, in a total reaction volume of 20 µL. To exclude the possibility of amplification errors resulting from HotStar Taq DNA polymerase, PCR and direct sequencing were repeated whenever a sequence variation was identified.
genotyping of common spp1 promoter polymorphism
Sequence variations with a prevalence of >5% were genotyped in cases and controls by dHPLC or restriction fragment length polymorphism analysis. dHLPC analysis was used for the detection of the polymorphisms c.155_156insG, c.66T>G, and c.244_245 insTG. Samples initially genotyped as heterozygous were reanalyzed without mixing the sample with a known homozygous control. This approach allows for determination of whether the sample was heterozygous or homozygous for the variant. Genotyping for the polymorphisms c.1776T>C, c.1748A>G, and c.443C>T was performed by using mismatch PCR for construction of restriction sites, followed by digestion with an appropriate restriction enzyme. For the genotyping of c.616G>T, allele-specific PCRs were established. Primer sequences for mismatch and allele-specific PCR, annealing temperatures, sizes of the PCR products, and restriction enzymes are listed in Table 2 of the online Data Supplement. The restriction enzymes were obtained from New England Biolabs, and digestions were performed according to the manufacturers instructions.
statistical analysis and power calculations
All polymorphisms were tested for confirmation with Hardy-Weinberg expectations in both cohorts analyzed. Allele and haplotype frequencies were compared between cases and controls using the Fisher exact test. The association of each polymorphism with PXE was measured by the odds ratio (OR) and 95% CI. Multiple testing correction was performed using the Bonferroni method. P <0.05 was considered significant after Bonferroni correction. All tests were executed with GraphPad Prism 4.0 (GraphPad Software). Power calculations were performed using the program developed by Skol et al. (29).
linkage disequilibrium structure and identification of haplotype blocks
Determination of linkage disequilibrium (LD) and haplotype blocks and frequencies were performed by using 2 different validated programs, Haploview 3.2 and PHASE v2.1.1 (30)(31), and comparing them. Haplotype blocks were defined according to the "spine of LD" setting in Haploview software, which is on the basis of each end marker of a block having a D' value of >0.8 (30).
| Results |
|---|
|
|
|---|
5% (c.1776T>C, c.1625A>G, and c.1282A>G) were not analyzed further. The remaining 6 SPP1 promoter polymorphisms were genotyped in an age- and sex-matched control cohort.
|
association analysis of individual sequence variations
Comparison of the allelic frequencies of the detected SPP1 promoter polymorphisms between PXE patients and healthy controls revealed the 4 variants c.1748A>G, c.155_156insG, c.66T>G, and c.244_245insTG to be significantly more frequent in the PXE group (P <0.05 each) (Table 2
). The association with PXE was significant for allele c.1748G (frequency in PXE patients vs controls: 0.333 vs 0.188, P = 0.0020), allele c.155_156GG (0.382 vs 0.204, P = 0.0002), allele c.66G (0.296 vs 0.177, P = 0.0102), and allele c.244_245TGTG (0.328 vs 0.199, P = 0.0066). After correction for multiple testing according to the Bonferroni method, only the 3 SPP1 promoter polymorphisms c.1748A>G, c.155_156insG, and c.244_245insTG remained significantly associated with PXE (Pcorrected < 0.05) (Table 2
). The presence of the disease-associated allele conferred an OR of 2.16 (95% CI 1.343.48) for allele c.1748G, 2.41 (1.513.82) for allele c.155_156GG, and 1.97 (1.233.15) for allele c.244_245TGTG.
To calculate the power of our sample size and genetic models, we used the algorithm of Skol et al. (29). The results indicate that our sample has a sufficient size to detect associations with a power of 79%90% for a relative risk of
2.0 for the most frequent polymorphisms (minor allele frequency 25%45%, significance
= 5%), assuming an additive or multiplicative model for PXE. Assuming a dominant or recessive model, our study had a power of 37%59% or 16%45%, respectively. For the less abundant sequence variants (minor allele frequency <8%), our study had a power ranging from 26%42%, assuming a multiplicative, additive, or dominant model if the relative risk exceeds 2.0; assuming a recessive model for PXE, the power of our study was 5% for these rare SPP1 promoter polymorphisms.
determination of ld structure and haplotype blocks
LD and haplotype blocks were evaluated in PXE patients and controls using Haploview 3.2 and PHASE v2.1.1 (30)(31). Significant LD was observed among the 6 SPP1 promoter polymorphisms, with a minor allele frequency of >5% as shown in Fig. 1
. Determination of LD structure revealed that the 3 variants c.1748A>G, c.155_156insG, and c.244_245insTG were in perfect linkage disequilibrium (D' >0.9), as determined by other groups (24).
|
The Haploview results showed reconstruction of 8 different haplotype blocks in the pooled sample of cases and controls, but only 5 had a frequency of >2% (Table 3
). Comparison of haplotype block frequencies among PXE patients and controls showed a decreased frequency of the major haplotype A (Table 3
) in the PXE group (Pcorrected = 0.0345, OR 1.8, 95% CI 1.192.71). The haplotype showing the highest frequency in the PXE patients (haplotype B in Table 3
) carries the alleles c.1748G, c.155_156GG, and c.244_245TGTG, which were individually increased in the PXE patients (Pcorrected < 0.05 each) (Table 2
).
|
We repeated analysis using the PHASE algorithm to compare haplotype patterns in PXE patients and controls. The PHASE program generated 14 different haplotype blocks (data not shown). The 5 most frequent haplotype patterns and the estimated frequencies were the same as determined by the Haploview algorithm. The remaining 9 haplotypes had frequencies <1%. The overall haplotype frequency distribution was significantly different between patients and controls (P = 0.04).
association of spp1 promoter polymorphisms with clinical features
Allelic frequencies of the 3 SPP1 promoter variants were analyzed in subgroups of the PXE patient groups to evaluate an association with the clinical features shown in Table 1
. We found no significant association between age, age at PXE onset, and number or kind of organs involved.
| Discussion |
|---|
|
|
|---|
A number of genes are likely to contribute to PXE susceptibility. Previously we showed that polymorphisms in the xylosyltransferase I and II genes (XYLT1 and XYLT2) result in a severe disease course of PXE (33). In this study, we analyzed SPP1, an interesting candidate gene with important functions in the regulation of biological calcification, for sequence variations in the proximal promoter region. Altered expression of SPP1 has been found in the dermis of PXE patients analyzed by immunoelectron microscopy (19). Recent studies have described strong SPP1 expression due to the Dyscalc1 locus in mice suffering from dystrophic cardiac calcification, with abcc6 gene as a potential candidate gene (20)(21). The abcc6 gene has just been excluded, however, emphasizing that only the N-terminal part of the gene was screened for sequence variations and that most PXE-associated ABCC6 mutations were found in the C-terminal part (5)(6)(34). The results of our investigation suggest that German PXE patients represent a different distribution of SPP1 promoter polymorphisms than an age- and sex-matched control cohort. The c.1748G, c.155_156insG, and c.244_245insTG alleles appear to be significantly more common in PXE patients. Haplotype analysis revealed 1 haplotype to be significantly reduced among PXE patients, whereas another haplotype, bearing the disease-associated alleles, was more often found in the PXE group (although not statistically significant). Conclusively, these 3 SPP1 promoter polymorphisms and the haplotype combining these disease-associated alleles could be interpreted as a genetic risk pattern for PXE.
Until now, no functional studies have been carried out with the SPP1 promoter polymorphism c.1748A>G. The polymorphic variant c.244_245insTG did not have a major effect on regulation of SPP1 gene expression (35). Analysis of the SPP1 promoter sequence revealed putative transcription factor binding sites for SP1 around c.66, for CBFA1/RUNX2 around c.155, and for MYT1 zinc finger factor at c.443 (Fig. 2
) (23). The polymorphism c.155_156insG generates a RUNX2-binding site. Another RUNX2-binding site was found 14 bp downstream and is highly conserved between species, and the RUNX2 factor was shown to bind better to the c.155_156GG allele than to the c.155_156G allele (23). RUNX2-binding sites are very important for regulation of SPP1 expression in bone tissue (36). A constitutive expression of RUNX2, combined with a supplementation of glucocorticoid hormones, resulted in a strong upregulation of SPP1 expression and finally in a biological matrix mineralization of primary dermal fibroblasts (37). Reporter gene expression experiments with the SPP1 promoter polymorphisms c.443C>T, c.155_156insG, and c.66T>G revealed significantly altered expression of SPP1 (23)(25). All haplotype combinations constructed by these 3 sequence variants resulted in increased reporter gene expression. The strongest expression was conferred by the G-insertion in position c.155 in combination with the c.66T allele (25). These results suggest that haplotypes should be preferred for disease association studies instead of single variants.
|
Although we identified polymorphisms and haplotype patterns in the SPP1 promoter that may play a significant role in the pathogenesis of PXE by altering gene expression, the mechanisms by which these variants predispose patients to PXE are unknown. The higher expression of SPP1 found in the dermis of PXE could result from the genetic risk pattern we identified in this study or could just accompany the calcification process. It is still uncertain whether a higher expression of SPP1 prevents or promotes the mineralization process. Phosphoproteins such as SPP1 were shown to initiate mineral formation when bound to collagen and to inhibit crystal growth when in solution (38). Therefore, SPP1 could induce calcification when bound to the scaffold of elastic fibers, as shown by Baccarani-Contri et al. (19). Whether SPP1 is bound to the elastin protein or to another component of the elastic fiber, maybe hyaluronic acid, is still unknown. Thus, the identification of the physiological substrate of MRP6 will shed light on the role of SPP1 and other extracellular matrix proteins in PXE.
There are some limitations to our study. Although the size of the PXE patient cohort was small, the power of the present study was adequate to detect an association of SPP1 promoter polymorphisms and susceptibility to PXE reliably. However, it cannot be totally excluded that relationships of smaller magnitude were missed in our analysis. Especially for the less frequent variants with a minor allele frequency <8%, the power did not exceed 80%, indicating a possibility of false-negative results. Our cohort size was not large enough to detect associations when assuming a dominant or recessive disease model. Another retrospective study analyzing the association of SPP1 promoter polymorphisms in PXE patients is now necessary to determine whether these polymorphisms are indeed a genetic risk factor for PXE.
We have not yet analyzed the effect of these SPP1 promoter variants on the expression of SPP1, for example in plasma or serum, for several reasons. The currently available ELISA assays are suitable for determination of SPP1 concentrations only in plasma samples, which were not available for the PXE patients and controls in this study. Despite several publications that report SPP1 concentrations in serum, the results from currently used ELISA systems should be interpreted with caution because of highly variable and contradictory SPP1 concentrations in plasma and serum samples, even among healthy controls and between identical ELISA systems (39).
In summary, we have observed an association between SPP1 promoter variants and PXE. Our findings add significant support to the role of proteins actively involved in regulating calcification processes in PXE and underscore the importance of the analysis of genegeneenvironment interactions in understanding the development of complex phenotypes such as PXE. Understanding the whole genetic risk pattern or patterns of PXE may provide new insights into the pathogenesis of the disease and eventually provide opportunities for its treatment, which are limited at present (1). For example, regulating the activity of the SPP1 signaling pathway to modulate its effect on calcification may be possible. Further studies are required to analyze the exact role of SPP1 in the pathophysiology of PXE.
| Acknowledgments |
|---|
Financial disclosures: None declared.
Acknowledgements: We are grateful to all the PXE patients and their relatives, whose cooperation made this study possible. Furthermore, we thank Peter Hof, chairman of the Selbsthilfegruppe für PXE Erkrankte Deutschlands e.V., and the members of the clinical ambulance for PXE at the Bethesda hospital in Freudenberg, Germany. We thank Alexandra Adam and Marlen Ewald for their excellent technical assistance and Sarah L. Kirkby for her linguistic advice.
| Footnotes |
|---|
| References |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |