Physical models of plant development.

The definition of shape in multicellular organisms is a major issue of developmental biology. It is well established that morphogenesis relies on genetic regulation. However, cells, tissues, and organism behaviors are also bound by the laws of physics, which limit the range of possible deformations organisms can undergo but also define what organisms must do to achieve specific shapes. Besides experiments, theoretical models and numerical simulations of growing tissues are powerful tools to investigate the link between genetic regulation and mechanics. Here, we provide an overview of the main mechanical models of plant morphogenesis developed so far, from subcellular scales to whole tissues. The common concepts and discrepancies between the various models are discussed.

Phyllotaxis as geometric canalization during plant development.

Why living forms develop in a relatively robust manner, despite various sources of internal or external variability, is a fundamental question in developmental biology. Part of the answer relies on the notion of developmental constraints: at any stage of ontogenesis, morphogenetic processes are constrained to operate within the context of the current organism being built. One such universal constraint is the shape of the organism itself, which progressively channels the development of the organism toward its final shape. Here, we illustrate this notion with plants, where strikingly symmetric patterns (phyllotaxis) are formed by lateral organs. This Hypothesis article aims first to provide an accessible overview of phyllotaxis, and second to argue that the spiral patterns in plants are progressively canalized from local interactions of nascent organs. The relative uniformity of the organogenesis process across all plants then explains the prevalence of certain patterns in plants, i.e. Fibonacci phyllotaxis.

Phenotyping and modeling of root hydraulic architecture reveal critical determinants of axial water transport.

Water uptake by roots is a key adaptation of plants to aerial life. Water uptake depends on root system architecture (RSA) and tissue hydraulic properties that, together, shape the root hydraulic architecture. This work investigates how the interplay between conductivities along radial (e.g. aquaporins) and axial (e.g. xylem vessels) pathways determines the water transport properties of highly branched RSAs as found in adult Arabidopsis (Arabidopsis thaliana) plants. A hydraulic model named HydroRoot was developed, based on multi-scale tree graph representations of RSAs. Root water flow was measured by the pressure chamber technique after successive cuts of a same root system from the tip toward the base. HydroRoot model inversion in corresponding RSAs allowed us to concomitantly determine radial and axial conductivities, providing evidence that the latter is often overestimated by classical evaluation based on the Hagen-Poiseuille law. Organizing principles of Arabidopsis primary and lateral root growth and branching were determined and used to apply the HydroRoot model to an extended set of simulated RSAs. Sensitivity analyses revealed that water transport can be co-limited by radial and axial conductances throughout the whole RSA. The number of roots that can be sectioned (intercepted) at a given distance from the base was defined as an accessible and informative indicator of RSA. The overall set of experimental and theoretical procedures was applied to plants mutated in ESKIMO1 and previously shown to have xylem collapse. This approach will be instrumental to dissect the root water transport phenotype of plants with intricate alterations in root growth or transport functions.

Benchmarking of deep learning algorithms for 3D instance segmentation of confocal image datasets.

Segmenting three-dimensional (3D) microscopy images is essential for understanding phenomena like morphogenesis, cell division, cellular growth, and genetic expression patterns. Recently, deep learning (DL) pipelines have been developed, which claim to provide high accuracy segmentation of cellular images and are increasingly considered as the state of the art for image segmentation problems. However, it remains difficult to define their relative performances as the concurrent diversity and lack of uniform evaluation strategies makes it difficult to know how their results compare. In this paper, we first made an inventory of the available DL methods for 3D cell segmentation. We next implemented and quantitatively compared a number of representative DL pipelines, alongside a highly efficient non-DL method named MARS. The DL methods were trained on a common dataset of 3D cellular confocal microscopy images. Their segmentation accuracies were also tested in the presence of different image artifacts. A specific method for segmentation quality evaluation was adopted, which isolates segmentation errors due to under- or oversegmentation. This is complemented with a 3D visualization strategy for interactive exploration of segmentation quality. Our analysis shows that the DL pipelines have different levels of accuracy. Two of them, which are end-to-end 3D and were originally designed for cell boundary detection, show high performance and offer clear advantages in terms of adaptability to new data.

Cortical tension overrides geometrical cues to orient microtubules in confined protoplasts.

