Circadian Control of Metabolism by the Clock Component TOC1

Photosynthesis in chloroplasts during the day and mitochondrial respiration during the night execute nearly opposing reactions that are coordinated with the internal cellular status and the external conditions. Here, we describe a mechanism by which the Arabidopsis clock component TIMING OF CAB EXPRESSION1 (TOC1) contributes to the diurnal regulation of metabolism. Proper expression of TOC1 is important for sustaining cellular energy and for the diel and circadian oscillations of sugars, amino acids and tricarboxylic acid (TCA) cycle intermediates. TOC1 binds to the promoter of the TCA-related gene FUMARASE 2 to repress its expression at night, which results in decreased fumarate accumulation in TOC1 over-expressing plants and increased in toc1-2 mutant. Genetic interaction studies confirmed that over-expression of FUMARASE 2 in TOC1 over-expressing plants alleviates the molecular and physiological energy-deprivation phenotypes of TOC1 over-expressing plants. Thus, we propose that the tandem TOC1-FUMARASE 2 is one of the mechanisms that contribute to the regulation of plant metabolism during the day and night.

The circadian clock influences the timing of metabolism in animals, and thus, clock disruption results in metabolic disorders (Xu et al., 2008;Dibner and Schibler, 2015;Gill et al., 2015;Brown, 2016). In plants, metabolism relies in essence on two energy-related organelles: chloroplasts and mitochondria (Sweetlove and Fernie, 2013). Chloroplasts convert light energy into sugars, whereas mitochondria turn the nutrients into energy to cover for the cellular energetic demands (Millar et al., 2011). Photorespiration also supplies a large amount of NADH to mitochondria (Lim et al., 2020). Photosynthesis in chloroplasts and mitochondrial respiration are closely connected  but execute nearly opposing reactions. Thus, these activities must be precisely coordinated in sync with the internal cellular status and the external conditions. Fumarate is one of the intermediates of the TCA cycle (Nunes-Nesi et al., 2013) that accumulates to high amounts in Arabidopsis leaves due to the activity of two fumarase (fumarate hydratase) (FUM) enzymes (FUM1 and FUM2) (Chia et al., 2000;Araújo et al., 2011). Arabidopsis FUM2 activity accounts for the majority of the fumarate accumulation in leaves (Chia et al., 2000;Nunes-Nesi et al., 2007). FUM2 expression and fumarate accumulation follow a diurnal oscillation (Chia et al., 2000;Nunes-Nesi et al., 2007;Pracharoenwattana et al., 2010), with higher accumulation during the light period. Fumarate and starch have been proposed as alternative carbon sinks for photosynthate (Chia et al., 2000;Tschoep et al., 2009;Pracharoenwattana et al., 2010).
The connection of the circadian clock with chloroplasts and photosynthetic activity has been firmly established (McClung et al., 2000;Dodd et al., 2005;Haydon et al., 2013;Noordally et al., 2013;Atkins and Dodd, 2014). Diel or circadian rhythms in metabolism-related gene expression or protein accumulation have been also reported (McClung et al., 2000;Michalecka et al., 2003;Blasing et al., 2005;Elhafez et al., 2006;Gibon et al., 2006;Lee et al., 2010;Flis et al., 2019). The use of single, double and triple clock mutant plants has shown altered gene expression and disrupted accumulation of organic acids and other metabolites (Fukushima et al., 2009;Sanchez-Villarreal et al., 2013;Flis et al., 2019) indicating that the circadian clock is important for proper metabolism. Here, we combine a wide-range of approaches including transcriptional, metabolomic, chromatin immunoprecipitation (ChIP) and genetic interaction studies to show that TOC1 function is important in the control of cellular metabolism at least in part through the oscillatory regulation of fumarase expression. We propose that TOC1 might function as a rheostat in the regulation of metabolism in sync with the day and night cycles.

