Computational image features of immune architecture is associated with clinical benefit and survival in gynecological cancers across treatment modalities

Background We present a computational approach (ArcTIL) for quantitative characterization of the architecture of tumor-infiltrating lymphocytes (TILs) and their interplay with cancer cells from digitized H&E-stained histology whole slide images and evaluate its prognostic role in three different gynecological cancer (GC) types and across three different treatment types (platinum, radiation and immunotherapy). Methods In this retrospective study, we included 926 patients with GC diagnosed with ovarian cancer (OC), cervical cancer, and endometrial cancer with available digitized diagnostic histology slides and survival outcome information. ArcTIL features quantifying architecture and spatial interplay between immune cells and the rest of nucleated cells (mostly comprised cancer cells) were extracted from the cell cluster graphs of nuclei within the tumor epithelial nests, surrounding stroma and invasive tumor front compartments on H&E-stained slides. A Cox proportional hazards model, incorporating ArcTIL features was fit on the OC training cohort (N=51), yielding an ArcTIL signature. A unique threshold learned from the training set stratified the patients into a low and high-risk group. Results The seven feature ArcTIL classifier was found to significantly correlate with overall survival in chemotherapy and radiotherapy-treated validation cohorts and progression-free survival in an immunotherapy-treated validation cohort. ArcTIL features relating to increased density of TILs in the epithelium and invasive tumor front were found to be associated with better survival outcomes when compared with those patients with an increased TIL density in the stroma. A statistically significant association was found between the ArcTIL signature and signaling pathways for blood vessel morphogenesis, vasculature development, regulation of cell differentiation, cell-substrate adhesion, biological adhesion, regulation of vasculature development, and angiogenesis. Conclusions This study reveals that computationally-derived features from the spatial architecture of TILs and tumor cells are prognostic in GCs treated with chemotherapy, radiotherapy, and checkpoint blockade and are closely associated with central biological processes that impact tumor progression. These findings could aid in identifying therapy-refractory patients and further enable personalized treatment decision-making.