The Effect of Individual Musculoskeletal Conditions on Depression: Updated Insights From an Irish Longitudinal Study on Aging

Few longitudinal studies have systematically investigated whether or how individual musculoskeletal conditions (IMCs) convey risks for negative psychological health outcomes, and approaches to assess such risk in the older population are lacking. In this Irish nationally representative longitudinal prospective study of 6,715 individuals aged 50 and above, machine learning algorithms and various models, including mediation models, were employed to elaborate the underlying mechanisms of IMCs leading to depression and to develop an IMC-induced negative psychological risk (IMCPR) classification approach. Resultantly, arthritis [odds ratio (95% confidence interval): 2.233 (1.700–2.927)], osteoporosis [1.681 (1.133–2.421)], and musculoskeletal chronic pain [MCP, 2.404 (1.838–3.151)] were found to increase the risk of depression after 2 years, while fracture and joint replacement did not. Interestingly, mediation models further demonstrated that arthritis per se did not increase the risk of depression; such risk was augmented only when arthritis-induced restrictions of activities (ARA) existed [proportion of mediation: 316.3% (ARA of usual), 213.3% (ARA of social and leisure), and 251.3% (ARA of sleep)]. The random forest algorithm attested that osteoarthritis, not rheumatoid arthritis, contributed the most to depressive symptoms. Moreover, bone mineral density was negatively associated with depressive symptoms. Systemic pain contributed the most to the increased risk of depression, followed by back, knee, hip, and foot pain (mean Gini-Index: 3.778, 2.442, 1.980, 1.438, and 0.879, respectively). Based on the aforementioned findings, the IMCPR classification approach was developed using an interpretable machine learning model, which stratifies participants into three grades. Among the IMCPR grades, patients with a grade of “severe” had higher odds of depression than those with a “mild” [odds ratio (95% confidence interval): 4.055 (2.907–5.498)] or “moderate” [3.584 (2.101–5.883)] grade. Females with a “severe” grade had higher odds of depression by 334.0% relative to those with a “mild” grade, while males had a relative risk of 258.4%. In conclusion, the present data provide systematic insights into the IMC-induced depression risk and updated the related clinical knowledge. Furthermore, the IMCPR classification approach could be used as an effective tool to evaluate this risk.

Introduction

Human life expectancy worldwide has increased at an unprecedented rate during the last 200 years, and age can no longer be shaped by natural selection (1, 2). Population aging escalates the global burden of diseases, which makes explorations for reducing the length and severity of late-life morbidity a major goal (2). Psychiatric disorders contribute significantly to the global disease burden, and depression is the most common psychiatric disorder, which increases the disease burden (3, 4). Depression remains a major health problem in humans (5). Globally, over 300 million people and more than 25% of older people are compromised by depression, which severely disrupts psychosocial functioning and diminishes their quality of life (6–8). This disorder is deemed a leading cause of disability worldwide by the World Health Organization, and it often leads to tragic consequences, including suicide (9–11). Indeed, major depressive disorder as a prime cause of years lived with disability (YLD) contributes to 4.2% of total YLD worldwide (4); an existing depression longitudinal study also discovered that the overall long-term suicide risk varied from 5.6 to 6.8%, and the suicide risk of men with major depression reached as high as 20% (12). In addition, in the context of global turbulence (especially considering the recent outbreak of COVID-19), the socioeconomic burden of this debilitating disorder will continue to increase significantly (7, 13). From 2011 to 2030, the cumulative impact of mental-associated material expenditures on the global economy was estimated at US$16.3 trillion (13). This sobering and salient evidence indicates that it is extremely important to explore the factors that affect depression to control or prevent this disease.

Musculoskeletal conditions contribute even more to the global burden of disease among older people than do neurological and mental disorders (7.5 vs. 6.6%) (14); unfortunately, the two disorders are often entangled. Globally, over 1.3 billion people suffer from musculoskeletal disorders, which are a substantial contributor to non-fatal health loss (4). For example, osteoporosis escalates the propensity of fragility fractures (15), for which a fracture of the hip or spine carries a 12-month excess mortality rate of up to 20% (16). Osteoarthritis, as a major cause of disability worldwide, often results in pain, comorbidity, and mortality (17, 18). YLD was augmented by 54% between 1990 and 2015 owing to low back pain (19). The worst part is that musculoskeletal conditions confer mental disorders, including depression (20–23), which seriously threatens the quality of life of older people. However, how musculoskeletal conditions affect psychiatric health remains elusive.

To systematically interpret the roles and underlying mechanisms of musculoskeletal conditions on psychopathology risk, we relied on various models, including the mediation model and random forest algorithm (24–26), focusing on four classes of associations: arthritis-induced restriction of activities, arthritis subtypes, bone mineral density, and severe musculoskeletal chronic pain (MCP). After these works were completed, we re-examined the musculoskeletal conditions affecting psychopathological health and derived an IMCPR classification approach for negative psychological risk stratification of older people with musculoskeletal diseases.

Materials and Methods

STROBE recommendations were employed to guide this reporting (27).

Participants

We drew from The Irish Longitudinal Study on Aging (TILDA), a longitudinal study with two waves of population-representative data from Irish participants aged 50 years and above and those that were sampled using geographic-clusters-based RANSAM sampling system (28, 29). The study targeted individuals aged 50 and older living in the Republic of Ireland in residential accommodation (excluding those living in nursing homes and other institutions). Full informed consent was given by 8,504 participants, and Trinity College Dublin Research Ethics Committee issues the Ethical approval.

From October 2009 to February 2011, participants were contacted for a Wave 1 interview (30), and a Wave 2 interview was conducted from April 2012 to January 2013 (31). The response rate adjusted for eligibility is 62.0%. The cohorts for this study were recruited obeying the process shown in Supplementary Figure 1.

Measures

Demographic data (i.e., gender and age, place of residence, smoking status, and marital status) were collected via a report at Wave 1.

