Combined computational modeling and experimental analysis integrating chemical and mechanical signals suggests possible mechanism of shoot meristem maintenance.

Stem cell maintenance in multilayered shoot apical meristems (SAMs) of plants requires strict regulation of cell growth and division. Exactly how the complex milieu of chemical and mechanical signals interact in the central region of the SAM to regulate cell division plane orientation is not well understood. In this paper, simulations using a newly developed multiscale computational model are combined with experimental studies to suggest and test three hypothesized mechanisms for the regulation of cell division plane orientation and the direction of anisotropic cell expansion in the corpus. Simulations predict that in the Apical corpus, WUSCHEL and cytokinin regulate the direction of anisotropic cell expansion, and cells divide according to tensile stress on the cell wall. In the Basal corpus, model simulations suggest dual roles for WUSCHEL and cytokinin in regulating both the direction of anisotropic cell expansion and cell division plane orientation. Simulation results are followed by a detailed analysis of changes in cell characteristics upon manipulation of WUSCHEL and cytokinin in experiments that support model predictions. Moreover, simulations predict that this layer-specific mechanism maintains both the experimentally observed shape and structure of the SAM as well as the distribution of WUSCHEL in the tissue. This provides an additional link between the roles of WUSCHEL, cytokinin, and mechanical stress in regulating SAM growth and proper stem cell maintenance in the SAM.

Cytokinin stabilizes WUSCHEL by acting on the protein domains required for nuclear enrichment and transcription.

Concentration-dependent transcriptional regulation and the spatial regulation of transcription factor levels are poorly studied in plant development. WUSCHEL, a stem cell-promoting homeodomain transcription factor, accumulates at a higher level in the rib meristem than in the overlying central zone, which harbors stem cells in the shoot apical meristems of Arabidopsis thaliana. The differential accumulation of WUSCHEL in adjacent cells is critical for the spatial regulation and levels of CLAVATA3, a negative regulator of WUSCHEL transcription. Earlier studies have revealed that DNA-dependent dimerization, subcellular partitioning and protein destabilization control WUSCHEL protein levels and spatial accumulation. Moreover, the destabilization of WUSCHEL may also depend on the protein concentration. However, the roles of extrinsic spatial cues in maintaining differential accumulation of WUS are not understood. Through transient manipulation of hormone levels, hormone response patterns and analysis of the receptor mutants, we show that cytokinin signaling in the rib meristem acts through the transcriptional regulatory domains, the acidic domain and the WUSCHEL-box, to stabilize the WUS protein. Furthermore, we show that the same WUSCHEL-box functions as a degron sequence in cytokinin deficient regions in the central zone, leading to the destabilization of WUSCHEL. The coupled functions of the WUSCHEL-box in nuclear retention as described earlier, together with cytokinin sensing, reinforce higher nuclear accumulation of WUSCHEL in the rib meristem. In contrast a sub-threshold level may expose the WUSCHEL-box to destabilizing signals in the central zone. Thus, the cytokinin signaling acts as an asymmetric spatial cue in stabilizing the WUSCHEL protein to lead to its differential accumulation in neighboring cells, which is critical for concentration-dependent spatial regulation of CLAVATA3 transcription and meristem maintenance. Furthermore, our work shows that cytokinin response is regulated independently of the WUSCHEL function which may provide robustness to the regulation of WUSCHEL concentration.

Cell-Based Model of the Generation and Maintenance of the Shape and Structure of the Multilayered Shoot Apical Meristem of Arabidopsis thaliana.

