Abstract
Multiple executive functions (EFs) are associated with trait anxiety and depression symptoms, but it is unclear how genetic and/or environmental factors account for these associations and whether they are explained by general variance underlying multiple EFs (i.e., common EF). In this study, 1,207 male twins completed seven EF tasks and measures of trait anxiety and depression symptoms at an average age of 62 years. The common EF factor was associated with both anxiety (r = −.24) and depression (r = −.32) symptoms. Anxiety and depression had a near-complete genetic overlap (rg = .96). Genetic influences shared with depression accounted for 83% of the phenotypic correlation between anxiety and common EF. For depression, the phenotypic correlation with common EF was explained by genetic influences shared with anxiety (69%) but also nonshared environmental influences distinct from anxiety (20%). These results suggest that genetic influences on common EF play a role in the observed genetic overlap between anxiety and depression.
A large body of research has examined the association between individual differences in anxiety and depression symptoms and executive functions (EFs), cognitive abilities that control and regulate behavior (Miyake & Friedman, 2012). According to attentional control theory, for example, subclinical levels of anxiety are associated with poor EFs in multiple contexts, including prepotent response inhibition, task-set shifting, and the updating of working memory (WM) representations (WM updating; Berggren & Derakshan, 2013; Eysenck & Derakshan, 2011; Eysenck, Derakshan, Santos, & Calvo, 2007). Associations between EFs and depression symptoms have also been observed in the general population (Franz et al., 2011; Holmes & Pizzagalli, 2007).
Despite the breadth of evidence for associations between EFs and tendencies to experience anxiety and depression symptoms, these findings have yet to be integrated into the unity and diversity framework of EF that emphasizes its genetic/environmental variance components and the general processes underlying many EF abilities (i.e., common EF; Friedman et al., 2016; Miyake & Friedman, 2012). The strong correlation between trait anxiety and depression symptoms is explained largely by overlapping genetic influences (Eley & Stevenson, 1999; Jardine, Martin, Henderson, & Rao, 1984; Kendler, Heath, Martin, & Eaves, 1987), and the associations between EFs and other clinically relevant constructs such as substance use, conduct disorder, and attention problems also appear to be due to genetic influences (Gustavson, Stallings, et al., 2017; Young et al., 2009). However, little research has examined whether genetic influences account for the associations between EFs and trait anxiety or depression symptoms (Franz et al., 2011; Friedman, du Pont, Corley, & Hewitt, 2018; Routledge et al., 2017), especially the overlapping genetic influences between anxiety and depression. Moreover, the associations between anxiety and separate EFs may reflect an underlying association with common EF, rather than distinct associations with individual EF processes. The current study examined associations between EFs and anxiety and depression symptoms in a genetically informative study of community-dwelling middle-age male twins in the Vietnam Era Twin Study of Aging (VETSA).
Common Variance in Anxiety and Depression Symptoms
Anxiety and depression exhibit substantial overlap in the general population (Caci, Bayle, Dossios, Robert, & Boyer, 2003; Endler, Cox, Parker, & Bagby, 1992). Much of the phenotypic correlation between levels of anxiety and depression symptoms in the general population is explained by genetic influences, regardless of whether symptoms are measured with diagnostic interviews (Hink et al., 2013; Krueger, McGue, & Iacono, 2001; Tackett et al., 2013) or with scales assessing trait tendencies (Eley & Stevenson, 1999; Jardine et al., 1984; Kendler et al., 1987). For example, a genetic correlation of rg = 1.0 between trait anxiety and depression symptoms was observed in children and adolescents (Eley & Stevenson, 1999), and smaller but substantial genetic correlations were observed in adults (rg = .79–.88; Jardine et al., 1984).
Nonshared environmental influences were also correlated in both studies, but these associations were moderate (re = .43–.63), suggesting that environmental factors also differentially influence anxiety and depression. These nonshared environmental influences refer to factors that make twins within a family less correlated, in contrast with shared environmental influences that makes twins more similar. Thus, both can be correlated across traits. These results are consistent with evidence from clinical studies suggesting that generalized anxiety disorder and major depressive disorder have nearly identical genetic influences but only some overlapping nonshared environmental influences (Kendler & Myers, 2014; Kendler, Neale, Kessler, Heath, & Eaves, 1992; Lyons et al., 1998). Although many studies have examined how anxiety and depression symptoms are associated with cognition, few have considered the genetic and environmental overlap between cognition and anxiety/depression symptoms (Franz et al., 2011; Friedman et al., 2018; Routledge et al., 2017), especially EFs.
Unity and Diversity of Executive Functions
Existing research on EFs and anxiety or depression has focused on the three correlated EFs discussed by Miyake et al. (2000): inhibition, shifting, and WM updating. However, in the more recent unity and diversity model of EFs (Friedman & Miyake, 2017; Miyake & Friedman, 2012), variation in these well-known EF abilities can be reconceptualized into different latent factors: a general factor explaining variation across all EF tasks (common EF) and EF-specific factors that explain additional variance in certain processes not captured by common EF (e.g., shifting-specific, WM-specific). By reconceptualizing variance in inhibition, shifting, and WM updating, researchers aim to better specify the cognitive/biological underpinnings and correlates of the variance unifying different aspects of EFs.
Our recent work on the present sample at an earlier wave of assessment identified evidence for the common EF factor in middle age (Gustavson, Panizzon, Franz, et al., 2018). We also identified a WM-specific factor accounting for additional variation in WM span tasks not captured by common EF, but not a shifting-specific factor. As yet, there is no evidence for an inhibition-specific factor in adolescence, early adulthood, or middle age, suggesting that inhibition is entirely encapsulated by common EF (Friedman et al., 2016; Gustavson, Stallings, et al., 2017; Ito et al., 2015).
Importantly, genetic influences account for a large portion of variance in common EF and the specific factors. Heritability estimates of common EF range from 46% to 81% in young and middle adulthood (Friedman et al., 2016; Gustavson, Panizzon, Franz, et al., 2018) and are essentially 100% in adolescence (Engelhardt, Briley, Mann, Harden, & Tucker-Drob, 2015; Friedman et al., 2008). Furthermore, genetic influences tend to explain most of the phenotypic correlations between common EF and other cognitive and clinical constructs throughout the life span (Engelhardt et al., 2016; Gustavson, Panizzon, Franz, et al., 2018; Gustavson, Stallings, et al., 2017; Herd et al., 2014).
Executive Functions, Anxiety, Depression, and Aging
Theoretical reviews of attentional control theory have demonstrated that higher levels of subclinical anxiety symptoms are associated with poorer inhibition, shifting, and WM updating (Berggren & Derakshan, 2013; Eysenck & Derakshan, 2011; Eysenck et al., 2007). Associations between EFs and subclinical depression symptoms have received far less attention, and the findings are somewhat inconsistent. Studies of college students have yielded mixed results regarding associations between depression symptoms and EFs (Holmes & Pizzagalli, 2007; Wagner, Alloy, & Abramson, 2015). Several studies have also observed associations between EFs and trait anxiety above and beyond depression and/or showed that depression symptoms were not associated with their candidate EF task in a model with trait anxiety (Edwards, Edwards, & Lyvers, 2015; Gustavson, Altamirano, Johnson, Whisman, & Miyake, 2017; Gustavson & Miyake, 2016).
Existing studies have focused on these associations in adolescence and young adulthood. When older adults are included, the samples have had generally wide age ranges (Bunce, Handley, & Gaines, 2008; Routledge et al., 2017). However, at an earlier wave of assessment of the present sample (ages 51–60 years), levels of depression were weakly but significantly associated with performance on multiple EF tasks (rs = −.07 to –.12; Franz et al., 2011) and mostly explained by shared genetic influences. A more recent twin study of adults (ages 18–62 years) revealed that a combined measure of anxiety and depression symptoms (42-item Depression Anxiety Stress Scales) was associated with worse inhibition (β = −0.14) and shifting (β = −0.12; Routledge et al., 2017). Although nonsignificant, at least half of both associations were accounted for by genetic influences, with evidence for a significant environmental association for inhibition only.
Cognitive abilities begin to decline in late middle age, so this may be an important time to examine the associations among EFs, trait anxiety, and depression symptoms. For example, anxiety and depression are associated with increased risk for dementia (Byers & Yaffe, 2011; Gallacher et al., 2009; Petkus, Reynolds, et al., 2016). Additionally, EFs may be especially susceptible to age-related decline in middle age (de Frias, Dixon, & Strauss, 2006), and decline in EFs may underlie age-related decline in other cognitive abilities (Hasher & Zacks, 1988; Lustig, Hasher, & Zacks, 2007; Salthouse, Atkinson, & Berish, 2003). In fact, longitudinal analyses suggest that trait anxiety is associated with decline in multiple cognitive abilities, including processing speed and attention, with evidence for some bidirectionality in these associations (Petkus, Reynolds, Wetherell, Kremen, & Gatz, 2017). However, the Petkus et al. (2017) study did not have EF or WM measures other than digit span. It is surprising that associations between anxiety, depression, and EFs have received so little attention in the aging literature, given their relevance to cognitive aging. Furthermore, anxiety and depression have been associated with dorsolateral prefrontal cortex structure and function (Engels et al., 2010; Herrington et al., 2010), an area also associated with EFs (Banich, 2009), suggesting that these constructs are linked through neurobiological factors that likely include genetic influences (Stollstorff et al., 2013).
The Current Study
We tested the hypothesis that anxiety and depression symptoms are associated with worse common EF ability. We also hypothesized that the association between common EF and anxiety/depression would be accounted for primarily by shared genetic influences. Specifically, we expected that most of the genetic influences on trait anxiety would be shared with depression symptoms and that their genetic associations with common EF would be primarily explained by genetic influences common to anxiety and depression rather than those unique to either one. This contrasts with an alternative (but potentially complementary) hypothesis that there are some unique genetic or environmental influences that also account for these associations. If so, these findings could shed light on the genetic/environmental factors that differentiate anxiety from depression. We also examined whether the WM-specific factor would be associated with anxiety or depression, as these associations have yet to be explicitly studied. We tested these hypotheses in a twin study of community-dwelling middle-age men who were not selected for their levels of anxiety or depression.
Method
Participants
Analyses were based on 1,207 men from 563 full twin pairs (330 monozygotic pairs, 233 dizygotic pairs, 5 pairs with unknown zygosity, 71 twins unpaired) who participated in Wave 2 of the VETSA (age: M = 61.72 years, SD = 2.44). VETSA participants were recruited from the Vietnam Era Twin Registry from a previous study (Tsuang, Bar, Harley, & Lyons, 2001). They are similar to American men in their age group with respect to health and lifestyle factors (e.g., M years of education = 13.81, SD = 2.11, range = 5–20), and nearly 80% of the sample did not serve in combat or in Vietnam (Kremen, Franz, & Lyons, 2013; Kremen et al., 2006). Most of the sample identify as White, non-Hispanic (88.1%), with additional individuals identifying themselves as Black (5.1%), more than one race (2.7%), American Indian (< 1%), or Pacific Islander (< 1%). Hispanic individuals accounted for 3.0% of the sample.
Measures
The neuropsychological tasks comprising the EF battery are widely used and have been described in previous work on this sample at an earlier wave of assessment (Gustavson, Panizzon, Franz, et al., 2018; Kremen, Moore, Franz, Panizzon, & Lyons, 2014; Kremen et al., 2006). After creating each of the dependent measures described below, we adjusted each measure for age by creating residualized scores. Moreover, because 193 individuals completed the EF tasks for the first time in this second wave of assessment, we also adjusted scores on each task to account for practice effects for those individuals who had completed the tests twice (Elman et al., 2018). These practice effect adjustments were based on procedures outlined by Rönnlund, Nyberg, Bäckman, and Nilsson (2005), in which returning participants are compared with matched replacement participants and compared with dropouts to account for attrition effects.
Inhibition
Inhibition was assessed with the Stroop task (Golden & Freshwater, 2002; Stroop, 1935) and the AX–Continuous Performance Test (Braver et al., 2001; Kremen et al., 2011). For the Stroop, the dependent measure was a residualized score for the total number of words identified during the color-word condition after adjusting for performance on the two baseline conditions (word only and color only). For the AX–Continuous Performance Test, the dependent measure was an arcsine-transformed signal detection index (d′) based on the hit rate for AX trials minus the false-alarm rate for BX trials. Additionally, d′ values of less than 0 were set to 0 to reduce the tail of the distribution. On this test, the target stimulus is an “X” that is immediately preceded by an “A”. BX trials refer to trials on which an “X” is preceded by a non-“A” letter. Because 70% of trials are AX, BX trials require inhibition of a prepotent tendency to press the target button in response to an “X”.
Working memory span
WM span was assessed with the letter-number sequencing subtest of the Wechsler Memory Scale–III (Wechsler, 1997), the reading span task (Daneman & Carpenter, 1980), and the combined forward and backward digit span subtests of the Wechsler Memory Scale–III. The dependent measures for letter-number sequencing and digit span were the total number of trials passed. The dependent measure for the reading span was the total number of correct words recalled across all 15 trials.
Shifting
Shifting was assessed using the trail making test and the category switching trial of the verbal fluency test, both from the Delis-Kaplan Executive Function System (Delis, Kaplan, & Kramer, 2001). The dependent measure for the trail making test was a residualized score for the time taken to complete the switching trial after adjusting for the time on the single-task trials (number sequencing and letter sequencing). The dependent measure of category switching was a residualized score for category switching accuracy (the number of times a participant correctly switched categories) after adjusting for the number of correct responses on the category fluency trials (animals and boys’ names).
Anxiety symptoms
Levels of trait anxiety were measured with the trait form of the State-Trait Anxiety Inventory (STAI; Spielberger, 1983). The trait form of the STAI has 20 items that are responded to on a scale from 1 to 4 and asks “how you generally feel.” Reliability, computed using Cronbach’s alpha, was high (α = .93). Additionally, this scale was log-transformed to improve the distributional characteristics. For use in additional analyses, we assessed state anxiety using the log-transformed state form of the STAI, which also has 20 items responded to on a scale from 1 to 4 and assesses anxiety symptoms on the basis of “how you feel right now, that is, at this moment.” State anxiety was assessed on the day of testing; trait anxiety was assessed via a mailed survey approximately 4 weeks prior to testing.
Depression symptoms
Depression symptoms were measured with the Center for Epidemiologic Studies–Depression scale (CESD; Radloff, 1977), which was also assessed via a mailed survey approximately 4 weeks prior to testing. The scale has 20 items, each responded to on a scale from 0 to 3, and asks about the experience of symptoms “during the past week.” Internal consistency was high (α = .89), and this scale was also log-transformed to improve distributional characteristics.
General cognitive ability
For use in additional analyses, general cognitive ability was measured with the Armed Forces Qualification Test (Bayroff & Anderson, 1963), a 100-question multiple-choice paper-and-pencil test that is highly correlated with traditional measures of IQ and general cognitive ability (Lyons et al., 2009). Percentile scores were transformed on the basis of military norms to normalize the distribution.
Data analysis
Phenotypic analyses were conducted using Mplus Version 7.2 (Muthén & Muthén, 2010), and genetic/environmental analyses were conducted using the structural equation modeling package OpenMx in the R programming environment (Version 3.3.1; Boker et al., 2001; R Core Team, 2016). Both programs account for missing observations using full-information maximum likelihood. Mplus was preferred for phenotypic analyses to fit multiple regression analyses involving latent variables (i.e., controlling for general cognitive ability and state anxiety), and OpenMx was preferred for genetic analyses to obtain likelihood-based 95% confidence intervals (CIs).
Model fit was determined using chi-square tests, the root-mean-square error of approximation (RMSEA), and the comparative fit index (CFI). Good-fitting models had chi-square values less than 2 times the degrees of freedom, RMSEA values less than 0.06, and CFI values greater than .95 (Hu & Bentler, 1998). Significance of individual parameters was established with chi-square difference tests (by fixing those parameters to zero) or using likelihood-based 95% CIs. For phenotypic models, chi-square tests and standard errors were adjusted for the within-family clustering of the data using scaling factors provided in the output (i.e., by using the “type = complex” command in Mplus; Satorra & Bentler, 2001).
Genetic analyses were based on the classical assumptions in twin designs. Additive genetic influences (A) are assumed to correlate at 1.0 for monozygotic twins and at 0.5 for dizygotic twins because monozygotic twins share 100% of their alleles identical by descent and dizygotic twins share, on average, 50% of their alleles identical by descent. Shared environmental influences (C) are assumed to correlate at 1.0 for both types of twins. Nonshared environmental influences (E), which also include measurement error, are not correlated in either monozygotic or dizygotic twins. We also assume that means and variances are identical across twin pairs (Twin 1 vs. Twin 2) and across zygosity (monozygotic twins vs. dizygotic twins). These standard assumptions for univariate twin models extend to multivariate analyses and to analyses in which phenotypic correlations are decomposed into their genetic (rg), shared environmental (rc), and nonshared environmental (re) components.
Using an approach developed in Gustavson, Panizzon, Franz, et al. (2018), we first fit the confirmatory common pathway model of EFs alone. This model, displayed in Figure 1, parses variance in the seven neuropsychological EF tasks into common EF and WM-specific factors, and each of these factors is decomposed into its genetic, shared environmental, and nonshared environmental components. Next, we combined this model of EFs with the measures of trait anxiety and depression symptoms using a Cholesky decomposition, which is a common first step in twin analyses. The Cholesky is a simple decomposition of the genetic and environmental variances and covariances among common EF, WM-specific, trait anxiety, and depression symptoms and can be used to compute the genetic/environmental correlations between all constructs before examining more theoretically driven models. Finally, we fit a confirmatory factor model in which we directly modeled the associations between common EF and genetic/environmental influences, capturing shared and unique variance in anxiety and depression. This model is displayed in Figure 2 and is described in further detail below.

