Erkan G?k?e,Murat Beyhan

Erkan G?k?e,Murat Beyhan,Department of Radiology,Tokat Gaziosmanpasa University,Faculty of Medicine,Tokat 60100,Turkey

Abstract BACKGROUND Increased use of functional magnetic resonance imaging(MRI)methods such as diffusion-weighted imaging(DWI)and dynamic contrast-enhanced(DCE)MRI consisting of sequential contrast series,allows us to obtain more information on the microstructure,cellularity,interstitial distance,and vascularity of tumors,which has increased the discrimination power for benign and malignant salivary gland tumors(SGTs).In the last few years,quantitative DCE MRI data containing T1 perfusion parameters(Ktrans,Kep and Ve),were reported to contribute to the differentiation of benign or malignant subtypes in SGTs.AIM To evaluate the diagnostic efficacy of DWI and semiquantitative and quantitative perfusion MRI parameters in SGTs.METHODS Diffusion MRI[apparent diffusion coefficient(ADC)value]with a 1.5 T MR machine,semiquantitative perfusion MRI[time intensity curve(TIC)pattern],and quantitative perfusion MRI examinations(Ktrans,Kep and Ve)of 73 tumors in 67 patients with histopathological diagnosis performed from 2017 to 2021 were retrospectively evaluated.In the ADC value and semiquantitative perfusion MRI measurements,cystic components of the tumors were not considered,and the region of interest(ROI)was manually placed through the widest axial section of the tumor.TIC patterns were divided into four groups:Type A = Tpeak >120 s;type B = Tpeak ≤ 120 s,washout ratio(WR)≥ 30%;type C = Tpeak ≤ 120 s,WR <30%;and type D = flat TIC.For the quantitative perfusion MRI analysis,a 3D ROI was placed in the largest solid component of the tumor,and the Ktrans,Kep and Ve values were automatically generated.RESULTS The majority of SGTs were located in the parotid glands(86.3%).Of all the SGTs,68.5% were benign and 31.5% were malignant.Significant differences were found for ADC values among pleomorphic adenomas(PMAs),Warthin's tumors(WTs),and malignant tumors(MTs)(P <0.001).PMAs had type A and WTs had type B TIC pattern while the vast majority of MTs and other benign tumors(OBTs)(54.5% and 45.5%,respectively)displayed type C TIC pattern.PMAs showed no washout,while the highest mean WR was observed in WTs(59% ± 11%).Ktrans values of PMAs,WTs,OBTs,and MTs were not significantly different.Kep values of PMAs and WTs were significantly different from those of OBTs and MTs.Mean Ve value of WTs was significantly different from those of PMAs,OBTs,and MTs(P <0.001).CONCLUSION The use of quantitative DCE parameters along with diffusion MRI and semiquantitative contrastenhanced MRI in SGTs could improve the diagnostic accuracy.

Key Words:Diffusion-weighted imaging;Dynamic contrast-enhanced imaging;Magnetic resonance imaging;Perfusion imaging;Salivary gland tumor;Tumor

INTRODUCTION