Plant Material, Growing Conditions and Luminescence Assays
Arabidopsis thaliana lines (Columbia, Col-0) were surface sterilized dried on sterile filter paper in a laminar flow cabinet and sown on plates containing Murashige and Skoog (MS) agar medium supplemented or not (as specified for each experiment) with 3% (w/v) sucrose. Following 48 h of stratification at 4 • C in the dark, plates were transferred to environmentallycontrolled chambers (Inkoa Sistemas) (Huang et al., 2012), TOC1 RNAi (Más et al., 2003a), toc1-2 (NASC, N2107710) and TMG-YFP/toc1-2 (Huang et al., 2012) plants were described elsewhere. Plant transformation was performed by floral dipping protocols using Agrobacterium tumefaciens (GV2260)-mediated transfer (Clough and Bent, 1998). In vivo luminescence assays were performed as previously described (Takahashi et al., 2015). Briefly, 6 day-old seedlings synchronized under LD cycles were transferred to 96-well plates and released into LL conditions. Analyses were performed with a LB960 luminometer (Berthold Technologies) using the Microwin software (Mikrotek Laborsysteme).

Plasmid Construction
The FUM2-ox construct was generated by PCR-mediated amplification of the FUM2 coding sequence (CDS) followed by subcloning into the pENTR/D-TOPO vector (Invitrogen). The resulting vector containing the FUM2 CDS was used to transform chemically competent Escherichia coli cells (one shot TOP10 cells, Gateway R ). The FUM2 CDS was introduced in the plant destination vector pGWB514 (35S pro, C-3xHA) (Nakagawa et al., 2007a,b) by homologous recombination using the LR reaction (Gateway R ). Expression vectors (containing the FUM2 CDS fused to 4x-MYC tag and expressed under the control of the 35S promoter) were also used to transform chemically competent E. coli cells (one shot TOP10, Gateway R ). WT and TOC1-ox plants were transformed with the FUM2-ox construct to generate the single FUM2-ox and the double FUM2-ox/TOC1ox plants, respectively.

RNA Extraction and Reverse Transcription Quantitative PCR
About 20 seedlings of 12-14 days after stratification (das) were collected every 4 h over a diurnal or circadian cycle. RNA was purified using the Maxwell 16 LEV simply RNA Tissue kit (Promega). RNA was incubated with RNase-free TURBO DNase (Ambion) to eliminate genomic DNA contamination. Single strand cDNA was synthesized using 1 µg of RNA using iScript TM Reverse Transcription Supermix for RT-Q-PCR (Biorad) or AffinityScript Q-PCR cDNA Synthesis Kit (Agilent). For qPCR analysis, cDNAs were diluted fivefold with nuclease-free water and qPCR was performed with 10% of diluted cDNA with Brilliant III Ultra-Fast SYBR Green qPCR Master Mix (Agilent) or iTag Universal SYBR Green Supermix (Bio-Rad) in a 96well CFX96 Touch Real-Time PCR detection system (Bio-Rad). The ISOPENTENYL PYROPHOSPHATE:DIMETHYLALLYL PYROPHOSPHATE ISOMERASE 2 (IPP2) gene was used as control (Huang et al., 2012). Data were analyzed using the second derivative maximum method. Resulting Cp values were converted into relative expression values using the comparative Ct method (Livak and Schmittgen, 2001).

Adenylate Measurements by HPLC
The ADP and ATP content was measured by High Performance Liquid Chromatography (HPLC) with 20 mg of fresh weight plant materials (Zhang et al., 2020). Briefly, 0.2 ml of 0.1 M HCL was mixed with the plant materials on ice. 15 µl extract/standard (different concentration) was mixed with 77 µl CP Buffer [62 mM, citric acid monohydrate and 76 mM (Na) 2 HPO 4 × 2H 2 O] and 8 µl 45% chloroacetaldehyde and incubated at 80 • C for 10 min. The mixture was then centrifuged at 16,000×g at 20 • C for 30 min. 90 µl of the supernatant was then measured by the HPLC [Hyperclone C18 (ODS) column (Phenomenex)]. The result was calculated using the standard ADP and ATP gradient.

