Introduction

Antidepressant drugs are the current standard treatment for major depressive disorder (MDD) and, among these, selective reuptake inhibitors (SSRIs) are the most commonly prescribed. However, their efficacy is variable and incomplete: 60–70% of depressed patients do not experience remission and 30–40% do not show a significant response.1 One of the main reasons for such incomplete efficacy is the poor comprehension of their mechanisms of action. A new hypothesis, the undirected susceptibility to change hypothesis, predicts that SSRI treatment does not drive changes in mood per se but, by increasing brain plasticity, creates a window of opportunity for a change in mood that is driven by the quality of the living conditions.2 In particular, the serotonin increase induced by SSRIs enhances brain plasticity and thus renders the individual more susceptible to the influence of the living conditions. The main consequence is the lack of a univocal outcome of SSRI administration: in a favorable environment treatment leads to a reduction of symptoms; by contrast, in a stressful environment it may lead to a worse prognosis. As a further consequence, SSRIs are expected to amplify the influence of living conditions on mood in a dose-dependent manner. Such hypothesis is supported by preclinical data3, 4, 5 showing that fluoxetine treatment leads to an improvement of the depression-like phenotype when administered in an enriched condition, while it leads to a worsening when administered in a stressful condition. In addition, SSRI treatment consequences on selected end points, such as vulnerability to obesity, have been shown to be dependent on the quality of the environment.6, 7 Finally, a number of clinical studies have shown that the environment moderates the effects of antidepressant treatment.1, 8, 9, 10, 11, 12, 13 However, the approaches used so far in these studies do not allow for assessing the effect of SSRIs on the susceptibility to the influence of the living conditions.

The main aim of the present study was to test whether SSRI treatment amplifies the influence of the living conditions on mood in a dose-dependent manner. We thus exploited the data collected in the framework of the Sequenced Treatment Alternatives to Relieve Depression (STAR*D) study. We considered a subpopulation of 591 patients having similar MDD severity and overlapping treatment history, and analyzed the efficacy of treatment between weeks 4 and 6 according to the dose received—either 20 or 40 mg per day of citalopram. Longer or different treatment periods or other patients’ groups could not be considered without losing data validity because of the limitations imposed by the STAR*D trial design. The socioeconomic characteristics included in the analysis are proxy of the quality of the patient's life environment and widely considered reliable indicators of socioeconomic status.14 Our prediction was that the 40 mg per day dose, compared to the 20 mg per day dose, would amplify the influence of sociodemographic characteristics on patients’ mood, since it should increase plasticity to a higher degree and thus lead to greater susceptibility to the environment. We therefore predicted citalopram to affect susceptibility to the living conditions in a dose-dependent manner. In addition, since we hypothesized that the environment drives the change in MDD induced by SSRIs, we predicted that patients living in conditions associated to a high quality of life should show a more effective response to treatment.

Materials and methods

Study organization

This paper is based on the STAR*D (ClinicalTrials.gov, number NCT00021528) study data. The design details of STAR*D are described accurately elsewhere.15 In brief, STAR*D was a multisite, prospective, randomized, multistep clinical trial conducted in the United States of America aimed to determine which of several treatments would be most effective for outpatients with nonpsychotic MDD.16 The study was conducted at 18 primary care and 23 psychiatric care centers and enrolled 4041 nonpsychotic MDD patients, aged 18–75, with a baseline score 14, on the 17-item Hamilton Depression Rating Scale (HAM-D17).

MDD symptom severity in the STAR*D clinical trial was measured longitudinally using the 16-item Quick Inventory of Depressive Symptomatology (QIDS). The QIDS is a briefer version of the more commonly used 30-item Inventory of Depressive Symptomatology—IDS. The QIDS is available in the clinician and self-rated version and has been designed to assess the severity of depressive symptoms through the evaluation of all the criterion symptom domains designated by the American Psychiatry Association Diagnostic and Statistical Manual of Mental Disorders—4th edition (DSM-IV) to diagnose a major depressive episode. Each item was scored on a scale from 0 to 3 points: 0 indicating no problem and 3 indicating severe problem. Total score ranges from 0, that is, not depressed, to 27, that is, most depressed. For further details on the validity, reliability, generalizability, scoring and interpretation, see http://www.ids-qids.org/.

