Low Predictive Value of FRAX Adjusted by Trabecular Bone Score for Osteoporotic Fractures in Korean Women: A Community-Based Cohort Study.

BACKGROUND: The value of the Fracture Risk Assessment Tool (FRAX) and the trabecular bone score (TBS) for assessing osteoporotic fracture risk has not been fully elucidated in Koreans. We conducted this study to clarify the predictive value of FRAX adjusted by TBS for osteoporotic fractures in Korean women.
METHODS: After screening 7,192 eligible subjects from the Ansung cohort, 1,165 women aged 45 to 76 years with available bone mineral density (BMD) and TBS data were enrolled in this study. We assessed their clinical risk factors for osteoporotic fractures and evaluated the predictive value of FRAX with or without BMD and TBS.
RESULTS: During the mean follow-up period of 7.5 years, 99 (8.5%) women suffered major osteoporotic fractures (MOFs) and 28 (2.4%) experienced hip fractures. FRAX without BMD, BMD-adjusted FRAX, and TBS-adjusted FRAX were significantly associated with the risk of MOFs (hazard ratio [HR] per percent increase, 1.08; 95% confidence interval [CI], 1.03 to 1.14; HR, 1.09; 95% CI, 1.03 to 1.15; and HR, 1.07; 95% CI, 1.02 to 1.13, respectively). However, BMD-adjusted FRAX and TBS-adjusted FRAX did not predict MOFs better than FRAX without BMD based on the Harrell's C statistic. FRAX probabilities showed limited value for predicting hip fractures. The cut-off values of FRAX without BMD, FRAX with BMD, and FRAX with BMD adjusted by TBS for predicting MOFs were 7.2%, 5.0%, and 6.7%, respectively.
CONCLUSION: FRAX with BMD and TBS adjustment did not show better predictive value for osteoporotic fractures in this study than FRAX without adjustment. Moreover, the cut-off values of FRAX probabilities for treatment might be lower in Korean women than in other countries.

INTRODUCTION

Osteoporosis is defined as a fragile bone status characterized by low bone mass, bone microarchitecture deterioration, and reduced bone strength, which increase the risk of fragility fractures [1]. Previous studies reported that mortality significantly increased after osteoporotic vertebral and hip fractures [2-4]. In Korea, the risk of experiencing osteoporotic fractures during the remaining lifetime in persons older than 50 years was estimated to be 23.8% for men and 59.5% for women [5]. Thus, the precise prediction of osteoporotic fractures would be helpful to reduce their consequent burden.

Bone mineral density (BMD) measured using dual-energy X-ray absorptiometry (DXA) is known to be a major determinant of osteoporotic fracture risk [6]. However, BMD provides little information about bone quality, including the bone microarchitecture. Furthermore, severe multiple vertebral compression fractures, degenerative changes, or aortic calcification can cause a false increase in spinal BMD [7]. Moreover, previous studies reported that most fragility fractures occurred in people with normal or osteopenic BMD values [8,9]. Recently, the lumbar spine trabecular bone score (TBS) was introduced as a measurement of bone texture that is correlated with fracture risk independently of BMD [10-12]. Several studies have reported that a low TBS was correlated with an increased fracture risk [13-15].

Clinical risk factors such as parental history of hip fracture, exposure to systemic glucocorticoids, history of prior fragility fracture, current smoking, high alcohol intake, and the presence of rheumatoid arthritis play a pivotal role in the assessment of osteoporotic fracture risk [16-20]. Although several fracture risk prediction tools based on clinical risk factors have been created, the Fracture Risk Assessment Tool (FRAX), which was developed in 2008 based on large prospective cohort studies, is widely used worldwide [21]. However, the risk factors and incidence of osteoporotic fractures vary widely among populations [22] and the incidence of osteoporotic fractures in Korea is known to be lower than that in Western countries [23,24]. The calculated FRAX probability was also lower in Koreans than in Western populations [25,26]. Thus, population-specific data on the treatment threshold for FRAX are required to predict the risk of fracture. Moreover, the value of the recently introduced TBS-adjusted FRAX for predicting osteoporotic fractures has not been validated in Koreans.

