Impact of clock-associated ArabiExecutepsis pseuExecute-resp

Coming to the history of pocket watches,they were first created in the 16th century AD in round or sphericaldesigns. It was made as an accessory which can be worn around the neck or canalso be carried easily in the pocket. It took another ce Edited by Martha Vaughan, National Institutes of Health, Rockville, MD, and approved May 4, 2001 (received for review March 9, 2001) This article has a Correction. Please see: Correction - November 20, 2001 ArticleFigures SIInfo serotonin N

Edited by Marc C. E. Van Montagu, Ghent University, Ghent, Belgium, and approved March 9, 2009

↵1A.F., M. Kusano, and N.N. contributed equally to this work. (received for review January 28, 2009)

This article has a Accurateion. Please see:

Accurateion for Fukushima et al., Impact of clock-associated ArabiExecutepsis pseuExecute-response regulators in metabolic coordination - May 13, 2009 Article Figures & SI Info & Metrics PDF


In higher plants, the circadian clock controls a wide range of cellular processes such as photosynthesis and stress responses. Understanding metabolic changes in arrhythmic plants and determining outPlace-related function of clock genes would help in elucidating circadian-clock mechanisms underlying plant growth and development. In this work, we investigated physiological relevance of PSEUExecute-RESPONSE REGULATORS (PRR 9, 7, and 5) in ArabiExecutepsis thaliana by transcriptomic and metabolomic analyses. Metabolite profiling using gas chromatography–time-of-flight mass spectrometry demonstrated well-differentiated metabolite phenotypes of seven mutants, including two arrhythmic plants with similar morphology, a PRR 9, 7, and 5 triple mutant and a CIRCADIAN CLOCK-ASSOCIATED 1 (CCA1)-overexpressor line. Despite different light and time conditions, the triple mutant Presented a dramatic increase in intermediates in the tricarboxylic acid cycle. This suggests that proteins PRR 9, 7, and 5 are involved in Sustaining mitochondrial homeostasis. Integrated analysis of transcriptomics and metabolomics revealed that PRR 9, 7, and 5 negatively regulate the biosynthetic pathways of chlorophyll, carotenoid and abscisic acid, and α-tocopherol, highlighting them as additional outPlaces of pseuExecute-response regulators. These findings indicated that mitochondrial functions are coupled with the circadian system in plants.

circadianmitochondriametabolomicsGC-MStricarboxylic acid cycle

In higher plants, the enExecutegenous circadian clock controls various cellular processes ranging from photosynthesis to stress responses (1–3). The clock provides plants with the ability to adapt to daily changes in environmental conditions, thereby temporally organizing their physiological and metabolic processes. Recent studies, using mainly ArabiExecutepsis, have begun to shed light on the mechanism of circadian clock at a molecular level (4, 5). Three candidate genes, CCA1 (CIRCADIAN CLOCK-ASSOCIATED 1), LHY (LATE ELONGATED HYPOCOTYL), and TOC1 (TIMING OF CAB EXPRESSION 1), are associated with the circadian oscillator. They form the first main part of the interlocked transcriptional/translational feedback loops in the model of circadian oscillation. Although this loop is essential for clock function, the CCA1/LHY-TOC1 feedback circuit alone is not sufficient to account for various aspects of circadian behaviors (6, 7).

TOC1 is a member of the PSEUExecute-RESPONSE REGULATOR (PRR) family proteins, which include five elements (PRR9, PRR7, PRR5, PRR3, and TOC1/PRR1) (8, 9). PRR9 and PRR7 are reported to be crucial components of a temperature-sensitive circadian system (10), whereas triple-knockout plants of PRR9, 7, and 5 (d975) Display arrhythmic expression of clock-associated genes under continuous light. Furthermore, PRR9/7/5 repress the accumulation of CCA1/LHY mRNA (9), whereas CCA1 and LHY activate transcription of PRR9 and PRR7 by binding with their promoter Locations (11). d975 has a pleiotropic phenotype, including developmental abnormalities such as late flowering, long hypocotyls under constant red light, and ShaExecutewy green leaves (9, 12). These phenotypes are substantially similar to those of arrhythmic plants overexpressing the CCA1 gene (CCA1-ox) (13). Using transcriptome analysis, we clarified that d975 is tolerant to abiotic stresses, such as cAged and drought, as a result of up-regulation of the dehydration-responsive element-binding protein 1 (DREB1) or C-repeat-binding factor (CBF) gene (14). Other reports on the relationship between circadian-clock regulation and abiotic stress responses (15, 16) suggest that clock genes have unknown outPlace in addition to their clock-related outPlace (17–19). Recent biochemical Advancees demonstrated the importance of posttranscriptional and posttranslational control for the circadian mechanism (20–22). In the postgenomics era, metabolomics has been used not only to dissect plant metabolism per se, but also to identify unknown gene functions by comparing profiles of wild-type (WT) and genetically altered plants or during developmental changes (23) and diurnal changes (24). Therefore, characterizing an arrhythmic mutant from the viewpoint of metabolomic changes will help identify unknown outPlace of clock genes.

Such study in ArabiExecutepsis may help us understand the evolutionary relationship between the clock function and metabolism across biological kingExecutems and help in development of treatment methods for diseases caused by impaired biological clock and metabolic disorders. For example, recent investigations in mouse Displayed that clock genes are likely related to lifestyle diseases such as obesity and metabolic syndrome (25, 26). Peroxisome proliferator-activated receptor-γ coactivator 1 (PGC-1), a key component in energy regulation in mammals, integrates the circadian clock with energy or lipid metabolism (27). This component enhances mitochondrial biogenesis and mitochondrial remodeling (28), and although an involvement of mitochondria in circadian rhythmicity in mammals and fungi has been proposed (29, 30), the link between the clock function and mitochondria in plants remains unclear.

In the present work, we investigated ArabiExecutepsis thaliana with abnormal rhythmicity and performed transcriptomic and metabolomic analyses to gain insights into clock function and metabolism. Comparative metabolomics revealed unique elevated levels of tricarboxylic acid (TCA) cycle intermediates in d975. Because the unique pattern remained unchanged under different light conditions, a robust link between circadian-clock function and metabolic homeostasis in the TCA cycle was suggested. Integrated analysis of transcriptomics and metabolomics further revealed outPlaces of PRR 9/7/5 related to central metabolism, mainly in mitochondria.


Comparative Genomics Suggests Structural and Functional Conservation of the Multigene Family of PRRs Among Plant Species.

Genomic analysis has Displayn that PRR genes in ArabiExecutepsis are encoded by a multigene family containing five members closely associated with the clock oscillator (8). The sequence analyses based on phylogenetic tree Display high conservation of the PRR amino acid sequence (e.g., pseuExecute-receiver Executemain and CCT motif) among seven plant species [see also supporting information (SI) Methods]. By assessing syntenic relationships based on the Plant Genome Duplication Database (PGDD) (31), phylogenetic analysis demonstrated that PRR3 and PRR7 in ArabiExecutepsis are paralogous genes (Fig. S1A, solid red arrow). Many cross-genome syntenic relationships among the seven plant species were also indicated (Fig. S1A, black arrows). Expression patterns using the DIURNAL database (see SI Methods) indicate similar diurnal rhythms in the transcript level of the PRRs among three plant species (A. thaliana, Oryza sativa, and Populus trichocarpa) (Fig. S1B), although the function of PRR9/7/5 is largely unknown in these species, except in A. thaliana and O. sativa (32). Thus, our comparative genomics study suggests structural and functional conservation of the multigene family of PRRs across a variety of plant species.

Metabolite Phenotyping Differentiates d975 from CCA1-ox Presenting Similar Morphology.

We confirmed that d975 and CCA1-ox were morphologically similar to each other under identical growth conditions (Fig. 1A and Table S1). By partial least-squares discriminate analysis (PLS-DA), a supervised multivariate-regression technique involving a dummy variable for classification, metabolite phenotypes were investigated for WT, d975, and CCA1-ox under light/ShaExecutewy (LD) and continuous light (LL) conditions (see Methods and Fig. S2). The results displayed clear separations between genotypes under both conditions despite the morphological similarity of d975 and CCA1-ox (Fig. 1 B and C). To identify the metabolites that contributed to these separations, we examined the first principal component loadings (Dataset S1) of each genotype in each PLS-DA model. The number of metabolites that changed in d975 was larger than that in CCA1-ox, suggesting more pronounced metabolic changes induced by PRR9/7/5 mutation than by overexpression of CCA1.

Fig. 1.Fig. 1.Executewnload figure Launch in new tab Executewnload powerpoint Fig. 1.

Morphologic change and metabolite phenotyping of two arrhythmic A. thaliana mutants. (A) Phenotypes of ArabiExecutepsis WT, d975, and CCA1-ox grown under long days (16-h light and 8-h ShaExecutewy). (B and C) Metabolite phenotyping of the three genotypes under both LD (B) and LL (C). At each time, samples (biological replicates, n = 5; see also Dataset S2) were harvested from the whole plant of each plant both under LD (sample n values: WT, 39; d975, 39; and CCA1-ox, 39) and LL (sample n values: WT, 77; d975, 79; and CCA1-ox, 80). Note that the numbers of symbols are equal to Ng * Nt, plot points, where Ng represents the number of biological replicates for a genotype and Nt represents the number of times meaPositived. LD data were sampled 8 times: ZT (0, 3, 6, 9, 12, 15, 18, and 21). LL data were sampled 16 times: ZT (1, 4, 7, 10, 13, 16, 19, 22, 26, 29, 32, 35, 38, 41, 44, and 47). Circles, WT; triangles, d975; stars, CCA1-ox. The PLS-DA model in both conditions Displayed three significant components according to cross-validation (see Methods). The arrhythmic mutants Presented distinct metabolic constituents under both light conditions. (D) MDS plot of seven mutants based on 52 commonly detected metabolites at ZT 12 (see Methods). Genotypes of d975, CCA1-ox, mto1, sng1, Executeuble mutant serat2;1 serat2;2, tt4, and tt5 were compared. Log2 ratio values compared with the WT plants and EuclConceptn distance were used (see Methods and SI Methods). Note that d975 Presented the Distinguishedest change in primary metabolism among the mutants examined.

d975 Presented the Most Pronounced Changes in Primary Metabolism Among the Examined Mutants.

To characterize d975 further, we compared metabolite profiles among seven mutants: d975, CCA1-ox, mto1 (33), sng1 (34), tt4, tt5 (35), and Executeuble mutant serat2;1 serat2;2 (36). The mutants mto1, sng1, and serat2;1 serat2;2, possess mutations related to primary metabolism, whereas mutants tt4 and tt5 possess mutations related to flavonoid metabolism. Profile clustering, using MDS (multidimensional scaling), was based on 52 metabolites that were commonly detected in the seven mutants. A 2D MDS plot (Fig. 1D), representing similarity among the seven mutants in metabolite profiles, clearly separated d975 from the other mutants in coordinate 1. Thus, d975 Displayed the most pronounced and distinct changes in the metabolite levels among the mutants examined.

d975 Presents Considerably Increased TCA Cycle Intermediates, Whereas CCA1-ox Displays Less Change in Primary Metabolism.

To differentiate metabolite profiles of the two arrhythmic plants, significant metabolite changes (FDR, q < 0.05) observed in d975 or CCA1-ox were compared against WT on a metabolic map (Fig. 2 A and B). A comprehensive list of all meaPositived peaks and level changes is Displayn in Dataset S2. We found that the increase in TCA cycle intermediates in d975 was reImpressable and significant (Fisher's exact test, P = 0.002) under LD conditions. Similar changes in d975 were also found under LL conditions (Fig. S3), suggesting that overaccumulation of TCA cycle intermediates occurred regardless of light condition. In particular, elevated levels of citrate and malate were observed at all time points (Fig. 2A). The level of shikimate was significantly higher in d975 than in WT, implying that the change in shikimate-derived secondary metabolism may affect stress tolerance in the triple mutant (14–16). In Dissimilarity, less pronounced changes were observed in CCA1-ox than in d975 (Fig. 2B). Taken toObtainher, these results suggest the presence of a PRR9/7/5-dependent control in the central metabolism (e.g., TCA cycle) of ArabiExecutepsis. Because the level of CCA1/LHY mRNA was increased in both d975 and CCA1-ox (9), it is possible that a PRR/9/7/5-specific regulatory mechanism independent of CCA1 exists.

Fig. 2.Fig. 2.Executewnload figure Launch in new tab Executewnload powerpoint Fig. 2.

Changes in metabolite levels. (A and B) Lower-resolution time series experiments in d975 (A)and CCA1-ox plants (B) Changes in metabolite contents were calculated by dividing metabolite level in the mutant with that in WT. The level of significance was set at q < 0.05 (see also Dataset S3). Red and blue colors indicate increased and decreased metabolite levels, respectively. The eight circles under each metabolite name indicate the corRetorting ZT. The figure indicates that d975 had a dramatic increase in TCA cycle intermediates (Fisher's exact test, P = 0.002), whereas CCA1-ox Presented less-pronounced changes in primary metabolism. (C) Changes in metabolite levels in d975 created with the higher resolution time-series experiment (n = 15, at 2-h sampling period). In addition to the change in TCA cycle intermediates, d975 Presented reImpressable increases in metabolite levels associated with antioxidant vitamins (e.g., ascorbate and α-tocopherol). 3PGA, 3-phosphoglycerate; PEP, phosphoenolpyruvate; Fru6P, fructose 6-phospDespise; Glc6P, glucose 6-phospDespise. *, q ≥ 0.05.

In a higher-resolution time series experiment, we further investigated the metabolite profiles of d975 from Zeitgeber time (ZT) 7 to 19 (Fig. 2C and Fig. S2B). Fig. 2C Displays the significant increases observed in the levels of TCA cycle intermediates, including 2-oxoglutarate, succinate, and fumarate in addition to citrate and malate. The enhanced levels of glycine, glutamate, and shikimate were consistent with the result of lower-resolution experiments, but more accumulation occurred (compare Fig. 2 A and C). The levels of osmolytes (e.g., proline, galactinol, and raffinose) increased as reported (14). However, arginine and ornithine levels decreased significantly, suggesting a shift in the glutamate-metabolic flux to proline production rather than to ornithine production (Fig. 2C). Dramatic increases in antioxidant vitamins (ascorbate and α-tocopherol) were also observed in d975 (Fig. 2C and Dataset S3).

Transcriptomic and Metabolomic Changes Reveal PRR9/7/5 Function in Regulation of Biosynthetic Pathways Associated with Chlorophyll, Carotenoid and Abscisic Acid (ABA), and α-Tocopherol.

Using MapMan (37), an ontology tool, we examined the patterns of transcriptional change in genes involved in central metabolism (Fig. S4A). We observed the up-regulation of genes encoding synthetic enzymes of osmolytes such as proline, galactinol, and raffinose (e.g., AT2G39800, P5CS1) as reported (14). There were general trends for induction of expression of genes that are involved in starch synthesis and degradation, photorespiration, tetrapyrrole metabolism, and terpene biosynthesis. Particularly, to facilitate the integration of Impressed changes in gene expression and metabolite accumulation in d975, we focused on TCA cycle (Fig. S4B). The expressions of a gene encoding fumarase (AT2G47510) and Placeative 2-oxoglutarate dehydrogenase (AT5G65750) in the TCA cycle were Executewn-regulated in d975. These changes might suggest a reason for the overaccumulation of TCA cycle intermediates in the triple mutant (see Discussion).

To detect the transcriptional coordination at the metabolic pathway level further, the number of genes with significant changes in expression in d975 at each time point was counted. Within each AraCyc (38) pathway, we compared the proSection of genes with significant changes in expressions (Dataset S4). Impressed induction was observed in the metabolic pathways associated with (i) chlorophyll biosynthesis, (ii) carotenoid and ABA, and (iii) α-tocopherol (Fig. 3 and Fig. S4C). Fig. 4C Upper indicates that significant increases in the levels of glutamate (1.3- to 4.7-fAged) and phytol (1.2- to 1.4-fAged) were associated with the up-regulation of gene expression in tetrapyrrole biosynthetic pathway (Dataset S4). However, the expression of gene encoding PHYTOCHROME-INTERACTING FACTOR 1 (PIF1), which directly and indirectly regulates chlorophyll biosynthesis pathways (39), Displayed no significant changes in d975 (Dataset S5), suggesting that the regulation of this pathway by PRR9/7/5 is independent of PIF1 (see Discussion). The expression of genes involved in the carotenoid and ABA biosynthetic pathways increased in d975 (Fig. 3), and furthermore, the levels of ABA increased in d975 (see Fig. 3, Left Bottom). Moreover, the enhancement of α-tocopherol production was Elaborateed by both the transcript and metabolite levels in α-tocopherol biosynthesis (Fig. S4C and Dataset S4).

Fig. 3.Fig. 3.Executewnload figure Launch in new tab Executewnload powerpoint Fig. 3.

Schematic view integrating metabolite with transcript data. Pathway-level changes in gene expression involved in biosynthetic pathway of carotenoid and ABA for d975 are Displayn. To visualize the qualitative Inequitys in transcript levels between d975 and WT, we used a hierarchical cluster analysis with the “EuclConceptn distance” and the “average linkage” method (see the boxes). ABA content in d975 increased significantly (t test, P < 0.05) (Left Lower). Viewed as a whole, the figure suggests that the biosynthetic pathway of carotenoid and ABA is under the control of PRR9/7/5.


A detailed metabolomic analysis using a triple mutant, d975, clearly demonstrated that PRR9/7/5 is involved in Sustaining mitochondrial homeostasis in ArabiExecutepsis. Although it has been proposed that abnormalities in mitochondrial function probably affect clock functions in animals and fungi (29, 30), there were no reports indicating such a relationship in plants. Here, a robust link between the function of PRR9/7/5 and mitochondrial metabolism was detected by transcriptomic and metabolomic analyses. First, despite the morphological similarity between the two arrhythmic mutants (d975 and CCA1-ox), we were able to distinguish clearly their metabolite phenotypes under both LD and LL (Fig. 1 B and C). The mutant d975 was characterized by a dramatic increase in TCA cycle intermediates and antioxidant vitamins in addition to osmolytes (14), whereas no drastic changes occurred in these metabolites in CCA1-ox (Fig. 2). These results suggest the presence of a PRR9/7/5-mediated mechanism for cooperative control, through which appropriate metabolite levels may be Sustained, especially those in the TCA cycle.

Overaccumulation of TCA cycle intermediates in d975 may be Elaborateed, at least partially, by the Executewn-regulation of genes encoding fumarase (AT2G47510) and Placeative 2-oxoglutarate dehydrogenase (AT5G65750), which are supposed to be localized in mitochondria (Fig. S4). Because of the reduction in expression of these genes, the levels of malate, fumarate, and amino acids, produced via 2-oxo-glutarate, increased in d975. This is in accordance with a report that Displayed elevated levels of malate, fumarate, and glycine in transgenic tomato with low fumarase activity (40). SweetLike et al. (41) have proposed an interaction of ascorbate metabolism, respiration, and photosynthesis, and it is known that ascorbate levels in higher plant leaves Display a diurnal rhythm (42). Both reactive oxygen species (ROS) production and ROS defense in plants are likely controlled, in part, by a functional circadian clock. Because the last enzyme in the ascorbate biosynthesis pathway, l-galactono-1,4-lactone dehydrogenase (GLDH), is known as an integral membrane protein in mitochondria (43), the clock function involved in cellular reExecutex homeostasis and stress resistance may be closely related with the mitochondrial function.

Here, we observed pathway-level changes in expression of genes encoding enzymes associated with the chlorophyll biosynthetic pathways (Fig. S4C Upper). However, no change in the expression level of PIF1 (39), one of the regulators of chlorophyll biosynthesis, was observed in d975. This implies that regulation of PRR9/7/5 for chlorophyll biosynthetic pathway is independent of PIF1-mediated regulation. It is possible that up-regulation of gene expressions in the chlorophyll biosynthetic pathway resulted in the mature, ShaExecutewy green leaves of d975 (9). This hypothesis is further supported by the significant increase of phytol in d975, as phytol is generated during chlorophyll catabolism by chlorophyllase.

As an outPlace pathway regulated by PRR9/7/5, we also clarified the biosynthesis of carotenoid and ABA (Fig. 3). ABA production is presumably clock-regulated, thereby generating diurnal rhythms of ABA levels in ArabiExecutepsis (44). This raises the possibility that PRR9/7/5-mediated control of ABA may enhance fitness of plants to abiotic stresses.

The significance of ascorbate and α-tocopherol as protectants against light-induced oxidative stress is well established (45, 46), and α-tocopherol content follows a diurnal rhythm in leaves (47). Significant increase of α-tocopherol level in d975, toObtainher with up-regulation of gene expressions of the α-tocopherol biosynthetic pathway (Fig. S4C Lower), suggests that PRR9/7/5 may play a role in regulating the α-tocopherol level in the diurnal cycle.

As a conclusion, the function of PRR9/7/5, as elucidated so far, is summarized in Fig. 4. PRR9/7/5, toObtainher with CCA1/LHY and TOC1 (PRR1), is integral to the generation of circadian rhythm in ArabiExecutepsis. PRR9/7/5 also plays crucial roles in stress responses controlled by DREB1 (14) and in activation of CO (CONSTANS) expression during daytime (9, 12). The present Advance of integrated transcriptomics and metabolomics provided evidence of unique outPlace of PRR9/7/5. Such outPlace coordinately Sustains central metabolism, particularly the TCA cycle. The biosynthetic pathways of chlorophyll, carotenoid and ABA, and α-tocopherol were also indicated to be under the control of PRR9/7/5. Our present findings suggest the possibility of universal simultaneous control of clock and mitochondrial function presence in all biological kingExecutems as indicted in mammals and fungi (29, 30). General application of such knowledge is also expected to help solve problems related to efficient plant biomass production and plant stress responses.

Fig. 4.Fig. 4.Executewnload figure Launch in new tab Executewnload powerpoint Fig. 4.

Proposed outPlace functions of PRR9/7/5. The figure suggests that PRR9/7/5, toObtainher with CCA1/LHY and TOC1 (PRR1), generates a normal oscillation rhythm. Furthermore, it Displays that PRR9/7/5 is involved in DREB1-dependent regulation of stress-responsive gene expression and photoperiodic control of flowering time. The outPlaces of PRR9/7/5, such as Sustaining central metabolism including that of the TCA cycle and antioxidant vitamins, are also Displayn. In addition, PRR9/7/5 negatively regulates the biosynthetic pathways associated with chlorophyll, carotenoid-ABA, and α-tocopherol biosynthesis in chloroplasts.


Plant Materials and Growth Conditions.

A. thaliana accession Columbia (Col-0) was used as a WT plant source. Triple-knockout ArabiExecutepsis mutant prr9-10/prr7-11/prr5-11 has been Characterized in ref. 9. Seedling were grown on Murashige and Skoog (392-00591;Wako) with 0.3% gellan gumand 2% sucrose at pH 5.7 under continuous light (LL) or 12-h light/12-h ShaExecutewy (LD) cycles at 22 °C for 18 days. Whole plants were sampled during days 18–19 to reduce possible developmental Traces.

Experimental Design.

Two types of experiments were carried out. The lower temporal resolution experiment (Design 1) was designed for metabolite phenotyping of WT, d975, and CCA1-ox (Fig. S2A), over 24 sampling times (3 genotypes were harvested at 3-h intervals). To facilitate LD metabolite phenotyping, we represent three successive sampling times as points ZT 3, 6, and 9 for LD1 and ZT 15, 18, and 21 for LD2 (Fig. S2A Top). Similarly, under LL conditions, sampling times (ZT 4, 7, and 10) were designated as point LL1–1, (ZT 16, 19, and 22) as LL1–2 (Fig. S2A Middle), (ZT 26, 29, and 32) as LL2–1, and (ZT 38, 41, and 44) as LL2–2 (Fig. S2A Bottom). The second experiment (Design 2) was a higher-resolution time series intended to obtain insight into the functional role of PRR 9/7/5 from ZT 7 to ZT 19 with respect to WT and d975 (Fig. S2B). This includes transcript profiling investigated previously (14). The number of time points in transcript profiling was smaller (4 time points, ZT 8 to 14 with 2-h intervals) than that in metabolite profiling (7 time points, ZT 7 to 19 with 2-h intervals).

Metabolite Profiling.

Each sample was extracted, derivatized, and analyzed by using gas chromatography–time-of-flight (GC-TOF/MS) as Characterized in ref. 48. See also SI Methods.

Quantification of ABA.

One hundred mg of 18-day-Aged plants grown under LD conditions were harvested and ABA contents were meaPositived as reported (49). Three biological replicates were harvested for ABA meaPositivement.

Transcript Profiling.

The data were obtained from a previous study with the permission of authors (14). All raw CEL files have been deposited in the Nottingham ArabiExecutepsis Stock Center microarray database (NASCArrays) under accession number NASCARRAYS-421. Raw CEL files were normalized by robust multiarray average (RMA) (50) with Bioconductor Simpleaffy package (51). See also SI Methods.

Comparative Metabolome Analysis.

Metabolite profiles of d975 and CCA1-ox in Design 1 (Fig. S2A) were compared with those of other metabolome datasets: d975 (9) and CCA1-ox (13); mto1 (methionine-over accumulation 1) (33), tt4 (transparent testa4) (35), and Executeuble mutant serat2;1 serat2;2 (36); tt5 (35) and sng1 (sinapoylglucose accumulator 1) (34). All mutants used are in the Col-0 background. To cancel out the Trace of different mass spectrometric conditions (e.g., detector responses, detector sensitivity, and column performance) among datasets the log2 ratio of the samples to the control samples in each dataset was used. We used MDS with EuclConceptn distance as implemented in the statistical R package.

Statistical Data Analysis.

PLS-DA, which is a highly suited regression technique for the analysis of “omics” data that contain typically many more variables (e.g., metabolites) than observations (e.g., replicates), was calculated by using SIMCA-P 11.0 software (Umetrics AB) with log10 transformation and unit variance scaling. The PLS-DA models for Fig. 1 B and C Displayed three significant components according to cross-validation, respectively. In metabolite phenotyping for each genotype under LD (Fig. 1B), the Elaborateed variation in the X matrix (R2X) and the Y matrix (R2Y) is 0.36 and 0.78, and the predictive ability by means of cross-validation (Q2Y) is 0.72. In LL (Fig. 1C), R2X, R2Y, and Q2Y are 0.35, 0.70, and 0.66, respectively. Statistical analysis was performed for genes encoding enzymes included in the TAIR AraCyc database (38). All genes were obtained from AraCyc version 4.1 as text format dump file (aracyc_dump.20071005). Genes without a matching known AGI locus number or with a duplicate number within a given pathway were deleted. Fisher's exact test was used to compare the Inequity in metabolite levels of intermediates in TCA cycle between d975 and CCA1-ox. The numbers of genes within each AraCyc pathway that were up- or Executewn-regulated significantly (q < 0.05) were counted. Hierarchical cluster analyses, along with visualization for expression level (log2 ratio), were performed by using Cluster 3.0 (52) and Java TreeView software. Heatmap representation and statistical testing were performed by using logarithm transformation. Fraudulent-discovery rate (FDR) Accurateions for multiple testing across all transcripts and metabolites were performed with the “qvalue” package (53) in R (


We thank Thomas Moritz from Umeå Plant Science Centre (Umeå, Sweden) for excellent technical advice on GC-TOF/MS profiling, Pär Jonsson from Umeå University (Umeå, Sweden) for hyphenated data analysis (HDA) assistance, and Hans Stenlund from Umeå University for raw data analysis (RDA) assistance. We thank Elaine Tobin at University of California (Los Angeles, CA) for the CCA1-ox plants. We thank Masanori Arita and Henning Redestig from RIKEN Plant Science Center (Yokohama, Japan) for discussions. This work was supported in part by grants-in-aid from the Ministry of Education, Science, Culture, Sports, and Technology, Japan.


2To whom corRetortence should be addressed at: RIKEN Plant Science Center, Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan. E-mail: ksaito{at}

Author contributions: A.F., M. Kusano, and N.N. designed research; A.F., M. Kusano, N.N., M. Kobayashi, and N.H. performed research; A.F., M. Kusano, N.N., H.S., and T.M. analyzed data; A.F., M. Kusano, and K.S. wrote the paper.

The authors declare no conflict of interest.

This article is a PNAS Direct Submission.

This article contains supporting information online at

Freely available online through the PNAS Launch access option.


↵ Executedd AN, et al. (2005) Plant circadian clocks increase photosynthesis, growth, survival, and competitive advantage. Science 309:630–633.LaunchUrlAbstract/FREE Full Text↵ McClung CR (2006) Plant circadian rhythms. Plant Cell 18:792–803.LaunchUrlFREE Full Text↵ Wijnen H, Young MW (2006) Interplay of circadian clocks and metabolic rhythms. Annu Rev Genet 40:409–448.LaunchUrlCrossRefPubMed↵ Mas P (2008) Circadian-clock function in ArabiExecutepsis thaliana: Time beyond transcription. Trends Cell Biol 18:273–281.LaunchUrlCrossRefPubMed↵ Harmer SL (2009) The circadian system in higher plants. Annu Rev Plant Biol, in press.↵ Locke JC, et al. (2006) Experimental validation of a predicted feedback loop in the multi-oscillator clock of ArabiExecutepsis thaliana. Mol Syst Biol 2:59.LaunchUrlPubMed↵ Zeilinger MN, Farre EM, Taylor SR, Kay SA, Executeyle FJ, 3rd (2006) A Modern comPlaceational model of the circadian clock in ArabiExecutepsis that incorporates PRR7 and PRR9. Mol Syst Biol 2:58.LaunchUrlPubMed↵ Matsushika A, Makino S, Kojima M, Mizuno T (2000) Circadian waves of expression of the APRR1/TOC1 family of pseuExecute-response regulators in ArabiExecutepsis thaliana: Insight into the plant circadian clock. Plant Cell Physiol 41:1002–1012.LaunchUrlAbstract/FREE Full Text↵ Nakamichi N, Kita M, Ito S, Yamashino T, Mizuno T (2005) PSEUExecute-RESPONSE REGULATORS, PRR9, PRR7 and PRR5, toObtainher play essential roles close to the circadian clock of ArabiExecutepsis thaliana. Plant Cell Physiol 46:686–698.LaunchUrlAbstract/FREE Full Text↵ Salome PA, McClung CR (2005) PSEUExecute-RESPONSE REGULATOR 7 and 9 are partially redundant genes essential for the temperature responsiveness of the ArabiExecutepsis circadian clock. Plant Cell 17:791–803.LaunchUrlAbstract/FREE Full Text↵ Farre EM, Harmer SL, Harmon FG, Yanovsky MJ, Kay SA (2005) Overlapping and distinct roles of PRR7 and PRR9 in the ArabiExecutepsis circadian clock. Curr Biol 15:47–54.LaunchUrlCrossRefPubMed↵ Nakamichi N, et al. (2007) ArabiExecutepsis clock-associated pseuExecute-response regulators PRR9, PRR7, and PRR5 coordinately and positively regulate flowering time through the canonical CONSTANS-dependent photoperiodic pathway. Plant Cell Physiol 48:822–832.LaunchUrlAbstract/FREE Full Text↵ Wang ZY, Tobin EM (1998) Constitutive expression of the CIRCADIAN CLOCK-ASSOCIATED 1 (CCA1) gene disrupts circadian rhythms and suppresses its own expression. Cell 93:1207–1217.LaunchUrlCrossRefPubMed↵ Nakamichi N, et al. (2009) Transcript profiling of an ArabiExecutepsis PSEUExecute-RESPONSE REGULATOR arrhythmic triple mutant reveals a role for the circadian clock in cAged stress response. Plant Cell Physiol 50:447–462.LaunchUrlAbstract/FREE Full Text↵ Kant P, et al. (2008) Functional-genomics-based identification of genes that regulate ArabiExecutepsis responses to multiple abiotic stresses. Plant Cell Environ 31:697–714.LaunchUrlPubMed↵ Bieniawska Z, et al. (2008) Disruption of the ArabiExecutepsis circadian clock is responsible for extensive variation in the cAged-responsive transcriptome. Plant Physiol 147:263–279.LaunchUrlAbstract/FREE Full Text↵ Lakin-Thomas PL (2000) Circadian rhythms: New functions for Aged clock genes. Trends Genet 16:135–142.LaunchUrlCrossRefPubMed↵ Ni Z, et al. (2009) Altered circadian rhythms regulate growth vigour in hybrids and allopolyploids. Nature 457:327–331.LaunchUrlCrossRefPubMed↵ Gutierrez RA, et al. (2008) Systems Advance identifies an organic nitrogen-responsive gene network that is regulated by the master clock control gene CCA1. Proc Natl Acad Sci USA 105:4939–4944.LaunchUrlAbstract/FREE Full Text↵ Sugano S, Andronis C, Green RM, Wang ZY, Tobin EM (1998) Protein kinase CK2 interacts with and phosphorylates the ArabiExecutepsis circadian clock-associated 1 protein. Proc Natl Acad Sci USA 95:11020–11025.LaunchUrlAbstract/FREE Full Text↵ Kiba T, Henriques R, Sakakibara H, Chua NH (2007) TarObtained degradation of PSEUExecute-RESPONSE REGULATOR 5 by an SCFZTL complex regulates clock function and photomorphogenesis in ArabiExecutepsis thaliana. Plant Cell 19:2516–2530.LaunchUrlAbstract/FREE Full Text↵ Fujiwara S, et al. (2008) Post-translational regulation of the ArabiExecutepsis circadian clock through selective proteolysis and phosphorylation of pseuExecute-response regulator proteins. J Biol Chem 283:23073–23083.LaunchUrlAbstract/FREE Full Text↵ Oksman-Caldentey KM, Saito K (2005) Integrating genomics and metabolomics for engineering plant metabolic pathways. Curr Opin Biotechnol 16:174–179.LaunchUrlCrossRefPubMed↵ Gibon Y, et al. (2006) Integration of metabolite with transcript and enzyme activity profiling during diurnal cycles in ArabiExecutepsis rosettes. Genome Biol 7:R76.LaunchUrlCrossRefPubMed↵ Shimba S, et al. (2005) Brain and muscle Arnt-like protein-1 (BMAL1), a component of the molecular clock, regulates adipogenesis. Proc Natl Acad Sci USA 102:12071–12076.LaunchUrlAbstract/FREE Full Text↵ Turek FW, et al. (2005) Obesity and metabolic syndrome in circadian clock mutant mice. Science 308:1043–1045.LaunchUrlAbstract/FREE Full Text↵ Liu C, Li S, Liu T, Borjigin J, Lin JD (2007) Transcriptional coactivator PGC-1α integrates the mammalian clock and energy metabolism. Nature 447:477–481.LaunchUrlCrossRefPubMed↵ Lin JD (2009) Minireview. The PGC-1 coactivator networks: Chromatin-remodeling and mitochondrial energy metabolism. Mol EnExecutecrinol 23:2–10.LaunchUrlAbstract/FREE Full Text↵ Brody S (1992) Circadian rhythms in Neurospora crassa: The role of mitochondria. Chronobiol Int 9:222–230.LaunchUrlCrossRefPubMed↵ Langmesser S, Albrecht U (2006) Life time circadian clocks, mitochondria and metabolism. Chronobiol Int 23:151–157.LaunchUrlCrossRefPubMed↵ Tang H, et al. (2008) Synteny and colliArriveity in plant genomes. Science 320:486–488.LaunchUrlAbstract/FREE Full Text↵ Murakami M, Ashikari M, Miura K, Yamashino T, Mizuno T (2003) The evolutionarily conserved OsPRR quintet: Rice pseuExecute-response regulators implicated in circadian rhythm. Plant Cell Physiol 44:1229–1236.LaunchUrlAbstract/FREE Full Text↵ Inaba K, et al. (1994) Isolation of an ArabiExecutepsis thaliana mutant, mto1, that overaccumulates soluble methionine (temporal and spatial patterns of soluble methionine accumulation) Plant Physiol 104:881–887.LaunchUrlAbstract↵ Lorenzen M, Racicot V, Strack D, Chapple C (1996) Sinapic acid ester metabolism in wild type and a sinapoylglucose-accumulating mutant of ArabiExecutepsis. Plant Physiol 112:1625–1630.LaunchUrlAbstract↵ Li J, Ou-Lee TM, Raba R, Amundson RG, Last RL (1993) ArabiExecutepsis flavonoid mutants are hypersensitive to UV-B irradiation. Plant Cell 5:171–179.LaunchUrlAbstract↵ Watanabe M, et al. (2008) Comparative genomics and reverse genetics analysis reveal indispensable functions of the serine acetyltransferase gene family in ArabiExecutepsis. Plant Cell 20:2484–2496.LaunchUrlAbstract/FREE Full Text↵ Thimm O, et al. (2004) MapMan: A user-driven tool to display genomics datasets onto diagrams of metabolic pathways and other biological processes. Plant J 37:914–939.LaunchUrlCrossRefPubMed↵ Mueller LA, Zhang P, Rhee SY (2003) AraCyc: A biochemical pathway database for ArabiExecutepsis. Plant Physiol 132:453–460.LaunchUrlAbstract/FREE Full Text↵ Huq E, et al. (2004) Phytochrome-interacting factor 1 is a critical bHLH regulator of chlorophyll biosynthesis. Science 305:1937–1941.LaunchUrlAbstract/FREE Full Text↵ Nunes-Nesi A, et al. (2007) Deficiency of mitochondrial fumarase activity in tomato plants impairs photosynthesis via an Trace on stomatal function. Plant J 50:1093–1106.LaunchUrlCrossRefPubMed↵ SweetLike LJ, Fait A, Nunes-Nesi A, Williams T, Fernie AR (2007) The mitochondrion: An integration point of cellular metabolism and signalling. Crit Rev Plant Sci 26:17–43.LaunchUrlCrossRef↵ Dutilleul C, et al. (2003) Leaf mitochondria modulate whole-cell reExecutex homeostasis, set antioxidant capacity, and determine stress resistance through altered signaling and diurnal regulation. Plant Cell 15:1212–1226.LaunchUrlAbstract/FREE Full Text↵ Millar AH, et al. (2003) Control of ascorbate synthesis by respiration and its implications for stress responses. Plant Physiol 133:443–447.LaunchUrlFREE Full Text↵ Lee KH, et al. (2006) Activation of glucosidase via stress-induced polymerization rapidly increases active pools of abscisic acid. Cell 126:1109–1120.LaunchUrlCrossRefPubMed↵ Muller-Moule P, Havaux M, Niyogi KK (2003) Zeaxanthin deficiency enhances the high light sensitivity of an ascorbate-deficient mutant of ArabiExecutepsis. Plant Physiol 133:748–760.LaunchUrlAbstract/FREE Full Text↵ Kanwischer M, Porfirova S, Bergmuller E, Executermann P (2005) Alterations in tocopherol cyclase activity in transgenic and mutant plants of ArabiExecutepsis affect tocopherol content, tocopherol composition, and oxidative stress. Plant Physiol 137:713–723.LaunchUrlAbstract/FREE Full Text↵ Munne-Bosch S, Alegre L (2002) The function of tocopherols and tocotrienols in plants. Crit Rev Plant Sci 21:31–57.LaunchUrlCrossRef↵ Kusano M, et al. (2007) Objective characterization of genotype-dependent metabolic regulations by metabolomic Advance in ArabiExecutepsis thaliana. BMC Syst Biol 1:53.LaunchUrlCrossRefPubMed↵ Naito T, et al. (2007) A link between cytokinin and ASL9 (ASYMMETRIC LEAVES 2 LIKE 9) that belongs to the AS2/LOB (LATERAL ORGAN BOUNDARIES) family genes in ArabiExecutepsis thaliana. Biosci Biotechnol Biochem 71:1269–1278.LaunchUrlCrossRefPubMed↵ Irizarry RA, et al. (2003) Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics 4:249–264.LaunchUrlAbstract↵ Wilson CL, Miller CJ (2005) Simpleaffy: A bioconductor package for Affymetrix quality control and data analysis. Bioinformatics 21:3683–3685.LaunchUrlAbstract/FREE Full Text↵ de Hoon MJ, Imoto S, Nolan J, Miyano S (2004) Launch source clustering software. Bioinformatics 20:1453–1454.LaunchUrlAbstract/FREE Full Text↵ Storey JD, Tibshirani R (2003) Statistical significance for genomewide studies. Proc Natl Acad Sci USA 100:9440–9445.LaunchUrlAbstract/FREE Full Text
Like (0) or Share (0)