Faculty Bibliography

Diffusion MRI is being used increasingly in studies of the brain and other parts of the body for its ability to provide quantitative measures that are sensitive to changes in tissue microstructure. However, inter-scanner and inter-protocol differences are known to induce significant measurement variability, which in turn jeopardises the ability to obtain 'truly quantitative measures' and challenges the reliable combination of different datasets. Combining datasets from different scanners and/or acquired at different time points could dramatically increase the statistical power of clinical studies, and facilitate multi-centre research. Even though careful harmonisation of acquisition parameters can reduce variability, inter-protocol differences become almost inevitable with improvements in hardware and sequence design over time, even within a site. In this work, we present a benchmark diffusion MRI database of the same subjects acquired on three distinct scanners with different maximum gradient strength (40, 80, and 300 mT/m), and with 'standard' and 'state-of-the-art' protocols, where the latter have higher spatial and angular resolution. The dataset serves as a useful testbed for method development in cross-scanner/cross-protocol diffusion MRI harmonisation and quality enhancement. Using the database, we compare the performance of five different methods for estimating mappings between the scanners and protocols. The results show that cross-scanner harmonisation of single-shell diffusion data sets can reduce the variability between scanners, and highlight the promises and shortcomings of today's data harmonisation techniques.

Tissue microstructure modeling of diffusion MRI signal is an active research area striving to bridge the gap between macroscopic MRI resolution and cellular-level tissue architecture. Such modeling in neuronal tissue relies on a number of assumptions about the microstructural features of axonal fiber bundles, such as the axonal shape (e.g., perfect cylinders) and the fiber orientation dispersion. However, these assumptions have not yet been validated by sufficiently high-resolution 3-dimensional histology. Here, we reconstructed sequential scanning electron microscopy images in mouse brain corpus callosum, and introduced a random-walker (RaW)-based algorithm to rapidly segment individual intra-axonal spaces and myelin sheaths of myelinated axons. Confirmed by a segmentation based on human annotations initiated with conventional machine-learning-based carving, our semi-automatic algorithm is reliable and less time-consuming. Based on the segmentation, we calculated MRI-relevant estimates of size-related parameters (inner axonal diameter, its distribution, along-axon variation, and myelin g-ratio), and orientation-related parameters (fiber orientation distribution and its rotational invariants; dispersion angle). The reported dispersion angle is consistent with previous 2-dimensional histology studies and diffusion MRI measurements, while the reported diameter exceeds those in other mouse brain studies. Furthermore, we calculated how these quantities would evolve in actual diffusion MRI experiments as a function of diffusion time, thereby providing a coarse-graining window on the microstructure, and showed that the orientation-related metrics have negligible diffusion time-dependence over clinical and pre-clinical diffusion time ranges. However, the MRI-measured inner axonal diameters, dominated by the widest cross sections, effectively decrease with diffusion time by ~ 17% due to the coarse-graining over axonal caliber variations. Furthermore, our 3d measurement showed that there is significant variation of the diameter along the axon. Hence, fiber orientation dispersion estimated from MRI should be relatively stable, while the "apparent" inner axonal diameters are sensitive to experimental settings, and cannot be modeled by perfectly cylindrical axons.

Development of therapies for neurological disorders depends on our ability to non-invasively diagnose and monitor the progression of underlying pathologies at the cellular level. Physics and physiology limit the resolution of human MRI to be orders of magnitude coarser than cell dimensions. Here we identify and quantify the MRI signal coming from within micrometer-thin axons in human white matter tracts in vivo, by utilizing the sensitivity of diffusion MRI to Brownian motion of water molecules restricted by cell walls. We study a specific power-law scaling of the diffusion MRI signal with the diffusion weighting, predicted for water confined to narrow axons, and quantify axonal water fraction and orientation dispersion.

This work evaluates the accuracy and precision of the Diffusion parameter EStImation with Gibbs and NoisE Removal (DESIGNER) pipeline, developed to identify and minimize common sources of methodological variability including: thermal noise, Gibbs ringing artifacts, Rician bias, EPI and eddy current induced spatial distortions, and motion-related artifacts. Following this processing pipeline, iterative parameter estimation techniques were used to derive diffusion parameters of interest based on the diffusion tensor and kurtosis tensor. We evaluated accuracy using a software phantom based on 36 diffusion datasets from the Human Connectome project and tested the precision by analyzing data from 30 healthy volunteers scanned three times within one week. Preprocessing with both DESIGNER or a standard pipeline based on smoothing (instead of noise removal) improved parameter precision by up to a factor of 2 compared to preprocessing with motion correction alone. When evaluating accuracy, we report average decreases in bias (deviation from simulated parameters) over all included regions for fractional anisotropy, mean diffusivity, mean kurtosis, and axonal water fraction of 9.7%, 8.7%, 4.2%, and 7.6% using DESIGNER compared to the standard pipeline, demonstrating that preprocessing with DESIGNER improves accuracy compared to other processing methods.

