help button home button Endocrine Society JCEM JCEM Call for Nominations for EIC
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

Journal of Clinical Endocrinology & Metabolism, doi:10.1210/jc.2005-0364
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
90/8/4829    most recent
Author Manuscript (PDF)
Right arrow Submit a related Letter to the Editor
Right arrow Purchase Article
Right arrow View Shopping Cart
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow Request Copyright Permission
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Garnero, P.
Right arrow Articles by Delmas, P. D.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Garnero, P.
Right arrow Articles by Delmas, P. D.
Related Collections
Right arrow Calcium and Bone Metabolism
Right arrow Female Endocrinology
The Journal of Clinical Endocrinology & Metabolism Vol. 90, No. 8 4829-4835
Copyright © 2005 by The Endocrine Society

Vitamin D Receptor Gene Polymorphisms Are Associated with the Risk of Fractures in Postmenopausal Women, Independently of Bone Mineral Density

P. Garnero, F. Munoz, O. Borel, E. Sornay-Rendu and P. D. Delmas

Institut National de la Santé et de la Recherche Médicale Research Unit 403 (P.G., F.M., O.B., E.S.-R., P.D.D.), 69003 Lyon, France; and Synarc (P.G.), 69003 Lyon, France

Address all correspondence and requests for reprints to: Dr. Patrick Garnero, Institut National de la Santé et de la Recherche Médicale Unit 403, Hôpital E Herriot, Pavillon F, 69437 Lyon Cedex 03, France. E-mail: patrick.garnero{at}synarc.com.


    Abstract
 Top
 Abstract
 Introduction
 Subjects and Methods
 Results
 Discussion
 References
 
Context: Osteoporosis is a systemic disease with a strong genetic component. Vitamin D receptor (VDR) gene polymorphisms explain only a small part of the genetic influence on the level of bone mineral density (BMD), whereas their effect on fracture remains uncertain.

Objective: The objective of this study was to investigate the relationships between VDR genotypes and fracture risk

Design: A prospective population-based cohort was studied.

Subjects: A total of 589 postmenopausal women (mean age, 62 yr) were followed prospectively during a median (interquartile) of 11 (1.1) yr.

Main Outcome Measure: The study measured incidents of vertebral and nonvertebral fractures.

Results: VDR allele B was significantly and dose dependently overrepresented in women who fractured, including 34 and 86 women with first incident vertebral and nonvertebral fragility fractures, respectively. This corresponded to an odds ratio of 1.5 (95% confidence interval, 0.95–2.40) for heterozygous carriers (bB, n = 286) and 2.10 (95% confidence interval, 1.16–3.79) for homozygous carriers (BB, n = 90) of the B allele, compared with women with the bb genotype (n = 213). VDR genotype groups did not differ for demographics, physical activity, grip strength, personal and maternal history of fracture, and calcium intake. The association was independent of BMD of the spine, hip, and radius, and of the BMD loss at the radius. The relationship between VDR polymorphisms and fracture risk was not altered after adjustment for baseline circulating levels of bone turnover markers, estradiol, dehydroepiandrosterone sulfate, SHBG, IGF-I, intact PTH, and 25 hydroxyvitamin D.

Conclusion: VDR genotypes are associated with the risk of fracture in postmenopausal women independently of BMD, rate of postmenopausal forearm BMD loss, bone turnover, and endogenous hormones. The mechanisms by which VDR genotypes influence bone strength remain to be determined.


    Introduction
 Top
 Abstract
 Introduction
 Subjects and Methods
 Results
 Discussion
 References
 
OSTEOPOROSIS IS A systemic disease characterized by decreased bone mass and microarchitectural deterioration of bone tissue, leading to increased bone fragility and fracture. Although this disease is clearly multifactorial involving several environmental influences, genetic factors play a major role in its pathogenesis. Twin and pedigree studies have shown that genetic influences account for 50–80% of the interindividual variability of bone mineral density (BMD) in young adults (1, 2, 3). Conversely, the magnitude of the influence of genetic factors on the rate of postmenopausal bone loss in postmenopausal women remains controversial (4, 5), with one study finding a significant genetic effect for spine, but not hip, BMD loss (5), and another reporting no significant influence (4).

Because a low BMD in postmenopausal women, which results from a low peak BMD and or a faster rate of bone loss, is a strong determinant of fracture risk, most genetic studies in osteoporosis have been dedicated to identifying gene polymorphisms that could explain the interindividual variance in these two components. During the past 10 yr, since the initial report from Morrison et al. (6), several genetic polymorphisms in candidate genes for osteoporosis, such as vitamin D receptor (VDR), estrogen receptor, PTH receptor, calcitonin receptor, IL-1 receptor antagonist, IL-6, TGF-ß, IGF-I, type I collagen, osteocalcin, apolipoprotein E, methylenetetrahydrofolate reductase, peroxisome proliferators-activated receptor, osteoprotegerin, and GnRH have been investigated (for reviews see Refs.7 and 8). More recently, associations were reported between BMD and gene polymorphisms for lipoprotein receptor-related protein 5 (9, 10, 11), cathepsin K (12), and pituitary glutaminyl cyclase (13). However, their relation with BMD and/or postmenopausal bone loss was inconsistent between studies, and when present, the association was of low magnitude.

