Association Between Dietary Selenium Intake and Bone Mineral Density: A Cross-sectional Study Using National Health and Nutrition Examination Survey

DOI: https://doi.org/10.21203/rs.3.rs-1506021/v1

Abstract

Assessments of whether dietary selenium intake is related to bone health are scarce, with few relevant studies limited by its small sample. The aim of present study was to investigated the association between dietary selenium intake and bone mineral density (BMD) levels in different sites, including total femur, femur neck, trochanter and intertrochanter in femur and lumbar spine in the National Health and Nutrition Examination Survey (NHANES) database. Generalized linear models for the association between dietary selenium intake and BMD and generalized additive model for the dose-response relationship were used. A total of 21939 participants were included, and the mean age was 40.68 ± 22.42 years, and 51.28% were male. In multivariable adjustment model, participants had highest quintiles of dietary selenium intake (Q5) were associated with increased BMD levels in total femur (β=0.014, 95CI%: 0.008, 0.020, P<0.001), femur neck (β=0.010, 95CI%: 0.004, 0.016, P=0.001), trochanter (β=0.011, 95CI%: 0.005, 0.017, P<0.001), intertrochanter (β=0.017, 95CI%: 0.010, 0.025, P<0.001) and lumbar spine (β=0.013, 95CI%: 0.005, 0.020, P<0.001) compared with those in quintiles 1 (Q1). The dose-response relationship showed the inverted U-shape relationship between dietary selenium relationship and BMD levels (P for non-linearity <0.05). Participants tended to have increased levels of BMD as the dietary selenium intake increased when dietary selenium was below the turning point, and then a negative relation was observed when dietary was higher than the turning point. Our study indicated that higher dietary selenium intake was associated with increased BMD levels in total femur, femur neck, trochanter, intertrochanter, and lumbar spine, and these relationships were nonlinear. Future high-quality, prospective longitudinal studies are needed to confirm these findings.

Introduction

Osteoporosis is a systemic bone disease that is common worldwide. Approximately 10.2 million people aged 50 years and over had osteoporosis in 2010 in the US, and the number is expected to reach 13.5 million by 2030 [1, 2]. Osteoporosis, characterized by reduced bone mass, low bone mineral density (BMD), and bone microstructure deterioration [3, 4], increases the risk of all-cause mortality, including cardiovascular- and cancer-related mortality [57]. The etiology of reduced bone mass and the development of osteoporosis is related to multiple factors, including genetic, environmental, and dietary factors [8, 9]. Furthermore, oxidative stress has been implicated as a causative factor for many disease states, including the diminished BMD in osteoporosis. Therefore, including antioxidant-rich whole plant foods in the diet and adopting lifestyle modifications to maintain appropriate antioxidant levels may increase BMD and reduce the risk of brittleness-related fractures [10].

The element selenium is an essential micronutrient that forms the so-called “selenoproteins” when incorporated into the polypeptide chain of proteins. Selenoproteins and Se-dependent enzymes are involved in many crucial biological mechanisms, such as antioxidant and anti-inflammatory pathways [11, 12], intracellular redox regulation [13], and thyroid hormone metabolism [14]. A previous observational study demonstrated that the serum Se concentration was positively associated with bone outcomes, including BMD and fracture risk [15]. In a cross-sectional study, the relationship between serum selenium and the risk of osteoporosis-related fracture was nonlinear, but a strong positive correlation was evident between osteoporosis-related fracture risk and relatively high selenium exposure [16]. At present, few studies have examined the effects of dietary selenium intake on bone health, and those that have are limited by small sample sizes.

One study suggested that antioxidant intake, including Se, β-carotene, vitamin C, and vitamin E, was inversely associated with the risk of osteoporotic hip fracture in an older population of smokers [17]. However, a prospective population-based cohort study indicated that sodium selenite supplementation did not affect bone turnover markers or physical performance in postmenopausal women with osteopenia [18].

This study investigated the association between dietary selenium intake and BMD in various bones, including the femur, femur neck, trochanter and intertrochanter in the femur, and lumbar spine. Data on dietary intake was drawn from the cross-sectional National Health and Nutrition Examination Survey (NHANES). We examined the nonlinear dose–response association between dietary selenium intake and BMD.

Methods

The protocol of the present study was approved by the Institutional Review Board of the National Center for Health Statistics in United States. All participants gave their written informed consent. All authors declared that all methods in this study were carried out in accordance with relevant guidelines and regulations.

Study population

We used the data from NHANES, which is a population-based cross-sectional survey designed to assess the health, nutritional status, and potential risk factors of the civilian, noninstitutionalized population of the US. The consecutive surveys are conducted by the National Center for Health Statistics (NCHS) of the Centers for Disease Control (CDC) via in-person interviews, physical examinations and laboratory tests in a mobile examination center (MEC). Approximately, 5000 individuals at 15 geographic sites are selected by a multistage, stratified probability sampling design every 2 years. NHANES was approved by the National Center for Health Statistics Institutional Review Board, and each participants provided written informed consent. We extracted and aggregated data from four cycles of NHANES (2005–2010, 2013-2014). Participants with missing information on dietary selenium intake and BMD were excluded. Lastly, data on 21 939 participants were available for analysis of total femur, femur neck, trochanter and intertrochanter in femur, and 18 116 participants were available for analysis of lumbar spine. The detailed description of the NHANES was published elsewhere[6,19].

Assessment of dietary selenium intake

Dietary selenium intake and other food components, including dietary fiber and calcium intake were obtained from 24-hour dietary recall interviews which was conducted in the Mobile Examination Center (MEC). During the interview, participants were asked to recall the details of food and beverages consumed in 24-hours period before the interview. For each participant, nutrient intake from each food or beverage were estimated. The dietary selenium intake was calculated as microgram per day (μg/day).

Bone mineral density measurement 

