MR diffusion kurtosis imaging for cancer diagnosis: A meta-analysis of the diagnostic accuracy of quantitative kurtosis value and diffusion coefficient

Optimization of quasi-diffusion magnetic resonance imaging for quantitative accuracy and time-efficient acquisition

PURPOSE

Quasi-diffusion MRI (QDI) is a novel quantitative technique based on the continuous time random walk model of diffusion dynamics. QDI provides estimates of the diffusion coefficient, <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mspace/> <mml:msub><mml:mi>D</mml:mi> <mml:mrow><mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {D}_{1,2} $$</mml:annotation></mml:semantics> </mml:math> in mm²  s^-1 and a fractional exponent, <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> , defining the non-Gaussianity of the diffusion signal decay. Here, the b-value selection for rapid clinical acquisition of QDI tensor imaging (QDTI) data is optimized.

METHODS

Clinically appropriate QDTI acquisitions were optimized in healthy volunteers with respect to a multi-b-value reference (MbR) dataset comprising 29 diffusion-sensitized images arrayed between <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>b</mml:mi> <mml:mo>=</mml:mo> <mml:mn>0</mml:mn></mml:mrow> <mml:annotation>$$ b=0 $$</mml:annotation></mml:semantics> </mml:math> and 5000 s mm^-2 . The effects of varying maximum b-value ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> </mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}} $$</mml:annotation></mml:semantics> </mml:math> ), number of b-value shells, and the effects of Rician noise were investigated.

RESULTS

QDTI measures showed <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> </mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}} $$</mml:annotation></mml:semantics> </mml:math> dependence, most significantly for <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> in white matter, which monotonically decreased with higher <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> </mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}} $$</mml:annotation></mml:semantics> </mml:math> leading to improved tissue contrast. Optimized 2 b-value shell acquisitions showed small systematic differences in QDTI measures relative to MbR values, with overestimation of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mspace/> <mml:mspace/> <mml:msub><mml:mi>D</mml:mi> <mml:mrow><mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ \kern0.50em {D}_{1,2} $$</mml:annotation></mml:semantics> </mml:math> and underestimation of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> in white matter, and overestimation of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>D</mml:mi> <mml:mrow><mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {D}_{1,2} $$</mml:annotation></mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> anisotropies in gray and white matter. Additional shells improved the accuracy, precision, and reliability of QDTI estimates with 3 and 4 shells at <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> <mml:mo>=</mml:mo> <mml:mn>5000</mml:mn></mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}}=5000 $$</mml:annotation></mml:semantics> </mml:math>  s mm^-2 , and 4 b-value shells at <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> <mml:mo>=</mml:mo> <mml:mn>3960</mml:mn></mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}}=3960 $$</mml:annotation></mml:semantics> </mml:math>  s mm^-2 , providing minimal bias in <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>D</mml:mi> <mml:mrow><mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {D}_{1,2} $$</mml:annotation></mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> compared to the MbR.

CONCLUSION

A highly detailed optimization of non-Gaussian dMRI for in vivo brain imaging was performed. QDI provided robust parameterization of non-Gaussian diffusion signal decay in clinically feasible imaging times with high reliability, accuracy, and precision of QDTI measures.

Diagnostic Performance of Diffusion Kurtosis Imaging for Benign and Malignant Breast Lesions: A Systematic Review and Meta-Analysis

Purpose

Magnetic resonance imaging (MRI) has a high sensitivity for differentiating between malignant and non-malignant breast lesions but is sometimes limited due to its low specificity. Here, we performed a meta-analysis to evaluate the diagnostic performance of mean kurtosis (MK) and mean diffusivity (MD) values in magnetic resonance diffusion kurtosis imaging (DKI) for benign and malignant breast lesions.

Methods

Original articles on relevant topics, published from 2010 to 2019, in PubMed, EMBASE, and WanFang databases were systematically reviewed. According to the purpose of the study and the characteristics of DKI reported, the diagnostic performances of MK and MD were evaluated, and meta-regression was conducted to explore the source of heterogeneity.

Results

Fourteen studies involving 1,099 (451 benign and 648 malignant) lesions were analyzed. The pooled sensitivity, pooled specificity, positive likelihood ratio, and negative likelihood ratio for MD were 0.84 (95% confidence interval (CI), 0.81-0.87), 0.83 (95% CI, 0.79-0.86), 4.44 (95% CI, 3.54-5.57), and 0.18 (95% CI, 0.13-0.26), while those for MK were 0.89 (95% CI, 0.86-0.91), 0.86 (95% CI, 0.82-0.89), 5.72 (95% CI, 4.26-7.69), and 0.13 (95% CI, 0.09-0.19), respectively. The overall area under the curve (AUC) was 0.91 for MD and 0.95 for MK.

Conclusions

Analysis of the data from 14 studies showed that MK had a higher pooled sensitivity, pooled specificity, and diagnostic performance for differentiating between breast lesions, compared with MD.

Multi-parameter diffusion and perfusion magnetic resonance imaging and radiomics nomogram for preoperative evaluation of aquaporin-1 expression in rectal cancer

PURPOSE

The overexpression of aquaporin-1 (AQP1) is associated with poor prognosis in rectal cancer. This study aimed to explore the value of multi-parameter diffusion and perfusion MRI and radiomics models in predicting AQP1 high expression.

METHODS

This prospective study was performed from July 2019 to February 2021, which included rectal cancer participants after preoperative rectal MRI, with diffusion-weighted imaging, intravoxel incoherent motion (IVIM), diffusion kurtosis imaging (DKI), and dynamic contrast-enhanced (DCE) sequences. Radiomic features were extracted from MR images, and immunohistochemical tests assessed AQP1 expression. Selected quantitative MRI and radiomic features were analyzed. Receiver operating characteristic (ROC) curves evaluated the predictive performance. The nomogram performance was evaluated by its calibration, discrimen, and clinical utility. The intraclass correlation coefficient evaluated the interobserver agreement for the MRI features.

RESULTS

110 participants with the age of 60.7 ± 12.5 years been enrolled in this study. The apparent diffusion coefficient (ADC), IVIM_D, DKI_diffusivity, and DCE_Ktrans were significantly higher in participants with high AQP1 expression than in those with low expression (P < 0.05). ADC (b = 1000, 2000, and 3000 s/mm²), IVIM_D, DKI_diffusivity, and DCE_Ktrans were positively correlated (r = 0.205, 0.275, 0.37, 0.235, 0.229, and 0.227, respectively; P < 0.05), whereas DKI_Kurtosis was negatively correlated (r = - 0.22, P = 0.021) with AQP1 expression. ADC (b = 3000 s/mm²), IVIM_D, DKI_ diffusivity, DKI_Kurtosis, and DCE_Ktrans had moderate diagnostic efficiencies for high AQP1 expression (AUC = 0.715, 0.636, 0.627, 0.633, and 0.632, respectively; P < 0.05). The radiomic features had excellent predictive efficiency for high AQP1 expression (AUC = 0.967 and 0.917 for training and validation). The model-based nomogram had C-indexes of 0.932 and 0.851 for the training and validation cohorts, which indicated good fitting to the calibration curves (p > 0.05).

CONCLUSION

Diffusion and perfusion MRI can indicate the aquaporin-1 expression in rectal cancer, and radiomic features can enhance the predictive efficiency for high AQP1 expression. A nomogram for high aquaporin-1 expression will improve clinical decision-making.

Radiomics Analysis Based on Diffusion Kurtosis Imaging and T2 Weighted Imaging for Differentiation of Pancreatic Neuroendocrine Tumors From Solid Pseudopapillary Tumors

Objective

To develop and validate a radiomics model of diffusion kurtosis imaging (DKI) and T2 weighted imaging for discriminating pancreatic neuroendocrine tumors (PNETs) from solid pseudopapillary tumors (SPTs).

Materials and Methods

Sixty-six patients with histopathological confirmed PNETs (n = 31) and SPTs (n = 35) were enrolled in this study. ROIs of tumors were manually drawn on each slice at T2WI and DWI (b = 1,500 s/mm²) from 3T MRI. Intraclass correlation coefficients were used to evaluate the interobserver agreement. Mean diffusivity (MD) and mean kurtosis (MK) were derived from DKI. The least absolute shrinkage and selection operator regression were used for feature selection.

Results

