- 1Department of Women's and Children's Health, Karolinska Institute, Stockholm, Sweden
- 2Clinical Epidemiology Division, Department of Medicine, Karolinska Institute, Stockholm, Sweden
- 3Department of Psychology, Faculty of Health Sciences, University of Deusto, Bilbao, Bizkaia, Spain
Objective: Structural brain volumetric differences have been investigated previously in very preterm children. However, children born extremely preterm, at the border of viability, have been studied to a lesser degree. Our group previously analyzed children born extremely preterm at term using voxel-based morphometry. In this study, we aimed to examine regional gray and white matter differences for children born extremely preterm derived from the same cohort during childhood. We also aimed to explore the effect of perinatal risk factors on brain volumes in the same group.
Methods: At 10 years of age, 51 children born extremely preterm (before 27 weeks and 0 days) and 38 term-born controls with high-quality 3.0 Tesla magnetic resonance images were included. Statistical analyses using voxel-based morphometry were conducted on images that were normalized using age-specific templates, modulated, and smoothed. Analyses were also performed in stratified groups of children born extremely preterm in the absence or presence of perinatal risk factors that have previously been shown to be associated with volumetric differences at term.
Results: We found volumetric decreases in gray and white matter in the temporal lobes, gray matter decreases in the precuneus gyri, and white matter decreases in the anterior cingulum for children born extremely preterm (all p < 0.001, and pfwe< 0.05). Gray and white matter increases were predominantly observed in the right posterior cingulum and occipital lobe (all p < 0.001, and pfwe< 0.05). Of the examined perinatal risk factors, intraventricular hemorrhage grades I-II compared with no intraventricular hemorrhage and patent ductus arteriosus ligation compared with no treated patent ductus arteriosus or patent ductus arteriosus treated with ibuprofen led to volumetric differences at 10 years of age (all p < 0.001, and pfwe< 0.05).
Conclusions: Children born extremely preterm exhibit volumetric alterations in a pattern overlapping that previously found at term, where many regions with differences are the main hubs of higher order networks. Some, but not all, risk factors known to be associated with structural alterations at term were associated with alterations at 10 years of age.
1. Introduction
In Sweden, 3.3 of 1,000 children are born before 27 weeks and 0 days and thus defined as extremely preterm (EPT) (1). Survival rates are high and have improved during the last decade (1). However, these children still have a high risk for adverse neurological outcomes. At the age of 6.5 years, 33% of children born EPT in Sweden during 2004–2007 had moderate-to-severe cognitive disability (2). Magnetic resonance imaging (MRI) has enabled us to further investigate the brain morphology of children born preterm, but the mechanisms underlying developmental challenges are still not fully understood.
Children born very preterm (before 32 weeks and 0 days) and EPT have reduced brain volumes at a global level during childhood compared with term-born children, and these volumetric reductions have been associated with adverse cognitive outcomes (3, 4). We have previously reported that global gray matter (GM) volume was reduced more than global white matter (WM) volume at term for children born EPT and that the relationship was inverted at 10 years of age (4, 5).
Voxel-based morphometry (VBM) encompasses a voxelwise comparison of brain tissue volumes between patient groups and provides information about where volumetric differences are localized (6). We have previously applied this method at term for children born EPT compared with term-born controls (5). The children born EPT exhibited a pattern of widespread regional volumetric alterations (5). Reduced GM and WM brain volumes were found in the temporal lobes and subcortical structures, whereas enlarged brain volumes were found in regions involved in visual processing (5).
To our knowledge, there are no previous studies investigating only children born EPT with VBM beyond infancy, but very preterm adolescents and young adults have been examined (7–9). Reduced volumes and growth of the temporal lobes are consistent findings for children born very preterm compared with term-born controls both during childhood, adolescence, and adulthood (7, 9). Furthermore, various volumetric alterations within the higher order brain networks are common findings that defy straightforward summary (7, 8). The increased volumes in occipital regions reported in infancy for children born EPT have not been found in assessments of older children born very preterm (8).
It remains to be elucidated whether these patterns of GM and WM alterations found at term are distributed similarly at an older age for children born EPT and whether volumetric differences for children born EPT in childhood are different from those previously observed in children born very preterm.
We have previously examined the effect of perinatal risk factors on brain volumes in children born EPT at term and found that those with intraventricular hemorrhage (IVH) grades I-II, patent ductus arteriosus (PDA) ligation, PDA treated with ibuprofen, and the lowest gestational age (GA) were associated with volumetric differences (5).
Low-grade IVH has been considered to not affect neurodevelopment, but at the time of writing, the body of literature suggests that even low-grade IVH can result in adverse neurodevelopmental outcomes, but to a lesser degree than high-grade IVH (10, 11). Whether PDA ligation affects neurodevelopment or not remains unclear. Studies have indicated adverse neurodevelopmental outcomes for children who underwent surgery for PDA (12), and others show no effect of PDA ligation on neurodevelopment (13). Other groups that have investigated the brain volumes of very preterm children have reported associations between brain volume, PDA ligation, and high-grade IVH (14, 15). There is a lack of studies that focus on the most immature children born EPT in relation to brain volumes and perinatal risk factors beyond infancy.
The main aim of this study was to examine the regional brain volume differences in a cohort of children born EPT compared with term-born controls at 10 years of age using VBM.
The secondary aim was to investigate the effect of perinatal risk factors on regional brain volume at 10 years of age. The risk factors examined were IVH grades I-II, PDA (treated with ibuprofen or ligation), and grade of prematurity.
2. Materials and methods
2.1. Participants
In this population-based study, 128 children born EPT (defined as birth before GA 27 weeks and 0 days) in Stockholm between 1 January 2004 and 31 March 2007, who survived to term (GA 40 weeks and 0 days) were invited for an MRI examination at ~10 years of age. The cohort is presented in Figure 1 and includes the same children born EPT that were previously described to investigate global brain volumes (4). In total, four children had moved away from the Stockholm area before term MRI, and 14 children were lost to follow-up. Children with cystic periventricular leukomalacia, IVH grade III, periventricular hemorrhagic infarction (previously labeled IVH grade IV), severe white matter abnormalities as defined by a published scoring system (16), focal brain lesions, or severe medical conditions were excluded. In all, six children born EPT were excluded because of severe medical conditions, detailed as follows: one due to Trisomi 21, one had acquired brain lesions, one had hemophagocytic lymphohistiocytosis, and three had congenital malformations (myelomeningocele, esophageal atresia, and vascular malformation). Seven children were excluded due to an interrupted MRI examination, inaccurate segmentation, or low-quality MRI (defined as images with motion artifacts or blurry contrast between GM and WM), leaving 51 children born EPT with high-quality MRI images. In total, nine of the included children born EPT were from multiple pregnancies, where three mothers had two children each included from a multiple pregnancy. Drop-out analyses were performed to compare the 51 children born EPT included with the 41 children excluded due to low quality on MRI or declined participation.
Figure 1. Flow chart of children born extremely preterm and term-born controls. MRI, magnetic resonance imaging.
Healthy, singleton, term-born children were randomly selected from the Swedish Medical Birth Registry and matched to children born EPT for place of birth, sex, day of birth, and maternal country of birth at age 2.5 years as part of the National Extremely Preterm Infants in Sweden Study (EXPRESS) (17). The Stockholm-based, full-term children were invited for an MRI examination at ~10 years of age, and of the 77 children who received an invite, 38 had high-quality MRI data and were included in this study. For a summary of the study group, see Figure 1.
This study was approved by the Stockholm Ethics Review Board, and written informed consent was obtained from the parents.
2.2. Perinatal risk factors and clinical characteristics
Perinatal information and clinical characteristics were gathered from medical records. Ultrasound scans were performed by a trained neonatologist to detect IVH during the first 3 days after birth, once a week until 27 weeks of gestation, and then every other week until term. Grading was performed according to the Papile classification (18). The definition of PDA in this study was the need for either PDA ligation or treatment with ibuprofen. The treatment decision was at the discretion of the neonatologist in charge (19). In general, children born EPT with a hemodynamically significant PDA according to the preceding echocardiographic evaluation were primarily treated with ibuprofen. The PDA ligation was then indicated if pharmacological therapy failed or when contraindications to ibuprofen were present. The definition of small for gestational age (SGA) was a birthweight < 2 SD below the mean. Bronchopulmonary dysplasia was defined as the need for oxygen at 36 weeks' gestation, and necrotizing enterocolitis was defined using Bell criteria (20).
The effects of perinatal risk factors were investigated by creating subgroups based on the presence or absence of the examined risk factors. The subgroups were as follows: IVH grades I-II or no IVH, PDA ligation or no treated PDA, PDA ligation or PDA treated with ibuprofen, PDA treated with ibuprofen or no treated PDA, and birth ≤GA 25 + 6 weeks or birth ≥GA 26 + 0 weeks. Risk factors were selected because they have been shown to be associated with altered brain volumes at term in children from the same cohort (5).
2.3. MRI acquisition at 10 years of age, preprocessing, and segmentation
Images were acquired using a SIGMA HDx General Electric 3.0-T MRI system (GE Healthcare, Milwaukee, WI, USA), and MRI examinations were performed without sedation. The protocol included a sagittal 3D T1-weighted image with the BRAVO SPGR sequence with a time to inversion of 400 ms, a field of view of 240 × 240 mm2, a flip angle of 12°, a slice thickness of 1.0 mm, and a voxel size of 1 × 0.938 × 0.938 mm3. All images were assessed by a neuroradiologist.
The T1-weighted images were reoriented in the plane of the anterior commissures, and the skull was stripped using the Brain Extraction Tool from the FMRIB Software Libraries version 5.0.1 (FMRIB Laboratory, University of Oxford, England, UK) (21). Remaining non-brain tissue after automatic brain extraction was removed through manually editing. The criteria to perform manual editing were based on anatomical references and guided by an experienced brain researcher (NP).
The pre-processed images were then automatically segmented into GM, WM, and cerebrospinal fluid (CSF) using the default settings in unified segmentation (22), carried out on the SPM8 software (Wellcome Department, University College, London, UK) using MATLAB version 7.5 (MathWorks, Natick, MA, USA). Age-specific templates from the Template-O-Matic (TOM) toolbox were used (23). The images were spatially non-linearly normalized using DARTEL (24). The segmented images were finally modulated with Jacobian modulation to preserve the amount of volume and subsequently smoothed using a 6 mm Gaussian kernel at full width at half maximum. The aims of smoothing were to improve the signal-to-noise ratio, ensure that the imaging data satisfy a Gaussian distribution, and overcome alignment problems.
Images were inspected before and after each processing step to assess quality. High-quality images were defined as images without motion artifacts, or blurry contrast between GM and WM. Out of the children born EPT with complete MRI data, one child was excluded due to motion artifacts, one had poor contrast between GM and WM, two had braces that generated artifacts, and one image was excluded after unified segmentation since it did not accurately segment brain tissues (Figure 1). Out of the term-born controls with complete MRI data, four children were excluded due to poor contrast between GM and WM and three children were excluded after inspection following unified segmentation due to incorrect classification of brain tissues (Figure 1).
The Easy Volume toolbox was used to calculate the intracranial volume (ICV), total WM, and total GM, and the method has been previously described (4).
2.4. Voxel-based morphometry analysis
Segmented, warped, smoothed, and modulated tissues were used for statistical analyses using a two-sample t-test design, adjusting for covariates within the VBM toolbox, using the SPM12 software (Wellcome Department, University College, London, UK).
In the first VBM analyses, children born EPT were compared with term-born controls. Both reduced and increased brain volumes for children born EPT were analyzed by inverting the contrast. The analyses were adjusted for sex, age at scan, and ICV. As analyses were adjusted for ICV, they are interpreted as relative volume differences.
Analyses were repeated with only singleton children, and the results are shown in Supplementary material.
Thereafter, the brain volumes of children born EPT with and without perinatal risk factors were compared using VBM. Both reduced and increased brain volumes for children born EPT with the examined perinatal risk factor were analyzed by inverting the contrast. In these analyses, GA, ICV, age at scan, and sex were used as covariates, except when children born EPT ≤GA 25 + 5 weeks or ≥GA 26 + 0 weeks were compared where only sex, age at scan, and ICV were adjusted for. When children born EPT treated with PDA ligation were compared with children with no treated PDA or with PDA treated with ibuprofen, we added IVH grades I-II as a covariate. When IVH grades I-II were compared with no IVH, PDA ligation was added as a covariate. However, since these adjustments did not change the results, these covariates were not included in the final analyses.
The initial cluster-forming height threshold was set at p < 0.001. All results presented in tables and figures were then family-wise error (FWE) corrected (p = 0.05) at cluster level using the Gaussian random field theory to account for multiple comparisons. Absolute threshold masking was set at 0.1 to avoid possible edge effects between tissue types (25). We report the corrected cluster-level p-values and their corresponding Montreal Neurological Institute (MNI) peak coordinates. The Harvard-Oxford cortical and subcortical structural atlases (26) and a conversion of the original Talairach Atlas registered to MNI space were used for anatomical orientation (27).
2.5. Statistical analysis
All the statistical analyses of clinical characteristics and total brain volumes were performed using SPSS, version 25 (IBM Corp, Armonk, NY, USA). Variables were tested for normality and homogeneity prior to analysis. Normally distributed continuous variables were compared using the two-sample t-test, and non-normally distributed variables were compared using the Mann–Whitney U-test. Categorical variables were compared using Pearson's chi-squared or Fisher's exact test, as appropriate. The ICV, total GM, and total WM volumes were compared between children born EPT with and without perinatal risk factors. These total brain volumes were first univariately analyzed and then analyzed as dependent variables in a multivariate general linear model with group status (presence or absence of a perinatal risk factor) as independent variables, adjusting for sex, age at scan, and GA. A statistical significance level of p < 0.05 was used.
3. Results
The clinical characteristics of the children born EPT and term-born controls are summarized in Table 1. Drop-out analyses of the children born EPT included in the study compared with the children declining participation or excluded due to low MRI quality are presented in Supplementary Table 1, and the included children born EPT had a higher GA and a lower rate of bronchopulmonary dysplasia than the children not included.
In Supplementary Tables 2a–e, the clinical characteristics and total brain volumes of the children born EPT with and without perinatal risk factors are presented. Supplementary Table 2a describes the children born EPT with IVH grades I-II, Supplementary Tables 2b–d describes the groups treated for PDA, and Supplementary Table 2e compares children born EPT born ≤GA 25 weeks and 6 days with those born ≥GA 26 weeks and 0 days. Of the 51 included children born EPT, there were 16 treated with PDA ligation, and 14 out of these were previously treated with ibuprofen. There were 20 children treated with only ibuprofen. Of the analyzed perinatal risk factors, only the children born EPT with PDA ligation had reduced ICV compared with children born EPT with no treated PDA (mean difference 81.8 cm3 and 95% confidence interval 2.2–161.4, p = 0.044) and children born EPT with PDA treated with ibuprofen (mean difference 80.0 cm3 and 95% confidence interval 3.3–156.6, p = 0.041), but this did not remain significant after adjusting for sex, age at scan, and GA, see Supplementary Tables 2b, c. No significant differences in ICV, total GM, or total WM brain volumes were found for children born with IVH grades I-II compared with children with no IVH (Supplementary Table 2a), PDA treated with ibuprofen compared with no treated PDA (Supplementary Table 2d), or ≤GA 25 weeks and 6 days compared with ≥GA 26 weeks and 0 days (Supplementary Table 2e). IVH grades I-II were more common in children treated with PDA ligation, and PDA ligation was more common in children with IVH grades I-II (Supplementary Tables 2a–c).
3.1. Regional volumetric differences for children born EPT compared with term-born controls
The children born EPT exhibited significant regional volumetric differences with increases and decreases in both GM and WM when adjusted for sex, age at scan, and ICV (Table 2, Figure 2).
Table 2. Regional volumetric differences between children born EPT (n = 51) and term-born controls (n = 38) analyzed with voxel-based morphometry at 10 years of age, adjusted for intracranial volume, sex, and age at scan.
Figure 2. Sagittal, coronal, and axial views of regional volumetric differences for children born extremely preterm compared with term-born controls. (A) Gray matter decreased, (B) white matter decreased, (C) gray matter increased, (D) white matter increased. The colored bars represent t-score. Differences are mapped on the render of one child born extremely preterm. Analyses are adjusted for sex, age at scan, and intracranial volume. Display orientation right = right. Results are presented with a threshold of p < 0.001, with family-wise error correction p < 0.05 at the cluster level.
The most prominent GM reductions were found in the temporal lobes. The large clusters had their peek MNI coordinates in the superior temporal gyri and were distributed bilaterally. The volumetric reduction extended medially to the insula and inferiorly to the middle temporal gyri. However, there were also significant GM reductions in children born EPT in the precuneus cortex bilaterally within the parietal lobe.
The WM reductions for the children born EPT were localized adjacent to the GM reductions and were most prominent in the middle temporal gyri. The large WM cluster with peak coordinates in the right middle temporal gyrus extended into the brainstem and also included the right thalamus, while the left WM cluster was only localized in the temporal lobe. In addition, reduced WM volumes were found in the anterior cingulum and cerebellum.
The most prominent volumetric increases for children born EPT compared with term-born controls for GM were localized in the cingulate gyrus with its peak in the right posterior cingulate cortex. Furthermore, there were GM increases in the lateral occipital lobe with bilateral distribution.
The WM increases were adjacent to the GM increases and involved the cingulum as well as occipital regions. All reported results were significant at p < 0.001 and cluster-level pfwe < 0.05.
Analyses were repeated with only the n=42 singleton children born EPT and added as Supplementary Table 3. As the regions with differences did not change, we chose to keep all children born EPT in the final analyses.
3.2. Regional volumetric differences for children born EPT with and without perinatal risk factors
All significant regional volumetric differences from analyses using VBM for children born EPT with and without perinatal risk factors are summarized in Table 3 and Figure 3. All reported results were significant at p < 0.001, and cluster-level pfwe < 0.05. Of the investigated perinatal risk factors, PDA ligation and IVH grades I-II resulted in volumetric differences.
Table 3. Volumetric differences for children born EPT with and without perinatal risk factors, results were adjusted for sex, gestational age, intracranial volume, and age at scan.
Figure 3. Regional volumetric differences for children born EPT with and without perinatal risk factors. (A) Gray matter decreased for IVH grades I-II compared with no IVH, (B) white matter decreased for PDA ligation compared with no treated PDA, (C) gray matter increased for PDA ligation compared with no treated PDA, (D) gray matter increased for PDA ligation compared with PDA treated with ibuprofen, (E) white matter increased for PDA ligation compared with PDA treated with ibuprofen. The color bar represents t-score. Differences are displayed on the render of one child born EPT. Analyses are adjusted for gestational age, sex, age at scan, and intracranial volume. Display orientation: right = right. EPT, extremely preterm; IVH, intraventricular hemorrhage; PDA, patent ductus arteriosus. Results are presented with a threshold of p < 0.001, with family-wise error correction p < 0.05 at the cluster level.
3.2.1. Children born EPT with IVH grades I–II compared with children born EPT with no IVH
Children born EPT with IVH grades I-II had reduced GM volume in the left hippocampus compared with children born EPT without IVH when adjusted for sex, age at scan, GA, and ICV (Table 3 and Figure 3). Children with IVH grades I-II did not have any enlarged regions compared with children without IVH.
3.2.2. Children born EPT treated with PDA ligation compared with children born EPT with no treated PDA
Children treated with PDA ligation had a reduced WM volume in the frontal lobe within the premotor area compared with children with no treated PDA (Table 3 and Figure 3). Multiple enlarged GM regions were observed in the occipital lobe (Table 3 and Figure 3). Adjustments for sex, age at scan, GA, and ICV were made.
3.2.3. Children born EPT treated with PDA ligation compared with children born EPT treated with ibuprofen
GM increases were observed in the occipital lobe when children treated with PDA ligation were compared with children treated with ibuprofen (Table 3 and Figure 3), when adjusted for sex, age at scan, GA, and ICV.
3.2.4. Children born EPT ≤ GA 25 weeks and 6 days compared with ≥GA 26 weeks and 0 days, and children born EPT treated with ibuprofen compared with no treated PDA
Children born EPT ≤ GA 25 weeks and 6 days compared with ≥GA 26 weeks and 0 days did not exhibit any regional volumetric differences after adjustments for sex, age at scan, and ICV. Likewise, treatment with ibuprofen compared with no treated PDA did not result in any regional volumetric differences when adjusted for sex, age at scan, GA, and ICV.
4. Discussion
This study demonstrated widespread differences in regional GM and WM volumes for children born EPT compared with term-born controls at 10 years of age. GM decreases were mainly localized bilaterally in the temporal and precuneus gyri, while the GM increases were prominent in the posterior cingulate gyrus on the right side and occipital lobe. The differences in WM and GM followed similar patterns, for example, when the WM decreases were localized in the temporal lobes and the WM increases were found in the occipital lobe and the right posterior cingulate gyrus. Regions with volumetric differences largely overlapped with those previously found at term (5). Of the examined perinatal risk factors known to affect regional brain volumes at term for children born EPT (5), PDA ligation compared with no treated PDA or PDA treated with ibuprofen and IVH grades I-II compared with no IVH resulted in brain volume differences at 10 years of age.
In typically developing children, GM has an “inverted U-shaped” growth pattern (28). The early increase in GM mirrors the rapid synaptogenesis during the third trimester and first years of life, followed by reduced speed of GM growth during childhood, which is associated with synaptic pruning and the myelination at the border between the cortical GM and WM (29–31). The timing of these processes varies for different regions of the brain (29, 30). The WM volume follows a more linear pattern of growth and increases relative to ICV up to the fourth decade of life (32). The increased WM volume during childhood reflects the myelination and maturation of brain connections (28, 29).
One possible explanation for the volumetric differences between children born EPT and term-born controls found in this study could be associated with the early stressor of abrupt extrauterine exposure, potentially modifying the fine-tuned processes occurring in the third trimester, involving premyelinating oligodendrocytes, synaptogenesis, migration, and subplate neurons, thus leading to deviations from typical brain trajectories (31).
Patterns of reduced brain tissue in some regions, accompanied by increased amounts of brain tissue in other regions, are expected because changes in one region often lead to a cascade of alterations in other regions (33).
As expected, WM differences were found adjacent to the GM differences. WM and GM are all components of the same neural circuit and are highly interdependent because there are inseparable connections between neurons, glial cells, and myelin (29).
We found more prominent volumetric increases for GM than for WM volumes when children born EPT were compared with term-born controls. Our group has previously shown that the total GM volume relative to ICV is increased at 10 years of age in the same cohort, while the relative total WM volume is decreased. The findings in this study provide a clearer image of the locations of these volumetric differences (4).
4.1. Regional volumetric differences for children born EPT compared with term-born controls
This study suggests selective vulnerability in brain regions, possibly related to the developmental stage in different areas of the brain, as previously suggested (7).
There were bilateral GM volume reductions in the temporal lobes for children born EPT compared with term-born controls. This pattern of volumetric changes matches the alterations previously demonstrated at term (5). Reduced GM volumes in the temporal gyri are supported by previous research investigating older children born very preterm (7, 8, 34), and the temporal lobes have also been shown to have reduced cortical thickness in children born very preterm (35). During the second half of pregnancy, there is a rapid increase in the synaptic density in the temporal lobe, which could be a vulnerable area for exposure to the extrauterine environment (36, 37). There is rapid in utero growth of subplate thickness in the superior temporal gyri during the third trimester, and in this study, we observed that this region was particularly affected (38). The temporal gyri are involved in multiple cognitive processes and are thus central to emotional regulation (7). It is understood that children born preterm are at risk for both cognitive and emotional deficits (2, 39). The superior temporal gyri are also critical components of language processing (40–42). Even in the absence of major disabilities, preterm children are more prone to experiencing difficulties in language functions than term-born children (42). Volumetric differences are associated with neurodevelopmental outcomes, such as cognition and language (3).
Interestingly, many of the regions where we found volumetric differences were the main hubs within higher order networks, which are important for sensory integration and socio-emotional regulation, and these hubs are present before birth (43). We know that children born EPT have a high rate of behavioral problems (44, 45). It has been proposed that the main cortical hubs are vulnerable to early adverse events, such as extremely preterm birth (43, 45).
We found reduced GM volume in the precuneus cortex bilaterally in children born EPT and enlarged GM and WM volumes in the posterior cingulate gyrus, which are important hubs within the default mode network (46). We also found reduced WM in the anterior cingulum, an important hub within the salience network, governing the ability to segregate external and internal stimuli and switch the activation between neural networks to guide behavior and attention (47). The observed volumetric alterations in regions within higher order networks are consistent with previous literature on children born very preterm (7, 8). Reduced WM in the anterior cingulum at term was demonstrated previously, indicating a stable alteration throughout childhood (5). We also found WM reductions in the brainstem, which have been observed previously in children derived from this cohort at term and in very preterm children at older assessments (5, 7).
Multiple instances of GM increase were found in the visual areas of the occipital lobe in children born EPT compared with term-born controls, and this was also found at term (5, 48). The enlarged volumes could be due to defect pruning programs or compensatory mechanisms with early onset (33). However, these findings at term may suggest that extrauterine visual exposure in children born EPT has led to increased volumes (5). The underlying mechanisms for these relative volumetric increases and how they relate to visual function are an issue for future investigation. It is important to note that as all the analyses are adjusted for ICV, the volumetric increases reflect the relative size of the specific region in the brain.
Previous research has suggested links between early cerebellar growth and cortical development at older ages in preterm populations (49, 50). Future studies could further investigate the association of early cerebellar growth and the final volumes, with special attention drawn to the regions affected in this study.
4.2. Brain volumes for children born EPT with and without perinatal risk factors
Of the examined perinatal risk factors, IVH grades I-II and PDA ligation were associated with brain volume differences.
For children born EPT with IVH grades I-II, the hippocampal volume was reduced on the left side when compared to children born EPT without IVH when analyses were adjusted for sex, GA, age at scan, and ICV. The results remained significant when adjusted for PDA ligation. Recent findings have shown that even low-grade IVH is associated with adverse cognitive outcomes compared with children with no IVH (10, 11). The ICV, total GM, and total WM did not differ for children born EPT with IVH grades I-II when compared to children born with no IVH; thus, the volumetric alterations were limited. Reduced hippocampal volume in less immature children born very preterm with any grade of IVH compared with preterm children without IVH has been previously reported (15). The anatomical proximity of the hippocampus to the ventricles could lead to a particular vulnerability to perinatal events that cause hypoxia and ischemia, such as IVH (51). Reduced hippocampal volumes have been associated with memory performance (15, 51, 52). Reduced hippocampal volume was only significant on the left side of the brain in this study, and preceding research using the same methodology has shown that reduced hippocampal volumes on the left side are associated with memory outcome (52).
At term, IVH grades I-II were associated with both global and regional volumetric differences (5). However, at the regional level, the areas of reduced GM volume were not consistent at term and at 10 years of age, which could indicate adaptational growth and altered brain trajectories (5).
Children born EPT treated with PDA ligation had reduced ICV compared with children with no treated PDA or PDA treated with ibuprofen, but after adjustments for GA, age at scan, and sex, the differences were not significant. When examined with VBM, the children born EPT treated with PDA ligation had reduced WM volume in the premotor cortex within the frontal lobe after adjustments for GA, sex, and ICV, and the results remained significant after adjustments for IVH grades I–II. Surgical closure of the PDA has been associated with neurosensory impairments, including a higher risk of cerebral palsy (12). Importantly, it is not possible to distinguish the risk of the surgery from the risk of hypooxygenation to the brain from a long-lasting hemodynamically significant ductus (14). Enlarged brain volumes were found in the occipital lobe of children treated with PDA ligation compared with children born EPT with no treated PDA or PDA treated with ibuprofen. These volumetric differences could be a sign of reorganization; however, this needs to be elucidated further in future research.
At term age, the children born EPT treated with PDA ligation had reduced brain volumes in multiple brain regions for both GM and WM, including the temporal gyri and cerebellum, when compared to children with no treated PDA and children treated with ibuprofen for PDA, as previously reported (5). Thus, this risk factor was associated with more prominent alterations at term than at 10 years of age, and there may be favorable compensatory growth during childhood.
No volumetric differences were found for children born EPT treated with ibuprofen for PDA compared with children with no treated PDA, while at term, the children born EPT treated with ibuprofen for PDA had enlarged volumes (5).
Interestingly, no global or regional volumetric differences were found for children born ≤GA 25 weeks and 6 days compared with children born ≥GA 26 weeks and 0 days. The effect of GA on brain volume is known from the literature (7, 53). However, IVH grades I-II and PDA ligation had larger impacts on regional brain volumes than prematurity per se in our homogeneous sample of children born EPT without major brain lesions at 10 years of age.
4.3. Strengths and limitations
A strength of this study is its population-based study design and its focus on children specifically born EPT before 27 gestational weeks without major brain lesions. Furthermore, this study extends beyond early childhood, assesses the regional brain volume growth, and evaluates the effect of perinatal risk factors on regional brain volumes.
This study had several limitations. First, because the study population was limited, the stratified groups were small when examining perinatal risk factors. For the same reason, we could not perform separate analyses on children with IVH grade I or IVH grade II only. Likewise, we could not perform sub-analyses of children with unilateral or bilateral IVH or if the IVH was on the right or left side. We could also not stratify children treated with only PDA ligation and children treated with PDA ligation and ibuprofen as the numbers were too few.
Furthermore, PDA treatment protocols and the operating technique for ligation vary between centers and countries, which may limit the generality of the results.
Also, the stratification for GA was made unevenly with children born from 22 weeks and 0 days up to 25 weeks and 6 days in one group, and only children from 26 weeks and 0 days up to 26 weeks and 6 days in the other group. This choice was made as there were too few surviving children with the lowest gestational ages, and to homogenize the stratification with previous analyses done at term.
We are aware that software and pipelines for image processing are constructed for adult brains. We did however use an age-specific atlas to increase the robustness of data averaging. Also, the VBM methodology has previously been used in multiple studies in similar populations (8).
Some of the volumetric differences were large and expanded over more than one brain region, which hampers the anatomical precision. In future, research could further investigate these regions with other methodologies, such as TBSS, or structural covariance.
Another limitation was the inclusion of children born EPT from multiple pregnancies. However, the VBM analyses were repeated with only singleton children and the results remained similar and were added as Supplementary material.
5. Conclusion
Children born EPT had altered regional brain volumes at 10 years of age compared with term-born controls. Surgically ligated PDA and IVH grades I-II were associated with alterations in regional brain volume growth. The results emphasize the need to optimize brain growth during the neonatal period and beyond and to limit the exposure to preventable perinatal risk factors. Longitudinal studies, including cognitive performance and language data, would help to determine the significance of these volumetric brain alterations.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by the Regional Ethical Review Board in Stockholm. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.
Author contributions
HK, JB, NP, LB, and UÅ: conceptualization. HK, LB, LF, NP, and UÅ: methodology. HK, NP, and LF: formal analysis. HK: writing—original draft. HK, JB, LB, LF, NP, and UÅ: writing—review and editing. HK and NP: visualization. UÅ and NP: supervision. UÅ: funding acquisition, data curation, and resources. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by grants from the Swedish Medical Research Council (grant number, 2017-03043) (UÅ), the regional agreement on medical training and clinical research (grant number, ALF SLL 20170243) between Stockholm County Council and the Karolinska Institutet (UÅ), the Swedish Order of Freemasons in Stockholm (NP), the Swedish Medical Society (UÅ), the Swedish Brain Foundation (grant number, FO2019-0045) (UÅ), the Philipson Foundation (UÅ), the Crown Princess Lovisa's Foundation (NP), the Sällskapet Barnavård, and the Erik and Edith Fernström Stiftelse (NP); Region Stockholm (clinical postdoctoral appointment) (JB); Region Stockholm (combined clinical residency and Ph.D. training program) (HK). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Acknowledgments
We are grateful to the research nurse, Lena Swartling-Schlinzig, who organized the MRI examinations in the cohort at 10 years of age. We are also grateful to all participating children and parents.
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/fneur.2023.1148781/full#supplementary-material
Abbreviations
EPT, extremely preterm; VBM, voxel-based morphometry; IVH, intraventricular hemorrhage; MRI, magnetic resonance imaging; PDA, patent ductus arteriosus; GM, gray matter; WM, white matter; GA, gestational age; ICV, intracranial volume.
References
1. Norman M, Hallberg B, Abrahamsson T, Bjorklund LJ, Domellof M, Farooqi A, et al. Association between year of birth and 1-year survival among extremely preterm infants in sweden during 2004–2007 and 2014–2016. JAMA. (2019) 321:1188–99. doi: 10.1001/jama.2019.2021
2. Serenius F, Ewald U, Farooqi A, Fellman V, Hafström M, Hellgren K, et al. Neurodevelopmental outcomes among extremely preterm infants 6.5 years after active perinatal care in Sweden. JAMA Pediatr. (2016) 170:954–63. doi: 10.1001/jamapediatrics.2016.1210
3. Monson BB, Anderson PJ, Matthews LG, Neil JJ, Kapur K, Cheong JL, et al. Examination of the pattern of growth of cerebral tissue volumes from hospital discharge to early childhood in very preterm infants. JAMA Pediatr. (2016) 170:772–9. doi: 10.1001/jamapediatrics.2016.0781
4. Kvanta H, Bolk J, Strindberg M, Jiménez-Espinoza C, Broström L, Padilla N, et al. Exploring the distribution of grey and white matter brain volumes in extremely preterm children, using magnetic resonance imaging at term age and at 10 years of age. PLoS ONE. (2021) 16:e0259717. doi: 10.1371/journal.pone.0259717
5. Padilla N, Alexandrou G, Blennow M, Lagercrantz H, Aden U. Brain growth gains and losses in extremely preterm infants at term. Cereb Cortex. (2015) 25:1897–905. doi: 10.1093/cercor/bht431
6. Ashburner J, Friston KJ. Voxel-based morphometry–the methods. Neuroimage. (2000) 11(6 Pt 1):805–21. doi: 10.1006/nimg.2000.0582
7. Nosarti C, Giouroukou E, Healy E, Rifkin L, Walshe M, Reichenberg A, et al. Grey and white matter distribution in very preterm adolescents mediates neurodevelopmental outcome. Brain. (2008) 131(Pt 1):205–17. doi: 10.1093/brain/awm282
8. Zhou L, Zhao Y, Liu X, Kuang W, Zhu H, Dai J, et al. Brain gray and white matter abnormalities in preterm-born adolescents: a meta-analysis of voxel-based morphometry studies. PLoS ONE. (2018) 13:e0203498. doi: 10.1371/journal.pone.0203498
9. Karolis VR, Froudist-Walsh S, Kroll J, Brittain PJ, Tseng CJ, Nam KW, et al. Volumetric grey matter alterations in adolescents and adults born very preterm suggest accelerated brain maturation. Neuroimage. (2017) 163:379–89. doi: 10.1016/j.neuroimage.2017.09.039
10. Argyropoulou MI, Astrakas LG, Xydis VG, Drougia A, Mouka V, Goel I, et al. Is Low-grade intraventricular hemorrhage in very preterm infants an innocent condition? structural and functional evaluation of the brain reveals regional neurodevelopmental abnormalities. AJNR Am J Neuroradiol. (2020) 41:542–7. doi: 10.3174/ajnr.A6438
11. Mukerji A, Shah V, Shah PS. Periventricular/intraventricular hemorrhage and neurodevelopmental outcomes: a meta-analysis. Pediatrics. (2015) 136:1132–43. doi: 10.1542/peds.2015-0944
12. Kabra NS, Schmidt B, Roberts RS, Doyle LW, Papile L, Fanaroff A. Neurosensory impairment after surgical closure of patent ductus arteriosus in extremely low birth weight infants: results from the trial of indomethacin prophylaxis in preterms. J Pediatr. (2007) 150:229–34.e1. doi: 10.1016/j.jpeds.2006.11.039
13. Youn Y, Seo YM, Yum SK, Sung IK. Patent ductus arteriosus ligation on neurodevelopmental outcomes at corrected 2 years. Ital J Pediatr. (2019) 45:110. doi: 10.1186/s13052-019-0698-7
14. Lemmers PM, Benders MJ, D'Ascenzo R, Zethof J, Alderliesten T, Kersbergen KJ, et al. Patent ductus arteriosus and brain volume. Pediatrics. (2016) 137:e20153090. doi: 10.1542/peds.2015-3090
15. Fernández de Gamarra-Oca L, Zubiaurre-Elorza L, Junqué C, Solana E, Soria-Pastor S, Vázquez É, et al. Reduced hippocampal subfield volumes and memory performance in preterm children with and without germinal matrix-intraventricular hemorrhage. Sci Rep. (2021) 11:2420. doi: 10.1038/s41598-021-81802-7
16. Inder TE, Wells SJ, Mogridge NB, Spencer C, Volpe JJ. Defining the nature of the cerebral abnormalities in the premature infant: a qualitative magnetic resonance imaging study. J Pediatr. (2003) 143:171–9. doi: 10.1067/S0022-3476(03)00357-3
17. Serenius F, Källén K, Blennow M, Ewald U, Fellman V, Holmström G, et al. Neurodevelopmental outcome in extremely preterm infants at 2.5 years after active perinatal care in Sweden. JAMA. (2013) 309:1810–20. doi: 10.1001/jama.2013.3786
18. Papile LA, Burstein J, Burstein R, Koffler H. Incidence and evolution of subependymal and intraventricular hemorrhage: a study of infants with birth weights less than 1,500 Gm. J Pediatr. (1978) 92:529–34. doi: 10.1016/S0022-3476(78)80282-0
19. Gudmundsdottir A, Broström L, Skiöld B, Källén K, Serenius F, Norman M, et al. The type and timing of patent ductus arteriosus treatment was associated with neurodevelopment when extremely preterm infants reached 6.5 years. Acta Paediatr. (2021) 110:510–20. doi: 10.1111/apa.15452
20. Gregory KE, Deforge CE, Natale KM, Phillips M, Van Marter LJ. Necrotizing enterocolitis in the premature infant: neonatal nursing assessment, disease pathogenesis, and clinical presentation. Adv Neonat Care. (2011) 11:155–64. doi: 10.1097/ANC.0b013e31821baaf4
21. Smith SM. Fast robust automated brain extraction. Hum Brain Mapp. (2002) 17:143–55. doi: 10.1002/hbm.10062
22. Ashburner J, Friston KJ. Unified segmentation. Neuroimage. (2005) 26:839–51. doi: 10.1016/j.neuroimage.2005.02.018
23. Wilke M, Holland SK, Altaye M, Gaser C. Template-O-matic: a toolbox for creating customized pediatric templates. Neuroimage. (2008) 41:903–13. doi: 10.1016/j.neuroimage.2008.02.056
24. Ashburner J. A fast diffeomorphic image registration algorithm. Neuroimage. (2007) 38:95–113. doi: 10.1016/j.neuroimage.2007.07.007
25. Ledoux AA, Boyer P, Phillips JL, Labelle A, Smith A, Bohbot VD. Structural hippocampal anomalies in a schizophrenia population correlate with navigation performance on a wayfinding task. Front Behav Neurosci. (2014) 8:88. doi: 10.3389/fnbeh.2014.00088
26. Desikan RS, Ségonne F, Fischl B, Quinn BT, Dickerson BC, Blacker D, et al. An automated labeling system for subdividing the human cerebral cortex on mri scans into gyral based regions of interest. Neuroimage. (2006) 31:968–80. doi: 10.1016/j.neuroimage.2006.01.021
27. Lancaster JL, Tordesillas-Gutiérrez D, Martinez M, Salinas F, Evans A, Zilles K, et al. Bias between mni and talairach coordinates analyzed using the ICBM-152 brain template. Hum Brain Mapp. (2007) 28:1194–205. doi: 10.1002/hbm.20345
28. Morita T, Asada M, Naito E. Contribution of neuroimaging studies to understanding development of human cognitive brain functions. Front Hum Neurosci. (2016) 10:464. doi: 10.3389/fnhum.2016.00464
29. Lenroot RK, Giedd JN. Brain development in children and adolescents: insights from anatomical magnetic resonance imaging. Neurosci Biobehav Rev. (2006) 30:718–29. doi: 10.1016/j.neubiorev.2006.06.001
30. Tamnes CK, Ostby Y, Fjell AM, Westlye LT, Due-Tonnessen P, Walhovd KB. Brain maturation in adolescence and young adulthood: regional age-related changes in cortical thickness and white matter volume and microstructure. Cereb Cortex. (2010) 20:534–48. doi: 10.1093/cercor/bhp118
31. Dimitrova R, Arulkumaran S, Carney O, Chew A, Falconer S, Ciarrusta J, et al. Phenotyping the preterm brain: characterizing individual deviations from normative volumetric development in two large infant cohorts. Cereb Cortex. (2021) 31:3665–77. doi: 10.1093/cercor/bhab039
32. Coupe P, Catheline G, Lanuza E, Manjon JV. Towards a unified analysis of brain maturation and aging across the entire lifespan: a MRI analysis. Hum Brain Mapp. (2017) 38:5501–18. doi: 10.1002/hbm.23743
33. Johnson MH, Jones EJ, Gliga T. Brain adaptation and alternative developmental trajectories. Dev Psychopathol. (2015) 27:425–42. doi: 10.1017/S0954579415000073
34. Thompson DK, Matthews LG, Alexander B, Lee KJ, Kelly CE, Adamson CL, et al. Tracking regional brain growth up to age 13 in children born term and very preterm. Nat Commun. (2020) 11:696. doi: 10.1038/s41467-020-14334-9
35. Lax ID, Duerden EG, Lin SY, Mallar Chakravarty M, Donner EJ, Lerch JP, et al. Neuroanatomical consequences of very preterm birth in middle childhood. Brain Struct Funct. (2013) 218:575–85. doi: 10.1007/s00429-012-0417-2
36. Huttenlocher PR, Dabholkar AS. Regional differences in synaptogenesis in human cerebral cortex. J Compar Neurol. (1997) 387:167–78. doi: 10.1002/(SICI)1096-9861(19971020)387:2<167::AID-CNE1>3.0.CO;2-Z
37. Kesler SR, Vohr B, Schneider KC, Katz KH, Makuch RW, Reiss AL, et al. Increased temporal lobe gyrification in preterm children. Neuropsychologia. (2006) 44:445–53. doi: 10.1016/j.neuropsychologia.2005.05.015
38. Corbett-Detig J, Habas PA, Scott JA, Kim K, Rajagopalan V, McQuillen PS, et al. 3d Global and regional patterns of human fetal subplate growth determined in utero. Brain Struct Funct. (2011) 215:255–63. doi: 10.1007/s00429-010-0286-5
39. Linsell L, Johnson S, Wolke D, Morris J, Kurinczuk JJ, Marlow N. Trajectories of behavior, attention, social and emotional problems from childhood to early adulthood following extremely preterm birth: a prospective cohort study. Eur Child Adolesc Psychiatry. (2019) 28:531–42. doi: 10.1007/s00787-018-1219-8
40. Eklof E, Martensson GE, Aden U, Padilla N. Reduced structural brain asymmetry during neonatal life is potentially related to autism spectrum disorders in children born extremely preterm. Autism Res. (2019) 12:1334–43. doi: 10.1002/aur.2169
41. Friederici AD. The brain basis of language processing: from structure to function. Physiol Rev. (2011) 91:1357–92. doi: 10.1152/physrev.00006.2011
42. Stipdonk LW, Franken MJP, Dudink J. Language outcome related to brain structures in school-aged preterm children: a systematic review. PLoS ONE. (2018) 13:e0196607. doi: 10.1371/journal.pone.0196607
43. Oldham S, Fornito A. The development of brain network hubs. Dev Cogn Neurosci. (2019) 36:100607. doi: 10.1016/j.dcn.2018.12.005
44. Marleau I, Vona M, Gagner C, Luu TM, Beauchamp MH. [Formula: see text] social cognition, adaptive functioning, and behavior problems in preschoolers born extremely preterm. Child Neuropsychol. (2021) 27:96–108. doi: 10.1080/09297049.2020.1797656
45. Padilla N, Eklof E, Martensson GE, Bolte S, Lagercrantz H, Aden U. Poor brain growth in extremely preterm neonates long before the onset of autism spectrum disorder symptoms. Cereb Cortex. (2017) 27:1245–52. doi: 10.1093/cercor/bhv300
46. Danielson NB, Guo JN, Blumenfeld H. The default mode network and altered consciousness in epilepsy. Behav Neurol. (2011) 24:55–65. doi: 10.1155/2011/912720
47. Menon V, Uddin LQ. Saliency, switching, attention and control: a network model of insula function. Brain Struct Funct. (2010) 214:655–67. doi: 10.1007/s00429-010-0262-0
48. 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) 185:654–63. doi: 10.1016/j.neuroimage.2018.07.021
49. Lee W, Al-Dossary H, Raybaud C, Young JM, Morgan BR, Whyte HE, et al. Longitudinal cerebellar growth following very preterm birth. J Magn Reson Imaging. (2016) 43:1462–73. doi: 10.1002/jmri.25098
50. Limperopoulos C, Chilingaryan G, Sullivan N, Guizard N, Robertson RL, du Plessis AJ. Injury to the premature cerebellum: outcome is related to remote cortical development. Cereb Cortex. (2014) 24:728–36. doi: 10.1093/cercor/bhs354
51. Strahle JM, Triplett RL, Alexopoulos D, Smyser TA, Rogers CE, Limbrick DD Jr, et al. Impaired hippocampal development and outcomes in very preterm infants with perinatal brain injury. Neuroimage Clin. (2019) 22:101787. doi: 10.1016/j.nicl.2019.101787
52. Giménez M, Junqué C, Narberhaus A, Caldú X, Salgado-Pineda P, Bargalló N, et al. Hippocampal gray matter reduction associates with memory deficits in adolescents with history of prematurity. Neuroimage. (2004) 23:869–77. doi: 10.1016/j.neuroimage.2004.07.029
53. Lemola S, Oser N, Urfer-Maurer N, Brand S, Holsboer-Trachsler E, Bechtel N, et al. Effects of gestational age on brain volume and cognitive functions in generally healthy very preterm born children during school-age: a voxel-based morphometry study. PLoS ONE. (2017) 12:e0183519. doi: 10.1371/journal.pone.0183519
Keywords: magnetic resonance imaging, extremely preterm, voxel-based morphometry, intraventricular hemorrhage, patent ductus arteriosus, gray matter volume, white matter volume
Citation: Kvanta H, Bolk J, Broström L, Fernández de Gamarra-Oca L, Padilla N and Ådén U (2023) Extreme prematurity and perinatal risk factors related to extremely preterm birth are associated with complex patterns of regional brain volume alterations at 10 years of age: a voxel-based morphometry study. Front. Neurol. 14:1148781. doi: 10.3389/fneur.2023.1148781
Received: 20 January 2023; Accepted: 28 April 2023;
Published: 19 May 2023.
Edited by:
Wenbin Guo, Central South University, ChinaReviewed by:
Dustin Scheinost, Yale University, United StatesGauthier Loron, Reims University, France
Copyright © 2023 Kvanta, Bolk, Broström, Fernández de Gamarra-Oca, Padilla and Ådén. 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: Hedvig Kvanta, aGVkdmlnLmt2YW50YSYjeDAwMDQwO2tpLnNl