Automatically Quantifying Spatial Heterogeneity of Immune and Tumor Hypoxia Environment and Predicting Disease Free Survival for Patients with Rectal Cancer

Abstract Immunohistochemistry (IHC) remains the gold standard for evaluating protein expression in tumor microenvironment analysis. This approach hinders robust correlation analyses between spatial heterogeneity in the tumor microenvironment and clinical outcomes like disease-free survival (DFS). To address these challenges, we developed an automated pipeline for quantitative IHC feature extraction. Our method integrates deep learning-based tumor segmentation with computational detection of invasive margins at varying distances. Deconvolution algorithms quantify diaminobenzidine (DAB) staining intensity across the tumor body and the invasive margin. Spatial heterogenetic DAB density patterns were subsequently analyzed for DFS correlation. Using 104 patient samples (57 training/47 validation) stained for CD3, CD8, CD31, and HIF-1α, we identified two prognostic feature categories (CD3/CD8 co-localization positive areas within the 0.25mm peripheral zone extending outward from the tumor-invasive front and HIF1-α-positive areas within a 0.75mm peripheral zone extending outward from the tumor-invasive front). Immune-related features demonstrated C-indices of 0.726 (training) and 0.626 (validation), while hypoxia-associated markers showed C-indices of 0.714 and 0.656, respectively. Integration of these features with pTNM staging enhanced DFS stratification compared to pTNM staging alone, improving C-indices from 0.702 to 0.819 (training) and 0.668 to 0.853 (validation). This automated pipeline addresses critical limitations in traditional IHC analysis by enabling: 1) Objective quantification of spatial DAB heterogeneity. 2) Identification of biologically interpretable prognostic features. 3) Enhanced predictive performance over conventional staging systems. Our findings suggest this methodology could standardize IHC-based prognostic assessments and inform personalized treatment strategies. Further validation in multicenter cohorts is warranted to confirm clinical applicability.