Simultaneous magnetic resonance diffusion and pseudo-diffusion tensor imaging

Diffusion-derived intravoxel-incoherent motion anisotropy relates to CSF and blood flow

This study investigates the feasibility of multi-b-value, multi-directional diffusion MRI for assessing the anisotropy of the cerebral pseudo-diffusion (D*)-tensor. We examine D*-tensor's potential to (1) reflect CSF and blood flow, and (2) detect microvascular architectural alterations in cerebral small vessel disease (cSVD) and aging.

METHODS

Multi-b-value diffusion MRI was acquired in 32 gradient directions for 11 healthy volunteers, and in six directions for 29 patients with cSVD and 14 controls at 3 T. A physics-informed neural network was used to estimate intravoxel incoherent motion (IVIM)-DTI model parameters, including the parenchymal slow diffusion (D-)tensor and the pseudo-diffusion (D*)-tensor, from which the fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD), and radial diffusivity (RD) were derived. Comparisons of D*-tensor metrics were made between lateral, third, and fourth ventricles and between the middle cerebral arteries and superior sagittal sinus. Group differences in D*-tensor metrics in normal-appearing white matter were analyzed using multivariable linear regression, correcting for age and sex.

RESULTS

D*-anisotropy aligned well with CSF flow and arterial blood flow. FA(D*), MD(D*), AD(D*), and RD(D*) were highest in the third, moderate in the fourth, and lowest in the lateral ventricles. The arteries showed higher MD(D*), AD(D*), and RD(D*) than the sagittal sinus. Higher FA(D*) in the normal-appearing white matter was related to cSVD diagnosis and older age, suggesting microvascular architecture alterations.

CONCLUSION

Multi-b-value, multi-directional diffusion analysis using the IVIM-DTI model enables assessment of the cerebral microstructure, fluid flow, and microvascular architecture, providing information on neurodegeneration, glymphatic waste clearance, and the vasculature in one measurement.

SpinDoctor-IVIM: A virtual imaging framework for intravoxel incoherent motion MRI

Intravoxel incoherent motion (IVIM) imaging is increasingly recognised as an important tool in clinical MRI, where tissue perfusion and diffusion information can aid disease diagnosis, monitoring of patient recovery, and treatment outcome assessment. Currently, the discovery of biomarkers based on IVIM imaging, similar to other medical imaging modalities, is dependent on long preclinical and clinical validation pathways to link observable markers derived from images with the underlying pathophysiological mechanisms. To speed up this process, virtual IVIM imaging is proposed. This approach provides an efficient virtual imaging tool to design, evaluate, and optimise novel approaches for IVIM imaging. In this work, virtual IVIM imaging is developed through a new finite element solver, SpinDoctor-IVIM, which extends SpinDoctor, a diffusion MRI simulation toolbox. SpinDoctor-IVIM simulates IVIM imaging signals by solving the generalised Bloch-Torrey partial differential equation. The input velocity to SpinDoctor-IVIM is computed using HemeLB, an established Lattice Boltzmann blood flow simulator. Contrary to previous approaches, SpinDoctor-IVIM accounts for volumetric microvasculature during blood flow simulations, incorporates diffusion phenomena in the intravascular space, and accounts for the permeability between the intravascular and extravascular spaces. The above-mentioned features of the proposed framework are illustrated with simulations on a realistic microvasculature model.

Integrated intravoxel incoherent motion tensor and diffusion tensor brain MRI in a single fast acquisition

The acquisition of intravoxel incoherent motion (IVIM) data and diffusion tensor imaging (DTI) data from the brain can be integrated into a single measurement, which offers the possibility to determine orientation-dependent (tensorial) perfusion parameters in addition to established IVIM and DTI parameters. The purpose of this study was to evaluate the feasibility of such a protocol with a clinically feasible scan time below 6 min and to use a model-selection approach to find a set of DTI and IVIM tensor parameters that most adequately describes the acquired data. Diffusion-weighted images of the brain were acquired at 3 T in 20 elderly participants with cerebral small vessel disease using a multiband echoplanar imaging sequence with 15 b-values between 0 and 1000 s/mm² and six non-collinear diffusion gradient directions for each b-value. Seven different IVIM-diffusion models with 4 to 14 parameters were implemented, which modeled diffusion and pseudo-diffusion as scalar or tensor quantities. The models were compared with respect to their fitting performance based on the goodness of fit (sum of squared fit residuals, chi² ) and their Akaike weights (calculated from the corrected Akaike information criterion). Lowest chi² values were found using the model with the largest number of model parameters. However, significantly highest Akaike weights indicating the most appropriate models for the acquired data were found with a nine-parameter IVIM-DTI model (with isotropic perfusion modeling) in normal-appearing white matter (NAWM), and with an 11-parameter model (IVIM-DTI with additional pseudo-diffusion anisotropy) in white matter with hyperintensities (WMH) and in gray matter (GM). The latter model allowed for the additional calculation of the fractional anisotropy of the pseudo-diffusion tensor (with a median value of 0.45 in NAWM, 0.23 in WMH, and 0.36 in GM), which is not accessible with the usually performed IVIM acquisitions based on three orthogonal diffusion-gradient directions.

