1
|
Graham EJ, Elhadad N, Albers D. Reduced model for female endocrine dynamics: Validation and functional variations. Math Biosci 2023; 358:108979. [PMID: 36792027 DOI: 10.1016/j.mbs.2023.108979] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2022] [Revised: 01/19/2023] [Accepted: 02/07/2023] [Indexed: 02/15/2023]
Abstract
A normally functioning menstrual cycle requires significant crosstalk between hormones originating in ovarian and brain tissues. Reproductive hormone dysregulation may cause abnormal function and sometimes infertility. The inherent complexity in this endocrine system is a challenge to identifying mechanisms of cycle disruption, particularly given the large number of unknown parameters in existing mathematical models. We develop a new endocrine model to limit model complexity and use simulated distributions of unknown parameters for model analysis. By employing a comprehensive model evaluation, we identify a collection of mechanisms that differentiate normal and abnormal phenotypes. We also discover an intermediate phenotype-displaying relatively normal hormone levels and cycle dynamics-that is grouped statistically with the irregular phenotype. Results provide insight into how clinical symptoms associated with ovulatory disruption may not be detected through hormone measurements alone.
Collapse
Affiliation(s)
- Erica J Graham
- Mathematics Department, Bryn Mawr College, Bryn Mawr, PA 19010, USA.
| | - Noémie Elhadad
- Department of Biomedical Informatics, Columbia University, New York, NY 10032, USA
| | - David Albers
- Pediatrics Department, University of Colorado Denver-Anschutz Medical Campus, Aurora, CO 80045, USA
| |
Collapse
|
2
|
Grindstaff JL, Beaty LE, Ambardar M, Luttbeg B. Integrating theoretical and empirical approaches for a robust understanding of endocrine flexibility. J Exp Biol 2022; 225:274311. [PMID: 35258612 PMCID: PMC8987727 DOI: 10.1242/jeb.243408] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
There is growing interest in studying hormones beyond single 'snapshot' measurements, as recognition that individual variation in the endocrine response to environmental change may underlie many rapid, coordinated phenotypic changes. Repeated measures of hormone levels in individuals provide additional insight into individual variation in endocrine flexibility - that is, how individuals modulate hormone levels in response to the environment. The ability to quickly and appropriately modify phenotype is predicted to be favored by selection, especially in unpredictable environments. The need for repeated samples from individuals can make empirical studies of endocrine flexibility logistically challenging, but methods based in mathematical modeling can provide insights that circumvent these challenges. Our Review introduces and defines endocrine flexibility, reviews existing studies, makes suggestions for future empirical work, and recommends mathematical modeling approaches to complement empirical work and significantly advance our understanding. Mathematical modeling is not yet widely employed in endocrinology, but can be used to identify innovative areas for future research and generate novel predictions for empirical testing.
Collapse
Affiliation(s)
| | - Lynne E Beaty
- School of Science, Penn State Erie - The Behrend College, Erie, PA 16563, USA
| | - Medhavi Ambardar
- Department of Biological Sciences, Fort Hays State University, Hays, KS 67601, USA
| | - Barney Luttbeg
- Department of Integrative Biology, Oklahoma State University, OK 74078, USA
| |
Collapse
|
3
|
Clément F, Crépieux P, Yvinec R, Monniaux D. Mathematical modeling approaches of cellular endocrinology within the hypothalamo-pituitary-gonadal axis. Mol Cell Endocrinol 2020; 518:110877. [PMID: 32569857 DOI: 10.1016/j.mce.2020.110877] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/27/2020] [Revised: 05/13/2020] [Accepted: 05/19/2020] [Indexed: 01/26/2023]
Abstract
The reproductive neuroendocrine axis, or hypothalamo-pituitary-gonadal (HPG) axis, is a paragon of complex biological system involving numerous cell types, spread over several anatomical levels communicating through entangled endocrine feedback loops. The HPG axis exhibits remarkable dynamic behaviors on multiple time and space scales, which are an inexhaustible source of studies for mathematical and computational biology. In this review, we will describe a variety of modeling approaches of the HPG axis from a cellular endocrinology viewpoint. We will in particular investigate the questions raised by some of the most striking features of the HPG axis: (i) the pulsatile secretion of hypothalamic and pituitary hormones, and its counterpart, the cell signaling induced by frequency-encoded hormonal signals, and (ii) the dual, gametogenic and glandular function of the gonads, which relies on the tight control of the somatic cell populations ensuring the proper maturation and timely release of the germ cells.
Collapse
Affiliation(s)
- Frédérique Clément
- Inria, Centre de Recherche Inria Saclay-Île-de-France, Palaiseau, France.
| | - Pascale Crépieux
- INRAE, UMR85, Unité Physiologie de la Reproduction et des Comportements, F-37380, Nouzilly, France; CNRS, UMR7247, F-37380, Nouzilly, France; Université de Tours, F-37041, Tours, France
| | - Romain Yvinec
- INRAE, UMR85, Unité Physiologie de la Reproduction et des Comportements, F-37380, Nouzilly, France; CNRS, UMR7247, F-37380, Nouzilly, France; Université de Tours, F-37041, Tours, France
| | - Danielle Monniaux
- INRAE, UMR85, Unité Physiologie de la Reproduction et des Comportements, F-37380, Nouzilly, France; CNRS, UMR7247, F-37380, Nouzilly, France; Université de Tours, F-37041, Tours, France
| |
Collapse
|
4
|
Churilov AN, Milton J, Salakhova ER. An integrate-and-fire model for pulsatility in the neuroendocrine system. CHAOS (WOODBURY, N.Y.) 2020; 30:083132. [PMID: 32872840 DOI: 10.1063/5.0010553] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/13/2020] [Accepted: 07/08/2020] [Indexed: 06/11/2023]
Abstract
A model for pulsatility in neuroendocrine regulation is proposed which combines Goodwin-type feedback control with impulsive input from neurons located in the hypothalamus. The impulsive neural input is modeled using an integrate-and-fire mechanism; namely, inputs are generated only when the membrane potential crosses a threshold, after which it is reset to baseline. The resultant model takes the form of a functional-differential equation with continuous and impulsive components. Despite the impulsive nature of the inputs, realistic hormone profiles are generated, including ultradian and circadian rhythms, pulsatile secretory patterns, and even chaotic dynamics.
Collapse
Affiliation(s)
- Alexander N Churilov
- Faculty of Mathematics and Mechanics, St. Petersburg State University, Universitetsky av. 28, Stary Peterhof, 198504 St. Petersburg, Russia
| | - John Milton
- Keck Science Department, The Claremont Colleges, 925 North Mills Ave., Claremont, California 91711, USA
| | - Elvira R Salakhova
- Faculty of Mathematics and Mechanics, St. Petersburg State University, Universitetsky av. 28, Stary Peterhof, 198504 St. Petersburg, Russia
| |
Collapse
|
5
|
Jessus C, Munro C, Houliston E. Managing the Oocyte Meiotic Arrest-Lessons from Frogs and Jellyfish. Cells 2020; 9:E1150. [PMID: 32392797 PMCID: PMC7290932 DOI: 10.3390/cells9051150] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2020] [Revised: 05/03/2020] [Accepted: 05/05/2020] [Indexed: 12/11/2022] Open
Abstract
During oocyte development, meiosis arrests in prophase of the first division for a remarkably prolonged period firstly during oocyte growth, and then when awaiting the appropriate hormonal signals for egg release. This prophase arrest is finally unlocked when locally produced maturation initiation hormones (MIHs) trigger entry into M-phase. Here, we assess the current knowledge of the successive cellular and molecular mechanisms responsible for keeping meiotic progression on hold. We focus on two model organisms, the amphibian Xenopus laevis, and the hydrozoan jellyfish Clytia hemisphaerica. Conserved mechanisms govern the initial meiotic programme of the oocyte prior to oocyte growth and also, much later, the onset of mitotic divisions, via activation of two key kinase systems: Cdk1-Cyclin B/Gwl (MPF) for M-phase activation and Mos-MAPkinase to orchestrate polar body formation and cytostatic (CSF) arrest. In contrast, maintenance of the prophase state of the fully-grown oocyte is assured by highly specific mechanisms, reflecting enormous variation between species in MIHs, MIH receptors and their immediate downstream signalling response. Convergence of multiple signalling pathway components to promote MPF activation in some oocytes, including Xenopus, is likely a heritage of the complex evolutionary history of spawning regulation, but also helps ensure a robust and reliable mechanism for gamete production.
Collapse
Affiliation(s)
- Catherine Jessus
- Laboratoire de Biologie du Développement - Institut de Biologie Paris Seine, LBD - IBPS, Sorbonne Université, CNRS, F-75005 Paris, France
| | - Catriona Munro
- Laboratoire de Biologie du Développement de Villefranche-sur-mer (LBDV), Sorbonne Université, CNRS, 06230 Villefranche-sur-mer, France;
- Inserm, Center for Interdisciplinary Research in Biology, Collège de France, PSL Research University, CNRS, 75005 Paris, France
| | - Evelyn Houliston
- Laboratoire de Biologie du Développement de Villefranche-sur-mer (LBDV), Sorbonne Université, CNRS, 06230 Villefranche-sur-mer, France;
| |
Collapse
|
6
|
D'Occhio MJ, Baruselli PS, Campanile G. Influence of nutrition, body condition, and metabolic status on reproduction in female beef cattle: A review. Theriogenology 2019; 125:277-284. [DOI: 10.1016/j.theriogenology.2018.11.010] [Citation(s) in RCA: 50] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2018] [Revised: 11/06/2018] [Accepted: 11/14/2018] [Indexed: 12/19/2022]
|
7
|
Zhang Q, Li J, Middleton A, Bhattacharya S, Conolly RB. Bridging the Data Gap From in vitro Toxicity Testing to Chemical Safety Assessment Through Computational Modeling. Front Public Health 2018; 6:261. [PMID: 30255008 PMCID: PMC6141783 DOI: 10.3389/fpubh.2018.00261] [Citation(s) in RCA: 40] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2018] [Accepted: 08/21/2018] [Indexed: 12/18/2022] Open
Abstract
Chemical toxicity testing is moving steadily toward a human cell and organoid-based in vitro approach for reasons including scientific relevancy, efficiency, cost, and ethical rightfulness. Inferring human health risk from chemical exposure based on in vitro testing data is a challenging task, facing various data gaps along the way. This review identifies these gaps and makes a case for the in silico approach of computational dose-response and extrapolation modeling to address many of the challenges. Mathematical models that can mechanistically describe chemical toxicokinetics (TK) and toxicodynamics (TD), for both in vitro and in vivo conditions, are the founding pieces in this regard. Identifying toxicity pathways and in vitro point of departure (PoD) associated with adverse health outcomes requires an understanding of the molecular key events in the interacting transcriptome, proteome, and metabolome. Such an understanding will in turn help determine the sets of sensitive biomarkers to be measured in vitro and the scope of toxicity pathways to be modeled in silico. In vitro data reporting both pathway perturbation and chemical biokinetics in the culture medium serve to calibrate the toxicity pathway and virtual tissue models, which can then help predict PoDs in response to chemical dosimetry experienced by cells in vivo. Two types of in vitro to in vivo extrapolation (IVIVE) are needed. (1) For toxic effects involving systemic regulations, such as endocrine disruption, organism-level adverse outcome pathway (AOP) models are needed to extrapolate in vitro toxicity pathway perturbation to in vivo PoD. (2) Physiologically-based toxicokinetic (PBTK) modeling is needed to extrapolate in vitro PoD dose metrics into external doses for expected exposure scenarios. Linked PBTK and TD models can explore the parameter space to recapitulate human population variability in response to chemical insults. While challenges remain for applying these modeling tools to support in vitro toxicity testing, they open the door toward population-stratified and personalized risk assessment.
Collapse
Affiliation(s)
- Qiang Zhang
- Department of Environmental Health, Rollins School of Public Health, Emory University, Atlanta, GA, United States
| | - Jin Li
- Unilever, Safety and Environmental Assurance Centre, Colworth Science Park, Sharnbrook, United Kingdom
| | - Alistair Middleton
- Unilever, Safety and Environmental Assurance Centre, Colworth Science Park, Sharnbrook, United Kingdom
| | - Sudin Bhattacharya
- Biomedical Engineering, Michigan State University, East Lansing, MI, United States
| | - Rory B Conolly
- Integrated Systems Toxicology Division, National Health and Environmental Effects Research Laboratory, United States Environmental Protection Agency, Durham, NC, United States
| |
Collapse
|
8
|
Coupling a reproductive function model to a productive function model to simulate lifetime performance in dairy cows. Animal 2018; 13:570-579. [PMID: 30037359 DOI: 10.1017/s1751731118001830] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
Reproductive success is a key component of lifetime performance in dairy cows but is difficult to predict due to interactions with productive function. Accordingly, this study introduces a dynamic model to simulate the productive and reproductive performance of a cow during her lifetime. The cow model consists of an existing productive function model (GARUNS) which is coupled to a new reproductive function model (RFM). The GARUNS model simulates the individual productive performance of a dairy cow throughout her lifespan. It provides, with a daily time step, changes in BW and composition, fetal growth, milk yield and composition and food intake. Genetic-scaling parameters are incorporated to scale individual performance and simulate differences within and between breeds. GARUNS responds to the discrete event signals 'conception' and 'death' (of embryo or fetus) generated by RFM. In turn, RFM responds to the GARUNS outputs concerning the cow's energetic status: the daily total processed metabolizable energy per kg BW (TPEW) and the net energy balance (EB). Reproductive function model models the reproductive system as a compartmental system transitioning between nine competence stages: prepubertal (PRPB), anestrous (ANST), anovulatory (ANOV), pre-ovulating (PREO), ovulating (OVUL), post-ovulating (PSTO), luteinizing (LUTZ), luteal (LUTL) and gestating (GEST). The transition from PRPB to ANST represents the start of reproductive activity at puberty. The cyclic path through ANST, PREO, OVUL, PSTO, LUTZ and LUTL forms the regime of ovulatory cycles, whereas ANOV and GEST are transient stages that interrupt this regime. Anovulatory refers explicitly to a stage in which ovulation cannot occur (i.e. interrupted cyclicity), whereas ANST is a pivotal stage within ovulatory cycles. Reproductive function model generates estradiol and progesterone hormonal profiles consistent with reference profiles derived from literature. Cyclicity is impacted by the GARUNS output EB and clearance of estradiol is impacted by TPEW. A farming system model was designed to describe different farm protocols of heat detection, insemination, feeding (amount and energy density), drying-off and culling. Results of model simulation (10 000 simulations of individual cows over 5000 days lifetime period, with randomly drawn genetic-scaling parameters and standard diet) are consistent with literature for reproductive performance. This model allows simulation of deviations in reproductive trajectories along physiological stages of the cow reproductive cycle. It thus provides the basis for evaluation of the relative importance of different factors affecting fertility at individual cow and herd levels across different breeds and management environments.
Collapse
|