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

The dot-compartment revealed? Diffusion MRI with ultra-strong gradients and spherical tensor encoding in the living human brain

The so-called “dot-compartment” is conjectured in diffusion MRI to represent small spherical spaces, such as cell bodies, in which the diffusion is restricted in all directions. Previous investigations inferred its existence from data acquired with directional diffusion encoding which does not permit a straightforward separation of signals from ‘sticks’ (axons) and signals from ‘dots’. Here we combine isotropic diffusion encoding with ultra-strong diffusion gradients (240 ​mT/m) to achieve high diffusion-weightings with high signal to noise ratio, while suppressing signal arising from anisotropic water compartments with significant mobility along at least one axis (e.g., axons). A dot-compartment, defined to have apparent diffusion coefficient equal to zero and no exchange, would result in a non-decaying signal at very high b-values (b≳7000s/mm2). With this unique experimental setup, a residual yet slowly decaying signal above the noise floor for b-values as high as 15000s/mm2 was seen clearly in the cerebellar grey matter (GM), and in several white matter (WM) regions to some extent. Upper limits of the dot-signal-fraction were estimated to be 1.8% in cerebellar GM and 0.5% in WM. By relaxing the assumption of zero diffusivity, the signal at high b-values in cerebellar GM could be represented more accurately by an isotropic water pool with a low apparent diffusivity of 0.12 μm2/ms and a substantial signal fraction of 9.7%. The T2 of this component was estimated to be around 61ms. This remaining signal at high b-values has potential to serve as a novel and simple marker for isotropically-restricted water compartments in cerebellar GM.

•

The “dot-compartment” is conjectured in diffusion MRI to represent e.g. cell bodies.

•

We combine isotropic encoding with ultra-strong gradients to study the dot-compartment.

•

A slowly decaying signal for high b-values was seen in cerebellar GM.

•

An apparent diffusivity of 0.12 and signal fraction of 9.7% were estimated.

•

The signal could serve as a novel and simple marker for spherical compartments.

Diffusion Tensor MRI to Distinguish Progressive Supranuclear Palsy from α-Synucleinopathies

Background The differential diagnosis of progressive supranuclear palsy (PSP) and Lewy body disorders, which include Parkinson disease and dementia with Lewy bodies, is often challenging due to the overlapping symptoms. Purpose To develop a diagnostic tool based on diffusion tensor imaging (DTI) to distinguish between PSP and Lewy body disorders at the individual-subject level. Materials and Methods In this retrospective study, skeletonized DTI metrics were extracted from two independent data sets: the discovery cohort from the Swedish BioFINDER study and the validation cohort from the Penn Frontotemporal Degeneration Center (data collected between 2010 and 2018). Based on previous neuroimaging studies and neuropathologic evidence, a combination of regions hypothesized to be sensitive to pathologic features of PSP were identified (ie, the superior cerebellar peduncle and frontal white matter) and fractional anisotropy (FA) was used to compute an FA score for each individual. Classification performances were assessed by using logistic regression and receiver operating characteristic analysis. Results In the discovery cohort, 16 patients with PSP (mean age ± standard deviation, 73 years ± 5; eight women, eight men), 34 patients with Lewy body disorders (mean age, 71 years ± 6; 14 women, 20 men), and 44 healthy control participants (mean age, 66 years ± 8; 26 women, 18 men) were evaluated. The FA score distinguished between clinical PSP and Lewy body disorders with an area under the curve of 0.97 ± 0.04, a specificity of 91% (31 of 34), and a sensitivity of 94% (15 of 16). In the validation cohort, 34 patients with PSP (69 years ± 7; 22 women, 12 men), 25 patients with Lewy body disorders (70 years ± 7; nine women, 16 men), and 32 healthy control participants (64 years ± 7; 22 women, 10 men) were evaluated. The accuracy of the FA score was confirmed (area under the curve, 0.96 ± 0.04; specificity, 96% [24 of 25]; and sensitivity, 85% [29 of 34]). Conclusion These cross-validated findings lay the foundation for a clinical test to distinguish progressive supranuclear palsy from Lewy body disorders. © RSNA, 2019 Online supplemental material is available for this article. See also the editorial by Shah in this issue.

In vivo demonstration of microscopic anisotropy in the human kidney using multidimensional diffusion MRI

PURPOSE

To demonstrate the feasibility of multidimensional diffusion MRI to probe and quantify microscopic fractional anisotropy (µFA) in human kidneys in vivo.

METHODS

Linear tensor encoded (LTE) and spherical tensor encoded (STE) renal diffusion MRI scans were performed in 10 healthy volunteers. Respiratory triggering and image registration were used to minimize motion artefacts during the acquisition. Kidney cortex-medulla were semi-automatically segmented based on fractional anisotropy (FA) values. A model-free analysis of LTE and STE signal dependence on b-value in the renal cortex and medulla was performed. Subsequently, µFA was estimated using a single-shell approach. Finally, a comparison of conventional FA and µFA is shown.

RESULTS

The hallmark effect of µFA (divergence of LTE and STE signal with increasing b-value) was observed in all subjects. A statistically significant difference between LTE and STE signal was found in the cortex and medulla, starting from b = 750 s/mm2 and b = 500 s/mm2 , respectively. This difference was maximal at the highest b-value sampled (b = 1000 s/mm2 ) which suggests that relatively high b-values are required for µFA mapping in the kidney compared to conventional FA. Cortical and medullary µFA were, respectively, 0.53 ± 0.09 and 0.65 ± 0.05, both respectively higher than conventional FA (0.19 ± 0.02 and 0.40 ± 0.02).

CONCLUSION

The feasibility of combining LTE and STE diffusion MRI to probe and quantify µFA in human kidneys is demonstrated for the first time. By doing so, we show that novel microstructure information-not accessible by conventional diffusion encoding-can be probed by multidimensional diffusion MRI. We also identify relevant technical limitations that warrant further development of the technique for body MRI.

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.

