Introduction
Depression and self-injurious thoughts and behaviors (SITBs), such as suicidal ideation (SI), suicide attempts, and non-suicidal self-injury (NSSI), are of great concern due to related morbidity and mortality (Fox et al., Reference Fox, Franklin, Ribeiro, Kleiman, Bentley and Nock2015; Nock, Joiner, Gordon, Lloyd-Richardson, & Prinstein, Reference Nock, Joiner, Gordon, Lloyd-Richardson and Prinstein2006). Though depression and SITBs occur across the lifespan, they frequently emerge in adolescence, a time of widespread change in the brain, body, and social landscape (Baetens et al., Reference Baetens, Claes, Onghena, Grietens, Van Leeuwen, Pieters and Griffith2014; Giedd et al., Reference Giedd, Lalonde, Celano, White, Wallace, Lee and Lenroot2009). Distal and proximal biological traits have been concurrently associated with NSSI. Alteration in the HPA axis and brain structure and function in frontolimbic regions implicated in stress and pain processing, cognitive control, and social/self representation [particularly the amygdala and medial prefrontal cortex (mPFC)] have been identified as potential biomarkers for NSSI risk (Casey, Heller, Gee, & Cohen, Reference Casey, Heller, Gee and Cohen2019; Kaess et al., Reference Kaess, Hooley, Klimes-Dougan, Koenig, Plener, Reichl and Cullen2021). Importantly, frontolimbic connections strengthen during adolescence (Gee et al., Reference Gee, Humphreys, Flannery, Goff, Telzer, Shapiro and Tottenham2013), with the mPFC increasingly providing inhibitory control of the amygdala (reviewed by Kaess et al., Reference Kaess, Hooley, Klimes-Dougan, Koenig, Plener, Reichl and Cullen2021). Though research has associated adolescent stress responses with psychopathology (Kaess et al., Reference Kaess, Hooley, Klimes-Dougan, Koenig, Plener, Reichl and Cullen2021; Liu, Cheek, & Nestor, Reference Liu, Cheek and Nestor2016; Stikkelbroek, Bodden, Kleinjan, Reijnders, & Baar, Reference Stikkelbroek, Bodden, Kleinjan, Reijnders and van Baar2016), the role of concordance between components of the stress systems remains unclear. Recent work indicates that discordance between these systems may be crucial in identifying individuals at greatest risk of experiencing depression and SITBs (Bendezú, Thai, Wiglesworth, Cullen, & Klimes-Dougan, Reference Bendezú, Thai, Wiglesworth, Cullen and Klimes-Dougan2022).
Stress experience, expression, and physiology links to depression and SITBs
Stress experience
An individual's perception of stress provides insights to their appraisals of themselves and stressors in their surrounding environments, which may be crucial to determining clinical outcomes (Cohen, Kamarck, & Mermelstein, Reference Cohen, Kamarck and Mermelstein1983). Indeed, perceived stress has been positively associated with depressive symptoms (Cristóbal-Narváez, Haro, & Koyanagi, Reference Cristóbal-Narváez, Haro and Koyanagi2020), SI (Ang & Huan, Reference Ang and Huan2006; Cole et al., Reference Cole, Wingate, Tucker, Rhoades-Kerswill, O'Keefe and Hollingsworth2015; Kleiman, Turner, Chapman, & Nock, Reference Kleiman, Turner, Chapman and Nock2018), and NSSI (Kim et al., Reference Kim, Cushman, Weissman, Puzia, Wegbreit, Tone and Dickstein2015; Renteria et al., Reference Renteria, Cazares, Baltz, Schreiner, Yalcinbas, Steinkellner and Gremel2021).
Stress expression
The expression of distress has been linked to both adaptive and pathological outcomes. Expression conveys one's needs and potentially facilitates the reception of social support, thus representing adaptive response to distress (Ackerman, Abe, & Izard, Reference Ackerman, Abe, Izard, Mascolo and Griffin1998; Calhoun et al., Reference Calhoun, Patterson, Bendezú, Helms, Owens, Rudolph and Prinstein2021; Gerin, Pieper, Levy, & Pickering, Reference Gerin, Pieper, Levy and Pickering1992; Kirschbaum, Klauer, Filipp, & Hellhammer, Reference Kirschbaum, Klauer, Filipp and Hellhammer1995). However, heightened stress expression in adolescents has been longitudinally linked to greater risk for NSSI and SI (Kilic et al., Reference Kilic, Bendezú, Nick, Calhoun, Rudolph, Hastings, Nock and Prinstein2023). Downplaying expression may reflect emotional competence and efforts to attain social goals (Dickerson & Kemeny, Reference Dickerson and Kemeny2004; Mauss & Robinson, Reference Mauss and Robinson2009). However, the suppression of expressed stress responses could be deleterious, as it is associated with exaggerated physiological stress responses (Schwerdtfeger & Kohlmann, Reference Schwerdtfeger, Kohlmann, Hentschel, Smith, Draguns and Ehlers2004; Schwerdtfeger & Rathner, Reference Schwerdtfeger and Rathner2016).
Stress physiology
When stress is encountered, multiple physiological systems are engaged, including the hypothalamic–pituitary–adrenal (HPA) axis, which elicits the production of cortisol, increasing arousal and mobilizing resources to navigate the stressor (Herman, Ostrander, Mueller, & Figueiredo, Reference Herman, Ostrander, Mueller and Figueiredo2005). Patterns of heightened cortisol reactivity have been found in adolescents with depression (Lopez-Duran, Kovacs, & George, Reference Lopez-Duran, Kovacs and George2009), though those with a history of NSSI consistently demonstrate blunted cortisol reactivity (Kaess et al., Reference Kaess, Hooley, Klimes-Dougan, Koenig, Plener, Reichl and Cullen2021). Adolescents with a history of SI may have heightened cortisol reactivity (Giletta et al., Reference Giletta, Calhoun, Hastings, Rudolph, Nock and Prinstein2015); however, blunted reactivity is associated with both past and future suicide attempts (Eisenlohr-Moul et al., Reference Eisenlohr-Moul, Miller, Giletta, Hastings, Rudolph, Nock and Prinstein2018). Alterations in HPA reactivity can at times be adaptive; heightened reactivity allows for greater sensitivity to future threat, while attenuated reactivity could protect against the deleterious effects of chronically elevated cortisol (Herman et al., Reference Herman, Ostrander, Mueller and Figueiredo2005). However, these changes may increase vulnerability to psychopathology.
Interplay between stress response systems
While experience, expression, and physiological (EEP) responses to stressors are independently linked to depression, SI and attempt, and NSSI, emerging research indicates that the concordance between these systems may be implicated in psychopathology (Bendezú et al., Reference Bendezú, Thai, Wiglesworth, Cullen and Klimes-Dougan2022). It is theorized that concordance across systems indicates healthy capacity for self-regulation, while cross-system discordance indicates disruption at the system level and thus increased risk for psychopathology (Mauss & Robinson, Reference Mauss and Robinson2009; Selye, Reference Selye1936). That said, evidence also shows that HPA reactivity acts as a moderator of stress experience and depressive symptoms [e.g. stronger positive associations observed in the presence of heightened HPA activation (Owens et al., Reference Owens, Helms, Rudolph, Hastings, Nock and Prinstein2019)]. As such, heightened responsivity within and across systems may index heightened risk of psychopathology.
The investigation of interplay between stress response systems has primarily employed variable-centered methods, which assess for differences on individual stress metrics between groups (Aldao, Reference Aldao2013). Such approaches fail to capture individual-level nuance or patterns of interaction among systems, which may be critical for linking stress to outcomes of interest (e.g. Campbell & Ehlert, Reference Campbell and Ehlert2012). As a result, some recent work has attempted to examine this cross-system interplay using multi-level, person-centered frameworks (e.g. Bendezú et al., Reference Bendezú, Calhoun, Patterson, Findley, Rudolph, Hastings and Prinstein2021; Simon, Rab, Goldstein, Magal, and Admon, Reference Simon, Rab, Goldstein, Magal and Admon2022). Bendezú et al. (Reference Bendezú, Calhoun, Patterson, Findley, Rudolph, Hastings and Prinstein2021) utilized multi-trajectory modeling of adolescent girls’ negative affect, positive affect, and cortisol responses to the Trier Social Stress Test (TSST) and identified four profiles. Hyporesponsive girls, those who exhibited higher negative affect, lower positive affect, and lower cortisol production, were at increased risk for SITBs 3 to 9 months later. Notably, low cortisol on its own was not significantly associated with these risky outcomes, suggesting a need to consider multi-system patterns in context.
Recent work has employed person-centered analyses of patterns of EEP responses, permitting the identification of varying degrees of correspondence across experience, expression, and physiology systems (Bendezú et al., Reference Bendezú, Thai, Wiglesworth, Cullen and Klimes-Dougan2022). Bendezú et al. (Reference Bendezú, Thai, Wiglesworth, Cullen and Klimes-Dougan2022) examined stress responses to the TSST in adolescents diagnosed with depression and healthy controls and found three relatively concordant EEP groups and one discordant group: High Experience, High Expression, and Low Physiology (Hexperi–Hexpres–Lphysio). Similar to the Hyporesponsive profile identified in (Bendezú et al., Reference Bendezú, Calhoun, Patterson, Findley, Rudolph, Hastings and Prinstein2021), the Hexperi–Hexpres–Lphysio discordant group was associated with greater depressive symptoms, SI, and NSSI engagement. Bendezú et al. (Reference Bendezú, Thai, Wiglesworth, Cullen and Klimes-Dougan2022) also found lower frontolimbic resting-state functional connectivity (RSFC) in Hexperi–Hexpres–Lphysio compared to the concordant groups. These results suggest that discordance in the stress response system may incur risk for depression and SITBs, and that such risk may be explained by deficits in neurologic mechanisms involved in self-regulation. Further work is needed to examine this theory.
The current study: aims and hypotheses
The present study aimed to further characterize the relationships among EEP correspondence, psychopathology, and frontolimbic functioning in a sample of adolescents assigned female at birth with a range of NSSI behaviors. In Aim 1, we identified participant EEP groups in response to a laboratory-based social stressor. Guided by Bendezú et al. (Reference Bendezú, Thai, Wiglesworth, Cullen and Klimes-Dougan2022), we hypothesized that groups with varying levels of concordance would emerge. In Aim 2, we examined whether the EEP groups differed in psychopathology. We hypothesized that more discordant groups would report more severe depressive symptoms, SI and attempt, and NSSI engagement relative to the concordant groups. Finally, in Aim 3, we examined whether groups differed in neurobiological signatures associated with stress processing: mPFC-amygdala RSFC and activation of amygdala and mPFC in response to emotional stimuli. We hypothesized that more discordant groups would exhibit decreased fronto-limbic RSFC and increased amygdala and decreased mPFC activation relative to the concordant groups.
Method
Procedure
Adolescents (12–17 years) assigned female sex at birth (postmenarcheal) were recruited through flyers posted in the community, digital marketing, and local clinics and hospitals for a longitudinal study focused on NSSI and brain development. The study was oversampled for participants who engaged in NSSI behaviors, and participants represent a range of severity in NSSI engagement. Specifically, we aimed to recruit four equally sized groups of NSSI severity: no NSSI, mild NSSI, moderate NSSI, and severe NSSI (methods for assigning these groups are described in the online Supplement). Participants were excluded at intake for current substance use disorder, lifetime bipolar disorder or psychotic disorder, intellectual or developmental disability, major medical illness that would impact brain structure or function, or contraindications to undergoing magnetic resonance imaging. The present study focused on the first assessment wave, which included three visits consisting of (1) consent, clinical assessment, and self-report measures, (2) modified TSST, and (3) neuroimaging. This study was approved by the Institutional Review Board at University of Minnesota. See Başgöze and colleagues (Reference Başgöze, Mirza, Silamongkol, Hill, Falke, Thai and Cullen2021a) for additional details.
Trier social stress test- modified (TSST)
EEP responses were captured in the context of a modified TSST, during which participants were asked to give an introduction speech and perform verbal arithmetic calculations while being observed by two examiners trained to remain neutral and to avoid giving reassurance or feedback (based on a validated measure of the procedure used by (Yim, Granger, & Quas, Reference Yim, Granger and Quas2010; Yim, Quas, Rush, Granger, & Skoluda, Reference Yim, Quas, Rush, Granger and Skoluda2015); adapted from the classic TSST (Kirschbaum, Pirke, & Hellhammer, Reference Kirschbaum, Pirke and Hellhammer1993). This procedure has been repeatedly used to reliably elicit stress responses in female youth (e.g. Klimes-Dougan et al., Reference Klimes-Dougan, Begnel, Almy, Thai, Schreiner and Cullen2019). Participants were scheduled to complete the TSST between 1 and 4 pm and were instructed to not eat, drink, smoke, or brush their teeth within the hour prior to providing saliva samples.
Experience of stress
To rate the experience of stress, participants were asked ‘How stressful was…’ the (1) preparation, (2) speech task, (3) math task, and (4) post-test period (e.g. ‘now’). Responses were provided on a scale from 1 (least stressful) to 5 (most stressful). This rating scale or similar rating scales have been previously employed in other studies (e.g. Gunnar et al., Reference Gunnar, Reid, Donzella, Miller, Gardow, Tsakonas and Bendezú2021; Klimes-Dougan et al., Reference Klimes-Dougan, Begnel, Almy, Thai, Schreiner and Cullen2019).
Expression of stress
Evaluators went through a training process where they rated videos of TSST administrations to ensure reliability in ratings of expressions of stress before engaging in the ratings for this study. Evaluators were provided with a list of behavioral expressions indicative of stress, which included fidgeting, blushing, freezing, and verbal references to stress, to scaffold their responses. This rating system has been previously employed in other studies (Klimes-Dougan et al., Reference Klimes-Dougan, Begnel, Almy, Thai, Schreiner and Cullen2019). Two examiners independently rated observed expressions of stress during (1) the speech and (2) arithmetic portions of the TSST from 1 (not stressed at all) to 6 (discontinued the procedure because they were showing signs of distress), and their scores were averaged together to create separate composite mean expressed stress scores for the speech and math tasks. Similar to past studies, intraclass correlations indicated that examiner ratings demonstrated moderate agreement (ICCspeech = 0.701, p < 0.001; ICCmath = 0.651, p < 0.001).Footnote † Footnote 1
Physiological response to stress
Physiological response was captured via five salivary cortisol samples collected during the procedure: CORT1 pre-task, CORT2 immediately after the speech and math section (+15 min), CORT3 at +30 min, CORT4 at+45 min, and CORT5 at +60 min. Passive drool techniques were used and samples were stored in a −25 °C freezer until they were shipped to Universität Trier in Trier, Germany for analysis. Assay methods were consistent with Dressendörfer, Kirschbaum, Rohde, Stahl, and Strasburger (Reference Dressendörfer, Kirschbaum, Rohde, Stahl and Strasburger1992). Cortisol values that were >3 SDs from the grand mean in either direction (CORT1 (n = 3), CORT2 (n = 2), CORT3 (n = 3), CORT5 (n = 2)) were winsorized prior to analyses.
Clinical measures
Depressive symptoms
The Beck Depression Inventory, second edition (BDI; Beck, Steer, & Brown, Reference Beck, Steer and Brown1996) was used to measure current depressive symptoms. Participants rated the frequency at which they experienced 18 depressive symptoms in the past two weeks. Sum scores were computed; higher scores indicating greater depressive symptoms (skewness = 0.7, kurtosis = 2.5). When the BDI was included as a covariate in follow up analyses of our suicide ideation, suicide attempt, and NSSI models, 17-item sum scores (excluding one item that assesses suicide ideation) were used.
Suicidal ideation (SI)
The Beck Scale for Suicidal Ideation (BSSI; Beck, Kovacs, & Weissman, Reference Beck, Kovacs and Weissman1979) measured current suicidal and self-injurious thoughts, desires, and behavior patterns. Responses to the first five items, which were asked of all participants and only measure SI, were summed, with higher scores indicating greater SI (skewness = 1.7, kurtosis = 2.9).
Nonsuicidal self-injury (NSSI) and suicide attempt (SA)
The Self-Injurious Thoughts and Behaviors Interview (SITBI; Nock, Holmberg, Photos, and Michel, Reference Nock, Holmberg, Photos and Michel2007) is a clinician-administered measure of lifetime thoughts of and engagement in SITBs. The number of suicide attempts and episodes of NSSI engagement were used for each participant. Lifetime report was used given that prior history of behavior/ideation is a significant risk factor for future engagement (Fox et al., Reference Fox, Franklin, Ribeiro, Kleiman, Bentley and Nock2015). A description of the operationalization of NSSI severity and the outcomes of these analyses can be found in the online Supplement. Additionally, sensitivity analysis for past year NSSI engagement was also conducted.
Neuroimaging
Brain scans were conducted using a Siemens 3 Tesla Prisma scanner and a 32-channel receive-only head coil, using the Human Connectome Project (Glasser et al., Reference Glasser, Coalson, Robinson, Hacker, Harwell, Yacoub and Van Essen2016) multiband sequences. Structural scans were acquired using T1-weighted, multiecho Magnetization-Prepared RApid Gradient-Echo sequence and a T2-weighted Sampling Perfection with Application optimized Contrasts using different flip angle Evolution sequence. A series of functional MRI (fMRI) scans were obtained consisting of whole brain T2*-weighted functional volumes with 2 mm isotropic voxel resolution using the HCP multiband echo planar imaging sequence. Participants first underwent a resting-state scan during which they were instructed to stay awake, focus their eyes on a fixation cross, and ‘not think about anything in particular.’ Participants also completed an emotion-face matching task in which they match either emotionally expressive faces (characterized by fear or anger) or neutral shapes (Hariri, Tessitore, Mattay, Fera, & Weinberger, Reference Hariri, Tessitore, Mattay, Fera and Weinberger2002). The emotion-face matching task was conducted via E-prime software (Schneider, Eschman, & Zuccolotto, Reference Schneider, Eschman and Zuccolotto2012), and stimuli were projected onto a screen inside the bore of the scanner. The online Supplement contains greater information about neuroimaging data processing.
Resting state analyses
The CIFTI-space gray-ordinate-wise time series were used to create an average time series for each of the Glasser and Harvard-Oxford parcellations. Then, the average time series from bilateral amygdala and mPFC were extracted, cross-correlated, and Fisher's z-transformed to yield z-scores representing the RSFC between amygdala and mPFC for right and left hemispheres. The mPFC (parcellations 10v, 10r, p32, s32, 9 m, p32pr, a32pr, d32, 8BM) regions were identified using the Glasser multi-modal parcellation guide (Glasser et al., Reference Glasser, Coalson, Robinson, Hacker, Harwell, Yacoub and Van Essen2016). Values were averaged across hemispheres to obtain a bilateral value.
Task fMRI analyses
FSL FMRI Expert Analysis Tool (Woolrich, Ripley, Brady, & Smith, Reference Woolrich, Ripley, Brady and Smith2001) was used to conduct a regression analysis measuring neural activation at each gray ordinate of the brain during the emotion matching task. We included two explanatory variables (emotion face and shape conditions) of the Hariri task (Hariri et al., Reference Hariri, Tessitore, Mattay, Fera and Weinberger2002). The produced contrast of emotion > shape z-score statistical maps in matrix form were projected back to CIFTI space for further analyses. Bilateral average z-scores for gray ordinates within amygdala and mPFC for the emotion > shape contrast were extracted for further analyses.
Analytic approach
Aim 1
Following Bendezú et al. (Reference Bendezú, Thai, Wiglesworth, Cullen and Klimes-Dougan2022), multi-trajectory modeling (Nagin, Jones, Passos, & Tremblay, Reference Nagin, Jones, Passos and Tremblay2018) was used to group adolescents based on the extent to which they exhibited similar stress EEP trajectories in response to the TSST. The PROC TRAJ procedure SAS 9.4; (SAS Institute, 2010) was utilized with the MULTGROUPS option. A nonsignificant Little's (Little, Reference Little1988) MCAR test, Χ 2 (123) = 63.36, p > 0.25, supported the use of Full-Information-Maximum likelihood to handle missing data. For expression, a linear polynomial function was estimated.Footnote 2 For experience and physiology, quadratic and cubic polynomial functions were initially estimated, respectively. At each step of model specification (e.g. one-group, two-group), non-significant highest order polynomial functions for each stress index were systematically eliminated. Non-significant cubic and quadratic terms were systematically removed from trajectories at each step, while linear parameters were retained irrespective of significance (Andruff, Carraro, Thompson, Gaudreau, & Louvet, Reference Andruff, Carraro, Thompson, Gaudreau and Louvet2009; Helgeson, Snyder, & Seltman, Reference Helgeson, Snyder and Seltman2004; Louvet, Gaudreau, Menaut, Genty, & Deneuve, Reference Louvet, Gaudreau, Menaut, Genty and Deneuve2009). The log Bayes factor approximation [2loge(B10)] was utilized at each step as a fit index. Given our sample size and recommendations from the PROC TRAJ procedure developers (Jones & Nagin, Reference Jones and Nagin2007; Nagin, Reference Nagin2005), we limited model specification to four groups. Following specification, we evaluated multi-trajectory modeling adequacy via average posterior probability (AvePPj > 0.70), odds of correct classification (OCCj > 5.00), and the ratio of the probability of subgroup assignment to the proportion of adolescents assigned to subgroups ([Probj/Propj]≈1) (Nagin et al., Reference Nagin, Jones, Passos and Tremblay2018). After adequacy evaluation, Wald tests were used to distinguish the groups, delineating how intercept and polynomial functions for each trajectory were relatively ‘higher’ or ‘lower’ (e.g. baseline, reactivity) across groups.
Aims 2 and 3
A series of regression analyses were employed to determine whether differences were observed between EEP groups in clinical and neurobiological measures of interest. In keeping with Bendezú et al. (Reference Bendezú, Thai, Wiglesworth, Cullen and Klimes-Dougan2022), we compared each group to the Hexperi–Hexpres–Lphysio group in all analyses. Linear regressions were run for linear dependent variables (e.g. depressive symptoms, SI, RSFC, regional activation) and negative binomial regressions were run for count dependent variables (e.g. NSSI, SA). EEP group membership, which was dummy coded when entered into the models (with Hexperi–Hexpres–Lphysio left out so that it would be the default reference group), was used as the primary predictor. Age, minority race or ethnicity (as opposed to white),Footnote 3 family income (<$ 60 000 v. ⩾$ 60 000), and HPA medication status (Granger, Hibel, Fortunato, & Kapelewski, Reference Granger, Hibel, Fortunato and Kapelewski2009) were included as covariates in all analyses. Follow up analyses were conducted for SI, SA, and NSSI controlling for depressive symptoms. We checked assumptions of normality for linear outcomes prior to analysis.
Results
Our final sample included 109 adolescent participants with complete EEP data (see Table 1 for descriptives of the study sample; sample sizes for each analysis are reported in the corresponding results table).
Participant demographic and clinical descriptives disaggregated by profiles
Note: Lexperi–Lexpres–Lphysio, Low experience, low expression, low physiology profile group; Hexperi–Hexpres–Hphysio, High experience, high expression, high physiology profile group; Hexperi–Lexpres–Mphysio, High experience, low expression, moderate physiology profile group; Hexperi–Hexpres–Lphysio, High experience, High expression, low physiology profile group; BDI, Beck Depression Inventory; BSSI, Beck Scale for Suicide Ideation; NSSI, non-suicidal self-injury; SA, suicide attempt; NSSI History, SA History, The number and percentage of individuals within each group that report any history of NSSI or suicide attempt; NSSI Incidents, SA Incidents, The average and standard deviation of the lifetime number of engagements in NSSI or SA reported by each group.
Aim 1: stress experience–expression–physiology groups
Multi-trajectory modeling parameter estimates, adequacy indices, and trajectory distinction analysis results are shown in Table 2. Model specification revealed four groups (Fig. 1): two- to one-group comparison [2loge(B10) = 263.66], three- to two-group comparison [2loge(B10) = 56.88], four- to three-group comparison [2loge(B10) = 10.46]. As per Nagin (Reference Nagin2005), the final model fit the data well. With multi-trajectory modeling, subgroup trajectories may differ in intercept but not response patterning (or vice versa). In these cases, we used that significantly differing aspect to help inform our characterization of the subgroups.
Parameter estimates (standard errors) and model adequacy indices for subgroups
Note: AvePPj, Average posterior probability; OCCj, Odds of correct classification; Probj, Probability of group assignment; Propj, Proportion of children assigned to each group; Ratio, Ratio of Probj to Propj; Upper-case and lower-case superscripts denote significant differences in intercept and polynomial estimates, respectively, within the same Trier Social Stress Test response index.
†p < 0.10 *p < 0.05.
Experience, expression, and physiological trajectories of subgroups.
Note: Physiology response was indexed by salivary cortisol (μg/dL). TSST, Trier Social Stress Test.
Three subgroups emerged whose stress response groups were generally consistent with those identified in Bendezú et al. (Reference Bendezú, Thai, Wiglesworth, Cullen and Klimes-Dougan2022): Low Experience–Low Expression–Low Physiology (Lexperi–Lexpres–Lphysio, n = 41), High Experience–High Expression–Low Physiology (Hexperi–Hexpres–Lphysio; n = 28), High Experience–High Expression–High Physiology (Hexperi–Hexpres–Hphysio, n = 19). A fourth subgroup emerged that was unique to our study sample; High Experience–Low Expression–Moderate Physiology (Hexperi–Lexpres–Mphysio, n = 25). There were no significant differences between these four groups in age, proportion of white v. Minority race/ethnicity youth, or gross family income (all p's > 0.05).
Aim 2: clinical correlates of stress EEP groups
As shown in Table 3, compared to the Lexperi–Lexpres–Lphysio and the Hexperi–Lexpres–Mphysio subgroups, the Hexperi–Hexpres–Lphysio group exhibited higher depressive symptoms and SI severity. Compared to the Hexperi–Hexpres–Hphysio and the Hexperi–Lexpres–Mphysio subgroups, those in the Hexperi–Hexpres–Lphysio group had higher odds of reporting greater lifetime SAs (Table 3). While there were no group differences in the episode number (severity) or in dichotomous indexes of lifetime NSSI engagement, there were some results consistent with predictions when restricting the NSSI engagement to past year. That is, those in the Hexperi–Hexpres–Lphysio group demonstrated significantly higher odds of greater NSSI engagements compared to those in the Hexperi–Lexpres–Mphysio group (Table S1).
Stress Experience–Expression–Physiology profiles and clinical outcomes
Note: b (se), unstandardized beta (standard error); B, standardized beta; BDI, Beck Depression Inventory; BSSI, Beck Scale for Suicide Ideation; NSSI, nonsuicidal self-injury; OR, odds ratio; 95% CI, 95% confidence interval; Lexperi–Lexpres–Lphysio, Low experience, low expression, low physiology profile group; Hexperi–Hexpres–Hphysio, High experience, high expression, high physiology profile group; Hexperi–Lexpres–Mphysio, High experience, low expression, moderate physiology profile group; Gross Family Income, below median (<$60 000) or above median (≥$60 000 = 1) (below = 0, above = 1); Suicide Attempt and NSSI are negative binomial count regression models; aReference group for effects is High Experience-High Expression-Low Physiology; bResult changes in significance when covarying for suicidal ideation; * p < 0.05.
Results did change in some cases in follow up analyses of SI severity, SA, and NSSI that controlled for depressive symptoms (Table S2). Specifically, depressive symptoms were the only significant predictor of SI severity when incorporated into the models. While group differences in lifetime SAs between the Hexperi–Hexpres–Hphysio and Hexperi–Hexpres–Lphysio groups were maintained, differences between the Hexperi–Lexpres–Mphysio and Hexperi–Hexpres–Lphysio were no longer significant.
Aim 3: neural correlates of stress EEP groups
Parameter estimates for our regression models are shown in Table 4. There were no significant differences between the Hexperi–Hexpres–Lphysio group and the other three groups with respect to amygdala–mPFC RSFC nor amygdala or mPFC activation to the emotional face-matching task.
Stress Experience–Expression–Physiology profiles and neural outcomes
Note: b (se), unstandardized beta (standard error); B, standardized beta; RSFC, resting state functional connectivity; mPFC, medial prefrontal cortex; Lexperi–Lexpres–Lphysio, Low experience, low expression, low physiology profile group; Hexperi–Hexpres–Hphysio, High experience, high expression, high physiology profile group; Hexperi–Lexpres–Mphysio, High experience, low expression, moderate physiology profile group; Gross Family Income, below median (<$60 000) or above median (≥$60 000 = 1) (below = 0, above = 1).
aReference group for effects is High experience-High expression-Low physiology; * p < 0.05.
Discussion
The present study extends evidence that discordance within EEP responses to stress are related to psychopathology and that these associations are more nuanced than can be ascertained using variable-centered approaches. Leveraging a unique sample of female adolescents at heightened risk for SITBs, we examined multi-level data (self-reported stress and psychopathology, observer reported stress, salivary cortisol, and neuroimaging) to assess whether within-person patterns across stress response indices were related to psychopathology and neural indices of activation and self-regulation. In taking a person-centered analytical approach, we identified four EEP groups: Hexperi–Hexpres–Lphysi, Hexperi–Hexpres–Hphysi, Hexperi–Lexpres–Mphysi, Lexperi–Lexpres–Lphysi. These EEP groups differed in their reports of depression, SI, suicide attempt, and NSSI (only for episodes reported in the past year), helping to further our understanding of the relationship between EEP discordance and risk for psychopathology.
The links between person-centered groups and psychopathology are illuminating; adolescents in the Hexperi–Hexpres–Lphysio group reported more severe depression and SI compared to those in the Lexperi–Lexpres–Lphysio and Hexperi–Lexpres–Mphysio groups and reported more frequent engagement in suicide attempts compared to those in the Lexperi–Lexpres–Lphysio and Hexperi–Hexpres–Hphysio groups. In a sample in which about two-thirds of participants had a history of NSSI engagement, no differences between groups were found in lifetime NSSI engagements, likely in part due to the pervasive shared risk factor across groups. Nevertheless, more proximal indexes of NSSI (past year) revealed group differences. Such patterns of elevated risk among youth in the Hexperi–Hexpres–Lphysio group are broadly consistent with previous research, which has demonstrated that high psychological stress reactivity confers risk for SI for youth with low physiological stress responsivity (Eisenlohr-Moul et al., Reference Eisenlohr-Moul, Miller, Giletta, Hastings, Rudolph, Nock and Prinstein2018; Owens et al., Reference Owens, Helms, Rudolph, Hastings, Nock and Prinstein2019). Moreover, as shown here, a person-centered approach is critically important in clarifying differences for those who exhibit low physiological reactivity. Though youth in the Lexperi–Lexpres–Lphysio and Hexperi–Hexpres–Lphysi groups exhibited similar physiological reactivity, the groups were qualitatively distinct with respect to reports of psychopathology. Thus, it is likely that research that considers only one indicator of stress response in isolation obscures system-level patterns of vulnerability. Further, this person-centered approach highlighted nuance within youth who exhibited high distress during the TSST. While youth in the Hexperi–Hexpres–Hphysio group did not significantly differ from those in the Hexperi–Hexpres–Lphysio group with respect to depressive symptoms, suicide ideation, and NSSI engagement, these concordantly high responding youth did report significantly less frequent engagement in suicide attempts than those in the discordant Hexperi–Hexpres–Lphysio group, potentially highlighting some protective features of this concordant response profile.
Importantly, the present study provided a partial-replication and extension of the work from Bendezú et al. (Reference Bendezú, Thai, Wiglesworth, Cullen and Klimes-Dougan2022) in a unique sample enriched for SITBs. Despite differences in samples and outcomes, the overall replication of EEP groups and associations with risk in an independent sample that is more focused on transdiagnostic risk factors for maladaptive outcomes underscores the robustness of these patterns. Conceptual replication is critical for scientific advances (Pashler & Wagenmakers, Reference Pashler and Wagenmakers2012). This is particularly true when employing person-centered approaches, as this approach carries with it concerns about whether identified groups are statistical artifacts or a result of over-fitting to one study sample. Our replication helps to allay those concerns, and tentatively suggests that EEP discordance may represent a reliable source of meaningful stress response heterogeneity that can be dependably captured using multi-level data and person-centered techniques such as multi-trajectory modeling.
Our study did not produce a full replication of the results from Bendezú et al. (Reference Bendezú, Thai, Wiglesworth, Cullen and Klimes-Dougan2022). In our study, no significant differences were found between the identified EEP groups in neurobiological function. These differences in results may reflect differences in characteristics of the study samples, of which there are multiple. First, the present transdiagnostic study oversampled for a history of NSSI while Bendezú et al. (Reference Bendezú, Thai, Wiglesworth, Cullen and Klimes-Dougan2022) employed a case-control design based on depressive disorder diagnosis. Though these clinical constructs frequently overlap, there are significant differences in the neurobiological and psychosocial mechanisms at play (Başgöze, Wiglesworth, Carosella, Klimes-Dougan, & Cullen, Reference Başgöze, Wiglesworth, Carosella, Klimes-Dougan and Cullen2021b). Second, the present sample was made up entirely of individuals assigned female sex at birth. Bendezú et al. (Reference Bendezú, Thai, Wiglesworth, Cullen and Klimes-Dougan2022) sample was 66% female, with adolescents in the Mexperi–Mexpres–Mphysio subgroup more likely to be male than female, potentially accounting for the absence of this subgroup in the current study. Additionally, the present sample was younger (mean age of 14.6 years) than that of Bendezú et al. (Reference Bendezú, Thai, Wiglesworth, Cullen and Klimes-Dougan2022) (mean age 16.5 years). As adolescence is a time of vast neural change, particularly in frontolimbic regions (Gee et al., Reference Gee, Humphreys, Flannery, Goff, Telzer, Shapiro and Tottenham2013; Power, Fair, Schlaggar, & Petersen, Reference Power, Fair, Schlaggar and Petersen2010), our hypothesized neural associations may emerge at later stages of development.
Additionally, while not considered in past research, concordance in EEP did not appear to be related to stress responsivity of the amygdala and mPFC. These findings contrast prior work that showed that left amygdala activation is related to the experience of stress during the TSST (Başgöze et al., Reference Başgöze, Mirza, Silamongkol, Hill, Falke, Thai and Cullen2021a). While it is possible that this relationship may emerge if measurements of EEP and amygdala/mPFC activation and connectivity were collected using synchronous evaluations, these results indicate a need for more work to understand the interplay between physiological and neural systems as related to stress system functioning.
Limitations
Our results should be considered in the context of the study limitations. First, while multi-trajectory modeling analyses have been used with smaller samples (Simon et al., Reference Simon, Rab, Goldstein, Magal and Admon2022), a larger sample would increase statistical power of our analyses. Second, there are considerations regarding the representativeness of our study sample and the generalizability of our results. This sample was majority white and higher-than-average socioeconomic status. Thus, our study failed to capture variation which may be expected due to diversity in racial and ethnic backgrounds, such as experiences of minority stress and historical oppression (Blair & Raver, Reference Blair and Raver2016; Flentje et al., Reference Flentje, Obedin-Maliver, Lubensky, Dastur, Neilands and Lunn2020). Additionally, the sample consisted solely of individuals assigned female sex at birth; it is important to note that stress processing differs as a function of biological sex (Taylor et al., Reference Taylor, Lewis, Gruenewald, Gurung, Updegraff and Klein2002). As the majority of the present findings replicated those from a sample of both male and female youth, the results may still be generalizable. The present sample also excluded adolescents with substance use disorder, bipolar disorder, and psychotic symptoms. These exclusion criteria were intended to decrease the variability in biological metrics due to factors unrelated to NSSI, though they may limit the generalizability of findings for all engaging in NSSI. Third, SA and NSSI history as well as experience and expression of stress during the TSST were captured using single item measures. Greater validity could be achieved through more thorough probing of different aspects of these phenomena in the form of multi-item assessments (e.g. Ammerman, Burke, Jacobucci, and McClure, Reference Ammerman, Burke, Jacobucci and McClure2021). Fourth, this study did not capture adequate metrics of, or analyze, menstrual phase, nicotine use, or engagement in physical activity, all of which have been shown to impact HPA activity and may contribute to error in group assignment. Finally, we are unable to ascertain whether these findings index correlates or causes of risk for pathology. Considering longitudinal patterns of clinical constructs of interest will be an important next step to delineate the developmental course of EEP groups (Bendezú et al., Reference Bendezú, Calhoun, Patterson, Findley, Rudolph, Hastings and Prinstein2021).
Clinical implications
Our findings suggest that examining concordance across multiple levels of stress system functioning can give insight into potentially heightened risk for depression and suicidal thoughts and behaviors, illuminating points of intervention. Based on the present results, blunted physiological reactivity in the presence of heightened experience and expression of stress may represent a risk factor for psychopathology such as depression and suicidal thoughts and behaviors. Adolescents displaying EEP discordance, especially in high-risk contexts, could be targeted with interventions supporting coping and stress appraisal, which could theoretically support multi-level integration and recalibration across systems (Calvete, Orue, Fernández-González, Chang, & Little, Reference Calvete, Orue, Fernández-González, Chang and Little2020; Kallianta et al., Reference Kallianta, Katsira, Tsitsika, Vlachakis, Chrousos, Darviri and Bacopoulou2021). Thus far, the limited work that has focused on modulating the cortisol response in youth has employed a vast array of measurement and treatment modalities and produced mixed findings (Slopen, McLaughlin, & Shonkoff, Reference Slopen, McLaughlin and Shonkoff2014).
Conclusion
The present study shows evidence that multi-level, person-centered analyses of stress response patterns offer unique insight into how interplay across elements of the stress response system within groups of individuals may confer risk for, or offer protection against, psychopathology. Importantly, results indicate that the discordant group with high expression and experienced stress and low cortisol output to stress demonstrated elevated depression, suicidal thoughts, and suicide attempts. Findings advance the EEP literature by demonstrating that unique patterns of stress responses are related to nuanced patterns of risk. Thus, interventions that attend to these heightened signs of distress may be particularly useful for identifying and supporting adolescents at risk for depression and SITBs.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S0033291723002003.
Financial support
Funding for this project provided by the National Institute for Mental Health R01MH107394. Author AW is supported by the National Science Foundation Graduate Research Fellowship Program.
Ethical standard
The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008.
Open access
