Determinants of trabecular bone score and prevalent vertebral fractures in women with fragility fractures: a cross-sectional sub-study of NoFRACT.

Determinants of trabecular bone score (TBS) and vertebral fractures assessed semiquantitatively (SQ1-SQ3) were studied in 496 women with fragility fractures. TBS was associated with age, parental hip fracture, alcohol intake and BMD, not SQ1-SQ3 fractures. SQ1-SQ3 fractures were associated with age, prior fractures, and lumbar spine BMD, but not TBS.
INTRODUCTION: Trabecular bone score (TBS) and vertebral fractures assessed by semiquantitative method (SQ1-SQ3) seem to reflect different aspects of bone strength. We therefore sought to explore the determinants of and the associations between TBS and SQ1-SQ3 fractures.
METHODS: This cross-sectional sub-study of the Norwegian Capture the Fracture Initiative included 496 women aged ≥ 50 years with fragility fractures. All responded to a questionnaire about risk factors for fracture, had bone mineral density (BMD) of femoral neck and/or lumbar spine assessed, TBS calculated, and 423 had SQ1-SQ3 fracture assessed.
RESULTS: Mean (SD) age was 65.6 years (8.6), mean TBS 1.27 (0.10), and 33.3% exhibited SQ1-SQ3 fractures. In multiple variable analysis, higher age (βper SD = - 0.26, 95% CI: - 0.36,- 0.15), parental hip fracture (β = - 0.29, 95% CI: - 0.54,- 0.05), and daily alcohol intake (β = - 0.43, 95% CI - 0.79, - 0.08) were associated with lower TBS. Higher BMD of femoral neck (βper SD = 0.34, 95% CI 0.25-0.43) and lumbar spine (βper SD = 0.40, 95% CI 0.31-0.48) were associated with higher TBS. In multivariable logistic regression analyses, age (ORper SD = 1.94, 95% CI 1.51-2.46) and prior fragility fractures (OR = 1.71, 95% CI 1.09-2.71) were positively associated with SQ1-SQ3 fractures, while lumbar spine BMD (ORper SD = 0.75 95% CI 0.60-0.95) was negatively associated with SQ1-SQ3 fractures. No association between TBS and SQ1-SQ3 fractures was found.
CONCLUSION: Since TBS and SQ1-SQ3 fractures were not associated, they may act as independent risk factors, justifying the use of both in post-fracture risk assessment.

Introduction

In recent years, trabecular bone score (TBS) and vertebral fracture assessment (VFA) have become established methods in clinical use, providing supplemental information on bone strength and future fracture risk [1-3]. Identification of the determinants of these indices is important for further understanding the pathophysiology of fracture risk and for identification of modifiable factors to prevent future fractures.

TBS is a grey-level textural index of trabecular bone microarchitecture obtained from dual-energy absorptiometry (DXA) images of the lumbar spine [4]. It is, however, still subject to discussion, which bone properties TBS actually reflects [5]. TBS predicts fragility fractures independently of bone mineral density (BMD) [2, 6–8]. The use of TBS for adjustment of 10-year probability of a hip or other major osteoporotic fracture (MOF) calculated by the Fracture Risk Assessment Tool (FRAX) is endorsed [9]. The retrospective Manitoba study is the largest study on determinants of TBS (n = 29,407 women)[10]. They found that TBS was negatively associated with glucocorticoid use, prior major fracture, rheumatoid arthritis, chronic obstructive pulmonary disease, high alcohol intake and higher body mass index (BMI) and positively associated with femoral neck and lumbar spine BMD, and 7–11% of the variation in TBS was explained by BMD [10]. Genetic factors explained approximately 45% of the variance in TBS in healthy Vietnamese subjects, whereas sex, age, and height accounted for about 28% of the total variance in TBS [11].

