Next Article in Journal
Clinical Significance of COVID-19 and Diabetes: In the Pandemic Situation of SARS-CoV-2 Variants including Omicron (B.1.1.529)
Previous Article in Journal
Application of Wearable Sensors Technology for Lumbar Spine Kinematic Measurements during Daily Activities following Microdiscectomy Due to Severe Sciatica
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Feasibility of Intravoxel Incoherent Motion (IVIM) and Dynamic Contrast-Enhanced Magnetic Resonance Imaging (DCE-MRI) in Differentiation of Benign Parotid Gland Tumors

1
2nd Department of Radiology, Medical University of Gdansk, 80-214 Gdansk, Poland
2
Department of Otolaryngology, Medical University of Gdansk, 80-214 Gdansk, Poland
*
Author to whom correspondence should be addressed.
Biology 2022, 11(3), 399; https://doi.org/10.3390/biology11030399
Submission received: 2 February 2022 / Revised: 1 March 2022 / Accepted: 2 March 2022 / Published: 4 March 2022
(This article belongs to the Section Medical Biology)

Abstract

:

Simple Summary

Multiparametric magnetic resonance imaging gains recognition in the diagnostic algorithm of salivary gland tumors, providing data from different diffusion-weighted imaging (DWI) models and dynamic contrast enhanced sequences (DCE-MRI). The aim of this prospective study is to identify quantitative intravoxel incoherent motion and DCE-MRI parameters of the most frequent benign parotid tumors, pleomorphic adenomas and Warthin tumors, to compare their utility and diagnostic accuracy. With a precise pre-treatment diagnostic tool, patients can potentially avoid unnecessary diagnostic procedures and be offered optimal treatment.

Abstract

Aim: The aim of this prospective study is to identify quantitative intravoxel incoherent motion and dynamic contrast-enhanced magnetic resonance imaging parameters of the most frequent benign parotid tumors, compare their utility and diagnostic accuracy. Methods: The study group consisted of 52 patients with 64 histopathologically confirmed parotid focal lesions. Parametric maps representing apparent diffusion coefficient (ADC), pure diffusion coefficient (D), pseudo-diffusion coefficient (D*), perfusion fraction (FP) and transfer constant (Ktrans), reflux constant (Kep), extra-vascular extra-cellular volume fraction (Ve), and initial area under curve in 60 s (iAUC) have been obtained from multiparametric MRI. Results: Statistically significant (p < 0.001) inter-group differences were found between pleomorphic adenomas (PA) and Warthin tumors (WT) in all tested parameters but iAUC. Receiver operating characteristic curves were constructed to determine the optimal cut-off levels of the most significant parameters allowing differentiation between WT and PA. The Area Under the Curve (AUC) values and thresholds were for ADC: 0.931 and 1.05, D: 0.896 and 0.9, Kep: 0.964 and 1.1 and Ve: 0.939 and 0.299, respectively. Lesions presenting with a combination of ADC, D, and Ve values superior to the cut-off and Kep values inferior to the cut-off are classified as pleomorphic adenomas. Lesions presenting with combination of ADC, D, and Ve values inferior to the cut-off and Kep values superior to the cut-off are classified as Warthin tumors. Conclusions: DWI, IVIM and quantitative analysis of DCE-MRI derived parameters demonstrated distinctive features of PAs and WT and as such they seem feasible in differentiation of benign parotid gland tumors.

1. Introduction

Three to six percent of tumors within the head and neck arise in the salivary glands [1], the majority within the parotid glands. Up to 80% of tumors within the parotids are benign with the two most common types being Warthin tumors (WTs) and pleomorphic adenomas (PAs). Despite both being benign lesions, they present with different clinical course if left untreated. Malignant transformation is reported in 1.9–23.3% of untreated PAs, whereas it is less than 1% for WTs [1,2]. Thus, WTs can be surveilled if asymptomatic, while surgery is a method of treatment for pleomorphic adenomas. The aim of surgery is to obtain a complete resection with clear histopathological margins, as recurrence rate after enucleation in up to 50% [2,3]. Additionally, special attention should be given to lesions within the deep lobe as they are not only difficult to diagnose in cytology due to accessibility issues but also demand special caution during operation in order to preserve surrounding structures, especially facial nerves. Hence, a precise pre-treatment diagnosis is essential.
Reported overall sensitivity of salivary gland fine needle aspiration cytology (FNAC), which usually serves as a first-line diagnostic procedure, ranges from 86% to 100%, with specificity from 90% to 100% [4]. Sensitivity and specificity of FNAC, however, depend on technical experience of both person performing the biopsy and the evaluating pathologist, quality of cytologic preparations, morphological heterogeneity of the lesions, presence of cystic components and lesion location (with limited access when the lesion is in the deep lobe) [4].
Multiparametric magnetic resonance imaging (mpMRI) gains recognition in the diagnostic algorithm of salivary gland tumors, providing data from different diffusion-weighted imaging (DWI) models and dynamic contrast enhanced sequences (DCE-MRI) [5,6,7,8,9]. Diffusion-weighted imaging reflects the cell density of the tissue. The most used quantitative parameter of DWI, apparent diffusion coefficient (ADC), relates to the microscopic diffusion of the water molecules and is assessed in a mono-exponential model; its role seems to be well established in radiological clinical practice [5]. The intravoxel incoherent motion (IVIM) model proposed by LeBihan separates the perfusion-related parameters (pseudo-diffusion coefficient D* and perfusion fraction FP) and the diffusion-related parameters (pure diffusion coefficient D) by employing a bi-exponential analysis [10,11,12]. Previous studies proved the utility of IVIM in head and neck tumors’ differentiation [13,14,15,16,17,18], evaluating lymph nodes [19,20,21,22], monitoring outcome and/or treatment response [23,24], in the assessment of radiation-induced changes [25], and evaluating inflammatory processes, such as Sjögren’s syndrome [26,27]. Additionally, the potential role of IVIM analysis in parotid gland oncology seems to be increasing [5]. Semi-quantitative analysis of the contrast enhancement pattern in DCE-MRI reflected as time intensity curves is widely employed in salivary gland tumor differentiation [5,28]. The utility of quantitative data derived from DCE-MRI such as Ktrans (transfer constant), Kep (reflux constant), Ve (extra-vascular extra-cellular volume fraction), and iAUC (initial area under curve in 60 s) is now under ongoing research [29,30,31,32,33].
The aim of this prospective study is to identify quantitative IVIM and DCE-MRI parameters of the most frequent benign parotid tumors and compare the utility and diagnostic accuracy of those parameters. We believe that assessment of well-established features such as ADC values and time intensity curves together with novel quantitative parameters derived from IVIM and DCE-MRI will reveal tumor biology and hence provide significant information on tumor characteristics.

