Weak functional connectivity in the human fetal brain prior to preterm birth.

It has been suggested that neurological problems more frequent in those born preterm are expressed prior to birth, but owing to technical limitations, this has been difficult to test in humans. We applied novel fetal resting-state functional MRI to measure brain function in 32 human fetuses in utero and found that systems-level neural functional connectivity was diminished in fetuses that would subsequently be born preterm. Neural connectivity was reduced in a left-hemisphere pre-language region, and the degree to which connectivity of this left language region extended to right-hemisphere homologs was positively associated with the time elapsed between fMRI assessment and delivery. These results provide the first evidence that altered functional connectivity in the preterm brain is identifiable before birth. They suggest that neurodevelopmental disorders associated with preterm birth may result from neurological insults that begin in utero.

Developmental problems are significantly more common in children born preterm. Epidemiological and meta-analytic studies indicate children born preterm are three times more likely to develop autism, attention deficit/hyperactivity disorders, and emotional disorders12, five times more likely to manifest neurological abnormalities3, and three to four times more likely to experience school failure4. Studies of the neuropathology of brain injury in the premature infant have identified alterations in cerebral white matter, grey matter, and projection fibers of the brain in infants born preterm56. Neural irregularities have been observed at the systems level as well, with altered wide-scale network connectional architecture identified in infants born preterm78910. These studies report diminished coherence of activity measured across brain circuits, indicative of weaker connectivity in individuals born preterm from infancy1112 through adulthood1314. Given the robust linkage between child neurological disorders and atypical large-scale neural connectivity1516 (see also refs 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37), theory has evolved that neurodevelopmental impairment following preterm birth may stem from alterations in neural connectivity38, and recent evidence provides nascent support for this claim1039.

A critical target for current research is to separate intrauterine from extrauterine influences on injurious brain development, and to isolate the earliest indicators of change in connectional architecture of the preterm brain. While much has been gained from examination of the postnatal preterm brain, examinations of human brain networks at or after preterm birth are confounded by potential insults conferred both by the absence of neuroprotective elements4041 and addition of neurotoxic influences, which are inherent conditions of early delivery4243. There has been an emphasis on hypoxia-ischemia and infection/inflammation as upstream etiologies of preterm brain injury, but there are several less frequently considered factors that also influence brain development, including extrauterine respiratory complications and the effects of fetal deprivation from maternal hormones and nutritional factors41. Studies of functional neural connectivity prior to preterm birth are needed to isolate processes that begin in the womb. If functional connectivity is altered in the preterm brain in utero, the untoward influences of extrauterine factors cannot be the source of those differences. The recent development of resting-state fMRI (rs-fMRI) methodology for the human fetus offers the first opportunity to investigate altered functional connectivity prior to birth4445.

In this study we utilized rs-fMRI to measure neural connectivity in utero in 32 human fetuses, 14 of which were subsequently born preterm. We employed an intrinsic connectivity distribution analysis of rs-fMRI data to map synchrony in MRI signals over time across the fetal brain. We then investigated the relationship between connectivity and gestational age at fMRI assessment and delivery. Lastly, post-hoc seed analysis was performed to explore the regions in which specific connections were most responsible for changes intrinsic connectivity distribution (ICD)46 value differences between preterm and full-term groups.

Results

Description of research cohort

Thirty-six women, mean age 25.3 years, SD, 5.6, underwent MRI between their 22nd and 36th week of pregnancy. Half of participating pregnant women were high-risk for early delivery and gave birth prior to the 37th week of pregnancy, mean GA preterm birth = 32 weeks, range = 24–35 weeks. The comparison group was case control matched based on GA at time of MRI and gender of the fetus (Table 1). Four high-risk participants were excluded due to fetal intrauterine growth restriction, or IUGR, which can influence neural network connectivity47. The final study sample consisted of 14 pregnancies that ended in preterm delivery between 24 and 35 weeks, and 18 uncomplicated pregnancies. Birth and placental pathology outcomes of the preterm sample are provided in Supplemental Table 1.

