Developmental timing in plants.

Plants exhibit reproducible timing of developmental events at multiple scales, from switches in cell identity to maturation of the whole plant. Control of developmental timing likely evolved for similar reasons that humans invented clocks: to coordinate events. However, whereas clocks are designed to run independently of conditions, plant developmental timing is strongly dependent on growth and environment. Using simplified models to convey key concepts, we review how growth-dependent and inherent timing mechanisms interact with the environment to control cyclical and progressive developmental transitions in plants.

The hidden diversity of vascular patterns in flower heads.

Vascular systems are intimately related to the shape and spatial arrangement of the plant organs they support. We investigate the largely unexplored association between spiral phyllotaxis and the vascular system in Asteraceae flower heads. We imaged heads of eight species using synchrotron-based X-ray micro-computed tomography and applied original virtual reality and haptic software to explore head vasculature in three dimensions. We then constructed a computational model to infer a plausible patterning mechanism. The vascular system in the head of the model plant Gerbera hybrida is qualitatively different from those of Bellis perennis and Helianthus annuus, characterized previously. Cirsium vulgare, Craspedia globosa, Echinacea purpurea, Echinops bannaticus, and Tanacetum vulgare represent variants of the Bellis and Helianthus systems. In each species, the layout of the main strands is stereotypical, but details vary. The observed vascular patterns can be generated by a common computational model with different parameter values. In spite of the observed differences of vascular systems in heads, they may be produced by a conserved mechanism. The diversity and irregularities of vasculature stand in contrast with the relative uniformity and regularity of phyllotactic patterns, confirming that phyllotaxis in heads is not driven by the vasculature.

Diversity of woodland strawberry inflorescences arises from heterochrony regulated by TERMINAL FLOWER 1 and FLOWERING LOCUS T.

A vast variety of inflorescence architectures have evolved in angiosperms. Here, we analyze the diversity and development of the woodland strawberry (Fragaria vesca) inflorescence. Contrary to historical classifications, we show that it is a closed thyrse: a compound inflorescence with determinate primary monopodial axis and lateral sympodial branches, thus combining features of racemes and cymes. We demonstrate that this architecture is generated by 2 types of inflorescence meristems differing in their geometry. We further show that woodland strawberry homologs of TERMINAL FLOWER 1 (FvTFL1) and FLOWERING LOCUS T (FvFT1) regulate the development of both the racemose and cymose components of the thyrse. Loss of functional FvTFL1 reduces the number of lateral branches of the main axis and iterations in the lateral branches but does not affect their cymose pattern. These changes can be enhanced or compensated by altering FvFT1 expression. We complement our experimental findings with a computational model that captures inflorescence development using a small set of rules. The model highlights the distinct regulation of the fate of the primary and higher-order meristems, and explains the phenotypic diversity among inflorescences in terms of heterochrony resulting from the opposite action of FvTFL1 and FvFT1 within the thyrse framework. Our results represent a detailed analysis of thyrse architecture development at the meristematic and molecular levels.

Phyllotaxis without symmetry: what can we learn from flower heads?

Phyllotaxis is commonly considered in the context of circular meristems or receptacles, yet non-circular (fasciated) structures also give rise to new primordia and organs. Here we investigate phyllotactic patterns in fasciated flower heads in the Asteraceae plant family. We begin by surveying the phenomenon of fasciation. We then show that phyllotactic patterns in fasciated heads can be generated by removing the inessential assumption of circularity from the previously published model of gerbera heads. To characterize these patterns, we revisit the conceptual framework in which phyllotactic patterns are commonly described. We note that some notions, in particular parastichies and parastichy numbers, maintain their significance in non-circular phyllotaxis, whereas others, in particular the divergence angle, need to be extended or lose their role. These observations highlight a number of open problems related to phyllotaxis in general, which may be elucidated by studies of fasciated heads.

Computational Models of Auxin-Driven Patterning in Shoots.

Auxin regulates many aspects of plant development and behavior, including the initiation of new outgrowth, patterning of vascular systems, control of branching, and responses to the environment. Computational models have complemented experimental studies of these processes. We review these models from two perspectives. First, we consider cellular and tissue-level models of interaction between auxin and its transporters in shoots. These models form a coherent body of results exploring different hypotheses pertinent to the patterning of new outgrowth and vascular strands. Second, we consider models operating at the level of plant organs and entire plants. We highlight techniques used to reduce the complexity of these models, which provide a path to capturing the essence of studied phenomena while running simulations efficiently.

