Impact of the calculation algorithm on biexponential fitting of diffusion-weighted MRI in upper abdominal organs

ssVERDICT: Self-supervised VERDICT-MRI for enhanced prostate tumor characterization

PURPOSE

Demonstrating and assessing self-supervised machine-learning fitting of the VERDICT (vascular, extracellular and restricted diffusion for cytometry in tumors) model for prostate cancer.

METHODS

We derive a self-supervised neural network for fitting VERDICT (ssVERDICT) that estimates parameter maps without training data. We compare the performance of ssVERDICT to two established baseline methods for fitting diffusion MRI models: conventional nonlinear least squares and supervised deep learning. We do this quantitatively on simulated data by comparing the Pearson's correlation coefficient, mean-squared error, bias, and variance with respect to the simulated ground truth. We also calculate in vivo parameter maps on a cohort of 20 prostate cancer patients and compare the methods' performance in discriminating benign from cancerous tissue via Wilcoxon's signed-rank test.

RESULTS

In simulations, ssVERDICT outperforms the baseline methods (nonlinear least squares and supervised deep learning) in estimating all the parameters from the VERDICT prostate model in terms of Pearson's correlation coefficient, bias, and mean-squared error. In vivo, ssVERDICT shows stronger lesion conspicuity across all parameter maps, and improves discrimination between benign and cancerous tissue over the baseline methods.

CONCLUSION

ssVERDICT significantly outperforms state-of-the-art methods for VERDICT model fitting and shows, for the first time, fitting of a detailed multicompartment biophysical diffusion MRI model with machine learning without the requirement of explicit training labels.

Application of Intravoxel Incoherent Motion in Clinical Liver Imaging: A Literature Review

Intravoxel incoherent motion (IVIM) modeling is a widely used double-exponential model for describing diffusion-weighted imaging (DWI) signal, with a slow component related to pure molecular diffusion and a fast component associated with microcirculatory perfusion, which compensates for the limitations of traditional DWI. IVIM is a noninvasive technique for obtaining liver pathological information and characterizing liver lesions, and has potential applications in the initial diagnosis and treatment monitoring of liver diseases. Recent studies have demonstrated that IVIM-derived parameters are useful for evaluating liver lesions, including nonalcoholic fatty liver disease (NAFLD), liver fibrosis and liver tumors. However, the results are not stable. Therefore, it is necessary to summarize the current applications of IVIM in liver disease research, identify existing shortcomings, and point out the future development direction. In this review, we searched for studies related to hepatic IVIM-DWI applications over the past two decades in the PubMed database. We first introduce the fundamental principles and influential factors of IVIM, and then discuss its application in NAFLD, liver fibrosis, and focal hepatic lesions. It has been found that IVIM is still unstable in ensuring the robustness and reproducibility of measurements in the assessment of liver fibrosis grade and liver tumors differentiation, due to inconsistent and substantial overlap in the range of IVIM-derived parameters for different fibrotic stages. In the end, the future direction of IVIM-DWI in the assessment of liver diseases is discussed, emphasizing the need for further research on the stability of IVIM-derived parameters, particularly perfusion-related parameters, in order to promote the clinical practice of IVIM-DWI. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY: Stage 3.

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.

Calculation of intravoxel incoherent motion parameter maps using a kernelized total difference-based method

Quantitative analysis of diffusion-weighted magnetic resonance imaging (DW-MRI) has been explored for many clinical applications since its development. In particular, the intravoxel incoherent motion (IVIM) model for DW-MRI has been commonly utilized in various organs. However, because of the presence of excessive noise, the IVIM parameter maps obtained from pixel-wise fitting are often unreliable. In this study, we propose a kernelized total difference-based curve-fitting method to estimate the IVIM parameters. Simulated DW-MRI data at five signal-to-noise ratios (i.e., 10, 20, 30, 50, and 100) and real abdominal DW-MRI data acquired on a 1.5-T MRI scanner with nine b-values (i.e., 0, 10, 25, 50, 100, 200, 300, 400, and 500 s/mm2) and six diffusion-encoding gradient directions were used to evaluate the performance of the proposed method. The results were compared with those obtained by three existing methods: trust-region reflective (TRR) algorithm, Bayesian probability (BP), and deep neural network (DNN). Our simulation results showed that the proposed method outperformed the other three comparing methods in terms of root-mean-square error. Moreover, the proposed method could preserve small details in the estimated IVIM parameter maps. The experimental results showed that, compared with the TRR method, the proposed method as well as the BP (and DNN) method could reduce the overestimation of the pseudodiffusion coefficient and improve the quality of IVIM parameter maps. For all studied abdominal organs except the pancreas, both the proposed method and the BP method could provide IVIM parameter estimates close to the reference values; the former had higher precision. The kernelized total difference-based curve-fitting method has the potential to improve the reliability of IVIM parametric imaging.

Intravoxel incoherent motion magnetic resonance imaging to assess early tumor response to radiation therapy: Review and future directions

Early prediction of radiation response by imaging is a dynamic field of research and it can be obtained using a variety of noninvasive magnetic resonance imaging methods. Recently, intravoxel incoherent motion (IVIM) has gained interest in cancer imaging. IVIM carries both diffusion and perfusion information, making it a promising tool to assess tumor response. Here, we briefly introduced the basics of IVIM, reviewed existing studies of IVIM in various type of tumors during radiotherapy in order to show whether IVIM is a useful technique for an early assessment of radiation response. 31/40 studies reported an increase of IVIM parameters during radiotherapy compared to baseline. In 27 studies, this increase was higher in patients with good response to radiotherapy. Future directions including implementation of IVIM on MR-Linac and its limitation are discussed. Obtaining new radiologic biomarkers of radiotherapy response could open the way for a more personalized, biology-guided radiation therapy.

Repeatability and reproducibility comparisons of liver IVIM imaging with free-breathing or respiratory-triggered sequences

For liver intravoxel incoherent motion (IVIM) data acquisition, respiratory-triggering (RT) MRI is commonly used, and there are strong motivations to shorten the scan duration. For the same scan duration, more b values or higher numbers of excitations can be allowed for free-breathing (FB) imaging than for RT. We studied whether FB can be used to replace RT when careful IVIM image acquisition and image processing are conducted. MRI data of 22 healthy participants were acquired using a 3.0 T scanner. Diffusion imaging was based on a single-shot spin-echo-type echo-planar sequence and 16 b values of 0, 2, 4, 7, 10, 15, 20, 30, 46, 60, 72, 100, 150, 200, 400, and 600 s/mm2 . Each subject attended two scan sessions with an interval of 10-20 days. For each scan session, a subject was scanned twice, first with RT and then with FB. The mean image acquisition time was 5.4 min for FB and 10.8 min for RT. IVIM parameters were calculated with bi-exponential model segmented fitting with a threshold b value of 60 s/mm2 , and fitting started from b = 2 s/mm2 . There was no statistically significant difference between IVIM parameters measured with FB imaging or RT imaging. Perfusion fraction ICC (intraclass correlation coefficient) for FB imaging and RT imaging in the same scan session was 0.824. For perfusion fraction, wSD (within-subject standard deviation), BA (Bland-Altman) difference, BA 95% limit, and ICC were 0.022, 0.0001, -0.0635~0.0637, and 0.687 for FB and 0.031, 0.0122, -0.0723~0.0967, and 0.611 for RT. For Dslow (×10-3  s/mm2 ), wSD, BA difference, BA 95% limit, and ICC were 0.057, 0.0268, -0.1258~0.1793, and 0.471 for FB and 0.073, -0.0078, -0.2170-0.2014, and <0.4 for RT. The Dfast coefficient of variation was 0.20 for FB imaging and 0.28 for RT imaging. All reproducibility indicators slightly favored FB imaging.

Radiologists' Rating for Comparative Qualitative Assessment of Intravoxel Incoherent Motion Using Novel Analysis Methods

OBJECTIVE

The objective was to assess qualitative interpretability and quantitative precision and reproducibility of intravoxel incoherent motion ( IVIM) parametric images evaluated using novel IVIM analysis methods for diagnostic accuracy.

METHODS

Intravoxel incoherent motion datasets of 55 patients (male/female = 41:14; age = 17.8 ± 5.5 years) with histopathology-proven osteosarcoma were analyzed. Intravoxel incoherent motion parameters-diffusion coefficient ( D ), perfusion fraction ( f ), and perfusion coefficient ( D* )-were estimated using 5 IVIM analysis methods-(i) biexponential (BE) model, (ii) BE-segmented fitting 2-parameter (BESeg-2), (iii) BE-segmented fitting 1-parameter (BESeg-1), (iv) BE model with total variation penalty function (BE + TV), and (v) BE model with Huber penalty function (BE + HPF). Qualitative scoring in a 5-point Likert scale (uninterpretable: 1; poor: 2; fair: 3; good: 4; excellent: 5) was performed by 2 radiologists for 4 criteria: (a) tumor shape and margin, (b) morphologic correlation, (c) noise suppression, and (d) overall interpretability. Interobserver agreement was evaluated using Spearman rank-order correlation ( rs ). Precision and reproducibility were evaluated using within-subject coefficient of variation (wCV) and between-subject coefficient of variation (bCV).

RESULTS

BE + TV and BE + HPF produced significantly ( P < 10 -3 ) higher qualitative scores for D (fair-good [3.3-3.8]) than BE (poor [2.3]) and for D* (poor-fair [2.2-2.7]) and f (fair-good [3.2-3.8]) than BE, BESeg-2, and BESeg-1 ( D* : uninterpretable-poor [1.3-1.9] and f : poor-fair [1.5-3]). Interobserver agreement for qualitative scoring was rs = 0.48-0.59, P < 0.009. BE + TV and BE + HPF showed significantly ( P < 0.05) improved reproducibility in estimating D (wCV: 24%-31%, bCV: 21%-31% improvement) than the BE method and D* (wCV: 4%-19%, bCV: 5%-19% improvement) and f (wCV: 25%-49%, bCV: 25%-47% improvement) than BE, BESeg-2, and BESeg-1 methods.

