Faculty Bibliography

PURPOSE/OBJECTIVE:) in tumors, both in animals and in humans. METHODS:) by adjusting the diffusion gradient strength. The tDKI data at each diffusion time t were used for a weighted linear least-squares fit method to estimate the diffusion-time-dependent diffusivity, D(t), and diffusional kurtosis, K(t). RESULTS:median and IQR of the two breast cancers were 70 (94) and 106 (92) ms. CONCLUSION/CONCLUSIONS:The results of this proof-of-concept study substantiate the feasibility of using tDKI to measure cellular-interstitial water exchange time without using an exogenous contrast agent.

PURPOSE/OBJECTIVE:To develop and evaluate a neural network-based method for Gibbs artifact and noise removal. METHODS:A convolutional neural network (CNN) was designed for artifact removal in diffusion-weighted imaging data. Two implementations were considered: one for magnitude images and one for complex images. Both models were based on the same encoder-decoder structure and were trained by simulating MRI acquisitions on synthetic non-MRI images. RESULTS:Both machine learning methods were able to mitigate artifacts in diffusion-weighted images and diffusion parameter maps. The CNN for complex images was also able to reduce artifacts in partial Fourier acquisitions. CONCLUSIONS:The proposed CNNs extend the ability of artifact correction in diffusion MRI. The machine learning method described here can be applied on each imaging slice independently, allowing it to be used flexibly in clinical applications.

Background Functional MRI improves preoperative planning in patients with brain tumors, but task-correlated signal intensity changes are only 2%-3% above baseline. This makes accurate functional mapping challenging. Marchenko-Pastur principal component analysis (MP-PCA) provides a novel strategy to separate functional MRI signal from noise without requiring user input or prior data representation. Purpose To determine whether MP-PCA denoising improves activation magnitude for task-based functional MRI language mapping in patients with brain tumors. Materials and Methods In this Health Insurance Portability and Accountability Act-compliant study, MP-PCA performance was first evaluated by using simulated functional MRI data with a known ground truth. Right-handed, left-language-dominant patients with brain tumors who successfully performed verb generation, sentence completion, and finger tapping functional MRI tasks were retrospectively identified between January 2017 and August 2018. On the group level, for each task, histograms of z scores for original and MP-PCA denoised data were extracted from relevant regions and contralateral homologs were seeded by a neuroradiologist blinded to functional MRI findings. Z scores were compared with paired two-sided t tests, and distributions were compared with effect size measurements and the Kolmogorov-Smirnov test. The number of voxels with a z score greater than 3 was used to measure task sensitivity relative to task duration. Results Twenty-three patients (mean age ± standard deviation, 43 years ± 18; 13 women) were evaluated. MP-PCA denoising led to a higher median z score of task-based functional MRI voxel activation in left hemisphere cortical regions for verb generation (from 3.8 ± 1.0 to 4.5 ± 1.4; P < .001), sentence completion (from 3.7 ± 1.0 to 4.3 ± 1.4; P < .001), and finger tapping (from 6.9 ± 2.4 to 7.9 ± 2.9; P < .001). Median z scores did not improve in contralateral homolog regions for verb generation (from -2.7 ± 0.54 to -2.5 ± 0.40; P = .90), sentence completion (from -2.3 ± 0.21 to -2.4 ± 0.37; P = .39), or finger tapping (from -2.3 ± 1.20 to -2.7 ± 1.40; P = .07). Individual functional MRI task durations could be truncated by at least 40% after MP-PCA without degradation of clinically relevant correlations between functional cortex and functional MRI tasks. Conclusion Denoising with Marchenko-Pastur principal component analysis led to higher task correlations in relevant cortical regions during functional MRI language mapping in patients with brain tumors. © RSNA, 2020 Online supplemental material is available for this article.

BACKGROUND:Monte Carlo simulations of diffusion are commonly used as a model validation tool as they are especially suitable for generating the diffusion MRI signal in complicated tissue microgeometries. New method: Here we describe the details of implementing Monte Carlo simulations in three-dimensional (3d) voxelized segmentations of cells in microscopy images. Using the concept of the corner reflector, we largely reduce the computational load of simulating diffusion within and exchange between multiple cells. Precision is further achieved by GPU-based parallel computations. RESULTS:Our simulation result of diffusion in the exemplified white matter axons segmented from a mouse brain demonstrates its power and value in validations of biophysical models. Furthermore, we provide the theoretical background for implementing a discretized diffusion process, and consider the finite-step effects of the particle-membrane reflection and permeation events, needed for efficient simulation of interactions with irregular boundaries, variable diffusion coefficient, and exchange. Comparison with existing methods: To our knowledge, this is the first Monte Carlo pipeline for MR signal simulations in a substrate composed of numerous realistic cells, accounting for their permeable and irregularly-shaped membranes. CONCLUSIONS:effects and flow.

