Introduction
Osteoporotic vertebral compression fracture (OVCF) is a common and clinically important manifestation of osteoporosis in older adults and is associated with pain, disability, progressive spinal deformity, impaired mobility, and reduced quality of life [
1,
2]. In spine practice, OVCF represents a fracture event and a marker of vertebral skeletal fragility that influences sagittal alignment, functional status, treatment decisions, and future fracture risks, even when vertebral fractures are clinically silent or asymptomatic [
3,
4].
Bone mineral density (BMD) measured by dual-energy X-ray absorptiometry (DXA) is the current reference standard for osteoporosis diagnosis based on the World Health Organization T-score threshold of ≤−2.5 [
5]. However, BMD explains only 60%–70% of bone strength, as vertebral strength also depends on bone quality, including trabecular microarchitecture, geometry, and tissue material properties [
6]. This limitation is particularly relevant in vertebral compression fractures because many fragility fractures occur despite nonosteoporotic BMD, and prevalent vertebral fractures are frequently observed among patients with clinical risk factors and nonosteoporotic lumbar BMD [
7].
Lumbar spine DXA interpretation in older adults is further complicated by degenerative changes, osteophytes, vertebral deformity, and overlying abdominal aortic calcification (AAC), which may artifactually elevate lumbar BMD and mask spinal fragility [
8,
9]. This diagnostic gap is clinically important because patients with radiographic OVCF may be classified as nonosteoporotic using conventional BMD thresholds, potentially delaying the recognition of bone fragility and appropriate osteoporosis care.
Trabecular bone score (TBS), a texture-based index derived from lumbar spine DXA images, provides complementary information on trabecular microarchitectural integrity not captured by BMD [
6]. Low TBS indicates degraded bone microarchitecture and predicts osteoporotic fracture risk independently of BMD [
10]. TBS has also been integrated into fracture-risk assessment models, including the Fracture Risk Assessment Tool (FRAX), to enhance future risk stratification beyond BMD alone [
11]. In patients with nonosteoporotic BMD, TBS may help detect prevalent vertebral fractures and occult vertebral fragility that may be missed by BMD-based assessment [
7,
9].
AAC is a radiographic marker of vascular calcification that may reflect the systemic bone–vascular interplay underlying skeletal fragility [
12,
13]. AAC can be readily identified on lateral lumbar radiographs using validated semiquantitative scoring methods [
12]. A greater AAC burden has been linked to lower BMD and elevated vertebral and nonvertebral fracture risk and lower TBS independent of BMD [
14,
15]. However, whether AAC burden provides additional clinical insight into BMD–TBS discordance and radiographically confirmed OVCF remains uncertain.
Although BMD, TBS, and AAC have each been investigated in relation to fracture risk, few studies have evaluated these measures together within a matched framework focused on radiographically confirmed OVCF. The present study addresses this gap by concurrently assessing the lumbar spine, femoral neck, and total hip BMD, TBS, and high-burden AAC in an Asian spine-practice cohort. This approach enabled us to evaluate whether TBS and high-burden AAC provide additional clinical context beyond conventional DXA-derived BMD assessment.
This study aimed to examine the relationship between TBS and OVCF after accounting for BMD at conventional DXA sites, including the lumbar spine, femoral neck, and total hip. Secondary objectives were to investigate BMD–TBS discordance, defined as nonosteoporotic BMD at all measured DXA sites with low TBS, and to evaluate whether high-burden AAC was associated with OVCF after adjusting for BMD and TBS.
Materials and Methods
Study design and participants
This retrospective matched case-control study was performed at Chulabhorn Hospital, Thailand, in accordance with the Declaration of Helsinki and was approved by the Human Research Ethics Committee, Chulabhorn Research Institute (IRB No. 055/2569). Informed consent was waived because only existing de-identified clinical and imaging data obtained during routine care were analyzed.
The study population consisted of patients aged ≥50 years who underwent DXA with lumbar spine and hip BMD measurements, with TBS derived from the same lumbar spine DXA scan, and were evaluated with available vertebral imaging between January 2025 and March 2026. Vertebral imaging included vertebral fracture assessment (VFA) or lateral lumbar spine radiography. The study was conducted and reported in accordance with the STROBE (Strengthening the Reporting of Observational Studies in Epidemiology) guidelines.
Case and control definitions
Cases were defined as patients with prevalent radiographic OVCF, defined as a Genant grade 2 or 3 vertebral compression fracture on VFA or lateral spine radiography in individuals aged ≥50 years, after excluding those with documented bone malignancy. OVCF was considered a radiographic marker of vertebral fragility rather than a BMD-based diagnosis. Grade 1 vertebral deformities were not used to define the case status.
Controls were selected from patients without radiographic vertebral fracture on VFA or lateral spine radiography and were matched to cases by sex and age (within ±5 years). All patients were required to have undergone complete lumbar spine BMD, TBS, and vertebral imaging needed for AAC assessment. Femoral neck and total hip BMD were additionally required for primary cohort analyses integrating hip BMD.
Participants with nonevaluable vertebral imaging, missing TBS, spinal instrumentation, or major artifacts interfering with BMD, TBS, or AAC assessment, documented metastatic or bone malignancy, or end-stage renal disease requiring hemodialysis, were excluded. Documented active osteoporosis treatment prior to the index DXA scan was reviewed when available; however, since treatment status could not be consistently verified, it was not employed as a strict exclusion criterion.
Matching procedure and analytic cohorts
Eligible controls were identified from the DXA and vertebral imaging database. For the primary analytic cohort, only OVCF cases with at least four eligible sex- and age-matched controls were included. When more than four controls were available, the four with the smallest absolute age difference were selected. The primary cohort included 47 OVCF cases and 188 controls using a strict 1:4 matching ratio.
A sensitivity cohort was subsequently assembled using all eligible OVCF cases and matched controls, yielding 51 OVCF cases and 223 controls. Because femoral neck and total hip BMD were available only for the primary cohort, analyses incorporating hip BMD, including the lowest all-site BMD T-score and any-site osteoporosis, were limited to the primary cohort.
DXA-derived BMD and TBS assessment
All DXA examinations were performed using a GE Lunar iDXA densitometer (GE Healthcare, Madison, WI, USA). Lumbar spine BMD was obtained from the posteroanterior L1–L4 DXA. Following the International Society for Clinical Densitometry recommendations, vertebrae with compression fractures, implants, major artifacts, local structural abnormalities, or a T-score difference >1.0 from adjacent vertebrae were excluded from analysis. BMD classification was not based on any single vertebra.
Femoral neck and total hip BMD were measured using the same index DXA. The lowest all-site BMD T-score was defined as the lowest among the lumbar spine, femoral neck, and total hip. TBS was derived from the index lumbar DXA acquisition and calculated from the same evaluable vertebral levels used for lumbar spine BMD analysis.
Vertebral fracture and AAC assessment
Vertebral fracture status was determined using VFA images or lateral lumbar spine radiographs. When both imaging modalities were available, lateral radiographs were preferentially used for vertebral fracture confirmation and AAC assessment. Vertebral fractures were graded independently by two orthopedic reviewers using the semiquantitative Genant method. Briefly, grade 1 represents a mild deformity (<25% vertebral height loss), grade 2 a moderate fracture (25%–40% height loss), and grade 3 a severe fracture (>40% height loss) [
16]. Only grade 2 or 3 deformities were classified as OVCF, and disagreements were resolved by consensus. For exploratory fracture severity analyses, OVCF cases were further characterized according to the maximum Genant grade, number of fractured vertebral levels, presence of multiple vertebral fractures, and a Genant-based fracture-burden score. The Genant-based fracture-burden score was calculated as the sum of Genant grades across all fractured vertebral levels; for example, a single grade 3 fracture yielded a score of 3, whereas graded 2 and 3 fractures at two vertebral levels yielded a score of 5. Inter-rater reliability for Genant grading and the Genant-based fracture-burden score was evaluated using weighted kappa and intraclass correlation coefficient (ICC), respectively.
AAC was assessed on lateral vertebral imaging using the Kauppila 24-point semiquantitative scoring system [
12]. Scores were assigned independently by two orthopedic reviewers blinded to BMD and TBS values; however, complete blinding to vertebral fracture status was not always possible because fracture grading and AAC scoring were performed on the same lateral images. The final AAC-24 score was defined as the mean of the two reviewers’ scores when the absolute difference was ≤3 points, whereas larger differences were resolved by consensus. Inter-rater reliability for AAC-24 was evaluated using the ICC for absolute agreement. AAC-24 was analyzed as a continuous and a categorical variable. High-burden AAC was prespecified as an AAC-24 score ≥9, representing extensive calcification burden in prior AAC-24 studies [
17,
18], and was not intended to represent a universal clinical treatment threshold.
Phenotype definitions
Low TBS was defined as TBS ≤1.230, corresponding to the commonly used degraded trabecular microarchitecture in TBS classification studies [
19–
21]. As population-specific TBS thresholds and fracture-based validation for Thai patients and broader Asian spine-practice cohorts remain limited, this threshold was applied uniformly to all participants. Any-site osteoporosis was defined as a T-score ≤−2.5 at the lumbar spine, femoral neck, or total hip. Nonosteoporotic BMD was defined as T-scores >−2.5 at all measured DXA sites. BMD–TBS discordance was defined as nonosteoporotic BMD at all measured DXA sites in the presence of low TBS. An exploratory combined phenotype was defined by the coexistence of nonosteoporotic BMD, low TBS, and high-burden AAC. This phenotype was evaluated descriptively and considered hypothesis-generating because it represented a secondary exploratory subgroup rather than a predefined primary exposure. Exploratory fracture-burden variables were defined as described above.
Statistical analysis
Baseline characteristics are presented as mean±standard deviation or number with percentage. Formal between-group comparisons were not emphasized because inferential analyses accounted for the matched study design.
Matched conditional logistic regression was employed to evaluate associations with OVCF in the primary 1:4 matched cohort, with models stratified by matched case-control sets. Continuous BMD and TBS variables were standardized and expressed per 1 standard deviation increase. Regression models examined lumbar spine BMD and TBS; the lowest all-site BMD T-score and TBS; and clinical-threshold variables, including low TBS, high-burden AAC, and BMD represented either as the continuous lowest all-site BMD T-score or as the conventional any-site osteoporosis threshold. The clinical-threshold model incorporating the continuous lowest all-site BMD T-score was prioritized as the principal model because it retained densitometric information across conventional DXA sites, whereas the model using any-site osteoporosis served as a threshold-based clinical comparator. Odds ratios (ORs)and 95% confidence intervals (CIs) were reported.
The sample size was determined by the number of eligible patients with radiographically confirmed OVCF and complete DXA, TBS, AAC, and hip BMD data during the study period; no formal a priori sample size calculation was performed. To evaluate statistical precision, a post hoc Wald-based precision assessment was conducted using the 95% CI of the low-TBS estimate from the principal clinical-threshold model.
Subgroup analyses described the prevalence of low TBS, high-burden AAC, and the combined phenotype among participants with nonosteoporotic BMD at all measured DXA sites. The combined exploratory phenotype was analyzed descriptively. Sensitivity analyses repeated the primary lumbar spine BMD-based regression models in the all-valid matched cohort. Exploratory case-only analyses examined associations between imaging markers and Genant-based fracture burden using the Mann-Whitney U test, Fisher’s exact test, and Spearman rank correlation, as appropriate. A two-sided p-value <0.05 was considered statistically significant. Analyses were performed using R ver. 4.3.3 (The R Foundation for Statistical Computing, Vienna, Austria) and the survival package (ver. 3.5.8).
Discussion
This matched case-control study provides clinically oriented relevant evidence on BMD–TBS discordance in radiographically confirmed OVCF within an Asian spine-practice cohort. Building on prior evidence linking TBS and AAC to fracture risk, we evaluated lumbar spine, femoral neck, and total hip BMD, TBS, and high-burden AAC within a single matched framework focused on Genant-confirmed OVCF. Three findings are most relevant. First, low TBS remained independently associated with OVCF after accounting for the lowest all-site BMD T-score. Second, more than half of OVCF cases did not meet the osteoporotic threshold at any measured DXA site, and low TBS was approximately twice as common among these patients as in matched controls. Third, although high-burden AAC was more prevalent among OVCF cases, its association was attenuated after accounting for all-site BMD and low TBS, supporting its interpretation as a contextual vascular marker rather than a standalone skeletal fragility predictor. These findings are relevant to spine practice because radiographic OVCF may reflect vertebral fragility despite nonosteoporotic BMD at conventional DXA sites.
A major challenge in vertebral fracture evaluation is that BMD alone does not fully reflect skeletal fragility. Although a T-score ≤−2.5 remains the diagnostic threshold for osteoporosis, many fragility fractures occur in individuals whose BMD is in the osteopenic or even normal BMD range [
5,
22–
24]. This limitation is particularly relevant in the spine, where degenerative changes, vertebral deformity, and vascular calcification may artifactually elevate areal BMD and obscure underlying vertebral fragility [
6,
9]. The present study extends these observations by demonstrating that this diagnostic gap was not confined to lumbar spine BMD alone: 53.2% of OVCF cases did not meet the osteoporotic threshold even after incorporating femoral neck and total hip BMD. This finding supports the clinical concern that conventional DXA thresholds may fail to identify vertebral fragility in a substantial proportion of patients with radiographically confirmed OVCF.
Hip BMD plays an important role in osteoporosis diagnosis and fracture-risk assessment; however, because femoral neck and total hip BMD are measured at nonvertebral skeletal sites, they may not fully reflect local spinal fragility. TBS, derived from the lumbar DXA texture, offers complementary information by indirectly assessing the trabecular microarchitecture [
6,
25,
26]. Prior cohort studies and meta-analyses have demonstrated that lower TBS predicts vertebral and major osteoporotic fractures independently of BMD and clinical risk factors [
11,
20,
27]. Among patients with nonosteoporotic BMD, TBS has also been shown to improve the detection of prevalent vertebral fractures beyond BMD alone [
7]. The present findings extend these observations to an all-site BMD framework: low TBS remained significantly associated with OVCF after adjusting for the lowest all-site BMD T-score, and this association persisted even when BMD was represented categorically as any-site osteoporosis. These results support TBS as a clinically useful adjunct for identifying vertebral fragility that may not be detected using conventional BMD thresholds [
20,
21]. Methodologically, TBS and lumbar BMD were derived from the same evaluable vertebral levels, minimizing regional mismatch. The association of low TBS in the model using the conventional any-site osteoporosis threshold further supports its clinical relevance when BMD is interpreted categorically. However, the model incorporating the continuous lowest all-site BMD T-score was emphasized because it preserves more densitometric information across conventional DXA sites.
The AAC findings should be interpreted cautiously because the observed associations varied across statistical models. In the primary cohort, high-burden AAC was more common among OVCF cases, but its association with OVCF was attenuated after adjustment for the lowest all-site BMD T-score and low TBS. In contrast, high-burden AAC remained associated with OVCF in the lumbar spine BMD-based sensitivity models. This discrepancy may be explained by several non-mutually exclusive factors. First, the AAC–BMD relationship may be site-dependent. Because AAC is anatomically adjacent to the lumbar spine, it may influence lumbar DXA interpretation [
9,
14], whereas femoral neck and total hip BMD provide non-spinal skeletal information that is less affected by local aortic calcification. Second, lumbar spine-only models may not fully account for skeletal fragility because lumbar BMD can be artifactually increased by degenerative change, vertebral deformity, or vascular calcification [
6,
9]. Incorporating hip BMD may therefore minimize residual confounding and attenuate the apparent AAC association. Third, the relatively smaller number of OVCF cases may have contributed to sampling variability and limited precision across model specifications. Taken together, these findings suggest that high-burden AAC should not be interpreted as a direct measure of bone quality or as a standalone skeletal predictor in this cohort.
This interpretation is consistent with the bone–vascular axis concept, which proposes that vascular calcification and skeletal deterioration may share overlapping biological pathways [
13]. Previous observational evidence has linked greater AAC burden to lower BMD and increased vertebral and nonvertebral fracture risk [
14], while aortic vascular calcification has been inversely associated with TBS in patients undergoing dialysis [
15]. In the present hospital-based cohort of non-dialysis patients, AAC was best interpreted as complementary vascular information that may coexist with skeletal fragility and degraded trabecular microarchitecture, rather than as a universal threshold for severe AAC or clinical intervention.
Methodologically, OVCF was conservatively defined using Genant grade 2 or 3 fractures rather than by diagnostic codes or mild grade 1 deformities, which have lower reproducibility [
16,
28,
29]. In exploratory case-only analyses, the Genant-based fracture-burden score was not significantly associated with BMD, TBS, AAC-24 score, low TBS, high-burden AAC, or BMD–TBS discordance. Because these analyses were restricted to patients with confirmed OVCF and were not powered to evaluate fracture-burden gradients, they should be interpreted descriptively. Prior evidence also indicates that TBS and VFA may capture complementary skeletal risk dimensions [
30].
This study has several strengths, including sex- and age-matched design, radiographic confirmation of OVCF, inclusion of a hospital-based Asian cohort, good-to-excellent inter-rater reliability, and alignment of TBS analysis with the same evaluable vertebral levels used for lumbar BMD. Nevertheless, several limitations should be acknowledged. First, the sample size was determined by the number of eligible OVCF cases with complete imaging and DXA-derived data during the study period, and no formal a priori sample size calculation was performed. A post hoc Wald-based precision assessment indicated that the study was primarily powered to detect relatively large associations, corresponding to an OR of approximately 2.9 or greater for a binary exposure with precision comparable to that of low TBS. Consequently, smaller independent associations may have been missed, particularly for high-burden AAC, interaction testing, and the combined exploratory phenotype. Second, the retrospective case-control design and the timing of DXA, TBS, and AAC assessments preclude causal or temporal inferences. Accordingly, low TBS observed in OVCF patients should not be interpreted as evidence that degraded trabecular microarchitecture necessarily preceded fracture occurrence. Third, osteoporosis treatment before DXA assessment was not excluded and may have influenced BMD and TBS differently. Prior antiresorptive or anabolic therapy may preserve or increase BMD despite persistent microarchitectural degradation, potentially increasing apparent BMD–TBS discordance. Conversely, effective treatment could mitigate fracture risk and attenuate associations with OVCF. Fourth, residual confounding by body mass index, comorbidities, medication exposure, fall risk, and other clinical fracture-risk factors cannot be excluded. Finally, the TBS cutoff of ≤1.230 was based on widely used international classification thresholds, but population-specific normative data and fracture-based validation for Thai patients and broader Asian spine-practice cohorts remain limited. In addition, the single-center, hospital-based design may limit the generalizability of the findings.
Clinically, these findings support the integrated assessment of BMD, TBS, and vertebral imaging when assessing vertebral fragility in older spine patients, particularly in those with radiographic OVCF despite nonosteoporotic BMD at conventional DXA sites.