Amygdala-medial prefrontal cortex connectivity relates to stress and mental health in early childhood.

Early life stress has been associated with disrupted functional connectivity between the amygdala and medial prefrontal cortex (mPFC), but it is unknown how early in development stress-related differences in amygdala-mPFC connectivity emerge. In a resting-state functional connectivity (rs-FC) analysis with 79 four- to seven-year-old children, we found a significant correlation between more adverse experiences and weaker amygdala-mPFC rs-FC. We also found that weaker amygdala-mPFC rs-FC was associated with higher levels of aggressive behavior and attention problems. These findings suggest that the impact of stress on emotional circuitry is detectable in early childhood and that this impact is associated with mental health difficulties. Connectivity in this circuit may be useful as a marker for mental health risk and for tracking the efficacy of early interventions.

Introduction

Childhood adversity increases the risk for mental health disorders (Shonkoff ; Green ; McLaughlin ). One pathway by which adversity may confer mental health risk is through changes to the amygdala, a brain region central to threat detection and emotional salience (Phelps and LeDoux, 2005), and medial prefrontal cortex (mPFC), a brain region thought to modulate amygdala activity in support of emotion regulation (Quirk and Beer, 2006; Etkin ; McLaughlin ). Disruptions in frontoamygdala functional connectivity have been observed in a wide array of mental health disorders, including depression (Kaiser ), anxiety (Kim ; Hamm ), post-traumatic stress disorder (Gilboa ; Bryant ), and personality and conduct disorders (Marsh ). Thus, characterizing the link between early adversity and atypical development of the frontoamygdala network may provide insight into the emergence of psychopathology.

Adversity has been linked to amygdala volume (increases: Tottenham ; Evans ; decreases: Hanson ), reduced cortical thickness and volume in mPFC (Hanson ; Kelly ), and increased amygdala reactivity during emotion processing and regulation tasks (McCrory ; Dannlowski ; Marusak ; Swartz ). Low socioeconomic status (SES), which elevates the risk of experiencing adversity (Brooks-Gunn and Duncan, 1997), is also associated with amygdala volume (Noble ; Evans ), mPFC structure (Noble ) and amygdala reactivity (Kim ).

In addition to its associations with the structure and function of the amygdala and mPFC, adversity has also been linked to the communication between these regions. Causal experimental designs in animals have shown that stress alters amygdala–PFC interactions (Sánchez ; Milad ). In correlational studies in humans, early life stress has been linked to altered amygdala–mPFC functional connectivity during emotional tasks (Gee ; Javanbakht ; Marusak ; Herringa ). Stress-related changes have also been observed in resting-state functional connectivity (rs-FC). Such changes are of particular interest because rs-FC is not biased by a task—it is thought to be a dynamic measure of the history of co-activation between brain regions (Guerra-Carrillo ; Gabard-Durnam ; Poldrack, 2017). In adults, exposure to childhood trauma is consistently negatively correlated with amygdala–mPFC rs-FC (Herringa ; Birn ; Fan ). Early childhood cortisol levels also show a negative relationship with amygdala–mPFC rs-FC in adulthood (Burghy ). Altered amygdala–mPFC rs-FC has been observed in older children and adolescents, between the ages of 9 and 15, but the directionality of the finding is reversed: two studies, one on negative life events (Pagliaccio ) and one on trauma exposure (Thomason ), observed a positive correlation between stress and amygdala–mPFC rs-FC. It is possible that stress interacts with ongoing developmental changes, leading to inconsistencies in correlation direction.

Amygdala–mPFC rs-FC not only reflects individual differences in early life stress exposure it also mediates the relationship between early life stress and mental health symptoms (Burghy ; Herringa ; Pagliaccio ). These findings suggest that early life stress alters frontoamygdala circuitry, conferring risk for mental health problems (VanTieghem and Tottenham, 2017). However, because studies to date have focused on participants from late childhood through adulthood, it is unclear how early such alterations emerge. Little is known about the effects of stress on amygdala–mPFC functional connectivity in early childhood, a time during which the brain is perhaps most plastic and amenable to intervention (Fox ). Further, because many previous studies have focused on severe cases of adversity and trauma, it is unknown whether more typical variation in early stress exposure relates to differences in amygdala–mPFC circuitry and mental health in early childhood.

