Abstract

Photosynthetic starch reserves that accumulate in Arabidopsis leaves during the day decrease approximately linearly with time at night to support metabolism and growth. We find that the rate of decrease is adjusted to accommodate variation in the time of onset of darkness and starch content, such that reserves last almost precisely until dawn. Generation of these dynamics therefore requires an arithmetic division computation between the starch content and expected time to dawn. We introduce two novel chemical kinetic models capable of implementing analog arithmetic division. Predictions from the models are successfully tested in plants perturbed by a night-time light period or by mutations in starch degradation pathways. Our experiments indicate which components of the starch degradation apparatus may be important for appropriate arithmetic division. Our results are potentially relevant for any biological system dependent on a food reserve for survival over a predictable time period.

https://doi.org/10.7554/eLife.00669.001

eLife digest

Many organisms build up reserves of food when it is available so that metabolism and growth can continue when food is no longer available. Plants, for example, use energy from the sun to produce starch during the day, which is then used as a source of energy during the night. In some plants the amount of starch increases linearly with time during the day, and declines linearly with time during the night, so that the plant contains very little starch when the sun rises the following day. Although a plant is likely to starve if it cannot store or consume starch effectively, very little is known about the mechanisms that plants use to ensure that they store enough starch and do not use it up too quickly.

Plants are able to track the time to dawn using their internal circadian clock so, as Scialdone et al. point out, if they can also track how much starch they have stored, they might somehow be dividing the amount of starch by the time to dawn to work out the rate at which starch can be consumed so that it lasts until sunrise. But could plants actually perform such calculations?

To gain some insight into this puzzle, Scialdone et al. constructed mathematical models in which information about the size of the starch store and the time until dawn was encoded in the concentrations of two kinds of molecules (called S for starch and T for time). In one model, they propose that the S molecules stimulate starch consumption, and the T molecules prevent this from happening, with the rate of starch consumption being equal to the concentration of S molecules divided by the concentration of T molecules. These models were able to reproduce the results of previous experiments, including experiments in which dawn arrived unexpectedly early or unexpectedly late.

Scialdone et al. then performed experiments which confirmed predictions of the models for the pattern of starch consumption in plants lacking relevant enzymes. These experiments also revealed that an enzyme called PWD may be the point at which results of the division computation are integrated into the starch consumption pathway. More generally, this work shows that sophisticated arithmetic computations can be important in biology. Moreover, whereas computers rely on digital logic, Scialdone et al. show that arithmetic computations can also be performed by exploiting the analogue dynamics that take place between molecules.

https://doi.org/10.7554/eLife.00669.002

Introduction

Organisms must control the rate of consumption of their stored food reserves to prevent starvation during periods when food acquisition is not possible. A classic example of this requirement is provided by the response to the light/dark cycle in plants. During the day, plants utilize solar energy for carbon assimilation through photosynthesis. During the night, when solar energy is unavailable, plants utilize stored carbohydrate—usually starch—to allow continued metabolism and growth. In the model plant Arabidopsis thaliana this phenomenon is essential for productivity: mutants with defects in either the accumulation or the degradation of starch have reduced productivity and exhibit symptoms of starvation (Usadel et al., 2008; Yazdanbakhsh et al., 2011). The leaf starch content of Arabidopsis increases approximately linearly with time during the day: more than half of the carbon assimilated via photosynthesis may be stored as semi-crystalline starch granules inside chloroplasts. At night, starch content decreases approximately linearly with time such that 95% of starch is utilized by dawn (Gibon et al., 2004; Graf et al., 2010). This pattern of utilization is extremely robust, and is achieved even when darkness comes unexpectedly early (Graf et al., 2010). It is also likely to be optimal for the efficient utilization of carbohydrate over the light/dark cycle (Gibon et al., 2004; Smith and Stitt, 2007; Graf and Smith, 2011; Stitt and Zeeman, 2012; Feugier and Satake, 2013). However, despite the high importance for plant productivity of precise control of starch degradation, the way in which such dynamics are generated is very poorly understood. One intriguing possibility is the existence of a mechanism that dynamically measures the starch content and the expected time to dawn, then arithmetically divides these two quantities to compute the appropriate starch degradation rate. Such a mechanism could ensure complete utilization of available starch reserves at dawn despite variation in both the starch content at the onset of darkness and the subsequent duration of darkness. Consistent with this idea, we have recently shown that computation of the appropriate starch degradation rate in a normal night requires the circadian clock (Graf et al., 2010), which indicates how information about the expected time to dawn is obtained.

Although levels of transcripts for starch-degrading enzymes undergo large daily changes, levels of the enzymes themselves do not (Smith et al., 2004; Lu et al., 2005; Kötting et al., 2005; Yu et al., 2005; and our unpublished data). Therefore, it is likely that an arithmetic division mechanism would control flux through starch degradation at a post-translational level. Post-translational control would also permit swifter modulation of the catalytic capacity of these abundant proteins than would be possible via transcription/translation. Accordingly, in this paper we propose appropriate analog division mechanisms that operate through post-translational chemical kinetics. More generally, we also consider analog chemical kinetic schemes that generate addition, subtraction and multiplication operations. To the best of our knowledge the starch degradation system constitutes the first concrete realization of such arithmetic operations in biology. In this context, we therefore introduce two mathematical models which can both implement arithmetic division between the starch content and the expected time to dawn. We then successfully test predictions from the two models by examining the pattern of starch degradation in abnormal light/dark cycles and in a range of mutant plants defective in components of the starch degradation apparatus. Finally, our experiments also indicate which components of the starch degradation apparatus may be important for the appropriate implementation of arithmetic division.

Results

We first investigated the robustness of a potential arithmetic division calculation to perturbations in both the numerator (starch content) and denominator (expected time to dawn). Previously, we showed that the rate of starch degradation is appropriately adjusted in response to an unexpectedly early night (imposition of darkness 8 hr after dawn on plants grown in 12-hr light/12-hr dark cycles) (Graf et al., 2010). We found that adjustment also occurs in response to an unexpectedly late night. Plants grown in 12-hr light/12-hr dark cycles were subjected to darkness at either 12 hr or 16 hr after dawn. In both cases starch content decreased approximately linearly with time during the night, but with different slopes such that starch reserves were almost exhausted by dawn (Figure 1A). We also found that similar adjustments could be performed in a cca1/lhy circadian clock mutant which has a free-running period of <24 hr (Alabadi et al., 2001). In 12-hr light/12-hr dark cycles, this mutant degrades its starch by approximately 21 hr after dawn, rather than the normal 24 hr (Graf et al., 2010; Figure 1B).When subjected to an unexpected early night, the starch degradation rate in the mutant was adjusted, such that starch reserves were again exhausted at around 21 hr after dawn (Figure 1B). Appropriate adjustments of the rate of starch degradation also occurred in wild-type plants in which environmental manipulations produced different starch contents at the end of the 12-hr light period. A subset of a uniform batch of plants was transferred to a reduced light level for a single light period, leading to a twofold reduction in the starch content at the onset of darkness. For these and control plants subjected to normal light levels, starch content decreased approximately linearly with time during the night, but with different slopes such that starch reserves were almost exhausted by dawn in both cases (Figure 1C). Appropriate adjustment of starch degradation also occurred in subsets of plants exposed to three different regimes of varying light intensity over a single light period that generated two different starch contents at the onset of darkness (Figure 1—figure supplement 1). To investigate whether this phenomenon is widespread among plants, we examined the model grass Brachypodium distachyon. Ancestors of Arabidopsis and Brachypodium diverged at least 140 Myr ago. Starch content in Brachypodium leaves increased through the light period. As in Arabidopsis, the approximately linear decrease of starch with time following either a normal night (12 hr after dawn) or an unexpectedly early night (8 hr after dawn) was such that starch reserves were almost depleted by dawn (Figure 1D).

Figure 1 with 1 supplement see all
Starch content levels from experiments with unexpected variation in either starch content at the onset of darkness or the time of onset of darkness.

(A) Starch turnover in Arabidopsis grown in 12-hr light/12-hr dark, then subject to unexpected early (8 hr, n = 6 individual rosettes, circles) normal (12 hr, n = 6, squares) or unexpected late (16 hr, n = 5, triangles) onset of darkness. (B) Starch turnover in Arabidopsis cca1/lhy mutant grown in 12-hr light/12-hr dark, then subject to unexpected early (9 hr, circles), or normal (12 hr, squares) onset of darkness (n = 6–10). (C) Starch turnover in Arabidopsis exposed to different daytime light levels: 90 µmol quanta m−2 s−1 (open squares) or 50 µmol quanta m−2 s−1 (filled squares) (both n = 5, previously all plants grown in 12-hr light/12-hr dark with 90 µmol quanta m−2 s−1). (D) Starch turnover in Brachypodium grown in 12-hr light/12-hr dark, then subject to unexpected early (8 hr, circles) or normal (12 hr, squares) onset of darkness (both n = 6). Error bars are standard error of the mean throughout.

