
94% of researchers rate our articles as excellent or good
Learn more about the work of our research integrity team to safeguard the quality of each article we publish.
Find out more
ORIGINAL RESEARCH article
Front. Oncol. , 04 March 2025
Sec. Head and Neck Cancer
Volume 15 - 2025 | https://doi.org/10.3389/fonc.2025.1505385
This article is part of the Research Topic Advancing Cancer Imaging Technologies: Bridging the Gap from Research to Clinical Practice View all 13 articles
Introduction: Accurate differentiation between pleomorphic adenomas (PA) and Warthin tumors (WT) in the parotid gland is challenging owing to overlapping imaging features. This study aimed to evaluate a nomogram combining dual-energy computed tomography (DECT) quantitative parameters and radiomics to enhance diagnostic precision.
Methods: This retrospective study included 120 patients with pathologically confirmed PA or WT, randomly divided into training and test sets (7:3). DECT features, including tumor CT values from 70 keV virtual monochromatic images (VMIs), iodine concentration (IC), and normalized IC (NIC), were analyzed. Independent predictors were identified via logistic regression. Radiomic features were extracted from segmented regions of interest and filtered using the K-best and least absolute shrinkage and selection operator. Radiomic models based on 70 keV VMIs and material decomposition images were developed using logistic regression (LR), support vector machine (SVM), and random forest (RF). The best-performing radiomics model was combined with independent DECT predictors to construct a model and nomogram. Model performance was assessed using ROC curves, calibration curves, and decision curve analysis (DCA).
Results: IC (venous phase), NIC (arterial phase), and NIC (venous phase) were independent DECT predictors. The DECT feature model achieved AUCs of 0.842 and 0.853 in the training and test sets, respectively, outperforming the traditional radiomics model (AUCs 0.836 and 0.834, respectively). The DECT radiomics model using arterial phase water-based images with LR showed improved performance (AUCs 0.883 and 0.925). The combined model demonstrated the highest discrimination power, with AUCs of 0.910 and 0.947. The combined model outperformed the DECT features and conventional radiomics models, with AUCs of 0.910 and 0.947, respectively (P<0.05). While the difference in AUC between the combined model and the DECT radiomics model was not statistically significant (P>0.05), it showed higher specificity, accuracy, and precision. DCA found that the nomogram gave the greatest net therapeutic effect across a broad range of threshold probabilities.
Discussion: The nomogram combining DECT features and radiomics offers a promising non-invasive tool for differentiating PA and WT in clinical practice.
Parotid gland tumors represent approximately 2–3% of all head and neck tumors, with approximately 80% being benign (1). Among them, pleomorphic adenoma (PA) is the most common, constituting a significant portion of parotid gland tumors of epithelial origin. Although typically benign, approximately 3.2% of PAs have the potential for malignant transformation (2, 3). In addition, PA can infiltrate surrounding tissues, forming pseudopods or satellite nodules, making complete excision using a traditional debulking surgery challenging. This results in a risk of recurrence, necessitating treatment via extracapsular dissection (ECD) (4, 5). Warthin tumors (WT) account for approximately 15% of primary parotid tumors, making them the second most common after PA (6). WT predominantly affects middle-aged men and smokers, with some cases involving both parotid glands. Although malignant transformation and recurrence are very rare (7), WT often triggers inflammatory reactions in surrounding tissues, complicating surgical intervention. Consequently, conservative treatment is typically preferred to avoid unnecessary surgical risks (8, 9). Therefore, differentiating between PA and WT is crucial because it directly influences clinical decision-making and individualized surgical planning.
Fine-needle aspiration Biopsy (FNAB) is the current “gold standard” for tumor diagnosis. However, it is an invasive procedure with risks, including potential damage to the facial nerve, which can lead to peripheral facial paralysis (10, 11). Ultrasound is effective in detecting tumors in the superficial lobes of the parotid gland, but visualizing deeper lesions can be challenging, and its accuracy is highly dependent on the operator’s experience (12). While magnetic resonance imaging (MRI) offers high tissue resolution, its use is limited by factors such as prolonged examination times and sensitivity to metal implants, making it less ideal as a primary diagnostic tool (13–16). Traditional computed tomography scans provide detailed information about the tumor’s anatomical location and lesion extent and can uniquely assess calcification and adjacent vascular conditions (17, 18). However, traditional CT relies on attenuation-based imaging, which lacks material-specific resolution.
Dual-energy CT (DECT) addresses this limitation by using dual-energy spectra acquired simultaneously during CT scanning to improve the analysis of physical and chemical properties. One common method for obtaining dual-energy spectra is the instantaneous switching between high and low tube voltages (140 kVp and 80 kVp) during data acquisition. DECT works by converting the absorption data of the same material at different energy levels into density data based on the variation in attenuation coefficients of different substances exposed to varying X-ray energy levels. This process produces material decomposition (MD) images (19). Given that iodine is the primary material used in contrast agents for enhanced scanning, measuring iodine concentration on water-iodide MD images allows for a detailed visualization of iodine distribution, which can reflect blood perfusion at lesion sites (20).
Radiomics involves using machine learning techniques to extract information from medical images that are not visible to the naked eye, capturing the microscopic characteristics of lesions. This allows for predicting the biological behavior of lesions, as well as patient prognosis and treatment outcomes. As a novel biomarker, radiomics is valuable in differentiating tumor types by assessing their heterogeneity (21). Radiomics has shown great promise in improving noninvasive individualized treatment and monitoring by differentiating parotid tumor types, evaluating the integrity of benign tumor envelopes, grading malignant tumors, and differentiating lymphoid lesions, and these studies use a variety of imaging techniques such as CT, MRI, US, and PET-CT (22–26). However, the current issue is the small number and varying quality of studies in this direction, with none of them split at several time points, prospectively verified, or cost-effective. The majority of the studies did not disclose public data, and one of the issues was the difficulty in reproducing the study (27). As a result, the clinical application needs to be further investigated. In this work, we included the energy spectrum CT substance separation image radiomics features and compared them to the performance of traditional CT histology, with the goal of providing evidence for the validity of radiomics in the parotid gland using a novel method.
The integration of radiomics with DECT has been increasingly applied to disease identification and prediction across various organ systems (28–30). However, its application to the parotid gland remains limited. In one study, Ajmi et al. used virtual monochromatic images from DECT combined with radiomics analysis to differentiate between PA and WT. They found that combining radiomics features from monochromatic images at multiple energy levels significantly improved diagnostic accuracy compared to using a radiomics model based on a single energy level (65 keV) (31). However, the study did not explore the potential value of radiomics features derived from MD images. We hypothesized that radiomics based on MD images from DECT, which quantitatively analyze iodine concentration differences between lesions, along with microscopic image analysis, could provide a novel approach for identifying parotid tumors in clinical practice. Therefore, this study aimed to assess the value of DECT MD image-based radiomics and evaluate the effectiveness of a nomogram that combines radiomics with DECT-independent predictors in differentiating the two most common parotid tumors, playing an adjunctive role in the development of patient treatment strategies using a non-invasive modality.
This retrospective study, which did not interfere with the clinical diagnosis or treatment of any patients, was approved by the Ethics Committee of the First Affiliated Hospital of Dalian Medical University (approval number: PJ-KS-KY-2024-230; date: April 3, 2024), with a waiver for written informed consent. We collected data from 128 patients diagnosed with PA and WT of the parotid gland who presented to our institution and underwent DECT enhancement between October 2021 and October 2023. The inclusion criteria were as follows: (1) primary parotid tumors with all types confirmed by surgical pathology; (2) no prior treatment before the CT examination; and (3) high-quality CT images without factors affecting the measurement of DECT parameters. The exclusion criteria were as follows: (1) history of chemotherapy or radiotherapy for parotid tumors; (2) poorly visible lesions on imaging; (3) incomplete clinical data; and (4) severe metallic or motion artifacts in images. After applying these criteria, we excluded one patient with extensive metal (denture) artifacts obscuring the tumor, two with poorly visualized lesions, and five with missing clinical data. This resulted in a total of 120 patients being included in the study, comprising 66 cases of PA and 54 cases of WT. The patients were randomly assigned to training and test sets in a 7:3 ratio, with 84 cases in the training set and 36 cases in the test set. The patient selection process is illustrated in Figure 1.
All patients underwent non-contrast and contrast-enhanced scans from the superior orbital rim to the level of the aortic arch using a GE Discovery CT 750 HD scanner. The scanning parameters were as follows: instantaneous switching between high (140 kVp) and low (80 kVp) tube voltages, a tube current of 600 mA, a gantry rotation speed of 0.6 seconds per rotation, a collimator width of 64 × 0.625 mm, and a slice thickness and spacing of 2.5 mm. The contrast agent used was iohexol (300 mg I/mL) administered at a dose of 100–200 mL via an elbow vein injection. DECT scans in the arterial phase (AP) and venous phase (VP) were initiated at 30 s and 60 s, respectively, after contrast injection. Iodine- and water-based MD sequences of the non-contrast (NP), arterial, and venous phases were reconstructed using Gemstone Spectral Imaging (GSI) post-processing software. The images were then imported into the picture archiving and communication system (PACS) for analysis.
The study workflow is outlined in Figure 2. Image processing and data measurements were performed using GSI Viewer 4.5 software (GE Healthcare). Key features such as the number, location, maximum diameter of the lesions, cystic changes, calcifications, and lymph node status were recorded and analyzed (Table 1). As illustrated in Figure 3, we examined 70 keV virtual monochromatic images (VMIs) and MD images (iodine-based and water-based). A radiologist with three years of experience in head and neck imaging, blinded to all pathological findings, measured the DECT quantitative parameters by placing regions of interest (ROIs) within the tumor and the ipsilateral common carotid artery. The ROIs were positioned to encompass the entire tumor while avoiding calcifications and prominent enhanced vessels. Measurements were centered at the level of maximum lesion display, averaged across three consecutive slices. Lesion attenuation (in Hounsfield units, HU) was obtained from the 70 keV VMIs for the NP, AP, and VP. Iodine concentrations for the tumor and common carotid artery were assessed using iodine-based MD images. The degree of enhancement and normalized iodine concentration (NIC) were calculated using the following formulas:
Figure 2. Overview of the study workflow of this study. DECT, dual-energy computed tomography; IC, iodine concentration; LASSO, least absolute shrinkage and selection operator; NIC, normalized iodine concentration; ROI, region of interest; VOI, volume of interest.
Figure 3. Representative dual-energy CT (DECT) material decomposition (MD) images. (white arrows, A-D) PA in a 67-year-old woman. (red arrows, E-H) WT in a 71-year-old man. (A, E) Iodine-based MD images in the arterial phase. (B, F) Water-based MD images in the arterial phase. (C, G) Iodine-based MD images in the venous phase. (D, H). Water-based MD images in the venous phase.
The largest lesion was selected for analysis in cases with multiple parotid lesions. Differences between WT and PA were analyzed using univariate logistic regression. Significant variables (p<0.05) were included in a multivariate logistic regression model to identify independent predictors of DECT features and calculate odds ratios (ORs) with 95% confidence intervals (CIs).
Lesions were manually segmented using the open-source software 3D Slicer (https://www.slicer.org). ROIs were delineated slice by slice to generate the volume of interest (VOI), carefully avoiding adjacent tissues and blood vessels. Layers heavily affected by artifacts, such as those near the mandible, were excluded from the segmentation. To assess intra- and inter-observer consistency, 30 patients (18 PA, 12 WT) were randomly selected, and their VOIs were re-segmented after one month by the same radiologist to evaluate intra-observer agreement. Another radiologist independently repeated the segmentation to assess inter-observer agreement. The reproducibility of radiomic features was quantified using the intraclass correlation coefficient (ICC).
The segmented VOIs were imported into the United Imaging Intelligence workstation (Shanghai, China). All sequences were aligned and normalized in terms of window width and level. A total of 2,048 features were extracted from each VOI, categorized into (1) Shape-based (2D/3D) features; (2) First-order histogram features; and (3) Second-order texture features, including gray level run length matrix (GLRLM), gray level co-occurrence matrix (GLCM), gray level dependence matrix (GLDM), gray level size zone matrix (GLSZM), and neighboring gray-tone difference matrix (NGTDM). These first- and second-order features were extracted not only from the original image but also from derived images processed with 14 filters, such as Mean, Bilateral, and Speckle Noise. The algorithms used for feature extraction follow the Image Biomarker Standardization Initiative (IBSI) guidelines (32).
Features with intra- and inter-observer ICC values greater than 0.75 were retained and standardized using z-score normalization. The optimal regularization parameters were determined through five-fold cross-validation in the training set. Univariate analysis was conducted using the K-best method with F-value tests to identify significant differences between classifications, selecting the top 10 most important features. These features were further refined using the least absolute shrinkage and selection operator (LASSO) analysis based on the feature coefficients at the maximum alpha value. The final model was trained using the selected features, and its performance and generalization ability were evaluated with the test set.
Predictive models were developed using logistic regression (LR), support vector machine (SVM), and random forest (RF) algorithms based on independent DECT predictors and filtered radiomic features. The diagnostic performance of each model was assessed by comparing the area under the curve (AUC) of receiver operating characteristic (ROC), as well as sensitivity, specificity, accuracy, and precision. Significant DECT features were utilized to construct the DECT-based model. Among the radiomics model, the best-performing model from the 70 keV non-contrast and contrast-enhanced sequences was selected as the conventional radiomics model. In contrast, the top-performing model based on MD images was selected as the DECT radiomics model. A combined model was developed by calculating the Rad-score from the features of the best MD model and integrating it with independent DECT predictors. A nomogram was constructed for this combined model. Calibration curves were generated to determine the goodness-of-fit between predicted and actual values across the four models. Decision curve analysis (DCA) was performed to assess the net benefit rates of each model, providing insights into their clinical utility.
Statistical analyses were conducted using SPSS software (version 26.0, IBM) and MedCalc software (version 22.009, https://www.medcalc.org). Quantitative parameters following a normal distribution are presented as mean ± standard deviation (x ± s) and were compared using the Student’s t-test. Non-normally distributed data are expressed as median and interquartile range [M (P25, P75)] and were analyzed using the Mann–Whitney U test to assess between-group differences. Qualitative data were analyzed using the chi-square test or Fisher’s exact test, depending on the data characteristics. Univariate and multivariate logistic regression analyses were employed to evaluate all parameters, with significant variables sequentially filtered to identify independent predictors. The Delong test was used to compare the differences in AUC values between the models, with a p-value of less than 0.05 considered statistically significant.
A total of 120 patients (78 males, 42 females) were included in this study, ranging in age from 24 to 78 years. Among them, 66 had PA and 54 WT. Detailed clinical data for both groups are presented in Table 1. Statistically significant differences were observed between PA and WT patients in terms of age, sex, smoking history, number of tumors, presence of cystic lesions, maximum lesion diameter, and ipsilateral lymph node status (all p<0.05).
As presented in Supplementary Table S1, no significant differences in DECT features were observed between the training and test sets. Within the training set, significant differences were observed between PA and WT in attenuation (NP and AP), enhancement degree (AP and VP), IC, and NIC (all p<0.05). Specifically, NP attenuation, AP attenuation, AP enhancement degree, IC, and NIC were significantly higher in WT than in PA. However, the VP enhancement degree, IC, and NIC were notably lower in WT (Table 2).
Logistic regression analysis (Table 3) identified IC (VP), NIC (AP), and NIC (VP) as independent predictors for distinguishing PA from WT. When these features were combined to construct the model, no significant performance differences were observed among the LR, SVM, and RF models (Supplementary Table S2 and Supplementary Figure S1). Given its relatively better performance in the test set, the LR model was ultimately selected as the DECT feature model.
From the 70 keV VMIs across the non-contrast, arterial, and venous phases, 784, 398, and 536 highly reproducible features (ICC>0.75) were extracted, respectively. Subsequent K-best and LASSO analyses retained two, four, and four non-zero coefficient features for each phase, respectively. Supplementary Table S3 provides details on the selected features and their corresponding LASSO coefficients. Diagnostic models were constructed using LR, SVM, and RF algorithms based on these selected features, with model performance detailed in Supplementary Table S4 and illustrated in Supplementary Figures S2-S4.
Among the models, the SVM model for the AP images demonstrated superior performance as the representative conventional radiomics model, achieving an AUC of 0.834 (95% CI, 0.692–0.977). The LASSO filtering process for the AP radiomics features is shown in Figure 4A.
Figure 4. LASSO coefficient variation with hyperparameters for 10 features screened using the K-best method in AP 70kev images (A) and AP iodine-based MD images (B). AP, arterial phase; LASSO, least absolute shrinkage and selection operator; VP, venous phase.
In the Iodine-based image analysis, we extracted 540, 884, and 674 features with ICCs greater than 0.75 in the non-contrast, arterial, and venous phases, respectively. After screening, the top four features were selected for each phase, resulting in AUCs ranging from 0.767 to 0.907 for the training models and 0.638 to 0.703 for the test models. For the water-based image analysis, we extracted 647, 639, and 548 features with ICCs greater than 0.75 across the three phases. The final selection included two, three, and one key features, achieving AUCs ranging from 0.843 to 0.944 in the training models and from 0.823 to 0.925 in the test models.
No significant differences were found among the LR, SVM, and RF models for any sequence. The LR model of the arterial-phase water-based images, which demonstrated the highest diagnostic accuracy in the test set, was selected as the final DECT radiomics model (Supplementary Table S5, Supplementary Figures S5-S10) with an AUC of 0.925 (95% CI, 0.844–0.999). The LASSO selection process for the arterial-phase water-based images is shown in Figure 4B, while Figure 5 illustrates Rad-score differences among patients.
Figure 5. Rad-scores for all included samples in the training (A) and test (B) sets. Label 1 represents pleomorphic adenoma (PA), and Label 2 represents Warthin’s tumor (WT).
The Rad-score formula is as follows:
The nomogram for the combined model was constructed using IC (VP), NIC (AP), NIC (VP), and Rad scores derived from the arterial-phase water-based images (Figure 6A). The performance of each model is shown in Supplementary Table S6 and Supplementary Figure S11. Heat maps indicate that the Rad score contributed the highest weight among the features (Figures 6B–D). A comparative analysis of the DECT features, traditional radiomics, and DECT radiomics is presented in Table 4, summarizing the combined model’s performance metrics. The ROC curves for all four models are depicted in Figure 7. The integrated model achieved AUCs of 0.91 in the training set and 0.947 in the test set, clearly outperforming both the DECT feature and traditional radiomics models. Although the AUC of the composite model was not significantly different from that of the DECT radiomics model, it exhibited greater specificity, accuracy, and precision. Calibration curves confirmed good alignment between predicted and actual outcomes for the training and test data (Figure 7). The DCA (Figure 8) demonstrates that the combined nomogram offers a higher net benefit across a wide range of threshold probabilities.
Figure 6. Radiomics nomogram was developed in the training set by combining IC(VP), NIC(AP), NIC(VP), and Rad-score. (B, C) Correlation heat maps for the combined model. (D) Unsupervised hierarchical clustering of the combined model’s features. Feature 1: IC(VP); Feature 2: NIC(AP); Feature 3: NIC(VP); Feature 4: Rad-score. AP, arterial phase; IC, iodine concentration; NIC, normalized iodine concentration; VP, venous phase.
Figure 7. Receiver operator characteristics (ROC) and calibration curves for the DECT features model, traditional radiomics model, DECT radiomics model, and nomogram in the training (A, C) and test (B, D) sets. DECT, DECT, dual-energy CT.
Figure 8. Decision curve analysis (DCA) for the classification of PA and WT in training (A) and test (B) sets. The nomogram provides the greatest net benefit compared to the other three datasets. PA, pleomorphic adenoma; PA, pleomorphic adenoma WT, Warthin’s tumor.
Age, sex, smoking history, and cystic degeneration all showed significant differences between adenolymphomas and pleomorphic adenomas (P<0.001). These findings are in line with earlier research that found adenolymphomas are more common in middle-aged and older men and that nicotine from tobacco is a significant pathogenetic factor in their development. Because of their faster growth rate, adenolymphomas are more likely to experience cystic degeneration and necrosis (33). The imaging presentations of various parotid gland tumors often overlap, making accurate diagnosis using conventional CT and MRI challenging. This diagnostic process heavily relies on the clinician’s expertise, limiting reliability. The diagnostic sensitivity of salivary gland tumors for FNAC varies greatly, ranging from 71 to 93%, according to several earlier studies (34). Aspiration is an invasive test, and there is a risk of facial nerve damage, which is a drawback that is hard to totally prevent. Previous studies have shown promising results by combining DECT with radiomics for conditions such as pancreatic and lung tumors (28–30). However, to our knowledge, the application of MD techniques combined with radiomics for parotid tumor classification has not been explored. In our study, we developed an MD image-based radiomics model. We constructed a nomogram by integrating DECT quantitative parameters with MD radiomics features to distinguish PA from WT in the parotid gland. Our results indicate that MD radiomics provides superior diagnostic accuracy compared to DECT features and traditional radiomics alone. The combined model showed further improved diagnostic accuracy and clinical utility over the standalone DECT feature, traditional radiomics, and DECT radiomics models. These results suggest that the nomogram could serve as a valuable tool for the noninvasive preoperative evaluation of PA and WT. Because the two cancers are treated differently, using the nomograms to distinguish between them helps with treatment decision-making.
Additionally, our study underscored the utility of DECT quantitative parameters in distinguishing between PA and WT. Specifically, IC (VP), NIC (AP), and NIC (VP) were identified as independent predictors. The IC obtained from iodine-based DECT images reflects tumor blood perfusion. We found that the NIC (AP) of WT was higher than that of PA, whereas IC (VP) and NIC (VP) were significantly lower, consistent with previous studies (35, 36). This phenomenon may be related to the pathological structures of the two tumors. WT is a well-vascularized benign tumor with a dense microvascular network and leaky blood vessels, resulting in rapid contrast agent uptake during the arterial phase and quick washout during the venous phase (37). In contrast, PA contains epithelial and mesenchymal components with abundant intercellular mucin, leading to slower contrast uptake due to limited vascularity (18). The combined AUCs of the three independent predictors for the training and test sets were 0.842 and 0.853, respectively, demonstrating the feasibility of using DECT quantitative parameters to differentiate PA from WT.
The 70 keV VMI can simulate the attenuation characteristics of conventional 120 kVp scans in head and neck imaging, making them suitable alternatives for radiomic analysis in this region (38). In our study, radiomic analysis using 70 keV images demonstrated superior diagnostic efficacy during the arterial phase compared to other phases, with a sensitivity and specificity of 0.812 and 0.75, respectively, in the test set. Previous studies have highlighted the potential of traditional CT radiomics in distinguishing PA from WT (23, 39). Feng et al. identified the arterial phase as the most informative for differentiating these tumors, with their findings generally aligning with our results (22). Similarly, Jung et al. found no significant difference between quantitative enhancement analysis and radiomics in distinguishing PA from WT (40). This aligns with our results, where there was no significant difference between the two methods, although traditional CT radiomics performed slightly worse than the DECT feature model.
Radiomics enhances the analysis of MD images in DECT by quantifying the distribution of substances such as iodine and water into voxel relationships, offering a more precise and scientific. In our study, after analyzing radiomic features from both iodine- and water-based images across different phases, we found that the water-based images in the arterial phase yielded the highest diagnostic performance. The DECT radiomics model significantly improved the differential diagnosis of PA and WT compared to traditional CT radiomics. Water-based imaging, achieved through post-processing techniques that remove iodine signals, reflects the relative water distribution in tissues. We hypothesize that variations in the amount and distribution of water within parotid malignancies led to the maximum diagnostic efficacy of water-based images, which are based on the imaging principle of eliminating iodine and emphasizing water. The fact that Michela Gabelloni et al. were able to identify pleomorphic adenomas and adenolymphomas of the parotid gland using radiomics taken from T2WI sequences in MRI with similarly high sensitivity, specificity, and accuracy supports the significant difference in water content between the two benign tumors (41). This discrepancy might be explained by the fact that adenolymphomas, which have closely packed cells and many mesenchymal components, including tiny blood vessels, do not have a significant water content, while the majority of pleomorphic adenomas have a high and widely distributed water content due to their high mucus component content (42). This approach, often referred to as virtual non-enhanced imaging, closely mirrors the resolution and contrast of conventional non-enhanced image (43, 44). Notably, in the present study, the radiomic features derived from water-based images in the arterial phase outperformed those from the non-enhanced phase. We hypothesize that the reason for this phenomenon may be that virtual non-contrast images are derived from enhanced images that are consistent with both the position taken by the patient during scanning and the level of the tumor display, whereas real non-enhanced images are difficult to match exactly, leading to differences in final model performance. This also suggests that virtual non-enhanced images could serve as a viable substitute for conventional non-enhanced scans, thereby reducing redundant imaging and minimizing patient radiation exposure. Interestingly, the radiomic features from iodine-based images across all phases were less effective and even inferior to those of traditional CT radiomics. This may be attributed to the fact that iodine-based images primarily reflect iodine distribution, neglecting the contributions of other substances. As a result, the contrast between tissues and lesions is limited, leading to suboptimal image resolution and weaker diagnostic performance. And it also indicates that iodine might have interfered with the contrast-enhanced image histology model, making it perform marginally worse overall than the virtual non-enhanced images.
Studies have demonstrated that LASSO-based feature selection offers significant advantages in the differential diagnosis of parotid tumors (45). In our study, three key radiomic features—firstorder-90Percentile, glcm-Joint Energy, and glrlm-Low Gray Level Run Emphasis—were identified from arterial-phase water-based images as the most robust discriminators. These features were extracted from images processed using boxmean and speckle noise filters, which enhance edge definition and texture details. The boxmean filter implements fast rectangular homogenization, while the speckle noise filter introduces noise proportional to pixel intensity, thereby refining texture analysis. The first-order features capture voxel intensity distributions within a ROI, with the 90Percentile representing the voxel intensity at the 90th percentile. GLCM features, such as Joint Energy, describe the spatial relationships between voxel pairs, where higher energy values indicate greater uniformity in gray-level distribution. GLRLM, such as Low Gray Level Run Emphasis measures the prevalence of consecutive low-intensity pixels, reflecting textural uniformity and heterogeneity (46). These features can indicate tumor heterogeneity through the uniformity of the texture. Consistent with prior studies, WT exhibits higher texture heterogeneity owing to its complex cellular makeup, which includes epithelial components and lymphoid stroma, resembling the patterns often seen in malignant tumors (37, 47). Furthermore, the higher Low Gray Level Run Emphasis in WT than in PA suggests a greater presence of low-intensity voxels, likely due to the cystic and necrotic regions commonly found in WT. This observation aligns with our findings, which show that cystic lesions are more frequent in WT than in PA, as noted in the clinical characteristics of our patient cohort.
LR is one of the most commonly used machine learning algorithms, applying the principles of linear regression of continuous variables and using a logistic function to predict categorical outcomes. In contrast, RF algorithm, which relies on ensemble learning via decision trees, often underperforms when applied to small datasets owing to attribute perturbation, leading to large generalization errors (48, 49). SVMs are generalized linear classifiers that binarily classify data in a supervised learning manner and can use slack variables to avoid the influence of linearly indistinguishable data on the results; thus, a smaller set of samples achieves the greatest possible effectiveness (50, 51). In our study, although the differences in performance among the LR, SVM, and RF models were not significant, the LR model exhibited slightly higher AUC values across both the DECT features and radiomics models. This suggests that the classification problem at hand is primarily linear, explaining the marginally lower performance of the nonlinear RF and SVM models. To our knowledge, this is the first study comparing the effectiveness of multiple machine-learning models in distinguishing between PA and WT parotid tumors. A previous study by Yu et al. indicated that radiomic features from enhanced images often follow nonlinear patterns, favoring SVM and RF models over LR (45). However, Zhiying et al. found that the XGBoost algorithm outperformed SVM in classifying radiomics features derived from T1WI and T2WI MRI sequences (52).
Our study has several limitations. First, it is based on a single-center dataset with a small sample size, which restricts the robustness of the test results. Multi-center studies with larger datasets are needed to validate these findings. Second, we compared only three algorithms—LR, RF, and SVM—while other machine learning techniques, such as k-nearest neighbor (KNN) or XGBOOST, were not evaluated. Third, our focus was limited to MD technology in radiomics, excluding other DECT techniques, such as virtual monochromatic imaging, which warrant further investigation. Furthermore, there are few prior studies on radiomics in the parotid gland, the evidence is mixed, and there are variations in scanning parameters between clinical practice and data collection. As a result, it is challenging to create a consistent standard, and there are some restrictions when applying it in a clinical setting. In the future, we will increase the sample size for additional validation. Finally, we analyzed only two common benign parotid tumors, excluding other tumor types. Future studies should investigate multiple parotid tumors to provide a more comprehensive analysis.
DECT quantitative features are valuable in discriminating between the two most common parotid tumors. Additionally, water-based radiomics models derived from arterial phase images outperform traditional radiomics models. Our novel nomogram, which integrates DECT quantitative features with arterial phase water-based radiomics, further improves diagnostic accuracy. As a non-invasive preoperative tool, this model provides more precise guidance for individualized clinical treatment strategies.
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
The studies involving humans were approved by Ethics Committee of the First Affiliated Hospital of Dalian Medical University (approval number: PJ-KS-KY-2024-230; date: April 3, 2024). The studies were conducted in accordance with the local legislation and institutional requirements. The ethics committee/institutional review board waived the requirement of written informed consent for participation from the participants or the participants’ legal guardians/next of kin because the study was retrospective, did not involve any clinical trials on the subjects, did not contain patient-sensitive information, did not involve commercial interests, and therefore exempted the subjects from informed consent.
ZG: Conceptualization, Formal analysis, Investigation, Methodology, Visualization, Writing – original draft. JL: Resources, Writing – review & editing. YH: Validation, Writing – review & editing. SC: Writing – review & editing. LW: Funding acquisition, Resources, Writing – review & editing.
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by Liaoning Medical Education Research Program, China, Grant number: 2022-N005-05.
We would like to express our gratitude to all those who helped us in conducting this study. We thank all subjects who participated in this study; the Department of Radiology of the First Affiliated Hospital of Dalian Medical University for creating a favorable research environment for our study; and Editage (www.editage.cn) for English language editing.
Author JL was employed by GE HealthCare.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declare that no Generative AI was used in the creation of this manuscript.
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2025.1505385/full#supplementary-material
AP, Arterial Phase; AUC, Area Under the Curve; CI, Confidence Interval; CT, Computed Tomography; DCA, Decision Curve Analysis; DECT, Dual-Energy Computed Tomography; GLCM, Gray-Level Co-occurrence Matrix; GLRLM, Gray-Level Run Length Matrix; IC, Iodine Concentration; KNN, k-Nearest Neighbor; LASSO, Least Absolute Shrinkage and Selection Operator; LR, Logistic Regression; MD, Material Decomposition; NIC, Normalized Iodine Concentration; PA, Pleomorphic Adenoma; RF, Random Forest; ROC, Receiver Operating Characteristic; ROI, Region of Interest; SVM, Support Vector Machine; T1WI, T1-Weighted Imaging; T2WI, T2-Weighted Imaging; VMI, Virtual Monochromatic Image; VP, venous phase.
1. Guerra G, Testa D, Montagnani S, Tafuri D, Salzano FA, Rocca A, et al. Surgical management of pleomorphic adenoma of parotid gland in elderly patients: Role of morphological features. Int J Surg. (2014) 12 Supplement 2:S12–6. doi: 10.1016/j.ijsu.2014.08.391
2. Hernandez-Prera JC, Skálová A, Franchi A, Rinaldo A, Vander Poorten V, Zbären P, et al. Pleomorphic adenoma: The great mimicker of Malignancy. Histopathology. (2021) 79:279–90. doi: 10.1111/his.14322
3. Levyn H, Subramanian T, Eagan A, Katabi N, Lin O, Badillo ND, et al. Risk of carcinoma in pleomorphic adenomas of the parotid. JAMA Otolaryngol Head Neck Surg. (2023) 149:1034–41. doi: 10.1001/jamaoto.2023.3212
4. Iro AK, Agaimy A, Müller SK, Sievert M, Iro H, Mantsopoulos K. Satellite nodules in pleomorphic adenomas of the parotid gland: A nightmare for less invasive parotid surgery? Oral Oncol. (2021) 115:105218. doi: 10.1016/j.oraloncology.2021.105218
5. Quer M, Vander Poorten V, Takes RP, Silver CE, Boedeker CC, de Bree R, et al. Surgical options in benign parotid tumors: A proposal for classification. Eur Arch Otorhinolaryngol. (2017) 274:3825–36. doi: 10.1007/s00405-017-4650-4
6. Aoki R, Tanaka T. Pathogenesis of Warthin’s tumor: Neoplastic or non-neoplastic? Cancers (Basel). (2024) 16:912. doi: 10.3390/cancers16050912
7. Patel DK, Morton RP. Demographics of benign parotid tumours: Warthin’s tumour versus other benign salivary tumours. Acta Oto-Laryngol. (2016) 136:83–6. doi: 10.3109/00016489.2015.1081276
8. Psychogios G, Bohr C, Constantinidis J, Canis M, Vander Poorten V, Plzak J, et al. Review of surgical techniques and guide for decision making in the treatment of benign parotid tumors. Eur Arch Otorhinolaryngol. (2021) 278:15–29. doi: 10.1007/s00405-020-06250-x
9. Quer M, Hernandez-Prera JC, Silver CE, Casasayas M, Simo R, Vander Poorten V, et al. Current trends and controversies in the management of Warthin tumor of the parotid gland. Diagnostics (Basel). (2021) 11:1467. doi: 10.3390/diagnostics11081467
10. Fisher R, Ronen O. Cytologic diagnosis of parotid gland Warthin tumor: Systematic review and meta-analysis. Head Neck. (2022) 44:2277–87. doi: 10.1002/hed.27099
11. Witt BL, Schmidt RL. Ultrasound-guided core needle biopsy of salivary gland lesions: A systematic review and meta-analysis. Laryngoscope. (2014) 124:695–700. doi: 10.1002/lary.24339
12. Rubini A, Guiban O, Cantisani V, D’Ambrosio F. Multiparametric ultrasound evaluation of parotid gland tumors: B-mode and color Doppler in comparison and in combination with contrast-enhanced ultrasound and elastography. A case report of a misleading diagnosis. J Ultrasound. (2021) 24:337–41. doi: 10.1007/s40477-020-00469-4
13. Luna LP, Coffey W 3rd, Alvin MD, Shanechi AM, Sankaran N, Rodriguez EF, et al. Parotid Warthin’s tumor: Novel MR imaging score as diagnostic indicator. Clin Imaging. (2022) 81:9–14. doi: 10.1016/j.clinimag.2021.09.005
14. Maraghelli D, Pietragalla M, Cordopatri C, Nardi C, Peired AJ, Maggiore G, et al. Magnetic resonance imaging of salivary gland tumours: Key findings for imaging characterisation. Eur J Radiol. (2021) 139:109716. doi: 10.1016/j.ejrad.2021.109716
15. Friedman E, Cai Y, Chen B. Imaging of major salivary gland lesions and disease. Oral Maxillofac Surg Clin North Am. (2023) 35:435–49. doi: 10.1016/j.coms.2023.02.007
16. Abdel Razek AAK, Mukherji SK. State-of-the-art imaging of salivary gland tumors. Neuroimaging Clin N Am. (2018) 28:303–17. doi: 10.1016/j.nic.2018.01.009
17. Zuo H. The clinical characteristics and CT findings of parotid and submandibular gland tumours. J Oncol. (2021) 2021:8874100. doi: 10.1155/2021/8874100
18. Kato H, Kawaguchi M, Ando T, Mizuta K, Aoki M, Matsuo M. Pleomorphic adenoma of salivary glands: Common and uncommon CT and MR imaging features. Jpn J Radiol. (2018) 36:463–71. doi: 10.1007/s11604-018-0747-y
19. Fernández-Pérez GC, Fraga Piñeiro C, Oñate Miranda M, Díez Blanco M, Mato Chaín J, Collazos Martínez MA. Dual-energy CT: Technical considerations and clinical applications. Radiol (Engl Ed). (2022) 64:445–55. doi: 10.1016/j.rxeng.2022.06.003
20. Sananmuang T, Agarwal M, Maleki F, Muthukrishnan N, Marquez JC, Chankowsky J, et al. Dual energy computed tomography in head and neck imaging: Pushing the Envelope. Neuroimaging Clin N Am. (2020) 30:311–23. doi: 10.1016/j.nic.2020.04.003
21. Mayerhoefer ME, Materka A, Langs G, Häggström I, Szczypiński P, Gibbs P, et al. Introduction to radiomics. J Nucl Med. (2020) 61:488–95. doi: 10.2967/jnumed.118.222893
22. Feng B, Wang Z, Cui J, Li J, Xu H, Yu D, et al. Distinguishing parotid polymorphic adenoma and Warthin tumor based on the CT radiomics nomogram: A multicenter study. Acad Radiol. (2023) 30:717–26. doi: 10.1016/j.acra.2022.06.017
23. Zheng M, Chen Q, Ge Y, Yang L, Tian Y, Liu C, et al. Development and validation of CT-based radiomics nomogram for the classification of benign parotid gland tumors. Med Phys. (2023) 50:947–57. doi: 10.1002/mp.16042
24. Li S, Su X, Ning Y, Zhang S, Shao H, Wan X, et al. CT based intratumor and peritumoral radiomics for differentiating complete from incomplete capsular characteristics of parotid pleomorphic adenoma: A two-center study. Discovery Oncol. (2023) 14:76. doi: 10.1007/s12672-023-00665-8
25. Zhang MH, Hasse A, Carroll T, Pearson AT, Cipriani NA, Ginat DT. Differentiating low and high grade mucoepidermoid carcinoma of the salivary glands using CT radiomics. Gland Surg. (2021) 10:1646–54. doi: 10.21037/gs-20-830
26. Zheng YM, Xu WJ, Hao DP, Liu XJ, Gao CP, Tang GZ, et al. A CT-based radiomics nomogram for differentiation of lympho-associated benign and Malignant lesions of the parotid gland. Eur Radiol. (2021) 31:2886–95. doi: 10.1007/s00330-020-07421-4
27. Aringhieri G, Fanni SC, Febi M, Colligiani L, Cioni D, Neri E. The role of radiomics in salivary gland imaging: A systematic review and radiomics quality assessment. Diagnostics. (2022) . 12:3002. doi: 10.3390/diagnostics12123002
28. Yao Q, Liu M, Yuan K, Xin Y, Qiu X, Zheng X, et al. Radiomics nomogram based on dual-energy spectral CT imaging to diagnose low bone mineral density. BMC Musculoskelet Disord. (2022) 23:424. doi: 10.1186/s12891-022-05389-4
29. Wu L, Cen C, Yue X, Chen L, Wu H, Yang M, et al. A clinical-radiomics nomogram based on dual-layer spectral detector CT to predict cancer stage in pancreatic ductal adenocarcinoma. Cancer Imaging. (2024) 24:55. doi: 10.1186/s40644-024-00700-z
30. Xu H, Zhu N, Yue Y, Guo Y, Wen Q, Gao L, et al. Spectral CT-based radiomics signature for distinguishing Malignant pulmonary nodules from benign. BMC Cancer. (2023) 23:91. doi: 10.1186/s12885-023-10572-4
31. Al Ajmi E, Forghani B, Reinhold C, Bayat M, Forghani R. Spectral multi-energy CT texture analysis with machine learning for tissue classification: An investigation using classification of benign parotid tumours as a testing paradigm. Eur Radiol. (2018) 28:2604–11. doi: 10.1007/s00330-017-5214-0
32. Zwanenburg A, Vallières M, Abdalah MA, Aerts HJWL, Andrearczyk V, Apte A, et al. The image biomarker standardization initiative: Standardized quantitative radiomics for high-throughput image-based phenotyping. Radiology. (2020) 295:328–38. doi: 10.1148/radiol.2020191145
33. Gao M, Hao Y, Huang MX, Ma DQ, Chen Y, Luo HY, et al. Salivary gland tumours in a northern Chinese population: a 50-year retrospective study of 7190 cases. Int J Oral Maxillofac Surg. (2017) 46:343–9. doi: 10.1016/j.ijom.2016.09.021
34. Mayer M, Alfarra MM, Möllenhoff K, Engels M, Arolt C, Quaas A, et al. The impact of lesion-specific and sampling-related factors on success of salivary gland fine-needle aspiration cytology. Head Neck Pathol. (2025) 19:1. doi: 10.1007/s12105-024-01741-3
35. Li L, Zhao Y, Luo D, Yang L, Hu L, Zhao X, et al. Diagnostic value of single-source dual-energy spectral computed tomography in differentiating parotid gland tumors: Initial results. Quant Imaging Med Surg. (2018) 8:588–96. doi: 10.21037/qims.2018.07.07
36. Wang Y, Hu H, Ban X, Jiang Y, Su Y, Yang L, et al. Evaluation of quantitative dual-energy computed tomography parameters for differentiation of parotid gland tumors. Acad Radiol. (2024) 31:2027–38. doi: 10.1016/j.acra.2023.08.024
37. Sučić M, Ljubić N, Perković L, Ivanović D, Pažanin L, Sučić Radovanović T, et al. Cytopathology and diagnostics of Warthin’s tumour. Cytopathology. (2020) 31:193–207. doi: 10.1111/cyt.12830
38. Matsumoto K, Jinzaki M, Tanami Y, Ueno A, Yamada M, Kuribayashi S. Virtual monochromatic spectral imaging with fast kilovoltage switching: Improved image quality as compared with that obtained with conventional 120-kVp CT. Radiology. (2011) 259:257–62. doi: 10.1148/radiol.11100978
39. Zhang D, Li X, Lv L, Yu J, Yang C, Xiong H, et al. Improving the diagnosis of common parotid tumors via the combination of CT image biomarkers and clinical parameters. BMC Med Imaging. (2020) 20:38. doi: 10.1186/s12880-020-00442-x
40. Jung YJ, Han M, Ha EJ, Choi JW. Differentiation of salivary gland tumors through tumor heterogeneity: A comparison between pleomorphic adenoma and Warthin tumor using CT texture analysis. Neuroradiology. (2020) 62:1451–8. doi: 10.1007/s00234-020-02485-x
41. Gabelloni M, Faggioni L, Attanasio S, Vani V, Goddi A, Colantonio S, et al. Can magnetic resonance radiomics analysis discriminate parotid gland tumors? A pilot study. Diagnostics. (2020) 10:900. doi: 10.3390/diagnostics10110900
42. Speight PM. Update on diagnostic difficulties in lesions of the minor salivary glands. Head Neck Pathol. (2007) 1:55–60. doi: 10.1007/s12105-007-0010-9
43. Greffier J, Villani N, Defez D, Dabli D, Si-Mohamed S. Spectral CT imaging: Technical principles of dual-energy CT and multi-energy photon-counting CT. Diagn Interv Imaging. (2023) 104:167–77. doi: 10.1016/j.diii.2022.11.003
44. So A, Nicolaou S. Spectral computed tomography: Fundamental principles and recent developments. Korean J Radiol. (2021) 22:86–96. doi: 10.3348/kjr.2020.0144
45. Yu Q, Wang A, Gu J, Li Q, Ning Y, Peng J, et al. Multiphasic CT-based radiomics analysis for the differentiation of benign and Malignant parotid tumors. Front Oncol. (2022) 12:913898. doi: 10.3389/fonc.2022.913898
46. Scapicchio C, Gabelloni M, Barucci A, Cioni D, Saba L, Neri E. A deep look into radiomics. Radiol Med. (2021) 126:1296–311. doi: 10.1007/s11547-021-01389-x
47. Hu Z, Guo J, Feng J, Huang Y, Xu H, Zhou Q. Value of T2-weighted-based radiomics model in distinguishing Warthin tumor from pleomorphic adenoma of the parotid. Eur Radiol. (2023) 33:4453–63. doi: 10.1007/s00330-022-09295-0
48. Ganaie MA, Tanveer M, Suganthan PN, Snasel V. Oblique and rotation double random forest. Neural Netw. (2022) 153:496–517. doi: 10.1016/j.neunet.2022.06.012
49. Hu J, Szymczak S. A review on longitudinal data analysis with random forest. Brief Bioinform. (2023) 24:bbad002. doi: 10.1093/bib/bbad002
50. Yan Y, Wang Y, Lei Y. Micro learning support vector machine for pattern classification: A high-speed algorithm. Comput Intell Neurosci. (2022) 2022:4707637. doi: 10.1155/2022/4707637
51. Pushpa B, Baskaran B, Vivekanandan S, Gokul P. Liver fat analysis using optimized support vector machine with support vector regression. Technol Health Care. (2023) 31:867–86. doi: 10.3233/THC-220254
Keywords: parotid tumor, dual-energy computed tomography, radiomics, machine learning, combined nomogram, identification
Citation: Gong Z, Li J, Han Y, Chen S and Wang L (2025) Nomogram combining dual-energy computed tomography features and radiomics for differentiating parotid warthin tumor from pleomorphic adenoma: a retrospective study. Front. Oncol. 15:1505385. doi: 10.3389/fonc.2025.1505385
Received: 02 October 2024; Accepted: 27 January 2025;
Published: 04 March 2025.
Edited by:
Abhishek Mahajan, The Clatterbridge Cancer Centre, United KingdomReviewed by:
Lorenzo Faggioni, University of Pisa, ItalyCopyright © 2025 Gong, Li, Han, Chen and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lijun Wang, d2FuZ2xqMzQ1QDE2My5jb20=
Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.
Research integrity at Frontiers
Learn more about the work of our research integrity team to safeguard the quality of each article we publish.