The presence of a vertebral fracture doubles the risk of subsequent fracture at any given BMD [12]. VFA of lateral spine x-rays is therefore recommended as a part of the post-fracture risk assessment [13]. With higher resolution and image quality of modern DXA equipment, VFA of spine images obtained by DXA has become easily available and provides reliable information on vertebral fracture status in the thoracolumbar spine (T4–L4) at a much lower radiation dose compared to conventional lateral spine x-rays. Determinants as age, BMD, height loss, and prior non-vertebral fractures are associated with prevalent vertebral fractures in a large population-based cohort study [14]. In the Rotterdam Study, incident vertebral fractures were reported to be strongly associated with prevalent vertebral fractures, early menopause, current smoking, use of walking aids, and low BMD [15]. In addition, the population-based European Prospective Osteoporosis Study (EPOS) showed that low BMD, low BMI, and late menarche were associated with higher incidence of vertebral fractures while use of hormonal replacement therapy (HRT) was protective [16]. No associations between lifestyle factors such as smoking, alcohol intake, physical activity, or milk consumption and incident vertebral fracture were found in that study. A heritability of vertebral fractures (vertebral height reduction > 25%) of more than 43% has been demonstrated [17].

Studies on the association between TBS and vertebral fractures diverge. Low TBS has been shown to be associated with a higher prevalence of vertebral fractures on VFA in women [18, 19], but there are also studies contradicting this notion [20]. In the prospective Manitoba study, TBS was associated with incident vertebral fractures in women but not in men [21]. Moreover, TBS was not associated with incident vertebral fractures in the men in the MrOs study [22]. In a previous publication of results from this sub-study of Norwegian Capture the Fracture Initiative (NoFRACT), we reported that more than half of the patients who had sustained a fragility fracture exhibited low TBS (≤ 1.23), prevalent fracture on VFA or both [23]. Only 14% of these patients exhibited osteoporosis at the femoral neck, which increased to 28% when using the site with lowest BMD T-score. In the patients with osteopenia at the femoral neck, 34% exhibited low TBS, 33% vertebral fracture on VFA, but only 13% exhibited both low TBS and prevalent vertebral fracture. This is suggesting that TBS and VFA captured different aspects of bone strength. To better understand the potential pathophysiology underlying this diversity, identification of the determinants of both traits in the same cohort is of interest. Keeping in mind that heritability is a major determinant, investigation of clinical phenotypes and modifiable environmental risk factors may enable identification and treatment of individuals at risk. To our knowledge, this has not been studied in individuals with fractures. The objectives of this study were therefore (i) to explore the determinants of TBS and prevalent vertebral fractures on VFA, (ii) to explore whether prevalent vertebral fractures are determinants of TBS, and (iii) to explore whether TBS is a determinant of vertebral fractures in a cohort of women with fragility fractures.

Material and methods

Study subjects

NoFRACT is a multicenter study at seven hospitals in Norway with 34,976 persons with fragility fractures enrolled by Jan 2019 [24]. The main aim is to investigate the effect of introducing a standardized intervention program consisting of a Fracture Liaison Service model of care for identification, assessment, and treatment of osteoporosis in patients with fragility fractures on the rate of subsequent fractures. Both women and men aged 50 years and older, who have recently sustained a clinical fragility fracture, are eligible to the intervention, with the exception of patients with fracture in fingers, toes, skull, and face.

This cross-sectional sub-study (NCT02608801) of NoFRACT (NCT02536898) included patients at the University Hospital of North Norway, Tromsø, from 1 Oct 2015 to 31 Dec 2017 and at Drammen Hospital in the south-eastern part of Norway from 1 Jan 2016 to 31 Dec 2017 [23, 25]. Of all patients of 50 years of age and above, attending to these hospitals with a fragility fracture, more than 90% (n = 2682) were identified and offered fracture risk assessment. For elderly in-patients with fractures of hip, vertebrae, a total of two or more fragility fractures, or 10-year probability of MOF ≥ 20% calculated using FRAX, the treatment decision was made without the need of DXA scan (n = 1235). The participants were recruited among those who were referred to DXA (n = 1447), of whom 839 consented to participate in the study, 675 women and 164 men [23]. We included a total of 496 women who all had responded to a questionnaire on risk factors for fracture, had valid measurement of TBS, and BMD of the femoral neck and lumbar spine, and 423 of them had VFA performed. No patients were excluded due to conditions known to affect bone metabolism, such as chronic kidney disease, use of anti-osteoporosis drugs, HRT, or premenopausal status. The study was approved by The Regional Committee for Medical and Health Research Ethics (REK 2014/2260).

Variables