2. Materials and Methods

2.1. Study Group

Eighty-eight consecutive patients referred to The Department of Otolaryngology with major salivary gland nodules revealed in physical examination and/or in ultrasound underwent parametric magnetic resonance examination of the parotid glands in our institution between July 2020 and October 2021. Inclusion criteria were: 18+ years of age, the presence of at least one focal lesion within parotid glands other than of inflammatory origin/cyst and the absence of contraindications to magnetic resonance imaging/intravenous contrast agent administration. FNAC performed less than two weeks prior to MRI was an exclusion criterium; in case of one patient biopsy was performed two days after mpMRI. Patients with a history of radiotherapy in the head and neck region were not enrolled.
After the initial analysis, 31 patients were excluded from further evaluation due to the presence of inflammatory lesions (6 pts), the presence of cysts (2 pts), submandibular gland lesions (2 pts), deviation from the examination protocol (5 pts), and a lack of post-operative histopathology examination results (8 pts). In eight patients, mpMRI did not confirm the presence of a focal lesion. We have further excluded 5 patients with malignant lesions due to a small number of lesions as well as group heterogeneity (1 muco-epidermoid carcinoma, 1 adenoid cystic carcinoma, 1 acinic cell carcinoma, 1 lymph node metastasis, 1 diffuse large B-cell lymphoma).
Final study cohort consisted of fifty-two patients with 64 histopathologically confirmed focal lesions within parotid glands: 33 women (mean age 57.67, age range 21–67 years old) and 19 men (mean age 53.05, age range 24–80 years old)—Table 1. Nine patients had multiple lesions (from 2 to 3 foci). All the patients underwent surgery and histopathology examination results, which were treated as a reference standard, revealed 31 Warthin tumors, 27 pleomorphic adenomas, 3 myoepitheliomas and 3 basal cell adenomas (Table 1).

2.2. Imaging Protocol

All examinations were acquired with a 1.5 T scanner (Siemens Magnetom Sola, Erlangen, Germany) with a multi-array head-coil according to the examination protocol detailed in Table 2. Diffusion-weighted sequence was performed with 13 different b-values (0, 10, 20, 30, 50, 80, 100, 300, 400, 500, 800, 1000, 2000 s/mm2). Dynamic contrast-enhanced sequence was obtained with 36 consecutive scans with a temporal resolution of 6 s after intravenous contrast agent (CA) injection (Gadobutrol; Gadovist, Bayer Schering Pharma, Leverkusen, Germany) at a dose of 0.1 mL/kg body weight and flow rate of 3 mL/s followed by a flush of 20 mL saline solution at the same rate. CA was administered at the start of the first post-contrast scan with the use of an automatic injection system. No adverse reactions were reported.
Informed consent was obtained from all patients involved in the study. The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Independent Ethics Committee for Scientific Research (NKBBN/26/2021).

2.3. Image Analysis

Images have been independently evaluated by two radiologists with over 5 years of experience in head and neck radiology (KM—observer 1, TN—observer 2) blinded to the final histopathological results. Tumors have been identified on axial T1-weighted Turbo Spin Echo images. Initial morphological evaluation was followed by assessment of ADC values and time intensity curves (TICs). In order to assess apparent diffusion coefficient values, regions of interest (ROIs) were drawn manually on ADC maps by each observer independently over the entire tumor at an axial cross-section where the tumor maximum diameter was previously measured, carefully avoiding cystic areas and post biopsy hemorrhage (in case of two patients). ROIs were then copied within dynamic contrast-enhanced images in order to construct TICs. With concordance to the literature [34] TIC patterns were categorized as follows: type A—progressive contrast enhancement with delayed time to peak enhancement (Tpeak over 120 s) and a poor washout (less than 10% at 5 min); type B—early enhancement (Tpeak prior to 120 s) and high washout pattern (WR over 30% at 5 min); and type C—early enhancement (Tpeak prior to 120 s) and low washout pattern (WR of less than 30% at 5 min). A flat curve without enhancement has not been observed within our cohort. The range and average ROI size were as follows: observer 1: range 0.11–7.34 cm2, mean 1.7 cm2, observer 2: range 0.10–7.03 cm2, mean 1.65 cm2, respectively.

2.3.1. Intravoxel Incoherent Motion (IVIM) Analysis

IVIM analysis was performed with commercially available software syngo.via Frontier MR Body Diffusion Toolbox (Siemens, Erlangen, Germany) with an assumption of bi-exponential signal evolution model and according to implemented algorithm 3 for full non-linear fitting of extracted IVIM-related features. ROIs were manually set within the same areas as those used in previous ADC/TICs analysis into each of constructed IVIM maps. In addition, correlating ROIs were set in the contralateral normal parotid gland parenchyma for reference. The following parameters were assessed and further statistically analyzed: D (pure diffusion coefficient), D* (pseudo-diffusion coefficient) and FP (perfusion fraction).

2.3.2. Quantitative DCE-MRI Analysis

Analyses were performed with a vendor-provided software syngo.via Tissue 4D (Siemens, Erlangen, Germany). We based our analysis on the Tofts model generating pixel-based maps from TICs with constant T1 (assumed T1 of 1300 ms; value obtained in our study as well as supported by literature [34]). Arterial input function (AIF) was set to intermediate [35], providing with an optimal median chi2 value at processing 0.021 (range 0.006–0.11). For quantitative analysis, the same 2 observers independently copied and pasted ROIs to each of permeability maps, obtaining values of Ktrans (transfer constant), Kep (reflux constant), Ve (extra-vascular extra-cellular volume fraction), and iAUC (initial area under curve in 60 s). As previously, correlating ROIs were set in the contralateral parotid gland for reference.

2.3.3. Statistical Analysis

