Basic Study Open Access
Copyright ©The Author(s) 2019. Published by Baishideng Publishing Group Inc. All rights reserved.
World J Gastroenterol. Aug 28, 2019; 25(32): 4727-4738
Published online Aug 28, 2019. doi: 10.3748/wjg.v25.i32.4727
Integrative analysis of the inverse expression patterns in pancreas development and cancer progression
Hong-Liang Zang, Guo-Min Huang, Xiao-Feng Tian, Department of Hepatobiliary and Pancreatic Surgery, China-Japan Union Hospital of Jilin University, Changchun 130033, Jilin Province, China
Hai-Ying Ju, Department of Hematology, Jilin Province Blood Center, Changchun 130000, Jilin Province, China
ORCID number: Hong-Liang Zang (0000-0002-9570-6158); Guo-Min Huang (0000-0002-7473-7780); Hai-Ying Ju (0000-0001-6428-803X); Xiao-Feng Tian (0000-0001-7146-5313).
Author contributions: Zang HL performed the majority of experiments and analyzed the data; Huang GM performed the molecular investigations; Ju HY designed and coordinated the research; Tian XF wrote the paper.
Conflict-of-interest statement: The authors declare no conflict of interest.
Data sharing statement: No additional data are available.
Open-Access: This article is an open-access article which was selected by an in-house editor and fully peer-reviewed by external reviewers. It is distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/
Corresponding author: Xiao-Feng Tian, MSc, Attending Doctor, Department of Hepatobiliary and Pancreatic Surgery, China-Japan Union Hospital of Jilin University, No.126 Sendai Street, Changchun 130033, Jilin Province, China. chuie03755@163.com
Telephone: +86-431-84995999
Received: May 5, 2019
Peer-review started: May 5, 2019
First decision: June 10, 2019
Revised: July 5, 2019
Accepted: July 19, 2019
Article in press: July 19, 2019
Published online: August 28, 2019

Abstract
BACKGROUND

As the malignant tumor, pancreatic cancer with a meager 5-years survival rate has been widely concerning. However, the molecular mechanisms that result in malignant transformation of pancreatic cells remain elusive.

AIM

To investigate the gene expression profiles in normal or malignant transformed pancreas development.

METHODS

MaSigPro and ANOVA were performed on two pancreas development datasets downloaded from the Gene Expression Omnibus database. Six pancreatic cancer datasets collected from TCGA database were used to establish differentially expressed genes related to pancreas development and pancreatic cancer. Moreover, gene clusters with highly similar interpretation patterns between pancreas development and pancreatic cancer progression were established by self-organizing map and singular value decomposition. Additionally, the hypergeometric test was performed to compare the corresponding interpretation patterns. Abnormal regions of metabolic pathway were analyzed using the Sub-pathway-GM method.

RESULTS

This study established the continuously upregulated and downregulated genes at different stages in pancreas development and progression of pancreatic cancer. Through analysis of the differentially expressed genes, we established the inverse and consistent direction development-cancer pattern associations. Based on the application of the Subpathway-GM analysis, we established 17 significant metabolic sub-pathways that were closely associated with pancreatic cancer. Of note, the most significant metabolites sub-pathway was related to glycerophospholipid metabolism.

CONCLUSION

The inverse and consistent direction development-cancer pattern associations were established. There was a significant correlation in the inverse patterns, but not consistent direction patterns.

Key Words: Pancreatic cancer, Pancreas development, Inverse pattern, Metabolites subpathway

Core tip: Here we analyzed the differentially expressed genes during pancreas development and pancreatic cancer progression. We noticed that genes that up-regulated in tumorigenesis were conversely suppressed in the development of pancreas. Reciprocally, upregulated gene expression pattern during pancreas development were negatively correlated with pancreatic cancer progression. Additionally, 17 significant metabolic sub-pathways, especially glycerophospholipid metabolism, were identified, which is highly correlated with pancreas cancer development.



INTRODUCTION

As one of the most malignant tumors, pancreatic cancer results in more than 0.4 million deaths per year[1]. Generally, it is used to refer to pancreatic adenocarcinoma, which is the most common and devastating in all types of pancreatic cancer. According to the Union for International Cancer Control classification criteria, patients with pancreatic cancer can be divided into four main stages, including stage I, II, III and IV[2]. The overall five-year survival rate of the pancreas is less than 5% and the median survival period is six months after diagnosis, which is the lowest survival rate among all types of cancers[3,4]. Early diagnosis of biomarkers and effective treatments can help us effectively conduct pancreatic cancer research.

