Diffusional kurtosis imaging (DKI) incorporation into an intravoxel incoherent motion (IVIM) MR model to measure cerebral hypoperfusion induced by hyperventilation challenge in healthy subjects

High-fidelity intravoxel incoherent motion parameter mapping using locally low-rank and subspace modeling

PURPOSE

Intravoxel incoherent motion (IVIM) is a quantitative magnetic resonance imaging (MRI) method used to quantify perfusion properties of tissue non-invasively without contrast. However, clinical applications are limited by unreliable parameter estimates, particularly for the perfusion fraction (f) and pseudodiffusion coefficient (D*). This study aims to develop a high-fidelity reconstruction for reliable estimation of IVIM parameters. The proposed method is versatile and amenable to various acquisition schemes and fitting methods.

METHODS

To address current challenges with IVIM, we adapted several advanced reconstruction techniques. We used a low-rank approximation of IVIM images and temporal subspace modeling to constrain the magnetization dynamics of the bi-exponential diffusion signal decay. In addition, motion-induced phase variations were corrected between diffusion directions and b-values, facilitating the use of high SNR real-valued diffusion data. The proposed method was evaluated in simulations and in vivo brain acquisitions in six healthy subjects and six individuals with a history of SARS-CoV-2 infection and compared with the conventionally reconstructed magnitude data. Following reconstruction, IVIM parameters were estimated voxel-wise.

RESULTS

Our proposed method reduced noise contamination in simulations, resulting in a 60%, 58.9%, and 83.9% reduction in the NRMSE for D, f, and D*, respectively, compared to the conventional reconstruction. In vivo, anisotropic properties of D, f, and D* were preserved with the proposed method, highlighting microvascular differences in gray matter between individuals with a history of COVID-19 and those without (p = 0.0210), which wasn't observed with the conventional reconstruction.

CONCLUSION

The proposed method yielded a more reliable estimation of IVIM parameters with less noise than the conventional reconstruction. Further, the proposed method preserved anisotropic properties of IVIM parameter estimates and demonstrated differences in microvascular perfusion in COVID-affected subjects, which weren't observed with conventional reconstruction methods.

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.

Magnetic Resonance Image Compilation Was Used in Conjunction with Prostate PI-RADS v2.1 Score Has Diagnostic Relevance for Benign and Malignant Prostate Lesions

Objective

To assess the diagnostic usefulness of magic in conjunction with PI-RADS v2.1 for prostate cancer malignant foci.

Methods

A total of 202 lesions (97 transitional zone lesions and 105 peripheral zone lesions) from 198 people were investigated retrospectively using traditional MRI and magic images. Each lesion has a unique pathological consequence. Lesions T1, T2, and PD values were employed as magic observation markers. The locations of the lesions were aggregated, and the paired t-test and receiver operating characteristic curve (ROC) were employed to find the indices with statistical significance in separating benign from malignant prostatic nodules (+1 point) and (−1 point) respectively. Draw a ROC curve and compare it to the PI-RADS v2.1 score using the magic positive and negative indices as well as the PI-RADS v2.1 score. By comparing the ROC curves scored separately, the diagnostic efficiency of the two scoring approaches for benign and malignant prostate lesions was investigated.

Results

T2 value has the highest diagnostic efficiency among the magic observation indices. T2 value of 77 ms for transitional zone lesions and T2 value of 89 ms for peripheral zone lesions are positive indices, whereas T2 value >77 ms and T2 value >89 ms are negative indexes. PI-RADS v2.1 combines one score and magic. In the transitional zone, the sensitivity, specificity, positive predictive value, and negative predictive value of the two scoring methods were 57.52, 87.70, 76.70, and 74.6 percent and 82.50, 73.68, 95.5, and 74.7 percent, respectively, and the AUC values were 0.735 and 0.846, respectively (P = 0.004); in the peripheral zone, the AUC values were 86.15 percent, 68.42 percent, 82.4.

Conclusions

Magic T2 value is a favorable sign for diagnosing benign and malignant prostate cancers when used in conjunction with PI-RADS v2.1. The end product exceeds PI-RADS v2.1 on its own, which is more useful in identifying benign and malignant prostate lesions, decreasing unnecessary puncture and alleviating patient pain.

The Histogram Analysis of Intravoxel Incoherent Motion-Kurtosis Model in the Diagnosis and Grading of Prostate Cancer-A Preliminary Study

Objectives

This study was conducted in order to explore the value of histogram analysis of the intravoxel incoherent motion-kurtosis (IVIM-kurtosis) model in the diagnosis and grading of prostate cancer (PCa), compared with monoexponential model (MEM).