CONCLUSIONS

BE + TV and BE + HPF demonstrated qualitatively and quantitatively improved IVIM parameter estimation and may be considered for clinical use further.

Highly accelerated multi-shot intravoxel incoherent motion diffusion-weighted imaging in brain enabled by parametric POCS-based multiplexed sensitivity encoding

Recently, intravoxel incoherent motion (IVIM) diffusion-weighted imaging (DWI) has also been demonstrated as an imaging tool for applications in neurological and neurovascular diseases. However, the use of single-shot diffusion-weighted echo-planar imaging for IVIM DWI acquisition leads to suboptimal data quality: for instance, geometric distortion and deteriorated image quality at high spatial resolution. Although the recently commercialized multi-shot acquisition methods, such as multiplexed sensitivity encoding (MUSE), can attain high-resolution and high-quality DWI with signal-to-noise ratio (SNR) performance superior to that of the conventional parallel imaging method, the prolonged scan time associated with multi-shot acquisition is impractical for routine IVIM DWI. This study proposes an acquisition and reconstruction framework based on parametric-POCSMUSE to accelerate the four-shot IVIM DWI with 70% reduction of total scan time (13 min 8 s versus 4 min 8 s). First, the four-shot IVIM DWI scan with 17 b values was accelerated by acquiring only one segment per b value except for b values of 0 and 600 s/mm2 . Second, an IVIM-estimation scheme was integrated into the parametric-POCSMUSE to enable joint reconstruction of multi-b images from under-sampled four-shot IVIM DWI data. In vivo experiments on both healthy subjects and patients show that the proposed framework successfully produced multi-b DW images with significantly higher SNRs and lower reconstruction errors than did the conventional acceleration method based on parallel imaging. In addition, the IVIM quantitative maps estimated from the data produced by the proposed framework showed quality comparable to that of fully sampled MUSE-reconstructed images, suggesting that the proposed framework can enable highly accelerated multi-shot IVIM DWI without sacrificing data quality. In summary, the proposed framework can make multi-shot IVIM DWI feasible in a routine MRI examination, with reasonable scan time and improved geometric fidelity.

Towards detection of early response in neoadjuvant chemotherapy of breast cancer using Bayesian intravoxel incoherent motion

Introduction

The early identification of good responders to neoadjuvant chemotherapy (NACT) holds a significant potential in the optimal treatment of breast cancer. A recent Bayesian approach has been postulated to improve the accuracy of the intravoxel incoherent motion (IVIM) model for clinical translation. This study examined the prediction and early sensitivity of Bayesian IVIM to NACT response.

Materials and methods

Seventeen female patients with breast cancer were scanned at baseline and 16 patients were scanned after Cycle 1. Tissue diffusion and perfusion from Bayesian IVIM were calculated at baseline with percentage change at Cycle 1 computed with reference to baseline. Cellular proliferative activity marker Ki-67 was obtained semi-quantitatively with percentage change at excision computed with reference to core biopsy.

Results

The perfusion fraction showed a significant difference (p = 0.042) in percentage change between responder groups at Cycle 1, with a decrease in good responders [-7.98% (-19.47-1.73), n = 7] and an increase in poor responders [10.04% (5.09-28.93), n = 9]. There was a significant correlation between percentage change in perfusion fraction and percentage change in Ki-67 (p = 0.042). Tissue diffusion and pseudodiffusion showed no significant difference in percentage change between groups at Cycle 1, nor was there a significant correlation against percentage change in Ki-67. Perfusion fraction, tissue diffusion, and pseudodiffusion showed no significant difference between groups at baseline, nor was there a significant correlation against Ki-67 from core biopsy.

Conclusion

The alteration in tumour perfusion fraction from the Bayesian IVIM model, in association with cellular proliferation, showed early sensitivity to good responders in NACT.

Clinical trial registration

https://clinicaltrials.gov/ct2/show/NCT03501394, identifier NCT03501394.

Intravoxel Incoherent Motion Quantification Dependent on Measurement SNR and Tissue Perfusion: A Simulation Study

Background

The intravoxel incoherent motion (IVIM) model extracts both functional and structural information of a tissue using motion-sensitizing gradients.

Objective

The Objective of the present work is to investigate the impact of signal to noise ratio (SNR) and physiologic conditions on the validity of IVIM parameters.

Material and Methods

This study is a simulation study, modeling IVIM at a voxel, and also done 10,000 times for every single simulation. Complex noises with various standard deviations were added to signal in-silico to investigate SNR effects on output validity. Besides, some blood perfusion situations for different tissues were considered based on their physiological range to explore the impacts of blood fraction at each voxel on the validity of the IVIM outputs. Coefficient variation (CV) and bias of the estimations were computed to assess the validity of the IVIM parameters.

Results

This study has shown that the validity of IVIM output parameters highly depends on measurement SNR and physiologic characteristics of the studied organ.

Conclusion

IVIM imaging could be useful if imaging parameters are correctly selected for each specific organ, considering hardware limitations.

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.

Probing Renal Microstructure and Function with Advanced Diffusion MRI: Concepts, Applications, Challenges, and Future Directions

Diffusion measurements in the kidney are affected not only by renal microstructure but also by physiological processes (i.e., glomerular filtration, water reabsorption, and urine formation). Because of the superposition of passive tissue diffusion, blood perfusion, and tubular pre-urine flow, the limitations of the monoexponential apparent diffusion coefficient (ADC) model in assessing pathophysiological changes in renal tissue are becoming apparent and motivate the development of more advanced diffusion-weighted imaging (DWI) variants. These approaches take advantage of the fact that the length scale probed in DWI measurements can be adjusted by experimental parameters, including diffusion-weighting, diffusion gradient directions and diffusion time. This forms the basis by which advanced DWI models can be used to capture not only passive diffusion effects, but also microcirculation, compartmentalization, tissue anisotropy. In this review, we provide a comprehensive overview of the recent advancements in the field of renal DWI. Following a short introduction on renal structure and physiology, we present the key methodological approaches for the acquisition and analysis of renal DWI data, including intravoxel incoherent motion (IVIM), diffusion tensor imaging (DTI), non-Gaussian diffusion, and hybrid IVIM-DTI. We then briefly summarize the applications of these methods in chronic kidney disease and renal allograft dysfunction. Finally, we discuss the challenges and potential avenues for further development of renal DWI. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY: Stage 2.

Estimation of intravoxel incoherent motion (IVIM) parameters in vertebral bone marrow: a comparative study of five algorithms

OBJECTIVES

To access the performances of different algorithms for quantification of Intravoxel incoherent motion (IVIM) parameters D, f, [Formula: see text] in Vertebral Bone Marrow (VBM).

MATERIALS AND METHODS

Five algorithms were studied: four deterministic algorithms (the One-Step and three segmented methods: Two-Step, Three-Step, and Fixed-[Formula: see text] algorithm) based on the least-squares (LSQ) method and a Bayesian probabilistic algorithm. Numerical simulations and quantification of IVIM parameters D, f, [Formula: see text] in vivo in vertebral bone marrow, were done on six healthy volunteers. The One-way repeated-measures analysis of variance (ANOVA) followed by Bonferroni's multiple comparison test (p value = 0.05) was applied.

RESULTS

In numerical simulations, the Bayesian algorithm provided the best estimation of D, f, [Formula: see text] compared to the deterministic algorithms. In vivo VBM-IVIM, the values of D and f estimated by the Bayesian algorithm were close to those of the One-Step method, in contrast to the three segmented methods.

DISCUSSION

The comparison of the five algorithms indicates that the Bayesian algorithm provides the best estimation of VBM-IVIM parameters, in both numerical simulations and in vivo data.

Quantitative tumor burden imaging parameters of the spleen at MRI for predicting treatment response in patients with acute leukemia

Objectives

To study the value of standardized volume and intravoxel incoherent motion (IVIM) parameters of the spleen based on tumor burden for predicting treatment response in newly diagnosed acute leukemia (AL).

Methods

Patients with newly diagnosed AL were recruited and underwent abdominal IVIM diffusion-weighted imaging within one week before the first induction chemotherapy. Quantitative parameters of magnetic resonance imaging (MRI) included the standardized volume (representing volumetric tumor burden) and IVIM parameters (standard apparent diffusion coefficient [sADC]; pure diffusion coefficient [D]; pseudo-diffusion coefficient [D∗]; and pseudo-perfusion fraction [f], representing functional tumor burden) of the spleen. Clinical biomarkers of tumor burden were collected. Patients were divided into complete remission (CR) and non-CR groups according to the treatment response after the first standardized induction chemotherapy, and the MRI and clinical parameters were compared between the two groups. The correlations of MRI parameters with clinical biomarkers were analyzed. Multivariate logistic regression was performed to determine the independent predictors for treatment response. Receiver operating characteristic curves were used to analyze the predicted performance.

Results

76 AL patients (CR: n = 43; non-CR: n = 33) were evaluated. Standardized spleen volume, sADC, D, f, white blood cell counts, and lactate dehydrogenase were significantly different between CR and non-CR groups (all p < 0.05). Standardized spleen volume, sADC, and D were correlated with white blood cell and lactate dehydrogenase, and f was correlated with lactate dehydrogenase (all p < 0.05). Standardized spleen volume (hazard ratio = 4.055, p = 0.042), D (hazard ratio = 0.991, p = 0.027), and f (hazard ratio = 1.142, p = 0.008) were independent predictors for treatment response, and the combination of standardized spleen volume, D, and f showed more favorable discrimination (area under the curve = 0.856) than individual predictors.

