Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
J Bone Miner Res. Author manuscript; available in PMC 2013 April 1.
Published in final edited form as:
PMCID: PMC3375687

Heritability of Prevalent Vertebral Fracture and Volumetric Bone Mineral Density and Geometry at the Lumbar Spine in Three Generations of the Framingham Study


Genetic factors likely contribute to the risk for vertebral fractures; however, there are few studies on the genetic contributions to vertebral fracture (VFrx), vertebral volumetric bone mineral density (vBMD) and geometry. Also the heritability (h2) for VFrx and its genetic correlation with phenotypes contributing to VFrx risk have not been established. This study aims to estimate the h2 of vertebral fracture, vBMD and cross-sectional-area (CSA) derived from quantitative computed tomography (QCT) scans, and to estimate the extent to which they share common genetic association in adults of European ancestry from three generations of Framingham Heart Study (FHS) families. Members of the FHS families were assessed for VFrx by lateral radiographs or QCT lateral scout views at 13 vertebral levels (T4-L4) using Genant’s semi-quantitative (SQ) scale (grades 0–3). Vertebral fracture was defined as having at least 25% reduction in height of any vertebra. We also analyzed QCT scans at the L3 level for integral (In.BMD) and trabecular (Tb.BMD) vBMD and cross-sectional area (CSA). Heritability estimates were calculated, and bivariate genetic correlation analysis was performed, adjusting for various covariates. For VFrx, we analyzed 4,099 individuals (148 VFrx cases) including 2,082 women and 2,017 men from 3 generations. Estimates of crude and multivariable-adjusted h2 were 0.43 to 0.69 (P< 1.1×10−2). 3,333 individuals including 1,737 men and 1,596 women from 2 generations had VFrx status and QCT-derived vBMD and CSA information. Estimates of crude and multivariable-adjusted h2 for vBMD and CSA ranged from 0.27 to 0.51. In a bivariate analysis, there was a moderate genetic correlation between VFrx and multivariable-adjusted In.BMD (−0.22) and Tb.BMD (−0.29). Our study suggests vertebral fracture, vertebral vBMD and CSA in adults of European ancestry are heritable, underscoring the importance of further work to identify the specific variants underlying genetic susceptibility to vertebral fracture, bone density and geometry.

Keywords: vertebral fracture, bone mineral density, heritability, QCT


Vertebral fractures (VFrx), comprise almost half of the 1.5 million new osteoporotic fractures annually in the United States (1). Prevalent vertebral fracture increases the risk of future fracture at the spine and other skeletal sites (2) and increases the risk of death in older adults (3). Studies have shown that both volumetric bone density (vBMD, g/cm3) and vertebral cross-sectional area (CSA) are important contributors to bone strength and VFrx risk (4), and low QCT-assessed vBMD has been associated with VFrx (5).

Whereas DXA-measured areal BMD (aBMD: g/cm2) is widely used clinically to assess fracture risk and is considered the “gold standard” for the diagnosis of osteoporosis, the conventional 2-D DXA does not provide volumetric bone density, and can be influenced by artifacts such as aortic calcification and degenerative changes in the spine. In contrast, quantitative computed tomography (QCT) provides 3-D measures of vBMD that offers the opportunity to examine specific compartments such as trabecular bone and to better estimate cross-sectional geometry (6,7).

While VFrx, vBMD and CSA are affected by many environmental factors, genetic factors also have been implicated. However, to our knowledge, there have been no previously published estimate of heritability for VFrx. The objective of this study was to estimate the heritability of VFrx, vBMD and CSA and to explore the extent to which VFrx and volumetric bone measurements share common genetic components.


Study Samples

Participants in this study include members of the original, 2nd, and 3rd generation cohorts in the Framingham Heart Study (FHS). The analysis of prevalent vertebral fracture includes 684 members of the original cohort who obtained lateral spine radiographs in years 1992–1993 (8) and 3448 members of the 2nd and 3rd generation cohorts who underwent CT examinations including lateral scout views in years 2002–2005 (9). The vBMD and CSA measurements were derived from the same QCT images in members of the 2nd and 3rd generation cohorts. All participants signed informed consent before enrollment and clinical examinations, and the study was approved by the Institutional Review Boards at Boston University Medical Center and Hebrew Rehabilitation Center.

Semi-Quantitative (SQ) Measurements for Vertebral Fracture

