Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Oncol., 25 September 2025

Sec. Head and Neck Cancer

Volume 15 - 2025 | https://doi.org/10.3389/fonc.2025.1650943

Clinical value of intratumoral and peritumoral CT radiomics models for discriminating benign and malignant parotid gland tumors

Cong Zhang,Cong Zhang1,2Naijing Shi,Naijing Shi1,2Yiru WangYiru Wang2Mohan HaoMohan Hao2Jinwu Ren*Jinwu Ren1*
  • 1Department of Medical Imaging, The First Central Hospital of Baoding, Baoding, Hebei, China
  • 2Graduate School, Chengde Medical University, Chengde, Hebei, China

Objective: To evaluate the utility of combining unenhanced and contrast-enhanced CT intratumoral and peritumoral radiomic features with clinical variables for distinguishing benign from malignant parotid gland tumors.

Methods: We retrospectively collected clinical and imaging data from 171 patients with pathologically confirmed parotid gland tumors treated at Baoding First Central Hospital between June 2019 and June 2025 (101 benign, 70 malignant). Tumor ROIs were manually delineated slice-by-slice on non-contrast, arterial-phase and venous-phase CT images, and peritumoral regions were automatically expanded by 1–4 mm. The cohort was randomly split into training and test sets at a 7:3 ratio. After extraction and selection of radiomic features, multiple models were constructed for intratumoral, various peritumoral ranges (1–4 mm) and intratumoral+peritumoral combinations. Model performance was evaluated by ROC curves, the optimal radiomics model was selected and integrated with the clinical model to produce a combined model, and a nomogram was subsequently developed.

Results: The AUC values of the intratumoral, peritumoral (1–4 mm) and intratumoral+peritumoral models in the training set were 0.966, 0.953, 0.927, 0.983, 0.947, 0.959, 0.956, 0.909 and 0.976, respectively; in the test set the AUCs were 0.797, 0.766, 0.791, 0.714, 0.710, 0.805, 0.836, 0.778 and 0.753, respectively. According to the DeLong test, in the training set the differences between intratumor+peritumor 3mm vs. peritumor 3mm and between intratumor+peritumor 3mm vs. intratumor+peritumor 4mm were statistically significant (p = 0.022 and p = 0.026, respectively); in the test set, differences among the models were not statistically significant (P > 0.05). From this, it can be seen the combined intratumoral + 2 mm peritumoral radiomics model demonstrated superior diagnostic performance compared to models based exclusively on either intratumoral or peritumoral features. Consequently, this model was designated as the optimal radiomic signature and was integrated with independent clinical risk factors—specifically symptomatology and tumor margin status—to construct a combined clinical–radiomics predictive model. In the training and test sets, the AUC values of the radiomics model were 0.956 and 0.836, respectively, while those of the clinical model were 0.774 and 0.703. The combined model achieved AUC values of 0.974 and 0.844, demonstrating significantly superior diagnostic performance compared to the standalone clinical or radiomics models, along with the highest clinical utility. According to the Delong test, in the training set the differences between the clinical model and the combined model, and between the clinical model and the radiomics model, were statistically significant (p = 0.000 and p = 0.000, respectively); in the test set, differences among the models were not statistically significant (P > 0.05).

Conclusion: A multiphase CT radiomics approach that fuses intratumoral features with a 2 mm peritumoral zone robustly distinguishes benign from malignant parotid gland tumors. Integration with key clinical predictors further enhances diagnostic accuracy, supporting clinical translation of the combined model for noninvasive tumor characterization.

1 Introduction

Salivary gland neoplasms represent one of the most common tumor entities within the head and neck region, of which parotid gland tumors (PGTs) account for approximately 80%, with nearly 20% exhibiting malignant behavior (1). Surgical excision remains the gold-standard treatment for PGTs; however, divergent histopathological subtypes necessitate tailored operative techniques and carry distinctly different prognoses (2). While benign parotid tumors generally confer favorable postoperative outcomes, malignant lesions demonstrate aggressive invasion, with a substantial risk of local recurrence and distant metastasis (3). Consequently, accurate preoperative discrimination between benign and malignant parotid lesions is of paramount clinical importance for guiding individualized treatment strategies.

Fine-needle aspiration biopsy (FNAB) is widely regarded as a reliable method for histological characterization of salivary gland masses (4). Nevertheless, FNAB is inherently invasive and may be complicated by hemorrhage, facial nerve injury, or acute inflammatory reactions (5, 6). Thus, noninvasive modalities capable of robustly distinguishing tumor malignancy are highly desirable. Conventional imaging techniques—including ultrasound, magnetic resonance imaging (MRI), and computed tomography (CT)—play pivotal roles in the diagnostic workup. Ultrasonography, as a first-line tool, effectively delineates cystic versus solid components and assesses lesion margins, yet its diagnostic yield is heavily operator-dependent and limited in evaluating deep lobe involvement. MRI offers superior soft-tissue contrast, and diffusion-weighted imaging (DWI) can aid in differentiating benign from malignant histology (7); however, its utility may be constrained by cost and contraindications such as metallic implants. CT reliably evaluates tumor margins, especially in the deep lobe, but demonstrates limited specificity for histopathological subtype discrimination (8). Although CT imaging involves ionizing radiation—especially salient in multiphasic scans and when contrast enhancement requires iodinated contrast agents—CT offers higher spatial resolution compared with MRI, which facilitates clearer delineation of tumor margins and thereby allows more precise definition of tumor extent and surrounding tissues.