Simultaneous Quantification of Anisotropic Microcirculation and Microstructure in Peripheral Nerve

Peripheral nerve injury is a significant public health challenge, and perfusion in the nerve is a potential biomarker for assessing the injury severity and prognostic outlook. Here, we applied a novel formalism that combined intravoxel incoherent motion (IVIM) and diffusion tensor imaging (DTI) to simultaneously characterize anisotropic microcirculation and microstructure in the rat sciatic nerve. Comparison to postmortem measurements revealed that the in vivo IVIM-DTI signal contained a fast compartment (2.32 ± 0.04 × 10−3 mm2/s mean diffusivity, mean ± sem, n = 6, paired t test p < 0.01) that could be attributed to microcirculation in addition to a slower compartment that had similar mean diffusivity as the postmortem nerve (1.04 ± 0.01 vs. 0.96 ± 0.05 × 10−3 mm2/s, p > 0.05). Although further investigation and technical improvement are warranted, this preliminary study demonstrates both the feasibility and potential for applying the IVIM-DTI methodology to peripheral nerves for quantifying perfusion in the presence of anisotropic tissue microstructure.

Intravoxel Incoherent Motion Magnetic Resonance Imaging in Skeletal Muscle: Review and Future Directions

Throughout the body, muscle structure and function can be interrogated using a variety of noninvasive magnetic resonance imaging (MRI) methods. Recently, intravoxel incoherent motion (IVIM) MRI has gained momentum as a method to evaluate components of blood flow and tissue diffusion simultaneously. Much of the prior research has focused on highly vascularized organs, including the brain, kidney, and liver. Unique aspects of skeletal muscle, including the relatively low perfusion at rest and large dynamic range of perfusion between resting and maximal hyperemic states, may influence the acquisition, postprocessing, and interpretation of IVIM data. Here, we introduce several of those unique features of skeletal muscle; review existing studies of IVIM in skeletal muscle at rest, in response to exercise, and in disease states; and consider possible confounds that should be addressed for muscle-specific evaluations. Most studies used segmented nonlinear least squares fitting with a b-value threshold of 200 sec/mm2 to obtain IVIM parameters of perfusion fraction (f), pseudo-diffusion coefficient (D*), and diffusion coefficient (D). In healthy individuals, across all muscles, the average ± standard deviation of D was 1.46 ± 0.30 × 10-3  mm2 /sec, D* was 29.7 ± 38.1 × 10-3  mm2 /sec, and f was 11.1 ± 6.7%. Comparisons of reported IVIM parameters in muscles of the back, thigh, and leg of healthy individuals showed no significant difference between anatomic locations. Throughout the body, exercise elicited a positive change of all IVIM parameters. Future directions including advanced postprocessing models and potential sequence modifications are discussed. LEVEL OF EVIDENCE: 2 TECHNICAL EFFICACY: Stage 2.

D* from diffusion MRI reveals a correspondence between ventricular cerebrospinal fluid volume and flow in the ischemic rodent model

