Circadian Clock Genes Are Correlated with Prognosis and Immune Cell Infiltration in Colon Adenocarcinoma

Abstract

Background: Colon adenocarcinoma (COAD) is a malignancy with a high incidence and is associated with poor quality of life. Dysfunction of circadian clock genes and disruption of normal rhythms are associated with the occurrence and progression of many cancer types. However, studies that systematically describe the prognostic value and immune-related functions of circadian clock genes in COAD are lacking.

Methods: Genomic data obtained from The Cancer Genome Atlas (TCGA) database was analyzed for expression level, mutation status, potential biological functions, and prognostic performance of core circadian clock genes in COAD. Their correlations with immune infiltration and TMB/MSI score were analyzed by Spearman's correlation analysis. Pearson's correlation analysis was performed to analyze their associations with drug sensitivity. Lasso Cox regression analysis was performed to construct a prognosis signature. Moreover, an mRNA-miRNA-lncRNA regulatory axis was also detected by ceRNA network.

Results: In COAD tissues, the mRNA levels of CLOCK, CRY1, and NR1D1 were increased, while the mRNA levels of ARNTL, CRY2, PER1, PER3, and RORA were decreased. We also summarized the relative genetic mutation variation landscape. GO and KEGG pathway analyses demonstrated that these circadian clock genes were primarily correlated with the regulation of circadian rhythms and glucocorticoid receptor signaling pathways. COAD patients with high CRY2, NR1D1, and PER2 expression had worse prognosis. A prognostic model constructed based on the 9 core circadian clock genes predicted the COAD patients' overall survival with medium to high accuracy. A significant association between prognostic circadian clock genes and immune cell infiltration was found. Moreover, the lncRNA KCNQ1OT1/hsa-miRNA-32-5p/PER2/CRY2 regulatory axis in COAD was also detected through a mRNA-miRNA-lncRNA network.

Conclusion: Our results identified CRY2, NR1D1, and PER2 as potential prognostic biomarkers for COAD patients and correlated their expression with immune cell infiltration. The lncRNA KCNQ1OT1/hsa-miRNA-32-5p/PER2/CRY2 regulatory axis was detected in COAD and might play a vital role in the occurrence and progression of COAD.