In this context, there exists an urgent need for a more precise, noninvasive preoperative assessment method to differentiate benign from malignant parotid gland tumors. Radiomics—a quantitative imaging analysis technique that extracts high-dimensional features from medical images—has emerged as a promising approach to capture tumor microenvironment and heterogeneity, potentially reflecting underlying molecular and cellular characteristics (9). Recent studies have applied radiomics in head and neck malignancies (10), glioblastoma (11), breast cancer (12), and hepatocellular carcinoma (13), demonstrating strong associations between radiomic signatures and tumor histology, grade, and prognosis. Preliminary investigations into parotid tumors have leveraged radiomic features for lesion classification (1416); however, most efforts have focused solely on intratumoral regions, overlooking the peritumoral microenvironment. Evidence suggests that peritumoral tissues may harbor crucial information regarding tumor invasiveness and heterogeneity (17). To date, few studies have systematically evaluated the peritumoral zone in radiomic analyses for parotid lesions, and optimal peritumoral margin size remains undetermined, often selected based on anecdotal experience or extrapolation from other tumor types.

Accordingly, the present study aims to develop and compare intratumoral and peritumoral CT radiomics models—across multiple peritumoral margins—derived from unenhanced and contrast-enhanced scans to distinguish benign from malignant parotid gland tumors. By integrating comprehensive radiomic profiling of both tumor and surrounding tissue, this work seeks to furnish clinicians with a more accurate and holistic noninvasive tool for preoperative decision-making.

2 Materials and methods

2.1 Patient cohort

Clinical and imaging data were retrospectively collected for 171 patients treated at The First Central Hospital of Baoding between June 2019 and June 2025 whose postoperative histopathology confirmed benign parotid gland tumors (BPT) or malignant parotid gland tumors (MPT). Inclusion criteria were: 1.histopathological diagnosis of primary parotid gland neoplasm; 2.availability of complete clinical records and high-quality CT scans, including unenhanced and two‐phase contrast-enhanced examinations; and 3.absence of significant motion or foreign‐body artifacts on CT images. Exclusion criteria comprised:1. prior parotid surgery or radiotherapy/chemotherapy; 2.history of invasive preoperative procedures such as biopsy; and 3. tumor diameter < 1 cm precluding reliable lesion segmentation (Figure 1). The cohort was randomly divided in a 7:3 ratio into a training set (n = 119; BPT = 68, MPT = 51) and a testing set (n = 52; BPT = 33, MPT = 19). The study protocol adhered to the Declaration of Helsinki and was approved by the Institutional Ethics Committee of The First Central Hospital of Baoding, with waiver of informed consent (Approval No. Kuai [2025]018).

Figure 1
Flowchart detailing a study on parotid tumors. It includes clinical and imaging data for 171 patients, with 101 having benign parotid tumors (BPT) and 70 malignant parotid tumors (MPT). Inclusion criteria: confirmed pathology, complete imaging data, and acceptable CT image quality. Exclusion criteria: history of parotid surgery, invasive procedures, or small tumors. The population is divided into a training set (119 patients: 68 BPT, 51 MPT) and a test set (52 patients: 33 BPT, 19 MPT).

Figure 1. Flowchart of patient selection in this study.

2.2 Imaging acquisition

All examinations were performed using a Philips Brilliance i CT 128-slice scanner, covering bilateral parotid regions in the supine position. Scanning parameters were: 120 kV tube voltage, automatic tube current modulation, 1 mm slice thickness with 1 mm interslice interval, pitch 0.984:1, and gantry rotation time 0.6 s. Contrast enhancement was achieved via peripheral intravenous injection of iopamidol ( 300mg I/mL; 80–100 mL) at 3 mL/s using a dual-head power injector. Arterial-phase images were acquired 30–50 s after injection, and venous-phase images at 65–70 s. Patients were instructed to remove dentures and minimize swallowing during image acquisition.

2.3 Imaging and clinical data assessment

CT imaging features were independently evaluated in a blinded fashion by two radiologists (each with > 5 years of head and neck diagnostic experience) unaware of the clinical information and pathological results. In this study, two radiologists independently evaluated the imaging features of all 171 patients. The initial discordance rate between the two readers was approximately 12% (about 21 cases). All discrepancies were resolved successfully through consensus discussion between the two radiologists, and a final unanimous assessment was reached. Imaging features included tumor location (left/right), margin definition (well-defined/ill-defined), shape (regular/irregular), deep-lobe involvement (determined by the imaginary line connecting the most dorsal point of the posterior facial vein to the most dorsal point of the ipsilateral mastoid process), calcification (present/absent), attenuation homogeneity (homogeneous/heterogeneous), lymph node enlargement (present/absent), maximum tumor diameter, peak enhancement phase (arterial/venous), and enhancement uniformity (homogeneous/heterogeneous). Clinical variables comprised sex, age, smoking history, alcohol consumption, and presence of clinical symptoms.

2.4 Image segmentation

