Two-Dimensional Correlation of Isotropic and Directional Diffusion Using NMR

Nonparametric distributions of tensor-valued Lorentzian diffusion spectra for model-free data inversion in multidimensional diffusion MRI

Magnetic resonance imaging (MRI) is the method of choice for noninvasive studies of micrometer-scale structures in biological tissues via their effects on the time- and frequency-dependent (restricted) and anisotropic self-diffusion of water. While new designs of time-dependent magnetic field gradient waveforms have enabled disambiguation between different aspects of translational motion that are convolved in traditional MRI methods relying on single pairs of field gradient pulses, data analysis for complex heterogeneous materials remains a challenge. Here, we propose and demonstrate nonparametric distributions of tensor-valued Lorentzian diffusion spectra, or "D(ω) distributions," as a general representation with sufficient flexibility to describe the MRI signal response from a wide range of model systems and biological tissues investigated with modulated gradient waveforms separating and correlating the effects of restricted and anisotropic diffusion.

Probing brain tissue microstructure with MRI: principles, challenges, and the role of multidimensional diffusion-relaxation encoding

Diffusion MRI uses the random displacement of water molecules to sensitize the signal to brain microstructure and to properties such as the density and shape of cells. Microstructure modeling techniques aim to estimate these properties from acquired data by separating the signal between virtual tissue 'compartments' such as the intra-neurite and the extra-cellular space. A key challenge is that the diffusion MRI signal is relatively featureless compared with the complexity of brain tissue. Another challenge is that the tissue microstructure is wildly different within the gray and white matter of the brain. In this review, we use results from multidimensional diffusion encoding techniques to discuss these challenges and their tentative solutions. Multidimensional encoding increases the information content of the data by varying not only the b-value and the encoding direction but also additional experimental parameters such as the shape of the b-tensor and the echo time. Three main insights have emerged from such encoding. First, multidimensional data contradict common model assumptions on diffusion and T2 relaxation, and illustrates how the use of these assumptions cause erroneous interpretations in both healthy brain and pathology. Second, many model assumptions can be dispensed with if data are acquired with multidimensional encoding. The necessary data can be easily acquired in vivo using protocols optimized to minimize Cramér-Rao lower bounds. Third, microscopic diffusion anisotropy reflects the presence of axons but not dendrites. This insight stands in contrast to current 'neurite models' of brain tissue, which assume that axons in white matter and dendrites in gray matter feature highly similar diffusion. Nevertheless, as an axon-based contrast, microscopic anisotropy can differentiate gray and white matter when myelin alterations confound conventional MRI contrasts.

Multidimensional encoding of restricted and anisotropic diffusion by double rotation of the q vector

Diffusion NMR and MRI methods building on the classic pulsed gradient spin-echo sequence are sensitive to many aspects of translational motion, including time and frequency dependence ("restriction"), anisotropy, and flow, leading to ambiguities when interpreting experimental data from complex heterogeneous materials such as living biological tissues. While the oscillating gradient technique specifically targets frequency dependence and permits control of the sensitivity to flow, tensor-valued encoding enables investigations of anisotropy in orientationally disordered materials. Here, we propose a simple scheme derived from the "double-rotation" technique in solid-state NMR to generate a family of modulated gradient waveforms allowing for comprehensive exploration of the 2D frequency-anisotropy space and convenient investigation of both restricted and anisotropic diffusion with a single multidimensional acquisition protocol, thereby combining the desirable characteristics of the oscillating gradient and tensor-valued encoding techniques. The method is demonstrated by measuring multicomponent isotropic Gaussian diffusion in simple liquids, anisotropic Gaussian diffusion in a polydomain lyotropic liquid crystal, and restricted diffusion in a yeast cell sediment.

Mapping the human connectome using diffusion MRI at 300 mT/m gradient strength: Methodological advances and scientific impact

Tremendous efforts have been made in the last decade to advance cutting-edge MRI technology in pursuit of mapping structural connectivity in the living human brain with unprecedented sensitivity and speed. The first Connectom 3T MRI scanner equipped with a 300 mT/m whole-body gradient system was installed at the Massachusetts General Hospital in 2011 and was specifically constructed as part of the Human Connectome Project. Since that time, numerous technological advances have been made to enable the broader use of the Connectom high gradient system for diffusion tractography and tissue microstructure studies and leverage its unique advantages and sensitivity to resolving macroscopic and microscopic structural information in neural tissue for clinical and neuroscientific studies. The goal of this review article is to summarize the technical developments that have emerged in the last decade to support and promote large-scale and scientific studies of the human brain using the Connectom scanner. We provide a brief historical perspective on the development of Connectom gradient technology and the efforts that led to the installation of three other Connectom 3T MRI scanners worldwide - one in the United Kingdom in Cardiff, Wales, another in continental Europe in Leipzig, Germany, and the latest in Asia in Shanghai, China. We summarize the key developments in gradient hardware and image acquisition technology that have formed the backbone of Connectom-related research efforts, including the rich array of high-sensitivity receiver coils, pulse sequences, image artifact correction strategies and data preprocessing methods needed to optimize the quality of high-gradient strength diffusion MRI data for subsequent analyses. Finally, we review the scientific impact of the Connectom MRI scanner, including advances in diffusion tractography, tissue microstructural imaging, ex vivo validation, and clinical investigations that have been enabled by Connectom technology. We conclude with brief insights into the unique value of strong gradients for diffusion MRI and where the field is headed in the coming years.

Multidimensional MRI for Characterization of Subtle Axonal Injury Accelerated Using an Adaptive Nonlocal Multispectral Filter

Multidimensional MRI is an emerging approach that simultaneously encodes water relaxation (T₁ and T₂) and mobility (diffusion) and replaces voxel-averaged values with subvoxel distributions of those MR properties. While conventional (i.e., voxel-averaged) MRI methods cannot adequately quantify the microscopic heterogeneity of biological tissue, using subvoxel information allows to selectively map a specific T₁-T₂-diffusion spectral range that corresponds to a group of tissue elements. The major obstacle to the adoption of rich, multidimensional MRI protocols for diagnostic or monitoring purposes is the prolonged scan time. Our main goal in the present study is to evaluate the performance of a nonlocal estimation of multispectral magnitudes (NESMA) filter on reduced datasets to limit the total acquisition time required for reliable multidimensional MRI characterization of the brain. Here we focused and reprocessed results from a recent study that identified potential imaging biomarkers of axonal injury pathology from the joint analysis of multidimensional MRI, in particular voxelwise T₁-T₂ and diffusion-T₂ spectra in human Corpus Callosum, and histopathological data. We tested the performance of NESMA and its effect on the accuracy of the injury biomarker maps, relative to the co-registered histological reference. Noise reduction improved the accuracy of the resulting injury biomarker maps, while permitting data reduction of 35.7 and 59.6% from the full dataset for T₁-T₂ and diffusion-T₂ cases, respectively. As successful clinical proof-of-concept applications of multidimensional MRI are continuously being introduced, reliable and robust noise removal and consequent acquisition acceleration would advance the field towards clinically-feasible diagnostic multidimensional MRI protocols.

Comparative analysis of signal models for microscopic fractional anisotropy estimation using q-space trajectory encoding

