Skip to main content
  • Research article
  • Open access
  • Published:

Variability in the prevalence of depression among adults with chronic pain: UK Biobank analysis through clinical prediction models

Abstract

Background

The prevalence of depression among people with chronic pain remains unclear due to the heterogeneity of study samples and definitions of depression. We aimed to identify sources of variation in the prevalence of depression among people with chronic pain and generate clinical prediction models to estimate the probability of depression among individuals with chronic pain.

Methods

Participants were from the UK Biobank. The primary outcome was a “lifetime” history of depression. The model’s performance was evaluated using discrimination (optimism-corrected C statistic) and calibration (calibration plot).

Results

Analyses included 24,405 patients with chronic pain (mean age 64.1 years). Among participants with chronic widespread pain, the prevalence of having a “lifetime” history of depression was 45.7% and varied (25.0–66.7%) depending on patient characteristics. The final clinical prediction model (optimism-corrected C statistic: 0.66; good calibration on the calibration plot) included age, BMI, smoking status, physical activity, socioeconomic status, gender, history of asthma, history of heart failure, and history of peripheral artery disease. Among participants with chronic regional pain, the prevalence of having a “lifetime” history of depression was 30.2% and varied (21.4–70.6%) depending on patient characteristics. The final clinical prediction model (optimism-corrected C statistic: 0.65; good calibration on the calibration plot) included age, gender, nature of pain, smoking status, regular opioid use, history of asthma, pain location that bothers you most, and BMI.

Conclusions

There was substantial variability in the prevalence of depression among patients with chronic pain. Clinically relevant factors were selected to develop prediction models. Clinicians can use these models to assess patients’ treatment needs. These predictors are convenient to collect during daily practice, making it easy for busy clinicians to use them.

Peer Review reports

Background

Chronic pain is one of the leading causes of disability, affecting more than 30% of people worldwide [1]. Depression is also a leading cause of disability, affecting approximately 5% of adults worldwide [2, 3]. It is generally understood that chronic pain and depression are commonly co-morbid disorders [1, 4]. Indeed, research suggests that chronic pain increases the risk of depression, and depression increases the risk of chronic pain [5, 6]. However, the prevalence of depression among people living with chronic pain remains unclear [1, 7]. Previous studies have reported prevalence estimates ranging from about 15% to 85% [8,9,10]. There are several possible reasons for such variation in the prevalence of depression among people with chronic pain reported across studies. Firstly, measures of depression and definitions of depression vary considerably across studies. For example, some studies measure current depression, while others measure lifetime depression [10]. Secondly, the extent of pain varied (e.g., regional vs widespread pain). Thirdly, the demographic and health characteristics of the populations sampled varied. For example, people with chronic pain who are female, have additional chronic health conditions, or have a lower socioeconomic status are thought to be at higher risk of depression [11].

A clinical prediction model could calculate the risk of a particular endpoint for individual patients by combining multiple predictors, which could be a useful way to accurately estimate the probability that patients with chronic pain suffer from depression based on their individual characteristics [12]. Recently published methodological papers have provided a framework for the development of valid clinical prediction models [13,14,15].

The selection of the appropriate dataset is important for the development of a valid clinical prediction model. Among potentially suitable datasets, we selected the UK Biobank dataset for the following reasons. Firstly, at its baseline visit, the UK Biobank recruited about 0.5 million participants across the UK, which provided a large sample size to start a study. Secondly, the “experience of pain” questionnaire (2019–2020) provides a comprehensive assessment of chronic pain, including regional or widespread pain, neuropathic or non-neuropathic pain, and pain location that bothers you most. Thirdly, the validity of the measurement of depression in the “online mental health self-assessment” questionnaire (2016–2017) is supported by a dual approach that includes both secondary care record linkage (i.e., diagnosis by a professional) and self-reporting of symptoms [16]. Using this dataset, we aimed to develop and internally validate clinical prediction models of depression among individuals with chronic pain.

Methods

Study sample

This study used data from the UK Biobank. UK Biobank is a large-scale biomedical database, which recruited approximately 500,000 people in the UK at its initial enrollment (from 13 March 2006 to 1 October 2010). Part of these participants received follow-up surveys. For example, about 157,000 participants received the “online mental health self-assessment” questionnaire from 13 July 2016 to 27 July 2017, and about 167,000 participants received the “experience of pain” questionnaire from 9 January 2019 to 18 April 2020 [17]. More details about the UK Biobank can be found in the registry online protocol: http://www.ukbiobank.ac.uk. The North West Multi-centre Ethics Committee granted ethical approval to access data from the UK Biobank, and all participants provided written informed consent.