Phyllotactic patterning of gerbera flower heads.

Phyllotaxis, the distribution of organs such as leaves and flowers on their support, is a key attribute of plant architecture. The geometric regularity of phyllotaxis has attracted multidisciplinary interest for centuries, resulting in an understanding of the patterns in the model plants Arabidopsis and tomato down to the molecular level. Nevertheless, the iconic example of phyllotaxis, the arrangement of individual florets into spirals in the heads of the daisy family of plants (Asteraceae), has not been fully explained. We integrate experimental data and computational models to explain phyllotaxis in Gerbera hybrida We show that phyllotactic patterning in gerbera is governed by changes in the size of the morphogenetically active zone coordinated with the growth of the head. The dynamics of these changes divides the patterning process into three phases: the development of an approximately circular pattern with a Fibonacci number of primordia near the head rim, its gradual transition to a zigzag pattern, and the development of a spiral pattern that fills the head on the template of this zigzag pattern. Fibonacci spiral numbers arise due to the intercalary insertion and lateral displacement of incipient primordia in the first phase. Our results demonstrate the essential role of the growth and active zone dynamics in the patterning of flower heads.

Latent Space Phenotyping: Automatic Image-Based Phenotyping for Treatment Studies.

Association mapping studies have enabled researchers to identify candidate loci for many important environmental tolerance factors, including agronomically relevant tolerance traits in plants. However, traditional genome-by-environment studies such as these require a phenotyping pipeline which is capable of accurately measuring stress responses, typically in an automated high-throughput context using image processing. In this work, we present Latent Space Phenotyping (LSP), a novel phenotyping method which is able to automatically detect and quantify response-to-treatment directly from images. We demonstrate example applications using data from an interspecific cross of the model C4 grass Setaria, a diversity panel of sorghum (S. bicolor), and the founder panel for a nested association mapping population of canola (Brassica napus L.). Using two synthetically generated image datasets, we then show that LSP is able to successfully recover the simulated QTL in both simple and complex synthetic imagery. We propose LSP as an alternative to traditional image analysis methods for phenotyping, enabling the phenotyping of arbitrary and potentially complex response traits without the need for engineering-complicated image-processing pipelines.

Phyllotaxis: is the golden angle optimal for light capture?

Phyllotactic patterns are some of the most conspicuous in nature. To create these patterns plants must control the divergence angle between the appearance of successive organs, sometimes to within a fraction of a degree. The most common angle is the Fibonacci or golden angle, and its prevalence has led to the hypothesis that it has been selected by evolution as optimal for plants with respect to some fitness benefits, such as light capture. We explore arguments for and against this idea with computer models. We have used both idealized and scanned leaves from Arabidopsis thaliana and Cardamine hirsuta to measure the overlapping leaf area of simulated plants after varying parameters such as leaf shape, incident light angles, and other leaf traits. We find that other angles generated by Fibonacci-like series found in nature are equally optimal for light capture, and therefore should be under similar evolutionary pressure. Our findings suggest that the iterative mechanism for organ positioning itself is a more likely target for evolutionary pressure, rather than a specific divergence angle, and our model demonstrates that the heteroblastic progression of leaf shape in A. thaliana can provide a potential fitness benefit via light capture.

Skin Patterning in Psoriasis by Spatial Interactions between Pathogenic Cytokines.

Disorders of human skin manifest themselves with patterns of lesions ranging from simple scattered spots to complex rings and spirals. These patterns are an essential characteristic of skin disease, yet the mechanisms through which they arise remain unknown. Here we show that all known patterns of psoriasis, a common inflammatory skin disease, can be explained in terms of reaction-diffusion. We constructed a computational model based on the known interactions between the main pathogenic cytokines: interleukins IL-17 and IL-23, and tumor necrosis factor TNF-α. Simulations revealed that the parameter space of the model contained all classes of psoriatic lesion patterns. They also faithfully reproduced the growth and evolution of the plaques and the response to treatment by cytokine targeting. Thus the pathogenesis of inflammatory diseases, such as psoriasis, may be readily understood in the framework of the stimulatory and inhibitory interactions between a few diffusing mediators.

