Using positional information to provide context for biological image analysis with MorphoGraphX 2.0

Stochastic gene expression in auxin signaling in the floral meristem of Arabidopsis thaliana

Cells display striking stochasticity in gene expression, which plays an important role in development, disease, and regeneration. Previous studies have found stochastic gene expression in bacteria, yeast, and from constitutive promoters in Arabidopsis. However, most promoters are non-constitutive. Stochastic gene expression from non-constitutive promoters in a multicellular organism, especially those with key developmental roles, remains largely uncharacterized. Here, we report stochastic expression of auxin responsive genes in the Arabidopsis floral meristem, using promoter reporters of DR5, ARABIDOPSIS HISTIDINE PHOSPHOTRANSFER PROTEIN6 (AHP6), and DNA BINDING WITH ONE FINGER5.8 (DOF5.8). We find highly variable DR5 expression patterns among younger meristems. Such variability is largely attributed to stochastic expression of DR5, mainly influenced by cell-intrinsic molecular noise. Expression of AHP6 and DOF5.8 is also noisy, although their noise is lower and has distinct spatiotemporal patterns unlike DR5. Finally, we propose spatial averaging as a mechanism that buffers cellular gene expression noise, allowing the formation of robust global expression patterns. Our study reveals stochastic gene expression downstream of auxin signaling, a key developmental player. Thus, stochastic gene expression from non-constitutive promoters, including those involved in hormone signaling, is an ordinary part of multicellular life.

Single-cell transcriptomics reveal how root tissues adapt to soil stress

Land plants thrive in soils showing vastly different properties and environmental stresses1. Root systems can adapt to contrasting soil conditions and stresses, yet how their responses are programmed at the individual cell scale remains unclear. Using single-cell RNA sequencing and spatial transcriptomic approaches, we showed major expression changes in outer root cell types when comparing the single-cell transcriptomes of rice roots grown in gel versus soil conditions. These tissue-specific transcriptional responses are related to nutrient homeostasis, cell wall integrity and defence in response to heterogeneous soil versus homogeneous gel growth conditions. We also demonstrate how the model soil stress, termed compaction, triggers expression changes in cell wall remodelling and barrier formation in outer and inner root tissues, regulated by abscisic acid released from phloem cells. Our study reveals how root tissues communicate and adapt to contrasting soil conditions at single-cell resolution.

Camelot: a computer-automated micro-extensometer with low-cost optical tracking

BACKGROUND

Plant growth and morphogenesis is a mechanical process controlled by genetic and molecular networks. Measuring mechanical properties at various scales is necessary to understand how these processes interact. However, obtaining a device to perform the measurements on plant samples of choice poses technical challenges and is often limited by high cost and availability of specialized components, the adequacy of which needs to be verified. Developing software to control and integrate the different pieces of equipment can be a complex task.

RESULTS

To overcome these challenges, we have developed a computer automated micro-extensometer combined with low-cost optical tracking (Camelot) that facilitates measurements of elasticity, creep, and yield stress. It consists of three primary components: a force sensor with a sample attachment point, an actuator with a second attachment point, and a camera. To monitor force, we use a parallel beam sensor, commonly used in digital weighing scales. To stretch the sample, we use a stepper motor with a screw mechanism moving a stage along linear rail. To monitor sample deformation, a compact digital microscope or a microscope camera is used. The system is controlled by MorphoRobotX, an integrated open-source software environment for mechanical experimentation. We first tested the basic Camelot setup, equipped with a digital microscope to track landmarks on the sample surface. We demonstrate that the system has sufficient accuracy to measure the stiffness in delicate plant samples, the etiolated hypocotyls of Arabidopsis, and were able to measure stiffness differences between wild type and a xyloglucan-deficient mutant. Next, we placed Camelot on an inverted microscope and used a C-mount microscope camera to track displacement of cell junctions. We stretched onion epidermal peels in longitudinal and transverse directions and obtained results similar to those previously published. Finally, we used the setup coupled with an upright confocal microscope and measured anisotropic deformation of individual epidermal cells during stretching of an Arabidopsis leaf.

CONCLUSIONS

The portability and suitability of Camelot for high-resolution optical tracking under a microscope make it an ideal tool for researchers in resource-limited settings or those pursuing exploratory biomechanics work.

CLAVATA signalling shapes barley inflorescence by controlling activity and determinacy of shoot meristem and rachilla

The large variety of inflorescence architectures evolved in grasses depends on shape, longevity and determinacy of meristems directing growth of the main and lateral axes. The CLAVATA pathway is known to regulate meristem size and inflorescence architecture in grasses. However, how individual meristem activities are determined and integrated to generate specific inflorescences is not yet understood. We found that activity of distinct meristems in the barley inflorescence is controlled by a signalling pathway comprising the receptor-like kinase Hordeum vulgare CLAVATA1 (HvCLV1) and the secreted CLAVATA3/EMBRYO-SURROUNDING REGION RELATED (CLE)-family peptide FON2-LIKE CLE PROTEIN1 (HvFCP1). HvFCP1 and HvCLV1 interact to promote spikelet formation, but restrict inflorescence meristem and rachilla proliferation. Hvfcp1 or Hvclv1 mutants generate additional rows of spikelets and supernumerary florets from extended rachilla activity. HvFCP1/HvCLV1 signalling coordinates meristem activity through regulation of trehalose-6-phosphate levels. Our discoveries outline a path to engineer inflorescence architecture via specific regulation of distinct meristem activities.

Polarity-guided uneven mitotic divisions control brassinosteroid activity in proliferating plant root cells

Brassinosteroid hormones are positive regulators of plant organ growth, yet their function in proliferating tissues remains unclear. Here, through integrating single-cell RNA sequencing with long-term live-cell imaging of the Arabidopsis root, we reveal that brassinosteroid activity fluctuates throughout the cell cycle, decreasing during mitotic divisions and increasing during the G1 phase. The post-mitotic recovery of brassinosteroid activity is driven by the intrinsic polarity of the mother cell, resulting in one daughter cell with enhanced brassinosteroid signaling, while the other supports brassinosteroid biosynthesis. The coexistence of these distinct daughter cell states during the G1 phase circumvents a negative feedback loop to facilitate brassinosteroid production while signaling increases. Our findings uncover polarity-guided, uneven mitotic divisions in the meristem, which control brassinosteroid hormone activity to ensure optimal root growth.

Growth directions and stiffness across cell layers determine whether tissues stay smooth or buckle

From smooth shapes to buckles, nature exhibits organs of various shapes and forms. How cells grow to produce smooth shaped leaves and sepals remain unclear. Here, we show that growth along the longitudinal axis during early developmental stages and comparable stiffness across both epidermal layers of Arabidopsis sepals are essential for smoothness, as seen in the wild type. We identified a mutant (as2-7D) with ectopic expression of ASYMMETRIC LEAVES 2 (AS2) on the outer epidermis. Our analysis reveals that ectopic AS2 expression causes the outer epidermis of as2-7D sepals to buckle during early stages of sepal development. We show that buckling of the outer epidermis occurs due to conflicting cell growth directions and unequal tissue stiffness across the epidermal layers. Overexpression of cyclin-dependent kinase (CDK) inhibitor Kip-related protein 1 (KRP1) in as2-7D restores sepal smoothness by aligning the growth directions of the outer epidermal cells along the longitudinal axis, while also increasing the overall stiffness of the outer epidermis. Furthermore, buckling is associated with the convergence of auxin efflux transporter protein PIN-FORMED 1 (PIN1) to generate outgrowth in the sepals at later stages, suggesting that buckling can initiate outgrowths. Our findings suggest that in addition to molecular cues influencing tissue mechanics, tissue mechanics can also modulate molecular signals, giving rise to well-defined shapes.

