Margaritis Voliotis

Stress can disrupt menstrual cycles, impair fertility and cause reproductive disfunction. The posterodorsal medial amygdala (MePD) integrates stress signals and regulates the gonadotropin-releasing hormone (GnRH) pulse generator through a dense network of GABA and Urocortin-3 (UCN3) neurons, yet the mechanisms underlying the circuitry remain poorly understood. Here, we combine in vivo mini-endoscopic calcium imaging, optogenetics, clustering analysis, and computational modeling to investigate the MePD circuitry in female mice. We uncover two anti-correlated GABA subpopulations in the MePD that are involved in the response to restraint stress and UCN3 neuron stimulation. Computational modeling suggests that mutual inhibition between these GABA groups drives their anti-correlated activity and predicts how these interactions shape downstream responses to stimulation of GABA and UCN3 neurons. In vivo optogenetics confirms that GABA neurons are critical for transmitting UCN3 signals to regulate luteinizing hormone (LH) pulse frequency. Together, our findings reveal amygdala GABAergic circuit mechanisms that mediate stress effects on reproductive health, linking emotional processing and neuroendocrine control.

The gonadotrophin-releasing hormone (GnRH) pulse generator is a critical neural oscillator that governs reproductive function through the pulsatile release of luteinising hormone (LH) and follicle stimulating hormone (FSH). Early electrophysiological studies, notably by Ernst Knobil, identified multiunit activity (MUA) volleys in the mediobasal hypothalamus that aligned with LH pulses, suggesting a neural basis for GnRH pulsatility. Although GnRH neurons exhibit some intrinsic secretory rhythmicity in vitro , their isolated electrophysiological signatures have proven inconsistent. Recent advances, including GCaMP fibre photometry in freely behaving mice, have revealed a precise correlation between episodic GnRH distal processes and LH pulses. However, it is now well established that a neural oscillator comprising hypothalamic kisspeptin neurones co-expressing neurokinin B and dynorphin, collectively termed the KNDy network, represents the core construct of the GnRH pulse generator. Understanding the dynamics of this network and its modulation by external inputs such as stress, metabolic cues, and circadian rhythms is essential. Computational modelling provides a systematic framework for integrating experimental data with mechanistic and predictive analyses to decode the GnRH pulse generator dynamics.

Stress can disrupt menstrual cycles, cause infertility, and lead to other reproductive disorders. The posterodorsal medial amygdala (MePD) processes stress signals and regulates the gonadotropin-releasing hormone (GnRH) pulse generator through GABAergic inhibitory projections to the hypothalamus. However, how stress is processed in the MePD - especially involving its dense GABA and Urocortin-3 (UCN3) neurons - remains poorly understood. In this study, we combine in vivo GRadient-INdex (GRIN) lens mini-endoscopic calcium imaging (to track neuronal activity), optogenetics, clustering analysis, and computational modeling to investigate the MePD circuitry. Our findings reveal two anti-correlated GABA subpopulations in the MePD that dictate responses to both UCN3 neuron stimulation and restraint stress. Our computational modeling suggests that mutual inhibition between these GABA groups drives the anti-correlated activity and predicts how these interactions shape downstream responses to stimulation of GABA and UCN3 neurons. We test these predictions using optogenetics and confirm that GABA neurons are critical for the transmission of UCN3 signals to regulate luteinizing hormone (LH) pulse frequency. Our study is the first to show how GABA neurons in the amygdala mediate stress effects on reproductive health, uncovering key neural mechanisms linking emotional and reproductive functions.

Technical limitations have prevented understanding of how growth factor signals are encoded in distinct activity patterns of the phosphoinositide 3-kinase (PI3K)/AKT pathway, and how this is altered by oncogenic pathway mutations. We introduce a kinetic, single-cell framework for precise calculations of PI3K-specific information transfer for different growth factors. This features live-cell imaging of PI3K/AKT activity reporters and multiplexed CyTOF measurements of PI3K/AKT and RAS/ERK signaling markers over time. Using this framework, we found that the PIK3CA ^H1047R oncogene was not a simple, constitutive activator of the pathway as often presented. Dose-dependent expression of PIK3CA ^H1047R in human cervical cancer and induced pluripotent stem cells corrupted the fidelity of growth factor-induced information transfer, with preferential amplification of epidermal growth factor receptor (EGFR) signaling responses compared to insulin-like growth factor 1 (IGF1) and insulin receptor signaling. PIK3CA ^H1047R did not only shift these responses to a higher mean but also enhanced signaling heterogeneity. We conclude that oncogenic PIK3CA ^H1047R corrupts information transfer in a growth factor-dependent manner and suggest new opportunities for tuning of receptor-specific PI3K pathway outputs for therapeutic benefit.