BMD measurement was obtained using the method of dual-energy x-ray absorptiometry (DXA) with Hologic QDR 4500A fan-beam densitometers (Hologic Inc., Bedford, MA, USA)[6], which is an internationally accepted standard-of-care screening tool used to assess fragility-fracture risk[20]. The DXA examinations were administered by trained and certified radiology technologists. DXA scans were administered to eligible survey participants 8 years of age and older. Pregnant female, self-reported history of radiographic contrast material in the past 7 days, and measured weighted over 300 pounds were ineligible for the DXA examination.

Other covariates

The covariates involved age (years), sex (male and female), race (non-Hispanic white, non-Hispanic black, Mexican American, or other race/ethnicity), education (under high school, high school, or above high school), and family income (Under $20,000, $20,000-55000, or $20,000 and over) were obtained from in-person household interviews. Body mass index (BMI) was calculated as weight divided by height squared (kg/m2), and was classified into underweight (<18.5 kg/m2), normal weight (18.5-24.9 kg/m2), overweight (25-29.9 kg/m2), or obese (≥30 kg/m2). Furthermore, leisure time physical activity (<500 MET/week, 500-999 MET/week, or ≥1000 MET/week), smoke status (yes or no), alcohol use (yes or no), diabetes (yes or no), and hypertension (yes or no) were ascertained using questionnaires, which were self-reported by the participants.

Statistical analysis

Descriptive characteristics were summarized using means ± standard deviation (SD) for continuous variables and proportions for categorical variables. Data analyses accounted for the masked variance and used the recommended weighting methodology. Dietary selenium intake was divided into quintiles, and the differences between quintiles of selenium in participants characteristics were determined using one-way ANOVA or χ2 tests. We used generalized linear models to assess the association between dietary selenium intake and BMD. Three models were used in current analysis: model 1 was crude model without adjustment for potential confounders; model 2 was adjusted for age and sex; model 3 was further adjusted for race, education, family income, body mass index, leisure time physical activity, smoke status, alcohol use, dietary fiber intake, calcium intake, diabetes, and hypertension. We also assess the dose-response relationship between dietary selenium intake and BMD by using a generalized additive model (GAM). Moreover, to test the robustness of our findings, we used subgroup analysis by age (<50 vs. ≥50 years) and sex. All statistical analysis were performed using R software (version 3.5.3) and EmpowerStats (R) (www.empowerstats.com, X&Y solutions, Inc. Boston MA). Statistical significance was indicated by a two-sided P value <0.05.

Results

Characteristics of the participants included in present study are shown in Table 1. Of the 21939 participants, the mean age was 40.68 ± 22.42 years, and 51.28% were male. Approximately, 43% of participants were non-Hispanic white, and 44.54% of participants had education level less than high school. Compare with participants in quintiles 1 (Q1), those in higher quintiles (Q2-Q5) tended to be younger, male, non-Hispanic white, and have higher educated, higher family income, more smoke and alcohol use, more leisure time physical activity, dietary fiber and calcium intake, and lower prevalence of diabetes and hypertension.

Table 1

Characteristics of the participants

 

Total population (n = 21939)

Dietary selenium intake (Quintile)

Q1 (n = 4376)

Q2 (n = 4395)

Q3 (n = 4369)

Q4 (n = 4403)

Q5 (n = 4396)

P value

Age, years

40.68 ± 22.42

41.87 ± 24.24

41.35 ± 23.66

40.29 ± 22.86

40.67 ± 21.67

39.21 ± 19.28

< 0.001

Male, %

11251 (51.28)

1509 (34.48)

1802 (41.00)

2141 (49.00)

2523 (57.30)

3276 (74.52)

< 0.001

Race, n (%)

           

< 0.001

Non-Hispanic white

9510 (43.35)

1782 (40.72)

1916 (43.59)

1864 (42.66)

2005 (45.54)

1943 (44.20)

 

Black

4684 (21.35)

1088 (24.86)

945 (21.50)

916 (20.97)

857 (19.46)

878 (19.97)

 

Mexican American

4549 (20.73)

877 (20.04)

926 (21.07)

943 (21.58)

885 (20.10)

918 (20.88)

 

Other Hispanic

1920 (8.75)

398 (9.10)

371 (8.44)

391 (8.95)

382 (8.68)

378 (8.60)

 

Other race/ethnicity

1276 (5.82)

231 (5.28)

237 (5.39)

255 (5.84)

274 (6.22)

279 (6.35)

 

Education, n (%)

           

< 0.001

Under high school

9764 (44.54)

2212 (50.62)

2049 (46.63)

1978 (45.32)

1812 (41.19)

1713 (38.99)

 

High school

4055 (18.50)

799 (18.28)

799 (18.18)

782 (17.92)

816 (18.55)

859 (19.55)

 

Above high school

8102 (36.96)

1359 (31.10)

1546 (35.18)

1605 (36.77)

1771 (40.26)

1821 (41.45)

 

Family income, %

           

< 0.001

Under $20,000

4830 (22.83)

1122 (26.71)

1021 (24.00)

956 (22.71)

845 (19.87)

886 (20.89)

 

$20,000–55,000

8239 (38.94)

1706 (40.61)

1676 (39.40)

1664 (39.52)

1631 (38.35)

1562 (36.83)

 

$55,000 and over

8090 (38.23)

1373 (32.68)

1557 (36.60)

1590 (37.77)

1777 (41.78)

1793 (42.28)

 

Height, cm

164.73 ± 12.72

161.32 ± 12.01

162.00 ± 12.60

163.77 ± 12.73

166.19 ± 12.44

170.35 ± 11.65

< 0.001

Weight, cm

73.41 ± 21.62

70.08 ± 20.60

70.51 ± 21.40

72.20 ± 21.67

75.17 ± 21.64

79.04 ± 21.47

< 0.001

Body mass index, kg/m2

26.66 ± 6.26

26.56 ± 6.31

26.43 ± 6.40

26.50 ± 6.30

26.83 ± 6.15

26.96 ± 6.13

0.013

Smoke status, %

7532 (47.54)

1406 (46.40)

1433 (46.26)

1399 (45.47)

