-
Multiscale population dynamics in reproductive biology: singular perturbation reduction in deterministic and stochastic models
Authors:
Celine Bonnet,
Keltoum Chahour,
Frédérique Clément,
Marie Postel,
Romain Yvinec
Abstract:
In this study, we describe different modeling approaches for ovarian follicle population dynamics, based on either ordinary (ODE), partial (PDE) or stochastic (SDE) differential equations, and accounting for interactions between follicles. We put a special focus on representing the population-level feedback exerted by growing ovarian follicles onto the activation of quiescent follicles. We take ad…
▽ More
In this study, we describe different modeling approaches for ovarian follicle population dynamics, based on either ordinary (ODE), partial (PDE) or stochastic (SDE) differential equations, and accounting for interactions between follicles. We put a special focus on representing the population-level feedback exerted by growing ovarian follicles onto the activation of quiescent follicles. We take advantage of the timescale difference existing between the growth and activation processes to apply model reduction techniques in the framework of singular perturbations. We first study the linear versions of the models to derive theoretical results on the convergence to the limit models. In the nonlinear cases, we provide detailed numerical evidence of convergence to the limit behavior. We reproduce the main semi-quantitative features characterizing the ovarian follicle pool, namely a bimodal distribution of the whole population, and a slope break in the decay of the quiescent pool with aging.
△ Less
Submitted 4 March, 2019;
originally announced March 2019.
-
Stochastic nonlinear model for somatic cell population dynamics during ovarian follicle activation
Authors:
Frédérique Clément,
Frédérique Robin,
Romain Yvinec
Abstract:
In mammals, female germ cells are sheltered within somatic structures called ovarian follicles, which remain in a quiescent state until they get activated, all along reproductive life. We investigate the sequence of somatic cell events occurring just after follicle activation, starting by the awakening of precursor somatic cells, and their transformation into proliferative cells. We introduce a no…
▽ More
In mammals, female germ cells are sheltered within somatic structures called ovarian follicles, which remain in a quiescent state until they get activated, all along reproductive life. We investigate the sequence of somatic cell events occurring just after follicle activation, starting by the awakening of precursor somatic cells, and their transformation into proliferative cells. We introduce a nonlinear stochastic model accounting for the joint dynamics of the two cell types, and allowing us to investigate the potential impact of a feedback from proliferative cells onto precursor cells. To tackle the key issue of whether cell proliferation is concomitant or posterior to cell awakening, we assess both the time needed for all precursor cells to awake, and the corresponding increase in the total cell number with respect to the initial cell number. Using the probabilistic theory of first passage times, we design a numerical scheme based on a rigorous Finite State Projection and coupling techniques to compute the mean extinction time and the cell number at extinction time. We find that the feedback term clearly lowers the number of proliferative cells at the extinction time. We calibrate the model parameters using an exact likelihood approach. We carry out a comprehensive comparison between the initial model and a series of submodels, which helps to select the critical cell events taking place during activation, and suggests that awakening is prominent over proliferation.
△ Less
Submitted 9 December, 2020; v1 submitted 4 March, 2019;
originally announced March 2019.
-
Advances in computational modeling approaches in pituitary gonadotropin signaling
Authors:
Romain Yvinec,
Pascale Crépieux,
Eric Reiter,
Anne Poupon,
Frédérique Clément
Abstract:
Pituitary gonadotropins play an essential and pivotal role in the control of human and animal reproduction within the hypothalamic-pituitary-gonadal (HPG) axis. The computational modeling of pituitary gonadotropin signaling encompasses phenomena of different natures such as the dynamic encoding of gonadotropin secretion, and the intracellular cascades triggered by gonadotropin binding to their cog…
▽ More
Pituitary gonadotropins play an essential and pivotal role in the control of human and animal reproduction within the hypothalamic-pituitary-gonadal (HPG) axis. The computational modeling of pituitary gonadotropin signaling encompasses phenomena of different natures such as the dynamic encoding of gonadotropin secretion, and the intracellular cascades triggered by gonadotropin binding to their cognate receptors, resulting in a variety of biological outcomes. We overview historical and ongoing issues in modeling and data analysis related to gonadotropin secretion in the field of both physiology and neuro-endocrinology. We mention the different mathematical formalisms involved, their interest and limits. We discuss open statistical questions in signal analysis associated with key endocrine issues. We also review recent advances in the modeling of the intracellular pathways activated by gonadotropins, which yields promising development for innovative approaches in drug discovery. The greatest challenge to be tackled in computational modeling of pituitary gonadotropin signaling is the embedding of gonadotropin signaling within its natural multi-scale environment, from the single cell level, to the organic and whole HPG level. The development of modeling approaches of G protein-coupled receptor signaling, together with multicellular systems biology may lead to unexampled mechanistic understanding with critical expected fallouts in the therapeutic management of reproduction.
△ Less
Submitted 1 August, 2018;
originally announced August 2018.
-
Analysis and calibration of a linear model for structured cell populations with unidirectional motion : Application to the morphogenesis of ovarian follicles
Authors:
Frédérique Clément,
Frédérique Robin,
Romain Yvinec
Abstract:
We analyze a multi-type age dependent model for cell populations subject to unidirectional motion, in both a stochastic and deterministic framework. Cells are distributed into successive layers; they may divide and move irreversibly from one layer to the next. We adapt results on the large-time convergence of PDE systems and branching processes to our context, where the Perron-Frobenius or Krein-R…
▽ More
We analyze a multi-type age dependent model for cell populations subject to unidirectional motion, in both a stochastic and deterministic framework. Cells are distributed into successive layers; they may divide and move irreversibly from one layer to the next. We adapt results on the large-time convergence of PDE systems and branching processes to our context, where the Perron-Frobenius or Krein-Rutman theorem can not be applied. We derive explicit analytical formulas for the asymptotic cell number moments, and the stable age distribution. We illustrate these results numerically and we apply them to the study of the morphodynamics of ovarian follicles. We prove the structural parameter identifiability of our model in the case of age independent division rates. Using a set of experimental biological data, we estimate the model parameters to fit the changes in the cell numbers in each layer during the early stages of follicle development.
△ Less
Submitted 14 December, 2017;
originally announced December 2017.
-
DynPeak : An algorithm for pulse detection and frequency analysis in hormonal time series
Authors:
Alexandre Vidal,
Qinghua Zhang,
Claire Médigue,
Stéphane Fabre,
Frédérique Clément
Abstract:
The endocrine control of the reproductive function is often studied from the analysis of luteinizing hormone (LH) pulsatile secretion by the pituitary gland. Whereas measurements in the cavernous sinus cumulate anatomical and technical difficulties, LH levels can be easily assessed from jugular blood. However, plasma levels result from a convolution process due to clearance effects when LH enters…
▽ More
The endocrine control of the reproductive function is often studied from the analysis of luteinizing hormone (LH) pulsatile secretion by the pituitary gland. Whereas measurements in the cavernous sinus cumulate anatomical and technical difficulties, LH levels can be easily assessed from jugular blood. However, plasma levels result from a convolution process due to clearance effects when LH enters the general circulation. Simultaneous measurements comparing LH levels in the cavernous sinus and jugular blood have revealed clear differences in the pulse shape, the amplitude and the baseline. Besides, experimental sampling occurs at a relatively low frequency (typically every 10 min) with respect to LH highest frequency release (one pulse per hour) and the resulting LH measurements are noised by both experimental and assay errors. As a result, the pattern of plasma LH may be not so clearly pulsatile. Yet, reliable information on the InterPulse Intervals (IPI) is a prerequisite to study precisely the steroid feedback exerted on the pituitary level. Hence, there is a real need for robust IPI detection algorithms. In this article, we present an algorithm for the monitoring of LH pulse frequency, basing ourselves both on the available endocrinological knowledge on LH pulse (shape and duration with respect to the frequency regime) and synthetic LH data generated by a simple model. We make use of synthetic data to make clear some basic notions underlying our algorithmic choices. We focus on explaining how the process of sampling affects drastically the original pattern of secretion, and especially the amplitude of the detectable pulses. We then describe the algorithm in details and perform it on different sets of both synthetic and experimental LH time series. We further comment on how to diagnose possible outliers from the series of IPIs which is the main output of the algorithm.
△ Less
Submitted 22 December, 2011;
originally announced December 2011.
-
Endogenous circannual rhythm in LH secretion: insight from signal analysis coupled with mathematical modelling
Authors:
Alexandre Vidal,
Claire Médigue,
Benoit Malpaux,
Frédérique Clément
Abstract:
In sheep as in many vertebrates, the seasonal pattern of reproduction is timed by the annual photoperiodic cycle, characterized by seasonal changes in the daylength. The photoperiodic information is translated into a circadian profile of melatonin secretion. After multiple neuronal relays (within the hypothalamus), melatonin impacts GnRH (gonadotrophin releasing hormone) secretion that in turn c…
▽ More
In sheep as in many vertebrates, the seasonal pattern of reproduction is timed by the annual photoperiodic cycle, characterized by seasonal changes in the daylength. The photoperiodic information is translated into a circadian profile of melatonin secretion. After multiple neuronal relays (within the hypothalamus), melatonin impacts GnRH (gonadotrophin releasing hormone) secretion that in turn controls ovarian cyclicity. The pattern of GnRH secretion is mirrored into that of LH (luteinizing hormone) secretion, whose plasmatic level can be easily measured. We addressed the question of whether there exists an endogenous circannual rhythm in a tropical sheep population that exhibits clear seasonal ovarian activity when ewes are subjected to temperate latitudes. We based our analysis on LH time series collected in the course of 3 years from ewes subjected to a constant photoperiodic regime. Due to intra- and inter- animal variability and unequal sampling times, the existence of an endogenous rhythm is not straightforward. We have used time-frequency signal processing methods to extract hidden rhythms from the data. To further investigate the LF (low frequency) and HF (high frequency) components of the signals, we have designed a mathematical model of LH plasmatic level accounting for the effect of experimental sampling times. The model enables us to confirm the existence of an endogenous circannual rhythm, to investigate the action mechanism of photoperiod on the pulsatile pattern of LH secretion (control of the interpulse interval) and to conclude that the HF component is mainly due to the experimental sampling protocol.
△ Less
Submitted 10 March, 2009;
originally announced March 2009.
-
Bifurcation-based parameter tuning in a model of the GnRH pulse and surge generator
Authors:
Frédérique Clément,
Alexandre Vidal
Abstract:
We investigate a model of the GnRH pulse and surge generator, with the definite aim of constraining the model GnRH output with respect to a physiologically relevant list of specifications. The alternating pulse and surge pattern of secretion results from the interaction between a GnRH secreting system and a regulating system exhibiting fast-slow dynamics. The mechanisms underlying the behavior o…
▽ More
We investigate a model of the GnRH pulse and surge generator, with the definite aim of constraining the model GnRH output with respect to a physiologically relevant list of specifications. The alternating pulse and surge pattern of secretion results from the interaction between a GnRH secreting system and a regulating system exhibiting fast-slow dynamics. The mechanisms underlying the behavior of the model are reminded from the study of the Boundary-Layer System according to the "dissection method" principle. Using singular perturbation theory, we describe the sequence of bifurcations undergone by the regulating (FitzHugh-Nagumo) system, encompassing the rarely investigated case of homoclinic connexion. Basing on pure dynamical considerations, we restrict the space of parameter search for the regulating system and describe a foliation of this restricted space, whose leaves define constant duration ratios between the surge and the pulsatility phase in the whole system. We propose an algorithm to fix the parameter values to also meet the other prescribed ratios dealing with amplitude and frequency features of the secretion signal. We finally apply these results to illustrate the dynamics of GnRH secretion in the ovine species and the rhesus monkey.
△ Less
Submitted 10 September, 2008;
originally announced September 2008.
-
Multi-scale modeling of follicular ovulation as a reachability problem
Authors:
Nki Echenim,
Frederique Clement,
Michel Sorine
Abstract:
During each ovarian cycle, only a definite number of follicles ovulate, while the others undergo a degeneration process called atresia. We have designed a multi-scale mathematical model where ovulation and atresia result from a hormonal controlled selection process. A 2D-conservation law describes the age and maturity structuration of the follicular cell population. In this paper, we focus on th…
▽ More
During each ovarian cycle, only a definite number of follicles ovulate, while the others undergo a degeneration process called atresia. We have designed a multi-scale mathematical model where ovulation and atresia result from a hormonal controlled selection process. A 2D-conservation law describes the age and maturity structuration of the follicular cell population. In this paper, we focus on the operating mode of the control, through the study of the characteristics of the conservation law. We describe in particular the set of microscopic initial conditions leading to the macroscopic phenomenon of either ovulation or atresia, in the framework of backwards reachable sets theory.
△ Less
Submitted 6 February, 2007; v1 submitted 20 July, 2006;
originally announced July 2006.
-
Mathematical modeling of the GnRH pulse and surge generator
Authors:
Frederique Clement,
Jean-Pierre Francoise
Abstract:
We propose a mathematical model allowing for the alternating pulse and surge pattern of GnRH (Gonadotropin Releasing Hormone) secretion. The model is based on the coupling between two systems running on different time scales. The faster system corresponds to the average activity of GnRH neurons, while the slower one corresponds to the average activity of regulatory neurons. The analysis of the s…
▽ More
We propose a mathematical model allowing for the alternating pulse and surge pattern of GnRH (Gonadotropin Releasing Hormone) secretion. The model is based on the coupling between two systems running on different time scales. The faster system corresponds to the average activity of GnRH neurons, while the slower one corresponds to the average activity of regulatory neurons. The analysis of the slow/fast dynamics exhibited within and between both systems allows to explain the different patterns (slow oscillations, fast oscillations and periodical surge) of GnRH secretion. Specifications on the model parameter values are derived from physiological knowledge in terms of amplitude, frequency and plateau length of oscillations. The behavior of the model is finally illustrated by numerical simulations reproducing natural ovarian cycles and either direct or indirect actions of ovarian steroids on GnRH secretion.
△ Less
Submitted 27 October, 2006; v1 submitted 15 November, 2005;
originally announced November 2005.