In the current study, we investigated whether exposure to stressful life events, such as the death of a family member, parental conflict or a serious accident, is associated with weakened amygdala–mPFC rs-FC in children between the ages of 4 and 7 years old. We also tested whether weakened amygdala–mPFC rs-FC is related to mental health symptoms that reflect early developmental risk for broad adult deficits in emotion regulation (Holtmann ; Bellani ). We focused on the dysregulation profile of the Child Behavior Checklist, which includes the ‘AAA’ symptoms: aggressive behavior, attention problems and anxiety/depression. This profile was first identified by a meta-analysis on the mental health symptoms associated with bipolar disorder in childhood (Mick ), and since has been associated with Axis I disorders more broadly (Biederman ; Mbekou ; Uchida ). Longitudinal studies have found that elevated AAA symptoms in childhood predict impaired psychosocial functioning in adulthood, including mood and anxiety disorders, attention deficit hyperactivity disorder (ADHD), disruptive behavior disorders and substance abuse (Meyer ; Althoff ; Holtmann ). Finally, because socioeconomic status (SES) has been linked to early life stress (Brooks-Gunn and Duncan, 1997), amygdala–mPFC connectivity (Javanbakht ) and mental health symptoms (Lorant ; Reiss, 2013), we recruited a socioeconomically diverse sample and conducted our analyses with and without controlling for SES. To our knowledge, this work represents the first investigation of the links between normative stress exposure, amygdala–mPFC functional connectivity and mental health symptoms in young children.

Methods

Participants

This study was approved by the Committee on the Use of Humans as Experimental Subjects at the Massachusetts Institute of Technology (MIT). All parents provided informed consent, and participants provided verbal assent. Participants were recruited from the Greater Boston Area as part of two larger studies, one study on executive function development, and one study on an intervention (only data prior to the intervention were included). Recruitment for the first study occurred through postings on parent forums, magazine ads, community family events, and Head Start programs. Recruitment for the second study occurred through schools. In an initial screening interview for both studies, parents were asked whether their child had a medical diagnosis of a neurological or psychiatric disorder. For the executive function study, children were excluded if they had a medical diagnosis. The parent of one child reported that their child may have ADHD, but was not taking medication or receiving other treatments, and so the child was included in this study. For the intervention study, children with diagnoses were included in the MRI portion of the study if their diagnoses did not prevent them from participating in the intervention. Two children with Autism Spectrum Disorder participated in this study, but their data were not included for the present analyses.

Resting-state scans were completed for 108 participants. Seventy-nine participants were included in the final sample. Demographic data of this final sample, which was 50.6% female, are summarized in Table 1. Participants were excluded for falling asleep during the scan (1 participant), missing the Life Events Scale for Children (3 participants) or a diagnosis of Autism Spectrum Disorder (2 participants). Participants were also excluded due to average head displacement (3 translations, 3 rotations) of >1 mm (23 participants). The participants who were excluded for motion (n = 23) had slightly, but not significantly, lower maternal education (Exc. median: 16, Inc. median: 16, t = 1.88, P = 0.06), family income (Exc. median: $63 K, Inc. median: $85.5 K, U = 690.0, P = 0.09) and were slightly younger (Exc. median: 6.0, Inc. median: 6.2, t = 1.9, P = 0.06) than included participants (n = 79). Excluded participants did not differ from included participants on exposure to stressful life events (Exc. median: 0, Inc. median: 0, U = 853.5, P = 0.89).

Table 1.

Descriptive statistics of key variables

	n	Mean (s.d.)	Median	Range	Shapiro–Wilk W (P)