MD and MK had a moderate diagnostic performancewith the area under curve (AUC) of 0.71 and 0.65, respectively. A radiomics model, which incorporated sex and age of patients and radiomics signature of the tumor, showed excellent discrimination performance with AUC of 0.97 and 0.86 in the primary and validation cohort. Moreover, the new model had better diagnostic performance than that of MD (P = 0.023) and MK (P = 0.004), and showed excellent differentiation with a sensitivity of 95.00% and specificity of 91.67% in primary cohort, and the sensitivity of 90.91% and specificity of 81.82% in the validation cohort. The accuracy of radiomics analysis, radiologist 1, and radiologist 2 for diagnosing SPTs and PNETs were 92.42, 77.27, and 78.79%, respectively. The accuracy of radiomics analysis was significantly higher than that of subjective diagnosis (P < 0.05).

Conclusions

Radiomics model could improve the diagnostic accuracy of SPTs and PNETs and contribute to determining an appropriate treatment strategy for pancreatic tumors.

Higher-order diffusion MRI characterization of mesorectal lymph nodes in rectal cancer

PURPOSE

Mesorectal lymph node staging plays an important role in treatment decision making. Here, we explore the benefit of higher-order diffusion MRI models accounting for non-Gaussian diffusion effects to classify mesorectal lymph nodes both 1) ex vivo at ultrahigh field correlated with histology and 2) in vivo in a clinical scanner upon patient staging.

METHODS

The preclinical investigation included 54 mesorectal lymph nodes, which were scanned at 16.4 T with an extensive diffusion MRI acquisition. Eight diffusion models were compared in terms of goodness of fit, lymph node classification ability, and histology correlation. In the clinical part of this study, 10 rectal cancer patients were scanned with diffusion MRI at 1.5 T, and 72 lymph nodes were analyzed with Apparent Diffusion Coefficient (ADC), Intravoxel Incoherent Motion (IVIM), Kurtosis, and IVIM-Kurtosis.

RESULTS

Compartment models including restricted and anisotropic diffusion improved the preclinical data fit, as well as the lymph node classification, compared to standard ADC. The comparison with histology revealed only moderate correlations, and the highest values were observed between diffusion anisotropy metrics and cell area fraction. In the clinical study, the diffusivity from IVIM-Kurtosis was the only metric showing significant differences between benign (0.80 ± 0.30 μm² /ms) and malignant (1.02 ± 0.41 μm² /ms, P = .03) nodes. IVIM-Kurtosis also yielded the largest area under the receiver operating characteristic curve (0.73) and significantly improved the node differentiation when added to the standard visual analysis by experts based on T₂ -weighted imaging.

CONCLUSION

Higher-order diffusion MRI models perform better than standard ADC and may be of added value for mesorectal lymph node classification in rectal cancer patients.

Diffusion kurtosis imaging and conventional diffusion weighted imaging to assess electrochemotherapy response in locally advanced pancreatic cancer

Background The aim of the study was to evaluate diagnostic performance of functional parameters derived by conventional mono-exponential approach of diffusion weighted imaging (DWI) and by diffusion kurtosis imaging (DKI) in the assessment of pancreatic tumours treated with electrochemotherapy (ECT). Patients and methods Twenty-one consecutive patients with locally advanced pancreatic adenocarcinoma subjected to ECT were enrolled in a clinical approved trial. Among twenty-one enrolled patients, 13/21 (61.9%) patients were subjected to MRI before and after ECT. DWI was performed with a 1.5 T scanner; a free breathing axial single shot echo planar DWI pulse sequence parameters were acquired using seven b value = 0, 50, 100, 150, 400, 800, 1000 s/mm2. Apparent diffusion coefficient by conventional mono-exponential approach and mean of diffusion coefficient (MD) and mean of diffusional kurtosis (MK) by DKI approach were derived from DWI. Receiver operating characteristic (ROC) analysis was performed and sensitivity, specificity, positive and negative predictive value were calculated. Results Among investigated diffusion parameters, only the MD derived by DKI showed a significant variation of values between pre and post treatment (p = 0.02 at Wilcoxon test) and a significant statistically difference for percentage change between responders and not responders (p = 0.01 at Kruskal Wallis test). MD had a good diagnostic performance with a sensitivity of 80%, a specificity of 100% and area under ROC of 0.933. Conclusions MD derived by DKI allows identifying responders and not responders patients subject to ECT treatment. MD had higher diagnostic performance to assess ECT response compared to conventional DWI derived parameters.