Conclusion

Standardized volume, D, and f of the spleen could be used to predict treatment response in newly diagnosed AL, and the combination of morphological and functional parameters would further improve the predicted performance. IVIM parameters of the spleen may be viable indicators for evaluating functional tumor burden in AL.

Deep learning intravoxel incoherent motion modeling: Exploring the impact of training features and learning strategies

PURPOSE

The development of advanced estimators for intravoxel incoherent motion (IVIM) modeling is often motivated by a desire to produce smoother parameter maps than least squares (LSQ). Deep neural networks show promise to this end, yet performance may be conditional on a myriad of choices regarding the learning strategy. In this work, we have explored potential impacts of key training features in unsupervised and supervised learning for IVIM model fitting.

METHODS

Two synthetic data sets and one in-vivo data set from glioma patients were used in training of unsupervised and supervised networks for assessing generalizability. Network stability for different learning rates and network sizes was assessed in terms of loss convergence. Accuracy, precision, and bias were assessed by comparing estimations against ground truth after using different training data (synthetic and in vivo).

RESULTS

A high learning rate, small network size, and early stopping resulted in sub-optimal solutions and correlations in fitted IVIM parameters. Extending training beyond early stopping resolved these correlations and reduced parameter error. However, extensive training resulted in increased noise sensitivity, where unsupervised estimates displayed variability similar to LSQ. In contrast, supervised estimates demonstrated improved precision but were strongly biased toward the mean of the training distribution, resulting in relatively smooth, yet possibly deceptive parameter maps. Extensive training also reduced the impact of individual hyperparameters.

CONCLUSION

Voxel-wise deep learning for IVIM fitting demands sufficiently extensive training to minimize parameter correlation and bias for unsupervised learning, or demands a close correspondence between the training and test sets for supervised learning.

Effects of fitting methods, high b-values and image quality on diffusion and perfusion quantification and reproducibility in the calf

PURPOSES

The study aimed to optimize diffusion-weighted imaging (DWI) image acquisition and analysis protocols in calf muscles by investigating the effects of different model-fitting methods, image quality, and use of high b-value and constraints on parameters of interest (POIs). The optimized modeling methods were used to select the optimal combinations of b-values, which will allow shorter acquisition time while achieving the same reliability as that obtained using 16 b-values.

METHODS

Test-retest baseline and high-quality DWI images of ten healthy volunteers were acquired on a 3T MR scanner, using 16 b-values, including a high b-value of 1200 s/mm², and structural T1-weighted images for calf muscle delineation. Three and six different fitting methods were used to derive ADC from monoexponential (ME) model and D_d, f_p, and D_p from intravoxel incoherent motion (IVIM) model, with or without the high b-value. The optimized ME and IVIM models were then used to determine the optimal combinations of b-values, obtainable with the least number of b-values, using the selection criteria of coefficient of variance (CV) ≤10% for all POIs.

RESULTS

The find minimum multivariate algorithm was more flexible and yielded smaller fitting errors. The 2-steps fitting method, with fixed D_d, performed the best for IVIM model. The inclusion of high b-value reduced outliers, while constraints improved 2-steps fitting only.

CONCLUSIONS

The optimal numbers of b-values for ME and IVIM models were nine and six b-values respectively. Test-retest reliability analyses showed that only ADC and D_d were reliable for calf diffusion evaluation, with CVs of 7.22% and 4.09%.

Multiparametric renal magnetic resonance imaging: A reproducibility study in renal allografts with stable function

Monitoring renal allograft function after transplantation is key for the early detection of allograft impairment, which in turn can contribute to preventing the loss of the allograft. Multiparametric renal MRI (mpMRI) is a promising noninvasive technique to assess and characterize renal physiopathology; however, few studies have employed mpMRI in renal allografts with stable function (maintained function over a long time period). The purposes of the current study were to evaluate the reproducibility of mpMRI in transplant patients and to characterize normal values of the measured parameters, and to estimate the labeling efficiency of Pseudo-Continuous Arterial Spin Labeling (PCASL) in the infrarenal aorta using numerical simulations considering experimental measurements of aortic blood flow profiles. The subjects were 20 transplant patients with stable kidney function, maintained over 1 year. The MRI protocol consisted of PCASL, intravoxel incoherent motion, and T1 inversion recovery. Phase contrast was used to measure aortic blood flow. Renal blood flow (RBF), diffusion coefficient (D), pseudo-diffusion coefficient (D*), flowing fraction ( f ), and T1 maps were calculated and mean values were measured in the cortex and medulla. The labeling efficiency of PCASL was estimated from simulation of Bloch equations. Reproducibility was assessed with the within-subject coefficient of variation, intraclass correlation coefficient, and Bland-Altman analysis. Correlations were evaluated using the Pearson correlation coefficient. The significance level was p less than 0.05. Cortical reproducibility was very good for T1, D, and RBF, moderate for f , and low for D*, while medullary reproducibility was good for T1 and D. Significant correlations in the cortex between RBF and f (r = 0.66), RBF and eGFR (r = 0.64), and D* and eGFR (r = -0.57) were found. Normal values of the measured parameters employing the mpMRI protocol in kidney transplant patients with stable function were characterized and the results showed good reproducibility of the techniques.

An unsupervised convolutional neural network method for estimation of intravoxel incoherent motion parameters

Objective. Intravoxel incoherent motion (IVIM) imaging obtained by fitting a biexponential model to multiple b-value diffusion-weighted magnetic resonance imaging (DW-MRI) has been shown to be a promising tool for different clinical applications. Recently, several deep neural network (DNN) methods were proposed to generate IVIM imaging. Approach. In this study, we proposed an unsupervised convolutional neural network (CNN) method for estimation of IVIM parameters. We used both simulated and real abdominal DW-MRI data to evaluate the performance of the proposed CNN-based method, and compared the results with those obtained from a non-linear least-squares fit (TRR, trust-region reflective algorithm) and a feed-forward backward-propagation DNN-based method. Main results. The simulation results showed that both the DNN- and CNN-based methods had lower coefficients of variation than the TRR method, but the CNN-based method provided more accurate parameter estimates. The results obtained from real DW-MRI data showed that the TRR method produced many biased IVIM parameter estimates that hit the upper and lower parameter bounds. In contrast, both the DNN- and CNN-based methods yielded less biased IVIM parameter estimates. Overall, the perfusion fraction and diffusion coefficient obtained from the DNN- and CNN-based methods were close to literature values. However, compared with the CNN-based method, both the TRR and DNN-based methods tended to yield increased pseudodiffusion coefficients (55%–180%). Significance. Our preliminary results suggest that it is feasible to estimate IVIM parameters using CNN.

Reproducibility of apparent diffusion coefficient measurement in normal prostate peripheral zone at 1.5T MRI

PURPOSE

The purpose of this study was to quantify the influence of factors of variability on apparent diffusion coefficient (ADC) estimation in the normal prostate peripheral zone (PZ).

MATERIALS AND METHODS

Fifty healthy volunteers underwent in 2017 (n = 17) or 2020 (n = 33) two-point (0, 800 s/mm²) prostate diffusion-weighted imaging in the morning on 1.5 T scanners A and B from different manufacturers. Additional five-point (50, 150, 300, 500, 800 s/mm²) acquisitions were performed on scanner B in the morning and evening. ADC was measured in PZ at midgland using ADC maps reconstructed with various b-value combinations. ADC distributions from 2017 and 2020 were compared using Wilcoxon rank sum test. ADC obtained in the same volunteers were compared using Bland Altman methodology. The 95% confidence interval upper limit of the repeatability/reproducibility coefficient defined the lowest detectable ADC difference.

RESULTS

Forty-nine participants with a mean age of 24.6 ± 3.8 [SD] years (range: 21-37 years) were finally included. ADC distributions from 2017 and 2020 were not significantly different and were combined. Despite high individual variability, there was no significant bias (10 × 10^-6 mm²/s, P = 0.58) between ADC measurements made on both scanners. On scanner B, differences in lowest b-values chosen within the 0-500 s/mm² range for two-point ADC computation induced significant biases (56-109 × 10^-6 mm²/s, P < 0.0001). ADC was significantly lower in the morning (bias: 33 × 10^-6 mm²/s, P = 0.006). The number of b-values had little influence on ADC values. The lowest detectable ADC difference varied from 85 × 10^-6 to 311 × 10^-6 mm²/s across scanners, b-value combinations and periods of the day.

CONCLUSIONS

The MRI scanner, the lowest b-value used and the period of the day induce substantial variability in ADC computation.

Bi-exponential fitting excluding b=0 data improves the scan-rescan stability of liver IVIM parameter measures and particularly so for the perfusion fraction

BACKGROUND

A prerequisite to translating intravoxel incoherent motion (IVIM) imaging into meaningful clinical applications is sufficient scan-rescan reproducibility. This study aims to confirm the hypothesis that IVIM data fitting by not using b=0 images will improve the stability of liver IVIM measurement.

METHODS

Healthy volunteers' liver IVIM images were prospectively acquired using a 1.5-T magnet or a 3.0 T with 16 b-values. Repeatability study subjects were scanned twice during the same session, resulted in 35 paired scans for 35 subjects (11 men, mean age: 41.82 years, range: 32-60 years; 24 women, mean age: 42.67 years, range: 20-71 years). IVIM analysis was performed with full-fitting and segmented-fitting with a threshold b-value of 60 s/mm2, and fitting started from b=0 s/mm2 or from b=2 s/mm2. Reproducibility study subjects were scanned and then rescanned with an interval of 5-18 days, resulted in 20 paired scans for 11 subjects (4 men, mean age: 26.25 years, range: 25-27 years; 7 women, mean age: 25.57 years, range: 24-27 years). IVIM analysis was performed with segmented-fitting with a threshold b-value of 50 s/mm2, and fitting started from b=0 s/mm2 or from b=3 s/mm2.