Hypothalamic kisspeptin (Kiss1) neurons are vital for pubertal development and reproduction. Arcuate nucleus Kiss1 (Kiss1^ARH) neurons are responsible for the pulsatile release of gonadotropin-releasing hormone (GnRH). In females, the behavior of Kiss1^ARH neurons, expressing Kiss1, neurokinin B (NKB), and dynorphin (Dyn), varies throughout the ovarian cycle. Studies indicate that 17β-estradiol (E2) reduces peptide expression but increases Slc17a6 (Vglut2) mRNA and glutamate neurotransmission in these neurons, suggesting a shift from peptidergic to glutamatergic signaling. To investigate this shift, we combined transcriptomics, electrophysiology, and mathematical modeling. Our results demonstrate that E2 treatment upregulates the mRNA expression of voltage-activated calcium channels, elevating the whole-cell calcium current that contributes to high-frequency burst firing. Additionally, E2 treatment decreased the mRNA levels of canonical transient receptor potential (TPRC) 5 and G protein-coupled K⁺ (GIRK) channels. When Trpc5 channels in Kiss1^ARH neurons were deleted using CRISPR/SaCas9, the slow excitatory postsynaptic potential was eliminated. Our data enabled us to formulate a biophysically realistic mathematical model of Kiss1^ARH neurons, suggesting that E2 modifies ionic conductances in these neurons, enabling the transition from high-frequency synchronous firing through NKB-driven activation of TRPC5 channels to a short bursting mode facilitating glutamate release. In a low E2 milieu, synchronous firing of Kiss1^ARH neurons drives pulsatile release of GnRH, while the transition to burst firing with high, preovulatory levels of E2 would facilitate the GnRH surge through its glutamatergic synaptic connection to preoptic Kiss1 neurons.

Hypothalamic kisspeptin (Kiss1) neurons are vital for pubertal development and reproduction. Arcuate nucleus Kiss1 (Kiss1ARH) neurons are responsible for the pulsatile release of Gonadotropin-releasing Hormone (GnRH). In females, the behavior of Kiss1ARH neurons, expressing Kiss1, Neurokinin B (NKB), and Dynorphin (Dyn), varies throughout the ovarian cycle. Studies indicate that 17β-estradiol (E2) reduces peptide expression but increases Vglut2 mRNA and glutamate neurotransmission in these neurons, suggesting a shift from peptidergic to glutamatergic signaling. To investigate this shift, we combined transcriptomics, electrophysiology, and mathematical modeling. Our results demonstrate that E2 treatment upregulates the mRNA expression of voltage-activated calcium channels, elevating the whole-cell calcium current and that contribute to high-frequency burst firing. Additionally, E2 treatment decreased the mRNA levels of Canonical Transient Receptor Potential (TPRC) 5 and G protein-coupled K+ (GIRK) channels. When TRPC5 channels in Kiss1ARH neurons were deleted using CRISPR, the slow excitatory postsynaptic potential (sEPSP) was eliminated. Our data enabled us to formulate a biophysically realistic mathematical model of the Kiss1ARH neuron, suggesting that E2 modifies ionic conductances in Kiss1ARH neurons, enabling the transition from high frequency synchronous firing through NKB-driven activation of TRPC5 channels to a short bursting mode facilitating glutamate release. In a low E2 milieu, synchronous firing of Kiss1ARH neurons drives pulsatile release of GnRH, while the transition to burst firing with high, preovulatory levels of E2 would facilitate the GnRH surge through its glutamatergic synaptic connection to preoptic Kiss1 neurons.

Normal reproductive function and fertility rely on the rhythmic secretion of gonadotropin-releasing hormone (GnRH), which is driven by the hypothalamic GnRH pulse generator. A key regulator of the GnRH pulse generator is the posterodorsal subnucleus of the medial amygdala (MePD), a brain region that is involved in processing external environmental cues, including the effect of stress. However, the neuronal pathways enabling the dynamic, stress-triggered modulation of GnRH secretion remain largely unknown. Here, we employ in silico modelling in order to explore the impact of dynamic inputs on GnRH pulse generator activity. We introduce and analyse a mathematical model representing MePD neuronal circuits composed of GABAergic and glutamatergic neuronal populations, integrating it with our GnRH pulse generator model. Our analysis dissects the influence of excitatory and inhibitory MePD projections’ outputs on the GnRH pulse generator’s activity and reveals a functionally relevant MePD glutamatergic projection to the GnRH pulse generator, which we probe with in vivo optogenetics. Our study sheds light on how MePD neuronal dynamics affect the GnRH pulse generator activity and offers insights into stress-related dysregulation.