Three‐phase CT images were exported from the PACS system and resampled to isotropic voxels of 1 mm × 1 mm × 1 mm; all image intensities were then normalized. A radiologist blinded to the pathological results manually delineated the tumor region of interest (ROI) on axial unenhanced, arterial‐phase, and venous‐phase images using ITK‐SNAP software, carefully excluding bone and vessels. The intratumoral ROI was subsequently expanded outward by 1 mm, 2 mm, 3 mm, and 4 mm to generate peritumoral ROIs (workflow shown in Figure 2); these were labeled Intra, Peri1mm, Peri2mm, Peri3mm, and Peri4mm. ROI consistency was maintained by first outlining the arterial‐phase tumor ROI, then applying the same boundaries to the unenhanced and venous‐phase ROIs, followed by manual adjustment to ensure accuracy. To assess segmentation reproducibility, a subset of 30 patients was randomly selected one month later for repeat delineation, and the intraclass correlation coefficient (ICC) was calculated.

Figure 2
CT scan images labeled A to E show the progression of a highlighted region on the left side in red, with a green outline appearing from panels B to E, indicating a focus area.

Figure 2. Intratumoral ROI (A), intratumoral plus peritumoral 1mm ROI (B), intratumoral plus peritumoral 2mm ROI (C), intratumoral plus peritumoral 3mm ROI (D), and intratumoral plus peritumoral 4mm ROI (E).

2.5 Feature extraction and selection

Radiomic features were extracted using PyRadiomics v3.0.1. Features with ICC ≥ 0.75 were retained and standardized by z-score normalization. Univariate analysis (Student’s t-test or Mann–Whitney U test) and Pearson correlation filtering were performed, followed by least absolute shrinkage and selection operator (LASSO) logistic regression with ten-fold cross-validation to optimize the penalty parameter and select the most robust predictors.

2.6 Model construction

Clinical–imaging predictors significantly associated with BPT and MPT were identified by univariate and multivariate logistic regression to build the clinical model. Given prior evidence that CT-based machine learning using a support vector machine (SVM) classifier achieves high accuracy for parotid tumor discrimination (18), SVM algorithms were employed to develop nine radiomics models: Intra; Peri1mm; Peri2mm; Peri3mm; Peri4mm; Intra + Peri1mm; Intra + Peri2mm; Intra + Peri3mm; and Intra + Peri4mm. Receiver operating characteristic (ROC) curves were plotted and area under the curve (AUC) values calculated to evaluate each model’s performance in differentiating benign from malignant parotid lesions. The optimal radiomics model was then combined with the clinical model to construct a joint model, which was visualized as a nomogram.

2.7 Statistical analysis

Clinical and imaging characteristics of the training and validation cohorts were analyzed using Python 3.7.0 with statsmodels v0.13. Continuous variables were compared by independent‐samples t-test or Mann–Whitney U test, and categorical variables by χ² test; P < 0.05 was considered statistically significant. Diagnostic performance of each model was assessed by ROC analysis, reporting AUC, accuracy, sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV). Delong’s test was used to compare AUCs between models. Calibration performance was evaluated via calibration curves, and clinical utility was appraised using decision curve analysis (DCA). The overall radiomics workflow is depicted in Figure 3.

Figure 3
Flowchart depicting a medical imaging process. It starts with CT imaging acquisition, showing a patient in a scanner and a series of scans. Next is VOI segmentation with highlighted scan regions. Feature selection includes various graphs and data charts. Nomogram plotting has an axis chart. Performance evaluation contains ROC curves and a heatmap for diagnostic accuracy.

Figure 3. Radiomics analysis workflow of this study.

3 Results

3.1 Clinical data and imaging features

This study enrolled 171 patients, of whom 101 were pathologically confirmed as having benign parotid tumors (BPT) and 70 as malignant parotid tumors (MPT); the cohort was randomly split 7:3 into a training set of 119 patients (68 BPT, 51 MPT) and a test set of 52 patients (33 BPT, 19 MPT). The detailed clinical and imaging characteristics are summarized in Table 1. There were no statistically significant differences between BPT and MPT groups in age, sex, smoking history, alcohol history, maximum tumor diameter, shape, location, distribution, calcification, attenuation homogeneity, peak enhancement phase, or enhancement uniformity (all P > 0.05). Furthermore, univariate and multivariate logistic regression analyses identified clinical symptoms and tumor margin status as independent predictors of malignancy (P < 0.05; see Table 2), upon which the clinical model was built, yielding AUCs of 0.774 and 0.703 in the training and test sets, respectively.

Table 1
www.frontiersin.org

Table 1. Baseline characteristics of patients in training and validation set.

Table 2
www.frontiersin.org

Table 2. Univariable and multivariable logistic regression analysis of factors.

3.2 Radiomics model development