RESULTS

Fitting without b=0 data generally improved the repeatability and reproducibility for both PF and Dslow, and particularly so for PF. For with b=0 data segmented fitting repeatability, PF had within-subject standard deviation of 0.019, bland-Atman 75% agreement limit of -31.52% to 28.35%, and ICC of 0.647, while these values were 0.009, -20.78% to 16.86%, and 0.837 for without b=0 analysis. Though the repeatability and reproducibility for Dfast generally also improved, they remained suboptimal. Measurement stability was better for repeatability than for reproducibility.

CONCLUSIONS

Scan-rescan repeatability and reproducibility of liver IVIM parameters can be improved by fitting without b=0 data, which is particularly so for PF.

Perfusion-Diffusion Ratio: A New IVIM Approach in Differentiating Solid Benign and Malignant Primary Lesions of the Liver

Introduction

In order to improve the efficacy of intravoxel incoherent motion (IVIM) parameters in characterising specific tissues, a new concept is introduced: the perfusion–diffusion ratio (PDR), which expresses the relationship between the signal S(b) decline rate as a result of IVIM and the rate of signal S(b) decline due to diffusion. The aim of this study was to investigate this novel approach in the differentiation of solid primary liver lesions. Material and Methods. Eighty-three patients referred for liver MRI between August 2017 and January 2020 with a suspected liver tumour were prospectively examined with the standard liver MRI protocol extended by DWI-IVIM sequence. Patients with no liver lesions, haemangiomas, or metastases were excluded. The final study population consisted of 34 patients with primary solid liver masses, 9 with FNH, 4 with regenerative nodules, 10 with HCC, and 11 with CCC. The PDR coefficient was introduced, defined as the ratio of the rate of signal S(b) decrease due to the IVIM effect to the rate of signal S(b) decrease due to the diffusion process, for b = 0.

Results

No significant differences were found between benign and malignant lesions in the case of IVIM parameters (f, D, or D^∗) and ADC. Significant differences were observed only for PDR, with lower values for malignant lesions (p = 0.03). The ROC analysis yielded an AUC value for PDR equal to 0.74, with a cut-off value of 5.06, sensitivity of 81%, specificity of 77%, and accuracy of 79%.

Conclusion

PDR proved to be more effective than IVIM parameters and ADC in the differentiation of solid benign and malignant primary liver lesions.

Evaluation of response in patients with hepatocellular carcinoma treated with intratumoral dendritic cell vaccination using intravoxel incoherent motion (IVIM) MRI and histogram analysis

BACKGROUND

Immunotherapy of hepatocellular carcinoma (HCC) is an emerging method with promising results. Immunotherapy can have an antitumor effect without affecting tumor size, calling for functional imaging methods for response evaluation.

PURPOSE

To evaluate the response to intratumoral injections with the immune primer ilixadencel in HCCs with diffusion-weighted magnetic resonance imaging (DW-MRI) using intravoxel incoherent motion (IVIM) and histogram analysis.

MATERIAL AND METHODS

A total of 17 patients with advanced HCC were treated with intratumoral injections with ilixadencel on three occasions 2-5 weeks apart. The patients were examined with IVIM before each injection as well as approximately three months after the first injection.

RESULTS

The 10th percentile of perfusion-related parameter D* decreased significantly after the first and second intratumoral injections of ilixadencel compared to baseline (P < 0.05). There was a non-significant trend of lower median region of interest f (perfusion fraction) before injection 2 compared to baseline (P = 0.07). There were significant correlations between the 10th percentile and median of D at baseline and change in tumor size after three months (r = 0.79, P < 0.01 and r = 0.72, P < 0.05, respectively).

CONCLUSION

DW-MRI with IVIM and histogram analysis revealed significant reductions of D* early after treatment as well as an association between D at baseline and smaller tumor growth at three months. The lower percentiles (10th and 50th) were found more important. Further research is needed to confirm our preliminary findings of reduced perfusion after ilixadencel vaccinations, suggesting a treatment effect on HCC.

Associations between angiogenic factors and intravoxel incoherent motion-derived parameters in diffusion-weighted magnetic resonance imaging of breast cancer

Intravoxel incoherent motion (IVIM) diffusion-weighted magnetic resonance imaging (MRI) can be used to estimate perfusion-related parameters, but these parameters may differ, based on the curve-fitting algorithm used for IVIM. Microvessel density (MVD) and vascular endothelial growth factor (VEGF) status are used as angiogenic factors in breast cancer. We aimed to investigate the relationship between MVD, VEGF, and intravoxel incoherent motion (IVIM)-derived parameters, obtained by 4 curve-fitting algorithms, in patients with invasive breast cancers.This retrospective study investigated IVIM-derived parameters, D (ie, tissue diffusivity), D∗ (ie, pseudodiffusivity), and f (ie, perfusion fraction), of 55 breast cancers, using 10 b values (range, 0-800 s/mm2) and 4 curve-fitting algorithms: algorithm 1, linear fitting of D and f first, followed by D∗; algorithm 2, linear fitting of D and f and nonlinear fitting of D∗; algorithm 3, linear fitting of D and f, linear fitting of D∗, and ignoring D contribution for low b values; and algorithm 4, full nonlinear fitting of D, f, and D∗. We evaluated whole-tumor histograms of D, f, and D∗ for their association with MVD and VEGF.D∗10, D∗25, D∗50, D∗mean, D∗75, D∗90, f10, and f25, derived using algorithm 3, were associated with VEGF expression (P = .043, P = 0.012, P = .019, P = .024, P = .044, P = .041, P = .010, and P = .005, respectively). However, no correlation existed between MVD and IVIM-derived parameters.Perfusion-related IVIM parameters obtained by curve-fitting algorithm 3 may reflect VEGF expression.

Usefulness of diffusion derived vessel density computed from a simplified IVIM imaging protocol: An experimental study with rat biliary duct blockage induced liver fibrosis

OBJECTIVES

Liver vessel density can be evaluated by DDVD (diffusion derived vessel density): DDVD_(b0b1) = Sb0/ROIarea0 - Sb1/ROIarea1, where Sb0 and Sb1 refer to the liver signal when b is 0 or 1 s/mm². Sb1 and ROIarea1 may be replaced by other b-values. With a rat biliary duct ligation (BDL) model, this study assesses the usefulness of liver DDVD computed from a simplified IVIM imaging protocol using b = 25 and b = 50 to replace b = 1 s/mm², alone and in combination with other IVIM parameters.

METHODS

Male Sprague-Dawley rats were used. The rat number was 5, 5, 5, and 3 respectively, for the timepoints of 7, 14, 21, 28 days post-BDL surgery. 12 rats had partial biliary duct recanalization performed after the rats had BDL for 7 days and then again followed-up for a mean of 14 days. Liver diffusion MRIs were acquired at 3.0 T with a b-value distribution of 0, 25, 50, 75, 100, 150, 300, 700, 1000 s/mm². DDVD_mean (control rats n = 6) was the mean of DDVD_(b0b25) and DDVD_(b0b50). IVIM fitting started from b = 0 s/mm² with segmented fitting and a threshold b of 50 s/mm² (n = 5 for control rats). Three 3-D spaces were constructed using a combination of the four diffusion parameters.

RESULTS

The control rats and BDL rats (n = 18) had a liver DDVD_mean of 84.0 ± 26.2 and 44.7 ± 14.4 au/pixel (p < 0.001). All 3-D spaces totally separated healthy livers and all fibrotic livers (n = 30, BDL rats and recanalization rats). The mean relative distance between healthy liver cluster and fibrotic liver cluster was 0.331 for PF, D_slow, and D_fast; 0.381 for PF, D_fast, and DDVD_mean; and 0.384 for PF, D_slow, and DDVD_mean.

CONCLUSION

A combination of PF, D_slow, and D_fast allows total separation of healthy livers and fibrotic livers and the integration of DDVD improved the separation.

Improved unsupervised physics-informed deep learning for intravoxel incoherent motion modeling and evaluation in pancreatic cancer patients

PURPOSE

Earlier work showed that IVIM-NETorig , an unsupervised physics-informed deep neural network, was faster and more accurate than other state-of-the-art intravoxel-incoherent motion (IVIM) fitting approaches to diffusion-weighted imaging (DWI). This study presents a substantially improved version, IVIM-NEToptim , and characterizes its superior performance in pancreatic cancer patients.

METHOD

In simulations (signal-to-noise ratio [SNR] = 20), the accuracy, independence, and consistency of IVIM-NET were evaluated for combinations of hyperparameters (fit S0, constraints, network architecture, number of hidden layers, dropout, batch normalization, learning rate), by calculating the normalized root-mean-square error (NRMSE), Spearman's ρ, and the coefficient of variation (CVNET ), respectively. The best performing network, IVIM-NEToptim was compared to least squares (LS) and a Bayesian approach at different SNRs. IVIM-NEToptim 's performance was evaluated in an independent dataset of 23 patients with pancreatic ductal adenocarcinoma. Fourteen of the patients received no treatment between two repeated scan sessions and nine received chemoradiotherapy between the repeated sessions. Intersession within-subject standard deviations (wSD) and treatment-induced changes were assessed.

RESULTS