Similar to the tumor progression, organ development can also be divided into some stages in a time-dependent manner[5]. Moreover, the development is strictly controlled by multiple signaling pathways and transcription factors[6]. Some researchers believe that cancer is the problem of developmental biology. It has been reported that serval genes and pathways affected the normal or malignant transformed pancreas development[7]. Pancreatic and duodenal homeobox 1 (PDX1), which is exclusively expressed in the pancreas, is essential for pancreas development. Recently research shows that dysfunction of PDX1 also promotes pancreatic cancer development and progression[8,9].

We therefore investigated the relationship between pancreas development and pancreatic cancer progression by analyzing two datasets related with pancreas development and six datasets related with pancreatic cancer. Through bioinformatics analysis, we established differentially gene expression profiles and multiple patterns that were consistent/inverse in the development-cancer patterns. Moreover, we found that there was significant negative correlation between inverse development and cancer. Of note, we identified 17 metabolic sub-pathways that were highly related with pancreatic cancer development.

MATERIALS AND METHODS
Pancreas development datasets

To further understand the pancreas development, we analyzed two databases (accession: GSE42094, GSE96697), which obtained from the Gene Expression Omnibus (GEO) database. 16 samples in 6 development stages data were obtained from GSE42094, including undifferentiated human embryonic stem cells, stage 1 (S1), S2, S3, S4, and S5. In the GSE96697 database, according to the glycoprotein 2 (GP2) and CDH1 interpretation, we classified the 7-weeks pancreas development data into three stages, including an early, middle and late stage. We established that high interpretation of GP2 was early stage, while low interpretation was the middle stage. Additionally, CDH1 low interpretation of CDH1 was established as a late stage.

Pancreatic cancer datasets

