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

ORIGINAL RESEARCH article

Front. Oncol., 03 September 2025

Sec. Thoracic Oncology

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

This article is part of the Research TopicRadiomics and Artificial Intelligence in Oncology ImagingView all 17 articles

Incremental diagnostic value of tumor habitat radiomics for risk stratification in thymic epithelial tumors

Yiqiao WangYiqiao Wang1Zhe ShiZhe Shi2Qinliang SunQinliang Sun3Tianzuo Wang*Tianzuo Wang2*Zhizhen Ran*Zhizhen Ran1*Jinling Zhang*Jinling Zhang1*
  • 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 (13). 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 (68). 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) (911). 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 (1417). 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 (1921). 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 (2426). 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.

Figure 1
Flowchart depicting a process divided into four stages: Segmentation, Feature Extraction, Feature Selection, and Modeling & Evaluation. Under Segmentation are images of a tumor and its regions. Feature Extraction includes first-order, shape, and texture features shown as graphs. Feature Selection shows ICC, feature importance, and RFECV charts. Modeling & Evaluation includes XGBoost, ROC, and calibration plots.

Figure 1. Flowchart of the research.

Figure 2
Flowchart depicting a study with 229 patients with pathologically confirmed TETs, where 9 were excluded due to recurrent lesions, poor-quality images, or undetermined WHO classification. The remaining 220 patients were split into a training set of 176 (66 high-risk, 110 low-risk) and a test set of 44 (16 high-risk, 28 low-risk).

Figure 2. Workflow of patient inclusion.

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.

Table 1
www.frontiersin.org

Table 1. Baseline characteristics.

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.

Figure 3
Elbow plot showing within-cluster sum of squared errors (SSE) against the number of clusters (k) from one to ten. A significant bend, indicating the optimal number of clusters, is highlighted at k equals three.

Figure 3. Line plot of finding optimal k.

Figure 4
Bar chart showing mean SHAP values for different features in two categories: whole (gray) and subregion (cyan). The top features include lbp-3D-m1_glszm_ZonePercentage and original_shape_Sphericity, with lbp-3D-m1_glszm_ZonePercentage having the highest value. The x-axis represents the mean SHAP value.

Figure 4. SHAP bar plot (combined model).

Figure 5
Heatmap visualization showing SHAP values for features across instances in two panels, A and B. Each panel displays feature impact on model output, with red indicating positive influence and blue indicating negative influence. Key features include lbp-3D-m1_glszm_ZonePercentage and original_shape_Sphericity. Graphs above heatmaps plot f(x) values. Color bars on the right indicate SHAP value ranges. Panel A covers instances from 0 to 176, while Panel B covers instances from 0 to 44.

Figure 5. SHAP heatmap (combined model): training set (A), test set (B).

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

Table 2
www.frontiersin.org

Table 2. Performance comparison of models.

Figure 6
Three charts labeled a, b, and c display the performance of different prediction models. Chart a shows ROC curves comparing true positive rates to false positive rates. Chart b depicts calibration plots of fraction of positives versus mean predicted probability. Chart c presents decision curves for standardized net benefit versus threshold probability. The models, Radiomics, Habitat, Combined, and All, are represented by brown, orange, cyan, and gray lines, respectively. A legend on the right identifies the model colors.

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

34. Truhn D, Schrading S, Haarburger C, Schneider H, Merhof D, and Kuhl C. Radiomic versus Convolutional Neural Networks Analysis for Classification of Contrast-enhancing Lesions at Multiparametric Breast MRI. Radiology. (2019) 290:290–7. doi: 10.1148/radiol.2018181352

PubMed Abstract | Crossref Full Text | Google Scholar

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), France

Reviewed by:

Catharina Lisson, Ulm University Medical Center, Germany
YuChuan 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==

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.