The diagrams present the 3 multivariate models fitted to the data: correlated factors solution (DSM-5 conceptualization) (A), 1-factor independent pathway model (unidimensional conceptualization) (B), and 2-factor independent pathway model (bidimensional distress and fear conceptualization) (C). The figure is for illustrative purposes only; only the genetic and nonshared environmental associations are shown. The diagram in part A illustrates only the genetic and nonshared environmental correlations between depression and anxiety subscales. Ac and Ac1 indicate additive genetic influences acting via a common factor on all variables; Ac2, additive genetic influences acting via a common factor on 3 fear variables; As, additive genetic influences acting on a specific variable; Ec, nonshared environmental influences acting via a common factor on all variables; and Es, nonshared environmental influences acting on a specific variable.
eAppendix. Recruitment Details
eTable 1. Multivariate Model Fit Statistics in Adolescence and Early Adulthood Excluding Siblings
eTable 2. Multivariate Genetic Analyses at 15 and 17 Years, Inclusive of Physical Injury Variable
eTable 3. Univariate Results
eTable 4. Longitudinal Phenotypic Continuity of Anxiety Subscales Within and Across Anxiety Measures
eTable 5. Multivariate Model Fit Statistics in Adolescence and Early Adulthood, Submodel Comparisons With Drop C and A
Waszczuk MA, Zavos HMS, Gregory AM, Eley TC. The Phenotypic and Genetic Structure of Depression and Anxiety Disorder Symptoms in Childhood, Adolescence, and Young Adulthood. JAMA Psychiatry. 2014;71(8):905-916. doi:10.1001/jamapsychiatry.2014.655
Copyright 2014 American Medical Association. All Rights Reserved. Applicable FARS/DFARS Restrictions Apply to Government Use.
The DSM-5 classifies mood and anxiety disorders as separate conditions. However, some studies in adults find a unidimensional internalizing factor that underpins anxiety and depression, while others support a bidimensional model where symptoms segregate into distress (depression and generalized anxiety) and fear factors (phobia subscales). However, little is known about the phenotypic and genetic structure of internalizing psychopathology in children and adolescents.
To investigate the phenotypic associations between depression and anxiety disorder symptom subscales and to test the genetic structures underlying these symptoms (DSM-5–related, unidimensional and bidimensional) across 3 developmental stages: childhood, adolescence, and early adulthood.
Design, Setting, and Participants
Two population-based prospective longitudinal twin/sibling studies conducted in the United Kingdom. The child sample included 578 twins (mean age, approximately 8 and 10 years at waves 1 and 2, respectively). The adolescent and early adulthood sample included 2619 twins/siblings at 3 waves (mean age, 15, 17, and 20 years at each wave).
Main Outcomes and Measures
Self-report symptoms of depression and anxiety disorders.
Phenotypically, when controlling for other anxiety subscales, depression symptoms were only associated with generalized anxiety disorder symptoms in childhood (r = 0.20-0.21); this association broadened to panic and social phobia symptoms in adolescence (r = 0.17-0.24 and r = 0.14-0.16, respectively) and all anxiety subscales in young adulthood (r = 0.06-0.19). The genetic associations were in line with phenotypic results. In childhood, anxiety subscales were influenced by a single genetic factor that did not contribute to genetic variance in depression symptoms, suggesting largely independent genetic influences on anxiety and depression. In adolescence, genetic influences were significantly shared between depression and all anxiety subscales in agreement with DSM-5 conceptualization. In young adulthood, a genetic internalizing factor influencing depression and all anxiety subscales emerged, alongside a small significant genetic fear factor.
Conclusions and Relevance
These results provide preliminary evidence for different phenotypic and genetic structures of internalizing disorder symptoms in childhood, adolescence, and young adulthood, with depression and anxiety becoming more associated from adolescence. The results inform molecular genetics research and transdiagnostic treatment approaches. The findings affirm the need to continue examining the classification of mood and anxiety disorders in diagnostic systems.
The publication of the DSM-5 has been central to the debate regarding the classification of depression and anxiety disorders.1 Depression and anxiety commonly co-occur2- 4 and are rarely diagnosed in isolation.2,5,6 They share multiple risk factors7 including substantial genetic overlap.8- 12 These observations argue against diagnosis-specific etiology of depression and anxiety. However, anxiety is heterogeneous1 and because of the age changes in internalizing disorders,5,7,13 it remains unclear whether all anxiety types are equally associated with depression across development.14- 18 To improve diagnostic classification, the current study investigated the etiologic structure of internalizing disorder symptoms in childhood, adolescence, and early adulthood.
Most studies investigating the structure of internalizing disorders and symptoms focus on adults. Some studies provide support for a unidimensional internalizing liability factor that underpins anxiety and depression6,19- 23 in line with evidence of shared genetic effects on several different types of anxiety disorders and depression.10,24,25 Another influential conceptualization proposes a bidimensional hierarchical model in which generalized anxiety disorder and depression form a distress factor, while the remaining anxiety disorders form a fear factor.26- 30 These 2 factors may be underpinned by separate genetic influences.31 Importantly, fear and distress are generally highly correlated with each other, thus the 2 conceptualizations are not mutually exclusive.
To our knowledge, few studies to date have used a developmental approach to investigating the structure of internalizing disorder symptoms to test whether the structure is consistent at different developmental stages. Phenotypic studies in children and adolescents provide mixed conclusions. Some support a unidimensional internalizing factor,32- 34 others identify the distress and fear dimensions,35,36 and 1 study found that depression and anxiety disorders generally cluster into DSM-related categories.37 Twin and family studies largely provide evidence for the shared etiology of mood and anxiety disorder symptoms in young people in line with the unidimensional conceptualization.38- 43 Importantly, most of these studies encompass broad age ranges spanning childhood and adolescence, thus the associations at specific developmental stages remain unknown.
Age effects are essential to consider given that depression and anxiety disorders are characterized by different ages at onset2,5,44 and have developmentally dynamic etiologies. Environmental influences tend to decrease, while heritability increases with age and genetic innovation and attenuation take place at multiple stages.12,13,45- 50 Furthermore, depression may differ substantially pre-adolescence and postadolescence,51- 57 with 1 study finding that only the latter shares genetic influences with anxiety disorders.40 Thus, it is plausible that despite continuing comorbidity of internalizing problems, the genetic structure changes during development.5
The present analyses examined these important taxonomic issues by using a genetically informed design to investigate the structure of internalizing psychopathology cross-sectionally at multiple ages: childhood, adolescence, and early adulthood. To our knowledge, this is the first study to combine 5 waves of phenotypic and genetic data on depression symptoms and 4 anxiety subscales—generalized anxiety disorder, panic, separation anxiety, and social phobia symptoms—to address this question from a developmental perspective. The genetic structures of internalizing symptoms were investigated using 3 alternative models based on previous research: DSM-5–related structure and unidimensional and bidimensional (fear and distress) models. Because of the mixed findings and broad age ranges of previous research, the current study tested alternative models in an exploratory manner.
The analyses used data from 2 longitudinal twin studies: waves 1 and 2 from the Emotions, Cognitions, Heredity and Outcome Study (ECHO, child twin sample) and waves 2 through 4 from the Genesis 12-19 Study (adolescent/young adult twin and sibling sample). Full recruitment details are provided elsewhere58,59 (eAppendix in the Supplement). The studies were given ethical approval by the research ethics committees of the Institute of Psychiatry, King’s College London, South London and Maudsley NHS Trust and of Goldsmiths University, London. Written informed consent was obtained from parents of children younger than 16 years and from adolescents older than 16 years. Sample characteristics are presented in Table 1.
Child participants completed the Children’s Depression Inventory,62 a 27-item self-report questionnaire examining affective, cognitive, and behavioral signs of current depression. Adolescents and young adults completed the Short Mood and Feelings Questionnaire,63 a 13-item self-report measure assessing how often depressive symptoms occurred in the previous 2 weeks. Responses were summed to give total depression scores. Both measures demonstrate good reliability and validity.62,63
Children’s anxiety disorder symptoms were measured using the Screen for Child Anxiety Related Emotional Disorders.64 Children indicated how often in the last 3 months they experienced symptoms described by 41 questionnaire items. The adolescents completed the Spence Children’s Anxiety Scale,65 a 38-item self-report questionnaire tapping common anxiety symptoms. Adults completed the Revised Symptoms of Anxiety Scale,66 an age-appropriate version of the Revised Child Anxiety and Depression Scale,67 consisting of 36 self-report items designed to assess DSM-IV anxiety and depressive disorder symptoms. Responses were summed to create 4 DSM-IV–related anxiety subscale scores: generalized anxiety, panic/somatic symptoms, separation anxiety, and social anxiety. All measures have sound psychometric properties.64- 67
The internal consistencies and descriptive statistics of all measures are presented in Table 1.
Descriptive statistics were conducted using Stata (StataCorp).68 The associations between depression and anxiety subtypes were explored using full and partial correlations. For example, to investigate the unique association between depression and generalized anxiety symptoms, the scores on all other anxiety scales were controlled. This tested associations over and above the relationships with other variables that might confound the association owing to high covariance.
The twin design compares the similarity between monozygotic (sharing 100% of their genes) and dizygotic (sharing on average 50% of their segregating genes) twin pairs. Relative differences in within-pair correlations allow estimations of the influences of additive genetics, shared environment, and nonshared environment. Quantitative genetic methods are described comprehensively elsewhere.69
Models were fitted using OpenMx70 within R,71 a structural equation modeling package for the analysis of genetically informative data. Sampling weights were incorporated into child analyses, although they did not influence the results in a manner that would alter interpretation.72 The weight controls for biases due to selection criteria. Lower weights were assigned to individuals from categories overrepresented in the sample and higher weights to individuals from categories underrepresented relative to the population distribution. As is standard in model fitting analysis, variables were regressed for age and sex,73 and any with skew greater than 1 were transformed.
Univariate genetic analyses were conducted on all variables at each wave. Owing to sample size, sex differences were only examined in G1912. Scalar sex differences that examine whether males and females showed differences in variance were tested. A scalar model was fitted in twin modeling analyses for all variables except for social phobia (for which there was no difference in variance between males and females). Quantitative sex differences were tested to see whether males and females differ in magnitude of genetic and environmental influences but such differences were not found.
Three multivariate models that test different genetic structures underpinning associations between depression and anxiety subscales were fitted. They are discussed in the following order: DSM-5–related, unidimensional, and bidimensional (fear and distress) structures. The first model was a correlated factors solution (Figure, A), which is in line with the DSM-5 conceptualization in which each disorder is classified independently but expected to correlate with other disorders. This model includes additive genetic, shared environmental, and nonshared environmental influences on each of the scales and tests whether the correlation between them is due to correlations among the genetic and environmental factors that influence each of them. Each set of influences is allowed to correlate with one another. As such, the correlation among the variables can be mediated via genetic or environmental routes.
The second model was a 1-factor independent pathway model (Figure, B). This model reflects the unidimensional conceptualization by allowing internalizing disorder symptoms to share common genetic and environmental influences. It tests whether there is a single set of common etiologic factors that influence depression and all anxiety subscales, accounting for their correlations, in addition to variable-specific factors. The model includes 1 set of common additive genetic, shared environmental, and nonshared environmental factors that influence each of the measured variables.
The third model was a 2-factor independent pathway model (Figure, C). This model is similar to the 1-factor independent pathway model but contains a second common genetic factor loading on the anxiety symptoms hypothesized to belong to the fear factor. This model reflects the bidimensional conceptualization and tests whether there are 2 common genetic factors (distress and fear) and 1 common nonshared environmental factor that influences all variables, accounting for their correlation, in addition to variable specific factors.
Models were fitted using raw data maximum likelihood. The core fit statistic was minus twice the log likelihood of the observations. This is not an overall measure of fit but provides a relative measure of fit because differences in minus twice the log likelihood between models are distributed as χ2. Therefore, to examine the overall fit of the genetic model, we compared the minus twice the log likelihood with that of a saturated model (one which fully describes data using the maximum number of free parameters, estimating variances, covariances, and means for the raw data to get a baseline index of fit). The fit of submodels was assessed by χ2 difference tests, the Akaike Information Criterion (AIC), and the Bayesian Information Criterion (BIC) (AIC = χ2 – 2df; BIC = χ2 – kln[n]), with lower χ2 values and more negative AIC and BIC values suggesting a better fit. If the difference between the AIC of 2 models was less than 10, the more parsimonious model was selected.74 Independent pathway models are nested in the correlated factors solution, and the 1-factor independent pathway model is nested in the 2-factor independent pathway model. Information about the precision of parameter estimates was obtained by likelihood-based confidence intervals. The analyses were repeated excluding siblings to narrow the age ranges (eTable 1 in the Supplement) and including an additional anxiety subscale: fear of physical injury (only available at the 2 adolescent time points; eTable 2 in the Supplement).
The results focused on the association between depression and the different anxiety subscales. The phenotypic and genetic associations among the anxiety subscales are presented elsewhere.72,75
Full correlations at all ages showed that depression symptoms were significantly associated with all anxiety subscales (Table 2). In childhood and adolescence, depression symptoms showed significantly stronger correlations with generalized anxiety symptoms (r = 0.36 to 0.60) than with all other subscales except for panic symptoms (r = 0.28 to 0.57).
Partial correlations that controlled for all other variables within times are shown in Table 2. In childhood, when controlling for concurrent associations, depression symptoms were only significantly associated with generalized anxiety symptoms (r = 0.21 and 0.20). At 15 years, partial correlations revealed that depression symptoms were significantly associated with 3 anxiety subscales: generalized anxiety (r = 0.19), panic (r = 0.24), and social phobia (r = 0.14) symptoms. At a mean age of 17 years and in young adulthood, depression symptoms were significantly associated with all anxiety subscales even when controlling for concurrent associations.
Univariate analyses revealed that genetic influences on depression and anxiety symptoms were generally small to moderate, shared environmental influences were small and nonsignificant, and nonshared environmental influences were large (eTable 3 in the Supplement). Multivariate model fitting results are presented in Table 3. Shared environmental influences were nonsignificant and were dropped from the models without a significant deterioration of the fit in adolescence and young adulthood; fit statistics and parameter estimates are therefore presented for models with additive genetic and nonshared environmental influences.
In childhood, the most restrictive 1-factor independent pathway model was the best fitting model (Table 4). The common genetic factor accounted for most of the genetic influences on all anxiety subscales but did not contribute to genetic variance in depression symptoms, which instead was influenced by unique genetic influences. There were moderate to large unique nonshared environmental influences on each symptom.
In adolescence, the least restrictive model, the correlated factors solution, showed the best fit to the data in line with DSM-5 conceptualization (Table 5). Genetic correlations were mostly large. Depression symptoms generally had higher genetic correlations with generalized anxiety (r = 0.71 and 0.74), panic (r = 0.78 and 0.61), and social phobia (r = 0.66 and 0.53) than with separation anxiety (r = 0.52 and 0.15) symptoms. Nonshared environmental correlations were generally moderate. Genetic influences explained a substantial proportion of the phenotypic correlation between depression and anxiety subscales (36% to 100%).
In young adulthood, a 2-factor independent pathway model showed the best fit to the data in line with a bidimensional conceptualization (Table 6). The first common genetic factor loaded significantly on all variables and accounted for most of the genetic variance. The second common genetic factor, specified to load on the fear variables, showed small but significant contributions to panic, separation anxiety, and social phobia symptoms. In addition, depression and generalized anxiety symptoms had significant unique genetic influences. The common nonshared environmental factor loaded significantly on all variables but there were also significant unique nonshared environmental influences on each variable.
To our knowledge, this study is the first to investigate the phenotypic and genetic structure of internalizing disorder symptoms at 3 developmental stages. The results provide preliminary evidence for developmental differences in the associations between depression and multiple anxiety disorder symptoms, advancing the search for an evidence-based conceptualization of internalizing disorders in diagnostic manuals.
We observed different etiologic structures of internalizing disorder symptoms at 3 developmental phases, with common genetic vulnerability across depression and anxiety disorder symptoms only emerging in adolescence. Specifically, in childhood, when controlling for concurrent associations, only the generalized anxiety disorder symptoms were associated with depression. Furthermore, childhood depression was influenced by separate genetic factors from the anxiety subscales. In adolescence, comorbidity began to increase—partial correlations revealed that at 15 years of age, depression was associated with 3 anxiety disorder subscales: generalized anxiety disorder, panic, and social phobia symptoms. At this developmental stage, the etiologic structure reflected the DSM-5 conceptualization of distinct but correlated disorders in contrast to previous studies that found support for unidimensional or bidimensional latent factor structures in young people.32- 36 These age differences may be explained by anxiety emerging in childhood, while depression peaks in adolescence,2,36and are in agreement with previous studies finding that depression pre-adolescence and postadolescence may differ substantially,51- 57 which could be explained by significant new genetic influences coming online after puberty.12,13,45,46,48
In young adulthood, these associations broadened even further, and depression was significantly correlated with all anxiety disorder symptom scales. Genetic analyses provided support for both unidimensional19- 22 and bidimensional26- 30 conceptualizations of internalizing psychopathology. The 2 genetic factors representing distress and fear emerged, although the genetic fear factor had a relatively small influence on the fear symptoms. The current results add to a debate as to whether generalized anxiety disorder ought to be classified together with depression,14- 17 and they suggest that at most ages, generalized anxiety disorder symptoms are no more closely related to depression than other anxiety subtypes. The exception is childhood, where the generalized anxiety disorder symptom subscale was the only one associated with depression, although this association was not underpinned by shared genes.
While genetic influences accounted for comorbidity, in agreement with the generalist genes hypothesis,76 the nonshared environment was largely symptom specific across development, accounting for most of the unique variance that makes each disorder symptom a discrete condition. These results carry implications for the molecular genetic studies of depression and anxiety, which in turn may inform clinical interventions.77- 79 The results provide preliminary support for broadening phenotypic definitions in linkage or association studies, as including adult cases with a variety of internalizing disorders underpinned by an overarching genetic internalizing factor would lead to increasing power to detect shared susceptibility loci.80 Conversely, the difference in the genetic results pre-adolescence and postadolescence also provides a preliminary argument for narrowing the phenotypic definitions by age.81
A key clinical implication of our findings was the support for transdiagnostic treatment approaches for anxiety and depression disorders, which are designed to target common elements of several disorders in 1 protocol.82- 87 The developmental pattern of the data suggests that while disorder-specific treatment may be more appropriate for pediatric patients, treatment focused on a range of symptoms common to internalizing disorders may be more appropriate for older patients. The evidence for a shared genetic etiologic factor is in agreement with the findings that internalizing disorders respond to similar interventions and therapies.23,86- 91
The genetically informative, representative samples and multiple time points were strengths of the current study. However, a number of limitations are noteworthy. First, the child sample was smaller than the adolescent/adult sample. Although considered large for phenotypic analyses, the child sample had reduced power to examine sex differences or shared environmental influences, and parameter estimates had large confidence intervals. Replication in larger pediatric twin samples is essential. However, because of the internal replication of results across the 2 time points, interpretations seem broadly applicable for childhood. Second, the inclusion of siblings in the Genesis 12-19 Study meant there were large age ranges in adolescence and early adulthood. However, 72% of the participants were twins, and additional analyses excluding siblings suggest that the results are applicable to tighter age ranges. Third, to inform understanding of comorbidity of internalizing disorders in clinical settings, the results should be replicated in clinical samples with comorbid diagnoses and using lifetime diagnostic interviews. However, internalizing symptoms are important markers of psychopathology92- 94 and because common mental disorders are quantitative traits,95 there is evidence that differently defined internalizing problems have the same etiology.8,96,97 Fourth, our study included self-report measures, allowing comparisons across waves. While studies have shown that young children can accurately report on their own internalizing symptoms,98,99 including parent-report measures at these waves may have strengthened our findings. Last, there are limitations inherent to the twin design, discussed comprehensively elsewhere.100 These have minimal and contrasting effects on parameter estimates that should therefore be taken as indicative rather than absolute.
Our results suggest that the phenotypic and genetic structure of internalizing disorder symptoms may differ across development. Depression and anxiety seem to be somewhat distinct in childhood but become more associated and share most of their genetic etiology from adolescence, with an overarching internalizing genetic factor emerging in early adulthood. The results have multiple implications for further research, taxonomy, and clinical practice. They affirm the need to continue examining developmental differences in the etiology of mood and anxiety disorders to ensure that the diagnostic conceptualization of psychopathology is age appropriate.
Corresponding Author: Thalia C. Eley, PhD, MRC Social, Genetic, and Developmental Psychiatry Centre, PO Box 80, Institute of Psychiatry, King’s College London, De’Crespigny Park, London, SE5 8AF, England (email@example.com).
Submitted for Publication: November 25, 2013; final revision received February 21, 2014; accepted March 24, 2014.
Published Online: June 11, 2014. doi:10.1001/jamapsychiatry.2014.655.
Author Contributions: Ms Waszczuk had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. Ms Waszczuk and Dr Zavos were the joint lead authors.
Study concept and design: All authors.
Acquisition, analysis, or interpretation of data: Waszczuk, Zavos, Eley.
Drafting of the manuscript: Waszczuk.
Critical revision of the manuscript for important intellectual content: Zavos, Gregory, Eley.
Statistical analysis: Waszczuk, Zavos.
Obtained funding: Gregory, Eley.
Study supervision: Zavos, Eley.
Conflict of Interest Disclosures: None reported.
Funding/Support: Both the Emotions, Cognitions, Heredity and Outcome and Genesis 12-19 studies were supported by a Medical Research Council Training Fellowship and a Career Development Award given to Prof Eley. Waves 1 through 3 of the Genesis 12-19 Study were also funded by the W. T. Grant Foundation, the University of London Central Research fund, and wave 4 was supported by the Economic and Social Research Council (RES-000-22-2206), the Institute of Social Psychiatry, and a Leverhulme Research Fellowship awarded to Dr Gregory. Ms Waszczuk was supported by a PhD studentship awarded by the Alexander von Humboldt Foundation. This study presents independent research partly funded by the National Institute for Health Research (NIHR) Biomedical Research Centre at South London and Maudsley NHS Foundation Trust and King’s College London.
Role of the Sponsor: The funders 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.
Disclaimer: The views expressed in this article are those of the authors and not necessarily those of the NHS, the NIHR, or the Department of Health.
Additional Contributions: We thank the families for their participation as well as numerous staff and students from the Social, Genetic, and Developmental Psychiatry Centre, Institute of Psychiatry, London and Goldsmiths, University of London.