Fluctuations of the transcription factor ATML1 generate the pattern of giant cells in the Arabidopsis sepal

Self-organization underlies developmental robustness in plants

Development is a self-organized process that builds on cells and their interactions. Cells are heterogeneous in gene expression, growth, and division; yet how development is robust despite such heterogeneity is a fascinating question. Here, we review recent progress on this topic, highlighting how developmental robustness is achieved through self-organization. We will first discuss sources of heterogeneity, including stochastic gene expression, heterogeneity in growth rate and direction, and heterogeneity in division rate and precision. We then discuss cellular mechanisms that buffer against such noise, including Paf1C- and miRNA-mediated denoising, spatiotemporal growth averaging and compensation, mechanisms to improve cell division precision, and coordination of growth rate and developmental timing between different parts of an organ. We also discuss cases where such heterogeneity is not buffered but utilized for development. Finally, we highlight potential directions for future studies of noise and developmental robustness.

OsMADS6-OsMADS32 and REP1 control palea cellular heterogeneity and morphogenesis in rice

Precise regulation of cell proliferation and differentiation is vital for organ morphology. Rice palea, serving as sepal, comprises two distinct regions: the marginal region (MRP) and body of palea (BOP), housing heterogeneous cell populations, which makes it an ideal system for studying organ morphogenesis. We report that the transcription factor (TF) REP1 promotes epidermal cell proliferation and differentiation in the BOP, resulting in hard silicified protrusion cells, by regulating the cyclin-dependent kinase gene, OsCDKB1;1. Conversely, TFs OsMADS6 and OsMADS32 are expressed exclusively in the MRP, where they limit cell division rates by inhibiting OsCDKB2;1 expression and promote endoreduplication, yielding elongated epidermal cells. Furthermore, reciprocal inhibition between the OsMADS6-OsMADS32 complex and REP1 fine-tunes the balance between cell division and differentiation during palea morphogenesis. We further show the functional conservation of these organ identity genes in heterogeneous cell growth in Arabidopsis, emphasizing a critical framework for controlling cellular heterogeneity in organ morphogenesis.

A gradient of the HD-Zip regulator Woolly regulates multicellular trichome morphogenesis in tomato

Homeodomain (HD) proteins regulate embryogenesis in animals such as the fruit fly (Drosophila melanogaster), often in a concentration-dependent manner. HD-leucine zipper (Zip) IV family genes are unique to plants and often function in the L1 epidermal cell layer. However, our understanding of the roles of HD-Zip IV family genes in plant morphogenesis is limited. In this study, we investigated the morphogenesis of tomato (Solanum lycopersicum) multicellular trichomes, a type of micro-organ in plants. We found that a gradient of the HD-Zip IV regulator Woolly (Wo) coordinates spatially polarized cell division and cell expansion in multicellular trichomes. Moreover, we identified a TEOSINTE BRANCHED1, CYCLOIDEA, and PROLIFERATING CELL NUCLEAR ANTIGEN BINDING FACTOR (TCP) transcription factor-encoding gene, SlBRANCHED2a (SlBRC2a), as a key downstream target of Wo that regulates the transition from cell division to cell expansion. High levels of Wo promote cell division in apical trichome cells, whereas in basal trichome cells, Wo mediates a negative feedback loop with SlBRC2a that forces basal cells to enter endoreduplication. The restricted high and low activities of Wo pattern the morphogenesis of tomato multicellular trichomes. These findings provide insights into the functions of HD-Zip IV genes during plant morphogenesis.

Genome-wide characterisation of HD-Zip transcription factors and functional analysis of PbHB24 during stone cell formation in Chinese white pear (Pyrus bretschneideri)

BACKGROUND

The homodomain-leucine zipper (HD-Zip) is a conserved transcription factor family unique to plants that regulate multiple developmental processes including lignificaion. Stone cell content is a key determinant negatively affecting pear fruit quality, which causes a grainy texture of fruit flesh, because of the lignified cell walls.

RESULTS

In this study, a comprehensive bioinformatics analysis of HD-Zip genes in Chinese white pear (Pyrus bretschneideri) (PbHBs) was performed. Genome-wide identification of the PbHB gene family revealed 67 genes encoding PbHB proteins, which could be divided into four subgroups (I, II, III, and IV). For some members, similar intron/exon structural patterns support close evolutionary relationships within the same subgroup. The functions of each subgroup of the PbHB family were predicted through comparative analysis with the HB genes in Arabidopsis and other plants. Cis-element analysis indicated that PbHB genes might be involved in plant hormone signalling and external environmental responses, such as light, stress, and temperature. Furthermore, RNA-sequencing data and quantitative real-time PCR (RT-qPCR) verification revealed the regulatory roles of PbHB genes in pear stone cell formation. Further, co-expression network analysis revealed that the eight PbHB genes could be classified into different clusters of co-expression with lignin-related genes. Besides, the biological function of PbHB24 in promoting stone cell formation has been demonstrated by overexpression in fruitlets.

CONCLUSIONS

This study provided the comprehensive analysis of PbHBs and highlighted the importance of PbHB24 during stone cell development in pear fruits.

Sepal shape variability is robust to cell size heterogeneity in Arabidopsis

How organisms produce organs with robust shapes and sizes is still an open question. In recent years, the Arabidopsis sepal has been used as a model system to study this question because of its highly reproducible shape and size. One interesting aspect of the sepal is that its epidermis contains cells of very different sizes. Previous reports have qualitatively shown that sepals with more or less giant cells exhibit comparable final size and shape. Here, we investigate this question using quantitative approaches. We find that a mixed population of cell size modestly contribute to the normal width of the sepal but is not essential for its shape robustness. Furthermore, in a mutant with increased cell and organ growth variability, the change in final sepal shape caused by giant cells is exaggerated but the shape robustness is not affected. This formally demonstrates that sepal shape variability is robust to cell size heterogeneity.

HD-Zip proteins modify floral structures for self-pollination in tomato