Salivary gland tumors(SGTs)account for about 2.0%-6.5% of all head and neck tumors.Approximately 70% of them originate from the parotid glands,and a small number have submandibular,sublingual,and minor salivary gland origins.While the majority of tumors from the parotid glands are benign,malignancies are more common in those located in other glands.Preoperative characterization of SGTs is important for treatment planning.The choice of surgery method for SGTs is closely associated with the histology of the tumor.Diagnosis is mostly based on combined evaluation of clinical features and findings from physical examinations,imaging and cytological observations.Fine-needle aspiration biopsy(FNAB)is the most commonly used method for cytological examinations but complex pathologies can result in false positives and false negatives in malignant tumors(MTs)[1,2].Conventional magnetic resonance imaging(MRI)is very useful for identifying the tumor location,morphology,extension,and its association with the nerves and inner structure[1-3].However,diagnosing MTs and benign tumors(BTs)by conventional MRI can be difficult due to overlapping findings[1,2,4].In recent years,an increase has been reported in diagnostic accuracy in SGTs for distinguishing between MTs and BTs with the use of diffusion-weighted imaging(DWI)and dynamic contrast-enhanced(DCE)MRI techniques[1,5-9].DCE MRI is used to track an exogenous,paramagnetic contrast agent in tissues and has been a powerful tool in the characterization of tumor hemodynamics[1,3,10,11].As a semiquantitative method in DCE MRI,patterns have been established by measuring time to peak(Tpeak)and washout ratio(WR)on the time intensity curve(TIC)[1,3].Tpeakis closely related to microvessel count while WR reflects the stromal cellularity grade.On quantitative DCE MRI,on the other hand,perfusion parameters such as Ktrans[volume transfer constant between blood plasma and the extracellular extravascular space(EES)],Kep(flux rate constant between the EES and plasma),and Ve(EES fractional volume)are used[1,3].Although there are many studies dealing with diffusion and semiquantitative DCE MRIs in SGTs,the number of quantitative MRI studies is limited[12-15].In the present study,the diagnostic value of diffusion MRI and semiquantitative and quantitative perfusion MRI parameters was evaluated in SGTs.

MATERIALS AND METHODS

Patients

The study was conducted retrospectively following approval by the local ethics committee(20-KAEK-105).A total of 67 patients with tumors originating from or involving the salivary glands were included.The study included patients who had swelling of the face or in the salivary glands,who were subjected to MRI,diffusion MRI,and perfusion MRI examinations at our hospital between April 2017 and February 2021 and who were diagnosed histopathologically after FNAB,Tru-cut biopsy,or surgical removal.For this study,patients whose neck and maxillofacial MRI examination reports included the description of a mass in the salivary glands were surveyed in picture archiving and communication systems.A total of 33 patients were excluded:2 patients with intra-lesion hemorrhage due to FNAB before the MRI examination,16 patients who had contrast-enhanced MRI but did not have perfusion MRI series,and 15 patients whose diagnosis was not confirmed histopathologically.Thus,a total of 73 MTs and BTs,which originated from major and minor salivary glands in 67 patients,were included in the study(Figure 1).

MRI scanning and measurements

MRI was performed on a 1.5 T superconducting MRI system[General Electric(GE)Signa Explorer Software Version 25;GE Healthcare,Milwaukee,WI,United States,2016]with head and neck array coils.Routine MRI sequences included axial T1-weighted[time to repetition/time to echo(TR/TE),456 ms/8.1 ms],in phase axial T2-weighted(TR/TE,3711 ms/82.8 ms),sagittal T2-weighted(TR/TE,4499 ms/88.2 ms),and coronal T2-weighted(TR/TE,4380 ms/84.6 ms).DCE MRI was performed with a T1-weighted 3D fast spoiled gradient echo sequence[TR/TE/time to inversion,3.8 ms/1.3 ms/15 ms;flip angle,20°].The contrast agent Gd-DTPA(Dotarem,Guerbet,France)was injected after the fourth dynamic sequence acquisition at a rate of 2.0 mL/sviathe right antecubital vein.The contrast agent was administered at a dose of 0.2 mmoL/kg body weight.Immediately after the injection of the contrast agent,a 20 mL saline flush was carried out at the same injection rate.In total,18-21 dynamic sequence acquisitions with 30 dynamic images per sequence were performed with total scanning time ranging from 3 min and 11 s to 5 min and 24 s.The location,morphology,and internal structure of the tumor were evaluated by conventional MRI(Figures 2A and B,3A and B,4A and B,5A and B).

DWI was performed using a multislice echo-planar single-shot spin-echo sequence,on the axial plane(TR/TE = 5476 ms/95.4 ms,field of view = 26 cm,matrix = 96 × 128,section thickness = 4-5 mm,and interslice gap = 4 mm).Three diffusion gradients were applied sequentially in the x,y,and z directions withbvalues of 0 and 1000 s/mm2(Figure 2B).The acquisition time varied from 60 s to 120 s.The apparent diffusion coefficient(ADC)maps were generated automatically.

"GE Advantage Windows Workstation 4.7" was used to determine ADC values on diffusion MRI and to perform measurements in semiquantitative and quantitative perfusion MRI.Image analysis and region of interest(ROI)measurements were carried out on a consensus basis by two neuroradiologists(Erkan G?k?e and Murat Beyhan with more than 12 and 7 years of work experience,respectively)who were not aware of the clinical status of the patients.On ADC value measurements,cystic components of the tumors were not considered,the ROI was manually placed through the widest axial section of the tumor,and the ADC value was determined as mm2/s(Figures 2C,3C,4C and 5C).Semiquantitative analysis of DCE MRI was based on TIC(Figures 2D and E,3D and E,4D and E,5D and E).Tpeakwas measured as the time from the point where the lesion began to show contrast enhancement to the point with the highest level of contrast enhancement.TICs were evaluated in four different categories based on Yabuuchiet al[10]:Type A = Tpeak>120 s;type B = Tpeak≤ 120 s,WR ≥ 30%;type C = Tpeak≤ 120 s,WR <30%;and type D = flat TIC.To confirm the accuracy of TIC and perfusion biomarker analyses,ROIs were drawn in a way to avoid the vascular and cystic parts of the tumors.Quantitative perfusion DCE MRI parameters were measured using the Tofts kinetic model[16].For quantitative perfusion MRI analysis,a 3D ROI was placed in the largest solid component of the tumor,and the Ktrans,Kepand Vevalues were generated automatically(Figures 2F-H,3F-H,4F-H,5F-H).

Figure 1 Patient inclusion and exclusion flowchart.

Figure 2 A 72-year-old male patient with a Warthin’s tumor in the right parotid gland.A:On axial plane T2-weighted image,a mildly hypointense(compared to the gland),smooth-contoured mass localized in the center of the parotid gland is observed;B:On the diffusion-weighted image,the mass appears to be hyperintense;C:ADC value was 0.8 × 10-3 mm2/s on the apparent diffusion coefficient map;D:The mass is hyperperfused on color-coded perfusion image;E:Type B time intensity curve shows a washout ratio of 75%;F,G,and H:Ktrans,Kep,and Ve values on quantitative perfusion images were 0.617 min-1,6.438 min-1,and 0.073,respectively.ADC:Apparent diffusion coefficient.

Figure 3 A 38-year-old male patient with a pleomorphic adenoma in the right submandibular gland.A:On sagittal plane T2-weighted image,a hyperintense(compared to the gland),smooth,slightly lobule-contoured mass is observed;B:On contrast-enhanced axial plane T1-weighted image,intense contrastenhancement is observed in the mass;C:The mass is hyperintense on the apparent diffusion coefficient map due to facilitated diffusion(ADC value:1.7 × 10-3 mm2/s);D:The mass is hypoperfused on color coded perfusion image;E:The tumor has type A time intensity curve;F,G,and H:Ktrans,Kep,and Ve values on quantitative perfusion images were 0.211 min-1,0.383 min-1,and 0.593,respectively.ADC:Apparent diffusion coefficient.

Figure 4 A 76-year-old female patient with diffuse large B cell lymphoma and a mass in the left submandibular gland region.A:On axial T2-weighted image,a smooth-contoured mass with homogeneous internal structure and an intensity similar to that of the submandibular gland is observed;B:Axial contrast-enhanced image shows the intense homogeneous contrast enhancement of the mass;C:On the apparent coefficient mapping image,the mass features prominent diffusion restriction(ADC value:0.6 × 10-3 mm2/s);D:The mass is hyperperfused on color coded perfusion image;E:Type B time intensity curve shows a 48% washout ratio;F,G,and H:Ktrans,Kep,and Ve values on quantitative perfusion images were 1.058 min-1,3.391 min-1,and 0.356,respectively.ADC:Apparent diffusion coefficient.

Statistical analysis

Statistical analyses were performed using SPSS 18.0 software(IBM,Chicago,IL,United States)and MedCalc statistical software version 20.009(MedCalc software bvba,Ostend,Belgium).For each parameter,the conformity of the groups to the normal distribution was evaluated by the Shapiro-Wilk test,and the Levene test was used to evaluate the homogeneity of variances.Data are expressed as the mean ± SD or frequency and percent.One-way ANOVA was used for the groups with a normal distribution for comparison of the groups,and Bonferroni correction was applied in multiple comparisons.The Kruskal-Wallis test was used to compare the groups that did not fit the normal distribution,and

Bonferroni correction was applied in multiple comparisons.The area under curve(AUC),sensitivity,specificity,positive predictive value(PPV),negative predictive value(NPV),and cut-off values of diagnostic parameters were calculated for each tumor group by performing receiver operating characteristic curve(ROC)analysis.

RESULTS

The age range of the 67 patients(40 male and 27 female)included in the study was 12-93 years(mean age = 56.9 ± 15.8 years).One patient had three lesions in the salivary glands,while four patients had two and the remaining 62 had one lesion.Thus,73 lesions in 67 patients were evaluated.The majority of the lesions(86.3%)were located in the parotid glands,while a small number(4.1%)originated from minor salivary glands.The locations,numbers,and frequencies of SGTs are shown in Table 1.Approximately two-thirds of the lesions(68.5%)were benign(Figures 2 and 3),and one-third(31.5%)was malignant(Figures 4 and 5).Warthin's tumors(WTs)(36.0%)were the most common BTs,followed by pleomorphic adenomas(PMAs)(28.0%).Of the MTs,squamous cell cancer(47.8%),adenoid cystic cancers(13.0%),and malignant lymphomas(13.0%)were the most common.The numbers of benign and malignant SGTs are provided in Table 2.The ADC values of PMAs were significantly higher than those of WTs,other benign tumors(OBTs),and MTs(P<0.001).However,there was no significant difference in ADC values for OBTs,WTs,and MTs.Significant differences were not found for ADC values of all BTs and MTs.The mean ADC values of SGTs are shown in Table 3.

Figure 5 A 79-year-old male patient with squamous cell cancer in the left parotid gland.A:On axial plane T1-weighted image,an irregularly contoured,hypointense mass involving skin and subcutaneous tissues is observed;B:On axial contrast-enhanced T1-weighted image,the mass shows an intense heterogeneous contrast enhancement;C:ADC value on the apparent diffusion coefficient map was 1.1 × 10-3 mm2/s;D:The mass is heterogeneously hyperperfused on color coded perfusion image;E:Type C time intensity curve shows a 10% washout ratio;F,G,and H:Ktrans,Kep,and Ve values on quantitative perfusion images were 0.993 min-1,1.659 min-1,and 0.612,respectively.ADC:Apparent diffusion coefficient.

Table 1 Locations,numbers,and frequencies of salivary gland tumors

An evaluation of Tpeakvalues of semiquantitative perfusion MRI parameters revealed that PMAs reached Tpeaksignificantly later(mean Tpeak= 202.74 ± 21.48 s)than WTs,OBTs,and MTs while the difference between OBTs and MTs for Tpeakvalues was not significant.WTs reached Tpeaksignificantly earlier than other tumors.With regard to WR,no washout was observed in PMAs.WTs had the highest mean WR value(59% ± 11%),which was significantly different from the mean WR values of MTs and OBTs.PMAs had type A and WTs had type B TIC pattern,while the majority of MTs and OBTs(54.5% and 45.5%,respectively)exhibited type C TIC pattern.Semiquantitative DCE MRI parameters of SGTs are provided in Table 4.

For quantitative perfusion MRI parameters,Ktransvalues of PMAs,WTs,OBTs,and MTs were not significantly different.The Kepvalue of WTs,on the other hand,was significantly higher than those of other tumors(P<0.001).For Vevalue,WTs and OBTs differed significantly from PMAs and MTs(P<0.001).An evaluation of all BTs and MTs showed significant differences for Kepand Vevalues(P<0.05)but not for Ktransvalues.Quantitative DCE MRI parameters of SGTs are shown in Table 3.

The results of ROC analysis and cut-off values used for the parameters of DWI,semiquantitative and quantitative MRI of PMAs,WTs,and malignant SGTs are given in Table 5.

Table 2 Numbers of benign and malignant salivary gland tumors

Table 3 Mean ADC values and quantitative dynamic contrast-enhanced magnetic resonance imaging parameters of different histopathologic salivary gland tumors

Table 4 Semiquantitative dynamic contrast-enhanced magnetic resonance imaging parameters of different histopathologic salivary gland tumors

Table 5 Comparison of pleomorphic adenomas,Warthin’s tumors,and malignant salivary gland tumors by ROC analysis using ADC,Tpeak,Ktrans,Kep,and Ve values

DISCUSSION

In recent years,diffusion MRI has been an essential complement to conventional sequences in the radiological evaluation of SGTs[1,3,6-10,13-15,17,18].Diffusion MRI allows us to evaluate the cellularity in tissues and the changes that physiological processes create on microstructural features.As malignant or benign SGTs include a highly heterogeneous group,their ADC values could also be highly variable.In cell-rich tumors such as WT and lymphoma,ADC values are low,but tumors containing hetero-geneous components such as PMA have higher ADC values[1,13].In many DWI studies involving SGTs,ADC values were reported to be useful in distinguishing BTs and MTs[6-8,17,19-21].However,there are also studies reporting that DWI was not sufficient to make this distinction but ADC values could be useful in distinguishing some subtypes of MTs or BTs[10,22-24].An evaluation of mean ADC values of all BTs and MTs in the present study showed that ADC values of BTs(0.98 × 10-3± 0.43 mm2/s)and MTs(0.95 × 10-3± 0.31 mm2/s)were similar and did not differ significantly.However,when specific tumoral subgroups were evaluated,significant differences were found in the mean ADC values among PMAs,WTs,and MTs(P<0.001).In a ROC analysis using an ADC cut-off value of >1.1 × 10-3mm2/s for PMAs,the AUC,sensitivity,and specificity were 97.7%,100%,and 89.8%,respectively.A ROC analysis of WTs using an ADC cut-off value of ≤ 0.8 × 10-3mm2/s,on the other hand,resulted in AUC,sensitivity,and specificity values of 74.2%,94.4%,and 58.2%,respectively.These values were 54.1%,78.3%,and 44.0%,respectively,for MTs with an ADC cut-off value of >0.7 × 10-3mm2/s.In the present study,the mean ADC value of malignant lymphomas was 0.56 × 10-3± 0.05 mm2/s,which was well below the average ADC value of all MTs.This finding indicated that diffusion MRI could be more useful in distinguishing the subgroups within both BTs and MTs than contributing to a more general distinction between MTs and BTs.