The aspiration of imaging tissue microstructure with MRI is to uncover micrometer-scale tissue features within millimeter-scale imaging voxels, in vivo. This kind of super-resolution has fueled a paradigm shift within the biomedical imaging community. However, what feels like an on-going revolution in MRI, has been conceptually experienced in physics decades ago; from this point of view, our current developments can be seen as Thomas Kuhn's "normal science" stage of progress. While the concept of model-based quantification below the nominal imaging resolution is not new, its possibilities in neuroscience and neuroradiology are only beginning to be widely appreciated. This disconnect calls for communicating the progress of tissue microstructure MR imaging to its potential users. Here, a number of recent research developments are outlined in terms of the overarching concept of coarse-graining the tissue structure over an increasing diffusion length. A variety of diffusion models and phenomena are summarized on the phase diagram of diffusion MRI, with the unresolved problems and future directions corresponding to its unexplored domains.

To study axonal microstructure with diffusion MRI, axons are typically modeled as straight impermeable cylinders, whereby the transverse diffusion MRI signal can be made sensitive to the cylinder's inner diameter. However, the shape of a real axon varies along the axon direction, which couples the longitudinal and transverse diffusion of the overall axon direction. Here we develop a theory of the intra-axonal diffusion MRI signal based on coarse-graining of the axonal shape by 3-dimensional diffusion. We demonstrate how the estimate of the inner diameter is confounded by the diameter variations (beading), and by the local variations in direction (undulations) along the axon. We analytically relate diffusion MRI metrics, such as time-dependent radial diffusivity D_⊥(t) and kurtosis K_⊥(t), to the axonal shape, and validate our theory using Monte Carlo simulations in synthetic undulating axons with randomly positioned beads, and in realistic axons reconstructed from electron microscopy images of mouse brain white matter. We show that (i) In the narrow pulse limit, the inner diameter from D_⊥(t) is overestimated by about twofold due to a combination of axon caliber variations and undulations (each contributing a comparable effect size); (ii) The narrow-pulse kurtosis K_⊥|_{t→∞} deviates from that in an ideal cylinder due to caliber variations; we also numerically calculate the fourth-order cumulant for an ideal cylinder in the wide pulse limit, which is relevant for inner diameter overestimation; (iii) In the wide pulse limit, the axon diameter overestimation is mainly due to undulations at low diffusion weightings b; and (iv) The effect of undulations can be considerably reduced by directional averaging of high-b signals, with the apparent inner diameter given by a combination of the axon caliber (dominated by the thickest axons), caliber variations, and the residual contribution of undulations.

MRI provides a unique non-invasive window into the brain, yet is limited to millimeter resolution, orders of magnitude coarser than cell dimensions. Here, we show that diffusion MRI is sensitive to the micrometer-scale variations in axon caliber or pathological beading, by identifying a signature power-law diffusion time-dependence of the along-fiber diffusion coefficient. We observe this signature in human brain white matter and identify its origins by Monte Carlo simulations in realistic substrates from 3-dimensional electron microscopy of mouse corpus callosum. Simulations reveal that the time-dependence originates from axon caliber variation, rather than from mitochondria or axonal undulations. We report a decreased amplitude of time-dependence in multiple sclerosis lesions, illustrating the potential sensitivity of our method to axonal beading in a plethora of neurodegenerative disorders. This specificity to microstructure offers an exciting possibility of bridging across scales to image cellular-level pathology with a clinically feasible MRI technique.

