Distinct tumor microenvironment landscapes of rectal cancer for prognosis and prediction of immunotherapy response.

Abstract

Tumor microenvironment (TME) affects the progression of rectal cancer (RC), and the clinical relevance of its immune elements was widely reported. Here we aim to delineate the complete TME landscape, including non-immune features, to improve our understanding of RC heterogeneity and provide a better strategy for precision medicine. Single-cell analysis of GSE161277 using Seurat and Cellcall was performed to identify cell-cell interactions. The ssGSEA was employed to quantify the TME elements in TCGA patients, which were further clustered into subtypes by hclust. WGCNA and LASSO were combined to construct a degenerated signature for prognosis, and its performance was validated in two GEO datasets. We proposed a subtyping strategy based on the abundance of both immune and non-immune components, which divided all RC patients into 4 subtypes (Immune-, Canonical-, Dormant- and Stem-like). Different subtypes exhibited distinct mutation landscapes, biological features, immune characteristics, immunotherapy responses and prognoses. Next, WGCNA and LASSO regression were combined to construct a 10-gene signature based on differentially expressed genes among different subtypes. Subgroups divided by this signature also exhibited different clinical parameters and responses to immune checkpoint blockades. Diverse machine learning algorithms were applied to achieve higher accuracy for survival prediction and a nomogram was further established in combination with M stage and age to provide an accurate and visual prediction of prognosis. We identified four TME-based RC subtypes with distinct biological and clinical features. Based on those subtypes, we also proposed a degenerated 10-gene signature to predict the prognosis and immunotherapy response.