Microscopic diffusion anisotropy imaging using diffusion-weighted MRI and multidimensional diffusion encoding is a promising method for quantifying clinically and scientifically relevant microstructural properties of neural tissue. Several methods for estimating microscopic fractional anisotropy (µFA), a normalized measure of microscopic diffusion anisotropy, have been introduced but the differences between the methods have received little attention thus far. In this study, the accuracy and precision of µFA estimation using q-space trajectory encoding and different signal models were assessed using imaging experiments and simulations. Three healthy volunteers and a microfibre phantom were imaged with five non-zero b-values and gradient waveforms encoding linear and spherical b-tensors. Since the ground-truth µFA was unknown in the imaging experiments, Monte Carlo random walk simulations were performed using axon-mimicking fibres for which the ground truth was known. Furthermore, parameter bias due to time-dependent diffusion was quantified by repeating the simulations with tuned waveforms, which have similar power spectra, and with triple diffusion encoding, which, unlike q-space trajectory encoding, is not based on the assumption of time-independent diffusion. The truncated cumulant expansion of the powder-averaged signal, gamma-distributed diffusivities assumption, and q-space trajectory imaging, a generalization of the truncated cumulant expansion to individual signals, were used to estimate µFA. The gamma-distributed diffusivities assumption consistently resulted in greater µFA values than the second order cumulant expansion, 0.1 greater when averaged over the whole brain. In the simulations, the generalized cumulant expansion provided the most accurate estimates. Importantly, although time-dependent diffusion caused significant overestimation of µFA using all the studied methods, the simulations suggest that the resulting bias in µFA is less than 0.1 in human white matter.

Mapping complex cell morphology in the grey matter with double diffusion encoding MR: A simulation study

This paper investigates the impact of cell body (namely soma) size and branching of cellular projections on diffusion MR imaging (dMRI) and spectroscopy (dMRS) signals for both standard single diffusion encoding (SDE) and more advanced double diffusion encoding (DDE) measurements using numerical simulations. The aim is to investigate the ability of dMRI/dMRS to characterize the complex morphology of brain cells focusing on these two distinctive features of brain grey matter. To this end, we employ a recently developed computational framework to create three dimensional meshes of neuron-like structures for Monte Carlo simulations, using diffusion coefficients typical of water and brain metabolites. Modelling the cellular structure as realistically connected spherical soma and cylindrical cellular projections, we cover a wide range of combinations of sphere radii and branching order of cellular projections, characteristic of various grey matter cells. We assess the impact of spherical soma size and branching order on the b-value dependence of the SDE signal as well as the time dependence of the mean diffusivity (MD) and mean kurtosis (MK). Moreover, we also assess the impact of spherical soma size and branching order on the angular modulation of DDE signal at different mixing times, together with the mixing time dependence of the apparent microscopic anisotropy (μA), a promising contrast derived from DDE measurements. The SDE results show that spherical soma size has a measurable impact on both the b-value dependence of the SDE signal and the MD and MK diffusion time dependence for both water and metabolites. On the other hand, we show that branching order has little impact on either, especially for water. In contrast, the DDE results show that spherical soma size has a measurable impact on the DDE signal's angular modulation at short mixing times and the branching order of cellular projections significantly impacts the mixing time dependence of the DDE signal's angular modulation as well as of the derived μA, for both water and metabolites. Our results confirm that SDE based techniques may be sensitive to spherical soma size, and most importantly, show for the first time that DDE measurements may be more sensitive to the dendritic tree complexity (as parametrized by the branching order of cellular projections), paving the way for new ways of characterizing grey matter morphology, non-invasively using dMRS and potentially dMRI.

Nuclear Magnetic Resonance Microimaging for the Qualitative Assessment of Root Canal Treatment: An Ex Vivo Preliminary Study

AIM

To assess the potential ability of nuclear magnetic resonance micro-imaging (mMRI) to visualize and identify soft tissue debris and unfilled spaces inside radicular canals in endodontic treated extracted teeth, for understanding the causes of treatment failure. Toward this goal, multi-parametric mMRI and cone beam computed tomography (CBCT) were compared.

METHODOLOGY

A non-recoverable root treated human tooth was extracted due to endodontic failure and excessive mobility. It was examined with both CBCT and mMRI: CBCT was performed with 0.125 mm voxel size (GXCB-500, Kavo-Gendex, Brea, CA, USA) and mMRI was performed with a spectrometer operating at 9.4T magnetic field (Bruker Avance-400, Bruker, Billerica, MA, USA). The mMRI images were obtained with a microimaging probe. Relaxation times (T1 and T2) and diffusion-weighted acquisition sequences were used to obtain multi-parametric maps of the extracted tooth (slice thickness of 200 µm and in plane resolution of 30 × 30 µm2).

RESULTS

T1 and T2 maps identified unfilled spaces around and close to Gutta-percha cones instead of CBCT images that were not able to highlight this aspect. T1, T2 and apparent diffusion coefficient (ADC) assumed different values in dentine and in voids, characterized by different dimensions. Moreover, they were able to discriminate between infiltrations of water only and deposits of biological material. Because Gutta-percha cones are constituted of hard, non-porous material, they do not provide a signal and in mMRI images appear as zones of noise.

CONCLUSIONS

Unlike the CBCT exam, mMRI can detect soft tissue debris and unfilled spaces inside radicular canals. Therefore, this in vitro study showed the potential of mMRI to evaluate the quality of the root canal treatment, suggesting its potential benefit in determining the causes of endodontic failure, without the use of ionizing radiation.

Toward nonparametric diffusion- T1 characterization of crossing fibers in the human brain

PURPOSE

To estimate T 1 for each distinct fiber population within voxels containing multiple brain tissue types.

METHODS

A diffusion- T 1 correlation experiment was carried out in an in vivo human brain using tensor-valued diffusion encoding and multiple repetition times. The acquired data were inverted using a Monte Carlo algorithm that retrieves nonparametric distributions P ( D , R 1 ) of diffusion tensors and longitudinal relaxation rates R 1 = 1 / T 1 . Orientation distribution functions (ODFs) of the highly anisotropic components of P ( D , R 1 ) were defined to visualize orientation-specific diffusion-relaxation properties. Finally, Monte Carlo density-peak clustering (MC-DPC) was performed to quantify fiber-specific features and investigate microstructural differences between white matter fiber bundles.

RESULTS

Parameter maps corresponding to P ( D , R 1 ) 's statistical descriptors were obtained, exhibiting the expected R 1 contrast between brain tissue types. Our ODFs recovered local orientations consistent with the known anatomy and indicated differences in R 1 between major crossing fiber bundles. These differences, confirmed by MC-DPC, were in qualitative agreement with previous model-based works but seem biased by the limitations of our current experimental setup.

CONCLUSIONS

Our Monte Carlo framework enables the nonparametric estimation of fiber-specific diffusion- T 1 features, thereby showing potential for characterizing developmental or pathological changes in T 1 within a given fiber bundle, and for investigating interbundle T 1 differences.

Magic DIAMOND: Multi-fascicle diffusion compartment imaging with tensor distribution modeling and tensor-valued diffusion encoding

