Hostname: page-component-78c5997874-g7gxr Total loading time: 0 Render date: 2024-11-03T12:10:33.118Z Has data issue: false hasContentIssue false

Resting-state connectivity subtype of comorbid PTSD and alcohol use disorder moderates improvement from integrated prolonged exposure therapy in Veterans

Published online by Cambridge University Press:  30 April 2021

Daniel M. Stout*
Affiliation:
Center of Excellence for Stress and Mental Health, VA San Diego Healthcare System, San Diego, CA, USA Research Service, VA San Diego Healthcare System, San Diego, CA, USA Department of Psychiatry, University of California San Diego, San Diego, CA, USA
Katia M. Harlé
Affiliation:
Research Service, VA San Diego Healthcare System, San Diego, CA, USA Department of Psychiatry, University of California San Diego, San Diego, CA, USA
Sonya B. Norman
Affiliation:
Center of Excellence for Stress and Mental Health, VA San Diego Healthcare System, San Diego, CA, USA Research Service, VA San Diego Healthcare System, San Diego, CA, USA Department of Psychiatry, University of California San Diego, San Diego, CA, USA National Center for PTSD, White River Junction, Vermont, USA
Alan N. Simmons
Affiliation:
Center of Excellence for Stress and Mental Health, VA San Diego Healthcare System, San Diego, CA, USA Research Service, VA San Diego Healthcare System, San Diego, CA, USA Department of Psychiatry, University of California San Diego, San Diego, CA, USA
Andrea D. Spadoni
Affiliation:
Center of Excellence for Stress and Mental Health, VA San Diego Healthcare System, San Diego, CA, USA Research Service, VA San Diego Healthcare System, San Diego, CA, USA Department of Psychiatry, University of California San Diego, San Diego, CA, USA
*
Author for correspondence: Daniel M. Stout, E-mail: dastout@health.ucsd.edu
Rights & Permissions [Opens in a new window]

Abstract

Background

Posttraumatic stress disorder (PTSD) and alcohol use disorder (AUD) are highly comorbid and are associated with significant functional impairment and inconsistent treatment outcomes. Data-driven subtyping of this clinically heterogeneous patient population and the associated underlying neural mechanisms are highly needed to identify who will benefit from psychotherapy.

Methods

In 53 comorbid PTSD/AUD patients, resting-state functional magnetic resonance imaging was collected prior to undergoing individual psychotherapy. We used a data-driven approach to subgroup patients based on directed connectivity profiles. Connectivity subgroups were compared on clinical measures of PTSD severity and heavy alcohol use collected at pre- and post-treatment.

Results

We identified a subgroup of patients associated with improvement in PTSD symptoms from integrated-prolonged exposure therapy. This subgroup was characterized by lower insula to inferior parietal cortex (IPC) connectivity, higher pregenual anterior cingulate cortex (pgACC) to posterior midcingulate cortex connectivity and a unique pgACC to IPC path. We did not observe any connectivity subgroup that uniquely benefited from integrated-coping skills or subgroups associated with change in alcohol consumption.

Conclusions

Data-driven approaches to characterize PTSD/AUD subtypes have the potential to identify brain network profiles that are implicated in the benefit from psychological interventions – setting the stage for future research that targets these brain circuit communication patterns to boost treatment efficacy.

Type
Original Article
Creative Commons
This is a work of the US Government and is not subject to copyright protection within the United States. Published by Cambridge University Press
Copyright
Copyright © United States Department of Veterans Affairs, 2021

Introduction

Rates of comorbid posttraumatic stress disorder and alcohol use disorder (PTSD/AUD) are extremely high (Blanco et al., Reference Blanco, Xu, Brady, Perez-Fuentes, Okuda and Wang2013). PTSD/AUD comorbidity is associated with greater symptom severity and functional impairment than either disorder alone (Norman, Haller, Hamblen, Southwick, & Pietrzak, Reference Norman, Haller, Hamblen, Southwick and Pietrzak2018). Although several studies show that gold-standard trauma-focused psychotherapies are efficacious interventions for some patients with comorbid PTSD/AUD (Norman et al., Reference Norman, Trim, Haller, Davis, Myers, Colvonen and Mayes2019; Roberts, Roberts, Jones, & Bisson, Reference Roberts, Roberts, Jones and Bisson2015; Ruglass et al., Reference Ruglass, Lopez-Castro, Papini, Killeen, Back and Hien2017), they are associated with high-dropout rates and a significant number of individuals maintaining elevated symptoms after treatment (Krystal et al., Reference Krystal, Davis, Neylan, Raskind, Schnurr, Stein and Huang2017; Steenkamp, Reference Steenkamp2016). One potential reason for modest outcomes is that PTSD and AUD are highly heterogeneous – suggesting that they are not distinct disorders but are comprised of sets of neurobiological mechanisms that have both unique and shared pathophysiology. Under this framework, there is a need to construct neurobiologically-grounded subtypes focused on dysfunctional brain networks that capture the heterogeneity (Feczko et al., Reference Feczko, Miranda-Dominguez, Marr, Graham, Nigg and Fair2019; Huys, Maia, & Frank, Reference Huys, Maia and Frank2016; Yamada et al., Reference Yamada, Hashimoto, Yahata, Ichikawa, Yoshihara, Okamoto and Kawato2017). Subtyping patients on shared brain network features instead of clinical symptoms could inform the mechanistic understanding of PTSD/AUD pathophysiology and provide clinically-useful biomarkers that are associated with treatment outcomes (Milham, Craddock, & Klein, Reference Milham, Craddock and Klein2017; Woo, Chang, Lindquist, & Wager, Reference Woo, Chang, Lindquist and Wager2017).

Mapping brain network communication has been instrumental in characterizing the pathophysiology of PTSD and AUD in addition to understanding how network connectivity patterns relate to treatment success (Sheynin et al., Reference Sheynin, Duval, King, Angstadt, Phan, Simon and Liberzon2020; Zhutovsky et al., Reference Zhutovsky, Thomas, Olff, van Rooij, Kennis, van Wingen and Geuze2019). A common method to measure brain network communication is through the use of functional connectivity (Abi-Dargham & Horga, Reference Abi-Dargham and Horga2016). Such work has revealed that PTSD and AUD are commonly observed to involve dysfunctional neural connectivity within and between several regions implicated in affective reactivity, cognition, and emotion regulation (Akiki, Averill, & Abdallah, Reference Akiki, Averill and Abdallah2017; Koch et al., Reference Koch, van Zuiden, Nawijn, Frijling, Veltman and Olff2016; Sheynin & Liberzon, Reference Sheynin and Liberzon2017). For example, increased pregenual anterior cingulate cortex (pgACC) connectivity with the amygdala is associated with higher PTSD symptoms after completing a trauma-focused psychotherapy (Cisler et al., Reference Cisler, Sigel, Kramer, Smitherman, Vanderzee, Pemberton and Kilts2015) and increased connectivity of the insula with the posterior midcingulate cortex (pMCC) is observed in PTSD patients who no longer met the criteria for PTSD post-treatment (Simmons, Norman, Spadoni, & Strigo, Reference Simmons, Norman, Spadoni and Strigo2013). In AUD, decreased connectivity between regions of the prefrontal cortex (PFC) and the subcortical network consisting of reward regions (e.g. ventral striatum) is associated with relapse (Camchong, Stenger, & Fein, Reference Camchong, Stenger and Fein2013) and increased connectivity between the PFC, striatum, and the insula is associated with residential treatment program completion (Kohno, Dennis, McCready, & Hoffman, Reference Kohno, Dennis, McCready and Hoffman2017).

