A mathematical model for understanding synergistic regulations and paradoxical feedbacks in the shoot apical meristem

Intercellular Communication in Shoot Meristems

The shoot meristem of land plants maintains the capacity for organ generation throughout its lifespan due to a group of undifferentiated stem cells. Most meristems are shaped like a dome with a precise spatial arrangement of functional domains, and, within and between these domains, cells interact through a network of interconnected signaling pathways. Intercellular communication in meristems is mediated by mobile transcription factors, small RNAs, hormones, and secreted peptides that are perceived by membrane-localized receptors. In recent years, we have gained deeper insight into the underlying molecular processes of the shoot meristem, and we discuss here how plants integrate internal and external inputs to control shoot meristem activities.

An updated model of shoot apical meristem regulation by ERECTA family and CLAVATA3 signaling pathways in Arabidopsis

The shoot apical meristem (SAM) gives rise to the aboveground organs of plants. The size of the SAM is relatively constant due to the balance between stem cell replenishment and cell recruitment into new organs. In angiosperms, the transcription factor WUSCHEL (WUS) promotes stem cell proliferation in the central zone of the SAM. WUS forms a negative feedback loop with a signaling pathway activated by CLAVATA3 (CLV3). In the periphery of the SAM, the ERECTA family receptors (ERfs) constrain WUS and CLV3 expression. Here, we show that four ligands of ERfs redundantly inhibit the expression of these two genes. Transcriptome analysis confirmed that WUS and CLV3 are the main targets of ERf signaling and uncovered new ones. Analysis of promoter reporters indicated that the WUS expression domain mostly overlaps with the CLV3 domain and does not shift along the apical-basal axis in clv3 mutants. Our three-dimensional mathematical model captured gene expression distributions at the single-cell level under various perturbed conditions. Based on our findings, CLV3 regulates cellular levels of WUS mostly through autocrine signaling, and ERfs regulate the spatial expression of WUS, preventing its encroachment into the peripheral zone.

Initiation of aboveground organ primordia depends on combined action of auxin, ERECTA family genes, and PINOID

Leaves and flowers are produced by the shoot apical meristem (SAM) at a certain distance from its center, a process that requires the hormone auxin. The amount of auxin and the pattern of its distribution in the initiation zone determine the size and spatial arrangement of organ primordia. Auxin gradients in the SAM are formed by PIN-FORMED (PIN) auxin efflux carriers whose polar localization in the plasma membrane depends on the protein kinase PINOID (PID). Previous work determined that ERECTA (ER) family genes (ERfs) control initiation of leaves. ERfs are plasma membrane receptors that enable cell-to-cell communication by sensing extracellular small proteins from the EPIDERMAL PATTERNING FACTOR/EPF-LIKE (EPF/EPFL) family. Here, we investigated whether ERfs regulate initiation of organs by altering auxin distribution or signaling in Arabidopsis (Arabidopsis thaliana). Genetic and pharmacological data suggested that ERfs do not regulate organogenesis through PINs while transcriptomics data showed that ERfs do not alter primary transcriptional responses to auxin. Our results indicated that in the absence of ERf signaling the peripheral zone cells inefficiently initiate leaves in response to auxin signals and that increased accumulation of auxin in the er erecta-like1 (erl1) erl2 SAM can partially rescue organ initiation defects. We propose that both auxin and ERfs are essential for leaf initiation and that they have common downstream targets. Genetic data also indicated that the role of PID in initiation of cotyledons and leaves cannot be attributed solely to regulation of PIN polarity and PID is likely to have other functions in addition to regulation of auxin distribution.

Nonmodular oscillator and switch based on RNA decay drive regeneration of multimodal gene expression

Periodic gene expression dynamics are key to cell and organism physiology. Studies of oscillatory expression have focused on networks with intuitive regulatory negative feedback loops, leaving unknown whether other common biochemical reactions can produce oscillations. Oscillation and noise have been proposed to support mammalian progenitor cells’ capacity to restore heterogenous, multimodal expression from extreme subpopulations, but underlying networks and specific roles of noise remained elusive. We use mass-action-based models to show that regulated RNA degradation involving as few as two RNA species—applicable to nearly half of human protein-coding genes—can generate sustained oscillations without explicit feedback. Diverging oscillation periods synergize with noise to robustly restore cell populations’ bimodal expression on timescales of days. The global bifurcation organizing this divergence relies on an oscillator and bistable switch which cannot be decomposed into two structural modules. Our work reveals surprisingly rich dynamics of post-transcriptional reactions and a potentially widespread mechanism underlying development, tissue regeneration, and cancer cell heterogeneity.

Identification, visualization, statistical analysis and mathematical modeling of high-feedback loops in gene regulatory networks

Background

Feedback loops in gene regulatory networks play pivotal roles in governing functional dynamics of cells. Systems approaches demonstrated characteristic dynamical features, including multistability and oscillation, of positive and negative feedback loops. Recent experiments and theories have implicated highly interconnected feedback loops (high-feedback loops) in additional nonintuitive functions, such as controlling cell differentiation rate and multistep cell lineage progression. However, it remains challenging to identify and visualize high-feedback loops in complex gene regulatory networks due to the myriad of ways in which the loops can be combined. Furthermore, it is unclear whether the high-feedback loop structures with these potential functions are widespread in biological systems. Finally, it remains challenging to understand diverse dynamical features, such as high-order multistability and oscillation, generated by individual networks containing high-feedback loops. To address these problems, we developed HiLoop, a toolkit that enables discovery, visualization, and analysis of several types of high-feedback loops in large biological networks.

Results

HiLoop not only extracts high-feedback structures and visualize them in intuitive ways, but also quantifies the enrichment of overrepresented structures. Through random parameterization of mathematical models derived from target networks, HiLoop presents characteristic features of the underlying systems, including complex multistability and oscillations, in a unifying framework. Using HiLoop, we were able to analyze realistic gene regulatory networks containing dozens to hundreds of genes, and to identify many small high-feedback systems. We found more than a 100 human transcription factors involved in high-feedback loops that were not studied previously. In addition, HiLoop enabled the discovery of an enrichment of high feedback in pathways related to epithelial-mesenchymal transition.

Conclusions

HiLoop makes the study of complex networks accessible without significant computational demands. It can serve as a hypothesis generator through identification and modeling of high-feedback subnetworks, or as a quantification method for motif enrichment analysis. As an example of discovery, we found that multistep cell lineage progression may be driven by either specific instances of high-feedback loops with sparse appearances, or generally enriched topologies in gene regulatory networks. We expect HiLoop’s usefulness to increase as experimental data of regulatory networks accumulate. Code is freely available for use or extension at https://github.com/BenNordick/HiLoop.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04405-z.

ERECTA family signaling constrains CLAVATA3 and WUSCHEL to the center of the shoot apical meristem

The shoot apical meristem (SAM) is a reservoir of stem cells that gives rise to all post-embryonic above-ground plant organs. The size of the SAM remains stable over time owing to a precise balance of stem cell replenishment versus cell incorporation into organ primordia. The WUSCHEL (WUS)/CLAVATA (CLV) negative feedback loop is central to SAM size regulation. Its correct function depends on accurate spatial expression of WUS and CLV3 A signaling pathway, consisting of ERECTA family (ERf) receptors and EPIDERMAL PATTERNING FACTOR LIKE (EPFL) ligands, restricts SAM width and promotes leaf initiation. Although ERf receptors are expressed throughout the SAM, EPFL ligands are expressed in its periphery. Our genetic analysis of Arabidopsis demonstrated that ERfs and CLV3 synergistically regulate the size of the SAM, and wus is epistatic to ERf genes. Furthermore, activation of ERf signaling with exogenous EPFLs resulted in a rapid decrease of CLV3 and WUS expression. ERf-EPFL signaling inhibits expression of WUS and CLV3 in the periphery of the SAM, confining them to the center. These findings establish the molecular mechanism for stem cell positioning along the radial axis.

HAM Gene Family and Shoot Meristem Development

Land plants develop highly diversified shoot architectures, all of which are derived from the pluripotent stem cells in shoot apical meristems (SAMs). As sustainable resources for continuous organ formation in the aboveground tissues, SAMs play an important role in determining plant yield and biomass production. In this review, we summarize recent advances in understanding one group of key regulators - the HAIRY MERISTEM (HAM) family GRAS domain proteins - in shoot meristems. We highlight the functions of HAM family members in dictating shoot stem cell initiation and proliferation, the signaling cascade that shapes HAM expression domains in shoot meristems, and the conservation and diversification of HAM family members in land plants. We also discuss future directions that potentially lead to a more comprehensive view of the HAM gene family and stem cell homeostasis in land plants.