Complementary approaches to understanding the plant circadian clock

Circadian clocks are oscillatory genetic networks that help organisms adapt to the 24-hour day/night cycle. The clock of the green alga Ostreococcus tauri is the simplest plant clock discovered so far. Its many advantages as an experimental system facilitate the testing of computational predictions. We present a model of the Ostreococcus clock in the stochastic process algebra Bio-PEPA and exploit its mapping to different analysis techniques, such as ordinary differential equations, stochastic simulation algorithms and model-checking. The small number of molecules reported for this system tests the limits of the continuous approximation underlying differential equations. We investigate the difference between continuous-deterministic and discrete-stochastic approaches. Stochastic simulation and model-checking allow us to formulate new hypotheses on the system behaviour, such as the presence of self-sustained oscillations in single cells under constant light conditions. We investigate how to model the timing of dawn and dusk in the context of model-checking, which we use to compute how the probability distributions of key biochemical species change over time. These show that the relative variation in expression level is smallest at the time of peak expression, making peak time an optimal experimental phase marker. Building on these analyses, we use approaches from evolutionary systems biology to investigate how changes in the rate of mRNA degradation impacts the phase of a key protein likely to affect fitness. We explore how robust this circadian clock is towards such potential mutational changes in its underlying biochemistry. Our work shows that multiple approaches lead to a more complete understanding of the clock.

💡 Research Summary

The paper presents a comprehensive, multi‑method study of the circadian clock in the green alga Ostreococcus tauri, which is currently the simplest known plant‑like clock. The authors first encode the core transcription‑translation feedback loop of this organism in Bio‑PEPA, a stochastic process algebra that allows a single formal model to be automatically translated into three distinct analytical frameworks: ordinary differential equations (ODEs), stochastic simulation algorithms (SSA), and probabilistic model‑checking.

Using the ODE translation, the authors obtain a deterministic description of the system that reproduces the overall period and amplitude observed experimentally. However, because Ostreococcus cells contain only a few hundred copies of each mRNA and protein, the continuous approximation underlying ODEs is stretched to its limits. The deterministic model fails to capture the large intrinsic noise and phase variability that arise when molecule numbers are low, especially during the trough phases of the oscillation.

To address this shortcoming, the authors apply Gillespie’s SSA to the same Bio‑PEPA specification. The stochastic simulations explicitly model each reaction event, thereby preserving the discrete nature of the system. Remarkably, the simulations predict that even under constant light (LL) conditions—traditionally considered arrhythmic for many organisms—single Ostreococcus cells can maintain self‑sustained oscillations. This hypothesis, which emerges only when stochasticity is taken into account, suggests that population‑averaged measurements may mask persistent rhythmicity at the single‑cell level.

The third analytical route employs the PRISM model checker. Here, the authors formalise dawn and dusk as timed logical events and compute the time‑dependent probability distributions of key species such as CCA1, TOC1, and their mRNA transcripts. Model‑checking results reveal that the coefficient of variation of these species is minimal at the time of peak expression, indicating that peak time is the most reliable experimental phase marker. Moreover, the probabilistic analysis quantifies how the distributions shift in response to light cues, providing a detailed picture of the clock’s responsiveness that is difficult to obtain from deterministic or purely stochastic simulations alone.

Building on these quantitative insights, the authors explore evolutionary robustness by varying the degradation rate of clock‑related mRNA (k_deg). Through systematic parameter sweeps and sensitivity analyses, they show that within a biologically plausible range, changes in k_deg cause only modest shifts in the phase of a downstream protein that is likely to affect fitness. This phase stability implies that the Ostreococcus clock is buffered against mutational perturbations in its biochemical parameters, a property that could have been selected for during evolution. When k_deg exceeds a critical threshold, however, the phase shifts dramatically, potentially compromising the organism’s adaptation to the day‑night cycle.

Overall, the study demonstrates that integrating continuous‑deterministic, discrete‑stochastic, and formal verification approaches yields a richer, more nuanced understanding of a biological oscillator than any single method could provide. The deterministic ODE analysis offers a quick overview of average dynamics, the stochastic SSA uncovers noise‑driven phenomena and single‑cell behaviours, and model‑checking supplies rigorous, quantitative statements about temporal probability distributions and system robustness. By applying this triple‑pronged strategy to the Ostreococcus clock, the authors not only validate the utility of Bio‑PEPA as a unifying modelling language but also generate new biological hypotheses—such as persistent LL oscillations and optimal phase markers—that can be tested experimentally. Their work thus sets a methodological benchmark for future studies of circadian systems and other complex biochemical networks.