Intravoxel incoherent motion MRI as a means to measure in vivo perfusion: A review of the evidence

Potential predictive value of IVIM MR for xerostomia in nasopharyngeal carcinoma

BACKGROUND AND PURPOSE

Xerostomia, caused by radiation-induced parotid damage, is the most commonly reported radiotherapy (RT) complication for nasopharyngeal carcinoma (NPC). The purpose of this study was to evaluate the value of intravoxel incoherent motion (IVIM) MR in monitoring radiation-induced parotid gland damage and predicting the risk of xerostomia.

METHODS

Fifty-four NPC patients were enrolled and underwent at least three IVIM MR scans: before (pre-RT), after 5 fractions of (5th-RT), halfway through (mid-RT), and after RT (post-RT). The degree of xerostomia patients was assessed before each MR examination. Furthermore, the time when patients first reported xerostomia symptoms was recorded. The changes in IVIM parameters throughout RT, as well as the relationships between IVIM parameters and xerostomia, were analysed.

RESULT

All IVIM parameters increased significantly from pre-RT to post-RT (p < 0.001). The rates of D, D* and f increase increased significantly from pre-RT to mid-RT (p < 0.001), indicating that cell necrosis mainly occurs in the first half of RT. In multivariate analysis, N3 (p = 0.014), pre-D (p = 0.007) and pre-D* (p = 0.003) were independent factors influencing xerostomia. D and f were significantly higher at 5th-RT than at pre-RT (both p < 0.05). IVIM detected parotid gland injury at 5th-RT at an average scanning time of 6.18 ± 1.07 days, earlier than the 11.94 ± 2.61 days when the patient first complained of xerostomia according to the RTOG scale (p < 0.001).

CONCLUSIONS

IVIM MR can dynamically monitor radiation-induced parotid gland damage and assess it earlier and more objectively than RTOG toxicity. Moreover, IVIM can screen people at risk of more severe xerostomia early.

Validation of MRI assessment of foot perfusion for improving treatment of patients with peripheral artery disease

INTRODUCTION

Information on tissue perfusion in the foot is important when treating patients with chronic limb-threatening ischemia. This study aims to test the reliability of different magnetic resonance sequences when measuring perfusion in the foot.

METHODS