To define chronic pain, we selected the “experience of pain” questionnaire (2019–2020) rather than the baseline visit (2006–2010) for the following reasons. Firstly, the number of pain types in the “experience of pain” questionnaire was much higher than at the baseline visit (i.e., 15 in the “experience of pain” questionnaire compared with 8 in the baseline visit). Secondly, the “experience of pain” questionnaire collected a number of additional pain-related variables (e.g., neuropathic pain or not, and the pain area that bothers you the most). To match the measurement time of chronic pain and depression, the analysis sample was restricted to participants who reported having pain for more than 5 years in the “experience of pain” questionnaire (2019–2020) and completed the “online mental health self-assessment” questionnaire (2016–2017). Based on the International Classification of Diseases 11th Revision definitions for chronic pain and the data availability of UK Biobank, chronic pain was classified as widespread pain (through the question “have you experienced pain or discomfort all over the body?”) and regional pain (i.e., leg pain, chest pain, feet pain, hand pain, arm pain, knee pain, hip pain, stomach or abdominal pain, back pain, neck or shoulder pain, facial pain, and headache) [18].

Although previous literature suggested that multisite pain is strongly related to mood disorders and played an important role in the development of chronic pain, UK Biobank has created a new question, “the pain area that bothers you the most,” in consideration of the fact that many people have multiple pains [19, 20]. Therefore, we included the pain area that bothers you the most as one of the predictors. We also collected the nature of pain (neuropathic and non-neuropathic pain) as one pain-related characteristic [21]. Details for defining pain can be found in Supplementary A.

Outcomes

We followed the framework that the UK Biobank team proposed to define the depression [16]. The primary outcome was a “lifetime” history of depression rather than present depression, because many mental disorders (e.g., depression) can fluctuate. By including those with a “lifetime” history, we are more likely to more comprehensively capture those with the condition. The dual approach was used to define a “lifetime” history of depression, which included both secondary care record linkage (i.e., diagnosed by a professional) and self-report of symptoms through the Composite International Diagnostic Interview-Short Form (CIDI-SF), depression module, lifetime version. CIDI-SF is a simplified version of its full version CIDI [22] which is a fully structured diagnostic interview, and one previous validation study showed CIDI-SF had comparable accuracy for diagnosing major depressive episodes when compared to CIDI [23]. Two reasons justified the choice of the dual approach: firstly, traditional full-version diagnostic interview is too expensive to be implemented in a cohort with a large sample size (e.g., UK Biobank). Secondly, secondary care record linkage can fail to identify patients with less severe illnesses as these patients are less likely to seek help from the professional compared with patients with more severe illnesses [24] Through this dual approach, all participants were classified as having no “lifetime” history of depression, having a “lifetime” history of subthreshold depressive symptoms, and having a “lifetime” history of depression.

Following the framework that the UK Biobank team proposed, the secondary outcome was present depression [16]. It is worth noting that the UK Biobank team identified present depression among participants with a history of depression, but did not provide clear justification for this approach. Readers should be aware of this point when interpreting the results. Present depression was defined through the Patient Health Questionnaire 9-question version (PHQ-9). PHQ-9 is a validated tool that included nine short screening questionnaires and is widely used in screening for depression [25].

The detailed algorithms and the corresponding R code to define the above outcome were provided by the official group, as available at https://data.mendeley.com/datasets/kv677c2th4/3.

Covariates

Previous systematic reviews have identified factors that are known to increase risk of depression [11, 26, 27]. Based on these findings and data availability in the UK Biobank and in daily practice, we consider the following variables as covariates: demographic characteristics (age, gender, ethnicity, and Townsend deprivation score which reflected socioeconomic status), body mass index (BMI), lifestyle behaviors (smoking status, alcohol consumption, and physical activity), comorbidities as identified in the recent international consensus on the definition of multimorbidity for research purposes (i.e., stroke, coronary artery disease, heart failure, peripheral artery disease, diabetes, Addison’s disease, cystic fibrosis, chronic obstructive pulmonary disease, asthma, Parkinson’s disease, epilepsy, multiple sclerosis, paralysis, solid organ cancers, hematological cancers, metastatic cancers, dementia, schizophrenia, connective tissue disease, chronic liver disease, inflammatory bowel disease, chronic kidney disease, end-stage kidney disease, and HIV/AIDS) [28], and regular opioid use. For participants with chronic regional pain, nature of pain, and pain location that bothers you most were also added. Definition details could be found in Supplementary B. Other pain severity-related variables were not included as predictors due to the concerns with the potential measurement bias. For example, pain intensity was measured through the question “Thinking about the last 24 hours, how would you rate your pain on a 0-10 scale, where 0 is ‘no pain’ and 10 is ‘pain’ as bad as it could be,” which may not align with the timeline of when patients completed the mental health questionnaire.

Statistical analysis

Baseline characteristics for participants with chronic pain were shown by depression status. Overall and subgroup prevalence of having: (1) a “lifetime” history of depression among participants with chronic widespread pain; (2) a “lifetime” history of depression among participants with chronic regional pain; (3) present depression among participants with chronic widespread pain; (4) present depression among participants with chronic regional pain were provided. Subgroup analyses were performed based on the “one covariate at a time” principle by each of the variables mentioned in the covariates section. Wald statistic was used to assess whether the prevalence differed by each covariate [29].