Materials and Methods

Thirty patients were included in this study. Single-shot echo-planar imaging (SS-EPI) diffusion-weighted images (b-values of 0, 20, 50, 100, 200, 500, 1,000, 1,500, 2,000 s/mm²) were acquired. The pathologies were confirmed by in-bore MR-guided biopsy. The postprocessing and measurements were processed using the software tool Matlab R2015b for the IVIM-kurtosis model and MEM. Regions of interest (ROIs) were drawn manually. Mean values of D, D*, f, K, ADC, and their histogram parameters were acquired. The values of these parameters in PCa and benign prostatic hyperplasia (BPH)/prostatitis were compared. Receiver operating characteristic (ROC) curves were used to investigate the diagnostic efficiency. The Spearman test was used to evaluate the correlation of these parameters and Gleason scores (GS) of PCa.

Results

For the IVIM-kurtosis model, D (mean, 10th, 25th, 50th, 75th, 90th), D* (90th), and f (10th) were significantly lower in PCa than in BPH/prostatitis, while D (skewness), D* (kurtosis), and K (mean, 75th, 90th) were significantly higher in PCa than in BPH/prostatitis. For MEM, ADC (mean, 10th, 25th, 50th, 75th, 90th) was significantly lower in PCa than in BPH/prostatitis. The area under the ROC curve (AUC) of the IVIM-kurtosis model was higher than MEM, without significant differences (z = 1.761, P = 0.0783). D (mean, 50th, 75th, 90th), D* (mean, 10th, 25th, 50th, 75th), and f (skewness, kurtosis) correlated negatively with GS, while D (kurtosis), D* (skewness, kurtosis), f (mean, 75th, 90th), and K (mean, 75th, 90th) correlated positively with GS. The histogram parameters of ADC did not show correlations with GS.

Conclusion

The IVIM-kurtosis model has potential value in the differential diagnosis of PCa and BPH/prostatitis. IVIM-kurtosis histogram analysis may provide more information in the grading of PCa than MEM.

Crossed cerebellar diaschisis after acute ischemic stroke detected by intravoxel incoherent motion magnetic resonance imaging

PURPOSE

To study the value of 3.0 T magnetic resonance imaging with intravoxel incoherent motion (IVIM) in the diagnosis of the crossed cerebellar diaschisis (CCD) after the unilateral supratentorial acute ischemic stroke.

METHODS

Seventy-four patients with acute ischemic stroke who underwent intravoxel incoherent motion (IVIM), arterial spin labeling (ASL), and conventional magnetic resonance imaging (MRI) scanning were enrolled. Intravoxel incoherent motion-derived perfusion-related parameters including fast diffusion coefficient (D*), slow diffusion coefficient (D), vascular volume fraction (f), and arterial spin-labeling-derived cerebral blood flow (CBF) of bilateral cerebellum were measured.

RESULTS

In the CCD-positive group, D*, D, and CBF values of the contralateral cerebellum decreased compared with those of the ipsilesional cerebellum (P < 0.05), whereas f significantly increased (P < 0.05). A positive correlation was detected between the slow diffusion coefficient-based asymmetry index (AI-D) and the cerebral blood flow-based asymmetry index (AI-CBF) (r = 0.515, P < 0.01), whereas the vascular volume fraction-based asymmetry index (AI-f) had a negative correlation with the cerebral blood flow-based asymmetry index (AI-CBF) (r =  - 0.485, P < 0.01). Furthermore, the area under the receiver operating characteristic (ROC) curve value of AI-D and AI-f was 0.81 and 0.76, respectively.

CONCLUSIONS

The IVIM is feasible for the detection of CCD. This technique might provide opportunities to further investigate the pathophysiology of CCD.

Diffusion kurtosis imaging: correlation analysis of quantitative model parameters with molecular features in advanced lung adenocarcinoma

Background:

Due to development of magnetic resonance-based functional imaging, it is easier to detect micro-structural alterations of tumor tissues. The aim of this study was to conduct a preliminary evaluation of the correlation of non-Gaussian diffusion kurtosis imaging (DKI) parameters with expression of molecular markers (epidermal growth factor receptor [EGFR]; anaplastic lymphoma kinase [ALK]; Ki-67 protein) in patients with advanced lung adenocarcinoma, using routine diffusion-weighted imaging as the reference standard.

Methods:

Data from patients with primary lung adenocarcinoma diagnosed at Cancer Hospital, Chinese Academy of Medical Sciences (CHCAMS) from 2016 to 2019 were collected for retrospective analysis. The pathologic and magnetic resonance imaging data of 96 patients who met the inclusion criteria were included in this study. Specifically, the K_app and D_app parameters measured from the DKI model; apparent diffusion coefficient (ADC) value from the diffusion-weighted imaging model; and the EGFR, ALK, and Ki-67 biomarkers detected by immunohistochemistry and/or molecular biology techniques after biopsy or surgery were evaluated. The relations between quantitative parameters (ADC, K_app, D_app) and pathologic outcomes (EGFR, ALK, and Ki-67 expression) were analyzed by Spearman correlation test.

Results:

Of the 96 lung adenocarcinoma lesions (from 96 patients), the number of EGFR- and ALK-positive and high Ki-67 expressing lesions were 53, 12, and 83, respectively. The K_app values were significantly higher among patients with EGFR-positive mutations (0.81 ± 0.12 vs. 0.66 ± 0.10, t = 6.41, P < 0.001), ALK rearrangement-negative (0.76 ± 0.12 vs. 0.60 ± 0.15, t = 4.09, P < 0.001), and high Ki-67 proliferative index (PI) (0.76 ± 0.12 vs. 0.58 ± 0.13, t = 4.88, P < 0.001). The D_app values were significantly lower among patients with high Ki-67 PI (3.19 ± 0.69 μm²/ms vs. 4.20 ± 0.83 μm²/ms, t = 4.80, P < 0.001) and EGFR-positive mutations (3.11 ± 0.73 μm²/ms vs. 3.59 ± 0.77 μm²/ms, t = 3.12, P = 0.002). The differences in mean D_app (3.73 ± 1.26 μm²/ms vs. 3.26 ± 0.68 μm²/ms, t = 1.96, P = 0.053) or ADC values ([1.34 ± 0.81] × 10⁻³ mm²/s vs. [1.33 ± 0.41] × 10⁻³ mm²/s, t = 0.07, P = 0.941) between the groups with or without ALK rearrangements were not statistically significant. The ADC values were significantly lower among patients with EGFR-positive mutation ([1.19 ± 0.37] × 10⁻³ mm²/s vs. [1.50 ± 0.53] × 10⁻³ mm²/s, t = 3.38, P = 0.001) and high Ki-67 PI ([1.28 ± 0.39] × 10⁻³ mm²/s vs. [1.67 ± 0.77] × 10⁻³ mm²/s, t = 2.88, P = 0.005). K_app was strongly positively correlated with EGFR mutations (r = 0.844, P = 0.008), strongly positively correlated with Ki-67 PI (r = 0.882, P = 0.001), and strongly negatively correlated with ALK rearrangements (r = −0.772, P = 0.001). D_app was moderately correlated with EGFR mutations (r = −0.650, P = 0.024) or Ki-67 PI (r = −0.734, P = 0.012). ADC was moderately correlated with Ki-67 PI (r = −0.679, P = 0.033).

Conclusions:

The K_app value of DKI parameters was strongly correlated with different expression of EGFR, ALK, and Ki-67 in advanced lung adenocarcinoma. The results potentially indicate a surrogate measure of the status of different molecular markers assessed by non-invasive imaging tools.

Application of diffusion kurtosis tensor MR imaging in characterization of renal cell carcinomas with different pathological types and grades

Background

To probe the feasibility and reproducibility of diffusion kurtosis tensor imaging (DKTI) in renal cell carcinoma (RCC) and to apply DKTI in distinguishing the subtypes of RCC and the grades of clear cell RCC (CCRCC).

Methods

Thirty-eight patients with pathologically confirmed RCCs [CCRCC for 30 tumors, papillary RCC (PRCC) for 5 tumors and chromophobic RCC (CRCC) for 3 tumors] were involved in the study. Diffusion kurtosis tensor MR imaging were performed with 3 b-values (0, 500, 1000s/mm²) and 30 diffusion directions. The mean kurtosis (MK), axial kurtosis (Ka), radial kurtosis (Kr) values and mean diffusity (MD) for RCC and contralateral normal parenchyma were acquired. The inter-observer agreements of all DKTI metrics of contralateral renal cortex and medulla were evaluated using Bland-Altman plots. Statistical comparisons with DKTI metrics of 3 RCC subtypes and between low-grade (Furman grade I ~ II, 22 cases) and high-grade (Furman grade III ~ IV, 8 cases) CCRCC were performed with ANOVA test and Student t test separately. Receiver operating characteristic (ROC) curve analyses were used to compare the diagnostic efficacy of DKTI metrics for predicting nuclear grades of CCRCC. Correlations between DKTI metrics and nuclear grades were also evaluated with Spearman correlation analysis.

Results

Inter-observer measurements for each metric showed great reproducibility with excellent ICCs ranging from 0.81 to 0.87. There were significant differences between the DKTI metrics of RCCs and contralateral renal parenchyma, also among the subtypes of RCC. MK and Ka values of CRCC were significantly higher than those of CCRCC and PRCC. Statistical difference of the MK, Ka, Kr and MD values were also obtained between CCRCC with high- and low-grades. MK values were more effective for distinguishing between low- and high- grade CCRCC (area under the ROC curve: 0.949). A threshold value of 0.851 permitted distinction with high sensitivity (90.9%) and specificity (87.5%).

Conclusion

Our preliminary results suggest a possible role of DKTI in differentiating CRCC from CCRCC and PRCC. MK, the principle DKTI metric might be a surrogate biomarker to predict nuclear grades of CCRCC.

Trial registration

ChiCTC, ChiCTR-DOD-17010833, Registered 10 March, 2017, retrospectively registered, http://www.chictr.org.cn/showproj.aspx?proj=17559.

Optimal Model Mapping for Intravoxel Incoherent Motion MRI

In general, only one diffusion model would be applied to whole field-of-view voxels in the intravoxel incoherent motion-magnetic resonance imaging (IVIM-MRI) study. However, the choice of the applied diffusion model can significantly influence the estimated diffusion parameters. The quality of the diffusion analysis can influence the reliability of the perfusion analysis. This study proposed an optimal model mapping method to improve the reliability of the perfusion parameter estimation in the IVIM study. Six healthy volunteers (five males and one female; average age of 38.3 ± 7.5 years). Volunteers were examined using a 3.0 Tesla scanner. IVIM-MRI of the brain was applied at 17 b-values ranging from 0 to 2,500 s/mm². The Gaussian model, the Kurtosis model, and the Gamma model were found to be optimal for the CSF, white matter (WM), and gray matter (GM), respectively. In the mean perfusion fraction (f_p) analysis, the GM/WM ratios were 1.16 (Gaussian model), 1.80 (Kurtosis model), 1.94 (Gamma model), and 1.54 (Optimal model mapping); in the mean pseudo diffusion coefficient (D^*) analysis, the GM/WM ratios were 1.18 (Gaussian model), 1.19 (Kurtosis model), 1.56 (Gamma model), and 1.24 (Optimal model mapping). With the optimal model mapping method, the estimated f_p and D^* were reliable compared with the conventional methods. In addition, the optimal model maps, the associated products of this method, may provide additional information for clinical diagnosis.

Isotropically weighted intravoxel incoherent motion brain imaging at 7T

Perfusion magnetic resonance imaging (MRI) is a promising non-invasive technique providing insights regarding the brain's microvascular architecture in vivo. The scalar perfusion metrics can be used for quantitative diagnostics of various brain abnormalities, in particular, in the stroke cases and tumours. However, conventional MRI-based perfusion approaches such as dynamic contrast-enhanced perfusion imaging or arterial spin labelling have a few weaknesses, for instance, contrast agent deposition, low signal-to-noise ratio, limited temporal and spatial resolution, and specific absorption rate constraints. As an alternative, the intravoxel incoherent motion (IVIM) approach exploits an extension of diffusion MRI in order to estimate perfusion parameters in the human brain. Application of IVIM imaging at ultra-high field MRI might employ the advantage of a higher signal-to-noise ratio, and thereby the use of higher spatial and temporal resolutions. In the present work, we demonstrate an application of recently developed isotropic diffusion weighted sequences to the evaluation of IVIM parameters at an ultra-high 7T field. The used sequence exhibits high immunity to image degrading factors and allows one to acquire the data in a fast and efficient way. Utilising the bi-exponential fitting model of the signal attenuation, we performed an extensive analysis of the IVIM scalar metrics obtained by a isotropic diffusion weighted sequence in vivo and compared results with a conventional pulsed gradient sequence at 7T. In order to evaluate a possible metric bias originating from blood flows, we additionally used a truncated b-value protocol (b-values from 100 to 200 s/mm² with the step 20 s/mm²) accompanied to the full range (b-values from 0 to 200 s/mm²). The IVIM scalar metrics have been assessed and analysed together with a large and middle vessel density atlas of the human brain. We found that the diffusion coefficients and perfusion fractions of the voxels consisting of large and middle vessels have higher values in contrast to other tissues. Additionally, we did not find a strong dependence of the IVIM metrics on the density values of the vessel atlas. Perspectives and limitations of the developed isotropic diffusion weighted perfusion are presented and discussed.

Characterization of active and infiltrative tumorous subregions from normal tissue in brain gliomas using multiparametric MRI

BACKGROUND

Targeted localized biopsies and treatments for diffuse gliomas rely on accurate identification of tissue subregions, for which current MRI techniques lack specificity.

PURPOSE

To explore the complementary and competitive roles of a variety of conventional and quantitative MRI methods for distinguishing subregions of brain gliomas.

STUDY TYPE

Prospective.

POPULATION

Fifty-one tissue specimens were collected using image-guided localized biopsy surgery from 10 patients with newly diagnosed gliomas.

FIELD STRENGTH/SEQUENCE

Conventional and quantitative MR images consisting of pre- and postcontrast T1 w, T2 w, T2 -FLAIR, T2 -relaxometry, DWI, DTI, IVIM, and DSC-MRI were acquired preoperatively at 3T.

ASSESSMENT

Biopsy specimens were histopathologically attributed to glioma tissue subregion categories of active tumor (AT), infiltrative edema (IE), and normal tissue (NT) subregions. For each tissue sample, a feature vector comprising 15 MRI-based parameters was derived from preoperative images and assessed by a machine learning algorithm to determine the best multiparametric feature combination for characterizing the tissue subregions.

STATISTICAL TESTS

For discrimination of AT, IE, and NT subregions, a one-way analysis of variance (ANOVA) test and for pairwise tissue subregion differentiation, Tukey honest significant difference, and Games-Howell tests were applied (P < 0.05). Cross-validated feature selection and classification methods were implemented for identification of accurate multiparametric MRI parameter combination.

RESULTS

After exclusion of 17 tissue specimens, 34 samples (AT = 6, IE = 20, and NT = 8) were considered for analysis. Highest accuracies and statistically significant differences for discrimination of IE from NT and AT from NT were observed for diffusion-based parameters (AUCs >90%), and the perfusion-derived parameter as the most accurate feature in distinguishing IE from AT. A combination of "CBV, MD, T2 _ISO, FLAIR" parameters showed high diagnostic performance for identification of the three subregions (AUC ∼90%).

DATA CONCLUSION

Integration of a few quantitative along with conventional MRI parameters may provide a potential multiparametric imaging biomarker for predicting the histopathologically proven glioma tissue subregions.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 3 J. Magn. Reson. Imaging 2018;48:938-950.

Measurement of the perfusion fraction in brain tumors with intravoxel incoherent motion MR imaging: validation with histopathological vascular density in meningiomas

OBJECTIVE

To evaluate the quantification performance of the perfusion fraction (f) measured with intravoxel incoherent motion (IVIM) MR imaging in a comparison with the histological vascular density in meningiomas.

METHODS

29 consecutive patients with meningioma (59.0 ± 16.8 years old, 8 males and 21 females) who underwent a subsequent surgical resection were examined with both IVIM imaging and a histopathological analysis. IVIM imaging was conducted using a single-shot SE-EPI sequence with 13 b-factors (0, 10, 20, 30, 50, 80, 100, 200, 300, 400, 600, 800, 1000 s mm^-2) at 3T. The perfusion fraction (f) was calculated by fitting the IVIM bi-exponential model. The 90-percentile f-value in the tumor region-of-interest (ROI) was defined as the maximum f-value (f-max). Histopathological vascular density (%Vessel) was measured on CD31-immunostainted histopathological specimens. The correlation and agreement between the f-values and %Vessel was assessed.

RESULTS

The f-max (15.5 ± 5.5%) showed excellent agreement [intraclass correlation coefficient (ICC) = 0.754] and a significant correlation (r = 0.69, p < 0.0001) with the %Vessel (12.9 ± 9.4%) of the tumors. The Bland-Altman plot analysis showed excellent agreement between the f-max and %Vessel (bias, -2.6%; 95% limits of agreement, from -16.0 to 10.8%). The f-max was not significantly different among the histological subtypes of meningioma.

CONCLUSION

An excellent agreement and a significant correlation were observed between the f-values and %Vessel. The f-value can be used as a noninvasive quantitative imaging measure to directly assess the vascular volume fraction in brain tumors. Advances in knowledge: The f-value measured by IVIM imaging showed a significant correlation and an excellent agreement with the histological vascular density in the meningiomas. The f-value can be used as a noninvasive and quantitative imaging measure to directly assess the volume fraction of capillaries in brain tumors.