Diffusion tensor imaging provides increased sensitivity to microstructural tissue changes compared to conventional anatomical imaging but also presents limited specificity. To tackle this problem, the DIAMOND model subdivides the voxel content into diffusion compartments and draws from diffusion-weighted data to estimate compartmental non-central matrix-variate Gamma distributions of diffusion tensors. It models each sub-voxel fascicle separately, resolving crossing white-matter pathways and allowing for a fascicle-element (fixel) based analysis of microstructural features. Alternatively, specific features of the intra-voxel diffusion tensor distribution can be selectively measured using tensor-valued diffusion-weighted acquisition schemes. However, the impact of such schemes on estimating brain microstructural features has only been studied in a handful of parametric single-fascicle models. In this work, we derive a general Laplace transform for the non-central matrix-variate Gamma distribution, which enables the extension of DIAMOND to tensor-valued encoded data. We then evaluate this "Magic DIAMOND" model in silico and in vivo on various combinations of tensor-valued encoded data. Assessing uncertainty on parameter estimation via stratified bootstrap, we investigate both voxel-based and fixel-based metrics by carrying out multi-peak tractography. We demonstrate using in silico evaluations that tensor-valued diffusion encoding significantly improves Magic DIAMOND's accuracy. Most importantly, we show in vivo that our estimated metrics can be robustly mapped along tracks across regions of fiber crossing, which opens new perspectives for tractometry and microstructure mapping along specific white-matter tracts.

A new framework for MR diffusion tensor distribution

The ability to characterize heterogeneous and anisotropic water diffusion processes within macroscopic MRI voxels non-invasively and in vivo is a desideratum in biology, neuroscience, and medicine. While an MRI voxel may contain approximately a microliter of tissue, our goal is to examine intravoxel diffusion processes on the order of picoliters. Here we propose a new theoretical framework and efficient experimental design to describe and measure such intravoxel structural heterogeneity and anisotropy. We assume that a constrained normal tensor-variate distribution (CNTVD) describes the variability of positive definite diffusion tensors within a voxel which extends its applicability to a wide range of b-values while preserving the richness of diffusion tensor distribution (DTD) paradigm unlike existing models. We introduce a new Monte Carlo (MC) scheme to synthesize realistic 6D DTD numerical phantoms and invert the MR signal. We show that the signal inversion is well-posed and estimate the CNTVD parameters parsimoniously by exploiting the different symmetries of the mean and covariance tensors of CNTVD. The robustness of the estimation pipeline is assessed by adding noise to calculated MR signals and compared with the ground truth. A family of invariant parameters and glyphs which characterize microscopic shape, size and orientation heterogeneity within a voxel are also presented.

Computing and visualising intra-voxel orientation-specific relaxation-diffusion features in the human brain

Diffusion MRI techniques are used widely to study the characteristics of the human brain connectome in vivo. However, to resolve and characterise white matter (WM) fibres in heterogeneous MRI voxels remains a challenging problem typically approached with signal models that rely on prior information and constraints. We have recently introduced a 5D relaxation-diffusion correlation framework wherein multidimensional diffusion encoding strategies are used to acquire data at multiple echo-times to increase the amount of information encoded into the signal and ease the constraints needed for signal inversion. Nonparametric Monte Carlo inversion of the resulting datasets yields 5D relaxation-diffusion distributions where contributions from different sub-voxel tissue environments are separated with minimal assumptions on their microscopic properties. Here, we build on the 5D correlation approach to derive fibre-specific metrics that can be mapped throughout the imaged brain volume. Distribution components ascribed to fibrous tissues are resolved, and subsequently mapped to a dense mesh of overlapping orientation bins to define a smooth orientation distribution function (ODF). Moreover, relaxation and diffusion measures are correlated to each independent ODF coordinate, thereby allowing the estimation of orientation-specific relaxation rates and diffusivities. The proposed method is tested on a healthy volunteer, where the estimated ODFs were observed to capture major WM tracts, resolve fibre crossings, and, more importantly, inform on the relaxation and diffusion features along with distinct fibre bundles. If combined with fibre-tracking algorithms, the methodology presented in this work has potential for increasing the depth of characterisation of microstructural properties along individual WM pathways.

Multidimensional correlation MRI

Multidimensional correlation spectroscopy is emerging as a novel MRI modality that is well suited for microstructure and microdynamic imaging studies, especially of biological specimens. Conventional MRI methods only provide voxel-averaged and mostly macroscopically averaged information; these methods cannot disentangle intra-voxel heterogeneity on the basis of both water mobility and local chemical interactions. By correlating multiple MR contrast mechanisms and processing the data in an integrated manner, correlation spectroscopy is able to resolve the distribution of water populations according to their chemical and physical interactions with the environment. The use of a non-parametric, phenomenological representation of the multidimensional MR signal makes no assumptions about tissue structure, thereby allowing the study of microscopic structure and composition of complex heterogeneous biological systems. However, until recently, vast data requirements have confined these types of measurement to non-localized NMR applications and prevented them from being widely and successfully used in conjunction with imaging. Recent groundbreaking advancements have allowed this powerful NMR methodology to be migrated to MRI, initiating its emergence as a promising imaging approach. This review is not intended to cover the entire field of multidimensional MR; instead, it focuses on pioneering imaging applications and the challenges involved. In addition, the background and motivation that have led to multidimensional correlation MR development are discussed, along with the basic underlying mathematical concepts. The goal of the present work is to provide the reader with a fundamental understanding of the techniques developed and their potential benefits, and to provide guidance to help refine future applications of this technology.

Accuracy and precision of statistical descriptors obtained from multidimensional diffusion signal inversion algorithms

In biological tissues, typical MRI voxels comprise multiple microscopic environments, the local organization of which can be captured by microscopic diffusion tensors. The measured diffusion MRI signal can, therefore, be written as the multidimensional Laplace transform of an intravoxel diffusion tensor distribution (DTD). Tensor-valued diffusion encoding schemes have been designed to probe specific features of the DTD, and several algorithms have been introduced to invert such data and estimate statistical descriptors of the DTD, such as the mean diffusivity, the variance of isotropic diffusivities, and the mean squared diffusion anisotropy. However, the accuracy and precision of these estimations have not been assessed systematically and compared across methods. In this article, we perform and compare such estimations in silico for a one-dimensional Gamma fit, a generalized two-term cumulant approach, and two-dimensional and four-dimensional Monte-Carlo-based inversion techniques, using a clinically feasible tensor-valued acquisition scheme. In particular, we compare their performance at different signal-to-noise ratios (SNRs) for voxel contents varying in terms of the aforementioned statistical descriptors, orientational order, and fractions of isotropic and anisotropic components. We find that all inversion techniques share similar precision (except for a lower precision of the two-dimensional Monte Carlo inversion) but differ in terms of accuracy. While the Gamma fit exhibits infinite-SNR biases when the signal deviates strongly from monoexponentiality and is unaffected by orientational order, the generalized cumulant approach shows infinite-SNR biases when this deviation originates from the variance in isotropic diffusivities or from the low orientational order of anisotropic diffusion components. The two-dimensional Monte Carlo inversion shows remarkable accuracy in all systems studied, given that the acquisition scheme possesses enough directions to yield a rotationally invariant powder average. The four-dimensional Monte Carlo inversion presents no infinite-SNR bias, but suffers significantly from noise in the data, while preserving good contrast in most systems investigated.

Diffusion tensor distribution imaging of an in vivo mouse brain at ultrahigh magnetic field by spatiotemporal encoding

