Quantification of microcirculatory parameters by joint analysis of flow-compensated and non-flow-compensated intravoxel incoherent motion (IVIM) data

Comparison of methods for intravoxel incoherent motion parameter estimation in the brain from flow-compensated and non-flow-compensated diffusion-encoded data

PURPOSE

Joint analysis of flow-compensated (FC) and non-flow-compensated (NC) diffusion MRI (dMRI) data has been suggested for increased robustness of intravoxel incoherent motion (IVIM) parameter estimation. For this purpose, a set of methods commonly used or previously found useful for IVIM analysis of dMRI data obtained with conventional diffusion encoding were evaluated in healthy human brain.

METHODS

Five methods for joint IVIM analysis of FC and NC dMRI data were compared: (1) direct non-linear least squares fitting, (2) a segmented fitting algorithm with estimation of the diffusion coefficient from higher b-values of NC data, (3) a Bayesian algorithm with uniform prior distributions, (4) a Bayesian algorithm with spatial prior distributions, and (5) a deep learning-based algorithm. Methods were evaluated on brain dMRI data from healthy subjects and simulated data at multiple noise levels. Bipolar diffusion encoding gradients were used with b-values 0-200 s/mm2 and corresponding flow weighting factors 0-2.35 s/mm for NC data and by design 0 for FC data. Data were acquired twice for repeatability analysis.

RESULTS

Measurement repeatability as well as estimation bias and variability were at similar levels or better with the Bayesian algorithm with spatial prior distributions and the deep learning-based algorithm for IVIM parametersD $$ D $$ andf $$ f $$ , and for the Bayesian algorithm only forv d $$ {v}_d $$ , relative to the other methods.

CONCLUSION

A Bayesian algorithm with spatial prior distributions is preferable for joint IVIM analysis of FC and NC dMRI data in the healthy human brain, but deep learning-based algorithms appear promising.

Spatial profiling of in vivo diffusion-weighted MRI parameters in the healthy human kidney

OBJECTIVE

Diffusion-weighted MRI is a technique that can infer microstructural and microcirculatory features from biological tissue, with particular application to renal tissue. There is extensive literature on diffusion tensor imaging (DTI) of anisotropy in the renal medulla, intravoxel incoherent motion (IVIM) measurements separating microstructural from microcirculation effects, and combinations of the two. However, interpretation of these features and adaptation of more specific models remains an ongoing challenge. One input to this process is a whole organ distillation of corticomedullary contrast of diffusion metrics, as has been explored for other renal biomarkers.

MATERIALS AND METHODS

In this work, we probe the spatial dependence of diffusion MRI metrics with concentrically layered segmentation in 11 healthy kidneys at 3 T. The metrics include those from DTI, IVIM, a combined approach titled "REnal Flow and Microstructure AnisotroPy (REFMAP)", and a multiply encoded model titled "FC-IVIM" providing estimates of fluid velocity and branching length.

RESULTS

Fractional anisotropy decreased from the inner kidney to the outer kidney with the strongest layer correlation in both parenchyma (including cortex and medulla) and medulla with Spearman correlation coefficients and p-values (r, p) of (0.42, <0.001) and (0.37, <0.001), respectively. Also, dynamic parameters derived from the three models significantly decreased with a high correlation from the inner to the outer parenchyma or medulla with (r, p) ranges of (0.46-0.55, <0.001).

CONCLUSIONS

These spatial trends might find implications for indirect assessments of kidney physiology and microstructure using diffusion MRI.

Velocity-compensated intravoxel incoherent motion of the human calf muscle

PURPOSE

To determine whether intravoxel incoherent motion (IVIM) describes the blood perfusion in muscles better, assuming pseudo diffusion (Bihan Model 1) or ballistic motion (Bihan Model 2).

METHODS

IVIM parameters were measured in 18 healthy subjects with three different diffusion gradient time profiles (bipolar with two diffusion times and one with velocity compensation) and 17 b-values (0-600 s/mm2) at rest and after muscle activation. The diffusion coefficient, perfusion fraction, and pseudo-diffusion coefficient were estimated with a segmented fit in the gastrocnemius medialis (GM) and tibialis anterior (TA) muscles.

RESULTS

Velocity-compensated gradients resulted in a decreased perfusion fraction (6.9% ± 1.4% vs. 4.4% ± 1.3% in the GM after activation) and pseudo-diffusion coefficient (0.069 ± 0.046 mm2/s vs. 0.014 ± 0.006 in the GM after activation) compared to the bipolar gradients with the longer diffusion encoding time. Increased diffusion coefficients, perfusion fractions, and pseudo-diffusion coefficients were observed in the GM after activation for all gradient profiles. However, the increase was significantly smaller for the velocity-compensated gradients. A diffusion time dependence was found for the pseudo-diffusion coefficient in the activated muscle.

CONCLUSION

Velocity-compensated diffusion gradients significantly suppress the IVIM effect in the calf muscle, indicating that the ballistic limit is mostly reached, which is supported by the time dependence of the pseudo-diffusion coefficient.

Advanced Diffusion-Weighted MRI for Cancer Microstructure Assessment in Body Imaging, and Its Relationship With Histology

Diffusion-weighted magnetic resonance imaging (DW-MRI) aims to disentangle multiple biological signal sources in each imaging voxel, enabling the computation of innovative maps of tissue microstructure. DW-MRI model development has been dominated by brain applications. More recently, advanced methods with high fidelity to histology are gaining momentum in other contexts, for example, in oncological applications of body imaging, where new biomarkers are urgently needed. The objective of this article is to review the state-of-the-art of DW-MRI in body imaging (ie, not including the nervous system) in oncology, and to analyze its value as compared to reference colocalized histology measurements, given that demonstrating the histological validity of any new DW-MRI method is essential. In this article, we review the current landscape of DW-MRI techniques that extend standard apparent diffusion coefficient (ADC), describing their acquisition protocols, signal models, fitting settings, microstructural parameters, and relationship with histology. Preclinical, clinical, and in/ex vivo studies were included. The most used techniques were intravoxel incoherent motion (IVIM; 36.3% of used techniques), diffusion kurtosis imaging (DKI; 16.7%), vascular, extracellular, and restricted diffusion for cytometry in tumors (VERDICT; 13.3%), and imaging microstructural parameters using limited spectrally edited diffusion (IMPULSED; 11.7%). Another notable category of techniques relates to innovative b-tensor diffusion encoding or joint diffusion-relaxometry. The reviewed approaches provide histologically meaningful indices of cancer microstructure (eg, vascularization/cellularity) which, while not necessarily accurate numerically, may still provide useful sensitivity to microscopic pathological processes. Future work of the community should focus on improving the inter-/intra-scanner robustness, and on assessing histological validity in broader contexts. LEVEL OF EVIDENCE: NA TECHNICAL EFFICACY: Stage 2.

Low b-value Diffusion Tensor Imaging to Analyze the Dynamics of Cerebrospinal Fluid: Resolving Intravoxel Pseudorandom Motion into Ordered and Disordered Motions

PURPOSE

Analysis of cerebrospinal fluid (CSF) dynamics may be beneficial for understanding the mechanisms and diagnosis of several neurological diseases. Low b-value diffusion tensor imaging (low-b DTI) is useful for observing the slow and complex motion of the CSF. Theoretically, a mathematical framework suggests that low-b DTI provides the variance of the pseudorandom motion of the CSF. Furthermore, low-b DTI could provide comprehensive information on fluid dynamics. Accordingly, we proposed an analysis technique that resolves intravoxel pseudorandom motion into ordered (linear) and disordered (random) motions based on the mathematical framework.

METHODS

The proposed analysis technique helps measure low-b DTI with multiple diffusion times and linearly fits its mean diffusivity (MD) with the diffusion time to obtain two parameters, double-slope Vv and y-intersect Dr, which represent the variance of the velocity distribution of linear motion and the diffusion coefficient of random motion, respectively. Seven healthy subjects were scanned to evaluate the proposed technique and investigate fluid dynamics in several representative ROIs.

