Abstract
This study examines additive, curvilinear, and interactive relations of anxiety and depression with several subjective indicators of intrapersonal (i.e., hope, self-compassion, shame) and interpersonal (i.e., social connectedness, quality of social relationships) functioning in a sample of adults (N = 547, Mage = 43.37 ± 12.02, female = 56.88%) seeking treatment for psychological difficulties. Results of complementary analyses were largely consistent with the hypothesis that increasing levels of anxiety and depression would correspond with worse psychosocial functioning, although nonlinear relations indicated that the effect of depression progressively attenuated at higher levels of symptom severity. Whereas the findings generally supported additive effects of anxiety and depression, the hypothesis that there would be synergistic effects of anxiety and depression was not supported. Supplementary group comparisons revealed that the functional implications of subsyndromal combinations of anxiety and depression may be comparable to those associated with symptoms that meet more traditional standards (i.e., syndromal or dimensional definitions) of comorbid anxiety–depression. The findings offer further insight into the complex relations of anxiety and depression with psychosocial functioning and emphasize the importance of detecting and offering appropriate treatments for anxiety and depression symptoms that coexist at subsyndromal levels.
Introduction
Mental illness is one of the largest contributors to the global burden of disease and accounts for 32.4% of disability worldwide (Vigo et al., 2016). Of the mental illness categories identified among the top 20 contributors to the global burden of disease (see Vos et al., 2015), those ranked highest are major depression (second) and anxiety disorders (ninth). Global prevalence estimates (e.g., Baxter et al., 2014) of major depressive disorder (4.4%) and anxiety disorders (4%) underscore the magnitude of the burden that each of these categories of mental illness carries. Yet, anxiety and depression are highly comorbid (Kroenke et al., 2016). A recent cross-national epidemiological analysis involving 74,045 adults from 24 countries estimated the 12-month prevalence of major depressive disorder comorbid with at least one 12-month anxiety disorder to be 41.6% (Kessler et al., 2015). Simultaneous co-occurrence may reflect a core transdiagnostic mechanism that underlies the cognitive vulnerabilities of anxiety and depression (Hong & Cheung, 2015). Psychometric analysis has found support for a bifactor structure of anxiety and depression, which suggests an underlying dimension may be common to both psychopathologies (Kroenke et al., 2016). An alternative view is that anxiety represents a vulnerability phenotype that increases the risk of stress-induced depression (Weger & Sandi, 2018). This proposition is consistent with epidemiological evidence that has found a majority of individuals with anxious major depressive disorder tend to report earlier onset of anxiety disorders than major depressive disorder (Kessler et al., 2015).
Alongside progress that has been made toward understanding the causal determinants of comorbid anxiety–depression, researchers have sought to determine whether the debilitating effects of either disorder occurring in isolation differs from co-occurring anxiety and depression. In a review of neurobiological evidence from 24 studies, Ionescu, Niciu, Mathews, et al. (2013) assert that anxiety and depression comorbidity is neurobiologically different from depression without anxiety and is associated with distinct psychosocial outcomes. Evidence from a range of studies (e.g., Fusar-Poli et al., 2014; Saris et al., 2017; Zhou et al., 2017) suggests that relations of comorbid anxiety and depression with diminished psychological and interpersonal functioning are amplified when compared to the effects of anxiety or depression alone.
Although anxiety and depression comorbidity has generally been linked to less favorable psychosocial functioning, there are several reasons for further examining anxiety and depression in combination. First, existing research in this area has often relied on modeling approaches that emphasize linear additive effects of anxiety and depression (e.g., Brown et al., 2010), thereby overlooking the salience of nonlinear and multiplicative effects that could provide useful insight into the complex interrelations of anxiety and depression. Second, researchers have typically applied syndromal (i.e., concurrent formal diagnosis of major depressive disorder and one or more anxiety disorder) or dimensional (i.e., diagnosis of major depressive disorder based on formal diagnostic criteria, along with concurrent anxiety symptoms that meet a given threshold on psychometrically validated measures) definitions to classify comorbid anxiety and depression (Ionescu, Niciu, Henter, et al., 2013). While a finite approach to classifying comorbidity is advantageous because it aligns more closely with formal diagnostic criteria that may be of clinical significance, overlapping symptoms of anxiety and depression frequently coexist at subsyndromal levels and can have similar effects on psychosocial functioning compared to symptoms that meet diagnosable levels (Wanders et al., 2016). Therefore, consideration of subsyndromal levels of anxiety and depression may be important for early detection of individuals at risk of symptom progression, formulating prognoses, and effectively treating specific configurations of anxiety and depression. In the current study, a multi-analytic approach is used to examine relations of anxiety and depression with several subjective markers capturing overarching intrapersonal and interpersonal domains of psychosocial functioning.
Anxiety, depression, and subjective indicators of psychosocial functioning
At the intrapersonal level, criterion variables of interest included hope, self-compassion, and shame. The concept of hope refers to a future-oriented thought process involving a person’s self-perceived capacity to find practical means (i.e., pathways) of reaching their goals and the motivation (i.e., agency) to attain them (Snyder, 2005). According to hope theory (for an overview, see Snyder et al., 2002), the pathways and agency components represent distinct processes that are both necessary ingredients for successfully yielding hopeful thinking. An abundance of research (for reviews, see Gallagher & Lopez, 2018) has linked hope with lower levels of internalizing symptoms (e.g., anxiety, depression), although evidence suggests that hope and psychological well-being have a dynamic influence on each other (Gum, 2018). For example, symptoms consistent with anxiety (e.g., persistent worry) or depression (e.g., apathy) may undermine goal-oriented thought processes (Kashdan et al., 2006; Thimm et al., 2013) that could promote adaptive goal pursuits and lead to improved mental health (Snyder et al., 2002).
Self-compassion is a positive self-orientation that involves being kind, loving, and understanding toward the self during instances of difficulty, inadequacy, and disappointment (Neff, 2003a). Conceptualizations of self-compassion assert that it buffers against the negative effects of self-denigrating cognitive-emotional experiences (e.g., self-criticism) that ordinarily accompany failure, thereby supporting psychological well-being (Neff, 2003b). Self-compassion may also be facilitated when lower levels of psychopathology (e.g., anxiety, depression) are experienced, with more severe psychopathology a potential barrier to taking a compassionate stance toward the self (MacBeth & Gumley, 2012). In support of theoretical links between self-compassion and mental health, findings of various studies indicate that self-compassion is inversely associated with symptoms of anxiety and depression (for a review, see Barnard & Curry, 2011).
Compared to the supportive role that hope and self-compassion play in sustaining well-being, shame is often associated with impaired psychological functioning. Shame is a dysphoric state involving perceived inferiority about the global self (Leach, 2017). Particularly when shame reflects a view that the whole self is unalterably flawed, it can have a profound impact on self-worth and debilitating mental health consequences (for an overview, see Cibich et al., 2016). Although experienced shame may heighten vulnerability to anxiety and depression, the onset of such internalizing symptoms might also be interpreted as evidence of a flawed self and therefore precipitate feelings of shame (for meta-analytic reviews, see Cândea & Szentagotai-Tăta, 2018; Kim et al., 2011).
Interpersonal functioning was assessed via two criterion variables, namely, social connectedness and perceived quality of social relationships. Social connectedness is an attribute reflecting one’s sense of self in relation to the social world (Lee & Robbins, 2000). In contrast to individuals with higher social connectedness, those on the lower end of the continuum have difficulty identifying with and feeling close to others, experience discomfort in social settings, and feel misunderstood or isolated (Lee et al., 2001). As people are inherently social beings, these kinds of relational challenges may evoke distress and contribute to psychological maladjustment, including depression and anxiety (see Lee & Robbins, 1998; Williams & Galliher, 2006). Similarly, mental health is also affected by the quality of a person’s social relationships. Findings of numerous studies link depression and anxiety with impairments on a range of behavioral (e.g., social support) and affective (e.g., belongingness) indicators of social functioning (e.g., Cacioppo et al., 2010; Saris et al., 2017), the effects of which can further undermine the quality of relationships a person is able to initiate and maintain (see Trompenaars et al., 2007).
The present study
Extending previous research, this study uses complementary analytical approaches to examine the independent and joint effects of anxiety and depression on indicators of intrapersonal and interpersonal functioning. It was expected that (a) outcomes resonating with worse intrapersonal (i.e., lower hope and self-compassion, higher shame) and interpersonal (i.e., lower social connectedness and quality of social relationships) functioning would be found among individuals with more severe symptoms of anxiety and depression, and (b) anxiety and depression would be synergistically associated with psychosocial functioning, such that relations of depression with lower quality psychosocial functioning would be amplified at higher levels of anxiety.
Methods
Participants
The sample (N = 547) comprised of individuals who had elected to attend a short-term residential psychotherapeutic treatment program (see Table 1, for sample characteristic details). Participants ranged between 19 and 73 years of age (Mage = 43.37, standard deviationage = 12.02), the majority of which were females (56.88%), self-identified as White (92.31%), and were married or in a domestic partnership (51.55%). Most had completed a bachelor’s degree or higher (72.76%), were currently employed (78.90%), and affiliated religiously with Christianity (63.79%).
Sample characteristics (N = 547).
Note. SD: standard deviation.
Measures
Anxiety
The generalized anxiety disorder seven-item (GAD-7; Spitzer et al., 2006) scale was used to measure the severity of generalized anxiety symptomology experienced by participants during the past two weeks (e.g., “Feeling nervous, anxious or on edge”). A four-point scale (0 = not at all; 3 = nearly every day) is used to rate each item. A recent meta-analysis (Plummer et al., 2015) found that sensitivity (.74) and specificity (.83) for detecting generalized anxiety disorder at a cut-off score of ≥10 were comparable to those reported in the initial validation study. A cut-off score of ≥10 on the GAD-7 has also moderate sensitivity (.66 to .74) and specificity (.80 to .81) for detecting other anxiety-related disorders, including posttraumatic stress disorder, panic disorder, and social anxiety disorder (Kroenke et al., 2007). Categories typically used to delineate severity of anxiety on the GAD-7 (see Spitzer et al., 2006) include minimal (≤4), mild (5–9), moderate (10–14), and severe (≥15). Internal consistency of the GAD-7 in the current sample was ωt = .92.
Depression
Depression symptoms were measured using the Patient Health Questionnaire-9 (PHQ-9; Kroenke et al., 2001). The nine-item scale assesses the frequency of depressive symptomology experienced during the past two weeks (e.g., “Feeling down, depressed, or hopeless”). Participants use a four-point scale (0 = not at all; 3 = nearly every day) to respond to the items. Several meta-analyses have revealed an optimal sensitivity and specificity criterion of the PHQ-9 for detecting major depressive disorder is ≥10 (see Levis et al., 2019; Mitchell et al., 2016). Commonly used categories to classify depression severity on the PHQ-9 (see Kroenke et al., 2001) include minimal (≤4), mild (5–9), moderate (10–14), moderately severe (15–19), and severe (≥20). In this sample, internal consistency of the PHQ-9 was ωt = .90.
Hope
Participants completed the Adult State Hope Scale (ASHS; Snyder et al., 1996). The ASHS contains six items designed to measure goal-directed thinking (e.g., “I can think of many ways to reach my current goals”), which are evenly distributed across two subscales (i.e., agency and pathways). Items are rated on an eight-point scale (1 = definitely false; 8 = definitely true), with subscale scores added for a total hope score (Snyder et al., 1996). The convergent, concurrent, and discriminant validity of the ASHS has been supported in several studies (Curry et al., 1997; Martin-Krumm et al., 2015; Snyder et al., 1996). Estimated internal consistency of total ASHS scores reported in prior studies was ≥.71 (Martin-Krumm et al., 2015; Snyder et al., 1996). In this study, the internal consistency of the total ASHS was ωt = .90.
Self-compassion
The Self-Compassion Scale (SCS; Neff, 2003b) measured the tendency of participants to be caring, kind, and understanding toward themselves during moments of difficulty, inadequacy, or failure. The 26-item scale consists of three sets of subscales with contrasting poles: self-kindness (five items) and self-judgment (five items), common humanity (four items) and isolation (four items), and mindfulness (four items) and overidentified (four items). Participants use a five-point scale (1 = almost never; 5 = almost always) to rate each item (e.g., “I’m kind to myself when I’m experiencing suffering”). Evidence from clinical and nonclinical samples of various countries supports the factorial, convergent, and criterion validity of the SCS (Neff, 2003b; Neff et al., 2017, 2019). Internal consistency values of total SCS scores range from .91 to .95 (Neff, 2003b; Neff et al., 2017), and test–retest reliability estimates over various intervals (e.g., three weeks, four weeks) were ≥.78 (Castilho et al., 2015; Neff, 2003b). A total self-compassion score was used in this study, which was generated by reverse scoring the negative subscale items, calculating subscale means, and adding together the means of all six subscales (ωt = .95).
Shame
Shame was measured using the 25-item Experience of Shame Scale (ESS; Andrews et al., 2002). The ESS captures various aspects of characterological (12 items), behavioral (nine items), and bodily shame (four items). Using a four-point scale (1 = not at all; 4 = very much), participants rated the items (e.g., “Have you felt ashamed of the sort of person you are?”) with reference to the past month. Findings of several studies support the concurrent, convergent, and incremental predictive validity of the ESS (Andrews et al., 2002; Vizin et al., 2016). Estimated internal consistency and test–retest reliability of total ESS scores over an 11-week interval were .92 and .83, respectively (Andrews et al., 2002; Cheung et al., 2004). In this study, total scores were derived by summing responses to all 25 items (ωt = .96).
Social connectedness
The Social Connectedness Scale—Revised (SCS-R; Lee et al., 2001) was used to assess participants’ sense of interpersonal closeness. The measure consists of 20 items (e.g., “I feel close to people”) that are rated on a six-point scale (1 = strongly disagree; 6 = strongly agree). After reverse scoring the 10 negatively phrased items, responses are summed for a total score. Research supports the convergent and discriminant validity of the SCS-R (Lee et al., 2001, 2008). Prior studies have reported internal consistency of scores for the SCS-R that range from .92 to .94 (Grieve et al., 2013; Lee et al., 2001, 2008). In the current sample, estimated internal consistency of the SCS-R was ωt = .95.
Social relationship quality
Participants completed the three items comprising the social relationships domain of the World Health Organization Quality of Life—BREF scale (WHOQOL-BREF; The WHOQOL Group, 1998). The items (e.g., “How satisfied are you with your personal relationships?”) measure perceived quality of interpersonal relationships, each of which is rated on a five-point scale (1 = not at all; 5 = completely). Scale scores are calculated by multiplying the mean of the three items by four. The social relationships domain of the WHOQOL-BREF has demonstrated appropriate convergent and discriminant validity in samples from a variety of countries (Hawthorne et al., 2006; Skevington & McCrate, 2012). The estimated internal consistency of scale scores reported in previous research ranges from .66 to .68 (Skevington et al., 2004; The WHOQOL Group, 1998). Internal consistency of the social relationships domain in this study was ωt = .70.
Physical health quality
Drawing on evidence linking internalizing symptoms and physical health problems (see Scott et al., 2016), the physical health domain of the WHOQOL-BREF scale (The WHOQOL Group, 1998) was used to measure perceived quality of physical health for inclusion as a covariate in the primary analyses. A five-point scale (1 = not at all; 5 = completely) is used to rate each of the seven items (e.g., “How well are you able to get around?”), two of which are reverse coded. Scale scores are derived by multiplying the mean of the seven items by four. Evidence supports the convergent and discriminant validity of the physical health domain (Hawthorne et al., 2006; Skevington & McCrate, 2012). Internal consistency values for the scale were ≥.82 across a number of studies (Skevington et al., 2004; The WHOQOL Group, 1998). In this study, internal consistency for this index was ωt = .80.
Procedure
This study utilizes baseline data from an ongoing longitudinal intervention project assessing the efficacy of a short-term residential multimodal experiential psychotherapy program. Data were collected between January 2018 and February 2019. Ethical approval to conduct this project was granted by an independent review board (protocol 2017/11/23). Upon arrival at the program premises, program attendees were invited to participate in the study. Interested individuals were directed to a large group room where they were presented with information about the nature of the project and their participation in it. Those who agreed to participate provided written informed consent and then self-completed the battery of measures. Participants first completed a set of sociodemographic items followed by a randomized ordering of the measures. 1
Statistical approach
Data preparation and exploratory analyses
The data were initially screened for missing values. A negligible quantity (0.13%) of item-level data was missing. They were imputed using a nonparametric random forest technique (Stekhoven & Bühlmann, 2012) based on 10,000 iterations (proportion falsely classified = .22). Standardized values of the primary study variables were screened for univariate outliers (z > |3.29|, p < .001), three of which were found (z ≥ |3.51|). Mahalanobis distance, χ2(8) = 26.12, p < .001, detected three multivariate outliers (D2 ≥ 29.34). An a priori decision was taken to omit all outliers before proceeding with subsequent analyses. Omega total (ωt) was used to estimate the internal consistency reliability of all measures (Revelle & Zinbarg, 2009), a procedure that is favored over alternative approaches (e.g., alpha) because it makes fewer assumptions about the internal structure of a scale (Peters, 2014). Internal consistency values for all measures were ≥.70. Univariate skewness (maximum = |0.59|) and kurtosis (maximum = |0.88|) did not reveal any preliminary concerns with normality. Pearson correlations were used to explore bivariate relations among the primary study variables.
Polynomial regression of anxiety and depression
To examine additive and synergistic effects of anxiety and depression, separate polynomial regression models were performed by regressing each of the outcome variables (Z) included in this study (i.e., hope, self-compassion, shame, social connectedness, and quality of social relationships) on anxiety (X), depression (Y), and the interaction term between anxiety and depression (X × Y). Regression models were estimated with quadratic terms (i.e., squared value) of anxiety (X2) and depression (Y2) included, as an interaction effect involving highly correlated variables (as is the case with anxiety and depression) may mask an underlying nonlinear effect (Matuschek & Kliegl, 2018) and lead to spurious conclusions about synergistic effects (i.e., an interaction effect may be observed when there is no true interaction). Anxiety and depression were mean centered prior to model estimation (Schielzeth, 2010). Age, sex, race/ethnicity, and quality of physical health were included as covariates in each model. When evidence of a quadratic effect emerged, the two-line approach (see Simonsohn, 2018) was used to test whether the quadratic function conforms to a U-shape association with an outcome. The approach uses an algorithm to identify an optimal breakpoint along the predictor and interrupted regressions to estimate average slopes, one for low values (i.e., below the breakpoint) and another for high values (i.e., above the breakpoint) of the predictor. The procedure provides an indication of whether the signs of the average slopes derived for low and high values of a predictor vary, thereby benefiting the interpretation of nonlinear effects. For each regression model, visual inspection of quantile-quantile (QQ) plots and residual plots produced via a Wallyplot technique (for an overview, see Ekstrøm, 2014) indicated that the residuals appeared approximately normal and homoscedastic in distribution. There were also no concerns with multicollinearity (i.e., all variance inflation factor values ≤4.14).
Anxiety–depression group comparisons
To further explore relations of anxiety and depression in combination, comparative analyses were used to identify whether distinct subgroups of participants with syndromal or subsyndromal combinations of anxiety and depression symptoms differed on indicators of psychosocial functioning. Applying a linear approach to classifying scores on the PHQ-9 and GAD-7, participants were designated into anxiety–depression groups that varied by the severity of anxiety and depression symptoms. Participants with symptoms that met the thresholds for both anxiety (≥10) and depression (≥10) were classified into the anxiety–depression group (n = 222, 41.04%). Those with symptoms meeting the threshold for depression (≥10) and subsyndromal symptoms of anxiety (5–9) comprised the depression with subsyndromal anxiety group (n = 49, 9.06%). Conversely, the anxiety with subsyndromal depression group (n = 62, 11.46%) consisted of participants with symptoms meeting the threshold for anxiety (≥10) and subsyndromal symptoms of depression (5–9). Subsyndromal symptoms of anxiety (5–9) and depression (5–9) reflected participants with subsyndromal anxiety–depression (n = 78, 14.42%). Because we were most interested in comparing subsyndromal combinations of anxiety and depression to symptomology that aligned with standards commonly used to classify comorbidity, participants who did not meet any of the aforementioned criteria (n = 130, 24.03%) were excluded. The four groups formed the anxiety–depression group variable used in subsequent analysis.
A multivariate analysis of covariance (MANCOVA) was performed to test for mean differences between the anxiety–depression groups on the criterion variables in this study (i.e., hope, self-compassion, shame, social connectedness, and quality of social relationships). The MANCOVA was run while controlling for the same covariates as those that were included in the polynomial regression analyses. For outcome variables that reached statistical significance, follow-up analysis of covariance (ANCOVA) tests were performed. Prior to interpreting effects, Wallyplots for normality of residuals (i.e., QQ plots) and homogeneity of variance (i.e., residual plots) were inspected. No issues concerning normality and homogeneity of residuals were identified. ANCOVAs that yielded significant effects for the anxiety–depression group variable were followed by planned post hoc contrasts in which the subsyndromal anxiety–depression group was compared to each of the other three groups. To correct for multiple comparisons, a Bonferroni adjustment was applied to post hoc comparisons for each criterion variable (i.e., p values were multiplied by three).
Results
Preliminary analyses indicated that the relation between anxiety and depression was large in effect size (r = .70, p < .001). Relations of anxiety and depression with the outcome variables in this study were in the expected direction and ranged from moderate to large in effect size (r = |.33 to .54|, all p values < .001; see Table 2).
Descriptive statistics, internal consistency estimates, and zero-order correlations among study variables.
Note. All correlations are statistically significant at p < .001. SD: standard deviation.
Polynomial regression of anxiety and depression
Polynomial regression model results are reported in Table 3. Although there was no evidence of a linear association between anxiety and social relationship quality (p = .078), linear relations of anxiety (all p values ≤ .018) and depression (all p values < .001) with all other criterion variables were in the expected direction and reached statistical significance. Depression (all p values ≤ .007), but not anxiety (all p values ≥ .088), evidenced curvilinear relations with all criterion variables. For hope, self-compassion, social connectedness, and social relationship quality, curvilinear effects were positive. Two-line tests indicated that at lower levels of depression, relations with these criterion variables were negative (z = −3.98 to −5.19, all p values ≤ .001). With the exception of social connectedness (z = 2.00, p = .046), relations of these variables with depression were absent at higher values of depression (z = −1.05 to 0.76, all p values ≥ .294). The curvilinear effect for shame was negative, with positive relations evidenced at lower values of depression (z = 4.53, p < .001) but not at higher values of depression (z = −1.37, p = .172). The results did not support a multiplicative effect of anxiety and depression on any of the criterion variables (all p values ≥ .277).
Polynomial regression summary statistics.
Note. Anxiety and depression variables mean centered for all analyses. SE: standard error; CI: confidence interval.
Female = 0, male = 1.
Other race/ethnicity = 0, White = 1.
*p < .05. **p < .01. ***p < .001.
Anxiety–depression group comparisons
The results of the MANCOVA, ANCOVA, and post hoc analyses are reported in Table 4. The MANCOVA yielded a significant effect for anxiety–depression group. Follow-up ANCOVAs supported significant differences in the criterion variables as a function of anxiety–depression group. Post hoc comparisons revealed better outcomes on hope, self-compassion, shame, and quality of social relationships among participants in the subsyndromal anxiety–depression group compared to the anxiety–depression group (all p values < .001). The subsyndromal anxiety–depression group also scored significantly higher on hope than the depression with subsyndromal anxiety group (p = .008). However, the subsyndromal anxiety–depression group did not differ significantly from the depression with subsyndromal anxiety group or the anxiety with subsyndromal depression group on any other criterion variable (all p values ≥ .063). The general pattern of findings suggests that insalubrious outcomes are more likely as the severity of mixed anxiety–depression symptoms increases. However, aspects of psychosocial functioning may be similarly impaired among people with subsyndromal combinations of anxiety and depression compared to those with symptoms that coincide with syndromal definitions of comorbid anxiety–depression.
Summary statistics for MANCOVA, ANCOVA, and post hoc comparison analyses.
Note. Post hoc comparisons (using a Bonferroni adjustment) were only performed following statistically significant (p < .05) univariate analyses. SD: standard deviation; MANCOVA: multivariate analysis of covariance; ANCOVA: analysis of covariance.
Female = 0, male = 1.
Other race/ethnicity = 0, White = 1.
*p < .05. **p < .01. ***p < .001. #p < .05, compared with subsyndromal anxiety and depression group.
Discussion
In this study, a multi-analytic approach was used to examine relations of combinations of anxiety and depression with subjective intrapersonal and interpersonal indicators of psychosocial functioning. The findings of the complementary analyses were largely consistent with the expectation that increasing levels of anxiety and depression would correspond with worse psychosocial functioning, although nonlinear relations indicated that the effect of depression progressively attenuated at higher levels of symptom severity. Whereas the findings generally supported additive effects of anxiety and depression, the hypothesis that there would be synergistic effects of anxiety and depression was not supported. That is, there was no evidence that links between depression and lower quality psychosocial functioning were amplified at higher levels of anxiety.
Our findings resonate with an abundance of research that has found concurrent anxiety and depression symptoms may be more debilitating compared to the effects of either set of symptoms occurring in isolation (e.g., Lamers et al., 2011; Wanders et al., 2016). Specifically, higher levels of anxiety and depression were independently associated with unfavorable outcomes on a majority of criterion variables. Independent effects resonate with the tripartite model of anxiety and depression (see Clark & Watson, 1991; Watson et al., 1995), which posits that each share a predominant, nonspecific component of general distress, but are differentiated by elements that are specific to anxiety (i.e., physiological hyperarousal) and depression (i.e., lack of positive affect). Because anxiety and depression were simultaneously included in the regression analyses, the general pattern of findings may reflect the unique contributions of features that are specific to anxiety and depression after partialing out the underlying component of general distress that is common to both (Brown et al., 2010).
Nonlinear trends revealed that relations of depression with worse psychosocial functioning tended to dissipate as the severity of depressive symptoms increases. Whereas the effect of unique depressive symptoms may be stronger at lower levels, the negative implications of such symptoms for psychosocial functioning tend to stabilize as symptoms progressively become more severe. Although curvilinear trends may reflect ceiling or floor effects originating from skewed distributions on criterion variables (Le et al., 2011), such methodological artifacts are unlikely given skewness values were within appropriate limits and nonlinear relations were consistent across all outcomes. Considered within the context of the tripartite model, perhaps negative effects of depression-specific symptoms (i.e., low positive affect) on psychosocial functioning weaken progressively as the severity of those symptoms increases. More generally, these findings indicate that the complexities involved in identifying the psychosocial implications of concurrent anxiety and depression go beyond linear combinations of symptoms.
Group comparisons largely reaffirmed prior evidence that has found subsyndromal symptoms of mixed anxiety–depression may be associated with similar levels of impairment when compared to people with symptoms that correspond with common definitions of anxiety–depression comorbidity (see Hettema et al., 2015). These findings advocate an inclusive approach to assessing comorbid symptoms which extends beyond the narrow boundaries of finite approaches that are frequently used. This perspective aligns with a nonreductionist view of mental health and criticisms leveled against the generalized application of diagnostic criteria without accounting for broader personal and contextual factors that might affect a person’s subjective experience (Castiglioni & Laudisa, 2015). As subsyndromal anxiety and depression have each been linked to an increased likelihood of a subsequent depressive or anxiety disorder (Karsten et al., 2011), a nonreductionist approach to comorbid anxiety–depression will likely benefit early detection of individuals at risk of symptom progression and promote prevention research initiatives targeting subsyndromal combinations of anxiety and depression symptoms.
By extending our focus beyond linear effects, the findings of this study offer an improved understanding of how anxiety and depression symptoms may be nonlinearly associated with intrapersonal and interpersonal indicators of functioning. The findings also underscore the importance of integrating curvilinear effects into models in which interaction effects involving highly correlated variables are estimated (Matuschek & Kliegl, 2018), thereby safeguarding against spurious conclusions (i.e., type I error) about interaction effects (Ganzach, 1997) and enhancing the validity of statistical conclusions (García-Pérez, 2012). Nevertheless, the findings ought to be interpreted with consideration of selected methodological limitations. First, this study is based on cross-sectional data, which prevents interpretations of causality and directionality from being made. It is possible that indicators of psychosocial functioning included in this study as criterion variables precede anxiety and depression in the causal order or that relations among the variables are bidirectional. Longitudinal studies are needed to ascertain how various concurrent combinations of anxiety and depression—both at syndromal and subsyndromal levels—develop, change over time, and are affected by treatment. Second, the sample consisted of a relatively homogenous group of participants, the majority of which were of nonminority race/ethnic status and highly educated. Caution should be applied in generalizing the findings to other populations, as risk factors for, prevalence of, and outcomes linked to anxiety and depression can vary by sociodemographic (e.g., race/ethnicity) characteristics (see Watkins et al., 2015; Zhao et al., 2009). Third, supplementary comparative analyses relied on a linear approach to generate anxiety–depression subgroups using scores on the PHQ-9 and GAD-7. Although this procedure has widely been applied in lieu of formal diagnoses (Mitchell et al., 2016), subgroups were derived using subjective responses to the PHQ-9 and GAD-7 rather than clinical assessments.
Conclusions
In summary, the general pattern of findings suggests that the implications of comorbid anxiety–depression are complex and vary based on linear (i.e., anxiety) and nonlinear (i.e., depression) trends in the progression of symptom severity. The findings also highlight the importance of detecting concurrent anxiety–depression that may exist at subsyndromal levels and offering customized treatments that appropriately align with the severity of impairment linked to combinations of anxiety and depression that do not meet more traditional criteria of comorbidity.