With the use of SPSS Statistica software (IBM Corporation, Armonk, NY, USA) a non-linear (non-Gaussian) distribution of variables (ADC, D, D*, Fp, Ktrans, Kep, iAUC) was confirmed. Further non-parametric tests were used for analysis of the variables. Inter-observer agreement levels were measured with intra-class correlation coefficient (ICC). ICCs were regarded as follows: values less than 0.5 indicate poor reliability, values between 0.5 and 0.75—moderate reliability, values between 0.75 and 0.9—good reliability and values greater than 0.90 indicate excellent reliability [36]. With good and excellent ICCs median values and interquartile range were calculated for each variable. Kruskal–Wallis test was used for evaluation of independent samples relations between Warthin tumors, pleomorphic adenomas and healthy parotids. Significance level was obtained at 0.05. Values of significance for multiple tests was corrected by Bonferonni method. In addition, receiver operating characteristic (ROC) curves were constructed to determine the optimal cut-off levels of the most significant parameters allowing differentiation between WT and PA. Myoepitheliomas and basal cell adenomas have been excluded from statistical analysis due to a small number of lesions.

3. Results

ICCs of tumor derived features presented with excellent inter-rater reliability except for D*—good agreement; p value < 0.001—Table 3.
Values of magnetic resonance derived features of the most common benign parotid gland tumors and of contralateral healthy parotid parenchyma are gathered in Table 4. The ADC values of WTs were significantly lower than that of PAs. Pleomorphic adenomas showed significantly higher D and Ve values than Warthin tumors (p < 0.001). D* values and Kep values of Warthin tumors were significantly higher than that of PAs (p < 0.001).
Kruskal–Wallis test for independent samples proved that statistically significant (p < 0.001) inter-group differences exist between PAs and WTs in all tested parameters but iAUC, pleomorphic adenomas differ from healthy parotid gland in all aspects while the values of D* and FP of Warthin tumors and healthy parotid glands are very similar—Table 5. The significance level was obtained at 0.05.
ROC curves were constructed to determine the optimal cut-off levels of the most significant parameters allowing differentiation between WT and PA—Figure 1. The area under the curve (AUC) values and thresholds were for ADC: 0.931 and 1.05, D: 0.896 and 0.9, Kep: 0.964 and 1.1 and Ve: 0.939 and 0.299, respectively.
Lesions presenting with combination of ADC, D, and Ve values superior to the cut-off and Kep values inferior to the cut-off are classified as pleomorphic adenomas. A representative case of PA is presented in Figure 2.
Lesions presenting with combination of ADC, D, Ve values inferior to the cut-off and Kep values superior to the cut-off are classified as Warthin tumors. Figure 3 highlights the imaging features of WTs.
AUC for the rest of analyzed features presented as follows: Ktrans 0.692, D* 0.763, FP 0.773 and iAUC 0.576. The correlation between the area under the ROC curve and diagnostic accuracy is interpreted as 0.9–1.0 excellent, 0.8–0.9 very good, 0.7–0.8 good, 0.6–0.7 sufficient, 0.5–0.6 bad, respectively; with a value less than 0.5 representing that the test is not useful. Due to bad diagnostic accuracy value of AUC for the iAUC further analysis was not performed for that selected parameter.
Lesions presenting with combination of ADC, D, and Ve values superior to cut-off and Kep values inferior to cut-off are classified as pleomorphic adenomas.
Lesions presenting with combination of ADC, D, and Ve values inferior to cut-off and Kep values superior to cut-off are classified as Warthin tumors.
Diagnostic accuracy of IVIM- and DCE-derived parameters are shown in Table 6.
When a combination of parameters (ADC, Kep and D) underwent analysis, we obtained diagnostic accuracy of 94.8% (95%CI 85.59–98.91), sensitivity of 92.59% (95%CI 75.71–99.09), specificity of 96.77% (95%CI 83.30–99,92), PPV 96.22% (95%CI 78.49–98.25) and NPV of 93.64% (95%CI 79.49–98.25).
Diagnostic accuracy increased to 96.55% (95%CI 88.09–99.58) when we analyzed a combination of ADC, Kep and D*; sensitivity—96.30% (95%CI 81.03–99.91), specificity—96.77% (95%CI 83.30–99.92), PPV 96.36% (95%CI 79.36–99.45) and NPV of 96.72% (95%CI 81.13–99.51).
However, we did not observe a further increase in diagnostic accuracy after including the fourth parameter (Ve) in the analysis. Accuracy in fact dropped to 91.33% (95%CI 80.95–97.11%) with sensitivity of 85.19% (95%CI 66.27–95.81), specificity of 96.77% (95%CI 83.30–99.92), PPV 95.90% (95%CI 77.19–99.39) and NPV of 88.05% (95%CI 74.84–94.80).

4. Discussion