https://doi.org/10.7554/eLife.00669.003

Overall, these results demonstrate that the control of starch degradation at night to achieve almost complete consumption at the expected time of dawn can accommodate unexpected variation in the time of onset of darkness, starch content at the start of the night, and patterns of starch accumulation during the preceding day. Although the rate of degradation is different in a circadian clock mutant with an altered period from that in the wild-type, the capacity to adjust starch degradation in response to an unexpectedly early night is not compromised. It is also likely that the mechanism underlying this control is present in evolutionarily-distant species of plants.

Robust computation of the appropriate starch degradation rate in response to perturbations of both the expected time to dawn and starch content is clearly consistent with the implementation of arithmetic division. Since it is conceptually unclear how such a computation might be performed, we turned to mathematical modeling to generate possible mechanisms.

The fact that starch exists as large polymers that are mostly inaccessible within granule matrices means that a separate measure is likely required to provide information about the total amount of starch present at any given time. Hence, we assume the existence of a soluble molecule S whose concentration is proportional to the amount of starch in a granule. Since plants are able to adjust the rate of starch degradation according to variations in two independent quantities (the expected time to dawn and the amount of starch present), two separate species of molecule are clearly required. Therefore, we further assume the existence of a molecule T whose concentration encodes information about the expected time to dawn.

In our first model the T molecule concentration is proportional to the expected time to dawn, except during a period after dawn when its value must be reset (Figure 2B). To compute the appropriate degradation rate, the S and T molecule concentrations must therefore be arithmetically divided. We propose that computations of this form can be carried out most simply using analog chemical kinetics. As shown in Figure 2A, it is straightforward to perform addition, subtraction and multiplication operations. Subtraction can be implemented through efficient sequestration and multiplication by a two-species chemical reaction. Division is slightly more intricate, but can be implemented using the model introduced in Figure 2D (other conceptually similar models are discussed in the ‘Materials and methods’). Here, S molecules associate with the starch granule surface, where they permit starch degradation (presumably in combination with other elements, as shown in Figure 2D) and where the S molecules can also be degraded. T molecules inhibit S molecule and starch degradation by binding to S and causing its detachment from the granule surface. In this way, it can be seen intuitively that a division-like operation can be implemented (for rigorous calculations, see ‘Materials and methods’). In Figure 2F,H,J, we show the best fit of this model to the data from Figure 1A–C with good results.

Chemical kinetic models capable of implementing analog arithmetic operations.

(A) Pictorial summaries of schemes for analog implementation of addition, subtraction and multiplication between the concentrations of two molecules S and T. Square brackets indicate concentrations. (B) and (C) Schematic behavior of the stromal concentrations of S and T molecules ([SC] and [TC] respectively), in (B) first and (C) second arithmetic division models. In the first model, the T molecule tracks the time to expected dawn after a reset-time tr. In the second model the T molecule concentration increases with time proportionally to 1/(expected time to dawn) between tr1 and tr2. (D) and (E) Pictorial summaries of (D) first and (E) second analog arithmetic division models (not all reactions shown in pictures, for full details see ‘Materials and methods’). In the reaction schemes, molecules not attached to the starch granule surface have a ‘C’ subscript. The blue disk represents components of the starch degradation apparatus potentially activated by the S molecule in the first model, and by the ST complex in the second model. Best fits (full lines) of first (F), (H), and (J) and second (G), (I), and (K) arithmetic division models to Arabidopsis data from Figure 1A–C.

https://doi.org/10.7554/eLife.00669.005

A second distinct possibility also exists for computing the appropriate degradation rate. We now assume that the T molecule concentration increases as the expected time of dawn approaches, before being reset. If this increase is such that the T molecule concentration is approximately proportional to 1/(expected time to dawn) (Figure 2C) then the appropriate degradation rate can be computed by multiplying the S and T concentrations. This is implemented by the reaction scheme shown in Figure 2E: S molecules associate with the starch granule surface, where they recruit T molecules from the stroma. The resulting molecular complex permits the degradation of starch (again presumably in combination with other elements, as shown in Figure 2E) and of the S molecule itself. The output of this second model is very similar to that of the first model (fits to data of Figure 1A–C shown in Figure 2G,I,K) in that a division computation is still performed, but now the timing information is encoded differently in the T molecule concentration. One potential difficulty is the need to pre-compute the reciprocal of the expected time to dawn. A simple possible scheme to achieve this goal is outlined in the ‘Materials and methods’. Of course, a combination of the two above models is also possible, involving both multiplication and division by factors dependent on the expected time to dawn, such that overall an appropriate division computation is still performed. However, the additional complexity required for such implementations makes such a combined model less likely.

Importantly we assume that the granule surface area does not limit the reaction rate as the granule shrinks, consistent with an approximately linear decrease of starch content with time. Accordingly, in both models the degradation reactions (for both the starch and the S molecules) occur only in region(s) of granule surface of overall fixed area as each granule shrinks. This would be the case if one or more of the additional components required for starch degradation (illustrated in Figure 2D,E) is present at a fixed number on the granule surface as the granule shrinks. Clearly the assumption of a non-limiting surface area cannot remain true if the granule shrinks to very small volumes, as could happen at the end of the night. However, our models still fit the experimental data well even at these times (see Figure 2F–K).

Taken together, our modeling results show that arithmetic division (as well as other arithmetic operations) can be implemented simply using analog post-translational chemical kinetics. Furthermore, output from both models fits our experimental data well. We note that there is some variation in the fitted parameters for both models (see Tables 3–6), which arises due to variation between experimental data sets even for a single genotype. Such variation is widely observed for measurements of primary metabolites over the light/dark cycle and could arise from batch to batch variation in expression levels of degradation or clock components.

As the fit of the two models to the data is equally good, distinguishing between them is currently challenging. However, we can test two critical predictions common to both models. The first prediction is that the degradation rate is continuously computed via arithmetic division during the night. Such a scheme clearly has advantages in its flexibility and potential to recover from unexpected perturbations. To test whether this prediction is correct, we interrupted a normal night with a period of light, ending 5 hr before the expected time of dawn. During this night-time light period starch accumulated to levels very similar to those present at the end of the day-time light period. Comparison of the rate of starch degradation following the night-time light period with the rate before this period, and with the rate that would have been expected over a normal 12-hr night, allowed a robust assessment of whether the degradation rate had been reset. We confirmed that the night-time light period did not re-entrain the circadian clock, by monitoring expression of the clock gene LHY (see ‘Materials and methods’ and Figure 3—figure supplement 1). For three independent experiments we found that the rates of starch degradation immediately following the night-time light period (between 19 hr and 21 hr after dawn) were significantly greater than both the corresponding actual rates before this period (between 12 hr and 14 hr after dawn) and the corresponding rates that would have been expected during a normal 12-hr night, see insets in Figure 3A–C (p = 1.6 × 10−3 and 3.8 × 10−5, respectively, details of the statistical analysis in the ‘Materials and methods’). The latter rates are those that would have ensured the complete depletion of the starch content measured at 12 hr at the expected time of dawn (24 hr). We obtained the same result by changing the duration and the starting time of the night-time light period (Figure 3D). This increase in starch degradation rate is consistent with our prediction that the rate is continuously computed during the night, rather than set only once at the first onset of darkness. In Figure 3 we show the best fits of the first model to the data, with good results. As before, the second model produced very similar fits (see Figure 3—figure supplement 2).

Figure 3 with 2 supplements see all
Starch content levels from experiments incorporating night-time light period.

Arabidopsis plants grown in 12-hr light/12-hr dark were subjected to onset of darkness at 12 hr, followed by an unexpected period of light, followed by extended darkness. (A)–(C) Three data sets (n = 12 individual rosettes, except n = 10 for C), in which the unexpected period of light was between 14 hr and 19 hr after dawn. (D) In the fourth dataset (n = 12) the period of light was between 16 hr and 20 hr after dawn. Full lines are best fits to the first division model. The second model produces very similar fits (see Figure 3—figure supplement 2). The insets show the respective starch degradation rates computed from the 12-hr and 14-hr experimental time points (dark grey bars) compared to those computed from the 19-hr and 21-hr experimental time points in panels (AC) or the 20-hr and 22-hr time points in panel (D) (light grey bars). The white bars are the expected starch degradation rates in a normal 12 hr night, that is rates that would have ensured the complete depletion of the starch content measured at 12 hr at the time of expected dawn (24 hr). Error bars are standard error of the mean throughout.

https://doi.org/10.7554/eLife.00669.006

The second prediction from the models concerns the effects on starch degradation of perturbations in components of the arithmetic division or degradation apparatus. To reproduce the normal pattern—in which starch is almost completely degraded by the expected time of dawn—our models require that a factor involved in specifying the relative degradation rates of the starch and S molecules must be fine-tuned to one (see ‘Materials and methods’). If instead this factor is <1, so that the starch is degraded too slowly, then rather than complete degradation the models predict that only a certain percentage of the starch will be consumed during the night regardless of the starch content at the end of the preceding light period (for full calculations, see ‘Materials and methods’). Accordingly, we predict that perturbations to parts of the arithmetic division or degradation apparatus might result in degradation of only a fraction of starch during the night. Interestingly, several mutants defective in proteins involved in, or related to, starch degradation exhibit approximately linearly decreasing starch content with time during the night, but have an elevated starch content at the end of the night. These include beta-amylase mutants bam3 and bam4, the debranching enzyme mutant isa3, and mutants lacking phosphoglucan water dikinase (pwd, also called gwd3), a glucan phosphate phosphatase (sex4) and a glucan phosphate phosphatase-like protein (lsf1) (Smith, 2012). We showed previously that this abnormal pattern of starch turnover was rapidly regained in lsf1 mutants that were transferred to normal light/dark cycles after starch was reduced to very low levels by prolonged darkness (Comparot-Moss et al., 2010). We re-analyzed these data and also performed a similar experiment on the sex4 mutant. For both mutants, we found that the fraction of end-of-light period starch content degraded during the night was approximately the same on successive nights (around 30% for the sex4 mutant and 45% for the lsf1 mutant), regardless of the starch content at the end of the respective preceding light period. This resulted in a progressive increase in the end-of-night starch content to an approximately constant value over 3–4 days after return to normal light/dark cycles (Figure 4A,B). Thus the pattern of starch degradation in these mutants is as predicted from the models for a situation in which the above factor is incorrectly set to a value <1.

Figure 4 with 1 supplement see all
Starch content levels in mutant Arabidopsis plants defective in components of the starch degradation apparatus.

(A) Starch content in wild-type (WT) plants and lsf1 and sex4 mutant plants during four days of 12-hr light/12-hr dark following 5 days of continuous darkness, where plants were transferred back into the light (at time 0 hr on the x-axis) 132 hr after the end of the previous light period (n = 6 individual rosettes). Data for wild-type and lsf1 plants are from (Comparot-Moss et al., 2010). (B) The percentage of starch degraded during each of the four nights in (A). (C)–(E) Starch content in lsf1, sex4 and pwd mutant plants grown in 12-hr light/12-hr dark cycles then subject to unexpected early (8 hr, circles) or normal (12 hr, squares) onset of darkness (n = 5). The continuous and dashed lines are linear fits to the normal and early night datasets respectively. (F) For each of the labeled genotypes, R is the ratio between the starch degradation rates (each normalized by their respective end-of-light period starch content and as determined from the linear fits) during the normal and early nights. The dashed line shows the expected value of R for wild-type (WT) plants, that is, ratio of rates that would ensure the complete depletion of the starch content in all cases at the time of expected dawn (24 hr). See ‘Materials and methods’ for details about the linear fitting and the calculation of R. Error bars are standard error of the mean throughout. Figure 4—figure supplement 1 shows the datasets used to calculate R for WT, bam3, bam4 and isa3.

https://doi.org/10.7554/eLife.00669.009

We next attempted to gain experimental insight into how the arithmetic division mechanism modulates flux through the starch degradation pathway. In our models the computation of the starch degradation rate is the result of interactions between the molecules encoding the information about starch content and time to dawn (S and T) and the starch degradation apparatus. If this idea is correct, we expect that mutants lacking components of the degradation apparatus involved in these interactions may also lack the ability to adjust the rate of starch degradation in response to variation in starch content or time of onset of darkness.

To look for a mutant in this category, we imposed an unexpectedly early night on six mutants, each lacking a protein involved in starch degradation (lsf1, sex4, bam3, bam4, isa3, pwd). All six mutants show approximately linear decreases in starch content with time during the night, but they have higher end-of-night starch contents than wild-type plants. This mutant collection covers the majority of currently-known components of the chloroplastic starch degradation apparatus. For five of the six mutants, the rates of starch degradation (as determined from linear fits) were lower during an unexpectedly early night (starting 8 hr after dawn) than in a normal night (starting 12 hr after dawn), as is the case in wild-type plants (Figure 4C,D and Figure 4—figure supplement 1). To quantify the adjustment of the starch degradation rates in the mutants, we calculated the ratio R between the degradation rates (normalized by the respective end-of-light period starch content) during the normal and the early night. For wild-type plants, the expected value of R is 16/12 ≈ 1.33, that is, the ratio between the expected lengths of the early and the normal night. In Figure 4F we show that the values of R for five of the six mutants are consistent with the wild-type value.

However for the sixth mutant, pwd, the rate of starch degradation was not adjusted in response to an unexpectedly early night (Figure 4E,F), as we found an R which was significantly different from the wild-type value (p = 0.01, details of the statistical analysis in the ‘Materials and methods’). This finding indicates that PWD—phosphoglucan, water dikinase—may be a node at which information about expected time to dawn and starch content is integrated to set an appropriate rate of starch degradation.

PWD contributes to a cycle of phosphorylation and dephosphorylation of glucosyl residues within starch polymers that is essential for normal starch degradation. Two enzymes—glucan water dikinase (GWD) and PWD—add phosphate groups to starch polymers at the granule surface, and two further enzymes SEX4 and LSF2—remove them. All four enzymes are necessary for normal rates of starch degradation, and loss of GWD almost completely prevents degradation at night (Smith, 2012). The phosphate groups are thought to disrupt the ordered packing of the starch molecules at the granule surface, allowing access to starch degrading enzymes (beta-amylases and isoamylase 3). Subsequent removal of the phosphate groups is essential for full degradation of the starch polymers because the phosphates block the action of the exo-acting beta-amylases (Baunsgaard et al., 2005; Kötting et al., 2005; Ritte et al., 2006; Hejazi et al., 2009). The cycle as a whole is thus an attractive candidate for integration of factors that modulate flux through starch degradation.

Because the phosphorylation/dephosphorylation cycle modifies the granule surface, it is a potential candidate not only for flux modulation but also for the storage of information about starch content. To discover whether phosphate groups may provide information about starch content, we measured the amount of granule-bound phosphate (on the 6-position of glucosyl residues, representing about 80% of the total phosphate) over the light/dark cycle. If phosphorylation simply tracks starch polymer synthesis during the day, phosphate levels per unit mass of starch will be constant and will thus contain no information about starch content. Surprisingly, however, we found a large increase/decrease in the level of phosphate per unit mass of starch over the light/dark cycle (Figure 5). This discovery implies that the S molecule may be a modulator of activities of enzymes of the phosphorylation/dephosphorylation cycle, thus generating a daily pattern of change in the accessibility of the granule surface to hydrolytic enzymes that approximately tracks starch content.

Daily change in starch phosphate content (measured as glucose 6-phosphate, G6P) in Arabidopsis leaves.

Results are normalized by total amount of glucose (Glc) in starch at each time point. Starch was extracted from rosettes of 26-day-old plants. n = 3 pools of 10 rosettes except at 24 hr time point, with n = 2 pools of 15 rosettes. Error bars represent the range (i.e., error bar edges correspond to highest and lowest values measured).

https://doi.org/10.7554/eLife.00669.011

Discussion

Overall, our work provides a new framework and perspective for understanding the control of reserve utilization in plants at night. Our experiments provide strong support for an implementation of arithmetic division in night-time starch degradation. We used mathematical modeling to generate two simple mechanisms capable of analog implementation of such an operation. Predictions from the models were then verified and a potential point identified at which the division computation could be integrated into the starch degradation pathway. We also showed that the phosphorylation state of the starch granule surface could provide information about starch content through the day.