Mechanical interactions between tissue layers underlie plant morphogenesis

Differential growth between tissues generates mechanical conflicts influencing organogenesis in plants. Here we use the anther, the male floral reproductive organ, as a model system to understand how cell dynamics and tissue-scale mechanics control 3D morphogenesis of a complex shape. Combining deep live-cell imaging, growth analysis, osmotic treatments, genetics and mechanical modelling, we show that localized expansion of internal cells actively drives anther lobe outgrowth, while the epidermis stretches in response. At later stages, mechanical load is transferred to the sub-epidermal layer (endothecium), contributing to proper organ shape. We propose the concept of 'inflation potential', encapsulating mechanical and anatomical features causing differential growth. Our data emphasize the active mechanical role of inner tissue in controlling both organ shape acquisition and cell dynamics in outer layers.

Symplastic guard cell connections buffer pressure fluctuations to promote stomatal function in grasses

Stomata regulate plant gas exchange via repeated turgor-driven changes of guard cell shape, thereby adjusting pore apertures. Grasses, which are among the most widespread plant families on the planet, are distinguished by their unique stomatal structure, which is proposed to have significantly contributed to their evolutionary and agricultural success. One component of their structure, which has received little attention, is the presence of a discontinuous adjoining cell wall of the guard cell pair. Here, we demonstrate the presence of these symplastic connections in a range of grasses and use finite element method simulations to assess hypotheses for their functional significance. Our results show that opening of the stomatal pore is maximal when the turgor pressure in dumbbell-shaped grass guard cells is equal, especially under the low pressure conditions that occur during the early phase of stomatal opening. By contrast, we demonstrate that turgor pressure differences have less effect on the opening of kidney-shaped guard cells, characteristic of the majority of land plants, where guard cell connections are rarely or not observed. Our data describe a functional mechanism based on cellular mechanics, which plausibly facilitated a major transition in plant evolution and crop development.

A Method to Visualize Cell Proliferation of Arabidopsis thaliana: A Case Study of the Root Apical Meristem

Plant growth and development rely on a delicate balance between cell proliferation and cell differentiation. The root apical meristem (RAM) of Arabidopsis thaliana is an excellent model to study the cell cycle due to the coordinated relationship between nucleus shape and cell size at each stage, allowing for precise estimation of the cell cycle duration. In this study, we present a method for high-resolution visualization of RAM cells. This is the first protocol that allows for simultaneous high-resolution imaging of cellular and nuclear stains, being compatible with DNA replication markers such as EdU, including fluorescent proteins (H2B::YFP), SYTOX DNA stains, and the cell wall stain SR2200. This protocol includes a clarification procedure that enables the acquisition of high-resolution 3D images, suitable for detailed subsequent analysis.

Spatial ploidy inference using quantitative imaging

Polyploidy (whole-genome multiplication) is a common yet under-surveyed property of tissues across multicellular organisms. Polyploidy plays a critical role during tissue development, following acute stress, and during disease progression. Common methods to reveal polyploidy involve either destroying tissue architecture by cell isolation or by tedious identification of individual nuclei in intact tissue. Therefore, there is a critical need for rapid and high-throughput ploidy quantification using images of nuclei in intact tissues. Here, we present iSPy (Inferring Spatial Ploidy), a new unsupervised learning pipeline that is designed to create a spatial map of nuclear ploidy across a tissue of interest. We demonstrate the use of iSPy in Arabidopsis, Drosophila, and human tissue. iSPy can be adapted for a variety of tissue preparations, including whole mount and sectioned. This high-throughput pipeline will facilitate rapid and sensitive identification of nuclear ploidy in diverse biological contexts and organisms.

WOX1 controls leaf serration development via temporally restricting BRASSINAZOLE RESISTANT 1 and CUP SHAPED COTYLEDON 3 expression in Arabidopsis

Leaves have evolved shape diversity, ranging from simple leaves with a smooth margin to complicated shapes with toothed/serrated, lobed, and dissected leaves with leaflets. In the model plant Arabidopsis thaliana with simple leaves producing a serrated margin, boundary regulatory factor genes CUP SHAPED COTYLEDON 2 (CUC2) and CUC3 play important roles in promoting leaf initiation and maintenance of serration. Stem cell-related WUSCHEL-RELATED HOMEOBOX1 (WOX1) and PRESSED FLOWER (PRS)/WOX3 are also essential for leaf margin morphogenesis, but the role of WOX1 and PRS as well as the relationships between CUCs and WOXs for tooth development are unclear. In this study, we found that WOX1, but not PRS, prevents overproduction of number of teeth and excessive tooth size by limiting CUC3 expression to a moderate level in a temporally regulated manner. We also revealed that BRASSINAZOLE RESISTANT 1 (BZR1), a known regulator of plant development including boundary regions, is involved in WOX1 negative regulation of tooth development by repressing CUC3 expression during the initiation/early stage of tooth development. WOX1 parallelly limits BZR1 and CUC3 expression from the late stage of the first two teeth, while it restricts CUC3 activity in a BZR1-dependent manner from the initiation/early stage of subsequently developed teeth. This study uncovers a new mechanism for WOX1 in fine-tuning the leaf margin geometry.

Evolution of a SHOOTMERISTEMLESS transcription factor binding site promotes fruit shape determination

In animals and plants, organ shape is primarily determined during primordium development by carefully coordinated growth and cell division1-3. Rare examples of post-primordial change in morphology (reshaping) exist that offer tractable systems for the study of mechanisms required for organ shape determination and diversification. One such example is morphogenesis in Capsella fruits whose heart-shaped appearance emerges by reshaping of the ovate spheroid gynoecium upon fertilization4. Here we use whole-organ live-cell imaging and single-cell RNA sequencing (scRNA-seq) analysis to show that Capsella fruit shape determination is based on dynamic changes in cell growth and cell division coupled with local maintenance of meristematic identity. At the molecular level, we reveal an auxin-induced mechanism that is required for morphological alteration and ultimately determined by a single cis-regulatory element. This element resides in the promoter of the Capsella rubella SHOOTMERISTEMLESS5 (CrSTM) gene. The CrSTM meristem identity factor positively regulates its own expression through binding to this element, thereby providing a feed-forward loop at the position and time of protrusion emergence to form the heart. Independent evolution of the STM-binding element in STM promoters across Brassicaceae species correlates with those undergoing a gynoecium-to-fruit shape change. Accordingly, genetic and phenotypic studies show that the STM-binding element is required to facilitate the shape transition and suggest a conserved molecular mechanism for organ morphogenesis.