Treatment

At Level 1, all participants were treated with citalopram, a SSRI, for a minimum of 8 weeks and were encouraged to complete 12 weeks to maximize benefit. All participants started treatment with a dose of 20 mg per day citalopram, with clinical visits at 2, 4, 6, 9 and 12 weeks. To ensure satisfactory dosing for an appropriate period of time, treatment was conducted using measurement-based care. This included flexible dosing recommendations based on symptom and side effect at each treatment visit.17 According to the Guidelines for discontinuing participants from the randomized treatment study, patient drop out could be due to a number of reasons including participant request and decision by clinicians to discontinue the study is in the best interest of the participant.

Dose adjustments were decided according to the QIDS-Clinician-rate (QIDS-C16) score: QIDS-C165, continue current dose; QIDS-C16=6–8, continue or increase current dose according to the clinician assessment; QIDS-C169, increase current dose. In addition, if the reduction in baseline symptom severity was found to be <20% at week 4, the initial dose (20 mg per day) was raised to 40 mg per day, assuming tolerable side effects. At week 4, participants with intolerable side effects could move to the next treatment level.17

Selection criteria

In the present study, data from level 1 of STAR*D were considered. Only patients having similar MDD severity and overlapping treatment history were included in the analysis. In particular, we selected patients treated with a dose of 20 mg per day for the first 4 weeks following enrollment, showing a QIDS-SR16 score equal to 6–10 (mild depression) on week 4 and whose information concerning the QIDS score on week 6 was available (Supplementary Figure S1). Longer or different treatment periods or other patients’ groups could not be considered without losing data validity because of the limitations imposed by the STAR*D trial design. In particular, according to the recommendations of the clinical trial design, patients with high QIDS score that indicates moderate, severe and very severe depression are compulsorily assigned to a dose increase. Including patients with all QIDS scores in the analysis would lead to bring severe depressed and potentially non-responding patients in the high dose group, erroneously leading to the conclusion that the high dose is less effective. The subset of patients considered in the present study show no meaningful clinical difference in depressive symptomatology (measured as QIDS-SR16 score at enrollment, at week 4 and as difference between week 4 and enrollment) according to the dose assignment (see Results section and Table 1). We therefore analyzed treatment efficacy between weeks 4 and 6 according to the dose received, either 20 or 40 mg per day of citalopram.

Table 1 Baseline sociodemographic characteristic of the evaluable sample by citalopram dose

Sociodemographic characteristics

The efficacy of the two dosing regimens was analyzed in relation to a number of sociodemographic characteristics considered as proxy of the socioeconomic status and the quality of the living environment of the patient:14 sex, race, marital status, employment status, insurance status, education, income, experience of traumatic events and drug abuse. In addition, we analyzed the following descriptors of the onset and progression of the psychopathology: age at onset of first major depressive episode (MDE), number of MDE, difference in QIDS-SR16 score between enrollment and week 4. Education was shown as years of schooling completed; we considered two levels of education: ⩾college (16 years). Income was analyzed subdividing the examined population into three monthly income classes: the low-income group was set at $1000, while middle and high income groups at $1000–$2500 and >$2500, respectively. Employment status was analyzed considering only two conditions: employed (self-employed, part-time employed and full-time employed) and unemployed (unemployed and retired). For marital status we considered: never married, married (married or cohabiting) and no more married (separated, divorced and widowed). Finally, we considered only two ethnicities, Caucasian (white) and non-Caucasian (all other ethnicities).

Outcome measures

MDD symptom severity was measured using the QIDS-SR16. Remission was defined as a QIDS-SR16 score 5. In order to determine the effects of citalopram treatment, according to sociodemographic characteristics, we considered three variables: (i) percent of patients showing an improvement, measured as a reduction 1 in QIDS-SR16 score between week 4 and 6; (ii) percent of patients achieving remission, measured as the attainment of a QIDS-SR16 score 5 at week 6; and (iii) variation in the QIDS-SR16 score, measured between weeks 4 and 6.

Data analysis