Prediction models (through logistic regression) to estimate the probability of depression for individuals with chronic pain were developed. The choice of logistic regression was based on its ease of understanding and communication, as well as its ability to handle binary outcomes [30]. To ensure precise predictions and prevent overfitting, the maximum number of candidate predictor parameters was estimated based on the criteria proposed (details in Supplementary C) by Riley et al. [31]. To minimize the influence of sparse data from binary predictors, we excluded predictors if the number of events in one level of the predictor was less than 10. If the remaining predictors were still more than the estimated maximum number, we excluded predictors with an insignificant Wald statistic. Considering most covariates have a small quantity of missing data (details in Table 1), a single imputation through the transcan function (i.e., a nonlinear additive transformation and imputation function) was used [29].

Table 1 Baseline characteristics for participants with chronic pain stratified by depression status

The modeling strategy we used was adapted from Harrell’s Regression Modeling Strategies (detailed in Fig. 1) [29] The full model, including all pre-specified predictors without variable selection, was considered the gold standard. However, clinicians may have insufficient resources (e.g., time) to collect all these predictors. Thus, the simplified model may be needed in daily practice. One significant benefit of Harrell’s simplified model is that it offers varying degrees of parsimony to clinicians based on their specific needs. This is achieved by estimating the contribution of each predictor. In our study, we provide two examples. Firstly, the simplified model (reported as equations and nomograms) has at least 95% of the performance compared with the full model. Secondly, we assume that the clinician only wants to collect the three most important predictors.

Fig. 1
figure 1

Summary of the modeling strategy

Model performance was assessed by the discrimination (through optimism-corrected C statistic) and calibration (through calibration plot) [12]. Optimism is defined as a bias due to overfitting. The bootstrap method is a class of resampling methods that samples a sub-dataset from the original one with replacement. The estimate of the optimism equals the C statistic from the original sample minus the C statistic from the bootstrap sample. In our study, this process was repeated 1000 times to get an average optimism. The final reported optimism-corrected C statistic equals the C statistic from the original sample minus the average optimism [29]. In addition, the C-statistic with the 95% confidence interval using 10-fold cross-validation was provided. We checked whether two continuous variables (age and BMI) should be modeled through splines and the results showed that they can be analyzed through the original form. Based on clinical knowledge and other literature, we assessed the potential interaction between age and ethnicity and the results showed that we did not need to include this interaction term in the model [32]. Details for modeling could be found in Supplementary D.

For chronic regional pain, although one prediction model may not work well for different categories, we did not develop a clinical prediction model for each category as the sample size may be insufficient. To explore the robustness of the prediction model for the overall chronic regional pain, we performed an additional analysis by evaluating its model performance for each category of chronic regional pain.

We reported this study based upon the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) and Transparent reporting of a multivariable prediction model for individual prognosis or diagnosis (TRIPOD) statement [33, 34]. All statistical analyses were performed in R, version 4.2.2 (R Group for Statistical Computing).

Results

Of the UK Biobank participants, 24,405 participants with chronic pain were included: 912 (3.7%) had present depression, 7952 (32.6%) had a “lifetime” history of depression, 5317 (21.8%) had a “lifetime” history of subthreshold depressive symptoms, and 11,137 (45.6%) had no “lifetime” history of depression. Figure 2 shows the selection process. Table 1 reports the participants’ baseline characteristics. Among included participants, 9228 (37.8%) were men, the mean (SD) age was 64.1 (7.5) years, and 23,706 (97.1%) were white. Univariate associations of the covariates with depression outcomes could be found in Supplementary E.

Fig. 2
figure 2

Flowchart of study participants

Primary outcome

Among participants with chronic widespread pain, the prevalence of having a “lifetime” history of depression was 45.7% (1716/3757) (Table 2). Subgroup analyses revealed that the prevalence ranged from 25.0 to 66.7% (Table 2). 26 predictors were included in the initial full prediction model (Supplementary F). The final simplified model (Supplementary G) with nine predictors (age, BMI, smoking status, physical activity, Townsend deprivation score, gender, history of asthma, history of heart failure, and history of peripheral artery disease) was built with its equation in Supplementary H and the nomogram in Fig. 3. The prediction model showed moderate discrimination (optimism-corrected C statistic was 0.66; C statistic from the 10-fold cross-validation: 0.67, 95% confidence interval [CI] 0.65 to 0.69) and good calibration (on the calibration plot) (Supplementary I). Age (as age increases by one year, the odds of having a “lifetime” history of depression decreases: odds ratio [OR] 0.94, 95% CI 0.93 to 0.95), gender (compared to females, males were less likely to have a “lifetime” history of depression: OR 0.56, 95% CI 0.47 to 0.65), and BMI (as the value of BMI increase by one, the odds of having a “lifetime” history of depression also increases: OR 1.02, 95% CI 1.01 to 1.03) were the three most important predictors.

