Utility of a Hybrid IVIM-DKI Model to Predict the Development of Distant Metastasis in Head and Neck Squamous Cell Carcinoma Patients

Patterns of Distant Metastasis in Head and Neck Cancer in a Tertiary Care Centre

Objectives/Hypothesis To find out the frequency and location of distant metastasis in head and neck malignancies. Our study also aims to find out the most common site leading to distant metastasis and the management of these distant metastasis cases. Methods 1558 patients treated for head and neck malignancy between 2017 and 2021 were retrospectively reviewed. The frequency and proportions were used to produce descriptive statistics. Results The highest number of head and neck malignancy cases were reported in the oral cavity which included 943 cases (60.52%). Patients with distant metastasis (M1) accounted for 4.73 percent of all cases (n = 90). Nasopharyngeal malignancy cases showed the highest M1 frequency (29.03%), whereas oral cavity patients had the lowest frequency (2.75%). The most common site of distant metastasis was in the lung (64%) followed by bone (18%) and the liver (11%). CT scan of the neck and thorax was the most commonly used diagnostic modality. The most common histopathological finding was squamous cell carcinoma (85%). Multimodality treatment was employed for most of the detected cases. Conclusion Distant metastasis at presentation is rare in head and neck cancer. The rate of distant metastasis in the present study was 4.73%, with the lung being the most common site. The overall survival of these patients depends on a variety of factors and more studies are needed in this regard.

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.

Registration and quantification network (RQnet) for IVIM-DKI analysis in MRI

PURPOSE

A deep learning method is proposed for aligning diffusion weighted images (DWIs) and estimating intravoxel incoherent motion-diffusion kurtosis imaging parameters simultaneously.

METHODS

We propose an unsupervised deep learning method that performs 2 tasks: registration and quantification for intravoxel incoherent motion-diffusion kurtosis imaging analysis. A common registration method in diffusion MRI is based on minimizing dissimilarity between various DWIs, which may result in registration errors due to different contrasts in different DWIs. We designed a novel unsupervised deep learning method for both accurate registration and quantification of various diffusion parameters. In order to generate motion-simulated training data and test data, 17 volunteers were scanned without moving their heads, and 4 volunteers moved their heads during the scan in a 3 Tesla MRI. In order to investigate the applicability of the proposed method to other organs, kidney images were also obtained. We compared the registration accuracy of the proposed method, statistical parametric mapping, and a deep learning method with a normalized cross-correlation loss. In the quantification part of the proposed method, a deep learning method that considered the diffusion gradient direction was used.

RESULTS

Simulations and experimental results showed that the proposed method accurately performed registration and quantification for intravoxel incoherent motion-diffusion kurtosis imaging analysis. The registration accuracy of the proposed method was high for all b values. Furthermore, quantification performance was analyzed through simulations and in vivo experiments, where the proposed method showed the best performance among the compared methods.

CONCLUSION

The proposed method aligns the DWIs and accurately quantifies the intravoxel incoherent motion-diffusion kurtosis imaging parameters.

Intravoxel incoherent motion imaging combined with diffusion kurtosis imaging to assess the response to radiotherapy in a rabbit VX2 malignant bone tumor model

Purpose

To combine intravoxel incoherent motion (IVIM) imaging and diffusion kurtosis imaging (DKI) parameters for the evaluation of radiotherapy response in rabbit VX2 malignant bone tumor model.

Material and methods

Forty-seven rabbits with bone tumor were prospectively enrolled and divided into pre-treatment, considerable effect and slight effect group. Treatment response was evaluated using IVIM-DKI. IVIM-based parameters (tissue diffusion [Dt], pseudo-diffusion [Dp], perfusion fraction [fp]), and DKI-based parameters (mean diffusion coefficient [MD] and mean kurtosis [MK]) were calculated for each animal. Corresponding changes in MRI parameters before and after radiotherapy in each group were studied with one-way ANOVA. Correlations of diffusion parameters of IVIM and DKI model were computed using Pearson’s correlation test. A diagnostic model combining different diffusion parameters was established using binary logistic regression, and its ROC curve was used to evaluate its diagnostic performance for determining considerable and slight effect to malignant bone tumor.

Results

After radiotherapy, Dt and MD increased, whereas fp and MK decreased (p <  0.05). The differences in Dt, fp, MD, and MK between considerable effect and slight effect groups were statistically significant (p <  0.05). A combination of Dt, fp, and MK had the best diagnostic performance for differentiating considerable effect from slight effect (AUC = 0.913, p <  0.001).

Conclusions

