Participants
Healthy participants for this study were recruited from Growing Up in Singapore Towards Healthy Outcomes (GUSTO), a longitudinal, population-based birth cohort. Pregnant women aged 18 years or older were recruited at their first trimester of pregnancy from two main public hospitals in Singapore between June 2009 and December 201086. Mother-child dyads were followed throughout pregnancy and beyond. This study followed the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) reporting guidelines for cohort studies87. To minimize potential birth complication effects, we only included children with gestational age at birth ≥ 34 weeks, birth weight ≥ 2000 g, and a 5-minute Apgar score ≥ 8 in our analysis. Furthermore, children with a neuropsychiatric diagnosis or on psychoactive mediations were also excluded. All mother-child dyads utilized in this study were healthy with no pre-identified neurodevelopmental problems.
Maternal Mental Health
Prenatal maternal mental health measures were collected during week 26 of the pregnancy. Parents completed three mental health questionnaires: BDI, STAI and EPDS.
The Beck Depression Inventory, Second Edition (BDI)88 is a list of 21 items describing common depressive symptoms. Each item contains four or seven statements describing varying severity of a common depressive symptom. Participants selected the statement that best described how they felt for the past 2 weeks. The EPDS89 has 10 items of depressive symptoms, and participants indicated how much each item described how they were feeling for past 7 days on a 4-point Likert scale. The STAI19 consisted of 40 items that are associated with the presence or absence of anxiety. For the first 20 items, participants responded to how much each item described how they felt right now on a 4-point Likert scale; for the next 20 items, they responded to how much the item described how they generally felt.
Previous work by Desiree et al.20 identified a bi-factor model by combining all the items from these three questionnaires, with a specific positive mental health factor and a general negative mental health factor. We used this positive mental health factor as our measure of positive maternal mental health and the general negative factor as our measure of negative maternal mental health.
Language
At 24 months of age, the Bayley Scales of Infant and Toddler Development, Third Edition90, was used to assess the child’s development in the domains of cognition, language, motor skills, socioemotional behaviours, and adaptability. For this study, we combined the child’s receptive and expressive communication language scores into a composite language score for further analysis.
At 9 years of age, the Wechsler Individual Achievement Test, 3rd Edition (WIAT-3) was administered. WIAT-3 assesses the child’s development in terms of spelling, oral fluency, pseudoword decoding and numerical ability. We used the oral fluency score as a measure of language ability at later timepoints.
Executive Function
Behaviour Rating Inventory of Executive Function, 2nd Edition (BRIEF-2; parent-report form)91 was administered at age 7 years (n = 620). BRIEF-2 is a rating scale that assesses everyday behaviours reflecting executive functions across the school-age span (ages 5–18) and contains 63 items within nine clinical scales and three indexes. The Shift and Emotional Control scales comprise the Emotional Regulation Index, the Inhibit and Self-Monitor scales comprise the Behavior Regulation Index, and the Initiate, Working Memory, Plan/Organize, Task Monitor, and Organization of Materials scales comprise the Cognitive Regulation Index. For this study, we combined the 3 indexes into a total score as a singular measure of executive function for further analysis.
Childhood Depression
Children completed the Children’s Depression Inventory 2 (CDI-II)55, a modified version of the BDI for children. This evaluation was carried out once per child, at either their year 8.5, year 9 or year 10 visit. This variability in evaluation time was due to experimental limitations. The CDI-II consists of 28 questions scored on a 3-point Likert scale. The questions were divided into Emotional and Functional subscales and the sum of all the responses provide the scores for the corresponding subscales and total CDI-II score. A total score of 20 or above (13.3% of girls and 13.5% of boys) indicates clinical levels of depressive symptoms. The CDI-II is stable across ages, applicable for both depressed and normal children92,93 and has good test-retest reliability and high internal consistency94. In this study, the total score was utilized as the outcome measure.
Correlation Analysis
We calculated the simple Pearson’s correlation between both positive maternal mental health and negative maternal mental health with all the outcome measures. We also calculated the Pearson’s correlation between our candidate mediators, language and executive function, and childhood depression.
Serial Mediation Analysis
We carried out mediation analysis using the Lavaan Toolbox95 in R. Specifically, we built a serial mediation model, with language at 2 years as the first mediator, executive function at 7 years as the second mediator, and childhood depressive symptoms at 8.5–10 years as the outcome. Prenatal negative maternal mental health and positive maternal mental health were used as independent predictors in the serial mediation model. The bias-corrected significance of the mediation pathways were estimated using bootstrapping with 10,000 iterations.
MRI Data Acquisition
The MRI data was acquired via a 3-Tesla (Magnetom Prisma; Siemens, Germany). The mean age of participants at time of scan 7.46, with a standard deviation of 0.14 years. We acquired 3D T1-weighted MPRAGE images with the following imaging parameters: repetition time = 2000 ms, echo time = 2.08 ms, field of view = 192×192 mm2, matrix size = 192×192, number of slices = 192, slice thickness = 1 mm. We also acquired one run of resting-state fMRI data: 5.24 min (120 timepoints), repetition time = 2620 ms, echo time = 27 ms, flip angle = 90 degrees, 3 mm isotropic voxels, matrix = 64×64, 48 interleaved axial slices with no gap.
MRI Data Processing
Preliminary pre-processing of the data was performed using the default MNI pipeline in the CONN toolbox (release 20b)96 and SPM1297 in Matlab R2020a. For each session, all scans were co-registered and resampled to the first scan using b-spline interpolation in the SPM12 functional realignment and unwarp procedure98. Sinc interpolation using SPM1299 was applied for slice time correction. We flagged scan volumes with framewise displacement above 0.9 mm or global BOLD signal changes above 5 standard deviations as potential outliers for subsequent confound regression. Both functional and anatomical data were normalized into standard MNI space and segmented into grey matter, white matter and CSF tissue using SPM 12 unified segmentation and normalization procedure100 with 4th order spline interpolation. A kernel of 6 mm full width half maximum (FWHM) was used for gaussian smoothing of the functional data.
Denoising of the functional data was carried out in two steps. First, ordinary least squares regression was used to remove potential noise confounds (cerebral blood flow, subject motion parameters, identified outliers and constant and first-order linear session effects)101. Next, temporal band-pass filtering was used to remove signals below 0.008 Hz and above 0.09 Hz102. Subjects with greater than 6 mm maximum motion, 0.6 mm mean motion or more than 20% outlier scans were excluded from further analysis.
Functional Connectivity Matrices Analysis
Cortical regions were parcellated into 17 functional networks using the atlas of Yeo et al.103 which were used as our regions of interest (ROIs). For each ROI, we calculated the average ROI timeseries by averaging the time course across all voxels in the ROI. We then calculated Fisher-transformed bivariate correlation coefficients between the timeseries for each pair of ROIs to generate the functional connectivity matrix for that subject. Next, we calculated the Pearson’s correlation between each element of the functional connectivity matrix with the corresponding negative maternal mental health and positive maternal mental health total score across all subjects to obtain a matrix of the correlations of functional connectivity fluctuations with negative maternal mental health and positive maternal mental health respectively.
Validation using an independent cohort
A partial validation of the study was carried out using the Singapore Preconception Study of Long-Term Maternal and Child Outcomes (S-PRESTO) cohort56 as an independent replication cohort. As the children in the SPRESTO study were all below 7 years old at the time of this study, the CDI-II was not administered for them and we were only able to test the correlations between positive and negative maternal mental health and the respective mediators. Due to limited overlap of measures used in both cohorts, we utilized slightly different measures of language and executive function. However, this could also be viewed as a strength as a replication would suggest the observed effects arise from the underlying language and executive function abilities and is not instrument dependent.
We applied the same bifactor structure to the three maternal mental health questionnaires obtained from the S-PRESTO mothers at pregnancy week 26 to extract out the equivalent general negative mental health factor and positive maternal mental health factors. We used the Peabody Picture Vocabulary Test, 5th Ed (PPVT-5)104 as the measure of language ability. The Behaviour Rating Inventory of Executive Function, Preschool Version (BRIEF-P)105 was used as the measure of executive function. Both tests were administered during the same visit, with the child’s age at test between 4.5 to 6.3 years. Due to the wide range of ages, we utilized the standardized scores for both measures to adjust for age-related changes.