RESUMEN
BACKGROUND. Tumor growth processes result in spatial heterogeneity, with the development of tumor subregions (i.e., habitats) having unique biologic characteristics. OBJECTIVE. The purpose of our study was to develop and validate a habitat model combining tumor and peritumoral radiomic features on chest CT for predicting invasiveness of lung adenocarcinoma. METHODS. This retrospective study included 1156 patients (mean age, 57.5 years; 464 men, 692 women), from three centers and a public dataset, who underwent chest CT before lung adenocarcinoma resection (variable date ranges across datasets). Patients from one center formed training (n = 500) and validation (n = 215) sets; patients from the other sources formed three external test sets (n = 249, 113, 79). For each patient, a single nodule was manually segmented on chest CT. The nodule segmentation was combined with an automatically generated 4-mm peritumoral region into a whole-volume volume of interest (VOI). A gaussian mixture model (GMM) identified voxel clusters with similar first-order energy across patients. GMM results were used to divide each patient's whole-volume VOI into multiple habitats, which were defined consistently across patients. Radiomic features were extracted from each habitat. After feature selection, a habitat model was developed for predicting invasiveness, with the use of pathologic assessment as a reference. An integrated model was constructed, combining features extracted from habitats and whole-volume VOIs. Model performance was evaluated, including in subgroups based on nodule density (pure ground-glass, part-solid, and solid). The code for habitat imaging and model construction is publicly available (https://github.com/Shangyoulan/Habitat/). RESULTS. Invasive cancer was diagnosed in 626 of 1156 patients. GMM identified four as the optimal number of voxel clusters and thus of per-patient tumor habitats. The habitat model had an AUC of 0.932 in the validation set and 0.881, 0.880, and 0.764 in the three external test sets. The integrated model had an AUC of 0.947 in the validation set and 0.936, 0.908, and 0.800 in the three external test sets. In the three external test sets combined, across nodule densities, AUCs for the habitat model were 0.836-0.869 and for the integrated model were 0.846-0.917. CONCLUSION. Habitat imaging combining tumoral and peritumoral radiomic features could help predict lung adenocarcinoma invasiveness. Prediction is improved when combining information on tumor subregions and the tumor overall. CLINICAL IMPACT. The findings may aid personalized preoperative assessments to guide clinical decision-making in lung adenocarcinoma.