1605 (49.08)

1689 (50.15)

< 0.001

Alcohol use, %

4503 (20.53)

612 (13.99)

742 (16.88)

858 (19.64)

972 (22.08)

1319 (30.00)

< 0.001

Leisure time physical activity, MET/week

           

0.045

< 500

10029 (54.68)

2118 (59.16)

2032 (57.29)

1996 (55.49)

1980 (53.35)

1903 (48.73)

 

500–999

2251 (12.27)

395 (11.03)

463 (13.05)

439 (12.20)

474 (12.77)

480 (12.29)

 

≥ 1000

6060 (33.04)

1067 (29.80)

1052 (29.66)

1162 (32.30)

1257 (33.87)

1522 (38.98)

 

Dietary fiber intake, g

15.67 ± 9.68

10.56 ± 7.23

13.43 ± 7.83

15.31 ± 8.39

17.18 ± 8.85

21.85 ± 11.60

< 0.001

Calcium intake, mg

936.97 ± 601.26

553.73 ± 363.11

768.28 ± 397.13

895.80 ± 446.07

1061.54 ± 529.48

1403.27 ± 792.23

< 0.001

Diabetes, %

1852 (8.59)

417 (9.70)

374 (8.67)

375 (8.75)

373 (8.63)

313 (7.23)

0.002

Hypertension, %

5733 (32.21)

1264 (36.96)

1193 (34.30)

1118 (32.51)

1134 (31.31)

1024 (26.65)

< 0.001

Bone mass density

             

Total femur, gm/cm2

0.95 ± 0.17

0.91 ± 0.17

0.93 ± 0.17

0.94 ± 0.17

0.97 ± 0.17

1.00 ± 0.17

< 0.001

Femur neck, gm/cm2

0.83 ± 0.16

0.81 ± 0.16

0.81 ± 0.16

0.83 ± 0.16

0.84 ± 0.15

0.87 ± 0.16

< 0.001

Trochanter, gm/cm2

0.72 ± 0.14

0.70 ± 0.14

0.70 ± 0.14

0.72 ± 0.14

0.74 ± 0.14

0.76 ± 0.14

< 0.001

Intertrochanter, gm/cm2

1.11 ± 0.21

1.07 ± 0.20

1.08 ± 0.21

1.11 ± 0.20

1.13 ± 0.20

1.18 ± 0.20

< 0.001

Lumbar spine, gm/cm2

0.97 ± 0.19

0.95 ± 0.19

0.95 ± 0.20

0.96 ± 0.19

0.98 ± 0.19

1.01 ± 0.17

< 0.001

 

The unadjusted and multivariable adjusted associations between dietary selenium intake and BMD, including total femur, femur neck, trochanter, intertrochanter and lumbar spine are shown in Table 2. In crude model, compared with the lowest quintiles of dietary selenium intake, higher dietary selenium intake was linked with increased BMD levels (all P-trend < 0.001). After adjustment for potential confounders, the results were also consistent. In multivariable model, participants had highest quintiles of dietary selenium intake (Q5) were associated with increased BMD levels in total femur (β = 0.014, 95CI%: 0.008, 0.020, P<0.001), femur neck (β = 0.010, 95CI%: 0.004, 0.016, P = 0.001), trochanter (β = 0.011, 95CI%: 0.005, 0.017, P<0.001), intertrochanter (β = 0.017, 95CI%: 0.010, 0.025, P<0.001) and lumbar spine (β = 0.013, 95CI%: 0.005, 0.020, P<0.001) compared with those in quintiles 1 (Q1).

Table 2

The association between dietary selenium intake and bone mass density

Dietary selenium intake (Quintile)

Model 1

Model 2

Model 3

β (95% CI)

P value

β (95% CI)

P value

β (95% CI)

P value

Total femur, gm/cm2

           

Q1

Reference

 

Reference

 

Reference

 

Q2

0.011 (0.004, 0.018)

0.002

0.005 (-0.002, 0.012)

0.150

0.005 (-0.001, 0.010)

0.091

Q3

0.028 (0.021, 0.035)

< 0.001

0.015 (0.008, 0.021)

< 0.001

0.007 (0.002, 0.013)

0.011

Q4

0.051 (0.044, 0.058)

< 0.001

0.030 (0.023, 0.037)

< 0.001

0.012 (0.006, 0.018)

< 0.001

Q5

0.090 (0.083, 0.097)

< 0.001

0.052 (0.045, 0.060)

< 0.001

0.014 (0.008, 0.020)

< 0.001

P-trend

< 0.001

 

< 0.001

 

< 0.001

 

Femur neck, gm/cm2

           

Q1

Reference

 

Reference

 

Reference

 

Q2

0.005 (-0.001, 0.012)

0.117

0.001 (-0.006, 0.007)

0.838

0.002 (-0.003, 0.007)

0.505

Q3

0.020 (0.013, 0.027)

< 0.001

0.009 (0.003, 0.015)

0.006

0.005 (-0.001, 0.010)

0.078

Q4

0.035 (0.029, 0.042)

< 0.001

0.020 (0.014, 0.027)

< 0.001

0.008 (0.003, 0.014)

0.004

Q5

0.066 (0.060, 0.073)

< 0.001

0.039 (0.032, 0.045)

< 0.001

0.010 (0.004, 0.016)

0.001

P-trend

< 0.001

 

< 0.001

 

< 0.001

 

Trochanter, gm/cm2

           

Q1

Reference

 

Reference

 

Reference

 

Q2

0.009 (0.003, 0.015)

0.003

0.003 (-0.002, 0.009)

0.228

0.003 (-0.002, 0.008)

0.189

Q3

0.021 (0.015, 0.026)

< 0.001

0.009 (0.003, 0.014)

0.002

0.003 (-0.001, 0.008)

0.172

Q4

0.040 (0.034, 0.045)

< 0.001

0.021 (0.016, 0.027)

< 0.001

0.009 (0.004, 0.014)

< 0.001

Q5