Among these candidate genes, the first and by far most extensively investigated is that coding for the VDR and, more specifically, the BsmI restriction fragment length polymorphism (RFLP). The VDR plays a role in regulating calcium homoeostasis through binding and nuclear translocation of 1{alpha}25(OH)2D3, affecting bone resorption, and increasing calcium absorption (14). The most recent metaanalysis of the influence of VDR on BMD and BMD change in women indicated a significant but very modest effect on spine BMD in postmenopausal women; individuals carrying the BB genotype have a 2.4% lower BMD than the Bb/bb genotype groups (15). No significant association could be detected with spine BMD in premenopausal women, with pre- and postmenopausal femoral neck BMD and with BMD change at the femoral neck (15). However, analyzing the relationship between BMD and VDR polymorphism may not be sufficient to assess the role of VDR in osteoporosis. Indeed, although a low BMD is a major determinant of increased fracture risk, about half of patients with incident fractures have baseline BMD above the WHO diagnostic threshold of osteoporosis (T score ≤ –2.5) (16, 17, 18, 19). Bone fragility also depends on the morphology, the architecture, and remodeling of bone as well as on the quality (properties) of the bone matrix. Some of these BMD-independent factors may be assessed to a certain extent by biochemical markers of bone turnover, bone ultrasound properties and hip geometry for which a genetic contribution has been reported (20, 21, 22, 23). Interestingly, a recent twin study performed in 6570 women showed that the genes involved in the genetic influence of wrist fracture are probably, for a vast majority, different from those implicated in the determination of forearm BMD (24). Thus, analyzing gene polymorphisms directly in relation to fracture risk may be more clinically relevant than studying the genetic influence on BMD.

A few studies have analyzed the associations between VDR BsmI genotypes and fracture. Most studies were cross-sectional, comparing the genotype distribution between women with prevalent vertebral or hip fracture and controls. They usually included a limited number of fracture cases and, like the BMD studies, have generated conflicting results (25, 26, 27, 28, 29). The association between VDR genotypes and the risk of incident fractures was investigated in three prospective case control or population studies and also yielded conflicting results (30, 31, 32). The potential mechanisms by which VDR genotypes may influence fracture risk, including bone turnover, rate of BMD changes, and circulating hormones were not investigated in these studies.

The aim of our study was to investigate the relationships between VDR BsmI polymorphisms and the risk of incident vertebral and nonvertebral fractures in postmenopausal women followed prospectively for up to 12 yr. We also analyzed the potential influence of BMD, rate of BMD change, bone turnover, circulating hormones, and clinical risk factors for fracture on this relationship.


    Subjects and Methods
 Top
 Abstract
 Introduction
 Subjects and Methods
 Results
 Discussion
 References
 
We investigated the postmenopausal women of the Os des Femmes de Lyon (OFELY) cohort, a prospective study of the determinants of bone loss that has been previously described in detail (33). The cohort of this study comprises 1039 Caucasian women, 31–89 yr of age, including 671 postmenopausal women, recruited between February 1992 and December 1993 from the regional section of a large health insurance company (Mutuelle Génerale de l’Education Nationale, Lyon, France) of the Rhône district (i.e. Lyon and its surroundings in France). Written informed consent was obtained from each woman, and the study was approved by the local ethical committee. We followed the women from this cohort during a median of 11 yr. In this analysis we investigated the 671 subjects who were menopausal at the baseline visit. DNA was not available for 82 women, resulting in a final group of 589 women.

Fracture assessment

Nonvertebral fractures. Prior fractures were those that occurred after the age of 40 yr, identified by self-reporting during the baseline questionnaire. Incident fractures were reported during each annual follow-up. For women who did not come to the clinical center, a letter was sent every year to identify the occurrence of fractures. All peripheral fractures were confirmed by radiographs or by a surgical report. Only low-trauma fractures, i.e. those occurring with falls from standing height or less, were taken into account, and we excluded fractures of fingers, toes, skull, and face.

Vertebral fractures. Lateral x-ray films of the thoracic and lumbar spine were obtained at baseline for 97% of women, and on average after 7.5 yr ± 1.6 (from 2.8–10.3) of follow-up for 86% of them. All prevalent and incident vertebral fractures were identified by the semiquantitative method of Genant et al. (34) by a trained rheumatologist. A vertebra was classified as fractured on the baseline radiograph if any vertical height (anterior, middle, and/or posterior) was reduced by more than 20%. An incident fracture was defined by a decrease of 20% or more and at least 4 mm in any vertebral height of one or more thoracic or lumbar vertebrae between follow-up and baseline x-ray films. We excluded vertebral fractures that occurred because of major trauma and vertebral deformities due to other causes than osteoporosis such as osteoarthritis and Scheuermann’s disease.

Clinical evaluation and physical examination

Women completed a questionnaire at the initial screening visit and at each annual follow-up as described previously (35). It included medical history, medication use such as hormone replacement therapy, calcium and vitamin D supplementation, calcium intake, reproductive characteristics, age of menopause, and family history of fragility fracture. The occurrence of falls during the past 12 months was recorded. Physical activity was expressed by a score calculated from sport or recreation and job and home activities as described previously (35). The grip strength was measured at baseline by a hand dynamometer (Vigorimeter; Martin, Tuttlingen, Germany) on the left and right hand, using the maximum of two readings for each hand.

Bone densitometry

BMD was measured at the lumbar spine, total hip, whole body (BMD and bone mineral content), and mid, distal, and ultradistal radius. The mid area is composed mainly (~95%) of cortical bone, the distal area comprises both cortical and trabecular bone (~25%), and the ultradistal contains more trabecular bone. Subsequent BMD measurements at the radius were performed at yr 2 and then annually. Forearm BMD was measured by pencil beam dual-energy x-ray absorptiometry (DXA) with a QDR 2000 device between 1992 and 1999 and with a QDR 1000+ between 1999 and 2003 (Hologic, Waltham, MA). A cross-calibration between the two scanners has been done. A control phantom was scanned every day, and all DXA measurements were performed by the same experienced operator. The in vitro long-term precision error of DXA over the 11-yr study was 0.41%.

A valid assessment of the rate of BMD change at the other skeletal sites (hip, spine, whole body), measured with a multidetector mode (fan beam) was not obtained in this cohort because of technical problems that specifically affected this mode. A drift in phantom values was observed over time that could not be reliably corrected.

Biochemistry

Blood samples were collected between 0800 and 0930 h after an overnight fast at baseline. Serum samples were stored frozen at –70 C until assayed.