One of the central problems in animal and plant developmental biology is deciphering how chemical and mechanical signals interact within a tissue to produce organs of defined size, shape, and function. Cell walls in plants impose a unique constraint on cell expansion since cells are under turgor pressure and do not move relative to one another. Cell wall extensibility and constantly changing distribution of stress on the wall are mechanical properties that vary between individual cells and contribute to rates of expansion and orientation of cell division. How exactly cell wall mechanical properties influence cell behavior is still largely unknown. To address this problem, a novel, subcellular element computational model of growth of stem cells within the multilayered shoot apical meristem (SAM) of Arabidopsis thaliana is developed and calibrated using experimental data. Novel features of the model include separate, detailed descriptions of cell wall extensibility and mechanical stiffness, deformation of the middle lamella, and increase in cytoplasmic pressure generating internal turgor pressure. The model is used to test novel hypothesized mechanisms of formation of the shape and structure of the growing, multilayered SAM based on WUS concentration of individual cells controlling cell growth rates and layer-dependent anisotropic mechanical properties of subcellular components of individual cells determining anisotropic cell expansion directions. Model simulations also provide a detailed prediction of distribution of stresses in the growing tissue which can be tested in future experiments.

DNA-dependent homodimerization, sub-cellular partitioning, and protein destabilization control WUSCHEL levels and spatial patterning.

The homeodomain transcription factor WUSCHEL (WUS) promotes stem cell maintenance in inflorescence meristems of Arabidopsis thaliana WUS, which is synthesized in the rib meristem, migrates and accumulates at lower levels in adjacent cells. Maintenance of WUS protein levels and spatial patterning distribution is not well-understood. Here, we show that the last 63-aa stretch of WUS is necessary for maintaining different levels of WUS protein in the rib meristem and adjacent cells. The 63-aa region contains the following transcriptional regulatory domains: the acidic region, the WUS-box, which is conserved in WUS-related HOMEOBOX family members, and the ethylene-responsive element binding factor-associated amphiphilic repression (EAR-like) domain. Our analysis reveals that the opposing functions of WUS-box, which is required for nuclear retention, and EAR-like domain, which participates in nuclear export, are necessary to maintain higher nuclear levels of WUS in cells of the rib meristem and lower nuclear levels in adjacent cells. We also show that the N-terminal DNA binding domain, which is required for both DNA binding and homodimerization, along with the homodimerization sequence located in the central part of the protein, restricts WUS from spreading excessively and show that the homodimerization is critical for WUS function. Our analysis also reveals that a higher level of WUS outside the rib meristem leads to protein destabilization, suggesting a new tier of regulation in WUS protein regulation. Taken together our data show that processes that influence WUS protein levels and spatial distribution are highly coupled to its transcriptional activity.

Threshold-dependent transcriptional discrimination underlies stem cell homeostasis.

Transcriptional mechanisms that underlie the dose-dependent regulation of gene expression in animal development have been studied extensively. However, the mechanisms of dose-dependent transcriptional regulation in plant development have not been understood. In Arabidopsis shoot apical meristems, WUSCHEL (WUS), a stem cell-promoting transcription factor, accumulates at a higher level in the rib meristem and at a lower level in the central zone where it activates its own negative regulator, CLAVATA3 (CLV3). How WUS regulates CLV3 levels has not been understood. Here we show that WUS binds a group of cis-elements, cis- regulatory module, in the CLV3-regulatory region, with different affinities and conformations, consisting of monomers at lower concentration and as dimers at a higher level. By deleting cis elements, manipulating the WUS-binding affinity and the homodimerization threshold of cis elements, and manipulating WUS levels, we show that the same cis elements mediate both the activation and repression of CLV3 at lower and higher WUS levels, respectively. The concentration-dependent transcriptional discrimination provides a mechanistic framework to explain the regulation of CLV3 levels that is critical for stem cell homeostasis.

WUSCHEL protein movement mediates stem cell homeostasis in the Arabidopsis shoot apex.

WUSCHEL (WUS) is a homeodomain transcription factor produced in cells of the niche/organizing center (OC) of shoot apical meristems. WUS specifies stem cell fate and also restricts its own levels by activating a negative regulator, CLAVATA3 (CLV3), in adjacent cells of the central zone (CZ). Here we show that the WUS protein, after being synthesized in cells of the OC, migrates into the CZ, where it activates CLV3 transcription by binding to its promoter elements. Using a computational model, we show that maintenance of the WUS gradient is essential to regulate stem cell number. Migration of a stem cell-inducing transcription factor into adjacent cells to activate a negative regulator, thereby restricting its own accumulation, is a theme that is unique to plant stem cell niches.