The dependence of the diffusion MRI signal on the diffusion time t is a hallmark of tissue microstructure at the scale of the diffusion length. Here we measure the time-dependence of the mean diffusivity D(t) and mean kurtosis K(t) in cortical gray matter and in 25 gray matter sub-regions, in 10 healthy subjects. Significant diffusivity and kurtosis time-dependence is observed for t=21.2-100 ms, and is characterized by a power-law tail ∼t^-ϑ with dynamical exponent ϑ. To interpret our measurements, we systematize the relevant scenarios and mechanisms for diffusion time-dependence in the brain. Using effective medium theory formalisms, we derive an exact relation between the power-law tails in D(t) and K(t). The estimated power-law dynamical exponent ϑ≃1/2 in both D(t) and K(t) is consistent with one-dimensional diffusion in the presence of randomly positioned restrictions along neurites. We analyze the short-range disordered statistics of synapses on axon collaterals in the cortex, and perform one-dimensional Monte Carlo simulations of diffusion restricted by permeable barriers with a similar randomness in their placement, to confirm the ϑ=1/2 exponent. In contrast, the Kärger model of exchange is less consistent with the data since it does not capture the diffusivity time-dependence, and the estimated exchange time from K(t) falls below our measured t-range. Although we cannot exclude exchange as a contributing factor, we argue that structural disorder along neurites is mainly responsible for the observed time-dependence of diffusivity and kurtosis. Our observation and theoretical interpretation of the t^-1/2 tail in D(t) and K(t) alltogether establish the sensitivity of a macroscopic MRI signal to micrometer-scale structural heterogeneities along neurites in human gray matter in vivo.

Beta amyloid (Aβ) accumulation is the earliest pathological marker of Alzheimer's disease (AD), but early AD pathology also affects white matter (WM) integrity. We performed a cross-sectional study including 44 subjects (23 healthy controls and 21 mild cognitive impairment or early AD patients) who underwent simultaneous PET-MR using 18F-Florbetapir, and were categorized into 3 groups based on Aβ burden: Aβ- [mean mSUVr ≤1.00], Aβi [1.00 < mSUVr <1.17], Aβ+ [mSUVr ≥1.17]. Intergroup comparisons of diffusion MRI metrics revealed significant differences across multiple WM tracts. Aβi group displayed more restricted diffusion (higher fractional anisotropy, radial kurtosis, axonal water fraction, and lower radial diffusivity) than both Aβ- and Aβ+ groups. This nonmonotonic trend was confirmed by significant continuous correlations between mSUVr and diffusion metrics going in opposite direction for 2 cohorts: pooled Aβ-/Aβi and pooled Aβi/Aβ+. The transient period of increased diffusion restriction may be due to inflammation that accompanies rising Aβ burden. In the later stages of Aβ accumulation, neurodegeneration is the predominant factor affecting diffusion.

Multi-parametric quantitative MRI (qMRI) of the spinal cord is a promising non-invasive tool to probe early microstructural damage in neurological disorders. It is usually performed in vivo by combining acquisitions with multiple signal readouts, which exhibit different thermal noise levels, geometrical distortions and susceptibility to physiological noise. This ultimately hinders joint multi-contrast modelling and makes the geometric correspondence of parametric maps challenging. We propose an approach to overcome these limitations, by implementing state-of-the-art microstructural MRI of the spinal cord with a unified signal readout in vivo (i.e. with matched spatial encoding parameters across a range of imaging contrasts). We base our acquisition on single-shot echo planar imaging with reduced field-of-view, and obtain data from two different vendors (vendor 1: Philips Achieva; vendor 2: Siemens Prisma). Importantly, the unified acquisition allows us to compare signal and noise across contrasts, thus enabling overall quality enhancement via multi-contrast image denoising methods. As a proof-of-concept, here we provide a demonstration with one such method, known as Marchenko-Pastur (MP) Principal Component Analysis (PCA) denoising. MP-PCA is a singular value (SV) decomposition truncation approach that relies on redundant acquisitions, i.e. such that the number of measurements is large compared to the number of components that are maintained in the truncated SV decomposition. Here we used in vivo and synthetic data to test whether a unified readout enables more efficient MP-PCA denoising of less redundant acquisitions, since these can be denoised jointly with more redundant ones. We demonstrate that a unified readout provides robust multi-parametric maps, including diffusion and kurtosis tensors from diffusion MRI, myelin metrics from two-pool magnetisation transfer, and T1 and T2 from relaxometry. Moreover, we show that MP-PCA improves the quality of our multi-contrast acquisitions, since it reduces the coefficient of variation (i.e. variability) by up to 17% for mean kurtosis, 8% for bound pool fraction (myelin-sensitive), and 13% for T1, while enabling more efficient denoising of modalities limited in redundancy (e.g. relaxometry). In conclusion, multi-parametric spinal cord qMRI with unified readout is feasible and provides robust microstructural metrics with matched resolution and distortions, whose quality benefits from multi-contrast denoising methods such as MP-PCA.