In plant cells, cortical microtubules (CMTs) generally control morphogenesis by guiding cellulose synthesis. CMT alignment has been proposed to depend on geometrical cues, with microtubules aligning with the cell long axis in silico and in vitro. Yet, CMTs are usually transverse in vivo, i.e., along predicted maximal tension, which is transverse for cylindrical pressurized vessels. Here, we adapted a microwell setup to test these predictions in a single-cell system. We confined protoplasts laterally to impose a curvature ratio and modulated pressurization through osmotic changes. We find that CMTs can be longitudinal or transverse in wallless protoplasts and that the switch in CMT orientation depends on pressurization. In particular, longitudinal CMTs become transverse when cortical tension increases. This explains the dual behavior of CMTs in planta: CMTs become longitudinal when stress levels become low, while stable transverse CMT alignments in tissues result from their autonomous response to tensile stress fluctuations.

Sugar availability suppresses the auxin-induced strigolactone pathway to promote bud outgrowth.

Apical dominance occurs when the growing shoot tip inhibits the outgrowth of axillary buds. Apically-derived auxin in the nodal stem indirectly inhibits bud outgrowth via cytokinins and strigolactones. Recently, sugar deprivation was found to contribute to this phenomenon. Using rose and pea, we investigated whether sugar availability interacts with auxin in bud outgrowth control, and the role of cytokinins and strigolactones, in vitro and in planta. We show that sucrose antagonises auxin's effect on bud outgrowth, in a dose-dependent and coupled manner. Sucrose also suppresses strigolactone inhibition of outgrowth and the rms3 strigolactone-perception mutant is less affected by reducing sucrose supply. However, sucrose does not interfere with the regulation of cytokinin levels by auxin and stimulates outgrowth even with optimal cytokinin supply. These observations were assembled into a computational model in which sucrose represses bud response to strigolactones, largely independently of cytokinin levels. It quantitatively captures our observed dose-dependent sucrose-hormones effects on bud outgrowth and allows us to express outgrowth response to various combinations of auxin and sucrose levels as a simple quantitative law. This study places sugars in the bud outgrowth regulatory network and paves the way for a better understanding of branching plasticity in response to environmental and genotypic factors.

Coupling water fluxes with cell wall mechanics in a multicellular model of plant development.

The growth of plant organs is a complex process powered by osmosis that attracts water inside the cells; this influx induces simultaneously an elastic extension of the walls and pressure in the cells, called turgor pressure; above a threshold, the walls yield and the cells grow. Based on Lockhart's seminal work, various models of plant morphogenesis have been proposed, either for single cells, or focusing on the wall mechanical properties. However, the synergistic coupling of fluxes and wall mechanics has not yet been fully addressed in a multicellular model. This work lays the foundations of such a model, by simplifying as much as possible each process and putting emphasis on the coupling itself. Its emergent properties are rich and can help to understand plant morphogenesis. In particular, we show that the model can display a new type of lateral inhibitory mechanism that amplifies growth heterogeneities due e.g to cell wall loosening.

Cytokinin signalling inhibitory fields provide robustness to phyllotaxis.

How biological systems generate reproducible patterns with high precision is a central question in science. The shoot apical meristem (SAM), a specialized tissue producing plant aerial organs, is a developmental system of choice to address this question. Organs are periodically initiated at the SAM at specific spatial positions and this spatiotemporal pattern defines phyllotaxis. Accumulation of the plant hormone auxin triggers organ initiation, whereas auxin depletion around organs generates inhibitory fields that are thought to be sufficient to maintain these patterns and their dynamics. Here we show that another type of hormone-based inhibitory fields, generated directly downstream of auxin by intercellular movement of the cytokinin signalling inhibitor ARABIDOPSIS HISTIDINE PHOSPHOTRANSFER PROTEIN 6 (AHP6), is involved in regulating phyllotactic patterns. We demonstrate that AHP6-based fields establish patterns of cytokinin signalling in the meristem that contribute to the robustness of phyllotaxis by imposing a temporal sequence on organ initiation. Our findings indicate that not one but two distinct hormone-based fields may be required for achieving temporal precision during formation of reiterative structures at the SAM, thus indicating an original mechanism for providing robustness to a dynamic developmental system.

Simulating Turgor-Induced Stress Patterns in Multilayered Plant Tissues.