Cleistogamy is a type of self-pollination that relies on the formation of a stigma-enclosing floral structure. We identify three homeodomain-leucine zipper IV (HD-Zip IV) genes that coordinately promote the formation of interlocking trichomes at the anther margin to unite neighboring anthers, generating a closed anther cone and cleistogamy (flower morphology necessitating strict self-pollination). These HD-Zip IV genes also control style length by regulating the transition from cell division to endoreduplication. The expression of these HD-Zip IV genes and their downstream gene, Style 2.1, was sequentially modified to shape the cleistogamy morphology during tomato evolution and domestication. Our results provide insights into the molecular basis of cleistogamy in modern tomato and suggest targets for improving fruit set and preventing pollen contamination in genetically modified crops.

SHR and SCR coordinate root patterning and growth early in the cell cycle

Precise control of cell division is essential for proper patterning and growth during the development of multicellular organisms. Coordination of formative divisions that generate new tissue patterns with proliferative divisions that promote growth is poorly understood. SHORTROOT (SHR) and SCARECROW (SCR) are transcription factors that are required for formative divisions in the stem cell niche of Arabidopsis roots1,2. Here we show that levels of SHR and SCR early in the cell cycle determine the orientation of the division plane, resulting in either formative or proliferative cell division. We used 4D quantitative, long-term and frequent (every 15 min for up to 48 h) light sheet and confocal microscopy to probe the dynamics of SHR and SCR in tandem within single cells of living roots. Directly controlling their dynamics with an SHR induction system enabled us to challenge an existing bistable model3 of the SHR-SCR gene-regulatory network and to identify key features that are essential for rescue of formative divisions in shr mutants. SHR and SCR kinetics do not align with the expected behaviour of a bistable system, and only low transient levels, present early in the cell cycle, are required for formative divisions. These results reveal an uncharacterized mechanism by which developmental regulators directly coordinate patterning and growth.

HD-Zip IV transcription factors: Drivers of epidermal cell fate integrate metabolic signals

The leaf epidermis comprises the outermost layer of cells that protect plants against environmental stresses such as drought, ultraviolet radiation, and pathogen attack. Research over the past decades highlights the role of class IV homeodomain leucine-zipper (HD-Zip IV) transcription factors (TFs) in driving differentiation of various epidermal cell types, such as trichomes, guard cells, and pavement cells. Evolutionary origins of this family in the charophycean green algae and HD-Zip-specific gene expression in the maternal genome provide clues to unlocking their secrets which include ties to cell cycle regulation. A distinguishing feature of these TFs is the presence of a lipid binding pocket that integrates metabolic information with gene expression. Identities of metabolic partners are beginning to emerge, uncovering feedback loops to maintain epidermal cell specification. Discoveries of associated molecular mechanisms are revealing fascinating links to phospholipid and sphingolipid metabolism and mechanical signaling.

Giant cells: multiple cells unite to survive

Multinucleated Giant Cells (MGCs) are specialized cells that develop from the fusion of multiple cells, and their presence is commonly observed in human cells during various infections. However, MGC formation is not restricted to infections alone but can also occur through different mechanisms, such as endoreplication and abortive cell cycle. These processes lead to the formation of polyploid cells, eventually resulting in the formation of MGCs. In Entamoeba, a protozoan parasite that causes amoebic dysentery and liver abscesses in humans, the formation of MGCs is a unique phenomenon and not been reported in any other protozoa. This organism is exposed to various hostile environmental conditions, including changes in temperature, pH, and nutrient availability, which can lead to stress and damage to its cells. The formation of MGCs in Entamoeba is thought to be a survival strategy to cope with these adverse conditions. This organism forms MGCs through cell aggregation and fusion in response to osmotic and heat stress. The MGCs in Entamoeba are thought to have increased resistance to various stresses and can survive longer than normal cells under adverse conditions. This increased survival could be due to the presence of multiple nuclei, which could provide redundancy in case of DNA damage or mutations. Additionally, MGCs may play a role in the virulence of Entamoeba as they are found in the inflammatory foci of amoebic liver abscesses and other infections caused by Entamoeba. The presence of MGCs in these infections suggests that they may contribute to the pathogenesis of the disease. Overall, this article offers valuable insights into the intriguing phenomenon of MGC formation in Entamoeba. By unraveling the mechanisms behind this process and examining its implications, researchers can gain a deeper understanding of the complex biology of Entamoeba and potentially identify new targets for therapeutic interventions. The study of MGCs in Entamoeba serves as a gateway to exploring the broader field of cell fusion in various organisms, providing a foundation for future investigations into related cellular processes and their significance in health and disease.

Increased gene expression variability hinders the formation of regional mechanical conflicts leading to reduced organ shape robustness

To relate gene networks and organ shape, one needs to address two wicked problems: i) Gene expression is often variable locally, and shape is reproducible globally; ii) gene expression can have cascading effects on tissue mechanics, with possibly counterintuitive consequences for the final organ shape. Here, we address such wicked problems, taking advantage of simpler plant organ development where shape only emerges from cell division and elongation. We confirm that mutation in VERNALIZATION INDEPENDENCE 3 (VIP3), a subunit of the conserved polymerase-associated factor 1 complex (Paf1C), increases gene expression variability in Arabidopsis. Then, we focused on the Arabidopsis sepal, which exhibits a reproducible shape and stereotypical regional growth patterns. In vip3 sepals, we measured higher growth heterogeneity between adjacent cells. This even culminated in the presence of negatively growing cells in specific growth conditions. Interestingly, such increased local noise interfered with the stereotypical regional pattern of growth. We previously showed that regional differential growth at the wild-type sepal tip triggers a mechanical conflict, to which cells resist by reinforcing their walls, leading to growth arrest. In vip3, the disturbed regional growth pattern delayed organ growth arrest and increased final organ shape variability. Altogether, we propose that gene expression variability is managed by Paf1C to ensure organ robustness by building up mechanical conflicts at the regional scale, instead of the local scale.

Spatial consistency of cell growth direction during organ morphogenesis requires CELLULOSE SYNTHASE INTERACTIVE1

