==== Front Neuroimage ClinNeuroimage ClinNeuroImage : Clinical2213-1582Elsevier S2213-1582(18)30234-110.1016/j.nicl.2018.07.022Regular ArticleNeurodevelopmental changes in the relationship between stress perception and prefrontal-amygdala functional circuitry Wu Jingsong a1Geng Xiujuan b1Shao Robin cd1Wong Nichol M.L. cTao Jing aChen Lidian cld@fjtcm.edu.cne⁎⁎⁎Chan Chetwyn C.H. chetwyn.chan@polyu.edu.hkf⁎⁎Lee Tatia M.C. tmclee@hku.hkbcd⁎a Rehabilitation Medicine College, Fujian University of Traditional Chinese Medicine, Fuzhou, Chinab State Key Laboratory of Brain and Cognitive Sciences, The University of Hong Kong, Hong Kongc Laboratory of Neuropsychology, The University of Hong Kong, Hong Kongd Laboratory of Social Cognitive Affective Neuroscience, The University of Hong Kong, Hong Konge Fujian University of Traditional Chinese Medicine, Fuzhou, Chinaf Applied Cognitive Neuroscience Laboratory, Department of Rehabilitation Sciences, The Hong Kong Polytechnic University, Hong Kong⁎ Correspondence to: T.M.C. Lee, Laboratory of Neuropsychology, The University of Hong Kong, Pokfulam Road, Hong Kong. tmclee@hku.hk⁎⁎ Correspondence to: C. C.H. Chan, Applied Cognitive Neuroscience Laboratory, Department of Rehabilitation Sciences, The Hong Kong Polytechnic University, Hung Hom, Hong Kong. chetwyn.chan@polyu.edu.hk⁎⁎⁎ Correspondence to: L. Chen, Fujian University of Traditional Chinese Medicine, No. 1 Huatuo Road Shangjie Minhou, Fuzhou 350122, China. cld@fjtcm.edu.cn1 These authors contributed equally to this work. 23 7 2018 2018 23 7 2018 20 267 274 15 1 2018 27 3 2018 21 7 2018 © 2018 The Authors2018This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).Our brain during distinct developmental phases may show differential responses to perceived psychological stress, yet existing research specifically examining neurodevelopmental changes in stress processing is scarce. To fill in this research gap, this functional magnetic resonance imaging (fMRI) study examined the relationship between perceived stress and resting-state neural connectivity patterns among 67 healthy volunteers belonging to three age groups (adolescents, young adults and adults), who were supposed to be at separate neurodevelopmental phases and exhibit different affect regulatory processes in the brain. While the groups showed no significant difference in self-reported general perceived stress levels, the functional connectivity between amygdala and ventromedial prefrontal cortex (vmPFC) was positively and negatively correlated with perceived stress in adolescents and young adults respectively, while no significant correlations were observed in adults. Furthermore, among adolescents, the causal functional interaction between amygdala and vmPFC exhibited bottom-up connectivity, and that between amygdala and subgenual anterior cingulate cortex exhibited top-down connectivity, both of which changed to bilateral directions, i.e. both bottom-up and top-down connections, in both young adults and adults, supporting the notion that the amygdala and prefrontal cortical circuitries undergo functional reorganizations during brain development. These novel findings have important clinical implications in treating stress-related affective disorders in young individuals. Highlights • Age moderates the relationship between prefrontal-amygdala circuitry and perceived stress. • The VMPFC-amygdala connectivity were distinct in different age groups. • The VMPFC-amygdala connectivity was positively related to stress in adolescents. • The VMPFC-amygdala connectivity was negatively related to stress in young adults. • The Ventral PFC-amygdala connectivity was bi-directional in adults. Keywords Amygdala, brain developmentDynamic causal modelingResting-state functional connectivityStressVentromedial prefrontal cortex ==== Body 1 Introduction Perceived psychological stress is a global, subjective evaluation of the impact of adversities in life that an individual experiences at a given point in time or over a course of time. Such evaluation depends on an individual's appraisal and the resources available for coping with the life adversity experienced (Lazarus and Folkman, 1984). A large body of evidence indicates that high perceived stress predisposes onset of mental illnesses such as depression (Hammen, 2015), anxiety disorders (Glynn et al., 2008), substance use (Rice and Van Arsdale, 2010), and posttraumatic stress disorder (PTSD) (Lagana and Reger, 2009). Notably, high levels of perceived stress and high vulnerabilities to developing affect-related disorders co-occur during adolescence, a period in which individuals undergo rapid and critical psychological and neurobiological developments (Gogtay et al., 2004; Schraml et al., 2011). Yet, existing research specifically investigating the neural correlates of subjective stress during different brain developmental phases is lacking. Such research can advance our understanding of the developmental-dependent neurobiological mechanisms that predispose and precipitate stress-related affective disorders. Stress-related neural processes generally implicate the functioning of reciprocal brain networks involving the limbic and prefrontal cortical (PFC) systems (Lupien et al., 2009; Morey et al., 2012). The amygdala, which is involved in fear processing and emotional memory formation and expression (Becker et al., 2012), is thought to activate the hypothalamus-pituitary-adrenal (HPA) axis in response to threat and challenge. After activation of the HPA axis, the hippocampus and frontal cortex are triggered in order to reduce the HPA axis activation to a homeostatic state (Lupien et al., 2009). Extensive evidence indicates a pivotal role of the amygdala and PFC, especially the ventral medial PFC (vmPFC), for stress-related processing (Gee et al., 2013a, Gee et al., 2013b; McEwen et al., 2016; Tottenham and Galván, 2016). These networks regulate stress-related cognitive, emotional and behavioral processes (Raio et al., 2013; Arnsten et al., 2015). Furthermore, imaging genetic studies suggest that the genotype-related variations in the function of amygdala-vmPFC circuitry implicate a system-level mechanism underlying normal emotional reactivity and genetic susceptibility for affective disorders such as depression, which is closely related to stress (Ressler and Mayberg, 2007; Tost and Meyer-Lindenberg, 2012). The subgenual anterior cingulate cortex that has extensive connections with limbic and paralimbic structures is also considered to play a key role in stress response and emotional regulation (Phillips et al., 2008; Pechtel and Pizzagalli, 2013). Importantly, the functional characteristics of neural circuitries in relate to stress perception and regulation likely vary across different developmental stages. Adolescence is a developmental period that is particularly sensitive to psychosocial adversity, and the affect regulatory networks undergo considerable reorganizations during adolescence and early adulthood (Lupien et al., 2009; Birn et al., 2014; Casey et al., 2008). Generally, limbic structures such as the amygdala and striatum, which are involved in affective reactivity processes, mature earlier than prefrontal areas including the vmPFC and dorsolateral PFC that monitor and regulate the subcortical limbic functions (Steinberg, 2008). Anatomically, the volume of the amygdala starts to stabilize around late adolescence (Wierenga et al., 2014), whereas the ventromedial prefrontal cortical thickness exhibits continuous change throughout adolescence to young adulthood (Shaw et al., 2008). Such asynchronous development of the amygdala and the vmPFC, with the former developing earlier than the latter, renders the amygdala–vmPFC functional coupling less stable in adolescents than in adults, and may cause significant momentary changes in affective processing in the former age group (John and Gross, 2004; Johnson et al., 2016). Also, white-matter fibers connecting amygdala and cortical regions, especially the medial PFC, continue to mature during adolescence until adulthood (Cunningham et al., 2002), which provides the anatomical basis for the functional connectivity between amygdala and vmPFC to continue developing during brain maturation (Gee et al., 2013a, Gee et al., 2013b). Furthermore, while animal studies show bidirectional connectivity between amygdala and orbitofrontal and medial PFC (Barbas et al., 2003), the amygdala to vmPFC projections were found to emerge earlier than the vmPFC to amygdala projections (Bouwmeester et al., 2002), suggesting that the top-down control from medial PFC to subcortical limbic regions may not fully mature until later in development. This cross-sectional neuroimaging study examined the relationship between perceived stress and resting-state functional and effective connectivity in three developmental stages, namely adolescence (high school students), young adults (undergraduate students), and adults. While we acknowledge that considerable individual differences exist within each age group, all individuals in a given age group shared common psychosocial environments (school vs. university vs. society), which aligned with and contributed to the neurodevelopmental profiles of those individuals (Lamm et al., 2014; Simmonds et al., 2014). As such, we considered each of the 3 groups as representing, at least to some extents, distinct and relatively homogeneous psychosocial stages of brain development, and findings on these samples would provide a comprehensive cross-sectional view of the relationships between the maturing brain and stress processing. Following the assumption that prefrontal-limbic networks that are not yet fully matured would be more sensitive and responsive to stress experience (Lupien et al., 2009; Birn et al., 2014), we hypothesized: (1) The association between perceived stress and resting-state functional connectivity between the amygdala and the PFC would differ among the three age groups. Specifically, the adolescent group would exhibit more positive coupling between the PFC-amygdala connectivity strength and levels of perceived stress compared to the other age groups; (2) Based on the assumption that the regulatory function of PFC on amygdala does not fully mature until adulthood, we hypothesized that the causal interactions (i.e. effective connectivity) between the amygdala and the PFC would be different across the developmental stages. Specifically, the PFC-amygdala functional connectivity would be primarily directed from the latter to the former in adolescents, and become bi-directional as people mature in adulthood. 2 Materials and methods 2.1 Participants The study was approved by the Fujian Traditional Chinese Medicine University Institutional Review Board. Prior to the study, written informed consent was obtained from the participant after a detailed description of the study was given. All participants were recruited from the community via advertisement and were screened based on the inclusion criteria. Sixty-seven healthy subjects were included who were assessed as being right-handed by using the Edinburgh Handedness Inventory (Oldfield, 1971), had normal or corrected-to-normal vision and hearing, reported no neurological, neurodegenerative or psychiatric diseases, no history of substance abuse, and were suitable to enter a magnetic resonance image (MRI) scanner. They then completed the Test trait Nonverbal Intelligence, third edition (TONI-III), and the Chinese version of the Perceived Stress Scale (PSS; (Yang and Huang, 2003)). Subjects were divided into three age groups including adolescence (n = 24; 11 females; mean age = 17.02 ± 0.94 years, range = 14.90 to 17.93 years), young adults (n = 22; 9 females; mean age = 19.55 ± 0.43 years, range = 19.00 to 20.61 years), and adults (n = 21; 9 females; mean age = 35.21 ± 4.19 years, range = 30.02 to 45.24 yeas). The three age groups were matched for gender (p = .94) and IQ measured by TONI III (p = .25) (Table 1). While the age difference between the adolescent and young adult groups may be considered small, it is worth noting that i) as stated above, these two groups represent distinct psychosocial and neurodevelopmental phases, and past studies have revealed differences between similarly-aged groups in both regional activations and white-matter structure of affective circuitries (Lamm et al., 2014; Simmonds et al., 2014); ii) existing whole-brain developmental mapping studies show that prefrontal cortical regions undergo considerable changes between late adolescence and early adulthood (Gogtay et al., 2004); and iii) The PSS measure generally lacks validation in children and young adolescent populations.Table 1 Demographic and psychometric characteristics of the participants. Table 1Variables Adolescent Young Adult Adult P value N 24 22 21 Age 17.02 ± 0.94 19.55 ± 0.43 35.21 ± 4.19 <0.0001** Range (14.90–17.93) (19.00–20.61) (30.02–45.24) Gender (m/f) 13/11 13/9 12/9 0.94 IQ(TONI-III) 105.71 ± 14.20 106.91 ± 12.93 111.90 ± 10.92 0.25 Range (84–135) (88–135) (96–135) PSS 25.21 ± 5.87 23.14 ± 6.18 21.24 ± 7.84 0.14 Range (16–37) (13–38) (9–33) Note: IQ, TONI-III; PSS, perceived stress scale. All participants underwent an MRI scan wherein they were instructed to keep their heads still and to remain awake. Anatomical imaging (5 min) and resting state imaging (6 min) were acquired and used in this study. Due to a failure in the acquisition of one subject's functional data, the number of functional images for later analysis was 66 instead of 67. 2.2 The stress measure The PSS is a widely used psychological instrument for measuring nonspecific perceived stress, and it is sensitive and valid for measuring subjective stress levels in both clinical and non-clinical samples (Cohen et al., 1983). In this context, ‘stress’ is defined and measured as situations that are perceived as overwhelming, exhausting, uncontrollable and/or unmanageable. Given the questionnaire primarily focuses on assessing the respondent's stress-related feelings, experiences and emotions, high scores on the PSS reflect greater subjective rather than objective stress experience, which is particularly pertinent to our investigation on the functioning of the neural affect regulatory circuitries. We used the 14-item PSS rated on a 5-point scale from never (0) to very often (4), with higher scores representing elevations in perceived stress. In this study, we focused on detecting the normal perceived stress levels of the participants. Sample items included “Normally, how often have you felt nervous and ‘stressed’?” and “Normally, how often have you found that you could not cope with all the things that you had to do?” The PSS showed a good internal consistency in the current sample with the Cronbach's alpha = 0.812. 2.3 Imaging data acquisition and processing Resting-state fMRI data were acquired on a 3-T GE MRI scanner (GE Medical Systems, Erlangen, USA) with an 8-channel GE head coil. High-resolution T1-weighted anatomical images were acquired using the magnetization-prepared rapid gradient echo (MPRAGE) sequence (TR/TE = 2000 ms /1.75 ms; flip angle = 15°; slice thickness = 1.0 mm; FOV = 240 × 240 × 160 mm3; resolution matrix = 256 × 256; voxel size = 0.94 × 0.94 × 1 mm3) for coregistration with the functional images. Using a T2*-weighted gradient-echo-planar imaging sequence, 180 functional volumes were acquired (slice number/TR/TE/flip angle = 40/2000 ms/30 ms/90°, matrix = 64 × 64, FOV = 225 × 225 × 140 mm3, voxel size = 3.52 × 3.52 × 3.5 mm3). Participants were instructed to stay awake with their eyes open and focus on the cross in the middle of the screen for the entire duration of the resting-state scanning, which lasted for a total of 6 min. 2.4 Functional connectivity estimation Resting-state functional MRI (rs-fMRI) data were preprocessed using Conn (www.nitrc.org/projects/conn), including the following steps: slice-timing correction, motion correction, scrubbing, normalization, nuisance regression, quadratic detrending, band-pass temporal filtering (0.008 < f < 0.09 Hz), and smoothing with an 8-mm Gaussian kernel. Outliers for scrubbing that could be due to movement or disequilibrium effect were defined using the threshold of time intensity and motion (global signal z value <3, and motion mm threshold <0.5 mm). The nuisance regression included the following regressors: 6 motion parameters, the first 5 principle components of the signals in CSF and white matter, and outlier indicators from the scrubbing step. After scrubbing, the average frame displacement in the three age groups were: 0.073, 0.064 and 0.081 mm in adolescents, young adults and adults respectively, and no significant differences were detected using one-way ANOVA (p = .226). The average number of outlier volumes to be scrubbed off in each age group were 8.46, 7.33 and 10.47 mm in adolescents, young adults and adults respectively, and no significant differences were detected using one-way ANOVA (p = .206). The resting-state functional connectivity (rsFC) maps of the amygdala were analyzed using SPM12 and MATLAB (The MathWorks, Inc., Natick, MA, version 2015a). In particular, the right amygdala has been indicated to be crucial in processing negative affective stimuli such as those related to fear (Baker and Kim, 2004; Pegna et al., 2005), and is closely associated with stress processing (Mothersill and Donohoe, 2016). Furthermore, genotypes related to anxiety trait and stress responsiveness, such as the 5-HTTLPR, showed more consistent relations with functions of the right amygdala (Munafo et al., 2008). Therefore the right amygdala was chosen as the seed region-of-interest (ROI). The rsFC of the left amygdala was also computed for sake of completeness. The amygdala ROI was defined following the subcortical delineation in the Harvard-Oxford atlas as provided by FSL (Smith et al., 2004) (see Supplementary Fig. 1 for a depiction of amygdala ROI overlain on the group-level rsFC analysis mask and on a representative participant's unsmoothed image). We included both lateral and mesial amygdala based on existing evidence that 1) the lateral and central amygdala nuclei are tightly interconnected and collectively constitute a coherent affect processing system, with affect-related information passing from the lateral to the mesial portion of amygdala (Maier, 2015); 2) the lateral and central amygdala are respectively involved in acquisition and consolidation of stress-related aversive memory (Haglund et al., 2007); 3) the vmPFC exerts influence on both lateral and mesial amygdala in regulating negative affect and stress (Etkin et al., 2011) and 4) both lateral and mesial amygdala show functional activations to negative affective and stress-related processes (Schiller et al., 2008; Kerr et al., 2012). A functional connectivity map was obtained by computing the correlation between each voxel's time course and the mean time course of the amygdala ROI. 2.5 Statistical analyses on the relationships between age groups, PSS and functional connectivity The effect of age groups on PSS was tested using one-way ANOVA in SPSS v.22 (IBM Corporation, Armonk, NY), controlling for gender and IQ. To test whether different age groups exhibited distinct associations between functional connectivity measures and perceived stress, a second-level full-factorial analysis in SPM12 was performed which tested the following effects: Age as the between-subject factor with three levels (adolescents, young adults and adults); PSS and the interaction of PSS and age as the covariate of interest; and Gender and IQ as nuisance covariates. We controlled for gender and IQ in both behavioral and imaging analyses in light of previous literature indicating potential modulation of stress perception and processing by those factors (e.g. Eisenegger et al., 2011; Wu et al., 2013; Lambert et al., 2014;). F-contrast was constructed to specifically examine the Age by PSS interactive effect. To account for Type 1 errors, we conducted FWE correction with uncorrected p < .001 at voxel level and p < .05 at cluster-level. As we were primarily interested in the connectivity between the amygdala and the PFC, we applied a searching mask defined by the entire PFC by combining 11 cortical areas defined by the Harvard-Oxford atlas as provided by FSL (see Supplementary Fig. 2). For completeness, we also examined the main gender, age, PSS and IQ effects on right amygdala rsFC using F test. Moreover, we also tested possible gender effect on perceived stress (Xu et al., 2015; Burghy et al., 2016) and whether the effect was moderated by different age groups. This was tested in a moderation linear regression model setting age as the moderator on the relationship between gender and PSS, controlling for IQ, using the PROCESS macro (Hayes, 2013) in SPSS v.22 (IBM Corporation, Armonk, NY). The effects were evaluated using bootstrapping (5000 times) as suitable for relatively small samples (Preacher et al., 2007). 2.6 Effective connectivity analyses Causal interactions between volumes of interest (VOIs) at rest were estimated for each individual by spectral dynamic causal modeling (sDCM) using SPM12. The sDCM analysis is widely applied for rsFC studies (Friston et al., 2013). It was conducted in the following steps: selection of the amygdala seed region and the region(s) exhibiting significant interactive effect of PSS by age on connectivity as VOIs; extraction of the first eigenvariate of the time courses in the seed VOI defined by the same mask used in the rsFC analysis, and other VOI(s) defined by a 6 mm sphere centered at the peak coordinates of the significant cluster(s); specification and estimation of the whole model space; and performing a family-level Bayesian model selection (BMS) analysis by grouping subsets of models based on some common feature shared by the models (Penny et al., 2010). We primarily focused on exploring the directionality of the functional couplings between the VOIs, therefore the models were grouped into 1) a family defined by connections from the amygdala seed to the significant clusters, 2) a family defined by connections from the clusters to the amygdala, and 3) a family defined by bidirectional connections between the amygdala seed and the clusters. The family with a posterior probability >95% was regarded as the winning family (Stephan et al., 2010). 3 Results 3.1 PSS across ages One-way ANOVA showed that the three age groups did not differ significantly in general PSS (p = .14). Post-hoc pair-wise t-tests showed no significant PSS differences between any pair of age groups: p = .653 between adolescents and young adults; p = .183 between adolescents and adults; and p = .734 between young adults and adults. 4 rsFC and stress across ages Using the right amygdala as seed, an interaction effect of PSS by age in rsFC was observed in two clusters in the subgenual anterior cingulate cortex (sgACC) with p (FWE corrected) = 0.007 and in the vmPFC with p (FWE corrected) = 0.037 (Fig. 1 and Table 2). To further examine the correlation patterns within the sgACC and vmPFC clusters in each age group, we extracted the average rsFC values in the clusters showing significant interactive effect of PSS by age from each individual subject. The partial correlation analyses were conducted between the extracted rsFC values of each significant cluster and the PSS in each age group, controlling for gender and IQ, using SPSS. Bonferroni correction was used to account for the six comparisons (three age groups and two clusters) (i.e. pcorrected = p*6). In the adolescent group, there was a significant positive correlation between PSS and the extracted rsFC of the amygdala–vmPFC connectivity with r = 0.570, pcorrected = 0.048, and between PSS and the rsFC of the amygdala-sgACC connectivity with r = 0.715, pcorrected < 0.001 (Fig. 2). In young adults, a significant negative correlation was found between PSS and the rsFC of the amygdala–vmPFC circuit, with r = −0.648, pcorrected = 0.018, while the negative correlation between PSS and the rsFC of the amygdala-sgACC circuit did not survive multiple correction (r = −0.546, pcorrected = 0.096). In adults, no significant correlation was observed between PSS and amygdala-vmPFC rsFC(r = 0.310, p = .196), or between PSS and amygdala-sgACC rsFC (r = −0.386, p = .103). No significant cluster was observed for the PSS-by-age interactive effect on the rsFC of the left amygdala (p > .1).Fig. 1 Interaction effect of PSS by age in rsFC with the right amygdala as seed. Both the subgenual cluster (top panel) and the vmPFC cluster (bottom panel) show significant interactive effect. Color bar denotes F statistics. NOTE: rsFC – resting-state functional connectivity; vmPFC – ventral medial prefrontal cortex. Fig. 1Table 2 Clusters showing significant perceived stress by age interactive effect on the resting-state functional connectivity with right amygdala as the seed. Table 2Anatomical regions Cluster Size (# voxels) MNI Coordinates (mm) F-Score Corrected p X y z Subgenual ACC, rectal gyrus, BA25 297 −4 14 −24 20.04 0.007 vmPFC, Rectal Gyrus, BA11 152 −2 48 −26 14.90 0.037 Fig. 2 Scatter plots between the PSS and the rsFC (averaged rsFC over the cluster showing significant interactive effect) in three age groups of adolescents, young adults, and adults. Note: rsFC – resting-state functional connectivity; vmPFC – ventral medial prefrontal cortex; sgACC – subgenual ACC. Fig. 2 We also tested the main effect of age, PSS and IQ on amygdala rsFC, respectively. There was no significant main age or IQ effect on either left or right amygdala rsFC. There was one cluster showing significant main PSS effect on right amygdala rsFC (see Supplementary Fig. 3), but no PSS main effect was observed on left amygdala rsFC (p (FWE) > 0.1). 4.1 Gender effect There was no significant effect of gender on PSS (b = −0.0489, SE = 1.7059, p = .9772), and the age moderation on the association between gender and PSS was not significant (b = 2.8121, SE = 2.1518, p = .1961). The main effect of gender on the rsFC of right amygdala was significant in a cluster localized in the anterior cingulate (X/Y/Z = 2, 28, 32, cluster size = 155 voxels, Fmax = 22.90, p (FWE corrected) = 0.007), with males showing more positive rsFC than females, consistent with existing literature (Gong et al., 2011) (Supplementary Fig. 4). However, the interactive effect of gender and PSS was not significant (p (FWE corrected>0.1). No significant gender or gender-by-PSS effect was observed for the left amygdala rsFC. 4.2 Effective connectivity As there was significant PSS by age interactive effect on the rsFC in the clusters of vmPFC and sgACC, we focused on detecting the directionality of the effective connectivities between the amygdala seed and those two clusters in each of the three age groups. The amygdala seed was defined as one VOI, and two 6-mm spheres were defined as the other two VOIs, centered at the peak coordinate (−4, 14, −24) for the sgACC cluster, and at the peak coordinate (−2, 48, −26) for the vmPFC cluster. The time courses of the three VOIs were extracted. A total of 26–1 = 63 models were estimated for each individual, which were classified into families based on connection directionality. In adolescents, the winning families were defined by effective connectivity from the amygdala to vmPFC with probability = .985, and from the sgACC to amygdala with probability = .956 (Fig. 3a), showing strong evidence for unidirectional connectivity pattern in both cases. In the young adult group, the clearly winning families were bidirectional connectivity between the amygdala and vmPFC (probability = 1), and bidirectional connectivity between the amygdala and sgACC (probability = 1) (Fig. 3b). In the adult group, again, the clearly winning families showed bidirectional connectivity between the amygdala and vmPFC (probability = 1), and bidirectional connectivity between the amygdala and sgACC (probability = 1) (Fig. 3b).Fig. 3 Effective connectivity estimated using spectral dynamic causal modeling analyses in adolescent group (a), young adult and adult groups (b). Fig. 3 5 Discussion By examining the correlations between perceived stress and brain functional networks in adolescents, young adults and adults, we revealed distinct neural functional correlates with perceived stress during different developmental periods. In adolescents, the functional connectivity between right amygdala and both the vmPFC and sgACC correlated positively with perceived stress, with the directionality being from the right amygdala to vmPFC and from the sgACC to right amygdala. In young adults, the right amygdala-vmPFC functional connectivity negatively correlated with perceived stress and showed bidirectional patterns, while no significant correlations between the amygdala connectivity and perceived stress was observed in adults. Based on these results, we propose that the amygdala-vmPFC network continues to mature and engages in different patterns of stress regulatory processes across different developmental periods. The amygdala is the first relay station in the processing of psychosocial stress stimuli, and it continues to develop during adolescence (Ostby et al., 2009). With abundant HPA axis hormone receptors, the amygdala is particularly responsive to stress and may be more sensitive to stress exposures during development than other stress-related regions. Animal studies have shown that the basolateral amygdala exhibits dendrite expansion in response to chronic stress (Boyle, 2013). Neuroimaging studies in humans have also found increased amygdala activity in associate with higher stress levels or greater stress disorder severity (Weber et al., 2013; Savic, 2015). After the amygdala activates the HPA axis in response to adverse stimuli, the medial PFC is then triggered to reduce the HPA axis activation to a homeostatic state (Lupien et al., 2009). Existing neuroimaging evidence indicates both structural (Kim and Richardson, 2009) and functional (Hare et al., 2008) connections between the amygdala and the medial PFC in humans. This circuit not only contributes to various fundamental aspects of emotional behavior including learning and regulation, but is also involved in stress processing and many stress-related disorders at different life stages. For example, weaker connectivity between the amygdala and the vmPFC has been shown to be associated with increased trait anxiety (Kim et al., 2011) and chronic stress (Jovanovic et al., 2011) in the adult population. Furthermore, adolescents with a history of childhood maltreatment (Herringa et al., 2013) or trauma (Pagliaccio et al., 2015) exhibit weaker connectivity between the amygdala and the vmPFC (Nooner et al., 2013). Weaker amygdala and vmPFC connectivity in adolescence has also been shown to mediate the associations between stress and anxiety symptoms (Pagliaccio et al., 2015) and PTSD symptoms (Cisler et al., 2013). Together with our findings of significant correlations between PSS and the amygdala–vmPFC rsFC, the right amygdala–vmPFC connectivity appears to have great implications for stress-related mental health. The connectivity of the amygdala–vmPFC circuit is immature in childhood and continues to develop during adolescence until adulthood (Cunningham et al., 2002; Pattwell et al., 2016). Neuroimaging studies have shown that the amygdala and medial PFC connectivity changes during brain maturation (Gee et al., 2013a, Gee et al., 2013b). The inconsistent functional couplings between the amygdala and the medial PFC during development may be associated with changes in emotional processing (John and Gross, 2004). The functional couplings between amygdala and medial PFC to emotional stimuli has been shown to shift from positive connectivity strength in childhood with average age of 10 years to negative connectivity in adolescence aged around 18 years (Gee et al., 2013a, Gee et al., 2013b). Our results also revealed a shift in the relationships between the perceived stress and the amygdala-vmPFC connectivity during brain maturation, i.e., positive correlation in adolescents and negative correlation in young adults. Furthermore, our effective connectivity analysis showed that the right amygdala-vmPFC connectivity was directed from the former to the latter structure in adolescents, but was bidirectional in young adults and adults. Previous studies revealed that the amygdala has strong bidirectional projections with the PFC (Barbas et al., 2003), and the vmPFC is thought to regulate the activity of the amygdala (Motzkin et al., 2015). In adults, the amygdala and the vmPFC are functionally coupled and work together in stress perception and regulation. Under non-stressed situations, the top-down regulation of the PFC on amygdala provides control over thoughts, emotions and behaviors (Cunningham et al., 2002). Our findings of a negative correlation between PSS and amygdala–vmPFC connectivity in young adults, with a bidirectional connectivity between the amygdala and the vmPFC, further support the top-down regulation model. The lack of correlation between PSS and the amygdala–vmPFC connectivity in well-matured adult brains may suggest that the fully developed circuit in adulthood is more stable, and may thus be less sensitive to stress. On the other hand, the amygdala and the vmPFC have asynchronous developmental trajectories. The amygdala matures relatively early both structurally and functionally (Ulfig et al., 2003; Hare et al., 2008). In contrast, the vmPFC shows protracted structural and functional developments (Gogtay et al., 2004). Tracing studies have also revealed that amygdala-to-vmPFC projections emerge earlier than vmPFC-to-amygdala projections (Bouwmeester et al., 2002), and these connections continue to develop through adolescence in rodents (Cunningham et al., 2002). Neuroimaging studies also reported that effective amygdala–vmPFC connectivity in emotional regulation tasks increases with age in childhood (Perlman and Pelphrey, 2011). Thus, the positive correlation between the PSS and the amygdala–to-vmPFC connectivity observed in adolescents may suggest robustly activated amygdala to stress that provides strong bottom-up affective inputs to the vmPFC, but insufficient top-down control of the yet immature vmPFC on the amygdala. Similar to the amygdala-vmPFC circuitry, the amgydala-sgACC rsFC was also shown to be positively correlated with PSS in adolescents, negatively (although not surviving multiple-test correction) correlated with PSS in young adults, and not correlated in adults. The amygdala-sgACC circuitry showed weaker functional connectivity during emotional processing among the short allele carriers of the 5′ promoter region (5-HTTLPR) of the human serotonin transporter gene, who were shown to possess higher risk for affective disorders (Lotrich and Pollock, 2004), compared to the long allele carriers (Pezawas et al., 2005). In the groups of young adults and adults, the amygdala-sgACC pathway showed bidirectional couplings, similar to the amygdala-vmPFC connectivity. However, the connectivity was directed from sgACC to amygdala in adolescents. Studies have revealed that sgACC is associated with automatic regulation of stress response, whereas regions in ventral and dorsal PFC are involved in voluntary emotional regulation (Phillips et al., 2008). Our findings of the bottom-up projection of amygdala to vmPFC and top-down projection of sgACC-to-amygdala in adolescents not only indicate divergent roles of various areas in medial PFC in stress processing, but also suggest that the automatic emotional regulation starts with a top-down model and develops to a model with bidirectional connections, whereas the voluntary emotional regulation starts with a bottom-up direction and matures to a bidirectional connectivity pattern after adolescence. The current findings should be considered in view of the following limitations. The present study used a cross-sectional design and was thus limited in delineating the complete longitudinal brain developmental trajectories. Future longitudinal studies are needed to provide more definitive evidence on the transitional relationship between perceived stress and prefrontal-amygdala connectivity patterns across different developmental stages. Also, while our imaging results suggested less-developed top-down affect regulatory functions between the vmPFC and amygdala, we obtained no significant differences across age groups on self-reported perceived stress. Two reasons may explain such null finding. First, all our participants were psychologically healthy, thus their levels of stress experience may tend to be more uniform and normalized (Yang and Huang, 2003). Second, the PSS focuses solely on subjective stress experience and provides no objective stress measures. It is thus possible that the adolescent age group generally experienced less objective stress than the older age groups, resulting in comparable levels of perceived stress. Future research could complement the current findings by including clinical samples and employ laboratory stress-manipulation techniques to provide more direct evidence on the implication of vmPFC-amygdala connectivity patterns on affect/stress regulation. 6 Conclusions Our findings support the notion that the neurodevelopmental trajectories from adolescence to adulthood contribute to changes in the relationship between perceived stress and PFC-limbic network functioning. Specifically, the amygdala-vmPFC circuitry was critically involved in regulation of stress processing, the nature and extents of which differed across adolescents, young adults and adults. Findings from this study help to further our understanding of the complex interplay between brain development, stress perception and PFC-limbic network functioning, with clinical implications in understanding and treatment of stress-related affective disorders in young individuals. Appendix A Supplementary data Supplementary material Image 1 Acknowledgements The project was supported by The University of Hong Kong May Endowed Professorship in Neuropsycholog, the KKHO international Charitable Foundation and National Rehabilitation Research Center of Traditional Chinese Medicine, Fujian Provincial Rehabilitation Industrial Institution, and Fujian Rehabilitation Tech Co-Innovation Center (Grant No. X2015003-Collaboration). The funders had no role in study design, data collection, data analysis, data interpretation, or writing of the report. All of the authors had access to the data, and all of the authors agreed to submit the paper for publication. Conflict of interest The authors declare no conflict of interest. Appendix A Supplementary data to this article can be found online at https://doi.org/10.1016/j.nicl.2018.07.022. ==== Refs References Arnsten A.F. Raskind M.A. Taylor F.B. Connor D.F. The effects of stress exposure on prefrontal cortex: translating basic research into successful treatments for post-traumatic stress disorder Neurobiology of stress 1 2015 89 99 25436222 Baker K.B. Kim J.J. Amygdalar lateralization in fear conditioning: evidence for greater involvement of the right amygdala Behav. Neurosci. 118 2004 15 23 14979779 Barbas H. Saha S. Rempel-Clower N. Ghashghaei T. Serial pathways from primate prefrontal cortex to autonomic areas may influence emotional expression BMC Neurosci. 4 2003 25 14536022 Becker B. Mihov Y. Scheele D. Kendrick K.M. Feinstein J.S. Matusch A. Aydin M. Reich H. Urbach H. Oros-Peusquens A.-M. Fear processing and social networking in the absence of a functional amygdala Biol. Psychiatry 72 2012 70 77 22218285 Birn R.M. Patriat R. Phillips M.L. Germain A. Herringa R.J. Childhood maltreatment and combat posttraumatic stress differentially predict fear-related fronto-subcortical connectivity Depress Anxiety 31 2014 880 892 25132653 Bouwmeester H. Smits K. Van Ree J.M. Neonatal development of projections to the basolateral amygdala from prefrontal and thalamic structures in rat J. Comp. Neurol. 450 2002 241 255 12209853 Boyle L.M. A neuroplasticity hypothesis of chronic stress in the BasolateralAmygdala Yale J. Biol. & Med. 86 2013 117 125 23766733 Burghy C.A. Fox M.E. Cornejo M.D. Experience-driven differences in childhood cortisol predict affect-relevant brain function and coping in adolescent monozygotic twins Sci. Rep. 6 2016 37081 27872489 Casey B.J. Getz S. Galvan A. The adolescent brain Developmental Rev. 28 2008 62 77 Cisler J.M. Steele J.S. Smitherman S. Lenow J.K. Kilts C.D. Neural processing correlates of assaultive violence exposure and PTSD symptoms during implicit threat processing: a network-level analysis among adolescent girls Psychiatry Res. Neuroimaging 214 2013 238 246 Cohen S. Kamarck T. Mermelstein R. A global measure of perceived stress Journal of health and social behavior 1983 385 396 6668417 Cunningham M.G. Bhattacharyya S. Benes F.M. Amygdalo-cortical sprouting continues into early adulthood: implications for the development of normal and abnormal function during adolescence J. Comp. Neurol. 453 2002 116 130 12373778 Eisenegger C. Haushofer J. Fehr E. The role of testosterone in social interaction Trends Cogn. Sci. 15 6 2011 263 271 21616702 Etkin A. Egner T. Kalisch R. Emotional processing in anterior cingulate and medial prefrontal cortex Trend Cogn Sci 15 2 2011 85 93 Friston K.J. Kahan J. Biswal B. Razi A. A DCM for resting state fMRI NeuroImage 94 2013 396 407 24345387 Gee D.G. Gabard-Durnam L.J. Flannery J. Goff B. Humphreys K.L. Telzer E.H. Hare T.A. Bookheimer S.Y. Tottenham N. Early developmental emergence of human amygdala–prefrontal connectivity after maternal deprivation Proc. Natl. Acad. Sci. 110 2013 15638 15643 24019460 Gee D.G. Humphreys K.L. Flannery J. Goff B. Telzer E.H. Shapiro M. Hare T.A. Bookheimer S.Y. Tottenham N. A developmental shift from positive to negative connectivity in human amygdala–prefrontal circuitry J. Neurosci. 33 2013 4584 4593 23467374 Glynn L.M. Schetter C.D. Hobel C.J. Sandman C.A. Pattern of perceived stress and anxiety in pregnancy predicts preterm birth Health Psychol. 27 2008 43 18230013 Gogtay N. Giedd J.N. Lusk L. Hayashi K.M. Greenstein D. Vaituzis A.C. Nugent T.F. Herman D.H. Clasen L.S. Toga A.W. Dynamic mapping of human cortical development during childhood through early adulthood Proc. Natl. Acad. Sci. U. S. A. 101 21 2004 8174 8179 15148381 Gong G. He Y. Evans A.C. Brain connectivity: gender makes a difference Neuroscientist 17 5 2011 Oct 575 591 21527724 Haglund M.E. Nestadt P.S. Cooper N.S. Southwick S.M. Charney D.S. Psychobiological mechanisms of resilience: relevance to prevention and treatment of stress-related psychopathology Dev. Psychopathol. 19 3 2007 889 920 17705907 Hammen C.L. Stress and depression: old questions, new approaches Curr. Opin. Psychol. 4 2015 80 85 Hare T.A. Tottenham N. Galvan A. Voss H.U. Glover G.H. Casey B.J. Biological substrates of emotional reactivity and regulation in adolescence during an emotional go-nogo task Biol. Psychiatry 63 2008 927 934 18452757 Hayes A.F. Introduction to Mediation, Moderation, and Conditional Process Analysis: A Regression-Based Approach 2013 Guilford Press Herringa R.J. Birn R.M. Ruttle P.L. Burghy C.A. Stodola D.E. Davidson R.J. Essex M.J. Childhood maltreatment is associated with altered fear circuitry and increased internalizing symptoms by late adolescence Proc. Natl. Acad. Sci. 110 2013 19119 19124 24191026 John O.P. Gross J.J. Healthy and unhealthy emotion regulation: personality processes, individual differences, and life span development J. Pers. 72 2004 1301 1334 15509284 Johnson C.M. Loucks F.A. Peckler H. Thomas A.W. Janak P.H. Wilbrecht L. Long-range orbitofrontal and amygdala axons show divergent patterns of maturation in the frontal cortex across adolescence Develop. Cognitive Neurosci. 18 2016 113 120 Jovanovic H. Perski A. Berglund H. Savic I. Chronic stress is linked to 5-HT 1A receptor changes and functional disintegration of the limbic networks NeuroImage 55 2011 1178 1188 21211567 Kerr D.L. McLaren D.G. Mathy R.M. Nitschke J.B. Controllability modulates the anticipatory response in the human ventromedial prefrontal cortex Front. Psychol. 3 2012 Dec 557 23550176 Kim J.H. Richardson R. The effect of the μ-opioid receptor antagonist naloxone on extinction of conditioned fear in the developing rat Learn. Mem. 16 2009 161 166 19223599 Kim M.J. Gee D.G. Loucks R.A. Davis F.C. Whalen P.J. Anxiety dissociates dorsal and ventral medial prefrontal cortex functional connectivity with the amygdala at rest Cereb. Cortex 21 2011 1667 1673 21127016 Lagana L. Reger S.L. A pilot study on perceived stress and PTSD symptomatology in relation to four dimensions of older women's physical health Aging and Mental Health 13 2009 885 893 19888708 Lambert K.G. Hyer M.M. Rzucidlo A.A. Bergeron T. Landis T. Bardi M. Contingency-based emotional resilience: effort-based reward training and flexible coping lead to adaptive responses to uncertainty in male rats Front Behav. Neurosci. 2014 8 24478655 Lamm C. Benson B.E. Guyer A.E. Perez-Edgar K. Fox N.A. Pine D.S. Ernst M. Longitudinal study of striatal activation to reward and loss anticipation from mid-adolescence into late adolescence/early adulthood Brain Cogn. 89 2014 51 60 24485273 Lazarus R.S. Folkman S. Coping and adaptation The Handbook of Behavioral Medicinem 1984 282 325 Lotrich F.E. Pollock B.G. Meta-analysis of serotonin transporter polymorphisms and affective disorders Psychiatr. Genet. 14 2004 121 129 15318024 Lupien S.J. McEwen B.S. Gunnar M.R. Heim C. Effects of stress throughout the lifespan on the brain, behaviour and cognition Nat. Rev. Neurosci. 10 2009 434 445 19401723 Maier S.F. Behavioral control blunts reactions to contemporaneous and future adverse events: medial prefrontal cortex plasticity and a corticostriatal network Neurobiol. Stress 1 2015 12 22 25506602 McEwen B.S. Nasca C. Gray J.D. Stress effects on neuronal structure: Hippocampus, amygdala, and prefrontal cortex Neuropsychopharmacology 41 2016 3 23 26076834 Morey R.A. Gold A.L. Labar K.S. Beall S.K. Brown V.M. Haswell C.C. Nasser J.D. Wagner H.R. McCarthy G. Amygdala volume changes in posttraumatic stress disorder in a large case-controlled veterans group Arch. Gen. Psychiatry 69 2012 1169 1178 23117638 Mothersill O. Donohoe G. Neural effects of social environmental stress - an activation likelihood estimation meta-analysis Psychol. Med. 46 2016 2015 2023 27216635 Motzkin J.C. Philippi C.L. Oler J.A. Kalin N.H. Baskaya M.K. Koenigs M. Ventromedial prefrontal cortex damage alters resting blood flow to the bed nucleus of stria terminalis Cortex; a journal devoted to the study of the nervous system and behavior 64 2015 281 288 25569763 Munafo M.R. Brown S.M. Hariri A.R. Serotonin transporter (5-HTTLPR) genotype and amygdala activation: a meta-analysis Biol. Psychiatry 63 2008 852 857 17949693 Nooner K.B. Mennes M. Brown S. Castellanos F.X. Leventhal B. Milham M.P. Colcombe S.J. Relationship of trauma symptoms to amygdala-based functional brain changes in adolescents J. Trauma. Stress. 26 2013 784 787 24343754 Oldfield R.C. The assessment and analysis of handedness: the Edinburgh inventory Neuropsychologia 9 1971 97 113 5146491 Ostby Y. Tamnes C.K. Fjell A.M. Westlye L.T. Due-Tonnessen P. Walhovd K.B. Heterogeneity in subcortical brain development: a structural magnetic resonance imaging study of brain maturation from 8 to 30 years J. Neurosci. 29 2009 11772 11782 19776264 Pagliaccio D. Luby J.L. Bogdan R. Agrawal A. Gaffrey M.S. Belden A.C. Botteron K.N. Harms M.P. Barch D.M. Amygdala functional connectivity, HPA axis genetic variation, and life stress in children and relations to anxiety and emotion regulation J. Abnorm. Psychol. 124 2015 817 26595470 Pattwell S.S. Liston C. Jing D. Ninan I. Yang R.R. Witztum J. Murdock M.H. Dincheva I. Bath K.G. Casey B.J. Dynamic changes in neural circuitry during adolescence are associated with persistent attenuation of fear memories Nature communications 2016 7 Pechtel P. Pizzagalli D.A. Disrupted reinforcement learning and maladaptive behavior in women with a history of childhood sexual abuse: a high-density event-related potential study JAMA psychiatry 70 2013 499 507 23487253 Pegna A.J. Khateb A. Lazeyras F. Seghier M.L. Discriminating emotional faces without primary visual cortices involves the right amygdala Nat. Neurosci. 8 2005 24 25 15592466 Penny W.D. Stephan K.E. Daunizeau J. Rosa M.J. Friston K.J. Schofield T.M. Leff A.P. Comparing families of dynamic causal models PLoS Comput. Biol. 6 2010 e1000709 Perlman S.B. Pelphrey K.A. Developing connections for affective regulation: age-related changes in emotional brain connectivity J. Exp. Child Psychol. 108 2011 607 620 20971474 Pezawas L. Meyer-Lindenberg A. Drabant E.M. Verchinski B.A. Munoz K.E. Kolachana B.S. Egan M.F. Mattay V.S. Hariri A.R. Weinberger D.R. 5-HTTLPR polymorphism impacts human cingulate-amygdala interactions: a genetic susceptibility mechanism for depression Nat. Neurosci. 8 2005 828 834 15880108 Phillips M.L. Ladouceur C.D. Drevets W.C. A neural model of voluntary and automatic emotion regulation: implications for understanding the pathophysiology and neurodevelopment of bipolar disorder Mol. Psychiatry 13 829 2008 833 857 Preacher K.J. Rucker D.D. Hayes A.F. Addressing moderated mediation hypotheses: theory, methods, and prescriptions Multivar. Behav. Res. 42 2007 185 227 Raio C.M. Orederu T.A. Palazzolo L. Shurick A.A. Phelps E.A. Cognitive emotion regulation fails the stress test Proc. Natl. Acad. Sci. U. S. A. 110 2013 15139 15144 23980142 Ressler K.J. Mayberg H.S. Targeting abnormal neural circuits in mood and anxiety disorders: from the laboratory to the clinic Nat. Neurosci. 10 2007 1116 1124 17726478 Rice K.G. Van Arsdale A.C. Perfectionism, perceived stress, drinking to cope, and alcohol-related problems among college students J. Couns. Psychol. 57 2010 439 Savic I. Structural changes of the brain in relation to occupational stress Cereb. Cortex 25 2015 Schiller D. Levy I. Niv Y. Ledoux J.E. Phelps E.A. From fear to safety and back: reversal of fear in the human brain J. Neurosci. 28 45 2008 11517 11525 18987188 Schraml K. Perski A. Grossi G. Simonsson-Sarnecki M. Stress symptoms among adolescents: the role of subjective psychosocial conditions, lifestyle, and self-esteem J. Adolesc. 34 5 2011 Oct 1 987 996 21147499 Shaw P. Kabani N.J. Lerch J.P. Eckstrand K. Lenroot R. Gogtay N. Greenstein D. Clasen L. Evans A. Rapoport J.L. Giedd J.N. Wise S.P. Neurodevelopmental trajectories of the human cerebral cortex J. Neurosci. 28 2008 3586 3594 18385317 Simmonds D.J. Hallquist M.N. Asato M. Luna B. Developmental stages and sex differences of white matter and behavioral development through adolescence: a longitudinal diffusion tensor imaging (DTI) study NeuroImage 92 2014 356 368 24384150 Smith S.M. Jenkinson M. Woolrich M.W. Beckmann C.F. Behrens T.E. Johansen-Berg H. Bannister P.R. De Luca M. Drobnjak I. Flitney D.E. Niazy R.K. Saunders J. Vickers J. Zhang Y. De Stefano N. Brady J.M. Matthews P.M. Advances in functional and structural MR image analysis and implementation as FSL NeuroImage 23 Suppl. 1 2004 S208 S219 15501092 Steinberg L. A social neuroscience perspective on adolescent risk-taking Developmental Rev. 28 2008 78 106 Stephan K.E. Penny W.D. Moran R.J. den Ouden H.E. Daunizeau J. Friston K.J. Ten simple rules for dynamic causal modeling NeuroImage 49 2010 3099 3109 19914382 Tost H. Meyer-Lindenberg A. Puzzling over schizophrenia: schizophrenia, social environment and the brain Nat. Med. 18 2012 211 213 22310688 Tottenham N. Galván A. Stress and the adolescent brain: amygdala-prefrontal cortex circuitry and ventral striatum as developmental targets Neurosci. Biobehav. Rev. 70 2016 217 227 27473936 Ulfig N. Setzer M. Bohl J. Ontogeny of the human amygdala Ann. N. Y. Acad. Sci. 985 2003 22 33 12724145 Weber M. Killgore W.D. Rosso I.M. Britton J.C. Schwab Z.J. Weiner M.R. Simon N.M. Pollack M.H. Rauch S.L. Voxel-based morphometric gray matter correlates of posttraumatic stress disorder J. Anxiety Disorders 27 2013 413 419 Wierenga L. Langen M. Ambrosino S. Van D.S. Oranje B. Durston S. Typical development of basal ganglia, hippocampus, amygdala and cerebellum from age 7 to 24 NeuroImage 96 2014 67 72 24705201 Wu G. Feder A. Cohen H. Kim J.J. Calderon S. Charney D.S. Mathé A.A. Understanding resilience Front. Behav. Neurosci. 7 2013 Xu X. Bao H. Strait K. Spertus J.A. Lichtman J.H. D’onofrio G. … Bueno H. Sex differences in perceived stress and early recovery in young and middle-aged patients with acute myocardial infarction Circulation 131 7 2015 614 25679303 Yang T.Z. Huang H.T. An epidemiological study on stress among urban residents in social transition period Chin. J. Epidemiol. 24 2003 760 764