Center for Epidemiological Studies Depression Scale (CES-D) was employed for depressive symptoms assessment. The CES-D is a 20-item scale assessing four major factors: depressed affect, absence of positive affect or anhedonia, somatic activity, or inactivity and interpersonal challenges (32). The response values are four-point Likert scales ranging from 0 to 3, yielding a total possible score of 60, with higher scores indicating greater depressive symptoms. When applied to the elder population, the CES-D scale works very well; an accumulated value of 16 and above was confirmed a sensitivity of 92% and a specificity of 87% for diagnosis of late-life depression (LLD) (33); the reliability coefficient was 0.85 to 0.91 for older people (34).

Individual musculoskeletal conditions were measured at Wave 1. Parameters collected includes arthritis (subtypes and degree of arthritis-induced restriction of activities), osteoporosis (T-score and stiffness), fracture, joint replacement, and MCP (most severe body part of MCP).

Osteoporosis increases the risk of fracture (35). An osteoporosis participant with a newly occurred fracture at Wave 2 (without fracture at Wave 1, but with fracture at Wave 2) was considered to be suffering from osteoporosis-related fracture.

Objective measures of physical health were collected in the fixed health assessment centers. These measures include T-score and stiffness (measured via heel ultrasound), triglyceride (TG), low-density lipoprotein cholesterol (LDL), high-density lipoprotein cholesterol (HDL), total cholesterol (CHOL), and body mass index (BMI).

Covariates were collected from theoretical, practical, or previous empirical evidence of association with musculoskeletal disorders or depression (36–39). Age, gender, place of residence, smoking status, marital status, chronic illnesses self-reported (i.e., an abnormal heart rhythm, congestive heart failure, angina, myocardial infarction, and diabetes), BMI, TG, LDL, HDL, CHOL, cognitive score, use of medication [i.e., benzodiazepines, antihypertensives (ATC C09, C08, C07, C03, C02), and anticholinergics (urological, antidepressant, other)], and alcohol were included.

Statistical Analyses

Binomial logistic regression quantified the associations [i.e., odds ratios (OR)] between arthritis, osteoporosis, fracture, joint replacement, MCP, and probable LLD at baseline and Wave 2. Hosmer–Lemeshow test was conducted, and Nagelkerke R² was subsequently calculated to evaluate the goodness of fit of models (40). Likelihood ratio tests assessed covariate significance (41).

For arthritis, arthritis-induced restriction of activities (ARA) was considered the potential mediators for mediation analysis (26). Kruskal–Wallis one-way ANOVA test (H-test, a pairwise comparison using Dwass–Steel–Crichtlow–Fligner test, and p-value adjustment using Benjamini and Hochberg method) revealed the potential difference of depressive symptoms among the different degrees of ARA (42). Stratified analyses were conducted to quantify associations between no arthritis, arthritis with/without ARA, and LLD. Gini Index calculated by random-forest algorithm (1,000 grew decision trees) quantified the contribution of subtypes (26).

For osteoporosis, osteoporosis-related fracture as a potential mediator participated in mediation analysis. Spearman's rank correlation coefficient assessed the correlations between bone mineral density (T-score and stiffness) and depressive symptoms. For MCP, the random-forest algorithm was performed to quantify contributions of the most severe body part of MCP.

Potential differences of depressive symptoms in musculoskeletal conditions were investigated by the Wilcoxon test, and H-test measured the gender-related potential difference.

IMCPR classification approach was constructed by decision tree regression referring to existing work (43), and musculoskeletal conditions associated with LLD were included for analysis. Associations between IMCPR grades and LLD were quantified (also considering gender-related associations).

R (version 3.6.3) was employed for statistical analyses. Packages used in this study are as follows: Hosmer–Lemeshow test (ResourceSelection package); likelihood ratio test (lmtest package); Wilcoxon test, Kruskal–Wallis one-way ANOVA test, Dwass–Steel–Crichtlow–Fligner test, Benjamini and Hochberg method used for p-value adjustment, and corresponding effect sizes calculation (ggstatsplot package); Random-Forest algorithm (randomForest package); and decision tree regression (rpart package).

Results

Demographics of the Participants at Baseline

The demographics of the Irish participants aged 50 and above are shown in Table 1. Among the musculoskeletal conditions, four (arthritis, osteoporosis, fracture, and MCP) had statistical significance between the no-depression and depression groups. Specifically, they were arthritis (p < 0.001, V = 0.083, 95% CI: 0.057–0.109, small effect size), osteoporosis (p < 0.001, V = 0.046, 95% CI: 0.014–0.073, small effect size), fracture (p = 0.007, V = 0.031, 95% CI: 0.004–0.057), and MCP (p < 0.001, V = 0.170, 95% CI: 0.146–0.195, small effect size).

TABLE 1

Table 1. Demographics of Irish participants aged 50 and above on baseline.

Preliminary Cross-Sectional and Longitudinal Results

To assess the effects of musculoskeletal conditions on the development of depression, binomial logistic regression was conducted to quantify the associations (i.e., ORs) between the musculoskeletal conditions and probable LLD at baseline and Wave 2 (Figure 1).

FIGURE 1

Figure 1. Preliminary cross-sectional and longitudinal results. Model 1 adjusted for gender and age. Model 2 adjusted for Model 1 and place of residence, smoking status, and marital status. Model 3 adjusted for Model 2 and an abnormal heart rhythm, congestive heart failure, angina, myocardial infarction, BMI, diabetes, TRIG, LDL, HDL, CHOL, cognitive score, use of medication [i.e., benzodiazepines, antihypertensives (ATC C09, C08, C07, C03, C02), and anticholinergics (urological, antidepressant, other)], and alcohol. OR, Odds ratios; CI, confidence intervals; p, p-value of binomial logistic regression result; P_H−L, p-value of Hosmer–Lemeshow test of goodness of fit; P_L−R, p-value of likelihood ratio test; M1:M2, likelihood ratio test conducted between model 1 and model 2; M2:M3, likelihood ratio test conducted between model 2 and model 3.

