Identification of prognostic genes and tumor-infiltrating immune cells in the tumor microenvironment of esophageal squamous cell carcinoma and esophageal adenocarcinoma

Qilin Huai; Wei Guo; Liankui Han; Demiao Kong; Liang Zhao; Peng Song; Yue Peng; Shugeng Gao

doi:10.21037/tcr-20-3078

Identification of prognostic genes and tumor-infiltrating immune cells in the tumor microenvironment of esophageal squamous cell carcinoma and esophageal adenocarcinoma

Transl Cancer Res. 2021 Apr;10(4):1787-1803. doi: 10.21037/tcr-20-3078.

Authors

Qilin Huai^#^{1

2}, Wei Guo^#³, Liankui Han^#², Demiao Kong², Liang Zhao³, Peng Song³, Yue Peng³, Shugeng Gao³

Affiliations

¹ Department of Graduate School, Zunyi Medical University, Zunyi, China.
² Department of Thoracic Surgery, Guizhou Provincial People's Hospital, Guiyang, China.
³ Department of Thoracic Surgery, National Cancer Center/National Clinical Research Center for Cancer/Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China.

^# Contributed equally.

Abstract

Background: Esophageal cancer (EC) is a highly aggressive malignancy that is classified as esophageal squamous cell carcinoma (ESCC) and esophageal adenocarcinoma (EAC). Infiltrating stromal/immune cells, a major component of the tumor immune microenvironment (TIME), have prognostic significance in various cancers.

Methods: In this study we investigated genes and immune factors in the tumor microenvironment (TME) of ESCC and EAC that can serve as prognostic biomarkers. Stromal and immune scores were calculated using the Estimation of Stromal and Immune Cells in Malignant Tumor Tissues Using Expression Data (ESTIMATE) algorithm based on gene expression profiles of patient-derived tumor tissues in The Cancer Genome Atlas database. The correlation between ESTIMATE scores and survival rates in EC were analyzed. A comparison of high and low stromal/immune score groups revealed multiple differentially expressed genes (DEGs) as candidate prognostic genes; their role in immune-related biological processes was evaluated by functional and protein-protein interaction (PPI) network analyses, and the genes were validated using Gene Expression Omnibus datasets. Additionally, 22 tumor-infiltrating immune cell (TIIC) subsets were analyzed using the CIBERSORT algorithm.

Results: Median stromal score was higher whereas immune score was lower in ESCC than in EAC (both P<0.01). Stromal score was lower in female as compared to male ESCC patients (P<0.05), and was significantly correlated with T stage (P<0.05). In EAC, median immune score was higher in female as compared to male patients (P<0.05) and was correlated with tumor-node-metastasis stage (P<0.05). The identified DEGs were mainly involved in lymphocyte (especially T-lymphocyte) activation and carbohydrate binding. Moreover, the levels of infiltrating resting-stage dendritic cells, CD8+ T cells, naïve B cells, activated mast cells, and resting memory CD4+ T cells were significantly correlated with EC prognosis (P<0.05).

Conclusions: The immune microenvironment of ESCC and EAC are quite different. We have found genes with prognostic value in multiple tumor databases.

Keywords: Esophageal cancer (EC); The Cancer Genome Atlas (TCGA); immune infiltration; prognosis; tumor immune microenvironment (TIME).