J Cancer 2020; 11(16):4841-4850. doi:10.7150/jca.46156 This issue
1. School of Bioscience and Bioengineering, South China University of Technology, Guangzhou, Guangdong, People's Republic of China.
2. The Second School of Clinical Medicine, Southern Medical University, Guangzhou, People's Republic of China.
3. Department of General Surgery, Guangdong Provincial People's Hospital, Guangdong Academy of Medical Sciences, Guangzhou, Guangdong, People's Republic of China.
4. School of Medicine, South China University of Technology, Guangzhou, Guangdong, People's Republic of China.
5. Shantou University Medical College, Shantou, Guangdong, People's Republic of China.
*These authors have contributed equally to this work.
Low expressions of PRKACB are related to the occurrence of various human malignancies. However, the prognostic value of PRKACB expression in colorectal cancer (CRC) patients remains controversial. In this analysis, PRKACB expression in CRC tumors was evaluated across the GEO, TCGA, and Oncomine databases, and a PRKACB survival analysis was performed based on the TCGA profile. We detected PRKACB in 7 GEO series (GSE110225, GSE32323, GSE44076, GSE9348, GSE41328, GSE21510, GSE68468) and TCGA spectra (all P <0.05). A meta-analysis performed in the Oncomine database revealed that PRKACB was significantly up-regulated in neoplastic tissues compared to normal tissues (all P <0.05). A Kaplan-Meier analysis demonstrated that lower PRKACB expression in tumors was significantly associated with poorer overall survival (OS) in patients with CRC (P <0.05). A subgroup analysis showed that low expression of PRKACB correlated with poor 1-, 3-, and 5-year OS (all P <0.05). Furthermore, in males (P = 0.0083), whites (P = 0.0463), and non-mucinous adenocarcinoma patients (P = 0.0108), the down-regulation of PRKACB expression was more significant for the OS prognostic value. Conclusion: PRKACB is down-regulated in tumors and associated with worsening OS in CRC patients.
Keywords: PRKACB, colorectal carcinoma, survival, therapeutic target
Colorectal cancer (CRC) is amongst the most prevalent digestive tract malignancies and it is the world's second most deadly cancer with almost 900 000 deaths recorded annually . Despite formidable advances in imaging techniques, surgery and multimodal therapy, the overall survival for patients with advanced CRC remains low. A study showed that the incidence of CRC in China has constantly been on the rise over the past three decades and equally predicted a further increase in the near future . According to the 2015 cancer statistics , it is estimated that about 2,376,300 new CRC cases and 191,000 CRC related deaths occurred in China, accounting for nearlyone-tenth of the global CRC burden. Even though outstanding radiological, surgical and multimodal therapeutic advances have recently been mad, the overall survival rates for patients with late stage CRC remain substantially low at ∼8-9% . A precise estimation of the prognosis plays an important role in the diagnostic and therapeutic management of CRC patients. Thus, identifying reliable and practical prognostic biomarkers and revealing novel CRC treatment targets is urgently required [5,6].
The protein kinase cAMP-dependent catalytic subunit β (PRKACB) encoding the cAMP-dependent protein kinase catalytic subunit β (PRKACB) is a member of the serine/threonine protein kinase family . It serves as a key regulatory point and is involved in cell proliferation, differentiation, regulation of apoptosis, and is closely related to physiological and pathological processes such as cell growth, gene expression, tumor proliferation and metastasis . Recently, the presence of the PRKACA and PRKACB fusion genes has been detected in various cancers such as bile duct cancers, fibrolamellar hepatocellular carcinoma, and pancreatic cancers, and there is a possibility that this might be implicated in the pathogenesis of cancer [9-11]. Some researches declared that PRKACB mutations could cause adrenal and bile duct tumors [12,13], and that high expressions of PRKACB is related to drug resistance in patients diagnosed with breast cancer . Previous studies have confirmed that the c-MYC was directly activated by PRKACB, which then induced neoplasia [15,16]. Conversely, another study suggested that the PRKACB might act as a tumor suppressor gene in non-small cell lung cancers . Nonetheless, the exact role played by PRKACB in colorectal cancer is still unknown. Thus, to evaluating the prognostic value of the PRKACB in CRC patients is fundamental.
To help elucidate the possible relationship between the PRKACB expression and CRC patient outcomes, we identified the PRKACB expression in GEO, Oncomine and TCGA databases and performed a survival analysis based on the TCGA profile, with hopes of providing useful insights into the development of colorectal cancer.
GEO microarray series (GSE110225, GSE32323, GSE44076, GSE9348, GSE41328, GSE21510, GSE68468) containing CRC tumor and non-tumor samples were obtained from the National Center for Biotechnology Information's (NCBI) Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/). Platforms and samples of GEO series were summarized in Table 1.
All of the publicly available colorectal cancer RNA-Seq data information were downloaded from The Cancer Genome Atlas (TCGA) official website (https://cancergenome.nih.gov/) before December 17, 2019, through the GDC Data Transfer Tool . TCGA RNAseq data was comprised of 568 tumor samples and 44 non-tumor samples.
Raw CEL files of the microarray from each GEO dataset were normalized by the quantile method of Robust Multichip Analysis (RMA) from the R affy package  and the normalized gene expression levels were presented as log2-transformed values by RMA. Gene expression of PRKACB was determined by tumor and non-tumor samples comparison from the R Limma package .
Studies comparing PRKACB between tumor and non-tumor samples in colorectal cancer were selected with a threshold of p-Value ≤ 1E-4, fold change ≥ 2 and top 10% gene rank in the Oncomine database (https://www.oncomine.org/).
The clinical data and PRKACB RNA Seq V2 data of 382 colorectal cancer patients (TCGA, Provisional) were downloaded from the cBioPortal database [23,24] (http://www.cBioPortal.org/), and the patients' information were matched with the sample ID by VLOOKUP index in EXCEL. After deleting some missing data, the PRKACB expression values were ranked from top to bottom and the median was taken as the cutoff point, these patients were then classified into a low expression group and high expression group to analyze the correlation between PRKACB expression with survival rates and clinical pathological characteristics.
Details of GEO series included in this analysis
|GEO series||Contributor(s), Year||Tumor||Non-tumor||Platform|
|GSE110225||Vlachavas E, 2018||30||30||[HG-U133A] Affymetrix Human Genome U133A Array[HG-U133_Plus_2] Affymetrix Human Genome U133 Plus 2.0 Array|
|GSE32323||Ahmed K et al., 2012||17||17||[HG-U133_Plus_2] Affymetrix Human Genome U133 Plus 2.0 Array|
|GSE44076||Solé X et al., 2014||98||98||[HG-U219] Affymetrix Human Genome U219 Array|
|GSE9348||Hong Y et al., 2010||70||12||[HG-U133_Plus_2] Affymetrix Human Genome U133 Plus 2.0 Array|
|GSE41328||Lin G et al., 2012||5||5||[HG-U133_Plus_2] Affymetrix Human Genome U133 Plus 2.0 Array|
|GSE21510||Tsukamoto S, 2011||19||25||[HG-U133_Plus_2] Affymetrix Human Genome U133 Plus 2.0 Array|
|GSE68468||NA, 2015||44||44||[HG-U133A] Affymetrix Human Genome U133A Array|
RNA-Seq data information of 568 tumor samples and 44 non-tumor samples were downloaded from TCGA official website and the edgeR package [21,22] was used to compare the mRNA expression of the tumor and non-tumor samples to screen for differential expression genes (DEGs). The P-value <0.05 and |logFC|>1 were chosen as cut-off criteria. The Spearman coefficients of the DEGs and PRKACB were calculated, whilst DEGs with p-value<0.05 were defined as PRKACB-related genes.
The Database for Annotation, Visualization, and Integrated Discovery (DAVID, http://david.ncifcrf.gov) (version 6.8)  is an online functional annotation tool that we applied for Gene Ontology (GO) enrichment analysis, including biological process (BP), cellular components (CC), and molecular function (MF). The DAVID database was also used to a perform pathway enrichment analysis. P-value <0.05 was considered as the threshold.
STRING (http://string-db.org) (version 11.0)  is an online biological database for the prediction of known and unknown protein interaction relationships. The PRKACB-related genes were uploaded to the STRING website to analyze the interactions between those proteins. The minimum required interaction score was set to 0.400 (medium confidence) and the protein nodes undergoing no interaction with other proteins were removed. Next, the PPI pairs were inputted into the Cytoscape software (http://www.cytoscape.org) (version 3.7.1)  to construct a PPI network and the top 10 hub genes were identified in accordance with the Cytoscape plug-in (degrees ranking of cytoHubba).
The mRNA expression of 7 complete samples from the Colorectal Adenocarcinoma (TCGA, Provisional) database was used to conduct an analysis by the Co-expression online analysis function in the cBioPort database (http://www.cBioPortal.org/) . P value <0.05 was considered as the threshold. Genes with Spearman correlation coefficients and PRKACB expression greater than 0.9 were screened and uploaded to Cytoscape software (http://www.cytoscape.org) (version 3.7.1)  to draw a gene co-expression network.
The details of the GEO series included in this analysis were summarized in Table 1. As illustrated in Figure 1, the expression of PRKACB in tumor samples was all significantly lower than non-tumor samples in GSE110225, GSE32323, GSE44076, GSE9348, GSE4128, GSE21510, GSE68468 and TCGA datasets (all P < 0.01, Figure 1).
mRNA expression levels of PRKACB between tumor and non-tumor samples in CRC patients on the GEO database series including GSE110225 (A), GSE32323 (B), GSE44076 (C), GSE9348 (D), GSE4128 (E), GSE21510 (F), GSE68468 (G), TCGA database (H).
For validation, we performed a meta-analysis of PRKACB expression in 15 analyses with a threshold set as p-Value ≤1E-4, fold change ≥ 2 and top 10% gene rank in the Oncomine database. As displayed in Figure 2 and compared with its value in normal tissues, PRKACB was significantly downregulated in rectal adenocarcinoma (P<0.0001, Figure 2B), rectal carcinoma (P<0.0001, Figure 2G), colon carcinoma (P<0.0001, Figure 2E-G), colon adenoma (P<0.0001, Figure 2E), colorectal carcinoma (P<0.0001, Figure 2C-D), and colorectal adenocarcinoma (P<0.0001, Figure 2D). Notwithstanding, the difference of PRKACB expression in rectal carcinoma has no value (P = 0.2021, Figure 2F). Besides, the PRKACB was also downregulated in cecal carcinoma (P<0.0005, Figure 2F) and cecal adenocarcinoma (P<0.0001, Figure 2G).
The results outlined that low expression of PRKACB in tumor tissues was considerably associated with poor disease-free survival (log rank P = 0.0025, Figure 3A) and overall survival (log rank P = 0.0186, Figure 3B) in patients with CRC. Moreover, a subgroup analysis revealed that the downregulation of PRKACB in tumor tissue was a risk factor for reduced 1 year (log rank P = 0.0100, HR = 9.273 (2.684-32.04), Figure 3C), 3 years (log rank P = 0.0041, HR = 3.477 (1.635-7.393), Figure 3D), and 5 years (log rank P = 0.0083, HR = 2.677 (1.367-5.244), Figure 3C) OS in CRC patients.
By the same token, we performed subgroup survival analyses in different populations. As shown in Figure 4, downregulation of PRKACB was associated with poorer survival in males and white patients, while no significant differences were found in females and black patients. Moreover, low PRKACB levels significantly contributed to worse OS in CRC patients without non-mucinous cancers (HR=2.557(1.328-4.924), log-rank P=0.0108, Figure 4E) and the down-regulation of PRKACB was observed in patients with stage III-IV colorectal cancer (HR = 2.931 (1.357-6.333), log-rank P = 0.0145, Figure 4G), but not identified as risk factors for patients with stage I-II colorectal cancer (P> 0.05, Figure 4H).
As delineated in Table 2, there are more male cases in the PRKACB low group (64.85% vs. 46.99%, P = 0.001) and PRKACB low group patients were significantly older than those in the PRKACB high group (66yr vs. 61.5yr, P = 0.022). However, BMI, tumor stage, lymph node stage, metastasis stage, AJCC stage, lymphovascular invasion, perineural invasion, vascular invasion, race category, person neoplasm status were not significantly different in PRKACB expression (P> 0.05).
The KEGG pathway enrichment of PRKACB interactive genes showed that neuroactive ligand-receptor interactions, pancreatic secretions, bile secretions, mineral absorption, salivary secretions, cAMP signaling pathways, glutamatergic synapse, GABAergic synapses, retrograde endocannabinoid signaling, circadian entrainment, etc. were the most enriched pathways (Figure 5A). Additionally, GO analysis results proved that PRKACB-related genes were significantly enriched in the complement activation classical pathway, proteolysis, complement activation, etc. at BP levels; extracellular region, extracellular space, plasma membrane, etc. at CC levels and antigen-binding, serine-type endopeptidase activity and hormone activity at MF levels (Figure 5B-D).
A total of 127 genes were filtered into the target genes PPI network complex, containing 100 nodes and 654 edges, 10 hub Genes (PRKACB, ATP2B2, MAPT, PHLPP2, ABCCB, GRIN2A, MYLK, GRIA1, BCHE, ADCY2) were screened according to Cytoscape 3.7.1 and its plug-in (Ranking degree of cytoHubba).
Identification of PRKACB co-expressed genes was completed by the cBioPortal database's online analysis function. A total of 182 genes with Spearman correlation coefficients greater than 0.9 (highly correlated) were selected and visualized via Cytoscape (version 3.7.1) (Figure 7). We found that FAM167A, NRIP3, RASL11B, ST13P4, TMEM99 completely correlated with PRKACB (Spearman's Correlation = 1) whereas ALPP, C3ORF70, JUND, and ZBTB7A were completely negative correlated with PRKACB (Spearman's Correlation = -1).
Over these past years, more and more data indicate that PRKACB is involved in the cancerization process of malignant tumors in different systems. PRKACB mutations and gene fusion of PRKACB can lead to adrenal, bile duct, liver, and pancreatic cancers [12,13], meanwhile down-regulation of PRKACB expression may be associated with poor survival in patients with non-small cell lung cancer. Makondi et al. found that targeting PRKACB may increase the responsiveness of colorectal tumors to irinotecan treatment , and Feng et al. also found pathological expressions of PRKACB in human colorectal cancer tissues infected with nucleobacteria . Consistent with previous studies , our study showed that PRKACB had a significant discrepancy of expression in normal and tumor tissues of colorectal cancer, and we also depicted that down-regulation of PRKACB expression is a risk factor for declining 1, 3, and 5-year survival rates in patients with CRC. This was equally associated with poor patient-free survival and overall survival. Consequently, we can confirm that PRKACB plays an important role in patients with colorectal cancer and affects the prognosis of patients.
Gene co-expression network analysis established that PRKACB had a completely positive correlation with FAM167A, NRIP3, RASL11B, ST13P4, TMEM99, and a completely negative correlation with ALPP, C3ORF70, JUND, and ZBTB7A. It has been confirmed by studies that NRIP3 and RASL11B play a role in suppressing cancer proliferation in breast cancer and renal cell carcinoma [31,32]. Meanwhile ALPP, JUND, ZBTB7A in gastric cancer, prostate cancer, breast cancer [33-35] and other cancers promote the progress of cancer. Thence, we can boldly speculate that PRKACB plays a synergistic role with these tumor suppressors in inhibiting tumor growth. On the contrary, oncogenes antagonize the tumor-suppressive effect of PRKACB.
Many studies have confirmed that microRNAs can regulate the expression of PRKACB [29,36-38], affect cell proliferation, adhesion, and metastasis, and eventually bring about the development of tumors. Also, studies have found that antisense oligodeoxynucleotides targeted to protein kinase subunits induce growth arrest, apoptosis, and differentiation in a variety of cancer cell lines both in vitro and in vivo . Similarly, we were able to regulate the development of colon cancer by targeting PRKACB, providing an important theoretical basis for the development of colon cancer chemotherapeutic drugs and new anti-EMT therapies.
Characteristics of CRC patients between PRKACB high and low groups
|High groups (n=166)||Low groups (n=165)||P|
|Age, median(IQR), years✱||61.5(21.75)||66(17)||0.022|
|BMI, kg/m2||<18.5 (n=4)||1||3||0.836|
|Tumor Stage||T1 (n=10)||4||6||0.849|
|Lymph Node Stage||N0 (n=188)||91||97||0.513|
|Metastasis Stage||M0 (n=228)||111||117||0.143|
|AJCC stage||I (n=52)||25||27||0.213|
|Lymphovascular invasion||YES (n=89)||43||46||0.910|
|Perineural Invasion||YES (n=51)||30||21||0.196|
|Vascular invasion||YES (n=61)||31||30||0.674|
|Race Category||WHITE (n=232)||112||120||0.479|
|BLACK OR AFRICAN AMERICAN (n=56)||28||28|
|AMERICAN INDIAN OR ALASKA NATIVE (n=1)||1||0|
|Person Neoplasm Status||TUMOR FREE (n=221)||115||106||0.845|
|WITH TUMOR (n=71)||36||35|
Comparison of PRKACB mRNA expression levels across 15 analyses in the Oncomine database. Meta-analysis of PRKACB expression in 15 analyses (A); PRKACB levels in Gaedcke Colorectal (B), Hong Colorectal (C), Skrzypczak Colorectal (D), Skrzypczak Colorectal 2 (E), Gaedcke Colorectal (F) and, TCGA Colorectal (G). Note: In the Kaiser Colon study, Colon Adenocarcinoma (n=41), Colon Mucinous Adenocarcinoma (n=13), Colon Signet Ring Cell Adenocarcinoma (n=2), Colon Small Cell Carcinoma (n=2), Rectosigmoid Adenocarcinoma (n=10), Rectosigmoid Mucinous Adenocarcinoma (n=2) were all merged into Colon Carcinoma; Rectal Adenocarcinoma (n=8), Rectal Mucinous Adenocarcinoma (n=4), Rectal Signet Ring Cell Adenocarcinoma (n=1) were merged into Rectal Carcinoma; In the TCGA Colorectal study, Colon Adenocarcinoma (n=101), Colon Mucinous Adenocarcinoma (n=22), Rectosigmoid Adenocarcinoma (n=3), Rectosigmoid Mucinous Adenocarcinoma (n=1) were merged into Colon Carcinoma; Rectal Adenocarcinoma (n=60), Rectal Mucinous Adenocarcinoma (n=6) were equally merged into Rectal Carcinoma.
Disease-free survival (A) and overall survival (B) of CRC patients grouped by PRKACB median cutoff in TCGA database; 1-year (C), 3-year (D) and 5-year (E) overall survivals comparison between high and low PRKACB groups.
Subgroup analyses of overall survival comparison in different populations [gender (A, B), race (C, D) and cancer type (E, F)] and different stage (G, H) with PRKACB median cutoffs in CRC patients.
The protein kinase cAMP-dependent catalytic subunit β (PRKACB) encoding the cAMP-dependent protein kinase catalytic subunit β (PRKACB) is a member of the serine/threonine protein kinase family . Activation of PKA catalytic activity is initiated by any signal causing an increase in intracellular cAMP concentration. The traditional view of PKA activation is that two cAMP molecules bind to each R subunit, causing a conformational change in the R subunit dimer, and the C subunits are released and become catalytically active through exposure of their active sites . PKs catalyze the transfer of phosphate groups onto Ser, Thr, or Tyr residues of target proteins. Phosphorylation of substrates represents a key regulatory mechanism in all eukaryotic cells41 and the various PKs target different substrates with a multitude of biological effects. In addition, studies have shown that PKA-mediated tamoxifen resistance in breast cancer is caused by 305 serine (S305) phosphorylation of Erα . Similarly, phosphorylation of PRKACB has also been found in neurological diseases . From a physiologic aspect, PRKACB may mainly regulate protein phosphorylation/dephosphorylation through the cAMP pathway. Unfortunately, we are yet to conduct experimental studies exploring the potential carcinogenic mechanism of PRKACB in the development of liver cancer. Even after taking into account previous reports, we are pragmatic on suggesting the hypothesis that low expression of PRKACB leads to poorer prognoses in CRC patients.
KEGG and GO biological function enrichment analyses of PRKACB related genes. KEGG signal pathway enrichment analysis (A); Biological process enrichment analysis (B), Cell component enrichment analysis (C) and molecular function enrichment analysis (D).
Hub genes of the PPI network. The darker the color, the bigger the degrees.
Construction of gene co-expression networks. Blue represents genes that are negatively related to PRKACB, and red points out genes that are positively related to PRKACB. The darker the color, the stronger the correlation.
This work was supported by grants from the Science and Technology Planning Project of Guangdong Province, China (No. 2017A030223006, 2016A020215128), the Science and Technology Planning Project of Guangzhou, China (No. 201704020077), the Research Fund of Guangdong General Hospital (No. y012015338), the Yuexiu Science and Information Center of Guangzhou Scientific Foundation (No. 2012‑GX‑046), the Second Batch of Scientific Research Projects of Dengfeng Plan (NO. DFJH201913), the Research Fund of CSCO-Roche Oncology (NO. Y-2019Roche-190) and the Research Fund of CSCO-Hansoh Oncology (NO. Y-HS2019/2-050).
The authors have declared that no competing interest exists.
1. Bray F, Ferlay J, Soerjomataram I. et al. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68:394-424
2. Zhang L, Cao F, Zhang GY. et al. Trends in and Predictions of Colorectal Cancer Incidence and Mortality in China From 1990 to 2025. Frontiers In Oncology. 2019;9:1-9
3. Chen W, Zheng R, Baade PD. et al. Cancer statistics in China, 2015. CA Cancer J Clin. 2016;66:115-132
4. Meyerhardt JA, Mayer RJ. Systemic therapy for colorectal cancer. N Engl J Med. 2005;352:87-476
5. Zhang F, Zhang Y, Zhao W. et al. Metabolomics for biomarker discovery in the diagnosis, prognosis, survival and recurrence of colorectal cancer: a systematic review. Oncotarget. 2017;8:35460-35472
6. Oh HH, Joo YE. Novel biomarkers for the diagnosis and prognosis of colorectal cancer. Intest Res. 2020;18:168-183
7. Zhu H, Lu W, Laurent C. et al. Genes encoding catalytic subunits of protein kinase A and risk of spina bifida. Birth Defects Res A Clin Mol Teratol. 2005;73:591-596
8. Skalhegg BS, Tasken K. Specificity in the cAMP/PKA signaling pathway. differential expression, regulation, and subcellular localization of subunits of PKA. Front Biosci. 1997;2:d331-342
9. Nakamura H, Arai Y, Totoki Y. et al. Genomic spectra of biliary tract cancer. Nat Genet. 2015;47:1003-1010
10. Singhi AD, Wood LD, Parks E. et al. Recurrent Rearrangements in PRKACA and PRKACB in Intraductal Oncocytic Papillary Neoplasms of the Pancreas and Bile Duct. Gastroenterology. 2019;158:573-582
11. Simon EP, Freije CA, Farber BA. et al. Transcriptomic characterization of fibrolamellar hepatocellular carcinoma. Proc Natl Acad Sci U S A. 2015;112:E5916-E5925
12. Espiard S, Knape MJ, Bathon K. et al. Activating PRKACB somatic mutation in cortisol-producing adenomas. JCI Insight. 2018;3:1-11
13. Loilome W, Juntana S, Namwat N. et al. PRKAR1A is overexpressed and represents a possible therapeutic target in human cholangiocarcinoma. Int J Cancer. 2011;129:34-44
14. Zwart W, Griekspoor A, Berno V. et al. PKA-induced resistance to tamoxifen is associated with an altered orientation of ERalpha towards co-activator SRC-1. Embo j. 2007;26:3534-3544
15. Wu KJ, Mattioli M, Morse HC 3rd. et al. c-MYC activates protein kinase A (PKA) by direct transcriptional activation of the PKA catalytic subunit beta (PKA-Cbeta) gene. Oncogene. 2002;21:7872-7882
16. Padmanabhan A, Li X, Bieberich CJ. Protein kinase A regulates MYC protein through transcriptional and post-translational mechanisms in a catalytic subunit isoform-specific manner. J Biol Chem. 2013;288:14158-14169
17. Chen Y, Gao Y, Tian Y. et al. PRKACB is downregulated in non-small cell lung cancer and exogenous PRKACB inhibits proliferation and invasion of LTEP-A2 cells. Oncol Lett. 2013;5:1803-1808
18. Wang Z, Jensen MA, Zenklusen JC. A Practical Guide to The Cancer Genome Atlas (TCGA). Methods Mol Biol. 2016;1418:111-141
19. Gautier L, Cope L, Bolstad BM. et al. affy-analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004;20:307-315
20. Ritchie ME, Phipson B, Wu D. et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43:e47-e47
21. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139-140
22. Nikolayeva O, Robinson MD. edgeR for differential RNA-seq and ChIP-seq analysis: an application to stem cell biology. Methods Mol Biol. 2014;1150:45-79
23. Gao J, Aksoy BA, Dogrusoz U. et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6:pl1-pl1
24. Cerami E, Gao J, Dogrusoz U. et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012;2:401-404
25. Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44-57
26. Szklarczyk D, Morris JH, Cook H. et al. The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. 2017;45:D362-D368
27. Killcoyne S, Carter GW, Smith J. et al. Cytoscape: a community-based framework for network modeling. Methods Mol Biol. 2009;563:219-239
28. Makondi PT, Chu CM, Wei PL. et al. Prediction of novel target genes and pathways involved in irinotecan-resistant colorectal cancer. PLoS One. 2017;12:1-18
29. Feng YY, Zeng DZ, Tong YN. et al. Alteration of microRNA-4474/4717 expression and CREB-binding protein in human colorectal cancer tissues infected with Fusobacterium nucleatum. PLoS One. 2019;14:1-13
30. Wu S, Wu F, Jiang Z. Identification of hub genes, key miRNAs and potential molecular mechanisms of colorectal cancer. Oncol Rep. 2017;38:2043-2050
31. He H, Dai J, Zhuo R. et al. Study on the mechanism behind lncRNA MEG3 affecting clear cell renal cell carcinoma by regulating miR-7/RASL11B signaling. J Cell Physiol. 2018;233:9503-9515
32. Dong-zhi C. Overexpression of NRIP3 gene inhibited metastasis in lymph node negative breast cancer. Journal of Shandong university (health science). 2019;50:83-85
33. Mao A, Chen M, Qin Q. et al. ZBTB7A promotes migration, invasion and metastasis of human breast cancer cells through NF-kappaB-induced epithelial-mesenchymal transition in vitro and in vivo. J Biochem. 2019;166:485-493
34. Liu S, Mao Q, Xue W. et al. High expression of ALPPL2 is associated with poor prognosis in gastric cancer. Hum Pathol. 2019;86:49-56
35. Elliott B, Millena AC, Matyunina L. et al. Essential role of JunD in cell proliferation is mediated via MYC signaling in prostate cancer cells. Cancer Lett. 2019;448:155-167
36. Irwandi RA, Khonsuphap P, Limlawan P. et al. miR-302a-3p regulates RANKL expression in human mandibular osteoblast-like cells. J Cell Biochem. 2018;119:4372-4381
37. Wang L, Liu J, Wang Q. et al. MicroRNA-200a-3p Mediates Neuroprotection in Alzheimer-Related Deficits and Attenuates Amyloid-Beta Overproduction and Tau Hyperphosphorylation via Coregulating BACE1 and PRKACB. Front Pharmacol. 2019;10:1-12
38. Wang J, Liu X, Wu H. et al. CREB up-regulates long non-coding RNA, HULC expression through interaction with microRNA-372 in liver cancer. Nucleic Acids Res. 2010;38:5366-5383
39. Nesterova MV, Johnson NR, Stewart T. et al. CpG immunomer DNA enhances antisense protein kinase A RIalpha inhibition of multidrug-resistant colon carcinoma growth in nude mice: molecular basis for combinatorial therapy. Clin Cancer Res. 2005;11:5950-5955
40. Sjoberg TJ, Kornev AP, Taylor SS. Dissecting the cAMP-inducible allosteric switch in protein kinase A RIalpha. Protein Sci. 2010;19:1213-1221
41. Manning G, Whyte DB, Martinez R. et al. The protein kinase complement of the human genome. Science. 2002;298:1912-1934
42. Michalides R, Griekspoor A, Balkenende A. et al. Tamoxifen resistance by a conformational arrest of the estrogen receptor alpha after PKA activation in breast cancer. Cancer Cell. 2004;5:597-605
Corresponding author: Xueqing Yao, MD, Ph.D., Department of General Surgery, Guangdong Provincial People′s Hospital, Guangdong Academy of Medical Sciences, Guangzhou 510080, China; Email: syyaoxueqingedu.cn.