RESULTS

The obtained data showed the validity of the technique, repeatability, and consistency across the subjects in ROIs, such as the lateral ventricle (LV), third ventricle (3V), fourth ventricle (4V), and Sylvian fissure (SF). The obtained parameters Vv and Dr highlighted different characteristics of fluid dynamics in the representative ROIs: low Vv and low Dr in the LV, high Vv and moderate Dr in the 3V, and moderate Vv and moderate Dr in the 4V and SF.

CONCLUSION

The proposed analysis technique will facilitate a comprehensive investigation of the complex dynamics of the CSF using resolved parameters representing ordered and disordered motions.

Extra-axonal contribution to double diffusion encoding-based pore size estimates in the corticospinal tract

OBJECTIVE

To study the origin of compartment size overestimation in double diffusion encoding MRI (DDE) in vivo experiments in the human corticospinal tract. Here, the extracellular space is hypothesized to be the origin of the DDE signal. By exploiting the DDE sensitivity to pore shape, it could be possible to identify the origin of the measured signal. The signal difference between parallel and perpendicular diffusion gradient orientation can indicate if a compartment is regular or eccentric in shape. As extracellular space can be considered an eccentric compartment, a positive difference would mean a high contribution to the compartment size estimates.

MATERIALS AND METHODS

Computer simulations using MISST and in vivo experiments in eight healthy volunteers were performed. DDE experiments using a double spin-echo preparation with eight perpendicular directions were measured in vivo. The difference between parallel and perpendicular gradient orientations was analyzed using a Wilcoxon signed-rank test and a Mann-Whitney U test.

RESULTS

Simulations and MR experiments showed a statistically significant difference between parallel and perpendicular diffusion gradient orientation signals ([Formula: see text]).

CONCLUSION

The results suggest that the DDE-based size estimate may be considerably influenced by the extra-axonal compartment. However, the experimental results are also consistent with purely intra-axonal contributions in combination with a large fiber orientation dispersion.

Investigation of time-dependent diffusion in extra-axial brain tumors using oscillating-gradient spin-echo

Oscillating gradient spin-echo (OGSE) sequences provide access to short diffusion times and may provide insight into micro-scale internal structures of pathologic lesions based on an analysis of changes in diffusivity with differing diffusion times. We hypothesized that changes in diffusivity acquired with a shorter diffusion time may permit elucidation of properties related to the internal structure of extra-axial brain tumors. This study aimed to investigate the utility of changes in diffusivity between short and long diffusion times for characterizing extra-axial brain tumors. In total, 12 patients with meningothelial meningiomas, 13 patients with acoustic neuromas, and 11 patients with pituitary adenomas were scanned with a 3 T magnetic resonance imaging (MRI) scanner with diffusion-weighted imaging (DWI) using OGSE and pulsed gradient spin-echo (PGSE) (effective diffusion times [Δ_eff]: 6.5 ms and 35.2 ms) with b-values of 0 and 1000 s/mm². Relative percentage changes between shorter and longer diffusion times were calculated using region-of-interest (ROI) analysis of brain tumors on λ₁, λ₂, λ₃, and mean diffusivity (MD) maps. The diffusivities of PGSE, OGSE, and relative percentage changes were compared among each tumor type using a multiple comparisons Steel-Dwass test. The mean (standard deviation) MD at Δ_eff of 6.5 ms was 1.07 ± 0.23 10^-3 mm²/s, 1.19 ± 0.18 10^-3 mm²/s, 1.19 ± 0.21 10^-3 mm²/s for meningothelial meningiomas, acoustic neuromas, and pituitary adenomas, respectively. The mean (standard deviation) MD at Δ_eff of 35.2 ms was 0.93 ± 0.22 10^-3 mm²/s, 1.07 ± 0.19 10^-3 mm²/s, 0.82 ± 0.21 10^-3 mm²/s for meningothelial meningiomas, acoustic neuromas, and pituitary adenomas, respectively. The mean (standard deviation) of the relative percentage change was 15.7 ± 4.4%, 12.4 ± 8.2%, 46.8 ± 11.3% for meningothelial meningiomas, acoustic neuromas, and pituitary adenomas, respectively. Compared to meningiomas and acoustic neuromas, pituitary adenoma exhibited stronger diffusion time-dependence with diffusion times between 6.5 ms and 35.2 ms (P < 0.05). In conclusion, differences in diffusion time-dependence may be attributed to differences in the internal structures of brain tumors. DWI with a short diffusion time may provide additional information on the microstructure of each tumor and contribute to tumor diagnosis.

Metabolic activity diffusion imaging (MADI): II. Noninvasive, high-resolution human brain mapping of sodium pump flux and cell metrics

We introduce a new ¹ H₂ O magnetic resonance approach: metabolic activity diffusion imaging (MADI). Numerical diffusion-weighted imaging decay simulations characterized by the mean cellular water efflux (unidirectional) rate constant (k_io ), mean cell volume (V), and cell number density (ρ) are produced from Monte Carlo random walks in virtual stochastically sized/shaped cell ensembles. Because of active steady-state trans-membrane water cycling (AWC), k_io reflects the cytolemmal Na⁺ , K⁺ ATPase (NKA) homeostatic cellular metabolic rate (^c MR_NKA ). A digital 3D "library" contains thousands of simulated single diffusion-encoded (SDE) decays. Library entries match well with disparate, animal, and human experimental SDE decays. The V and ρ values are consistent with estimates from pertinent in vitro cytometric and ex vivo histopathological literature: in vivo V and ρ values were previously unavailable. The library allows noniterative pixel-by-pixel experimental SDE decay library matchings that can be used to advantage. They yield proof-of-concept MADI parametric mappings of the awake, resting human brain. These reflect the tissue morphology seen in conventional MRI. While V is larger in gray matter (GM) than in white matter (WM), the reverse is true for ρ. Many brain structures have k_io values too large for current, invasive methods. For example, the median WM k_io is 22s^-1 ; likely reflecting mostly exchange within myelin. The k_io •V product map displays brain tissue ^c MR_NKA variation. The GM activity correlates, quantitatively and qualitatively, with the analogous resting-state brain ¹⁸ FDG-PET tissue glucose consumption rate (^t MR_glucose ) map; but noninvasively, with higher spatial resolution, and no pharmacokinetic requirement. The cortex, thalamus, putamen, and caudate exhibit elevated metabolic activity. MADI accuracy and precision are assessed. The results are contextualized with literature overall homeostatic brain glucose consumption and ATP production/consumption measures. The MADI/PET results suggest different GM and WM metabolic pathways. Preliminary human prostate results are also presented.

Metabolic activity diffusion imaging (MADI): I. Metabolic, cytometric modeling and simulations

Evidence mounts that the steady-state cellular water efflux (unidirectional) first-order rate constant (k_io [s^-1 ]) magnitude reflects the ongoing, cellular metabolic rate of the cytolemmal Na⁺ , K⁺ -ATPase (NKA), ^c MR_NKA (pmol [ATP consumed by NKA]/s/cell), perhaps biology's most vital enzyme. Optimal ¹ H₂ O MR k_io determinations require paramagnetic contrast agents (CAs) in model systems. However, results suggest that the homeostatic metabolic k_io biomarker magnitude in vivo is often too large to be reached with allowable or possible CA living tissue distributions. Thus, we seek a noninvasive (CA-free) method to determine k_io in vivo. Because membrane water permeability has long been considered important in tissue water diffusion, we turn to the well-known diffusion-weighted MRI (DWI) modality. To analyze the diffusion tensor magnitude, we use a parsimoniously primitive model featuring Monte Carlo simulations of water diffusion in virtual ensembles comprising water-filled and -immersed randomly sized/shaped contracted Voronoi cells. We find this requires two additional, cytometric properties: the mean cell volume (V [pL]) and the cell number density (ρ [cells/μL]), important biomarkers in their own right. We call this approach metabolic activity diffusion imaging (MADI). We simulate water molecule displacements and transverse MR signal decays covering the entirety of b-space from pure water (ρ = V = 0; k_io undefined; diffusion coefficient, D₀ ) to zero diffusion. The MADI model confirms that, in compartmented spaces with semipermeable boundaries, diffusion cannot be described as Gaussian: the nanoscopic D (D_n ) is diffusion time-dependent, a manifestation of the "diffusion dispersion". When the "well-mixed" (steady-state) condition is reached, diffusion becomes limited, mainly by the probabilities of (1) encountering (ρ, V), and (2) permeating (k_io ) cytoplasmic membranes, and less so by D_n magnitudes. Importantly, for spaces with large area/volume (A/V; claustrophobia) ratios, this can happen in less than a millisecond. The model matches literature experimental data well, with implications for DWI interpretations.

