- 1Liyuan Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
- 2Affiliated Hospital of Jining Medical University, Jining, China
Gastric cancer (GC) is a prevalent form of cancer worldwide and has a high death rate, with less than 40% of patients surviving for 5 years. GC demonstrates a vital characteristic of evading regulatory cell death (RCD). However, the extent to which RCD patterns are clinically significant in GC has not been well investigated. The study created a regulatory cell death index (RCDI) signature by employing 101 machine-learning algorithms. These algorithms were based on the expression files of 1292 GC patients from 6 multicenter cohorts. RCDI is a reliable and robust determinant of the likelihood of surviving in general. Furthermore, the precision of RCDI surpasses that of the 20 signatures that have been previously disclosed. The presence of RCDI signature is closely linked to immunological characteristics, such as the infiltration of immune cells, the presence of immunotherapy markers, and the activation of immune-related functions. This suggests that there is a higher level of immune activity in cases with RCDI signature. Collectively, the use of RCDI has the potential to be a strong and encouraging method for enhancing the clinical results of individual individuals with GC.
1 Introduction
Gastric cancer (GC) is a widespread malignancy with a rapidly increasing incidence each year (Thrift and El-Serag, 2020). GC is the fifth most common cancer worldwide in terms of new cases and the fourth leading cause of death globally, according to the 2020 Global Cancer Epidemiology statistics. In China, GC is the third most prevalent type of cancer, with approximately 480,000 new cases. It is responsible for 12.4% of all cancer-related fatalities in the country (Sung et al., 2021). This is directly associated with the high prevalence of GC patients in China and the significant variability in tumor biology and clinical characteristics. Forecasting the outcome of GC poses a significant obstacle in the current clinical approach to GC treatment.
Escape from cell death is a crucial trait exhibited by tumors (Tan et al., 2021). Regulatory cell death (RCD) is a form of cell death that happens when signal transduction modules are triggered to maintain the stability of the internal environment (Tang et al., 2019). Studies have found that numerous pathways involved in regulatory cell death signaling play a role in the development and advancement of GC. Studying these pathways will help to progress the development of diagnostic and treatment methods for GC.
Rayes et al. (2019) verified that the process of neutrophil-induced NETosis contributed to facilitating the spread of GC cells. Hao et al. (2017) reported that suppressing the expression of cysteine dioxygenase 1 in GC cells can replenish the level of glutathione (GSH), enhance the activity of glutathione peroxidase 4 (GPX4), hinder the generation of reactive oxygen species (ROS), and impede iron-induced cell death. Moreover, elevated concentrations of arachidonic acid and adrenic acid in GC facilitate the production of polyunsaturated fatty acids (PUFA), which trigger lipid peroxidation and accelerate iron-mediated cell death (Lee et al., 2020). Guo et al. (2019) found that celastrol, which has anticancer properties, elevated the levels of phosphorylation of RIP1 and RIP3. This resulted in the necrotic death of GC cells. The studies collectively indicate that RCD plays a crucial role in both the formation and advancement of GC.
However, the majority of the previously mentioned research solely examines the influence of a specific RCD mode on GC. Currently, there is a lack of extensive understanding of the interactions among RCD patterns in GC. Additionally, there is a limited quantity of particular studies on the functional aspects of these processes in GC. To tackle these regions of low comprehension, we have developed a new metric called the regulatory cell death index (RCDI). This metric is specifically formulated to predict the efficacy and prognosis of treatment strategies for GC. We identified heterogeneity in patients with GC and assessed their clinical prospects based on the RCDI signature, which provides valuable guidance for selecting the most effective treatment. The specific process of the study is shown in Figure 1.
2 Materials and methods
2.1 Data acquisition and processing
The gene expression data, consisting of 375 GC samples and 32 samples of para-cancer tissue, along with the clinical data of GC patients, were obtained from the TCGA database (https://www.tcga.org) (Cancer Genome Atlas Research Network, 2014). The bulk-seq datasets including GSE84437 (Yoon et al., 2020), GSE13861 (Cho et al., 2011), GSE15459 (Tao et al., 2011), GSE26253 (Oh et al., 2018), GSE26901 (Oh et al., 2018), and GSE57303 dataset (Qian et al., 2014) as well as scRNA-seq GSE183904 dataset (Kumar et al., 2022) were obtained from the Gene Expression Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo/).
Ultimately, the study encompassed 6 unique cohorts of patients, totaling 1,292 individuals diagnosed with GC, in order to establish the prognostic signature. The criteria for inclusion were as follows: 1) The individuals were diagnosed with GC. 2) The patients possess comprehensive survival statistics (Refer to Supplementary Table S1 for comprehensive clinical parameters). The GSE84437 cohort is designated as the training cohort. The validation cohort consists of the following datasets: GSE13861, GSE15459, GSE26253, GSE26901, and GSE57303 cohort.
PRJEB25780 dataset including the immunotherapy information and RNAseq data of 78 GC patients was obtained from a published study (Kim et al., 2018) (Refer to Supplementary Table S2 for details).
2.2 Functional enrichment analysis
Gene Ontology (GO) enrichment analysis is a widely used method for conducting large-scale investigations on the functional enrichment of biological processes (BPs), molecular functions (MFs), and cell components (CCs). KEGG enrichment analysis, a widely used bioinformatics technique, is employed to examine gene pathways and functional enrichment in a certain gene set. The above enrichment analysis was performed based on the “clusterProfiler” R package.
2.3 Consensus clustering
The “ConsensusClusterPlus” R package (Wilkerson and Hayes, 2010) was employed to find the cluster of RCD-related genes based on the expression data of RCD-related genes in the GSE84437 cohort and all samples. Subsequently, the optimal number of clusters was determined using the consensus score matrix, CDF curve, PAC score, and Nbclust.
2.4 Integration of machine learning algorithms
To enhance the precision and consistency of the RCDI signature, the study incorporated ten machine-learning algorithms into our analysis. These algorithms encompass Lasso (Friedman et al., 2010), CoxBoost (Binder et al., 2009), random survival forest (RSF) (Rigatti, 2017), elastic network (Enet) (Friedman et al., 2010), Ridge (Friedman et al., 2010), Stepwise Cox (Núñez et al., 2011), partial least squares regression for Cox (plsRcox) (Lê Cao et al., 2011), supervised principal components (SuperPC) (Bair and Tibshirani, 2004), generalized boosted regression modeling (GBM) (Guo and Chang, 2022), and survival support vector machine (survival-SVM) (Van Belle et al., 2011). Among these algorithms, Lasso, stepwise Cox, CoxBoost, and RSF have shown feature selection capabilities. Therefore, we integrated these algorithms to produce a consensus model. A total of 101 algorithm combinations were performed to construct prediction models using the 10-fold cross-validation technique.
2.5 Collection of biomarkers in cancer immunotherapy
The relationship between the RCDI signature and immune cell infiltration in tumor immune microenvironment (TIME) was investigated based on the TIMER algorithm (Li et al., 2017), CIBERSORT algorithm (Newman et al., 2019), quantiseq algorithm (Finotello and Trajanoski, 2018), MCPcounter algorithm (Becht et al., 2016), and EPIC algorithm (Racle et al., 2017).
Furthermore, seven published immunotherapeutic biomarkers were enrolled. The “easier” package (Lapuente-Santana et al., 2021) was used to calculate Cytotoxic activity (CYT) (Rooney et al., 2015), IFNy signature (IFNy) (Ayers et al., 2017), Roh immune score (Roh_IS) (Roh et al., 2017), chemokine signature (chemokines) (Messina et al., 2012), Davoli immune signature (Davoli_IS) (Cabrita et al., 2020), extended immune signature (Ayers_expIS) (Ayers et al., 2017). TIDE scores were retrieved from the TIDE database (http://tide.dfci.harvard.edu/).
2.6 Drug sensitivity analysis
The “pRRophetic” R package was applied to predict the therapeutic response of GC patients to common drugs, and the value of the RCDI signature in guiding the selection of drugs for GC patients was assessed based on the IC50 values in different RCDI score groupings.
2.7 Single-cell analysis
The scRNA-seq data of GSE183904 were analyzed using the R package “Seurat” (4.0.3). Cells with <300 genes, or >20% mitochondrial genes were excluded. Subsequently, this study performed Principal Component Analysis (PCA) using the first 1,500 highly variable genes. The first 15 principal components were then chosen to construct the t-Distributed Stochastic Neighbor Embedding (t-SNE) plot. Clustree graphically represents the Seurat resolution concerning cluster clusters, providing guidance on the number of cell clusters and selecting the most suitable level of resolution. The “FindAllMarkers” function retrieves the hallmark genes associated with each cell cluster.
Cell clusters were annotated using reference data from the Human Cell Atlas and were subsequently refined based on specific cell biomarkers including Epithelial cells (EPCAM, KRT18, and MUC1) Endothelial cells (VWF) Fibroblasts (LUM) Plasma cells (MZB1) NK/T cells (CD2 and NKG7) B cells (CD79A and MS4A1) Macro/Mono (CD14 and CD68) and Mast cells (MS4A2).
The “infercnv” R package was utilized to analyze the copy number variation (CNV) of epithelial cells. Predictions were made about the existence of malignant epithelial cells and normal epithelial cells based on the CNV score. The “irGSEA” R package was employed to conduct gene set enrichment analysis (GSEA) on the scRNAseq dataset using the “AUCell”, “UCell”, “singscore”, and “ssgsea” algorithms.
2.8 The protein expressions of prognostic genes
The protein expression of the prognostic hub genes between GC and normal tissues was validated using immunohistochemistry (IHC) provided by the Human Protein Atlas database (HPA, https://www.proteinatlas.org/).
2.9 Statistical analysis
The R package “limma” was employed to extract differentially expressed genes (DEGs). Statistical differences between groups were determined by Student’s t-test for normally distributed variables, and for non-normally distributed variables, statistical differences between groups were determined by the Wilcoxon test. The statistical studies were conducted using the R project (version 4.3.3).
3 Results
3.1 Variant landscape of RCD-related genes in GC patients
In this study, we collected a total of 20 RCD patterns and 2013 key regulatory genes (Refer to Supplementary Table S3 for details) from the existing published articles (Qin et al., 2023; Liu et al., 2023a; Liu et al., 2023b), MSigDB Database (http://software.broadinstitute.org/gsea/msigdb/index.jsp), KEGG database (https://www.genome.jp/kegg), and Gene cards database (https://www.genecards.org/). We removed 449 duplicate gene symbols, resulting in 1564 RCD-related genes for subsequent analysis. (Figure 2A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1447970/fphar-15-1447970-HTML-r1/image_m/fphar-15-1447970-g002.jpg)
Figure 2. Variant landscape of RCD-related genes in GC patients. (A) The upset plot displaying diverse RCD patterns and key regulatory genes. (B) Volcano plot of the RCD-related DEGs. Points with labels are obvious DEGs which adjusted. p-value < 0.05 and |log2FC| > 1. (C) GO enrichment analyses based on the DEGs. (D) KEGG enrichment analyses based on the DEGs. (E) An oncoplot of RCD-related genes in the TCGA-STAD cohort. (F) CNV values of RCD-related genes in the TCGA-STAD cohort.
A total of 285 genes with significant differential expression (adjusted p-value < 0.05 and |log2FC| > 1) were identified in the TCGA-STAD cohort (Figure 2B). Furthermore, the DEGs are associated with many cell death pathways and signal transduction pathways linked to cancer, as demonstrated by the KEGG and GO enrichment studies (Figures 2C,D). The TCGA-STAD cohort was used to assess the variation in RCD-related genes.
The findings indicated that approximately 85.14% (361 out of 424) of individuals with GC exhibited genetic alterations. Figure 2E displayed the top 20 mutations in RCD-related genes, with TP53 exhibiting the highest mutation frequency at 46%. The examination of CNV status revealed frequent alterations in RCD-related genes. Analysis revealed that WWOX had the highest degree of CNV deletion, whereas MYC displayed the most pronounced CNV amplification (Figure 2F).
3.2 Identification of RCD-related genes associated with GC
The “ConsensusClusterPlus” R package was employed to conduct consensus cluster analysis on a set of 1,564 RCD-related genes. This study utilized a uniform clustering approach to divide the GC data into k clusters, with k ranging from 2 to 7. The cumulative distribution function (CDF) curve of the consensus score matrix and the proportion of ambiguous clustering (PAC) statistics for the fuzzy clustering ratio suggest that the optimal number of clusters is attained when k = 2 (Figures 3A,B). Two distinct subcategories of RCD, namely, Cluster A and Cluster B, were identified, (Figure 3C). Figure 3D illustrates the heterogeneity in RCD-related pathway activity between two clusters based on the ssGSEA algorithm. A total of 859 RCD-related genes that showed differential expression were identified when comparing the two RCD subtypes (p-value < 0.05) (Figure 3E).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1447970/fphar-15-1447970-HTML-r1/image_m/fphar-15-1447970-g003.jpg)
Figure 3. Identification of RCD-related genes associated with GC. (A) Consensus clustering model with cumulative distribution function (CDF) for k = 2–7 (k means cluster count). (B) Relative change in the area under the CDF curve for k = 2–7. (C) The consensus score matrix of all samples when k = 2. A higher consensus score between two samples indicates they are more likely to be grouped into the same cluster in different iterations. (D) The heatmap displaying the heterogeneity in RCD-related pathway activity between two clusters based on the ssGSEA algorithm. (E) The volcano plot of the RCD-related DEGs. Points with labels are obvious DEGs with adjusted. p-value < 0.05.
3.3 Construction and validation of the RCDI signature
The best-performing predictive signature was determined as the signature with the greatest mean C-index in five external validation cohorts, due to overfitting in the training cohort (Figure 4A). The findings indicated that the Lasso + RSF algorithm combination demonstrated the highest average C-index (0.647), making it the optimal combination of algorithms for developing the RCDI signature.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1447970/fphar-15-1447970-HTML-r1/image_m/fphar-15-1447970-g004.jpg)
Figure 4. Construction and validation of the RCDI signature. (A) A total of 101 combinations of machine learning algorithms for the RCDI signatures via a 10-fold cross-validation framework based on the GSE84437 cohort. The C-index of each model was calculated across validation datasets, including GSE13861, GSE15459, GSE26253, GSE26901, and GSE57303 cohort. (B) The coefficients of the most useful prognostic genes based on the lasso algorithm in the GSE84437 cohort. (C) The importance of the 17 most valuable genes based on the RSF algorithm in the GSE84437 cohort. (D) Kaplan-Meier survival curve of OS between patients in high-RCDI score group and low-RCDI score group in GSE13861, GSE15459, GSE26253, GSE26901, GSE57303, and GSE84437 cohort. (E) ROC analysis of RCDI in GSE13861, GSE15459, GSE26253, GSE26901, GSE57303, and GSE84437 cohort.
Lasso algorithm was used to screen out the most valuable genes (Figure 4B). RSF algorithm was further used to filtrate the most reliable model (Figure 4C). Log-rank score test for splitting survival trees was conducted. First, the x-variable x was assumed to be ordered as
where
where
Afterwards, the RCDI score for each sample was calculated. Kaplan-Meier analysis and assessment of prognostic performance were conducted. Demonstrating a substantial difference in survival time between the low-RCDI score and high-RCDI score groups in all six cohorts (Figure 4D). The accuracy and reliability of the RCDI signature predicting 1-, 2-, 3-, 4- and 5-year survival of GC patients was supported by empirical evidence that the area under the curve (AUC) values exceeded 0.65 in multiple distinct cohorts (Figure 4E).
3.4 The relationship between RCDI signature and TIME characteristics
To assess the role of RCDI signature in GC TIME, we assessed the relationship between RCDI signature and immune infiltrating cells (Figure 5A). Based on the TIMER algorithm, CIBERSORT algorithm, quantiseq algorithm, MCPcounter, and EPIC algorithm, the RCDI signature was correlated with the majority of tumor immune infiltrating cells. The immune heterogeneity between the high-RCDI score group and the low-RCDI score group was demonstrated.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1447970/fphar-15-1447970-HTML-r1/image_m/fphar-15-1447970-g005.jpg)
Figure 5. The relationship between RCDI signature and TIME characteristics. (A) Heatmap displaying the correlation between the RCDI signature and immune infiltrating cells. (B) Box plot displaying the correlation between the RCDI signature and immunotherapy response in the immunotherapy cohort. (C) ROC curves of the RCDI signature to predict the benefits of immunotherapy in the immunotherapy cohort. (D) Box plot displaying the correlation between the RCDI signature and immune modulators. (E) Heatmap displaying the correlation between the RCDI signature and 13 immune-related processes.
We also evaluated the predictive value of RCDI signature in immunotherapy using immunotherapy-related datasets. We found that in the PRJEB25780 cohort, The RCDI score of responders was significantly lower than that of non-responders (Figure 5B). The receiver operating characteristic (ROC) analysis showed that the RCDI signature exhibited a superior ability to predict the efficacy of immunotherapy-based treatment (Figure 5C).
In addition, the study also evaluated the relationship between RCDI signature and known immune modulators (CYT, IFNy, Davoli_IS, Roh_IS, Ayers_expIS, chemokines, and TIDE) (Figure 5D). The values of most of the immune modulators (CYT, IFNy, Davoli_IS, Roh_IS, Ayers_expIS, and chemokines) were significantly higher in the low RCDI scores group. Meanwhile, the TIDE score was significantly lower in the low RCDI scores group. A low TIDE score indicates less likelihood of immune escape and a better response to immunotherapy. Based on the ssGSEA algorithm, the RCDI signature was significantly correlated with most immune-related processes (Figure 5E).
These findings indicated that GC patients with lower RCDI scores may experience more favorable outcomes from immunotherapy treatment.
3.5 Assessment and clinical application of the RCDI signature
Subsequently, we further performed univariate Cox regression for each characteristic of all cohorts. The results suggest that the RCDI signature may be utilized as an independent prognostic indication of unfavorable outcomes (Figure 6A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1447970/fphar-15-1447970-HTML-r1/image_m/fphar-15-1447970-g006.jpg)
Figure 6. Assessment and clinical application of the RCDI signature. (A) Forest plot displaying associations between RCDI signature and other clinical features and the prognosis of GC based on univariate Cox regression analysis. (B) C-index comparison of RCDI signature and 20 previously published signatures. (C) Box plot displaying the correlation between the RCDI signature and chemotherapy drugs. (D) A nomogram was established to predict the prognostic of GC patients based on age, gender, and RCDI score.
Machine learning-based gene expression analysis can be utilized to forecast the onset of diseases, aiding in the early detection of diseases and facilitating research on novel treatments. Over the past few years, there has been a proliferation of disease prediction models related to GC. We searched literature about GC that focused on disease prediction models. After excluding papers that had an unclear formula for the prediction model and lacked gene expression data in both the training and validation groups, a total of 19 prediction models for GC-related diseases were collected (Supplementary Table S4). These traits encompass a range of metabolic and cell death mechanisms, such as cuproptosis, ferroptosis, and autophagy. The C-index of each signature in both the training cohort and the validation cohorts is computed and then compared with the C-index of the RCDI signature. The superiority of the RCDI signature is evident in comparison to the majority of signatures in each cohort (Figure 6B).
In addition, the drug sensitivity analysis revealed that GC patients with high RCDI scores saw a significant increase in their sensitivity to commonly used chemotherapeutic medicines for GC (Figure 6C). It is suggested that the RCDI signature has a potential guiding effect on the treatment of GC patients. Ultimately, to facilitate clinical application, a nomogram was created, integrating the factors of age, gender, and RCDI score (Figure 6D).
3.6 Dissection of tumor microenvironment based on RCD patterns
To identify the optimal resolution for unsupervised clustering, the effectiveness of 15 distinct resolution values was assessed by employing “clustree” R package (Supplementary Figure S1A). The authors selected a resolution of 0.5 for the initial distinction of cell kinds, as indicated by the pre-assigned notes. The tSNE plot identified 17 separate cell clusters, with each cluster being assigned a unique color (Supplementary Figure S1B). A total of 8 known cell types were identified based on specific cell marker genes (Supplementary Figures S1C,D).
Then, epithelial cells were isolated and subsequently reaggregated into 14 clusters using the “Seurat” and “clustree” R packages (Supplementary Figures S2A,B). In order to identify malignant epithelial cells in GC, the “infercnv” R package was employed to analyze the CNV in each epithelial cell on a wide scale and quantify it as the CNV score, with fibroblasts and endothelial cells as reference (Supplementary Figure S2C). The results showed that the CNV scores of clusters 0, 1, 3–6, 8–9, and 12–13 were significantly higher than the reference cells and there was no significant difference in CNV scores between clusters 2, 7, 10, and 11 with reference cells (Supplementary Figure S2D). Therefore, clusters 0, 1, 3–6, 8–9, and 12–13 are categorized as malignant epithelial cells, whereas clusters 2, 7, 10, and 11 are categorized as normal epithelial cells (Supplementary Figure S2E). To summarize, a total of 9 unique cell types were effectively recognized and recorded (Figure 7A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1447970/fphar-15-1447970-HTML-r1/image_m/fphar-15-1447970-g007.jpg)
Figure 7. Heterogeneity of multiple RCD patterns in cell subsets in GC patients. (A) The t-SNE plot displaying the composition of cells in the microenvironment of GC. (B) The boxplot displaying the RCDI score in different cell types. (C) Histogram displaying the number of RCD patterns in different cell types. (D) The heatmap displaying the RCD patterns in different cell types.
Furthermore, we explored the detailed distribution of RCDI scores in GC patients using single-cell RNA transcriptome data. We consistently found that malignant cells had a higher RCDI score than other cells (Figure 7B). This further explains why the higher RCDI score group had a worse prognosis.
Ultimately, the study preliminarily described the differences in RCD-related pathway activity of various cell subtypes in the GC microenvironment (Figures 7C,D).
4 Discussion
GC is a highly aggressive tumor characterized by a high rate of invasion, recurrence, and a bleak prognosis. In order to enhance the survival rate, it is necessary to conduct early screening of the factors associated with poor prognosis and provide tailored individual treatment. However, to achieve this, there is a need for more sensitive and accurate biomarkers. The current categorization of GC primarily relies on the histological features of tumor cells. Furthermore, numerous studies have demonstrated that molecular markers can also uncover clinical significance and predictive worth. As our knowledge of GC biology grows, researchers are exploring emerging prognostic criteria and novel molecular markers to improve prognosis accuracy and tailor treatment strategies.
The study employed scRNAseq and bulk-RNAseq to characterize the pathway activity profiling of RCD patterns in GC patients. In addition, 101 distinct machine-learning algorithm combinations were utilized to develop a stable RCDI signature, which was derived by analyzing comprehensive bulk-RNAseq datasets. The stability and reliability of the RCDI signature are ensured by utilizing the advantages of each algorithm and adopting the set learning technique. By conducting validation on several datasets, the signature demonstrated exceptional performance in forecasting the prognosis of GC patients.
Furthermore, the majority of genes enrolled in the RCDI signature have been confirmed to be involved in the progression of GC. The research of Yanrong et al. (Cui et al. 2020a) study suggested that NOTCH3 is a prognostic factor correlating with immune tolerance in GC. Furthermore, as a crucial target of miR-491-5p/miR-875-5p, NOTCH3 has been shown to promote gastric carcinogenesis by upregulating PHLDB2 expression and activating the Akt pathway (Kang et al., 2021). Song and Zhou (2021) discovered that HOXA10 facilitates the process of epithelial-mesenchymal transition, which contributes to the spread of GC. This is achieved, in part, through the regulation of the TGFB2/Smad/METTL3 signal transduction pathway. Zhang et al. (2022a) verified that the increased expression of PDK4 promoted the proliferation, migration, and invasion capacity of GC cells. Besides, Buckley et al. (2022) discovered that individuals with recessive mutations in the BRCA2 gene have a higher likelihood of developing GC. This finding suggests that GC could be considered part of the whole spectrum of cancer risks associated with the BRCA1/2 genes. Tang et al. (Buckley et al. 2022) reported that NOX4-driven ROS formation regulates the proliferation and apoptosis of GC cells through the GLI1 pathway. Furthermore, Tu et al. (2023) demonstrated that HTR2B Regulates Lipid Metabolism to Inhibit Ferroptosis in GC. Wang et al. (2022) reported that APAF1-binding long noncoding RNA promotes tumor growth and multidrug resistance in GC by blocking apoptosome assembly. Kwon et al. (2011) reported the frequent upregulation of LAMC2 in GC due to promoter demethylation. Similarly, Ii et al. (2011) discussed the co-expression of Laminin [β3 and γ2] chains and the epigenetic inactivation of the Laminin [α3] chain in GC. Peng et al. (37016377) reported FYN/TOPK/HSPB1 axis facilitates the proliferation and metastasis of GC. The research of Xu et al. (2020) suggested the downregulation of GLS2 has been linked to its role as a tumor suppressor gene in GC. In conclusion, these studies suggest that programmed cell death plays a critical role in the development and progression of GC, and RCDI signature could serve as a biomarker for assessing GC. Moreover, multiple studies have demonstrated that measuring the levels of TUBB3 in the serum can assist in determining the appropriate chemotherapy drugs for patients with advanced gastric cancer (Yu et al., 2012; Luo et al., 2015; Huang et al., 2013; Di Bartolomeo et al., 2021). TICAM1, also known as TRIF, is a protein involved in the TRIF-IFN-I pathway, which plays a role in Helicobacter-induced gastric cancer (Bali et al., 2024). Additionally, TICAM1 has been identified as a promising target for immune therapy in gastric cancer (Cui et al., 2020b). Zhang et al. (2022b) demonstrated that CD19+ CD24hi CD38hi regulatory B cells were higher significantly in patients with gastric cancer than in the healthy group. Ci et al. (2024) reported that AP1G1 has been shown to promote gastric cancer stem cell (GCSC) characteristics, indicating its involvement in shaping the properties of these cancer cells. Additionally, AP1G1 has been implicated in the inhibition of gastric cancer growth and tumor development, suggesting its potential as a therapeutic target for gastric cancer treatment (Ma et al., 2021). Furthermore, the downregulation of AP1G1 expression has been associated with the suppression of apoptosis in gastric cancer cells, further emphasizing the significance of AP1G1 in gastric cancer progression (Ci et al., 2024). Moreover, microRNAs have been implicated in gastric cancer, with studies showing the expression of ELL3 was associated with exosomal and non-exosomal microRNAs (Feng et al., 2019). Chen et al. (2020) demonstrated that PTPN22 was associated with the T stage and pathological grade of STAD. Zhang et al. (2003) delved into the roles of TPA, also known as PLAT, in the apoptosis of gastric cancer cells. They found that TPA inhibits PKB activity and causes its degradation, leading to apoptosis in gastric cancer cells.
Furthermore, it has been verified that NOX4 (Tang et al., 2018) and TGFB2 (Zhang and Li, 2016) exhibit elevated expression levels in tumor tissues relative to normal tissues. Furthermore, we obtained immunohistochemistry (IHC) staining images from GC and healthy stomach tissue associated with remaining RCD genes from the HPA database. The results indicated that there were variations in the protein expression levels of the remaining RCD genes between GC and healthy stomach tissues (Supplementary Figure S3; Supplementary Table S5).
The tumor immune microenvironment is a critical factor in predicting the responsiveness of GC to immunotherapy and evaluating its prognosis. The study conducted a comprehensive examination of immune infiltration and discovered that GC patients in the high-RCDI score cohort displayed a notable abundance of immune cells, such as NK cells, macrophages, and dendritic cells. Moreover, the association between immunotherapy indicators and immunotherapy datasets, along with the RCDI signature, suggests that persons with lower RCDI scores in GC are likely to have better immunotherapy outcomes. Our research findings suggest that RCDI has the potential to be a valuable biomarker for predicting genomic patterns and evaluating the efficacy of immunotherapy in patients with GC.
A significant limitation of this study is the lack of in vitro or in vivo experiments to directly confirm our results. Although bioinformatics analysis and computational methodologies are widely used, experimental validation remains an essential aspect of scientific research. By performing experiments, we can get valuable insights into the functional implications of observed patterns and improve the trustworthiness of our findings. Therefore, future research must give priority to conducting targeted experiments to verify and build upon the results of RCDI.
5 Conclusion
In general, our findings suggest that the RCDI signature can serve as a valuable tool for guiding treatment decisions and improving patient outcomes.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
JS: Conceptualization, Methodology, Software, Supervision, Validation, Writing–original draft. LR: Software, Validation, Writing–original draft. SZ: Methodology, Software, Writing–original draft. YZ: Funding acquisition, Investigation, Writing–review and editing. YS: Methodology, Software, Supervision, Visualization, Writing–review and editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Acknowledgments
All authors would like to express our sincere thanks for sharing the online databases.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2024.1447970/full#supplementary-material
References
Ayers, M., Lunceford, J., Nebozhyn, M., Murphy, E., Loboda, A., Kaufman, D. R., et al. (2017). IFN-γ–related mRNA profile predicts clinical response to PD-1 blockade. J. Clin. Investigation 127 (8), 2930–2940. doi:10.1172/JCI91190
Bair, E., and Tibshirani, R. (2004). Semi-supervised methods to predict patient survival from gene expression data. PLoS Biol. 2 (4), e108. doi:10.1371/journal.pbio.0020108
Bali, P., Lozano-Pope, I., Hernandez, J., Estrada, M. V., Corr, M., Turner, M. A., et al. (2024). TRIF-IFN-I pathway in Helicobacter-induced gastric cancer in an accelerated murine disease model and patient biopsies. iScience 27 (4), 109457. doi:10.1016/j.isci.2024.109457
Becht, E., Giraldo, N. A., Lacroix, L., Buttard, B., Elarouci, N., Petitprez, F., et al. (2016). Estimating the population abundance of tissue-infiltrating immune and stromal cell populations using gene expression. Genome Biol. 17, 218. doi:10.1186/s13059-016-1070-5
Binder, H., Allignol, A., Schumacher, M., and Beyersmann, J. (2009). Boosting for high-dimensional time-to-event data with competing risks. Bioinformatics 25 (7), 890–896. doi:10.1093/bioinformatics/btp088
Buckley, K. H., Niccum, B. A., Maxwell, K. N., and Katona, B. W. (2022). Gastric cancer risk and pathogenesis in BRCA1 and BRCA2 carriers. Cancers (Basel) 14 (23), 5953. doi:10.3390/cancers14235953
Cabrita, R., Lauss, M., Sanna, A., Donia, M., Skaarup Larsen, M., Mitra, S., et al. (2020). Author Correction: tertiary lymphoid structures improve immunotherapy and survival in melanoma. Nature 580 (7801), E1. doi:10.1038/s41586-020-2155-6
Cancer Genome Atlas Research Network (2014). Comprehensive molecular profiling of lung adenocarcinoma. Nature 511 (7511), 543–550. doi:10.1038/nature13385
Chen, J., Zhao, X., Yuan, Y., and Jing, J. J. (2020). The expression patterns and the diagnostic/prognostic roles of PTPN family members in digestive tract cancers. Cancer Cell. Int. 20, 238. doi:10.1186/s12935-020-01315-7
Cho, J. Y., Lim, J. Y., Cheong, J. H., Park, Y. Y., Yoon, S. L., Kim, S. M., et al. (2011). Gene expression signature-based prognostic risk score in gastric cancer. Clin. Cancer Res. 17 (7), 1850–1857. doi:10.1158/1078-0432.CCR-10-2180
Ci, Y., Zhang, Y., and Zhang, X. (2024). Methylated lncRNAs suppress apoptosis of gastric cancer stem cells via the lncRNA-miRNA/protein axis. Cell. Mol. Biol. Lett. 29 (1), 51. doi:10.1186/s11658-024-00568-8
Cui, L., Wang, X., and Zhang, D. (2020b). TLRs as a promise target along with immune checkpoint against gastric cancer. Front. Cell. Dev. Biol. 8, 611444. doi:10.3389/fcell.2020.611444
Cui, Y., Li, Q., Li, W., Wang, Y., Lv, F., Shi, X., et al. (2020a). NOTCH3 is a prognostic factor and is correlated with immune tolerance in gastric cancer. Front. Oncol. 10, 574937. doi:10.3389/fonc.2020.574937
Di Bartolomeo, M., Raimondi, A., Cecchi, F., Catenacci, D. V. T., Schwartz, S., Sellappan, S., et al. (2021). Association of high TUBB3 with resistance to adjuvant docetaxel-based chemotherapy in gastric cancer: translational study of ITACA-S. Tumori 107 (2), 150–159. doi:10.1177/0300891620930803
Feng, Y., Hang, W., Sang, Z., Li, S., Xu, W., Miao, Y., et al. (2019). Identification of exosomal and non-exosomal microRNAs associated with the drug resistance of ovarian cancer. Mol. Med. Rep. 19 (5), 3376–3392. doi:10.3892/mmr.2019.10008
Finotello, F., and Trajanoski, Z. (2018). Quantifying tumor-infiltrating immune cells from transcriptomics data. Cancer Immunol. Immunother. 67 (7), 1031–1040. doi:10.1007/s00262-018-2150-z
Friedman, J., Hastie, T., and Tibshirani, R. (2010). Regularization paths for generalized linear models via coordinate descent. J. Stat. Softw. 33 (1), 1–22. doi:10.18637/jss.v033.i01
Guo, C. Y., and Chang, K. H. (2022). A novel algorithm to estimate the significance level of a feature interaction using the extreme gradient boosting machine. Int. J. Environ. Res. Public Health 19 (4), 2338. doi:10.3390/ijerph19042338
Guo, D., Zhang, W., Yang, H., Bi, J., Xie, Y., Cheng, B., et al. (2019). Celastrol induces necroptosis and ameliorates inflammation via targeting biglycan in human gastric carcinoma. Int. J. Mol. Sci. 20 (22), 5716. doi:10.3390/ijms20225716
Hao, S., Yu, J., He, W., Huang, Q., Zhao, Y., Liang, B., et al. (2017). Cysteine dioxygenase 1 mediates erastin-induced ferroptosis in human gastric cancer cells. Neoplasia 19 (12), 1022–1032. doi:10.1016/j.neo.2017.10.005
Huang, J., Hu, H., Xie, Y., Tang, Y., Liu, W., and Zhong, M. (2013). Effect of TUBB3, TS and ERCC1 mRNA expression on chemoresponse and clinical outcome of advanced gastric cancer by multiplex branched-DNA liquid chip technology. Zhong Nan Da Xue Xue Bao Yi Xue Ban 38 (6), 582–589. doi:10.3969/j.issn.1672-7347.2013.06.006
Ii, M., Yamamoto, H., Taniguchi, H., Adachi, Y., Nakazawa, M., Ohashi, H., et al. (2011). Co-expression of laminin β3 and γ2 chains and epigenetic inactivation of laminin α3 chain in gastric cancer. Int. J. Oncol. 39 (3), 593–599. doi:10.3892/ijo.2011.1048
Kang, W., Zhang, J., Huang, T., Zhou, Y., Wong, C. C., Chan, R. C. K., et al. (2021). NOTCH3, a crucial target of miR-491-5p/miR-875-5p, promotes gastric carcinogenesis by upregulating PHLDB2 expression and activating Akt pathway. Oncogene 40 (9), 1578–1594. doi:10.1038/s41388-020-01579-3
Kim, S. T., Cristescu, R., Bass, A. J., Kim, K. M., Odegaard, J. I., Kim, K., et al. (2018). Comprehensive molecular characterization of clinical responses to PD-1 inhibition in metastatic gastric cancer. Nat. Med. 24 (9), 1449–1458. doi:10.1038/s41591-018-0101-z
Kumar, V., Ramnarayanan, K., Sundar, R., Padmanabhan, N., Srivastava, S., Koiwa, M., et al. (2022). Single-cell Atlas of lineage States, tumor microenvironment, and subtype-specific expression programs in gastric cancer. Cancer Discov. 12 (3), 670–691. doi:10.1158/2159-8290.CD-21-0683
Kwon, O. H., Park, J. L., Kim, M., Kim, J. H., Lee, H. C., Kim, H. J., et al. (2011). Aberrant up-regulation of LAMB3 and LAMC2 by promoter demethylation in gastric cancer. Biochem. Biophys. Res. Commun. 406 (4), 539–545. doi:10.1016/j.bbrc.2011.02.082
Lapuente-Santana, Ó., van Genderen, M., Hilbers, P. A. J., Finotello, F., and Eduati, F. (2021). Interpretable systems biomarkers predict response to immune-checkpoint inhibitors. Patterns (N Y) 2 (8), 100293. doi:10.1016/j.patter.2021.100293
Lê Cao, K. A., Boitard, S., and Besse, P. (2011). Sparse PLS discriminant analysis: biologically relevant feature selection and graphical displays for multiclass problems. BMC Bioinforma. 12, 253. doi:10.1186/1471-2105-12-253
Lee, J. Y., Nam, M., Son, H. Y., Hyun, K., Jang, S. Y., Kim, J. W., et al. (2020). Polyunsaturated fatty acid biosynthesis pathway determines ferroptosis sensitivity in gastric cancer. Proc. Natl. Acad. Sci. U. S. A. 117 (51), 32433–32442. doi:10.1073/pnas.2006828117
Li, T., Fan, J., Wang, B., Traugh, N., Chen, Q., Liu, J. S., et al. (2017). TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 77 (21), e108–e110. doi:10.1158/0008-5472.CAN-17-0307
Liu, J., Zhang, M., Sun, Q., Qin, X., Gao, T., Xu, Y., et al. (2023a). Construction of a novel MPT-driven necrosis-related lncRNAs signature for prognosis prediction in laryngeal squamous cell carcinoma. Environ. Sci. Pollut. Res. Int. 30 (31), 77210–77225. doi:10.1007/s11356-023-26996-1
Liu, X., Nie, L., Zhang, Y., Yan, Y., Wang, C., Colic, M., et al. (2023b). Actin cytoskeleton vulnerability to disulfide stress mediates disulfidptosis. Nat. Cell. Biol. 25 (3), 404–414. doi:10.1038/s41556-023-01091-2
Luo, Y., Li, Z., Cui, S., Shen, C., Zhao, J., Wu, M., et al. (2015). Retraction note: joint detection of ERCC1, TUBB3, and TYMS guidance selection of docetaxel, 5-fluorouracil and cisplatin (DDP) individual chemotherapy in advanced gastric cancer patients. Eur. J. Med. Res. 20 (1), 29. doi:10.1186/s40001-015-0134-4
Ma, H., Wu, J., Zhou, M., Wu, J., Wu, Z., Lin, L., et al. (2021). Inhibition of glutamine uptake improves the efficacy of cetuximab on gastric cancer. Integr. Cancer Ther. 20, 15347354211045349. doi:10.1177/15347354211045349
Messina, J. L., Fenstermacher, D. A., Eschrich, S., Qu, X., Berglund, A. E., Lloyd, M. C., et al. (2012). 12-Chemokine gene signature identifies lymph node-like structures in melanoma: potential for patient selection for immunotherapy? Sci. Rep. 2, 765. doi:10.1038/srep00765
Newman, A. M., Steen, C. B., Liu, C. L., Gentles, A. J., Chaudhuri, A. A., Scherer, F., et al. (2019). Determining cell-type abundance and expression from bulk tissues with digital cytometry. Nat. Biotechnol. 37 (7), 773–782. doi:10.1038/s41587-019-0114-2
Núñez, E., Steyerberg, E. W., and Núñez, J. (2011). Regression modeling strategies. Rev. Española Cardiol. 64 (6), 501–507. doi:10.1016/j.recesp.2011.01.019
Oh, S. C., Sohn, B. H., Cheong, J. H., Kim, S. B., Lee, J. E., Park, K. C., et al. (2018). Clinical and genomic landscape of gastric cancer with a mesenchymal phenotype. Nat. Commun. 9 (1), 1777. doi:10.1038/s41467-018-04179-8
Qian, Z., Zhu, G., Tang, L., Wang, M., Zhang, L., Fu, J., et al. (2014). Whole genome gene copy number profiling of gastric cancer identifies PAK1 and KRAS gene amplification as therapy targets. Genes Chromosomes Cancer 53 (11), 883–894. doi:10.1002/gcc.22196
Qin, H., Abulaiti, A., Maimaiti, A., Abulaiti, Z., Fan, G., Aili, Y., et al. (2023). Integrated machine learning survival framework develops a prognostic model based on inter-crosstalk definition of mitochondrial function and cell death patterns in a large multicenter cohort for lower-grade glioma. J. Transl. Med. 21 (1), 588. doi:10.1186/s12967-023-04468-x
Racle, J., de Jonge, K., Baumgaertner, P., Speiser, D. E., and Gfeller, D. (2017). Simultaneous enumeration of cancer and immune cell types from bulk tumor gene expression data. eLife 6, e26476. doi:10.7554/eLife.26476
Rayes, R. F., Mouhanna, J. G., Nicolau, I., Bourdeau, F., Giannias, B., Rousseau, S., et al. (2019). Primary tumors induce neutrophil extracellular traps with targetable metastasis promoting effects. JCI Insight 5 (16), e128008. doi:10.1172/jci.insight.128008
Roh, W., Chen, P. L., Reuben, A., Spencer, C. N., Prieto, P. A., Miller, J. P., et al. (2017). Integrated molecular analysis of tumor biopsies on sequential CTLA-4 and PD-1 blockade reveals markers of response and resistance. Sci. Transl. Med. 9 (379), eaah3560. doi:10.1126/scitranslmed.aah3560
Rooney, M. S., Shukla, S. A., Wu, C. J., Getz, G., and Hacohen, N. (2015). Molecular and genetic properties of tumors associated with local immune cytolytic activity. Cell 160 (1–2), 48–61. doi:10.1016/j.cell.2014.12.033
Song, C., and Zhou, C. (2021). HOXA10 mediates epithelial-mesenchymal transition to promote gastric cancer metastasis partly via modulation of TGFB2/Smad/METTL3 signaling axis. J. Exp. Clin. Cancer Res. 40 (1), 62. doi:10.1186/s13046-021-01859-0
Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 71 (3), 209–249. doi:10.3322/caac.21660
Tan, Y., Chen, Q., Li, X., Zeng, Z., Xiong, W., Li, G., et al. (2021). Pyroptosis: a new paradigm of cell death for fighting against cancer. J. Exp. Clin. Cancer Res. 40 (1), 153. doi:10.1186/s13046-021-01959-x
Tang, C. T., Lin, X. L., Wu, S., Liang, Q., Yang, L., Gao, Y. J., et al. (2018). NOX4-driven ROS formation regulates proliferation and apoptosis of gastric cancer cells through the GLI1 pathway. Cell. Signal 46, 52–63. doi:10.1016/j.cellsig.2018.02.007
Tang, D., Kang, R., Berghe, T. V., Vandenabeele, P., and Kroemer, G. (2019). The molecular machinery of regulated cell death. Cell. Res. 29 (5), 347–364. doi:10.1038/s41422-019-0164-5
Tao, J., Deng, N. T., Ramnarayanan, K., Huang, B., Oh, H. K., Leong, S. H., et al. (2011). CD44-SLC1A2 gene fusions in gastric cancer. Sci. Transl. Med. 3 (77), 77ra30. doi:10.1126/scitranslmed.3001423
Thrift, A. P., and El-Serag, H. B. (2020). Burden of gastric cancer. Clin. Gastroenterol. Hepatol. 18 (3), 534–542. doi:10.1016/j.cgh.2019.07.045
Tu, R. H., Wu, S. Z., Huang, Z. N., Zhong, Q., Ye, Y. H., Zheng, C. H., et al. (2023). Neurotransmitter receptor HTR2B regulates lipid metabolism to inhibit ferroptosis in gastric cancer. Cancer Res. 83 (23), 3868–3885. doi:10.1158/0008-5472.CAN-23-1012
Van Belle, V., Pelckmans, K., Van Huffel, S., and Suykens, J. A. K. (2011). Improved performance on high-dimensional survival data by application of Survival-SVM. Bioinformatics 27 (1), 87–94. doi:10.1093/bioinformatics/btq617
Wang, Q., Chen, C., Xu, X., Shu, C., Cao, C., Wang, Z., et al. (2022). APAF1-Binding long noncoding RNA promotes tumor growth and multidrug resistance in gastric cancer by blocking apoptosome assembly. Adv. Sci. (Weinh) 9 (28), e2201889. doi:10.1002/advs.202201889
Wilkerson, M. D., and Hayes, D. N. (2010). ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking. Bioinformatics 26 (12), 1572–1573. doi:10.1093/bioinformatics/btq170
Xu, L., Zhou, D., Li, F., and Ji, L. (2020). Glutaminase 2 functions as a tumor suppressor gene in gastric cancer. Transl. Cancer Res. 9 (8), 4906–4913. doi:10.21037/tcr-20-2246
Yoon, S. J., Park, J., Shin, Y., Choi, Y., Park, S. W., Kang, S. G., et al. (2020). Deconvolution of diffuse gastric cancer and the suppression of CD34 on the BALB/c nude mice model. BMC Cancer 20 (1), 314. doi:10.1186/s12885-020-06814-4
Yu, J., Gao, J., Lu, Z., Li, Y., and Shen, L. (2012). Serum levels of TUBB3 correlate with clinical outcome in Chinese patients with advanced gastric cancer receiving first-line paclitaxel plus capecitabine. Med. Oncol. 29 (5), 3029–3034. doi:10.1007/s12032-012-0292-y
Zhang, B., Wu, Q., Ye, X. F., Liu, S., Lin, X. F., and Chen, M. C. (2003). Roles of PLC-gamma2 and PKCalpha in TPA-induced apoptosis of gastric cancer cells. World J. Gastroenterol. 9 (11), 2413–2418. doi:10.3748/wjg.v9.i11.2413
Zhang, W., and Li, Y. (2016). miR-148a downregulates the expression of transforming growth factor-β2 and SMAD2 in gastric cancer. Int. J. Oncol. 48 (5), 1877–1885. doi:10.3892/ijo.2016.3437
Zhang, Y., Wu, J., Zhang, H., and Wu, C. (2022b). The regulation between CD4+CXCR5+ follicular helper T (tfh) cells and CD19+CD24hiCD38hi regulatory B (breg) cells in gastric cancer. J. Immunol. Res. 2022, 9003902. doi:10.1155/2022/9003902
Keywords: gastric cancer, regulatory cell death, immune microenvironment, prognostic signature, immunotherapy
Citation: Sun J, Rao L, Zhou S, Zeng Y and Sun Y (2024) Unraveling the regulatory cell death pathways in gastric cancer: a multi-omics study. Front. Pharmacol. 15:1447970. doi: 10.3389/fphar.2024.1447970
Received: 12 June 2024; Accepted: 26 August 2024;
Published: 09 September 2024.
Edited by:
Donatella Delle Cave, National Research Council (CNR), ItalyReviewed by:
Jung-Mao Hsu, China Medical University, TaiwanPratik Bhojnagarwala, Wistar Institute, United States
Copyright © 2024 Sun, Rao, Zhou, Zeng and Sun. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yulan Zeng, 1989LY0551@hust.edu.cn; Yalu Sun, 18487344536@163.com