Prevalent vertebral fracture was assessed from T4 to L4 using the semi-quantitative grading system (10). In the original cohort, lateral radiographs (8) were used and in the 2nd and 3rd generation cohorts, scout images (9) were independently evaluated by trained radiologists. The SQ readings at each vertebral level were graded as 0 (no fracture), 1 (mild fracture: 20–25% reduction in any vertebral height), 2 (moderate fracture: 25–40% reduction) and 3 (severe fracture: > 40% reduction). In this study, prevalent fracture was defined as SQ grade ≥ 2 at any level, although we also repeated the analyses considering grade 1 and higher vertebral fractures. We also repeated the vertebral fracture heritability analyses in the subset of individuals who had both vertebral fracture assessments and vBMD measures.

Volumetric Bone Density and Vertebral Cross-sectional Area

QCT images of the lumbar spine obtained for 2nd and 3rd generation cohorts were analyzed using previously published software (11). The integral vBMD (In.BMD, g/cm3) was measured for the entire vertebral body including both cortical and trabecular compartments. Trabecular vBMD (Tb.BMD, g/cm3) was determined within an elliptical region inserted into the anterior vertebral body. CSA (cm2) was calculated from a central 10 mm thick slice perpendicular to the axis of the vertebral body. Note that although lumbar levels 2 to 4 were targeted for analysis on the QCT image, in many instances the L2 level was not included due to technical issues. Also, we excluded images due to poor quality or severe scoliosis. Lumbar level 3 (L3) was present in the highest proportion of participants (around 98%) such that for analyses of vBMD and CSA we only present results for L3.

Statistical Analysis

Descriptive statistics were generated using SAS v9.1 (SAS, Inc, Cary, NC). We used SOLAR statistical software to estimate heritability (h2) (12). We calculated h2 estimates for VFrx, vBMD and CSA, after adjustment for generation (model 1), for age, age2 and gender (model 2), and additionally for height (model 3). Age2 was included in the model based on the observation of non-linear relationship between prevalence of VFrx and age in our data. Using model 2 covariates, we also estimated h2 of VFrx additionally adjusting for either CSA, In.BMD, or Tb.BMD (Models 4a–4c). Bivariate analyses were also conducted in SOLAR to examine the genetic correlation (RhoG) among QCT-derived measurements of vBMDs, CSA and VFrx, adjusting for generation, age, age2 and gender.


Sample Characteristics

The characteristics of the sample are presented in Table 1. We analyzed data from 4,099 individuals including 2,082 women and 2,017 men for VFrx analysis while there were only 3,333 individuals including 1,737 men and 1,596 women available for vBMD analysis. Among 148 VFrx cases with SQ grade ≥ 2, 48 were present in men and 100 in women, with ages ranging from 34 to 96 years. The Original Cohort participants (age: 72 ~ 96) had the highest proportion of VFrx cases, 13.5% vs. less than 3% in Generation 2 (age: 40~87) and less than 1% in Generation3 (age: 31~72).

Table 1
Characteristics of participants from three generations of the Framingham Heart Study

Heritability Estimates and Bivariate Correlation Analysis

Both volumetric BMD measures (In.BMD, Tb.BMD) and CSA showed significant and strong h2 in crude analyses and in models adjusting for covariates including age, age2, gender and with/without height. The h2 in covariate-adjusted models consistently remained around 0.41 (SE = 0.04) and 0.43 (0.04) for Tb.BMD and In.BMD, respectively. Depending on covariates included, h2 for CSA ranged from 0.27 to 0.43 (Table 2). The h2 for VFrx with grade ≥2 was 0.43 to 0.53, similar to the values observed for vBMD in all models. The h2 for VFrx increased (approximately 0.69) after additional adjustment for either vBMD or CSA measurement. Using only individuals who were included in the vBMD analyses, the heritability was similar (approximately 70%) as observed using the larger sample. When we repeated the analysis using only VFrx with grade ≥1, the heritability was much lower (19% to 27%) (intra-reader agreement for grade ≥1, kappa =0.56 ~0.59, was not as good as agreement for grade ≥2, kappa=0.68~0.72). Bivariate analyses showed a negative genetic correlation between In.BMD and VFrx (RhoG= −0.22, p-value=0.047) and a negative genetic correlation between Tb.BMD and VFrx (RhoG = 0.29, p-value = 0.015).

Table 2
Heritability of vertebral fracture and vertebral volumetric bone mineral density and cross-sectional area