Of the 4041 participants, a sample of 591 patients was identified according to the selection criteria (Supplementary Figure S1). Summary statistics are presented as percentages for discrete variables, and as means and s.d. for continuous variables. Logistic regression models were used to assess the association of all sociodemographic characteristics with QIDS-SR16 score improvement and remission rate, within each dose group, and to compute odds ratios. The Wald test was used to compute the overall significance of the categorical sociodemographic factors presenting more than two levels. The use of the odds ratios as a measure of risk is debatable and its approximation to the relative risk (RR) is acceptable only when the risk of the event (in the present study, the improvement or remission) in the control group (in the present study, the favorable sociodemographic condition) is very low, that is, under the so-called 'rare disease assumption'.18 For this reason, we additionally computed the percent of improvement or remission in the favorable (that is, associated to a high quality of life) and in the unfavorable (that is, associated to a low quality of life) sociodemographic condition, and the RR of unfavorable versus favorable condition, with the corresponding 95% confidence interval. Finally, the preventive fraction (1−RR) was computed to estimate the effect size and thus the magnitude of influence of the living conditions on the outcome. As in our study the outcome measured is a beneficial one, that is improvement or remission, values of 1−RR>0 indicate that the unfavorable condition decreases the probability of the beneficial outcome compared to the favorable condition, and vice versa. The variation in the QIDS-SR16 score between weeks 4 and 6 was analyzed with analyses of variance including sociodemographic characteristics as between-subject factors. Separate analyses of variance were performed for each sociodemographic characteristic within the two treatment doses. Post hoc comparisons were performed using the Tukey’s test. The variation in percent of improvement, stationary and worsening according to sociodemographic characteristics has been analyzed independently in the two dose groups with X2-test.

Results

Baseline characteristics

Of 591 patients who comprised the evaluable sample, 357 (60.4%) were treated with 20 mg of citalopram per day and 234 (39.6%) with 40 mg of citalopram per day between weeks 4 and 6. Females comprised two-thirds of the sample and minority representation was robust. The two dose groups did not show any meaningful clinical difference before receiving different citalopram doses (Supplementary Table S1). At enrollment, their QIDS-SR16 scores were almost overlapping (20 mg per day dose group: mean 14.14, s.d. 3.78; 40 mg per day dose group: mean 14.7, s.d. 3.79). The two groups did not differ also in QIDS-SR16 score at week 4 (20 mg per day dose group: mean 7.93, s.d. 1.39; 40 mg per day dose group: mean 8.16, s.d. 1.44). Finally, as shown in Table 1, no difference between the 20 and 40 mg per day dose groups was found in response to treatment as both showed an almost overlapping reduction of QIDS-SR score during the first 4 weeks of treatment, when all patients received the 20 mg per day dose. Table 1 summarizes further baseline characteristics of the two groups.

Sociodemographic characteristics and dosage associated with improvement

According to our hypothesis, the two citalopram dosages amplified the influence of sociodemographic characteristics on the percent of patients showing an improvement in a dose-dependent manner. In the 20 mg per day dose group, patients' response was not significantly affected by sociodemographic characteristics. By contrast, in the 40 mg per day dose group, sociodemographic characteristics were associated to significantly different outcomes. In the latter group, a higher rate of improvement was associated with having a working employment status, more than 16 years of education, high income and a private insurance (Table 2). In order to appreciate the influence of the environment in each dose group, its magnitude was measured as the preventive fraction (1−RR). The comparison of the preventive fraction, indicating the percent of patients showing a positive outcome when exposed to a favorable environment, between the two dose groups shows that the influence of sociodemographic characteristics on treatment outcome was overall markedly higher in the 40 mg per day dose group. In particular, the magnitude of the influence of the employment status (unemployed versus employed) and education (⩾college) on improvement rate was, respectively, 5- and 37-fold in the 40 compared to the 20 mg per day dose group. In addition, the magnitude of the influence on improvement rate of income was 6- (high versus low) and 8-fold (high versus middle) and of insurance was 22- (no versus public) and 8-fold (private versus public) in the 40 compared to the 20 mg per day group (Supplementary Table S2).

Table 2 Baseline sociodemographic characteristics associated with percentage of patients improved by citalopram dosea

Rate of remission

As for the rate of improvement, the two citalopram dosages produced a different amplification of the influence of the sociodemographic characteristics on the percent of patients achieving remission. In particular, in the 40 mg per day group, a significantly higher rate of remission was found to be associated with being employed and having longer education (Table 3). The magnitude of the influence of the employment status (unemployed versus employed) and education (⩾college) on remission rate was, respectively, six- and twofold in the 40 compared to the 20 mg per day group (Supplementary Table S3).