Characterization of parotid focal lesions is crucial, as treatment differs accordingly to histopathological type of the tumor. Clinical findings provide limited information and FNAC is not always conclusive. Thus, pre-operative magnetic resonance imaging gains recognition in a pre-treatment diagnostic algorithm. Interpreted in conjunction with conventional morphological sequences, dynamic contrast-enhanced MRI and diffusion-weighted MRI have high potential for the characterization of parotid tumors [37], according to research comparable to that of FNAC [7,38]. Nowadays, novel advanced DWI analysis techniques are under ongoing research with the aim to identify characteristic features of different histopathological types of parotid tumors in order to further increase diagnostic accuracy [39]. In our study DWI, IVIM and quantitative analysis of DCE-MRI derived parameters demonstrated distinctive features of pleomorphic adenomas and Warthin tumors as confirmed with a Kruskal–Wallis test for independent samples. The test proved that statistically significant (p < 0.001) inter-group differences exist between pleomorphic adenomas and Warthin tumors in all tested parameters but iAUC.
Pleomorphic adenomas presented with mean ADC value of 1.32 × 10−3 mm2/s which is in concordance to previous studies reporting a cut-off ADC value between 1.267 × 10−3 mm2/s and 1.60 × 10−3 mm2/s [40,41]. The ADC values of Warthin tumors were significantly lower than that of PAs. In our study the ADC cut-off value distinguishing between WTs and PAs was established at the level of 1.05 × 10−3 mm2/s (AUC 0.931) with diagnostic accuracy of 94.8%, sensitivity 88.9%, specificity 100%, PPV 100% and NPV 91%. A similar ADC value was obtained by Nada et al. with a cut-off 1.16 ± 0.31 (SD) × 10−3 mm2/s [42].
IVIM, allowing separate evaluation of perfusion and diffusion-related features of salivary gland tumors, may help discriminate parotid gland tumors with good and excellent accuracy. Warthin tumors presented with significantly lower D values than pleomorphic adenomas with a threshold value equal or less than 0.9 × 10−3 mm2/s. Previously, Sumi et al. [43] in their prospective study including 12 pleomorphic adenomas and 8 Warthin tumors evaluated the D cut-off value at the level of 1.1 × 10−3 mm2/s allowing differentiation of WT from PAs with 100% accuracy, specificity and sensitivity. As the authors imply these results may be attributed to large areas of densely packed small lymphocytes in Warthin tumors, which greatly limit pure diffusion. D* values of WTs were significantly higher than in case of PAs. D* reflects tumor vascularity being proportional to the average blood velocity and mean capillary segment length [8].
The role of quantitative parameters of DCE-MRI reflecting the microvascular density of the lesion is less studied than IVIM. Warthin tumors showed significantly higher Kep values than pleomorphic adenomas (the median value was 2.49 min−1). This is in concordance with results obtained by Yabuuchi et al. [30] In their study on the added value of permeability imaging Kep, representing a rate constant between the extravascular space and the blood plasma, proved to be the only significant feature derived from DCE-MRI. Further the authors interpolated the results on TIC pattern, speculating that high Kep values correspond to a high wash-out rate, typically seen in case of WTs [30]. Additionally, in our study 93.5% (29 tumors out of 31 total WTs) of Warthin tumors showed a type B time intensity curve, characterized by an early enhancement (Tpeak prior to 120 s) and high washout pattern (WR over 30% at 5 min). In addition, the Ve values of WTs was significantly lower that than of pleomorphic adenomas. Taking that into consideration, we agree with Yabuuchi et al. that further analysis of Kep and Ve parameters may prove them feasible for precise evaluation of the contrast agent movement between extravascular extracellular space and blood pool [28]. Furthermore, in the light of recent research by Xu et al., Ve seems to have a discriminative potential between benign and malignant lesions [39].
The ROC analysis allowed us to propose the following thresholds allowing to differentiate WTs from PAs: ADC 1.05 (AUC 0.931), D 0.9 (AUC 0.896), Kep 1.1 (AUC 0.964) and Ve 0.299 (AUC 0.939). Lesions presenting with combination of ADC, D, and Ve values superior to the cut-off and Kep values inferior to the cut-off are classified as pleomorphic adenomas. Lesions presenting with combination of ADC, D, and Ve values inferior to the cut-off and Kep values superior to the cut-off are classified as Warthin tumors.
Similarly to Patella et al. [33], we found that a combination of analyzed parameters increases diagnostic accuracy and its indices, however the optimal combination included no more than three features (with the best values obtained with ADC, Kep and D or ADC, Kep and D*).
The authors are aware of the following limitations. Firstly, it was a single-center study. Secondly, the sample size calculation has not been conducted as we recruited all consecutive patients within certain time period about inclusion and exclusion criteria. Additionally, due to a small number of lesions as well as group heterogeneity malignant lesions needed to be excluded from the study. However, we are now collecting cases for a validation study and researching parotid malignancies. At this point we have not yet analyzed different histopathological subgroups of PAs and WTs; however, as this study is a part of an ongoing prospective research, we will hopefully be able to investigate if and at to what extend the analyzed features may be dependent on pathological subtype of the tumor.

5. Conclusions

In our study IVIM and quantitative DCE-MRI derived parameters demonstrated distinctive features of PAs and WTs and as such they seem feasible in differentiation of benign parotid gland tumors in a multiparametric magnetic resonance approach. With a precise pre-treatment diagnostic tool, patients can be offered optimal treatment.

Author Contributions

