<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Plant Sci.</journal-id>
<journal-title>Frontiers in Plant Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Plant Sci.</abbrev-journal-title>
<issn pub-type="epub">1664-462X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fpls.2017.01859</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Plant Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Photosynthetic Entrainment of the Circadian Clock Facilitates Plant Growth under Environmental Fluctuations: Perspectives from an Integrated Model of Phase Oscillator and Phloem Transportation</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name><surname>Ohara</surname> <given-names>Takayuki</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="author-notes" rid="fn001"><sup>&#x0002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/456114/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Satake</surname> <given-names>Akiko</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/27234/overview"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Graduate School of Environmental Science, Hokkaido University</institution>, <addr-line>Sapporo</addr-line>, <country>Japan</country></aff>
<aff id="aff2"><sup>2</sup><institution>Department of Biology, Faculty of Science, Kyushu University</institution>, <addr-line>Fukuoka</addr-line>, <country>Japan</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Katrin Kahlen, Hochschule Geisenheim University, Germany</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Norihito Nakamichi, Nagoya University, Japan; Matt Jones, University of Essex, United Kingdom</p></fn>
<fn fn-type="corresp" id="fn001"><p>&#x0002A;Correspondence: Takayuki Ohara <email>t.ohara.63920&#x00040;ees.hokudai.ac.jp</email></p></fn>
<fn fn-type="other" id="fn002"><p>This article was submitted to Plant Biophysics and Modeling, a section of the journal Frontiers in Plant Science</p></fn></author-notes>
<pub-date pub-type="epub">
<day>30</day>
<month>10</month>
<year>2017</year>
</pub-date>
<pub-date pub-type="collection">
<year>2017</year>
</pub-date>
<volume>8</volume>
<elocation-id>1859</elocation-id>
<history>
<date date-type="received">
<day>31</day>
<month>08</month>
<year>2017</year>
</date>
<date date-type="accepted">
<day>11</day>
<month>10</month>
<year>2017</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2017 Ohara and Satake.</copyright-statement>
<copyright-year>2017</copyright-year>
<copyright-holder>Ohara and Satake</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract>
<p>Plants need to avoid carbon starvation and resultant growth inhibition under fluctuating light environments to ensure optimal growth and reproduction. As diel patterns of carbon metabolism are influenced by the circadian clock, appropriate regulation of the clock is essential for plants to properly manage their carbon resources. For proper adjustment of the circadian phase, higher plants utilize environmental signals such as light or temperature and metabolic signals such as photosynthetic products; the importance of the latter as phase regulators has been recently elucidated. A mutant of <italic>Arabidopsis thaliana</italic> that is deficient in phase response to sugar has been shown, under fluctuating light conditions, to be unable to adjust starch turnover and to realize carbon homeostasis. Whereas, the effects of light entrainment on growth and survival of higher plants are well studied, the impact of phase regulation by sugar remains unknown. Here we show that endogenous sugar entrainment facilitates plant growth. We integrated two mathematical models, one describing the dynamics of carbon metabolism in <italic>A. thaliana</italic> source leaves and the other growth of sink tissues dependent on sucrose translocation from the source. The integrated model predicted that sugar-sensitive plants grow faster than sugar-insensitive plants under constant as well as changing photoperiod conditions. We found that sugar entrainment enables efficient carbon investment for growth by stabilizing sucrose supply to sink tissues. Our results highlight the importance of clock entrainment by both exogenous and endogenous signals for optimizing growth and increasing fitness.</p>
</abstract>
<kwd-group>
<kwd>starch degradation</kwd>
<kwd>sucrose</kwd>
<kwd>homeostasis</kwd>
<kwd>mathematical model</kwd>
<kwd>dynamics</kwd>
<kwd>phloem transportation</kwd>
<kwd><italic>Arabidopsis thaliana</italic></kwd>
</kwd-group>
<counts>
<fig-count count="8"/>
<table-count count="0"/>
<equation-count count="9"/>
<ref-count count="46"/>
<page-count count="12"/>
<word-count count="7187"/>
</counts>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>Introduction</title>
<p>Plants are inevitably exposed to daily and seasonal variations in light environments. To continuously grow in fluctuating environments, it is crucial for plants to stably supply carbon resources for respiration and growth. Plant growth in the day is supported by the supply of photosynthates, particularly soluble sugars such as sucrose that are transported from photosynthetic leaves (source tissues) to sink tissues (e.g., roots). Plants grow even in nighttime using carbon resources accumulated during the preceding daytime. <italic>Arabidopsis thaliana</italic>, a model plant, partitions a large fraction of assimilated carbon into insoluble starch, which is degraded at night to produce sucrose (Smith and Stitt, <xref ref-type="bibr" rid="B38">2007</xref>; Stitt and Zeeman, <xref ref-type="bibr" rid="B40">2012</xref>). Because early exhaustion of starch results in carbon starvation and ensuing growth inhibition (Graf et al., <xref ref-type="bibr" rid="B19">2010</xref>; Yazdanbakhsh et al., <xref ref-type="bibr" rid="B45">2011</xref>), careful management of starch metabolism is essential to cope with daily and seasonal fluctuations of light conditions.</p>
<p>In <italic>A. thaliana</italic>, starch amount increases during the day at an almost constant rate and decreases almost linearly at night (Caspar et al., <xref ref-type="bibr" rid="B5">1985</xref>; Gibon et al., <xref ref-type="bibr" rid="B17">2004</xref>; Smith et al., <xref ref-type="bibr" rid="B39">2004</xref>; Lu et al., <xref ref-type="bibr" rid="B27">2005</xref>). Plants in shorter photoperiods accumulate starch more rapidly during the day and degrade it more slowly at night than in longer photoperiods (Lu et al., <xref ref-type="bibr" rid="B27">2005</xref>). Plants also adjust the rate of starch degradation immediately in response to an unexpectedly early or late onset of night (Lu et al., <xref ref-type="bibr" rid="B27">2005</xref>; Graf et al., <xref ref-type="bibr" rid="B19">2010</xref>; Scialdone et al., <xref ref-type="bibr" rid="B35">2013</xref>). The circadian clock underlying the approximately 24-h cycle of biological processes is implicated in the control of starch metabolism (Graf and Smith, <xref ref-type="bibr" rid="B18">2011</xref>). Wild type <italic>A. thaliana</italic> (Ws) exhausts starch reserves about 24 h after the last dawn even under non-24 h light/dark cycles (T-cycles) (Graf et al., <xref ref-type="bibr" rid="B19">2010</xref>), indicating that the timing of starch exhaustion is programmed by the circadian clock. In the circadian clock mutant <italic>cca1/lhy</italic>, in which the functional clock has a period of about 17 h (Locke et al., <xref ref-type="bibr" rid="B26">2005</xref>), the depletion of starch occurs prematurely under a 24-h T-cycle but coincides with dawn under a 17-h T-cycle (Graf et al., <xref ref-type="bibr" rid="B19">2010</xref>). These studies suggest that coordination of the internal timing of starch turnover with environmental cycles is necessary to avoid carbon starvation.</p>
<p>Phase adjustment of the circadian clock to external stimuli such as light or temperature is fundamental for synchronizing biological processes with environments (Johnson et al., <xref ref-type="bibr" rid="B22">2003</xref>). In addition to signals from the external environment, endogenous signals such as photosynthates are also important regulators of the circadian phase in <italic>A. thaliana</italic>. We previously reported that phase adjustment by sugar is necessary for plants to flexibly regulate carbon metabolism in fluctuating light environments (Seki et al., <xref ref-type="bibr" rid="B37">2017</xref>). We developed a phase oscillator model describing phase regulation of the circadian clock by sucrose. This model predicted that phase adjustment of the circadian clock by sucrose is crucial for homeostatic regulation of carbon resources. These theoretical predictions were confirmed by physiological experiments using the mutant <italic>pseudoresponse regulator 7&#x02013;11</italic> (<italic>prr7&#x02013;11</italic>), the circadian clock of which does not show clear phase response to sucrose pulse (Haydon et al., <xref ref-type="bibr" rid="B21">2013</xref>).</p>
<p>Whereas, clock entrainment by exogenous signals such as light has been shown to be advantageous for competition and survival in several organisms (Woelfle et al., <xref ref-type="bibr" rid="B43">2004</xref>; Dodd et al., <xref ref-type="bibr" rid="B11">2005</xref>), the advantages of clock entrainment by photosynthetic products remain elusive. Here we theoretically evaluate the effects of clock entrainment by sugar on plant growth, a good proxy for plant fitness (Younginger et al., <xref ref-type="bibr" rid="B46">2017</xref>). We extended the phase oscillator model for the circadian clock (Seki et al., <xref ref-type="bibr" rid="B37">2017</xref>) by incorporating growth dynamics of sink tissues, including shoot apical meristems and roots. Growth dynamics are described by modeling phloem transportation of sucrose from source to sink tissues (Seki et al., <xref ref-type="bibr" rid="B36">2015</xref>; Satake et al., <xref ref-type="bibr" rid="B34">2016</xref>). We demonstrate that plant growth is facilitated by endogenous sugar entrainment under long photoperiods because the entrainment enables the stable supply of sucrose from source to sink tissues irrespective of light fluctuation. In short photoperiods, however, the effect of sugar entrainment on growth is negligibly small. Our results provide important theoretical evidence that circadian-phase adjustment by endogenous signals is advantageous for plant growth.</p>
</sec>
<sec id="s2">
<title>Model and methods</title>
<p>To investigate the effect of clock regulation by sugar on plant growth, we integrated two previously developed models, one describing the dynamics of starch and sucrose metabolism in source leaves (Seki et al., <xref ref-type="bibr" rid="B37">2017</xref>) and the other growth of sink tissues dependent on phloem transportation of sucrose from source tissues (Seki et al., <xref ref-type="bibr" rid="B36">2015</xref>; Satake et al., <xref ref-type="bibr" rid="B34">2016</xref>). Photosynthetic products in the source leaf are partitioned into sucrose and starch (Seki et al., <xref ref-type="bibr" rid="B37">2017</xref>; Figure <xref ref-type="fig" rid="F1">1A</xref>). Sucrose in the source leaf is loaded into the phloem, moves through the phloem tube, and is unloaded at sink tissues where it is used for respiration and growth (Seki et al., <xref ref-type="bibr" rid="B36">2015</xref>; Satake et al., <xref ref-type="bibr" rid="B34">2016</xref>; Figures <xref ref-type="fig" rid="F1">1B&#x02013;D</xref>). These two models were coupled by incorporating a term that represents sucrose translocation from source to sink. We explain the detailed structures of each component of our new integrated model in the following sections. We consider two sink tissues, the shoot and root apical meristems (SAM and RAM, respectively), because plant growth and development mainly occur in these organs (Figure <xref ref-type="fig" rid="F1">1C</xref>). Our model can be extended to the structure including multiple sinks in the complex phloem network as studied previously (Seki et al., <xref ref-type="bibr" rid="B36">2015</xref>; Satake et al., <xref ref-type="bibr" rid="B34">2016</xref>).</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p>Model describing sugar dynamics in source and sink tissues and translocation of sucrose through the phloem tube. <bold>(A)</bold> In a source leaf, the circadian clock and carbon metabolism are reciprocally regulated. <bold>(B)</bold> Translocated sucrose is used for growth and respiration in the sink. Each structural component is expressed by a cylinder of radius <italic>r</italic><sub>X</sub> and height <italic>l</italic><sub>X</sub>, where X is either G (source), T (tube), or Y (sink). <bold>(C)</bold> Tubal structure of the model. <italic>p</italic><sub><italic>i</italic></sub> represents hydrostatic pressure at the apex of component <italic>i</italic> (<italic>i</italic> &#x003F5; {0,&#x02026;, 5}). Xylem adjacent to the phloem tube is also schematized. SAM, shoot apical meristem; RAM, root apical meristem. <bold>(D)</bold> Fluxes in the model. Red arrows (<italic>J</italic><sub><italic>i</italic></sub>) indicate phloem sap flow. Light blue arrows (<italic>w</italic><sub><italic>i</italic></sub>) indicate pure water flow due to osmosis, which occurs at the region represented by light blue circles.</p></caption>
<graphic xlink:href="fpls-08-01859-g0001.tif"/>
</fig>
<sec>
<title>Sugar dynamics in source leaves and sink tissues</title>
<sec>
<title>Source leaves</title>
<p>In source leaves (Figure <xref ref-type="fig" rid="F1">1A</xref>), carbon is assimilated by photosynthesis at a rate <italic>a</italic> during the light period. The length of the light period is given by &#x003C4;<sub>L</sub>. A fraction &#x003B3; of total photoassimilates is partitioned into starch (<italic>C</italic>) for storage and a fraction 1 &#x02212; &#x003B3; is partitioned into sucrose (<italic>S</italic><sub>G</sub>). Sucrose is consumed for respiration at a rate <italic>h</italic><sub>G</sub> and for transportation at a rate &#x003B7;<sub><italic>G</italic></sub>. Starch is degraded into sucrose at a rate &#x003B2;, which is assumed to be under the control of the circadian clock and thus is a function of the phase &#x003D5; of the circadian oscillator (Seki et al., <xref ref-type="bibr" rid="B37">2017</xref>). These processes are formalized as follows:</p>
<disp-formula id="E1"><label>(1)</label><mml:math id="M1"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mfrac><mml:mrow><mml:mi>d</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac><mml:msub><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mi>G</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>a</mml:mi><mml:mi>L</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>-</mml:mo><mml:mi>&#x003B3;</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>&#x0002B;</mml:mo><mml:mi>&#x003B2;</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>&#x003D5;</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:msup><mml:mrow><mml:mi>C</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x003BA;</mml:mi></mml:mrow></mml:msup><mml:mo>-</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>h</mml:mi></mml:mrow><mml:mrow><mml:mi>G</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003B7;</mml:mi></mml:mrow><mml:mrow><mml:mi>G</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:msub><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mi>G</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E2"><label>(2)</label><mml:math id="M2"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mfrac><mml:mrow><mml:mi>d</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac><mml:mi>C</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>a</mml:mi><mml:mi>L</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mi>&#x003B3;</mml:mi><mml:mo>-</mml:mo><mml:mi>&#x003B2;</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>&#x003D5;</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:msup><mml:mrow><mml:mi>C</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x003BA;</mml:mi></mml:mrow></mml:msup><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>L</italic>(<italic>t</italic>) indicates the light condition (defined as 1 under light and 0 under dark) and &#x003BA; is a constant. Starch reserve is accumulated during the light period at a rate determined by balance between <italic>a</italic>&#x003B3; and &#x003B2;<italic>C</italic><sup>&#x003BA;</sup>. The starch degradation rate &#x003B2;(&#x003D5;) is assumed to show diel oscillation due to regulation by the circadian clock. We assume that &#x003B2;(&#x003D5;) shows a peak at dawn and a trough at the subjective dusk &#x003D5;<sup>&#x0002A;</sup> because a previous study showed that this oscillation pattern is ideal to minimize fluctuations in the sucrose supply to sinks (Figure <xref ref-type="supplementary-material" rid="SM1">S1A</xref>; Seki et al., <xref ref-type="bibr" rid="B37">2017</xref>).</p>
<p>In the previous model (Seki et al., <xref ref-type="bibr" rid="B37">2017</xref>), transportation of sucrose from source leaves was assumed to be constant. In our new model, sucrose transportation is assumed to occur based on the pressure-flow hypothesis (M&#x000FC;nch, <xref ref-type="bibr" rid="B30">1930</xref>) with an assumption that the flux in phloem obeys the Hagen&#x02013;Poiseuille law (Supplementary Material, section 1). Sucrose dynamics are now described both in the source leaves and phloem tubes. At the source leaves, sucrose is loaded into the adjacent phloem tube at a rate &#x003B7;<sub>G</sub> (Figure <xref ref-type="fig" rid="F1">1A</xref>). At the phloem tube, loaded sucrose is transported to sink tissues (see Supplementary Material, section 1 for detailed explanation). The dynamics of sucrose concentration at the phloem tube adjacent to the source (<italic>g</italic><sub>0</sub>) is given by</p>
<disp-formula id="E9"><label>(3)</label><mml:math id="M14"><mml:mfrac><mml:mi>d</mml:mi><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac><mml:msub><mml:mi>g</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow></mml:mfrac><mml:mo>&#x0007B;</mml:mo><mml:msub><mml:mi>&#x003B7;</mml:mi><mml:mi>G</mml:mi></mml:msub><mml:msub><mml:mi>S</mml:mi><mml:mi>G</mml:mi></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:msub><mml:mi>g</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:msub><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mo>&#x02212;</mml:mo><mml:msub><mml:mi>J</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>+</mml:mo></mml:msub><mml:mo>&#x02212;</mml:mo><mml:msub><mml:mi>g</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:msub><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:msub><mml:mi>J</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>+</mml:mo></mml:msub><mml:mo>&#x0007D;</mml:mo><mml:mo>,</mml:mo></mml:math></disp-formula>
<p>where <italic>g</italic><sub>1</sub>(<italic>t</italic>)[&#x02013;<italic>J</italic><sub>1</sub>(<italic>t</italic>)]<sub>&#x0002B;</sub> &#x02212; <italic>g</italic><sub>0</sub>(<italic>t</italic>)[<italic>J</italic><sub>1</sub>(<italic>t</italic>)]<sub>&#x0002B;</sub> describes the rate of sucrose change due to flux <italic>J</italic><sub>1</sub>(<italic>t</italic>) at phloem tube 1 (Figures <xref ref-type="fig" rid="F1">1C,D</xref>). <italic>g</italic><sub>1</sub> and <italic>V</italic><sub>0</sub> represent sucrose concentration at tube 1 and the volume of tube 0, respectively (Figure <xref ref-type="fig" rid="F1">1C</xref>).</p>
<p>Similar to the previous model (Seki et al., <xref ref-type="bibr" rid="B37">2017</xref>), the phase &#x003D5; of the circadian oscillator is modeled by</p>
<disp-formula id="E3"><label>(4)</label><mml:math id="M3"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mfrac><mml:mrow><mml:mi>d</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac><mml:mi>&#x003D5;</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>&#x003C9;</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mi>L</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>&#x003D5;</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:msub><mml:mrow><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:mi>L</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mo>&#x0007E;</mml:mo></mml:mover></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mi>S</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>&#x003D5;</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:msub><mml:mrow><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:mi>S</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mo>&#x0007E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>G</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where &#x003C9; is the angular frequency of the oscillator. The term <italic>Z</italic><sub>L</sub>(&#x003D5;)<italic>f</italic><sub>L</sub>(<inline-formula><mml:math id="M4"><mml:mover accent="true"><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mo>&#x0007E;</mml:mo></mml:mover></mml:math></inline-formula>) represents the effect of light stimulus, which is assumed to reset the phase to 0 at dawn and to &#x003C4;<sub>L</sub> at dusk (Seki et al., <xref ref-type="bibr" rid="B37">2017</xref>). <italic>Z</italic><sub>S</sub>(&#x003D5;) is a phase response curve (PRC) to a sucrose pulse showing phase advance until &#x003D5; &#x0003D; &#x003D5;<sup>&#x0002A;</sup> and phase delay thereafter as previously determined (Figure <xref ref-type="supplementary-material" rid="SM1">S1B</xref>; Seki et al., <xref ref-type="bibr" rid="B37">2017</xref>). Sugar input <italic>f</italic><sub>S</sub> (Supplementary Material, section 2) is defined by the Hill function of the rate of change in sucrose level (i.e., <inline-formula><mml:math id="M5"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mo>&#x0007E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>G</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mi>d</mml:mi><mml:msub><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mi>G</mml:mi></mml:mrow></mml:msub><mml:mo>/</mml:mo><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:math></inline-formula>). The functions of <italic>Z</italic><sub>S</sub>(&#x003D5;) and <italic>f</italic><sub>S</sub> have been demonstrated to be optimal for minimizing sucrose fluctuation (Seki et al., <xref ref-type="bibr" rid="B37">2017</xref>). We assume that the phase shift of the circadian clock by sugar takes place only in the light period because the phase shift at night did not improve sucrose homeostasis (Seki et al., <xref ref-type="bibr" rid="B37">2017</xref>).</p>
</sec>
<sec>
<title>Sink tissues</title>
<p>Translocated sucrose is unloaded into the sinks at a rate &#x003B7;<sub>Y</sub> from the adjacent phloem tubes with sucrose concentration <italic>g</italic><sub><italic>i</italic></sub> (Figure <xref ref-type="fig" rid="F1">1B</xref>). Sucrose in the sink (<italic>S</italic><sub>Y</sub>) is consumed for respiration and growth at rates <italic>h</italic><sub>Y</sub> and &#x003B1;(<italic>S</italic><sub>Y</sub>), respectively. These processes are formalized as follows:</p>
<disp-formula id="E4"><label>(5)</label><mml:math id="M6"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mfrac><mml:mrow><mml:mi>d</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac><mml:msub><mml:mrow><mml:mi>g</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>g</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mo>-</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>J</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mo>&#x0002B;</mml:mo></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>g</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>J</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mo>&#x0002B;</mml:mo></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003B7;</mml:mi></mml:mrow><mml:mrow><mml:mtext>Y</mml:mtext></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mo>}</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E5"><label>(6)</label><mml:math id="M7"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mfrac><mml:mrow><mml:mi>d</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac><mml:msub><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mtext>Y</mml:mtext></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003B7;</mml:mi></mml:mrow><mml:mrow><mml:mtext>Y</mml:mtext></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mi>g</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>h</mml:mi></mml:mrow><mml:mrow><mml:mtext>Y</mml:mtext></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mtext>Y</mml:mtext></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>-</mml:mo><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mtext>Y</mml:mtext></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>V</italic><sub><italic>i</italic></sub> is the volume of component <italic>i</italic> (<italic>i</italic> &#x0003D; 4 for SAM and <italic>i</italic> &#x0003D; 5 for RAM; <italic>i</italic> &#x0003D; 1, 2, and 3 for connecting tubes between the source and sink; Figure <xref ref-type="fig" rid="F1">1C</xref>). The first and second terms in the right-hand side of Equation (5) represent the solution inflow and outflow, respectively (see Supplementary Material, section 1). The function for sucrose consumption rate for growth &#x003B1;(<italic>S</italic><sub>Y</sub>) will be explained in the next subsection. To clarify the effect of sugar entrainment on plant growth, we simplified the model by assuming that the sugar and growth dynamics of SAM and RAM are identical. Therefore, the values of &#x003B7;<sub>Y</sub> and <italic>h</italic><sub>Y</sub> as well as the parameters in &#x003B1;(<italic>S</italic><sub>Y</sub>) are the same in the two sinks.</p>
</sec>
</sec>
<sec>
<title>Growth of sink tissues</title>
<sec>
<title>Growth dynamics</title>
<p>We formalized the growth kinetics of the sink tissue based on sucrose supply because a strong correlation between growth rate and sucrose supply has been reported in both the light and dark periods for <italic>A. thaliana</italic> (Sulpice et al., <xref ref-type="bibr" rid="B41">2014</xref>; Mengin et al., <xref ref-type="bibr" rid="B28">2017</xref>). When sugar supply is sufficient, growth is promoted by the target of rapamycin (TOR) kinase, the expression level of which correlates with <italic>A. thaliana</italic> shoot and root growth (Deprost et al., <xref ref-type="bibr" rid="B10">2007</xref>; Lastdrager et al., <xref ref-type="bibr" rid="B24">2014</xref>). On the contrary, Snf1-related kinase 1 (SnRK1) inhibits growth in response to low carbon availability (Baena-Gonz&#x000E1;lez and Sheen, <xref ref-type="bibr" rid="B2">2008</xref>; Lastdrager et al., <xref ref-type="bibr" rid="B24">2014</xref>). In addition, growth rate is likely to be saturated as sucrose supply is increased (Sulpice et al., <xref ref-type="bibr" rid="B41">2014</xref>). Given these empirical findings, we assume that growth-related sucrose consumption rate is an increasing and saturating function of sucrose supply (Figure <xref ref-type="fig" rid="F2">2A</xref>):</p>
<disp-formula id="E6"><label>(7)</label><mml:math id="M8"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mtext>Y</mml:mtext></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003B1;</mml:mi></mml:mrow><mml:mrow><mml:mtext>Max</mml:mtext></mml:mrow></mml:msub><mml:mfrac><mml:mrow><mml:msubsup><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mtext>Y</mml:mtext></mml:mrow><mml:mrow><mml:mi>&#x000F1;</mml:mi></mml:mrow></mml:msubsup></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>K</mml:mi></mml:mrow><mml:mo>&#x0007E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>&#x000F1;</mml:mi></mml:mrow></mml:msup><mml:mo>&#x0002B;</mml:mo><mml:msubsup><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mtext>Y</mml:mtext></mml:mrow><mml:mrow><mml:mi>&#x000F1;</mml:mi></mml:mrow></mml:msubsup></mml:mrow></mml:mfrac><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where &#x003B1;<sub>Max</sub>, &#x000F1;, and <inline-formula><mml:math id="M9"><mml:mover accent="true"><mml:mrow><mml:mi>K</mml:mi></mml:mrow><mml:mo>&#x0007E;</mml:mo></mml:mover></mml:math></inline-formula> are constants. These values are estimated using the published data of fresh biomass in <italic>A. thaliana</italic> (Caspar et al., <xref ref-type="bibr" rid="B5">1985</xref>) as explained later. The rate of increase in sink fresh biomass (<italic>W</italic><sub>Y</sub>) is then described by</p>
<disp-formula id="E7"><label>(8)</label><mml:math id="M10"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mfrac><mml:mrow><mml:mi>d</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac><mml:msub><mml:mrow><mml:mi>W</mml:mi></mml:mrow><mml:mrow><mml:mtext>Y</mml:mtext></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>&#x003BB;</mml:mi><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mtext>Y</mml:mtext></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:msub><mml:mrow><mml:mi>W</mml:mi></mml:mrow><mml:mrow><mml:mtext>Y</mml:mtext></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where &#x003BB; is the conversion rate of sucrose for growth. Therefore, growth rate is represented by &#x003BB;&#x003B1;(<italic>S</italic><sub>Y</sub>). Since dry weight and fresh weight of <italic>A. thaliana</italic> Col-0 display qualitatively similar increase patterns (Caspar et al., <xref ref-type="bibr" rid="B5">1985</xref>; Christophe et al., <xref ref-type="bibr" rid="B8">2008</xref>), our model can also be applicable to the analysis of dry biomass by appropriate scaling. We tested the robustness of our results by using an alternative formalization of growth rate as a linear increasing function of sucrose (Supplementary Material, section 3; Figures <xref ref-type="supplementary-material" rid="SM1">S4</xref>, <xref ref-type="supplementary-material" rid="SM1">S5</xref>).</p>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p><bold>(A)</bold> Sucrose consumption rate &#x003B1; for growth of the sink (SAM or RAM). The unit for sucrose is &#x003BC;molC<sub>6</sub> g<sup>&#x02212;1</sup>FW. <bold>(B)</bold> Time evolution of growth of the wild type (black) and the starchless <italic>pgm</italic> mutant (red) grown in a 12-h photoperiod. In <bold>(B)</bold>, lines represent the sum of the simulated growth of the two sinks and circles represent the published growth data (Caspar et al., <xref ref-type="bibr" rid="B5">1985</xref>).</p></caption>
<graphic xlink:href="fpls-08-01859-g0002.tif"/>
</fig>
</sec>
<sec>
<title>Parameter estimation of the sucrose consumption rate for growth</title>
<p>The parameter values of the sucrose consumption rate &#x003B1;(<italic>S</italic><sub>Y</sub>) in Equation (7) were estimated by fitting the simulated growth curves to published data measuring fresh weight of <italic>A. thaliana</italic> wild type (Col-0) and the starchless <italic>phosphoglucomutase</italic> (<italic>pgm</italic>) mutant grown in a 12-h photoperiod (Caspar et al., <xref ref-type="bibr" rid="B5">1985</xref>; Figure <xref ref-type="fig" rid="F2">2B</xref>). Because the number of the data is limited (data at four time points per genotype), we combine the data of both genotypes for the parameter estimation. The subjective dusk &#x003D5;<sup>&#x0002A;</sup> was set to 12 h because both the wild type and <italic>pgm</italic> were likely to be completely entrained to a 12-h photoperiod in the experiment. We simulated the growth of the wild type (see next subsection) and starchless mutant using an initial value of 0.0005 at time <italic>t</italic> &#x0003D; 0 (corresponding to the first observation of plant fresh weight in the experiment). Because the <italic>pgm</italic> mutant does not accumulate starch, the carbon partitioning rate for starch (&#x003B3; in Equations 1 and 2) was set to 0 for the mutant. The parameter values &#x003B1;<sub>Max</sub>, &#x000F1;, and <inline-formula><mml:math id="M11"><mml:mover accent="true"><mml:mrow><mml:mi>K</mml:mi></mml:mrow><mml:mo>&#x0007E;</mml:mo></mml:mover></mml:math></inline-formula> in Equation (7) were estimated by minimizing a following cost function:</p>
<disp-formula id="E8"><label>(9)</label><mml:math id="M12"><mml:mtable columnalign='left'><mml:mtr><mml:mtd><mml:mi>P</mml:mi><mml:mo stretchy='false'>(</mml:mo><mml:msub><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mtext>Max</mml:mtext></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mi>&#x000F1;</mml:mi><mml:mo>,</mml:mo><mml:mover accent='true'><mml:mi>K</mml:mi><mml:mo>&#x002DC;</mml:mo></mml:mover><mml:mo stretchy='false'>)</mml:mo><mml:mo>=</mml:mo><mml:mstyle displaystyle='true'><mml:munder><mml:mo>&#x02211;</mml:mo><mml:mi>X</mml:mi></mml:munder><mml:mrow><mml:mstyle displaystyle='true'><mml:munderover><mml:mo>&#x02211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mn>4</mml:mn></mml:munderover><mml:mrow><mml:msup><mml:mrow><mml:mo>&#x0007B;</mml:mo><mml:msubsup><mml:mi>W</mml:mi><mml:mi>i</mml:mi><mml:mrow><mml:mi>E</mml:mi><mml:mi>X</mml:mi><mml:mi>P</mml:mi></mml:mrow></mml:msubsup><mml:mo>&#x02212;</mml:mo><mml:mn>2</mml:mn><mml:msub><mml:mi>W</mml:mi><mml:mrow><mml:mi>Y</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msub><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mtext>Max</mml:mtext></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mi>&#x000F1;</mml:mi><mml:mo>,</mml:mo><mml:mover accent='true'><mml:mi>K</mml:mi><mml:mo>&#x002DC;</mml:mo></mml:mover></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x0007D;</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mstyle></mml:mrow></mml:mstyle><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <inline-formula><mml:math id="M13"><mml:msubsup><mml:mrow><mml:mi>W</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>E</mml:mi><mml:mi>X</mml:mi><mml:mi>P</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> and <italic>W</italic><sub>Y,<italic>i</italic></sub> are fresh weight at the <italic>i</italic>th time point in the experiment and the growth at corresponding times in the simulation, respectively (Figure <xref ref-type="fig" rid="F2">2B</xref>). <italic>X</italic> corresponds to the wild type and <italic>pgm</italic> mutant. As fresh weight is likely to be measured from whole plants in the experiment (Caspar et al., <xref ref-type="bibr" rid="B5">1985</xref>), we used the sum (2<italic>W</italic><sub>Y,<italic>i</italic></sub>) of the growth of SAM (<italic>W</italic><sub>Y,<italic>i</italic></sub>) and RAM (also <italic>W</italic><sub>Y,<italic>i</italic></sub>) for the parameter estimation.</p>
</sec>
<sec>
<title>Simulation conditions</title>
<p>To evaluate the effect of clock entrainment by sugar on growth, we simulated the growth dynamics of sugar-sensitive (wild type) and sugar-insensitive (mutant) plants in constant photoperiods (ranging from 8 to 16 h) as well as under changing photoperiod conditions. The wild type adjusts the phase of the circadian clock by sugar as formalized in Equation (4), while the sugar-insensitive mutant lacks this response to sugar [i.e., <italic>f</italic><sub>S</sub>(<inline-formula><mml:math id="M15"><mml:msub><mml:mover accent='true'><mml:mi>S</mml:mi><mml:mo>&#x002DC;</mml:mo></mml:mover><mml:mtext>G</mml:mtext></mml:msub></mml:math></inline-formula>) in Equation (4) is zero]. Both plants respond to light signals in the same manner. Although phase regulation by sugar reduces fluctuation of carbon resources, sucrose dynamics still deviate from homeostasis in both the wild type and mutant unless these plants precisely predict the timing of dusk (i.e., &#x003D5;<sup>&#x0002A;</sup> &#x0003D; &#x003C4;<sub>L</sub>; Seki et al., <xref ref-type="bibr" rid="B37">2017</xref>). To investigate the effect of sucrose homeostasis on growth, we also consider an ideal plant that can maintain perfect homeostasis in a steady state in any photoperiod (hereafter termed the homeostatic plant). We previously determined the optimal function for starch degradation rate to minimize sucrose fluctuation in a given photoperiod (Seki et al., <xref ref-type="bibr" rid="B37">2017</xref>). The homeostatic plant is assumed to possess these functions in each photoperiod and is simulated by setting &#x003D5;<sup>&#x0002A;</sup> being equal to &#x003C4;<sub>L</sub> in any condition (&#x003C4;<sub>L</sub> &#x003F5;(0, 24]) in contrast to the fixed value of &#x003D5;<sup>&#x0002A;</sup> in the wild type and sugar-insensitive mutant (summarized in a table in Figure <xref ref-type="supplementary-material" rid="SM1">S1</xref>). Comparison of growth dynamics among these three plants enables us to address the differential effects of endogenous sugar entrainment and sucrose homeostasis on growth. Because SAM and RAM grow at similar rates, we present the growth dynamics of SAM as the sink growth.</p>
<p>Under the constant photoperiod conditions, plant growth was simulated over 10 days at various photoperiods (8, 9, 10, &#x02026;, 16 h). We then compared the growth increment over 10 days among plant types by calculating the differences for the following pairs: (wild type&#x02014;mutant), (homeostatic&#x02014;mutant), and (homeostatic&#x02014;wild type). Under the changing photoperiod conditions, plants grown for 5 days in an 8-h or 16-h photoperiod were transferred to a 16-h or 8-h photoperiod and grown for an additional 5 days (Seki et al., <xref ref-type="bibr" rid="B37">2017</xref>). For the wild type and sugar-insensitive mutant, we mainly analyzed the plants with the subjective dusk &#x003D5;<sup>&#x0002A;</sup> of 10 h. In addition, we also simulated 10-day growth of the plants with various values of &#x003D5;<sup>&#x0002A;</sup> (8, 9, 11, 12 h) in constant photoperiods to investigate the influence of &#x003D5;<sup>&#x0002A;</sup> on growth dynamics. Other parameter values are listed in Table <xref ref-type="supplementary-material" rid="SM1">S1</xref>. Numerical integration of the ordinary differential equations was performed with the fourth-order Runge-Kutta method using Mathematica (version 10; Wolfram Research).</p>
</sec>
</sec>
</sec>
<sec sec-type="results" id="s3">
<title>Results</title>
<p>When the plants are grown in a constant photoperiod, the growth of the sink (SAM or RAM) of the all three types increases as the photoperiod is lengthened (Figure <xref ref-type="fig" rid="F3">3A</xref>). The sugar-sensitive wild type and homeostatic plant grow significantly faster than the sugar-insensitive mutant under long days (Figure <xref ref-type="fig" rid="F3">3B</xref>). When photoperiod is shorter than 12 h, the growth difference between the mutant and the others almost disappears. Growth of the wild type is slower than that of the homeostatic plant in long photoperiods (Figure <xref ref-type="fig" rid="F3">3B</xref>), suggesting that the minimization of sucrose fluctuation is the most effective strategy for efficient growth under these conditions. We confirmed the robustness of our results under different parameter values for carbon metabolism (Figure <xref ref-type="supplementary-material" rid="SM1">S2</xref>) as well as for the phloem tube network (Figure <xref ref-type="supplementary-material" rid="SM1">S3</xref>).</p>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption><p><bold>(A)</bold> Growth of the sink (SAM or RAM) of the mutant, wild type, and homeostatic plant over 10 days under constant photoperiod conditions. <bold>(B)</bold> Growth difference between the wild type and mutant (upper), between the homeostatic and mutant (middle), and between the homeostatic and wild type (lower).</p></caption>
<graphic xlink:href="fpls-08-01859-g0003.tif"/>
</fig>
<p>Under a long day (16 L/ 8 D), sucrose levels in the source and sink (SAM or RAM) are highly variable in the sugar-insensitive mutant, moderately variable in the wild type, and almost constant in the homeostatic plant (Figures <xref ref-type="fig" rid="F4">4A,B</xref>). Diel patterns of growth reflect sucrose profiles in the sink, revealing a substantially lower growth rate during the evening in the mutant and a moderately lower growth rate during the evening in the wild type compared to the homeostatic plant (Figure <xref ref-type="fig" rid="F4">4C</xref>). The difference in evening growth rate is the major reason for differential growth among the three plant types in a long photoperiod. The mutant accumulates the largest amount of starch, with some amount unused even at the end of night (Figure <xref ref-type="fig" rid="F4">4D</xref>), indicating inefficient translocation of carbon. In contrast, the wild type accumulates less starch than the mutant and invests a larger amount of photoassimilate for growth (Figure <xref ref-type="fig" rid="F4">4D</xref>).</p>
<fig id="F4" position="float">
<label>Figure 4</label>
<caption><p>Predicted profiles of <bold>(A)</bold> sucrose in the source, <bold>(B)</bold> sucrose in the sink (SAM or RAM), <bold>(C)</bold> growth rate of the sink (SAM or RAM), <bold>(D)</bold> starch, <bold>(E)</bold> phase, and <bold>(F)</bold> starch degradation rate of the plants grown in a 16-h photoperiod. The unit for sucrose and starch is &#x003BC;molC<sub>6</sub> g<sup>&#x02212;1</sup>FW. White background, light period; Gray background, dark period.</p></caption>
<graphic xlink:href="fpls-08-01859-g0004.tif"/>
</fig>
<p>The different sucrose, starch, and growth profiles between the wild type and sugar-insensitive mutant are caused by the clock plasticity in response to sugar. The decrease in sucrose levels at dawn (Figure <xref ref-type="fig" rid="F4">4A</xref>) is sensed by the wild type as a negative sugar signal, driving a phase delay of the circadian oscillator (Figure <xref ref-type="fig" rid="F4">4E</xref>). This phase delay in the morning increases the starch degradation rate during the light period (Figure <xref ref-type="fig" rid="F4">4F</xref>), resulting in elevation of the sucrose level (Figure <xref ref-type="fig" rid="F4">4A</xref>) and decrease of the starch level (Figure <xref ref-type="fig" rid="F4">4D</xref>). Such a phase shift in the circadian oscillator does not occur in the mutant (Figure <xref ref-type="fig" rid="F4">4E</xref>). Excessively high sucrose in the mutant around dawn (Figure <xref ref-type="fig" rid="F4">4B</xref>) does not significantly contribute to growth (Figure <xref ref-type="fig" rid="F4">4C</xref>) due to the saturating property of the growth-related sucrose consumption rate &#x003B1; (Figure <xref ref-type="fig" rid="F2">2A</xref>), while the sucrose decrease around dusk substantially reduces growth rate since at this low level &#x003B1; is almost linearly dependent on sucrose.</p>
<p>Under a short day (8 L/16 D; Figure <xref ref-type="fig" rid="F5">5</xref>), the wild type and sugar-insensitive mutant show the similar sugar and growth dynamics as reported previously (Seki et al., <xref ref-type="bibr" rid="B37">2017</xref>). The elevation of sucrose at dawn (Figure <xref ref-type="fig" rid="F5">5A</xref>) gives rise to a positive sugar signal and resultant phase advance in the wild type (Figure <xref ref-type="fig" rid="F5">5E</xref>). However, this phase shift does not substantially change the starch degradation rate (Figure <xref ref-type="fig" rid="F5">5F</xref>), so the sugar and growth dynamics are almost the same for the wild type and mutant (Figures <xref ref-type="fig" rid="F5">5A&#x02013;D</xref>). The homeostatic plant displays distinct growth dynamics (Figure <xref ref-type="fig" rid="F5">5C</xref>). Nevertheless, the difference in 10-day growth among plants is very small in the short photoperiod (Figure <xref ref-type="fig" rid="F3">3B</xref>). Since all three plants fix less carbon per day compared to under long day conditions, sucrose concentration is low both in the source and sink (Figures <xref ref-type="fig" rid="F5">5A,B</xref>). This leads to a slow growth rate due to low supply of sucrose. Because all three plants are limited by sucrose under short day conditions, the growth difference among them is indistinguishably small.</p>
<fig id="F5" position="float">
<label>Figure 5</label>
<caption><p>Predicted profiles of <bold>(A)</bold> sucrose in the source, <bold>(B)</bold> sucrose in the sink (SAM or RAM), <bold>(C)</bold> growth rate of the sink (SAM or RAM), <bold>(D)</bold> starch, <bold>(E)</bold> phase, and <bold>(F)</bold> starch degradation rate of the plants grown in an 8-h photoperiod. Details are as in the legend to Figure <xref ref-type="fig" rid="F4">4</xref>.</p></caption>
<graphic xlink:href="fpls-08-01859-g0005.tif"/>
</fig>
<p>We examined the effect of the timing of the subjective dusk &#x003D5;<sup>&#x0002A;</sup> on growth dynamics of the wild type (Figures <xref ref-type="fig" rid="F6">6</xref>, <xref ref-type="fig" rid="F7">7</xref>). Changes in &#x003D5;<sup>&#x0002A;</sup> do not markedly alter the values of 10-day growth of the sink and its difference among plant types (Figure <xref ref-type="fig" rid="F6">6</xref>). On the other hand, growth rate is strongly dependent on the timing of both the internal and external dusk. When the external dusk occurs later than the internal dusk (i.e., the length of the light period &#x003C4;<sub>L</sub> is larger than &#x003D5;<sup>&#x0002A;</sup>), the growth rate peaks around dawn and decreases around dusk (Figures <xref ref-type="fig" rid="F4">4C</xref>, <xref ref-type="fig" rid="F7">7</xref>). As the photoperiod is shortened, amplitude of growth rate decreases and eventually becomes almost zero. When the value of &#x003C4;<sub>L</sub> is smaller than &#x003D5;<sup>&#x0002A;</sup>, the opposite growth pattern is observed, with a peak around dusk and a trough around dawn (Figures <xref ref-type="fig" rid="F5">5C</xref>, <xref ref-type="fig" rid="F7">7</xref>).</p>
<fig id="F6" position="float">
<label>Figure 6</label>
<caption><p><bold>(A)</bold> 10-day growth of the sink (SAM or RAM) of the mutant, wild type, and homeostatic plant, and <bold>(B)</bold> growth differences among the plants at various values of the subjective dusk &#x003D5;<sup>&#x0002A;</sup>.</p></caption>
<graphic xlink:href="fpls-08-01859-g0006.tif"/>
</fig>
<fig id="F7" position="float">
<label>Figure 7</label>
<caption><p>Comparison of growth rate of the sink (SAM or RAM) of the wild type at various values of the subjective dusk &#x003D5;<sup>&#x0002A;</sup>. In each panel, growth dynamics in various photoperiods is represented. Different colors correspond to different light conditions (upper left panel).</p></caption>
<graphic xlink:href="fpls-08-01859-g0007.tif"/>
</fig>
<p>When the plants are transferred from long to short or short to long days, the growth of the wild type and homeostatic plant are consistently higher than that of the sugar-insensitive mutant (Figure <xref ref-type="fig" rid="F8">8A</xref>). All three plants are able to restore normal sugar and growth dynamics in about 2 days after the transfer, although the dynamics are complex immediately after the photoperiod change (Figure <xref ref-type="fig" rid="F8">8B</xref>). Therefore, the growth differences under these conditions may reflect the results under the constant photoperiod conditions that growth of the mutant is inferior to the others especially under a long day (Figure <xref ref-type="fig" rid="F3">3</xref>). Note that the dynamics of other variables such as sucrose in source were reported in a previous study (Seki et al., <xref ref-type="bibr" rid="B37">2017</xref>).</p>
<fig id="F8" position="float">
<label>Figure 8</label>
<caption><p><bold>(A)</bold> Growth of the sink (SAM or RAM) of the mutant, wild type, and homeostatic plant transferred from an 8-h to 16-h photoperiod (left) and from a 16-h to 8-h photoperiod (right) and the growth difference among the plants. Color codes as in Figure <xref ref-type="fig" rid="F3">3A</xref>. <bold>(B)</bold> Predicted profiles of sucrose in the sink and growth rate around the photoperiod change. Details are as in the legend to Figure <xref ref-type="fig" rid="F4">4</xref>.</p></caption>
<graphic xlink:href="fpls-08-01859-g0008.tif"/>
</fig>
</sec>
<sec sec-type="discussion" id="s4">
<title>Discussion</title>
<p>Our results provide the first theoretical evidence that clock entrainment by photosynthetic products improves the fitness of higher plants. The sugar-sensitive wild type is predicted to grow faster than the sugar-insensitive mutant under long day conditions (Figure <xref ref-type="fig" rid="F3">3</xref>) because the phase shift of the circadian oscillator by sugar signals enables efficient sugar allocation for growth, while the mutant accumulates carbon as insoluble starch unusable for growth. The growth of the wild type is also higher than the mutant under changing photoperiod conditions (Figure <xref ref-type="fig" rid="F8">8A</xref>). Since the wild type and mutant are assumed to possess the same entrainment property to light, the lower growth in the mutant stems solely from the lack of sugar-induced phase adjustment. These results suggest that clock entrainment by endogenous sugar signals, in addition to entrainment by exogenous light signals, optimizes plant growth in nature, where the day length gradually changes.</p>
<p>Our results also provide important information about the internal timing of plants. Our model predicts that the growth pattern of the wild type displays a maximum around dawn and a minimum around dusk in a 16-h photoperiod (Figure <xref ref-type="fig" rid="F4">4C</xref>) or more generally under conditions where the value of &#x003C4;<sub>L</sub> (length of the light period) is larger than &#x003D5;<sup>&#x0002A;</sup> (timing of the subjective dusk) (Figure <xref ref-type="fig" rid="F7">7</xref>). When the value of &#x003C4;<sub>L</sub> is smaller than &#x003D5;<sup>&#x0002A;</sup>, the pattern is reversed (Figure <xref ref-type="fig" rid="F5">5C</xref> where &#x003C4;<sub>L</sub> &#x0003D; 8 h and &#x003D5;<sup>&#x0002A;</sup> &#x0003D; 10 h; Figure <xref ref-type="fig" rid="F7">7</xref>). Yazdanbakhsh et al. (<xref ref-type="bibr" rid="B45">2011</xref>) reported, in both 16-h and 8-h photoperiods, a similar diel pattern in root elongation growth of <italic>A. thaliana</italic> (Col-0), which possibly correlates with fresh biomass change since both growth measures consider dry matter production and water content. These findings imply that &#x003D5;<sup>&#x0002A;</sup> of this accession is about 8 h. This idea is consistent with the PRC of the same accession to a sucrose pulse showing a transition from phase advance to delay, which corresponds to the phase &#x003D5;<sup>&#x0002A;</sup>, at 6&#x02013;10 h zeitgeber time (Haydon et al., <xref ref-type="bibr" rid="B21">2013</xref>; Seki et al., <xref ref-type="bibr" rid="B37">2017</xref>). These data suggest that a relatively short light period is subjectively anticipated by this accession. Other ecotypes of <italic>A. thaliana</italic> that anticipate later subjective dusk (i.e., large &#x003D5;<sup>&#x0002A;</sup>) will display the inversed growth pattern in a short photoperiod and later phase advance-to-delay transition of the sucrose-pulsed PRC.</p>
<p>Numerous mathematical models on growth of various plants have been developed (e.g., Thornley and Johnson, <xref ref-type="bibr" rid="B42">1990</xref>; Chew et al., <xref ref-type="bibr" rid="B7">2014</xref>, <xref ref-type="bibr" rid="B6">2017</xref>; Feller et al., <xref ref-type="bibr" rid="B13">2015</xref>; Barillot et al., <xref ref-type="bibr" rid="B3">2016a</xref>,<xref ref-type="bibr" rid="B4">b</xref>). Among them, two multiscale models of <italic>A. thaliana</italic> include the circadian clock sub-model and are able to quantitatively predict growth (Chew et al., <xref ref-type="bibr" rid="B7">2014</xref>, <xref ref-type="bibr" rid="B6">2017</xref>). Because the phase of the clock is dynamically regulated by sugar signals in our model and such phase adjustment is not considered in the previous models (Chew et al., <xref ref-type="bibr" rid="B7">2014</xref>, <xref ref-type="bibr" rid="B6">2017</xref>), only our model enables us to evaluate the impact of phase regulation by sugar on plant growth. Our model correctly predicts the empirical finding that wild type <italic>A. thaliana</italic> (Col-0) and the mutant <italic>prr7&#x02013;11</italic> in which phase response to sugar is abolished grow similarly in a 12-h photoperiod (Chew et al., <xref ref-type="bibr" rid="B6">2017</xref>; Figure <xref ref-type="fig" rid="F3">3</xref>). We suggest that this is due to the relatively weak impact of sugar entrainment on growth, as our modeling results show that the sugar-induced phase shift has less effect on growth dynamics in shorter photoperiods (Figure <xref ref-type="fig" rid="F5">5</xref>). These considerations are reminiscent of the empirical finding that response of <italic>A. thaliana</italic> to induced increases in trehalose-6-phosphate, a potential signal metabolite of sucrose status, was more intense in a 16-h photoperiod than in a 12-h photoperiod (Figueroa et al., <xref ref-type="bibr" rid="B15">2016</xref>). Our model also predicts significantly faster growth of the wild type than the sugar-insensitive mutant under long day conditions (Figure <xref ref-type="fig" rid="F3">3</xref>), which will be experimentally tested using wild type <italic>A. thaliana</italic> and <italic>prr7&#x02013;11</italic>. Wild type is likely to grow faster than <italic>prr7&#x02013;11</italic> under long day conditions because the latter accumulates more starch than the former in a 16-h photoperiod (Seki et al., <xref ref-type="bibr" rid="B37">2017</xref>) which implies that <italic>prr7&#x02013;11</italic> can utilize less sucrose for growth than wild type. However, it should be noted that mutation of <italic>PRR7</italic> affects not only sugar signaling but also light signaling in the plant circadian system (Kaczorowski and Quail, <xref ref-type="bibr" rid="B23">2003</xref>; Farr&#x000E9; et al., <xref ref-type="bibr" rid="B12">2005</xref>), and thus the growth characteristics of <italic>prr7&#x02013;11</italic> can be influenced by variations in light as well as sugar entrainment. Nevertheless, sugar entrainment is likely to be a dominant factor because our previous study has theoretically shown that the deficient in sugar signaling only is sufficient to correctly predict patterns of starch turnover of <italic>prr7&#x02013;11</italic> (Seki et al., <xref ref-type="bibr" rid="B37">2017</xref>).</p>
<p>There are several possibilities to expand our model. We used a constant rate <italic>a</italic> for carbon capture, which was determined in a previous study based on experiments where nutrient levels were controlled (Gibon et al., <xref ref-type="bibr" rid="B17">2004</xref>; Feugier and Satake, <xref ref-type="bibr" rid="B14">2013</xref>). In higher plants, however, it has been demonstrated that sugar accumulation downregulates photosynthesis, possibly through decreased activity of Rubisco (Araya et al., <xref ref-type="bibr" rid="B1">2006</xref>; Ribeiro et al., <xref ref-type="bibr" rid="B32">2012</xref>; Quentin et al., <xref ref-type="bibr" rid="B31">2013</xref>; Lobo et al., <xref ref-type="bibr" rid="B25">2015</xref>). Formalizing the carbon capture rate as a decreasing function of sucrose concentration could reduce the growth of the sugar-insensitive mutant under long day conditions due to the very high sucrose level around dawn (Figures <xref ref-type="fig" rid="F4">4A,B</xref>). Regarding circadian clock properties, we considered a circadian oscillator with a period of 24 h (i.e., &#x003C9; &#x0003D; 1). A recent study has reported a positive correlation between the free-running period of natural populations of <italic>Mimulus guttatus</italic> and the latitude of their geographic origin (Greenham et al., <xref ref-type="bibr" rid="B20">2017</xref>). In <italic>A. thaliana</italic>, a similar but weak correlation has been found (Michael et al., <xref ref-type="bibr" rid="B29">2003</xref>), although advantages of such variation in the period remain elusive. Growth simulations of plants with fast- and slow-running clocks could potentially reveal the adaptive value of natural variation in the circadian period in terms of growth optimization. Moreover, we can consider the possibility that two or more clock genes are involved in the phase response to sugar, where each gene can have differential responsiveness. Although the crucial role of <italic>PRR7</italic> in the response has been established (Haydon et al., <xref ref-type="bibr" rid="B21">2013</xref>), other clock genes could participate in it; for instance, <italic>CCA1</italic> will be a good candidate, mutation of which disrupts the dependency of the circadian period on sucrose concentration in growth media (Haydon et al., <xref ref-type="bibr" rid="B21">2013</xref>). To examine this possibility, the single phase oscillator will be expanded to coupled two or multiple oscillators model, in which each phase oscillator is defined with different PRC to sucrose and sugar input function (<italic>Z</italic><sub>S</sub> and <italic>f</italic><sub>S</sub> in Equation (4), respectively). Alternative way is to use the clock gene-regulatory network model (Fogelmark and Troein, <xref ref-type="bibr" rid="B16">2014</xref>; De Caluw&#x000E9; et al., <xref ref-type="bibr" rid="B9">2016</xref>) with the explicit formalization of the interaction between clock genes and sugar signals.</p>
<p>We conclude that plants optimize growth by monitoring nutrient status and utilizing endogenous sugar signals as circadian-phase regulators. Photosynthetic products act not only as direct growth substances but also as feedback signals to the clock to realize the efficient carbon-usage for growth. It is plausible that plants as sessile organisms utilize signals from metabolism because they seem more controllable than environmental signals. In this sense, the phase regulation by endogenous signals is potentially even more important for plants than animals that also utilize metabolic signals as zeitgebers (Woller et al., <xref ref-type="bibr" rid="B44">2016</xref>). Our data strongly support the concept that the circadian clock improves the fitness of organisms by forming complex feedback loops with the signaling pathways it controls (Sanchez and Kay, <xref ref-type="bibr" rid="B33">2016</xref>).</p>
</sec>
<sec id="s5">
<title>Author contributions</title>
<p>TO and AS developed the model and wrote the manuscript. TO analyzed the model.</p>
<sec>
<title>Conflict of interest statement</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
</sec>
</body>
<back>
<ack><p>We thank Prof. Alex A. R. Webb, Dr. Camila Caldana, and Dr. Motohide Seki for fruitful discussions and useful comments on the manuscript.</p>
</ack>
<sec sec-type="supplementary-material" id="s6">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fpls.2017.01859/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fpls.2017.01859/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Presentation1.pdf" id="SM1" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Araya</surname> <given-names>T.</given-names></name> <name><surname>Noguchi</surname> <given-names>K.</given-names></name> <name><surname>Terashima</surname> <given-names>I.</given-names></name></person-group> (<year>2006</year>). <article-title>Effects of carbohydrate accumulation on photosynthesis differ between sink and source leaves of <italic>Phaseolus vulgaris</italic> L</article-title>. <source>Plant Cell Physiol.</source> <volume>47</volume>, <fpage>644</fpage>&#x02013;<lpage>652</lpage>. <pub-id pub-id-type="doi">10.1093/pcp/pcj033</pub-id><pub-id pub-id-type="pmid">16540483</pub-id></citation>
</ref>
<ref id="B2">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Baena-Gonz&#x000E1;lez</surname> <given-names>E.</given-names></name> <name><surname>Sheen</surname> <given-names>J.</given-names></name></person-group> (<year>2008</year>). <article-title>Convergent energy and stress signaling</article-title>. <source>Trends Plant Sci.</source> <volume>13</volume>, <fpage>474</fpage>&#x02013;<lpage>482</lpage>. <pub-id pub-id-type="doi">10.1016/j.tplants.2008.06.006</pub-id><pub-id pub-id-type="pmid">18701338</pub-id></citation>
</ref>
<ref id="B3">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Barillot</surname> <given-names>R.</given-names></name> <name><surname>Chambon</surname> <given-names>C.</given-names></name> <name><surname>Andrieu</surname> <given-names>B.</given-names></name></person-group> (<year>2016a</year>). <article-title>CN-Wheat, a functional&#x02013;structural model of carbon and nitrogen metabolism in wheat culms after anthesis. I. Model description</article-title>. <source>Ann. Bot.</source> <volume>118</volume>, <fpage>997</fpage>&#x02013;<lpage>1013</lpage>. <pub-id pub-id-type="doi">10.1093/aob/mcw143</pub-id></citation>
</ref>
<ref id="B4">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Barillot</surname> <given-names>R.</given-names></name> <name><surname>Chambon</surname> <given-names>C.</given-names></name> <name><surname>Andrieu</surname> <given-names>B.</given-names></name></person-group> (<year>2016b</year>). <article-title>CN-Wheat, a functional&#x02013;structural model of carbon and nitrogen metabolism in wheat culms after anthesis. II. Model evaluation</article-title>. <source>Ann. Bot.</source> <volume>118</volume>, <fpage>1015</fpage>&#x02013;<lpage>1031</lpage>. <pub-id pub-id-type="doi">10.1093/aob/mcw144</pub-id></citation>
</ref>
<ref id="B5">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Caspar</surname> <given-names>T.</given-names></name> <name><surname>Huber</surname> <given-names>S. C.</given-names></name> <name><surname>Somerville</surname> <given-names>C.</given-names></name></person-group> (<year>1985</year>). <article-title>Alterations in growth, photosynthesis, and respiration in a starchless mutant of <italic>Arabidopsis thaliana</italic> (L.) deficient in chloroplast phosphoglucomutase activity</article-title>. <source>Plant Physiol</source>. <volume>79</volume>, <fpage>11</fpage>&#x02013;<lpage>17</lpage>. <pub-id pub-id-type="doi">10.1104/pp.79.1.11</pub-id><pub-id pub-id-type="pmid">16664354</pub-id></citation>
</ref>
<ref id="B6">
<citation citation-type="other"><person-group person-group-type="author"><name><surname>Chew</surname> <given-names>Y. H.</given-names></name> <name><surname>Seaton</surname> <given-names>D. D.</given-names></name> <name><surname>Mengin</surname> <given-names>V.</given-names></name> <name><surname>Flis</surname> <given-names>A.</given-names></name> <name><surname>Mugford</surname> <given-names>S. T.</given-names></name> <name><surname>Smith</surname> <given-names>A. M.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>Linking circadian time to growth rate quantitatively via carbon metabolism</article-title>. Biorxiv <pub-id pub-id-type="doi">10.1101/105437</pub-id></citation>
</ref>
<ref id="B7">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chew</surname> <given-names>Y. H.</given-names></name> <name><surname>Wenden</surname> <given-names>B.</given-names></name> <name><surname>Flis</surname> <given-names>A.</given-names></name> <name><surname>Mengin</surname> <given-names>V.</given-names></name> <name><surname>Taylor</surname> <given-names>J.</given-names></name> <name><surname>Davey</surname> <given-names>C. L.</given-names></name> <etal/></person-group>. (<year>2014</year>). <article-title>Multiscale digital Arabidopsis predicts individual organ and whole-organism growth</article-title>. <source>Proc. Natl. Acad. Sci. U.S.A.</source> <volume>111</volume>, <fpage>4127</fpage>&#x02013;<lpage>4136</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.1410238111</pub-id><pub-id pub-id-type="pmid">25197087</pub-id></citation>
</ref>
<ref id="B8">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Christophe</surname> <given-names>A.</given-names></name> <name><surname>Letort</surname> <given-names>V.</given-names></name> <name><surname>Hummel</surname> <given-names>I.</given-names></name> <name><surname>Coum&#x000E8;de</surname> <given-names>P.-H.</given-names></name> <name><surname>de Reffye</surname> <given-names>P.</given-names></name> <name><surname>Lecoeur</surname> <given-names>J.</given-names></name></person-group> (<year>2008</year>). <article-title>A model-based analysis of the dynamics of carbon balance at the whole-plant level in <italic>Arabidopsis thaliana</italic></article-title>. <source>Funct. Plant Biol.</source> <volume>35</volume>, <fpage>1147</fpage>&#x02013;<lpage>1162</lpage>. <pub-id pub-id-type="doi">10.1071/FP08099</pub-id></citation>
</ref>
<ref id="B9">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>De Caluw&#x000E9;</surname> <given-names>J.</given-names></name> <name><surname>Xiao</surname> <given-names>Q.</given-names></name> <name><surname>Hermans</surname> <given-names>C.</given-names></name> <name><surname>Verbruggen</surname> <given-names>N.</given-names></name> <name><surname>Leloup</surname> <given-names>J.-C.</given-names></name> <name><surname>Gonze</surname> <given-names>D.</given-names></name></person-group> (<year>2016</year>). <article-title>A compact model for the complex plant circadian clock</article-title>. <source>Front. Plant Sci.</source> <volume>7</volume>:<fpage>74</fpage>. <pub-id pub-id-type="doi">10.3389/fpls.2016.00074</pub-id><pub-id pub-id-type="pmid">26904049</pub-id></citation>
</ref>
<ref id="B10">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Deprost</surname> <given-names>D.</given-names></name> <name><surname>Yao</surname> <given-names>L.</given-names></name> <name><surname>Sormani</surname> <given-names>R.</given-names></name> <name><surname>Moreau</surname> <given-names>M.</given-names></name> <name><surname>Leterreux</surname> <given-names>G.</given-names></name> <name><surname>Nicolai</surname> <given-names>M.</given-names></name> <etal/></person-group>. (<year>2007</year>). <article-title>The Arabidopsis TOR kinase links plant growth, yield, stress resistance and mRNA translation</article-title>. <source>EMBO Rep.</source> <volume>8</volume>, <fpage>864</fpage>&#x02013;<lpage>870</lpage>. <pub-id pub-id-type="doi">10.1038/sj.embor.7401043</pub-id><pub-id pub-id-type="pmid">17721444</pub-id></citation>
</ref>
<ref id="B11">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dodd</surname> <given-names>A. N.</given-names></name> <name><surname>Salathia</surname> <given-names>N.</given-names></name> <name><surname>Hall</surname> <given-names>A.</given-names></name> <name><surname>K&#x000E9;vei</surname> <given-names>E.</given-names></name> <name><surname>T&#x000F3;th</surname> <given-names>R.</given-names></name> <name><surname>Nagy</surname> <given-names>F.</given-names></name> <etal/></person-group>. (<year>2005</year>). <article-title>Plant circadian clocks increase photosynthesis, growth, survival, and competitive advantage</article-title>. <source>Science</source> <volume>309</volume>, <fpage>630</fpage>&#x02013;<lpage>633</lpage>. <pub-id pub-id-type="doi">10.1126/science.1115581</pub-id><pub-id pub-id-type="pmid">16040710</pub-id></citation>
</ref>
<ref id="B12">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Farr&#x000E9;</surname> <given-names>E. M.</given-names></name> <name><surname>Harmer</surname> <given-names>S. L.</given-names></name> <name><surname>Harmon</surname> <given-names>F. G.</given-names></name> <name><surname>Yanovsky</surname> <given-names>M. J.</given-names></name> <name><surname>Kay</surname> <given-names>S. A.</given-names></name></person-group> (<year>2005</year>). <article-title>Overlapping and distinct roles of PRR7 and PRR9 in the <italic>Arabidopsis circadian</italic> clock</article-title>. <source>Curr. Biol.</source> <volume>15</volume>, <fpage>47</fpage>&#x02013;<lpage>54</lpage>. <pub-id pub-id-type="doi">10.1016/j.cub.2004.12.067</pub-id><pub-id pub-id-type="pmid">15649364</pub-id></citation>
</ref>
<ref id="B13">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Feller</surname> <given-names>C.</given-names></name> <name><surname>Favre</surname> <given-names>P.</given-names></name> <name><surname>Janka</surname> <given-names>A.</given-names></name> <name><surname>Zeeman</surname> <given-names>S. C.</given-names></name> <name><surname>Gabriel</surname> <given-names>J.-P.</given-names></name> <name><surname>Reinhardt</surname> <given-names>D.</given-names></name></person-group> (<year>2015</year>). <article-title>Mathematical modeling of the dynamics of shoot-root interactions and resource partitioning in plant growth</article-title>. <source>PLoS ONE</source> <volume>10</volume>:<fpage>e0127905</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0127905</pub-id><pub-id pub-id-type="pmid">26154262</pub-id></citation>
</ref>
<ref id="B14">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Feugier</surname> <given-names>F. G.</given-names></name> <name><surname>Satake</surname> <given-names>A.</given-names></name></person-group> (<year>2013</year>). <article-title>Dynamical feedback between circadian clock and sucrose availability explains adaptive response of starch metabolism to various photoperiods</article-title>. <source>Front. Plant Sci.</source> <volume>3</volume>:<fpage>305</fpage>. <pub-id pub-id-type="doi">10.3389/fpls.2012.00305</pub-id><pub-id pub-id-type="pmid">23335931</pub-id></citation>
</ref>
<ref id="B15">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Figueroa</surname> <given-names>C. M.</given-names></name> <name><surname>Feil</surname> <given-names>R.</given-names></name> <name><surname>Ishihara</surname> <given-names>H.</given-names></name> <name><surname>Watanabe</surname> <given-names>M.</given-names></name> <name><surname>K&#x000F6;lling</surname> <given-names>K.</given-names></name> <name><surname>Krause</surname> <given-names>U.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Trehalose 6&#x02013;phosphate coordinates organic and amino acid metabolism with carbon availability</article-title>. <source>Plant J.</source> <volume>85</volume>, <fpage>410</fpage>&#x02013;<lpage>423</lpage>. <pub-id pub-id-type="doi">10.1111/tpj.13114</pub-id><pub-id pub-id-type="pmid">26714615</pub-id></citation>
</ref>
<ref id="B16">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fogelmark</surname> <given-names>K.</given-names></name> <name><surname>Troein</surname> <given-names>C.</given-names></name></person-group> (<year>2014</year>). <article-title>Rethinking transcriptional activation in the <italic>Arabidopsis circadian</italic> clock</article-title>. <source>PLoS Comput. Biol.</source> <volume>10</volume>:<fpage>e1003705</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pcbi.1003705</pub-id><pub-id pub-id-type="pmid">25033214</pub-id></citation>
</ref>
<ref id="B17">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gibon</surname> <given-names>Y.</given-names></name> <name><surname>Bl&#x000E4;sing</surname> <given-names>O. E.</given-names></name> <name><surname>Palacios-Rojas</surname> <given-names>N.</given-names></name> <name><surname>Pankovic</surname> <given-names>D.</given-names></name> <name><surname>Hendriks</surname> <given-names>J. H. M.</given-names></name> <name><surname>Fisahn</surname> <given-names>J.</given-names></name> <etal/></person-group>. (<year>2004</year>). <article-title>Adjustment of diurnal starch turnover to short days: depletion of sugar during the night leads to a temporary inhibition of carbohydrate utilization, accumulation of sugars and post-translational activation of ADP-glucose pyrophosphorylase in the following light period</article-title>. <source>Plant J.</source> <volume>39</volume>, <fpage>847</fpage>&#x02013;<lpage>862</lpage>. <pub-id pub-id-type="doi">10.1111/j.1365-313X.2004.02173.x</pub-id><pub-id pub-id-type="pmid">15341628</pub-id></citation>
</ref>
<ref id="B18">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Graf</surname> <given-names>A.</given-names></name> <name><surname>Smith</surname> <given-names>A. M.</given-names></name></person-group> (<year>2011</year>). <article-title>Starch and the clock: the dark side of plant productivity</article-title>. <source>Trends Plant Sci.</source> <volume>16</volume>, <fpage>169</fpage>&#x02013;<lpage>175</lpage>. <pub-id pub-id-type="doi">10.1016/j.tplants.2010.12.003</pub-id><pub-id pub-id-type="pmid">21216654</pub-id></citation>
</ref>
<ref id="B19">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Graf</surname> <given-names>A.</given-names></name> <name><surname>Schlereth</surname> <given-names>A.</given-names></name> <name><surname>Stitt</surname> <given-names>M.</given-names></name> <name><surname>Smith</surname> <given-names>A. M.</given-names></name></person-group> (<year>2010</year>). <article-title>Circadian control of carbohydrate availability for growth in Arabidopsis plants at night</article-title>. <source>Proc. Natl. Acad. Sci. U.S.A.</source> <volume>107</volume>, <fpage>9458</fpage>&#x02013;<lpage>9463</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.0914299107</pub-id><pub-id pub-id-type="pmid">20439704</pub-id></citation>
</ref>
<ref id="B20">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Greenham</surname> <given-names>K.</given-names></name> <name><surname>Lou</surname> <given-names>P.</given-names></name> <name><surname>Puzey</surname> <given-names>J. R.</given-names></name> <name><surname>Kumar</surname> <given-names>G.</given-names></name> <name><surname>Arnevik</surname> <given-names>C.</given-names></name> <name><surname>Farid</surname> <given-names>H.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>Geographic variation of plant circadian clock function in natural and agricultural settings</article-title>. <source>J. Biol. Rhythms</source> <volume>32</volume>, <fpage>26</fpage>&#x02013;<lpage>34</lpage>. <pub-id pub-id-type="doi">10.1177/0748730416679307</pub-id><pub-id pub-id-type="pmid">27920227</pub-id></citation>
</ref>
<ref id="B21">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Haydon</surname> <given-names>M. J.</given-names></name> <name><surname>Mielczarek</surname> <given-names>O.</given-names></name> <name><surname>Robertson</surname> <given-names>F. C.</given-names></name> <name><surname>Hubbard</surname> <given-names>K. E.</given-names></name> <name><surname>Webb</surname> <given-names>A. A. R.</given-names></name></person-group> (<year>2013</year>). <article-title>Photosynthetic entrainment of the <italic>Arabidopsis thaliana</italic> circadian clock</article-title>. <source>Nature</source> <volume>502</volume>, <fpage>689</fpage>&#x02013;<lpage>692</lpage>. <pub-id pub-id-type="doi">10.1038/nature12603</pub-id><pub-id pub-id-type="pmid">24153186</pub-id></citation>
</ref>
<ref id="B22">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Johnson</surname> <given-names>C. H.</given-names></name> <name><surname>Elliott</surname> <given-names>J. A.</given-names></name> <name><surname>Foster</surname> <given-names>R.</given-names></name></person-group> (<year>2003</year>). <article-title>Entrainment of circadian programs</article-title>. <source>Chronobiol. Int.</source> <volume>20</volume>, <fpage>741</fpage>&#x02013;<lpage>774</lpage>. <pub-id pub-id-type="doi">10.1081/CBI-120024211</pub-id><pub-id pub-id-type="pmid">14535352</pub-id></citation>
</ref>
<ref id="B23">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kaczorowski</surname> <given-names>K. A.</given-names></name> <name><surname>Quail</surname> <given-names>P. H.</given-names></name></person-group> (<year>2003</year>). <article-title>Arabidopsis pseudo-response regulator7 is a signaling intermediate in phytochrome-regulated seedling deetiolation and phasing of the circadian clock</article-title>. <source>Plant Cell</source> <volume>15</volume>, <fpage>2654</fpage>&#x02013;<lpage>2665</lpage>. <pub-id pub-id-type="doi">10.1105/tpc.015065</pub-id><pub-id pub-id-type="pmid">14563930</pub-id></citation>
</ref>
<ref id="B24">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lastdrager</surname> <given-names>J.</given-names></name> <name><surname>Hanson</surname> <given-names>J.</given-names></name> <name><surname>Smeekens</surname> <given-names>S.</given-names></name></person-group> (<year>2014</year>). <article-title>Sugar signals and the control of plant growth and development</article-title>. <source>J. Exp. Bot.</source> <volume>65</volume>, <fpage>799</fpage>&#x02013;<lpage>807</lpage>. <pub-id pub-id-type="doi">10.1093/jxb/ert474</pub-id><pub-id pub-id-type="pmid">24453229</pub-id></citation>
</ref>
<ref id="B25">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lobo</surname> <given-names>A. K. M.</given-names></name> <name><surname>Martins</surname> <given-names>M. O.</given-names></name> <name><surname>Neto</surname> <given-names>M. C. L.</given-names></name> <name><surname>Machado</surname> <given-names>E. C.</given-names></name> <name><surname>Ribeiro</surname> <given-names>R. V.</given-names></name> <name><surname>Silveira</surname> <given-names>J. A. G.</given-names></name></person-group> (<year>2015</year>). <article-title>Exogenous sucrose supply changes sugar metabolism and reduces photosynthesis of sugarcane through the down-regulation of Rubisco abundance and activity</article-title>. <source>J. Plant Physiol.</source> <volume>179</volume>, <fpage>113</fpage>&#x02013;<lpage>121</lpage>. <pub-id pub-id-type="doi">10.1016/j.jplph.2015.03.007</pub-id><pub-id pub-id-type="pmid">25863283</pub-id></citation>
</ref>
<ref id="B26">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Locke</surname> <given-names>J. C. W.</given-names></name> <name><surname>Southern</surname> <given-names>M. M.</given-names></name> <name><surname>Kozma-Bogn&#x000E1;r</surname> <given-names>L.</given-names></name> <name><surname>Hibberd</surname> <given-names>V.</given-names></name> <name><surname>Brown</surname> <given-names>P. E.</given-names></name> <name><surname>Turner</surname> <given-names>M. S.</given-names></name> <etal/></person-group>. (<year>2005</year>). <article-title>Extension of a genetic network model by iterative experimentation and mathematical analysis</article-title>. <source>Mol. Syst. Biol.</source> <volume>1</volume>:<fpage>0013</fpage>. <pub-id pub-id-type="doi">10.1038/msb4100018</pub-id><pub-id pub-id-type="pmid">16729048</pub-id></citation>
</ref>
<ref id="B27">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lu</surname> <given-names>Y.</given-names></name> <name><surname>Gehan</surname> <given-names>J. P.</given-names></name> <name><surname>Sharkey</surname> <given-names>T. D.</given-names></name></person-group> (<year>2005</year>). <article-title>Daylength and circadian effects on starch degradation and maltose metabolism</article-title>. <source>Plant Physiol</source>. <volume>138</volume>, <fpage>2280</fpage>&#x02013;<lpage>2291</lpage>. <pub-id pub-id-type="doi">10.1104/pp.105.061903</pub-id><pub-id pub-id-type="pmid">16055686</pub-id></citation>
</ref>
<ref id="B28">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mengin</surname> <given-names>V.</given-names></name> <name><surname>Pyl</surname> <given-names>E.-T.</given-names></name> <name><surname>Moraes</surname> <given-names>T. A.</given-names></name> <name><surname>Sulpice</surname> <given-names>R.</given-names></name> <name><surname>Krohn</surname> <given-names>N.</given-names></name> <name><surname>Encke</surname> <given-names>B.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>Photosynthate partitioning to starch in <italic>Arabidopsis thaliana</italic> is insensitive to light intensity but sensitive to photoperiod due to a restriction on growth in the light in short photoperiods</article-title>. <source>Plant Cell Environ.</source> <volume>40</volume>, <fpage>2608</fpage>&#x02013;<lpage>2627</lpage>. <pub-id pub-id-type="doi">10.1111/pce.13000</pub-id><pub-id pub-id-type="pmid">28628949</pub-id></citation>
</ref>
<ref id="B29">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Michael</surname> <given-names>T. P.</given-names></name> <name><surname>Salom&#x000E9;</surname> <given-names>P. A.</given-names></name> <name><surname>Yu</surname> <given-names>H. J.</given-names></name> <name><surname>Spencer</surname> <given-names>T. R.</given-names></name> <name><surname>Sharp</surname> <given-names>E. L.</given-names></name> <name><surname>McPeek</surname> <given-names>M. A.</given-names></name> <etal/></person-group>. (<year>2003</year>). <article-title>Enhanced fitness conferred by naturally occurring variation in the circadian clock</article-title>. <source>Science</source> <volume>302</volume>, <fpage>1049</fpage>&#x02013;<lpage>1053</lpage>. <pub-id pub-id-type="doi">10.1126/science.1082971</pub-id><pub-id pub-id-type="pmid">14605371</pub-id></citation>
</ref>
<ref id="B30">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>M&#x000FC;nch</surname> <given-names>E.</given-names></name></person-group> (<year>1930</year>). <source>Die Stoffbewegung in der Pflanze</source>. <publisher-loc>Jena</publisher-loc>: <publisher-name>Fischer</publisher-name>.</citation>
</ref>
<ref id="B31">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Quentin</surname> <given-names>A. G.</given-names></name> <name><surname>Close</surname> <given-names>D. C.</given-names></name> <name><surname>Hennen</surname> <given-names>L. M. H. P.</given-names></name> <name><surname>Pinkard</surname> <given-names>E. A.</given-names></name></person-group> (<year>2013</year>). <article-title>Down-regulation of photosynthesis following girdling, but contrasting effects on fruit set and retention, in two sweet cherry cultivars</article-title>. <source>Plant Physiol. Biochem.</source> <volume>73</volume>, <fpage>359</fpage>&#x02013;<lpage>367</lpage>. <pub-id pub-id-type="doi">10.1016/j.plaphy.2013.10.014</pub-id><pub-id pub-id-type="pmid">24189522</pub-id></citation>
</ref>
<ref id="B32">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ribeiro</surname> <given-names>R. V.</given-names></name> <name><surname>Machado</surname> <given-names>E. C.</given-names></name> <name><surname>Habermann</surname> <given-names>G.</given-names></name> <name><surname>Santos</surname> <given-names>M. G.</given-names></name> <name><surname>Oliveira</surname> <given-names>R. F.</given-names></name></person-group> (<year>2012</year>). <article-title>Seasonal effects on the relationship between photosynthesis and leaf carbohydrates in orange trees</article-title>. <source>Funct. Plant Biol.</source> <volume>39</volume>, <fpage>471</fpage>&#x02013;<lpage>480</lpage>. <pub-id pub-id-type="doi">10.1071/FP11277</pub-id></citation>
</ref>
<ref id="B33">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sanchez</surname> <given-names>S. E.</given-names></name> <name><surname>Kay</surname> <given-names>S. A.</given-names></name></person-group> (<year>2016</year>). <article-title>The plant circadian clock: from a simple timekeeper to a complex developmental manager</article-title>. <source>Cold Spring Harbor Perspect. Biol.</source> <volume>8</volume>:<fpage>a027748</fpage>. <pub-id pub-id-type="doi">10.1101/cshperspect.a027748</pub-id><pub-id pub-id-type="pmid">27663772</pub-id></citation>
</ref>
<ref id="B34">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Satake</surname> <given-names>A.</given-names></name> <name><surname>Seki</surname> <given-names>M.</given-names></name> <name><surname>Iima</surname> <given-names>M.</given-names></name> <name><surname>Teramoto</surname> <given-names>T.</given-names></name> <name><surname>Nishiura</surname> <given-names>Y.</given-names></name></person-group> (<year>2016</year>). <article-title>Florigen distribution determined by a source-sink balance explains the diversity of inflorescence structures in Arabidopsis</article-title>. <source>J. Theor. Biol.</source> <volume>395</volume>, <fpage>227</fpage>&#x02013;<lpage>237</lpage>. <pub-id pub-id-type="doi">10.1016/j.jtbi.2016.01.035</pub-id><pub-id pub-id-type="pmid">26845309</pub-id></citation>
</ref>
<ref id="B35">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Scialdone</surname> <given-names>A.</given-names></name> <name><surname>Mugford</surname> <given-names>S. T.</given-names></name> <name><surname>Feike</surname> <given-names>D.</given-names></name> <name><surname>Skeffington</surname> <given-names>A.</given-names></name> <name><surname>Borrill</surname> <given-names>P.</given-names></name> <name><surname>Graf</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2013</year>). <article-title>Arabidopsis plants perform arithmetic division to prevent starvation at night</article-title>. <source>Elife</source> <volume>2</volume>:<fpage>e00669</fpage>. <pub-id pub-id-type="doi">10.7554/eLife.00669</pub-id><pub-id pub-id-type="pmid">23805380</pub-id></citation>
</ref>
<ref id="B36">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Seki</surname> <given-names>M.</given-names></name> <name><surname>Feugier</surname> <given-names>F. G.</given-names></name> <name><surname>Song</surname> <given-names>X.-J.</given-names></name> <name><surname>Ashikari</surname> <given-names>M.</given-names></name> <name><surname>Nakamura</surname> <given-names>H.</given-names></name> <name><surname>Ishiyama</surname> <given-names>K.</given-names></name> <etal/></person-group>. (<year>2015</year>). <article-title>A mathematical model of phloem sucrose transport as a new tool for designing rice panicle structure for high grain yield</article-title>. <source>Plant Cell Physiol.</source> <volume>56</volume>, <fpage>605</fpage>&#x02013;<lpage>619</lpage>. <pub-id pub-id-type="doi">10.1093/pcp/pcu191</pub-id><pub-id pub-id-type="pmid">25516572</pub-id></citation>
</ref>
<ref id="B37">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Seki</surname> <given-names>M.</given-names></name> <name><surname>Ohara</surname> <given-names>T.</given-names></name> <name><surname>Hearn</surname> <given-names>T. J.</given-names></name> <name><surname>Frank</surname> <given-names>A.</given-names></name> <name><surname>da Silva</surname> <given-names>V. C. H.</given-names></name> <name><surname>Caldana</surname> <given-names>C.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>Adjustment of the <italic>Arabidopsis circadian</italic> oscillator by sugar signaling dictates the regulation of starch metabolism</article-title>. <source>Sci. Rep.</source> <volume>7</volume>:<fpage>8305</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-017-08325-y</pub-id></citation>
</ref>
<ref id="B38">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Smith</surname> <given-names>A. M.</given-names></name> <name><surname>Stitt</surname> <given-names>M.</given-names></name></person-group> (<year>2007</year>). <article-title>Coordination of carbon supply and plant growth</article-title>. <source>Plant Cell Environ.</source> <volume>30</volume>, <fpage>1126</fpage>&#x02013;<lpage>1149</lpage>. <pub-id pub-id-type="doi">10.1111/j.1365-3040.2007.01708.x</pub-id><pub-id pub-id-type="pmid">17661751</pub-id></citation>
</ref>
<ref id="B39">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Smith</surname> <given-names>S. M.</given-names></name> <name><surname>Fulton</surname> <given-names>D. C.</given-names></name> <name><surname>Chia</surname> <given-names>T.</given-names></name> <name><surname>Thorneycroft</surname> <given-names>D.</given-names></name> <name><surname>Chapple</surname> <given-names>A.</given-names></name> <name><surname>Dunstan</surname> <given-names>H.</given-names></name> <etal/></person-group>. (<year>2004</year>). <article-title>Diurnal changes in the transcriptome encoding enzymes of starch metabolism provide evidence for both transcriptional and posttranscriptional regulation of starch metabolism in Arabidopsis leaves</article-title>. <source>Plant Physiol</source>. <volume>136</volume>, <fpage>2687</fpage>&#x02013;<lpage>2699</lpage>. <pub-id pub-id-type="doi">10.1104/pp.104.044347</pub-id><pub-id pub-id-type="pmid">15347792</pub-id></citation>
</ref>
<ref id="B40">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Stitt</surname> <given-names>M.</given-names></name> <name><surname>Zeeman</surname> <given-names>S. C.</given-names></name></person-group> (<year>2012</year>). <article-title>Starch turnover: pathways, regulation and role in growth</article-title>. <source>Curr. Opin. Plant Biol.</source> <volume>15</volume>, <fpage>282</fpage>&#x02013;<lpage>292</lpage>. <pub-id pub-id-type="doi">10.1016/j.pbi.2012.03.016</pub-id><pub-id pub-id-type="pmid">22541711</pub-id></citation>
</ref>
<ref id="B41">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sulpice</surname> <given-names>R.</given-names></name> <name><surname>Flis</surname> <given-names>A.</given-names></name> <name><surname>Ivakov</surname> <given-names>A. A.</given-names></name> <name><surname>Apelt</surname> <given-names>F.</given-names></name> <name><surname>Krohn</surname> <given-names>N.</given-names></name> <name><surname>Encke</surname> <given-names>B.</given-names></name> <etal/></person-group>. (<year>2014</year>). <article-title>Arabidopsis coordinates the diurnal regulation of carbon allocation and growth across a wide range of photoperiods</article-title>. <source>Mol. Plant</source> <volume>7</volume>, <fpage>137</fpage>&#x02013;<lpage>155</lpage>. <pub-id pub-id-type="doi">10.1093/mp/sst127</pub-id><pub-id pub-id-type="pmid">24121291</pub-id></citation>
</ref>
<ref id="B42">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Thornley</surname> <given-names>J. H. M.</given-names></name> <name><surname>Johnson</surname> <given-names>I. R.</given-names></name></person-group> (<year>1990</year>). <source>Plant and Crop Modeling: A Mathematical Approach to Plant and Crop Physiology</source>. <publisher-loc>Oxford</publisher-loc>: <publisher-name>Clarendon Press</publisher-name>.</citation>
</ref>
<ref id="B43">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Woelfle</surname> <given-names>M. A.</given-names></name> <name><surname>Ouyang</surname> <given-names>Y.</given-names></name> <name><surname>Phanvijhitsiri</surname> <given-names>K.</given-names></name> <name><surname>Johnson</surname> <given-names>C. H.</given-names></name></person-group> (<year>2004</year>). <article-title>The adaptive value of circadian clocks: an experimental assessment in cyanobacteria</article-title>. <source>Curr. Biol.</source> <volume>14</volume>, <fpage>1481</fpage>&#x02013;<lpage>1486</lpage>. <pub-id pub-id-type="doi">10.1016/j.cub.2004.08.023</pub-id><pub-id pub-id-type="pmid">15324665</pub-id></citation>
</ref>
<ref id="B44">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Woller</surname> <given-names>A.</given-names></name> <name><surname>Duez</surname> <given-names>H.</given-names></name> <name><surname>Staels</surname> <given-names>B.</given-names></name> <name><surname>Lefranc</surname> <given-names>M.</given-names></name></person-group> (<year>2016</year>). <article-title>A mathematical model of the liver circadian clock linking feeding and fasting cycles to clock function</article-title>. <source>Cell Rep.</source> <volume>17</volume>, <fpage>1087</fpage>&#x02013;<lpage>1097</lpage>. <pub-id pub-id-type="doi">10.1016/j.celrep.2016.09.060</pub-id><pub-id pub-id-type="pmid">27760313</pub-id></citation>
</ref>
<ref id="B45">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yazdanbakhsh</surname> <given-names>N.</given-names></name> <name><surname>Sulpice</surname> <given-names>R.</given-names></name> <name><surname>Graf</surname> <given-names>A.</given-names></name> <name><surname>Stitt</surname> <given-names>M.</given-names></name> <name><surname>Fisahn</surname> <given-names>J.</given-names></name></person-group> (<year>2011</year>). <article-title>Circadian control of root elongation and C partitioning in <italic>Arabidopsis thaliana</italic></article-title>. <source>Plant Cell Environ.</source> <volume>34</volume>, <fpage>877</fpage>&#x02013;<lpage>894</lpage>. <pub-id pub-id-type="doi">10.1111/j.1365-3040.2011.02286.x</pub-id><pub-id pub-id-type="pmid">21332506</pub-id></citation>
</ref>
<ref id="B46">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Younginger</surname> <given-names>B. S.</given-names></name> <name><surname>Sirov&#x000E1;</surname> <given-names>D.</given-names></name> <name><surname>Cruzan</surname> <given-names>M. B.</given-names></name> <name><surname>Ballhorn</surname> <given-names>D. J.</given-names></name></person-group> (<year>2017</year>). <article-title>Is biomass a reliable estimate of plant fitness?</article-title> <source>Appl. Plant Sci.</source> <volume>5</volume>:<fpage>1600094</fpage>. <pub-id pub-id-type="doi">10.3732/apps.1600094</pub-id><pub-id pub-id-type="pmid">28224055</pub-id></citation>
</ref>
</ref-list>
<fn-group>
<fn fn-type="financial-disclosure"><p><bold>Funding.</bold> TO was supported by Grant-in-Aid for JSPS Fellows Number 16J04492. AS was supported by JSPS KAKENHI Number BR150301.</p>
</fn>
</fn-group>
</back>
</article>