Precision of region of interest-based tri-exponential intravoxel incoherent motion quantification and the role of the Intervoxel spatial distribution of flow velocities

PURPOSE

The purpose of this work was to obtain precise tri-exponential intravoxel incoherent motion (IVIM) quantification in the liver using 2D (b-value and first-order motion moment [M1 ]) IVIM-DWI acquisitions and region of interest (ROI)-based fitting techniques.

METHODS

Diffusion MRI of the liver was performed in 10 healthy volunteers using three IVIM-DWI acquisitions: conventional monopolar, optimized monopolar, and optimized 2D (b-M1 ). For each acquisition, bi-exponential and tri-exponential full, segmented, and over-segmented ROI-based fitting and a newly proposed blood velocity SDdistribution (BVD) fitting technique were performed to obtain IVIM estimates in the right and left liver lobes. Fitting quality was evaluated using corrected Akaike information criterion. Precision metrics (test-retest repeatability, inter-reader reproducibility, and inter-lobar agreement) were evaluated using Bland-Altman analysis, repeatability/reproducibility coefficients (RPCs), and paired sample t-tests. Precision was compared across acquisitions and fitting methods.

RESULTS

High repeatability and reproducibility was observed in the estimations of the diffusion coefficient (Dtri  = [1.03 ± 0.11] × 10-3  mm2 /s; RPCs ≤ 1.34 × 10-4  mm2 /s), perfusion fractions (F1  = 3.19 ± 1.89% and F2  = 16.4 ± 2.07%; RPCs ≤ 2.51%), and blood velocity SDs (Vb,1  = 1.44 ± 0.14 mm/s and Vb,2  = 3.62 ± 0.13 mm/s; RPCs ≤ 0.41 mm/s) in the right liver lobe using the 2D (b-M1 ) acquisition in conjunction with BVD fitting. Using these methods, significantly larger (p < 0.01) estimates of Dtri and F1 were observed in the left lobe in comparison to the right lobe, while estimates of Vb,1 and Vb,2 demonstrated high interlobar agreement (RPCs ≤ 0.45 mm/s).

CONCLUSIONS

The 2D (b-M1 ) IVIM-DWI data acquisition in conjunction with BVD fitting enables highly precise tri-exponential IVIM quantification in the right liver lobe.

Reduction of the cardiac pulsation artifact and improvement of lesion conspicuity in flow‐compensated diffusion images in the liver—A quantitative evaluation of postprocessing algorithms

PURPOSE

To enhance image quality of flow-compensated diffusion-weighted liver MRI data by increasing the lesion conspicuity and reducing the cardiac pulsation artifact using postprocessing algorithms.

METHODS

Diffusion-weighted image data of 40 patients with liver lesions had been acquired at 1.5 T. These data were postprocessed with 5 different algorithms (weighted averaging, p-mean, percentile, outlier exclusion, and exception set). Four image properties of the postprocessed data were evaluated for optimizing the algorithm parameters. These properties were the lesion to tissue contrast-to-noise ratio (CNR), the reduction of the cardiac pulsation artifact, the data consistency, and the vessel darkness. They were combined into a total quality score ( Q total , $$ {Q}_{\mathrm{total}}, $$ set to 1 for the trace-weighted reference image), which was used to rate the image quality objectively.

RESULTS

The weighted averaging algorithm performed best according to the total quality score ( Q total = 1.111 ± 0.067 $$ {Q}_{\mathrm{total}}=1.111\pm 0.067 $$ ). The further ranking was outlier exclusion algorithm ( Q total = 1.086 ± 0.061 $$ {Q}_{\mathrm{total}}=1.086\pm 0.061 $$ ), p-mean algorithm ( Q total = 1.045 ± 0.049 $$ {Q}_{\mathrm{total}}=1.045\pm 0.049 $$ ), percentile algorithm ( Q total = 1.012 ± 0.049 $$ {Q}_{\mathrm{total}}=1.012\pm 0.049 $$ ), and exception set algorithm ( Q total = 0.957 ± 0.027 $$ {Q}_{\mathrm{total}}=0.957\pm 0.027 $$ ). All optimized algorithms except for the exception set algorithm corrected the pulsation artifact and increased the lesion CNR. Changes in Q total $$ {Q}_{\mathrm{total}} $$ were significant for all optimized algorithms except for the percentile algorithm. Liver ADC was significantly reduced (except for the exception set algorithm), particularly in the left lobe.

CONCLUSION

Postprocessing algorithms should be used for flow-compensated liver DWI. The proposed weighted averaging algorithm seems to be suited best to increase the image quality of artifact-corrupted flow-compensated diffusion-weighted liver data.

VERDICT MRI for radiation treatment response assessment in neuroendocrine tumors

Noninvasive methods to study changes in tumor microstructure enable early assessment of treatment response and thus facilitate personalized treatment. The aim of this study was to evaluate the diffusion MRI model, Vascular, Extracellular and Restricted Diffusion for Cytometry in Tumors (VERDICT), for early response assessment to external radiation treatment and to compare the results with those of more studied sets of parameters derived from diffusion-weighted MRI data. Mice xenografted with human small intestine tumors were treated with external radiation treatment, and diffusion MRI experiments were performed on the day before and up to 2 weeks after treatment. The diffusion models VERDICT, ADC, IVIM, and DKI were fitted to MRI data, and the treatment response of each tumor was calculated based on pretreatment tumor growth and post-treatment tumor volume regression. Linear regression and correlation analysis were used to evaluate each model and their respective parameters for explaining the treatment response. VERDICT analysis showed significant changes from day -1 to day 3 for the intracellular and extracellular volume fraction, as well as the cell radius index (p < 0.05; Wilcoxon signed-rank test). The strongest correlation between the diffusion model parameters and the tumor treatment response was seen for the ADC, kurtosis-corrected diffusion coefficient, and intracellular volume fraction on day 3 (τ = 0.47, 0.52, and -0.49, respectively, p < 0.05; Kendall rank correlation coefficient). Of all the tested models, VERDICT held the strongest explanatory value for the tumor treatment response on day 3 (R²  = 0.75, p < 0.01; linear regression). In conclusion, VERDICT has potential for early assessment of external radiation treatment and may provide further insights into the underlying biological effects of radiation on tumor tissue. In addition, the results suggest that the time window for assessment of treatment response using dMRI may be narrow.

Flow-compensated diffusion encoding in MRI for improved liver metastasis detection