In addition to diffusion MRI,the parameters of semiquantitative DCE MRI(TIC patterns)have also been frequently used in recent years for the differential diagnosis of SGTs[9,12,20,25,26].On DCE MRI,TIC is obtained from signal intensity changes before the contrast agent administration,during the transition of contrast agent from the capillary bed to extravascular-intercellular distance,and during the washing of contrast agent from the tissue[1,18].TIC patterns are correlated with tumor cellularity and vascularity[1,5,18,27].PMAs have progressive contrast-enhancement due to low microvessel count and cellularity-stromal grade,and their washout patterns are mostly negative and,to a lesser degree,in the form of a plateau[1,27].In the present study,type A TIC pattern(curve pattern with progression towards the late phases)was observed in all PMAs.Tpeakvalues ranged from 161.80 s to 251.70 s.The average Tpeakvalue of PMAs(202.74 ± 21.48 s)was significantly longer compared to the Tpeakvalues of all other SGTs.In ROC analysis of PMAs using a cut-off value of Tpeak>120 s,AUC,sensitivity,specificity,PPV,and NPV were 94.7%,100%,89.8%,70.0%,and 100%,respectively.WTs feature rapid contrast enhancement and washout due to their high microvessel count and cellularity-stromal grade.In the present study,type B TIC pattern(Tpeak≤ 120 s,WR ≥ 30%)was observed in all WTs.Tpeakvalues ranged from 10.80 s to 46.40 s,while WR varied from 31% to 75%.The mean Tpeakof WTs(20.26 ± 11.72 s)was significantly shorter than that of other parotid lesions.The average WR value of WTs(59.33% ± 10.99%)was significantly higher than that of any other tumors.In ROC analysis of WTs using a cut-off value of WR >43%,AUC,specificity,and PPV were quite high(98.1%,94.4%,and 90.9%,respectively).Due to their high microvessel count and lower cellularity-stromal grade,MTs have rapid enhancement but their washouts tend to be slower than those of WTs[1,27].In the present study,the mean Tpeakvalue of MTs(60.60 ± 55.78 s)was significantly shorter than the Tpeakvalue of PMAs.The mean WR value of MTs(18.48 ± 18.38%)was significantly lower than that of WTs,but was not different from the mean WR value of OBTs.In ROC analysis of MTs with cut-off values of Tpeak≤ 120 s and WR ≤ 49%,sensitivities were quite high(91.3% and 100%,respectively)but specificities were quite low(36.0% and 30.0%,respectively).A survey of semiquantitative DCE MRI studies in the literature showed that PMAs generally had type A pattern,while WTs had type B and MTs had type C TIC patterns[4,5,26].TIC patterns are considered to have a higher diagnostic accuracy in distinguishing subgroups in SGTs compared to their power to distinguish all BTs from MTs.However,it was mentioned that TIC patterns had higher specificity especially in PMAs and WTs while their specificity in MTs was lower[4,18,25,26].In their study with all SGTs,Lamet al[26]showed that all MTs except lymphomas showed type C TIC pattern(Tpeak<150 s and WR <30%),while 70% of lymphomas had type B TIC pattern(Tpeak<150 s and WR ≥ 30%).Similar to the findings of Lamet al[26],66.7% of lymphomas in the present study showed type B TIC pattern.However,unlike their findings,some other MTs showed types A,B,and D TIC patterns.There are also studies in the literature reporting that all WTs had type B TIC pattern[4,10,12].In accordance with their findings,100% of WTs in the present study featured type B TIC pattern.Subtypes of SGTs in the present study generally had similar TIC patterns to those reported in the literature.

The literature contains several studies on quantitative DCE perfusion MRI parameters(Ktrans,Kep,and Ve)in SGTs[3,14,15,28].In these studies,mean Ktransvalues for PMAs ranged from 0.101 ± 0.069 min?1to 0.217 ± 0.036 min?1,mean Kepvalues from 0.245 ± 0.160 min?1to 0.567 ± 0.048 min?1,and mean Vevalues from 0.360 to 0.590 ± 0.478,while mean Ktransvalues for WTs varied between 0.105 min?1± 0.064 min?1and 0.464 ± 0.036 min?1,mean Kepvalues between 0.729 ± 0.112 min?1and 2.299 ± 1.312 min?1,and mean Vevalues between 0.1439 ± 0.093 and 0.272 ± 0.013.For MTs,mean Ktransvalues varied from 0.130± 0.035 min?1to 0.327 ± 0.030 min?1;mean Kepvalues from 0.463 ± 0.103 min?1to 0.784 ± 0.064 min?1;and mean Vevalues from 0.264 ± 0.119 to 0.445 ± 0.025.In all of these studies in the literature,the Ktransvalues of PMAs were lower than those of other SGTs[3,14,15,28].Xuet al[3]found that the mean Ktransvalue of PMAs was slightly different from that of WTs(P= 0.05).Yabuuchiet al[14]found no significant differences among Ktransvalues of other SGTs.Huanget al[15]found that the Ktransvalues of PMAs were significantly lower than those of other SGTs.Similar to the results of Yabuuchiet al[14],in our study,mean Ktransvalue of PMAs was the lowest among all SGTs,but it was not significantly different from those of other tumors.In the studies by Xuet al[3],Yabuuchiet al[14],and Huanget al[15],the mean Kepvalue was the lowest in PMAs and highest in WTs.Kepvalues of PMAs,WTs,and MTs in the studies of both Xuet al[3]and Yabuuchiet al[14]were significantly different.However,in the study by Huanget al[15],the Kepvalue of only WTs was significantly different from those of other tumors.In another study by Huanget al[28],significant differences were found in Kepvalues between WTs and PMAs,and between WTs and OBTs.Similar to the results of those studies in the literature,the mean Kepvalue in the present study was the lowest in PMAs and highest in WTs,and Kepvalues of PMAs and WTs were significantly different from those of other tumors[3,13,14,28].Xuet al[3],Yabuuchiet al[14],and Huanget al[15]found that mean Vevalues of WTs were significantly lower than those of other tumors.Similar to the results of their studies,the mean Vevalue of WTs in the present study was significantly lower than those of other tumors[3,14].In another study by Huanget al[28],unlike other studies,the Vevalue of WTs and the Vevalues of PMAs and OBTs were found to be significantly lower.In ROC analysis using a cut-off value of Kep≥ 2.44 min?1for WTs,the AUC,sensitivity,and specificity were 97.3%,100%,and 85.5%,respectively.On the other hand,in ROC analysis using a cut-off value of Ve≤ 0.17,quite high AUC,sensitivity,and specificity values(95.8%,100%,and 90.9%,respectively)were obtained.High Kepand low Vevalues in WTs are explained by the limited extravascular and extracellular space in these tumors.As many studies in the literature and the present study revealed,ADC and TIC patterns of WTs could overlap with those of MTs[11].However,similar to the findings of the studies in the literature,the present study showed that quantitative perfusion MRI parameters Kepand Vecould contribute greatly to distinguishing WTs from MTs[3,14,15].Nevertheless,our findings need to be verified by future quantitative perfusion MRI studies performed with larger series.

There are some limitations in this study.First,the parameters(number of dynamic series,acquisition time,etc.)varied on perfusion MRI series due to the retrospective nature of the study.Second,most of the tumors in our study were benign SGTs,and the number of MTs in the primary salivary gland was relatively low,which may have resulted in an overestimation of the diagnostic accuracy.Third,the manual definition of ROI might have increased the variability in quantitative measurements.Although the cystic-necrotic components of the lesions were excluded from the ROI in our study,contamination of these areas can lead to significant changes in quantitative values,even if it is small in manual measurements.Fourth,for the measurements of ADC values and semiquantitative and quantitative DCE perfusion MRI parameters,interobserver agreement could not be evaluated in the study as the measurements were made by two observers with consensus.

CONCLUSION

Combined use of quantitative DCE MRI along with diffusion MRI and semiquantitative DCE MRI could help radiologists in the differential diagnosis of different subtypes of SGTs by providing higher diagnostic accuracy.

ARTICLE HIGHLIGHTS

Research background

Conventional magnetic resonance imaging(MRI)provides more data than other radiological modalities in determining the extent of tumor spread in salivary gland tumors(SGTs)and assessing its relationship to vascular and neural structures,but falls short of distinguishing subtypes of SGTs.As the malignant or benign nature of SGTs affects the treatment protocol,it is important to differentiate between malignant(MTs)and benign tumors(BTs)noninvasively with high diagnostic accuracy.