The index fractures leading to inclusion were fragility fractures of the forearm (n = 196), ankle (n = 90), proximal humerus (n = 68), hip (n = 36), vertebrae (thoracic and lumbar) (n = 23), and other sites (n = 83). No women had rib or cervical fractures as index fracture. The vertebral fractures that led to inclusion were diagnosed on x-ray, CT, or MRI, not by VFA. Information on number of fractures after the age of 50, number of falls during the last 12 months before inclusion, parental history of hip fractures (yes vs. no), comorbidity, medication, postmenopausal status (yes vs. no), age at menopause, number of children, total months of breastfeeding, smoking (yes vs. no), daily alcohol intake (yes vs. no), physical activity (h/week), and dairy products (units/day) were collected through a questionnaire.

Height and weight was measured without shoes and heavy clothing, and BMI was calculated as weight per square meter height. BMD was measured at the lumbar spine (L1–L4), femoral neck, and total hip at both sides and using DXA Prodigy Pro in Tromsø and iDXA Pro in Drammen (both GE Lunar, Madison, WI, USA). Phantom Quality Assurance (QA) of the DXA equipment was performed daily. All fractured lumbar vertebrae were excluded, and BMD T-scores at femoral neck and total hip were calculated using the Third National Health and Nutrition Examination Survey with reference data of female Caucasians aged 20–29 years [26]. iNsight software (MediMaps, Geneva, Switzerland) version 3.0.1 was used to calculate TBS values from the DXA scans used for lumbar spine BMD (L1–L4), and fractured vertebrae were excluded. European (Medimaps) reference population was used.

Lateral images of the thoracolumbar spine (T4–L4) were obtained from DXA scans with the patient in a lateral decubitus position with lumbar support and hips flexed 90 degrees. An experienced physician (TTB) performed VFA and grading of the fracture severity using Genant’s method for the semiquantitative (SQ) vertebral fracture scoring [27]. This combines the visual identification of fracture of vertebral bodies (height loss of the anterior, middle, posterior, or the whole vertebra) and grading of the fracture by percentage of height loss. A SQ score of 0 (SQ0) (< 20% height loss) is a non-fractured vertebra, SQ1 (20–25% height loss) a mild fracture, SQ2 (25–40% height loss) a moderate fracture, and SQ3 (≥ 40% height loss) a severe fracture. We classified the presence of one or more SQ1, SQ2, or SQ3 fractures as SQ1–SQ3 fractures, one or more SQ2 or SQ3 fractures as SQ2–SQ3 fractures, and presence of at least one SQ3 fracture as ≥ SQ3. The inter-observer agreement of SQ1–SQ3 fractures tested against another experienced clinician showed a κ of 0.84 (95% CI 0.70, 0.98) [23].

Statistical analyses

Characteristics of the cohort were calculated as mean ± standard deviation (SD) for the continuous variables and number with percentages (%) for categorical variables. Continuous variables were checked for normality using quantile-quantile (QQ) plot. All these variables were normally distributed except breastfeeding which was log-transformed in further analyses. Scatterplots were performed between continuous variables and visually checked for linearity. Univariable linear regression analyses were performed to investigate associations between the outcome variable TBS and clinical relevant determinants (age, BMI, history of prior fractures after the age of 50, falls during the last 12 months before inclusion, parental history of hip fractures, comorbidities, use of medications, number of children, breastfeeding, currently smoking, daily alcohol consumption, physical activity, consumption of dairy products, prevalence of SQ1–SQ3 fractures, femoral neck and lumbar spine BMD). Only determinants significant at p-level < 0.10 were retained and included in multivariable models. Multiple linear regression analyses were performed, and because of potential multi-collinearity between TBS, femoral neck, and lumbar spine BMD, we tested different models for each of the traits to explore the attributional variance in the outcome variable by the introduced determinants. Non-significant determinants were removed one by one until the exposure variables with statistical significant association remained (p < 0.05). Determinants that had been removed were reintroduced one by one to re-check for significance. Results are presented as β coefficients with 95% confidence intervals (CIs), p values, and explained variance (R2). The same procedure was performed with femoral neck and lumbar spine BMD as outcome variables. Univariable logistic regression analyses were performed to explore associations between SQ1–SQ3 fractures (yes vs. no) as the outcome variable and the determinants used in the linear regression models (listed above). Variables with significant association at p level < 0.10 were included in the multivariable logistic regression analyses. We tested different models including TBS, femoral neck, and lumbar spine BMD as determinants to evaluate the association with SQ1–SQ3 fractures. Determinants without significant association with the outcome SQ1–SQ3 fractures were removed one by one until the final model contained only the exposure variables with significant association (p < 0.05). Then, the removed variables were reintroduced one by one to re-check for significance. Results are presented by odds ratio (OR) with 95% CI. Evaluation of the predictive accuracy of the models was assessed by calibration and discrimination. Calibration was evaluated by the Hosmer and Lemeshow goodness-of-fit test. A statistically non-significant Hosmer and Lemeshow result (p > 0.05) suggests that the fit of the model is acceptable. Discrimination of SQ1–SQ3 fractures was evaluated by analysis of the area under the receiver operating characteristic (ROC) curve. We defined acceptable discriminatory capability as an area under the ROC curve greater than 0.7. Standardized regression coefficients (βper SD) or odds ratio (ORper SD) with 95% CI was used to facilitate the comparison of the strength of the associations between the continuous exposure variables and the outcome variables. The number of determinants in the multiple linear regression models did not exceed 10% of the number of observations. The number of determinants in the multiple logistic regression models never exceeded 14, as there were 141 SQ1–SQ3 fractures (maximum 10 events per determinant). All analyses were performed using Stata v15 (Version 15, StataCorp LP, TX, USA).