For the cross-sectional studies, the prevalence of probable LLD was 8.55% (n = 429). Among all of the musculoskeletal conditions, arthritis, osteoporosis, and MCP were confirmed to be significantly associated with probable LLD. The crude model attested that arthritis, osteoporosis, and MCP were associated with 76.2% (OR = 1.762, 95% CI: 1.434–2.160; p < 0.001), 55.4% (OR = 1.554, 95% CI: 1.158–2.055; p = 0.003), and 195.1% (OR = 2.951, 95% CI: 2.412–3.618; p < 0.001) higher odds of probable LLD. Sensitivity analyses further substantiated the robustness of these associations. The OR, as the quantitative indicator of associations between arthritis, osteoporosis, MCP, and depression, fluctuated a little but was generally stable (Figure 1).

For the longitudinal studies, the prevalence of probable LLD was 4.92% (n = 226). As we expected, arthritis, osteoporosis, and MCP were reaffirmed to be significantly associated with probable LLD. Arthritis, osteoporosis, and MCP were associated with 123.3% (OR = 2.233, 95% CI: 1.700–2.927; p < 0.001), 68.1% (OR = 1.681, 95% CI: 1.133–2.421; p = 0.007), and 140.4% (OR = 2.404, 95% CI: 1.838–3.151; p < 0.001) higher odds of probable LLD in the crude model. In further sensitivity analyses, the robustness of these associations was confirmed (Figure 1).

Comprehensively considering the results of the goodness-of-fit test and the likelihood ratio test (Figure 1), model 1, model 2, and model 3 were deemed different during the sensitivity analyses. Thus, the robustness of the associations between arthritis, osteoporosis, MCP, and LLD was validated using sensitivity analyses (44).

Arthritis to LLD

To determine whether the arthritis-induced restriction of activities (ARA) explained the increased depression risk, we tested a series of mediation models. Interestingly, our results were somewhat contrary to the results observed in other works, which simply concluded that arthritis was significantly linked to enhanced risks for depression (45, 46); however, the causal pathways between arthritis and depression are far from simple. The present results newly discovered and elaborated one of the potential causal pathways between arthritis and depression, as shown in Figure 2A. In cross-sectional mediation analyses, the direct effects of arthritis attenuated the risk for depression, while ARAs presented completed mediation effects that mediated 316.3% (ARA of usual), 213.3% (ARA of social and leisure), and 251.3% (ARA of sleep). The results of the longitudinal studies were very similar. Arthritis-induced direct effects attenuated or did not affect the depression risk, and ARAs demonstrated mediation effects that mediated 112.0% (ARA of usual), 132.3% (ARA of social and leisure), and 74.6% (ARA of sleep). Further analyses of potential differences in depressive symptoms corroborated the results of the mediation analyses. No statistically significant difference was detected between participants without arthritis and arthritis participants without ARA. In contrast, depressive symptoms gradually increased with the increasing degree of ARA (p < 0.001, middle effect sizes), regardless of the cross-sectional (Figure 2C) or longitudinal design (Figure 2D). Stratified analyses using binomial logistic regression reaffirmed that no significant association existed between arthritis without ARA and probable depression (Figure 2E). Considering that the participants without arthritis and the arthritis participants without ARA had no significant difference in depressive symptoms and that their group means were similar, the two were classified together for further stratified analyses. Arthritis with ARA was associated with 141.1% (OR = 2.411, 95% CI: 2.011–2.884 for cross-sectional) and 126.8% (OR = 2.268, 95% CI: 1.755–2.911 for longitudinal) higher odds of probable depression (vs. no arthritis and arthritis without ARA). When further restricted to participants with arthritis, arthritis with ARA was associated with 200.8% (OR = 3.008, 95% CI: 2.106–4.414 for cross-sectional) and 66.1% (OR = 1.661, 95% CI: 1.104–2.564 for longitudinal) higher odds of probable depression.

FIGURE 2

Figure 2. Further analyses for restriction of activities, and subtypes in arthritis. (A) Results of cross-sectional and longitudinal mediation analysis for arthritis, in which ARA (usual), ARA (social and leisure), and ARA (sleep) were considered as potential mediators. (B) Random-Forest algorithm-derived Mean Gini Index weighed the contribution of arthritis subtypes to outcomes of depression. (C) Cross-sectional potential difference of depressive symptoms among participants with/without arthritis (stratified by degree of ARA). (D) Longitudinal potential difference of depressive symptoms among participants with/without arthritis (stratified by degree of ARA). (E) Stratified analyses using binomial logistic regression quantified associations between arthritis (No ARA), arthritis (ARA), and depression at baseline and Wave 2. ARA, Restriction of activities caused by arthritis; ACME, average causal mediation effects; ADE, average direct effects; Prop. mediated, proportion mediated; χ², chi-square; n_obs, number of observation; û, group mean; $\widehat{{\epsilon }^2},$effect size of ε²; p, p-value. p-value adjusted < 0.05 was represented by “*,” p-value adjusted < 0.01 was represented by “**,” p-value adjusted < 0.001 was represented by “***”.

Subsequently, the contributions of arthritis subtypes to depressive symptoms were investigated in Wave 2. As a result of the random forest displayed in Figure 2B, osteoarthritis contributed the most to the increased depression risk, followed by two or more arthritis, rheumatoid arthritis, and other arthritis (the mean Gini index was 1.089, 0.860, 0.713, and 0.352 for longitudinal studies).

Osteoporosis to LLD