Model of the genetic (A), shared environmental (C), and nonshared environmental (E) influences on common EF- and WM-specific factors. Circles and ellipses indicate latent variables and squares indicate measured variables. The variance explained by each latent factor can be computed by squaring the factor loading on that factor. Significant factor loadings are displayed with black text and lines (p < .05). EF = executive function; WM = working memory; AX-CPT = AX–Continuous Performance Test.

Confirmatory model of the genetic and environmental associations between common EF, trait anxiety, and depression symptoms. Variance in trait anxiety and depression symptoms are decomposed into genetic/environmental influences common to both anxiety and depression (AA+D, EA+D) and genetic/environmental influences unique to anxiety (AAnx, EAnx) or depression (ADep, EDep). Genetic/environmental influences on working-memory-specific factor and shared environmental influences on all constructs are not pictured here, but all shared environmental influences were nonsignificant. Factor loadings and residual genetic/environmental influences on the EF tasks are also not shown here but were similar to those displayed in Figure 1. Significant factor loadings and correlations are displayed in bold with black lines (p < .05). EF = executive function; STAI = State-Trait Anxiety Inventory; CESD = Center for Epidemiologic Studies–Depression scale.
Results
Preliminary analyses
Descriptive statistics
Descriptive statistics for all measures are summarized in Table 1, shown before adjusting for age on each of the dependent measures. The full phenotypic correlation matrix between all measures is displayed in the Appendix. As would be expected in a middle-age community-dwelling sample, the large majority were at normal levels of anxiety and depression. For example, only 172 individuals (14.3% of the sample) were at risk for clinical depression (i.e., CESD score > 15).
Descriptive Statistics for All Measures in the Study
Note: Descriptive statistics for all dependent measures are shown before adjusting for age (for all measures) but after transformation (for measures of anxiety and depression). The final two columns display the cross-twin correlations for monozygotic (rMZ) and dizygotic (rDZ) twins. Boldface type indicates significant cross-twin correlations (p < .05). AX-CPT = AX–Continuous Performance Test; STAI = State-Trait Anxiety Inventory; CESD = Center for Epidemiologic Studies–Depression scale; AFQT = Armed Forces Qualification Test percentile score (transformed).
Unity and diversity model of EF
The genetic and environmental model of EFs alone, displayed in Figure 1, was largely consistent with the same genetic/environmental decomposition reported for these latent variables in the first wave of the VETSA (Gustavson, Panizzon, Franz, et al., 2018). Common EF accounted for variation and covariation in all seven neuropsychological EF tasks, and the WM-specific factor accounted for residual variation and covariation across the three WM span tasks not already captured by common EF. Common EF was explained primarily by genetic influences (heritability of a2 = .77, or the square of the factor loading of .88, 95% CI = [.37, .90]) but also by some nonshared environmental influences (e2 = .19, 95% CI = [.08, .29]) and nonsignificant shared environmental influences (c2 = .04, 95% CI = [.00, .39]). The WM-specific factor was also explained by nonshared environmental influences (e2 = .29, 95% CI = [.08, .50]). Although shared environmental influences accounted for almost all of the remaining variation in the WM-specific factor (c2 = .79), these influences were not significant because of the low power to detect the difference between genetic and shared environmental influences on this residual factor (Martin, Eaves, Kearsey, & Davies, 1978). Both the genetic and shared environmental influences on the WM-specific factor could not be removed from the model without a significant drop in fit, χ2diff(2) = 18.74, p < .001.
Phenotypic analyses
We examined the phenotypic associations between the two latent EF factors and the measures of anxiety and depression using a correlational measurement model, displayed in Table 2a, which fit the data well, χ2(21) = 35.20, p = .027, RMSEA = 0.024, CFI = .992. Confirming our hypothesis, lower common EF was significantly associated with greater trait anxiety (r = −.24) and depression symptoms (r = −.32), both χ2diff(1) > 14.05, ps < .001. In contrast, a low WM-specific factor was not significantly associated with trait anxiety (r = −.03) or depression symptoms (r = .03), both χ2diff(1) < .07, p > .785.
Latent Variable Correlations Between Executive Function and Anxiety/Depression Symptoms
Note: The final column displays the proportion of variation in each trait that is explained by genetic influences (a2), shared environmental influences (c2), or nonshared environmental influences (e2) in the Cholesky decomposition. Genetic/environmental influences differ slightly from those reported in the main text because they were derived from the Cholesky decomposition rather than the confirmatory model (see Fig. 2). Factor loadings on the EF latent variables are not shown but were similar to those displayed in Figure 1. Common EF and WM-specific factor were fixed to be uncorrelated. Boldface type indicates significant correlations and genetic/environmental variance components (p < .05). EF = executive function; WM = working memory.
Genetic and environmental models of anxiety, depression, and executive functions
The genetic/environmental correlations among EF factors, trait anxiety, and depression symptoms derived from the Cholesky decomposition are displayed in Tables 2b to 2d. This basic model fit the data well, χ2(313) = 400.14, p < .001, RMSEA = 0.021, CFI = .965.
To directly test the hypothesis that the associations between common EF and anxiety/depression were due to those genetic influences underlying the commonality between anxiety and depression, we fit a confirmatory model based on one outlined by Loehlin (1996). In this model—displayed in Figure 2—there are genetic and environmental factors that capture the common variance between anxiety and depression symptoms (AA+D, EA+D; with equated factor loadings) and that correlate with the genetic/environmental influences on common EF (AEF, EEF). Additionally, there are factors that capture the remaining genetic/environmental influences on anxiety (AAnx, EAnx) and depression symptoms (ADep, EDep), some of which are also correlated with common EF. 1 This model also fit the data well, χ2(321) = 405.81, p < .001, RMSEA = 0.020, CFI = .966, and did not fit worse than the Cholesky decomposition, χ2diff(8) = 5.67, p = .684.
Although not displayed in Figure 2, the WM-specific factor was included in this model but restricted to not correlate with anxiety or depression, consistent with the lack of phenotypic correlations. Additionally, because shared environmental influences accounted for little variance in any construct (see below), these influences were estimated (CA+D, CAnx, CDep) but not allowed to correlate with those for common EF (CEF) and are also not displayed in Figure 2. In both genetically informative models, genetic/environmental correlations were not estimated directly but derived algebraically from statistically equivalent covariance matrices.
Genetic influences
As shown in Figure 2, genetic influences explained a total of 36% of the variation in trait anxiety (a2 = .36, 95% CI = [.12, .48]) and 28% of the variation in depression symptoms (a2 = .28, 95% CI = [.04, .40]). All of the genetic influences on depression were explained by the common factor, and there were only small and nonsignificant genetic influences unique to anxiety (AAnx; a2 = .08, 95% CI = [.00, .12]). Thus, anxiety and depression have a near-complete genetic overlap.
Consistent with our hypothesis, there was a significant genetic correlation between common EF (AEF) and the genetic influences shared between anxiety and depression (AA+D; rg = −.44, 95% CI = [−1.0, −0.25]). This correlation accounted for 83% of the phenotypic correlation between common EF and trait anxiety and 69% of the phenotypic correlation between common EF and depression symptoms. In contrast, the genetic correlation between common EF and those genetic influences unique to anxiety (AAnx) was nonsignificant (rg = −.44, 95% CI = [−1.0, 1.0]), explaining only 3% of the phenotypic correlation.
Shared environmental influences
Although not pictured in Figure 2, shared environmental influences accounted for a small and nonsignificant amount of variation in both anxiety (c2 = .04, 95% CI = [.00, .25]) and depression (c2 = .04, 95% CI = [.00, .25]) in the confirmatory model. Although nonsignificant, these influences were explained entirely by the common factor for both anxiety and depression (CA+D), with no influences estimated on the specific shared environmental factors (CAnx, CDep).
As described above, we did not estimate shared environmental correlations between common EF and anxiety/depression symptoms. If a shared environmental correlation between common EF (CEF) and the influences shared between anxiety and depression (CA+D) were added to this model, it would be nonsignificant (rg = −1.0, 95% CI = [−1.0, 1.0]) and would explain less than 1% of the phenotypic correlation. In this alternate model including a shared environmental correlation, the significant genetic correlation between common EF and genetic influences shared between anxiety and depression (AA+D) would remain similar in magnitude but would be nonsignificant (rg = −.38, 95% CI = [−1.0, 0.06]), most likely because of a loss of power (Martin et al., 1978).
Nonshared environmental influences
As shown in Figure 2, nonshared environmental influences explained most of the variation in trait anxiety (e2 = .59, 95% CI = [.52, .68]) and depression symptoms (e2 = .68, 95% CI = [.60, .76]). Most of these influences were explained by the common factor (EA+D), explaining 38% of the total variance in either trait, but there were also influences unique to anxiety or depression (explaining the remaining 19% and 26% of variance, respectively).
The nonshared environmental correlation between common EF (EEF) and the influences shared between anxiety and depression (EA+D) was nonsignificant (re = −.13, 95% CI = [−.35, .09]), explaining 14% and 12% of the phenotypic correlation between common EF and anxiety or depression, respectively. Interestingly, there was a significant nonshared environmental correlation between common EF and environmental influences specific to depression (EDep; re = −.27, 95% CI = [−.52, −.04]), explaining the remaining 20% of their phenotypic correlation.
Additional analyses
We conducted two additional phenotypic multiple regression analyses to address potential confounds. First, there is a considerable phenotypic correlation between common EF and measures of general cognitive ability or intelligence (Friedman et al., 2008; Gustavson, Panizzon, Franz, et al., 2018), as evidenced by the correlation of r = .71 between common EF and general cognitive ability in these data. Therefore, it was possible that some of the associations between EF and anxiety or depression symptoms described here were explained by a broader association between general cognitive ability and anxiety and/or depression symptoms rather than common EF. After controlling for general cognitive ability using multiple regression procedures, the phenotypic association between common EF and trait anxiety remained significant (β = −0.22, p = .016), as did the association between common EF and depression symptoms (β = −0.33, p = .001). Moreover, although general cognitive ability was negatively correlated with both trait anxiety (r = −.21, p < .001) and depression symptoms (r = −.25, p < .001), these associations were no longer significant after controlling for the EF latent factors (β = −0.05, p = .499 for anxiety, β = −0.01, p = .910 for depression). Thus, if anything, the association between general cognitive ability and anxiety/depression seems to be explained by variance in common EF, rather than vice versa.
Second, although we focused here on the association between EFs and trait anxiety, it is possible that these associations are actually explained by state anxiety. Again using multiple regression, the phenotypic association between common EF and trait anxiety remained significant after controlling for state anxiety (β = −0.21, p = .001). State anxiety was correlated with both common EF (r = −.16, p = .002) and the WM-specific factor (r = −.14, p = .022), but it was not associated with common EF in the multiple regression model (β = −0.05, p = .366). Additional genetic analyses involving state anxiety, displayed in the Appendix (Fig. A1), revealed that state anxiety was not significantly explained by genetic influences common to trait anxiety and depression. Rather, it was primarily explained by nonshared environmental influences distinct from those on trait anxiety or depression (e2 = .64). There were no significant genetic or environmental correlations between common EF and state anxiety, all χ2diff(1) < 3.51, p > .061. Therefore, the association between EFs and trait anxiety described above does not appear to be accounted for by people’s current level of anxiety.
Discussion
We examined the associations between EFs and anxiety and depression symptoms in a community-dwelling sample of middle-age male twins. First, we confirmed that the associations between EFs and anxiety/depression symptoms were due to common variance spanning multiple EFs (i.e., common EF). Second, we showed that these associations were mostly explained by genetic influences. Specifically, the correlation between genetic influences on common EF (AEF) and genetic influences shared between anxiety and depression (AA+D) accounted for most of the phenotypic correlation between common EF and anxiety (83%) or depression (69%). Additionally, there was a significant correlation between the nonshared environmental influences on common EF and those environmental influences unique to depression (EDep), accounting for 20% of their phenotypic correlation. Therefore, anxiety and depression were not associated with EFs entirely through the same mechanisms; depression was associated with common EF through both genetic and environmental influences, whereas anxiety was associated primarily through genetic influences.
These findings have some direct implications for theoretical models, such as attentional control theory, that aim to explain the associations between EFs and subclinical anxiety (Eysenck & Derakshan, 2011; Eysenck et al., 2007). Three of attentional control theory’s hypotheses state that levels of anxiety are associated with inhibition, shifting, and WM updating using the framework inspired by Miyake et al. (2000), but these hypotheses have yet to be applied to the more recent unity and diversity model of EFs (Miyake & Friedman, 2012). Our results suggest that it may be more parsimonious to view this as an association between anxiety and the general processes underlying EFs across multiple subdomains (i.e., common EF).
In contrast, the WM-specific factor was uncorrelated with either anxiety or depression symptoms. Meta-analytic estimates have shown that there are significant associations between trait anxiety and WM span (Moran, 2016). Our results may appear to be in conflict with the meta-analytic results, but it should be noted that the WM-specific factor captures only the variance in WM span that is independent of common EF. The remaining variance in WM is captured by common EF, which is associated with trait anxiety. Models that directly quantify these sources of variance may be useful in further understanding how different aspects of WM processes are associated with other constructs. Although there was no shifting-specific factor in this sample, future work should investigate the extent to which anxiety is related to common EF versus shifting-specific processes as well, as there is some debate about whether associations between anxiety and shifting are specific to certain situations (Edwards et al., 2015; Gustavson, Altamirano, et al., 2017).
Previous work has suggested that almost all of the genetic influences on trait anxiety are shared with those on depression symptoms and that nonshared environmental influences are moderately correlated but are also unique to each construct (Eley & Stevenson, 1999; Jardine et al., 1984; Kendler et al., 1987). Similar results were observed here. All of the genetic variance in depression was accounted for by the shared anxiety/depression genetic factor (AA+D), and only a small and nonsignificant portion of the variance in anxiety was accounted for by unique genetic influences (AAnx). Importantly, some of these shared genetic influences between anxiety and depression were the same as those genetic influences on common EF, and there was no evidence that the unique influences on anxiety were associated with common EF. Although genetic influences accounted for much of the covariance between common EF and anxiety/depression, it should be noted that most of the genetic variance in common EF is independent of anxiety or depression (i.e., only about 20% of the total genetic variance in common EF was shared with genetic variance on anxiety and/or depression).
Interestingly, the environmental influences unique to depression symptoms (EDep) were negatively associated with common EF (re = −.27), explaining 20% of their phenotypic correlation. It is unclear what environmental influences explain this association; it is clear only that they are distinct from the environmental influences on trait anxiety and are influences that differentiate twins (as opposed to shared environmental influences that make twins more similar). Moreover, we could have tested an alternate but statistically equivalent model in which common EF was negatively correlated with environmental influences common to anxiety/depression (ECom) but positively correlated with influences unique to anxiety (EAnx). We found it more likely that there are environmental influences unique to common EF and depression and not anxiety, rather than environmental influences that act in opposite directions on common EF and anxiety (completely cancelling each other out). Nevertheless, it will be important to examine these associations further as they may shed light on the situations that distinguish depression from anxiety. Given these results, it is also not surprising that despite their almost complete genetic overlap, anxiety and depression can have some independent associations with cognition and/or dementia (e.g., controlling for one another; Burton, Campbell, Jordan, Strauss, & Mallen, 2013; Petkus et al., 2017).
In addition to the novel genetic findings, our findings extend work on the association between EFs and anxiety and depressive disorders to show that these associations exist even in the normal range of anxiety and depression symptoms. These results are consistent with associations between EFs and generalized anxiety disorder or major depressive disorder (Caspi et al., 2014; Snyder, 2013; Snyder, Miyake, & Hankin, 2015) and evidence that much of these associations are explained by the shared variance between anxiety and depression (Caspi et al., 2014; Martel et al., 2017). Nevertheless, future work on clinical-level associations with EFs may benefit from similar multivariate models that emphasize common EF (e.g., Martel et al., 2017).
These results are also relevant to recent efforts to use cognitive control training as an intervention for anxiety and/or depression. There is mounting evidence that interventions that target EF processes (e.g., attentional control, WM) are successful at reducing anxiety and depression symptoms in individuals with elevated symptoms (Bettis et al., 2017; Calkins, McMorran, Siegle, & Otto, 2015; Owens, Koster, & Derakshan, 2013; Sari, Koster, Pourtois, & Derakshan, 2016). Our findings that the associations between EFs and anxiety/depression symptoms are driven by common EF suggest that interventions may be most successful if they incorporate EF tasks from multiple domains (inhibition, shifting, WM), rather than modular approaches that target only one process at a time (e.g., WM span only). Moreover, given the substantial overlap between anxiety and depression, EF training for anxious individuals may improve their depression symptoms (and vice versa), but existing studies have typically considered only one construct or the other.
Strengths and limitations
Examining associations at the level of latent variables was important for being able to draw conclusions. By isolating variance in common EF and WM-specific latent variables, we controlled for task-specific variance on each EF task, including measurement error. As displayed in the Appendix, there were only weak correlations between individual EF tasks and trait anxiety or depressive measures (rs = −.01 to −.19), consistent with previous estimates (Franz et al., 2011; Routledge et al., 2017). Phenotypic correlations between common EF and trait anxiety (r = −.26) or depression symptoms (r = −.33) were more substantial at the latent variable level, suggesting that the true associations may have been underestimated in prior studies.
These results also add to our growing understanding of the genetic influences on common EF. Previous work using twin studies in adolescence and young adulthood has shown that the associations between common EF and other clinically relevant constructs tend to be explained by genetic influences more so than environmental influences (Gustavson, Stallings, et al., 2017; Young et al., 2009). This was not as surprising in adolescence when nearly 100% of the individual differences are explained by genetic influences (Engelhardt et al., 2015; Friedman et al., 2016), but by midlife, environmental influences also explain a considerable amount of the variance (Gustavson, Panizzon, Franz, et al., 2018). Nevertheless, genetic influences accounted for most of the phenotypic association between common EF and trait anxiety and much of the phenotypic correlation between common EF and depression symptoms, suggesting that underlying genetic risk factors continue to account for its association with other constructs even in middle age.
Our results are cross-sectional and therefore cannot directly inform the causal associations between these constructs. Poor EFs may predispose individuals to experience uncontrollable negative thoughts (e.g., worry and rumination) and lead to anxiety or depression, but it is also possible that prolonged anxiety or depression episodes result in changes to the brain that affect cognitive control abilities involved in EFs. A longitudinal analysis has provided evidence for bidirectionality of associations between anxiety and cognition, including processing speed, attention, and memory, suggesting that these associations may reinforce each other over time (Petkus et al., 2017), although EFs were not assessed in that study. Furthermore, the genetic variance in anxiety may increase in late life (past age 70 years), suggesting that these findings in middle age may differ with increasing age (Petkus, Gatz, Reynolds, Kremen, & Wetherell, 2016). To the extent that the associations among common EF, anxiety, and depression are explained by genetic influences, shared genetic risk factors may be an underlying cause of both cognition and psychopathology.
Finally, although we examined individual differences in EFs at the level of latent variables, our assessment of anxiety and depression was more limited. The trait scale of the STAI does not assess anxiety symptoms as comprehensively as scales such as the 42-item Depression Anxiety Stress Scales or Index of Depression and Anxiety Symptoms (Lovibond & Lovibond, 1995; Watson et al., 2007), and it has been criticized for assessing depression in addition to anxiety (Bieling, Antony, & Swinson, 1998; Caci et al., 2003). The CESD was designed to measure depression symptoms only in the past week. However, the CESD scores were correlated with CESD scores assessed in Wave 1 of VETSA approximately 6 years earlier (r = .58, N = 1,000; Franz et al., 2011), suggesting that they captured some relatively stable tendencies to experience depression symptoms. The genetic association between common EF and depression symptoms also suggests that this association is explained in part by stable influences, especially considering the perfectly stable genetic influences on EFs over the past 6 years in this sample (rg = 1.0; Gustavson, Panizzon, Elman, et al., 2018). Nevertheless, it will be useful to further examine these associations with multiple measures of anxiety and depression symptoms, including those that assess tendencies for depression on longer time scales, to examine how strongly these associations are driven by current symptoms versus stable personality traits and to model their shared and unique variance at the level of latent variables.
Conclusion
We found that EFs are associated with anxiety and depression symptoms through general sources of variance in common EF and the genetic variance underlying anxiety and depression. These associations were observed in late middle age, a period during which many cognitive abilities begin to decline. Given that both anxiety and depression symptoms are relevant to age-related cognitive decline and dementia (Byers & Yaffe, 2011; Gallacher et al., 2009; Petkus et al., 2017; Petkus, Reynolds, et al., 2016) and that the decline in EFs may underlie age-related decline in other cognitive abilities (Lustig et al., 2007; Persad, Abeles, Zacks, & Denburg, 2002), these results suggest that it is important to consider EFs when examining the dynamic association between anxiety, depression, and cognitive decline.
Footnotes
Appendix
Phenotypic Correlations Between All Measures of the Study
| Measure | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
|---|---|---|---|---|---|---|---|---|---|---|
| 1. Stroop | ||||||||||
| 2. AX-CPT |
|
|||||||||
| 3. Trail making test |
|
|
||||||||
| 4. Category switching |
|
|
|
|||||||
| 5. Letter-number |
|
|
|
|
||||||
| 6. Reading span |
|
|
|
|
|
|||||
| 7. Digit span |
|
|
|
|
|
|
||||
| 8. STAI–trait | −0.05 |
|
|
−0.01 |
|
|
|
|||
| 9. CESD |
|
|
|
−0.02 |
|
|
|
|
||
| 10. STAI–state | −0.02 | −0.06 |
|
−0.05 |
|
|
|
|
|
|
| 11. AFQT |
|
|
|
|
|
|
|
|
|
|
Note: Boldface type indicates significant correlations (p < .05). AX-CPT = AX–Continuous Performance Test; STAI = State-Trait Anxiety Inventory; CESD = Center for Epidemiologic Studies–Depression scale; AFQT = Armed Forces Qualification Test.
Acknowledgements
The content of this article is the responsibility of the authors and does not represent official views of the National Institute on Aging, National Institutes of Health, or the Veterans Administration. Numerous organizations provided invaluable assistance in the conduct of the Vietnam Era Twin Registry, including U.S. Department of Veterans Affairs; U.S. Department of Defense; National Personnel Records Center, National Archives and Records Administration; Internal Revenue Service; National Opinion Research Center; National Research Council, National Academy of Sciences; and Institute for Survey Research, Temple University. The authors gratefully acknowledge the continued cooperation of the twins and the efforts of many staff members.
Action Editor
Kelly L. Klump served as action editor for this article.
Author Contributions
D. E. Gustavson developed the study concept and analyzed the data under the supervision of W. S. Kremen. Testing and data collection were performed by staff of the Vietnam Era Longitudinal Twin Study of Aging (principal investigators were W. S. Kremen and M. J. Lyons). All the authors contributed to the study design, provided critical revisions, and approved the final manuscript for submission.
Declaration of Conflicting Interests
The author(s) declared that there were no conflicts of interest with respect to the authorship or the publication of this article.
Funding
This research was supported by National Institute on Aging Grants AG050595, AG022381, and AG047903. The Cooperative Studies Program of the U.S. Department of Veterans Affairs provided financial support for development and maintenance of the Vietnam Era Twin Registry.
