Faculty Bibliography

PURPOSE: To study and reduce the effect of Gibbs ringing artifact on computed diffusion parameters. METHODS: We reduce the ringing by extrapolating the k-space of each diffusion weighted image beyond the measured part by selecting an adequate regularization term. We evaluate several regularization terms and tune the regularization parameter to find the best compromise between anatomical accuracy of the reconstructed image and suppression of the Gibbs artifact. RESULTS: We demonstrate empirically and analytically that the Gibbs artifact, which is typically observed near sharp edges in magnetic resonance images, has a significant impact on the quantification of diffusion model parameters, even for infinitesimal diffusion weighting. We find the second order total generalized variation to be a good choice for the penalty term to regularize the extrapolation of the k-space, as it provides a parsimonious representation of images, a practically full suppression of Gibbs ringing, and the absence of staircasing artifacts typical for total variation methods. CONCLUSIONS: Regularized extrapolation of the k-space data significantly reduces truncation artifacts without compromising spatial resolution in comparison to the default option of window filtering. In particular, accuracy of estimating diffusion tensor imaging and diffusion kurtosis imaging parameters improves so much that unconstrained fits become possible. Magn Reson Med, 2015. (c) 2015 Wiley Periodicals, Inc.

There is a need for accurate quantitative non-invasive biomarkers to monitor myelin pathology in vivo and distinguish myelin changes from other pathological features including inflammation and axonal loss. Conventional MRI metrics such as T2, magnetization transfer ratio and radial diffusivity have proven sensitivity but not specificity. In highly coherent white matter bundles, compartment-specific white matter tract integrity (WMTI) metrics can be directly derived from the diffusion and kurtosis tensors: axonal water fraction, intra-axonal diffusivity, and extra-axonal radial and axial diffusivities. We evaluate the potential of WMTI to quantify demyelination by monitoring the effects of both acute (6weeks) and chronic (12weeks) cuprizone intoxication and subsequent recovery in the mouse corpus callosum, and compare its performance with that of conventional metrics (T2, magnetization transfer, and DTI parameters). The changes observed in vivo correlated with those obtained from quantitative electron microscopy image analysis. A 6-week intoxication produced a significant decrease in axonal water fraction (p<0.001), with only mild changes in extra-axonal radial diffusivity, consistent with patchy demyelination, while a 12-week intoxication caused a more marked decrease in extra-axonal radial diffusivity (p=0.0135), consistent with more severe demyelination and clearance of the extra-axonal space. Results thus revealed increased specificity of the axonal water fraction and extra-axonal radial diffusivity parameters to different degrees and patterns of demyelination. The specificities of these parameters were corroborated by their respective correlations with microstructural features: the axonal water fraction correlated significantly with the electron microscopy derived total axonal water fraction (rho=0.66; p=0.0014) but not with the g-ratio, while the extra-axonal radial diffusivity correlated with the g-ratio (rho=0.48; p=0.0342) but not with the electron microscopy derived axonal water fraction. These parameters represent promising candidates as clinically feasible biomarkers of demyelination and remyelination in the white matter.

The presence of micrometer-level restrictions leads to a decrease of diffusion coefficient with diffusion time. Here we investigate this effect in human white matter in vivo. We focus on a broad range of diffusion times, up to 600 ms, covering diffusion length scales up to about 30 mum. We perform stimulated echo diffusion tensor imaging on 5 healthy volunteers and observe a relatively weak time-dependence in diffusion transverse to major fiber tracts. Remarkably, we also find notable time-dependence in the longitudinal direction. Comparing models of diffusion in ordered, confined and disordered media, we argue that the time-dependence in both directions can arise due to structural disorder, such as axonal beads in the longitudinal direction, and the random packing geometry of fibers within a bundle in the transverse direction. These time-dependent effects extend beyond a simple picture of Gaussian compartments, and may lead to novel markers that are specific to neuronal fiber geometry at the micrometer scale.

