The demographic characteristics of the participants
The demographic characteristics of the 9,016 participants based on lumbar spine BMD quartiles are presented in Table 1. Participants in the highest quartile of BMD were more likely to be male, Mexican American, more educated, nondrinkers, and nonsmokers. Among the subgroups, no significant correlation was found between lumbar BMD and physical activity or TC (P > 0.05). Besides, compared to those with lower BMD levels, those with higher BMD levels had higher levels of HbA1c, TG, creatinine, uric acid, BMI, weight, WC, WHTR, and VAI, whereas were less likely to have higher levels of total calcium, phosphorus, and HDL-C (P < 0.001).
Table 1
Baseline characteristics of the research population based on lumbar spine BMD quartiles.
Characteristics
|
Quartiles of lumbar-spine BMD
|
|
|
Q1
(≤ 0.919)
|
Q2
(0.919-1.019)
|
Q3
(1.019-1.115)
|
Q4
(≥ 1.15)
|
P value
|
Age (years)
|
41.16 ± 20.99
|
38.96 ± 17.04
|
37.69 ± 15.79
|
40.22 ± 15.44
|
<0.001
|
Gender
|
|
|
|
|
|
Male
|
1037 (46.42%)
|
1089 (48.16%)
|
1153 (51.20%)
|
1220 (53.73%)
|
<0.001
|
Female
|
1196 (53.58%)
|
1172 (51.84%)
|
1109 (48.80%)
|
1051 (46.27%)
|
|
Race/ethnicity (%)
|
|
|
|
|
<0.001
|
Non-Hispanic white
|
277 (12.41%)
|
264 (11.65%)
|
213 (9.46%)
|
154 (6.77%)
|
|
Non-Hispanic black
|
361 (16.17%)
|
313 (13.84%)
|
233 (10.35%)
|
222 (9.74%)
|
|
Mexican American
|
1461 (65.40%)
|
1508 (66.71%)
|
1556 (69.14%)
|
1512 (66.64%)
|
|
Other race
|
134 (6.02%)
|
176 (7.81%)
|
249 (11.06%)
|
383 (16.85%)
|
|
Level of education (%)
|
|
|
|
|
<0.001
|
Less Than High School Grade
|
339 (15.17%)
|
410 (18.15%)
|
340 (15.11%)
|
311 (13.70%)
|
|
High School Grad
|
410 (18.37%)
|
416 (18.45%)
|
406 (18.05%)
|
506 (22.29%)
|
|
More Than High School Grade
|
891 (39.92%)
|
1073 (47.75%)
|
1184 (52.54%)
|
1266 (55.74%)
|
|
Missing
|
593 (26.54%)
|
362 (15.97%)
|
321 (14.30%)
|
188 (8.28%)
|
|
Alcohol consumption (%)
|
|
|
|
|
<0.001
|
Nondrinker
|
1108 (49.61%)
|
1378 (60.94%)
|
1499 (66.60%)
|
1625 (71.55%)
|
|
Low to moderate drinker
|
454 (20.32%)
|
417 (18.44%)
|
373 (16.55%)
|
362 (15.96%)
|
|
Heavy drinker
|
671 (30.07%)
|
466 (20.62%)
|
379 (16.85%)
|
284 (12.49%)
|
|
Smoking behavior (%)
|
|
|
|
|
<0.001
|
Current
|
968 (43.33%)
|
488 (21.57%)
|
392 (17.41%)
|
233 (10.28%)
|
|
Past
|
556 (24.90%)
|
733 (32.43%)
|
884 (39.26%)
|
969 (42.64%)
|
|
Never
|
709 (31.77%)
|
1040 (46.00%)
|
975 (43.33%)
|
1069 (47.08%)
|
|
Physical activity (%)
|
|
|
|
|
0.156
|
Low
|
344 (15.41%)
|
423 (18.73%)
|
398 (17.66%)
|
440 (19.39%)
|
|
Moderate
|
445 (19.95%)
|
344 (15.20%)
|
449 (19.95%)
|
422 (18.58%)
|
|
High
|
151 (6.74%)
|
126 (5.56%)
|
152 (6.74%)
|
141 (6.19%)
|
|
Not recorded
|
1293 (57.91%)
|
1368 (60.52%)
|
1252 (55.65%)
|
1268 (55.84%)
|
|
HbA1c (%)
|
5.48 ± 0.74
|
5.43 ± 0.76
|
5.45 ± 0.79
|
5.53 ± 0.86
|
<0.001
|
Total cholesterol (mmol/L)
|
4.95 ± 1.17
|
4.97 ± 1.05
|
4.91 ± 1.05
|
4.94 ± 1.02
|
0.348
|
TG (mmol/L)
|
1.54 ± 1.26
|
1.60 ± 1.34
|
1.59 ± 1.45
|
1.76 ± 1.62
|
<0.001
|
Total calcium (mmol/L)
|
2.38 ± 0.09
|
2.37 ± 0.09
|
2.36 ± 0.09
|
2.35 ± 0.09
|
<0.001
|
Creatinine (µmol/L)
|
70.24 ± 28.46
|
74.24 ± 23.80
|
75.80 ± 29.19
|
78.72 ± 27.58
|
<0.001
|
Uric acid (µmol/L)
|
300.34 ± 74.24
|
318.05 ± 80.91
|
319.31 ± 80.24
|
329.82 ± 81.44
|
<0.001
|
Phosphorus (mmol/L)
|
1.31 ± 0.83
|
1.24 ± 0.19
|
1.23 ± 0.19
|
1.21 ± 0.18
|
<0.001
|
HDL-C (mmol/L)
|
1.44 ± 0.43
|
1.38 ± 0.43
|
1.34 ± 0.38
|
1.31 ± 0.39
|
<0.001
|
BMI (kg/m2)
|
24.48 ± 5.27
|
26.36 ± 5.33
|
27.64 ± 5.61
|
29.16 ± 5.52
|
<0.001
|
Weight (kg)
|
66.66 ± 16.95
|
74.55 ± 17.20
|
80.25 ± 17.82
|
86.03 ± 17.57
|
<0.001
|
Height (cm)
|
164.56 ± 9.78
|
167.92 ± 9.52
|
170.27 ± 9.70
|
171.71 ± 17.57
|
<0.001
|
WC (cm)
|
87.21 ± 14.85
|
91.66 ± 14.59
|
94.71 ± 14.88
|
98.28 ± 14.16
|
<0.001
|
WHTR
|
0.53 ± 0.09
|
0.55 ± 0.09
|
0.56 ± 0.09
|
0.57 ± 0.09
|
<0.001
|
VAI
|
3.05 ± 4.13
|
3.61 ± 5.53
|
3.72 ± 6.56
|
4.38 ± 6.23
|
<0.001
|
Mean ± SD for continuous variables: the P value was calculated by the weighted linear regression model. (%) for categorical variables: the P value was calculated by the weighted chi-square test.
BMD, bone density mineral; TG, triglycerides; HDL-C, high-density lipoprotein cholesterol; BMI, body mass index; WC, waist circumference; WHTR, waist-to-height ratio; VAI, visceral adiposity index.
Correlation between VAI and BMD in all parts of the body
The adjusted correlations between VAI and BMD in all parts of the body were presented in Table 2. We found that VAI was significantly positively correlated with BMD levels in the unadjusted model (β = 0.002, 95%CI: 0.002–0.003, P<0.001), and the positive association was still present in model 2 (β = 0.002, 95%CI: 0.002–0.003, P<0.001) and model 3 (β = 0.002, 95%CI: 0.001–0.002, P<0.001) after adjusting for confounders. In the multivariable model 3, every 1 increase in VAI was associated with a 0.002 increase in lumbar BMD. Furthermore, this positive correlation remained significant when the BMD measurement site was changed to the total femur, femoral neck. After transforming VAI into a categorical variable, participants in the highest VAI quartile had a 0.049g/cm2 higher BMD than those in the lowest VAI quartile. And the trend remained significant among different BMD quartile groups (P < 0.01). The subgroup analyses stratified by gender and age are reported in Table 3. In both genders, participants with higher VAI levels had significantly higher BMD than those with lower levels. Meanwhile, BMD was increased with the increase of VAI levels in different age subgroups (P < 0.01).
Table 2
Associations between VAI and BMD in all parts of the body.
|
Lumbar-spine BMD
|
Total femur BMD
|
Femoral neck BMD
|
Model 1
β (95% CI)
|
0.002 (0.002, 0.003) <0.001
|
0.001 (0.001, 0.002)
<0.001
|
0.002 (0.002, 0.003) <0.001
|
Model 2
β (95% CI)
|
0.002 (0.002, 0.003) <0.001
|
0.001 (0.001, 0.002) <0.001
|
0.001 (0.001, 0.002) <0.001
|
Model 3
β (95% CI)
|
0.002 (0.001, 0.002) <0.001
|
0.001 (0.001, 0.002) <0.001
|
0.001 (0.001, 0.002) <0.001
|
Logistic regression models:
Model 1: no covariates were adjusted.
Model 2 was adjusted for demographic factors, including gender, age, race.
Model 3 was adjusted for gender, age, race, education, physical activity, smoking behavior, alcohol consumption, HbA1c, total cholesterol, triglycerides, total calcium, creatinine, phosphorus, uric acid.
Table 3
Subgroup analyses stratified by gender and race.
|
Model 1
β (95% CI)
|
Model 2
β (95% CI)
|
Model 3
β (95% CI)
|
VAI
|
0.002 (0.002, 0.003) <0.001
|
0.002 (0.002, 0.003) <0.001
|
0.002 (0.001, 0.002) <0.001
|
VAI (quartile)
|
Q1
|
Reference
|
Reference
|
Reference
|
Q2
|
0.017 (0.008, 0.025) <0.001
|
0.022 (0.014, 0.031) <0.001
|
0.014 (0.006, 0.021) <0.001
|
Q3
|
0.035 (0.026, 0.043) <0.001
|
0.044 (0.036, 0.053) <0.001
|
0.030 (0.022, 0.038) <0.001
|
Q4
|
0.052 (0.044, 0.061) <0.001
|
0.065 (0.056, 0.074) <0.001
|
0.049 (0.040, 0.057) <0.001
|
P for trend
|
<0.001
|
<0.001
|
<0.001
|
Subgroup analysis stratified by gender
|
Male
|
0.002 (0.001, 0.003) <0.001
|
0.002 (0.001, 0.003) <0.001
|
0.001 (0.000, 0.002) <0.001
|
Female
|
0.002 (0.001, 0.003) <0.001
|
0.003 (0.002, 0.004) <0.001
|
0.003 (0.002, 0.004) <0.001
|
Subgroup analysis stratified by age
|
Age ≤18
|
0.006 (0.002, 0.009)
0.002
|
0.010 (0.006, 0.013)
<0.001
|
0.002 (0.001, 0.002) <0.001
|
Age 19-45
|
0.001 (-0.000, 0.001)
0.122
|
0.001 (0.001, 0.002)
<0.001
|
0.001 (0.000, 0.002)
0.033
|
Age 46-59
|
0.002 (0.001, 0.003) <0.001
|
0.002 (0.001, 0.003) <0.001
|
0.002 (0.001, 0.003) <0.001
|
Age ≥60
|
0.005 (0.003, 0.006) <0.001
|
0.003 (0.001, 0.004) <0.001
|
0.006 (0.001, 0.010)
0.010
|
Logistic regression models:
Model 1: no covariates were adjusted.
Model 2 was adjusted for demographic factors, including gender, age, race.
Model 3 was adjusted for gender, age, race, education, physical activity, smoking behavior, alcohol consumption, HbA1c, total cholesterol, triglycerides, total calcium, creatinine, phosphorus, uric acid.
The subgroup analyses stratified by gender and age
Through multiple linear regression analysis, we demonstrated a significant relationship between VAI and lumbar BMD. Therefore, we used smooth curve fittings and generalized additive models to further investigate the correlation between VAI and lumbar BMD. As shown in Figure 2-4, VAI was positively correlated with lumbar BMD. The correlation between VAI and lumbar BMD was an inverted U-shaped curve in the age ≤18 groups, with a point of inflection at 4.86. As shown in Table 4, for individuals with a VAI <4.86, every 1 unit increase in VAI was associated with a 0.025 g/cm2 higher lumbar BMD (95%CI: 0.015–0.035); In contrast, for individuals with a VAI >4.86, a 1 unit increase in VAI was associated with a 0.027 g/cm2 decrease in lumbar BMD (95%CI: -0.040– -0.014).
Table 4
Threshold effect analysis of VAI on lumbar BMD in the age ≤18 groups using the two-piecewise linear regression model.
Lumbar bone mineral density
|
Adjusted β (95% CI), P value
|
Age ≤18
|
|
Fitting by standard linear model
|
0.006 (0.002, 0.009) 0.002
|
Fitting by two-piecewise linear model
|
|
Inflection point
|
4.86
|
VAI <4.86
|
0.025 (0.015, 0.035) <0.001
|
VAI >4.86
|
-0.027 (-0.040, -0.014) <0.001
|
Log likelihood ratio
|
<0.001
|
Gender, race, education, physical activity, smoking behavior, alcohol consumption, HbA1c, total cholesterol, total calcium, creatinine, phosphorus, and uric acid were adjusted.
VAI, visceral adiposity index; CI, confidence interval.