A high-resolution gene expression map of the Arabidopsis shoot meristem stem cell niche.

The shoot apical meristem (SAM) acts as a reservoir for stem cells. The central zone (CZ) harbors stem cells. The stem cell progenitors differentiate in the adjacent peripheral zone and in the rib meristem located just beneath the CZ. The SAM is further divided into distinct clonal layers: the L1 epidermal, L2 sub-epidermal and L3 layers. Collectively, SAMs are complex structures that consist of cells of different clonal origins that are organized into functional domains. By employing fluorescence-activated cell sorting, we have generated gene expression profiles of ten cell populations that belong to different clonal layers as well as domains along the central and peripheral axis. Our work reveals that cells in distinct clonal layers exhibit greater diversity in gene expression and greater transcriptional complexity than clonally related cell types in the central and peripheral axis. Assessment of molecular functions and biological processes reveals that epidermal cells express genes involved in pathogen defense: the L2 layer cells express genes involved in DNA repair pathways and telomere maintenance, and the L3 layers express transcripts involved in ion balance and salt tolerance besides photosynthesis. Strikingly, the stem cell-enriched transcriptome comprises very few hormone-responsive transcripts. In addition to providing insights into the expression profiles of hundreds of transcripts, the data presented here will act as a resource for reverse genetic analysis and will be useful in deciphering molecular pathways involved in cell type specification and their functions.

Plant stem cell maintenance involves direct transcriptional repression of differentiation program.

In animal systems, master regulatory transcription factors (TFs) mediate stem cell maintenance through a direct transcriptional repression of differentiation promoting TFs. Whether similar mechanisms operate in plants is not known. In plants, shoot apical meristems serve as reservoirs of stem cells that provide cells for all above ground organs. WUSCHEL, a homeodomain TF produced in cells of the niche, migrates into adjacent cells where it specifies stem cells. Through high-resolution genomic analysis, we show that WUSCHEL represses a large number of genes that are expressed in differentiating cells including a group of differentiation promoting TFs involved in leaf development. We show that WUS directly binds to the regulatory regions of differentiation promoting TFs; KANADI1, KANADI2, ASYMMETRICLEAVES2 and YABBY3 to repress their expression. Predictions from a computational model, supported by live imaging, reveal that WUS-mediated repression prevents premature differentiation of stem cell progenitors, being part of a minimal regulatory network for meristem maintenance. Our work shows that direct transcriptional repression of differentiation promoting TFs is an evolutionarily conserved logic for stem cell regulation.

WUSCHEL mediates stem cell homeostasis by regulating stem cell number and patterns of cell division and differentiation of stem cell progenitors.

Plant stem cell populations, unlike their animal counterparts, do not use cell migration and oriented cell divisions to maintain their size, and therefore require a precise coordination between self-renewing divisions of stem cells, and rates of cell division and differentiation among stem cell progenitors. Shoot apical meristems (SAMs) of higher plants harbor a set of stem cells within the central zone (CZ) that divide infrequently. Stem cell daughters that are displaced towards the surrounding peripheral zone (PZ) divide at a faster rate and enter into differentiation at specific locations to form leaves or flowers. The relative ratios of cells in the CZ and the PZ are maintained, despite a constant displacement of cells from the CZ into the PZ, and subsequent allocation of cells within the PZ to form organ primordia. The mechanisms that mediate this homeostatic balance are not well understood. A homeodomain transcription factor WUSCHEL, expressed in the rib meristem (RM), located beneath the CZ, has been shown to provide nonautonomous cues for stem cell specification. By employing transient spatial manipulation and live imaging, we show that an elevated level of WUS not only induces expansion of the CZ, but also results in increased cell division rates in cells of the PZ; conversely, decreases in WUS level lead to a smaller CZ and are associated with a reduction in cell division rate. Moreover, low levels of WUS lead to enlarged organ primordia, by elevating the responsiveness of the PZ cells to the plant hormone auxin. This reveals a function of WUS in mediating the balance between differentiating and non-differentiating cells of the PZ. Regulation of stem cell numbers, growth and differentiation patterns by a single transcription factor forms a interconnected and self-correcting feedback loop to provide robustness to stem cell homeostasis in a dynamic cellular environment.