Age (years)	79	6.06 (0.96)	6.21	4.05–7.96	0.98 (0.16)
Family income (thousands of dollars)	78	106.6 (78.8)	85.5	5–350	0.93 (<0.001)
Maternal education (years)	78	16.0 (2.8)	16.0	12–20	0.99 (0.92)
Stressful life events (LES-C)	79	2.8 (5.1)	0	0–26	0.71 (<0.001)
Aggressive behavior	38	53.5 (6.7)	50	50–77	0.69 (<0.001)
Attention problems	38	53.1 (5.1)	51	50–67	0.71 (<0.001)
Anxious/depressed	38	52.7 (4.8)	50.5	50–72	0.70 (<0.001)
Somatic complaints	38	53.0 (4.5)	50.0	50–64	0.71 (<0.001)

Note: Results from Shapiro–Wilk tests for normality are shown in the rightmost column, and reveal that only age and maternal education do not significantly deviate from a normal distribution.

Questionnaires

Parents were asked to report total annual income and maternal education. Income data was missing for one child, and maternal education data was missing for one child. Parents completed a modified form of the Life Events Scale for Young Children (LES-C) (Coddington, 1972), which asked parents to report whether specific events had happened to their child within the last 12 months and how stressful the child found these events, from 0 to 4. Examples of items include: ‘Your child had a serious accident or illness’, ‘A family member or close relative died’ and ‘You separated or got divorced from your partner’. LES-C score was calculated by summing the stress ratings for all events. In one of the studies (n = 38), parents also completed the Child Behavior Checklist (CBCL) to identify behavioral and emotional problems (Achenbach and Rescorla, 2000, 2001). Depending on the age of the participant, either the preschool (CBCL/1½–5) or the school-aged (CBCL/6–18) version was administered. Across the two versions of the CBCL, there were four subscales that were consistent: the AAA subscales (Aggressive Behavior, Attention Problems, Anxious/Depressed) and the Somatic Complaints subscale. As described earlier, there is mounting evidence that the dysregulation profile described by the AAA subscales is associated with later psychopathology (Mick ; Meyer ; Althoff ; Holtmann ; Bellani ; Biederman ; Mbekou ; Uchida ). However, because there has also been some research linking Somatic Complaints to anxiety, mood disorders and ADHD (Kasius ; Egger ), all four subscales were analyzed in relation to amygdala connectivity. Descriptive statistics on the LES-C and the four subscales of the CBCL are shown in Table 1. Correlations among these variables and demographic measures are presented in Table 2.

Table 2.

Correlations between demographic and questionnaire measures

	1. Age	2. Family Income	3. Maternal Education	4. LES-C	5. Aggressive Behavior	6. Attention Problems	7. Anxious/ Depressed	8. Somatic Complaints
1. Age (years)	–
2. Family income (thousands of dollars)	0.23*	–
	0.05
	78
3. Maternal education (years)	0.11	0.66***	–
	0.33	<0.001
	78	77
4. Stressful life events (LES-C)	−0.09	−0.38***	−0.27*	–
	0.43	<0.001	0.02
	79	78	78
5. Aggressive behavior	0.12	−0.34*	−0.07	0.29	–
	0.49	0.04	0.66	0.08
	38	37	37	38
6. Attention problems	0.02	−0.25	−0.07	0.10	0.36*	–
	0.92	0.13	0.68	0.57	0.03
	38	37	37	38	38
7. Anxious/depressed	0.09	−0.28	0.02	0.20	0.47**	0.38*	–
	0.61	0.09	0.89	0.24	<0.001	0.02
	38	37	37	38	38	38
8. Somatic complaints	0.19	−0.05	0.04	<0.001	0.22	0.36*	0.43**	–
	0.25	0.77	0.82	0.98	0.18	0.03	0.01
	38	37	37	38	38	38	38

Note: Top values are Spearman rho statistics, middle values are P values, bottom values are ns.

P < 0.05,

P < 0.01,

P < 0.001.

Neuroimaging data acquisition

Imaging was performed at the Athinoula A. Martinos Imaging Center at MIT. Prior to the scanning session, participants were acclimated to the scanning environment with a mock scanner that simulates typical MRI noises. Participants practiced keeping still in the mock scanner, and were given feedback (paused video) if they moved their heads. During the MRI scan, a researcher stayed in the scanner room with the participant to reassure the child and to touch the child’s foot if the child moved.