In simulations (SNR = 20), IVIM-NEToptim outperformed IVIM-NETorig in accuracy (NRMSE(D) = 0.177 vs 0.196; NMRSE(f) = 0.220 vs 0.267; NMRSE(D*) = 0.386 vs 0.393), independence (ρ(D*, f) = 0.22 vs 0.74), and consistency (CVNET (D) = 0.013 vs 0.104; CVNET (f) = 0.020 vs 0.054; CVNET (D*) = 0.036 vs 0.110). IVIM-NEToptim showed superior performance to the LS and Bayesian approaches at SNRs < 50. In vivo, IVIM-NEToptim showed significantly less noisy parameter maps with lower wSD for D and f than the alternatives. In the treated cohort, IVIM-NEToptim detected the most individual patients with significant parameter changes compared to day-to-day variations.

CONCLUSION

IVIM-NEToptim is recommended for accurate, informative, and consistent IVIM fitting to DWI data.

Echo time dependence of biexponential and triexponential intravoxel incoherent motion parameters in the liver

PURPOSE

Intravoxel incoherent motion (IVIM) studies are performed with different acquisition protocols. Comparing them requires knowledge of echo time (TE) dependencies. The TE-dependence of the biexponential perfusion fraction f is well-documented, unlike that of its triexponential counterparts f₁ and f₂ and the biexponential and triexponential pseudodiffusion coefficients D^* , D 1 ∗ , and D 2 ∗ . The purpose was to investigate the TE-dependence of these parameters and to check whether the triexponential pseudodiffusion compartments are associated with arterial and venous blood.

METHODS

Fifteen healthy volunteers (19-58 y; mean: 24.7 y) underwent diffusion-weighted imaging of the abdomen with 24 b-values (0.2-800 s/mm² ) at TEs of 45, 60, 75, and 90 ms. Regions of interest (ROIs) were manually drawn in the liver. One set of bi- and triexponential IVIM parameters per volunteer and TE was determined. The TE-dependence was assessed with the Kruskal-Wallis test.

RESULTS

TE-dependence was observed for f (P < .001), f₁ (P = .001), and f₂ (P < .001). Their median values at the four measured TEs were: f: 0.198/0.240/0.274/0.359, f₁ : 0.113/0.139/0.146/0.205, f₂ : 0.115/0.155/0.182/0.194. D, D^* , D 1 ∗ , and D 2 ∗ showed no significant TE-dependence. Their values were: diffusion coefficient D (10^-4 mm² /s): 9.45/9.63/9.75/9.41, biexponential D^* (10^-2 mm² /s): 5.26/5.52/6.13/5.82, triexponential D 1 ∗ (10^-2 mm² /s): 1.73/2.91/2.25/2.51, triexponential D 2 ∗ (mm² /s): 0.478/1.385/0.616/0.846.

CONCLUSION

f₁ and f₂ show similar TE-dependence as f, ie, increase with rising TE; an effect that must be accounted for when comparing different studies. The diffusion and pseudodiffusion coefficients might be compared without TE correction. Because of the similar TE-dependence of f₁ and f₂ , the triexponential pseudodiffusion compartments are most probably not associated to venous and arterial blood.

Continuous diffusion spectrum computation for diffusion-weighted magnetic resonance imaging of the kidney tubule system

Background

The use of rigid multi-exponential models (with a priori predefined numbers of components) is common practice for diffusion-weighted MRI (DWI) analysis of the kidney. This approach may not accurately reflect renal microstructure, as the data are forced to conform to the a priori assumptions of simplified models. This work examines the feasibility of less constrained, data-driven non-negative least squares (NNLS) continuum modelling for DWI of the kidney tubule system in simulations that include emulations of pathophysiological conditions.

Methods

Non-linear least squares (LS) fitting was used as reference for the simulations. For performance assessment, a threshold of 5% or 10% for the mean absolute percentage error (MAPE) of NNLS and LS results was used. As ground truth, a tri-exponential model using defined volume fractions and diffusion coefficients for each renal compartment (tubule system: D_tubules , f_tubules ; renal tissue: D_tissue , f_tissue ; renal blood: D_blood , f_blood ;) was applied. The impact of: (I) signal-to-noise ratio (SNR) =40-1,000, (II) number of b-values (n=10-50), (III) diffusion weighting (b-range_small =0-800 up to b-range_large =0-2,180 s/mm²), and (IV) fixation of the diffusion coefficients D_tissue and D_blood was examined. NNLS was evaluated for baseline and pathophysiological conditions, namely increased tubular volume fraction (ITV) and renal fibrosis (10%: grade I, mild) and 30% (grade II, moderate).

Results

NNLS showed the same high degree of reliability as the non-linear LS. MAPE of the tubular volume fraction (f_tubules ) decreased with increasing SNR. Increasing the number of b-values was beneficial for f_tubules precision. Using the b-range_large led to a decrease in MAPE _ftubules compared to b-range_small. The use of a medium b-value range of b=0-1,380 s/mm² improved f_tubules precision, and further b_max increases beyond this range yielded diminishing improvements. Fixing D_blood and D_tissue significantly reduced MAPE _ftubules and provided near perfect distinction between baseline and ITV conditions. Without constraining the number of renal compartments in advance, NNLS was able to detect the (fourth) fibrotic compartment, to differentiate it from the other three diffusion components, and to distinguish between 10% vs. 30% fibrosis.

Conclusions

This work demonstrates the feasibility of NNLS modelling for DWI of the kidney tubule system and shows its potential for examining diffusion compartments associated with renal pathophysiology including ITV fraction and different degrees of fibrosis.

Bayesian inference using hierarchical and spatial priors for intravoxel incoherent motion MR imaging in the brain: Analysis of cancer and acute stroke

The intravoxel incoherent motion (IVIM) model allows to map diffusion (D) and perfusion-related parameters (F and D*). Parameter estimation is, however, error-prone due to the non-linearity of the signal model, the limited signal-to-noise ratio (SNR) and the small volume fraction of perfusion in the in-vivo brain. In the present work, the performance of Bayesian inference was examined in the presence of brain pathologies characterized by hypo- and hyperperfusion. In particular, a hierarchical and a spatial prior were combined. Performance was compared relative to conventional segmented least squares regression, hierarchical prior only (non-segmented and segmented data likelihoods) and a deep learning approach. Realistic numerical brain IVIM simulations were conducted to assess errors relative to ground truth. In-vivo, data of 11 central nervous system cancer patients and 9 patients with acute stroke were acquired. The proposed method yielded reduced error in simulations for both the cancer and acute stroke scenarios compared to other methods across the whole investigated SNR range. The contrast-to-noise ratio of the proposed method was better or on par compared to the other techniques in-vivo. The proposed Bayesian approach hence improves IVIM parameter estimation in brain cancer and acute stroke.

Assessment of Imaging Modalities Against Liver Biopsy in Nonalcoholic Fatty Liver Disease: The Amsterdam NAFLD-NASH Cohort

Background

Noninvasive diagnostic methods are urgently required in disease stratification and monitoring in nonalcoholic fatty liver disease (NAFLD). Multiparametric magnetic resonance imaging (MRI) is a promising technique to assess hepatic steatosis, inflammation, and fibrosis, potentially enabling noninvasive identification of individuals with active and advanced stages of NAFLD.

Purpose

To examine the diagnostic performance of multiparametric MRI for the assessment of disease severity along the NAFLD disease spectrum with comparison to histological scores.

Study Type

Prospective, cohort.

Population

Thirty‐seven patients with NAFLD.

Field Strength/Sequence

Multiparametric MRI at 3.0 T consisted of magnetic resonance (MR) spectroscopy (MRS) with multi‐echo stimulated‐echo acquisition mode, magnitude‐based and three‐point Dixon using a two‐dimensional multi‐echo gradient echo, MR elastography (MRE) using a generalized multishot gradient‐recalled echo sequence and intravoxel incoherent motion (IVIM) using a multislice diffusion weighted single‐shot echo‐planar sequence.

Assessment

Histological steatosis grades were compared to proton density fat fraction measured by MRS (PDFF_MRS), magnitude‐based MRI (PDFF_MRI‐M), and three‐point Dixon (PDFF_Dixon), as well as FibroScan® controlled attenuation parameter (CAP). Fibrosis and disease activity were compared to IVIM and MRE. FibroScan® liver stiffness measurements were compared to fibrosis levels. Diagnostic performance of all imaging parameters was determined for distinction between simple steatosis and nonalcoholic steatohepatitis (NASH).

Statistical Tests

Spearman's rank test, Kruskal–Wallis test, Dunn's post‐hoc test with Holm‐Bonferroni P‐value adjustment, receiver operating characteristic curve analysis. A P‐value <0.05 was considered statistically significant.

Results

Histological steatosis grade correlated significantly with PDFF_MRS (r_s = 0.66, P < 0.001), PDFF_MRI‐M (r_s = 0.68, P < 0.001), and PDFF_Dixon (r_s = 0.67, P < 0.001), whereas no correlation was found with CAP. MRE and IVIM diffusion and perfusion significantly correlated with disease activity (r_s = 0.55, P < 0.001, r_s = −0.40, P = 0.016, r_s = −0.37, P = 0.027, respectively) and fibrosis (r_s = 0.55, P < 0.001, r_s = −0.46, P = 0.0051; r_s = −0.53, P < 0.001, respectively). MRE and IVIM diffusion had the highest area‐under‐the‐curve for distinction between simple steatosis and NASH (0.79 and 0.73, respectively).

Data Conclusion

Multiparametric MRI is a promising method for noninvasive, accurate, and sensitive distinction between simple hepatic steatosis and NASH, as well as for the assessment of steatosis and fibrosis severity.

Level of Evidence

2

Technical Efficacy

2

Intravoxel Incoherent Motion Model in Differentiating the Pathological Grades of Esophageal Carcinoma: Comparison of Mono-Exponential and Bi-Exponential Fit Model

PURPOSE

