Xinyuan Li1,2,3,4
Muhammad Naveed Iqbal Qureshi3,4
David P. Laplante5
Guillaume Elgbeili2
Sherri Lee Jones6
Suzanne King2,6*
Pedro Rosa-Neto3,4,7- 1Integrated Program in Neuroscience, McGill University, Montreal, QC, Canada
- 2Mental Health and Society Division, Douglas Mental Health University Institute, Montreal, QC, Canada
- 3Translational Neuroimaging Laboratory, McGill University Research Centre for Studies in Aging, Montreal, QC, Canada
- 4Montreal Neurological Institute, McGill University, Montreal, QC, Canada
- 5Centre for Child Development and Mental Health, Lady Davis Institute-Jewish General Hospital, Montreal, QC, Canada
- 6Department of Psychiatry, McGill University, Montreal, QC, Canada
- 7Department of Neurology and Neurosurgery, McGill University, Montreal, QC, Canada
Background: Studies have shown that prenatal maternal stress alters volumes of the amygdala and hippocampus, and alters functional connectivity between the amygdala and prefrontal cortex. However, it remains unclear whether prenatal maternal stress (PNMS) affects volumes and functional connectivity of these structures at their subdivision levels.
Methods: T1-weighted MRI and resting-state functional MRI were obtained from 19-year-old young adult offspring with (n = 39, 18 male) and without (n = 65, 30 male) exposure to PNMS deriving from the 1998 ice storm. Volumes of amygdala nuclei, hippocampal subfields and prefrontal subregions were computed, and seed-to-seed functional connectivity analyses were conducted.
Results: Compared to controls, young adult offspring exposed to disaster-related PNMS had larger volumes of bilateral whole amygdala, driven by the lateral, basal, central, medial, cortical, accessory basal nuclei, and corticoamygdaloid transition; larger volumes of bilateral whole hippocampus, driven by the CA1, HATA, molecular layer, fissure, tail, CA3, CA4, and DG; and larger volume of the prefrontal cortex, driven by the left superior frontal. Inversely, young adult offspring exposed to disaster-related PNMS had lower functional connectivity between the whole amygdala and the prefrontal cortex (driven by bilateral frontal poles, the left superior frontal and left caudal middle frontal); and lower functional connectivity between the hippocampal tail and the prefrontal cortex (driven by the left lateral orbitofrontal).
Conclusion: These results suggest the possibility that effects of disaster-related PNMS on structure and function of subdivisions of offspring amygdala, hippocampus and prefrontal cortex could persist into young adulthood.
1. Introduction
The fetal programming hypothesis proposes that fetal exposure to an adverse intrauterine environment induces long-lasting changes in the offspring (Barker, 2003). Studies have shown that pre-natal maternal depression (Li et al., 2022) and disaster-related prenatal maternal stress (PNMS) (Lafortune et al., 2021) are associated with offspring cognitive, socio-emotional and behavioral development that may persist into young adulthood. This hypothesis also applies to brain development (Lautarescu et al., 2020). A recent review suggests that pre-natal maternal depression, anxiety and stressful life events are associated with atypical volumes and functional connectivity of widespread brain regions in the offspring (Lautarescu et al., 2020).
Prenatal maternal stress has been associated with heightened stress reactivity in the offspring (Yong Ping et al., 2015), which may be explained by alterations in the development of the fetal brain. In the brain, the amygdala, hippocampus and prefrontal cortex contain a very high density of glucocorticoid receptors, and are thus highly sensitive to elevated levels of glucocorticoids (Herman et al., 2005). The amygdala, hippocampus and prefrontal cortex are implicated in the regulation of the hypothalamic pituitary adrenal (HPA) axis: the amygdala is implicated in activation of glucocorticoid secretion whereas the hippocampus and prefrontal cortex are largely inhibitory to glucocorticoid secretion (Herman et al., 2005). It has been documented that the amygdala, hippocampus and prefrontal cortex are concurrently vulnerable to stress exposure (McEwen et al., 2016; Henigsberg et al., 2019; Merino et al., 2021). Moreover, the three regions have been reported to work in concert to regulate the stress response: the amygdala detects potential danger; the hippocampus encodes environmental information associated with the stressor and the prefrontal cortex modulates associations between cues and stressor (Godoy et al., 2018).
The amygdala is a key structure involved in the regulation of fear (Ressler, 2010). Structural magnetic resonance imaging (MRI) research has shown that pre-natal maternal cortisol (Buss et al., 2012), depression (Wen et al., 2017), anxiety (Acosta et al., 2019), and disaster-related PNMS (Jones et al., 2019) are associated with increased whole amygdala volume in young children (Buss et al., 2012; Wen et al., 2017; Acosta et al., 2019) and in adolescents (Jones et al., 2019). The amygdala, however, is a heterogeneous structure composed of multiple nuclei that play distinct functional roles (Sah et al., 2003). For instance, the lateral and basal nuclei are involved in fear acquisition, and the central and medial nuclei are involved in fear expression and execution of fear responses (Ressler, 2010; Duvarci and Pare, 2014). To date, there is scant literature on associations between PNMS and offspring amygdala nuclei volumes.
The hippocampus plays a central role in learning and memory (Wingenfeld and Wolf, 2014). Three studies have shown that pre-natal maternal psychological distress (Lehtola et al., 2020; Moog et al., 2021) and anxiety (Wu et al., 2020) are associated with decreased whole hippocampal volume in fetuses (Wu et al., 2020) or in neonates (Lehtola et al., 2020; Moog et al., 2021); disaster-related PNMS has been associated with increased whole hippocampal volume in adolescents (Cao-Lei et al., 2021), while two other studies have reported no associations between pre-natal maternal cortisol (Buss et al., 2012) or stressful life events (Marečková et al., 2018) and whole hippocampal volume in children (Buss et al., 2012) or in young adults (Marečková et al., 2018). The mixed findings may be attributed to differences in types of stress and offspring age at scanning. As with the amygdala, the hippocampus has been divided into multiple subfields with distinct functional roles (Fogwe et al., 2018). For example, rodent evidence showed that PNMS was associated with a reduction of offspring dendritic arborization and synaptic density in the CA1 and CA3 subfields but not in others (Moog et al., 2021). There is little research on associations between PNMS and hippocampal subfield volumes in human offspring, however.
There is an emerging literature on PNMS and offspring brain functional connectivity. The prefrontal cortex modulates cognitive control functions and is mainly involved in working memory, self-regulatory, and goal-directed behaviors (McEwen and Morrison, 2013). It may act as a hub system that integrates sensory, affective, social, and memory-related information from the amygdala and the hippocampus to coordinate behavioral and peripheral physiological responses according to contextual demands that are appraised as stressful (Ginty et al., 2019). One study showed that pre-natal maternal anxiety was associated with reduced prefrontal gray matter density in child offspring (Buss et al., 2010).
Resting-state functional MRI studies on humans reported associations between PNMS and offspring whole amygdala-prefrontal resting-state functional connectivity (rs-FC). Specifically, pre-natal stressful life events were associated with decreased whole amygdala-medial prefrontal rs-FC in infants (Humphreys et al., 2020). Pre-natal maternal depression was associated with increased rs-FC from the whole amygdala to the anterior cingulate cortex and orbitofrontal cortex in infants (Qiu et al., 2015). Pre-natal maternal depression was associated with increased rs-FC from the whole amygdala to the ventromedial (Qiu et al., 2015) and the dorsal prefrontal cortex (Posner et al., 2016) in infants. The discrepancies noted here paint a potentially complex picture that depends on types of stress and on prefrontal subregions. Compared to the abundant research on the effects of PNMS on whole amygdala-prefrontal rs-FC, the influence of PNMS on whole hippocampus-prefrontal rs-FC has been much less explored. It remains unclear on associations between PNMS and functional connectivity of amygdala nuclei-prefrontal subregions and hippocampal subfields-prefrontal subregions.
In contrast to pre-natal maternal depression, anxiety, and stressful life events, studying PNMS deriving from sudden-onset natural disasters has several advantages, primarily because natural disasters are independent events that are not confounded by parental characteristics. In January 1998, five continuous days of freezing rain produced an ice storm in southern Quebec, Canada, that resulted in the failure of the regional power grid depriving millions of residents of electricity. In June 1998, we launched the world’s first prospective longitudinal disaster-related PNMS cohort: Project Ice Storm. This project recruited women who were pregnant during the crisis or became pregnant within 3 months following the ice storm, and comprehensively assessed three aspects of PNMS in each woman: objective hardship, subjective distress and cognitive appraisal of the crisis. The measurements of Project Ice Storm offspring began from age 6 months, and continued approximately every 2 years thereafter, till 19 years old.
Here, we studied 19-year-old young adult offspring from Project Ice Storm and typically developing controls from two datasets: the Autism Brain Imaging Data Exchange (ABIDE) and the Attention Deficit Hyperactivity Disorder-200 (ADHD-200). Our first goal was to examine group differences between young adult offspring exposed to disaster-related PNMS and controls in volumes of amygdala nuclei, hippocampal subfields and prefrontal subregions. Our secondary goal was to examine the group differences in rs-FC of amygdala nuclei-prefrontal subregions and hippocampal subfields-prefrontal subregions. Finally, we aimed to determine the extent to which the severity of disaster-related PNMS was associated with any volume and rs-FC showing between-group differences.
2. Materials and methods
2.1. Participants
2.1.1. Ice Storm
In June 1998, the initial Project Ice Storm cohort consisted of 176 women. At age 19, there were 39 young adult offspring (18M/21F) who underwent structural MRI and resting-state functional MRI scans. As shown in Supplementary Table 1, no significant differences were detected on the three aspects of PNMS between the current sample of 39 families and the 137 families who did not participate in MRI scanning at age 19; the difference was that the current sample of 39 families had significantly higher socioeconomic status. Among the 39 mothers, 30.8% (12/39) were middle class; 46.2% (18/39) were upper middle class; and 23.1% (9/39) were upper class. When the ice storm peaked on 9 January 1998, 28.2% (11/39) were within 3-month of conception; 30.8% (12/39) were in the 1st trimester of pregnancy; 23.1% (9/39) were in the 2nd trimester; and 17.9% (7/39) were in the 3rd trimester.
All phases of Project Ice Storm were approved by the Research Ethics Board of Douglas Mental Health University Institute. All the participants provided written informed consent.
2.1.2. Controls
Control participants were obtained from the ABIDE1 and the ADHD-2002. All control participants were typically developing participants who were characterized by absence of Autism Spectrum Disorder and Attention-Deficit/Hyperactivity Disorder diagnoses, as well as by absence of major neurological or psychiatric disorders. Despite no measurements of PNMS in controls, no natural disasters were recorded between 1996 and 1999 in the locations where the participants were recruited. As such, we can assume that the control participants were not systematically exposed to a population-level natural disaster. Control participants were selected to match Ice Storm participants on age, sex, handedness, intelligence quotient and in-scanner head motion (Supplementary Table 2).
2.2. Three aspects of prenatal maternal stress
2.2.1. Objective hardship
In June 1998, 5 months after the onset of the ice storm, the severity of maternal objective hardship experienced by pregnant women was assessed according to four dimensions of disaster exposure: Threat (e.g., injuries), Loss (e.g., loss of personal income), Scope (e.g., duration without electricity), and Change (e.g., temporary shelter) (Bromet and Dew, 1995; Laplante et al., 2007). Each dimension was scored on a scale of 0 (no exposure) to 8 (high exposure). A total score, referred to as Storm32, was calculated by summing scores across all four dimensions (Laplante et al., 2007). The test-retest reliability of Storm32 (assessed in the same women 6 years later) was satisfactory (r = 0.79) (St-Hilaire et al., 2015).
2.2.2. Subjective distress
In June 1998, maternal subjective distress was assessed using a validated 22-item French version (Brunet et al., 2003) of the Impact of Event Scale–Revised (IES-R) (Weiss et al., 1997), the gold-standard screening for post-traumatic stress disorder (PTSD). A total score of 33 is a cut-off for probable PTSD (Creamer et al., 2003). The scale rates the severity of symptoms in the preceding 7 days in three dimensions relevant to PTSD: intrusive thoughts, hyperarousal, and avoidance. Each dimension was scored of 0 (not at all) to 4 (extremely). The scale has satisfactory test-retest reliability for the total score (r = 0.76) (Brunet et al., 2003). Log-transformed values of the total score were used in the analyses due to skewed distribution.
2.2.3. Cognitive appraisal
In June 1998, maternal cognitive appraisal of the crisis was assessed using the following question: “Overall, what were the consequences of the ice storm on you and your family?.” Response options were rated as three options: “negative” (“–1”), “neutral” (“0”), and “positive” (“1”). We have shown that this measure has predictive validity by significantly correlating with child outcomes [e.g., BMI and central adiposity (Cao-Lei et al., 2016), C-peptide (Cao-Lei et al., 2018), and DNA methylation (Cao-Lei et al., 2015)].
2.3. MRI data acquisition
2.3.1. Ice Storm
MR images were acquired using a 3T Siemens MAGNETOM Trio TIM Syngo MRI scanner, with a 12-channel head coil. Anatomical images were obtained using a T1-weighted (T1w) Magnetization Prepared Rapid Gradient Echo sequence: 192 slices; Repetition Time (TR) = 2,400 s; Echo Time (TE) = 2.43 ms; slice thickness = 1 mm; flip angle = 8°; matrix = 256 × 256. Resting-sate functional images were acquired using a T2*-weighted echo-planar imaging sequence: 42 slices; TR = 2,600 ms; TE = 30 ms; flip angle = 90°; slice thickness = 3.4 mm; FoV = 218 mm, matrix = 64 × 64. Throughout the 5:01 min resting-state functional MRI scan, participants were instructed to lie still with their eyes open.
2.3.2. Controls
The scanning parameters of controls are described in Supplementary Table 3.
2.4. MRI data pre-processing
One Ice Storm participant was excluded from resting-state functional MRI pre-processing due to scan artifacts (Supplementary Figure 1). fMRIPrep 1.5.7 (Esteban et al., 2019) was used for pre-processing. The T1w image was corrected for intensity non-uniformity with N4BiasFieldCorrection (Tustison et al., 2010), distributed with ANTs 2.2.0 (Avants et al., 2008), and used as T1w-reference throughout the workflow. The T1w-reference was then skull-stripped with a Nipype implementation of the antsBrainExtraction.sh workflow (from ANTs), using OASIS30ANTs as target template. Brain tissue segmentation of gray matter, white matter and cerebrospinal fluid was performed on the brain-extracted T1w using fast FSL 5.0.9 (Zhang et al., 2001). Volume-based spatial normalization to the Montreal Neurological Institute (MNI) space was performed through non-linear registration with ants. Registration (ANTs 2.2.0), using brain-extracted versions of both T1w reference and the T1w template. For each of the blood-oxygen-level-dependent (BOLD) runs found per subject, the following pre-processing was performed. First, a reference volume and its skull-stripped version were generated using a custom methodology of fMRIPrep. A deformation field to correct for susceptibility distortions was estimated based on fMRIPrep’s fieldmap-less approach. The deformation field is that resulting from co-registering the BOLD reference to the same-subject T1w-reference with its intensity inverted (Huntenburg, 2014; Wang et al., 2017). Registration is performed with ants. Registration (ANTs 2.2.0), and the process regularized by constraining deformation to be non-zero only along the phase-encoding direction, and modulated with an average fieldmap template (Treiber et al., 2016). Based on the estimated susceptibility distortion, a corrected echo-planar imaging reference was calculated for a more accurate co-registration with the anatomical reference. The BOLD reference was then co-registered to the T1w reference using flirt FSL 5.0.9 (Jenkinson and Smith, 2001) with the boundary-based registration (Greve and Fischl, 2009) cost-function. Co-registration was configured with nine degrees of freedom to account for distortions remaining in the BOLD reference. Head-motion parameters with respect to the BOLD reference (transformation matrices, and six corresponding rotation and translation parameters) are estimated before any spatiotemporal filtering using mcflirt FSL 5.0.9 (Jenkinson et al., 2002). The BOLD time-series were resampled onto their original, native space by applying a single, composite transform to correct for head-motion and susceptibility distortions. The BOLD time-series were then resampled into the MNI space. Automatic removal of motion artifacts using independent component analysis (ICA-AROMA) (Pruim et al., 2015), was performed on the spatially normalized, pre-processed BOLD on MNI space time-series after removal of non-steady state volumes and spatial smoothing with an isotropic, Gaussian kernel of 6 mm FWHM. In addition, we conducted white matter and cerebrospinal fluid signal removal from the BOLD time series and temporally bandpass filtering (>0.01 Hz).
2.5. FreeSurfer segmentation
fMRIPrep pre-processed T1w brain regions in MNI space were segmented using FreeSurfer 7.1.1 (Fischl, 2012) and its library tool recon-all. The segmentation of the amygdala and the hippocampus was performed using segmentHA_T1.sh (Iglesias et al., 2015; Saygin et al., 2017). The amygdala nuclei and hippocampal subfields are shown in Figure 1: there were 9 amygdala nuclei, and 19 hippocampal subfields. AFNI 3dcalc was used to combine the hippocampal head and body: the 19 subfields were regrouped into 12 subfields. Eight prefrontal subregions were obtained based on the Desikan-Killiany Atlas (Figure 2). The quality of the segmentation was visually inspected with FreeView by X.L. The anterior amygdaloid area was absent for one Ice Storm participant and three controls (Supplementary Figure 2).
Figure 1. Segmentation of the hippocampus and the amygdala of one Ice Storm participant. The columns from left to right represent axial, coronal, and sagittal views, respectively. The hippocampal subfields and the amygdala nuclei (left and right) were labeled with different colors. The second-row zooms in the right hippocampus and the right amygdala without any labels as the reference. The third-row zooms in the right hippocampus and the right amygdala with labels as the reference.
Figure 2. Segmentation of the prefrontal cortex of one Ice Storm participant. The columns from left to right represent axial, coronal, and sagittal views, respectively. The prefrontal subregions (left and right) were labeled with different colors.
2.6. Between-group volumetric differences
Using IBM SPSS Statistics 22 (IBM Corp, 2019), an analysis of covariance (ANCOVA) model, controlling for TR, TE, and sex, was conducted to compare Ice Storm participants with controls on volumetric differences of: (1) whole amygdala; (2) whole hippocampus; (3) 9 amygdala nuclei; (4) 12 hippocampal subfields, and (5) 16 prefrontal subregions. The ANCOVA results were false discovery rate (FDR)-corrected for 18 comparisons based on 9 amygdala nuclei (left and right), 24 comparisons based on 12 hippocampal subfields (left and right), and 16 comparisons based on the eight prefrontal subregions (left and right). We used partial eta squared (η2p) as effect sizes (0.01 = small; 0.06 = medium; 0.14 = large) (Cohen, 1988; Lakens, 2013). We conducted sensitivity analyses by omitting the five participants with the 1.5T scanner. The sensitivity analyses indicated that there were no changes to the significance levels in the ANCOVA results.
2.7. Seed-to-seed rs-FC and between-group rs-FC differences
At the individual level, using the whole amygdala or whole hippocampus segmented in MNI space (left and right, separately) as the seed, Pearson correlation coefficients were calculated between the average BOLD time courses extracted from the whole amygdala or whole hippocampus and the average BOLD time courses of the prefrontal subregions. Using the 9 amygdala nuclei or 12 hippocampal subfields (left and right, separately) segmented in MNI space as seeds, Pearson correlation coefficients were calculated between the average BOLD time courses extracted from the 9 amygdala nuclei or 12 hippocampal subfields and the average BOLD time courses of the prefrontal subregions. Resultant seed-to-seed Pearson correlation coefficients were converted to normally distributed z-values.
Using the CONN functional connectivity toolbox 19b (Whitfield-Gabrieli and Nieto-Castanon, 2012), ANCOVA controlling for TR, TE, in-scanner eye status and sex, was applied to test differences between Ice Storm participants and controls in: (1) whole amygdala-prefrontal rs-FC; (2) whole hippocampus-prefrontal rs-FC; (3) amygdala nuclei-prefrontal rs-FC, and (4) hippocampal subfield-prefrontal rs-FC. Significance thresholds were set to p < 0.001 for uncorrected, and p < 0.05 for FDR-correction. Due to the segmentation failure of the anterior amygdaloid area, one Ice Storm participant and three controls were excluded from the rs-FC analyses of the whole amygdala and the anterior amygdaloid area. Finally, the ANCOVA results were FDR-corrected for 18 comparisons based on 9 amygdala nuclei (left and right), and 24 comparisons based on 12 hippocampal subfields (left and right). We conducted sensitivity analyses by omitting the five participants with the 1.5T scanner. The sensitivity analyses indicated that there were no changes to the significance levels in the ANCOVA results.
2.8. Association between PNMS and volume and rs-FC within the Ice Storm group
For volume and rs-FC showing significant FDR-corrected between-group differences, the following linear regressions within the Ice Storm group were conducted with IBM SPSS Statistics 22 (IBM Corp, 2019), controlling for sex, to determine: (1) associations between the three aspects of PNMS and volume and rs-FC and; (2) associations between volume and rs-FC and interactions of the three aspects of PNMS (objective hardship × subjective distress, objective hardship × cognitive appraisal, subjective distress × cognitive appraisal). Significant interactions at the uncorrected level (p < 0.05) were probed with PROCESS macro (Hayes and Preacher, 2013) to identify regions of significance. The above regression analyses were FDR-corrected for multiple comparisons based on three aspects of PNMS and the number of brain regions in which volume and rs-FC showing between-group differences.
3. Results
3.1. Larger volumes of the amygdala, hippocampus, and prefrontal cortex
3.1.1. Amygdala volume
Ice Storm participants had significantly larger bilateral whole amygdala volumes compared to controls (left: p = 0.001; right: p = 0.001, Table 1). Ice Storm participants had larger volumes, that survived FDR correction, of bilateral lateral, basal, central, medial, cortical, accessory basal nuclei, and corticoamygdaloid transition compared to controls; left central and right accessory basal nuclei showed large effect sizes (Table 1 and Figure 3). The volumes of anterior amygdaloid area and paralaminar nucleus did not differ between the two groups (Table 1).
Table 1. Ice Storm participants had larger volumes of whole amygdala, amygdala nuclei, whole hippocampus, hippocampal subfields, and prefrontal subregions.
Figure 3. Effects size estimates measured by partial eta squared (η2p) for the left (blue) and right (orange) whole amygdala and amygdala nuclei. η2p = 0.01 indicates a small effect; η2p = 0.06 indicates a medium effect, and η2p = 0.14 indicates a large effect. ***q < 0.001, **q < 0.01, *q < 0.05.
3.1.2. Hippocampal volume
Ice Storm participants had significantly larger bilateral whole hippocampal volumes compared to controls (left: p = 0.006; right: p = 0.004, Table 1). Ice Storm participants had larger volumes, that survived FDR correction, of the right CA1, the right HATA, the right molecular layer, the left fissure, bilateral tail, bilateral CA3, bilateral CA4, and bilateral DG, but smaller right fimbria volume compared to controls; right CA3 showed large effect size (Table 1 and Figure 4).
Figure 4. Effects size estimates measured by partial eta squared (η2p) for the left (blue) and right (orange) whole hippocampus and hippocampal subfields. η2p = 0.01 indicates a small effect; η2p = 0.06 indicates a medium effect, and η2p = 0.14 indicates a large effect. ***q < 0.001, **q < 0.01, *q < 0.05.
3.1.3. Prefrontal cortex volume
Ice Storm participants had larger left superior frontal volume, that survived FDR correction, compared to controls (p = 0.001, q = 0.016), while other remaining prefrontal subregions did not show significant between-group volumetric difference (Table 1).
3.2. Associations between PNMS and volumes of amygdala, hippocampus, and prefrontal cortex within the ice storm group (none survived FDR correction)
3.2.1. Amygdala volume
The linear regression results were FDR-corrected for 42 comparisons based on three aspects of PNMS and 14 amygdala nuclei showing between-group volumetric differences. No main effects of PNMS were observed for amygdala nuclei volumes. An interaction between objective hardship and subjective distress was observed for the right medial nucleus volume (p = 0.010): larger right medial nucleus volume was seen with a combination of higher objective hardship and higher subjective distress (Supplementary Figure 3). However, this interaction did not survive FDR correction (q = 0.420). An interaction of subjective distress and cognitive appraisal was observed for the right medial nucleus volume (p = 0.028): larger right medial nucleus volume was seen with a combination of higher subjective distress and more negative cognitive appraisal (Supplementary Figure 4). However, this interaction did not survive FDR correction (q = 0.588). An interaction of objective hardship and subjective distress was observed for the right cortical nucleus volume (p = 0.049): larger right cortical nucleus volume was seen with a combination of higher objective hardship and higher subjective distress (Supplementary Figure 5). However, this interaction did not survive FDR correction (q = 0.686).
3.2.2. Hippocampal volume
The linear regression results were FDR-corrected for 39 comparisons based on three aspects of PNMS and 13 hippocampal subfields showing between-group volumetric differences. The lower the maternal subjective distress, the larger the young adult offspring’s right hippocampal tail volume (beta = −0.437, p = 0.006, Supplementary Figure 6). However, this main effect did not survive FDR correction (q = 0.234). Similarly, we observed an interaction between subjective distress and cognitive appraisal for the left hippocampal tail volume (p = 0.046): larger left hippocampal tail volume was seen with a combination of lower subjective distress and more positive cognitive appraisal (Supplementary Figure 7). However, this interaction did not survive FDR correction (q = 0.322). In addition, we observed an interaction between objective hardship and subjective distress on the right CA4 volume (p = 0.049): larger right CA4 volume was seen with a combination of higher objective hardship and higher subjective distress (Supplementary Figure 8). However, this interaction did not survive FDR correction (q = 0.322). We also observed an interaction of subjective distress and cognitive appraisal on the right CA4 volume (p = 0.034). When we probed this interaction, there were no regions of significance (Supplementary Figure 9), and this interaction did not survive FDR correction (q = 0.322). Likewise, we found an interaction between objective hardship and subjective distress on the right DG volume (p = 0.041): larger right DG volume was seen with a combination of higher objective hardship and higher subjective distress (Supplementary Figure 10). However, this interaction did not survive FDR correction (q = 0.322). We also observed an interaction between subjective distress and cognitive appraisal on the right DG volume (p = 0.049), but there were no regions of significance (Supplementary Figure 11), and this interaction did not survive FDR correction (q = 0.322).
3.2.3. Prefrontal cortex volume
No main effects nor interactions were observed between PNMS and the left superior frontal volume.
3.3. Lower amygdala-prefrontal rs-FC
Compared to controls, Ice Storm participants exhibited lower rs-FC between the right whole amygdala and (1) the left caudal middle frontal (p = 0.003, q = 0.015); (2) the left superior frontal (p < 0.001, q = 0.003); (3) the left frontal pole (p < 0.001, q < 0.001), and (4) the right frontal pole (p = 0.005, q = 0.021) (Table 2). For the right amygdala nuclei, Ice Storm participants had lower rs-FC (1) from the right lateral nucleus to the left frontal pole; (2) from the right basal nucleus to the left superior frontal and the left caudal middle frontal; (3) from the right accessory basal nucleus to the right frontal pole and the left superior frontal; and (4) from the right paralaminar nucleus to bilateral superior frontal and bilateral caudal middle frontal (Supplementary Table 4 and Supplementary Figure 12).
Table 2. Ice Storm participants had lower amygdala-prefrontal rs-FC and lower hippocampus-prefrontal rs-FC compared to controls.
3.4. Lower hippocampus-prefrontal rs-FC
Ice Storm participants exhibited significantly lower FDR-corrected rs-FC between the right hippocampal tail and the left lateral orbitofrontal than controls (p < 0.001, q < 0.001) (Table 2). In addition, Ice Storm participants had lower rs-FC from the left parasubiculum to bilateral superior frontal and the right frontal pole compared to controls (shown in Supplementary Table 4 and Supplementary Figure 12).
3.5. Associations between PNMS and whole amygdala-prefrontal rs-FC within the Ice Storm group (none survived FDR correction)
No main effects of PNMS were observed for whole amygdala-prefrontal rs-FC. An interaction between objective hardship and cognitive appraisal was observed for rs-FC between the right whole amygdala and the right frontal pole (p = 0.022): lower rs-FC between the right whole amygdala and the right frontal pole was seen with a combination of higher objective hardship and more negative cognitive appraisal (Supplementary Figure 13). However, this interaction did not survive FDR correction (q = 0.264) for 12 comparisons based on three aspects of PNMS and four whole amygdala-prefrontal rs-FC showing between-group differences.
3.6. Associations between PNMS and hippocampal tail-prefrontal rs-FC within the ice storm group (none survived FDR correction)
The higher the maternal objective hardship, the lower rs-FC between the right hippocampal tail and the left lateral orbitofrontal (beta = −0.361, p = 0.028; Supplementary Figure 14). This main effect did not survive FDR correction (q = 0.084) for three comparisons based on three aspects of PNMS.
4. Discussion
To our knowledge, this is the first study examining volumes and rs-FC of the amygdala, hippocampus and prefrontal cortex in young adult offspring exposed in utero to varying levels of disaster-related PNMS. Primarily, we found that, compared to controls, young adult offspring exposed to disaster-related PNMS had larger volumes of the amygdala, hippocampus, and prefrontal cortex but lower amygdala-prefrontal connectivity and lower hippocampus-prefrontal connectivity. In addition, within the Ice Storm group, we found several associations between the severity of objective hardship or subjective distress, or interactions between PNMS variables, that explained variance in volume or functional connectivity, although none survived FDR corrections.
We observed that Ice Storm participants had larger bilateral whole amygdala volumes than controls. Larger amygdala may not be adaptive given that our previous Project Ice Storm findings at age 11½ indicated that larger right whole amygdala volume was associated with more severe externalizing behaviors (Jones et al., 2019). Our current results are inconsistent with a previous study finding that pre-natal stressful life events (e.g., break-up or divorce from partner, consideration of abortion, violence, serious illness or death in the family, financial difficulties) were associated with decreased whole amygdala volume in young adulthood (Mareckova et al., 2021). This discrepancy from our finding might be mainly attributed to different operationalizations of stress exposure; most of the stressors in the previous study (Mareckova et al., 2021) could have been brought on by parental characteristics that increase propensity to create stressful life events and that also confer genetic risk to the offspring that might be seen in amygdala development, whereas exposure to a natural disaster is an “independent” stressor that is outside of the control of the individual. In that sense, Project Ice Storm could be considered a natural experiment, and the current comparison between Ice Storm-exposed participants and controls involves random assignment to groups. Another possible explanation is that the stressors captured in the previous study (Mareckova et al., 2021) were all from the first half of pregnancy while Project Ice Storm included exposures from 3-months preconception to the very end of pregnancy. Although our previous Ice Storm results at age 11½ found that larger right whole amygdala volume was associated higher maternal subjective distress during 2nd and 3rd trimesters of pregnancy in boys, and higher maternal objective hardship in girls (Jones et al., 2019), at age 19 we found no associations between the severity of PNMS and bilateral whole amygdala volumes. It is possible, then, that the mere exposure to the ice storm in utero is the active ingredient in the effect on whole amygdala volume that persists into young adulthood.
Regarding amygdala nuclei, we found that, compared to controls, Ice Storm participants exhibited larger volumes of bilateral lateral, basal, central, medial, cortical, accessory basal nuclei, and corticoamygdaloid transition, but not of the anterior amygdaloid area or paralaminar nucleus. The absence of volumetric increases in the latter nuclei might be, in part, due to their low density of binding receptors for stress-inducing hormones and neurotransmitters (LeDoux, 2007). Evidence from rodents showed that glucocorticoids can directly bind to glucocorticoid receptors in the basolateral complex, consisting of lateral, basal, accessory basal and paralaminar nuclei, in which brief stress exposure triggers an increase in the spine density (Roozendaal et al., 2009). The basolateral complex is the largest and is the main input site of the amygdala (deCampo and Fudge, 2012). Central and medial nuclei, the main output area, are involved in processing glucocorticoid signaling and regulating autonomic, behavioral and hormonal response to stress via efferent projections to hypothalamus and the bed nucleus of the stria terminals. The superficial complex, including cortical, corticoamygdaloid transition and anterior amygdaloid area, are the major targets of olfactory projections (Sah et al., 2003), and are involved in selective social processing of the sensory inputs (Goossens et al., 2009). The volumetric increases in these nuclei may reflect stress-sensitive neurodevelopment in utero.
Unlike larger amygdala volume, Ice Storm participants exhibited lower rs-FC between the prefrontal cortex and the right whole amygdala, driven by the right basolateral complex. It is evident that basal and lateral nuclei receive fearful inputs from the visual cortex and project to prefrontal cortex (Ressler, 2010). Upon acute stress, the amygdala overactivates while the medial prefrontal cortex deactivates (Herman et al., 2005). The medial prefrontal cortex is largely inhibitory to HPA axis secretion, which in turn deactivates the amygdala (Bremner, 2006; Dedovic et al., 2009). Our finding suggests inadequate integration of amygdala fear acquisition into the prefrontal processing circuit in response to disaster-related PNMS.
As for the hippocampus, we found that, compared to controls, Ice Storm participants had larger bilateral whole hippocampal volumes. In typically developing general populations, hippocampal volume has been reported to increase with age from childhood to young adulthood (Xu et al., 2020). Our results of larger bilateral whole hippocampal volumes may be explained, in part, by the “predictive adaptive response hypothesis” (Gluckman et al., 2005), which indicates that exposure to stress in the intra-utero environment biases protective stress responses to better adapt the development of an organism to the ex-utero environment (Sandman et al., 2012; Glover et al., 2017) by overdevelopment of the hippocampus. Further, we identified hippocampal subfield-specific volumetric alterations: larger volumes of the HATA, CA1, molecular layer, fissure, tail, CA3, CA4, and DG, but smaller fimbria volume. The HATA is tightly co-located and interconnected with the amygdala at the cellular level, and is involved in memory processing (Li et al., 2021); the CA1 is active in pattern completion while the CA3 is active in pattern separation, and the DG is involved in both processes (Postel et al., 2019). Among the hippocampal subfields, the vulnerability of the CA3 has been most-frequently reported possibly due to its role as the main target for glucocorticoids (Cai et al., 2007; Jia et al., 2010; Marečková et al., 2018; Zhang et al., 2019). In line with this, the CA3 in our study was the most vulnerable with the largest between-group effect size (η2p = 0.168). In contrast, the smaller right fimbria volume possibly indicates inconsistent regulatory roles of the hippocampal subfields in response to PNMS. The hippocampal tail is relevant for spatial information and negative emotion (Li et al., 2021). Our study extends prior research showing larger right hippocampal tail volume, but lower rs-FC between the right hippocampal tail and the left lateral orbitofrontal in young adult offspring of mothers exposed to natural disaster-related pre-natal stress. Interestingly, the whole hippocampus rs-FC did not differ between Ice Storm participants and controls, suggesting that subfield rs-FC, instead of whole hippocampus rs-FC, better underlies the neuropathology of PNMS.
In our study, Ice Storm participants had larger left superior frontal volume than controls. The superior frontal, located at the superior part of the prefrontal cortex, has been reported to be involved in motor control tasks, working memory, and higher cognitive processing (Li et al., 2013). We found that, compared to controls, Ice Storm participants had lower rs-FC between the right whole amygdala and the left superior frontal extending to the left caudal middle frontal and bilateral frontal poles, and lower rs-FC between the right hippocampal tail and the left lateral orbitofrontal cortex, which consistently point to decreased prefrontal integration in amygdala and hippocampal circuitry of pre-natal stress. Further, these results suggest that the prefrontal integration of the amygdala and hippocampal information is subregion-specific. It has been reported that the left superior frontal gyrus mainly contributes to working memory (du Boisgueheneuc et al., 2006) and the left caudal middle frontal gyrus mainly contributes to self-initiated elaborative strategies (Husa et al., 2017); frontal poles are mainly involved in action selection (Kovach et al., 2012). We propose that these subregions contribute differently to the processes of high-level cognitive functions via integrating information from the amygdala. Lateral orbitofrontal cortex mainly contributes to decision-making by combining prior with current information (Nogueira et al., 2017). In addition, the connectivity between lateral orbitofrontal cortex and hippocampal tail may process memory consolidation of information.
Within the Ice Storm group, we tested associations between the severity of PNMS and amygdala nuclei volumes at the uncorrected level. These overall trends extend our previous findings of PNMS influences on whole amygdala volume at age 11½ (Jones et al., 2019) to specific amygdala nuclei at age 19. We found that higher objective hardship, higher subjective distress and more negative cognitive appraisal interacted to predict larger medial nucleus volume, and that higher objective hardship and higher subjective distress interacted to predict larger volumes of medial and cortical nuclei. However, these findings should be interpreted with caution given that none of these interactions survived FDR correction.
Similar to the whole amygdala, no robustly significant associations were observed between the severity of PNMS and the whole hippocampal volume, in line with previous reports of no associations between pre-natal stressful life events and the whole hippocampal volume in young adult offspring (Favaro et al., 2015; Marečková et al., 2018). We also tested the associations between the severity of PNMS and hippocampal subfield volumes; none of the results survived FDR correction. The results from Project Ice Storm at age 11½ indicated that (1) greater maternal objective hardship was associated with increased right whole hippocampal volume in girls; (2) COMT genotype moderated the association between maternal objective hardship and right whole hippocampal volume in boys; and (3) COMT genotype moderated the association between maternal subjective distress and right whole hippocampal volume in girls (Cao-Lei et al., 2021). Some differences between age 11½ and the current study are that the age 11½ analyses were conducted with whole hippocampal volume only, and that the age 11½ analyses were conducted with boys and girls separately. In comparison, the current study further included analyses at the level of hippocampal subfields. In this sense, future studies may be needed to identify whether genetics and sex moderate associations between PNMS and hippocampal subfield volumes.
Within the Ice Storm group, we tested associations between the severity of PNMS and hippocampal subfield volumes at the uncorrected level. These results are worth mentioning, despite their failure to survive FDR corrections, in order to inform future hypotheses. Unlike our findings suggesting that higher PNMS is consistently associated with larger amygdala nuclei, these findings suggest that PNMS predicts larger hippocampal subfield volumes in a stress- and subfield- specific manner (Herman et al., 2005): higher subjective distress predicted smaller tail volume while higher objective hardship and higher subjective distress interacted to predict larger volumes of CA4 and DG. We speculate that the observed inconsistency might be attributed to much higher heterogeneity for hippocampal subfields than for amygdala nuclei. Within the Ice Storm group, at the uncorrected level, we observed that a combination of higher maternal objective hardship and more negative maternal cognitive appraisal jointly predicted lower offspring amygdala-frontal rs-FC whereas higher maternal objective hardship predicted lower offspring hippocampus-lateral orbitofrontal rs-FC, indicating that prefrontal regulatory circuits of amygdala and hippocampus might be stress- and subregion- specific. As noted, however, these findings should be interpreted with caution given that the aforementioned main effects and interactions failed to survive FDR correction.
Despite the observed group differences, and the ability of the PNMS to predict amygdala and hippocampal structure and function within the Ice Storm group, one must be cautious in making causal interpretations. It is possible that the results reflect the direct effects of in utero exposure to PNMS on the fetal brain that endures to age 19 years. It may also be the case, however, that the PNMS created a physiological vulnerability to post-natal environmental stressors that could be responsible for indirect effects on brain 19 years later.
There were some limitations. First, our sample was small which limited the ability to go deeper in our analyses to test moderations by sex and timing in utero. The sample size also limits our ability to determine the extent to which reactivity to any post-natal environmental stressors might have mediated associations between PNMS and offspring brain at age 19. Second, since MRI data for Ice Storm and control groups were collected from different scanning sites, we minimized this limitation by controlling for the key scanning parameters. Third, the socioeconomic status of the families in this study is higher than the median of the population from which they were recruited in 1998, such that the current results may not generalize to lower class populations.
In conclusion, our findings highlight long-lasting effects of maternal stress before birth on the volumes and rs-FC of the amygdala, the hippocampus and the prefrontal cortex 19 years after exposure, in which their vulnerability might be stress- and subdivision- specific.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
All phases of Project Ice Storm were approved by the Research Ethics Board of Douglas Mental Health University Institute. All the participants provided written informed consent.
Author contributions
XL contributed to the design, control MRI data acquisition, all analyses, and the draft of the manuscript. MQ contributed to the volume and rs-FC analysis methodology. DL contributed to the statistical analysis and MRI data acquisition. GE contributed to the statistical analyses. SJ contributed to the Ice Storm MRI data acquisition. SK and PR-N contributed to the supervision. All authors have approved the final manuscript.
Funding
Project Ice Storm was supported by grants to SK and colleagues by the Canadian Institutes of Health Research (MOP-57849, MOP-111177, and MOP-125892). XL was supported by the China Scholarship Council (201906170056).
Acknowledgments
We thank the families of Project Ice Storm participants. For controls recruited in our study, we thank the contributors at the nine centres from the ABIDE (first and second releases; https://fcon_1000.projects.nitrc.org/indi/abide/) and the ADHD-200 (https://fcon_1000.projects.nitrc.org/indi/adhd200/) for their efforts in the collection, organization and sharing of their datasets. We also thank Sulantha Mathotaarachchi and Vladimir S. Fonov from PR-N Lab for programming support.
Conflict of interest
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnhum.2023.1094039/full#supplementary-material
Footnotes
- ^ https://fcon_1000.projects.nitrc.org/indi/abide/
- ^ https://fcon_1000.projects.nitrc.org/indi/adhd200/
References
Acosta, H., Tuulari, J., Scheinin, N., Hashempour, N., Rajasilta, O., Lavonius, T., et al. (2019). Maternal pregnancy-related anxiety is associated with sexually dimorphic alterations in amygdala volume in 4-year-old children. Front. Behav. Neurosci. 13:175.
Avants, B., Epstein, C., Grossman, M., and Gee, J. (2008). Symmetric diffeomorphic image registration with cross-correlation: Evaluating automated labeling of elderly and neurodegenerative brain. Med. Image Anal. 12, 26–41. doi: 10.1016/j.media.2007.06.004
Barker, D. (2003). The developmental origins of adult disease. Eur. J. Epidemiol. 18, 733–736. doi: 10.1023/A:1025388901248
Bremner, J. (2006). Traumatic stress: Effects on the brain. Dialogues Clin. Neurosci. 8, 445–461. doi: 10.31887/DCNS.2006.8.4/jbremner
Bromet, E., and Dew, M. (1995). Review of psychiatric epidemiologic research on disasters. Epidemiol. Rev. 17, 113–119. doi: 10.1093/oxfordjournals.epirev.a036166
Brunet, A., St-Hilaire, A., Jehel, L., and King, S. (2003). Validation of a French version of the impact of event scale-revised. Can. J. Psychiatry 48, 56–61. doi: 10.1177/070674370304800111
Buss, C., Davis, E., Muftuler, L., Head, K., and Sandman, C. (2010). High pregnancy anxiety during mid-gestation is associated with decreased gray matter density in 6-9-year-old children. Psychoneuroendocrinology 35, 141–153. doi: 10.1016/j.psyneuen.2009.07.010
Buss, C., Davis, E., Shahbaba, B., Pruessner, J., Head, K., and Sandman, C. (2012). Maternal cortisol over the course of pregnancy and subsequent child amygdala and hippocampus volumes and affective problems. Proc. Natl. Acad. Sci. U.S.A. 109, E1312–E1319. doi: 10.1073/pnas.1201295109
Cai, Q., Zhu, Z., Li, H., Fan, X., Jia, N., Bai, Z., et al. (2007). Prenatal stress on the kinetic properties of Ca2+ and K+ channels in offspring hippocampal CA3 pyramidal neurons. Life Sci. 80, 681–689. doi: 10.1016/j.lfs.2006.10.016
Cao-Lei, L., Dancause, K., Elgbeili, G., Laplante, D., Szyf, M., and King, S. (2016). Pregnant women’s cognitive appraisal of a natural disaster affects their children’s BMI and central adiposity via DNA methylation: Project ice storm. Early Hum. Dev. 103, 189–192. doi: 10.1016/j.earlhumdev.2016.09.013
Cao-Lei, L., Dancause, K., Elgbeili, G., Laplante, D., Szyf, M., and King, S. (2018). DNA methylation mediates the effect of maternal cognitive appraisal of a disaster in pregnancy on the child’s C-peptide secretion in adolescence: Project ice storm. PLoS One 13:e0192199. doi: 10.1371/journal.pone.0192199
Cao-Lei, L., Elgbeili, G., Massart, R., Laplante, D., Szyf, M., and King, S. (2015). Pregnant women’s cognitive appraisal of a natural disaster affects DNA methylation in their children 13 years later: Project ice storm. Transl. Psychiatry 5:e515. doi: 10.1038/tp.2015.13
Cao-Lei, L., Yogendran, S., Dufoix, R., Elgbeili, G., Laplante, D., and King, S. (2021). Prenatal maternal stress from a natural disaster and hippocampal volumes: Gene-by-environment interactions in young adolescents from project ice storm. Front. Behav. Neurosci. 15:706660. doi: 10.3389/fnbeh.2021.706660
Cohen, J. (1988). Chapter 8. The analysis of variance and covariance. Statistical power analysis for the behavioral sciences. New York, NY: Routledge, 273–406. doi: 10.1016/B978-0-12-179060-8.50013-X
Creamer, M., Bell, R., and Failla, S. (2003). Psychometric properties of the impact of event scale - revised. Behav. Res. Ther. 41, 1489–1496. doi: 10.1016/j.brat.2003.07.010
deCampo, D., and Fudge, J. (2012). Where and what is the paralaminar nucleus? A review on a unique and frequently overlooked area of the primate amygdala. Neurosci. Biobehav. Rev. 36, 520–535. doi: 10.1016/j.neubiorev.2011.08.007
Dedovic, K., D’Aguiar, C., and Pruessner, J. (2009). What stress does to your brain: A review of neuroimaging studies. Can. J. Psychiatry 54, 6–15. doi: 10.1177/070674370905400104
du Boisgueheneuc, F., Levy, R., Volle, E., Seassau, M., Duffau, H., Kinkingnehun, S., et al. (2006). Functions of the left superior frontal gyrus in humans: A lesion study. Brain 129(Pt 12), 3315–3328. doi: 10.1093/brain/awl244
Duvarci, S., and Pare, D. (2014). Amygdala microcircuits controlling learned fear. Neuron 82, 966–980. doi: 10.1016/j.neuron.2014.04.042
Esteban, O., Markiewicz, C., Blair, R., Moodie, C., Isik, A., Erramuzpe, A., et al. (2019). fMRIPrep: A robust preprocessing pipeline for functional MRI. Nat. Methods 16, 111–116. doi: 10.1038/s41592-018-0235-4
Favaro, A., Tenconi, E., Degortes, D., Manara, R., and Santonastaso, P. (2015). Neural correlates of prenatal stress in young women. Psychol. Med. 45, 2533–2543. doi: 10.1017/S003329171500046X
Fogwe, L., Reddy, V., and Mesfin, F. (2018). Neuroanatomy, hippocampus. Tampa, FL: StatPearls Publishing.
Ginty, A., Kraynak, T., Kuan, D., and Gianaros, P. (2019). Ventromedial prefrontal cortex connectivity during and after psychological stress in women. Psychophysiology 56:e13445. doi: 10.1111/psyp.13445
Glover, L., Schoenfeld, T., Karlsson, R., Bannerman, D., and Cameron, H. (2017). Ongoing neurogenesis in the adult dentate gyrus mediates behavioral responses to ambiguous threat cues. PLoS Biol. 15:e2001154. doi: 10.1371/journal.pbio.2001154
Gluckman, P., Hanson, M., and Spencer, H. (2005). Predictive adaptive responses and human evolution. Trends Ecol. Evol. 20, 527–533. doi: 10.1016/j.tree.2005.08.001
Godoy, L., Rossignoli, M., Delfino-Pereira, P., Garcia-Cairasco, N., and de Lima Umeoka, E. (2018). A comprehensive overview on stress neurobiology: Basic concepts and clinical implications. Front. Behav. Neurosci. 12:127. doi: 10.3389/fnbeh.2018.00127
Goossens, L., Kukolja, J., Onur, O., Fink, G., Maier, W., Griez, E., et al. (2009). Selective processing of social stimuli in the superficial amygdala. Hum. Brain Mapp. 30, 3332–3338. doi: 10.1002/hbm.20755
Greve, D., and Fischl, B. (2009). Accurate and robust brain image alignment using boundary-based registration. Neuroimage 48, 63–72. doi: 10.1016/j.neuroimage.2009.06.060
Hayes, A., and Preacher, K. (2013). “Conditional process modeling: Using structural equation modeling to examine contingent causal processes,” in Structural equation modeling: A second course, eds G. R. Hancock and R. O. Mueller (Charlotte, NC: Information Age Publishing), 219–266.
Henigsberg, N., Kalember, P., Petrović, Z., and Šečić, A. (2019). Neuroimaging research in posttraumatic stress disorder - Focus on amygdala, hippocampus and prefrontal cortex. Prog. Neuropsychopharmacol. Biol. Psychiatry 90, 37–42. doi: 10.1016/j.pnpbp.2018.11.003
Herman, J., Ostrander, M., Mueller, N., and Figueiredo, H. (2005). Limbic system mechanisms of stress regulation: Hypothalamo-pituitary-adrenocortical axis. Prog. Neuropsychopharmacol. Biol. Psychiatry 29, 1201–1213. doi: 10.1016/j.pnpbp.2005.08.006
Humphreys, K., Camacho, M., Roth, M., and Estes, E. (2020). Prenatal stress exposure and multimodal assessment of amygdala-medial prefrontal cortex connectivity in infants. Dev. Cogn. Neurosci. 46:100877. doi: 10.1016/j.dcn.2020.100877
Huntenburg, J. (2014). Evaluating nonlinear coregistration of BOLD EPI and T1w images. Berlin: Freie Universität Berlin.
Husa, R., Gordon, B., Cochran, M., Bertolin, M., Bond, D., and Kirchhoff, B. (2017). Left caudal middle frontal gray matter volume mediates the effect of age on self-initiated elaborative encoding strategies. Neuropsychologia 106, 341–349. doi: 10.1016/j.neuropsychologia.2017.10.004
Iglesias, J., Augustinack, J., Nguyen, K., Player, C., Player, A., Wright, M., et al. (2015). A computational atlas of the hippocampal formation using ex vivo, ultra-high resolution MRI: Application to adaptive segmentation of in vivo MRI. Neuroimage 115, 117–137. doi: 10.1016/j.neuroimage.2015.04.042
Jenkinson, M., Bannister, P., Brady, M., and Smith, S. (2002). Improved optimization for the robust and accurate linear registration and motion correction of brain images. Neuroimage 17, 825–841. doi: 10.1006/nimg.2002.1132
Jenkinson, M., and Smith, S. (2001). A global optimisation method for robust affine registration of brain images. Med. Image Anal. 5, 143–156. doi: 10.1016/S1361-8415(01)00036-6
Jia, N., Yang, K., Sun, Q., Cai, Q., Li, H., Cheng, D., et al. (2010). Prenatal stress causes dendritic atrophy of pyramidal neurons in hippocampal CA3 region by glutamate in offspring rats. Dev. Neurobiol. 70, 114–125. doi: 10.1002/dneu.20766
Jones, S., Dufoix, R., Laplante, D., Elgbeili, G., Patel, R., Chakravarty, M., et al. (2019). Larger amygdala volume mediates the association between prenatal maternal stress and higher levels of externalizing behaviors: Sex specific effects in project ice storm. Front. Hum. Neurosci. 13:144. doi: 10.3389/fnhum.2019.00144
Kovach, C., Daw, N., Rudrauf, D., Tranel, D., O’Doherty, J., and Adolphs, R. (2012). Anterior prefrontal cortex contributes to action selection through tracking of recent reward trends. J. Neurosci. 32, 8434–8442. doi: 10.1523/JNEUROSCI.5468-11.2012
Lafortune, S., Laplante, D., Elgbeili, G., Li, X., Lebel, S., Dagenais, C., et al. (2021). Effect of natural disaster-related prenatal maternal stress on child development and health: A meta-analytic review. Int. J. Environ. Res. Public Health 18:8332. doi: 10.3390/ijerph18168332
Lakens, D. (2013). Calculating and reporting effect sizes to facilitate cumulative science: A practical primer for t-tests and ANOVAs. Front. Psychol. 4:863. doi: 10.3389/fpsyg.2013.00863
Laplante, D., Zelazo, P., Brunei, A., and King, S. (2007). Functional play at 2 years of age: Effects of prenatal maternal stress. Infancy 12, 69–93. doi: 10.1111/j.1532-7078.2007.tb00234.x
Lautarescu, A., Craig, M., and Glover, V. (2020). Prenatal stress: Effects on fetal and child brain development. Int. Rev. Neurobiol. 150, 17–40. doi: 10.1016/bs.irn.2019.11.002
Lehtola, S., Tuulari, J., Scheinin, N., Karlsson, L., Parkkola, R., Merisaari, H., et al. (2020). Newborn amygdalar volumes are associated with maternal prenatal psychological distress in a sex-dependent way. Neuroimage Clin. 28:102380. doi: 10.1016/j.nicl.2020.102380
Li, L., Pan, N., Zhang, L., Lui, S., Huang, X., Xu, X., et al. (2021). Hippocampal subfield alterations in pediatric patients with post-traumatic stress disorder. Soc. Cogn. Affect. Neurosci. 16, 334–344. doi: 10.1093/scan/nsaa162
Li, W., Qin, W., Liu, H., Fan, L., Wang, J., Jiang, T., et al. (2013). Subregions of the human superior frontal gyrus and their connections. Neuroimage 78, 46–58. doi: 10.1016/j.neuroimage.2013.04.011
Li, X., Laplante, D., Paquin, V., Lafortune, S., Elgbeili, G., and King, S. (2022). Effectiveness of cognitive behavioral therapy for perinatal maternal depression, anxiety and stress: A systematic review and meta-analysis of randomized controlled trials. Clin. Psychol. Rev. 92:102129. doi: 10.1016/j.cpr.2022.102129
Mareckova, K., Marecek, R., Andryskova, L., Brazdil, M., and Nikolova, Y. (2021). Impact of prenatal stress on amygdala anatomy in young adulthood: Timing and location matter. Biol. Psychiatry Cogn. Neurosci. Neuroimaging 7, 231–238. doi: 10.1016/j.bpsc.2021.07.009
Marečková, K., Mareček, R., Bencurova, P., Klánová, J., Dušek, L., and Brázdil, M. (2018). Perinatal stress and human hippocampal volume: Findings from typically developing young adults. Sci. Rep. 8:4696. doi: 10.1038/s41598-018-23046-6
McEwen, B., and Morrison, J. (2013). The brain on stress: Vulnerability and plasticity of the prefrontal cortex over the life course. Neuron 79, 16–29. doi: 10.1016/j.neuron.2013.06.028
McEwen, B., Nasca, C., and Gray, J. (2016). Stress effects on neuronal structure: Hippocampus, amygdala, and prefrontal cortex. Neuropsychopharmacology 41, 3–23. doi: 10.1038/npp.2015.171
Merino, E., Raya-Salom, D., Teruel-Martí, V., Adell, A., Cervera-Ferri, A., and Martínez-Ricós, J. (2021). Effects of acute stress on the oscillatory activity of the hippocampus-amygdala-prefrontal cortex network. Neuroscience 476, 72–89. doi: 10.1016/j.neuroscience.2021.09.009
Moog, N., Nolvi, S., Kleih, T., Styner, M., Gilmore, J., Rasmussen, J., et al. (2021). Prospective association of maternal psychosocial stress in pregnancy with newborn hippocampal volume and implications for infant social-emotional development. Neurobiol. Stress 15:100368. doi: 10.1016/j.ynstr.2021.100368
Nogueira, R., Abolafia, J., Drugowitsch, J., Balaguer-Ballester, E., Sanchez-Vives, M., and Moreno-Bote, R. (2017). Lateral orbitofrontal cortex anticipates choices and integrates prior with current information. Nat. Commun. 8:14823. doi: 10.1038/ncomms14823
Posner, J., Cha, J., Roy, A., Peterson, B., Bansal, R., Gustafsson, H., et al. (2016). Alterations in amygdala-prefrontal circuits in infants exposed to prenatal maternal depression. Transl. Psychiatry 6:e935. doi: 10.1038/tp.2016.146
Postel, C., Viard, A., André, C., Guénolé, F., de Flores, R., Baleyte, J., et al. (2019). Hippocampal subfields alterations in adolescents with post-traumatic stress disorder. Hum. Brain Mapp. 40, 1244–1252. doi: 10.1002/hbm.24443
Pruim, R., Mennes, M., van Rooij, D., Llera, A., Buitelaar, J., and Beckmann, C. F. (2015). ICA-AROMA: A robust ICA-based strategy for removing motion artifacts from fMRI data. Neuroimage 112, 267–277. doi: 10.1016/j.neuroimage.2015.02.064
Qiu, A., Anh, T., Li, Y., Chen, H., Rifkin-Graboi, A., Broekman, B., et al. (2015). Prenatal maternal depression alters amygdala functional connectivity in 6-month-old infants. Transl. Psychiatry 5:e508. doi: 10.1038/tp.2015.3
Ressler, K. (2010). Amygdala activity, fear, and anxiety: Modulation by stress. Biol. Psychiatry 67, 1117–1119. doi: 10.1016/j.biopsych.2010.04.027
Roozendaal, B., McEwen, B., and Chattarji, S. (2009). Stress, memory and the amygdala. Nat. Rev. Neurosci. 10, 423–433. doi: 10.1038/nrn2651
Sah, P., Faber, E., Lopez De Armentia, M., and Power, J. (2003). The amygdaloid complex: Anatomy and physiology. Physiol. Rev. 83, 803–834. doi: 10.1152/physrev.00002.2003
Sandman, C., Davis, E., and Glynn, L. (2012). Prescient human fetuses thrive. Psychol. Sci. 23, 93–100. doi: 10.1177/0956797611422073
Saygin, Z., Kliemann, D., Iglesias, J., van der Kouwe, A., Boyd, E., Reuter, M., et al. (2017). High-resolution magnetic resonance imaging reveals nuclei of the human amygdala: Manual segmentation to automatic atlas. Neuroimage 155, 370–382. doi: 10.1016/j.neuroimage.2017.04.046
St-Hilaire, A., Steiger, H., Liu, A., Laplante, D., Thaler, L., Magill, T., et al. (2015). A prospective study of effects of prenatal maternal stress on later eating-disorder manifestations in affected offspring: Preliminary indications based on the Project Ice Storm cohort. Int. J. Eat. Disord. 48, 512–516. doi: 10.1002/eat.22391
Treiber, J., White, N., Steed, T., Bartsch, H., Holland, D., Farid, N., et al. (2016). Characterization and correction of geometric distortions in 814 diffusion weighted images. PLoS One 11:e0152472. doi: 10.1371/journal.pone.0152472
Tustison, N., Avants, B., Cook, P., Zheng, Y., Egan, A., Yushkevich, P., et al. (2010). N4ITK: Improved N3 bias correction. IEEE Trans. Med. Imaging 29, 1310–1320. doi: 10.1109/TMI.2010.2046908
Wang, S., Peterson, D., Gatenby, J., Li, W., Grabowski, T., and Madhyastha, T. (2017). Evaluation of field map and nonlinear registration methods for correction of susceptibility artifacts in diffusion MRI. Front. Neuroinform. 11:17. doi: 10.3389/fninf.2017.00017
Weiss, D., Marmar, C., Wilson, J., and Keane, T. (1997). Assessing psychological trauma and PTSD. New York, NY: The Guilford Press.
Wen, D., Poh, J., Ni, S., Chong, Y., Chen, H., Kwek, K., et al. (2017). Influences of prenatal and postnatal maternal depression on amygdala volume and microstructure in young children. Transl. Psychiatry 7:e1103. doi: 10.1038/tp.2017.74
Whitfield-Gabrieli, S., and Nieto-Castanon, A. (2012). Conn: A functional connectivity toolbox for correlated and anticorrelated brain networks. Brain Connect. 2, 125–141. doi: 10.1089/brain.2012.0073
Wingenfeld, K., and Wolf, O. (2014). Stress, memory, and the hippocampus. Front. Neurol. Neurosci. 34, 109–120. doi: 10.1159/000356423
Wu, Y., Lu, Y., Jacobs, M., Pradhan, S., Kapse, K., Zhao, L., et al. (2020). Association of prenatal maternal psychological distress with fetal brain growth, metabolism, and cortical maturation. JAMA Netw. Open 3:e1919940. doi: 10.1001/jamanetworkopen.2019.19940
Xu, Q., Zuo, C., Liao, S., Long, Y., and Wang, Y. (2020). Abnormal development pattern of the amygdala and hippocampus from childhood to adulthood with autism. J. Clin. Neurosci. 78, 327–332. doi: 10.1016/j.jocn.2020.03.049
Yong Ping, E., Laplante, D., Elgbeili, G., Hillerer, K., Brunet, A., O’Hara, M., et al. (2015). Prenatal maternal stress predicts stress reactivity at 2½ years of age: The Iowa flood study. Psychoneuroendocrinology 56, 62–78. doi: 10.1016/j.psyneuen.2015.02.015
Zhang, H., Li, Q., Shang, Y., Xiao, X., Xu, X., Zhang, J., et al. (2019). Effect of prenatal stress on neural oscillations in developing hippocampal formation. Prog. Neuropsychopharmacol. Biol. Psychiatry 89, 456–464. doi: 10.1016/j.pnpbp.2018.10.019
Keywords: volume, resting-state functional connectivity, prenatal maternal stress, amygdala nuclei, hippocampal subfields, prefrontal cortex
Citation: Li X, Qureshi MNI, Laplante DP, Elgbeili G, Jones SL, King S and Rosa-Neto P (2023) Neural correlates of disaster-related prenatal maternal stress in young adults from Project Ice Storm: Focus on amygdala, hippocampus, and prefrontal cortex. Front. Hum. Neurosci. 17:1094039. doi: 10.3389/fnhum.2023.1094039
Received: 09 November 2022; Accepted: 11 January 2023;
Published: 01 February 2023.
Edited by:
Chuanliang Han, Chinese Academy of Sciences (CAS), ChinaReviewed by:
Meijia Li, Vrije Universiteit Brussel, BelgiumXianglian Jia, Stanford University, United States
Jiahua Xu, Chinese Institute for Brain Research, China
Copyright © 2023 Li, Qureshi, Laplante, Elgbeili, Jones, King and Rosa-Neto. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Suzanne King, 
suzanne.king@mcgill.ca