Plasticity of resting state brain networks in recovery from stress

Introduction

When the homeostatic mechanisms are disrupted, namely through prolonged stress exposure, maladaptive responses take place and trigger inappropriate functional responses. It is well-established that prolonged stress has deleterious impact in multiple biological systems, including the central nervous system. In fact, prolonged stress exposure impairs spatial working memory, perceptual attention, behavioral flexibility, and decision making both in rodents and in humans (Joels et al., 2004; Cerqueira et al., 2005; Dias-Ferreira et al., 2009; Soares et al., 2012; Yuen et al., 2012), which has been associated with structural and functional changes of several brain regions. Importantly, some of these maladaptive structural and functional responses to increased chronic stress were shown to be reversible (Sousa et al., 1998; Heine et al., 2004; Cerqueira et al., 2005; Goldwater et al., 2009; Bian et al., 2012; Soares et al., 2012), including evidence showing that as trait positive affect may potentiate recovery and adaptive response (Papousek et al., 2010). However, certain stress effects and specific structural and functional changes may endure after this recovery period (Joels et al., 2004; Gourley et al., 2013). Of notice, most stress recovery studies were performed in rodent models.

A growing field of functional magnetic resonance imaging (fMRI) has provided new insights into the functional connectivity across different brain regions. Indeed, resting state fMRI is being widely used to assess brain regional interactions that comprise the resting state networks (RSNs) (De Luca et al., 2006; Fox and Raichle, 2007), both during resting periods and task-induced deactivations. Moreover, alterations in the normal patterns of RSNs have been associated with several disease states and neuropsychiatric disorders (Zhang and Raichle, 2010; Meda et al., 2012; Sripada et al., 2012), including stress exposure (Soares et al., 2013). Indeed, we previously reported that stressed participants had an hyperactivation pattern of the default mode (DMN), dorsal attention (DAN), ventral attention (VAN), sensorimotor (SMN), and primary visual (VN) networks, paralleled by structural constriction of the DMN brain regions (Soares et al., 2013).

The existence of plastic events in the RSNs after recovery from chronic stress is, however, largely unknown. Indeed, Vaisvaser et al. (2013), using an acute social stress model, examined stress-induced responses in the RSNs and cortisol levels before stress, immediately after the acute stress exposure and 2 h later. The authors found a “recovery” pattern of the DMN connectivity after stress exposure in two of the central hubs of the DMN (seed ROIs at the posterior cingulate cortex and hippocampus), but not in the amygdala-hippocampal disconnectivity that was sustained at 2 h post-stress. Moreover, this increased connectivity was inversely correlated with cortisol levels (Vaisvaser et al., 2013). These results suggest that even acute psychosocial stressors are associated with a prolonged post-stress DMN connectivity response in specific brain regions. This study used only an acute stress model and studies addressing how RSNs respond to chronic stress and identifying specific networks that are associated to an efficient recovery are absent. Therefore, the present study examined the effects of chronic stress on the RSNs following recovery and investigated region-specific changes during successful recovery from chronic stress exposure.

Materials and Methods

Participants, Psychological Tests, and Cortisol Measurements

The participants included in this study were 6 stress participants submitted to prolonged psychological stress exposure (3 males, 3 females; mean age, 23.83 ± 0.37), the same 6 stress participants, 6 weeks after the end of the exposure to stress and 6 controls (3 males, 3 females; mean age, 24.33 ± 1.24). Control participants included a cohort of medical students under their normal academic activities, whereas the stress group included participants that had just finished their long period of preparation for the medical residence selection exam (stress group). Participants responded to a laterality test and to a self-administered questionnaire regarding stress assessment (Perceived Stress Scale—PSS) (Cohen et al., 1983). Participants were further assessed with the Hamilton anxiety scale—HAS (Hamilton, 1959) and the Hamilton depression scale—HDS (Hamilton, 1967) by a certified psychologist. Upon filling of the questionnaires, and immediately before the imaging acquisitions, participants collected saliva samples with the help of Salivette (Sarstedt, Germany) collection devices. Collection took place between 9 and 5 p.m. in all participants. Samples were stored at −20°C until the biologically active, free fraction of the stress hormone cortisol was analyzed using an immunoassay (IBL, Hamburg).

Ethics Statement

The study was conducted in accordance with the principles expressed in the Declaration of Helsinki and was approved by the Ethics Committee of Hospital de Braga (Portugal). The study goals and tests were explained to all participants and all gave informed written consent.

Data Acquisition