Actin-mediated avoidance of tricellular junction influences global topology at the Arabidopsis shoot apical meristem

Division plane orientation contributes to cell shape and topological organization, playing a key role in morphogenesis, but the precise physical and molecular mechanism influencing these processes remains largely obscure in plants. In particular, it is less clear how the placement of the new walls occurs in relation to the walls of neighboring cells. Here, we show that genetic perturbation of the actin cytoskeleton results in more rectangular cell shapes and higher incidences of four-way junctions, perturbing the global topology of cells in the shoot apical meristem of Arabidopsis thaliana. Actin mutants also exhibit changes in the expansion rate of the new versus the maternal cell wall after division, affecting the evolution of internal angles at tricellular junctions. Further, the increased width of the preprophase band in the actin mutant contributes to inaccuracy in the placement of the new cell wall. Computational simulation further substantiates this hypothesis and reproduces the observed cell shape defects.

Topological analysis of 3D digital ovules identifies cellular patterns associated with ovule shape diversity

Tissue morphogenesis remains poorly understood. In plants, a central problem is how the 3D cellular architecture of a developing organ contributes to its final shape. We address this question through a comparative analysis of ovule morphogenesis, taking advantage of the diversity in ovule shape across angiosperms. Here, we provide a 3D digital atlas of Cardamine hirsuta ovule development at single cell resolution and compare it with an equivalent atlas of Arabidopsis thaliana. We introduce nerve-based topological analysis as a tool for unbiased detection of differences in cellular architectures and corroborate identified topological differences between two homologous tissues by comparative morphometrics and visual inspection. We find that differences in topology, cell volume variation and tissue growth patterns in the sheet-like integuments and the bulbous chalaza are associated with differences in ovule curvature. In contrast, the radialized conical ovule primordia and nucelli exhibit similar shapes, despite differences in internal cellular topology and tissue growth patterns. Our results support the notion that the structural organization of a tissue is associated with its susceptibility to shape changes during evolutionary shifts in 3D cellular architecture.

Age-associated growth control modifies leaf proximodistal symmetry and enabled leaf shape diversification

Biological shape diversity is often manifested in modulation of organ symmetry and modification of the patterned elaboration of repeated shape elements.1,2,3,4,5 Whether and how these two aspects of shape determination are coordinately regulated is unclear.5,6,7 Plant leaves provide an attractive system to investigate this problem, because they often show asymmetries along the proximodistal (PD) axis of their blades, along which they can also produce repeated marginal outgrowths such as serrations or leaflets.1 One aspect of leaf shape diversity is heteroblasty, where the leaf form in a single genotype is modified with progressive plant age.8,9,10,11 In Arabidopsis thaliana, a plant with simple leaves, SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 9 (SPL9) controls heteroblasty by activating CyclinD3 expression, thereby sustaining proliferative growth and retarding differentiation in adult leaves.12,13 However, the precise significance of SPL9 action for leaf symmetry and marginal patterning is unknown. By combining genetics, quantitative shape analyses, and time-lapse imaging, we show that PD symmetry of the leaf blade in A. thaliana decreases in response to an age-dependent SPL9 expression gradient, and that SPL9 action coordinately regulates the distribution and shape of marginal serrations and overall leaf form. Using comparative analyses, we demonstrate that heteroblastic growth reprogramming in Cardamine hirsuta, a complex-leafed relative of A. thaliana, also involves prolonging the duration of cell proliferation and delaying differentiation. We further provide evidence that SPL9 enables species-specific action of homeobox genes that promote leaf complexity. In conclusion, we identified an age-dependent layer of organ PD growth regulation that modulates leaf symmetry and has enabled leaf shape diversification.

Robust organ size in Arabidopsis is primarily governed by cell growth rather than cell division patterns

Organ sizes and shapes are highly reproducible, or robust, within a species and individuals. Arabidopsis thaliana sepals, which are the leaf-like organs that enclose flower buds, have consistent size and shape, indicating robust development. Cell growth is locally heterogeneous due to intrinsic and extrinsic noise. To achieve robust organ shape, fluctuations in cell growth must average to an even growth rate, which requires that fluctuations are uncorrelated or anti-correlated in time and space. Here, we live image and quantify the development of sepals with an increased or decreased number of cell divisions (lgo mutant and LGO overexpression, respectively), a mutant with altered cell growth variability (ftsh4), and double mutants combining these. Changes in the number of cell divisions do not change the overall growth pattern. By contrast, in ftsh4 mutants, cell growth accumulates in patches of over- and undergrowth owing to correlations that impair averaging, resulting in increased organ shape variability. Thus, we demonstrate in vivo that the number of cell divisions does not affect averaging of cell growth, preserving robust organ morphogenesis, whereas correlated growth fluctuations impair averaging.

Modulation of cell differentiation and growth underlies the shift from bud protection to light capture in cauline leaves

Plant organs have evolved into diverse shapes for specialized functions despite emerging as simple protrusions at the shoot apex. Cauline leaves serve as photosynthetic organs and protective structures for emerging floral buds. However, the growth patterns underlying this dual function remain unknown. Here, we investigate the developmental dynamics shaping Arabidopsis (Arabidopsis thaliana) cauline leaves underlying their functional diversification from other laminar organs. We show that cauline leaves display a significant delay in overall elongation compared with rosette leaves. Using live imaging, we reveal that their functional divergence hinges on early modulation of the timing of cell differentiation and cellular growth rates. In contrast to rosette leaves and sepals, cell differentiation is delayed in cauline leaves, fostering extended proliferation, prolonged morphogenetic activity, and growth redistribution within the organ. Notably, cauline leaf growth is transiently suppressed during the early stages, keeping the leaf small and unfolded during the initiation of the first flowers. Our findings highlight the unique developmental timing of cauline leaves, underlying their shift from an early protective role to a later photosynthetic function.

An optimized live imaging and multiple cell layer growth analysis approach using Arabidopsis sepals

Arabidopsis thaliana sepals are excellent models for analyzing growth of entire organs due to their relatively small size, which can be captured at a cellular resolution under a confocal microscope. To investigate how differential growth of connected cell layers generate unique organ morphologies, it is necessary to live-image deep into the tissue. However, imaging deep cell layers of the sepal (or plant tissues in general) is practically challenging. Image processing is also difficult due to the low signal-to-noise ratio of the deeper tissue layers, an issue mainly associated with live imaging datasets. Addressing some of these challenges, we provide an optimized methodology for live imaging sepals, and subsequent image processing. For live imaging early-stage sepals, we found that the use of a bright fluorescent membrane marker, coupled with increased laser intensity and an enhanced Z- resolution produces high-quality images suitable for downstream image processing. Our optimized parameters allowed us to image the bottommost cell layer of the sepal (inner epidermal layer) without compromising viability. We used a 'voxel removal' technique to visualize the inner epidermal layer in MorphoGraphX image processing software. We also describe the MorphoGraphX parameters for creating a 2.5D mesh surface for the inner epidermis. Our parameters allow for the segmentation and parent tracking of individual cells through multiple time points, despite the weak signal of the inner epidermal cells. While we have used sepals to illustrate our approach, the methodology will be useful for researchers intending to live-image and track growth of deeper cell layers in 2.5D for any plant tissue.