Bone formation was evaluated with human RIA for serum-intact osteocalcin, bone-specific alkaline phosphatase, and intact procollagen type I N-terminal propeptide as previously described (36). The intra- and interassay coefficients of variation for all three assays are less than 10%. Bone resorption was evaluated by measuring serum ß isomerized C-terminal cross-linking telopeptide of type I collagen (ßCTX) and urinary N-terminal cross-linking telopeptide of type I collagen by ELISA (36). Intra- and interassay coefficients of variation are lower than 10% for all assays. Urinary {alpha} and ß isomerized CTX were also measured to derive the urinary {alpha} to ß CTX ratio (37).

Serum total estradiol, estrone, testosterone, and dehydroepiandrosterone (DHEA) sulfate were measured by tritiated RIA after diethylether extraction or dilution (for DHEA sulfate) as previously described (36). The total assay precision is below 10% for all assays, and the limit of detection of estradiol is 3 pg/ml. SHBG was measured by RIA, serum-intact PTH by immunochemiluminometric assay, and 25 hydroxyvitamin D [25 (OH) D] level was determined by a radio-binding assay as previously described (36). Serum IGF-1 was measured by RIA after acid ethanol extraction.

VDR genotyping

Genomic DNA was extracted from whole blood according to standard methods and used for PCR of the sequences flanking the BsmI, TaqI, and ApaI restriction sites with a Corbett FTS-1 Thermal sequencer (Corbett Research, Mortlake, New South Wales, Australia). The PCR products were then digested by BsmI (New England Biolabs, Inc., Beverly, MA) or ApaI or TaqI (Promega Inc., Madison, WI), and genotypes were analyzed by agarose gel electrophoresis as previously described (33).

Statistical analyses

Differences between genotype groups in baseline clinical, biochemical, and BMD parameters were assessed by ANOVA (continuous variables) or Pearson {chi}2 (categorical variables) after logarithmic transformation of bone marker and hormone levels. Radius BMD changes were calculated from the mean of the individual slopes divided by the intercept and were expressed in percentage per year. Odds ratios (OR) with 95% confidence interval were obtained by logistic regression analyses to estimate the relative risk of fractures by genotypes of the risk allele, with no copies of the risk allele as the reference group. First, we calculated crude ORs and then, secondly we adjusted for potential confounding factors (age, BMD, prior fractures, bone markers, hormones). The population attributable risk (PAR), which corresponds to the proportion of women with incident fracture in the population that is attributable to the risk factor assuming a causal relationship, was calculated as population attributable risk = p * (OR – 1)/[1 + p * (OR – 1)] where p is prevalence of the risk factor. All statistical analyses were performed on SAS software (version 8; SAS Institute, Cary, NC).


    Results
 Top
 Abstract
 Introduction
 Subjects and Methods
 Results
 Discussion
 References
 
Baseline characteristics

VDR genotype distribution was in Hardy-Weinberg equilibrium, with no overrepresentation of the BB genotypes that represented 15% of the population. The BsmI, TaqI, and ApaI genotypes were highly associated as previously reported (33). Indeed, the genotype TT was found with bb in 98% of cases, and genotype tt was associated with genotype BB in 97% of cases. Similarly, genotype aa was found with bb in 60% of women and genotype AA was found with BB in 99% of cases. As shown in Table 1Go, women in the three BsmI genotype groups were similar in age, time since menopause, height, weight, hormone replacement use, and calcium intake. There were also no significant differences in these parameters between TaqI and ApaI genotype groups (data not shown).


View this table:
[in this window]
[in a new window]
 
TABLE 1. Baseline characteristics of postmenopausal women in relation to VDR BsmI genotypes

 
During a median (interquartile) 11 (1.1) yr follow-up, first incident fractures were recorded in 120 women, including 34 women with vertebral fractures and 86 with nonvertebral fractures.

Association between VDR genotypes and incident fracture

When women were grouped according to VDR BsmI genotypes, all fragility fractures were significantly overrepresented in women carrying the B allele (Table 2Go). Allele B was dose dependently associated with increased risk of all fragility fractures; women with genotype bB had a relative risk of fracture (95% confidence interval) of 1.50 (0.95–2.40) compared with women in the bb genotype, and those carrying the BB genotype had a relative risk of fracture of 2.10 (1.16–3.79). When nonvertebral and vertebral fractures were analyzed separately, VDR genotypes were significantly associated with an increased risk of nonvertebral fractures, but not with incident vertebral fractures (Table 2Go). When all fractures, i.e. not only those resulting from low trauma, were analyzed (n = 136), VDR genotypes were also associated with fracture risk, women with bB and BB genotype having a relative risk of fracture (95% confidence interval) of 1.44 (0.83–3.15) and 2.48 (1.34–4.57), respectively compared with women with genotype bb.


View this table:
[in this window]
[in a new window]
 
TABLE 2. Number of women with incident fractures according to VDR BsmI genotypes

 
As expected from the high association between TaqI and BsmI polymorphisms, TaqI genotypes were also associated with fracture risk with 15.9, 21.5, and 27.6% of all incident fragility fractures occurring in genotype TT, Tt, and tt, respectively (P {chi}2 = 0.06). In contrast, there was no significant association between ApaI polymorphisms and incidence of fracture (P {chi}2 = 0.45).

We then performed haplotype analysis by combining genotypes for BsmI, ApaI, and TaqI polymorphisms. Among the different combinations, genotypes bbAATT (0.85%), BbAATT (0.17%), BBAATT (0.17%), BBAATt (0.34%), BbAaTT (0.51%), bbAaTt (0.85%), BBAatt (0.17%), and BbaaTT (0.17%) accounted each for less than 1% of the population. Thus, we excluded from this haplotype analysis the 19 women containing one of these genotypes. Table 3Go shows the distribution of fractures in women grouped according to the remaining common genotypes. Because there was an overrepresentation of fractures in women carrying the "BAt" haplotype, we grouped women according to carrier status for this VDR haplotype as homozygous carriers (BBAAtt) and heterozygous carriers (including genotypes BbAaTt and BbAATt) of the risk haplotype and compared with women not carrying the haplotype (including genotypes bbaaTT and bbAaTT). There was no significant difference in the variables shown on Table 1Go between women grouped according to these VDR genotypes (data not shown). Women homozygous for the risk haplotype had a relative risk (95% confidence interval) of all fragility fractures of 2.0 (1.1–3.7) compared with women not carrying the haplotype, whereas heterozygous carriers had a relative risk of 1.4 (0.8–2.4).