Magnetic resonance (MR) diffusion-weighted imaging (DWI) is often used to detect focal liver lesions (FLLs), though DWI image quality can be limited in the left liver lobe owing to the pulsatile motion of the nearby heart. Flow-compensated (FloCo) diffusion encoding has been shown to reduce this pulsation artifact. The purpose of this prospective study was to intra-individually compare DWI of the liver acquired with conventional monopolar and FloCo diffusion encoding for assessing metastatic FLLs in non-cirrhotic patients. Forty patients with known or suspected multiple metastatic FLLs were included and measured at 1.5 T field strength with a conventional (monopolar) and a FloCo diffusion encoding EPI sequence (single refocused; b-values, 50 and 800 s/mm²). Two board-certified radiologists analyzed the DWI images independently. They issued Likert-scale ratings (1 = worst, 5 = best) for pulsation artifact severity and counted the difference of lesions visible at b = 800 s/mm² separately for small and large FLLs (i.e., < 1 cm or > 1 cm) and separately for left and right liver lobe. Differences between the two diffusion encodings were assessed with the Wilcoxon signed-rank test. Both readers found a reduction in pulsation artifact in the liver with FloCo encoding (p < 0.001 for both liver lobes). More small lesions were detected with FloCo diffusion encoding in both liver lobes (left lobe: six and seven additional lesions by readers 1 and 2, respectively; right lobe: five and seven additional lesions for readers 1 and 2, respectively). Both readers found one additional large lesion in the left liver lobe. Thus, flow-compensated diffusion encoding appears more effective than monopolar diffusion encoding for the detection of liver metastases.

b value and first-order motion moment optimized data acquisition for repeatable quantitative intravoxel incoherent motion DWI

PURPOSE

To design a b value and first-order motion moment (M₁ ) optimized data acquisition for repeatable intravoxel incoherent motion (IVIM) quantification in the liver.

METHODS

Cramer-Rao lower bound optimization was performed to determine optimal monopolar and optimal 2D samplings of the b-M₁ space based on noise performance. Monte Carlo simulations were used to evaluate the bias and variability in estimates obtained using the proposed optimal samplings and conventional monopolar sampling. Diffusion MRI of the liver was performed in 10 volunteers using 3 IVIM acquisitions: conventional monopolar, optimized monopolar, and b-M₁ -optimized gradient waveforms (designed based on the optimal 2D sampling). IVIM parameter maps of diffusion coefficient, perfusion fraction, and blood velocity SD were obtained using nonlinear least squares fitting. Noise performance (SDs), stability (outlier percentage), and test-retest or scan-rescan repeatability (intraclass correlation coefficients) were evaluated and compared across acquisitions.

RESULTS

Cramer-Rao lower bound and Monte Carlo simulations demonstrated improved noise performance of the optimal 2D sampling in comparison to monopolar samplings. Evaluating the designed b-M₁ -optimized waveforms in healthy volunteers, significant decreases (p < 0.05) in the SDs and outlier percentages were observed for measurements of diffusion coefficient, perfusion fraction, and blood velocity SD in comparison to measurements obtained using monopolar samplings. Good-to-excellent repeatability (intraclass correlation coefficients ≥ 0.77) was observed for all 3 parameters in both the right and left liver lobes using the b-M₁ -optimized waveforms.

CONCLUSIONS

2D b-M₁ -optimized data acquisition enables repeatable IVIM quantification with improved noise performance. 2D acquisitions may advance the establishment of IVIM quantitative biomarkers for liver diseases.

Evaluating diffusion dispersion across an extended range of b-values and frequencies: Exploiting gap-filled OGSE shapes, strong gradients, and spiral readouts

Purpose

To address the long echo times and relatively weak diffusion sensitization that typically limit oscillating gradient spin‐echo (OGSE) experiments, an OGSE implementation combining spiral readouts, gap‐filled oscillating gradient shapes providing stronger diffusion encoding, and a high‐performance gradient system is developed here and utilized to investigate the tradeoff between b‐value and maximum OGSE frequency in measurements of diffusion dispersion (i.e., the frequency dependence of diffusivity) in the in vivo human brain. In addition, to assess the effects of the marginal flow sensitivity introduced by these OGSE waveforms, flow‐compensated variants are devised for experimental comparison.

Methods

Using DTI sequences, OGSE acquisitions were performed on three volunteers at b‐values of 300, 500, and 1000 s/mm² and frequencies up to 125, 100, and 75 Hz, respectively; scans were performed for gap‐filled oscillating gradient shapes with and without flow sensitivity. Pulsed gradient spin‐echo DTI acquisitions were also performed at each b‐value. Upon reconstruction, mean diffusivity (MD) maps and maps of the diffusion dispersion rate were computed.

Results

The power law diffusion dispersion model was found to fit best to MD measurements acquired at b = 1000 s/mm² despite the associated reduction of the spectral range; this observation was consistent with Monte Carlo simulations. Furthermore, diffusion dispersion rates without flow sensitivity were slightly higher than flow‐sensitive measurements.

Conclusion

The presented OGSE implementation provided an improved depiction of diffusion dispersion and demonstrated the advantages of measuring dispersion at higher b‐values rather than higher frequencies within the regimes employed in this study.

On the generalizability of diffusion MRI signal representations across acquisition parameters, sequences and tissue types: Chronicles of the MEMENTO challenge

Diffusion MRI (dMRI) has become an invaluable tool to assess the microstructural organization of brain tissue. Depending on the specific acquisition settings, the dMRI signal encodes specific properties of the underlying diffusion process. In the last two decades, several signal representations have been proposed to fit the dMRI signal and decode such properties. Most methods, however, are tested and developed on a limited amount of data, and their applicability to other acquisition schemes remains unknown. With this work, we aimed to shed light on the generalizability of existing dMRI signal representations to different diffusion encoding parameters and brain tissue types. To this end, we organized a community challenge - named MEMENTO, making available the same datasets for fair comparisons across algorithms and techniques. We considered two state-of-the-art diffusion datasets, including single-diffusion-encoding (SDE) spin-echo data from a human brain with over 3820 unique diffusion weightings (the MASSIVE dataset), and double (oscillating) diffusion encoding data (DDE/DODE) of a mouse brain including over 2520 unique data points. A subset of the data sampled in 5 different voxels was openly distributed, and the challenge participants were asked to predict the remaining part of the data. After one year, eight participant teams submitted a total of 80 signal fits. For each submission, we evaluated the mean squared error, the variance of the prediction error and the Bayesian information criteria. The received submissions predicted either multi-shell SDE data (37%) or DODE data (22%), followed by cartesian SDE data (19%) and DDE (18%). Most submissions predicted the signals measured with SDE remarkably well, with the exception of low and very strong diffusion weightings. The prediction of DDE and DODE data seemed more challenging, likely because none of the submissions explicitly accounted for diffusion time and frequency. Next to the choice of the model, decisions on fit procedure and hyperparameters play a major role in the prediction performance, highlighting the importance of optimizing and reporting such choices. This work is a community effort to highlight strength and limitations of the field at representing dMRI acquired with trending encoding schemes, gaining insights into how different models generalize to different tissue types and fiber configurations over a large range of diffusion encodings.

Signal to noise and b-value analysis for optimal intra-voxel incoherent motion imaging in the brain

Intravoxel incoherent motion (IVIM) is a method that can provide quantitative information about perfusion in the human body, in vivo, and without contrast agent. Unfortunately, the IVIM perfusion parameter maps are known to be relatively noisy in the brain, in particular for the pseudo-diffusion coefficient, which might hinder its potential broader use in clinical applications. Therefore, we studied the conditions to produce optimal IVIM perfusion images in the brain. IVIM imaging was performed on a 3-Tesla clinical system in four healthy volunteers, with 16 b values 0, 10, 20, 40, 80, 110, 140, 170, 200, 300, 400, 500, 600, 700, 800, 900 s/mm², repeated 20 times. We analyzed the noise characteristics of the trace images as a function of b-value, and the homogeneity of the IVIM parameter maps across number of averages and sub-sets of the acquired b values. We found two peaks of noise of the trace images as function of b value, one due to thermal noise at high b-value, and one due to physiological noise at low b-value. The selection of b value distribution was found to have higher impact on the homogeneity of the IVIM parameter maps than the number of averages. Based on evaluations, we suggest an optimal b value acquisition scheme for a 12 min scan as 0 (7), 20 (4), 140 (19), 300 (9), 500 (19), 700 (1), 800 (4), 900 (1) s/mm².

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/mm² 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  mm² /sec, D* was 29.7 ± 38.1 × 10^-3  mm² /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.

