Modelling Cellular Kinematics in Self-Similar Plant Growth

Biophysics Internship - 2025

Cover image credit from Figure 3 of (Prusinkiewicz & Runions, 2012).

Some plant structures exhibit self-similar growth, that is, they maintain stable geometries as their cells grow, divide, and flow within them. This offers an analytically tractable system for exploring their growth kinematics and the coupling between genetic regulation and morphoelasticity in plant development. After first trying to examine the PDEs governing the system analytically we then used a custom numerical solver for morphodynamics, namely Organism, to create models to emulate these forms, capturing cellular growth and proliferation while preserving overall geometry.

Specifically, we looked at emulating the Shoot Apical Meristem (SAM) and the apical hook as shown in Figure 1. We did this by performing a map between Lagrangian, curvilinear and Eulerian coordinate systems. A parabolic map was used to model the SAM and is given by \(\begin{equation} \begin{split} s &\mapsto x = sr \\ r &\mapsto y = \frac{1}{2}(r^2 - s^2) \,, \label{ODE} \end{split} \end{equation}\) where $(s,r)$ span the curvilinear space and $(x,y)$ span the Eulerian space. The apical hook required the more complicated transformation \(\begin{equation} \begin{split} \theta(s) &= \int\limits_0^s d\tilde{s} \, \kappa(\tilde{s}) \\ x(s,r) &= \int\limits_0^s d\tilde{s} \,\sin(\theta(\tilde{s})) + r\cos(\theta(s)) \\ y(s,r) &= \int\limits_0^s d\tilde{s} \, \cos(\theta(\tilde{s})) - r\sin(\theta(s)) \end{split} \end{equation}\) where $\kappa=\kappa(\tilde{s})$ is the curvature distribution–we used a gamma distribution.

Figure 1. (Left) Illustration of the SAM (Moore J, 2025). (Right) Image of the apical hook (Huang et al., 2020).

The resultant models are shown in Figure 2. The cells flow within a fixed envelope while growing and dividing. By keeping the curvilinear space the same (up to length scales) and solely varying the map to the Eulerian space, completely different plant structures can be replicated. Gene regulatory networks can then be placed on these structures and coupled to growth mechanics in the orthogonal, $r$, direction.

Figure 2. Models of the SAM (left) and apical hook (Right).

References

2025

  1. moore2025vcu.png
    VCU BIOL 152: Introduction to Biological Sciences II
    Jennings D Moore J
    2025
    Richmond (VA): Virginia Commonwealth University

2020

  1. huang2020salicylic.png
    Salicylic acid suppresses apical hook formation via NPR1-mediated repression of EIN3 and EIL1 in Arabidopsis
    P Huang, Z Dong, P Guo, and 5 more authors
    The Plant Cell, 2020

2012

  1. prusinkiewicz2012computational.png
    Computational models of plant development and form
    P Prusinkiewicz and A Runions
    New Phytologist, 2012