Exploring the key genetic association between chronic pancreatitis and pancreatic ductal adenocarcinoma through integrated bioinformatics

Abstract

Background: Pancreatic ductal adenocarcinoma (PDAC) develops rapidly and has a poor prognosis. It has been demonstrated that pancreatic ductal adenocarcinoma and chronic pancreatitis (CP) have a close connection. However, the underlying mechanisms for chronic pancreatitis transforming into pancreatic ductal adenocarcinoma are still unclear. The purpose of this study was to identify real hub genes in the development of chronic pancreatitis and pancreatic ductal adenocarcinoma. Methods: RNA-seq data of chronic pancreatitis and pancreatic ductal adenocarcinoma were downloaded from the Gene Expression Omnibus (GEO) database. Weighted gene co-expression network analysis (WGCNA) was performed to construct a gene co-expression network between chronic pancreatitis and pancreatic ductal adenocarcinoma. GEO2R and a Venn diagram were used to identify differentially expressed genes. Then visualized networks were constructed with ClueGO, and modules of PPI network were calculated by MCODE plugin. Further validation of the results was carried out in two additional cohorts. Analyses of CEL-coexpressed genes and regulators including miRNAs and transcription factors were performed by using the corresponding online web tool. Finally, the influence of CEL in the tumor immune microenvironment (TIME) was assessed by immune contextual analysis. Results: With the help of WGCNA and GEO2R, four co-expression modules and six hub genes were identified, respectively. ClueGO enrichment analysis and MCODE cluster analysis revealed that the dysfunctional transport of nutrients and trace elements might contribute to chronic pancreatitis and pancreatic ductal adenocarcinoma development. The real hub gene CEL was identified with a markedly low expression in pancreatic ductal adenocarcinoma in external validation sets. According to the miRNA-gene network construction, hsa-miR-198 may be the key miRNA. A strong correlation exists between CEL and TIME after an evaluation of the influence of CEL in TIME. Conclusion: Our study revealed the dysfunctional transport of nutrients and trace elements may be common pathogenesis of pancreatic ductal adenocarcinoma and chronic pancreatitis. Examination on these common pathways and real hub genes may shed light on the underlying mechanism.

Keywords: CEL; WCGNA; bioinformatics; chronic pancreatitis; differential gene analysis; pancreatic ductal adenocarcinoma.