IntroductionCuproptosis is a novel identified regulated cell death (RCD), which is correlated with the development, treatment response and prognosis of cancer. However, the potential role of cuproptosis-related genes (CRGs) in the tumor microenvironment (TME) of gastric cancer (GC) remains unknown.MethodsTranscriptome profiling, somatic mutation, somatic copy number alteration and clinical data of GC samples were downloaded from the Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) database to describe the alterations of CRGs from genetic and transcriptional fields. Differential, survival and univariate cox regression analyses of CRGs were carried out to investigate the role of CRGs in GC. Cuproptosis molecular subtypes were identified by using consensus unsupervised clustering analysis based on the expression profiles of CRGs, and further analyzed by GO and KEGG gene set variation analyses (GSVA). Genes in distinct molecular subtypes were also analyzed by GO and KEGG gene enrichment analyses (GSEA). Differentially expressed genes (DEGs) were screened out from distinct molecular subtypes and further analyzed by GO enrichment analysis and univariate cox regression analysis. Consensus clustering analysis of prognostic DEGs was performed to identify genomic subtypes. Next, patients were randomly categorized into the training and testing group at a ratio of 1:1. CRG Risk scoring system was constructed through logistic least absolute shrinkage and selection operator (LASSO) cox regression analysis, univariate and multivariate cox analyses in the training group and validated in the testing and combined groups. Real-time quantitative polymerase chain reaction (RT-qPCR) was used to evaluate the expression of key Risk scoring genes. Sensitivity and specificity of Risk scoring system were examined by using receiver operating characteristic (ROC) curves. pRRophetic package in R was used to investigate the therapeutic effects of drugs in high- and low- risk score group. Finally, the nomogram scoring system was developed to predict patients’ survival through incorporating the clinicopathological features and CRG Risk score.ResultsMost CRGs were up-regulated in tumor tissues and showed a relatively high mutation frequency. Survival and univariate cox regression analysis revealed that LIAS and FDX1 were significantly associated with GC patients’ survival. After consensus unsupervised clustering analysis, GC patients were classified into two cuproptosis molecular subtypes, which were significantly associated with clinical features (gender, age, grade and TNM stage), prognosis, metabolic related pathways and immune cell infiltration in TME of GC. GO enrichment analyses of 84 DEGs, obtained from distinct molecular subtypes, revealed that DEGs primarily enriched in the regulation of metabolism and intracellular/extracellular structure in GC. Univariate cox regression analysis of 84 DEGs further screened out 32 prognostic DEGs. According to the expression profiles of 32 prognostic DEGs, patients were re-classified into two gene subtypes, which were significantly associated with patients’ age, grade, T and N stage, and survival of patients. Nest, the Risk score system was constructed with moderate sensitivity and specificity. A high CRG Risk score, characterized by decreased microsatellite instability-high (MSI-H), tumor mutation burden (TMB) and cancer stem cell (CSC) index, and high stromal and immune score in TME, indicated poor survival. Four of five key Risk scoring genes expression were dysregulated in tumor compared with normal samples. Moreover, CRG Risk score was greatly related with sensitivity of multiple drugs. Finally, we established a highly accurate nomogram for promoting the clinical applicability of the CRG Risk scoring system.DiscussionOur comprehensive analysis of CRGs in GC demonstrated their potential roles in TME, clinicopathological features, and prognosis. These findings may improve our understanding of CRGs in GC and provide new perceptions for doctors to predict prognosis and develop more effective and personalized therapy strategies.
The development of cancer was determined by not only the intrinsic properties of cancer cells, but also the communication between cancer cells and tumor microenvironment (TME). We applied ESTIMATE and CIBERSORT algorithms to calculate the immune/stromal component and tumor-infiltrating immune cells (TICs) in TME of BC. The results showed that immune component in TME predicted patients’ survival and associated with progression of BC. Differentially expressed genes (DEGs) were primarily enriched in immune-related activities. Finally, CCL19 was acquired which shared the leading nodes in PPI network and was associated with patients’ survival. High expression of CCL19 predicted better prognosis and participated in progression of BC. Genes in CCL19 up-regulated group were enriched in immune-related activities and these functions might depend on the communications between CCL19 and multiple TICs in TIME. In conclusion, CCL19 functioned as a potential prognostic biomarker and a modulator of TIME in BC through communicating with various TICs.
BackgroundTumor immune microenvironment (TIME) plays a significant role in the initiation and progression of bladder urothelial carcinoma (BLCA). However, there are only a few researches regarding the association between immune-related genes and tumor-infiltrating immune cells (TICs) in TIME of BLCA.MethodsWe calculated the proportion of immune/stromal component and TICs of 414 BLCA samples and 19 normal samples downloaded from TCGA database with the help of ESTIMATE and CIBERSORT algorithms. Differentially expressed genes (DEGs) were obtained from the comparison between Stromal and Immune Score and further analyzed by GO and KEGG enrichment analysis, as well as PPI network and COX regression analysis. CXCL12 was overlapping among the above analyses. Single gene analysis of CXCL12 was carried out through difference analysis, paired analysis and GSEA. The association between CXCL12 and TICs was assessed by difference analysis and correlation analysis.ResultsImmune and stromal component in TIME of BLCA were associated with patients’ clinicopathological characteristics. 284 DEGs were primarily enriched in immune-associated activities, among which CXCL12 was the most significant gene sharing the leading nodes in PPI network and being closely related with patients’ survival. Single gene analysis and immunohistochemistry revealed that CXCL12 was down-regulated in BLCA samples and significantly related with the clinicopathological characteristics of patients. Further analysis suggested that CXCL12 was involved in the immune-associated activities probably through its close cross-talk with TICs.ConclusionsCXCL12 down-regulation could be a potential biomarker to predict the unbalanced immune status of TIME of BLCA, which might provide an extra insight for the immunotherapy of BLCA.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.