Abstract:Background Esophageal squamous cell carcinoma (ESCC) is the most common type of esophageal cancer and the seventh most prevalent cause of cancer-related death worldwide. Tumor microenvironment (TME) has been confirmed to play an crucial role in ESCC progression, prognosis, and the response to immunotherapy. There is a need for predictive biomarkers of TME-related processes to better prognosticate ESCC outcomes. Aim To identify a novel gene signature linked with the TME to predict the prognosis of ESCC. Methods We calculated the immune/stromal scores of 95 ESCC samples from The Cancer Genome Atlas (TCGA) using the ESTIMATE algorithm, and identified differentially expressed genes (DEGs) between high and low immune/stromal score patients. The key prognostic genes were further analyzed by the intersection of protein–protein interaction (PPI) networks and univariate Cox regression analysis. Finally, a risk score model was constructed using multivariate Cox regression analysis. We evaluated the associations between the risk score model and immune infiltration via the CIBERSORT algorithm. Moreover, we validated the signature using the Gene Expression Omnibus (GEO) database. Within the ten gene signature, five rarely reported genes were further validated with quantitative real time polymerase chain reaction (qRT-PCR) using an ESCC tissue cDNA microarray. Results A total of 133 up-regulated genes were identified as DEGs. Ten prognostic genes were selected based on intersection analysis of univariate COX regression analysis and PPI, and consisted of C1QA, C1QB, C1QC, CD86, C3AR1, CSF1R, ITGB2, LCP2, SPI1, and TYROBP (HR>1, p<0.05). The expression of 9 of these genes in the tumor samples were significantly higher compared to matched adjacent normal tissue based on the GEO database (p<0.05). Next, we assessed the ability of the ten-gene signature to predict the overall survival of ESCC patients, and found that the high-risk group had significantly poorer outcomes compared to the low-risk group using univariate and multivariate analyses in the TCGA and GEO cohorts (HR=2.104, 95% confidence interval:1.343-3.295, p=0.001; HR=1.6915, 95% confidence interval:1.053-2.717, p=0.0297). Additionally, receiver operating characteristic (ROC) curve analysis demonstrated a relatively sensitive and specific profile for the signature (1-, 2-, 3-year AUC=0.672, 0.854, 0.81). To identify the basis for these differences in the TME, we performed correlation analyses and found a significant positive correlation with M1 and M2 macrophages and CD8+ T cells, as well as a strong correlation to M2 macrophage surface markers. A nomogram based on the risk score and select clinicopathologic characteristics was constructed to predict overall survival of ESCC patients. For validation, qRT-PCR of an ESCC patient cDNA microarray was performed, and demonstrated that C1QA, C3AR1, LCP2, SPI1, and TYROBP were up-regulated in tumor samples and predict poor prognosis. Conclusion This study established and validated a novel 10-gene signature linked with M2 macrophages and poor prognosis in ESCC patients. Importantly, we identified C1QA, C3AR1, LCP2, SPI1, and TYROBP as novel M2 macrophage-correlated survival biomarkers. These findings may identify potential targets for therapy in ESCC patients.

Integrated genomic profiling and modelling for risk stratification in patients with advanced oesophagogastric adenocarcinoma

Genomic analysis of oesophageal carcinoma (EC) to identify recurrent mutations in histone methyltransferases as a distinctive subset.

Comprehensive Profiling of MDM2/TP53 Genomic Aberration in Chinese Patients with Diverse Malignancies.

An Integrated Model of Clinical Information and Gene Expression for Prediction of Survival in Ovarian Cancer Patients

The genomic landscape of oesophagogastric junctional adenocarcinoma

Genomic and transcriptomic alterations associated with drug vulnerabilities and prognosis in adenocarcinoma at the gastroesophageal junction

Gene expression profiling of mucinous ovarian tumors and comparison with upper and lower gastrointestinal tumors identifies markers associated with adverse outcomes.

Whole-genome sequencing of 508 patients identifies key molecular features associated with poor prognosis in esophageal squamous cell carcinoma

Abstract 3890: Comprehensive molecular and immunological characterization of early-onset esophagogastric cancer

Multiplex profiling of peritoneal metastases from gastric adenocarcinoma identified novel targets and molecular subtypes that predict treatment response

Revisiting MET: Clinical Characteristics and Treatment Outcomes of Patients with Locally Advanced or Metastatic, MET-Amplified Esophagogastric Cancers

Developing a prognosis and chemotherapy evaluating model for colon adenocarcinoma based on mitotic catastrophe-related genes

Molecular Phenotypes and Potential Therapeutic Targets in Esophagogastric Adenocarcinoma Unearthed by Whole Genome and Transcriptome Analyses

Multicentre validation of CT grey-level co-occurrence matrix features for overall survival in primary oesophageal adenocarcinoma

Characteristics and response to Next‐generation sequencing‐guided therapy in locally advanced or metastatic esophageal cancer

Enhancer demethylation-regulated gene score identified molecular subtypes, inspiring immunotherapy or CDK4/6 inhibitor therapy in oesophageal squamous cell carcinoma

Development and Validation of a Prognostic Gene Signature Correlated With M2 Macrophage Infiltration in Esophageal Squamous Cell Carcinoma

The Tumour Immune Microenvironment Drives Survival Outcomes and Therapeutic Response in an Integrated Molecular Analysis of Gastric Adenocarcinoma

Cell states and neighborhoods in distinct clinical stages of primary and metastatic esophageal adenocarcinoma

Prognostic value of genomic mutation signature associated with immune microenvironment in southern Chinese patients with esophageal squamous cell carcinoma