Dataset Information

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

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.

SUBMITTER: Yao J

PROVIDER: S-EPMC8677679 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Objective: The transcription factor forkhead box protein O1 (FOXO1) is critical for regulating cytokine and chemokine secretion. However, its function in the tumor microenvironment (TME) remains largely unexplored. In this study, we characterized the prognostic value of FOXO1 and the interaction between tumor-derived FOXO1 and M2 macrophages in esophageal squamous cell carcinoma (ESCC). Methods: FOXO1 expression and macrophage infiltration in clinical samples and mouse models were quantified using quantitative real-time polymerase chain reaction (qRT-PCR) and immunohistochemistry staining. Western blotting, qRT-PCR, and enzyme-linked immunosorbent assay were used to evaluate chemokine ligand 20 (CCL20) and colony stimulating factor 1 (CSF-1) expression in FOXO1(+) and FOXO1(-) tumor cells. Macrophage phenotypes were determined using qRT-PCR, flow cytometry, and RNA sequencing. Transcriptional activity was measured using chromatin immunoprecipitation (ChIP)-qPCR. Tumor viability was investigated using XTT proliferation and foci formation assays. Results: FOXO1 upregulation in tumor tissues was found to drive the polarization of M0 macrophages and infiltration of M2 macrophages into the TME, resulting in worse prognosis in ESCC patients. CSF-1, a vital factor inducing M0-to-M2 polarization, was upregulated via a FOXO1-mediated mechanism. RNA sequencing results corroborated that the FOXO1-induced macrophages exhibited similar molecular signatures to the IL4-stimulated M2 macrophages. The transwell assays showed that FOXO1 promoted the migration of M2 macrophages via CCL20 secretion, which could be inhibited using an anti-CCL20 antibody. FOXO1(+) tumor-induced M2 macrophages promoted tumor proliferation via the FAK-PI3K-AKT pathway and the PI3K inhibitor could effectively impede the oncogenical process. Conclusions: FOXO1 facilitated M0-to-M2 polarization and the recruitment of M2 macrophages in the TME via the transcriptional modulation of CCL20 and CSF-1. Our data deciphered the FOXO1-dependent mechanism in M2 macrophage infiltration in the TME of ESCC, which has implications for the development of novel prognostic and therapeutic targets to optimize the current treatment against ESCC.

Project description:Esophageal carcinoma (EC) is a common malignant cancer worldwide. Esophageal squamous cell carcinoma (ESCC), the main type of EC, is difficult to treat because of the widespread morbidity, high fatality rates, and low quality of life caused by postoperative complications and no specific molecular target. In this study, we screened genes to establish a prognostic model for ESCC. The transcriptome expression profiles of 81 ESCC tissues and 340 normal esophageal mucosal epithelium tissues were obtained from The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) cohorts. The transcriptome expression datasets of 19 esophageal squamous carcinoma cell lines were downloaded from Cancer Cell Line Encyclopedia (CCLE). The R software Limma package was used to identify 6,231 differentially expressed genes and 647 differentially expressed immune-related genes between normal and ESCC tissues. Gene functional analysis was performed using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). Weighted gene co-expression network analysis (WGCNA) was used to screen out 18 immune-related prognostic genes. We then established the prognostic and risk signature using these genes, and the patients were divided into low-risk and high-risk groups. Compared with high-risk group patients, the low-risk group patients had longer overall survival. M1 macrophages and resting dendritic cells were differentially distributed between the low-risk and high-risk groups and were related to patient survival. We also examined the functional immune cell and immune molecule levels in low-risk and high-risk group patients, with significant differences in the tumor microenvironment between the two groups. To further verify the accuracy of the prognostic risk model, we performed area under the ROC curve (AUC) analysis. The AUC value was 0.931 for the prognostic risk, which was better than the microsatellite instability (MSI) and Tumor Immune Dysfunction and Exclusion (TIDE) scores. In conclusion, we found 18 immune-related prognostic genes related to the occurrence of ESCC and established a prognostic model for predicting disease severity.