Characterisation of microvessel blood velocity and segment length in the brain using multi-diffusion-time diffusion-weighted MRI

Multi-diffusion-time diffusion-weighted MRI can probe tissue microstructure, but the method has not been widely applied to the microvasculature. At long diffusion-times, blood flow in capillaries is in the diffusive regime, and signal attenuation is dependent on blood velocity (v) and capillary segment length (l). It is described by the pseudo-diffusion coefficient (D*=vl/6) of intravoxel incoherent motion (IVIM). At shorter diffusion-times, blood flow is in the ballistic regime, and signal attenuation depends on v, and not l. In theory, l could be estimated using D* and v. In this study, we compare the accuracy and repeatability of three approaches to estimating v, and therefore l: the IVIM ballistic model, the velocity autocorrelation model, and the ballistic approximation to the velocity autocorrelation model. Twenty-nine rat datasets from two strains were acquired at 7 T, with b-values between 0 and 1000 smm^-2 and diffusion times between 11.6 and 50 ms. Five rats were scanned twice to assess scan-rescan repeatability. Measurements of l were validated using corrosion casting and micro-CT imaging. The ballistic approximation of the velocity autocorrelation model had lowest bias relative to corrosion cast estimates of l, and had highest repeatability.

Characterization and correction of cardiovascular motion artifacts in diffusion-weighted imaging of the pancreas

PURPOSE

To assess the effects of cardiovascular-induced motion on conventional DWI of the pancreas and to evaluate motion-robust DWI methods in a motion phantom and healthy volunteers.

METHODS

3T DWI was acquired using standard monopolar and motion-compensated gradient waveforms, including in an anatomically accurate pancreas phantom with controllable compressive motion and healthy volunteers (n = 8, 10). In volunteers, highly controlled single-slice DWI using breath-holding and cardiac gating and whole-pancreas respiratory-triggered DWI were acquired. For each acquisition, the ADC variability across volunteers, as well as ADC differences across parts of the pancreas were evaluated.

RESULTS

In motion phantom scans, conventional DWI led to biased ADC, whereas motion-compensated waveforms produced consistent ADC. In the breath-held, cardiac-triggered study, conventional DWI led to heterogeneous DW signals and highly variable ADC across the pancreas, whereas motion-compensated DWI avoided these artifacts. In the respiratory-triggered study, conventional DWI produced heterogeneous ADC across the pancreas (head: 1756 ± 173 × 10^-6 mm² /s; body: 1530 ± 338 × 10^-6 mm² /s; tail: 1388 ± 267 × 10^-6 mm² /s), with ADCs in the head significantly higher than in the tail (P < .05). Motion-compensated ADC had lower variability across volunteers (head: 1277 ± 102 × 10^-6 mm² /s; body: 1204 ± 169 × 10^-6 mm² /s; tail: 1235 ± 178 × 10^-6 mm² /s), with no significant difference (P ≥ .19) across the pancreas.

CONCLUSION

Cardiovascular motion introduces artifacts and ADC bias in pancreas DWI, which are addressed by motion-robust DWI.

Effect of compression garments on muscle perfusion in delayed-onset muscle soreness: A quantitative analysis using intravoxel incoherent motion MR perfusion imaging

The aim of this prospective cohort study was to evaluate the effect of compression garments under resting conditions and after the induction of delayed-onset muscle soreness (DOMS) by MR perfusion imaging using intravoxel incoherent motion (IVIM). Magnetic resonance imaging of both lower legs of 16 volunteers was performed before and after standardized eccentric exercises that induced DOMS. A compression garment (21-22 mmHg) was worn during and for 6 h after exercise on one randomly selected leg. IVIM MR imaging, represented as total muscle perfusion D*f, perfusion fraction f and tissue diffusivity D, were compared between baseline and directly, 30 min, 6 h and 48 h after exhausting exercise with and without compression. Creatine kinase levels and T2-weighted images were acquired at baseline and after 48 h. DOMS was induced in the medial head of the gastrocnemius muscle (MGM) in all volunteers. Compression garments did not show any significant effect on IVIM perfusion parameters at any time point in the MGM or the tibialis anterior muscle (p > 0.05). Microvascular perfusion in the MGM increased significantly in both the compressed and noncompressed leg between baseline measurements and those taken directly after and 30 min after the exercise: the relative median f increased by 31.5% and 24.7% in the compressed and noncompressed leg, respectively, directly after the exercise compared with the baseline value. No significant change in tissue perfusion occurred 48 h after the induction of DOMS compared with baseline. It was concluded that compression garments (21-22 mmHg) do not alter microvascular muscle perfusion at rest, nor do they have any significant effect during the regeneration phase of DOMS. In DOMS, only a short-term effect of increased muscle perfusion (30 min after exercise) was observed, with normalization occurring during regeneration after 6-48 h. The normalization of perfusion independently of compression after 6 h may have implications for diagnostic and therapeutic strategies and for the better understanding of pathophysiological pathways in DOMS.

Association of Head Impact Exposure with White Matter Macrostructure and Microstructure Metrics

Prior studies have reported white matter abnormalities associated with a history of cumulative concussion and/or repetitive head impacts (RHI) in contact sport athletes. Growing evidence suggests these abnormalities may begin as more subtle changes earlier in life in active younger athletes. We investigated the relationship between prior concussion and contact sport exposure with multi-modal white matter microstructure and macrostructure using magnetic resonance imaging. High school and collegiate athletes (n = 121) completed up to four evaluations involving neuroimaging. Linear mixed-effects models examined associations of years of contact sport exposure (i.e., RHI proxy) and prior concussion across multiple metrics of white matter, including total white matter volume, diffusion tensor imaging (DTI) metrics, diffusion kurtosis imaging (DKI) metrics, and quantitative susceptibility mapping (QSM). A significant inverse association between cumulative years of contact sport exposure and QSM was observed, F(1, 237.77) = 4.67, p = 0.032. Cumulative contact sport exposure was also associated with decreased radial diffusivity, F(1, 114.56) = 5.81, p = 0.018, as well as elevated fractional anisotropy, F(1, 115.32) = 5.40, p = 0.022, and radial kurtosis, F(1, 113.45) = 4.03, p = 0.047. In contrast, macroscopic white matter volume was not significantly associated with cumulative contact sport exposure (p > 0.05). Concussion history was not significantly associated with QSM, DTI, DKI, or white matter volume (all, p > 0.05). Cumulative contact sport exposure is associated with subtle differences in white matter microstructure, but not gross white matter macrostructure, in young active athletes. Longitudinal follow-up is required to assess the progression of these findings to determine their contribution to potential adverse effects later in life.

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.

Double diffusion encoding and applications for biomedical imaging

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

Motion-compensated gradient waveforms for tensor-valued diffusion encoding by constrained numerical optimization

PURPOSE

Diffusion-weighted MRI is sensitive to incoherent tissue motion, which may confound the measured signal and subsequent analysis. We propose a "motion-compensated" gradient waveform design for tensor-valued diffusion encoding that negates the effects bulk motion and incoherent motion in the ballistic regime.

METHODS

Motion compensation was achieved by constraining the magnitude of gradient waveform moment vectors. The constraint was incorporated into a numerical optimization framework, along with existing constraints that account for b-tensor shape, hardware restrictions, and concomitant field gradients. We evaluated the efficacy of encoding and motion compensation in simulations, and we demonstrated the approach by linear and planar b-tensor encoding in a healthy heart in vivo.

RESULTS

The optimization framework produced asymmetric motion-compensated waveforms that yielded b-tensors of arbitrary shape with improved efficiency compared with previous designs for tensor-valued encoding, and equivalent efficiency to previous designs for linear (conventional) encoding. Technical feasibility was demonstrated in the heart in vivo, showing vastly improved data quality when using motion compensation. The optimization framework is available online in open source.

CONCLUSION