We aimed to evaluate the predictive value of the TBS-adjusted FRAX for osteoporotic fractures, and to suggest cut-off values of FRAX probabilities for osteoporosis treatment in this community-based cohort study of Korean women.

METHODS

Study population

The Ansung cohort was established for the Korean Genome and Epidemiology Study in 2001 to investigate the prevalence and determinants of chronic diseases in Korea. This is an ongoing prospective study involving biennial examinations. In Ansung, sampling was based on a clustering sampling method using mailing, door-to-door visits, and telephone requests within five randomly selected local townships (termed myeon in Korea) in the district. After screening 7,192 eligible subjects, a total of 5,018 subjects (2,239 men and 2,279 women) aged between 40 and 69 years were recruited in 2001. Finally, a total of 1,165 women aged 45 to 76 years who voluntarily consented to undergo BMD and TBS measurements during 2007 to 2008 were enrolled because TBS has been available since 2007. The period encompassed by the study extended until 2015 to 2016, and the mean follow-up period was 7.5 years. Forty subjects died during this period. Clinical risk factors for osteoporotic fractures were investigated using self-administered questionnaires, which gathered information on age, menopausal status, history of prior osteoporotic fracture, parental history of hip fragility fracture, smoking status, alcohol consumption, and comorbidities. Smoking status was categorized as never, previous, or current. Alcohol intake history was assessed including the types and quantity of alcohol. All data were collected by trained staff using standardized questionnaires through face-to-face interviews. The questionnaires collected information on all types of arthritis without distinguishing rheumatoid arthritis from osteoarthritis.

The study was approved by the Institutional Review Board of Ajou University Hospital and written informed consent was obtained from all patients (AJIRB-CRO-07-012, AJIRB-MED-SUR-13-016, AJIRB-BMR-SUR-15-179, AJIRB-BMR-SUR-16-121).

Assessment of osteoporotic fractures

The subjects were contacted every 2 years and asked whether they had a new fracture through the self-questionnaire. We included only incident osteoporotic fractures of the spine, hip, humerus, and wrist as major osteoporotic fractures (MOFs). Osteoporotic fractures were defined as fragility fractures that occurred without major trauma or fractures due to a fall from a height lower than the subject’s height. We confirmed osteoporotic fractures only through the self-reported questionnaire because the time of occurrence of vertebral fractures observed on radiography without any medical history was unclear, and we could not assess osteoporotic fractures immediately during the follow-up interval.

Measurements of BMD and TBS

Lumbar spine BMD and femoral neck BMD were measured using DXA (Lunar Prodigy, General Electric Medical Systems, Chicago, IL, USA) and analyzed (Encore Software version 11.0, Minneapolis, MN, USA) according to the manufacturer’s instructions. When the lumbar spine BMD (L1–L4) was not appropriate because of a compression fracture or severe sclerosis, we used the mean BMD value of the remaining lumbar spine. All TBS values were retrospectively measured with iNsight software version 2.0.0.1 (MedImaps, Geneva, Switzerland), using spine DXA files from the database to ensure that the investigators were blinded to all clinical parameters and outcomes. The software used the raw DXA image of the anteroposterior spine with the same region of interest as that used for the BMD measurements. The instruments were calibrated using anthropomorphic phantoms.

FRAX and TBS-adjusted FRAX calculations