Gene expression map of the Arabidopsis shoot apical meristem stem cell niche.

Despite the central importance of stem cells in plant growth and development, the molecular signatures associated with them have not been revealed. Shoot apical meristems (SAMs) harbor a small set of stem cells located at the tip of each plant and they are surrounded by several million differentiating cells. This imposes a major limitation in isolating pure populations of stem cells for genomic analyses. We have developed a system to isolate pure populations of distinct cell types of the SAMs, including stem cells. We have used this system to profile gene expression from 4 different cell samples of SAMs. The cell sample-specific gene expression profiling has resulted in a high-resolution gene expression map to reveal gene expression networks specific to individual spatial domains of SAMs. We demonstrate that the cell sample-specific expression profiling is sensitive in identifying rare transcripts expressed in a few specific subsets of cells of SAMs. Our extensive RNA in situ analysis reveals that the expression map can be used as a predictive tool in analyzing the spatial expression patterns of genes and it has led to the identification of unique gene expression patterns within the SAMs. Furthermore, our work reveals an enrichment of DNA repair and chromatin modification pathways in stem cells suggesting that maintenance of genome stability and flexible chromatin may be crucial for stem cell function. The gene expression map should guide future reverse genetics experiments, high-resolution analyses of cell-cell communication networks and epigenetic modifications.

Concentration-dependent transcriptional switching through a collective action of cis-elements.

Gene expression specificity of homeobox transcription factors has remained paradoxical. WUSCHEL activates and represses CLAVATA3 transcription at lower and higher concentrations, respectively. We use computational modeling and experimental analysis to investigate the properties of the cis-regulatory module. We find that intrinsically each cis-element can only activate CLAVATA3 at a higher WUSCHEL concentration. However, together, they repress CLAVATA3 at higher WUSCHEL and activate only at lower WUSCHEL, showing that the concentration-dependent interactions among cis-elements regulate both activation and repression. Biochemical experiments show that two adjacent functional cis-elements bind WUSCHEL with higher affinity and dimerize at relatively lower levels. Moreover, increasing the distance between cis-elements prolongs WUSCHEL monomer binding window, resulting in higher CLAVATA3 activation. Our work showing a constellation of optimally spaced cis-elements of defined affinities determining activation and repression thresholds in regulating CLAVATA3 transcription provides a previously unknown mechanism of cofactor-independent regulation of transcription factor binding in mediating gene expression specificity.

Automated tracking of stem cell lineages of Arabidopsis shoot apex using local graph matching.

Shoot apical meristems (SAMs) of higher plants harbor stem-cell niches. The cells of the stem-cell niche are organized into spatial domains of distinct function and cell behaviors. A coordinated interplay between cell growth dynamics and changes in gene expression is critical to ensure stem-cell homeostasis and organ differentiation. Exploring the causal relationships between cell growth patterns and gene expression dynamics requires quantitative methods to analyze cell behaviors from time-lapse imagery. Although technical breakthroughs in live-imaging methods have revealed spatio-temporal dynamics of SAM-cell growth patterns, robust computational methods for cell segmentation and automated tracking of cells have not been developed. Here we present a local graph matching-based method for automated-tracking of cells and cell divisions of SAMs of Arabidopsis thaliana. The cells of the SAM are tightly clustered in space which poses a unique challenge in computing spatio-temporal correspondences of cells. The local graph-matching principle efficiently exploits the geometric structure and topology of the relative positions of cells in obtaining spatio-temporal correspondences. The tracker integrates information across multiple slices in which a cell may be properly imaged, thus providing robustness to cell tracking in noisy live-imaging datasets. By relying on the local geometry and topology, the method is able to track cells in areas of high curvature such as regions of primordial outgrowth. The cell tracker not only computes the correspondences of cells across spatio-temporal scale, but it also detects cell division events, and identifies daughter cells upon divisions, thus allowing automated estimation of cell lineages from images captured over a period of 72 h. The method presented here should enable quantitative analysis of cell growth patterns and thus facilitating the development of in silico models for SAM growth.