To further understand the development and progression of pancreatic cancer, we analyzed five pancreatic cancer datasets which obtained from TCGA (https://portal.gdc.cancer.gov/), ICGC (https://icgc.org/) and GEO database[10-12]. All of these datasets contained five or six clinical stages and available clinical information. The accession numbers in GEO database were GSE62452, GSE79668, and GSE102238. Detail information of the datasets was listed in Table 1.

Table 1 Information of the pancreatic cancer datasets.
Dataset nameTCGAGSE62452PACA-AU (ICGC)PACA-CA (ICGC)GSE79670GSE102238
Histological typePDACPDACPDACPDACPDACPDAC
Samples14569892025048
Gender Female, Male68, 7731, 3841, 4894, 10819, 3121, 27
StageI (IA, IB)3, 90, 40, 60, 113, 60, 8
StageII (IIA, IIB)23, 10410, 3620, 4728, 1554, 3119, 15
StageIII3139563
StageIV367303
PlatformRNA-seqAffymetrixRNA-seqRNA-seqRNA-seqAgilent
Metabolite highly related to PACA

According to the previous studies, we established the differentially abundant metabolites and converted them to KEGG compound IDs[13-16]. These metabolites were obtained from the differentially abundant between control and pancreatic cancer. A total of 60 pancreatic cancer-associated differentially abundant metabolites were established.

Identification of differentially expressed genes in time-course datasets

We found the gene expression profile with significant differences between experimental groups in time course datasets using maSigPro approach. MaSigPro were applied to establish experimental groups using the dummy variables with a two-regression step, and it is an R package for significance analysis time-course microarray experiments. Briefly, after adjust, a global-scale regression model with all the defined variables, differently expressed genes was established. We applied a variable selection strategy to investigate the virtual difference between groups to establish the different profiles. According to the patient clinical stages, we could establish the pancreatic cancer interpretation profiles as “time-series” datasets. We selected 0.05 as false discovery rate (FDR) significance threshold for the analysis in this study without a special request. In GSE96697 dataset, due to the fewer development stage, we used ANOVA method to get differentially expressed genes (DEGs).

Identification of gene cluster with consistent expression patterns

To the investigation of the pancreas development and pancreatic cancer development and progression associated interpretation patterns, we performed an initial screening on the interpretation matrix of DEGs. We pre-processed the DEGs using a self-organizing map (SOM) strategy and performed the pattern recognition using the singular value decomposition (SVD) strategy. We conducted a comprehensive SOM-SVD analysis under the guidance of as previous described (paper). Briefly, the analysis could conclude into three steps. Initially, we conducted the SOM transformation. Secondly, SOM output was decomposed using SVD strategy. FDR, which is used for significant neuron assessment, is set at 0.10 to select the pancreas development and pancreatic cancer progression pattern genes. The identified genes were analyzed using component plane presentation-integrated SOM for gene clustering. The highly similar expression pattern clusters were used for pancreas development and pancreatic cancer progression integrative analysis.

Function for Geneset enrichment analysis

In our study, we annotated the identified interpretation patterns using the Database for Annotation, Visualization and Integrated Discovery, which are comprehensive sets of functional annotation tools. The identified interpretation patterns were clustered into more than 40 Gene Ontology-Biological Process (GO-BP) terms. R package iSubpathwayMiner was applied for subpathway enrichment analysis.

Comparison of gene expression pattern between development and cancer

We identified the up-regulated interpretation patterns and down-regulated interpretation patterns in the pancreas development as dev-Up and dev-Dw. To analyze the associations between pancreas development and pancreatic progression, we use the hypergeometric test in our study, as previously described[17]. Briefly, we compared the dev-Up vs can-Up, dev-Up vs can-Dw, dev-Dw vs can-Up and dev-Dw vs can-Up to generated the development and cancer dataset pairs. Moreover, we considered that P < 0.05 was the significant correlation between pancreas development and pancreatic cancer progression. The p-value was calculated according to the following formula:

Figure formula
Figure formula 

In this formula, N-development and N-cancer indicate the number of pancreas development and pancreatic cancer regulation genes. M indicates the whole genome number, while k indicates the number of the common gene.

RESULTS
Identification of genes differentially expressed in pancreas development

To investigate the transcriptional pattern during pancreas development, we interrogated the interpretation patterns in GSE42094, GSE96697, which are the datasets of pancreas development. Using the maSigPro method, we identified 3069 DEGs at different time points of pancreas development in GSE42094 dataset. We then employed SOM-SVD strategy to select the topology-preserving DEGs, according to the interpretation matrices (Figure 1A). Reciprocally, the selected genes in a topology-preserving selection further confirmed the alteration in “time-series” processes. The entire genes were automatically selected in this method. Then, the matrices with 1257 genes, which obtained from SOM-SVD analysis, were clustered into four gene clusters (cluster 1-4). As shown in Figure 1A, the cluster 2 and cluster 4 contained genes that were transiently upregulated in the early stage of development, then decreased gradually along with development. We therefore named these gene sets as continuous down-regulated expression patterns. In contrast, genes with low expression level and were increased gradually in the latter stage of development, which were observed in cluster 1 and cluster 3. They were thus identified as continuously up-regulated expression patterns.

Figure 1
Figure 1 Continuous differential expression patterns in pancreas development. A: Clustering of human pancreas development genes. In the heat map, green indicates down-regulated, and red indicates up-regulated. In the line graphs, lines represent the tendency of the cluster changes; B: Gene Ontology-Biological Process annotations of the continuously up- and down-regulated genes.

For the GSE96697 data set, we used the ANOVA method to analyze the expression data set and identified 3078 differential expressed genes. The K-means clustering method was used to establish patterns of interpretation of DEG sets, and we established six clusters. Similar to the previous continuous adjustment pattern recognition, we found that cluster 4 and cluster 6 are the interpretation modes of down-regulation, and cluster 2 is the interpretation mode of up-regulation (Figure 2). From this analysis, we found that 641 and 616 genes in the GSE42094 and GSE96697 data sets were continuously up-regulated, while 1059 and 1052 genes were continuously down-regulated. To investigate the biological characteristics of genes related with pancreatic development, functional enrichment analysis was used to aggregate genes in a consistent up- or down-regulation model. We found that ''lipid digestion' and ''cholesterol homeostasis'' were enriched during pancreas development and continued to up-regulate the mode of interpretation. "Cell proliferation," "mitotic nuclear division" is annotated in an interpretive pattern that is continuously down-regulated (Figure 1B).

Figure 2
Figure 2 The expression patterns identified from GSE96697. Among these results, up- and downregulated patterns (red, green color labeled) were defined and further analyzed in our study.
Exploration of genes that continuously regulated in pancreatic cancer progression

To measure the status of genes related with pancreas development in the progression of pancreatic cancer, we analyzed six pancreatic cancer data sets described in Materials and Methods. We considered each clinical stage of pancreatic cancer as the point in time to determine progression patterns, similar to pancreatic developmental analysis. Thus, we can determine the various interpretation patterns in these pancreatic cancer data sets (Figure 3A). Comparative analysis was used to study the mode of interpretation of continuous regulation between pancreatic development and pancreatic cancer progression. We established the interpretation mode of the up-regulation of the tumor and the interpretation mode of the down-regulation according to the following criteria. Exceeding [(n-1)/ 2 + 1] the interpretation level of the adjacent stage changes with the same trend, n repeats the number of stages in each data set. The absolute slope is more significant than 0.05. We established six gene clusters in the up-regulated pattern and six of the down-regulated patterns (Figure 3B).

Figure 3
Figure 3 The expression patterns identified from six pancreatic cancer datasets. A: Among these results, up- and downregulated patterns (red, blue color labeled) were defined and further analyzed in our study; B: Pie diagram showed the up- and down-regulation genes in each pancreatic cancer dataset.
Comparative analysis of gene expression between pancreas development and pancreatic cancer progression

To investigate the relationship between pancreatic development and pancreatic cancer progression, we performed a hypergeometric test. As shown in Figure 4A, retro-pancreatic development and pancreatic cancer patterns, including dev-Up vs can-Dw and dev-Dw vs can-Up. On the other hand, we found a weak relationship between inconsistent development and cancer patterns. Specifically, we did not find any significant correlation between any two dev-Dw and can-Dw datasets. The inverse interpretation patterns were clustered into GO-BP terms to analyze the biological function, including dev-Up vs can-Dw and Dev-Dw vs can-Up patterns (Figure 4B). The results showed that 141 genes with dev-Up vs can-Dw were mainly associated with immune-related BP terms, including ‘‘T Cell Proliferation” and ‘‘Innate Immune Response In Mucosa’’. Furthermore, 202 genes with dev-Dw vs can-Up were exclusively involved in proliferation-related BP terms, including “Cell Division” and ‘‘DNA Replication Initiation’’. Collectively, cell proliferation activity, as one of the essential characteristics in the malignant tumor, was gradually enhanced along with cancer progression, which was consistent with previous studies.

Figure 4
Figure 4 Integrated analysis of expression patterns in pancreas development and pancreatic cancer progression. A: Comparison of interpretation patterns between pancreas development and cancer progression. Color in each cell indicated the p-value; B: Gene Ontology-Biological Process annotations of the inverse interpretation patterns, including dev-up vs can-dw and dev-dw vs can-up.
Identification of Metabolic Sub-pathways Associated with pancreatic cancer

A total of 343 genes were established in the inverse interpretation patterns. To further establish pancreatic cancer associated with metabolic subpathways, we found 60 unique differentially abundant metabolites, which might be associated with pancreatic cancer progression. After integrating the analysis of 343 genes and 60 metabolites, Sub-pathway-GM strategy was employed to establish the critical, abnormal regions were identified in each metabolic pathway. Subsequently, we set FDR < 0.01 as a threshold for further analysis of 343 genes and 60 differential metabolite pathways and identified 17 significant metabolic sub-avenues (Table 2). Among these established sub-pathways, the most significant was "Glycerophospholipid metabolism" (path:00564_1) which was critical for lipid metabolism. Our data thus demonstrate that activation of metabolite pathways, especially lipid metabolism, are crucial for pancreatic cancer development.

Table 2 17 significant metabolic subpathways in pancreatic cancer.
Pathway IDPathway nameAnn molecule ratioAnn bg ratioP valueFDR
path:00564_1Glycerophospholipid metabolism9/37847/250513.03E-085.15E-07
path:00330_1Arginine and proline metabolism5/37818/250515.55E-064.72E-05
path:00020_1Citrate cycle (TCA cycle)3/3789/250510.0002676130.001516472
path:00630_1Glyoxylate and dicarboxylate metabolism3/37811/250510.0005139850.002184435
path:00140_4Steroid hormone biosynthesis4/37832/250510.0013137070.003783719
path:00430_1Taurine and hypotaurine metabolism2/3784/250510.001335430.003783719
path:00061_3Fatty acid biosynthesis2/3785/250510.0022034990.004682436
path:00650_1Butanoate metabolism2/3785/250510.0022034990.004682436
path:00565_3Ether lipid metabolism3/37821/250510.0037045220.00699743
path:00260_2Glycine, serine and threonine metabolism2/3787/250510.0045355970.007455604
path:00140_7Steroid hormone biosynthesis3/37823/250510.0048242140.007455604
path:00230_1Purine metabolism4/37858/250510.0114126180.016167876
path:00480_1Glutathione metabolism3/37838/250510.0194664350.025456107
path:00230_2Purine metabolism3/37843/250510.0269572070.032733751
path:00071_1Fatty acid metabolism2/37819/250510.0327858660.037157315
path:00240_1Pyrimidine metabolism2/37848/250510.1634417220.17365683
path:00310_1Lysine degradation1/37821/250510.273427560.27342756
DISCUSSION

Accumulative studies indicate that metabolism is substantial for cancer initiation and progression. Alterations of gene-related with metabolism in tumors provide increased energy for cancer cell proliferation even under nutritional deficient or hypoxia condition[18]. Integrative analysis of metabolic pathway and metabolites facilitates us to better understanding of the underlying mechanism and potential drug-targets of pancreatic cancer[13,19].

Following analysis of the pancreas development database and pancreatic cancer database, we established the 202 genes with dev-Dw and can-Up, which were mainly associated with cell proliferation. Consistent to previous studies, our results confirm that the uncontrolled proliferative activity of cancer cell is the most remarkable hallmarks of carcinogenesis, increased along with tumor progression[22]. Conversely, the proliferative capacity is decreased along with normal pancreas development[20-22]. DNA replication is occurred accompanied by cell proliferation. Due to the infinite hyperplasia of cancer cell, the DNA replication is continuously up-regulated in pancreatic cancer[23,24]. Due to the limited proliferation of healthy organs, DNA replication continuously downregulates pancreatic development[25]. We established 141 genes with dev-Up, can-Dw, including ‘‘T Cell Proliferation” and ‘‘Innate Immune Response In Mucosa''. As a dynamic process, the carcinogenesis, associated with immunoediting[26]. Along with cancer progressed, the host immunosurveillance were suppressed, which in turn led to cancer immune escape. Interestingly, though analysis of metabolism pathways, we found that genes related with steroid hormone biosynthesis were dysregulated during cancer development. Steroid can elicits immunosuppressive effects and restricts T cell-mediated cancer eradication[27]. We therefore hypothesized that enhanced steroid production were the primarily cause for cancer immunoediting.

After analysis of the pancreas development database and pancreatic cancer database, we established that Glycerophospholipid metabolism, as an essential subpathway in lipid metabolism, was the most significant sub-pathway. Many pieces of research have been shown that lipid metabolism was strongly linked with pancreas cancer. The pancreatic lipase, as a lipolytic enzyme, was thought to be one of the predictors for prognosis and cancer-specific mortality in pancreas cancer[13]. Pancreatic lipase 1 and 2 is a significant lipase for lipid hydrolysis in pancreatic cancer patients, which is significantly reduced compared to healthy controls. Additionally, the inversed expressed genes and metabolites in the glycerophospholipid metabolism are associated with pancreatic cancer. SLC44A4, closely associated with acetylcholine synthesis and transport, was markedly upregulated in advanced and undifferentiated epithelial tumors, especially in prostate and pancreatic cancer[28,29]. Choline could decrease the risk of developing pancreatic cancer[29]. Besides, choline and phosphocholine have been reported to be associated with other cancers, including breast cancer, ovarian tumor[30-31]. Some researchers believe that increased choline and phosphocholine are the critical aspects of tumor metabolism and tumor cell migration[32]. Besides lipid metabolism, recent studies also reveal that activation of glycolysis and citrate cycle pathway promotes cancer development[33].

Above all, here we identified a series of gene related with metabolism via bioinformatics analysis, which are crucial for cancer development. We believe that our findings may provide potential targets for the treatment or prognosis of pancreatic cancer.

ARTICLE HIGHLIGHTS
Research background

Pancreatic diseases remain as one of the most feared and clinically challenging diseases to treat despite continual improvements in therapies.

Research motivation

To develop agents into a targeted drug for explicitly killing cancer cells.

Research objectives

To explore the molecular interpretation patterns of pancreas development and cancer progression.

Research methods

This study used the ANOVA method, self-organizing map-singular value decomposition analysis, enrichment analysis, and hypergeometric test.

Research results

The results investigate continuously dysregulated interpretation patterns in pancreas development and pancreatic cancer.

Research conclusions

Integrative analysis of continuously dysregulated interpretation patterns to establish the inverse interpretation in metabolites and gene levels. Through integrating the genes with metabolites, some key abnormal regions of metabolic pathways have been established.

Research perspectives

With the increase of human disease database, a larger-scale integrative analysis is needed for the correlation with pancreas development and cancer. We believe the more convince underlying mechanism and potential drug development targets could be supposed by the larger-scale development and integrative cancer analysis in further. Also, this method could be used for other diseases investigation.

ACKNOWLEDGEMENTS

We appreciate Dr. Xi Huang for supports in editing the manuscript.

Footnotes

Manuscript source: Unsolicited manuscript

Specialty type: Gastroenterology and hepatology

Country of origin: China

Peer-review report classification

Grade A (Excellent): A

Grade B (Very good): B, B

Grade C (Good): C, C

Grade D (Fair): 0

Grade E (Poor): 0

P-Reviewer: Ahmad Z, Bolshakova GB, Löhr JM, Misiakos EP, Zavras N S-Editor: Yan JP L-Editor: Filipodia E-Editor: Qi LL

References
1.  GBD 2015 Mortality and Causes of Death Collaborators. Global, regional, and national life expectancy, all-cause mortality, and cause-specific mortality for 249 causes of death, 1980-2015: A systematic analysis for the Global Burden of Disease Study 2015. Lancet. 2016;388:1459-1544.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 4353]  [Cited by in F6Publishing: 4031]  [Article Influence: 503.9]  [Reference Citation Analysis (1)]
2.  Cascinu S, Falconi M, Valentini V, Jelic S; ESMO Guidelines Working Group. Pancreatic cancer: ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-up. Ann Oncol. 2010;21:v55-v58.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 107]  [Cited by in F6Publishing: 121]  [Article Influence: 8.6]  [Reference Citation Analysis (0)]
3.  Hidalgo M. Pancreatic cancer. N Engl J Med. 2010;362:1605-1617.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 2078]  [Cited by in F6Publishing: 2117]  [Article Influence: 151.2]  [Reference Citation Analysis (2)]
4.  Bond-Smith G, Banga N, Hammond TM, Imber CJ. Pancreatic adenocarcinoma. BMJ. 2012;344:e2476.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 109]  [Cited by in F6Publishing: 119]  [Article Influence: 9.9]  [Reference Citation Analysis (0)]
5.  Jennings RE, Berry AA, Strutt JP, Gerrard DT, Hanley NA. Human pancreas development. Development. 2015;142:3126-3137.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 186]  [Cited by in F6Publishing: 176]  [Article Influence: 22.0]  [Reference Citation Analysis (0)]
6.  Quilichini E, Haumaitre C. Implication of epigenetics in pancreas development and disease. Best Pract Res Clin Endocrinol Metab. 2015;29:883-898.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 24]  [Cited by in F6Publishing: 21]  [Article Influence: 2.3]  [Reference Citation Analysis (0)]
7.  Fleming AK, Storz P. Protein kinase C isoforms in the normal pancreas and in pancreatic disease. Cell Signal. 2017;40:1-9.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 11]  [Cited by in F6Publishing: 12]  [Article Influence: 1.7]  [Reference Citation Analysis (0)]
8.  Vinogradova TV, Sverdlov ED. PDX1: A Unique Pancreatic Master Regulator Constantly Changes Its Functions during Embryonic Development and Progression of Pancreatic Cancer. Biochemistry (Mosc). 2017;82:887-893.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 13]  [Cited by in F6Publishing: 14]  [Article Influence: 2.0]  [Reference Citation Analysis (0)]
9.  Roy N, Takeuchi KK, Ruggeri JM, Bailey P, Chang D, Li J, Leonhardt L, Puri S, Hoffman MT, Gao S, Halbrook CJ, Song Y, Ljungman M, Malik S, Wright CV, Dawson DW, Biankin AV, Hebrok M, Crawford HC. PDX1 dynamically regulates pancreatic ductal adenocarcinoma initiation and maintenance. Genes Dev. 2016;30:2669-2683.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 76]  [Cited by in F6Publishing: 65]  [Article Influence: 8.1]  [Reference Citation Analysis (0)]
10.  Yang S, He P, Wang J, Schetter A, Tang W, Funamizu N, Yanaga K, Uwagawa T, Satoskar AR, Gaedcke J, Bernhardt M, Ghadimi BM, Gaida MM, Bergmann F, Werner J, Ried T, Hanna N, Alexander HR, Hussain SP. A Novel MIF Signaling Pathway Drives the Malignant Character of Pancreatic Cancer by Targeting NR3C2. Cancer Res. 2016;76:3838-3850.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 128]  [Cited by in F6Publishing: 186]  [Article Influence: 23.3]  [Reference Citation Analysis (0)]
11.  Kirby MK, Ramaker RC, Gertz J, Davis NS, Johnston BE, Oliver PG, Sexton KC, Greeno EW, Christein JD, Heslin MJ, Posey JA, Grizzle WE, Vickers SM, Buchsbaum DJ, Cooper SJ, Myers RM. RNA sequencing of pancreatic adenocarcinoma tumors yields novel expression patterns associated with long-term survival and reveals a role for ANGPTL4. Mol Oncol. 2016;10:1169-1182.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 63]  [Cited by in F6Publishing: 56]  [Article Influence: 7.0]  [Reference Citation Analysis (0)]
12.  He W, Wu J, Shi J, Huo YM, Dai W, Geng J, Lu P, Yang MW, Fang Y, Wang W, Zhang ZG, Habtezion A, Sun YW, Xue J. IL22RA1/STAT3 Signaling Promotes Stemness and Tumorigenicity in Pancreatic Cancer. Cancer Res. 2018;78:3293-3305.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 70]  [Cited by in F6Publishing: 77]  [Article Influence: 12.8]  [Reference Citation Analysis (0)]
13.  Zhang G, He P, Tan H, Budhu A, Gaedcke J, Ghadimi BM, Ried T, Yfantis HG, Lee DH, Maitra A, Hanna N, Alexander HR, Hussain SP. Integration of metabolomics and transcriptomics revealed a fatty acid network exerting growth inhibitory effects in human pancreatic cancer. Clin Cancer Res. 2013;19:4983-4993.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 186]  [Cited by in F6Publishing: 227]  [Article Influence: 20.6]  [Reference Citation Analysis (0)]
14.  Budhu A, Terunuma A, Zhang G, Hussain SP, Ambs S, Wang XW. Metabolic profiles are principally different between cancers of the liver, pancreas and breast. Int J Biol Sci. 2014;10:966-972.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 21]  [Cited by in F6Publishing: 23]  [Article Influence: 2.3]  [Reference Citation Analysis (0)]
15.  Battini S, Faitot F, Imperiale A, Cicek AE, Heimburger C, Averous G, Bachellier P, Namer IJ. Metabolomics approaches in pancreatic adenocarcinoma: Tumor metabolism profiling predicts clinical outcome of patients. BMC Med. 2017;15:56.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 44]  [Cited by in F6Publishing: 38]  [Article Influence: 5.4]  [Reference Citation Analysis (0)]
16.  An Y, Cai H, Yang Y, Zhang Y, Liu S, Wu X, Duan Y, Sun D, Chen X. Identification of ENTPD8 and cytidine in pancreatic cancer by metabolomic and transcriptomic conjoint analysis. Cancer Sci. 2018;109:2811-2821.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 7]  [Cited by in F6Publishing: 9]  [Article Influence: 1.5]  [Reference Citation Analysis (0)]
17.  Zhang C, Li C, Xu Y, Feng L, Shang D, Yang X, Han J, Sun Z, Li Y, Li X. Integrative analysis of lung development-cancer expression associations reveals the roles of signatures with inverse expression patterns. Mol Biosyst. 2015;11:1271-1284.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 6]  [Cited by in F6Publishing: 6]  [Article Influence: 0.8]  [Reference Citation Analysis (0)]
18.  Le A, Rajeshkumar NV, Maitra A, Dang CV. Conceptual framework for cutting the pancreatic cancer fuel supply. Clin Cancer Res. 2012;18:4285-4290.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 47]  [Cited by in F6Publishing: 49]  [Article Influence: 4.5]  [Reference Citation Analysis (0)]
19.  Hanahan D, Weinberg RA. Hallmarks of cancer: The next generation. Cell. 2011;144:646-674.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 39812]  [Cited by in F6Publishing: 43204]  [Article Influence: 3323.4]  [Reference Citation Analysis (4)]
20.  Ames BN, Shigenaga MK, Gold LS. DNA lesions, inducible DNA repair, and cell division: Three key factors in mutagenesis and carcinogenesis. Environ Health Perspect. 1993;101:35-44.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 144]  [Cited by in F6Publishing: 145]  [Article Influence: 4.7]  [Reference Citation Analysis (0)]
21.  Fürthauer M, González-Gaitán M. Endocytosis, asymmetric cell division, stem cells and cancer: Unus pro omnibus, omnes pro uno. Mol Oncol. 2009;3:339-353.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 22]  [Cited by in F6Publishing: 23]  [Article Influence: 1.5]  [Reference Citation Analysis (0)]
22.  Knudsen ES, Balaji U, Mannakee B, Vail P, Eslinger C, Moxom C, Mansour J, Witkiewicz AK. Pancreatic cancer cell lines as patient-derived avatars: Genetic characterisation and functional utility. Gut. 2018;67:508-520.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 63]  [Cited by in F6Publishing: 66]  [Article Influence: 11.0]  [Reference Citation Analysis (0)]
23.  Cancer Genome Atlas Research Network. Electronic address: andrew_aguirre@dfci.harvard.edu.; Cancer Genome Atlas Research Network. Integrated Genomic Characterization of Pancreatic Ductal Adenocarcinoma. Cancer Cell. 2017;32:185-203.e13.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 1247]  [Cited by in F6Publishing: 1162]  [Article Influence: 166.0]  [Reference Citation Analysis (0)]
24.  Cicenas J, Kvederaviciute K, Meskinyte I, Meskinyte-Kausiliene E, Skeberdyte A, Cicenas J. KRAS, TP53, CDKN2A, SMAD4, BRCA1, and BRCA2 Mutations in Pancreatic Cancer. Cancers (Basel). 2017;9:pii: E42.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 124]  [Cited by in F6Publishing: 162]  [Article Influence: 23.1]  [Reference Citation Analysis (0)]
25.  Anderson RM, Bosch JA, Goll MG, Hesselson D, Dong PD, Shin D, Chi NC, Shin CH, Schlegel A, Halpern M, Stainier DY. Loss of Dnmt1 catalytic activity reveals multiple roles for DNA methylation during pancreas development and regeneration. Dev Biol. 2009;334:213-223.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 121]  [Cited by in F6Publishing: 117]  [Article Influence: 7.8]  [Reference Citation Analysis (0)]
26.  Aroldi F, Zaniboni A. Immunotherapy for pancreatic cancer: Present and future. Immunotherapy. 2017;9:607-616.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 31]  [Cited by in F6Publishing: 33]  [Article Influence: 5.5]  [Reference Citation Analysis (0)]
27.  Patel S, Choudhary M, Chandra RK, Bhardwaj AK, Tripathi MK. Sex steroids exert a suppressive effect on innate and cell mediated immune responses in fresh water teleost, Channa punctatus. Dev Comp Immunol. 2019;100:103415.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 4]  [Cited by in F6Publishing: 5]  [Article Influence: 1.0]  [Reference Citation Analysis (0)]
28.  Mattie M, Raitano A, Morrison K, Morrison K, An Z, Capo L, Verlinsky A, Leavitt M, Ou J, Nadell R, Aviña H, Guevara C, Malik F, Moser R, Duniho S, Coleman J, Li Y, Pereira DS, Doñate F, Joseph IB, Challita-Eid P, Benjamin D, Stover DR. The Discovery and Preclinical Development of ASG-5ME, an Antibody-Drug Conjugate Targeting SLC44A4-Positive Epithelial Tumors Including Pancreatic and Prostate Cancer. Mol Cancer Ther. 2016;15:2679-2687.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 25]  [Cited by in F6Publishing: 29]  [Article Influence: 3.6]  [Reference Citation Analysis (0)]
29.  Song P, Rekow SS, Singleton CA, Sekhon HS, Dissen GA, Zhou M, Campling B, Lindstrom J, Spindel ER. Choline transporter-like protein 4 (CTL4) links to non-neuronal acetylcholine synthesis. J Neurochem. 2013;126:451-461.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 28]  [Cited by in F6Publishing: 30]  [Article Influence: 2.7]  [Reference Citation Analysis (0)]
30.  Aboagye EO, Bhujwalla ZM. Malignant transformation alters membrane choline phospholipid metabolism of human mammary epithelial cells. Cancer Res. 1999;59:80-84.  [PubMed]  [DOI]  [Cited in This Article: ]
31.  Ackerstaff E, Glunde K, Bhujwalla ZM. Choline phospholipid metabolism: A target in cancer cells? J Cell Biochem. 2003;90:525-533.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 256]  [Cited by in F6Publishing: 264]  [Article Influence: 12.6]  [Reference Citation Analysis (0)]
32.  Zhan B, Wen S, Lu J, Shen G, Lin X, Feng J, Huang H. Identification and causes of metabonomic difference between orthotopic and subcutaneous xenograft of pancreatic cancer. Oncotarget. 2017;8:61264-61281.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 25]  [Cited by in F6Publishing: 25]  [Article Influence: 3.6]  [Reference Citation Analysis (0)]
33.  Ren JG, Seth P, Ye H, Guo K, Hanai JI, Husain Z, Sukhatme VP. Citrate Suppresses Tumor Growth in Multiple Models through Inhibition of Glycolysis, the Tricarboxylic Acid Cycle and the IGF-1R Pathway. Sci Rep. 2017;7:4537.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 62]  [Cited by in F6Publishing: 79]  [Article Influence: 11.3]  [Reference Citation Analysis (0)]