A total of 1,834 radiomic features per ROI were extracted using PyRadiomics v3.0.1, comprising 360 first-order features, 14 shape features, and 1,460 texture features (440 gray-level co-occurrence matrix, 280 gray-level dependence matrix, 320 gray-level run-length matrix, 320 gray-level size-zone matrix, and 14 neighborhood gray-tone difference matrix features), resulting in 5,502 features across the three CT phases. After LASSO selection, 20, 33, 39, 34, and 60 features were retained for the Intra, Peri1 mm, Peri2 mm, Peri3 mm, and Peri4 mm ROIs respectively (see Figure 4). In both the training and test cohorts, the combined intratumoral + 2 mm peritumoral model achieved AUCs of 0.956 and 0.836, outperforming all other radiomics models (Table 3). The ROC curves for all nine radiomics models are displayed in Figure 5; Thus, the intratumoral + 2 mm peritumoral model emerged as the optimal radiomics signature, According to the Delong test, in the training set the differences between intratumor+peritumor 3mm vs. peritumor 3mm and between intratumor+peritumor 3mm vs. intratumor+peritumor 4mm were statistically significant (p = 0.022 and p = 0.026, respectively); in the test set, differences among the models were not statistically significant (P > 0.05). Figure 6 illustrates the Delong test results for each model.

Figure 4
Panel A shows a regularization plot with coefficients against Lambda, displaying various colored lines converging toward zero. Panel B illustrates a mean squared error (MSE) plot against Lambda, featuring a U-shaped curve with error bars. Panel C is a horizontal bar chart showing feature names and their coefficient values, with bars extending in both positive and negative directions.

Figure 4. Process of extracting radiomics features. (A) The vertical axis of LASSO regression represents the coefficients corresponding to 5502 features, and the horizontal axis represents the adjustment parameter λ; (B) The optimal λ (λ=0.0168) is selected through mRMR; (C) 56 IntraPeri optimal features have been extracted. IntraPeri, Intratumoral + 2mm peritumoral combined model; LASSO, least absolute shrinkage and selection operator; mRMR, maximum relevance-minimum redundancy.

Table 3
www.frontiersin.org

Table 3. The predictive performance of the nine models.

Figure 5
Two ROC curves labeled A and B compare different SVM models' performance in medical imaging analysis. Curve A shows higher performance with Intra_SVM achieving an AUC of 0.966, while B shows lower overall performance with Intra_SVM at an AUC of 0.797. The legend lists various model specifications and corresponding AUC values, highlighting performance variations.

Figure 5. ROC curves of nine radiomic models (A) training set; (B) test set.

Figure 6
Two colorful heatmaps labeled “Cohort train Delong” and “Cohort test Delong” compare support vector machine (SVM) models with different peri- and intra-peritumoral measurements. A color scale from blue to red indicates values 0.0 to 0.9.

Figure 6. Summary of DeLong test results for each model.

3.3 Nomogram construction

The optimal radiomics model (Intra + Peri2 mm) was combined with the clinical model to generate a nomogram (joint model) as detailed in Table 4 and Figure 7. In the training and test sets, the radiomics model AUCs were 0.956 and 0.836, respectively, while the clinical model yielded AUCs of 0.774 and 0.703. The combined model demonstrated AUC values of 0.974 and 0.844, surpassing those of the individual radiomics and clinical models, and indicating superior diagnostic performance. According to the Delong test, in the training set the differences between the clinical model and the combined model, and between the clinical model and the radiomics model, were statistically significant (p = 0.000 and p = 0.000, respectively); in the test set, differences among the models were not statistically significant (P > 0.05). Delong’s comparisons among the three models are shown in Figure 8. The nomogram assigns substantial weight to the radiomics score in malignancy risk estimation (Figure 9). Calibration curves for all three models in both cohorts confirmed excellent goodness-of-fit (Figure 10), with Hosmer–Lemeshow P values of 0.392 and 0.435 for the training and test sets, respectively, indicating no significant deviation from perfect calibration. Decision curve analysis (Figure 11) demonstrated that the joint model achieved the highest net clinical benefit for parotid tumor discrimination in both cohorts.

Table 4
www.frontiersin.org

Table 4. Diagnostic efficacy of clinical, radiomics, and nomogram models in both training and testing cohorts.

Figure 7
Two ROC curve charts labeled A and B. Chart A is titled “Cohort train ROC” and chart B is titled “Cohort test ROC”. Both charts plot sensitivity against 1-specificity. Chart A shows three lines: clinic (magenta), radiomic (dotted black), and combined (cyan), with AUCs of 0.774, 0.956, and 0.974 respectively. Chart B shows similar lines with AUCs of 0.703, 0.836, and 0.844. The diagonal line represents random classification.

Figure 7. ROC curves demonstrating the discriminatory performance of radiomics, clinical, and integrated models in parotid gland tumor classification (A) training set; (B) test set.

Figure 8
Heatmaps labeled “Cohort train Delong” and “Cohort test Delong.” Each features squares representing Radiomic and Clinic data. Values are shown in colors from red to blue on a scale indicating intensity. Image A shows Radiomic (0.000) and Combined (0.000) for Clinic, Radiomic (0.163); Image B shows Radiomic (0.125) and Combined (0.077) for Clinic, Radiomic (0.692).

Figure 8. Delong test results comparing the performance of clinical, radiomics, and combined models.

Figure 9
Nomogram with five sections labeled Points, Symptom, Margin, Radiomic, Total Points, and Risk. Each section has a numerical scale, allowing for scoring and risk assessment based on input values.

Figure 9. Nomogram of the combined model.