Diffusion tensor distribution (DTD) imaging builds on principles from diffusion, solid-state and low-field NMR spectroscopies, to quantify the contents of heterogeneous voxels as nonparametric distributions, with tensor "size", "shape" and orientation having direct relations to corresponding microstructural properties of biological tissues. The approach requires the acquisition of multiple images as a function of the magnitude, shape and direction of the diffusion-encoding gradients, leading to long acquisition times unless fast image read-out techniques like EPI are employed. While in previous in vivo human brain studies performed at 3 T this proved a viable option, porting these measurements to very high magnetic fields and/or to heterogeneous organs induces B0 - and B1 -inhomogeneity artifacts that challenge the limits of EPI. To overcome such challenges, we demonstrate here that high spatial resolution DTD of mouse brain can be carried out at 15.2 T with a surface-cryoprobe, by relying on SPatiotemporal ENcoding (SPEN) imaging sequences. These new acquisition and data-processing protocols are demonstrated with measurements on in vivo mouse brain, and validated with synthetic phantoms designed to mimic the diffusion properties of white matter, gray matter and cerebrospinal fluid. While still in need of full extensions to 3D mappings and of scanning additional animals to extract more general physiological conclusions, this work represents another step towards the model-free, noninvasive in vivo characterization of tissue microstructure and heterogeneity in animal models, at ≈0.1 mm resolutions.

Double diffusion encoding and applications for biomedical imaging

Diffusion Magnetic Resonance Imaging (dMRI) is one of the most important contemporary non-invasive modalities for probing tissue structure at the microscopic scale. The majority of dMRI techniques employ standard single diffusion encoding (SDE) measurements, covering different sequence parameter ranges depending on the complexity of the method. Although many signal representations and biophysical models have been proposed for SDE data, they are intrinsically limited by a lack of specificity. Advanced dMRI methods have been proposed to provide additional microstructural information beyond what can be inferred from SDE. These enhanced contrasts can play important roles in characterizing biological tissues, for instance upon diseases (e.g. neurodegenerative, cancer, stroke), aging, learning, and development. In this review we focus on double diffusion encoding (DDE), which stands out among other advanced acquisitions for its versatility, ability to probe more specific diffusion correlations, and feasibility for preclinical and clinical applications. Various DDE methodologies have been employed to probe compartment sizes (Section 3), decouple the effects of microscopic diffusion anisotropy from orientation dispersion (Section 4), probe displacement correlations, study exchange, or suppress fast diffusing compartments (Section 6). DDE measurements can also be used to improve the robustness of biophysical models (Section 5) and study intra-cellular diffusion via magnetic resonance spectroscopy of metabolites (Section 7). This review discusses all these topics as well as important practical aspects related to the implementation and contrast in preclinical and clinical settings (Section 9) and aims to provide the readers a guide for deciding on the right DDE acquisition for their specific application.

Direct and specific assessment of axonal injury and spinal cord microenvironments using diffusion correlation imaging

We describe a practical two-dimensional (2D) diffusion MRI framework to deliver specificity and improve sensitivity to axonal injury in the spinal cord. This approach provides intravoxel distributions of correlations of water mobilities in orthogonal directions, revealing sub-voxel diffusion components. Here we use it to investigate water diffusivities along axial and radial orientations within spinal cord specimens with confirmed, tract-specific axonal injury. First, we show using transmission electron microscopy and immunohistochemistry that tract-specific axonal beading occurs following Wallerian degeneration in the cortico-spinal tract as direct sequelae to closed head injury. We demonstrate that although some voxel-averaged diffusion tensor imaging (DTI) metrics are sensitive to this axonal injury, they are non-specific, i.e., they do not reveal an underlying biophysical mechanism of injury. Then we employ 2D diffusion correlation imaging (DCI) to improve discrimination of different water microenvironments by measuring and mapping the joint water mobility distributions perpendicular and parallel to the spinal cord axis. We determine six distinct diffusion spectral components that differ according to their microscopic anisotropy and mobility. We show that at the injury site a highly anisotropic diffusion component completely disappears and instead becomes more isotropic. Based on these findings, an injury-specific MR image of the spinal cord was generated, and a radiological-pathological correlation with histological silver staining % area was performed. The resulting strong and significant correlation (r=0.70,p < 0.0001) indicates the high specificity with which DCI detects injury-induced tissue alterations. We predict that the ability to selectively image microstructural changes following axonal injury in the spinal cord can be useful in clinical and research applications by enabling specific detection and increased sensitivity to injury-induced microstructural alterations. These results also encourage us to translate DCI to higher spatial dimensions to enable assessment of traumatic axonal injury, and possibly other diseases and disorders in the brain.

The present and the future of microstructure MRI: From a paradigm shift to normal science

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 ongoing 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.

Transferring principles of solid-state and Laplace NMR to the field of in vivo brain MRI

Magnetic resonance imaging (MRI) is the primary method for noninvasive investigations of the human brain in health, disease, and development but yields data that are difficult to interpret whenever the millimeter-scale voxels contain multiple microscopic tissue environments with different chemical and structural properties. We propose a novel MRI framework to quantify the microscopic heterogeneity of the living human brain as spatially resolved five-dimensional relaxation-diffusion distributions by augmenting a conventional diffusion-weighted imaging sequence with signal encoding principles from multidimensional solid-state nuclear magnetic resonance (NMR) spectroscopy, relaxation-diffusion correlation methods from Laplace NMR of porous media, and Monte Carlo data inversion. The high dimensionality of the distribution space allows resolution of multiple microscopic environments within each heterogeneous voxel as well as their individual characterization with novel statistical measures that combine the chemical sensitivity of the relaxation rates with the link between microstructure and the anisotropic diffusivity of tissue water. The proposed framework is demonstrated on a healthy volunteer using both exhaustive and clinically viable acquisition protocols.

Retaining information from multidimensional correlation MRI using a spectral regions of interest generator

Multidimensional correlation magnetic resonance imaging (MRI) is an emerging imaging modality that is capable of disentangling highly heterogeneous and opaque systems according to chemical and physical interactions of water within them. Using this approach, the conventional three dimensional MR scalar images are replaced with spatially resolved multidimensional spectra. The ensuing abundance in microstructural and chemical information is a blessing that incorporates a real challenge: how does one distill and refine it into images while retaining its significant components? In this paper we introduce a general framework that preserves the spectral information from spatially resolved multidimensional data. Equal weight is given to significant spectral components at the single voxel level, resulting in a summarized image spectrum. This spectrum is then used to define spectral regions of interest that are utilized to reconstruct images of sub-voxel components. Using numerical simulations we first show that, contrary to the conventional approach, the proposed framework preserves spectral resolution, and in turn, sensitivity and specificity of the reconstructed images. The retained spectral resolution allows, for the first time, to observe an array of distinct [Formula: see text]-[Formula: see text]-[Formula: see text] components images of the human brain. The robustly generated images of sub-voxel components overcome the limited spatial resolution of MRI, thus advancing multidimensional correlation MRI to fulfilling its full potential.

Direction-averaged diffusion-weighted MRI signal using different axisymmetric B-tensor encoding schemes

Purpose

It has been shown, theoretically and in vivo, that using the Stejskal‐Tanner pulsed‐gradient, or linear tensor encoding (LTE), and in tissue exhibiting a “stick‐like” diffusion geometry, the direction‐averaged diffusion‐weighted MRI signal at high b‐values ( 7000<b<10000s/mm2) follows a power‐law, decaying as 1/b. It has also been shown, theoretically, that for planar tensor encoding (PTE), the direction‐averaged diffusion‐weighted MRI signal decays as 1/b. We aimed to confirm this theoretical prediction in vivo. We then considered the direction‐averaged signal for arbitrary b‐tensor shapes and different tissue substrates to look for other conditions under which a power‐law exists.

