Identification of TIMELESS and RORA as key clock molecules of non-small cell lung cancer and the comprehensive analysis.

Background: The incidence rate of non-small cell lung cancer (NSCLC) has been increasing worldwide, and the correlation of circadian rhythm disruption with a raised risk of cancer and worse prognosis has been shown by accumulating evidences recently. On the other hand, drug resistance and the impact of tumor heterogeneity have been inevitable in NSCLC therapy. These both lead to an urgent need to identify more useful prognostic and predictive markers for NSCLC diagnosis and treatment, especially on the aspect of circadian clock genes.

Methods: The expression of the main clock genes in cancer was probed with TIMER and Oncomine databases. The prognostic value of key clock genes was probed systematically with the Kaplan-Meier estimate and Cox regression on samples from TCGA database. RT-qPCR was performed on patient tissue samples to further validate the results from databases. The functional enrichment analysis was performed using the "ClusterProfiler" R package, and the correlation of key clock genes with tumor mutation burden, immune checkpoint, and immune infiltration levels were also assessed using multiple algorithms including TIDE, TIMER2.0, and XCELL.

Results: TIMELESS was significantly upregulated in lung tissue of clinical lung cancer patients as well as TCGA and Oncomine databases, while RORA was downregulated. Multivariate Cox regression analysis indicated that TIMELESS (P = 0.004, HR = 1.21 [1.06, 1.38]) and RORA (P = 0.047, HR = 0.868 [0.755, 0.998]) has a significant correlation with overall survival in NSCLC. Genes related to TIMELESS were enriched in the cell cycle and immune system, and the function of RORA was mainly focused on oncogenic signaling pathways or glycosylation and protein activation. Also, TIMELESS was positively correlated with tumor mutation burden while RORA was negatively correlated with it. TIMELESS and RORA were also significantly correlated with immune checkpoint and immune infiltration levels in NSCLC. Additionally, TIMELESS showed a significant positive relationship with lipid metabolism.

Conclusions: TIMELESS and RORA were identified as key clock genes in NSCLC, and were independent prognostic factors for overall survival in NSCLC. The function of them were assessed in many aspects, indicating the strong potential of the two genes to serve as biomarkers for NSCLC progression and prognosis.