A combination of IVIM- and DKI-based parameters allowed the non-invasive assessment of cellular, vascular, and microstructural changes in malignant bone tumors after radiotherapy, and holds great potential for monitoring the efficacy of tumor radiotherapy.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40644-022-00488-w.

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.

IVIM-DKI for differentiation between prostate cancer and benign prostatic hyperplasia: comparison of 1.5 T vs. 3 T MRI

OBJECTIVE

To implement an advanced spatial penalty-based reconstruction to constrain the intravoxel incoherent motion (IVIM)-diffusion kurtosis imaging (DKI) model and investigate whether it provides a suitable alternative at 1.5 T to the traditional IVIM-DKI model at 3 T for clinical characterization of prostate cancer (PCa) and benign prostatic hyperplasia (BPH).

MATERIALS AND METHODS

Thirty-two patients with biopsy-proven PCa were recruited for MRI examination (n = 16 scanned at 1.5 T, n = 16 scanned at 3 T). Diffusion-weighted imaging (DWI) with 13 b values (b = 0 to 2000 s/mm² up to 3 averages, 1.5 T: TR = 5.774 s, TE = 81 ms and 3 T: TR = 4.899 s, TE = 100 ms), T2-weighted, and T1-weighted imaging were used on the 1.5 T and 3 T MRI scanner, respectively. The IVIM-DKI signal was modeled using the traditional IVIM-DKI model and a novel model in which the total variation (TV) penalty function was combined with the traditional model to optimize non-physiological variations. Paired and unpaired t-tests were used to compare intra-scanner and scanner group differences in IVIM-DKI parameters obtained using the novel and the traditional models. Analysis of variance with post hoc test and receiver operating characteristic (ROC) curve analysis were used to assess the ability of parameters obtained using the novel model (at 1.5 T) and the traditional model (at 3 T) to characterize prostate lesions.

RESULTS

IVIM-DKI modeled using novel model with TV spatial penalty function at 1.5 T, produced parameter maps with 50-78% lower coefficient of variation (CV) than traditional model at 3 T. Novel model estimated higher D with lower D*, f and k values at both field strengths compared to traditional model. For scanner differences, the novel model at 1.5 T estimated lower D* and f values as compared to traditional model at 3 T. At 1.5 T, D and f values were significantly lower with k values significantly higher in tumor than BPH and healthy tissue. D (AUC: 0.98), f (AUC: 0.82), and k (AUC: 0.91) parameters estimated using novel model showed high diagnostic performance in cancer lesion detection at 1.5 T.

DISCUSSION

In comparison with the IVIM-DKI model at 3 T, IVIM-DKI signal modeled with the TV penalty function at 1.5 T showed lower estimation errors. The proposed novel model can be utilized for improved detection of prostate lesions.

Microstructural white matter alterations in Alzheimer's disease and amnestic mild cognitive impairment and its diagnostic value based on diffusion kurtosis imaging: a tract-based spatial statistics study

This prospective study aimed to explore the white matter microstructural alterations in Alzheimer's disease (AD) and amnestic mild cognitive impairment (aMCI) using the Tract-based Spatial Statistics (TBSS) method of diffusion kurtosis imaging (DKI).Diffusion images were collected from 45 AD patients, 42 aMCI patients, and 35 healthy controls (HC). The differences of DKI-derived parameters, including kurtosis fractional anisotropy (KFA), mean kurtosis (MK), fractional anisotropy (FA), and mean diffusivity (MD), were compared across the three groups using the TBSS method. Correlation between the altered DKI-derived parameters and the Mini-Mental State Examination (MMSE) and Montreal Cognitive Assessment (MoCA) scores were analyzed. A receiver operating characteristic curve (ROC) was used to evaluate the diagnostic performance of different white matter parameters with the strongest correlations. As a result, compared with the HC group, KFA values decreased significantly in the aMCI group. Compared with both the HC and aMCI groups, the FA, KFA, and MK values decreased significantly and the MD value increased significantly in the AD group. FA, MD, KFA, and MK values of many white matter fiber tracts were significantly correlated with MMSE and MoCA scores. The area under the ROC curve (AUC) for the splenium of corpus callosum KFA values were highest for the diagnosis of aMCI and AD patients. In conclusion, the compactness and complexity of white matter microstructures were reduced in AD and aMCI patients. DKI can provide information about the severity of AD progression, and KFA might be more sensitive for the detection of white matter microstructural alterations.

Preoperatively Grading Rectal Cancer with the Combination of Intravoxel Incoherent Motions Imaging and Diffusion Kurtosis Imaging

Purpose

To combine Intravoxel Incoherent Motions (IVIM) imaging and diffusion kurtosis imaging (DKI) which can aid in the quantification of different biological inspirations including cellularity, vascularity, and microstructural heterogeneity to preoperatively grade rectal cancer.

Methods

A total of 58 rectal patients were included into this prospective study. MRI was performed with a 3T scanner. Different combinations of IVIM-derived and DKI-derived parameters were performed to grade rectal cancer. Pearson correlation coefficients were applied to evaluate the correlations. Binary logistic regression models were established via integrating different DWI parameters for screening the most sensitive parameter. Receiver operating characteristic analysis was performed for evaluating the diagnostic performance.

Results

For individual DWI-derived parameters, all parameters except the pseudodiffusion coefficient displayed the capability of grading rectal cancer (p < 0.05). The better discrimination between high- and low-grade rectal cancer was achieved with the combination of different DWI-derived parameters. Similarly, ROC analysis suggested the combination of D (true diffusion coefficient), f (perfusion fraction), and K_app (apparent kurtosis coefficient) yielded the best diagnostic performance (AUC = 0.953, p < 0.001). According to the result of binary logistic analysis, cellularity-related D was the most sensitive predictor (odds ratio: 9.350 ± 2.239) for grading rectal cancer.

Conclusion

The combination of IVIM and DKI holds great potential in accurately grading rectal cancer as IVIM and DKI can provide the quantification of different biological inspirations including cellularity, vascularity, and microstructural heterogeneity.

An optimal acquisition and post-processing pipeline for hybrid IVIM-DKI in head and neck

Purpose

To optimize the diffusion‐weighting b values and postprocessing pipeline for hybrid intravoxel incoherent motion diffusion kurtosis imaging in the head and neck region.

Methods

Optimized diffusion‐weighting b value sets ranging between 5 and 30 b values were constructed by optimizing the Cramér‐Rao lower bound of the hybrid intravoxel incoherent motion diffusion kurtosis imaging model. With this model, the perfusion fraction, pseudodiffusion coefficient, diffusion coefficient, and kurtosis were estimated. Sixteen volunteers were scanned with a reference b value set and 3 repeats of the optimized sets, of which 1 with volunteers swallowing on purpose.

The effects of (1) b value optimization and number of b values, (2) registration type (none vs. intervolume vs. intra‐ and intervolume registration), and (3) manual swallowing artifact rejection on the parameter precision were assessed.

Results

The SD was higher in the reference set for perfusion fraction, diffusion coefficient, and kurtosis by a factor of 1.7, 1.5, and 2.3 compared to the optimized set, respectively. A smaller SD (factor 0.7) was seen in pseudodiffusion coefficient. The sets containing 15, 20, and 30 b values had comparable repeatability in all parameters, except pseudodiffusion coefficient, for which set size 30 was worse. Equal repeatability for the registration approaches was seen in all parameters of interest. Swallowing artifact rejection removed the bias when present.

Conclusion

To achieve optimal hybrid intravoxel incoherent motion diffusion kurtosis imaging in the head and neck region, b value optimization and swallowing artifact image rejection are beneficial. The optimized set of 15 b values yielded the optimal protocol efficiency, with a precision comparable to larger b value sets and a 50% reduction in scan time.

Prediction of distant metastases in patients with squamous cell carcinoma of head and neck using DWI and DCE-MRI

BACKGROUND

The primary purpose was to evaluate the prognostic potential of diffusion imaging (DWI) and dynamic contrast enhanced magnetic resonance imaging (DCE-MRI) in predicting distant metastases in squamous cell carcinoma of head and neck (HNSCC) patients. The secondary aim was to examine differences in DWI and DCE-MRI-derived parameters on the basis of human papilloma virus (HPV) status, differentiation grade, and nodal stage of HNSCC.

METHODS

Fifty-six patients underwent pretreatment DWI and DCE-MRI. Patients were divided into groups who subsequently did (n = 12) or did not develop distant metastases (n = 44). Median values of apparent diffusion coefficient (ADC), volume transfer constant (K^trans ), and mean intracellular water-lifetime (τ_i ) and volume were computed from metastatic lymph nodes and were compared between two groups. Prognostic utility of HPV status, differentiation grading, and nodal staging was also evaluated both in isolation or in combination with MRI parameters in distinguishing patients with and without distant metastases. Additionally, MRI parameters were compared between two groups based on dichotomous HPV status, differentiation grade, and nodal stage.

RESULTS

