Abstract
During the second-to-third trimester, the neuronal pathways of the fetal brain experience rapid development, resulting in the complex architecture of the interwired network at birth. While diffusion MRI-based tractography has been employed to study the prenatal development of structural connectivity network (SCN) in preterm neonatal and postmortem fetal brains, the in utero development of SCN in the normal fetal brain remains largely unknown. In this study, we utilized in utero dMRI data from human fetuses of both sexes between 26 and 38 gestational weeks to investigate the developmental trajectories of the fetal brain SCN, focusing on intrahemispheric connections. Our analysis revealed significant increases in global efficiency, mean local efficiency, and clustering coefficient, along with significant decrease in shortest path length, while small-worldness persisted during the studied period, revealing balanced network integration and segregation. Widespread short-ranged connectivity strengthened significantly. The nodal strength developed in a posterior-to-anterior and medial-to-lateral order, reflecting a spatiotemporal gradient in cortical network connectivity development. Moreover, we observed distinct lateralization patterns in the fetal brain SCN. Globally, there was a leftward lateralization in network efficiency, clustering coefficient, and small-worldness. The regional lateralization patterns in most language, motor, and visual-related areas were consistent with prior knowledge, except for Wernicke's area, indicating lateralized brain wiring is an innate property of the human brain starting from the fetal period. Our findings provided a comprehensive view of the development of the fetal brain SCN and its lateralization, as a normative template that may be used to characterize atypical development.
Significance Statement
We studied the normal development of intrahemispheric corticocortical structural connectivity networks (SCNs) of the fetal brain from 26 to 38 gestational weeks using in utero diffusion MRI data. Graph theory-based analysis revealed significant enhancement in network efficiency and clustering, as well as persistent small-worldness with age, revealing balanced integration and segregation in the fetal brain SCN during the studied period, supported by regional developmental patterns. Leftward lateralization in network efficiency, clustering coefficient, and small-worldness was observed. Regional lateralization patterns in most language, motor, and visual-related areas were consistent with prior knowledge. We also summarized the challenges of investigating the fetal brain SCN development and provided suggestions for future studies.
Introduction
The fetal brain neuronal pathways rapidly develop during the second-to-third trimester. Histological studies demonstrated the growth of thalamocortical fibers (Kostovic and Jovanov-Milosevic, 2006), the ingrowth of cortical axons, and the emergence of corticocortical fibers with primary gyrification (Takahashi et al., 2014). These fiber pathways interconnect different cortical regions, forming complex brain networks (Sporns et al., 2005). Altered fetal brain structural connectivity are linked with neurodevelopmental abnormalities (Jakab et al., 2015), underscoring the importance of characterizing normative structural connectivity development.
Diffusion magnetic resonance imaging (dMRI) is a powerful neuroimaging tool for the noninvasive characterization of neuronal microstructure and brain injuries (Assaf et al., 2019). It also enables the reconstruction of white matter (WM) fiber pathways via tractography (Jeurissen et al., 2019). Several early studies have explored fetal brain WM development using ex vivo dMRI tractography of brain specimens (Huang et al., 2006, 2009; Takahashi et al., 2012). Recent advances in fast imaging and processing techniques have facilitated in utero studies of fetal brain WM fiber bundles (Hooker et al., 2020; Jaimes et al., 2020; Wilson et al., 2021) and spatiotemporal atlases depicting their development (Khan et al., 2019; Chen et al., 2022).
In addition, dMRI-based fiber tractography allows quantifications of whole-brain structural wiring patterns. By defining brain regions as nodes and traced streamlines as edges, a structural connectivity network (SCN) can be established. SCN topology can be quantitatively measured by graph theory metrics. Global efficiency (Eglob) quantifies the exchange of information across the entire network, reflecting network integration. Local efficiency (Eloc), clustering coefficient (CC), and modularity reflect network segregation. Eloc also quantifies a network's resistance to localized failures. Small-worldness (SW) characterizes a network with both high clustering and short paths, indicating balanced integration and segregation. Details on network measures are provided below, Network properties. Previous studies have revealed integration and segregation processes in prenatal brain SCN using preterm neonatal (Brown et al., 2014; Zhao et al., 2019a) and postmortem fetal brain data (Song et al., 2017), highlighting increases in Eglob, Eloc, and clustering coefficient (CC), and the emergence of small-world and rich-club organizations (Ball et al., 2014; Brown et al., 2014; Song et al., 2017; Zhao et al., 2019a). Fetal brain functional connectivity network (FCN) studies also indicated segregation patterns such as increasing CC (Cao et al., 2017), the existence of network communities (Turk et al., 2019), and integration patterns involving strengthened intermodule connections (Thomason et al., 2014). Complex functional topologies such as small-world, rich-club, and functional gradients were also observed (Fransson et al., 2011; Thomason et al., 2014; Turk et al., 2019; Moore et al., 2023).
While some studies used in utero fetal brain dMRI data to construct SCNs (Jakab et al., 2015; Hunt et al., 2018), they focused on specific conditions or comparative analyses, lacking descriptions of temporal changes or integration–segregation tendencies in the normal fetal brain SCN. In this work, we hypothesized the existence of a balanced integration and segregation process during the in utero development of fetal brain SCN and tested this hypothesis by measuring the developmental trajectories of network properties indicative of these tendencies.
Exploring the asymmetry of the developing brain network is another intriguing area. Brain structural asymmetry offers insights into hemisphere-specific functions (Toga and Thompson, 2003). Previous research suggested stronger left hemisphere connectivity in neonatal brains (Ratnarajah et al., 2013), especially relevant for language and motor skills. Resting-state functional MRI highlighted stronger connections in the left language and auditory networks in fetuses between 18 and 37 gestational weeks (W), while the visual and attention networks exhibited rightward lateralization (Kasprian et al., 2011). However, the lateralization patterns of dMRI-based fetal brain SCN were not investigated.
This study aims to establish SCN using in utero fetal brain dMRI data from 26W to 38W. We focused on intrahemispheric corticocortical SCN and seek to assess the coexistence of network integration and segregation, characterize regional variations, and delineate early lateralization patterns in the fetal brain SCN.
Materials and Methods
Experimental design
In this study, we constructed SCN of the fetal brain using in utero dMRI data acquired from 26W to 38W (see below, dMRI data processing and Building connectome). Graph theory analysis was used to evaluate the network properties of the fetal brain SCN, such as the nodal strength and betweenness, global and local efficiency, modularity, shortest path length, clustering coefficient, and small-worldness (see below, Network properties). We then performed correlation analyses between these measures and GA to study the developmental changes in fetal brain SCN. In addition, we investigated the asymmetry of SCN and its developmental pattern (see below, Statistical analysis).
Subjects and data acquisition
We enrolled 171 healthy pregnant women at Shandong Provincial Hospital after obtaining approval from the Institutional Review Board. The enrollment criteria included the following: (1) gestational age (GA) between 24 and 38 weeks; (2) no maternal comorbidities; and (3) fetuses with no clinical, ultrasound, or MRI evidence of brain abnormalities. Exclusion criteria were as follows: (1) pregnant women with gestational diabetes mellitus or hypertensive disorder complicating pregnancy, (2) multiple pregnancies, and (3) fetuses with chromosomal or genetic abnormalities diagnosed by amniocentesis. The scans were performed on a 3 T MR scanner (MAGNETOM Skyra, Siemens Healthineers) using an 18-channel body coil. The pregnant women were in the supine position when the fetuses were being scanned. We adopted a diffusion-weighted echoplanar imaging sequence with a b-value of 600 s/mm2, 30 gradient directions, eight nondiffusion-weighted images, 1.73 mm × 1.73 mm in-plane resolution in axial slices, 4 mm slice thickness, and two averages. Other acquisition parameters were as follows: field of view (FOV), 260 × 260 mm2; repetition time (TR), 3,900 ms; and echo time (TE), 80 ms. GeneRalized Auto-calibrating Partial Parallel Acquisition (GRAPPA) parallel imaging with an acceleration factor of 2 was adopted. Acquisition time for each subject was 3–4 min. Quality checks were performed before and after data processing by visual inspections to exclude images with poor data quality. This ensured that only high-quality images were used in the study. The data inclusion and exclusion flowchart was shown in Figure 1A.
dMRI data processing
The dMRI data processing was conducted using MRtrix3 and SVRTK software. The raw data underwent a preprocessing pipeline in MRtrix3 (https://www.mrtrix.org/), including denoising, volume-to-volume alignment, correction of eddy-current–induced distortions, and bias field correction (Tournier et al., 2019). The images were further processed using SVRTK (https://github.com/SVRTK/SVRTK; Deprez et al., 2020). The pipeline first performed iterative slice-to-volume registration (SVR) and superresolution reconstruction of all 2D image slices, regardless of gradient directions, to generate a motion-free mean diffusion-weighted (DW) 3D volume at 1.2 mm isotropic resolution (Kuklisova-Murgasova et al., 2012), which was used to initialize the spherical harmonics model. The original DW slices along with their slice-specific motion parameters resulting from SVR were then fitted to reconstruct all DW volumes with the corrected gradient tables (Deprez et al., 2020).
Then the FODs were calculated using constrained spherical deconvolution in MRtrix3 (Tournier et al., 2007). Probabilistic tractography was performed based on the FODs using a signal cutoff of 0.08 and a maximum angle threshold of 30° (Tournier et al., 2010; Chen et al., 2022). Given the considerable variability in fetal brain sizes during the studied period, we randomly selected 10 seeds per voxel to initiate tracking, which ensured a uniform seeding density across different brains. Subsequently, we applied SIFT2 to the resulting tractograms to match the distribution of streamlines to the underlying FOD signals and mitigate false positives (Smith et al., 2015).
Building connectome
Thirty-nine cortical regions in each hemisphere were defined based on the CRL fetal brain parcellations (Gholipour et al., 2014). First, the b0 image of each subject was nonlinearly registered to the CRL fetal brain T2W atlas at the corresponding GA using the SyN registration algorithm in Advanced Normalization Tools (https://github.com/ANTsX/ANTs). Next, the parcellation in the atlas space was back-transformed to the individual subject spaces using the inverse transformations from the previous step. The nodes were restricted to the cortical plate, by combining the CRL parcellation with manually delineated tissue segmentation (Xu et al., 2022).
The edge strength between every two nodes i and j was defined by the weighted sum of the streamlines. Specifically, the total number of streamlines between noes i and j was denoted as Nij and Wn was the weight of the nth streamline derived by SIFT2. Since larger ROIs inherently contain more streamlines than smaller ones, which could introduce a bias in favor of nodes with larger volumes exhibiting stronger connectivity, we normalized the number of streamlines by the sum of volumes of the two nodes, which is particularly important in characterizing the developing SCN with growing ROI volumes. Thus, the connectivity between node i and j was expressed as follows:
Network analysis
Network properties
We calculated the global properties of each intrahemispheric network. Global efficiency (Eglob), mean local efficiency (Eloc), shortest path length (Lp), nodal strength, and nodal betweenness were calculated in MATLAB using the GRETNA 2.0 toolbox (https://www.nitrc.org/projects/gretna/). Network efficiency assesses how effectively information flows in the network. Eglob and Eloc focus on the overall network and local groups, respectively. Lp measures the shortest route between nodes in a network. The network modularity evaluates how a network divides into distinct communities, which was calculated using the Louvain community detection algorithm with a resolution coefficient (gamma) of 2, evaluating how the networks divide into distinct communities. We adopted the algorithm proposed by Zhang and Horvath (2005) to calculate the clustering coefficient (CC) generalized for weighted networks, which was defined as follows:
Statistical analysis
Statistical analyses were performed in MATLAB using the GRETNA 2.0 toolbox. To investigate the age dependence of fetal brain SCN, we used Pearson's correlation analysis to evaluate the change of network properties, nodal strength, nodal betweenness, and edge strength with GA. To study the asymmetry of fetal brain SCN, we used paired t tests to compare the network properties, nodal strength, nodal betweenness, and edge strength between the left and right hemispheres of each subject. Subject GA was defined as a covariate when making comparisons. To quantify the degree of laterality patterns, the laterality index (LI) of each measurement was calculated by the following:
Results
Subjects
Out of 171 in utero dMRI scans acquired from 24W to 38W, 13 were excluded due to excessive signal dropouts from intraslice head motion. The remaining 158 scans underwent dMRI preprocessing and reconstruction, FOD estimation, and tractography. Twenty-six scans with signal bias that could not be corrected and 18 with poor fiber tracking results due to noisy reconstructed images were subsequently excluded. Ultimately, a total of 114 scans from 26 to 38 weeks were used for SCN analysis (Fig. 1). The sex of each fetus was unrevealed in compliance with the data privacy policy.⇓⇓
A, Flowchart of data inclusion and exclusion. B, The gestational age distribution of included fetal subjects.
The pipeline for building fetal brain structural connectivity network. Thirty-nine cortical regions in each hemisphere were defined based on the CRL atlas. Whole-brain probabilistic tractography was performed followed by SIFT2 streamline reweighting. The edge strength was calculated by the weighted sum of streamlines between every two nodes normalized by the number of voxels in nodal ROIs.
Change of network properties with age. Eglob, Eloc, and CC increased significantly and Lp decreased significantly over GA in both hemispheres. SW was found above 1 throughout the late-second-to-third gestational trimester and remained stable. The modularity had no significant correlation with GA. Eglob (t = 2.77; FDR-corrected p = 0.013), Eloc (t = 3.41; FDR-corrected p = 0.003), CC (t = 3.65; FDR-corrected p = 0.002), and SW (t = 2.60; FDR-corrected p = 0.016) were significantly lateralized toward the left hemisphere. Lp showed significant rightward lateralization (t = −2.37; FDR-corrected p = 0.023). The network modularity had no significant asymmetry (*p < 0.05, **p < 0.01, ***p < 0.001).
Network properties
Eglob, Eloc, and CC increased significantly while Lp decreased significantly over GA in both hemispheres. SW was found above 1 throughout the late-second-to-third gestational trimester and remained stable. The modularity had no significant correlation with GA (Fig. 3).
Regional connectivity patterns
Nodal strength exhibited widespread significant increase with GA in (FDR-corrected p < 0.05; Fig. 4A), with the correlation coefficient ranging from r = 0.24 (in the left Fusiform) to r = 0.91 (in the left medial orbital part of the superior frontal gyrus, ORBmed). We identified the nodes with r values exceeding the mean plus one standard deviation of all nodal r values as those with the most significant nodal strength increase. Notably, the bilateral ORBmed, anterior cingulate gyrus (ACG), and middle cingulate gyrus (MCG) and the right calcarine and MOG displayed the most significant growth in nodal strength (with r ranging from 0.77 to 0.91). Conversely, the degree of the right parahippocampal gyrus showed a negative correlation with GA (r = −0.31; FDR-corrected p = 0.001).
Regional development of fetal brain SCN. A, Cortical regions with a significant change in nodal strength with GA. B, Cortical regions with a significant change in nodal betweenness with GA. C, Edges with significant increase in strength over GA.
A significant increase in nodal betweenness was mainly observed in the parietal, middle temporal, and cingulate cortices (Fig. 4B), with the correlation coefficient ranging from r = 0.22 (in the right olfactory) to r = 0.47 (in the right anterior cingulum gyrus). The bilateral insula, cuneus, parahippocampal gyrus, and the right rolandic operculum (ROL) exhibited a significant decrease in nodal betweenness over time (with r ranging from −0.24 to −0.51).
Edges were observed to be significantly strengthened over GA, with the most significant increase occurring in the short-ranged connections in the inferior frontal cortex and around the central sulcus (Fig. 4C). A few long-ranged edges were found to be negatively correlated with GA (Fig. 6B), which however, may suffer from errors induced by the overestimation of WM connectivity in early GAs (discussed in Discussions, Challenges in in utero fetal brain SCN analysis).
Lateralization
In terms of global network metrics, Eglob (t = 2.77; FDR-corrected p = 0.013), Eloc (t = 3.41; FDR-corrected p = 0.003), CC, and SW (t = 2.60; FDR-corrected p = 0.016) were significantly lateralized toward the left hemisphere. Lp showed significant rightward lateralization (t = −2.37; FDR-corrected p = 0.023). The network modularity displayed no significant asymmetry. The LI of the global network properties did not exhibit any significant correlation with GA.
Regarding nodal betweenness, the orbital part of the inferior frontal gyrus (ORBinf), the olfactory, the angular gyrus, ROL, and cuneus exhibited significant lateralization to the left (Fig. 5A), with the ORBinf displaying increasing leftward lateralization and the olfactory showing decreasing leftward lateralization with GA (Fig. 5D). In contrast, the rectus, posterior cingulate, inferior and superior occipital, and middle and superior temporal gyri (STG), had significantly higher betweenness in the right hemisphere (Fig. 5A). The LI of the supramarginal gyrus increased significantly with GA (Fig. 5D) but this area had no significant lateralization in general.
Regional lateralization of fetal brain SCN. A–C, Lateralization of nodal betweenness, nodal strength, and edge strength. Inferior frontal and anterior-to-superior temporal nodes exhibited the most significant leftward lateralization while the posterior temporal and the occipital cortices had rightward lateralization. D–F, Correlation between LI of the regional properties (nodal betweenness, nodal strength, and edge strength) and GA. ANG, angular gyrus; CAL, calcarine; CUN, cuneus; HES, Heschl's gyrus; IFGoperc, opercular part of the inferior frontal gyrus; IFGtri, triangular part of the inferior frontal gyrus; IOG, inferior occipital gyrus; ITG, inferior temporal gyrus; MCG, middle cingulate gyrus; MFG, middle frontal gyrus; MOG, middle occipital gyrus; MTG, middle temporal gyrus; OLF, olfactory; ORBinf, orbital part of the inferior frontal gyrus; ORBmed, medial orbital part of the superior frontal gyrus; PCG, posterior cingulate gyrus; PCL, paracentral lobule; PCUN, precuneus; ROL, rolandic operculum; SFG, superior frontal gyrus; SFGmed, medial part of the superior frontal gyrus; SMA, supplementary motor area; SMG, supramarginal gyrus; SOG, superior occipital gyrus; STG, superior temporal gyrus; TPOsup, superior temporal pole.
The nodal strength, which reflects the strength of connectivity for each node in the network, demonstrated leftward lateralization in temporal regions including the superior temporal pole, ROL, inferior temporal, and Heschl's gyri. Frontal regions with leftward asymmetry involved the ORBinf, ORBmed, the triangular part of the inferior frontal gyrus (IFGtri), the middle frontal gyrus, and the supplementary motor area (SMA). The precuneus and cuneus also showed leftward asymmetry (Fig. 5B). The LI of nodal strength had a significant negative correlation with GA in ROL (Fig. 5E), which means it had increasing leftward lateralization. Nodes showing significant rightward lateralization primarily located in the occipital, parietal, and limbic cortices, including the calcarine, cingulate, supramarginal, inferior and superior occipital gyri, and paracentral lobule (Fig. 5B).
The distribution of lateralized network edges, which also rely on the strength of connectivity, showed similar spatial patterns as nodal strength. Edges with leftward asymmetry were observed in all brain lobes, while rightward-lateralized edges were mostly distributed in the posterior brain (Fig. 5C). LI of a few edges showed a significant correlation with GA. LI of the edges between ROL and the frontal cortex, between SMA and the paracentral lobule, as well as between the inferior temporal gyrus and middle occipital gyrus showed a negative correlation with GA. LI of the edge between the inferior occipital and middle temporal gyri, between the supramarginal gyrus and STG, and between the opercular part of the inferior frontal gyrus and the medial part of the superior frontal gyrus had a positive correlation with GA (Fig. 5F).
Discussion
In this study, we collected in utero dMRI to evaluate the normal developmental trajectory of the fetal brain SCN. The observed increases in Eglob, Eloc, and CC, decrease in Lp, and persistent SW from 26W to 38W reflected balanced network integration and segregation, supported by region-specific changes in nodal strength, betweenness, and edge strength. Eglob, Eloc, CC, and SW exhibited significant leftward lateralization, while Lp showed rightward asymmetry. The laterality of nodal strength, betweenness, and edge strength varied across brain regions. Our results provided valuable insights into fetal brain SCN development during the second-to-third trimester.
The normal development of fetal brain SCN
Development of network topology
The increasing Eglob and Eloc with GA indicated enhanced global and local capacities. The decreasing Lp suggested nodes becoming more closely connected. Increase in CC indicated enhancing local clustering. The small-worldness persisted throughout the studied period, reflecting both high clustering and short paths in the network, maintaining high capacity for both local and global information transfer, consistent with previous findings in fetuses (Marami et al., 2017; Song et al., 2017). Our observations validated the balance between network integration and segregation during in utero development (Zhao et al., 2019b), a distinction from postnatal network integration (Huang et al., 2015).
Regional reorganization
The regional analyses revealed general strengthened network connectivity. The growth rate of nodal strength, implied by their correlations with GA, exhibited posterior-to-anterior and medial-to-lateral gradients, which were the most significant in the bilateral ORBmed, ACG, and MCG, along with the right calcarine and SOG. The rapid growth of the calcarine and SOG supported prior findings that the primary visual network strengthened significantly in preterm infants (Cao et al., 2017). The ORBmed is a neonatal structural hub (van den Heuvel and Sporns, 2013) and becomes a functional hub in 2-year-old infants, indicating the emerging default mode network (Gao et al., 2011). The ACG and MCG were reported to be sensitive to infant crying (Mutschler et al., 2016), similar to their social functioning in adults (Fujiwara et al., 2007). ACG also becomes a structural hub in toddlerhood (Huang et al., 2015). These postnatal patterns might be supported by the rapid prenatal structural development we observed.
The strengthening of short-ranged edges coincided with the emergence of corticocortical u-fibers during the second-to-third trimester (Takahashi et al., 2012), possibly contributing to the segregation pattern revealed by increasing Eloc and CC. A few long-ranged edges, negatively correlated with GA, may stem from the overestimations of fiber connectivity in early GAs (Fig. 6B).
Unique challenges in fetal brain SCN analyses using in utero dMRI. A, Colored FA maps and whole-brain tractography in the 26W and 36W fetal brains. The bilateral frontal cortex “touched” with each other on the 26W colored FA and tractography due to insufficient resolution. B, Edges with significant negative correlation with GA. C, The cingulum bundle on the colored FA maps and in tractography. It was blurred with cortical fiber structures on the 26W fetal brain MRI (yellow arrow) and appeared thicker in tractography.
Nodal betweenness notably increased in bilateral frontal, temporal, and cingulate cortices, suggesting their increasingly central roles in the brain, influencing the network's shortest paths. These observations collectively implied intricate reorganization within the fetal brain SCN, characterized by strengthened local connections and a restructured centrality of nodes, maintaining an efficient small-world network organization.
Lateralization
Global lateralization patterns
The Eglob, Eloc, CC, and SW showed significant leftward lateralization and Lp was significantly lower in the left hemisphere, indicating that the left fetal brain hemisphere exhibits closer connections between nodes and within neighborhoods. These observations suggested higher information-transferring capacity in the left hemisphere of the fetal brain, similar to patterns in infants and right-handed adults (Ratnarajah et al., 2013; Caeyenberghs and Leemans, 2014). Our findings suggested that the global laterality patterns are innate and have emerged during in utero neurodevelopment, independent of postnatal experience.
Regional lateralization and possible functional explanations
Regional analysis revealed varied lateralization patternsin the fetal brain. As task-based experiments in utero are challenging, we lack specific knowledge about regional fetal brain functions. Nevertheless, we can still compare the observed structural asymmetries with known functional distributions in later life. Specifically, we compared the observed patterns in the language, motor, and visual-related brain regions with previous studies. Overall, our findings highlighted strong similarities between fetal brain structural asymmetry and the known lateralization patterns in both structural and functional aspects.
Language
Leftward asymmetries in nodal strength or betweenness were observed in language-related regions including the angular, Heschl's, inferior temporal gyri, the superior temporal pole, ROL, and IFGtri, which comprises Broca's area. The betweenness of ORBinf, the degree of ROL, and edges connected with ROL exhibited increasing leftward lateralization. These observations suggested leftward structural lateralization in language-related regions during fetal stage, strengthening with age, similar to findings in neonates (Ratnarajah et al., 2013) and adults (Caeyenberghs and Leemans, 2014). It also coincided with resting-state functional patterns that regions potentially develop into Broca's area became more left-lateralized in older fetuses (Thomason et al., 2014).
Contrarily, the STG, partially overlapping with Wernicke's speech area, displayed rightward-lateralized nodal betweenness. Its leftward functional lateralization was observed in neonates (Perani et al., 2011). However, studies highlighted that cortical folding, associated with local corticocortical fiber connections (Takahashi et al., 2012), progresses earlier and faster in the right STG compared with the left (Dubois et al., 2008; Glasel et al., 2011; Habas et al., 2012; Williams et al., 2021; Xu et al., 2022). This morphological disparity likely induces preferential connectivity underlying the right STG, leading to the observed rightward lateralization.
Motor
Nodal strength displayed leftward asymmetry in SMA, part of the motor cortex (Nachev et al., 2008). The edge between SMA and another motor-related region, the paracentral lobule, showed leftward change with GA (Fig. 4F), possibly leading to the leftward asymmetry of motor-related regions in later development (Ratnarajah et al., 2013; Caeyenberghs and Leemans, 2014).
Visual
Rightward lateralization was observed in primary visual cortices and visuospatial-related regions including the calcarine, cingulate, middle temporal, superior, and inferior occipital gyri. This agreed with the rightward functional asymmetry in the primary and association visual cortices reported in neonates (Williams et al., 2023). Among these regions, the nodal strength of the right calcarine and SOG had faster growth with GA compared with others, resembling the increasing nodal efficiency with age in visuospatial-related regions in the right hemisphere of adults (Caeyenberghs and Leemans, 2014).
To summarize, our observations provided insights into the structural lateralization of the fetal brain and their potential implications for functional laterality in language, motor, and visual-related regions. The specific relationship between structural and functional connectivity needs further investigation.
Challenges in in utero fetal brain SCN analysis
Our study provided preliminary insights into the developmental features of fetal brain SCN using in utero dMRI-based tractography. However, unique features of in utero fetal brain dMRI data present considerable challenges that warrant attention in future investigations. This section summarizes these challenges, particularly in relation to limited image quality and the rapid fetal brain development, and provides suggestions for future investigations.
Immaturity of fetal brain white matter
The fetal brain exhibits tissue heterogeneity and ongoing maturation processes, impacting dMRI properties and complicating tractography. In fetal brains younger than 30W, the highly anisotropic radial glial cells (RGCs) can be estimated by tractography but are not true WM fibers (Kostovic et al., 2002). The RGCs guide neuronal migration and gradually diminish as the process concludes (Sidman and Rakic, 1973; McKinstry et al., 2002). The weakened RGC-related connectivity with GA does not signify a decline in WM connectivity (Fig. 6B), warranting cautious interpretation of fetal brain connectivity from dMRI tractography.
Limited image quality
Acquiring high-quality in utero fetal brain dMRI data is hindered by the unpredictable fetal motion and insufficient resolution relative to the small brain sizes, introducing artifacts and partial volume effects. In our results, for example, the bilateral frontal cortex “touched” with each other (Fig. 6A), and the cingulum bundle merged with adjacent fibers (Fig. 6C). These effects led to overestimations of structural connectivity in early GAs, especially in the interhemispheric regions. To ensure reliability of our findings, we excluded interhemispheric data and performed separate calculations within each hemisphere. This approach allowed investigations on specific scientific questions such as lateralization but overlooked interhemispheric connections. Therefore, future studies aiming for comprehensive fetal brain SCN analyses with higher image quality should integrate interhemispheric connections. Additionally, including subcortical structures can enrich our understanding of the complete fetal brain connectome.
Moreover, employing low imaging resolutions might overlook short-ranged connections, potentially underestimating network measures such as local efficiency. Future studies need sufficient imaging resolution to address this bias.
Handling fetal brains of different GAs
Another inherent challenge in fetal brain imaging studies lies in the considerable difference in brain sizes across different GAs. Employing a uniform imaging resolution, while ensuring standardized comparisons, presents challenges in capturing detailed structures in smaller fetal brains. Using adaptive imaging resolutions for different brains could potentially establish a common “effective resolution,” yet it might yield disparate observations that require careful interpretation. Balancing between keeping methodological consistency and accommodating anatomical variation is crucial.
As for tractography, our study applied a common threshold for all subjects to eliminate erroneous connections, similar to Song et al. (2017). While prevalent studies adopted uniform tracking strategies for fetal brains of different GAs (Hooker et al., 2020; Machado-Rivas et al., 2021), different seeding densities and signal cutoffs across GAs can offer insights into development-related SCN changes. Moreover, the choices of these values lack gold standards and should be rigorously tested when dealing with different datasets.
Conclusion
In this study, we utilized in utero dMRI data collected from 26 to 38 weeks of gestation to investigate the developmental trajectories of the fetal brain structural connectivity network. Our analysis revealed significant increases in both global efficiency and mean local efficiency during the studied period. Small-worldness was found as early as 26 weeks and remained stable until birth. These developmental patterns reflected balanced network integration and segregation in the fetal brain. Short-ranged corticocortical connectivity strengthened significantly, and the development of nodal strength displayed a posterior-to-anterior and medial-to-lateral developmental sequence. The network efficiency and small-worldness showed significant lateralization toward the left hemisphere. The identified lateralization patterns in language, motor, and visual-related regions were consistent with known lateralization patterns in brain structure and function. An exception was observed in Wernicke's area, which displayed a contrasting rightward structural lateralization compared with its known leftward functional asymmetry. Overall, our findings provide novel insights into the development of the fetal brain network and its lateralization.
Footnotes
This work was supported by the Ministry of Science and Technology of the People’s Republic of China (2021ZD0200202), the National Natural Science Foundation of China (81971606, 82122032), the Science and Technology Department of Zhejiang Province (202006140, 2022C03057), and the Young Investigator Program of United Imaging (UIH-QNJJ-2021001).
The authors declare no competing financial interests.
- Correspondence should be addressed to Dan Wu at danwu.bme{at}zju.edu.cn or Guangbin Wang at wgb7932596{at}hotmail.com.