Our gradient waveform design is both more flexible and efficient than previous methods, facilitating tensor-valued diffusion encoding in tissues in which motion would otherwise confound the signal. The proposed design exploits asymmetric encoding times, a single refocusing pulse or multiple refocusing pulses, and integrates compensation for concomitant gradient effects throughout the imaging volume.

Probing the ballistic microcirculation in placenta using flow-compensated and non-compensated intravoxel incoherent motion imaging

PURPOSE

Intravoxel incoherent motion (IVIM) imaging is widely used to evaluate microcirculatory flow, which consists of diffusive and ballistic flow components. We proposed a joint use of flow-compensated (FC) and non-compensated (NC) diffusion gradients to probe the fraction and velocity of ballistic flow in the placenta.

METHODS

Forty pregnant women were included in this study and scanned on a 1.5T clinical scanner. FC and NC diffusion MRI (dMRI) sequences were achieved using a pair of identical or mirrored bipolar gradients. A joint FC-NC model was established to estimate the fraction (f_b ) and velocity (v_b ) of the ballistic flow. Conventional IVIM parameters (f, D, and D*) were obtained from the FC and NC data, separately. The v_b and f·D*, as placental flow velocity measurements, were correlated with the umbilical-artery Doppler ultrasound indices and gestational ages.

RESULTS

The ballistic flow component can be observed from the difference between the FC and NC dMRI signal decay curves. v_b fitted from the FC-NC model showed strong correlations with umbilical-artery impedance indices, the systolic-to-diastolic (SD) ratio and pulsatility index (PI), with correlation coefficients of 0.65 and 0.62. The f·D* estimated from the NC data positively correlated with SD and PI, while the FC-based f·D* values showed weak negative correlations. Significant gestational-age dependence was also found in the flow velocity measurements.

CONCLUSION

Our results demonstrated the feasibility of using FC and NC dMRI to noninvasively measure ballistic flow velocity in the placenta, which may be used as a new marker to evaluate placenta microcirculation.

Differentiation of high-grade and low-grade intra-axial brain tumors by time-dependent diffusion MRI

INTRODUCTION

Oscillating gradient spin-echo (OGSE) sequences enable acquisitions with shorter diffusion times. There is growing interest in the effect of diffusion time on apparent diffusion coefficient (ADC) values in patients with cancer. However, little evidence exists regarding its usefulness for differentiating between high-grade and low-grade brain tumors. The purpose of this study is to investigate the utility of changes in the ADC value between short and long diffusion times in distinguishing low-grade and high-grade brain tumors.

MATERIAL AND METHODS

Eleven patients with high-grade brain tumors and ten patients with low-grade brain tumors were scanned using a 3 T magnetic resonance imaging with diffusion-weighted imaging (DWI) using OGSE and PGSE (effective diffusion time [Δ_eff]: 6.5 ms and 35.2 ms) and b-values of 0 and 1000 s/mm². Using a region of interest (ROI) analysis of the brain tumors, we measured the ADC for two Δeff (ADC_Δeff) values and computed the subtraction ADC (ΔADC = ADC_{6.5 ms} - ADC_{35.2 ms}) and the relative ADC (ΔADC = (ADC_{6.5 ms} - ADC_{35.2 ms}) / ADC_{35.2 ms} × 100). The maximum values for the subtraction ADC (ΔADC_max) and the relative ADC (rADC_max) on the ROI were compared between low-grade and high-grade tumors using the Wilcoxon rank-sum test. A P-value <.05 was considered significant. The ROIs were also placed in the normal white matter of patients with high- and low-grade brain tumors, and ΔADC_max values were determined.

RESULTS

High-grade tumors had significantly higher ΔADC_max and rADC_max than low-grade tumors. The ΔADC_max values of the normal white matter were lower than the ΔADC_max of high- and low-grade brain tumors.

CONCLUSION

The dependence of ADC values on diffusion time between 6.5 ms and 35.2 ms was stronger in high-grade tumors than in low-grade tumors, suggesting differences in internal tissue structure. This finding highlights the importance of reporting diffusion times in ADC evaluations and might contribute to the grading of brain tumors using DWI.

Limits to flow detection in phase contrast MRI

Pulsed gradient spin echo (PGSE) complex signal behavior becomes dominated by attenuation rather than oscillation when displacements due to flow are similar or less than diffusive displacements. In this "slow-flow" regime, the optimal displacement encoding parameter q for phase contrast velocimetry depends on the diffusive length scale q s l o w = 1 / l D = 1 / 2 D Δ rather than the velocity encoding parameter v _enc = π/(qΔ). The minimum detectable mean velocity using the difference between the phase at +q _slow and -q _slow is 〈 v m i n 〉 = 1 / SNR D / Δ . These theories are then validated and applied to MRI by performing PGSE echo planar imaging experiments on water flowing through a column with a bulk region and a beadpack region at controlled flow rates. Velocities as slow as 6 μm/s are detected with velocimetry. Theories, MRI experimental protocols, and validation on a controlled phantom help to bridge the gap between porous media NMR and pre-clinical phase contrast and diffusion MRI.

Double Diffusion Encoding for Probing Radiation-Induced Microstructural Changes in a Tumor Model: A Proof-of-Concept Study With Comparison to the Apparent Diffusion Coefficient and Histology

BACKGROUND

Microstructure analyses are gaining interest in cancer MRI as an alternative to the conventional apparent diffusion coefficient (ADC), of which the determinants remain unclear.

PURPOSE

To assess the sensitivity of parameters calculated from a double diffusion encoding (DDE) sequence to changes in a tumor's microstructure early after radiotherapy and to compare them with ADC and histology.

STUDY TYPE

Cohort study on experimental tumors.

ANIMAL MODEL

Sixteen WAG/Rij rats grafted with one rhabdomyosarcoma fragment in each thigh. Thirty-one were imaged at days 1 and 4, of which 17 tumors received a 20 Gy radiation dose after the first imagery.

FIELD STRENGTH/SEQUENCE

3T. Diffusion-weighted imaging, DDE with flow compensated, and noncompensated measurements.

ASSESSMENTS

1) To compare, after irradiation, DDE-derived parameters (intracellular fraction, cell size, and cell density) to their histological counterparts (fraction of stained area, minimal Feret diameter, and nuclei count, respectively). 2) To compare percentage changes in DDE-derived parameters and ADC. 3) To evaluate the evolution of DDE-derived parameters describing perfusion.

STATISTICAL TESTS

Wilcoxon rank sum test.

RESULTS

1) Intracellular fraction, cell size, and cell density were respectively lower (-24%, P < 0.001), higher (+7.5%, P < 0.001) and lower (-38%, P < 0.001) in treated tumors as compared to controls. Fraction of stained area, minimal Feret diameter, and nuclei count were respectively lower (-20%, P < 0.001), higher (+28%, P < 0.001), and lower (-34%, P < 0.001) in treated tumors. 2) The magnitude of ADC's percentage change due to irradiation (16.4%) was superior to the one of cell size (8.4%, P < 0.01) but inferior to those of intracellular fraction (35.5%, P < 0.001) and cell density (42%, P < 0.001). 3) After treatment, the magnitude of the vascular fraction's decrease was higher than the increase of flow velocity (33.3%, vs. 13.3%, P < 0.001).

DATA CONCLUSION

The DDE sequence allows quantitatively monitoring the effects of radiotherapy on a tumor's microstructure, whereas ADC only reveals global changes.

EVIDENCE LEVEL

2.

TECHNICAL EFFICACY

Stage 4. J. Magn. Reson. Imaging 2020;52:941-951.

Gibbs point field model quantifies disorder in microvasculature of U87-glioblastoma

