Evren Özarslan

Resolution of complex tissue microarchitecture using the diffusion orientation transform (DOT)

This article describes an accurate and fast method for fiber orientation mapping using multidirectional diffusion-weighted magnetic resonance (MR) data. This novel approach utilizes the Fourier transform relationship between the water displacement probabilities and diffusion-attenuated MR signal expressed in spherical coordinates. The radial part of the Fourier integral is evaluated analytically under the assumption that MR signal attenuates exponentially. The values of the resulting functions are evaluated at a fixed distance away from the origin. The spherical harmonic transform of these functions yields the Laplace series coefficients of the probabilities on a sphere of fixed radius. Alternatively, probability values can be computed nonparametrically using Legendre polynomials. Orientation maps calculated from excised rat nervous tissue data demonstrate this techniques ability to accurately resolve crossing fibers in anatomical regions such as the optic chiasm. This proposed methodology has a trivial extension to multiexponential diffusion-weighted signal decay. The developed methods will improve the reliability of tractography schemes and may make it possible to correctly identify the neural connections between functionally connected regions of the nervous system.

Generalized scalar measures for diffusion MRI using trace, variance, and entropy

This paper details the derivation of rotationally invariant scalar measures from higher‐rank diffusion tensors (DTs) and functions defined on a unit sphere. This was accomplished with the use of an expression that generalizes the evaluation of the trace operator to tensors of arbitrary rank, and even to functions whose domains are the unit sphere. It is shown that the mean diffusivity is invariant to the selection of tensor rank for the model used. However, this rank invariance does not apply to the anisotropy measures. Therefore, a variance‐based, general anisotropy measure is proposed. Also an information theoretical parametrization of anisotropy is introduced that is frequently more consistent with the meaning attributed to anisotropy. We accomplished this by associating anisotropy with the amount of orientational information present in the data, regardless of the imaging technique used. Using a simplified model of fibrous tissue, we simulated anisotropy values with varying orientational complexity and tensor models. Simulations suggested that a lower‐rank tensor model may produce artificially low anisotropy values in voxels with complex structure. This was confirmed with a spin‐echo experiment performed on an excised rat brain. Magn Reson Med 53:866–876, 2005.

Structural insights from high-resolution diffusion tensor imaging and tractography of the isolated rat hippocampus

The hippocampus is a critical structure for learning and memory formation injured by diverse neuropathologies such as epilepsy or Alzheimers disease. Recently, clinical investigations have attempted to use diffusion tensor MRI as a more specific surrogate marker for hippocampal damage. To first better understand the tissue architecture of healthy hippocampal regions, this study characterized 10 rat hippocampi with diffusion tensor imaging (DTI) at 50-microm in-plane image resolution using a 14.1-T magnet. Chemical fixation of the dissected and straightened rat hippocampus provided a simple, effective way to reduce partial volume effects when segmenting hippocampal regions and improved mean signal-to-noise per unit time (e.g. 50.6+/-4.4 at b=1250 s/mm2 in 27 min). Contrary to previous reports that water diffusion is homogeneous throughout the nervous system, statistically different mean diffusivities were observed (e.g. 0.238+/-0.054 and 0.318+/-0.084 microm2/ms for the molecular and granule cell layers respectively) (ANOVA, P<0.05). Different hippocampal subregions had lower fractional anisotropy than uniformly fibrous structures like corpus callosum because of their complex architecture. DTI-derived color fiber orientation maps and tractography demonstrated most components of the trisynaptic intrahippocampal pathway (e.g. orientations in stratum lacunosum-moleculare were dominated by perforant and Schaffer fibers) and also permitted some assessment of connectivity in the rat hippocampus.

Conventions and nomenclature for double diffusion encoding NMR and MRI

Stejskal and Tanners ingenious pulsed field gradient design from 1965 has made diffusion NMR and MRI the mainstay of most studies seeking to resolve microstructural information in porous systems in general and biological systems in particular. Methods extending beyond Stejskal and Tanners design, such as double diffusion encoding (DDE) NMR and MRI, may provide novel quantifiable metrics that are less easily inferred from conventional diffusion acquisitions. Despite the growing interest on the topic, the terminology for the pulse sequences, their parameters, and the metrics that can be derived from them remains inconsistent and disparate among groups active in DDE. Here, we present a consensus of those groups on terminology for DDE sequences and associated concepts. Furthermore, the regimes in which DDE metrics appear to provide microstructural information that cannot be achieved using more conventional counterparts (in a model‐free fashion) are elucidated. We highlight in particular DDEs potential for determining microscopic diffusion anisotropy and microscopic fractional anisotropy, which offer metrics of microscopic features independent of orientation dispersion and thus provide information complementary to the standard, macroscopic, fractional anisotropy conventionally obtained by diffusion MR. Finally, we discuss future vistas and perspectives for DDE. Magn Reson Med 75:82–87, 2016.

Q-space trajectory imaging for multidimensional diffusion MRI of the human brain.