We found in this sample of European ancestry that genetic factors contribute significantly to VFrx, lumbar vertebral vBMD and CSA. Risk factors for VFrx, including integral vBMD and trabecular vBMD of the L3 vertebra, and L3 CSA, demonstrated moderate to high heritability (27% to 51%), while VFrx, defined as at least 25% height reduction in any vertebral level, was highly heritable (>43%). Spinal BMD was genetically correlated with VFrx suggesting that there are common genetic variants affecting VFrx and its risk factors. Although we observed a higher heritability when examining only vertebral fractures of grades two and higher, the number of such vertebral fractures was relatively small (n=148), therefore the results should be confirmed by family studies with larger numbers of VFrx.

Studying the heritability of VFrx, spinal vBMD and CSA is important since no one has yet investigated the genetic variants associated with VFrx. These results provide a strong justification to pursue the genetic contributors to these morbid and mortal fractures. This study highlights three major findings using spine images from the large family-based cohort, the Framingham Study. First, VFrx are heritable with estimates of h2 between 43–53% depending on the adjustment model. This heritability estimate is comparable to other types of the osteoporotic fracture (13,14), though lower than that for hip fractures occurring earlier in life (65%) (15) which is similar to the VFrx heritability estimate (around 68%) with adjustment for vBMD or CSA. Second, QCT-derived measurements of volumetric bone density and geometry are moderately-to-highly heritable (Tb.BMD ~41%, In.BMD ~43%, and CSA ~37%). Similarly, DXA-derived BMD and bone geometry are under strong genetic control with 50–70% heritability. Importantly, adjusting CSA for age, sex, and body size noticeably increased h2 of CSA, implying that genetic factors may operate independently of effects on body size. We also observed that the h2 of CSA remained the same (37.5%) if models included weight adjustment (the result is not shown).

Third, there were moderate genetic correlations between VFrx and multivariable-adjusted L3 In.BMD (RhoG = −0.22, p=0.047) and Tb.BMD (RhoG= −0.29, p=0.015). Although these genetic correlations are higher than both phenotypic and environmental correlations (RhoE), they are modest in magnitude. Genetic contributions to a risk factor, vBMD and CSA, may differ from the genetic contribution to the ultimate disease phenotype, VFrx in this study. This is especially true for a “proxy phenotype,” such as DXA-derived aBMD, and a complex event, such as VFrx. Similarly, Deng and colleagues (16) estimated that less than 1% of additive genetic variance is shared between aBMD and fractures at the hip, which strengthens the argument that there are other (non-skeletal) factors that contribute to the risk of fractures. Our estimates based on single vertebral level L3 are higher than those previously reported (with genetic correlation 0.26 for In.BMD and 0.07 for CSA), but the common genetic contribution to VFrx shared with the vBMD or CSA is still modest. Neither vBMD nor CSA alone can individually serve as a perfect surrogate of the skeleton’s ability to withstand the forces that produce fractures. Similarly, they are not perfect proxies for genetic studies of VFrx either.

There are several caveats worth noting. We conducted this analysis using up to 4,099 available samples of European ancestry in the Framingham Study. Results may not be generalizable to other ethnic groups. Moreover, the prevalence of osteoporotic fractures varies across different ethnicities. Research has shown the highest fracture rate in women of European ancestry while the mortality rate due to fracture is higher in African Americans (17). Studying prevalent VFrx is another limitation and future studies should examine incident VFrx. The bivariate analyses were conducted using SOLAR, which treats the phenotype as quantitative in the correlation. This may lead to inflated type I error for our VFrx measurement which was dichotomous.

In conclusion, our study suggests that prevalent vertebral fracture, vertebral vBMD and CSA in Caucasian adults are heritable and some of the heritability is shared between the traits, underscoring the importance of further work to identify the specific variants underlying genetic susceptibility to vertebral fracture and its risk factors.


This work was supported by NIH R01AR053986, R01AR/AG041398, T32 AG023480, K01 AR053118, and by the National Heart, Lung, and Blood Institute (NHLBI), Framingham Heart Study (NIH/NHLBI Contract N01-HC-25195).



All authors state that they have no conflicts of interest.

Authors’ roles: Study design and data collection: DK, MLB, LAC, EJS, DPK. Data analysis: CTL, YZ. Data interpretation: CTL, DK, MLB, LAC, DPK. Drafting manuscript: CTL, DK, DPK. Revising manuscript content: CTL, DK, YZ, YHH, HKG, KEB, TFL, EJS, SD, MLB, LAC, DPK. Approving final version of manuscript: CTL, DK, YZ, YHH, HKG, KEB, TFL, EJS, SD, MLB, LAC, DPK. CTL, DK and DPK take responsibility for the integrity of the data analysis.