Sixteen healthy volunteers had their right foot scanned in a test/retest study with six different magnetic resonance sequences (BOLD, multi-echo gradient echo (mGRE), 2D and 3D pCASL, PASL FAIR, and DWI with intravoxel incoherent motion (IVIM) with quantitative measurements of perfusion. For five sequences, cuff-induced ischemia followed by a hyperactive response was measured. Images of the feet were segmented into angiosomes and perfusion data were extracted from the five angiosomes.

RESULTS

BOLD, PASL FAIR, mGRE, and DWI with IVIM had low mean differences between the first and second scans, while the results of 2D and 3D pCASL had the highest differences. Based on a paired t-test, BOLD, and FAIR were able to distinguish between perfusion and no perfusion in all angiosomes with p-values below 0.01. This was not the case with 2D and 3D pCASL with p-values above 0.05 in all angiosomes. The mGRE could not distinguish between perfusion and no perfusion in the lateral side of the foot.

CONCLUSION

BOLD, mGRE, pASL FAIR, and DWI with IVIM seem to give more robust results compared to 2D and 3D pCASL. Further studies on patients with peripheral artery disease should explore if the sequences can have clinical relevance when assessing tissue ischemia and results of revascularization.

IMPLICATIONS FOR PRACTICE

This study provides knowledge that could be used to improve the diagnosis of patient with chronic limb-threatening ischemia to explore tissue perfusion.

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.

IDEAL-IQ combined with intravoxel incoherent motion diffusion-weighted imaging for quantitative diagnosis of osteoporosis

BACKGROUND

Osteoporosis (OP) is a common chronic metabolic bone disease characterized by decreased bone mineral content and microstructural damage, leading to increased fracture risk. Traditional methods for measuring bone density have limitations in accurately distinguishing vertebral bodies and are influenced by vertebral degeneration and surrounding tissues. Therefore, novel methods are needed to quantitatively assess changes in bone density and improve the accurate diagnosis of OP.

METHODS

This study aimed to explore the applicative value of the iterative decomposition of water and fat with echo asymmetry and least-squares estimation-iron (IDEAL-IQ) sequence combined with intravoxel incoherent motion diffusion-weighted imaging (IVIM-DWI) for the diagnosis of osteoporosis. Data from 135 patients undergoing dual-energy X-ray absorptiometry (DXA), IDEAL-IQ, and IVIM-DWI were prospectively collected and analyzed. Various parameters obtained from IVIM-DWI and IDEAL-IQ sequences were compared, and their diagnostic efficacy was evaluated.

RESULTS

Statistically significant differences were observed among the three groups for FF, R2*, f, D, DDC values, and BMD values. FF and f values exhibited negative correlations with BMD values, with r=-0.313 and - 0.274, respectively, while R2*, D, and DDC values showed positive correlations with BMD values, with r = 0.327, 0.532, and 0.390, respectively. Among these parameters, D demonstrated the highest diagnostic efficacy for osteoporosis (AUC = 0.826), followed by FF (AUC = 0.713). D* exhibited the lowest diagnostic performance for distinguishing the osteoporosis group from the other two groups. Only D showed a significant difference between genders. The AUCs for IDEAL-IQ, IVIM-DWI, and their combination were 0.74, 0.89, and 0.90, respectively.

CONCLUSIONS

IDEAL-IQ combined with IVIM-DWI provides valuable information for the diagnosis of osteoporosis and offers evidence for clinical decisions. The superior diagnostic performance of IVIM-DWI, particularly the D value, suggests its potential as a more sensitive and accurate method for diagnosing osteoporosis compared to IDEAL-IQ. These findings underscore the importance of integrating advanced imaging techniques into clinical practice for improved osteoporosis management and highlight the need for further research to explore the full clinical implications of these imaging modalities.

MRI for Differentiation between HPV-Positive and HPV-Negative Oropharyngeal Squamous Cell Carcinoma: A Systematic Review

Human papillomavirus (HPV) is an important risk factor for oropharyngeal squamous cell carcinoma (OPSCC). HPV-positive (HPV+) cases are associated with a different pathophysiology, microstructure, and prognosis compared to HPV-negative (HPV-) cases. This review aimed to investigate the potential of magnetic resonance imaging (MRI) to discriminate between HPV+ and HPV- tumours and predict HPV status in OPSCC patients. A systematic literature search was performed on 15 December 2022 on EMBASE, MEDLINE ALL, Web of Science, and Cochrane according to PRISMA guidelines. Twenty-eight studies (n = 2634 patients) were included. Five, nineteen, and seven studies investigated structural MRI (e.g., T1, T2-weighted), diffusion-weighted MRI, and other sequences, respectively. Three out of four studies found that HPV+ tumours were significantly smaller in size, and their lymph node metastases were more cystic in structure than HPV- ones. Eleven out of thirteen studies found that the mean apparent diffusion coefficient was significantly higher in HPV- than HPV+ primary tumours. Other sequences need further investigation. Fourteen studies used MRI to predict HPV status using clinical, radiological, and radiomics features. The reported areas under the curve (AUC) values ranged between 0.697 and 0.944. MRI can potentially be used to find differences between HPV+ and HPV- OPSCC patients and predict HPV status with reasonable accuracy. Larger studies with external model validation using independent datasets are needed before clinical implementation.

Intravoxel incoherent motion as a surrogate marker of perfused vascular density in rat brain

Intravoxel incoherent motion (IVIM) MRI has emerged as a valuable technique for the assessment of tissue characteristics and perfusion. However, there is limited knowledge about the relationship between IVIM-derived measures and changes at the level of the vascular network. In this study, we investigated the potential use of IVIM MRI as a noninvasive tool for measuring changes in cerebral vascular density. Variations in quantitative immunohistochemical measurements of the vascular density across different regions in the rat brain (cortex, corpus callosum, hippocampus, thalamus, and hypothalamus) were related to the pseudo-diffusion coefficient D* and the flowing blood fraction f in healthy Wistar rats. We assessed whether region-wise differences in the vascular density are reflected by variations in the IVIM measurements and found a significant positive relationship with the pseudo-diffusion coefficient (p < 0.05, β = 0.24). The effect of cerebrovascular alterations, such as blood-brain barrier (BBB) disruption on the perfusion-related IVIM parameters, is not well understood. Therefore, we investigated the effect of BBB disruption on the IVIM measures in a rat model of metabolic and vascular comorbidities (ZSF1 obese rat) and assessed whether this affects the relationship between the cerebral vascular density and the noninvasive IVIM measurements. We observed increased vascular permeability without detecting any differences in diffusivity, suggesting that BBB leakage is present before changes in the tissue integrity. We observed no significant difference in the relationship between cerebral vascular density and the IVIM measurements in our model of comorbidities compared with healthy normotensive rats.

Magnetic Resonance Perfusion Imaging of Prostate

Magnetic resonance (MR) perfusion imaging, both with and without exogenous contrast agents, has the potential to assess tissue perfusion and vascularity in prostate cancer. Dynamic contrast-enhanced (DCE) MRI is an important element of the clinical non-invasive multiparametric MRI, which can be used to differentiate benign from malignant lesions, to stage tumors, and to monitor response to therapy. The arterial spin labeled (ASL) and intravoxel incoherent motion (IVIM) diffusion-weighted MRI have the advantage of quantitative perfusion measurements without the concerns of gadolinium-based contrast agent safety and retention issues. The adoption of these non-contrast techniques in clinical practice needs more research and clinical evaluation.

Clinical Interpretation of Intravoxel Incoherent Motion Perfusion Imaging in the Brain

Intravoxel incoherent motion (IVIM) perfusion imaging extracts information on blood motion in biological tissue from diffusion-weighted MR images. The method is attractive from a clinical stand point, because it measures in essence local quantitative perfusion, without intravenous contrast injection. Currently, the clinical interpretation of IVIM perfusion maps focuses on the IVIM perfusion fraction maps, but improvements in image quality of the IVIM pseudo-diffusion maps, using advanced postprocessing tools involving artificial intelligence, could lead to an increased interest in this parameters, as it could provide additional local perfusion information in the clinical setting, not otherwise available with other perfusion techniques.

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

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

Intravoxel incoherent motion and diffusion kurtosis imaging at 3T MRI: Application to ischemic stroke

BACKGROUND AND PURPOSE

The DKI-IVIM model that incorporates DKI (diffusional kurtosis imaging) into the IVIM (Intravoxel Incoherent Motion) concept was investigated to assess its utility for both enhanced diffusion characterization and perfusion measurements in ischemic stroke at 3 T.

METHODS

Fifteen stroke patients (71 ± 11 years old) were enrolled and DKI-IVIM analysis was performed using 9 b-values from 0 to 1500 s/mm² chosen with the Cramer-Rao-Lower-Bound optimization approach. Pseudo-diffusion coefficient D*, perfusion fraction f, blood flow-related parameter fD*, the diffusion coefficient D and an additional parameter, the kurtosis, K were determined in the ischemic lesion and controlateral normal tissue based on a region of interest approach. The apparent diffusion coefficient (ADC) and arterial spin labelling (ASL) cerebral blood flow (CBF) parameters were also assessed and parametric maps were obtained for all parameters.

RESULTS

Significant differences were observed for all diffusion parameters with a significant decrease for D (p < 0.0001), ADC (p < 0.0001), and a significant increase for K (p < 0.0001) in the ischemic lesions of all patients. f decreased significantly in these regions (p = 0.0002). The fD* increase was not significant (p = 0.56). The same significant differences were found with a motion correction except for fD* (p = 0.47). CBF significantly decreased in the lesions. ADC was significantly positively correlated with D (p < 0.0001) and negatively with K (p = 0.0002); K was also negatively significantly correlated with D (p = 0.01).

CONCLUSIONS

DKI-IVIM model enables for simultaneous cerebral perfusion and enhanced diffusion characterization in an acceptable clinically acquisition time for the ischemic stroke diagnosis with the additional kurtosis factor estimation, that may better reflect the microstructure heterogeneity.

Comprehensive Brain Tumour Characterisation with VERDICT-MRI: Evaluation of Cellular and Vascular Measures Validated by Histology

The aim of this work was to extend the VERDICT-MRI framework for modelling brain tumours, enabling comprehensive characterisation of both intra- and peritumoural areas with a particular focus on cellular and vascular features. Diffusion MRI data were acquired with multiple b-values (ranging from 50 to 3500 s/mm2), diffusion times, and echo times in 21 patients with brain tumours of different types and with a wide range of cellular and vascular features. We fitted a selection of diffusion models that resulted from the combination of different types of intracellular, extracellular, and vascular compartments to the signal. We compared the models using criteria for parsimony while aiming at good characterisation of all of the key histological brain tumour components. Finally, we evaluated the parameters of the best-performing model in the differentiation of tumour histotypes, using ADC (Apparent Diffusion Coefficient) as a clinical standard reference, and compared them to histopathology and relevant perfusion MRI metrics. The best-performing model for VERDICT in brain tumours was a three-compartment model accounting for anisotropically hindered and isotropically restricted diffusion and isotropic pseudo-diffusion. VERDICT metrics were compatible with the histological appearance of low-grade gliomas and metastases and reflected differences found by histopathology between multiple biopsy samples within tumours. The comparison between histotypes showed that both the intracellular and vascular fractions tended to be higher in tumours with high cellularity (glioblastoma and metastasis), and quantitative analysis showed a trend toward higher values of the intracellular fraction (fic) within the tumour core with increasing glioma grade. We also observed a trend towards a higher free water fraction in vasogenic oedemas around metastases compared to infiltrative oedemas around glioblastomas and WHO 3 gliomas as well as the periphery of low-grade gliomas. In conclusion, we developed and evaluated a multi-compartment diffusion MRI model for brain tumours based on the VERDICT framework, which showed agreement between non-invasive microstructural estimates and histology and encouraging trends for the differentiation of tumour types and sub-regions.

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.

Placental perfusion using intravoxel incoherent motion MRI combined with Doppler findings in differentiating between very low birth weight infants and small for gestational age infants

INTRODUCTION

Abnormalities in placental morphology and function can lead to small for gestational age infants (SGA) and very low birth weight infants (VLBWI). In this study, we explored the value of intravoxel incoherent motion (IVIM) histogram parameters, MRI morphological parameters, and Doppler findings of the placenta in differentiating between VLBWI and SGA.

METHODS

Thirty-three pregnant women who were diagnosed with SGA and met the inclusion criteria were enrolled in this retrospective study and divided into two groups: 22 with non-VLBWI and 11 with VLBWI. The IVIM histogram parameters (perfusion fraction (f), true diffusion coefficient (D), and pseudo-diffusion coefficient (D*)), MRI morphological parameters, and Doppler findings were compared between groups. The diagnostic efficiency was compared by receiver operating characteristic (ROC) curve analysis.

RESULTS

The D_mean, D_90th, D*_90th, f_max, and placental area of the VLBWI group were significantly lower than those of the non-VLBWI group (p < 0.05). The values of umbilical artery pulsatility index, resistance index (RI), and peak systolic velocity/end-diastolic velocity were significantly higher in the VLBWI group than in the non-VLBWI group (p < 0.05). D_90th, placental area, and umbilical artery RI had the highest areas under the curve (AUCs) of the ROC curves, at 0.787, 0.785, and 0.762, respectively. A combined predictive model (D_90th, placental area, and umbilical artery RI) improved the performance in differentiating between VLBWI and SGA compared with the single model (AUC = 0.942).

DISCUSSION

IVIM histogram (D_90th) and MRI morphological (placental area) parameters and a Doppler finding (umbilical artery RI) may be sensitive indicators for differentiating between VLBWI and SGA.

Loss of aquaporin-4 results in glymphatic system dysfunction via brain-wide interstitial fluid stagnation

The glymphatic system is a fluid transport network of cerebrospinal fluid (CSF) entering the brain along arterial perivascular spaces, exchanging with interstitial fluid (ISF), ultimately establishing directional clearance of interstitial solutes. CSF transport is facilitated by the expression of aquaporin-4 (AQP4) water channels on the perivascular endfeet of astrocytes. Mice with genetic deletion of AQP4 (AQP4 KO) exhibit abnormalities in the brain structure and molecular water transport. Yet, no studies have systematically examined how these abnormalities in structure and water transport correlate with glymphatic function. Here, we used high-resolution 3D magnetic resonance (MR) non-contrast cisternography, diffusion-weighted MR imaging (MR-DWI) along with intravoxel-incoherent motion (IVIM) DWI, while evaluating glymphatic function using a standard dynamic contrast-enhanced MR imaging to better understand how water transport and glymphatic function is disrupted after genetic deletion of AQP4. AQP4 KO mice had larger interstitial spaces and total brain volumes resulting in higher water content and reduced CSF space volumes, despite similar CSF production rates and vascular density compared to wildtype mice. The larger interstitial fluid volume likely resulted in increased slow but not fast MR diffusion measures and coincided with reduced glymphatic influx. This markedly altered brain fluid transport in AQP4 KO mice may result from a reduction in glymphatic clearance, leading to enlargement and stagnation of fluid in the interstitial space. Overall, diffusion MR is a useful tool to evaluate glymphatic function and may serve as valuable translational biomarker to study glymphatics in human disease.

Prediction of cerebral hyperperfusion following carotid endarterectomy using intravoxel incoherent motion magnetic resonance imaging

OBJECTIVES

One of the risk factors for cerebral hyperperfusion following carotid endarterectomy (CEA) is a chronic reduction in cerebral perfusion pressure due to internal carotid artery (ICA) stenosis, which is clinically detected as increased cerebral blood volume (CBV). The perfusion fraction (f) is one of the intra-voxel incoherent motion (IVIM) parameters obtained using magnetic resonance (MR) imaging that theoretically reflects CBV. The present study aimed to determine whether preoperative IVIM-f on MR imaging predicts development of cerebral hyperperfusion following CEA.

MATERIALS AND METHODS

Sixty-eight patients with unilateral ICA stenosis (≥ 70%) underwent preoperative diffusion-weighted 3-T MR imaging, and IVIM-f maps were generated from these data. Quantitative brain perfusion single-photon emission computed tomography (SPECT) was performed before and immediately after CEA. Regions-of-interest (ROIs) were automatically placed in the bilateral middle cerebral artery territories in all images using a three-dimensional stereotactic ROI template, and affected-to-contralateral ratios in the ROIs were calculated on IVIM-f maps.

RESULTS

Nine patients (13%) exhibited postoperative hyperperfusion (cerebral blood flow increases of ≥ 100% compared with preoperative values in the ROIs on brain perfusion SPECT). Only high IVIM-f ratios were significantly associated with the occurrence of postoperative hyperperfusion (95% confidence interval, 253.8-6774.2; p = 0.0031) on logistic regression analysis. The sensitivity, specificity, and positive and negative predictive values of the IVIM-f ratio to predict the occurrence of postoperative hyperperfusion were 100%, 81%, 45%, and 100%, respectively.

CONCLUSIONS

Preoperative IVIM-f on MR imaging can predict development of cerebral hyperperfusion following CEA.

Prediction of short-term treatment outcome of nasopharyngeal carcinoma based on voxel incoherent motion imaging and arterial spin labeling quantitative parameters

Purpose

To evaluate the early response of chemoradiotherapy (CRT) in nasopharyngeal carcinoma (NPC) based on intravoxel incoherent motion diffusion-weighted imaging (IVIM-DWI) and three-dimensional pseudo-continuous arterial spin labeling (3D pCASL).

Materials and methods

Forty patients diagnosed with NPC were recruited and divided into complete remission (CR) and partial remission (PR) group after CRT. All patients underwent IVIM and ASL and the related parameters was obtained. These parameters include pure diffusion coefficient (D), pseudo-diffusion coefficient (D*), perfusion fraction (f), average blood flow ( BF_avg), minimum blood flow (BF_min), and maximum blood flow (BF_max). Student's t test was used to compare the difference in ASL and IVIM derived parameters between CR and PR. The Areas under curve (AUC) of the receiver operating characteristic (ROC) was used to analyze the diagnostic performance of each parameter of ASL and IVIM to the treatment outcome.

Results

the D value of IVIM in CR group was lower than that of the PR group ( P = 0.014),. Among the parameters of ASL, the BF_avg and BF_max of the CR group were higher than those of the PR group(p = 0.004,0.013), but the BF_min had no statistical significance in the two groups（P = 0.54）. AUC of D, BF_avg, and BF_max is about 0.731, 0.753, and 0.724, respectively, all of their combined AUC diagnosis was 0.812.

Conclusion

The early response of NPC after CRT can predict by IVIM's diffusion parameters and ASL-related blood flow parameters.

Microvascular Dysfunction Associates With Outcomes in Hypertrophic Cardiomyopathy: Insights From the Intravoxel Incoherent Motion MRI

BACKGROUND

Although intravoxel incoherent motion (IVIM) MRI has emerged as an in vivo marker of tissue diffusion and perfusion, its prognostic value in patients with hypertrophic cardiomyopathy (HCM) remains unclear.

PURPOSE

To investigate whether IVIM-MRI derived parameters are associated with outcomes in patients with HCM.

STUDY TYPE

Prospective cohort.

SUBJECTS

A total of 112 patients (51.72 ± 17.13 years) with suspected or known HCM.

FIELD STRENGTH/SEQUENCE

Single-shot echo planar IVIM imaging, balanced steady-state free precession, and phase-sensitive inversion-recovery late gadolinium enhancement (LGE) sequences at 3 T.

ASSESSMENT

All patients were followed up of 29.3 ± 12.3 months for combined major adverse cardiac events (MACE) including cardiac death, aborted sudden death, heart transplantation, and rehospitalization for heart failure. The CVI42 imaging platform was used to assess morphological and functional MRI indices and to quantify LGE. The Body Diffusion Toolbox was used to derive pseudo diffusion (D*), water molecular diffusion (D) and perfusion fraction (f).

STATISTICAL TESTS

Univariable and stepwise multivariable Cox model analyses were used to investigate the association between variables and composite endpoints. Kaplan-Meier curves were constructed to assess event-free survival, and the event rates were compared by the log-rank test.

RESULTS

A total of 19 patients reached endpoints. Patients with MACE showed a significantly impaired D* value, lower f value, and more extensive LGE than those without MACE (all, P < 0.05), while there was no significant difference in D value (P = 0.285). In the Cox regression models, D* value (hazard ratio [HR] 0.93; 95% CI: 0.88-0.98) and f value (HR 0.65; 95% CI: 0.45-0.92) were independent predictors for MACE. Moreover, in Kaplan-Meier survival analysis, the incidence of MACE was significantly higher in patients with decreased D* value and f value.

CONCLUSIONS

Impaired D* and f values derived from IVIM-MRI are associated with adverse outcomes in patients with HCM.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY STAGE: 2.

What happens to brain outside the thermal ablation zones? An assessment of needle-based therapeutic ultrasound in survival swine

BACKGROUND

In stereotactic radiosurgery, isodose lines must be considered to determine how surrounding tissue is affected. In thermal ablative therapy, such as laser interstitial thermal therapy (LITT), transcranial MR-guided focused ultrasound (tcMRgFUS), and needle-based therapeutic ultrasound (NBTU), how the surrounding area is affected has not been well studied.

OBJECTIVE

We aimed to quantify the transition zone surrounding the ablation core created by magnetic resonance-guided robotically-assisted (MRgRA) delivery of NBTU using multi-slice volumetric 2-D magnetic resonance thermal imaging (MRTI) and subsequent characterization of the resultant tissue damage using histopathologic analysis.

METHODS

Four swine underwent MRgRA NBTU using varying duration and wattage for treatment delivery. Serial MRI images were obtained, and the most representative were overlaid with isodose lines and compared to brain tissue acquired postmortem which underwent histopathologic analysis. These results were also compared to predicted volumes using a finite element analysis model. Contralateral brain tissue was used for control data.

RESULTS

Intraoperative MRTI thermal isodose contours were characterized and comprehensively mapped to post-operative MRI images and qualitatively compared with histological tissue sections postmortem. NBTU 360° ablations induced smaller lesion volumes (33.19 mm³; 120 s, 3 W; 30.05 mm³, 180 s, 4 W) versus 180° ablations (77.20 mm³, 120 s, 3 W; 109.29 mm³; 180 s; 4 W). MRTI/MRI overlay demonstrated the lesion within the proximal isodose lines. The ablation-zone was characterized by dense macrophage infiltration and glial/neuronal loss as demonstrated by glial fibrillary acidic protein (GFAP) and neurofilament (NF) absence and avid CD163 staining. The transition-zone between lesion and normal brain demonstrated decreased macrophage infiltration and measured ∼345 microns (n - 3). We did not detect overt hemorrhages or signs of edema in the adjacent spared tissue.

CONCLUSION

We successfully performed MRgRA NBTU ablation in swine and demonstrated minimal histologic changes extended past the ablation-zone. The lesion was characterized by macrophage infiltration and glial/neuronal loss which decreased through the transition-zone.

Longitudinal Correlations Between Intravoxel Incoherent Motion (IVIM) and Dynamic Contrast-Enhanced (DCE) MRI During Radiotherapy in Prostate Cancer Patients

Purpose

Intravoxel incoherent motion (IVIM) is a promising technique that can acquire perfusion information without the use of contrast agent, contrary to the more established dynamic contrast-enhanced (DCE) technique. This is of interest for treatment response monitoring, where patients can be imaged on each treatment fraction. In this study, longitudinal correlations between IVIM- and DCE parameters were assessed in prostate cancer patients receiving radiation treatment.

Materials and Methods

20 prostate cancer patients were treated on a 1.5 T MR-linac with 20 x 3 or 3.1 Gy. Weekly IVIM and DCE scans were acquired. Tumors, the peripheral zone (PZ), and the transition zone (TZ) were delineated on a T₂-weighted scan acquired on the first fraction. IVIM and DCE scans were registered to this scan and the delineations were propagated. Median values from these delineations were used for further analysis. The IVIM parameters D, f, D^* and the product fD^* were calculated. The Tofts model was used to calculate the DCE parameters K^trans, k_ep and v_e. Pearson correlations were calculated for the IVIM and DCE parameters on values from the first fraction for each region of interest (ROI). For longitudinal analysis, the repeated measures correlation coefficient was used to determine correlations between IVIM and DCE parameters in each ROI.

Results

When averaging over patients, an increase during treatment in all IVIM and DCE parameters was observed in all ROIs, except for D in the PZ and TZ. No significant Pearson correlations were found between any pair of IVIM and DCE parameters measured on the first fraction. Significant but low longitudinal correlations were found for some combinations of IVIM and DCE parameters in the PZ and TZ, while no significant longitudinal correlations were found in the tumor. Notably in the TZ, for both f and fD^*, significant longitudinal correlations with all DCE parameters were found.

Conclusions

The increase in IVIM- and DCE parameters when averaging over patients indicates a measurable response to radiation treatment with both techniques. Although low, significant longitudinal correlations were found which suggests that IVIM could potentially be used as an alternative to DCE for treatment response monitoring.

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.

Value of IVIM in Differential Diagnoses between Benign and Malignant Solitary Lung Nodules and Masses: A Meta-analysis

Purpose

This study aims to evaluate the accuracy of intravoxel incoherent motion diffusion-weighted imaging (IVIM-DWI) in distinguishing malignant and benign solitary pulmonary nodules and masses.

Methods

Studies investigating the diagnostic accuracy of IVIM-DWI in lung lesions published through December 2020 were searched. The standardized mean differences (SMDs) of the apparent diffusion coefficient (ADC), tissue diffusivity (D), pseudo-diffusivity (D*), and perfusion fraction (f) were calculated. The sensitivity, specificity, area under the curve (AUC), publication bias, and heterogeneity were then summarized, and the source of heterogeneity and the reliability of combined results were explored by meta-regression and sensitivity analysis.

Results

A total of 16 studies including 714 malignant and 355 benign lesions were included. Significantly lower ADC, D, and f values were found in malignant pulmonary lesions compared to those in benign lesions. The D value showed the best diagnostic performance (sensitivity = 0.90, specificity = 0.71, AUC = 0.91), followed by ADC (sensitivity = 0.84, specificity = 0.75, AUC = 0.88), f (sensitivity = 0.70, specificity = 0.62, AUC = 0.71), and D^* (sensitivity = 0.67, specificity = 0.61, AUC = 0.67). There was an inconspicuous publication bias in ADC, D, D* and f values, moderate heterogeneity in ADC, and high heterogeneity in D, D*, and f values. Subgroup analysis suggested that both ADC and D values had a significant higher sensitivity in “nodules or masses” than that in “nodules.”

Conclusions

The parameters derived from IVIM-DWI, especially the D value, could further improve the differential diagnosis between malignant and benign solitary pulmonary nodules and masses.

Systematic Review Registration:https://www.crd.york.ac.uk/PROSPERO/#myprospero, identifier: CRD42021226664

Predictors of Outcome after (Chemo)Radiotherapy for Node-Positive Oropharyngeal Cancer: The Role of Functional MRI

The prognosis of a subset of patients with locally advanced oropharyngeal cancer (LA-OPC) is still poor despite improvements in patient selection and treatment. Identifying specific patient- and tumor-related factors can help to select those patients who need intensified treatment. We aimed to assess the role of historical risk factors and novel magnetic resonance imaging (MRI) biomarkers in predicting outcomes in these patients. Patients diagnosed with LA-OPC were studied with diffusion-weighted imaging (DWI) and dynamic-contrast enhanced MRI at baseline and at the 10th radiotherapy (RT) fraction. Clinical information was collected as well. The endpoint of the study was the development of disease progression, locally or distantly. Of the 97 patients enrolled, 68 were eligible for analysis. Disease progression was recorded in 21 patients (11 had loco-regional progression, 10 developed distant metastases). We found a correlation between N diameter and disease control (p = 0.02); features such as p16 status and extranodal extension only showed a trend towards statistical significance. Among perfusion MRI features, higher median values of K_ep both in primary tumor (T, p = 0.016) and lymph node (N, p = 0.003) and lower median values of v_e (p = 0.018 in T, p = 0.004 in N) correlated with better disease control. K_ep P90 and N diameter were identified by MRMR algorithm as the best predictors of outcome. In conclusion, the association of non-invasive MRI biomarkers and patients and tumor characteristics may help in predicting disease behavior and patient outcomes in order to ensure a more customized treatment.

Quantification of Tumor Hypoxia through Unsupervised Modelling of Consumption and Supply Hypoxia MR Imaging in Breast Cancer

Simple Summary

Hypoxia in solid tumors is common in most solid cancers and is associated with treatment resistance to both chemo- and radiation-therapy. There is also reason to believe that hypoxia is an important determinant of metastic disease. Identifying hypoxia in solid tumors is important in treatment planning and decision making. In 2018 Hompland et al. proposed a method, based on quantifying consumption and supply of oxygen from diffusion weighted magnetic resonance imaging, to estimate the hypoxic fraction of a solid tumor. The method was based on training model parameters on a known hypoxia state in prostate cancer. In the present study we verified the validity of the consumption and supply concept in breast cancer. Furthermore, we developed and validated a new approach to the concept that does not require a ground truth to train the parameters.

Abstract

The purpose of the present study is to investigate if consumption and supply hypoxia (CSH) MR-imaging can depict breast cancer hypoxia, using the CSH-method initially developed for prostate cancer. Furthermore, to develop a generalized pan-cancer application of the CSH-method that doesn’t require a hypoxia reference standard for training the CSH-parameters. In a cohort of 69 breast cancer patients, we generated, based on the principles of intravoxel incoherent motion modelling, images reflecting cellular density (apparent diffusion coefficient; ADC) and vascular density (perfusion fraction; fp). Combinations of the information in these images were compared to a molecular hypoxia score made from gene expression data, aiming to identify a way to apply the CSH-methodology in breast cancer. Attempts to adapt previously proposed models for prostate cancer included direct transfers and model parameter rescaling. A novel approach, based on rescaling ADC and fp data to give more nuanced response in the relevant physiologic range, was also introduced. The new CSH-method was validated in a prostate cancer cohort with known hypoxia status. The proposed CSH-method gave estimates of hypoxia that was strongly correlated to the molecular hypoxia score in breast cancer, and hypoxia as measured in pathology slices stained with pimonidazole in prostate cancer. The generalized approach to CSH-imaging depicted hypoxia in both breast and prostate cancers and requires no model training. It is easy to implement using readily available technology and encourages further investigation of CSH-imaging in other cancer entities and in other settings, with the goal being to overcome hypoxia-induced resistance to treatment.

A Clustering Approach to Improve IntraVoxel Incoherent Motion Maps from DW-MRI Using Conditional Auto-Regressive Bayesian Model

The Intra-Voxel Incoherent Motion (IVIM) model allows to estimate water diffusion and perfusion-related coefficients in biological tissues using diffusion weighted MR images. Among the available approaches to fit the IVIM bi-exponential decay, a segmented Bayesian algorithm with a Conditional Auto-Regressive (CAR) prior spatial regularization has been recently proposed to produce more reliable coefficient estimation. However, the CAR spatial regularization can generate inaccurate coefficient estimation, especially at the interfaces between different tissues. To overcome this problem, the segmented CAR model was coupled in this work with a k-means clustering approach, to separate different tissues and exclude voxels from other regions in the CAR prior specification. The proposed approach was compared with the original Bayesian CAR method without clustering and with a state-of-the-art Bayesian approach without CAR. The approaches were tested and compared on simulated images by calculating the estimation error and the coefficient of variation (CV). Furthermore, the proposed method was applied to some illustrative real images of oncologic patients. On simulated images, the proposed innovation reduced the average error of 47%, 21% and 58% for D, f and D*, respectively, compared to the state-of-the-art Bayesian approach, and of 48% and 34% for D and f, respectively, compared to the original CAR, while it achieved the same error for D*. The clustering approach was also able to consistently reduce the CV for each coefficient. On real images, the novel approach did not alter the IVIM maps obtained by the original CAR method, with the advantage of reducing their typical blotchy appearance at the boundaries. The proposed approach represents a valuable improvement over the state-of-the-art Bayesian CAR method and provides more reliable IVIM coefficient estimation, and is less sensitive to bias and inconsistency at tissue/tissue and tissue/background interfaces.

Tumor Microenvironment Modifications Recorded With IVIM Perfusion Analysis and DCE-MRI After Neoadjuvant Radiotherapy: A Preclinical Study

PURPOSE

Neoadjuvant radiotherapy (NeoRT) improves tumor local control and facilitates tumor resection in many cancers. Some clinical studies demonstrated that both timing of surgery and RT schedule influence tumor dissemination, and subsequently patient overall survival. Previously, we developed a pre-clinical model demonstrating the impact of NeoRT schedule and timing of surgery on metastatic spreading. We report on the impact of NeoRT on tumor microenvironment by MRI.

METHODS

According to our NeoRT model, MDA-MB 231 cells were implanted in the flank of SCID mice. Tumors were locally irradiated (PXI X-Rad SmART) with 2x5Gy and then surgically removed at different time points after RT. Diffusion-weighted (DW) and Dynamic contrast enhancement (DCE) MRI images were acquired before RT and every 2 days between RT and surgery. IntraVoxel Incoherent Motion (IVIM) analysis was used to obtain information on intravascular diffusion, related to perfusion (F: perfusion factor) and subsequently tumor vessels perfusion. For DCE-MRI, we performed semi-quantitative analyses.

RESULTS

With this experimental model, a significant and transient increase of the perfusion factor F [50% of the basal value (n=16, p<0.005)] was observed on day 6 after irradiation as well as a significant increase of the WashinSlope with DCE-MRI at day 6 (n=13, p<0.05). Using immunohistochemistry, a significant increase of perfused vessels was highlighted, corresponding to the increase of perfusion in MRI at this same time point. Moreover, Tumor surgical resection during this peak of vascularization results in an increase of metastasis burden (n=10, p<0.05).

CONCLUSION

Significant differences in perfusion-related parameters (F and WashinSlope) were observed on day 6 in a neoadjuvant radiotherapy model using SCID mice. These modifications are correlated with an increase of perfused vessels in histological analysis and also with an increase of metastasis spreading after the surgical procedure. This experimental observation could potentially result in a way to personalize treatment, by modulating the time of surgery guided on MRI functional data, especially tumor perfusion.

The Role of Patient- and Treatment-Related Factors and Early Functional Imaging in Late Radiation-Induced Xerostomia in Oropharyngeal Cancer Patients

Simple Summary

In the present prospective study, we assessed the role of various Magnetic Resonance Imaging biomarkers combined with self-assessed xerostomia questionnaires and patient- and treatment-related factors, in predicting xerostomia at 12 months after chemoradiotherapy for oropharyngeal squamous cell carcinoma. We hypothesized that the integration of pre-treatment imaging biomarkers, which addresses the tissue heterogeneity and individual variations among patients, could improve the accuracy of conventional prediction models that are based only on dose information, ultimately providing a better understanding of the pathophysiological mechanisms underlying radiation induced salivary dysfunction. The implementation of multifactorial models, driven by machine learning algorithms, may improve prediction accuracy of radiation-induced toxicity and tailor individual treatment options for patients.

Abstract

The advent of quantitative imaging in personalized radiotherapy (RT) has offered the opportunity for a better understanding of individual variations in intrinsic radiosensitivity. We aimed to assess the role of magnetic resonance imaging (MRI) biomarkers, patient-related factors, and treatment-related factors in predicting xerostomia 12 months after RT (XER₁₂) in patients affected by oropharyngeal squamous cell carcinoma (OSCC). Patients with locally advanced OSCC underwent diffusion-weighted imaging (DWI) and dynamic-contrast enhanced MRI at baseline; DWI was repeated at the 10th fraction of RT. The Radiation Therapy Oncology Group (RTOG) toxicity scale was used to evaluate salivary gland toxicity. Xerostomia-related questionnaires (XQs) were administered weekly during and after RT. RTOG toxicity ≥ grade 2 at XER₁₂ was considered as endpoint to build prediction models. A Decision Tree classification learner was applied to build the prediction models following a five-fold cross-validation. Of the 89 patients enrolled, 63 were eligible for analysis. Thirty-six (57.1%) and 21 (33.3%) patients developed grade 1 and grade 2 XER₁₂, respectively. Including only baseline variables, the model based on DCE-MRI and V65 (%) (volume of both glands receiving doses ≥ 65 Gy) had a fair accuracy (77%, 95% CI: 66.5–85.4%). The model based on V65 (%) and XQ-Int_mid (integral of acute XQ scores from the start to the middle of RT) reached the best accuracy (81%, 95% CI: 71–88.7%). In conclusion, non-invasive biomarkers from DCE-MRI, in combination with dosimetric variables and self-assessed acute XQ scores during treatment may help predict grade 2 XER₁₂ with a fair to good accuracy.

Monitoring Bevacizumab-Induced Tumor Vascular Normalization by Intravoxel Incoherent Motion Diffusion-Weighted MRI

BACKGROUND

Accurate monitoring of tumor blood vessel normalization progression is beneficial to accurate treatment of patients. At present, there is a lack of safe and noninvasive monitoring methods.

PURPOSE

To serial monitor the vascular normalization time window of tumor antiangiogenesis treatment through intravoxel incoherent motion diffusion-weighted imaging (IVIM-DWI) and histopathological methods.

STUDY TYPE

Exploratory animal study.

POPULATION

Sixty rat C6 glioma models were randomly and equally divided into the control groups (N = 30) and bevacizumab treatment groups (N = 30). Twenty-five for magnetic resonance imaging (MRI) and five for electron microscope testing in each group.

FIELD STRENGTH/SEQUENCE

T1-weighted imaging (T1WI), T2WI with a fast spin echo sequence and IVIM-DWI with a spin-echo echo-planar imaging sequence at 3 T.

ASSESSMENT

IVIM-DWI quantitative parameters (f, D, D*, and fD*) were obtained on days 0, 2, 4, 6, and 8 after bevacizumab treatment. After MRI, the microvessel density (MVD), pericyte coverage, and hypoxia-inducible factor-1α (HIF-1α) were assessed. Electron microscope observation was performed at each time point.

STATISTICAL TESTS

One-way analysis of variance and Student's t-tests were used to compare differences within and between groups. Spearman's correlation coefficient (r) assess the correlation between IVIM and pathological parameters. The intragroup correlation coefficient was determined to assess the repeatability of each IVIM parameter.

RESULTS

The IVIM-DWI perfusion parameters (f and fD*) of the treated group were higher than the control group on days 2 and 4. Compared to the control group, MVD decreased on days 2 and pericyte coverage increased on days 4 in the treatment group. Electron microscopy showed that the tight junctions of the treatment group were prolonged on days 2-4. In the control group, f had the highest correlation with MVD (r = 0.689). In the treated group, f had a good correlation with pericyte coverage (r = 0.557), HIF-1α had a moderately positive correlation with f (r = 0.480) and fD*(r = 0.447).

DATA CONCLUSION

The vascular normalization time window of bevacizumab treatment of glioma was days 2-4 after antiangiogenesis treatment, which could be monitored noninvasively by IVIM-DWI.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 3.

Four functional magnetic resonance imaging techniques for skeletal muscle exploration, a systematic review

BACKGROUND

The study of muscle health has become more relevant lately, due to global aging and a higher incidence of musculoskeletal pathologies. Current exploration techniques, such as electromyography, do not provide accurate spatial information.

OBJECTIVE

The objective of this work is to perform a systematic review of the literature to synthesize the contributions that can offer functional MRI techniques commonly used in neuroimaging, applied to skeletal muscle: Blood Oxygen Level Dependent (BOLD), IntraVoxel Incoherent Motion (IVIM), Arterial Spin Labeling (ASL) and Dynamic Contrast Enhanced (DCE).

EVIDENCE ACQUISITION

Web of Science and Medline databases were searched, over the last 10 years, focused on the use of BOLD, ASL, IVIM or DCE in skeletal muscle.

EVIDENCE SYNTHESIS

59 articles were included after applying the selection criteria. 37 studies were performed in healthy subjects, and 22 in patients with different pathologies: in peripheral arterial disease, systemic sclerosis, diabetes, osteoporosis, adolescent idiopathic scoliosis, and dermatomyositis. Reference values in healthy subjects still vary in some cases.

CONCLUSION

The studies show the feasibility of implementing functional MRI through BOLD, ASL, IVIM or DCE imaging in several muscles and their possible utility in different pathologies. A synthesis of how to implement such exploration is given here.

CLINICAL IMPACT

These four techniques are based on sequences already present in clinical MRI scanners, therefore, their use for functional muscle exploration does not require additional investment. These techniques allow visualization and quantification of parameters associated with the vascular health of the muscles and represent interesting support for musculoskeletal exploration.

A Fast Protocol for Multiparametric Characterisation of Diffusion in the Brain and Brain Tumours

Multi-parametric tissue characterisation is demonstrated using a 4-minute protocol based on diffusion trace acquisitions. Three diffusion regimes are covered simultaneously: pseudo-perfusion, Gaussian, and non-Gaussian diffusion. The clinical utility of this method for fast multi-parametric mapping for brain tumours is explored. A cohort of 17 brain tumour patients was measured on a 3T hybrid MR-PET scanner with a standard clinical MRI protocol, to which the proposed multi-parametric diffusion protocol was subsequently added. For comparison purposes, standard perfusion and a full diffusion kurtosis protocol were acquired. Simultaneous amino-acid (¹⁸F-FET) PET enabled the identification of active tumour tissue. The metrics derived from the proposed protocol included perfusion fraction, pseudo-diffusivity, apparent diffusivity, and apparent kurtosis. These metrics were compared to the corresponding metrics from the dedicated acquisitions: cerebral blood volume and flow, mean diffusivity and mean kurtosis. Simulations were carried out to assess the influence of fitting methods and noise levels on the estimation of the parameters. The diffusion and kurtosis metrics obtained from the proposed protocol show strong to very strong correlations with those derived from the conventional protocol. However, a bias towards lower values was observed. The pseudo-perfusion parameters showed very weak to weak correlations compared to their perfusion counterparts. In conclusion, we introduce a clinically applicable protocol for measuring multiple parameters and demonstrate its relevance to pathological tissue characterisation.

Assessment of Impaired Cerebrovascular Reactivity in Chronic Cerebral Ischemia using Intravoxel Incoherent Motion Magnetic Resonance Imaging

BACKGROUND

The severity of chronic cerebral ischemia can be assessed using cerebrovascular reactivity (CVR) to acetazolamide (ACZ) challenge, which is measured by single-photon emission computed tomography (SPECT); however, this is an invasive method. We investigated whether intravoxel incoherent motion (IVIM) magnetic resonance imaging (MRI) can assess impaired CVR in preoperative patients with chronic cerebral ischemia and compared it to SPECT-CVR.

METHODS

Forty-seven patients with unilateral cervical carotid artery stenosis underwent diffusion-weighted MRI with 11 b-values in the range of 0-800 s/mm² and cerebral perfusion SPECT with the ACZ challenge. The perfusion fraction (f) and diffusion coefficient (D) of the IVIM parameters were calculated using a bi-exponential model. The f and D values and these ratios of the ipsilateral middle cerebral artery territory against the contralateral side were compared with the CVR values of the affected side calculated from the SPECT data.

RESULTS

The IVIM-f and D values in the affected side were significantly higher than those in the unaffected side (median: 7.74% vs. 7.45%, p = 0.027; 0.816 vs. 0.801 10^-3mm²/s, p < 0.001; respectively). However, there were no significant correlations between the f or D values and SPECT-CVR values in the affected side. In contrast, the f ratio showed a moderate negative correlation with the SPECT-CVR values (r = -0.40, p = 0.006) and detected impaired CVR (< 18.4%) with a sensitivity/specificity of 0.71/0.90.

CONCLUSION

The IVIM perfusion parameter, f, can noninvasively assess impaired CVR with high sensitivity and specificity in patients with unilateral cervical carotid artery stenosis.

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/mm2, 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/mm2.

Intravoxel incoherent motion diffusion-weighted imaging in the characterization of Alzheimer's disease

OBJECTIVES

Alzheimer's disease (AD) is the most common type of dementia, and characterizing brain changes in AD is important for clinical diagnosis and prognosis. This study was designed to evaluate the classification performance of intravoxel incoherent motion (IVIM) diffusion-weighted imaging in differentiating between AD patients and normal control (NC) subjects and to explore its potential effectiveness as a neuroimaging biomarker.

METHODS

Thirty-one patients with probable AD and twenty NC subjects were included in the prospective study. IVIM data were subjected to postprocessing, and parameters including the apparent diffusion coefficient (ADC), slow diffusion coefficient (D_s), fast diffusion coefficient (D_f), perfusion fraction (fp) and D_f*fp were calculated. The classification model was developed and confirmed with cross-validation (group A/B) using Support Vector Machine (SVM). Correlations between IVIM parameters and Mini-Mental State Examination (MMSE) scores in AD patients were investigated using partial correlation analysis.

RESULTS

Diffusion MRI revealed significant region-specific differences that aided in differentiating AD patients from controls. Among the analyzed regions and parameters, the D_f of the right precuneus (PreR) (ρ = 0.515; P = 0.006) and the left cerebellum (CL) (ρ = 0.429; P = 0.026) demonstrated significant associations with the cognitive function of AD patients. An area under the receiver operating characteristics curve (AUC) of 0.84 (95% CI: 0.66, 0.99) was calculated for the validation in dataset B after the prediction model was trained on dataset A. When the datasets were reversed, an AUC of 0.90 (95% CI: 0.75, 1.00) was calculated for the validation in dataset A, after the prediction model trained in dataset B.

CONCLUSION

IVIM imaging is a promising method for the classification of AD and NC subjects, and IVIM parameters of precuneus and cerebellum might be effective biomarker for the diagnosis of AD.

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.

A Multi-Variate framework to assess reliability and discrimination power of Bayesian estimation of Intravoxel Incoherent Motion parameters

PURPOSE

To propose a multivariate multi-step framework for a systematic assessment of the estimation reliability and discriminability of Intravoxel Incoherent Motion (IVIM) model parameters.

METHODS

Monte-Carlo simulations were generated on a range of SNRs and in different IVIM combinations considering: i) a dense discretization with 24 b-values; ii) a discretization with 9 b-values. A state-of-the-art Bayesian fitting method was adopted. The framework assessed: i) the best model between mono- and bi-exponential, through the BIC index; ii) the fitting accuracy; iii) the power in discriminating two different IVIM parameters distributions of estimated coefficients, using a multivariate test. Exemplificative oncologic cases were also presented.

