Association between plasma total homocysteine level within normal range and bone mineral density in adults.

BACKGROUND: Growing evidence indicates that homocysteine is a noteworthy marker for general health status. However, research regarding plasma total homocysteine (tHcy) levels and bone mineral density (BMD) is sparse and controversial. Hence, we aimed to investigate the association between plasma tHcy level within normal range and lumbar BMD in adults.
METHODS: In this cross-sectional study, using the National Health and Nutrition Examination Survey database, data on 10748 adults aged between 30 and 85 years were analyzed. The weighted multiple logistic regression analyses were conducted to evaluate the association between plasma tHcy level and lumbar BMD. The fitted smoothing curves were performed to explore potential non-linear relationships. When non-linearity was detected, we further calculated the inflection point using a recursive algorithm and constructed a weighted two-piecewise linear regression model.
RESULTS: After adjusting for all the covariates, the association between plasma tHcy and lumbar BMD was different in various age groups (adults aged 30-49 years: β = -0.0004, 95% CI -0.0025, 0.0018; adults aged 50-69 years: β = 0.0001, 95% CI -0.0025, 0.0026; adults aged 70-85 years: β = 0.0050, 95% CI 0.0008, 0.0092). In the subgroup analysis stratified by gender, this association also differed based on gender. There was a negative trend in females (aged 30-49 years: β = -0.0022, 95% CI -0.0054, 0.0011; aged 50-69 years: β = -0.0028, 95% CI -0.0062, 0.0007), and a positive trend in males (aged 30-49 years: β = 0.0018, 95% CI -0.0012, 0.0048; aged 50-69 years: β = 0.0027, 95% CI -0.0009, 0.0063) in both 30-49 years group and 50-69 years group. In the 70-85 years group, this association was significantly positive in males (β = 0.0136, 95% CI 0.0068, 0.0204), but was not significantly different in females (β = 0.0007, 95% CI -0.0046, 0.0060).
CONCLUSION: The correlation between plasma tHcy level within the normal range and lumbar BMD differs by age and gender.

Introduction

Homocysteine (Hcy) is a non-essential amino acid derived from methionine metabolism. Growing evidence indicates that Hcy is a noteworthy marker for general health status, and increased plasma total homocysteine (tHcy) level is regarded as an independent risk factor for various human diseases, including cardiovascular diseases, chronic kidney diseases, neurological disorders, and bone tissue damages [1].

In the clinical diagnosis of osteoporosis and prediction of the risk of osteoporotic fracture, bone mineral density (BMD) is a commonly used quantitative indicator [2]. Nevertheless, bone turnover markers (BTMs) are preferable as dynamic indices for evaluating the status of bone turnover and progression after drug administration during the clinical treatment of osteoporosis. Plasma tHcy level, as a bone matrix-related marker, has been mentioned in the clinical guideline for the diagnosis and treatment of osteoporosis [3]. However, the mechanisms underlying the relationship between plasma tHcy level and BMD have not yet been unraveled. Besides, data regarding plasma tHcy level and BMD are relatively sparse, and inconsistent results regarding their relationship, including inverse [4-11], mixed [12-14], and no associations [15-21], have been reported. Therefore, we used the National Health and Nutrition Examination Survey (NHANES) database to investigate the relationship of plasma tHcy level within the normal range with lumbar BMD in adults aged between 30 and 85 years.

Methods

Study participants

The data of this study were obtained from the NHANES (1999–2006), which was an ongoing survey conducted by the Center for Disease Control and Prevention (CDC). CDC used a multistage, complex clustered probability design to select a representative sample of non-institutionalized United States civilians. The survey data are made available on the internet for researchers. A total of 11901 participants aged between 30 and 85 years remained after the exclusion of 1903 subjects with missing plasma tHcy data, 1405 subjects with missing lumbar BMD data, and 1236 subjects with cancer. We further excluded 1153 subjects with plasma tHcy levels outside the normal physiological range of 5–15 μmol/L [22], resulting in a final study population of 10748 people. All protocols were approved by the research ethical review board of the National Center for Health Statistics (NCHS), and informed consent forms were obtained from all participants.

Study variables

