Next Article in Journal
High-Throughput Screening Identified Compounds Sensitizing Tumor Cells to Glucose Starvation in Culture and VEGF Inhibitors In Vivo
Next Article in Special Issue
Preoperative Imaging Evaluation after Downstaging of Pancreatic Ductal Adenocarcinoma: A Multi-Center Study
Previous Article in Journal
Extracellular Vesicles as Transmitters of Hypoxia Tolerance in Solid Cancers
Previous Article in Special Issue
Prognostic Role of High-Grade Tumor Budding in Pancreatic Ductal Adenocarcinoma: A Systematic Review and Meta-Analysis with a Focus on Epithelial to Mesenchymal Transition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Integrative Data Mining and Omics-Based Translational Model for the Identification and Validation of Oncogenic Biomarkers of Pancreatic Cancer

1
College of Pharmacy, Seoul National University, Seoul 08826, Korea
2
Department of Biomedical Sciences, College of Medicine, Inha University, 3-ga, Sinheung-dong, Jung-gu, Incheon 400-712, Korea
3
School of Medicine, Vietnam National University, Ho Chi Minh 70000, Vietnam
4
Department of Statistics, Seoul National University, Seoul 08826, Korea
5
Department of Medicine, College of Medicine, Inha University, 3-ga, Sinheung-dong, Jung-gu, Incheon 400-712, Korea
6
Division of Life and Pharmaceutical Sciences, Ewha Womans University, Seoul 120-750, Korea
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Cancers 2019, 11(2), 155; https://doi.org/10.3390/cancers11020155
Submission received: 14 December 2018 / Revised: 18 January 2019 / Accepted: 21 January 2019 / Published: 29 January 2019
(This article belongs to the Special Issue Advances in Pancreatic Cancer Research)

Abstract

:
Substantial alterations at the multi-omics level of pancreatic cancer (PC) impede the possibility to diagnose and treat patients in early stages. Herein, we conducted an integrative omics-based translational analysis, utilizing next-generation sequencing, transcriptome meta-analysis, and immunohistochemistry, combined with statistical learning, to validate multiplex biomarker candidates for the diagnosis, prognosis, and management of PC. Experiment-based validation was conducted and supportive evidence for the essentiality of the candidates in PC were found at gene expression or protein level by practical biochemical methods. Remarkably, the random forests (RF) model exhibited an excellent diagnostic performance and LAMC2, ANXA2, ADAM9, and APLP2 greatly influenced its decisions. An explanation approach for the RF model was successfully constructed. Moreover, protein expression of LAMC2, ANXA2, ADAM9, and APLP2 was found correlated and significantly higher in PC patients in independent cohorts. Survival analysis revealed that patients with high expression of ADAM9 (Hazard ratio (HR)OS = 2.2, p-value < 0.001), ANXA2 (HROS = 2.1, p-value < 0.001), and LAMC2 (HRDFS = 1.8, p-value = 0.012) exhibited poorer survival rates. In conclusion, we successfully explore hidden biological insights from large-scale omics data and suggest that LAMC2, ANXA2, ADAM9, and APLP2 are robust biomarkers for early diagnosis, prognosis, and management for PC.

1. Introduction

Pancreatic cancer (PC) is the seventh leading cause of cancer-related death worldwide, with less than seven percent of patients surviving after five years from the initial diagnosis [1]. Early diagnosis at a curable stage of PC is exceptionally challenging. Resection surgery may increase the 5-year survival rate to 20%, but only 10–20% of the patients are eligible for this procedure [2]. The symptoms of PC are usually not apparent due to its location deep inside the abdomen, and there are no current biomarkers with sufficient sensitivity and specificity for accurate diagnosis [3]. The cut-off level of carbohydrate antigen (CA) 19-9 for differentiating malignant and benign tumors is still controversial due to unconvincing diagnostic values. Additionally, the sensitivity of CA 19-9 drops dramatically in early and small tumors [4]. The lack of effective screening strategies and the paucity of powerful systematic therapies are also significant causes of PC fatality [5]. Collectively, novel methods for early detection and effective management of PC are urgently needed.
Only a small fraction of an enormous number of cancer biomarkers has been successfully implemented in the clinical practice. Cancer is characterized by tumor heterogeneity and a variety of molecular biosignatures; this calls for a more effective strategy of simultaneously combining a set of biomarkers [6]. A study has been reported on the improvement in diagnostic sensitivity of cancer using panels of combined biomarkers rather than individual molecules [7]. In PC, novel approaches applying methylated DNAs, circulating tumor DNAs, miRNAs, proteins, and metabolites have been developed and used with CA 19-9 in an attempt to enhance diagnostic and prognostic performance [8,9,10]. For instance, the sensitivity and specificity of the combined CA 19-9 and serum cytokines panel for differentiating PC from the normal group were 85.7% and 92.3%, respectively, and those of the combined CA19-9, ICAM-1, and OPG panel were 88% and 90%, respectively [11,12]. However, it is worth mentioning that several confounding factors may affect the establishment of a cancer biomarker panel. The samples for research on PC usually stem from patients with the advanced disease, as most patients are diagnosed at late stages, whereas biomarkers are expected to be for early, curable stages [13]. Additional clinical parameters, notably jaundice, are also worth considering when recruiting PC samples for biomarker studies [14].
Recent large-scale gene expression data studies have been conducted to shed light on molecular characterization and to provide a system-level view of human cancers, with a minor focus on PC [15,16]. Remarkably, a compendium including 2516 genes that showed changes in expression levels associated with PC was investigated, which accounted for nearly 13% of the known human coding genes [17]. However, the huge data load makes it impossible to translate the results into clinical applications. In this context, a promising and practical set derived from this database is needed. In the current study, we applied high-throughput next-generation sequencing (NGS)-based gene expression analysis, transcriptome meta-analysis, various data mining techniques with the previously published compendium and databases, and immunohistochemistry (IHC) on various cohorts for candidate validation to introduce multiplex signatures capable of differentiating PC tissue from the normal tissue. A powerful statistical learning technique and easy-to-comprehend explanation approach, using local interpretable model-agnostic explanations algorithm, were also utilized to enhance the diagnostic practice. The impact of the proposed biomarkers on patients’ survival was examined using The Cancer Genome Atlas (TCGA) cohort. Ultimately, an exceptional panel of biomarkers was introduced and initially validated. We also aimed to systematically assess the most promising therapeutic targets for PC. Finally, the associated biological and functionally altered signaling cascades of the introduced biosignatures in PC carcinogenesis have been discussed.

2. Results

2.1. Integrative Screening Approach Suggests 23 Novel Biomarker Candidates for the Diagnosis and Management of Pancreatic Cancer

We conducted an effect size meta-analysis from two datasets to detect differentially expressed genes. Dataset GSE16515 included 36 cancerous samples and 16 adjacent normal tissues. Dataset GSE28735 contained 45 matching pairs of pancreatic cancer and adjacent non-tumor tissues. Table S1 lists the characteristics of all utilized microarray data. We identified 3915 significantly upregulated and 4007 significantly downregulated genes in the PC tissue compared to normal tissue. It is worthy to note that we targeted upregulated genes in an effort to identify oncogenic biomarkers for the diagnosis, prognosis, and treatment of PC. It turned out that less than 7% (266 out of 3915) of overexpressed genes demonstrated a combined effect size (cES) greater than or equal to 1.5, the selection criterion for potential candidates (Figure 1a). In parallel, we utilized an exhaustive compendium of PC biomarkers collected by Harsha et al. as the filtering resource [17]. The genes coding for secretory proteins and membrane proteins overexpressed in PC, precursor lesions, the stroma associated with PC, and other subtypes of PC, but not in pancreatitis, were selected. Of note, we tried to limit our study to genes that are only overexpressed in PC but not in chronic pancreatitis, which are also those having potential to distinguish PC from chronic pancreatitis, since it may be difficult to distinguish chronic pancreatitis from PC in some cases, as there have been no excellent biomarkers to distinguish these two conditions [17,18]. Given these criteria, 95 secretory and 114 membrane protein-coding genes from the compendium fulfilled the criteria. Based on the Venn diagram analysis of the candidates from the meta-analysis and the compendium, five secretory (AGR2, GAPDH, LAMC2, MMP11, and TAGLN2), 12 membrane (CD82, CLDN18, EPHA2, EZR, FXYD3, GPRC5A, ITGA2, ITGB6, MET, MST1R, NQO1, and SLC2A1), and six “dual” (ADAM9, ANXA2, APLP2, CDH3, MSLN, and SERPINB5) protein-coding genes were selected for further analysis (Figure 1b,c).

2.2. Literature- and Experiment-Based Validation of the 23 Selected Candidates

