Radiomics prognostic analysis of PET/CT images in a multicenter head and neck cancer cohort: investigating ComBat strategies, sub-volume characterization, and automatic segmentation

Purpose This study aimed to investigate the impact of several ComBat harmonization strategies, intra-tumoral sub-volume characterization, and automatic segmentations for progression-free survival (PFS) prediction through radiomics modeling for patients with head and neck cancer (HNC) in PET/CT images. Methods The HECKTOR MICCAI 2021 challenge set containing PET/CT images and clinical data of 325 oropharynx HNC patients was exploited. A total of 346 IBSI-compliant radiomic features were extracted for each patient’s primary tumor volume defined by the reference manual contours. Modeling relied on least absolute shrinkage Cox regression (Lasso-Cox) for feature selection (FS) and Cox proportional-hazards (CoxPH) models were built to predict PFS. Within this methodological framework, 8 different strategies for ComBat harmonization were compared, including before or after FS, in feature groups separately or all features directly, and with center or clustering-determined labels. Features extracted from tumor sub-volume clustering were also investigated for their prognostic additional value. Finally, 3 automatic segmentations (2 threshold-based and a 3D U-Net) were also compared. All results were evaluated with the concordance index ( C -index). Results Radiomics features without harmonization, combined with clinical factors, led to models with C -index values of 0.69 in the testing set. The best version of ComBat harmonization, i.e., after FS, for feature groups separately and relying on clustering-determined labels, achieved a C -index of 0.71. The use of features extracted from tumor sub-volumes further improved the C -index to 0.72. Models that relied on the automatic segmentations yielded close but slightly lower prognostic performance (0.67–0.70) compared to reference contours. Conclusion A standard radiomics pipeline allowed for prediction of PFS in a multicenter HNC cohort. Applying a specific strategy of ComBat harmonization improved the performance. The extraction of intra-tumoral sub-volume features and automatic segmentation could contribute to the improvement and automation of prognosis modeling, respectively.