Quantitative measurement of cerebrospinal fluid (CSF) flow and volume and longitudinal monitoring of CSF dynamics provide insights into the compensatory characteristics of post-stroke CSF. In this study, we compared the MRI pseudo-diffusion index (D*) of live and sacrificed rat brains to confirm the effect of ventricular CSF flow on diffusion signals. We observed the relationship between the CSF peak velocities and D* through Monte Carlo (MC) simulations to further understand the source of D* contrast. We also determined the dominant CSF flow using D* in three directions. Finally, we investigated the dynamic evolutions of ventricular CSF flow and volume in a stroke rat model (n = 8) from preoperative to up to 45 days after surgery and determined the correlation between ventricular CSF volume and flow. MC simulations showed a strong positive correlation between the CSF peak velocity and D* (r = 0.99). The dominant CSF flow variations in the 3D ventricle could be measured using the maximum D* map. A longitudinal positive correlation between ventricular CSF volume and D* was observed in the lateral (r = 0.74) and ventral-third (r = 0.81) ventricles, respectively. The directional D* measurements provide quantitative CSF volume and flow information, which would provide useful insights into ischemic stroke with diffusion MRI.

Crossed cerebellar diaschisis after acute ischemic stroke detected by intravoxel incoherent motion magnetic resonance imaging

PURPOSE

To study the value of 3.0 T magnetic resonance imaging with intravoxel incoherent motion (IVIM) in the diagnosis of the crossed cerebellar diaschisis (CCD) after the unilateral supratentorial acute ischemic stroke.

METHODS

Seventy-four patients with acute ischemic stroke who underwent intravoxel incoherent motion (IVIM), arterial spin labeling (ASL), and conventional magnetic resonance imaging (MRI) scanning were enrolled. Intravoxel incoherent motion-derived perfusion-related parameters including fast diffusion coefficient (D*), slow diffusion coefficient (D), vascular volume fraction (f), and arterial spin-labeling-derived cerebral blood flow (CBF) of bilateral cerebellum were measured.

RESULTS

In the CCD-positive group, D*, D, and CBF values of the contralateral cerebellum decreased compared with those of the ipsilesional cerebellum (P < 0.05), whereas f significantly increased (P < 0.05). A positive correlation was detected between the slow diffusion coefficient-based asymmetry index (AI-D) and the cerebral blood flow-based asymmetry index (AI-CBF) (r = 0.515, P < 0.01), whereas the vascular volume fraction-based asymmetry index (AI-f) had a negative correlation with the cerebral blood flow-based asymmetry index (AI-CBF) (r =  - 0.485, P < 0.01). Furthermore, the area under the receiver operating characteristic (ROC) curve value of AI-D and AI-f was 0.81 and 0.76, respectively.

CONCLUSIONS

The IVIM is feasible for the detection of CCD. This technique might provide opportunities to further investigate the pathophysiology of CCD.

Low b-value diffusion tensor imaging for measuring pseudorandom flow of cerebrospinal fluid

PURPOSE

Cerebrospinal fluid (CSF) plays an important role in the clearance system of the brain. Recently, low b-value diffusion tensor imaging (low-b DTI) has been reported to be useful in the observation of CSF flow; however, the precise flow property observed by low-b DTI has not been fully investigated. Accordingly, a mathematical framework of low-b DTI is proposed for investigating CSF and clarifying its pseudorandom flow.

THEORY

The framework will show that the limit of the diffusion tensor as b-value decreases to zero approximately represents the covariance of the velocity distribution of the CSF's pseudorandom flow.

METHODS

The low b-value diffusion tensor (DT_L ) of whole-brain CSF was obtained using diffusion-weighted echo-planar imaging. Seven healthy volunteers were scanned for intersubject analysis; three of the volunteers was consecutively scanned for repeatability analysis. Obtained DT_L was visually assessed by ellipsoid-representation map and was statistically evaluated by calculating mean diffusivity (MD) and fractional anisotropy (FA) in regions of interest (ROIs) representing intensive pseudorandom flow.

RESULTS

Obtained DT_L consistently shows large and anisotropic diffusivity in some segments of CSF, typically the ROIs around the foramen of Monro, the aqueduct, the prepontine cistern, the middle cerebral artery, and the Sylvian fissure throughout the study. The statistical analysis shows high repeatability and consistently high MD and FA in all the ROIs for all the volunteers.

CONCLUSION

From the viewpoint of the proposed framework, the high and anisotropic DT_L in the ROIs indicates large covariance of velocity distribution, which represents intensive pseudorandom flows of CSF.

A Bayesian approach to diffusional kurtosis imaging

PURPOSE

Diffusional kurtosis metrics show high performance for detecting pathological changes and are therefore expected to be disease biomarkers. Kurtosis maps, however, tend to be noisy. The maps' visual quality is crucial for disease diagnosis, even when kurtosis is being used quantitatively. A Bayesian method was proposed to curtail the large statistical error inherent in kurtosis estimation while maintaining potential application to biomarkers.

THEORY

Gaussian priors are determined from first-step estimations implemented using the least-square method (LSM). The likelihood-function variance is determined from the residuals of the estimation. Although the proposed approach is similar to a regularized LSM, regularization parameters do not have to be artificially adjusted. An appropriate balance between denoising and preventing false shrinkages of metric dispersions is automatically achieved.

METHODS

Map qualities achieved using the conventional and proposed methods were compared. The receiver-operating characteristic analysis was performed for glioma-grade differentiation using simulated low- and high-grade glioma DWI datasets. Noninferiority of the proposed method was tested for areas under the curves (AUCs).

RESULTS

The noisier the conventional maps, the better the proposed Bayesian method improved them. Noninferiority of the proposed method was confirmed by AUC tests for all kurtosis-related metrics. Superiority of the proposed method was also established for several metrics.

CONCLUSIONS

The proposed approach improved noisy kurtosis maps while maintaining their performances as biomarkers without increasing data acquisition requirements or arbitrarily choosing LSM regularization parameters. This approach may enable the use of higher-order terms in diffusional kurtosis imaging (DKI) fitting functions by suppressing overfitting, thereby improving the DKI-estimation accuracy.

Investigation of intravoxel incoherent motion tensor imaging for the characterization of the in vivo human heart

PURPOSE

To investigate intravoxel incoherent motion (IVIM) tensor imaging of the in vivo human heart and elucidate whether the estimation of IVIM tensors is affected by the complexity of pseudo-diffusion components in myocardium.

METHODS

The cardiac IVIM data of 10 healthy subjects were acquired using a diffusion weighted spin-echo echo-planar imaging sequence along 6 gradient directions with 10 b values (0~400 s/mm² ). The IVIM data of left ventricle myocardium were fitted to the IVIM tensor model. The complexity of myocardial pseudo-diffusion components was reduced through exclusion of low b values (0 and 5 s/mm² ) from the IVIM curve-fitting analysis. The fractional anisotropy, mean fraction/mean diffusivity, and Westin measurements of pseudo-diffusion tensors (f_p and D*) and self-diffusion tensor (D), as well as the angle between the main eigenvector of f_p (or D*) and that of D, were computed and compared before and after excluding low b values.

RESULTS

The fractional anisotropy values of f_p and D* without low b value participation were significantly higher (P < .001) than those with low b value participation, but an opposite trend was found for the mean fraction/diffusivity values. Besides, after removing low b values, the angle between the main eigenvector of f_p (or D*) and that of D became small, and both f_p and D* tensors presented significant decrease of spherical components and significant increase of linear components.

CONCLUSION

The presence of multiple pseudo-diffusion components in myocardium indeed influences the estimation of IVIM tensors. The IVIM tensor model needs to be further improved to account for the complexity of myocardial microcirculatory network and blood flow.

Pattern recognition analysis of directional intravoxel incoherent motion MRI in ischemic rodent brains

This study aimed to demonstrate a reliable automatic segmentation method for independently separating reduced diffusion and decreased perfusion areas in ischemic stroke brains using constrained nonnegative matrix factorization (cNMF) pattern recognition in directional intravoxel incoherent motion MRI (IVIM-MRI). First, the feasibility of cNMF-based segmentation of IVIM signals was investigated in both simulations and in vivo experiments. The cNMF analysis was independently performed for S₀ -normalized and scaled (by the difference between the maximum and minimum) IVIM signals, respectively. Segmentations of reduced diffusion (from S₀ -normalized IVIM signals) and decreased perfusion (from scaled IVIM signals) areas were performed using the corresponding cNMF pattern weight maps. Second, Monte Carlo simulations were performed for directional IVIM signals to investigate the relationship between the degree of vessel alignment and the direction of the diffusion gradient. Third, directional IVIM-MRI experiments (x, y and z diffusion-gradient directions, 20 b values at 7 T) were performed for normal (n = 4), sacrificed (n = 1, no flow) and ischemic stroke models (n = 4, locally reduced flow). The results showed that automatic segmentation of the hypoperfused lesion using cNMF analysis was more accurate than segmentation using the conventional double-exponential fitting. Consistent with the simulation, the double-exponential pattern of the IVIM signals was particularly strong in white matter and ventricle regions when the directional flows were aligned with the applied diffusion-gradient directions. cNMF analysis of directional IVIM signals allowed robust automated segmentation of white matter, ventricle, vascular and hypoperfused regions in the ischemic brain. In conclusion, directional IVIM signals were simultaneously sensitive to diffusion and aligned flow and were particularly useful for automatically segmenting ischemic lesions via cNMF-based pattern recognition.

Quantitative histomorphometry of capillary microstructure in deep white matter

White matter lesions represent a major risk factor for dementia in elderly people. Magnetic Resonance Imaging (MRI) studies have demonstrated cerebral blood flow reduction in age-related white matter lesions, indicating that vascular alterations are involved in developing white matter lesions. Hypoperfusion and changes in capillary morphology are generally linked to dementia. However, a quantitative study describing these microvascular alterations in white matter lesions is missing in the literature; most previous microvascular studies being on the cortex. The aim of this work is to identify and quantify capillary microstructural changes involved in the appearance of deep subcortical lesions (DSCL). We characterize the distribution of capillary diameter, thickness, and density in the deep white matter in a population of 75 elderly subjects, stratified into three equal groups according to DSCL: Control (subject without DSCL), Lesion (sample presenting DSCL), and Normal Appearing White Matter (NAWM, the subject presented DSCL but not at the sampled tissue location). Tissue samples were selected from the Cognitive Function and Aging Study (CFAS), a cohort representative of an aging population, from which immunohistochemically-labeled histological images were acquired. To accurately estimate capillary diameters and thicknesses from the 2D histological images, we also introduce a novel semi-automatic method robust to non-perpendicular incidence angle of capillaries into the imaging plane, and to non-circular deformations of capillary cross sections. Subjects with DSCL presented a significant increase in capillary wall thickness, a decrease in the diameter intra-subject variability (but not in the mean), and a decrease in capillary density. No significant difference was observed between controls and NAWM. Both capillary wall thickening and reduction in capillary density contribute to the reduction of cerebral blood flow previously reported for white matter lesions. The obtained distributions provide reliable statistics of capillary microstructure useful to inform the modeling of human cerebral blood flow, for instance to define microcirculation models for their estimation from MRI or to perform realistic cerebral blood flow simulations.

Enhancing intravoxel incoherent motion parameter mapping in the brain using k-b PCA

Intravoxel incoherent motion (IVIM) imaging of diffusion and perfusion parameters in the brain using parallel imaging suffers from local noise amplification. To address the issue, signal correlations in space and along the diffusion encoding dimension are exploited jointly using a constrained image reconstruction approach. IVIM imaging was performed on a clinical 3 T MR system with diffusion weighting along six gradient directions and 16 b-values encoded per direction across a range of 0-900 s/mm² . Data were collected in 11 subjects, retrospectively undersampled in k-space with net factors ranging from 2 to 6 and reconstructed using CG-SENSE and the proposed k-b PCA approach. Results of k-b PCA and CG-SENSE from retrospectively undersampled data were compared with those from the fully sampled reference. In addition, prospective single-shot k-b undersampling was implemented and data were acquired in five additional volunteers. IVIM parameter maps were derived using a segmented least-squares method. The proposed k-b PCA method outperformed CG-SENSE in terms of reconstruction errors for effective undersampling factors of 3 and beyond. Undersampling artifacts were effectively removed with k-b PCA up to sixfold undersampling. At net sixfold undersampling, relative errors (compared with the fully sampled reference) of image magnitude and IVIM parameters (D, f and D^* ) were (median ± interquartile range): 3.5 ± 3.7 versus 25.3 ± 25.8%, 2.7 ± 3.6 versus 14.2 ± 20.4%, 15.1 ± 26.1 versus 96.6 ± 67.4% and 14.8 ± 26.6 versus 100 ± 195.1% for k-b PCA versus CG-SENSE, respectively. Acquisition with sixfold prospective undersampling yielded average IVIM parameters in the brain of 0.79 ± 0.18 × 10^-3  mm² /s for D, 7.35 ± 7.27% for f and 7.11 ± 2.39 × 10^-3  mm² /s for D^* . Constrained reconstruction using k-b PCA improves IVIM parameter mapping from undersampled data when compared with CG-SENSE reconstruction. Prospectively undersampled single-shot echo planar imaging acquisition was successfully employed using k-b PCA, demonstrating a reduction of image artifacts and noise relative to parallel imaging.