Published papers and data submitted to public repositories were extensively mined for literature-based validation. The overexpression of the 23 candidates in pancreatic ductal adenocarcinoma (PDAC) was experimentally identified and validated by multiple techniques, such as microarray, quantitative proteomics, and antibody-based methods. According to the results, each gene was overexpressed with a change greater than two-fold at least once in PC or PC cell lines at the mRNA or protein level. For instance, ANXA2 was preferably reported to be overexpressed in PDAC tissue in the mRNA analysis using DNA Microarray and in the protein analyses using IHC, isotope-coded affinity tag, western blot, mass spectrometry, two-dimensional electrophoresis, and label-free quantitation mass spectrometry.
We further conducted an NGS gene expression analysis on a cohort containing six PDAC samples and six controls for experiment-based validation. Despite the heterogeneity and other confounding factors, the NGS result was promising, showing 17 out of 23 consistently overexpressed genes with a cut-off of two-fold change (log FC = 1.00), and 14 out of 23 genes using log FC = 1.50. Remarkably, SERPINB5 was more than 250 times overexpressed in the PC tissue compared to normal tissue (log FC = 8.35). LAMC2, CDH3, and ANXA2 exhibited the altered expression difference with log FC of 3.91, 2.72, and 1.16, respectively. Collectively, the experimental and in silico analyses largely supported our 23 selected candidates as reliable biomarkers for differentiating PC tissue from normal tissue.
Genetic alterations occur commonly in patients with curable pancreatic pre-malignant lesions, making the biomarker discovery at early stages indispensable [19]. The altered expression of the 23 genes in pancreatic intraductal papillary-mucinous neoplasm (IPMN) and intraductal papillary-mucinous carcinoma (IPMC) was also explored. Among the selected candidates, AGR2, LAMC2, TAGLN2, MSLN, SERPINB5, CLDN18, EZR, FXYD3, GPRC5A, ITGA2, MST1R, NQO1, and SLC2A1 were upregulated in pancreatic neoplasms compared to normal tissue, while CDH3, EPHA2, GPRC5A, and SLC2A1 were upregulated in PC compared to pancreatic neoplasms, providing clues that they may be involved in PC initiation and development.
The list of techniques for gene expression, identification, and validation, cES, and gene expression in premalignant pancreatic lesions are summarized in Table 1.

2.3. ADAM9, ANXA2, ITGA2, MET, and LAMC2 Exhibit Substantial Impact on the Survival of PC Patients

The Kaplan-Meier (KM) survival analysis and the multivariate Cox regression analysis were conducted based on the available data from TCGA pancreatic cancer cohort. The KM survival analysis revealed that PC patients with high expression levels of 17 of the candidate genes exhibited a poorer overall survival (OS) or disease-free survival (DFS). Particularly, patients overexpressing ADAM9, ANXA2, ITGA2, or MET, had both poorer OS and DFS rates than those with low gene expression. Figure 2 shows the KM plots of the OS for ADAM9, ANXA2, ITGA2, and MET expression. Of note, the HR was greater than 2 for both analyses, except for ITGA2 that had an HR of 1.9 for the DFS analysis. LAMC2 exhibited a significantly poor DFS (HR = 1.8, p-value = 0.012). A multivariate Cox regression analysis of the TCGA cohort was adopted to evaluate the impact of the 23 genes in patient survival. The analysis revealed that LAMC2, ADAM9, ANXA2, CDH3, SERPINB5, EPHA2, GPRC5A, ITGA2, ITGB6, and MET were likely associated with the worse outcome of PC with FDR < 0.05. The detailed information regarding Cox regression survival analysis of these genes is described in Table 2 and Figure S1.

2.4. Supervised Machine Learning Classification Demonstrates that ADAM9, ANXA2, LAMC2, and APLP2 Are Important in Differentiating PC Tissue from Normal Tissue

The role of 11 secretory protein-coding genes for differentiating PC tissue from normal tissue was further tested. We visualized the expression profiling of the 11 genes by means of boxplots. For instance, the expression levels of two promising candidates, CDH3 and LAMC2, were significantly higher in PDAC compared to normal tissue (Figure 3a,b). Boxplots of the remaining potential genes can be found in Figure S2. As illustrated in Figure 3c, two classes of PC and normal tissues showed a relatively separated tendency in the principal component analysis. However, there was still an overlapped region between them, making grouping of those cases more challenging. A heatmap analysis further provided a general overview of relative differences in gene expression between PC and normal tissue (Figure 3d). As expected, in general, the 11 candidates exhibited a higher level of expression in PC tissues than in normal tissue. The dataset was then randomly divided into the training set (70%) and test set (30%) for supervised RF model calibration and validation. RF classification was conducted to build the prediction model based on 11 potential genes. The number of randomly chosen variables at each split ranged from one to 11, and the data were cross-validated tenfold to find the best tuning parameter, which was repeated five times. The optimal RF model (mtry of two) in our trial revealed an excellent performance with an area under the ROC curve of 0.96 in the training set. The model validated in the test set also exhibited promising results with only one misclassified case, reporting accuracy of 0.96, sensitivity of 1.00, and specificity of 0.93 (Figure 3e). Of note, LAMC2, ANXA2, APLP2, and ADAM9 were the most important variables of the model, although other genes facilitated the decision (Figure 3f).
Local Interpretable Model-Agnostic Explanations (LIME) algorithm was applied to explain the “decision rules” of our black-box model. The first case was classified into the normal group, with a probability of 0.96 based on the gene expression rules: ANXA2 ≤ 9.77, APLP2 ≤ 8.31, ADAM9 ≤ 6.74, 5.31 < LAMC2 ≤ 6.07, and 3.61 < SERPINB5 ≤ 4.46. Another case was certainly predicted as PC based on the rules: LAMC2 > 7.82, ANXA2 > 10.62, ADAM9 > 8.09, MMP11 > 6.58, and AGR2 > 9.68, reaching the probability of 1.00. We deduced from the typical cases that a higher expression level of ANXA2 (cut-off threshold = 9.77), ADAM9 (cut-off threshold = 6.74), LAMC2 (cut-off threshold = 6.07), AGR2 (cut-off threshold = 9.68), MMP11 (cut-off threshold = 6.58), APLP2 (cut-off threshold = 8.31), and SERPINB5 (cut-off threshold = 4.46) was an indication of PC in this investigation. The visualization of the above explanations of the four representative cases is given in Figure S3.

2.5. Immunohistochemistry Analysis Reveals High Protein Expression Levels of LAMC2, ADAM9, ANXA2, and APLP2 in Pancreatic Cancer

We assessed the expression of LAMC2, ADAM9, and ANXA2 at the protein level using tissue microarray samples of 86 patients and normal controls. They were chosen due to their roles in the diagnostic model, prognostic impact, and profound validated evidence. Typical immunostaining is shown in Figure 4a, and the cohort characteristics are given in Table S2. Individual proteins showed moderate to strong immunostaining in most PC tissues, and weak or negative immunostaining in normal tissues; the differences are statistically significant (Figure 4b, p-value < 0.001). LAMC2, ADAM9, and ANXA2 exhibited moderate correlations with their protein levels according to Spearman’s pairwise correlation (Figure 4c, Figures S4 and S5). The protein expression of LAMC2, ADAM9, and ANXA2 between PC and adjacent tissues (eight paired samples) presented similar results, although only ADAM9 showed statistical significance (Figure S5, p-value = 0.031). The protein expression of APLP2 was also investigated in a cohort comprising 64 PC and 64 adjacent normal tissues (Table S3). As shown in Figure S6, the expression of APLP2 was significantly higher in PC tissues (p-value = 0.0139), making them valuable for further investigation. The Pearson’s correlation analysis of gene expression between LAMC2ADAM9 (R = 0.71, p-value < 0.001), ANXA2LAMC2 (R = 0.75, p-value < 0.001), and ANXA2ADAM9 (R = 0.78, p-value < 0.001) revealed strong correlations among them (Figure 4d).

2.6. Functional Analysis and Public Data Mining to Create a Comprehensive Picture of the Biological Processes in Pancreatic Cancer

Cancer Hallmarks Analytics Tool (CHAT) was employed to clarify the involvement of our candidates in the biological processes leading to PC according to the current ten hallmarks of cancer. Surprisingly, 18 of the 23 genes were implicated in at least five cancer hallmarks, suggesting that they may be highly associated with PC carcinogenesis and progression. Finally, a pathway enrichment analysis was performed based on the 397 upregulated genes with a cES greater than 1.5, derived from the meta-analysis of five datasets, to provide solid evidence on the PC related biological processes. As a result, various KEGG-annotated pathways were enriched in the PC tissue compared to normal tissue, including prominent cancer-related pathways, such as cell cycle and intrinsic molecular processes of PC and others, such as prostate cancer and chronic myeloid leukemia. The results indicate both unique and overlapping features of PC and other cancers, which would be helpful for further studies. In the protein–protein interaction network of 397 upregulated genes, strong connections among the nodes were observed (score > 0.900). There were two major clusters formed by pathways in cancer and cell cycles in the network (Figure S7). Pathway enrichment, protein–protein networks with known and predicted interactions between the nodes, and CHAT analysis can be found in Tables S4 and S5, and Figure S7.