Results

The mean age of the 496 women was 65.6 years (± 8.6) and 193 (42.1%) women had sustained one or more fractures prior to the index fracture after the age of 50 (Table 1). One in three exhibited SQ1–SQ3 fractures on VFA. Mean TBS was 1.27 and BMD at femoral neck, total hip, and lumbar spine was 0.795 g/cm2, 0.835 g/cm2, and 1.027 g/cm2, respectively. The mean number of falls during the last 12 months prior to inclusion was 1.3, about one in four women reported a history of parental hip fracture, and one in ten reported daily alcohol intake. The women had on average given birth to 2.2 children (including the 7.8% of the women who were nulliparous). The median duration of total time of breastfeeding was 10 months (range 0–96 months).

Table 1

Characteristics of 496 women 50 years and older with a recent fragility fracture

	Number	Total
Age, years	496	65.6 ± 8.6
Body mass index, kg/m2	496	26.1 ± 4.0
History of prior fracture after age of 50	459	193 (42.1)
Number of falls during the last 12 months	489	1.3 ± 0.7
Parental history of hip fracture	410	100 (24.4)
Ulcus/gastritis	476	56 (11.8)
Asthma/COLD	489	53 (10.8)
Type 1 and type 2 diabetes	495	23 (4.6)
Rheumatoid arthritis	488	22 (4.5)
Myocardial infaction	496	13 (2.6)
Cancer	495	71 (14.3)
Hypothyreosis	490	54 (11.2)
Hyperthyreosis	490	13 (2.7)
Use of antacids	483	76 (15.7)
Use of statins	488	119 (24.4)
Use of oral glucocorticoids	486	28 (5.8)
Menopausal estrogen supplementation	436	29 (6.7)
Use of AOD at baseline	446	41 (9.2)
Postmenopausal status	447	429 (96.0)
Age at menopause, years	399	48.7 ± 4.9
Number of children	432	2.2 ± 1.2
Nullipara	432	34 (7.8)
Breastfeeding, monthsa	325	10 (0 – 96)
No breastfeeding	325	45 (13.9)
Currently smoking	479	71 (14.8)
Daily alcohol intake	487	48 (9.8)
Physical activity, h/week	454	2.8 ± 1.8
Intake of dairy products, units/day	480	2.3 ± 1.1
Trabecular bone score lumbar spine	496	1.27 ± 0.10
SQ1-SQ3 fracture	423	141 (33.3)
Femoral neck BMD, g/cm2	482	0.795 ± 0.105
Total hip BMD, g/cm2	482	0.835 ± 0.118
Lumbar spine BMD, g/cm2	496	1.027 ± 0.166
Lowest BMD any site, T-score	496	− 2.1 ± 0.9

COLD chronic obstructive lung disease, AOD anti-osteoporosis drugs, SQ vertebral fractures assessed by semiquantitative method, BMD bone mineral density

Values are mean ± SD or n (%), except amedian (range)

Determinants of TBS