RESULTS

The bi-exponential fitting was reliable for perfusion fraction higher than 5%, with high accuracy in D estimation, acceptable error for f, but high uncertainty in D*. The discrimination of two distributions is generally feasible if differences in D values (at least 0.3 x10^-3 mm²/s) are present; in the case of similar D values, a minimal difference of 5% in f can be discriminated just in case of balanced sample size and dense b-values discretization, whereas the impact of D* is quite negligible. These results were also supported by clinical examples.

CONCLUSIONS

IVIM model is generally accurate in estimating diffusion, but uncertainties related to perfusion estimation are not negligible and compromise the discrimination power when different populations should be differentiated. The proposed framework should be adopted as interpretative guidelines to better understand when IVIM model applied on real data can provide reliable findings.

Simulation of intravoxel incoherent perfusion signal using a realistic capillary network of a mouse brain

PURPOSE

To simulate the intravoxel incoherent perfusion magnetic resonance magnitude signal from the motion of blood particles in three realistic vascular network graphs from a mouse brain.

METHODS

In three networks generated from the cortex of a mouse scanned by two-photon laser microscopy, blood flow in each vessel was simulated using Poiseuille's law. The trajectories, flow speeds and phases acquired by a fixed number of simulated blood particles during a Stejskal-Tanner bipolar pulse gradient scheme were computed. The resulting magnitude signal was obtained by integrating all phases and the pseudo-diffusion coefficient D* was estimated by fitting an exponential signal decay. To better understand the anatomical source of the intravoxel incoherent motion (IVIM) perfusion signal, the above was repeated restricting the simulation to various types of vessel.