2.7. Drug–Gene Interaction and miRNA–Gene for Further Development of Novel Therapeutics

Increasing evidence suggests that aberrantly expressed miRNAs are attractive targets for preventing PC initiation, progression, metastasis, and chemoresistance. Given their enormous potential as therapeutic tools for cancer management, we checked whether the transcripts of our 18 membrane protein-coding genes were the targets of miRNAs based on experimentally validated evidence. It turned out that nine candidates exhibited a strong connection with miRNAs (ADAM9, ANXA2, MSLN, SERPINB5, EZR, GPRC5A, ITGA2, MET, and SLC2A1). The Pearson’s correlation coefficient between each candidate and the miRNA was also recorded. Finally, we input our genes in curated resources combined by the drug-gene interaction database (DGIdb) to examine if they coded for a drug-metabolizing enzyme or a drug transporter (so-called druggable genes), and then determined possible drug-gene interactions. As an example of the results, ADAM9 was positively correlated with hsa-miR-126-3p (the correlation coefficient, 0.315; p-value < 0.05) with strong evidence from the reporter assay, western blot, and qPCR. According to DrugBank, the drug Ilomastat was detected to inhibit the cell-surface protein ADAM9. The detailed information regarding miRNA–gene interaction and drug–gene interaction is presented in Table 3.

3. Discussion

PC is a complicated and highly heterogeneous disease, yet it takes approximately two decades for a primary pancreatic tumor to develop from initiation to cancer death, opening a large window for an effective screening and early detection program [20]. In that context, a plethora of biomarkers have been studied, but most of them have failed to be translated into clinical practice. Robust biomarkers, even identified by advanced, reliable technologies, have not exhibited a better diagnostic and prognostic ability than conventional biomarkers in large-scale clinical trials [21]. In recent years, advances in computational techniques have facilitated and enhanced cancer omics-based research [22]. The current study utilizes an unbiased data-driven approach together with the integration of omics data, machine learning, and conventional immunohistochemistry approaches to suggest and validate potential biomarkers for PC diagnostics and management. Our approach covers a broad spectrum of cancer biology, of which LAMC2, ANXA2, ADAM9, and APLP2 are the most important features of the tested diagnostic model, and their implementation in clinical settings could be further examined. The roles of LAMC2, ANXA2, ADAM9, and APLP2 in cancer biology have been documented [23,24,25,26,27]. However, our study has been a pioneer in suggesting the combination of LAMC2, ANXA2, ADAM9, and APLP2 and statistical learning algorithms to properly detect PDAC in early stages. Of note, the expression at the gene and protein level of LAMC2, ANXA2, and ADAM9 are positively correlated to PC. PC patients overexpressing LAMC2 potentially have a higher chance of getting distant metastasis and poorer prognosis [28]. The elevation of LAMC2 in serum was also considered an indicator for differentiating healthy people from an early stage of PDAC with outstanding performance [29]. Similarly, the high expression of stromal ANXA2 was significantly correlated with short disease-free survival and overall survival, while the high expression of ADAM9 was associated with poor tumor differentiation and short overall survival in PDAC patients [26]. We also screened for the presence of somatic mutations of our candidates using the most comprehensive global database for the annotation of somatic variants in human cancer, the Catalogue of Somatic Mutations in Cancer (COSMIC) [30]. COSMIC has systematically curated published cancer mutation data and provides in-depth knowledge on cancer genes, including their involvement in cancer genesis and mutation mechanisms. It turned out that our candidates were not frequently mutated genes. Further studies integrating the suggested biosignature and conventional biomarkers, such as CA 19-9, are needed to validate the performance of the panel. Nevertheless, our investigation demonstrates that with the power of machine learning and omics-based big data, scientists could gather patient tumor data and clinical parameters together with biomarker screening using the combined panel to explore the pattern of alterations in multiple PDAC patient tumors, thereby increasing the diagnostic and prognostic accuracy.
Out of the 23 genes in our panel, 18 were involved in at least five cancer hallmarks, indicating the relation of these genes and associated products with cancer initiation and progression. However, with respect to panel, the roles of immune destruction, cellular energetics, and replicative immortality in PC have not been studied frequently compared with the other hallmarks. More studies regarding these aspects are needed to fully explore the potential of our multiplex panels in PC diagnosis and treatment. In recent years, extensive efforts have been made to shed light on the association between immune destructions, cancer, and associated regulators [31,32]. The immune hallmark can be summarized as the ability to survive in the chronically infiltrated microenvironment, to escape the immune recognition, and to suppress the immune response [33]. Despite their great impact on cancer treatment and prevention, these processes may also give some clues for developing new screening and diagnostic approaches. Of note, the deregulating cellular energetics or cancer metabolism has been recently recognized as a cancer hallmark [34]. Investigations in the past 20 years and the successful implementation of metabolic tumor imaging into clinical practice have once again emphasized that altered metabolic processes are critical for tumor survival and proliferation [35]. Among them, the alterations of glycolysis, amino acids, and lipid metabolism are the most important changes of cancer cellular energetics [36]. Multiple metabolic enzymes and pathways could be investigated as attractive targets not only for anticancer therapeutics but also for developing novel biomarkers [37].
It was expected that miRNA would diminish drug resistance of tumors through various regulatory mechanisms [38,39]. In our investigation, ADAM9 and ANXA2 were attractive candidates for miRNA-based anti-PC therapies. They exhibited strong interactions with miR-126, miR-33a, and miR-125a (ADAM9); and miR-155 and miR-206 (ANXA2). In PC, miR-126 regulated ADAM9 by targeting the epithelial-mesenchymal transition, the initial step of the metastatic cascade [40]. ANXA2 and KRAS were also direct downstream targets of miR-206 that modulated cancer cell growth, metastasis, and lymphangiogenesis in PDAC [41]. In addition to miRNAs, various drug molecules are worthy of more attention. Lately, in an effort to implement evidence-based personalized treatment recommendations into clinical settings, some software and databases have been developed to assist the decision of PC treatment. A recent evidence mining software analyzing next-generation sequencing data suggests pharmacogenomic biomarkers and has potential to predict PC treatment efficacy and toxicity [42]. Further studies are necessary to fully explore the potential of ADAM9 and ANXA2, with high level of evidence of clinical validity, among other candidates, as drug candidates in PC management with regard to their associated miRNAs.
The current study set out to introduce an integrative data mining model for translational research to facilitate the discovery and validation of new biomarkers. In particular, a comprehensive literature and experimental data mining approach to examine potential roles of overexpressed genes as biomarkers for diagnosis and prognosis of pancreatic cancer were conducted. Thus, as a limitation, down-regulated genes were neglected in all of our analyses, which might result in the missing of valuable candidates for further validation. In addition, a follow-up study integrating four proposed candidates and statistical learning is needed to examine their actual clinical validity.

4. Materials and Methods

4.1. Ethical Approval and Consent of Patients

The Institutional Review Board (IRB) of Seoul National University reviewed and exempted our protocol for a full ethical approval application when using human gene expression data from the public database without personal information, survival data, and commercial tissue micro-array (SNU 16-11-014 and SNU 17-05-043). The IRB of Inha University approved the protocol of the NGS gene expression analysis (INHAUH 2016-08-008-001). Every patient signed an informed consent before any procedures were conducted. The experiments followed the Ethical Principles of the World Medical Association Declaration of Helsinki.

4.2. Antibodies