Higher age (βper SD = − 0.26, 95% CI − 0.36, − 0.15), a history of parental hip fracture (β = − 0.29, 95% CI − 0.54, − 0.05), and daily alcohol intake (yes vs. no) (β = − 0.43, 95% CI − 0.79, − 0.08) were associated with a lower TBS (Table 2; Fig. 1). Higher age at menopause (βper SD = 0.11, 95% CI 0.01, 0.21), higher BMD of the femoral neck (βper SD = 0.34, 95% CI 0.25, 0.43), and lumbar spine (βper SD = 0.40, 95% CI 0.31, 0.48) were associated with higher TBS. In models additionally including SQ1–SQ3 fractures, femoral neck or lumbar spine BMD, age at menopause were no longer associated with TBS. SQ1–SQ3 fractures were not associated with TBS in models including age and daily alcohol consumption. Replacing SQ1–SQ3 fractures with SQ2–SQ3 fractures or SQ3 fractures in multivariable analyses did not change the results. The model including the significant determinant age, parental hip fracture, daily alcohol consumption, number of children, and lumbar spine BMD explained 28% of the variance in TBS. Femoral neck BMD explained 8% and lumbar spine BMD 18% of the variance in TBS.

Table 2

Association between risk factors for fracture and 1 standard deviation (SD) difference in trabecular bone score (TBS)

	Univariable analyses		Multivariable modela		Modela+ SQ1-SQ3		Modela+ femoral neck BMD		Modela+ lumbar spine BMD
	β (95% CI)	p	β (95% CI)	p	β (95% CI)	p	β (95% CI)		β (95% CI)	p
Age per SD	− 0.29 (− 0.37, − 0.20)	< 0.001	− 0.26 (− 0.36, − 0.15)	< 0.001	− 0.25 (− 0.35,− 0.15)	< 0.001	− 0.16 (− 0.24, − 0.07)	0.001	− 0.24 (− 0.32, − 0.15)	< 0.001
BMI per SD	0.07 (0.02, 0,16)	0.116
History of prior fracture	− 0.19 (0.37, 0.01)	0.038
Parental hip fracture	− 0.25 (− 0.47, − 0.03)	0.028	− 0.29 (− 0.54, − 0.05)	0.017					− 0.28 (− 0.48, − 0.07)	0.008
Using AOD at baseline	− 0.32 (− 0.64, − 0.00)	0.047
Age at menopause per SD	0.08 (− 0.02, 0.18)	0.099	0.11 (0.01, 0.21)	0.036
Number of children	0.08 (− 0 00, 0.16)	0.058							0.09 (0.01, 0.17)	0.036
Daily alcohol consumption	− 0.33 (− 0.63, − 0.03)	0.030	− 0.43 (− 0.79, − 0.08)	0.017	− 0.47 (− 0.76, − 0.17)	0.002	− 0.32 (− 0.60, − 0.04)	0.027	− 0.41 (− 0.70, − 0.12)	0.006
SQ1− SQ3 fractures	− 0.31 (− 0.50, − 0.11)	0.003	–	–	− 0.14 (− 0.34, 0.07)	0.184	–		–	–
Femoral neck BMD per SD	0.41 (0.33, 0.50)	< 0.001	–	–	–	–	0.34 (0.25, 0.43)	< 0.001	–	–
Lumbar spine BMD per SD	0.44 (0.36, 0.52)	< 0.001	–	–	–	–	–	–	0.40 (0.31, 0.48)	< 0.001
R2	–		0.11		0.09		0.19		0.29

SQ vertebral fractures assessed by semiquantitative method, BMD bone mineral density, β beta coefficient in linear regression analysis, CI confidence interval

aIncluding age, BMI (body mass index), prior fracture, parental hip fracture, using AOD at baseline, age at menopause, number of children and daily alcohol intake. Only significant results are shown

Fig. 1

Associations between femoral neck and lumbar spine bone mineral density (BMD), trabecular bone score (TBS) and vertebral fractures on VFA (SQ1-SQ3) with attributed variance of their determinants. BMI body mass index

Determinants of SQ1–SQ3 fractures

Age (ORper SD = 1.94, 95% CI 1.51–2.46) and a history of prior fractures (OR = 1.71, 95% CI 1.09–2.71) were positively associated with SQ1–SQ3 fractures (Table 3). Lumbar spine BMD (ORper SD = 0.75 95% CI 0.60–0.95) was negatively associated with SQ1–SQ3 fractures in models including age and prior fracture. SQ1–SQ3 fractures were neither associated with TBS nor femoral neck BMD in models including age and prior fractures. In analyses with SQ2–SQ3 fractures and SQ3 fractures as outcome variables, the results were similar.

Table 3