The intertwining between mechanics and developmental biology is extensively studied at the shoot apical meristem of land plants. Indeed, plant morphogenesis heavily relies on mechanics; tissue deformations are fueled by turgor-induced forces, and cell mechanosensitivity plays a major regulatory role in this dynamics. Since measurements of forces in growing meristems are still out of reach, our current knowledge relies mainly on theoretical and numerical models. So far, these modeling efforts have been mostly focusing on the epidermis, where aerial organs are initiated. In many models, the epidermis is assimilated to its outermost cell walls and described as a thin continuous shell under pressure, thereby neglecting the inner walls. There is, however, growing experimental evidence suggesting a significant mechanical role of these inner walls. The aim of this work is to investigate the influence of inner walls on the mechanical homeostasis of meristematic tissues. To this end, we simulated numerically the effect of turgor-induced loading, both in realistic flower buds and in more abstract structures. These simulations were performed using finite element meshes with subcellular resolution. Our analysis sheds light on the mechanics of growing plants by revealing the strong influence of inner walls on the epidermis mechanical stress pattern especially in negatively curved regions. Our simulations also display some strong and unsuspected features, such as a correlation between stress intensity and cell size, as well as differential response to loading between epidermal and inner cells. Finally, we monitored the time evolution of the mechanical stresses felt by each cell and its descendants during the early steps of flower morphogenesis.

Regulation of plant cell wall stiffness by mechanical stress: a mesoscale physical model.

A crucial question in developmental biology is how cell growth is coordinated in living tissue to generate complex and reproducible shapes. We address this issue here in plants, where stiff extracellular walls prevent cell migration and morphogenesis mostly results from growth driven by turgor pressure. How cells grow in response to pressure partly depends on the mechanical properties of their walls, which are generally heterogeneous, anisotropic and dynamic. The active control of these properties is therefore a cornerstone of plant morphogenesis. Here, we focus on wall stiffness, which is under the control of both molecular and mechanical signaling. Indeed, in plant tissues, the balance between turgor and cell wall elasticity generates a tissue-wide stress field. Within cells, mechano-sensitive structures, such as cortical microtubules, adapt their behavior accordingly and locally influence cell wall remodeling dynamics. To fully apprehend the properties of this feedback loop, modeling approaches are indispensable. To that end, several modeling tools in the form of virtual tissues have been developed. However, these models often relate mechanical stress and cell wall stiffness in relatively abstract manners, where the molecular specificities of the various actors are not fully captured. In this paper, we propose to refine this approach by including parsimonious biochemical and biomechanical properties of the main molecular actors involved. Through a coarse-grained approach and through finite element simulations, we study the role of stress-sensing microtubules on organ-scale mechanics.

Systems analysis of auxin transport in the Arabidopsis root apex.

Auxin is a key regulator of plant growth and development. Within the root tip, auxin distribution plays a crucial role specifying developmental zones and coordinating tropic responses. Determining how the organ-scale auxin pattern is regulated at the cellular scale is essential to understanding how these processes are controlled. In this study, we developed an auxin transport model based on actual root cell geometries and carrier subcellular localizations. We tested model predictions using the DII-VENUS auxin sensor in conjunction with state-of-the-art segmentation tools. Our study revealed that auxin efflux carriers alone cannot create the pattern of auxin distribution at the root tip and that AUX1/LAX influx carriers are also required. We observed that AUX1 in lateral root cap (LRC) and elongating epidermal cells greatly enhance auxin's shootward flux, with this flux being predominantly through the LRC, entering the epidermal cells only as they enter the elongation zone. We conclude that the nonpolar AUX1/LAX influx carriers control which tissues have high auxin levels, whereas the polar PIN carriers control the direction of auxin transport within these tissues.

Integration of hormonal signaling networks and mobile microRNAs is required for vascular patterning in Arabidopsis roots.

As multicellular organisms grow, positional information is continually needed to regulate the pattern in which cells are arranged. In the Arabidopsis root, most cell types are organized in a radially symmetric pattern; however, a symmetry-breaking event generates bisymmetric auxin and cytokinin signaling domains in the stele. Bidirectional cross-talk between the stele and the surrounding tissues involving a mobile transcription factor, SHORT ROOT (SHR), and mobile microRNA species also determines vascular pattern, but it is currently unclear how these signals integrate. We use a multicellular model to determine a minimal set of components necessary for maintaining a stable vascular pattern. Simulations perturbing the signaling network show that, in addition to the mutually inhibitory interaction between auxin and cytokinin, signaling through SHR, microRNA165/6, and PHABULOSA is required to maintain a stable bisymmetric pattern. We have verified this prediction by observing loss of bisymmetry in shr mutants. The model reveals the importance of several features of the network, namely the mutual degradation of microRNA165/6 and PHABULOSA and the existence of an additional negative regulator of cytokinin signaling. These components form a plausible mechanism capable of patterning vascular tissues in the absence of positional inputs provided by the transport of hormones from the shoot.