RESULTS

The characteristics of the three microvascular networks were respectively vessel lengths (mean ± std. dev.) 67.2 ± 53.6 μm, 59.8 ± 46.2 μm and 64.5 ± 50.9 μm, diameters 6.0 ± 3.5 μm, 5.7 ± 3.6 μm and 6.1 ± 3.7 μm and simulated blood velocity 0.9 ± 1.7 μm/ms, 1.4 ± 2.5 μm/ms and 0.7 ± 2.1 μm/ms. Exponential fitting of the simulated signal decay as a function of b-value resulted in the following D*-values [10-3 mm2 /s]: 31.7, 40.4 and 33.4. The signal decay for low b-values was the largest in the larger vessels, but the smaller vessels and the capillaries accounted for more of the total volume of the networks.

CONCLUSION

This simulation improves the theoretical understanding of the IVIM perfusion estimation method by directly linking the MR IVIM perfusion signal to an ultra-high resolution measurement of the microvascular network and a realistic blood flow simulation.

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.

Measuring Perfusion: Intravoxel Incoherent Motion MR Imaging

The signal acquired in vivo using a diffusion-weighted MR imaging (DWI) sequence is influenced by blood motion in the tissue. This means that perfusion information from a DWI sequence can be obtained in addition to thermal diffusion, if the appropriate sequence parameters and postprocessing methods are applied. This is commonly regrouped under the denomination intravoxel incoherent motion (IVIM) perfusion MR imaging. Of relevance, the perfusion information acquired with IVIM is essentially local, quantitative and acquired without intravenous injection of contrast media. The aim of this work is to review the IVIM method and its clinical applications.