Odds ratio (OR) of presence of SQ1-SQ3 fractures (yes vs. no) with possible risk factors for fracture

	Univariable analyses		Multivariable modela		Modela+ TBS		Modela+ femoral neck BMD		Modela+ lumbar spine BMD
	OR (95% CI)	p	OR (95% CI)	p	OR (95% CI)	p	OR (95% CI)	p	OR (95% CI)	p
Age per SD	2.02 (1.60,2.56)	< 0.001	1.94 (1.51, 2.46)	< 0.001	1.86 (1.44, 2.40)	< 0.001	1.87 (1.43, 2.45)	< 0.001	1.93 (1.49, 2.48)	< 0.001
BMI per SD	1.08 (0.89, 1.32)	0.437
History of prior fracture	2.37 (1.54, 3.64)	< 0.001	1.71 (1.09, 2.71)	0.021	1.71(1.08, 2.70)	0.022	1.66 (1.04, 2.65)	0.033	1.76 (1.11, 2.79)	0.017
Cerebrovascular event	4.20 (1.24, 14.19)	0.021
TBS per SD	0.73(0.59, 0.90)	0.003	–	–	0.81 (0.64, 1.03)	0.084	–		–	–
Femoral neck BMD per SD	0.70 (0.56, 0.87)	0.002	–	–	–	–	0.91 (0.71, 1.16)	0.452	–	–
Lumbar spine BMD per SD	0.70 (0.56, 0.87)	0.002	–	–	–	–	–		0.75 (0.60, 0.95)	0.017

TBS trabecular bone score, BMD bone mineral density, CI confidence interval

aIncluding age, BMI (body mass index), prior fracture, and cerebrovascular accident. Only significant or important results are shown

Determinants of femoral neck BMD

Higher age (βper SD = − 0.35, 95% CI: − 0.43, − 0.27) was associated with lower femoral neck BMD, while higher BMI (βper SD = 0.28, 95% CI: 0.20, 0.36) and higher TBS (βper SD = 0.32, 95% CI: 0.25, 0.40) were associated with higher femoral neck BMD (Table 4). Age, BMI, myocardial infarction, and TBS explained 32% of the variance in femoral neck BMD. TBS explained 10% of the variance in femoral neck BMD.

Table 4

Association between risk factors for fracture and 1 standard deviation (SD) difference in femoral neck bone mineral density (BMD)

	Univariable analyses		Multivariable modela		Modela+ TBS		Modela+ SQ1-SQ3		Modela+ SQ1-SQ3 + TBS
	β (95% CI)	p	β (95% CI)	p	β (95% CI)	p	β (95% CI)	p	β (95% CI)	p
Age per SD	− 0.38 (− 0.46, − 0.30)	< 0.001	− 0.35(− 0.43,− 0.27)	< 0.001	− 0.26 (− 0.33, − 0.17)	< 0.001	− 0.30 (− 0.40, − 0.21)	< 0.001	− 0.26 (− 0.33, − 0.18)	< 0.001
BMI per SD	0.31 (0.23, 0.40)	< 0.001	0.28 (0.20, 0.36)	< 0.001	0.27 (0.20, 0.35)	< 0.001	0.28 (0.20, 0.37)	< 0.001	0.27 (0.20, 0.35)	< 0.001
Prior fracture	− 0.30 (− 0.49, − 0.12)	0.002
Asthma/COLD	− 0.34 (− 0.62, − 0.05)	0.021
Myocardial infarction	− 0.62 (− 1.19, − 0.05)	0.035			− 0.50 (− 0.98, − 0.03)	0.039	− 0.60(− 1.14, − 0.06)	0.029	− 0.50 (− 0.98, − 0.03)	0.039
Using AOD at baseline	− 0.51 (− 0.83, − 0.19)	0.002					− 0.34 (− 0.66, − 0.02)	0.040
TBS per SD	0.41 (0.33, 0.49)	< 0.001	–	–	0.32 (0.25, 0.40)	< 0.001	-	–	0.29 (0.20, 0.38)	< 0.001
SQ1-SQ3 fractures	− 0.33 (− 0.53, − 0.13)	0.001	–	–	–	–	− 0.20 (− 0.39, 0.00)	0.051	− 0.12 (− 0.30, 0.05)	0.173
Lumbar spine BMD per SD	0.57 (0.49, 0.64)	< 0.001	–	–	–	–	–	–	–	–
R2	–		0.22		0.32		0.25		0.31