To compare the diagnostic efficiency of the mono-exponential model and bi-exponential model deriving from intravoxel incoherent motion diffusion-weighted imaging (IVIM-DWI) in differentiating the pathological grade of esophageal squamous cell carcinoma (ESCC).

METHODS

Fifty-four patients with ESCC were divided into three groups of poorly-differentiated (PD), moderately-differentiated (MD), and well-differentiated (WD), and underwent the IVIM-DWI scan. Mono-exponential (Dmono, D*mono, and fmono) and bi-exponential fit parameters (Dbi, D*bi, and fbi) were calculated using the IVIM data for the tumors. Mean parameter values of three groups were compared using a one-way ANOVA followed by post hoc tests. The receiver operating characteristic curve was drawn for differentiating pathological grade of ESCC. Correlations between pathological grades and IVIM parameters were analyzed.

RESULTS

There were significant differences in fmono and fbi among the PD, MD and WD ESCC groups (all p<0.05). The fmono were 0.32 ± 0.07, 0.23 ± 0.08, and 0.16 ± 0.05, respectively, and the fbi were 0.35 ± 0.08, 0.26 ± 0.10, and 0.18 ± 0.07, respectively. There was a significant difference in the Dmono between the WD and the PD group (1.48 ± 0.51* 10-3 mm2/s versus 1.05 ± 0.44*10-3 mm2/s, p<0.05), but there was no significant difference between the WD and MD groups, MD and PD groups (all p>0.05). The D*mono, Dbi, and D*bi showed no significant difference among the three groups (all p>0.05). The area under the curve (AUC) of Dmono, fmono and fbi in differentiating WD from PD ESCC were 0.764, 0.961 and 0.932, and the sensitivity and specificity were 92.9% and 60%, 92.9% and 90%, 85.7% and 100%, respectively. The AUC of fmono and fbi in differentiating MD from PD ESCC were 0.839 and 0.757, and the sensitivity and specificity were 78.6% and 80%, 85.7% and 70%, respectively. The AUC of fmono and fbi in differentiating MD from WD ESCC were 0.746 and 0.740, and the sensitivity and specificity were 65% and 85%, 80% and 60%, respectively. The pathologically differentiated grade was correlated with all IVIM parameters (all p<0.05).

CONCLUSIONS

The mono-exponential IVIM model is superior to the bi-exponential IVIM model in differentiating pathological grades of ESCC, which may be a promising imaging method to predict pathological grades of ESCC.

Evidence of Tri-Exponential Decay for Liver Intravoxel Incoherent Motion MRI: A Review of Published Results and Limitations

Diffusion weighted imaging (DWI) and intravoxel incoherent motion (IVIM) have been explored to assess liver tumors and diffused liver diseases. IVIM reflects the microscopic translational motions that occur in voxels in magnetic resonance (MR) DWI. In biologic tissues, molecular diffusion of water and microcirculation of blood in the capillary network can be assessed using IVIM DWI. The most commonly applied model to describe the DWI signal is a bi-exponential model, with a slow compartment of diffusion linked to pure molecular diffusion (represented by the coefficient D_slow), and a fast compartment of diffusion, related to microperfusion (represented by the coefficient D_fast). However, high variance in D_fast estimates has been consistently shown in literature for liver IVIM, restricting its application in clinical practice. This variation could be explained by the presence of another very fast compartment of diffusion in the liver. Therefore, a tri-exponential model would be more suitable to describe the DWI signal. This article reviews the published evidence of the existence of this additional very fast diffusion compartment and discusses the performance and limitations of the tri-exponential model for liver IVIM in current clinical settings.

Quantification of intravoxel incoherent motion with optimized b-values using deep neural network

PURPOSE

To develop a framework for quantifying intravoxel incoherent motion (IVIM) parameters, where a neural network for quantification and b-values for diffusion-weighted imaging are simultaneously optimized.

METHOD

A deep neural network (DNN) method is proposed for accurate quantification of IVIM parameters from multiple diffusion-weighted images. In addition, optimal b-values are selected to acquire the multiple diffusion-weighted images. The proposed framework consists of an MRI signal generation part and an IVIM parameter quantification part. Monte-Carlo (MC) simulations were performed to evaluate the accuracy of the IVIM parameter quantification and the efficacy of b-value optimization. In order to analyze the effect of noise on the optimized b-values, simulations were performed with five different noise levels. For in vivo data, diffusion images were acquired with the b-values from four b-values selection methods for five healthy volunteers at 3T MRI system.

RESULTS

Experiment results showed that both the optimization of b-values and the training of DNN were simultaneously performed to quantify IVIM parameters. We found that the accuracies of the perfusion coefficient (D_p ) and perfusion fraction (f) were more sensitive to b-values than the diffusion coefficient (D) was. Furthermore, when the noise level changed, the optimized b-values also changed. Therefore, noise level has to be considered when optimizing b-values for IVIM quantification.

CONCLUSION

The proposed scheme can simultaneously optimize b-values and train DNN to minimize quantification errors of IVIM parameters. The trained DNN can quantify IVIM parameters from the diffusion-weighted images obtained with the optimized b-values.

Repeatability of IVIM biomarkers from diffusion-weighted MRI in head and neck: Bayesian probability versus neural network

Purpose

The intravoxel incoherent motion (IVIM) model for DWI might provide useful biomarkers for disease management in head and neck cancer. This study compared the repeatability of three IVIM fitting methods to the conventional nonlinear least‐squares regression: Bayesian probability estimation, a recently introduced neural network approach, IVIM‐NET, and a version of the neural network modified to increase consistency, IVIM‐NET_mod.

Methods

Ten healthy volunteers underwent two imaging sessions of the neck, two weeks apart, with two DWI acquisitions per session. Model parameters (ADC, diffusion coefficient Dt, perfusion fraction fp, and pseudo‐diffusion coefficient Dp) from each fit method were determined in the tonsils and in the pterygoid muscles. Within‐subject coefficients of variation (wCV) were calculated to assess repeatability. Training of the neural network was repeated 100 times with random initialization to investigate consistency, quantified by the coefficient of variance.

Results

The Bayesian and neural network approaches outperformed nonlinear regression in terms of wCV. Intersession wCV of Dt in the tonsils was 23.4% for nonlinear regression, 9.7% for Bayesian estimation, 9.4% for IVIM‐NET, and 11.2% for IVIM‐NET_mod. However, results from repeated training of the neural network on the same data set showed differences in parameter estimates: The coefficient of variances over the 100 repetitions for IVIM‐NET were 15% for both Dt and fp, and 94% for Dp; for IVIM‐NET_mod, these values improved to 5%, 9%, and 62%, respectively.

Conclusion

Repeatabilities from the Bayesian and neural network approaches are superior to that of nonlinear regression for estimating IVIM parameters in the head and neck.

Analysis of Renal Diffusion-Weighted Imaging (DWI) Using Apparent Diffusion Coefficient (ADC) and Intravoxel Incoherent Motion (IVIM) Models

Analysis of renal diffusion-weighted imaging (DWI) data to derive markers of tissue properties requires careful consideration of the type, extent, and limitations of the acquired data. Alongside data quality and general suitability for quantitative analysis, choice of diffusion model, fitting algorithm, and processing steps can have consequences for the precision, accuracy, and reliability of derived diffusion parameters. Here we introduce and discuss important steps for diffusion-weighted image processing, and in particular give example analysis protocols and pseudo-code for analysis using the apparent diffusion coefficient (ADC) and intravoxel incoherent motion (IVIM) models. Following an overview of general principles, we provide details of optional steps, and steps for validation of results. Illustrative examples are provided, together with extensive notes discussing wider context of individual steps, and notes on potential pitfalls.This publication is based upon work from the COST Action PARENCHIMA, a community-driven network funded by the European Cooperation in Science and Technology (COST) program of the European Union, which aims to improve the reproducibility and standardization of renal MRI biomarkers. This analysis protocol chapter is complemented by two separate chapters describing the basic concepts and experimental procedure.

Confounding factors in multi-parametric q-MRI protocol: A study of bone marrow biomarkers at 1.5 T

OBJECT

The MRI tissue characterization of vertebral bone marrow includes the measurement of proton density fat fraction (PDFF), T₁ and T₂* relaxation times of the water and fat components (T_1W, T_1F, T₂*_W, T₂*_F), IVIM diffusion D, perfusion fraction f and pseudo-diffusion coefficient D*. However, the measurement of these vertebral bone marrow biomarkers (VBMBs) is affected with several confounding factors. In the current study, we investigated these confounding factors including the regional variation taking the example of variation between the anterior and posterior area in lumbar vertebrae, B₁ inhomogeneity and the effect of fat suppression on f.

MATERIALS AND METHODS

A fat suppressed diffusion-weighted sequence and two 3D gradient multi-echo sequences were used for the measurements of the seven VBMBs. A turbo flash B₁ map sequence was used to estimate B₁ inhomogeneities and thus, to correct flip angle for T₁ quantification. We introduced a correction to perfusion fraction f measured with fat suppression, namely f_PDFF.

RESULTS

A significant difference in the values of PDFF, f and f_PDFF, T_1F, T₂*_W and D was observed between the anterior and posterior region. Although, little variations of flip angle were observed in this anterior-posterior direction in one vertebra but larger variations were observed in head-feet direction from L1 to L5 vertebrae.

DISCUSSION

The regional difference in PDFF, f_PDFF and T₂*_W can be ascribed to differences in the trabecular bone density and vascular network within vertebrae. The regional variation of VBMBs shows that care should be taken in reproducing the same region-of-interest location along a longitudinal study. The same attention should be taken while measuring f in fatty environment, and measuring T₁. Furthermore, the MRI-protocol presented here allows for measurements of seven VBMBs in less than 6 min and is of interest for longitudinal studies of bone marrow diseases.