We measured and compared heritability estimates for measures of functional brain connectivity extracted using the Enhancing Neuroimaging Genetics through Meta-Analysis (ENIGMA) rsfMRI analysis pipeline in two cohorts: the genetics of brain structure (GOBS) cohort and the HCP (the Human Connectome Project) cohort. These two cohorts were assessed using conventional (GOBS) and advanced (HCP) rsfMRI protocols, offering a test case for harmonization of rsfMRI phenotypes, and to determine measures that show consistent heritability for in-depth genome-wide analysis. The GOBS cohort consisted of 334 Mexican-American individuals (124M/210F, average age = 47.9 ± 13.2 years) from 29 extended pedigrees (average family size = 9 people; range 5-32). The GOBS rsfMRI data was collected using a 7.5-min acquisition sequence (spatial resolution = 1.72 × 1.72 × 3 mm³ ). The HCP cohort consisted of 518 twins and family members (240M/278F; average age = 28.7 ± 3.7 years). rsfMRI data was collected using 28.8-min sequence (spatial resolution = 2 × 2 × 2 mm³ ). We used the single-modality ENIGMA rsfMRI preprocessing pipeline to estimate heritability values for measures from eight major functional networks, using (1) seed-based connectivity and (2) dual regression approaches. We observed significant heritability (h² = 0.2-0.4, p < .05) for functional connections from seven networks across both cohorts, with a significant positive correlation between heritability estimates across two cohorts. The similarity in heritability estimates for resting state connectivity measurements suggests that the additive genetic contribution to functional connectivity is robustly detectable across populations and imaging acquisition parameters. The overarching genetic influence, and means to consistently detect it, provides an opportunity to define a common genetic search space for future gene discovery studies.

Biophysical modeling of macroscopic diffusion-weighted MRI (dMRI) signal in terms of microscopic cellular parameters holds the promise of quantifying the integrity of white matter. Unfortunately, even fairly simple multi-compartment models of proton diffusion in the white matter do not provide a unique, biophysically plausible solution. Here we report a nontrivial diffusion MRI signal dependence on echo time (TE) in human white matter in vivo. We demonstrate that such TE dependence originates from compartment-specific T2 values and that it is a promising "orthogonal measure" able to break the degeneracy in parameter estimation, and to yield important relaxation metrics robustly. We thereby enable the precise estimation of the intra- and extra-axonal water T2 relaxation values, which is precluded by a limited signal-to-noise ratio when using multi-echo relaxometry alone.

For many pathologies, early structural tissue changes occur at the cellular level, on the scale of micrometers or tens of micrometers. Magnetic resonance imaging (MRI) is a powerful non-invasive imaging tool used for medical diagnosis, but its clinical hardware is incapable of reaching the cellular length scale directly. In spite of this limitation, microscopic tissue changes in pathology can potentially be captured indirectly, from macroscopic imaging characteristics, by studying water diffusion. Here we focus on water diffusion and NMR relaxation in the human prostate, a highly heterogeneous organ at the cellular level. We present a physical picture of water diffusion and NMR relaxation in the prostate tissue, that is comprised of a densely-packed cellular compartment (composed of stroma and epithelium), and a luminal compartment with almost unrestricted water diffusion. Transverse NMR relaxation is used to identify fast and slow T

PURPOSE:Diffusion kurtosis imaging (DKI) is an advanced magnetic resonance imaging modality that is known to be sensitive to changes in the underlying microstructure of the brain. Image voxels in diffusion weighted images, however, are typically relatively large making them susceptible to partial volume effects, especially when part of the voxel contains cerebrospinal fluid. In this work, we introduce the "Diffusion Kurtosis Imaging with Free Water Elimination" (DKI-FWE) model that separates the signal contributions of free water and tissue, where the latter is modeled using DKI. THEORY AND METHODS:A theoretical study of the DKI-FWE model, including an optimal experiment design and an evaluation of the relative goodness of fit, is carried out. To stabilize the ill-conditioned estimation process, a Bayesian approach with a shrinkage prior (BSP) is proposed. In subsequent steps, the DKI-FWE model and the BSP estimation approach are evaluated in terms of estimation error, both in simulation and real data experiments. RESULTS:Although it is shown that the DKI-FWE model parameter estimation problem is ill-conditioned, DKI-FWE was found to describe the data significantly better compared to the standard DKI model for a large range of free water fractions. The acquisition protocol was optimized in terms of the maximally attainable precision of the DKI-FWE model parameters. The BSP estimator is shown to provide reliable DKI-FWE model parameter estimates. CONCLUSION:The combination of the DKI-FWE model with BSP is shown to be a feasible approach to estimate DKI parameters, while simultaneously eliminating free water partial volume effects. Magn Reson Med 80:802-813, 2018. © 2018 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 NonCommercial License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited and is not used for commercial purposes.

We develop a general analytical and numerical framework for estimating intra- and extra-neurite water fractions and diffusion coefficients, as well as neurite orientational dispersion, in each imaging voxel. By employing a set of rotational invariants and their expansion in the powers of diffusion weighting, we analytically uncover the nontrivial topology of the parameter estimation landscape, showing that multiple branches of parameters describe the measurement almost equally well, with only one of them corresponding to the biophysical reality. A comprehensive acquisition shows that the branch choice varies across the brain. Our framework reveals hidden degeneracies in MRI parameter estimation for neuronal tissue, provides microstructural and orientational maps in the whole brain without constraints or priors, and connects modern biophysical modeling with clinical MRI.