Next Article in Journal
Documenting Pharmacogenomic Test Results in Electronic Health Records: Practical Considerations for Primary Care Teams
Next Article in Special Issue
The Contribution of “Individual Participant Data” Meta-Analyses of Psychotherapies for Depression to the Development of Personalized Treatments: A Systematic Review
Previous Article in Journal
The Potential of Gamma Secretase as a Therapeutic Target for Cardiac Diseases
Previous Article in Special Issue
Machine Learning-Based Behavioral Diagnostic Tools for Depression: Advances, Challenges, and Future Directions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Patient Stratification Approach to Identifying the Likelihood of Continued Chronic Depression and Relapse Following Treatment for Depression

1
Centre for Outcomes Research and Effectiveness (CORE), Research Department of Clinical, Educational & Health Psychology, University College London, 1-19 Torrington Place, London WC1E 7HB, UK
2
Department of Psychiatry, University of California, Los Angeles, CA 90095, USA
3
Statistical Science, University College London, 1-19 Torrington Place, London WC1E 7HB, UK
4
Department of Psychology, School of Arts and Sciences, 425 S. University Avenue, Philadelphia, PA 19104, USA
5
Centre for Academic Mental Health, Population Health Sciences, Bristol Medical School, University of Bristol, Oakfield House, Bristol BS8 2BN, UK
6
Centre for Academic Primary Care, Population Health Sciences, Bristol Medical School, University of Bristol, Canynge Hall, Bristol BS8 2PS, UK
7
Department of Health Sciences, University of York, ARRC Building, Heslington, York YO10 5DD, UK
8
Department of Psychology, Vanderbilt University, Nashville, TN 37235, USA
9
Primary Care Research Centre, Primary Care, Population Sciences and Medical Education, Faculty of Medicine, University of Southampton, Aldermoor Health Centre, Southampton SO16 5ST, UK
10
Department of Psychology, University of Exeter, Sir Henry Wellcome Building for Mood Disorders Research, Perry Road, Exeter EX4 4QG, UK
11
Institute of Health Research, University of Exeter College of Medicine and Health, Exeter EX1 2LU, UK
12
Department of Health and Caring Sciences, Western Norway University of Applied Sciences, Inndalsveien 28, 5063 Bergen, Norway
13
Division of Psychiatry, University College London, Maple House, London W1T 7NF, UK
14
Camden & Islington NHS Foundation Trust, St. Pancras Hospital, 4 St. Pancras Way, London NW1 0PE, UK
15
iCope—Camden & Islington Psychological Therapies Services—Camden & Islington NHS Foundation Trust, St. Pancras Hospital, 4 St. Pancras Way, London NW1 0PE, UK
*
Author to whom correspondence should be addressed.
J. Pers. Med. 2021, 11(12), 1295; https://doi.org/10.3390/jpm11121295
Submission received: 17 September 2021 / Revised: 15 November 2021 / Accepted: 17 November 2021 / Published: 4 December 2021

Abstract

:
Background: Subgrouping methods have the potential to support treatment decision making for patients with depression. Such approaches have not been used to study the continued course of depression or likelihood of relapse following treatment. Method: Data from individual participants of seven randomised controlled trials were analysed. Latent profile analysis was used to identify subgroups based on baseline characteristics. Associations between profiles and odds of both continued chronic depression and relapse up to one year post-treatment were explored. Differences in outcomes were investigated within profiles for those treated with antidepressants, psychological therapy, and usual care. Results: Seven profiles were identified; profiles with higher symptom severity and long durations of both anxiety and depression at baseline were at higher risk of relapse and of chronic depression. Members of profile five (likely long durations of depression and anxiety, moderately-severe symptoms, and past antidepressant use) appeared to have better outcomes with psychological therapies: antidepressants vs. psychological therapies (OR (95% CI) for relapse = 2.92 (1.24–6.87), chronic course = 2.27 (1.27–4.06)) and usual care vs. psychological therapies (relapse = 2.51 (1.16–5.40), chronic course = 1.98 (1.16–3.37)). Conclusions: Profiles at greater risk of poor outcomes could benefit from more intensive treatment and frequent monitoring. Patients in profile five may benefit more from psychological therapies than other treatments.

1. Introduction

Depression is amongst the most burdensome diseases across the globe [1]. It is highly prevalent, affecting approximately 320 million people annually [2,3], results in significant impairment for prolonged periods of time [4], and is typically thought to follow a ‘relapsing-remitting’ course, with multiple episodes throughout life [5]. Much of our understanding regarding the course of depression comes from general population studies of adults [6,7], but clinical samples are needed to ensure greater relevance to healthcare settings. This is especially important as the course of depression is quite different for those that seek treatment compared to those that do not, with the rate of relapse much higher in clinical samples [8,9].
Treatments for depression including continued use of antidepressant medications can be effective in reducing the risk of relapse and of chronic depression [10]. Although many patients do not want long-term treatment with antidepressants [11], there is an increasing trend for them to be maintained on antidepressants for two years or more [12,13], and the use of antidepressants has greatly increased over the last three decades [14]. Although pharmacogenetics may hold promise for personalizing antidepressant treatments [15,16], it is becoming increasingly common that antidepressants are prescribed indefinitely for all patients with a history of previous relapses, to mitigate the risks of chronic illness and relapse [11,12,17,18]. This is despite evidence of the waning effectiveness of medications over time [19,20,21,22], and even increased risk of relapse and other harmful outcomes for those remaining on medications after several years [17,23]. Some psychological therapies offer protection against relapse which is equivalent to that offered by continuation of antidepressants at six months post-treatment [24,25,26]. Psychological therapies can also be used after the acute phase of treatment, i.e., when a patient reaches remission, these can have a preventive effect [27,28], reducing the risk of relapse by approximately 15–29% up to two years post-treatment [29]. Further, combining antidepressants and psychotherapies or augmenting one with the other when remission has not been achieved with monotherapy improves the chances of avoiding relapse [27,30,31]. However, the benefits of such treatments are not universal, the costs of either treatment can be prohibitive, and psychological treatments in particular are not readily available throughout the world [29,32]. Coupled with the prevalence of depression, this makes it impractical to offer prophylactic treatment to all patients with depression. As such, identifying subgroups of patients for whom further treatment might reduce the likelihood of continued chronic depression or the risk of relapse could have important implications in terms of clinical outcomes and reduced healthcare costs.
Identifying subgroups of patients by better understanding the complex relationships between pre-treatment patient characteristics can allow for the development of ‘clinical phenotypes’ [33], informing prognosis and supporting treatment decision making. One way of doing so is via latent profile analysis (LPA), which is part of the latent variable mixture modelling family of analytic approaches. Based on probability theory, these models identify homogeneous subpopulations by identifying stratified groups of patients [34,35] based on the relationships between available indicator variables (both continuous and categorical variables can be used in the LPA). These methods have been applied to a range of clinical presentations, including in Alzheimer’s disease [36], Crohn’s disease [37], and depression and anxiety [38], where they have also been used to predict likely treatment outcomes for people receiving psychological therapy in routine care settings [39,40]. However, such methods have not been applied to consider the longer-term prognosis of patients after the end of acute-phase treatment for depression, including continued chronic depression or relapse. A number of studies have attempted to develop models to predict relapses to depression [41,42,43]. However, these have been hampered by methodological problems, including small sample sizes, inappropriate handling of missing data, and a lack of validation, and they have largely failed to make accurate predictions [44]. Other studies to develop more accurate models are planned [45], but no findings are available yet, and the clinical utility of any such models are yet to be established. Patient stratification techniques such as LPA offer an alternative means of predicting outcomes which may have greater face-validity for clinicians and patients because the characteristics of each profile can be described in more accessible terms [40]. This is often not possible with regression based or supervised machine learning models because the reasons for any individual being classified as having a particular degree of risk for the outcome of interest may be hard to determine and often harder still to explain when there are multiple predictor variables [46,47]. Furthermore, variable-centred approaches such as regression modelling require a different model, and therefore a different set of associations, by outcomes, whereas patient stratification approaches, following the identification of profiles, can allow exploration of associations between profiles and multiple different outcomes [48]. The application of latent profiling methods to stratify patients and identify individual risks of poor longer-term prognosis might therefore have important clinical utility.
The aims of this study were to identify the stratified profiles of adults treated for depression in primary care settings, based on pre-treatment characteristics, and explore differences in the course of illness and risk of relapse between patients in these profiles. The associations between the types of treatment received on outcomes within each profile were also explored.

2. Materials and Methods

