Imaging diffusional variance by MRI [public] : The role of tensor-valued diffusion encoding and tissue heterogeneity

Diffusion MRI provides a non-invasive probe of tissue microstructure. We recently proposed a novel method for diffusion-weighted imaging, so-called q-space trajectory encoding, that facilitates tensor-valued diffusion encoding. This method grants access to b-tensors with multiple shapes and enables us to probe previously unexplored aspects of the tissue microstructure. Specifically, we can disentangle diffusional heterogeneity that originates from isotropic and anisotropic tissue structures; we call this diffusional variance decomposition (DIVIDE).In Paper I, we investigated the statistical uncertainty of the total diffusional variance in the healthy brain. We found that the statistical power was heterogeneous between brain regions which needs to be taken into account when interpreting results.In Paper II, we showed how spherical tensor encoding can be used to separate the total diffusional variance into its isotropic and anisotropic components. We also performed initial validation of the parameters in phantoms, and demonstrated that the imaging sequence could be implemented on a high-performance clinical MRI system. In Paper III and V, we explored DIVIDE parameters in healthy brain tissue and tumor tissue. In healthy tissue, we found that diffusion anisotropy can be probed on the microscopic scale, and that metrics of anisotropy on the voxel scale are confounded by the orientation coherence of the microscopic structures. In meningioma and glioma tumors, we found a strong association between anisotropic variance and cell eccentricity, and between isotropic variance and variable cell density. In Paper IV, we developed a method to optimize waveforms for tensor-valued diffusion encoding, and in Paper VI we demonstrated that whole-brain DIVIDE is technically feasible at most MRI systems in clinically feasible scan times

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

Purpose: Asymmetric gradient waveforms are attractive for diffusion encoding due to their superior efficiency, however, the asymmetry may cause a residual gradient moment at the end of the encoding. Depending on the experiment setup, this residual moment may cause significant signal bias and image artifacts. The purpose of this study was to develop an asymmetric gradient waveform design for tensor-valued diffusion encoding that is not affected by concomitant gradient. Methods: The Maxwell index was proposed as a scalar invariant that captures the effect of concomitant gradients and was constrained in the numerical optimization to 100 (mT/m)$^2$ms to yield Maxwell-compensated waveforms. The efficacy of this design was tested in an oil phantom, and in a healthy human brain. For reference, waveforms from literature were included in the analysis. Simulations were performed to investigate if the design was valid for a wide range of experiments and if it could predict the signal bias. Results: Maxwell-compensated waveforms showed no signal bias in oil or in the brain. By contrast, several waveforms from literature showed gross signal bias. In the brain, the bias was large enough to markedly affect both signal and parameter maps, and the bias could be accurately predicted by theory. Conclusion: Constraining the Maxwell index in the optimization of asymmetric gradient waveforms yields efficient tensor-valued encoding with concomitant gradients that have a negligible effect on the signal. This waveform design is especially relevant in combination with strong gradients, long encoding times, thick slices, simultaneous multi-slice acquisition and large/oblique FOVs

Diffusional exchange versus microscopic kurtosis from CTI: two conflicting interpretations of the same data

Correlation tensor imaging (CTI) is a new diffusion MRI framework that utilises double diffusion encoding (DDE) to resolve isotropic, anisotropic and microscopic kurtosis sources. Microscopic kurtosis in CTI is provided by the contrast between SDE and parallel DDE signals at the same b-value. Multi-Gaussian exchange (MGE) is a diffusion MRI framework that utilises DDE to measure exchange. The highest exchange sensitivity in MGE is obtained by contrasting SDE and DDE signals at the same b-value. CTI and MGE can thus be applied to analyse the same data but provide conflicting interpretations of that data. We perform Monte Carlo simulations in different geometries with varying levels of exchange to determine which approach is more compatible with the data. Simulations reveal that in all microstructures considered, CTI microscopic kurtosis drastically increases when exchange is introduced. Furthermore, in microstructures that are well-described by the multi-Gaussian assumption, CTI-estimated microscopic kurtosis increases with both the exchange rate and the mixing time, despite fulfilment of the long-mixing-time condition of CTI. Increasing the exchange rate by a factor of 2 positively biases CTI microscopic kurtosis by approximately the same factor. At a modest exchange rate of 10 /s, varying the mixing time from 12 to 100 ms increases CTI microscopic kurtosis by at least a factor of 3. To address this problem, we propose a heuristic approach to combine CTI and MGE to estimate intra-compartmental kurtosis unconfounded by exchange and demonstrate its feasibility using numerical simulations

Gradient waveform design for tensor-valued encoding in diffusion MRI

Diffusion encoding along multiple spatial directions per signal acquisition can be described in terms of a b-tensor. The benefit of tensor-valued diffusion encoding is that it unlocks the "shape of the b-tensor" as a new encoding dimension. By modulating the b-tensor shape, we can control the sensitivity to microscopic diffusion anisotropy which can be used as a contrast mechanism; a feature that is inaccessible by conventional diffusion encoding. Since imaging methods based on tensor-valued diffusion encoding are finding an increasing number of applications we are prompted to highlight the challenge of designing the optimal gradient waveforms for any given application. In this review, we first establish the basic design objectives in creating field gradient waveforms for tensor-valued diffusion MRI. We also survey additional design considerations related to limitations imposed by hardware and physiology, potential confounding effects that cannot be captured by the b-tensor, and artifacts related to the diffusion encoding waveform. Throughout, we discuss the expected compromises and tradeoffs with an aim to establish a more complete understanding of gradient waveform design and its impact on accurate measurements and interpretations of data.Comment: Invited review, submitted in May 2020 to the Journal of Neuroscience Methods. 46 pages, 9 figures, 35 equation

Spherical convolutional neural networks can improve brain microstructure estimation from diffusion MRI data

Diffusion magnetic resonance imaging is sensitive to the microstructural properties of brain tissue. However, estimating clinically and scientifically relevant microstructural properties from the measured signals remains a highly challenging inverse problem that machine learning may help solve. This study investigated if recently developed rotationally invariant spherical convolutional neural networks can improve microstructural parameter estimation. We trained a spherical convolutional neural network to predict the ground-truth parameter values from efficiently simulated noisy data and applied the trained network to imaging data acquired in a clinical setting to generate microstructural parameter maps. Our network performed better than the spherical mean technique and multi-layer perceptron, achieving higher prediction accuracy than the spherical mean technique with less rotational variance than the multi-layer perceptron. Although we focused on a constrained two-compartment model of neuronal tissue, the network and training pipeline are generalizable and can be used to estimate the parameters of any Gaussian compartment model. To highlight this, we also trained the network to predict the parameters of a three-compartment model that enables the estimation of apparent neural soma density using tensor-valued diffusion encoding

Variability in diffusion kurtosis imaging: Impact on study design, statistical power and interpretation.

Diffusion kurtosis imaging (DKI) is an emerging technique with the potential to quantify properties of tissue microstructure that may not be observable using diffusion tensor imaging (DTI). In order to help design DKI studies and improve interpretation of DKI results, we employed statistical power analysis to characterize three aspects of variability in four DKI parameters; the mean diffusivity, fractional anisotropy, mean kurtosis, and radial kurtosis. First, we quantified the variability in terms of the group size required to obtain a statistical power of 0.9. Second, we investigated the relative contribution of imaging and post-processing noise to the total variance, in order to estimate the benefits of longer scan times versus the inclusion of more subjects. Third, we evaluated the potential benefit of including additional covariates such as the size of the structure when testing for differences in group means. The analysis was performed in three major white matter structures of the brain: the superior cingulum, the corticospinal tract, and the mid-sagittal corpus callosum, extracted using diffusion tensor tractography and DKI data acquired in a healthy cohort. The results showed heterogeneous variability across and within the white matter structures. Thus, the statistical power varies depending on parameter and location, which is important to consider if a pathogenesis pattern is inferred from DKI data. In the data presented, inter-subject differences contributed more than imaging noise to the total variability, making it more efficient to include more subjects rather than extending the scan-time per subject. Finally, strong correlations between DKI parameters and the structure size were found for the cingulum and corpus callosum. Structure size should thus be considered when quantifying DKI parameters, either to control for its potentially confounding effect, or as a means of reducing unexplained variance

Assessment of Precision and Accuracy of Brain White Matter Microstructure using Combined Diffusion MRI and Relaxometry

Joint modeling of diffusion and relaxation has seen growing interest due to its potential to provide complementary information about tissue microstructure. For brain white matter, we designed an optimal diffusion-relaxometry MRI protocol that samples multiple b-values, B-tensor shapes, and echo times (TE). This variable-TE protocol (27 min) has as subsets a fixed-TE protocol (15 min) and a 2-shell dMRI protocol (7 min), both characterizing diffusion only. We assessed the sensitivity, specificity and reproducibility of these protocols with synthetic experiments and in six healthy volunteers. Compared with the fixed-TE protocol, the variable-TE protocol enables estimation of free water fractions while also capturing compartmental $T_2$ relaxation times. Jointly measuring diffusion and relaxation offers increased sensitivity and specificity to microstructure parameters in brain white matter with voxelwise coefficients of variation below 10%.Comment: 8 figure

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-minute b-tensor encoding protocol for diffusion MRI-based assessment of the microscopic anisotropy and tissue heterogeneity in a wide range of intracranial tumors. Methods: B-tensor encoding was performed in 42 patients with intracranial tumors (gliomas, meningiomas, adenomas, 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 faster 3-minute 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 \pm 0.06$ versus $0.45\pm0.08, p = 0.003)$. Metastases had higher tissue heterogeneity $(MK_I = 0.57\pm0.07)$ than both the glioblastomas $(0.44\pm0.06, p < 0.001)$ and meningiomas $(0.46\pm0.06, p = 0.03)$. Conclusion: Evaluation of the microscopic anisotropy and tissue heterogeneity in intracranial tumor patients is feasible in clinically relevant times frames.Comment: Submitted to Magnetic Resonance in Medicin

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