Expression patterns and prognostic values of the cyclin-dependent kinase 1 and cyclin A2 gene cluster in pancreatic adenocarcinoma

Abstract

Objective

Pancreatic adenocarcinoma (PAAD) is one of the most lethal malignant tumors worldwide. Various studies based on cell lines, preclinical mouse models, and human tissue samples have shown that cell cycle-associated proteins are involved in the tumorigenesis and progression of PAAD.

Methods

Herein, we analyzed the relationships between CDK1 and CCNA2 gene expression and prognosis in patients with pancreatic cancer, using information from the Oncomine, cBioportal, Kaplan–Meier Plotter, and GEPIA databases.

Results

Expression levels of CDK1 and CCNA2 were significantly higher in PAAD compared with control tissues, and were associated with more advanced tumor stage. Survival analyses using the Kaplan–Meier Plotter database further confirmed that increased expression levels of CDK1 and CCNA2 were associated with a poor prognosis in patients with pancreatic cancer.

Conclusions

The results of this study suggest that CDK1 and CCNA2 may be potential therapeutic targets and prognostic biomarkers in patients with PAAD.

Keywords

CDK1 CCNA2 pancreatic adenocarcinoma gene expression prognosis database

Introduction

Pancreatic adenocarcinoma (PAAD) is one of the most lethal malignant tumors worldwide. According to GLOBOCAN 2012 estimates, more than 331,000 people die from pancreatic cancer each year, accounting for 4.0% of all cancer-related deaths. PAAD is ranked as the seventh leading cause of cancer-related deaths in both sexes.¹,² The incidence of pancreatic cancer is positively correlated with the level of economic development, and is higher in developed countries compared with developing countries. The 5-year survival rate for pancreatic cancer is estimated to be <5%.³ Despite extensive studies aimed at improving surgery and targeted therapies, the overall survival (OS) rate for patients with pancreatic cancer is extremely low, and the treatment of recurrent pancreatic cancer still presents huge challenges. There is therefore an urgent need to explore the mechanisms of pancreatic cancer and to identify new tumor biomarkers with high specificity and sensitivity.

Cyclin-dependent kinase (CDK) 1 is the only essential CDK⁴,⁵ required to facilitate the G2-M phase transition. It regulates G1 phase progression and G1-S transition.⁶ Alterations in CDK1 gene activity often lead to uncontrolled proliferation of tumor cells, which is a hallmark of malignant tumors. Cyclin A2 (CCNA2) is expressed in most human tissues and is a highly conserved member of the cyclin protein family. It plays a key role in controlling the G1-to-S and G2-to-M cell cycle transitions, in addition to being a key regulator of the hematopoietic lineage and in embryonic cells.⁷ Although dysregulated expression levels of CDK1 and CCNA2 genes and their links to patient prognosis and features have been reported in some pancreatic cancer studies,⁸ their expression profiles and prognostic relevance still remain unclear. In this study, we systematically evaluated the expression levels, potential functions, and prognostic relevance of CDK1 and CCNA2 in PAAD using a large number of data sets available in public databases.

Materials and methods

Ethics statement

This study was approved by the institutional review board of Shandong Provincial Third Hospital (Permit Number: No. 102) and was conducted according to the ethical guidelines of the Declaration of Helsinki. All the datasets were extracted from online public databases, for which written informed consent had been obtained.

Oncomine analysis

Oncomine (www.oncomine.org) is a public database that can be used to analyze gene expression in various cancers. It uses Student’s t-tests for two-class differential expression analyses. In the current study, we entered the names of the genes of interest into Oncomine to analyze their expression levels in pancreatic cancer. The expression levels of CDK1 and CCNA2 in cancer specimens were compared with those in normal controls. The cut-off of P-value was set at 0.05.

GEPIA dataset

