Deterministic and stochastic models of Arabidopsis Thaliana flowering

Emrah Haspolat, Benoit Huard, Maia Angelova

Research output: Contribution to journalArticlepeer-review

5 Citations (Scopus)
39 Downloads (Pure)


Experimental studies of the flowering of Arabidopsis Thaliana have shown that
a large complex gene regulatory network (GRN) is responsible for its regulation.
This process has been mathematically modelled with deterministic differential
equations by considering the interactions between gene activators and inhibitors
[26, 28]. However, due to complexity of the model, the properties of the network
and the roles of the individual genes cannot be deducted from the numerical
solution the published work offers. Here, we propose simplifications of the model, based on decoupling of the original GRN to motifs, described with three and two differential equations. A stable solution of the original model is sought by linearisation of the original model which contributes to further investigation of the role of the individual genes to the flowering. Furthermore, we study the role of noise by introducing and investigating two types of stochastic elements into the model.

The deterministic and stochastic non-linear dynamic models of Arabidopsis flowering time are considered by following the deterministic delayed model introduced in [26]. Steady state regimes and stability of the deterministic original
model are investigated analytically and numerically. By decoupling some concentrations, the system was reduced to emphasise the role played by the transcription factor Suppressor of Overexpression of Constants1 (SOC1) and the important floral meristem identity genes, Leafy (LFY ) and Apetala1 (AP1). Two-dimensional motifs, based on the dynamics of LFY and AP1, are obtained from the reduced network and parameter ranges ensuring flowering are determined. Their stability analysis shows that LFY and AP1 are regulating each other for flowering, matching experimental findings. New sufficient conditions of mean square stability in the stochastic model are obtained using a stochastic Lyapunov approach. Our numerical simulations demonstrate that the reduced models of Arabidopsis flowering time, describing specific motifs of the GRN, can capture the essential behaviour of the full system and also introduce the conditions of
flowering initiation. Additionally, they show that stochastic effects can change the behaviour of the stability region through a stability switch. This study thus contributes to a better understanding of the role of LFY and AP1 in Arabidopsis
Original languageEnglish
Pages (from-to)277-311
Number of pages35
JournalBulletin of Mathematical Biology
Issue number1
Early online date8 Nov 2018
Publication statusPublished - 15 Jan 2019


Dive into the research topics of 'Deterministic and stochastic models of Arabidopsis Thaliana flowering'. Together they form a unique fingerprint.

Cite this