Maxwell-compensated design of asymmetric gradient waveforms for tensor-valued diffusion encoding

PURPOSE

Diffusion encoding with asymmetric gradient waveforms is appealing because the asymmetry provides superior efficiency. However, concomitant gradients may cause a residual gradient moment at the end of the waveform, which can cause significant signal error and image artifacts. The purpose of this study was to develop an asymmetric waveform designs for tensor-valued diffusion encoding that is not sensitive to concomitant gradients.

METHODS

The "Maxwell index" was proposed as a scalar invariant to capture the effect of concomitant gradients. Optimization of "Maxwell-compensated" waveforms was performed in which this index was constrained. Resulting waveforms were compared to waveforms from literature, in terms of the measured and predicted impact of concomitant gradients, by numerical analysis as well as experiments in a phantom and in a healthy human brain.

RESULTS

Maxwell-compensated waveforms with Maxwell indices below 100 (mT/m)2 ms showed negligible signal bias in both numerical analysis and experiments. By contrast, several waveforms from literature showed gross signal bias under the same conditions, leading to a signal bias that was large enough to markedly affect parameter maps. Experimental results were accurately predicted by theory.

CONCLUSION

Constraining the Maxwell index in the optimization of asymmetric gradient waveforms yields efficient diffusion encoding that negates the effects of concomitant fields while enabling arbitrary shapes of the b-tensor. This waveform design is especially useful in combination with strong gradients, long encoding times, thick slices, simultaneous multi-slice acquisition, and large FOVs.

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.

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.

Peripheral nerve stimulation limits of a high amplitude and slew rate magnetic field gradient coil for neuroimaging

PURPOSE

To establish peripheral nerve stimulation (PNS) thresholds for an ultra-high performance magnetic field gradient subsystem (simultaneous 200-mT/m gradient amplitude and 500-T/m/s gradient slew rate; 1 MVA per axis [MAGNUS]) designed for neuroimaging with asymmetric transverse gradients and 42-cm inner diameter, and to determine PNS threshold dependencies on gender, age, patient positioning within the gradient subsystem, and anatomical landmarks.

METHODS

The MAGNUS head gradient was installed in a whole-body 3T scanner with a custom 16-rung bird-cage transmit/receive RF coil compatible with phased-array receiver brain coils. Twenty adult subjects (10 male, mean ± SD age = 40.4 ± 11.1 years) underwent the imaging and PNS study. The tests were repeated by displacing subject positions by 2-4 cm in the superior-inferior and anterior-posterior directions.

RESULTS

The x-axis (left-right) yielded mostly facial stimulation, with mean ΔG_min = 111 ± 6 mT/m, chronaxie = 766 ± 76 µsec. The z-axis (superior-inferior) yielded mostly chest/shoulder stimulation (123 ± 7 mT/m, 620 ± 62 µsec). Y-axis (anterior-posterior) stimulation was negligible. X-axis and z-axis thresholds tended to increase with age, and there was negligible dependency with gender. Translation in the inferior and posterior directions tended to increase the x-axis and z-axis thresholds, respectively. Electric field simulations showed good agreement with the PNS results. Imaging at MAGNUS gradient performance with increased PNS threshold provided a 35% reduction in noise-to-diffusion contrast as compared with whole-body performance (80 mT/m gradient amplitude, 200 T/m/sec gradient slew rate).

CONCLUSION

The PNS threshold of MAGNUS is significantly higher than that for whole-body gradients, which allows for diffusion gradients with short rise times (under 1 msec), important for interrogating brain microstructure length scales.

MK-curve - Characterizing the relation between mean kurtosis and alterations in the diffusion MRI signal

Diffusion kurtosis imaging (DKI) is a diffusion MRI (dMRI) technique to quantify brain microstructural properties. While DKI measures are sensitive to tissue alterations, they are also affected by signal alterations caused by imaging artifacts such as noise, motion and Gibbs ringing. Consequently, DKI often yields output parameter values (e.g. mean kurtosis; MK) that are implausible. These include implausible values that are outside of the range dictated by physics/biology, and visually apparent implausible values that form unexpected discontinuities, being too high or too low comparing with their neighborhood. These implausible values will introduce bias into any following data analyses (e.g. between-population statistical computation). Existing studies have attempted to correct implausible DKI parameter values in multiple ways; however, these approaches are not always effective. In this study, we propose a novel method for detecting and correcting voxels with implausible values to enable improved DKI parameter estimation. In particular, we focus on MK parameter estimation. We first characterize the relation between MK and alterations in the dMRI signal including diffusion weighted images (DWIs) and the baseline (b0) images. This is done by calculating MK for a range of synthetic DWI or b0 for each voxel, and generating curves (MK-curve) representing how alterations to the input dMRI signals affect the resulting output MK. We find that voxels with implausible MK values are more likely caused by artifacts in the b0 images than artifacts in DWIs with higher b-values. Accordingly, two characteristic b0 values, which define a range of synthetic b0 values that generate implausible MK values, are identified on the MK-curve. Based on this characterization, we propose an automatic approach for detection of voxels with implausible MK values by comparing a voxel's original b0 signal to the identified two characteristic b0 values, along with a correction strategy to replace the original b0 in each detected implausible voxel with a synthetic b0 value computed from the MK-curve. We evaluate the method on a DKI phantom dataset and dMRI datasets from the Human Connectome Project (HCP), and we compare the proposed correction method with other previously proposed correction methods. Results show that our proposed method is able to identify and correct most voxels with implausible DKI parameter values as well as voxels with implausible diffusion tensor parameter values.

Diffusion MRI microstructural models in the cervical spinal cord - Application, normative values, and correlations with histological analysis