Extracellular matrices contain fibril-like polymers often organized in parallel arrays. Although their role in morphogenesis has been long recognized, it remains unclear how the subcellular control of fibril synthesis translates into organ shape. We address this question using the Arabidopsis sepal as a model organ. In plants, cell growth is restrained by the cell wall (extracellular matrix). Cellulose microfibrils are the main load-bearing wall component, thought to channel growth perpendicularly to their main orientation. Given the key function of CELLULOSE SYNTHASE INTERACTIVE1 (CSI1) in guidance of cellulose synthesis, we investigate the role of CSI1 in sepal morphogenesis. We observe that sepals from csi1 mutants are shorter, although their newest cellulose microfibrils are more aligned compared to wild-type. Surprisingly, cell growth anisotropy is similar in csi1 and wild-type plants. We resolve this apparent paradox by showing that CSI1 is required for spatial consistency of growth direction across the sepal.

Enhancer activation via TCP and HD-ZIP and repression by Dof transcription factors mediate giant cell-specific expression

Proper cell-type identity relies on highly coordinated regulation of gene expression. Regulatory elements such as enhancers can produce cell type-specific expression patterns, but the mechanisms underlying specificity are not well understood. We previously identified an enhancer region capable of driving specific expression in giant cells, which are large, highly endoreduplicated cells in the Arabidopsis thaliana sepal epidermis. In this study, we use the giant cell enhancer as a model to understand the regulatory logic that promotes cell type-specific expression. Our dissection of the enhancer revealed that giant cell specificity is mediated primarily through the combination of two activators and one repressor. HD-ZIP and TCP transcription factors are involved in the activation of expression throughout the epidermis. High expression of HD-ZIP transcription factor genes in giant cells promoted higher expression driven by the enhancer in giant cells. Dof transcription factors repressed the activity of the enhancer such that only giant cells maintained enhancer activity. Thus, our data are consistent with a conceptual model whereby cell type-specific expression emerges from the combined activities of three transcription factor families activating and repressing expression in epidermal cells.

Anisotropic cell growth at the leaf base promotes age-related changes in leaf shape in Arabidopsis thaliana

Plants undergo extended morphogenesis. The shoot apical meristem (SAM) allows for reiterative development and the formation of new structures throughout the life of the plant. Intriguingly, the SAM produces morphologically different leaves in an age-dependent manner, a phenomenon known as heteroblasty. In Arabidopsis thaliana, the SAM produces small orbicular leaves in the juvenile phase, but gives rise to large elliptical leaves in the adult phase. Previous studies have established that a developmental decline of microRNA156 (miR156) is necessary and sufficient to trigger this leaf shape switch, although the underlying mechanism is poorly understood. Here we show that the gradual increase in miR156-targeted SQUAMOSA PROMOTER BINDING PROTEIN-LIKE transcription factors with age promotes cell growth anisotropy in the abaxial epidermis at the base of the leaf blade, evident by the formation of elongated giant cells. Time-lapse imaging and developmental genetics further revealed that the establishment of adult leaf shape is tightly associated with the longitudinal cell expansion of giant cells, accompanied by a prolonged cell proliferation phase in their vicinity. Our results thus provide a plausible cellular mechanism for heteroblasty in Arabidopsis, and contribute to our understanding of anisotropic growth in plants.

Cell Cycle Dynamics during Stomatal Development: Window of MUTE Action and Ramification of Its Loss-of-Function on an Uncommitted Precursor

Plants develop in the absence of cell migration. As such, cell division and differentiation need to be coordinated for functional tissue formation. Cellular valves on the plant epidermis, stomata, are generated through a stereotypical sequence of cell division and differentiation events. In Arabidopsis, three master regulatory transcription factors, SPEECHLESS (SPCH), MUTE and FAMA, sequentially drive initiation, proliferation and differentiation of stomata. Among them, MUTE switches the cell cycle mode from proliferative asymmetric division to terminal symmetric division and orchestrates the execution of the single symmetric division event. However, it remains unclear to what extent MUTE regulates the expression of cell cycle genes through the symmetric division and whether MUTE accumulation itself is gated by the cell cycle. Here, we show that MUTE directly upregulates the expression of cell cycle components throughout the terminal cell cycle phases of a stomatal precursor, not only core cell cycle engines but also check-point regulators. Time-lapse live imaging using the multicolor Plant Cell Cycle Indicator revealed that MUTE accumulates up to the early G2 phase, whereas its successor and direct target, FAMA, accumulate at late G2 through terminal mitosis. In the absence of MUTE, meristemoids fail to differentiate and their G1 phase elongates as they reiterate asymmetric divisions. Together, our work provides the framework of cell cycle and master regulatory transcription factors to coordinate a single symmetric cell division and suggests a mechanism for the eventual cell cycle arrest of an uncommitted stem-cell-like precursor at the G1 phase.

A HD-ZIP transcription factor specifies fates of multicellular trichomes via dosage-dependent mechanisms in tomato

Hair-like structures are shared by most living organisms. The hairs on plant surfaces, commonly referred to as trichomes, form diverse types to sense and protect against various stresses. However, it is unclear how trichomes differentiate into highly variable forms. Here, we show that a homeodomain leucine zipper (HD-ZIP) transcription factor named Woolly controls the fates of distinct trichomes in tomato via a dosage-dependent mechanism. The autocatalytic reinforcement of Woolly is counteracted by an autoregulatory negative feedback loop, creating a circuit with a high or low Woolly level. This biases the transcriptional activation of separate antagonistic cascades that lead to different trichome types. Our results identify the developmental switch of trichome formation and provide mechanistic insights into the progressive fate specification in plants, as well as a path to enhancing plant stress resistance and the production of beneficial chemicals.

Large-scale analysis and computer modeling reveal hidden regularities behind variability of cell division patterns in Arabidopsis thaliana embryogenesis

Noise plays a major role in cellular processes and in the development of tissues and organs. Several studies have examined the origin, the integration or the accommodation of noise in gene expression, cell growth and elaboration of organ shape. By contrast, much less is known about variability in cell division plane positioning, its origin and links with cell geometry, and its impact on tissue organization. Taking advantage of the first-stereotyped-then-variable division patterns in the embryo of the model plant Arabidopsis thaliana, we combined 3D imaging and quantitative cell shape and cell lineage analysis together with mathematical and computer modeling to perform a large-scale, systematic analysis of variability in division plane orientation. Our results reveal that, paradoxically, variability in cell division patterns of Arabidopsis embryos is accompanied by a progressive reduction of heterogeneity in cell shape topology. The paradox is solved by showing that variability operates within a reduced repertoire of possible division plane orientations that is related to cell geometry. We show that in several domains of the embryo, a recently proposed geometrical division rule recapitulates observed variable patterns, suggesting that variable patterns emerge from deterministic principles operating in a variable geometrical context. Our work highlights the importance of emerging patterns in the plant embryo under iterated division principles, but also reveal domains where deviations between rule predictions and experimental observations point to additional regulatory mechanisms.

Meeting report - Cell size and growth: from single cells to the tree of life

In April 2022, The Company of Biologists hosted their first post-pandemic in-person Workshop at Buxted Park Country House in the Sussex countryside. The Workshop, entitled 'Cell size and growth: from single cells to the tree of life', gathered a small group of early-career and senior researchers with expertise in cell size spanning a broad range of organisms, including bacteria, yeast, animal cells, embryos and plants, and working in fields from cell biology to ecology and evolutionary biology. The programme made ample room for fruitful discussions and provided a much-needed opportunity to discuss the most recent findings relating to the regulation of cell size and growth, identify the emerging challenges for the field, and build a community after the pandemic.

Context-specific functions of transcription factors controlling plant development: From leaves to flowers

Plant development is regulated by transcription factors that often act in more than one process and stage of development. Yet the molecular mechanisms that govern the functional diversity and specificity of these proteins remains far from understood. Flower development provides an ideal context to study these mechanisms since the development of distinct floral organs depends on similar but distinct combinations of transcriptional regulators. Recent work also highlights the importance of leaf polarity regulators as additional key factors in flower initiation, floral organ morphogenesis, and possibly floral organ positioning. A detailed understanding of how these factors work in combination will enable us to address outstanding questions in flower development including how distinct shapes and positions of floral organs are generated. Experimental approaches and computer-based modeling will be required to characterize gene-regulatory networks at the level of single cells.

Cell type-specific dynamics underlie cellular growth variability in plants

Coordination of growth, patterning and differentiation is required for shaping organs in multicellular organisms. In plants, cell growth is controlled by positional information, yet the behavior of individual cells is often highly heterogeneous. The origin of this variability is still unclear. Using time-lapse imaging, we determined the source and relevance of cellular growth variability in developing organs of Arabidopsis thaliana. We show that growth is more heterogeneous in the leaf blade than in the midrib and petiole, correlating with higher local differences in growth rates between neighboring cells in the blade. This local growth variability coincides with developing stomata. Stomatal lineages follow a specific, time-dependent growth program that is different from that of their surroundings. Quantification of cellular dynamics in the leaves of a mutant lacking stomata, as well as analysis of floral organs, supports the idea that growth variability is mainly driven by stomata differentiation. Thus, the cell-autonomous behavior of specialized cells is the main source of local growth variability in otherwise homogeneously growing tissue. Those growth differences are buffered by the immediate neighbors of stomata and trichomes to achieve robust organ shapes.

A hierarchical transcriptional network activates specific CDK inhibitors that regulate G2 to control cell size and number in Arabidopsis

How cell size and number are determined during organ development remains a fundamental question in cell biology. Here, we identified a GRAS family transcription factor, called SCARECROW-LIKE28 (SCL28), with a critical role in determining cell size in Arabidopsis. SCL28 is part of a transcriptional regulatory network downstream of the central MYB3Rs that regulate G2 to M phase cell cycle transition. We show that SCL28 forms a dimer with the AP2-type transcription factor, AtSMOS1, which defines the specificity for promoter binding and directly activates transcription of a specific set of SIAMESE-RELATED (SMR) family genes, encoding plant-specific inhibitors of cyclin-dependent kinases and thus inhibiting cell cycle progression at G2 and promoting the onset of endoreplication. Through this dose-dependent regulation of SMR transcription, SCL28 quantitatively sets the balance between cell size and number without dramatically changing final organ size. We propose that this hierarchical transcriptional network constitutes a cell cycle regulatory mechanism that allows to adjust cell size and number to attain robust organ growth.

Fifteen compelling open questions in plant cell biology

As scientists, we are at least as excited about the open questions-the things we do not know-as the discoveries. Here, we asked 15 experts to describe the most compelling open questions in plant cell biology. These are their questions: How are organelle identity, domains, and boundaries maintained under the continuous flux of vesicle trafficking and membrane remodeling? Is the plant cortical microtubule cytoskeleton a mechanosensory apparatus? How are the cellular pathways of cell wall synthesis, assembly, modification, and integrity sensing linked in plants? Why do plasmodesmata open and close? Is there retrograde signaling from vacuoles to the nucleus? How do root cells accommodate fungal endosymbionts? What is the role of cell edges in plant morphogenesis? How is the cell division site determined? What are the emergent effects of polyploidy on the biology of the cell, and how are any such "rules" conditioned by cell type? Can mechanical forces trigger new cell fates in plants? How does a single differentiated somatic cell reprogram and gain pluripotency? How does polarity develop de-novo in isolated plant cells? What is the spectrum of cellular functions for membraneless organelles and intrinsically disordered proteins? How do plants deal with internal noise? How does order emerge in cells and propagate to organs and organisms from complex dynamical processes? We hope you find the discussions of these questions thought provoking and inspiring.

Cell biology of the leaf epidermis: Fate specification, morphogenesis, and coordination

As the outermost layer of plants, the epidermis serves as a critical interface between plants and the environment. During leaf development, the differentiation of specialized epidermal cell types, including stomatal guard cells, pavement cells, and trichomes, occurs simultaneously, each providing unique and pivotal functions for plant growth and survival. Decades of molecular-genetic and physiological studies have unraveled key players and hormone signaling specifying epidermal differentiation. However, most studies focus on only one cell type at a time, and how these distinct cell types coordinate as a unit is far from well-comprehended. Here we provide a review on the current knowledge of regulatory mechanisms underpinning the fate specification, differentiation, morphogenesis, and positioning of these specialized cell types. Emphasis is given to their shared developmental origins, fate flexibility, as well as cell cycle and hormonal controls. Furthermore, we discuss computational modeling approaches to integrate how mechanical properties of individual epidermal cell types and entire tissue/organ properties mutually influence each other. We hope to illuminate the underlying mechanisms coordinating the cell differentiation that ultimately generate a functional leaf epidermis.

