Differentiation of Low- and High-Grade Gliomas Using High b-Value Diffusion Imaging with a Non-Gaussian Diffusion Model

Value of fractional-order calculus (FROC) model diffusion-weighted imaging combined with simultaneous multi-slice (SMS) acceleration technology for evaluating benign and malignant breast lesions

BACKGROUND

This study explores the diagnostic value of combining fractional-order calculus (FROC) diffusion-weighted model with simultaneous multi-slice (SMS) acceleration technology in distinguishing benign and malignant breast lesions.

METHODS

178 lesions (73 benign, 105 malignant) underwent magnetic resonance imaging with diffusion-weighted imaging using multiple b-values (14 b-values, highest 3000 s/mm2). Independent samples t-test or Mann-Whitney U test compared image quality scores, FROC model parameters (D,, ), and ADC values between two groups. Multivariate logistic regression analysis identified independent variables and constructed nomograms. Model discrimination ability was assessed with receiver operating characteristic (ROC) curve and calibration chart. Spearman correlation analysis and Bland-Altman plot evaluated parameter correlation and consistency.

RESULTS

Malignant lesions exhibited lower D, and ADC values than benign lesions (P < 0.05), with higher values (P < 0.05). In SSEPI-DWI and SMS-SSEPI-DWI sequences, the AUC and diagnostic accuracy of D value are maximal, with D value demonstrating the highest diagnostic sensitivity, while value exhibits the highest specificity. The D and combined model had the highest AUC and accuracy. D and ADC values showed high correlation between sequences, and moderate. Bland-Altman plot demonstrated unbiased parameter values.

CONCLUSION

SMS-SSEPI-DWI FROC model provides good image quality and lesion characteristic values within an acceptable time. It shows consistent diagnostic performance compared to SSEPI-DWI, particularly in D and values, and significantly reduces scanning time.

Comparison of continuous-time random walk and fractional order calculus models in characterizing breast lesions using histogram analysis

OBJECTIVE

To compare the diagnostic performance of different mathematical models for DWI and explore whether parameters reflecting spatial and temporal heterogeneity can demonstrate better diagnostic accuracy than the diffusion coefficient parameter in distinguishing benign and malignant breast lesions, using whole-tumor histogram analysis.

METHODS

This retrospective study was approved by the institutional ethics committee and included 104 malignant and 42 benign cases. All patients underwent breast magnetic resonance imaging (MRI) with a 3.0 T MR scanner using the simultaneous multi-slice (SMS) readout-segment ed echo-planar imaging (rs-EPI). Histogram metrics of Mono- apparent diffusion coefficient (ADC), CTRW, and FROC-derived parameters were compared between benign and malignant breast lesions, and the diagnostic performance of each diffusion parameter was evaluated. Statistical analysis was performed using Mann-Whitney U test and receiver operating characteristic (ROC) curve.

RESULTS

The DFROC-median exhibited the highest AUC for distinguishing benign and malignant breast lesions (AUC = 0.965). The temporal heterogeneity parameter αCTRW-median generated a statistically higher AUC compared to the spatial heterogeneity parameter βCTRW-median (AUC = 0.850 and 0.741, respectively; p = 0.047). Finally, the combination of median values of CTRW parameters displayed a slightly higher AUC than that of FROC parameters, with no significant difference however (AUC = 0.971 and 0.965, respectively; p = 0.172).

CONCLUSIONS

The diffusion coefficient parameter exhibited superior diagnostic performance in distinguishing breast lesions when compared to the temporal and spatial heterogeneity parameters.

Multi-b-value DWI to evaluate the synergistic antiproliferation and anti-heterogeneity effects of bufalin plus sorafenib in an orthotopic HCC model

BACKGROUND

Multi-b-value diffusion-weighted imaging (DWI) with different postprocessing models allows for evaluating hepatocellular carcinoma (HCC) proliferation, spatial heterogeneity, and feasibility of treatment strategies. We assessed synergistic effects of bufalin+sorafenib in orthotopic HCC-LM3 xenograft nude mice by using intravoxel incoherent motion (IVIM), diffusion kurtosis imaging (DKI), a stretched exponential model (SEM), and a fractional-order calculus (FROC) model.

METHODS

Twenty-four orthotopic HCC-LM3 xenograft mice were divided into bufalin+sorafenib, bufalin, sorafenib treatment groups, and a control group. Multi-b-value DWI was performed using a 3-T scanner after 3 weeks' treatment to obtain true diffusion coefficient Dt, pseudo-diffusion coefficient Dp, perfusion fraction f, mean diffusivity (MD), mean kurtosis (MK), distributed diffusion coefficient (DDC), heterogeneity index α, diffusion coefficient D, fractional order parameter β, and microstructural quantity μ. Necrotic fraction (NF), standard deviation (SD) of hematoxylin-eosin staining, and microvessel density (MVD) of anti-CD31 staining were evaluated. Correlations of DWI parameters with histopathological results were analyzed, and measurements were compared among four groups.

RESULTS

In the final 22 mice, f positively correlated with MVD (r = 0.679, p = 0.001). Significantly good correlations of MK (r = 0.677), α (r = -0.696), and β (r= -0.639) with SD were observed (all p < 0.010). f, MK, MVD, and SD were much lower, while MD, α, β, and NF were higher in bufalin plus sorafenib group than control group (all p < 0.050).

CONCLUSION

Evaluated by IVIM, DKI, SEM, and FROC, bufalin+sorafenib was found to inhibit tumor proliferation and angiogenesis and reduce spatial heterogeneity in HCC-LM3 models.

RELEVANCE STATEMENT

Multi-b-value DWI provides potential metrics for evaluating the efficacy of treatment in HCC.

KEY POINTS

• Bufalin plus sorafenib combination may increase the effectiveness of HCC therapy. • Multi-b-value DWI depicted HCC proliferation, angiogenesis, and spatial heterogeneity. • Multi-b-value DWI may be a noninvasive method to assess HCC therapeutic efficacy.

Ultra high b-value diffusion weighted imaging enables better molecular grading stratification over histological grading in adult-type diffuse glioma

PURPOSE

Accurate preoperative radiological staging of adult-type diffuse glioma is crucial for effective prognostic stratification and selection of appropriate therapeutic interventions. The purpose of this study was to compare the effectiveness of apparent diffusion coefficient (ADC) maps generated from ultrahigh b-value diffusion-weighted imaging (DWI) for molecular grading with that for histological grading of adult-type diffuse glioma, and to evaluate the correlation between these ADC maps and molecular and histological biomarkers.

METHODS

This study retrospectively enrolled forty adult-type diffuse glioma patients, diagnosed using the 2021 WHO classification criteria. Preoperative imaging data, including multiple b-value DWI and conventional magnetic resonance imaging, were collected. Tumors were graded using both histological and molecular criteria. Histogram analysis was conducted to generate 14 parameters for each tumor. Receiver operating characteristic curves and the area under the curve (AUC) were used to evaluate tumor grading and molecular status differentiation. Analysis of histological biomarkers was performed by calculating the Pearson and Spearman correlation coefficients of continuous and hierarchical variables, respectively.

RESULTS

The intensity-related parameters for molecular grading were found to be superior to those for histological grading for the identification of WHO grade 4 (WHO4) adult-type diffuse glioma. The AUC of both grading systems increased with increasing b-values, with ADC8000-based histogram parameters showing the best results (molecular grading, square root: AUC = 0.897; histological grading, median: AUC = 0.737). The intensity-related parameters could also differentiate molecular WHO4 gliomas from histologically lower-grade gliomas (ADC8000-based square root: AUC = 0.919), and different ADC8000-based kurtosis was observed between molecular and histological WHO4 gliomas (AUC = 0.833). Significant correlations between the Ki-67 index and molecular status prediction for IDH, CDKN2A, and EGFR were also demonstrated.

CONCLUSION

The histogram parameters derived from high b-value ADC maps were found to be more effective for differentiating molecular grades of WHO4 adult-type diffuse glioma than for differentiating histological grades.

Accelerating High b-Value Diffusion-Weighted MRI Using a Convolutional Recurrent Neural Network (CRNN-DWI)

PURPOSE

To develop a novel convolutional recurrent neural network (CRNN-DWI) and apply it to reconstruct a highly undersampled (up to six-fold) multi-b-value, multi-direction diffusion-weighted imaging (DWI) dataset.

METHODS

A deep neural network that combines a convolutional neural network (CNN) and recurrent neural network (RNN) was first developed by using a set of diffusion images as input. The network was then used to reconstruct a DWI dataset consisting of 14 b-values, each with three diffusion directions. For comparison, the dataset was also reconstructed with zero-padding and 3D-CNN. The experiments were performed with undersampling rates (R) of 4 and 6. Standard image quality metrics (SSIM and PSNR) were employed to provide quantitative assessments of the reconstructed image quality. Additionally, an advanced non-Gaussian diffusion model was employed to fit the reconstructed images from the different approaches, thereby generating a set of diffusion parameter maps. These diffusion parameter maps from the different approaches were then compared using SSIM as a metric.

RESULTS

Both the reconstructed diffusion images and diffusion parameter maps from CRNN-DWI were better than those from zero-padding or 3D-CNN. Specifically, the average SSIM and PSNR of CRNN-DWI were 0.750 ± 0.016 and 28.32 ± 0.69 (R = 4), and 0.675 ± 0.023 and 24.16 ± 0.77 (R = 6), respectively, both of which were substantially higher than those of zero-padding or 3D-CNN reconstructions. The diffusion parameter maps from CRNN-DWI also yielded higher SSIM values for R = 4 (>0.8) and for R = 6 (>0.7) than the other two approaches (for R = 4, <0.7, and for R = 6, <0.65).

CONCLUSIONS

CRNN-DWI is a viable approach for reconstructing highly undersampled DWI data, providing opportunities to reduce the data acquisition burden.

Assessment of Prognostic Factors and Molecular Subtypes of Breast Cancer With a Continuous-Time Random-Walk MR Diffusion Model: Using Whole Tumor Histogram Analysis

BACKGROUND

The continuous-time random-walk (CTRW) diffusion model to evaluate breast cancer prognosis is rarely reported.

PURPOSE

To investigate the correlations between apparent diffusion coefficient (ADC) and CTRW-specific parameters with prognostic factors and molecular subtypes of breast cancer.

STUDY TYPE

Retrospective.

POPULATION

One hundred fifty-seven women (median age, 50 years; range, 26-81 years) with histopathology-confirmed breast cancer.

FIELD STRENGTH/SEQUENCE

Simultaneous multi-slice readout-segmented echo-planar imaging at 3.0T.

ASSESSMENT

The histogram metrics of ADC, anomalous diffusion coefficient (D), temporal diffusion heterogeneity (α), and spatial diffusion heterogeneity (β) were calculated for whole-tumor volume. Associations between histogram metrics and prognostic factors (estrogen receptor [ER], progesterone receptor [PR], human epidermal growth factor receptor 2 [HER2], and Ki-67 proliferation index), axillary lymph node metastasis (ALNM), and tumor grade were assessed. The performance of histogram metrics, both alone and in combination, for differentiating molecular subtypes (HER2-positive, Luminal or triple negative) was also assessed.

STATISTICAL TESTS

Comparisons were made using Mann-Whitney test between different prognostic factor statuses and molecular subtypes. Receiver operating characteristic curve analysis was used to assess the performance of mean and median histogram metrics in differentiating the molecular subtypes. A P value <0.05 was considered statistically significant.

RESULTS

The histogram metrics of ADC, D, and α differed significantly between ER-positive and ER-negative status, and between PR-positive and PR-negative status. The histogram metrics of ADC, D, α, and β were also significantly different between the HER2-positive and HER2-negative subgroups, and between ALNM-positive and ALNM-negative subgroups. The histogram metrics of α and β significantly differed between high and low Ki-67 proliferation subgroups, and between histological grade subgroups. The combination of αmean and βmean achieved the highest performance (AUC = 0.702) to discriminate the Luminal and HER2-positive subtypes.

DATA CONCLUSION

Whole-tumor histogram analysis of the CTRW model has potential to provide additional information on the prognosis and intrinsic subtyping classification of breast cancer.

EVIDENCE LEVEL

4 TECHNICAL EFFICACY: Stage 2.

Optimal selection of b-values for differential diagnosis of mediastinal lymph nodes using diffusion-weighted imaging

This study proposed to investigate the optimal selection of b-values in diffusion-weighted imaging for distinguishing malignant from benign mediastinal lymph nodes. Diffusion-weighted imaging with six b-values was performed on 35 patients at 1.5 T. Image quality score, signal-to-noise ratio, and relative contrast ratio of lymph node to chest muscle were compared between the diffusion-weighted images with a b-value up to 800 and 1000 s/mm². Using a lower and an upper b-value in the range of 0-1000 s/mm², eight apparent diffusion coefficient maps were obtained from a mono-exponential model. Receiver operating characteristic analysis was employed to evaluate the performance of the apparent diffusion coefficients for distinguishing malignant from benign mediastinal lymph nodes by using the area under the curve as a criterion. The mean image quality score and the relative contrast ratio showed no difference between b-values of 800 and 1000 s/mm². In the receiver operating characteristic analysis, the areas under the curve of apparent diffusion coefficient with b-value pairs of (0, 800), (0, 1000), and (50, 800) s/mm² were significantly higher than those from the other b-value pairs. No significant difference was observed among the three b-value pairs. Apparent diffusion coefficient obtained from b-value pairs of (0, 800), (0, 1000), and (50, 800) s/mm² showed superior diagnostic performance compared to the other b-value combinations. Based on several practical considerations, the b-value pair of (50, 800) s/mm² is recommended for differential diagnosis of mediastinal lymph nodes.

Characterization of breast lesions using multi-parametric diffusion MRI and machine learning

Objective. To investigate quantitative imaging markers based on parameters from two diffusion-weighted imaging (DWI) models, continuous-time random-walk (CTRW) and intravoxel incoherent motion (IVIM) models, for characterizing malignant and benign breast lesions by using a machine learning algorithm.Approach. With IRB approval, 40 women with histologically confirmed breast lesions (16 benign, 24 malignant) underwent DWI with 11b-values (50 to 3000 s/mm2) at 3T. Three CTRW parameters,Dm,α, andβand three IVIM parametersDdiff,Dperf, andfwere estimated from the lesions. A histogram was generated and histogram features of skewness, variance, mean, median, interquartile range; and the value of the 10%, 25% and 75% quantiles were extracted for each parameter from the regions-of-interest. Iterative feature selection was performed using the Boruta algorithm that uses the Benjamin Hochberg False Discover Rate to first determine significant features and then to apply the Bonferroni correction to further control for false positives across multiple comparisons during the iterative procedure. Predictive performance of the significant features was evaluated using Support Vector Machine, Random Forest, Naïve Bayes, Gradient Boosted Classifier (GB), Decision Trees, AdaBoost and Gaussian Process machine learning classifiers.Main Results. The 75% quantile, and median ofDm; 75% quantile off;mean, median, and skewness ofβ;kurtosis ofDperf; and 75% quantile ofDdiffwere the most significant features. The GB differentiated malignant and benign lesions with an accuracy of 0.833, an area-under-the-curve of 0.942, and an F1 score of 0.87 providing the best statistical performance (p-value < 0.05) compared to the other classifiers.Significance. Our study has demonstrated that GB with a set of histogram features from the CTRW and IVIM model parameters can effectively differentiate malignant and benign breast lesions.

Differentiation of benign and malignant breast lesions using diffusion-weighted imaging with a fractional-order calculus model

PURPOSE

To assess the feasibility of using three diffusion parameters (D, β, and μ) derived from fractional-order calculus (FROC) diffusion model for improving the differentiation between benign and malignant breast lesions.

METHOD

In this prospective study, 103 patients with breast lesions were enrolled. All subjects underwent diffusion-weighted imaging (DWI) with 12b values. Inter-observer agreement with respect to quantification of parameters by two radiologists was assessed using intraclass coefficient. Conventional apparent diffusion coefficient (ADC) and three FROC model parameters D, β, and μ were compared between the benign lesion and malignant lesion groups using the Mann-Whitney U test. Then, a comprehensive prediction model was created by using binary logistic regression. Receiver operating characteristic (ROC) curve analysis was performed to evaluate the diagnostic performance of the parameters using histopathological diagnosis as the reference standard.

RESULTS

The FROC parameters and ADC all exhibited significant differences between benign lesions and malignant lesions (P＜0.001). Among the individual parameters, the sensitivity of μ was higher than ADC (95.92% for μ vs 91.84% for ADC), and the specificity of β was higher than ADC (72.22% for β vs 70.37% for ADC). The combination of ADC and FROC parameters (D and β) generated the largest area under the ROC curve (0.841) when compared with individual parameters, indicating an improved performance for differentiating benign lesions from malignant lesions.

CONCLUSIONS

This study demonstrated the feasibility of using the FROC diffusion model to improve the accuracy of identifying malignant breast lesions.

Preoperative assessment of microvascular invasion of hepatocellular carcinoma using non-Gaussian diffusion-weighted imaging with a fractional order calculus model: A pilot study

PURPOSE

To assess the clinical potential of a set of new diffusion parameters (D, β, and μ) derived from fractional order calculus (FROC) diffusion model in predicting microvascular invasion (MVI) of hepatocellular carcinoma (HCC).

MATERIALS AND METHODS

Between January 2019 to November 2020, a total of 63 patients with HCC were enrolled in this study. Diffusion-weighted images were acquired by using ten b-values (0-2000 s/mm²). The FROC model parameters including diffusion coefficient (D), fractional order parameter (β), a microstructural quantity (μ) together with a conventional apparent diffusion coefficient (ADC) were calculated. Intraclass coefficients were calculated for assessing the agreement of parameters quantified by two radiologists. The differences of these values between the MVI-positive and MVI-negative HCC groups were compared by using independent sample t-test or the Mann-Whitney U test. Then the parameters showing significant differences between subgroups, including the β and D, were integrated to develop a comprehensive predictive model via binary logistic regression. The diagnostic performance was evaluated by receiver operating characteristic (ROC) analysis.

RESULTS

Among all the studied diffusion parameters, significant differences were found in D, β, and ADC between the MVI-positive and MVI-negative groups. MVI-positive HCCs showed significantly higher β values (0.65 ± 0.17 vs. 0.51 ± 0.13, P = 0.001), along with lower D values (0.84 ± 0.11 μm²/ms vs. 1.03 ± 0.13 μm²/ms, P < 0.001) and lower ADC values (1.38 ± 0.46 μm²/ms vs. 2.09 ± 0.70 μm²/ms, P < 0.001) than those of MVI-negative HCCs. According to the ROC analysis, the combination of D and β demonstrated the largest area under the ROC curve (0.920) compared with individual parameters (D: 0.912; β: 0.733; and ADC: 0.831) for differentiating MVI-positive from MVI-negative HCCs.

CONCLUSIONS

The FROC parameters can be used as noninvasive quantitative imaging markers for preoperatively predicting the MVI status of HCCs.

Incorporating multiple magnetic resonance diffusion models to differentiate low- and high-grade adult gliomas: a machine learning approach

BACKGROUND

Accurate grading of gliomas is a challenge in imaging diagnosis. This study aimed to evaluate the performance of a machine learning (ML) approach based on multiparametric diffusion-weighted imaging (DWI) in differentiating low- and high-grade adult gliomas.

METHODS

A model was developed from an initial cohort containing 74 patients with pathology-confirmed gliomas, who underwent 3 tesla (3T) diffusion magnetic resonance imaging (MRI) with 21 b values. In all, 112 histogram features were extracted from 16 parameters derived from seven diffusion models [monoexponential, intravoxel incoherent motion (IVIM), diffusion kurtosis imaging (DKI), fractional order calculus (FROC), continuous-time random walk (CTRW), stretched-exponential, and statistical]. Feature selection and model training were performed using five randomly permuted five-fold cross-validations. An internal test set (15 cases of the primary dataset) and an external cohort (n=55) imaged on a different scanner were used to validate the model. The diagnostic performance of the model was compared with that of a single DWI model and DWI radiomics using accuracy, sensitivity, specificity, and the area under the curve (AUC).

RESULTS

Seven significant multiparametric DWI features (two from the stretched-exponential and FROC models, and three from the CTRW model) were selected to construct the model. The multiparametric DWI model achieved the highest AUC (0.84, versus 0.71 for the single DWI model, P<0.05), an accuracy of 0.80 in the internal test, and both AUC and accuracy of 0.76 in the external test.

CONCLUSIONS

Our multiparametric DWI model differentiated low- (LGG) from high-grade glioma (HGG) with better generalization performance than the established single DWI model. This result suggests that the application of an ML approach with multiple DWI models is feasible for the preoperative grading of gliomas.

Value of non-Gaussian diffusion imaging with a fractional order calculus model combined with conventional MRI for differentiating histological types of cervical cancer

OBJECTIVE

This study aimed to evaluate the value of a fractional order calculus (FROC) model combined with conventional magnetic resonance imaging (MRI) for differentiating cervical adenocarcinoma (CAC) from squamous cell carcinoma (SCC).

METHODS

Diffusion-weighted imaging (DWI) with 9 b values (0-2000s/mm²) was carried out in 57 cervical cancer patients. Diffusion coefficient (D), fractional order parameter (β), and microstructural quantity (μ) together with apparent diffusion coefficient (ADC) were calculated and compared between the CAC and SCC groups. Conventional MRI features included T₂WI signal intensity (SI), unenhanced-T₁WI SI, enhanced-T₁WI SI, and ∆T₁WI SI, which were also compared between the two groups. Receiver operating characteristic (ROC) analysis was employed to assess the performance of FROC parameters, ADC, and conventional MRI features in differentiating CAC from SCC.

RESULTS

β was significantly lower in the CAC group than in the SCC group (0.682 ± 0.054 vs. 0.723 ± 0.084, P = 0.035), while D and μ were not significantly different between the two groups (D, P = 0.171; μ, P = 0.127). There was no significant difference in the ADC value between the two groups (P = 0.053). In conventional MRI features, enhanced-T₁WI SI was significantly higher in the SCC group than in the CAC group (985.78 ± 130.83 vs. 853.92 ± 149.65, P = 0.002). The area under the curve (AUC) of β, ADC, and enhanced-T₁WI SI was 0.700, 0.683, and 0.799, respectively. The combination of β, ADC, and enhanced-T₁WI SI revealed optimal diagnostic performance in differentiating CAC from SCC (AUC = 0.930), followed by β + enhanced-T₁WI SI (AUC = 0.869), ADC+ enhanced-T₁WI SI (AUC = 0.817), and β + ADC (AUC = 0.761).

CONCLUSION

The FROC model can serve as a noninvasive and quantitative imaging technique for differentiating CAC from SCC. β combined with ADC and enhanced-T₁WI SI had the highest diagnostic efficiency.

Staging Chronic Hepatitis B Related Liver Fibrosis with a Fractional Order Calculus Diffusion Model

RATIONALE AND OBJECTIVES

Accurately staging liver fibrosis is of great clinical significance. We aimed to evaluate the clinical potential of the non-Gaussian fractional order calculus (FROC) diffusion model in staging liver fibrosis.

MATERIALS AND METHODS

A total of 82 patients with chronic hepatitis B (CHB) were included in this prospective study. Diffusion weighted imaging (DWI)-derived parameters including the diffusion coefficient (D), fractional order parameter (β) and microstructural quantity (μ) sourced from FROC-DWI, and apparent diffusion coefficient (ADC) derived from mono-exponential DWI, as well as the aspartate aminotransferase-to-platelet ratio index (APRI) and fibrosis-4 (FIB-4) were calculated. Their correlations with fibrosis stages and the diagnostic efficacy in predicting liver fibrosis were assessed and compared.

RESULTS

D (r = -0.667), β (r = -0.671), μ (r = -0.481), and ADC (r = -0.665) displayed significant correlations with fibrosis stages (p < 0.001). D, β and ADC (p < 0.01) were independently associated with fibrosis; and compared to inflammatory activity, fibrosis was the independent factor significantly correlated with D, β and ADC (p < 0.001). There were no significant differences between the area under curves of D, β, μ or their combinations and ADC for predicting different fibrosis stages (p > 0.05). The diagnostic performance of the combined index with four diffusion metrics was better than D, β, μ or ADC used alone (p < 0.05) as well as APRI or FIB-4 (p < 0.01) in fibrosis staging.

CONCLUSION

FROC-DWI was valuable in staging liver fibrosis in patients with CHB, but there were no significant differences between the FROC-DWI parameters and the classical ADC. However, the combined DWI-derived index including D, β, μ and ADC offered the best diagnostic efficacy and may serve as a reliable tool for fibrosis evaluation, superior to APRI and FIB-4.

Cervical Carcinoma: Evaluation Using Diffusion MRI With a Fractional Order Calculus Model and its Correlation With Histopathologic Findings

Objective

The objective of the study is to investigate the feasibility of using the fractional order calculus (FROC) model to reflect tumor subtypes and histological grades of cervical carcinoma.

Methods

Sixty patients with untreated cervical carcinoma underwent multi-b-value diffusion-weighted imaging (DWI) at 3.0T magnetic resonance imaging (MRI). The mono-exponential and the FROC models were fitted. The differences in the histological subtypes and grades were evaluated by the Mann–Whitney U test. Receiver operating characteristic (ROC) analyses were performed to assess the diagnostic performance and to determine the best predictor for both univariate analysis and multivariate analysis. Differences between ROC curves were tested using the Hanley and McNeil test, while the sensitivity, specificity, and accuracy were compared using the McNemar test. P-value <0.05 was considered as significant difference. The Bonferroni corrections were applied to reduce problems associated with multiple comparisons.

Results

Only the parameter β, derived from the FROC model could differentiate cervical carcinoma subtypes (P = 0.03) and the squamous cell carcinoma (SCC) lesions exhibited significantly lower β than that in the adenocarcinoma (ACA) lesions. All the individual parameters, namely, ADC, β, D, and μ derived from the FROC model, could differentiate low-grade cervical carcinomas from high-grade ones (P = 0.022, 0.009, 0.004, and 0.015, respectively). The combination of all the FROC parameters showed the best overall performance, providing the highest sensitivity (81.2%) and AUC (0.829).

Conclusion

The parameters derived from the FROC model were able to differentiate the subtypes and grades of cervical carcinoma.

Machine Learning Based on Diffusion Kurtosis Imaging Histogram Parameters for Glioma Grading

Glioma grading plays an important role in surgical resection. We investigated the ability of different feature reduction methods in support vector machine (SVM)-based diffusion kurtosis imaging (DKI) histogram parameters to distinguish glioma grades. A total of 161 glioma patients who underwent magnetic resonance imaging (MRI) from January 2017 to January 2020 were included retrospectively. The patients were divided into low-grade (n = 61) and high-grade (n = 100) groups. Parametric DKI maps were derived, and 45 features from the DKI maps were extracted semi-automatically for analysis. Three feature selection methods [principal component analysis (PCA), recursive feature elimination (RFE) and least absolute shrinkage and selection operator (LASSO)] were used to establish the glioma grading model with an SVM classifier. To evaluate the performance of SVM models, the receiver operating characteristic (ROC) curves of SVM models for distinguishing glioma grades were compared with those of conventional statistical methods. The conventional ROC analysis showed that mean diffusivity (MD) variance, MD skewness and mean kurtosis (MK) C50 could effectively distinguish glioma grades, particularly MD variance. The highest classification distinguishing AUC was found using LASSO at 0.904 ± 0.069. In comparison, classification AUC by PCA was 0.866 ± 0.061, and 0.899 ± 0.079 by RFE. The SVM-PCA model with the lowest AUC among the SVM models was significantly better than the conventional ROC analysis (z = 1.947, p = 0.013). These findings demonstrate the superiority of DKI histogram parameters by LASSO analysis and SVM for distinguishing glioma grades.

Grading of clear cell renal cell carcinoma using diffusion MRI with a fractional order calculus model

BACKGROUND

The fractional order calculus (FROC) model has been developed to describe restrained motion of water molecules as well as microstructural heterogeneity, providing a novel tool for non-invasive tumor grading.

PURPOSE

To evaluate the role of the FROC model in characterizing clear cell renal cell carcinoma (ccRCC) grades.

MATERIAL AND METHODS

A total of 59 patients diagnosed with ccRCC were included in this prospective study. The diffusion metrics derived from the mono-exponential model (apparent diffusion coefficient [ADC]), intra-voxel incoherent motion [IVIM] model [D, D*, f], and FROC model [D_froc, β, μ]) were calculated and compared between low- and high-grade ccRCCs. Binary logistic regression analysis was performed to establish the diagnostic models. Receiver operating characteristic (ROC) analysis and DeLong test were performed to evaluate and compare the diagnostic performance of metrics in grading ccRCC.

RESULTS

All the metrics except D* and f exhibited statistical differences between low- and high-grade ccRCCs. ROC analysis showed individual FROC parameters, μ, D_froc, and β, outperformed ADC and IVIM parameters in grading ccRCC. For single parameter, μ demonstrated the highest AUC value, sensitivity, and diagnostic accuracy in discriminating the two ccRCC groups while β exhibited the optimal specificity. Importantly, the combination of D_froc, μ, and β could further improve the diagnostic performance.

CONCLUSION

The FROC parameters were superior to ADC and IVIM parameters in grading ccRCC, indicating the great potential of the FROC model in distinguishing low- and high-grade ccRCCs.

The relationship between diffusion heterogeneity and microstructural changes in high-grade gliomas using Monte Carlo simulations

PURPOSE

Diffusion-weighted imaging (DWI) may aid accurate tumor grading. Decreased diffusivity and increased diffusion heterogeneity measures have been observed in high-grade gliomas using the non-monoexponential models for DWI. However, DWI measures concerning tissue characteristics in terms of pathophysiological and structural changes are yet to be established. Thus, this study aims to investigate the relationship between the diffusion measurements and microstructural changes in the presence of high-grade gliomas using a three-dimensional Monte Carlo simulation with systematic changes of microstructural parameters.

METHODS

Water diffusion was simulated in a microenvironment along with changes associated with the presence of high-grade gliomas, including increases in cell density, nuclear volume, extracellular volume (VF_ex), and extracellular tortuosity (λ_ex), and changes in membrane permeability (P_mem). DWI signals were simulated using a pulsed gradient spin-echo sequence. The sequence parameters, including the maximum gradient strength and diffusion time, were set to be comparable to those of clinical scanners and advanced human MRI systems. The DWI signals were fitted using the gamma distribution and diffusional kurtosis models with b-values up to 6000 and 2500 s/mm², respectively.

RESULTS

The diffusivity measures (apparent diffusion coefficients (ADC), D_gamma of the gamma distribution model and D_app of the diffusional kurtosis model) decreased with increases in cell density and λ_ex, and a decrease in P_mem. These diffusivity measures increased with increases in nuclear volume and VF_ex. The diffusion heterogeneity measures (σ_gamma of the gamma distribution model and K_app of the diffusional kurtosis model) increased with increases in cell density or nuclear volume at the low P_mem, and a decrease in P_mem. Increased σ_gamma was also associated with an increase in VF_ex.

CONCLUSION

Among simulated microstructural changes, only increases in cell density at low P_mem or decreases in P_mem corresponded to both the decreased diffusivity and increased diffusion heterogeneity measures. The results suggest that increases in cell density at low P_mem or decreases in P_mem may be associated with the diffusion changes observed in high-grade gliomas.

Correlation of dual energy computed tomography electron density measurements with cerebral glioma grade

OBJECTIVE

To correlate dual energy computed tomography electron density measurements with histopathological cerebral glioma grading to determine whether it can be used as a non-invasive predictor of cerebral glioma grade.

MATERIALS AND METHODS

Fifty patients with suspected cerebral gliomas on imaging scheduled to undergo resection were included. We tested our hypothesis that with increasing glioma grade, increased tumor cellularity should translate into increased electron density and if a statistically significant difference between electron density of low-grade gliomas and high-grade gliomas is seen, we may have a clinical use of dual energy computed tomography as a non-invasive tool to predict cerebral glioma grade.A pre-operative dual energy computed tomography scan of the brain was performed, and electron density measurements calculated from the solid part of the tumor. Obtaining a ratio with electron density of contralateral normal brain parenchyma normalized these values. The minimum, maximum and mean electron density and their normalized values recorded between high-grade gliomas and low-grade gliomas were compared for presence of statistical significance.

RESULTS

A statistically significant difference was found between all six parameters recorded (minimum electron density and normalized values, mean electron density and normalized values, maximum electron density and normalized values) between low-grade gliomas and high-grade gliomas. The predictivity ranged from 75% (for minimum electron density and maximum normalized values) to 81.25% (for mean normalized values). All six parameters were found to have statistically significant positive correlation with Ki-67 index.

CONCLUSION

Dual energy computed tomography electron density measurements in cerebral gliomas are predictive of pre-operative differentiation of low-grade gliomas from high-grade gliomas and show a linear, statistically significant positive correlation with Ki-67 index.

Predicting the aggressiveness of peripheral zone prostate cancer using a fractional order calculus diffusion model

PURPOSE

To evaluate the performance of parameters D, β, μ from the Fractional Order Calculus (FROC) model at differentiating peripheral zone (PZ) prostate cancer (PCa) MATERIAL AND METHODS: 75 patients who underwent targeted MRI-guided TRUS prostate biopsy within 6 months of MRI were reviewed retrospectively. Regions of interest (ROI) were placed on suspicious lesions on MRI scans. ROIs were then correlated to pathological results based on core biopsy location. The final tumor count is a total: 23 of GS 6 (3 + 3), 36 of GS 7 (3 + 4), 18 of GS 7 (4 + 3), and 19 of GS ≥ 8. Diffusion-weighted imaging (DWI) scans were fitted into the FROC and monoexponential model to calculate ADC and FROC parameters: anomalous diffusion coefficient D, intravoxel diffusion heterogeneity β, and spatial parameter μ. The performance of FROC parameters and ADC at differentiating PCa grade was evaluated with receiver operating characteristic (ROC) analysis.

RESULTS

In differentiating low (GS 6) vs. intermediate (GS 7) risk PZ PCa, combination of (D, β) provides the best performance with AUC of 0.829 with significance of p = 0.018 when compared to ADC (AUC of 0.655). In differentiating clinically significant (GS 6) vs. clinically significant (GS ≥ 7) PCa, combination of (D, β, μ) provides highest AUC of 0.802 when compared to ADC (AUC of 0.671) with significance of p = 0.038. Stratification of intermediate (GS 7) and high (GS ≥ 8) risk PCa with FROC did not reach a significant difference when compared to ADC.

CONCLUSION

Combination of FROC parameters shows greater performance than ADC at differentiating low vs. intermediate risk and clinically insignificant vs. significant prostate cancers in peripheral zone lesions. The FROC diffusion model holds promise as a quantitative imaging technique for non-invasive evaluation of PZ PCa.

Evaluation of a fractional-order calculus diffusion model and bi-parametric VI-RADS for staging and grading bladder urothelial carcinoma

OBJECTIVES

To evaluate the feasibility of high b-value diffusion-weighted imaging (DWI) for distinguishing non-muscle-invasive bladder cancer (NMIBC) from muscle-invasive bladder cancer (MIBC) and low- from high-grade bladder urothelial carcinoma using a fractional-order calculus (FROC) model as well as a combination of FROC DWI and bi-parametric Vesical Imaging-Reporting and Data System (VI-RADS).

METHODS

Fifty-eight participants with bladder urothelial carcinoma were included in this IRB-approved prospective study. Diffusion-weighted images, acquired with 16 b-values (0-3600 s/mm²), were analyzed using the FROC model. Three FROC parameters, D, β, and μ, were used for delineating NMIBC from MIBC and for tumor grading. A receiver operating characteristic (ROC) analysis was performed based on the individual FROC parameters and their combinations, followed by comparisons with apparent diffusion coefficient (ADC) and bi-parametric VI-RADS based on T2-weighted images and DWI.

RESULTS

D and μ were significantly lower in the MIBC group than in the NMIBC group (p = 0.001 for each), and D, β, and μ all exhibited significantly lower values in the high- than in the low-grade tumors (p ≤ 0.011). The combination of D, β, and μ produced the highest specificity (85%), accuracy (78%), and the area under the ROC curve (AUC, 0.782) for distinguishing NMIBC and MIBC, and the best sensitivity (89%), specificity (86%), accuracy (88%), and AUC (0.892) for tumor grading, all of which outperformed the ADC. The combination of FROC parameters with bi-parametric VI-RADS improved the AUC from 0.859 to 0.931.

CONCLUSIONS

High b-value DWI with a FROC model is useful in distinguishing NMIBC from MIBC and grading bladder tumors.

KEY POINTS

• Diffusion parameters derived from a FROC diffusion model may differentiate NMIBC from MIBC and low- from high-grade bladder urothelial carcinomas. • Under the condition of a moderate sample size, higher AUCs were achieved by the FROC parameters D (0.842) and μ (0.857) than ADC (0.804) for bladder tumor grading with p ≤ 0.046. • The combination of the three diffusion parameters from the FROC model can improve the specificity over ADC (85% versus 67%, p = 0.031) for distinguishing NMIBC and MIBC and enhance the performance of bi-parametric VI-RADS.

Comparison of mono-exponential, bi-exponential, kurtosis, and fractional-order calculus models of diffusion-weighted imaging in characterizing prostate lesions in transition zone

PURPOSE

To compare various models of diffusion-weighted imaging including mono-exponential, bi-exponential, diffusion kurtosis (DK) and fractional-order calculus (FROC) models in diagnosing prostate cancer (PCa) in transition zone (TZ) and distinguish the high-grade PCa [Gleason score (GS) ≥ 7] lesions from the total of low-grade PCa (GS ≤ 6) lesions and benign prostatic hyperplasia (BPH) in TZ.

METHODS

80 Patients with 103 lesions were included in this study. Nine metrics [including apparent diffusion coefficient (ADC) derived from mono-exponential model, slow diffusion coefficient (D_s), fast diffusion coefficient (D_f),, and f (the fraction of fast diffusion) from bi-exponential model; mean diffusivity (MD) and mean kurtosis (MK) from DK model; diffusion coefficient (D), fractional-order derivative in space (β), and spatial metric (μ) from FROC model] were calculated. Comparisons between BPH and PCa lesions as well as between clinically significant PCa (CsPCa) (GS ≥ 7, n = 31) and clinically insignificant lesions (Cins) (GS ≤ 6 and BPH, n = 72) of these metrics were conducted. Mann-Whitney U-test and receiver operating characteristic (ROC) analysis were used for statistical evaluations.

RESULTS

The areas under the ROC curve (AUC) values of β derived from FROC model were 0.778 and 0.853 in differentiating PCa from BPH and in differentiating CS (GS ≥ 7) from Cins (GS ≤ 6 and BPH), both were the highest compared to other metrics. The AUC value of β was significantly higher than that of ADC (P = 0.009) in differentiating CS from Cins, while the differentiation between BPH and PCa did not reach the statistical significance when comparing with ADC (P = 0.089).

CONCLUSION

Although no significant difference was found in distinguishing PCa from BPH, the metric β derived from FROC model was superior to other diffusion metrics in differentiation between CS and Cins in TZ.

Quartile histogram assessment of glioma malignancy using high b-value diffusion MRI with a continuous-time random-walk model

The purpose of this study is to investigate the feasibility of using a continuous-time random-walk (CTRW) diffusion model, together with a quartile histogram analysis, for assessing glioma malignancy by probing tissue heterogeneity as well as cellularity. In this prospective study, 91 patients (40 females, 51 males) with histopathologically proven gliomas underwent MRI at 3 T. The cohort included 42 grade II (GrII), 19 grade III (GrIII) and 29 grade IV (GrIV) gliomas. Echo-planar diffusion-weighted imaging was conducted using 17 b-values (0-4000 s/mm² ). Three CTRW model parameters, including an anomalous diffusion coefficient D_m , and two parameters related to temporal and spatial diffusion heterogeneity α and β, respectively, were obtained. The mean parameter values within the tumor regions of interest (ROIs) were computed by utilizing the first quartile of the histograms as well as the full ROI for comparison. A Bonferroni-Holm-corrected Mann-Whitney U-test was used for the group comparisons. Individual and combinations of the CTRW parameters were evaluated for the characterization of gliomas with a receiver operating characteristic analysis. All first-quartile mean CTRW parameters yielded significant differences (p-values < 0.05) between pair-wise comparisons of GrII (D_m : 1.14 ± 0.37 μm² /ms; α: 0.904 ± 0.03, β: 0.913 ± 0.06), GrIII (D_m : 0.88 ± 0.21 μm² /ms; α: 0.888 ± 0.01, β: 0.857 ± 0.06) and GrIV gliomas (D_m : 0.73 ± 0.22 μm² /ms; α: 0.878 ± 0.01; β: 0.791 ± 0.07). The highest sensitivity, specificity, accuracy and area-under-the-curve of using the combinations of the first-quartile parameters were 84.2%, 78.5%, 75.4% and 0.76 for GrII and GrIII classification; 86.2%, 89.4%, 75% and 0.76 for GrIII and GrIV classification; and 86.2%, 85.7%, 84.5% and 0.90 for GrII and GrIV classification, respectively. Quartile-based analysis produced higher accuracy and area-under-the-curve than the full ROI-based analysis in all classifications. The CTRW diffusion model, together with a quartile-based histogram analysis, offers a new way for probing tumor structural heterogeneity at a subvoxel level, and has potential for in vivo assessment of glioma malignancy to complement histopathology.

Altered Motor Excitability in Patients With Diffuse Gliomas Involving Motor Eloquent Areas: The Impact of Tumor Grading

BACKGROUND

Diffuse gliomas have an increased biological aggressiveness across the World Health Organization (WHO) grading system. The implications of glioma grading on the primary motor cortex (M1)-corticospinal tract (CST) excitability is unknown.

OBJECTIVE

To assess the excitability of the motor pathway with navigated transcranial magnetic stimulation (nTMS).

METHODS

Retrospective cohort study of patients admitted for surgery with diffuse gliomas within motor eloquent areas. Demographic, clinical, and nTMS-related variables were collected. The Cortical Excitability Score (CES 0 to 2 according to the number of abnormal interhemispheric resting motor threshold (RMT) ratios) was calculated for patients where bilateral upper and lower limb mapping was performed.

RESULTS

A total of 45 patients were included: 9 patients had a low-grade glioma and 36 patients had a high-grade glioma. The unadjusted analysis revealed an increase in the latency of the motor evoked potential of the lower limb with an increase of the WHO grade (P = .038). The adjusted analysis confirmed this finding (P = .013) and showed a relation between the increase in the WHO and a decreased RMT (P = .037) of the motor evoked responses in the lower limb. When CES was calculated, an increase in the score was related with an increase in the WHO grade (unadjusted analysis-P = .0001; adjusted analysis-P = .001) and in isocitrate dehydrogenase (IDH) wild-type tumors (unadjusted analysis-P = .020).

CONCLUSION

An increase in the WHO grading system and IDH wild-type tumors are associated with an abnormal excitability of the motor eloquent areas in patients with diffuse gliomas.

In vivo assessment of Lauren classification for gastric adenocarcinoma using diffusion MRI with a fractional order calculus model

OBJECTIVES

To evaluate the performance of a fractional order calculus (FROC) diffusion model for imaging-based assessment of Lauren classification in gastric adenocarcinoma.

METHODS

In this study, 43 patients (15 females, 28 males) with gastric adenocarcinoma underwent MRI at 1.5 T. According to pathology-based Lauren classification, 10 patients had diffuse-type, 20 had intestinal-type, and 13 had mixed-type lesions. The diffuse and mixed types were combined as diffuse-and-mixed type to be differentiated from the intestinal type using diffusion MRI. Diffusion-weighted images were acquired by using eleven b-values (0-2000 s/mm²). Three FROC model parameters comprising diffusion coefficient D, intravoxel diffusion heterogeneity β, and a microstructural quantity μ, together with a conventional apparent diffusion coefficient (ADC), were estimated. The mean parameter values in the tumour were computed by using a percentile histogram analysis. Individual or linear combinations of the mean parameters in the tumour were used to differentiate the diffuse-and-mixed type from the intestinal type using descriptive statistics and receiver operating characteristic (ROC) analyses.

RESULTS

Significant differences were observed between diffuse-and-mixed-type and intestinal-type lesions in D (0.99 ± 0.20 μm²/ms vs. 1.11 ± 0.23 μm²/ms; p = 0.036), β (0.37 ± 0.08 vs. 0.43 ± 0.11; p = 0.043), μ (7.92 ± 2.79 μm vs. 9.87 ± 1.52 μm; p = 0.038), and ADC (0.81 ± 0.34 μm²/ms vs. 0.96 ± 0.19 μm²/ms; p = 0.033). Among the individual parameters, μ produced the largest area under the ROC curve (0.739). The combinations of (D, β, μ) and (β and μ) produced the best overall performance with a sensitivity of 0.739, specificity of 0.750, accuracy of 0.744, and area under the curve of 0.793 (95% confidence interval: 0.657-0.929).

CONCLUSION

Diffusion MRI with the FROC model holds promise for non-invasive assessment of Lauren classification for gastric adenocarcinoma.

KEY POINTS

• High b-value diffusion MRI with a FROC model that is sensitive to tissue microstructures can differentiate the diffuse-and-mixed type from intestinal type of gastric adenocarcinoma. • The combination of FROC parameters produced the best result for distinguishing the diffuse-and-mixed type from the intestinal type with an area under the receiver operating characteristic curve of 0.793. • The FROC model parameters, individually or conjointly, hold promise for repeated, non-invasive evaluations of gastric adenocarcinoma at various time points throughout disease progression or regression to complement conventional Lauren classification.

Multi-attribute decision-making method based on a novel distance measure of linguistic intuitionistic fuzzy sets

Linguistic intuitionistic fuzzy sets can qualitatively rather than quantitatively express data in the form of membership degree. But quantitative tools are required to handle qualitative information. Therefore, an improved linguistic scale function, which can more accurately manifest the subjective feelings of decision-makers, is employed to deal with linguistic intuitionistic information. Subsequently, due to some commonly used distance measures do not comprehensively evaluate the information of linguistic intuitionistic fuzzy sets, an improved distance measure of linguistic intuitionistic fuzzy sets is designed. It considers the cross-evaluation information to get more realistic reasoning results. In addition, a new similarity measure defined by nonlinear Gaussian diffusion model is proposed, which can provide different response scales for different information between various schemes. The properties of these measures are also studied in detail. On this basis, a method in linguistic intuitionistic fuzzy environment is developed to handle multi-attribute decision-making problems. Finally, an illustrative example is given to demonstrate the effectiveness of the proposed method and the influence of the parameters is analyzed.

Differentiation of salivary gland tumor using diffusion-weighted imaging with a fractional order calculus model

OBJECTIVE

To evaluate the feasibility of using imaging parameters (D, β and μ) obtained from fractional order calculus (FROC) diffusion model to differentiate salivary gland tumors.

METHODS

15 b-value (0-2000 s/mm²) diffusion-weighted imaging (DWI) was scanned in 62 patients with salivary gland tumors (47 benign and 15 malignant). Diffusion coefficient D, fractional order parameter β (which correlates with tissue heterogeneity) and a microstructural quantity μ of the solid portion within the tumor were calculated, and compared between benign and malignant groups, or among pleomorphic adenoma (PA), Warthin's tumor (WT), and malignant tumor (MT) groups. Performance of FROC parameters for differentiation was assessed using receiver operating characteristic analysis.

RESULTS

None of the FROC parameters exhibited significant differences between benign and malignant group (D, p = 0.150; β, p = 0.967; μ, p = 0.693). WT showed significantly lower D (p < 0.001) and β (p < 0.001), while higher μ (p = 0.001) than PA. Combination of D, β and μ showed optimal diagnostic performance (area under the curve, AUC, 0.998). MT showed significantly lower D (p = 0.001) and β (p = 0.025) than PA, while no significant difference was found on μ (p = 0.064). Combination of D and β showed optimal diagnostic performance (AUC, 0.933). Significant difference was found on β (p = 0.027) between MT and WT, while not on D (p = 0.806) and μ (p = 0.789). Setting a βof 0.615 as the cut-off value, optimal diagnostic performance could be obtained (AUC = 0.806).

CONCLUSION

A non-Gaussian FROC diffusion model can serve as a noninvasive and quantitative imaging technique for differentiating salivary gland tumors.

ADVANCES IN KNOWLEDGE

(1) PA showed higher D and β and lower μ than WT. (2) PA had higher D and β than MT. (3) WT demonstrated lower β than MT. (4) β, as a new FROC parameter, could offer an added value to the differentiation.

High-Spatial-Resolution Diffusion MRI in Parkinson Disease: Lateral Asymmetry of the Substantia Nigra

Background Motor symptoms in Parkinson disease (PD) have exhibited lateral asymmetry, suggesting asymmetric neuronal loss in the substantia nigra (SN). Diffusion MRI may be able to help confirm tissue microstructural alterations in the substantia nigra to probe for the presence of asymmetry. Purpose To investigate lateral asymmetry in the SN of patients with PD by using diffusion MRI with both Gaussian and non-Gaussian models. Materials and Methods In this cross-sectional study conducted from March 2015 to March 2017, 27 participants with PD and 27 age-matched healthy control (HC) participants, all right handed, underwent MRI at 3.0 T. High-spatial-resolution diffusion images were acquired with a reduced field of view by using seven b values up to 3000 sec/mm². A continuous-time random-walk (CTRW) non-Gaussian diffusion model was used to produce anomalous diffusion coefficient (D_m) and temporal (α) and spatial (β) diffusion heterogeneity indexes followed by a Gaussian diffusion model to yield an apparent diffusion coefficient (ADC). Individual or linear combinations of diffusion parameters in the SN were unilaterally and bilaterally compared between the PD and HC groups. Results In the bilateral comparison between the PD and HC groups, differences were observed in β (0.67 ± 0.06 [standard deviation] vs 0.64 ± 0.04, respectively; P = .016), ADC (0.48 μm²/msec ± 0.08 vs 0.53 μm²/msec ± 0.06, respectively; P = .03), and the combination of CTRW parameters (P = .02). In the unilateral comparison, differences were observed in all diffusion parameters on the left SN (P < .03), but not on the right (P > .20). In a receiver operating characteristic (ROC) analysis to delineate left SN abnormality in PD, the combination of D_m, α, and β produced the best sensitivity (sensitivity, 0.78); the combination of D_m and β produced the best specificity (specificity, 0.85); and the combination of α and β produced the largest area under the ROC curve (area under the ROC curve, 0.73). Conclusion These results suggest that quantitative diffusion MRI is sensitive to brain tissue changes in participants with Parkinson disease and provide evidence of substantia nigra lateral asymmetry in this disease. © RSNA, 2019 Online supplemental material is available for this article.

A fractional motion diffusion model for a twice-refocused spin-echo pulse sequence

The purpose of this study was to develop an analytical expression for a fractional motion (FM) diffusion model to characterize diffusion-induced signal attenuation in a twice-refocused spin-echo (TRSE) sequence that is resilient to eddy currents, and to demonstrate its applicability to human brain imaging in vivo. Based on the FM theory, which provides a unified statistical description for Langevin motions, the diffusion-weighted (DW) MR signal was measured with a TRSE sequence that balances the concomitant gradients. The analytical expression was fitted to a set of DW images acquired with 14 b-values (0-4000 s/mm² ) from a total of 10 healthy human subjects at 3 T, yielding three FM parameter maps based on anomalous diffusion coefficient D_{φ, ψ} , diffusion increment variance φ, and diffusion correlation ψ, respectively. These parameters were used to characterize different brain regions in gray matter (GM), white matter (WM), and cerebrospinal fluid. The analytical expression for the TRSE-based FM model accurately described diffusion signal attenuation in healthy brain tissues at high b-values. TRSE's robustness against eddy currents was illustrated by comparing results from an expression for a conventional Stejskal-Tanner sequence. The TRSE-based FM model also produced consistent GM-WM contrast (p < 0.01) across all brain regions studied, whereas the consistency was not observed with the Stejskal-Tanner-based FM model. This new analytical expression is expected to enable further investigations to probe tissue structures by exploiting anomalous diffusion properties without being hindered by eddy-current perturbations at high b-values.

Diffusion MRI of cancer: From low to high b-values

Following its success in early detection of cerebral ischemia, diffusion-weighted imaging (DWI) has been increasingly used in cancer diagnosis and treatment evaluation. These applications are propelled by the rapid development of novel diffusion models to extract biologically valuable information from diffusion-weighted MR signals, and significant advances in MR hardware that has enabled image acquisition with high b-values. This article reviews recent technical developments and clinical applications in cancer imaging using DWI, with a special emphasis on high b-value diffusion models. The article is organized in four sections. First, we provide an overview of diffusion models that are relevant to cancer imaging. The model parameters are discussed in relation to three tissue properties-cellularity, vascularity, and microstructures. An emphasis is placed on characterization of microstructural heterogeneity, given its novelty and close relevance to cancer. Second, we illustrate diffusion MR clinical applications in each of the following three categories: 1) cancer detection and diagnosis; 2) cancer grading, staging, and classification; and 3) cancer treatment response prediction and evaluation. Third, we discuss several practical issues, including selection of image acquisition parameters, reproducibility and reliability, motion management, image distortion, etc., that are commonly encountered when applying DWI to cancer in clinical settings. Lastly, we highlight a few ongoing challenges and provide some possible future directions, particularly in the area of establishing standards via well-organized multicenter clinical trials to accelerate clinical translation of advanced DWI techniques to improving cancer care on a large scale. Level of Evidence: 5 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2019;49:23-40.

White matter structural differences in OSA patients experiencing residual daytime sleepiness with high CPAP use: a non-Gaussian diffusion MRI study

OBJECTIVES

To investigate factors associated with residual sleepiness in patients who were highly adherent to continuous positive airway pressure (CPAP). Nocturnal inactivity, comorbidities, concomitant medications, and, in particular, white matter (WM) differences using diffusion magnetic resonance imaging (MRI) were explored using a continuous-time random-walk (CTRW) model.

METHODS

Twenty-seven male patients (30-55 years of age) with obstructive sleep apnea (OSA) received CPAP as the only treatment (CPAP ≥ 6 h/night) for at least 30 days. Based on the Psychomotor Vigilance Task (PVT) results, participants were divided into a non-sleepy group (lapses ≤ 5; n = 18) and a sleepy group (lapses > 5; n = 9). Mean nocturnal inactivity (sleep proxy) was measured using actigraphy for one week. Diffusion-weighted imaging (DWI) with high b-values, as well as diffusion tensor imaging (DTI), was performed on a 3 T MRI scanner. The DWI dataset was analyzed using the CTRW model that yielded three parameters: temporal diffusion heterogeneity α, spatial diffusion heterogeneity β, and an anomalous diffusion coefficient D_m. The differences in α, β, and D_m between the two groups were investigated by a whole-brain analysis using tract-based spatial statistics (TBSS), followed by a regional analysis on individual fiber tracts using a standard parcellation template. Results from the CTRW model were compared with those obtained from DTI. The three CTRW parameters were also correlated with the clinical assessment scores, Epworth Sleepiness Scale (ESS), PVT lapses, and PVT mean reaction time (MRT) in specific fiber tracts.

RESULTS

There were no differences between groups in mean sleep duration, comorbidities, and the number or type of medications, including alerting and sedating medications. In the whole-brain DWI analysis, the sleepy group showed higher α (17.27% of the WM voxels) and D_m (17.14%) when compared to the non-sleepy group (P < 0.05), whereas no significant difference in β was observed. In the regional fiber analysis, the sleepy and non-sleepy groups showed significant differences in α, β, or their combinations in a total of 12 fiber tracts; whereas similar differences were not observed in DTI parameters, when age was used as a covariate. Additionally, moderate to strong correlations between the CTRW parameters (α, β, or D_m) and the sleepiness assessment scores (ESS, PVT lapses, or PVT MRT) were observed in specific fiber tracts (|R| = 0.448-0.654, P = 0.0003-0.019).

CONCLUSIONS

The observed differences in the CTRW parameters between the two groups indicate that WM alterations can be a possible mechanism to explain reversible versus residual sleepiness observed in OSA patients with identical high level of CPAP use. The moderate to strong correlations between the CTRW parameters and the clinical scores suggest the possibility of developing objective and quantitative imaging markers to complement clinical assessment of OSA patients.

Anisotropy of anomalous diffusion improves the accuracy of differentiating low- and high-grade cerebral gliomas

BACKGROUND

Anomalous diffusion model has been introduced and shown to be beneficial in clinical applications. However, only the directionally averaged values of anomalous diffusion parameters were investigated, and the anisotropy of anomalous diffusion remains unexplored. The aim of this study was to demonstrate the feasibility of using anisotropy of anomalous diffusion for differentiating low- and high-grade cerebral gliomas.

METHODS

Diffusion MRI images were acquired from brain tumor patients and analyzed using the fractional motion (FM) model. Twenty-two patients with histopathologically confirmed gliomas were selected. An anisotropy metric for the FM-related parameters, including the Noah exponent (α) and the Hurst exponent (H), was introduced and their values were statistically compared between the low- and high-grade gliomas. Additionally, multivariate logistic regression analysis was performed to assess the combination of the anisotropy metric and the directionally averaged value for each parameter. The diagnostic performances for grading gliomas were evaluated using a receiver operating characteristic (ROC) analysis.

RESULTS

The Hurst exponent H was more anisotropic in high-grade than in low-grade gliomas (P = 0.015), while no significant difference was observed for the anisotropy of α. The ROC analysis revealed that larger areas under the ROC curves were produced for the combination of α (1) and the combination of H (0.813) compared with the directionally averaged α (0.979) and H (0.594), indicating an improved performance for tumor differentiation.

CONCLUSION

The anisotropy of anomalous diffusion can provide distinctive information and benefit the differentiation of low- and high-grade gliomas. The utility of anisotropic anomalous diffusion may have an improved effect for investigating pathological changes in tissues.

Texture analysis- and support vector machine-assisted diffusional kurtosis imaging may allow in vivo gliomas grading and IDH-mutation status prediction: a preliminary study

We sought to investigate, whether texture analysis of diffusional kurtosis imaging (DKI) enhanced by support vector machine (SVM) analysis may provide biomarkers for gliomas staging and detection of the IDH mutation. First-order statistics and texture feature extraction were performed in 37 patients on both conventional (FLAIR) and mean diffusional kurtosis (MDK) images and recursive feature elimination (RFE) methodology based on SVM was employed to select the most discriminative diagnostic biomarkers. The first-order statistics demonstrated significantly lower MDK values in the IDH-mutant tumors. This resulted in 81.1% accuracy (sensitivity = 0.96, specificity = 0.45, AUC 0.59) for IDH mutation diagnosis. There were non-significant differences in average MDK and skewness among the different tumour grades. When texture analysis and SVM were utilized, the grading accuracy achieved by DKI biomarkers was 78.1% (sensitivity 0.77, specificity 0.79, AUC 0.79); the prediction accuracy for IDH mutation reached 83.8% (sensitivity 0.96, specificity 0.55, AUC 0.87). For the IDH mutation task, DKI outperformed significantly the FLAIR imaging. When using selected biomarkers after RFE, the prediction accuracy achieved 83.8% (sensitivity 0.92, specificity 0.64, AUC 0.88). These findings demonstrate the superiority of DKI enhanced by texture analysis and SVM, compared to conventional imaging, for gliomas staging and prediction of IDH mutational status.

Non-Gaussian diffusion imaging with a fractional order calculus model to predict response of gastrointestinal stromal tumor to second-line sunitinib therapy

PURPOSE

To demonstrate the clinical value of a non-Gaussian diffusion model using fractional order calculus (FROC) for early prediction of the response of gastrointestinal stromal tumor to second-line sunitinib targeted therapy.

METHODS

Fifteen patients underwent sunitinib treatment after imatinib resistance. Diffusion-weighted imaging with multiple b-values was performed before treatment (baseline) and 2 weeks (for early prediction of response) after initiating sunitinib treatment. Conventional MRI images at 12 weeks were used to determine the good and poor responders according to the modified Choi criteria for MRI. Diffusion coefficient D, fractional order parameter β (which correlates to intravoxel tissue heterogeneity), and a microstructural quantity µ were calculated using the FROC model. The FROC parameters and the longest diameter of the lesion, as well as their changes after 2 weeks of treatment, were compared between the good and poor responders. Additionally, the pretreatment FROC parameters were individually combined with the change in D (ΔD) using a logistic regression model to evaluate response to sunitinib treatment with a receiver operating characteristic analysis.

RESULTS

Forty-two good-responding and 32 poor-responding lesions were identified. Significant differences were detected in pretreatment β (0.67 versus 0.74, P = 0.011) and ΔD (45.7% versus 12.4%, P = 0.001) between the two groups. The receiver operating characteristic analysis showed that ΔD had a significantly higher predictive power than the tumor size change (area under the curve: 0.725 versus 0.580; 0.95 confidence interval). When ΔD was combined with pretreatment β, the area under the curve improved to 0.843 with a predictive accuracy of 75.7% (56 of 74).

CONCLUSIONS

The non-Gaussian FROC diffusion model showed clinical value in early prediction of gastrointestinal stromal tumor response to second-line sunitinib targeted therapy. The pretreatment FROC parameter β can increase the predictive accuracy when combined with the change in diffusion coefficient during treatment. Magn Reson Med 79:1399-1406, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

On a fractional order calculus model in diffusion weighted breast imaging to differentiate between malignant and benign breast lesions detected on X-ray screening mammography

OBJECTIVE

To evaluate a fractional order calculus (FROC) model in diffusion weighted imaging to differentiate between malignant and benign breast lesions in breast cancer screening work-up using recently introduced parameters (βFROC, DFROC and μFROC).

MATERIALS AND METHODS

This retrospective analysis within a prospective IRB-approved study included 51 participants (mean 58.4 years) after written informed consent. All patients had suspicious screening mammograms and indication for biopsy. Prior to biopsy, full diagnostic contrast-enhanced MRI examination was acquired including diffusion-weighted-imaging (DWI, b = 0,100,750,1500 s/mm2). Conventional apparent diffusion coefficient Dapp and FROC parameters (βFROC, DFROC and μFROC) as suggested further indicators of diffusivity components were measured in benign and malignant lesions. Receiver operating characteristics (ROC) were calculated to evaluate the diagnostic performance of the parameters.

RESULTS

29/51 patients histopathologically revealed malignant lesions. The analysis revealed an AUC for Dapp of 0.89 (95% CI 0.80-0.98). For FROC derived parameters, AUC was 0.75 (0.60-0.89) for DFROC, 0.59 (0.43-0.75) for βFROC and 0.59 (0.42-0.77) for μFROC. Comparison of the AUC curves revealed a significantly higher AUC of Dapp compared to the FROC parameters DFROC (p = 0.009), βFROC (p = 0.003) and μFROC (p = 0.001).

CONCLUSION

In contrast to recent description in brain tumors, the apparent diffusion coefficient Dapp showed a significantly higher AUC than the recently proposed FROC parameters βFROC, DFROC and μFROC for differentiating between malignant and benign breast lesions. This might be related to the intrinsic high heterogeneity within breast tissue or to the lower maximal b-value used in our study.

A fractional motion diffusion model for grading pediatric brain tumors

Objectives

To demonstrate the feasibility of a novel fractional motion (FM) diffusion model for distinguishing low- versus high-grade pediatric brain tumors; and to investigate its possible advantage over apparent diffusion coefficient (ADC) and/or a previously reported continuous-time random-walk (CTRW) diffusion model.

Materials and methods

With approval from the institutional review board and written informed consents from the legal guardians of all participating patients, this study involved 70 children with histopathologically-proven brain tumors (30 low-grade and 40 high-grade). Multi-b-value diffusion images were acquired and analyzed using the FM, CTRW, and mono-exponential diffusion models. The FM parameters, D_fm, φ, ψ (non-Gaussian diffusion statistical measures), and the CTRW parameters, D_m, α, β (non-Gaussian temporal and spatial diffusion heterogeneity measures) were compared between the low- and high-grade tumor groups by using a Mann-Whitney-Wilcoxon U test. The performance of the FM model for differentiating between low- and high-grade tumors was evaluated and compared with that of the CTRW and the mono-exponential models using a receiver operating characteristic (ROC) analysis.

Results

The FM parameters were significantly lower (p < 0.0001) in the high-grade (D_fm: 0.81 ± 0.26, φ: 1.40 ± 0.10, ψ: 0.42 ± 0.11) than in the low-grade (D_fm: 1.52 ± 0.52, φ: 1.64 ± 0.13, ψ: 0.67 ± 0.13) tumor groups. The ROC analysis showed that the FM parameters offered better specificity (88% versus 73%), sensitivity (90% versus 82%), accuracy (88% versus 78%), and area under the curve (AUC, 93% versus 80%) in discriminating tumor malignancy compared to the conventional ADC. The performance of the FM model was similar to that of the CTRW model.

Conclusions

Similar to the CTRW model, the FM model can improve differentiation between low- and high-grade pediatric brain tumors over ADC.

•

The fractional motion (FM) diffusion model was applied to pediatric brain tumors.

•

The FM model parameters can be sensitive to tissue microstructures.

•

The FM model outperforms the mono-exponential diffusion model.

•

The FM model performs similarly to the continuous-time random-walk (CTRW) model.

•

Our results challenge those from recent biophysics studies in cell cultures.