Multi-compartment tissue modeling using diffusion magnetic resonance imaging has proven valuable in the brain, offering novel indices sensitive to the tissue microstructural environment in vivo on clinical MRI scanners. However, application, characterization, and validation of these models in the spinal cord remain relatively under-studied. In this study, we apply a diffusion "signal" model (diffusion tensor imaging, DTI) and two commonly implemented "microstructural" models (neurite orientation dispersion and density imaging, NODDI; spherical mean technique, SMT) in the human cervical spinal cord of twenty-one healthy controls. We first provide normative values of DTI, SMT, and NODDI indices in a number of white matter ascending and descending pathways, as well as various gray matter regions. We then aim to validate the sensitivity and specificity of these diffusion-derived contrasts by relating these measures to indices of the tissue microenvironment provided by a histological template. We find that DTI indices are sensitive to a number of microstructural features, but lack specificity. The microstructural models also show sensitivity to a number of microstructure features; however, they do not capture the specific microstructural features explicitly modelled. Although often regarded as a simple extension of the brain in the central nervous system, it may be necessary to re-envision, or specifically adapt, diffusion microstructural models for application to the human spinal cord with clinically feasible acquisitions - specifically, adjusting, adapting, and re-validating the modeling as it relates to both theory (i.e. relevant biology, assumptions, and signal regimes) and parameter estimation (for example challenges of acquisition, artifacts, and processing).

Linear, planar and spherical tensor-valued diffusion MRI data by free waveform encoding in healthy brain, water, oil and liquid crystals

Recently, several biophysical models and signal representations have been proposed for microstructure imaging based on tensor-valued, or multidimensional, diffusion MRI. The acquisition of the necessary data requires non-conventional pulse sequences, and data is therefore not available to the wider diffusion MRI community. To facilitate exploration and development of analysis techniques based on tensor-valued diffusion encoding, we share a comprehensive data set acquired in a healthy human brain. The data encompasses diffusion weighted images using linear, planar and spherical diffusion tensor encoding at multiple b-values and diffusion encoding directions. We also supply data acquired in several phantoms that may support validation. The data is hosted by GitHub: https://github.com/filip-szczepankiewicz/Szczepankiewicz_DIB_2019.

Resolving degeneracy in diffusion MRI biophysical model parameter estimation using double diffusion encoding

PURPOSE

Biophysical tissue models are increasingly used in the interpretation of diffusion MRI (dMRI) data, with the potential to provide specific biomarkers of brain microstructural changes. However, it has been shown recently that, in the general Standard Model, parameter estimation from dMRI data is ill-conditioned even when very high b-values are applied. We analyze this issue for the Neurite Orientation Dispersion and Density Imaging with Diffusivity Assessment (NODDIDA) model and demonstrate that its extension from single diffusion encoding (SDE) to double diffusion encoding (DDE) resolves the ill-posedness for intermediate diffusion weightings, producing an increase in accuracy and precision of the parameter estimation.

METHODS

We analyze theoretically the cumulant expansion up to fourth order in b of SDE and DDE signals. Additionally, we perform in silico experiments to compare SDE and DDE capabilities under similar noise conditions.

RESULTS

We prove analytically that DDE provides invariant information non-accessible from SDE, which makes the NODDIDA parameter estimation injective. The in silico experiments show that DDE reduces the bias and mean square error of the estimation along the whole feasible region of 5D model parameter space.

CONCLUSIONS

DDE adds additional information for estimating the model parameters, unexplored by SDE. We show, as an example, that this is sufficient to solve the previously reported degeneracies in the NODDIDA model parameter estimation.

Searching for the neurite density with diffusion MRI: Challenges for biophysical modeling

In vivo mapping of the neurite density with diffusion MRI (dMRI) is a high but challenging aim. First, it is unknown whether all neurites exhibit completely anisotropic ("stick-like") diffusion. Second, the "density" of tissue components may be confounded by non-diffusion properties such as T2 relaxation. Third, the domain of validity for the estimated parameters to serve as indices of neurite density is incompletely explored. We investigated these challenges by acquiring data with "b-tensor encoding" and multiple echo times in brain regions with low orientation coherence and in white matter lesions. Results showed that microscopic anisotropy from b-tensor data is associated with myelinated axons but not with dendrites. Furthermore, b-tensor data together with data acquired for multiple echo times showed that unbiased density estimates in white matter lesions require data-driven estimates of compartment-specific T2 values. Finally, the "stick" fractions of different biophysical models could generally not serve as neurite density indices across the healthy brain and white matter lesions, where outcomes of comparisons depended on the choice of constraints. In particular, constraining compartment-specific T2 values was ambiguous in the healthy brain and had a large impact on estimated values. In summary, estimating neurite density generally requires accounting for different diffusion and/or T2 properties between axons and dendrites. Constrained "index" parameters could be valid within limited domains that should be delineated by future studies.

Diffusion MRI in the brain - Theory and concepts

Over the past two decades, diffusion MRI has become an essential tool in neuroimaging investigations. This is due to its sensitivity to the motion of water molecules as they diffuse through the microstructural environment, allowing diffusion MRI to be used as a 'probe' of tissue microstructure. Furthermore, this sensitivity is strongly direction-dependent, notably in brain white matter, due to the alignment of structures that restrict or hinder the motion of water molecules, notably axonal membranes. This provides a means of inferring the orientation of fibres in vivo, and by use of appropriate fibre-tracking algorithms, of delineating the path of white matter tracts in the brain. The ability to perform so-called tractography in humans in vivo non-invasively is unique to diffusion MRI, and is now used in applications such as neurosurgery planning and more broadly within investigations of brain connectomics. This review describes the theory and concepts of diffusion MRI and describes its most important areas of application in the brain, with a strong focus on tractography.

Monte Carlo Simulations of Water Exchange Through Myelin Wraps: Implications for Diffusion MRI