Investigating DWI changes in white matter of meningioma patients treated with proton therapy

PURPOSE

To evaluate changes in diffusion and perfusion-related properties of white matter (WM) induced by proton therapy, which is capable of a greater dose sparing to organs at risk with respect to conventional X-ray radiotherapy, and to eventually expose early manifestations of delayed neuro-toxicities.

METHODS

Apparent diffusion coefficient (ADC) and IVIM parameters (D, D* and f) were estimated from diffusion-weighted MRI (DWI) in 46 patients affected by meningioma and treated with proton therapy. The impact on changes in diffusion and perfusion-related WM properties of dose and time, as well as the influence of demographic and pre-treatment clinical information, were investigated through linear mixed-effects models.

RESULTS

Decreasing trends in ADC and D were found for WM regions hit by medium-high (30-40 Gy(RBE)) and high (>40 Gy(RBE)) doses, which are compatible with diffusion restriction due to radiation-induced cellular injury. Significant influence of dose and time on median ADC changes were observed. Also, D* showed a significant dependency on dose, whereas f consistently showed no dependency on dose and time. Age, gender and surgery extent were also found to affect changes in ADC.

CONCLUSIONS

These results overall agree with those from studies conducted on cohorts of mixed proton and X-ray radiotherapy patients. Future work should focus on relating our findings with clinical information of co-morbidities and thus exploiting such or more advanced imaging data to build normal tissue complication probability models to better integrate clinical and dose information.

Quantitative Magnetic Resonance Imaging for Biological Image-Guided Adaptive Radiotherapy

MRI-guided radiotherapy systems have the potential to bring two important concepts in modern radiotherapy together: adaptive radiotherapy and biological targeting. Based on frequent anatomical and functional imaging, monitoring the changes that occur in volume, shape as well as biological characteristics, a treatment plan can be updated regularly to accommodate the observed treatment response. For this purpose, quantitative imaging biomarkers need to be identified that show changes early during treatment and predict treatment outcome. This review provides an overview of the current evidence on quantitative MRI measurements during radiotherapy and their potential as an imaging biomarker on MRI-guided radiotherapy systems.

Associations between IVIM histogram parameters and histopathology in rectal cancer

PURPOSE

Histogram analysis can better reflect tumor heterogeneity than conventional imaging analysis. The present study analyzed possible correlations between histogram analysis parameters derived from Intravoxel-incoherent imaging (IVIM) and histopathological features in rectal cancer (RC).

METHODS

Seventeen patients with histopathologically proven rectal adenocarcinomas were retrospectively acquired. In all cases, pelvic MRI was performed. Diffusion weighted imaging was obtained using a multi-slice single-shot echo-planar imaging sequence with b values of 0, 50, 200, 500 and 1000 s/mm². Simplified IVIM analysis was performed using the IntelliSpace portal, version 10 and the following images were generated: f (perfusion fraction) map, D (true diffusion coefficient) map, and ADC map utilizing all b-values. Histogram based analysis of signal intensities was performed for every IVIM map using an in-house matlab tool. Histopathology was investigated using Ki 67 specimens with calculation of Ki 67-index and cellularity. CD31 stained specimens were used for calculation of microvessel density (MVD).

RESULTS

There were statistically significant correlations between Ki 67 index and mode derived from ADC as well as entropy from f, r=-0.50, p=.04 and r=-0.55, p=.02, respectively. MVD correlated well with parameters derived from f.

CONCLUSION

IVIM histogram analysis parameters can reflect histopathology in RC. ADC and D values are associated with proliferation potential. Perfusion fraction f is associated with MVD.

Differentiating atypical hemangiomas and vertebral metastases: a field-of-view (FOV) and FOCUS intravoxel incoherent motion (IVIM) diffusion-weighted imaging (DWI) study

PURPOSE

Some atypical vertebral hemangiomas (VHs) may mimic metastases on routine MRI and can result in misdiagnosis and ultimately to additional imaging, biopsy and unnecessary costs. The purpose of this study is to assess the utility of intravoxel incoherent motion (IVIM) diffusion-weighted imaging (DWI) on account of field-of-view optimized and constrained undistorted single shot (FOCUS) in distinguishing atypical VHs and vertebral metastases.

METHODS

A total of 25 patients with vertebral metastases and 25 patients with atypical VHs were confirmed by clinical follow-up or pathology. IVIM-DWI imaging was performed at different b values (0, 30, 50, 100, 150, 200, 400, 600, 800, 1000 mm²/s). IVIM parameters [the true diffusion coefficient (D), pseudodiffusion coefficient (D*), standard apparent diffusion coefficient (ADC), and perfusion fraction (f)] were calculated and compared between two groups by using Student's t test. A receiver operating characteristic analysis was performed.

