Next Article in Journal
Radiation in Early-Stage Breast Cancer: Moving beyond an All or Nothing Approach
Previous Article in Journal
Combining Molecular Subtypes with Multivariable Clinical Models Has the Potential to Improve Prediction of Treatment Outcomes in Prostate Cancer at Diagnosis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of a Two-Gene Signature and Establishment of a Prognostic Nomogram Predicting Overall Survival in Diffuse-Type Gastric Cancer

1
Digestive Medicine Center, The Seventh Affiliated Hospital of Sun Yat-sen University, Shenzhen 518107, China
2
Department of Thoracic Surgery, The Third Affiliated Hospital of Sun Yat-sen University, Guangzhou 510630, China
3
Division of Hematology/Oncology, Department of Pediatrics, The Seventh Affiliated Hospital of Sun Yat-sen University, Shenzhen 518107, China
4
Gastrointestinal Surgery Center, The First Affiliated Hospital of Sun Yat-sen University, Guangzhou 510080, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Curr. Oncol. 2023, 30(1), 171-183; https://doi.org/10.3390/curroncol30010014
Submission received: 28 November 2022 / Revised: 18 December 2022 / Accepted: 20 December 2022 / Published: 23 December 2022
(This article belongs to the Section Gastrointestinal Oncology)

Abstract

:
Background: It is widely acknowledged that the molecular biological characteristics of diffuse-type gastric cancer are different from intestinal-type gastric cancer. Notwithstanding that significant progress in high-throughput sequencing technology has been made, there is a paucity of effective prognostic biomarkers for diffuse gastric cancer for clinical practice. Methods: We downloaded four GEO datasets (GSE22377, GSE38749, GSE47007 and GSE62254) to establish and validate a prognostic two-gene signature for diffuse gastric cancer. The TGCA-STAD dataset was used for external validation. The optimal gene signature was established by using Cox regression analysis. Receiver operating characteristic (ROC) methodology was used to find the best prognostic model. Gene set enrichment analysis was used to analyze the possible signaling pathways of the two genes (MEF2C and TRIM15). Results: A total of four differently expressed genes (DEGs) (two upregulated and two downregulated) were identified. After a comprehensive analysis, two DEGs (MEF2C and TRIM15) were utilized to construct a prognostic model. A prognostic prediction model was constructed according to T stage, N stage, M stage and the expression of MEF2C and TRIM15. The area under the time-dependent receiver operator characteristic was used to evaluate the performance of the prognosis model in the GSE62254 dataset. Conclusions: We demonstrated that MEF2C and TRIM15 might be key genes. We also established a prognostic nomogram based on the two-gene signature that yielded a good performance for predicting overall survival in diffuse-type gastric cancer.

1. Introduction

Gastric cancer (GC) is the third most common cause of cancer-related death globally, with a low 5-year survival rate [1]. Indeed, in recent years, with the establishment of multidisciplinary team (MDT) care, GC treatment has substantially improved. Albeit multiple therapeutic approaches are available currently, it is difficult to accurately determine the optimal treatment for an individual gastric cancer (GC) patient due to clinical and genetic heterogeneity [2]. The Lauren classification is an internationally recognized histopathological classification that sorts GC into three subtypes: intestinal type, diffuse type and mixed type [3]. However, several studies have shown significant heterogeneity in biological behavior among the three subtypes of gastric cancers [4,5,6,7]. Importantly, diffuse-type GC has been documented to exhibit more aggressive behavior and a poorer prognosis, accounting for one-third of all GC patients [5]. This subtype of GC is more common in women and younger patients and is more genetic-related [8]. Nowadays, the AJCC TNM staging system is widely used to evaluate the prognosis of GC patients [9,10]. Nevertheless, multiple studies reported that the AJCC TNM staging system sometimes exhibited poor accuracy in predicting the prognosis of cancer patients, partly due to staging migration [11]. With the advancement of cancer molecular biology, cancer-specific gene signatures have been harnessed to develop new prediction tools. However, there are limited prognosis-related genes associated with diffuse-type GC. Accordingly, an effective prediction model is urgently needed to help oncologists evaluate GC patient prognosis in the clinic.
In this study, we obtained three GC datasets with Lauren classification data from the GEO database to identify differentially expressed genes (DEGs). Then, the DEGs were validated with external datasets. Univariate regression analysis and Kaplan–Meier analysis were applied to select the overall survival-related DEGs. A prognostic nomogram incorporating prognostic gene signatures and clinical survival-related characteristics was established to predict overall survival. In summary, we uncovered new gene signatures for diffuse-type GC and established a nomogram that exhibited good performance in predicting the overall survival of diffuse-type GC patients.

2. Materials and Methods

2.1. Data Collection and Processing

The gene expression profile data (GSE22377, GSE38749, GSE47007 and GSE62254) were searched and downloaded from the Gene Expression Omnibus (GEO) (https://www.ncbi.nlm.nih.gov/geo/ (accessed on 1 June 2022)) database. All included datasets met the following criteria: (1) histopathological information with detailed records of the Lauren classification available; (2) more than ten samples included; (3) gene expression profile in RNA level; (4) datasets were described in English language. Exclusion criteria included not diffuse- and intestinal-type gastric cancer; normal tissue of diffuse-type gastric cancer patients. In addition, data of 83 diffuse-type gastric cancer patients were extracted from TCGA-STAD dataset (https://portal.gdc.cancer.gov (accessed on 1 June 2022)) for external validation.
Probes were matched to the gene symbols using the annotation files provided by the manufacturer. The median expression value was calculated if multiple probes matched a single gene. Robust multi-array average (RMA)-normalized data were log2-transformed for further analysis. Limma package [12] in R/Bioconductor software (version 3.6.1, R Foundation for Statistical Computing, Vienna, Austria) was applied to screen the DEGs between diffuse-type GC and intestinal-type GC in GSE22377, GSE38749 and GSE47007 datasets. | log2Fold change (FC) | ≥ 1 and p value < 0.05 were set as the thresholds for significant differential expression.

2.2. Validation of Expression Level of DEGs

The screened DEGs were validated with a large dataset (GSE62254). Diffuse-type and intestinal-type histopathological data were available from the dataset. We analyzed the expression of DEGs between these two subtypes. Finally, a boxplot was constructed to visualize the expression level.

2.3. Kaplan–Meier Plotter Analysis

Kaplan–Meier Plotter (https://kmplot.com/analysis/ (accessed on 1 June 2022)) is an online website that integrates RNA-sequencing data of GEO dataset and survival information [13]. We used Gastric Cancer section of Kaplan–Meier Plotter (https://kmplot.com/analysis/index.php?p=service&cancer=gastric (accessed on 1 June 2022)) to explore the prognostic value of DEGs in diffuse-type gastric cancer [14].

2.4. Clinical Correlation Analysis and Biological Process Prediction

We extracted the clinical information from the GSE62254 dataset. SPSS 24.0 was used to perform a Chi-square test between every single gene and clinical pathological characteristics. A p value less than 0.05 was considered statistically significant. To understand the biological process of the identified prognostic gene signatures, gene set enrichment analysis (GSEA) was performed by using a Java GSEA desktop application (downloaded from http://www.broad.mit.edu/gsea (accessed on 1 June 2022)). The GSE62254 samples were divided into high- and low-expression groups according to the median value. The GSE62254 dataset was analyzed with the GTM file (c2.KEGG.v6.2) to identify enriched KEGG pathways. Four files containing expression datasets, gene sets, phenotype labels and chip platforms were required for running GSEA. |NES| > 1 and FDR < 0.25 were considered statistically significant.

2.5. The Establishment of the Predictive Nomogram

After testing for collinearity, prognostic gene signatures and relevant clinical parameters were included to establish a prognostic nomogram via a stepwise Cox regression model to predict the 1-, 3- and 5-year overall survival of diffuse-type gastric cancer patients in the GSE62254 dataset. A time-dependent ROC curve, Harrell’s concordance index and a calibration curve were utilized to assess the performance of the prognostic nomogram. Decision curve analysis was used to evaluate the net benefit of the program compared with TNM staging alone.
Based on the prognostic nomogram, the diffuse-type gastric cancer samples from the GSE62254 dataset were assigned to high-risk and low-risk score groups according to the median risk score. Kaplan–Meier analysis was performed to demonstrate the relationship between risk score and overall survival time by using the “survival” package. A log-rank test was used to distinguish the differences between groups.

2.6. External Validation of Two-Gene Signature-Based Nomogram by TCGA Dataset

To further confirm the prediction value of the two-gene signature nomogram, we performed ROC analysis to show the predictive performance of TNM staging and the nomogram-based model. Kaplan–Meier analysis and the log-rank test were applied to demonstrate the survival difference between the high-risk group and low-risk group. Decision curve analysis was also utilized to quantify the clinical benefits of the nomogram at different threshold probabilities. The above analyses of external validation were performed by using the TCGA-STAD dataset.

3. Results

3.1. Identification of Diffuse-Type Gastric Cancer-Specific Gene Signatures

The flowchart of the screening process used in our study to identify diffuse-type gastric cancer gene signatures is shown in Figure 1. The details of the GEO datasets included in this study are displayed in Table 1. A total of 991 (GSE22377), 166 (GSE38749) and 171 (GSE47007) DEGs were identified between the diffuse-type and intestinal-type gastric cancer datasets. Two genes (COL4A3, MEF2C) were highly expressed in diffuse-type gastric cancer, whereas two (TRIM15, MMP12) were lowly expressed in diffuse-type gastric cancer.

3.2. Validation of the Expression Level of Four Differentiated Expressed Genes

In this study, the expression level of the four DEGs identified was validated in a large dataset (GSE62254). Two upregulated and one downregulated gene were identified (Figure 2B).

3.3. Clinical Correlation Analysis of Three DEGs

Detailed clinical information of 134 patients from the GSE62254 dataset was extracted. A Chi-square test was used to evaluate the relationship between the three DEGs and clinical pathological characteristics (Table 2). In brief, the expression level of COL4A3 and TRIM15 were significantly correlated with the T stage and age, respectively, while the expression level of MEF2C was significantly correlated with age, T stage and TNM stage.

3.4. Kaplan–Meier Analysis and Evaluation of Prognostic Factors in Diffuse-Type Gastric Cancer

The survival information of COL4A3, MEF2C and TRIM15 was freely obtained in Kaplan–Meier Plotter. In this study, we assessed the difference between the expression level of the three DEGs and overall survival in diffuse-type gastric cancer. Notably, in the results of Kaplan–Meier Plotter, we used the best cutoff value of COL4A3, MEF2C and TRIM15 expression to divide diffuse-type gastric cancer patients into a high expression- and low expression-group, respectively. It was found that the high expression of COL4A3 and MEF2C and the low expression of TRIM15 were associated with worse OS for diffuse-type gastric cancer patients (Figure 3). However, the false discovery rate (FDR) was 50%, 50% and over 50% for the survival difference of COL4A3, MEF2C and TRIM15, respectively. These FDR values were high. Thus, it made it difficult for us to evaluate the prognostic value of these genes.
A univariate Cox regression analysis was performed to evaluate the prognostic value and identify the risk factors in diffuse-type gastric cancer. The results of the univariate Cox regression analysis demonstrated that T stage (p < 0.01), N stage (p < 0.01), M stage (p < 0.001) and the expression level of MEF2C (p < 0.01) and TRIM15 (p < 0.05) were significantly correlated with overall survival in diffuse-type gastric cancer (Table 3).

3.5. Establishment of the Prognostic Nomogram of Diffuse-Type Gastric Cancer

The clinical information of the 134 diffuse-type gastric cancer patients from the GSE62254 dataset was used to construct a prognostic nomogram for predicting 1-, 3-, 5-year overall survival based on a stepwise Cox regression model (Figure 4A). T stage, N stage, M stage and the expression of TRIM15 and MEF2C were parameters included in the nomogram. The calibration curves showed a good consistency between the actual and the nomogram-predicted 1-, 3- and 5-year overall survival probabilities (Figure 4B). The risk score was calculated as Formula (1):
e^((−0.52425 + [0.621×Expression value of MEF2C] + [(−1.2188)×Expression value of TRIM15] + βT + βN + βM))
When T stage is T2, T3 or T4, the value of βT is 0, 0.3812 or −0.2084, respectively.
When N stage is N0, N1, N2 or N3, the value of βN is 0, 0.8741, 1.4293 or 2.7244, respectively.
When M stage is M0 or M1, the value of βM is 0 or 0.994, respectively.
Figure 4. Construction and internal validation of nomogram. (A) The predictive nomogram was established with T-stage, N-stage, M-stage and expression level of MEF2C and TRIM15 in diffuse-type gastric cancer. (B) A comparison of predictive and actual outcome for 1-, 3- and 5-year survival probabilities in the nomogram is demonstrated in calibration curves for internal validation.
Figure 4. Construction and internal validation of nomogram. (A) The predictive nomogram was established with T-stage, N-stage, M-stage and expression level of MEF2C and TRIM15 in diffuse-type gastric cancer. (B) A comparison of predictive and actual outcome for 1-, 3- and 5-year survival probabilities in the nomogram is demonstrated in calibration curves for internal validation.
Curroncol 30 00014 g004

3.6. Evaluation the Predictive Performance of Nomogram and External Validation of Nomogram

The AUC of the predicted 1-, 3- and 5-year overall survival were 0.82, 0.84 and 0.87, respectively (Figure 5A). When the seventh AJCC TNM stage was used, the AUC values for the 1-, 3- and 5-year overall survival predictions were 0.79, 0.79 and 0.84, respectively (Figure 5A). As seen in Figure 4A, the calculated overall score could estimate the survival prognosis (1-, 3- and 5-year survival probabilities), and the C-index of our nomogram model was 0.766 (95% CI = 0.711 ~ 0.821). The Kaplan–Meier analysis showed a significant difference in prognostic value between the high-risk and low-risk groups (Figure 5B). Decision curve analysis showed the effectiveness of the nomogram was better than the seventh AJCC TNM staging system.
For external validation, using the TCGA-STAD dataset, the results of ROC curves, Kaplan–Meier analysis and decision curve analysis were similar with the training set (GSE62254). The AUC of predicted 1-, 3- and 5-year overall survival were 0.61, 0.69 and 0.77, which were better than the TNM stating system (Figure 5D) in the TCGA-STAD dataset. Also, the patients in the low-risk group had a more favorable overall prognosis than the high-risk group in the validation set (Figure 5E). Similarly, decision curve analysis also demonstrated the net benefit of the nomogram was better than the TNM staging system in the TCGA-STAD dataset (Figure 5F). Therefore, the nomogram showed better discriminatory ability than the seventh AJCC TNM classification.

3.7. Scatter Point

The patients were divided into two groups according to the scoring of the nomogram (Figure 6A). We utilized a scatter plot to reveal the relationship between the level of risk score and the overall survival of diffuse-type GC. As a result, the high-risk group exhibited significantly poorer overall survival (Figure 6B).

3.8. Gene Set Enrichment Analysis (GSEA)

To further investigate the molecular mechanism of MEF2C and TRIM15 in diffuse-type gastric cancer, GSEA was conducted. The results showed high expression MEF2C enriched in the MAPK signaling pathway (Figure 7A) and low expression TRIM15 enriched in glycosaminoglycan biosynthesis chondroitin sulfate (Figure 7B). In addition, the high-risk group was compared to the low-risk group in GSEA. We demonstrated that the high-risk group was significantly enriched in the top 13 biological pathways, including “Regulation of actin cytoskeleton”, “Focal adhesion”, “Vasopressin regulated water reabsorption”, “Dilated cardiomyopathy”, “Glycosaminoglycan biosynthesis chondroitin sulfate”, “Calcium signaling pathway”, “Notch signaling pathway”, “ECM receptor interaction”, “Phosphatidylinositol signaling system”, “MAPK signaling pathway”, ”Hypertrophic cardiomyopathy”, “Vascular smooth muscle contraction”, “FcγR-Mediated Phagocytosis” (Supplementary Materials Figure S1).

4. Discussion

Compared to intestinal-type gastric cancer, diffuse-type gastric cancer exhibits a more aggressive phenotype with a relatively poor prognosis and a 5-year overall survival rate of 32.1% [15]. It is widely acknowledged that the treatment failure of diffuse-type gastric cancer is due to drug resistance and disease progression, including tumor recurrence and metastasis. The prognostic model is important to clinicians to provide individualized treatment by determining which patients would benefit most from a particular or a combination of treatment approaches, including radical surgery, adjuvant chemotherapy, neoadjuvant chemotherapy, targeted molecular medicine or immunotherapy. However, today the big problem is that the prognostic model based on clinical characteristics and histopathological characteristics is not accurate [16,17]. Accordingly, it is important to develop a novel prognostic model for improving patient management by stratifying patients according to their characteristics.
In the present study, we constructed a nomogram that incorporated a two-gene (MEF2C and TRIM15) signature and clinicopathological parameters to assist clinicians in determining the prognosis of individual diffuse-type GC patients. The sensitivity and specificity of our prognostic model were more satisfying than the TNM staging system (Figure 5). Gene set enrichment analysis showed that MEF2C and TRIM15 were closely related to invasive and metastasis signaling pathways in diffuse-type gastric cancer. The MAPK signaling pathway was the most significant in the high MEF2C expression group of diffuse-type gastric cancer (Figure 7A). Moreover, the glycosaminoglycan biosynthesis chondroitin sulfate signaling pathway was the most enriched in the low TRIM15 expression group of diffuse-type gastric cancer (Figure 7B). MEF2C upregulation and TRIM15 downregulation in diffuse-type gastric cancer were related to poor prognoses.
In recent years, multiple gene signatures, mRNAs or non-coding RNAs have been used to evaluate the prognosis of gastric cancer patients [18,19,20]. Nevertheless, rare studies have focused on the Lauren subtype-specific gene signature to evaluate the prognosis of diffuse-type gastric cancer. Moreover, few studies have sought to combine the TNM stage with the multi-gene signature to assess the prognosis of diffuse-type gastric cancer. A previous study reported a three-gene signature to predict the prognosis of diffuse-type gastric cancer [21]. However, the prognostic model only considered the three-gene expression level, but lacked the clinical parameters of the diffuse-type gastric cancer patients. The TNM staging system only considers tumor invasion depth, lymph node metastasis and distant metastasis. The biological characteristics of the tumors, such as the immune infiltration status, drug response and intracellular signal pathways are not reflected in the TNM staging system. However, the genomic sequence of the tumor is an effective tool to uncover heterogeneous malignance [22,23]. Several tumor biomarkers can help guide treatment decisions, including Human Epidermal Growth Factor Receptor-2 (HER2), Programmed Cell Death-Ligand 1 (PDL1) and Vascular Endothelial Growth Factor Receptor (VEGFR) [24,25,26]. Accordingly, in the current study, we identified risk factors, including age, T stage, N stage, M stage, the expression level of TRIM15 and MEF2C and established a prognostic model. Finally, a nomogram integrating a two-gene signature and clinicopathologic features was constructed and yielded an accurate prediction of overall survival. Through ROC curves, Kaplan–Meier analysis and decision curve analysis of the external validation in the TCGA-STAD dataset, as a supplement to AJCC staging, our two-gene signature and nomogram demonstrated a similar predictive performance with the training set (Figure 5). Our predictive nomogram will exhibit a potential value of diffuse-type GC in future clinical practice. Similarly, several previous studies integrated clinical features and risk scores based on the expression level of risk genes into a novel prognostic nomogram [27,28,29]. The predictive value of their integrated nomograms was also better than using the risk factor alone. These studies and our present study have a certain reference significance for future clinical research.
MEF2C and TRIM15 have previously been reported to be associated with gastric cancer. Interestingly, Myocyte Enhancer Factor 2C (MEF2C) has been documented in pathways of organelle biogenesis and maintenance and transcriptional misregulation in cancer, which involved DNA-binding transcription factor activity and protein heterodimerization activity. MEF2C has been associated with DNA methylation and enhanced PD-L1 expression in gastric cancer [30,31]. Recent studies have also shown that MEF2C plays an important role in myocilin mediating cancer-induced muscle wasting and cachexia in cancer patients [32] and regulates chemotherapeutic resistance [33] and the disease progression of acute myeloid leukemia [34]. TRIM15 is a member of the tripartite motif (TRIM) family. The protein encoded by TRIM15 has a TRIM motif, including three zinc-binding domains, a RING, a B-box type 1, a B-box type 2 and a coiled-coil region. However, the biological function of TRIM15 remains unknown. Our GSEA results showed that TRIM15 was correlated with the glycosaminoglycan biosynthesis chondroitin sulfate signaling pathway in diffuse-type gastric cancer. Importantly, a recent study has found that the expression of TRIM15 is an independent risk factor of prognosis in gastric cancer patients [35]. However, the roles of the MEF2C and TRIM15 genes in diffuse-type gastric cancer are still unclear. Our current study disclosed that MEF2C and TRIM15 could promote invasion and metastasis through cancer-related signaling pathways: the MEF2C-activated MAPK signaling pathway and the TRIM15-activated glycosaminoglycan biosynthesis chondroitin sulfate signaling pathway (Figure 7). Moreover, our studies revealed poor prognoses associated with upregulated MEF2C and downregulated TRIM15 expression in diffuse-type gastric cancer. Accordingly, the present research revealed the roles of these two genes in diffuse-type gastric cancer and established a risk model to complement the AJCC staging system to improve the outcomes of diffuse-type gastric cancer. However, further in vivo studies are needed to explore the molecular mechanism underlying the oncological function of these two genes in diffuse-type gastric cancer.
This study contains several limitations. First, our study was based on RNA sequence data rather than proteomics, which could have affected the accuracy of our prediction model. Accordingly, the expression of these two genes should be analyzed in another study with a large sample of diffuse-type gastric cancer patients to validate the predictive performance of our model. Furthermore, it may be hard to promote the utilization of multi-genome sequencing during clinical practice due to its high price and practicability. With the development of sequencing technology and precision medicine, the identified two-gene signature will be clinically feasible. Moreover, our predictive model should be externally validated with another large sample of diffuse-type gastric cancer patients.

5. Conclusions

In summary, we constructed a nomogram that incorporated a two-gene signature and clinicopathological parameters to assist clinicians in determining the prognosis of individual GC patients. Our nomogram is simple to use and can be harnessed to provide optimal treatment and make medical decisions. To the best of our knowledge, the two-gene prognostic signature described and the nomogram constructed have not been reported previously. The current study provides a new perspective of the molecular mechanisms underlying prognosis prediction in diffuse-type gastric cancer. In addition, MEF2C and TRIM15 were obtained by a pooled analysis of multiple datasets and are accordingly highly reliable. Importantly, these two genes may be potential molecular targets for the treatment of diffuse-type gastric cancer.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/curroncol30010014/s1, Figure S1: GSEA showing the top 13 enriched biological pathways of high-risk group.

Author Contributions

Conceptualization, S.C. and J.X.; methodology, S.Y.; software, H.W. (Huabin Wang); validation, G.L., X.J. and J.L.; formal analysis, J.X.; data curation, J.Z., H.W. (Huijin Wang), H.L. and H.W. (Han Wang); writing—original draft preparation, S.C.; writing—review and editing, J.X. and S.Y.; supervision, C.Z. and Y.H. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Sanming Project of Medicine in Shenzhen (No. SZSM201911010) and the Guangdong Province Key Laboratory of Digestive Cancer Research (No. 2021B1212040006).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All dataset analysis for this study can be found in the Gene Expression Omnibus (GEO) database, the Cancer Genome Atlas (TCGA) database and the website of Kaplan–Meier Plotter.

Acknowledgments

The authors would like to acknowledge the Gene Expression Omnibus (GEO) database and the Cancer Genome Atlas (TCGA) database.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Smyth, E.C.; Nilsson, M.; Grabsch, H.I.; van Grieken, N.C.; Lordick, F. Gastric cancer. Lancet 2020, 396, 635–648. [Google Scholar] [CrossRef] [PubMed]
  2. Van Cutsem, E.; Sagaert, X.; Topal, B.; Haustermans, K.; Prenen, H. Gastric cancer. Lancet 2016, 388, 2654–2664. [Google Scholar] [CrossRef] [PubMed]
  3. Lauren, P. The Two Histological Main Types of Gastric Carcinoma: Diffuse and So-Called Intestinal-Type Carcinoma. An Attempt at a Histo-Clinical Classification. Acta Pathol. Microbiol. Scand. 1965, 64, 31–49. [Google Scholar] [CrossRef] [PubMed]
  4. Aoyama, T.; Hutchins, G.; Arai, T.; Sakamaki, K.; Miyagi, Y.; Tsuburaya, A.; Ogata, T.; Oshima, T.; Earle, S.; Yoshikawa, T.; et al. Identification of a high-risk subtype of intestinal-type Japanese gastric cancer by quantitative measurement of the luminal tumor proportion. Cancer Med. 2018, 7, 4914–4923. [Google Scholar] [CrossRef] [PubMed]
  5. Chen, Y.C.; Fang, W.L.; Wang, R.F.; Liu, C.A.; Yang, M.H.; Lo, S.S.; Wu, C.W.; Li, A.F.; Shyr, Y.M.; Huang, K.H. Clinicopathological Variation of Lauren Classification in Gastric Cancer. Pathol. Oncol. Res. 2016, 22, 197–202. [Google Scholar] [CrossRef]
  6. Lauren, P. Histogenesis of intestinal and diffuse types of gastric carcinoma. Scand. J. Gastroenterol. Suppl. 1991, 180, 160–164. [Google Scholar]
  7. Qiu, M.Z.; Cai, M.Y.; Zhang, D.S.; Wang, Z.Q.; Wang, D.S.; Li, Y.H.; Xu, R.H. Clinicopathological characteristics and prognostic analysis of Lauren classification in gastric adenocarcinoma in China. J. Transl. Med. 2013, 11, 58. [Google Scholar] [CrossRef] [Green Version]
  8. Lee, J.Y.; Gong, E.J.; Chung, E.J.; Park, H.W.; Bae, S.E.; Kim, E.H.; Kim, J.; Do, Y.S.; Kim, T.H.; Chang, H.S.; et al. The Characteristics and Prognosis of Diffuse-Type Early Gastric Cancer Diagnosed during Health Check-Ups. Gut Liver 2017, 11, 807–812. [Google Scholar] [CrossRef]
  9. Fang, C.; Wang, W.; Deng, J.Y.; Sun, Z.; Seeruttun, S.R.; Wang, Z.N.; Xu, H.M.; Liang, H.; Zhou, Z.W. Proposal and validation of a modified staging system to improve the prognosis predictive performance of the 8th AJCC/UICC pTNM staging system for gastric adenocarcinoma: A multicenter study with external validation. Cancer Commun. (Lond) 2018, 38, 67. [Google Scholar] [CrossRef] [Green Version]
  10. Son, T.; Sun, J.; Choi, S.; Cho, M.; Kwon, I.G.; Kim, H.I.; Cheong, J.H.; Choi, S.H.; Noh, S.H.; Woo, Y.; et al. Multi-institutional validation of the 8th AJCC TNM staging system for gastric cancer: Analysis of survival data from high-volume Eastern centers and the SEER database. J. Surg. Oncol. 2019, 120, 676–684. [Google Scholar] [CrossRef]
  11. Wang, Y.; Zhang, J.; Guo, S.; Dong, Z.; Meng, X.; Zheng, G.; Yang, D.; Zheng, Z.; Zhao, Y. Implication of lymph node staging in migration and different treatment strategies for stage T2N0M0 and T1N1M0 resected gastric cancer: A SEER population analysis. Clin. Transl. Oncol. 2019, 21, 1499–1509. [Google Scholar] [CrossRef] [PubMed]
  12. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015, 43, e47. [Google Scholar] [CrossRef] [PubMed]
  13. Lanczky, A.; Gyorffy, B. Web-Based Survival Analysis Tool Tailored for Medical Research (KMplot): Development and Implementation. J. Med. Internet Res. 2021, 23, e27633. [Google Scholar] [CrossRef]
  14. Szasz, A.M.; Lanczky, A.; Nagy, A.; Forster, S.; Hark, K.; Green, J.E.; Boussioutas, A.; Busuttil, R.; Szabo, A.; Gyorffy, B. Cross-validation of survival associated biomarkers in gastric cancer using transcriptomic data of 1,065 patients. Oncotarget 2016, 7, 49322–49333. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Kao, Y.C.; Fang, W.L.; Wang, R.F.; Li, A.F.; Yang, M.H.; Wu, C.W.; Shyr, Y.M.; Huang, K.H. Clinicopathological differences in signet ring cell adenocarcinoma between early and advanced gastric cancer. Gastric Cancer 2019, 22, 255–263. [Google Scholar] [CrossRef] [Green Version]
  16. Kwa, M.; Makris, A.; Esteva, F.J. Clinical utility of gene-expression signatures in early stage breast cancer. Nat. Rev. Clin. Oncol. 2017, 14, 595–610. [Google Scholar] [CrossRef]
  17. Wang, K.; Li, H.L.; Xiong, Y.F.; Shi, Y.; Li, Z.Y.; Li, J.; Zhang, X.; Li, H.Y. Development and validation of nomograms integrating immune-related genomic signatures with clinicopathologic features to improve prognosis and predictive value of triple-negative breast cancer: A gene expression-based retrospective study. Cancer Med. 2019, 8, 686–700. [Google Scholar] [CrossRef]
  18. Ghafouri-Fard, S.; Taheri, M. Long non-coding RNA signature in gastric cancer. Exp. Mol. Pathol. 2020, 113, 104365. [Google Scholar] [CrossRef]
  19. Frilling, A.; Modlin, I.M.; Kidd, M.; Russell, C.; Breitenstein, S.; Salem, R.; Kwekkeboom, D.; Lau, W.Y.; Klersy, C.; Vilgrain, V.; et al. Recommendations for management of patients with neuroendocrine liver metastases. Lancet Oncol. 2014, 15, e8–e21. [Google Scholar] [CrossRef]
  20. Li, X.; Zhang, Y.; Zhang, Y.; Ding, J.; Wu, K.; Fan, D. Survival prediction of gastric cancer by a seven-microRNA signature. Gut 2010, 59, 579–585. [Google Scholar] [CrossRef]
  21. Bao, B.; Zheng, C.; Yang, B.; Jin, Y.; Hou, K.; Li, Z.; Zheng, X.; Yu, S.; Zhang, X.; Fan, Y.; et al. Identification of Subtype-Specific Three-Gene Signature for Prognostic Prediction in Diffuse Type Gastric Cancer. Front. Oncol. 2019, 9, 1243. [Google Scholar] [CrossRef] [PubMed]
  22. Hong, M.; Tao, S.; Zhang, L.; Diao, L.T.; Huang, X.; Huang, S.; Xie, S.J.; Xiao, Z.D.; Zhang, H. RNA sequencing: New technologies and applications in cancer research. J. Hematol. Oncol. 2020, 13, 166. [Google Scholar] [CrossRef]
  23. Dentro, S.C.; Leshchiner, I.; Haase, K.; Tarabichi, M.; Wintersinger, J.; Deshwar, A.G.; Yu, K.; Rubanova, Y.; Macintyre, G.; Demeulemeester, J.; et al. Characterizing genetic intra-tumor heterogeneity across 2,658 human cancer genomes. Cell 2021, 184, 2239–2254 e2239. [Google Scholar] [CrossRef] [PubMed]
  24. Meric-Bernstam, F.; Johnson, A.M.; Dumbrava, E.E.I.; Raghav, K.; Balaji, K.; Bhatt, M.; Murthy, R.K.; Rodon, J.; Piha-Paul, S.A. Advances in HER2-Targeted Therapy: Novel Agents and Opportunities Beyond Breast and Gastric Cancer. Clin. Cancer Res. 2019, 25, 2033–2041. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Brody, R.; Zhang, Y.; Ballas, M.; Siddiqui, M.K.; Gupta, P.; Barker, C.; Midha, A.; Walker, J. PD-L1 expression in advanced NSCLC: Insights into risk stratification and treatment selection from a systematic literature review. Lung Cancer 2017, 112, 200–215. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Malekan, M.; Ebrahimzadeh, M.A. Vascular Endothelial Growth Factor Receptors [VEGFR] as Target in Breast Cancer Treatment: Current Status in Preclinical and Clinical Studies and Future Directions. Curr. Top. Med. Chem. 2022, 22, 891–920. [Google Scholar] [CrossRef]
  27. Qing, X.; Xu, W.; Liu, S.; Chen, Z.; Ye, C.; Zhang, Y. Molecular Characteristics, Clinical Significance, and Cancer Immune Interactions of Angiogenesis-Associated Genes in Gastric Cancer. Front. Immunol. 2022, 13, 843077. [Google Scholar] [CrossRef]
  28. Bai, Y.; Wei, C.; Zhong, Y.; Zhang, Y.; Long, J.; Huang, S.; Xie, F.; Tian, Y.; Wang, X.; Zhao, H. Development and Validation of a Prognostic Nomogram for Gastric Cancer Based on DNA Methylation-Driven Differentially Expressed Genes. Int. J. Biol. Sci. 2020, 16, 1153–1165. [Google Scholar] [CrossRef]
  29. Nie, K.; Zheng, Z.; Wen, Y.; Shi, L.; Xu, S.; Wang, X.; Zhou, Y.; Fu, B.; Li, X.; Deng, Z.; et al. Construction and validation of a TP53-associated immune prognostic model for gastric cancer. Genomics 2020, 112, 4788–4795. [Google Scholar] [CrossRef]
  30. Liu, D.; Ma, X.; Yang, F.; Xiao, D.; Jia, Y.; Wang, Y. Discovery and validation of methylated-differentially expressed genes in Helicobacter pylori-induced gastric cancer. Cancer Gene Ther. 2020, 27, 473–485. [Google Scholar] [CrossRef]
  31. Menyhart, O.; Pongor, L.S.; Gyorffy, B. Mutations Defining Patient Cohorts With Elevated PD-L1 Expression in Gastric Cancer. Front. Pharmacol. 2018, 9, 1522. [Google Scholar] [CrossRef] [PubMed]
  32. Judge, S.M.; Deyhle, M.R.; Neyroud, D.; Nosacka, R.L.; D’Lugos, A.C.; Cameron, M.E.; Vohra, R.S.; Smuder, A.J.; Roberts, B.M.; Callaway, C.S.; et al. MEF2c-Dependent Downregulation of Myocilin Mediates Cancer-Induced Muscle Wasting and Associates with Cachexia in Patients with Cancer. Cancer Res. 2020, 80, 1861–1874. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Brown, F.C.; Still, E.; Koche, R.P.; Yim, C.Y.; Takao, S.; Cifani, P.; Reed, C.; Gunasekera, S.; Ficarro, S.B.; Romanienko, P.; et al. MEF2C Phosphorylation Is Required for Chemotherapy Resistance in Acute Myeloid Leukemia. Cancer Discov. 2018, 8, 478–497. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Tarumoto, Y.; Lu, B.; Somerville, T.D.D.; Huang, Y.H.; Milazzo, J.P.; Wu, X.S.; Klingbeil, O.; El Demerdash, O.; Shi, J.; Vakoc, C.R. LKB1, Salt-Inducible Kinases, and MEF2C Are Linked Dependencies in Acute Myeloid Leukemia. Mol. Cell 2018, 69, 1017–1027 e1016. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Zhou, W.; Chen, H.; Ruan, Y.; Zeng, X.; Liu, F. High Expression of TRIM15 Is Associated with Tumor Invasion and Predicts Poor Prognosis in Patients with Gastric Cancer. J. Investig. Surg. 2020, 34, 853–861. [Google Scholar] [CrossRef]
Figure 1. Flowchart showing the process of establishing the gene signature and prognostic nomogram of diffuse-type gastric cancer in this study.
Figure 1. Flowchart showing the process of establishing the gene signature and prognostic nomogram of diffuse-type gastric cancer in this study.
Curroncol 30 00014 g001
Figure 2. (A) The Venn diagram shows the DEGs with upregulated genes and downregulated genes in diffuse-type GC. (B) Expression levels of four DEGs between diffuse-type and intestinal-type gastric cancer were validated in GSE62254 dataset.
Figure 2. (A) The Venn diagram shows the DEGs with upregulated genes and downregulated genes in diffuse-type GC. (B) Expression levels of four DEGs between diffuse-type and intestinal-type gastric cancer were validated in GSE62254 dataset.
Curroncol 30 00014 g002
Figure 3. Validation of prognostic value of three genes in diffuse-type gastric cancer (n = 241) by Kaplan–Meier Plotter.
Figure 3. Validation of prognostic value of three genes in diffuse-type gastric cancer (n = 241) by Kaplan–Meier Plotter.
Curroncol 30 00014 g003
Figure 5. (A) Receiver operating characteristic (ROC) curves were used to compare the predictive ability of nomogram model and TNM stage for 1-, 3- and 5-year survival probabilities in training set. (B) The Kaplan–Meier curve shows the high-risk group with worse prognosis in training set. (C) Decision curve analysis of nomogram in training set for OS. (D) The ROC curves of nomogram in external validation set (TCGA-STAD). (E) Kaplan–Meier analysis of high-risk group and low-risk group in external validation set (TCGA-STAD). (F) Decision curve analysis of nomogram in external validation set.
Figure 5. (A) Receiver operating characteristic (ROC) curves were used to compare the predictive ability of nomogram model and TNM stage for 1-, 3- and 5-year survival probabilities in training set. (B) The Kaplan–Meier curve shows the high-risk group with worse prognosis in training set. (C) Decision curve analysis of nomogram in training set for OS. (D) The ROC curves of nomogram in external validation set (TCGA-STAD). (E) Kaplan–Meier analysis of high-risk group and low-risk group in external validation set (TCGA-STAD). (F) Decision curve analysis of nomogram in external validation set.
Curroncol 30 00014 g005
Figure 6. The association between risk score, OS and survival status in the predictive nomogram is shown in scatter plot. (A) Risk score ranking of patients. (B) Distribution of risk score and associated overall survival.
Figure 6. The association between risk score, OS and survival status in the predictive nomogram is shown in scatter plot. (A) Risk score ranking of patients. (B) Distribution of risk score and associated overall survival.
Curroncol 30 00014 g006
Figure 7. Gene set enrichment analysis was used to showed MEF2C high expression group (A) and TRIM15 low expression group (B) enriched signaling pathways.
Figure 7. Gene set enrichment analysis was used to showed MEF2C high expression group (A) and TRIM15 low expression group (B) enriched signaling pathways.
Curroncol 30 00014 g007
Table 1. Details of GEO datasets in this study.
Table 1. Details of GEO datasets in this study.
DatasetsPlatformCase No.Sample SizeSubmitterApplication
GSE22377(G-U133_Plus_2)
Affymetrix Human
Genome U133 Plus 2.0 Array
4324 intestinal
19 diffuse
Förster S, MDC, Berlin, GermanyIdentification
of DEGs
GSE38749(HG-U133_Plus_2)
Affymetrix Human
Genome U133 Plus 2.0 Array
154 intestinal
10 diffuse
1 mixed
Pasini FS, Faculdade
de Medicina da USP, São Paulo, Brazil
Identification
of DEGs
GSE47007(HG_U95Av2)
Affymetrix Human
Genome U95 Version 2 Array
3018 intestinal
12 diffuse
Sasaki H, National
Cancer Center Reseach Institute, Tokyo, Japan
Identification
of DEGs
GSE62254(HG-U133_Plus_2)
Affymetrix Human
Genome U133 Plus 2.0 Array
300146 intestinal
134 diffuse
17 mixed
3 indeterminate
Nebozhyn M,
Merck, Inc., PA, USA
Validation
Table 2. Stratified analysis of COL4A3, MEF2C and TRIM15 for diffuse-type gastric cancer patients (n = 134) in terms of prognosis.
Table 2. Stratified analysis of COL4A3, MEF2C and TRIM15 for diffuse-type gastric cancer patients (n = 134) in terms of prognosis.
Clinical
Features
COL4A3 Expressionp ValueMEF2C Expressionp ValueTRIM15 Expressionp Value
HighLow HighLow HighLow
Age≤6032360.15145230.016 *31370.014 *
>602343 3036 4422
GenderMale28320.23438220.12231290.366
Female2747 3737 4430
Depth of invasionT218460.014 *25390.001 **4024
T33228 4317 3129
T455 73 46
Lymph nodeN0440.178530.33440.678
metastasisN11637 2528 3320
N21823 2714 2219
N31715 1814 1616
Distant metastasisM045690.37760540.06364500.925
M11010 155 119
AJCC stageI140.59230.008 **410.51
II826 1123 2113
III2326 3316 2722
IV2323 2917 2323
* p value < 0.05; ** p value < 0.01.
Table 3. The univariate Cox regression analysis between three genes and other clinical characteristics and OS in diffuse-type gastric cancer.
Table 3. The univariate Cox regression analysis between three genes and other clinical characteristics and OS in diffuse-type gastric cancer.
CharacteristicsHR95% CIp Value
Age1.0050.987–1.0230.617
GenderFemale1
Male0.77120.495–1.2030.252
Depth of InvasionT21
T31.9061.181–3.0780.008 **
T42.7691.301–5.8930.008 **
Lymph Node of MetastasisN01
N11.4880.346–6.3980.593
N23.1770.753–13.4110.115
N311.0172.608–46.5390.001 **
Distant MetastasisM01
M13.6772.172–6.224<0.0001 ***
Gene ExpressionCOL4A3 expression1.2810.685–2.3950.438
MEF2C expression3.971.512–10.4300.005 **
TRIM15 expression0.28610.104–0.7900.015 *
* p value < 0.05; ** p value < 0.01; *** p value < 0.001; HR, hazard ratio; CI, confidence interval.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, S.; Xu, J.; Yin, S.; Wang, H.; Liu, G.; Jin, X.; Zhang, J.; Wang, H.; Wang, H.; Li, H.; et al. Identification of a Two-Gene Signature and Establishment of a Prognostic Nomogram Predicting Overall Survival in Diffuse-Type Gastric Cancer. Curr. Oncol. 2023, 30, 171-183. https://doi.org/10.3390/curroncol30010014

AMA Style

Chen S, Xu J, Yin S, Wang H, Liu G, Jin X, Zhang J, Wang H, Wang H, Li H, et al. Identification of a Two-Gene Signature and Establishment of a Prognostic Nomogram Predicting Overall Survival in Diffuse-Type Gastric Cancer. Current Oncology. 2023; 30(1):171-183. https://doi.org/10.3390/curroncol30010014

Chicago/Turabian Style

Chen, Songyao, Jiannan Xu, Songcheng Yin, Huabin Wang, Guangyao Liu, Xinghan Jin, Junchang Zhang, Huijin Wang, Han Wang, Huan Li, and et al. 2023. "Identification of a Two-Gene Signature and Establishment of a Prognostic Nomogram Predicting Overall Survival in Diffuse-Type Gastric Cancer" Current Oncology 30, no. 1: 171-183. https://doi.org/10.3390/curroncol30010014

APA Style

Chen, S., Xu, J., Yin, S., Wang, H., Liu, G., Jin, X., Zhang, J., Wang, H., Wang, H., Li, H., Liang, J., He, Y., & Zhang, C. (2023). Identification of a Two-Gene Signature and Establishment of a Prognostic Nomogram Predicting Overall Survival in Diffuse-Type Gastric Cancer. Current Oncology, 30(1), 171-183. https://doi.org/10.3390/curroncol30010014

Article Metrics

Back to TopTop