Diffusion magnetic resonance imaging (dMRI) yields parameters sensitive to brain tissue microstructure. A structurally important aspect of this microstructure is the myelin wrapping around the axons. This paper investigated the forward problem concerning whether water exchange via the spiraling structure of the myelin can meaningfully contribute to the signal in dMRI. Monte Carlo simulations were performed in a system with intra-axonal, myelin, and extra-axonal compartments. Diffusion in the myelin was simulated as a spiral wrapping the axon, with a custom number of wraps. Exchange (or intra-axonal residence) times were analyzed for various number of wraps and axon diameters. Pulsed gradient sequences were employed to simulate the dMRI signal, which was analyzed using different methods. Diffusional kurtosis imaging analysis yielded the radial diffusivity (RD) and radial kurtosis (RK), while the two-compartment Kärger model yielded estimates the intra-axonal volume fraction ( ν _ic ) and exchange time ( τ ). Results showed that τ was on the sub-second level for geometries with axon diameters below 1.0 μ m and less than eight wraps. Otherwise, exchange was negligible compared to typical experimental durations, with τ of seconds or longer. In situations where exchange influenced the signal, estimates of RK and ν _ic increased with the number of wraps, while RD decreased. τ estimates from simulated signals were in agreement with predicted ones. In conclusion, exchange through spiraling myelin permits sub-second τ for small diameters and low number of wraps. Such conditions may arise in the developing brain or in neurodegenerative disease, and thus the results could aid the interpretation of dMRI studies.

Estimation of diffusion, perfusion and fractional volumes using a multi-compartment relaxation-compensated intravoxel incoherent motion (IVIM) signal model

Compartmental diffusion MRI models that account for intravoxel incoherent motion (IVIM) of blood perfusion allow for estimation of the fractional volume of the microvascular compartment. Conventional IVIM models are known to be biased by not accounting for partial volume effects caused by free water and cerebrospinal fluid (CSF), or for tissue-dependent relaxation effects. In this work, a three-compartment model (tissue, free water and blood) that includes relaxation terms is introduced. To estimate the model parameters, in vivo human data were collected with multiple echo times (TE), inversion times (TI) and b-values, which allowed a direct relaxation estimate alongside estimation of perfusion, diffusion and fractional volume parameters. Compared to conventional two-compartment models (with and without relaxation compensation), the three-compartment model showed less effects of CSF contamination. The proposed model yielded significantly different volume fractions of blood and tissue compared to the non-relaxation-compensated model, as well as to the conventional two-compartment model, suggesting that previously reported parameter ranges, using models that do not account for relaxation, should be reconsidered.

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.

Advances in computational and statistical diffusion MRI

Computational methods are crucial for the analysis of diffusion magnetic resonance imaging (MRI) of the brain. Computational diffusion MRI can provide rich information at many size scales, including local microstructure measures such as diffusion anisotropies or apparent axon diameters, whole-brain connectivity information that describes the brain's wiring diagram and population-based studies in health and disease. Many of the diffusion MRI analyses performed today were not possible five, ten or twenty years ago, due to the requirements for large amounts of computer memory or processor time. In addition, mathematical frameworks had to be developed or adapted from other fields to create new ways to analyze diffusion MRI data. The purpose of this review is to highlight recent computational and statistical advances in diffusion MRI and to put these advances into context by comparison with the more traditional computational methods that are in popular clinical and scientific use. We aim to provide a high-level overview of interest to diffusion MRI researchers, with a more in-depth treatment to illustrate selected computational advances.

Challenges in diffusion MRI tractography - Lessons learned from international benchmark competitions

Diffusion MRI (dMRI) fiber tractography has become a pillar of the neuroimaging community due to its ability to noninvasively map the structural connectivity of the brain. Despite widespread use in clinical and research domains, these methods suffer from several potential drawbacks or limitations. Thus, validating the accuracy and reproducibility of techniques is critical for sound scientific conclusions and effective clinical outcomes. Towards this end, a number of international benchmark competitions, or "challenges", has been organized by the diffusion MRI community in order to investigate the reliability of the tractography process by providing a platform to compare algorithms and results in a fair manner, and evaluate common and emerging algorithms in an effort to advance the state of the field. In this paper, we summarize the lessons from a decade of challenges in tractography, and give perspective on the past, present, and future "challenges" that the field of diffusion tractography faces.

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.

Building connectomes using diffusion MRI: why, how and but

Why has diffusion MRI become a principal modality for mapping connectomes in vivo? How do different image acquisition parameters, fiber tracking algorithms and other methodological choices affect connectome estimation? What are the main factors that dictate the success and failure of connectome reconstruction? These are some of the key questions that we aim to address in this review. We provide an overview of the key methods that can be used to estimate the nodes and edges of macroscale connectomes, and we discuss open problems and inherent limitations. We argue that diffusion MRI-based connectome mapping methods are still in their infancy and caution against blind application of deep white matter tractography due to the challenges inherent to connectome reconstruction. We review a number of studies that provide evidence of useful microstructural and network properties that can be extracted in various independent and biologically relevant contexts. Finally, we highlight some of the key deficiencies of current macroscale connectome mapping methodologies and motivate future developments.

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.

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.

Tensor-valued diffusion encoding for diffusional variance decomposition (DIVIDE): Technical feasibility in clinical MRI systems

Microstructure imaging techniques based on tensor-valued diffusion encoding have gained popularity within the MRI research community. Unlike conventional diffusion encoding-applied along a single direction in each shot-tensor-valued encoding employs diffusion encoding along multiple directions within a single preparation of the signal. The benefit is that such encoding may probe tissue features that are not accessible by conventional encoding. For example, diffusional variance decomposition (DIVIDE) takes advantage of tensor-valued encoding to probe microscopic diffusion anisotropy independent of orientation coherence. The drawback is that tensor-valued encoding generally requires gradient waveforms that are more demanding on hardware; it has therefore been used primarily in MRI systems with relatively high performance. The purpose of this work was to explore tensor-valued diffusion encoding on clinical MRI systems with varying performance to test its technical feasibility within the context of DIVIDE. We performed whole-brain imaging with linear and spherical b-tensor encoding at field strengths between 1.5 and 7 T, and at maximal gradient amplitudes between 45 and 80 mT/m. Asymmetric gradient waveforms were optimized numerically to yield b-values up to 2 ms/μm2. Technical feasibility was assessed in terms of the repeatability, SNR, and quality of DIVIDE parameter maps. Variable system performance resulted in echo times between 83 to 115 ms and total acquisition times of 6 to 9 minutes when using 80 signal samples and resolution 2×2×4 mm3. As expected, the repeatability, signal-to-noise ratio and parameter map quality depended on hardware performance. We conclude that tensor-valued encoding is feasible for a wide range of MRI systems-even at 1.5 T with maximal gradient waveform amplitudes of 33 mT/m-and baseline experimental design and quality parameters for all included configurations. This demonstrates that tissue features, beyond those accessible by conventional diffusion encoding, can be explored on a wide range of MRI systems.