Figure 10
Two calibration plots labeled A and B. Plot A shows cohort train calibration, displaying mean predicted probability versus fraction of positives. Plot B shows cohort test calibration with similar axes. Both plots include lines for clinic, radiomic, and combined data compared to a perfectly calibrated line. Clinic data is in blue, radiomic in orange, and combined in green. Dotted lines indicate perfect calibration.

Figure 10. calibration curves of clinical, radiomics, and nomogram models (A) training set; (B) test set.

Figure 11
Two decision curve analysis graphs comparing net benefit against threshold probability. Graph A, labeled “Cohort train DCA,” shows clinic, radiomic, and combined strategies alongside “Treat all” and “Treat none” baselines. Graph B, labeled “Cohort test DCA,” presents similar data for the test cohort with respective color-coded lines.

Figure 11. DCA curves of clinical, radiomics, and nomogram models (A) training set; (B) test set.

4 Discussion

Benign and malignant parotid gland tumors exhibit markedly distinct biological behaviors, leading to divergent treatment strategies and prognosis. For benign lesions, superficial or partial parotidectomy is typically preferred; by contrast, malignant tumors are associated with higher rates of recurrence and metastasis and often warrant more extensive resection, such as total or radical parotidectomy (19). Consequently, accurate preoperative differentiation between benign and malignant tumors is critical for optimizing therapeutic decision-making and patient outcomes. Conventional imaging assessments rely heavily on the radiologist’s expertise and subjective interpretation, and overlapping radiographic features among different histological subtypes further limit diagnostic precision (20). As a result, the accuracy of preoperative characterization of parotid tumors remains suboptimal (21).

Radiomics—an emerging artificial intelligence–driven methodology—enables noninvasive extraction of high-dimensional quantitative features from routine medical images, thereby revealing intrinsic tumor biology (22). Prior studies have predominantly focused on intratumoral radiomic signatures (23, 24), with comparatively little attention paid to the peritumoral microenvironment. Xu et al. demonstrated the utility of tumor-based radiomic features for distinguishing benign from malignant parotid lesions (25). Zheng et al. constructed an MRI-based radiomics nomogram for preoperative differentiation of benign and malignant parotid tumors. By extracting texture features from T1WI and fat-suppressed T2WI and integrating clinical factors (such as deep-lobe invasion and peritumoral tissue infiltration), the model demonstrated high discriminative performance in both training and validation cohorts (AUC > 0.93). The pathological basis is that malignant tumors often demonstrate infiltrative growth and high tissue heterogeneity—attributes that can be quantitatively captured by radiomic features—thus surpassing visual assessment and aiding clinical decision-making (26). Faggioni et al. demonstrated that MRI-based radiomic features—particularly T2-weighted Skewness and pcsT1-weighted GLCM_InverseVariance—can effectively discriminate pleomorphic adenoma (PA) from Warthin tumor (WT), exhibiting high specificity. One study found that on contrast-enhanced T1-weighted images, GLCM_InverseVariance was the most discriminative feature (AUC=0.90), reflecting the relatively homogeneous enhancement pattern of WT; on T2-weighted images, skewness had the highest specificity (88%), with high skewness values associated with intralesional cystic change, necrosis, or hypervascularity typical of WT. WT commonly contains cystic components and lymphoid stroma with complex architecture, whereas PA is relatively more homogeneous; these histopathologic differences are quantitatively captured by radiomic features, enabling noninvasive and accurate preoperative discrimination (27). Evidence suggests that the peritumoral zone harbors critical biological information—such as angiogenesis, lymphovascular invasion, and stromal reactions—that may drive tumor progression (28, 29). These pathophysiological changes can likewise be captured through radiomic analysis. Indeed, peritumoral radiomics has shown promise in predicting non–small cell lung cancer subtypes (30), epithelial ovarian cancer extra-pelvic metastasis (31), and malignant brain tumor classification (32). However, its application to parotid gland tumors remains underexplored.

In this study, the diagnostic efficacy of intratumoral and peritumoral radiomic features for discriminating benign and malignant parotid gland tumors was systematically evaluated, revealing that the peritumoral microenvironment may harbor critical biological cues for tumor characterization. Likewise, the selection of peritumoral margin size significantly influences radiomic predictive performance (33), and the systematic identification of an optimal expansion zone can further enhance model accuracy. Although optimal peritumoral ranges vary by anatomical site, definitive evidence for parotid tumors has been lacking. By constructing models at 1 mm, 2 mm, 3 mm, and 4 mm peritumoral expansions, we sought to determine the ideal margin for parotid neoplasms and confirmed that peritumoral heterogeneity provides substantial complementary value in malignancy prediction. We integrated multiphase CT features—unenhanced, arterial, and venous—into nine distinct radiomics models and found the 2 mm expansion to be optimal, yielding AUCs of 0.927 and 0.791 in the training and test cohorts, respectively, outperforming the 1 mm, 3 mm, and 4 mm models. Furthermore, the combined intratumoral + 2 mm peritumoral model demonstrated superior discrimination (AUC = 0.956 in training, 0.836 in testing), with high accuracy and specificity. Consistent with Shen et al.’s findings that a 2 mm peritumoral margin maximizes AUC when compared to 5 mm (34), our results corroborate the primacy of the 2 mm zone. Prior research indicates that support vector machine (SVM)–based classifiers outperform other machine learning algorithms (18) and, owing to their robust generalizability, suitability for small datasets, and interpretability, we employed SVM throughout. Moreover, the Peri1 mm model also demonstrated robust diagnostic performance; likewise, the Peri3 mm (AUC = 0.983) and Peri4 mm (AUC = 0.976) models achieved high AUCs in the training cohort, yet their discriminative accuracy declined markedly in the test cohort. Moreover, we observed that standalone peritumoral models underperformed relative to intratumoral models, potentially reflecting our moderately sized cohort and the predominance of larger tumors, which inherently contain abundant discriminative information, thereby limiting the incremental value of peritumoral features.

