Circadian Genes MBOAT2/CDA/LPCAT2/B4GALT5 in the Metabolic Pathway Serve as New Biomarkers of PACA Prognosis and Immune Infiltration

SIMPLE SUMMARY: In this study, we found 299 circadian genes that were differentially expressed in pancreatic cancer (PACA). The expression levels of MBOAT2, CDA, LPCAT2 and B4GALT5, which were enriched in metabolic pathways, were not only highly connected with overall survival and other clinical parameters in PACA patients, but were also closely correlated with the level of immune cell infiltration. In conclusion, there was a close link between circadian genes, cancer marker pathways and immune infiltration. These findings could lay the foundation for the use of personalized chronotherapy and immune checkpoint inhibitors in the treatment of pancreatic cancer. ABSTRACT: Pancreatic cancer (PACA) is a highly malignant tumor with a poor prognosis. Recent studies have discovered substantial differences in the expression levels of several circadian genes in PACA samples compared to normal samples. The goal of this research was to find differentially expressed rhythm genes (DERGs) in PACA samples and determine their role in the development of PACA. A total of 299 DERGs were identified in PACA, including 134 downregulated genes and 165 upregulated genes. DERGs were significantly abundant in the metabolic pathway and immune response pathways, according to GO and KEGG analyses. Survival analyses showed that PACA patients who had higher expression levels of MBOAT2/CDA/LPCAT2/B4GALT5 had shorter overall survival times. Using cell assay verification, the mRNA levels of MBOAT2/CDA/LPCAT2/B4GALT5 in Patu-8988 and PNAC-1 cells were found to be significantly higher than those in HPDE6-C7 cells, which was in line with previous studies on PACA patient data. Through conducting univariate Cox analysis, it was determined that MBOAT2/CDA/LPCAT2/B4GALT5 expression, age and grade were all high-risk factors. The MBOAT2/CDA/LPCAT2/B4GALT5 genes were independently correlated with overall survival, according to the multivariate Cox analysis. The proportion of immune cells in PACA and normal samples significantly changed, according to the immune infiltration analysis. Furthermore, MBOAT2/CDA/LPCAT2/B4GALT5 expression levels were significantly related to the level of immune cell infiltration. The protein–protein interaction network of the MBOAT2/CDA/LPCAT2/B4GALT5 genes included 54 biological nodes and 368 interacting genes. In conclusion, the finding of these DERGs adds to the investigation of the molecular processes underlying the onset and progression of PACA. In the future, DERGs may serve as prognostic and diagnostic biomarkers as well as drug targets for chronotherapy in PACA patients.