Antidepressants are first-line treatments for major depressive disorder (MDD).1–3 However, many patients fail to respond, with response rates averaging around 50% in clinical trials,Reference Undurraga and Baldessarini4 and it is important to identify these patients as soon as possible to minimise the duration of ineffective treatment and the time until response. Clinical guidelines currently recommend 4–8 weeks of treatment before considering a change in treatment in patients who show no improvement,1–3 although the evidence base for this recommendation is limited. Antidepressant effects can be detected within the first treatment week,Reference Taylor, Freemantle, Geddes and Bhagwagar5 and numerous studies show that early improvement is associated with later response or remission.Reference Stassen, Angst, Hell, Scharfetter and Szegedi6–Reference Wagner, Engel, Engelmann, Herzog, Dreimüller and Müller10 However, these studies disagree on whether lack of early improvement justifies a change in management. For instance, one study found that only 4% of participants with minimal improvement after 2 weeks reached remission by week 4,Reference Szegedi, Jansen, Van Willigenburg, Van Der Meulen, Stassen and Thase7 but in another study, 44% of participants without early improvement still responded after 12 weeks of treatment.Reference Uher, Mors, Rietschel, Rajewska-Rager, Petrovic and Zobel8 On average, most studies indicate that at least 20–30% of participants without early improvement attain response or remission after 4–12 weeks of treatment, which is reduced from the overall probability of around 50%, but is not negligible.Reference Gorwood, Bayle, Vaiva, Courtet, Corruble and Llorca9, Reference Wagner, Engel, Engelmann, Herzog, Dreimüller and Müller10 Conversely, many early improvers do not achieve such outcomes. Hence, better predictive models are desirable.
One possibility to extend these models is to examine individual symptoms, rather than only the total depression score. There are meaningful differences between symptoms (e.g. regarding risk factors and disabilityReference Fried and Nesse11), and severity of specific symptoms is associated with prognosis.Reference Uher, Perlis, Henigsberg, Zobel, Rietschel and Mors12, Reference Chekroud, Zotti, Shehzad, Gueorguieva, Johnson and Trivedi13 Previous studies have found that response can be predicted by early improvement in several symptoms, including depressed mood, somatic symptoms and loss of insight.Reference Tokuoka, Takahashi, Ozeki, Kuga, Yoshikawa and Tsuji14–Reference Farabaugh, Bitran, Witte, Alpert, Chuzi and Clain19 However, these studies did not investigate whether improvement in individual symptoms is more informative than improvement in the total score alone – something we therefore investigated in this study. We also examined whether there are interactions between early-improving symptoms, gender and age. Finally, we examined whether individual symptoms are differentially predictive for response to different antidepressant classes.
Method
Data source and trial selection
We requested individual patient data from Clinical Study Data Request,20 a data-sharing platform providing data from (among others) trials of antidepressants developed by participating sponsors (GlaxoSmithKline and Lilly). We examined second-generation antidepressants (SGAs, defined as selective serotonin reuptake inhibitors (SSRIs), serotonin–norepinephrine reuptake inhibitors (SNRIs) or other antidepressants approved after 1987) because older antidepressants are considered second-line options. However, we also included trials of new chemical entities that have never been approved for treatment of MDD, if an approved SGA was used as an active comparator. We included randomised, placebo- or active-comparator-controlled double-blind trials for MDD in adults, which had a minimum duration of 6 weeks and used the Hamilton Rating Scale for Depression (HRSD) at baseline, week 2 and either week 6 (±1) or week 12 (±1) (or both). We excluded trials that specifically included only participants with additional symptoms (e.g. MDD with pain).
Patient population
We only included participants assigned to placebo or SGAs. No eligible trials included participants assigned to non-SSRI/SNRI SGAs (e.g. mirtazapine), so our final sample consisted of participants assigned to placebo, SSRIs or SNRIs. We took a complete-case approach, only including participants who had valid HRSD scores at baseline, week 2 and week 6 or 12. Week 2 visits took place on day 14 (±7 days), week 6 visits on day 42 (±14 days) and week 12 visits on day 84 (±14). If a participant had multiple visits within the eligible time frame, we selected the visit closest to the intended visit day or, if eligible visits were equally close (e.g. day 35 and day 49), we randomly selected one of the visits.
Training and test data
We randomly split the data into an 80% training set and 20% test set, stratified by treatment group (placebo, SSRI or SNRI). The training set was used for model discovery and cross-validation, and prediction accuracy was assessed in the test set.
Outcomes and predictors
Our primary outcome was response (≥50% reduction in HRSD (17-item version) score) at week 6.Reference Leucht, Fennema, Engel, Kaspers-Janssen, Lepping and Szegedi21 Secondary outcomes were remission (score of ≤7 on the HRSD 17-item version) at week 6, and response and remission at week 12. We chose response as our primary outcome rather than remission because we believed it to be a more realistically achievable outcome by week 6, because it is less dependent on baseline severity and because it is commonly used as the primary outcome in other meta-analyses.Reference Undurraga and Baldessarini4, Reference Cipriani, Furukawa, Salanti, Chaimani, Atkinson and Ogawa22
Improvement in symptoms was calculated from the baseline and week 2 HRSD items and dichotomised into no improvement (no change or worsening) or improvement (improvement in item score of 1 or more) for each individual symptom. Baseline HRSD items were dichotomised into absent (score of 0) or present (score of ≥1). Early improvement in the total HRSD score was dichotomised into no improvement (<20% improvement) or improvement (≥20% improvement), consistent with other studies,Reference Gorwood, Bayle, Vaiva, Courtet, Corruble and Llorca9 whereas the baseline HRSD score was standardised. As demographic variables, we included age (standardised) and gender.
Statistical analysis
Our primary analyses only included antidepressant-treated participants. For variable selection, we used least absolute shrinkage and selection operator (lasso) logistic regression,Reference Tibshirani23 implemented in the glmnet package (version 2.0–5) for R software (version 3.3.0). For each outcome (response and remission at week 6 and 12), we built four models: a baseline model, including only baseline HRSD score, HRSD items, age and gender; a total improvement model, including these baseline variables and early improvement in the total HRSD score; an item improvement model, including all these variables and early improvement in the 17 HRSD items; and an item interactions model, including all of the above and all two-way interactions among early-improving items, age and gender. The optimal regularisation penalty (λ) was determined by ten-fold cross-validation. We favoured sparser models by choosing the largest λ whose deviance was within one s.e. of the minimal deviance.Reference Friedman, Hastie and Tibshirani24
From each lasso model, we selected all variables with non-zero coefficients to build a mixed-effects logistic regression model with a random intercept for trial, using the lme4 package (version 1.1–12). Hence, we built four separate mixed-effects models for each outcome.
Model performance
The prediction accuracy of each mixed-effects model was assessed in the test set by determining the area under the receiver operating characteristic curve (AUC) (R package pROC, version 1.8). The model with the highest AUC was considered the best model. We also determined the accuracy, sensitivity, specificity, positive predictive value (PPV) and negative predictive value (NPV) for each model by assigning participants with a model-predicted probability of response/remission of ≥50% to the response/remission group.
Secondary analyses and post hoc analyses
We performed secondary analyses in the total group of antidepressant- and placebo-treated participants. In these analyses, we included treatment group (placebo versus SSRI versus SNRI) as a predictor to examine whether associations between early-improving items and outcome were dependent on antidepressant class (suggesting a drug-specific mechanism).
In our main and secondary analyses, we dichotomised early improvement, for comparability with other studies. However, we conducted additional post hoc analyses in which baseline item scores and early improvement (change from baseline in the total score and the individual items) were included as continuous variables.
To examine the influence of taking a complete-case approach, we performed single imputation of participants with a valid baseline and week 2 visit, using the mice package for R. We used predictive mean matching to predict HRSD score at week 6 or 12, using as predictors all variables also included in our main analyses. We repeated our main analyses in the imputed data-set.
Ethics statement
Ethical approval was not required for this study as it used only de-identified patient data. The study is registered with Clinicaltrials.gov, under the identifier NCT02934035.
Results
Trials and patients
We requested and received data for 32 trials. However, two trials proved to be ineligible (no week 6 or 12 visit). The remaining 30 trials investigated duloxetine (15 trials), paroxetine (13 trials) or new chemical entities (two trials). Thirteen trials also included other SGAs (escitalopram, fluoxetine, paroxetine or venlafaxine). These 30 trials included 10 365 participants, of whom 8242 had a week 6 visit. The ten trials with a duration of ≥12 weeks included 4487 participants, of whom 3103 had a week 12 visit. Sample characteristics are shown in Table 1. Supplementary Table 1, available at https://doi.org/10.1192/bjp.2018.122, provides further details about the individual trials.
Sample characteristics
HRSD, Hamilton Rating Scale for Depression; SNRI, serotonin–noradrenaline reuptake inhibitor; SSRI, selective serotonin reuptake inhibitor.
Variable selection
Variable selection was performed in the training data-set. Detailed information about the variables selected by the lasso regressions are provided in Supplementary Tables 2–5. In brief, all improvement models selected early improvement in the total score. The item improvement models generally selected most of the early-improving HRSD items. However, items 3 (suicide) and 15 (hypochondria) were never selected, whereas items 1 (depressed mood), 2 (guilt), 4 (early insomnia), 7 (work and activities), 10 (psychological anxiety) and 13 (general somatic symptoms) were always selected. Baseline HRSD items were selected infrequently. All item interaction models selected a number of interactions among early-improving symptoms; additional interactions between early-improving symptoms and age or gender were only selected by the item interaction model for remission at week 12.
Model performance
Model performance was assessed in the test data-set. For response at week 6, the baseline model performed quite poorly (AUC 0.60). The total improvement model performed significantly better (AUC 0.73), and the item improvement and item interactions model performed similarly (AUC 0.77) and significantly better than the total improvement model. For remission at week 6 and response and remission at week 12, the patterns were similar, although model performance was worse for the week 12 outcomes (see Table 2 and Supplementary Figures 1–4 for the receiver operating characteristic curves).
Model performance in the test data-set
AUC, area under the (receiver operating characteristic) curve; N/A, not applicable (undefined); NPV, negative predictive value; PPV, positive predictive value.
The accuracy, sensitivity, specificity, PPV and NPV of each model are also given in Table 2. There were only minor differences between the three early improvement models. At week 6, 51% of antidepressant-treated participants in the test set responded and 33% remitted. The total improvement model predicted non-response for 38% of participants; the associated NPV was 0.74, indicating that 26% of these participants were false negatives who did actually respond by week 6 (Fig. 1). Conversely, of participants who were predicted to respond, 67% actually responded. Based on the most parsimonious model with the highest AUC, the item improvement model, 29% of predicted non-responders actually did respond by week 6 (NPV 0.71), whereas 70% of predicted responders actually responded.
Actual probability of response at week 6 according to participants' predicted outcome (non-response versus response) for each model. Results are based on the test data-set. The dashed line indicates the baseline probability of response. The models predicted non-response for 42% (baseline), 38% (total improvement), 46% (item improvement) and 47% (item interactions) of participants.
For remission at week 6, the total improvement model identified a large majority group (85% of participants) with a slightly reduced probability of remission (28%) and a minority group with an increased probability of remission (64%). Results were similar for the item improvement model (24% and 66% probability of remission, respectively) (Supplementary Figure 5).
By week 12, 69% of participants in the test set responded and 51% remitted. The total improvement model predicted non-response for 13% of participants, but these participants still had a 38% probability of response, whereas predicted non-responders according to the item improvement model still had a 43% probability of response. For remission, 45% of participants were predicted non-remitters according to the total improvement model, and these had a 36% probability of remission compared with a 34% probability of remission according to the item improvement model (see Supplementary Figures 6 and 7).
Post hoc, we also used the predicted probability of responding or remitting to divide participants into quintiles and examined each quintile's actual probability of response or remission (Supplementary Figures 8–11). This more fine-grained approach suggested that the improvement models could identify a risk group with poor outcomes at week 6, but prediction was less accurate and not much better than the baseline model by week 12.
Secondary analyses
Treatment group (placebo versus SSRI versus SNRI) was a significant predictor of response and remission at both week 6 and week 12. However, models that only included a main effect for treatment group performed as well in the test data-set as models with interactions between group and other variables, indicating no evidence for different associations between early-improving symptoms and response or remission depending on antidepressant class (Supplementary Table 6).
Post hoc analyses
Because the total improvement model performed nearly as well as models that included individual items, we performed additional analyses with continuous change from baseline, as dichotomising a continuous variable might affect model performance. For response at week 6, the lasso regressions for the total improvement, item improvement and item interaction models all selected the same variables (baseline score and change from baseline). The AUC of this model in the test data-set was 0.79. For remission at week 6 and response and remission at week 12, the lasso regressions did select individual items for the item improvement and/or the item interactions model. However, these models had similar AUCs as the (more parsimonious) total improvement model. The AUC for the total improvement model was 0.79 for remission at week 6, 0.71 for response at week 12 and 0.75 for remission at week 12 (Supplementary Table 7).
Figure 2 depicts the probability of response or remission as a function of the percentage change from baseline in the total HRSD score at week 2. The probability of response at week 6 was 91% for the few participants (163 out of 6058 participants; 3%) who improved ≥80% by week 2, decreasing to 17% for participants who showed any early worsening (573 participants; 9%). At week 12, however, even participants who showed early worsening still had a 39% probability of responding.
Participants with missing data for week 6 or week 12 outcomes were largely comparable to participants that did have valid data, but were less likely to show early improvement (Supplementary Table 8). In the imputed data-sets, model performance was comparable to our main analyses (Supplementary Table 9), although overall probabilities of response and remission were slightly lower.
Proportion of participants who responded or remitted according to the percentage improvement from baseline. Error bars (dashed lines) indicate the 95% CIs.
Discussion
Principal findings
In this individual patient data meta-analysis, we investigated whether early improvement in individual HRSD symptoms could predict response and remission better than early improvement in the total score alone. Consistent with previous literature, we found that patients without early improvement were less likely to respond or remit.Reference Gorwood, Bayle, Vaiva, Courtet, Corruble and Llorca9, Reference Wagner, Engel, Engelmann, Herzog, Dreimüller and Müller10 In our main analyses, a model with individual symptoms did perform better than a model that only included total improvement. However, the difference was relatively small, and post hoc analyses examining continuous change from baseline did not confirm an added benefit from examining individual symptoms. There was also no evidence that interactions between age, gender and symptoms improved model performance. Our secondary analyses also found no evidence that associations between early-improving symptoms and outcome differed between placebo, SSRIs and SNRIs. Taken together, our results show that early improvement is a non-specific predictor of response and remission, regardless of treatment type.
Although our results confirm the predictive value of early improvement, this value was still relatively limited, especially for longer-term outcomes. Some authors have suggested that non-improvers have virtually no chance of attaining remission and that these patients' treatment should be adapted,Reference Szegedi, Jansen, Van Willigenburg, Van Der Meulen, Stassen and Thase7 but our results indicate that these patients can still achieve response and remission. Further, although it was possible to identify a small group of participants with a very slim chance of attaining remission by week 6 (7%) by looking specifically at risk quintiles, even participants in the highest-risk quintile still had a 21% chance of remitting and a 48% chance of responding by week 12. This decrease in predictive ability for longer-term outcomes suggests that clinicians can anticipate relatively poor outcomes in early non-improvers by week 6, but that a patient's eventual outcome after 12 weeks of treatment cannot be predicted with any certainty after 2 weeks of treatment.
These findings are consistent with results from another large, 12-week trial (the Genome-Based Therapeutic Drugs for Depression study).Reference Uher, Mors, Rietschel, Rajewska-Rager, Petrovic and Zobel8 Another study found that the probability that non-responders would respond within the next 2 weeks was stable throughout the first 12 weeks of treatment, at around 15%.Reference Posternak, Baer, Nierenberg and Fava25 Hence, although lack of early improvement can identify a group of high-risk patients that may need to be monitored more intensively or benefit from additional counselling, a degree of caution in adapting pharmacological treatment may be warranted, especially because the evidence for alternative treatment strategies is scarce. Switching antidepressants is no more effective than continuing the same antidepressant,Reference Bschor, Kern, Henssler and Baethge26 and there is also little evidence in favour of the effectiveness of early dose escalation, although escalation is clearly associated with reduced tolerability.Reference Ruhé, Huyser, Swinkels and Schene27, Reference Dold, Bartova, Rupprecht and Kasper28 This may, however, be different for non-SSRI/SNRI antidepressants like mirtazapine,Reference Ueno, Nakajima, Suzuki, Abe, Sato and Mimura29 or for escalating from very low dosages to standard dosages.Reference Hieronymus, Nilsson and Eriksson30 Other strategies, such as augmentation, may be more successful but also result in decreased tolerability.Reference Zhou, Ravindran A, Qin, Del Giovane, Li and Bauer31 Such strategies might be appropriate for some patients without early improvement, for instance if a fast response is essential because of suicidality, but are likely to be premature for many patients. Given the limited predictive accuracy of models based on symptoms alone, inclusion of a broader set of predictors (e.g. psychiatric history, comorbidity or adverse events) may be necessary to achieve better predictions.
Previous research has indicated that symptoms are not interchangeable, and that the depression sum score could obscure important information.Reference Fried and Nesse11 Several studies have also found that early improvement in specific symptoms is associated with good outcomes,Reference Tokuoka, Takahashi, Ozeki, Kuga, Yoshikawa and Tsuji14–Reference Farabaugh, Bitran, Witte, Alpert, Chuzi and Clain19 in seeming contrast to our work. However, none of these studies included early improvement in the total score, so their findings are not directly comparable to ours. Furthermore, a variety of symptoms were found to be predictive, including general somatic symptoms, gastrointestinal symptoms, insomnia, depressed mood, agitation, loss of interest, feeling slowed down and others, which also suggests that the association between early-improving symptoms and outcome is not particularly specific. Our lasso regressions also tended to select most of the HRSD items rather than a few specific items, although some items were consistently not selected (suicidality and hypochondriasis). These results suggest that, with regard to early improvement, individual symptoms do not add meaningful predictive information to the sum score (especially when taken as continuous). This may be because symptoms are actually more or less interchangeable in this regard, and early improvement in any symptom is associated with good outcomes, or because symptoms are correlated and tend to improve together. However, it could also be related to the reliability of individual items. Single items are more strongly affected by random error than multi-item scales, for which the random error can balance out, degrading the predictive ability of a symptom. Furthermore, as our outcomes were derived from the HRSD sum score, they are inherently dependent on improvement in all individual items, although this would not a priori exclude differences in predictive ability, particularly if the probability or time course of improvement differs.
Our post hoc analyses show that the association between early improvement and outcome is gradual. Although a cut-off, such as ≥20% improvement, may be easier to use in clinical practice, there is no major difference between patients on either side of this cut-off value. The likelihood of response or remission does, however, seem to plateau as the percentage improvement by week 2 drops below around 10%. For instance, the likelihood of response by week 6 is only 17% for patients who deteriorate early in treatment, and the likelihood of remission is only 13%. By week 12, however, around 39% of patients who deteriorate early in treatment have responded and 26% have remitted, which suggests that good outcomes are still possible for these patients (although less likely), if a longer period until remission can be tolerated. These results may therefore offer some guidance to clinicians who are faced with patients showing variable degrees of early improvement and need to decide between continuing, switching or intensifying treatment.
Strengths and limitations
Among the strengths of our study is our large sample size, achieved through combining individual patient data. We used a rigorous approach to building predictive models, including using lasso to prevent over-fitting and using separate test data to examine model performance. We also examined multiple outcomes (response and remission) and both a short and a longer time frame (6 and 12 weeks).
A limitation of our study is that we did not take dosing schedules into account. One study has found that early improvement was more predictive when rapid rather than slow dose escalation was used.Reference Szegedi, Jansen, Van Willigenburg, Van Der Meulen, Stassen and Thase7 However, there is only limited evidence for a dose-response relationship for SGAs,Reference Hieronymus, Nilsson and Eriksson30, Reference Adli, Baethge, Heinz, Langlitz and Bauer32, Reference Jakubovski, Varigonda, Freemantle, Taylor and Bloch33 and dose escalation usually also continues beyond 2 weeks in clinical practice.
We also took a complete-case approach because we were interested in predicting outcomes in patients who are receiving treatment. Our results therefore do not apply to participants who discontinue their medication and drop out of the trial. However, a sensitivity analysis, in which we imputed missing data for participants who did have a baseline and week 2 visit, but did not have week 6 or week 12 visits, yielded similar model performance, although overall response and remission rates were reduced slightly (e.g. from 51.3% to 49.8% for week 6 response).
An additional limitation is that our data were derived from clinical trials with strict inclusion and exclusion criteria. Hence, the study population represents only a subset of treatment-seeking patients, and participants may, on average, have better outcomes than patients seen in clinical practice.Reference Wisniewski, Rush, Nierenberg, Gaynes, Warden and Luther34 Further research is therefore necessary to confirm that our results generalise to the broader patient population, including those with extensive comorbidity or chronic depression.
We also specifically examined the ability of early improvement in individual symptoms to predict response or remission. Other clinically relevant outcomes, such as suicidality or functional impairment, were not examined, and it is conceivable that individual symptoms do add predictive information for such outcomes. Similarly, we specifically examined symptoms derived from the clinician-rated HRSD because the HRSD is most commonly used in depression trials. Other instruments, however, may cover a different set of symptoms,Reference Gullion and Rush35 and there could also be differences between self-reported scales and clinician-rated instruments.Reference Uher, Perlis, Placentino, Dernovšek, Henigsberg and Mors36
Finally, we chose the threshold of ≥50% probability to assign participants to the response or remission category. This is a reasonable cut-off with the advantage of being independent of the data. However, a different cut-off could increase the NPV (at the cost of PPV). In principle, this might identify a group of participants with a lower probability of response or remission. However, because of decreasing specificity, this group would become progressively smaller as NPV increases, which would limit clinical applicability. In post hoc analyses, we examined risk quintiles, which suggested that a small group of participants with poor outcomes at week 6 could be identified, but predictive accuracy was reduced for week 12 outcomes. Similar results were obtained when examining continuous early improvement, which suggests that this is the upper bound of predictive accuracy that can be achieved on the basis of symptoms alone.
Conclusions
Our results show that a model with only early improvement in the total score is about as predictive as models that also contain individual symptoms. Hence, clinicians need not focus on specific symptoms, but can gain as much information about the likelihood of response or remission from improvement in the total score alone, particularly if improvement is interpreted as a continuous measure. However, although absence of early improvement identifies a group of patients that is at increased risk of poor outcomes and that may need to be monitored more closely, it does not rule out later response or remission with certainty. Therefore, adapting antidepressant treatment because of limited improvement in the first 2 weeks would be premature for many patients.
Supplementary material
Supplementary material is available online at https://doi.org/10.1192/bjp.2018.122.
Open access
eLetters
No eLetters have been published for this article.