Fragility fractures are the most common and severe complications of osteoporosis (47–49); therefore, we tested whether they mediated the increased risk for LLD. Mediation analysis found that osteoporosis-associated fractures that newly occurred at Wave 2 did not mediate the association between osteoporosis and LLD (Supplementary Table 1). T-score and stiffness as indicators for osteoporosis were negatively correlated with depressive symptoms, regardless of baseline (Figures 3A,D) or Wave 2 (Figures 3B,E).

FIGURE 3

Figure 3. Further analyses for osteoporosis and MCP. (A) Association between T-score and depressive symptoms at baseline. (B) Association between T-score and depressive symptoms at Wave 2. (C) Mean Gini Index quantified the contributions of different body parts with the most severe MCP. (D) Association between stiffness and depressive symptoms at baseline. (E) Association between stiffness and depressive symptoms at Wave 2. (F) Association between the most severe body part of MCP and depressive symptoms at Wave 2.

MCP to LLD

We investigated whether the difference in the most severe body part of MCP affects the increased risk for LLD at Wave 2. As a result of the random forest presented in Figure 3C, systemic pain contributed the most to the increased risk for LLD, followed by back, knees, hips, and feet (the mean Gini index was 3.778, 2.442, 1.980, 1.438, and 0.879, respectively). It was significantly associated with depressive symptoms (Figure 3F).

Potential Difference of Depressive Symptoms at Wave 2

As we expected, participants with arthritis, osteoporosis, or MCP tended to report more severe depressive symptoms (arthritis: p < 0.001, $\widehat{r}$ = 0.109, 95% CI: 0.082–0.130, small effect size; osteoporosis: p = 0.010, $\widehat{r}$ = 0.040, 95% CI: 0.014–0.068, small effect size; MCP: p < 0.001, $\widehat{r}$ = 0.169, 95% CI: 0.146–0.194, small effect size) (Supplementary Figure 2). Among the three disorders, females were always more likely to report more severe depressive symptoms than males (Supplementary Figure 3). These resultant data indicate that gender-related differences in depressive symptoms were non-negligible.

IMC-Induced Negative Psychological Risk Classification Approach

Musculoskeletal disorders (arthritis, osteoporosis, and MCP) significantly associated with LLD were included for decision tree regression. The derived IMCPR classification approach, including arthritis and MCP, could stratify older people into three grades of mild, moderate, and severe (Figure 4A). Details of the IMCPR classification approach are given in Supplementary Table 2. Within the tree regression analysis of the overall cohort, the most significant factor conveying risk for LLD was arthritis. The prevalence of probable LLD was 3.521% (mild), 6.203% (moderate), and 12.891% (severe). Analyses for depressive symptoms reported significant differences among the three grades (p < 0.001, $\widehat{{\epsilon }^2}$ = 0.035, 95% CI: 0.029–0.047, small effect size) (Figure 4B), and the significant difference remained when further considering gender (p < 0.001, $\widehat{{\epsilon }^2}=$ 0.042, 95% CI: 0.032–0.055, small effect size) (Figure 4C). Of particular note, females always tended to report more severe depressive symptoms than males. Overall, participants who were classified as “severe” had 305.5%/123.9% higher odds of LLD than “mild”/“moderate” and “moderate” participants, who possessed 81.2% higher odds of LLD than “mild” participants (Figure 4D). Gender differences must be recognized. The results of the stratified analyses for females agreed with this tendency of OR distribution, but for males, it partly did not (Figure 4D). Males had higher odds of LLD only when they were graded as “severe.” Specifically, females who were classified as “severe” had higher odds of depression by 334.0% relative to those classified as “mild,” while males who were classified as “severe” only had higher odds of depression by 258.4% relative to those classified as “mild”.

FIGURE 4

Figure 4. IMC-induced negative psychological risk classification approach. (A) IMCPR classification approach was derived using tree regression. (B) Difference of depressive symptoms among participants with different IMCPR grades. (C) Potential gender-related difference of depressive symptoms among participants with different IMCPR grades. (D) Stratified analyses using binomial logistic regression derived odds ratios (OR) and 95% confidence intervals (CI) as indicators of association between IMCPR grades and depressive disorder, in overall, female, or male cohort. NA, No arthritis.

Discussion

This prospective study on data from a population–representative sample of older people living in the same context represents, to our knowledge, the first attempt to systematically investigate whether musculoskeletal conditions and surgery are associated with increased risks for negative psychological health outcomes and to further explore the underlying mediators. It is also the first large longitudinal study using machine learning that sought to derive a classification approach to assess the potential augmented psychiatric risks of older people who suffered from musculoskeletal disorders.

We investigated not only whether but also how individual musculoskeletal conditions that commonly occur in older people, such as arthritis, osteoporosis, fracture, joint replacement operations, and MCP, could confer an enhanced risk for LLD. One important finding revealed by the present study was that it is not the original arthritis itself but rather the arthritis-induced restriction of individual activities that plays a prominent role in the development of psychiatric disorders in these elderly populations. The present results updated the recognition generalized in other works, which simply deemed arthritis as a hazardous factor resulting in depression (45, 46); however, the causal pathways between arthritis and depression were underappreciated. The resultant data found that arthritis did not convey a risk for negative psychological health outcomes unless it restricted individual activities (usual, social and leisure, or sleep). Arthritis sufferers with ARA augmented the odds of probable LLD after 2 years by 126.8%, while arthritis sufferers without ARA were not compromised by such risks. Restriction of activities completely mediated the association between arthritis and depression (ARA: 316.3% of usual, 213.3% of social and leisure, and 251.3% of sleep), and depressive symptoms gradually augmented with the escalating degree of ARA.