Table 1

Summary of participant and fMRI data characteristics by group.

	Fetuses (n = 14) later born prematurely	Fetuses (n = 18) later born at term
Fetal GA at scan	29.6 (3.7)	30.3 (3.2)
Fetal GA at delivery*	32.4 (4.2)	38.8 (1.1)
Maternal depressive symptoms, m (SD)	15.3 (12.1)	12.4 (8.5)
Maternal anxiety symptoms, m (SD)	36.9 (12.7)	34.1 (8.2)
Maternal stress, m (SD)
Penn State Worry Questionnaire	19.9 (9.0)	14.9 (5.9)
Perceived Stress Scale*	46.1 (14.8)	40.2 (8.7)
Fetal gestational age at MRI, m (SD)	29.7 (3.7)	30.25 (3.3)
Maternal IQ, m (SD)	82.8 (13.8)	78.8 (13.1)
Ethnicity, n (%)
Caucasian	0	1 (5.6)
African-American	13 (92.9)	14 (77.8)
Asian	0	1 (5.6)
Other	1 (7.1)	2 (11.1)
Education, n (%)
No GED/High-school diploma	3 (21.4)	6 (33.3)
GED/High-school diploma	3 (21.4)	5 (27.8)
Some college	6 (42.9)	6 (33.3)
2-yr college degree	1 (7.1)	1 (5.6)
4-yr college degree	1 (7.1)	0
Annual Income, n (%)
<$10,000	9 (64.3)	6 (35.5)
$10,000–$20,000	2 (14.3)	7 (41.2)
$20,000–$30,000	2 (14.3)	3 (17.6)
$30,000–$40,000	0	1 (5.9)
>$40 000	1 (7.1)	0
Motion during resting-state scan
Translational mean movement, m (SD)	0.4 (0.1)	0.4 (0.1)
Rotational mean movement, m (SD)	0.7 (0.3)	0.7 (0.4)
Translational RMS, m (SD)	0.2 (0.1)	0.3 (0.1)
Rotational RMS, m (SD)	0.0 (0.0)	0.0 (0.0)
fMRI data characteristics
SAR, m (SD)	0.26 (0.1)	0.3 (0.1)
Number of fMRI frames analyzed, m (SD)	171 (81.2)	185.4 (50.4)
Proportion frames retained after exclusion of periods of movement, % (SD)	0.6 (0.2)	0.5 (0.1)

Chi-square tests compared race/ethnicity, education, and income between groups. Two-sample independent t-tests compared all other variables. As planned, fetal GA at delivery* was significantly different, p < 001. All other comparisons were non-significant, using two-tailed p < 0.05. A single trend was observed for the Perceived Stress Scale** at p = 0.066; p’s > 0.1 for all other comparisons. Depressive and anxiety symptoms were measured using the Center for Epidemiologic Studies Depression Scale and State-Trait Anxiety Inventory, respectively. Intelligence Quotient (IQ) was measured using verbal and matrix reasoning subtests of the Wechsler Abbreviated Scale of Intelligence. Gestational age (GA) reported in weeks; translational (x, y, z) movement reported in mm; rotational, in degrees; Specific Absorbtion Rate (SAR) in units of watts per kilogram (W/kg). Abbreviations: standard deviation, SD; mean, m; root-mean-square, RMS (head position change).

Studies that demonstrate postnatal differences in brain development as a function of sex48, maternal prenatal stress49, and socioeconomic status50 introduce the possibility that intrauterine functional connectivity may differ as a result of such factors. We examined the possibility that these or other key indicators may differ between our study groups. No differences were observed between study groups in fetal age or gender, maternal demographic or IQ measures, nor fMRI data quality measures (Table 1). From this we infer that group differences do not result from features of the social and psychological context that may differ between groups or from data quality, but rather the altered brain development in fetuses later born preterm.

Diminished connectivity in fetuses later born preterm

ICD comparisons of preterm- and term-born fetuses revealed a substantial area of the left hemisphere, proximal to what will later become Broca’s area, where connectivity was greater in term-born fetuses (Fig. 1). In contrast, there were no areas in which connectivity was greater in those later born preterm. Observed effects agree with evidence of atypical brain connectivity in preterm born neonates obtained using similar MRI systems level analyses111251, and animal electrophysiological and cellular models of prematurity and inflammation that report diminished dendritic branching, arborization, oligodendrocyte maturation, and cell-to-cell communication in affected animals52. A follow up analysis examined the possibility that preterm labor unduly contributed to group differences, and found that even after removing 4 fetuses born within 7 days of the scan, the effect remained significant (Supplemental Fig. 1). This effect also remained significant when repeated in a multiple linear regression model with age, sex, and motion variables as covariates (Supplementary Table 2). Differences in connectivity in this putative pre-language region support the notion that language impairments and related altered connectivity often observed in individuals born preterm535455 may arise from altered brain development that commences in utero.

Figure 1

Comparison of preterm (PT) - and term (T) -born fetuses using voxel-level connectivity.

Single group maps for (A) term and (B) preterm-born fetuses, showing putative hub regions in the fetal brain at 29.6 weeks. Significant differences (p < 0.05, corrected) between preterm- and term-born fetuses were observed in eventual left hemisphere language regions in the frontal lobe as shown on (C) surface rendering, and (D) axial slices. No regions exhibited increased ICD for preterm-born fetuses compared to term born fetuses. For the preterm-born fetuses, ICD values extracted from a seed centered on the peak difference (27 mm left, 18 mmanterior, 5 mm inferior of the center of image data) were significantly correlated with both gestational age (GA; panel (E)) at scan (r = 0.69, p = 0.003, df = 12) and (F) at birth (r = 0.51, p = 0.03, df = 12) using one-tailed Spearman rank correlation.

Neural connectivity relates to age at MRI and length of gestation

We reasoned that if this fetal pre-language region is a critical hub of connectional group differences, then variation in connectivity of this region may relate to individual differences in the preterm-born sample. Correspondence between functional connectivity, or ICD, of this region, and both gestational age at the time of scan and gestational age at delivery were evaluated. Generally, postnatal evaluation of the human premature brain has pointed to early developmental delays, some of which improve with age, either through restorative or compensatory processes56. We therefore expected that gestational age at time of scan may relate to the level of connectivity in this region, with older preterm-born fetuses demonstrating connectivity values more similar to the term-born controls. Results supported this idea, with positive correspondence between rank age at scan and values in the ICD peak region in the preterm group, r = 0.69, p = 0.003 (Fig. 1E). An analogous examination of the association between neural functional connectivity and gestational age at delivery revealed, again, a significant positive correlation in preterm fetuses within the ICD peak region, r = 0.51, p = 0.03 (Fig. 1F). This latter observation suggests that more significant impairments observed in those born extremely preterm may begin with altered neural functional connectivity before birth. Together, these findings demonstrate that the extent to which this pre-language region is connected with other brain regions, rather than functioning independently, is related to both gestational age at scan and length of gestation. Those destined to have longer gestational duration, were more similar in functional connectivity to those born at term.

Altered proto-linguistic functional connectivity

To better understand differences between groups, we performed follow-up region of interest seed connectivity analysis using this putative pre-language, ICD peak region as an area of interest. Signal from this region was correlated with signal intensity in every other grey matter voxel to identify brain regions that demonstrated similar signal properties over time. In the preterm-born fetal group, we observed diminished signal correlation in ipsilateral posterior superior temporal gyrus, extending into the posterior pre-language brain network (Fig. 2; Supplementary Table 2), when compared to term-born control cases. Previously, we have shown in typically developing fetuses that a large decrease in modularity between these regions occurs during the third trimester57, indicative of increased crosstalk between these regions with advancing age. Our current results indicate that functional organization of the larger pre-language system, encompassing both Broca’s and Wernicke’s areas, may be reduced in fetuses who subsequently are born preterm.