A pre-registered protocol for this study and analytic plan is available through the Open Science Framework (DOI 10.17605/OSF.IO/Y52HG) (https://osf.io/y52hg/ accessed on 29 July 2021).

2.1. Participants

The Depression in General Practice (Dep-GP) individual patient dataset [49] was used to derive the sample for the current analysis. The overall aims of the Dep-GP project are to identify prognostic factors associated with outcomes from treatment for depression in primary care. The dataset contains individual patient data from the participants of 12 randomised controlled trials (RCTs) of treatments for depression. Information regarding the dataset and protocol for the methods of identifying relevant studies for the Dep-GP dataset are available elsewhere [49,50]. Importantly, all trials included in the Dep-GP dataset were of adults with depression who had been recruited from general practices and used the most common comprehensive diagnostic and screening measure of depression and anxiety disorders used in primary care RCTs—the Revised Clinical Interview Schedule (CIS-R) [51].

Inclusion Criteria

For this analysis, individual patient data were sought from RCTs of adults with depression recruited in primary care included in the Dep-GP database, that: used the CIS-R at baseline, had an endpoint at 3 to 4 months post-baseline, and collected data on depressive symptoms at follow-up periods after the 3 to 4 month post-baseline endpoint. Seven studies in the Dep-GP database met these criteria (presented in Table 1).

2.2. Measures

The relevant measures in the seven included studies were:
CIS-R [51]: establishes the nature and severity of common mental disorder symptoms and provides primary and secondary diagnoses according to ICD-10 criteria. CIS-R is made up of 14 symptom subsections. Five subscales contribute to a depression score, covering core features of depression, depressive thoughts, fatigue, concentration and forgetfulness, and sleep; the remaining subscales cover anxiety disorders and related symptoms (generalized anxiety, worry, irritability, somatic anxiety/hypochondriasis, panic, obsessions, compulsions, phobias (split into agoraphobia, social phobia, and specific phobias). Each section includes information on the duration of the problems assessed, and there is an additional subsection that covers general health and disabilities, none of which contributes to either the depression or anxiety total scores. The total scores are calculated by summing the scores of the contributing individual subscales.
Beck Depression Inventory (BDI-II) [52]: This 21-item assessment is used to measure depressive symptoms; each item is scored 0–3 with a maximum score obtainable of 63. A cut-off of ≥14 is used to indicate significant symptoms of depression and ≤10 is used for remission.
Patient Health Questionnaire (PHQ-9) [53]: This is a 9-item depression screening measure. Items are scored 0–3; a cut-off of ≥10 is used to indicate “caseness”, that is symptoms of depression likely to be commensurate with reaching diagnostic criteria for a major depressive episode. A cut-off of ≤9 is used for remission.
Edinburgh Postnatal Depression Scale (EPDS) [54]. This 10-item measure assesses symptoms of depression among women in the post-natal period. Each item is scored 0–3 and the maximum obtainable score is 30, with scores of ≥13 indicative of caseness. A cut-off of ≤12 is used for remission.
In addition, socio-demographics including age, gender, ethnicity, and employment status were also available.

2.3. Outcomes

Two primary outcomes were included in this analysis, with participants included in the analyses of one outcome only, depending on whether they were or were not in remission at 3 to 4 months:
(1)
Relapse: This outcome is for patients who were in remission at 3 to 4 months and is defined as not in remission at either 6 to 8 or 9 to 12 months on the primary depressive symptom measure used (either BDI-II, PHQ-9, or EPDS) (See Table 1 for primary outcomes).
(2)
Continued chronic course of depression: This outcome is for patients not in remission at 3 to 4 months and is defined as caseness at all subsequent time-points on the primary symptom measure used in the study.

2.4. Data Analysis

2.4.1. Latent Profile Analysis (LPA)

The selection of variables to include in the profiling (termed indicator variables, which can be continuous or categorical) was determined from previous analyses, identifying prognostic factors in the dataset and findings from LPA in similar settings. Baseline depression severity, as well as the duration of depressive illness, have both been identified as strong prognostic factors of treatment outcome, as have initial anxiety symptom severity levels and the duration of comorbid anxiety disorders [50]. In addition to these, the age of participant as well as the employment status of participant were considered as indicators, given their role in discriminating profiles in previous LPA [39]. Seven patient factors, identified at baseline, were included in the LPA. These factors, and how they were transformed for the current analysis, are listed below:
(a)
The sum of the scores on the depression sub-scales of the CIS-R;
(b)
The sum score of the anxiety subscales of the CIS-R:
(c)
Duration of depression made binary: less than one year, or one year, or greater;
(d)
Duration of non-depressive (i.e., anxiety) symptoms measured on the CIS-R made binary: less than one year, or one year or greater;
(e)
History of antidepressant use (yes/no);
(f)
Age at baseline;
(g)
Employment status (employed or unemployed).
A number of metrics were used to identify the optimal LPA solution, in line with established guidance and previous studies employing this analytic approach [39,55]. The following metrics were considered together, with the decision regarding the optimal solution based on a balance between evidence provided by the Vuong-Lo-Mendell-Rubin Likelihood Ratio test (VLMR-LRT) [56], the Bayesian Information Criteria (BIC), sample-size adjusted BIC (SABIC), and entropy values of the model. The VLMR-LRT tests the K model (specified model with K profiles) against the model with one less profile (K-1 model), and a p-value less than 0.05 is taken to indicate the K model is a better fit to the data (a p-value >0.05 indicates little evidence that the K model provides a better fit than the K-1 model). Lower BIC and SABIC values are considered to indicate better model fit, and higher entropy values indicate better classification accuracy of the model. LPA was conducted in Mplus V8 [57].

Association between Profiles and Outcomes

Once the best fitting profile solution was identified, all individuals were allocated to the profiles to which they had the highest probability of membership. Logistic regression models were constructed for both primary outcomes, with profile included as an independent categorical variable to explore differences in outcomes between profiles. All regression analyses were performed in Stata16 [58]; odds ratios (OR) and 95% confidence intervals (95% CI) are presented.
First, differences in the likelihood of relapse or of continued chronic course of depression between the profiles were explored by constructing logistic regression models with allocated treatment (within the RCT, 15 levels in total) included as a covariate (variables used in profiling were not included as covariates in regression models). The difference in each outcome across the profiles as a whole was initially explored through an omnibus test of the variance explained by the model (i.e., an F-test). Because there was evidence of such a difference, the above analyses were conducted comparing the outcomes for participants in each profile relative to both the largest profile as well as another profile which had the average likelihood of remission at 3 to 4 months. The adjusted probability of each outcome was presented for each identified profile, with forest plots constructed to display these differences.
Second, the likelihood of relapse or continued chronic course of depression was compared when participants were randomised to (1) psychological treatments, (2) antidepressant treatment, or (3) ‘Treatment as usual’ or usual care (TAU) arms in their study. Two of the interventions listed in Table 1 (Listening intervention and Physical activity), made up a fourth group labelled ‘Other’, but this was not used in this set of analyses owing to very small numbers of available participants. These groupings were chosen to maximise power whilst providing an informative comparison between intervention types. An omnibus test for the difference in each outcome across all profiles whilst fitting an interaction between profile and treatment type was initially conducted. As there was evidence of differences, analyses were conducted within profiles, to estimate differences in the likelihood of each outcome when different treatments were received for members of each profile, in line with previous analyses looking at treatment differences between identified profiles [40].

2.4.2. Missing Data

Missing data on latent profile indicators was handled in Mplus using Full Information Maximum-Likelihood (FIML), following previous analyses using these approaches [39]. For the analysis exploring the association between outcomes and the identified profiles, individuals with missing outcome data (i.e., lost to follow up) had their outcome data imputed; imputation was also conducted for any participants missing baseline data. Missing baseline and outcome data were imputed using multiple imputation with chained equations (MICE) in Stata 16.0 [58]. This approach uses regression models to impute missing values. A number of imputed datasets (in this study we used 50) were produced to reflect the uncertainty/variability in the imputation process. Where data were not reasonably able to be log transformed to meet normality assumptions, predictive mean matching (PMM) via a k-nearest neighbours approach was used [59]; in this study, we used k = 10 nearest neighbours. Linear regression was used for approximately normally distributed continuous variables, logistic regression models for binary variables, and ordinal and multinomial regression models for ordered and unordered categorical variables, respectively. All imputation models were built using data on baseline and outcome variables following conventions [60]. Only variables with less than 50% missing data were imputed (although in practice no variables had more than 50% missing data; see extended data for degrees of missing by variable [49]), and only participants with missing data on fewer than 50% of the variables were included (this resulted in two participants being removed from the IPCRESS study). Sensitivity analyses including only individuals with complete data were performed to check for differences between imputed and observed model results.

2.5. Ethical Considerations and Trial Registrations

All included studies were granted ethical approvals and all participants gave informed consent (see Supplementary Materials D for details). No additional NHS ethical approval was required for this study: HRA reference 712/86/32/81.

3. Results

3.1. Descriptive Statistics

The included sample is presented in Table 2. The mean age was 43 years (sd = 13.8), with 73% of participants identifying as female and nearly 50% in employment. Just under 50% of the sample reported experiencing depression for over 1 year at baseline, and 73% reported significant anxiety symptoms for over a year at baseline.

3.2. Latent Profile Analysis

The model fit statistics for the profile-solutions are presented in Supplementary Online Table S1. Overall, considering the model fit and classification accuracy statistics, the seven-profile model solution was selected. There was little evidence that an eight-profile solution provided a better fit to the data than the seven-profile model (VLMR-LRT p-value = 0.434), whereas there was evidence that the seven-profile solution was a better fit than the six-profile model (seven-profile VLMR-LRT p-value = 0.019). The AIC and BIC values continued decreasing with the increase in the number of profiles, and the entropy value was highest for the seven-profile model, providing further support that the seven-profile model was the best fit for the data. Each participant was allocated to the profile to which they had the highest probability of membership [40]. The characteristics of individual profiles identified are described below and presented in Table 3.
Profile 1: This was the smallest profile identified (making up 5% of the sample) and was characterised by lower baseline depression and anxiety severity scores compared to all other profiles, whilst having a high likelihood of a duration of anxiety disorders over one year long at baseline.
Profile 2: This profile was the youngest group identified (mean = 29.43 years) and was characterised by patients that were very likely to have had depression for less than one year at baseline and were very likely to be unemployed.
Profile 3: Patients in this profile were similar to Profile 2, in that they had severity scores below the study average and were likely to have suffered from depression for less than a year. However, they were older than Profile 2 members on average (mean age = 40) and were very likely to be employed and have a history of antidepressant use.
Profile 4: This was one of the largest profiles (19% of the sample) and had the second highest initial symptom severity scores (after Profile 7), although the duration of depression for patients in this profile was likely to be less than one year long at baseline.
Profile 5: Members of this profile appeared to be more chronic sufferers of depression as they were highly likely to have been suffering from both depression and anxiety for more than one year, and over three quarters had a history of antidepressant use.
Profile 6: This profile had the oldest average age (62 years), were very likely to be unemployed, and very likely to have both a history of using antidepressants and having suffered from anxiety for over a year.
Profile 7: This group presented with the highest average depression and anxiety scores and a very high likelihood of suffering from depression and anxiety for over a year and having a history of antidepressant use. This was the largest profile identified (22% of the sample).

3.3. Associations between Profiles and Outcomes

Following the identification of the profiles, differences in outcomes between the profiles were explored. Participants were split into those who were in remission at the 3 to 4 months (study) endpoint and were included in the ‘relapse’ outcome analysis, or not in remission and therefore were included in the analysis of the continued chronic course outcome. These categories were mutually exclusive, so participants could not be included in both. Figure 1 presents the proportion of individuals in remission at 3 to 4 months, those who relapsed, and those who followed a continued chronic course of depression, by profile. Profile 1 (low severity) were the most likely to be in remission, not relapse, and not follow a continued chronic course, whereas Profile 7 were least likely to be in remission. Profile 3 (likely first episode, employed) had a near 50% likelihood of remission, a 20% chance of relapse for those previously in remission, and a 50% chance of continued chronic course (if not in remission at 3 to 4 months).
The differences in the likelihood of relapse for individuals who were in remission at 3 to 4 months between profiles is presented in Table 4. There was evidence that the likelihood of relapse was not the same across all profiles (omnibus test p < 0.001). When using Profile 7 (the largest profile) (“high severity”) as the reference group (left side of table), Profile 1 (OR = 0.30 (95% CI = 0.13; 0.68)), Profile 2 (OR = 0.51 (95% CI = 0.27; 0.97)), Profile 3 (OR = 0.38 (95% CI = 0.22; 0.65)), and Profile 5 (OR = 0.44 (95% CI = 0.27; 0.73)) were less likely to relapse. When using Profile 3 (the profile with a 50% likelihood of remission at 3 to 4 months) (“first episode, employed”) as the reference group (right side of Table 3), Profile 4 (OR = 1.82 (95% CI = 1.11; 2.99)), Profile 6 (OR = 1.78 (95% CI = 1.04; 3.04)), and Profile 7 (OR = 2.63 (95% CI = 1.53; 5.43)) were more likely to relapse. The probability of relapse, adjusted by treatment received, and 95% CIs are presented in Figure 2. Sensitivity analyses only conducted with individuals who had complete data at the follow up time periods are presented in Supplementary Online Table S2 and show that all findings were replicated, apart from the decreased odds of relapse for Profile 2 when compared to Profile 7.
Associations between profile and differences in the likelihood of continued chronic course of depression for participants not in remission (i.e., still scoring in the clinical range) at 3 to 4 months between profiles are presented in Table 5. There was evidence that the likelihood of a continued chronic course was not the same across all profiles (omnibus test p < 0.001). Compared to Profile 7, all other profiles were less likely to follow a continued chronic course, and when compared to Profile 3, Profile 4 (OR = 1.84 (95% CI = 1.24; 2.72)), Profile 6 (OR = 1.77 (95% CI = 1.11; 2.82)), and Profile 7 (OR = 3.15 (95% CI = 2.11; 4.71)), were more likely to follow a continued chronic course. The probability of a continued chronic course with 95% CIs are presented in Figure 3. Sensitivity analyses conducted only with individuals that had complete data at the follow up time periods replicated these findings (presented in Supplementary Online Table S3).

3.4. Association between Treatment Type Received and Outcome by Profile

Finally, differences in the likelihood of relapse and of a continued chronic course were compared within each profile when individuals received different types of treatment. We compared the odds of relapse and continued chronic course of depression between individuals who were randomised to: (1) psychological treatments, (2) antidepressant medication(s), or (3) Treatment as usual.
Differences in the odds of relapse between randomised treatments are presented in Table 6. There was evidence that differences between the profiles for each treatment type were not the same (omnibus test p < 0.001). There was little evidence of differences in the likelihood of later relapse between treatment types for individuals from any profile, except for Profile 5, where relapse was more likely for individuals who received treatment as usual (OR = 2.51 (95% CI = 1.16; 5.40)) or antidepressants (OR = 2.92 (95% CI = 1.24; 6.87)) compared to psychological treatments. Table 7 presents the results of logistic regression models comparing the likelihood of chronic course between treatments, and again there was evidence of differences between profile across treatment type received (omnibus test p < 0.001). Findings were similar as for relapse; there was evidence that for those in Profile 5, both TAU (OR = 1.98 (95% CI = 1.16; 3.37)) and antidepressants (OR = 2.27 (95% CI = 1.27; 4.06)) were associated with higher odds of a continued chronic compared to psychological treatment. There was also limited evidence that those in Profile 4 were more likely to relapse with treatment as usual compared to psychological therapy (OR = 2.12 (95% CI = 1.00; 4.51)).
Sensitivity analyses conducted with observed (non-imputed) data are presented in Supplementary Online Tables S4 (relapse) and S5 (continued chronic course). Findings for differences for Profile 5 were replicated. Similar to the primary analyses, the odds of relapse were higher when TAU was received compared to psychological treatment for Profile 4 (OR = 2.65 (95% CI = 1.22; 5.78)) and for Profile 6 when antidepressants were received compared to psychological treatments (OR = 2.91 (95% CI = 1.14; 7.42)). These associations were stronger in the observed data compared to those in the primary analyses using imputed data.

4. Discussion

This study used latent profile analysis to identify stratified subgroups of patients seeking treatment for depression in primary care. Seven distinct profiles were identified, representing a range of subgroups varying in relation to initial symptom severity, comorbidity, duration of illness, age, employment status, and treatment history. Profiles included a low severity group, to a younger first episode group, an older unemployed group with long duration of illness, and a group with very severe symptoms and longer-term depressive illness. Significant differences in both the likelihood of later relapse in individuals who were previously in remission, and of continued chronic course of illness in those who were not in remission, were observed between profiles. The likelihood of relapse ranged from 17.5% (Profile 1) to 45.7% (Profile 7) between profiles, whilst continued chronic depression ranged from 41.5% (Profile 1) to 77.5% (Profile 7). The profiles could be ranked from highest to lowest likelihood in the same order for each outcome. Psychological treatments were associated with a lower likelihood of relapse and continued chronic course for members of Profile 5 (longer duration of depression and anxiety group), compared to both treatment as usual and antidepressant medications.

4.1. Limitations

Whilst this study demonstrates the potential utility of stratification methods in the personalisation of treatment for depression, there are a number of limitations impacting the generalisability of findings. Firstly, this study only included data from seven RCTs, and all were conducted in UK primary care settings, perhaps because the use of the same measure of baseline symptoms and diagnoses (the CIS-R) was an inclusion criterion. The use of the CIS-R was required for data harmonisation purposes, but alternative measures could be used in addition as part of future analyses. The use of RCT data also risks selection biases; for example, there was a lack of ethnic diversity and older mean age of the current sample relative to some other clinical populations with depression in primary care [32]. However, all of the studies were pragmatic trials, increasing the likelihood that the present sample is somewhat representative of other depressed patients presenting in primary care across the world. The fact that all study participants were recruited in primary care also offers an improvement on much of the extant literature, where previously there has often been a lack of information regarding where participants were recruited from [61]. Further, primary care is one of the most common routes into treatment for adults with depression [3,62], so the results here may be generalizable to large proportions of patients with depression. However, the observed profiles may not generalize as well to secondary or inpatient care settings, or to patients seen in different treatment contexts.
The length of follow-up period varied between the studies, and potential changes in depression status between the time points could not be accounted for with the current dataset. This may include further treatments, or complications with antidepressant medication withdrawal, neither of which we could explore in the current dataset. More frequent measurement of depression symptoms following the end of treatment could allow for the identification of the clinical profiles of patients who are likely to relapse earlier, which would further their utility in supporting clinical decision making.
Although data came from participants of RCTs, we analysed data post-treatment and across studies, losing the benefits of randomization in controlling for confounding. Furthermore, a number of potential confounding factors were not available in the dataset, and some of these factors have also been associated with increased risk of relapse. For example, childhood maltreatment, neuroticism, rumination, and interpersonal stress have all been linked to increased risk of relapse [6]. As well as using these factors in future analyses, the utility of profiling approaches to identify at-risk groups may also be enhanced with subgroup specific factors not available here. For example, although gender is not independently prognostic for depression treatment outcomes or course [6,63], and does not typically discriminate between latent profiles of patients with depression [39,40], data on pregnancy or menopause may have helped further refine the latent profiles identified in this study [64,65,66]. In addition, none of the randomised treatments included relapse prevention specific psychological treatments (e.g., mindfulness based cognitive therapy or continuation CBT) [29], a number of common antidepressant medications were not represented in the studies included here, and none included starting combination therapies or augmentation strategies for those not in remission at the primary end-point, so consideration of differential benefits between types of treatment were somewhat limited here. It is also noteworthy that some participants randomized to treatment as usual at baseline may have received a course of antidepressants during the trial. A further limitation to the analysis was due to the focus on continued chronic course of depression and relapse as outcomes. The relapse outcome required participants to have been in remission 3 to 4 months after starting treatment, and therefore individuals who did not remit (over 60% of the sample) were excluded for analyses of this outcome. Whilst unavoidable, this loss of sample will have impacted the available power for analyses. Linked to this, despite the inclusion of seven RCTs, the final sample sizes for comparisons between profiles were small for some estimates, and especially when further grouping by treatment received, limiting the available power to detect differences within profiles. This highlights the need for sufficient sample sizes to explore the further utility of stratification methods to inform clinical practice for adults with depression, and the actual clinical utility of these profiles can only be assessed through a prospective study. Further, there were multiple comparisons made, and it might be argued that adjustments could have been performed to account for the likelihood of Type 1 errors. However, following Rothman [67], and because all analyses followed a pre-registered analysis plan, this was not deemed necessary for this exploratory analysis. Finally, the use of these outcomes might have introduced additional bias in the assessment of differences within profiles across treatment types. If it were the case that one treatment led to a much greater chance of being in remission at 3 to 4 months than another, comparing outcomes after this time-point would have led to considerable selection bias. However, as the results were consistent across the two primary outcomes, this type of selection bias seems to be unlikely in the present study.

4.2. Implications

Findings from this study highlight the potential use of patient stratification methods to identify subgroups of individuals suffering from depression who are more or less likely to follow different courses of illness. Some of these individuals, such as those in Profiles 6 and 7, were at particularly high risk of continued chronic illness, or of later relapse. Conversely, patients in Profiles 1 and 3 were at relatively lower risk of continued chronic depression or relapse. The identification of these groups, potentially at the point of pre-treatment assessment, or a first consultation with a clinician might help inform treatment planning to mitigate these long-term risks. This might involve consideration of more or less regular reviews with a clinician, more aggressive treatments or the augmentation of treatments to prevent chronic depression [30,68], or more conservative management and the use of relapse prevention focussed psychological therapies. Relapse prevention focussed cognitive behavioural therapy, mindfulness based cognitive therapy, or interpersonal psychotherapy have all proven efficacious in reducing the risk of relapse up to one year post-treatment [29]. Information on the risks of chronic depression and relapse might also inform treatment decisions on continuation or discontinuation of medications [69,70]. Additionally, the present study has also provided an analysis which, if replicated in a larger sample, might suggest that psychological therapies could be more beneficial for mitigating the risks of both relapse and continued chronic depression for individuals in Profile 5, relative to treatment as usual or treatment with antidepressant medications. Members of this profile were highly likely to have been suffering from both depression and anxiety for more than one year, and to have a history of antidepressant use, and made up 20% of the analytic sample. Further, in sensitivity analyses, there appeared to be a benefit of psychological therapies compared to treatment as usual and antidepressants for those in Profiles 4 and 6, respectively. If these findings were replicated on external data, it might support calls to increase access to psychological therapies [32] and could potentially inform personalization of care [71]. The risks of both continued chronic depression and relapse were quite high across this study sample. Given the limitations noted above, to better understand how these phenomena can be mitigated and the ways in which treatment can be optimally personalized to do so, it is particularly important that treatment studies continue to collect data to capture these common outcomes, and in routine care services, follow-up appointments and relapse prevention interventions are normalized and integrated into usual practices [32,72].

4.3. Conclusions

The potential utility of patient stratification approaches for supporting clinical decision making have been demonstrated in areas of physical healthcare, and to a lesser extent in predicting response in mental health treatments, but research has not considered the value of stratification for predicting long-term outcomes. This novel study presents a patient stratification approach that identified subgroups of participants who were either more or less likely to either relapse or follow a continued chronic course of depressive illness following treatment for depression. Seven profiles were identified using baseline severity, duration of illness, history of antidepressant use, employment status, and age. Profiles included young people presenting with their first episode, older people with more chronic illness, and individuals with severe symptoms and chronic illness pre-treatment. Substantial differences in the likelihood of relapse and following a continued chronic course between profiles were observed, with some profiles at particular risk of both later relapse and of chronic course of illness. Treatment type was not consistently associated with different outcomes in most profiles, perhaps due to the numbers of participants in stratified profile-by-treatment-groups being low overall, making it difficult to confidently ascertain differences. However, members of Profile 5 appeared to have better outcomes with psychological therapies compared to either antidepressants or treatment as usual, and there was some evidence that members of Profiles 4 and 6 may also potentially have had better outcomes from psychological therapy. This study demonstrates the potential utility of patient identification approaches for ascertaining the likelihood of different long-term outcomes, thereby supporting clinical decision making. The findings might inform subgroups of patients for whom the provision of addition treatments are likely to be needed due to poorer initial response, as well as interventions to reduce the risk of relapse for those who do initially respond.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/jpm11121295/s1, Table S1: Model fit statistics, Table S2: Associations with relapse for each profile compared to two reference groups, the largest overall profile (profile 7) and the profile with approximately equal numbers in remission and not in remission at 3-to-4 months (profile 3) (observed data), Table S3: Associations with continued chronic course of depression for each profile compared two reference groups, the largest overall profile (profile 7) and the profile with approximately equal numbers in remission and not in remission at 3-to-4 months (profile 3) (observed data), Table S4: Association between treatment types and relapse within each profile (observed data), Table S5: Association between treatment type and continued chronic depression within each profile (observed data), Table S6: Ethical approval and Trial Registration details of the included studies from the Dep-GP IPD database.

Author Contributions

R.S. and J.E.J.B. conceived of the project and applied for funding to support this work, along with G.L., S.P., S.D.H., R.J.D., S.G., T.K., E.W. and G.A. All of the above, along with Z.D.C., D.S., N.W., D.K., S.B., E.L. and D.R., wrote the initial protocol document and the plan for the current study. R.S. was responsible for the data analysis, with support from J.E.J.B. and consultation from all other authors. R.S. wrote the original manuscript with support from J.E.J.B. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Welcome Trust through a Clinical Research Fellowship to J.E.J.B. (201292/Z/16/Z), the Royal College of Psychiatrists through a grant to J.E.J.B., R.S. and S.P., MQ Foundation (for ZC: MQDS16/72), the Higher Education Funding Council for England and the National Institute of Health Research (NIHR) University College London Hospitals Biomedical Research Centre (S.P.), University College London (R.S., G.A., G.L.), Vanderbilt University (S.D.H.), University of Southampton (T.K.), University of Exeter (E.W.) and University of York (S.G.), NIHR Biomedical Research Centre at the University Hospitals Bristol and Weston NHS Foundation Trust and the University of Bristol (N.W.). The views expressed are those of the authors and not necessarily those of the NIHR or the Department of Health and Social Care).

Institutional Review Board Statement

All included studies were granted ethical approvals by NHS Research Ethics Committees (see Supplementary Table S6). No additional NHS ethical approval was required for this study: HRA reference 712/86/32/81.

Informed Consent Statement

Written informed consent was obtained from all participants.

Data Availability Statement

Requests for sharing of the IPD used in this study can be made to the corresponding author, any sharing of data will be subject to obtaining appropriate agreements from the chief investigators or data custodians for each individual trial dataset used here.

Acknowledgments

The authors would like to thank Ian R White and Katherine Clarke for their contribution to the general protocol used to develop the Dep-GP dataset which was used in this study.

Conflicts of Interest

The authors declare no conflict of interest..

References

  1. James, S.L.; Abate, D.; Abate, K.H.; Abay, S.M.; Abbafati, C.; Abbasi, N.; Abbastabar, H.; Abd-Allah, F.; Abdela, J.; Abdelalim, A. Global, regional, and national incidence, prevalence, and years lived with disability for 354 diseases and injuries for 195 countries and territories, 1990–2017: A systematic analysis for the Global Burden of Disease Study 2017. Lancet 2018, 392, 1789–1858. [Google Scholar] [CrossRef] [Green Version]
  2. Vos, T.; Allen, C.; Arora, M.; Barber, R.M.; Bhutta, Z.A.; Brown, A.; Carter, A.; Casey, D.C.; Charlson, F.J.; Chen, A.Z.; et al. Global, regional, and national incidence, prevalence, and years lived with disability for 310 diseases and injuries, 1990–2015: A systematic analysis for the Global Burden of Disease Study 2015. Lancet 2016, 388, 1545–1602. [Google Scholar] [CrossRef] [Green Version]
  3. Thornicroft, G.; Chatterji, S.; Evans-Lacko, S.; Gruber, M.; Sampson, N.; Aguilar-Gaxiola, S.; Al-Hamzawi, A.; Alonso, J.; Andrade, L.; Borges, G.; et al. Undertreatment of people with major depressive disorder in 21 countries. Br. J. Psychiatry 2017, 210, 119–124. [Google Scholar] [CrossRef]
  4. Kessler, R.C.; Berglund, P.; Demler, O.; Jin, R.; Koretz, D.; Merikangas, K.R.; Rush, A.J.; Walters, E.E.; Wang, A.; Rovner, B.; et al. The epidemiology of major depressive disorder Results for the National Comorbidity Survey Replication (NCS-R). JAMA 2003, 4, 186–187. [Google Scholar] [CrossRef]
  5. Burcusa, S.L.; Iacono, W.G. Risk for recurrence in depression. Clin. Psychol. Rev. 2007, 27, 959–985. [Google Scholar] [CrossRef] [Green Version]
  6. Buckman, J.E.J.; Underwood, A.; Clarke, K.; Saunders, R.; Hollon, S.D.; Fearon, P.; Pilling, S. Risk factors for relapse and recurrence of depression in adults and how they operate: A four-phase systematic review and meta-synthesis. Clin. Psychol. Rev. 2018, 64, 13–38. [Google Scholar] [CrossRef]
  7. Moffitt, T.E.; Caspi, A.; Taylor, A.; Kokaua, J.; Milne, B.J.; Polanczyk, G.; Poulton, R. How common are common mental disorders? Evidence that lifetime prevalence rates are doubled by prospective versus retrospective ascertainment. Psychol. Med. 2010, 40, 899–909. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Eaton, W.W.; Shao, H.; Nestadt, G.; Lee, B.H.; Bienvenu, O.J.; Zandi, P. Population-Based Study of First Onset and Chronicity in Major Depressive Disorder. Arch. Gen. Psychiatry 2008, 65, 513. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Ali, S.; Rhodes, L.; Moreea, O.; McMillan, D.; Gilbody, S.; Leach, C.; Lucock, M.; Lutz, W.; Delgadillo, J. How durable is the effect of low intensity CBT for depression and anxiety? Remission and relapse in a longitudinal cohort study. Behav. Res. Ther. 2017, 94, 1–8. [Google Scholar] [CrossRef]
  10. Geddes, J.R.; Carney, S.M.; Davies, C.; Furukawa, T.A.; Kupfer, D.J.; Frank, E.; Goodwin, G.M. Relapse prevention with antidepressant drug treatment in depressive disorders: A systematic review. Lancet 2003, 361, 653–661. [Google Scholar] [CrossRef]
  11. McHugh, R.K.; Whitton, S.W.; Peckham, A.D.; Welge, J.A.; Otto, M.W. Patient preference for psychological vs pharmacologic treatment of psychiatric disorders: A meta-analytic review. J. Clin. Psychiatry 2013, 74, 595–602. [Google Scholar] [CrossRef]
  12. Thase, M.E. Preventing Relapse and Recurrence of Depression: A Brief Review of Therapeutic Options. CNS Spectr. 2006, 11, 12–21. [Google Scholar] [CrossRef] [PubMed]
  13. Reimherr, F.W.; Amsterdam, J.D.; Quitkin, F.M.; Rosenbaum, J.F.; Fava, M.; Zajecka, J.; Beasley, C.M.; Michelson, D.; Roback, P.; Sundell, K. Optimal length of continuation therapy in depression: A prospective assessment during long-term fluoxetine treatment. Am. J. Psychiatry 1998, 155, 1247–1253. [Google Scholar] [CrossRef] [PubMed]
  14. Walker, A.J.; Croker, R.; Curtis, H.J.; MacKenna, B.; Goldacre, B. Trends in antidepressant prescribing in England. The Lancet Psychiatry 2021, 8, 278–279. [Google Scholar] [CrossRef]
  15. Simon, G.E.; Perlis, R.H. Personalized Medicine for Depression: Can We Match Patients with Treatments? Am. J. Psychiatry 2010, 167, 1445–1455. [Google Scholar] [CrossRef] [PubMed]
  16. Uher, R.; Perroud, N.; Ng, M.Y.M.; Hauser, J.; Henigsberg, N.; Maier, W.; Mors, O.; Placentino, A.; Rietschel, M.; Souery, D.; et al. Genome-Wide Pharmacogenetics of Antidepressant Response in the GENDEP Project. Am. J. Psychiatry 2010, 167, 555–564. [Google Scholar] [CrossRef]
  17. Fava, G.A.; Ruini, C.; Sonino, N. Treatment of Recurrent Depression: A Sequential Psychotherapeutic and Psychopharmacological Approach. CNS Drugs 2003, 17, 1109–1117. [Google Scholar] [CrossRef] [PubMed]
  18. Lewis, G.; Marston, L.; Duffy, L.; Freemantle, N.; Gilbody, S.; Hunter, R.; Kendrick, T.; Kessler, D.; Mangin, D.; King, M.; et al. Maintenance or Discontinuation of Antidepressants in Primary Care. N. Engl. J. Med. 2021, 385, 1257–1267. [Google Scholar] [CrossRef] [PubMed]
  19. Andrews, G.; Sanderson, K.; Corry, J.; Lapsley, H.M. Using epidemiological data to model efficiency in reducing the burden of depression. J. Ment. Health Policy Econ. 2000, 3, 175–186. [Google Scholar] [CrossRef] [PubMed]
  20. Pies, R. Are antidepressants effective in the acute and long-term treatment of depression? Sic et non. Innov. Clin. Neurosci. 2012, 9, 31–40. [Google Scholar]
  21. Amsterdam, J.D.; Williams, D.; Michelson, D.; Adler, L.A.; Dunner, D.L.; Nierenberg, A.A.; Reimherr, F.W.; Schatzberg, A.F. Tachyphylaxis after repeated antidepressant drug exposure in patients with recurrent major depressive disorder. Neuropsychobiology 2009, 59, 227–233. [Google Scholar] [CrossRef]
  22. Fava, G.A.; Ruini, C.; Rafanelli, C.; Grandi, S. Cognitive Behavior Approach to Loss of Clinical Effect During Long-Term Antidepressant Treatment: A Pilot Study. Am. J. Psychiatry 2002, 159, 2094–2095. [Google Scholar]
  23. Andrews, P.W.; Thomson, J.A.; Amstadter, A.; Neale, M.C. Primum non nocere: An evolutionary analysis of whether antidepressants do more harm than good. Front. Psychol. 2012, 3, 117. [Google Scholar] [CrossRef] [Green Version]
  24. Cuijpers, P.; Hollon, S.D.; van Straten, A.; Bockting, C.; Berking, M.; Andersson, G. Does cognitive behaviour therapy have an enduring effect that is superior to keeping patients on continuation pharmacotherapy? A meta-analysis. BMJ Open 2013, 3, e002542. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Hollon, S.D.; DeRubeis, R.J.; Shelton, R.C.; Amsterdam, J.D.; Salomon, R.M.; O’Reardon, J.P.; Lovett, M.L.; Young, P.R.; Haman, K.L.; Freeman, B.B.; et al. Prevention of Relapse Following Cognitive Therapy vs Medications in Moderate to Severe Depression. Arch. Gen. Psychiatry 2005, 62, 417–422. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Dobson, K.S.; Hollon, S.D.; Dimidjian, S.; Schmaling, K.B.; Kohlenberg, R.J.; Gallop, R.J.; Rizvi, S.L.; Gollan, J.K.; Dunner, D.L.; Jacobson, N.S. Randomized trial of behavioral activation, cognitive therapy, and antidepressant medication in the prevention of relapse and recurrence in major depression. J. Consult. Clin. Psychol. 2008, 76, 468–477. [Google Scholar] [CrossRef]
  27. Klein, D.N.; Santiago, N.J.; Vivian, D.; Blalock, J.A.; Kocsis, J.H.; Markowitz, J.C.; McCullough, J.P., Jr.; Rush, A.J.; Trivedi, M.H.; Arnow, B.A.; et al. Cognitive-behavioral analysis system of psychotherapy as a maintenance treatment for chronic depression. J. Consult. Clin. Psychol. 2004, 72, 681–688. [Google Scholar] [CrossRef]
  28. Bockting, C.L.H. Breaking the rhythm of depression: Cognitive behavior therapy and relapse prevention for depression. Psihol. Teme 2010, 19, 273–287. [Google Scholar]
  29. Clarke, K.; Mayo-Wilson, E.; Kenny, J.; Pilling, S. Can non-pharmacological interventions prevent relapse in adults who have recovered from depression? A systematic review and meta-analysis of randomised controlled trials. Clin. Psychol. Rev. 2015, 39, 58–70. [Google Scholar] [CrossRef]
  30. Denton, W.H.; Wittenborn, A.K.; Golden, R.N. Augmenting antidepressant medication treatment of depressed women with emotionally focused therapy for couples: A randomized pilot study. J. Marital Fam. Ther. 2012, 38, 23–38. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Breedvelt, J.J.F.; Warren, F.C.; Segal, Z.; Kuyken, W.; Bockting, C.L. Continuation of Antidepressants vs Sequential Psychological Interventions to Prevent Relapse in Depression. JAMA Psychiatry 2021, 78, 868–875. [Google Scholar] [CrossRef]
  32. Clark, D.M. Realizing the Mass Public Benefit of Evidence-Based Psychological Therapies: The Iapt Program. Annu. Rev. Clin. Psychol. 2018, 7, 159–183. [Google Scholar] [CrossRef] [Green Version]
  33. Calfee, C.S.; Delucchi, K.; Parsons, P.E.; Thompson, B.T.; Ware, L.B.; Matthay, M.A. Subphenotypes in acute respiratory distress syndrome: Latent class analysis of data from two randomised controlled trials. Lancet. Respir. Med. 2014, 2, 611–620. [Google Scholar] [CrossRef] [Green Version]
  34. Berlin, K.S.; Williams, N.A.; Parra, G.R. An introduction to latent variable mixture modeling (Part 1): Overview and cross-sectional latent class and latent profile analyses. J. Pediatr. Psychol. 2014, 39, 174–187. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Lazarsfield, P.F.; Henry, N.W. Latent Structure Analysis; Houghton Mifflin: Boston, MA, USA, 1968. [Google Scholar]
  36. Scheltens, N.M.E.; Galindo-Garre, F.; Pijnenburg, Y.A.L.; Van Der Vlies, A.E.; Smits, L.L.; Koene, T.; Teunissen, C.E.; Barkhof, F.; Wattjes, M.P.; Scheltens, P.; et al. The identification of cognitive subtypes in Alzheimer’s disease dementia using latent class analysis. J. Neurol. Neurosurg. Psychiatry 2016, 87, 235–243. [Google Scholar] [CrossRef]
  37. Cleynen, I.; John, J.M.M.; Henckaerts, L.; van Moerkercke, W.; Rutgeerts, P.; van Steen, K.; Vermeire, S. Molecular reclassification of Crohn’s disease by cluster analysis of genetic variants. PLoS ONE 2010, 5, e12952. [Google Scholar] [CrossRef] [PubMed]
  38. Unick, G.J.; Snowden, L.; Hastings, J. Heterogeneity in comorbidity between major depressive disorder and generalized anxiety disorder and its clinical consequences. J. Nerv. Ment. Dis. 2009, 197, 215–224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Saunders, R.; Cape, J.; Fearon, P.; Pilling, S. Predicting treatment outcome in psychological treatment services by identifying latent profiles of patients. J. Affect. Disord. 2016, 197, 107–115. [Google Scholar] [CrossRef]
  40. Saunders, R.; Buckman, J.E.J.; Pilling, S. Latent variable mixture modelling and individual treatment prediction. Behav. Res. Ther. 2020, 124, 103505. [Google Scholar] [CrossRef]
  41. Klein, N.S.; Holtman, G.A.; Bockting, C.L.H.; Heymans, M.W.; Burger, H. Development and validation of a clinical prediction tool to estimate the individual risk of depressive relapse or recurrence in individuals with recurrent depression. J. Psychiatr. Res. 2018, 104, 1–7. [Google Scholar] [CrossRef] [PubMed]
  42. Van Loo, H.M.; Aggen, S.H.; Gardner, C.O.; Kendler, K.S. Multiple risk factors predict recurrence of major depressive disorder in women. J. Affect. Disord. 2015, 180, 52–61. [Google Scholar] [CrossRef] [Green Version]
  43. Wang, J.L.; Patten, S.; Sareen, J.; Bolton, J.; Schmitz, N.; Macqueen, G. Development and validation of a prediction algorithm for use by health professionals in prediction of recurrence of major depression. Depress. Anxiety 2014, 31, 451–457. [Google Scholar] [CrossRef]
  44. Moriarty, A.S.; Meader, N.; Snell, K.I.; Riley, R.D.; Paton, L.W.; Chew-Graham, C.A.; Gilbody, S.; Churchill, R.; Phillips, R.S.; Ali, S.; et al. Prognostic models for predicting relapse or recurrence of major depressive disorder in adults. Cochrane Database Syst. Rev. 2021, 2021, 1–57. [Google Scholar] [CrossRef]
  45. Moriarty, A.S.; Paton, L.W.; Snell, K.I.E.; Riley, R.D.; Buckman, J.E.J.; Gilbody, S.; Chew-Graham, C.A.; Ali, S.; Pilling, S.; Meader, N.; et al. The development and validation of a prognostic model to PREDICT Relapse of depression in adult patients in primary care: Protocol for the PREDICTR study. Diagnostic Progn. Res. 2021, 5, 12. [Google Scholar] [CrossRef]
  46. Dwyer, D.B.; Falkai, P.; Koutsouleris, N. Machine Learning Approaches for Clinical Psychology and Psychiatry. Annu. Rev. Clin. Psychol. 2018, 14, 91–118. [Google Scholar] [CrossRef] [PubMed]
  47. Hart, A.; Wyatt, J. Evaluating black-boxes as medical decision aids: Issues arising from a study of neural networks. Informatics Heal. Soc. Care 1990, 15, 229–236. [Google Scholar] [CrossRef] [PubMed]
  48. Kent, P.; Stochkendahl, M.J.; Christensen, H.W.; Kongsted, A. Could the clinical interpretability of subgroups detected using clustering methods be improved by using a novel two-stage approach? Chiropr. Man. Therap. 2015, 23, 1–25. [Google Scholar] [CrossRef] [Green Version]
  49. Buckman, J.E.J.; Saunders, R.; Cohen, Z.D.; Clarke, K.; Ambler, G.; DeRubeis, R.J.; Gilbody, S.; Hollon, S.D.; Kendrick, T.; Watkins, E.; et al. What factors indicate prognosis for adults with depression in primary care? A protocol for meta-analyses of individual patient data using the Dep-GP database. Wellcome Open Res. 2020, 4, 69. [Google Scholar] [CrossRef]
  50. Buckman, J.E.J.; Saunders, R.; Cohen, Z.D.; Barnett, P.; Clarke, K.; Ambler, G.; DeRubeis, R.J.; Gilbody, S.; Hollon, S.D.; Kendrick, T.; et al. The contribution of depressive ‘disorder characteristics’ to determinations of prognosis for adults with depression: An individual patient data meta-analysis. Psychol. Med. 2021, 51, 1068–1081. [Google Scholar] [CrossRef]
  51. Lewis, G.; Pelosi, A.J.; Araya, R.; Dunn, G. Measuring psychiatric disorder in the community: A standardized assessment for use by lay interviewers. Psychol. Med. 1992, 22, 465–486. [Google Scholar] [CrossRef] [PubMed]
  52. Beck, A.T.; Steer, R.A.; Brown, G.K. Manual for the Beck Depression Inventory-II; Psychology Corporation: San Antonio, TX, USA, 1996; pp. 1–82. [Google Scholar]
  53. Kroenke, K.; Spitzer, R.L.; Williams, J.B.W. The PHQ-9: Validity of a Brief Depression Severity Measure. J. Gen. Intern. Med. 2001, 16, 606–613. [Google Scholar] [CrossRef]
  54. Cox, J.L.; Holden, J.M.; Sagovsky, R. Detection of Postnatal Depression: Development of the 10-item Edinburgh Postnatal Depression Scale. Br. J. Psychiatry 1987, 150, 782–786. [Google Scholar] [CrossRef] [Green Version]
  55. Nylund, K.L.; Asparouhov, T.; Muthén, B.O. Deciding on the number of classes in latent class analysis and growth mixture modeling: A Monte Carlo simulation study. Struct. Equ. Model. 2007, 14, 535–569. [Google Scholar] [CrossRef]
  56. Lo, Y.; Mendell, N.R.; Rubin, D.B. Testing the number of components in a normal mixture. Biometrika 2001, 88, 767–778. [Google Scholar] [CrossRef]
  57. Muthén, L.; Muthén, B. Mplus Version 7 User’s Guide; Muthén & Muthén: Los Angeles, CA, USA, 2012. [Google Scholar]
  58. Stata Statistical Software: Release 16; Stata Press: College Station, TX, USA, 2019.
  59. Horton, N.J.; Lipsitz, S.R. Multiple Imputation in Practice: Comparison of Software Packages for Regression Models With Missing Variables. Am. Stat. 2001, 55, 244–254. [Google Scholar] [CrossRef]
  60. Royston, P.; White, I.R. Multiple Imputation by Chained Equations (MICE): Implementation in Stata. J. Stat. Softw. 2011, 45, 1–20. [Google Scholar] [CrossRef] [Green Version]
  61. Dal-Ré, R.; Janiaud, P.; Ioannidis, J.P.A. Real-world evidence: How pragmatic are randomized controlled trials labeled as pragmatic? BMC Med. 2018, 16, 49. [Google Scholar] [CrossRef] [PubMed]
  62. McManus, S.; Bebbington, P.; Jenkins, R.; Brugha, T.S. Mental Health and Wellbeing in England: Adult Psychiatric Morbidity Survey 2014; NHS Digital: Leeds, UK, 2016. [Google Scholar]
  63. Buckman, J.E.J.; Saunders, R.; Stott, J.; Arundell, L.-L.; O’Driscoll, C.; Davies, M.R.; Eley, T.C.; Hollon, S.D.; Kendrick, T.; Ambler, G.; et al. Role of age, gender and marital status in prognosis for adults with depression: An individual patient data meta-analysis. Epidemiol. Psychiatr. Sci. 2021, 30, e42. [Google Scholar] [CrossRef] [PubMed]
  64. Faravelli, C.; Alessandra Scarpato, M.; Castellini, G.; Lo Sauro, C. Gender differences in depression and anxiety: The role of age. Psychiatry Res. 2013, 210, 1301–1303. [Google Scholar] [CrossRef] [PubMed]
  65. Epperson, C.N.; Amin, Z.; Naftolin, F.; Cappiello, A.; Czarkowski, K.A.; Stiklus, S.; Anderson, G.M.; Krystal, J.H. The resistance to depressive relapse in menopausal women undergoing tryptophan depletion: Preliminary findings. J. Psychopharmacol. 2007, 21, 414–420. [Google Scholar] [CrossRef]
  66. Joffe, H.; Cohen, L.S. Estrogen, serotonin, and mood disturbance: Where is the therapeutic bridge? Biol. Psychiatry 1998, 44, 798–811. [Google Scholar] [CrossRef]
  67. Rothman, K.J. No adjustments are needed for multiple comparisons. Epidemiology 1990, 1, 43–46. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Carvalho, A.F.; Cavalcante, J.L.; Castelo, M.S.; Lima, M.C.O. Augmentation strategies for treatment-resistant depression: A literature review. J. Clin. Pharm. Ther. 2007, 32, 415–428. [Google Scholar] [CrossRef] [PubMed]
  69. Berwian, I.M.; Walter, H.; Seifritz, E.; Huys, Q.J.M. Predicting relapse after antidepressant withdrawal—A systematic review. Psychol. Med. 2017, 47, 426–437. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Cromarty, P.; Jonsson, J.; Moorhead, S.; Freeston, M.H. Cognitive behaviour therapy for withdrawal from antidepressant medication: A single case series. Behav. Cogn. Psychother. 2011, 39, 77–97. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Cohen, Z.D.; DeRubeis, R.J. Treatment Selection in Depression. Annu. Rev. Clin. Psychol. 2018, 14, 209–236. [Google Scholar] [CrossRef]
  72. Clark, D.M.; Layard, R.; Smithies, R.; Richards, D.A.; Suckling, R.; Wright, B. Improving access to psychological therapy: Initial evaluation of two UK demonstration sites. Behav. Res. Ther. 2009, 47, 910–920. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Outcomes between profiles.
Figure 1. Outcomes between profiles.
Jpm 11 01295 g001
Figure 2. Adjusted probability of relapse by profile (probability and 95% CIs presented).
Figure 2. Adjusted probability of relapse by profile (probability and 95% CIs presented).
Jpm 11 01295 g002
Figure 3. Adjusted probability of continued chronic course of depression by profile (probability and 95% CIs presented).
Figure 3. Adjusted probability of continued chronic course of depression by profile (probability and 95% CIs presented).
Jpm 11 01295 g003
Table 1. Description of included studies from the Dep-GP IPD dataset.
Table 1. Description of included studies from the Dep-GP IPD dataset.
StudyNSampleInterventionsPrimary Outcome MeasureFollow Up Time Points (Months from Baseline)
CADET581Adults ≥18, ICD-10 depressive episodeCollaborative Care + TAU vs. TAUPHQ-94 and 12
COBALT469Adults 18–75 with treatment resistant depression, scoring ≥14 BDI-IITAU vs. CBT + TAUPHQ-9 (& BDI-II at baseline)3, 6, and 9
IPCRESS299Adults scoring ≥14 BDI-II and GP confirmed diagnosis of depressioniCBT vs. TAUBDI-II4 and 6
MIR480Adults ≥18 taking SSRIs or SNRIs at adequate dose for ≥6 weeks, and scored ≥14 on BDI-IISSRIs or SNRIs + Mirtazapine vs. SSRIs or SNRIs + Pill PlaceboBDI-II (also PHQ-9)3, 6, and 12
REEACT691Adults with PHQ-9 > 10 presenting to GP with depressionMoodgym vs. Beating the Blues vs. TAUPHQ-94 and 12
RESPOND220Women meeting criteria for MDD within 6-months post-partumADM vs. listening interventionEPDS4 and 10
TREAD361Adults 18–69 who met diagnostic criteria for MDD and scored ≥14 on BDI-IITAU vs. Physical Activity + TAUBDI-II4, 7, and 12
Notes: MDD: Major depressive disorder; TAU: Treatment as usual; iCBT: internet-delivered CBT; ADM: Antidepressant medication.
Table 2. Descriptive statistics of the sample.
Table 2. Descriptive statistics of the sample.
Continuous VariablesNM (sd)
Depression severity score303713.8 (3.4)
Anxiety severity score303714.0 (6.6)
Age303742.6 (13.8)
Categorical VariablesCategoryN (%)
Employed (N (%))Yes1511 (49.8%)
No1524 (50.2%)
Missing2 (0.1%)
Duration of depressionLess than 2 weeks42 (1.4%)
2 weeks to 6 months937 (30.9%)
6 to 12 months564 (18.6%)
1 and 2 years457 (15.1%)
More than 2 years1037 (34.2%)
Duration of anxietyLess than 2 weeks6 (0.2%)
2 weeks to 6 months342 (11.3%)
6 to 12 months299 (9.9%)
1 and 2 years320 (10.5%)
More than 2 years1385 (45.6%)
Missing685 (22.6%)
History of antidepressant useYes2069 (68.1%)
No965 (31.8%)
Missing3 (0.1%)
SexMale868 (28.6%)
Female2167 (71.4%)
Missing2 (0.1%)
Ethnicity groupWhite2840 (93.5%)
Black and Minority Ethnic195 (6.4)
Missing2 (0.1%)
Treatment typeTreatment as usual1076 (35.4%)
Psychological interventions1079 (35.5%)
Antidepressant medication589 (19.4%)
Other293 (9.7%)
Table 3. Descriptive statistics for baseline variables in each identified profile.
Table 3. Descriptive statistics for baseline variables in each identified profile.
Profile 1 (5%)Profile 2 (10%)Profile 3 (15%)Profile 4 (19%)Profile 5 (20%)Profile 6 (10%)Profile 7 (22%)
Depression severity *6.1 (2.0)12.4 (3.3)11.8 (2.1)16.0 (1.9)13.0 (2.3)12.2 (2.4)16.8 (2.0)
Anxiety severity *6.2 (3.7)10.3 (4.1)9. 8 (4.2)18.4 (4.7)11.1 (4.4)10.1 (4.5)20.8 (4.9)
Age *41.8 (13.0)29.1 (8.5)40.6 (10.8)37.7 (11.6)41.6 (11.4)63.3 (7.2)46.3 (12.2)
Depression 1+ year §Yes48 (34.5%)5 (1.6%)8 (1.8%)15 (2.6%)607 (99.7%)187 (62.1%)624 (94.8%)
No91 (65.5%)311 (98.4%)433 (98.2%)558 (97.4%)2 (0.3%)114 (37.9%)34 (5.2%)
Anxiety + 1 year §Yes26 (18.7%)65 (20.6%)117 (26.5%)216 (37.7%)467 (76.7%)227 (75.4%)587 (89.2%)
No12 (8.6%)184 (58.2%)154 (34.9%)262 (45.7%)4 (0.7%)26 (8.6%)0
Missing101 (72.7%)67 (21.2%)170 (38.6%)95 (16.6%)138 (22.7%)48 (15.6%)71 (10.8%)
History ADM use §Yes73 (52.5%)144 (45.6%)256 (58.1%)398 (69.5%)455 (74.7%)233 (77.4%)510 (77.5%)
No65 (46.8%)172 (54.4%)184 (41.7%)174 (30.4%)154 (25.3%)68 (22.6%)148 (22.5%)
Missing1 (0.7%)01 (0.2%)1 (0.2%)000
Employed §Yes98 (70.5%)0441 (100%)329 (57.4%)430 (70.6%)0213 (32.4%)
No41 (29.5%)316 (100%)0243 (42.4%)178 (29.2%)301 (100%)445 (67.6%)
Missing0001 (0.2%)1 (0.2%)00
Note: * Means and standard deviations presented; § N and % presented; ADM: antidepressant medication.
Table 4. Associations with relapse for each profile compared to two reference groups: the largest overall profile (Profile 7) and the profile with approximately equal numbers in remission and not in remission at 3 to 4 months (Profile 3).
Table 4. Associations with relapse for each profile compared to two reference groups: the largest overall profile (Profile 7) and the profile with approximately equal numbers in remission and not in remission at 3 to 4 months (Profile 3).
Reference = Profile 7 Reference = Profile 3
OR95% CI OR95% CI
Profile 10.30(0.13; 0.68)Profile 10.79(0.37; 1.70)
Profile 20.51(0.27; 0.97)Profile 21.35(0.73; 2.50)
Profile 30.38(0.22; 0.65)Profile 3RefRef
Profile 40.69(0.41; 1.16)Profile 41.82(1.11; 2.99)
Profile 50.44(0.27; 0.73)Profile 51.16(0.70; 1.93)
Profile 60.68(0.40; 1.15)Profile 61.78(1.04; 3.04)
Profile 7RefRefProfile 72.63(1.53; 5.43)
Note: Odds ratios and confidence intervals are adjusted for the randomised treatment in each RCT.
Table 5. Associations with continued chronic course of depression for each profile compared to two reference groups: the largest overall profile (Profile 7) and the profile with approximately equal numbers in remission and not in remission at 3 to 4 months (Profile 3).
Table 5. Associations with continued chronic course of depression for each profile compared to two reference groups: the largest overall profile (Profile 7) and the profile with approximately equal numbers in remission and not in remission at 3 to 4 months (Profile 3).
Reference = Profile 7 Reference = Profile 3
OR95% CI OR95% CI
Profile 10.23(0.12; 0.46)Profile 10.74(0.67; 1.47)
Profile 20.41(0.25; 0.68)Profile 21.29(0.75; 2.19)
Profile 30.32(0.21; 0.47)Profile 3Ref.Ref.
Profile 40.58(0.41; 0.83)Profile 41.84(1.24; 2.72)
Profile 50.39(0.28; 0.54)Profile 51.24(0.83; 1.83)
Profile 60.56(0.38; 0.84)Profile 61.77(1.11; 2.82)
Profile 7Ref.Ref.Profile 73.15(2.11; 4.71)
Note: Odds ratios and confidence intervals are adjusted for the randomised treatment in each RCT.
Table 6. Association between treatment types and relapse within each profile.
Table 6. Association between treatment types and relapse within each profile.
TAU (vs. Psychological Interventions)Antidepressants (vs. Psychological Interventions)Antidepressants (vs. TAU)
OR95% CIOR95% CIOR95% CI
Profile 1 (n = 87)2.18(0.52; 9.11)4.31(0.57; 32.84)1.98(0.29; 13.68)
Profile 2 (n = 118)2.39(0.76; 7.54)1.02(0.29; 3.58)0.42(0.13; 1.37)
Profile 3 (n = 207)1.15(0.53; 2.46)1.23(0.31; 4.94)1.07(0.26; 4.38)
Profile 4 (n = 178)2.12(1.00; 4.51)0.55(0.14; 2.20)0.26(0.06; 1.07)
Profile 5 (n = 215)2.51(1.16; 5.40)2.92(1.24; 6.87)1.16(0.52; 2.59)
Profile 6 (n = 128)1.92(0.76; 4.83)2.40(0.96; 5.96)1.25(0.50; 3.12)
Profile 7 (n = 120)1.36(0.59; 3.16)1.91(0.70; 5.21)1.40(0.53; 3.70)
Note: Treatment in parentheses is the reference category.
Table 7. Association between treatment type and a continued chronic course of depression within each profile.
Table 7. Association between treatment type and a continued chronic course of depression within each profile.
TAU (vs. Psychological Interventions)Antidepressants (vs. Psychological Interventions)Antidepressants (vs. TAU)
OR95% CIOR95% CIOR95% CI
Profile 1 (n = 50)2.35(0.57; 9.61)2.51(0.45; 14.07)1.07(0.19; 5.95)
Profile 2 (n = 115)1.12(0.41; 3.04)0.61(0.16; 2.30)0.54(0.16; 1.89)
Profile 3 (n = 187)1.70(0.81; 3.54)0.86(0.32; 2.29)0.50(0.19; 1.34)
Profile 4 (n = 303)0.93(0.50; 1.72)1.21(0.49; 2.98)1.30(0.55; 3.04)
Profile 5 (n = 365)1.98(1.16; 3.37)2.27(1.27; 4.06)1.15(0.64; 2.06)
Profile 6 (n = 162)1.18(0.52; 2.68)1.94(0.84; 4.47)1.64(0.73; 3.70)
Profile 7 (n = 509)1.22(0.75; 2.00)1.60(0.84; 3.03)1.31(0.70; 2.43)
Note: Treatment listed in parentheses is the reference category.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Saunders, R.; Cohen, Z.D.; Ambler, G.; DeRubeis, R.J.; Wiles, N.; Kessler, D.; Gilbody, S.; Hollon, S.D.; Kendrick, T.; Watkins, E.; et al. A Patient Stratification Approach to Identifying the Likelihood of Continued Chronic Depression and Relapse Following Treatment for Depression. J. Pers. Med. 2021, 11, 1295. https://doi.org/10.3390/jpm11121295

AMA Style

Saunders R, Cohen ZD, Ambler G, DeRubeis RJ, Wiles N, Kessler D, Gilbody S, Hollon SD, Kendrick T, Watkins E, et al. A Patient Stratification Approach to Identifying the Likelihood of Continued Chronic Depression and Relapse Following Treatment for Depression. Journal of Personalized Medicine. 2021; 11(12):1295. https://doi.org/10.3390/jpm11121295

Chicago/Turabian Style

Saunders, Rob, Zachary D. Cohen, Gareth Ambler, Robert J. DeRubeis, Nicola Wiles, David Kessler, Simon Gilbody, Steve D. Hollon, Tony Kendrick, Ed Watkins, and et al. 2021. "A Patient Stratification Approach to Identifying the Likelihood of Continued Chronic Depression and Relapse Following Treatment for Depression" Journal of Personalized Medicine 11, no. 12: 1295. https://doi.org/10.3390/jpm11121295

APA Style

Saunders, R., Cohen, Z. D., Ambler, G., DeRubeis, R. J., Wiles, N., Kessler, D., Gilbody, S., Hollon, S. D., Kendrick, T., Watkins, E., Richards, D., Brabyn, S., Littlewood, E., Sharp, D., Lewis, G., Pilling, S., & Buckman, J. E. J. (2021). A Patient Stratification Approach to Identifying the Likelihood of Continued Chronic Depression and Relapse Following Treatment for Depression. Journal of Personalized Medicine, 11(12), 1295. https://doi.org/10.3390/jpm11121295

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop