Construction of store-operated calcium entry-related gene signature for predicting prognosis and indicates immune microenvironment infiltration in stomach adenocarcinomas

Abstract

Abstract Gastric adenocarcinoma (STAD) is the most prevalent malignancy of the human digestive system and the fourth leading cause of cancer-related death. Calcium pools, especially Ca2+ entry (SOCE) for storage operations, play a crucial role in maintaining intracellular and extracellular calcium balance, influencing cell activity, and facilitating tumor progression. Nevertheless, the prognostic and immunological value of SOCE in STAD has not been systematically studied. The objective of this study was to develop a risk model for SOCE signature and to explore its correlation with clinical characteristics, prognosis, tumor microenvironment (TME), as well as response to immunotherapy, chemotherapy, and targeted drugs. We used the TCGA, GEO (GSE84437 and GSE159929), cBioPortal and TIMER databases to acquire mRNA expression data for STAD, along with patient clinical indicators, single-cell sequencing data, genomic variation information, and correlations of immune cell infiltration. An analysis of SOCE genes based on tumor vs. normal tissue comparisons, pan-cancer dimension, single-cell sequencing, DNA mutation, and copy number variation revealed that SOCE genes significantly impact the survival of STAD patients and are differentially involved in the immune response. SOCE co-expressed genes were identified by Pearson analysis, and subsequently protein-protein interaction (PPI) and gene function enrichment analysis indicated that coexpressed genes were associated with multicellular pathways. Based on TCGA and GSE84437 datasets, a multifactor Cox proportional hazard regression analysis was conducted to identify SOCE co-expressed genes associated with overall survival in STAD patients. Several mRNA prognostic genes, including PTPRJ, GPR146, LTBP3, FBLN1, EFEMP2, ADAMTS7 and LBH, were identified, which could be used as effective prognostic predictors for STAD patients. In both training and test datasets, receiver operating characteristic (ROC) curves were utilized to evaluate and illustrate the predictive capability of this characteristic in forecasting overall survival of STAD patients. The qPCR and human protein atlas (HPA) were employed to assess mRNA expression and protein levels. Furthermore, the ESTIMATE, TIMER, CIBERSORT, QUANTISEQ, MCPCOUNTER, xCell and EPIC algorithms were utilized to perform immune score and analyze immune cell infiltration. It effectively revealed the difference in prognosis and immune cell infiltration in TME between high-risk and low-risk groups based on the risk signature associated with SOCE. Notably, significant differences in tumor immune dysfunction and rejection (TIDE) scores between the two groups, suggesting that patients in the low-risk group may exhibit a more favorable response to ICIS treatment. The GDSC database and R packages for predictive analysis were utilized to analyze responses to chemotherapy and targeted drugs in high-risk and low-risk groups. In summary, the 7-gene signature associated with SOCE serves as a significant biomarker for evaluating the TME and predicting the prognosis of STAD patients. In addition, it may provide valuable insights for developing effective immunotherapy and chemotherapy regiments for patients with STAD.

Keywords

• Data mining
• Prognosis
• Biological maker
• Expression difference
• Functional enrichment analysis