Figure 2

Comparison of preterm- and term-born fetuses using seed connectivity.

(A) To investigate which specific connection may be driving our voxel-level ICD results, we performed follow-up seed connectivity using a cubic seed (shown in red) centered at the peak coordinate of between-group differences (27 mm left, 18 mm anterior, 5 mm inferior of the center of image data). (B) Significant differences (p < 0.05, corrected) between preterm- and term-born fetuses were observed in regions eventually becoming primary auditory cortex and Wernicke’s area as shown on (B) surface rendering and (C) axial slices. No regions exhibited increased seed connectivity for preterm-born fetuses compared to term-born fetuses.

Remaining gestational duration linked to contralateral connectivity

We also examined connectivity of the area of peak differences between groups, the anterior pre-language region, within the preterm fetal group to ascertain whether connectivity profiles were different in fetuses who would be born soon after their MRI date, versus those born further out from the time of MRI. Whole-brain regression from the peak ICD region revealed that connectivity to right-hemisphere pre-language homologues was diminished in those fetuses that would be born soon after MRI, defined as duration of time between scan and delivery. This effect remained significant even when controlling for age at scan and age at delivery (Fig. 3). This is striking given knowledge from our prior work that connectivity to contralateral hemisphere neural homologues increases with advancing fetal age44. This new observation suggests that fetuses at risk for impending preterm delivery had weaker left to right hemisphere connectivity in this lateral fronto-parietal region, the region that differentiated preterm versus term born whole-brain functional connectivity. This altered connectivity between hemispheres in utero may set the foundation for altered cross-hemisphere connectivity observed in infants51, children55, and adolescents5354 born preterm.

Figure 3

Whole brain regression of time between scan and birth for fetuses (n = 14) subsequently born preterm.

Using the same seed connectivity as in Fig. 2, we investigated whether specific connections to the seed region were associated with increased time between fetal MRI and birth. (A) Connectivity between the left hemisphere proto-language seed and right hemisphere language homologues in the parietal lobes were significantly correlated (p < 0.05, corrected) with increased time between fetal MRI and birth. No regions exhibited a significant negative correlation. This association suggests that cross-hemisphere connectivity between language regions may be predictive of longer in utero development for those at risk of preterm birth. Connectivity values averaged from a seed centered on the peak difference and linear fit from spearman rank correlation are extracted for each subject and used to visualize the observed effect in (B). Using partial spearman correlation, these correlations remained significant after controlling for gestational age (GA) at scan (r = 0.78, 0 = 0.002, df = 11) and GA at birth (r = 0.83, p = 0.001, df = 11), suggesting a unique effect of time between fetal MRI and birth.

Underdetermined role of inflammation in altered neural functional connectivity

Conditions that increase risk for prematurity may also be those that enhance risk for deficient neural functional connectivity of the fetus. Intra-amniotic infection increases risk for both preterm delivery and white matter injury5859. The association between intra-amniotic infection and neural connectivity is as yet untested in humans, but is indirectly supported by postnatal data demonstrating links between white matter injury and neural functional connectivity measured in premature infants60. Placental pathology reports obtained for 28 participants in our study sample were evaluated for presence or absence of inflammatory lesions. Pathology confirmed high prevalence of inflammatory lesions, either acute or chronic, manifested as focal or generalized acute chorioamnionitis, chronic chorioamnionitis, and funisitis in the study sample (Supplemental Table 1). Indicators of inflammation were present in placental histology reports of all women in the preterm group. As expected there were significantly fewer indications of inflammatory lesions in the comparison group, χ2(2, N = 28) = 8.089, p = 0.007. These clinical variables may relate to neural functional connectivity, but due to small sample size and potential for type II error, it was not appropriate to examine those potential associations here.