The Intra + Peri2 mm model demonstrated a pivotal role in enhancing the predictive performance for parotid gland tumor characterization. By integrating intratumoral and 2 mm peritumoral radiomic signatures, the IntraPeri2 mm model achieved AUCs of 0.956 and 0.836 in the training and test cohorts, respectively. Consequently, we infer that intratumoral and peritumoral features are complementary, jointly conferring superior accuracy in distinguishing benign from malignant parotid lesions. Decision curve analysis further confirmed that the IntraPeri2 mm model yields a meaningful net clinical benefit.

Furthermore, this study identified clinical symptoms and tumor margin status as independent predictors of malignancy, and combined these clinical‐imaging variables with the optimal radiomic signature to construct a joint model. The resulting combined model attained AUCs of 0.974 in the training set and 0.844 in the test set—outperforming both the clinical‐only and radiomics‐only models—and was visualized via a nomogram. The integration of these clinical and radiomic parameters enables more accurate preoperative stratification of parotid gland tumors, thereby holding significant implications for individualized clinical decision-making.

5 Study limitations

This investigation is subject to several constraints. First, its retrospective design inherently predisposes it to selection bias. Second, as a single-center study with a modest and imbalanced sample, external validation was not feasible; the limited cohort may constrain the robustness and generalizability of our findings, underscoring the need for prospective, multicenter trials with larger, more heterogeneous populations. Third, manual delineation of tumor ROIs introduces intra- and interobserver variability; the implementation and validation of automated or semi-automated segmentation algorithms are warranted to enhance reproducibility and mitigate observer-dependent discrepancies.

6 Conclusion

Multiphase CT-derived radiomic signatures encompassing both intratumoral and peritumoral regions demonstrate high discriminatory power for benign versus malignant parotid gland tumors. Integration of these quantitative imaging biomarkers with independent clinical predictors yields a composite model with superior diagnostic accuracy and clinical net benefit. This noninvasive framework offers a promising avenue for precision stratification and individualized management of parotid neoplasms.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding author.

Ethics statement

The studies involving humans were approved by Ethics Committee of Baoding First Central Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.

Author contributions

CZ: Data curation, Methodology, Investigation, Validation, Writing – review & editing, Software, Writing – original draft. NS: Writing – review & editing, Writing – original draft, Data curation. YW: Validation, Supervision, Writing – original draft, Writing – review & editing. MH: Writing – original draft, Writing – review & editing, Data curation. JR: Writing – review & editing, Writing – original draft.

Funding

The author(s) declare that no financial support was received for the research and/or publication of this article.

Acknowledgments

The authors gratefully acknowledge the expert contributions of the Radiology, Oral Surgery, and Pathology Department staff at The First Central Hospital of Baoding for their invaluable clinical insights and technical assistance.

Conflict of interest

The 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.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

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.

References

1. Zhan KY, Khaja SF, Flack AB, and Day TA. Benign parotid tumors. Otolaryngol Clin North Am. (2016) 49:327–42. doi: 10.1016/j.otc.2015.10.005

PubMed Abstract | Crossref Full Text | Google Scholar

2. Stoia S, Băciuț G, Lenghel M, Badea R, Csutak C, Rusu GM, et al. Cross-sectional imaging and cytologic investigations in the preoperative diagnosis of parotid gland tumors – an updated literature review. Bosn J Basic Med Sci. (2020) 21:19–32. doi: 10.17305/bjbms.2020.5028

PubMed Abstract | Crossref Full Text | Google Scholar

3. Maahs GS, Oppermann PDO, Maahs LGP, MaChado Filho G, and Ronchi AD. Parotid gland tumors: a retrospective study of 154 patients. Braz J Otorhinolaryngol. (2015) 81:301–6. doi: 10.1016/j.bjorl.2015.03.007

PubMed Abstract | Crossref Full Text | Google Scholar

4. Haldar S, Sinnott JD, Tekeli KM, Turner SS, and Howlett DC. Biopsy of parotid masses: review of current techniques. World J Radiol. (2016) 8:501. doi: 10.4329/wjr.v8.i5.501

PubMed Abstract | Crossref Full Text | Google Scholar

5. Bahar G, Dudkiewicz M, Feinmesser R, Joshua B, Braslavsky D, Popovtzer A, et al. Acute parotitis as a complication of fine-needle aspiration in Warthin’s tumor. A unique finding of a 3-year experience with parotid tumor aspiration. Otolaryngol Neck Surg. (2006) 134:646–9. doi: 10.1016/j.otohns.2005.10.050

