Spatiotemporal tissue maturation of thalamocortical pathways in the human fetal brain
Figures

Methods pipeline to estimate and quantify thalamocortical tracts development.
(Top row) (a) Whole-brain connectomes generated for each gestational week template. (b) Atlas-defined masks of the thalamus and cortical areas were used to extract white matter pathways of interest from the connectomes. (c) These pathways were transformed to the native fetal diffusion space, (d) the values were sampled along the tract. (f) Whole-tract average diffusion metrics were calculated or (g) values sampled along the tract were aligned to an atlas of transient fetal compartments.

Tractography of thalamocortical pathways in different gestational week templates across the second to third trimester.
Tracts project to five different cortical areas, the prefrontal cortex, primary motor cortex, primary sensory cortex, posterior parietal cortex, and primary visual cortex, coloured according to the anterior-posterior axis.

Tract density imaging parcellation of thalamus at different fetal ages.
(a) A schematic of expected cortical connectivity arrangement across the thalamus, based on Morel’s parcellation of the adult thalamic nuclei. (b) Axial slices of thalamic parcellation, thresholded for the top 20% of voxels, colour-coded according to streamline connectivity of different tracts at 23 weeks, (c) 27 weeks, (d) 31 weeks, and (e) 35 weeks.

Diffusion metric age trajectories for each tract.
(a) Distribution of age among the fetal cohort (n=140) in gestational weeks. (b) Whole-tract average tissue (top) and fluid fractions (bottom) for each subject in the left (orange) and right (blue) hemisphere, plotted against gestational age (GA) of the subject, best fit by second-order polynomials (AT = anterior thalamic radiation, OR = optic radiation, PP = posterior parietal tract, TS = thalamic-sensory tract, TM = thalamic-motor tract).

Whole-tract average of diffusion tensor metrics in the fetal cohort (n=140), fractional anisotropy (FA), and mean diffusivity (MD), in thalamocortical tracts across gestational age (every other week shown).

Microstructural composition of fetal compartments traversed by developing thalamic white matter.
Tracts were overlayed on the atlas of fetal compartments (examples highlight the difference between fetal brain structure in early prenatal [25 weeks] on far left, and late prenatal [35 weeks] on far right). Tissue fraction trends (top row) and fluid fraction trends (bottom row), normalised to 1, between the thalamus and cortex (thalamocortical tract axis) for the (a) thalamic-motor tract, (b) optic radiation, and (c) anterior thalamic radiation. Subjects were grouped by age, and average trajectories plotted for early prenatal (22–25.5 weeks), mid prenatal (26–31.5 weeks), late prenatal (32–36 weeks). Error bars represent the standard deviation among all subjects in each group. Atlas-derived tissue boundaries are marked on the trajectories to reveal the changing tissue properties of each layer between early, mid, and late prenatal development (cortical spinal fluid = CSF, cortical plate = CP, subplate = SP, intermediate zone = IZ, ventricular zone = VZ, deep grey matter = GM, immature white matter = WM).

Microstructural composition of fetal compartments traversed by developing thalamic white matter.
(a) Thalamic-sensory tract and (b) posterior parietal tract. Tracts were overlayed on the atlas of fetal compartments (examples highlight the difference between fetal brain structure in early prenatal [25 weeks] on far left, and late prenatal [35 weeks] on far right). Tissue fraction trends (top row) and fluid fraction trends (bottom row) normalised to 1, between the thalamus and cortex (thalamocortical tract axis). Subjects were grouped by age, and average trajectories plotted for early prenatal (22–25.5 weeks), mid prenatal (26–31.5 weeks), late prenatal (32–36 weeks). Error bars represent the standard deviation among all subjects in each group. Atlas-derived tissue boundaries are marked on the trajectories to reveal the changing tissue properties of each layer between early, mid, and late prenatal development (cortical spinal fluid = CSF, cortical plate = CP, subplate = SP, intermediate zone = IZ, ventricular zone = VZ, deep grey matter = GM, immature white matter = WM).

Sociodemographic and neonatal follow-up information.
Including index of multiple deprivation (IMD), London. Gestational age at birth, birth weight, and ethnicity.

Examples of QC scoring scale for subjects, between 0 and 3, with 0 (=failure, e.g. subject moved out of the field of view) to 3 (=high quality).

Correlations between spherical harmonics and radial decomposition (SHARD)-derived indices describing subject motion and gestational age (signal-to-noise ratio [SNR], rotation and translation).

Signal amplitude decay in tissue vs. fluid response functions at different b-shells (left).
Polar plots of the distinct orientation density functions (ODFs) between the corpus callosum (white matter, tissue compartment) and ventricle (cortical spinal fluid [CSF], fluid compartment) after multi-shell multi-tissue constrained spherical deconvolution (MSMT-CSD) is applied to the data (right).

Tissue fraction trajectories along thalamocortical tract axis in each gestational week, prior to grouping into early, mid, and late prenatal (every other week shown).
