Published online Mar 15, 2024. doi: 10.4251/wjgo.v16.i3.919
Peer-review started: September 7, 2023
First decision: December 5, 2023
Revised: December 16, 2023
Accepted: February 2, 2024
Article in press: February 2, 2024
Published online: March 15, 2024
Processing time: 186 Days and 18.5 Hours
Treatment options for patients with gastric cancer (GC) continue to improve, but the overall prognosis is poor. The use of PD-1 inhibitors has also brought benefits to patients with advanced GC and has gradually become the new standard treatment option at present, and there is an urgent need to identify valuable biomarkers to classify patients with different characteristics into subgroups.
To determined the effects of differentially expressed immune-related genes (DEIRGs) on the development, prognosis, tumor microenvironment (TME), and treatment response among GC patients with the expectation of providing new biomarkers for personalized treatment of GC populations.
Gene expression data and clinical pathologic information were downloaded from The Cancer Genome Atlas (TCGA), and immune-related genes (IRGs) were searched from ImmPort. DEIRGs were extracted from the intersection of the differentially-expressed genes (DEGs) and IRGs lists. The enrichment pathways of key genes were obtained by analyzing the Kyoto Encyclopedia of Genes and Genomes (KEGGs) and Gene Ontology (GO) databases. To identify genes asso
We collected 412 GC and 36 adjacent tissue samples, and identified 3627 DEGs and 1311 IRGs. A total of 482 DEIRGs were obtained. GO analysis showed that DEIRGs were mainly distributed in immunoglobulin complexes, receptor ligand activity, and signaling receptor activators. KEGG pathway analysis showed that the top three DEIRGs enrichment types were cytokine-cytokine receptors, neuroactive ligand receptor interactions, and viral protein interactions. We ultimately obtained an immune-related signature based on 10 genes, including 9 risk genes (LCN1, LEAP2, TMSB15A mRNA, DEFB126, PI15, IGHD3-16, IGLV3-22, CGB5, and GLP2R) and 1 protective gene (LGR6). Kaplan-Meier survival analysis, receiver operating characteristic curve analysis, and risk curves confirmed that the risk model had good predictive ability. Multivariate COX analysis showed that age, stage, and risk score were independent prognostic factors for patients with GC. Meanwhile, patients in the low-risk group had higher tumor mutation burden and immunophenotype, which can be used to predict the immune checkpoint inhibitor response. Both cytotoxic T lymphocyte antigen4+ and programmed death 1+ patients with lower risk scores were more sensitive to immunotherapy.
In this study a new prognostic model consisting of 10 DEIRGs was constructed based on the TME. By providing risk factor analysis and prognostic information, our risk model can provide new directions for immunotherapy in GC patients.
Core Tip: We ultimately obtained an Immune Related Signature based on 10 genes, including 9 risk genes (LCN1, LEAP2, TMSB15A mRNA, DEFB126, PI15, IGHD3-16, IGLV3-22, CGB5, GLP2R) and 1 protective gene (LGR6). Kaplan-Meier survival analysis, ROC analysis and risk curve confirmed that the risk model has good predictive ability. Multivariate COX analysis showed that age, stage and risk score were independent prognostic factors. Patients in the low-risk group had higher tumor mutation burden and immunophenotype score, which can be used to predict immune checkpoint inhibitor response. Both cytotoxic T lymphocyte antigen4+ and programmed death 1+ patients with lower risk scores were more sensitive to immunotherapy.