The FRAX without BMD (including age, sex, body weight, height, history of prior osteoporotic fracture, parental history of hip fragility fracture, current smoking, arthritis, and alcohol consumption >3 units/day), FRAX with BMD, and TBS-adjusted FRAX probabilities for MOFs and hip fractures were calculated using the tool for Koreans provided on the FRAX website (https://www.sheffield.ac.uk/FRAX/index.aspx).

Statistical analysis

Data were given as mean±standard deviation (SD) or number (%). We performed the chi-square test for categorical variables and the Student’s t test for continuous variables. Hazard ratios (HRs) with 95% confidence intervals (CIs) for the time to incident fracture were obtained from Cox proportional hazard models using Stata version 14.2 (StataCorp., College Station, TX, USA). The Harrell’s C statistic was used to compare the predictive capability of FRAX without BMD, FRAX with BMD, and TBS-adjusted FRAX probabilities. To evaluate the additional contribution of TBS to the prediction of osteoporotic fractures in osteopenic subjects, we conducted a subgroup analysis. The cut-off value of each FRAX probability was calculated using the ‘maximally selected rank statistics’ (maxstat) method through the ‘maxstat’ and ‘rpart’ packages of R version 4.3.2 (R Foundation for Statistical Computing, Vienna, Austria) [27]. The details of this method were described by Hothorn and Lausen [27,28]. Briefly, this method assesses all of the data points from the continuous FRAX probability data and establishes a cut-off point where the standardized statistics take their maximum significance for the separation of subjects at a high risk for osteoporotic fracture.

RESULTS

The baseline characteristics of the study subjects are shown in Table 1. The mean age±SD was 61.6±8.7 years. The study subjects comprised 1,072 (92.0%) postmenopausal women. In total, 123 women (10.6%) had a history of prior osteoporotic fracture. Their lumbar spine BMD, femoral neck BMD, and lumbar spine TBS were 0.957±0.179, 0.800±0.133, and 1.363± 0.097 g/cm2, respectively (Table 1). The BMD-adjusted FRAX probability for MOF and hip fracture was 5.4%±2.9% and 1.1%±1.6%, respectively. The TBS-adjusted FRAX probability was 5.4%±3.2% for MOF and 1.1%±1.6% for hip fracture.

Table 1

Baseline Characteristics of the Study Subjects

Characteristic	Women (n=1,165)
Age, yr	61.6±8.7
Height, cm	152.3±6.1
Weight, kg	57.8±8.7
BMI, kg/m2	24.9±3.4
Menopause	1,072 (92.0)
Previous history of osteoporotic fracture	123 (10.6)
Parental history of hip fracture	0 (0.0)
Current smoking	18 (1.5)
≥3 Alcoholic drinks per day	17 (1.5)
Arthritis	28 (2.4)
Lumbar spine T-score	−1.289±1.494
Lumbar spine BMD, g/cm2	0.957±0.179
Femur neck T-score	−0.846±1.106
Femur neck BMD, g/cm2	0.800±0.133
L1-4 TBS	1.363±0.097
FRAX probability for MOF without BMD, %	6.4±3.5
FRAX probability for hip fracture without BMD, %	1.9±1.9
FRAX probability for MOF with BMD, %	5.4±2.9
FRAX probability for hip fracture with BMD, %	1.1±1.6
FRAX probability for MOF with BMD and TBS, %	5.4±3.2
FRAX probability for hip fracture with BMD and TBS, %	1.1±1.6

Values are expressed as mean±standard deviation or number (%).

BMI, body mass index; BMD, bone mineral density; L1-4, lumbar spine level 1–4; TBS, trabecular bone score; MOF, major osteoporotic fracture; FRAX, Fracture Risk Assessment Tool.

During the median follow-up period of 7.5 years, 99 (8.5%) women experienced MOFs. The most common fracture site was the wrist (n=31, 2.7%), followed by the hip (n=28, 2.4%), spine (n=21, 1.8%), and humerus (n=19, 1.6%). The women who experienced MOFs were significantly older and were more likely to have a prior history of osteoporotic fracture than those who did not. Furthermore, the women with MOFs showed significantly lower lumbar spine BMD, femoral neck BMD, and TBS than those without MOFs. The FRAX probabilities for MOF and hip fracture were significantly higher in the women who experienced MOFs than in those who did not (Table 2).

Table 2

Comparison of Clinical Risk Factors and Bone-Related Parameters According to the Occurrence of Major Osteoporotic Fractures

Variable	Women without MOF (n=1, 066)	Women with MOF (n=99)	P value
Age, yr	61.3±8.8	64.6±7.5	<0.001
Height, cm	152.4±6.1	151.4±6.2	0.153
Weight, kg	57.9±8.7	57.4±8.9	0.590
BMI, kg/m2	24.9±3.4	25.0±3.4	0.857
Menopause	976 (91.6)	96 (97.0)	0.077
Previous history of osteoporotic fracture	99 (9.3)	24 (24.2)	<0.001
Current smoking	16 (1.5)	2 (2.0)	0.661
≥3 Alcoholic drinks per day	16 (1.5)	1 (1.0)	1.000
Arthritis	25 (2.3)	3 (3.0)	0.726
Lumbar spine T-score	−1.249±1.498	−1.723±1.381	0.003
Lumbar spine BMD, g/cm2	0.961±0.180	0.904±0.166	0.003
Femoral neck T-score	−0.817±1.117	−1.162±0.936	0.003
Femoral neck BMD, g/cm2	0.802±0.134	0.760±0.112	0.025
L1-4 TBS	1.366±0.097	1.332±0.088	0.001
FRAX probability for MOF without BMD, %	6.3±3.4	7.9±3.7	<0.001
FRAX probability for hip fracture without BMD, %	1.8±1.9	2.5±1.9	0.001
FRAX probability for MOF with BMD, %	5.3±2.9	6.4±3.0	<0.001
FRAX probability for hip fracture with BMD, %	1.1±1.5	1.5±1.8	0.014
FRAX probability for MOF with BMD and TBS, %	5.3±3.2	6.5±3.3	<0.001
FRAX probability for hip fracture with BMD and TBS, %	1.1±1.6	1.5±1.8	0.012

Values are expressed as mean±standard deviation or number (%).

MOF, major osteoporotic fracture; BMI, body mass index; BMD, bone mineral density; L1-4, lumbar spine level 1–4; TBS, trabecular bone score; FRAX, Fracture Risk Assessment Tool.

In Cox proportional hazard regression models, the significant risk factors for MOFs were age, prior osteoporotic fracture, lumbar spine BMD, femur neck BMD, and lumbar spine TBS. Among these, a history of a prior osteoporotic fracture was the strongest risk factor for MOFs (HR, 3.09; 95% CI, 1.69 to 5.64). The risk for MOFs significantly increased by 63% as age increased by 10 years. All types of FRAX probabilities were shown to be related to the risk of MOFs (Table 3). However, in the multivariate analysis, only a history of a prior osteoporotic fracture (HR, 2.30; 95% CI, 1.22 to 4.36) and the BMD-adjusted FRAX probability (HR, 1.08; 95% CI, 1.01 to 1.15) remained significant risk factors for MOFs (Table 3).

Table 3

Cox Regression Models of Age, BMI, Previous History of Fracture, BMD, TBS, and FRAX Probability for Major Osteoporotic Fracture

Variable	Unadjusted	Adjusted
Age (per 10-yr increase)	1.63 (1.19–2.24)	1.39 (0.87–2.22)
BMI (per SD decrease)	0.95 (0.89–1.03)	0.94 (0.87–1.02)
Menopause	1.67 (0.52–5.33)	0.67 (0.17–2.57)
Previous osteoporotic fracture	3.09 (1.69–5.64)	2.30 (1.22–4.36)
LS BMD (per SD decrease)	1.20 (1.00–1.44)	1.05 (0.79–1.38)
FN BMD (per SD decrease)	1.32 (1.03–1.69)	0.99 (0.66–1.47)
L1-4 TBS (per SD decrease)	1.43 (1.11–1.84)	1.16 (0.83–1.62)
FRAX probability without BMD, %	1.10 (1.04–1.17)	1.08 (1.01–1.15)
FRAX probability with BMD, %	1.10 (1.03–1.18)	1.06 (0.98–1.15)
FRAX probability with BMD and TBS, %	1.09 (1.02–1.15)	NA

Values are expressed as hazard ratio (95% confidence interval). Clinical risk factors (age, BMI, menopause, and previous osteoporotic fracture), LS BMD, FN BMD, and L1-4 TBS were adjusted by each of the other factors. FRAX probability without BMD and FRAX probability with BMD were adjusted by TBS. Unadjusted hazard ratios were calculated by univariate Cox proportional hazard model analysis. Adjusted hazard ratios were calculated by multivariate Cox proportional hazard model analysis.

BMI, body mass index; BMD, bone mineral density; TBS, trabecular bone score; FRAX, Fracture Risk Assessment Tool; SD, standard deviation; LS, lumbar spine; FN, femoral neck; L1-4, lumbar spine level 1–4; NA, not available.

No variables were found to significantly elevate the risk of hip fracture owing to the low incidence of hip fractures. The HRs (95% CIs) of FRAX probabilities without BMD, with BMD, and with BMD adjusted by TBS for hip fractures were 1.00 (95% CI, 0.99 to 1.02), 1.01 (95% CI, 0.99 to 1.03), and 1.01 (95% CI, 0.99 to 1.03), respectively (data not shown).

We compared the predictive power of each model by using the Harrell’s C statistic (Table 4). The FRAX probability with BMD and the TBS-adjusted FRAX probability for MOFs did not predict MOFs better than the FRAX probability without BMD. As the FRAX probabilities were very low in this study, we calculated their cut-off values by using the maxstat method for MOF prediction. The cut-off values of FRAX probabilities without BMD, with BMD, and with BMD adjusted by TBS were 7.2%, 5.0%, and 6.7%, respectively (all P<0.001).

Table 4

Harrell’s C Statistics and Cut-off Values of FRAX Probability without BMD, with BMD, and with BMD and TBS for Predicting Major Osteoporotic Fractures

Variable	Harrell C	Cut-off value, %
FRAX probability without BMD	0.717a	7.2
FRAX probability with BMD	0.711b	5.0
FRAX probability with BMD and TBS	0.697c	6.7

FRAX, Fracture Risk Assessment Tool; BMD, bone mineral density; TBS, trabecular bone score.

P value between a and b = 0.356; P value between b and c=0.137; P value between c and a=0.936.

DISCUSSION

In the present study, we demonstrated that clinical risk factors such as age and prior osteoporotic fractures were significant predictors of MOFs. However, BMD- or TBS-adjusted FRAX probabilities did not show improved predictive value for MOFs compared with the clinical risk factor-based FRAX probability alone in community-dwelling Korean women. In addition, the cut-off value of the FRAX probability without BMD for predicting MOFs was 7.2%, which is much lower than the FRAX probability for MOFs of 20% that is generally used as a threshold for osteoporosis treatment in Western countries.

The present study demonstrated that the FRAX probability without BMD for MOFs in women was comparable to the BMD- or TBS-adjusted FRAX probability. In previous studies assessing the utility of FRAX, the FRAX probability with BMD predicted osteoporotic fractures as well as a multiple logistic regression model combining clinical risk factors including age, body weight, femoral neck BMD, and prior fracture history in women [29-31]. However, most studies that showed good predictive values of FRAX for osteoporotic fractures were performed in Caucasian populations [32-34]. In fact, insufficient data have been published for validating the predictive power of FRAX probabilities in the Korean population. Despite the limited evidence from previous Korean studies regarding the feasibility of FRAX, prior research has indicated that FRAX probabilities in Koreans underestimated the risk of MOFs compared with those in other ethnicities [25,26,33]. Our results imply that FRAX probabilities for Koreans should be readjusted using appropriately weighted risk factors influencing osteoporotic fractures.

Previous studies have shown that TBS may be of additional value in predicting osteoporotic fractures. A Manitoba cohort study demonstrated that spine TBS and BMD were independent risk factors for osteoporotic fractures, and that the combination of these two measurements was superior to either measurement alone in postmenopausal women [14]. According to another study from the United States that validated the predictive value of FRAX in older men, TBS was independently associated with incident major osteoporotic and hip fractures, and the predictive power of FRAX with BMD after TBS adjustment was significantly improved [35]. Moreover, a meta-analysis of studies including 17,809 men and women from 14 prospective population-based cohorts showed that the gradient of risk (GR) of TBS for MOFs was 1.44 (95% CI, 1.35 to 1.53) when adjusted for age and time since baseline. After an additional adjustment for FRAX probability with BMD, TBS remained a significant independent predictor of osteoporotic fractures (GR, 1.32; 95% CI, 1.24 to 1.41). Adjusting the FRAX probability by TBS, derived from McCloskey et al. [34], resulted in a small increase in the GR (1.76 [95% CI, 1.65 to 1.87] vs. 1.70 [95% CI, 1.60 to 1.81]) [15]. However, we demonstrated that a low TBS increased the risk for osteoporotic fractures, but failed to show an additional prognostic role of TBS-adjusted FRAX in predicting MOFs in Korean women. Ripamonti et al. [36] reported that TBS was not associated with spine fragility fractures in postmenopausal women with osteoporosis, although TBS—but not spine BMD—was a significant risk factor for spine fragility fractures in postmenopausal women without osteoporosis. In addition, Mirzaei et al. [37] compared the predictive value of three different models using BMD, TBS, and the combination of BMD and TBS for spine fragility fractures in postmenopausal women, and they showed no clinical benefit of TBS in postmenopausal women. Adding TBS to BMD or FRAX neither improved the predictive value of BMD for vertebral fractures nor changed the treatment decisions made based on FRAX [37]. In this study, no additional benefit of TBS adjustment was found for the predictive power of the FRAX with BMD model for MOFs and hip fractures. These results could be attributed to the distinctly lower FRAX probabilities in Koreans than in Caucasian populations or even in Japanese people. The narrow distribution of FRAX probabilities in this study can also explain the small value of TBS-adjusted FRAX for fracture risk assessment. No definitive evidence exists to explain the low FRAX probabilities recorded in this study; however, relatively young subjects were included and both the BMD and TBS values were high in subjects regardless of whether they experienced MOFs. Moreover, the clinical risk factors assessed through the self-reported questionnaires might be causally linked to the low FRAX probabilities. Consequentially, we also observed that FRAX probabilities were attenuated after BMD and TBS adjustment in our study.

In the present study, the cut-off value of FRAX probabilities for MOF was 7.2% and that for hip fracture could not be calculated due to the low incidence of hip fracture. The cut-off point of FRAX 10-year probability requiring treatment in United States is 20% for MOF and 3% for hip fracture. In Japan, the cut-off point used as an intervention threshold is 15% for MOF and 3% for hip fracture. Compared to these values, our study showed very low cut-off points for MOF. A recent study of the FRAX-based intervention threshold in seven Latin American counties reported that the cut-off points of FRAX calculated with BMD according to age varied from country to country. The cut-off values of FRAX for MOF ranged from 1.5% to 27.5% in Argentina, from 3.8% to 25.2% in Brazil, from 1.6% to 20.0% in Chile, from 0.6% to 10.2% in Colombia, from 0.9% to 13.6% in Ecuador, from 2.6% to 20.0% in Mexico, and from 0.7% to 22.0% in Venezuela at the ages of 40 and 90 years, respectively [24]. As previously described, our study subjects were younger (with ages ranging from 45 to 76 years) than the Latin American study subjects, and we did not calculate cut-off points for each age group. Moreover, the low FRAX probabilities, which may have been caused by the possibility of underestimating clinical risk factors, and the relatively high BMD and TBS values in the Korean population might lead to a relatively low intervention threshold of FRAX probabilities.

To the best of our knowledge, this is the first large-prospective cohort study conducted in Koreans for almost a decade to verify the ability of TBS-adjusted FRAX to predict osteoporotic fracture. Another strength of this study is that it suggested cut-off values for each FRAX probability as intervention thresholds.

There are also several limitations of our study. First, its short follow-up period of less than 10 years could be a limitation because FRAX algorithms calculate 10-year osteoporotic fracture probability. As a result, the incidence of MOF might have been unrealistically low in this study. Second, the total number of study subjects was too small to support generalized conclusions. Third, we excluded men (n=907) due to the extremely low incidence of MOFs in men (n=36, 4.0%), which might have led to selection bias. Fourth, another reason why the incidence of MOFs was low in our study might have been that cases of morphometric vertebral compression fractures with no directly relevant medical history were excluded due to the unclear onset of fractures. This exclusion might have caused the incidence of osteoporotic fractures to be underestimated. Lastly, the lack of specific clinical information about underlying illnesses affecting bone health due to inaccuracies in the questionnaires might cause bias. For instance, it was not possible to discriminate rheumatoid arthritis from osteoarthritis. Further large prospective cohort studies with long-term follow-up are warranted to confirm the feasibility of FRAX in the Korean population.

In conclusion, the FRAX probability with BMD adjusted for TBS did not predict MOFs and hip fracture better than the FRAX probability with BMD or the FRAX probability without BMD in Korean women. Lower cut-off values of FRAX probabilities appear to be needed to determine whether osteoporosis treatment is indicated in Korean women.