The ultimate promise of diffusion MRI (dMRI) models is specificity to neuronal microstructure, which may lead to distinct clinical biomarkers using noninvasive imaging. While multi-compartment models are a common approach to interpret water diffusion in the brain in vivo, the estimation of their parameters from the dMRI signal remains an unresolved problem. Practically, even when q space is highly oversampled, nonlinear fit outputs suffer from heavy bias and poor precision. So far, this has been alleviated by fixing some of the model parameters to a priori values, for improved precision at the expense of accuracy. Here we use a representative two-compartment model to show that fitting fails to determine the five model parameters from over 60 measurement points. For the first time, we identify the reasons for this poor performance. The first reason is the existence of two local minima in the parameter space for the objective function of the fitting procedure. These minima correspond to qualitatively different sets of parameters, yet they both lie within biophysically plausible ranges. We show that, at realistic signal-to-noise ratio values, choosing between the two minima based on the associated objective function values is essentially impossible. Second, there is an ensemble of very low objective function values around each of these minima in the form of a pipe. The existence of such a direction in parameter space, along which the objective function profile is very flat, explains the bias and large uncertainty in parameter estimation, and the spurious parameter correlations: in the presence of noise, the minimum can be randomly displaced by a very large amount along each pipe. Our results suggest that the biophysical interpretation of dMRI model parameters crucially depends on establishing which of the minima is closer to the biophysical reality and the size of the uncertainty associated with each parameter

Dynamic contrast enhanced (DCE) MRI has emerged as a reliable and diagnostically useful functional imaging technique. DCE protocol typically lasts 3-15 minutes and results in a time series of N volumes. For automated analysis, it is important that volumes acquired at different times be spatially coregistered. We have recently introduced a novel 4D, or volume time series, coregistration tool based on a user-specified target volume of interest (VOI). However, the relationship between coregistration accuracy and target VOI size has not been investigated. In this study, coregistration accuracy was quantitatively measured using various sized target VOIs. Coregistration of 10 DCE-MRI mouse head image sets were performed with various sized VOIs targeting the mouse brain. Accuracy was quantified by measures based on the union and standard deviation of the coregistered volume time series. Coregistration accuracy was determined to improve rapidly as the size of the VOI increased and approached the approximate volume of the target (mouse brain). Further inflation of the VOI beyond the volume of the target (mouse brain) only marginally improved coregistration accuracy. The CPU time needed to accomplish coregistration is a linear function of N that varied gradually with VOI size. From the results of this study, we recommend the optimal size of the VOI to be slightly overinclusive, approximately by 5 voxels, of the target for computationally efficient and accurate coregistration.

Diffusion MRI combined with biophysical modeling allows for the description of a white matter (WM) fiber bundle in terms of compartment specific white matter tract integrity (WMTI) metrics, which include intra-axonal diffusivity (Daxon), extra-axonal axial diffusivity (De||), extra-axonal radial diffusivity (De upper left and right quadrants), axonal water fraction (AWF), and tortuosity (alpha) of extra-axonal space. Here we derive these parameters from diffusion kurtosis imaging to examine their relationship to concentrations of global WM N-acetyl-aspartate (NAA), creatine (Cr), choline (Cho) and myo-Inositol (mI), as measured with proton MR spectroscopy (1H-MRS), in a cohort of 25 patients with mild traumatic brain injury (MTBI). We found statistically significant (p<0.05) positive correlations between NAA and Daxon, AWF, alpha, and fractional anisotropy; negative correlations between NAA and De, upper left and right quadrants and the overall radial diffusivity (D upper left and right quadrants). These correlations were supported by similar findings in regional analysis of the genu and splenium of the corpus callosum. Furthermore, a positive correlation in global WM was noted between Daxon and Cr, as well as a positive correlation between De|| and Cho, and a positive trend between De|| and mI. The specific correlations between NAA, an endogenous probe of the neuronal intracellular space, and WMTI metrics related to the intra-axonal space, combined with the specific correlations of De|| with mI and Cho, both predominantly present extra-axonally, corroborate the overarching assumption of many advanced modeling approaches that diffusion imaging can disentangle between the intra- and extra-axonal compartments in WM fiber bundles. Our findings are also generally consistent with what is known about the pathophysiology of MTBI, which appears to involve both intra-axonal injury (as reflected by a positive trend between NAA and Daxon) as well as axonal shrinkage, demyelination, degeneration, and/or loss (as reflected by correlations between NAA and De upper left and right quadrants, AWF, and alpha).

Interpreting brain diffusion MRI measurements in terms of neuronal structure at a micrometer level is an exciting unresolved problem. Here we consider diffusion transverse to a bundle of fibers, and show theoretically, as well as using Monte Carlo simulations and measurements in a phantom made of parallel fibers mimicking axons, that the time dependent diffusion coefficient approaches its macroscopic limit slowly, in a (lnt)/t fashion. The logarithmic singularity arises due to short range disorder in the fiber packing. We identify short range disorder in axonal fibers based on histological data from the splenium, and argue that the time dependent contribution to the overall diffusion coefficient from the extra-axonal water dominates that of the intra-axonal water. This dominance may explain the bias in measuring axon diameters in clinical settings. The short range disorder is also reflected in the linear frequency dependence of the diffusion coefficient measured with oscillating gradients, in agreement with recent experiments. Our results relate the measured diffusion to the mesoscopic structure of neuronal tissue, uncovering the sensitivity of diffusion metrics to axonal arrangement within a fiber tract, and providing an alternative interpretation of axonal diameter mapping techniques.

White matter microstructural changes during the first three years of healthy brain development are characterized using two different models developed for limited clinical diffusion data: White Matter Tract Integrity (WMTI) metrics from Diffusional Kurtosis Imaging (DKI) and Neurite Orientation Dispersion and Density Imaging (NODDI). Both models reveal a non-linear increase in intra-axonal water fraction and in tortuosity of the extra-axonal space as a function of age, in the genu and splenium of the corpus callosum and the posterior limb of the internal capsule. The changes are consistent with expected behavior related to myelination and asynchrony of fiber development. The intra- and extracellular axial diffusivities as estimated with WMTI do not change appreciably in normal brain development. The quantitative differences in parameter estimates between models are examined and explained in the light of each model's assumptions and consequent biases, as highlighted in simulations. Finally, we discuss the feasibility of a model with fewer assumptions.

The purpose of this work was to carry out diffusion tensor imaging (DTI) at multiple diffusion times Td in skeletal muscle in normal subjects and chronic exertional compartment syndrome (CECS) patients and analyze the data with the random permeable barrier model (RPBM) for biophysical specificity. Using an institutional review board approved HIPAA-compliant protocol, seven patients with clinical suspicion of CECS and eight healthy volunteers underwent DTI of the calf muscle in a Siemens MAGNETOM Verio 3 T scanner at rest and after treadmill exertion at four different Td values. Radial diffusion values lambdarad were computed for each of seven different muscle compartments and analyzed with RPBM to produce estimates of free diffusivity D0 , fiber diameter a, and permeability kappa. Fiber diameter estimates were compared with measurements from literature autopsy reference for several compartments. Response factors (post/pre-exercise ratios) were computed and compared between normal controls and CECS patients using a mixed-model two-way analysis of variance. All subjects and muscle compartments showed nearly time-independent diffusion along and strongly time-dependent diffusion transverse to the muscle fibers. RPBM estimates of fiber diameter correlated well with corresponding autopsy reference. D0 showed significant (p < 0.05) increases with exercise for volunteers, and a increased significantly (p < 0.05) in volunteers. At the group level, response factors of all three parameters showed trends differentiating controls from CECS patients, with patients showing smaller diameter changes (p = 0.07), and larger permeability increases (p = 0.07) than controls. Time-dependent diffusion measurements combined with appropriate tissue modeling can provide enhanced microstructural specificity for in vivo tissue characterization. In CECS patients, our results suggest that high-pressure interfiber edema elevates free diffusion and restricts exercise-induced fiber dilation. Such specificity may be useful in differentiating CECS from other disorders or in predicting its response to either physical therapy or fasciotomy