Prognostic model construction and immune microenvironment analysis of esophageal cancer based on gene expression data...

Accumulating evidence suggests that microRNA-target genes are closely related to tumorigenesis and progression. This study aims to screen the intersection of differentially expressed mRNAs (DEmRNAs) and the target genes of differentially expressed microRNAs (DEmiRNAs), and to construct a prognostic gene model of esophageal cancer (EC).

Gene expression, microRNA expression, somatic mutation, and clinical information data of EC from The Cancer Genome Atlas (TCGA) database were used. The intersection of DEmRNAs and the target genes of DEmiRNAs predicted by the Targetscan database and microRNA Data Integration Portal (mirDIP) database were screened. The screened genes were used to construct a prognostic model of EC. Then, the molecular and immune signatures of these genes were explored. Finally, the GSE53625 dataset from the Gene Expression Omnibus (GEO) database was further used as a validation cohort to confirm the prognostic value of the genes.

Six genes on the grounds of the intersection of DEmiRNAs target genes and DEmRNAs were identified as prognostic genes, including ARHGAP11A, H1.4, HMGB3, LRIG1, PRR11, and COL4A1. Based on the median risk score calculated for these genes, EC patients were divided into a high-risk group (n=72) and a low-risk group (n=72). Survival analysis showed that the high-risk group had a significantly shorter survival time than the low-risk group (TCGA and GEO, P<0.001). The nomogram evaluation showed high reliability in predicting the 1-year, 2-year, and 3-year survival probability of EC patients. Compared to low-risk group, higher expression level of M2 macrophages was found in high-risk group of EC patient (P<0.05), while STAT3 checkpoints showed attenuated expression level in high-risk group.

A panel of differential genes was identified as potential EC prognostic biomarkers and showed great clinical significance in EC prognosis.