ADAM9 antibody was purchased from Invitrogen (PA5-25959, Waltham, MA, USA), Annexin A2 (ANXA2) antibody was purchased from Cell Signaling (#8235, Beverly, MA, USA), and LAMC2 and APLP2 antibodies were purchased from Sigma (AMAb91098 and HPA039319, St. Louis, MO, USA).

4.3. Pancreatic Cancer and Normal Control Cohorts

4.3.1. NGS Sample Collection

Tissue samples were obtained from six PC patients, and six paired controls (stage IIa/IIb, median age of 66.5, range from 55 to 75, male:female of 2:1) using standard protocols of the Inha University Hospital. All samples were formalin-fixed paraffin-embedded and stored at 4 °C.

4.3.2. Microarray Sample Collection

The Affymetrix-based datasets GSE16515, GSE28735, GSE15471, GSE18670, and GSE41368, were included in the gene expression meta-analysis. After that, we divided the datasets into two groups with similar sample sizes: set A including, GSE16515 and GSE28735 for variable selection purposes; and set B, including GSE15471, GSE18670, and GSE41368 for statistical learning-based validation of the proposed set of novel biomarkers. Dataset GSE19650 was also included for the examination of the multistep PC development.

4.3.3. Tissue-Array Sample Collection

Tissue microarrays containing PC tissues, adjacent normal tissues, and normal tissues (HPan-A150CS-02 and PA1001b) were purchased from US Biomax Inc. (Rockville, MD, USA).

4.4. RNA Extraction from Formalin-Fixed Paraffin-Embedded Tissue

Formalin-fixed paraffin-embedded tissue was cut using a razor blade (20 µm), transferred to a 1.5 mL tube, and filled with mineral oil. The tube was incubated at 70 °C for 3 min to melt the paraffin. The sample was then centrifuged for 1 min at 16,000× g. The pellet was washed with 1 mL of ethanol and air dried. Next, 150 µL protease K digestion buffer containing 500 µg/mL protease K was added, and the mixture was incubated for 3 h. 1 mL of Trizol was added to the sample and incubated at room temperature for at least 5 min to dissociate nucleoprotein complexes. Subsequently, 0.2 mL chloroform was added, and the mixture was vigorously vortexed for 15 s and centrifuged at 12,000× g for 10 min at 4 °C. The aqueous phase was transferred to a fresh tube. Total RNA was precipitated by mixing with an equal volume of isopropyl alcohol at −20 °C for at least 1 h, followed by 10 min of centrifugation at 12,000× g at 4 °C. Finally, RNA pellet was washed with 70% ethanol, briefly air-dried, and resuspended in RNase-free water.

4.5. RNA Data Processing and Analysis

To construct cDNA libraries, 1 µg of total RNA was used with TruSeq RNA library kit (Illumina, San Diego, CA, USA). The protocol included polyA-selected RNA extraction, RNA fragmentation, random hexamer-primed reverse transcription, and 100-nt paired-end sequencing by Illumina HiSeq2500 (Illumina). The libraries were quantified using qPCR and qualified using an Agilent Technologies 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). To estimate the expression level, the RNA-seq reads were mapped to the genome of Homo sapiens. The reference genome sequence of Homo sapiens and the annotation data were downloaded from the UCSC website (http://genome.uscs.edu). The transcript counts at the gene level were calculated [43]. Data normalization and statistical analysis were then conducted using edgeR. We determined the significant result by adjusting the absolute value of fold change ≥ 2.

4.6. Microarray Data Processing and Gene Expression Meta-Analysis

We normalized raw gene expression data using a robust multi-array analysis implemented in the Bioconductor affy package [44]. Database for Annotation, Visualization, and Integrated Discovery (DAVID) version 6.8 was applied to get the Entrez IDs. The Empirical Bayes (ComBat) cross-study normalization method was applied to remove batch effects. The random effects model was chosen as the statistical meta-analysis method based on the result of Cochran’s Q test. The analysis was implemented in the comprehensive web-based tool NetworkAnalyst [45].

4.7. Immunohistochemistry (IHC) Experiments

The procedures of IHC were carried out as described previously [46]. Two experts independently assessed staining scores. The immunostaining intensity was graded as 0, 1, 2, or 3 for the tissue with stainless cells, <10% stained cells, 10–50% stained cells, and > 50% stained cells, respectively. We excluded the tissues that were not suitable for analysis (e.g., not adenocarcinoma or tissue loss cores). IHC scoring was visualized by Aperio ScanScope digital slide scanners, and analyzed by the vendor’s software (Leica, Wetzlar, Germany). The Mann-Whitney and Wilcoxon matched-pairs signed rank tests were conducted to test two-group differences among PC tissues, adjacent normal tissues, and normal tissues using GraphPad Prism 6 (San Diego, CA, USA) where applicable.

4.8. Unsupervised and Supervised Machine Learning Algorithms

4.8.1. Variable Selection

We treated the overexpressed genes coding for secretory and membrane proteins separately because of their distinct clinical roles. Only those coding for secretory proteins and not upregulated in pancreatitis from the compendium were selected for a further classification.

4.8.2. Data Exploration and Visualization

The principal component analysis (PCA) was incorporated to facilitate data visualization and to detect potential outliers prior to the class assignment analysis using FactoMineR version 1.35 and factoextra version 1.0.4 [47]. A heatmap was also applied to highlight the differences in gene expression between the two comparison groups by means of MetaboAnalyst 3.0 [48].

4.8.3. Random Forests Classification Model and Explanation

Classification and Regression Training (caret) package version 6.0.77 was used for data splitting and supervised classification analysis [49]. Set B was divided into the training set and the test set with a ratio of 7:3. A random forests (RF) model using 500 trees was optimized on the training set. The number of variables for splitting at each tree node (mtry) was tuned. Five-time repeated 10-fold cross-validation sampling iteration was used during the training process to select the optimal model. A seed number was given prior to data splitting and model training to achieve reproducible results. We used the explainer local interpretable model-agnostic explanations (LIME) to interpret the rationale behind the black-box RF model’s predictions [50].

4.9. Correlation Analysis

Pearson’s correlation analysis and Spearman’s rank correlation analysis were conducted for gene expression level and protein expression score, respectively.

4.10. Kaplan-Meier Plots and Cox Regression Analysis

The overall survival and disease-free survival were investigated using the Kaplan-Meier method with the log-rank test. We set the high and low gene expression level groups by the median value. The overall survival plot and disease-free survival plot were obtained with the hazard ratios (HR) and the 95% confidence interval information. The whole process was implemented using the web-based tool GEPIA [51]. Additionally, multivariate Cox regression analysis of TCGA PC patients was also extracted using Oncolnc to determine the potential prognostic biomarkers [52].

4.11. Database Mining, Pathway Enrichment, STRING, and CHAT Analyses

Kyoto Encyclopedia of Genes, Genomes (KEGG) pathway enrichment, and OmicsNet and STRING protein-protein interaction analyses and visualizations were conducted using the significantly overexpressed genes. The Database of Human Pancreatic Cancer (HPCDb) and the Pancreatic Expression Database (PED) were used as resources for literature data mining [53,54]. A public literature review of a reported gene list in PC was also conducted under the guidance of Cancer Hallmarks Analytics Tool (CHAT) [55].

4.12. Drug-Gene Interaction and miRNA–Gene Interaction for Further Developing Novel Therapeutics

To examine the connection between our membrane protein-coding genes with miRNA, we employed the miRTarBase [56]. The miRTarBase offers data regarding experimentally validated miRNA–target interactions. The Pearson’s correlation coefficient between each candidate and miRNA was also recorded. Subsequently, the drug-gene interaction database was investigated to examine if our genes coded for a drug-metabolizing enzyme or a drug transporter [57].

4.13. Statistical Significance Level

R statistic 3.4.2 was used to implement the analysis except otherwise stated. A p-value of 0.05 was used as the cut-off for significance in statistical tests. False discovery rate (FDR) (Benjamini-Hochberg method), when applicable (e.g., in NGS differential analysis), with a threshold of 0.05, was utilized for all multiple hypothesis tests.

4.14. Availability of Data and Materials

NGS data generated and analyzed during the current study are available in the Gene Expression Omnibus (GEO) with the assession code of GSE119224. Other datasets analyzed during the current study are available in the GEO (GSE16515, GSE28735, GSE15471, GSE18670, GSE41368, and GSE19650). All data generated or analyzed during this study are included in this published article and its supplementary information files.

5. Conclusions

Recent advances in high-throughput multi–omic technologies along with advanced bioinformatics and machine learning techniques have provided a new paradigm for the discovery of human disease biomarkers [10,46,58,59]. This perspective should remark the continued employment of multi-omic approaches and machine learning techniques under the guidance of big data and biological mechanisms in cancer biomarker research. In view of this, we successfully demonstrated an integrative method to introduce and validate a panel of LAMC2, ANXA2, ADAM9, and ALPL2 for innovating PC diagnosis, prognosis, and management. Future studies are expected to take a step forward and concentrate on validating biomarkers by integrating multi-omics data in an epidemiological and clinical context-dependent manner.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6694/11/2/155/s1. Figure S1: Overall survival and disease-free survival plots of the 23 candidates. Figure S2: Boxplots of the remaining nine secretory protein-coding genes; p-value < 0.0001 in all analyses. Figure S3: LIME of four representative cases. Figure S4: Spearman’s pairwise correlation analysis for ANXA2, ADAM9, and LAMC2. Figure S5: Protein expression level of three promising candidates in eight matched PC and normal samples. Figure S6: Protein expression level of APLP2 in 64 matched PC and normal samples. p-value < 0.05. Figure S7: Protein-protein interaction network of 397 upregulated genes in the meta-analysis. Table S1: Characteristics of the included microarray gene expression datasets. Table S2: Characteristics of the cohort used for assessing the expression of LAMC2, ADAM9, and ANXA2. Table S3: Characteristics of the cohort used for assessing the expression of APLP2. Table S4: Significantly enriched pathways of 397 upregulated genes from the meta-analysis of five datasets. Table S5: Summary of the CHAT analysis.

Author Contributions

Conceptualization, S.W.K., S.-S.H. and N.P.L.; methodology, N.P.L., K.H.J., S.P., N.H.A., J.E.M., H.M.K. and J.H.L.; software, N.P.L., N.H.A., T.D.N., S.P. and S.L.; formal analysis, N.P.L., K.H.J., N.H.A., H.H.Y., T.D.N., S.P., S.J.Y., J.E.M., H.M.K. and J.H.L.; investigation, S.W.K., S.-S.H., S.L., J.M.K. and J.L.; resources, S.W.K. and S.-S.H.; data curation, N.P.L., K.H.J., N.H.A., H.H.Y., T.D.N. and S.P.; writing—original draft preparation, N.P.L. and K.H.J.; writing—review and editing, all authors; visualization, N.P.L., K.H.J., N.H.A., H.H.Y., T.D.N. and S.P.; supervision, S.W.K., S.-S.H., S.L., J.L. and J.M.K.; funding acquisition, S.W.K. and S.-S.H.

Funding

This research was funded by National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIP) (NRF-2018R1A5A2024425), the National Research Foundation of Korea (NRF) grant (2018R1A2A1A05077263), and by BK21 Plus Program in 2018.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Balachandran, V.P.; Luksza, M.; Zhao, J.N.; Makarov, V.; Moral, J.A.; Remark, R.; Herbst, B.; Askan, G.; Bhanot, U.; Senbabaoglu, Y.; et al. Identification of unique neoantigen qualities in long-term survivors of pancreatic cancer. Nature 2017, 551, 512–516. [Google Scholar] [CrossRef] [PubMed]
  2. Raimondi, S.; Maisonneuve, P.; Lowenfels, A.B. Epidemiology of pancreatic cancer: An overview. Nat. Rev. Gastroenterol. Hepatol. 2009, 6, 699–708. [Google Scholar] [CrossRef] [PubMed]
  3. Kleeff, J.; Korc, M.; Apte, M.; La Vecchia, C.; Johnson, C.D.; Biankin, A.V.; Neale, R.E.; Tempero, M.; Tuveson, D.A.; Hruban, R.H.; et al. Pancreatic cancer. Nat. Rev. Dis. Primers 2016, 2, 16022. [Google Scholar] [CrossRef] [PubMed]
  4. Decker, G.A.; Batheja, M.J.; Collins, J.M.; Silva, A.C.; Mekeel, K.L.; Moss, A.A.; Nguyen, C.C.; Lake, D.F.; Miller, L.J. Risk factors for pancreatic adenocarcinoma and prospects for screening. Gastroenterol. Hepatol. 2010, 6, 246–254. [Google Scholar]
  5. Semaan, A.; Maitra, A. Pancreatic cancer in 2017: Rebooting pancreatic cancer knowledge and treatment options. Nat. Rev. Gastroenterol. Hepatol. 2018, 15, 76–78. [Google Scholar] [CrossRef] [PubMed]
  6. Sant, G.R.; Knopf, K.B.; Albala, D.M. Live-single-cell phenotypic cancer biomarkers-future role in precision oncology? NPJ Precis. Oncol. 2017, 1, 21. [Google Scholar] [CrossRef] [PubMed]
  7. Russell, M.R.; Graham, C.; D’Amato, A.; Gentry-Maharaj, A.; Ryan, A.; Kalsi, J.K.; Ainley, C.; Whetton, A.D.; Menon, U.; Jacobs, I.; et al. A combined biomarker panel shows improved sensitivity for the early detection of ovarian cancer allowing the identification of the most aggressive type II tumours. Br. J. Cancer 2017, 117, 666–674. [Google Scholar] [CrossRef] [Green Version]
  8. Lee, H.S.; Jang, C.Y.; Kim, S.A.; Park, S.B.; Jung, D.E.; Kim, B.O.; Kim, H.Y.; Chung, M.J.; Park, J.Y.; Bang, S.; et al. Combined use of CEMIP and CA 19-9 enhances diagnostic accuracy for pancreatic cancer. Sci. Rep. 2018, 8, 3383. [Google Scholar] [CrossRef]
  9. Chan, A.; Prassas, I.; Dimitromanolakis, A.; Brand, R.E.; Serra, S.; Diamandis, E.P.; Blasutig, I.M. Validation of biomarkers that complement CA19.9 in detecting early pancreatic cancer. Clin. Cancer Res. Off. J. Am. Assoc. Cancer Res. 2014, 20, 5787–5795. [Google Scholar] [CrossRef]
  10. Long, N.P.; Yoon, S.J.; Anh, N.H.; Nghi, T.D.; Lim, D.K.; Hong, Y.J.; Hong, S.-S.; Kwon, S.W. A systematic review on metabolomics-based diagnostic biomarker discovery and validation in pancreatic cancer. Metabolomics 2018, 14, 109. [Google Scholar] [CrossRef]
  11. Zeh, H.J.; Winikoff, S.; Landsittel, D.P.; Gorelik, E.; Marrangoni, A.M.; Velikokhatnaya, L.; Winans, M.T.; Lee, K.; Moser, A.; Bartlett, D.; et al. Multianalyte profiling of serum cytokines for detection of pancreatic cancer. Cancer Biomark. Sect. A Dis. Markers 2005, 1, 259–269. [Google Scholar] [CrossRef]
  12. Brand, R.E.; Nolen, B.M.; Zeh, H.J.; Allen, P.J.; Eloubeidi, M.A.; Goldberg, M.; Elton, E.; Arnoletti, J.P.; Christein, J.D.; Vickers, S.M.; et al. Serum biomarker panels for the detection of pancreatic cancer. Clin. Cancer Res. Off. J. Am. Assoc. Cancer Res. 2011, 17, 805–816. [Google Scholar] [CrossRef] [PubMed]
  13. Goymer, P. Early detection for pancreatic cancer. Nat. Rev. Cancer 2008, 8, 408. [Google Scholar] [CrossRef]
  14. Sener, S.F.; Fremgen, A.; Menck, H.R.; Winchester, D.P. Pancreatic cancer: A report of treatment and survival trends for 100,313 patients diagnosed from 1985–1995, using the National Cancer Database. J. Am. Coll. Surg. 1999, 189, 1–7. [Google Scholar] [CrossRef]
  15. Hu, J.; Locasale, J.W.; Bielas, J.H.; O’Sullivan, J.; Sheahan, K.; Cantley, L.C.; Heiden, M.G.V.; Vitkup, D. Heterogeneity of tumor-induced gene expression changes in the human metabolic network. Nat. Biotechnol. 2013, 31, 522. [Google Scholar] [CrossRef] [PubMed]
  16. Peng, X.; Chen, Z.; Farshidfar, F.; Xu, X.; Lorenzi, P.L.; Wang, Y.; Cheng, F.; Tan, L.; Mojumdar, K.; Du, D.; et al. Molecular Characterization and Clinical Relevance of Metabolic Expression Subtypes in Human Cancers. Cell Rep. 2018, 23, 255–269.e4. [Google Scholar] [CrossRef] [PubMed]
  17. Harsha, H.C.; Kandasamy, K.; Ranganathan, P.; Rani, S.; Ramabadran, S.; Gollapudi, S.; Balakrishnan, L.; Dwivedi, S.B.; Telikicherla, D.; Selvan, L.D.N.; et al. A Compendium of Potential Biomarkers of Pancreatic Cancer. PLoS Med. 2009, 6, e1000046. [Google Scholar] [CrossRef]
  18. Zhang, X.; Shi, S.; Zhang, B.; Ni, Q.; Yu, X.; Xu, J. Circulating biomarkers for early diagnosis of pancreatic cancer: Facts and hopes. Am. J. Cancer Res. 2018, 8, 332–353. [Google Scholar] [PubMed]
  19. Oldfield, L.E.; Connor, A.A.; Gallinger, S. Molecular Events in the Natural History of Pancreatic Cancer. Trends Cancer 2017, 3, 336–346. [Google Scholar] [CrossRef]
  20. Yachida, S.; Jones, S.; Bozic, I.; Antal, T.; Leary, R.; Fu, B.; Kamiyama, M.; Hruban, R.H.; Eshleman, J.R.; Nowak, M.A.; et al. Distant metastasis occurs late during the genetic evolution of pancreatic cancer. Nature 2010, 467, 1114. [Google Scholar] [CrossRef]
  21. Diamandis, E.P. The failure of protein cancer biomarkers to reach the clinic: Why, and what can be done to address the problem? BMC Med. 2012, 10, 87. [Google Scholar] [CrossRef] [PubMed]
  22. Ding, L.; Wendl, M.C.; McMichael, J.F.; Raphael, B.J. Expanding the computational toolbox for mining cancer genomes. Nat. Rev. Genet. 2014, 15, 556. [Google Scholar] [CrossRef] [PubMed]
  23. Garg, M.; Braunstein, G.; Koeffler, H.P. LAMC2 as a therapeutic target for cancers. Expert Opin. Ther. Targets 2014, 18, 979–982. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Mauri, P.; Scarpa, A.; Nascimbeni, A.C.; Benazzi, L.; Parmagnani, E.; Mafficini, A.; Della Peruta, M.; Bassi, C.; Miyazaki, K.; Sorio, C. Identification of proteins released by pancreatic cancer cells by multidimensional protein identification technology: A strategy for identification of novel cancer markers. FASEB J. Off. Publ. Fed. Am. Soc. Exp. Biol. 2005, 19, 1125–1127. [Google Scholar] [CrossRef] [PubMed]
  25. Zheng, L.; Jaffee, E.M. Annexin A2 is a new antigenic target for pancreatic cancer immunotherapy. Oncoimmunology 2012, 1, 112–114. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Grützmann, R.; Lüttges, J.; Sipos, B.; Ammerpohl, O.; Dobrowolski, F.; Alldinger, I.; Kersting, S.; Ockert, D.; Koch, R.; Kalthoff, H.; et al. ADAM9 expression in pancreatic cancer is associated with tumour type and is a prognostic factor in ductal adenocarcinoma. Br. J. Cancer 2004, 90, 1053. [Google Scholar] [CrossRef] [PubMed]
  27. Pandey, P.; Sliker, B.; Peters, H.L.; Tuli, A.; Herskovitz, J.; Smits, K.; Purohit, A.; Singh, R.K.; Dong, J.; Batra, S.K.; et al. Amyloid precursor protein and amyloid precursor-like protein 2 in cancer. Oncotarget 2016, 7, 19430–19444. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Takahashi, S.; Hasebe, T.; Oda, T.; Sasaki, S.; Kinoshita, T.; Konishi, M.; Ochiai, T.; Ochiai, A. Cytoplasmic expression of laminin γ2 chain correlates with postoperative hepatic metastasis and poor prognosis in patients with pancreatic ductal adenocarcinoma. Cancer 2002, 94, 1894–1901. [Google Scholar] [CrossRef] [PubMed]
  29. Kosanam, H.; Prassas, I.; Chrystoja, C.C.; Soleas, I.; Chan, A.; Dimitromanolakis, A.; Blasutig, I.M.; Ruckert, F.; Gruetzmann, R.; Pilarsky, C.; et al. Laminin, gamma 2 (LAMC2): A promising new putative pancreatic cancer biomarker identified by proteomic analysis of pancreatic adenocarcinoma tissues. Mol. Cell. Proteom. 2013, 12, 2820–2832. [Google Scholar] [CrossRef] [PubMed]
  30. Tate, J.G.; Bamford, S.; Jubb, H.C.; Sondka, Z.; Beare, D.M.; Bindal, N.; Boutselakis, H.; Cole, C.G.; Creatore, C.; Dawson, E.; et al. COSMIC: The Catalogue of Somatic Mutations In Cancer. Nucleic Acids Res. 2019, 47, D941–D947. [Google Scholar] [CrossRef]
  31. Coussens, L.M.; Werb, Z. Inflammation and cancer. Nature 2002, 420, 860. [Google Scholar] [CrossRef] [PubMed]
  32. Bach, D.-H.; Lee, S.K. Long noncoding RNAs in cancer cells. Cancer Lett. 2018, 419, 152–166. [Google Scholar] [CrossRef] [PubMed]
  33. Cavallo, F.; De Giovanni, C.; Nanni, P.; Forni, G.; Lollini, P.L. 2011: The immune hallmarks of cancer. Cancer Immunol. Immunother. CII 2011, 60, 319–326. [Google Scholar] [CrossRef]
  34. Hanahan, D.; Weinberg, R.A. Hallmarks of cancer: The next generation. Cell 2011, 144, 646–674. [Google Scholar] [CrossRef] [PubMed]
  35. Cairns, R.A.; Mak, T.W. The current state of cancer metabolism. Nat. Rev. Cancer 2016, 16, 613. [Google Scholar] [CrossRef]
  36. Phan, L.M.; Yeung, S.C.; Lee, M.H. Cancer metabolic reprogramming: Importance, main features, and potentials for precise targeted anti-cancer therapies. Cancer Biol. Med. 2014, 11, 1–19. [Google Scholar] [CrossRef]
  37. Muñoz-Pinedo, C.; El Mjiyad, N.; Ricci, J.E. Cancer metabolism: Current perspectives and future directions. Cell Death Dis. 2012, 3, e248. [Google Scholar] [CrossRef] [PubMed]
  38. Bach, D.-H.; Hong, J.-Y.; Park, H.J.; Lee, S.K. The role of exosomes and miRNAs in drug-resistance of cancer cells. Int. J. Cancer 2017, 141, 220–230. [Google Scholar] [CrossRef]
  39. Bach, D.-H.; Park, H.J.; Lee, S.K. The Dual Role of Bone Morphogenetic Proteins in Cancer. Mol. Ther. Oncol. 2018, 8, 1–13. [Google Scholar] [CrossRef]
  40. Hamada, S.; Satoh, K.; Fujibuchi, W.; Hirota, M.; Kanno, A.; Unno, J.; Masamune, A.; Kikuta, K.; Kume, K.; Shimosegawa, T. MiR-126 acts as a tumor suppressor in pancreatic cancer cells via the regulation of ADAM9. Mol. Cancer Res. MCR 2012, 10, 3–10. [Google Scholar] [CrossRef] [PubMed]
  41. Keklikoglou, I.; Hosaka, K.; Bender, C.; Bott, A.; Koerner, C.; Mitra, D.; Will, R.; Woerner, A.; Muenstermann, E.; Wilhelm, H.; et al. MicroRNA-206 functions as a pleiotropic modulator of cell proliferation, invasion and lymphangiogenesis in pancreatic adenocarcinoma by targeting ANXA2 and KRAS genes. Oncogene 2015, 34, 4867–4878. [Google Scholar] [CrossRef] [PubMed]
  42. Malgerud, L.; Lindberg, J.; Wirta, V.; Gustafsson-Liljefors, M.; Karimi, M.; Moro, C.F.; Stecker, K.; Picker, A.; Huelsewig, C.; Stein, M.; et al. Bioinformatory-assisted analysis of next-generation sequencing data for precision medicine in pancreatic cancer. Mol. Oncol. 2017, 11, 1413–1429. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef]
  44. Gautier, L.; Cope, L.; Bolstad, B.M.; Irizarry, R.A. Affy—Analysis of Affymetrix GeneChip data at the probe level. Bioinformatics 2004, 20, 307–315. [Google Scholar] [CrossRef] [PubMed]
  45. Xia, J.; Gill, E.E.; Hancock, R.E.W. NetworkAnalyst for statistical, visual and network-based meta-analysis of gene expression data. Nat. Protoc. 2015, 10, 823. [Google Scholar] [CrossRef]
  46. Long, N.P.; Jung, K.H.; Yoon, S.J.; Anh, N.H.; Nghi, T.D.; Kang, Y.P.; Yan, H.H.; Min, J.E.; Hong, S.S.; Kwon, S.W. Systematic assessment of cervical cancer initiation and progression uncovers genetic panels for deep learning-based early diagnosis and proposes novel diagnostic and prognostic biomarkers. Oncotarget 2017, 8, 109436–109456. [Google Scholar] [CrossRef]
  47. Lê, S.; Josse, J.; Husson, F. FactoMineR: An R Package for Multivariate Analysis. J. Stat. Softw. 2008, 25, 18. [Google Scholar] [CrossRef]
  48. Xia, J.; Sinelnikov, I.V.; Han, B.; Wishart, D.S. MetaboAnalyst 3.0—Making metabolomics more meaningful. Nucleic Acids Res. 2015, 43, W251–W257. [Google Scholar] [CrossRef]
  49. Kuhn, M. Building Predictive Models in R Using the caret Package. J. Stat. Softw. 2008, 28, 26. [Google Scholar] [CrossRef]
  50. Pedersen, T.L.; Benesty, M. Lime: Local Interpretable Model-Agnostic Explanations. R Package Version 0.4.1. 2018. Available online: https://CRAN.R-project.org/package=lime (accessed on 27 January 2019).
  51. Tang, Z.; Li, C.; Kang, B.; Gao, G.; Li, C.; Zhang, Z. GEPIA: A web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017, 45, W98–W102. [Google Scholar] [CrossRef]
  52. Anaya, J. OncoLnc: Linking TCGA survival data to mRNAs, miRNAs, and lncRNAs. PeerJ Comput. Sci. 2016, 2, e67. [Google Scholar] [CrossRef]
  53. Tan, Y.; Tan, Y.; Lu, L.; Zhang, H.; Sun, C.; Liang, Y.; Zeng, J.; Yang, X.; Li, D.; Zou, H. HPCDb: An integrated database of pancreatic cancer. bioRxiv 2017. [Google Scholar] [CrossRef]
  54. Dayem Ullah, A.Z.; Cutts, R.J.; Ghetia, M.; Gadaleta, E.; Hahn, S.A.; Crnogorac-Jurcevic, T.; Lemoine, N.R.; Chelala, C. The pancreatic expression database: Recent extensions and updates. Nucleic Acids Res. 2014, 42, D944–D949. [Google Scholar] [CrossRef] [PubMed]
  55. Baker, S.; Ali, I.; Silins, I.; Pyysalo, S.; Guo, Y.; Hogberg, J.; Stenius, U.; Korhonen, A. Cancer Hallmarks Analytics Tool (CHAT): A text mining approach to organize and evaluate scientific literature on cancer. Bioinformatics 2017, 33, 3973–3981. [Google Scholar] [CrossRef] [PubMed]
  56. Chou, C.H.; Shrestha, S.; Yang, C.D.; Chang, N.W.; Lin, Y.L.; Liao, K.W.; Huang, W.C.; Sun, T.H.; Tu, S.J.; Lee, W.H.; et al. miRTarBase update 2018: A resource for experimentally validated microRNA-target interactions. Nucleic Acids Res. 2018, 46, D296–D302. [Google Scholar] [CrossRef] [PubMed]
  57. Cotto, K.C.; Wagner, A.H.; Feng, Y.-Y.; Kiwala, S.; Coffman, A.C.; Spies, G.; Wollam, A.; Spies, N.C.; Griffith, O.L.; Griffith, M. DGIdb 3.0: A redesign and expansion of the drug–gene interaction database. Nucleic Acids Res. 2018, 46, D1068–D1073. [Google Scholar] [CrossRef] [PubMed]
  58. Lee, D.-K.; Long, N.P.; Jung, J.; Kim, T.J.; Na, E.; Kang, Y.P.; Kwon, S.W.; Jang, J. Integrative lipidomic and transcriptomic analysis of X-linked adrenoleukodystrophy reveals distinct lipidome signatures between adrenomyeloneuropathy and childhood cerebral adrenoleukodystrophy. Biochem. Biophys. Res. Commun. 2018. [Google Scholar] [CrossRef]
  59. Long, N.P.; Park, S.; Anh, N.H.; Min, J.E.; Yoon, S.J.; Kim, H.M.; Nghi, T.D.; Lim, D.K.; Park, J.H.; Lim, J.; et al. Efficacy of Integrating a Novel 16-Gene Biomarker Panel and Intelligence Classifiers for Differential Diagnosis of Rheumatoid Arthritis and Osteoarthritis. J. Clin. Med. 2019, 8, 50. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Screening process of innovative biomarker candidates for pancreatic cancer diagnosis and treatment. (a) Combined effect size distributions in transcriptome meta-analysis. (b) Selection of membrane protein-coding genes. (c) Selection of secretory protein-coding genes. The left circle of the Venn diagram (orange) represents the candidates from meta-analysis and the right circle of the Venn diagram (blue) represents the candidates curated from [17].