Orientationally-averaged diffusion-attenuated magnetic resonance signal for locally-anisotropic diffusion

Diffusion-attenuated MR signal for heterogeneous media has been represented as a sum of signals from anisotropic Gaussian sub-domains to the extent that this approximation is permissible. Any effect of macroscopic (global or ensemble) anisotropy in the signal can be removed by averaging the signal values obtained by differently oriented experimental schemes. The resulting average signal is identical to what one would get if the micro-domains are isotropically (e.g., randomly) distributed with respect to orientation, which is the case for “powdered” specimens. We provide exact expressions for the orientationally-averaged signal obtained via general gradient waveforms when the microdomains are characterized by a general diffusion tensor possibly featuring three distinct eigenvalues. This extends earlier results which covered only axisymmetric diffusion as well as measurement tensors. Our results are expected to be useful in not only multidimensional diffusion MR but also solid-state NMR spectroscopy due to the mathematical similarities in the two fields.

Eddy current nulled constrained optimization of isotropic diffusion encoding gradient waveforms

PURPOSE

Isotropic diffusion encoding efficiently encodes additional microstructural information in combination with conventional linear diffusion encoding. However, the gradient-intensive isotropic diffusion waveforms generate significant eddy currents, which cause image distortions. The purpose of this study is to present a method for designing isotropic diffusion encoding waveforms with intrinsic eddy current nulling.

METHODS

Eddy current nulled gradient waveforms were designed using a constrained optimization waveform for a 3T GE Premier MRI system. Encoding waveforms were designed for a variety of eddy current null times and sequence timings to evaluate the achievable b-value. Waveforms were also tested with both eddy current nulling and concomitant field compensation. Distortion reduction was tested in both phantoms and the in vivo human brain.

RESULTS

The feasibility of isotropic diffusion encoding with intrinsic correction of eddy current distortion and signal bias from concomitant fields was demonstrated. The constrained optimization algorithm produced gradient waveforms with the specified eddy current null times. The reduction in the achievable diffusion weighting was dependent on the number of eddy current null times. A reduction in the eddy current-induced image distortions was observed in both phantoms and in vivo human subjects.

CONCLUSION

The proposed framework allows the design of isotropic diffusion-encoding sequences with reduced image distortion.

Effects of nongaussian diffusion on "isotropic diffusion" measurements: An ex-vivo microimaging and simulation study

Designing novel diffusion-weighted pulse sequences to probe tissue microstructure beyond the conventional Stejskal-Tanner family is currently of broad interest. One such technique, multidimensional diffusion MRI, has been recently proposed to afford model-free decomposition of diffusion signal kurtosis into terms originating from either ensemble variance of isotropic diffusivity or microscopic diffusion anisotropy. This ability rests on the assumption that diffusion can be described as a sum of multiple Gaussian compartments, but this is often not strictly fulfilled. The effects of nongaussian diffusion on single shot isotropic diffusion sequences were first considered in detail by de Swiet and Mitra in 1996. They showed theoretically that anisotropic compartments lead to anisotropic time dependence of the diffusion tensors, which causes the measured isotropic diffusivity to depend on gradient frame orientation. Here we show how such deviations from the multiple Gaussian compartments assumption conflates orientation dispersion with ensemble variance in isotropic diffusivity. Second, we consider additional contributions to the apparent variance in isotropic diffusivity arising due to intracompartmental kurtosis. These will likewise depend on gradient frame orientation. We illustrate the potential importance of these confounds with analytical expressions, numerical simulations in simple model geometries, and microimaging experiments in fixed spinal cord using isotropic diffusion encoding waveforms with 7.5 ms duration and 3000 mT/m maximum amplitude.

Review of diffusion MRI studies in chronic white matter diseases

Diffusion MRI studies characterizing the changes in white matter (WM) due to vascular cognitive impairment, which includes all forms of small vessel disease are reviewed. We reviewed the usefulness of diffusion methods in discriminating the affected WM regions and its relation to cognitive impairment. These studies were categorized based on the diffusion MRI techniques used. The most common method was the diffusion tensor imaging, whereas other methods included diffusion weighted imaging, diffusion kurtosis imaging, intravoxel incoherent motion, and studies based on diffusion tractography. The diffusion measures showed correlation with cognitive scores and disease progression, with mean diffusivity being the most robust parameter. Future studies should focus on incorporating multi-compartment and higher order diffusion models, which can handle the presence of multiple and crossing fibers inside a voxel.

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.

Isotropically weighted intravoxel incoherent motion brain imaging at 7T