Why plants make puzzle cells, and how their shape emerges.

The shape and function of plant cells are often highly interdependent. The puzzle-shaped cells that appear in the epidermis of many plants are a striking example of a complex cell shape, however their functional benefit has remained elusive. We propose that these intricate forms provide an effective strategy to reduce mechanical stress in the cell wall of the epidermis. When tissue-level growth is isotropic, we hypothesize that lobes emerge at the cellular level to prevent formation of large isodiametric cells that would bulge under the stress produced by turgor pressure. Data from various plant organs and species support the relationship between lobes and growth isotropy, which we test with mutants where growth direction is perturbed. Using simulation models we show that a mechanism actively regulating cellular stress plausibly reproduces the development of epidermal cell shape. Together, our results suggest that mechanical stress is a key driver of cell-shape morphogenesis.

The use of plant models in deep learning: an application to leaf counting in rosette plants.

Deep learning presents many opportunities for image-based plant phenotyping. Here we consider the capability of deep convolutional neural networks to perform the leaf counting task. Deep learning techniques typically require large and diverse datasets to learn generalizable models without providing a priori an engineered algorithm for performing the task. This requirement is challenging, however, for applications in the plant phenotyping field, where available datasets are often small and the costs associated with generating new data are high. In this work we propose a new method for augmenting plant phenotyping datasets using rendered images of synthetic plants. We demonstrate that the use of high-quality 3D synthetic plants to augment a dataset can improve performance on the leaf counting task. We also show that the ability of the model to generate an arbitrary distribution of phenotypes mitigates the problem of dataset shift when training and testing on different datasets. Finally, we show that real and synthetic plants are significantly interchangeable when training a neural network on the leaf counting task.

Crops In Silico: Generating Virtual Crops Using an Integrative and Multi-scale Modeling Platform.

Multi-scale models can facilitate whole plant simulations by linking gene networks, protein synthesis, metabolic pathways, physiology, and growth. Whole plant models can be further integrated with ecosystem, weather, and climate models to predict how various interactions respond to environmental perturbations. These models have the potential to fill in missing mechanistic details and generate new hypotheses to prioritize directed engineering efforts. Outcomes will potentially accelerate improvement of crop yield, sustainability, and increase future food security. It is time for a paradigm shift in plant modeling, from largely isolated efforts to a connected community that takes advantage of advances in high performance computing and mechanistic understanding of plant processes. Tools for guiding future crop breeding and engineering, understanding the implications of discoveries at the molecular level for whole plant behavior, and improved prediction of plant and ecosystem responses to the environment are urgently needed. The purpose of this perspective is to introduce Crops in silico (cropsinsilico.org), an integrative and multi-scale modeling platform, as one solution that combines isolated modeling efforts toward the generation of virtual crops, which is open and accessible to the entire plant biology community. The major challenges involved both in the development and deployment of a shared, multi-scale modeling platform, which are summarized in this prospectus, were recently identified during the first Crops in silico Symposium and Workshop.

A common developmental program can produce diverse leaf shapes.

Eudicot leaves have astoundingly diverse shapes. The central problem addressed in this paper is the developmental origin of this diversity. To investigate this problem, we propose a computational model of leaf development that generalizes the largely conserved molecular program for the reference plants Arabidopsis thaliana, Cardamine hirsuta and Solanum lycopersicum. The model characterizes leaf development as a product of three interwoven processes: the patterning of serrations, lobes and/or leaflets on the leaf margin; the patterning of the vascular system; and the growth of the leaf blade spanning the main veins. The veins play a significant morphogenetic role as a local determinant of growth directions. We show that small variations of this model can produce diverse leaf shapes, from simple to lobed to compound. It is thus plausible that diverse shapes of eudicot leaves result from small variations of a common developmental program.

A Model of Differential Growth-Guided Apical Hook Formation in Plants.

