Intravoxel incoherent motion magnetic resonance imaging: basic principles and clinical applications

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.

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.

Multiparametric MRI of Solid Renal Masses: Principles and Applications of Advanced Quantitative and Functional Methods for Tumor Diagnosis and Characterization

Solid renal masses (SRMs) are increasingly detected and encompass both benign and malignant masses, with renal cell carcinoma (RCC) being the most common malignant SRM. Most patients with SRMs will undergo management without a priori pathologic confirmation. There is an unmet need to noninvasively diagnose and characterize RCCs, as significant variability in clinical behavior is observed and a wide range of differing management options exist. Cross-sectional imaging modalities, including magnetic resonance imaging (MRI), are increasingly used for SRM characterization. Multiparametric (mp) MRI techniques can provide insight into tumor biology by probing different physiologic/pathophysiologic processes noninvasively. These include sequences that probe tissue microstructure, including intravoxel incoherent motion diffusion-weighted imaging (IVIM-DWI) and T₁ relaxometry; oxygen metabolism (blood oxygen level dependent [BOLD-MRI]); as well as vascular flow and perfusion (dynamic contrast-enhanced MRI [DCE-MRI] and arterial spin labeling [ASL]). In this review, we will discuss each mpMRI method in terms of its principles, roles, and discuss the results of human studies for SRM assessment. Future validation of these methods may help to enable a personalized management approach for patients with SRM in the emerging era of precision medicine. EVIDENCE LEVEL: 5. TECHNICAL EFFICACY: 2.

Application of advanced diffusion models from diffusion weighted imaging in a large cohort study of breast lesions

BACKGROUND

To evaluate multiple parameters in multiple b-value diffusion-weighted imaging (DWI) in characterizing breast lesions and predicting prognostic factors and molecular subtypes.

METHODS

In total, 504 patients who underwent 3-T magnetic resonance imaging (MRI) with T1-weighted dynamic contrast-enhanced (DCE) sequences, T2-weighted sequences and multiple b-value (7 values, from 0 to 3000 s/mm²) DWI were recruited. The average values of 13 parameters in 6 models were calculated and recorded. The pathological diagnosis of breast lesions was based on the latest World Health Organization (WHO) classification.

RESULTS

Twelve parameters exhibited statistical significance in differentiating benign and malignant lesions. alpha demonstrated the highest sensitivity (89.5%), while sigma demonstrated the highest specificity (77.7%). The stretched-exponential model (SEM) demonstrated the highest sensitivity (90.8%), while the biexponential model demonstrated the highest specificity (80.8%). The highest AUC (0.882, 95% CI, 0.852-0.912) was achieved when all 13 parameters were combined. Prognostic factors were correlated with different parameters, but the correlation was relatively weak. Among the 6 parameters with significant differences among molecular subtypes of breast cancer, the Luminal A group and Luminal B (HER2 negative) group had relatively low values, and the HER2-enriched group and TNBC group had relatively high values.

CONCLUSIONS

All 13 parameters, independent or combined, provide valuable information in distinguishing malignant from benign breast lesions. These new parameters have limited meaning for predicting prognostic factors and molecular subtypes of malignant breast tumors.

Differential value of diffusion kurtosis imaging and intravoxel incoherent motion in benign and malignant solitary pulmonary lesions

Objective

To investigate the diagnostic value of diffusion kurtosis imaging (DKI) and intravoxel incoherent motion (IVIM) whole-lesion histogram parameters in differentiating benign and malignant solitary pulmonary lesions (SPLs).

Materials and Methods

Patients with SPLs detected by chest CT examination and with further routine MRI, DKI and IVIM-DWI functional sequence scanning data were recruited. According to the pathological results, SPLs were divided into a benign group and a malignant group. Independent samples t tests (normal distribution) or Mann‒Whitney U tests (nonnormal distribution) were used to compare the differences in DKI (Dk, K), IVIM (D, D*, f) and ADC whole-lesion histogram parameters between the benign and malignant SPL groups. The receiver operating characteristic (ROC) curve was used to evaluate the diagnostic efficiency of the histogram parameters and determine the optimal threshold. The area under the curve (AUC) of each histogram parameter was compared by the DeLong method. Spearman rank correlation was used to analyze the correlation between histogram parameters and malignant SPLs.

Results

Most of the histogram parameters for diffusion-related values (Dk, D, ADC) of malignant SPLs were significantly lower than those of benign SPLs, while most of the histogram parameters for the K value of malignant SPLs were significantly higher than those of benign SPLs. DKI (Dk, K), IVIM (D) and ADC were effective in differentiating benign and malignant SPLs and combined with multiple parameters of the whole-lesion histogram for the D value, had the highest diagnostic efficiency, with an AUC of 0.967, a sensitivity of 90.00% and a specificity of 94.03%. Most of the histogram parameters for the Dk, D and ADC values were negatively correlated with malignant SPLs, while most of the histogram parameters for the K value were positively correlated with malignant SPLs.

Conclusions

DKI (Dk, K) and IVIM (D) whole-lesion histogram parameters can noninvasively distinguish benign and malignant SPLs, and the diagnostic performance is better than that of DWI. Moreover, they can provide additional information on SPL microstructure, which has important significance for guiding clinical individualized precision diagnosis and treatment and has potential clinical application value.

Perspective: Disentangling the effects of tES on neurovascular unit

Transcranial electrical stimulation (tES) can modulate the neurovascular unit, including the perivascular space morphology, but the mechanisms are unclear. In this perspective article, we used an open-source "rsHRF toolbox" and an open-source functional magnetic resonance imaging (fMRI) transcranial direct current stimulation (tDCS) data set to show the effects of tDCS on the temporal profile of the haemodynamic response function (HRF). We investigated the effects of tDCS in the gray matter and at three regions of interest in the gray matter, namely, the anodal electrode (FC5), cathodal electrode (FP2), and an independent site remote from the electrodes (PZ). A "canonical HRF" with time and dispersion derivatives and a finite impulse response (FIR) model with three parameters captured the effects of anodal tDCS on the temporal profile of the HRF. The FIR model showed tDCS onset effects on the temporal profile of HRF for verum and sham tDCS conditions that were different from the no tDCS condition, which questions the validity of the sham tDCS (placebo). Here, we postulated that the effects of tDCS onset on the temporal profile of HRF are subserved by the effects on neurovascular coupling. We provide our perspective based on previous work on tES effects on the neurovascular unit, including mechanistic grey-box modeling of the effects of tES on the vasculature that can facilitate model predictive control (MPC). Future studies need to investigate grey-box modeling of online effects of tES on the neurovascular unit, including perivascular space, neurometabolic coupling, and neurovascular coupling, that can facilitate MPC of the tES dose-response to address the momentary ("state") and phenotypic ("trait") factors.

Advanced Techniques in Head and Neck Cancer Imaging: Guide to Precision Cancer Management

Precision treatment requires precision imaging. With the advent of various advanced techniques in head and neck cancer treatment, imaging has become an integral part of the multidisciplinary approach to head and neck cancer care from diagnosis to staging and also plays a vital role in response evaluation in various tumors. Conventional anatomic imaging (CT scan, MRI, ultrasound) remains basic and focuses on defining the anatomical extent of the disease and its spread. Accurate assessment of the biological behavior of tumors, including tumor cellularity, growth, and response evaluation, is evolving with recent advances in molecular, functional, and hybrid/multiplex imaging. Integration of these various advanced diagnostic imaging and nonimaging methods aids understanding of cancer pathophysiology and provides a more comprehensive evaluation in this era of precision treatment. Here we discuss the current status of various advanced imaging techniques and their applications in head and neck cancer imaging.

Heterogeneity evaluation of multi-high b-value apparent diffusion coefficient on cerebral ischemia in MCAO rat

Purpose

To assess brain damage in a rat model of cerebral ischemia based on apparent diffusion coefficient (ADC) data obtained from multi-high b-values and evaluate the relationship between Aquaporin 4 (AQP4) expression and ADC.

Methods

Thirty eight male Sprague-Dawley rats were randomized into two groups: (1) sham controls (n = 6) and (2) cerebral ischemia (successful model, n = 19). All rats underwent diffusion-weighted imaging (DWI) with both standard b-values and multi-high b-values (2,500-4,500 s/mm²) using a 3.0-T device. Standard ADC (ADC_st) maps and multi-high b-value ADCs (ADC_mh) were calculated, respectively. Aquaporin 4 expression was quantified using Western blot. Relative values of ADC_st and ADC_mh, AQP4 expression were compared between the sham group and the ischemia group. Correlations between ADC values and AQP4 expression were evaluated.

Results

At 0.5 h after suture insertion, the value of ADC_mh on the lesion was obviously decreased, and there was no difference in lesion volume when compared with ADC_st. After reperfusion, besides similar regions where ADC_st values decreased, we also found additional large values on ADC_mh within the cortex of the ipsilateral side or surrounding the lesion. The lesion evolution of the large value on ADC_mh was quite different from other indicators. But the total ADC_mh values were still significantly associated with ADC_st. The AQP4 protein expression level was appreciably increased after middle cerebral artery occlusion (MCAO), but there was no correlation between AQP4 expression either with ADC_mh or ADC_st.

Conclusion

We found the large values on ADC_mh during the progression of cerebral infarction is varied, but there was no correlation between ADC_mh values and AQP4 expression. ADC_mh may indicate the heterogeneity of ischemia lesions, but the underlying pathological basis should be further explored.

A comparative study of functional MRI in predicting response of regional nodes to induction chemotherapy in patients with nasopharyngeal carcinoma

Purpose

To identify and compare the value of functional MRI (fMRI) in predicting the early response of metastatic cervical lymph nodes (LNs) to induction chemotherapy (IC) in nasopharyngeal carcinoma (NPC) patients.

Methods

This prospective study collected 94 metastatic LNs from 40 consecutive NPC patients treated with IC from January 2021 to May 2021. Conventional diffusion-weighted imaging, diffusion kurtosis imaging, intravoxel incoherent motion, and dynamic contrast-enhanced magnetic resonance imaging were performed before and after IC. The parameter maps apparent diffusion coefficient (ADC), mean diffusion coefficient (MD), mean kurtosis (MK), Dslow, Dfast, perfusion fraction (PF), Ktrans, Ve, and Kep) of the metastatic nodes were calculated by the Functool postprocessing software. All LNs were classified as the responding group (RG) and non-responding group (NRG) according to Response Evaluation Criteria in Solid Tumors 1.1. The fMRI parameters were compared before and after IC and between the RG and the NRG. The significant parameters are fitted by logistic regression analysis to produce new predictive factor (PRE)–predicted probabilities. Logistic regression analysis and receiver operating characteristic (ROC) curves were performed to further identify and compare the efficacy of the parameters.

Results

After IC, the mean values of ADC, MD, and Dslow significantly increased, while MK, Dfast, and Ktrans values decreased dramatically, while no significant difference was detected in Ve and Kep. Compared with NRG, PF-pre and Ktrans-pre values in the RG were higher statistically. The areas under the ROC for the pretreatment PF, Ktrans, and PRE were 0.736, 0.722, and 0.810, respectively, with the optimal cutoff value of 222 × 10-4, 934 × 10-3/min, and 0.6624, respectively.

Conclusions

The pretreatment fMRI parameters PF and Ktrans showed promising potential in predicting the response of the metastatic LNs to IC in NPC patients.

Clinical Trial Registration

This study was approved by the ethics board of the Chinese PLA General Hospital, and registered on 30 January 2021, in the Chinese Clinical Trial Registry; http://www.chictr.org.cn/showproj.aspx?proj=121198, identifier (ChiCTR2100042863).

Comparison of the pre-treatment functional MRI metrics' efficacy in predicting Locoregionally advanced nasopharyngeal carcinoma response to induction chemotherapy

Background

Functional MRI (fMRI) parameters analysis has been proven to be a promising tool of predicting therapeutic response to induction chemotherapy (IC) in nasopharyngeal carcinoma (NPC). The study was designed to identify and compare the value of fMRI parameters in predicting early response to IC in patients with NPC.

Methods

This prospective study enrolled fifty-six consecutively NPC patients treated with IC from January 2021 to May 2021. Conventional diffusion weighted imaging (DWI), diffusion kurtosis imaging (DKI), intravoxel incoherent motion (IVIM) and dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) protocols were performed before and after IC. Parameters maps (ADC, MD, MK, D_slow, D_fast, PF, K^trans, V_e and K_ep) of the primary tumor were calculated by the Functool post-processing software. The participants were classified as responding group (RG) and non-responding group (NRG) according to Response Evaluation Criteria in Solid Tumors 1.1. The fMRI parameters were compared before and after IC and between RG with NRG. Logistic regression analysis and ROC were performed to further identify and compare the efficacy of the parameters.

Results

After IC, the mean values of ADC(p < 0.001), MD(p < 0.001), D_slow(p = 0.001), PF(p = 0.030) and V_e(p = 0.003) significantly increased, while MK(p < 0.001), D_fast(p = 0.009) and K_ep(p = 0.003) values decreased dramatically, while no significant difference was detected in K^trans(p = 0.130). Compared with NRG, ADC-pre(p < 0.001), MD-pre(p < 0.001) and D_slow-pre(p = 0.002) values in RG were lower, while MK-pre(p = 0.017) values were higher. The areas under the ROC curves for the ADC-pre, MD-pre, MK-pre, D_slow-pre and PRE were 0.885, 0.855, 0.809, 0.742 and 0.912, with the optimal cutoff value of 1210 × 10^− 6 mm²/s, 1010 × 10^− 6 mm²/s, 832 × 10^− 6, 835 × 10^− 6 mm²/s and 0.799 respectively.

Conclusions

The pretreatment conventional DWI (ADC), DKI (MD and MK), and IVIM (D_slow) values derived from fMRI showed a promising potential in predicting the response of the primary tumor to IC in NPC patients.

Trial registration

This study was approved by ethics board of the Chinese PLA General Hospital, and registered on January 30, 2021, in Chinese Clinical Trial Registry (ChiCTR2100042863).

Supplementary Information

The online version contains supplementary material available at 10.1186/s40644-021-00428-0.

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.