Figure 1. Screening process of innovative biomarker candidates for pancreatic cancer diagnosis and treatment. (a) Combined effect size distributions in transcriptome meta-analysis. (b) Selection of membrane protein-coding genes. (c) Selection of secretory protein-coding genes. The left circle of the Venn diagram (orange) represents the candidates from meta-analysis and the right circle of the Venn diagram (blue) represents the candidates curated from [17].
Cancers 11 00155 g001
Figure 2. Kaplan–Meier (KM) plots of the overall survival of four promising prognostic candidates. (a) KM plot of ADAM9. (b) KM plot of ANXA2. (c) KM plot of ITGA2. (d) KM plot of MET.
Figure 2. Kaplan–Meier (KM) plots of the overall survival of four promising prognostic candidates. (a) KM plot of ADAM9. (b) KM plot of ANXA2. (c) KM plot of ITGA2. (d) KM plot of MET.
Cancers 11 00155 g002
Figure 3. Data exploration and diagnostic performance of 11 biomarker candidates in the Random Forests model. (a) The gene expression of CDH3 is higher in PDAC than in normal controls (n = 96), ****: p < 0.0001. (b) The gene expression of LAMC2 is higher in PDAC than in normal controls (n = 96). (c) Principal component analysis of PDAC versus normal controls. (d) Heatmap analysis of 11 biomarker candidates. (e) ROC curve of the random forests model in the test set (n = 28). (f) The importance scores of 11 biomarker candidates in the random forests model.
Figure 3. Data exploration and diagnostic performance of 11 biomarker candidates in the Random Forests model. (a) The gene expression of CDH3 is higher in PDAC than in normal controls (n = 96), ****: p < 0.0001. (b) The gene expression of LAMC2 is higher in PDAC than in normal controls (n = 96). (c) Principal component analysis of PDAC versus normal controls. (d) Heatmap analysis of 11 biomarker candidates. (e) ROC curve of the random forests model in the test set (n = 28). (f) The importance scores of 11 biomarker candidates in the random forests model.
Cancers 11 00155 g003
Figure 4. Immunohistochemical analysis for determination of LAMC2, ADAM9, and ANXA2 (Annexin A2) gene expression. (a) Scoring system for the three proteins in pancreatic tissues. (b) IHC scores of LAMC2, ADAM9, and ANXA2 in pancreatic tumors and normal controls (n = 86). ***: p < 0.001; ****: p < 0.0001. (c) Pairwise correlation of LAMC2, ADAM9, and ANXA2 expression. (d) Pairwise correlation of LAMC2, ADAM9, and ANXA9 at gene expression level (log2 of transcripts per million) in TGCA PC, TCGA normal pancreas, and Genotype-Tissue Expression Project (GTEx) derived normal pancreas.
Figure 4. Immunohistochemical analysis for determination of LAMC2, ADAM9, and ANXA2 (Annexin A2) gene expression. (a) Scoring system for the three proteins in pancreatic tissues. (b) IHC scores of LAMC2, ADAM9, and ANXA2 in pancreatic tumors and normal controls (n = 86). ***: p < 0.001; ****: p < 0.0001. (c) Pairwise correlation of LAMC2, ADAM9, and ANXA2 expression. (d) Pairwise correlation of LAMC2, ADAM9, and ANXA9 at gene expression level (log2 of transcripts per million) in TGCA PC, TCGA normal pancreas, and Genotype-Tissue Expression Project (GTEx) derived normal pancreas.
Cancers 11 00155 g004aCancers 11 00155 g004b
Table 1. Literature- and experiment-based validation of the 23 candidates.
Table 1. Literature- and experiment-based validation of the 23 candidates.
Gene SymbolEntrez IDRNA AlterationProtein AlterationNGS ResultsMeta-AnalysisNeoplasm # vs. NormalcyPDAC vs. Neoplasm #
Log FCp-ValueFDRcESp-Value
AGR210551↑RT-PCR 1,2, ↑SAGE 1,2, ↑DM 1↑IHC 1, ↑LFQ-MS 1, ↑WB 23.671.55 × 10−157.37 × 10−131.930X
GAPDH2597↑DM 1, ↑SAGE 1, NB 2↑ICAT 1, ↑WB 1,2, ↓SILAC-TMS 21.486.08 × 10−41.27 × 10−21.682.02 × 10−13XX
LAMC23918↑DM 1,2↑IHC 13.911.15 × 10−71.17 × 10−52.370X
MMP114320↑DM 1, ↑ISH 1, ↑SAGE 1,2, ↑NB 1,2↑IHC 1, ↑WB 2NANANA1.627.90 × 10−14XX
TAGLN28407NA↑SILAC-TMS 1, ↑ICAT 1NANANA1.731.64 × 10−5X
ADAM98754↑DM 1,2, ↑RT-PCR 1,2↑IHC 1NANANA1.711.41 × 10−5XX
ANXA2302↑DM 1↑IHC 1, ↑ICAT 1, ↑WB 1,2, ↑MS 1, ↑TDE 1, ↑LFQ-MS 11.162.26 × 10−33.35 × 10−21.638.33 × 10−9XX
APLP2334↑DM 1,2↑SILAC-TMS 1,2, ↑SILAC 2NANANA1.531.11 × 10−5XX
CDH31001↑DM 1,2, ↑RT-PCR 1,2↑IHC 12.728.70 × 10−64.50 × 10−42.020X
MSLN10232↓DM 1, ↑DM 1,2, ↑RT-PCR 1,2, ↑SAGE 1,2, ↑ISH 1↑IHC 15.211.89 × 10−158.77 × 10−132.034.88 × 10−4X
SERPINB55268↑DM 1,2, ↑NB 2↑IHC 1,2, ↑SILAC-TMS 18.356.54 × 10−302.96 × 10−261.960X
CD823732↑NB 1, ↑ISH 1, ↑DM 1↑IHC 11.394.55 × 10−41.03 × 10−21.627.90 × 10−14XX
CLDN1851208↑DM 1↑IHC 12.579.43 × 10−52.97 × 10−31.647.60 × 10−9X
EPHA21969↑DM 1,2↑SILAC-TMS 1, ↑WB 21.694.95 × 10−51.78 × 10−31.562.99 × 10−13X
EZR7430NA↑IHC 1, SILAC-TMS 1NANANA1.680X
FXYD35349↑DM 1,2, ↑SAGE 1, ↑NB 1, ↑ISH 1NA4.714.03 × 10−215.21 × 10−181.910X
GPRC5A9052↑DM 1,2, ↑SAGE 1,2↑IHC 13.982.08 × 10−82.61 × 10−61.911.25 × 10−8
ITGA23673↑DM 1,2↑IHC 1, ↑WB 12.941.20 × 10−68.72 × 10−52.042.02 × 10−13X
ITGB63694↑DM 1↑IHC 12.564.55 × 10−62.71 × 10−41.792.78 × 10−8XX
MET4233↑DM 1, ↑RT-PCR 1,2NANANANA1.901.97 × 10−7XX
MST1R4486↑DM 1,2, ↑RT-PCR 2↑IHC 1, ↑WB 22.849.03 × 10−76.81 × 10−52.060X
NQO11728↑DM 1,2↑IHC 1, ↑WB 23.711.88 × 10−93.04 × 10−72.151.63 × 10−5X
SLC2A16513↑DM 1,2↑IHC 13.286.34 × 10−101.16 × 10−72.030
# Data obtained from the dataset GSE19650. 1 Evidence in PDAC tissue. 2 Evidence in cancer cell lines. Abbreviations: cES: combined effect size; DM: DNA microarray; FDR: false discovery rate; ICAT: isotope-coded affinity tag; IHC: immunohistochemistry; IPMN: intraductal papillary mucinous neoplasia; ISH: in situ hybridization; LFQ-MS: label-free quantitation mass spectrometry; NB: northern blot; NGS: next-generation sequencing; PDAC: pancreatic ductal adenocarcinoma; RT-PCR: reverse transcription-polymerase chain reaction; SAGE: serial analysis of gene expression; SILAC-TMS: stable isotope labeling by/with amino acids in cell culture-tandem mass spectrometry; TDE: two dimensional electrophoresis; WB: western blot; NA: Not applicable/available.
Table 2. Combined effect sizes, Cox regression analysis, and alterations at the mRNA and protein level of the 23 candidates.
Table 2. Combined effect sizes, Cox regression analysis, and alterations at the mRNA and protein level of the 23 candidates.
Gene SymbolEntrez IDCox Regression
Cox CoefficientFDR
AGR2105510.050.76
GAPDH25970.210.19
LAMC239180.420.02
MMP1143200.170.28
TAGLN284070.190.22
ADAM987540.410.02
ANXA23020.320.04
APLP23340.230.14
CDH310010.370.03
MSLN102320.260.10
SERPINB552680.380.02
CD8237320.110.51
CLDN18512080.090.63
EPHA219690.330.05
EZR74300.280.08
FXYD353490.170.30
GPRC5A90520.370.03
ITGA236730.410.02
ITGB636940.490.01
MET42330.660.00
MST1R44860.220.17
NQO117280.140.43
SLC2A165130.270.08
Table 3. Drug–gene and miRNA–gene interactions of 18 membrane proteins.
Table 3. Drug–gene and miRNA–gene interactions of 18 membrane proteins.
Gene SymbolmiRNA–Gene InteractionDruggable GenomeDrug–Gene Interaction 3
miRNA 1 That Targets the GeneValidation MethodsExpression Profile in PC (Correlation Coefficient 2)
Strong EvidenceLess Strong Evidence
ADAM9hsa-miR-126-3pRA, WB, qPCRMA0.315YesIlomastat
hsa-miR-33a-5pRA, WB, qPCR 0.373
hsa-miR-125a-5pqPCR 0.320
ANXA2hsa-miR-155-5pRA, WB, qPCRMA−0.599YesNA
hsa-miR-206RA, WB, qPCR −0.312
APLP2NA YesNA
CDH3NA YesNA
MSLNhsa-miR-21-5pRA 0.715YesAmatuximab
SERPINB5hsa-miR-21-5pRA, WB, qPCR MA, NGS0.709YesNA
hsa-miR-103a-3pRA, WB, qPCR MA, NGS0.501
CD82NA YesNA
CLDN18NA YesClaudiximab
EPHA2NA YesDasatinib, Dorsomorphin, Regorafenib, Vandetanib
EZRhsa-miR-183-5pRA, WB, qPCR 0.424NANA
hsa-miR-204-5pRA, WB, qPCRMA−0.542
hsa-miR-205-5pRAMA, NGS0.358
FXYD3NA NoNA
GPRC5Ahsa-miR-103a-3pRA, WB, qPCRNGS0.514YesNA
ITGA2hsa-miR-16-5pqPCRNGS, pSILAC0.509YesAbciximab, CHEMBL36326, Eptifibatide, Tirofiban, Vatelizumab
ITGB6NA YesDI17E6, Intetumumab, STX-100
METhsa-miR-34c-5pRA, WB, qPCRMA, NGS0.337YesARRY-300, ABT-700, AMG-337, AMG-208, Amuvatinib, Altiratinib, Amoxicillin, Alectinib, BMS-698769, BMS-777607, BMS-794833, BMS-817378, BPI-9016, Crizotinib, Clofibrate, Cabozantinib, Cabozantinib S-Malate, Capmatinib, Crenolanib, Emibetuzumab, EMD-1204831, Foretinib, Golvatinib, JNJ-38877605, LY-2875358, MK-8033, MGCD-265, Merestinib, MK-2461, Onartuzumab, PF-04217903, PHA-665752, Pyrazinamide, Rilonacept, SGX-523, Savolitinib, Tepotinib, Tivantinib, Tanespimycin, SAR-125844, TAS-115
hsa-miR-199a-3pRA, WB, qPCRMA0.365
hsa-miR-34a-5pRA, WB, qPCR 0.478
hsa-miR-23b-3pRA, WB, qPCR 0.447
hsa-miR-27a-3pRA, WB, qPCR 0.444
hsa-miR-27b-3pRA, WB, qPCR 0.458
hsa-miR-31-5pRA, qPCR 0.594
hsa-miR-34a-3pWB 0.380
MST1RNA YesBMS-777607, Foretinib, MGCD-265, MK-2461, MK-8033
NQO1NA YesApaziquone, Dicumarol, Vatiquinone
SLC2A1hsa-miR-22-3pRA, WB, qPCRNGSs0.331YesNA
1 Only miRNAs that have strong evidence were included, 2 p-value of the correlation coefficient was less than 0.05. 3 Only defined interactions were included. Abbreviation: MA: microarray; NGS: next-generation sequencing; pSILAC: pulsed stable isotope labeling by amino acids in cell culture; qPCR: quantitative polymerase chain reaction; RA: reporter assay; WB: western blot.