TBS trabecular bone score, SQ vertebral fractures assessed by semiquantitative method, β beta coefficients in linear regression analysis, CI confidence interval, COLD chronic obstructive lung disease, R explained variance

aIncluding age, BMI (body mass index), prior fracture, asthma/COLD, myocardial infarction and previous AOD (anti-osteoporosis drugs). Only significant results are shown

Determinants of lumbar spine BMD

SQ1–SQ3 fractures (βper SD = − 0.25, 95% CI − 0.43, − 0.07) were associated with a lower lumbar spine BMD, while higher age (βper SD = 0.20, 95% CI 0.12, 0.29), higher BMI (βper SD = 0.20, 95% CI 0.12, 0.29), diabetes (β = 0.61, 95% CI: 0.21, 1.00) and high TBS (βper SD = 0.43, 95% CI 0.34, 0.51) were associated with higher lumbar spine BMD (Table 5). TBS and SQ1–SQ3 fractures explained 17% and 1% of the variance in lumbar spine BMD, respectively.

Table 5

Association between risk factors for fracture and 1 standard deviation (SD) difference in lumbar spine bone mineral density (BMD)

	Univariable analyses		Multivariablemodela		Modela+ TBS		Modela+ SQ1-SQ3		Modela+ SQ1-SQ3 + TBS
	β (95% CI)	p	β (95% CI)	p	β (95% CI)	p	β (95% CI)	p	β (95% CI)	p
Age per SD	− 0.14 (− 0.22, − 0.05)	0.003	− 0.12 (− 0.21, − 0.03)	0.006					0.09 (0.00, 0.18)	0.045
BMI per SD	0.24 (0.15, 0.32)	< 0.001	0.22 (0.13, 0.30)	< 0.001	0.19 (0.12, 0.27)	< 0.001	0.24 (0.15, 0.33)	< 0.001	0.20 (0.12, 0.29)	< 0.001
Fallsb per SD	− 0.10 (− 0.22, 0.02)	0.088	− 0.12 (− 0.23, − 0.00)	0.049
Diabetes	0.60 (0.19, 1.02)	0.005	0.51 (0.10, 0.92)	0.014	0.65 (0.28, 1.02)	0.001	0.44 (0.02, 0.87)	0.042	0.61 (0.21, 1.00)	0.003
Rheumatoid arthritis	0.40 (− 0.06, 0.80)	0.092
Smoking	0.23 (− 0.03, 0.48)	0.079
TBS per SD	0.44 (0.36, 0.52)	< 0.001	–		0.44 (0.36, 0.51)	< 0.001	–	–	0.43 (0.34, 0.51)	< 0.001
SQ1-SQ3 fractures	− 0.32 (− 0.52, − 0.12)	0.001	–		–	–	− 0.33 (− 0.52, − 0.14)	0.001	− 0.25 (− 0.43, − 0.07)	0.008
FN BMD per SD	0.54 (0.47, 0.61)	< 0.001	–		–	–	–	–	–	–
R2	–		0.09		0.25		0.10		0.26

TBS trabecular bone score, SQ vertebral fractures assessed by semiquantitative method, β beta coefficients in linear regression analysis, CI confidence interval, FN femoral neck

aIncluding age, BMI (body mass index), falls, diabetes, rheumatoid arthritis and smoking. Only significant results are shown

bNumber of falls the last 12 months prior to inclusion

Discussion

No association between TBS and SQ1–SQ3 fractures could be demonstrated in this cohort of women with prevalent fragility fractures. We found that higher age, a history of parental hip fracture, and daily alcohol consumption were associated with lower TBS, while higher femoral neck and lumbar spine BMD were associated with higher TBS. Age and prior fractures were positively associated with SQ1–SQ3 fractures, whereas lumbar spine BMDs were negatively associated with SQ1–SQ3 fractures.

