Abstract
PURPOSE: Pediatric adrenocortical tumors (pACTs) are rare and clinically heterogeneous. Existing risk stratification systems rely on fixed thresholds and linear assumptions, which may limit their prognostic accuracy-particularly for nonmetastatic, locally advanced cases. We aimed to develop an interpretable machine learning (ML) model for individualized survival prediction using only routine clinical features. METHODS: We retrospectively analyzed 97 patients with pACT from the German Pediatric Oncology Hematology-Malignant Endocrine Tumors Registry (1997-2024). An Extreme Gradient Boosting Cox proportional hazards model was trained using 4 features-tumor volume, distant metastases, pathologic T stage, and resection status-identified via systematic feature evaluation across 11 737 model combinations. Performance was assessed using a stratified 80/20 train-test split, 500 bootstrap iterations, and Harrell's concordance index (C-index). SHapley Additive exPlanations (SHAP) were used for interpretability. RESULTS: The model achieved strong prognostic performance (test-set C-index: 0.925; bootstrap mean: 0.891, 95% confidence interval: 0.817-0.946). SHAP analysis confirmed the dominant influence of metastatic status, followed by tumor volume, T stage, and resection status. The model uncovered nonlinear and additive effects, including a SHAP- and bootstrap-guided tumor volume cut-off (190 mL, 95% confidence interval 127-910 mL) that only slightly differed from conventional thresholds. Stratification remained robust in subgroups, including nonmetastatic patients with advanced local disease. CONCLUSION: This interpretable ML model enables individualized survival prediction in pACT using only routine clinical data. It offers a clinically accessible and clinically meaningful complement to existing scoring systems, particularly in patients with ambiguous risk profiles who may benefit from more personalized management.