In the elderly subjects, the findings of arthritis collected in the present study are in line with previous work (50–52), demonstrating that arthritis may contribute to depressive disorders/symptoms over time. The present study provided further supporting evidence to clinicians who may transfer the focus from arthritis itself to the arthritis-induced restriction of activities as a significant trigger for medical interventions among elderly sufferers. The assessment scale of the degree of ARA is enclosed in this study for potential medical applications. Moreover, rheumatoid arthritis is generally considered a critical factor in the development of depression in the reported studies (45, 53, 54), which is contrary to our findings. The random forest algorithm confirmed that among the subtypes of arthritis, osteoarthritis contributed the most to the increased risk for LLD (the Gini Index was 1.089 for longitudinal studies), which may be explained by osteoarthritis being a common and disabling condition that affects the daily functioning and independence of older people (55, 56). It was reported that osteoarthritis, as the most prevalent subtype of arthritis in older people, contributes more to YLD than rheumatoid arthritis (4). Obviously, osteoarthritis and its underlying mechanisms leading to an increased risk for psychological disorders warrant further investigation.

Our findings highlight the importance of arthritis-induced restriction of activities, which completely mediated the association between arthritis and depression. Activity restrictions are strongly associated with depression in older people (57). Unfortunately, arthritis is always inseparable from the restriction of activities (56). Given the high incidence of arthritis in older people, it is remarkably important to emphasize the symptoms that can actually lead to depression.

Osteoporosis and MCP often convey a risk for depressive disorder (58, 59), which was corroborated in the present study on population–representative samples. In addition, we found that osteoporosis-related fracture, as the most common complication of osteoporosis, did not mediate the enhanced risk for LLD. More importantly, it is noteworthy that such severe, acute trauma should not be considered a potential risk factor for LLD according to our findings. T-score and stiffness, as diagnostic clinical parameters for osteopenia and osteoporosis (60), were negatively associated with depressive symptoms. This result indicated that depressive symptoms were augmented with the increasing severity of osteoporosis (decreased bone mineral density), which could serve as a potential indicator for increased risk for LLD in the development of osteoporosis. For MCP, we attested that different severe body parts of MCP convey different degrees of depressive symptoms, which might be conducive to guiding stratified interventions for older people with MCP.

The most important aspect of the present study lies in the IMC-induced negative psychological risk classification approach for the prediction and quantification of increased risks for depression. The IMCPR classification approach was derived by decision tree regression from all musculoskeletal disorders that were confirmed to be associated with an augmented risk for probable LLD in our study. Our results demonstrated that older people with a higher grade of IMCPR suffered an augmented risk for worse psychological health outcomes. In detail, a patient with a grade of “severe” had 305.5/123.9% higher odds of depression than those with a “mild”/“moderate” grade.

These findings might offer insight into the clinical treatments of older people who suffer from musculoskeletal disorders and highlight the importance of refined stratified psychological interventions to ensure the quality of life of the older population. Older people classified at lower grades of IMCPR would not require excessive psychological intervention, while those classified as higher grades demand comprehensive interventions. Musculoskeletal disorders, which generally compromise older people, exacerbate the global disease burden and consume enormous socioeconomic costs (4, 14). Aided by our IMCPR classification approach, older people with higher odds of LLD could be regularly identified; such a classification approach could not only reduce unnecessary psychological interventions but also potentialize precision medicine and optimize the redistribution of social-medical resources.

Gender-related differences in depressive disorder/symptoms were also taken into consideration in the present study. Among musculoskeletal disorders, females always tend to report higher odds of depression or poorer depressive symptoms than males. Such gender-related differences remained in the IMCPR classification. Females had two- to three-fold higher rates of depression than males, as previously reported (61), of which several explanations have been offered: (1) Sexual dimorphism at the transcriptional level mediates gender-specific stress susceptibility (62). (2) Gender-related differences in neurotransmitter, endocrine, and metabolic systems result in differences in mood susceptibility (61). Our results from the stratified analyses revealed that, as we hypothesized, females had higher odds of probable LLD than males in the older population who were classified at the same IMCPR grade. Females with a “severe” grade had higher odds of LLD by 334.0% relative to those with a “mild” grade, while males had higher odds of LLD by only 258.4%. The proposed IMCPR classification approach considering gender-related differences has huge potential for clinical application because it was confirmed as outstanding for negative psychological risk stratification referring to different genders.

We acknowledge that this study has limitations. Chiefly, limited to the methods and original data of large-scale cohort research (TILDA), both activity restriction and the IMCPR classification approach are based on self-reported symptoms, and the present IMCPR classification approach demands further validation in clinical practice as already planned. For the same reason, we could not evaluate the effects of frailty score, mobility scores, step count, and vitamin D status. These indicators will be adopted in the follow-up studies after further optimized design with the inclusion of more extensive and detailed source data. Furthermore, the genetic mechanism of ARA-mediated effects requires further experimental research.

Overall, the present study, based on machine learning algorithms and various models, including a mediation model, discovered and elaborated the role and the underlying associations of musculoskeletal conditions in the development of depression in the elderly population and further fostered a novel IMCPR approach for depressive risk stratification. We found that arthritis does not convey a risk for negative psychological health outcomes unless it restricts individual activities (usual, social and leisure, or sleep). Furthermore, osteoarthritis, not rheumatoid arthritis, contributes the most to depressive symptoms. Among the musculoskeletal disorders, arthritis with activity restriction was the leading contributor to depression. Bone mineral density was negatively associated with depressive symptoms, and the depressive symptoms varied with the body part with the most severe MCP. The most important finding is that the IMCPR classification approach that we developed could specifically identify older people who require psychological intervention, thus optimizing social–medical resource redistribution and refining precision medicine by reducing the related global medical consumption due to musculoskeletal diseases.

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.

Ethics Statement

The studies involving human participants were reviewed and approved by Trinity College Dublin Research Ethics Committee. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

WJ: conceptualization and visualization. WJ and ZL: methodology and writing—original draft preparation. WJ, ZL, and YZ: data curation. WJ, ZL, and ZC: formal analysis. MG: writing—review and editing, project administration, and funding acquisition. All authors contributed to the article and approved the submitted version.

Funding

