DyCeModel: a tool for 1D simulation for distribution of plant hormones controlling tissue patterning

To study the mechanisms of growth and development, it is necessary to analyze the dynamics of the tissue patterning regulators in time and space and to take into account their effect on the cellular dynamics within a tissue. Plant hormones are the main regulators of the cell dynamics in plant tissues; they form gradients and maxima and control molecular processes in a concentration-dependent manner. Here, we present DyCeModel, a software tool implemented in MATLAB for one-dimensional simulation of tissue with a dynamic cellular ensemble, where changes in hormone (or other active substance) concentration in the cells are described by ordinary differential equations (ODEs). We applied DyCeModel to simulate cell dynamics in plant meristems with different cellular structures and demonstrated that DyCeModel helps to identify the relationships between hormone concentration and cellular behaviors. The tool visualizes the simulation progress and presents a video obtained during the calculation. Importantly, the tool is capable of automatically adjusting the parameters by fitting the distribution of the substance concentrations predicted in the model to experimental data taken from the microscopic images. Noteworthy, DyCeModel makes it possible to build models for distinct types of plant meristems with the same ODEs, recruiting specific input characteristics for each meristem. We demonstrate the tool's efficiency by simulation of the effect of auxin and cytokinin distributions on tissue patterning in two types of Arabidopsis thaliana stem cell niches: the root and shoot apical meristems. The resulting models represent a promising framework for further study of the role of hormone-controlled gene regulatory networks in cell dynamics.

CisCross: A gene list enrichment analysis to predict upstream regulators in Arabidopsis thaliana

Having DNA-binding profiles for a sufficient number of genome-encoded transcription factors (TFs) opens up the perspectives for systematic evaluation of the upstream regulators for the gene lists. Plant Cistrome database, a large collection of TF binding profiles detected using the DAP-seq method, made it possible for Arabidopsis. Here we re-processed raw DAP-seq data with MACS2, the most popular peak caller that leads among other ones according to quality metrics. In the benchmarking study, we confirmed that the improved collection of TF binding profiles supported a more precise gene list enrichment procedure, and resulted in a more relevant ranking of potential upstream regulators. Moreover, we consistently recovered the TF binding profiles that were missing in the previous collection of DAP-seq peak sets. We developed the CisCross web service (https://plamorph.sysbio.ru/ciscross/) that gives more flexibility in the analysis of potential upstream TF regulators for Arabidopsis thaliana genes.

Salicylic Acid in Root Growth and Development

In plants, salicylic acid (SA) is a hormone that mediates a plant’s defense against pathogens. SA also takes an active role in a plant’s response to various abiotic stresses, including chilling, drought, salinity, and heavy metals. In addition, in recent years, numerous studies have confirmed the important role of SA in plant morphogenesis. In this review, we summarize data on changes in root morphology following SA treatments under both normal and stress conditions. Finally, we provide evidence for the role of SA in maintaining the balance between stress responses and morphogenesis in plant development, and also for the presence of SA crosstalk with other plant hormones during this process.

3D Analysis of Mitosis Distribution Pattern in the Plant Root Tip with iRoCS Toolbox

The protocol allows to define and characterize mitosis distribution patterns in the plant root meristem. The method does not require genetic markers, which makes it applicable to plants of different non-transgenic genotypes, including ecotypes, mutants, and non-model plant species. Computer analysis of the mitosis distribution in three dimensions with iRoCS Toolbox identifies statistically significant changes in proliferation activity within specific root tissues and cell lineages.

Cell Dynamics in WOX5-Overexpressing Root Tips: The Impact of Local Auxin Biosynthesis

Root stem cell niche functioning requires the formation and maintenance of the specific "auxin-rich domain" governed by directional auxin transport and local auxin production. Auxin maximum co-localizes with the WOX5 expression domain in the quiescent center that separates mitotically active proximal and distal root meristems. Here we unravel the interconnected processes happening under WOX5 overexpression by combining in vivo experiments and mathematical modeling. We showed that WOX5-induced TAA1-mediated auxin biosynthesis is the cause, whereas auxin accumulation, PIN transporters relocation, and auxin redistribution between proximal and distal root meristems are its subsequent effects that influence the formation of the well-described phenotype with an enlarged root cap. These findings helped us to clarify the role of WOX5, which serves as a local QC-specific regulator that activates biosynthesis of non-cell-autonomous signal auxin to regulate the distal meristem functioning. The mathematical model with WOX5-mediated auxin biosynthesis and auxin-regulated cell growth, division, and detachment reproduces the columella cells dynamics in both wild type and under WOX5 dysregulation.

3D analysis of mitosis distribution highlights the longitudinal zonation and diarch symmetry in proliferation activity of the Arabidopsis thaliana root meristem

To date CYCB1;1 marker and cortex cell lengths have been conventionally used to determine the proliferation activity of the Arabidopsis root meristem. By creating a 3D map of mitosis distribution we showed that these markers overlooked that stele and endodermis save their potency to divide longer than the cortex and epidermis. Cessation of cell divisions is not a random process, so that mitotic activity within the endodermis and stele shows a diarch pattern. Mitotic activity of all root tissues peaked at the same distance from the quiescent center (QC); however, different tissues stopped dividing at different distances, with cells of the protophloem exiting the cell cycle first and the procambial cells being the last. The robust profile of mitotic activity in the root tip defines the longitudinal zonation in the meristem with the proliferation domain, where all cells are able to divide; and the transition domain, where the cell files cease to divide. 3D analysis of cytokinin deficient and cytokinin signaling mutants showed that their proliferation domain is similar to that of the wild type, but the transition domain is much longer. Our data suggest a strong inhibitory effect of cytokinin on anticlinal cell divisions in the stele.

[A model study of the role of proteins CLV1, CLV2, CLV3, and WUS in regulation of the structure of the shoot apical meristem]

In order to elucidate the role of proteins CLV1, CLV2, CLV3, and WUS in the mechanism underlying the maintenance of compartmental structure (spatial arrangement of the zones of biosynthesis of marker proteins) of the shoot apical meristem, a model of such mechanism was developed. Computational experiments led to biologically plausible solutions only when synthesis of substance W in a space between the organizing center and meristem apex was limited by the mechanism based on interaction of CLV3 with membrane receptor CLV1/CKV2 and lower boundary of the zone of W synthesis was determined by isoline of the corresponding threshold level of substance Y concentration. The model of the "reaction-diffusion" type formalizing the role proteins CLV1, CLV2, CLV3, and WUS can describe the basis of the mechanism underlying regulation of the compartmental structure of the shoot apical meristem and positioning of the organizing center in a certain site of the cell ensemble of such meristem.