Table 2 Overall and subgroup prevalence of depression among participants with chronic pain
Fig. 3
figure 3

Nomogram for estimating the probability of having a “lifetime” history of depression for individuals with chronic widespread pain. Gender: male—1 and female—0. History of one comorbidity: yes—1 and no—0. Instructions for the use of the nomogram: (1) locate the answer for each predictor, (2) draw a straight line upward to the point axis and record the score, (3) calculate the total score for all predictors and locate the score in the total points axis, (4) draw a straight line downward to the probability axis to estimate the individual’s probability of having a “lifetime” history of depression

Among participants with chronic regional pain, the prevalence of having a “lifetime” history of depression was 30.2% (6235/20,648) (Table 2). Subgroup analyses revealed that the prevalence ranged from 21.4 to 70.6% (Table 2). Thirty predictors were included in the initial full prediction model (Supplementary F). The final simplified model (Supplementary G) with eight predictors (age, gender, nature of pain, smoking status, regular opioid use, history of asthma, pain location that bothers you most, and BMI) was built with its equation in Supplementary H and the nomogram in Fig. 4. The prediction model showed moderate discrimination (optimism-corrected C statistic was 0.65; C statistic from the 10-fold cross-validation: 0.66, 95% CI 0.65 to 0.66) and good calibration (on the calibration plot) (Supplementary I). Age (as age increases by one year, the odds of having a “lifetime” history of depression decreases: OR 0.96, 95% CI 0.96 to 0.96), gender (compared to females, males were less likely to have a “lifetime” history of depression: OR 0.53, 95% CI 0.50 to 0.57), and nature of pain (compared with patients with non-neuropathic pain, patients with neuropathic pain were more likely to have a “lifetime” history of depression: OR 1.47, 95% CI 1.36 to 1.58) were the three most important predictors.

Fig. 4
figure 4

Nomogram for estimating the probability of having a “lifetime” history of depression for individuals with chronic regional pain. Gender: male—1 and female—0. Nature of pain: neuropathic pain—1 and non-neuropathic pain—0. Regular opioid use: yes—1 and no—0. History of asthma: yes—1 and no—0. Pain location that bothers you most: arm pain—a, back pain—b, chest pain—c, facial pain—d, feet pain—e, hand pain—f, headache—g, hip pain—h, knee pain—i, leg pain—j, neck or shoulder pain—k, and stomach or abdominal pain—l. Instructions for the use of the nomogram: (1) locate the answer for each predictor, (2) draw a straight line upward to the point axis and record the score, (3) calculate the total score for all predictors and locate the score in the total points axis, (4) draw a straight line downward to the probability axis to estimate the individual’s probability of having a “lifetime” history of depression

Secondary outcome

Among participants with chronic widespread pain, the prevalence of having present depression was 10.5% (396/3757) (Table 2). Subgroup analyses revealed that the prevalence ranged from 4.5 to 33.3% (Table 2). In total, 13 predictors were included in the initial full prediction model (Supplementary F). The final simplified model (Supplementary G) with seven predictors (age, BMI, smoking status, physical activity, Townsend deprivation score, history of peripheral artery disease, and history of chronic kidney disease) was built with its equation in Supplementary H and the nomogram in Supplementary J. The prediction model showed moderate discrimination (optimism-corrected C statistic was 0.75; C statistic from the 10-fold cross-validation: 0.76, 95% CI 0.74 to 0.79) and good calibration (on the calibration plot) (Supplementary I). Age (as age increases by one year, the odds of having present depression decreases: OR 0.91, 95% CI 0.90 to 0.93), BMI (as the value of BMI increases by one, the odds of having present depression also increases: OR 1.04, 95% CI 1.02 to 1.06), and smoking status (compared to current smokers, both former [OR 0.62, 95% CI 0.44 to 0.86] and never [OR 0.47, 95% CI 0.34 to 0.65] smokers were less likely to have present depression) were the three most important predictors.

Among participants with chronic regional pain, the prevalence of having present depression was 2.5% (516/20,648) (Table 2). Subgroup analyses revealed that the prevalence ranged from 1.4 to 26.7% (Table 2). In total, 17 predictors were included in the initial full prediction model (Supplementary F). The final simplified model (Supplementary G) with 10 predictors (age, BMI, nature of pain, pain location that bothers you most, Townsend deprivation score, regular opioid use, physical activity, smoking status, history of diabetes, and history of chronic obstructive pulmonary disease) was built with its equation in Supplementary H and the nomogram in Supplementary J. The prediction model showed moderate discrimination (optimism-corrected C statistic was 0.74; C statistic from the 10-fold cross-validation: 0.75, 95% CI 0.73 to 0.77) and good calibration (on the calibration plot) (Supplementary I). Age (as age increases by one year, the odds of having present depression decrease: OR 0.93, 95% CI 0.92 to 0.94), BMI (as the value of BMI increases by one, the odds of having present depression also increases: OR 1.06, 95% CI 1.04 to 1.07), and nature of pain (compared with patients with non-neuropathic pain, patients with neuropathic pain were more likely to have present depression: OR 1.71, 95% CI 1.40 to 2.10) were the three most important predictors.