The unexpected finding of no association between TBS and SQ1–SQ3 fractures has previously been shown in a study of elderly Swedish women with prevalent vertebral fractures [20] and in men in the Manitoba study and MrOs study with incident vertebral fractures [21, 22]. There are other studies showing an association between TBS and SQ1–SQ3 fractures in women [18, 19]. We do not know why these findings diverge between studies, but several factors may be of importance. Firstly, genetic factors explain a large part of the variation in TBS and vertebral fractures, and genetic factors may vary significantly between countries and the cohorts studied. Scandinavian women are a genetically homogenous group, and they have the highest prevalence of vertebral fractures in Europe [28]. Secondly, TBS is derived from the same images as lumbar spine BMD. Lumbar spine BMD accounted for a larger contribution to the variance of TBS (20%) than femoral neck BMD (10%) in our cohort. SQ1–SQ3 fractures were negatively associated with lumbar spine BMD, and only about 1% of the variance in lumbar spine BMD was explained by SQ1–SQ3 fractures, suggesting a weak association. No association between SQ1–SQ3 fractures and femoral neck BMD was found. We might infer that in our cohort of women with fragility fractures, SQ1–SQ3 fractures reflect reduced bone strength that is not captured by BMD nor TBS.

The determinants of TBS identified in our study differ somewhat from what was found in the population-based Canadian Manitoba study [10]. We found a negative association of TBS with age and daily alcohol intake, but no associations with BMI, glucocorticoid use, rheumatoid arthritis, or chronic obstructive pulmonary disease. This could be due to the smaller sample size in our study compared to the Manitoba study (496 vs. 29,407) or that the source of the data differed. We collected self-reported information on lifestyle, comorbidities, and medication, with limitations regarding self-reported data (i.e. possible information bias). In the Manitoba study, this information was retrieved from register-based data of physician billing claims, hospital discharge abstracts, and provincal retail pharmacy database, with uncertainities regarding the validity of registered information (i.e. whether the discharge record reflects the disease at hand). Both studies reported a positive association between femoral neck and lumbar spine BMD and TBS. The variance in TBS explained by femoral neck and lumbar spine BMD in the Manitoba women was 7% and 11%, respectively, compared to 20% and 28% in the NoFRACT women. This might be explained by the inclusion of healthier subjects in the Manitoba cohort where only 14% had a previous major fracture, compared to the NoFRACT cohort, where all the women had a clinical fragility fracture. In addition, patients with BMI > 37 kg/m2 were not excluded in the Manitoba study, which may have diluted the association between BMD and TBS.

We found higher OR for SQ1–SQ3 fractures in patients with a history of a prior fracture, increasing age and decreasing lumbar spine BMD, in line with other studies [29-31]. No association between SQ1–SQ3 fractures and body weight or BMI was found in adjusted models, which is in accordance with some studies [29-31], but contrary to others [32, 33]. Studies on incident vertebral fractures, however, show similar results. In the EPOS study, including 3402 women, late menarche (after 16 years of age) was associated with an increased risk of incident vertebral fracture, whereas HRT, increasing body weight, and BMI were protective [16]. Lifestyle factors such as smoking, alcohol intake, physical activity, or milk consumption showed no association with incident vertebral fractures, in line with our findings in patients with prevalent vertebral fractures. In the prospective population-based Rotterdam Study including 2467 women, age, low BMD, prevalent vertebral fractures, early menopause (< 45 years of age), currently smoking, and walking aid use were associated with incident vertebral fractures in women [15]. Prevalent and incident vertebral fractures are not comparable, although the risk factors should be similar.

To the best of our knowledge, this is the first study on determinants of TBS in a cohort of women with fragility fractures. However, there are some limitations. The study sample was perhaps too small to show associations with the main outcome variables (TBS, SQ1–SQ3 fractures) and certain risk factors, such as glucocorticoid use, rheumatoid arthritis, and chronic obstructive pulmonary disease. On the other hand, we showed that more of the variance in TBS could be explained by BMD than in the Manitoba study, maybe due to exclusion of obese patients or the selection of patients who all had suffered at least one fragility fracture. As a result of the cross-sectional design, it was not possible to determine the temporal nature of any observed associations.

In conclusion, no statistical significant association between TBS and SQ1–SQ3 fractures was found. Higher age, a history of parental hip fracture, and daily alcohol consumption were associated with lower TBS, while higher femoral neck and lumbar spine BMD were associated with higher TBS. Age and a history of prior fractures were positively associated with SQ1–SQ3 fractures, while lumbar spine BMD was negatively associated with SQ1–SQ3 fractures. Modifiable risk factors such as daily alcohol intake and low BMD are important targets in fracture prevention strategies. Since TBS and SQ1–SQ3 fractures are not associated, each of them may act as independent risk factors for fracture, justifying the use of both in post-fracture risk assessment.