Differential cell growth enables flexible organ bending in the presence of environmental signals such as light or gravity. A prominent example of the developmental processes based on differential cell growth is the formation of the apical hook that protects the fragile shoot apical meristem when it breaks through the soil during germination. Here, we combined in silico and in vivo approaches to identify a minimal mechanism producing auxin gradient-guided differential growth during the establishment of the apical hook in the model plant Arabidopsis thaliana Computer simulation models based on experimental data demonstrate that asymmetric expression of the PIN-FORMED auxin efflux carrier at the concave (inner) versus convex (outer) side of the hook suffices to establish an auxin maximum in the epidermis at the concave side of the apical hook. Furthermore, we propose a mechanism that translates this maximum into differential growth, and thus curvature, of the apical hook. Through a combination of experimental and in silico computational approaches, we have identified the individual contributions of differential cell elongation and proliferation to defining the apical hook and reveal the role of auxin-ethylene crosstalk in balancing these two processes.

Auxin-driven patterning with unidirectional fluxes.

The plant hormone auxin plays an essential role in the patterning of plant structures. Biological hypotheses supported by computational models suggest that auxin may fulfil this role by regulating its own transport, but the plausibility of previously proposed models has been questioned. We applied the notion of unidirectional fluxes and the formalism of Petri nets to show that the key modes of auxin-driven patterning-the formation of convergence points and the formation of canals-can be implemented by biochemically plausible networks, with the fluxes measured by dedicated tally molecules or by efflux and influx carriers themselves. Common elements of these networks include a positive feedback of auxin efflux on the allocation of membrane-bound auxin efflux carriers (PIN proteins), and a modulation of this allocation by auxin in the extracellular space. Auxin concentration in the extracellular space is the only information exchanged by the cells. Canalization patterns are produced when auxin efflux and influx act antagonistically: an increase in auxin influx or concentration in the extracellular space decreases the abundance of efflux carriers in the adjacent segment of the membrane. In contrast, convergence points emerge in networks in which auxin efflux and influx act synergistically. A change in a single reaction rate may result in a dynamic switch between these modes, suggesting plausible molecular implementations of coordinated patterning of organ initials and vascular strands predicted by the dual polarization theory.

Modelling biomechanics of bark patterning in grasstrees.

BACKGROUND AND AIMS: Bark patterns are a visually important characteristic of trees, typically attributed to fractures occurring during secondary growth of the trunk and branches. An understanding of bark pattern formation has been hampered by insufficient information regarding the biomechanical properties of bark and the corresponding difficulties in faithfully modelling bark fractures using continuum mechanics. This study focuses on the genus Xanthorrhoea (grasstrees), which have an unusual bark-like structure composed of distinct leaf bases connected by sticky resin. Due to its discrete character, this structure is well suited for computational studies. METHODS: A dynamic computational model of grasstree development was created. The model captures both the phyllotactic pattern of leaf bases during primary growth and the changes in the trunk's width during secondary growth. A biomechanical representation based on a system of masses connected by springs is used for the surface of the trunk, permitting the emergence of fractures during secondary growth to be simulated. The resulting fracture patterns were analysed statistically and compared with images of real trees. KEY RESULTS: The model reproduces key features of grasstree bark patterns, including their variability, spanning elongated and reticulate forms. The patterns produced by the model have the same statistical character as those seen in real trees. CONCLUSIONS: The model was able to support the general hypothesis that the patterns observed in the grasstree bark-like layer may be explained in terms of mechanical fractures driven by secondary growth. Although the generality of the results is limited by the unusual structure of grasstree bark, it supports the hypothesis that bark pattern formation is primarily a biomechanical phenomenon.

Genetic control of plant development by overriding a geometric division rule.

Formative cell divisions are critical for multicellular patterning. In the early plant embryo, such divisions follow from orienting the division plane. A major unanswered question is how division plane orientation is genetically controlled, and in particular whether this relates to cell geometry. We have generated a complete 4D map of early Arabidopsis embryogenesis and used computational analysis to demonstrate that several divisions follow a rule that uses the smallest wall area going through the center of the cell. In other cases, however, cell division clearly deviates from this rule, which invariably leads to asymmetric cell division. By analyzing mutant embryos and through targeted genetic perturbation, we show that response to the hormone auxin triggers a deviation from the "shortest wall" rule. Our work demonstrates that a simple default rule couples division orientation to cell geometry in the embryo and that genetic regulation can create patterns by overriding the default rule.

A division in PIN-mediated auxin patterning during organ initiation in grasses.