Conceptualization, K.M. and B.M.; methodology, K.M.; software, A.G.; validation, K.M., B.M., T.N. and E.S.; formal analysis, K.M.; investigation, K.M., T.N. and A.G.; resources, K.M.; data curation, K.M.; writing—original draft preparation K.M. and A.G.; writing—review and editing, K.M., B.M. and E.S.; visualization, K.M.; supervision, B.M. and E.S.; project administration, K.M.; funding acquisition, K.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received part funding from an internal grant acknowleged by Medical University of Gdansk, Poland to K.M. The grant “Young Creator of Science” (number 664/323/61/71-1208) is conducted within the “Excellence Initiative—Research University” Programme.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Independent Ethics Committee for Scientific Research of Medical University of Gdansk, Poland. Protocol code: NKBBN/26/2021. Date of approval 17 February 2021.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Data available on request due to restrictions (privacy). The data presented in this study are available on request from the corresponding author. The data are not publicly available due to above mentioned restrictions.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Barnes, L.; Eveson, J.W.; Reichart, P.; Sidransky, D. World Health Organization Classification of Tumors—Pathology and Genetics of Head and Neck Tumors; IARC Press: Lyon, France, 2005; pp. 209–281. [Google Scholar]
  2. Witt, R.L.; Eisele, D.W.; Morton, R.P.; Nicolai, P.; Poorten, V.V.; Zbären, P. Etiology and management of recurrent parotid pleomorphic adenoma. Laryngoscope 2015, 125, 888–893. [Google Scholar] [CrossRef] [PubMed]
  3. Riad, M.A.; Abdel-Rahman, H.; Ezzat, W.F.; Adly, A.; Dessouky, O.; Shehata, M. Variables related to recurrence of pleomorphic adenomas: Outcome of parotid surgery in 182 cases. Laryngoscope 2011, 121, 1467–1472. [Google Scholar] [CrossRef]
  4. Faquin, W.C.; Rossi, E.D.; Baloch, Z.; Barkan, G.A.; Foschini, M.P.; Kurtycz, D.F.I.; Pusztaszeri, M.; Vielh, P. The Milan System for Reporting Salivary Gland Cytopathology; Springer International Publishing AG: Cham, Switzerland, 2018; pp. 1–10. [Google Scholar]
  5. Coudert, H.; Mirafzal, S.; Dissard, A.; Boyer, L.; Montoriol, P.F. Multiparametric magnetic resonance imaging of parotid tumors: A systematic review. Diagn. Interv. Imaging 2021, 102, 121–130. [Google Scholar] [CrossRef] [PubMed]
  6. Noij, D.P.; Martens, R.M.; Marcus, J.T.; de Bree, R.; Leemans, C.R.; Castelijns, J.A.; de Jong, M.C.; de Graaf, P. Intravoxel incoherent motion magnetic resonance imaging in head and neck cancer: A systematic review of the diagnostic and prognostic value. Oral Oncol. 2017, 68, 81–91. [Google Scholar] [CrossRef] [PubMed]
  7. Liang, Y.Y.; Xu, F.; Guo, Y.; Wang, J. Diagnostic accuracy of magnetic resonance imaging techniques for parotid tumors, a systematic review and meta-analysis. Clin. Imaging 2018, 52, 36–43. [Google Scholar] [CrossRef] [PubMed]
  8. Mikaszewski, B.; Markiet, K.; Smugała, A.; Stodulski, D.; Szurowska, E.; Stankiewicz, C. Diffusion-weighted MRI in the differential diagnosis of parotid malignancies and pleomorphic adenomas: Can the accuracy of dynamic MRI be enhanced? Oral Surg. Oral Med. Oral Pathol. Oral Radiol. 2017, 124, 95–103. [Google Scholar] [CrossRef]
  9. Mikaszewski, B.; Markiet, K.; Smugała, A.; Stodulski, D.; Szurowska, E.; Stankiewicz, C. Clinical and demographic data improve diagnostic accuracy of dynamic contrast-enhanced and diffusion-weighted MRI in differential diagnostics of parotid gland tumors. Oral Oncol. 2020, 111, 104932. [Google Scholar] [CrossRef]
  10. Le Bihan, D.; Breton, E.; Lallemand, D.; Aubin, M.L.; Vignaud, J.; Laval-Jeantet, M. Separation of diffusion and perfusion in intravoxel incoherent motion MR imaging. Radiology 1988, 168, 497–505. [Google Scholar] [CrossRef]
  11. Le Bihan, D.; Turner, R. The capillary network: A link between IVIM and classical perfusion. Magn. Reason. Med. 1992, 27, 171–178. [Google Scholar] [CrossRef]
  12. Le Bihan, D. What can we see with IVIM MRI? Neuroimage 2019, 187, 56–67. [Google Scholar] [CrossRef]
  13. Ai, Q.Y.; King, A.D.; Chan, J.S.M.; Chen, W.; Chan, K.C.A.; Woo, J.K.S.; Zee, B.C.Y.; Chan, A.T.C.; Poon, D.M.C.; Ma, B.B.Y.; et al. Distinguishing early-stage nasopharyngeal carcinoma from benign hyperplasia using intravoxel incoherent motion diffusion-weighted MRI. Eur. Radiol. 2019, 29, 5627–5634. [Google Scholar] [CrossRef] [PubMed]
  14. Martínez Barbero, J.P.; Rodríquez Jiménez, I.; Martin Noguerol, T.; Luna Alcalá, A. Utility of MRI diffusion techniques in the evaluation of tumors of the head and neck. Cancers 2013, 5, 875–889. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Fujima, N.; Yoshida, D.; Sakashita, T.; Homma, A.; Tsukahara, A.; Tha, K.K.; Kudo, K.; Shirato, H. Intravoxel incoherent motion diffusion-weighted imaging in head and neck squamous cell carcinoma: Assessment of perfusion-related parameters compared to dynamic contrast-enhanced MRI. Magn. Reason. Imaging 2014, 32, 1206–1213. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Marzi, S.; Piludu, F.; Vidiri, A. Assessment of diffusion parameters by intravoxel incoherent motion MRI in head and neck squamous cell carcinoma. NMR Biomed. 2013, 26, 1806–1814. [Google Scholar] [CrossRef] [PubMed]
  17. Fujima, N.; Yoshida, D.; Sakashita, T.; Homma, A.; Tsukahara, A.; Shimizu, Y.; Tha, K.K.; Kudo, K.; Shirato, H. Prediction of the treatment outcome using intravoxel incoherent motion and diffusional kurtosis imaging in nasal or sinonasal squamous cell carcinoma patients. Eur. Radiol. 2017, 27, 956–965. [Google Scholar] [CrossRef] [Green Version]
  18. Ma, G.; Xu, X.Q.; Zhu, L.N.; Jiang, J.S.; Su, G.Y.; Hu, H.; Bu, S.S.; Wu, F.Y. Intravoxel Incoherent Motion Magnetic Resonance Imaging for Assessing Parotid Gland Tumors: Correlation and Comparison with Arterial Spin Labeling Imaging. Korean J. Radiol. 2021, 22, 243–252. [Google Scholar] [CrossRef]
  19. Hauser, T.; Essig, M.; Jensen, A.; Laun, F.B.; Münter, M.; Maier-Hein, K.H.; Stieltjes, B. Prediction of treatment response in head and neck carcinomas using IVIM-DWI: Evaluation of lymph node metastasis. Eur. J. Radiol. 2014, 83, 783–787. [Google Scholar] [CrossRef]
  20. Hejduk, B.; Bobek-Billewicz, B.; Rutkowski, T.; Hebda, A.; Zawadzka, A.; Jurkowski, M.K. Application of Intravoxel Incoherent Motion (IVIM) Model for Differentiation Between Metastatic and Non-Metastatic Head and Neck Lymph Nodes. Pol. J. Radiol. 2017, 82, 506–510. [Google Scholar] [CrossRef] [Green Version]
  21. Marzi, S.; Piludu, F.; Forina, C.; Sanguineti, G.; Covello, R.; Spriano, G.; Vidiri, A. Correlation study between intravoxel incoherent motion MRI and dynamic contrast-enhanced MRI in head and neck squamous cell carcinoma: Evaluation in primary tumors and metastatic nodes. Magn. Reason. Imaging 2017, 37, 1–8. [Google Scholar] [CrossRef]
  22. Liang, L.; Luo, X.; Lian, Z.; Chen, W.; Zhang, B.; Dong, Y.; Liang, C.; Zhang, S. Lymph node metastasis in head and neck squamous carcinoma: Efficacy of intravoxel incoherent motion magnetic resonance imaging for the differential diagnosis. Eur. J. Radiol. 2017, 90, 159–165. [Google Scholar] [CrossRef]
  23. Paudyal, R.; Oh, J.H.; Riaz, N.; Venigalla, P.; Li, J.; Hatzoglou, V.; Leeman, J.; Nunez, D.A.; Lu, Y.; Deasy, J.O.; et al. Intravoxel incoherent motion diffusion-weighted MRI during chemoradiation therapy to characterize and monitor treatment response in human papillomavirus head and neck squamous cell carcinoma. J. Magn. Reason. Imaging 2017, 45, 1013–1023. [Google Scholar] [CrossRef]
  24. Ding, Y.; Hazle, J.D.; Mohamed, A.S.; Frank, S.J.; Hobbs, B.P.; Colen, R.R.; Gunn, G.B.; Wang, J.; Kalpathy-Cramer, J.; Garden, A.S.; et al. Intravoxel incoherent motion imaging kinetics during chemoradiotherapy for human papillomavirus-associated squamous cell carcinoma of the oropharynx: Preliminary results from a prospective pilot study. NMR Biomed. 2015, 28, 1645–1654. [Google Scholar] [CrossRef] [PubMed]
  25. Marzi, S.; Farneti, A.; Vidiri, A.; Di Giuliano, F.; Marucci, L.; Spasiano, F.; Terrenato, I.; Sanguineti, G. Radiation-induced parotid changes in oropharyngeal cancer patients: The role of early functional imaging and patient-/treatment-related factors. Radiat Oncol. 2018, 13, 189. [Google Scholar] [CrossRef] [PubMed]
  26. Chu, C.; Zhou, N.; Zhang, H.; Dou, X.; Li, M.; Liu, S.; Zhu, Y.; Chen, W.; Chan, Q.; He, J.; et al. Correlation between intravoxel incoherent motion MR parameters and MR nodular grade of parotid glands in patients with Sjögren’s syndrome: A pilot study. Eur. J. Radiol. 2017, 86, 241–247. [Google Scholar] [CrossRef] [PubMed]
  27. Su, G.Y.; Xu, X.Q.; Wang, Y.Y.; Hu, H.; Shen, J.; Hong, X.N.; Shi, H.B.; Wu, F.Y. Feasibility study of using intravoxel incoherent motion mri to detect parotid gland abnormalities in early-stage Sjögren syndrome patients. J. Magn. Reason. Imaging 2016, 43, 1455–1461. [Google Scholar] [CrossRef] [PubMed]
  28. Yabuuchi, H.; Fukuya, T.; Tajima, T.; Hachitanda, Y.; Tomita, K.; Koga, M. Salivary gland tumors: Diagnostic value of gadolinium-enhanced dynamic MR imaging with histopathologic correlation. Radiology 2003, 226, 345–354, Erratum in Radiology 2003, 227, 909. [Google Scholar] [CrossRef]
  29. Sumi, M.; Nakamura, T. Head and neck tumors: Assessment of perfusion-related parameters and diffusion coefficients based on the intravoxel incoherent motion model. AJNR Am. J. Neuroradiol. 2013, 34, 410–416. [Google Scholar] [CrossRef] [Green Version]
  30. Yabuuchi, H.; Kamitani, T.; Sagiyama, K.; Yamasaki, Y.; Hida, T.; Matsuura, Y.; Hino, T.; Murayama, Y.; Yasumatsu, R.; Yamamoto, H. Characterization of parotid gland tumors: Added value of permeability MR imaging to DWI and DCE-MRI. Eur. Radiol. 2020, 30, 6402–6412. [Google Scholar] [CrossRef]
  31. Xu, X.Q.; Choi, Y.J.; Sung, Y.S.; Yoon, R.G.; Jang, S.W.; Park, J.E.; Heo, Y.J.; Baek, J.H.; Lee, J.H. Intravoxel Incoherent Motion MR Imaging in the Head and Neck: Correlation with Dynamic Contrast-Enhanced MR Imaging and Diffusion-Weighted Imaging. Korean J. Radiol. 2016, 17, 641–649. [Google Scholar] [CrossRef] [Green Version]
  32. Sumi, M.; Nakamura, T. Head and neck tumours: Combined MRI assessment based on IVIM and TIC analyses for the differentiation of tumors of different histological types. Eur. Radiol. 2014, 24, 223–231. [Google Scholar] [CrossRef]
  33. Patella, F.; Franceschelli, G.; Petrillo, M.; Sansone, M.; Fusco, R.; Pesapane, F.; Pompili, G.; Ierardi, A.M.; Saibene, A.M.; Moneghini, L.; et al. A multiparametric analysis combining DCE-MRI- and IVIM -derived parameters to improve differentiation of parotid tumors: A pilot study. Future Oncol. 2018, 14, 2893–2903. [Google Scholar] [CrossRef] [PubMed]
  34. ITIS Foundation. Available online: https://itis.swiss/virtual-population/tissue-properties/database/relaxation-times/ (accessed on 1 November 2021).
  35. Parker, G.J.; Roberts, C.; Macdonald, A.; Buonaccorsi, G.A.; Cheung, S.; Buckley, D.L.; Jackson, A.; Watson, Y.; Davies, K.; Jayson, G.C. Experimentally-derived functional form for a population-averaged high-temporal-resolution arterial input function for dynamic contrast-enhanced MRI. Magn. Reason. Med. 2006, 56, 993–1000. [Google Scholar] [CrossRef] [PubMed]
  36. Koo, T.K.; Li, M.Y. A Guideline of Selecting and Reporting Intraclass Correlation Coefficients for Reliability Research. J. Chiropr. Med. 2017, 15, 155–163, Erratum in J. Chiropr. Med. 2017, 16, 346. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Friedman, E.; Patino, M.O.; Abdel Razek, A.A.K. MR Imaging of Salivary Gland Tumors. Magn. Reson. Imaging Clin. N. Am. 2022, 30, 135–149. [Google Scholar] [CrossRef]
  38. Mikaszewski, B.; Markiet, K.; Smugała, A.; Stodulski, D.; Szurowska, E.; Stankiewicz, C. An algorithm for preoperative differential diagnostics of parotid tumours on the basis of their dynamic and diffusion-weighted magnetic resonance images: A retrospective analysis of 158 cases. Folia Morphol. 2018, 77, 29–35. [Google Scholar] [CrossRef] [Green Version]
  39. Xu, Z.; Chen, M.; Zheng, S.; Chen, S.; Xiao, J.; Hu, Z.; Lu, L.; Yang, Z.; Lin, D. Differential diagnosis of parotid gland tumours: Application of SWI combined with DWI and DCE-MRI. Eur. J. Radiol. 2022, 146, 110094. [Google Scholar] [CrossRef]
  40. Matsusue, E.; Fujihara, Y.; Matsuda, E.; Tokuyasu, Y.; Nakamoto, S.; Nakamura, K.; Ogawa, T. Differentiating parotid tumors by quantitative signal intensity evaluation on MR imaging. Clin. Imaging 2017, 46, 37–43. [Google Scholar] [CrossRef]
  41. Elmokadem, A.H.; Abdel Khalek, A.M.; Abdel Wahab, R.M.; Tharwat, N.; Gaballa, G.M.; Elata, M.A.; Amer, T. Diagnostic Accuracy of Multiparametric Magnetic Resonance Imaging for Differentiation Between Parotid Neoplasms. Can. Assoc. Radiol. J. 2019, 70, 264–272. [Google Scholar] [CrossRef] [Green Version]
  42. Nada, A.A.; Youssef, A.A.; Basmy, A.E.; Amin, A.W.; Shokry, A.M. Diffusion-weighted imaging of the parotid gland: Can the apparent diffusion coefficient discrim- inate between normal and abnormal parotid gland tissues? Erciyes Med. J. 2017, 39, 125–130. [Google Scholar] [CrossRef] [Green Version]
  43. Sumi, M.; Van Cauteren, M.; Sumi, T.; Obara, M.; Ichikawa, Y.; Nakamura, T. Salivary gland tumors: Use of intravoxel incoherent motion MR imaging for assessment of diffusion and perfusion for the differentiation of benign from malignant tumors. Radiology 2012, 263, 770–777. [Google Scholar] [CrossRef]
