- 1West China Hospital, Sichuan University, Chengdu, Sichuan, China
- 2West China School of Nursing, Sichuan University, Chengdu, Sichuan, China
Background: The association between dietary selenium intake and arthritis, rheumatoid arthritis (RA), and osteoarthritis (OA) is inconsistent in previous studies and remain unclear. To investigate their relationship, this study was performed.
Methods: Data from the National Health and Nutrition Examination Survey (2003–2016) were downloaded and further analyzed. Dietary Se intake was classified according to quartiles with quartile 1 (Q1) having the lowest intake and quartile 4 (Q4) having the highest intake. Weighted logistic regression was used to investigate the association between dietary selenium intake and arthritis, RA, and OA. Subgroup analyses were performed to verify the findings. To further examine the non-linear relationship between dietary selenium intake and OA, restricted cubic spline (RCS) was adopted.
Results: In the crude model, the highest level of dietary selenium intake was siginificantly associated with decreased risks of arthritis (OR: 0.40, 95% CI: 0.37, 0.44) and rheumatoid arthritis (OR: 0.47, 95% CI: 0.40, 0.54), respectively. In the fully adjusted model, dietary selenium intake was not associated with risk of arthritis and RA (all P > 0.05). Conversely, the risk of OA was noted for participants with higher selenium intake (odds ratio of quartile 4 = 1.33, 95% CI = 1.07–1.65, P < 0.05). In the subgroup analyses, participants with diabetes had a higher risk of OA when ingested high selenium levels than those without diabetes (P < 0.001). The results of RCS showed that significant overall trends were found between dietary selenium intake and osteoarthritis (P for overall < 0.05). However, non-linear association was not detected in this association (P for non-linear > 0.05).
Conclusion: Using data from NHANES, this study discloses that high dietary selenium intake might be associated with risk of OA. However, the generalization of conclusion needs further examination because of the limitation of dietary questionnaire survey.
1. Background
As the most common joint disease, osteoarthritis (OA) is prevalent globally, especially in the aged. As summarized by Sun et al. (1), the prevalence of OA is 25.03% for lumbar, 21.51% for knee, 20.46% for cervical vertebra in the middle-aged and elderly Chinese. This figure is higher in women than in men. OA is characterized by cartilage degradation, subchondral bone, and synovium lesions, leading to obvious joint pain and finally, loss of function (2). According to the report from Cross et al. (3), OA is the leading course of global disability and the years of life lived with disability increase from 10.5 million in 1990 to 17.1 million in 2010. The resultant socioeconomic burden and adverse impact for patients are heavy (4).
Joint replacement is an effective way to treat the end-stage OA with destroyed cartilage (2). However, the joint function after surgery can be poor and worse still, the lifespan of replaced joint is limited. Patients may face reoperation, consequently bringing more costs and pain. Therefore, identifying risk factors seems requisite to suppress the high incidence in the aged. Previous studies have disclosed that genetics, obesity, unhealthy diet can increase the risk of OA (5–7). Among these risk factors, the role of trace element selenium in the onset of OA is gradually noted. In China, Wang et al. (8). recruited 1,032 subjects aged ≥50 years and found that lower selenium is associated with higher risks of OA. Kurz et al. (9) also reported that diet selenium supplementation could prevent the development of mechanically induced OA, which echoed the findings from Wang et al. However, in Turkey, Yazar et al. (10) identified no association between synovial fluid and plasma selenium and OA. This insignificant finding is replicated in a dog model of OA (11). The inconsistent results may be attributed to sample size, ethnic disparity, and economic status. Therefore, more evidence with more is required to clarify the unclear association between selenium intake and OA.
To clarify the association between selenium intake and OA, we used the data from National Health and Nutrition Examination Survey (NHANES) in Americans, which is a well-designed cohort with adequate sample size (12). In addition, we also examined the non-linear relationship between selenium intake and OA, providing more evidence to determine the threshold effects of selenium.
2. Methods
2.1. Study population
A representative, non-institutionalized sample of the U.S. population has been collected on dietary habits and health condition every two years through the NHANES since 1999 as a national survey. It combines detailed in-person interviews, physical examinations, computer-based questionnaires, and laboratory tests to collect a wide range of quantitative and qualitative information (13). A detailed description of the NHANES survey methods can be found on the website (http://www.cdc.gov/nchs/nhanes/index.htm). A written informed consent was obtained from all participants before any data collection was conducted by the National Center for Health Statistics ethics review board.
Based on questionnaire information regarding arthritis and dietary nutrients intake, we used data from seven independent waves of NHANES (2003–2016). Participants were limited to those over 20 years of age, being non-pregnant, and having completed data of dietary selenium intake and arthritis outcomes.
2.2. Assessment of arthritis
The arthritis was assessed via the following question from NHANES codebook: “Has a doctor or other health professional ever told you that you had arthritis?”. Response options were “Yes” or “No.” Rheumatoid arthritis and osteoarthritis were assessed via the following question: “Which type of arthritis was it?”, and response options were “Rheumatoid arthritis,” “Osteoarthritis,” “Psoriatic arthritis,” “Other,” “Refused,” and “Don't know.” Percentage of psoriatic arthritis was quite low (64 cases, 0.91%), thus we did not consider the association between dietary selenium intake and psoriatic arthritis.
2.3. Assessment of dietary selenium intake
NHANES collected data on food intake over two non-consecutive days, with the first interview conducted in person and the second by phone. In the current analysis, we used the mean food intake over 2 days. Participants provided details of the meals they consumed in the past 24 h, and nutrients were estimated using the Food and Nutrient Database for Dietary Studies published by the USDA (14).
2.4. Covariates
In this study, the following variables were considered as covariates: Age (years), sex, education (above high school, high school, or less than high school), Races (Mexican American, non-Hispanic, non-Hispanic White, other Hispanic, and other races), Alcohol intake [never, former, current (Heavy, mild, and moderate)], smoking status (former, never, now), BMI, physical activity [metabolic equivalent (MET), min/week], poverty income ratio, diabetes, and cardiovascular disease (CVD). Poverty income ratio was calculated based on previous study, from <1.3 (Low), 1.3–3.5 (Median), and >3.5 (High) (15). Participants were considered to have diabetes or CVD: self-reported doctor's diagnosis of diabetes or CVD, currently taking medicine for controlling blood glucose or combating heart disease.
2.5. Statistical analysis
Dietary Se intakes were classified according to quartiles with quartile 1 (Q1) having the lowest intake and quartile 4 (Q4) having the highest intake. In this study, sample characteristics were presented as means and standard deviations or as percentages, and the means and standard deviations of subjects were compared using one-way ANOVA analysis or Chi-square tests. A complex multistage probability sampling process was accounted for by survey weights, sample strata, and sample clusters. A simple linear scaling of the 2-year weights created seven-wave weights (2003–2016). The association of arthritis, rheumatoid arthritis, and osteoarthritis with quartiles of Se intake was examined using logistic regression models adjusted for multivariables. The independent association between arthritis types and dietary selenium intake was examined using three models: model 1, no adjusted variable; model 2, adjusted for age, sex; model 3, adjusted for age, sex, races, education, alcohol intake, smoking status, poverty income ratio, MET, CVD, BMI and diabetes. Odds ratios (OR) and their associated 95% confidence intervals (CIs) were used to estimate the strength of the association for multivariate models.
Additional analyses were also conducted. First, selenium quartiles were transformed to continuous variables (values range from 1 to 4) and included into logistic regression model to find the potential linear trends in the association. Second, subgroup analyses were adopted stratified by age, sex, races, education, alcohol intake, smoke, poverty income ratio, CVD and diabetes. Non-linear association test was utilized to identify the modification of association between different sub population. Third, restricted cubic spline was plotted to find the dose-response trends between dietary selenium intake and arthritis risk, four knots were defined at 25, 50, 75, and 95% percentiles of dietary selenium intake, and reference value was located at the median value. The STATA version 17.0 were used to conduct all statistical analyses. A P-value <0.05 was considered statistically significant.
3. Results
As shown in Table 1, 26,620 participants with dietary selenium intake classified by quartile were included in this study. Participants in the quartile 4 group tend to be younger, male, non-Hispanic White, mild-drinker, non-smoker, obese, engaging in more physical activity, locating at highest poverty income ratio, receiving above high school education, and having no diabetes, CVD, arthritis, rheumatoid arthritis, and osteoarthritis (P < 0.05).
Compared to the reference level (Q1), the association between dietary selenium intake and arthritis was significantly negative in the null model [OR = 0.74 (95% CI: 0.69–0.79) for Q2, OR = 0.58 (95% CI: 0.54–0.62) for Q3, OR = 0.40 (95% CI: 0.37–0.44) Table 2], and significant linear trend was found at P for trend < 0.001. However, this association was insignificant after adjusting for covariates with sex and age (P > 0.05), and the OR of quartile 4 was 1.07 (95% CI: 0.95–1.20, Table 2). Meanwhile, this association was consistent when adjusted for other covariates, and the OR of quartile 4 was 1.06 (95% CI: 0.94–1.19, Table 2).
Similarly, the association between dietary selenium intake and rheumatoid arthritis was non-significant as comparison to quartile 1 in the full model. In the model 1, high selenium intake was associated with the decreased risk of rheumatoid arthritis [OR = 0.83 (95% CI: 0.73–0.95) for Q2, OR = 0.61 (95% CI: 0.53–0.70) for Q3, OR = 0.47 (95% CI: 0.40–0.54), Table 3], and the association was gradually increased (P for trend < 0.001). In the model 2 and model 3, high selenium intake was not associated with risk of rheumatoid arthritis, and the OR of quartile 4 was 0.95 (95% CI: 0.77–1.18) and 0.94 (95% CI: 0.76–1.18; Table 3).
The positive association between dietary selenium intake and osteoarthritis was found in the full model (Table 4). In the model 1, dietary selenium intake was related to the decreased risk of osteoarthritis [OR = 0.76 (95% CI: 0.66–0.87) for Q2, OR = 0.67 (95% CI: 0.58–0.77) for Q3, OR = 0.55 (95% CI: 0.48–0.64), Table 4]. However, the inverse association was found in model 2 and model 3. The OR of quartile 4 was 1.33 (95% CI: 1.07–1.65) for model 2, and 1.28 (95% CI: 1.03–1.59) for model 3, respectively (Table 4).
There was no significant difference between dietary selenium intake and osteoarthritis stratified with subgroup, except for diabetes (P for interaction > 0.05, Table 5). In detailed, the results indicated that this association was not supported in the subgroup analysis stratified by age, sex, races, education, alcohol intake, smoking, poverty income ratio, and CVD (Table 5). Besides, participants with diabetes had a greater risk of osteoarthritis when ingested high selenium levels than those without diabetes (P < 0.001), though this effect size was insignificant in each group (Table 5).
The restricted cubic spline showed that significant overall trends were found between dietary selenium intake and osteoarthritis (P for overall <0.05). However, non-linear association was not detected in this association (P for non-linear >0.05), which demonstrated that threshold value of dietary selenium intake for osteoarthritis needs further to research (Figure 1).
Figure 1. The RCS plot between dietary selenium intake and osteoarthritis. RCS, restricted cubic spline.
4. Discussion
In this study, we explored the association between selenium intake and OA in a large sample from NHANES. The risky role of high dietary selenium intake is noted for OA in Americans. Additionally, the non-linear relationship between selenium and OA is not detected, indicating the constantly rising risk with the increase of dietary selenium intake.
The association between selenium and OA is discordant in different studies. In a population-based cross-sectional study by Wang et al. (8), lower plasma selenium is found to be associated with higher risk of OA. This may be explained by the dysregulation of redox homeostasis triggered by selenium deficiency (16). The oxidative stress subsequently impairs cartilage and further leads to OA. However, as indicated by Perri et al. (17), low dietary selenium intake was common in the elderly. And the study by Wang et al. (8) was carried out among participants over 50 years old and the elderly are prone to OA, which may bias the findings. The evidence from Kurz et al. (9) in a rat model also supports this finding. They found that dietary selenium intake could prevent the progress of OA. The protective effect of high plasma selenium was also reported in a Mendelian randomization study by Qu et al. (18). But the protective effect was only observed in women, not in men, indicating the sex-specific association. In our study, high dietary selenium intake is found to increase the risk of OA, which is less reported in previous studies. Therefore, an excessive selenium intake is not recommended for the general population.
The mechanism linking selenium to OA still remains unknown. According to previous studies, metabolism associated inflammation, instead of oxidative stress, may be responsible for the adverse effect. One previous study based on mice model revealed that excessive intake of selenium could lead to the over-expression of GPx1, which further interfered with insulin signaling pathway (19). Transgenic mice overexpressing GPx1 developed hyperglycemia, and obesity, which may then trigger systemic inflammation and consequently OA (20, 21). Besides, some Se-rich food items, like red meat, processed meats, cheese and milk products has been suggested those appeared to increase risks of OA (22, 23). The high-fat and high-cholesterol food intake might raise body weight and promote OA risks.
Whether selenium supplementation is beneficial to the general population is controversial. Previous prospective studies have shown that high selenium is associated with low overall mortality (24, 25). The national survey in U.S. among 13,887 adults suggested that non-linear association was detected between increased serum selenium concentration and all-cause or cancer motility, but not in cardiovascular motility. Briefly, the inverse association was identified at selenium levels below 130 ng/ml, whereas moderate positive association was examined at selenium levels above 150 ng/ml (24). However, recent analysis from U.S. population suggested that higher serum selenium concentration could decrease the odds of all-cause and cardiovascular motility (26). In consistent results revealed that selenium intake requirement might be different across various characteristics, and corresponding threshold needs further investigation. The recommended average daily selenium intake is 60 μg for men and 53 μg for women (27). In our results (Figure 1), dietary selenium intake <100 μg don't increase the risk of OA. Hence, additional dietary selenium intake for patients with selenium deficiency are still safe and will not increase the risk of OA.
This study has some strengths and limitations. The main merit is the adequate sample size in Americans, allowing for higher statistical power. Second, sampling weights were considered in the analyses, reducing the bias from oversampling. The weighted dataset is more representative. Third, we also detected the non-linear relationship between dietary selenium intake and OA risk, providing more evidence for the threshold effect. This study also has some limitations. The principal demerit is the cross-sectional design of this study, which could not identify the causal association and avoid the bias from confounding factors. Future cohort studies should be performed to verify the findings. Additionally, the diagnosis of OA is mainly based on self-report, reducing the accuracy. Further verification using more accurate method like MRI or CT should be considered. Last but not least, we obtained dietary Se intake using questionnaire survey, which has some limitation in measurement accuracy, because of regional heterogeneity of the food production and bias from self-reported. Hence, further investigation should include biomarkers such as total serum selenium or selenoprotein P concentrations as supplement.
In conclusion, this study discloses that high dietary selenium intake is associated with risk of osteoarthritis. However, this risk is not detected for participants with dietary selenium intake <100 μg/day. For the general population, regular dietary selenium intake should be carefully considered.
5. Conclusions
Using data from NHANES, this study discloses that high dietary selenium intake might be associated with risk of osteoarthritis. However, current conclusion should be cautious because of limitation of questionnaire survey.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
XD and YT wrote, revised and reviewed the manuscript, drafted the study design, and supervised all processes. All authors contributed to the article and approved the submitted version.
Acknowledgments
Thanks to the NHANES team for collecting and sharing datasets.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Sun X, Zhen X, Hu X, Li Y, Gu S, Gu Y, et al. Osteoarthritis in the middle-aged and elderly in China: prevalence and influencing factors. Int J Environ Res Public Health. (2019) 16:4701. doi: 10.3390/ijerph16234701
2. Glyn-Jones S, Palmer AJ, Agricola R, Price AJ, Vincent TL, Weinans H, et al. Osteoarthritis. Lancet. (2015) 386:376–87. doi: 10.1016/S0140-6736(14)60802-3
3. Cross M, Smith E, Hoy D, Nolte S, Ackerman I, Fransen M, et al. The global burden of hip and knee osteoarthritis: estimates from the global burden of disease 2010 study. Ann Rheum Dis. (2014) 73:1323–30. doi: 10.1136/annrheumdis-2013-204763
4. Hiligsmann M, Cooper C, Arden N, Boers M, Branco JC, Luisa Brandi M, et al. Health economics in the field of osteoarthritis: an expert's consensus paper from the European Society for Clinical and Economic Aspects of Osteoporosis and Osteoarthritis (ESCEO). Semin Arthritis Rheum. (2013) 43:303–13. doi: 10.1016/j.semarthrit.2013.07.003
5. Ratneswaran A, Kapoor M. Osteoarthritis year in review: genetics, genomics, epigenetics Osteoarthritis Cartilage. (2021) 29:151–60. doi: 10.1016/j.joca.2020.11.003
6. Kulkarni K, Karssiens T, Kumar V, Pandit H. Obesity and osteoarthritis. Maturitas. (2016) 89:22–8. doi: 10.1016/j.maturitas.2016.04.006
7. Zhang Y, Ye M, Zhao Y, Xiong Y, Shen S, Yu Q, et al. Higher dietary Se intake is associated with the risk of new-onset fracture: a National Longitudinal Study for 20 years. Front Nutr. (2021) 8:719147. doi: 10.3389/fnut.2021.719147
8. Wang N, Xie M, Lei G, Zeng C, Yang T, Yang Z, et al. A cross-sectional study of association between plasma selenium levels and the prevalence of osteoarthritis: data from the xiangya osteoarthritis study. J Nutr Health Aging. (2022) 26:197–202. doi: 10.1007/s12603-022-1739-2
9. Kurz B, Jost B, Schünke M. Dietary vitamins and selenium diminish the development of mechanically induced osteoarthritis and increase the expression of antioxidative enzymes in the knee joint of STR/1N mice. Osteoarthritis Cartilage. (2002) 10119–26. doi: 10.1053/joca.2001.0489
10. Yazar M, Sarban S, Kocyigit A, Isikan UE. Synovial fluid and plasma selenium, copper, zinc, and iron concentrations in patients with rheumatoid arthritis and osteoarthritis. Biol Trace Elem Res. (2005) 106:123–32. doi: 10.1385/BTER:106:2:123
11. De Oliveira El-Warrak A, Rouma M, Amoroso A, Boysen SR, Chorfi Y. Measurement of vitamin A, vitamin E, selenium, and L-lactate in dogs with and without osteoarthritis secondary to ruptured cranial cruciate ligament. Can Vet J. (2012) 53:1285–8.
12. Zhang Y, He Q, Zhang W, Xiong Y, Shen S, Yang J, et al. Non-linear associations between visceral adiposity index and cardiovascular and cerebrovascular diseases: results from the NHANES (1999-2018). Front Cardiovasc Med. (2022) 9:908020. doi: 10.3389/fcvm.2022.908020
13. Johnson CL, Paulose-Ram R, Ogden CL, Carroll MD, Kruszon-Moran D, Dohrmann SM, et al. National health and nutrition examination survey: analytic guidelines, 1999-2010. Vital Health Stat 2. (2013) 161:1–24.
14. Ahluwalia N, Dwyer J, Terry A, Moshfegh A, Johnson C. Update on NHANES dietary data: focus on collection, release, analytical considerations, and uses to inform public policy. Adv Nutr. (2016) 7:121–34. doi: 10.3945/an.115.009258
15. Ogden CL, Carroll MD, Fakhouri TH, Hales CM, Fryar CD, Li X, et al. Prevalence of obesity among youths by household income and education level of head of household - United States 2011-2014. MMWR Morb Mortal Wkly Rep. (2018) 67:186–9. doi: 10.15585/mmwr.mm6706a3
16. Kang D, Lee J, Wu C, Guo X, Lee BJ, Chun JS, et al. The role of selenium metabolism and selenoproteins in cartilage homeostasis and arthropathies. Exp Mol Med. (2020) 52:1198–208. doi: 10.1038/s12276-020-0408-y
17. Perri G, Mendonça N, Jagger C, Walsh J, Eastell R, Mathers JC, et al. Dietary selenium intakes and musculoskeletal function in very old adults: analysis of the newcastle 85+ study. Nutrients. (2020) 12:2068. doi: 10.3390/nu12072068
18. Qu Z, Yang F, Hong J, Wang W, Li S, Jiang G, et al. Causal relationship of serum nutritional factors with osteoarthritis: a Mendelian randomization study. Rheumatology. (2021) 60:2383–90. doi: 10.1093/rheumatology/keaa622
19. Labunskyy VM, Lee BC, Handy DE, Loscalzo J, Hatfield DL, Gladyshev VN. Both maximal expression of selenoproteins and selenoprotein deficiency can promote development of type 2 diabetes-like phenotype in mice. Antioxid Redox Signal. (2011) 14:2327–36. doi: 10.1089/ars.2010.3526
20. Mcclung JP, Roneker CA, Mu W, Lisk DJ, Langlais P, Liu F, et al. Development of insulin resistance and obesity in mice overexpressing cellular glutathione peroxidase. Proc Natl Acad Sci U S A. (2004) 101:8852–7. doi: 10.1073/pnas.0308096101
21. Veronese N, Cooper C, Reginster JY, Hochberg M, Branco J, Bruyère O, et al. Type 2 diabetes mellitus and osteoarthritis. Semin Arthritis Rheum. (2019) 49:9–19. doi: 10.1016/j.semarthrit.2019.01.005
22. Jin J, Li J, Gan Y, Liu J, Zhao X, Chen J, et al. Red meat intake is associated with early onset of rheumatoid arthritis: a cross-sectional study. Sci Rep. (2021) 11:5681. doi: 10.1038/s41598-021-85035-6
23. Datta P, Zhang Y, Parousis A, Sharma A, Rossomacha E, Endisha H, et al. High-fat diet-induced acceleration of osteoarthritis is associated with a distinct and sustained plasma metabolite signature. Sci Rep. (2017) 7:8205. doi: 10.1038/s41598-017-07963-6
24. Bleys J, Navas-Acien A, Guallar E. Serum selenium levels and all-cause, cancer, and cardiovascular mortality among US adults Arch Intern Med. (2008) 168:404–10. doi: 10.1001/archinternmed.2007.74
25. Akbaraly NT, Arnaud J, Hininger-Favier I, Gourlet V, Roussel AM, Berr C. Selenium and mortality in the elderly: results from the EVA study. Clin Chem. (2005) 51:2117–23. doi: 10.1373/clinchem.2005.055301
26. Qiu Z, Geng T, Wan Z, Lu Q, Guo J, Liu L, et al. Serum selenium concentrations and risk of all-cause and heart disease mortality among individuals with type 2 diabetes. Am J Clin Nutr. (2022) 115:53–60. doi: 10.1093/ajcn/nqab241
Keywords: dietary selenium intake, osteoarthritis, NHANES, risk factor, arthritis
Citation: Deng X and Tan Y (2023) A national cross-sectional analysis of selenium intake and risk of osteoarthritis: NHANES 2003–2016. Front. Public Health 10:1047605. doi: 10.3389/fpubh.2022.1047605
Received: 23 September 2022; Accepted: 22 December 2022;
Published: 09 January 2023.
Edited by:
Kayo Masuko, Sanno Medical Center, JapanReviewed by:
Lutz Schomburg, Charité Universitätsmedizin Berlin, GermanyMohd Aleemuddin Quamri, National Institute of Unani Medicine (NIUM), India
Copyright © 2023 Deng and Tan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yongqiong Tan, MTEyNTg2NTQ3OCYjeDAwMDQwO3FxLmNvbQ==