Figure 1. Full multivariate model for 6 anxiety disorders in men and women. Latent variables are indicated by circles and observed variables by rectangles. A1 and A2 are the first and second common additive genetic factors; C1 and C2, the first and second common shared (familial) environmental factors; and E1 and E2, the first and second common unique (individual) environmental factors. ASP, CSP, and ESP refer to disorder-specific additive genetic, shared environmental, and unique environmental factors, respectively. Path coefficients on the arrows, representing standardized regression coefficients, indicate the strength of the influence of the latent variable the arrow originates from on the observed variable it ends on. GAD indicates generalized anxiety disorder.
Figure 2. Best-fitting model for 6 anxiety disorders. The model contains 2 common additive genetic factors (A1 and A2), 1 common shared environmental factor (C1), and 1 common unique environmental factor (E1). The only disorder-specific factors that were significant and were retained in the model were an agoraphobia-specific genetic factor, ASP, plus the unique environmental disorder-specific factors, Esp. GAD indicates generalized anxiety disorder.
Hettema JM, Prescott CA, Myers JM, Neale MC, Kendler KS. The Structure of Genetic and Environmental Risk Factors for Anxiety Disorders in Men and Women. Arch Gen Psychiatry. 2005;62(2):182–189. doi:10.1001/archpsyc.62.2.182
The anxiety disorders exhibit high levels of lifetime comorbidity with one another. Understanding the underlying causes of this comorbidity can provide insight into the etiology of the disorders and inform classification and treatment.
To explain anxiety disorder comorbidity by examining the structure of the underlying genetic and environmental risk factors.
Lifetime diagnoses for 6 anxiety disorders (generalized anxiety disorder, panic disorder, agoraphobia, social phobia, animal phobia, and situational phobia) were obtained during personal interviews from a population-based twin registry. Multivariate structural equation modeling that allowed for sex differences was performed.
General community sample.
More than 5000 members of male-male and female-female twin pairs from the Virginia Adult Twin Study of Psychiatric and Substance Use Disorders.
Main Outcome Measures
Parameter estimates for best-fitting model.
The full model, which contained 2 common genetic, shared environmental, and unique environmental factors plus disorder-specific factors, could be constrained to equality across male and female study participants. In the best-fitting model, the genetic influences on anxiety were best explained by 2 additive genetic factors common across the disorders. The first loaded most strongly in generalized anxiety disorder, panic disorder, and agoraphobia, whereas the second loaded primarily in the 2 specific phobias. Social phobia was intermediate in that it was influenced by both genetic factors. A small role for shared environmental influences was observed owing to a single common factor that accounted for less than 12% of the total variance for any disorder. Unique environmental influences could be explained by a single common factor plus disorder-specific effects.
The underlying structure of the genetic and environmental risk factors for the anxiety disorders is similar between men and women. Genes predispose to 2 broad groups of disorders dichotomized as panic-generalized-agoraphobic anxiety vs the specific phobias. The remaining associations between the disorders are largely explained by a unique environmental factor shared across the disorders and, to a lesser extent, a common shared environmental factor.
Comorbidity among the anxiety disorders has proved to be the rule rather than the exception,1- 3 such that it is unlikely for an individual to experience a single anxiety disorder in isolation.4 The etiology of comorbidity can be loosely divided into origins of causation (disorder A causes disorder B) or shared risk factors (both disorders are etiologically related to the same underlying causes).5 It is the latter that this study seeks to elucidate regarding the underlying risk factors that increase comorbidity among the anxiety disorders.
Familial aggregation that largely results from genetic risk factors has been documented for all of the major anxiety disorders.6 Although twin studies have suggested that genetic risk factors shared between anxiety disorders may account for excess rates of comorbidity, the results of family studies have been inconsistent. In one of the earliest twin studies to examine multiple anxiety disorders, Torgersen7 reported a 2:1 monozygotic-dizygotic concordance ratio for all anxiety disorders but found no monozygotic pairs in which twins had the same disorder, most likely because of the limited sample size. Prior twin studies in the current sample have found genetic and environmental factors common across the phobic disorders in men8 and women9 and among major depression, generalized anxiety disorder (GAD), panic disorder, and specific phobias in both sexes.10,11 Similarly, analyses from a large twin data set, the Vietnam Era Twin Registry, suggested significant genetic correlation between GAD and panic disorder12 and between these syndromes and posttraumatic stress disorder.13 Contrary to these findings, most family studies14- 19 of the anxiety disorders reported relative specificity in their familial aggregation.
Like genetic factors, the experience of stressful life events, either shared with other family members or unique to the individual, is etiologically related to the development of anxiety disorders and is likely to have nonspecific effects across disorders. For example, many forms of childhood adversity are associated with adult psychiatric outcomes with little specificity,20 including the anxiety disorders.21 Childhood sexual abuse has been found to elevate the risk for adult GAD and panic disorder, as well as other psychiatric and substance use disorders.22,23 Similarly, childhood parental loss or separation has been linked to various forms of adult psychopathologic characteristics.24- 26
In this study, we used multivariate twin analysis to examine the latent structure of the genetic and environmental risk factors that underlie 6 anxiety disorders. We attempted to answer the following questions: (1) How similar is the pattern of genetic and environmental risk factors across the anxiety disorders in men and women? (2) Can one identify shared risk factor domains (genetic, common familial environment, or unique individual environment) that account for the phenotypic correlations (comorbidity) among the anxiety disorders? (3) If shared risk factors exist, how are they structured (ie, how many separate factors are necessary to best describe patterns of comorbidity)? (4) What are the relative roles of risk factors that are common across the anxiety disorders vs those that are disorder specific?
The sample in this study derives from the population-based Virginia Adult Twin Study of Psychiatric and Substance Use Disorders.27 Female-female (FF) twin pairs, from birth years 1934 to 1974, became eligible if both members previously responded to a mailed questionnaire between 1987 and 1988, the response rate to which was approximately 64%. They were approached for 4 subsequent waves of personal interviews from 1988 to 1997, with individual cooperation rates ranging from 85% to 92%. Between waves 1 and 2, from which the current analyses derive, the pairwise cooperation rate was approximately 90% and the time between the 2 interviews was approximately 13 to 15 months. The male-male and male-female (MMMF) twin pairs, covering the birth years 1940 to 1974, were ascertained in a separate study, with an initial cooperation rate of 72.4%, and were approached for 2 waves of interviews from 1993 until 1998. Zygosity was determined by a combination of standard questions,28 photographs, and DNA analysis.27 The mean (SD) ages of the FF and MMMF samples at their final interviews were, respectively, 36.6 (8.1) years and 36.8 (9.1) years. Interviewers had a master’s degree in a mental health–related field or a bachelor’s degree in this area plus 2 years of clinical experience. At each wave, members of a twin pair were interviewed by different individuals who were blind to clinical information about the co-twin. In this report, we examine only same-sex twin pairs and will therefore refer to male-male pairs as the MM sample.
The phobias were diagnosed using an adaptation of DSM-III criteria,29 which required the presence of 1 or more of 22 fears that the respondent recognized as unreasonable and that, in the judgment of the interviewer, objectively interfered with the respondent’s life.8 They were grouped according to standard nomenclature as agoraphobia, social phobia, or specific phobia, the last having been further subtyped by the nature of the respondent’s fears (animal, situational, or blood injury). Since a low prevalence of GAD and panic disorder had been problematic in previous analyses,30,31 modified DSM-III-R32 diagnostic criteria were used for these conditions. We adopted a broad diagnostic approach to these 2 disorders, reducing the minimum duration from 6 months to 1 month for GAD and requiring a history of panic attacks that met at least 2 symptomatic criteria, peaking within 30 minutes for panic disorder. We have shown that these approaches reflect the same continuum of liability as the fully syndromal disorders.30,31
The diagnostic data used in these analyses came from various waves of the 2 samples. In the FF sample, all diagnoses came from the first wave (2156 participants) except for panic disorder, which was derived from wave 2 (1999 participants) approximately 12 to 15 months later. The pairwise cooperation rate was approximately 90% between the 2 interviews. For the MM sample, all diagnoses were obtained from the second wave (2939 participants). Comorbidity in these analyses refers to the lifetime occurrence of 2 or more disorders within an individual.
We have outlined elsewhere our approach to multivariate genetic analysis.9,33 The models described herein are based on a liability threshold model, the strengths and limitations of which have been previously discussed.33,34 Like traditional factor analysis, multivariate genetic analysis seeks to explain covariation among multiple variables with a few factors. However, traditional exploratory factor analysis is purely descriptive, whereas multivariate genetic analysis provides insight into the sources of resemblance.
Using the software package Mx,35 we fit models by the method of maximum likelihood to data from all individual twins, including those with missing data or without an interviewed co-twin. This method can reduce the impact of cooperation bias if the data are “missing at random” as described by Little and Rubin.36
For these analyses, we used independent pathway models that contained 2 genetic common (A1, A2), 2 shared environmental common (C1, C2), and 2 unique environmental common (E1, E2) factors in addition to disorder-specific factors. Two factors were used based on indications from other analyses that phobias may be etiologically distinct from the other anxiety disorders11 and because of computational limitations in trying to model more than 2 factors. Also, to keep computer run-times tractable, we restricted our analyses to the simultaneous modeling of 6 disorders, including data from only 2 of 3 available specific phobia subtypes for these analyses.
We fit this model to both the MM and FF data sets, first allowing all parameters, including thresholds, to vary by sex. We then constrained all of the path estimates (but not the thresholds) to be equal in the 2 sexes. We did this for all of the path estimates en masse to test the global structural equivalence across sexes rather than test each individual pathway. We fit subsequent, simpler submodels with the goal of obtaining an overall best-fitting model that best explained the data by taking maximal account of risk factors shared between the anxiety disorders. To minimize the overall number of tests and therefore type I error, we compared submodels created by simplifications across classes of factors (common vs disorder specific) rather than testing the significance of each of the individual path loadings. Twice the difference in log likelihood between a higher order and submodel yields a statistic that is asymptotically distributed as χ2 with the df equal to the difference in their number of parameters. We used the Akaike information criterion (AIC)37 for model selection. The lower its value, the better the balance between explanatory power and parsimony.
Table 1 gives the lifetime prevalence rates for the 6 anxiety disorders in our sample of male and female twins. Comorbidity among the disorders is expressed as both tetrachoric correlations and odds ratios (taking into account the clustered nature of twin data) separately for male and female participants. Tetrachoric correlation refers to the correlation between the liability for the 2 disorders within a participant, as determined by the liability threshold model. All of the values except 1 are significant at P<.05. The overall pattern is similar across the sexes. The strongest associations are for agoraphobia with social phobia, GAD, and panic, with more modest associations between these 3 and the specific phobias, which are reflected in the factor structure described herein.
The results of our model-fitting process are given in Table 2. Models 1 and 2 are full 2-factor independent pathway models, which we will refer to as 2-2-2 models, with the digits indicating, respectively, the number of common (across anxiety disorders) genetic, shared (within family) environmental, and unique (individual) environmental factors (Figure 1). These models also include factors that affect the individual anxiety disorders in a disorder-specific fashion. Model 1 allows for the effects of the various factors to differ between men and women, whereas model 2 constrains these to be equal across the sexes. Model 2 provides a better balance of parsimony and explanatory power by lower AIC (by 52.3 units) and is not significantly different from model 1 by a χ2 difference test (Δχ245 = 37.74, P = .77), suggesting that the pattern of genetic and environmental risk factors does not differ significantly between men and women, as had been seen in prior univariate analyses of these disorders. All subsequent models will therefore ignore sex differences in the effects of the latent factors on the disorders.
Models 3a and 3b attempted to simplify the factor structure by removing the effects of disorder-specific additive genetics or shared environment. Both provided an improvement in parsimony via lower AIC compared with model 2. However, model 4a, in which both were removed simultaneously, provided the most parsimonious model with the full 2-2-2 structure for the common factors. However, given the substantial disorder-specific genetic loading on agoraphobia in the full model, we tested its significance separately by including it in model 4b. A comparison between models 4a and 4b shows a significant effect of this factor (Δχ21 = 4.01, P = .04), which was retained in subsequent submodels.
Models 5a (1-2-2), 5b (2-1-2), and 5c (2-2-1) tested the hypotheses that 2 vs 1 additive genetic, shared environmental, or unique environmental common factors, respectively, are required to best explain the pattern of comorbidity between the anxiety disorders. Model 5a, in which the effects of A2 were constrained to 0 but the effects of C2 and E2 were retained, provided a slightly improved AIC but a nearly significant deterioration in fit as indicated by the χ2 difference (P = .11). Model 5c (2-2-1) provided the lowest AIC of this class of models, suggesting that the unique environmental risk factors can be effectively explained using 1 common factor (E1) plus disorder-specific factors.
Models 6a and 6b compared the 2 next simpler submodels (1-2-1 vs 2-1-1) with model 5c. The removal of the effects of A2 in model 6a produced a significant deterioration in fit (P = .04) and a higher (less negative) AIC, suggesting that 2 additive genetic factors common to all the anxiety disorders best explain the data. Conversely, the improvement in AIC provided by model 6b compared with 5c suggests that 1 shared environmental risk factor common across disorders (vs 2) can adequately explain the data.
Further simplifications to the underlying risk structure beyond model 6b did not produce better fit to the data. Models 7a and 7b, which tested the 2-0-1 and 2-1-0 structures against model 6b, provided worse fit by AIC and could be rejected by the χ2 difference test. Model 8, which again tested the significance of the A2 factor but within the simplified factor structure background, produced a severe deterioration in fit by the χ2 difference statistic and thus could also be rejected with confidence.
Figure 1 and Figure 2 depict the path estimates for the full model for men and women (model 2) and the best-fitting model (model 6b). As can be seen, the overall pattern of loadings is similar between these models, although model 6b provides a substantially improved balance of parsimony and explanatory power and is more easily interpretable. The rest of this section will focus on these latter results. The proportions of variance in liability for the 6 anxiety disorders due to the factors included in this best-fitting model are given in Table 3.
Figure 2 shows that A1 substantially affects GAD, panic, agoraphobia, and to a lesser extent social phobia, whereas A2 primarily affects the 2 specific phobia subtypes with less substantial loadings in GAD, panic, and social phobia. This is reflected in Table 3, where A1 accounts for only approximately 1% of the total variance of the animal and situational phobias, and A2 accounts for less than 5% of the variance of the first 4 disorders. Agoraphobia alone was found to have a significant proportion of its genetic risk not shared with the other disorders. The total genetic proportion of variance (that is, heritability) was estimated to be approximately 25% to 35% for all disorders except social phobia, in which it accounted for only 10% of the total variance in our data.
The common shared environmental factor, C1, has loadings of 0.33 or less, accounting for approximately 10% or less of the variance for any of the disorders (Table 3). The common unique environmental factor, E1, loaded most strongly on agoraphobia but had substantial loadings across all of the disorders except social phobia. Disorder-specific effects were generally higher than those arising from E1 except for agoraphobia.
We applied multivariate structural equation modeling to lifetime diagnoses of 6 anxiety disorders assessed in male and female twins to examine the pattern of genetic and environmental risk factors that underlie their observed comorbidity. In this effort, we sought to answer a series of 4 questions, the last 3 of which will be consolidated for the purposes of discussion.
We tested whether the impact of genetic and environmental risk factors on the anxiety disorders and their patterns of association differed between men and women by constraining the path loadings to be equal across the sexes. This produced a nonsignificant deterioration in model fit and a more parsimonious explanation of the data. This is consistent with prior analyses of the individual anxiety disorders in this sample. This finding suggests that despite the nearly 2-fold greater prevalence rate in women, the same underlying liability structure broadly accounts for patterns of comorbidity among the anxiety disorders, independent of sex.
We sought to account for the observed comorbidity among the anxiety disorders by starting with a full model that contained 2 factors common to all of the anxiety disorders in each of the risk factor domains (additive genetic, shared environmental, and unique environmental effects), with residual variance explained by disorder-specific factors. By testing successively simpler submodels, we could most parsimoniously explain this comorbidity with a model that contained 4 factors common across the 6 disorders: 2 genetic factors, a single shared environmental factor, and a single unique environmental factor. In this best-fitting model, residual effects were limited (except for agoraphobia-specific genetics) to disorder-specific unique environmental factors that include the effects of measurement error.
The first genetic factor had highest loadings on GAD, panic, agoraphobia, and to a lesser extent, social phobia. A shared genetic diathesis among these disorders is not surprising heuristically and may help to explain their common response to antidepressant medications. However, residual differences among these disorders may be mediated by the relative impact of the second genetic factor and differential effects of the environmental risk factors. One might hypothesize a set of genes that increases the liability of developing a common intermediate phenotype, such as an anxious personality trait, which then interacts with modifier genes or stressful life events to surpass the threshold for the development of 1 or more disorders. Similarly, the finding that the specific phobias load primarily on a second genetic factor uncorrelated with the first suggests that their genetic etiology may be largely distinct from the other disorders. Interestingly, specific phobias do not share the same treatment response to medications as the other disorders studied.
These results have potential implications for the conceptual unity of the anxiety disorders within psychiatric nomenclature. Although heterogeneous in their symptom presentation, age at onset, course of illness, and biological correlates, the results of this study suggest that etiologically, the anxiety disorders possess a relatively simple genetic architecture. Our analyses place agoraphobia and, to a lesser extent, social phobia closer etiologically to GAD and panic disorder than to the specific phobias. Given the way that agoraphobia is defined in relation to panic disorder, this is not surprising. For social phobia, the decision to place it in a separate category from the specific phobias is in part supported by these analyses. In addition, although GAD and panic disorder were created as separate disorders out of the former category of anxiety neurosis, their risk structure is remarkably similar.
Because of computational limitations, we could simultaneously model only 6 disorders and thus did not include all available specific phobia subtypes. To examine the effect of omitting blood-injury phobia on the present results and also understand how it fits into the risk structure, we analyzed a model that substituted it for social phobia. The results place blood-injury phobia closer in its genetic risk structure to agoraphobia than to the specific phobias (details available on request from the authors) and suggest that it is the genetic risk factors for situational and animal phobia that are clearly distinguishable from those of the other anxiety disorders.
It is interesting to compare these findings with another recent multivariate analysis in this sample that compared the risk structure of internalizing vs externalizing disorders.11 A secondary analysis in that report focused on a subset of 5 internalizing disorders: major depression, GAD, panic disorder, and situational and animal phobias. That study also reported a broadly similar 2-factor genetic architecture in which the specific phobias loaded primarily in a factor that was distinct from the other disorders. This suggests that the overall factor structure is robust to the specific set of internalizing disorders examined within our data set and is therefore unlikely a result obtained by chance. One potential difference was that in the present study we used an orthogonal rotation of the factor structure with uncorrelated latent factors for our modeling, whereas the former analysis applied an oblique rotation that allowed a nonzero correlation between the common factor pairs. The former is conceptually clearer but may not be the best representation of reality. When we allowed the factors in the present study’s full model to be correlated via an oblique rotation instead, the correlation between the 2 common genetic factors was estimated at less than 10%, which supports their relative independence. However, the 2 unique environmental common factors, E1 and E2, were found to have approximately 50% correlation. Although these estimates differ somewhat from the prior analysis, they are unlikely significantly different, and any real differences probably arise from using different sets of phenotypes. We cannot compare the 2 studies’ best-fitting models, however, because our prior study focused on the broad factor structure and did not explore results from more parsimonious submodels.
The findings of this analysis have important implications for studies that seek to find the genes that increase liability for the anxiety disorders. For monogenic illnesses, genetic studies have traditionally emphasized the importance of selecting a “pure” phenotype to maximize the power of detecting the effect of the causative gene. To do so may be more difficult (and likely less appropriate) for complex, multifactorial disorders that exhibit nonmendelian inheritance and high rates of comorbidity. Gene-finding studies for the anxiety disorders may benefit from a multivariate approach using a broadened phenotype that includes a range of disorders with shared genetic risk factors. For example, a linkage study38 that used such a broadened phenotype of panic disorder together with other anxiety disorders has reported significant linkage to chromosome 9q31, whereas a family study39 that included panic disorder, bladder problems, severe headaches, mitral valve prolapse, and thyroid conditions found significant linkage on chromosomes 13q. Further developments in gene finding across disorders may be expected from multivariate linkage analysis.40
The results of this analysis should be interpreted in the context of several potential limitations. First, these results are predicated on the basic tenets that underlie the twin model, such as the equal environment assumption, which have been discussed in detail elsewhere.41 Second, we used broadened diagnostic definitions for GAD and panic disorder to offset modeling limitations that arise owing to their low prevalence. Although univariate analyses of these disorders have suggested that these approaches reflect the same continuum of liability as the fully syndromal disorders,30,31 it is unclear whether this is true of the multivariate distribution of these liabilities. Third, lifetime history of each anxiety disorder was assessed at 1 time point, which potentially confounds effects of unique environment and measurement error, reducing the estimates of additive genetic and shared environmental parameters accordingly. We have found, for phobias in our FF sample, that improving the diagnostic reliability by reducing error via multiple, sequential assessments substantially increased the estimate of additive genetic effects.42 Fourth, we assessed panic disorder in the female twins in a separate interview from the other disorders and thus treated it differently from the male twins. To assess whether ascertainment bias that resulted from this difference affected our results, we performed a logistic regression to determine whether panic symptoms at the first interview predicted cooperation at the second, controlling for other relevant predictors such as zygosity and age. We found no significant effects, suggesting little impact from this difference. Fifth, we could not analyze the entire spectrum of anxiety disorders as classified in the DSM because we did not assess obsessive-compulsive disorder or posttraumatic stress disorder in this sample. Sixth, our twin model is based on the assumption that shared risk factors underlie phenotypic comorbidity. As mentioned previously, other mechanisms are possible5 but were not tested. Furthermore, our multivariate model is based on the hypothesis that the 3 risk factor domains have different effects on the pattern of comorbidity (independent pathway model), as opposed to the hypothesis that genes and environment act on disorder comorbidity via a common pathway.43 We tested this common pathway model, but it did not produce results that fit the data as well as the model presented herein (details available on request from the authors). Finally, because the sample was made up entirely of white participants, the results may not generalize to other ethnic groups.
The underlying structure of genetic and environmental risk factors for the anxiety disorders does not significantly differ between men and women. Two genetic factors predispose to 2 broad groups of disorders dichotomized as generalized and panic anxiety plus agoraphobia vs the specific phobias, with social phobia intermediate between these. The remaining associations between the disorders are largely explained by a unique environmental factor shared across the disorders and, to a lesser extent, a common shared environmental factor. In the most parsimonious model, the only disorder-specific factors that contribute significantly are from an agoraphobia-specific genetic factor and unique environmental effects for each disorder. This suggests that the underlying patterns of liability can be broadly understood as follows: individual genetic factors derive from 2 sets of genes that increase risk for the 2 classes of panic-generalized-agoraphobic anxiety vs specific phobias. Risk across all of the anxiety disorders may be further increased by life experiences either shared with other family members or unique to the individual, with varying impact (loadings) depending on the disorder. Add to this a set of unique environmental factors that increase the liability for one anxiety disorder independent of another. Although these last factors appear to carry the largest proportion of variance, this is likely an overestimate due to confounding by measurement error, with a compensatory underestimate of the true effects of genes and/or shared environment.
Correspondence: John M. Hettema, MD, PhD, Department of Psychiatry, PO Box 980126, Richmond, VA 23298-0126 (firstname.lastname@example.org).
Submitted for Publication: December 19, 2003; final revision received June 14, 2004; accepted July 9, 2004.
Funding/Support: This study was supported by National Institutes of Health (Bethesda, Md) grants MH-40828, MH-65322, MH-20030, DA-11287, MH/AA/DA-49492, and K08 MH-66277-1. The Mid-Atlantic Twin Registry (MATR), directed by Judy Silberg, PhD, and Lindon Eaves, PhD, has received support from the National Institutes of Health; the Carman Trust, Richmond, Va; the WM Keck Foundation, Los Angeles, Calif; the John Templeton Foundation, Radnor, Pa; and the Robert Wood Johnson Foundation, Princeton, NJ.
Acknowledgment: We acknowledge the contribution of the Virginia Twin Registry, now part of the MATR, for ascertainment of participants for this study.