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.

Le manuscrit disponible sur cette page "Clement_Vidal.ps" est une version pré-printWe 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 ofInternational audienceWe 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

In this work we study mixed mode oscillations in a model of secretion of GnRH (Gonadotropin Releasing Hormone). The model is a phantom burster consisting of two feedforward coupled FitzHugh-Nagumo systems, with three time scales. The forcing system (Regulator) evolves on the slowest scale and acts by moving the slow null-cline of the forced system (Secretor). There are three modes of dynamics: pulsatility (transient relaxation oscillation), surge (quasi steady state) and small oscillations related to the passage of the slow null-cline through a fold point of the fast null-cline. We derive a variety of reductions, taking advantage of the mentioned features of the system. We obtain two results; one on the local dynamics near the fold in the parameter regime corresponding to the presence of small oscillations and the other on the global dynamics, more specifically on the existence of an attracting limit cycle. Our local result is a rigorous characterization of small canards and sectors of rotation in the case of folded node with an additional time scale, a feature allowing for a clear geometric argument. The global result gives the existence of an attracting unique limit cycle, which, in some parameter regimes, remains attracting and unique even during passages through a canard explosion.

No abstract