GEPIA (Gene Expression Profiling Interactive Analysis) is an online database that combines TCGA and GTEx data to compile RNA-seq results from 9,736 tumor and 8,587 normal samples, to allow assessment of differential expression profiles, patient outcomes, and other analyses.⁹ We used these public bioinformatics platforms to evaluate the expression profiles of CDK1 and CCNA2 in pancreatic cancer.

UALCAN database

The UALCAN database (http://ualcan.path.uab.edu) is an interactive web resource including RNA-seq and clinical data for 31 cancer types from the TCGA database. It can be used to compare expression of the gene of interest between tumor and normal samples, and to analyze the relationships between gene expression and clinicopathologic parameters.¹⁰ In this study, we used UALCAN to analyze the mRNA expression levels of CDK1 and CCNA2 in pancreatic cancer tissues and their association with clinicopathologic parameters. The significance threshold was P < 0.05.

RNA extraction and quantitative reverse transcription-polymerase chain reaction (PCR)

We verified the mRNA expression levels of CDK1 and CCNA2 in 10 pairs of PAAD tissues and adjacent non-tumor tissues collected during surgery from patients at Shandong Provincial Third Hospital, between June and November 2019. Written informed consent was obtained from all participating patients. Approximately 20 mg of tissue was obtained from each sample and subjected to RNA extraction using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s protocols. All PCR reactions were performed using an ABI Prism 5700 Sequence Detection System (PerkinElmer Applied Biosystems, Foster City, CA, USA), according to the manufacturer’s instructions.

Human Protein Atlas

The Human Protein Atlas database (https://www.proteinatlas.org) contains immunohistochemical expression data for nearly 20 common kinds of tumors, including 12 individual tumor subtypes for each tumor.¹¹ The database allows researchers to identify the expression patterns of certain proteins in particular types of tumors. We directly compared CDK1 and CCNA2 expression levels between normal and pancreatic cancer tissues based on immunohistochemistry images.

Prognostic value analyses

We assessed the prognostic relevance of CDK1 and CCNA2 expression in pancreatic cancer using the OncoLnc (www.oncolnc.org) and GEPIA databases. For each gene, the median expression value was used to stratify patients into low- and high-expressing groups, and the OS of pancreatic cancer patients was assessed, with corresponding 95% confidence intervals and P-values. P < 0.05 was considered statically significant.

Functional pathway enrichment analyses

The Gene Ontology (GO) resource allows the identification of attributes that are significantly overrepresented in a set of genes, including cellular component, molecular function, and biological process attributes.¹² In addition, the Kyoto Encyclopedia of Genes and Genomes (KEGG) database is used to understand the high-level functions and utilities of biological systems. The online STRING (https://string-db.org/) tool is used to provide investigators with systematic and comprehensive functional annotation tools to reveal the biological meaning behind an extensive list of genes. In our study, GO and KEGG enrichment analyses were conducted for CDK1 and CCNA2 using STRING. The significance threshold was P < 0.05.

Results

CDK1 and CCNA2 gene expression in various cancers

We compared CDK1 and CCNA2 gene expression levels in all cancer and normal tissues using the Oncomine database. CDK1 mRNA levels were significantly elevated in pancreatic cancer in two datasets, and CCNA2 was significantly upregulated in one pancreatic cancer dataset (P < 0.05) (Figure 1a). CDK1 and CCNA2 were expressed in most cancers, indicating that these genes may be cancer-promoting genes. We further queried the Oncomine database to analyze the expression levels of CDK1 and CCNA2 in several pancreatic cancer studies, and found that the transcription levels of both genes were significantly higher in pancreatic cancer compared with normal pancreatic tissues (P < 0.05) (Figure 1b).

Figure 1.

Transcription levels of CDK1 and CCNA2 in different types of cancers (Oncomine). (a) Expression levels of CDK1 and CCNA2 in different types of cancers. Red, over-expression; blue, under-expression. (b) Expression levels of CDK1 and CCNA2 in several pancreatic cancer studies. Red, over-expression. Significance threshold P < 0.05.

Relationship between CDK1 and CCNA2 expression levels and clinicopathological findings in patients with PAAD

Expression levels of CDK1 and CCNA2 were significantly higher in PAAD compared with normal control tissues, according to GEPIA (P < 0.05) (Figure 2a, b). There was also a clear correlation between CDK1 and CCNA2 expression and PAAD tumor grade, with more advanced grades associated with higher expression of CDK1 and CCNA2 (Figure 2c). We also examined CDK1 and CCNA2 protein expression patterns in PAAD using the Human Protein Atlas, and confirmed that CDK1 and CCNA2 protein levels were elevated in PAAD relative to normal controls (Figure 3). We further verified the mRNA expression levels of CDK1 and CCNA2 in 10 pairs of PAAD and adjacent non-tumor tissues collected during surgery from patients at Shandong Provincial Third Hospital (aged 55–70 years, 5 male and 5 female), and showed that CDK1 and CCNA2 mRNA levels were significantly higher in PAAD compared with normal tissues (P < 0.05) (Appendix 1). Taken together, our results showed that transcriptional and proteomic expression levels of CDK1 and CCNA2 were increased in patients with PAAD.

Figure 2.

Correlation between CDK1/CCNA2 expression and tumor stage in patients with pancreatic adenocarcinoma (PAAD). Expression profiles of CDK1 (a) and CCNA2 (b) in patients with PAAD. (c) Correlation between CDK1 and CCNA2 expression levels and tumor stage in patients with PAAD (UALCAN database). X-axis indicates PAAD tumor stage; y-axis indicates expression level of CDK1/CCNA2. *P < 0.05 compared with normal control. PAAD, pancreatic adenocarcinoma; T, tumor; N, normal.

Figure 3.

CDK1 and CCNA2 expression in pancreatic adenocarcinoma (PAAD) and normal tissues from the Human Protein Atlas. CDK1 and CCNA2 proteins were highly expressed in PAAD but not in normal pancreatic tissues, according to immunohistochemistry. PAAD, pancreatic adenocarcinoma.

Prognostic values of CDK1 and CCNA2 in PAAD patients

We evaluated the prognostic values of CDK1 and CCNA2 in PAAD patients using the GEPIA and OncoLnc databases, and showed that increased expression levels of both genes were significantly associated with poor OS (P = 0.001 and P = 0.025, respectively) (Figure 4).

Figure 4.

Prognostic values of CDK1 and CCNA2 expression in patients with pancreatic adenocarcinoma. Overall survival curves for patients with high (red) and low (blue) CDK1 and CCNA2 expression levels were plotted using the GEPIA and OncoLnc databases, with a threshold of P < 0.05. TPM, transcripts per million.

Functional enrichment of CDK1 and CCNA2 in PAAD patients

We explored the functions of CDK1 and CCNA2 by GO and KEGG analyses. According to GO analyses, CDK1 and CCNA2 and their neighboring genes were primarily enriched for signal transduction involved in DNA damage checkpoint, histone phosphorylation, cellular protein metabolic process, protein phosphorylation, G2/M transition of mitotic cell cycle, mitotic cell cycle checkpoint, nuclear division, negative regulation of mitotic cell cycle, regulation of cell cycle, and cell division (Figure 5a). Enriched molecular functions included transcription regulation by protein domain-specific binding, enzyme binding, protein binding, enzyme regulator activity, ATP binding, protein kinase binding, cyclin binding, histone kinase activity, anaphase-promoting complex binding, and cyclin-dependent protein serine/threonine kinase activity (Figure 5b). Cellular component annotations for these genes included perinuclear region of cytoplasm, mitotic spindle, nucleoplasm, cyclin B1-CDK1 complex, transferase complex, cytosol, microtubule organizing center, spindle, chromosomal region, and condensed chromosome (Figure 5c). The KEGG pathways for these two genes are shown in Table 1. Among these pathways, cell cycle, oocyte meiosis, cellular senescence, p53 signaling pathway, viral carcinogenesis, human T-lymphotropic virus type-I infection, FoxO signaling pathway, Epstein–Barr virus infection, phosphoinositide 3-kinase-Akt signaling pathway, and pathways in cancer were involved in the tumor development and pathogenesis of PAAD (Figure 6a, b).

Figure 5.

Functional Enrichment Analysis of CDK1 and CCNA2 in patients with pancreatic adenocarcinoma. GO enrichment analysis predicted the functions of the target genes from three aspects: biological processes (a), cellular components (b), and molecular functions (c).

Table 1.

KEGG pathway analysis of CDK1 and CCNA2 in pancreatic cancer.

Pathway ID	Pathway name	Gene count	False discovery rate	Genes
Hsa04110	Cell cycle	12	1.61E-25	BUB1, BUB1B, CCNA2, CCNB1, CCNB2, CCNE1, CDC20, CDK1, CDK2, MAD2L1, PLK1, PTTG1
Hsa04114	Oocyte meiosis	10	6.36E-20	BUB1, CCNB1, CCNB2, CCNE1, CDC20, CDK1, CDK2, MAD2L1, PLK1, PTTG1
Hsa04218	Cellular senescence	6	1.49E-09	CCNA2, CCNB1, CCNB2, CCNE1, CDK1, CDK2
Hsa04115	p53 signaling pathway	5	2.24E-09	CCNB1, CCNB2, CCNE1, CDK1, CDK2
Hsa05203	Viral carcinogenesis	5	2.23E-07	CCNA2, CCNE1, CDC20, CDK1, CDK2
Hsa05166	HTLV-I infection	5	8.71E-07	BUB1B, CCNB2, CDC20, MAD2L1, PTTG1
Hsa04068	FoxO signaling pathway	4	2.86E-06	CCNB1, CCNB2, CDK2, PLK1
Hsa05169	Epstein–Barr virus infection	3	0.00047	CCNA2, CDK1, CDK2
Hsa04151	PI3K-Akt signaling pathway	2	0.0239	CCNE1, CDK2
hsa05200	Pathways in cancer	2	0.0467	CCNE1, CDK2

HTLV, human T-lymphotropic virus; PI3K, phosphoinositide 3-kinase.

Figure 6.

p53 signal pathway and cell cycle pathway regulated by CDK1 and CCNA2 in pancreatic adenocarcinoma (PAAD). p53 signal pathway (a) and cell cycle pathway (b) regulated by the CDK1 and CCNA2 gene cluster in PAAD.

Discussion

PAAD remains an important unresolved global health problem, largely because of a lack of early diagnosis and effective treatment targets.¹³ Stratifying patients based on biological markers and exploring the prognostic roles of such biomarkers can help to identify novel therapeutic approaches and appropriate treatment strategies. In the current study, we investigated the mRNA expression patterns and prognostic values of CDK1 and CCNA2 in PAAD. These results will support efforts to improve the diagnosis, prognostic analysis, and treatment of patients with PAAD.¹⁴,¹⁵

CDK1 is a cell cycle regulator associated with tumor growth.⁴ CDK1 can execute all the events required for cell cycle division, indicating that it is sufficient for the proliferation of mammalian cells. CDK1 expression was shown to be significantly increased in colorectal cancer, hepatocellular carcinoma, and human melanoma,^16–18 while another study found that up-regulation of CDK1 expression was associated with a poor prognosis in patients with breast cancer and epithelial ovarian cancer.^19–21 Furthermore, CDK1 inhibitors significantly reduced the proportion of cells in the S and G2/M stages in a pancreatic ductal adenocarcinoma (PDAC) tumor cell model, and inhibition of CDK1 expression also showed promising anti-cancer activity in pancreatic cancer cells.²² All these findings were consistent with the current results. We found that CDK1 expression was significantly increased in patients with pancreatic cancer, and its up-regulation was significantly associated with a poor prognosis. CCNA2 is a member of the highly conserved cyclin family, and has been found to be up-regulated in a variety of tumors, indicating potential roles in tumor progression and transformation.²³ High expression of CCNA2 was associated with increased proliferation of hepatocellular and breast cancer cells, which might be helpful in monitoring the efficacy of chemotherapy.²⁴,²⁵ Bioinformatics analysis showed that over-expression of CCNA2 was associated with the development of PDAC, which was consistent with our results. We demonstrated that over-expression of CCNA2 was significantly associated with a poor prognosis in patients with pancreatic cancer. Based on the above findings, we concluded that the CDK1 and CCNA2 gene cluster may provide a novel combined prognostic biomarker for PAAD.

This study had several limitations. First, the data used in this study were provided by various public databases, and were not generated by us. Second, the results were web-based and were not verified by biological experiments. Further mechanistic studies based on our findings are therefore needed.

In conclusion, this study provides a comprehensive bioinformatics analysis of the roles of CDK1 and CCNA2 in the progression and development of PAAD. The results suggest that the CDK1 and CCNA2 gene cluster may provide potential therapeutic targets for PAAD, and as well as biomarkers for identifying high-risk populations among pancreatic cancer patients. These results might also contribute to further studies of the therapeutic potential of various cell cycle inhibitors.

Supplemental Material

sj-pdf-1-imr-10.1177_0300060520930113 - Supplemental material for Expression patterns and prognostic values of the cyclin-dependent kinase 1 and cyclin A2 gene cluster in pancreatic adenocarcinoma

Supplemental material, sj-pdf-1-imr-10.1177_0300060520930113 for Expression patterns and prognostic values of the cyclin-dependent kinase 1 and cyclin A2 gene cluster in pancreatic adenocarcinoma by Peng Jiang, Ming Zhang, Liangliang Gui and Kai Zhang in Journal of International Medical Research

Footnotes

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Declaration of conflicting interest

The authors declare that there is no conflict of interest.

Funding

This research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors.

ORCID iD

Kai Zhang

Supplemental material

Supplemental material for this article is available online.

Appendix 1.

CDK1 and CCNA2 mRNA expression levels in 10 pairs of pancreatic adenocarcinoma and adjacent non-tumor tissues.

References

Ferlay

Soerjomataram

Dikshit

, et al. Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer 2015; 136: E359–E386. DOI: 10.1002/ijc.29210.

Hidalgo

Cascinu

Kleeff

, et al. Addressing the challenges of pancreatic cancer: future directions for improving outcomes. Pancreatology 2015; 15: 8–18. DOI: 10.1016/j.pan.2014.10.001.

Vincent

Herman

Schulick

, et al. Pancreatic cancer. Lancet 2011; 378: 607–620. DOI: 10.1016/S0140-6736(10)62307-0.

Santamaria

Barriere

Cerqueira

, et al. Cdk1 is sufficient to drive the mammalian cell cycle. Nature 2007; 448: 811–815. DOI: 10.1038/nature06046.

Malumbres

Barbacid

Cell cycle, CDKs and cancer: a changing paradigm.

Nat Rev Cancer 2009; 9: 153–166. DOI: 10.1038/nrc2602.

Enserink

Kolodner

RD.

An overview of Cdk1-controlled targets and processes.

Cell Div 2010; 5: 11. DOI: 10.1186/1747-1028-5-11.

Arsic

Bendris

Peter

, et al. A novel function for Cyclin A2: control of cell invasion via RhoA signaling. J Cell Biol 2012; 196: 147–162. DOI: 10.1083/jcb.201102085.

Dong

Huang

Zhang

, et al. Overexpression of BUB1B, CCNA2, CDC20, and CDK1 in tumor tissues predicts poor survival in pancreatic ductal adenocarcinoma. Biosci Rep 2019; 39: BSR20182306. DOI: 10.1042/BSR20182306.

Tang

Kang

, et al. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res 2017; 45: W98–W102. DOI: 10.1093/nar/gkx247.

10.

Chandrashekar

Bashel

Balasubramanya

SAH

, et al. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia 2017; 19: 649–658. DOI: 10.1016/j.neo.2017.05.002.

11.

Asplund

Edqvist

Schwenk

, et al. Antibodies for profiling the human proteome-The Human Protein Atlas as a resource for cancer research. Proteomics 2012; 12: 2067–2077. DOI: 10.1002/pmic.201100504.

12.

Ashburner

Ball

Blake

, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet 2000; 25: 25–29. DOI: 10.1038/75556.

13.

Seufferlein

Bachet

Van Cutsem

, et al. Pancreatic adenocarcinoma: ESMO-ESDO Clinical Practice Guidelines for diagnosis, treatment and follow-up. Ann Oncol 2012; 23: vii33–vii40. DOI: 10.1093/annonc/mds224.

14.

Dunne

Hezel

AF.

Genetics and biology of pancreatic ductal adenocarcinoma.

Hematol Oncol Clin North Am 2015; 29: 595–608. DOI: 10.1016/j.hoc.2015.04.003.

15.

Ying

Dey

Yao

, et al. Genetics and biology of pancreatic ductal adenocarcinoma. Genes Dev 2016; 30: 355–385. DOI: 10.1101/gad.275776.115.

16.

Thorenoor

Faltejskova-Vychytilova

Hombach

, et al. Long non-coding RNA ZFAS1 interacts with CDK1 and is involved in p53-dependent cell cycle control and apoptosis in colorectal cancer. Oncotarget 2016; 7: 622–637. DOI: 10.18632/oncotarget.5807.

17.

Ravindran Menon

Luo

Arcaroli

, et al. CDK1 interacts with Sox2 and promotes tumor initiation in human melanoma. Cancer Res 2018; 78: 6561–6574. DOI: 10.1158/0008-5472.CAN-18-0330.

18.

Wang

Chok

, et al. Blocking CDK1/PDK1/beta-Catenin signaling by CDK1 inhibitor RO3306 increased the efficacy of sorafenib treatment by targeting cancer stem cells in a preclinical model of hepatocellular carcinoma. Theranostics 2018; 8: 3737–3750. DOI: 10.7150/thno.25487.

19.

Chen

Xie

, et al. Association study of germline variants in CCNB1 and CDK1 with breast cancer susceptibility, progression, and survival among Chinese Han women. PLoS One 2013; 8: e84489. DOI: 10.1371/journal.pone.0084489.

20.

Yang

Cho

Shin

, et al. Accumulation of cytoplasmic Cdk1 is associated with cancer growth and survival rate in epithelial ovarian cancer. Oncotarget 2016; 7: 49481–49497. DOI: 10.18632/oncotarget.10373.

21.

Huang

Wang

, et al. The expression of CDK1 is associated with proliferation and can be a prognostic factor in epithelial ovarian cancer. Tumour Biol 2015; 36: 4939–4948. DOI: 10.1007/s13277-015-3141-8.

22.

Piao

Zhu

Sun

, et al. High expression of CDK1 and BUB1 predicts poor prognosis of pancreatic ductal adenocarcinoma. Gene2019; 701: 15–22. DOI: 10.1016/j.gene.2019.02.081.

23.

Uhlen

Oksvold

Fagerberg

, et al. Towards a knowledge-based Human Protein Atlas. Nat Biotechnol 2010; 28: 1248–1250. DOI: 10.1038/nbt1210-1248.

24.

Yang

Gong

Wang

, et al. Waltonitone inhibits proliferation of hepatoma cells and tumorigenesis via FXR-miR-22-CCNA2 signaling pathway. Oncotarget 2016; 7: 75165–75175. DOI: 10.18632/oncotarget.12614.

25.

Gao

Han

, et al. CCNA2 is a prognostic biomarker for ER+ breast cancer and tamoxifen resistance. PLoS One 2014; 9: e91771. DOI: 10.1371/journal.pone.0091771.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

0.08 MB