Research motivation

In recent years,advanced MRI techniques such as diffusion-weighted imaging(DWI)and semiquantitative MRI have been increasingly used in the radiological evaluation of SGTs.However,various studies on quantitative dynamic contrast-enhanced(DCE)perfusion MRI parameters(Ktrans,Kep,and Ve)in SGTs are limited.Therefore,in this study,the effectiveness of advanced MRI applications,including all three methods,in the diagnosis of SGTs was evaluated in light of the literature.

Research objectives

To determine the diagnostic efficiency of DWI and DCE(semiquantitative perfusion)MRI and quantitative perfusion MRI parameters in SGTs.

Research methods

Apparent diffusion coefficient(ADC)values of SGTs on DWI were measured with manually inserted regions of interest,excluding the cystic components of the tumors.Time intensity curve(TIC)patterns were created for semiquantitative perfusion MRI based on Tpeakand washout ratios(WRs)of tumors.On quantitative DCE MRI,perfusion parameters such as Ktrans[volume transfer constant between blood plasma and extracellular extravascular space(EES)],Kep(flux rate constant between the EES and plasma),and Ve(EES fractional volume)were used.

Research results

The ADC values of pleomorphic adenomas(PMAs)were significantly higher than those of Warthin's tumors(WTs),other benign tumors(OBTs),and MTs(P<0.001).However,there was no significant difference in ADC values for OBTs,WTs,and MTs.PMAs had type A and WTs had type B TIC pattern while the vast majority of MTs and OBTs(54.5% and 45.5%,respectively)displayed type C TIC pattern.PMAs showed no washout,while the highest mean WR was observed in WTs.For quantitative perfusion MRI parameters,the Kepvalue of WTs was significantly higher than those of other tumors(P<0.001).For the Vevalue,WTs and OBTs differed significantly from PMAs and MTs(P<0.001).Ktransvalues of PMAs,WTs,OBTs,and MTs were not significantly different.

Research conclusions

DWI and semiquantitative and quantitative perfusion MRI,which provide more information on the microstructure,cellularity,interstitial distance,and vascularity of tumors,have increased the discrimination power for subtypes of SGTs.

Research perspectives

Although there is some overlap in the findings of the subtypes of SGTs obtained by advanced MRI methods,the combined use of DWI and semiquantitative and quantitative perfusion MRI will increase the power for distinguishing subtypes of SGTs.

ACKNOWLEDGEMENTS

We thank Demir O and Gürp?nar AB for their help with the statistical analyses.

FOOTNOTES

Author contributions:G?k?e E designed the study;Beyhan M supervised the study;G?k?e E and Beyhan M participated in literature research and manuscript preparation,and read and approved the final version.

Institutional review board statement:This study was reviewed and approved by the Ethics Committee of the Tokat Gaziosmanpasa University Faculty of Medicine(20-KAEK-105).

Informed consent statement:Patients were not required to give informed consent for the study as figures from picture archiving and communication systems were studied retrospectively.

Conflict-of-interest statement:All the authors report no relevant conflicts of interest for this article.

Data sharing statement:No additional data are available.

Open-Access:This article is an open-access article that was selected by an in-house editor and fully peer-reviewed by external reviewers.It is distributed in accordance with the Creative Commons Attribution NonCommercial(CC BYNC 4.0)license,which permits others to distribute,remix,adapt,build upon this work non-commercially,and license their derivative works on different terms,provided the original work is properly cited and the use is noncommercial.See:https://creativecommons.org/Licenses/by-nc/4.0/

Country/Territory of origin:Turkey

ORCID number:Erkan G?k?e 0000-0003-3947-2972;Murat Beyhan 0000-0002-8630-4632.

S-Editor:Liu XF

L-Editor:Wang TQ

P-Editor:Liu XF