A circadian clock gene-related signature for predicting prognosis and its association with sorafenib response in hepatocellular carcinoma

Background: Hepatocellular carcinoma (HCC), one of the highest causes of cancer-associated death, has effective treatments, especially for patients with advanced HCC. Circadian rhythm participates in several important physiological functions, and its chronic disruption results in many disordered diseases, including cancer. However, the role of circadian rhythm in the overall survival (OS) of patients with HCC remains unclear. Methods: We investigated the expression, copy number variation (CNV), and mutation profiles of core circadian clock genes in normal and tumor tissues. We developed and validated a messenger RNA signature (mRNASig) based on prognostic circadian clock genes. A set of bioinformatic tools were applied for functional annotation and tumor-associated microenvironment (TME) analysis. Results: Core circadian clock genes were disrupted in terms of the transcription and CNV of HCC samples. The mRNASig, including NPAS2, NR1D1, PER1, RORC, and TIMELESS, was constructed. We divided patients with HCC into high-risk group and low-risk group based on the median value of the risk score. The high-risk group had a poorer prognosis than the low-risk group. The high-risk group was associated with malignant processes (e.g., proliferation, oncogenic pathway, DNA repair), metabolism, and tumor mutational burden (TMB). Surprisingly, the low-risk group was associated with enriched angiogenesis and was linked to enhanced response to sorafenib. Moreover, the high-risk group showed poor infiltration of CD8 T cells and natural killer cells accompanied by higher expression of CTLA4, PDCD1, TIGIT, and TIM3. Additionally, the mRNASig was associated with TMB. Conclusions: The mRNASig based on core circadian clock genes is a potential prognostic signature and therapeutic strategy and is significantly associated with the malignant biology of HCC.