Methods

We considered the signal decay for high b‐values for encoding geometries ranging from 2‐dimensional PTE, through isotropic or spherical tensor encoding to LTE. When a power‐law behavior was suggested, this was tested using in silico simulations and, when appropriate, in vivo using ultra‐strong (300 mT/m) gradients.

Results

Our in vivo results confirmed the predicted 1/b power law for PTE. Moreover, our analysis showed that using an axisymmetric b‐tensor a power‐law only exists under very specific conditions: (a) “stick‐like” tissue geometry and purely LTE or purely PTE waveforms; and (b) "pancake‐like" tissue geometry and a purely LTE waveform.

Conclusions

A complete analysis of the power‐law dependencies of the diffusion‐weighted signal at high b‐values has been performed. Only three specific forms of encoding result in a power‐law dependency, pure linear and pure PTE when the tissue geometry is “stick‐like” and pure LTE when the tissue geometry is "pancake‐like". The different exponents of these encodings could be used to provide independent validation of the presence of different tissue geometries in vivo.

Probing Surface-to-Volume Ratio of an Anisotropic Medium by Diffusion NMR with General Gradient Encoding

Since the seminal paper by Mitra et al., diffusion MR has been widely used in order to estimate surface-to-volume ratios. In this paper, we generalize Mitra's formula for arbitrary diffusion encoding waveforms, including recently developed q-space trajectory encoding sequences. We show that the surface-to-volume ratio can be significantly misestimated using the original Mitra's formula without taking into account the applied gradient profile. In order to obtain more accurate estimation in anisotropic samples, we propose an efficient and robust optimization algorithm to design diffusion gradient waveforms with prescribed features. Our results are supported by Monte Carlo simulations.

Tensor-valued diffusion MRI in under 3 minutes: an initial survey of microscopic anisotropy and tissue heterogeneity in intracranial tumors

PURPOSE

To evaluate the feasibility of a 3-minutes protocol for assessment of the microscopic anisotropy and tissue heterogeneity based on tensor-valued diffusion MRI in a wide range of intracranial tumors.

METHODS

B-tensor encoding was performed in 42 patients with intracranial tumors (gliomas, meningiomas, adenomas, and metastases). Microscopic anisotropy and tissue heterogeneity were evaluated by estimating the anisotropic kurtosis (MK_A ) and isotropic kurtosis (MK_I ), respectively. An extensive imaging protocol was compared with a 3-minutes protocol.

RESULTS

The fast imaging protocol yielded parameters with characteristics in terms of bias and precision similar to the full protocol. Glioblastomas had lower microscopic anisotropy than meningiomas (MK_A = 0.29 ± 0.06 vs. 0.45 ± 0.08, P = 0.003). Metastases had higher tissue heterogeneity (MK_I = 0.57 ± 0.07) than both the glioblastomas (0.44 ± 0.06, P < 0.001) and meningiomas (0.46 ± 0.06, P = 0.03).

CONCLUSION

Evaluation of the microscopic anisotropy and tissue heterogeneity in intracranial tumor patients is feasible in clinically relevant times frames.

Multiple dimensions for random walks

Current trends in diffusion NMR and MRI methods development are reviewed. While great efforts are still directed towards further improving the spectral, spatial, and relaxation rate resolution of basic diffusion measurements, recent improvements in magnetic field gradient technology on whole-body scanners have enabled an exciting line of research involving MRI implementations of advanced diffusion NMR methods with motion-encoding gradient waveforms designed for multidimensional separation and correlation of properties like short-time diffusivity, restriction, anisotropy, flow, and exchange, thereby opening up for highly specific characterization of microstructure and heterogeneity in healthy and diseased tissues in a clinical setting.

In vivo magnetic resonance imaging and spectroscopy. Technological advances and opportunities for applications continue to abound

Over the past decades, the field of in vivo magnetic resonance (MR) has built up an impressive repertoire of data acquisition and analysis technologies for anatomical, functional, physiological, and molecular imaging, the description of which requires many book volumes. As such it is impossible for a few authors to have an authoritative overview of the field and for a brief article to be inclusive. We will therefore focus mainly on data acquisition and attempt to give some insight into the principles underlying current advanced methods in the field and the potential for further innovation. In our view, the foreseeable future is expected to show continued rapid progress, for instance in imaging of microscopic tissue properties in vivo, assessment of functional and anatomical connectivity, higher resolution physiologic and metabolic imaging, and even imaging of receptor binding. In addition, acquisition speed and information content will continue to increase due to the continuous development of approaches for parallel imaging (including simultaneous multi-slice imaging), compressed sensing, and MRI fingerprinting. Finally, artificial intelligence approaches are becoming more realistic and will have a tremendous effect on both acquisition and analysis strategies. Together, these developments will continue to provide opportunity for scientific discovery and, in combination with large data sets from other fields such as genomics, allow the ultimate realization of precision medicine in the clinic.

Multidimensional diffusion MRI with spectrally modulated gradients reveals unprecedented microstructural detail

Characterization of porous media is essential in a wide range of biomedical and industrial applications. Microstructural features can be probed non-invasively by diffusion magnetic resonance imaging (dMRI). However, diffusion encoding in conventional dMRI may yield similar signatures for very different microstructures, which represents a significant limitation for disentangling individual microstructural features in heterogeneous materials. To solve this problem, we propose an augmented multidimensional diffusion encoding (MDE) framework, which unlocks a novel encoding dimension to assess time-dependent diffusion specific to structures with different microscopic anisotropies. Our approach relies on spectral analysis of complex but experimentally efficient MDE waveforms. Two independent contrasts to differentiate features such as cell shape and size can be generated directly by signal subtraction from only three types of measurements. Analytical calculations and simulations support our experimental observations. Proof-of-concept experiments were applied on samples with known and distinctly different microstructures. We further demonstrate substantially different contrasts in different tissue types of a post mortem brain. Our simultaneous assessment of restriction size and shape may be instrumental in studies of a wide range of porous materials, enable new insights into the microstructure of biological tissues or be of great value in diagnostics.

Diffusion tensor distribution imaging

Conventional diffusion MRI yields voxel-averaged parameters that suffer from ambiguities for heterogeneous anisotropic materials such as brain tissue. Using principles from solid-state NMR spectroscopy, we have previously introduced the shape of the diffusion encoding tensor as a separate acquisition dimension that disentangles isotropic and anisotropic contributions to the observed diffusivities, thereby allowing for unconstrained data inversion into diffusion tensor distributions with "size," "shape," and orientation dimensions. Here we combine our recent non-parametric data inversion algorithm and data acquisition protocol with an imaging pulse sequence to demonstrate spatial mapping of diffusion tensor distributions using a previously developed composite phantom with multiple isotropic and anisotropic components. We propose a compact format for visualizing two-dimensional arrays of the distributions, new scalar parameters quantifying intra-voxel heterogeneity, and a binning procedure giving maps of all relevant parameters for each of the components resolved in the multidimensional distribution space.

Multimodal integration of diffusion MRI for better characterization of tissue biology

The contrast in diffusion-weighted MR images is due to variations of diffusion properties within the examined specimen. Certain microstructural information on the underlying tissues can be inferred through quantitative analyses of the diffusion-sensitized MR signals. In the first part of the paper, we review two types of approach for characterizing diffusion MRI signals: Bloch's equations with diffusion terms, and statistical descriptions. Specifically, we discuss expansions in terms of cumulants and orthogonal basis functions, the confinement tensor formalism and tensor distribution models. Further insights into the tissue properties may be obtained by integrating diffusion MRI with other techniques, which is the subject of the second part of the paper. We review examples involving magnetic susceptibility, structural tensors, internal field gradients, transverse relaxation and functional MRI. Integrating information provided by other imaging modalities (MR based or otherwise) could be a key to improve our understanding of how diffusion MRI relates to physiology and biology.

Imaging brain microstructure with diffusion MRI: practicality and applications

This article gives an overview of microstructure imaging of the brain with diffusion MRI and reviews the state of the art. The microstructure-imaging paradigm aims to estimate and map microscopic properties of tissue using a model that links these properties to the voxel scale MR signal. Imaging techniques of this type are just starting to make the transition from the technical research domain to wide application in biomedical studies. We focus here on the practicalities of both implementing such techniques and using them in applications. Specifically, the article summarizes the relevant aspects of brain microanatomy and the range of diffusion-weighted MR measurements that provide sensitivity to them. It then reviews the evolution of mathematical and computational models that relate the diffusion MR signal to brain tissue microstructure, as well as the expanding areas of application. Next we focus on practicalities of designing a working microstructure imaging technique: model selection, experiment design, parameter estimation, validation, and the pipeline of development of this class of technique. The article concludes with some future perspectives on opportunities in this topic and expectations on how the field will evolve in the short-to-medium term.

Quantifying brain microstructure with diffusion MRI: Theory and parameter estimation

We review, systematize and discuss models of diffusion in neuronal tissue, by putting them into an overarching physical context of coarse-graining over an increasing diffusion length scale. From this perspective, we view research on quantifying brain microstructure as occurring along three major avenues. The first avenue focusses on transient, or time-dependent, effects in diffusion. These effects signify the gradual coarse-graining of tissue structure, which occurs qualitatively differently in different brain tissue compartments. We show that transient effects contain information about the relevant length scales for neuronal tissue, such as the packing correlation length for neuronal fibers, as well as the degree of structural disorder along the neurites. The second avenue corresponds to the long-time limit, when the observed signal can be approximated as a sum of multiple nonexchanging anisotropic Gaussian components. Here, the challenge lies in parameter estimation and in resolving its hidden degeneracies. The third avenue employs multiple diffusion encoding techniques, able to access information not contained in the conventional diffusion propagator. We conclude with our outlook on future directions that could open exciting possibilities for designing quantitative markers of tissue physiology and pathology, based on methods of studying mesoscopic transport in disordered systems.

A Feature Extraction Method for P2P Botnet Detection Using Graphic Symmetry Concept

A DDoS (Distributed Denial of Service) attack makes use of a botnet to launch attacks and cause node congestion of wireless sensor networks, which is a common and serious threat. Due to the various kinds of features required in a Peer-to-Peer (P2P) botnet for DDoS attack detection via current machine learning methods and the failure to effectively detect encrypted botnets, this paper extracts the data packet size and the symmetric intervals in flow according to the concept of graphic symmetry. Combined with flow information entropy and session features, the frequency domain features can be sorted so as to obtain features with better correlations, which solves the problem of multiple types of features required for detection. Information entropy corresponding to the flow size can distinguish an encrypted botnet. This method is implemented through machine learning techniques. Experimental results show that the proposed method can detect the P2P botnet for DDoS attack and the detection accuracy is higher than that of traditional feature detection.

Microscopic anisotropy misestimation in spherical-mean single diffusion encoding MRI

Purpose

Microscopic fractional anisotropy (µFA) can disentangle microstructural information from orientation dispersion. While double diffusion encoding (DDE) MRI methods are widely used to extract accurate µFA, it has only recently been proposed that powder‐averaged single diffusion encoding (SDE) signals, when coupled with the diffusion standard model (SM) and a set of constraints, could be used for µFA estimation. This study aims to evaluate µFA as derived from the spherical mean technique (SMT) set of constraints, as well as more generally for powder‐averaged SM signals.

Methods

SDE experiments were performed at 16.4 T on an ex vivo mouse brain (Δ/δ = 12/1.5 ms). The µFA maps obtained from powder‐averaged SDE signals were then compared to maps obtained from DDE‐MRI experiments (Δ/τ/δ = 12/12/1.5 ms), which allow a model‐free estimation of µFA. Theory and simulations that consider different types of heterogeneity are presented for corroborating the experimental findings.

Results

µFA, as well as other estimates derived from powder‐averaged SDE signals produced large deviations from the ground truth in both gray and white matter. Simulations revealed that these misestimations are likely a consequence of factors not considered by the underlying microstructural models (such as intercomponent and intracompartmental kurtosis).

Conclusion

Powder‐averaged SMT and (2‐component) SM are unable to accurately report µFA and other microstructural parameters in ex vivo tissues. Improper model assumptions and constraints can significantly compromise parameter specificity. Further developments and validations are required prior to implementation of these models in clinical or preclinical research.

Measuring non-parametric distributions of intravoxel mean diffusivities using a clinical MRI scanner

We measure spectra of water mobilities (i.e., mean diffusivities) from intravoxel pools in brain tissues of healthy subjects with a non-parametric approach. Using a single-shot isotropic diffusion encoding (IDE) preparation, we eliminate signal confounds caused by anisotropic diffusion, including microscopic anisotropy, and acquire in vivo diffusion-weighted images (DWIs) over a wide range of diffusion sensitizations. We analyze the measured IDE signal decays using a regularized inverse laplace transform (ILT) to derive a probability distribution of mean diffusivities of tissue water in each voxel. Based on numerical simulations we assess the sensitivity and accuracy of our ILT analysis and optimize an experimental protocol for use with clinical MRI scanners. In vivo spectra of intravoxel mean diffusivities measured in healthy subjects generally show single-peak distributions throughout the brain parenchyma, with small differences in peak location and shape among white matter, cortical and subcortical gray matter, and cerebrospinal fluid. Mean diffusivity distributions (MDDs) with multiple peaks are observed primarily in voxels at tissue interfaces and are likely due to partial volume contributions. To quantify tissue-specific MDDs with improved statistical power, we average voxel-wise normalized MDDs in corresponding regions-of-interest (ROIs). This non-parametric, rotation-invariant assessment of isotropic diffusivities of tissue water may reflect important microstructural information, such as cell packing and cell size, and active physiological processes, such as water transport and exchange, which may enhance biological specificity in the clinical diagnosis and characterization of ischemic stroke, cancer, neuroinflammation, and neurodegenerative disorders and diseases.

Characterizing intra-axonal water diffusion with direction-averaged triple diffusion encoding MRI