4D-MRI driven MR-guided online adaptive radiotherapy for abdominal stereotactic body radiation therapy on a high field MR-Linac: Implementation and initial clinical experience

•

This is the first clinical use of 4D-MRI in an online radiation therapy workflow.

•

Eleven patients received SBRT on an Elekta MR-Linac using ATP and ATS workflows.

•

A parallel contour editing approach was successfully utilized in the ATS workflow.

Background and purpose

In this report, we describe our implementation and initial clinical experience using 4D-MRI driven MR-guided online adaptive radiotherapy (MRgOART) for abdominal stereotactic body radiotherapy (SBRT) on the Elekta Unity MR-Linac.

Materials and methods

Eleven patients with abdominal malignancies were treated with free-breathing SBRT in three to five fractions on a 1.5 T MR-Linac. Online adaptive plans were generated using Adapt-To-Position (ATP) or Adapt-To-Shape (ATS) workflows based on motion averaged or mid-position images derived from a pre-beam 4D-MRI. A high performance server positioned on the local MR-Linac machine network was utilized for 4D-MR image reconstruction. A parallel contour editing approach was employed in the ATS workflow. Intravoxel incoherent motion (IVIM) and T2 mapping sequences were acquired during adaptive planning in both ATP and ATS workflows for treatment response monitoring. Adaptive plans were delivered under real-time cine image motion monitoring.

Results

The shortest 4D-MRI time-to-image was the motion averaged image, followed by mid position and respiratory binned images. In this cohert of patients, 50% of treatments utilized the ATS workflow; the remaining treatments utilized the ATP workflow. Mid-position images were utilized as daily planning images for two of the eleven patients. The mean daily adaptive plan secondary dose calculation and ArcCheck 3D Gamma passing rates were 97.5% (92.1–100.0%) and 99.3% (96.2–100.0%), respectively. The median overall treatment times for abdominal SBRT was 46 and 62 min for ATP and ATS workflows, respectively.

Conclusion

We have successfully implemented and utilized a 4D-MRI driven MRgOART process with ATP and ATS workflows for free-breathing abdominal SBRT on a 1.5 T Elekta Unity MR-Linac.

Reliable estimation of brain intravoxel incoherent motion parameters using denoised diffusion-weighted MRI

In this study, we evaluate whether diffusion-weighted magnetic resonance imaging (DW-MRI) data after denoising can provide a reliable estimation of brain intravoxel incoherent motion (IVIM) perfusion parameters. Brain DW-MRI was performed in five healthy volunteers on a 3 T clinical scanner with 12 different b-values ranging from 0 to 1000 s/mm² . DW-MRI data denoised using the proposed method were fitted with a biexponential model to extract perfusion fraction (PF), diffusion coefficient (D) and pseudo-diffusion coefficient (D*). To further evaluate the accuracy and precision of parameter estimation, IVIM parametric images obtained from one volunteer were used to resimulate the DW-MRI data using the biexponential model with the same b-values. Rician noise was added to generate DW-MRI data with various signal-to-noise ratio (SNR) levels. The experimental results showed that the denoised DW-MRI data yielded precise estimates for all IVIM parameters. We also found that IVIM parameters were significantly different between gray matter and white matter (P < 0.05), except for D* (P = 0.6). Our simulation results show that the proposed image denoising method displays good performance in estimating IVIM parameters (both bias and coefficient of variation were <12% for PF, D and D*) in the presence of different levels of simulated Rician noise (SNR_b=0 = 20-40). Simulations and experiments show that brain DW-MRI data after denoising can provide a reliable estimation of IVIM parameters.

Simultaneously spatial and temporal Higher-Order Total Variations for noise suppression and motion reduction in DCE and IVIM

In many applications based on kinetic evaluation analysis and model fitting, quantitative mapping retrieved on data series from modalites such as MRI is completed on a voxel-by-voxel basis, where motion and low signal to noise ratio (SNR) would considerably degenerate the reliability of estimations. The coherence of image series in space and time can be used as prior knowledge to mitigate this occurrence. In this study, spatial and temporal higher-order total variations (HOTVs) are applied on a data series of MRI signal (e.g. dynamic contrast-enhanced (DCE) MRI and intravoxel incoherent motion (IVIM) MRI) to exploit the coherence of signal in space and time to minimize the variabilities caused by motion as well as improving quality of images with low SNR while retaining the physical details of original data properly. Simultaneously applying spatial and temporal HOTVs on images is non-trivial in implementation since it is a non-smooth optimization problem with multiple regularizers. Therefore, we use the proximal gradient method as well as a primal-dual split proximal mechanism to address the problem properly. In addition to increase the reliability of quantitative parametric map estimations, this preprocessing procedure can be included into many existing map estimation algorithms and pipelines effortlessly. We demonstrate our method on the parametric maps estimation for DCE MRI and IVIM MRI.

A novel bayesian approach with conditional autoregressive specification for intravoxel incoherent motion diffusion-weighted MRI

The Intra-Voxel Incoherent Motion (IVIM) model is largely adopted to estimate slow and fast diffusion coefficients of water molecules in biological tissues, which are used in cancer applications. The most reported fitting approach is a voxel-wise segmented non-linear least square, whereas Bayesian approaches with a direct fit, also considering spatial regularization, were proposed too. In this work a novel segmented Bayesian method was proposed, also in combination with a spatial regularization through a Conditional Autoregressive (CAR) prior specification. The two segmented Bayesian approaches, with and without CAR specification, were compared with two standard least-square and a direct Bayesian fitting methods. All approaches were tested on simulated images and real data of patients with head-and-neck and rectal cancer. Estimation accuracy and maps noisiness were quantified on simulated images, whereas the coefficient of variation and the goodness of fit were evaluated for real data. Both versions of the segmented Bayesian approach outperformed the standard methods on simulated images for pseudo-diffusion (D^∗ ) and perfusion fraction (f), whilst the segmented least-square fitting remained the less biased for the diffusion coefficient (D). On real data, Bayesian approaches provided the less noisy maps, and the two Bayesian methods without CAR generally estimated lower values for f and D^∗ coefficients with respect to the other approaches. The proposed segmented Bayesian approaches were superior, in terms of estimation accuracy and maps quality, to the direct Bayesian model and the least-square fittings. The CAR method improved the estimation accuracy, especially for D^∗ .

Optimal acquisition scheme for flow-compensated intravoxel incoherent motion diffusion-weighted imaging in the abdomen: An accurate and precise clinically feasible protocol

PURPOSE

Flow-compensated (FC) diffusion-weighted MRI (DWI) for intravoxel-incoherent motion (IVIM) modeling allows for a more detailed description of tissue microvasculature than conventional IVIM. The long acquisition time of current FC-IVIM protocols, however, has prohibited clinical application. Therefore, we developed an optimized abdominal FC-IVIM acquisition with a clinically feasible scan time.

METHODS

Precision and accuracy of the FC-IVIM parameters were assessed by fitting the FC-IVIM model to signal decay curves, simulated for different acquisition schemes. Diffusion-weighted acquisitions were added subsequently to the protocol, where we chose the combination of b-value, diffusion time and gradient profile (FC or bipolar) that resulted in the largest improvement to its accuracy and precision. The resulting two optimized FC-IVIM protocols with 25 and 50 acquisitions (FC-IVIMopt25 and FC-IVIMopt50 ), together with a complementary acquisition consisting of 50 diffusion-weighting (FC-IVIMcomp ), were acquired in repeated abdominal free-breathing FC-IVIM imaging of seven healthy volunteers. Intersession and intrasession within-subject coefficient of variation of the FC-IVIM parameters were compared for the liver, spleen, and kidneys.

RESULTS

Simulations showed that the performance of FC-IVIM improved in tissue with larger perfusion fraction and signal-to-noise ratio. The scan time of the FC-IVIMopt25 and FC-IVIMopt50 protocols were 8 and 16 min. The best in vivo performance was seen in FC-IVIMopt50 . The intersession within-subject coefficients of variation of FC-IVIMopt50 were 11.6%, 16.3%, 65.5%, and 36.0% for FC-IVIM model parameters diffusivity, perfusion fraction, characteristic time and blood flow velocity, respectively.

CONCLUSIONS

We have optimized the FC-IVIM protocol, allowing for clinically feasible scan times (8-16 min).

On the identification of the blood vessel confounding effect in intravoxel incoherent motion (IVIM) Diffusion-Weighted (DW)-MRI in liver: An efficient sparsity based algorithm

IntraVoxel Incoherent Motion (IVIM) Diffusion-Weighted Magnetic Resonance Imaging (DW-MRI) is of great interest for evaluating tissue diffusion and perfusion and producing parametric maps in clinical applications for liver pathologies. However, the presence of macroscopic blood vessels (not capillaries) in a given Region of Interest (ROI) results in a confounding effect that bias the quantification of tissue perfusion. Therefore, it is necessary to identify those voxels affected by blood vessels. In this paper, an efficient algorithm for an automatic identification of blood vessels in a given ROI is proposed. It relies on the sparsity of the spatial distribution of blood vessels. This sparsity prior can be easily incorporated using the all-voxel IVIM-MRI model introduced in this paper. In addition to the identification of blood vessels, the proposed algorithm provides a quantification of blood vessels, tissue diffusion and tissue perfusion of all voxels in a given ROI, in one single step. Besides, two strategies are proposed in this paper to deal with the nonnegativity of the model parameters. The efficiency of the proposed algorithm compared to the Non-Negative Least Square (NNLS)-based method, recently introduced to deal with the confounding blood vessel effect in the IVIM-MRI model, is confirmed using both realistic and real DW-MR images.

Intravoxel incoherent motion imaging for diagnosing and staging the liver fibrosis and inflammation

PURPOSE

