Vijay Chickarmane

Multiple feedback loops through cytokinin signaling control stem cell number within the Arabidopsis shoot meristem

A central unanswered question in stem cell biology, both in plants and in animals, is how the spatial organization of stem cell niches are maintained as cells move through them. We address this question for the shoot apical meristem (SAM) which harbors pluripotent stem cells responsible for growth of above-ground tissues in flowering plants. We find that localized perception of the plant hormone cytokinin establishes a spatial domain in which cell fate is respecified through induction of the master regulator WUSCHEL as cells are displaced during growth. Cytokinin-induced WUSCHEL expression occurs through both CLAVATA-dependent and CLAVATA-independent pathways. Computational analysis shows that feedback between cytokinin response and genetic regulators predicts their relative patterning, which we confirm experimentally. Our results also may explain how increasing cytokinin concentration leads to the first steps in reestablishing the shoot stem cell niche in vitro.

Transcriptional dynamics of the embryonic stem cell switch

Recent ChIP experiments of human and mouse embryonic stem cells have elucidated the architecture of the transcriptional regulatory circuitry responsible for cell determination, which involves the transcription factors OCT4, SOX2, and NANOG. In addition to regulating each other through feedback loops, these genes also regulate downstream target genes involved in the maintenance and differentiation of embryonic stem cells. A search for the OCT4–SOX2–NANOG network motif in other species reveals that it is unique to mammals. With a kinetic modeling approach, we ascribe function to the observed OCT4–SOX2–NANOG network by making plausible assumptions about the interactions between the transcription factors at the gene promoter binding sites and RNA polymerase (RNAP), at each of the three genes as well as at the target genes. We identify a bistable switch in the network, which arises due to several positive feedback loops, and is switched on/off by input environmental signals. The switch stabilizes the expression levels of the three genes, and through their regulatory roles on the downstream target genes, leads to a binary decision: when OCT4, SOX2, and NANOG are expressed and the switch is on, the self-renewal genes are on and the differentiation genes are off. The opposite holds when the switch is off. The model is extremely robust to parameter changes. In addition to providing a self-consistent picture of the transcriptional circuit, the model generates several predictions. Increasing the binding strength of NANOG to OCT4 and SOX2, or increasing its basal transcriptional rate, leads to an irreversible bistable switch: the switch remains on even when the activating signal is removed. Hence, the stem cell can be manipulated to be self-renewing without the requirement of input signals. We also suggest tests that could discriminate between a variety of feedforward regulation architectures of the target genes by OCT4, SOX2, and NANOG.

Variability in the Control of Cell Division Underlies Sepal Epidermal Patterning in Arabidopsis thaliana

Live cell imaging and computational modeling explains how variability in the timing of cell division generates a characteristic pattern of cell sizes during development.

Cytokinin signaling as a positional cue for patterning the apical–basal axis of the growing Arabidopsis shoot meristem

The transcription factor WUSCHEL (WUS) acts from a well-defined domain within the Arabidopsis thaliana shoot apical meristem (SAM) to maintain a stem cell niche. A negative-feedback loop involving the CLAVATA (CLV) signaling pathway regulates the number of WUS-expressing cells and provides the current paradigm for the homeostatic maintenance of stem cell numbers. Despite the continual turnover of cells in the SAM during development, the WUS domain remains patterned at a fixed distance below the shoot apex. Recent work has uncovered a positive-feedback loop between WUS function and the plant hormone cytokinin. Furthermore, loss of function of the cytokinin biosynthetic gene, LONELY GUY (LOG), results in a wus-like phenotype in rice. Herein, we find the Arabidopsis LOG4 gene is expressed in the SAM epidermis. We use this to develop a computational model representing a growing SAM to suggest the plausibility that apically derived cytokinin and CLV signaling, together, act as positional cues for patterning the WUS domain within the stem cell niche. Furthermore, model simulations backed by experimental data suggest a previously unknown negative feedback between WUS function and cytokinin biosynthesis in the Arabidopsis SAM epidermis. These results suggest a plausible dynamic feedback principle by which the SAM stem cell niche is patterned.

A Computational Model for Understanding Stem Cell, Trophectoderm and Endoderm Lineage Determination

Background Recent studies have associated the transcription factors, Oct4, Sox2 and Nanog as parts of a self-regulating network which is responsible for maintaining embryonic stem cell properties: self renewal and pluripotency. In addition, mutual antagonism between two of these and other master regulators have been shown to regulate lineage determination. In particular, an excess of Cdx2 over Oct4 determines the trophectoderm lineage whereas an excess of Gata-6 over Nanog determines differentiation into the endoderm lineage. Also, under/over-expression studies of the master regulator Oct4 have revealed that some self-renewal/pluripotency as well as differentiation genes are expressed in a biphasic manner with respect to the concentration of Oct4. Methodology/Principal Findings We construct a dynamical model of a minimalistic network, extracted from ChIP-on-chip and microarray data as well as literature studies. The model is based upon differential equations and makes two plausible assumptions; activation of Gata-6 by Oct4 and repression of Nanog by an Oct4–Gata-6 heterodimer. With these assumptions, the results of simulations successfully describe the biphasic behavior as well as lineage commitment. The model also predicts that reprogramming the network from a differentiated state, in particular the endoderm state, into a stem cell state, is best achieved by over-expressing Nanog, rather than by suppression of differentiation genes such as Gata-6. Conclusions The computational model provides a mechanistic understanding of how different lineages arise from the dynamics of the underlying regulatory network. It provides a framework to explore strategies of reprogramming a cell from a differentiated state to a stem cell state through directed perturbations. Such an approach is highly relevant to regenerative medicine since it allows for a rapid search over the host of possibilities for reprogramming to a stem cell state.