For large diffusion weightings, the direction-averaged diffusion MRI (dMRI) signal from white matter is typically dominated by the contribution of water confined to axons. This fact can be exploited to characterize intra-axonal diffusion properties, which may be valuable for interpreting the biophysical meaning of diffusion changes associated with pathology. However, using just the classic Stejskal-Tanner pulse sequence, it has proven challenging to obtain reliable estimates for both the intrinsic intra-axonal diffusivity and the intra-axonal water fraction. Here we propose to apply a modification of the Stejskal-Tanner sequence designed for achieving such estimates. The key feature of the sequence is the addition of a set of extra diffusion encoding gradients that are orthogonal to the direction of the primary gradients, which corresponds to a specific type of triple diffusion encoding (TDE) MRI sequence. Given direction-averaged dMRI data for this TDE sequence, it is shown how the intra-axonal diffusivity and the intra-axonal water fraction can be determined by applying simple, analytic formulae. The method is illustrated with numerical simulations, which suggest that it should be accurate for b-values of about 4000 s/mm² or higher.

On modeling

Mapping tissue microstructure with MRI holds great promise as a noninvasive window into tissue organization at the cellular level. Having originated within the realm of diffusion NMR in the late 1970s, this field is experiencing an exponential growth in the number of publications. At the same time, model-based approaches are also increasingly incorporated into advanced MRI acquisition and reconstruction techniques. However, after about two decades of intellectual and financial investment, microstructural mapping has yet to find a single commonly accepted clinical application. Here, we suggest that slow progress in clinical translation may signify unresolved fundamental problems. We outline such problems and related practical pitfalls, as well as review strategies for developing and validating tissue microstructure models, to provoke a discussion on how to bridge the gap between our scientific aspirations and the clinical reality. We argue for recalibrating the efforts of our community toward a more systematic focus on fundamental research aimed at identifying relevant degrees of freedom affecting the measured MR signal. Such a focus is essential for realizing the truly revolutionary potential of noninvasive three-dimensional in vivo microstructural mapping.

Rotationally-invariant mapping of scalar and orientational metrics of neuronal microstructure with diffusion MRI

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.

Liquid crystal phantom for validation of microscopic diffusion anisotropy measurements on clinical MRI systems

PURPOSE

To develop a phantom for validating MRI pulse sequences and data processing methods to quantify microscopic diffusion anisotropy in the human brain.

METHODS

Using a liquid crystal consisting of water, detergent, and hydrocarbon, we designed a 0.5-L spherical phantom showing the theoretically highest possible degree of microscopic anisotropy. Data were acquired on the Connectome scanner using echo-planar imaging signal readout and diffusion encoding with axisymmetric b-tensors of varying magnitude, anisotropy, and orientation. The mean diffusivity, fractional anisotropy (FA), and microscopic FA (µFA) parameters were estimated.

RESULTS

The phantom was observed to have values of mean diffusivity similar to brain tissue, and relaxation times compatible with echo-planar imaging echo times on the order of 100 ms. The estimated values of µFA were at the theoretical maximum of 1.0, whereas the values of FA spanned the interval from 0.0 to 0.8 as a result of varying orientational order of the anisotropic domains within each voxel.

CONCLUSIONS

The proposed phantom can be manufactured by mixing three widely available chemicals in volumes comparable to a human head. The acquired data are in excellent agreement with theoretical predictions, showing that the phantom is ideal for validating methods for measuring microscopic diffusion anisotropy on clinical MRI systems. Magn Reson Med 79:1817-1828, 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-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.

Magnetic resonance temporal diffusion tensor spectroscopy of disordered anisotropic tissue

Molecular diffusion measured with diffusion weighted MRI (DWI) offers a probe for tissue microstructure. However, inferring microstructural properties from conventional DWI data is a complex inverse problem and has to account for heterogeneity in sizes, shapes and orientations of the tissue compartments contained within an imaging voxel. Alternative experimental means for disentangling the signal signatures of such features could provide a stronger link between the data and its interpretation. Double diffusion encoding (DDE) offers the possibility to factor out variation in compartment shapes from orientational dispersion of anisotropic domains by measuring the correlation between diffusivity in multiple directions. Time dependence of the diffusion is another effect reflecting the dimensions and distributions of barriers. In this paper we extend on DDE with a modified version of the oscillating gradient spin echo (OGSE) experiment, giving a basic contrast mechanism closely linked to both the temporal diffusion spectrum and the compartment anisotropy. We demonstrate our new method on post mortem brain tissue and show that we retrieve the correct temporal diffusion tensor spectrum in synthetic data from Monte Carlo simulations of random walks in a range of disordered geometries of different sizes and shapes.

Multidimensional correlation of nuclear relaxation rates and diffusion tensors for model-free investigations of heterogeneous anisotropic porous materials

Despite their widespread use in non-invasive studies of porous materials, conventional MRI methods yield ambiguous results for microscopically heterogeneous materials such as brain tissue. While the forward link between microstructure and MRI observables is well understood, the inverse problem of separating the signal contributions from different microscopic pores is notoriously difficult. Here, we introduce an experimental protocol where heterogeneity is resolved by establishing 6D correlations between the individual values of isotropic diffusivity, diffusion anisotropy, orientation of the diffusion tensor, and relaxation rates of distinct populations. Such procedure renders the acquired signal highly specific to the sample's microstructure, and allows characterization of the underlying pore space without prior assumptions on the number and nature of distinct microscopic environments. The experimental feasibility of the suggested method is demonstrated on a sample designed to mimic the properties of nerve tissue. If matched to the constraints of whole body scanners, this protocol could allow for the unconstrained determination of the different types of tissue that compose the living human brain.

Double diffusion encoding MRI for the clinic

PURPOSE

The purpose of this study is to develop double diffusion encoding (DDE) MRI methods for clinical use. Microscopic diffusion anisotropy measurements from DDE promise greater specificity to changes in tissue microstructure compared with conventional diffusion tensor imaging, but implementation of DDE sequences on whole-body MRI scanners is challenging because of the limited gradient strengths and lengthy acquisition times.

METHODS

A custom single-refocused DDE sequence was implemented on a 3T whole-body scanner. The DDE gradient orientation scheme and sequence parameters were optimized based on a Gaussian diffusion assumption. Using an optimized 5-min DDE acquisition, microscopic fractional anisotropy (μFA) maps were acquired for the first time in multiple sclerosis patients.

RESULTS

Based on simulations and in vivo human measurements, six parallel and six orthogonal diffusion gradient pairs were found to be the minimum number of diffusion gradient pairs necessary to produce a rotationally invariant measurement of μFA. Simulations showed that optimal precision and accuracy of μFA measurements were obtained using b-values between 1500 and 3000 s/mm² . The μFA maps showed improved delineation of multiple sclerosis lesions compared with conventional fractional anisotropy and distinct contrast from T₂ -weighted fluid attenuated inversion recovery and T₁ -weighted imaging.

CONCLUSION

The μFA maps can be measured using DDE in a clinical setting and may provide new opportunities for characterizing multiple sclerosis lesions and other types of tissue degeneration. Magn Reson Med 80:507-520, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Distinguishing neuronal from astrocytic subcellular microstructures using in vivo Double Diffusion Encoded 1H MRS at 21.1 T