To evaluate the diagnostic accuracy of intravoxel incoherent motion (IVIM) model parameters for the diagnosis and staging of liver fibrosis and inflammation in patients with chronic hepatitis B.

METHODS

Fifty-four patients with chronic hepatitis B and 42 healthy volunteers were included in the study. All subjects were examined by 3 T magnetic resonance imaging. Diffusion-weighted imaging was undertaken with sixteen b values. IVIM parameters [D (true diffusion coefficient), D* (pseudo-diffusion coefficient), f (perfusion fraction)] were calculated. Histological evaluation of biopsy samples was considered the reference standard for the staging of liver fibrosis and inflammation. Differences in IVIM parameters between patient and control groups were analyzed. In the patient group, fibrosis stage and inflammation grade groups were analyzed with respect to IVIM parameters. The correlation was assessed between IVIM parameters and Ishak-modified scale of fibrosis stages and inflammation grades.

RESULTS

The D was significantly lower in the patient group than the control group, p = 0.038 with Cohen's d effect size of 0.452. D was significantly different between fibrosis stage levels. D values decreased in fibrosis stages from the minimal to moderate to marked fibrosis. Fibrosis grades significantly negatively correlated with D and D* values, p = 0.001, and 0.021, respectively. In addition, inflammation grades negatively correlated with f values, p = 0.047.

CONCLUSION

D values measured with IVIM imaging may help to diagnose liver fibrosis. IVIM imaging could be an alternative to liver biopsy for the staging of liver fibrosis.

Comparison of tri-exponential decay versus bi-exponential decay and full fitting versus segmented fitting for modeling liver intravoxel incoherent motion diffusion MRI

OBJECTIVES

To determine whether bi- or tri-exponential models, and full or segmented fittings, better fit the intravoxel incoherent motion (IVIM) imaging signal of healthy livers.

METHODS

Diffusion-weighted images were acquired with a 3 T scanner using a respiratory-triggered echo-planar sequence and 16 b-values (0-800 s/mm² ). Eighteen healthy volunteers had their livers scanned twice in the same session, and then once in another session. Liver parenchyma region-of-interest-based measurements were processed with bi-exponential and tri-exponential models, with both full fitting and segmented fitting (threshold b-value = 200 s/mm² ).

RESULTS

With the signal of all scans averaged, bi-exponential model full fitting showed D_slow  = 1.14 × 10^-3  mm² /s, D_fast  = 193.6 × 10^-3  mm² /s, and perfusion fraction (PF) = 16.9%, and segmented fitting showed D_slow  = 0.98 × 10^-3  mm² /s, D_fast  = 42.2 × 10^-3  mm² /s, and PF = 23.3%. IVIM parameters derived from the tri-exponential model were similar for full fitting and segmented fitting, with slow (D'_slow  = 0.98 × 10^-3  mm² /s; F'_slow  = 76.4 or 76.6%), fast (D'_fast  = 15.1 or 15.4 × 10^-3  mm² /s; F'_fast  = 11.8 or 11.7%) and very fast (D'_Vfast  = 445.0 or 448.8 × 10^-3  mm² /s; F'_Vfast  = 11.8 or 11.7%) diffusion compartments. The tri-exponential model provided an overall better fit than the bi-exponential model. For the bi-exponential model, full fitting provided a better fit at very low and low b-values compared with segmented fitting, with the latter tending to underestimate D_fast ; however, the segmented method demonstrated lower error in signal prediction for high b-values. Compared with full fitting, tri-exponential segmented fitting offered better scan-rescan reproducibility.

CONCLUSION

For healthy liver, tri-exponential modeling is preferred to bi-exponential modeling. For the bi-exponential model, segmented fitting underestimates D_fast , but offers a more accurate estimation of D_slow .

Consensus-based technical recommendations for clinical translation of renal diffusion-weighted MRI

Objectives

Standardization is an important milestone in the validation of DWI-based parameters as imaging biomarkers for renal disease. Here, we propose technical recommendations on three variants of renal DWI, monoexponential DWI, IVIM and DTI, as well as associated MRI biomarkers (ADC, D, D*, f, FA and MD) to aid ongoing international efforts on methodological harmonization.

Materials and methods

Reported DWI biomarkers from 194 prior renal DWI studies were extracted and Pearson correlations between diffusion biomarkers and protocol parameters were computed. Based on the literature review, surveys were designed for the consensus building. Survey data were collected via Delphi consensus process on renal DWI preparation, acquisition, analysis, and reporting. Consensus was defined as ≥ 75% agreement.

Results

Correlations were observed between reported diffusion biomarkers and protocol parameters. Out of 87 survey questions, 57 achieved consensus resolution, while many of the remaining questions were resolved by preference (65–74% agreement). Summary of the literature and survey data as well as recommendations for the preparation, acquisition, processing and reporting of renal DWI were provided.

Discussion

The consensus-based technical recommendations for renal DWI aim to facilitate inter-site harmonization and increase clinical impact of the technique on a larger scale by setting a framework for acquisition protocols for future renal DWI studies. We anticipate an iterative process with continuous updating of the recommendations according to progress in the field.

Electronic supplementary material

The online version of this article (10.1007/s10334-019-00790-y) contains supplementary material, which is available to authorized users.

Optimization of b-value schemes for estimation of the diffusion coefficient and the perfusion fraction with segmented intravoxel incoherent motion model fitting

PURPOSE

Intravoxel incoherent motion (IVIM) modeling for estimation of the diffusion coefficient (D) and perfusion fraction (f) is increasingly popular, but no consensus on standard protocols exists. This study provides a framework for optimization of b-value schemes for reduced estimation uncertainty of D and f from segmented model fitting.

THEORY

Analytical expressions for uncertainties of D and f from segmented model fitting were derived as Cramer-Rao lower bounds (CRLBs).

METHODS

Optimized b-value schemes were obtained for 3 to 12 acquisitions and in the limit of infinitely many acquisitions through constrained minimization of the CRLBs, with b-values constrained to be 0 or 200 to 800 s/mm² . The optimized b-value scheme with eight acquisitions was compared with b-values linearly distributed in the allowed range using simulations and in vivo liver data from seven healthy volunteers.

RESULTS

All optimized b-value schemes contained exactly three unique b-values regardless of the total number of acquisitions (0, 200, and 800 s/mm² ) with repeated acquisitions distributed approximately as 1:2:2. Compared with linearly distributed b-values, the variability of estimates of D and f was reduced by approximately 30% as seen both in simulations and in repeated in vivo measurements.

CONCLUSION

The uncertainty of IVIM D and f estimates can be reduced by the use of optimized b-value schemes.

Principal component analysis fosr fast and model-free denoising of multi b-value diffusion-weighted MR images

Despite the utility of tumour characterisation using quantitative parameter maps from multi-b-value diffusion-weighted MRI (DWI), clinicians often prefer the use of the image with highest diffusion-weighting (b-value), for instance for defining regions of interest (ROIs). However, these images are typically degraded by noise, as they do not utilize the information from the full acquisition. We present a principal component analysis (PCA) approach for model-free denoising of DWI data. PCA-denoising was compared to synthetic MRI, where a diffusion model is fitted for each voxel and a denoised image at a given b-value is generated from the model fit. A quantitative comparison of systematic and random errors was performed on data simulated using several diffusion models (mono-exponential, bi-exponential, stretched-exponential and kurtosis). A qualitative visual comparison was also performed for in vivo images in six healthy volunteers and three pancreatic cancer patients. In simulations, the reduction in random errors from PCA-denoising was substantial (up to 55%) and similar to synthetic MRI (up to 53%). Model-based synthetic MRI denoising resulted in substantial (up to 29% of signal) systematic errors, whereas PCA-denoising was able to denoise without introducing systematic errors (less than 2%). In vivo, the signal-to-noise ratio (SNR) and sharpness of PCA-denoised images were superior to synthetic MRI, resulting in clearer tumour boundaries. In the presence of motion, PCA-denoising did not cause image blurring, unlike image averaging or synthetic MRI. Multi-b-value MRI can be denoised model-free with our PCA-denoising strategy that reduces noise to a level similar to synthetic MRI, but without introducing systematic errors associated with the synthetic MRI method.

Threshold Isocontouring on High b-Value Diffusion-Weighted Images in Magnetic Resonance Mammography

OBJECTIVES

Motivated by the similar appearance of malignant breast lesions in high b-value diffusion-weighted imaging (DWI) and positron emission tomography, the purpose of this work was to evaluate the applicability of a threshold isocontouring approach commonly used in positron emission tomography to analyze DWI data acquired from female human breasts with minimal interobserver variability.

METHODS

Twenty-three female participants (59.4 ± 10.0 years) with 23 lesions initially classified as suggestive of cancers in x-ray mammography screening were subsequently imaged on a 1.5-T magnetic resonance imaging scanner. Diffusion-weighted imaging was performed prior to biopsy with b values of 0, 100, 750, and 1500 s/mm. Isocontouring with different threshold levels was performed on the highest b-value image to determine the voxels used for subsequent evaluation of diffusion metrics. The coefficient of variation was computed by specifying 4 different regions of interest drawn around the lesion. Additionally, a receiver operating statistical analysis was performed.

RESULTS

Using a relative threshold level greater than or equal to 0.85 almost completely suppresses the intra-individual and inter-individual variability. Among 4 studied diffusion metrics, the diffusion coefficients from the intravoxel incoherent motion model returned the highest area under curve value of 0.9. The optimal cut-off diffusivity was found to be 0.85 μm/ms with a sensitivity of 87.5% and specificity of 90.9%.

CONCLUSION

Threshold isocontouring on high b-value maps is a viable approach to reliably evaluate DWI data of suspicious focal lesions in magnetic resonance mammography.