Discussion

Defining the nascent architecture of the preterm brain provides a basis for understanding the underlying etiology of neurological impairments that can accompany prematurity. Using recent developments in fetal resting-state fMRI we examined neural functional connectivity in 32 human fetuses and found reduced connectional integrity in fetuses that would subsequently be born preterm, particularly in regions of the left hemisphere that later support language processing. Strength of functional connectivity in this region was related to gestational age at delivery, such that those born closer to expected due date demonstrated connectivity profiles more similar to those in the case-matched control group. Thus, observed differences in neural functional connectivity were related to proximal health outcomes in the preterm group.

These results demonstrate that neurological connectivity differences associated with human preterm birth begin in utero, prior to the potentially injurious experiences of early delivery. These constitute the first human data to suggest that disabilities frequently accompanying extreme prematurity, such as autism and ADHD, may derive from pre-existing intrauterine neurological conditions, especially given that these disorders have neuroconnectional bases61626364. We found that diminished functional connectivity was also linked to individual differences in the preterm group, suggesting that maturation of brain circuitry may be associated with the circumstances of an individual pregnancy. We observed that both gestational age at time of scan and gestational age at delivery positively relate to the degree of connectivity, such that fetuses with increased intrauterine course demonstrate connectivity values more similar to the term-born controls. This pattern builds on prior multi-level postnatal neuroscience findings demonstrating that brain circuits flexibly adjust across the life-span and in response to injury, environmental programming, and/or disease6566. Similarly, fetuses born preterm often overcome initial functional deficits during postnatal or early childhood development. Our results indicated that left-hemisphere pre-lingual regions and cross-hemispheric connections were most significantly impacted, fitting with the language deficits that frequently manifest in those born extremely preterm6768. This observation may also imply centrality of lingual brain networks in organization of emergent functional brain systems.

Knowledge that functional connectivity differs before birth in those born preterm encourages evaluation of prenatal intervention. Potential efficacy of prenatal neurobehavioral therapy is supported by growing evidence from non-human animals and from humans that learning occurs prior to birth. Markham and colleagues have demonstrated that birds learn about the calls of conspecifics while still in the egg, and that this embryonic programming alters neural cellular activation after birth and also characteristics of later vocalizations69. A landmark study by Partanen and colleagues presented evidence of experienced-based neuroplasticity in human fetuses. This group exposed fetuses to varied quantities of word-like sounds and found that, after birth, those exposed to intrauterine stimuli showed behavioral learning and enhanced brain activity that scaled with quantity of exposure70. They concluded that prenatal exposure to complex sounds may lead to the development of a more effective neural network for information processing after birth. These studies substantiate presence of cross-species learning-induced neural plasticity in utero. Understanding that extrauterine environmental stimuli not only induce behavioral response in the fetus, but also change brain function in lasting ways, supports the possibility that prenatal behavioral intervention may serve as an effective therapy for those at risk.

While this study presents a novel comparison of human brain function before birth in fetuses destined either for term or preterm delivery, considerations about nascent fetal fMRI methodology merit discussion. Challenges inherent in the methodology include small fetal head volume, influences of physiological signals originating from mother and fetus, limited constraints over motion, and variation in orientation71. These concerns are not entirely unique to fetal fMRI, as partial volume effects, physiological noise, motion, and image registration are broad concerns for the larger MR field, and many strategies exist for addressing these72737475, several having been applied in the present study. However, the fetus represents an extreme case for these crucial areas and considerable work remains to be done to bring fetal fMRI to its full potential, including such things as development of specialized tools and atlases for fetal MRI. Until that time, current best practices include normalization to age-appropriate fetal templates76, use of subject specific anatomical segmentation, and use of data-driven, rather than spatially constrained, analytic approaches such as independent components analyses77, spatial-spectral parcellation5778, ICD46 (employed here), and multivariate distance-based analyses79 methods.

