Abstract
Stress is implicated in psychosis etiology and exacerbation, but pathogenesis toward brain network alterations in schizophrenia remain unclear. White matter connects limbic and prefrontal regions responsible for stress response regulation, and white matter tissues are also vulnerable to glucocorticoid aberrancies. Using a novel psychological stressor task, we studied cortisol stress responses over time and white matter microstructural deficits in schizophrenia spectrum disorder (SSD). Cortisol was measured at baseline, 0-, 20-, and 40-min after distress induction by a psychological stressor task in 121 SSD patients and 117 healthy controls (HC). White matter microstructural integrity was measured by 64-direction diffusion tensor imaging. Fractional anisotropy (FA) in white matter tracts were related to cortisol responses and then compared to general patterns of white matter tract deficits in SSD identified by mega-analysis. Differences between 40-min post-stress and baseline, but not acute reactivity post-stress, was significantly elevated in SSD vs HC, time × diagnosis interaction F2.3,499.9 = 4.1, p = 0.013. All SSD white matter tracts were negatively associated with prolonged cortisol reactivity but all tracts were positively associated with prolonged cortisol reactivity in HC. Individual tracts most strongly associated with prolonged cortisol reactivity were also most impacted in schizophrenia in general as established by the largest schizophrenia white matter study (r = −0.56, p = 0.006). Challenged with psychological stress, SSD and HC mount similar cortisol responses, and impairments arise in the resolution timeframe. Prolonged cortisol elevations are associated with the white matter deficits in SSD, in a pattern previously associated with schizophrenia in general.
Similar content being viewed by others
Introduction
Schizophrenia spectrum disorder (SSD) is a severe mental illness characterized by psychosis, negative symptoms, and cognitive impairment. Its etiology remains unknown but likely includes stress-related pathophysiology, which is linked to the onset and exacerbation of psychosis, symptom severity [1], and prognosis [2]. Stress and adversity are also associated with elevated and dysregulated cortisol [3, 4]; and prolonged elevation of glucocorticoid levels is associated with brain functional and structural changes [5]. However, evidence bridging prolonged elevations of glucocorticoids, an abnormal stress response, and structural brain changes in SSD are limited.
White matter (WM) may be particularly vulnerable to pathogenic effects of stress, as oligodendrocytes and myelination are adversely affected by glucocorticoids [6]. Chronically elevated glucocorticoids impair WM reorganization following injury while promoting gliosis [7, 8]. Cortisol may also indirectly influence WM integrity through pro-inflammatory cytokines [9, 10] leading to neuronal and oligodendrocyte cellular and receptor level dysfunction and damaging myelination patterns in pre-clinical studies [8, 11]. In clinical studies, chronic stress is associated with WM lesions [12] and reduced fractional anisotropy (FA), an index of WM microstructural integrity using diffusion tensor imaging (DTI) [13]. Importantly, reduced FA is consistently associated with SSD patients [14], evident even in antipsychotic-naive first-episode patients [15] and in non-ill, first-degree relatives [16]. These WM changes are present in a regionally specific pattern, with the most affected WM tracts at the interhemispheric callosal and frontal fibers such as the corpus callosum (CC) and the anterior corona radiate (ACR), while other regions like the corticospinal tract (CST) are not [17, 18]. Although not traditionally thought of as tracts with direct stress responsiveness, these WM tracts connect prefrontal cortical and limbic structures that are responsible for mediating HPA-axis responses to stress [19, 20]. Tracts more specifically associated with limbic regional connectivity such as the fornix and the cingulate gyrus associated cingulum (CGC) are also among the more affected WM tracts in SSD [21, 22]. Abnormal cortisol reactivity in SSD could also represent a communication failure in these WM structures [23].
Despite extensive studies, the underlying mechanism of the decreased WM integrity in SSD remains elusive. We hypothesized that aberrant stress-induced glucocorticoid responses in SSD are partly responsible for these reported WM deficits. Furthermore, as discussed above, not all WM tracts are equally affected in SSD [18, 24]. If aberrant stress responses in SSD are contributory to WM deficits, we expect that regional tracts more vulnerable in SSD and/or more closely associated with limbic functions are also more strongly associated with the stress response abnormalities. Therefore, we hypothesized that inadequate or prolonged resolution of cortisol levels in response to stress would be associated with WM alterations in a tract-specific manner involving regions, in part, responsible for limbic communications previously mentioned.
We have previously examined the first hypothesis in 30 patients and 33 controls, although the sample size was limited and we did not find a statistically significant stress by diagnosis interaction [14]. With a larger cohort of participants, we re-examined the hypothesis and additionally examined the hypothesis that abnormal glucocorticoid response to psychological stress is associated with schizophrenia-related regional vulnerability patterns in the WM.
Materials and methods
Participants
Participants with SSD (n = 121) and HC (n = 117) completed the stress challenge task and DTI brain imaging; the sample included the 30 patients and 33 controls in the preliminary report [14]. SSD patients (n = 121: 84 males; age = 36.6 ± 13.3 years [SD]) were recruited from outpatient clinics at the Maryland Psychiatric Research Center and neighboring mental health outpatient clinics. HC (n = 117: 60 males; age = 35.4 ± 13.6 years) were recruited through media advertisements. Time interval between the psychological stress test and imaging session was 6.0 ± 8.0 days, with no significant group difference (t = −1.3, p = 0.2). The Structured Clinical Interview for DSM-IV or −5 was utilized to obtain diagnoses of either schizophrenia or schizoaffective disorder, categorized as SSD. Controls had no current DSM severe mental illness diagnoses, however past single episode depression was not exclusionary. Controls had no family history of psychosis in the prior two generations. Current or past major medical and neurological illnesses, history of head injury with cognitive sequelae, intellectual disability, substance dependence within the past 6 months, or current substance abuse were exclusionary for all participants. Substance use disorder was based on SCID interview, while urine toxicology screening was performed for verification purposes when indicated to confirm current substance use and allowed for excluding participants based on results. Clinical and cognitive assessments were performed using: Brief Psychiatric Rating Scale (BPRS) and Brief Negative Symptom Scale (BNSS) [25]; Digit Symbol Substitution test (for processing speed) from the Wechsler Adult Intelligence Scale-3 [26] and Digit Sequencing task (for working memory) from the Brief Assessment of Cognition in Schizophrenia [27]. Antipsychotic regimens for participants were converted into CPZ equivalents: 9 on both typical and atypical, 69 on atypical, 10 on typical, 19 on clozapine, one on clozapine and typical, seven on clozapine and atypical, and six not on any antipsychotic. Demographic and clinical characteristics are summarized in Table 1. Protocols were approved by the University of Maryland IRB and written informed consent to each participant in the study was collected.
Psychological stressor tasks
The psychological stressor tasks include the Paced Auditory Serial Addition Task (PASAT) [28] and the Mirror-Tracing Persistence Task [29], order randomly assigned and completed in series, to study an individual’s distress tolerability [30]. Procedures were followed as previously established [31]. Briefly, in the PASAT, a computer mouse was used to select the correct sum of consecutive numbers presented briefly on a computer screen. The MTPT asked participants to trace a dot along the outline of a star-shaped image on the computer using the mouse. Response errors were met with a loud, ~90 decibel, aversive buzzer noise. Learning sessions were performed prior to the experimental one. Speed and accuracy of responses were measured, and an algorithm automatically titrated the task presentation to provide similar challenges across participants. Participants were incentivized to continue, but could terminate the sessions at any point, coded in binary as: 1 = completed one or both tests, “distress tolerant”; 0 = quit both tests, “distress intolerant”.
Participants provided saliva at four timepoints: prior to task initiation (baseline), immediately (t = 0), 20-, and 40-min post-task completion/quitting (Fig. 1A). Thirty minutes before the baseline sample collection and after the post 0-min collection timepoint, participants sat comfortably, watched a scenic video, and/or read leisurely. All testing sessions were held approximately between 1200 and 1600 h. Participants refrained from eating, drinking, or smoking for 1-h before testing. Saliva samples were placed on ice and immediately stored after the test at −80 °C until assay. Prior to assay, samples were thawed and centrifuged at 10,000 g for 10-min. Cortisol was assayed using a commercial enzyme immunoassay kit (Salimetrics), following the manufacturer recommended protocol. Intra-assay coefficient of variance (CV) was 7.77% and inter-assay CV was 4.18% in our lab. Salivary cortisol levels (µg/dL) are highly correlated with serum levels and typically peak at 30 min following the initiation of a stressor task [32]. Physiological cortisol reactivities were calculated as the difference between the concentration obtained at baseline and 0-, 20-, and 40-min post-task completion, as the immediate, acute, and prolonged cortisol reactivity, respectively, in agreement with previously published studies [14, 33, 34].
Diffusion tensor imaging (DTI)
All MR examinations were performed at the University of Maryland Center for Brain Imaging Research, using a Siemens 3T TRIO MRI (Erlangen, Germany) system equipped with a 32-channel phase array head coil. DTI data were collected using a single-shot, echo-planar, single refocusing spin-echo, T2-weighted sequence, with ‘GeneRalized Autocalibrating Partially Parallel Acquisitions’, acceleration factor 2 [35], yielding voxel dimensions 1.7 × 1.7 × 3.0 mm, acquisition time ~8 min. The sequence parameters were: echo time/repetition time = 87/8000 ms, field-of-view = 200 mm, axial slice orientation with 50 slices and no gaps, five b = 0 images and 64 isotropically distributed diffusion-weighted directions with b = 700 s/mm2. ENIGMA-DTI pipeline (https://www.nitrc.org/projects/enigma_dti) was used for tract-based spatial statistical analysis of diffusion anisotropy [36] with protocol demonstrating excellent reproducibility [37, 38]. FA images were created using previously published protocols [14, 39] from the JHU atlas [40] and then nonlinearly aligned to a group-wise, minimal-deformation target (MDT) brain using the FLIRT method [36]. Briefly, the diffusion tensor was fit to the motion and eddy current diffusion data. RMSDIFF [41] was used to estimate the root mean square (RMS) movement distance between diffusion sensitized and b = 0 images. The RMS difference was calculated by comparing two 4 × 4 transformation matrixes: a transformation matrix from each frame to the first b0 image and an identity matrix which served as the no-movement reference. The advantage of the RMS distance includes both translation and rotation effects, providing an index of motion for the whole brain [42]. Motion statistics are displayed in Table 1 and included in all covariate analysis. All data passed quality assurance control of 2.5 mm accumulated motion during the scan.
The QA/QC parameters were derived empirically as an adequate threshold in a test-retest examination of subject for motion, as described previously by our group and others [37]. The group’s MDT brain was identified by warping all individual brain images in the group to each other. Group-average FA image was used to create a group-wise skeleton of 24 bilaterally-averaged WM tracts, thresholded at FA = 0.20 level to eliminate non-WM voxels, and FA values were projected onto the group-wise skeleton to account for residual misalignment among individual WM tracts.
Statistical analysis
Group differences on demographics and clinical measures were examined with ANOVA or χ2. The primary measure of cortisol reactivity in response to the stressor challenge was analyzed by repeated measures (RM)-ANOVA, where the four timepoints were the repeated measures, and diagnosis was the between-subject factor. Post-hoc analyses for significant diagnosis x time interaction were then performed between individuals and baseline. Greenhouse–Geisser corrections were used to report RM-ANOVA. Correlation analyses between WM FA and psychological stressor task measures were corrected for motion, age, sex, and current cigarette smoking status, all features known to affect WM FA, and reported as partial r-values. Fisher’s r-to-z transformation was performed to compare group difference in correlation coefficients. All significance was two-tailed. False discovery rate (FDR) was used for multiple comparisons of 24 WM tracts with a q-value set at <0.05.
For WM regional vulnerability assessment, the ENIGMA study provided a meta-analysis of the case/control effects associated with schizophrenia in 23 major WM regions (Cohen’s d) [17] that represent the expected pattern of disease effects across different WM tracts as established by the largest DTI study in schizophrenia [18], and can be found in Supplementary Table 1. Correlational coefficients from the correlation analyses between each WM tract and psychological stressor task measures in the current study was then correlated with the ENIGMA tract effect sizes, which investigated whether tracts associated with the psychological stress are more likely tracts also affected by schizophrenia in general as determined by the ENIGMA study. The abbreviations of the WM tracts are also included in Supplementary Table 2.
Results
Behavioral and glucocorticoid response to the psychological stressor
SSD and HC participants were similar in age (t = −0.72; p = 0.47), but not sex (69% vs 52% male, respectively; χ2 = 7.78; p < 0.01) or current smoking status (43.8% vs 32.5%, respectively; χ2 = 3.23; p = 0.048) (Table 1).
SSD participants were significantly more likely to be distress intolerant than HC (32% compared to 18%, respectively; χ2 = 5.80; p = 0.016) (Fig. 1B), in agreement with our preliminary report [31]. At baseline, cortisol levels were not significantly different based on group (F1,221 = 0.37, p = 0.55). SSD and HC similarly increased cortisol immediately following psychological stressor, followed by declines from 0 to 40 min post-stress, with group differences in the final timepoint (Fig. 1A). RM-ANOVA demonstrated a significant time (F2.3,499.9 = 4.4, p = 0.005) and time x diagnosis (F2.3,499.9 = 4.1, p = 0.013) interaction, but no significant diagnosis main effect (F1,217 = 0.018, p = 0.9).
To identify the key timepoint for this interaction effect, we performed post-hoc analysis with a series of RM-ANOVA tests comparing baseline with each timepoint. Between baseline and immediately post-stress showed a significant time effect (F1,236 = 7.83, p = 0.006) without diagnosis (F1,236 = 0.61, p = 0.43) or time × diagnosis interactions (F1,236 = 0.024, p = 0.88). Between baseline and 20-min post-stress, there was no significant time effect (F1,231 = 0.48, p = 0.49), diagnosis (F1,231 = 0.03, p = 0.86), or time × diagnosis interaction (F1,231 = 1.58, p = 0.21). Between baseline and 40-min post-stress, however, there was a significant time (F1,223 = 10.43, p = 0.001) and time × diagnosis interaction (F1,223 = 4.30, p = 0.039), but no diagnosis main effect (F1,223 = 0.26, p = 0.61). We directly compared this cortisol reactivity interaction (post-stress 40-min minus baseline) by group differences, and SSD had significantly higher (less reduced) relative cortisol response (F1,221 = 4.7, p = 0.03) (Fig. 1A), with no change once including sex and smoking as covariates (F1,218 = 4.84, p = 0.03). However, 3 cortisol reactivity values were >3 standard deviation (SD) from the mean (n = 2 HC, n = 1 SSD, 2 slightly above 3 SD, one at 4 SD) and removing them did change the significant group differences (t = −1.4, p = 0.15) although still in the same trend. As no outliers were present in the individual cortisol timepoints and the cortisol reactivity outliers were ≤4 SD, these samples were retained for further analysis. Distress intolerance was also added as a covariate (distress intolerance vs. non-distress intolerance) because it was significantly different between groups (Fig. 1B), in agreement with our preliminary analysis [14], but was not significant (F1,218 = 0.03, p = 0.85), and the prolonged cortisol reactivity effect remained (F1,218 = 5.9, p = 0.02). Distress intolerance then subdivided groups (Fig. 1C) but no significant diagnosis by distress intolerance (DI vs. DT) effect (p = 0.7–0.9) or distress intolerance main effect (p = 0.33–0.52) emerged at any of the timepoints, suggesting cortisol levels primarily tracked with diagnosis. Subsequent analyses focused on prolonged cortisol reactivity.
As a complementary analysis, differences of cortisol levels from its peak to the 40-min timepoint (peak minus 40 min) were also compared between groups, and controls (0.089 ± 0.0054) showed significantly larger decline in cortisol level compared to SSD (0.059 ± 0.0054) (t = 2.6, p = 0.009), again suggesting a slower decline from the peak levels in SSD.
White matter and the prolonged cortisol response
Prolonged cortisol reactivity in SSD was significantly correlated with whole-brain averaged WM FA (Suppl. Fig. 1A), (r = 0.22, p = 0.02). This significant association was reduced if excluding the three cortisol reactivity values that fell between 3 and 4 standard deviations of the mean (r = −0.1, p = 0.2) (Supplementary Fig. 1B). Tract-level correlational analyses was also performed, and several major WM tracts in controls (Fig. 2A) and SSD (Fig. 2B) were observed to be significantly correlated with cortisol reactivity. These include the SLF (r = 0.26, p = 0.006), IFO (r = 0.25, p = 0.01), and CGC (0.2, p = 0.04) in the HC group and the Fx (r = −0.23, p = 0.01) in the SSD group. In comparing group tract-level correlation coefficients, correlations most significantly differed in SLF (z = 3.03, p = 0.002) and CGC (z = 2.71, p = 0.006), as well as trend level significance (not surviving FDR correction) in BCC (z = 2.12, p = 0.03), superior fronto-occipital fasciculus (SFO) (z = 2.1, p = 0.03), and Fx (z = 2.0, p = 0.04).
We also used whole-brain FA as a covariate for further analysis (Supplementary Fig. 2). The general pattern was altered but HC still demonstrated positive correlations between the prolonged cortisol reactivity and WM FA in most of the tracts, while SSD demonstrated negative correlations in more of the tracts. However, no tract showed significant correlations after covarying out whole-brain average FA. Tracts were then rank ordered and compared to patterns generated with versus without including whole-brain average FA as a covariate. Their ranks were highly correlated in both HC (r = 0.85, p ≤ 0.001) and SSD (r = 0.92, p < 0.001). Lastly, a Fisher’s r-to-z transformation for the correlation coefficient differences was performed, revealing significant differences remained in 2 of the 6 tracts that were significant both with and without covarying whole-brain FA: CGC (z = 2.42, p = 0.01) and SLF (z = 2.8, p = 0.005) between HC and SSD. FA values are in Supplementary Table 1 along with effect sizes to compare with ENIGMA findings. A correlation analysis between the effect sizes of the tracts in our sample and those of the ENIGMA found that the tract-wise effect sizes of the current sample were significantly correlated to those reported in ENIGMA (r = 0.44, p = 0.03).
Prolonged cortisol reactivity and white matter deficit patterns in schizophrenia
As the ENIGMA-Schizophrenia DTI study provided the rank of WM tract effect sizes (Cohen’s d) affected by schizophrenia in general, we studied how they may be related to prolonged cortisol reactivity (Fig. 3). Correlations between prolonged cortisol reactivity and different WM tracts were significantly associated with ENIGMA tract effect sizes in SSD (r = −0.56, p = 0.006) (Fig. 3F). This appears specific to the prolonged cortisol reactivity, as immediate (Fig. 3D) and acute (Fig. 3E) cortisol reactivities did not show statistically significant associations. This relationship was unobserved in HC (Fig. 3A–C).
Prolonged cortisol response and clinical characteristics
Distress intolerance was not significantly associated with whole-brain WM FA (r = 0.14, p = 0.16) or prolonged cortisol reactivity (r = 0.02, p = 0.86) in SSD. In SSD, prolonged cortisol reactivity was not associated with CPZ (r = 0.05, p = 0.64), age (r = 14, p = 0.16), sex (r = −0.02, p = 0.84), smoking status (r = 0.01, p = 0.92), BPRS (r = −0.03, p = 0.7), BNSS (r = −0.06, p = 0.5) or working memory (r = −0.06, p = 0.56). Prolonged cortisol, however, was significantly associated with processing speed (r = −0.23, p = 0.02), although this would not be significant after corrected for multiple comparisons. In HC, no significant correlations were found between prolonged cortisol response and age (r = 0.02, p = 0.9), sex (r = −0.05, p = −0.6), smoking (r = 0.02, p = 0.8), working memory (r = 0.04, p = 0.7), or processing speed (r = −0.04, p = 0.7). Furthermore, the correlation coefficients between cortisol reactivity and processing speed were significantly different between SSD and controls (z = −2.09, p = 0.036).
Discussion
Our study found that the laboratory-administered psychological stressor task successfully induced cortisol secretion in SSD and HC similarly. At 40 min post-stress, SSD demonstrated a significantly elevated cortisol level compared to that of HC. This prolonged cortisol reactivity was significantly correlated with WM FA. Nearly all tracts were oppositely associated with prolonged cortisol reactivity between groups: positively associated in HC but negatively associated in SSD. These findings strongly mapped onto the schizophrenia-vulnerable tracts previously established by mega-analysis.
Results demonstrate that SSD patients have abnormally elevated cortisol responses in the prolonged timeframe following acute challenge, despite similar baseline and initial post-stress cortisol levels compared with HC. Cortisol levels in HC demonstrate the expected pattern of daytime diurnal decline [34] post-stress and drop rapidly after the peak, and reach well below the baseline level at 40 min, while in SSD the cortisol levels at 40 min post-stress essentially only return approximately to the baseline level (Fig. 1). Cortisol response to stress in SSD is controversial in the literature. Earlier diathesis-stress models [43] have employed the Trier Social Stress Test [44] to examine this hypothesis, although studies produced variable, inconclusive results [45,46,47]. Consistent with the diathesis-stress model, SSD demonstrates increased basal morning [48] and decreased overall daily cortisol levels [48], flattened awakening responses [49], and sensitivity to early life psychosocial stressors [50, 51]. Meta-analysis suggests SSD displays overall similar subjective responses to stress ratings and no group differences in cortisol response compared to healthy controls [52], consistent with our immediate or acute timepoint findings. Our results instead support a prolonged elevation of cortisol responses after psychological distress in support of the diathesis-stress model for SSD. Impairments in recovery and/or sustained cortisol release after psychological stressors may serve a key component of the diathesis-stress dysfunction in SSD, reconciling contradictory reports in the field.
Linking reduced FA and elevated prolonged cortisol responses in SSD here supports the association of chronic stress with WM lesions [12] and microstructural damage [13]. Regulating increased cortisol from stress demands integrates prefrontal, amygdalar, and hippocampal [19] regions, impaired in SSD [53]. Prolonged cortisol levels could represent failure of upstream HPA-axis structures to provide feedback inhibition [20]. Poor resolution of the initial cortisol response or prolonged release of additional cortisol after psychological stress, repeatedly occurring over time, may cumulatively lead to WM alterations in SSD shown here and across literature. Alternatively, WM deficits could be the original abnormalities impairing associative fibers such as the fornix, cingulum, frontal corona radiata, and the body of the CC connecting limbic, prefrontal, and the lateral hypothalamus [23, 54], leading to inadequately resolved cortisol responses to psychological stress. We found there to be significant group-level differences on comparing correlation coefficients between SSD and HC that were not present within the same group as significantly correlated WM tract and cortisol levels, for instance between the SFO or BCC. We interpret this finding in that there was still a relationship between WM and prolonged cortisol reactivity in these tracts, although the relationship was not strong enough to be statistically significant. However, the difference in this relationship was still significant between groups based on disease presence/absence.
Whole-brain average FA has been used as a covariate for tract-specific effects in some [17, 55] but not other [22, 56, 57] major DTI studies to account for a common variance in the microstructural integrity. Covarying whole-brain average FA markedly reduced and often eliminated significant group difference findings on specific tracts [17], as we have seen here, which are expected as whole-brain FA is an average of the values of multiple tracts. Meanwhile, we observed that the order representing the strength of tract-specific correlation coefficients between different regions and cortisol reactivity were largely retained with versus without the whole-brain FA covariate. The patient-control differences for their correlation coefficients remained significant in the CGC and SLF.
Understanding why only the prolonged, but not acute cortisol reactivity, was significantly correlated with FA may be clinically important. Pre-clinical models show that prolonged glucocorticoid exposure inhibits oligodendrocyte growth globally [6] and reduces myelin in a dose-dependent manner [7, 8]. Rats with maladaptive stress responses showed abnormally elevated cortisol levels and reduced FA following stress [55]. Meanwhile, acute cortisol elevations serve adaptive stress response mechanisms [56] crucial in neural plasticity [58] and activates key protein mediators for oligodendrocyte function [59]. Our findings between early and prolonged cortisol responses and WM appears consistent with these reports.
Opposite patterns of WM tracts and prolonged cortisol reactivity in SSD versus HC may reflect the mechanisms discussed above. Chronically elevated cortisol has been associated with widespread reductions in FA across WM tracts particularly in the CC and cingulate cingulum, in agreement with our findings presented here, and is observed in otherwise non-psychiatric patients recovered from hypercortisolemic states imparted from Cushing’s syndrome [5]. Oligodendrocyte damage is found ubiquitously throughout gray matter and WM regions following prolonged cortisol exposure [6]. Therefore, negative correlations across all tracts in SSD may suggest a reduced physiologic compensatory return of elevated cortisol levels back to baseline that, over time, impacts microstructure health globally. Indeed, in HC, this cortisol response is adequately regulated, and thus positive correlations between all WM tracts and the 40-min cortisol response may indicate physiologically adaptive cortisol effects on WM in response to stress.
WM deficits in SSD are regionally specific [17], although mechanisms of tract-level vulnerabilities remain unclear. Here, prolonged cortisol response was the only stressor-related measure associated with WM FA in SSD, also in a regionally specific pattern, and the two tract-specific patterns were significantly correlated (Fig. 3). These results suggest WM microstructural changes seen broadly in SSD may be related to abnormally prolonged cortisol responses from psychological stress. What is driving the relationship is the strong effects from tracts such as corona radiata, callosal, and fornix fibers that are critical connections between the limbic and the prefrontal functions [23, 54], but weak effects from motor tracts such as PLIC and CST. Similar deficit patterns were found in ENIGMA-schizophrenia samples [17, 53] and linked to core cognitive deficits, particularly processing speed [24]. Therefore, although the nominally significant correlation between processing speed and prolonged cortisol levels was not corrected for multiple comparisons, the fact that processing speed was the only measure among all the clinical measures we had explored that was nominally associated with the prolonged cortisol levels is intriguing. Findings from the current study provide new evidence that the regionally specific SSD FA patterns represented by the ENIGMA study results may in part be a consequence of unresolved stress effects on WM.
Our findings should be considered with several limitations. Antipsychotic medications may impact cortisol secretion; however, current dose was not significantly associated with our findings. Sleep has known effects on rhythmicity of cortisol secretion [60], which was not recorded and will be important to factor in future studies. We also noted that some of the tracts, such as the SFO, were significantly associated with the prolonged cortisol reactivity but did not belong to the list of those suspected to primarily mediate cortisol’s effects, and thus additional research is needed to the underlying mechanism of their associations. Past substance use disorder was unaccounted for in the current sample, which may impart dysfunction on the stress response and WM impairments. The experimental paradigm would benefit from more frequent timepoint gatherings to better understand the cortisol level changes with refinement. This study is cross-sectional, preventing causal interpretability. The positive correlation between psychological stress-induced cortisol levels and WM FA in healthy controls was not hypothesized. Although it is consistent with an adaptive mechanism, this idea requires further experimental testing. FA was the only index used for WM assessment as it is most sensitive to SSD deficits compared to other diffusion parameters [17]. The DTI acquisition protocol was made backward compatible with the legacy imaging protocol executed on the previous version of the Siemens 3T scanner, including weaker gradient and image reconstruction systems, with further details supplied in the methods section. Diurnal variability in neuroendocrine responses may exist despite many controlled aspects of the testing protocol. Roles as a patient vs control in clinical settings, motivational deficits, symptoms like paranoia, or the presence of a research staff may also act as confounders. In our distress task, we minimized staff roles by automating testing sequences.
In conclusion, prolonged cortisol elevations following psychological stressors may provide new perspectives to understand how stress response mechanisms are involved in SSD. Tract-specific SSD WM deficits are important stress-related factors. These observations may guide efforts identifying interventions where targeted management could extend a neuroprotective effect.
Funding and disclosure
Support was received from NIH grants R01MH116948, R01MH112180, P50MH103222, and the University of Maryland/Sheppard Pratt Psychiatry Residency PSTP Program. LEH has received or is planning to receive research funding or consulting fees from Mitsubishi, Your Energy Systems LLC, Neuralstem, Taisho, Heptares, Pfizer, Luye Pharma, Sound Pharma, Takeda, and Regeneron. All other authors declare no financial interests that could represent a conflict of interest.
References
Ritsner M, Maayan R, Gibel A, Strous RD, Modai I, Weizman A. Elevation of the cortisol/dehydroepiandrosterone ratio in schizophrenia patients. Eur Neuropsychopharmacol. 2004;14:267–73.
Schilling EA, Aseltine RH Jr, Gore S. Adverse childhood experiences and mental health in young adults: a longitudinal survey. BMC Public Health. 2007;7:30.
Fisher PA, Van Ryzin MJ, Gunnar MR, Mitigating HPA. axis dysregulation associated with placement changes in foster care. Psychoneuroendocrinology. 2011;36:531–9.
van Andel HW, Jansen LM, Grietens H, Knorth EJ, van der Gaag RJ. Salivary cortisol: a possible biomarker in evaluating stress and effects of interventions in young foster children? Eur Child Adolesc Psychiatry. 2014;23:3–12.
van der Werff SJ, Andela CD, Nienke Pannekoek J, Meijer OC, van Buchem MA, Rombouts SA, et al. Widespread reductions of white matter integrity in patients with long-term remission of Cushing’s disease. Neuroimage Clin. 2014;4:659–67.
Alonso G. Prolonged corticosterone treatment of adult rats inhibits the proliferation of oligodendrocyte progenitors present throughout white and gray matter regions of the brain. Glia. 2000;31:219–31.
Zia MT, Vinukonda G, Vose LR, Bhimavarapu BB, Iacobas S, Pandey NK, et al. Postnatal glucocorticoid-induced hypomyelination, gliosis, and neurologic deficits are dose-dependent, preparation-specific, and reversible. Exp Neurol. 2015;263:200–13.
Miguel-Hidalgo JJ, Carter K, Deloach PH, Sanders L, Pang Y. Glucocorticoid-Induced Reductions of Myelination and Connexin 43 in Mixed Central Nervous System Cell Cultures Are Prevented by Mifepristone. Neuroscience 2019;411:255–69.
Sudheimer KD, O’Hara R, Spiegel D, Powers B, Kraemer HC, Neri E, et al. Cortisol, cytokines, and hippocampal volume interactions in the elderly. Front Aging Neurosci. 2014;6:153.
Xu Y, Fang Z, Wu C, Xu H, Kong J, Huang Q, et al. The long-term effects of adolescent social defeat stress on oligodendrocyte lineage cells and neuroinflammatory mediators in mice. Neuropsychiatr Dis Treat. 2020;16:1321–30.
Woolley CS, Gould E, McEwen BS. Exposure to excess glucocorticoids alters dendritic morphology of adult hippocampal pyramidal neurons. Brain Res. 1990;531:225–31.
Johansson L, Skoog I, Gustafson DR, Olesen PJ, Waern M, Bengtsson C, et al. Midlife psychological distress associated with late-life brain atrophy and white matter lesions: a 32-year population study of women. Psychosom Med. 2012;74:120–5.
Sun Y, Wang Z, Ding W, Wan J, Zhuang Z, Zhang Y, et al. Alterations in white matter microstructure as vulnerability factors and acquired signs of traffic accident-induced PTSD. PLoS One. 2013;8:e83473.
Nugent KL, Chiappelli J, Sampath H, Rowland LM, Thangavelu K, Davis B, et al. Cortisol reactivity to stress and its association with white matter integrity in adults with schizophrenia. Psychosom Med. 2015;77:733–42.
Cheung V, Cheung C, McAlonan GM, Deng Y, Wong JG, Yip L, et al. A diffusion tensor imaging study of structural dysconnectivity in never-medicated, first-episode schizophrenia. Psychol Med. 2008;38:877–85.
Camchong J, Lim KO, Sponheim SR, Macdonald AW. Frontal white matter integrity as an endophenotype for schizophrenia: diffusion tensor imaging in monozygotic twins and patients’ nonpsychotic relatives. Front Hum Neurosci. 2009;3:35.
Kelly S, Jahanshad N, Zalesky A, Kochunov P, Agartz I, Alloza C, et al. Widespread white matter microstructural differences in schizophrenia across 4322 individuals: results from the ENIGMA Schizophrenia DTI Working Group. Mol Psychiatry. 2018;23:1261–69.
Kochunov P, Huang J, Chen S, Li Y, Tan S, Fan F, et al. White matter in schizophrenia treatment resistance. Am J Psychiatry. 2019;176:829–38.
Dedovic K, Duchesne A, Andrews J, Engert V, Pruessner JC. The brain and the stress axis: the neural correlates of cortisol regulation in response to stress. Neuroimage. 2009;47:864–71.
Ulrich-Lai YM, Herman JP. Neural regulation of endocrine and autonomic stress responses. Nat Rev Neurosci. 2009;10:397–409.
Abdul-Rahman MF, Qiu A, Sim K. Regionally specific white matter disruptions of fornix and cingulum in schizophrenia. PLOS ONE. 2011;6:e18652.
Koshiyama D, Fukunaga M, Okada N, Morita K, Nemoto K, Usui K, et al. White matter microstructural alterations across four major psychiatric disorders: mega-analysis study in 2937 individuals. Mol Psychiatry. 2020;25:883–95.
Saunders RC, Aggleton JP. Origin and topography of fibers contributing to the fornix in macaque monkeys. Hippocampus. 2007;17:396–411.
Kochunov P, Coyle TR, Rowland LM, Jahanshad N, Thompson PM, Kelly S, et al. Association of white matter with core cognitive deficits in patients with schizophrenia. JAMA Psychiatry. 2017;74:958–66.
Kirkpatrick B, Strauss GP, Nguyen L, Fischer BA, Daniel DG, Cienfuegos A, et al. The brief negative symptom scale: psychometric properties. Schizophr Bull. 2011;37:300–05.
WechslerDW. Adult Intelligence Scale. 3rd ed. San Antonio, TX: Psychological Corporation; 1997.
Keefe RS, Goldberg TE, Harvey PD, Gold JM, Poe MP, Coughenour L. The Brief Assessment of Cognition in Schizophrenia: reliability, sensitivity, and comparison with a standard neurocognitive battery. Schizophr Res. 2004;68:283–97.
Lejuez C, Kahler C, Brown R. A modified computer version of the Paced Auditory Serial Addition Task (PASAT) as a laboratory-based stressor. Behav Therapist. 2003;26:290–93.
Renna M, Chin S, Seeley S, Fresco D, Heimberg R, Mennin D. The use of the mirror tracing persistence task as a measure of distress tolerance in generalized anxiety disorder. J Ration Emot Cognit Beh Ther. 2017;36:1–9.
Daughters SB, Lejuez CW, Strong DR, Brown RA, Breen RB, Lesieur HR. The relationship among negative affect, distress tolerance, and length of gambling abstinence attempt. J Gambl Stud. 2005;21:363–78.
Chiappelli J, Pocivavsek A, Nugent KL, Notarangelo FM, Kochunov P, Rowland LM, et al. Stress-induced increase in kynurenic acid as a potential biomarker for patients with schizophrenia and distress intolerance. JAMA Psychiatry. 2014;71:761–8.
Kirschbaum C, Wust S, Hellhammer D. Consistent sex differences in cortisol responses to psychological stress. Psychosom Med. 1992;54:648–57.
Cullen AE, Zunszain PA, Dickson H, Roberts RE, Fisher HL, Pariante CM, et al. Cortisol awakening response and diurnal cortisol among children at elevated risk for schizophrenia: relationship to psychosocial stress and cognition. Psychoneuroendocrinology. 2014;46:1–13.
Ljubijankić N, Popović-Javorić R, Sćeta S, Sapcanin A, Tahirović I, Sofić E. Daily fluctuation of cortisol in the saliva and serum of healthy persons. Bosn J Basic Med Sci. 2008;8:110–5.
Griswold MA, Jakob PM, Heidemann RM, Nittka M, Jellus V, Wang J, et al. Generalized autocalibrating partially parallel acquisitions (GRAPPA). Magn Reson Med. 2002;47:1202–10.
Jahanshad N, Kochunov PV, Sprooten E, Mandl RC, Nichols TE, Almasy L, et al. Multi-site genetic analysis of diffusion images and voxelwise heritability analysis: a pilot project of the ENIGMA-DTI working group. NeuroImage. 2013;81:455–69.
Acheson A, Wijtenburg SA, Rowland LM, Winkler A, Mathias CW, Hong LE, et al. Reproducibility of tract-based white matter microstructural measures using the ENIGMA-DTI protocol. Brain Behav. 2017;7:e00615.
McGuire SA, Wijtenburg SA, Sherman PM, Rowland LM, Ryan M, Sladky JH, et al. Reproducibility of quantitative structural and physiological MRI measurements. Brain Behav. 2017;7:e00759.
Savransky A, Chiappelli J, Rowland LM, Wisner K, Shukla DK, Kochunov P, et al. Fornix structural connectivity and allostatic load: empirical evidence from schizophrenia patients and healthy controls. Psychosom Med. 2017;79:770–76.
Wakana S, Jiang H, Nagae-Poetscher LM, van Zijl PC, Mori S. Fiber tract-based atlas of human white matter anatomy. Radiology. 2004;230:77–87.
Smith SM, Jenkinson M, Woolrich MW, Beckmann CF, Behrens TE, Johansen-Berg H, et al. Advances in functional and structural MR image analysis and implementation as FSL. NeuroImage. 2004;23:S208–19.
Kochunov P, Lancaster JL, Glahn DC, Purdy D, Laird AR, Gao F, et al. Retrospective motion correction protocol for high-resolution anatomical MRI. Hum Brain Mapp. 2006;27:957–62.
Walker EF, Diforio D. Schizophrenia: a neural diathesis-stress model. Psychol Rev. 1997;104:667–85.
Kirschbaum C, Pirke KM, Hellhammer DH. The ‘Trier Social Stress Test’-a tool for investigating psychobiological stress responses in a laboratory setting. Neuropsychobiology. 1993;28:76–81.
Brenner K, Liu A, Laplante DP, Lupien S, Pruessner JC, Ciampi A, et al. Cortisol response to a psychosocial stressor in schizophrenia: blunted, delayed, or normal? Psychoneuroendocrinology. 2009;34:859–68.
Jansen LM, Gispen-de Wied CC, Gademan PJ, De Jonge RC, van der Linden JA, Kahn RS. Blunted cortisol response to a psychosocial stressor in schizophrenia. Schizophr Res. 1998;33:87–94.
van Venrooij JA, Fluitman SB, Lijmer JG, Kavelaars A, Heijnen CJ, Westenberg HG, et al. Impaired neuroendocrine and immune response to acute stress in medication-naive patients with a first episode of psychosis. Schizophr Bull. 2012;38:272–9.
Walker EF, Trotman HD, Pearce BD, Addington J, Cadenhead KS, Cornblatt BA, et al. Cortisol levels and risk for psychosis: initial findings from the North American prodrome longitudinal study. Biol Psychiatry. 2013;74:410–7.
Day FL, Valmaggia LR, Mondelli V, Papadopoulos A, Papadopoulos I, Pariante CM, et al. Blunted cortisol awakening response in people at ultra high risk of developing psychosis. Schizophr Res. 2014;158:25–31.
Seitz R, Vracotas N, Bechard-Evans L, King S, Abadi S, Joober R, et al. The Trier Social Stress Test in first episode psychosis patients: Impact of perceived stress, protective factors and childhood trauma. Psychoneuroendocrinology. 2019;105:155–63.
Lange C, Huber CG, Fröhlich D, Borgwardt S, Lang UE, Walter M. Modulation of HPA axis response to social stress in schizophrenia by childhood trauma. Psychoneuroendocrinology. 2017;82:126–32.
Lange C, Deutschenbaur L, Borgwardt S, Lang UE, Walter M, Huber CG. Experimentally induced psychosocial stress in schizophrenia spectrum disorders: a systematic review. Schizophrenia Res. 2017;182:4–12.
van Erp TG, Hibar DP, Rasmussen JM, Glahn DC, Pearlson GD, Andreassen OA, et al. Subcortical brain volume abnormalities in 2028 individuals with schizophrenia and 2540 healthy controls via the ENIGMA consortium. Mol Psychiatry. 2016;21:585.
Poletti CE, Creswell G. Fornix system efferent projections in the squirrel monkey: an experimental degeneration study. J Comp Neurol. 1977;175:101–28.
Penke L, Maniega SM, Murray C, Gow AJ, Valdés Hernández MC, Clayden JD, et al. a general factor of brain white matter integrity predicts information processing speed in healthy older people. J Neurosci. 2010;30:7569–74.
Villalón-Reina JE, Martínez K, Qu X, Ching CRK, Nir TM, Kothapalli D, et al. Altered white matter microstructure in 22q11.2 deletion syndrome: a multisite diffusion tensor imaging study. Mol Psychiatry. 2020;25:2818–31.
Favre P, Pauling M, Stout J, Hozer F, Sarrazin S, Abé C, et al. Widespread white matter microstructural abnormalities in bipolar disorder: evidence from mega- and meta-analyses across 3033 individuals. Neuropsychopharmacology. 2019;44:2285–93.
Diamond DM, Bennett MC, Fleshner M, Rose GM. Inverted-U relationship between the level of peripheral corticosterone and the magnitude of hippocampal primed burst potentiation. Hippocampus. 1992;2:421–30.
Hinds LR, Chun LE, Woodruff ER, Christensen JA, Hartsock MJ, Spencer RL. Dynamic glucocorticoid-dependent regulation of Sgk1 expression in oligodendrocytes of adult male rat brain by acute stress and time of day. PLOS ONE. 2017;12:e0175075.
Leproult R, Copinschi G, Buxton O, Van Cauter E. Sleep loss results in an elevation of cortisol levels the next evening. Sleep. 1997;20:865–70.
Author information
Authors and Affiliations
Contributions
Author contributions included conceptual and study design (ELG, LEH), data collection (JC, PK, ZBM), and interpretation (ELG, JC, MDK, LEH). All authors (ELG, JC, MDK, XD, ZBM, BMA, HO, JS, SL, SV, YM, HB, SC, YT, PK, LEH) participated in drafting, editing, and revising the manuscript.
Corresponding author
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Goldwaser, E.L., Chiappelli, J., Kvarta, M.D. et al. White matter in prolonged glucocorticoid response to psychological stress in schizophrenia. Neuropsychopharmacol. 46, 2312–2319 (2021). https://doi.org/10.1038/s41386-021-01077-4
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41386-021-01077-4
This article is cited by
-
Abnormal cortisol profile during psychosocial stress among patients with schizophrenia in a Chinese population
Scientific Reports (2022)
-
Cortical complexity in world trade center responders with chronic posttraumatic stress disorder
Translational Psychiatry (2021)