Using quantitative methods to understand leaf epidermal development

As the interface between plants and the environment, the leaf epidermis provides the first layer of protection against drought, ultraviolet light, and pathogen attack. This cell layer comprises highly coordinated and specialised cells such as stomata, pavement cells and trichomes. While much has been learned from the genetic dissection of stomatal, trichome and pavement cell formation, emerging methods in quantitative measurements that monitor cellular or tissue dynamics will allow us to further investigate cell state transitions and fate determination in leaf epidermal development. In this review, we introduce the formation of epidermal cell types in Arabidopsis and provide examples of quantitative tools to describe phenotypes in leaf research. We further focus on cellular factors involved in triggering cell fates and their quantitative measurements in mechanistic studies and biological patterning. A comprehensive understanding of how a functional leaf epidermis develops will advance the breeding of crops with improved stress tolerance.

Sculpting the surface: Structural patterning of plant epidermis

Plant epidermis are multifunctional surfaces that directly affect how plants interact with animals or microorganisms and influence their ability to harvest or protect from abiotic factors. To do this, plants rely on minuscule structures that confer remarkable properties to their outer layer. These microscopic features emerge from the hierarchical organization of epidermal cells with various shapes and dimensions combined with different elaborations of the cuticle, a protective film that covers plant surfaces. Understanding the properties and functions of those tridimensional elements as well as disentangling the mechanisms that control their formation and spatial distribution warrant a multidisciplinary approach. Here we show how interdisciplinary efforts of coupling modern tools of experimental biology, physics, and chemistry with advanced computational modeling and state-of-the art microscopy are yielding broad new insights into the seemingly arcane patterning processes that sculpt the outer layer of plants.

Quantitative imaging of RNA polymerase II activity in plants reveals the single-cell basis of tissue-wide transcriptional dynamics

The responses of plants to their environment are often dependent on the spatiotemporal dynamics of transcriptional regulation. While live-imaging tools have been used extensively to quantitatively capture rapid transcriptional dynamics in living animal cells, the lack of implementation of these technologies in plants has limited concomitant quantitative studies in this kingdom. Here, we applied the PP7 and MS2 RNA-labelling technologies for the quantitative imaging of RNA polymerase II activity dynamics in single cells of living plants as they respond to experimental treatments. Using this technology, we counted nascent RNA transcripts in real time in Nicotiana benthamiana (tobacco) and Arabidopsis thaliana. Examination of heat shock reporters revealed that plant tissues respond to external signals by modulating the proportion of cells that switch from an undetectable basal state to a high-transcription state, instead of modulating the rate of transcription across all cells in a graded fashion. This switch-like behaviour, combined with cell-to-cell variability in transcription rate, results in mRNA production variability spanning three orders of magnitude. We determined that cellular heterogeneity stems mainly from stochasticity intrinsic to individual alleles instead of variability in cellular composition. Together, our results demonstrate that it is now possible to quantitatively study the dynamics of transcriptional programs in single cells of living plants.

Cell identity specification in plants: lessons from flower development

Multicellular organisms display a fascinating complexity of cellular identities and patterns of diversification. The concept of 'cell type' aims to describe and categorize this complexity. In this review, we discuss the traditional concept of cell types and highlight the impact of single-cell technologies and spatial omics on the understanding of cellular differentiation in plants. We summarize and compare position-based and lineage-based mechanisms of cell identity specification using flower development as a model system. More than understanding ontogenetic origins of differentiated cells, an important question in plant science is to understand their position- and developmental stage-specific heterogeneity. Combinatorial action and crosstalk of external and internal signals is the key to cellular heterogeneity, often converging on transcription factors that orchestrate gene expression programs.

Modelling the spatial crosstalk between two biochemical signals explains wood formation dynamics and tree-ring structure

In conifers, xylogenesis during a growing season produces a very characteristic tree-ring structure: large, thin-walled earlywood cells followed by narrow, thick-walled latewood cells. Although many factors influence the dynamics of differentiation and the final dimensions of xylem cells, the associated patterns of variation remain very stable from one year to the next. While radial growth is characterized by an S-shaped curve, the widths of xylem differentiation zones exhibit characteristic skewed bell-shaped curves. These elements suggest a strong internal control of xylogenesis. It has long been hypothesized that much of this regulation relies on a morphogenetic gradient of auxin. However, recent modelling studies have shown that while this hypothesis could account for the dynamics of stem radial growth and the zonation of the developing xylem, it failed to reproduce the characteristic tree-ring structure. Here, we investigated the hypothesis of regulation by a crosstalk between auxin and a second biochemical signal, by using computational morphodynamics. We found that, in conifers, such a crosstalk is sufficient to simulate the characteristic features of wood formation dynamics, as well as the resulting tree-ring structure. In this model, auxin controls cell enlargement rates while another signal (e.g. cytokinin, tracheary element differentiation inhibitory factor) drives cell division and auxin polar transport.

A Quarter Century History of ATML1 Gene Research

The cloning of the ATML1 gene, encoding an HD-ZIP class IV transcription factor, was first reported in 1996. Because ATML1 mRNA was preferentially detected in the shoot epidermis, cis-regulatory sequences of ATML1 have been used to drive gene expression in the outermost cells of the shoot apical meristem and leaves, even before the function of ATML1 was understood. Later studies revealed that ATML1 is required for developmental processes related to shoot epidermal specification and differentiation. Consistent with its central role in epidermal development, ATML1 activity has been revealed to be restricted to the outermost cells via several regulatory mechanisms. In this review, we look back on the history of ATML1 research and provide a perspective for future studies.

Arabidopsis sepals: A model system for the emergent process of morphogenesis

During development, Arabidopsis thaliana sepal primordium cells grow, divide and interact with their neighbours, giving rise to a sepal with the correct size, shape and form. Arabidopsis sepals have proven to be a good system for elucidating the emergent processes driving morphogenesis due to their simplicity, their accessibility for imaging and manipulation, and their reproducible development. Sepals undergo a basipetal gradient of growth, with cessation of cell division, slow growth and maturation starting at the tip of the sepal and progressing to the base. In this review, I discuss five recent examples of processes during sepal morphogenesis that yield emergent properties: robust size, tapered tip shape, laminar shape, scattered giant cells and complex gene expression patterns. In each case, experiments examining the dynamics of sepal development led to the hypotheses of local rules. In each example, a computational model was used to demonstrate that these local rules are sufficient to give rise to the emergent properties of morphogenesis.