CLAVATA3 mediated simultaneous control of transcriptional and post-translational processes provides robustness to the WUSCHEL gradient.

Regulation of the homeodomain transcription factor WUSCHEL concentration is critical for stem cell homeostasis in Arabidopsis shoot apical meristems. WUSCHEL regulates the transcription of CLAVATA3 through a concentration-dependent activation-repression switch. CLAVATA3, a secreted peptide, activates receptor kinase signaling to repress WUSCHEL transcription. Considering the revised regulation, CLAVATA3 mediated repression of WUSCHEL transcription alone will lead to an unstable system. Here we show that CLAVATA3 signaling regulates nuclear-cytoplasmic partitioning of WUSCHEL to control nuclear levels and its diffusion into adjacent cells. Our work also reveals that WUSCHEL directly interacts with EXPORTINS via EAR-like domain which is also required for destabilizing WUSCHEL in the cytoplasm. We develop a combined experimental and computational modeling approach that integrates CLAVATA3-mediated transcriptional repression of WUSCHEL and post-translational control of nuclear levels with the WUSCHEL concentration-dependent regulation of CLAVATA3. We show that the dual control by the same signal forms a seamless connection between de novo WUSCHEL synthesis and sub-cellular partitioning in providing robustness to the WUSCHEL gradient.

Live-imaging stem-cell homeostasis in the Arabidopsis shoot apex.

A precise spatio-temporal regulation of growth and differentiation is crucial to maintain a stable population of stem cells in the shoot apical meristems (SAMs) of higher plants. The real-time and simultaneous observations of dynamics of cell identity transitions, growth patterns, and signaling machinery involved in cell-cell communication is crucial to gain a mechanistic view of stem-cell homeostasis. In this article, I review recent advances in understanding the regulatory dynamics of stem-cell maintenance in Arabidopsis thaliana and discuss future challenges involved in transforming the static maps of genetic interactions into a dynamic framework representing functional molecular and cellular interactions in living SAMs.

Fluorescence activated cell sorting of shoot apical meristem cell types.

Growing tips of plants harbor a set of stem cells in structures called shoot apical meristems (SAMs) which provide cells for development of aboveground biomass. Despite a periodic differentiation of stem cell progenitors into leaves, the stem cell pool remains constant over time. Genetic analysis has revealed molecular pathways involved in stem-cell specification, cell division patterns, and organ differentiation. Stem cells within SAMs are few in number, which imposes a limitation to the experimental approaches that can be used for deciphering the gene regulatory networks that underlie cell fate transitions. Here, I provide detailed experimental protocols for the protoplasting and subsequent purification through cell sorting of SAM cells, which allows genome-wide analyses of gene expression patterns at a single cell-type resolution.

Gene expression analysis of shoot apical meristem cell types.

Shoot apical meristems (SAMs) of higher plants harbor a set of stem-cells and provide cells for the development of all the above-ground biomass of plants. Most of the important pattern formation events such as maintenance of stem-cell identity, specification and differentiation of leaf/flower primordia, and temporal control of the transition from vegetative to reproductive program are determined in SAMs. Genetic analysis has revealed molecular and hormonal pathways involved in stem-cell maintenance, organ differentiation, and flowering time. However, limited information is available as to how different pathways interact with each other to function as a network in specifying different cell types and their function. Deciphering gene networks that underlie cell fate transitions requires new approaches aimed at assaying genome-scale expression patterns of genes at a single cell-type resolution. Here we provide details of experimental methods involved in protoplasting of SAM cells, generating cell type-specific gene expression profiles, and analysis platforms for identifying and inferring gene networks.

