Modelling the widespread effects of TOC1 signalling on the plant circadian clock and its outputs
© Pokhilko et al.; licensee BioMed Central Ltd. 2013
Received: 10 October 2012
Accepted: 8 March 2013
Published: 19 March 2013
24-hour biological clocks are intimately connected to the cellular signalling network, which complicates the analysis of clock mechanisms. The transcriptional regulator TOC1 (TIMING OF CAB EXPRESSION 1) is a founding component of the gene circuit in the plant circadian clock. Recent results show that TOC1 suppresses transcription of multiple target genes within the clock circuit, far beyond its previously-described regulation of the morning transcription factors LHY (LATE ELONGATED HYPOCOTYL) and CCA1 (CIRCADIAN CLOCK ASSOCIATED 1). It is unclear how this pervasive effect of TOC1 affects the dynamics of the clock and its outputs. TOC1 also appears to function in a nested feedback loop that includes signalling by the plant hormone Abscisic Acid (ABA), which is upregulated by abiotic stresses, such as drought. ABA treatments both alter TOC1 levels and affect the clock’s timing behaviour. Conversely, the clock rhythmically modulates physiological processes induced by ABA, such as the closing of stomata in the leaf epidermis. In order to understand the dynamics of the clock and its outputs under changing environmental conditions, the reciprocal interactions between the clock and other signalling pathways must be integrated.
We extended the mathematical model of the plant clock gene circuit by incorporating the repression of multiple clock genes by TOC1, observed experimentally. The revised model more accurately matches the data on the clock’s molecular profiles and timing behaviour, explaining the clock’s responses in TOC1 over-expression and toc1 mutant plants. A simplified representation of ABA signalling allowed us to investigate the interactions of ABA and circadian pathways. Increased ABA levels lengthen the free-running period of the clock, consistent with the experimental data. Adding stomatal closure to the model, as a key ABA- and clock-regulated downstream process allowed to describe TOC1 effects on the rhythmic gating of stomatal closure.
The integrated model of the circadian clock circuit and ABA-regulated environmental sensing allowed us to explain multiple experimental observations on the timing and stomatal responses to genetic and environmental perturbations. These results crystallise a new role of TOC1 as an environmental sensor, which both affects the pace of the central oscillator and modulates the kinetics of downstream processes.
KeywordsCircadian rhythms Biological clocks Gene regulatory networks Mathematical model Systems biology
Circadian clocks allow most eukaryotes and some prokaryotes to anticipate the environmental day/night cycle, through rhythmic modulation of multiple physiological processes . Circadian clocks are characterised by free running rhythms with a period of ~24 h even in the absence of any environmental cues, such as in constant light conditions. Their timing is synchronised with the environmental day-night cycle principally by responses to light and temperature, but clock circuits also respond to many additional stimuli. The signalling pathways from such stimuli are often found to be rhythmically controlled by the clock, forming nested feedback loops that modify the circadian oscillator. The gene circuit of the higher plant clock consists of interlocked morning and evening loops operating through negative feedback mechanisms . TOC1 is an important component of the plant clock since its mutation or overexpression dramatically change the properties of the central oscillator .
Although the importance of TOC1 in the plant circadian clock has long being recognized, the precise role of TOC1 in the clock has been a matter of debate [2–7]. Based on indirect genetic evidence it was first proposed that TOC1 activated the expression of the key transcription factors LHY and CCA1, although several observations on clock mutants remained paradoxical. Recent experimental and computational data [2, 3, 7] have demonstrated that TOC1 functions as a repressor rather than an activator. The studies also show that TOC1 repression is not limited to CCA1 and LHY but to nearly all of the genes at the core of the oscillator . The observed reduction in the expression of multiple clock genes by up-regulation of TOC1 could not be described by our previous model, where only LHY and CCA1 were affected by TOC1 [2, 3]. Here, we extend our previous model by including the repression of multiple clock genes by TOC1 and explore TOC1 effect on the clock. In particular, we analysed clock dynamics in various TOC1-misexpressed lines and addressed the long-standing paradox of simultaneous reduction of LHY and CCA1 amplitudes in both toc1 mutant and TOC1-overexpressed (TOC1-ox) plants [7, 8].
The recently discovered pervasive negative effect of TOC1 on gene expression [3, 7], the overlap between microarrays of the TOC1- and ABA-regulated genes and the gated induction of TOC1 by the stress hormone abscisic acid ABA  suggests that TOC1 might integrate environmental and circadian information to regulate downstream physiology at specific time of the day . Here we combined the clock model with key processes regulated by ABA and explored mutual interactions between clock and ABA pathways. The model describes the regulation of stomatal aperture, one of the most important outputs of ABA signalling, which is gated by the clock . Being a sensor of water availability in plants, ABA level increases under dry conditions, which accelerates closing of stomata in the afternoon [9–11]. In agreement with the microarray studies, TOC1 mis-expressing plants were shown to have altered ABA-dependent stomata closure and affected responses to drought conditions . To model the observed gating of stomatal opening by the clock, we included the interaction of ABA with the ABAR protein (also known as the H subunit of Mg-chelatase, CHLH or GUN5) – a key component of the ABA signalling pathway, arguably postulated as ABA receptor . Besides being one of the downstream targets of TOC1, ABAR was shown to be a key component in the gating of ABA signalling by the clock [9, 13]. We used the model to explore the dynamics of the opening of stomatal pore under various genetic and environmental perturbations, such as mutations of clock genes and by changing the light and water conditions, and showed that the model provides good match to existing experimental data. Our results show how the dual roles of TOC1 are integrated, combining timing and environmental information to modulate downstream physiology.
Results and discussion
Extending the clock model by inhibition of gene expression by TOC1
The latest model of the Arabidopsis oscillator  was extended based on our recent data on the inhibition of expression of multiple clock genes by TOC1 and induction of TOC1 through ABA signalling. The model is described by a system of 32 ODEs. The basic structure of the plant clock was kept from our previous P2011 model . Briefly, the key components of the clock are the LHY and CCA1 morning genes and the evening complex (EC) genes. LHY and CCA1 are similarly expressed around dawn and act together to regulate multiple targets. They are described by single variable LHY/CCA1 as before. EC (EVENING COMPLEX) genes LUX, ELF3 (EARLY FLOWERING 3) and ELF4 are expressed around dusk and form the EC protein complex, which suppresses expression of multiple genes at night (Figure 1; ). LHY/CCA1 expression is suppressed by the PRR proteins in a morning loop, while the EC is negatively regulated by the ubiquitin E3 ligase COP1 (CONSTITUTIVE PHOTOMORPHOGENIC 1), which targets ELF3 protein to degradation by proteasome and by GI protein, in the evening loop. The morning and evening loops are further interlocked through the suppression of evening gene expression by LHY/CCA1 and suppression of PRR9 and TOC1 (PRR1) expression by the EC. The model also includes the F box protein ZTL (ZEITLUPE), which negatively regulates the level of TOC1 protein. Light entrains the clock similarly to the P2011 model through several mechanisms, which are supported by experimental data . The main mechanisms include acute activation of LHY/CCA1, PRR9 and GI transcription (eqs. 1, 5, 27) immediately after dawn; stabilization of PRR proteins in presence of light (eqs. 6, 8, 10) and light-dependent regulation of the EC by COP1 and GI proteins (eqs. 24, 21, 22).
The TOC1-related reactions were substantially extended. Firstly, we added multiple reactions of inhibition of clock gene expression by TOC1. This includes direct inhibition of expression of LHY/CCA1 as in the P2011 model (eq. 1), and also PRR9, PRR7, PRR5, LUX, ELF4 and GI (eqs. 5, 7, 9, 18, 14, 27), which is based on our data  (Figure 1A). The details of TOC1 interactions with other regulators (e.g. LHY/CCA1, EC, protein P) to modulate gene expression are largely unknown, so we assumed that TOC1 acts as a non-competitive inhibitor of gene expression Secondly, we added physiologically relevant activation of TOC1 expression by the ABA signalling pathway (Figure 1B, ), described in section 2 below. This provided an additional level of input to the clock through ABA, which is induced by stress. Thirdly, we added regulation of stomata aperture as a clock output, which is directly related to both ABA and TOC1 signalling (Figure 1B), as presented in section 2 below.
Where stand for dimensionless concentrations of mRNA and protein, respectively. The time unit is an hour. Index “i” labels the molecular components as follows:
ABAR ABAR mRNA
AR ABA complex with ABAR protein
COP1c COP1 cytoplasmic protein
COP1d COP1 nuclear protein, day
COP1n COP1 nuclear protein, night
E3 ELF3 mRNA
E3c ELF3 cytoplasmic protein
E3n ELF3 nuclear protein
E4 ELF4 mRNA and protein
E34 nuclear protein complex ELF3-ELF4
EC nuclear protein complex ELF4-ELF3-LUX
EGc cytoplasmic protein complex ELF3-GI
EGn nuclear protein complex ELF3-GI
G GI mRNA
Gc GI cytoplasmic protein
Gn GI nuclear protein
Gn_tot total amount of GI protein in nucleus
L LHY/CCA1 mRNA and protein
Lmod LHY/CCA1 modified protein
Ltot total amount of LHY/CCA1 protein
LUX LUX mRNA and protein
P5 PRR5 mRNA and protein
P7 PRR7 mRNA and protein
P9 PRR9 mRNA and protein
PP2C PP2C active protein
SnRK2 SnRK2 active protein
T TOC1 mRNA and protein
ZG cytoplasmic protein complex GI-ZTL
ZTL ZTL protein
Symbol “s” corresponds to stomata aperture in relative units (the maximum of s is equal to 1). A0 is the total ABA level in relative units. ABA level for most of calculations was set to 1, unless it is stated otherwise (ABA0 = 1).
Where dawn and dusk are the phases of dawn and dusk (normally dawn = 0); T is the duration of twilight (we used T = 0.05 h); tanh and floor – standard functions of hyperbolic tangent and rounding operation. This representation of periodic input corresponds to the Input Signal Step Function, which now has convenient software support in SBML . Mutations of the clock genes were simulated by decreasing the rate of transcription of the corresponding gene to zero. TOC1-ox plants were modelled by adding the background activation of TOC1 translation (parameter n16 in eq. 12; n16 = 0 for wild type).
Modelling of ABA signalling and its relation with the clock
To model the mutual connections between the clock and ABA signalling, we extended the clock model by including ABA-induced reactions.
The first group of reactions of ABA signalling is related with the inactivation of protein phosphatase PP2C (protein phosphatase of 2C type) by the ABA complex (; Figure 1B). Based on the data showing relatively small diurnal changes in ABA concentration in plants [10, 16], we assumed that ABA level is constant in our model and treated it as a parameter ABA0. Next we included the potential ABA receptor ABAR, which is a necessary component of the regulation of both TOC1 expression and stomatal aperture by ABA . ABAR expression has a strong circadian pattern, peaking in the morning . Our analysis of the ABAR promoter in Arabidopsis revealed 2 CCA1 binding sites (CBS) AAATCT , in addition to the known TOC1-binding sites [3, 9], suggesting a possible mechanism for the regulation of ABAR transcription by the clock. This was described in our model through activation of ABAR expression by LHY/CCA1 and its inhibition by TOC1 ([3, 17], eq. 29). To simplify the model and reduce the number of equations, we used a quasi-steady state approximation for the ABA complex with ABAR, assuming that ABAR protein level follows the cognate mRNA. PP2C activity was described through its inhibition by this ABA complex (eq. 30).
The second important group of reactions of ABA signalling is related with the activation of SNF1 (Sucrose-Nonfermenting Kinase1)-related protein kinase SnRK2 by ABA . This activation happens through double-negative interactions: the ABA complex inactivates PP2C and PP2C inactivates SnRK2. So we included in the model an equation for SnRK2 activity, which is inhibited by PP2C (eq. 31). After activation by ABA, SnRK2 mediates multiple reactions of the ABA pathway, including induction of transcription of multiple target genes which possess ABRE elements (ABA-responsive elements) in their promoters . As TOC1 was reported to be induced by ABA , we analysed the TOC1 promoter and identified 2 ABRE elements. This suggested that a mechanism of TOC1 induction by ABA is similar to other genes and mediated by SnRK2 (Figure 1B). Therefore we included in the model the activation of TOC1 expression by SnRK2, assuming that this activation competes with inhibition of TOC1 by LHY/CCA1 protein (, eq. 11). An analogous mechanism of ABA signalling through the ABA/PP2C/SnRK2 pathway is involved in the acceleration of stomatal closing by ABA , which we included into the model. The equation for stomatal aperture describes stomatal opening in a phenomenological way, which allowed us qualitatively compare the model behaviour with physiological data. Additionally to the inhibition of stomatal opening by ABA signalling through active SnRK2, we included the known activation of stomatal opening in the presence of light (, eq. 32).
So finally the clock model was extended by including eqs. (29)-(32) for ABAR, PP2C, SnRK2 and stomata aperture (variable “s”). Thus we connected the clock with a simple version of ABA signalling, which enabled us to simulate TOC1 and stomatal regulation by ABA. The effective parameters of the ABA-related reactions were optimized together with other parameters the clock model against multiple perturbations as described below.
Solving the system of ODEs
The extension of P2011 by adding newly discovered processes of the inhibition of the PRR9, PRR7, PRR5, LUX, ELF4 and GI clock genes by TOC1 (Figure 1A) resulted in the addition of 7 new parameters to P2011. Our model was also extended compared to P2011 by adding reactions describing ABA-related processes (Figure 1B), which introduced 14 more new parameters. Finally, the model describes the dynamics of stomatal opening, which was absent from P2011, resulting in 5 more new parameters. Therefore, the higher level of biological complexity of the model and its extension to hormonal signalling through ABA inevitably increased the parameter space, adding 26 parameters compared to P2011. This allowed our model to describe experimental data that were not described by P2011, such as the inhibition of multiple clock genes by TOC1, stimulation of TOC1 expression by ABA signalling, changing the clock period by ABA and the dynamics of stomata presented in the Results below. In addition, the new structure of the clock, which now incorporates ABA signalling to the clock through TOC1 (Figure 1A) allowed us to improve the description of the data on the clock kinetics compared to P2011, as shown in the Results below.
Parameter values were either constrained based on experimental data or fitted to multiple time-series data sets, similarly to the P2011 model. 49 out of the total 133 parameters were constrained based on the available experimental data. Most of these parameters have the same values as in P2011 because they were constrained by the same data. Other parameters were fitted to multiple datasets similarly to P2011, but their values may be different from P2011 because the structures of the models are different. The optimal set of parameters is presented in Additional file 2: Table S1. The new model, under the optimal parameter values, retains most of its properties from P2011, keeping the good fit to data from various conditions, such as diurnal cycles of the wild type and lhy/cca1, lhy/cca1/gi and elf3 mutants, as well as continuous light or darkness in the wild type and multiple mutants. For example, the period values for the simulated free-running clock in wild type and mutant plants were 24.5 h for wild type under constant light conditions, 26.6 h for wild type under constant darkness and 17.7 h, 21.4 h, 30.1 h, 30.7 h, 21.2 h for lhy/cca1, toc1, ztl, prr79 and gi mutants in constant light, which match the experimental observations [8, 19–22]. Additionally lhy/cca1/gi, elf3, elf4 and lux mutants were arrhythmic in constant light, in agreement with the data [8, 23–25]. We also tested the sensitivity of the model to the variations of the newly introduced parameters of the inhibition of target genes by TOC1 and ABA signalling. Additional file 1: Figure S2 shows the relative changes in the amplitude of LHY expression and the clock period in constant light conditions under 10% changes of each parameter. We observed less than 3% changes in the period and less than 12% change in the amplitude of LHY mRNA under variations of these parameters, which shows that the extended model retains its robustness to parameter perturbations from P2011 (Additional file 1: Figure S2). Thus, the optimal set of parameters demonstrated good correspondence to multiple datasets and provides robust behaviour in the model. However, given the significant complexity of the model, which is related to the complexity of the biological processes involved, and the number of unknown parameters, we cannot exclude the existence of other parameter sets that might equally describe the data. The detailed study of this large parameter space lies outside the scope of this paper. The model with the current parameter set matches well enough to the data to make it a useful instrument for the exploration of the possible mechanisms of the observed biological phenomenon.
Modulation of circadian dynamics via clock gene repression by TOC1
Over-expression of TOC1 also resulted in a lower level of LHY/CCA1 expression, which is due to the direct suppression of LHY/CCA1 expression by the increased level of the repressor, TOC1 (Figure 2A). This effect of TOC1 was present already in the P2011 model. However, the P2011 model failed to describe the longer period of the TOC1-overexpressor line (TOC1-ox) under constant light, which was observed in the data . The discrepancy in the earlier model was caused by increased levels of the PRR inhibitors as well as TOC1, which resulted in a reduction of the LHY/CCA1 level. This led to a shortening of the period in the simulated TOC1-ox lines in P2011 rather than a longer period. The current model describes the experimentally observed prolongation of the period due to delay in LHY/CCA1 in TOC1-ox compared to the wild type because repression by TOC1 reduces the level of the other PRRs (Figure 2C). The simulated TOC1-ox in Figure 2A has a 1.5 h longer period than the wild type. The period difference between the TOC1-ox and wild type increases with the increase in the level of TOC1 expression in our model, in agreement with the data [21, 27]. For example, 3-fold increase of the TOC1 overexpression level compared to the level shown in Figure 2A results in 5 h longer period of the TOC1-ox compared the wild type (while it has 1 h shorter period in P2011). The large prolongation of the TOC1-ox period (5 h) corresponds to the experimental observations for the TOC1-ox lines with high overexpression levels of TOC1. The oscillations are damping in TOC1-ox in our model simulations (Figure 2A) similarly to the data , to the point of the loss of rhythmicity of the simulated TOC1-ox at very high levels of TOC1-overexpression, which also agrees with the data . The model suggests that the period increase is due to the simultaneous up-regulation of TOC1 and down-regulation of the other PRR expression in the TOC1-ox lines, which allows a longer interval of LHY/CCA1 expression than in the wild type.
Figure 2C shows schematically the mechanisms of regulation of LHY/CCA1 expression by the wave of PRR inhibitors in lines with different levels of TOC1: In the toc1 mutant, increased PRR levels reduce LHY/CCA1 levels, while the absence of TOC1, the last inhibitor, shortens the period. In TOC1-ox, LHY/CCA1 expression is repressed mainly by TOC1 protein, which declines when LHY is raising (LHY and TOC1 are anti-phasic) and thus allowing oscillations to occur but with longer period than in wild type (Figure 2A). Thus the model offers an explanation of the low LHY and CCA1 levels in the toc1 mutant  that is consistent with the biochemical data showing that TOC1 is a repressor [3, 7]. The mechanism depends upon the inhibition of multiple clock genes by TOC1, particularly the PRRs (Figure 2C).
Additionally to affecting the expression of the morning genes LHY/CCA1, PRR9, PRR7, PRR5, TOC1 negatively regulates the evening genes LUX (LUX ARRHYTHMO), ELF4 (EARLY FLOWERING 4) and GI (GIGANTEA) in the experimental data  and in our model. This is in clear contrast with the P2011 model, where TOC1 affected only LHY/CCA1 expression directly. The new model suggests that increased levels of evening gene expression result in the arrhythmia of the lhy/cca1/toc1 triple mutant under constant light conditions, as opposed to the damped oscillation of the lhy/cca1 double mutant. This agrees with experimental observations .
TOC1-mediated effects of ABA on circadian timing
Altogether, the extension of the P2011 model including the widespread repression of core components by TOC1 (Figure 1A) and the inclusion of the regulatory function of ABA on TOC1 expression allowed us to improve the description of multiple datasets on the clock kinetics in wild type and mutant plants, and suggest mechanisms to explain previously paradoxical data.
Gating of ABA signalling by the clock
The revised clock model explains clock gene profiles of mutants with altered TOC1 expression, resulting from the pervasive transcriptional repression by TOC1. The model was extended to explore the mutual connections between TOC1 targets and ABA signalling, resulting in downregulation of clock gene expression and lengthening of the free running period under high ABA conditions. Moreover, TOC1’s participation in the rhythmic gating of ABA signalling can explain well-characterised rhythms in physiological processes downstream of ABA, such as the regulation of stomata aperture. The clock effects on stress signalling are significant at the whole-plant level, as plant survival under dry conditions was reduced in TOC1-ox plants and increased in toc1 mutants . Moreover, the close overlap between TOC1- and ABA-regulated genes in transcriptome profiles  suggests that other ABA-mediated stress responses are similarly responsive to TOC1. Integrating such physiological response pathways, even if they are represented in a simplified way, with the clock model will allow to expand our knowledge on the temporal regulation of plant physiology in its normal context, the unavoidable day-night cycle.
The system of ordinary differential equations was solved using MATLAB, integrated with the stiff solver ode15s (The MathWorks UK, Cambridge). The SBML version of the model will be available upon publication from the Biomodels database  (accession number BIOMD0000000445) and the Plant Systems Modelling portal (http://www.plasmo.ed.ac.uk). A MATLAB version of the model is available from the authors upon request.
LATE ELONGATED HYPOCOTYL
CIRCADIAN CLOCK ASSOCIATED 1
PRR7, PRR5, PSEUDO-RESPONSE REGULATORs 9, 7, 5
TIMING OF CAB EXPRESSION 1
EARLY FLOWERING 3
EARLY FLOWERING 4
CONSTITUTIVE PHOTOMORPHOGENIC 1
ABA receptor (magnesium-chelatase H subunit)
ABA-responsive elements on gene promoters
Protein phosphatase of 2C type
SNF1 (Sucrose-Nonfermenting Kinase1)-related protein kinase
- 12 L:
12D: 12 h light/12 h dark daily cycle of light.
This work was supported by the European Commission FP7 Collaborative Project TiMet (project 245143). SynthSys is a Centre for Integrative and Systems Biology supported by BBSRC and EPSRC award D019621. Work in P.M. laboratory is supported by grants from the Ramón Areces Foundation, from the Spanish Ministry of Science and Innovation (MICINN) (BIO2010-16483) and from EUROHORCS (European Heads Of Research Councils) and the European Science Foundation (ESF) through the EURYI Award.
- Zhang EE, Kay SA: Clocks not winding down: unravelling circadian networks. Nat Rev Mol Cell Biol. 2010, 11: 764-776. 10.1038/nrm2995.PubMedView ArticleGoogle Scholar
- Pokhilko A, Fernandez AP, Edwards KD, Southern MM, Halliday KJ, Millar AJ: The clock gene circuit in Arabidopsis includes a repressilator with additional feedback loops. Mol Syst Biol. 2012, 8: 574-587.PubMedPubMed CentralView ArticleGoogle Scholar
- Huang W, Perez-Garcia P, Pokhilko A, Millar AJ, Antoshechkin I, Riechmann JL, Mas P: Mapping the core of the arabidopsis circadian clock defines the network structure of the oscillator. Science. 2012, 336: 75-79. 10.1126/science.1219075.PubMedView ArticleGoogle Scholar
- Alabadi D, Oyama T, Yanovsky MJ, Harmon FG, Mas P, Kay SA: Reciprocal regulation between TOC1 and LHY/CCA1 within the Arabidopsis circadian clock. Science. 2001, 293: 880-883. 10.1126/science.1061320.PubMedView ArticleGoogle Scholar
- Pruneda-Paz JL, Breton G, Para A, Kay SA: A functional genomics approach reveals CHE as a component of the Arabidopsis circadian clock. Science. 2009, 323: 1481-1485. 10.1126/science.1167206.PubMedPubMed CentralView ArticleGoogle Scholar
- Baudry A, Ito S, Song YH, Strait AA, Kiba T, Lu S, Henriques R, Pruneda-Paz JL, Chua NH, Tobin EM: F-box proteins FKF1 and LKP2 act in concert with ZEITLUPE to control Arabidopsis clock progression. Plant Cell. 2010, 22: 606-622. 10.1105/tpc.109.072843.PubMedPubMed CentralView ArticleGoogle Scholar
- Gendron JM, Pruneda-Paz JL, Doherty CJ, Gross AM, Kang SE, Kay SA: Arabidopsis circadian clock protein, TOC1, is a DNA-binding transcription factor. Proc Natl Acad Sci USA. 2012, 109: 21152-21157. 10.1073/pnas.1210799110.PubMedPubMed CentralView ArticleGoogle Scholar
- Locke JC, Southern MM, Kozma-Bognar L, Hibberd V, Brown PE, Turner MS, Millar AJ: Extension of a genetic network model by iterative experimentation and mathematical analysis. Mol Syst Biol. 2005, 1: 0013-PubMedView ArticleGoogle Scholar
- Legnaioli T, Cuevas J, Mas P: TOC1 functions as a molecular switch connecting the circadian clock with plant responses to drought. EMBO J. 2009, 28: 3745-3757. 10.1038/emboj.2009.297.PubMedPubMed CentralView ArticleGoogle Scholar
- Correia MJ, Pereira JS, Chaves MM, Rodrigues ML, Pacheco CA: ABA xylem concentrations determine maximum daily leaf conductance of field-grown Vitis vinifera L. plants. Plant Cell Environ. 1995, 18: 511-521. 10.1111/j.1365-3040.1995.tb00551.x.View ArticleGoogle Scholar
- Raghavendra AS, Gonugunta VK, Christmann A, Grill E: ABA perception and signalling. Trends Plant Sci. 2010, 15: 395-401. 10.1016/j.tplants.2010.04.006.PubMedView ArticleGoogle Scholar
- Shen YY, Wang XF, Wu FQ, Du SY, Cao Z, Shang Y, Wang XL, Peng CC, Yu XC, Zhu SY: The Mg-chelatase H subunit is an abscisic acid receptor. Nature. 2006, 443: 823-826. 10.1038/nature05176.PubMedView ArticleGoogle Scholar
- Wu FQ, Xin Q, Cao Z, Liu ZQ, Du SY, Mei C, Zhao CX, Wang XF, Shang Y, Jiang T: The magnesium-chelatase H subunit binds abscisic acid and functions in abscisic acid signaling: new evidence in Arabidopsis. Plant Physiol. 2009, 150: 1940-1954. 10.1104/pp.109.140731.PubMedPubMed CentralView ArticleGoogle Scholar
- Adams R, Clark A, Yamaguchi A, Hanlon N, Tsorman N, Ali S, Lebedeva G, Goltsov A, Sorokin A, Akman OE: SBSI: an extensible distributed software infrastructure for parameter estimation in systems biology. Bioinformatics. 2013, 29: 664-665. 10.1093/bioinformatics/btt023.PubMedPubMed CentralView ArticleGoogle Scholar
- Pokhilko A, Ramos JA, Holtan H, Maszle DR, Khanna R, Millar AJ: Ubiquitin ligase switch in plant photomorphogenesis: A hypothesis. J Theor Biol. 2011, 270: 31-41. 10.1016/j.jtbi.2010.11.021.PubMedPubMed CentralView ArticleGoogle Scholar
- Neales TF, Masia A, Zhang J, Davies WJ: The effects of partially drying part of the root system of helianthus annuus on the abscisic acid content of the roots, xylem sap and leaves. J Exp Botan. 1989, 40: 1113-1120. 10.1093/jxb/40.10.1113.View ArticleGoogle Scholar
- Harmer SL, Kay SA: Positive and negative factors confer phase-specific circadian regulation of transcription in Arabidopsis. Plant Cell. 2005, 17: 1926-1940. 10.1105/tpc.105.033035.PubMedPubMed CentralView ArticleGoogle Scholar
- Dodd AN, Parkinson K, Webb AAR: Independent circadian regulation of assimilation and stomatal conductance in the ztl-1 mutant of Arabidopsis. New Phytol. 2004, 162: 63-70. 10.1111/j.1469-8137.2004.01005.x.View ArticleGoogle Scholar
- Farre EM, Harmer SL, Harmon FG, Yanovsky MJ, Kay SA: Overlapping and distinct roles of PRR7 and PRR9 in the Arabidopsis circadian clock. Curr Biol. 2005, 15: 47-54. 10.1016/j.cub.2004.12.067.PubMedView ArticleGoogle Scholar
- Mas P, Kim WY, Somers DE, Kay SA: Targeted degradation of TOC1 by ZTL modulates circadian function in Arabidopsis thaliana. Nature. 2003, 426: 567-570. 10.1038/nature02163.PubMedView ArticleGoogle Scholar
- Mas P, Alabadi D, Yanovsky MJ, Oyama T, Kay SA: Dual role of TOC1 in the control of circadian and photomorphogenic responses in Arabidopsis. Plant Cell. 2003, 15: 223-236. 10.1105/tpc.006734.PubMedPubMed CentralView ArticleGoogle Scholar
- Strayer C, Oyama T, Schultz TF, Raman R, Somers DE, Mas P, Panda S, Kreps JA, Kay SA: Cloning of the Arabidopsis clock gene TOC1, an autoregulatory response regulator homolog. Science. 2000, 289: 768-771. 10.1126/science.289.5480.768.PubMedView ArticleGoogle Scholar
- McWatters HG, Bastow RM, Hall A, Millar AJ: The ELF3 zeitnehmer regulates light signalling to the circadian clock. Nature. 2000, 408: 716-720. 10.1038/35047079.PubMedView ArticleGoogle Scholar
- McWatters HG, Kolmos E, Hall A, Doyle MR, Amasino RM, Gyula P, Nagy F, Millar AJ, Davis SJ: ELF4 is required for oscillatory properties of the circadian clock. Plant Physiol. 2007, 144: 391-401. 10.1104/pp.107.096206.PubMedPubMed CentralView ArticleGoogle Scholar
- Hazen SP, Schultz TF, Pruneda-Paz JL, Borevitz JO, Ecker JR, Kay SA: LUX ARRHYTHMO encodes a Myb domain protein essential for circadian rhythms. Proc Natl Acad Sci USA. 2005, 102: 10387-10392. 10.1073/pnas.0503029102.PubMedPubMed CentralView ArticleGoogle Scholar
- Martin-Tryon EL, Kreps JA, Harmer SL: GIGANTEA acts in blue light signaling and has biochemically separable roles in circadian clock and flowering time regulation. Plant Physiol. 2007, 143: 473-486.PubMedPubMed CentralView ArticleGoogle Scholar
- Makino S, Matsushika A, Kojima M, Yamashino T, Mizuno T: The APRR1/TOC1 quintet implicated in circadian rhythms of Arabidopsis thaliana: I. Characterization with APRR1-overexpressing plants. Plant Cell Physiol. 2002, 43: 58-69. 10.1093/pcp/pcf005.PubMedView ArticleGoogle Scholar
- Ding Z, Doyle MR, Amasino RM, Davis SJ: A complex genetic interaction between Arabidopsis thaliana TOC1 and CCA1/LHY in driving the circadian clock and in output regulation. Genetics. 2007, 176: 1501-1510. 10.1534/genetics.107.072769.PubMedPubMed CentralView ArticleGoogle Scholar
- Hanano S, Domagalska MA, Nagy F, Davis SJ: Multiple phytohormones influence distinct parameters of the plant circadian clock. Genes Cells. 2006, 11: 1381-1392. 10.1111/j.1365-2443.2006.01026.x.PubMedView ArticleGoogle Scholar
- Webb AAR: Stomatal rhythms. Biological rhythms and photoperiodism in plants. Edited by: Lumsden PJ, Millar AJ. 1998, Oxford: Bios Scientific Publications, 66-79.Google Scholar
- Le Novere N, Bornstein B, Broicher A, Courtot M, Donizelli M, Dharuri H, Li L, Sauro H, Schilstra M, Shapiro B: BioModels Database: a free, centralized database of curated, published, quantitative kinetic models of biochemical and cellular systems. Nucleic Acids Res. 2006, 34: D689-691. 10.1093/nar/gkj092.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.