Participants were scanned on a clinical approved Siemens Magnetom Avanto 1.5 T (Siemens Medical Solutions, Erlangen, Germany) on Hospital de Braga using the Siemens 12-channel receive-only head coil. The imaging sessions, including one structural T1, one resting state functional, and two task related functional acquisitions, were conducted in the same day and the Siemens Auto Align scout protocol was used to minimize variations in head positioning. For structural analysis, a T1 high-resolution anatomical sequence, 3D MPRAGE (magnetization prepared rapid gradient echo) was performed with the following scan parameters: repetition time (TR) = 2.4 s, echo time (TE) = 3.62 ms, 160 sagittal slices with no gap, field-of-view (FoV) = 234 mm, flip angle (FA) = 8°, in-plane resolution = 1.2 × 1.2 mm² and slice thickness = 1.2 mm. During resting-state fMRI acquisition, using gradient echo T2* weighted echo-planar images (EPIs), participants were instructed to keep the eyes closed and to think about nothing in particular. The imaging parameters were: 100 volumes, TR = 3 s, TE = 50 ms, FA = 90°, in-plane resolution = 3.4 × 3.4 mm², 30 interleaved slices, slice thickness = 5 mm, imaging matrix 64 × 64 and FoV = 220 mm. fMRI paradigm acquisition was acquired using: TR = 2 s, TE = 20 ms, FA = 90°, in-plane resolution and slice thickness 3.3 mm, 38 ascending interleaved axial slices with no gap and FoV = 212 mm. The functional paradigm acquisitions were previously described (Soares et al., 2012) and the paradigm was presented using the fully integrated fMRI system IFIS-SA.

Image Pre-Processing

Before any data processing and analysis, all acquisitions were visually inspected and confirmed that they were not affected by critical head motion and that participants had no brain lesions.