Scanning was performed using a Siemens MAGNETOM Trio Tim 3 T MRI scanner with a 32 channel coil. A whole-brain, high-resolution, T1-weighted multi-echo structural scan (MPRAGE) was collected (acquisition parameters: TR = 2530 ms, TE = 1.64 ms/3.5 ms/5.36 ms/7.22 ms, flip angle = 7°, voxel size = 1 mm isotropic, matrix size = 192×192, 176 sagittal slices, FOV = 192 mm). This sequence was optimized for participants with high motion (Tisdall ). A 5 minute T2*-weighted gradient echo resting-state scan was also collected (acquisition parameters: TR = 2500 ms, TE = 30 ms, flip angle = 90°, voxel size = 3.0×3.0×3.3 mm, matrix size = 64×64, 41 axial slices, 120 volumes, FOV = 192 mm) with a real-time motion correction technique called Prospective Acquisition Correction (PACE), which adjusts slice orientation and position as the volumes are collected (Thesen ). The first four volumes of each scan were automatically discarded to allow time for scanner magnetization to reach equilibrium. Participants looked at a fixation cross throughout the scan.

Structural analyses

Structural data were used in registration (see below), to identify a bilateral amygdala seed, and to calculate bilateral amygdala volume. The structural images were examined and rated based on quality (1 = best quality to 4 = worst quality, mean = 2.20, s.d. = 0.68) by two researchers who were blind to any other data about the participants. Lower quality structural scans were associated with higher scores on the LES-C (rs(77) = 0.23, P = 0.046, 95% CI [0.03, 0.44]) and on the Aggressive Behavior subscale of the CBCL (rs(36) = 0.32, P = 0.03, 95% CI [0.07, 0.63]). No relations were found between scan quality and the Attention Problems (rs(36) = 0.22, P = 0.19, 95% CI [−0.07, 0.52]), Anxious/Depressed (rs(36) = 0.04, P = 0.81, 95% CI [−0.28, 0.35]) or Somatic Complaints (rs(36) = −0.17, P = 0.30, 95% CI [−0.47, 0.11]) subscales. Scan quality was included as a covariate in structural analyses, and controlled for in whole-brain resting-state analyses.

Cortical reconstruction and volumetric segmentation of the structural images were performed using the standard FreeSurfer processing stream (Dale ). The surfaces were visually inspected for errors and manually edited by researchers who were blind to participant information. Cortical reconstructions were checked for accuracy after editing, and were deemed accurate by blinded researchers for all participants included in the functional connectivity analyses. The amygdala segmentations were also visually inspected for quality, and were found to be accurate for all participants.

Resting-state analyses

