Introduction
The relationship between depression and mortality has been frequently studied. Indeed, three systematic reviews in the general population (Wulsin et al. Reference Wulsin, Vaillant and Wells1999; Cuijpers et al. Reference Cuijpers, Vogelzangs, Twisk, Kleiboer, Li and Penninx2014) and in older people (Schulz et al. Reference Schulz, Drayer and Rollman2002) cover almost all the available literature on the issue, with compelling evidence of an increased risk of dying linked to depression. In addition, depressive symptoms have also been associated with mortality in community-dwelling elderly (Sun et al. Reference Sun, Schooling, Chan, Ho and Lam2011; White et al. Reference White, Zaninotto, Walters, Kivimaki, Demakakos, Shankar, Kumari, Gallacher and Batty2015). In residential and nursing homes, some (Rovner et al. Reference Rovner, German, Brant, Clark, Burton and Folstein1991; O'Connor & Vallerand, Reference O'Connor and Vallerand1998; Barca et al. Reference Barca, Engedal, Laks and Selbaek2010; Kane et al. Reference Kane, Yochim and Lichtenberg2010; Drageset et al. Reference Drageset, Eide and Ranhoff2013) but not all cohort studies (Cohen-Mansfield et al. Reference Cohen-Mansfield, Marx, Lipson and Werner1999; Parmelee et al. Reference Parmelee, Katz and Lawton1992; Cuijpers, Reference Cuijpers2001; Sutcliffe et al. Reference Sutcliffe, Burns, Challis, Mozley, Cordingley, Bagley and Huxley2007) have found a positive association of depression and depressive symptoms with mortality. Most of these studies were conducted over short follow-up periods and none has properly examined the potential role of detection in the association between depression and mortality, though one study did provide some related data (Rovner et al. Reference Rovner, German, Brant, Clark, Burton and Folstein1991). In a previous report, we found that 42% of depressed residents were undetected, corresponding to 18% of the whole institutionalised population (Damián et al. Reference Damián, Pastor-Barriuso and Valderrama-Gama2010). Accordingly, the aims of this study were to measure the association of depression with long-term all-cause mortality and evaluate a potential effect modification according to its detection status in a representative population sample residing in facilities for older people.
Method
Study population
This cohort study used mortality follow-up data from a baseline survey conducted from June 1998 through June 1999 in a stratified cluster sample of residents aged 65 years or older in residential and nursing homes of Madrid, Spain. We initially selected 25 public/subsidised and 30 private institutions with probability proportional to their size, and then sampled ten men and ten women from each public/subsidised facility, and 5 men and 5 women from each private facility. Due to refusal or prolonged absence, 85 residents could not be included, leading to an overall response rate of 89% (715 out of 800 sample residents). Thirty nine subjects were randomly substituted with residents of the same facility and sex, yielding a total of 754 baseline interviews. In the present study we excluded subjects with severe cognitive impairment since they were likely unable to remember or report effective depressive symptoms.
The Carlos III Institutional Review Board approved the study. Informed consent was obtained verbally and documented from all study participants or their next of kin.
Baseline data collection
Structured questionnaires were administered by trained geriatricians to all residents, their main caregivers, and the facility physicians to collect baseline data on sociodemographic characteristics, medical history, depressive symptoms, functional dependency, and cognitive status.
Physician's diagnosis of depression was ascertained by interviewing facility physicians (or nurses in 8% of residents) with access to medical history. Antidepressants used for the preceding 7 days (code N06A of the World Health Organization Anatomical Therapeutic Chemical Classification) were noted by reviewing medical records. A ten-item version of the Geriatric Depression Scale (GDS) was administered verbally to each resident to assess depressive symptoms over the previous week without focusing on physical complaints (D'Ath et al. Reference D'Ath, Katona, Mullan, Evans and Katona1994). This self-report measure has shown adequate diagnostic accuracy in the institutional setting, including residents with mild to moderate cognitive impairment. The overall discrimination of the ten-item GDS was 0.86 for identifying clinical major or minor depression among nursing home patients, with sensitivity rates of 75–86% and specificity rates of 70–77% for the optimal cut-off score of 4 or higher (Shah et al. Reference Shah, Phongsathorn, Bielawska and Katona1996; Jongenelis et al. Reference Jongenelis, Pot, Eisses, Gerritsen, Derksen, Beekman, Kluiter and Ribbe2005). Although the ten-item GDS has not been validated for Spanish population, translation and adaption of longer 15- and 30-item versions to the Spanish language were straightforward and showed similar psychometric properties to those of the original questionnaires (Izal & Montorio, Reference Izal and Montorio1993; Martínez de la Iglesia et al. Reference Martínez de la Iglesia, Onís Vilches, Dueñas-Herrero, Albert Colomer, Aguado Taberné and Luque Luque2002). (See supplementary material.)
Depression was defined as a physician's diagnosis of the condition, use of antidepressants, or a score of 4 or higher on the ten-item GDS. Residents with depression were considered to be detected if they had a physician's diagnosis or had been prescribed antidepressants, irrespective of depressive symptoms; whereas depressed residents were deemed to be undetected if they scored 4 or higher on the ten-item GDS, without documented diagnosis or treatment.
Chronic conditions, including cancer, obstructive pulmonary disease, arrhythmias, hypertension, ischaemic heart disease, congestive heart failure, peripheral arterial disease, stroke, diabetes, anemia, Alzheimer's disease, other dementias, Parkinson's disease, epilepsy, anxiety disorders and arthritis, were gathered from the physician's interview. Functional dependency in performing basic activities of daily living was assessed by residents or their main caregivers (if assigned, 45%) using the modified Barthel index (Shah et al. Reference Shah, Vanclay and Cooper1989). Residents were classified as functionally independent (100 points), mild to moderate dependency (61–99 points), and severe to total dependency (0–60 points) (Shah et al. Reference Shah, Vanclay and Cooper1989).
Cognitive status was evaluated using both the Short Portable Mental Status Questionnaire (SPMSQ, 0–10 errors) (Pfeiffer, Reference Pfeiffer1975), which was adapted to the institutional setting, and the Minimum Data Set Cognition Scale (MDS-COGS, 0–10 points) (Hartmaier et al. Reference Hartmaier, Sloane, Guess and Koch1994; Gruber-Baldini et al. Reference Gruber-Baldini, Zimmerman, Mortimore and Magaziner2000), which obtained an assessment from the main caregivers based on selected Minimum Data Set questions. This latter scale was used because the SPMSQ could only be administered to 61% of residents for logistic reasons. Residents were classified as normal cognition (≤2 education-adjusted SPMSQ errors and ≤1 MDS-COGS points), mild to moderate cognitive impairment (3–7 SPMSQ errors and ≤8 MDS-COGS points, or ≤7 SPMSQ errors and 2–8 MDS-COGS points), and severe cognitive impairment (≥8 SPMSQ errors or ≥9 MDS-COGS points).
Mortality ascertainment
Mortality was ascertained by requesting data on residents’ vital status to the participating facilities and through linkage to the Spanish National Death Index, which includes all deaths registered in Spain since 1987 (Ministerio de Sanidad Servicios Sociales e Igualdad, 2014). Residents contributed follow-up time from their 1998–1999 baseline interview or age 75 years, whichever occurred later, until death, age 105 years, or 15 September 2013, whichever occurred first. We restricted follow-up to the age interval 75–104 years because there were only 33 and 2 deaths in residents aged 65–74 and ≥105 years, respectively.
Statistical analysis
Since depressive symptoms may reflect frailty conditions near death (reverse causation bias), early deaths within the first year from the baseline interview were analyzed separately from the remaining follow-up. Hazard ratios for mortality and 95% confidence interval (CI) for residents with detected and undetected depression compared with non-depressed residents were estimated using Cox proportional hazards models with age as the time scale (Thiebaut & Benichou, Reference Thiebaut and Benichou2004) and adjusted for sex, type of facility (public, subsidised or private), number of chronic conditions other than depression (0–1, 2–3 or ≥4), and functional dependency (independent, mild/moderate or severe/total). We tested for effect modifications by including interactions of detected and undetected depression with each of the above baseline covariates. Non-parametric survival curves for residents with no, detected, and undetected depression were estimated as the baseline survival functions from a depression-stratified Cox model adjusted to the overall weighted percentages of baseline confounders. To allow for non-proportional hazards over age, we fitted age-dependent Cox models with interactions of detected and undetected depression with a restricted quadratic spline function of age with knots at 80, 90 and 100 years (Hess, Reference Hess1994).
We also fitted a Weibull model with different location and scale parameters for each depression group and adjusted to the overall weighted percentages of baseline confounders. This model provided similar survival curves to those obtained by non-parametric methods, but allowed estimating the absolute years of life gained and lost for residents with detected and undetected depression compared with non-depressed residents as the difference in percentiles from their estimated Weibull distributions, with 95% CI derived from delta methods (Cox et al. Reference Cox, Chu, Schneider and Munoz2007).
Sensitivity analyses were performed to evaluate the robustness of the results to alternative definitions of depression and further adjustments for potential confounders. Specifically, we replicated the analyses excluding antidepressants from the definition of depression, as well as using a stricter cut-off score of 5 or higher on the ten-item GDS for the presence of depressive symptoms. In addition, although residents with severe cognitive impairment were excluded from the analyses, the association between depression and mortality was further adjusted for residual differences in cognitive functioning (unimpaired or mild/moderate impairment) in the subset of residents with available data.
Due to the complex design and the different selection probabilities (residents in public/subsidised facilities and men were oversampled) all analyses were weighted to the underlying population distribution and accounted for the effect of stratification and clustering on point and interval estimates (using Taylor-linearised variance estimation). Analyses were performed using Stata, version 13.1 (Stata Corp., College Station, Texas) and R, version 2.15 (R Foundation for Statistical Computing, Vienna, Austria).
Results
From the 754 participants in the baseline survey, we excluded 54 residents due to insufficient information on their mortality status at the end of follow-up, 33 residents who died before 75 years of age and 76 residents with severe cognitive impairment. Thus, the final sample included 591 residents who were followed up for mortality from 75 to 104 years of age. During 3575 person-years of follow-up (median/maximum of 4.8/15.2 years), 499 participants died, including 70 deaths within the first year of follow-up. Mortality rates were higher in men, residents in public/subsidised facilities, and subjects with mild to moderate cognitive impairment, and increased steadily with age and degree of functional dependency (Table 1).
Population distribution and mortality rates by age interval and baseline characteristics of residents in nursing homes of Madrid, Spain, 1998–1999 to 2013
* Unweighted counts and weighted percentages. Figures across age intervals add up to more than the overall sample size because subjects may contribute to different age intervals during follow-up.
† Weighted mortality rates per 1000 person-years and 95% CI.
The baseline prevalences of detected and undetected depression were 25.9 and 18.8%, respectively, adding up to an overall prevalence of depression of 44.7%. The estimates of the association of depression-related variables with mortality during the first year of follow-up were imprecise and mostly compatible with a null effect (Table 2). For the remaining long-term follow-up period, baseline physician-reported depression and antidepressant use were associated with 39 and 32% reduced risks of death in multivariate-adjusted models, respectively. In contrast, mortality risk increased by 24% among residents with baseline depressive symptoms (GDS score ≥4) compared with those without symptoms (Table 2). As a result of these opposite associations, the composite depression variable showed virtually no effect on mortality (hazard ratio 0.87, 95% CI: 0.70–1.08). However, when depression was disaggregated by its detection status, residents with detected depression registered a 37% lower mortality risk than those without depression, whereas residents with undetected depression displayed a 35% higher, not statistically significant, mortality risk (Table 2). The hazard ratio for mortality comparing undetected with detected depression was 2.14 (95% CI: 1.34–3.44). Accordingly, the adjusted survival curves display consistently lower mortality over time in residents with detected depression compared with undepressed residents (p for homogeneity of survival curves = 0.005) and substantially higher mortality in residents with undetected depression compared with those with no depression (p = 0.06) and detected depression (p = 0.002) (Fig. 1). The hazard ratio for detected depression v. no depression remained fairly stable below 1 over the entire age range (p for proportional hazards = 0.57), unlike that for undetected depression (p = 0.004), which displayed a markedly increased risk of death at ages below 85 years and a dilution in the strength of the association at older ages, mainly due to the differential depletion of residents with undetected depression (Fig. 2).
Adjusted non-parametric survival curves from age 75 for residents with no depression, detected depression, and undetected depression in nursing homes of Madrid, Spain, 1998–1999 to 2013.
Age-dependent hazard ratios for mortality (solid curves) and 95% CI (dashed curves) in residents with detected and undetected depression compared with residents without depression in nursing homes of Madrid, Spain, 1998–1999 to 2013. The histograms represent the weighted numbers of person-years by age interval among residents with no depression (shaded bars), detected depression (left white bars) and undetected depression (right white bars).
Association of baseline depression-related variables with mortality risk by follow-up period among residents in nursing homes of Madrid, Spain, 1998–1999 to 2013
* Unweighted counts and weighted percentages.
† Hazard ratios and 95% CI were obtained from Cox models splitting each subject's exposure time into the first year and the remaining follow-up and adjusted for age, sex, type of facility, number of chronic conditions other than depression and functional dependency.
‡ Score of 4 or higher on the ten-item GDS.
§ Detected depression was defined as a physician's diagnosis of depression or antidepressant use, and undetected depression as a score of 4 or higher on the ten-item GDS without documented diagnosis or treatment.
Using a Weibull model, which showed a very good fit to the non-parametric survival curves (appreciable in the left panel of Fig. 3), it was estimated that the lower, median and upper quartiles of life expectancy at age 75 for non-depressed residents were 80.6, 84.6 and 89.6 years, respectively. By the same quartiles, residents with detected depression gained 0.4, 1.8 and 4.0 years of life, whereas those with undetected depression lost 4.3, 6.3 and 8.0 years, respectively (right panel of Fig. 3). The median life expectancy decreased by 8.1 years (95% CI: 2.3–13.9 years) in residents with undetected depression compared with those detected.
Years of life gained and lost due to detected and undetected depression among residents in nursing homes of Madrid, Spain, 1998–1999 to 2013. The left panel shows the parametric survival curves (smooth lines) obtained from a Weibull model with different location and scale parameters for each depression group, together with the corresponding non-parametric survival curves (step functions) for comparison. The right panel shows the absolute years of life gained and lost (solid curves) and their 95% CI (dashed curves) for residents with detected and undetected depression compared with non-depressed residents, as estimated from the Weibull model.
In subgroup analyses, the long-term beneficial effect of detected depression on mortality was higher in women (p for interaction = 0.002), private facilities (p = 0.09) and residents with mild-to-moderate cognitive impairment (p = 0.003). The harmful effect of undetected depression was quite homogeneous across subgroups defined by baseline covariates (Fig. 4).
Hazard ratios for mortality (squares with area inversely proportional to the variance) and 95% CI (horizontal lines) comparing detected and undetected depression with no depression in pre-specified subgroups of residents in nursing homes of Madrid, Spain, 1998–1999 to 2013.
In sensitivity analyses, we excluded antidepressants from the definition of depression and results were strengthened, with hazard ratios for long-term mortality comparing detected and undetected depression with no depression of 0.69 (95% CI: 0.51–0.94) and 1.50 (95% CI: 1.11–2.03), respectively. These associations did not change appreciably after adjustment for antidepressant use (hazard ratios 0.72 and 1.47). Also, we used a more stringent cut-off score of 5 or higher on the ten-item GDS, which resulted in a slight attenuation of the association between undetected depression and mortality (hazard ratio 1.24, 95% CI: 0.78–1.98). Finally, further adjustment for cognitive impairment (unimpaired or mild/moderate) in 78.9% of residents with available data yielded similar but less precise results (hazard ratios 0.59, 95% CI: 0.41–0.83, for detected depression and 1.43, 95% CI: 0.93–2.20, for undetected depression).
Discussion
The association between depression and mortality has been consistent in most studies of community-dwelling older people but less clear in the case of the institutionalised population. However, this study found somewhat conflicting results, with a null overall association. In addition, physician-reported depression was associated with lower mortality. In-depth examination of the potential confounding effect of additional variables, such as number of medicines and several diseases, led to no appreciable change in the findings. Nonetheless when our analysis focused on depressive symptoms, we found an increase in risk in residents with symptoms. The analyses were repeated using a stricter GDS cut point of ≥5 (over 10), thus increasing specificity and decreasing sensitivity. This led some individuals to change from undetected to non-depressed with an expected dilution of the association – which in no way suffices to change our initial conclusion. A higher specificity can achieve less false positives but at the expense of increasing the probability of false negatives. In line with our results false negatives are a worse consequence because subjects can remain unattended and we show this can increase their mortality risk.
When the depressed group was broken down into detected and undetected, undetected individuals displayed a higher mortality than did those with no depression. In addition, it is nevertheless striking that the risk would be lower among those with detected depression than it was in the non-depressed group. Although we are unable to furnish an adequate explanation for this, we believe that depression detection might be part of a more general mechanism which links closer care and supervision, greater attention, responsiveness, etc. to lower mortality. In view of the fact that residents for whom detection was easier would share a better prognosis, adjustment was made for the main prognostic variables but the results remained unchanged. Detection may prove difficult in older persons, particularly in nursing homes where the quiet resident can be overlooked (Fenton et al. Reference Fenton, Raskin, Gruber-Baldini, Menon, Zimmerman, Kaup, Loreck, Ruskin and Magaziner2004). One might also speculate that, once residents have been assigned a depression diagnosis (and a treatment) and have assimilated it, this could trigger motivational changes in their attitude towards better self-care, compliance with treatments (depression and other conditions), and increased participation in activities (preventive, rehabilitative and social). In brief, detection could help break the vicious cycle of depression leading to negative conduct in self-care/activities/behaviours, which in turn leads to illness and ultimately to a return to a depressed mood.
We wish to highlight the importance of including the non-depressed group in the analytical framework not only because its absence yields a notable loss of information, but also because the non-depressed state is a part of the natural history involving detection: non-depressed, depressed but undetected and depressed but detected. In addition, the finding of non-depressed group showing a somewhat higher risk compared with detected, though admittedly incidental, is worthy of further investigation because it might be closely related to relevant aspects of care in these facilities.
For the purposes of this study, subjects using antidepressants were deemed to be depressed and detected. Admittedly, this choice is open to discussion because antidepressant use may be indicated for conditions other than depression; albeit this condition was the predominant reason for prescription in our study (63.1% of residents with prescribed antidepressants had physician's diagnosis of depression). Thus, we performed sensitivity analyses excluding antidepressants from depression definition and findings were even more marked.
According to our results it is possible that detection of depression be preventive only in women. This is in agreement with a finding on subthreshold depression in older community residents (Hybels et al. Reference Hybels, Pieper and Blazer2002). Men are usually more reluctant to admit or assimilate the diagnosis, leading to worse resources for resilience or coping (Hinton et al. Reference Hinton, Zweifach, Oishi, Tang and Unutzer2006; Murray et al. Reference Murray, Banerjee, Byng, Tylee, Bhugra and Macdonald2006). In fact, the general pattern (of lower mortality in individuals with detected depression) was not confirmed in men, among whom depression detection did not seem to make a notable difference (Fig. 4).
Among the most comparable studies, Drageset et al. (Reference Drageset, Eide and Ranhoff2013) recently reported depression scores associated with 5-year mortality in a sample of cognitively intact residents from 30 nursing homes in Norway. Another of the few nursing-home-based studies to measure mortality with a long follow-up (even if undertaken in a single facility) failed to find an association between depressed affect, as evaluated by nurses, and mortality, though some association was observed in the short term (1 year) (Cohen-Mansfield et al. Reference Cohen-Mansfield, Marx, Lipson and Werner1999). It should be noted that we found no association in the 1-year follow-up in the present study. The effects of detection/non-detection are discernible with the passage of time, and indeed ours is the only study to give an adequate description of how the hazard ratio changes with time. In primary care, Bogner et al. (Reference Bogner, Morales, Reynolds, Cary and Bruce2012) found that persistent depressive symptoms were related to mortality and demonstrate that repeated assessment of depression may be particularly important in preventing excess mortality associated to depression in older adults. We speculate that results even more favourable would be expected in nursing home population.
Potential pathways in the association of depression and mortality are multiple and complex. Some of them are discussed below.
Rovner et al. (Reference Rovner, German, Brant, Clark, Burton and Folstein1991) suggested some mechanisms, including insomnia, poor nutrition due to anorexia, immobility from psychomotor retardation and altered immune function. Cuijpers (Reference Cuijpers2001) also mentioned immune function and added interference with patients’ motivation towards recovery. Somewhat similar factors are mentioned in Kane et al. (Reference Kane, Yochim and Lichtenberg2010) (poor self-care, less physical activity and immune system functioning). Kiecolt-Glaser & Glaser (Reference Kiecolt-Glaser and Glaser2002) review the evidence on the role of immune function, adding a simple but important message, the interaction of depression with ageing on enhancing risks for morbidity and mortality. This would be particularly relevant in frail elders.
Strengths and limitations
Our study was conducted on a representative sample of individuals living in institutions (enhancing external validity), was based on a substantial sample size (limiting random error), and collected information on most of the relevant variables (limiting confounding and allowing for assessment of effect modification). Even so, some degree of depression misclassification is possible, in that we did not follow a structured assessment based on DSM-IV or ICD-10 criteria. Nonetheless we think that the group scoring 4 or more in the GDS-10 is likely to include not only a majority of depression cases as determined by that clinical criteria, but also a number of subclinical or subthreshold cases (i.e., a condition in which a person has depressive symptoms, but does not meet the criteria for a depressive disorder), likewise found to be serious and chronic (Beekman et al. Reference Beekman, Geerlings, Deeg, Smit, Schoevers, de Beurs, Braam, Penninx and Van Tilburg2002) and associated with mortality as well (Cuijpers et al. Reference Cuijpers, Vogelzangs, Twisk, Kleiboer, Li and Penninx2013).
Thus, taking into account our results, perhaps patients scoring 4 or more GDS-10 symptoms should merit the label of depression. At any rate, one of the main messages derived from our findings is that depressive symptoms, not declared as depression by physicians, may lead to higher mortality. As for the outcome, we believe that some deaths might not have been identified, something that would eventually generate non-differential misclassification and, in general, lead to dilution of measured associations. Finally, we recommend that these results be generalised only to subjects without severe cognitive impairment.
Conclusions
The long-term all-cause mortality risk associated with depression depends on the condition's detection, with a better prognosis in those cases where it is detected. Had detection played a causal role, the absolute numbers of years of life lost due to undetected depression that our study estimates may be remarkable. Regarding implications on the care of nursing homes, it would seem imperative to adopt a proactive approach to the detection and treatment of depression because, whatever its causal role vis-à-vis mortality, what is undeniable is that depression leads to avoidable suffering. We endorse the recommendation that GDS be used for subjects who are not cognitively impaired, followed by initiation of treatment and monitoring of the outcome (Boyle et al. Reference Boyle, Roychoudhury, Beniak, Cohn, Bayer and Katz2004), and the Cornell Scale for Depression in Dementia for subjects with cognitive impairment (Alexopoulos et al. Reference Alexopoulos, Abrams, Young and Shamoian1988).
Supplementary Material
To view supplementary material for this article, please visit http://dx.doi.org/10.1017/S2045796015001171
Acknowledgements
Ana Villaverde-Hueso, pharmacist at the Rare Diseases Research Institute (Carlos III Institute of Health), carried out the identification and classification of medications.
Financial Support
This work was supported by the Spanish Health Research Fund (Fondo de Investigación Sanitaria) (PI96/0201).
Conflict of Interest
Authors declare no conflict of interest.
Ethical Standards
The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008.
Open access