Perfusion magnetic resonance imaging (MRI) is a promising non-invasive technique providing insights regarding the brain's microvascular architecture in vivo. The scalar perfusion metrics can be used for quantitative diagnostics of various brain abnormalities, in particular, in the stroke cases and tumours. However, conventional MRI-based perfusion approaches such as dynamic contrast-enhanced perfusion imaging or arterial spin labelling have a few weaknesses, for instance, contrast agent deposition, low signal-to-noise ratio, limited temporal and spatial resolution, and specific absorption rate constraints. As an alternative, the intravoxel incoherent motion (IVIM) approach exploits an extension of diffusion MRI in order to estimate perfusion parameters in the human brain. Application of IVIM imaging at ultra-high field MRI might employ the advantage of a higher signal-to-noise ratio, and thereby the use of higher spatial and temporal resolutions. In the present work, we demonstrate an application of recently developed isotropic diffusion weighted sequences to the evaluation of IVIM parameters at an ultra-high 7T field. The used sequence exhibits high immunity to image degrading factors and allows one to acquire the data in a fast and efficient way. Utilising the bi-exponential fitting model of the signal attenuation, we performed an extensive analysis of the IVIM scalar metrics obtained by a isotropic diffusion weighted sequence in vivo and compared results with a conventional pulsed gradient sequence at 7T. In order to evaluate a possible metric bias originating from blood flows, we additionally used a truncated b-value protocol (b-values from 100 to 200 s/mm² with the step 20 s/mm²) accompanied to the full range (b-values from 0 to 200 s/mm²). The IVIM scalar metrics have been assessed and analysed together with a large and middle vessel density atlas of the human brain. We found that the diffusion coefficients and perfusion fractions of the voxels consisting of large and middle vessels have higher values in contrast to other tissues. Additionally, we did not find a strong dependence of the IVIM metrics on the density values of the vessel atlas. Perspectives and limitations of the developed isotropic diffusion weighted perfusion are presented and discussed.

Advances in microstructural diffusion neuroimaging for psychiatric disorders

Understanding the neuropathological underpinnings of mental disorders such as schizophrenia, major depression, and bipolar disorder is an essential step towards the development of targeted treatments. Diffusion MRI studies utilizing the diffusion tensor imaging (DTI) model have been extremely successful to date in identifying microstructural brain abnormalities in individuals suffering from mental illness, especially in regions of white matter, although identified abnormalities have been biologically non-specific. Building on DTI's success, in recent years more advanced diffusion MRI methods have been developed and applied to the study of psychiatric populations, with the aim of offering increased sensitivity to subtle neurological abnormalities, as well as improved specificity to candidate pathologies such as demyelination and neuroinflammation. These advanced methods, however, usually come at the cost of prolonged imaging sequences or reduced signal to noise, and they are more difficult to evaluate compared with the more simplified approach taken by the now common DTI model. To date, a limited number of advanced diffusion MRI methods have been employed to study schizophrenia, major depression and bipolar disorder populations. In this review we survey these studies, compare findings across diverse methods, discuss the main benefits and limitations of the different methods, and assess the extent to which the application of more advanced diffusion imaging approaches has led to novel and transformative information with regards to our ability to better understand the etiology and pathology of mental disorders.

Microstructural imaging of the human brain with a 'super-scanner': 10 key advantages of ultra-strong gradients for diffusion MRI

The key component of a microstructural diffusion MRI 'super-scanner' is a dedicated high-strength gradient system that enables stronger diffusion weightings per unit time compared to conventional gradient designs. This can, in turn, drastically shorten the time needed for diffusion encoding, increase the signal-to-noise ratio, and facilitate measurements at shorter diffusion times. This review, written from the perspective of the UK National Facility for In Vivo MR Imaging of Human Tissue Microstructure, an initiative to establish a shared 300 mT/m-gradient facility amongst the microstructural imaging community, describes ten advantages of ultra-strong gradients for microstructural imaging. Specifically, we will discuss how the increase of the accessible measurement space compared to a lower-gradient systems (in terms of Δ, b-value, and TE) can accelerate developments in the areas of 1) axon diameter distribution mapping; 2) microstructural parameter estimation; 3) mapping micro-vs macroscopic anisotropy features with gradient waveforms beyond a single pair of pulsed-gradients; 4) multi-contrast experiments, e.g. diffusion-relaxometry; 5) tractography and high-resolution imaging in vivo and 6) post mortem; 7) diffusion-weighted spectroscopy of metabolites other than water; 8) tumour characterisation; 9) functional diffusion MRI; and 10) quality enhancement of images acquired on lower-gradient systems. We finally discuss practical barriers in the use of ultra-strong gradients, and provide an outlook on the next generation of 'super-scanners'.

A fractional motion diffusion model for a twice-refocused spin-echo pulse sequence

The purpose of this study was to develop an analytical expression for a fractional motion (FM) diffusion model to characterize diffusion-induced signal attenuation in a twice-refocused spin-echo (TRSE) sequence that is resilient to eddy currents, and to demonstrate its applicability to human brain imaging in vivo. Based on the FM theory, which provides a unified statistical description for Langevin motions, the diffusion-weighted (DW) MR signal was measured with a TRSE sequence that balances the concomitant gradients. The analytical expression was fitted to a set of DW images acquired with 14 b-values (0-4000 s/mm² ) from a total of 10 healthy human subjects at 3 T, yielding three FM parameter maps based on anomalous diffusion coefficient D_{φ, ψ} , diffusion increment variance φ, and diffusion correlation ψ, respectively. These parameters were used to characterize different brain regions in gray matter (GM), white matter (WM), and cerebrospinal fluid. The analytical expression for the TRSE-based FM model accurately described diffusion signal attenuation in healthy brain tissues at high b-values. TRSE's robustness against eddy currents was illustrated by comparing results from an expression for a conventional Stejskal-Tanner sequence. The TRSE-based FM model also produced consistent GM-WM contrast (p < 0.01) across all brain regions studied, whereas the consistency was not observed with the Stejskal-Tanner-based FM model. This new analytical expression is expected to enable further investigations to probe tissue structures by exploiting anomalous diffusion properties without being hindered by eddy-current perturbations at high b-values.

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.

Relevance of time-dependence for clinically viable diffusion imaging of the spinal cord

Purpose

Time‐dependence is a key feature of the diffusion‐weighted (DW) signal, knowledge of which informs biophysical modelling. Here, we study time‐dependence in the human spinal cord, as its axonal structure is specific and different from the brain.

Methods

