|
|
||||||||
elj2
Kocjan2
Faculty of Pharmacy, Chair of Clinical Biochemistry, University of Ljubljana, Askerceva cesta 7, SI-1000 Ljubljana, Slovenia 1 Rheumatology Section, Molecular Medicine Centre, Western General Hospital, University of Edinburgh, Edinburgh EH4 2XU, UK 2 Department of Endocrinology, Diabetes and Metabolic Diseases, University Medical Centre, Zaloska cesta 2, SI-1000 Ljubljana, Slovenia
(Correspondence should be addressed to J Marc Email: janja.marc{at}ffa.uni-lj.si)
| Abstract |
|---|
|
|
|---|
B ligand, is one of the candidate genes for the genetic susceptibility to osteoporosis. As variations in the TNFSF11 gene promoter could alter its expression, the aim of the study was to evaluate the functional influence of three polymorphisms in the promoter and to investigate their association with bone mineral density (BMD) and biochemical markers in postmenopausal women. A total of 404 postmenopausal women were genotyped for the presence of TNFSF11 gene promoter polymorphisms –290C>T, –643C>T and –693G>C. Two common haplotypes, CCG and TTC, which occur in 44.3 and 49.3% of subjects respectively, were subjected to functional analysis. Amplified fragments were cloned into pGL3-basic reporter plasmid, which was co-transfected with pRL-TK plasmid into HEK293 cells. Dual luciferase reporter assay was performed. BMD and biochemical markers were measured. Reporter gene analysis showed significantly higher luciferase activity in CCG than in TTC haplotype (P=0.018). Both showed association with lumbar spine BMD (BMD-ls; P=0.005 and 0.007 for TTC and CCG respectively), whereas in femoral neck there was no association with BMD. In postmenopausal osteoporosis, association with BMD-ls was established in –290C>T, –643C>T and –693G>C (P values: 0.001, 0.041 and 0.013 respectively). Association with femoral neck BMD was shown in –693G>C (P=0.049). No association was found with biochemical markers in any of the groups. Our results suggest that in postmenopausal osteoporosis, TNFSF11 gene promoter polymorphisms –290C>T, –643C>T and –693G>C play a functional role in the genetic regulation of BMD.
| Introduction |
|---|
|
|
|---|
Genetic ablation experiments have shown that members of the tumour necrosis factor (TNF) ligand and receptor (TNFR) superfamilies are essential factors of osteoclast biology and bone metabolism, namely receptor activator of nuclear factor-
B ligand (RANKL; Lacey et al. 1998), RANK (Anderson et al. 1997) and osteoprotegerin (OPG; Simonet et al. 1997). The TNF superfamily member 11 (TNFSF11) gene encodes for RANKL protein and is located on human chromosome 13q14 (Walsh & Choi 2003).
RANKL expression is regulated by various hormones (glucocorticoids, vitamin D, oestrogen and parathyroid hormone), cytokines (transforming growth factor-β, TNF-
, interleukins 1, 4, 6, 11 and 17, and prostaglandin E2) and various mesenchymal transcription factors (TFs) such as cbfa-1 and peroxisome proliferators-activated receptor
(Hofbauer & Heufelder 2001). The osteoblastic protein of the Wnt signalling pathway, β-catenin, also regulates RANKL and OPG expression in mice and is therefore involved in the process of osteoclastogenesis (Glass et al. 2005, Holmen et al. 2005).
As RANKL is the final effector of the osteoclastogenesis pathway, it is definitely one of the candidate genes for the regulation of susceptibility to osteoporosis (Ralston & Crombrugghe 2006). In a 2 kb 5'-flanking region of human TNFSF11 gene, several motifs with significant sequence similarity to TFs (>85%) have already been identified (Roccisana et al. 2004). SNPs in the promoter region could influence the binding of TFs and gene expression, and finally impact on protein synthesis. A preliminary study of TNFSF11 gene promoter SNPs in 115 postmenopausal women showed the association of –290C>T SNP with femoral neck bone mineral density (BMD-fn), whereas SNPs –643C>T and –693G>C showed no association with BMD (Mencej et al. 2006). Their genetic influence on BMD has been evaluated, but no functional analysis of the studied SNPs has been performed. Therefore, we investigated whether the studied SNPs are the functional cause of the associated phenotype or, possibly, just in linkage disequilibrium (LD) with another site. Thus, the aims of the present study were: 1) to assess the functional role of certain TNFSF11 gene promoter haplotypes in transcriptional regulation of TNFSF11 gene expression and 2) to explore the potential associations of TNFSF11 gene promoter haplotypes and SNPs with BMD and biochemical markers of bone turnover in osteoporotic and non-osteoporotic postmenopausal women.
| Materials and methods |
|---|
|
|
|---|
We evaluated 404 postmenopausal women, aged 42–91 years, who were referred to the outpatient departments of the University Medical Centre, Ljubljana, or the General Hospital, Celje, for BMD measurement. Each patient was examined clinically and routine biochemical tests were performed to exclude systemic and metabolic bone diseases other than primary osteoporosis. None had previously taken any drug known to influence bone metabolism. Osteoporosis was defined according to the World Health Organization criteria (WHO 1994). The study was approved by the Ethical Committee of the Republic of Slovenia and written informed consent was obtained from each patient participating in the study.
Genotyping
Genomic DNA was isolated from peripheral blood leukocytes using standard methodology. To detect –643C>T (rs9533156) and –693G>C (rs9533155) SNPs, TNFSF11 gene promoter fragment (–888 to –337 bp) was amplified by PCR, followed by restriction fragment length polymorphism analysis as described previously (Mencej et al. 2006). Genotyping of –290C>T (rs9525641) SNP was performed by the TaqMan allelic discrimination method, using the TaqMan SNP genotyping assay (Applied Biosystems, Foster City, CA, USA) under the following conditions: 10 min at 95 °C, followed by 43 cycles of 15 s at 95 °C and 1 min at 60 °C. Reactions were performed on ABI Prism 7000 sequence detection system.
BMD measurement
BMD measurements at the lumbar spine (L1–L4), total hip, and femoral neck were performed using dual-energy X-ray absorptiometry (DXA; QDR-4500, Hologic, Inc., Waltham, MA, USA) in Ljubljana and Celje. Cross-calibration study of the precision of measurements between the centres had previously been performed. A correction factor was not considered necessary.
Biochemical markers of bone turnover
Serum-free soluble RANKL (sRANKL) was measured by an enzyme immunoassay (sRANKL ELISA, Biomedica, Vienna, Austria). Serum C-terminal cross-linking telopeptides of type I collagen (CTX) were measured by enzyme immunoassay (Serum CrossLaps ELISA, Nordic Bioscience Diagnostics A/S, Herlev, Denmark). Osteocalcin (OC) in heparinised plasma was measured by a solid-phase, two-site chemiluminescent enzyme-labelled immunometric assay (Immulite Osteocalcin, Diagnostic Product Corporation, Los Angeles, CA, USA).
Luciferase reporter assay
The selected TNFSF11 gene promoter sequence spans nucleotides –1715 to +94 relative to the transcription start site. Human genomic DNA from two individuals with the two most common haplotypes, CCG and TTC, was amplified by PCR, using the primers: 5'-ATTATGGTACCGTCAAGGAGCAGGGAGAGAATG-3' and 5'-TCGGAGATCTCGCTTCGGAGCTCTCCTC-3'. Primers were designed according to the GenBank sequence AF544022 with the addition of the KpnI and BglII recognition sequence at the 5'-end. PCR was carried out using Taq polymerase (Qiagen) under the following conditions: 94 °C for 4 min, followed by 33 cycles of 94 °C for 30 s, 58 °C for 1 min and 72 °C for 2 min. The resulting PCR products were digested with KpnI and BglII (New England Biolabs, Beverly, MA, USA) and inserted into pGL3-basic luciferase reporter vector (Promega) digested at the same sites.
The HEK293 cells were maintained in
-minimal essential medium (
-MEM) supplemented with L-glutamine (1% v/v), fetal calf serum (10% v/v) and penicillin/streptomycin (1% v/v) (Gibco). For transfection, cells were plated at a density of 100 000 cells/well in 24-well tissue culture plates. After 24 h, the cells were transfected in six replicates with a mixture of Fugene HD transfection reagent (Roche Applied Science) at a ratio (µl reagent:µl DNA) of 3:2,
-MEM and 500 ng DNA/well, including 100 ng/well of pRP-TK control reporter vector. Cells were harvested 48 h post-transfection, and the luciferase assay performed using Dual-Luciferase Reporter Assay System (Promega). Luminescence was measured using a BIO-TEK Synergy HT multidetection microplate reader (Fisher Scientific, Pittsburgh, PA, USA).
Statistical analysis
Hardy–Weinberg equilibrium was tested for each SNP and group of participants using the
2 test. In each studied SNP, frequency distributions of genotypes in osteoporotic and non-osteoporotic postmenopausal women were compared using the
2 test. The standardized measures of LD, denoted as D' and r2, were assessed using the EMLD software (University of Texas, Houston, TX, USA), which calculates pairwise LD based on SNP genotype data from unrelated individuals (epi.mdanderson.org/
qhuang/Software/pub.htm). For two independent loci, the difference (D) is that between the actual and expected gametic frequencies and is usually expressed as a standardized difference (D'). D measures the statistical association of alleles in forming gametes and is related to the Pearson correlation coefficient (r). D' is used to assess the probability of historical recombination, whereas r2 is the most relevant measure in the association studies (Mueller 2004). PHASE program was used to estimate haplotypes from genotype data for individual participants (Stephens et al. 2001). For the purpose of statistical analysis, we used only haplotype data in which the probability of correct haplotype assignment by PHASE in an individual participant was estimated to be 95% or greater. Kolmogorov–Smirnov normality test was conducted before association analysis and data transformation was performed where appropriate. In haplotype analysis, subjects were coded according to whether they had two copies, one copy or no copies of the haplotype under study. Differences in BMD and biochemical markers between the genotype and haplotype subgroups were tested using ANOVA. For comparison of luciferase activities, the two-sample t-test was used. Statistical analyses were carried out using SPSS version 14.0 (SPSS Inc., Chicago, IL, USA) and P values <0.05 were considered statistically significant.
Possible TF binding sites at the sites of our nucleotide changes in TNFSF11 gene promoter were searched using TFSerach ver 1.3 (Heinemeyer et al. 1998).
| Results |
|---|
|
|
|---|
|
Evidence of strong LD was obtained within each of the three possible pairs: –290C>T and –643C>T, D'=0.96, r2=0.87; –290C>T and –693G>C, D'=0.93, r2=0.84; –643C>T and –693G>C, D'=0.92, r2=0.78.
In the group of postmenopausal women, the haplotypes, defined by –290C>T, –643C>T and –693G>C TNFSF11 gene promoter SNPs, were distributed as follows: TTC (49.3%), CTC (0.6%), TCC (2.6%), CCC (1.5%), TTG (1.2%), CTG (0.1%), TCG (0.4%) and CCG (44.3%). All eight possible haplotypes were found, but TTC and CCG accounted for 93.6% of the alleles observed.
To study whether different haplotypes have functional relevance, transient transfections with two different plasmids, containing these two most common haplotypes were performed for transcriptional activity evaluation. Forty-eight-hours post-transfection Renilla and firefly luciferase activities were measured and expressed as a ratio of firefly to Renilla activity. The ratio was 2.966 and 2.543 in the CCG and TTC haplotypes respectively, the difference between the ratios being statistically significant (P=0.018).
Analysis of TF binding motifs revealed numerous possible TF binding sites in the TNFSF11 gene promoter, two of which include two of the studied SNPs (Fig. 1). Polymorphism –290C>T lies in a potential TF binding motif for heat shock factors (95.3% homology) and –693G>C in ADR1 (87.7% homology).
|
|
|
| Discussion |
|---|
|
|
|---|
As the RANKL/RANK/OPG system plays a crucial role in osteoclast development, the TNFSF11 gene is certainly a candidate for the genetic regulation of bone mass and remodelling. In men, SNP rs9594782 has shown an association with BMD (Hsu et al. 2006). In Korean postmenopausal women, SNPs rs12721445 and rs2277438 in TNFSF11 gene have been studied and the latter was related with BMD, but only in combination with G1181C SNP in the OPG gene (Kim et al. 2007). In a study on Slovenian postmenopausal women, three polymorphic nucleotide variations (–290C>T, –643C>T and –693G>C) were studied in the TNFSF11 gene promoter. In –290C>T, an association was shown with BMD-fn in postmenopausal women (Mencej et al. 2006). Therefore, we focused our research on the SNPs –290C>T, –643C>T and –693G>C, which have already been studied in Caucasian postmenopausal women. All three studied SNPs in the TNFSF11 gene promoter region were found to be in strong LD. However, it is possible that the three studied SNPs are in LD with some other nearby sequence variation that is the actual functional cause of the association with BMD. To evaluate the functional role of the studied SNPs, luciferase reporter gene experiments were performed. Comparison of the two most common inferred haplotypes, CCG and TTC, showed a small but significant difference in their transcriptional activities, CCG being more active. Thus, higher serum RANKL concentrations are expected in the CCG haplotype.
To confirm the positive results from the functional analysis study, association analysis was performed. In 182 osteoporotic women, all studied SNPs in the TNFSF11 gene promoter were associated with BMD-ls, and one of them, –693G>C, also with BMD-fn. However, no association with BMD was found in 222 non-osteoporotic women. In the previous study on 115 postmenopausal women, only –290C>T was associated with BMD-fn. The present study was conducted on a larger number of participants and the homogeneity of groups was improved, as participants were involved according to their menopausal status and the presence of osteoporosis.
One of the focuses of the research was haplotype analysis. All eight possible haplotypes were inferred in our group. The two most common ones, TTC and CCG, showed an association with BMD-ls. Women carrying two alleles with TTC haplotype had higher BMD-ls than those without. The values of BMD-ls were similar for osteoporotic women with either one or two copies of the TTC allele. The CCG haplotype was also associated with BMD-ls, but BMD-ls values were lowest in individuals with two copies of CCG haplotype and highest in those with no CCG allele. The CCG and TTC haplotypes were not associated with BMD-fn, although women with two CCG alleles had lower BMD-fn.
Biochemical markers of bone turnover, which respond to changes more rapidly than BMD, did not show any association with the studied genotypes or haplotypes in any of the groups. Serum RANKL concentrations in osteoporotic women also showed no association with any of the studied genotypes or haplotypes. Higher transcriptional activity was observed in the CCG haplotype and, as expected, much higher RANKL concentrations were measured in 45 serum samples from osteoporotic women with two CCG alleles. As RANKL induces osteoclastogenesis, higher levels of RANKL could lead to increased bone loss and lower BMD. Indeed, osteoporotic women with two copies of CCG haplotype had lower BMD than those with two copies of TTC allele.
The transcriptional activity of the TNFSF11 gene could be influenced, if one or all SNPs were to lie in a TF binding site. So far, several putative TFs binding motifs with significant homology (>85%) have been identified in the human TNFSF11 gene promoter region including heat shock factor-responsive elements, CdxA, SRY, MZF-1, Nkx-2, NF-E2, AML-1a, C/EBP-β, ADR1, AP-4, GATA-1, GATA-2 and c-Myb (Roccisana et al. 2004). Our analysis showed that the two studied SNPs in the TNFSF11 gene promoter are located in possible binding sites for known TFs. Heat shock factors have already been functionally evaluated further upstream in human TNFSF11 gene promoter (Roccisana et al. 2004). However, no functional analysis of the ADR1 TF binding site has been performed on human TNFSF11 gene promoter.
In conclusion, the results of our study support the proposal that the TNFSF11 gene promoter region SNPs –290C>T, –643C>T and –693G>C have a functional role in osteoporosis development, due to the difference in transcription activity between two common haplotypes. Haplotype CCG had higher transcriptional activity and was associated with lower BMD-ls. Therefore, the studied haplotypes could be considered as risk factors for genetic susceptibility to postmenopausal osteoporosis, since a statistically significant association with BMD-ls has been demonstrated in all three, as well as in the two most common haplotype forms TTC and CCG. These findings were further supported by the association of –693G>C polymorphism with BMD-fn. These results provide important additional information on the functional and genetic role of TNFSF11 gene in bone density. However, the genetic effects of the TNFSF11 gene promoter SNPs need further confirmation on a larger group of subjects and the next step in the explanation of their functional role will require identification of the binding TF(s).
| Acknowledgements |
|---|
ak for their technical assistance. We thank Prof. Roger Pain for advice on the English language. We also thank the European Calcified Tissue Society for the support on functional analysis work. This work was supported by grant P3-0298 provided by Ministry of Education, Science and Sport of the Republic of Slovenia. The functional analysis work was supported by European Calcified Tissue Society Exchange Scholarship Grant. There is no conflict of interest that would prejudice its impartiality. | References |
|---|
|
|
|---|
Glass DA II, Bialek P, Ahn JD, Starbuck M, Patel MS, Clevers H, Taketo MM, Long F, McMahon AP, Lang RA et al. 2005 Canonical Wnt signaling in differentiated osteoblasts controls osteoclast differentiation. Developmental Cell 8 751–764.[CrossRef][Web of Science][Medline]
Heinemeyer T, Wingender E, Reuter I, Hermjakob H, Kel AE, Kel OV, Ignatieva EV, Ananko EA, Podkolodnaya OA, Kolpakov FA et al. 1998 Databases on transcriptional regulation: TRANSFAC, TRRD, and COMPEL. Nucleic Acids Research 26 364–370.
Hofbauer LC & Heufelder AE 2001 Role of receptor activator of nuclear factor-
B ligand and osteoprotegrin in bone cell biology. Journal of Molecular Medicine 79 243–253.[CrossRef][Web of Science][Medline]
Holmen SL, Zylstra CR, Mukherjee A, Sigler RE, Faugere MC, Bouxsein ML, Deng L, Clemens TL & Williams BO 2005 Essential role of beta-catenin in postnatal bone acquisition. Journal of Biological Chemistry 280 21162–21168.
Hsu Y, Niu T, Terwedow H, Xu X, Feng Y, Li Z, Brain JD, Rosen CJ, Laird N & Xu X 2006 Variation in genes involved in the RANKL/RANK/OPG bone remodelling pathway are associated with bone mineral density at different skeletal sites in men. Human Genetics 118 568–577.[CrossRef][Web of Science][Medline]
International HapMap Consortium A second generation human haplotype map of over 3.1 million SNPs Nature 449 2007 851–861.[CrossRef][Medline]
Kim JG, Kim JH, Kim JY, Ku SY, Jee BC, Suh CS, Kim SH & Choi YM 2007 Association between osteoprotegerin (OPG), receptor activator of nuclear factor-
B (RANK), and RANK ligand (RANKL) gene polymorphisms and circulating OPG, soluble RANKL levels, and bone mineral density in Korean postmenopausal women. Menopause 14 913–918.[CrossRef][Web of Science][Medline]
Knight JC 2003 Functional implications of genetic variation in non-coding DNA for disease susceptibility and gene regulation. Clinical Science 104 493–501.[CrossRef][Web of Science][Medline]
Lacey DL, Timms E, Tan HL, Kelley MJ, Dunstan CR, Burgess T, Elliott R, Colombero A, Elliott G, Scully S et al. 1998 Osteoprotegrin ligand is a cytokine that regulates osteoclast differentiation and activation. Cell 93 165–176.[CrossRef][Web of Science][Medline]
Mencej S, Prezelj J, Kocijancic A, Ostanek B & Marc J 2006 Association of TNFSF11 gene promoter polymorphisms with bone mineral density in postmenopausal women. Maturitas 55 219–226.[CrossRef][Web of Science][Medline]
Mueller JC 2004 Linkage disequilibrium for different scales and applications. Briefings in Bioinformatics 5 355–364.
Ralston SH & Crombrugghe B 2006 Genetic regulation of bone mass and susceptibility to osteoporosis. Genes and Development 20 2492–2506.
Roccisana JL, Kawanabe N, Kajiya H, Koide M, Roodman GD & Reddy SV 2004 Functional role for heat shock factors in the transcriptional regulation of human RANK ligand expression in stromal/osteoblast cells. Journal of Biological Chemistry 12 10500–10507.
Simonet WS, Lacey DL, Dunstan CR, Kelley M, Chang MS, Luthy R, Nguyen HQ, Wooden S, Bennett L, Boone T et al. 1997 Osteoprotegrin: a novel secreted protein involved in the regulation of bone density. Cell 89 309–319.[CrossRef][Web of Science][Medline]
Stephens M, Smith NJ & Donnelly P 2001 A new statistical method for haplotype reconstruction from population data. American Journal of Human Genetics 68 978–989.[CrossRef][Web of Science][Medline]
Troen BR 2003 Molecular mechanisms underlying osteoclast formation and activation. Experimental Gerontology 38 605–614.[CrossRef][Web of Science][Medline]
Walsh MC & Choi Y 2003 Biology of TRANCE axis. Cytokine and Growth Factor Reviews 14 251–263.[CrossRef][Web of Science][Medline]
World Health Organisation Assessment of fracture risk and its application to screening for postmenopausal osteoporosis: report of a WHO Study Group World Health Organization Technical Report Series 843 1994 1–129.[Medline]
Received in final form 13 March 2008
Accepted 28 March 2008
Made available online as an Accepted Preprint 28 March 2008
This article has been cited by other articles:
![]() |
S. Mencej-Bedrac, J. Prezelj, T. Kocjan, K. Teskac, B. Ostanek, M. Smelcer, and J. Marc The combinations of polymorphisms in vitamin D receptor, osteoprotegerin and tumour necrosis factor superfamily member 11 genes are associated with bone mineral density J. Mol. Endocrinol., March 1, 2009; 42(3): 239 - 247. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |