Introduction
Depression is a common and disabling disorder in later life (Kok & Reynolds, Reference Kok and Reynolds2017), and has a poorer prognosis compared with depression in younger individuals (Jeuring et al., Reference Jeuring, Stek, Huisman, Oude Voshaar, Naarding, Collard and Comijs2018; Schaakxs et al., Reference Schaakxs, Comijs, Lamers, Kok, Beekman and Penninx2018). Among other factors, this worse prognosis of late-life depression (LLD) might be explained by aging-related physiological changes such as the occurrence of metabolic syndrome (MetS) and inflammation (Dowlati et al., Reference Dowlati, Herrmann, Swardfager, Liu, Sham, Reim and Lanctôt2010; Howren, Lamkin, & Suls, Reference Howren, Lamkin and Suls2009; Koponen, Jokelainen, Keinänen-Kiukaanniemi, Kumpusalo, & Vanhala, Reference Koponen, Jokelainen, Keinänen-Kiukaanniemi, Kumpusalo and Vanhala2008; Köhler et al., Reference Köhler, Freitas, Maes, de Andrade, Liu, Fernandes and Carvalho2017; Marijnissen et al., Reference Marijnissen, Smits, Schoevers, van den Brink, Holewijn, Franke and Oude Voshaar2013; Repousi, Masana, Sanchez-Niubo, Haro, & Tyrovolas, Reference Repousi, Masana, Sanchez-Niubo, Haro and Tyrovolas2018).
MetS and depression have a bidirectional association. MetS predicts the onset and persistence of depression (Marijnissen et al., Reference Marijnissen, Vogelzangs, Mulder, van den Brink, Comijs and Oude Voshaar2017; Vogelzangs et al., Reference Vogelzangs, Beekman, van Reedt Dortland, Schoevers, Giltay, de Jonge and Penninx2014a, Vogelzangs, Comijs, Oude Voshaar, Stek, & Penninx, Reference Vogelzangs, Comijs, Oude Voshaar, Stek and Penninx2014b), and vice versa depression the onset of MetS (Pan et al., Reference Pan, Keum, Okereke, Sun, Kivimaki, Rubin and Hu2012). The association between MetS and depression may be explained by immune-inflammatory alterations, which are present in both conditions (Capuron et al., Reference Capuron, Su, Miller, Bremner, Goldberg, Vogt and Vaccarino2008; Martinac et al., Reference Martinac, Babić, Bevanda, Vasilj, Glibo, Karlović and Jakovljević2017). Meta-analyses consistently demonstrate an association between depression and inflammation, also in the absence of systemic inflammatory diseases (Capuron et al., Reference Capuron, Su, Miller, Bremner, Goldberg, Vogt and Vaccarino2008; Dowlati et al., Reference Dowlati, Herrmann, Swardfager, Liu, Sham, Reim and Lanctôt2010; Howren et al., Reference Howren, Lamkin and Suls2009; Köhler et al., Reference Köhler, Freitas, Maes, de Andrade, Liu, Fernandes and Carvalho2017). A proposed explanation for the association between inflammation and depression is the ‘cytokine hypothesis’, in which depression is considered to be the result of stress-related increased production of pro-inflammatory cytokines (Schiepers, Wichers, & Maes, Reference Schiepers, Wichers and Maes2005; Smith, Reference Smith1991), which eventually leads to less availability of serotonin and over-activation of the hypothalamic, pituitary, and adrenal axis (Maes, Leonard, Myint, Kubera, & Verkerk, Reference Maes, Leonard, Myint, Kubera and Verkerk2011; Wichers & Maes, Reference Wichers and Maes2004). This association, however, is less consistent in older populations than in younger populations (Vogelzangs et al., Reference Vogelzangs, Beekman, van Reedt Dortland, Schoevers, Giltay, de Jonge and Penninx2014a, Reference Vogelzangs, Comijs, Oude Voshaar, Stek and Penninx2014b). Moreover, metabolic dysregulations without inflammation have also been found in depressed older patients (Capuron et al., Reference Capuron, Su, Miller, Bremner, Goldberg, Vogt and Vaccarino2008). Nonetheless, depressed older persons with inflammation tend to have a worse prognosis (Gallagher, Kiss, Lanctot, & Herrmann, Reference Gallagher, Kiss, Lanctot and Herrmann2017) and higher depression symptom levels over time (Au, Smith, Gariépy, & Schmitz, Reference Au, Smith, Gariépy and Schmitz2015) than those without inflammation.
The inconsistent results in older samples might be explained by the etiological heterogeneity among LLD patients. Data-driven cluster analyses applied to MetS and inflammatory markers may identify etiologically distinct subgroups of patients with LLD. Until now, subtyping depression using data-driven clustering techniques has predominantly been based on patients' depressive symptom profiles (Marquand, Wolfers, Mennes, Buitelaar, & Beckmann, Reference Marquand, Wolfers, Mennes, Buitelaar and Beckmann2016; van Loo, de Jonge, Romeijn, Kessler, & Schoevers, Reference van Loo, de Jonge, Romeijn, Kessler and Schoevers2012). In population-based studies, latent class analyses (LCA) have mostly yielded different classes (or subtypes) reflecting different levels of depressive symptom severity (e.g. Mezuk & Kendler, Reference Mezuk and Kendler2012; Ten Have et al. Reference Ten Have, Lamers, Wardenaar, Beekman, de Jonge, van Dorsselaer and de Graaf2016), whereas only a few studies yielded classes with different symptom patterns (e.g. Wanders et al. Reference Wanders, van Loo, Vermunt, Meijer, Hartman, Schoevers and de Jonge2016; Wardenaar, Wanders, Ten Have, de Graaf, & de Jonge, Reference Wardenaar, Wanders, Ten Have, de Graaf and de Jonge2017). In clinical samples of depressed patients, LCAs have mostly yielded classes with different symptom profiles, such as melancholic severe depression and atypical severe depression, in both younger (Lamers et al., Reference Lamers, de Jonge, Nolen, Smit, Zitman, Beekman and Penninx2010; Li et al., Reference Li, Aggen, Shi, Gao, Li, Tao and Kendler2014) and older patients (Veltman et al., Reference Veltman, Lamers, Comijs, de Waal, Stek, van der Mast and Rhebergen2017). Nonetheless, the atypical subtype was only associated with higher levels of inflammatory parameters among younger (Lamers et al., Reference Lamers, Vogelzangs, Merikangas, de Jonge, Beekman and Penninx2013), but not among older patients (Veltman et al., Reference Veltman, Lamers, Comijs, Stek, van der Mast and Rhebergen2018). An explanation might be that (atypical) depression early in life is closely associated with metabolic-inflammatory dysregulation, whereas in later life atypical depression has a more heterogeneous etiology. To our knowledge, the only study that applied LCA on somatic biomarkers was conducted in depressed younger (18–65 years) patients and showed that the overweight subgroup had the highest levels of psychopathology (Beijers et al., Reference Beijers, Wardenaar, Bosker, Lamers, van Grootheest, de Boer and Schoevers2018).
The first objective of the present study was to classify depressed older patients according to their patterns of metabolic and inflammatory dysregulations. We hypothesized to identify a subgroup free of metabolic-inflammatory dysregulation and one or more subgroups with a specific metabolic-inflammatory profile. Our second objective was to compare the prognosis of depression over a 2-year follow-up between identified subgroups. We hypothesized to find a worse prognosis for subgroups characterized by metabolic-inflammatory dysregulation.
Methods
Study design and sample
The study is embedded in the Netherlands Study of Depression in Older persons (NESDO), an ongoing multi-site cohort study designed to examine the course and consequences of depressive disorders in older persons (⩾ 60 years) (Comijs et al., Reference Comijs, van Marwijk, van der Mast, Naarding, Oude Voshaar, Beekman and Smit2011, Reference Comijs, Nieuwesteeg, Kok, van Marwijk, van der Mast, Naarding and Stek2015). In brief, the cohort consists of 378 depressed and 132 non-depressed older persons aged from 60 to 93 years, recruited from mental health care institutes and general practitioners from 2007 until 2010. Participants with a primary diagnosis of dementia, a Mini Mental State Examination (MMSE) score of <18 (out of 30 points) (Folstein, Folstein, & McHugh, Reference Folstein, Folstein and McHugh1975), a primary psychotic or bipolar disorder or insufficient command of the Dutch language were excluded. Inclusion was based on the diagnosis of depression according to the Diagnostic and Statistical Manual of Mental Disorders-IV-R (DSM-IV-R) criteria and was defined as a past 6-months major depressive disorder (MDD) (95%), dysthymic disorder (26.5%) or past-month minor depression (two to four depressive symptoms lasting at least 2 weeks) (5%). The comparison group of non-depressed participants was recruited at the same general practices that recruited the depressed participants and had no history of depressive disorder.
Baseline assessments included written questionnaires, a face-to-face interview, a medical examination, cognitive tests and the collection of fasten blood samples in the morning. Information was gathered about mental health outcomes, demographic characteristics and psychosocial, biological, cognitive and genetic determinants. At 2-year follow-up, all measures open to change were administered again (Comijs et al., Reference Comijs, Nieuwesteeg, Kok, van Marwijk, van der Mast, Naarding and Stek2015). Furthermore, the course of depression was assessed every 6-months using postal questionnaires.
The ethical review boards of the participating institutes approved this study and all participants provided written informed consent.
For the present study, we selected only the depressed patient group (n = 378). A total of 14 patients were excluded because of a fever in the week prior to blood withdrawal, as this may bias the level of inflammatory markers. Of the remaining 364 patients, 273 patients participated in the 2-year follow-up visit.
Assessment of depression
At baseline and 2-year follow-up, the Composite International Diagnostic Interview (CIDI; WHO version 2.1) was used to assess whether criteria were met for a diagnosis of depression or dysthymia, according to the DSM-IV-R criteria. The CIDI is a structured clinical interview with high validity for depressive and anxiety disorders (Kessler et al., Reference Kessler, Birnbaum, Shahly, Bromet, Hwang, McLaughlin and Stein2010; Wittchen et al., Reference Wittchen, Robins, Cottler, Sartorius, Burke and Regier1991). Questions were added to diagnose current minor depression according to the research criteria of the DSM-IV-R (Comijs et al., Reference Comijs, van Marwijk, van der Mast, Naarding, Oude Voshaar, Beekman and Smit2011).
Severity of depression was measured with the well-validated 30-item Inventory of Depressive Symptomatology-Self Report (IDS-SR; Hegeman et al., Reference Hegeman, Wardenaar, Comijs, de Waal, Kok and van der Mast2012; Rush, Gullion, Basco, Jarrett, & Trivedi, Reference Rush, Gullion, Basco, Jarrett and Trivedi1996). The course of depression severity was assessed using the repeated IDS-SR scores taken at baseline and every 6-months until 2-year follow-up (five IDS-SR scores per participant).
Patients who, after 2-year follow-up, no longer met the DSM-IV-TR criteria for depression (either 6-months MDD or dysthymic disorder, or past-month minor depression) were considered remitted from depression.
Biomarkers
Several metabolic and inflammatory biomarkers were used as input variables for the LCA to identify biomarker-based subgroups of patients.
Metabolic dysregulation
Seven included markers of metabolic dysregulation were waist circumference, levels of triglycerides, high-density lipoprotein (HDL), low-density lipoprotein (LDL), glucose level, as well as systolic and diastolic blood pressure. Waist circumference was assessed in centimeters halfway between lower rib margin and the iliac crest following normal expiration upon light clothing. Triglyceride, HDL cholesterol, LDL cholesterol, total cholesterol and glucose levels were measured in morning fasting blood samples using routine standardized laboratorial methods. Blood pressure was averaged over two readings measured in a supine position using an electric Omron sphygmomanometer. All variables were considered as continuous measures in the LCA.
Body mass index (BMI) and total cholesterol were not included as these markers were strongly (r > 0.60) correlated with other markers (r = 0.79 for BMI and waist circumference; r = 0.92 for total- and LDL-cholesterol). The correlations between included variables varied between −0.42 and 0.59.
Inflammation
Four inflammation markers, i.e. high-sensitive C-reactive protein (hsCRP), interleukin-6 (IL-6), Growth Differentiation Factor-15 (GDF-15) and neutrophil gelatinase-associated lipocalin (NGAL), were assessed in fasting blood samples obtained in the morning between 8 and 9 a.m. after an overnight fasting, and kept frozen at −80 °C. These markers have previously been observed to be associated with depression (Dowlati et al., Reference Dowlati, Herrmann, Swardfager, Liu, Sham, Reim and Lanctôt2010; Howren et al., Reference Howren, Lamkin and Suls2009; Köhler et al., Reference Köhler, Freitas, Maes, de Andrade, Liu, Fernandes and Carvalho2017; Naudé et al., Reference Naudé, Eisel, Comijs, Groenewold, De Deyn, Bosker and Oude Voshaar2013; Teunissen, Durieux-Lu, Blankenstein, Oude Voshaar, & Comijs, Reference Teunissen, Durieux-Lu, Blankenstein, Oude Voshaar and Comijs2016). All markers were assessed using well-validated assays. hsCRP levels were measured in duplicate by using an immunoturbidimetric assay, IL-6 and NGAL by using an enzyme-linked immunosorbent assay and GDF-15 using an automated assay on Abbott Architect. Intra- and inter-assay coefficients of variation were 2% and 2% for hsCRP, 8% and 12% for IL-6 and 3% and 5% for NGAL. The intra-assay coefficient of variation for GDF-15 was 3.8%.
Patient characteristics
Several additional measurements were performed to compare any identified LCA classes.
Metabolic dysregulation
In addition to the MetS variables used in the LCA, we assessed the presence of MetS (yes/no), diabetes mellitus (yes/no), BMI and total cholesterol at baseline. The MetS was defined according to the National Cholesterol Education Programme (NCEP)-Adult Treatment Panel III guidelines (Grundy et al., Reference Grundy, Cleeman, Merz, Brewer, Clark and Hunninghake2004).
Clinical characteristics
Key clinical characteristics of interest were: age of onset of depressive disorder (CIDI question), depression severity based on IDS-SR sum score, melancholic and/or atypical features [presence of DSM-defined melancholic features were assessed based on the IDS-items according to Kahn's algorithm (Khan et al., Reference Khan, Carrithers, Preskorn, Lear, Wisniewski, John Rush and Fava2006) and DSM-defined atypical features based on IDS-items according to Novick's algorithm (Novick et al., Reference Novick, Stewart, Wisniewski, Cook, Manev and Nierenberg2005)]. Furthermore, mood, motivation and somatic-affective symptom dimensions according to the IDS subscale scores (Hegeman et al., Reference Hegeman, Wardenaar, Comijs, de Waal, Kok and van der Mast2012), diagnosis of dysthymia according to DSM-IV-R criteria, a current anxiety disorder (CIDI), level of anxiety based on the Beck Anxiety Inventory (BAI) sum scores and the use of antidepressants (yes/no, as well as differentiated into categories SSRI/TCA/MAO/other).
Demographics
Standard demographics included age, sex, educational level (in years) and marital status (partner yes/no).
Life-style and physical health indicators
Indicators of life style and physical health included smoking (yes/no), use of alcohol, physical activity, global cognitive functioning and number of chronic diseases.
The number of alcoholic drinks per day was based on the first two items of the Alcohol Use Disorders Identification Test (AUDIT) (Aalto, Alho, Halme, & Seppä, Reference Aalto, Alho, Halme and Seppä2011; Babor, Kranzler, & Lauerman, Reference Babor, Kranzler and Lauerman1989). Physical activity was measured with the last-seven-days short-form (8-items) of the International Physical Activity Questionnaire (IPAQ) and classified as minimal, moderate or high (Craig et al., Reference Craig, Marshall, Sjöström, Bauman, Booth, Ainsworth and Oja2003). Global cognitive functioning was assessed by the MMSE (Folstein et al., Reference Folstein, Folstein and McHugh1975). The number of chronic diseases was assessed by self-report questions with high accuracy compared to general practitioner information (Kriegsman, Penninx, van Eijk, Boeke, & Deeg, Reference Kriegsman, Penninx, van Eijk, Boeke and Deeg1996). We asked for cardiac disease (including myocardial infarction), peripheral atherosclerosis, stroke, diabetes mellitus, chronic obstructive pulmonary disease (asthma, chronic bronchitis or pulmonary emphysema), arthritis (rheumatoid arthritis or osteoarthritis), cancer or any other disease.
Statistical analysis
First, LCAs were performed to identify biomarker-based patient subgroups. Models with increasing numbers of classes were fit to the data and the most optimal model was selected based on the lowest values for the Bayesian information criterion (BIC) and Akaike information criterion (AIC). In addition, model interpretability (i.e. sufficient qualitative differentiation between classes) was considered when choosing the number of classes. The LCAs were run using a robust maximum likelihood estimator (MLR) in Mplus 5.0 (Muthén & Muthén, Reference Muthén and Muthén2007). To avoid a model solution at a local maximum, each model was run with multiple random starts. After identification of the optimal model, each patient was allocated to a class based on their highest posterior class probability.
After class-assignment, the subgroups were compared with respect to demographic and clinical characteristics at baseline. Next, we examined the association of class-membership at baseline (independent variable) with remission of LLD at 2-year follow-up (dependent variable) using logistic regression. Finally, the course of depressive symptoms over time as assessed with the 6-monthly IDS-SR was examined in the different classes using mixed model analysis (Twisk, de Boer, de Vente, & Heymans, Reference Twisk, de Boer, de Vente and Heymans2013). Models with random coefficients for intercept and/or slope per subject were compared, and the following elements were subsequently tested to determine the best fitting model, using the likelihood ratio test: (1) a general linear change in depression severity in all classes, (2) a general quadratic – i.e. parabolic – change over time, (3) stable differences in depression severity between the classes and (4) differences in depression course between the classes. All models were adjusted for covariates known to be associated with the course of depression, i.e. age, sex, level of education, cognitive functioning (MMSE), number of chronic diseases, antidepressant drug use and lifestyle characteristics (i.e. use of alcohol, smoking and physical activity).
These analyses were conducted in SPSS, version 25. We considered p values less than 0.05 as significant.
Results
Characteristics of the study population
The 364 depressed older patients had a mean age of 70.7 (standard deviation 7.3; range 60–90) years and 66.2% were females (Tables 1 and 2).
Table 1. Metabolic and inflammatory differences between subgroups at baseline
a Total sample consisting of six classes.
b p value from ANOVA or χ2.
c Hochberg post hoc analysis was chosen because of differences in class-size.
d Log-transformation was performed, the values listed are the retransformed mean and s.d. of the log values.
e Parameters not included in the LCA.
Table 2. Demographic and clinical differences between subgroups at baseline
Class A = healthy subgroup, class B = metabolic and inflammatory dysregulation subgroup, class C = severe inflammation, class D = mild inflammation.
a Total sample consisting of six classes.
Latent class analyses (LCA)
The best fitting LCA model according to the AIC and BIC was the six-class model, which also had good interpretability. Two classes were not used in the subsequent analyses, since these classes were too small and based on an extreme value of only one specific parameter, i.e. hsCRP (class 5, three patients) and IL-6 (class 6, two patients).
Furthermore, a sensitivity analysis, performed by repeating the LCA restricted to the subsample of patients with a current (6-month) MDD (n = 345), resulted in a comparable class solution.
Table 1 shows the metabolic and inflammation parameters of the four selected classes (A through D). Class A (n = 181, 49.7%) consisted of a relatively healthy subgroup of depressed patients with low scores across MetS and inflammation markers and is referred to as ‘Healthy’ (healthy solely referring to their metabolic-inflammatory status). This depressed subgroup was characterized by comparatively lower measures of waist circumference, triglyceride, glucose and BMI, a lower frequency of diabetes and the MetS, and higher levels of HDL cholesterol whereas inflammation markers were not elevated.
In class B (n = 137, 37.6%), 75.2% of the patients had MetS, reflected by higher waist circumference, triglyceride levels and systolic and diastolic blood pressure; as well as slightly higher levels of inflammatory markers. Therefore, this subgroup is referred to as ‘Metabolic and inflammatory dysregulation’.
Class C (n = 27, 7.4%) was characterized by higher levels of inflammation markers, mostly hsCRP and IL-6. This subgroup is referred to as ‘Severe inflammation’.
Similarly, class D (n = 14, 3.8%) showed increased inflammation levels, although less severe than the previous group, and was characterized by elevated measures of GDF-15 and NGAL. This subgroup is referred to as ‘Mild inflammation’.
Characteristics of LCA identified subgroups
Table 2 shows the demographical and clinical characteristics of the four depressed subgroups (A through D, n = 359). The healthy subgroup (class A) predominantly consisted of physically active females with few chronic diseases and less severe levels of depression and anxiety than the other classes. In contrast, the metabolic and inflammatory dysregulation subgroup (class B) was less physically active, suffered from more chronic diseases and had more severe levels of depression and anxiety than the other subgroups. The higher depression severity in this class was based on higher scores on the mood and motivation items compared to the other classes. There was no difference in frequency of atypical or melancholic depression subtype between the classes. However, this metabolic and inflammatory dysregulation subgroup more frequently suffered from dysthymia. The severe inflammation subgroup (class C) did not distinguish itself from the other subgroups on demographic or clinical parameters. The mild inflammation subgroup (class D) was the oldest and less physically active subgroup with a higher age of depression onset than the metabolic and inflammatory dysregulation subgroup (class B) and the severe inflammation (class C) subgroup, and was the only predominantly male group.
Remission of late-life depression at 2-year follow-up
Follow-up data were missing for 86/359 patients (24.0%). In total 23/86 patients were deceased at 2-year follow-up, the rest was lost in follow-up because of different reasons. The dropout-rate differed significantly between the depressed subgroups (χ2 = 10.0, df = 3, p = 0.018). The mild inflammation subgroup (class D) had the highest dropout-rate [7/14 (50.0%), of which three deceased), followed by the severe inflammation subgroup [class C, 8/27 (29.6%), of which three deceased], the metabolic and inflammatory dysregulation subgroup [class B, 38/137 (27.7%), of which 11 deceased] and the healthy subgroup [class A, 33/181 (18.2%), of which six deceased].
Patients who no longer participated at 2-year follow-up differed from the other patients in several baseline metabolic, inflammation, demographic and clinical parameters. Drop-outs had a higher waist circumference [96.4 (13.7) v. 92.6 (12.9); t = −2.4, df = 356, p = 0.019), a higher BMI [27.2 (4.7) v. 26.0 (4.3); t = −2.2, df = 357, p = 0.026], a higher level of hsCRP [2.31 (2.87) v. 1.71 (3.12); t = −2.1, df = 346, p = 0.037], a higher level of GDF-15 [1080.19 (1.71) v. 856.25 (1.56); t = −3.9, df = 349, p ⩽ 0.001] and were less physically active (low 45.2%, moderate 27.4%, high 27.4% v. low 26.5%, moderate 41.3%, high 32.2%; χ2 = 10.9, df = 2, p = 0.004) than patients whom did not drop-out.
At 2-year follow-up 141/273 (51.6%) patients did no longer meet the DSM-IV-TR criteria for depression (either 6-month MDD or dysthymic disorder, or past-month minor depression). There was a significant difference in remission-rates between the four depressed subgroups (χ2 = 12.2, df = 3, p = 0.007). The healthy subgroup (class A) had the highest remission rate (90/148, 60.8%), followed by the mild inflammation subgroup (class D, 4/7, 57.1%) and the severe inflammation subgroup (class C, 9/19, 47.4%). The metabolic and inflammatory dysregulation subgroup had the lowest remission rate (class B, 38/99, 38.4%).
Logistic regression confirmed that class membership predicted 2-year remission (Wald = 11.9, df = 3, p = 0.008), but the association was no longer significant after adjustment for covariates (Wald = 6.3, df = 3, p = 0.098), although the difference between the two largest subgroups remained. Patients in the metabolic and inflammatory dysregulation subgroup (class B) had lower odds of achieving remission compared to patients in the healthy subgroup (class A), in both the unadjusted [odds ratio (OR) 0.40 (95% confidence interval (CI) 0.24–0.67), p = 0.001] and adjusted analysis [OR 0.49 (95% CI 0.26–0.91), p = 0.025].
Course of late-life depression over 2 years
As shown in Table 3, the successive inclusion of a linear effect of time, a quadratic effect of time, and a main effect of class improved the model fit of the mixed model exploring the severity of depressive symptoms over time as indicated by the likelihood ratio tests at each step. A further inclusion of either an interaction effect of class with time (step 4A) or with time square (step 4B), did not improve the model fit any further. This means that no significant differences were found between the classes in depression course. At all steps, a random intercept and random slope model proved to be the best fitting model, as indicated by the likelihood ratio test.
Table 3. Model selection and optimal mixed model analysis of course of depression severity over 2-year follow-up
a Adjusted for age, sex, years of education, smoking, alcohol use, physical activity, somatic comorbidity, global cognitive functioning and use of antidepressants.
b No other significant differences between classes (p < 0.05) were found.
The optimal model showed that for all patient classes, depression severity is best estimated to decline linearly from baseline by (after adjustment for covariates) 6.0 points per year on the IDS-SR (95% CI 4.1–7.8), but this decline tapers off at a quadratic rate of 1.8 (1.0–2.7) points per squared number of years since baseline, resulting in an expected net decline of 4.3 points over the first year and 0.6 over the second year. Furthermore, the metabolic and inflammatory dysregulation subgroup (class B) consistently scored 4.8 (2.4–7.2) points above the healthy subgroup (class A), from baseline on, and this difference is only partly explained by the factors known to influence the course of depression which were included as covariates. After adjustment for these covariates, the stable difference between the two subgroups remained, but was reduced to 3.1 (0.5–5.6).
Because almost all change in depression severity occurred in the first year of follow-up, we also analyzed the course of depression over this period separately. Here too the models including both a random intercept and slope proved to be the best fitting models, at all steps. All consecutive additions of (1) a linear effect of time, (2) a quadratic effect of time, (3) a main effect of class and (4) an interaction effect of class with either time (4A) or time square (4B) improved the fit of the model (see Table 4). Because in step 4 the addition of the interaction between class and time improved the fit of the model more than the addition of interaction between class and time square (as indicated by the reduction in log likelihood), this model was chosen as best fitting model for the course of depression severity over the first follow-up year.
Table 4. Model selection and optimal mixed model analysis of course of depression severity over the first year of follow-up
a Adjusted for age, sex, years of education, smoking, alcohol use, physical activity, somatic comorbidity, global cognitive functioning and use of antidepressants.
b Difference between class A and class B in unadjusted analysis (t = 3.163; p = 0.002), not in adjusted analysis (t = 1.590; p = 0.113).
c No other significant differences between classes (p < 0.05) were found.
Table 4 and Fig. 1 show that for the severe inflammation subgroup (class C, reference group in Table 4) it can be seen, that – after adjustment for covariates – a sharp linear decline from baseline in depression severity occurred (at a rate of 6.9 points per year on the IDS-SR; 95% CI 0.9–12.8), but this decline soon leveled off at a quadratic rate of 7.8 (4.0–11.6) points per squared number of years since baseline, in such a way that after about five months a maximum decline of 1.5 points from the baseline level was reached, after which the severity rose again to above baseline level after 1 year. The course of depression in the severe inflammation subgroup (class C) did not differ significantly from that seen in the mild inflammation subgroup (class D). In the metabolic and inflammatory dysregulation (class B) and healthy (class A) subgroups, the decline in depression severity over the first follow-up year was significantly larger than in the severe inflammation subgroup (class C), with a difference of 5.1 (<−0.1 to 10.1) and 6.5 (1.6–11.5) points per year, respectively. Finally, a difference was seen in the unadjusted analysis between the metabolic and inflammatory dysregulation (class B) and healthy (class A) subgroups, with the latter scoring 4.4 (1.7–7.2) points lower on average at baseline, but this difference was no longer significant after correction for covariates.
Fig. 1. Estimated course of depression severity over the first year (adjusted for age, sex, years of education, smoking, alcohol use, physical activity, somatic comorbidity, global cognitive functioning and use of antidepressants). Class A: ‘Healthy’, class B: ‘Metabolic and inflammatory dysregulation’, class C: ‘Severe inflammation’, class D: ‘Mild inflammation’.
Discussion
Main findings
Among 364 depressed older patients, LCA identified six classes according to metabolic-inflammatory parameters. The largest class consisted of 181 (49.7%) depressed older persons showing minimal metabolic-immune dysregulation. The other five groups did show deviations from this health pattern in their metabolic and/or inflammatory marker levels. The largest effected subgroup consisted of patients with an overall, but relatively mild level of metabolic-inflammatory dysregulation (n = 137, 37.6%). Smaller subgroups were characterized by increased levels of hsCRP and IL-6 (severe inflammation, n = 27, 7.4%), by increased levels of GDF-15 and NGAL (mild inflammation, n = 14, 3.8%), an isolated extremely elevated level of hsCRP (n = 3), or an isolated extremely elevated level of GDF-15 (n = 2). The three larger subgroups with low grade inflammation showed the highest levels of psychopathology compared to the healthy subgroup. This was most pronounced for the metabolic inflammatory dysregulation subgroup, which had the highest level of depression severity. In line with our hypothesis, patients in the subgroup with both metabolic and inflammatory dysregulation consistently scored higher on depression severity over the follow-up period than patients without metabolic-inflammatory dysregulation, and were less often remitted at 2 year follow-up. Patients with severe inflammation showed the least decline in depression severity in the first year.
Metabolic and inflammatory characterized depression
To our knowledge, only one study applied LCA to identify depression subgroups within depressed younger (18–65 years) patients using somatic biomarkers, including metabolic and inflammatory markers. Three subgroups were identified that were labeled as ‘lean’, ‘average’ and ‘overweight’ based on their average BMI value. However, this labeling was somewhat arbitrary, as higher BMI values in that study were also associated with abnormalities in the other somatic parameters. In line with the associations between MetS, depression (severity) and the role of inflammation, the overweight subgroup had the highest levels of psychopathology (Beijers et al., Reference Beijers, Wardenaar, Bosker, Lamers, van Grootheest, de Boer and Schoevers2018). In contrast to our findings, however, the course of psychopathology did not differ between these groups (Beijers et al., Reference Beijers, Wardenaar, Bosker, Lamers, van Grootheest, de Boer and Schoevers2018).
While our identified subgroup with both metabolic and inflammatory dysregulation scored significantly worse on all metabolic and inflammatory parameters compared to the healthy subgroup (except for LDL- and total cholesterol), the average level of inflammatory markers did not exceed the cut-off for low-grade inflammation applied in cardiovascular research (e.g. a hsCRP ⩾3; Pearson et al., Reference Pearson, Mensah, Alexander, Anderson, Cannon and Criqui2003).
Additionally, the metabolic and inflammatory dysregulation subgroup more often suffered from dysthymia and had higher levels of anxiety at baseline, and consistently showed significantly higher depression severity levels over the follow-up period compared to the healthy subgroup. Thus, the small elevation in inflammation in this subgroup has great clinical consequences. Possibly, pathways involved in depression are highly sensitive to small changes in inflammation markers resulting in depression.
Although at younger age metabolic-inflammatory dysregulation seems to be related to atypical depression and the somatic-affective symptoms of depression, we did not find these associations in later-life. First, we did not find an association between depression subtypes (atypical/ melancholic) and a specific class. In adult depressed patients an association between atypical depression and metabolic-inflammatory dysregulation was found (Lamers et al., Reference Lamers, Vogelzangs, Merikangas, de Jonge, Beekman and Penninx2013; Milaneschi et al., Reference Milaneschi, Lamers, Peyrot, Baune, Breen and Dehghan2017; Vogelzangs et al., Reference Vogelzangs, Beekman, van Reedt Dortland, Schoevers, Giltay, de Jonge and Penninx2014a, Reference Vogelzangs, Comijs, Oude Voshaar, Stek and Penninx2014b), nevertheless previous findings in older depressed patients were not consistent (Veltman et al., Reference Veltman, Lamers, Comijs, Stek, van der Mast and Rhebergen2018; Vogelzangs et al., Reference Vogelzangs, Beekman, van Reedt Dortland, Schoevers, Giltay, de Jonge and Penninx2014a, Reference Vogelzangs, Comijs, Oude Voshaar, Stek and Penninx2014b). Therefore, distinction between atypical and melancholic depression in later life does not seem to be driven by metabolic-inflammatory dysregulation. Furthermore, in contrast to previous studies, higher depression severity levels were not associated with somatic depression symptoms (Capuron et al., Reference Capuron, Su, Miller, Bremner, Goldberg, Vogt and Vaccarino2008; Dantzer, O'Connor, Freund, Johnson, & Kelley, Reference Dantzer, O'Connor, Freund, Johnson and Kelley2008; Duivis, Vogelzangs, Kupper, de Jonge, & Penninx, Reference Duivis, Vogelzangs, Kupper, de Jonge and Penninx2013; Marijnissen et al., Reference Marijnissen, Smits, Schoevers, van den Brink, Holewijn, Franke and Oude Voshaar2013).
Inflammatory depression
The severe inflammation subgroup was characterized by higher levels of hsCRP and IL-6 and the mild inflammation subgroup was characterized by elevated levels of NGAL and GDF-15. Patients in both subgroups were more obese and more often suffered from MetS compared to the healthy subgroup. Nonetheless, the level of inflammatory dysregulation in both subgroups was higher than can be expected based on solely the metabolic dysregulation. Therefore, the specific elevated levels of hsCRP and IL-6 in the severe inflammation subgroup may be due to an underlying primary inflammatory disease. Likewise, the elevated levels of NGAL and GDF-15 in the mild inflammation subgroup may represent the aging process, which can be accompanied by inflammatory dysregulation also known as ‘inflammaging’ (Franceschi et al., Reference Franceschi, Zaikin, Gordleeva, Ivanchenko, Bonifazi, Storci and Bonafè2018), as this subgroup was the oldest, most physically inactive subgroup with the highest age of onset.
Depression prognosis according to metabolic and inflammatory status
The healthy subgroup had a higher rate of depression remission at 2-year follow-up compared to the metabolic and inflammatory dysregulation subgroup. Furthermore, only the healthy subgroup continued to show a decrease in depression symptoms after the first half-year follow-up, whereas the other subgroups showed a deterioration of depression symptoms in the second half-year of follow-up. An explanation might be that patients were included after they had enlisted themselves for the start of depression treatment. Therefore, depression treatment might lose its effectiveness after the first half-year in patients with metabolic and/or inflammatory dysregulation. Chronic inflammation might thus represent a continuous vulnerability that averts long-term stabilization. A meta-analysis found that treatment targeting the inflammatory dysregulation in depression as an add-on to antidepressant treatment has beneficial effects on the depression of patients with MDD (Köhler-Forsberg et al., Reference Köhler-Forsberg, Nicolaisen Lydholm, Hjorthoj, Nordentoft, Mors and Benros2019). Examples of this treatment are non-steroidal anti-inflammatory drugs (NSAIDs) and statins. Nonetheless further research is necessary to understand the mechanism responsible for this anti-depressant effect. It is unclear whether the positive effects of anti-inflammatory treatment are caused by direct involvement in the depression cascade, or that the effect may be mediated by effects on somatic disease. Furthermore, effects were small, which suggests that only specific subgroups might benefit from augmentation with anti-inflammatory drugs. The subgroups with inflammatory dysregulation identified by our study might therefore be of particular interest for further research on anti-inflammatory therapy in the treatment of depression. These findings are in correspondence with previous studies in which both metabolic dysregulation and inflammation were found to be associated with worse depression outcomes in younger patients and cardiac patients (Capuron et al., Reference Capuron, Su, Miller, Bremner, Goldberg, Vogt and Vaccarino2008; Duivis et al., Reference Duivis, Vogelzangs, Kupper, de Jonge and Penninx2013; Howren et al., Reference Howren, Lamkin and Suls2009; Marazziti, Rutigliano, Baroni, Landi, & Dell'Osso, Reference Marazziti, Rutigliano, Baroni, Landi and Dell'Osso2014).
Strengths and limitations
Strengths of this prospective study are the substantial sample size with availability of several clinical, and immune and metabolic parameters for LLD patients. Nonetheless, the severe and mild inflammation subgroups were relatively small subgroups with higher missings at 2-year follow-up. This could have led to attrition bias and limited statistical power to identify small but possibly meaningful differences between the classes. Moreover, the two extremely small classes, based on an extreme value of only one specific parameter, could not be further explored reliably. These subgroups might represent relatively rare cases with a specific underlying condition associated with a specific inflammatory profile. Finally, the subgroups may differ in the treatment received during the follow-up period. The data collected in the NESDO study did not allow us to control the analyses for these potential differences. However, at baseline the four subgroups did not differ in the medication treatment received.
Conclusion
Nearly half of the patients with LLD showed metabolic-inflammatory dysregulation that was associated with more severe depression and worse depression prognosis. Future studies should examine whether the largest group with an overall, but relatively mild level of dysregulation would benefit from augmentation with anti-inflammatory drugs like NSAIDs or statins.
Author contributions
All authors had substantial contributions to conception and design or analysis and interpretation of the data. All authors had substantial contributions to drafting the article and revising it critically for important intellectual content. All authors have approved the final version.
Financial support
The infrastructure for NESDO is funded through the Fonds NutsOhra, Stichting tot Steun VCVGZ, NARSAD The Brain and Behaviour Research Fund, and the participating universities and mental health care organizations (VU University Medical Center, Leiden University Medical Center, University Medical Center Groningen, Radboud University Nijmegen Medical Center, and GGZ inGeest, GGNet, GGZ Nijmegen, and Parnassia).
Conflict of interest
The authors have no financial of any other kind of personal conflicts with this paper.
Sponsor's role
The authors maintained full independence in the conduct of this work. The sponsors had no role in design, methods, subject recruitment, data collections, analysis or preparation of the manuscript.
Impact statement
We certify that this work is novel clinical research. The potential impact of this research on clinical care and health policy includes the following: as nearly half of the patients with late-life depressions suffer from metabolic-inflammatory dysregulation which is also associated with more severe depression and worse prognosis, a metabolic-inflammatory targeted treatment may benefit depressed older patients.
Open access
