Ultrasound-Based Habitat Radiomics for Differentiating Phyllodes Tumor from Breast Fibroadenoma.
Authors
Affiliations (7)
Affiliations (7)
- Department of Ultrasound, the Affiliated Hospital of Qingdao University, Qingdao 266003, China (C.X., F.C., S.L., C.Z., P.L.); Department of Interventional Ultrasound, Fifth Medical Center of Chinese PLA General Hospital, Beijing, China (C.X., W.D., Z.L., J.Y., P.L.).
- Department of Interventional Ultrasound, Fifth Medical Center of Chinese PLA General Hospital, Beijing, China (C.X., W.D., Z.L., J.Y., P.L.).
- Department of Ultrasound, the Affiliated Hospital of Qingdao University, Qingdao 266003, China (C.X., F.C., S.L., C.Z., P.L.).
- Department of Ultrasound, Ruijin Hospital, Shanghai Jiaotong University School of Medicine, Shanghai, China (X.J., J.Z.); Faculty of Medical Imaging Technology, College of Health Science and Technology, Shanghai Jiao Tong University School of Medicine, Shanghai, China (X.J., J.Z.).
- Department of Ultrasound, The First Affiliated Hospital of Anhui Medical University, Anhui, China (J.L.); Department of Ultrasound, Anhui Public Health Clinical Center, Anhui, China (J.L.).
- Senior Department of General Surgery, Chinese PLA General Hospital, Beijing 100853, China (J.S.).
- Department of Ultrasound, the Affiliated Hospital of Qingdao University, Qingdao 266003, China (C.X., F.C., S.L., C.Z., P.L.); Department of Interventional Ultrasound, Fifth Medical Center of Chinese PLA General Hospital, Beijing, China (C.X., W.D., Z.L., J.Y., P.L.). Electronic address: [email protected].
Abstract
Phyllodes tumors (PTs) and fibroadenomas (FAs) have overlapping ultrasound (US) appearances. This study aimed to develop and evaluate an interpretable machine-learning model integrating US radiomics, habitat features, and clinical variables. This retrospective multicenter study included 1950 pathologically confirmed lesions (1029 FAs and 921 PTs) divided into training (n=918), validation (n=474), and test (n=558) cohorts. Nineteen voxelwise features were used for K-means habitat segmentation. The number of clusters was selected from K=2-10 in the training cohort and evaluated across 50 random seeds. Whole-lesion, habitat, and combined habitat-radiomics (HR) signatures were developed using a train-only feature-selection pipeline, and ExtraTrees was selected as the classifier. The final combined model was built by probability-level late fusion of the ExtraTrees HR model and a logistic clinical model based on age, lesion size, and orientation. K=3 yielded the highest Calinski-Harabasz index (200,927.01). Clustering was stable across 50 random seeds, with mean adjusted Rand indices of 0.9942, 0.9948, and 0.9941 in the training, validation, and test cohorts, respectively. Feature selection retained 39 whole-lesion features and 34 habitat features (habitats 1-3: 1, 14, and 19). The combined model achieved AUCs of 0.921 (95% CI, 0.904-0.939), 0.797 (95% CI, 0.754-0.835), and 0.844 (95% CI, 0.813-0.877) in the training out-of-fold, validation, and test cohorts. In this multicenter cohort, a model integrating US habitat radiomics with age, lesion size, and orientation differentiated PTs from FAs with consistent discrimination across data partitions. Shapley additive explanations analysis provided global and case-level explanations of the combined prediction.