Psychosocial stress negatively impacts reproductive function by inhibiting pulsatile luteinizing hormone (LH) secretion. The posterodorsal medial amygdala (MePD) is responsible in part for processing stress and modulating the reproductive axis. Activation of the neurokinin 3 receptor (NK3R) suppresses the gonadotropin‐releasing hormone (GnRH) pulse generator, under hypoestrogenic conditions, and NK3R activity in the amygdala has been documented to play a role in stress and anxiety. We investigate whether NK3R activation in the MePD is involved in mediating the inhibitory effect of psychosocial stress on LH pulsatility in ovariectomised female mice. First, we administered senktide, an NK3R agonist, into the MePD and monitored the effect on pulsatile LH secretion. We then delivered SB222200, a selective NK3R antagonist, intra‐MePD in the presence of predator odour, 2,4,5‐trimethylthiazole (TMT) and examined the effect on LH pulses. Senktide administration into the MePD dose‐dependently suppresses pulsatile LH secretion. Moreover, NK3R signalling in the MePD mediates TMT‐induced suppression of the GnRH pulse generator, which we verified using a mathematical model. The model verifies our experimental findings: (i) predator odour exposure inhibits LH pulses, (ii) activation of NK3R in the MePD inhibits LH pulses and (iii) NK3R antagonism in the MePD blocks stressor‐induced inhibition of LH pulse frequency in the absence of ovarian steroids. These results demonstrate for the first time that NK3R neurons in the MePD mediate psychosocial stress‐induced suppression of the GnRH pulse generator.

The hypothalamus is the central regulator of reproductive hormone secretion. Pulsatile secretion of gonadotropin releasing hormone (GnRH) is fundamental to physiological stimulation of the pituitary gland to release luteinizing hormone (LH) and follicle stimulating hormone (FSH). Furthermore, GnRH pulsatility is altered in common reproductive disorders such as polycystic ovary syndrome (PCOS) and hypothalamic amenorrhea (HA). LH is measured routinely in clinical practice using an automated chemiluminescent immunoassay method and is the gold standard surrogate marker of GnRH. LH can be measured at frequent intervals (e.g., 10 minutely) to assess GnRH/LH pulsatility. However, this is rarely done in clinical practice because it is resource intensive, and there is no open-access, graphical interface software for computational analysis of the LH data available to clinicians. Here we present hormoneBayes, a novel open-access Bayesian framework that can be easily applied to reliably analyze serial LH measurements to assess LH pulsatility. The framework utilizes parsimonious models to simulate hypothalamic signals that drive LH dynamics, together with state-of-the-art (sequential) Monte-Carlo methods to infer key parameters and latent hypothalamic dynamics. We show that this method provides estimates for key pulse parameters including inter-pulse interval, secretion and clearance rates and identifies LH pulses in line with the widely used deconvolution method. We show that these parameters can distinguish LH pulsatility in different clinical contexts including in reproductive health and disease in men and women (e.g., healthy men, healthy women before and after menopause, women with HA or PCOS). A further advantage of hormoneBayes is that our mathematical approach provides a quantified estimation of uncertainty. Our framework will complement methods enabling real-time in-vivo hormone monitoring and therefore has the potential to assist translation of personalized, data-driven, clinical care of patients presenting with conditions of reproductive hormone dysfunction.

Hypothalamic kisspeptin (Kiss1) neurons are vital for pubertal development and reproduction. Arcuate nucleus Kiss1 (Kiss1 ^ARH ) neurons are responsible for the pulsatile release of Gonadotropin-releasing Hormone (GnRH). In females, the behavior of Kiss1 ^ARH neurons, expressing Kiss1, Neurokinin B (NKB), and Dynorphin (Dyn), varies throughout the ovarian cycle. Studies indicate that 17β-estradiol (E2) reduces peptide expression but increases Vglut2 mRNA and glutamate neurotransmission in these neurons, suggesting a shift from peptidergic to glutamatergic signaling. To investigate this shift, we combined transcriptomics, electrophysiology, and mathematical modeling. Our results demonstrate that E2 treatment upregulates the mRNA expression of voltage-activated calcium channels, elevating the whole-cell calcium current and that contribute to high-frequency burst firing. Additionally, E2 treatment decreased the mRNA levels of Canonical Transient Receptor Potential (TPRC) 5 and G protein-coupled K ⁺ (GIRK) channels. When TRPC5 channels in Kiss1 ^ARH neurons were deleted using CRISPR, the slow excitatory postsynaptic potential (sEPSP) was eliminated. Our data enabled us to formulate a biophysically realistic mathematical model of the Kiss1 ^ARH neuron, suggesting that E2 modifies ionic conductances in Kiss1 ^ARH neurons, enabling the transition from high frequency synchronous firing through NKB-driven activation of TRPC5 channels to a short bursting mode facilitating glutamate release. In a low E2 milieu, synchronous firing of Kiss1 ^ARH neurons drives pulsatile release of GnRH, while the transition to burst firing with high, preovulatory levels of E2 would facilitate the GnRH surge through its glutamatergic synaptic connection to preoptic Kiss1 neurons.