The hormone auxin plays a crucial role in plant morphogenesis. In the shoot apical meristem, the PIN-FORMED1 (PIN1) efflux carrier concentrates auxin into local maxima in the epidermis, which position incipient leaf or floral primordia. From these maxima, PIN1 transports auxin into internal tissues along emergent paths that pattern leaf and stem vasculature. In Arabidopsis thaliana, these functions are attributed to a single PIN1 protein. Using phylogenetic and gene synteny analysis we identified an angiosperm PIN clade sister to PIN1, here termed Sister-of-PIN1 (SoPIN1), which is present in all sampled angiosperms except for Brassicaceae, including Arabidopsis. Additionally, we identified a conserved duplication of PIN1 in the grasses: PIN1a and PIN1b. In Brachypodium distachyon, SoPIN1 is highly expressed in the epidermis and is consistently polarized toward regions of high expression of the DR5 auxin-signaling reporter, which suggests that SoPIN1 functions in the localization of new primordia. In contrast, PIN1a and PIN1b are highly expressed in internal tissues, suggesting a role in vascular patterning. PIN1b is expressed in broad regions spanning the space between new primordia and previously formed vasculature, suggesting a role in connecting new organs to auxin sinks in the older tissues. Within these regions, PIN1a forms narrow canals that likely pattern future veins. Using a computer model, we reproduced the observed spatio-temporal expression and localization patterns of these proteins by assuming that SoPIN1 is polarized up the auxin gradient, and PIN1a and PIN1b are polarized to different degrees with the auxin flux. Our results suggest that examination and modeling of PIN dynamics in plants outside of Brassicaceae will offer insights into auxin-driven patterning obscured by the loss of the SoPIN1 clade in Brassicaceae.

An intracellular partitioning-based framework for tissue cell polarity in plants and animals.

Tissue cell polarity plays a major role in plant and animal development. We propose that a fundamental building block for tissue cell polarity is the process of intracellular partitioning, which can establish individual cell polarity in the absence of asymmetric cues. Coordination of polarities may then arise through cell-cell coupling, which can operate directly, through membrane-spanning complexes, or indirectly, through diffusible molecules. Polarity is anchored to tissues through organisers located at boundaries. We show how this intracellular partitioning-based framework can be applied to both plant and animal systems, allowing different processes to be placed in a common evolutionary and mechanistic context.

The interplay between inflorescence development and function as the crucible of architectural diversity.

BACKGROUND: Most angiosperms present flowers in inflorescences, which play roles in reproduction, primarily related to pollination, beyond those served by individual flowers alone. An inflorescence's overall reproductive contribution depends primarily on the three-dimensional arrangement of the floral canopy and its dynamics during its flowering period. These features depend in turn on characteristics of the underlying branching structure (scaffold) that supports and supplies water and nutrients to the floral canopy. This scaffold is produced by developmental algorithms that are genetically specified and hormonally mediated. Thus, the extensive inflorescence diversity evident among angiosperms evolves through changes in the developmental programmes that specify scaffold characteristics, which in turn modify canopy features that promote reproductive performance in a particular pollination and mating environment. Nevertheless, developmental and ecological aspects of inflorescences have typically been studied independently, limiting comprehensive understanding of the relations between inflorescence form, reproductive function, and evolution. SCOPE: This review fosters an integrated perspective on inflorescences by summarizing aspects of their development and pollination function that enable and guide inflorescence evolution and diversification. CONCLUSIONS: The architecture of flowering inflorescences comprises three related components: topology (branching patterns, flower number), geometry (phyllotaxis, internode and pedicel lengths, three-dimensional flower arrangement) and phenology (flower opening rate and longevity, dichogamy). Genetic and developmental evidence reveals that these components are largely subject to quantitative control. Consequently, inflorescence evolution proceeds along a multidimensional continuum. Nevertheless, some combinations of topology, geometry and phenology are represented more commonly than others, because they serve reproductive function particularly effectively. For wind-pollinated species, these combinations often represent compromise solutions to the conflicting physical influences on pollen removal, transport and deposition. For animal-pollinated species, dominant selective influences include the conflicting benefits of large displays for attracting pollinators and of small displays for limiting among-flower self-pollination. The variety of architectural components that comprise inflorescences enable diverse resolutions of these conflicts.