This work describes a new diffusion MR framework for imaging and modeling of microstructure that we call q-space trajectory imaging (QTI). The QTI framework consists of two parts: encoding and modeling. First we propose q-space trajectory encoding, which uses time-varying gradients to probe a trajectory in q-space, in contrast to traditional pulsed field gradient sequences that attempt to probe a point in q-space. Then we propose a microstructure model, the diffusion tensor distribution (DTD) model, which takes advantage of additional information provided by QTI to estimate a distributional model over diffusion tensors. We show that the QTI framework enables microstructure modeling that is not possible with the traditional pulsed gradient encoding as introduced by Stejskal and Tanner. In our analysis of QTI, we find that the well-known scalar b-value naturally extends to a tensor-valued entity, i.e., a diffusion measurement tensor, which we call the b-tensor. We show that b-tensors of rank 2 or 3 enable estimation of the mean and covariance of the DTD model in terms of a second order tensor (the diffusion tensor) and a fourth order tensor. The QTI framework has been designed to improve discrimination of the sizes, shapes, and orientations of diffusion microenvironments within tissue. We derive rotationally invariant scalar quantities describing intuitive microstructural features including size, shape, and orientation coherence measures. To demonstrate the feasibility of QTI on a clinical scanner, we performed a small pilot study comparing a group of five healthy controls with five patients with schizophrenia. The parameter maps derived from QTI were compared between the groups, and 9 out of the 14 parameters investigated showed differences between groups. The ability to measure and model the distribution of diffusion tensors, rather than a quantity that has already been averaged within a voxel, has the potential to provide a powerful paradigm for the study of complex tissue architecture.

Fiber orientation mapping using generalized diffusion tensor imaging

Generalized diffusion tensor imaging uses tensors of arbitrary ranks to model the angular variations in the diffusivities measured by magnetic resonance imaging (MRI) methods. However, a diffusivity profile alone is not readily capable of producing distinct fiber orientations. In this work, we show how it is possible to get the displacement probability profile for water molecules from the higher rank diffusion tensors and validate the technique via simulations of one, two and three fiber systems. Finally, we present fiber orientation results for an image from an excised rat brain.

Clinical feasibility of using mean apparent propagator (MAP) MRI to characterize brain tissue microstructure.

Diffusion tensor imaging (DTI) is the most widely used method for characterizing noninvasively structural and architectural features of brain tissues. However, the assumption of a Gaussian spin displacement distribution intrinsic to DTI weakens its ability to describe intricate tissue microanatomy. Consequently, the biological interpretation of microstructural parameters, such as fractional anisotropy or mean diffusivity, is often equivocal. We evaluate the clinical feasibility of assessing brain tissue microstructure with mean apparent propagator (MAP) MRI, a powerful analytical framework that efficiently measures the probability density function (PDF) of spin displacements and quantifies useful metrics of this PDF indicative of diffusion in complex microstructure (e.g., restrictions, multiple compartments). Rotation invariant and scalar parameters computed from the MAP show consistent variation across neuroanatomical brain regions and increased ability to differentiate tissues with distinct structural and architectural features compared with DTI-derived parameters. The return-to-origin probability (RTOP) appears to reflect cellularity and restrictions better than MD, while the non-Gaussianity (NG) measures diffusion heterogeneity by comprehensively quantifying the deviation between the spin displacement PDF and its Gaussian approximation. Both RTOP and NG can be decomposed in the local anatomical frame for reference determined by the orientation of the diffusion tensor and reveal additional information complementary to DTI. The propagator anisotropy (PA) shows high tissue contrast even in deep brain nuclei and cortical gray matter and is more uniform in white matter than the FA, which drops significantly in regions containing crossing fibers. Orientational profiles of the propagator computed analytically from the MAP MRI series coefficients allow separation of different fiber populations in regions of crossing white matter pathways, which in turn improves our ability to perform whole-brain fiber tractography. Reconstructions from subsampled data sets suggest that MAP MRI parameters can be computed from a relatively small number of DWIs acquired with high b-value and good signal-to-noise ratio in clinically achievable scan durations of less than 10min. The neuroanatomical consistency across healthy subjects and reproducibility in test-retest experiments of MAP MRI microstructural parameters further substantiate the robustness and clinical feasibility of this technique. The MAP MRI metrics could potentially provide more sensitive clinical biomarkers with increased pathophysiological specificity compared to microstructural measures derived using conventional diffusion MRI techniques.

Measurement Tensors in Diffusion MRI: Generalizing the Concept of Diffusion Encoding

In traditional diffusion MRI, short pulsed field gradients (PFG) are used for the diffusion encoding. The standard Stejskal-Tanner sequence uses one single pair of such gradients, known as single-PFG (sPFG). In this work we describe how trajectories in q-space can be used for diffusion encoding. We discuss how such encoding enables the extension of the well-known scalar b-value to a tensor-valued entity we call the diffusion measurement tensor. The new measurements contain information about higher order diffusion propagator covariances not present in sPFG. As an example analysis, we use this new information to estimate a Gaussian distribution over diffusion tensors in each voxel, described by its mean (a diffusion tensor) and its covariance (a 4th order tensor).

Fast orientation mapping from HARDI

This paper introduces a new, accurate and fast method for fiber orientation mapping using high angular resolution diffusion imaging (HARDI) data. The approach utilizes the Fourier relationship between the water displacement probabilities and diffusion attenuated magnetic resonance (MR) signal expressed in spherical coordinates. The Laplace series coefficients of the water displacement probabilities are evaluated at a fixed distance away from the origin. The computations take under one minute for most three-dimensional datasets. We present orientation maps computed from excised rat optic chiasm, brain and spinal cord images. The developed method will improve the reliability of tractography schemes and make it possible to correctly identify the neural connections between functionally connected regions of the nervous system.

Analysis of the effects of noise, DWI sampling, and value of assumed parameters in diffusion MRI models

This study was a systematic evaluation across different and prominent diffusion MRI models to better understand the ways in which scalar metrics are influenced by experimental factors, including experimental design (diffusion‐weighted imaging [DWI] sampling) and noise.