Brain-Based Biotypes of Psychiatric Vulnerability in the Acute Aftermath of Trauma
Abstract
Objective:
Major negative life events, such as trauma exposure, can play a key role in igniting or exacerbating psychopathology. However, few disorders are diagnosed with respect to precipitating events, and the role of these events in the unfolding of new psychopathology is not well understood. The authors conducted a multisite transdiagnostic longitudinal study of trauma exposure and related mental health outcomes to identify neurobiological predictors of risk, resilience, and different symptom presentations.
Methods:
A total of 146 participants (discovery cohort: N=69; internal replication cohort: N=77) were recruited from emergency departments within 72 hours of a trauma and followed for the next 6 months with a survey, MRI, and physiological assessments.
Results:
Task-based functional MRI 2 weeks after a motor vehicle collision identified four clusters of individuals based on profiles of neural activity reflecting threat reactivity, reward reactivity, and inhibitory engagement. Three clusters were replicated in an independent sample with a variety of trauma types. The clusters showed different longitudinal patterns of posttrauma symptoms.
Conclusions:
These findings provide a novel characterization of heterogeneous stress responses shortly after trauma exposure, identifying potential neuroimaging-based biotypes of trauma resilience and psychopathology.
The diathesis-stress model of psychopathology has remained one of the most well-supported theories addressing the causes of mental disorders. In combination with predisposing factors, antecedent stressors increase risk for the onset and recurrence of depression ( 1), schizophrenia ( 2), insomnia ( 3), and posttraumatic stress disorder (PTSD) ( 4). However, stress and its severity or chronicity alone cannot account for the wide variety of different types of mental health outcomes that can follow major stressful life events ( 5), theoretically driven by existing individual differences ( 6). These variations and their biological bases are not well captured by existing definitions of psychiatric disorders. In the present study, our objective was to discover brain-based profiles to map heterogeneity following a stressor in a nationwide longitudinal study of trauma exposure and subsequent mental health outcomes, the Advancing Understanding of Recovery After Trauma (AURORA) study ( 7).
Neuroimaging is an attractive tool for mapping symptoms to biology. Previous efforts to account for heterogeneity have often explored brain-wide patterns of activation or connectivity to identify “biotypes,” subtypes of a particular form of psychopathology that differ in their neurophysiological features ( 8– 10). The identification of such subgroups may in turn improve our understanding of variance in outcomes and response to treatment. However, previous work has defined the neuroimaging features of interest on the basis of their association with either a specific symptom type or the response to treatment. These backward inferences constrain the solution to features that already have high relevance to a diagnostic category, potentially excluding features that contribute to atypical symptom profiles and raising concerns related to overfitting ( 11). Taking a complementary approach, we constructed a forward inference model, examining neuroimaging profiles in the acute posttrauma period and then investigating their association with the emergence of later symptoms. The goal was to identify posttrauma biotypes with relevance to overall stress vulnerability and resilience but not specific to a particular diagnosis or symptom.
Neural models of stress vulnerability involve hyperreactivity of regions involved in threat detection and the fear response, such as the amygdala, insula, and dorsal anterior cingulate cortex (dACC) ( 12– 15). In addition, both chronic depression and PTSD emerging following a stressful event appear to be preceded by low reward reactivity in affective-evaluative regions, including the nucleus accumbens (NAcc), amygdala, and orbitofrontal cortex (OFC) ( 16– 20). Finally, reduced pretrauma or early posttrauma engagement of regions involved in inhibition, including the ventromedial prefrontal cortex (vmPFC) and hippocampus, is also predictive of greater subsequent PTSD and lower resilience ( 12, 21– 24). Therefore, an early posttrauma profile of co-occurring high threat reactivity, low reward reactivity, and low inhibition would likely be predictive of later chronic symptoms of PTSD and depression. However, it is not yet clear whether these features co-occur within particular individuals or groups. Building a brain-based model of individual differences in the response to major stressors is critical for efforts to construct effective intervention and prevention strategies for stress-related psychiatric disorders.
Here, we collected functional MRI (fMRI) scans in a regionally diverse cohort of civilian trauma survivors 2 weeks posttrauma ( 7). Participants reported on symptoms of psychopathology through mobile surveys over the first 6 months posttrauma. fMRI-based phenotypes used in the biotyping analysis were motivated by previous longitudinal studies of stress exposure and included all brain regions previously linked with vulnerability to poststress psychopathology within the domains of threat responsivity, reward responsivity, and inhibition/impulsivity as described above. Participants engaged in fMRI tasks that were simple in their design and interpretation and that have been widely used to probe threat ( 12), reward ( 25), and response inhibition ( 26). Multivariate profiles of regional activation were entered into a hierarchical clustering analysis to identify brain-based groupings of individuals in the early posttrauma period, indicative of distinct biotypes. We predicted that fMRI-based clusters would be associated with different patterns of subsequent posttrauma symptoms across PTSD, dissociation, anxiety, depression, and impulsivity. Finally, to better understand whether the clusters overlap with widely known biomarkers of chronic posttraumatic pathology, such as deficits in fear inhibition ( 27) and extinction ( 28), we tested whether the clusters differed in these features in a fear-potentiated startle paradigm conducted on the same day as the fMRI scan.
METHODS
Participants
Participants were recruited from emergency departments as part of a multisite longitudinal study of adverse neuropsychiatric sequelae of trauma ( 7). Twenty-two emergency departments within the Northeast, Southern, Mid-Atlantic, and Midwest regions of the United States enrolled patients within 72 hours of trauma exposure. All participants were ages 18–75, able to speak and read English, oriented to time and place, and physically able to use a smartphone, and they had possessed a smartphone for more than 1 year. Potential participants were excluded if they had a solid organ injury greater than grade 1 or a significant hemorrhage, required a chest tube or general anesthesia, or were likely to be admitted for >72 hours. MRI scans and psychophysiology data were collected a mean of 18 days (SD=6) later at a laboratory visit at McLean Hospital (Belmont, Mass.), Emory University (Atlanta), Temple University (Philadelphia), or Wayne State University (Detroit), which were each located in proximity to multiple enrolling sites. Written informed consent was obtained as approved by each site’s institutional review board.
Data collection for the AURORA study is ongoing. The discovery cohort included an initial sample that was restricted to motor vehicle-related traumas for participants with at least 8 weeks of follow-up data by March 2019 (94 patients completed the MRI visit). Data for the replication cohort were separate from the discovery cohort because these data had not yet been released by the time of the initial analysis. A second freeze and release of the survey data was broadened to include all trauma types with at least 8 weeks of follow-up data by mid-October 2019. Unique participants in this second freeze made up the replication cohort (additional participants, N=108). After quality control, 69 participants in the discovery cohort and 77 in the replication cohort were retained for analyses. The study participants’ demographic characteristics are presented in Table 1.
Cohort | ||||||
---|---|---|---|---|---|---|
Discovery (N=69) | Replication (N=77) | Group Comparison | ||||
Characteristic | N | % | N | % | χ 2 | p |
Site | 1.81 | 0.61 | ||||
1 | 5 | 7.2 | 2 | 2.6 | ||
2 | 28 | 40.7 | 32 | 41.6 | ||
3 | 19 | 27.5 | 24 | 31.2 | ||
4 | 17 | 24.6 | 19 | 24.7 | ||
Female | 51 | 73.9 | 48 | 62.3 | 2.23 | 0.14 |
Race/ethnicity | 6.62 | 0.09 | ||||
Hispanic/Latino | 11 | 15.9 | 12 | 15.6 | ||
White | 20 | 29.0 | 34 | 44.2 | ||
Black/African American | 37 | 53.6 | 27 | 35.1 | ||
Asian, Hawaiian, American Indian, or other | 1 | 1.4 | 4 | 5.2 | ||
Education | 20.02 | 0.05 | ||||
Less than high school | 3 | 4.3 | 3 | 3.9 | ||
High school diploma/GED | 26 | 37.7 | 16 | 20.8 | ||
Some college/associate’s degree | 29 | 42.0 | 29 | 37.6 | ||
Bachelor’s degree | 8 | 11.6 | 19 | 24.7 | ||
Master’s degree | 3 | 4.3 | 8 | 10.4 | ||
Doctorate | 0 | 0.0 | 2 | 2.6 | ||
Employment status | 2.15 | 0.71 | ||||
Employed | 48 | 69.6 | 48 | 62.3 | ||
Retired | 1 | 1.4 | 3 | 3.9 | ||
Homemaker | 1 | 1.4 | 1 | 1.3 | ||
Student | 2 | 2.9 | 4 | 5.2 | ||
Unemployed, disabled, or other | 9 | 13.0 | 14 | 18.2 | ||
Data missing | 8 | 11.6 | 7 | 9.1 | ||
Yearly family income | 4.49 | 0.48 | ||||
<$19,000 | 17 | 24.6 | 17 | 22.1 | ||
$19,001–$35,000 | 18 | 26.1 | 19 | 24.7 | ||
$35,001–$50,000 | 10 | 14.5 | 13 | 16.9 | ||
$50,001–$75,000 | 4 | 5.8 | 9 | 11.7 | ||
$75,001–$100,000 | 3 | 4.3 | 7 | 9.1 | ||
>$100,000 | 9 | 13.0 | 5 | 6.5 | ||
Data missing | 8 | 11.6 | 7 | 9.1 | ||
Trauma type | 42.82 | 9.51 × 10 –7 | ||||
Motor vehicle collision | 69 | 100.0 | 41 | 53.2 | ||
Physical assault | 0 | 0.0 | 15 | 19.5 | ||
Sexual assault | 0 | 0.0 | 2 | 2.6 | ||
Incident causing traumatic stress exposure to many people | 0 | 0.0 | 1 | 1.3 | ||
Nonmotorized collision (e.g., bicycle, skateboard) | 0 | 0.0 | 6 | 7.8 | ||
Fall <10 feet | 0 | 0.0 | 4 | 5.2 | ||
Burn | 0 | 0.0 | 1 | 1.3 | ||
Animal-related | 0 | 0.0 | 3 | 3.9 | ||
Other | 0 | 0.0 | 4 | 5.2 | ||
Mean | SD | Mean | SD | t | p | |
Age (years) | 33.6 | 12.1 | 35.3 | 13.7 | –0.80 | 0.42 |
Injury severity scale | 2.5 | 2.4 | 2.3 | 1.8 | 0.80 | 0.42 |
Patient-rated chance of dying | 5.6 | 3.6 | 4.8 | 3.3 | 1.48 | 0.14 |
Demographic and trauma-related characteristics of trauma survivors in the Advancing Understanding of Recovery After Trauma study
Demographic Variables and Psychiatric Assessment
Trauma severity was measured using an injury severity score, as well as participants’ subjective ratings of their chances of dying. Assessments of pretrauma risk factors included a general physical health status assessment, childhood maltreatment assessment, and demographic variables. Assessments of posttrauma outcomes, including PTSD symptoms, depression symptoms, dissociative symptoms, anxiety symptoms, and impulsivity, were assessed for the pretrauma period (queried in the emergency department) and at 2 weeks (days 7–21), 8 weeks (days 46–67), 3 months (days 77–104), and 6 months (days 168–195) posttrauma. Measures and scoring details are summarized in the online supplement.
MRI
Acquisition.
Brain imaging data were acquired on four separate Siemens 3-T MRI scanners using the two-dimensional echo-planar blood-oxygen-level-dependent sequence for functional scans and a magnetization-prepared rapid acquisition gradient echo T 1-weighted image for structural scans. Site-specific sequence parameters are presented in Table S5 in the online supplement.
fMRI tasks.
The three fMRI tasks ( Figure 1) included a threat task designed to probe reactivity to social threat cues ( 12), an inhibition task, which was a modified version of Liebenluft’s stop-signal task ( 26), and a reward task, which was a short version of Delgado’s monetary reward task ( 25).
Preprocessing and analysis.
Full preprocessing information is reported in the online supplement. Functional images were preprocessed with fMRIPrep, version 1.2.2 ( 29). Echo-planar imaging scans were coregistered to the T 1-weighted images, then spatially realigned, slice-time corrected, and normalized to the 2009 ICBM-152 template. Volume-wise motion and other sources of artifact were corrected using ICA-AROMA ( 30). To handle cases in which motion was likely too high for effective independent component analysis correction, we also implemented an overall motion threshold for any run with >15% of volumes showing ≥1-mm framewise displacement. Images were then smoothed with a 6-mm kernel. Site-by-site quality metrics are plotted in Figure S3 in the online supplement.
The final sample was restricted to participants with good-quality data across all three fMRI tasks (threat, inhibition, and reward). We did not interpolate any data point because the goal of the clustering analysis was to identify existing patterns of activation across the three tasks. Participants were excluded for falx calcification (discovery cohort, N=0; replication cohort, N=5); discontinuing the scan before completing all three tasks (discovery cohort, N=4; replication cohort, N=7); superthreshold head motion on one or more tasks (discovery cohort, N=11; replication cohort, N=5); technical reasons, such as problems with stimulus display on one or more tasks (discovery cohort, N=2; replication cohort, N=4); or low behavioral performance on either the inhibition or reward task (<75% of trials receiving a button press, indicating sleepiness or low effort; discovery cohort, N=8; replication cohort, N=8). The analysis therefore included 69 participants in the discovery cohort and 77 in the replication cohort.
Statistical modeling of the fMRI data and region of interest definitions are detailed in the online supplement. Regions of interest were defined anatomically and included the left and right amygdala, insula, subgenual anterior cingulate cortex (sgACC), dACC (threat: fearful > neutral faces), nucleus accumbens (NAcc), OFC, amygdala (reward: monetary gain > loss), hippocampus, and vmPFC (inhibition: no-go < go).
Fear-potentiated startle.
Participants completed fear-acquisition and extinction tasks on the same day as the MRI scan. Details on the data acquisition and paradigm are presented in the online supplement.
Clustering Analysis
Analyses were conducted in R, version 3.6.3, with RStudio, version 1.2.1335. All tests were two-tailed and used a significance threshold of 0.05, with family-wise error correction as noted in the Results section. Clustering was conducted on data from the regions of interest extracted from the three fMRI tasks, using hierarchical agglomerated clustering, with the cluster package, version 2.1.0, following Ward’s criterion (agnes function). This is a bottom-up method designed to preserve the existing structure of the data without imposing assumptions of linearity, appropriate for exploratory analysis. The optimal number of clusters was determined using silhouette ( 31) and distance ( 32) methods. Nonparametric bootstrapping using the fpc package, version 2.2.5, was applied to the cluster solutions, with 1,000 iterations. After the initial hierarchical clustering, the data were randomly resampled with replacement. In each bootstrap, clustering was performed on the resampled data, and the new cluster most similar to each original cluster was identified by saving the maximum Jaccard coefficient (indexing similarity) for each old-new comparison ( 33). This was repeated, and a mean permuted Jaccard coefficient was computed across all the bootstraps by cluster. Permuted Jaccard coefficient therefore represents the proportion of individuals from each original cluster solution that were again clustered together in the permuted data. A permuted Jaccard coefficient of 0.6–0.75 indicates stable clusters, >0.75 represents high stability, and <0.50 is thought to indicate cluster instability ( 33); clusters were considered reconstituted on any bootstrap with a permuted Jaccard coefficient >0.60.
The replication was assessed quantitatively using a train and test approach. We trained a simple k-nearest-neighbors (knn function) model ( 34) with the class, version 7.3, package using the discovery cohort data, labeled using the cluster labels from the hierarchical clustering solution. We applied this knn model (“test”) in the replication cohort to obtain a new set of labels. We then compared these new labels to the de novo hierarchical clustering of the replication cohort in caret, version 6.0 ( 35).
Analysis of Posttrauma Outcomes by Cluster
Because clustering produced some small cells within-cohort (Ns as low as 11), we combined the clusters that replicated between the discovery and replication cohorts for further characterization. Cluster assignments from the initial cluster solutions were retained, rather than reclustering in a combined data set.
Chi-square tests (categorical variables) or one-way analyses of variance (ANOVAs) (continuous variables) were used to assess whether demographic factors or trauma-related factors differed between the cluster groups.
Given the multiple overlapping adverse mental health outcomes of trauma, we used multivariate analysis of variance (MANOVA) to test whether the pattern of subsequent mental health outcomes varied across the cluster groups. The outcome was a vector of standardized scores for PTSD symptoms, depression symptoms, dissociation, anxiety, and impulsivity. Predictors included cluster, assessment time point (time-invariant term for the 2-week, 8-week, 3-month, and 6-month posttrauma assessments; linear and quadratic terms), cluster-by-time point interaction, cohort, and a random effect for participant. Post hoc tests separating each outcome type were conducted using linear mixed models in the lme4 package, with the same set of predictors used in the MANOVA. To test whether the fMRI-based clustering provided incremental information above and beyond pretrauma symptom levels, we conducted secondary models including pretrauma symptoms that participants reported in the emergency department. Initial AURORA study findings indicated that among sociodemographic risk factors, pretrauma symptom levels were the strongest predictor of later PTSD and depression symptom severity ( 36, 37).
Finally, we tested whether a dimensional model of the 2-week fMRI data outperformed cluster assignment in predicting posttrauma outcomes using the first three principal components from the principal component analysis of the nine regions of interest in the combined discovery and replication data set as dimensional predictors of later outcomes. These models used the same structure as the cluster-based analyses. Model fit for dimensional and cluster-based models were directly compared.
Fear-potentiated startle during fear conditioning and extinction was also collected on the same day as the MRI scan (acquisition details are presented in the online supplement). ANOVAs tested whether the cluster groups varied in fear-potentiated startle responses during either acquisition or extinction as a function of cluster, cohort, block, and conditioned stimulus (CS) type and interactions between cluster, block, and CS type.
RESULTS
Covariance Among the fMRI Tasks and Regions of Interest
To assess for feature redundancy, we examined the covariance structure between the tasks and regions of interest. Regions of interest showed positive within-task covariance but not between tasks ( Figure 2A,B). The small to moderate correlations suggested that each task and region would contribute unique variance to a clustering analysis. Interestingly, similar regions were uncorrelated from one task to another; for example, participants’ amygdala reactivity to threat was not correlated with amygdala reactivity to reward (r=0.00, p=0.97). This suggests that the “crud” factor (everything correlates with everything) ( 38) was very low across this set of tasks.
Clustering of Individuals Using Task-Based fMRI 2 Weeks Posttrauma
Hierarchical clustering was first applied in the discovery sample, with 69 survivors of motor vehicle accidents. A four-cluster solution was identified ( Figure 2C; see also Figure S1A,B in the online supplement). Silhouette results suggested an optimal clustering with two groups (k=2) but with only a small decrement in width for k=4, whereas Hartigan’s distance metric showed an optimal gain in cluster cohesiveness at k=4. Examination of fMRI activation patterns ( Figure 2E,G) indicated that individuals in cluster 1 (permuted Jaccard coefficient=0.52) showed high reactivity to both threat and reward, with little engagement of regulatory regions in threat or inhibition. Thus, we classified the cluster 1 group as “reactive/disinhibited.” Individuals in cluster 2 (permuted Jaccard coefficient=0.54) showed threat responsivity predominated by the sgACC but low reward reactivity, and we classified this group as “low reward/high threat.” Individuals in cluster 3 (permuted Jaccard coefficient=0.52) showed no reactivity to threat, nor engagement of the vmPFC or hippocampus during inhibition, but very high responsivity to reward, and we classified this group as “high reward.” Finally, individuals in cluster 4 (permuted Jaccard coefficient=0.56) showed marked deactivation to threat in the amygdala, dACC, and insula, some activation of the hippocampus in the inhibition task, and little reactivity to reward, and we classified this group as “inhibited.”
The replication cohort included 77 participants with a variety of different trauma types, including interpersonal traumas. Here, the most favorable clustering solution included three groups, with agreement between the silhouette and distance metrics at k=3 ( Figure 2D; see also Figure S1C and D in the online supplement). The groups appeared to be consistent with cluster 1 (reactive/disinhibited; permuted Jaccard coefficient=0.53), cluster 2 (low reward/high threat; permuted Jaccard coefficient=0.73), and cluster 4 (inhibited; permuted Jaccard coefficient=0.55) from the discovery sample ( Figure 2F,H). There was a striking absence of a high reward-like phenotype; individuals who showed high reward reactivity also showed high threat reactivity.
The inclusion of higher-impact traumas may have pushed reward-responsive individuals toward higher threat reactivity. To test this, we combined both cohorts and examined effects of either injury severity or interpersonal violence on threat reactivity in the amygdala. Injury severity positively predicted amygdala reactivity (F=4.58, df=1, 144, p=0.03; see also Figure S2 in the online supplement), whereas interpersonal compared with noninterpersonal trauma did not (p=0.61). High reward was therefore likely subsumed under the reactive/disinhibited phenotype, related to higher-acuity traumas.
In the quantitative assessment of replication for clusters 1, 2, and 4, we assessed the extent to which a model trained on the clustering solution from the discovery cohort could predict the clustering solution within the replication cohort. The model trained on the discovery cohort data had 65.0% (95% CI=53.2, 75.5) accuracy in predicting the original hierarchical clustering-based labels in the replication cohort, compared with a 45.4% no-information rate (p=0.0005, kappa=0.45). This indicated that the clustering solution in the replication cohort could be recapitulated above and beyond chance levels using only the features of the discovery cohort solution.
The clusters were unrelated to demographic, health-related, trauma-related, or site-specific factors in follow-up testing (see the online supplement) and thus appeared to reflect covert neurocognitive features.
Prospective Trajectories of Mental Health Among the Four Clusters
We next assessed trauma-related outcomes across the clusters in a combined sample of 125 individuals from the clusters that replicated across both cohorts (see Figure S3 in the online supplement). The clusters showed different multivariate symptom profiles posttrauma (F=2.25, df=2, 948, p=0.013) ( Figure 3A). Although assessment time point was included as a factor in the model, there was no interaction of cluster by time point on the symptom profile (p=0.82). Follow-up tests were then performed for each symptom type separately. First, there was an effect of cluster on the longitudinal model of PTSD symptoms (Wald χ 2=6.47, p=0.039), with the highest symptoms in the reactive/disinhibited cluster ( Figure 3B). The effect of cluster was reduced when pretrauma PTSD symptoms were added as a predictor in the model (cluster effect: χ 2=5.10, p=0.078; pretrauma PTSD severity effect: χ 2=20.19, p<0.001).
Second, there was an effect of cluster on the longitudinal model of anxiety symptoms (Wald χ 2=6.23, p=0.044) that was higher in the reactive/disinhibited cluster ( Figure 3E). This effect held after including pretrauma anxiety symptoms as a predictor in the model (cluster effect: χ 2=6.07, p=0.048; pretrauma anxiety severity effect: χ 2=72.37, p<0.001), suggesting that cluster information provided unique predictive value above baseline symptoms. Finally, cluster did not predict depressive (p=0.19) or dissociative symptoms (p=0.86) or impulsivity (p=0.96) ( Figure 3C,D,F).
Cluster-Based Compared With Dimensional Models for Predicting Longitudinal Posttrauma Outcomes
We tested the utility of the discrete clusters against a dimensional model of the fMRI data for predicting longitudinal trajectories of stress-related symptoms. Dimensional fMRI predictors were continuous covariates reflecting threat reactivity, reward reactivity, and inhibition. Models with these covariates as predictors of later posttrauma symptoms showed negligible improvement in the model fit over the cluster-based models (see Table S3 in the online supplement). The individual fMRI dimensions were not linearly associated with any posttrauma outcome, with the exception of a negative association between the inhibition-related fMRI dimension and later dissociative symptoms (p=0.044). In models that included both clusters and dimensions competing for the variance in posttrauma outcomes, cluster assignment still predicted subsequent PTSD symptoms (p=0.012), whereas the three fMRI dimensions did not (all p values >0.05). For anxiety, neither cluster (p=0.080) nor dimensions (all p values >0.05) were significant in the head-to-head model. In summary, the dimensional model did not provide better predictive value than the cluster-based models.
Convergent Validity With Fear-Learning Phenotypes
On the day of the MRI scan, participants also completed a fear-potentiated startle paradigm that included fear conditioning and extinction. During fear conditioning, effects of CS (F=12.73, df=1, 468, p=0.0004) and the CS-by-block interaction (F=6.06, df=1, 468, p=0.003) suggested that fear conditioning occurred and that discrimination between the CS+ and CS− developed across acquisition. There was a significant cluster-by-block interaction (F=4.13, df=4, 468, p=0.003), such that the low reward/high threat cluster showed the highest fear-potentiated startle responses to both the CS+ and CS− at the beginning of fear conditioning, but this cluster was comparable to the other cluster groups by the end of the task ( Figure 3G). There was no cluster-by-CS interaction (p=0.82). During fear extinction, startle responses to the CS+ and CS− showed the expected decline over time (block effect: F=28.79, df=3, 623, p=0.02 × 10 −14), indicating the presence of extinction learning, but there was no interaction of CS by block (p=0.64), indicating no difference in the extinction pattern for CS+ compared with CS−. This was consistent with findings from previous studies of chronic PTSD using this startle paradigm ( 39). There was again an interaction of cluster by block, indicating different rates of extinction in the different cluster groups (F=2.35, df=6, 623, p=0.03). The low reward/high threat cluster showed the highest fear-potentiated startle responses to both the CS+ and CS− at the beginning of extinction, decreasing to become comparable to the other cluster groups by the end of the task ( Figure 3H).
Voxel-Wise Whole-Brain Comparison of Cluster Groups
Finally, to identify brain regions outside the primary regions of interest included in the clustering, we conducted whole-brain analysis in the combined sample of 125 participants ( Figure 4; see also Table S4 in the online supplement), comparing the three replicated cluster groups within the threat, reward, and inhibition fMRI tasks. The reactive/disinhibited cluster showed greater activation than the other two cluster groups in a mesopontine cluster overlapping with the median raphe nucleus and ventral tegmental area, as well as the hypothalamus, dACC, and insula in response to threat cues. In contrast, the low reward/high threat cluster showed greater activation in the left and right amygdala, hippocampus, and insula in response to threat cues. The reactive/disinhibited cluster also showed greater reactivity than the other two cluster groups in the amygdala, hippocampus, and rostral anterior cingulate cortex in response to reward. The inhibited cluster showed no region of greater activation compared with the other two cluster groups.
DISCUSSION
In a well-characterized cohort followed longitudinally in the aftermath of trauma, we identified participant clusters in a manner that was agnostic to standard diagnostic categories for posttrauma outcomes using fMRI across several neurocognitive dimensions of interest, including threat, reward, and inhibition. In the discovery cohort of motor vehicle accident survivors, four clusters were observed, and three were replicated in a cohort with a wider variety of index traumas. Given the timing at 2 weeks posttrauma, the clustering likely reflects a combination of traits that predate the trauma, as well as acute stress responses in the wake of the traumatic event. The clusters were not related to the demographic characteristics and background variables of the participants (e.g., gender, childhood trauma) but could still plausibly reflect pretrauma factors, such as genetics, family history, or temperament ( 40). Our findings confirmed the hypothesis that clusters may be associated with different posttrauma outcomes: different longitudinal patterns emerged over the first 6 months posttrauma, with the reactive/disinhibited cluster associated with subsequent heightened symptoms of PTSD/hypervigilance and anxiety. In addition, these findings represent an important step toward defining a neuroimaging-based longitudinal prediction model for stress-related resilience and vulnerability.
Our results suggest that an unsupervised forward-inference model is tractable for modeling heterogeneity in stress-related psychiatric outcomes, despite the lack of constraints on the model. We used very simple and transparent tools for clustering, uninformed by psychiatric symptoms or diagnoses, and found strong evidence for overlap (65%) in the cluster solutions arising from two fully independent hierarchical clustering solutions in different subcohorts. This level of overlap is consistent with cluster replication levels seen in larger brain-based clustering efforts, such as in the Adolescent Brain Cognitive Development study ( 41), and suggests that multivariate task-based fMRI data contain consistent information about individual differences. In fact, a major value of this study is its demonstration of task-based fMRI as a useful tool for mapping psychiatric heterogeneity. Task-based fMRI may not be needed to identify unitary biomarkers, such as for a diagnosis (e.g., PTSD) or a symptom (e.g., hyperarousal), where resting-state MRI is likely preferable for its rich information on neural circuit function and low barriers to translation. However, for the purposes of resolving heterogeneity, task fMRI shows clear strengths. For example, the signal from the amygdala was clearly uncorrelated in threat tasks compared with reward tasks ( Figure 2A,B), and this variability was important in resolving the clusters. Such information would not have been apparent from analyses of intrinsic network activity.
Interpretation of the Biotypes
The reactive/disinhibited cluster was the most interesting candidate as a risk-related biotype. Individuals in this cluster showed threat hyperreactivity, particularly in the insula and dACC, accompanied by high NAcc reward reactivity, as well as higher subsequent PTSD symptoms. This was partly contrary to previous findings showing that lower NAcc reward response predicted subsequent PTSD ( 16). Threat and reward reactivity have rarely been assessed concurrently in previous studies of trauma and related outcomes. However, preclinical findings indicate that interacting pathways regulate both threat and reward reactivity. For example, stress-related hyperactivity of the basolateral amygdala can directly influence NAcc function via direct efferent projections ( 42) and change reward-seeking behaviors ( 43). Participants in the reactive/disinhibited cluster also showed greater activation in the reticular nuclei (median raphe, ventral tegmental area) during the threat task. The reticular formation stimulates wakefulness and arousal ( 44, 45). The role of these nuclei in dopamine and serotonin synthesis may point toward tailored intervention opportunities for the future; early studies are exploring dopaminergic modulation to treat PTSD ( 46).
The fMRI features observed in the low reward/high threat cluster included moderate responsivity to threat dominated by sgACC activation, along with markedly low reactivity to reward. Reduced reward responsivity in the NAcc, amygdala, and OFC is characteristic of major depressive disorder ( 47), as is sgACC hyperreactivity to sadness-inducing stimuli. The heightened fear-potentiated startle shown by this group during early fear acquisition and extinction is consistent with patterns previously observed in comorbid PTSD and major depressive disorder ( 27), which is more common posttrauma than each disorder alone ( 48). This cluster group also showed greater threat reactivity in the amygdala, insula, and hippocampus in whole-brain analyses compared with the other cluster groups. Together, the findings suggest the possibility that PTSD-related symptom groups may be divided into a low reward/high threat group driven more by cortical function, and a reactive/disinhibited group, driven more by brainstem nuclei.
Participants in the inhibited cluster appeared to be most consistent with active coping, with low threat reactivity accompanied by relatively high vmPFC and hippocampus engagement during inhibition. Individual features of this pattern have previously been associated with resilience. For example, lower amygdala threat reactivity predicts lower future PTSD symptoms ( 12, 13). Similarly, greater vmPFC and hippocampal activation during inhibition has been associated with resilience ( 22, 26). Our findings indicate that some individuals show the combined profile in the acute posttrauma period. However, this was the smallest cluster, and it is possible that inhibition as an adaptation to the stress of the index trauma had not yet fully emerged by 2 weeks posttrauma. As data collection continues, the AURORA study will include additional neuroimaging 6 months posttrauma, allowing a window into the further development of these profiles.
Participants in the high reward cluster showed high reward reactivity, low threat reactivity, and low inhibition, a pattern suggesting preserved positive affect in the context of low top-down regulation. However, this pattern may only be observed when the emotional impact of trauma is relatively low, and this cluster was not observed in the replication cohort.
Limitations
To reduce participant burden, symptom assessments were abbreviated and based on self-report. This limited our ability to directly compare the outcomes to gold-standard assessments of psychiatric disorders. However, future work could apply these biotypes to archival data with interview-based assessments to extend our findings. Additionally, with sample sizes of 69 participants in the discovery cohort and 77 in the replication cohort, the study may be underpowered. However, we are encouraged that three clusters were replicated, indicating generalizability even across different types of trauma.
CONCLUSIONS
Neuroimaging phenotypes emerging in the early aftermath of trauma are associated with risk of or resilience to trauma-related psychopathology. Contrary to our initial predictions that heightened threat and blunted reward reactivity may reflect stress vulnerability, a cluster showing heightened reactivity to both threat and reward was associated with the subsequent maintenance of the highest levels of PTSD symptoms. Heightened reward reactivity in the early aftermath of a major stressor may be an underexplored risk mechanism for the development of stress-related disorders. The biotypes identified here, with further development to assess normative values and precision, may provide important information about targeted interventions to address different forms of future stress-related psychopathology.
1 : Diathesis-stress theories in the context of life stress research: implications for the depressive disorders . Psychol Bull 1991 ; 110 : 406 – 425 Crossref, Medline, Google Scholar
2 : The role of genes, stress, and dopamine in the development of schizophrenia . Biol Psychiatry 2017 ; 81 : 9 – 20 Crossref, Medline, Google Scholar
3 : Cognitive behavioral therapy for insomnia in posttraumatic stress disorder: a randomized controlled trial . Sleep (Basel) 2014 ; 37 : 327 – 341 Crossref, Medline, Google Scholar
4 : Stress Resilience. Amsterdam , Elsevier , 2020 , pp 197 – 207 Crossref, Google Scholar
5 : 636,120 ways to have posttraumatic stress disorder . Perspect Psychol Sci 2013 ; 8 : 651 – 662 Crossref, Medline, Google Scholar
6 : Understanding heterogeneity in PTSD: fear, dysphoria, and distress . Depress Anxiety 2014 ; 31 : 97 – 106 Crossref, Medline, Google Scholar
7 : The AURORA Study: a longitudinal, multimodal library of brain biology and function after traumatic stress exposure . Mol Psychiatry 2020 ; 25 : 283 – 296 Crossref, Medline, Google Scholar
8 : Resting-state connectivity biomarkers define neurophysiological subtypes of depression . Nat Med 2017 ; 23 : 28 – 38 Crossref, Medline, Google Scholar
9 : Biotypes of functional brain engagement during emotion processing differentiate heterogeneity in internalizing symptoms and interpersonal violence histories among adolescent girls . J Psychiatr Res 2020 ; 121 : 197 – 206 Crossref, Medline, Google Scholar
10 : Identification of distinct psychosis biotypes using brain-based biomarkers . Am J Psychiatry 2016 ; 173 : 373 – 384 Link, Google Scholar
11 : Evaluating the evidence for biotypes of depression: methodological replication and extension of . Neuroimage Clin 2019 ; 22 : 101796 Crossref, Medline, Google Scholar
12 : Amygdala reactivity and anterior cingulate habituation predict PTSD symptom maintenance after acute civilian trauma . Biol Psychiatry 2017 ; 81 : 1023 – 1029 Crossref, Medline, Google Scholar
13 : Human vulnerability to stress depends on amygdala’s predisposition and hippocampal plasticity . Proc Natl Acad Sci USA 2009 ; 106 : 14120 – 14125 Crossref, Medline, Google Scholar
14 : Amygdala response to negative stimuli predicts PTSD symptom onset following a terrorist attack . Depress Anxiety 2014 ; 31 : 834 – 842 Crossref, Medline, Google Scholar
15 : Cortical salience network activation precedes the development of delusion severity . Psychol Med 2016 ; 46 : 2741 – 2748 Crossref, Medline, Google Scholar
16 : Imbalanced neural responsivity to risk and reward indicates stress vulnerability in humans . Cereb Cortex 2013 ; 23 : 28 – 35 Crossref, Medline, Google Scholar
17 : The brain’s response to reward anticipation and depression in adolescence: dimensionality, specificity, and longitudinal predictions in a community-based sample . Am J Psychiatry 2015 ; 172 : 1215 – 1223 Link, Google Scholar
18 : Blunted ventral striatum development in adolescence reflects emotional neglect and predicts depressive symptoms . Biol Psychiatry 2015 ; 78 : 598 – 605 Crossref, Medline, Google Scholar
19 : Amygdala reward reactivity mediates the association between preschool stress response and depression severity . Biol Psychiatry 2018 ; 83 : 128 – 136 Crossref, Medline, Google Scholar
20 : Ventromedial prefrontal cortex processing during emotional evaluation in late-life depression: a longitudinal functional magnetic resonance imaging study . Biol Psychiatry 2008 ; 64 : 349 – 355 Crossref, Medline, Google Scholar
21 : Neural correlates of recovery from post-traumatic stress disorder: a longitudinal fMRI investigation of memory encoding . Neuropsychologia 2011 ; 49 : 1771 – 1778 Crossref, Medline, Google Scholar
22 : The role of the hippocampus in predicting future PTSD symptoms in recently traumatized civilians . Biol Psychiatry 2018 ; 84 : 106 – 115 Crossref, Medline, Google Scholar
23 : Neural indicators of interpersonal anger as cause and consequence of combat training stress symptoms . Psychol Med 2017 ; 47 : 1561 – 1572 Crossref, Medline, Google Scholar
24 : Hippocampal activation during contextual fear inhibition related to resilience in the early aftermath of trauma . Behav Brain Res 2021 ; 408 : 113282 Crossref, Medline, Google Scholar
25 : Savoring the past: positive memories evoke value representations in the striatum . Neuron 2014 ; 84 : 847 – 856 Crossref, Medline, Google Scholar
26 : Reduced neural activation during an inhibition task is associated with impaired fear inhibition in a traumatized civilian sample . Cortex 2013 ; 49 : 1884 – 1891 Crossref, Medline, Google Scholar
27 : Impaired fear inhibition is a biomarker of PTSD but not depression . Depress Anxiety 2010 ; 27 : 244 – 251 Crossref, Medline, Google Scholar
28 : Neurobiological basis of failure to recall extinction memory in posttraumatic stress disorder . Biol Psychiatry 2009 ; 66 : 1075 – 1082 Crossref, Medline, Google Scholar
29 : fMRIPrep: a robust preprocessing pipeline for functional MRI . Nat Methods 2019 ; 16 : 111 – 116 Crossref, Medline, Google Scholar
30 : Evaluation of ICA-AROMA and alternative strategies for motion artifact removal in resting state fMRI . Neuroimage 2015 ; 112 : 278 – 287 Crossref, Medline, Google Scholar
31 : Silhouettes: a graphical aid to the interpretation and validation of cluster analysis . J Comput Appl Math 1987 ; 20 : 53 – 65 Crossref, Google Scholar
32 : Clustering Algorithms. Hoboken, NJ , John Wiley and Sons , 1975 Google Scholar
33 : Cluster-wise assessment of cluster stability . Comput Stat Data Anal 2007 ; 52 : 258 – 271 Crossref, Google Scholar
34 : Pattern Recognition and Neural Networks. Cambridge, UK , Cambridge University Press , 2007 Google Scholar
35 : Building predictive models in R using the caret package . J Stat Softw 2008 ; 28 : 1 – 26 Crossref, Medline, Google Scholar
36 : Socio-demographic and trauma-related predictors of depression within eight weeks of motor vehicle collision in the AURORA study . Psychol Med (Online ahead of print, October 29, 2020) Google Scholar
37 : Socio-demographic and trauma-related predictors of PTSD within 8 weeks of a motor vehicle collision in the AURORA study . Mol Psychiatry (Online ahead of print, October 19, 2020) Google Scholar
38 : Crud (re)defined . Adv Methods Pract Psychol Sci 2020 ; 3 : 238 – 247 Crossref, Medline, Google Scholar
39 : Fear extinction in traumatized civilians with posttraumatic stress disorder: relation to symptom severity . Biol Psychiatry 2011 ; 69 : 556 – 563 Crossref, Medline, Google Scholar
40 : The complex neurobiology of resilient functioning after childhood maltreatment . BMC Med 2020 ; 18 : 32 Crossref, Medline, Google Scholar
41 : Decomposing complex links between the childhood environment and brain structure in school-aged youth . Dev Cogn Neurosci 2021 ; 48 : 100919 Crossref, Medline, Google Scholar
42 : Topographical organization of amygdaloid projections to the caudatoputamen, nucleus accumbens, and related striatal-like areas of the rat brain . Neuroscience 1991 ; 44 : 15 – 33 Crossref, Medline, Google Scholar
43 : Basolateral amygdala and stress-induced hyperexcitability affect motivated behaviors and addiction . Transl Psychiatry 2017 ; 7 : e1194 – e1194 Crossref, Medline, Google Scholar
44 : An ascending reticular activating system in the brain stem . AMA Arch Neurol Psychiatry 1952 ; 67 : 145 – 154, discussion 167–171 Crossref, Medline, Google Scholar
45 : Neuroanatomic connectivity of the human ascending arousal system critical to consciousness and its disorders . J Neuropathol Exp Neurol 2012 ; 71 : 531 – 546 Crossref, Medline, Google Scholar
46 : L-DOPA and consolidation of fear extinction learning among women with posttraumatic stress disorder . Transl Psychiatry 2020 ; 10 : 287 Crossref, Medline, Google Scholar
47 : Reduced caudate and nucleus accumbens response to rewards in unmedicated individuals with major depressive disorder . Am J Psychiatry 2009 ; 166 : 702 – 710 Link, Google Scholar
48 : Comorbidity of 9/11-related PTSD and depression in the World Trade Center Health Registry 10–11 years postdisaster . J Trauma Stress 2014 ; 27 : 680 – 688 Crossref, Medline, Google Scholar