The Abbott Imx was used to determine plasma tHcy level for NHANES 1999–2001, and the Abbott AxSym was used starting NHANES 2002. The certified radiology technologists used standard radiologic techniques and protocols to measure and analyze lumbar BMD on a Hologic QDR-4500A fan-beam densitometer. Age, gender, race, income-poverty ratio, education level, alcohol consumption, smoking behavior, physical activity, body mass index, blood urea nitrogen, total protein, total cholesterol, serum uric acid, serum folate, serum vitamin B12, serum phosphorus, and serum calcium were considered potential covariates and confounding factors, which were adjusted in the analytic models. The detailed information on the variables in this study can be found at www.cdc.gov/nchs/nhanes/.

Statistical analysis

The NHANES sample weights were used as recommended by the NCHS. Statistical analyses were performed using R (version 3.5.3) and the EmpowerStats software (http://www.empowerstats.com). P values <0.05 were considered statistically significant. The weighted linear regression models were used to analyze the difference between dichotomous variables, and the weighted chi-square tests were used for continuous variables. The weighted multivariate linear regression analyses were conducted for examining the association between plasma tHcy level and lumbar BMD. The smooth curve fittings were further performed to explore their potential non-linear relationships. When non-linearity was detected, we further calculated the inflection point using a recursive algorithm and constructed a weighted two-piecewise linear regression model.

Results

Table 1 shows the basic characteristics of 10748 participants included in the present study. Compared with the plasma tHcy quartile 1 group, the other quartile groups were older, more likely to be males, had higher blood urea nitrogen, total protein, total cholesterol, serum uric acid, serum phosphorus, and serum calcium, and had lower serum folate and vitamin B12.

Table 1

Description of 10748 participants included in the present study

Plasma total homocysteine quartiles (μmol/L)	Q1(5-6.86)	Q2(6.87-8.22)	Q3(8.23-9.95)	Q4(9.96-15)	P value
Age (years)	44.31 ± 10.71	47.24 ± 11.86	50.15 ± 12.76	55.27 ± 14.31	<0.001
Gender (%)					<0.001
Male	31.07	48.73	61.36	65.57
Female	68.93	51.27	38.64	34.43
Race (%)					<0.001
Non-Hispanic White	66.75	73.54	75.27	76.86
Non-Hispanic Black	10.76	10.33	10.58	10.82
Mexican American	9.35	6.82	5.02	3.95
Other race	13.13	9.31	9.13	8.37
Body mass index (kg/m2)	28.31 ± 6.70	28.86 ± 6.49	28.60 ± 5.90	28.81 ± 6.03	0.004
Education level (%)					<0.001
Less than high school	18.49	17.68	18.73	22.24
High school	22.56	25.37	26.40	28.40
More than high school	58.95	56.96	54.88	49.36
Income-poverty ratio	3.19 ± 1.56	3.27 ± 1.52	3.26 ± 1.54	3.04 ± 1.53	<0.001
Physical activity (%)					<0.001
Not walk very much	14.66	15.03	15.84	19.79
Walk a lot	27.89	27.38	26.34	27.87
Climb often	21.30	20.47	17.30	15.48
Heavy activity	30.60	31.00	32.64	28.28
Not recorded	5.54	6.13	7.88	8.58
Smoking behavior (%)					<0.001
None	57.23	51.12	46.64	41.37
Past	22.45	26.51	28.00	29.13
Current	20.32	22.36	25.35	29.50
Alcohol consumption (%)					<0.001
Non-drinker	33.22	32.40	29.99	36.10
Moderate alcohol use	41.36	35.15	33.15	30.16
High alcohol use	25.42	32.45	36.86	33.74
Blood urea nitrogen (mg/dL)	12.39 ± 3.81	12.99 ± 3.92	13.69 ± 4.32	14.57 ± 5.65	<0.001
Total protein (mg/dL)	7.25 ± 0.46	7.26 ± 0.46	7.29 ± 0.46	7.31 ± 0.50	<0.001
Total cholesterol (mg/dL)	201.1 ± 37.3	205.8 ± 40.2	210.1 ± 42.1	210.4 ± 45.9	<0.001
Serum uric acid (mg/dL)	4.81 ± 1.20	5.30 ± 1.30	5.59 ± 1.32	6.01 ± 1.39	<0.001
Serum folate (ng/mL)	15.75 ± 10.20	14.82 ± 15.76	14.07 ± 11.01	13.04 ± 9.07	<0.001
Serum vitamin B12( pg/mL)	611.6 ± 1758.3	534.8 ± 314.4	490.0 ± 240.1	450.8 ± 299.7	<0.001
Serum phosphorus (mg/dL)	3.63 ± 0.53	3.67 ± 0.56	3.68 ± 0.56	3.68 ± 0.57	0.002
Serum calcium (mg/dL)	9.39 ± 0.35	9.46 ± 0.36	9.52 ± 0.35	9.57 ± 0.38	<0.001
Lumbar bone mineral density (g/cm2)	1.05 ± 0.15	1.04 ± 0.15	1.04 ± 0.16	1.04 ± 0.17	0.054

Mean ± SD for continuous variables: P value was calculated by the weighted linear regression model.

% for categorical variables: P value was calculated by the weighted chi-square test

The results of different multivariate linear regression models for adults aged between 30–49 years, 50–69 years, and 70–85 years are shown in Tables 2, 3, and 4, respectively. After adjusting for all the covariates presented in Table 1, the association between plasma tHcy level and lumbar BMD differed by age (adults aged 30–49 years: β = −0.0004, 95% CI −0.0025, 0.0018; adults aged 50–69 years: β = 0.0001, 95% CI −0.0025, 0.0026; adults aged 70–85 years: β = 0.0050, 95% CI 0.0008, 0.0092). In the subgroup analysis stratified by gender, this association also differed based on gender. There was a negative trend in females (aged 30–49 years: β = −0.0022, 95% CI −0.0054, 0.0011; aged 50–69 years: β = −0.0028, 95% CI −0.0062, 0.0007), and a positive trend in males (aged 30–49 years: β = 0.0018, 95% CI −0.0012, 0.0048; aged 50–69 years: β = 0.0027, 95% CI −0.0009, 0.0063) in both 30–49 years group and 50–69 years group. In the 70–85 years group, this association was significantly positive in males (β = 0.0136, 95% CI 0.0068, 0.0204), but was not significantly different in females (β = 0.0007, 95% CI −0.0046, 0.0060). In the subgroup analysis stratified by race, this association again differed by age and gender (Tables 2, 3, and 4).

Table 2

Association of plasma total homocysteine (umol/L) with lumbar bone mineral density (g/cm2) in adults aged 30–49 years (n = 4977)

	Model 1β (95% CI)	Model 2β (95% CI)	Model 3β (95% CI)
Plasma total homocysteine	−0.0025 (−0.0046, -0.0005)	−0.0012 (−0.0033, 0.0009)	−0.0004 (−0.0025, 0.0018)
Stratified by gender
Male	0.0007 (−0.0022, 0.0036)	−0.0001 (−0.0030, 0.0028)	0.0018 (−0.0012, 0.0048)
Female	-0.0025 (-0.0057, 0.0007)	−0.0025 (−0.0057, 0.0006)	−0.0022 (−0.0054, 0.0011)
Stratified by race
Non-Hispanic White	−0.0045 (−0.0074–0.0016)	−0.0023 (−0.0054, 0.0008)	−0.0007 (−0.0040, 0.0026)
Non-Hispanic Black	0.0018 (−0.0028, 0.0064)	0.0038 (−0.0012, 0.0087)	0.0047 (−0.0005, 0.0100)
Mexican American	-0.0014 (-0.0059, 0.0030)	0.0028 (−0.0019, 0.0076)	0.0015 (−0.0034, 0.0064)
Other race/ethnicity	-0.0040 (-0.0107, 0.0027)	−0.0032 (−0.0106, 0.0041)	−0.0009 (−0.0088, 0.0071)

Model 1: no covariates were adjusted.

Model 2: age, gender, race were adjusted.

Model 3: age, gender, race, body mass index, education level, income-poverty ratio, physical activity, smoking behavior, alcohol consumption, blood urea nitrogen, total protein, total cholesterol, serum uric acid, serum folate, serum vitamin B12, serum phosphorus, and serum calcium were adjusted.

In the subgroup analysis stratified by gender or race, the model is not adjusted for the stratification variable itself

Table 3

Association of plasma total homocysteine (umol/L) with lumbar bone mineral density (g/cm2) in adults aged 50–69 years (n = 4123)

	Model 1β (95% CI)	Model 2β (95% CI)	Model 3β (95% CI)
Plasma total homocysteine	0.0040 (0.0016, 0.0063)	0.0002 (−0.0022, 0.0026)	0.0001 (−0.0025, 0.0026)
Stratified by gender
Male	0.0048 (0.0013, 0.0082)	0.0025 (−0.0009, 0.0059)	0.0027 (−0.0009, 0.0063)
Female	−0.0026 (−0.0059, 0.0007)	−0.0017 (−0.0050, 0.0016)	−0.0028 (−0.0062, 0.0007)
Stratified by race
Non-Hispanic White	0.0025 (−0.0008, 0.0058)	−0.0002 (−0.0036, 0.0032)	−0.0003 (−0.0040, 0.0034)
Non-Hispanic Black	0.0040 (−0.0013, 0.0093)	0.0005 (−0.0049, 0.0059)	0.0007 (−0.0052, 0.0066)
Mexican American	0.0050 (−0.0001, 0.0101)	0.0008 (−0.0045, 0.0062)	0.0018 (−0.0038, 0.0075)
Other race/ethnicity	0.0032 (−0.0050, 0.0115)	0.0028 (−0.0061, 0.0116)	0.0021 (−0.0073, 0.0115)

Model 1: no covariates were adjusted.

Model 2: age, gender, race/ethnicity were adjusted.

Model 3: age, gender, race, body mass index, education level, income-poverty ratio, physical activity, smoking behavior, alcohol consumption, blood urea nitrogen, total protein, total cholesterol, serum uric acid, serum folate, serum vitamin B12, serum phosphorus, and serum calcium were adjusted.

In the subgroup analysis stratified by gender or race, the model is not adjusted for the stratification variable itself

Table 4

Association of plasma total homocysteine (umol/L) with lumbar bone mineral density (g/cm2) in adults aged 70-85 years (n = 1648).

	Model 1β (95% CI)	Model 2β (95% CI)	Model 3β (95% CI)
Plasma total homocysteine	0.0103 (0.0061, 0.0144)	0.0053 (0.0013, 0.0093)	0.0050 (0.0008, 0.0092)
Stratified by gender
Male	0.0107 (0.0043, 0.0170)	0.0096 (0.0031, 0.0162)	0.0136 (0.0068, 0.0204)
Female	0.0015 (-0.0035, 0.0065)	0.0024 (-0.0027, 0.0075)	0.0007 (-0.0046, 0.0060)
Stratified by race
Non-Hispanic White	0.0097 (0.0046, 0.0149)	0.0050 (-0.0001, 0.0100)	0.0047 (-0.0006, 0.0100)
Non-Hispanic Black	0.0126 (0.0001, 0.0251)	0.0179 (0.0050, 0.0308)	0.0139 (0.0004, 0.0273)
Mexican American	0.0002 (-0.0160, 0.0163)	-0.0040 (-0.0228, 0.0149)	0.0040 (-0.0190, 0.0270)
Other race/ethnicity	0.0095 (0.0054, 0.0136)	0.0053 (0.0013, 0.0093)	0.0050 (0.0008, 0.0092)

Model 1: no covariates were adjusted.

Model 2: age, gender, race/ethnicity were adjusted.

Model 3: age, gender, race, body mass index, education level, income-poverty ratio, physical activity, smoking behavior, alcohol consumption, blood urea nitrogen, total protein, total cholesterol, serum uric acid, serum folate, serum vitamin B12, serum phosphorus, and serum calcium were adjusted.

In the subgroup analysis stratified by gender or race, the model is not adjusted for the stratification variable itself.

We further explored the potential non-linear relationship of plasma tHcy level and lumbar BMD using smooth curve fittings (Figs. 1, 2 and 3). Non-linear relationships of plasma tHcy levels with lumbar BMD in males aged 30–49 years and 50–69 years were detected. We further calculated the inflection points at 9.0 μmol/L for males aged 30–49 years and 10.0 μmol/L for 50–69 years (Table 5).

Fig. 1

The association between plasma total homocysteine and lumbar bone mineral density in adults aged 30–49 years, stratified by gender. Age, race, body mass index, education level, income-poverty ratio, physical activity, smoking behavior, alcohol consumption, blood urea nitrogen, total protein, total cholesterol, serum uric acid, serum folate, serum vitamin B12, serum phosphorus, and serum calcium were adjusted

Fig. 2.

The association between plasma total homocysteine and lumbar bone mineral density in adults aged 50–69 years, stratified by gender. Age, race, body mass index, education level, income-poverty ratio, physical activity, smoking behavior, alcohol consumption, blood urea nitrogen, total protein, total cholesterol, serum uric acid, serum folate, serum vitamin B12, serum phosphorus, and serum calcium were adjusted

Fig. 3

The association between plasma total homocysteine and lumbar bone mineral density in adults aged 70–85 years, stratified by gender. Age, race, body mass index, education level, income-poverty ratio, physical activity, smoking behavior, alcohol consumption, blood urea nitrogen, total protein, total cholesterol, serum uric acid, serum folate, serum vitamin B12, serum phosphorus, and serum calcium were adjusted

Table 5

Threshold effect analysis of plasma total homocysteine on lumbar bone mineral density using a two-piecewise linear regression model

Lumbar bone mineral density	Adjusted ß (95% CI)
Plasma total homocysteine
Males aged 30-49 years
Fitting by standard linear model	0.0018 −0.0012, 0.0048)
Fitting by two-piecewise linear model
Inflection point	9.0
Plasma total homocysteine <9.0 (umol/L)	0.0042 (−0.0015, 0.0099)
Plasma total homocysteine >9.0 (umol/L)	-0.0005 (−.0059, 0.0050)
Log-likelihood ratio	0.323
Males aged 50-69 years
Fitting by standard linear model	0.0027 (−0.0009, 0.0063)
Fitting by two-piecewise linear model
Inflection point	10.0
Plasma total homocysteine <10.0 (umol/L)	0.0098 (0.0032, 0.0163)
Plasma total homocysteine >10.0 (umol/L)	−0.0053 −0.0125, 0.0019)
Log-likelihood ratio	0.011