The TILDA study was supported by the Atlantic Philanthropies, Irish Life Plc, and the Irish Government. Additional funding for this work was provided by the Health Research Board of Ireland (grant number HRA_PHS/2012/30). Accessed via the Irish Social Science Data Archive—www.ucd.ie/issda. Additionally, this study was supported by grants from the National Natural Science Foundation of China (81371339 and 81870980). The funders were not involved in the study design, methods, subject recruitment, data collection, analysis, interpretation of data, writing of the paper, or submission for publication.

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.

Acknowledgments

We thank all individuals who took part in this study and the TILDA team who were responsible for the planning and administering of the study and making the datasets of TILDA available on their website (tilda.tcd.ie).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2021.697649/full#supplementary-material

Abbreviations

IMC, Individual musculoskeletal conditions; IMCPR, IMC-induced negative psychological risk; ARA, arthritis-induced restrictions of activities; YLDs, years lived with disability; MCP, musculoskeletal chronic pain; TILDA, The Irish Longitudinal Study on Aging; CES-D, Center for Epidemiological Studies Depression Scale; LLD, late-life depression; TG, triglyceride; LDL, low-density lipoprotein cholesterol; HDL, high-density lipoprotein cholesterol; CHOL, total cholesterol; BMI, body mass index; OR, odds ratios; CI, confidence intervals; P_H−L, p-value of Hosmer–Lemeshow test of goodness of fit; P_L−R, p-value of likelihood ratio test; M1:M2, likelihood ratio test conducted between model 1 and model 2; M2:M3, likelihood ratio test conducted between model 2 and model 3; ACME, average causal mediation effects; ADE, average direct effects; Prop. mediated, proportion mediated; χ², chi-square; n_obs, number of observation; û, group mean; $\widehat{{\epsilon }^2}$, effect size of ε²; p, p-value.

References

1. Oeppen J, Vaupel JW. Demography. Broken limits to life expectancy. Science. (2002) 296:1029–31. doi: 10.1126/science.1069675

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Partridge L, Deelen J, Slagboom PE. Facing up to the global challenges of ageing. Nature. (2018) 561:45–56. doi: 10.1038/s41586-018-0457-8

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Ferrari AJ, Charlson FJ, Norman RE, Patten SB, Freedman G, Murray CJ, et al. Burden of depressive disorders by country, sex, age, and year: findings from the global burden of disease study 2010. PLoS Med. (2013) 10:e1001547. doi: 10.1371/journal.pmed.1001547

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Vos T, Abajobir AA, Abate KH, Abbafati C, Abbas KM, Abd-Allah F, et al. Global, regional, and national incidence, prevalence, and years lived with disability for 328 diseases and injuries for 195 countries, 1990–2016: a systematic analysis for the global burden of disease study 2016. Lancet. (2017) 390:1211–59. doi: 10.1016/S0140-6736(17)32154-2

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Smith K. Mental health: a world of depression. Nature. (2014) 515:181. doi: 10.1038/515180a

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Briggs R, McCarroll K, O'Halloran A, Healy M, Kenny RA, Laird E. Vitamin D deficiency is associated with an increased likelihood of incident depression in community-dwelling older adults. J Am Med Dir Assoc. (2019) 20:517–23. doi: 10.1016/j.jamda.2018.10.006

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Gururajan A, Reif A, Cryan JF, Slattery DA. The future of rodent models in depression research. Nat Rev Neurosci. (2019) 20:686–701. doi: 10.1038/s41583-019-0221-6

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Malhi GS, Mann JJ. Depression. Lancet. (2018) 392:2299–312. doi: 10.1016/S0140-6736(18)31948-2

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Nock MK, Hwang I, Sampson NA, Kessler RC. Mental disorders, comorbidity and suicidal behavior: results from the national comorbidity survey replication. Mol Psychiatry. (2010) 15:868–76. doi: 10.1038/mp.2009.29

PubMed Abstract | CrossRef Full Text

10. World Health Organization. Depression. Geneva: WHO (2019).

Google Scholar

11. Steel N, Ford JA, Newton JN, Davis ACJ, Vos T, Naghavi M, et al. Changes in health in the countries of the UK and 150 english local authority areas 1990–2016: a systematic analysis for the global burden of disease study 2016. Lancet. (2018) 392:1647–61. doi: 10.1016/S0140-6736(18)32207-4

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Bradvik L, Mattisson C, Bogren M, Nettelbladt P. Long-term suicide risk of depression in the lundby cohort 1947-1997–severity and gender. Acta Psychiatr Scand. (2008) 117:185–91. doi: 10.1111/j.1600-0447.2007.01136.x

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Knapp M, Wong G. Economics and mental health: the current scenario. World Psychiatry. (2020) 19:3–14. doi: 10.1002/wps.20692

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Prince MJ, Wu F, Guo Y, Gutierrez Robledo LM, O'Donnell M, Sullivan R, et al. The burden of disease in older people and implications for health policy and practice. Lancet. (2015) 385:549–62. doi: 10.1016/S0140-6736(14)61347-7

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Rachner TD, Khosla S, Hofbauer LC. Osteoporosis: now and the future. Lancet. (2011) 377:1276–87. doi: 10.1016/S0140-6736(10)62349-5

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Center JR, Nguyen TV, Schneider D, Sambrook PN, Eisman JA. Mortality after all major types of osteoporotic fracture in men and women: an observational study. Lancet. (1999) 353:878–82. doi: 10.1016/S0140-6736(98)09075-8

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Loeser RF, Collins JA, Diekman BO. Ageing and the pathogenesis of osteoarthritis. Nat Rev Rheumatol. (2016) 12:412–20. doi: 10.1038/nrrheum.2016.65

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Tachmazidou I, Hatzikotoulas K, Southam L, Esparza-Gordillo J, Haberland V, Zheng J, et al. Identification of new therapeutic targets for osteoarthritis through genome-wide analyses of UK biobank data. Nat Genet. (2019) 51:230–6. doi: 10.1038/s41588-018-0327-1

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Hartvigsen J, Hancock MJ, Kongsted A, Louw Q, Ferreira ML, Genevay S, et al. What low back pain is and why we need to pay attention. Lancet. (2018) 391:2356–67. doi: 10.1016/S0140-6736(18)30480-X

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Buttgereit F, Smolen JS, Coogan AN, Cajochen C. Clocking in: chronobiology in rheumatoid arthritis. Nat Rev Rheumatol. (2015) 11:349–56. doi: 10.1038/nrrheum.2015.31

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Gerrits MM, van Oppen P, van Marwijk HW, van der Horst H, Penninx BW. The impact of chronic somatic diseases on the course of depressive and anxiety disorders. Psychother Psychosom. (2013) 82:64–6. doi: 10.1159/000338636

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Joaquim AF, Appenzeller S. Neuropsychiatric manifestations in rheumatoid arthritis. Autoimmun Rev. (2015) 14:1116–22. doi: 10.1016/j.autrev.2015.07.015

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Schwartz N, Temkin P, Jurado S, Lim BK, Heifets BD, Polepalli JS, et al. Chronic pain. Decreased motivation during chronic pain requires long-term depression in the nucleus accumbens. Science. (2014) 345:535–42. doi: 10.1126/science.1253994

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Imai K, Keele L, Yamamoto TJS. Identification, inference and sensitivity analysis for causal mediation effects. arXiv. (2010) 25:51–71. doi: 10.1214/10-STS321

