- 1Department of CT Diagnosis, Second Affiliated Hospital of Harbin Medical University, Harbin, China
- 2Department of Radiology, Heilongjiang Red Cross Sengong General Hospital, Harbin, China
- 3Department of Ultrasound, Second Affiliated Hospital of Harbin Medical University, Harbin, China
Purpose: To determine the incremental diagnostic value of habitat radiomics for risk stratification of thymic epithelial tumors (TETs) based on contrast-enhanced CT (CECT).
Methods: This retrospective study included 220 patients with pathologically confirmed TETs (82 high-risk [B2/B3/thymic carcinoma] and 138 low-risk [A/AB/B1]) who underwent preoperative CECT. Tumors were segmented into 3 subregions (habitats) using k-means clustering, and radiomic features were extracted from both whole-tumor and subregions. After feature selection (variance threshold, reproducibility evaluation, XGBoost-based importance ranking, and recursive feature elimination), three machine learning models were developed (1): a conventional radiomics model (2), a habitat radiomics model, and (3) a combined model integrating both feature sets. Model performance was evaluated using ROC analysis, net reclassification improvement (NRI), integrated discrimination improvement (IDI), calibration metrics, and decision curve analysis (DCA).
Results: The combined model demonstrated superior discrimination (AUC: 0.900) compared to the conventional (AUC: 0.819) and habitat (AUC: 0.734) radiomics models in the independent test set. Although DeLong’s test showed no statistically significant difference (p=0.161), the performance of combined model demonstrated incremental diagnostic value (NRI: 0.286; IDI: 0.209). Calibration and DCA confirmed its robustness and higher net benefit across decision thresholds. While the models’ training performance might suggest overfitting, their test results demonstrate generalizability.
Conclusions: The habitat radiomics approach enables accurate risk stratification prediction in TETs and demonstrates potential as a clinically valuable tool to augment the performance of conventional radiomics models in routine practice.
Introduction
Thymic epithelial tumors (TETs) represented the most prevalent primary neoplasms originating in the anterior mediastinum (1–3). According to the 2021 WHO classification, TETs were histologically stratified into six distinct subtypes: type A, AB, B1, B2, B3, and thymic carcinoma (TC) (4). The histopathological classification was one of the important prognostic factors, particularly in guiding postoperative therapeutic decision-making within the context of neoadjuvant treatment strategies (5). Clinically, these subtypes exhibited markedly divergent biological behaviors: types A, AB, and B1 TETs were associated with favorable prognoses and significantly lower recurrence rates, whereas types B2, B3, and TC demonstrated aggressive clinical courses with substantially elevated recurrence risks (6–8). TETs could be further classified into two prognostic groups: low-risk TETs (LRT, encompassing A, AB, and B1) and high-risk TETs (HRT, including B2, B3, and TC) (9–11). Reliable noninvasive risk stratification prior to therapeutic intervention was critical for prognostication and treatment optimization.
Among imaging modalities, contrast-enhanced computed tomography (CECT) was the first choice for evaluating TETs due to its cost-effectiveness, widespread accessibility, and capacity to delineate detailed tumor morphology and vascular enhancement patterns (12, 13). Previous studies had demonstrated the utility of CECT features in TETs risk stratification (14–17). However, models relying exclusively on conventional radiological signs exhibited limited predictive performance in discriminating high-risk from low-risk subtypes (18). Numerous studies had demonstrated the efficacy of CECT-based radiomics in accurately differentiating LRT from HRT (19–21). Nevertheless, most prior conventional radiomics investigations had predominantly focused on whole-tumor analysis as a single region of interest (ROI), with insufficient attention paid to tumor subregions exhibiting heterogeneous (22, 23). Recently, habitat imaging had emerged as a novel paradigm that specifically identifies intratumoral heterogeneous regions or cellular subpopulations (24–26). By analyzing tumor subregions, this approach enabled more precise characterization of spatial heterogeneity and enhanced delineation of intrinsic biological features, thereby holding significant promise for refining TETs risk stratification. To date, no studies had investigated the potential of habitat-based radiomics for the risk stratification of TETs. Therefore, we implemented habitat-specific radiomics analysis using CECT imaging for TETs risk stratification, quantitatively assessing its complementary benefits to standard whole-tumor radiomics.
Materials and methods
Patients
This study was approved by the institutional review board of the hospital with a waiver for informed (KY2025-157). Figure 1 outlined the study workflow. This retrospective study consecutively enrolled patients with pathologically confirmed TETs at the Affiliated Second Hospital of Harbin Medical University between December 2017 and December 2024, with Figure 2 illustrating the workflow of patient inclusion. The inclusion criteria were (1): histologically confirmed WHO classification of TETs (2); preoperative CECT performed within 1 month prior to surgery (3); no prior oncologic therapy (chemotherapy/radiotherapy). The exclusion criteria were (1): recurrent lesion (2); poor-quality image (motion artifacts or incomplete coverage) (3); undetermined WHO classification. The TETs were classified into LRT (A, AB, and B1) and HRT (B2, B3, and TC). Patients with TETs were randomly divided into training (80%) and test sets (20%) using stratified sampling while preserving the original WHO subtype distribution.
Image acquisition
Three CT scanners (GE Discovery CT750HD; GE LightSpeed VCT; PHILIPS iCT256) were used to perform CECT. The scanning parameters were listed in Supplementary Material 1. The contrast material about 80–100 mL was intravenously administered at a rate of 2.5 mL/s. All patients underwent supine-positioned CT examinations spanning from the pulmonary apex to the infradiaphragmatic region. Arterial phase was obtained by threshold trigger and the venous phase scanning was performed after 30s. Image reconstruction was consistently performed using a mediastinum-optimized convolution kernel to ensure accurate interpretation of mediastinal anatomy. The axial images in venous phase were exclusively selected for subsequent quantitative analysis to avoid potential influences from superior vena cava artifacts that were commonly observed in the arterial phase.
Segmentation
ROIs were manually delineated on axial slices (3D volume) through consensus-based delineation by two radiologists (each with >5 years of expertise in thoracic imaging) using ITK-SNAP software (version 4.2.2, www.itksnap.org). To ensure unbiased evaluation, both radiologists were blinded to all pathological findings and clinical outcomes.
To mitigate outlier effects, intensity values were normalized by truncating the histogram extremities (0.5th to 99.5th percentiles). All images were first spatially normalized to align with a common coordinate space, then isotropically resampled to 1×1×1 mm3 voxel dimensions to ensure uniform spatial resolution.
The tumor heterogeneity was quantitatively assessed by partitioning the images into distinct subregions (habitats) using the k-means clustering algorithm. The elbow method, based on the within-cluster sum of squared errors (SSE), was employed to identify the most plausible segmentation. To determine the optimal number of clusters, we evaluated values of k ranging from 2 to 10. The elbow point was determined by locating the value of k where the rate of decrease in within-cluster SSE sharply diminishes.
Feature extraction
To strictly prevent data leakage, the test set was reserved solely for independent model evaluation and never involved in any prior steps, including feature selection or model development. Radiomic feature extraction was performed on both the whole-tumor ROI and its subregions using PyRadiomics (version 3.1.0, https://pyradiomics.readthedocs.io/en/latest/), encompassing shape-based, first-order statistical, and texture-based features with all available filters applied. Feature definitions and computational methods were listed in the PyRadiomics documentation. For each patient, a total of 6,752 radiomic features were obtained, consisting of 1,688 features from the whole-tumor ROI and an additional 5,064 features (1,688 × 3) derived from three distinct subregions. The extracted features from each subregion were labeled with their corresponding cluster number. Feature normalization was subsequently performed using z-score standardization, where the mean and standard deviation (SD) calculated exclusively from the training set were applied to both training and test sets.
Feature selection
The feature selection pipeline was implemented in four sequential steps (1): Initial variance thresholding eliminated non-informative features (zero variance across the training set) (2); Reproducibility assessment using a randomly selected subset of 30 patients by two radiologists independently delineating ROIs, retaining only features demonstrating good reproducibility (ICC > 0.75) (3); Feature importance was assessed using an XGBoost classifier, with the top 300 most discriminative features selected according to the importance (4); Final selection through recursive feature elimination with 5-fold cross-validation (RFECV) using a XGBoost classifier to determine the minimal optimal feature subset maximizing AUC. To mitigate the class imbalance between the two groups, the support vector machine synthetic minority over-sampling technique (SVMSMOTE) was employed to augment the sample size of high-risk TETs through synthetic data generation.
Model development
The conventional radiomics feature and habitat radiomics feature were independently selected through the established pipeline. Subsequently, two machine learning models were developed using the XGBoost algorithm on the SVMSMOTE-balanced training set (1): a conventional radiomics model incorporating conventional imaging features, and (2) a habitat radiomics model capturing tumor subregional heterogeneity patterns. Following the same selection procedure, a combined model was constructed by using the selected features extracted from the whole tumor and subregions.
Statistics
Sex and group were reported as counts (percentages), while age was summarized as mean ± SD. The difference between training set and test set were assessed using chi-square test (or Fisher’s exact test) for categorical variables and independent t-test (or Mann-Whitney U test) for continuous variables. Model performance was evaluated through receiver operating characteristic (ROC) analysis with calculation of AUC, complemented by metrics including accuracy, sensitivity, and specificity. The incremental prognostic value of habitat features was quantitatively assessed using net reclassification improvement (NRI), integrated discrimination improvement (IDI), and DeLong’s test for AUC comparison. Calibration was verified through Brier score, Brier skill score and calibration curve analysis, while clinical utility was appraised via decision curve analysis (DCA). To enhance interpretability, SHAP (SHapley Additive exPlanations) values were computed to determine feature importance and explain model decisions. All statistical analyses were conducted in Python 3.11. A two-tailed p-values <0.05 considered statistically significant.
Results
Baseline characteristics of patients
Of 229 patients, 9 were excluded according to exclusion criteria. The workflow was shown in Figure 2. The eligible 220 patients with TETs, randomly divided into training (n=176, 80%) and test (n=44, 20%) sets. According to WHO-based risk stratification, 82 (37.3%) were HRT and 138 (62.7%) LRT. Detailed baseline characteristics are summarized in Table 1. There were no statistically significant differences (p>0.05) between the training and test sets in terms of age, gender, group and WHO classification.
Feature selection
Using the elbow method (Figure 3), we identified k = 3 as the optimal number of tumor subregions. From the initial 6,752 extracted features, 6,164 demonstrated non-zero variance, with subsequent reproducibility testing (ICC>0.75) retaining 1,946 robust features (1,429 whole-tumor and 517 subregion-derived). All ICC values were listed in Supplementary Material 2. RFECV distilled these to a 16-feature signature selected from the top 300 features for combined model, whose importance is quantified in Figure 4. The features developed (n=28) conventional and habitat (n=40) radiomics models were listed in Supplementary Material 3. Model interpretability was enhanced through SHAP value heatmap plots Figure 5, which delineate individualized feature impacts across the training and test sets.
Evaluation
The comparative performance analysis of the three predictive models was listed in Table 2. All three models attained an AUC of 1.000, along with accuracy, sensitivity, and specificity all reaching 1.000 on the training set. In the test set, the combined model achieved best performance compared to other models (AUC: 0.819, 0.734 vs 0.900; accuracy: 0.750, 0.705 vs 0.864). Although the models’ performance on the training set could indicate potential overfitting, their strong and consistent results on the independent test set suggest reasonable generalizability. The ROC curve was shown in Figure 6a. The comparison via DeLong’s test revealed no statistically significant difference between conventional radiomics and combined models (p=0.161). However, the combined model showed clinically meaningful incremental value, evidenced by NRI (0.286) and IDI (0.209).

Figure 6. Performance comparison in test set. ROC curves (a), calibration curve (b), and DCA curves (c).
Calibration curves (Figure 6b) revealed that the combined model show favorable agreement and had best alignment across the full probability range compared to other models. Quantitative evaluation revealed that the combined model achieved a significantly lower Brier score (0.121) compared to the conventional (0.193) and habitat radiomics models (0.255), corresponding to a Brier skill score of 37.2% relative to the conventional radiomics model as the reference. DCA (Figure 6c) demonstrated superior clinical utility of the combined model, which provided a higher net benefit across a wide range of threshold probabilities compared to other models.
Discussion
This study developed and validated a CECT-based radiomics model for the risk stratification of TETs, demonstrating favorable predictive performance. Radiomic features extracted from tumor habitat heterogeneity could improve the diagnostic accuracy compared to conventional radiomics models. Moreover, SHAP value analysis was employed to quantitatively visualize the feature importance in the combined prediction model.
CECT serves as the preferred non-invasive imaging modality for preoperative evaluation of TETs, providing critical information regarding tumor morphology and vascularity that guides therapeutic decision-making (27). Given the high recurrence rates and poor prognosis associated with HRT, comprehensive preoperative assessment are essential (28). This underscores the clinical importance of CECT-based risk stratification prior to intervention. The arterial phase imaging is susceptible to beam-hardening artifacts from high-concentration contrast media in the superior vena cava, which may compromise image quality and distort texture analysis of adjacent lesions (29). Therefore, our study specifically employed venous phase images for radiomic feature extraction and subsequent model development.
Radiomics represented a robust and well-validated quantitative imaging analysis methodology that facilitated high-throughput extraction of minable data from radiological images (30, 31). By leveraging advanced machine learning algorithms, this approach enabled the identification of clinically relevant imaging biomarkers and the development of predictive models for precise estimation of clinical outcomes and therapeutic endpoints. Prior investigations had demonstrated the superior performance of radiomics in TETs risk stratification, with enhanced discriminative capability compared to conventional clinical model (18, 29). In the current study, our conventional radiomics model demonstrated favorable performance, achieving an AUC of 0.819 and an accuracy of 0.750 the independent test set, with balanced sensitivity and specificity indicating minimal class prediction bias. Ren et al. reported a radiomics model for TET histologic subtyping with an AUC of 0.860 and accuracy of 0.750, while Zhou et al. achieved an AUC of 0.716 and accuracy of 0.736 in predicting histological risk categorization (18, 32). The slight variations of observed performance may be attributable to methodological differences in feature selection and algorithmic implementation. Our approach utilized RFECV, an iterative selection process that optimizes feature combinations while preserving synergistic predictive value and model generalizability. Furthermore, the XGBoost algorithm’s tree-based architecture could accept feature collinearity while autonomously capturing critical nonlinear relationships through hierarchical feature interactions.
The complex cellular and molecular milieu surrounding tumor cells - the tumor microenvironment - critically influenced oncogenesis, disease progression, and therapeutic outcomes (33). Most prior studies on the risk stratification of TETs had primarily focused on whole-tumor analysis, while largely neglecting the potential significance of intratumoral subregion (23). Our study implemented habitat imaging to explore the incremental value of characterizing spatial tumor heterogeneity in TETs. We systematically partitioned TETs into three distinct subregions and extracted habitat-specific features. Habitat analysis identified distinct subregions: Habitat_1/2 represented solid tumor components (Habitat_2 showing lower enhancement), while Habitat_3 in corresponded to necrotic foci or tumor margin. Among 16 features developed combined model, 7 were conventional radiomic features, and other 9 were habitat-derived features. Quantitative analysis of tumor subregional architecture provided enhanced characterization of intratumoral spatial heterogeneity and enabled more biologically faithful representation of tumor phenotypic features. The combined model achieved superior discriminative performance (AUC=0.900, accuracy=0.864) in the test set compared to either feature set alone. Although DeLong test revealed no statistically significant AUC difference between conventional and combined models (p>0.05), the numerical improvement and enhanced classification metrics substantiated the incremental value of habitat radiomics for TETs risk stratification, suggesting that spatial heterogeneity provided complementary biological information.
The radiomic feature selection for the combined predictive model identified three primary discriminators with significant prognostic value: Ibp-3D-m1_glszm_ZonePercentage, Ibp-3D-m2_glcm_DifferenceVariance, and original_shape_Sphericity. These biomarkers collectively characterize distinct tumor biological behaviors, where the high-risk cohort exhibited (1): decreased Sphericity values, indicative of irregular morphological patterns and invasive growth tendencies (2); elevated DifferenceVariance, reflecting marked intratumoral textural heterogeneity; and (3) reduced ZonePercentage, suggesting disordered spatial distribution of tumor zones. Furthermore, three additional habitat-specific features provided complementary pathophysiological insights: original_shape_Maximum2DDiameterSlice_1 demonstrated strong correlation with proliferative activity in solid tumor components. wavelet-LIH_glcm_SumSquares_2 precisely quantified textural heterogeneity within viable tumor regions, where heightened values corresponded to areas of cellular atypia and structural disorganization. wavelet-HHH_firstorder_Range_3 effectively captured the complexity of necrotic and marginal zones through density dispersion metrics, with elevated values indicating pathological processes including hemorrhagic transformation, dystrophic calcification, or residual tumor infiltration.
Unlike manually crafted radiomics features, deep learning autonomously extracted task-specific features directly from images. With rapid advancements in deep learning methodologies, these data-driven features had emerged as a powerful complement to conventional radiomics features in medical imaging (34). This advantage is exemplified in the work of Zhou et al, who systematically compared conventional radiomics, deep learning, and combined models for the risk stratification of TETs (32). Their findings demonstrated that the combined model significantly outperformed conventional radiomics (AUC improvement from 0.716 to 0.786; accuracy increase from 0.736 to 0.774), likely due to the complementary nature of hand-crafted radiomic features and deep learning-derived representations (32). Building upon the established risk stratification framework of Zhou et al., our study demonstrates that integrating habitat-specific radiomic features yields significant improvements in predictive performance compared to conventional approaches. The combined model exhibited enhanced diagnostic capability, with accuracy increasing from 0.750 to 0.864 and AUC improving from 0.819 to 0.900. This performance surpasses previously reported results (accuracy: 0.864 vs 0.774; AUC: 0.900 vs 0.786), representing a meaningful advancement in predictive modeling. Notably, our findings reveal that habitat feature incorporation provides greater performance enhancement in conventional radiomics models, suggesting these biologically relevant features may capture distinct and complementary tumor characteristics.
Despite its methodological rigor, this study had several limitations that should be acknowledged. First, the retrospective nature of the analysis might still introduce selection bias, particularly given the specific inclusion/exclusion criteria applied. Second, although SVMSMOTE was implemented to address class imbalance between high- and low-risk TETs (37.3% vs. 62.7%), the fundamental disproportion in WHO subtype prevalence might still influence model performance in real-world clinical settings. Third, the utilization of multiple CT scanners, while reflecting clinical reality, introduced inherent technical variability despite rigorous standardization protocols. Fourth, while interobserver variability in manual segmentation was mitigated through ICC-based filtering, residual subjectivity inherent in manual annotations may persist. Finally, although not employed in this study, the harmonization technique (e.g., ComBat) could further improve feature robustness in future multi-center investigations by explicitly addressing scanner-induced heterogeneity.
In conclusion, the present study establishes that CECT-based habitat radiomics offered significant improvements in risk stratification for TETs compared to conventional radiomic approaches, particularly within clinically relevant decision-making ranges. The enhanced predictive performance of our combined model substantiated the importance of characterizing intratumoral heterogeneity through advanced habitat analysis, demonstrating substantial potential for clinical translation.
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 authors.
Ethics statement
Ethical approval was not required for the study involving humans in accordance with the local legislation and institutional requirements. Written informed consent to participate in this study was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and the institutional requirements.
Author contributions
YW: Investigation, Writing – review & editing, Writing – original draft, Conceptualization, Resources. ZS: Writing – review & editing, Formal analysis, Supervision. QS: Formal analysis, Writing – review & editing. TW: Investigation, Writing – review & editing, Validation, Methodology, Software, Visualization, Data curation, Writing – original draft. ZR: Writing – review & editing, Validation, Formal analysis, Resources, Software, Investigation, Methodology, Visualization. JZ: Resources, Formal analysis, Writing – review & editing, Methodology, Data curation, Conceptualization, Project administration.
Funding
The author(s) declare that no financial support was received for the research and/or publication of this article.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2025.1630485/full#supplementary-material
References
1. Özçıbık Işık G and Turna A. Surgical treatment of thymic epithelial tumor and myasthenia gravis. Front Surg. (2024) 11:1467789. doi: 10.3389/fsurg.2024.1467789
2. Duan L-W, Li J, Liu S-Q, Wang C-C, Liu H-K, Hu S-L, et al. Screening of optimal candidates for postoperative radiotherapy among patients with resectable thymic tumors: a multicenter retrospective study. Lung Cancer Amst Neth. (2025) 203:108542. doi: 10.1016/j.lungcan.2025.108542
3. Rioja P, Ruiz R, Galvez-Nino M, Lozano S, Valdiviezo N, Olivera M, et al. Epidemiology of thymic epithelial tumors: 22-years experience from a single-institution. Thorac Cancer. (2020) 12:420–5. doi: 10.1111/1759-7714.13760
4. Kuhn E, Pescia C, Mendogni P, Nosotti M, and Ferrero S. Thymic epithelial tumors: an evolving field. Life. (2023) 13:314. doi: 10.3390/life13020314
5. Yu F, Gu Z, Zhang X, Xu N, Hao X, Wang C, et al. A re-examination of neoadjuvant therapy for thymic tumors: A long and winding road. Cancers. (2024) 16:1680. doi: 10.3390/cancers16091680
6. Park MS, Chung KY, Kim KD, Yang WI, Chung JH, Kim YS, et al. Prognosis of thymic epithelial tumors according to the new World Health Organization histologic classification. Ann Thorac Surg. (2004) 78:992–7. doi: 10.1016/j.athoracsur.2004.03.097
7. Lee GD, Kim HR, Choi SH, Kim Y-H, Kim DK, and Park S-I. Prognostic stratification of thymic epithelial tumors based on both Masaoka-Koga stage and WHO classification systems. J Thorac Dis. (2016) 8:901–10. doi: 10.21037/jtd.2016.03.53
8. D’Angelillo RM, Trodella L, Ramella S, Cellini N, Balducci M, Mantini G, et al. Novel prognostic groups in thymic epithelial tumors: assessment of risk and therapeutic strategy selection. Int J Radiat Oncol. (2008) 71:420–7. doi: 10.1016/j.ijrobp.2007.10.006
9. Lee HS, Oh JS, Park YS, Jang SJ, Choi IS, and Ryu J-S. Differentiating the grades of thymic epithelial tumor Malignancy using textural features of intratumoral heterogeneity via 18F-FDG PET/CT. Ann Nucl Med. (2016) 30:309–19. doi: 10.1007/s12149-016-1062-2
10. Nakajo M, Takeda A, Katsuki A, Jinguji M, Ohmura K, Tani A, et al. The efficacy of 18F-FDG-PET-based radiomic and deep-learning features using a machine-learning approach to predict the pathological risk subtypes of thymic epithelial tumors. Br J Radiol. (2022) 95:20211050. doi: 10.1259/bjr.20211050
11. Chiappetta M, Mendogni P, Cattaneo M, Evangelista J, Farina P, Pizzuto DA, et al. Is PET/CT able to predict histology in thymic epithelial tumours? A narrative review. Diagnostics. (2022) 13:98. doi: 10.3390/diagnostics13010098
12. Dai H, Lan B, Li S, Huang Y, Jiang G, and Tian J. Prognostic CT features in patients with untreated thymic epithelial tumors. Sci Rep. (2023) 13:2910. doi: 10.1038/s41598-023-30041-z
13. Zhou Q, Huang X, Xie Y, Liu X, Li S, and Zhou J. Role of quantitative energy spectrum CT parameters in differentiating thymic epithelial tumours and thymic cysts. Clin Radiol. (2022) 77:136–41. doi: 10.1016/j.crad.2021.10.017
14. Hu Y-C, Wu L, Yan L-F, Wang W, Wang S-M, Chen B-Y, et al. Predicting subtypes of thymic epithelial tumors using CT: new perspective based on a comprehensive analysis of 216 patients. Sci Rep. (2014) 4:6814. doi: 10.1038/srep06984
15. Gao R, Zhou J, Zhang J, Zhu J, Wang T, and Yan C. Quantitative CT parameters combined with preoperative systemic inflammatory markers for differentiating risk subgroups of thymic epithelial tumors. BMC Cancer. (2023) 23:1158. doi: 10.1186/s12885-023-11332-0
16. Yu C, Li T, Yang X, Xin L, Zhao Z, Yang Z, et al. The maximal contrast-enhanced range of CT for differentiating the WHO pathological subtypes and risk subgroups of thymic epithelial tumors. Br J Radiol. (2023) 96:20221076. doi: 10.1259/bjr.20221076
17. Yu C, Li T, Zhang R, Yang X, Yang Z, Xin L, et al. Dual-energy CT perfusion imaging for differentiating WHO subtypes of thymic epithelial tumors. Sci Rep. (2020) 10:5511. doi: 10.1038/s41598-020-62466-1
18. Ren C, Li M, Zhang Y, and Zhang S. Development and validation of a CT-texture analysis nomogram for preoperatively differentiating thymic epithelial tumor histologic subtypes. Cancer Imaging. (2020) 20:86. doi: 10.1186/s40644-020-00364-5
19. Blüthgen C, Patella M, Euler A, Baessler B, Martini K, von Spiczak J, et al. Computed tomography radiomics for the prediction of thymic epithelial tumor histology, TNM stage and myasthenia gravis. PloS One. (2021) 16:e0261401. doi: 10.1371/journal.pone.0261401
20. Rajamohan N, Goyal A, Kandasamy D, Bhalla AS, Parshad R, Jain D, et al. CT texture analysis in evaluation of thymic tumors and thymic hyperplasia: correlation with the international thymic Malignancy interest group (ITMIG) stage and WHO grade. Br J Radiol. (2021) 94:20210583. doi: 10.1259/bjr.20210583
21. Ozawa Y, Hara M, and Shibamoto Y. The impact of radiomics in predicting oncologic behavior of thymic epithelial tumors. Mediastinum. (2019) 3:27. doi: 10.21037/med.2019.06.03
22. Liu J, Yin P, Wang S, Liu T, Sun C, and Hong N. CT-based radiomics signatures for predicting the risk categorization of thymic epithelial tumors. Front Oncol. (2021) 11:628534. doi: 10.3389/fonc.2021.628534
23. Lu X-F and Zhu T-Y. Diagnostic performance of radiomics model for preoperative risk categorization in thymic epithelial tumors: a systematic review and meta-analysis. BMC Med Imaging. (2023) 23:115. doi: 10.1186/s12880-023-01083-6
24. Jia P, Li Y, Li H, Li Y, Qin H, Xie A, et al. Habitat radiomics assists radiologists in accurately diagnosing lymph node metastasis of adenocarcinoma of the esophagogastric junction. Insights Imaging. (2025) 16:90. doi: 10.1186/s13244-025-01969-9
25. Li S, Dai Y, Chen J, Yan F, and Yang Y. MRI-based habitat imaging in cancer treatment: current technology, applications, and challenges. Cancer Imaging Off Publ Int Cancer Imaging Soc. (2024) 24:107. doi: 10.1186/s40644-024-00758-9
26. Mu W, Liang Y, Hall LO, Tan Y, Balagurunathan Y, Wenham R, et al. (18)F-FDG PET/CT habitat radiomics predicts outcome of patients with cervical cancer treated with chemoradiotherapy. Radiology. (2020) 2:e190218. doi: 10.1148/ryai.2020190218
27. Koyasu S. Imaging of thymic epithelial tumors-a clinical practice review. Mediastinum Hong Kong China. (2024) 8:41. doi: 10.21037/med-23-66
28. Khandelwal A, Sholl LM, Araki T, Ramaiya NH, Hatabu H, and Nishino M. Patterns of metastasis and recurrence in thymic epithelial tumours: longitudinal imaging review in correlation with histological subtypes. Clin Radiol. (2016) 71:1010–17. doi: 10.1016/j.crad.2016.05.007
29. Chen X, Feng B, Xu K, Chen Y, Duan X, Jin Z, et al. Development and validation of a deep learning radiomics nomogram for preoperatively differentiating thymic epithelial tumor histologic subtypes. Eur Radiol. (2023) 33:6804–16. doi: 10.1007/s00330-023-09690-1
30. Inglese M, Conti A, and Toschi N. Radiomics across modalities: a comprehensive review of neurodegenerative diseases. Clin Radiol. (2025) 85:106921. doi: 10.1016/j.crad.2025.106921
31. Xu L, Bressem K, Adams L, Poddubnyy D, and Proft F. AI for imaging evaluation in rheumatology: applications of radiomics and computer vision-current status, future prospects and potential challenges. Rheumatol Adv Pract. (2025) 9:rkae147. doi: 10.1093/rap/rkae147
32. Zhou H, Bai HX, Jiao Z, Cui B, Wu J, Zheng H, et al. Deep learning-based radiomic nomogram to predict risk categorization of thymic epithelial tumors: A multicenter study. Eur J Radiol. (2023) 168:111136. doi: 10.1016/j.ejrad.2023.111136
33. Li Z, Li J, Bai X, Huang X, and Wang Q. Tumor microenvironment as a complex milieu driving cancer progression: a mini review. Clin Transl Oncol. (2024) 27:1943–52. doi: 10.1007/s12094-024-03697-w
Keywords: thymic epithelial tumors, habitat radiomics, risk stratification, computed tomography, machine learning
Citation: Wang Y, Shi Z, Sun Q, Wang T, Ran Z and Zhang J (2025) Incremental diagnostic value of tumor habitat radiomics for risk stratification in thymic epithelial tumors. Front. Oncol. 15:1630485. doi: 10.3389/fonc.2025.1630485
Received: 17 May 2025; Accepted: 25 July 2025;
Published: 03 September 2025.
Edited by:
Morgan Michalet, Institut du Cancer de Montpellier (ICM), FranceReviewed by:
Catharina Lisson, Ulm University Medical Center, GermanyYuChuan Hu, Air Force Medical University, China
Copyright © 2025 Wang, Shi, Sun, Wang, Ran and Zhang. 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: Tianzuo Wang, YWdudHd6QDE2My5jb20=; Zhizhen Ran, NDU1NzYwNTI4QHFxLmNvbQ==; Jinling Zhang, SmlubGluZ3poYW5nemlAMTYzLmNvbQ==