0.069 (0.063, 0.075)

< 0.001

0.037 (0.031, 0.042)

< 0.001

0.011 (0.005, 0.017)

< 0.001

P-trend

< 0.001

 

< 0.001

 

< 0.001

 

Intertrochanter, gm/cm2

           

Q1

Reference

 

Reference

 

Reference

 

Q2

0.013 (0.005, 0.022)

0.002

0.007 (-0.001, 0.015)

0.099

0.006 (0.000, 0.013)

0.048

Q3

0.034 (0.025, 0.042)

< 0.001

0.020 (0.011, 0.028)

< 0.001

0.011 (0.004, 0.017)

0.001

Q4

0.060 (0.052, 0.069)

< 0.001

0.037 (0.029, 0.046)

< 0.001

0.015 (0.008, 0.022)

< 0.001

Q5

0.106 (0.097, 0.114)

< 0.001

0.066 (0.057, 0.074)

< 0.001

0.017 (0.010, 0.025)

< 0.001

P-trend

< 0.001

 

< 0.001

 

< 0.001

 

Lumbar spine, gm/cm2

           

Q1

Reference

 

Reference

 

Reference

 

Q2

-0.001 (-0.009, 0.008)

0.892

-0.001 (-0.009, 0.008)

0.866

-0.001 (-0.008, 0.005)

0.654

Q3

0.008 (-0.000, 0.017)

0.061

0.010 (0.002, 0.019)

0.017

0.003 (-0.003, 0.010)

0.329

Q4

0.032 (0.023, 0.041)

< 0.001

0.031 (0.023, 0.039)

< 0.001

0.011 (0.004, 0.018)

0.001

Q5

0.057 (0.048, 0.066)

< 0.001

0.057 (0.048, 0.065)

< 0.001

0.013 (0.005, 0.020)

< 0.001

P-trend

< 0.001

 

< 0.001

 

< 0.001

 
Model 1: crude model; Model 2: Adjusted for age and sex; Model 3: Adjusted for age, sex (male and female), race (non-Hispanic white, non-Hispanic black, Mexican American, other race/ethnicity or missing), education (under high school, high school, above high school, or missing), family income (Under $20,000, $20,000-55000, $20,000 and over, or missing), body mass index (underweight, normal weight, overweight, obese, or missing), leisure time physical activity (<500 MET/week, 500–999 MET/week, ≥ 1000 MET/week, or missing), smoke status (yes, no, or missing), alcohol use (yes, no, or missing), dietary fiber intake, calcium intake, diabetes (yes, no, or missing), and hypertension (yes, no, or missing).


We further used GAM to examine the dose-response association, and found the inverted U-shape association between dietary selenium intake and BMD levels (Fig. 1). Overall, participants tended to have increased levels of BMD as the dietary selenium intake increased when dietary selenium intake was below the turning point, and then the BMD decrease as the dietary selenium intake increased when dietary selenium intake was higher than the turning point. Furthermore, we used a log likelihood ratio test showed the significant differences between linear model to and segmented regression model, indicating that the associations between dietary selenium intake and BMD levels of total femur (P < 0.001), femur neck (P = 0.003), trochanter (P = 0.005), intertrochanter (P < 0.001), and lumbar spine (P = 0.004) were non-linear.

In the subgroups stratified by age, the highest quintiles of dietary selenium intake were positively associated with increased BMD levels in total femur, femur neck, trochanter, intertrochanter, and lumbar spine (Table 3). By contrast, no statistical associations were observed between highest quintiles of dietary selenium intake and BMD level in female.

Table 3

The association between dietary selenium intake and bone mass density by different age and sex

 

Age<50

Age ≥ 50

Male

Female

β (95% CI)

P value

β (95% CI)

P value

β (95% CI)

P value

β (95% CI)

P value

Total femur, gm/cm2

               

Q1

Reference

 

Reference

 

Reference

 

Reference

 

Q2

-0.002 (-0.008, 0.005)

0.641

0.013 (0.004, 0.022) 2

0.003

0.006 (-0.003, 0.015)

0.186

0.002 (-0.005, 0.008)

0.620

Q3

0.004 (-0.003, 0.011)

0.221

0.013 (0.004, 0.022)

0.005

0.005 (-0.004, 0.014)

0.260

0.006 (-0.001, 0.013)

0.082

Q4

0.010 (0.003, 0.018)

0.005

0.017 (0.008, 0.027)

< 0.001

0.010 (0.001, 0.019)

0.029

0.012 (0.005, 0.020)

0.001

Q5

0.015 (0.007, 0.022)

< 0.001

0.020 (0.009, 0.030)

< 0.001

0.016 (0.007, 0.025)

< 0.001

0.004 (-0.005, 0.013)

0.356

P-trend

< 0.001

 

< 0.001

 

< 0.001

 

0.012

 

Femur neck, gm/cm2

               

Q1

Reference

 

Reference

 

Reference

 

Reference

 

Q2

0.000 (-0.007, 0.007)

0.959

0.006 (-0.002, 0.014)

0.156

0.002 (-0.006, 0.011)

0.599

0.000 (-0.006, 0.006)

0.960

Q3

0.005 (-0.002, 0.012)

0.158

0.009 (0.001, 0.018)

0.031

0.002 (-0.006, 0.010)

0.646

0.006 (-0.001, 0.012)

0.095

Q4

0.010 (0.002, 0.017)

0.009

0.013 (0.005, 0.022)

0.003

0.005 (-0.003, 0.014)

0.205

0.010 (0.003, 0.017)

0.005

Q5

0.012 (0.004, 0.019)

0.003

0.019 (0.009, 0.029)

< 0.001

0.012 (0.004, 0.021)

0.005

0.001 (-0.008, 0.009)

0.869

P-trend

< 0.001

 

< 0.001

 

0.002

 

0.055

 

Trochanter, gm/cm2

               

Q1

Reference

 

Reference

 

Reference

 

Reference

 

Q2

-0.002 (-0.008, 0.004)

0.484

0.010 (0.002, 0.017)

0.014

0.005 (-0.003, 0.013)

0.221

-0.000 (-0.006, 0.006)

0.976

Q3

0.001 (-0.006, 0.007)

0.848

0.008 (-0.000, 0.016)

0.059

0.002 (-0.006, 0.010)

0.610

0.002 (-0.004, 0.008)

0.559

Q4

0.007 (0.001, 0.013)

0.032

0.012 (0.004, 0.020)

0.004

0.007 (-0.001, 0.015)

0.078

0.008 (0.002, 0.015)

0.013

Q5

0.011 (0.004, 0.017)

0.003

0.015 (0.005, 0.024)

0.002

0.013 (0.004, 0.021)

0.003

0.004 (-0.004, 0.012)

0.344

P-trend

< 0.001

 

0.003

 

0.002

 

0.034

 

Intertrochanter, gm/cm2

               

Q1

Reference

 

Reference

 

Reference

 

Reference

 

Q2

-0.002 (-0.010, 0.006)

0.678

0.017 (0.007, 0.028)

0.001

0.007 (-0.003, 0.018)

0.163

0.004 (-0.004, 0.011)

0.372

Q3

0.007 (-0.001, 0.015)

0.092

0.017 (0.007, 0.028)

0.002

0.008 (-0.003, 0.018)

0.144

0.010 (0.001, 0.018)

0.024

Q4

0.012 (0.004, 0.021)

0.004

0.021 (0.010, 0.033)

< 0.001

0.012 (0.001, 0.022)

0.030

0.016 (0.007, 0.025)

< 0.001

Q5

0.018 (0.009, 0.027)

< 0.001

0.022 (0.010, 0.035)

< 0.001

0.018 (0.007, 0.029)

< 0.001

0.006 (-0.005, 0.017)

0.251

P-trend

< 0.001

 

< 0.001

 

< 0.001

 

0.006

 

Lumbar spine, gm/cm2

               

Q1

Reference

 

Reference

 

Reference

 

Reference

 

Q2

-0.006 (-0.014, 0.001)

0.072

0.005 (-0.008, 0.017)

0.469

-0.001 (-0.011, 0.010)

0.887

-0.004 (-0.012, 0.003)

0.270

Q3

0.002 (-0.006, 0.009)

0.648

0.004 (-0.009, 0.017)

0.555

-0.000 (-0.010, 0.010)

0.960

0.003 (-0.006, 0.011)

0.502

Q4

0.013 (0.005, 0.020)

0.001

0.006 (-0.007, 0.020)

0.370

0.006 (-0.004, 0.016)

0.232

0.014 (0.005, 0.023)

0.002

Q5

0.017 (0.009, 0.025)

< 0.001

0.003 (-0.012, 0.018)

0.705

0.016 (0.005, 0.026)

0.003

-0.002 (-0.013, 0.009)

0.759

P-trend

< 0.001

 

0.642

 

< 0.001

 

0.055

 
Model 1: crude model; Model 2: Adjusted for age and sex; Model 3: Adjusted for age, sex (male and female), race (non-Hispanic white, non-Hispanic black, Mexican American, other race/ethnicity or missing), education (under high school, high school, above high school, or missing), family income (Under $20,000, $20,000-55000, $20,000 and over, or missing), body mass index (underweight, normal weight, overweight, obese, or missing), leisure time physical activity (<500 MET/week, 500–999 MET/week, ≥ 1000 MET/week, or missing), smoke status (yes, no, or missing), alcohol use (yes, no, or missing), dietary fiber intake, calcium intake, diabetes (yes, no, or missing), and hypertension (yes, no, or missing).

Discussion

We analyzed data from a large-scale nationally representative population-based cross-sectional study (NHANES) and determined that higher dietary selenium intake was associated with increased BMD in the femur, femur neck, trochanter, intertrochanter, and lumbar spine, and the effect remained consistent across various age groups and genders. Furthermore, the results from the generalized additive model suggested a nonlinear inverted U-shape association between dietary selenium intake and BMD.

Severe Se deficiency is associated with Keshan disease, an endemic osteoarticular cardiomyopathy that is characterized by the selective necrosis of articular and growth plate chondrocytes [21]. Bone has the second-highest proportion of selenium (16%) in the body, only exceeded by skeletal muscles (27.5%) [22]. Few studies have examined the association between dietary selenium intake and bone health, leading to inconsistent epidemiological results. One study demonstrated that dietary selenium supplementation did not attenuate mammary tumorigenesis–mediated bone loss in a male mouse breast cancer model [23]. A recent randomized double-blinded controlled study by Walsh et al. reported that 200 µg/day selenite supplementation did not affect the musculoskeletal health of postmenopausal women [24]. However, in the present study, a higher dietary selenium intake did result in increased BMD. In line with our findings, a cross-sectional study that included 6267 participants demonstrated that compared with those in the lowest quartile of dietary selenium intake, those belonging to the fourth quartile exhibited a lower odds ratio for osteoporosis (OR: 0.47, 95% CI: 0.31–0.73) [25]. Zhang et al. observed that Se intake was negatively associated with the risk of osteoporotic hip fracture [17].

The biological mechanisms responsible for the effects of Se intake on BMD are uncertain. A previous study demonstrated that changes in the redox state can alter the bone remodeling process, which allows continuous bone regeneration through the coordination of the three major types of bone cells: osteoclasts, osteoblasts, and osteocytes [26]. Changes in reactive oxygen species (ROS) and/or antioxidant systems may be involved in the pathogenesis of bone loss. Osteoblast and osteocyte apoptosis induced by ROS leads to osteoclast formation and inhibits mineralization and osteogenesis. Excessive osteocyte apoptosis is associated with oxidative stress that leads to imbalanced osteoclast formation, which results in increased bone remodeling and bone loss [2628]. Moreover, Se plays a crucial role in antioxidant, immunological, and anti-inflammatory processes. The physiological function of the essential micronutrient Se is mainly mediated by selenoproteins [29], which have antioxidant activities and are known to maintain the redox cell balance, protect against oxidative stress caused by ROS, and regulate inflammation and osteocyte proliferation and differentiation [30]. Furthermore, it has been reported that interleukin-6 (IL-6) and other cytokines play a crucial role in the pathogenesis of osteoporosis [31]. Therefore, the anti-inflammatory effect of Se may be partly mediated by inhibiting the activity of IL-6 and cytokines [32, 25]. Another potential mechanism linking Se to bone health is the relationship of Se-dependent glutathione peroxidase to thyroid protection [33]. Therefore, Se deficiency may increase the level of thyroid hormone in the blood, leading to accelerated bone loss and osteoporosis [34].

Although limited data confirm the effects of dietary selenium supplementation on bone health, the accumulated evidence indicates a positive association between circulating Se concentrations and bone outcomes. A population-based cohort study conducted in five European cities demonstrated that higher Se levels were associated with increased hip BMD and decreased bone formation at the beginning of the research [15]. Other studies have indicated that Se deficiency can hinder bone growth and alter bone metabolism [35, 36]. In a survey conducted in the United States, increased serum Se concentrations were associated with increased femur BMD, decreased Fracture Risk Assessment Tool (FRAX) scores, and a reduced history of bone fractures [37]. A study, that used plasma Se and selenoprotein P as biomarkers demonstrated that an increase in Se content was associated with an increase in BMD in the lumbar spine and hip in European postmenopausal women [38]. In addition, low hair selenium levels have been reported to be associated with low lumbar and femoral BMD values in Korean adults [39]. The finding of a positive correlation between dietary selenium intake and blood selenium concentration [40, 41] suggests that Se supplementation may positively influence bone health in selenium -deficient patients.

Some observational studies have reported a U-shape relationship between serum Se and the risk of diabetes, coronary heart disease, anemia, and all-cause mortality [40, 4244]. Data from NHANES Ⅲ indicate an inverse association between serum Se and all-cause mortality at low selenium levels (< 130 ng/mL) and a modest increase in mortality at high Se levels (> 150 ng/mL) [45]. Similar to the findings of these studies, the results of our study indicate a positive relationship between dietary selenium intake and BMD when dietary selenium intake is below a certain threshold and a negative relationship when dietary selenium intake is higher than that threshold. This may be because selenium is an essential element with a narrow safety margin, and higher concentrations often lead to toxicity [42]. Furthermore, the regulation of selenium levels in the body mainly depends on its excretion rather than absorption. When dietary selenium intake is high enough to optimize the levels of selenium protein, any further intake is completely offset by excreta, allowing for only a slight increase in systemic selenium [46].

The present study has several strengths. To our knowledge, this is the first study to assess the nonlinear relationship between dietary selenium intake and BMD in the femur, femur neck, trochanter, intertrochanter, and lumbar spine. Our study was based on a large nationally representative survey, and BMD was measured in a reliable independent lab using established methods. Moreover, we adjusted numerous potentially confounding factors, including socioeconomic, lifestyle, and nutrient intake factors. Our study also has some limitations. First, as a cross-sectional study, inferring a causal association between dietary selenium intake and BMD is challenging. Second, we did not assess the association between dietary selenium intake and the risk of osteoporosis because of the limitations of the original data. Furthermore, bone remodeling is a continuous physiological process that entails bone resorption by osteoclasts and bone formation by osteoblasts. Future high-quality, prospective longitudinal studies are required to confirm the findings of this study.

In summary, our study suggests that higher dietary selenium intake is associated with increased BMD in the femur, femur neck, trochanter, intertrochanter, and lumbar spine. Furthermore, this study identified an inverted U-shape relationship between dietary selenium intake and BMD. Future high-quality, prospective longitudinal studies are required to confirm these findings.

Conclusion

In summary, our study suggests that higher dietary selenium intake was associated with increased BMD levels in total femur, femur neck, trochanter, intertrochanter, and lumbar spine. Furthermore, the inverted U-shape relationship between dietary selenium intake and BMD levels were observed. Future high-quality, prospective longitudinal studies are needed to confirm these causalities.

Declarations

Acknowledgements We thank the third National Health and Nutrition Examination Survey for the permission for use of the data of the NHANES 2005–2010, 2013-2014.

Authors' contributions GX analyzed and interpreted the data, and was a major contributor in writing the manuscript. RL was the main contributor in the revising the contents of the original report. All authors read and approved the final manuscript.

Funding The authors declare that no funds, grants, or other support were received during the preparation of this manuscript.

Availability of data and materials The datasets analyzed during the current study are available in the website of the NHANES: https://www.cdc.gov/nchs/index.htm.

Compliance with Ethical Standards

Ethics approval and consent to participate This study was approved by the Institutional Review Board of the National Center for Health Statistics. All participants gave their written informed consent. All authors declared that all methods in this study were carried out in accordance with relevant guidelines and regulations.

Consent for Publication The authors consent the publication.

Code Availability Not applicable.

Conflict of Interest The authors declare no competing interests.

Code Availability Not applicable.