Our analysis may also be relevant to a broader class of biological processes, where food reserves accumulated in advance of periods of predictable length without further food intake are just sufficient to permit survival to the expected end of the period. For example, migrating little stints (Calidris minuta) arriving at their Arctic summer breeding grounds after a 5000-km journey have sufficient remaining lipid reserves for an average of only 0.6 days (Tulp et al., 2009). During the 4-month fast period of egg-incubating male emperor penguins (Aptenodytes forsteri), lipid reserves are used such that they reach a critical depletion level at approximately the point at which the females are expected to return. Unexpected extension of the fast period leads to catabolism of protein and abandonment of offspring in favor of hunting for food (Groscolas and Robin, 2001). As in Arabidopsis leaves, the rate of reserve utilization in these examples can potentially be computed by arithmetically dividing the reserve levels by the anticipated time of fasting.

It is a longstanding idea that cells are able to use proteins to store and process information through networks of interactions (Bray, 1995). Understanding how such biochemical networks work and what kind of computations they perform is an ongoing challenge (see Deckard and Sauro, 2004; Lim et al., 2013). Our analysis here has underlined the utility of analog chemical kinetics in performing arithmetic computations in biology. Importantly, we have for the first time provided a concrete example of a biological system where such a computation is of fundamental importance. This contrasts to previous work where elegant theoretical implementations of arithmetic operations lacked specific biological applications (Cory and Perkins, 2008; Buisman et al., 2009). Analog chemical kinetic approaches may potentially also be useful for calculations in synthetic biology applications (Benenson, 2012), where they are likely to prove much simpler to implement than alternative schemes based on much more complex digital circuitry.

Materials and methods

Plant material

Request a detailed protocol

Arabidopsis thaliana (in the Col0 background, except for cca1/lhy and its wild-type which were in the Ws background) and Brachypodium distachyon (Bd21) plants were grown as in Graf et al. (2010) on soil in 12-hr light/12-hr dark with 200 µmol quanta m−2 s−1 illumination at a constant temperature of 20°C for 21 days.

Measurement of starch

Request a detailed protocol

Plants were harvested and extracted in dilute perchloric acid for analysis of starch, which was then quantified enzymatically as previously described (Graf et al., 2010).

Transcript analysis

Request a detailed protocol

RNA was extracted from plant material and qPCR was performed as in Graf et al. (2010). Oligonucleotide primer sequences were as follows:

PrimerSequence 5′ to 3′
LHY-F AT1G01060GACTCAAACACTGCCCAGAAGA
LHY-RAT1G01060CGTCACTCCCTGAAGGTGTATTT
ACT2-F AT3G18780ACTTTCATCAGCCGTTTTGA
ACT2-R AT3G18780ACGATTGGTTGAATATCATCAG

Starch phosphate measurements

Request a detailed protocol

Starch granules were prepared as in Ritte et al. (2000) and used immediately without drying. Granules were resuspended in two pellet volumes of water and boiled for 15 min then digested with 20 U amylogucosidase (Roche) and 2 U α-amylase (Megazyme) in 100 mM Na acetate pH 4.8 for 9 hr at 37°C. Glucose was assayed enzymatically following (Hargreaves and ap Rees, 1988) and glucose 6-phosphate was measured enzymatically using the fluorimetric assay of (Zhu et al., 2009).

Mathematical Modelling

Quantification of the starch content and the expected time to dawn

Request a detailed protocol

Starch is laid down as semi-crystalline granules in chloroplasts in the light. Granule surfaces are then subject to degradation during periods of darkness. We assume here that the computation of the degradation rate is performed autonomously inside each chloroplast. In order to correctly compute the appropriate starch degradation rate, knowledge of the total chloroplast starch content is required. We assume that there exists a molecule S whose number NStot is proportional to the total starch content ΔStot in a chloroplast at the end of a light period:

ΔStot=αNStot.

We also assume that the S molecule can exchange rapidly between the granule surfaces and the surrounding chloroplast compartment. Moreover, we propose the existence of a molecule T that encodes information about the expected time to dawn. The T molecule dynamics will presumably be controlled by the circadian clock and are assumed to be independent of the starch degradation process.

In the first model (Figure 2B,D), the T concentration in the compartment surrounding the granules, [TC], is proportional to the expected time to dawn Δt = tdayt during most of the light/dark cycle. Note, however, that this cannot be true at all times, as around the time of expected dawn, after [TC] has dropped to low levels, its value must be reset. So, if dawn occurs at t = 0, we assume that:

(1) [TC]={βtdaytrtrt0ttrβ(tdayt)trttday,

where β is a constant, tday is the period of the light/dark cycle and tr is the reset-time after which [TC] starts tracking the time to expected dawn (see Figure 2B).

These dynamics could be the result of production of T at a constant rate, with a comparatively low degradation rate for 0 ≤ t ≤ tr. During tr ≤ t ≤ tday, [TC] is assumed to decrease linearly with time, with a rate β. This could happen through an efficient sequestration of TC by another molecule which is produced at a rate β and is degraded at a comparatively low rate for tr ≤ t ≤ tday. These simple [TC] dynamics can reproduce the available experimental data on starch degradation during the hours after dawn. Experimentally, if the onset of darkness unexpectedly occurs only 6 hr after the previous dawn in plants previously grown in 12-hr light/12-hr dark conditions, the starch is degraded too quickly, while for an unexpected onset of darkness at 8 hr, the starch degradation occurs at the appropriate rate (Lu et al., 2005; Graf et al., 2010). This can be explained by the above model, with the T dynamics given in Equation 1 and assuming a reset time tr ≳ 8 hr, as we have [TC] < βΔt during the resetting period, and hence the starch degradation rate will initially be too high.

Nevertheless, if the early night condition is maintained each day for a few days in succession, acclimation occurs: the appropriate starch degradation rate is computed even after a light period as short as 4 hr. This suggests that tr can be reduced, thereby resetting [TC] more quickly, and hence allowing the arithmetic division mechanism to work effectively under the new light/dark cycle conditions.

In the second model (Figure 2C,E), we assume that for most of the light/dark cycle [TC] is proportional to 1/(tdayt). Again, this cannot be true around dawn, where [TC ] must be reset. For the sake of simplicity, we assume that during the reset period [TC] has a linear profile with time:

(2) [TC]={βtday+tr1tr2(t+tdaytr2tdaytr1+tr1ttdaytr2)0ttr1βtdayttr1ttr2βtday+tr1tr2(ttr2tdaytr1+tday+tr1ttdaytr2)tr2ttday,

where β is a constant, and tr1 and tr2 define the time interval in which [TC] is proportional to 1/Δt (see Figure 2C). In the second model, these dynamics can also explain why starch degradation is too rapid following a very early and unexpected onset of darkness. During the resetting period after dawn [TC] > β/Δt, and therefore the starch degradation rate will initially be too high.

One apparent difficulty with the second model is in finding a simple way to generate a T molecule concentration that scales as 1/Δt. In fact, kinetics of this kind are simple to produce. For example, if T promotes further production of the T molecule itself, with a quadratic nonlinearity, then, away from saturation, we have d[TC]dt=η[TC]2. It is straightforward to show that the solution to this equation is:

[TC](t)=[TC](tr1)1η(ttr1)[TC](tr1).

For [TC](t) = β/Δt as required, we then find [TC](tr1) = β/(tdaytr1) and η = 1/β.

Of course, more complex schemes could also generate a similar behavior; we simply wish to point out here that a 1/Δt behavior for [TC] is a plausible assumption.

Two models that appropriately calculate the starch degradation rate

We now describe in full mathematical detail the models introduced in the main text.

Model 1

Request a detailed protocol

In this model, division is implemented by dividing the S concentration by the T concentration, where the T molecule concentration is proportional to Δt for most of the light/dark cycle. Such a division can be implemented by the following set of reactions:

(3) rS:SCSrST1:S+TCSTrST2:STSC+TCrD1:SrD2:S+StarchS,

where the subscript ‘C’ refers to molecules in the compartment surrounding the granules. The reversible rS reaction describes exchange of the S molecules between the granule surface and surrounding compartment, with forward reaction parameter fS and backward rate bS. The T molecules can be recruited by the S molecules on the granule surface (reaction rST1, reaction parameter fST1) and the resulting complex can then dissociate leading to the detachment of S and T from the granule surface (reaction rST2, rate fST2). The S molecule on the granule can be degraded (reaction rD1, reaction parameter fD1) or it can permit starch degradation (reaction rD2, reaction parameter fD2).

In all the calculations in this section, we assume that bSfD1 and that the dynamics of the S and T molecules can be taken to be in quasi-steady-state. As described in the main text, we assume that the degradation reactions (for both the starch and the S molecule) can only occur in a region of overall fixed area Ad as each granule shrinks. We also assume that all granules in a given chloroplast are approximately equal in area and volume. The differential equations describing the dynamics of the total number of S molecules, NStot, and the total amount of starch in a chloroplast, ΔStot, are:

(4) dNStotdt=ndNSdt=nfD1[S]AddΔStotdt=ndΔSdt=nmSfD2[S]Ad=nμ=μtot,

where mS is the amount of starch degraded by an S molecule in a single degradation reaction, n is the number of granules in the chloroplast, μ is the starch degradation rate for an individual granule and μtot is the total starch degradation rate.

Using Equation 4, we find that dΔStotdt=μtot=mSfD2fD1dNStotdt. We assumed that at the beginning of the dark period ΔStot=αNStot. It is easy to show that, if α=mSfD2fD1, such a proportionality between the starch content and the number of S molecules continues to hold true throughout the degradation process. The total number of S molecules is NStot=[SC]V+n[S]A+n[ST]A, where V is the compartmental volume surrounding the granules, and A is the surface area of a granule. Using the above quasi-steady state assumption and the law of mass action, we obtain:

(5) [S]=kS[SC]1+kST1[TC],[ST]=fST1fST2kS[SC][TC]1+kST1[TC],

with kS = fS/bS, kST1 = fST1/bS. Assuming that nkS AV, bS/fST2 ≲1, then n[S] A, n[ST]A ≪ [SC]V, and hence:

(6) [SC]=NStotV=NStotnV0,

where we used the observation that the number of starch granules per unit volume inside the chloroplast is approximately constant, that is, the number of granules n is proportional to the chloroplast volume V (Crumpton-Taylor et al., 2012). For the T molecule, assuming the concentration in the compartment surrounding the granules is regulated by the circadian clock, we have (for much of the light/dark cycle, see Equation 1):

(7) [TC]=βΔt.

The starch degradation rate for an individual granule is then:

μ=mSfD2[S]Ad=mSfD2AdkS[SC]1+kST1[TC].

According to Equation 6, [SC] is proportional to the amount of starch in a granule, whereas from Equation 7, [TC] is proportional to the expected time to next dawn. Sufficiently far from the time of expected dawn when kST1 [TC] ≫ 1, we then have:

(8) μ=mSfD2AdkSNStotnV0kST1βΔt.

Hence, summing over all the granules, and using that ΔStot=αNStot=mSfD2fD1NStot (see above in this section), the total starch degradation rate is:

dΔStotdt=μtot=fD1AdkSkST1V0βΔStotΔt.

According to this equation, starch contents are completely depleted by the time of expected dawn, and for a degradation rate μtot=ΔStotΔt with normalization relation fD1AdkSkST1V0β=1, starch contents decrease linearly with time.

Model 2

Request a detailed protocol

In this model, division is now implemented by multiplying the T concentration by the S concentration, where the T molecule concentration is proportional to 1/Δt for most of the light/dark cycle. Such a multiplication can be implemented by the following set of reactions:

rS:SCSrST1:S+TCSTrST2:STS+TCrD1:STTCrD2:ST+StarchST.

Reversible reaction rS again describes the exchange of S molecules between the surrounding compartment and the granule surface (forward reaction parameter fS and backward reaction parameter bS). S can recruit T molecules from the surrounding compartment and form the complex ST (reaction rST1, reaction parameter fST). As the ST complex dissociates, T is released back in the stroma, while S remains attached to the granule (reaction rST2, reaction parameter bST). The ST complex permits S degradation, with T released from the granule surface (reaction rD1, reaction parameter fD1), as well as starch degradation (reaction rD2, reaction parameter fD2). Similar to our assumptions for Model 1, we assume here that bSTfD1, with a fixed degradation area Ad for the starch as well as for the S molecule as each granule shrinks. We also again assume that the dynamics of S and T can be taken to be in quasi-steady-state and that all granules in a given chloroplast have approximately equal areas and volumes.

The differential equations describing the dynamics of the total amount of starch in a chloroplast, ΔStot, and of the total number of S molecules, NStot, are now:

dNStotdt=ndNSdt=nfD1[ST]AddΔStotdt=ndΔSdt=nmSfD2[ST]Ad=nμ=μtot.

The above relations again ensure that if α=mSfD2fD1, then ΔStot=αNStot at all times during degradation. The total number of S molecules is NStot=[SC]V+n[S]A+n[ST]A, where V is the compartmental volume surrounding the granules, and where A is the surface area of a granule. Using the above quasi-steady state assumption and the law of mass action, we obtain:

[S]=kS[SC],[ST]=kSTkS[SC][TC],

with kS = fS/bS and kST = fST/bST. We assume that nkS AV, kST[TC] ≲ 1, giving n[S]A, n[ST]A ≪ [SC]V, and hence:

[SC]=NStotV=NStotnV0,

where we again used the observation that the number of starch granules per unit volume inside the chloroplast is approximately constant, that is, the number of granules n is proportional to the volume V. For the T molecules, we have (for much of the light/dark cycle, see Equation 2):

[TC]=β/Δt.

The starch degradation rate for an individual granule is therefore:

μ=mSfD2[ST]Ad=mSfD2AdkSTkS[SC][TC].

According to the above equation, [TC] is proportional to the reciprocal of the expected time to the next dawn. Thus, we have:

μ=mSfD2AdkSTkSβNStotnV0Δt.

Summing over all the granules, and using ΔStot=αNStot=mSfD2fD1NStot, the total starch degradation rate is

dΔStotdt=μtot=fD1AdβkSTkSV0ΔStotΔt.

According to this equation, starch contents are completely depleted by the time of expected dawn, and for a degradation rate μtot=ΔStotΔt with normalization relation fD1AdβkSTkSV0=1, starch contents decrease linearly with time.

Calculation of the starch content during the degradation process

In this section we describe how starch contents during the degradation process were calculated using the equations previously derived.

Model 1

Request a detailed protocol

The starch content as function of time is obtained by solving the equation dΔStotdt=μtot. For Model 1, from the previous section, we have:

μtot=mSfD2AdkSn[SC]1+kST1[TC].

As we showed, [SC] can, at all times, be approximated by [SC]=NStotnV0=ΔStotαnV0=fD1mSfD2nV0ΔStot; therefore, we find:

dΔStotdt=fD1AdkSV0ΔStot1+kST1[TC].

If we put γ=fD1AdkskST1V0β, we then have

dΔStotdt=γβkST1ΔStot1+kST1[TC].

To convert this equation into one for starch content ρS in mg g−1 FW, we use ρS = (ΔStot/n) ν, (ΔStot/n) being the amount of starch in a single granule in mg and ν the number of starch granules g−1 FW. This leads to

dρSdt=γρS(βkST1)1+([TC]/β).

Given Equation 1 for [TC], the above differential equation can be exactly solved analytically. If we define the functions:

f(1)(x,y)=((βkST1)1+tdaytrtry(βkST1)1+tdaytrtrx)γtrtdaytr,f(2)(x,y)=((βkST1)1+tdayy(βkST1)1+tdayx)γ,

and t* is the time at which degradation starts, with ρS(t*) = ρ0 the corresponding starch content, then, for 0ttr, the solution is:

ρS(t)={ρ0f(1)(t,t)tttrρ0f(1)(t,tr)f(2)(tr,t)trttday,

while, for trttday, the solution becomes:

ρS(t)={ρ0f(2)(t,t)tttdayρ0f(2)(t,tday)f(1)(0,ttday)tdayttday+tr.

According to these equations, for (βkST1)−1 sufficiently small, starch contents are almost completely depleted by the time of expected dawn, and for γ = 1, starch contents decrease linearly with time for most of the dark period, except around the time of expected dawn, when [TC] is being reset (see Figure 2B).

Model 2

Request a detailed protocol

Similarly, for Model 2, we have:

dρSdt=γρS[TC]β,

where γ=fD1AdβkSTkSV0 and [Tc] is defined by Equation 2. Then, given the three functions:

f(1)(x,y)=exp{γ2(yx)(2τ2y+xτ12)},f(2)(x,y)=(tdayytdayx)γ,f(3)(x,y)=exp{γ2(yx)(2τ3y+xτ12)},

with

1τ12=tr2tr1(tday+tr1tr2)(tdaytr1)(tdaytr2),1τ2=1(tday+tr1tr2)(tdaytr2tdaytr1+tr1tdaytr2),1τ3=1(tday+tr1tr2)(tday+tr1tdaytr2tr2tdaytr1),

if the initial condition is ρS(t*) = ρ0, and 0ttr1, the exact solution is:

ρS(t)={ρ0f(1)(t,t)tttr1ρ0f(1)(t,tr1)f(2)(tr1,t)tr1<ttr2ρ0f(1)(t,tr1)f(2)(tr1,tr2)f(3)(tr2,t)tr2<ttday.

For tr1ttr2, we have:

ρS(t)={ρ0f(2)(t,t)tttr2ρ0f(2)(t,tr2)f(3)(tr2,t)tr2ttdayρ0f(2)(t,tr2)f(3)(tr2,tday)f(1)(0,ttday)tdayttr1+tday,

and for tr2ttday:

ρS(t)={ρ0f(3)(t,t)tttdayρ0f(3)(t,tday)f(1)(0,ttday)tdayttr1+tdayρ0f(3)(t,tday)f(1)(0,tr1)f(2)(tr1,ttday)tr1+tdayttr2+tday.

From the previous equations, it is easily seen that for tr2 close to tday starch contents are almost completely depleted by the time of expected dawn, and for γ = 1, starch contents decrease linearly with time for tr1ttr2, although this linearity no longer holds around the time of expected dawn, when [TC] is being reset (see Figure 2C).

Parameters

For Model 1, the values of the following parameters are required: γ, (βkST1)−1, ρ0, tr; and for Model 2: γ, ρ0, tr1 and tr2.

In order to fit the data from the night-time light period experiments, we also considered the possibility that the T dynamics given by Equations 1 and 2 (see also Figure 2B,C) can be phase shifted by a time t0. The addition of t0 as an extra fitting parameter for the night-time light period experiments was justified by our data on transcripts of the clock gene LHY (Figure 3—figure supplement 1), showing that the night-time light period may induce a phase shift in the expression of some genes.

The phase shift of the T dynamics determines a shift in the time of expected dawn, that, in Model 1, is defined as the time at which [TC] falls to zero (see Equation 1 and Figure 2B) and which in Model 2 is the time at which [TC] would diverge without a reset (see Equation 2 and Figure 2C). Therefore, if the T dynamics are phase shifted by t0, in our models the starch degradation rates are adjusted in such a way so as to deplete the starch reserves at (24 + t0) hr after the previous dawn, instead of the normal 24 hr. We found this phenotype in cca1/lhy plants, which run out of starch earlier than 24 hr after the previous dawn. Hence, to reproduce the phenotype of cca1/lhy in our models, t0 was also used as a fitting parameter. A more extensive discussion about cca1/lhy can be found below in the ‘Mutant phenotypes’ section.

A full list of the system parameters is shown in Tables 1 and 2.

Table 1

Full list of the system parameters for Model 1

https://doi.org/10.7554/eLife.00669.012
SymbolDefinition
ρ0Starch content at the beginning of the dark period.
γNormalization variable.
(βkST1)−1kST1 is the ratio of the reaction parameter associated with reaction rST1 with the backward rate of reaction rS. β is the proportionality constant between [TC] and Δt.
trTime at which [TC] levels finish being reset at the beginning of the day (see Figure 2B).
t0Phase shifting parameter of the [TC] dynamics given by Equation 1. The next dawn is expected to come (24 + t0) hr after the previous one.
Table 2

Full list of the system parameters for Model 2

https://doi.org/10.7554/eLife.00669.013
SymbolDefinition
ρ0Starch content at the beginning of the dark period.
γNormalization variable.
tr1Time at which [TC] levels finish being reset at the beginning of the day (see Figure 2C).
tr2Time at which [TC] levels start being reset at the end of the day (see Figure 2C).
t0Phase shifting parameter of the [TC] dynamics given by Equation 2. The next dawn is expected to come (24 + t0) hr after the previous one.

Data fitting

Request a detailed protocol

The best fit to a given dataset was found by minimizing the function:

L(θ)=i=1N(ρStheory(ti,θ)ρSexp(ti)ρSexp(ti))2,

where ρSexp(ti) are the N mean values of starch contents measured at times ti. ρStheory(ti,θ) is the starch value predicted by the model at time ti with parameter values θ. The set of parameters θ that minimizes L corresponds to the maximum-likelihood estimates of the parameters under the assumption that the experimental measurements are normally distributed around the theoretical values with a constant relative error.

Note that, in the fits of the data from the cca1/lhy mutant plants for the early and the normal night, we did not consider the two data points closest in time to 24 hr from the datasets, as they are characterized by very low values of starch content and, therefore, are likely to be affected by higher relative errors compared to the other data points. For the same reason, for the linear fits shown in Figure 4—figure supplement 1A we did not consider the data point at t = 24 hr in the normal night datasets and the data points at t = 22 hr and t = 24 hr in the early night datasets.

A simulated annealing algorithm was used for the minimization of L(θ). The parameter values of the best fits of the models are given in Tables 3–6, along with the ranges in which the parameters were allowed to vary.

Table 3

The values of the parameters of the Model 1 best fits shown in Figure 2

https://doi.org/10.7554/eLife.00669.014
Model 1
WT early night (panel F)WT normal night (panel F)WT late night
(panel F)
cca1/lhy early night (panel H)cca1/lhy normal night (panel H)WT low light level (panel J)WT normal light level (panel J)
ρ0 (mg g−1 FW) within 10% of the measured value8.511.011.74.25.13.66.1
γ (0.7–3.0)1.81.81.91.21.21.51.3
(βkST1)−1 (1.5–5.0) hr5.05.05.01.61.72.11.5
tr (9.0–12.0) hr9.0Any value in the specified rangeAny value in the specified range11.711.7Any value in the specified rangeAny value in the specified range
t0 (−5.0–5.0) hr for the cca1/lhy data, t0 = 0 for WT0.00.00.0−4.2−2.50.00.0
  1. With each parameter the range used in the best fit search is indicated.

  2. WT: wild-type plants.

Table 4

The values of the parameters of the Model 2 best fits shown in Figure 2

https://doi.org/10.7554/eLife.00669.015
Model 2
WT early night (panel G)WT normal night (panel G)WT late night (panel G)cca1/lhy early night (panel I)cca1/lhy normal night (panel I)WT low light level (panel K)WT normal light level (panel K)
ρ0 (mg g−1 FW) within 10% of the measured value8.311.012.14.25.13.55.8
γ (0.7–3.0)1.11.01.11.31.21.20.9
tr1 (9.0–12.0) hr11.012.012.011.011.410.510.2
tr2 (20.0–23.0) hr20.120.020.021.821.421.422.5
t0 (−5.0–5.0) hr for the cca1/lhy data, t0 = 0 for WT0.00.00.0−2.4−1.30.00.0
  1. With each parameter the range used in the best fit search is indicated.

  2. WT: wild-type plants.

Table 5

The values of the parameters of the Model 1 best fits shown in Figure 3

https://doi.org/10.7554/eLife.00669.016
Model 1
Panel APanel BPanel CPanel D
ρ0 (mg g−1 FW) within 10% of the measured value5.49.07.74.2
γ (0.7–3.0)2.41.41.11.0
(βkST1)−1 (1.5–5.0) hr5.02.62.12.0
tr (9.0–12.0) hr9.09.010.312.0
t0 (−5.0—5.0) hr4.33.32.4−0.5
  1. With each parameter the range used in the best fit search is indicated.

Table 6

The values of the parameters of the Model 2 best fits shown in Figure 3—figure supplement 2

https://doi.org/10.7554/eLife.00669.017
Model 2
Panel APanel BPanel CPanel D
ρ0 (mg g−1 FW) within 10% of the measured value5.59.17.74.2
γ (0.7–3.0)1.51.30.70.7
tr1 (9.0–12.0) hr12.09.09.09.0
tr2 (20.0–23.0) hr20.020.922.122.1
t0 (−5.0–5.0) hr5.05.02.40.0
  1. With each parameter the range used in the best fit search is indicated.

Recalculation of the starch degradation rate after the night-time light period

Request a detailed protocol

By combining our results from the night-time light period experiments shown in Figure 3A–C, we can show that the evidence in favor of a re-calculation of the starch degradation rate after the night-time light period is statistically significant. In the experiments shown in Figure 3A–C, plants were subjected to 5 hr of light in the middle of the night (between 14 hr and 19 hr after the previous dawn). The starch content at each time point was measured by averaging the starch content of n = 10–12 individual rosettes and the standard error of the mean was also calculated.

The starch degradation rate before the night-time light period in the i-th experiment was measured as:

μBi=ρSi(12hr)ρSi(14hr)14hr12hr,

that is, as the difference between the starch content at 12 hr and 14 hr divided by the time interval.

Similarly, the starch degradation rate after the night-time light period for the i-th experiment was:

μAi=ρSi(19hr)ρSi(21hr)21hr19hr,

These rates were averaged over the three experiments, and the difference between the two averages was calculated:

μAμB=1Ni=1N(μAiμBi),

where N = 3 was the number of experiments. The error on this quantity was estimated by propagating standard errors of the mean. We found μA − μB = 0.43 ± 0.14 (mg g−1 FW)/hr which is significantly greater than 0 (one-tailed Welch’s t-test, p = 1.6 × 10−3).

We also calculated the average difference between μAi and

μnormali=ρSi(12hr)12hr,

where μnormali is the rate that would have been expected over a normal 12-hr night. We found that this average is

μAμnormal=1Ni=1N(μAiμnormali)=0.48±0.11(mgg1FW)/hr,

where the error was again estimated by propagating the standard errors of the mean.

This value is again significantly greater than 0 (one-tailed Welch’s t-test, p = 3.8 × 10−5).

These results provide strong evidence against the hypothesis that a fixed degradation rate is set only once at the first onset of darkness, and is instead compatible with our prediction that the rate is continuously re-computed throughout the night.

Mutant phenotypes

Request a detailed protocol

In the following, we show how the models can explain the mutant phenotypes discussed in the main text (see Figures 1B and 4 and Figure 4—figure supplement 1).

cca1/lhy

Request a detailed protocol

LHY and CCA1 are central components of the clock. The cca1/lhy mutant has a free-running period of significantly <24 hr under continuous light.

This mutant is characterized by too high a rate of starch degradation (see Figure 1B). Indeed, the starch reserve is exhausted around 21–22 hr after the previous dawn, instead of 24 hr as in the wild type. Interestingly, if these mutant plants are given an early night, the starch degradation rate is adjusted such that all the starch is again degraded around 21–22 hr after previous dawn (see Figure 1B).

The models we discussed previously can straightforwardly explain such a phenotype by assuming that in this mutant, the time of expected dawn is shifted to a time t < 24 hr after the previous dawn. There are different perturbations of the T dynamics that can produce this effect in our models, and still reproduce the data from cca1/lhy equally well. For instance, for cca1/lhy in Model 1, [TC] after being reset, could decrease more steeply than in the wild-type, and drop to zero at a time t < 24 hr after the previous dawn. Then [TC] could remain at low levels, before rising again around 24 hr. Another perturbation appropriate for both Models 1 and 2, would be to phase shift the T dynamics given by Equations 1 and 2, respectively by a time t0 < 0 hr, in such a way that the time of expected dawn becomes (24 + t0) hr < 24 hr, as discussed above in the ‘Parameters’ section. For the sake of simplicity, we fitted the cca1/lhy data by assuming that the latter perturbation takes place, and accordingly we used t0 as a fitting parameter.

lsf1, sex4, bam3, bam4, isa3

Request a detailed protocol

lsf1 and sex4 mutant plants were kept in the dark for 132 hr after the end of the previous light period, then exposed to normal 12-hr light/12-hr dark cycles and the starch content measured. As Figure 4A shows, during the days following the prolonged dark period the mutant plants failed to degrade their entire starch reserve at night. Instead, the total starch content degraded by the end of each night, expressed as a percentage of the starch content at the end of the respective preceding light period, was approximately constant and much lower than that of wild-type plants (Figure 4B).

Interestingly, the same mutants could also adjust their starch degradation rate in response to an unexpected early night. In order to quantitatively verify this observation, we performed linear fits of the data from the unexpected early and normal nights:

ρSearly(t)=μearly(ttearly)+ρS,0early,ρSnormal(t)=μnormal(ttnormal)+ρS,0normal,

where (μearly,ρS,0early) and (μnormal,ρS,0normal) are the fitting parameters, with tearly = 8 hr, tnormal = 12 hr the times of onset of darkness respectively for the unexpected early and normal night. From these fits, we calculated the ratio R between the starch degradation rates in the normal and the unexpected early night, normalized by the respective starch content at the time of onset of darkness:

R=μnormal/ρS,0normalμearly/ρS,0early.

If the degradation rate is adjusted as happens in wild-type plants, during the unexpected early night relative to a normal night, the same fraction f of the initial starch content should be degraded by the time of expected dawn. Therefore, μnormalfρS,0normaltdaytnormal and μearlyfρS,0earlytdaytearly; hence, Rtdaytearlytdaytnormal=16121.33. As Figure 4F shows, the values of R found for lsf1 and sex4 are compatible with 1.33 (dashed line), therefore supporting the hypothesis that these mutants are able to appropriately adjust their starch degradation rate in response to an unexpected early night. We also found the same type of rate adjustment in bam3, bam4 and isa3 mutants (see Figure 4F and Figure 4—figure supplement 1), which also failed to exhaust their starch reserves by the end of the night.

The models make a precise prediction about how such a phenotype can be produced. We will focus on Model 1, but analogous arguments also hold true for Model 2. We showed that for Model 1, in order to obtain a linearly decreasing starch content with time, and with the appropriate degradation rate to ensure complete starch depletion at the time of expected dawn, two normalization conditions must hold:

fD1AdkSkST1V0β=γ=1,mSfD2αfD1=χ=1.

The first condition ensures that the number of S molecules, and therefore the starch content, decreases linearly with time during the night. The second condition is needed to keep the number of S molecules and the amount of starch strictly proportional, so that both these quantities are fully degraded by the end of the night. We now assume that the second of the two conditions is not valid, and, in particular, that χ < 1. In this case, the starch degradation proceeds at a slower rate, breaking the proportionality between the number of S molecules and the amount of starch. By using Equations 4 and 8, and the initial condition at the start of the dark period, ΔStot=αNStot, we find that a fraction of starch approximately equal to χ only is degraded by the time of expected dawn. Such a phenotype is clearly compatible with the disruptions caused by the above mutations. Indeed by assuming that the second condition above is not valid, our models can be well fit to the full night-time starch profiles in all the above mutants (data not shown).

One way to perturb only the second of the two normalization conditions shown above (and not the first), would be to alter fD2, which is the starch degradation reaction parameter. This reasoning fits well with the known functions of the above genes, all of which play roles, directly or indirectly, in starch degradation.

pwd

Request a detailed protocol

The pwd mutant lacks the glucan water dikinase responsible for the addition of the phosphate groups to the 3-position of glucose moieties in starch. This mutation generated an approximately linear decrease of starch content with time during the night, but with a rate which was too low to ensure the complete utilization of the starch reserve by the time of expected dawn. Yet, as opposed to the previously discussed mutations, this mutant did not have the ability to adjust the starch degradation rate in response to an unexpected early night (see Figure 4E,F). Indeed, for this mutant we found that the ratio between the degradation rates (normalized by the respective end-of-light period starch content) during the normal and the early night was R = 1.10 ± 0.10, which is significantly less than the wild-type value (R = 16/12, one-tailed Z-test, p = 0.01).

This finding indicates that PWD is potentially a key node where information about starch content and expected time to dawn from the circadian clock are integrated to control starch degradation dynamics. This makes PWD an obvious target for future experiments that aim to elucidate the identities of the S and T molecules.

Other examples of models implementing arithmetic division

Request a detailed protocol

The models detailed above use two distinctly different methods of implementing the division operation. In one case, exemplified by Model 1, the S and T molecule concentrations are divided, with the T molecule concentration tracking the time to expected dawn. In Model 2, the S and T concentrations are multiplied, with the T molecule concentration tracking the reciprocal of the time to expected dawn. The precise implementation of these two methods of performing arithmetic division could, however, vary, with slightly different reaction schemes but the same underlying principles. To illustrate this point, we briefly outline other models related to Model 1 above, where the TC molecule again has the behavior given in Equation 1. We first consider:

rS:SCSrT:TCTrT2:T+TT2rST:S+TSTrST2:S+T2ST2rSCT2:ST2SC+T2rD1:STTrD2:ST+StarchST.

Here both S and T molecules can directly and reversibly associate to the granule surface (reactions rS and rT). Once bound to the surface, S and T can form a complex ST (reaction rST) that permits S molecule and starch degradation (through reactions rD1 and rD2 respectively). The T molecule can also hinder starch degradation by forming dimers (reaction rT2) that are able to associate with S molecules (reaction rST2) and induce them to detach from the granule surface (reaction rSCT2). By a similar analysis to that carried out for Model 1 and Model 2, and for similar reasons, it can be shown that this model can also implement arithmetic division between the S and T molecule concentrations.

More complex models can also easily accommodate additional molecules, which, for instance, could recruit S and T molecules to the granule surface and be part of the starch degradation apparatus. For example, in the set of reactions:

rM:MCMrS:SC+MSMrT:TC+MTMrSM:SMSC+MCrTM:TMTC+MCrD1:SMMCrD2:SM+StarchSM,

M molecules which reversibly bind to the granule surface (reaction rM), recruit the S and T molecules from the surrounding compartment (reactions rS and rT), to form the SM and TM complexes, which can then dissociate from the surface through the reactions rSM and rTM. The SM complex can permit S molecule and starch degradation (reactions rD1 and rD2). Since the T molecule hinders starch degradation by binding to M molecules and preventing them from binding to S, it can be shown that this model can also implement arithmetic division between the S and T molecule concentrations.

Validation of the use of a night-time light period

Request a detailed protocol

Exposure of Arabidopsis plants to light during the normal night could potentially lead to the early, light-induced expression of clock genes, and under some circumstances might re-entrain the clock. Such re-entrainment might affect the level of the T molecule. To discover whether our experimental conditions gave rise to such problems, we investigated the behavior of transcript levels of LHY, a central clock gene, by qPCR analysis. For the experiment shown in Figure 3B, LHY transcript levels were slightly elevated during the night-time light period, then peaked about 2 hr later than in plants that were not exposed to the light period (Figure 3—figure supplement 1). This result indicates that the night-time light period has only minor effects on the clock, which is not re-entrained: if this were the case the peak of LHY expression observed with the night-time light period would not be expected to occur at around 24 hr after the previous dawn.

References

    1. Smith AM
    (2012)
    Starch in the Arabidopsis plant
    Starch/Stärke 64:421–434.

Article and author information

Author details

  1. Antonio Scialdone

    Department of Computational and Systems Biology, John Innes Centre, Norwich, United Kingdom
    Contribution
    AS, Constructed and analyzed the models, statistically analyzed the results, wrote the paper
    Contributed equally with
    Sam T Mugford
    Competing interests
    The authors declare that no competing interests exist.
  2. Sam T Mugford

    Department of Metabolic Biology, John Innes Centre, Norwich, United Kingdom
    Contribution
    STM, Performed the Arabidopsis experiments, wrote the paper
    Contributed equally with
    Antonio Scialdone
    Competing interests
    The authors declare that no competing interests exist.
  3. Doreen Feike

    Department of Metabolic Biology, John Innes Centre, Norwich, United Kingdom
    Contribution
    DF, Performed the Arabidopsis experiments
    Competing interests
    The authors declare that no competing interests exist.
  4. Alastair Skeffington

    Department of Metabolic Biology, John Innes Centre, Norwich, United Kingdom
    Contribution
    AS, Performed the Arabidopsis experiments
    Competing interests
    The authors declare that no competing interests exist.
  5. Philippa Borrill

    1. Department of Metabolic Biology, John Innes Centre, Norwich, United Kingdom
    2. Department of Crop Genetics, John Innes Centre, Norwich, United Kingdom
    Contribution
    PB, Performed the Brachypodium experiments
    Competing interests
    The authors declare that no competing interests exist.
  6. Alexander Graf

    Department of Metabolic Biology, John Innes Centre, Norwich, United Kingdom
    Present address
    Laboratory of Plant Biotechnology, ETH Zurich, Zurich, Switzerland
    Contribution
    AG, Performed the Arabidopsis experiments
    Competing interests
    The authors declare that no competing interests exist.
  7. Alison M Smith

    Department of Metabolic Biology, John Innes Centre, Norwich, United Kingdom
    Contribution
    AMS, Conceived the study, wrote the paper
    For correspondence
    alison.smith@jic.ac.uk
    Competing interests
    The authors declare that no competing interests exist.
  8. Martin Howard

    Department of Computational and Systems Biology, John Innes Centre, Norwich, United Kingdom
    Contribution
    MH, Conceived the study, constructed and analyzed the models, wrote the paper
    For correspondence
    martin.howard@jic.ac.uk
    Competing interests
    The authors declare that no competing interests exist.

Funding

Biotechnology and Biological Sciences Research Council (UK) (BB/J004561/1)

  • Alison M Smith
  • Martin Howard

European Union FP7 (Collaborative project grant 245143 (TiMet))

  • Alison M Smith

The John Innes Foundation

  • Alison M Smith

The Gatsby Charitable Foundation

  • Alison M Smith

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Mark Stitt and members of the Howard group for discussions.

Reviewing Editor

  1. Detlef Weigel, Max Planck Institute for Developmental Biology, Germany

Version history

  1. Received: February 25, 2013
  2. Accepted: May 31, 2013
  3. Version of Record published: June 25, 2013 (version 1)

Copyright

© 2013, Scialdone et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 6,966
    Page views
  • 720
    Downloads
  • 121
    Citations

Article citation count generated by polling the highest count across the following sources: Scopus, Crossref, PubMed Central.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Antonio Scialdone
  2. Sam T Mugford
  3. Doreen Feike
  4. Alastair Skeffington
  5. Philippa Borrill
  6. Alexander Graf
  7. Alison M Smith
  8. Martin Howard
(2013)
Arabidopsis plants perform arithmetic division to prevent starvation at night
eLife 2:e00669.
https://doi.org/10.7554/eLife.00669

Share this article

https://doi.org/10.7554/eLife.00669

Further reading

    1. Plant Biology
    Ivan Kulich, Julia Schmid ... Jiří Friml
    Research Article

    Root gravitropic bending represents a fundamental aspect of terrestrial plant physiology. Gravity is perceived by sedimentation of starch-rich plastids (statoliths) to the bottom of the central root cap cells. Following gravity perception, intercellular auxin transport is redirected downwards leading to an asymmetric auxin accumulation at the lower root side causing inhibition of cell expansion, ultimately resulting in downwards bending. How gravity-induced statoliths repositioning is translated into asymmetric auxin distribution remains unclear despite PIN auxin efflux carriers and the Negative Gravitropic Response of roots (NGR) proteins polarize along statolith sedimentation, thus providing a plausible mechanism for auxin flow redirection. In this study, using a functional NGR1-GFP construct, we visualized the NGR1 localization on the statolith surface and plasma membrane (PM) domains in close proximity to the statoliths, correlating with their movements. We determined that NGR1 binding to these PM domains is indispensable for NGR1 functionality and relies on cysteine acylation and adjacent polybasic regions as well as on lipid and sterol PM composition. Detailed timing of the early events following graviperception suggested that both NGR1 repolarization and initial auxin asymmetry precede the visible PIN3 polarization. This discrepancy motivated us to unveil a rapid, NGR-dependent translocation of PIN-activating AGCVIII kinase D6PK towards lower PMs of gravity-perceiving cells, thus providing an attractive model for rapid redirection of auxin fluxes following gravistimulation.

    1. Plant Biology
    Daniel S Yu, Megan A Outram ... Simon J Williams
    Research Article

    Plant pathogens secrete proteins, known as effectors, that function in the apoplast or inside plant cells to promote virulence. Effector recognition by cell-surface or cytosolic receptors results in the activation of defence pathways and plant immunity. Despite their importance, our general understanding of fungal effector function and recognition by immunity receptors remains poor. One complication often associated with effectors is their high sequence diversity and lack of identifiable sequence motifs precluding prediction of structure or function. In recent years, several studies have demonstrated that fungal effectors can be grouped into structural classes, despite significant sequence variation and existence across taxonomic groups. Using protein X-ray crystallography, we identify a new structural class of effectors hidden within the secreted in xylem (SIX) effectors from Fusarium oxysporum f. sp. lycopersici (Fol). The recognised effectors Avr1 (SIX4) and Avr3 (SIX1) represent the founding members of the Fol dual-domain (FOLD) effector class, with members containing two distinct domains. Using AlphaFold2, we predicted the full SIX effector repertoire of Fol and show that SIX6 and SIX13 are also FOLD effectors, which we validated experimentally for SIX6. Based on structural prediction and comparisons, we show that FOLD effectors are present within three divisions of fungi and are expanded in pathogens and symbionts. Further structural comparisons demonstrate that Fol secretes effectors that adopt a limited number of structural folds during infection of tomato. This analysis also revealed a structural relationship between transcriptionally co-regulated effector pairs. We make use of the Avr1 structure to understand its recognition by the I receptor, which leads to disease resistance in tomato. This study represents an important advance in our understanding of Fol-tomato, and by extension plant–fungal interactions, which will assist in the development of novel control and engineering strategies to combat plant pathogens.