What is quantitative plant biology?

Quantitative plant biology is an interdisciplinary field that builds on a long history of biomathematics and biophysics. Today, thanks to high spatiotemporal resolution tools and computational modelling, it sets a new standard in plant science. Acquired data, whether molecular, geometric or mechanical, are quantified, statistically assessed and integrated at multiple scales and across fields. They feed testable predictions that, in turn, guide further experimental tests. Quantitative features such as variability, noise, robustness, delays or feedback loops are included to account for the inner dynamics of plants and their interactions with the environment. Here, we present the main features of this ongoing revolution, through new questions around signalling networks, tissue topology, shape plasticity, biomechanics, bioenergetics, ecology and engineering. In the end, quantitative plant biology allows us to question and better understand our interactions with plants. In turn, this field opens the door to transdisciplinary projects with the society, notably through citizen science.

Epidermal restriction confers robustness to organ shapes

The shape of comparable tissues and organs is consistent among individuals of a given species, but how this consistency or robustness is achieved remains an open question. The interaction between morphogenetic factors determines organ formation and subsequent shaping, which is ultimately a mechanical process. Using a computational approach, we show that the epidermal layer is essential for the robustness of organ geometry control. Specifically, proper epidermal restriction allows organ asymmetry maintenance, and the tensile epidermal layer is sufficient to suppress local variability in growth, leading to shape robustness. The model explains the enhanced organ shape variations in epidermal mutant plants. In addition, differences in the patterns of epidermal restriction may underlie the initial establishment of organ asymmetry. Our results show that epidermal restriction can answer the longstanding question of how cellular growth noise is averaged to produce precise organ shapes, and the findings also shed light on organ asymmetry establishment.

Transcriptomic Analysis and Specific Expression of Transcription Factor Genes in the Root and Sporophyll of Dryopteris fragrans (L.) Schott

Background: Dryopteris fragrans, which is densely covered with glandular trichomes, is considered to be one of the ferns with the most medicinal potential. The transcriptomes from selected tissues of D. fragrans were collected and analyzed for functional and comparative genomic studies. The aim of this study was to determine the transcriptomic characteristics of wild D. fragrans sporangium in tissues from the SR (root), SL (sporophyll), and TRL (sporophyll with glandular trichomes removed). Results: Cluster analysis identified genes that were highly expressed in an organ-specific manner according to read mapping, feature counting, and normalization. The functional map identified gene clusters that can uniquely describe the function of each tissue. We identified a group of three tissue-specific transcription factors targeting the SL, SR, and TRL. In addition, highly expressed transcription factors (TFs) were found in each tissue-specific gene cluster, where ERF and bHLH transcription factors were the two types showing the most distinct expression patterns between the three different tissues. The specific expression of transcription factor genes varied between the different types of tissues. The numbers of transcription factors specifically expressed in the roots and sporophylls were 60 and 30, respectively, while only seven were found for the sporophylls with glandular trichomes removed. The expression of genes known to be associated with the development of glandular trichomes in flowering plants, including MIXTA, ATML1, and MYB106, were also validated and are discussed. In particular, a unigene encoding MIXTA was identified and exhibited the highest expression level in SL in D. fragrans. Conclusions: This study is the first report of global transcriptomic analysis in different tissues of D. fragrans, and the first to discuss these findings in the context of the development of homologous glandular trichomes. These results set the stage for further research on the development, stress resistance, and secondary metabolism of D. fragrans glandular trichomes.

Does Gene Expression Noise Play a Functional Role in Plants?

Gene expression in individual cells can be surprisingly noisy. In unicellular organisms this noise can be functional; for example, by allowing a subfraction of the population to prepare for environmental stress. The role of gene expression noise in multicellular organisms has, however, remained unclear. In this review, we discuss how new techniques are revealing an unexpected level of variability in gene expression between and within genetically identical plants. We describe recent progress as well as speculate on the function of transcriptional noise as a mechanism for generating functional phenotypic diversity in plants.

Leaf Trichome Distribution Pattern in Arabidopsis Reveals Gene Expression Variation Associated with Environmental Adaptation

Gene expression varies stochastically even in both heterogenous and homogeneous cell populations. This variation is not simply useless noise; rather, it is important for many biological processes. Unicellular organisms or cultured cell lines are useful for analyzing the variation in gene expression between cells; however, owing to technical challenges, the biological relevance of this variation in multicellular organisms such as higher plants remain unclear. Here, we addressed the biological relevance of this variation between cells by examining the genetic basis of trichome distribution patterns in Arabidopsis thaliana. The distribution pattern of a trichome on a leaf is stochastic and can be mathematically represented using Turing’s reaction-diffusion (RD) model. We analyzed simulations based on the RD model and found that the variability in the trichome distribution pattern increased with the increase in stochastic variation in a particular gene expression. Moreover, differences in heat-dependent variability of the trichome distribution pattern between the accessions showed a strong correlation with environmental factors to which each accession was adapted. Taken together, we successfully visualized variations in gene expression by quantifying the variability in the Arabidopsis trichome distribution pattern. Thus, our data provide evidence for the biological importance of variations in gene expression for environmental adaptation.

Imaging flowers: a guide to current microscopy and tomography techniques to study flower development

Developmental biology relies heavily on our ability to generate three-dimensional images of live biological specimens through time, and to map gene expression and hormone response in these specimens as they undergo development. The last two decades have seen an explosion of new bioimaging technologies that have pushed the limits of spatial and temporal resolution and provided biologists with invaluable new tools. However, plant tissues are difficult to image, and no single technology fits all purposes; choosing between many bioimaging techniques is not trivial. Here, we review modern light microscopy and computed projection tomography methods, their capabilities and limitations, and we discuss their current and potential applications to the study of flower development and fertilization.

Can the French flag and reaction-diffusion models explain flower patterning? Celebrating the 50th anniversary of the French flag model