References

  1. Guo J, Huang Y, Bian S, Zhao C, Jin Y, Yu D, Wu X, Zhang D, Cao W, Jing F, Chen G (2018) Associations of urinary polycyclic aromatic hydrocarbons with bone mass density and osteoporosis in U.S. adults, NHANES 2005–2010. Environ Pollut 240:209–218. doi:10.1016/j.envpol.2018.04.108
  2. Wright NC, Looker AC, Saag KG, Curtis JR, Delzell ES, Randall S, Dawson-Hughes B (2014) The recent prevalence of osteoporosis and low bone mass in the United States based on bone mineral density at the femoral neck or lumbar spine. J Bone Miner Res 29 (11):2520–2526. doi:10.1002/jbmr.2269
  3. Hendrickx G, Boudin E, Van Hul W (2015) A look behind the scenes: the risk and pathogenesis of primary osteoporosis. Nat Rev Rheumatol 11 (8):462–474. doi:10.1038/nrrheum.2015.48
  4. Kanis JA, McCloskey EV, Johansson H, Oden A, Melton LJ, 3rd, Khaltaev N (2008) A reference standard for the description of osteoporosis. Bone 42 (3):467–475. doi:10.1016/j.bone.2007.11.001
  5. Bliuc D, Nguyen ND, Milch VE, Nguyen TV, Eisman JA, Center JR (2009) Mortality risk associated with low-trauma osteoporotic fracture and subsequent fracture in men and women. JAMA 301 (5):513–521. doi:10.1001/jama.2009.50
  6. Cai S, Fan J, Zhu L, Ye J, Rao X, Fan C, Zhong Y, Li Y (2020) Bone mineral density and osteoporosis in relation to all-cause and cause-specific mortality in NHANES: A population-based cohort study. Bone 141:115597. doi:10.1016/j.bone.2020.115597
  7. Lorentzon M, Nilsson AG, Johansson H, Kanis JA, Mellstrom D, Sundh D (2019) Extensive undertreatment of osteoporosis in older Swedish women. Osteoporos Int 30 (6):1297–1305. doi:10.1007/s00198-019-04872-4
  8. Rizzoli R, Bianchi ML, Garabedian M, McKay HA, Moreno LA (2010) Maximizing bone mineral mass gain during growth for the prevention of fractures in the adolescents and the elderly. Bone 46 (2):294–305. doi:10.1016/j.bone.2009.10.005
  9. Rizzoli R, Biver E, Brennan-Speranza TC (2021) Nutritional intake and bone health. Lancet Diabetes Endocrinol 9 (9):606–621. doi:10.1016/S2213-8587(21)00119-4
  10. Kimball JS, Johnson JP, Carlson DA (2021) Oxidative Stress and Osteoporosis. J Bone Joint Surg Am 103 (15):1451–1461. doi:10.2106/JBJS.20.00989
  11. Huang Z, Rose AH, Hoffmann PR (2012) The role of selenium in inflammation and immunity: from molecular mechanisms to therapeutic opportunities. Antioxid Redox Signal 16 (7):705–743. doi:10.1089/ars.2011.4145
  12. Mattmiller SA, Carlson BA, Sordillo LM (2013) Regulation of inflammation by selenium and selenoproteins: impact on eicosanoid biosynthesis. J Nutr Sci 2:e28. doi:10.1017/jns.2013.17
  13. Rayman MP (2012) Selenium and human health. Lancet 379 (9822):1256–1268. doi:10.1016/S0140-6736(11)61452-9
  14. Combs GF, Jr., Midthune DN, Patterson KY, Canfield WK, Hill AD, Levander OA, Taylor PR, Moler JE, Patterson BH (2009) Effects of selenomethionine supplementation on selenium status and thyroid hormone concentrations in healthy adults. Am J Clin Nutr 89 (6):1808–1814. doi:10.3945/ajcn.2008.27356
  15. Hoeg A, Gogakos A, Murphy E, Mueller S, Kohrle J, Reid DM, Gluer CC, Felsenberg D, Roux C, Eastell R, Schomburg L, Williams GR (2012) Bone turnover and bone mineral density are independently related to selenium status in healthy euthyroid postmenopausal women. J Clin Endocrinol Metab 97 (11):4061–4070. doi:10.1210/jc.2012-2121
  16. Galvez-Fernandez M, Grau-Perez M, Garcia-Barrera T, Ramirez-Acosta S, Gomez-Ariza JL, Perez-Gomez B, Galan-Labaca I, Navas-Acien A, Redon J, Briongos-Figuero LS, Duenas-Laita A, Perez-Castrillon JL, Tellez-Plaza M, Martin-Escudero JC (2021) Arsenic, cadmium, and selenium exposures and bone mineral density-related endpoints: The HORTEGA study. Free Radic Biol Med 162:392–400. doi:10.1016/j.freeradbiomed.2020.10.318
  17. Zhang J, Munger RG, West NA, Cutler DR, Wengreen HJ, Corcoran CD (2006) Antioxidant intake and risk of osteoporotic hip fracture in Utah: an effect modified by smoking status. Am J Epidemiol 163 (1):9–17. doi:10.1093/aje/kwj005
  18. Bonjour JP, Carrie AL, Ferrari S, Clavien H, Slosman D, Theintz G, Rizzoli R (1997) Calcium-enriched foods and bone mass growth in prepubertal girls: a randomized, double-blind, placebo-controlled trial. J Clin Invest 99 (6):1287–1294. doi:10.1172/JCI119287
  19. Johnson CL, Paulose-Ram R, Ogden CL, Carroll MD, Kruszon-Moran D, Dohrmann SM, Curtin LR (2013) National health and nutrition examination survey: analytic guidelines, 1999–2010. Vital Health Stat 2 (161):1–24
  20. CDC (2012) National Health and Nutrition Examination Survey, 2009–2010 Data documentation, Codebook, and Frequencies: dual Energy X-ray Absorptiometry.
  21. Hou J, Zhu L, Chen C, Feng H, Li D, Sun S, Xing Z, Wan X, Wang X, Li F, Guo X, Xiong P, Zhao S, Li S, Liu J, Sun D (2021) Association of selenium levels with the prevention and control of Keshan disease: A cross-sectional study. J Trace Elem Med Biol 68:126832. doi:10.1016/j.jtemb.2021.126832
  22. Zachara BA, Pawluk H, Bloch-Boguslawska E, Sliwka KM, Korenkiewicz J, Skok Z, Ryc K (2001) Tissue level, distribution, and total body selenium content in healthy and diseased humans in Poland. Arch Environ Health 56 (5):461–466. doi:10.1080/00039890109604483
  23. Yan L, Nielsen FH, Sundaram S, Cao J (2020) Dietary Selenium Supplementation Does Not Attenuate Mammary Tumorigenesis-Mediated Bone Loss in Male MMTV-PyMT Mice. Biol Trace Elem Res 194 (1):221–227. doi:10.1007/s12011-019-01767-7
  24. Walsh JS, Jacques RM, Schomburg L, Hill TR, Mathers JC, Williams GR, Eastell R (2021) Effect of selenium supplementation on musculoskeletal health in older women: a randomised, double-blind, placebo-controlled trial. Lancet Healthy Longev 2 (4):e212-e221. doi:10.1016/S2666-7568(21)00051-9
  25. Wang Y, Xie D, Li J, Long H, Wu J, Wu Z, He H, Wang H, Yang T, Wang Y (2019) Association between dietary selenium intake and the prevalence of osteoporosis: a cross-sectional study. BMC Musculoskelet Disord 20 (1):585. doi:10.1186/s12891-019-2958-5
  26. Domazetovic V, Marcucci G, Iantomasi T, Brandi ML, Vincenzini MT (2017) Oxidative stress in bone remodeling: role of antioxidants. Clin Cases Miner Bone Metab 14 (2):209–216. doi:10.11138/ccmbm/2017.14.1.209
  27. Henriksen K, Neutzsky-Wulff AV, Bonewald LF, Karsdal MA (2009) Local communication on and within bone controls bone remodeling. Bone 44 (6):1026–1033. doi:10.1016/j.bone.2009.03.671
  28. Robling AG, Bonewald LF (2020) The Osteocyte: New Insights. Annu Rev Physiol 82:485–506. doi:10.1146/annurev-physiol-021119-034332
  29. Rayman MP (2000) The importance of selenium to human health. Lancet 356 (9225):233–241. doi:10.1016/S0140-6736(00)02490-9
  30. Vescini F, Chiodini I, Palermo A, Cesareo R, De Geronimo V, Scillitani A, Gennari L, Falchetti A (2021) Selenium: A Trace Element for a Healthy Skeleton - A Narrative Review. Endocr Metab Immune Disord Drug Targets 21 (4):577–585. doi:10.2174/1871530320666200628030913
  31. Manolagas SC (1998) The role of IL-6 type cytokines and their receptors in bone. Ann N Y Acad Sci 840:194–204. doi:10.1111/j.1749-6632.1998.tb09563.x
  32. Duntas LH (2009) Selenium and inflammation: underlying anti-inflammatory mechanisms. Horm Metab Res 41 (6):443–447. doi:10.1055/s-0029-1220724
  33. Schomburg L, Kohrle J (2008) On the importance of selenium and iodine metabolism for thyroid hormone biosynthesis and human health. Mol Nutr Food Res 52 (11):1235–1246. doi:10.1002/mnfr.200700465
  34. Williams GR, Bassett JHD (2018) Thyroid diseases and bone health. J Endocrinol Invest 41 (1):99–109. doi:10.1007/s40618-017-0753-4
  35. Cao JJ, Gregoire BR, Zeng H (2012) Selenium deficiency decreases antioxidative capacity and is detrimental to bone microarchitecture in mice. J Nutr 142 (8):1526–1531. doi:10.3945/jn.111.157040
  36. Moreno-Reyes R, Egrise D, Neve J, Pasteels JL, Schoutens A (2001) Selenium deficiency-induced growth retardation is associated with an impaired bone metabolism and osteopenia. J Bone Miner Res 16 (8):1556–1563. doi:10.1359/jbmr.2001.16.8.1556
  37. Wu CC, Wang CK, Yang AM, Lu CS, Lin CY (2021) Selenium status is independently related to bone mineral density, FRAX score, and bone fracture history: NHANES, 2013 to 2014. Bone 143:115631. doi:10.1016/j.bone.2020.115631
  38. Al EAA, Parsian H, Fathi M, Faghihzadeh S, Hosseini SR, Nooreddini HG, Mosapour A (2018) ALOX12 gene polymorphisms and serum selenium status in elderly osteoporotic patients. Adv Clin Exp Med 27 (12):1717–1722. doi:10.17219/acem/75689
  39. Park KC, Kwon Y, Lee Y, Kim DK, Jang Y, Lee S (2020) Low selenium levels are associated with decreased bone mineral densities. J Trace Elem Med Biol 61:126534. doi:10.1016/j.jtemb.2020.126534
  40. Lin J, Shen T (2021) Association of dietary and serum selenium concentrations with glucose level and risk of diabetes mellitus: A cross sectional study of national health and nutrition examination survey, 1999–2006. J Trace Elem Med Biol 63:126660. doi:10.1016/j.jtemb.2020.126660
  41. Xia Y, Hill KE, Byrne DW, Xu J, Burk RF (2005) Effectiveness of selenium supplements in a low-selenium area of China. Am J Clin Nutr 81 (4):829–834. doi:10.1093/ajcn/81.4.829
  42. Xie B, Wang J, Zhang J, Chen M (2020) Dietary and serum selenium in coronary heart disease and all-cause mortality: An international perspective. Asia Pac J Clin Nutr 29 (4):827–838. doi:10.6133/apjcn.202012_29(4).0019
  43. Wang XL, Yang TB, Wei J, Lei GH, Zeng C (2016) Association between serum selenium level and type 2 diabetes mellitus: a non-linear dose-response meta-analysis of observational studies. Nutr J 15 (1):48. doi:10.1186/s12937-016-0169-6
  44. Zhou Q, Zhang B, Chen X, Chen Q, Hao L (2021) Association of serum selenium with anemia-related indicators and risk of anemia. Food Sci Nutr 9 (6):3039–3047. doi:10.1002/fsn3.2261
  45. Bleys J, Navas-Acien A, Guallar E (2008) Serum selenium levels and all-cause, cancer, and cardiovascular mortality among US adults. Arch Intern Med 168 (4):404–410. doi:10.1001/archinternmed.2007.74
  46. Burk RF, Hill KE (2015) Regulation of Selenium Metabolism and Transport. Annu Rev Nutr 35:109–134. doi:10.1146/annurev-nutr-071714-034250