The functional imaging data were preprocessed using Nipype, a Python-based framework for flexibly integrating neuroimaging analysis tools (Gorgolewski ). The software packages used in this preprocessing pipeline included FMRIB Software Library (FSL v5.0.8; Jenkinson ), FreeSurfer (v5.3.0; Dale ), Advanced Normalization Tools (ANTs v2.1.0; Avants ) and Nipype’s implementation of Artifact Detection Tools (ART; http://www.nitrc.org/projects/artifact_detect/).

Simultaneous realignment and slice timing correction was conducted using an algorithm implemented in Nipy (Roche, 2011). Outlier volumes in the functional data were defined using ART based on composite motion (>2 mm of head displacement between volumes) and global signal intensity (>3 s.d.s from the mean). Included participants had an average composite motion of <1 mm (mean = 0.35 mm, s.d. = 0.23 mm), and no participants had >15% outlier volumes (mean = 3.23%, s.d. = 3.69%).

The following confounds were regressed out of the functional data: 6 realignment parameters (3 translations, 3 rotations) and their first-order derivatives, outlier volumes flagged by ART (one nuisance regressor per outlier), composite motion, and linear and quadratic polynomials to detrend the data. Five principal components were also derived from segmentations of both cerebrospinal fluid (CSF) and white matter (WM), and regressed from the data, in order to help correct for physiological noise like heart rate and respiration (aCompCor; Behzadi ). FMRIB’s Automated Segmentation Tool (FAST; Zhang ) was used to generate the CSF and WM segmentations from the structural image.

The functional data were bandpass filtered (0.01–0.1 Hz), spatially smoothed with an isotropic 6 mm Gaussian kernel (FWHM), and normalized to the OASIS-30 Atropos template (in MNI152 2 mm space) in a two-step process. First, the median functional image was coregistered to the reconstructed surfaces using FreeSurfer’s bbregister (Greve and Fischl, 2009); next, the structural image was registered to the OASIS-30 Atropos MNI152 template using ANTs. The transformation matrices generated by these two steps were concatenated, allowing images to be transformed directly from functional to MNI space in a single interpolation step.

Bilateral amygdala was selected as the seed, and defined using FreeSurfer’s individually generated subcortical segmentations (Fischl ). The average time series of the bilateral amygdala seed was extracted from the unsmoothed functional data, and FSL’s GLM tool was used to generate whole-brain subject-level connectivity maps. Because confounds related to motion and physiological noise were already filtered out of the functional data during the preprocessing stream, these subject-level GLMs only contained the seed time series as a regressor.

Group-level analyses were conducted with FMRIB’s Local Analysis of Mixed Effects tool (FSL’s FLAME 1), within a mask limited to voxels with unanimous coverage across subjects (which covered the cortex but not inferior cerebellum). The following group-level analyses were conducted: (i) positive and negative group average, (ii) positive and negative correlations with age, (iii) gender differences, (iv) positive and negative correlations with LES-C, and (v) positive and negative correlations with CBCL subscales. All analyses control for number of outliers. We controlled for age and gender unless these variables were the variables of interest. We controlled for dataset (executive function study or intervention study) in analyses 1–4. Analysis 5 only involved data from the executive function study. Residuals of the whole-brain GLMs were normally distributed.

Z-statistic maps were corrected for multiple comparisons with parametric clusterwise inference using FSL’s cluster tool (relies on Gaussian Random Field Theory) at a cluster-defining threshold of z = 3.1 (P < 0.001), neighborhood size of 26, and an FWE-corrected threshold of P < 0.05, based on evidence from Eklund et al. that false positives are not well controlled at a less stringent threshold (Eklund ). Following recommendations for neuroimaging reporting (Poldrack ), uncorrected statistical maps are available on NeuroVault (https://neurovault.org/collections/ZGBHJQOZ/). To determine whether reduced functional connectivity between the amygdala and mPFC is associated with mental health symptoms, parameter estimates were extracted from the whole-brain results and correlated with CBCL scores. All analyses were conducted in Stata 14.1. Confidence intervals were bootstrapped with 1000 repetitions.

Summary of motion considerations

Because young children are more likely to move their heads in the scanner, and because motion is especially problematic for resting-state fMRI (Satterthwaite ; Power ), we adapted our data collection procedures to minimize motion and controlled for motion in our analyses. Here, we summarize these steps. First, children practiced keeping still in a mock scanner while watching a video. The video paused if the child moved >1 mm. The child practiced in the mock scanner for ∼15 minutes. During the scan, a researcher who had previously worked with the child stayed in the scanner room to reassure the child and to remind the child to stay still. Motion correction of the resting-state data was performed both prospectively and retrospectively: prospectively by using real-time motion correction during data collection, and retrospectively by performing simultaneous motion and slice timing correction during preprocessing. Participants were excluded if average composite motion was >1 mm. Motion outlier volumes in the resting-state data were defined based on composite motion >2 mm between volumes. No included participants had >15% outlier volumes. At the single-subject level, we regressed out motion parameters including realignment parameters and their first-order derivatives, outlier volumes and composite motion. Finally, we controlled for number of outlier volumes in all group analyses. We did not observe correlations between number of outliers and age, family income, maternal education, amygdala–mPFC rs-FC parameter estimates or the CBCL measures (P > 0.05).

Results

The amygdala showed widespread positive functional connectivity with mPFC, as well as with other hubs of the limbic and default mode networks (Yeo ), and with subcortical regions including the hippocampus, striatum and thalamus (Figure 1A, n = 79, z = 3.1, P < 0.05). No regions showed significant negative connectivity with the amygdala. Additionally, there were no significant age-related changes or gender differences in amygdala functional connectivity.

Fig. 1.

(A) Amygdala functional connectivity in the group average. Model controls for age, gender, outliers and dataset, and is corrected for multiple comparisons at z = 3.1, P <0.05. N = 79. (B) Negative correlation between amygdala connectivity and stressful life events. Peak voxel coordinates (x, y, z, in MNI space): −10, 58, −6. Maximum z-statistic = 4.75, cluster volume = 699 voxels. Model controls for age, gender, outliers and dataset, and is corrected for multiple comparisons at z = 3.1, P <0.05. N = 79. Scatterplot shows relationship between extracted parameter estimates (adjusted for age, gender, outliers and dataset) and exposure to stressful life events (total score on Life Events Scale for Children). mPFC, medial prefrontal cortex; rs-FC, resting-state functional connectivity.

Greater exposure to stressful life events was significantly correlated with weaker functional connectivity between bilateral amygdala and bilateral mPFC (Figure 1B, n = 79, z = 3.1, P < 0.05). There were no significant positive correlations between LES-C and amygdala functional connectivity. Because greater exposure to stressful life events (LES-C) was related to lower family income (rs(76) = −0.38, P = 0.0007, 95% CI [−0.58, −0.20]), and less strongly to lower maternal education (rs(76) =  −0.27, P = 0.02, 95% CI [−0.48, −0.07]) (Table 2), the analyses were repeated controlling for family income and maternal education, and the results were similar (Supplementary Figure S1). The results were also similar after controlling for structural scan quality (Supplementary Figure S2).

To better understand whether weaker amygdala–mPFC connectivity was related to susceptibility to mental health problems, we extracted parameter estimates from the whole-brain result. These parameter estimates did not significantly deviate from a normal distribution (W = 0.97, P = 0.09). Amygdala–mPFC rs-FC was not correlated with age (r(77) = −0.0005, P = 0.997, 95% CI [−0.24, 0.23]), did not differ by gender (t(77) = 1.38, P = 0.17, 95% CI [−0.23, 0.22]), did not differ between datasets (t(77) = 0.16, P = 0.87, 95% CI [−0.24, 0.23]), and was not related to outliers (rs(77) = −0.05, P = 0.68, 95% CI [−0.30, 0.19]). To be consistent with the whole-brain analyses, analyses of extracted parameter estimates controlled for age, gender and outliers.

Amygdala–mPFC rs-FC was negatively related to Aggressive Behavior (n = 38, z = −2.13, P = 0.03, ϐ = −6.60, 95% CI [−12.67, −0.53], Figure 2A) and to Attention Problems (n = 38, z = −2.17, P = 0.03, ϐ = −4.63, 95% CI [−8.80, −0.46], Figure 2B). No relationships were observed between amygdala–mPFC rs-FC and the Anxious/Depressed subscale (n = 38, z = −0.15, P = 0.89, ϐ = −0.31, 95% CI [−4.36, 3.74]) or the Somatic Complaints subscale (n = 38, z = 1.22, P = 0.22, ϐ = 2.23, 95% CI [−1.36, 5.82]). At the whole-brain level, amygdala functional connectivity with ventral mPFC (orbital PFC) was related to Aggressive Behavior (Supplementary Figure S3). Whole-brain analyses linking amygdala functional connectivity to the other CBCL subscales were non-significant.

Fig. 2.

Negative relationships between amygdala–mPFC rs-FC and mental health symptoms. Scatterplots show parameter estimates extracted from result in Figure 1B, adjusted for age, gender and outliers. (A) Relationship between amygdala–medial prefrontal cortex (mPFC) resting-state functional connectivity (rs-FC) and the Aggressive Behavior subscale of the Child Behavior Checklist (CBCL) (n = 38). (B) Relationship between amygdala–mPFC rs-FC and the Attention Problems subscale of the CBCL (n = 38).

Amygdala–mPFC rs-FC was positively related to family income (n = 78, z = 3.09, P = 0.002, ϐ = 70.46, 95% CI [25.80, 115.13]), and slightly but not significantly related to maternal education (n = 78, z = 1.73, P = 0.08, ϐ = 1.48, 95% CI [−0.20, 3.16]). Controlling for family income and maternal education, the relationships between amygdala–mPFC rs-FC and CBCL subscales were similar in directionality but non-significant (n = 38, Aggressive Behavior: z = −1.44, P = 0.15, ϐ= −4.49, 95% CI [−10.60, 1.62]); Attention Problems: z = −1.89, P = 0.06, ϐ = −3.82, 95% CI [−7.78, 0.15]; Anxious/Depressed: z = 0.53, P = 0.60, ϐ = 1.17, 95% CI [−3.20, 5.55]; Somatic Complaints: z = 3.58, P = 0.06, ϐ = 3.58, 95% CI [−0.20, 7.36]).

Bilateral amygdala volume was not related to LES-C, family income or maternal education (P > 0.2). Amygdala volume additionally was not related to the Aggressive Behavior, Attention Problems, Anxious/Depressed or Somatic Complaints subscales (P > 0.2). Relationships were not changed by including intracranial volume as a covariate. The absence of significant relationships between variables of interest and amygdala volume suggests that differences in amygdala volume are unlikely to explain differences in amygdala functional connectivity. Further, controlling for bilateral amygdala volume in the whole-brain analysis with LES-C did not change the results (Supplementary Figure S4).

Discussion

The present study examined relationships between frontoamygdala circuitry, normative stressful experiences and mental health in early childhood. Exposure to more stressful life events was associated with decreased connectivity between bilateral amygdala and mPFC in young children between the ages of four and seven. Further, decreased connectivity between bilateral amygdala and mPFC was associated with increased levels of aggressive behavior and attention problems. These findings show that the relation between stressful life events and amygdala–mPFC functional connectivity is present early in development, and that this connectivity is associated with symptoms of poorer mental health. Further, whereas related prior findings observed altered frontoamygdala connectivity in older children following trauma, the present study shows that variation in connectivity is also related to variation in normative life stressors.

Weakened functional coupling may be a sign of less effective communication, and perhaps less powerful regulatory control of mPFC over the amygdala. Our results are consistent with previous studies in adults demonstrating that childhood adversity is related to weaker amygdala–mPFC functional coupling at rest (Herringa ; Birn ; Fan ). Our findings complement work in older children and adolescents showing stress-related alterations to amygdala–mPFC rs-FC, although these results showed positive correlations between adversity and amygdala–mPFC rs-FC (Pagliaccio ; Thomason ). The directionality of the link between stress and frontoamygdala connectivity has varied in developmental studies, and the origin of this variability is unclear (Marusak ). It could be due to methodological variability (Marusak ) or could emerge from the interaction between stress and ongoing development. For example, there is evidence that stress is associated with faster maturation of frontoamygdala circuitry (Gee ; Herringa ; Thijssen ). In the 3 year age range of our sample, we did not observe measurable age-related changes in connectivity in this circuit. Our finding of overall positive functional connectivity between amygdala and mPFC was, however, consistent with prior studies (Roy ; Gabard-Durnam ). We did not find evidence for negative amygdala functional connectivity with any regions of the brain. Although previous studies have shown negative connectivity between amygdala and regions like superior frontal gyrus, middle frontal gyrus, and the parietal and occipital lobes, these studies primarily included adolescents and adults (Roy ; Gabard-Durnam ). It may be that robust negative amygdala functional connectivity does not emerge until later in development. For example, a study examining age-related changes in the functional connectivity of mPFC found a shift from positive to negative connectivity with dorsolateral PFC between childhood and adolescence (Chai ).

We also found that weaker amygdala–mPFC functional coupling was associated with more aggressive behavior and more attention problems, behaviors that predict risk for adult psychopathology (Meyer ; Althoff ; Holtmann ; Bellani ; Biederman ; Mbekou ; Uchida ). We did not find a relationship between amygdala–mPFC rs-FC and somatic complaints, and, unlike in studies of older participants (Burghy ; Herringa ; Pagliaccio ), we did not find a relationship with anxiety and depression. It is unclear whether our findings regarding symptoms of anxiety and depression differ from others because amygdala–mPFC connectivity is not related to symptoms of anxiety and depression in young children, or because parental report is less reliable than self-report on internalizing symptoms. Although our results add to the evidence that stress-related changes in amygdala connectivity confer risk of mental health problems, there is also evidence that such changes are protective (Gee ; Herringa ; Thijssen ). It is as yet unclear as to when variation in frontoamygdala connectivity is protective or problematic for mental health.

Consistent with many other studies (Brooks-Gunn and Duncan, 1997), we found a positive correlation between SES, specifically family income, and stressful life events, suggesting that lower SES is associated with greater stress exposure. However, after controlling for SES, the relationship between stressful life events and amygdala–mPFC rs-FC remained similar. This suggests that it is stress, and not other associated features of lower SES, that is related to frontoamygdala circuitry (McLaughlin ). Controlling for SES weakened relationships between amygdala–mPFC rs-FC and mental health symptoms, suggesting that other features of SES may be contributing to mental health. A larger sample size, perhaps with a weaker correlation between SES and adversity, is needed to fully understand how SES interacts with the ‘stress → amygdala–mPFC rs-FC → mental health’ pathway. In addition to changing the risk of stress exposure, SES could change the way in which stress affects the brain (e.g. via the presence or absence of social support) or could impact whether specific neural connectivity patterns give rise to mental health problems. For example, higher SES could be associated with greater compensatory activity in other brain networks.

Amygdala volume was not significantly related to stressful life events, SES, or mental health symptoms, suggesting that amygdala connectivity differences were unlikely to be driven by differences in volume. Previous studies have demonstrated both positive and negative associations between amygdala volume and early life stress (Tottenham ; Lupien ; Hanson ). It may be that volumetric changes are specific to subregions of the amygdala that we could not identify at the resolution of our data.

This study has a number of limitations. First, our measure of stress exposure, the LES-C, is subjective. It requires an accurate memory for stressful events taking place within the previous year, as well as parent assessments of how stressful those events were for the child. Because the scale reflects differences in how children responded to the events, as reported by the parent, it is possible that weaker amygdala–mPFC connectivity confers a predisposition to finding events more stressful. Second, because the relationships between LES-C and the CBCL subscales were not significant in this sample, we were not able to test whether amygdala–mPFC rs-FC mediates the relationship between stress exposure and mental health, and therefore, we cannot rule out other relationships among stress, frontoamygdala circuitry and mental health. Third, our data are cross-sectional, so it cannot address questions about stress-induced changes in developmental trajectory. Longitudinal data from early childhood (or even infancy) into adolescence, with information on the occurrence, duration and severity of stressors, would be useful for further elucidating the mechanisms by which stress exposure shapes mental health. Finally, it is important to note that we cannot infer a causal role of stress from our correlational results. Factors such as genetics and prenatal exposures may covary with stress exposure, amygdala–mPFC rs-FC and mental health.

In sum, we demonstrated that relationships between normative stress exposure and frontoamygdala connectivity emerge in early childhood. Previous studies had conducted brain imaging years after the occurrence of early life stress, an approach with few practical applications for intervention. Our results suggest that abnormal amygdala functional connectivity in young children could be a potential marker of latent risk for poor emotion regulation capacity, and may manifest as clinically relevant symptoms later in life. However, our results also suggest that frontoamygdala connectivity could serve as a potential marker of early intervention efficacy. Early intervention programs, such as the PATHS program for socioemotional learning (Domitrovich ; Fishbein ), have shown promise for improving cognitive and emotional skills in children who have experienced adversity (Raver, 2012). However, there are still considerable individual differences in responsivity to such programs (Bradshaw ). An important future direction is to understand whether neural measures, such as amygdala–mPFC connectivity, enhance the prediction of intervention outcomes in the short- and long-term (Gabrieli ).

Supplementary data

Supplementary data are available at SCAN online.

Conflict of interest. None declared.

Click here for additional data file.