Circadian Regulation Patterns With Distinct Immune Landscapes in Gliomas Aid in the Development of a Risk Model to Predict Prognosis and Therapeutic Response

Abstract

Circadian disruption in tumorigenesis has been extensively studied, but how circadian rhythm (CR) affects the formation of tumor microenvironment (TME) and the crosstalk between TME and cancer cells is largely unknown, especially in gliomas. Herein, we retrospectively analyzed transcriptome data and clinical parameters of glioma patients from public databases to explore circadian rhythm-controlled tumor heterogeneity and characteristics of TME in gliomas. Firstly, we pioneered the construction of a CR gene set collated from five datasets and review literatures. Unsupervised clustering was used to identify two CR clusters with different CR patterns on the basis of the expression of CR genes. Remarkably, the CR cluster-B was characterized by enriched myeloid cells and activated immune-related pathways. Next, we applied principal component analysis to construct a CRscore to quantify CR patterns of individual tumors, and the function of the CRscore in prognostic prediction was further verified by univariate and multivariate regression analyses in combination with a nomogram. The CRscore could not only be an independent factor to predict prognosis of glioma patients but also guide patients to choose suitable treatment strategies: immunotherapy or chemotherapy. A glioma patient with a high CRscore might respond to immune checkpoint blockade, whereas one with a low CRscore could benefit from chemotherapy. In this study, we revealed that circadian rhythms modulated tumor heterogeneity, TME diversity, and complexity in gliomas. Evaluating the CRscore of an individual tumor would contribute to gaining a greater understanding of the tumor immune status of each patient, enhancing the accuracy of prognostic prediction, and suggesting more effective treatment options.