Additional analyses

For the primary outcome (i.e., a “lifetime” history of depression), the results showed that the model developed for the overall chronic regional pain also worked well for all categories (optimism-corrected C statistics: 0.62 to 0.67) of chronic regional pain except for stomach pain (optimism-corrected C statistic: 0.59). For the secondary outcome (i.e., the present depression), the results showed that the model developed for the overall chronic regional pain worked well for all categories (optimism-corrected C statistics: 0.69 to 0.78) of chronic regional pain except for chest pain (optimism-corrected C statistics: 0.64), feet pain (optimism-corrected C statistics: 0.65), hand pain (optimism-corrected C statistics: 0.66), and headache (optimism-corrected C statistics: 0.60).

Discussion

Key results

We found that there was substantial variability in the prevalence of having a “lifetime” history of depression among patients with chronic pain. Among participants with chronic widespread pain, the prevalence of having a “lifetime” history of depression was 45.7%; subgroup analyses indicated that the prevalence ranged from 25.0 to 66.7%.

This study developed and evaluated clinical prediction models to estimate the probability of having a “lifetime” history of depression among patients with chronic pain. Among participants with chronic widespread pain, the final clinical prediction model consisted of nine predictors, including age, BMI, smoking status, physical activity, Townsend deprivation score, gender, history of asthma, history of heart failure, and history of peripheral artery disease. Among participants with chronic regional pain, the final clinical prediction model consisted of eight predictors, including age, gender, nature of pain, smoking status, regular opioid use, history of asthma, pain location that bothers you most, and BMI.

Comparison with previous studies

Using the terms “chronic pain,” “depression,” and “UK Biobank” in PubMed (from the inception to March 1, 2024), we found 13 studies including chronic pain and depression through the analysis of UK Biobank [19, 20, 35,36,37,38,39,40,41,42,43,44,45]. Of the 13 studies, five focused on genetic information [35, 36, 39, 41, 45], five were association analyses [19, 37, 38, 42, 44], one examined the role of coffee in the association between chronic pain and depression [43], one was a clinical prediction model for the development and spread of chronic pain [20], and one assessed risk factors for facial pain [40]. We also extended the search to clinical prediction models based on other datasets and found no other relevant studies. Therefore, this is the first study to develop prediction models that estimate individuals’ probability of experiencing depression among participants with chronic pain. Our models reported through TRIPOD guidelines, showed moderate discrimination and good calibration.

Although our study could not answer the question of bidirectional causality between chronic pain and depression, readers should bear in mind the complex interplay between chronic pain and depression when interpreting the results. Previous studies have reported the role of depression in the chronicity of pain, especially the nociplastic type of pain (fibromyalgia), and the role of chronic pain in the development of depression [20, 46]. Repeated measurements of both pain and depression could facilitate a deeper exploration into whether chronic pain predisposes patients to depression or vice versa [47].

Limitations

Several limitations should be noted. Firstly, the difference in the measurement time for chronic pain and depression status might bring bias. Although we restricted the analysis sample to those whose pain duration was more than 5 years, we could not totally exclude the influence of recall bias [48]. We also did not find formal analysis to assess the reliability of this retrospective way to define chronic pain, further studies should be performed to assess the accuracy of the estimate. Secondly, genetic information (e.g., polygenic risk scores) may add additional value in predicting depression among individuals with chronic pain, as previous studies have found the genetic relationship between pain and depression [49, 50]. However, the data we applied for this project did not include genetic data. This should be investigated in further studies. Thirdly, although we did not include the number of pain sites as one of the predictors considering the measurement issue in the relevant questionnaire in the UK Biobank, this variable may provide useful information, which should be collected in future studies with more accurate pain questionnaires. Fourthly, although an external validation would be beneficial, a suitable dataset for comparison with the UK Biobank was not found. Further validation studies that prospectively collect data with the comprehensive assessment of chronic pain and depression status are therefore still needed. Fifthly, participants in this study were from the UK, between the ages of 47 to 80, meaning results may not be generalizable to other countries or age groups. Finally, non-white patients were grouped into one category (i.e., the ethnicity in the analysis was treated as a binary variable: white and non-white) in this study to facilitate analysis. However, this way might mask the differences among these non-white patients, which should be explored in future studies.

Implications of clinical practice and future research

Results from this study can support clinicians in deciding upon treatment priorities for the patient. Importantly, the predictors included are easily collected by clinicians. To further enhance the model, future researchers should focus on improving the quality of the measurement instruments and look to objective assessment when possible. They should also consider other potentially important predictors to improve the predictive accuracy of the model, such as genetic information. Finally, external validation should take place. As Riley et al. mentioned in their new methodological paper, researchers should focus on the target population and setting in which the model is planned to be implemented, especially when the intended population or setting is different from the one in which the model was developed (e.g., UK Biobank) [14].