It has been 50 years since Lewis Wolpert introduced the French flag model proposing the patterning of different cell types based on threshold concentrations of a morphogen diffusing in the tissue. Sixty-seven years ago, Alan Turing introduced the idea of patterns initiating de novo from a reaction-diffusion network. Together these models have been used to explain many patterning events in animal development, so here we take a look at their applicability to flower development. First, although many plant transcription factors move through plasmodesmata from cell to cell, in the flower there is little evidence that they specify fate in a concentration-dependent manner, so they cannot yet be described as morphogens. Secondly, the reaction-diffusion model appears to be a reasonably good description of the formation of spots of pigment on petals, although additional nuances are present. Thirdly, aspects of both of these combine in a new fluctuation-based patterning system creating the scattered pattern of giant cells in Arabidopsis sepals. In the future, more precise imaging and manipulations of the dynamics of patterning networks combined with mathematical modeling will allow us to better understand how the multilayered complex and beautiful patterns of flowers emerge de novo.

Floral organ development goes live

The chance to watch floral organs develop live is not to be missed! Here, we outline reasons why quantitative, live-cell imaging is an important approach to study floral morphogenesis, and provide a basic workflow of how to get started. We highlight key advances in morphodynamics of lateral organ development, and discuss recent work that uses live confocal imaging to address the regulation of floral organ number, its robustness, and patterning mechanisms that exploit stochasticity.

MpFEW RHIZOIDS1 miRNA-Mediated Lateral Inhibition Controls Rhizoid Cell Patterning in Marchantia polymorpha

Lateral inhibition patterns differentiated cell types among equivalent cells during development in bacteria, metazoans, and plants. Tip-growing rhizoid cells develop among flat epidermal cells in the epidermis of the early-diverging land plant Marchantia polymorpha. We show that the majority of rhizoid cells develop individually, but some develop in linear, one-dimensional groups (chains) of between 2 and 7 rhizoid cells in wild-type plants. The distribution of rhizoid cells can be accounted for within a simple cellular automata model of lateral inhibition. The model predicted that in the absence of lateral inhibition, two-dimensional rhizoid cell groups (clusters) form. These can be larger than those formed with lateral inhibition. M. polymorpha rhizoid differentiation is positively regulated by the ROOT HAIR DEFECTIVE SIX-LIKE1 (MpRSL1) basic-helix-loop-helix (bHLH) transcription factor, which is directly repressed by the FEW RHIZOIDS1 (MpFRH1) microRNA (miRNA). To test if MpFRH1 miRNA acts during lateral inhibition, we generated loss-of-function (lof) mutants without the MpFRH1 miRNA. Two-dimensional clusters of rhizoids develop in Mpfrh1^lof mutants as predicted by the model for plants that lack lateral inhibition. Furthermore, two-dimensional clusters of up to 9 rhizoid cells developed in the Mpfrh1^lof mutants compared to a maximum number of 7 observed in wild-type groups. The higher steady-state levels of MpRSL1 mRNA in Mpfrh1^lof mutants indicate that MpFRH1-mediated lateral inhibition involves the repression of MpRSL1 activity. Together, the modeling and genetic data indicate that MpFRH1 miRNA mediates lateral inhibition by repressing MpRSL1 during pattern formation in the M. polymorpha epidermis.

A signal cascade originated from epidermis defines apical-basal patterning of Arabidopsis shoot apical meristems

In multicellular organisms, a long-standing question is how spatial patterns of distinct cell types are initiated and maintained during continuous cell division and proliferation. Along the vertical axis of plant shoot apical meristems (SAMs), stem cells are located at the top while cells specifying the stem cells are located more basally, forming a robust apical-basal pattern. We previously found that in Arabidopsis SAMs, the HAIRY MERISTEM (HAM) family transcription factors form a concentration gradient from the epidermis to the interior cell layers, and this gradient is essential for the stem cell specification and the apical-basal patterning of the SAMs. Here, we uncover that epidermis specific transcription factors, ARABIDOPSIS THALIANA MERISTEM LAYER 1 (ATML1) and its close homolog, define the concentration gradient of HAM in the SAM through activating a group of microRNAs. This study provides a molecular framework linking the epidermis-derived signal to the stem cell homeostasis in plants.

A concentration gradient of HAM transcription factors specifies apical-basal patterning in the Arabidopsis shoot apical meristem. Here, the authors show that epidermal expression of the ATML1 transcription factor defines this concentration gradient via activation of mobile micro RNA.

Stochastic pulsing of gene expression enables the generation of spatial patterns in Bacillus subtilis biofilms

Stochastic pulsing of gene expression can generate phenotypic diversity in a genetically identical population of cells, but it is unclear whether it has a role in the development of multicellular systems. Here, we show how stochastic pulsing of gene expression enables spatial patterns to form in a model multicellular system, Bacillus subtilis bacterial biofilms. We use quantitative microscopy and time-lapse imaging to observe pulses in the activity of the general stress response sigma factor σ^B in individual cells during biofilm development. Both σ^B and sporulation activity increase in a gradient, peaking at the top of the biofilm, even though σ^B represses sporulation. As predicted by a simple mathematical model, increasing σ^B expression shifts the peak of sporulation to the middle of the biofilm. Our results demonstrate how stochastic pulsing of gene expression can play a key role in pattern formation during biofilm development.

Stochastic pulsing of gene expression can generate phenotypic diversity in a genetically identical population of cells. Here, the authors show that stochastic pulsing in the expression of a sigma factor enables the formation of spatial patterns in a multicellular system, Bacillus subtilis bacterial biofilms.

Growth dynamics of the Arabidopsis fruit is mediated by cell expansion