Microvascular proliferation in glioblastoma multiforme is a biological key mechanism to facilitate tumor growth and infiltration and a main target for treatment interventions. The vascular architecture can be obtained by Single Plane Illumination Microscopy (SPIM) to evaluate vascular heterogeneity in tumorous tissue. We make use of the Gibbs point field model to quantify the order of regularity in capillary distributions found in the U87 glioblastoma model in a murine model and to compare tumorous and healthy brain tissue. A single model parameter Γ was assigned that is linked to tissue-specific vascular topology through Monte-Carlo simulations. Distributions of the model parameter Γ differ significantly between glioblastoma tissue with mean 〈Γ_G〉=2.1±0.4, as compared to healthy brain tissue with mean 〈Γ_H〉=4.9±0.4, suggesting that the average Γ-value allows for tissue differentiation. These results may be used for diagnostic magnetic resonance imaging, where it has been shown recently that Γ is linked to tissue-inherent relaxation parameters.

Motion-compensated b-tensor encoding for in vivo cardiac diffusion-weighted imaging

Motion is a major confound in diffusion-weighted imaging (DWI) in the body, and it is a common cause of image artefacts. The effects are particularly severe in cardiac applications, due to the nonrigid cyclical deformation of the myocardium. Spin echo-based DWI commonly employs gradient moment-nulling techniques to desensitise the acquisition to velocity and acceleration, ie, nulling gradient moments up to the 2nd order (M2-nulled). However, current M2-nulled DWI scans are limited to encode diffusion along a single direction at a time. We propose a method for designing b-tensors of arbitrary shapes, including planar, spherical, prolate and oblate tensors, while nulling gradient moments up to the 2nd order and beyond. The design strategy comprises initialising the diffusion encoding gradients in two encoding blocks about the refocusing pulse, followed by appropriate scaling and rotation, which further enables nulling undesired effects of concomitant gradients. Proof-of-concept assessment of in vivo mean diffusivity (MD) was performed using linear and spherical tensor encoding (LTE and STE, respectively) in the hearts of five healthy volunteers. The results of the M2-nulled STE showed that (a) the sequence was robust to cardiac motion, and (b) MD was higher than that acquired using standard M2-nulled LTE, where diffusion-weighting was applied in three orthogonal directions, which may be attributed to the presence of restricted diffusion and microscopic diffusion anisotropy. Provided adequate signal-to-noise ratio, STE could significantly shorten estimation of MD compared with the conventional LTE approach. Importantly, our theoretical analysis and the proposed gradient waveform design may be useful in microstructure imaging beyond diffusion tensor imaging where the effects of motion must be suppressed.

A mixed waveform protocol for reduction of the cardiac motion artifact in black-blood diffusion-weighted imaging of the liver

OBJECTIVE

Diffusion-weighted imaging (DWI) in the liver suffers from signal loss due to the cardiac motion artifact, especially in the left liver lobe. The purpose of this work was to improve the image quality of liver DWI in terms of cardiac motion artifact reduction and achievement of black-blood images in low b-value images.

MATERIAL AND METHODS

Ten healthy volunteers (age 20-31 years) underwent MRI examinations at 1.5 T with a prototype DWI sequence provided by the vendor. Two diffusion encodings (i.e. waveforms), monopolar and flow-compensated, and the b-values 0, 20, 50, 100, 150, 600 and 800 s/mm² were used. Two Likert scales describing the severity of the pulsation artifact and the quality of the black-blood state were defined and evaluated by two experienced radiologists. Regions of interest (ROIs) were manually drawn in the right and left liver lobe in each slice and combined to a volume of interest (VOI). The mean and coefficient of variation were calculated for each normalized VOI-averaged signal to assess the severity of the cardiac motion artifact. The ADC was calculated using two b-values once for the monopolar data and once with mixed data, using the monopolar data for the small and the flow-compensated data for the high b-value. A Wilcoxon rank sum test was used to compare the Likert scores obtained for monopolar and flow-compensated data.

RESULTS

At b-values from 20 to 150 s/mm², unlike the flow-compensated diffusion encoding, the monopolar encoding yielded black blood in all images with a negligible signal loss due to the cardiac motion artifact. At the b-values 600 and 800 s/mm², the flow-compensated encoding resulted in a significantly reduced cardiac motion artifact, especially in the left liver lobe, and in a black-blood state. The ADC calculated with monopolar data was significantly higher in the left than in the right liver lobe.

CONCLUSION

It is recommendable to use the following mixed waveform protocol: Monopolar diffusion encodings at small b-values and flow-compensated diffusion encodings at high b-values.

Impact of b-Value Sampling Scheme on Brain IVIM Parameter Estimation in Healthy Subjects

Purpose:

Intravoxel incoherent motion (IVIM) analysis has attracted the interest of the clinical community due to its close relationship with microperfusion. Nevertheless, there is no clear reference protocol for its implementation; one of the questions being which b-value distribution to use. This study aimed to stress the importance of the sampling scheme and to show that an optimized b-value distribution decreases the variance associated with IVIM parameters in the brain with respect to a regular distribution in healthy volunteers.

Methods:

Ten volunteers were included in this study; images were acquired on a 1.5T MR scanner. Two distributions of 16 b-values were used: one considered ‘regular’ due to its close association with that used in other studies, and the other considered ‘optimized’ according to previous studies. IVIM parameters were adjusted according to the bi-exponential model, using two-step method. Analysis was undertaken in ROI defined using in the Automated Anatomical Labeling atlas, and parameters distributions were compared in a total of 832 ROI.

Results:

Maps with fewer speckles were obtained with the ‘optimized’ distribution. Coefficients of variation did not change significantly for the estimation of the diffusion coefficient D but decreased by approximately 39% for the pseudo-diffusion coefficient estimation and by 21% for the perfusion fraction. Distributions of adjusted parameters were found significantly different in 50% of the cases for the perfusion fraction, in 80% of the cases for the pseudo-diffusion coefficient and 17% of the cases for the diffusion coefficient. Observations across brain areas show that the range of average values for IVIM parameters is smaller in the ‘optimized’ case.

Conclusion:

Using an optimized distribution, data are sampled in a way that the IVIM signal decay is better described and less variance is obtained in the fitted parameters. The increased precision gained could help to detect small variations in IVIM parameters.

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.

Effect of intravoxel incoherent motion on diffusion parameters in normal brain

At very low diffusion weighting the diffusion MRI signal is affected by intravoxel incoherent motion (IVIM) caused by dephasing of magnetization due to incoherent blood flow in capillaries or other sources of microcirculation. While IVIM measurements at low diffusion weightings have been frequently used to investigate perfusion in the body as well as in malignant tissue, the effect and origin of IVIM in normal brain tissue is not completely established. We investigated the IVIM effect on the brain diffusion MRI signal in a cohort of 137 radiologically-normal patients (62 male; mean age = 50.2 ± 17.8, range = 18 to 94). We compared the diffusion tensor parameters estimated from a mono-exponential fit at b = 0 and 1000 s/mm² versus at b = 250 and 1000 s/mm². The asymptotic fitting method allowed for quantitative assessment of the IVIM signal fraction f* in specific brain tissue and regions. Our results show a mean (median) percent difference in the mean diffusivity of about 4.5 (4.9)% in white matter (WM), about 7.8 (8.7)% in cortical gray matter (GM), and 4.3 (4.2)% in thalamus. Corresponding perfusion fraction f* was estimated to be 0.033 (0.032) in WM, 0.066 (0.065) in cortical GM, and 0.033 (0.030) in the thalamus. The effect of f* with respect to age was found to be significant in cortical GM (Pearson correlation ρ ​= ​0.35, p ​= ​3*10^-5) and the thalamus (Pearson correlation ρ = 0.20, p = 0.022) with an average increase in f* of 5.17*10^-4/year and 3.61*10^-4/year, respectively. Significant correlations between f* and age were not observed for WM, and corollary analysis revealed no effect of gender on f*. Possible origins of the IVIM effect in normal brain tissue are discussed.

Spectral Diffusion Analysis of Intravoxel Incoherent Motion MRI in Cerebral Small Vessel Disease

Background

Cerebral intravoxel incoherent motion (IVIM) imaging assumes two components. However, more compartments are likely present in pathologic tissue. We hypothesized that spectral analysis using a nonnegative least‐squares (NNLS) approach can detect an additional, intermediate diffusion component, distinct from the parenchymal and microvascular components, in lesion‐prone regions.