Conclusions

There was substantial variability in the prevalence of depression among patients with chronic pain. Clinically relevant factors were selected to develop prediction models. Clinicians can use these models to assess patients’ treatment needs. These predictors are convenient to collect during daily practice, making it easy for busy clinicians to use them.

Availability of data and materials

Access to data can be requested via application to the UK Biobank (https://www.ukbiobank.ac.uk/). The study was conducted under UK Biobank project number 80584.

Abbreviations

BMI:

Body mass index

CIDI-SF:

Composite International Diagnostic Interview-Short Form

PHQ-9:

Patient Health Questionnaire 9-question version

SD:

Standard deviation

STROBE:

Strengthening the Reporting of Observational Studies in Epidemiology

TRIPOD:

Transparent reporting of a multivariable prediction model for individual prognosis or diagnosis

References

  1. Cohen SP, Vase L, Hooten WM. Chronic pain: an update on burden, best practices, and new advances. Lancet. 2021;397(10289):2082–97.

    Article  PubMed  Google Scholar 

  2. Global burden of 369 diseases and injuries in 204 countries and territories, 1990-2019: a systematic analysis for the Global Burden of Disease Study 2019. Lancet. 2020;396(10258):1204-22.

  3. Organization WH. Depression. 2023.  https://www.who.int/health-topics/depression#tab=tab_1.

  4. Patten SB, Williams JV, Wang J. Mental disorders in a population sample with musculoskeletal disorders. BMC Musculoskelet Disord. 2006;7:37.

    Article  PubMed  PubMed Central  Google Scholar 

  5. McCarthy LH, Bigal ME, Katz M, Derby C, Lipton RB. Chronic pain and obesity in elderly people: results from the Einstein aging study. J Am Geriatr Soc. 2009;57(1):115–9.

    Article  PubMed  Google Scholar 

  6. Ohayon MM, Stingl JC. Prevalence and comorbidity of chronic pain in the German general population. J Psychiatr Res. 2012;46(4):444–50.

    Article  PubMed  Google Scholar 

  7. Bair MJ, Robinson RL, Katon W, Kroenke K. Depression and pain comorbidity: a literature review. Arch Intern Med. 2003;163(20):2433–45.

    Article  PubMed  Google Scholar 

  8. Cherif F, Zouari HG, Cherif W, Hadded M, Cheour M, Damak R. Depression prevalence in neuropathic pain and its impact on the quality of life. Pain Res Manag. 2020;2020:7408508.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Elbinoune I, Amine B, Shyen S, Gueddari S, Abouqal R, Hajjaj-Hassouni N. Chronic neck pain and anxiety-depression: prevalence and associated risk factors. Pan Afr Med J. 2016;24:89.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Kleykamp BA, Ferguson MC, McNicol E, Bixho I, Arnold LM, Edwards RR, et al. The prevalence of psychiatric and chronic pain comorbidities in fibromyalgia: an ACTTION systematic review. Semin Arthritis Rheum. 2021;51(1):166–74.

    Article  PubMed  Google Scholar 

  11. McCarron RM, Shapiro B, Rawles J, Luo J. Depression. Ann Intern Med. 2021;174(5):Itc65-itc80.

  12. Steyerberg EW, Moons KG, van der Windt DA, Hayden JA, Perel P, Schroter S, et al. Prognosis Research Strategy (PROGRESS) 3: prognostic model research. PLoS Med. 2013;10(2):e1001381.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Collins GS, Dhiman P, Ma J, Schlussel MM, Archer L, Van Calster B, et al. Evaluation of clinical prediction models (part 1): from development to external validation. BMJ. 2024;384:e074819.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Riley RD, Archer L, Snell KIE, Ensor J, Dhiman P, Martin GP, et al. Evaluation of clinical prediction models (part 2): how to undertake an external validation study. BMJ. 2024;384:e074820.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Riley RD, Snell KIE, Archer L, Ensor J, Debray TPA, van Calster B, et al. Evaluation of clinical prediction models (part 3): calculating the sample size required for an external validation study. BMJ. 2024;384:e074821.

    Article  PubMed  Google Scholar 

  16. Davis KAS, Coleman JRI, Adams M, Allen N, Breen G, Cullen B, et al. Mental health in UK Biobank - development, implementation and results from an online questionnaire completed by 157 366 participants: a reanalysis. BJPsych Open. 2020;6(2):e18.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Sudlow C, Gallacher J, Allen N, Beral V, Burton P, Danesh J, et al. UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS Med. 2015;12(3):e1001779.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Treede RD, Rief W, Barke A, Aziz Q, Bennett MI, Benoliel R, et al. Chronic pain as a symptom or a disease: the IASP Classification of Chronic Pain for the International Classification of Diseases (ICD-11). Pain. 2019;160(1):19–27.

    Article  PubMed  Google Scholar 

  19. Nicholl BI, Mackay D, Cullen B, Martin DJ, Ul-Haq Z, Mair FS, et al. Chronic multisite pain in major depression and bipolar disorder: cross-sectional study of 149,611 participants in UK Biobank. BMC Psychiatry. 2014;14:350.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Tanguay-Sabourin C, Fillingim M, Guglietti GV, Zare A, Parisien M, Norman J, et al. A prognostic risk score for development and spread of chronic pain. Nat Med. 2023;29(7):1821–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Bouhassira D, Attal N, Alchaar H, Boureau F, Brochet B, Bruxelle J, et al. Comparison of pain syndromes associated with nervous or somatic lesions and development of a new neuropathic pain diagnostic questionnaire (DN4). Pain. 2005;114(1–2):29–36.

    Article  PubMed  Google Scholar 

  22. Robins LN, Wing J, Wittchen HU, Helzer JE, Babor TF, Burke J, et al. The Composite International Diagnostic Interview. An epidemiologic Instrument suitable for use in conjunction with different diagnostic systems and in different cultures. Arch Gen Psychiatry. 1988;45(12):1069–77.

    Article  CAS  PubMed  Google Scholar 

  23. Kessler RC, Andrews G, Mroczek D, Ustun B, Wittchen H-U. The World Health Organization Composite International Diagnostic Interview short-form (CIDI-SF). Int J Methods Psychiatr Res. 1998;7(4):171–85.

    Article  Google Scholar 

  24. Newson RS, Karlsson H, Tiemeier H. Epidemiological fallacies of modern psychiatric research. Nord J Psychiatry. 2011;65(4):226–37.

    Article  PubMed  Google Scholar 

  25. Manea L, Gilbody S, McMillan D. Optimal cut-off score for diagnosing depression with the Patient Health Questionnaire (PHQ-9): a meta-analysis. CMAJ. 2012;184(3):E191-6.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Hammen C. Risk Factors for Depression: An Autobiographical Review. Annu Rev Clin Psychol. 2018;14:1–28.

    Article  PubMed  Google Scholar 

  27. Read JR, Sharpe L, Modini M, Dear BF. Multimorbidity and depression: a systematic review and meta-analysis. J Affect Disord. 2017;221:36–46.

    Article  PubMed  Google Scholar 

  28. Ho ISS, Azcoaga-Lorenzo A, Akbari A, Davies J, Khunti K, Kadam UT, et al. Measuring multimorbidity in research: Delphi consensus study. BMJ Med. 2022;1(1):e000247.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Jr FEH. Regression Modeling Strategies 2023. Available from: http://hbiostat.org/rmsc/.

  30. Harrell FE. Regression modeling strategies: with applications to linear models, logistic regression, and survival analysis. 1st edn. New York, NY: Springer-Verlag; 2001.

  31. Riley RD, Ensor J, Snell KIE, Harrell FE Jr, Martin GP, Reitsma JB, et al. Calculating the sample size required for developing a clinical prediction model. BMJ. 2020;368:m441.

    Article  PubMed  Google Scholar 

  32. Clift AK, Dodwell D, Lord S, Petrou S, Brady M, Collins GS, et al. Development and internal-external validation of statistical and machine learning models for breast cancer prognostication: cohort study. BMJ. 2023;381:e073800.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Collins GS, Reitsma JB, Altman DG, Moons KG. Transparent reporting of a multivariable prediction model for individual prognosis or diagnosis (TRIPOD): the TRIPOD statement. BMJ. 2015;350:g7594.

    Article  PubMed  Google Scholar 

  34. von Elm E, Altman DG, Egger M, Pocock SJ, Gøtzsche PC, Vandenbroucke JP. The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement: guidelines for reporting observational studies. Lancet. 2007;370(9596):1453–7.

    Article  Google Scholar 

  35. Chen Y, Sun Y, Wang L, Xu K, Wang DW. Genetic insights into associations of multisite chronic pain with common diseases and biomarkers using data from the UK Biobank. Br J Anaesth. 2024;132(2):372–82.

    Article  CAS  PubMed  Google Scholar 

  36. Farrell SF, Kho PF, Lundberg M, Campos AI, Rentería ME, de Zoete RMJ, et al. A Shared Genetic Signature for Common Chronic Pain Conditions and its Impact on Biopsychosocial Traits. J Pain. 2023;24(3):369–86.

    Article  CAS  PubMed  Google Scholar 

  37. Foster M, Niedzwiedz CL. Associations between multimorbidity and depression among breast cancer survivors within the UK Biobank cohort: a cross-sectional study. BMC Cancer. 2021;21(1):650.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Frank P, Batty GD, Pentti J, Jokela M, Poole L, Ervasti J, et al. Association between depression and physical conditions requiring Hospitalization. JAMA Psychiatry. 2023;80(7):690–9.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Johnston KJA, Adams MJ, Nicholl BI, Ward J, Strawbridge RJ, Ferguson A, et al. Genome-wide association study of multisite chronic pain in UK Biobank. PLoS Genet. 2019;15(6):e1008164.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Macfarlane TV, Beasley M, Macfarlane GJ. Self-reported facial pain in UK Biobank study: prevalence and associated factors. J Oral Maxillofac Res. 2014;5(3):e2.

    Article  PubMed  PubMed Central  Google Scholar 

  41. McIntosh AM, Hall LS, Zeng Y, Adams MJ, Gibson J, Wigmore E, et al. Genetic and environmental risk for chronic pain and the contribution of risk variants for major depressive disorder: a family-based mixed-model analysis. PLoS Med. 2016;13(8):e1002090.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Nicholl BI, Smith DJ, Cullen B, Mackay D, Evans J, Anderson J, et al. Ethnic differences in the association between depression and chronic pain: cross sectional results from UK Biobank. BMC Fam Pract. 2015;16:128.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Qin X, Li C, Wei W, He D, Zhao Y, Cai Q, et al. Assessing the association of coffee consumption on the relationship of chronic pain with depression and anxiety. Nutr Neurosci. 2024;27(3):196–206.

    Article  PubMed  Google Scholar 

  44. Thompson W, Swain S, Zhao SS, Kamps A, Coupland C, Kuo C, et al. Causal association between subtypes of osteoarthritis and common comorbidities: a Mendelian randomisation study. Osteoarthr Cartil Open. 2023;5(4):100414.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Zhang Z, Liu L, Zhang H, Li C, Chen Y, Zhang J, et al. The genetic structure of pain in depression patients: a genome-wide association study and proteome-wide association study. J Psychiatr Res. 2022;156:547–56.

    Article  PubMed  Google Scholar 

  46. IsHak WW, Wen RY, Naghdechi L, Vanle B, Dang J, Knosp M, et al. Pain and depression: a systematic review. Harvard Rev Psychiatry. 2018;26(6):352–63.

    Article  Google Scholar 

  47. VanderWeele TJ, Rothman KJ. Formal causal models. In Lash TL, VanderWeele TJ, Haneuse S, Rothman KJ (Eds.), Modern epidemiology (4th ed., pp. 33–51). Wolters Kluwer. 2021.

  48. Coughlin SS. Recall bias in epidemiologic studies. J Clin Epidemiol. 1990;43(1):87–91.

    Article  CAS  PubMed  Google Scholar 

  49. Meng W, Adams MJ, Reel P, Rajendrakumar A, Huang Y, Deary IJ, et al. Genetic correlations between pain phenotypes and depression and neuroticism. Eur J Hum Genet. 2020;28(3):358–66.

    Article  PubMed  Google Scholar 

  50. Tang B, Meng W, Hägg S, Burgess S, Jiang X. Reciprocal interaction between depression and pain: results from a comprehensive bidirectional Mendelian randomization study and functional annotation analysis. Pain. 2022;163(1):e40–8.

    Article  CAS  PubMed  Google Scholar 

Download references

Authors’ Twitter handles

MLF: @ProfManuelaF

Funding

MF is funded by the National Health and Medical Research Council of Australia. HZ is funded by the Taishan Scholars Program of Shandong Province-Young Taishan Scholars (tsqn201909197). SF is funded by the Taishan Scholars Program of Shandong Province- Pandeng Taishan Scholars (tspd20210320). SL is a National Health and Medical Research Council Emerging Leadership Fellow (GNT2017373).

Author information

Authors and Affiliations

Authors

Contributions

LC: conceptualization, methodology, formal analysis, writing — original draft. CEA: conceptualization, methodology, writing — original draft. BS: methodology, writing — review and editing. MRR: methodology, writing — review and editing. DBA: writing — review and editing. YC: writing — review and editing. DBP: writing — review and editing. JLH: writing — review and editing. SL: writing — review and editing. MB: writing — review and editing. PRB: writing — review and editing. NKA: writing — review and editing. PHF: writing — review and editing. HZ: resources, writing — review and editing. SF: resources, writing — review and editing, supervision. MLF: writing — review and editing, supervision. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Hengxing Zhou or Shiqing Feng.

Ethics declarations

Ethics approval and consent to participate

The North West Multi-centre Ethics Committee granted ethical approval to access data from the UK Biobank (IRAS project ID: 299116), and all participants provided written informed consent.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

 Supplementary A. Details for defining pain. Supplementary B. Details for covariates. Supplementary C. Sample size calculation. Supplementary D. Details for modeling. Supplementary E. Univariate associations of the covariates with depression outcomes. Supplementary F. Predictors in the initial full model. Supplementary G. Details for model approximation. Supplementary H. Equations of the final simplified model. Supplementary I. Model performance. Supplementary J. Nomograms for secondary outcomes.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, L., Ashton-James, C.E., Shi, B. et al. Variability in the prevalence of depression among adults with chronic pain: UK Biobank analysis through clinical prediction models. BMC Med 22, 167 (2024). https://doi.org/10.1186/s12916-024-03388-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12916-024-03388-x

Keywords