We run Monte Carlo simulations using a synthetic model of spinal cord white matter (WM) (large axons), and of brain WM (smaller axons). Furthermore, we study clinically feasible multi‐shell DW scans of the cervical spinal cord (b = 0; b = 711 s mm⁻²; b = 2855 s mm⁻²), obtained using three diffusion times (Δ of 29, 52 and 76 ms) from three volunteers.

Results

Both intra‐/extra‐axonal perpendicular diffusivities and kurtosis excess show time‐dependence in our synthetic spinal cord model. This time‐dependence is reflected mostly in the intra‐axonal perpendicular DW signal, which also exhibits strong decay, unlike our brain model. Time‐dependence of the total DW signal appears detectable in the presence of noise in our synthetic spinal cord model, but not in the brain. In WM in vivo, we observe time‐dependent macroscopic and microscopic diffusivities and diffusion kurtosis, NODDI and two‐compartment SMT metrics. Accounting for large axon calibers improves fitting of multi‐compartment models to a minor extent.

Conclusions

Time‐dependence of clinically viable DW MRI metrics can be detected in vivo in spinal cord WM, thus providing new opportunities for the non‐invasive estimation of microstructural properties. The time‐dependence of the perpendicular DW signal may feature strong intra‐axonal contributions due to large spinal axon caliber. Hence, a popular model known as “stick” (zero‐radius cylinder) may be sub‐optimal to describe signals from the largest spinal axons.

Development of a robust diffusion-MR elastography (dMRE) technique to mitigate intravoxel phase dispersion

Diffusion-magnetic resonance elastography (dMRE) is an emerging practical technique that can acquire diffusion magnetic resonance imaging and MRE simultaneously. However, a signal loss attributable to intravoxel phase dispersion (IVPD) interferes with the calculation of the apparent diffusion coefficient (ADC). This study presents an approach to dMRE that reduces the influence of IVPD by introducing a new pulse sequence. The existing and proposed techniques were performed using a phantom comprising five rods with different elasticities at 60 Hz vibration to investigate the accuracy of previous and proposed dMRE techniques. The measures of ADC and stiffness, obtained by using both dMRE techniques, were compared with conventional spin-echo (SE) diffusion and SE-MRE. Then, we evaluated those differences by using the mean of absolute differences (MAD) in each rod within the phantom. The results of the MAD of the stiffness from both dMRE techniques showed almost no difference. In contrast, the value of the ADC MAD (MAD ≒ 0.16 × 10^-3 mm²/s), obtained in the soft region within the phantom with the previous dMRE technique, was large. This value was about 2.7 times that of the value produced by the proposed dMRE technique. This difference must reflect the degree of influence of IVPD in both techniques. These results demonstrate that our dMRE technique is a robust method for addressing the signal loss attributable to IVPD.

Modeling white matter microstructure with fiber ball imaging

Fiber ball imaging (FBI) provides a means of calculating the fiber orientation density function (fODF) in white matter from diffusion MRI (dMRI) data obtained over a spherical shell with a b-value of about 4000 s/mm2 or higher. By supplementing this FBI-derived fODF with dMRI data acquired for two lower b-value shells, it is shown that several microstructural parameters may be estimated, including the axonal water fraction (AWF) and the intrinsic intra-axonal diffusivity. This fiber ball white matter (FBWM) modeling method is demonstrated for dMRI data acquired from healthy volunteers, and the results are compared with those of the white matter tract integrity (WMTI) method. Both the AWF and the intra-axonal diffusivity obtained with FBWM are found to be significantly larger than for WMTI, with the FBWM values for the intra-axonal diffusivity being more consistent with recent results obtained using isotropic diffusion weighting. An important practical advantage of FBWM is that the only nonlinear fitting required is the minimization of a cost function with just a single free parameter, which facilitates the implementation of efficient and robust numerical routines.

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.

Cumulant expansions for measuring water exchange using diffusion MRI

The rate of water exchange across cell membranes is a parameter of biological interest and can be measured by diffusion magnetic resonance imaging (dMRI). In this work, we investigate a stochastic model for the diffusion-and-exchange of water molecules. This model provides a general solution for the temporal evolution of dMRI signal using any type of gradient waveform, thereby generalizing the signal expressions for the Kärger model. Moreover, we also derive a general nth order cumulant expansion of the dMRI signal accounting for water exchange, which has not been explored in earlier studies. Based on this analytical expression, we compute the cumulant expansion for dMRI signals for the special case of single diffusion encoding (SDE) and double diffusion encoding (DDE) sequences. Our results provide a theoretical guideline on optimizing experimental parameters for SDE and DDE sequences, respectively. Moreover, we show that DDE signals are more sensitive to water exchange at short-time scale but provide less attenuation at long-time scale than SDE signals. Our theoretical analysis is also validated using Monte Carlo simulations on synthetic structures.

White matter biomarkers from diffusion MRI

As part of an issue celebrating 2 decades of Joseph Ackerman editing the Journal of Magnetic Resonance, this paper reviews recent progress in one of the many areas in which Ackerman and his lab has made significant contributions: NMR measurement of diffusion in biological media, specifically in brain tissue. NMR diffusion signals display exquisite sensitivity to tissue microstructure, and have the potential to offer quantitative and specific information on the cellular scale orders of magnitude below nominal image resolution when combined with biophysical modeling. Here, I offer a personal perspective on some recent advances in diffusion imaging, from diffusion kurtosis imaging to microstructural modeling, and the connection between the two. A new result on the estimation accuracy of axial and radial kurtosis with axially symmetric DKI is presented. I moreover touch upon recently suggested generalized diffusion sequences, promising to offer independent microstructural information. We discuss the need and some methods for validation, and end with an outlook on some promising future directions.

Advanced diffusion imaging for assessing normal white matter development in neonates and characterizing aberrant development in congenital heart disease

Background