An incoherent feed-forward loop involving bHLH transcription factors, Auxin and CYCLIN-Ds regulates style radial symmetry establishment in Arabidopsis

The bilateral-to-radial symmetry transition occurring during the development of the Arabidopsis thaliana female reproductive organ (gynoecium) is a crucial biological process linked to plant fertilization and seed production. Despite its significance, the cellular mechanisms governing the establishment and breaking of radial symmetry at the gynoecium apex (style) remain unknown. To fill this gap, we employed quantitative confocal imaging coupled with MorphoGraphX analysis, in vivo and in vitro transcriptional experiments, and genetic analysis encompassing mutants in two bHLH transcription factors necessary and sufficient to promote transition to radial symmetry, SPATULA (SPT) and INDEHISCENT (IND). Here, we show that defects in style morphogenesis correlate with defects in cell-division orientation and rate. We showed that the SPT-mediated accumulation of auxin in the medial-apical cells undergoing symmetry transition is required to maintain cell-division-oriented perpendicular to the direction of organ growth (anticlinal, transversal cell division). In addition, SPT and IND promote the expression of specific core cell-cycle regulators, CYCLIN-D1;1 (CYC-D1;1) and CYC-D3;3, to support progression through the G1 phase of the cell cycle. This transcriptional regulation is repressed by auxin, thus forming an incoherent feed-forward loop mechanism. We propose that this mechanism fine-tunes cell division rate and orientation with the morphogenic signal provided by auxin, during patterning of radial symmetry at the style.

A deep learning-based toolkit for 3D nuclei segmentation and quantitative analysis in cellular and tissue context

We present a new set of computational tools that enable accurate and widely applicable 3D segmentation of nuclei in various 3D digital organs. We have developed an approach for ground truth generation and iterative training of 3D nuclear segmentation models, which we applied to popular CellPose, PlantSeg and StarDist algorithms. We provide two high-quality models trained on plant nuclei that enable 3D segmentation of nuclei in datasets obtained from fixed or live samples, acquired from different plant and animal tissues, and stained with various nuclear stains or fluorescent protein-based nuclear reporters. We also share a diverse high-quality training dataset of about 10,000 nuclei. Furthermore, we advanced the MorphoGraphX analysis and visualization software by, among other things, providing a method for linking 3D segmented nuclei to their surrounding cells in 3D digital organs. We found that the nuclear-to-cell volume ratio varies between different ovule tissues and during the development of a tissue. Finally, we extended the PlantSeg 3D segmentation pipeline with a proofreading tool that uses 3D segmented nuclei as seeds to correct cell segmentation errors in difficult-to-segment tissues.

Tradeoff between speed and robustness in primordium initiation mediated by auxin-CUC1 interaction

Robustness is the reproducible development of a phenotype despite stochastic noise. It often involves tradeoffs with other performance metrics, but the mechanisms underlying such tradeoffs were largely unknown. An Arabidopsis flower robustly develops four sepals from four precisely positioned auxin maxima. The development related myb-like 1 (drmy1) mutant generates noise in auxin signaling that disrupts robustness in sepal initiation. Here, we find that increased expression of CUP-SHAPED COTYLEDON1 (CUC1), a boundary specification transcription factor, in drmy1 underlies this loss of robustness. CUC1 surrounds and amplifies stochastic auxin noise in drmy1 to form variably positioned auxin maxima and sepal primordia. Removing CUC1 from drmy1 provides time for noisy auxin signaling to resolve into four precisely positioned auxin maxima, restoring robust sepal initiation. However, removing CUC1 decreases the intensity of auxin maxima and slows down sepal initiation. Thus, CUC1 increases morphogenesis speed but impairs robustness against auxin noise. Further, using a computational model, we find that the observed phenotype can be explained by the effect of CUC1 in repolarizing PIN FORMED1 (PIN1), a polar auxin transporter. Lastly, our model predicts that reducing global growth rate improves developmental robustness, which we validate experimentally. Thus, our study illustrates a tradeoff between speed and robustness during development.

Spatiotemporally distinct responses to mechanical forces shape the developing seed of Arabidopsis

Organ morphogenesis depends on mechanical interactions between cells and tissues. These interactions generate forces that can be sensed by cells and affect key cellular processes. However, how mechanical forces, together with biochemical signals, contribute to the shaping of complex organs is still largely unclear. We address this question using the seed of Arabidopsis as a model system. We show that seeds first experience a phase of rapid anisotropic growth that is dependent on the response of cortical microtubule (CMT) to forces, which guide cellulose deposition according to shape-driven stresses in the outermost layer of the seed coat. However, at later stages of development, we show that seed growth is isotropic and depends on the properties of an inner layer of the seed coat that stiffens its walls in response to tension but has isotropic material properties. Finally, we show that the transition from anisotropic to isotropic growth is due to the dampening of cortical microtubule responses to shape-driven stresses. Altogether, our work supports a model in which spatiotemporally distinct mechanical responses control the shape of developing seeds in Arabidopsis.

Optical microscopic imaging, manipulation, and analysis methods for morphogenesis research

Morphogenesis is a developmental process of organisms being shaped through complex and cooperative cellular movements. To understand the interplay between genetic programs and the resulting multicellular morphogenesis, it is essential to characterize the morphologies and dynamics at the single-cell level and to understand how physical forces serve as both signaling components and driving forces of tissue deformations. In recent years, advances in microscopy techniques have led to improvements in imaging speed, resolution and depth. Concurrently, the development of various software packages has supported large-scale, analyses of challenging images at the single-cell resolution. While these tools have enhanced our ability to examine dynamics of cells and mechanical processes during morphogenesis, their effective integration requires specialized expertise. With this background, this review provides a practical overview of those techniques. First, we introduce microscopic techniques for multicellular imaging and image analysis software tools with a focus on cell segmentation and tracking. Second, we provide an overview of cutting-edge techniques for mechanical manipulation of cells and tissues. Finally, we introduce recent findings on morphogenetic mechanisms and mechanosensations that have been achieved by effectively combining microscopy, image analysis tools and mechanical manipulation techniques.

Development and application of a 3D image analysis strategy for focused ion beam - Scanning electron microscopy tomography of porous soft materials

In recent years, the potential of porous soft materials in various device technologies has increased in importance due to applications in fields, such as wearable electronics, medicine, and transient devices. However, understanding the 3-dimensional architecture of porous soft materials at the microscale remains a challenge. Herein, we present a method to structurally analyze soft materials using Focused Ion Beam - Scanning Electron Microscopy (FIB-SEM) tomography. Two materials, polymethyl methacrylate (PMMA) membrane and pine wood veneer were chosen as test-cases. FIB-SEM was successfully used to reconstruct the true topography of these materials in 3D. Structural and physical properties were subsequently deduced from the rendered 3D models. The methodology used segmentation, coupled with optimized thresholding, image processing, and reconstruction protocols. The 3D models generated pore size distribution, pore inter-connectivity, tortuosity, thickness, and curvature data. It was shown that FIB-SEM tomography provides both an informative and visual depiction of structure. To evaluate and validate the FIB-SEM reconstructions, porous properties were generated from the physical property analysis techniques, gas adsorption analysis using Brunauer-Emmett-Teller (BET) surface area analysis and mercury intrusion porosimetry (MIP) analysis. In general, the data obtained from the FIB-SEM reconstructions was well-matched with the physical data. RESEARCH HIGHLIGHTS: Porous specimens of both synthetic and biological nature, a poly(methyl methacrylate) membrane and a pine veneer respectively, are reconstructed via FIB-SEM tomography without resin-embedding. Different thresholding and reconstruction methods are explored whereby shadowing artifacts are present with the aid of free open-source software. Reconstruction data is compared to physical data: MIP, gas adsorption isotherms which are analyzed via BET and Barrett-Joyner-Halenda (BJH) analysis to yield a full picture of the materials.