CrossRef Full Text | Google Scholar

25. Liaw A, Wiener MJR. Classification and regression by randomForest. R News. (2002) 2:18–22. Available online at: https://cogns.northwestern.edu/cbmg/LiawAndWiener2002.pdf

Google Scholar

26. Tingley D, Yamamoto T, Hirose K, Keele L, Imai K. Mediation: R package for causal mediation analysis. J. Stat. Softw. (2014) 59:5. doi: 10.18637/jss.v059.i05

CrossRef Full Text | Google Scholar

27. von Elm E, Altman DG, Egger M, Pocock SJ, Gotzsche PC, Vandenbroucke JP, et al. [The strengthening the reporting of observational studies in epidemiology [STROBE] statement: guidelines for reporting observational studies]. Gac Sanit. (2008) 22:144–50. doi: 10.1157/13119325

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Whelan BJ, Savva GM. Design and methodology of the irish longitudinal study on ageing. J Am Geriatr Soc. (2013) 61 (Suppl. 2):S265–8. doi: 10.1111/jgs.12199

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Whelan BJJE, Review S. Ransam a random sample design for Ireland. Econ Soc Res Insti. (1979) 10:169–74.

Google Scholar

30. Barrett A, Burke H, Cronin H, Hickey A, Kamiya Y, Kenny RA, et al. Fifty Plus in Ireland 2011: First Results From the Irish Longitudinal Study on Ageing (TILDA). Dublin: Trinity College Dublin (2011). doi: 10.38018/TildaRe.2011-00. Available online at: https://repository.rcsi.com/articles/report/Fifty_plus_in_Ireland_2011_First_results_from_The_Irish_Longitudinal_Study_on_Ageing_TILDA_/10770389

CrossRef Full Text | Google Scholar

31. Nolan A, O'Regan C, Dooley C, Wallace D, Hever A, Cronin H, et al. The Over 50s in a Changing Ireland: Economic Circumstances, Health and Well-Being. Dublin: CARDI (2014). doi: 10.38018/TildaRe.2014-00

CrossRef Full Text

32. Radloff LS. The CES-D scale: a self-report depression scale for research in the general population. Appl Psychol Measure. (1977) 1:385–401. doi: 10.1177/014662167700100306

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Lyness JM, Noel TK, Cox C, King DA, Conwell Y, Caine EDJAoim. Screening for depression in elderly primary care patients: a comparison of the center for epidemiologic studies—depression scale and the geriatric depression scale. Arch Intern Med. (1997) 157:449–54. doi: 10.1001/archinte.157.4.449

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Sawyer Radloff L, Teri LJCG. Use of the center for epidemiological studies-depression scale with older adults. Clin Gerontol J Aging Ment Health. (1986) 5:119–36. doi: 10.1300/J018v05n01_06

PubMed Abstract | CrossRef Full Text

35. Farr JN, Khosla S. Skeletal changes through the lifespan–from growth to senescence. Nat Rev Endocrinol. (2015) 11:513–21. doi: 10.1038/nrendo.2015.89

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Kim SY, Lee JK, Oh DJ, Kong IG, Choi HG. Depression and incident hip fracture: a longitudinal follow-up study using a national sample cohort. Medicine. (2019) 98:e16268. doi: 10.1097/MD.0000000000016268

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Richard A, Rohrmann S, Lohse T, Eichholzer M. Is body weight dissatisfaction a predictor of depression independent of body mass index, sex and age? Results of a cross-sectional study. BMC Public Health. (2016) 16:863. doi: 10.1186/s12889-016-3497-8

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Wysokinski A, Strzelecki D, Kloszewska I. Levels of triglycerides, cholesterol, LDL, HDL and glucose in patients with schizophrenia, unipolar depression and bipolar disorder. Diabetes Metab Syndr. (2015) 9:168–76. doi: 10.1016/j.dsx.2015.04.004

PubMed Abstract | CrossRef Full Text

39. Yan XY, Huang SM, Huang CQ, Wu WH, Qin Y. Marital status and risk for late life depression: a meta-analysis of the published literature. J Int Med Res. (2011) 39:1142–54. doi: 10.1177/147323001103900402

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Subhash R, Jonah L, Solymos P. ResourceSelection: Resource Selection (Probability) Functions for Use-Availability Data. R package version 0.3–0 (2019).