Another topic requiring careful consideration in fetal fMRI is that we possess limited understanding of the physiological basis of blood-oxygen level dependent (BOLD) fMRI signals in the fetal brain. Current understanding is derived from assumptions regarding neurovascular coupling that stem from neurophysiological research performed in animals80. However, emergent data suggests some of these assumptions may be justified even in the immature human fetal brain. For example, studies are beginning to show that fetal functional connectivity MRI measures are congruous with what we know about neural development during this time. We and others have shown that intrahemispheric44 and long-range8182 fetal fMRI signals become more synchronized with age, which mirror known principles of fetal anatomical development8384. In addition, the few studies that have investigated BOLD signal in the antenatal period report positive BOLD contrast responses in preterm infants85 and in animal models of the preterm period86. Furthermore, murine studies of pre- and postnatal angioarchitechtonics support tight coupling of neural and vascular dynamics across early development87. Overall, it is well reasoned that signal covariation measured in fetal resting-state studies reflects the establishment of communicative architecture of the brain. While support for this is scarce and interpretive caution is advisable, available data do support this position.

In conclusion, we provide the first evidence that neural pathways are likely to be altered prior to preterm delivery. This discovery suggests that factors influencing early delivery may also impact development of the human brain, which has implications for life-long health. Future work will address sources such as infection and inflammation that may play a causal role in altering these parallel pathways, bringing us closer to understanding both the primary neurological injury and the optimal timing for early intervention.

Methods

Participants

Singleton pregnancies with normal brain anatomy assessed by ultrasound and MRI examination reporting no contraindications for MRI were eligible to participate. MRI T1 weighted images revealed no brain injury in any fetal cases at the time of prenatal MRI examination. Despite that no areas suggestive of brain lesions were observed, the presence of micro bleeds, or other more subtle forms of injury, cannot be completely excluded. The mean age of fetuses at the time of MRI was 29.9 weeks, SD, 3.4, post-menstrual gestational age (GA). Ultrasound (US) examination administered by study physician (E.H.-A.) was performed within 1 week of MRI examination to determine fetal GA. All women were native English speakers. All women provided written informed consent before undergoing MRI examination. Participation was approved by the Institutional Review Board of the National Institute of Child Health and Human Development (NICHD) and by the Human Investigation Committee of Wayne State University. All experiments were performed in accordance with relevant guidelines and regulations.

Image acquisition

Fetal MR exams were performed with a Siemens Verio 70-cm open-bore 3-T scanner with a light-weight (~550 g) abdominal 4-Channel Siemens Flex Coil. The MRI examination lasted 45 min. Images were collected with the following parameters and computed SAR values: (i) localizer [repetition time/echo time (TR/TE), 20/4.2 ms; 10-mm thickness; SAR = 0.21]; (ii) T2 anatomical (TR/TE, 3500/140 ms; 3-mm slice thickness; repeated 4–6 times; SAR = 0.53); (iii) echo planar imaging (EPI) BOLD (TR/TE, 2000/30 ms; 180 frames; 4-mm slice thickness; axial; repeated twice; SAR = 0.3); (iv) susceptibility weighted imaging (SWI) (TR/TE, 30/21.2 ms; 4-mm slice thickness; SAR = 0.06).

Connectivity preprocessing

Fetal fMRI data were censored75 for motion using criterion <1 mm frame-to-frame translation and <1.5 degrees rotation. 56% of data collected were retained after motion censoring. Censored data were motion corrected, normalized to common template space, and smoothed as previously described57 using both manual and automatic methods. Further connectivity analysis was performed using BioImage Suite88. Several covariates of no interest were regressed from the data including linear and quadratic drifts, 6 motion parameters, mean cerebral-spinal-fluid (CSF) signal, and mean white-matter signal. The data were temporally smoothed with a zero mean unit variance Gaussian filter (approximate cutoff frequency = 0.12 Hz). A gray-matter mask defined in template space was applied to the data so only gray matter voxels were used in further calculations.

Intrinsic functional connectivity

After preprocessing, intrinsic connectivity distribution (ICD; Fig. 4) was computed at the voxel level for each subject as described previously46. Similar to most voxel-based functional connectivity measures, ICD involves correlating the time-course for any voxel with every other time-course in the brain or brain hemisphere, and then a summary statistic based on the network theory measure degree was calculated. ICD models the entire distribution of correlation thresholds using a Wiebull distribution avoiding the need for choosing an arbitrary connectivity threshold. This parameterization is akin to modeling the change in network theory metric degree, as the threshold used to calculate degree is increased, with a stretched exponential. Specifically, the time-course for any gray matter voxel was correlated with every other voxel in the gray matter. A histogram of these correlations was constructed to estimate the distribution of connections to the current voxel. This distribution was converted to a survival function and the survival function was fitted with a stretched exponential with unknown variance. As variance controls the spread of the distribution of connections, a larger variance indicates a greater number of high correlation connections. Finally, this process is repeated for all voxels in the gray matter resulting in a whole-brain parametric image summarizing the connectivity of each tissue element. ICD was computed in each study group (Supplemental Fig. 2) and compared to identify brain regions where functional connectivity differed between groups.

Figure 4

Overview of ICD analysis.

For any voxel in the gray matter, the voxel’s time course is correlated with every other voxel’s time course. This procedure is repeated for every voxel, resulting in a voxel-by-voxel correlation matrix. From this matrix, a single row is extracted (representing all correlation to a voxel) and converted to a histogram to estimate the distribution of connectivity for that voxel. From this distribution, a survival function is constructed and parameterized with a stretch exponential with unknown parameters α and β. For our purposes, α controls the rate of decay of the survival function with a larger α indicating a slower decay and larger global connectivity. The α ICD maps can then be thresholded to reveal putative hubs of the fetal connectome.

To interrogate relative differences in connectivity, each participant’s map was normalized by subtracting the mean across all voxels and dividing by the standard deviation across all voxels. This z-score-like normalization does not change the underlying connectivity pattern but allows for investigation of relative differences in connectivity in the presence of large global differences in connectivity89.

Follow-up seed connectivity

Follow-up seed analysis was performed to explore (post-hoc) connectivity of the node where ICD connectivity differed most significantly between groups. Signal intensity was extracted from voxels centered on the ICD peak at location 27 mm left, 18 mm anterior, 5 mm inferior of the center of image data registered to a common 32-week fetal brain template76. The time course of the seed region in a given participant was then computed as the average time course across all voxels, comprising a 512 mm3 cube, in the seed region. This time course was correlated with the time course for every other voxel in the gray matter to create a map of r-values, reflecting seed-to-whole-brain connectivity. These r-values were transformed to z-values using Fisher’s transform yielding one map for each participant representing the strength of correlation to the seed region.

Statistical analysis

For imaging data, voxel-wise two sample independent t-tests were used to compare the connectivity data between groups. Voxel-wise Pearson’s correlation was used to assess association between connectivity and GA. Spearman rank correlation and partial Spearman rank correlation were used to assess the association between extracted connectivity values and GA. Results are shown at a cluster-level threshold of p < 0.05 family-wise error (FWE) correction as determined by AFNI’s 3dClustSim program. The cluster threshold was determined using cluster-forming threshold of p = 0.001, 10,000 iteration, smoothness estimated using a mixture of Exponential and Gaussian distributions (i.e. the –acf option), and 16.0.09 release of AFNI. For non-imaging data, two sample independent t-tests and chi-square tests were used with significances assessed at p < 0.05.

Additional Information

How to cite this article: Thomason, M. E. et al. Weak functional connectivity in the human fetal brain prior to preterm birth. Sci. Rep. 7, 39286; doi: 10.1038/srep39286 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.