Purpose

To investigate the presence of this intermediate diffusion component and its relation with cerebral small vessel disease (cSVD)‐related lesions.

Study Type

Prospective cross‐sectional study.

Population

Patients with cSVD (n = 69, median age 69.8) and controls (n = 39, median age 68.9).

Field Strength/Sequence

Whole‐brain inversion recovery IVIM acquisition at 3.0T.

Assessment

Enlarged perivascular spaces (PVS) were rated by three raters. White matter hyperintensities (WMH) were identified on a fluid attenuated inversion recovery (FLAIR) image using a semiautomated algorithm.

Statistical Tests

Relations between IVIM measures and cSVD‐related lesions were studied using the Spearman's rank order correlation.

Results

NNLS yielded diffusion spectra from which the intermediate volume fraction f_int was apparent between parenchymal diffusion and microvasular pseudodiffusion. WMH volume and the extent of MRI‐visible enlarged PVS in the basal ganglia (BG) and centrum semiovale (CSO) were correlated with f_int in the WMHs, BG, and CSO, respectively. f_int was 4.2 ± 1.7%, 7.0 ± 4.1% and 13.6 ± 7.7% in BG and 3.9 ± 1.3%, 4.4 ± 1.4% and 4.5 ± 1.2% in CSO for the groups with low, moderate, and high number of enlarged PVS, respectively, and increased with the extent of enlarged PVS (BG: r = 0.49, P < 0.01; CSO: r = 0.23, P = 0.02). f_int in the WMHs was 27.1 ± 13.1%, and increased with the WMH volume (r = 0.57, P < 0.01).

Data Conclusion

We revealed the presence of an intermediate diffusion component in lesion‐prone regions of cSVD and demonstrated its relation with enlarged PVS and WMHs. In tissue with these lesions, tissue degeneration or perivascular edema can lead to more freely diffusing interstitial fluid contributing to f_int.

Level of Evidence: 2

Technical Efficacy: Stage 2

J. Magn. Reson. Imaging 2020;51:1170–1180.

Multiple dimensions for random walks

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

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

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

Evidence of the diffusion time dependence of intravoxel incoherent motion in the brain

PURPOSE

To investigate the diffusion time (T_D ) dependence of intravoxel incoherent motion (IVIM) signals in the brain.

METHODS

A 3-compartment IVIM model was proposed to characterize 2 types of microcirculatory flows in addition to tissue water in the brain: flows that cross multiple vascular segments (pseudo-diffusive) and flows that stay in 1 segment (ballistic) within T_D . The model was first evaluated using simulated flow signals. Experimentally, flow-compensated (FC) pulsed-gradient spin-echo (PGSE) and oscillating-gradient spin-echo (OGSE) sequences were tested using a flow phantom and then used to examine IVIM signals in the mouse brain with T_D ranging from ~2.5 ms to 40 ms on an 11.7T scanner.

RESULTS

By fitting the model to simulated flow signals, we demonstrated the T_D dependency of the estimated fraction of pseudo-diffusive flow and the pseudo-diffusion coefficient (D*), which were dictated by the characteristic timescale of microcirculatory flow (τ). Flow phantom experiments validated that the OGSE and FC-PGSE sequences were not susceptible to the change in flow velocity. In vivo mouse brain data showed that both the estimated fraction of pseudo-diffusive flow and D* increased significantly as T_D increased.

CONCLUSION

We demonstrated that IVIM signals measured in the brain are T_D -dependent, potentially because more microcirculatory flows approach the pseudo-diffusive limit as T_D increases with respect to τ. Measuring the T_D dependency of IVIM signals may provide additional information on microvascular flows in the brain.

Multidimensional diffusion MRI with spectrally modulated gradients reveals unprecedented microstructural detail

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

On the Field Strength Dependence of Bi- and Triexponential Intravoxel Incoherent Motion (IVIM) Parameters in the Liver

BACKGROUND

Studies on intravoxel incoherent motion (IVIM) imaging are carried out with different acquisition protocols.

PURPOSE

To investigate the dependence of IVIM parameters on the B₀ field strength when using a bi- or triexponential model.

STUDY TYPE

Prospective.

STUDY POPULATION

20 healthy volunteers (age: 19-28 years).

FIELD STRENGTH/SEQUENCE

Volunteers were examined at two field strengths (1.5 and 3T). Diffusion-weighted images of the abdomen were acquired at 24 b-values ranging from 0.2 to 500 s/mm² .

ASSESSMENT

ROIs were manually drawn in the liver. Data were fitted with a bi- and a triexponential IVIM model. The resulting parameters were compared between both field strengths.

STATISTICAL TESTS

One-way analysis of variance (ANOVA) and Kruskal-Wallis test were used to test the obtained IVIM parameters for a significant field strength dependency.

RESULTS

At b-values below 6 s/mm² , the triexponential model provided better agreement with the data than the biexponential model. The average tissue diffusivity was D = 1.22/1.00 μm² /msec at 1.5/3T. The average pseudodiffusion coefficients for the biexponential model were D^* = 308/260 μm² /msec at 1.5/3T; and for the triexponential model D 1 * = 81.3/65.9 μm² /msec, D 2 * = 2453/2333 μm² /msec at 1.5/3T. The average perfusion fractions for the biexponential model were f = 0.286/0.303 at 1.5/3T; and for the triexponential model f₁ = 0.161/0.174 and f₂ = 0.152/0.159 at 1.5/3T. A significant B₀ dependence was only found for the biexponential pseudodiffusion coefficient (ANOVA/KW P = 0.037/0.0453) and tissue diffusivity (ANOVA/KW: P < 0.001).

DATA CONCLUSION

Our experimental results suggest that triexponential pseudodiffusion coefficients and perfusion fractions obtained at different field strengths could be compared across different studies using different B₀ . However, it is recommended to take the field strength into account when comparing tissue diffusivities or using the biexponential IVIM model. Considering published values for oxygenation-dependent transversal relaxation times of blood, it is unlikely that the two blood compartments of the triexponential model represent venous and arterial blood.

LEVEL OF EVIDENCE

1 Technical Efficacy Stage: 2 J. Magn. Reson. Imaging 2019;50:1883-1892.

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.

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.

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 robust deconvolution method to disentangle multiple water pools in diffusion MRI

The diffusion-weighted magnetic resonance imaging (dMRI) signal measured in vivo arises from multiple diffusion domains, including hindered and restricted water pools, free water and blood pseudo-diffusion. Not accounting for the correct number of components can bias metrics obtained from model fitting because of partial volume effects that are present in, for instance, diffusion tensor imaging (DTI) and diffusion kurtosis imaging (DKI). Approaches that aim to overcome this shortcoming generally make assumptions about the number of considered components, which are not likely to hold for all voxels. The spectral analysis of the dMRI signal has been proposed to relax assumptions on the number of components. However, it currently requires a clinically challenging signal-to-noise ratio (SNR) and accounts only for two diffusion processes defined by hard thresholds. In this work, we developed a method to automatically identify the number of components in the spectral analysis, and enforced its robustness to noise, including outlier rejection and a data-driven regularization term. Furthermore, we showed how this method can be used to take into account partial volume effects in DTI and DKI fitting. The proof of concept and performance of the method were evaluated through numerical simulations and in vivo MRI data acquired at 3 T. With simulations our method reliably decomposed three diffusion components from SNR = 30. Biases in metrics derived from DTI and DKI were considerably reduced when components beyond hindered diffusion were taken into account. With the in vivo data our method determined three macro-compartments, which were consistent with hindered diffusion, free water and pseudo-diffusion. Taking free water and pseudo-diffusion into account in DKI resulted in lower mean diffusivity and higher fractional anisotropy values in both gray and white matter. In conclusion, the proposed method allows one to determine co-existing diffusion compartments without prior assumptions on their number, and to account for undesired signal contaminations within clinically achievable SNR levels.