Class 1 indicates brief contact (77.0%); class 2, prolonged initial contact (12.8%); class 3, later reentry (7.1%); and class 4, persistent contact (3.1%). The model included the following covariates: sex, birth year, rural birthplace, age at initial diagnosis, calendar year at initial diagnosis, inpatient treatment at initial diagnosis, record of suicide attempts or self-harm, initial episode severity (mild [reference], moderate, severe without psychotic features, severe with psychotic features, or severity unspecified), and parental history of psychiatric diagnoses in the Danish Psychiatric Central Research Register (depression, bipolar disorder, psychotic illness, substance abuse, and anxiety or somatoform disorders). Data markers indicate predicted probabilities; error bars, 95% CIs.
eFigure 1. Sample Selection Pipeline
eFigure 2. Scree-like Plot of AIC and BIC Values
eFigure 3. Trajectory Patterns in the Sample Including Individuals With No Links to Parents in the Civil Register Who Otherwise Meet Criteria for Inclusion
eFigure 4. Trajectory Patterns in Sample Members With F32 as Their First Recorded MDD Diagnosis
eFigure 5. Trajectory Patterns in Members of the Sample With F33 as Their First Recorded MDD Diagnosis
eFigure 6. Trajectory Patterns in the Sample Including Individuals Who Also Received a Bipolar or Schizophrenia-Spectrum Diagnosis
eResults. Sensitivity Analyses
eTable 1. Characteristics of the Sample Including Individuals With No Links to Parents in the Civil Register Who Otherwise Meet Criteria for Inclusion
eTable 2. Characteristics of Members of the sample with F32 as Their First Recorded MDD Diagnosis
eTable 3. Characteristics of Members of the Sample With F33 as Their First Recorded MDD Diagnosis
eTable 4. Characteristics of the Sample Including Individuals Who Also Received a Bipolar or Schizophrenia-Spectrum Diagnosis
Customize your JAMA Network experience by selecting one or more topics from the list below.
Musliner KL, Munk-Olsen T, Laursen TM, Eaton WW, Zandi PP, Mortensen PB. Heterogeneity in 10-Year Course Trajectories of Moderate to Severe Major Depressive Disorder: A Danish National Register-Based Study. JAMA Psychiatry. 2016;73(4):346–353. doi:10.1001/jamapsychiatry.2015.3365
Evidence suggests that long-term trajectories of major depressive disorder (MDD) are heterogeneous. The Danish Psychiatric Central Research Register (DPCRR) provides a rare opportunity to examine patterns and correlates of long-term trajectories in a large sample of patients with moderate to severe MDD.
To characterize patterns and correlates of 10-year course trajectories of MDD in the DPCRR.
Design, Setting, and Participants
A cohort containing 11 640 individuals born in Denmark in 1955 or later with their first recorded MDD diagnosis in the DPCRR between 1995 and 2002 was established. Patients were followed for 10 years from the date of their initial MDD diagnosis. Data were obtained from Danish civil and psychiatric national registers in June 2013 and were analyzed from April 4, 2014, to December 17, 2015. Correlates of trajectory class membership were sex, characteristics of the first recorded MDD episode (ie, age, severity, inpatient treatment, and record of suicide attempt or self-harm), and psychiatric diagnoses in parents (ie, depression, bipolar disorder, schizophrenia-spectrum disorders, substance abuse, and anxiety or somatoform disorders).
Main Outcomes and Measures
The outcome variable was past-year contact at a psychiatric hospital with a main diagnosis of MDD during each of the 10 years following the initial MDD diagnosis. Trajectories were modeled using latent class growth analysis.
The sample included 11 640 individuals (7493 [64.4%] women) aged 18 to 48 years (mean [SD], 31.4 [7.3]) at their first recorded MDD diagnosis. Four trajectory classes were identified: brief contact (77.0%) (characterized by low probability of contact after 2 years); prolonged initial contact (12.8%) (characterized by high decreasing probability of contact during the first 5 years); later reentry (7.1%) (characterized by moderate probability of contact during the second 5 years); and persistent contact (3.1%) (characterized by high or moderate probability of contact throughout). Female sex (odds ratio [OR] range, 1.82-2.22), inpatient treatment (OR range, 1.40-1.50), and severity at first recorded MDD episode (OR range: moderate, 1.61-1.84; severe, 1.93-2.23; and psychotic, 2.73-3.07) were associated with more severe trajectories. Parental anxiety (OR, 1.34 [95% CI, 1.10-1.63]) and depression (OR, 1.63 [95% CI, 1.28-2.09]) were associated with the prolonged initial contact and later reentry classes, respectively. Parental schizophrenia was associated with the persistent contact class (OR range, 2.55-3.04).
Conclusions and Relevance
Most people treated for moderate to severe MDD in Danish psychiatric hospitals do not receive additional MDD treatment after 2 years; however, a minority receive specialty treatment for up to a decade. Observable heterogeneity in the course may be indicative of underlying etiologic differences.
Evidence suggests that trajectories of long-term illness course among individuals with major depressive disorder (MDD) are heterogeneous. Early work by Merikangas et al,1 Angst,2 and colleagues suggested that single-episode MDD was rare, occurring in fewer than 30% of cases; at the other extreme, an estimated 15% to 20% of patients with MDD exhibited a chronic course. More recently, Eaton et al3 examined the long-term course of MDD in 92 incident community cases from the Epidemiological Catchment Area follow-up. The investigators reported that 50% of the sample had a single episode, 35% had a course characterized by recovery and subsequent recurrence, and 15% experienced chronic MDD. Results from the Collaborative Depression Study,4,5 a longitudinal study of more than 400 patients with clinical MDD, showed that although most patients with MDD recovered from an episode within 1 year, 12% of the patients had not recovered after 5 years and 7% had not recovered after 10 years.
More recently, studies have used group-based trajectory models6-9 to examine heterogeneity in trajectories over time. This approach has many benefits: it is data driven, incorporates measurement error, and can be used to identify latent subgroups within a population as well as correlates of membership in a particular subgroup.7 To date, many studies have used these methods to examine trajectories of depressive symptoms in general population samples.10-34 Few studies, however, have examined heterogeneity in trajectories of illness course among clinical MDD cases,35-39 and those that did were limited by either small sample sizes,35,38 short follow-up periods,37,38 or infrequent outcome measurements.35,39
Characterizing heterogeneity in long-term MDD course trajectories is important for several reasons. First, individuals who follow different trajectories likely have different treatment needs. Second, differences in the patterns of long-term course may indicate differences in the underlying causes of depressive symptoms. Parsing the heterogeneity in the course of depression and identifying factors associated with different course trajectories have the potential to change how we understand and treat depression, and ultimately could improve our ability to intervene effectively and rationally to prevent or mitigate this disorder.
The Danish Psychiatric Central Research Register (DPCRR)40 contains records of all visits (inpatient, outpatient, and emergency) to psychiatric hospitals in Denmark from 1995 onward. Using information from the DPCRR and other Danish national registers, we identified a cohort of individuals first treated for MDD in a psychiatric hospital between 1995 and 2002 and followed up these individuals for 10 years after their first recorded MDD diagnosis. Our goals in this study were to examine heterogeneity in 10-year trajectories of the course of illness among patients in Denmark with moderate to severe MDD and to evaluate the extent to which sex, characteristics of the initial recorded MDD diagnosis, and psychiatric diagnoses in a parent are associated with different trajectory patterns.
Question: What are the patterns and correlates of 10-year trajectories of contact with the Danish specialty psychiatric treatment system for major depressive disorder (MDD)?
Findings: Latent class growth analysis identified 4 trajectory patterns: brief contact (77.0%), prolonged initial contact (12.8%), later reentry (7.1%), and persistent contact (3.1%).
Meaning: Most people treated for moderate to severe MDD in Danish psychiatric hospitals do not receive additional MDD treatment after 2 years; however, a few receive specialty treatment for up to 10 years, suggesting that the burden of MDD is distributed unevenly across cases.
Data for this study were obtained from the Danish national civil and psychiatric registers. The DPCRR40 was established in 1969 and initially recorded only inpatient visits to psychiatric hospitals in Denmark. Starting in 1995, the register began recording all outpatient and emergency visits as well. Among other variables, the DPCRR records the start and end dates of each visit as well as the diagnosis given at discharge by the treating psychiatrist. Diagnoses are based on the International Classification of Diseases, Eighth Revision (ICD-8)41 from 1970 to 1993 and on the International Statistical Classification of Diseases and Related Health Problems, Tenth Revision (ICD-10)42 from 1994 to the present. Visits to primary care physicians for psychiatric conditions are not included in the DPCRR. The Danish Civil Registration System43,44 contains demographic variables, such as sex, birthdate, and birthplace, for all Danish citizens and legal residents since 1968. The Danish Civil Registration System also includes reliable links to parents for individuals born in 1955 or later.44
A detailed description of the sample selection pipeline is presented in eFigure 1 in the Supplement. The source population for the study sample included all individuals born in Denmark in 1955 or later and treated for MDD (ICD-10 codes F32 and F33) in a psychiatric hospital between 1995 and 2002. To reduce prevalence bias, individuals with a record of inpatient MDD treatment in the DPCRR before 1995 (ICD-8 codes 296.09, 296.29, 298.09, and 300.49) were excluded. The sample was further limited to individuals whose first recorded MDD diagnosis occurred after age 18 years because child- and adolescent-onset MDD cases may have a different course trajectory than adult-onset cases. We also excluded individuals with a diagnosis of bipolar disorder, schizophrenia, or schizophrenia-related disorders45; those who died or emigrated within the 10-year period following their initial MDD diagnosis; individuals without links to parents; and persons with a sibling already in the study sample. The final sample contained 11 640 individuals.
This study was approved by the Danish Data Protection Agency. By law in Denmark, informed consent is not required for register-based studies.
Trajectories were modeled using latent class growth analyses6,7,9 in SAS proc traj (SAS, version 9.3; SAS Institute Inc).46,47 The primary response variable was past-year inpatient, outpatient, or emergency contact at a psychiatric hospital with a main diagnosis of MDD. We fit models with 1 to 7 classes, specifying quartic polynomial terms for trajectory slopes. If a polynomial term for a given class was not statistically significant, we ran the model without it; if the resulting model had a lower Bayesian information criterion,48 we used the later model for comparison with the other latent class growth analysis models.
We selected the 4-class model as the final model based on a combination of fit statistics, model precision, and clinical utility (Table 1). Models with additional latent classes had better Bayesian information criterion values (Table 1), but many of the added classes contained small proportions (1%-2%) of the sample and appeared to be minor variations on trajectory patterns already present in the 4-class model. We examined scree-like plots49,50 of the Bayesian information criterion and Akaike information criterion51,52 values (eFigure 2 in the Supplement) and confirmed that improvement in model fit was negligible after 4 classes.
After selecting the 4-class model, we refit it with the following correlates of trajectory class membership: sex, age at first recorded MDD diagnosis, inpatient treatment at first MDD diagnosis, previous record of suicide attempt or self-harm,53 severity of the initial MDD diagnosis (mild, moderate, severe without psychotic features, severe with psychotic features, and severity unspecified), and parental records of psychiatric diagnoses in the DPCRR (unipolar depression, bipolar disorder, schizophrenia and related disorders, substance abuse, and anxiety or somatoform disorders).45 In addition to the primary correlates of interest, we included birth year (1955-1984), calendar year of initial recorded MDD diagnosis (1995-2002), and urban vs rural birthplace in the models to account for secular trends and geographic variability in treatment access. Statistical significance was assessed at the Bonferroni-corrected α level of .017.
We conducted a series of sensitivity analyses (eResults, eTables 1-4, and eFigures 3-6 in the Supplement) to examine the effect of our inclusion or exclusion criteria on the findings. The results of these analyses suggested that none of the inclusion or exclusion criteria had a significant influence on the trajectory patterns or covariate associations.
As a secondary analysis to the main study hypotheses, we examined the differences in the frequencies of comorbid psychiatric diagnoses in the DPCRR during the 10-year follow-up period by trajectory class membership. Comorbid diagnoses included less severe forms of depression, such as dysthymia or recurrent brief depressive episodes (ICD-10 codes F34, F38, and F39), anxiety disorders (codes F40-F43), substance abuse (codes F10-F19), personality disorders (code F60), and other ICD-10 F code diagnoses (other than bipolar or schizophrenia-spectrum disorders). Statistical significance (α = .007) was assessed using χ2 and 2-tailed t tests. Data analysis was conducted from April 4, 2014, to December 17, 2015.
The study sample included 11 640 patients (7493 [64.4%] women) with a mean (SD) age at first recorded MDD diagnosis of 31.4 (7.3) years (Table 2). A total of 2893 individuals (24.9%) received inpatient treatment at their first MDD diagnosis, 1477 (12.7%) had a previous record of suicide attempt or self-harm, 1392 (12.0%) had severe symptoms without psychotic features, and 332 (2.9%) had psychotic symptoms (eg, delusions, hallucinations, or stupor).
The Figure shows the probabilities of contact with the specialty psychiatric treatment system for MDD each year during the 10-year follow-up period for each trajectory class. The largest class (class 1) included 77.0% of the sample (9244 patients with MDD based on maximum posterior probability of class membership) and was characterized by brief contact at a psychiatric hospital for MDD. Individuals in class 1 had a 12% probability of contact for MDD in the DPCRR between the first and second years following their first MDD diagnosis and a 1% probability during the remainder of the 10-year follow-up period. Class 2 (1335 [12.8%]) was characterized by prolonged initial contact. In this class, the probability of contact at a psychiatric hospital for MDD was high (72%) between the first and second year after the initial MDD diagnosis and then declined steadily to 2% to 3% in the second half of the follow-up period. Class 3 (725 [7.1%]) was characterized by later reentry into the specialty psychiatric treatment system. These individuals began with a 25% probability of contact for MDD between the first and second years after initial diagnosis, which decreased to 6% to 8% and then increased to 35% to 45% in the second half of the follow-up period. Class 4 (336 [3.1%]) was characterized by persistent contact. These individuals had a high probability of contact for MDD (80%-90%) throughout the first two-thirds of the follow-up period and a moderately decreasing probability thereafter.
The correlates of trajectory class membership are reported in Table 3. Reported as odds ratios (ORs) (95% CIs), female sex was associated with membership in the prolonged initial illness (1.82 [1.57-2.11]), later reentry (1.41 [1.16-1.71]), and persistent contact (2.22 [1.69-2.92]) classes compared with the brief contact class. Inpatient treatment at the first MDD diagnosis was associated with membership in the prolonged initial contact (1.50 [1.29-1.73]) and later reentry (1.40 [1.15-1.71]) classes compared with the brief contact class. Previous recorded suicide attempt or self-harm was associated with membership in the brief contact (1.39 [1.12-1.72]) and later reentry (1.55 [1.13-2.14]) classes compared with the prolonged initial contact class. There was a dose-response–like association between the severity of the first recorded MDD episode and trajectory class membership: as severity increased, so did the odds of membership in a more severe trajectory class compared with the brief contact class. The exception to this pattern was the later reentry class, for which the effect of psychotic symptoms was not significant.
Reported as ORs (95% CIs), having a parent with a diagnosis of an anxiety or somatoform disorder in the DPCRR was significantly associated with membership in the prolonged initial contact class compared with the brief contact class (1.34 [1.10-1.63]). Parental depression diagnosis was significantly associated with membership in the later reentry class compared with the brief contact (1.63 [1.28-2.09]) and prolonged initial contact (1.49 [1.10-2.03]) classes. Parental diagnosis of schizophrenia or related conditions was significantly associated with membership in the persistent contact class compared with all other classes (brief, 2.55 [1.62-4.02]; prolonged, 2.22 [1.30-3.81]; and later reentry, 3.04 [1.57-5.88]).
We examined the frequency and number of other psychiatric diagnoses received by participants during the 10-year period following their initial MDD diagnosis (Table 4). Overall, 5619 patients (48.3%) received at least 1 additional psychiatric diagnosis during the follow-up period. Predictably, the most common comorbid diagnosis was anxiety (3345 patients [28.7%]). Comorbid psychiatric diagnoses were least common among individuals in the brief contact class (4273 [46.2%] overall; 2304 [24.9%] between 2 and 10 years after initial MDD diagnosis) and most common among individuals in the later reentry class (484 [66.8%]). The elevated burden of comorbid diagnoses in the later reentry class resulted primarily from increased rates of anxiety, substance abuse, and personality disorders rather than milder depression diagnoses.
In this study, we examined heterogeneity in 10-year trajectories of moderate to severe MDD in the DPCRR using latent class growth analysis. The final model contained 4 classes: brief contact (77.0%), prolonged initial contact (12.8%), later reentry (7.1%), and persistent contact (3.1%) with the Danish psychiatric treatment system for MDD. Female sex, inpatient treatment at the first diagnosis, and severity of the initial episode were associated with membership in more severe trajectory classes compared with the brief contact class. Anxiety or somatoform diagnoses and depression diagnoses in parents were associated with membership in the prolonged initial contact and later reentry classes, respectively; schizophrenia and related diagnoses in parents were associated with membership in the persistent contact class.
Our results suggest that most individuals who receive treatment for MDD at a psychiatric hospital in Denmark do not receive additional specialized care for MDD for more than 2 years after their initial diagnosis. Of these persons, 75.1% do not have any contact at a psychiatric hospital with any diagnosis between 2 and 10 years after their initial MDD diagnosis, which suggests that they are functioning better than their counterparts in the other trajectory classes. However, the extent to which individuals in the brief contact group experience residual subclinical depression symptoms and/or receive subsequent care from their primary care physician is unclear.
A notable minority (3.1%) of the moderate to severe cases in our sample followed a trajectory characterized by persistent contact at a psychiatric hospital for MDD. This finding suggests that a large proportion of specialized MDD treatment goes to a small proportion of cases. The finding also suggests that the personal and public health burden of MDD is distributed unevenly across cases, which has important implications for public health and practice. An additional 7.1% of the sample had intermittent contact for MDD, but the high level of psychiatric comorbidities in this group (66.8% had ≥1 additional diagnosis) suggests that many of these patients remain in contact with the Danish psychiatric treatment system for other psychiatric and/or substance use disorders.
Of all the potential correlates, severity of the first diagnosis was most strongly associated with trajectory class membership: the more severe the first diagnosis, the higher the probability of a more severe 10-year trajectory. Inpatient treatment at the first visit was also associated with a more severe trajectory, but the effect was weaker. Female sex was most strongly associated with membership in trajectories characterized by prolonged contact. This is consistent with previous evidence that female sex is a risk factor both for developing depression54 and for experiencing recurrent episodes55,56; however, it may also reflect women’s greater willingness to seek specialized care compared with men. None of these variables demonstrated particular usefulness in distinguishing between the 3 more severe classes. Health care professionals might view these variables as nonspecific indicators of risk for a more severe 10-year course trajectory. The fact that previous recorded suicide attempt or self-harm was associated with membership in the later reentry group (comprising patients with a particularly large burden of comorbid anxiety, substance abuse, and personality disorders) suggests that a suicide attempt or self-harming behavior may be a marker for instability in the course of MDD or in the diagnostic category.
Perhaps the most interesting finding was that different psychiatric diagnoses in parents were associated with different MDD trajectory patterns in offspring. This finding raises the possibility that differences in observable MDD trajectories are influenced by genetic factors, which is consistent with the hypothesis that depression is an etiologically heterogeneous disorder. As such, course trajectory may represent an MDD phenotype that can be used to identify more etiologically homogeneous MDD cases.
Previous evidence for an association between parental psychopathology and the course of MDD is mixed. Lieb et al57 found that parental history of depression was associated with greater severity, chronicity, and recurrence of depression in offspring. However, in their review of factors associated with MDD recurrence, Hardeveld et al58 found that only 1 study59 of the 7 addressing this issue identified an association between family history of MDD and recurrence risk. Angst et al60 failed to find an association between parental history of depression or anxiety and MDD course patterns in the Zurich study. Rhebergen et al37 observed a strong effect (OR, 2.57; 95% CI, 0.75-8.76) of family history of depression on the odds of membership in the chronic, severe, 2-year trajectory class in the Netherlands Study of Depression and Anxiety; however, the association was not statistically significant.
To our knowledge, this is the first study to identify an association between parental history of schizophrenia spectrum disorders and chronic course trajectory in MDD. Recent evidence from the Psychiatric Genomics Consortium61 suggests that there is a genetic association between MDD and schizophrenia. Depressive symptoms are also common among patients with schizophrenia.62,63 The present findings suggest that a parental history of schizophrenia spectrum disorders may confer risk for a subtype of MDD characterized by chronic illness. This area of inquiry is in its infancy; further research on this topic is needed before any firm conclusions can be reached.
Several limitations should be considered when interpreting these results. First, the DPCRR does not include records of patients with MDD who never receive treatment or who receive treatment solely from their primary care physicians. As such, the results of this study may not be generalizable to milder cases of MDD. Although factors such as income, educational level, and geographic location influence treatment-seeking behavior to some degree, Denmark has nearly universal access to health care; psychiatric treatment is free, and the number of psychiatrists per capita is one of the highest in the world.64 Given the relative sparseness of structural barriers to treatment, the factor most likely to differentiate those who receive treatment in a specialty psychiatric hospital (and are therefore included in the DPCRR) from those who receive treatment solely from their primary care physicians or receive no treatment is the severity of the illness.
Second, because the DPCRR began recording outpatient visits in 1995, parents with psychiatric disorders who did not receive inpatient treatment before that time would have been misclassified as having no psychiatric history. This absence of data will have had the largest effect on our measures of parental anxiety and depression since these disorders are far more likely to be treated in outpatient settings than are bipolar disorder or schizophrenia. The associations between parental anxiety or depression and trajectory class membership should therefore be interpreted accordingly.
Finally, given the relative youth of the sample, our results may not be generalizable to late-onset MDD cases. In addition, since inclusion in the study sample was contingent on surviving for at least 10 years after the first MDD diagnosis, these results are not necessarily generalizable to patients with MDD who commit suicide.
Most individuals treated for MDD in psychiatric hospitals in Denmark exit the specialty psychiatric treatment system within a few years after their initial MDD diagnosis. However, a few require treatment for up to a decade. Different psychiatric disorders in parents are associated with different MDD trajectory patterns in their offspring, which suggests that observable heterogeneity in the course of MDD may reflect differences in the genetic underpinnings of the disorder.
Correction: This article was corrected on April 6, 2016, to fix a sentence in the Abstract.
Corresponding Author: Katherine L. Musliner, PhD, National Centre for Register-Based Research, Aarhus University, Fuglesangs Allé 4, Bldg K, 8210 Aarhus V, Denmark (email@example.com).
Submitted for Publication: October 14, 2015; final revision received and accepted December 17, 2015.
Published Online: March 2, 2016. doi:10.1001/jamapsychiatry.2015.3365.
Author Contributions: Dr Musliner had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.
Study concept and design: Musliner, Eaton.
Acquisition, analysis, or interpretation of data: All authors.
Drafting of the manuscript: Musliner.
Critical revision of the manuscript for important intellectual content: All authors.
Statistical analysis: Musliner, Zandi.
Obtained funding: Eaton, Mortensen.
Administrative, technical, or material support: Laursen, Eaton, Mortensen.
Study supervision: Laursen, Eaton, Mortensen.
Conflict of Interest Disclosures: None reported.
Funding/Support: This work was supported by an unrestricted grant from the Lundbeck Foundation: The Lundbeck Foundation Initiative for Integrative Psychiatric Research (iPSYCH). Dr Musliner received funding for doctoral training with the National Institute of Mental Health Psychiatric Epidemiology Institutional Training grant 2T32MH014592-36.
Role of the Funder/Sponsor: The funding sources had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
Additional Contributions: Karen Swartz, MD, Jeannie-Marie Leoutsakos, PhD (Johns Hopkins School of Medicine), Alden Gross, PhD, and Tamar Mendelson, PhD (Johns Hopkins Bloomberg School of Public Health) provided advice and input on the design and execution of the study. There was no financial compensation.