Table 3 Baseline sociodemographic characteristics associated with remission by citalopram dosea

Variation in the QIDS-SR16 score

The variation in the QIDS-SR16 score produced results in line with the previous ones, as it significantly differed according to sociodemographic variables only in the 40 mg per day dose group. In particular, in this group, being of Caucasian ethnicity (F (1,232)=5.334, P=0.0218), having a private insurance (F (2,185)=4.427, P=0.0132), a high income (F (2,225)=3.629, P=0.0281) or more years of education (F (1,187)=11.344, P=0.0009) was associated with a significant larger reduction of QIDS-SR16 score (Figure 1). Post hoc analysis revealed a significant reduction of the QIDS-SR16 comparing high versus low income (P<0.05) and having a private versus a public insurance (P<0.05).

Figure 1
figure 1

Variation in the QIDS-SR16 score, measured between weeks 4 and 6, in patients receiving either the 20 or 40 mg per day dose. Negative values indicate a reduction of the QIDS-SR16 score, that is, an improvement of the patient’s mental health condition. Only patients receiving the 40 mg per day dose showed a variation in the QIDS-SR16 significantly moderated by the sociodemographic characteristics. *Indicates difference with the first group; *P<0.05 and ***P<0.001. QIDS-SR16, Quick Inventory of Depressive Symptomatology Self-Report.

Discussion

The results of the present study show that citalopram amplifies the influence of the living conditions on mood in a dose-dependent manner as sociodemographic characteristics modify treatment response in the same direction in the two dose groups, but in the 40 mg per day dose group the effect is much larger and reaches statistical significance. In addition, the magnitude of the influence of the living conditions on mood is much greater—up to 37-fold—in the 40 compared to the 20 mg per day dose group. These results support the undirected susceptibility to change hypothesis,2 which predicts that, as SSRI increases the susceptibility to the environment, treatment outcome is more profoundly affected by the quality of the living conditions in patients receiving high dosages of SSRIs.

The STAR*D clinical trial provides a unique opportunity to investigate the role of citalopram as amplifier of the influence of the living conditions on mood. It has allowed to consider a subpopulation of patients with similar MDD severity and overlapping treatment history in order to analyze the amplification of the influence of the sociodemographic features induced by different dosages of citalopram. Given the clear ethical limitations to perform a clinical trial aimed at directly measuring the effects of the SSRIs in amplifying the beneficial, but especially the detrimental effects of the environment in patients, the STAR*D data set is the best alternative to test the undirected susceptibility to change hypothesis. In addition, the STAR*D clinical trial is the largest ecologically valid 'real world' study of outpatients with nonpsychotic major depressive disorder to date.1, 19

Sociodemographic characteristics modified treatment outcome in the same direction in the two dose groups, but these changes did not reach statistical significance in the 20 mg per day dose group. By contrast, in the 40 mg per day dose group, each one of five sociodemographic characteristics—income, education, ethnicity, insurance and employment—significantly affected treatment outcome. In addition, the magnitude of the effect of the sociodemographic characteristics on mood, measured as the percent of patients showing a positive outcome when exposed to a favorable environment, was dose-dependent. In particular, the influence of the living conditions was much greater in the 40 than in the 20 mg per day dose group: improvement rate went from a minimum of fivefold for employment status to a maximum of 37-fold for education while remission rate went from twofold for education to eightfold for income.