Fruit have evolved a sophisticated tissue and cellular architecture to secure plant reproductive success. Postfertilization growth is perhaps the most dramatic event during fruit morphogenesis. Several studies have proposed that fertilized ovules and developing seeds initiate signaling cascades to coordinate and promote the growth of the accompanying fruit tissues. This dynamic process allows the fruit to conspicuously increase its size and acquire its final shape and means for seed dispersal. All these features are key for plant survival and crop yield. Despite its importance, we lack a high-resolution spatiotemporal map of how postfertilization fruit growth proceeds at the cellular level. In this study, we have combined live imaging, mutant backgrounds in which fertilization can be controlled, and computational modeling to monitor and predict postfertilization fruit growth in Arabidopsis We have uncovered that, unlike leaves, sepals, or roots, fruit do not exhibit a spatial separation of cell division and expansion domains; instead, there is a separation into temporal stages with fertilization as the trigger for transitioning to cell expansion, which drives postfertilization fruit growth. We quantified the coordination between fertilization and fruit growth by imaging no transmitting tract (ntt) mutants, in which fertilization fails in the bottom half of the fruit. By combining our experimental data with computational modeling, we delineated the mobility properties of the seed-derived signaling cascades promoting growth in the fruit. Our study provides the basis for generating a comprehensive understanding of the molecular and cellular mechanisms governing fruit growth and shape.

Cell Size Control in Plants

The genetic control of the characteristic cell sizes of different species and tissues is a long-standing enigma. Plants are convenient for studying this question in a multicellular context, as their cells do not move and are easily tracked and measured from organ initiation in the meristems to subsequent morphogenesis and differentiation. In this article, we discuss cell size control in plants compared with other organisms. As seen from yeast cells to mammalian cells, size homeostasis is maintained cell autonomously in the shoot meristem. In developing organs, vacuolization contributes to cell size heterogeneity and may resolve conflicts between growth control at the cellular and organ levels. Molecular mechanisms for cell size control have implications for how cell size responds to changes in ploidy, which are particularly important in plant development and evolution. We also discuss comparatively the functional consequences of cell size and their potential repercussions at higher scales, including genome evolution.

Double or Nothing? Cell Division and Cell Size Control

Size is a fundamental property that must be tightly regulated to ensure that cells and tissues function efficiently. Dynamic size control allows unicellular organisms to adapt to environmental changes, but cell size is also integral to multicellular development, affecting tissue size and structure. Despite clear evidence for homeostatic cell size maintenance, we are only now beginning to understand cell size regulation in the actively dividing meristematic tissues of higher plants. We discuss here how coupled advances in live cell imaging and modelling are uncovering dynamic mechanisms for size control mediated at the cellular level. We argue that integrated models of cell growth and division will be necessary to predict cell size and fully understand multicellular growth and development.

Paf1c defects challenge the robustness of flower meristem termination in Arabidopsis thaliana

Although accumulating evidence suggests that gene regulation is highly stochastic, genetic screens have successfully uncovered master developmental regulators, questioning the relationship between transcriptional noise and intrinsic robustness of development. To identify developmental modules that are more or less resilient to large-scale genetic perturbations, we used the Arabidopsis polymerase II-associated factor 1 complex (Paf1c) mutant vip3, which is impaired in several RNA polymerase II-dependent transcriptional processes. We found that the control of flower termination was not as robust as classically pictured. In angiosperms, the floral female organs, called carpels, display determinate growth: their development requires the arrest of stem cell maintenance. In vip3 mutant flowers, carpels displayed a highly variable morphology, with different degrees of indeterminacy defects up to wild-type size inflorescence emerging from carpels. This phenotype was associated with variable expression of two key regulators of flower termination and stem cell maintenance in flowers, WUSCHEL and AGAMOUS. The phenotype was also dependent on growth conditions. Together, these results highlight the surprisingly plastic nature of stem cell maintenance in plants and its dependence on Paf1c.

Summary: Using a mutant with increased transcriptional noise, we reveal that stem cell maintenance is not as robust as anticipated in plants, even leading to major defects in essential developmental processes such as flower indeterminacy.

Connected through the force: mechanical signals in plant development

As multicellular organisms, plants acquire characteristic shapes through a complex set of biological processes known as morphogenesis. Biochemical signalling underlies much of development, as it allows cells to acquire specific identities based on their position within tissues and organs. However, as growing physical structures, plants, and their constituent cells, also experience internal and external physical forces that can be perceived and can influence key processes such as growth, polarity, and gene expression. This process, which adds another layer of control to growth and development, has important implications for plant morphogenesis. This review provides an overview of recent research into the role of mechanical signals in plant development and aims to show how mechanical signalling can be used, in concert with biochemical signals, as a cue allowing cells and tissues to coordinate their behaviour and to add robustness to developmental processes.

Volumetric finite-element modelling of biological growth

Differential growth is the driver of tissue morphogenesis in plants, and also plays a fundamental role in animal development. Although the contributions of growth to shape change have been captured through modelling tissue sheets or isotropic volumes, a framework for modelling both isotropic and anisotropic volumetric growth in three dimensions over large changes in size and shape has been lacking. Here, we describe an approach based on finite-element modelling of continuous volumetric structures, and apply it to a range of forms and growth patterns, providing mathematical validation for examples that admit analytic solution. We show that a major difference between sheet and bulk tissues is that the growth of bulk tissue is more constrained, reducing the possibility of tissue conflict resolution through deformations such as buckling. Tissue sheets or cylinders may be generated from bulk shapes through anisotropic specified growth, oriented by a polarity field. A second polarity field, orthogonal to the first, allows sheets with varying lengths and widths to be generated, as illustrated by the wide range of leaf shapes observed in nature. The framework we describe thus provides a key tool for developing hypotheses for plant morphogenesis and is also applicable to other tissues that deform through differential growth or contraction.

Brassicaceae flowers: diversity amid uniformity

The mustard family Brassicaceae, which includes the model plant Arabidopsis thaliana, exhibits morphological stasis and significant uniformity of floral plan. Nonetheless, there is untapped diversity in almost every aspect of floral morphology in the family that lends itself to comparative study, including organ number, shape, form, and color. Studies on the genetic basis of morphological diversity, enabled by extensive genetic tools and genomic resources and the close phylogenetic distance among mustards, have revealed a mosaic of conservation and divergence in numerous floral traits. Here I review the morphological diversity of the flowers of Brassicaceae and discuss studies addressing the underlying genetic and developmental mechanisms shaping floral diversity. To put flowers in the context of the floral display, I describe diversity in inflorescence morphology and the variation that exists in the structures preceding the floral organs. Reconstructing the floral morphospace in Brassicaceae coupled with next-generation sequencing data and unbiased approaches to interrogate gene function in species throughout the mustard phylogeny offers promising ways to understand how developmental mechanisms originate and diversify.