Reaction-diffusion patterns in plant tip morphogenesis: bifurcations on spherical caps

Size regulation of the lateral organ initiation zone and its role in determining cotyledon number in conifers

Introduction

Unlike monocots and dicots, many conifers, particularly Pinaceae, form three or more cotyledons. These are arranged in a whorl, or ring, at a particular distance from the embryo tip, with cotyledons evenly spaced within the ring. The number of cotyledons, n_c, varies substantially within species, both in clonal cultures and in seed embryos. n_c variability reflects embryo size variability, with larger diameter embryos having higher n_c. Correcting for growth during embryo development, we extract values for the whorl radius at each n_c. This radius, corresponding to the spatial pattern of cotyledon differentiation factors, varies over three-fold for the naturally observed range of n_c. The current work focuses on factors in the patterning mechanism that could produce such a broad variability in whorl radius. Molecularly, work in Arabidopsis has shown that the initiation zone for leaf primordia occurs at a minimum between inhibitor zones of HD-ZIP III at the shoot apical meristem (SAM) tip and KANADI (KAN) encircling this farther from the tip. PIN1-auxin dynamics within this uninhibited ring form auxin maxima, specifying primordia initiation sites. A similar mechanism is indicated in conifer embryos by effects on cotyledon formation with overexpression of HD-ZIP III inhibitors and by interference with PIN1-auxin patterning.

Methods

We develop a mathematical model for HD-ZIP III/KAN spatial localization and use this to characterize the molecular regulation that could generate (a) the three-fold whorl radius variation (and associated n_c variability) observed in conifer cotyledon development, and (b) the HD-ZIP III and KAN shifts induced experimentally in conifer embryos and in Arabidopsis.

Results

This quantitative framework indicates the sensitivity of mechanism components for positioning lateral organs closer to or farther from the tip. Positional shifting is most readily driven by changes to the extent of upstream (meristematic) patterning and changes in HD-ZIP III/KAN mutual inhibition, and less efficiently driven by changes in upstream dosage or the activation of HD-ZIP III. Sharper expression boundaries can also be more resistant to shifting than shallower expression boundaries.

Discussion

The strong variability seen in conifer n_c (commonly from 2 to 10) may reflect a freer variation in regulatory interactions, whereas monocot (n_c = 1) and dicot (n_c = 2) development may require tighter control of such variation. These results provide direction for future quantitative experiments on the positional control of lateral organ initiation, and consequently on plant phyllotaxy and architecture.

Pattern Propagation Driven by Surface Curvature

Pattern dynamics on curved surfaces are found everywhere in nature. The geometry of surfaces has been shown to influence dynamics and play a functional role, yet a comprehensive understanding is still elusive. Here, we report for the first time that a static Turing pattern on a flat surface can propagate on a curved surface, as opposed to previous studies, where the pattern is presupposed to be static irrespective of the surface geometry. To understand such significant changes on curved surfaces, we investigate reaction-diffusion systems on axisymmetric curved surfaces. Numerical and theoretical analyses reveal that both the symmetries of the surface and pattern participate in the initiation of pattern propagation. This study provides a novel and generic mechanism of pattern propagation that is caused by surface curvature, as well as insights into the general role of surface geometry.

Two-stage patterning dynamics in conifer cotyledon whorl morphogenesis

Background and Aims

Conifer embryos, unlike those of monocots or dicots, have variable numbers of cotyledons, even within the same species. Cotyledons form in a single whorl on a dome-shaped embryo. The closely spaced cotyledons are not found outside this ring, indicating a radial control on where they can form. Polar transport of the hormone auxin affects outgrowth of distinct cotyledons, but not the radial aspect of the whorl or the within-whorl spacing between cotyledons. A quantitative model of plant growth regulator patterning is needed to understand the dynamics of this complex morphogenetic process.

Methods

A two-stage reaction-diffusion model is developed for the spatial patterning of growth regulators on the embryo surface, with a radial pattern (P1) constraining the shorter-wavelength cotyledon pattern (P2) to a whorl. These patterns drive three-dimensional (3-D) morphogenesis by catalysing local surface growth.

Key Results

Growth driven by P2 generates single whorls across the experimentally observed range of two to 11 cotyledons, as well as the circularly symmetric response to auxin transport interference. These computations are the first corroboration of earlier theoretical proposals for hierarchical control of whorl formation. The model generates the linear relationship between cotyledon number and embryo diameter observed experimentally. This accounts for normal integer cotyledon number selection, as well as the less common cotyledon fusings and splittings observed experimentally. Flattening of the embryo during development may affect the upward outgrowth angle of the cotyledons.

Conclusions

Cotyledon morphogenesis is more complex geometrically in conifers than in angiosperms, involving 2-D patterning which deforms a surface in three dimensions. This work develops a quantitative framework for understanding the growth and patterning dynamics involved in conifer cotyledon development, and applies more generally to the morphogenesis of whorls with many primordia.

Mapping nano-scale mechanical heterogeneity of primary plant cell walls

Nanoindentation experiments are performed using an atomic force microscope (AFM) to quantify the spatial distribution of mechanical properties of plant cell walls at nanometre length scales. At any specific location on the cell wall, a complex (non-linear) force-indentation response occurs that can be deconvoluted using a unique multiregime analysis (MRA). This allows an unambiguous evaluation of the local transverse elastic modulus of the wall. Nanomechanical measurements on suspension-cultured cells (SCCs), derived from Italian ryegrass (Lolium multiflorum) starchy endosperm, show three characteristic modes of deformation and a spatial distribution of elastic moduli across the surface. 'Soft' and 'hard' domains are found across length scales between 0.1 µm and 3 µm, which is well above a typical pore size of the polysaccharide mesh. The generality and wider applicability of this mechanical heterogeneity is verified through in planta characterization on leaf epidermal cells of Arabidopsis thaliana and L. multiflorum The outcomes of this research provide a basis for uncovering and quantifying the relationships between local wall composition, architecture, cell growth, and/or morphogenesis.