Tradeoff Between Speed and Robustness in Primordium Initiation Mediated by Auxin-CUC1 Interaction

Robustness is the reproducible development of a phenotype despite stochastic noise. It often involves tradeoffs with other performance metrics, but the mechanisms underlying such tradeoffs were largely unknown. An Arabidopsis flower robustly develops four sepals from four precisely positioned auxin maxima. The development related myb-like 1 (drmy1) mutant generates noise in auxin signaling that disrupts robustness in sepal initiation. Here, we found that increased expression of CUP-SHAPED COTYLEDON1 (CUC1), a boundary specification transcription factor, in drmy1 underlies this loss of robustness. CUC1 surrounds and amplifies stochastic auxin noise in drmy1 to form variably positioned auxin maxima and sepal primordia. Removing CUC1 from drmy1 provides time for noisy auxin signaling to resolve into four precisely positioned auxin maxima, restoring robust sepal initiation. However, removing CUC1 decreases auxin maxima intensity and slows down sepal initiation. Thus, CUC1 increases morphogenesis speed but impairs robustness against auxin noise. Further, using a computational model, we found that the observed phenotype can be explained by the effect of CUC1 in repolarizing PIN FORMED1 (PIN1), a polar auxin transporter. Lastly, our model predicts that reducing global growth rate improves developmental robustness, which we validated experimentally. Thus, our study illustrates a tradeoff between speed and robustness during development.

Mechanical forces in plant tissue matrix orient cell divisions via microtubule stabilization

Plant morphogenesis relies exclusively on oriented cell expansion and division. Nonetheless, the mechanism(s) determining division plane orientation remain elusive. Here, we studied tissue healing after laser-assisted wounding in roots of Arabidopsis thaliana and uncovered how mechanical forces stabilize and reorient the microtubule cytoskeleton for the orientation of cell division. We identified that root tissue functions as an interconnected cell matrix, with a radial gradient of tissue extendibility causing predictable tissue deformation after wounding. This deformation causes instant redirection of expansion in the surrounding cells and reorientation of microtubule arrays, ultimately predicting cell division orientation. Microtubules are destabilized under low tension, whereas stretching of cells, either through wounding or external aspiration, immediately induces their polymerization. The higher microtubule abundance in the stretched cell parts leads to the reorientation of microtubule arrays and, ultimately, informs cell division planes. This provides a long-sought mechanism for flexible re-arrangement of cell divisions by mechanical forces for tissue reconstruction and plant architecture.

Two orthogonal differentiation gradients locally coordinate fruit morphogenesis

Morphogenesis requires the coordination of cellular behaviors along developmental axes. In plants, gradients of growth and differentiation are typically established along a single longitudinal primordium axis to control global organ shape. Yet, it remains unclear how these gradients are locally adjusted to regulate the formation of complex organs that consist of diverse tissue types. Here we combine quantitative live imaging at cellular resolution with genetics, and chemical treatments to understand the formation of Arabidopsis thaliana female reproductive organ (gynoecium). We show that, contrary to other aerial organs, gynoecium shape is determined by two orthogonal, time-shifted differentiation gradients. An early mediolateral gradient controls valve morphogenesis while a late, longitudinal gradient regulates style differentiation. Local, tissue-dependent action of these gradients serves to fine-tune the common developmental program governing organ morphogenesis to ensure the specialized function of the gynoecium.

Role of PMR4 and PDLP1 in priming of early acting penetration defense by resistance-inducing β-amino acids

R-β-homoserine (RBH) and β-aminobutyric acid (BABA) induce resistance against the oomycete Hyaloperonospora arabidopsidis (Hpa) in Arabidopsis, which is based on priming of multiple defense layers, including early acting penetration resistance at the cell wall. Here, we have examined the molecular basis of RBH- and BABA-primed defense by cell wall papillae against Hpa. Three-dimensional reconstruction of Hpa-induced papillae by confocal microscopy revealed no structural differences between control-, RBH-, and BABA-treated plants after Hpa challenge. However, mutations affecting POWDERY MILDEW RESISTANCE 4 or PLASMODESMATA LOCATED PROTEINs (PDLPs) only impaired BABA-induced penetration resistance and not RBH-induced penetration resistance. Furthermore, PDLP1 over-expression mimicked primed penetration resistance, while the intensity of GFP-tagged PDLP1 at germinating Hpa conidiospores was increased in BABA-primed plants but not RBH-primed plants. Our study reveals new regulatory layers of immune priming by β-amino acids and supports the notion that penetration resistance is a multifaceted defense layer that can be achieved through seperate pathways.

Growth and tension in explosive fruit

Exploding seed pods of the common weed Cardamine hirsuta have the remarkable ability to launch seeds far from the plant. The energy for this explosion comes from tension that builds up in the fruit valves. Above a critical threshold, the fruit fractures along its dehiscence zone and the two valves coil explosively, ejecting the seeds. A common mechanism to generate tension is drying, causing tissues to shrink. However, this does not happen in C. hirsuta fruit. Instead, tension is produced by active contraction of growing exocarp cells in the outer layer of the fruit valves. Exactly how growth causes the exocarp tissue to contract and generate pulling force is unknown. Here we show that the reorientation of microtubules in the exocarp cell cortex changes the orientation of cellulose microfibrils in the cell wall and the consequent cellular growth pattern. We used mechanical modeling to show how tension emerges through growth due to the highly anisotropic orientation of load-bearing cellulose microfibrils and their effect on cell shape. By explicitly defining the cell wall as multi-layered in our model, we discovered that a cross-lamellate pattern of cellulose microfibrils further enhances the developing tension in growing cells. Therefore, the interplay of cell wall properties with turgor-driven growth enables the fruit exocarp to generate sufficient tension to power explosive seed dispersal.

Cell-cycle-linked growth reprogramming encodes developmental time into leaf morphogenesis

How is time encoded into organ growth and morphogenesis? We address this question by investigating heteroblasty, where leaf development and form are modified with progressing plant age. By combining morphometric analyses, fate-mapping through live-imaging, computational analyses, and genetics, we identify age-dependent changes in cell-cycle-associated growth and histogenesis that underpin leaf heteroblasty. We show that in juvenile leaves, cell proliferation competence is rapidly released in a "proliferation burst" coupled with fast growth, whereas in adult leaves, proliferative growth is sustained for longer and at a slower rate. These effects are mediated by the SPL9 transcription factor in response to inputs from both shoot age and individual leaf maturation along the proximodistal axis. SPL9 acts by activating CyclinD3 family genes, which are sufficient to bypass the requirement for SPL9 in the control of leaf shape and in heteroblastic reprogramming of cellular growth. In conclusion, we have identified a mechanism that bridges across cell, tissue, and whole-organism scales by linking cell-cycle-associated growth control to age-dependent changes in organ geometry.

An optimized live imaging and growth analysis approach for Arabidopsis Sepals

Background

Arabidopsis thaliana sepals are excellent models for analyzing growth of entire organs due to their relatively small size, which can be captured at a cellular resolution under a confocal microscope [1]. To investigate how growth of different tissue layers generates unique organ morphologies, it is necessary to live-image deep into the tissue. However, imaging deep cell layers of the sepal is practically challenging, as it is hindered by the presence of extracellular air spaces between mesophyll cells, among other factors which causes optical aberrations. Image processing is also difficult due to the low signal-to-noise ratio of the deeper tissue layers, an issue mainly associated with live imaging datasets. Addressing some of these challenges, we provide an optimized methodology for live imaging sepals and subsequent image processing. This helps us track the growth of individual cells on the outer and inner epidermal layers, which are the key drivers of sepal morphogenesis.

Results

For live imaging sepals across all tissue layers at early stages of development, we found that the use of a bright fluorescent membrane marker, coupled with increased laser intensity and an enhanced Z- resolution produces high-quality images suitable for downstream image processing. Our optimized parameters allowed us to image the bottommost cell layer of the sepal (inner epidermal layer) without compromising viability. We used a 'voxel removal' technique to visualize the inner epidermal layer in MorphoGraphX [2, 3] image processing software. Finally, we describe the process of optimizing the parameters for creating a 2.5D mesh surface for the inner epidermis. This allowed segmentation and parent tracking of individual cells through multiple time points, despite the weak signal of the inner epidermal cells.

Conclusion

We provide a robust pipeline for imaging and analyzing growth across inner and outer epidermal layers during early sepal development. Our approach can potentially be employed for analyzing growth of other internal cell layers of the sepals as well. For each of the steps, approaches, and parameters we used, we have provided in-depth explanations to help researchers understand the rationale and replicate our pipeline.

Enabling neighbour labelling: using synthetic biology to explore how cells influence their neighbours

Cell-cell interactions are central to development, but exploring how a change in any given cell relates to changes in the neighbour of that cell can be technically challenging. Here, we review recent developments in synthetic biology and image analysis that are helping overcome this problem. We highlight the opportunities presented by these advances and discuss opportunities and limitations in applying them to developmental model systems.

Robust organ size in Arabidopsis is primarily governed by cell growth rather than cell division patterns

Organ sizes and shapes are highly reproducible, or robust, within a species and individuals. Arabidopsis thaliana sepals, which are the leaf-like organs that enclose flower buds, have consistent size and shape, which indicates robust development. Counterintuitively, variability in cell growth rate over time and between cells facilitates robust development because cumulative cell growth averages to a uniform rate. Here we investigate how sepal morphogenesis is robust to changes in cell division but not robust to changes in cell growth variability. We live image and quantitatively compare the development of sepals with increased or decreased cell division rate (lgo mutant and LGO overexpression, respectively), a mutant with altered cell growth variability (ftsh4), and double mutants combining these. We find that robustness is preserved when cell division rate changes because there is no change in the spatial pattern of growth. Meanwhile when robustness is lost in ftsh4 mutants, cell growth accumulates unevenly, and cells have disorganized growth directions. Thus, we demonstrate in vivo that both cell growth rate and direction average in robust development, preserving robustness despite changes in cell division.

Stiffness transitions in new walls post-cell division differ between Marchantia polymorpha gemmae and Arabidopsis thaliana leaves

Plant morphogenesis is governed by the mechanics of the cell wall-a stiff and thin polymeric box that encloses the cells. The cell wall is a highly dynamic composite material. New cell walls are added during cell division. As the cells continue to grow, the properties of cell walls are modulated to undergo significant changes in shape and size without breakage. Spatial and temporal variations in cell wall mechanical properties have been observed. However, how they relate to cell division remains an outstanding question. Here, we combine time-lapse imaging with local mechanical measurements via atomic force microscopy to systematically map the cell wall's age and growth, with their stiffness. We make use of two systems, Marchantia polymorpha gemmae, and Arabidopsis thaliana leaves. We first characterize the growth and cell division of M. polymorpha gemmae. We then demonstrate that cell division in M. polymorpha gemmae results in the generation of a temporary stiffer and slower-growing new wall. In contrast, this transient phenomenon is absent in A. thaliana leaves. We provide evidence that this different temporal behavior has a direct impact on the local cell geometry via changes in the junction angle. These results are expected to pave the way for developing more realistic plant morphogenetic models and to advance the study into the impact of cell division on tissue growth.

Perspectives on Scaffold Designs with Roles in Liver Cell Asymmetry and Medical and Industrial Applications by Using a New Type of Specialized 3D Bioprinter

Cellular asymmetry is an important element of efficiency in the compartmentalization of intracellular chemical reactions that ensure efficient tissue function. Improving the current 3D printing methods by using cellular asymmetry is essential in producing complex tissues and organs such as the liver. The use of cell spots containing at least two cells and basement membrane-like bio support materials allows cells to be tethered at two points on the basement membrane and with another cell in order to maintain cell asymmetry. Our model is a new type of 3D bioprinter that uses oriented multicellular complexes with cellular asymmetry. This novel approach is necessary to replace the sequential and slow processes of organogenesis with rapid methods of growth and 3D organ printing. The use of the extracellular matrix in the process of bioprinting with cells allows one to preserve the cellular asymmetry in the 3D printing process and thus preserve the compartmentalization of biological processes and metabolic efficiency.

Three-Dimensional Morphological Analysis Revealed the Cell Patterning Bases for the Sexual Dimorphism Development in the Liverwort Marchantia polymorpha

In land plants, sexual dimorphism can develop in both diploid sporophytes and haploid gametophytes. While developmental processes of sexual dimorphism have been extensively studied in the sporophytic reproductive organs of model flowering plants such as stamens and carpels of Arabidopsis thaliana, those occurring in gametophyte generation are less well characterized due to the lack of amenable model systems. In this study, we performed three-dimensional morphological analyses of gametophytic sexual branch differentiation in the liverwort Marchantia polymorpha, using high-depth confocal imaging and a computational cell segmentation technique. Our analysis revealed that the specification of germline precursors initiates in a very early stage of sexual branch development, where incipient branch primordia are barely recognizable in the apical notch region. Moreover, spatial distribution patterns of germline precursors differ between males and females from the initial stage of primordium development in a manner dependent on the master sexual differentiation regulator MpFGMYB. At later stages, distribution patterns of germline precursors predict the sex-specific gametangia arrangement and receptacle morphologies seen in mature sexual branches. Taken together, our data suggest a tightly coupled progression of germline segregation and sexual dimorphism development in M. polymorpha.

Interspersed expression of CUP-SHAPED COTYLEDON2 and REDUCED COMPLEXITY shapes Cardamine hirsuta complex leaf form

How genetically regulated growth shapes organ form is a key problem in developmental biology. Here, we investigate this problem using the leaflet-bearing complex leaves of Cardamine hirsuta as a model. Leaflet development requires the action of two growth-repressing transcription factors: REDUCED COMPLEXITY (RCO), a homeodomain protein, and CUP-SHAPED COTYLEDON2 (CUC2), a NAC-domain protein. However, how their respective growth-repressive actions are integrated in space and time to generate complex leaf forms remains unknown. By using live imaging, we show that CUC2 and RCO are expressed in an interspersed fashion along the leaf margin, creating a distinctive striped pattern. We find that this pattern is functionally important because forcing RCO expression in the CUC2 domain disrupts auxin-based marginal patterning and can abolish leaflet formation. By combining genetic perturbations with time-lapse imaging and cellular growth quantifications, we provide evidence that RCO-mediated growth repression occurs after auxin-based leaflet patterning and in association with the repression of cell proliferation. Additionally, through the use of genetic mosaics, we show that RCO is sufficient to repress both cellular growth and proliferation in a cell-autonomous manner. This mechanism of growth repression is different to that of CUC2, which occurs in proliferating cells. Our findings clarify how the two growth repressors RCO and CUC2 coordinate to subdivide developing leaf primordia into distinct leaflets and generate the complex leaf form. They also indicate different relationships between growth repression and cell proliferation in the patterning and post-patterning stages of organogenesis.

Tissues as networks of cells: towards generative rules of complex organ development

Network analysis is a well-known and powerful tool in molecular biology. More recently, it has been introduced in developmental biology. Tissues can be readily translated into spatial networks such that cells are represented by nodes and intercellular connections by edges. This discretization of cellular organization enables mathematical approaches rooted in network science to be applied towards the understanding of tissue structure and function. Here, we describe how such tissue abstractions can enable the principles that underpin tissue formation and function to be uncovered. We provide an introduction into biologically relevant network measures, then present an overview of different areas of developmental biology where these approaches have been applied. We then summarize the general developmental rules underpinning tissue topology generation. Finally, we discuss how generative models can help to link the developmental rule back to the tissue topologies. Our collection of results points at general mechanisms as to how local developmental rules can give rise to observed topological properties in multicellular systems.

Whole-mount smFISH allows combining RNA and protein quantification at cellular and subcellular resolution

Multicellular organisms result from complex developmental processes largely orchestrated through the quantitative spatiotemporal regulation of gene expression. Yet, obtaining absolute counts of messenger RNAs at a three-dimensional resolution remains challenging, especially in plants, owing to high levels of tissue autofluorescence that prevent the detection of diffraction-limited fluorescent spots. In situ hybridization methods based on amplification cycles have recently emerged, but they are laborious and often lead to quantification biases. In this article, we present a simple method based on single-molecule RNA fluorescence in situ hybridization to visualize and count the number of mRNA molecules in several intact plant tissues. In addition, with the use of fluorescent protein reporters, our method also enables simultaneous detection of mRNA and protein quantity, as well as subcellular distribution, in single cells. With this method, research in plants can now fully explore the benefits of the quantitative analysis of transcription and protein levels at cellular and subcellular resolution in plant tissues.

Brassinosteroid coordinates cell layer interactions in plants via cell wall and tissue mechanics

Growth coordination between cell layers is essential for development of most multicellular organisms. Coordination may be mediated by molecular signaling and/or mechanical connectivity between cells, but how genes modify mechanical interactions between layers is unknown. Here we show that genes driving brassinosteroid synthesis promote growth of internal tissue, at least in part, by reducing mechanical epidermal constraint. We identified a brassinosteroid-deficient dwarf mutant in the aquatic plant Utricularia gibba with twisted internal tissue, likely caused by mechanical constraint from a slow-growing epidermis. We tested this hypothesis by showing that a brassinosteroid mutant in Arabidopsis enhances epidermal crack formation, indicative of increased tissue stress. We propose that by remodeling cell walls, brassinosteroids reduce epidermal constraint, showing how genes can control growth coordination between layers by means of mechanics.

Grasses exploit geometry to achieve improved guard cell dynamics

Stomata are controllable micropores formed between two adjacent guard cells (GCs) that regulate gas flow across the plant surface.¹ Grasses, among the most successful organisms on the planet and the main food crops for humanity, have GCs flanked by specialized lateral subsidiary cells (SCs).^2,3,4 SCs improve performance by acting as a local pool of ions and metabolites to drive changes in turgor pressure within the GCs that open/close the stomatal pore.^4,5,6,7,8 The 4-celled complex also involves distinctive changes in geometry, having dumbbell-shaped GCs compared with typical kidney-shaped stomata.^2,4,9 However, the degree to which this distinctive geometry contributes to improved stomatal performance, and the underlying mechanism, remains unclear. To address this question, we created a finite element method (FEM) model of a grass stomatal complex that successfully captures experimentally observed pore opening/closure. Exploration of the model, including in silico and experimental mutant analyses, supports the importance of a reciprocal pressure system between GCs and SCs for effective stomatal function, with SCs functioning as springs to restrain lateral GC movement. Our results show that SCs are not essential but lead to a more responsive system. In addition, we show that GC wall anisotropy is not required for grass stomatal function (in contrast to kidney-shaped GCs¹⁰) but that a relatively thick GC rod region is needed to enhance pore opening. Our results demonstrate that a specific cellular geometry and associated mechanical properties are required for the effective functioning of grass stomata.

3D cellular morphometrics of ovule primordium development in Zea mays reveal differential division and growth dynamics specifying megaspore mother cell singleness

Introduction

Differentiation of spore mother cells marks the somatic-to-reproductive transition in higher plants. Spore mother cells are critical for fitness because they differentiate into gametes, leading to fertilization and seed formation. The female spore mother cell is called the megaspore mother cell (MMC) and is specified in the ovule primordium. The number of MMCs varies by species and genetic background, but in most cases, only a single mature MMC enters meiosis to form the embryo sac. Multiple candidate MMC precursor cells have been identified in both rice and Arabidopsis, so variability in MMC number is likely due to conserved early morphogenetic events. In Arabidopsis, the restriction of a single MMC per ovule, or MMC singleness, is determined by ovule geometry. To look for potential conservation of MMC ontogeny and specification mechanisms, we undertook a morphogenetic description of ovule primordium growth at cellular resolution in the model crop maize.

Methods

We generated a collection of 48 three-dimensional (3D) ovule primordium images for five developmental stages, annotated for 11 cell types. Quantitative analysis of ovule and cell morphological descriptors allowed the reconstruction of a plausible developmental trajectory of the MMC and its neighbors.

Results

The MMC is specified within a niche of enlarged, homogenous L2 cells, forming a pool of candidate archesporial (MMC progenitor) cells. A prevalent periclinal division of the uppermost central archesporial cell formed the apical MMC and the underlying cell, a presumptive stack cell. The MMC stopped dividing and expanded, acquiring an anisotropic, trapezoidal shape. By contrast, periclinal divisions continued in L2 neighbor cells, resulting in a single central MMC.