1. Riggs BL, Melton LJ. The Worldwide Problem of Osteoporosis - Insights Afforded by Epidemiology. Bone. 1995;17(5):S505–S511. [PubMed]
2. Klotzbuecher CM, Ross PD, Landsman PB, Abbott TA, Berger M. Patients with prior fractures have an increased risk of future fractures: A summary of the literature and statistical synthesis. Journal of Bone and Mineral Research. 2000;15 (4):721–739. [PubMed]
3. Ensrud KE, Thompson DE, Cauley JA, Nevitt MC, Kado DM, Hochberg MC, Santora AC, Black DM, Gr FITR. Prevalent vertebral deformities predict mortality and hospitalization in older women with low bone mass. Journal of the American Geriatrics Society. 2000;48(3):241–249. [PubMed]
4. Ebbesen EN, Thomsen JS, Beck-Nielsen H, Nepper-Rasmussen HJ, Mosekilde L. Lumbar vertebral body compressive strength evaluated by dual-energy X-ray absorptiometry, quantitative computed tomography, and ashing. Bone. 1999;25 (6):713–24. [PubMed]
5. Duan YB, Duboeuf F, Munoz F, Delmas PD, Seeman E. The fracture risk index and bone mineral density as predictors of vertebral structural failure. Osteoporosis International. 2006;17(1):54–60. [PubMed]
6. Goodpaster BH, Kelley DE, Thaete FL, He J, Ross R. Skeletal muscle attenuation determined by computed tomography is associated with skeletal muscle lipid content. Journal of Applied Physiology. 2000;89(1):104–110. [PubMed]
7. Goodpaster BH, Thaete FL, Kelley DE. Composition of skeletal muscle evaluated with computed tomography. In Vivo Body Composition Studies. 2000;904:18–24. [PubMed]
8. Samelson EJ, Hannan MT, Zhang YQ, Genant HK, Felson DT, Kiel DP. Incidence and risk factors for vertebral fracture in women and men: 25-year follow-up results from the population-based Framingham Study. Journal of Bone and Mineral Research. 2006;21(8):1207–1214. [PubMed]
9. Samelson EJ, Christiansen BA, Demissie S, Broe KE, Zhou Y, Meng CA, Yu W, Cheng X, O’Donnell CJ, Hoffmann U, Genant HK, Kiel DP, Bouxsein ML. Reliability of vertebral fracture assessment using multidetector CT lateral scout views: the Framingham Osteoporosis Study. Osteoporosis International. 2011;22(4):1123–1131. [PMC free article] [PubMed]
10. Genant HK, Wu CY, Vankuijk C, Nevitt MC. Vertebral Fracture Assessment Using a Semiquantitative Technique. Journal of Bone and Mineral Research. 1993;8(9):1137–1148. [PubMed]
11. Lang T, LeBlanc A, Evans H, Lu Y, Genant H, Yu A. Cortical and trabecular bone mineral loss from the spine and hip in long-duration spaceflight. J Bone Miner Res. 2004;19(6):1006–12. [PubMed]
12. Almasy L, Blangero J. Multipoint quantitative-trait linkage analysis in general pedigrees. American Journal of Human Genetics. 1998;62(5):1198–1211. [PubMed]
13. MacGregor AJ, Snieder H, Spector TD. Genetic factors and osteoporotic fractures in elderly people - Twin data support genetic contribution to risk of fracture. British Medical Journal. 2000;320(7250):1669–1670. [PMC free article] [PubMed]
14. Recker RR, Deng HW, Chen WM, Recker S, Stegman MR, Li JL, Davies KM, Zhou Y, Deng HY, Heaney R. Genetic determination of Colles’ fracture and differential bone mass in women with and without Colles’ fracture. Journal of Bone and Mineral Research. 2000;15(7):1243–1252. [PubMed]
15. Michaelsson K, Melhus H, Ferm H, Ahlbom A, Pedersen NL. Genetic liability to fractures in the elderly. Arch Intern Med. 2005;165(16):1825–30. [PubMed]
16. Deng HW, Mahaney MC, Willams JT, Li J, Conway T, Davies KM, Li JL, Deng HY, Recker RR. Relevance of the genes for bone mass variation to susceptibility to osteoporotic fractures and its implications to gene search for complex human diseases. Genetic Epidemiology. 2002;22(1):12–25. [PubMed]
17. Cauley JA. Defining Ethnic and Racial Differences in Osteoporosis and Fragility Fractures. Clinical Orthopaedics and Related Research. 2011;469(7):1891–1899. [PMC free article] [PubMed]