Share and Cite

MDPI and ACS Style

Long, N.P.; Jung, K.H.; Anh, N.H.; Yan, H.H.; Nghi, T.D.; Park, S.; Yoon, S.J.; Min, J.E.; Kim, H.M.; Lim, J.H.; et al. An Integrative Data Mining and Omics-Based Translational Model for the Identification and Validation of Oncogenic Biomarkers of Pancreatic Cancer. Cancers 2019, 11, 155. https://doi.org/10.3390/cancers11020155

AMA Style

Long NP, Jung KH, Anh NH, Yan HH, Nghi TD, Park S, Yoon SJ, Min JE, Kim HM, Lim JH, et al. An Integrative Data Mining and Omics-Based Translational Model for the Identification and Validation of Oncogenic Biomarkers of Pancreatic Cancer. Cancers. 2019; 11(2):155. https://doi.org/10.3390/cancers11020155

Chicago/Turabian Style

Long, Nguyen Phuoc, Kyung Hee Jung, Nguyen Hoang Anh, Hong Hua Yan, Tran Diem Nghi, Seongoh Park, Sang Jun Yoon, Jung Eun Min, Hyung Min Kim, Joo Han Lim, and et al. 2019. "An Integrative Data Mining and Omics-Based Translational Model for the Identification and Validation of Oncogenic Biomarkers of Pancreatic Cancer" Cancers 11, no. 2: 155. https://doi.org/10.3390/cancers11020155

APA Style

Long, N. P., Jung, K. H., Anh, N. H., Yan, H. H., Nghi, T. D., Park, S., Yoon, S. J., Min, J. E., Kim, H. M., Lim, J. H., Kim, J. M., Lim, J., Lee, S., Hong, S. -S., & Kwon, S. W. (2019). An Integrative Data Mining and Omics-Based Translational Model for the Identification and Validation of Oncogenic Biomarkers of Pancreatic Cancer. Cancers, 11(2), 155. https://doi.org/10.3390/cancers11020155

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop