Predicting age and clinical risk from the neonatal connectome.

The development of perinatal brain connectivity underpins motor, cognitive and behavioural abilities in later life. Diffusion MRI allows the characterisation of subtle inter-individual differences in structural brain connectivity. Individual brain connectivity maps (connectomes) are by nature high in dimensionality and complex to interpret. Machine learning methods are a powerful tool to uncover properties of the connectome which are not readily visible and can give us clues as to how and why individual developmental trajectories differ. In this manuscript we used Deep Neural Networks and Random Forests to predict demographic and neurodevelopmental characteristics from neonatal structural connectomes in a large sample of babies (n = 524) from the developing Human Connectome Project. We achieved an accurate prediction of post menstrual age (PMA) at scan in term-born infants (mean absolute error (MAE) = 0.72 weeks, r = 0.83 and p < 0.001). We also achieved good accuracy when predicting gestational age at birth in a cohort of term and preterm babies scanned at term equivalent age (MAE = 2.21 weeks, r = 0.82, p < 0.001). We subsequently used sensitivity analysis to obtain feature relevance from our prediction models, with the most important connections for prediction of PMA and GA found to predominantly involve frontal and temporal regions, thalami, and basal ganglia. From our models of PMA at scan for infants born at term, we computed a brain maturation index (predicted age minus actual age) of individual preterm neonates and found a significant correlation between this index and motor outcome at 18 months corrected age. Our results demonstrate the applicability of machine learning techniques in analyses of the neonatal connectome and suggest that a neural substrate of brain maturation with implications for future neurodevelopment is detectable at term equivalent age from the neonatal connectome.

Scattered slice SHARD reconstruction for motion correction in multi-shell diffusion MRI.

Diffusion MRI offers a unique probe into neural microstructure and connectivity in the developing brain. However, analysis of neonatal brain imaging data is complicated by inevitable subject motion, leading to a series of scattered slices that need to be aligned within and across diffusion-weighted contrasts. Here, we develop a reconstruction method for scattered slice multi-shell high angular resolution diffusion imaging (HARDI) data, jointly estimating an uncorrupted data representation and motion parameters at the slice or multiband excitation level. The reconstruction relies on data-driven representation of multi-shell HARDI data using a bespoke spherical harmonics and radial decomposition (SHARD), which avoids imposing model assumptions, thus facilitating to compare various microstructure imaging methods in the reconstructed output. Furthermore, the proposed framework integrates slice-level outlier rejection, distortion correction, and slice profile correction. We evaluate the method in the neonatal cohort of the developing Human Connectome Project (650 scans). Validation experiments demonstrate accurate slice-level motion correction across the age range and across the range of motion in the population. Results in the neonatal data show successful reconstruction even in severely motion-corrupted subjects. In addition, we illustrate how local tissue modelling can extract advanced microstructure features such as orientation distribution functions from the motion-corrected reconstructions.

Preterm birth alters the development of cortical microstructure and morphology at term-equivalent age.

INTRODUCTION: The dynamic nature and complexity of the cellular events that take place during the last trimester of pregnancy make the developing cortex particularly vulnerable to perturbations. Abrupt interruption to normal gestation can lead to significant deviations to many of these processes, resulting in atypical trajectory of cortical maturation in preterm birth survivors. METHODS: We sought to first map typical cortical micro- and macrostructure development using invivo MRI in a large sample of healthy term-born infants scanned after birth (n = 259). Then we offer a comprehensive characterization of the cortical consequences of preterm birth in 76 preterm infants scanned at term-equivalent age (37-44 weeks postmenstrual age). We describe the group-average atypicality, the heterogeneity across individual preterm infants, and relate individual deviations from normative development to age at birth and neurodevelopment at 18 months. RESULTS: In the term-born neonatal brain, we observed heterogeneous and regionally specific associations between age at scan and measures of cortical morphology and microstructure, including rapid surface expansion, greater cortical thickness, lower cortical anisotropy and higher neurite orientation dispersion. By term-equivalent age, preterm infants had on average increased cortical tissue water content and reduced neurite density index in the posterior parts of the cortex, and greater cortical thickness anteriorly compared to term-born infants. While individual preterm infants were more likely to show extreme deviations (over 3.1 standard deviations) from normative cortical maturation compared to term-born infants, these extreme deviations were highly variable and showed very little spatial overlap between individuals. Measures of regional cortical development were associated with age at birth, but not with neurodevelopment at 18 months. CONCLUSION: We showed that preterm birth alters cortical micro- and macrostructural maturation near the time of full-term birth. Deviations from normative development were highly variable between individual preterm infants.

Heterogeneity in Brain Microstructural Development Following Preterm Birth.

Preterm-born children are at increased risk of lifelong neurodevelopmental difficulties. Group-wise analyses of magnetic resonance imaging show many differences between preterm- and term-born infants but do not reliably predict neurocognitive prognosis for individual infants. This might be due to the unrecognized heterogeneity of cerebral injury within the preterm group. This study aimed to determine whether atypical brain microstructural development following preterm birth is significantly variable between infants. Using Gaussian process regression, a technique that allows a single-individual inference, we characterized typical variation of brain microstructure using maps of fractional anisotropy and mean diffusivity in a sample of 270 term-born neonates. Then, we compared 82 preterm infants to these normative values to identify brain regions with atypical microstructure and relate observed deviations to degree of prematurity and neurocognition at 18 months. Preterm infants showed strikingly heterogeneous deviations from typical development, with little spatial overlap between infants. Greater and more extensive deviations, captured by a whole brain atypicality index, were associated with more extreme prematurity and predicted poorer cognitive and language abilities at 18 months. Brain microstructural development after preterm birth is highly variable between individual infants. This poorly understood heterogeneity likely relates to both the etiology and prognosis of brain injury.

On the need for bundle-specific microstructure kernels in diffusion MRI.

Probing microstructure with diffusion magnetic resonance imaging (dMRI) on a scale orders of magnitude below the imaging resolution relies on biophysical modelling of the signal response in the tissue. The vast majority of these biophysical models of diffusion in white matter assume that the measured dMRI signal is the sum of the signals emanating from each of the constituent compartments, each of which exhibits a distinct behaviour in the b-value and/or orientation domain. Many of these models further assume that the dMRI behaviour of the oriented compartments (e.g. the intra-axonal space) is identical between distinct fibre populations, at least at the level of a single voxel. This implicitly assumes that any potential biological differences between fibre populations are negligible, at least as far as is measurable using dMRI. Here, we validate this assumption by means of a voxel-wise, model-free signal decomposition that, under the assumption above and in the absence of noise, is shown to be rank-1. We evaluate the effect size of signal components beyond this rank-1 representation and use permutation testing to assess their significance. We conclude that in the healthy adult brain, the dMRI signal is adequately represented by a rank-1 model, implying that biologically more realistic, but mathematically more complex fascicle-specific microstructure models do not capture statistically significant or anatomically meaningful structure, even in extended high-b diffusion MRI scans.

MRtrix3: A fast, flexible and open software framework for medical image processing and visualisation.

MRtrix3 is an open-source, cross-platform software package for medical image processing, analysis and visualisation, with a particular emphasis on the investigation of the brain using diffusion MRI. It is implemented using a fast, modular and flexible general-purpose code framework for image data access and manipulation, enabling efficient development of new applications, whilst retaining high computational performance and a consistent command-line interface between applications. In this article, we provide a high-level overview of the features of the MRtrix3 framework and general-purpose image processing applications provided with the software.

A framework for multi-component analysis of diffusion MRI data over the neonatal period.

We describe a framework for creating a time-resolved group average template of the developing brain using advanced multi-shell high angular resolution diffusion imaging data, for use in group voxel or fixel-wise analysis, atlas-building, and related applications. This relies on the recently proposed multi-shell multi-tissue constrained spherical deconvolution (MSMT-CSD) technique. We decompose the signal into one isotropic component and two anisotropic components, with response functions estimated from cerebrospinal fluid and white matter in the youngest and oldest participant groups, respectively. We build an orientationally-resolved template of those tissue components from data acquired from 113 babies between 33 and 44 weeks postmenstrual age, imaged as part of the Developing Human Connectome Project. These data were split into weekly groups, and registered to the corresponding group average templates using a previously-proposed non-linear diffeomorphic registration framework, designed to align orientation density functions (ODF). This framework was extended to allow the use of the multiple contrasts provided by the multi-tissue decomposition, and shown to provide superior alignment. Finally, the weekly templates were registered to the same common template to facilitate investigations into the evolution of the different components as a function of age. The resulting multi-tissue atlas provides insights into brain development and accompanying changes in microstructure, and forms the basis for future longitudinal investigations into healthy and pathological white matter maturation.

Modelling white matter with spherical deconvolution: How and why?

Since the realization that diffusion MRI can probe the microstructural organization and orientation of biological tissue in vivo and non-invasively, a multitude of diffusion imaging methods have been developed and applied to study the living human brain. Diffusion tensor imaging was the first model to be widely adopted in clinical and neuroscience research, but it was also clear from the beginning that it suffered from limitations when mapping complex configurations, such as crossing fibres. In this review, we highlight the main steps that have led the field of diffusion imaging to move from the tensor model to the adoption of diffusion and fibre orientation density functions as a more effective way to describe the complexity of white matter organization within each brain voxel. Among several techniques, spherical deconvolution has emerged today as one of the main approaches to model multiple fibre orientations and for tractography applications. Here we illustrate the main concepts and the reasoning behind this technique, as well as the latest developments in the field. The final part of this review provides practical guidelines and recommendations on how to set up processing and acquisition protocols suitable for spherical deconvolution.

Fibre-specific white matter reductions in Alzheimer's disease and mild cognitive impairment.

Alzheimer's disease is increasingly considered a large-scale network disconnection syndrome, associated with progressive aggregation of pathological proteins, cortical atrophy, and functional disconnections between brain regions. These pathological changes are posited to arise in a stereotypical spatiotemporal manner, targeting intrinsic networks in the brain, most notably the default mode network. While this network-specific disruption has been thoroughly studied with functional neuroimaging, changes to specific white matter fibre pathways within the brain's structural networks have not been closely investigated, largely due to the challenges of modelling complex white matter structure. Here, we applied a novel technique known as 'fixel-based analysis' to comprehensively investigate fibre tract-specific differences at a within-voxel level (called 'fixels') to assess potential axonal loss in subjects with Alzheimer's disease and mild cognitive impairment. We hypothesized that patients with Alzheimer's disease would exhibit extensive degeneration across key fibre pathways connecting default network nodes, while patients with mild cognitive impairment would exhibit selective degeneration within fibre pathways connecting regions previously identified as functionally implicated early in Alzheimer's disease. Diffusion MRI data from Alzheimer's disease (n = 49), mild cognitive impairment (n = 33), and healthy elderly control subjects (n = 95) were obtained from the Australian Imaging, Biomarkers and Lifestyle study of ageing. We assessed microstructural differences in fibre density, and macrostructural differences in fibre bundle morphology using fixel-based analysis. Whole-brain analysis was performed to compare groups across all white matter fixels. Subsequently, we performed a tract of interest analysis comparing fibre density and cross-section across 11 selected white matter tracts, to investigate potentially subtle degeneration within fibre pathways in mild cognitive impairment, initially by clinical diagnosis alone, and then by including amyloid status (i.e. a positive or negative amyloid PET scan). Our whole-brain analysis revealed significant white matter loss manifesting both microstructurally and macrostructurally in Alzheimer's disease patients, evident in specific fibre pathways associated with default mode network nodes. Reductions in fibre density and cross-section in mild cognitive impairment patients were only exhibited within the posterior cingulum when statistical analyses were limited to tracts of interest. Interestingly, these degenerative changes did not appear to be associated with high amyloid accumulation, given that amyloid-negative, but not positive, mild cognitive impairment subjects exhibited subtle focal left posterior cingulum deficits. The findings of this study demonstrated a stereotypical distribution of white matter degeneration in patients with Alzheimer's disease, which was in line with canonical findings from other imaging modalities, and with a network-based conceptualization of the disease.awx355media15726254535001.

Time-efficient and flexible design of optimized multishell HARDI diffusion.

PURPOSE: Advanced diffusion magnetic resonance imaging benefits from collecting as much data as is feasible but is highly sensitive to subject motion and the risk of data loss increases with longer acquisition times. Our purpose was to create a maximally time-efficient and flexible diffusion acquisition capability with built-in robustness to partially acquired or interrupted scans. Our framework has been developed for the developing Human Connectome Project, but different application domains are equally possible. METHODS: Complete flexibility in the sampling of diffusion space combined with free choice of phase-encode-direction and the temporal ordering of the sampling scheme was developed taking into account motion robustness, internal consistency, and hardware limits. A split-diffusion-gradient preparation, multiband acceleration, and a restart capacity were added. RESULTS: The framework was used to explore different parameters choices for the desired high angular resolution diffusion imaging diffusion sampling. For the developing Human Connectome Project, a high-angular resolution, maximally time-efficient (20 min) multishell protocol with 300 diffusion-weighted volumes was acquired in >400 neonates. An optimal design of a high-resolution (1.2 × 1.2 mm2 ) two-shell acquisition with 54 diffusion weighted volumes was obtained using a split-gradient design. CONCLUSION: The presented framework provides flexibility to generate time-efficient and motion-robust diffusion magnetic resonance imaging acquisitions taking into account hardware constraints that might otherwise result in sub-optimal choices. Magn Reson Med 79:1276-1292, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Investigating white matter fibre density and morphology using fixel-based analysis.

Voxel-based analysis of diffusion MRI data is increasingly popular. However, most white matter voxels contain contributions from multiple fibre populations (often referred to as crossing fibres), and therefore voxel-averaged quantitative measures (e.g. fractional anisotropy) are not fibre-specific and have poor interpretability. Using higher-order diffusion models, parameters related to fibre density can be extracted for individual fibre populations within each voxel ('fixels'), and recent advances in statistics enable the multi-subject analysis of such data. However, investigating within-voxel microscopic fibre density alone does not account for macroscopic differences in the white matter morphology (e.g. the calibre of a fibre bundle). In this work, we introduce a novel method to investigate the latter, which we call fixel-based morphometry (FBM). To obtain a more complete measure related to the total number of white matter axons, information from both within-voxel microscopic fibre density and macroscopic morphology must be combined. We therefore present the FBM method as an integral piece within a comprehensive fixel-based analysis framework to investigate measures of fibre density, fibre-bundle morphology (cross-section), and a combined measure of fibre density and cross-section. We performed simulations to demonstrate the proposed measures using various transformations of a numerical fibre bundle phantom. Finally, we provide an example of such an analysis by comparing a clinical patient group to a healthy control group, which demonstrates that all three measures provide distinct and complementary information. By capturing information from both sources, the combined fibre density and cross-section measure is likely to be more sensitive to certain pathologies and more directly interpretable.

Early development of structural networks and the impact of prematurity on brain connectivity.

Preterm infants are at high risk of neurodevelopmental impairment, which may be due to altered development of brain connectivity. We aimed to (i) assess structural brain development from 25 to 45 weeks gestational age (GA) using graph theoretical approaches and (ii) test the hypothesis that preterm birth results in altered white matter network topology. Sixty-five infants underwent MRI between 25+3 and 45+6 weeks GA. Structural networks were constructed using constrained spherical deconvolution tractography and were weighted by measures of white matter microstructure (fractional anisotropy, neurite density and orientation dispersion index). We observed regional differences in brain maturation, with connections to and from deep grey matter showing most rapid developmental changes during this period. Intra-frontal, frontal to cingulate, frontal to caudate and inter-hemispheric connections matured more slowly. We demonstrated a core of key connections that was not affected by GA at birth. However, local connectivity involving thalamus, cerebellum, superior frontal lobe, cingulate gyrus and short range cortico-cortical connections was related to the degree of prematurity and contributed to altered global topology of the structural brain network. The relative preservation of core connections at the expense of local connections may support more effective use of impaired white matter reserve following preterm birth.

Language ability in preterm children is associated with arcuate fasciculi microstructure at term.

In the mature human brain, the arcuate fasciculus mediates verbal working memory, word learning, and sublexical speech repetition. However, its contribution to early language acquisition remains unclear. In this work, we aimed to evaluate the role of the direct segments of the arcuate fasciculi in the early acquisition of linguistic function. We imaged a cohort of 43 preterm born infants (median age at birth of 30 gestational weeks; median age at scan of 42 postmenstrual weeks) using high b value high-angular resolution diffusion-weighted neuroimaging and assessed their linguistic performance at 2 years of age. Using constrained spherical deconvolution tractography, we virtually dissected the arcuate fasciculi and measured fractional anisotropy (FA) as a metric of white matter development. We found that term equivalent FA of the left and right arcuate fasciculi was significantly associated with individual differences in linguistic and cognitive abilities in early childhood, independent of the degree of prematurity. These findings suggest that differences in arcuate fasciculi microstructure at the time of normal birth have a significant impact on language development and modulate the first stages of language learning. Hum Brain Mapp 38:3836-3847, 2017. © 2017 Wiley Periodicals, Inc.

Periventricular Nodular Heterotopia: Detection of Abnormal Microanatomic Fiber Structures with Whole-Brain Diffusion MR Imaging Tractography.

Purpose To investigate whether it is possible in patients with periventricular nodular heterotopia (PVNH) to detect abnormal fiber projections that have only previously been reported in the histopathology literature. Materials and Methods Whole-brain diffusion-weighted (DW) imaging data from 14 patients with bilateral PVNH and 14 age- and sex-matched healthy control subjects were prospectively acquired by using 3.0-T magnetic resonance (MR) imaging between August 1, 2008, and December 5, 2012. All participants provided written informed consent. The DW imaging data were processed to generate whole-brain constrained spherical deconvolution (CSD)-based tractography data and super-resolution track-density imaging (TDI) maps. The tractography data were overlaid on coregistered three-dimensional T1-weighted images to visually assess regions of heterotopia. A panel of MR imaging researchers independently assessed each case and indicated numerically (no = 1, yes = 2) as to the presence of abnormal fiber tracks in nodular tissue. The Fleiss κ statistical measure was applied to assess the reader agreement. Results Abnormal fiber tracks emanating from one or more regions of heterotopia were reported by all four readers in all 14 patients with PVNH (Fleiss κ = 1). These abnormal structures were not visible on the tractography data from any of the control subjects and were not discernable on the conventional T1-weighted images of the patients with PVNH. Conclusion Whole-brain CSD-based fiber tractography and super-resolution TDI mapping reveals abnormal fiber projections in nodular tissue suggestive of abnormal organization of white matter (with abnormal fibers both within nodules and projecting to the surrounding white matter) in patients with bilateral PVNH. © RSNA, 2016.

Connectivity-based fixel enhancement: Whole-brain statistical analysis of diffusion MRI measures in the presence of crossing fibres.

In brain regions containing crossing fibre bundles, voxel-average diffusion MRI measures such as fractional anisotropy (FA) are difficult to interpret, and lack within-voxel single fibre population specificity. Recent work has focused on the development of more interpretable quantitative measures that can be associated with a specific fibre population within a voxel containing crossing fibres (herein we use fixel to refer to a specific fibre population within a single voxel). Unfortunately, traditional 3D methods for smoothing and cluster-based statistical inference cannot be used for voxel-based analysis of these measures, since the local neighbourhood for smoothing and cluster formation can be ambiguous when adjacent voxels may have different numbers of fixels, or ill-defined when they belong to different tracts. Here we introduce a novel statistical method to perform whole-brain fixel-based analysis called connectivity-based fixel enhancement (CFE). CFE uses probabilistic tractography to identify structurally connected fixels that are likely to share underlying anatomy and pathology. Probabilistic connectivity information is then used for tract-specific smoothing (prior to the statistical analysis) and enhancement of the statistical map (using a threshold-free cluster enhancement-like approach). To investigate the characteristics of the CFE method, we assessed sensitivity and specificity using a large number of combinations of CFE enhancement parameters and smoothing extents, using simulated pathology generated with a range of test-statistic signal-to-noise ratios in five different white matter regions (chosen to cover a broad range of fibre bundle features). The results suggest that CFE input parameters are relatively insensitive to the characteristics of the simulated pathology. We therefore recommend a single set of CFE parameters that should give near optimal results in future studies where the group effect is unknown. We then demonstrate the proposed method by comparing apparent fibre density between motor neurone disease (MND) patients with control subjects. The MND results illustrate the benefit of fixel-specific statistical inference in white matter regions that contain crossing fibres.

Identification and interpretation of microstructural abnormalities in motor pathways in adolescents born preterm.

There has been extensive interest in assessing the long-term effects of preterm birth on brain white matter microstructure using diffusion MRI. Our aim in this study is to explore diffusion MRI differences between adolescents born preterm and term born controls, with a specific interest in characterising how such differences are manifested in white matter regions containing predominantly single or crossing fibre populations. Probabilistic high angular resolution tractography together with large deformation spatial normalisation were used to objectively investigate diffusion tensor parameters at regular intervals along fibre tracts of 45 adolescents born before 33 weeks of gestation and 30 term-born typically developing adolescents. Diffusion parameters were significantly different between preterms and controls at several levels along the cortico-spinal, thalamo-cortical and transcallosal pathways. Within the predominantly single fibre regions of the corpus callosum and internal capsule, in the preterms mean diffusivity (MD) was found to be increased while fractional anisotropy (FA) was decreased compared to controls. In contrast, however, where these pathways traversed the centrum semiovale, FA and MD were both significantly increased. The major contributor to reduced FA in preterms in predominantly single fibre regions was the increased radial eigenvalue (i.e. increased radial diffusivity). In predominantly crossing-fibre regions, the tensor eigenvalues are not meaningful, and the observed increase in FA is likely to be due to a decrease in anisotropy in one of the contributing fibre bundles. Similar differences (although less pronounced) were observed after excluding preterms with radiological signs of preterm brain injury from the sample. In summary, white matter microstructure was found to be altered in motor pathways in adolescents born preterm. Disruption of white matter (WM) microstructure in a single fibre region with resulting higher radial diffusivity leads to lower FA, whereas selective disruption of one fibre population in a crossing fibre region is observed to lead to higher FA. These findings challenge the common simplistic interpretation of FA as a measure of WM tract integrity.

Quantification of track-weighted imaging (TWI): characterisation of within-subject reproducibility and between-subject variability.

Recently several novel image contrasts derived from whole-brain fibre tracking-data (tractograms) have been introduced. The novel contrasts of these track-weighted imaging (TWI) methods may provide important information for clinical neuroscience studies. However, before they can be used reliably to generate quantitative measures, it is important to characterise their within-subject reproducibility, and between-subject variability. In this work we compute the within-subject reproducibility (intra-scan, intra-session and inter-session), and between-subject variability of TWI for a number of different TWI contrasts across multiple subjects. The results are used in simple voxel-wise power calculations within illustrative regions of interest to provide guidelines for required sample sizes and observable effect sizes for individual subjects and between groups. It was found that the required sample sizes and observable effect sizes varied considerably between different TWI maps and for different ROIs. For some TWI contrast and ROI combinations, the power calculations yielded clinically practical values. These results provide important information concerning the potential usefulness and sensitivity of TWI maps for individual diagnosis, longitudinal studies and group comparisons, as well as for study designs.

Determination of the appropriate b value and number of gradient directions for high-angular-resolution diffusion-weighted imaging.

High-angular-resolution diffusion-weighted imaging (HARDI) is one of the most common MRI acquisition schemes for use with higher order models of diffusion. However, the optimal b value and number of diffusion-weighted (DW) directions for HARDI are still undetermined, primarily as a result of the large number of available reconstruction methods and corresponding parameters, making it impossible to identify a single criterion by which to assess performance. In this study, we estimate the minimum number of DW directions and optimal b values required for HARDI by focusing on the angular frequency content of the DW signal itself. The spherical harmonic (SH) series provides the spherical analogue of the Fourier series, and can hence be used to examine the angular frequency content of the DW signal. Using high-quality data acquired along 500 directions over a range of b values, we estimate that SH terms above l = 8 are negligible in practice for b values up to 5000 s/mm(2), implying that a minimum of 45 DW directions is sufficient to fully characterise the DW signal. l > 0 SH terms were found to increase as a function of b value, levelling off at b = 3000 s/mm(2), suggesting that this value already provides the highest achievable angular resolution. In practice, it is recommended to acquire more than the minimum of 45 DW directions to avoid issues with imperfections in the uniformity of the DW gradient directions and to meet signal-to-noise requirements of the intended reconstruction method.

Reduced structural connectivity in non-motor networks in children born preterm and the influence of early postnatal human cytomegalovirus infection.

Introduction: Preterm birth is increasingly recognized to cause lifelong functional deficits, which often show no correlate in conventional MRI. In addition, early postnatal infection with human cytomegalovirus (hCMV) is being discussed as a possible cause for further impairments. In the present work, we used fixel-based analysis of diffusion-weighted MRI to assess long-term white matter alterations associated with preterm birth and/or early postnatal hCMV infection. Materials and methods: 36 former preterms (PT, median age 14.8 years, median gestational age 28 weeks) and 18 healthy term-born controls (HC, median age 11.1 years) underwent high angular resolution DWI scans (1.5 T, b = 2 000 s/mm2, 60 directions) as well as clinical assessment. All subjects showed normal conventional MRI and normal motor function. Early postnatal hCMV infection status (CMV+ and CMV-) had been determined from repeated screening, ruling out congenital infections. Whole-brain analysis was performed, yielding fixel-wise metrics for fiber density (FD), fiber cross-section (FC), and fiber density and cross-section (FDC). Group differences were identified in a whole-brain analysis, followed by an analysis of tract-averaged metrics within a priori selected tracts associated with cognitive function. Both analyses were repeated while differentiating for postnatal hCMV infection status. Results: PT showed significant reductions of fixel metrics bilaterally in the cingulum, the genu corporis callosum and forceps minor, the capsula externa, and cerebellar and pontine structures. After including intracranial volume as a covariate, reductions remained significant in the cingulum. The tract-specific investigation revealed further reductions bilaterally in the superior longitudinal fasciculus and the uncinate fasciculus. When differentiating for hCMV infection status, no significant differences were found between CMV+ and CMV-. However, comparing CMV+ against HC, fixel metric reductions were of higher magnitude and of larger spatial extent than in CMV- against HC. Conclusion: Preterm birth can lead to long-lasting alterations of WM micro- and macrostructure, not visible on conventional MRI. Alterations are located predominantly in WM structures associated with cognitive function, likely underlying the cognitive deficits observed in our cohort. These observed structural alterations were more pronounced in preterms who suffered from early postnatal hCMV infection, in line with previous studies suggesting an additive effect.