PubMed Abstract | Crossref Full Text | Google Scholar

6. Kechagias N, Ntomouchtsis A, Valeri R, Patrikidou A, Kitikidou K, Xirou P, et al. Fine-needle aspiration cytology of salivary gland tumours: a 10-year retrospective analysis. Oral Maxillofac Surg. (2012) 16:35–40. doi: 10.1007/s10006-011-0291-8

PubMed Abstract | Crossref Full Text | Google Scholar

7. Sun Q, Ma C, Dong M, Jiang M, and Tao X. Effects of region of interest sizes on apparent diffusion coefficient measurements of pleomorphic adenoma, Warthin tumor, and normal parotid parenchyma. Quant Imaging Med Surg. (2019) 9:681–90. doi: 10.21037/qims.2019.04.11

PubMed Abstract | Crossref Full Text | Google Scholar

8. Dong Y, Lei G, Wang S, Zheng S, Ge Y, and Wei F. Diagnostic value of CT perfusion imaging for parotid neoplasms. Dentomaxillofac Radiol. (2014) 43:20130237. doi: 10.1259/dmfr.20130237

PubMed Abstract | Crossref Full Text | Google Scholar

9. Lambin P, Leijenaar RTH, Deist TM, Peerlings J, De Jong EEC, Van Timmeren J, et al. Radiomics: the bridge between medical imaging and personalized medicine. Nat Rev Clin Oncol. (2017) 14:749–62. doi: 10.1038/nrclinonc.2017.141

PubMed Abstract | Crossref Full Text | Google Scholar

10. Wu W, Ye J, Wang Q, Luo J, and Xu S. CT-based radiomics signature for the preoperative discrimination between head and neck squamous cell carcinoma grades. Front Oncol. (2019) 9:821. doi: 10.3389/fonc.2019.00821

PubMed Abstract | Crossref Full Text | Google Scholar

11. Chaddad A, Kucharczyk MJ, Daniel P, Sabri S, Jean-Claude BJ, Niazi T, et al. Radiomics in glioblastoma: current status and challenges facing clinical implementation. Front Oncol. (2019) 9:374. doi: 10.3389/fonc.2019.00374

PubMed Abstract | Crossref Full Text | Google Scholar

12. Xie T, Wang Z, Zhao Q, Bai Q, Zhou X, Gu Y, et al. Machine learning-based analysis of MR multiparametric radiomics for the subtype classification of breast cancer. Front Oncol. (2019) 9:505. doi: 10.3389/fonc.2019.00505

PubMed Abstract | Crossref Full Text | Google Scholar

13. Wu C, Chen J, Fan Y, Zhao M, He X, Wei Y, et al. Nomogram based on CT radiomics features combined with clinical factors to predict Ki-67 expression in hepatocellular carcinoma. Front Oncol. (2022) 12:943942. doi: 10.3389/fonc.2022.943942

PubMed Abstract | Crossref Full Text | Google Scholar

14. Zheng Y, Xu W, Hao D, Liu X, Gao C, Tang G, 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

PubMed Abstract | Crossref Full Text | Google Scholar

15. 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

PubMed Abstract | Crossref Full Text | Google Scholar

16. 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

PubMed Abstract | Crossref Full Text | Google Scholar

17. Jiang X, Wang J, Deng X, Xiong F, Zhang S, Gong Z, et al. The role of microenvironment in tumor angiogenesis. J Exp Clin Cancer Res. (2020) 39:204. doi: 10.1186/s13046-020-01709-5

PubMed Abstract | Crossref Full Text | Google Scholar

18. Zheng Y, Zhou D, Liu H, and Wen M. CT-based radiomics analysis of different machine learning models for differentiating benign and Malignant parotid tumors. Eur Radiol. (2022) 32:6953–64. doi: 10.1007/s00330-022-08830-3

PubMed Abstract | Crossref Full Text | Google Scholar

19. Moore MG, Yueh B, Lin DT, Bradford CR, Smith RV, and Khariwala SS. Controversies in the workup and surgical management of parotid neoplasms. Otolaryngol Neck Surg. (2021) 164:27–36. doi: 10.1177/0194599820932512

PubMed Abstract | Crossref Full Text | Google Scholar

20. Gökçe E. Multiparametric magnetic resonance imaging for the diagnosis and differential diagnosis of parotid gland tumors. J Magn Reson Imaging. (2020) 52:11–32. doi: 10.1002/jmri.27061

PubMed Abstract | Crossref Full Text | Google Scholar

21. Vogl TJ, Albrecht MH, Nour-Eldin N-EA, Ackermann H, Maataoui A, Stöver T, et al. Assessment of salivary gland tumors using MRI and CT: impact of experience on diagnostic accuracy. Radiol Med (Torino). (2018) 123:105–16. doi: 10.1007/s11547-017-0813-z

PubMed Abstract | Crossref Full Text | Google Scholar

22. Wei W, Liu Z, Rong Y, Zhou B, Bai Y, Wei W, et al. A computed tomography-based radiomic prognostic marker of advanced high-grade serous ovarian cancer recurrence: a multicenter study. Front Oncol. (2019) 9:255. doi: 10.3389/fonc.2019.00255

PubMed Abstract | Crossref Full Text | Google Scholar