View this table:
[in this window]
[in a new window]
 
TABLE 3. Number of postmenopausal women with fracture according to VDR BsmI, ApaI, and TaqI genotypes

 
Because this haplotype analysis gave association with fracture risk that was similar to that obtained with the BsmI polymorphism only, in the subsequent analyses, only BsmI genotypes were investigated.

VDR genotypes and risk factors for fractures

We investigated whether the relationships between VDR genotypes and fracture risk could be mediated by other risk factors for fracture that we identified in the OFELY cohort.

As previously reported (33), VDR genotypes were not associated with BMD at any of the skeletal sites investigated, including the total hip (Table 4Go). We thus investigated differences between VDR genotype groups in non-BMD risk factors for fractures.


View this table:
[in this window]
[in a new window]
 
TABLE 4. Risk factors for fractures in postmenopausal women in relation to VDR BsmI genotypes

 
Clinical factors

Six BMD-independent risk factors for all fractures have been previously identified in the OFELY cohort (35) including older age, number of falls in the previous year, lower grip strength, low physical activity, maternal history of fracture, and personal history of fractures. As shown on Table 1Go, there was no significant difference in any of these clinical risk factors between the VDR genotype groups.

Skeletal-related factors

In the OFELY cohort, we previously identified increased bone turnover, and more specifically, increased serum CTX and bone-specific alkaline phosphatase, as BMD-independent factors for all fractures (36). There was no significant difference in any of the biochemical markers of bone turnover investigated between VDR genotype groups (Table 4Go). We also previously reported that an increased urinary {alpha}/ß CTX was associated with an increased risk of all fractures, independently of BMD and of bone turnover rate (37). In this study, we found no difference in this ratio between VDR genotype groups (Table 4Go). As shown on Table 4Go, there was no significant difference in the rate of BMD change at the mid, distal, and ultradistal radius between VDR genotype groups.

Endogenous hormones

Among the various hormones we measured in the OFELY cohort at baseline, decreased serum levels of estradiol, DHEA sulfate, and IGF-I and increased concentrations of SHBG and PTH were associated with increased risk of all fractures independently of BMD (36). As shown on Table 4Go, there was no significant difference in the levels of these hormones and VDR genotypes, except for slightly higher PTH values in the bB and BB genotypes. Serum calcium and 25 (OH) D levels were also similar in the three genotype groups.

Multivariable model for the prediction of fracture risk

In multivariable logistic regression, none of the risk factors investigated above, taken individually or in combination, substantially altered the association between VDR genotypes and the risk of all fractures (data not shown). Because most of the variables above were associated with BMD and/or age and only PTH was significantly associated with VDR genotypes, we then performed a forward multiple variable model including age, hip BMD, and prevalent fractures as key established independent risk factors, and PTH and VDR polymorphisms. As shown on Table 5Go, the risk of fracture associated with genotype BB was of similar magnitude to that of low BMD and prior fractures, whereas serum PTH did not contribute significantly to the predictive model. As shown on Fig. 1Go, among women with a total hip BMD T-score more than –2.5 and genotype bb, only 11% had incident fracture, whereas 24% of women in the same BMD category, but carrying the BB genotype, had fracture. The highest fracture rate was observed in women with a BMD T-score less than –2.5 and carrying the BB genotype.


View this table:
[in this window]
[in a new window]
 
TABLE 5. Multiple variable logistic regression model for the prediction of all fractures

 


View larger version (20K):
[in this window]
[in a new window]
 
FIG. 1. Percentage of incident fractures according to total hip BMD levels and VDR BsmI genotypes. Low BMD was defined as a BMD T-score less than –2.5. Figures in the bars represent the number of incident fractures on the total number of women. *, P < 0.05 for allele dose association of VDR genotypes with incident fractures.

 

    Discussion
 Top
 Abstract
 Introduction
 Subjects and Methods
 Results
 Discussion
 References
 
In our study, we found that VDR genotypes were strongly associated with an increased risk of fractures independent of BMD levels, rate of radius BMD loss, and clinical and biochemical variables that have previously been shown to be associated with fracture risk.

Although VDR genotypes have been suggested previously to be associated with BMD or postmenopausal BMD loss, most recent metaanalyses indicated that, if any, the influence is limited to the spine and is very modest (15). In our population of postmenopausal women, we confirmed that VDR genotypes were not associated with BMD level at any skeletal sites including the spine, the hip, the forearm, and the whole body. The power of previous studies to analyze the association between VDR genotypes and rate of BMD changes was somewhat impaired because of the imprecision in the estimate of individual BMD loss, resulting from the short follow up that did not exceed 5 yr (15, 38, 39). To our knowledge, our study is the longest that has examined associations between VDR genotypes and rate of BMD changes measured over 11 yr at a very precise site, the radius. Thus, our data strongly suggest that VDR genotypes do not influence postmenopausal bone loss at the radius.