Adaptive geometric tessellation for 3D reconstruction of anisotropically developing cells in multilayer tissues from sparse volumetric microscopy images.

The need for quantification of cell growth patterns in a multilayer, multi-cellular tissue necessitates the development of a 3D reconstruction technique that can estimate 3D shapes and sizes of individual cells from Confocal Microscopy (CLSM) image slices. However, the current methods of 3D reconstruction using CLSM imaging require large number of image slices per cell. But, in case of Live Cell Imaging of an actively developing tissue, large depth resolution is not feasible in order to avoid damage to cells from prolonged exposure to laser radiation. In the present work, we have proposed an anisotropic Voronoi tessellation based 3D reconstruction framework for a tightly packed multilayer tissue with extreme z-sparsity (2-4 slices/cell) and wide range of cell shapes and sizes. The proposed method, named as the 'Adaptive Quadratic Voronoi Tessellation' (AQVT), is capable of handling both the sparsity problem and the non-uniformity in cell shapes by estimating the tessellation parameters for each cell from the sparse data-points on its boundaries. We have tested the proposed 3D reconstruction method on time-lapse CLSM image stacks of the Arabidopsis Shoot Apical Meristem (SAM) and have shown that the AQVT based reconstruction method can correctly estimate the 3D shapes of a large number of SAM cells.

Quantitative analysis of live-cell growth at the shoot apex of Arabidopsis thaliana: algorithms for feature measurement and temporal alignment.

Study of the molecular control of organ growth requires establishment of the causal relationship between gene expression and cell behaviors. We seek to understand this relationship at the shoot apical meristem (SAM) of model plant Arabidopsis thaliana. This requires the spatial mapping and temporal alignment of different functional domains into a single template. Live-cell imaging techniques allow us to observe real-time organ primordia growth and gene expression dynamics at cellular resolution. In this paper, we propose a framework for the measurement of growth features at the 3D reconstructed surface of organ primordia, as well as algorithms for robust time alignment of primordia. We computed areas and deformation values from reconstructed 3D surfaces of individual primordia from live-cell imaging data. Based on these growth measurements, we applied a multiple feature landscape matching (LAM-M) algorithm to ensure a reliable temporal alignment of multiple primordia. Although the original landscape matching (LAM) algorithm motivated our alignment approach, it sometimes fails to properly align growth curves in the presence of high noise/distortion. To overcome this shortcoming, we modified the cost function to consider the landscape of the corresponding growth features. We also present an alternate parameter-free growth alignment algorithm which performs as well as LAM-M for high-quality data, but is more robust to the presence of outliers or noise. Results on primordia and guppy evolutionary growth data show that the proposed alignment framework performs at least as well as the LAM algorithm in the general case, and significantly better in the case of increased noise.

Cytokinin signaling regulates pavement cell morphogenesis in Arabidopsis.

The puzzle piece-shaped Arabidopsis leaf pavement cells (PCs) with interdigitated lobes and indents is a good model system to investigate the mechanisms that coordinate cell polarity and shape formation within a tissue. Auxin has been shown to coordinate the interdigitation by activating ROP GTPase-dependent signaling pathways. To identify additional components or mechanisms, we screened for mutants with abnormal PC morphogenesis and found that cytokinin signaling regulates the PC interdigitation pattern. Reduction in cytokinin accumulation and defects in cytokinin signaling (such as in ARR7-over-expressing lines, the ahk3cre1 cytokinin receptor mutant, and the ahp12345 cytokinin signaling mutant) enhanced PC interdigitation, whereas over-production of cytokinin and over-activation of cytokinin signaling in an ARR20 over-expression line delayed or abolished PC interdigitation throughout the cotyledon. Genetic and biochemical analyses suggest that cytokinin signaling acts upstream of ROPs to suppress the formation of interdigitated pattern. Our results provide novel mechanistic understanding of the pathways controlling PC shape and uncover a new role for cytokinin signaling in cell morphogenesis.