Root system markup language: toward a unified root architecture description language.

The number of image analysis tools supporting the extraction of architectural features of root systems has increased in recent years. These tools offer a handy set of complementary facilities, yet it is widely accepted that none of these software tools is able to extract in an efficient way the growing array of static and dynamic features for different types of images and species. We describe the Root System Markup Language (RSML), which has been designed to overcome two major challenges: (1) to enable portability of root architecture data between different software tools in an easy and interoperable manner, allowing seamless collaborative work; and (2) to provide a standard format upon which to base central repositories that will soon arise following the expanding worldwide root phenotyping effort. RSML follows the XML standard to store two- or three-dimensional image metadata, plant and root properties and geometries, continuous functions along individual root paths, and a suite of annotations at the image, plant, or root scale at one or several time points. Plant ontologies are used to describe botanical entities that are relevant at the scale of root system architecture. An XML schema describes the features and constraints of RSML, and open-source packages have been developed in several languages (R, Excel, Java, Python, and C#) to enable researchers to integrate RSML files into popular research workflow.

Algorithmic height compression of unordered trees.

By nature, tree structures frequently present similarities between their sub-parts. Making use of this redundancy, different types of tree compression techniques have been designed in the literature to reduce the complexity of tree structures. A popular and efficient way to compress a tree consists of merging its isomorphic subtrees, which produces a directed acyclic graph (DAG) equivalent to the original tree. An important property of this method is that the compressed structure (i.e. the DAG) has the same height as the original tree, thus limiting partially the possibility of compression. In this paper we address the problem of further compressing this DAG in height. The difficulty is that compression must be carried out on substructures that are not exactly isomorphic as they are strictly nested within each-other. We thus introduced a notion of quasi-isomorphism between subtrees that makes it possible to define similar patterns along any given path in a tree. We then proposed an algorithm to detect these patterns and to merge them, thus leading to compressed structures corresponding to DAGs augmented with return edges. In this way, redundant information is removed from the original tree in both width and height, thus achieving minimal structural compression. The complete compression algorithm is then illustrated on the compression of various plant-like structures.

A computational framework for 3D mechanical modeling of plant morphogenesis with cellular resolution.

The link between genetic regulation and the definition of form and size during morphogenesis remains largely an open question in both plant and animal biology. This is partially due to the complexity of the process, involving extensive molecular networks, multiple feedbacks between different scales of organization and physical forces operating at multiple levels. Here we present a conceptual and modeling framework aimed at generating an integrated understanding of morphogenesis in plants. This framework is based on the biophysical properties of plant cells, which are under high internal turgor pressure, and are prevented from bursting because of the presence of a rigid cell wall. To control cell growth, the underlying molecular networks must interfere locally with the elastic and/or plastic extensibility of this cell wall. We present a model in the form of a three dimensional (3D) virtual tissue, where growth depends on the local modulation of wall mechanical properties and turgor pressure. The model shows how forces generated by turgor-pressure can act both cell autonomously and non-cell autonomously to drive growth in different directions. We use simulations to explore lateral organ formation at the shoot apical meristem. Although different scenarios lead to similar shape changes, they are not equivalent and lead to different, testable predictions regarding the mechanical and geometrical properties of the growing lateral organs. Using flower development as an example, we further show how a limited number of gene activities can explain the complex shape changes that accompany organ outgrowth.

Lateral root morphogenesis is dependent on the mechanical properties of the overlaying tissues.

In Arabidopsis, lateral root primordia (LRPs) originate from pericycle cells located deep within the parental root and have to emerge through endodermal, cortical, and epidermal tissues. These overlaying tissues place biomechanical constraints on the LRPs that are likely to impact their morphogenesis. This study probes the interplay between the patterns of cell division, organ shape, and overlaying tissues on LRP morphogenesis by exploiting recent advances in live plant cell imaging and image analysis. Our 3D/4D image analysis revealed that early stage LRPs exhibit tangential divisions that create a ring of cells corralling a population of rapidly dividing cells at its center. The patterns of division in the latter population of cells during LRP morphogenesis are not stereotypical. In contrast, statistical analysis demonstrated that the shape of new LRPs is highly conserved. We tested the relative importance of cell division pattern versus overlaying tissues on LRP morphogenesis using mutant and transgenic approaches. The double mutant aurora1 (aur1) aur2 disrupts the pattern of LRP cell divisions and impacts its growth dynamics, yet the new organ's dome shape remains normal. In contrast, manipulating the properties of overlaying tissues disrupted LRP morphogenesis. We conclude that the interaction with overlaying tissues, rather than the precise pattern of divisions, is most important for LRP morphogenesis and optimizes the process of lateral root emergence.

Systems analysis of shoot apical meristem growth and development: integrating hormonal and mechanical signaling.

The shoot apical meristem (SAM) is a small population of stem cells that continuously generates organs and tissues. This review covers our current understanding of organ initiation by the SAM in Arabidopsis thaliana. Meristem function and maintenance involves two major hormones, cytokinins and auxins. Cytokinins appear to play a major role in meristem maintenance and in controlling meristematic properties, such as cell proliferation. Self-organizing transport processes, which are still only partially understood, lead to the patterned accumulation of auxin at particular positions, where organs will grow out. A major downstream target of auxin-mediated growth regulation is the cell wall, which is a determinant for both growth rates and growth distribution, but feedbacks with metabolism and the synthetic capacity of the cytoplasm are crucial as well. Recent work has also pointed at a potential role of mechanical signals in growth coordination, but the precise mechanisms at work remain to be elucidated.

Multiscale systems analysis of root growth and development: modeling beyond the network and cellular scales.

Over recent decades, we have gained detailed knowledge of many processes involved in root growth and development. However, with this knowledge come increasing complexity and an increasing need for mechanistic modeling to understand how those individual processes interact. One major challenge is in relating genotypes to phenotypes, requiring us to move beyond the network and cellular scales, to use multiscale modeling to predict emergent dynamics at the tissue and organ levels. In this review, we highlight recent developments in multiscale modeling, illustrating how these are generating new mechanistic insights into the regulation of root growth and development. We consider how these models are motivating new biological data analysis and explore directions for future research. This modeling progress will be crucial as we move from a qualitative to an increasingly quantitative understanding of root biology, generating predictive tools that accelerate the development of improved crop varieties.

A data-driven integrative model of sepal primordium polarity in Arabidopsis.

Flower patterning is determined by a complex molecular network but how this network functions remains to be elucidated. Here, we develop an integrative modeling approach that assembles heterogeneous data into a biologically coherent model to allow predictions to be made and inconsistencies among the data to be found. We use this approach to study the network underlying sepal development in the young flower of Arabidopsis thaliana. We constructed a digital atlas of gene expression and used it to build a dynamical molecular regulatory network model of sepal primordium development. This led to the construction of a coherent molecular network model for lateral organ polarity that fully recapitulates expression and interaction data. Our model predicts the existence of three novel pathways involving the HD-ZIP III genes and both cytokinin and ARGONAUTE family members. In addition, our model provides predictions on molecular interactions. In a broader context, this approach allows the extraction of biological knowledge from diverse types of data and can be used to study developmental processes in any multicellular organism.

Functional-structural plant models: a growing paradigm for plant studies.

A number of research groups in various areas of plant biology as well as computer science and applied mathematics have addressed modelling the spatiotemporal dynamics of growth and development of plants. This has resulted in development of functional-structural plant models (FSPMs). In FSPMs, the plant structure is always explicitly represented in terms of a network of elementary units. In this respect, FSPMs are different from more abstract models in which a simplified representation of the plant structure is frequently used (e.g. spatial density of leaves, total biomass, etc.). This key feature makes it possible to build modular models and creates avenues for efficient exchange of model components and experimental data. They are being used to deal with the complex 3-D structure of plants and to simulate growth and development occurring at spatial scales from cells to forest areas, and temporal scales from seconds to decades and many plant generations. The plant types studied also cover a broad spectrum, from algae to trees. This special issue of Annals of Botany features selected papers on FSPM topics such as models of morphological development, models of physical and biological processes, integrated models predicting dynamics of plants and plant communities, modelling platforms, methods for acquiring the 3-D structures of plants using automated measurements, and practical applications for agronomic purposes.