Figure 1. ROC of the most significant parameters allowing differentiation between WT and PA. AUC values and thresholds were for (A) ADC. AUC 0.931, threshold 1.05. (B) Pure diffusion coefficient (D) AUC 0.896, threshold 0.9. (C) Kep. AUC 0.964, threshold 1.1. (D) Ve. AUC 0.939, threshold 0.299.
Figure 1. ROC of the most significant parameters allowing differentiation between WT and PA. AUC values and thresholds were for (A) ADC. AUC 0.931, threshold 1.05. (B) Pure diffusion coefficient (D) AUC 0.896, threshold 0.9. (C) Kep. AUC 0.964, threshold 1.1. (D) Ve. AUC 0.939, threshold 0.299.
Biology 11 00399 g001
Figure 2. Pleomorphic adenoma (arrow) in superficial lobe of the right parotid. (A) ADC map. Lesion presents with high signal intensity correlating with lack of diffusion restriction. (B) Time Intensity Curve showing continuous contrast agent uptake consistant with type A curve and typical for adenomas. (C) IVIM. D map (high values). (D) Permeability. Kep map (low values).
Figure 2. Pleomorphic adenoma (arrow) in superficial lobe of the right parotid. (A) ADC map. Lesion presents with high signal intensity correlating with lack of diffusion restriction. (B) Time Intensity Curve showing continuous contrast agent uptake consistant with type A curve and typical for adenomas. (C) IVIM. D map (high values). (D) Permeability. Kep map (low values).
Biology 11 00399 g002
Figure 3. Warthin tumor (arrow) in superficial lobe of the left parotid. (A) ADC map. Lesion presents with low signal intensity correlating with diffusion restriction. (B) Time intensity curve showing rapid contrast agent uptake with high-rate washout corresponding to type B curve typical for WT. (C) IVIM. D map (low values). (D) Permeability. Kep map (high values).
Figure 3. Warthin tumor (arrow) in superficial lobe of the left parotid. (A) ADC map. Lesion presents with low signal intensity correlating with diffusion restriction. (B) Time intensity curve showing rapid contrast agent uptake with high-rate washout corresponding to type B curve typical for WT. (C) IVIM. D map (low values). (D) Permeability. Kep map (high values).
Biology 11 00399 g003
Table 1. Final study cohort demographic data.
Table 1. Final study cohort demographic data.
Warthin TumorsPleomorphic AdenomasMyoepitheliomasBasal Cell Adenomas
SexFemales171932
Males14801
Total number of lesions312733
Mean BMIFemales29.226.325.324.9
Males32.423.7-26.5
SideRight1518-1
Left16932
LocationSuperficial lobe2515-2
Deep lobe6103-
Accessory lobe-2-1
Size * (mm)Females1717388
Males2015-19
LymphadenopathyFemales6111
Males72-1
Facial nerve palsyFemales0000
Males01 (post-operative)00
* Mean biggest diameter in axial plane.
Table 2. Imaging protocol.
Table 2. Imaging protocol.
SequenceOrientationTR (ms)TE (ms)SL/SPCFOVMatrix
T2 BladeSagittal40701443.0/3.6250∗250320∗320
T2 TIRMCoronal60001003.0/3.6229∗271540∗640
T1 TSETransverse555113.0/3.6195∗240416∗512
T2 TSE DIXONTransverse3050753.5/4.0218∗240442∗608
ep2d_DWITransverse8130654.0236∗250208∗220
ep2d_DWI ADCTransverse8130654.0236∗250208∗220
T1 Vibe DIXONTransverse7.72.41.0i189∗220163∗224
T1 Vibe CM Dyn 1–36Transverse4.21.63.5i260∗260154∗192
T1 Vibe DIXON CMTransverse7.72.41.0i189∗220163∗224
TIRM: turbo inversion recovery magnitude; Vibe: volumetric interpolated breath-hold examination; TSE: Turbo Spin Echo; DWI: diffusion-weighted imaging; ADC: apparent diffusion coefficient; CM: contrast medium; TR: repetition time; TE: echo time; SL: slice thickness; FOV: field of view; i: isotropic.
Table 3. Results of the ICC analysis.
Table 3. Results of the ICC analysis.
IVIMPermeability
DD*FPROIKtransKepVeiAUCROI
ICC0.9630.8270.9250.9970.9630.9900.9400.9670.995
95% CI0.957–0.9840.715–0.8950.876–0.9550.995–0.9980.940–0.9780.984–0.9940.901–0.9630.946–0.9800.991–0.998
p value < 0.001. ICC: intra-class correlation coefficient; CI: Confidence Interval; D: pure diffusion coefficient; D*: pseudo-diffusion coefficient; FP: fraction perfusion; Ktrans: transfer constant; Kep: reflux constant; Ve: extra-vascular extra-cellular volume fraction; iAUC: initial area under curve in 60 s.
Table 4. Magnetic resonance features of common benign parotid gland tumors.
Table 4. Magnetic resonance features of common benign parotid gland tumors.
Warthin TumorsPleomorphic AdenomasContralateral Parotid
Magnetic Resonance ParametersMedianInterquartile RangeMedianInterquartile RangeMedianInterquartile Range
ADC 10−3 mm2/s0.680.11.320.250.840.21
IVIMD 10−3 mm2/s0.630.171.190.530.710.12
D*10−3 mm2/s135.9145.547.0520.7117.7641.89
FP %24.88.534.414.523.68.3
PermeabilityKtrans min−10.410.290.170.340.130.14
Kep min−12.491.040.310.811.220.49
Ve0.160.10.620.760.120.09
iAUC0.270.110.240.380.140.15
No of lesionsNo of lesions
TICType A019
Type B291
Type C27
ADC: Apparent Diffusion Coefficient; IVIM: Incoherent Voxel Motion; D: pure diffusion coefficient; D*: pseudo-diffusion coefficient; FP: fraction perfusion; Ktrans: transfer constant; Kep: reflux constant; Ve: extra-vascular extra-cellular volume fraction; iAUC: initial area under curve in 60 s; TIC: Time Intensity Curve.
Table 5. Kruskal–Wallis test for independent samples.
Table 5. Kruskal–Wallis test for independent samples.
Sample 1/Sample 2 Pleomorphic Adenomas vs. Warthin TumorsPleomorphic Adenomas vs. Healthy ParotidsWarthin Tumors vs. Healthy Parotids
ADC<0.001<0.001<0.001
IVIMD<0.001<0.0010.049
D*<0.001<0.0010.922
FP<0.001<0.0010.52
PermeabilityKtrans0.0040.042<0.001
Kep<0.001<0.001<0.001
Ve<0.001<0.0010.004
iAUC0.319<0.001<0.001
ADC: apparent diffusion coefficient; IVIM: incoherent voxel motion; D: pure diffusion coefficient; D*: pseudo-diffusion coefficient; FP: fraction perfusion; Ktrans: transfer constant; Kep: reflux constant; Ve: extra-vascular extra-cellular volume fraction; iAUC: initial area under curve in 60 s.
Table 6. Diagnostic accuracy of IVIM- and DCE-derived parameters.
Table 6. Diagnostic accuracy of IVIM- and DCE-derived parameters.
Sensitivity %95% CISpecificity %95% CIPPV %95% CINPV %95% CIAccuracy %95% CI
ADC88.979.8–97.710088.8–100100NA91.077.7–96.794.885.6–98.9
IVIMD81.561.9–93.787.170.2–96.484.968.8–93.484.170.4–92.284.572.6–92.6
D*85.266.3–95.870.952.0–85.872.2459.5–82.284.468.0–93.277.764.8–87.5
FP48.228.7–68.193.678.6–99.286.962.1–96.467.158.3–74.872.258.9–83.2
PermeabilityKtrans44.425.5–64.710088.8–100100NA66.959.2–73.973.960.7–84.5
Kep88.970.8–97.796.883.3–99.996.177.9–99.490.177.1–96.693.183.2–98.1
Ve85.266.3–95.810088.8–100100NA88.475.5–94.993.183.2–98.1
iAUCNANANANANANANANANANA
ADC: apparent diffusion coefficient; IVIM: incoherent voxel motion; CI: confidence interval; D: pure diffusion coefficient; D*: pseudo-diffusion coefficient; FP: fraction perfusion; Ktrans: transfer constant; Kep: reflux constant; Ve: extra-vascular extra-cellular volume fraction; iAUC: initial area under curve in 60 s.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Markiet, K.; Glinska, A.; Nowicki, T.; Szurowska, E.; Mikaszewski, B. Feasibility of Intravoxel Incoherent Motion (IVIM) and Dynamic Contrast-Enhanced Magnetic Resonance Imaging (DCE-MRI) in Differentiation of Benign Parotid Gland Tumors. Biology 2022, 11, 399. https://doi.org/10.3390/biology11030399

