Abstract
Introduction
In recent years, the global trend of population aging has become increasingly severe. An epidemiological study revealed that over 11% of the population is currently aged 60 years or older, and this proportion is expected to double by 2050. 1 The middle-aged and older adults suffering from osteoporosis is rapidly growing. Osteoporosis is primarily characterized by reduced bone mass and diminished bone strength, making affected individuals more susceptible to fractures in various parts of the body. This condition significantly increases the risk of disability and mortality, imposing a substantial burden on both individuals and society. Currently, treatment strategies for osteoporosis include the use of bisphosphonates, calcitonin, prophylactic vitamin D supplementation, and increased calcium intake. 2 Given its high prevalence and mortality rates, identifying potential risk factors and facilitating early diagnosis of osteoporosis are urgent priorities.
Major depressive disorder (MDD) is one of the most common mental health conditions, affecting approximately 16% of adults in the United States. 3 As depression progresses, patients often experience sleep disturbances, physical exhaustion, and abnormal limb sensations, which can lead to a range of complications, including cardiovascular and endocrine disorders, and in severe cases, suicidal behaviors. Middle-aged and older adults, being a vulnerable population both physically and mentally as well as socially, are at a heightened risk of depression due to factors such as chronic illnesses and cognitive impairments. 4 Current treatments for depression in the older adults primarily involve the use of antidepressant medications, coupled with enhanced social support and attention to their psychological well-being. Early recognition of depressive symptoms can aid in predicting the severity of the condition and the development of potential complications.
Previous studies have identified a correlation between depressive symptoms and osteoporosis, and this relationship has gained medical consensus among clinicians. 5 This association may extend to broader psychological conditions, as evidenced by Godde et al 6 using multi-wave HRS data, which established that a history of psychological disorders was significantly associated with higher odds of osteoporosis diagnosis in middle-aged and older Americans. Globally, numerous cohort studies have observed a significant association between depression and osteoporosis, noting that the risk of osteoporosis increases with the severity of depression and the use of certain antidepressant medications. 5 Current research on the mechanisms linking depressive symptoms to osteoporosis primarily focuses on three pathways: First, the neuro-endocrine pathway, where chronic stress in depression activates the hypothalamic-pituitary-adrenal (HPA) axis, leading to sustained elevation of glucocorticoids that directly suppress osteoblast activity 7 ; second, inflammatory mechanisms wherein depression-associated proinflammatory cytokines upregulate osteoclast differentiation 8 ; and third, behavioral mediation pathways (e.g., physical inactivity, smoking) affecting the bone microenvironment. 9 However, most of these studies have been conducted on relatively small sample sizes, highlighting the need for large-scale cohort studies to further explore the relationship between depressive symptoms and osteoporosis. Building upon the “Stress-Bone Integration Model” theoretical framework, this study employs large-scale cohort data to examine how chronic psychological stress influences bone homeostasis through multi-target neuro-immune-endocrine mechanisms, while exploring the modifying effects of individual differences such as age and sex.
The National Health and Nutrition Examination Survey (NHANES) is a large-scale ongoing study, while the Health and Retirement Study (HRS) is characterized by its cohort-based design, both offering valuable opportunities to investigate the relationship between depression severity and osteoporosis. These public databases include depression questionnaire scores and osteoporosis survey data, benefiting from the advantage of large sample sizes. Through statistical analysis and computation of these two extensive datasets, we aim to provide robust evidence for the association between depressive symptoms and osteoporosis.
Methods
Study Population
This dual-database study adheres to the ethical standards of the relevant national human experimentation committees and the principles of the Declaration of Helsinki (revised in 2013). The data were derived from a national cross-sectional study (NHANES 2005-2006, 2007-2008, and 2009-2010 cycles; protocol details available at: https://www.cdc.gov/nchs/nhanes/about/erb.html?CDC_AAref_Val=https://www.cdc.gov/nchs/nhanes/irba98.htm) and a cross-sectional analysis of the 2012 wave from the Health and Retirement Study (HRS) cohort (institutional review documentation: https://hrs.isr.umich.edu/sites/default/files/biblio/HRS_IRB_Information-10-2017.pdf), with the latter treated as cross-sectional data to ensure methodological consistency in comparing depression-osteoporosis associations between datasets.The Institutional Review Board (IRB) of the National Center for Health Statistics (NCHS) approved the collection of NHANES data, and the establishment of HRS was authorized by the National Institute on Aging (NIA U01AG009740) and the Social Security Administration.
This study follows the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guidelines (Supplemental Material).
This study analyzed participants aged ≥50 years from both datasets, applying this age threshold to NHANES (2005-2010) to align with HRS’s (2012) design as a population-based study specifically targeting adults aged 50 and older. We selected the 2012 HRS wave for methodological consistency with NHANES’ cross-sectional design and because it contained the most complete CES-D 8 and osteoporosis data among recent waves. Included participants had complete data on depression (PHQ-9 for NHANES; CES-D 8 for HRS), osteoporosis status, and key covariates (demographics, behaviors, comorbidities). After excluding individuals with missing data, the analytic samples comprised 3612 from NHANES (from 31 034 initially eligible) and 4307 from HRS (from 20 395 initially eligible) (Figure 1). The Flow Chart of Study Population Selection
Depressive Symptoms
In NHANES and HRS, depressive symptoms were assessed using the PHQ-9 and CES-D 8 depression scoring scales, respectively, which evaluate a series of emotions experienced by participants in daily life. The total scores for PHQ-9 and CES-D 8 are 27 and 8, respectively. A score of 10 or higher on the PHQ-9 and 4 or higher on the CES-D 8 was used to assign a value of 1 to the variable. The accuracy of these two depression scales in identifying depressive symptoms has been validated by related studies.10,11
Osteoporosis
The osteoporosis data in both NHANES and HRS were derived from identical self-reported questionnaire items that served as the primary outcome measure in this study. Specifically, the 2012 HRS wave - which first incorporated osteoporosis assessment - and NHANES used parallel questions: (1) “Has a doctor ever told you that you have osteoporosis?” followed by (2) “Are you currently taking medication to treat osteoporosis?” Participants affirming either question were classified as having osteoporosis. This approach for outcome ascertainment reflects primary care screening practices where initial assessment relies on clinical history before confirmatory testing.
Covariates
Continuous variables included age and waist circumference (WC), which were analyzed as continuous in all primary analyses (baseline descriptions and generalized linear models [GLMs] of depression-osteoporosis relationships). For subgroup analyses only, age and WC were converted to categorical variables.
Statistical Analysis
Statistical analyses were conducted using R software (version 4.4.1; R Foundation for Statistical Computing).13Cross-sectional data from NHANES were analyzed to identify associations and observations, followed by a comparison using cohort data from HRS. Descriptive statistics included means ± SD (normality verified by Shapiro-Wilk tests) and frequencies (%). Group comparisons used t-tests (with Levene’s test for equal variances) and χ2 tests (with Fisher’s exact test when expected counts <5). All statistical assumptions were explicitly tested and met.
After excluding participants with missing data on key variables (depressive symptoms, osteoporosis diagnosis, or covariates), we utilized GLMs to investigate the association between depressive symptoms and osteoporosis in NHANES and HRS data. We assumed a binomial distribution for the dependent variable with a logit link function. Despite using different depression scales (PHQ-9 in NHANES, CES-D8 in HRS), methodological consistency was maintained through uniform adjustment for all specified covariates and application of the same GLM framework to both datasets. Model diagnostics included examination of residuals, Cook’s distances, and variance inflation factors (all VIFs < 2), confirming no substantial outliers, multicollinearity, or model misspecification. Had overdispersion been detected, we would have employed quasibinomial regression to adjust standard errors. Covariate adjustments were performed in both datasets using three models. Model 1 included no adjustments. Model 2 adjusted for age, sex, education level, and marital status. Model 3 further included lifestyle factors (smoking, alcohol consumption, and physical activity) and comorbid conditions such as hypertension, diabetes, and arthritis. The selection of these covariates was primarily based on their established associations with both depression and osteoporosis in prior epidemiological studies.14,15Based on these models, depressive symptoms were dichotomized as present (PHQ-9 ≥10 in NHANES; CES-D-8 ≥4 in HRS) or absent, serving as the independent variable, and their association with osteoporosis (binary: yes/no) as the dependent variable was assessed using GLMs.
Additionally, we hypothesized that the association of depressive symptoms with osteoporosis might vary by individual characteristics. To test this, we conducted pre-specified subgroup analyses stratified by sex, age (<65 vs ≥65 years), diabetes status (yes/no), arthritis status (present/absent), smoking status (current/former/never), physical activity level (low/moderate/high), and WC (<90 vs ≥90 cm),
12
as these factors may modify the depression-osteoporosis relationship. In all interaction analyses, the GLMs included both main effects (depressive symptoms and each subgroup variable) and their multiplicative interaction terms (depressive symptoms × subgroup variable) to ensure proper model specification. This approach allows the assessment of whether the association between depressive symptoms and osteoporosis differs across subgroups while accounting for baseline effects of each variable.Statistical significance was defined as
Results
Characteristics of Participants in the NHANES(2005-2010) and HRS(2012)
Abbreviations: WC, waist circumference;SD, standard deviation.
Continuous variables (e.g., age, waist circumference): Mean (SD) reported; group comparisons by two-tailed Student’s t-test (
Cross-cohort comparisons revealed descriptive differences: NHANES participants were younger on average (65.1 years) than HRS participants (66.7 years), had a higher proportion with education beyond high school (67.0% vs 45.6%), and reported lower arthritis prevalence (47.1% vs 59.8%). It should be explicitly noted that these cross-cohort comparisons are based on unadjusted descriptive data without formal statistical testing between cohorts, and therefore require cautious interpretation.
Stratified Generalized Linear Models of Depressive Symptom Scores and Osteoporosis Associations
Abbreviations: OR,odds ratio; CI,confidence intervals.
Model adjustments:Model 1,no covariate was adjusted. Model 2,age,sex,marital status,education were adjusted.Model 3,age,sex,marital status,education,drinking status,smoking status,waist circumference,combination conditions(hypertension, diabetes, arthritis) and physical activity were adjusted.
Statistical tests: Two-tailed Wald tests were used to calculate
Stratified Generalized Linear Models of Depression Diagnosis and Osteoporosis Associations
Abbreviations: OR, odds ratio; CI, confidence intervals.
Model adjustments: Model 1, no covariate was adjusted. Model 2, age,sex,marital status,education were adjusted.Model 3, age, sex, marital status, education, drinking status, smoking status, waist circumference,combination conditions(hypertension, diabetes, arthritis) and physical activity were adjusted.
Statistical tests: Two-tailed Wald tests were used to calculate
Generalized Linear Model of Associations Between Depressive Symptom Quartiles and Osteoporosis
Abbreviations: OR, odds ratio;CI, confidence interval.
Model adjustments: Model 1, no covariate was adjusted. Model 2, age, sex, marital status, education were adjusted.Model 3, age, sex, marital status, education, drinking status,smoking status,waist circumference,combination conditions(hypertension, diabetes, arthritis) and physical activity were adjusted.
Statistical tests: Trend
The results of the subgroup analyses are presented in Figure 2. In the NHANES cohort, the relationship between depressive symptoms and osteoporosis demonstrated a broad positive correlation across subgroups such as age, sex, and diabetes status. When stratified by age, the <65 age group showed a significant increase in osteoporosis prevalence with worsening depressive symptoms (OR = 1.062, 95% CI:1.028-1.097). Stratification by hypertension status revealed a significant increase in the odds ratio for individuals with hypertension (OR = 1.074, 95% CI:1.041-1.107). Additionally, individuals with arthritis exhibited a significantly higher prevalence of osteoporosis (OR = 1.056, 95% CI:1.026-1.086).In the HRS cohort, similar widespread positive correlations were observed in the subgroup analyses. Stratification by diabetes status revealed that among individuals with diabetes, the prevalence of osteoporosis increased with greater depressive severity (OR = 1.113, 95% CI:1.016-1.219).The interaction analyses presented in Figure 2 indicate no statistically significant effect modifications in either the NHANES or HRS cohorts. Forest Plot for Subgroup Analysis of Depressive Symptoms and Osteoporosis. Abbreviations: WC,waist circumference;OR,odds ratio;CI,confidence Interval. Model adjustments:Age,sex,marital status,education,drinking status,smoking status,waist circumference,combination conditions(hypertension, Diabetes, Arthritis) and Physical Activity were Adjusted. Error Bars: 95% Confidence Intervals for Adjusted ORs.Sample Sizes: NHANES (n = 3612), HRS (n = 4307).Statistics:Subgroup ORs and 
Discussions
The association between depressive symptoms and osteoporosis requires further validation. In this study, we quantified depressive symptoms using the PHQ-9 and CES-D 8 depression scoring scales to explore their relationship with osteoporosis. While previous studies have examined depression-osteoporosis relationships using NHANES data, 16 our study offers the following strengths: (1) implementing harmonized analytical protocols across both NHANES and HRS cohorts to enhance comparability, and (2) conducting pre-specified, theory-driven evaluations of potential moderators including key comorbidities (diabetes, arthritis) and demographic factors. Previous studies using HRS and NHANES data have developed predictive models for osteoporosis diagnosis, 17 identifying mental health issues as predictors associated with higher osteoporosis risk - a finding consistent with our results. However, our study specifically focuses on elucidating the depression-osteoporosis association and investigating how various factors modify this relationship, rather than simply identifying osteoporosis-related factors.Our findings reveal a positive correlation between depressive symptoms and osteoporosis. When depressive symptoms were dichotomized, the association between depression and osteoporosis remained positive. Trend analysis demonstrated a linear relationship between depressive symptoms and osteoporosis. Subgroup analyses further confirmed the widespread association across different populations and explained the impact of specific stratifications on this relationship.
Osteoporosis signifies a reduction in bone density, which can primarily be attributed to three factors: (1) inadequate bone acquisition, (2) abnormal increases in bone resorption, and (3) decreased efficiency in bone remodeling. 5 Depression may lead to functional abnormalities in inflammatory mediators, as evidenced by elevated levels of pro-inflammatory factors or inflammatory markers such as IL-1β, IL-6, TNF-α, and CRP in depressed patients. These increases can reduce bone density and heighten fracture risk.18,19In addition to the influence of pro-inflammatory factors, depressed patients often experience neurohormonal imbalances, further disrupting bone metabolism. Increased activation of the HPA axis in depressed individuals results in elevated cortisol secretion, which inhibits osteoblast formation and promotes bone resorption. 20 Beyond neuroendocrine mechanisms, reduced secretion of gonadal hormones (estrogen and testosterone) and vitamin D deficiency in depressed patients contribute to bone loss and decreased bone density.21,22Antidepressant medications may also exacerbate bone loss. In a database study by Rajha, H.E., 23 the use of conventional antidepressants was associated with a 44% increase in osteoporosis prevalence. Efendioglu, E.M., et al. 24 reported that older adults with depression using selective serotonin reuptake inhibitors (SSRIs) showed poorer outcomes with anti-osteoporotic treatments. Furthermore, a meta-analysis highlighted that antidepressant use, particularly widely prescribed SSRIs, increased the risk of hip fractures by 2.5 times. 25 Adverse effects of antidepressants, such as balance disorders and dizziness, significantly elevate the risk of falls and fractures among the older adults.
Depressive symptoms are generally associated with osteoporosis prevalence among middle-aged and older adults. The causes of worsening depressive symptoms in middle-aged and older populations are multifaceted.Firstly, during the aging process, serotonin (5-HT) activity undergoes significant changes. The number of 5-HT2A receptors decreases markedly during the transition from youth to middle age. 26 Additionally, an animal study found that the binding rate of 5-HT1A receptors declines in older adults, potentially explaining the reduced efficacy of antidepressants in older adults with depression. 27 Endocrine factors also play a crucial role in both late-life depression and osteoporosis. Beyond the previously mentioned gonadal hormones, alterations in levels of hormones such as growth hormone and parathyroid hormone can influence bone remodeling processes, either synergistically or independently.28,29Moreover, with advancing age, chronic diseases such as diabetes and arthritis gradually affect overall physical health and mental well-being, exacerbating the severity of depressive symptoms.
Diabetes is one of the most common chronic diseases among middle-aged and older adults. Intriguingly, our stratified analyses across both cohorts demonstrated that depressive symptoms were linked to osteoporosis specifically in individuals with diabetes. Diabetic patients often experience insulin resistance, which is associated with a higher prevalence of mental health issues. An animal study demonstrated that mice lacking insulin receptors in astrocytes exhibited impaired dopamine release mechanisms, leading to depressive-like behaviors. 30 While insulin and insulin-like growth factor-1 (IGF-1) may alleviate depressive symptoms, receptor antagonists for these pathways can obstruct these beneficial effects. 31 On the other hand, diabetic patients are often in a metabolically stressed state, with chronic activation of the HPA axis. This leads to hypercortisolism, which contributes to both insulin resistance and depressive symptoms.32,33 Additionally, structural and functional abnormalities in the brain, as well as systemic inflammation, further link diabetes to depression. 34 The effects of diabetes on bone density remain inconclusive. While many studies suggest that insulin plays a critical role in bone metabolism, osteoporosis is recognized as a complication of type 1 diabetes, whereas reduced bone density is not typically observed in type 2 diabetes. 35 However, Vestergaard, P., et al. 36 found that both types of diabetes were associated with a higher risk of fractures compared to healthy individuals. This may indicate that blood glucose levels are a critical factor, although the underlying mechanisms require further investigation.
The impact of arthritis on depressive symptoms and osteoporosis was also significant.Arthritis often negatively affects the mental health and musculoskeletal health of middle-aged and older adults.The immune regulatory mechanisms and pro-inflammatory factors associated with rheumatoid arthritis (RA) may influence relevant neurotransmitters, suggesting that RA and depression may share similar pathogenic mechanisms. 37 While the pathophysiological links between osteoarthritis (OA) and depression, beyond immune and inflammatory factors, remain unclear, their association has been well established by numerous studies.38,39Arthritis, as a musculoskeletal disease, is closely related to osteoporosis through various pathways, including the disease itself, treatment methods such as corticosteroid therapy, and its impact on patients’ lifestyle habits, such as prolonged sedentary behavior or bedrest.40,41
The sex-specific characteristics observed in the middle-aged and older adults in this study warrant attention. The influence of sex on depression and bone density remains inconclusive. Most studies suggest that postmenopausal women, due to hormonal changes, are more susceptible to both depressive symptoms and osteoporosis compared to men. However, Mussolino, M.E. et al. 42 reported that severe depression in younger men was significantly associated with reduced bone density, a phenomenon not observed in women. Meanwhile, a prospective cohort study by Whooley, M.A. et al. 43 on men aged 50 and older found no association between depression and osteoporosis.Depression may increase osteoporosis risk in older men by reducing bioavailable testosterone - a mechanism supported by documented links between: (1) testosterone-depressive symptoms inverse correlation, 44 and (2) testosterone-bone density positive association. 45
These discrepancies across studies may stem from differences in baseline characteristics of the populations studied, bone density measurement methods, and definitions of depression. To understand this sex-specific difference, we analyzed potential mechanisms. From a behavioral perspective, men exhibit high-risk behaviors such as smoking and alcohol consumption, which are more frequent compared to women and are recognized risk factors for osteoporosis. 46 Biologically, cortisol may play a critical role in this disparity. Given the involvement of the HPA axis in these associations, we draw a comparison to Cushing’s syndrome, which is characterized by elevated cortisol levels. Men with Cushing’s syndrome have a higher prevalence of osteoporosis than women, 47 and studies have shown that vertebral fractures occur more frequently in men with this condition. 48 While the mechanisms of depression and Cushing’s syndrome differ, and the degree of cortisol elevation is not identical, these findings offer another perspective on sex differences.Additionally, the role of neuropeptide Y (NPY) warrants consideration. NPY, widely distributed in the nervous system, plays a crucial role in maintaining homeostasis and regulates osteoblasts and chondrocytes through central and peripheral pathways, forming a neuro-bone metabolic axis. 49 NPY also plays a key role in the pathogenesis of depression, 50 with Ishida, H. et al finding that NPY injection into the hippocampal CA3 region produced antidepressant effects. 51 Furthermore, a meta-analysis indicated that women taking psychotropic medications exhibited higher levels of NPY compared to men. 52 This suggests that elevated NPY levels may mitigate the risk of bone density loss in women with depression, providing a possible explanation for the observed sex differences.
This study utilized WC as a clinically superior anthropometric measure for obesity evaluation, 12 particularly valuable in assessing central adiposity and metabolic risks where body mass index (BMI) data availability was limited.We observed a significant positive correlation between depressive symptoms and osteoporosis in populations with larger WC. As a key indicator of abdominal obesity, increased WC was closely associated with higher all-cause mortality and metabolic disease risks in middle-aged and older adults. 53 Recent genetic studies have demonstrated genomic colocalization between depression-associated loci and obesity-susceptibility genes. 54 For instance, the neuronal growth regulator 1 gene (NEGR1) influences energy metabolism through hypothalamic appetite regulation networks. 55 Mechanistically, obesity-related chronic low-grade inflammation promotes the release of proinflammatory cytokines and upregulates NLRP3 inflammasome expression, leading to impaired glucocorticoid receptor function, dysregulated cortisol metabolism, and HPA axis overactivation56-58 - these neuroendocrine abnormalities have been previously shown to correlate with decreased bone mineral density (BMD).Notably, beyond typical obesity, sarcopenic obesity represents a distinct phenotype where reduced skeletal muscle mass coupled with dysregulated lipid metabolism significantly disrupts bone homeostasis. 12 Currently, precision intervention strategies for such musculoskeletal comorbidities remain to be further explored.
In conclusion, the association between depression and osteoporosis in middle-aged and older adults requires attention to the following aspects: (1) (2) (3)
These findings highlight the complexity of the interaction between depression and osteoporosis and underscore the importance of a multifaceted approach to understanding and addressing this relationship.
This study included middle-aged and older adults aged ≥50 years from the NHANES and HRS databases, benefiting from a large sample size and confirming the association between depressive symptoms and osteoporosis. The findings contribute to osteoporosis research and prevention by emphasizing the importance of monitoring bone density in middle-aged and older adults with depression. The deepening of depressive symptoms may serve as a predictive factor for osteoporosis. The use of depression scales and self-reported disease diagnoses offers a quick and cost-effective approach. However, this study also has several limitations. First, as a cross-sectional study, it cannot establish a causal relationship between depressive symptoms and osteoporosis.We plan to address this limitation through future longitudinal cohort studies to better characterize temporal associations. Second, the lack of pre-study power analysis may have limited our ability to detect subtle associations. Our variable selection approach was less rigorous than the systematic filtering methods employed by Godde et al.
Conclusion
This study demonstrates a significant association between depressive symptoms (assessed using validated PHQ-9 and CES-D8 scales) and increased osteoporosis risk in middle-aged and older adults, with consistent findings across diverse subpopulations. While the cross-sectional design precludes causal conclusions, the observed dose-response relationship suggests depressive symptoms may adversely affect bone health. These findings highlight three key clinical needs: (1) integrating bone density monitoring into depression management, (2) incorporating depression screening in osteoporosis care, and (3) developing coordinated prevention strategies for these comorbid conditions. At the public health level, our results emphasize the importance of addressing mental health in aging populations through enhanced social support and targeted interventions. Future longitudinal studies should clarify the temporal relationship and potential bidirectional mechanisms between these conditions.
Supplemental Material
Supplemental Material - Depressive Symptoms and Osteoporosis in Middle-Aged and Older Adults: A Cross-Sectional Analysis of NHANES and HRS Data
Supplemental Material for Depressive Symptoms and Osteoporosis in Middle-Aged and Older Adults: A Cross-Sectional Analysis of NHANES and HRS Data by Wen Zhang, Yi Tang, Lei Chen, Zhongyi Zhang, Xinyu Hu, Kai Cheng, Jiaju Zhou, and Peijian Tong in Geriatric Orthopaedic Surgery & Rehabilitation
Footnotes
Author Note
Ethical Considerations
Author Contributions
Funding
Declaration of Conflicting Interests
Data Availability Statement
Supplemental Material
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
