Nonlinear compartmental modeling to monitor ovarian follicle population dynamics on the whole lifespan.

In this work, we introduce a compartmental model of ovarian follicle development all along lifespan, based on ordinary differential equations. The model predicts the changes in the follicle numbers in different maturation stages with aging. Ovarian follicles may either move forward to the next compartment (unidirectional migration) or degenerate and disappear (death). The migration from the first follicle compartment corresponds to the activation of quiescent follicles, which is responsible for the progressive exhaustion of the follicle reserve (ovarian aging) until cessation of reproductive activity. The model consists of a data-driven layer embedded into a more comprehensive, knowledge-driven layer encompassing the earliest events in follicle development. The data-driven layer is designed according to the most densely sampled experimental dataset available on follicle numbers in the mouse. Its salient feature is the nonlinear formulation of the activation rate, whose formulation includes a feedback term from growing follicles. The knowledge-based, coating layer accounts for cutting-edge studies on the initiation of follicle development around birth. Its salient feature is the co-existence of two follicle subpopulations of different embryonic origins. We then setup a complete estimation strategy, including the study of structural identifiability, the elaboration of a relevant optimization criterion combining different sources of data (the initial dataset on follicle numbers, together with data in conditions of perturbed activation, and data discriminating the subpopulations) with appropriate error models, and a model selection step. We finally illustrate the model potential for experimental design (suggestion of targeted new data acquisition) and in silico experiments.

Stochastic chemical kinetics of cell fate decision systems: From single cells to populations and back.

Stochastic chemical kinetics is a widely used formalism for studying stochasticity of chemical reactions inside single cells. Experimental studies of reaction networks are generally performed with cells that are part of a growing population, yet the population context is rarely taken into account when models are developed. Models that neglect the population context lose their validity whenever the studied system influences traits of cells that can be selected in the population, a property that naturally arises in the complex interplay between single-cell and population dynamics of cell fate decision systems. Here, we represent such systems as absorbing continuous-time Markov chains. We show that conditioning on non-absorption allows one to derive a modified master equation that tracks the time evolution of the expected population composition within a growing population. This allows us to derive consistent population dynamics models from a specification of the single-cell process. We use this approach to classify cell fate decision systems into two types that lead to different characteristic phases in emerging population dynamics. Subsequently, we deploy the gained insights to experimentally study a recurrent problem in biology: how to link plasmid copy number fluctuations and plasmid loss events inside single cells to growth of cell populations in dynamically changing environments.