Measuring cellular microstructures non-invasively and achieving specificity towards a cell-type population within an interrogated in vivo tissue, remains an outstanding challenge in brain research. Magnetic Resonance Spectroscopy (MRS) provides an opportunity to achieve cellular specificity via the spectral resolution of metabolites such as N-Acetylaspartate (NAA) and myo-Inositol (mI), which are considered neuronal and astrocytic markers, respectively. Yet the information typically obtained with MRS describes metabolic concentrations, diffusion coefficients or relaxation rates rather than microstructures. Understanding how these metabolites are compartmentalized is a challenging but important goal, which so far has been mainly addressed using diffusion models. Here, we present direct in vivo evidence for the confinement of NAA and mI within sub-cellular components, namely, the randomly oriented process of neurons and astrocytes, respectively. Our approach applied Relaxation Enhanced MRS at ultrahigh (21.1 T) field, and used its high ¹H sensitivity to measure restricted diffusion correlations for NAA and mI using a Double Diffusion Encoding (DDE) filter. While very low macroscopic anisotropy was revealed by spatially localized Diffusion Tensor Spectroscopy, DDE displayed characteristic amplitude modulations reporting on confinements in otherwise randomly oriented anisotropic microstructures for both metabolites. This implies that for the chosen set of parameters, the DDE measurements had a biased sensitivity towards NAA and mI sited in the more confined environments of neurites and astrocytic branches, than in the cell somata. These measurements thus provide intrinsic diffusivities and compartment diameters, and revealed subcellular neuronal and astrocytic morphologies in normal in vivo rat brains. The relevance of these measurements towards human applications—which could in turn help understand CNS plasticity as well as diagnose brain diseases—is discussed.

Magnetic resonance microdynamic imaging reveals distinct tissue microenvironments

Magnetic resonance imaging (MRI) provides a powerful set of tools with which to investigate biological tissues noninvasively and in vivo. Tissues are heterogeneous in nature; an imaging voxel contains an ensemble of different cells and extracellular matrix components. A long-standing challenge has been to infer the content of and interactions among these microscopic tissue components within a macroscopic imaging voxel. Spatially resolved multidimensional relaxation-diffusion correlation (REDCO) spectroscopy holds the potential to deliver such microdynamic information. However, to date, vast data requirements have mostly relegated these type of measurements to nuclear magnetic resonance applications and prevented them from being widely and successfully used in conjunction with imaging. By using a novel data acquisition and processing strategy in this study, spatially resolved REDCO could be performed in reasonable scanning times with excellent prospects for clinical applications. This new MR imaging framework-which we term "magnetic resonance microdynamic imaging (MRMI)"-permits the simultaneous noninvasive and model-free quantification of multiple subcellular, cellular, and interstitial tissue microenvironments within a voxel. MRMI is demonstrated with a fixed spinal cord specimen, enabling the quantification of microscopic tissue components with unprecedented specificity. Tissue components, such as axons, neuronal and glial soma, and myelin were identified on the basis of their multispectral signature within individual imaging voxels. These tissue elements could then be composed into images and be correlated with immunohistochemistry findings. MRMI provides novel image contrasts of tissue components and a new family of microdynamic biomarkers that could lead to new diagnostic imaging approaches to probe biological tissue alterations accompanied by pathological or developmental changes.

Imaging Local Diffusive Dynamics Using Diffusion Exchange Spectroscopy MRI

The movement of water between microenvironments presents a central challenge in the physics of soft matter and porous media. Diffusion exchange spectroscopy (DEXSY) is a powerful 2D nuclear magnetic resonance method for measuring such exchange, yet it is rarely used because of its long scan time requirements. Moreover, it has never been combined with magnetic resonance imaging (MRI). Using probability theory, we vastly reduce the required data, making DEXSY MRI feasible for the first time. Experiments are performed on a composite nerve tissue phantom with restricted and free water-exchanging compartments.

Ultrafast NMR diffusion measurements exploiting chirp spin echoes

Standard diffusion NMR measurements require the repetition of the experiment multiple times with varying gradient strength or diffusion delay. This makes the experiment time-consuming and restricts the use of hyperpolarized substances to boost sensitivity. We propose a novel single-scan diffusion experiment, which is based on spatial encoding of two-dimensional data, employing the spin-echoes created by two successive adiabatic frequency-swept chirp π pulses. The experiment is called ultrafast pulsed-field-gradient spin-echo (UF-PGSE). We present a rigorous derivation of the echo amplitude in the UF-PGSE experiment, justifying the theoretical basis of the method. The theory reveals also that the standard analysis of experimental data leads to a diffusion coefficient value overestimated by a few per cent. Although the overestimation is of the order of experimental error and thus insignificant in many practical applications, we propose that it can be compensated by a bipolar gradient version of the experiment, UF-BP-PGSE, or by corresponding stimulated-echo experiment, UF-BP-pulsed-field-gradient stimulated-echo. The latter also removes the effect of uniform background gradients. The experiments offer significant prospects for monitoring fast processes in real time as well as for increasing the sensitivity of experiments by several orders of magnitude by nuclear spin hyperpolarization. Furthermore, they can be applied as basic blocks in various ultrafast multidimensional Laplace NMR experiments. Copyright © 2016 John Wiley & Sons, Ltd.

Multidimensional diffusion MRI

Principles from multidimensional NMR spectroscopy, and in particular solid-state NMR, have recently been transferred to the field of diffusion MRI, offering non-invasive characterization of heterogeneous anisotropic materials, such as the human brain, at an unprecedented level of detail. Here we revisit the basic physics of solid-state NMR and diffusion MRI to pinpoint the origin of the somewhat unexpected analogy between the two fields, and provide an overview of current diffusion MRI acquisition protocols and data analysis methods to quantify the composition of heterogeneous materials in terms of diffusion tensor distributions with size, shape, and orientation dimensions. While the most advanced methods allow estimation of the complete multidimensional distributions, simpler methods focus on various projections onto lower-dimensional spaces as well as determination of means and variances rather than actual distributions. Even the less advanced methods provide simple and intuitive scalar parameters that are directly related to microstructural features that can be observed in optical microscopy images, e.g. average cell eccentricity, variance of cell density, and orientational order - properties that are inextricably entangled in conventional diffusion MRI. Key to disentangling all these microstructural features is MRI signal acquisition combining isotropic and directional dimensions, just as in the field of multidimensional solid-state NMR from which most of the ideas for the new methods are derived.

Comparative analysis of isotropic diffusion weighted imaging sequences

Visualisation of living tissue structure and function is a challenging problem of modern imaging techniques. Diffusion MRI allows one to probe in vivo structures on a micrometer scale. However, conventional diffusion measurements are time-consuming procedures, because they require several measurements with different gradient directions. Considerable time savings are therefore possible by measurement schemes that generate an isotropic diffusion weighting in a single shot. Multiple approaches for generating isotropic diffusion weighting are known and have become very popular as useful tools in clinical research. Thus, there is a strong need for a comprehensive comparison of different isotropic weighting approaches. In the present work we introduce two new sequences based on simple (co)sine modulations and compare their performance to established q-space magic-angle spinning sequences and conventional DTI, using a diffusion phantom assembled from microcapillaries and in vivo experiments at 7T. The advantages and disadvantages of all compared schemes are demonstrated and discussed.

Director orientations in lyotropic liquid crystals: diffusion MRI mapping of the Saupe order tensor

The macroscopic physical properties of a liquid crystalline material depend on both the properties of the individual crystallites and the details of their spatial arrangement. We propose a diffusion MRI method to estimate the director orientations of a lyotropic liquid crystal as a spatially resolved field of Saupe order tensors. The method relies on varying the shape of the diffusion-encoding tensor to disentangle the effects of voxel-scale director orientational order and the local diffusion anisotropy of the solvent. Proof-of-concept experiments are performed on water in lamellar and reverse hexagonal liquid crystalline systems with intricate patterns of director orientations.