Abstract
We have previously demonstrated that the pancreas can recover from chronic pancreatitis (CP) lesions in the cerulein-induced mouse model. To explore how pancreatic recovery is achieved at the molecular level, we used RNA-sequencing (seq) and profiled transcriptomes during CP transition to recovery. CP was induced by intraperitoneally injecting cerulein in C57BL/6 mice. Time-matched controls (CON) were given normal saline. Pancreata were harvested from mice 4 days after the final injections (designated as CP and CON) or 4 weeks after the final injections (designated as CP recovery (CPR) and control recovery (CONR)). Pancreatic RNAs were extracted for RNA-seq and quantitative (q) PCR validation. Using RNA-seq, we identified a total of 3,600 differentially expressed genes (DEGs) in CP versus CON and 166 DEGs in CPR versus CONR. There are 132 DEGs overlapped between CP and CPR and 34 DEGs unique to CPR. A number of selected pancreatic fibrosis-relevant DEGs were validated by qPCR. The top 20 gene sets enriched from DEGs shared between CP and CPR are relevant to extracellular matrix and cancer biology, whereas the top 10 gene sets enriched from DEGs specific to CPR are pertinent to DNA methylation and specific signaling pathways. In conclusion, we identified a distinct set of DEGs in association with extracellular matrix and cancer cell activities to contrast CP and CPR. Once during ongoing CP recovery, DEGs relevant to DNA methylation and specific signaling pathways were induced to express. The DEGs shared between CP and CPR and the DEGs specific to CPR may serve as the unique transcriptomic signatures and biomarkers for determining CP recovery and monitoring potential therapeutic responses at the molecular level to reflect pancreatic histological resolution.