Lower but not significantly different K^trans (0.51 ± 0.15 minute^-1 vs 0.60 ± 0.05 minute^-1 ) and not significantly different τ_i (0.13 ± 0.03 second vs 0.19 ± 0.02 second) were observed in patients who developed distant metastases than those who did not. Additionally, no significant differences in ADC or volume were found. τ_i, was the best parameter in discriminating two groups with moderate sensitivity (67%) and specificity (61.4%). Multivariate logistic regression analyses did not improve the overall prognostic performance for combination of all variables. A trend toward higher τ_i was observed in HPV-positive patients than those with HPV-negative patients. Also, a trend toward higher K^trans was observed in poorly differentiated HNSCCs than those with moderately differentiated HNSCCs.

CONCLUSION

Pretreatment DCE-MRI may be useful in predicting distant metastases in HNSCC.

Perfusion-driven Intravoxel Incoherent Motion (IVIM) MRI in Oncology: Applications, Challenges, and Future Trends

Recent developments in MR hardware and software have allowed a surge of interest in intravoxel incoherent motion (IVIM) MRI in oncology. Beyond diffusion-weighted imaging (and the standard apparent diffusion coefficient mapping most commonly used clinically), IVIM provides information on tissue microcirculation without the need for contrast agents. In oncology, perfusion-driven IVIM MRI has already shown its potential for the differential diagnosis of malignant and benign tumors, as well as for detecting prognostic biomarkers and treatment monitoring. Current developments in IVIM data processing, and its use as a method of scanning patients who cannot receive contrast agents, are expected to increase further utilization. This paper reviews the current applications, challenges, and future trends of perfusion-driven IVIM in oncology.

Intravoxel incoherent motion and ADC measurements for differentiating benign from malignant thyroid nodules: utilizing the most repeatable region of interest delineation at 3.0 T

Background

There is a growing need for a reproducible and effective imaging method for the quantitative differentiation of benign from malignant thyroid nodules. This study aimed to investigate the performances of intravoxel incoherent motion (IVIM) parameters and the apparent diffusion coefficient (ADC) in differentiating malignant from benign thyroid nodules derived from the most repeatable region of interest (ROI) delineation.

Methods

Forty-three patients with 46 pathologically confirmed thyroid nodules underwent diffusion-weighted imaging (DWI) with 8 b values. Two observers measured the intravoxel incoherent motion (IVIM) parameters (D, f and D*) and the apparent diffusion coefficient (ADC), ADC₆₀₀ and ADC₉₉₀ values using whole-lesion (W-L) ROI and IVIM parameters using single-section (S-S) ROI delineation. The intraclass correlation coefficients (ICCs) and Bland-Altman plots were used to evaluate the intra- and interobserver variability. The diagnostic performance of these parameters was evaluated by generating receiver operating characteristic (ROC) curves.

Results

The ICC values of all IVIM with W-L ROI delineation were higher than those with S-S ROI delineation, and excellent intra- and interobserver reproducibility was obtained. According to the Bland-Altman plots, the 95% limits of agreement of the IVIM parameters determined by the W-L ROIs revealed smaller absolute intra- and interobserver variability than those determined by S-S ROIs. The D and ADC₆₀₀ values obtained from the W-L ROIs were the most powerful parameters in differentiating benign from the malignant nodules [area under the ROC curve = 0.962 and 0.970, P = 0.771].

Conclusions

The W-L ROI of the thyroid was considered an effective method for obtaining IVIM measurements with excellent reproducibility for differentiating benign from malignant nodules.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Recent advances in MRI of the head and neck, skull base and cranial nerves: new and evolving sequences, analyses and clinical applications

MRI is an invaluable diagnostic tool in the investigation and management of patients with pathology of the head and neck. However, numerous technical challenges exist, owing to a combination of fine anatomical detail, complex geometry (that is subject to frequent motion) and susceptibility effects from both endogenous structures and exogenous implants. Over recent years, there have been rapid developments in several aspects of head and neck imaging including higher resolution, isotropic 3D sequences, diffusion-weighted and diffusion-tensor imaging as well as permeability and perfusion imaging. These have led to improvements in anatomic, dynamic and functional imaging. Further developments using contrast-enhanced 3D FLAIR for the delineation of endolymphatic structures and black bone imaging for osseous structures are opening new diagnostic avenues. Furthermore, technical advances in compressed sensing and metal artefact reduction have the capacity to improve imaging speed and quality, respectively. This review explores novel and evolving MRI sequences that can be employed to evaluate diseases of the head and neck, including the skull base.