Introduction
Depression is a common reason for consultation in primary care centers around the world [Reference Finley, Chan, Garrison, Korownyk, Kolber and Campbell1]. Despite there being many clinically recommended treatment options [2], most individuals remain untreated or are not able to access evidence-based treatments [Reference Moitra, Santomauro, Collins, Vos, Whiteford and Saxena3]. Therefore, integrating such treatments into primary care has become a major international healthcare priority [Reference Patel, Saxena, Lund, Thornicroft, Baingana and Bolton4]. Initiatives in England [Reference Clark5], Canada [Reference Naeem, Pikard, Rao, Ayub and Munshi6], Australia [Reference Cromarty, Drummond, Francis, Watson and Battersby7], Norway [Reference Knapstad, Nordgreen and Smith8], and Spain [Reference Cano-Vindel, Muñoz-Navarro, Moriana, Ruiz-Rodríguez, Medrano and González-Blanch9] to integrate psychological therapies into primary care have demonstrated that they are effective, cost-effective, and able to be utilized at scale [Reference Gaskell, Simmonds-Buckley, Kellett, Stockton, Somerville and Rogerson10].
Despite the efficacy of these treatments, a considerable proportion of patients do not improve or achieve a sustained recovery [Reference Saunders, Cohen, Ambler, DeRubeis, Wiles and Kessler11–Reference Cuijpers, Oud, Karyotaki, Noma, Quero and Cipriani14]. Even when patients do achieve recovery, a large number experience relapse or recurrence of depression [Reference Wardenaar, Conradi and de Jonge15, Reference Prieto-Vila, Estupiñá and Cano-Vindel16]. These phenomena are not new; it has been argued for many decades that large proportions of patients with depression will not get better with treatment and large proportions will either experience chronic difficulties or relapse in the years after treatment [Reference Judd17]. However, few studies have investigated the heterogeneity of the course of depression both during therapy and beyond, and this might hold promise for improving long-term outcomes. Modeling the trajectories across time could identify subgroups of patients at risk of poorer outcomes that might then be offered an alternative means of managing their depression, delivering more personalized care [Reference Khoury and Will18].
Prior studies that have modeled the trajectory of symptom change during treatment in primary care for people with depression or anxiety disorders have used routinely data collected from the IAPT services in the United Kingdom [Reference Skelton, Carr, Buckman, Davies, Goldsmith and Hirsch12, Reference Saunders, Buckman, Cape, Fearon, Leibowitz and Pilling13]. These studies found four distinct trajectories/subgroups of change in depressive symptoms: three of them characterized by improvement (small or large) at different moments of the treatment and one characterized by no response or chronicity. However, neither of these studies was able to investigate trajectories during follow-up in the months or years after therapy. A review of observational studies conducted in primary care investigated the proportions of patients in four predetermined subgroups based on their outcomes: recovery, late recovery, relapse, and chronicity [Reference Steinert, Hofmann, Kruse and Leichsenring19]. They found that the proportions of patients in each subgroup varied across the studies, with between 35 and 60% of participants experiencing some sort of recovery, between 10 and 17% having a chronic course, and between 7 and 65% experiencing a relapse. However, this review did not include studies that modeled the trajectories of change. Population studies with large samples and longitudinal follow-up have modeled the trajectories of depressive symptoms, identifying four common trajectories: constantly high or low depressive symptoms and trajectories characterized by increasing or decreasing depressive symptoms over time. Notably, the most prevalent trajectory is constantly low depressive symptoms, observed in more than 70% of the population [Reference Formánek, Csajbók, Wolfová, Kučera, Tom and Aarsland20–Reference Lim, Cheng, Kabir and Thorpe22].
One of the benefits of using methods to identify individuals following distinct trajectories of change is that the association between pre-treatment patient characteristics and the trajectory classes can be investigated. There are a number of important indicators associated with prognosis which can be easily measured in primary care before treatment for depression has commenced, foremost among them is the overall severity of depression [Reference Saunders, Cohen, Ambler, DeRubeis, Wiles and Kessler11–Reference Saunders, Buckman, Cape, Fearon, Leibowitz and Pilling13, Reference Buckman, Saunders, Cohen, Barnett, Clarke and Ambler23]. Other related factors are also associated with prognosis following treatment, including comorbid anxiety symptoms, panic, or somatization symptoms [Reference Skelton, Carr, Buckman, Davies, Goldsmith and Hirsch12, Reference Saunders, Buckman, Cape, Fearon, Leibowitz and Pilling13, Reference Buckman, Saunders, Cohen, Barnett, Clarke and Ambler23]. Other factors appear to be evidenced as prognostic factors, such as employment status [Reference Buckman, Saunders, Stott, Arundell, O’Driscoll and Davies24, Reference McDevitt-Petrovic, Shevlin and Kirby25] and marital status [Reference Wardenaar, Conradi and de Jonge15, Reference Buckman, Saunders, Stott, Cohen, Arundell and Eley26], and there is a group of factors that are not typically evaluated in many studies, but which have been found to be associated with worse prognosis: higher disability [Reference Saunders, Buckman, Cape, Fearon, Leibowitz and Pilling13, Reference Steinert, Hofmann, Kruse and Leichsenring19]; higher suicidal risk or behavior [Reference Wardenaar, Conradi and de Jonge15, Reference McDevitt-Petrovic, Shevlin and Kirby25], higher rumination, cognitive, and attentional biases [Reference Tait, Edmeade and Delgadillo27, Reference Buckman, Saunders, Fearon, Leibowitz and Pilling28]; and higher anhedonia [Reference Khazanov, Xu, Dunn, Cohen, DeRubeis and Hollon29].
Despite the relevance of this area of research, many of the prior studies have been conducted with a narrow group of the population with depression treated in primary care, having largely been conducted in just two countries (i.e., United Kingdom or Netherlands) based on the routine practice datasets of IAPT services or in a randomized controlled trial (RCT) with a small sample size in which a limited set of variables were measured, usually focusing only on clinical variables and disability. This may limit the generalizability to other populations or settings. Therefore, the aims of the current study were as follows: (1) to identify distinct trajectories of long-term changes (12 month follow-up) in depression symptoms in a large sample of patients from an RCT delivered in Spanish primary care settings and (2) to explore the association between baseline patient characteristics and specific trajectories following changes in depression symptoms during the 1-year follow-up.
Method
Participants
The data used for the current study were collected as part of the PsicAP clinical trial [Reference Cano-Vindel, Muñoz-Navarro, Moriana, Ruiz-Rodríguez, Medrano and González-Blanch9], in which patients were randomly allocated to receive either treatment as usual (TAU) and Transdiagnostic – cognitive behavioral treatment (TD-CBT) or TAU alone. In total, 1061 patients with emotional disorders (anxiety, depression, or somatization) were recruited across 22 primary care centers of the Spanish National Health System. TAU was delivered by the patient’s general practitioner (GP) and consisted of psychopharmacotherapy (antidepressants [ADM] or anxiolytics) and/or informal counseling. The TD-CBT was delivered by a clinical psychologist and involved seven group sessions [Reference Muthén and Muthén30]. Patients included in the trial were aged from 18 to 65 years, scored above the cut-off points on one or more of the screening scales for depression, anxiety, or somatoform disorder (PHQ-9 ≥ 10; GAD-7 ≥ 10; PHQ-15 ≥ 5, respectively) and were excluded if they reported: severe symptoms of depression (PHQ-9 ≥ 24); high level of disability (SDS ≥ 26); recent suicidal behavior; were already receiving psychological treatment; had difficulties understanding Spanish; had a diagnosis of substance dependence disorder, or a severe psychological disorder (i.e., personality disorders, eating disorders, bipolar disorder, or a psychotic condition) confirmed by an interview developed by a clinical psychologist. The study was conducted in accordance with the principles of the Declaration of Helsinki and following the Spanish Law on Data Protection (EUDRACT: 2013–001, 955–11). The study protocol was approved by the National Ethics Committee and the Spanish Agency of Medicines and Medical Devices (code: ISRCTN58437086) [Reference Cano-Vindel, Muñoz-Navarro, Wood, Limonero, Medrano and Ruiz-Rodríguez31].
For the purpose of the present study, the eligibility criteria for inclusion were that patients had to demonstrate clinically significant depression symptoms (PHQ-9 ≥ 5) at baseline assessment, and complete the pre-treatment, post-treatment, plus at least one follow-up (3, 6, 12 months) assessment to provide sufficient data for the modeling approach. There were no found statistical differences between the analytical sample for the current study and all participants who started the trial, except on age and marital status (see Supplementary Table 1). A flowchart of the sample is detailed in Supplementary Figure 1.
Data analysis
Trajectory class modeling
Growth mixture modeling (GMM) [Reference Muthén, Brown, Masyn, Jo, Khoo and Yang53] is a longitudinal structural equation modeling approach that aims to identify distinct subgroups of individuals in a sample that demonstrate similar patterns of response over time, employing pre-treatment, post-treatment, and follow-up timepoints (3, 6, and 12 months) PHQ-9 scores. To analyze the subgroups of latent classes, GMMs were performed modeling up to six classes of trajectories, identifying different slopes and intercepts for each number of latent classes, which usually had better fit to the data than the average trajectory [Reference Wardenaar54]. To determine the optimal number of classes, each model (k) was compared to the previous model (k – 1) on the following recommended model fit statistics: the Vuong-Lo-Mendell-Rubin likelihood ratio test (VLMR-LRT), where a p-value of < 0.05 indicates the k model is a better fit for the data than the k−1 model, the Akaike information criterion (AIC) and the Bayesian information criterion (BIC), for which the lowest value between models indicates better fit. Finally, the entropy value of each model was considered, where scores range from 0 to 1 to indicate the accuracy of classification into latent classes, with a value ≥ 0.8 indicating that at least 80% of the individuals were correctly classified into latent classes. Scores between 0.8 and 0.4 indicate medium accuracy and ≤ 0.4 low accuracy [Reference Clark and Muthén55].
GMM analysis was conducted in Mplus version 8.7 [Reference Muthén and Muthén30]. Missing PHQ-9 data were handled using full information maximum likelihood and the expectation maximization algorithm in Mplus [Reference Dempster, Laird and Rubin56].
Association of patient and treatment characteristics with trajectory class
Once the optimum class solution was observed, patients were allocated to the trajectory that they had the highest likelihood of membership before associations between patient characteristics and trajectory was assessed. In accordance with Rothman [Reference Rothman57], the multinomial logistic regression analyses were used without making adjustments for multiple comparisons in order to facilitate the exploration of potential associations between baseline patient characteristics (sociodemographic, clinical, cognitive-emotional, disability and quality of life, and the influence of the treatment received) with each of the identified trajectory classes (see Table 1 for list of variables). SPSS version 27 was used for these analyses [Reference Corp58].
Results
Descriptive statistics
A total sample of 483 patients of the original RCT [Reference Cano-Vindel, Muñoz-Navarro, Moriana, Ruiz-Rodríguez, Medrano and González-Blanch9] met the inclusion criteria for this analysis. Of those, 483 had completed pre-treatment and post-treatment assessments, 414 at 3 months, 361 at 6 months, and 316 at 12 months. Table 2 presents the characteristics of the sample.
Abbreviations: ERQ, Emotional Regulation Questionnaire; GAD-7, Generalized Anxiety Disorder-7; IACTA, Inventory of Cognitive Activity in Anxiety Disorders; MCQ, Metacognition Questionnaire; PHQ-9, Patient Health Questionnaire-9; PHQ-15, Patient Health Questionnaire-15; PSWQ, Penn State Worry Questionnaire; RRS, Rumination Response Scale; SD, standard deviation; SDS, Sheehan Disability Scale; TAU, treatment as usual; TDG-CBT, transdiagnostic group cognitive-behavioral therapy; WHOQOL, World Health Organization Quality of Life.
Trajectories of depressive symptoms
A four-class model was selected as the optimum solution according to the BIC criteria (Table 3). However, according to the VLMR-LRT p-value criterion, the best fit to the data was the two-class model because it is the only one with significant differences from the previous class. Nevertheless, according to Nylund et al. [Reference Nylund, Asparouhov and Muthén59], in case of discrepancies between criteria, the BIC should be chosen given that it is the most consistent indicator. Class allocation resulted in the following trajectory groups (Figure 1).
Abbreviations: Adj-BIC, sample size-adjusted Bayesian information criterion; AIC, Akaike information criterion; BIC, Bayesian information criterion; VLMR-LRT, Vuong-Lo–Mendell–Rubin likelihood ratio test.
Notes: Model fit indices for the favored model are in bold.
Class 1-recovery (n = 310; 64.18%): It is characterized by moderate symptoms at baseline, (PHQ-9mean = 12.6; SD = 4.5), a pronounced decrease at post-treatment (PHQ-9mean = 6.1; SD = 4), and a gradual continuation of symptom reduction at follow-up: 3 months (PHQ-9mean = 5.4; SD = 3.2), 6 months (PHQ-9mean = 5.2; SD = 3.6), and 12 months (PHQ-9mean = 4.9; SD = 3.4).
Class 2-late recovery (n = 49; 10.15%): It is characterized by moderate–severe symptoms at baseline (PHQ-9mean = 16.5; SD = 4.7), post-treatment (PHQ-9mean = 16; SD = 4.9), and 3-month follow-up assessments (PHQ-9mean = 16.2; SD = 4.1), although it shows a gradual symptom reduction was observed at 6 (PHQ-9mean = 10.1; SD = 4.5) and 12-month follow-up assessments (PHQ-9mean = 6.9; SD = 3.5).
Class 3-chronic (n = 58; 12.0%): It is characterized by moderate–severe symptoms at baseline (PHQ-9mean = 18.5; SD = 3.7) and maintaining similar levels throughout all assessments: post-treatment (PHQ-9mean = 18.4; SD = 5.1), 3 months (PHQ-9mean = 20.8; SD = 3.6), 6 months (PHQ-9mean = 20.5; SD = 3.5), and 12 months (PHQ-9mean = 19.3; SD = 4.4).
Class 4-relapse (n = 66; 13.7%): It is characterized by moderate–severe depressive symptoms at baseline (PHQ-9mean = 15.9; SD = 4.6), a decrease at post-treatment (PHQ-9mean = 10.6; SD = 4.2), near to the cut-off point (PHQ-9mean = 10), and a gradual increase of symptoms during the follow-up assessments: 3 months (PHQ-9mean = 10.8; SD = 3.3), 6 months (PHQ-9mean = 12.8; SD = 3.4), and 12 months (PHQ-9mean = 15.8; SD = 3.4).
Associations of baseline variables with trajectory class
The associations between baseline characteristics and depression trajectories were analyzed using multinomial regression models. The Recovery class was used as the reference group as it was the most common trajectory (Table 4).
Abbreviations: ERQ, Emotional Regulation Questionnaire; GAD-7, Generalized Anxiety Disorder-7; IACTA, Inventory of Cognitive Activity in Anxiety Disorders; MCQ, Metacognition Questionnaire; PHQ-9, Patient Health Questionnaire-9; PHQ-15, Patient Health Questionnaire-15; PSWQ, Penn State Worry Questionnaire; RRS, Rumination Response Scale; SD, standard deviation; SDS, Sheehan Disability Scale; TAU, Treatment as usual; TDG-CBT, Transdiagnostic group cognitive-behavioral therapy; WHOQOL, World Health Organization Quality of Life.
Note: * p < 0.05. ** p < 0.01.
The likelihood of being in Class 2-late recovery, relative to Class 1-recovery was higher in patients who had: received TAU alone (OR(95%CI) = 5.18(2.46–10.9)); higher baseline scores for depression (OR(95%CI) = 1.20(1.12–1.29)), generalized anxiety (OR(95%CI) = 1.12(1.01–1.24)), ruminative response (OR(95%CI) = 1.18(1.04–1.33)), and disability (OR(95%CI) = 1.07(1.03–1.12)). Besides, absence of suicidal thoughts (OR(95%CI) = .39(.18–.84)) was associated with a lower likelihood of following the late recovery trajectory respective to recovery.
The likelihood of being in trajectory Class 3-chronicity, compared to Class 1-recovery was higher in those patients that had: received TAU alone (OR(95%CI) = 15.42(6.08–39.14)); higher baseline scores for depression (OR(95%CI) = 1.35(1.24–1.46)), somatization (OR(95%CI) = 1.22(1.10–1.35)); generalized anxiety (OR(95%CI) = 1.20(1.07–1.35)), and emotional suppression (OR(95%CI) = 1.12(1.04–1.2)); and those patients with higher age (OR(95%CI) = 1.05(1.01–1.09)). Besides, higher scores on quality of life (OR(95%CI) = .30(.16–.56)), absence of panic disorder symptoms (OR(95%CI) = .39(.16–.96)) and not taking ADM (OR(95%CI) = .39(.17–.93)) were associated with a lower likelihood of following chronicity than recovery.
The likelihood of being in trajectory Class 4-relapse, compared to Class 1-recovery, was higher in patients with had: basic education level (OR(95%CI) = 2.28(1.09–4.77)); received TAU alone (OR(95%CI) = 2.03(1.11–3.7)); and higher baseline scores for depression (OR(95%CI) = 1.17(1.10–1.24)) and somatization (OR(95%CI) = 1.08(1.01–1.17)).
A binomial regression model was constructed to compare Class 2-late recovery and Class 4-relapse due to similar baseline severity and completely different trajectories across time using Class 4-relapse as the reference group. Higher likelihood of being in Class 2-late recovery was observed for those that had received TAU alone (OR(95%CI) = 2.55(1.1–5.93)) higher ruminative response scores (OR(95%CI) = 1.21(1.04–1.39)); and higher disability scores (OR(95%CI) = 1.06(1.01–1.12)). In addition, higher metacognitive response scores (OR(95%CI) = .87(.77–.99)) were associated with a lower likelihood of following late recovery than relapse (see Table 5).
Abbreviations: ERQ, Emotional Regulation Questionnaire; GAD-7, Generalized Anxiety Disorder-7; IACTA, Inventory of Cognitive Activity in Anxiety Disorders; MCQ, Metacognition Questionnaire; PHQ-9, Patient Health Questionnaire-9; PHQ-15, Patient Health Questionnaire-15; PSWQ, Penn State Worry Questionnaire; RRS, Rumination Response Scale; SD, standard deviation; SDS, Sheehan Disability Scale; TAU, Treatment as usual; TDG-CBT, Transdiagnostic group cognitive-behavioral therapy; WHOQOL, World Health Organization Quality of Life.
Note: * p < 0.05. ** p < 0.01.
Discussion
Main findings
This study identified four different trajectories of change in depressive symptoms during the 1-year follow-up, and baseline variables associated with each trajectory in patients with mild to moderate depressive symptoms seen in primary care. The overall trajectory across the sample was observed to be moderate symptoms at baseline, followed by a pronounced reduction of symptoms at post-treatment, followed by a gradual reduction during follow-up assessments. The four identified trajectories were named “recovery,” “late recovery,” “relapse,” and “chronicity.” The “recovery” trajectory was the most common across patients (64%), and there was a similar proportion following each of the remaining trajectories (10–12%). The four subgroups of patients differed in several baseline characteristics.
The four trajectories found here were similar to those identified in prior studies in primary care [Reference Wardenaar, Conradi and de Jonge15, Reference Steinert, Hofmann, Kruse and Leichsenring19], but differ in prevalence, with fewer patients belonging to the chronicity and relapse trajectories here than in prior studies. This could be due to the different amounts of follow-up time in each study, as more relapses tend to occur in studies with longer follow-up times [Reference Prieto-Vila, Estupiñá and Cano-Vindel16]. Additionally, it may be related to the characteristics of the patients. Our study only included patients with mild to moderate depressive symptoms, whereas in other studies, patients with more severe symptoms were also included, and it is well evidenced that severity is associated with chronicity and relapse [Reference Saunders, Cohen, Ambler, DeRubeis, Wiles and Kessler11, Reference Buckman, Saunders, Cohen, Barnett, Clarke and Ambler23]. Population-based studies have also found heterogeneity in depression symptoms trajectories and an association between baseline severity and chronicity [Reference Formánek, Csajbók, Wolfová, Kučera, Tom and Aarsland20–Reference Lim, Cheng, Kabir and Thorpe22].
It is noticeable that TAU was associated with an increased likelihood of following chronicity (OR(95%CI) = 15.42(6.08–39.14)), late recovery (OR(95%CI) = 5.18(2.46–10.9)) and relapse (OR(95%CI) = 2.03(1.11–3.7)) than recovery. Therefore, these results could suggest that the addition of TDG-CBT is one of the most important factors to achieve a sustained recovery over time. This is consistent with findings from meta-analyses of primary care studies that have indicated the effectiveness of psychological therapy [Reference Gaskell, Simmonds-Buckley, Kellett, Stockton, Somerville and Rogerson10] and better outcomes from psychological therapy compared to TAU [Reference Cuijpers, Oud, Karyotaki, Noma, Quero and Cipriani14]. It is also noteworthy that only 26% of the study sample took ADM, 43% of whom followed the chronicity trajectory. Previous studies showed that patients with prescribed ADM had poorer response to psychological treatment even after adjusting for baseline severity [Reference Skelton, Carr, Buckman, Davies, Goldsmith and Hirsch12, Reference Buckman, Saunders, Cohen, Barnett, Clarke and Ambler23, Reference González-Blanch, Muñoz-Navarro, Medrano, Moriana, Ruiz-Rodríguez and Cano-Vindel60].
Lower levels of comorbid generalized anxiety symptoms, panic disorder symptoms, and somatization were also associated with a greater likelihood of following the recovery trajectory class. This could be because of the common factors of emotional disorders [Reference Hofmann and Barlow61], although both anxiety symptoms and panic disorder have been found to be independently associated with depression treatment outcomes in primary care, even after accounting for depressive symptom severity [Reference Skelton, Carr, Buckman, Davies, Goldsmith and Hirsch12, Reference Saunders, Buckman, Cape, Fearon, Leibowitz and Pilling13, Reference Buckman, Saunders, Cohen, Barnett, Clarke and Ambler23, Reference Khazanov, Xu, Dunn, Cohen, DeRubeis and Hollon29].
Maladaptive emotional regulation strategies are linked to the development and maintenance of emotional disorders [Reference Aldao and Nolen-Hoeksema62], so one of the goals of transdiagnostic psychological therapies is helping patients to develop adaptive emotional regulation strategies [Reference Hofmann and Barlow61]. Previous studies of the PsicAP trial focused on the mediating effect of cognitive processes and emotional regulation strategies on treatment outcomes had shown that patients assigned to TAU+TDG-CBT had significant changes in emotional regulation strategies (worry, rumination, metacognitive beliefs, and emotional suppression) in post-treatment and in depressive symptoms compared to patients who had received TAU alone [Reference Barrio-Martínez, González-Blanch, Priede, Muñoz-Navarro, Medrano and Moriana63]. Similar results were also found for rumination in the UK primary care settings finding that a higher pre-treatment level of rumination was associated with a worse prognosis [Reference Tait, Edmeade and Delgadillo27]. In the present study, we found that higher baseline levels of emotional suppression were associated with the chronicity trajectory class and higher levels of rumination were associated with the late recovery trajectory.
Participants with a basic education were more likely to follow the relapse trajectory compared to the recovery trajectory, consistent with prior studies [Reference Saunders, Cohen, Ambler, DeRubeis, Wiles and Kessler11, Reference Prieto-Vila, Estupiñá and Cano-Vindel16]. However, this association may be influenced by confounding factors such as lower socioeconomic status, housing and other sociodemographic variables associated with depression treatment prognosis in primary care [Reference Buckman, Saunders, Stott, Cohen, Arundell and Eley26], none of which was possible to adjust for in this study. Other studies in primary care have found an association between being unemployed and chronicity [Reference Skelton, Carr, Buckman, Davies, Goldsmith and Hirsch12, Reference McDevitt-Petrovic, Shevlin and Kirby25]. However, in the present study, unemployment was not associated with any trajectory, which could be due to cultural or social reasons.
Differences between late recovery and relapse trajectories are noteworthy given that both groups started with similar baseline severity, but the trajectories are completely different. Patients who followed the late recovery trajectory were characterized by higher baseline scores in disability and rumination, and lower scores in metacognitive beliefs. They were also more likely to have been allocated to receive TAU alone. The comparison between trajectories such as these has not been carried out in previous studies, and the assessment of cognitive-emotional variables, which is an aspect in which both trajectories differ at baseline, are not usually employed either. It was expected that the relapse group would show higher rumination levels at baseline, like other studies [Reference Figueroa, Ruhé, Koeter, Spinhoven, Van der Does and Bockting64]. However, our findings here agree with those that had found that a higher rumination is related with a maintenance and residual symptomatology of depression [Reference Aldao and Nolen-Hoeksema62] which could explain that higher rumination levels are related with the trajectory of late recovery.
Strengths and limitations
This is the first study to investigate trajectories of change in depressive symptoms and associations between pre-treatment patient characteristics and likely trajectory class membership over treatment and up to 1-year follow-up in Spanish primary care services. Likewise, this study is one of the first investigations worldwide based on an RCT that considers longitudinal repeated measures and uses GMM models to identify subgroups with similar patterns in patient trajectories. In addition, this was the first study of psychological therapy to have assessed putative cognitive-emotional mechanisms associated with the onset and maintenance of depression and its resolution during therapy and to consider their associations with different trajectories.
However, there were several limitations to the study. First, a large proportion of patients dropped out following the end of treatment and so did not complete all follow-up assessments, which may have biased the results. It is noteworthy, though, that the rate of drop-out during the post-treatment and follow-up were similar in both treatment groups [Reference Cano-Vindel, Muñoz-Navarro, Moriana, Ruiz-Rodríguez, Medrano and González-Blanch9], and similar to the rate reported in other RCTs conducted in primary care according to the meta-analysis by Bortolotti et al. [Reference Bortolotti, Menchetti, Bellini, Montaguti and Berardi65]. Second, the study lacks information about the clinical history of the patients; therefore, it is likely that the depressive symptoms did not start at the beginning of the trial, and patients with a chronic course may have been more likely to take antidepressants prior to the trial starting. Moreover, the sample for this study consisted of middle-aged women. It was observed that the prevalence of females with depression in the sample was higher compared to the general population in Spain [66] or in primary care settings [Reference King, Nazareth, Levy, Walker, Morris and Weich67]. This difference could be attributed to the fact that this profile could be more inclined to participate in an RCT involving group therapy for emotional disorders. Alternatively, GPs may be more inclined to refer women than men to a study of this nature. Consequently, it is possible that this sample may not be fully representative of patients in primary care.
Conclusions
The present study contributes to a growing and promising area of knowledge on the study of the prognosis of depression and associated baseline characteristics in primary care settings.
These findings highlight the importance of adding TDG-CBT to TAU to improve the chances of recovering from depressive symptoms. Furthermore, it provides a useful tool for the personalization of depression treatment and the identification of subgroups at higher risk based on baseline characteristics. For example, patients exhibiting higher comorbidity, baseline severity of depression, or maladaptive emotional regulation skills were more likely to have a worse prognosis compared to patients with lower scores in these domains. This could have significant implications for the management of depression in primary care and the capacity of these public healthcare services to reduce the burden of depression. While this study offers an initial exploration of personalized treatments, further research, including replication studies and analyses of trajectories per treatment and related characteristics, is necessary before these findings can significantly influence care decisions.
Supplementary material
The supplementary material for this article can be found at http://doi.org/10.1192/j.eurpsy.2024.27.
Author contribution
M.P.-V.: Conceptualization, methodology, formal analysis, investigation, writing – original draft, writing – review and editing. A.C.-V.: Validation, investigation, resources, writing – review and editing, supervision, project administration, funding acquisition. C.G.-B., F.J.E.P., J.E.J.B., R.S.: Conceptualization, methodology, validation, investigation, formal analysis, writing – review and editing, supervision. R.M.-N., J.A.M., S.B.-M., P.R.-R., M.C.-G.: Validation, writing – review and editing.
Financial support
This research was supported by the Spanish Ministry of Science (RETOS grant PID2019-107243RB-C21 and FPI grant PRE2020-092381).
Competing interest
The authors declare none.
Comments
No Comments have been published for this article.