We found an association between VDR genotype and nonvertebral fracture, but not vertebral fractures. That might be a statistical power issue as there were relatively few incident vertebral fractures (n = 34) compared with nonvertebral ones (n = 86). Conversely, the difference in prediction might be related to the mechanism underlying the association between VDR genotype and nonvertebral fracture risk that is unclear as this stage. Two previous prospective case-control studies in postmenopausal women have reported conflicting results on the association between VDR genotypes and fracture risk. In the Nurses’ Health Study including postmenopausal women from 43–69 yr of age, VDR BB genotype was associated with a 2-fold increased risk of hip, but not distal radius fractures, which was not dose dependent (30). Conversely in the larger but older (>65 yr of age) population of the Study of Osteoporotic Fractures (31), no association was found between VDR genotypes and the risk of hip, vertebral, nonhip, and nonvertebral fractures. Finally, in a random sample of 1004 women aged 55–80 yr from the Rotterdam study, it was reported a significant dose-dependent relationship between VDR gene haplotypes and the risk of incident nonvertebral fractures occurring during a mean 3.8-yr follow-up (32). This relationship was independent of age, weight, and BMD. In a longer 7-yr follow-up of the same population, the authors subsequently reported a borderline significant association with incident vertebral fractures, but in contrast to their first report, no significant relationship with incident nonvertebral fractures (40). However, it should be mentioned that in the Rotterdam study, bat haplotypes constructed from direct analysis of BsmI, ApaI, and TaqI RFLPs were analyzed, and it was suggested by the authors that this may have improved the sensitivity to find significant association, although this has not been formally demonstrated. In our study, combining VDR BsmI, ApaI, and TaqI polymorphisms gave very similar results as the BsmI genotype only, although this analysis was not based on direct haplotyping. As reported in previous studies (15, 30), we found that B was the risk allele whereas others (32) reported that women carrying the b allele were at risk for fracture. These discrepancies may be due to technical errors leading to potential difference in genotype determination—although this is unlikely—or differences in population characteristics including interactions with environmental factors. Haplotyping of multiple VDR polymorphisms combined with functional analyses is likely to help in understanding these differences.

Our findings strongly suggest that VDR genotypes act directly and/or indirectly on bone fragility by mechanisms that are not reflected by the measurements of BMD or bone turnover. The originality and major strength of our study relates to the extensive clinical and biochemical investigation we performed on these women, which allowed us to adjust our analyses for several risk factors for fracture. Among the clinical parameters, grip strength and the number of falls in the previous years have been shown to be predictive of fractures independently of BMD in the OFELY cohort (35) and may also be associated with VDR polymorphisms as vitamin D regulates neuromuscular functions (41). Physical activity and calcium intake have also been suggested to interact with VDR genotypes on their effects on BMD or BMD changes (39, 42). We found no significant difference between genotypes in these physical, dietary, and lifestyle factors, and adjustment for these parameters did not weaken the association between VDR genotypes and fractures, suggesting that they are not likely to be involved in mediating the observed effect.

One of the major BMD-independent risk factors that has emerged in the recent years is an increased bone turnover which can be detected by sensitive biochemical markers. Among circulating hormones, low 17-ß estradiol, high SHBG, low DHEA sulfate, or low IGF-I have been shown to be associated with increased risk of hip, vertebral, and all fractures (36, 44) and genetic studies have reported a significant interaction between VDR and estrogen receptor {alpha} polymorphisms on BMD (45) or fracture risk (40). In our study, there was no significant association between VDR genotypes, biochemical markers of bone turnover, and the circulating levels of all these hormones. In addition, their inclusion in a multivariable model did not modify the effect of VDR polymorphisms on fracture risk.

An obvious potential mechanism by which VDR genotypes may influence bone fragility is through the calcium regulating system and the close relationships between the effects of vitamin D and PTH on bone metabolism. Although VDR alleles have been shown to be associated with differences in intestinal calcium absorption (46) and bone metabolism response to 1.25 (OH)2 D3 (47), the relationships between VDR genotypes and circulating levels of PTH, 25 (OH) D, or calcium have also been inconclusive both in healthy postmenopausal women and hemodialysis patients (48). In our study, we found significantly higher levels of serum-intact PTH in patients carrying the B allele, but no differences in serum 25 (OH) D and total calcium, in agreement with some, but not all studies (48). However, the difference between genotypes was very small, not dose dependent, and adjustment for circulating PTH levels did not modify the relationships between VDR genotypes and fracture risk. Consequently the biological relevance of this association between VDR genotypes and serum-intact PTH is unclear and could have resulted from chance findings.

Our study has strengths and some limitations. This is one of the largest and the longest prospective study on the relationships between VDR genotypes and fracture risk in postmenopausal women. As mentioned above, we were able to investigate most of the risk factors that could have confounded the association between VDR genotypes and fractures risk. However, we did not analyze the potential influence of VDR on bone geometry and bone matrix quality, two important determinants of skeletal fragility that may be genetically determined. It has been suggested that the effect of VDR on fracture risk may in part be mediated by its interaction with collagen I gene polymorphisms (32), which influence bone fragility mainly by changes in bone matrix quality rather than on BMD (49). However, we did not find any relationship between VDR genotypes and the urinary ratio of native to isomerized CTX, which has been proposed as an index of bone collagen matrix quality (37). In addition, no relationships between VDR genotypes and hip axis length or ultrasound parameters could be found in a twin study (50). We could not assess rate of BMD changes at the spine or the hip. We measured circulating levels of hormones that may not adequately reflect their local bone environment concentrations. We performed a single measurement of markers of bone turnover and hormones, which may not accurately reflect the bone turnover and hormone status because of the variability in the measurements. We did not perform direct haplotyping of BsmI, ApaI, and TaqI VDR RFLPs. The VDR polymorphisms we have investigated are probably nonfunctional and thus are unlikely to be directly involved in their association with fracture risk. However, they could be in linkage disequilibrium with one or more functional polymorphisms yet to be identified elsewhere in the VDR gene, which may explain the relationships we observed (43). Finally, VDR genotypes were not available for 82 women of the initial cohort. The baseline demographic characteristics of these women were similar and not significantly different from the rest of the cohort, except for an older age (66 vs. 62 yr), and consequently a slightly lower physical activity and higher prevalence of prior fracture. However, because of the absence of relationships between VDR genotypes and these parameters, this is not likely to have influenced our findings.

