https://orcid.org/0000-0003-1950-2345
Figures
Abstract
Age and sex differences in brain metabolite concentrations in early life are not well understood. We examined the associations of age and sex with brain metabolite levels in healthy neonates, and investigated the associations between neonatal brain metabolite concentrations and developmental outcomes. Forty-one infants (36–42 gestational weeks at birth; 39% female) of predominantly Hispanic/Latina mothers (mean 18 years of age) underwent MRI scanning approximately two weeks after birth. Multiplanar chemical shift imaging was used to obtain voxel-wise maps of N-acetylaspartate (NAA), creatine, and choline concentrations across the brain. The Bayley Scales of Infant and Toddler Development, a measure of cognitive, language, and motor skills, and mobile conjugate reinforcement paradigm, a measure of learning and memory, were administered at 4 months of age. Findings indicated that postmenstrual age correlated positively with NAA concentrations in multiple subcortical and white matter regions. Creatine and choline concentrations showed similar but less pronounced age related increases. Females compared with males had higher metabolite levels in white matter and subcortical gray matter. Neonatal NAA concentrations were positively associated with learning and negatively associated with memory at 4 months. Age-related increases in NAA, creatine, and choline suggest rapid development of neuronal viability, cellular energy metabolism, and cell membrane turnover, respectively, during early life. Females may undergo earlier and more rapid regional developmental increases in the density of viable neurons compared to males.
Citation: Merz EC, Monk C, Bansal R, Sawardekar S, Lee S, Feng T, et al. (2020) Neonatal brain metabolite concentrations: Associations with age, sex, and developmental outcomes. PLoS ONE 15(12): e0243255. https://doi.org/10.1371/journal.pone.0243255
Editor: Emma Duerden, Western University, CANADA
Received: August 7, 2020; Accepted: November 17, 2020; Published: December 17, 2020
Copyright: © 2020 Merz et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: The DOI to our publicly available data deposit is https://doi.gin.g-node.org/10.12751/g-node.k2o9ik/.
Funding: Funding for this research was provided by a National Institute of Health (NIH) grant (R01MH093677-01) awarded to C. Monk and B. Peterson. The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH or other funders. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Brain development during the perinatal period and in the first few years of life is increasingly being studied. A variety of magnetic resonance imaging (MRI) techniques have been used to gain insight into structural and functional brain development in the first years of life [1]. Developmental changes in brain metabolite concentrations very early in life are not well understood, however. Magnetic resonance spectroscopy (MRS) is a non-invasive imaging technique used to measure in vivo concentrations of brain metabolites, including N-acetylaspartate (NAA), an index of the density of viable neurons [2]. Although MRS has been applied to the study of the neonatal brain (S1 Table), the vast majority of these studies have focused primarily on preterm or clinical samples, such as those with traumatic brain injury or hypoxia-ischemia. In addition, few studies have examined sex differences in neurochemical metabolites in the brain during early life, though animal models have revealed sex differences in a range of early neurodevelopmental processes [3], and sexual dimorphisms in brain structure and function have been found in infancy [1, 4, 5].
The vast majority of MRS studies on neonatal populations have used single-voxel spectroscopy. In comparison to single-voxel spectroscopy, multi-voxel spectroscopic techniques, such as multiplanar chemical shift imaging (MPCSI) [6–8], have many advantages, including higher spatial resolution, fewer partial volume effects because more voxels are positioned in a single tissue type, and a larger total coverage area that permits regional localization of the effects of interest [9, 10]. The main goal of this study was to use MPCSI to examine age correlates and sex differences in metabolite concentrations across the brain in healthy newborns.
NAA, creatine, and choline
MRS studies of NAA, creatine (Cr), and choline (Cho) yield important information about neurochemical processes. NAA, present primarily in neurons [2, 11], contributes to signaling between neurons and oligodendrocytes and participates in myelin synthesis by oligodendrocytes [12]. Cr is involved in cellular energy metabolism and storage, and Cho is a marker of the structural integrity and turnover of cell membranes [2].
Age correlates of neonatal brain metabolite concentrations
NAA concentrations during the neonatal period have been reported to increase significantly with age in both white and gray matter, and in multiple brain regions that include the thalamus, basal ganglia, hippocampus, frontal cortex, and cerebellum (S1 Table). In a study of preterm neonates, NAA and Cr concentrations in the basal ganglia, centrum semiovale, and cerebellum increased with post-conception age (from 30–43 weeks), with NAA displaying the greatest increases and no significant age correlates detected for Cho [13]. One large study of a patient radiological database reported dramatic postconceptional increases in NAA and Cr concentrations with age in several white matter and cortical and subcortical gray matter regions in the first three postnatal months [14]. No prior studies have used high-spatial-resolution, whole-brain spectroscopic imaging techniques such as MPCSI to investigate age correlates of brain metabolite concentrations in healthy newborns.
Sex differences in neonatal brain metabolites
Animal studies have demonstrated sex differences in a variety of neurodevelopmental processes, ranging from the rate of neurogenesis to synaptic physiology [3]. Many of these differences originate prenatally, when exposure to gonadal hormones differs across fetuses [15]. Although several studies have reported sexual dimorphisms in brain structure in infancy [1, 16], sex differences in brain metabolite concentrations have not been assessed extensively in neonates or infants [13, 14]. One single-voxel MRS study of mostly preterm neonates reported no significant sex differences in NAA, Cr, or Cho in the basal ganglia, centrum semiovale, or cerebellum [13]. To date, the extent to which brain metabolite concentrations differ by sex in healthy neonates is largely unknown.
Neonatal brain metabolites and later developmental outcomes
Higher neonatal NAA concentrations have been associated with more optimal developmental outcomes in clinical populations (S2 Table) [17–19]. In a longitudinal study of preterm neonates, steeper increases in the NAA/Cho ratio in basal ganglia nuclei and white matter regions between 28 and 40 weeks postmenstrual age (PMA) were associated with better subsequent motor, language, and cognitive performance at 18 months corrected age [20]. Other studies, however, have not found significant associations (S2 Table) [21, 22]. Collectively, metabolite research to date has been performed primarily on infants with significant clinical diagnoses that affect the brain. Given that findings from these studies may not generalize to healthy infants, research is needed on typically developing neonates. Furthermore, most previous studies have reported metabolite ratios (e.g., NAA/Cr, NAA/Cho), which are difficult to interpret from a developmental perspective, as both the numerator and denominator of a ratio can vary with age [2, 23].
Moreover, single-voxel spectroscopy has been used in the vast majority of studies examining associations between age and neonatal brain metabolite concentrations (S1 Table), even though multi-voxel spectroscopy methods, such as MPCSI, have several advantages over single-voxel spectroscopy [10]. Single-voxel spectroscopy does not simultaneously measure metabolite levels across a large number of brain regions, and thus the regional specificity of the metabolic findings in neonates (e.g., anterior vs. parietal, temporal, or subcortical) and their tissue specificity (i.e., gray matter [GM] vs. white matter [WM]) is difficult to discern [24]. Multi-voxel spectroscopy methods, on the other hand, can be used to measure metabolites in contiguous voxels across most of the brain [6] and at a higher spatial resolution than is possible using single-voxel spectroscopy [6–8].
Current study
Many cognitive, emotional, and behavioral problems are thought to have their origins early in development, and many have prominent sex-specific differences in prevalence and symptom severity [25]; contributors to adverse health outcomes likely are represented in brain metabolite concentrations in early infancy [26–29]. Therefore, mapping the metabolic correlates of age and sex in the infant brain is important to understand the determinants of later health and illness. The goals of this study were to examine age and sex differences in NAA, Cr, and Cho levels in the healthy newborn brain and to investigate associations of neonatal brain metabolite levels with 4-month cognitive, language, and motor skills. MPCSI data were acquired approximately two weeks after birth in 41 infants (39% female; 36–42 gestational weeks at birth) of healthy adolescent or young adult mothers (14–20 years). The Bayley Scales of Infant and Toddler Development, Third Edition (Bayley-III) [30] and the mobile conjugate reinforcement paradigm, an assessment of learning and memory [31], were administered at 4 months of age. We hypothesized that postmenstrual age (PMA) at the time of MRI scanning would correlate positively with NAA concentrations across widespread regions of the brain. We also expected to detect sex differences in NAA concentrations, although the predicted direction and locations could not be specified [1]. Finally, we explored whether NAA concentrations shortly after birth were associated with developmental outcomes at 4 months of age. Cr and Cho were included as exploratory analyses that would aid interpretation of NAA findings.
Methods
Participants
Data used in this study came from a larger study of healthy, nulliparous pregnant adolescents and young adults who were actively engaged in prenatal care. They were recruited through the Department of Obstetrics and Gynecology at Columbia University Irving Medical Center (CUIMC) and Weill Cornell Medical College, and flyers posted in the local vicinity. Interested individuals were contacted by phone and screened for eligibility. Exclusionary criteria were smoking, use of recreational drugs, lack of fluency in English, and use of nitrates, steroids, beta blockers, triptans, or psychiatric medications.
Of the 154 mothers who enrolled in the study and were offered newborn MRI scans and infancy follow-up sessions, 82 newborns underwent MRI scans, and MPCSI data were acquired in 41 newborns. This pulse sequence was subsequently discontinued at our location following an upgrade that did not support the MPCSI pulse sequence, thus precluding scanning more infants. Mothers (93% Hispanic/Latina) were 14–20 years of age (70% were 18–19 years old). Infants (39% female) were born at 36–42 gestational weeks with typical birth weights; they ranged in PMA from 39–47 weeks at the time of MRI scan (Table 1). At 4 months, 27 and 26 infants had usable data for the Bayley-III and mobile conjugate reinforcement task, respectively. There were no major differences between those with and without developmental data at 4 months (S1 Text).
Procedures
Mothers brought their infants for MRI scans as newborns and for developmental assessments at age 4 months. Infant medical records were used to ascertain gestational age at birth, birth weight, pregnancy and birth complications, and sex. Written informed consent was provided by all mothers, and all procedures were approved by the Institutional Review Board of the New York State Psychiatric Institute.
MRI acquisition
Non-sedated infants underwent MRI scanning on a 3T whole body MR scanner (GE Health Care, Milwaukee, WI) equipped with an 8–channel head coil. MPCSI data were acquired in 6 axial oblique slices parallel to the AC-PC line, with the second bottom-most slice containing the AC-PC plane. Parameters for the MPCSI sequence were: TR = 2800 ms, TE = 144 ms, spectral width = 2000 Hz, number of complex data points = 512, FOV = 24 cm, slice thickness = 10.0 mm, slice spacing = 2.0 mm, number of phase encoding steps = 24 x 24. Water signal was suppressed using the CHESS sequence. Lipid signal was suppressed by placing eight angulated saturation bands around the brain. MPCSI data were spatially registered to a template brain using a localizer image of high in-plane resolution in the same orientation and slice locations as the MPCSI data. Those data were acquired using a 2D spin echo sequence with the following parameters: TR = 300 ms, TE = 10 ms, FOV = 24 cm, slice thickness = 10.0 mm, spacing = 2.0 mm, acquisition matrix = 256 x 128, image zero-padded to 256 x 256. We also acquired T2-weighted images for normalizing MRS data into a common template space. Parameters for the 2D spin echo sequence used to acquire the T2-weighted image were as follows: TR = 10,000 ms, TE = 129 ms, FOV = 100 cm, slice thickness = 1.0 mm, spacing = 1.0 mm, acquisition matrix = 192 x 192.
We used a long echo time (TE = 144 ms) because it flattened the spectral baseline and minimized the contamination of metabolite signals from macromolecules, thereby improving the accuracy of spectral fitting for metabolite quantification. However, the long TE limited metabolite measurement to NAA, Cr, and Cho, because these metabolites are present in greater abundance than other metabolites in the brain [32].
MPCSI processing
Full MPCSI processing procedures are provided in the S1 Text. A brief summary of these procedures is provided here. Rigorous quality assurance procedures were employed to ensure that the data used in analyses were not corrupted by noise or participant motion. More specifically, MRS data quality was assured by reconstructing the data, assessing it for excess noise, and examining the spectrum in each voxel for baseline distortions, signal contamination by lipid signal from the scalp, incorrect placement of suppression bands, or broadening of line width. Voxels with any of these artifacts or lipid contamination were rejected and were not processed further.
We processed signal from each coil of the 8-channel head coil separately before combining their processed MRS signals to generate the spectroscopic images [33]. The combined signal was then loaded into the software 3DiCSI (3D Interactive Chemical Shift Imaging) [34] to identify MRS voxels within the brain and save spectral data for those voxels [6]. We used model-based spectral fitting to model the spectrum in each voxel with Voigtian curves to the peaks for NAA, Cr, and Cho metabolites (S1 Text). Background noise was calculated as the standard deviation of the real part of the complex data in the regions free of signal from metabolites. A spectroscopic image for each metabolite was then generated as the ratio of the peak area to the background noise for each MRS voxel within the brain, accounting for variations in receiver and transmitter gain. Fig 1 shows a representative spectrum in a voxel of the MPCSI dataset.
Principal resonances were fit for NAA plus N-acetyl-aspartyl-glutamate (2.0 parts per million [ppm]; green), Cr plus phosphocreatine (3.0 ppm; yellow), Ch compounds (3.2 ppm; red), as well as contaminating lipids (blue). The quality of the spectra and signal-to-noise are high. The orange line at the top of the figure is the residual after fitting metabolite peaks to the MR spectrum.
We corrected the spectroscopic images for partial-volume effects and for the MPCSI point-spread function. For each MPCSI voxel and metabolite, we used linear regression to estimate the concentration of that metabolite in brain tissue and CSF using the levels of that metabolite in neighboring voxels and the proportions of brain tissue and CSF in each voxel. During spatial normalization, MPCSI data for each participant were coregistered into the coordinate space of a T2-weighted image of a template brain. The template brain was selected from our cohort of newborns using procedures that ensured the template brain was morphologically the most representative brain in the cohort (S1 Text).
The MPCSI saturation bands applied to suppress lipid signal from the scalp were not as precisely shaped as the scalp, and they unavoidably suppressed metabolite signals from several portions of cortical gray matter. Consequently, metabolite measures for many participants were available only in voxels of white matter and deep gray matter nuclei. We therefore presented results at only those voxels that had usable data from 75% or more of the infants.
Measures
PMA at Scan PMA at the time of MRI scan is the sum of gestational age at birth and postnatal age at the time of scan [35]. Gestational age at birth was calculated based on ultrasound examinations and the date of the last reported menstrual cycle from the medical record.
Cognitive, Language, and Motor Skills The Bayley-III [30] was administered at postnatal age 4 months to assess cognitive, language, and motor skills (S1 Text). Raw and age-standardized scores (M = 100, SD = 15) were calculated; raw scores were used in analyses. The Bayley-III has acceptable psychometric properties and norms based on a racially- and socioeconomically-representative sample [30, 36], although problems with Bayley-III reliability and validity have been noted [37–39].
Learning and Memory At 4-months postpartum, mothers brought their infants to the lab on two consecutive days for participation in the mobile conjugate reinforcement paradigm [31, 40] (S1 Text). After placing infants in a specially set up crib, each videotaped 15-min session began with a 3-min non-reinforcement phase (baseline), followed by a 9-min reinforcement phase (three 3-min learning blocks), and a final 3-min non-reinforcement phase (immediate retention). During periods of reinforcement, when infants kicked, a mobile hanging above the crib bounced, whereas during periods of non-reinforcement, the mobile did not bounce as a result of their kicking. On Day 1, the initial 3-min period of non-reinforcement (baseline) provides a measure of the infant's baseline kick rate, and the final 3-min period of non-reinforcement provides a measure of the infant's immediate retention. On Day 2, kicking during the initial period of non-reinforcement (Day 2 baseline) reflects the infant's long-term (24-hr) retention of the contingency.
Trained coders used the videotapes to count the number of times infants kicked per minute. A second trained coder independently coded 20% of the sessions. Inter-rater reliability was high (Spearman rank correlation = 0.95; range: 0.88–0.99). Ratio scores were calculated by dividing the mean kick rate during each learning and retention block by the mean kick rate during Day 1 baseline [40].
Statistical analyses
Descriptive statistics were computed using SAS (version 9.4). Multiple linear regression was applied voxel-wise throughout the brain, using an in-house developed program, separately for each of the three metabolite values (NAA, Cho, and Cr) entered as the dependent variable and with PMA at scan and sex entered simultaneously as independent variables. Similarly, NAA, Cho, and Cr concentrations were regressed voxel-wise onto Bayley-III and mobile conjugate reinforcement paradigm scores, with PMA and sex included as covariates.
We used the Benjamini–Yekutieli procedure [41] for False Discovery Rate (FDR) to control for false positives within each statistical map [42–44]; p-values surviving FDR correction at a corrected p < 0.05 level were color-coded as statistical parametric maps on the T2-template. This FDR procedure is a rigorous statistical procedure that limits false discovery in findings at a specified rate, which in our analyses is 0.05. This procedure will ensure that no more than 5% of the reported findings are generated by chance due to noise in the data. By permitting a prespecified fraction of false positives in our findings, this procedure provides significantly greater statistical power to detect true findings in the data. Furthermore, FDR is adaptive and adjusts the threshold for significant findings depending on the noise in the data. In contrast, other alternative familywise error rate procedures, such as Holm-Bonferroni correction, control for false positives across the entire family of hypotheses tested at a specified significance level α. Although Bonferroni correction will ensure that the probability of finding one or more false positive is smaller than α across all hypotheses, this procedure is highly conservative [45] providing low statistical power and therefore too often failing to detect true positives. Thus we opted to balance false positive and false negative findings through use of FDR correction.
Results
Descriptive statistics for the Bayley-III and mobile conjugate reinforcement paradigm are presented in Table 2. The infants in this study demonstrated typical patterns of learning and memory performance on the mobile conjugate reinforcement paradigm (see also [40]).
Postmenstrual age correlates positively with newborn brain metabolite concentrations
Controlling for infant sex, advancing PMA at scan was significantly associated with higher NAA concentrations in white matter (superior longitudinal fasciculus), superior frontal gyrus, and subcortical gray matter (thalamus, basal ganglia) (Fig 2). Similar results were evident for Cr and, to a lesser extent, Cho. Representative scatterplots are presented in Figs 3 and 4. Effect sizes for PMA at scan across voxels are presented in S1 Fig.
Findings for creatine and choline were weaker but in the same direction and in overlapping brain regions as the findings for NAA. Anatomical reference slices are presented in the top row. Maps indicating the voxels for which at least 75% of the sample had usable data are presented in the bottom row. FDR correction for multiple comparisons was done with FDR at p < .05. Warm colors (red, orange) indicate significant positive associations and cool colors (shades of blue) indicate significant inverse associations, with higher degrees of warmth/coolness corresponding to smaller FDR-corrected p-values as shown in the color bar. SLF, superior longitudinal fasciculus; BG, basal ganglia; Th, thalamus; SFG, superior frontal gyrus; Cd, caudate.
Postmenstrual age at scan (in weeks) was significantly positively associated with NAA concentrations in the (a) left and (c) right basal ganglia after controlling for sex. FDR-corrected p-values are shown in the color bar. Metabolite data were sampled from the regions indicated by the white dots (b). NAA is shown in arbitrary units (a.u.).
Postmenstrual age at scan (in weeks) was significantly positively associated with NAA concentrations in the (a) left caudate and (c) right superior longitudinal fasciculus after controlling for sex. FDR-corrected p-values are shown in the color bar. Metabolite data in the scatterplots were sampled from the regions indicated by the white dots (b). NAA is shown in arbitrary units (a.u.).
Female newborns have higher brain metabolite concentrations than male newborns
Females compared to males had greater NAA concentrations in white matter (corpus callosum), subcortical gray matter (thalamus), insula, and superior frontal gyrus. Females compared to males had greater Cho concentrations in the thalamus, insula, superior frontal gyrus, and premotor regions. Females compared to males had greater Cr concentrations in the thalamus, corpus callosum, and insula (Fig 5). Bar graphs for representative voxels are presented in Fig 6. Effect sizes for sex across voxels are presented in S2 Fig.
Anatomical reference slices are presented in the top row. Maps indicating the voxels for which at least 75% of the sample had usable data are presented in the bottom row. FDR correction for multiple comparisons was done with FDR at p < .05. Cool colors (shades of blue) indicate significantly higher metabolite concentrations in females compared to males. The magnitude of FDR-corrected p-values is color-coded as shown in the color bar. Th, thalamus; SFG, superior frontal gyrus; Ins, insula; PM, premotor cortex; CC, corpus callosum.
Females had higher NAA concentrations than males in the (a) insula and (c) superior frontal gyrus. FDR-corrected p-values are shown in the color bar. The error bars represent standard error (SE). Metabolite data in the scatterplots were sampled from the regions indicated by the white dots (b and d). NAA is shown in arbitrary units (a.u.).
Newborn NAA concentrations correlate with 4-month learning and memory
Neonatal NAA concentrations correlated significantly with performance on the mobile conjugate reinforcement paradigm at 4 months of age. Higher NAA concentrations in the superior longitudinal fasciculus (SLF), corpus callosum, insula, and thalamus were associated with greater learning. In contrast, lower NAA concentrations in the corpus callosum, SLF, insula, and basal ganglia were associated with greater immediate and long-term retention (Fig 7). Scatterplots for representative voxels are presented in Fig 8. Effect sizes for scores on this task across voxels are presented in S3 Fig. Newborn brain metabolite concentrations were not significantly associated with performance on the Bayley-III at 4 months.
NAA concentrations in the SLF and splenium of the corpus callosum were significantly positively associated with learning (b and c). NAA concentrations in the SLF and genu of the corpus callosum were significantly inversely associated with immediate and long-term retention (d and e). P-values are FDR-corrected and color coded if corrected p is < .05. Warm colors (red, orange) indicate significant positive associations and cool colors (shades of blue) indicate significant inverse associations, with higher degrees of warmth/coolness corresponding to smaller FDR-corrected p-values as shown in the color bar. SLF, superior longitudinal fasciculus; CCs, splenium of corpus callosum; CCg, genu of corpus callosum; Insp, posterior insula; Th, thalamus; Insa, anterior insula; BG, basal ganglia; MB, midbrain.
NAA concentrations in the SLF (b and d) were significantly positively associated with learning skills (a and c), while controlling for PMA and sex. NAA concentrations in a different portion of the SLF (f) were significantly inversely associated with long-term retention (e). FDR-corrected p-values are shown in the color bar. Metabolite data in the scatterplots were sampled from the regions indicated by the white dots. NAA is shown in arbitrary units (a.u.).
Discussion
The goals of this study were to examine postmenstrual age (PMA)- and sex-related variability in brain metabolite concentrations in healthy neonates and the associations of neonatal brain metabolite levels with 4-month developmental outcomes. Previous MRS studies with infants have focused primarily on preterm or clinical samples and have not examined sex differences. In addition, previous research has largely used single-voxel spectroscopy techniques (S1 Table), whereas multi-voxel spectroscopic techniques, such as MPCSI, have higher spatial resolution and can shed light on regional specificity of the associations [6–8]. This study is the first to use MPCSI to assess concentrations of NAA, Cr, and Cho in the healthy newborn brain.
Age correlates of neonatal NAA concentrations
As hypothesized, older PMA was significantly associated with higher NAA concentrations in anterior and posterior deep white matter and subcortical gray matter. NAA is synthesized in the mitochondria of neurons and stored largely in neuronal bodies and their axons. A portion of its signal may be attributable to N-acetyl-aspartyl-glutamate (NAAG), a prominent neuropeptide that regulates glutamate and dopamine release [46, 47]. Increasing NAA concentrations reflect greater mitochondrial energy metabolism, demands for lipid synthesis, or NAAG [12, 48]. Our findings of strong positive correlations of NAA concentrations with increasing PMA suggest rapid developmental increases in the density of neuronal tissue in these regions in early postnatal life [2].
PMA-related increases in NAA were found in the basal ganglia and thalamus, which undergo rapid morphological growth in early life [1, 49]. Sensorimotor portions of the basal ganglia, together with the sensorimotor cortex, drive motor learning and control. The thalamus relays basal ganglia output to the cortex and mediates information flow between cortical circuits. PMA-related increases in NAA in these structures during the perinatal period are consistent with the rapid rates of synaptogenesis and dendritic growth that characterize this developmental period in these nuclei (Gilmore et al., 2018).
PMA-related increases in NAA were also found in higher-order association areas and associative fiber bundles (e.g., frontal white matter, SLF). The SLF, which connects dorsal frontal with inferior and superior parietal cortices, is a component of neural networks that subserve numerous cognitive functions, including attention and self-regulation [50]. Increases in NAA concentrations may contribute to the rapid early developmental increases in organization (fractional anisotropy [FA]) of white matter tracts, including the SLF [1, 6, 51]. PMA-related increases in NAA in white matter are consistent with rapid myelination during the first year of life [1, 52, 53] and the lipid synthesis it demands [12, 54].
Sex differences in neonatal NAA concentrations
Females compared to males had higher NAA concentrations in cerebral white matter (insula, superior frontal gyrus, corpus callosum), premotor cortex, and subcortical gray matter (thalamus). These findings are consistent with an emerging literature showing infant sex differences in regional gray and white matter volumes [1, 4, 5, 49, 55] and white matter microstructure [1, 56, 57], and with previous theoretical and empirical work suggesting that the brain may develop earlier and more rapidly in females than in males [16, 58, 59]. For example, female newborns have a higher rate of myelination in the corpus callosum between 3 and 60 months [56], higher FA values across widespread white matter regions [57], and larger volume of the motor cortex [4].
Given that the sex differences we found were present shortly after birth, they are likely the product of hormonal and genetic influences on the fetal and newborn brain, rather than of differential environments and experiences in postnatal life [3, 60]. Differential exposure to gonadal steroid hormones and differential spatial and temporal profiles of receptors for those hormones in male and female fetuses produces sexual dimorphisms in brain structure and function [15, 61]. It has been suggested that testosterone may slow neural development, and some animal research has supported this conjecture [58]. These sex differences in prenatal brain development may lead to sex differences in how infants elicit and respond to stimulation from their environment during postnatal life [62–65].
Associations of neonatal NAA concentrations with 4-month developmental outcomes
Newborn NAA concentrations were significantly associated with 4-month developmental outcomes as measured by the mobile conjugate reinforcement paradigm. Greater NAA concentrations in the SLF and corpus callosum were associated with higher learning scores on the task. These findings are consistent with previous work in preterm samples that has associated higher neonatal NAA levels with higher scores on developmental measures [17–19, 66], although other studies have not found these associations (S2 Table) [21, 22, 67]. In contrast, lower NAA concentrations in these locations were associated with higher retention scores on the task. It is possible that these lower levels may signify a slower rate of development or more prolonged maturation that supports greater memory task performance, consistent with previous longitudinal MRI work on the neural underpinnings of general cognitive development [68].
Age correlates of neonatal creatine and choline concentrations
Cr concentrations also increased with advancing PMA at scan, albeit less dramatically and in fewer brain regions compared with NAA. Findings for Cho were similar but even more attenuated. This pattern of findings across metabolites is consistent with previous studies in which age correlates for NAA and Cr in neonates have been detected more frequently than for Cho (S1 Table) [13, 14, 69, 70]. Age-related increases in Cr may represent an increase in energy metabolism or the need for increased energy reserve with advancing age, whereas increases in Cho likely reflect an increased rate of cell membrane synthesis or turnover. The common pattern of findings, in which age correlates were in the same direction and overlapping regions may also be driven by a developmental increase in mass and density of neuronal tissue, which will necessarily entail an increase in all metabolites that support neuronal processes, including cell energetics and membrane turnover.
Limitations
This study had several limitations. Interpretation of age correlates from cross-sectional data should be regarded with caution [71], as definitive identification of developmental trajectories requires repeated measures over time. Metabolite concentrations were not available for most cortical gray matter over the surface convexities, because MPCSI saturation bands applied to suppress lipid signal from the scalp were not as precisely shaped as the scalp, and they therefore unavoidably suppressed metabolite signals from cortical gray matter. Attrition in follow-up also limited our sample size and statistical power for testing correlations of metabolite concentrations with 4-month developmental outcomes. In addition, all infants were born to adolescent or young adult mothers, who are often considered an at-risk population, though screening and exclusionary criteria ensured a healthy, reasonably representative sample of newborns for women in this age group. Nonetheless, findings from this study may not be generalizable to the children of older mothers. In addition, participant motion during image acquisition is a common concern in neuroimaging studies and could have affected the results reported here. However, the infants in this study were asleep throughout the acquisition of the MPCSI data, and their heads were immobilized. These acquisition procedures coupled with the rigorous quality assurance procedures employed in this study make it unlikely that participant motion during image acquisition corrupted the analyses or results of this study.
Results from this study shed light on typical patterns of early metabolic development across the brain, providing the basis for identifying regional developmental deviations that may signal risk for neurodevelopmental and psychiatric disorders. Sex differences in early metabolite concentrations, with females tending to have higher regional concentrations than males, may relate to the well-documented sex differences in risk for neurodevelopmental and psychiatric disorders.
Supporting information
S1 Table. Studies that investigated associations of age with neonatal brain N-acetylaspartate (NAA), creatine, and choline concentrations.
https://doi.org/10.1371/journal.pone.0243255.s001
(DOCX)
S2 Table. Studies that investigated associations between neonatal brain metabolite concentrations and later developmental outcomes.
https://doi.org/10.1371/journal.pone.0243255.s002
(DOCX)
S1 Fig. Strength of associations between postmenstrual age (PMA) at scan and brain metabolite concentrations.
The magnitude of association is displayed in terms of beta values from multiple linear regression applied voxel-wise throughout the brain separately for each of the three metabolite values (NAA, Ch, and Cr) entered as the dependent variable and PMA at scan and sex entered simultaneously as independent variables. Anatomical reference slices are presented in the top row. Maps indicating the voxels for which at least 75% of the sample had usable data are presented in the bottom row.
https://doi.org/10.1371/journal.pone.0243255.s003
(TIF)
S2 Fig. The magnitude of sex differences in brain metabolite concentrations in neonates.
The magnitude of the association is displayed in terms of beta values from multiple linear regression applied voxel-wise throughout the brain separately for each of the three metabolite values (NAA, Ch, and Cr) entered as the dependent variable and PMA at scan and sex entered simultaneously as independent variables. Anatomical reference slices are presented in the top row. Maps indicating the voxels for which at least 75% of the sample had usable data are presented in the bottom row.
https://doi.org/10.1371/journal.pone.0243255.s004
(TIF)
S3 Fig. The magnitude of associations between neonatal NAA concentrations and performance on the mobile conjugate reinforcement paradigm at 4 months.
Beta values are displayed from multiple linear regression applied voxel-wise throughout the brain for NAA entered as the dependent variable and ratio score as the independent variable, with PMA and sex included as covariates. Brain maps are presented for ratio scores on learning block 2, learning block 3, and the long-term retention block (Day 2 baseline) for the mobile conjugate reinforcement paradigm.
https://doi.org/10.1371/journal.pone.0243255.s005
(TIF)
S1 Text. Supporting information about the methods of the study.
https://doi.org/10.1371/journal.pone.0243255.s006
(DOCX)
Acknowledgments
The authors are grateful to the families who participated in this study for making it possible. The authors also thank their dedicated research assistants, Laura Kurzius, Willa Marquis, Colleen Doyle, Ashley Rainford, Carlo Nati, Zachary Toth, and Kirwan Walsh, for help with enrolling families and collecting data.
References
- 1. Gilmore JH, Knickmeyer RC, Gao W. Imaging structural and functional brain development in early childhood. Nat Rev Neurosci. 2018 Feb 16;19(3):123–37. pmid:29449712
- 2. Rae CD. A guide to the metabolic pathways and function of metabolites observed in human brain 1H magnetic resonance spectra. Neurochem Res. 2014 Jan;39(1):1–36. pmid:24258018
- 3. McCarthy MM, Nugent BM, Lenz KM. Neuroimmunology and neuroepigenetics in the establishment of sex differences in the brain. Nat Rev Neurosci. 2017 Aug;18(8):471–84. pmid:28638119
- 4. Knickmeyer RC, Wang J, Zhu H, Geng X, Woolson S, Hamer RM, et al. Impact of Sex and Gonadal Steroids on Neonatal Brain Structure. Cereb Cortex. 2014 Oct;24(10):2721–31. pmid:23689636
- 5. Knickmeyer RC, Xia K, Lu Z, Ahn M, Jha SC, Zou F, et al. Impact of Demographic and Obstetric Factors on Infant Brain Volumes: A Population Neuroscience Study. Cereb Cortex. 2017 Dec 1;27(12):5616–25. pmid:27797836
- 6. Hao X, Xu D, Bansal R, Dong Z, Liu J, Wang Z, et al. Multimodal Magnetic Resonance Imaging: The Coordinated Use of Multiple, Mutually Informative Probes to Understand Brain Structure and Function. Hum Brain Mapp. 2013 Feb;34(2):253–71. pmid:22076792
- 7. Nelson MB, O’Neil SH, Wisnowski JL, Hart D, Sawardekar S, Rauh V, et al. Maturation of Brain Microstructure and Metabolism Associates with Increased Capacity for Self-Regulation during the Transition from Childhood to Adolescence. J Neurosci. 2019 Oct 16;39(42):8362–75. pmid:31444243
- 8. O’Neill J, Bansal R, Goh S, Rodie M, Sawardekar S, Peterson BS. Parsing the Heterogeneity of Brain Metabolic Disturbances in Autism Spectrum Disorder. Biol Psychiatry. 2020 Jan 15;87(2):174–84. pmid:31427037
- 9. Meyerhoff DJ, Durazzo TC, Ende G. Chronic alcohol consumption, abstinence and relapse: brain proton magnetic resonance spectroscopy studies in animals and humans. Curr Top Behav Neurosci. 2013;13:511–40. pmid:21688208
- 10. Zhu H, Barker PB. MR Spectroscopy and Spectroscopic Imaging of the Brain. Methods Mol Biol Clifton NJ. 2011;711:203–26. pmid:21279603
- 11. Baslow MH. N-acetylaspartate in the vertebrate brain: metabolism and function. Neurochem Res. 2003 Jun;28(6):941–53. pmid:12718449
- 12. Moffett JR, Ross B, Arun P, Madhavarao CN, Namboodiri MAA. N-Acetylaspartate in the CNS: From Neurodiagnostics to Neurobiology. Prog Neurobiol. 2007 Feb;81(2):89–131. pmid:17275978
- 13. Tomiyasu M, Aida N, Endo M, Shibasaki J, Nozawa K, Shimizu E, et al. Neonatal Brain Metabolite Concentrations: An In Vivo Magnetic Resonance Spectroscopy Study with a Clinical MR System at 3 Tesla. PLoS ONE [Internet]. 2013 Nov 28 [cited 2018 Sep 17];8(11). Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3842974/
- 14. Blüml S, Wisnowski JL, Nelson MD, Paquette L, Gilles FH, Kinney HC, et al. Metabolic Maturation of the Human Brain From Birth Through Adolescence: Insights From In Vivo Magnetic Resonance Spectroscopy. Cereb Cortex. 2013 Dec;23(12):2944–55. pmid:22952278
- 15. Lombardo MV, Ashwin E, Auyeung B, Chakrabarti B, Taylor K, Hackett G, et al. Fetal Testosterone Influences Sexually Dimorphic Gray Matter in the Human Brain. J Neurosci. 2012 Jan 11;32(2):674–80. pmid:22238103
- 16. Lebel C, Deoni S. The development of brain white matter microstructure. NeuroImage. 2018 15;182:207–18. pmid:29305910
- 17. Bapat R, Narayana PA, Zhou Y, Parikh NA. Magnetic Resonance Spectroscopy at Term-Equivalent Age in Extremely Preterm Infants: Association With Cognitive and Language Development. Pediatr Neurol. 2014 Jul;51(1):53–9. pmid:24938140
- 18. Kendall GS, Melbourne A, Johnson S, Price D, Bainbridge A, Gunny R, et al. White matter NAA/Cho and Cho/Cr ratios at MR spectroscopy are predictive of motor outcome in preterm infants. Radiology. 2014 Apr;271(1):230–8. pmid:24475798
- 19. Lally PJ, Montaldo P, Oliveira V, Soe A, Swamy R, Bassett P, et al. Magnetic resonance spectroscopy assessment of brain injury after moderate hypothermia in neonatal encephalopathy: a prospective multicentre cohort study. Lancet Neurol. 2019 Jan;18(1):35–45. pmid:30447969
- 20. Chau V, Synnes A, Grunau RE, Poskitt KJ, Brant R, Miller SP. Abnormal brain maturation in preterm neonates associated with adverse developmental outcomes. Neurology. 2013 Dec 10;81(24):2082–9. pmid:24212394
- 21. Augustine EM, Spielman DM, Barnes PD, Sutcliffe TL, Dermon JD, Mirmiran M, et al. Can magnetic resonance spectroscopy predict neurodevelopmental outcome in very low birth weight preterm infants? J Perinatol. 2008 Sep;28(9):611–8. pmid:18615089
- 22. Taylor MJ, Vandewouw MM, Young JM, Card D, Sled JG, Shroff MM, et al. Magnetic resonance spectroscopy in very preterm-born children at 4 years of age: developmental course from birth and outcomes. Neuroradiology. 2018 Oct;60(10):1063–73. pmid:30105622
- 23. Karl A, Werner A. The use of proton magnetic resonance spectroscopy in PTSD research—meta-analyses of findings and methodological review. Neurosci Biobehav Rev. 2010 Jan;34(1):7–22. pmid:19559046
- 24. Durazzo TC, Meyerhoff DJ, Mon A, Abé C, Gazdzinski S, Murray DE. Chronic Cigarette Smoking in Healthy Middle-Aged Individuals Is Associated With Decreased Regional Brain N-acetylaspartate and Glutamate Levels. Biol Psychiatry. 2016 Mar 15;79(6):481–8. pmid:25979621
- 25. Ruigrok ANV, Salimi-Khorshidi G, Lai M-C, Baron-Cohen S, Lombardo MV, Tait RJ, et al. A meta-analysis of sex differences in human brain structure. Neurosci Biobehav Rev. 2014 Feb;39(100):34–50. pmid:24374381
- 26. Allen PJ. Creatine metabolism and psychiatric disorders: Does creatine supplementation have therapeutic value? Neurosci Biobehav Rev. 2012 May;36(5):1442–62. pmid:22465051
- 27. Cohen-Gilbert JE, Jensen JE, Silveri MM. Contributions of Magnetic Resonance Spectroscopy to Understanding Development: Potential Applications in the study of Adolescent Alcohol Use and Abuse. Dev Psychopathol. 2014 May;26(2):405–23. pmid:24621605
- 28. Hodes GE, Epperson CN. Sex Differences in Vulnerability and Resilience to Stress Across the Life Span. Biol Psychiatry. 2019 May 7;
- 29. Rubinow DR, Schmidt PJ. Sex differences and the neurobiology of affective disorders. Neuropsychopharmacology. 2019;44(1):111–28. pmid:30061743
- 30.
Bayley N. Bayley Scales of Infant and Toddler Development-Third Edition. San Antonio, TX: Psychological Corporation; 2006.
- 31. Rovee-Collier C, Hartshorn K, DiRubbo M. Long-term maintenance of infant memory. Dev Psychobiol. 1999 Sep;35(2):91–102. pmid:10461123
- 32. Goh S, Dong Z, Zhang Y, DiMauro S, Peterson BS. Mitochondrial Dysfunction as a Neurobiological Subtype of Autism Spectrum Disorder: Evidence From Brain Imaging. JAMA Psychiatry. 2014 Jun 1;71(6):665. pmid:24718932
- 33. Dong Z, Peterson B. The rapid and automatic combination of proton MRSI data using multi-channel coils without water suppression. Magn Reson Imaging. 2007 Oct;25(8):1148–54. pmid:17905247
- 34. Zhao Q, Patriotis P, Arias-Mendoza F, Stoyanova R, Brown T. 3D Interactive Chemical Shift Imaging: A Comprehensive Software Program for Data Analysis and Quantification. In 2007.
- 35. Engle WA, American Academy of Pediatrics Committee on Fetus and Newborn. Age terminology during the perinatal period. Pediatrics. 2004 Nov;114(5):1362–4. pmid:15520122
- 36. Albers CA, Grieve AJ. Review of Bayley Scales of Infant and Toddler Development—Third Edition. J Psychoeduc Assess. 2007;25(2):180–90.
- 37. Anderson PJ, De Luca CR, Hutchinson E, Roberts G, Doyle LW, Victorian Infant Collaborative Group. Underestimation of developmental delay by the new Bayley-III Scale. Arch Pediatr Adolesc Med. 2010 Apr;164(4):352–6. pmid:20368488
- 38. Anderson PJ, Burnett A. Assessing developmental delay in early childhood—concerns with the Bayley-III scales. Clin Neuropsychol. 2017 Feb;31(2):371–81. pmid:27687612
- 39. Moore T, Johnson S, Haider S, Hennessy E, Marlow N. Relationship between test scores using the second and third editions of the Bayley Scales in extremely preterm children. J Pediatr. 2012 Apr;160(4):553–8. pmid:22048046
- 40. Merz EC, McDonough L, Huang YL, Foss S, Werner E, Monk C. The mobile conjugate reinforcement paradigm in a lab setting. Dev Psychobiol. 2017 Jul;59(5):668–72. pmid:28436585
- 41. Yekutieli D, Benjamini Y. Resampling-based false discovery rate controlling multiple test procedures for correlated test statistics. J Stat Plan Inference. 1999 Dec 1;82(1):171–96.
- 42. Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J R Stat Soc Ser B Methodol. 1995;57(1):289–300.
- 43. Benjamini Y. Discovering the false discovery rate. J R Stat Soc Ser B Stat Methodol. 2010;72(4):405–16.
- 44. Chumbley J, Worsley K, Flandin G, Friston K. Topological FDR for neuroimaging. NeuroImage. 2010 Feb 15;49(4):3057–64. pmid:19944173
- 45. Chen S-Y, Feng Z, Yi X. A general introduction to adjustment for multiple comparisons. J Thorac Dis. 2017 Jun;9(6):1725–9. pmid:28740688
- 46. Pouwels PJ, Frahm J. Differential distribution of NAA and NAAG in human brain as determined by quantitative localized proton MRS. NMR Biomed. 1997 Apr;10(2):73–8. pmid:9267864
- 47. Pouwels PJ, Frahm J. Regional metabolite concentrations in human brain as determined by quantitative localized proton MRS. Magn Reson Med. 1998 Jan;39(1):53–60. pmid:9438437
- 48. Amaral A, Hadera MG, Kotter M, Sonnewald U. Oligodendrocytes Do Not Export NAA-Derived Aspartate In Vitro. Neurochem Res. 2017;42(3):827–37. pmid:27394419
- 49. Qiu A, Fortier MV, Bai J, Zhang X, Chong Y-S, Kwek K, et al. Morphology and microstructure of subcortical structures at birth: a large-scale Asian neonatal neuroimaging study. NeuroImage. 2013 Jan 15;65:315–23. pmid:23000785
- 50. Goh S, Bansal R, Xu D, Hao X, Liu J, Peterson BS. Neuroanatomical correlates of intellectual ability across the life span. Dev Cogn Neurosci. 2011 Jul;1(3):305–12. pmid:22436512
- 51. Wijtenburg SA, McGuire SA, Rowland LM, Sherman PM, Lancaster JL, Tate DF, et al. Relationship between fractional anisotropy of cerebral white matter and metabolite concentrations measured using 1H magnetic resonance spectroscopy in healthy adults. NeuroImage. 2013 Feb 1;0:161–8.
- 52. Dubois J, Dehaene-Lambertz G, Kulikova S, Poupon C, Hüppi PS, Hertz-Pannier L. The early development of brain white matter: a review of imaging studies in fetuses, newborns and infants. Neuroscience. 2014 Sep 12;276:48–71. pmid:24378955
- 53. Tau GZ, Peterson BS. Normal Development of Brain Circuits. Neuropsychopharmacology. 2010 Jan;35(1):147–68. pmid:19794405
- 54. Bjartmar C, Battistuta J, Terada N, Dupree E, Trapp BD. N-acetylaspartate is an axon-specific marker of mature white matter in vivo: a biochemical and immunohistochemical study on the rat optic nerve. Ann Neurol. 2002 Jan;51(1):51–8. pmid:11782984
- 55. Alexander B, Kelly CE, Adamson C, Beare R, Zannino D, Chen J, et al. Changes in neonatal regional brain volume associated with preterm birth and perinatal factors. NeuroImage. 2019 15;185:654–63. pmid:30016676
- 56. Deoni SCL, Dean DC, O’Muircheartaigh J, Dirks H, Jerskey BA. Investigating white matter development in infancy and early childhood using myelin water faction and relaxation time mapping. Neuroimage. 2012 Nov 15;63(3):1038–53. pmid:22884937
- 57. Kelly CE, Cheong JLY, Gabra Fam L, Leemans A, Seal ML, Doyle LW, et al. Moderate and late preterm infants exhibit widespread brain white matter microstructure alterations at term-equivalent age relative to term-born controls. Brain Imaging Behav. 2016 Mar;10(1):41–9. pmid:25739350
- 58. DiPietro JA, Voegtline KM. The gestational foundation of sex differences in development and vulnerability. Neuroscience. 2017 Feb 7;342:4–20. pmid:26232714
- 59. Wheelock MD, Hect JL, Hernandez-Andrade E, Hassan SS, Romero R, Eggebrecht AT, et al. Sex differences in functional connectivity during fetal brain development. Dev Cogn Neurosci. 2019 Mar 5;36:100632. pmid:30901622
- 60. McCarthy MM, Arnold AP, Ball GF, Blaustein JD, De Vries GJ. Sex Differences in the Brain: The Not So Inconvenient Truth. J Neurosci. 2012 Feb 15;32(7):2241–7. pmid:22396398
- 61. McCarthy MM. Is sexual differentiation of brain and behavior epigenetic? Curr Opin Behav Sci. 2019 Feb 1;25:83–8. pmid:31106255
- 62. Else-Quest NM, Hyde JS, Goldsmith HH, Van Hulle CA. Gender differences in temperament: a meta-analysis. Psychol Bull. 2006 Jan;132(1):33–72. pmid:16435957
- 63. Eriksson M, Marschik PB, Tulviste T, Almgren M, Pérez Pereira M, Wehberg S, et al. Differences between girls and boys in emerging language skills: evidence from 10 language communities. Br J Dev Psychol. 2012 Jun;30(Pt 2):326–43. pmid:22550951
- 64. Galsworthy MJ, Dionne G, Dale PS, Plomin R. Sex differences in early verbal and non-verbal cognitive development. Dev Sci. 2000 May 1;3(2):206–15.
- 65. Özçalişkan Ş, Goldin-Meadow S. Sex differences in language first appear in gesture. Dev Sci. 2010 Sep 1;13(5):752–60. pmid:20712741
- 66. Van Kooij BJM, Benders MJNL, Anbeek P, Van Haastert IC, De Vries LS, Groenendaal F. Cerebellar volume and proton magnetic resonance spectroscopy at term, and neurodevelopment at 2 years of age in preterm infants. Dev Med Child Neurol. 2012 Mar;54(3):260–6. pmid:22211363
- 67. Gadin E, Lobo M, Paul DA, Sem K, Steiner KV, Mackley A, et al. Volumetric MRI and MRS and early motor development of infants born preterm. Pediatr Phys Ther. 2012;24(1):38–44. pmid:22207464
- 68. Shaw P, Greenstein D, Lerch J, Clasen L, Lenroot R, Gogtay N, et al. Intellectual ability and cortical development in children and adolescents. Nature. 2006 Mar 30;440(7084):676–9. pmid:16572172
- 69. Kreis R, Hofmann L, Kuhlmann B, Boesch C, Bossi E, Hüppi PS. Brain metabolite composition during early human brain development as measured by quantitative in vivo 1H magnetic resonance spectroscopy. Magn Reson Med. 2002 Dec;48(6):949–58. pmid:12465103
- 70. Tanifuji S, Akasaka M, Kamei A, Araya N, Asami M, Matsumoto A, et al. Temporal brain metabolite changes in preterm infants with normal development. Brain Dev. 2017 Mar;39(3):196–202. pmid:27838187
- 71. Kraemer HC, Yesavage JA, Taylor JL, Kupfer D. How can we learn about developmental processes from cross-sectional studies, or can we? Am J Psychiatry. 2000 Feb;157(2):163–71. pmid:10671382