Age, race, body mass index, education level, income-poverty ratio, physical activity, smoking behavior, alcohol consumption, blood urea nitrogen, total protein, total cholesterol, serum uric acid, serum folate, serum vitamin B12, serum phosphorus, and serum calcium were adjusted

Discussion

In this population-based study of US adults, we found that (1) the association between plasma tHcy level and lumbar BMD was different in various age groups; (2) this association was also different by gender. To the best of our knowledge, this study is thus far the largest sample size study on the relationship between plasma tHcy level and BMD in a general adult population.

It was reported that Hcy could affect proper bone metabolism. In vitro, high Hcy levels can modulate the bone remodeling process by inhibiting osteoblastic differentiation; inducing apoptosis in osteoblasts, osteocytes, and bone marrow stromal cells; and increasing osteoclast activity [23]. In an experimental rat model with hyperhomocysteinemia induced by a 2.4% methionine-enriched diet for 12 weeks, there was a significant decrease in the bone formation marker osteocalcin and an increase in urinary N-terminal collagen I telopeptides compared with normal rats [24]. In the same rat model with hyperhomocysteinemia, a 19% reduction occurred in the bone volume at the femoral neck and a 45% reduction at the distal femur [25].

The relationship of the Hcy level with BMD has been studied in various populations, but the conclusions remain inconsistent. Some studies demonstrated that high Hcy levels were associated with lower BMD [4-11], whereas other studies found no significant correlation between them [15-21]. In particular, unlike all the mentioned research, some studies reported mixed correlations. The results of a cross-sectional survey of 446 post-menopausal women showed that tHcy levels were negatively associated with BMD of the total femur, but not of the femoral neck or lumbar spine [13]. The results of a cross-sectional study of 3337 healthy Korean adults suggested that the correlation of Hcy with BMD was different based on sex and age, and BMD of the lumbar spine and femur decreased as the tHcy levels increased in women aged <50 years, but no significant correlation was found in all age groups of men or other age groups of women [14]. The results of the Hordaland Homocysteine Study (3070 women and 2268 men, aged 47–50 and 71–75 years) suggest that plasma tHcy level is an independent risk factor for low BMD in women but not in men [12]. In our study, we conducted subgroup analyses to make a better use of the data following the STrengthening the Reporting of OBservational studies in Epidemiology (STROBE) guideline [26]. Further, we found that the association between plasma tHcy level and lumbar BMD differed by age and gender. Therefore, these conflicting conclusions may be attributed to the heterogeneity among studies, including study design, study size, and differences in participants’ selection, such as age and gender.

The size is the major strength of this study; we investigated plasma tHcy and lumbar BMD of 10748 samples of the multiracial population. Thus, subgroup analyses could be performed due to the large sample size. In addition, we used smooth curve fittings to explore potential non-linear relationships. However, some limitations should be noted. The cross-sectional nature of the study precludes any inferences about causality. Besides, because the participants with cancer were excluded, study results cannot be generalized to these special populations. Furthermore, the data of this study were analyzed with no exclusion of other diseases that may influence bone health. This makes the results more generalizable but may weaken the observed association.

In conclusion, we found that the correlation between plasma tHcy level within normal range and lumbar BMD differed according to age and gender. Our findings provide new insights to advance the research of the link between Hcy and bone health.