In summary, we found that VDR gene polymorphisms were dose dependently associated with the risk of all fractures independently of BMD, rate of postmenopausal BMD loss, and circulating levels of several biochemical variables. The biological mechanisms by which VDR genotypes influence bone strength remain to be determined.


    Footnotes
 
This work was partly supported by a contract from Institut National de la Santé et de la Recherche Médicale/Merck, Sharp & Dohme Chibret (the OFELY study).

First Published Online May 10, 2005

Abbreviations: BMD, Bone mineral density; CTX, C-terminal cross-linking telopeptide of type I collagen; DHEA, dehydroepiandrosterone; DXA, dual-energy x-ray absorptiometry; OFELY, Os des Femmes de Lyon; 25 (OH) D, 25 hydroxyvitamin D; OR, odds ratio; RFLP, restriction fragment length polymorphism; VDR, vitamin D receptor.

Received February 18, 2005.

Accepted April 28, 2005.


    References
 Top
 Abstract
 Introduction
 Subjects and Methods
 Results
 Discussion
 References
 

  1. Dequecker J, Nijs J, Verstraeten A, Geusens P, Gevers G 1987 Genetic determinants of bone mineral content at the spine and radius: a twin study. Bone 8:207–209[Medline]
  2. Pocock NA, Eisman JA, Hopper JL, Yeates PN, Sambrook PN, Ebert S 1987 Genetic determinants of bone mass in adults: a twin study. J Clin Invest 80:706–710
  3. Slemenda CW, Christian JC, Williams CJ, Norton JA, Johnston Jr CC 1991 Genetic determinants of bone mass in adult women: a reevaluation of the twin model and potential importance of gene interaction on heritability estimates. J Bone Miner Res 6:561–567[Medline]
  4. Christian JC, Yu PL, Slemenda CW, Johnston Jr CC 1989 Heritability of bone mass: a longitudinal study in aging male twins. Am J Hum Genet 44:429–433[Medline]
  5. Kelly PJ, Nguyen T, Hopper J, Pocock N, Sambrook P, Eisman J 1993 Changes in axial bone density with age: a twin study. J Bone Miner Res 8:11–17[Medline]
  6. Morrison NA, Qi JC, Tokita A Kelly PJ, Crofts L, Nguyen TV, Sambrook PN, Eisman JA 1994 Prediction of bone density from vitamin D receptor alleles. Nature 367:284–287[CrossRef][Medline]
  7. Hobson EE, Ralston SH 2001 Role of genetic factors in the pathophysiology and management of osteoporosis. Clin Endocrinol (Oxf) 54:1–9[CrossRef][Medline]
  8. Baldock PA, Eisman JA 2004 Genetic determinants of bone mass. Curr Opin Rheumatol 16:450–456[CrossRef][Medline]
  9. Ferrari SL, Deutsch S, Chodhury U, Chevalley T, Bonjour J-P, Dermitzakis T, Rizzoli R, Antonnonarakis E 2004 Polymorphisms in the low density-lipoprotein receptor-related protein 5 (LRP5) gene are associated with variation with vertebral bone mass, vertebral bone size, and stature in whites. Am J Hum Genet 74:866–875[CrossRef][Medline]
  10. Koay MA, Woon PY, Zhang Y, Miles LJ, Duncan EL, Ralston SH, Compston JE, Cooper C, Keen R, Langdahl BL, MacLelland A, O’Riordan J, Pols HA, Reid DM, Uitterlinden AG, Wass JAH, Brown MA 2004 Influence of LRP5 polymorphisms on normal BMD variation. J Bone Miner Res 19:1619–1627[CrossRef][Medline]
  11. Koller DL, Ichikawa S, Johnso ML, Lai D, Xuei X, Edenberg HJ, Conneally PM, Hui SL, Johnston CC, Peacock M, Foround T, Econs MJ 2005 Contribution of the LRP5 gene to normal variation in peak BMD in women. J Bone Miner Res 20:75–80[CrossRef][Medline]
  12. Giraudeau FS, McGinnis RE, Gray IC, O’Brien EJ, Doncaster KE, Spurr NK, Ralston SH, Reid DM, Wood H 2004 Characterization of common genetic variants in cathepsin K and testing for association with bone mineral density in a large cohort of perimenopausal women from Scotland. J Bone Miner Res 10:31–41
  13. Ezura Y, Katjita M, Ishida R, Yoshida S, Yoshida H, Suzuki T, Hosoi T, Inoue S, Shiraki M, Orimo H, Emi M 2004 Association of multiple nucleotide variations in the pituitary glutaminyl cyclase gene (QPCT) with low radial BMD in adult women. J Bone Miner Res 19:1296–1301[CrossRef][Medline]
  14. Pike JW, Yamamoto H, Shevde NK 2002 Vitamin D receptor-mediated gene regulation mechanisms and current concepts of vitamin D analog selectivity. Adv Ren Replace Ther 9:168–174[CrossRef][Medline]
  15. Thakkinstian A, D’Este C, Eisman J, Nguyen T, Attia J 2004 Meta-analysis of molecular association studies: vitamin D receptor gene polymorphisms and BMD as a case study. J Bone Miner Res 19:419–428[CrossRef][Medline]
  16. Siris ES, Miller PD, Barrett-Connor E, Faulkner KG, Wehren LE, Abbott TA, Berger ML, Santora AC, Sherwood LM 2001 Identification and fracture outcomes of undiagnosed low bone mineral density in postmenopausal women: results from the National Osteoporosis Risk Assessment. JAMA 286:2815–2822[Abstract/Free Full Text]
  17. Leslie WD, Metge C, Ward L 2003 Contribution of clinical risk factors to bone density-based absolute fracture risk assessment in postmenopausal women. Osteoporos Int 14:334–338[CrossRef][Medline]
  18. Miller PD, Siris ES, Barrett-Connor E, Faulkner KG, Wehren LE, Abbott TA, Chen YT, Berger ML, Santora AC, Sherwood LM 2002 Prediction of fracture risk in postmenopausal white women with peripheral bone densitometry: evidence from the National Osteoporosis Risk Assessment. J Bone Miner Res 17:2222–2230[CrossRef][Medline]
  19. Schuit SC E, van der Klift M, Weel AEAM, de Laet CEDH, Burger H, Seeman E, Hofman A, Uitterlinden AG, van Leeuwen JPTM, Pols HAP 2004 Fracture incidence and association with bone mineral density in elderly men and women: the Rotterdam Study. Bone 34:195–202[Medline]
  20. Kelly PJ, Hopper JL, Macaskill GT, Pocock NA, Sambrook PN, Eisman JA 1991 Genetic factors in bone turnover. J Clin Endocrinol Metab 72:808–813[Abstract]
  21. Tokita A, Kelly PJ, Nguyen TV, Qi JC, Morrison NA, Risteli L, Risteli J, Sambrook PN, Eisman JA 1994 Genetic influences on type I collagen synthesis and degradation: further evidence for genetic regulation of bone turnover. J Clin Endocrinol Metab 78:1461–1466[Abstract]
  22. Garnero P, Arden NK, Griffiths G, Delmas PD, Spector TD 1996 Genetic influence on bone turnover in postmenopausal twins. J Clin Endocrinol Metab 81:140–146[Abstract]
  23. Arden NK, Baker J, Hogg C, Baan K, Spector TD 1996 The heritability of bone mineral density, ultrasound of the calcaneus and hip axis length: a study of postmenopausal twins. J Bone Miner Res 11:530–534[Medline]
  24. Andrew T, Antionades L, Scurrah KJ, McGregor AJ, Spector TD 2005 Risk of wrist fracture in women is heritable and is influenced by genes that are largely independent of those influencing BMD. J Bone Miner Res 20:67–74[CrossRef][Medline]
  25. Looney JE, Yoon HK, Fischer M, Farley SM, Farley JR, Wergedal JE, Baylink DJ 1995 Lack of a high prevalence of the BB vitamin D receptor genotype in severely osteoporotic women. J Clin Endocrinol Metab 80:2158–2162[Abstract]
  26. Houston LA, Grant SFA, Reid DM, Ralston SH 1996 Vitamin D receptor polymorphism, bone mineral density, and osteoporotic vertebral fracture: studies in a UK population. Bone 18:249–252[Medline]
  27. Gomez C, Naves ML, Barrios Y, Diaz JB, Fernandez JL, Salido E, Torres A, Cannata JB 1999 Vitamin D receptor gene polymorphisms, bone mass, bone loss and prevalence of vertebral fracture: differences in postmenopausal women and men. Osteoporos Int 10:175–182[CrossRef][Medline]
  28. Aerssens J, Dequeker J, Peeters J, Breemans S, Broos P, Boonen S 2000 Polymorphisms of the VDR, ER and COLIA1 genes and osteoporotic hip fracture in elderly postmenopausal women. Osteoporos Int 11:583–591[CrossRef][Medline]
  29. Langdahl BL, Gravholt CH, Brixen K, Eriksen EF 2000 Polymorphisms in the vitamin D receptor gene and bone mass, bone turnover and osteoporotic fractures. Eur J Clin Invest 30:608–617[CrossRef][Medline]
  30. Feskanich D, Hunter DJ, Willett WC, Hankinson SE, Hollis BW, Hough HL, Kelsey KT, Colditz GA 1998 Vitamin D receptor genotype and the risk of bone fractures in women. Epidemiology 9:535–539[Medline]
  31. Ensrud KE, Stone K, Cauley JA, White C, Zmuda JM, Nguyen TV, Eisman JA, Cummings SR 1999 Vitamin D receptor gene polymorphisms and the risk of fractures in older women. For the Study of Osteoporotic Fractures Research Group. J Bone Miner Res 14:1637–1645[CrossRef][Medline]
  32. Uitterlinden AG, Weel AEAM, Burger H, Fang Y, Van Duijn CM Hofman A, Van Leeuwen JPTM, Pols HAP 2001 Interaction between the vitamin D receptor gene and collagen type I{alpha}1 gene in susceptibility for fracture. J Bone Miner Res 16:379–385[CrossRef][Medline]
  33. Garnero P, Borel O, Sornay-Rendu E, Delmas PD 1996 Vitamin D receptor gene polymorphisms are not related to bone turnover, rate of bone loss and bone mass in postmenopausal women: the OFELY study. J Bone Miner Res 11:827–834[Medline]
  34. Genant HK, Wu CY, Van Kuijk C, Nevitt MC 1993 Vertebral fracture assessment using a semi quantitative technique. J Bone Miner Res 8:1137–1148[Medline]
  35. Albrand G, Munoz F, Sornay-Rendu E, Duboeuf F, Delmas PD 2003 Independent predictors of all osteoporosis related fractures in healthy postmenopausal women. The OFELY study. Bone 32:78–91[Medline]
  36. Garnero P, Sornay-Rendu E, Claustrat B, Delmas PD 2000 Biochemical markers of bone turnover, endogenous hormones and the risk of fracture in postmenopausal women: the OFELY study. J Bone Miner Res 15:1526–1536[CrossRef][Medline]
  37. Garnero P, Cloos P, Sornay-Rendu E, Qvist P, Delmas PD 2002 Type I collagen racemization and isomerization and the risk of fracture in postmenopausal women: the OFELY study. J Bone Miner Res 17:826–833[CrossRef][Medline]
  38. Hansen TS, Abrahamsen B, Henriksen FL, Hermann AP, Jensen LB, Horder M, Gram J 1998 Vitamin D receptor alleles do not predict bone mineral density or bone loss in Danish perimenopausal women. Bone 22:571–575[Medline]
  39. Krall EA, Parry P, Lichter JB, Dawson-Hughes B 1995 Vitamin D receptor alleles and rates of bone loss: influences of years since menopause and calcium intake. J Bone Miner Res 10:978–984[Medline]
  40. Colin EM, Uitterlinden AJ, Meurs JBJ, Bergink AP, van de Klift M, Fang Y, Arp PP, Hofman A, Van Leeuwen JPTM, Pols HAP 2003 Interaction between vitamin D receptor genotype and estrogen receptor a genotype influences vertebral fracture risk. J Clin Endocrinol Metab 88:3777–3784[Abstract/Free Full Text]
  41. Pfeifer M, Begerow B, Minne HW 2002 Vitamin D and muscle function. Osteoporos Int 13:187–194[CrossRef][Medline]
  42. Ferrari SL, Rizzoli R, Slosman DO, Bonjour JP 1998 Do dietary calcium and age explain the controversy surrounding the relationship between bone mineral density and vitamin D receptor gene polymorphisms? J Bone Miner Res 13:363–370[CrossRef][Medline]
  43. Uitterlinden AG, Fang Y, Van Meurs JB, Pols HA, Van Leeuwen JP 2004 Genetics and biology of vitamin D receptor polymorphisms. Gene 338:143–156[CrossRef][Medline]
  44. Cummings SR, Browner WS, Bauer D, Stone K, Ensrud K, Jamal S, Ettinger B 1998 Endogenous hormones and the risk of hip and vertebral fractures among older women. N Engl J Med 339:733–738
  45. Willing M, Sowers M, Aron D, Clark MK, Burns T, Bunten C, Crutchfield M, D’Agostino D, Jannausch M 1998 Bone mineral density and its change in white women: estrogen and vitamin D receptor genotypes and their interaction. J Bone Miner Res 13:695–705[CrossRef][Medline]
  46. Gennari L, Becherini L, Masi L, Gonnelli S, Cepollaro C, Martini S, Mansani R, Brandi ML 1997 Vitamin D receptor genotypes and intestinal calcium absorption in postmenopausal women. J Clin Endocrinol Metab 82:1772–1775[Abstract/Free Full Text]
  47. Howard G, Nguyen T, Morrison N, Watanabe T, Sambrook P, Eisman J, Kelly PJ 1995 Genetic influences on bone density: physiological correlates of vitamin D receptor gene alleles in premenopausal women. J Clin Endocrinol Metab 80:2800–2805[Abstract]
  48. Bover J, Bosh RJ 1999 Vitamin D receptor polymorphisms as a determinant of bone mass and PTH secretion: from facts to controversies. Nephrol Dial Transplant 14:1066–1068[Free Full Text]
  49. Mann V, Hobson EE, Li B, Stewart TL, Grant SFA, Robins SP, Aspden RM, Ralston SHA 2001 COL1A1 Sp1 binding polymorphism predisposes to osteoporotic fracture by affecting bone density and quality. J Clin Invest 197:899–907
  50. Arden NK, Keen RW, Lanchbury JS, Spector TD 1996 Polymorphisms of the vitamin D receptor gene do not predict quantitative ultrasound of the calcaneus or hip axis length. Osteoporos Int 6:334–337[CrossRef][Medline]