It is worth noting that, in line with our hypothesis and with previous studies showing that a favorable and supportive environment increases antidepressant efficacy,1, 10, 11, 13, 20 the present results indicate not only that citalopram amplifies the influence of the living conditions on mood, but also that the quality of the living environment drives the change in mood. In particular, in the 40 mg per day dose group, where this change reaches statistical significance, improvement and remission were shown at significantly higher rates by patients living in conditions associated with a high quality of life, such as having a working employment status, more than 16 years of education and a high income. The exception concerns insurance where those individuals having public insurance showed the worse outcome, even compared to those having no insurance.21 This is concordant with previous studies reporting, for instance, that having public insurance predicts the highest attrition.17 In addition, higher rates of improvement and remission were associated with being of Caucasian ethnicity. The sociodemographic characteristics here found to affect treatment outcome have been previously shown to both determine rates of major depression morbidity and mortality22 and affect SSRI outcome.1 According to our hypothesis,2 even a worsening of symptomatology could be predicted when citalopram treatment is administered in an adverse environment. However, only a very limited worsening of the QIDS-SR score was expected in patients receiving the treatment while living in an unfavorable condition (Figure 1) because of the implementation of the STAR*D guidelines for discontinuing participants, which recommend that patients drop out of the study when clinicians consider that discontinuation is in the best interest of the participant. Nevertheless, an effect of treatment, both for better and for worse,23 can be appreciated when considering the percent of patients showing improvement or worsening according to the quality of the living conditions in the two dose groups (Supplementary Figure S2). Indeed, the 40, compared to 20 mg per day dose, leads to an increase not only of the percent of improvements in patients living in a favorable environment, but also of the percent of worsening in patients living in an unfavorable environment. For instance, in the 20 mg per day dose group, patients with a college degree showing an improvement were 62% and patients with a high school degree showing a worsening were 21%. In the 40 mg per day group, these percent rise up, respectively, to 71% and 35%. As further example, in the 20 mg per day dose group, patients with a high income showing an improvement were 58% and patients with a low income showing a worsening were 35% while, in the 40 mg per day group, these are, respectively, 68% and 40% (Supplementary Figure S2).

Taking into account the role of SSRIs as amplifier of the influence of living conditions on mood—and, consequently, the quality of the environment as a moderator of SSRI treatment outcome24—may explain the inconsistency of the findings in the literature about SSRI efficacy. Indeed, overlooking the environment as the key stratifying factor in treatment response might have led to the reported contradictory results in antidepressant outcome. The conceptual shift in considering the action of SSRIs from being the causative factor for recovering from major depression to acting as a permissive factor for the influence of the living conditions on the psychopathology allows to reconcile experimental and clinical data that apparently do not fit together.2 For instance, the theoretical framework currently available in the literature leads to the paradox that the same molecular mechanism of action has two opposite outcomes: high extracellular serotonin levels are beneficial when induced by SSRI administration, but confer a high risk to develop psychopathology when associated to the s variant of the serotonin-transporter-linked polymorphic region (5-HTTLPR). This discordant picture can be coherently interpreted in light of our results and the undirected susceptibility to change hypothesis, positing that high serotonin levels lead to increased plasticity and thus to high susceptibility to change, which may promote either an improvement or a worsening, according to the quality of the environment.2 It is worth noting that the effects of citalopram described in the present paper may represent only part of the action of SSRIs on mood as the main target of these drugs, that is, the serotoninergic system, has a high molecular complexity and is involved in a wide range of physiological functions.

The major limitation of the present study is the short time frame (2-week period) over which data of the STAR*D clinical trial have been extracted. This period had to be chosen because it is the only one allowing to consider patients with overlapping clinical history. Longer or different treatment periods or other patient groups would have not allowed to keep data validity. By contrast, the fact that coherent significant results have been found in such a short time frame suggests that the described phenomenon is robust. Further limitations include open treatment design, the use of a single antidepressant agent (citalopram) and the lack of placebo control. Although data analyses did not include correction for multiple comparisons, the overall consistence of the results support their reliability. It is worth noting that depression treatment disparities experienced by different ethnicities may result from stigma, clinician failure to engage with the patient, poor patient activation, treatment adherence and other factors, including biological differences.25

In conclusion, acknowledging the role of SSRIs as amplifier of the influence of the living conditions on mood represents a critical step in developing a personalized medicine approach aimed at better matching patients with treatment and avoiding potential harmful consequences. The control of the patients' living environment could be achieved by training patients to cope with harsh conditions, for instance, through cognitive behavioral therapy,26 or taking them in charge through appropriate specialized services, as it is unlikely that people can rapidly and effectively change their living milieu. The cost of this approach is limited as no new psychoactive molecules need to be developed, while the benefits for the patients could be substantial. Finally, the undirected susceptibility to change hypothesis may contribute building a new theoretical framework capable to integrate the 'chemical imbalance theory' with other hypotheses acknowledging the importance of social–psychological factors in MDD, as both approaches are needed to explain the mechanisms underlying the recovery from the disease.