GC-MS Based Metabolite Profiling
Metabolite analyses were performed essentially as previously described (Lisec et al., 2006). Briefly, seedlings were collected at the indicated time points, immediately frozen in liquid nitrogen, and stored at −80 • C until further analysis. Metabolite extraction was performed by rapid grinding in liquid nitrogen and immediate addition of the extraction buffer (1,400 µL of methanol plus 60 µL of 0.2 mg ribitol mL −1 water). The extraction, derivatization, standard addition, and sample injection were conducted as described (Lisec et al., 2006) with minor modifications in the equipment. An autosampler Gerstel Multi-Purpose system (Gerstel GmbH & Co., KG, Mülheim an der Ruhr, Germany) was used to inject the samples to a chromatograph coupled to a time-of-flight mass spectrometer (GC-MS) system, Leco Pegasus HT TOF-MS (LECO Corporation, St. Joseph, MI, United States). Metabolites were identified by comparison with database entries of authentic standards (Kopka et al., 2005;Schauer et al., 2005). Chromatograms and mass spectra were evaluated using the Chroma TOF 1.0 (Leco 1 ) and the TagFinder 4.0 software (Luedemann et al., 2012). The relative content of metabolites was calculated by normalization of signal intensity to that of ribitol, which was added as an internal standard, and by the fresh weight of the material. All data were also processed using the Xcalibur 4.0 software (Thermo Fisher Scientific, Waltham, MA, United States) to verify the metabolite identification and annotation. Identification and annotation of detected peaks followed the recommendations for reporting metabolite data .

Rosette Phenotyping
Whole rosettes of 35-day old seedlings growing on soil pots were detached and immediately weigh in a precision balance (OHAUS Pioneer PA114C Analytical Balance) to calculate fresh weight. Detached rosettes were immediately placed on top of a transilluminator and photographed (NIKON D7000) for further analysis. Rosette images were processed with ImageJ software (Schneider et al., 2012). First, images were split using the RGB color tool. Blue channel images were then converted to 8 bit black/white images. After the same scale was set for all samples, rosettes were defined using the wand tracing tool and the area and the perimeter was obtained.

Circadian Regulation of Nuclear-Encoded Mitochondrial-Related Transcriptome
To determine the oscillatory transcription of nuclear-encoded mitochondrial-related genes, we used the web tool DIURNAL (Mockler et al., 2007;Michael et al., 2008) and found that about 65% of the mitochondrial-related genes (Chrobok et al., 2016) displayed an oscillatory pattern of expression under diel conditions with spread peak phases but particularly enriched in the middle of the night at ZT17-22 ( Figure 1A). Time course assays by Reverse Transcription-Quantitative Polymerase Chain Reaction (RT-QPCR) confirmed the diurnal oscillation of selected genes belonging to mitochondrial-related metabolic pathways (Figures 1C,D and Supplementary Figures 1A-J).
To ascertain whether the oscillatory pattern of expression is controlled by the circadian clock, we also used the web tool DIURNAL (Mockler et al., 2007;Michael et al., 2008) and performed the analyses under LL conditions. Our studies revealed that around 42% of the mitochondrial-related genes showed a circadian oscillation ( Figure 1B). These results were also confirmed by RT-QPCR analyses ( Supplementary  Figures 1K-T). Furthermore, analyses of ACONITASE1 (ACO1) and ALTERNATIVE NAD(P)H DEHYDROGENASE 1 (NDA1) promoters fused to the LUCIFERASE also uncovered the rhythmic promoter activities under both diel and circadian conditions (Figures 1E,F). Altogether, the data from the web tool DIURNAL, our RT-QPCR and promoter activity analyses, all indicate the diel and circadian oscillation of some nuclearencoded mitochondrial-related genes.

TOC1 Regulates the Circadian Expression of Mitochondrial-Related Genes
Previous transcriptomic profiles showed that the clock components PRR9, 7, and 5 might be important in the control of metabolism (Fukushima et al., 2009). The study was conducted with plants growing in medium containing sucrose and excluded the founder member of the PRR family, TOC1, which plays distinct roles within the circadian system and in the control of particular outputs. Our analyses using TOC1 miss-expressing plants showed that the rhythmic expression of genes encoding proteins belonging to the mitochondrial electron transport chain (ETC), the tricarboxylic acid (TCA) cycle, as well as expression of NADH dehydrogenases, alternative oxidases (AOX) and mitochondrial transporters (TRANSP) was affected in TOC1 over-expressing and mutant plants grown in a medium without sucrose ( Supplementary Figures 2A-C). Peak accumulation appeared to be delayed in TOC1-ox and advanced in toc1-2. Comparative analyses of individual genes in the three genotypes confirmed the altered expression, showing that some genes were down-regulated in TOC1-ox and conversely, up-regulated in toc1-2 (Figures 2A-I). Less frequently, the amplitude of some genes was increased in TOC1-ox ( Figure 2I). As TOC1 functions as a general repressor (Gendron et al., 2012;Huang et al., 2012;Pokhilko et al., 2012), these genes might be indirectly regulated by the miss-expression of TOC1.

TOC1 Regulates the Diurnal and Circadian Accumulation of Central Metabolites and Cellular Energy
To examine metabolite accumulation in Wild-Type (WT) and TOC1 miss-expressing plants, we performed targeted metabolomics by Gas Chromatography/Mass Spectrometry (GC/MS) analyses in plants grown in a medium without sucrose (Supplementary Table 1). In WT plants, a number of metabolites displayed diurnal fluctuations with higher accumulation during the day at ZT7 than at night at ZT19 (Supplementary Figure 3A). Examples include sugars such as glucose, amino acids such as serine, or TCA intermediates such as succinic acid. The oscillatory pattern was altered in TOC1-ox for a subset of these metabolites, with increased accumulation either in the morning or in the evening (Supplementary Figure 3A and Figure 3B). Metabolites with high amplitude fluctuations in WT showed decreased accumulation in TOC1-ox ( Figure 3B and Supplementary Figure 3A) and were slightly increased in toc1-2 (Supplementary Figures 3A,B). Under LL conditions, metabolites also fluctuated with increased accumulation during the subjective day (e.g., glycine) or night (e.g., phenylalanine) ( Figure 3A). The rhythmic fluctuations were altered in TOC1ox and toc1-2 mutant plants (Figures 3A-C and Supplementary  Figure 3B). Amino acids accumulated either during the day or night in TOC1-ox, whereas sugars and fumaric acid displayed a clear reduction. Conversely, fructose, fumaric acid or the amino acid lysine were increased in toc1-2 ( Supplementary  Figures 3A,B). Thus, proper expression of TOC1 is important for the diel and circadian oscillations of sugars, amino acids and TCA intermediates.
We next examined whether a main output of metabolism, i.e., energy production in the form of ATP, was altered when the clock is not properly functioning. We used a highly sensitive method to quantify total adenine nucleotides (Zhang et al., 2020). Our results showed a clear variation in adenosine diphosphate (ADP) and adenosine triphosphate (ATP) content during the day compared to the night in WT plants (Figures 4A,B). The amplitude of the oscillation was considerably reduced in TOC1 miss-expressing plants suggesting that proper TOC1 expression is necessary for sustaining diurnal energy homeostasis (Figures 4A,B). Calculation of the ATP/ADP ratio, a measure of the energetic homeostasis, showed that the day-night differences in WT plants were reduced or lost in TOC1 miss-expressing plants ( Figure 4C). Consistent with the constitutive overexpression, the ATP/ADP ratio was reduced in TOC1-ox at both time points, whereas the effect of toc1-2 mutation was more evident during the night (Figure 4C). The results connect the circadian function of TOC1 with metabolism and cellular energy.
Our results and those previously reported (Fukushima et al., 2009) suggest a role for the PRR protein family modulating the rhythms of primary metabolites. However, it is unknown whether this connection is direct or occurs through other molecular clock components and/or clock outputs. As TOC1 functions as a general repressor (Gendron et al., 2012;Huang et al., 2012;Pokhilko et al., 2012), we compared metabolites miss-regulated in TOC1-ox and toc1-2 relative to WT, and focused on those downregulated in TOC1-ox and up-regulated in toc1-2. We reasoned that metabolites in this cluster might provide an indication of the possible direct regulation by TOC1 of genes involved in the generation of these metabolites. Our analyses uncovered the sugar fructose, the amino acid glycine, and the TCA intermediate fumaric acid as down-regulated in TOC1-ox and up-regulated in toc1-2 under both LD (Supplementary Figure 3C, dotted orange oval) and LL (Figure 3C, dotted orange oval).

TOC1 Controls the Diurnal and Circadian Expression of FUM2 and Directly Binds to Its Promoter
Fumaric acid is an intermediary of the TCA cycle synthesized by the activity of two fumarases (FUM1 and FUM2). As the diurnal and circadian oscillation of fumaric acid dampened low in TOC1-ox and high in toc1-2 ( Figure 5A and Supplementary Figure 4A), we focused on the possible regulation of the FUM genes by TOC1. Time course analyses showed that FUM1 expression followed a low amplitude oscillation ( Figure 5B) that was not manifestly altered in TOC1-ox (Supplementary Figure 4B). In contrast, FUM2 expression clearly oscillated with a peak of expression in the middle of the day under different photoperiods and under LL conditions (Supplementary  Figures 4C,D). The amplitude of FUM2 expression was reduced in TOC1-ox and increased in toc1-2 under entraining conditions ( Figure 5C). Under LL, FUM2 circadian peak of expression was abolished in TOC1-ox ( Figure 5D) and upregulated during the subjective day in toc1-2 (Supplementary Figure 4D). These results suggest that the transcriptional effects of TOC1-ox and toc1-2 correlate with the phenotypes of fumaric acid accumulation.
Our previous chromatin immunoprecipitation massive sequencing (ChIP-Seq) analyses have shown that FUM2 is one of the targets of TOC1 (Huang et al., 2012). To examine whether regulation of FUM2 gene expression by TOC1 occurs through direct binding to the FUM2 locus, we performed ChIP assays followed by Quantitative PCR (QPCR). Our assays using TOC1-ox plants showed amplification of the FUM2 promoter to a degree similar to that observed for a positive control (CCA1 promoter) (Figure 5E and Supplementary Figure 4E). The promoters of other metabolism-related genes (ACO1, DIC1, NDA1, and PDK), a negative control (ACTIN) or samples processed without antibody (-α) did not show significant enrichment under our experimental conditions (Figure 5E and Supplementary Figure 4E). We also performed ChIP assays with TOC1 Minigen (TMG) seedlings expressing the TOC1 genomic fragment under TOC1 promoter (Más et al., 2003a). Our analyses revealed a rhythmic binding with peak enrichment   around ZT15 (Figure 5F). Binding was observed in a FUM2 promoter region containing circadian-related motifs including the Evening Element-like motif (EEL and AATATCT), previously identified as a TOC1-binding motif, CCA1 Binding Site (CBS and AAAAATCT) and Morning Element (ME and CCACAC) (Michael et al., 2008) (Figure 5E). Together, the results are consistent with a direct binding of TOC1 to the FUM2 promoter to regulate its diurnal and circadian transcriptional expression that ultimately also correlates with fumaric acid abundance.

Over-Expression of FUM2 Restores the Metabolism-Related Phenotypes of TOC1-ox
If the phenotypes observed in TOC1-ox are due to the repression of FUM2, its over-expression in the TOC1-ox background should decrease the severity of the phenotypes. Thus, we next conducted genetic interaction studies by transforming TOC1-ox plants with a FUM2 over-expressing construct (FUM2-ox)  Figures 5A,B). The severity of the rosette size and fresh weight phenotypes of TOC1-ox was also significantly improved in the double FUM2/TOC1-ox plants compared to single TOC1-ox (Figures 6A,B and Supplementary Figure 5C). A similar phenotypic recovery was observed in two different double overexpressing lines ( Figure 6C and Supplementary Figure 5C). These results suggest that the reduced expression of FUM2 in TOC1-ox plants contributes to the observed phenotypes and that over-expression of FUM2 alleviates these phenotypes. The results also limit the possibility that the observed phenotypes of TOC1-ox plants are due to indirect pleiotropic effects.

(Supplementary
Growing conditions that mimic energy deprivation trigger alternative metabolic pathways including the induction of genes involved in the catabolism of Branched-Chain Amino Acids (BCAA) (Pedrotti et al., 2018). Notably, the expression of a number of amino acid catabolic genes was up-regulated in TOC1-ox during the night (Figures 7A,D,G). The up-regulation was more evident under shorter photoperiods (Supplementary Figure 5A). The genes were up-regulated in TOC1-ox although the plants were grown under normal cycles and not under energy deprivation conditions in which these genes are induced. Other genes encoding amino acid catabolic enzymes as well as sugar starvation markers were also induced in TOC1ox plants (Supplementary Figures 6B,C). Conversely, gene expression was down-regulated in FUM2-ox and in toc1-2 (Figures 7A,B,D,E,H). The up-regulation of BCCA catabolic gene expression in TOC1-ox was overcome in the double FUM2/TOC1-ox plants (Figures 7C,F,I), suggesting that the observed phenotypes in TOC1-ox are indeed mediated by the reduced expression of FUM2. The results support a redirected proteolytic metabolism due to the altered energy status of TOC1 miss-expressing plants.

DISCUSSION
Temporal compartmentalization of metabolism is pervasive in unicellular and multicellular organisms. Due to their sessile nature and highly compartmentalized cells (Sweetlove and Fernie, 2013), an accurate spatiotemporal distribution of metabolic processes is particularly essential for plants. Clockdependent mechanisms in cooperation with specific regulatory pathways allow plants to actively anticipate and prepare reactions in sync with internal metabolic signals and external environmental cues (Smith and Stitt, 2007;Sanchez-Villarreal et al., 2013). We have identified a molecular mechanism that connects the circadian clock with the regulation of metabolism in Arabidopsis thaliana.
Our results show that proper expression of TOC1 is important for the diel and circadian oscillations of sugars, amino acids and TCA intermediates. TCA intermediates, including fumarate were also found to accumulate in the triple d975 (PRR9, PRR7, and PRR5) mutant plants (Fukushima et al., 2009). Compared to TOC1-ox, the phenotypes of toc1-2 plants were less severe, suggesting a possible functional redundancy with other members of the PRR family (for instance PRR5). Overall, The white and gray areas represent day and night, respectively. Some data are repeated in different graphs to facilitate the comparisons among genotypes. Two (FUM2-ox) and three (WT, TOC1-ox, toc1-2, and FUM2/TOC1-ox) biological replicates were performed, with plants grown independently, with samples collected, processed and analyzed at different times.
the metabolite profiles were comparable to those previously reported (Flis et al., 2019) with some differences that might be explained by the different plant developmental stages and growing conditions. Our results show that TOC1 regulates the expression of FUM2. Consistently, FUM2 expression was down-regulated by induction of REVEILLE8 (Hsu et al., 2013), a previously described activator of TOC1 (Hsu et al., 2013;Xing et al., 2015;Ma et al., 2018). The altered fumarate accumulation in TOC1 miss-expressing plants suggests that the transcriptional regulation of FUM2 is responsible for the observed changing patterns of fumarate. A similar conclusion can be drawn by our genetic interaction studies using FUM2ox plants. As fumarate and starch serve as alternative carbon sinks for photosynthate (Chia et al., 2000;Tschoep et al., 2009;Pracharoenwattana et al., 2010), the regulation of fumarate accumulation (and other metabolites) by TOC1 provides a novel way for the circadian clock to control cellular metabolism in addition to the one previously described through the control of starch metabolism.
We also found that TOC1 directly binds to the FUM2 promoter to regulate its diurnal and circadian transcriptional expression. Although previous studies have provided evidence that multiple circadian clock outputs regulate diel patterns of C and N metabolism, the binding of TOC1 to the FUM2 promoter and the regulation of FUM2 expression provide a novel direct molecular mechanism connecting TOC1 with circadian metabolism. Our results also suggest that under normal LD cycles, the energy demand in TOC1 miss-expressing plants is not properly supplied so that alternative pathways for generating energy are triggered. Many of the enzymes of the BCAA pathway are located in mitochondria (Taylor et al., 2004) and the BCAA genes display rhythms during the diurnal cycle (Blasing et al., 2005;Lee et al., 2010), being transcriptionally repressed by light (Ishizaki et al., 2005(Ishizaki et al., , 2006. It is plausible that the disruption of energy homeostasis in TOC1-ox triggers an alarm response to obtain energy from other sources, in this case the catabolism of branched amino acids. Analyses of Arabidopsis tic mutant plants showed an altered starch, carbohydrate, amino acid and fumarate accumulation. Other TCA cycle intermediates were not significantly changed in the tic mutant (Sanchez-Villarreal et al., 2013) suggesting that the phenotypes are due to an altered balance of carbon and nitrogen usage.
In mammals, the circadian regulation of metabolism relies on both systemic as well as local signals (Brown, 2016) affecting basic pathways at a cellular level and metabolic homeostasis within the whole organism (Asher and Schibler, 2011). In plants, the mitochondrial proteome varies depending on the metabolic requirements of the different tissues (des Francs-Small et al., 1992;Bardel et al., 2002;Lee et al., 2008). Based on the circadian communication between shoots and roots in Arabidopsis (James et al., 2008;Takahashi et al., 2015;Chen et al., 2020), it would be interesting to dissect the possible cell-, tissue-, and organspecificity of metabolic regulation by the clock, and identify possible systemic signal(s) able to transmit energy status among different parts of the plant for improved productivity and fitness. Furthermore, metabolism-related organelles do not function in isolation within the cell but are tightly connected with other key cellular pathways including the cell cycle (Kianian and Kianian, 2014). As the circadian clock controls the timing of the cell cycle in plants (Fung-Uceda et al., 2018), it is possible that the circadian interaction with the cell cycle might also play a role regulating cellular metabolism in plants. Our results pave the way for further studies exploiting the chronobiology of metabolism to improve metabolic homeostasis and cellular energy in plants and animals.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

AUTHOR CONTRIBUTIONS
LC-C, TY, YZ, and MO performed the experiments. AF contributed with reagents, ideas, and comments. PM conceived the project, designed the experiments, and wrote the manuscript. All authors read, revised, and approved the manuscript.

FUNDING
The Mas laboratory is funded from the FEDER/Spanish Ministry of Economy and Competitiveness (PID2019-106653GB-I00/AEI/10.13039/501100011033), from the Ramon Areces Foundation and from the Generalitat de Catalunya (AGAUR) (2017 SGR 1211). PM laboratory also acknowledges financial support from the CERCA Program/Generalitat de Catalunya and by the Spanish Ministry of Economy and Competitiveness through the "Severo Ochoa Program for Centers of Excellence in R&D" (CEX2019-000902-S). MO is a recipient of a "Severo Ochoa" Internationalization Postdoctoral Program fellowship. This work was also supported by funding from the Max-Planck Society (AF, TY, and YZ) and the European Union's Horizon 2020 project PlantaSYST .