This article has been cited by other articles:


Home page
Endocr. Rev.Home page
P. Chavassieux, E. Seeman, and P. D. Delmas
Insights into Material and Structural Basis of Bone Fragility from Diseases Associated with Fractures: How Determinants of the Biomechanical Properties of Bone Are Compromised by Disease
Endocr. Rev., April 1, 2007; 28(2): 151 - 164.
[Abstract] [Full Text] [PDF]


Home page
J. Clin. Endocrinol. Metab.Home page
Y. Fang, J. B. J. van Meurs, F. Rivadeneira, N. M. van Schoor, J. P. T. van Leeuwen, P. Lips, H. A. P. Pols, and A. G. Uitterlinden
Vitamin D Receptor Gene Haplotype Is Associated with Body Height and Bone Size
J. Clin. Endocrinol. Metab., April 1, 2007; 92(4): 1491 - 1501.
[Abstract] [Full Text] [PDF]


Home page
IBMS BoneKEyHome page
S. Ferrari
Single Gene Mutations and Variations Affecting Bone Turnover and Strength: a Selective 2006 Update
IBMS BoneKEy, December 1, 2006; 3(12): 11 - 29.
[Abstract] [Full Text] [PDF]


Home page
J. Clin. Endocrinol. Metab.Home page
T. V. Nguyen, L. M. Esteban, C. P. White, S. F. Grant, J. R. Center, E. M. Gardiner, and J. A. Eisman
Contribution of the Collagen I {alpha}1 and Vitamin D Receptor Genes to the Risk of Hip Fracture in Elderly Women
J. Clin. Endocrinol. Metab., December 1, 2005; 90(12): 6575 - 6579.
[Abstract] [Full Text] [PDF]


Home page
Hum Mol GenetHome page
A. d'Alesio, M. Garabedian, J. P. Sabatier, G. Guaydier-Souquieres, C. Marcelli, A. Lemacon, O. Walrant-Debray, and F. Jehan
Two single-nucleotide polymorphisms in the human vitamin D receptor promoter change protein-DNA complex formation and are associated with height and vitamin D status in adolescent girls
Hum. Mol. Genet., November 15, 2005; 14(22): 3539 - 3548.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
90/8/4829    most recent
Author Manuscript (PDF)
Right arrow Submit a related Letter to the Editor
Right arrow Purchase Article
Right arrow View Shopping Cart
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow Request Copyright Permission
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Garnero, P.
Right arrow Articles by Delmas, P. D.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Garnero, P.
Right arrow Articles by Delmas, P. D.
Related Collections
Right arrow Calcium and Bone Metabolism
Right arrow Female Endocrinology


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Endocrinology Endocrine Reviews J. Clin. End. & Metab.
Molecular Endocrinology Recent Prog. Horm. Res. All Endocrine Journals