- Original article
- Open Access
Integrative analysis of gene expression profiles reveals specific signaling pathways associated with pancreatic duct adenocarcinoma
Cancer Communications volume 38, Article number: 13 (2018)
Pancreatic duct adenocarcinoma (PDAC) remains a major health problem because conventional cancer treatments are relatively ineffective against it. Microarray studies have linked many genes to pancreatic cancer, but the available data have not been extensively mined for potential insights into PDAC. This study attempted to identify PDAC-associated genes and signaling pathways based on six microarray-based profiles of gene expression in pancreatic cancer deposited in the gene expression omnibus database.
Pathway network methods were used to analyze core pathways in six publicly available pancreatic cancer gene (GSE71989, GSE15471, GSE16515, GSE32676, GSE41368 and GSE28735) expression profiles. Genes potentially linked to PDAC were assessed for potential impact on survival time based on data in The Cancer Genome Atlas and International Cancer Genome Consortium databases, and the expression of one candidate gene (CKS2) and its association with survival was examined in 102 patients with PDAC from our hospital. Effects of CKS2 knockdown were explored in the PDAC cell lines BxPC-3 and CFPAC-1.
The KEGG signaling pathway called “pathway in cancer” may play an important role in pancreatic cancer development and progression. Five genes (BIRC5, CKS2, ITGA3, ITGA6 and RALA) in this pathway were significantly associated with survival time in patients with PDAC. CKS2 was overexpressed in PDAC samples from our hospital, and higher CKS2 expression in these patients was associated with shorter survival time. CKS2 knockdown substantially inhibited PDAC cell proliferation in vitro.
Analysis integrating existing microarray datasets allowed identification of the “pathway in cancer” as an important signaling pathway in PDAC. This integrative approach may be powerful for identifying genes and pathways involved in cancer.
Pancreatic ductal adenocarcinoma (PDAC), which accounts for approximately 90% of all pancreatic tumors , is the ninth leading cause of cancer-related deaths in China . Due to the lack of effective methods for early diagnosis and therapy, pancreatic cancer has the second highest ratio of mortality to incidence (approximately 88%) among all types of cancer . Microarray platforms have produced numerous datasets of gene expression profiles in pancreatic cancer, which are publicly available through the gene expression omnibus (GEO) of the National Center for Biotechnology Information . Studies of such data have identified several genes associated with pancreatic cancer, but these datasets have yet to be mined extensively to identify genes and pathways linked to PDAC.
To identify additional genes that may be associated with PDAC, we mined six frequently cited profiles of genes expressed in pancreatic cancer that have been deposited in the GEO database. We used pathway network analysis to identify core pathways that may affect pancreatic cancer development and progression. Our analysis of genes differentially expressed in pancreatic cancer took into account expression quantitative trait loci , copy number variation [6, 7], and DNA methylation , all of which can affect gene expression levels.
Affymetrix microarray data
The following pancreatic cancer gene expression profiles were downloaded from the GEO database (http://www.ncbi.nlm.nih.gov/geo/): GSE71989, GSE15471, GSE16515, and GSE32676, which were obtained using an Affymetrix GeneChip® Human Genome U133 Plus 2.0 Array; as well as GSE41368 and GSE28735, which were obtained using an Affymetrix Human Gene 1.0 ST Array.
We downloaded the probe-level data from the GEO database as CEL files, then we used the GCBI on-line analysis tool (http://www.gcbi.com.cn) to normalize and log2-transform the data. Genes differentially expressed between tumor samples and non-tumor samples were filtered by requiring a fold difference of at least two and a maximum false discovery rate of 0.05. Pathway analysis was performed using DAVID software (https://david.ncifcrf.gov) as described .
Core pathway analysis
Pathways enriched for genes differentially expressed between tumor and non-tumor tissues were identified and connected in a pathway network to reveal relationships among the pathways. A pathway network was constructed using the GCBI on-line analysis tool, and for each pathway in the network, the number of upstream pathways (in-degree) and the number of downstream pathways (out-degree) were determined in order to assess the degree centrality of that pathway in the network. Greater degree centrality means that a pathway regulates or is regulated by more pathways, implying a more important role in the network.
CKS2 analysis based on database mining
Copy number variation of CKS2 in various types of cancer was analyzed using the cBioPortal database (http://www.cbioportal.org), and CKS2 overexpression in cancer was analyzed using the Oncomine database (http://www.omcomine.org). CKS2 expression levels in various cancer cell lines were analyzed using the CCLE database (https://portals.broadinstitute.org/ccle/home).
PDAC cell culture
Human PDAC cell lines BxPC-3 and CFPAC-1 were purchased from the Cell Bank of Type Culture Collection (Chinese Academy of Sciences, Shanghai Institute of Biochemistry and Cell Biology). The human immortalized pancreatic duct epithelial cell line HPDE6-C7 was purchased from Biotechnology Company. All cell lines were characterized by DNA fingerprinting and tested for mycoplasma and found to be free of infection.
BxPC-3 and HPDE6-C7 cells were cultured in RPMI-1640 medium (Gibco, Thermo Fisher Scientific, USA) supplemented with 10% fetal bovine serum (FBS; Hyclone, USA), while CFPAC-1 cells were cultured in IMDM (Gibco, Thermo Fisher Scientific) supplemented with 10% FBS. All cell lines were grown without antibiotics in an atmosphere of 5% CO2 and 99% relative humidity at 37 °C.
Survival analysis of key genes
SurvExpress  was used to analyze survival associated with key genes found in this study.
Surgical biopsies of tumor tissues and paired non-tumor tissues from patients definitively diagnosed with PDAC based on histology were obtained from the Department of Hepatobiliary Surgery at Sun Yat-sen Memorial Hospital (Guangzhou, China). Tissues were preserved in liquid nitrogen immediately following surgical removal, paraffin-embedded and stored properly before use in the present study. None of the patients had received anticancer treatment before biopsy collection. Clinicopathology and follow-up data for patients are shown in Additional file 1: Table S1.
Western blotting for CKS2
Total protein (20 μg) was extracted from cells or patient tissue samples, resolved by SDS-PAGE and transferred to PVDF membranes (Millipore, Germany). Membranes were incubated with antibodies against CKS2 (ab155078; Abcam, Cambridge, UK) or β-actin (ab8227; Abcam), then with horseradish peroxidase-conjugated secondary antibody from the Phototope-HRP Western blot detection kit (Cell Signaling Technology, USA). Protein bands were visualized using the components in this kit.
CKS2 knockdown in PDAC cells
Short interfering RNA (siRNA) targeting CKS2 and non-targeting control siRNA (Additional file 2: Table S2) were purchased from GenePharma (Suzhou, China). PDAC cells were transfected with siRNA using Lipofectamine 2000 (Life Technologies, USA) according to the manufacturer’s instructions.
Proliferation of PDAC cell lines
Cells were seeded into 96-well plates (2 × 103 cells in 100 µL per well) and cultured for certain periods. Then cell viability was measured using a CCK-8 assay (Dojindo, Japan). Samples were analyzed as six replicates per experiment, and three independent experiments were performed.
RNA extraction and qRT-PCR analysis
Total RNA from pancreatic cancer cell lines was extracted using TRIzol regent (Invitrogen, USA). First-strand cDNA was synthesized using the PrimeScript™ RT reagent Kit (TaKaRa, Japan) and amplified by qRT-PCR on an ABI Prism 7900 sequence detection system (Applied Biosystems) using SYBR Green. Sequences of gene-specific primers are listed in Additional file 3: Table S3.
All statistical analyses were performed using Graphpad Prism 5.0 and P < 0.05 was considered significant. Differences between two groups were assessed for significance using Student’s t test. Kaplan–Meier survival estimates were plotted, and inter-group differences were assessed using the log-rank test.
Genes differentially expressed in pancreatic cancer based on integrative analysis of GEO datasets
We examined six gene expression profiles from the GEO database (Additional file 4: Table S4) in order to identify genes differentially expressed between PDAC and normal tissues. One sample each from GSE15471 to GSE28735 were excluded from further analysis because they did not meet quality controls. Analysis of the remaining samples indicated a larger number of genes up-regulated in PDAC than down-regulated (Fig. 1).
The “pathway in cancer” is the top core pathway in pancreatic cancer
Genes differentially expressed between PDAC and non-tumor tissues in each of the six profiles were assembled into functional pathways. As expected, the pathways that were enriched in differentially expressed genes varied among the six profiles (data not shown). We identified 32 significant pathways common to the six profiles (Fig. 2a). We also used network analysis to identify pathways common to the six profiles that were likely to be most important in PDAC, which led to 13 common core pathways (Fig. 2b): pathway in cancer, p53 signaling pathway, focal adhesion, cytokine–cytokine receptor interaction, regulation of actin cytoskeleton, glycolysis/gluconeogenesis, systemic lupus erythematosus, bladder cancer, small cell lung cancer, complement and coagulation cascades, ECM-receptor interaction, axon guidance, and renal cell carcinoma. Interestingly, all 13 common core pathways were among the 32 common significant pathways (Fig. 2c). Details of the 13 common core pathways in each profile are shown in Table 1. The pathway in cancer (KEGG ID: 05200) was the top significant pathway, with a total degree centrality of 116. This suggests that the pathway in cancer may play an important role in PDAC development and progression.
Association of BIRC5, CKS2, ITGA3, ITGA6, and RALA with survival time in PDAC
To characterize the function of the pathway in cancer in PDAC, we analyzed the differentially expressed genes of this pathway in each of the six expression profiles. A total of 137 genes in this pathway were found to be differentially expressed in pancreatic cancer across all six datasets. Mining of The Cancer Genome Atlas (Internal ID: 811) and the International Cancer Genome Consortium database (Internal ID: 812) provided survival data connected with these 137 genes (Additional file 5: Table S5). We identified five genes that have been associated with patient survival time: baculoviral IAP repeat containing 5 (BIRC5), CDC28 protein kinase regulatory subunit 2 (CKS2), integrin subunit alpha 3 (ITGA3), integrin subunit alpha 6 (ITGA6), and RAS-like proto-oncogene A (RALA) (Fig. 3a). Patients with high expression of these five genes showed shorter survival. Consistent with these database results, we examined CKS2 expression in 102 patients with PDAC at our hospital and found that higher expression was associated with shorter survival (Fig. 3b).
Overexpression of CKS2 in PDAC
To determine which of these survival-related genes may be most important, we analyzed the fold-difference in the expression of each gene between PDAC and non-tumor tissues in the six expression profiles (Table 2). CKS2 showed the largest difference in expression, being up-regulated 4.82-fold in PDAC. Overexpression of CKS2 in PDAC compared with normal tissues has also been reported in the Oncomine database (Fig. 4), and data from the cbioportal database showed that CKS2 is amplified in many types of cancer (Fig. 5a). CKS2 expression levels positively correlate with copy number in 15 tumor cell lines (Fig. 5b).
CKS2 expression negatively correlates with DNA methylation
Tumor cells show dysregulation of DNA methylation levels and patterns, leading to significant changes in the overall level of DNA methylation in tumors . Data in The Cancer Genome Atlas for 149 PDAC tissue samples show that methylation levels at the CKS2 locus negatively correlate with CKS2 expression (r = − 0.368, P < 0.001; Fig. 5c).
CKS2 promotes PDAC cell proliferation
To characterize the potential functional role of CKS2, we analyzed genes whose expression parallels that of CKS2 based on RNA sequencing data from PDAC tissues in The Cancer Genome Atlas. A total of 325 genes co-expressed with CKS2 were identified (absolute r > 0.5; Additional file 6: Table S6), and these genes were enriched in the cell cycle pathway (Fig. 6a). Analysis of tumor and non-tumor tissues from our cohort of 102 patients with PDAC showed that CKS2 mRNA levels were higher in tumor tissues (Fig. 6b). CKS2 protein levels were higher in PDAC cell lines BxPC-3 and CFPAC-1 than in immortalized non-cancer pancreatic duct epithelial cell line HPDE6-C7 (Fig. 6c). Consistent with these results, CKS2 protein levels were higher in 20 samples of pancreatic cancer tissue than in the corresponding samples of adjacent normal tissue (Fig. 6d). Knockdown of CKS2 expression in the PDAC cell lines BxPC-3 and CFPAC-1 substantially suppressed proliferation (Fig. 6e, f).
In this study, we analyzed PDAC-associated genes and pathways in six frequently cited, published pancreatic cancer microarray datasets in the GEO database: GSE71989 [12, 13], GSE15471 , GSE16515 [15, 16], GSE32676 [17, 18], GSE41368 , and GSE28735 [18, 20]. We investigated differences in gene expression between pancreatic cancer tissues and normal tissues. Significant pathways enriched for differentially expressed genes were assembled into a pathway network, and core pathways showing high centrality degree in this network were chosen for further study. The pathway in cancer (KEGGID: 05200) was the top core pathway, and one gene in this pathway, CKS2, has been shown to be overexpressed in PDAC relative to normal tissue. CKS2 expression correlates positively with copy number and negatively with DNA methylation. CKS2 overexpression may be associated with shorter survival of patients with PDAC, probably due to CKS2 function in promoting proliferation of pancreatic cancer cells.
We used essentially all samples from these six GEO datasets to identify genes differentially expressed between PDAC tissue and normal tissue, although we excluded some samples to avoid the batch effect in previous studies . This also allowed us to pool the results from each dataset to perform integrative analysis. As part of this analysis, we created pathway networks , which were introduced to mine GEO data. Molecules in different signaling pathways interact with one another, so the pathways in a cell can be connected into a system. This is, to the best of our knowledge, the first analysis of PDAC-associated genes at the level of pathway networks. We obtained different pathway networks from the six GEO datasets, reflecting differences in samples, testing platform and other reasons (Additional file 7: Figure S1). If a pathway is important for pancreatic cancer development or progression, it should appear consistently across microarray experiments. We found 32 significant pathways common to all six GEO datasets, and within those, 13 core pathways, which are particularly likely to be important in the development and progression of pancreatic cancer. Of these core pathways, the pathway in cancer (KEGG ID: 05200) showed the highest degree centrality, suggesting that it is the most important pathway in pancreatic cancer. The pathway in cancer is an integrated signaling pathway that includes RAS-mediated pathways . KRAS is a well-known driver gene for pancreatic cancer, and mutation of this gene occurs in about 90% of patients . This suggests that our pathway network analysis can give reliable results.
We then analyzed whether all differentially expressed genes in the pathway in cancer across six GEO datasets are associated with PDAC survival, based on data in The Cancer Genome Atlas and International Cancer Genome Consortium databases. In both databases, expression of BIRC5, CKS2, ITGA3, ITGA6 and RALA were significantly associated with survival time. We did not identify any expression quantitative trait loci associated with differential expression of these five genes in pancreatic tissues, based on data in the GTEx database . This result may reflect the limited sample size. Of the five differentially expressed genes, CKS2 showed the greatest dysregulation in PDAC, which may be due to changes in copy number or/and DNA methylation (Fig. 5b, c).
CKS2, a gene encoding CDC28 protein kinase regulatory subunit, is located in chromosomal region 9q22.2. Previous studies have reported aberrant expression of this gene in cancers of the bladder , stomach , prostate , esophagus  and breast . Our database mining showed that many types of tumor overexpress CKS2 (Additional file 8: Figure S2 and Additional file 9: S3). Although CKS2 has been shown to play an important role in many types of cancer, its function in PDAC has not been reported. Our database mining showed that the genes associated with CKS2 are involved mainly in cell proliferation processes, such as the cell cycle and DNA replication. This implies that CKS2 knockdown may inhibit PDAC cell proliferation, which we confirmed directly in cultures of PDAC cell lines. These results are consistent with those obtained with other types of cancer.
BIRC5 is a member of the inhibitor of apoptosis gene family, which is overexpressed in many types of tumors [31, 32]. BIRC5 expression can be regulated by STAT3 , and it can be suppressed by curcumin, which induces apoptosis and ultimately inhibits PDAC cell proliferation . ITGA3 and ITGA6 are two members of the integrin alpha chain family of proteins. ITGA3 may promote proliferation or invasion in various types of cancer [35,36,37]. ITGA6 expression has been correlated with expression of migration-related genes , and it can promote the epithelial-mesenchymal transition  and tumor invasion . ITGA6 also plays a role in tumorigenesis . ITGA3 and ITGA6, both members of the integrin signaling pathway, are overexpressed in pancreatic cancer . RALA may participate in pancreatic cancer initiation . The details of how these genes participate in pancreatic cancer development and progression require further investigation.
Integrative analysis of six GEO microarray datasets of pancreatic cancer suggests that the pathway in cancer (KEGG ID: 05200) is an essential pathway in the development and progression of PDAC. High CKS2 expression is significantly associated with shorter survival time in patients with PDAC, which may reflect a function of CKS2 in promoting PDAC cell proliferation. These results demonstrate the potential of integrative data-mining from the GEO database in order to gain further insights into PDAC and other diseases.
cancer cell line encyclopedia
fetal bovine serum
gene expression omnibus
International Cancer Genome Consortium
pancreatic duct adenocarcinoma
short interfering RNA
The Cancer Genome Atlas
Hackeng WM, Hruban RH, Offerhaus GJ, Brosens LA. Surgical and molecular pathology of pancreatic neoplasms. Diagn Pathol. 2016;11:47.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2016. CA Cancer J Clin. 2016;66:7–30.
Chen W, Zheng R, Baade PD, Zhang S, Zeng H, Bray F, et al. Cancer statistics in China, 2015. CA Cancer J Clin. 2016;66:115–32.
Barrett T, Suzek TO, Troup DB, Wilhite SE, Ngau WC, Ledoux P, et al. NCBI GEO: mining millions of expression profiles–database and tools. Nucleic Acids Res. 2005;33:D562–6.
Yang S, Liu Y, Jiang N, Chen J, Leach L, Luo Z, et al. Genome-wide eQTLs and heritability for gene expression traits in unrelated individuals. BMC Genomics. 2014;15:13.
Stranger BE, Forrest MS, Dunning M, Ingle CE, Beazley C, Thorne N, et al. Relative impact of nucleotide and copy number variation on gene expression phenotypes. Science. 2007;315:848–53.
Gamazon ER, Stranger BE. The impact of human copy number variation on gene expression. Brief Funct Genomics. 2015;14:352–7.
Lee CJ, Evans J, Kim K, Chae H, Kim S. Determining the effect of DNA methylation on gene expression in cancer cells. Methods Mol Biol. 2014;1101:161–78.
da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.
Aguirre-Gamboa R, Gomez-Rueda H, Martinez-Ledesma E, Martinez-Torteya A, Chacolla-Huaringa R, Rodriguez-Barrientos A, et al. SurvExpress: an online biomarker validation tool and database for cancer gene expression data using survival analysis. PLoS ONE. 2013;8:e74250.
Wajed SA, Laird PW, DeMeester TR. DNA methylation: an alternative pathway to cancer. Ann Surg. 2001;234:10–20.
Giulietti M, Occhipinti G, Principato G, Piva F. Weighted gene co-expression network analysis reveals key genes involved in pancreatic ductal adenocarcinoma development. Cell Oncol. 2016;39:379–88.
Zhu T, Gao YF, Chen YX, Wang ZB, Yin JY, Mao XY, et al. Genome-scale analysis identifies GJB2 and ERO1LB as prognosis markers in patients with pancreatic cancer. Oncotarget. 2017;8:21281–9.
Badea L, Herlea V, Dima SO, Dumitrascu T, Popescu I. Combined gene expression analysis of whole-tissue and microdissected pancreatic ductal adenocarcinoma identifies genes specifically overexpressed in tumor epithelia. Hepatogastroenterology. 2008;55:2016–27.
Pei H, Li L, Fridley BL, Jenkins GD, Kalari KR, Lingle W, et al. FKBP51 affects cancer cell response to chemotherapy by negatively regulating Akt. Cancer Cell. 2009;16:259–66.
Yang D, Zhu Z, Wang W, Shen P, Wei Z, Wang C, et al. Expression profiles analysis of pancreatic cancer. Eur Rev Med Pharmacol Sci. 2013;17:311–7.
Donahue TR, Tran LM, Hill R, Li Y, Kovochich A, Calvopina JH, et al. Integrative survival-based molecular profiling of human pancreatic cancer. Clin Cancer Res. 2012;18:1352–63.
Li Z, Tian T, Lv F, Chang Y, Wang X, Zhang L, et al. Six1 promotes proliferation of pancreatic cancer cells via upregulation of cyclin D1 expression. PLoS ONE. 2013;8:e59203.
Frampton AE, Castellano L, Colombo T, Giovannetti E, Krell J, Jacob J, et al. MicroRNAs cooperatively inhibit a network of tumor suppressor genes to promote pancreatic tumor growth and progression. Gastroenterology. 2014;146(268–77):e18.
Zhang G, Schetter A, He P, Funamizu N, Gaedcke J, Ghadimi BM, et al. DPEP1 inhibits tumor cell invasiveness, enhances chemosensitivity and predicts clinical outcome in pancreatic ductal adenocarcinoma. PLoS ONE. 2012;7:e31507.
Chen C, Grennan K, Badner J, Zhang D, Gershon E, Jin L, et al. Removing batch effects in analysis of expression microarray data: an evaluation of six batch adjustment methods. PLoS ONE. 2011;6:e17238.
Dutta B, Wallqvist A, Reifman J. PathNet: a tool for pathway analysis using topological information. Source Code Biol Med. 2012;7:10.
Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2017;45:D353–61.
Ryan DP, Hong TS, Bardeesy N. Pancreatic adenocarcinoma. N Engl J Med. 2014;371:1039–49.
Consortium GT. Human genomics. The Genotype-Tissue Expression (GTEx) pilot analysis: multitissue gene regulation in humans. Science. 2015;348:648–60.
Chen R, Feng C, Xu Y. Cyclin-dependent kinase-associated protein Cks2 is associated with bladder cancer progression. J Int Med Res. 2011;39:533–40.
Kang MA, Kim JT, Kim JH, Kim SY, Kim YH, Yeom YI, et al. Upregulation of the cyclin kinase subunit CKS2 increases cell proliferation rate in gastric cancer. J Cancer Res Clin Oncol. 2009;135:761–9.
Stanbrough M, Bubley GJ, Ross K, Golub TR, Rubin MA, Penning TM, et al. Increased expression of genes converting adrenal androgens to testosterone in androgen-independent prostate cancer. Cancer Res. 2006;66:2815–25.
Kita Y, Nishizono Y, Okumura H, Uchikado Y, Sasaki K, Matsumoto M, et al. Clinical and biological impact of cyclin-dependent kinase subunit 2 in esophageal squamous cell carcinoma. Oncol Rep. 2014;31:1986–92.
van ‘t Veer LJ, Dai H, van de Vijver MJ, He YD, Hart AA, Mao M, et al. Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002;415:530–6.
Adida C, Berrebi D, Peuchmaur M, Reyes-Mugica M, Altieri DC. Anti-apoptosis gene, survivin, and prognosis of neuroblastoma. Lancet. 1998;351:882–3.
Tanaka K, Iwamoto S, Gon G, Nohara T, Iwamoto M, Tanigawa N. Expression of survivin and its relationship to loss of apoptosis in breast carcinomas. Clin Cancer Res. 2000;6:127–34.
Glienke W, Hausmann E, Bergmann L. Targeting STAT3 signaling in pancreatic cancer promotes antiapoptotic gene expression. Pancreas. 2011;40:323–4.
Glienke W, Maute L, Wicht J, Bergmann L. Curcumin inhibits constitutive STAT3 phosphorylation in human pancreatic cancer cell lines and downregulation of survivin/BIRC5 gene expression. Cancer Invest. 2010;28:166–71.
Akiyama SK, Olden K, Yamada KM. Fibronectin and integrins in invasion and metastasis. Cancer Metastasis Rev. 1995;14:173–89.
Koshizuka K, Hanazawa T, Kikkawa N, Arai T, Okato A, Kurozumi A, et al. Regulation of ITGA3 by the anti-tumor miR-199 family inhibits cancer cell migration and invasion in head and neck cancer. Cancer Sci. 2017;108:1681–92.
Sakaguchi T, Yoshino H, Yonemori M, Miyamoto K, Sugita S, Matsushita R, et al. Regulation of ITGA3 by the dual-stranded microRNA-199 family as a potential prognostic marker in bladder cancer. Br J Cancer. 2017;116:1077–87.
Golbert DC, Correa-de-Santana E, Ribeiro-Alves M, de Vasconcelos AT, Savino W. ITGA6 gene silencing by RNA interference modulates the expression of a large number of cell migration-related genes in human thymic epithelial cells. BMC Genomics. 2013;14(Suppl 6):S3.
Zhang K, Myllymaki SM, Gao P, Devarajan R, Kytola V, Nykter M, et al. Oncogenic K-Ras upregulates ITGA6 expression via FOSL1 to induce anoikis resistance and synergizes with alphaV-Class integrins to promote EMT. Oncogene. 2017;36:5681–94.
Brooks DL, Schwab LP, Krutilina R, Parke DN, Sethuraman A, Hoogewijs D, et al. ITGA6 is directly regulated by hypoxia-inducible factors and enriches for cancer stem cell activity and invasion in metastatic breast cancer models. Mol Cancer. 2016;15:26.
Kwon J, Lee TS, Lee HW, Kang MC, Yoon HJ, Kim JH, et al. Integrin alpha 6: a novel therapeutic target in esophageal squamous cell carcinoma. Int J Oncol. 2013;43:1523–30.
Van den Broeck A, Vankelecom H, Van Eijsden R, Govaere O, Topal B. Molecular markers associated with outcome and metastasis in human pancreatic cancer. J Exp Clin Cancer Res. 2012;31:68.
Lim KH, O’Hayer K, Adam SJ, Kendall SD, Campbell PM, Der CJ, et al. Divergent roles for RalA and RalB in malignant growth of human pancreatic carcinoma cells. Curr Biol. 2006;16:2385–94.
DL conceptualized and supervised this study. JL contributed to the study design; JL, WenleT and LP conducted data analysis. JL and QC performed the proliferation assay. JianZ, JialiangZ and XH performed sample preparation and association analysis. JL, CW and WenT were involved in manuscript preparation. All authors read and approved the final manuscript.
We thank Prof. Rufu Chen from Department of Hepatobiliary Surgery, Sun Yat-sen Memorial Hospital to provide PDAC and adjacent normal samples for this study.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated or analyzed during this study are included in this published article.
Consent for publication
We all consent for publication.
Ethics approval and consent to participate
Informed consent was obtained from all participants, and this study was approved by the Department of Pancreaticobiliary, Sun Yat-sen Memorial Hospital and Sun Yat-sen University Cancer Center. All experiments on the participants in this study were performed in accordance with the relevant guidelines and regulations.
This research was supported by the National Natural Science Foundation of China (No. 81572793, 81772586 and 81602461), a Grant (No. 2016A030313283) from the Natural Science Foundation of Guangdong Province and the Peking Union Medical College Doctor Innovation Fund (No. 2015-0710-17).
Characteristics of patients with pancreatic cancer in this study.
Short interfering RNA (siRNA) sequences used in this study.
Primers used for real-time PCR in this study.
Pancreatic cancer microarray datasets included in the study.
Hazard ratio (HR) of death for genes in the pathway in cancer that are differentially expressed in pancreatic cancer, based on data in two databases.
Genes associated with CKS2 expression level in pancreatic cancer tissues, based on data in The Cancer Genome Atlas.
Pathway networks of all GEO datasets. Pathway networks for the datasets GSE15471 (A), GSE16515 (B), GSE32676 (C), GSE71989 (D), GSE28735 (E) and GSE42368 (F) were drawn using the GCBI online tool as described in Methods. Yellow dots indicate involvement of up- and down-regulated signaling pathway genes; red dots, up-regulated signaling pathway genes; and blue dots, down-regulated signaling pathway genes. The arrow points from the upstream toward the downstream signaling pathway.
CKS2 expression in different types of cancer cell lines, based on data from the CCLE database (https://portals.broadinstitute.org/ccle/home).
CKS2 expression across multiple types of cancer and corresponding normal tissues. The numbers of studies reporting up- or down-regulation is indicated, respectively, in red or blue boxes. Color intensity reflects the best gene rank percentile for analyses within the box. The following settings were used for the analysis: P=0.05, fold change=all, gene rank=all, and data type=mRNA (http://www.oncomine.org).
About this article
Cite this article
Li, J., Tan, W., Peng, L. et al. Integrative analysis of gene expression profiles reveals specific signaling pathways associated with pancreatic duct adenocarcinoma. Cancer Commun 38, 13 (2018). https://doi.org/10.1186/s40880-018-0289-9
- Pancreatic cancer
- Pathway network