RESULTS

Quantitative analysis of standard ADC and D parameters showed significantly lower values in vertebral metastases when compared to atypical hemangiomas [ADC value: (0.70 ± 0.12) × 10^-3 mm²/s vs (1.14 ± 0.28) × 10^-3 mm²/s; D value: (0.47 ± 0.07) × 10^-3 mm²/s vs (0.76 ± 0.14) × 10^-3 mm²/s, all P < 0.01]. The sensitivity and specificity of D value were 93.8% and 92.3%, respectively.

CONCLUSION

The standard ADC value and D value may be used as an indicator to distinguish vertebral metastases from atypical VHs. FOCUS IVIM-derived parameters provide potential value in the quantitatively differentiating vertebral metastases from vertebral atypical hemangiomas.

Functional diagnosis of placenta accreta by intravoxel incoherent motion model diffusion-weighted imaging

OBJECTIVES

To investigate the diagnostic value of intravoxel incoherent motion (IVIM) DWI for placenta accreta by comparing diffusion and perfusion characteristics of placentas with accreta lesions (APs) with those of normal placentas (NPs).

METHODS

Twenty-five pregnant women with AP and 24 with NP underwent 3-T magnetic resonance examinations with IVIM-DWI. The perfusion percentage (f), pseudo-diffusion coefficient (D*), and diffusion coefficient (D) values were calculated from different ROIs: the entire-plane of the AP (AP-ROI) and NP (NP-ROI) and the implanted (IR-ROI) and non-implanted region (NIR-ROI) of the AP. The AP-ROIs and NP-ROIs were compared using covariance analysis; the IR-ROIs and NIR-ROIs were compared using the Wilcoxon signed-rank test. ROC curves were produced to evaluate the parameters for predicting placenta accreta.

RESULTS

The f and D* values for the AP-ROIs ([45.0 ± 7.63]%, [11.64 ± 2.15]mm²/s) were significantly higher than those for the NP-ROIs ([31.85 ± 5.96]%, [9.04 ± 3.13]mm²/s) (both p < 0.05); the IR-ROIs (54.8%, 14.03 mm²/s) were also significantly higher than the NIR-ROIs (37.4%, 11.4 mm²/s) (both p < 0.05). No significant differences were found between the D values of the AP-ROIs and NP-ROIs (p > 0.05) or of the IR-ROIs and NIR-ROIs (p > 0.05). The areas under the curve for f and D* of the ROC curves were 0.93 and 0.79, respectively.

CONCLUSIONS

These results suggest that the IVIM parameters f and D* can be used to quantitatively evaluate the higher perfusion of AP when compared with NP. Furthermore, IVIM may be a useful functional diagnostic technique to predict placenta accreta.

KEY POINTS

• Intravoxel incoherent motion (IVIM) may be a useful diagnostic technique to quantitatively estimate the perfusion of the placenta. • The perfusion percentage (f) and pseudo-diffusion coefficient (D*) values differed significantly between placentas with accreta lesions and normal placentas. • ROC curves showed that perfusion percentage (f) and pseudo-diffusion coefficient (D*) values could accurately predict placenta accreta.

Progress of intravoxel incoherent motion diffusion-weighted imaging in liver diseases

Traditional magnetic resonance (MR) diffusion-weighted imaging (DWI) uses a single exponential model to obtain the apparent diffusion coefficient to quantitatively reflect the diffusion motion of water molecules in living tissues, but it is affected by blood perfusion. Intravoxel incoherent motion (IVIM)-DWI utilizes a double-exponential model to obtain information on pure water molecule diffusion and microcirculatory perfusion-related diffusion, which compensates for the insufficiency of traditional DWI. In recent years, research on the application of IVIM-DWI in the diagnosis and treatment of hepatic diseases has gradually increased and has achieved considerable progress. This study mainly reviews the basic principles of IVIM-DWI and related research progress in the diagnosis and treatment of hepatic diseases.

Diffusion MR Imaging in the Head and Neck: Principles and Applications

Diffusion imaging is a functional MR imaging tool that creates tissue contrast representative of the random, microscopic translational motion of water molecules within human body tissues. Long considered a cornerstone MR imaging sequence for brain imaging, diffusion-weighted imaging (DWI) increasingly is used for head and neck imaging. This review reports the current state of diffusion techniques for head and neck imaging, including conventional DWI, DWI trace with apparent diffusion coefficient map, diffusion tensor imaging, intravoxel incoherent motion, and diffusion kurtosis imaging. This article describes background physics, reports supportive evidence and potential pitfalls, highlights technical advances, and details practical clinical applications.

Diffusion kurtosis imaging in head and neck cancer: A correlation study with dynamic contrast enhanced MRI

PURPOSE

To investigate the biophysical meaning of Diffusion Kurtosis Imaging (DKI) parameters via correlations with the perfusion parameters obtained from a long Dynamic Contrast Enhanced MRI scan, in head and neck (HN) cancer.

METHODS

Twenty two patients with newly diagnosed HN tumor were included in the present retrospective study. Some patients had multiple lesions, therefore a total of 26 lesions were analyzed. DKI was acquired using 5b values at 0, 500, 1000,1500 and 2000 s/mm². DCE-MRI was obtained with 130 dynamic volumes, with a temporal resolution of 5 s, to achieve a long scan time (>10 min). The apparent diffusion coefficient D_app and apparent diffusional kurtosis K_app were calculated voxel-by-voxel, removing the point at b value = 0 to eliminate possible perfusion effects on the parameter estimations. The transfer constants K^trans and K_ep, v_e, and the histogram-based entropy (En) and interquartile range (IQR) of each DCE-MRI parameter were quantified. Correlations between all variables were investigated by the Spearman's Rho correlation test.

RESULTS

Moderate relationships emerged between D_app and K_ep (Rho =  - 0.510, p = 0.009), and between D_app and v_e (Rho = 0.418, p = 0.038). En(K_ep) was significantly related to K_app (Rho = 0.407, p = 0.043), while IQR(K_ep) showed an inverse association with D_app (Rho = -0.422, p = 0.035).

CONCLUSIONS

A weak to intermediate correlation was found between DKI parameters and both K_ep and v_e. The kurtosis was associated to the intratumoral heterogeneity and complexity of the capillary permeability, expressed by En(K_ep).

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.

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^∗ .

Demonstration of asymmetric muscle perfusion of the back after exercise in patients with adolescent idiopathic scoliosis using intravoxel incoherent motion (IVIM) MRI

The purpose of this work was to quantify muscular perfusion patterns of back muscles after exercise in patients with adolescent idiopathic scoliosis (AIS) using intravoxel incoherent motion (IVIM) MR perfusion imaging. The paraspinal muscles of eight patients with AIS (Cobb angle 35 ± 10°, range [25-47°]) and nine healthy volunteers were scanned with a 1.5 T MRI, at rest and after performing a symmetric back muscle exercise on a Roman chair. An IVIM sequence with 16 b-values from 0 to 900 s/mm² was acquired, and the IVIM bi-exponential signal equation model was fitted in two steps. Perfusion asymmetries were evaluated using the blood flow related IVIM fD* parameter in regions of interest placed within the paraspinal muscles. Statistical significance was assessed using a Student t-test. The observed perfusion pattern after performing a Roman chair muscle exercise differed consistently in patients with AIS compared with healthy normal volunteers, and consisted of an asymmetrical increase in IVIM fD* [10^-3 mm² /s] above the lumbar convexity from 6.5 ± 5.8 to 28.8 ± 26.8 (p < 0.005), with no increase in the concavity (decrease from 6.5 ± 10.0 to 3.2 ± 1.5 (p = 0.19)), compared with a bilateral symmetric increase in the healthy volunteers (right, increase from 3.3 ± 2.1 to 10.1 ± 4.6 (p < 0.05); left, 6.7 ± 10.7 to 13.3 ± 7.0 (p < 0.05)). In conclusion, patients with AIS exhibit significant asymmetric muscle perfusion over the convexity of the scoliotic curvature after Roman chair exercise.