41. Zeileis A, Hothorn T. Diagnostic Checking in Regression Relationships. Magdeburg: Otto-von-Guericke-Universität Magdeburg (2002). Available online at: http://pkg.cs.ovgu.de/LNF/i386/5.10/R/LNFr-lmtest/reloc/R-2.10/library/lmtest/doc/lmtest-intro.pdf

Google Scholar

42. Patil I, Powell C. ggstatsplot:“ggplot2” Based Plots With Statistical Details. CRAN. (2018). Retrieved from https://cranr-project.org/web/packages/ggstatsplot

Google Scholar

43. Cohn SL, Pearson AD, London WB, Monclair T, Ambros PF, Brodeur GM, et al. The international neuroblastoma risk group (INRG) classification system: an INRG task force report. J Clin Oncol. (2009) 27:289–97. doi: 10.1200/JCO.2008.16.6785

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Wang Y, Pan Y, Zhao X, Li H, Wang D, Johnston SC, et al. Clopidogrel with aspirin in acute minor stroke or transient ischemic attack (CHANCE) trial: one-year outcomes. Circulation. (2015) 132:40–6. doi: 10.1161/CIRCULATIONAHA.114.014791

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Vallerand IA, Patten SB, Barnabe C. Depression and the risk of rheumatoid arthritis. Curr Opin Rheumatol. (2019) 31:279–84. doi: 10.1097/BOR.0000000000000597

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Albayrak Gezer I, Balkarli A, Can B, Bagcaci S, Kucuksen S, Kucuk A. Pain, depression levels, fatigue, sleep quality, and quality of life in elderly patients with rheumatoid arthritis. Turk J Med Sci. (2017) 47:847–53. doi: 10.3906/sag-1603-147

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Cauley JA, Chalhoub D, Kassem AM, Fuleihan Gel H. Geographic and ethnic disparities in osteoporotic fractures. Nat Rev Endocrinol. (2014) 10:338–51. doi: 10.1038/nrendo.2014.51

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Richards JB, Rivadeneira F, Inouye M, Pastinen TM, Soranzo N, Wilson SG, et al. Bone mineral density, osteoporosis, and osteoporotic fractures: a genome-wide association study. Lancet. (2008) 371:1505–12. doi: 10.1016/S0140-6736(08)60599-1

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Weycker D, Edelsberg J, Barron R, Atwood M, Oster G, Crittenden DB, et al. Predictors of near-term fracture in osteoporotic women aged >/=65 years, based on data from the study of osteoporotic fractures. Osteoporos Int. (2017) 28:2565–71. doi: 10.1007/s00198-017-4103-3

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Guglielmo D, Hootman JM, Boring MA, Murphy LB, Theis KA, Croft JB, et al. Symptoms of anxiety and depression among adults with arthritis - United States, 2015-2017. MMWR Morb Mortal Wkly Rep. (2018) 67:1081–7. doi: 10.15585/mmwr.mm6739a2

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Sambamoorthi U, Shah D, Zhao X. Healthcare burden of depression in adults with arthritis. Expert Rev Pharmacoecon Outcomes Res. (2017) 17:53–65. doi: 10.1080/14737167.2017.1281744

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Sturgeon JA, Finan PH, Zautra AJ. Affective disturbance in rheumatoid arthritis: psychological and disease-related pathways. Nat Rev Rheumatol. (2016) 12:532–42. doi: 10.1038/nrrheum.2016.112

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Kwiatkowska B, Klak A, Raciborski F, Maslinska M. The prevalence of depression and insomnia symptoms among patients with rheumatoid arthritis and osteoarthritis in Poland: a case control study. Psychol Health Med. (2019) 24:333–43. doi: 10.1080/13548506.2018.1529325

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Nerurkar L, Siebert S, McInnes IB, Cavanagh J. Rheumatoid arthritis and depression: an inflammatory perspective. Lancet Psychiatry. (2019) 6:164–73. doi: 10.1016/S2215-0366(18)30255-4

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Hunter DJ, Bierma-Zeinstra S. Osteoarthritis. Lancet. (2019) 393:1745–59. doi: 10.1016/S0140-6736(19)30417-9

PubMed Abstract | CrossRef Full Text

56. Lee JE, Martire LM, Zarit SH, Rovine MJ. Activity restriction and depressive symptoms in older couples. J Aging Health. (2017) 29:1251–67. doi: 10.1177/0898264316660413

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Mausbach BT, Chattillion EA, Moore RC, Roepke SK, Depp CA, Roesch S. Activity restriction and depression in medical patients and their caregivers: a meta-analysis. Clin Psychol Rev. (2011) 31:900–8. doi: 10.1016/j.cpr.2011.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Leino-Arjas P, Rajaleid K, Mekuria G, Nummi T, Virtanen P, Hammarstrom A. Trajectories of musculoskeletal pain from adolescence to middle age: the role of early depressive symptoms, a 27-year follow-up of the Northern Swedish cohort. Pain. (2018) 159:67–74. doi: 10.1097/j.pain.0000000000001065

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Williams LJ, Pasco JA, Jacka FN, Henry MJ, Dodd S, Berk M. Depression and bone metabolism. A review. Psychother Psychosom. (2009) 78:16–25. doi: 10.1159/000162297

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Ali K, El Said SM, Adly NN, Abdul-Rahman SA. The relation between calcaneus stiffness index as a measure of bone density and body mass index in an Egyptian cohort. J Multidiscip Healthc. (2019) 12:1085–90. doi: 10.2147/JMDH.S230730

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Duman RS. Sex-specific disease-associated modules for depression. Nat Med. (2017) 23:1015–7. doi: 10.1038/nm.4391

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Labonte B, Engmann O, Purushothaman I, Menard C, Wang J, Tan C, et al. Sex-specific transcriptional signatures in human depression. Nat Med. (2017) 23:1102–11. doi: 10.1038/nm.4386

PubMed Abstract | CrossRef Full Text | Google Scholar