A limitation of these investigations is that they do not consider the clinical and neural heterogeneity within these disorders. To date, the majority of neuroimaging studies exclude individuals with commonly comorbid conditions, limiting the generalizability and translation of their findings (Blanco et al., Reference Blanco, Wall, Lindquist, Rodríguez-Fernández, Franco, Wang and Olfson2016). Instead, data-driven approaches that characterize highly representative patients on functional connectivity patterns show promise in identifying neurobiological subtypes within broader heterogeneous diagnostic groups and in predicting treatment response (Drysdale et al., Reference Drysdale, Grosenick, Downar, Dunlop, Mansouri, Meng and Liston2017; Etkin et al., Reference Etkin, Maron-Katz, Wu, Fonzo, Huemer, Vértes and O'Hara2019; Maron-Katz et al., Reference Maron-Katz, Zhang, Narayan, Wu, Toll, Naparstek and Etkin2019).

Methods to improve the identification of the complex and directional neural patterns underlying psychiatric heterogeneity are needed (Reid et al., Reference Reid, Headley, Mill, Sanchez-Romero, Uddin, Marinazzo and Cole2019). Effective connectivity-derived biotypes have been shown to be superior in differentiating PTSD patients from healthy controls and better predict symptom severity than static connectivity measures (Jin et al., Reference Jin, Jia, Lanka, Rangaprakash, Li, Liu and Deshpande2017; Li et al., Reference Li, Zhu, Jiang, Jin, Zhang, Guo and Liu2014). One novel data-driven effective connectivity strategy called group iterative multiple model estimation with subgrouping (S-GIMME) has been shown to be particularly successful in deriving reliable and parsimonious directed connectivity and in separating clinically heterogeneous populations into distinct subtypes, even in samples as low as 25 (Beltz & Gates, Reference Beltz and Gates2017; Gates, Lane, Varangis, Giovanello, & Guskiewicz, Reference Gates, Lane, Varangis, Giovanello and Guskiewicz2017). A recent review of network modeling approaches shows that GIMME is receiving a growing interest as an analytical tool to understand the causal relationships between brain networks (Robinaugh, Hoekstra, Toner, & Borsboom, Reference Robinaugh, Hoekstra, Toner and Borsboom2020). S-GIMME applied to rs-fMRI has been fruitful in identifying subtypes in cognitive task ability (McCormick & Telzer, Reference McCormick and Telzer2018), accurately distinguishing depressed patients from healthy individuals (Price et al., Reference Price, Lane, Gates, Kraynak, Horner, Thase and Siegle2017b), and revealing major depressive disorder, attention-deficit/hyperactivity disorder, and transdiagnostic internalizing psychiatric disorder and personality subtypes (Chahal et al., Reference Chahal, Weissman, Hallquist, Robins, Hastings and Guyer2020; Lazarus, Sened, & Rafaeli, Reference Lazarus, Sened and Rafaeli2020; Price, Gates, Kraynak, Thase, & Siegle, Reference Price, Gates, Kraynak, Thase and Siegle2017a; Price et al., Reference Price, Beltz, Woody, Cummings, Gilchrist and Siegle2020). However, no studies to date have leveraged this data-driven approach to identify neural subtypes of comorbid PTSD/AUD individuals or to treatment response (Gilpin & Weiner, Reference Gilpin and Weiner2017).

To this end, we applied S-GIMME to rs-fMRI in a sample of 53 psychotherapy treatment-seeking Veterans diagnosed with comorbid PTSD and AUD participating in a randomized clinical trial examining treatment effectiveness of trauma-focused and non-trauma-focused therapies (Norman et al., Reference Norman, Trim, Haller, Davis, Myers, Colvonen and Mayes2019). Connectivity-based subtyping was conducted using seven regions of interest (ROIs) that were selected from a recent meta-analysis showing brain regions, including the cingulate cortex, insula, and the PFC, that distinguished PTSD and AUD patients from healthy controls (Klaming et al., Reference Klaming, Harlé, Infante, Bomyea, Kim and Spadoni2019). The aim of our investigation was to extend these findings by examining whether pre-treatment effective connectivity of these PTSD and AUD-relevant brain regions can be used to identify subgroups of patients and to test whether effective connectivity subgroups predict psychotherapeutic treatment improvement in Veterans.

Materials and methods

Participants

A total of 59 subjects (mean age = 40.23, s.d. = 10.73; range = 25–63; 47 men) from a parent psychotherapy treatment trial (NCT01601067) (Norman et al., Reference Norman, Trim, Haller, Davis, Myers, Colvonen and Mayes2019) comparing trauma-focused [integrated-prolonged exposure (I-PE)] to non-trauma-focused [integrated-coping skills (I-CS)] psychotherapy for comorbid PTSD/AUD agreed to complete neuroimaging prior to being assigned to a treatment arm (see Table 1 for sample characteristics). I-PE and I-CS therapy arms were delivered in 90 min individual sessions (12–16 sessions for both therapy arms; Norman et al., Reference Norman, Trim, Haller, Davis, Myers, Colvonen and Mayes2019). Subjects completed an average of 10.32 (s.d. = 4.33) therapy sessions [I-PEmean = 9.10, I-PESD = 4.35; I-CSmean = 11.79, I-CSSD = 3.89; t(50.67) = 2.37, p = 0.02]. Five subjects did not complete scanning due to technical and scheduling issues and one subject was removed due to do excessive motion artifact (>0.3 mm mean FD) measured via frame-wise displacement, leaving a total of 53 pretreatment resting-state scans. Thirty-eight subjects were retained at post-treatment. Subjects provided informed consent and were monetarily compensated for their participation. Study procedures were approved by the institutional review board of the Veteran's Administration San Diego Healthcare System.

Table 1. Descriptive statistics and sample characteristics

Note. Data reported are the raw descriptive statistics. No connectivity subgroup differences were found on descriptive measures. n pre, number of subjects at pre-treatment; n post, number of subjects at post-treatment; CAPS-5, Clinician Administered PTSD Scale for DSM-5.

Clinical measures

Clinician Administered PTSD Scale for DSM-5 (CAPS-5). The CAPS-5 (score range 0–80, with 0 indicating no PTSD symptoms and 80 indicating extreme ratings across all symptoms), a 30-item structured interview considered to be the criterion standard for PTSD, was the primary measure of PTSD symptoms and diagnosis. The CAPS-5 diagnosis has been shown to demonstrate strong interrater reliability (κ = 0.78), and severity scores had strong internal consistency (α = 0.88) (Weathers et al., Reference Weathers, Bovin, Lee, Sloan, Schnurr, Kaloupek and Marx2018).

Percentage of heavy drinking days (PHDD). Frequency and quantity of alcohol use were assessed using the Timeline Follow-Back, a calendar-assisted structured clinical interview that displays good psychometric properties (Sobell & Sobell, Reference Sobell, Sobell, Litten and Allen1992). The PHDD was calculated by dividing the number of days in which five or more drinks for men or five or more drinks for women were consumed by the total number of days in the reference period (i.e. 90 days at baseline and time from baseline assessment for post-treatment assessment).

fMRI acquisition and preprocessing

MRI acquisition and preprocessing are detailed in the Supplementary Materials.

Regions of interest

Five to 15 ROIs are recommended when using S-GIMME for identifying directed connectivity paths and sample subtyping (https://gimme.web.unc.edu/63-2/gimme-basics/). Therefore, we chose to take a hypothesis-driven approach for ROI selection based upon a recent meta-analysis from our group that identified seven key brain regions important for distinguishing between PTSD, AUD, their overlap, and healthy control populations (Klaming et al., Reference Klaming, Harlé, Infante, Bomyea, Kim and Spadoni2019). Our aim was to extend this work by determining whether effective connectivity of these regions can also identify subgroups of comorbid PTSD/AUD patients and whether these subgroups are clinically relevant. To this end, we extracted mean BOLD time-series from seven spherical ROIs with a 10 mm radius. ROIs were placed in three regions of the rostral cingulate cortex (see Shackman et al., Reference Shackman, Salomons, Slagter, Fox, Winter and Davidson2011 for detailed anatomical and functional description of the cingulate cortex): one in the pMCC (x = 0, y = −4, z = 34), one in the pgACC (x = 2, y = 42, z = 20), and one in the subgenual cingulate cortex (sgACC; x = 0, y = 6, z = −6). One ROI was placed in the left middle frontal gyrus (MFG; x = −46, y = 10, z = 38). One was placed in the inferior parietal cortex (IPC) near the supramarginal gyrus (x = 52, y = −46, z = 44). One in the left insula (x = −38, y = −4, z = 6) and one in the right putamen (x = 34, y = −4, z = 0). The extracted time-series from these seven ROIs were then entered into the connectivity clustering analysis described below.

Data-driven subgroup clustering via effective connectivity

S-GIMME in R (package: gimme) was used to identify whether clinically relevant subgroups of subjects can be derived based upon pretreatment resting-state effective connectivity patterns (Gates et al., Reference Gates, Lane, Varangis, Giovanello and Guskiewicz2017) (online Supplementary Fig. S1). S-GIMME is a unified structural equation modeling strategy that reliably recovers both the presence and direction of connectivity paths within heterogeneous individuals (Beltz & Gates, Reference Beltz and Gates2017; Gates, Molenaar, Iyer, Nigg, & Fair, Reference Gates, Molenaar, Iyer, Nigg and Fair2014). S-GIMME can be used to identify common pathways across all subjects but importantly can classify subjects within subgroups based upon unique connectivity profiles and estimate subject-level connectivity (Gates et al., Reference Gates, Lane, Varangis, Giovanello and Guskiewicz2017).

First, S-GIMME enters all the ROIs in a single step to mitigate the problem of variable order during the model-building procedure. Next, compared to a null model that assumes no connections exist, S-GIMME searches across individuals and adds lagged (lag = 1 TR) and contemporaneous directed paths that are common for the majority of subjects (group-level). S-GIMME uses modification indices to count individuals whose connectivity paths exist across individuals and only includes the path if it improves model fit for the majority of the subjects (75%) and must survive a conservative Bonferroni correction of 0.05/number of subjects (N = 53 in the current analysis). Paths not meeting these criteria are pruned. This process is iteratively repeated for each connectivity path and stops adding paths when no new path would improve the model-fit indices in the majority of the sample. If any remaining connection no longer improves model fit during this iterative process, then that path will be removed from the model. This pruning step can occur at the individual, sub-group, or group-level model estimation process (Beltz & Gates, Reference Beltz and Gates2017). This approach is robust to outliers and sign differences of parameter estimates (Gates et al., Reference Gates, Lane, Varangis, Giovanello and Guskiewicz2017). Next, S-GIMME uses Walktrap unsupervised community detection (Gates, Henry, Steinley, & Fair, Reference Gates, Henry, Steinley and Fair2016; Orman, Labatut, & Cherifi, Reference Orman, Labatut and Cherifi2012) on the similarity matrix derived from the individual-level estimates of the group-level connectivity paths to identify an optimal number of subgroups who have similar connectivity profiles. Once subgroups are identified, S-GIMME searches for unique connectivity paths for each data-driven subgroup. Finally, S-GIMME identifies individual-level connections using the group- and subgroup-derived temporal patterns as priors (not reported here).

S-GIMME stops adding new paths/connections once the model achieves good fit on at least two out of the following four model fit indices (Brown, Reference Brown2006): confirmatory fit index (CFI) ⩾0.95; non-normed fit index ⩾0.95; standardized root mean square residual (SRMR) ⩽0.05; root mean square error of approximation ⩽0.05 (Gates & Molenaar, Reference Gates and Molenaar2012). Stopping the S-GIMME path search based on model fit instead of stopping based on individual path significance ensures that extraneous connectivity paths are not added, a critical step to minimize complex models and model overfitting. A simulation study showed that these model fits can be achieved with low sample sizes (e.g. N = 25) (Gates et al., Reference Gates, Lane, Varangis, Giovanello and Guskiewicz2017). After identification of subgroups derived from GIMME, further characterization can be completed by comparing subgroups on clinical and treatment measures. Because the contemporaneous and lagged pathways were identical and highly correlated in our analysis (R 2s > 0.79), we only visualize the contemporaneous pathways for simplicity in the main report (Price, et al., Reference Price, Gates, Kraynak, Thase and Siegle2017a, Reference Price, Lane, Gates, Kraynak, Horner, Thase and Siegle2017b). Full connectivity maps are presented in online Supplementary Fig. S1b and S1d.

Subgroup characterization analysis

To examine connectivity subgroup differences in treatment outcomes, a linear mixed-effects model (lmer in the lme4 R package) was conducted on CAPS-5 scores. Connectivity subgroup, treatment group (I-PE and I-CS), and time (pre-treatment and post-treatment) were entered as fixed effects and the intercept was specified as random to account for repeated measurements per subject. We used the Satterthwaite's method to calculate the Type III sum of squares. An identical model was used to examine subgroup effects on PHDD. Post-hoc independent-samples t tests were corrected for multiple comparisons using the Benjamini and Hochberg (Reference Benjamini and Hochberg1995) false discovery rate (FDR).

Exploratory subgroup analysis

We conducted a series of exploratory analyses to further phenotype the connectivity-derived subgroups. We examined brain volume differences and additional clinical and behavioral measures. The aim of these analyses is exploratory and for generating future research questions (see Supplementary Materials).

Results

Subgroup connectivity profiles

Final model fit indices indicated excellent fit across individuals, CFI = 0.96 and SRMR = 0.03. S-GIMME identified two subgroups with common and unique network effective connectivity patterns. First, across all subjects, S-GIMME identified a series of seven contemporaneous directed connectivity pathways (Fig. 1; Table 2). These connectivity paths consisted of directed paths from the left insula to the right putamen, MFG, and the IPC, directed paths from the putamen to the pgACC and the sgACC, a path from the MFG to the pgACC, and one path from the pgACC to the pMCC. Although these paths are common across subgroups, subgroup differences in strength in these networks emerged in two out of the seven networks. Specifically, subgroup 2 had significantly lower insula to IPC connectivity (t = 3.32, p = 0.003, p FDR = 0.011) and greater pgACC to pMCC connectivity compared to subgroup 1 (t = −4.18, p = 0.0002, p FDR = 0.001; Fig. 2). Next, examination of subgroup connectivity revealed that subgroup 2 was further characterized by a pgACC to IPC path, while subgroup 1 had IPC to MFG and putamen to pMCC connectivity paths (Fig. 1b).

Fig. 1. S-GIMME identified common and unique effective connectivity. (a) S-GIMME group-level effective connectivity paths. Paths represent connections common to most subjects (red arrows). Green border on arrow represents path significantly greater for subgroup 1. Yellow border on arrow represents path significantly greater for subgroup 2. (b) Subgroup 1 was further characterized by two effective connectivity paths: putamen to pMCC, and IPC to MFG (green arrows). (c) Subgroup 2 was characterized by one additional effective connectivity path: pgACC to IPC (yellow arrows). Auto-regressive and lagged paths are not shown.

Fig. 2. Connectivity subgroup differences in group-level paths. Subgroup 1 had a greater insula to IPC path compared to subgroup 2. Subgroup 2 had a greater pgACC to pMCC path compared to subgroup 1. *ps < 0.004 (uncorrected) and psFDR < 0.02. See online version for color.

Table 2. Path estimates and Z-scores

Note. β Means, s.e.s, Z, and p values are computed as an average across subjects.

Clinical characterization of subgroups

We found a significant subgroup × therapy group × time interaction controlling for medication status [F (1,41.2) = 4.58, p = 0.038]. Within the I-CS therapy group, subgroup did not significantly predict PTSD symptom outcomes (subgroup main effect, p = 0.82, subgroup × time interaction, p = 0.63). In contrast, for patients undergoing I-PE, subgroup was associated with PTSD symptoms (subgroup × time interaction, p = 0.041). Post-hoc t tests revealed that before beginning I-PE, there were no subgroup differences (t = −0.17, p = 0.87). As shown in Fig. 3a, patients in subgroup 2 showed a greater reduction in PTSD symptoms post-treatment (t = 2.63, p = 0.022, p FDR = 0.044). Subgroups particularly moderated treatment outcomes on avoidance and hyperarousal symptoms (three-way interaction ps < 0.017; Supplementary Materials).

Fig. 3. Connectivity cluster moderated psychotherapy outcome for PTSD symptoms. (a) Clinician Administered PTSD Scale for DSM-5 (CAPS-5) total severity scores (estimated marginal means) for pre- and post-treatment assessments. Connectivity subgroup 2 showed a bigger change in PTSD symptoms than subgroup 1 in the I-PE groups: subgroup × time interaction, *p = 0.041. (b) Percentage of heavy drinking days (PHDD) for pre- and post-treatment assessments (estimated marginal means). Connectivity subgroup did not impact improvement in drinking symptoms assessed using the PHDD. See online version for color.

In light of the therapy group difference in the number of sessions attended (therapy group: t = 2.37, p = 0.02), and the qualitative difference in the number of sessions attended between connectivity subgroup (subgroup: t = −0.39, p = 0.70) and the qualitative difference in therapy drop-out status between connectivity subgroups and therapy groups (subgroup: χ2 = 1.34, p = 0.25; therapy group: χ2 = 1.97, p = 0.16), we added number of attended sessions and drop-out status as covariates into the model. The subgroup × therapy group × time interaction remained significant (p = 0.047).

To determine whether any of the seven paths that showed subgroup differences influenced treatment outcomes, we computed a multiple linear regression. We entered the seven directed connectivity path β values as the predictors of CAPS total score post-treatment with CAPS pretreatment, therapy group, connectivity cluster subgroup, and medication status as covariates. Results indicate that only the putamen to sgACC path significantly predicted CAPS post-treatment scores after controlling for CAPS at pre-treatment [β = 33.81, t = 2.43, p = 0.022 (uncorrected)].

Subgroup was unrelated to drinking behavior as measured with the PHDD (Fs < 1.09, ps > 0.30; Fig. 3b). See Supplementary Materials for the results associated with Therapy Group Type on the clinical symptoms (see also Norman et al., Reference Norman, Trim, Haller, Davis, Myers, Colvonen and Mayes2019). Connectivity subgroups identified did not significantly differ on age, sex, ethnicity, race, medication status, treatment group assignment (I-CS v. I-PE), number of treatment sessions attended, or treatment drop-out (ps > 0.25). None of the seven connectivity paths predicted drinking behavior (ps > 0.27). There were no significant differences between subgroups on MRI movement-related variables: TRs censor, DVARS, or FD (Fs < 1.21, ps > 0.23).

Exploratory results

We did find subgroup differences in volumes of a handful of regions, including the PFC and medial PFC (ps < 0.07; uncorrected). There were no differences between subgroups on diagnostic comorbidities (ps > 0.46; online Supplementary Tables S1 and S2).

Discussion

The aim of the current investigation was to identify neural connectivity subtypes of comorbid PTSD/AUD patients and to test whether the data-driven subtypes are associated with psychotherapy outcomes. We observed that the PTSD/AUD connectivity subtype 2 was associated with PTSD symptom improvement from I-PE. Furthermore, this connectivity subgroup was characterized by lower insula to IPC, higher pgACC to pMCC connectivity, and was further characterized by a pgACC to IPC connectivity path compared to subgroup 1. Both subgroups showed similar improvement from I-CS psychotherapy. Finally, we did not observe that connectivity subgroup moderated treatment outcomes for alcohol use. Our results provide novel evidence that data-driven approaches to characterize PTSD/AUD based on underlying neural connectivity patterns can detect clinically relevant subtypes.

There were five connectivity patterns that differentiated subgroup 2 from subgroup 1. Two primary regions played a particularly important role in subgroup differences. These regions included an ROI in the (1) IPC, and (2) two ROIs in the cingulate cortex: the pgACC and the pMCC. Our findings are consistent with previous work showing that PTSD and AUD are associated with aberrant activation and connectivity for each of these regions (Akiki et al., Reference Akiki, Averill and Abdallah2017; Gilpin & Weiner, Reference Gilpin and Weiner2017) and activity in these regions are associated with clinical outcomes (Malejko, Abler, Plener, & Straub, Reference Malejko, Abler, Plener and Straub2017; van Rooij, Kennis, Vink, & Geuze, Reference van Rooij, Kennis, Vink and Geuze2016; Zakiniaeiz, Scheinost, Seo, Sinha, & Constable, Reference Zakiniaeiz, Scheinost, Seo, Sinha and Constable2016). The subgroup connectivity profile differences among these regions may partially explain why subgroup 2 benefited more from I-PE than I-CS.

First, compared to subgroup 1, subgroup 2 had lower insula to IPC connectivity. This is consistent with data showing that compared to healthy controls, PTSD is associated with decreased IPC and insula connectivity (Zhang et al., Reference Zhang, Xie, Chen, Li, Guo and Chen2016). The IPC is highly connected anatomically and functionally across the brain (Wang et al., Reference Wang, Zhang, Rong, Wei, Zheng, Fox and Jiang2016; Zhang & Li, Reference Zhang and Li2014) and is a critical region implicated in attention, cognitive control, and emotion regulation (Etkin, Buchel, & Gross, Reference Etkin, Buchel and Gross2015). The insula is an integral hub in the salience network (Menon and Uddin, Reference Menon,V.2010), plays a prominent role in integrating cognitive and affectively salient information (Craig, Reference Craig2009; Gu et al., Reference Gu, Liu, Van Dam, Hof and Fan2013), and is important for interoceptive awareness (Critchley et al., Reference Critchley, Wiens, Rotshtein, Ohman and Dolan2004). Here, we show that the subgroup that benefited the most from I-PE had lower insula to IPC connectivity compared to the other subgroup, suggesting that this group may be better able to disengage attention and cognitive control from affective and interoceptive information that can interfere with goal-related behavior (Kuckertz et al., Reference Kuckertz, Amir, Boffa, Warren, Rindt, Norman and McLay2014).

Second, we found that subgroup 1 was further characterized by an IPC to MFG connectivity path. Experimental evidence indicates that higher levels of anxiety are associated with increased parietal and prefrontal activation that biases attention and working memory to potential threat (Balderston et al., Reference Balderston, Hale, Hsiung, Torrisi, Holroyd, Carver and Grillon2017; Stout, Shackman, Pedersen, Miskovich, & Larson, Reference Stout, Shackman, Pedersen, Miskovich and Larson2017). Our observation that subgroup 1 had unique IPC to MFG connectivity and did not improve as much from I-PE as subgroup 2 may suggest that subgroup 1 may be characterized by maladaptive utilization of working memory – a cognitive process important for facilitating extinction learning (Stout et al., Reference Stout, Acheson, Moore, Gur, Baker, Geyer and Risbrough2018). These results suggest that patients similar to subgroup 1 that have overactive pretreatment frontoparietal connectivity may see increased benefit from prolonged exposure therapy when an adjunctive treatment that targets the frontoparietal network is added (i.e. working memory training or repetitive transcranial magnetic stimulation) (Bomyea, Stein, & Lang, Reference Bomyea, Stein and Lang2015; Osuch et al., Reference Osuch, Benson, Luckenbaugh, Geraci, Post and McCann2009).

Third, subgroup 2 was characterized by a pgACC to IPC path. The pgACC is one of two primary regions in the anterior portion of the rostral cingulate cortex (Shackman et al., Reference Shackman, Salomons, Slagter, Fox, Winter and Davidson2011; Vogt, Nimchinsky, Vogt, & Hof, Reference Vogt, Nimchinsky, Vogt and Hof1995). The pgACC is well-positioned anatomically and functionally to exert influence on the flexible deployment of cognitive and affective control (Rolls, Reference Rolls2019; Shackman et al., Reference Shackman, Salomons, Slagter, Fox, Winter and Davidson2011; Tang et al., Reference Tang, Jbabdi, Zhu, Cottaar, Grisot, Lehman and Haber2019), is a critical node in the default mode network and internal cognitive processing (Raichle, Reference Raichle2015; Smallwood, Brown, Baird, & Schooler, Reference Smallwood, Brown, Baird and Schooler2012), and may be an important hub for dysfunction across a broad-spectrum of psychiatric disorders (Goodkind et al., Reference Goodkind, Eickhoff, Oathes, Jiang, Chang, Jones-Hagata and Etkin2015). Importantly, individuals with PTSD show dysfunctional pgACC connectivity (Kennis, Rademaker, van Rooij, Kahn, & Geuze, Reference Kennis, Rademaker, van Rooij, Kahn and Geuze2015) and decreased activation during affective interference tasks (Clausen et al., Reference Clausen, Francisco, Thelen, Bruce, Martin, McDowd and Aupperle2017). Clinically, increased IPC and pgACC connectivity is associated with improvements from cognitive behavioral therapy in depressed patients (Sambataro et al., Reference Sambataro, Doerig, Hänggi, Wolf, Brakowski, Holtforth and Spinelli2018). Furthermore, the pgACC resting-state connectivity results observed in the current report extend the structural neuroimaging findings showing the clinical importance of this region in prolonged exposure therapy for PTSD. For example, pretreatment pgACC cortical thickness and volume is associated with improved PTSD symptoms following prolonged exposure (Bryant et al., Reference Bryant, Felmingham, Whitford, Kemp, Hughes, Peduto and Williams2008). Collectively, our results in combination with the extant literature indicate that pgACC connectivity and structure play an integral role in psychological intervention improvement. Such improvement may be through the role of the pgACC in cognitive or emotion regulation of internally-based cognitions, such as rumination and traumatic memories (Palomero-Gallagher et al., Reference Palomero-Gallagher, Hoffstaedter, Mohlberg, Eickhoff, Amunts and Zilles2019).

The fourth and fifth subgroup connectivity difference involved the pMCC. The midcingulate cortex (MCC) is located on the dorsal portion of the rostral cingulate cortex and is divided into anterior and posterior (pMCC) subdivisions (Vogt et al., Reference Vogt, Nimchinsky, Vogt and Hof1995). A meta-analysis revealed that the pMCC and pgACC are consistently activated during negative affective states and pain (Shackman et al., Reference Shackman, Salomons, Slagter, Fox, Winter and Davidson2011). Furthermore, low pMCC connectivity with the amygdala prospectively predicts increased PTSD severity (Belleau et al., Reference Belleau, Ehret, Hanson, Brasel, Larson and DeRoon-Cassini2020), and PTSD patients show increased putamen activity when experiencing thermal pain, aversive images, and when experiencing an electric shock as an unconditioned stimulus during a fear learning task (Elman et al., Reference Elman, Upadhyay, Langleben, Albanese, Becerra and Borsook2018; Geuze et al., Reference Geuze, Westenberg, Jochims, de Kloet, Bohus, Vermetten and Schmahl2007; Linnman, Zeffiro, Pitman, & Milad, Reference Linnman, Zeffiro, Pitman and Milad2011). Here, we observed that subgroup 1 had weaker pgACC to pMCC and had a putamen to pMCC connectivity path, providing evidence for a circuit centered in the pMCC implicated in increased negative affect and pain that when less connected may also impair emotional processing – a necessary component of prolonged exposure therapy (Foa, Hembree, Rothbaum, & Rauch, Reference Foa, Hembree, Rothbaum and Rauch2019).

Biomarkers to aid diagnosis and treatment response are few, potentially requiring patients to cycle through multiple failed treatment trials, increasing suffering and resource expenditures. Neuroimaging-based biotyping has had recent success in predicting external measures of functioning and clinical outcomes using standard functional connectivity (Drysdale et al., Reference Drysdale, Grosenick, Downar, Dunlop, Mansouri, Meng and Liston2017; Etkin et al., Reference Etkin, Maron-Katz, Wu, Fonzo, Huemer, Vértes and O'Hara2019) and with effective connectivity measures (Li et al., Reference Li, Zhu, Jiang, Jin, Zhang, Guo and Liu2014; Price et al., Reference Price, Beltz, Woody, Cummings, Gilchrist and Siegle2020). Leveraging data-driven approaches to biotype clinical heterogeneity may better capture the complex symptom presentation of PTSD/AUD and thus be a more robust treatment response biomarker and aid in clinical decision making. Advancements in analytic strategies and methodological tools that allow for the measurement and manipulation (e.g. pharmacological or brain stimulation) of causal network relationships will be critical to enhance the clinical utility of neuroimaging (Reid et al., Reference Reid, Headley, Mill, Sanchez-Romero, Uddin, Marinazzo and Cole2019).

Limitations

There are limitations that will need to be addressed in future studies. First, our sample was restricted to a clinical sample of a mostly male Veteran PTSD/AUD patient population. We also did not include a healthy control group or separate PTSD and AUD groups to determine whether the network connectivity patterns observed here reflect disordered connectivity. Although our findings were specific to PTSD symptoms, it will be necessary to examine if the connectivity-derived subgroups and connectivity paths generalize to women and non-Veteran populations in addition to testing the specificity of effective connectivity-derived profiles in PTSD patients without AUD and in PTSD patients with other comorbidities. Second, our overall sample size is small, limiting our ability to probe how the within-subgroup connectivity paths may relate to clinical and behavioral outcomes and generalizability. However, our sample size is consistent with other recent neuroimaging investigations of psychotherapy outcomes in singly diagnosed individuals (sample size <45) (Bryant et al., Reference Bryant, Erlinger, Felmingham, Malhi, O'Donnell, Williams and Korgaonkar2020; Reggente et al., Reference Reggente, Moody, Morfini, Sheen, Rissman, O'Neill and Feusner2018; Zhutovsky et al., Reference Zhutovsky, Thomas, Olff, van Rooij, Kennis, van Wingen and Geuze2019) and S-GIMME simulation work shows robust model fits and a reduction of model overfitting even at small sample sizes as low as 25 (Gates et al., Reference Gates, Lane, Varangis, Giovanello and Guskiewicz2017). Nonetheless, further work with high-powered samples, multi-site clinical trials, and validation datasets will be required to determine the replicability and validity of the subtypes observed in our study and their treatment implications (Poldrack, Huckins, & Varoquaux, Reference Poldrack, Huckins and Varoquaux2020). Third, a limitation of effective connectivity analysis, and to multivariate analysis more broadly, is that variable selection can bias model fits when estimating directional relationships and the results are dependent on the variables entered into the model (Shi, Westerhuis, Rosén, Landberg, & Brunius, Reference Shi, Westerhuis, Rosén, Landberg and Brunius2019). Here, we used a targeted, hypothesis-driven approach for variable selection. The ROIs that we selected were based on prior research by our group that showed relevant brain regions for PTSD or AUD (Klaming et al., Reference Klaming, Harlé, Infante, Bomyea, Kim and Spadoni2019). A direction for future research will be to broaden the brain regions selected for subtyping to further identify clinically relevant neural circuits involved in benefiting from psychological interventions. Fourth, we did not find evidence of a biotype or connectivity path associated with change in alcohol use. It is unclear if AUD outcomes are separately detectable in the presence of PTSD symptomatology. The high baseline level of heavy alcohol use and low baseline variability in this measure may account for the absence of findings related to alcohol outcomes.

Conclusion

Our study provides novel evidence that the capacity to benefit from psychological interventions in patients with comorbid PTSD/AUD may involve balanced communication between key nodes from networks important for the integration of affect and cognition. Although our results do not provide direct clinical guidance, they do show that patients can be subtyped based on clinically relevant brain network relationships – providing a deeper understanding of the neural mechanisms contributing to symptom improvement that can be tested in future investigations. These results raise the possibility that recruiting clinically complex patient samples and leveraging data-driven network analyses may provide an unprecedented level of insight into the characterization of comorbid PTSD/AUD and identify brain network profiles that are associated with benefiting from standard treatments.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/S0033291721001513.

Acknowledgements

The views expressed in this article are those of the authors and do not necessarily reflect the position or policy of the Department of Veterans Affairs or the United States government.

Financial support

This work was supported by the VA Center of Excellence for Stress and Mental Health (all authors). DMS is supported by a VA Clinical Science Research and Development Career Development Award (IK2 CX001861). KMH is supported by a VA Clinical Science Research and Development Career Development Award (IK2 CX001584). SBN is supported by a VA Clinical Science Research and Development Merit Grant (I01 RX002412). ANS is supported by a VA Clinical Science Research Merit Award (I01 CX001542). ADS is supported by a VA Clinical Science Research and Development Merit Grant (I01 CX001762). The current work was supported by a VA Clinical Science Research and Development Career Development Award (IK2 CX000864) (ADS), with parent study support from a VA Clinical Science Research Merit Award (I01 CX000756) (SBN).

Conflict of interest

None.

Ethical standards

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.

References

Abi-Dargham, A., & Horga, G. (2016). The search for imaging biomarkers in psychiatric disorders. Nature Medicine, 22(11), 12481255. doi: 10.1038/nm.4190.CrossRefGoogle ScholarPubMed
Akiki, T. J., Averill, C. L., & Abdallah, C. G. (2017). A network-based neurobiological model of PTSD: Evidence from structural and functional neuroimaging studies. Current Psychiatry Reports, 19(11), 81. doi: 10.1007/s11920-017-0840-4.CrossRefGoogle ScholarPubMed
Balderston, N. L., Hale, E., Hsiung, A., Torrisi, S., Holroyd, T., Carver, F. W., … Grillon, C. (2017). Threat of shock increases excitability and connectivity of the intraparietal sulcus. ELife, 6, e23608. doi: 10.7554/eLife.23608.CrossRefGoogle ScholarPubMed
Belleau, E. L., Ehret, L. E., Hanson, J. L., Brasel, K. J., Larson, C. L., & DeRoon-Cassini, T. A. (2020). Amygdala functional connectivity in the acute aftermath of trauma prospectively predicts severity of posttraumatic stress symptoms. Neurobiology of Stress, 12, 100217. doi: 10.1016/j.ynstr.2020.100217.CrossRefGoogle ScholarPubMed
Beltz, A. M., & Gates, K. M. (2017). Network mapping with GIMME. Multivariate Behavioral Research, 52(6), 789804. doi: 10.1080/00273171.2017.1373014.CrossRefGoogle ScholarPubMed
Benjamini, Y. & Hochberg, Y. (1995). Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society. Series B (Methodological), 57, 289300. https://doi.org/10.1111/j.2517-6161.1995.tb02031.xCrossRefGoogle Scholar
Blanco, C., Wall, M. M., Lindquist, M. A., Rodríguez-Fernández, J. M., Franco, S., Wang, S., & Olfson, M. (2016). Generalizability of neuroimaging studies in 5 common psychiatric disorders based on the national epidemiologic survey on alcohol and related conditions (NESARC). The Journal of Clinical Psychiatry, 77(12), e1618e1625. doi: 10.4088/JCP.15m10264.CrossRefGoogle ScholarPubMed
Blanco, C., Xu, Y., Brady, K., Perez-Fuentes, G., Okuda, M., & Wang, S. (2013). Comorbidity of posttraumatic stress disorder with alcohol dependence among US adults: Results from National Epidemiological Survey on Alcohol and Related Conditions. Drug and Alcohol Dependence, 132(3), 630638. doi: 10.1016/j.drugalcdep.2013.04.016.CrossRefGoogle ScholarPubMed
Bomyea, J., Stein, M. B., & Lang, A. J. (2015). Interference control training for PTSD: A randomized controlled trial of a novel computer-based intervention. Journal of Anxiety Disorders, 34, 3342. doi: 10.1016/j.janxdis.2015.05.010.CrossRefGoogle ScholarPubMed
Brown, T. A. (2006). Confirmatory factor analysis for applied research. New York: The Guilford Press.Google Scholar
Bryant, R. A., Erlinger, M., Felmingham, K., Malhi, G. S., O'Donnell, M. L., Williams, L. M., … Korgaonkar, M. S. (2020). Differential neural predictors of treatment response for fear and dysphoric features of posttraumatic stress disorder. Depression and Anxiety, 37, 1026–1036. doi: 10.1002/da.23061.CrossRefGoogle ScholarPubMed
Bryant, R. A., Felmingham, K., Whitford, T. J., Kemp, A., Hughes, G., Peduto, A., & Williams, L. M. (2008). Rostral anterior cingulate volume predicts treatment response to cognitive-behavioural therapy for posttraumatic stress disorder. Journal of Psychiatry & Neuroscience, 33(2), 142146. Retrieved from https://pubmed.ncbi.nlm.nih.gov/18330460.Google ScholarPubMed
Camchong, J., Stenger, V. A., & Fein, G. (2013). Resting-state synchrony in short-term versus long-term abstinent alcoholics. Alcoholism, Clinical and Experimental Research, 37(5), 794803. doi: 10.1111/acer.12037.CrossRefGoogle ScholarPubMed
Chahal, R., Weissman, D. G., Hallquist, M. N., Robins, R. W., Hastings, P. D., & Guyer, A. E. (2020). Neural connectivity biotypes: Associations with internalizing problems throughout adolescence. Psychological Medicine, 111. Advance online publication. doi: 10.1017/S003329172000149X.Google ScholarPubMed
Cisler, J. M., Sigel, B. A., Kramer, T. L., Smitherman, S., Vanderzee, K., Pemberton, J., & Kilts, C. D. (2015). Amygdala response predicts trajectory of symptom reduction during trauma-focused cognitive-behavioral therapy among adolescent girls with PTSD. Journal of Psychiatric Research, 71, 3340. doi: 10.1016/j.jpsychires.2015.09.011.CrossRefGoogle ScholarPubMed
Clausen, A. N., Francisco, A. J., Thelen, J., Bruce, J., Martin, L. E., McDowd, J., … Aupperle, R. L. (2017). PTSD and cognitive symptoms relate to inhibition-related prefrontal activation and functional connectivity. Depression and Anxiety, 34(5), 427436. doi: 10.1002/da.22613.CrossRefGoogle ScholarPubMed
Craig, A.D. (2009). How do you feel – now? The anterior insula and human awareness. Nature Reviews Neuroscience, 10, 5970. https://doi.org/10.1038/nrn2555Google Scholar
Critchley, H. D., Wiens, S., Rotshtein, P., Ohman, A.,& Dolan, R. J. (2004). Neural systems supporting interoceptive awareness. Nature Neuroscience, 7(2), 189195. https://doi.org/10.1038/nn1176Google Scholar
Drysdale, A. T., Grosenick, L., Downar, J., Dunlop, K., Mansouri, F., Meng, Y., … Liston, C. (2017). Resting-state connectivity biomarkers define neurophysiological subtypes of depression. Nature Medicine, 23(1), 2838. doi: 10.1038/nm.4246.CrossRefGoogle ScholarPubMed
Elman, I., Upadhyay, J., Langleben, D. D., Albanese, M., Becerra, L., & Borsook, D. (2018). Reward and aversion processing in patients with post-traumatic stress disorder: Functional neuroimaging with visual and thermal stimuli. Translational Psychiatry, 8(1), 240. doi: 10.1038/s41398-018-0292-6.CrossRefGoogle ScholarPubMed
Etkin, A., Buchel, C., & Gross, J. J. (2015). The neural bases of emotion regulation. Nature Reviews. Neuroscience, 16(11), 693700. doi: 10.1038/nrn4044.CrossRefGoogle ScholarPubMed
Etkin, A., Maron-Katz, A., Wu, W., Fonzo, G. A., Huemer, J., Vértes, P. E., … O'Hara, R. (2019). Using fMRI connectivity to define a treatment-resistant form of post-traumatic stress disorder. Science Translational Medicine, 11(486), eaal3236. doi: 10.1126/scitranslmed.aal3236.CrossRefGoogle ScholarPubMed
Feczko, E., Miranda-Dominguez, O., Marr, M., Graham, A. M., Nigg, J. T., & Fair, D. A. (2019). The heterogeneity problem: Approaches to identify psychiatric subtypes. Trends in Cognitive Sciences, 23(7), 584601. doi: 10.1016/j.tics.2019.03.009.CrossRefGoogle ScholarPubMed
Foa, E., Hembree, E., Rothbaum, B. O., & Rauch, S. (2019). Prolonged exposure therapy for PTSD: Emotional processing of traumatic experiences therapist guide (2nd ed.). New York, NY: Oxford University Press. 10.1093/med: psych/9780195308501.001.0001.CrossRefGoogle Scholar
Gates, K. M., Henry, T., Steinley, D., & Fair, D. A. (2016). A Monte Carlo evaluation of weighted community detection algorithms. Frontiers in Neuroinformatics, 10, 45. doi: 10.3389/fninf.2016.00045.CrossRefGoogle ScholarPubMed
Gates, K. M., Lane, S. T., Varangis, E., Giovanello, K., & Guskiewicz, K. (2017). Unsupervised classification during time-series model building. Multivariate Behavioral Research, 52(2), 129148. doi: 10.1080/00273171.2016.1256187.CrossRefGoogle ScholarPubMed
Gates, K. M., Molenaar, P. C. M., Iyer, S. P., Nigg, J. T., & Fair, D. A. (2014). Organizing heterogeneous samples using community detection of GIMME-derived resting state functional networks. PLoS ONE, 9(3), e91322. doi: 10.1371/journal.pone.0091322.CrossRefGoogle ScholarPubMed
Gates, K. M., & Molenaar, P. C. M. (2012). Group search algorithm recovers effective connectivity maps for individuals in homogeneous and heterogeneous samples. NeuroImage, 63(1), 310319. doi: 10.1016/j.neuroimage.2012.06.026.CrossRefGoogle ScholarPubMed
Geuze, E., Westenberg, H. G. M., Jochims, A., de Kloet, C. S., Bohus, M., Vermetten, E., & Schmahl, C. (2007). Altered pain processing in Veterans with posttraumatic stress disorder. Archives of General Psychiatry, 64(1), 7685. doi: 10.1001/archpsyc.64.1.76.CrossRefGoogle ScholarPubMed
Gilpin, N. W., & Weiner, J. L. (2017). Neurobiology of comorbid post-traumatic stress disorder and alcohol-use disorder. Genes, Brain, and Behavior, 16(1), 1543. doi: 10.1111/gbb.12349.CrossRefGoogle ScholarPubMed
Goodkind, M., Eickhoff, S. B., Oathes, D. J., Jiang, Y., Chang, A., Jones-Hagata, L. B., … Etkin, A. (2015). Identification of a common neurobiological substrate for mental illness. JAMA Psychiatry, 72(4), 305315. doi: 10.1001/jamapsychiatry.2014.2206.CrossRefGoogle ScholarPubMed
Gu, X., Liu, X., Van Dam, N. T., Hof, P. R.& Fan, J. (2013). Cognition-emotion integration in the anterior insular cortex. Cerebral Cortex, 23, 20–27. https://doi.org/10.1093/cercor/bhr367Google Scholar
Huys, Q. J. M., Maia, T. V., & Frank, M. J. (2016). Computational psychiatry as a bridge from neuroscience to clinical applications. Nature Neuroscience, 19(3), 404413. doi: 10.1038/nn.4238.CrossRefGoogle ScholarPubMed
Jin, C., Jia, H., Lanka, P., Rangaprakash, D., Li, L., Liu, T., … Deshpande, G. (2017). Dynamic brain connectivity is a better predictor of PTSD than static connectivity. Human Brain Mapping, 38(9), 44794496. doi: 10.1002/hbm.23676.CrossRefGoogle ScholarPubMed
Kennis, M., Rademaker, A. R., van Rooij, S. J. H., Kahn, R. S., & Geuze, E. (2015). Resting state functional connectivity of the anterior cingulate cortex in veterans with and without post-traumatic stress disorder. Human Brain Mapping, 36(1), 99109. doi: 10.1002/hbm.22615.CrossRefGoogle ScholarPubMed
Klaming, R., Harlé, K. M., Infante, M. A., Bomyea, J., Kim, C., & Spadoni, A. D. (2019). Shared gray matter reductions across alcohol use disorder and posttraumatic stress disorder in the anterior cingulate cortex: A dual meta-analysis. Neurobiology of Stress, 10, 100132. doi: 10.1016/j.ynstr.2018.09.009.CrossRefGoogle ScholarPubMed
Koch, S. B. J., van Zuiden, M., Nawijn, L., Frijling, J. L., Veltman, D. J., & Olff, M. (2016). Aberrant resting-state brain activity in posttraumatic stress disorder: A meta-analysis and systematic review. Depression and Anxiety, 33(7), 592605. doi: 10.1002/da.22478.CrossRefGoogle ScholarPubMed
Kohno, M., Dennis, L. E., McCready, H., & Hoffman, W. F. (2017). Executive control and striatal resting-state network interact with risk factors to influence treatment outcomes in alcohol-use disorder. Frontiers in Psychiatry, 8, 182. doi: 10.3389/fpsyt.2017.00182.CrossRefGoogle ScholarPubMed
Krystal, J. H., Davis, L. L., Neylan, T. C., Raskind, M. A., Schnurr, P. P., Stein, M. B., … Huang, G. D. (2017). It is time to address the crisis in the pharmacotherapy of posttraumatic stress disorder: A consensus statement of the PTSD psychopharmacology working group. Biological Psychiatry, 82(7), e51e59. doi: 10.1016/j.biopsych.2017.03.007.CrossRefGoogle ScholarPubMed
Kuckertz, J. M., Amir, N., Boffa, J. W., Warren, C. K., Rindt, S. E. M., Norman, S., … McLay, R. (2014). The effectiveness of an attention bias modification program as an adjunctive treatment for post-traumatic stress disorder. Behaviour Research and Therapy, 63, 2535. doi: 10.1016/j.brat.2014.09.002.CrossRefGoogle ScholarPubMed
Lazarus, G., Sened, H., & Rafaeli, E. (2020). Subjectifying the personality state: Theoretical underpinnings and an empirical example. European Journal of Personality, 34(6), 1017–1036. doi: 10.1002/per.2278.CrossRefGoogle Scholar
Li, X., Zhu, D., Jiang, X., Jin, C., Zhang, X., Guo, L., … Liu, T. (2014). Dynamic functional connectomics signatures for characterization and differentiation of PTSD patients. Human Brain Mapping, 35(4), 17611778. doi: 10.1002/hbm.22290.CrossRefGoogle ScholarPubMed
Linnman, C., Zeffiro, T. A., Pitman, R. K., & Milad, M. R. (2011). An fMRI study of unconditioned responses in post-traumatic stress disorder. Biology of Mood & Anxiety Disorders, 1(1), 8. doi: 10.1186/2045-5380-1-8.CrossRefGoogle ScholarPubMed
Malejko, K., Abler, B., Plener, P. L., & Straub, J. (2017). Neural correlates of psychotherapeutic treatment of post-traumatic stress disorder: A systematic literature review. Frontiers in Psychiatry, 8, 85. Retrieved from https://www.frontiersin.org/article/10.3389/fpsyt.2017.00085.CrossRefGoogle ScholarPubMed
Maron-Katz, A., Zhang, Y., Narayan, M., Wu, W., Toll, R. T., Naparstek, S., … Etkin, A. (2019). Individual patterns of abnormality in resting-state functional connectivity reveal two data-driven PTSD subgroups. American Journal of Psychiatry, 177(3), 244253. doi: 10.1176/appi.ajp.2019.19010060.CrossRefGoogle ScholarPubMed
McCormick, E. M., & Telzer, E. H. (2018). Contributions of default mode network stability and deactivation to adolescent task engagement. Scientific Reports, 8(1), 18049. doi: 10.1038/s41598-018-36269-4.CrossRefGoogle ScholarPubMed
Menon,V., & Uddin, L.Q.. (2010). Saliency, switching, attention and control: a network model of insula function. Brain Structure & Function, 214, 655–667. https://doi.org/10.1007/s00429-010-0262-0.Google Scholar
Milham, M. P., Craddock, R. C., & Klein, A. (2017). Clinically useful brain imaging for neuropsychiatry: How can we get there? Depression and Anxiety, 34(7), 578587. doi: 10.1002/da.22627.CrossRefGoogle Scholar
Norman, S. B., Haller, M., Hamblen, J. L., Southwick, S. M., & Pietrzak, R. H. (2018). The burden of co-occurring alcohol use disorder and PTSD in U.S. Military veterans: Comorbidities, functioning, and suicidality. Psychology of Addictive Behaviors, 32(2), 224229. doi: 10.1037/adb0000348.CrossRefGoogle ScholarPubMed
Norman, S. B., Trim, R., Haller, M., Davis, B. C., Myers, U. S., Colvonen, P. J., … Mayes, T. (2019). Efficacy of integrated exposure therapy vs integrated coping skills therapy for comorbid posttraumatic stress disorder and alcohol use disorder: A randomized clinical trial. JAMA Psychiatry, 76(8), 791799. doi: 10.1001/jamapsychiatry.2019.0638.CrossRefGoogle ScholarPubMed
Orman, G. K., Labatut, V., & Cherifi, H. (2012). Comparative evaluation of community detection algorithms: A topological approach. Journal of Statistical Mechanics: Theory and Experiment, 2012(08), P08001. doi: 10.1088/1742-5468/2012/08/p08001.CrossRefGoogle Scholar
Osuch, E. A., Benson, B. E., Luckenbaugh, D. A., Geraci, M., Post, R. M., & McCann, U. (2009). Repetitive TMS combined with exposure therapy for PTSD: A preliminary study. Journal of Anxiety Disorders, 23(1), 5459. doi: 10.1016/j.janxdis.2008.03.015.CrossRefGoogle ScholarPubMed
Palomero-Gallagher, N., Hoffstaedter, F., Mohlberg, H., Eickhoff, S. B., Amunts, K., & Zilles, K. (2019). Human pregenual anterior cingulate cortex: Structural, functional, and connectional heterogeneity. Cerebral Cortex, 29(6), 25522574. doi: 10.1093/cercor/bhy124.CrossRefGoogle ScholarPubMed
Poldrack, R. A., Huckins, G., & Varoquaux, G. (2020). Establishment of best practices for evidence for prediction: A review. JAMA Psychiatry, 77(5), 534540. doi: 10.1001/jamapsychiatry.2019.3671.CrossRefGoogle ScholarPubMed
Price, R. B., Beltz, A. M., Woody, M. L., Cummings, L., Gilchrist, D., & Siegle, G. J. (2020). Neural connectivity subtypes predict discrete attentional-bias profiles among heterogeneous anxiety patients. Clinical Psychological Science, 8(3), 491505. doi: 10.1177/2167702620906149.CrossRefGoogle ScholarPubMed
Price, R. B., Gates, K., Kraynak, T. E., Thase, M. E., & Siegle, G. J. (2017a). Data-driven subgroups in depression derived from directed functional connectivity paths at rest. Neuropsychopharmacology, 42(13), 26232632. doi: 10.1038/npp.2017.97.CrossRefGoogle ScholarPubMed
Price, R. B., Lane, S., Gates, K., Kraynak, T. E., Horner, M. S., Thase, M. E., & Siegle, G. J. (2017b). Parsing heterogeneity in the brain connectivity of depressed and healthy adults during positive mood. Biological Psychiatry, 81(4), 347357. doi: 10.1016/j.biopsych.2016.06.023.CrossRefGoogle ScholarPubMed
Raichle, M. E. (2015). The brain's default mode network. Annual Review of Neuroscience, 38(1), 433447. doi: 10.1146/annurev-neuro-071013-014030.CrossRefGoogle ScholarPubMed
Reggente, N., Moody, T. D., Morfini, F., Sheen, C., Rissman, J., O'Neill, J., & Feusner, J. D. (2018). Multivariate resting-state functional connectivity predicts response to cognitive behavioral therapy in obsessive–compulsive disorder. Proceedings of the National Academy of Sciences, 115(9), 2222 LP2227. doi: 10.1073/pnas.1716686115.CrossRefGoogle ScholarPubMed
Reid, A. T., Headley, D. B., Mill, R. D., Sanchez-Romero, R., Uddin, L. Q., Marinazzo, D., … Cole, M. W. (2019). Advancing functional connectivity research from association to causation. Nature Neuroscience, 22(11), 17511760. doi: 10.1038/s41593-019-0510-4.CrossRefGoogle ScholarPubMed
Roberts, N. P., Roberts, P. A., Jones, N., & Bisson, J. I. (2015). Psychological interventions for post-traumatic stress disorder and comorbid substance use disorder: A systematic review and meta-analysis. Clinical Psychology Review, 38, 2538. doi: 10.1016/j.cpr.2015.02.007.CrossRefGoogle ScholarPubMed
Robinaugh, D. J., Hoekstra, R. H. A., Toner, E. R., & Borsboom, D. (2020). The network approach to psychopathology: A review of the literature 2008–2018 and an agenda for future research. Psychological Medicine, 50(3), 353366. doi: 10.1017/S0033291719003404.CrossRefGoogle Scholar
Rolls, E. T. (2019). The cingulate cortex and limbic systems for emotion, action, and memory. Brain Structure and Function, 224(9), 30013018. doi: 10.1007/s00429-019-01945-2.CrossRefGoogle ScholarPubMed
Ruglass, L. M., Lopez-Castro, T., Papini, S., Killeen, T., Back, S. E., & Hien, D. A. (2017). Concurrent treatment with prolonged exposure for co-occurring full or subthreshold posttraumatic stress disorder and substance use disorders: A randomized clinical trial. Psychotherapy and Psychosomatics, 86(3), 150161. doi: 10.1159/000462977.CrossRefGoogle ScholarPubMed
Sambataro, F., Doerig, N., Hänggi, J., Wolf, R. C., Brakowski, J., Holtforth, M. G., … Spinelli, S. (2018). Anterior cingulate volume predicts response to psychotherapy and functional connectivity with the inferior parietal cortex in major depressive disorder. European Neuropsychopharmacology, 28(1), 138148. doi: 10.1016/j.euroneuro.2017.11.008.CrossRefGoogle ScholarPubMed
Shackman, A. J., Salomons, T. V., Slagter, H. A., Fox, A. S., Winter, J. J., & Davidson, R. J. (2011). The integration of negative affect, pain and cognitive control in the cingulate cortex. Nature Reviews Neuroscience, 12(3), 154167. doi: 10.1038/nrn2994.CrossRefGoogle ScholarPubMed
Sheynin, J., Duval, E. R., King, A. P., Angstadt, M., Phan, K. L., Simon, N. M., … Liberzon, I. (2020). Associations between resting-state functional connectivity and treatment response in a randomized clinical trial for posttraumatic stress disorder. Depression and Anxiety, 37, 1037–1046. doi: 10.1002/da.23075.CrossRefGoogle Scholar
Sheynin, J., & Liberzon, I. (2017). Circuit dysregulation and circuit-based treatments in posttraumatic stress disorder. Neuroscience Letters, 649, 133138. doi: 10.1016/j.neulet.2016.11.014.CrossRefGoogle ScholarPubMed
Shi, L., Westerhuis, J. A., Rosén, J., Landberg, R., & Brunius, C. (2019). Variable selection and validation in multivariate modelling. Bioinformatics (Oxford, England), 35(6), 972980. doi: 10.1093/bioinformatics/bty710.Google ScholarPubMed
Simmons, A. N., Norman, S. B., Spadoni, A. D., & Strigo, I. A. (2013). Neurosubstrates of remission following prolonged exposure therapy in veterans with posttraumatic stress disorder. Psychotherapy and Psychosomatics, 82(6), 382389. doi: 10.1159/000348867.CrossRefGoogle ScholarPubMed
Smallwood, J., Brown, K., Baird, B., & Schooler, J. W. (2012). Cooperation between the default mode network and the frontal-parietal network in the production of an internal train of thought. Brain Research, 1428, 6070. doi: 10.1016/j.brainres.2011.03.072.CrossRefGoogle ScholarPubMed
Sobell, L. C., & Sobell, M. B. (1992). Timeline Follow-Back: A Technique for Assessing Self-Reported Alcohol Consumption. In Litten, R. Z., & Allen, J. P. (Eds.), Measuring Alcohol Consumption (pp. 4172). Totowa, NJ: Humana Press. 10.1007/978-1-4612-0357-5_3.CrossRefGoogle Scholar
Steenkamp, M. M. (2016). True evidence-based care for posttraumatic stress disorder in military personnel and veterans. JAMA Psychiatry, 73(5), 431432. doi: 10.1001/jamapsychiatry.2015.2879.CrossRefGoogle ScholarPubMed
Stout, D. M., Acheson, D. T., Moore, T. M., Gur, R. C., Baker, D. G., Geyer, M. A., & Risbrough, V. B. (2018). Individual variation in working memory is associated with fear extinction performance. Behaviour Research and Therapy, 102, 5259. doi: 10.1016/j.brat.2018.01.002.CrossRefGoogle ScholarPubMed
Stout, D. M., Shackman, A. J., Pedersen, W. S., Miskovich, T. A., & Larson, C. L. (2017). Neural circuitry governing anxious individuals’ mis-allocation of working memory to threat. Scientific Reports, 7(1), 8742. doi: 10.1038/s41598-017-08443-7.CrossRefGoogle ScholarPubMed
Tang, W., Jbabdi, S., Zhu, Z., Cottaar, M., Grisot, G., Lehman, J. F., … Haber, S. N. (2019). A connectional hub in the rostral anterior cingulate cortex links areas of emotion and cognitive control. ELife, 8, e43761. doi: 10.7554/eLife.43761.CrossRefGoogle ScholarPubMed
van Rooij, S. J. H., Kennis, M., Vink, M., & Geuze, E. (2016). Predicting treatment outcome in PTSD: A longitudinal functional MRI study on trauma-unrelated emotional processing. Neuropsychopharmacology, 41(4), 11561165. doi: 10.1038/npp.2015.257.CrossRefGoogle Scholar
Vogt, B. A., Nimchinsky, E. A., Vogt, L. J., & Hof, P. R. (1995). Human cingulate cortex: Surface features, flat maps, and cytoarchitecture. The Journal of Comparative Neurology, 359(3), 490506. doi: 10.1002/cne.903590310.CrossRefGoogle ScholarPubMed
Wang, J., Zhang, J., Rong, M., Wei, X., Zheng, D., Fox, P. T., … Jiang, T. (2016). Functional topography of the right inferior parietal lobule structured by anatomical connectivity profiles. Human Brain Mapping, 37(12), 43164332. doi: 10.1002/hbm.23311.CrossRefGoogle ScholarPubMed
Weathers, F. W., Bovin, M. J., Lee, D. J., Sloan, D. M., Schnurr, P. P., Kaloupek, D. G., … Marx, B. P. (2018). The Clinician-Administered PTSD Scale for DSM-5 (CAPS-5): Development and initial psychometric evaluation in military veterans. Psychological Assessment, 30(3), 383–395. doi: 10.1037/pas0000486.CrossRefGoogle ScholarPubMed
Woo, C.-W., Chang, L. J., Lindquist, M. A., & Wager, T. D. (2017). Building better biomarkers: Brain models in translational neuroimaging. Nature Neuroscience, 20(3), 365377. doi: 10.1038/nn.4478.CrossRefGoogle ScholarPubMed
Yamada, T., Hashimoto, R.-I., Yahata, N., Ichikawa, N., Yoshihara, Y., Okamoto, Y., … Kawato, M. (2017). Resting-state functional connectivity-based biomarkers and functional MRI-based neurofeedback for psychiatric disorders: A challenge for developing theranostic biomarkers. The International Journal of Neuropsychopharmacology, 20(10), 769781. doi: 10.1093/ijnp/pyx059.CrossRefGoogle ScholarPubMed
Zakiniaeiz, Y., Scheinost, D., Seo, D., Sinha, R., & Constable, R. T. (2016). Cingulate cortex functional connectivity predicts future relapse in alcohol dependent individuals. NeuroImage: Clinical, 13, 181187. doi: 10.1016/j.nicl.2016.10.019.CrossRefGoogle ScholarPubMed
Zhang, S., & Li, C.-S. R. (2014). Functional clustering of the human inferior parietal lobule by whole-brain connectivity mapping of resting-state functional magnetic resonance imaging signals. Brain Connectivity, 4(1), 5369. doi: 10.1089/brain.2013.0191.Google ScholarPubMed
Zhang, Y., Xie, B., Chen, H., Li, M., Guo, X., & Chen, H. (2016). Disrupted resting-state insular subregions functional connectivity in post-traumatic stress disorder. Medicine, 95(27), e4083. doi: 10.1097/MD.0000000000004083.CrossRefGoogle ScholarPubMed
Zhutovsky, P., Thomas, R. M., Olff, M., van Rooij, S. J. H., Kennis, M., van Wingen, G. A., & Geuze, E. (2019). Individual prediction of psychotherapy outcome in posttraumatic stress disorder using neuroimaging data. Translational Psychiatry, 9(1), 326. doi: 10.1038/s41398-019-0663-7.CrossRefGoogle ScholarPubMed
Figure 0

Table 1. Descriptive statistics and sample characteristics

Figure 1

Fig. 1. S-GIMME identified common and unique effective connectivity. (a) S-GIMME group-level effective connectivity paths. Paths represent connections common to most subjects (red arrows). Green border on arrow represents path significantly greater for subgroup 1. Yellow border on arrow represents path significantly greater for subgroup 2. (b) Subgroup 1 was further characterized by two effective connectivity paths: putamen to pMCC, and IPC to MFG (green arrows). (c) Subgroup 2 was characterized by one additional effective connectivity path: pgACC to IPC (yellow arrows). Auto-regressive and lagged paths are not shown.

Figure 2

Fig. 2. Connectivity subgroup differences in group-level paths. Subgroup 1 had a greater insula to IPC path compared to subgroup 2. Subgroup 2 had a greater pgACC to pMCC path compared to subgroup 1. *ps < 0.004 (uncorrected) and psFDR < 0.02. See online version for color.

Figure 3

Table 2. Path estimates and Z-scores

Figure 4

Fig. 3. Connectivity cluster moderated psychotherapy outcome for PTSD symptoms. (a) Clinician Administered PTSD Scale for DSM-5 (CAPS-5) total severity scores (estimated marginal means) for pre- and post-treatment assessments. Connectivity subgroup 2 showed a bigger change in PTSD symptoms than subgroup 1 in the I-PE groups: subgroup × time interaction, *p = 0.041. (b) Percentage of heavy drinking days (PHDD) for pre- and post-treatment assessments (estimated marginal means). Connectivity subgroup did not impact improvement in drinking symptoms assessed using the PHDD. See online version for color.

Supplementary material: File

Stout et al. supplementary material

Stout et al. supplementary material

Download Stout et al. supplementary material(File)
File 5.6 MB