Elucidating developmental trajectories of white matter (WM) microstructure is critically important for understanding normal development and regional vulnerabilities in several brain disorders. Diffusion Weighted Imaging (DWI) is currently the method of choice for in-vivo white matter assessment. A majority of neonatal studies use the standard Diffusion Tensor Imaging (DTI) model although more advanced models such as the Neurite Orientation Dispersion and Density Imaging (NODDI) model and the Gaussian Mixture Model (GMM) have been used in adult population. In this study, we compare the ability of these three diffusion models to detect regional white matter maturation in typically developing control (TDC) neonates and regional abnormalities in neonates with congenital heart disease (CHD).

Methods

Multiple b-value diffusion Magnetic Resonance Imaging (dMRI) data were acquired from TDC neonates (N = 16) at 38 to 47 gestational weeks (GW) and CHD neonates (N = 19) aged 37 weeks to 41 weeks. Measures calculated from the diffusion signal included not only Mean Diffusivity (MD) and Fractional Anisotropy (FA) derived from the standard DTI model, but also three advanced diffusion measures, namely, the fiber Orientation Dispersion Index (ODI), the isotropic volume fraction (Viso), and the intracellular volume fraction (Vic) derived from the NODDI model. Further, we used two novel measures from a non-parametric GMM, namely the Return-to-Origin Probability (RTOP) and Return-to-Axis Probability (RTAP), which are sensitive to axonal/cellular volume and density respectively. Using atlas-based registration, 22 white matter regions (6 projection, 4 association, and 1 callosal pathways bilaterally in each hemisphere) were selected and the mean value of all 7 measures were calculated in each region. These values were used as dependent variables, with GW as the independent variable in a linear regression model. Finally, we compared CHD and TDC groups on these measures in each ROI after removing age-related trends from both the groups.

Results

Linear analysis in the TDC population revealed significant correlations with GW (age) in 12 projection pathways for MD, Vic, RTAP, and 11 pathways for RTOP. Several association pathways were also significantly correlated with GW for MD, Vic, RTAP, and RTOP. The right callosal pathway was significantly correlated with GW for Vic. Consistent with the pathophysiology of altered development in CHD, diffusion measures demonstrated differences in the association pathways involved in language systems, namely the Uncinate Fasciculus (UF), the Inferior Fronto-occipital Fasciculus (IFOF), and the Superior Longitudinal Fasciculus (SLF). Overall, the group comparison between CHD and TDC revealed lower FA, Vic, RTAP, and RTOP for CHD bilaterally in the a) UF, b) Corpus Callosum (CC), and c) Superior Fronto-Occipital Fasciculus (SFOF). Moreover, FA was lower for CHD in the a) left SLF, b) bilateral Anterior Corona Radiata (ACR) and left Retrolenticular part of the Internal Capsule (RIC). Vic was also lower for CHD in the left Posterior Limb of the Internal Capsule (PLIC). ODI was higher for CHD in the left CC. RTAP was lower for CHD in the left IFOF, while RTOP was lower in CHD in the: a) left ACR, b) left IFOF and c) right Anterior Limb of the Internal Capsule (ALIC).

Conclusion

In this study, all three methods revealed the expected changes in the WM regions during the early postnatal weeks; however, GMM outperformed DTI and NODDI as it showed significantly larger effect sizes while detecting differences between the TDC and CHD neonates. Future studies based on a larger sample are needed to confirm these results and to explore clinical correlates.

Bayesian uncertainty quantification in linear models for diffusion MRI

Diffusion MRI (dMRI) is a valuable tool in the assessment of tissue microstructure. By fitting a model to the dMRI signal it is possible to derive various quantitative features. Several of the most popular dMRI signal models are expansions in an appropriately chosen basis, where the coefficients are determined using some variation of least-squares. However, such approaches lack any notion of uncertainty, which could be valuable in e.g. group analyses. In this work, we use a probabilistic interpretation of linear least-squares methods to recast popular dMRI models as Bayesian ones. This makes it possible to quantify the uncertainty of any derived quantity. In particular, for quantities that are affine functions of the coefficients, the posterior distribution can be expressed in closed-form. We simulated measurements from single- and double-tensor models where the correct values of several quantities are known, to validate that the theoretically derived quantiles agree with those observed empirically. We included results from residual bootstrap for comparison and found good agreement. The validation employed several different models: Diffusion Tensor Imaging (DTI), Mean Apparent Propagator MRI (MAP-MRI) and Constrained Spherical Deconvolution (CSD). We also used in vivo data to visualize maps of quantitative features and corresponding uncertainties, and to show how our approach can be used in a group analysis to downweight subjects with high uncertainty. In summary, we convert successful linear models for dMRI signal estimation to probabilistic models, capable of accurate uncertainty quantification.

Influence of the size and curvedness of neural projections on the orientationally averaged diffusion MR signal

Neuronal and glial projections can be envisioned to be tubes of infinitesimal diameter as far as diffusion magnetic resonance (MR) measurements via clinical scanners are concerned. Recent experimental studies indicate that the decay of the orientationally-averaged signal in white-matter may be characterized by the power-law, Ē(q) ∝ q^-1, where q is the wavenumber determined by the parameters of the pulsed field gradient measurements. One particular study by McKinnon et al. [1] reports a distinctively faster decay in gray-matter. Here, we assess the role of the size and curvature of the neurites and glial arborizations in these experimental findings. To this end, we studied the signal decay for diffusion along general curves at all three temporal regimes of the traditional pulsed field gradient measurements. We show that for curvy projections, employment of longer pulse durations leads to a disappearance of the q^-1 decay, while such decay is robust when narrow gradient pulses are used. Thus, in clinical acquisitions, the lack of such a decay for a fibrous specimen can be seen as indicative of fibers that are curved. We note that the above discussion is valid for an intermediate range of q-values as the true asymptotic behavior of the signal decay is Ē(q) ∝ q^-4 for narrow pulses (through Debye-Porod law) or steeper for longer pulses. This study is expected to provide insights for interpreting the diffusion-weighted images of the central nervous system and aid in the design of acquisition strategies.