We study a mathematical model describing the dynamics of a pluripotent stem cell population involved in the blood production process in the bone marrow. This model is a differential equation with a time delay. The delay describes the cell cycle duration and is uniformly distributed on an interval. We obtain stability conditions independent of the delay. We also show that the distributed delay can destabilize the entire system. In particularly, it is shown that Hopf bifurcations can occur.
This paper is devoted to the analysis of a mathematical model of blood cells production in the bone marrow (hematopoiesis). The model is a system of two age-structured partial differential equations. Integrating these equations over the age, we obtain a system of two nonlinear differential equations with distributed time delay corresponding to the cell cycle duration. This system describes the evolution of the total cell populations. By constructing a Lyapunov functional, it is shown that the trivial equilibrium is globally asymptotically stable if it is the only equilibrium. It is also shown that the nontrivial equilibrium, the most biologically meaningful one, can become unstable via a Hopf bifurcation. Numerical simulations are carried out to illustrate the analytical results. The study maybe helpful in understanding the connection between the relatively short cell cycle durations and the relatively long periods of peripheral cell oscillations in some periodic hematological diseases.
How far is neuroepithelial cell proliferation in the developing central nervous system a deterministic process? Or, to put it in a more precise way, how accurately can it be described by a deterministic mathematical model? To provide tracks to answer this question, a deterministic system of transport and diffusion partial differential equations, both physiologically and spatially structured, is introduced as a model to describe the spatially organized process of cell proliferation during the development of the central nervous system. As an initial step towards dealing with the three-dimensional case, a unidimensional version of the model is presented. Numerical analysis and numerical tests are performed. In this work we also achieve a first experimental validation of the proposed model, by using cell proliferation data recorded from histological sections obtained during the development of the optic tectum in the chick embryo.
We introduce and analyze several aspects of a new model for cell differentiation. It assumes that differentiation of progenitor cells is a continuous process. From the mathematical point of view, it is based on partial differential equations of transport type. Specifically, it consists of a structured population equation with a nonlinear feedback loop. This models the signaling process due to cytokines, which regulate the differentiation and proliferation process. We compare the continuous model to its discrete counterpart, a multi-compartmental model of a discrete collection of cell subpopulations recently proposed by Marciniak-Czochra et al. in 2009 to investigate the dynamics of the hematopoietic system. We obtain uniform bounds for the solutions, characterize steady state solutions, and analyze their linearized stability. We show how persistence or extinction might occur according to values of parameters that characterize the stem cells self-renewal. We also perform numerical simulations and discuss the qualitative behavior of the continuous model vis a vis the discrete one.
Stem cells can precisely and robustly undergo cellular differentiation and lineage commitment, referred to as stemness. However, how the gene network underlying stemness regulation reliably specifies cell fates is not well understood. To address this question, we applied a recently developed computational method, Random Circuit Perturbation (RACIPE), to a nine-component gene regulatory network (GRN) governing stemness, from which we identified fifteen robust gene states. Among them, four out of the five most probable gene states exhibit gene expression patterns observed in single mouse embryonic cells at 32-cell and 64-cell stages. These gene states can be robustly predicted by the stemness GRN but not by randomiz
Collective cell migration is crucial in many biological processes such as wound healing, tissue morphogenesis, and tumor progression. The leading front of a collective migrating epithelial cell layer often destabilizes into multicellular finger-like protrusions, each of which is guided by a leader cell at the fingertip. Here, we develop a subcellular-element-based model of this fingering instability, which incorporates leader cells and other related properties of a monolayer of epithelial cells. Our model recovers multiple aspects of the dynamics, especially the traction force patterns and velocity fields, observed in experiments on MDCK cells. Our model predicts the necessity of the leader cell and its minimal functions for the formation and maintenance of a stable finger pattern. Meanwhile, our model allows for an analysis of the role of supra-cellular actin cable on the leading front, predicting that while this observed structure helps maintain the shape of the finger, it is not required in order to form a finger. In addition, we also study the driving instability in the context of continuum active fluid model, which justifies some of our assumptions in the computational approach. In particular, we show that in our model no finger protrusions would emerge in a phenotypically homogenous active fluid and hence the role of the leader cell and its followers are often critical.