23. Qi J, Gao A, Ma X, Song Y, Zhao G, Bai J, et al. Differentiation of benign from Malignant parotid gland tumors using conventional MRI based on radiomics nomogram. Front Oncol. (2022) 12:937050. doi: 10.3389/fonc.2022.937050

PubMed Abstract | Crossref Full Text | Google Scholar

24. Piludu F, Marzi S, Ravanelli M, Pellini R, Covello R, Terrenato I, et al. MRI-based radiomics to differentiate between benign and Malignant parotid tumors with external validation. Front Oncol. (2021) 11:656918. doi: 10.3389/fonc.2021.656918

PubMed Abstract | Crossref Full Text | Google Scholar

25. Xu Y, Shu Z, Song G, Liu Y, Pang P, Wen X, et al. The role of preoperative computed tomography radiomics in distinguishing benign and Malignant tumors of the parotid gland. Front Oncol. (2021) 11:634452. doi: 10.3389/fonc.2021.634452

PubMed Abstract | Crossref Full Text | Google Scholar

26. Zheng Y, Li J, Liu S, Cui J, Zhan J, Pang J, et al. MRI-based radiomics nomogram for differentiation of benign and Malignant lesions of the parotid gland. Eur Radiol. (2021) 31:4042–52. doi: 10.1007/s00330-020-07483-4

PubMed Abstract | Crossref Full Text | Google Scholar

27. Faggioni L, Gabelloni M, De Vietro F, Frey J, Mendola V, Cavallero D, et al. Usefulness of MRI-based radiomic features for distinguishing Warthin tumor from pleomorphic adenoma: performance assessment using T2-weighted and post-contrast T1-weighted MR images. Eur J Radiol Open. (2022) 9:100429. doi: 10.1016/j.ejro.2022.100429

PubMed Abstract | Crossref Full Text | Google Scholar

28. Zhang X, Zhang G, Qiu X, Yin J, Tan W, Yin X, et al. Radiomics under 2D regions, 3D regions, and peritumoral regions reveal tumor heterogeneity in non-small cell lung cancer: a multicenter study. Radiol Med (Torino). (2023) 128:1079–92. doi: 10.1007/s11547-023-01676-9

PubMed Abstract | Crossref Full Text | Google Scholar

29. Kolesnikoff N, Chen C-H, and Samuel MS. Interrelationships between the extracellular matrix and the immune microenvironment that govern epithelial tumour progression. Clin Sci. (2022) 136:361–77. doi: 10.1042/CS20210679

PubMed Abstract | Crossref Full Text | Google Scholar

30. Tang X, Huang H, Du P, Wang L, Yin H, and Xu X. Intratumoral and peritumoral CT-based radiomics strategy reveals distinct subtypes of non-small-cell lung cancer. J Cancer Res Clin Oncol. (2022) 148:2247–60. doi: 10.1007/s00432-022-04015-z

PubMed Abstract | Crossref Full Text | Google Scholar

31. Wang X, Wei M, Chen Y, Jia J, Zhang Y, Dai Y, et al. Intratumoral and peritumoral MRI-based radiomics for predicting extrapelvic peritoneal metastasis in epithelial ovarian cancer. Insights Imaging. (2024) 15:281. doi: 10.1186/s13244-024-01855-w

PubMed Abstract | Crossref Full Text | Google Scholar

32. Liu D, Chen J, Ge H, Hu X, Yang K, Liu Y, et al. Differentiation of Malignant brain tumor types using intratumoral and peritumoral radiomic features. Front Oncol. (2022) 12:848846. doi: 10.3389/fonc.2022.848846

PubMed Abstract | Crossref Full Text | Google Scholar

33. Wu L, Gao C, Ye J, Tao J, Wang N, Pang P, et al. The value of various peritumoral radiomic features in differentiating the invasiveness of adenocarcinoma manifesting as ground-glass nodules. Eur Radiol. (2021) 31:9030–7. doi: 10.1007/s00330-021-07948-0

PubMed Abstract | Crossref Full Text | Google Scholar

34. Shen Q. Preoperative CT-based intra- and peri-tumoral radiomic models for differentiating benign and Malignant tumors of the parotid gland: a two-center study. Am J Cancer Res. (2024) 14:4445–58. doi: 10.62347/AXQW1100

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: parotid gland tumor, radiomics, machine learning, intratumoral, peritumoral

Citation: Zhang C, Shi N, Wang Y, Hao M and Ren J (2025) Clinical value of intratumoral and peritumoral CT radiomics models for discriminating benign and malignant parotid gland tumors. Front. Oncol. 15:1650943. doi: 10.3389/fonc.2025.1650943

Received: 20 June 2025; Accepted: 09 September 2025;
Published: 25 September 2025.

Edited by:

Alberto Rodriguez-Archilla, University of Granada, Spain

Reviewed by:

Lorenzo Faggioni, University of Pisa, Italy
Mohamed Mekky, Oral and Maxillofacial Department Faculty of Dentistry Kafrelsheik University, Egypt

Copyright © 2025 Zhang, Shi, Wang, Hao and Ren. 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: Jinwu Ren, cmp3emhhbmdsaWhvbmdAMTYzLmNvbQ==

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.