Abstract

Introduction. The somatic symptom component of depression is associated with increased hospitalisation and mortality and poorer health-related quality of life (HRQOL). However, the relationship of subsets of depression symptoms with frailty and outcomes is not known. This study aimed to (1) explore the relationship between the Clinical Frailty Scale (CFS) and components of depression and (2) their association with mortality, hospitalisation, and HRQOL in haemodialysis recipients. Methods. We conducted a prospective cohort study of prevalent haemodialysis recipients, with deep bio-clinical phenotyping including CFS and PHQ-9 somatic (fatigue, poor appetite, and poor sleep) and cognitive component scores. EuroQol EQ-5D summary index assessed HRQOL at the baseline. Electronic linkage to English national administration datasets ensured robust follow-up data for hospitalisation and mortality events. Findings. Somatic (β = 0.067; 95% C.I. 0.029 to 0.104; ) and cognitive (β = 0.062; 95% C.I. 0.034 to 0.089; P<0.001) components were associated with increased CFS scores. Both somatic (β = −0.062; 95% C.I. −0.104 to −0.021; P<0.001) and cognitive (β = 0.052; 95% C.I. −0.081 to −0.024; ) scores were associated with lower HRQOL. Somatic scores lost mortality association on addition of CFS to the multivariable model (HR1.06; 95% C.I. 0.977 to 1.14; ). Cognitive symptoms were not associated with mortality. Neither the component score was associated with hospitalisation on multivariable analyses. Conclusions. Both somatic and cognitive depression symptoms are associated with frailty and poorer HRQOL in haemodialysis recipients but were not associated with mortality or hospitalisation when adjusted for frailty. The risk profile of depression somatic scores may be related to overlap with symptoms of frailty.

1. Introduction

Frailty is a syndrome of increased vulnerability to poor resolution of homeostasis after a stressor event [1]. It is associated with poor patient outcomes including mortality, hospitalisation, and disability [2]. It is commonly defined using the Clinical Frailty Scale (CFS) [3], a simple global measure of frailty based upon activities of daily living after clinical interview. In prevalent haemodialysis cohorts, estimates of frailty prevalence range from 26 to 54% using the CFS [4, 5] and it is associated with mortality and hospitalisation [4].

Depression is under-recognised in haemodialysis populations [6] and is associated with increased mortality and hospitalisation [712]. The Patient Health Questionnaire-9 (PHQ-9) has been validated in dialysis recipients with 92% sensitivity and 92% specificity for depression [13]. However, our previous work has not been able to demonstrate links between PHQ-9 depression and these outcomes (submitted for peer-review). Work elsewhere has explored dividing depression scores into somatic and cognitive symptom components [14]. The somatic symptom subset of the Beck Depression Inventory (BDI) has been associated with mortality and hospitalisation in haemodialysis recipients in Dutch [15] and Jordanian [16] haemodialysis recipients. In the Dutch cohort, depression, somatic symptoms, and cognitive symptoms were all associated with lower quality of life [15].

Frailty has been associated with poorer quality of life in Brazilian [17] and UK [18] nondialysis CKD cohorts. McAdams-DeMarco and colleagues found that self-rated fair/poor quality of life was more likely in frail haemodialysis recipients and that frailty was associated with worsening quality of life over time [19]. A decline in quality of life in frail haemodialysis recipients was also observed in a small Canadian cohort [20].

Whilst there are reports of a reciprocal relationship between frailty and depression [21, 22], no such exploration of the relationship between frailty versus the somatic and cognitive components of depression has taken place. This may be important as the somatic symptoms of depression such as tiredness, poor sleep, and lack of appetite show considerable overlap with those of frailty. Therefore, the aims of this study were to (1) explore the relationship between frailty and components of depression and (2) assess the association between the components of depression with mortality, hospitalisation, and quality of life in haemodialysis recipients.

2. Materials and Methods

2.1. Study Design

Frailty Intervention Trial iN End-Stage patients on haemodialysiS (FITNESS) is a cohort multiple randomised controlled trial (cmRCT) [23] split into two stages, for which the full protocol has been published elsewhere [24]. The study protocol was subject to favourable opinion by the South Birmingham Research Ethics Committee (Ref: 17/WM/0381) and institutional review board assessment of University Hospitals Birmingham NHS Foundation Trust (RRK6082). Here, we report data from the first stage of the FITNESS project, a cohort study with extensive baseline phenotyping for frailty and other bio-clinical parameters. The study is reported in accordance with STROBE guidelines [25].

2.2. Study Setting

This study was performed in a single nephrology centre in Birmingham, England, which oversees one in-hospital dialysis unit and ten private-provider satellite units across urban and rural settings across the West Midlands, with consequent diversity of ethnic and socioeconomic groups. Patient eligibility was ascertained using hospital electronic patient records (EPR) and liaison with healthcare professionals at each dialysis unit. Eligible patients were contacted in person on dialysis. Study investigators provided written and verbal information to prospective participants and answered any queries. Sufficient time was allowed to consider the information, before willing patients gave written informed consent to participate.

2.3. Eligibility Criteria

Inclusion criteria included adults aged 18 and over, anyone receiving regular haemodialysis for at least 3 months’ duration, and the ability to give informed consent. Patients were excluded if they received inpatient care within 4 weeks of recruitment unless for the purposes of vascular access, to avoid potential confounding by frailty associated with hospitalisation.

2.4. Baseline Assessment

All study participants underwent baseline assessment during one of their usual scheduled dialysis sessions. Participants were not assessed on the first dialysis session after the weekend interval (i.e., Monday or Tuesday), to prevent confounding by the longer interdialytic interval before assessment. Where participants dialysed twice weekly, the dialysis session after the shortest interval was chosen for baseline assessment.

Study investigations are described in detail in our methodology study [24]. Briefly, prior to connection to dialysis, participants were invited to complete a timed 4 metre walk from standing and to test bilateral hand-grip strength via a dynamometer (Takei Grip D, Takei Scientific Instruments Co. Ltd., Japan). Montreal cognitive assessments (MoCAs [26]) were also performed prior to dialysis connection. When connected to dialysis, patients were clinically interviewed, collecting demographic, social and medical history data, alongside assessment of activities of daily living (ADL) disability, and frailty-specific questionnaires. Depression symptoms were assessed using the PHQ-9 questionnaire [13]. The Physical Activity Index was derived from the GP physical activity (GPPAQ) [27] questionnaires via a validated formula to give a global measure of participant physical activity. [27, 28] Electronic patient records were interrogated for data upon biochemistry, dialysis adequacy, comorbidities, prescribed medications, alongside dialysis vintage, and previous renal replacement therapies. Self-reported change in health (henceforth “health change”) was assessed with the question “How has your health changed in the last year?” with potential responses of “Better,” “The Same,” or “Worse” [29]. The English Index of Multiple Deprivation 2015 (IMD) was used to assess socioeconomic deprivation [30]. This is a national multiple deprivation model calculated from multiple socioeconomic data points at a postal code (Zip code) level. A composite score is obtained and split into quintiles of deprivation, with 1 representing the most deprived and 5 the least deprived area, respectively.

PHQ-9 somatic component scores included questions 3 (“Trouble falling or staying asleep or sleeping too much”), 4 (“Feeling tired or having little energy”), and 5 (“Poor appetite or overeating”). Cognitive component scores included questions 1 (“Little interest or pleasure in doing things”), 2 (“Feeling down, depressed, or hopeless”), 6 (“Feeling bad about yourself or that you are a failure or have let yourself or your family down”), 7 (“Trouble concentrating on things, such as reading the newspaper or watching television”), 8 (“Moving or speaking so slowly that other people could have noticed? Or the opposite-being so fidgety or restless that you have been moving around a lot more than usual”), and 9 (“Thoughts that you would be better off dead or of hurting yourself in some way”) [14].

The Clinical Frailty Scale was obtained by interpretation of ADL questionnaire responses, with possible responses of 1–9. A CFS of 1 represented very fit and 8 severely frail. A CFS score of 9 was attributed to those who were terminally ill but not overtly frail. A CFS score of ≥5 was considered frail.

2.5. Outcomes

Mortality and cause of death data were obtained by electronic record linkage to the Office of National Statistics (ONS), a UK-wide database of death certification. Hospitalisation data were obtained via Hospital Episode Statistics (HES), a clinical coding database containing all secondary care episodes in any English NHS hospital. Hospital admissions were defined as any hospital episode lasting ≥1 night. Transfers between hospitals were treated as one continuous admission, and length of stay in such episodes was calculated from admission at the initial hospital to discharge from the final hospital.

Health-related quality of life was assessed using Euroqol EQ-5D-3L. To allow global assessment of HRQOL, EQ-5D-3L responses were converted into a single-measure EQ-5D summary index via a standardised formula, validated in UK populations [31]. The EQ-5D summary index score ranges between 1 (if no HRQOL deficits reported) and −0.716 (if extreme deficit in every domain).

2.6. Recruitment

A power calculation was based upon US data demonstrating an association of frailty with an adjusted risk ratio of 2.24 for 1-year mortality and 1.56 for 1-year mortality/hospitalisation in haemodialysis recipients [32]. We assumed a nonfrail risk of mortality and mortality/hospitalisation to be 5% and 40%, respectively. Powered to 0.8 and with a confidence interval of 0.95, a sample size of 602 was therefore calculated for 1-year mortality and 150 patients for 1-year mortality/hospitalisation. Upon discussion and agreement with the sponsor, however, recruiting 602 participants was not felt to be feasible in this single centre. As such a revised target of 500 participants was set with follow-up beyond 1 year.

2.7. Statistics

Statistical analysis was performed using STATA 17 (StataCorp 2019, Stata Statistical Software: Release 17, College Station, TX: StataCorp LLC) and R version 4.0.4 (R Foundation for Statistical Computing, Vienna, Austria). Categorical data were presented as numbers and percentages, and continuous variables were reported as medians and interquartile ranges (IQRs).

Time-to-event outcomes were analysed with Cox’s proportional hazards model. The proportional hazard assumption was checked via interrogation of the log-negative-log plots of the within-group survivorship functions versus log time. Furthermore, we compared Kaplan–Meier (observed) with Cox (expected) survival curves for study variables (reported as hazard ratios (HRs) with 95% confidence intervals (CIs)).

Linear regression analyses explored the relationship between continuous frailty and PHQ-9 scores and between these same scores and EQ-5D summary index. The linearity assumption was checked by visually comparing plots of observed values by linear and LOWESS fit and by plotting observed versus predicted residual values. Linear regressions were performed unadjusted and as a series of adjusted models based upon a priori covariables selected for known or suspected relationship with the outcome of interest. Due to the number of covariables, adjusted models were constructed in a predetermined stepwise manner. For frailty, model 1 included depression, age, gender, and ethnicity. Model 2 added education level, social support, and IMD quintile. Model 3 added Charlson Index (CKD excluded), MoCA score, smoking status, self-rated change in health, and overall health. Model 4 added use of walking aids, Physical Activity Index, slow timed walk, and low grip strength.

Regressions for the EQ-5D summary index were performed unadjusted and subject to a separate set of a priori covariable models, based upon known or suspected relationship with HRQOL, and covariables found to significantly associate with frailty and/or PHQ-9 scores were also added to reduce confounding. Due to the large number of covariables identified, models were constructed in a predetermined stepwise approach. Model 1 included age, gender, ethnicity, education level, social support, IMD quintile, and employment status. Model 2 added to these haemodialysis vintage, Charlson comorbidity index (CKD omitted), haemoglobin, Kt/V, and current use of antidepressant medication. Model 3 added use of walking aids, slow walking (or inability to walk), and Physical Activity Index. Model 4 added EQ self-rated health today (continuous score from 0 to 100) and self-rated health change.

Count data were explored by negative binomial regression, death-censored and offset by length of follow-up, to give incidence rate ratios (IRRs). Negative binomial distribution was confirmed by interrogation of means and variances and visual inspection of observed versus expected distribution plots. Zero-truncated negative binomial regressions were performed for nights per admission, as by definition these could not equal zero.

We performed both unadjusted and adjusted negative binomial and Cox regressions. Covariables for adjusted analyses were selected a priori based upon a proven or suspected relationship with hospitalisation and/or mortality (age, sex, ethnicity (grouped into white, south Asian, black, and other ethnicities), body mass index, index of multiple deprivation, Charlson comorbidity index (CKD omitted), number of hospitalisation episodes, number of medications, smoking status, serum albumin, use of walking aids, dialysis vintage, self-reported change in health, and kidney transplant wait-listing). Furthermore, adjusted models were performed using the aforementioned covariables plus the addition of the CFS score. No transplant-listed participants died within 1 year of recruitment, so this covariable was omitted in the final logistic regression models for 1-year mortality.

A dummy variable was used to handle missing IMD quintile data. All other covariables were assumed missing at random as they had <1% data missing and were therefore handled via listwise deletion. Statistical significance was set at a value <0.05.

3. Results

3.1. Study Cohort Demographics

Figure 1 shows the PRISMA study flow of participant recruitment to the FITNESS study, with 485 prevalent haemodialysis patients with baseline frailty assessments and data linkage. Follow-up was 678 days (interquartile range: 531–812 days), with minimum potential follow-up of 365 days from recruitment. Baseline demographics are described in detail elsewhere [5]. Table 1 shows key demographics stratified by frailty status at study recruitment.

3.2. PHQ-9 Somatic and Affective Scores

The median PHQ-9 somatic and cognitive scores were 3 (IQR 1, 6) and 1 (IQR 0, 5), respectively. Scores ranged from 0 to 9 for the PHQ-9 somatic score and 0–18 for the PHQ-9 cognitive score.

3.3. Relationship with Frailty

Figures 2 and 3 show that both somatic and cognitive scores were positively associated with the CFS score on simple and all multiple linear regression models. Full results of the final models for somatic and cognitive scores are shown in Supplementary Tables 1 and 2.

3.4. Association with Mortality

PHQ-9 somatic component was associated with mortality on univariable analysis (HR 1.10; 95% C.I. 1.02, 1.17; ), but cognitive component was not (HR 1.00; 95% C.I. 0.950, 1.05; ). Figure 4 shows that the somatic component retained this association on multivariable analysis with CFS omitted but lost the significance on addition of CFS to the model. The affective component was not associated with mortality regardless of CFS inclusion. Full adjusted models are shown in Supplementary Tables 36.

3.5. Association with Hospital Admissions

The PHQ-9 somatic score was associated with increased rates of hospital admissions on univariable analysis (IRR 1.05; 95% C.I. 1.01, 1.10; ) but not upon multivariable analysis (IRR 1.00; 95% C.I. 0.96, 1.04; ).

PHQ-9 cognitive scores were not associated with rates of hospital admissions on either univariable (IRR 1.02; 95% C.I. 1.00, 1.06; ) or multivariable analyses (IRR 0.98; 95% C.I. 0.95, 1.01; ). Fully adjusted model results for somatic and cognitive scores are shown in Tables 2 and 3, respectively; of the a priori covariables in these models, age, CFS, the Charlson Index, number of previous admissions, and walking aid use were all associated with higher rates of admissions, while black ethnicity was associated with lower admission rates when included in both cognitive and somatic component score models. Supplementary Tables 7 and 8 show that omission of the CFS score from multivariable models did not significantly alter results.

3.6. Association with Quality of Life

Increases in both somatic and cognitive components of PHQ-9 were associated with lower EQ Summary Index scores on fractional regression, as shown in Table 4. Furthermore, increasing CFS scores were also associated with significant reductions in EQ Summary Index on all models, independent of somatic and cognitive component scores. Fully adjusted fractional regression model results are shown in Supplementary Tables 9 and 10.

4. Discussion

Both depression and frailty have been associated with mortality and hospitalisation in dialysis patients [712]. However, previous work in the FITNESS cohort has shown that depression was not significantly associated with either of these important negative outcomes [33]. The somatic and cognitive components of depression have also been associated with hospitalisation and mortality [15, 16], though their relationship with frailty has not been explored. Here, we show that both the somatic and cognitive components of depression are associated with increasing frailty. Only somatic symptoms were associated with increased mortality and hospitalisation in univariable analyses but lost these associations on multivariable analyses. Both cognitive and somatic symptoms were associated with lower health-related quality of life. Lack of association between somatic depression symptoms and hospitalisation or mortality contradicts reports elsewhere, but the association with poorer quality of life indicates that these symptoms remain important to haemodialysis recipients.

Work within the FITNESS cohort has shown that while PHQ-9 scores associate with frailty, they predict neither admissions nor mortality in prevalent haemodialysis recipients [33]. Schouten and colleagues showed a differential risk profile between the somatic and cognitive components of the BDI [15]. Somatic, but not cognitive, symptoms were associated with all-cause mortality, whereas both dimensions were associated with hospitalisation and poorer HRQOL [15]. Khalil et al. found a significant association between both somatic and cognitive component scores of the BDI with mortality and hospitalisation in a Jordanian haemodialysis cohort [16]. In the FITNESS cohort, however, the cognitive component did not associate with hospitalisation or mortality. Furthermore, while the somatic component did associate with mortality on multivariable analysis, it lost this association on addition of frailty into the model. We must exercise caution when comparing data across national and socio-cultural boundaries, obtained using different methodologies and heterogenous depression scores. However, neither the Dutch nor Jordanian cohorts included a frailty measure in their analyses [15, 16]. Given the overlap between symptoms of frailty and somatic depression symptoms, we may speculate that the somatic depressive symptoms may represent a surrogate marker for frailty in haemodialysis recipients. This may explain the loss of mortality association upon the addition of a de facto measure of frailty and a more powerful associate with negative outcomes. These novel findings indicate some of the complexity inherent to frailty assessment in heavily comorbid populations such as haemodialysis.

To our knowledge, FITNESS is amongst the first studies to explore relationships between frailty and depression symptom subsets in a large prospective haemodialysis cohort. Strengths of the study include the large cohort size, prospectively recruited, with diversity of population representative of our local populace [34]. The cohort is deeply phenotyped, allowing for a broad range of medical, social, and lifestyle factors to be included in our analyses. Furthermore, electronic data linkage ensures robust data capture of hospitalisation and mortality. However, we must advise caution in applying our data to non-English populations; validation of our findings elsewhere is required. We have adjusted for many potential confounders in our analyses, but we must be cautious about overfitting the models to our cohort. Covariables were added to our models in a stepwise manner to mitigate for this. As complexity of the models increased, the effect of the independent variable was attenuated, but we would argue that the inferences remained the same regardless of the model used. Limitations also include the single baseline data collection for frailty and depression phenotyping, both frailty and depression are dynamic states, and serial measurements would arguably improve both accuracy and analytical detail [10, 35]. The method of obtaining CFS was not subject to MDT discussion, which represents a deviation from the original CFS validation cohort [3]. However, we suggest that our approach is comparable to use of the CFS in clinical practice [36]. The EQ-5D Summary Index for HRQOL is validated in UK populations, but the relationship to other quality of life measures is not clear. Finally, as with all observational data, we report associations rather than causation, and we must be cautious when applying these findings to the individual haemodialysis recipient in clinical practice.

To conclude, both somatic and cognitive components of depression are associated with frailty and poorer HRQOL in haemodialysis recipients, but they are not associated with mortality or hospitalisation on fully adjusted models including frailty. These data may suggest that there is an overlap between frailty and depression in their associations with negative patient outcomes. Further work is warranted to better understand and distinguish individual versus cumulative contributions from overlapping comorbidities towards adverse outcomes in prevalent haemodialysis patients.

Data Availability

Data underlying this manuscript will be made available from the corresponding author upon reasonable request.

Additional Points

We performed a prospective cohort study investigating the relationship between symptom subsets of the PHQ-9 score with the Clinical Frailty Scale in prevalent haemodialysis recipients. Both somatic and cognitive component scores of PHQ-9 were associated with frailty and poorer HRQOL but were not associated with mortality or hospitalisation when adjusted for frailty.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by Queen Elizabeth Hospital Charity, Fund Number 17-3-886.

Supplementary Materials

Supplementary Table 1: multiple linear regression of the CFS score associated with the PHQ-9 somatic score. Fully adjusted model 4. Supplementary Table 2: multiple linear regression of the CFS score associated with the PHQ-9 cognitive score. Fully adjusted model 4. Supplementary Table 3: fully adjusted Cox regression model of mortality associated with the PHQ-9 somatic component score. CFS omitted. Supplementary Table 4: fully adjusted Cox regression of mortality associated with the PHQ-9 somatic component score. CFS included. Supplementary Table 5: fully adjusted Cox regression of mortality associated with the PHQ-9 cognitive component score. CFS omitted. Supplementary Table 6: fully adjusted Cox regression model of mortality associated with the PHQ-9 cognitive component score. CFS included. Supplementary Table 7: incidence rate ratios of hospital admissions associated with the PHQ-9 somatic component score. Adjusted model with CFS omitted. Supplementary Table 8: incidence rate ratios of hospital admissions associated with the PHQ-9 cognitive component score. Adjusted model with CFS omitted. Supplementary Table 9: fractional regression of association of PHQ-9 somatic component scores with EQ Summary Index: fully adjusted model 4. Supplementary Table 10: fractional regression of association of PHQ-9 cognitive component scores with EQ Summary Index: fully adjusted model 4. (Supplementary Materials)