An endogenous molecular clockwork drives various cellular pathways including metabolism and the cell cycle. Its dysregulation is able to prompt pathological phenotypes including cancer. Besides dramatic metabolic alterations, cancer cells display severe changes in the clock phenotype with likely consequences in tumor progression and treatment response. In this study, we use a comprehensive systems-driven approach to investigate the effect of clock disruption on metabolic pathways and its impact on drug response in a cellular model of colon cancer progression. We identified distinctive time-related transcriptomic and metabolic features of a primary tumor and its metastatic counterpart. A mapping of the expression data to a comprehensive genome-scale reconstruction of human metabolism allowed for the in-depth functional characterization of 24 h-oscillating transcripts and pointed to a clock-driven metabolic reprogramming in tumorigenesis. In particular, we identified a set of five clock-regulated glycolysis genes, ALDH3A2, ALDOC, HKDC1, PCK2, and PDHB with differential temporal expression patterns. These findings were validated in organoids and in primary fibroblasts isolated from normal colon and colon adenocarcinoma from the same patient. We further identified a reciprocal connection of HKDC1 to the clock in the primary tumor, which is lost in the metastatic cells. Interestingly, a disruption of the core-clock gene BMAL1 impacts on HKDC1 and leads to a time-dependent rewiring of metabolism, namely an increase in glycolytic activity, as well as changes in treatment response. This work provides novel evidence regarding the complex interplay between the circadian clock and metabolic alterations in carcinogenesis and identifies new connections between both systems with pivotal roles in cancer progression and response to therapy.

The mammalian circadian clock and the cell cycle are two major biological oscillators whose coupling influences cell fate decisions. In the present study, we use a model-driven experimental approach to investigate the interplay between clock and cell cycle components and the dysregulatory effects of RAS on this coupled system. In particular, we focus on the Ink4a/Arf locus as one of the bridging clock-cell cycle elements. Upon perturbations by the rat sarcoma viral oncogene (RAS), differential effects on the circadian phenotype were observed in wild-type and Ink4a/Arf knock-out mouse embryonic fibroblasts (MEFs), which could be reproduced by our modelling simulations and correlated with opposing cell cycle fate decisions. Interestingly, the observed changes can be attributed to in silico phase shifts in the expression of core-clock elements. A genome-wide analysis revealed a set of differentially expressed genes that form an intricate network with the circadian system with enriched pathways involved in opposing cell cycle phenotypes. In addition, a machine learning approach complemented by cell cycle analysis classified the observed cell cycle fate decisions as dependent on Ink4a/Arf and the oncogene RAS and highlighted a putative fine-tuning role of Bmal1 as an elicitor of such processes, ultimately resulting in increased cell proliferation in the Ink4a/Arf knock-out scenario. This indicates that the dysregulation of the core-clock might work as an enhancer of RAS-mediated regulation of the cell cycle. Our combined in silico and in vitro approach highlights the important role of the circadian clock as an Ink4a/Arf-dependent modulator of oncogene-induced cell fate decisions, reinforcing its function as a tumour-suppressor and the close interplay between the clock and the cell cycle network. In mammals, the circadian clock controls the punctual regulation of biological processes, which, in turn, affect physiology and behaviour, allowing for the synchronisation of internal time to environmental light-dark cycles. Malfunctions of the circadian clock are associated with pathological phenotypes including cancer. Given the range of molecular time-dependent processes, including metabolism, DNA repair, and the cell cycle, the clock is hypothesised to act as a tumour suppressor. With the help of mathematical modelling and whole-genome analysis combined with machine learning, we investigated the RAS-dependent dysregulation of the circadian clock. We find that the tumour-suppressor Ink4a/Arf acts as a key mediator of RAS oncogene-induced changes in the circadian system, thereby mediating the interplay between the clock and the cell cycle.

Cancer cells interrelate with the bordering host microenvironment that encompasses the extracellular matrix and a nontumour cellular component comprising fibroblasts and immune-competent cells. The tumour microenvironment modulates cancer onset and progression, but the molecular factors managing this interaction are not fully understood. Malignant transformation of a benign tumour is among the first crucial events in colorectal carcinogenesis. The role of tumour stroma fibroblasts is well-described in cancer, but less well-characterized in benign tumours. In the current work we utilized fibroblasts isolated from tubulovillous adenoma, which has high risk for malignant transformation, to study the interaction between benign tumour stroma and the circadian clock machinery. We explored the role of the biological clock in this interplay taking advantage of an experimental model, represented by the co-culture of colon cancer cells with normal fibroblasts or tumour-associated fibroblasts, isolated from human colorectal tumour specimens. When co-cultured with tumour-associated fibroblasts, colon cancer cells showed alterations in their circadian and metabolic parameters, with decreased apoptosis, increased colon cancer cell viability, and increased resistance to chemotherapeutic agents. In conclusion, the interactions among colon cancer cells and tumour-associated fibroblasts affect the molecular clockwork and seem to aggravate malignant cell phenotypes, suggesting a detrimental effect of this interplay on cancer dynamics.

Accumulating evidence points to a significant role of the circadian clock in the regulation of splicing in various organisms, including mammals. Both dysregulated circadian rhythms and aberrant pre-mRNA splicing are frequently implicated in human disease, in particular in cancer. To investigate the role of the circadian clock in the regulation of splicing in a cancer progression context at the systems-level, we conducted a genome-wide analysis and compared the rhythmic transcriptional profiles of colon carcinoma cell lines SW480 and SW620, derived from primary and metastatic sites of the same patient, respectively. We identified spliceosome components and splicing factors with cell-specific circadian expression patterns including SRSF1, HNRNPLL, ESRP1, and RBM 8A, as well as altered alternative splicing events and circadian alternative splicing patterns of output genes (e.g., VEGFA, NCAM1, FGFR2, CD44) in our cellular model. Our data reveals a remarkable interplay between the circadian clock and pre-mRNA splicing with putative consequences in tumor progression and metastasis.

Mounting evidence points to a role of the circadian clock in the temporal regulation of post-transcriptional processes in mammals, including alternative splicing (AS). In this study, we carried out a computational analysis of circadian and ultradian rhythms on the transcriptome level to characterise the landscape of rhythmic AS events in published datasets covering 76 tissues from mouse and olive baboon. Splicing-related genes with 24-h rhythmic expression patterns showed a bimodal distribution of peak phases across tissues and species, indicating that they might be controlled by the circadian clock. On the output level, we identified putative oscillating AS events in murine microarray data and pairs of differentially rhythmic splice isoforms of the same gene in baboon RNA-seq data that peaked at opposing times of the day and included oncogenes and tumour suppressors. We further explored these findings using a new circadian RNA-seq dataset of human colorectal cancer cell lines. Rhythmic isoform expression patterns differed between the primary tumour and the metastatic cell line and were associated with cancer-related biological processes, indicating a functional role of rhythmic AS that might be implicated in tumour progression. Our data shows that rhythmic AS events are widespread across mammalian tissues and might contribute to a temporal diversification of the proteome.