AMA Style

Markiet K, Glinska A, Nowicki T, Szurowska E, Mikaszewski B. Feasibility of Intravoxel Incoherent Motion (IVIM) and Dynamic Contrast-Enhanced Magnetic Resonance Imaging (DCE-MRI) in Differentiation of Benign Parotid Gland Tumors. Biology. 2022; 11(3):399. https://doi.org/10.3390/biology11030399

Chicago/Turabian Style

Markiet, Karolina, Anna Glinska, Tomasz Nowicki, Edyta Szurowska, and Boguslaw Mikaszewski. 2022. "Feasibility of Intravoxel Incoherent Motion (IVIM) and Dynamic Contrast-Enhanced Magnetic Resonance Imaging (DCE-MRI) in Differentiation of Benign Parotid Gland Tumors" Biology 11, no. 3: 399. https://doi.org/10.3390/biology11030399

APA Style

Markiet, K., Glinska, A., Nowicki, T., Szurowska, E., & Mikaszewski, B. (2022). Feasibility of Intravoxel Incoherent Motion (IVIM) and Dynamic Contrast-Enhanced Magnetic Resonance Imaging (DCE-MRI) in Differentiation of Benign Parotid Gland Tumors. Biology, 11(3), 399. https://doi.org/10.3390/biology11030399

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop