Construction of mRNA prognosis signature associated with differentially expressed genes in early stage of stomach adenocarcinomas based on TCGA and GEO datasets.

Stomach adenocarcinomas (STAD) are the most common malignancy of the human digestive system and represent the fourth leading cause of cancer-related deaths. As early-stage STAD are generally mild or asymptomatic, patients with advanced STAD have short overall survival. Early diagnosis of STAD has a considerable influence on clinical outcomes. The mRNA expression data and clinical indicators of STAD and normal tissues were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database. The gene expression differences were analyzed by R packages, and gene function enrichment analysis was performed. Kaplan-Meier method and univariate Cox proportional risk regression analysis were used to screen differential expressed genes (DEGs) related to survival of STAD patients. Multivariate Cox proportional risk regression analysis was used to further screen and determine the prognostic DEGs in STAD patients, and to construct a multigene prognostic prediction signature. The accuracy of predictive signature was tested by receiver operating characteristic (ROC) curve software package, and the nomogram of patients with STAD was drawn. Cox regression was used to investigate the correlation between multigene prognostic signature and clinical factors. The predictive performance of this model was compared with two other models proposed in previous studies using KM survival analysis, ROC curve analysis, Harrell consistency index and decision curve analysis (DCA). qRT-PCR and Western blot were used to verify the expression levels of prognostic genes. The pathways and functions of possible involvement of features were predicted using the GSEA method. A total of 569 early-stage specific DEGs were retrieved from TCGA-STAD dataset, including 229 up-regulated genes and 340 down-regulated genes. Enrichment analysis showed that the early-stage specific DEGs were associated with cytokine-cytokine receptor interaction, neuroactive ligand-receptor interaction, and calcium signaling pathway. Multiple Cox regression algorithm was used to identify 10 early-stage specific DEGs associated with overall survival (P < 0.01) of STAD patients, and a multi-mRNA prognosis signature was established. The patients were divided into high-risk group and low-risk group according to the risk score. In the training set, the prognostic signature was positively correlated with tumor size and stage (P < 0.05), survival curve (P < 0.001) and time-dependent ROC (AUC = 0.625). In the training dataset and test dataset, the both signatures had good predictive efficiencies. Cox regression and DCA analysis revealed that the prognostic signature was an independent factor and had a better predict effect than the conventional TNM stage classification method and the earlier published biomarkers on the prognosis of STAD patients. In this study, based on the early-stage specifically expressed genes, the prognostic signature constructed through TCGA and GEO datasets may become an indicator for clinical prognosis assessment of STAD and a new strategy for targeted therapy in the future.

Jiang F ，Lin H ，Yan H ，Sun X ，Yang J ，Dong M ... -  《-》

Identification of critical prognosis signature associated with lymph node metastasis of stomach adenocarcinomas.

Lymph node metastasis (LNM) is an important factor affecting the prognosis of patients with gastric adenocarcinoma (STAD), which is the most common malignancy of the human digestive system. Current detection techniques have limited sensitivity and specificity, and there is a lack of effective biomarkers to screen for LNM. Therefore, it is critical to screen for biomarkers that predict LNM in STAD. Gene expression differential analysis (false discovery rate < 0.05, |log2Fold change| ≥1.5) was performed on 102 LNM samples, 224 non-LNM samples, and 29 normal gastric tissue samples from The Cancer Genome Atlas (TCGA) STAD dataset, and 269 LNM-specific genes (DEGs) were obtained. Enrichment analysis showed that LNM-specific genes functioned mainly in cytokine-cytokine receptor interactions, calcium signaling, and other pathways. Ten DEGs significantly associated with overall survival in STAD patients were screened by multivariate Cox regression, and an LNM-based 10-mRNA prognostic signature was established (Logrank P < 0.0001). This 10-mRNA signature was well predicted in both the TCGA training set and the Gene Expression Omnibus validation dataset (GSE84437) and was associated with survival in patients with LNM or advanced-stage STAD. Using Kaplan-Meier survival, receiver operating characteristic curve, C-index analysis, and decision curve analysis, the 10-mRNA signature was found to be a more effective predictor of prognosis in STAD patients than the other two reported models (P < 0.0005). Protein-protein interaction network and gene set enrichment analysis of the 10-mRNA signature revealed that the signature may affect the expression of multiple biological pathways and related genes. Finally, the expression levels of prognostic genes in STAD tissues and cell lines were verified using qRT-PCR, Western blot, and the Human Protein Atlas database. Taken together, the prognostic signature constructed in this study may become an indicator for clinical prognostic assessment of LNM-STAD and provide a new strategy for future targeted therapy.

Wang X ，Zhang W ，Guo Y ，Zhang Y ，Bai X ，Xie Y ... -  《World Journal of Surgical Oncology》

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

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.

Zhang Z ，Wang C ，Shi W ，Wang Z ，Fu W ... -  《Scientific Reports》

A 9‑gene expression signature to predict stage development in resectable stomach adenocarcinoma.

Stomach adenocarcinoma (STAD) is a highly heterogeneous disease and is among the leading causes of cancer-related death worldwide. At present, TNM stage remains the most effective prognostic factor for STAD. Exploring the changes in gene expression levels associated with TNM stage development may help oncologists to better understand the commonalities in the progression of STAD and may provide a new way of identifying early-stage STAD so that optimal treatment approaches can be provided. The RNA profile retrieving strategy was utilized and RNA expression profiling was performed using two large STAD microarray databases (GSE62254, n = 300; GSE15459, n = 192) from the Gene Expression Omnibus (GEO) and the RNA-seq database within the Cancer Genome Atlas (TCGA, n = 375). All sample expression information was obtained from STAD tissues after radical resection. After excluding data with insufficient staging information and lymph node number, samples were grouped into earlier-stage and later-stage. Samples in GSE62254 were randomly divided into a training group (n = 172) and a validation group (n = 86). Differentially expressed genes (DEGs) were selected based on the expression of mRNAs in the training group and the TCGA group (n = 156), and hub genes were further screened by least absolute shrinkage and selection operator (LASSO) logistic regression. Receiver operating characteristic (ROC) curves were used to evaluate the performance of the hub genes in distinguishing STAD stage in the validation group and the GSE15459 dataset. Univariate and multivariate Cox regressions were performed sequentially. 22 DEGs were commonly upregulated (n = 19) or downregulated (n = 3) in the training and TCGA datasets. Nine genes, including MYOCD, GHRL, SCRG1, TYRP1, LYPD6B, THBS4, TNFRSF17, SERPINB2, and NEBL were identified as hub genes by LASSO-logistic regression. The model achieved discrimination in the validation group (AUC = 0.704), training-validation group (AUC = 0.743), and GSE15459 dataset (AUC = 0.658), respectively. Gene Set Enrichment Analysis (GSEA) was used to identify the potential stage-development pathways, including the PI3K-Akt and Calcium signaling pathways. Univariate Cox regression indicated that the nine-gene score was a significant risk factor for overall survival (HR = 1.28, 95% CI 1.08-1.50, P = 0.003). In the multivariate Cox regression, only SCRG1 was an independent prognostic predictor of overall survival after backward stepwise elimination (HR = 1.21, 95% CI 1.11-1.32, P < 0.001). Through a series of bioinformatics and validation processes, a nine-gene signature that can distinguish STAD stage was identified. This gene signature has potential clinical application and may provide a novel approach to understanding the progression of STAD.

Liu Z ，Liu H ，Wang Y ，Li Z ... -  《BMC GASTROENTEROLOGY》

Angiogenesis-related lncRNAs predict the prognosis signature of stomach adenocarcinoma.

Stomach adenocarcinoma (STAD), which accounts for approximately 95% of gastric cancer types, is a malignancy cancer with high morbidity and mortality. Tumor angiogenesis plays important roles in the progression and pathogenesis of STAD, in which long noncoding RNAs (lncRNAs) have been verified to be crucial for angiogenesis. Our study sought to construct a prognostic signature of angiogenesis-related lncRNAs (ARLncs) to accurately predict the survival time of STAD. The RNA-sequencing dataset and corresponding clinical data of STAD were acquired from The Cancer Genome Atlas (TCGA). ARLnc sets were obtained from the Ensemble genome database and Molecular Signatures Database (MSigDB, Angiogenesis M14493, INTegrin pathway M160). A ARLnc-related prognostic signature was then constructed via univariate Cox and multivariate Cox regression analysis in the training cohort. Survival analysis and Cox regression were performed to assess the performance of the prognostic signature between low- and high-risk groups, which was validated in the validation cohort. Furthermore, a nomogram that combined the clinical pathological characteristics and risk score conducted to predict the overall survival (OS) of STAD. In addition, ARLnc-mRNA coexpression pairs were constructed with Pearson's correlation analysis and visualized to infer the functional annotation of the ARLncs by gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. The expression of four ARLncs in STAD and their correlation with the angiogenesis markers, CD34 and CD105, were also validated by RT-qPCR in a clinical cohort. A prognostic prediction signature including four ARLncs (PVT1, LINC01315, AC245041.1, and AC037198.1) was identified and constructed. The OS of patients in the high-risk group was significantly lower than that of patients in the low-risk group (p < 0.001). The values of the time-dependent area under the curve (AUC) for the ARLnc signature for 1-, 3-, and 5- year OS were 0.683, 0.739, and 0.618 in the training cohort and 0.671, 0.646, and 0.680 in the validation cohort, respectively. Univariate and multivariate Cox regression analyses indicated that the ARLnc signature was an independent prognostic factor for STAD patients (p < 0.001). Furthermore, the nomogram and calibration curve showed accurate prediction of the survival time based on the risk score. In addition, 262 mRNAs were screened for coexpression with four ARLncs, and GO analysis showed that mRNAs were mainly involved in biological processes, including angiogenesis, cell adhesion, wound healing, and extracellular matrix organization. Furthermore, correlation analysis showed that there was a positive correlation between risk score and the expression of the angiogenesis markers, CD34 and CD105, in TCGA datasets and our clinical sample cohort. Our study constructed a prognostic signature consisting of four ARLnc genes, which was closely related to the survival of STAD patients, showing high efficacy of the prognostic signature. Thus, the present study provided a novel biomarker and promising therapeutic strategy for patients with STAD.

Han C ，Zhang C ，Wang H ，Li K ，Zhao L ... -  《BMC CANCER》