- Open Access
Identification of hub genes with diagnostic values in pancreatic cancer by bioinformatics analyses and supervised learning methods
© The Author(s). 2018
- Received: 7 August 2018
- Accepted: 24 October 2018
- Published: 14 November 2018
Pancreatic cancer is one of the most lethal tumors with poor prognosis, and lacks of effective biomarkers in diagnosis and treatment. The aim of this investigation was to identify hub genes in pancreatic cancer, which would serve as potential biomarkers for cancer diagnosis and therapy in the future.
Combination of two expression profiles of GSE16515 and GSE22780 from Gene Expression Omnibus (GEO) database was served as training set. Differentially expressed genes (DEGs) with top 25% variance followed by protein-protein interaction (PPI) network were performed to find candidate genes. Then, hub genes were further screened by survival and cox analyses in The Cancer Genome Atlas (TCGA) database. Finally, hub genes were validated in GSE15471 dataset from GEO by supervised learning methods k-nearest neighbor (kNN) and random forest algorithms.
After quality control and batch effect elimination of training set, 181 DEGs bearing top 25% variance were identified as candidate genes. Then, two hub genes, MMP7 and ITGA2, correlating with diagnosis and prognosis of pancreatic cancer were screened as hub genes according to above-mentioned bioinformatics methods. Finally, hub genes were demonstrated to successfully differ tumor samples from normal tissues with predictive accuracies reached to 93.59 and 81.31% by using kNN and random forest algorithms, respectively.
All the hub genes were associated with the regulation of tumor microenvironment, which implicated in tumor proliferation, progression, migration, and metastasis. Our results provide a novel prospect for diagnosis and treatment of pancreatic cancer, which may have a further application in clinical.
- Pancreatic cancer
- Bioinformatics analysis
- Differentially expressed genes
- Hub genes
Pancreatic cancer is one of the most lethal tumors due to the poor prognosis, and now it is the fourth or fifth most common causes of cancer mortality in developed countries . And it is estimated that by the year 2020, pancreatic cancer would move to the second leading cause of death . Although some advances in understanding the molecular mechanisms of pancreatic cancer have been achieved, there still exist difficulties in early diagnosis due to non-specific symptoms and lacking effective testing identification, making it usually found in its late stage . Until now, 1-year survival in pancreatic cancer patients is still not significantly improved , and the 5-year survival is less than 10% .
Numerous studies have focused on the investigation of biomarkers and molecular mechanisms of pancreatic cancers, and it is demonstrated that accumulated mutations in genes like oncogene Kras, and tumor-suppressor genes including P16 as well as TP53 resulted in the occurrence of pancreatic cancer . One study performed the whole-genome sequencing and copy number variation (CNV) analyses showed that several genes including TP53, SMAD4, CDKN2A, ARID1A, ROBO2, PREX2, and KDM6A were disrupt resulting from chromosomal rearrangements in pancreatic ductal adenocarcinomas patients . Molecular mechanisms researches demonstrated that overexpression of protein-coupled receptor GPR87 enhanced pancreatic cancer aggressiveness by activating NF-κB signaling pathway . Moreover, Zhong and colleges have found that functional P38 MAPK activity contributed to overall survival through suppressing JNK signaling in pancreatic cancer . In addition, aberrant expressions of some microRNAs have emerged as an important hallmark of cancer recently . It was reported that microRNA-21 was overexpressed in pancreatic cancer, and could serve as a potential predictor of survival . One study has found that miR-506 facilitated pancreatic cancer progression and chemoresistance via SPHK1/Akt/NF-κB signaling pathway . Another study demonstrated that suppressing microRNA-34 expression downregulated Bcl-2 and Notch1/2 in pancreatic cancer cells, as well as significantly inhibited cell growth and invasion, induced apoptosis and G1 and G2/M arrest in cell cycle, and sensitized the cells to chemotherapy and radiation .
However, traditional experimental methods as mentioned above could only identify single gene or a few genes at once, which limits large-scale investigation of hub genes and pathways in the systematic biology level. Development of microarray and sequencing technologies provides better methods for biomarker screening and molecular mechanism discovery in cancer research. Recent years with the accessibility of multi-omics database like Gene Expression Omnibus (GEO)  as well as The Cancer Genome Atlas (TCGA)  and so on, it is now possible to acquire multi-sample data and compare cancer profiles with normal profiles in multiple omics dimensions. On one hand, omics data in multiple dimensions leading to the system biology- and/or network-based approach, which could better understand the dysregulated molecular mechanisms in cancer development and progression . On the other hand, biology- and/or network-based method can not only identify critical genes but also can detect corresponding pathways and/or interactive network, which may provide better insights into molecular mechanisms investigation than dysregulated gene analysis individually . For example, Kras was proved to be the most frequently mutated gene in pancreatic ductal adenocarcinoma , and the mutation of Kras was a hallmark of pancreatic cancer . However, inhibitors targeting Kras gene were largely unsuccessful, while some omics-based strategies targeting Kras correlated pathways and interactive genes were proved to bear better therapeutic effects than targeting Kras individually .
To date, diagnosis of pancreatic cancer is mainly based on clinical signs and pathology confirmation. However, the specific symptoms and pathological imagines may only be detected unambiguously at the late stage of pancreatic cancer, which may lead to a limited therapies and poor prognosis. This raises an urgent need for the development of reliable biomarkers which can effectively differ tumor from normal tissues based on analyses of gene expression profiles. Herein, in order to identify novel diagnostic predictors and molecular markers, we integrated two microarray datasets from GEO database, and 11 candidate genes significantly differentially expressed between tumor and normal samples were screened by bioinformatics analyses. Then two hub genes, matrix metallopeptidase 7 (MMP7) and integrin, alpha 2 (ITGA2), were further identified by survival and cox analyses in TCGA database. These two hub genes were validated in another expression profile from GEO database, demonstrating that these hub genes can successfully differ normal tissues from tumor samples. The predictive accuracies of k-nearest neighbor (kNN) and random forest algorithms were almost 94% and almost 82%, respectively. Results in our study may provide an auxiliary evidence of pancreatic cancer diagnosis and therapy in the future.
Data collection and preprocessing
All the datasets were performed by Affymetrix Human Genome U133 Plus 2.0 Array (Affymetrix, Santa Clara, CA, USA). GSE16515 dataset included 36 malignant pancreatic samples and 16 normal pancreatic samples, while the corresponding numbers in GSE22780 dataset were 8 and 8. In order to obtain sample balance, combination of GSE16515 and GSE22780 was used as training set to determine hub genes. Besides, raw expression data of GSE15471 was downloaded from GEO, also performed by Affymetrix Human Genome U133 Plus 2.0 Array. It composed of 39 normal and 39 malignant pancreatic samples, and served as testing set.
Firstly, the quality of all the datasets were detected with “affyPLM” package in R, herein FitPLM weight, residual, relative log expression (RLE), normalized unscaled standard errors (NUSE), and RNA degradation images were evaluated. Then robust multiarray averaging (RMA) with “affy” package was used to do the background correction and normalization. Before subsequent hub gene selection in training set, empirical Bayes framework with “sva” package in R was used to adjust the batch effects between these two datasets.
In addition, we also downloaded RNA-sequencing data of pancreatic cancer from The Cancer Genome Atlas (TCGA) database (https://cancergenome.nih.gov/), and all the raw data were also converted into gene symbol expression matrix by R software and Perl software.
Differentially expressed genes screening
Herein, “limma” package was used to detect differentially expressed genes (DEGs) between malignant pancreatic samples and normal samples in training set with the threshold of adj.P value < 0.01 and absolute log2-based fold change > 1.
Candidate gene selection
Variance of every DEGs in different samples were calculated and sorted by descending order, and the top 25% results were selected. Then, 181 genes bearing top 25% variance were uploaded in Search Tool for the Retrieval of Interacting Genes (STRING) database (https://string-db.org/), and PPI network was constructed  by setting minimum required interaction score at 0.700. Then a plug-in Cytohubba in Cytoscape  was used to further screen candidate genes. Herein, degree algorithm was applied and the screening criterion was degree > 5.
Hub gene screening by survival and cox regression analyses in TCGA
Candidate genes were further screened by survival analysis and cox regression analysis in TCGA database with “survival” package. Genes with P value less than 0.05 both in survival analysis and cox analysis were further screened as hub genes.
Gene ontology annotation and pathway analyses of candidate genes
In order to depict the biological function of candidate genes, gene ontology (GO) biological process enrichments were performed through Database for Annotation, Visualization and Integrated Discovery (DAVID) (https://david.ncifcrf.gov/) [24, 25]. And the visualization of GO results was performed by “GOplot” package in R.
Validation of hub genes by supervised learning methods
In order to verify whether these hub genes were “real hub genes” to discriminate tumor and normal samples, kNN algorithm in “class” package and random forest algorithm in “randomForest” package were performed. The accuracy was used to evaluate the predictive results. Herein, random forest algorithm was rerun for 100 times, and the mean value of the accuracies was calculated finally.
Identification of DEGs
Determination of candidate genes
Selection of hub genes by survival and cox analyses
Functional annotation and pathway enrichment
Functional annotation of two hub genes ITGA2 and MMP7
Positive regulation of smooth muscle contraction
Response to l-ascorbic acid
Positive regulation of cell projection organization
Mesodermal cell differentiation
Extracellular matrix organization
Response to drug
Collagen catabolic process
Extracellular matrix disassembly
Prediction of pancreatic cancer by hub genes
Comparison of predictive accuracy resulted from different screening methods
Minimum required interaction score
Accuracy of kNN algorithm
Mean accuracy of random forest algorithm (rerun 100 times)
Method 1: 724 DGEs-181 candidate genes-genes bearing top 10 degrees in PPI-2 hub genes by survival analysis and cox analysis
Method 2: 724 DGEs-181 candidate genes-genes bearing top 10 degrees in PPI
ALB, EGF, FN1, ITGA2, COL1A2, SPARC, COL3A1, TIMP1, COL5A1, COL11A1, MMP7
Method 3: 724 DGEs-genes bearing top 10 degrees in PPI-2 hub genes by survival analysis and cox analysis
Method 4: 724 DGEs-genes bearing top 10 degrees in PPI
CCNB1, CCNA2, MAD2L1, TOP2A, UBE2C, CDC20, TTK, MELK, BUB1B, NDC80
Method 5: 724 DGEs-181 candidate genes-genes bearing top 10 degrees in PPI-1 hub genes by survival analysis and cox analysis
Method 6: 724 DGEs-181 candidate genes-genes bearing top 10 degrees in PPI
ALB, EGF, ITGA2, FN1, COL1A2, TIMP1, MMP1, COL3A1, PTGS2, CEL
Method7:724 DGEs-genes bearing top10 degrees in PPI-2 hub genes by survival analysis and cox analysis
Method8:724 DGEs-genes bearing top 10 degrees in PPI
ALB, GAPDH, EGF, TOP2A, CCNB1, NDC80, CCNA2, CDC20, UBE2C, BUB1B, MAD2L1, TTK, OIP5, KIF11
Compared with other cancers, the occurrence of pancreatic cancer is relatively rare; however, it is still a lethal disease with poor prognosis. Until now, there still lacks effective therapies against pancreatic cancer, and many novel therapies are in the experimental stage. Therefore, it is important to find some potential hub genes playing crucial roles in regulating cancer occurrence and progression, which may become key targets in the treatment of pancreatic cancer in the future. In addition, these hub genes effectively differing cancer tissues from normal samples may provide novel auxiliary evidence in pancreatic cancer diagnosis. It is demonstrated that pancreatic cancer results from the accumulation of acquired mutations, which may lead to the upregulation of some oncogenes and downregulation of some tumor-suppressing genes and genomic maintenance genes . Therefore, there might exist some DEGs between normal and tumor samples, and these DEGs may play important roles in regulating tumor occurrence, development, and progression. In the present study, two genes ITGA2 and MMP7 were screened from DEGs as hub genes by using a series of bioinformatics methods, and they could discriminate normal samples and tumor samples.
The matrix metalloproteinase (MMPs) is a family of enzymes, bearing the capability to cleave extracellular matrix substrates , as well as promotes the release of pro-TNF-α, Fas ligand, and some cytokines in various cancers cells . One previous study has experimentally demonstrated that genes in matrix metallopeptidase family, collagen family, and integrin family were upregulated in pancreatic cancer, and they may correlate with cancer activity and poor prognosis . MMPs also involved in proliferative, migrating, and differentiated processes in cells . The interaction between MMPs and extracellular ligand induced a series of signaling cascade, and thus led to the functional regulation of intracellular and extracellular activities. The expression of MMP7 has been reported to be upregulated in several kinds of cancer, including colon cancer , pancreatic cancer , breast cancer , gastric cancer , and esophageal cancer . One study has demonstrated that multiplex detection of pancreatic biomarkers CA19-9, MMP7, and MUC4 in sera samples were of high sensitivity, which may act as the critical biomarker in diagnosis of pancreatic cancer . Another study compared tumor tissues with healthy control samples revealed that MMP7 was highly predictive for advanced stage of pancreatic cancer, which strongly associated with N1 status, T3/T4 stage, moderate/poor differentiation, and perineural invasion . It has been reported that Stat3 was a critical factor to facilitate precursor formation and enforced MMP7 expression in pancreatic cancer cells, while MMP7 level was correlated with metastasis and survival in pancreatic cancer patients .
ITGA2 encoding by ITGA2 gene is the alpha subunit of the transmembrane receptor integrin, and it mainly exerts the adhesive roles in cell-cell interaction, also promotes the generation and adhesion of newly synthesized extracellular matrix [37, 38]. The polymorphisms of ITGA2 gene was related to the poor survival of nasopharyngeal carcinoma . ITGA2 gene was reported to play migrating roles in colon cancer cells , and it expressed in colorectal cancer with liver metastasis tissues but absent in normal tissue . In addition, epigenetic modifications such as DNA methylation were also important in tumorigenesis, and hypomethylation of ITGA2 with high gene expression was associated with poor survival in pancreatic cancer patients . One research has found that ITGA2 was overexpressed in a variety of gastric cancer patients mainly playing pro-survival roles, and the blockage of ITGA2 could induce apoptosis and inhibit cell migration in gastric cancer . Another research in gastric cancer revealed that HMGA2, FOXL2, and ITGA2 were increased in metastatic lymph nodes and distant metastases in gastric cancer, and suppressing the HMGA2-FOXL2-ITGA2 pathway could serve as a new strategy in further treatment in gastric cancer . The transcriptional co-activators yes-associated protein (YAP) was considered as oncogene in many types of cancer; ITGA2 stimulating YAP activity was associated with unfavorable survival of pancreatic cancer patients .
In order to validate whether these genes were real hub genes, another mRNA expression profile GSE15471 from GEO database was utilized as testing set. Herein, kNN and random forest algorithms were performed to detect whether these hub genes could successfully distinguish tumor tissues from normal samples. We can see from Table 2 that hub genes selected by method 1 in this study represented the highest accuracy reaching to 94% approximately with 2.56% false negative and 3.84% false positive. In the cases of differing from tumor and normal samples, reduction of false negative results was more important than the reduction of false positive result. Since false negative results may lead to wrongly diagnose pancreatic cancer as normal condition, it may result in the delay of timely treatment, and further lead to more serious progression of disease as well as more waste of medical resources and costs. Bedsides, random forest algorithm also represented highly predictive accuracy of 81.31% after rerun for 100 times of method 1. Therefore, method 1 bear highly predictive accuracies in both of the two methods, and it could be inferred that these two hub genes were real hub genes, which could successfully discriminate normal and tumor samples. Another interesting result could be found in Table 2 that selection of genes with top 25% variance obviously increased the predictive accuracy from 70 to 94% (method 1 vs. method 3).
In addition, we can choose different minimum required interaction score when constructing PPI network. Minimum required interaction score is a threshold providing a score for each interactive pair, which is computed as the joint probability from different evidence (e.g., protein interaction, fusion, co-expression, text mining). Higher score may represent more confident interaction while lower score may lead to more false positives . In order to elucidate whether setting different minimum required interaction score may have influence on hub gene selection, predictive accuracy was compared (shown in Table 2). It can be found that higher minimum required interaction score led to much higher predictive accuracies; method 5 bear the highest accuracy of 85% while the predictive accuracy of method 1 could reach to almost 95% by using kNN method.
Hub genes screened in this study were rational. Firstly, all the candidate genes and these two hub genes were closely correlated with the progression of tumor. As shown in GO enrichment, most of the candidate genes were implicated in the biological process of extracellular matrix, cell adhesion, cell proliferation, and signal transduction; they play important role in the progression of cancers. Moreover, both of the hub genes were implicated in the regulation of tumor microenvironment, including the regulation of tumor cells, stroma cells, extracellular matrix (ECM), and some extracellular molecules like cytokines as well as chemokines. It has been demonstrated that microenvironment was usually dysregulated and disorganized in cancer cells. Thus, disordered microenvironment may be favorable to tumor proliferation, progression, invasion and metastasis, and exert drug-hampering roles [46, 47], and now some treatment strategies have focused on the regulation of tumor microenvironment. Since pancreatic cancer was featured as uncontrolled and malignant invasion and migration, therefore we can infer that these hub genes implicated in tumor microenvironment might be core meditators in pancreatic cancer diagnosis and therapy. Secondly, two supervised learning methods were performed, and both of the predictive results of these two hub genes were good with lower false negative in discriminating tumor samples from normal samples.
However, there also exist some limitations in our study. Firstly, the number of samples in our study is not too much. According to the dataset screening criteria, three datasets were included in our study. There were 146 samples totally, of which 68 were training set and 78 were testing set. In the future, with more and more investigations about pancreatic cancer would be performed, more samples should be included. Secondly, in this study, we mainly focused on the genes in the pancreatic tissue not the genes from circulating tumor cells (CTC) nor circulating tumor DNA (ctDNA) in peripheral blood, since the genes in tissue are more accurate to analyze the important biomarkers. Moreover, the datasets about peripheral blood in GEO database are not enough to do the same research. In the future, the microarray analysis of DNA in peripheral blood of pancreatic cancer patients should be further proposed. Thirdly, in our study, all the hub genes were screened and validated only by bioinformatics method, and further exploration of the biological functions and molecular mechanisms of these hub genes both in vitro and in vivo are needed to be fulfilled.
In summary, we conducted a series of bioinformatics methods to find DEGs, further screened and validated hub genes. These two hub genes, ITGA2 and MMP7, may act as potential diagnostic and therapeutic biomarkers in pancreatic cancer patients. This study provides several useful hub genes for future in vitro and in vivo investigations of their molecular mechanisms in pancreatic cancer diagnosis and therapy. And profile data mining by bioinformatics analysis is an available method to find potential diagnostic and prognostic biomarkers systematically. Nevertheless, further molecular mechanisms investigations by biological experiments are still needed to be verified in pancreatic cancer cells.
We sincerely thank Professor Qiu Li and Doctor Pengfei Zhang (Department of Cancer Center, West China Hospital) for their good suggestions on pancreatic cancer diagnosis, treatment, and management.
This work was supported by Science & Technology Department of Sichuan Province funding project (No. 2016FZ0108), Health and Family Planning Commission of Sichuan Province project (No. 17PJ443), international cooperation project of Science and Technology Department of Sichuan Province (2016HH0069).
Availability of data and materials
The datasets used and/or analyzed during the current study are downloaded from the GEO database (http://www.ncbi.nlm.nih.gov/geo/), and the datasets used in this study were GSE16515, GSE22780, and GSE15471.
CL participated in study design, data preprocessing and analysis, as well as manuscript writing. XZ involved in data preprocessing, data analysis, and results double check. HY participated in supervised learning analysis, and manuscript writing. YG implicated in statistics analysis and language polishing. WZ designed the study and help to write the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Li C, Heidt DG, Dalerba P, Burant CF, Zhang L, Adsay V, Wicha M, Clarke MF, Simeone DM. Identification of pancreatic cancer stem cells. Cancer Res. 2007;67:1030–7.View ArticleGoogle Scholar
- Rahib L, Smith BD, Aizenberg R, Rosenzweig AB, Fleshman JM, Matrisian LM. Projecting cancer incidence and deaths to 2030: the unexpected burden of thyroid, liver, and pancreas cancers in the United States. Cancer Res. 2014;74:2913–21.View ArticleGoogle Scholar
- Korc M. Pancreatic cancer-associated stroma production. Am J Surg. 2007;194:S84–6.View ArticleGoogle Scholar
- Kleeff J, Korc M, Apte M, La Vecchia C, Johnson CD, Biankin AV, Neale RE, Tempero M, Tuveson DA, Hruban RH, Neoptolemos JP. Pancreatic cancer. Nat Rev Dis Primers. 2016;2:16022.View ArticleGoogle Scholar
- Neoptolemos JP, Palmer DH, Ghaneh P, Psarelli EE, Valle JW, Halloran CM, Faluyi O, O'Reilly DA, Cunningham D, Wadsley J, et al. Comparison of adjuvant gemcitabine and capecitabine with gemcitabine monotherapy in patients with resected pancreatic cancer (ESPAC-4): a multicentre, open-label, randomised, phase 3 trial. Lancet. 2017;389:1011–24.View ArticleGoogle Scholar
- Waddell N, Pajic M, Patch AM, Chang DK, Kassahn KS, Bailey P, Johns AL, Miller D, Nones K, Quek K, et al. Whole genomes redefine the mutational landscape of pancreatic cancer. Nature. 2015;518:495–501.View ArticleGoogle Scholar
- Wang L, Zhou W, Zhong Y, Huo Y, Fan P, Zhan S, Xiao J, Jin X, Gou S, Yin T, et al. Overexpression of G protein-coupled receptor GPR87 promotes pancreatic cancer aggressiveness and activates NF-kappaB signaling pathway. Mol Cancer. 2017;16:61.View ArticleGoogle Scholar
- Zhong Y, Naito Y, Cope L, Naranjo-Suarez S, Saunders T, Hong SM, Goggins MG, Herman JM, Wolfgang CL, Iacobuzio-Donahue CA. Functional p38 MAPK identified by biomarker profiling of pancreatic cancer restrains growth through JNK inhibition and correlates with improved survival. Clin Cancer Res. 2014;20:6200–11.View ArticleGoogle Scholar
- Khan MA, Zubair H, Srivastava SK, Singh S, Singh AP. Insights into the role of microRNAs in pancreatic cancer pathogenesis: potential for diagnosis, prognosis, and therapy. Adv Exp Med Biol. 2015;889:71–87.View ArticleGoogle Scholar
- Dillhoff M, Liu J, Frankel W, Croce C, Bloomston M. MicroRNA-21 is overexpressed in pancreatic cancer and a potential predictor of survival. J Gastrointest Surg. 2008;12:2171–6.View ArticleGoogle Scholar
- Li J, Wu H, Li W, Yin L, Guo S, Xu X, Ouyang Y, Zhao Z, Liu S, Tian Y, et al. Downregulated miR-506 expression facilitates pancreatic cancer progression and chemoresistance via SPHK1/Akt/NF-kappaB signaling. Oncogene. 2016;35:5501–14.View ArticleGoogle Scholar
- Ji Q, Hao XB, Zhang M, Tang WH, Meng Y, Li L, Xiang DB, DeSano JT, Bommer GT, Fan DM, et al. MicroRNA miR-34 inhibits human pancreatic cancer tumor-initiating cells. PLoS One. 2009;4:e6816.View ArticleGoogle Scholar
- Edgar R, Domrachev M, Lash AE. Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30:207–10.View ArticleGoogle Scholar
- Cancer Genome Atlas Research N, Weinstein JN, Collisson EA, Mills GB, Shaw KR, Ozenberger BA, Ellrott K, Shmulevich I, Sander C, Stuart JM. The Cancer Genome Atlas Pan-Cancer analysis project. Nat Genet. 2013;45:1113–20.View ArticleGoogle Scholar
- Rajamani D, Bhasin MK. Identification of key regulators of pancreatic cancer progression through multidimensional systems-level analysis. Genome Medicine. 2016;8:38.View ArticleGoogle Scholar
- Chuang HY, Lee E, Liu YT, Lee D, Ideker T. Network-based classification of breast cancer metastasis. Mol Syst Biol. 2007;3:140.View ArticleGoogle Scholar
- Sivakumar S, de Santiago I, Chlon L, Markowetz F. Master regulators of oncogenic KRAS response in pancreatic cancer: an integrative network biology analysis. PLoS Medicine / Public Library of Science. 2017;14:e1002223.Google Scholar
- Muzumdar MD, Chen PY, Dorans KJ, Chung KM, Bhutkar A, Hong E, Noll EM, Sprick MR, Trumpp A, Jacks T. Survival of pancreatic cancer cells lacking KRAS function. Nat Commun. 2017;8:1090.View ArticleGoogle Scholar
- Wolfgang CL, Herman JM, Laheru DA, Klein AP, Erdek MA, Fishman EK, Hruban RH. Recent progress in pancreatic cancer. CA Cancer J Clin. 2013;63:318–48.View ArticleGoogle Scholar
- Pei H, Li L, Fridley BL, Jenkins GD, Kalari KR, Lingle W, Petersen G, Lou Z, Wang L. FKBP51 affects cancer cell response to chemotherapy by negatively regulating Akt. Cancer Cell. 2009;16:259–66.View ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, Simonovic M, Roth A, Santos A, Tsafou KP, et al. STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015;43:D447–52.View ArticleGoogle Scholar
- Chin CH, Chen SH, Wu HH, Ho CW, Ko MT, Lin CY. cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol. 2014;8(Suppl 4):S11.View ArticleGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al. Gene ontology: tool for the unification of biology. The gene ontology consortium. Nat Genet. 2000;25:25–9.View ArticleGoogle Scholar
- Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA. DAVID: database for annotation, visualization, and integrated discovery. Genome Biol. 2003;4:P3.View ArticleGoogle Scholar
- Van Wart HE, Birkedal-Hansen H. The cysteine switch: a principle of regulation of metalloproteinase activity with potential applicability to the entire matrix metalloproteinase gene family. Proc Natl Acad Sci U S A. 1990;87:5578–82.View ArticleGoogle Scholar
- Jang B, Jung H, Choi S, Lee YH, Lee S-T. Oh E-S: Syndecan-2 cytoplasmic domain up-regulates matrix metalloproteinase-7 expression via the protein kinase Cgamma-mediated FAK/ERK signaling pathway in colon cancer. J Biol Chem. 2017;292:16321–32.View ArticleGoogle Scholar
- Moffitt RA, Marayati R, Flate EL, Volmar KE, Loeza SG, Hoadley KA, Rashid NU, Williams LA, Eaton SC, Chung AH, et al. Virtual microdissection identifies distinct tumor- and stroma-specific subtypes of pancreatic ductal adenocarcinoma. Nat Genet. 2015;47:1168–78.View ArticleGoogle Scholar
- Malemud CJ. Matrix metalloproteinases (MMPs) in health and disease: an overview. Front Biosci. 2006;11:1696–701.View ArticleGoogle Scholar
- Chen SH, Hung WC, Wang P, Paul C, Konstantopoulos K. Mesothelin binding to CA125/MUC16 promotes pancreatic cancer cell motility and invasion via MMP-7 activation. Sci Rep. 2013;3:1870.View ArticleGoogle Scholar
- Lin HY, Sun SM, Lu XF, Chen PY, Chen CF, Liang WQ, Peng CY. CCR10 activation stimulates the invasion and migration of breast cancer cells through the ERK1/2/MMP-7 signaling pathway. Int Immunopharmacol. 2017;51:124–30.View ArticleGoogle Scholar
- Xu J, E C, Yao Y, Ren S, Wang G, Jin H. Matrix metalloproteinase expression and molecular interaction network analysis in gastric cancer. Oncol Lett. 2016;12:2403–8.View ArticleGoogle Scholar
- Juchniewicz A, Kowalczuk O, Milewski R, Laudanski W, Dziegielewski P, Kozlowski M, Niklinski J. MMP-10, MMP-7, TIMP-1 and TIMP-2 mRNA expression in esophageal cancer. Acta Biochim Pol. 2017;64:295–9.View ArticleGoogle Scholar
- Banaei N, Foley A, Houghton JM, Sun YB, Kim B. Multiplex detection of pancreatic cancer biomarkers using a SERS-based immunoassay. Nanotechnology. 2017;28:455101.View ArticleGoogle Scholar
- Wang SC, Parekh JR, Porembka MR, Nathan H, D'Angelica MI, DeMatteo RP, Fong Y, Kingham TP, Jarnagin WR, Allen PJ. A pilot study evaluating serum MMP7 as a preoperative prognostic marker for pancreatic ductal adenocarcinoma patients. J Gastrointest Surg. 2016;20:899–904.View ArticleGoogle Scholar
- Fukuda A, Wang SC, JPt M, Folias AE, Liou A, Kim GE, Akira S, Boucher KM, Firpo MA, Mulvihill SJ, Hebrok M. Stat3 and MMP7 contribute to pancreatic ductal adenocarcinoma initiation and progression. Cancer Cell. 2011;19:441–55.View ArticleGoogle Scholar
- Bergelson JM, St John N, Kawaguchi S, Chan M, Stubdal H, Modlin J, Finberg RW. Infection by echoviruses 1 and 8 depends on the alpha 2 subunit of human VLA-2. J Virol. 1993;67:6847–52.PubMedPubMed CentralGoogle Scholar
- Graham KL, Halasz P, Tan Y, Hewish MJ, Takada Y, Mackow ER, Robinson MK, Coulson BS. Integrin-using rotaviruses bind alpha2beta1 integrin alpha2 I domain via VP4 DGE sequence and recognize alphaXbeta2 and alphaVbeta3 by using VP7 during cell entry. J Virol. 2003;77:9969–78.View ArticleGoogle Scholar
- Ban EZ, Lye MS, Chong PP, Yap YY, Lim SYC, Abdul Rahman H. Association of hOGG1 Ser326Cys, ITGA2 C807T, TNF-A -308G>a and XPD Lys751Gln polymorphisms with the survival of Malaysian NPC patients. PLoS One. 2018;13:e0198332.View ArticleGoogle Scholar
- Ferraro A, Boni T, Pintzas A. EZH2 regulates cofilin activity and colon cancer cell migration by targeting ITGA2 gene. PLoS ONE [electronic Resource]. 2014;9:e115276.View ArticleGoogle Scholar
- Yang Q, Bavi P, Wang JY, Roehrl MH. Immuno-proteomic discovery of tumor tissue autoantigens identifies olfactomedin 4, CD11b, and integrin alpha-2 as markers of colorectal cancer with liver metastases. J Proteome. 2017;168:53–65.View ArticleGoogle Scholar
- Nones K, Waddell N, Song S, Patch AM, Miller D, Johns A, Wu J, Kassahn KS, Wood D, Bailey P, et al. Genome-wide DNA methylation patterns in pancreatic ductal adenocarcinoma reveal epigenetic deregulation of SLIT-ROBO, ITGA2 and MET signaling. Int J Cancer. 2014;135:1110–8.View ArticleGoogle Scholar
- Chuang YC, Wu HY, Lin YL, Tzou SC, Chuang CH, Jian TY, Chen PR, Chang YC, Lin CH, Huang TH, et al. Blockade of ITGA2 induces apoptosis and inhibits cell migration in gastric Cancer. Biol Proced Online. 2018;20:10.Google Scholar
- Dong J, Wang R, Ren G, Li X, Wang J, Sun Y, Liang J, Nie Y, Wu K, Feng B, et al. HMGA2-FOXL2 Axis regulates metastases and epithelial-to-mesenchymal transition of chemoresistant gastric cancer. Clin Cancer Res. 2017;23:3461–73.View ArticleGoogle Scholar
- Rozengurt E, Sinnett-Smith J, Eibl G. Yes-associated protein (YAP) in pancreatic cancer: at the epicenter of a targetable signaling network associated with patient survival. Signal Transduct Target Ther. 2018;3:11.View ArticleGoogle Scholar
- Kessenbrock K, Plaks V, Werb Z. Matrix metalloproteinases: regulators of the tumor microenvironment. Cell. 2010;141:52–67.View ArticleGoogle Scholar
- Hanahan D, Coussens LM. Accessories to the crime: functions of cells recruited to the tumor microenvironment. Cancer Cell. 2012;21:309–22.View ArticleGoogle Scholar