Discussion

We propose a model where anisotropic ovule growth in maize drives L2 divisions and MMC elongation, coupling ovule geometry with MMC fate.

Deep machine learning for cell segmentation and quantitative analysis of radial plant growth

Plants produce the major part of terrestrial biomass and are long-term deposits of atmospheric carbon. This capacity is to a large extent due to radial growth of woody species - a process driven by cambium stem cells located in distinct niches of shoot and root axes. In the model species Arabidopsis thaliana, thousands of cells are produced by the cambium in radial orientation generating a complex organ anatomy enabling long-distance transport, mechanical support and protection against biotic and abiotic stressors. These complex organ dynamics make a comprehensive and unbiased analysis of radial growth challenging and asks for tools for automated quantification. Here, we combined the recently developed PlantSeg and MorphographX image analysis tools, to characterize tissue morphogenesis of the Arabidopsis hypocotyl. After sequential training of segmentation models on ovules, shoot apical meristems and adult hypocotyls using deep machine learning, followed by the training of cell type classification models, our pipeline segments complex images of transverse hypocotyl sections with high accuracy and classifies central hypocotyl cell types. By applying our pipeline on both wild type and phloem intercalated with xylem (pxy) mutants, we also show that this strategy faithfully detects major anatomical aberrations. Collectively, we conclude that our established pipeline is a powerful phenotyping tool comprehensively extracting cellular parameters and providing access to tissue topology during radial plant growth.

Brassinosteroid gene regulatory networks at cellular resolution in the Arabidopsis root

Brassinosteroids are plant steroid hormones that regulate diverse processes, such as cell division and cell elongation, through gene regulatory networks that vary in space and time. By using time series single-cell RNA sequencing to profile brassinosteroid-responsive gene expression specific to different cell types and developmental stages of the Arabidopsis root, we identified the elongating cortex as a site where brassinosteroids trigger a shift from proliferation to elongation associated with increased expression of cell wall-related genes. Our analysis revealed HOMEOBOX FROM ARABIDOPSIS THALIANA 7 (HAT7) and GT-2-LIKE 1 (GTL1) as brassinosteroid-responsive transcription factors that regulate cortex cell elongation. These results establish the cortex as a site of brassinosteroid-mediated growth and unveil a brassinosteroid signaling network regulating the transition from proliferation to elongation, which illuminates aspects of spatiotemporal hormone responses.

An optimized pipeline for live imaging whole Arabidopsis leaves at cellular resolution

BACKGROUND

Live imaging is the gold standard for determining how cells give rise to organs. However, tracking many cells across whole organs over large developmental time windows is extremely challenging. In this work, we provide a comparably simple method for confocal live imaging entire Arabidopsis thaliana first leaves across early development. Our imaging method works for both wild-type leaves and the complex curved leaves of the jaw-1D mutant.

RESULTS

We find that dissecting the cotyledons, affixing a coverslip above the samples and mounting samples with perfluorodecalin yields optimal imaging series for robust cellular and organ level analysis. We provide details of our complementary image processing steps in MorphoGraphX software for segmenting, tracking lineages, and measuring a suite of cellular properties. We also provide MorphoGraphX image processing scripts we developed to automate analysis of segmented images and data presentation.

CONCLUSIONS

Our imaging techniques and processing steps combine into a robust imaging pipeline. With this pipeline we are able to examine important nuances in the cellular growth and differentiation of jaw-D versus WT leaves that have not been demonstrated before. Our pipeline is approachable and easy to use for leaf development live imaging.

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.

The cellular basis for synergy between RCO and KNOX1 homeobox genes in leaf shape diversity

Leaves of seed plants provide an attractive system to study the development and evolution of form. Leaves show varying degrees of margin complexity ranging from simple, as in Arabidopsis thaliana, to fully dissected into leaflets in the closely related species Cardamine hirsuta. Leaflet formation requires actions of Class I KNOTTED1-LIKE HOMEOBOX (KNOX1) and REDUCED COMPLEXITY (RCO) homeobox genes, which are expressed in the leaves of C. hirsuta but not A. thaliana. Evolutionary studies indicate that diversification of KNOX1 and RCO genes was repeatedly associated with increased leaf complexity. However, whether this gene combination represents a developmentally favored avenue for leaflet formation remains unknown, and the cell-level events through which the combined action of these genes drives leaflet formation are also poorly understood. Here we show, through a genetic screen, that when a C. hirsuta RCO transgene is expressed in A. thaliana, then ectopic KNOX1 expression in leaves represents a preferred developmental path for leaflet formation. Using time-lapse growth analysis, we demonstrate that KNOX1 expression in the basal domain of leaves leads to prolonged and anisotropic cell growth. This KNOX1 action, in synergy with local growth repression by RCO, is instrumental in generating rachises and petiolules, the linear geometrical elements, that bear leaflets in complex leaves. Our results show how the combination of cell-level growth analyses and genetics can help us understand how evolutionary modifications in expression of developmentally important genes are translated into diverse leaf shapes.

Mechanical conflict caused by a cell-wall-loosening enzyme activates de novo shoot regeneration

Cellular heterogeneity is a hallmark of multicellular organisms. During shoot regeneration from undifferentiated callus, only a select few cells, called progenitors, develop into shoot. How these cells are selected and what governs their subsequent progression to a patterned organ system is unknown. Using Arabidopsis thaliana, we show that it is not just the abundance of stem cell regulators but rather the localization pattern of polarity proteins that predicts the progenitor's fate. A shoot-promoting factor, CUC2, activated the expression of the cell-wall-loosening enzyme, XTH9, solely in a shell of cells surrounding the progenitor, causing different mechanical stresses in these cells. This mechanical conflict then activates cell polarity in progenitors to promote meristem formation. Interestingly, genetic or physical perturbations to cells surrounding the progenitor impaired the progenitor and vice versa. These suggest a feedback loop between progenitors and their neighbors for shoot regeneration in the absence of tissue-patterning cues.

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.

The annotation and analysis of complex 3D plant organs using 3DCoordX

A fundamental question in biology concerns how molecular and cellular processes become integrated during morphogenesis. In plants, characterization of 3D digital representations of organs at single-cell resolution represents a promising approach to addressing this problem. A major challenge is to provide organ-centric spatial context to cells of an organ. We developed several general rules for the annotation of cell position and embodied them in 3DCoordX, a user-interactive computer toolbox implemented in the open-source software MorphoGraphX. 3DCoordX enables rapid spatial annotation of cells even in highly curved biological shapes. Using 3DCoordX, we analyzed cellular growth patterns in organs of several species. For example, the data indicated the presence of a basal cell proliferation zone in the ovule primordium of Arabidopsis (Arabidopsis thaliana). Proof-of-concept analyses suggested a preferential increase in cell length associated with neck elongation in the archegonium of Marchantia (Marchantia polymorpha) and variations in cell volume linked to central morphogenetic features of a trap of the carnivorous plant Utricularia (Utricularia gibba). Our work demonstrates the broad applicability of the developed strategies as they provide organ-centric spatial context to cellular features in plant organs of diverse shape complexity.