To achieve signal stabilization and allow participants to adjust to the scanner noise, the first 5 resting state fMRI volumes (15 s) were discarded. Data preprocessing was performed using SPM8 (Statistical Parametrical Mapping, version 8, http://www.fil.ion.ucl.ac.uk) analysis software. Images were firstly corrected for slice timing using first slice as reference and SPM8's Fourier phase shift interpolation. To correct for head motion, images were realigned to the mean image with a six-parameter rigid-body spatial transformation and estimation was performed at 0.9 quality, 4 mm separation, 5 mm FWHM smoothing kernel using 2nd degree B-Spline interpolation. No participants exceed head motion higher than 2 mm in translation or 1° in rotation. Images were then spatially normalized to the MNI (Montreal Neurological Institute) standard coordinate system using SPM8 EPI template and trilinear interpolation. Data were then re-sampled to 3 × 3 × 3 mm3 using sinc interpolation, smoothed to decrease spatial noise with a 8 mm, full-width at half-maximum (FWHM), Gaussian kernel, temporally band-pass filtered (0.01–0.08 Hz) and the linear trend was removed. The pre-processing of fMRI paradigm images was previously described (Soares et al., 2012).

Independent Component Analysis and Identification of RSN

Spatial independent component analysis was conducted for using the Group ICA 2.0d of fMRI Toolbox (GIFT, http://www.icatb.sourceforge.net) (Calhoun et al., 2001; Correa et al., 2005). Concisely, spatial ICA analysis is a fully data-driven approach that consists in extracting the non-overlapping spatial maps with temporally coherent time courses that maximize independence. The methodology employed by GIFT can be summarized in three main stages: dimensionality reduction, estimation of the group independent components, and back-reconstruction of each subject's corresponding independent components. The reduction of dimensionality of the functional data and computational load was performed with Principal Component Analysis (PCA) in the concatenated dataset over all subjects, independently of the groups. Then, 20 independent components were estimated, based on a good trade-off (clustering/splitting) between preserving the information in the data while reducing its size (Beckmann et al., 2005; Zuo et al., 2010), using the iterative Infomax algorithm. The ICASSO tool was used to assess the ICA reliability, and 20 computational runs were performed on the dataset, during which the components were being recomputed and compared across runs and the robustness of the results was ensured (Himberg et al., 2004). The previous steps result in the estimation of a mixing matrix with partitions, unique to each subject. The individual independent components were then back-reconstructed from the group-level components. This back-reconstruction step is accomplished by projecting each subject's data onto the inverse of the partition of the calculated matrix corresponding to that subject. The obtained independent components were expressed in t-statistic maps, which were finally converted to a z-statistic. Z-statistic describes the voxels that contributed more intensely to a specific independent component, providing a degree of functional connectivity within the network (Bartels and Zeki, 2005; Beckmann et al., 2005). The final components were visually inspected, sorted, and spatially correlated with resting state functional networks from (Shirer et al., 2012). Each subject's map corresponding to the best-fit component of each RSN was used to perform group statistical analyses.

RSN Deactivation During fMRI Task Analysis

The fMRI decision-making paradigm analyzed to investigate the task-induced deactivations consisted of two different event-related jittered design sessions. First session of valued actions with reward delivery and, after 30 min break, the second session consisted of the devalued actions with the outcome devaluation and extinction. Both sessions had 150 trials, each with 1.5 s for decision, 4 s with the choice highlighted, and 2 s for reward delivery, followed by the interstimulus interval with mean duration of 4 s [please see Soares et al. (2012), for further details].

fMRI paradigm was analyzed by creating a set of regressors at resting and decision making periods, which were convolved with the hemodynamic response function. In order to reliably map task-induced deactivations, we combined all the resting periods (resting baseline condition) and all the decision periods (decision condition), given that decision periods were equally demanding. The contrast used to assess task-induced deactivations was the resting baseline condition minus decision condition. Resulting functional patterns were masked with the previously described RSNs masks (Shirer et al., 2012).

Structural Analysis

Structural analysis based on segmentation of brain structures from T1 high-resolution anatomical data was performed using the freely available Freesurfer toolkit version 5.0 (http://surfer.nmr.mgh.harvard.edu). Intracranial volume (ICV) was used to correct the volumes and the processing pipeline was the same as previously described (Soares et al., 2012). DMN was defined by the summed volume of the angular gyrus of inferior parietal lobe, the posterior cingulate, the precuneus, and the frontopolar region (Raichle et al., 2001; Buckner et al., 2008). The summed volume of the middle frontal gyrus (dorsolateral and prefrontal region) and the posterior parietal region constituted the DAN (Seeley et al., 2007; Sridharan et al., 2008). VAN was constituted by the sum of the temporal-parietal junction and the ventral frontal cortex volumes (Fox et al., 2006). SMN was defined by the summed volume of the paracentral, precentral postcentral, and the cerebellum (Shirer et al., 2012). The summed volume of the cuneus, pericalcarine, and the lingual region constituted the primary VN (Shirer et al., 2012). Auditory network (AN) was defined by the summed volume of the temporal transverse and the temporal superior (Shirer et al., 2012).

Statistical Analyses

Results of the psychological scales, cortisol levels, and regional volumes were analyzed in the IBM SPSS Statistics software, v.19 (IBM, New York). Comparisons between the stress recovered and stress were done with paired samples t-test and between stress recovered and control with two-tailed independent-samples t-test. For all these comparisons significance level was set at 0.05. Values are presented as mean ± standard error of the mean.

Group analysis of the resting state fMRI and task induced deactivations were performed using the second level random effect analyses in SPM8. Initially, within group analyses were performed only to confirm the functional connectivity of the RSNs in the different groups, using one-sample t-tests. Therefore, between group analyses were implemented with directional two-sample t-tests. Functional results for all RSNs were considered significant at p < 0.05 corrected for multiple comparisons using a combination of an uncorrected height threshold of p < 0.025 with a minimum cluster size. The cluster size was determined over 1000 Monte Carlo simulations using AlphaSim program distributed with REST software tool (http://resting-fmri.sourceforge.net/). AlphaSim input parameters were the following: individual voxel probability threshold = 0.025, cluster connection radius = 3 mm, gaussian filter width (FWHM) = 8 mm, number of Monte Carlo simulations = 1000 and mask was set to the corresponding RSN template mask. Anatomical labeling was defined by a combination of visual inspection and Anatomical Automatic Labeling atlas (AAL) (Tzourio-Mazoyer et al., 2002).

Results

Physiological and Behavioral Results

Stress impact was confirmed in several parameters: PSS [Figure 1A; t₍₁₀₎ = 2.52, P < 0.05; Stress Group M = 35.50; SD = 2.59; Control Group M = 30.17; SD = 4.49] and in the HAS [Figure 1A; t₍₁₀₎ = 2.37, P < 0.05 Stress Group M = 11.00; SD = 7.95; Control Group M = 3.00; SD = 2.28], and depression scores [HAD, Figure 1A; t₍₁₀₎ = 3.65, P < 0.01; Stress Group M = 7.50; SD = 2.59; Control Group M = 3.17; SD = 1.33], but only by a trend when it regards to salivary cortisol levels [Figure 1B; t₍₁₀₎ = 1.69, P = 0.12; Stress Group M = 0.44; SD = 0.29; Control Group M = 0.23; SD = 0.10]. After a stress-free period of 6 weeks after the end of the stress exposure, we observed that all the psychological changes were restored [Figure 1A; SPP: t₍₅₎ = 3.72, P < 0.05; Stress Group M = 35.50; SD = 2.59; Stress-Recovered Group M = 30.00; SD = 3.03; anxiety score: t₍₅₎ = 2.86, P < 0.05, Stress Group M = 11.00; SD = 7.95; Stress-Recovered Group M = 2.17; SD = 2.71]; depression score: HAD: [t₍₄₎ = 4.84, P < 0.01; Stress Group M = 7.50; SD = 2.59; Stress-Recovered Group M = 2.5; SD = 0.35], except salivary cortisol levels [t₍₅₎ = 0.67, P = 0.53; Stress Group M = 0.44; SD = 0.29; Stress-Recovered Group M = 1.38; SD = 0.8]. Importantly, stress-recovered group did not differ with the control group in all psychological and salivary cortisol measures [PSS: t₍₁₀₎ = −0.08, P = 0.94; HAS: t₍₁₀₎ = −0.58, P = 0.58; HAD: t₍₁₀₎ = −0.85, P = 0.41; (Figure 1B) Cortisol: t₍₁₀₎ = 1.42, P = 0.19].

FIGURE 1

Figure 1. Clinical characteristics of the cohort. (A) Perceived Stress Scale (PSS), Hamilton Anxiety Scale (HAS); Hamilton Depression Scale (HAD) and (B) Salivary Cortisol levels in the two groups (stress, stress recovered, and control group); *P < 0.05.

Functional Connectivity Results

The ICA analysis revealed the typical spatial pattern of functional connectivity and deactivation in DMN, DAN, VAN, SMN, VN, and AN in all experimental conditions (results not shown).

RSNs in stress and stress—recovered groups

Increased resting functional connectivity was identified in DMN, VAN, and SMN and decreased connectivity in DAN and AN in the stress group when compared to stress recovered participants (Figure 2 and Table 1).

FIGURE 2

Figure 2. The recovery from stress in resting state networks (RSNs) at rest. The images depict areas in which stress participants display greater functional connectivity than stress recovered in the default mode network (DMN) (A), ventral attention network (VAN) (B) and sensorimotor network (SMN) (C) and lower functional connectivity in the dorsal attention network (DAN) (D) and auditory network (AN) (E). Results were extracted by independent component analysis and using paired t-tests, with results considered significant at a corrected for multiple comparisons p < 0.05 threshold.

TABLE 1

Table 1. Group differences (Stress vs. Stress recovered) at rest, in brain regions of the DMN, VAN, SMN, DAN, and AN maps (paired t-tests, corrected for multiple comparisons, p < 0.05).

Regarding DMN, stress group displayed increased functional connectivity mainly in the left cingulum, frontal medial orbitofrontal, right precuneus, and in the left lingual (Table 1). Increased functional connectivity was also found in VAN in stress group in the left parietal inferior and superior, right middle and superior frontal regions (Table 1) whereas in the SMN, increased functional connectivity was found in the left cerebellum (Table 1). In contrast, decreased functional connectivity was found in stress group in DAN, namely in the right parietal inferior, supramarginal, frontal inferior opercularis, and precentral regions (Table 1) as well as in the AN (left superior temporal region) (Table 1).

RSNs in stress—recovered and control groups

Regarding the functional connectivity comparison between stress-recovered and controls, we found that the former presented an increased functional connectivity in the DAN, SMN, and VN. Increased connectivity in the left superior occipital, bilateral superior parietal, right postcentral, left middle and superior frontal, bilateral inferior frontal opercularis and bilateral precentral was found in the DAN of stress-recovered compared to control group (Table 2). A differential pattern of functional connectivity was observed for the VAN that is, while the stress-recovered group presented higher connectivity in the left inferior parietal and bilateral angular, they presented decreased functional connectivity in the bilateral inferior parietal, left angular, bilateral middle frontal and left inferior frontal triangularis. Additionally, stress-recovered group showed decreased connectivity in the DMN in the right anterior cingulate, in the SMN in the bilateral precentral, left paracentral, right postcentral, and bilateral cerebellum and in the VN in the bilateral calcarine (Table 2) when compared to controls (Figure 3).

TABLE 2

Table 2. Group differences (Stress recovered vs. Controls) at rest, in brain regions of the DAN, VAN, SMN, VN, and DMNN maps (two sample t-tests, corrected for multiple comparisons, p < 0.05).

FIGURE 3

Figure 3. Comparison between stress recovered participants and controls in resting state networks (RSNs). The images show areas in which stress recovered participants display greater functional connectivity than controls in the dorsal attention network (DAN) (A), sensorimotor network (SMN) (B), and primary visual network (VN) (C). Lower functional connectivity was found in the default mode network (DMN) (D). Ventral attention network (VAN) (E) displays increased functional connectivity in different regions both in stress recovered (orange) and in control (blue) participants. Results were extracted by independent component analysis and using two-sample t-tests, with results considered significant at a corrected for multiple comparisons p < 0.05 threshold.

Task-induced deactivations in stress and stress—recovered groups

In task-induced deactivations, decreased deactivations in DMN, SMN, and AN were found in stress group when compared to stress-recovered participants (Figure 4 and Table 3). More specifically, decreased deactivations in the left medial frontal orbitofrontal and superior medial frontal were found in DMN of stress group (Table 3). In SMN, stress group presented lower functional deactivation in the left cerebellum (Table 3). The left superior temporal and rolandic operculum in AN were less deactivated in stress group compared to stress-recovered participants (Table 3). No significant region was found to display greater deactivation in stressed participants than in stress recovered in any of the studied RSNs.

FIGURE 4

Figure 4. The recovery from stress in resting state networks (RSNs) during task-induced deactivations. The images illustrate areas of decreased deactivation in stress group when compared to stressed recovered participants in the default mode network (DMN) (A), sensorimotor network (SMN) (B), and auditory network (AN) (C), extracted by general linear model analysis and using paired t-tests, with results considered significant at a corrected for multiple comparisons p < 0.05 threshold. Importantly, no areas of increased deactivation of these RSNs were found in stressed individuals when compared to stress recovered.

TABLE 3

Table 3. Group differences (Stress < Stress recovered) in brain regions of the DMN, SMN, and AD maps in task-induced deactivation (paired t-tests, corrected for multiple comparisons, p < 0.05).

Task-induced deactivations in stress—recovered and control groups

To test for the degree of plasticity in RSNs, we compared deactivation between stress-recovered participants and controls. In this comparison, we found decreased deactivations in DMN, both attention networks, and AN (Figure 5 and Table 4) in stress-recovered group. In DMN, stress-recovered group showed decreased deactivations in the left cuneus, anterior cingulate, right medial frontal orbitofrontal, fusiform and middle temporal and in the left inferior parietal in DAN (Table 4). In VAN, stress-recovered group showed lower deactivation in the left superior parietal and in AN in the bilateral superior temporal (Table 4). No significant region was found to display greater deactivation in stress recovered than in control participants in any of the studied RSNs.

FIGURE 5

Figure 5. Comparison between stress recovered participants and controls in resting state networks (RSNs) during task-induced deactivations. The images demonstrate areas of decreased deactivation in stress-recovered participants when compared to controls in the default mode network (DMN) (A), dorsal attention network (DAN) (B), ventral attention network (VAN) (C), and auditory network (AN) (D), extracted by general linear model analysis and using two-sample t-tests, with results considered significant at a corrected for multiple comparisons p < 0.05 threshold. Importantly, no areas of increased deactivation of these RSNs were found in stress recovered when compared to controls participants.

TABLE 4

Table 4. Group differences (Stress recovered < Controls) in brain regions of the DMN, DAN, VAN, and AN maps in task-induced deactivation (two sample t-tests, corrected for multiple comparisons, p < 0.05).

Expansion/Contraction Maps of the RSNs

Whole brain analysis for relative ICVs did not differ between experimental groups. We showed in a previous study (Soares et al., 2013) that exposure to stress triggered a significant reduction in total DMN volume (corrected for ICV) with specific contraction in the left pCC, and bilateral parietal inferior brain regions. Herein, however, we did not find any significant differences in the volume of any of the RSNs between stress participants before and after recovery from stress. No significant areas of expansion or constriction were found in the dorsal and ventral attention networks, SMN, AN, and primary VN between stress and stress-recovered participants (p = 0.99, p = 0.98, p = 0.87, p = 0.84, and p = 0.99, respectively) and between stress-recovered and control groups (p = 0.89, p = 0.54, p = 0.18, p = 0.47, and p = 0.87, respectively).

Discussion

In this study, we analzsed how the RSNs respond and change following recovery after chronic stress exposure. Our hypothesis was of a continuous recovery effect, in which the connectivity would be decreasing from stress toward the control group. Indeed, we observed a decreased resting functional connectivity in the DMN, VAN, and SMN after stress recovery. Additionally, decreased functional connectivity was also observed in the DAN, SMN, and VN networks in controls, when compared with stress-recovered group. However, only a specific brain region of the DMN (the right anterior cingulate cortex—ACC) showed increased functional connectivity in controls when compared with stress-recovered participants. Results of increased functional connectivity of the DMN at rest after chronic stress exposure are consistent with those previously reported (Soares et al., 2013). More recently, Vaisvaser et al. (2013) evidenced similar results using an acute social stress model.

In the current study, we explored further the plasticity of the RSNs after recovery from the impact of chronic stress-induced changes and showed for the first time that all RSNs, with the exception of the DAN and AN, displayed a functional recovery after the cessation of the exposure to stress. Notably, the comparison with controls allowed us to observe a return to the initial levels the functional connectivity of the DMN, VAN, and AN, but still a sustained pattern of increased functional connectivity of the DAN, SMN, and VN networks.

These results suggest that DAN, SMN, and VN are less plastic when recovering from the impact of stress exposure. The DAN network has been associated with top-down attention processes as inhibitory control, working memory, and response selection. These cognitive processes depend upon the prefrontal integrity (dorsal frontal regions), which are brain regions vulnerable to the effects of stress (Cerqueira et al., 2007). Indeed, animal studies evidenced stress-related prefrontal remodeling (e.g., selective atrophy of the prefrontal cortex, elimination of dendritic spines) after chronic stress exposure (Cerqueira et al., 2007; Gourley et al., 2013). This stress-related prefrontal structural reorganization has been associated with impaired perceptual attention, behavioral flexibility, and decision making in rodents and humans (Cerqueira et al., 2005; Dias-Ferreira et al., 2009; Soares et al., 2012; Yuen et al., 2012). Interestingly, studies analysing the recovery of posttraumatic stress disorder reported that an increased thickness of the dorsolateral prefrontal cortex was associated with greater symptomatic alleviation (Lyoo et al., 2011). The concomitant SMN and VN sustained increased functional connectivity are possibly associated with a motor and visual readiness state that is required for the stress response. This specific pattern of plasticity suggests that some RSNs may be a tool for monitoring effective anti-stress interventions, similar to that proposed to verify the effect of the treatments in several neuropsychiatric diseases (Achard and Bullmore, 2007).

Besides the functional plastic recovery in the connectivity of the RSNs at rest, we also observed a continuum in the pattern of deactivation—that is, there was an increased deactivation from stress toward the control group in all the RSNs. DMN deactivation has been associated with reallocation of attentional resources to cognitively demanding tasks (Hu et al., 2013). Moreover, task-induced RSNs deactivation is correlated with behavioral performance: for example, stronger DMN deactivation in a working memory task predicts better performance (Uddin et al., 2009; Mayer et al., 2010). Increased deactivation observed in our control group and in our stress-recovered participants (comparing with stress). Additionally, abnormal patterns of RSNs deactivation have been associated with several neuropsychiatric diseases (Pomarol-Clotet et al., 2008; Guerrero-Pedraza et al., 2012).

This study shows that while the functional remodeling of RSNs endures, the structural changes (volumetry) of the brain areas involving these networks is still absent after this period of recovery, as no significant areas of expansion or constriction were found in the networks between stress and stress recovered participants; however, in contrast to the difference previously reported in the volumetry of the DMN after stress exposure (Soares et al., 2013) we also did not find significant differences between stress recovered and controls. Difference in results may be related with the limited sample size; the fact that our groups did not differ in physiological cortisol levels and finally, because no direct comparisons were made between stress and control groups.

In summary, the present study contributes to better understand the plastic phenomena that occur in RSNs after the cessation of stress exposure. While we have previously shown the existence of stress-related impairments in the activation-deactivation of RSNs (Soares et al., 2013), here we demonstrate that a functional remodeling of the activation-deactivation pattern of the RSNs takes place following chronic-stress recovery. Although promising, our results should be interpreted with caution mainly due to the reduced size of our sample; therefore, future studies should try to replicate these observations in a larger sample, ideally using exactly the same participants in all conditions, as controls, stressed, and after recovery.

Author Contributions

José M. Soares and Adriana Sampaio contributed in literature search, figures, study design, data collection, data analysis, data interpretation, and writing. Paulo Marques and Luís M. Ferreira as contributed in data collection and data analysis. Nadine C. Santos, Fernanda Marques, Joana A. Palha, João J. Cerqueira, and Nuno Sousa contributed in study design, data interpretation, and writing.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We are thankful to all study participants. José M. Soares, Paulo Marques, and Nadine C. Santos are supported by fellowships of the project SwitchBox-FP7-HEALTH-2010-grant 259772-2; Fernanda Marques is supported by the fellowship SFRH/BPD/33379/2008 funded by the Fundação para a Ciência e Tecnologia (FCT, Portugal). The work was supported by SwitchBox-FP7-HEALTH-2010-grant 259772-2.

References

Achard, S., and Bullmore, E. (2007). Efficiency and cost of economical brain functional networks. PLoS Comput. Biol. 3:e17. doi: 10.1371/journal.pcbi.0030017

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Bartels, A., and Zeki, S. (2005). Brain dynamics during natural viewing conditions–a new guide for mapping connectivity in vivo. Neuroimage 24, 339–349. doi: 10.1016/j.neuroimage.2004.08.044

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Beckmann, C. F., Deluca, M., Devlin, J. T., and Smith, S. M. (2005). Investigations into resting-state connectivity using independent component analysis. Philos. Trans. R. Soc. Lond. B Biol. Sci. 360, 1001–1013. doi: 10.1098/rstb.2005.1634

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Bian, Y., Pan, Z., Hou, Z., Huang, C., Li, W., and Zhao, B. (2012). Learning, memory, and glial cell changes following recovery from chronic unpredictable stress. Brain Res. Bull. 88, 471–476. doi: 10.1016/j.brainresbull.2012.04.008

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Buckner, R. L., Andrews-Hanna, J. R., and Schacter, D. L. (2008). The brain's default network: anatomy, function, and relevance to disease. Ann. N.Y. Acad. Sci. 1124, 1–38. doi: 10.1196/annals.1440.011

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Calhoun, V. D., Adali, T., Pearlson, G. D., and Pekar, J. J. (2001). A method for making group inferences from functional MRI data using independent component analysis. Hum. Brain Mapp. 14, 140–151. doi: 10.1002/hbm.1048

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Cerqueira, J. J., Mailliet, F., Almeida, O. F., Jay, T. M., and Sousa, N. (2007). The prefrontal cortex as a key target of the maladaptive response to stress. J. Neurosci. 27, 2781–2787. doi: 10.1523/JNEUROSCI.4372-06.2007

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Cerqueira, J. J., Pego, J. M., Taipa, R., Bessa, J. M., Almeida, O. F., and Sousa, N. (2005). Morphological correlates of corticosteroid-induced changes in prefrontal cortex-dependent behaviors. J. Neurosci. 25, 7792–7800. doi: 10.1523/JNEUROSCI.1598-05.2005

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Cohen, S., Kamarck, T., and Mermelstein, R. (1983). A global measure of perceived stress. J. Health Soc. Behav. 24, 385–396. doi: 10.2307/2136404

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Correa, N., Adali, T., Li, Y., and Calhoun, V. (2005). Comparison of blind source separation algorithms for FMRI using a new Matlab toolbox: GIFT. Proc. IEEE Int. Conf. Acoust. Speech Signal Process. 5, 401–404. doi: 10.1109/ICASSP.2005.1416325

CrossRef Full Text

De Luca, M., Beckmann, C. F., De Stefano, N., Matthews, P. M., and Smith, S. M. (2006). fMRI resting state networks define distinct modes of long-distance interactions in the human brain. Neuroimage 29, 1359–1367. doi: 10.1016/j.neuroimage.2005.08.035

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Dias-Ferreira, E., Sousa, J. C., Melo, I., Morgado, P., Mesquita, A. R., Cerqueira, J. J., et al. (2009). Chronic stress causes frontostriatal reorganization and affects decision-making. Science 325, 621–625. doi: 10.1126/science.1171203

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Fox, M. D., Corbetta, M., Snyder, A. Z., Vincent, J. L., and Raichle, M. E. (2006). Spontaneous neuronal activity distinguishes human dorsal and ventral attention systems. Proc. Natl. Acad. Sci. U.S.A. 103, 10046–10051. doi: 10.1073/pnas.0604187103

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Fox, M. D., and Raichle, M. E. (2007). Spontaneous fluctuations in brain activity observed with functional magnetic resonance imaging. Nat. Rev. Neurosci. 8, 700–711. doi: 10.1038/nrn2201

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Goldwater, D. S., Pavlides, C., Hunter, R. G., Bloss, E. B., Hof, P. R., McEwen, B. S., et al. (2009). Structural and functional alterations to rat medial prefrontal cortex following chronic restraint stress and recovery. Neuroscience 164, 798–808. doi: 10.1016/j.neuroscience.2009.08.053

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Gourley, S. L., Swanson, A. M., and Koleske, A. J. (2013). Corticosteroid-induced neural remodeling predicts behavioral vulnerability and resilience. J. Neurosci. 33, 3107–3112. doi: 10.1523/JNEUROSCI.2138-12.2013

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Guerrero-Pedraza, A., McKenna, P. J., Gomar, J. J., Sarro, S., Salvador, R., Amann, B., et al. (2012). First-episode psychosis is characterized by failure of deactivation but not by hypo- or hyperfrontality. Psychol. Med. 42, 73–84. doi: 10.1017/S0033291711001073

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hamilton, M. (1959). The assessment of anxiety states by rating. Br. J. Med. Psychol. 32, 50–55. doi: 10.1111/j.2044-8341.1959.tb00467.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hamilton, M. (1967). Development of a rating scale for primary depressive illness. Br. J. Soc. Clin. Psychol. 6, 278–296. doi: 10.1111/j.2044-8260.1967.tb00530.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Heine, V. M., Maslam, S., Zareno, J., Joels, M., and Lucassen, P. J. (2004). Suppressed proliferation and apoptotic changes in the rat dentate gyrus after acute and chronic stress are reversible. Eur. J. Neurosci. 19, 131–144. doi: 10.1046/j.1460-9568.2003.03100.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Himberg, J., Hyvarinen, A., and Esposito, F. (2004). Validating the independent components of neuroimaging time series via clustering and visualization. Neuroimage 22, 1214–1222. doi: 10.1016/j.neuroimage.2004.03.027

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hu, Y., Chen, X., Gu, H., and Yang, Y. (2013). Resting-state glutamate and GABA concentrations predict task-induced deactivation in the default mode network. J. Neurosci. 33, 18566–18573. doi: 10.1523/JNEUROSCI.1973-13.2013

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Joels, M., Karst, H., Alfarez, D., Heine, V. M., Qin, Y., Van Riel, E., et al. (2004). Effects of chronic stress on structure and cell function in rat hippocampus and hypothalamus. Stress 7, 221–231. doi: 10.1080/10253890500070005

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Lyoo, I. K., Kim, J. E., Yoon, S. J., Hwang, J., Bae, S., and Kim, D. J. (2011). The neurobiological role of the dorsolateral prefrontal cortex in recovery from trauma. Longitudinal brain imaging study among survivors of the South Korean subway disaster. Arch. Gen. Psychiatry 68, 701–713. doi: 10.1001/archgenpsychiatry.2011.70

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Mayer, J. S., Roebroeck, A., Maurer, K., and Linden, D. E. (2010). Specialization in the default mode: task-induced brain deactivations dissociate between visual working memory and attention. Hum. Brain Mapp. 31, 126–139. doi: 10.1002/hbm.20850

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Meda, S. A., Gill, A., Stevens, M. C., Lorenzoni, R. P., Glahn, D. C., Calhoun, V. D., et al. (2012). Differences in resting-state functional magnetic resonance imaging functional network connectivity between schizophrenia and psychotic bipolar probands and their unaffected first-degree relatives. Biol. Psychiatry 71, 881–889. doi: 10.1016/j.biopsych.2012.01.025

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Papousek, I., Nauschnegg, K., Paechter, M., Lackner, H. K., Goswami, N., and Schulter, G. (2010). Trait and state positive affect and cardiovascular recovery from experimental academic stress. Biol. Psychol. 83, 108–115. doi: 10.1016/j.biopsycho.2009.11.008

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Pomarol-Clotet, E., Salvador, R., Sarro, S., Gomar, J., Vila, F., Martinez, A., et al. (2008). Failure to deactivate in the prefrontal cortex in schizophrenia: dysfunction of the default mode network? Psychol. Med. 38, 1185–1193. doi: 10.1017/S0033291708003565

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Raichle, M. E., Macleod, A. M., Snyder, A. Z., Powers, W. J., Gusnard, D. A., and Shulman, G. L. (2001). A default mode of brain function. Proc. Natl. Acad. Sci. U.S.A. 98, 676–682. doi: 10.1073/pnas.98.2.676

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Seeley, W. W., Menon, V., Schatzberg, A. F., Keller, J., Glover, G. H., Kenna, H., et al. (2007). Dissociable intrinsic connectivity networks for salience processing and executive control. J. Neurosci. 27, 2349–2356. doi: 10.1523/jneurosci.5587-06.2007

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Shirer, W. R., Ryali, S., Rykhlevskaia, E., Menon, V., and Greicius, M. D. (2012). Decoding subject-driven cognitive states with whole-brain connectivity patterns. Cereb. Cortex 22, 158–165. doi: 10.1093/cercor/bhr099

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Soares, J. M., Sampaio, A., Ferreira, L. M., Santos, N. C., Marques, F., Palha, J. A., et al. (2012). Stress-induced changes in human decision-making are reversible. Transl. Psychiatry 2, e131. doi: 10.1038/tp.2012.59.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Soares, J. M., Sampaio, A., Ferreira, L. M., Santos, N. C., Marques, P., Marques, F., et al. (2013). Stress impact on resting state brain networks. PLoS ONE 8:e66500. doi: 10.1371/journal.pone.0066500

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sousa, N., Almeida, O. F., Holsboer, F., Paula-Barbosa, M. M., and Madeira, M. D. (1998). Maintenance of hippocampal cell numbers in young and aged rats submitted to chronic unpredictable stress. Comparison with the effects of corticosterone treatment. Stress 2, 237–249. doi: 10.3109/10253899809167288

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sridharan, D., Levitin, D. J., and Menon, V. (2008). A critical role for the right fronto-insular cortex in switching between central-executive and default-mode networks. Proc. Natl. Acad. Sci. U.S.A. 105, 12569–12574. doi: 10.1073/pnas.0800005105

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sripada, R. K., King, A. P., Garfinkel, S. N., Wang, X., Sripada, C. S., Welsh, R. C., et al. (2012). Altered resting-state amygdala functional connectivity in men with posttraumatic stress disorder. J. Psychiatry Neurosci. 37, 241–249. doi: 10.1503/jpn.110069

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Tzourio-Mazoyer, N., Landeau, B., Papathanassiou, D., Crivello, F., Etard, O., Delcroix, N., et al. (2002). Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage 15, 273–289. doi: 10.1006/nimg.2001.0978

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Uddin, L. Q., Kelly, A. M., Biswal, B. B., Xavier Castellanos, F., and Milham, M. P. (2009). Functional connectivity of default mode network components: correlation, anticorrelation, and causality. Hum. Brain Mapp. 30, 625–637. doi: 10.1002/hbm.20531

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Vaisvaser, S., Lin, T., Admon, R., Podlipsky, I., Greenman, Y., Stern, N., et al. (2013). Neural traces of stress: cortisol related sustained enhancement of amygdala-hippocampal functional connectivity. Front. Hum. Neurosci. 7:313. doi: 10.3389/fnhum.2013.00313

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Yuen, E. Y., Wei, J., Liu, W., Zhong, P., Li, X., and Yan, Z. (2012). Repeated stress causes cognitive impairment by suppressing glutamate receptor expression and function in prefrontal cortex. Neuron 73, 962–977. doi: 10.1016/j.neuron.2011.12.033

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Zhang, D., and Raichle, M. E. (2010). Disease and the brain's dark energy. Nat. Rev. Neurol. 6, 15–28. doi: 10.1038/nrneurol.2009.198

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Zuo, X. N., Kelly, C., Adelstein, J. S., Klein, D. F., Castellanos, F. X., and Milham, M. P. (2010). Reliable intrinsic connectivity networks: test-retest evaluation using ICA and dual regression approach. Neuroimage 49, 2163–2177. doi: 10.1016/j.neuroimage.2009.10.080

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text