Identification and characterization of prognosis-related genes in the tumor microenvironment of esophageal squamous cell carcinoma

Int Immunopharmacol. 2021 Jul:96:107616. doi: 10.1016/j.intimp.2021.107616. Epub 2021 May 28.

Abstract

Background: Esophageal squamous cell carcinoma (ESCC) is the main pathological subtype of esophageal cancer with high incidence and mortality. Immune and stromal cells in the tumor microenvironment (TME) profoundly affect the development of ESCC.

Methods: In this study, we used the ESTIMATE algorithm to calculate the immune and stromal scores of ESCC samples in The Cancer Genome Atlas (TCGA) database. Next, we used the R package limma to identify differentially expressed genes (DEGs) from high- versus low-immune/stromal score groups and these DEGs were further utilized to analyze the functional annotations, protein-protein interaction (PPI) networks and overall survival of patients with ESCC. Finally, we identified the biological roles of core gene C3AR1 in the TME of ESCC using the TCGA database and in vitro experiments.

Results: We obtained the immune and stromal scores of ESCC samples and further evaluated the impact of these scores on the prognosis and clinical parameters of patients with ESCC. Next, we identified 410 DEGs from high- versus low-immune/stromal score groups and to gain better understanding of the biological functions and characteristics of DEGs. Among these DEGs, 69 were correlated with the overall survival of patients with ESCC and C3AR1 was identified as a core gene for the regulation of most genes in the network. We found that C3AR1 was positively correlated with M2 macrophages and immune inhibitory molecules (T-cell immunoglobulin and mucin domain 3 (TIM-3), programmed cell death-1 (PD-1)), but not with M1 macrophages. We also observed a higher expression of CD163 and CD206, which were the markers for M2 macrophages in the TLQP-21 TFA (the agonist of C3AR1)groups than in the control groups.

Conclusion: Based on the ESTIMATE algorithm, we obtained and characterized prognosis-related genes in the TME of ESCC samples from the TCGA database. We have further revealed that C3AR1 may cause an immunosuppressive microenvironment by affecting the polarization of macrophages to M2 phenotype and lead to the progression of ESCC, which indicates that C3AR1 may be a potential target for immunotherapy.

Keywords: Esophageal squamous cell carcinoma; Prognosis-related genes; TCGA; Tumor microenvironment.

MeSH terms

  • Antigens, CD / metabolism
  • Antigens, Differentiation, Myelomonocytic / metabolism
  • Biomarkers, Tumor / immunology
  • Computational Biology
  • Databases, Factual
  • Esophageal Neoplasms / genetics*
  • Esophageal Neoplasms / immunology*
  • Esophageal Neoplasms / mortality
  • Esophageal Squamous Cell Carcinoma / genetics*
  • Esophageal Squamous Cell Carcinoma / immunology*
  • Esophageal Squamous Cell Carcinoma / mortality
  • Female
  • Gene Expression Regulation, Neoplastic / immunology
  • Hepatitis A Virus Cellular Receptor 2 / metabolism
  • Humans
  • Kaplan-Meier Estimate
  • Macrophages / immunology
  • Macrophages / metabolism
  • Male
  • Membrane Glycoproteins / metabolism
  • Middle Aged
  • Prognosis
  • Programmed Cell Death 1 Receptor / metabolism
  • Protein Interaction Maps / genetics
  • Protein Interaction Maps / immunology
  • Receptors, Cell Surface / metabolism
  • Receptors, Complement / agonists
  • Receptors, Complement / immunology*
  • Receptors, Complement / metabolism
  • Receptors, Immunologic / metabolism
  • Tumor Microenvironment / genetics*
  • Tumor Microenvironment / immunology*

Substances

  • Antigens, CD
  • Antigens, Differentiation, Myelomonocytic
  • Biomarkers, Tumor
  • CD163 antigen
  • HAVCR2 protein, human
  • Hepatitis A Virus Cellular Receptor 2
  • MRC1 protein, human
  • Membrane Glycoproteins
  • PDCD1 protein, human
  • Programmed Cell Death 1 Receptor
  • Receptors, Cell Surface
  • Receptors, Complement
  • Receptors, Immunologic
  • complement C3a receptor