Computational modeling of the hematopoietic erythroid-myeloid switch reveals insights into cooperativity, priming, and irreversibility

Hematopoietic stem cell lineage choices are decided by genetic networks that are turned ON/OFF in a switch-like manner. However, prior to lineage commitment, genes are primed at low expression levels. Understanding the underlying molecular circuitry in terms of how it governs both a primed state and, at the other extreme, a committed state is of relevance not only to hematopoiesis but also to developmental systems in general. We develop a computational model for the hematopoietic erythroid-myeloid lineage decision, which is determined by a genetic switch involving the genes PU.1 and GATA-1. Dynamical models based upon known interactions between these master genes, such as mutual antagonism and autoregulation, fail to make the system bistable, a desired feature for robust lineage determination. We therefore suggest a new mechanism involving a cofactor that is regulated as well as recruited by one of the master genes to bind to the antagonistic partner that is necessary for bistability and hence switch-like behavior. An interesting fallout from this architecture is that suppression of the cofactor through external means can lead to a loss of cooperativity, and hence to a primed state for PU.1 and GATA-1. The PU.1–GATA-1 switch also interacts with another mutually antagonistic pair, –FOG-1. The latter pair inherits the state of its upstream master genes and further reinforces the decision due to several feedback loops, thereby leading to irreversible commitment. The genetic switch, which handles the erythroid-myeloid lineage decision, is an example of a network that implements both a primed and a committed state by regulating cooperativity through recruitment of cofactors. Perturbing the feedback between the master regulators and downstream targets suggests potential reprogramming strategies. The approach points to a framework for lineage commitment studies in general and could aid the search for lineage-determining genes.

Conservation analysis of large biochemical networks

MOTIVATION Large biochemical networks pose a unique challenge from the point of view of evaluating conservation laws. The computational problem in most cases exceeds the capability of available software tools, often resulting in inaccurate computation of the number and form of conserved cycles. Such errors have profound effects on subsequent calculations, particularly in the evaluation of the Jacobian which is a critical quantity in many other calculations. The goal of this paper is to outline a new algorithm that is computationally efficient and robust at extracting the correct conservation laws for very large biochemical networks. RESULTS We show that our algorithm can perform the conservation analysis of large biochemical networks, and can evaluate the correct conserved cycles when compared with other similar software tools. Biochemical simulators such as Jarnac and COPASI are successful at extracting only a subset of the conservation laws that our algorithm can. This is illustrated with examples for some large networks which show the advantages of our method.

Computational morphodynamics of plants: integrating development over space and time

The emerging field of computational morphodynamics aims to understand the changes that occur in space and time during development by combining three technical strategies: live imaging to observe development as it happens; image processing and analysis to extract quantitative information; and computational modelling to express and test time-dependent hypotheses. The strength of the field comes from the iterative and combined use of these techniques, which has provided important insights into plant development.

Bifurcation discovery tool

MOTIVATION Biochemical networks often yield interesting behavior such as switching, oscillation and chaotic dynamics. This article describes a tool that is capable of searching for bifurcation points in arbitrary ODE-based reaction networks by directing the user to regions in the parameter space, where such interesting dynamical behavior can be observed. RESULTS We have implemented a genetic algorithm that searches for Hopf bifurcations, turning points and bistable switches. The software is implemented as a Systems Biology Workbench (SBW) enabled module and accepts the standard SBML model format. The interface permits a user to choose the parameters to be searched, admissible parameter ranges, and the nature of the bifurcation to be sought. The tool will return the parameter values for the model for which the particular behavior is observed. AVAILABILITY The software, tutorial manual and test models are available for download at the following website: http:/www.sys-bio.org/ under the bifurcation link. The software is an open source and licensed under BSD.

Computational Morphodynamics: A Modeling Framework to Understand Plant Growth

Computational morphodynamics utilizes computer modeling to understand the development of living organisms over space and time. Results from biological experiments are used to construct accurate and predictive models of growth. These models are then used to make novel predictions that provide further insight into the processes involved, which can be tested experimentally to either confirm or rule out the validity of the computational models. This review highlights two fundamental challenges: (a) to understand the feedback between mechanics of growth and chemical or molecular signaling, and (b) to design models that span and integrate single cell behavior with tissue development. We review different approaches to model plant growth and discuss a variety of model types that can be implemented to demonstrate how the interplay between computational modeling and experimentation can be used to explore the morphodynamics of plant development.