- Research article
- Open Access
Expression and prognostic value of cell-cycle-associated genes in gastric adenocarcinoma
BMC Gastroenterology volume 18, Article number: 81 (2018)
Gastric carcinoma is a malignant disease, and gastric adenocarcinoma (GAC) is the most common histological type. Molecular profiling of GAC has been extensively performed, but few have focused on the clinical significance of gene clusters of the cell cycle.
We investigated the genetic profile of cell-cycle-associated genes in a GAC cohort. The mRNA expression and clinical data were downloaded from TCGA, according to cBioportal. We conducted a series of analyses to detect the relationships between these genes and GAC.
From all the patients, 5 clusters were identified based on mRNA expression of 122 cell-cycle-associated genes. Cluster 1 showed the worst prognosis and is characterized by extremely high expression of WEE2 and CCNE1. Comparison of the gene patterns showed that 16 genes expressed were distinctly varied between each cluster. In addition, investigations into the prognostic role of the 16 genes suggested that high expression of ESPL1 and MCM5 were significantly correlated with favorable outcomes. Moreover, we detected that ESPL1 and MCM5 gene expression were negatively correlated with GAC pathologic stage progression.
This study revealed a gene expression pattern of the cell cycle in different GAC subgroups, and suggested individual genes were associated with the clinical outcome and AJCC stages. These results suggest a novel prognostic strategy for GAC and provide information for patient stratification and trials of targeted therapies.
Gastric carcinoma (GC) remains the fifth most prevalent cancer and the third leading cause of cancer-related deaths worldwide [1, 2]. Gastric adenocarcinoma (GAC) is the main type of GC and is associated with poor survival rates . Although the incidence of GAC has been reduced over the past years, it imposes a critical issue globally . The characteristics of diverse histological (phenotypes) and genotypes manifest that GAC is a heterogeneous disease . Previous studies uncovered the genetic profiling of GAC by performing gene expression or DNA sequencing, but rare profiling led to alterations of gene panels associated with specific biological events [5,6,7]. The cell cycle process is a highly organized event and regulated duplication of genetic material and cell division, and aberrant cell cycle activity is a hallmark of cancer . A comprehensive molecular characterization study elicited frequent amplifications of cell cycle mediators (CCNE1, CCND1 and CDK6) in GAC; however, it remains unclear if specific patterns of cell-cycle-associated gene expression across different subtypes or stages bears any significance on patient outcome or is correlated with genetic alterations .
Tumor stages are generally used to characterize disease progression and determine metastasis and prognosis . Gene expression involved in the development of GC is also variant in different stages . In this study, we focused on mRNA expressing variations of genes associated with the cell cycle in different GAC subgroups, as well as its correlation to tumor stages. We identified significantly differentially expressed genes in GAC samples and the prognostic impact . Specific genes that impacted the outcomes showed a variant expression panel between different stages. Finally, we report a unique set of cell-cycle-associated genes in GAC that serve as a divider of biological characteristics, which revealed the potential for therapeutic strategy targeting cell-cycle-associated genes.
Samples and database
We obtained RNA-Seq data and the corresponding clinical records of 415 gastric adenocarcinoma patients (GAC) of TCGA from cBioPortal for Cancer Genomics (http://cbioportal.org) . We filtrated the data base on whether the mRNA z-score, tumor stage and overall survival were clearly recorded. Collectively, the data set included 228 samples for the clustering study and 206 samples for the stage study, respectively.
A panel of cell-cycle-associated genes was derived from the KEGG pathway database (http://www.kegg.jp/kegg/), as previously described . In total, 124 genes were listed; however, mRNA expression data of genes MYC and CCND1 were unobtainable, and the remaining 122 candidate genes were analyzed (Additional file 1: Table S1).
Hierarchical cluster analysis of 122 genes expressed in each sample was used to group samples with similar gene expression patterns. Samples with differential cell-cycle-associated gene expression between separate clusters were identified from the whole data set. The expression levels of GAC patients were shown as mRNA z-scores and grouped using the hierarchical clustering algorithm in the Gene Cluster 3.0 program . The heat-map of cluster indication and tumor stage divided pattern was generated using the Java Treeview program .
Prognostic implication analyses
To evaluate the relationship between the cell-cycle-associated genes and patient clinical outcome, we used GraphPad Prism 6 for Windows (GraphPad Software, Inc., California, US; Version 6.01, 2012) for construction and overall survival (OS) comparisons in different gastric adenocarcinoma clusters and different tumor stages. Additionally, OS difference analysis between low and high expression cohorts were conducted in GraphPad Prism 6.
Survival curves were plotted using a Kapla-Meier analysis and compared using the log-rank test in GraphPad Prism 6. Associations between clinical characteristics and clustering variable were analyzed by Fisher’s exact test and Pearson/Spearman correlation. Difference of gene expression between clusters was conducted by ANOVA analysis. Correlation between each element was conducted by regression analysis. Both results were performed in SPSS 19.0 (IBM, Inc., New York, US). A P value of less than 0.05 was considered statistically significant.
Expression profile of cell-cycle-associated genes in GAC
To understand biological functions and utilities of the biological system, we retrieved 124 cell-cycle-associated genes (PATHWAY: map04110) from the KEGG database. Additionally, we selected 122 genes that were detectable and sorted the mRNA expression values from the RNA-seq data of GAC (Additional file 2: Figure S1).
The expression levels of the 122 genes were calculated by mRNA z-scores compared to the expression distribution of each gene from tumors that were diploid for the genes in 415 GAC cases (RNA Seq V2 RSEM), based on TCGA data in GC. Following the clinical filtrated criteria, we achieved 228 cases with survival records to perform cluster analysis. Five clusters of GAC were distinguished based on the gene expression panels (Fig. 1a). The statistical analysis of the clinical and histological characteristics of each clusters was performed, and the result shows that the nodes pathologic stages and tumor stages was significantly different (The characteristic of histological types was excluded for the reason of the incompleteness of data) (Table 1). Moreover, correlational analysis showed the significant correlations of the three parameters and clustering situation (Table 1).
We compared the median survival and survival curves between each cluster further and the results showed that cluster 1 demonstrated the worst prognosis (12.35 months) compared to other clusters (26.31 months) (Fig. 1b). In contrast, clusters 2 and 4 showed favorable prognosis with a median survival of 30.88 months and 26.45 months, respectively (Table 2). Comparison of survival curves revealed a significant difference between outcomes of cluster 1 and cluster 4 (P = 0.0278) (Additional file 1: Table S2).
Variation of cell-cycle-associated gene expression in different clusters
Additionally, we examined the genes that expressed variation between each cluster. According to the filter criteria of significant difference (P < 0.01), there were 42 differentially expressed genes of cluster 1 versus the other clusters, 32 differentially expressed genes of cluster 1 versus cluster 2 and 38 differentially expressed genes of cluster 1 versus cluster 4 (Fig. 1c). A cross-reference of the three gene groups showed that 16 genes were included in the overlaps (Fig. 1d). A comparison of the 16 genes of mRNA expression demonstrated that WEE2 and CCNE1 were dramatically elevated in cluster 1 (P = 2.89085E-09 and P = 6.96046E-26, respectively) (Fig. 1e and Additional file 1: Table S3).
Expression of cell-cycle-associated genes were correlated with prognosis in GAC
To address the prognostic roles of genes with cluster-differential expression, we evenly divided the samples into two groups based on the mRNA z-scores from high values to low values. Survival curves according to each single gene are shown in Fig. 2. Comparison of the median survival and survival curves in the 16 genes between low-expression and high-expression groups revealed that ESPL1 and MCM5 were significantly associated with clinical outcome (P = 0.0448 and 0.0048, respectively) (Table 3). High expression of ESPL1 and MCM5 indicated a favorable median survival (18.33 vs 28.71 months and 19.94 vs 59.49 months, respectively). In addition, elevated CCNE2 and TTK expression also showed a trend of correlation with good prognosis, but with a non-significant difference (Table 3).
Expression of ESPL1 and MCM5 were associated with tumor stage progression
To detect the impacted factors of differential gene expression, we filtered out a sample set with intact pathologic records. According to investigating the expression of cluster 1-specific genes corresponding to each clinical manifestation, we detected obvious differences of the 16 genes expressed in each tumor stage (Fig. 3a). In addition, ESPL1 and MCM5 expression were negatively correlated with tumor stage progression (Fig. 3b, c), and regression analysis indicated that ESPL1 and MCM5 were correlated to tumor stages (Pearson correlation = − 0.25713; P = 9.54E-05 and Pearson correlation = − 0.13982; P = 0.023, respectively) (Additional file 1: Table S4).
Additionally, we investigated the correlation between cluster 1 specifically expressed genes and tumor stages. There were 15 genes, including ESPL1 and MCM5, that showed relativity (P < 0.05) with tumor stage changing (Fig. 3b) (Additional file 1: Table S4). A correlation matrix of each gene and stage showed that FZR1 was also weakly correlated with stages (Pearson correlation = − 0.20806; P = 0.001). What’s more, we detected that the expression of ESPL1 was positively correlated with the expression of FZR1, SFN, MCM5, CDC20 and CDC25C and negatively correlated with the expression of GADD45B and ANAPC13. MCM5 was positively correlated with the expression of CDC25C, ESPL1, FZR1 and CDC20 and negatively correlated with the expression of SMAD2, RBL2, GADD45B and ANAPC13 (Table 4) (Fig. 3c, d).
In summary, we studied the expression profiles of 122 cell-cycle-associated genes and sorted the GAC samples into 5 clusters according to mRNA z-score distribution. Sixteen specific genes showed differential expression between each cluster that were identified. We examined the prognostic roles of the 16 genes and detected that ESPL1 and MCM5 were significantly associated with overall survival. Moreover, we detected that the expression of ESPL1 and MCM5 were negatively correlated with tumor stage progression. These results suggested the significance of cell-cycle-related gene expression in the development and progression of GAC and provided potential targets for GAC therapy.
Uncontrolled tumor cell proliferation via aberrant expression of various cell cycle genes is one of the most essential features in multiple cancers. Therefore, cell-cycle-regulated genes are considerable targets in cancer therapy . Through this study of the mRNA expression of 122 cell-cycle-associated genes, we demonstrated a molecular profile, which defined five genomic clusters of GAC. Among them, cluster 1 manifested the worst survival and was characterized by dramatically elevated expression of WEE2 and CCNE1. WEE2 is an oocyte-specific protein tyrosine kinase that phosphorylates and inhibits CDK1, and acts as a key regulator of meiosis during both prophase I and metaphase II . In GAC, WEE2 might play a similar role by mediated the CDK1 phosphorylation. Amplification of CCNE1 (G1/S-specific cyclin-E1) is associated with poor outcome in breast, lung, and other solid cancers . However, prognosis implication analysis of single genes did not show any correlation of unfavorable outcome and highly expressed WEE2 and CCNE1 (P = 0.4349 and P = 0.2887, respectively). This discrepancy revealed that WEE2 and CCNE1 were indirectly regulated the GAC by triggering other molecular events.
In entire GAC cases, we detected that elevated ESPL1 and MCM5 expression were significantly associated favorable prognosis, and the expression variability were reversely consisted with stage progression. ESPL1, a protease (Separase) encoded gene, was reported overexpressing in mammary adenocarcinomas and related to tumor initiation and progression, but not mentioned in GAC . According to cleavage of cohesin complex subunit, ESPL1 might promote metaphase/anaphase transition during the cell cycle in GAC . MCM5, minichromosome maintenance complex component 5, which encoded protein is involved in the initiation of DNA replication and the effect might reduce the accumulation of genetic variation in GAC. The prognostic analysis of MCM5 in current study showed that the upregulated MCM5 was associated with a favorable prognosis. This result was in accordance with an previous immunohistochemical study, which suggested that patients with high MCM5 expression had significantly shorter survival times and MCM5 was associated with clinicopathological parameters in gastric adenocarcinoma . Besides, additional 13 genes specifically expressed in cluster 1 were also correlated to tumor stage. The regression analysis showed that ESPL1 and MCM5 was positively correlated (Person = 0.419389, P = 1.75479E-10) and both them were correlated to other cell cycle-specific genes expression, suggesting a co-regulation effect through genes expression changing in GAC.
Although previous studies revealed the significant role of cell cycle associated genes in gastric cancer, most works were focus on individual genes [23,24,25]. In current study, we discuss the panorama of cell cycle associated genes in GAC samples. Those results might appeal to the further studies of various cell cycle inhibitors with therapeutic potential. For example, evaluation of cell cycle derangement in thyroid tumors may serve as a useful tool for both DTC diagnosis and prognosis ; PLK1 and AK inhibitors display the potential for being employed in innovative therapeutic strategies for improving T-ALL patient outcome ; Inhibition of Aurora and Polo-like kinases suggest that targeting G2-M regulators may represent a novel approach for treatment of human .
Irrespective of tissue origin and adjuvant therapy, this work revealed the gene expression profile of cell cycle association in GAC. The impact and functions of distinctive genes need to be further investigated. We believe these results will facilitate the exploration of novel therapies, ultimately improving clinical outcome from this intractable disease.
American Joint Committee on Cancer
Differentiated thyroid carcinoma
Kyoto encyclopedia of genes and genomes
Whole transcriptome sequencing
T-cell acute lymphoblastic leukemia
The cancer genome atlas
Fuchs CS, Niedzwiecki D, Mamon HJ, Tepper JE, Ye X, Swanson RS, Enzinger PC, Haller DG, Dragovich T, Alberts SR, et al. Adjuvant Chemoradiotherapy With Epirubicin, Cisplatin, and Fluorouracil Compared With Adjuvant Chemoradiotherapy With Fluorouracil and Leucovorin After Curative Resection of Gastric Cancer: Results From CALGB 80101 (Alliance). J Clin Oncol Off J Am Soc Clin Oncol. 2017;35(32):JCO2017742130.
Van Cutsem E, Sagaert X, Topal B, Haustermans K, Prenen H. Gastric cancer. Lancet (London, England). 2016;388(10060):2654–64.
Torre LA, Siegel RL, Ward EM, Jemal A. Global Cancer incidence and mortality rates and trends--an update. Cancer Epidemiol Biomarkers Prev. 2016;25(1):16–27.
Ajani JA, Lee J, Sano T, Janjigian YY, Fan D, Song S. Gastric adenocarcinoma. Nat Rev Dis Prim. 2017;3:17036.
Tan IB, Ivanova T, Lim KH, Ong CW, Deng N, Lee J, Tan SH, Wu J, Lee MH, Ooi CH, et al. Intrinsic subtypes of gastric cancer, based on gene expression pattern, predict survival and respond differently to chemotherapy. Gastroenterology. 2011;141(2):476–85. 485 e471–411
Wang K, Kan J, Yuen ST, Shi ST, Chu KM, Law S, Chan TL, Kan Z, Chan AS, Tsui WY, et al. Exome sequencing identifies frequent mutation of ARID1A in molecular subtypes of gastric cancer. Nat Genet. 2011;43(12):1219–23.
Lei Z, Tan IB, Das K, Deng N, Zouridis H, Pattison S, Chua C, Feng Z, Guan YK, Ooi CH, et al. Identification of molecular subtypes of gastric cancer with different responses to PI3-kinase inhibitors and 5-fluorouracil. Gastroenterology. 2013;145(3):554–65.
Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144(5):646–74.
The Cancer Genome Atlas Research Network. Comprehensive molecular characterization of gastric adenocarcinoma. Nature. 2014;513 (7517):202–9.
Park JH, van Wyk H, Roxburgh CSD, Horgan PG, Edwards J, McMillan DC. Tumour invasiveness, the local and systemic environment and the basis of staging systems in colorectal cancer. Br J Cancer. 2017;116(11):1444–50.
Ge H, He X, Guo L, Yang X. Clinicopathological significance of HSP27 in gastric cancer: a meta-analysis. OncoTargets and therapy. 2017;10:4543–51.
Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, Jacobsen A, Byrne CJ, Heuer ML, Larsson E, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer discovery. 2012;2(5):401–4.
Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, Sun Y, Jacobsen A, Sinha R, Larsson E, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6(269):pl1.
Wang Y, Wang Y, Duan X, Wang Y, Zhang Z. Interleukin-1 receptor-associated kinase 1 correlates with metastasis and invasion in endometrial carcinoma. J Cell Biochem. 2017;119(3):2545–55.
de Hoon MJ, Imoto S, Nolan J, Miyano S. Open source clustering software. Bioinformatics (Oxford, England). 2004;20 (9):1453–4.
Saldanha AJ. Java Treeview--extensible visualization of microarray data. Bioinformatics (Oxford, England). 2004;20(17):3246–8.
Otto T, Sicinski P. Cell cycle proteins as promising targets in cancer therapy. Nat Rev Cancer. 2017;17(2):93–115.
Oh JS, Susor A, Conti M. Protein tyrosine kinase Wee1B is essential for metaphase II exit in mouse oocytes. Science (New York, NY). 2011;332(6028):462–5.
Etemadmoghadam D, Au-Yeung G, Wall M, Mitchell C, Kansara M, Loehrer E, Batzios C, George J, Ftouni S, Weir BA, et al. Resistance to CDK2 inhibitors is associated with selection of polyploid cells in CCNE1-amplified ovarian cancer. Clin Can Res. 2013;19(21):5960–71.
Mukherjee M, Ge G, Zhang N, Edwards DG, Sumazin P, Sharan SK, Rao PH, Medina D, Pati D. MMTV-Espl1 transgenic mice develop aneuploid, estrogen receptor alpha (ERalpha)-positive mammary adenocarcinomas. Oncogene. 2014;33(48):5511–22.
Kumar R. Separase: function beyond cohesion cleavage and an emerging oncogene. J Cell Biochem. 2017;118(6):1283–99.
Giaginis C, Giagini A, Tsourouflis G, Gatzidou E, Agapitos E, Kouraklis G, Theocharis S. MCM-2 and MCM-5 expression in gastric adenocarcinoma: clinical significance and comparison with Ki-67 proliferative marker. Dig Dis Sci. 2011;56(3):777–85.
Yang C, Wen Y, Li H, Zhang D, Zhang N, Shi X, Jiang B, Ma X, Yang P, Tang H, et al. Overexpression of minichromosome maintenance 2 predicts poor prognosis in patients with gastric cancer. Oncol Rep. 2012;27(1):135–42.
Zhao BW, Chen S, Li YF, Xiang J, Zhou ZW, Peng JS, Chen YB. Low expression of CDK10 correlates with adverse prognosis in gastric carcinoma. J Cancer. 2017;8(15):2907–14.
Baba K, Kitajima Y, Miyake S, Nakamura J, Wakiyama K, Sato H, Okuyama K, Kitagawa H, Tanaka T, Hiraki M, et al. Hypoxia-induced ANGPTL4 sustains tumour growth and anoikis resistance through different mechanisms in scirrhous gastric cancer cell lines. Sci Rep. 2017;7(1):11127.
Melck A, Masoudi H, Griffith OL, Rajput A, Wilkins G, Bugis S, Jones SJ, Wiseman SM. Cell cycle regulators show diagnostic and prognostic utility for differentiated thyroid cancer. Ann Surg Oncol. 2007;14(12):3403–11.
Sparta AM, Bressanin D, Chiarini F, Lonetti A, Cappellini A, Evangelisti C, Evangelisti C, Melchionda F, Pession A, Bertaina A, et al. Therapeutic targeting of Polo-like kinase-1 and Aurora kinases in T-cell acute lymphoblastic leukemia. Cell cycle (Georgetown, Tex). 2014;13(14):2237–47.
Markant SL, Esparza LA, Sun J, Barton KL, McCoig LM, Grant GA, Crawford JR, Levy ML, Northcott PA, Shih D, et al. Targeting sonic hedgehog-associated medulloblastoma through inhibition of aurora and polo-like kinases. Cancer Res. 2013;73(20):6310–22.
This work was funded by the National Natural Science Foundation of China (Nos. 81702355, 81372387, 81572820).
Availability of data and materials
The datasets analysed during the current study are available in the TCGA repository (http://www.cbioportal.org/).
Ethics approval and consent to participate
Data obtained from the TCGA open-access database was collected from tumors of patients who provided informed consent based on the guidelines from the TCGA Ethics, Law and Policy Group.
The authors declare that they no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Gene listed in KEGG cell cycle pathway (hsa04110). Table S2. Comparison of Survival Curves of each cluster by Log-rank (Mantel-Cox) test. Table S3. Comparison of individual gene expressions of each cluster. Table S4. Correlation between Cluster-specific genes expression and tumor stages. (DOCX 24 kb).
Figure S1. Networks of cell cycle associated genes. The interaction values between each connected gene were exported from String-db. (TIF 1159 kb).
About this article
Cite this article
Wang, D., Zhu, H., Guo, M. et al. Expression and prognostic value of cell-cycle-associated genes in gastric adenocarcinoma. BMC Gastroenterol 18, 81 (2018). https://doi.org/10.1186/s12876-018-0811-1
- Gastric adenocarcinoma
- Cell cycle
- Gene expression
- Tumor stages
- Prognostic implications