Modulation of pulsatile GnRH dynamics across the ovarian cycle via changes in the network excitability and basal activity of the arcuate kisspeptin network

Abstract
Editor's evaluation
Introduction
Results
Discussion
Materials and methods
Data availability
References
Article and author information
Metrics

Abstract

Pulsatile GnRH release is essential for normal reproductive function. Kisspeptin secreting neurons found in the arcuate nucleus, known as KNDy neurons for co-expressing neurokinin B, and dynorphin, drive pulsatile GnRH release. Furthermore, gonadal steroids regulate GnRH pulsatile dynamics across the ovarian cycle by altering KNDy neurons' signalling properties. However, the precise mechanism of regulation remains mostly unknown. To better understand these mechanisms, we start by perturbing the KNDy system at different stages of the estrous cycle using optogenetics. We find that optogenetic stimulation of KNDy neurons stimulates pulsatile GnRH/LH secretion in estrous mice but inhibits it in diestrous mice. These in vivo results in combination with mathematical modelling suggest that the transition between estrus and diestrus is underpinned by well-orchestrated changes in neuropeptide signalling and in the excitability of the KNDy population controlled via glutamate signalling. Guided by model predictions, we show that blocking glutamate signalling in diestrous animals inhibits LH pulses, and that optic stimulation of the KNDy population mitigates this inhibition. In estrous mice, disruption of glutamate signalling inhibits pulses generated via sustained low-frequency optic stimulation of the KNDy population, supporting the idea that the level of network excitability is critical for pulse generation. Our results reconcile previous puzzling findings regarding the estradiol-dependent effect that several neuromodulators have on the GnRH pulse generator dynamics. Therefore, we anticipate our model to be a cornerstone for a more quantitative understanding of the pathways via which gonadal steroids regulate GnRH pulse generator dynamics. Finally, our results could inform useful repurposing of drugs targeting the glutamate system in reproductive therapy.

Editor's evaluation

Successful reproduction requires luteinizing hormone that is secreted from the pituitary gland in pulses which vary over the female reproductive cycle. The pulses arise by patterned secretion of a hypothalamic 'releasing factor', the secretion of which is itself governed by a population of hypothalamic neurons that express the neuropeptide kisspeptin. The paper by Voliotis et al. combines novel experimental evidence from transgenic mice with an elegant mathematical model to analyze how the kisspeptin neurons generate the varying pulsatile patterns. Testing this model in a way that will either confirm or refute it, may prove challenging.

https://doi.org/10.7554/eLife.71252.sa0

Introduction

The dynamics of gonadotropin-releasing hormone (GnRH) secretion is critical for reproductive health. In female animals, GnRH secretion is tightly regulated across the ovarian cycle. Pulsatile secretion dominates most of the cycle, with frequency and amplitude modulated by the ovarian steroid feedback. Positive feedback from increasing estradiol levels triggers a preovulatory surge of GnRH/LH secretion (Christian and Moenter, 2010). Furthermore, there is ample evidence that ARC kisspeptin neurons are prime mediators of the ovarian steroid feedback on the pulsatile dynamics of GnRH/LH secretion (McQuillan et al., 2019), although the mechanisms remain unclear.

In vitro studies have shown that gonadal steroids have a dramatic effect on the electrophysiology of ARC kisspeptin neurons. For instance, spontaneous firing activity of ARC kisspeptin neurons from castrated mice appears elevated compared to intact animals (Ruka et al., 2016) and estradiol replacement attenuates ARC kisspeptin neuron activity in gonadectomised animals (Ruka et al., 2016; Wang et al., 2018). More recently, fibre photometry data from female mice show that the ARC kisspeptin neuronal population (KNDy network) pulses at a relatively constant frequency throughout the ovarian cycle apart from the estrous phase where the frequency is dramatically reduced (McQuillan et al., 2019). This slowdown of LH frequency is thought to be a direct consequence of the increasing progesterone levels associated with ovulation (McQuillan et al., 2019), although studies using the rhesus monkey show that raising pre-ovulatory estrogen levels are also important (O’Byrne et al., 1991). Studies in sheep indicate the inhibitory effects of progesterone are mediated through increased dynorphin signalling (Goodman et al., 2011; Moore et al., 2018); however, this is less clear in mice where ovarian steroids have a negative effect on Dyn mRNA levels (Navarro et al., 2009).

Perplexing is also the differential effect of various neuromodulators on LH secretion depending on the gonadal steroid background. For instance, N-methyl-D-aspartate (NMDA) robustly inhibits LH pulses in the ovariectomised monkey, whereas in the presence of estradiol this effect is reversed, and NMDA stimulates LH secretion (Reyes et al., 1990; Reyes et al., 1991). Similar reversal of action on LH dynamics depending on the underlying ovarian steroid milieu has been also documented for other neurotransmitter and neuropeptides in other species (Kalra and Kalra, 1983; Brann and Mahesh, 1992; Arias et al., 1993; Bonavera et al., 1994; Scorticati et al., 2004) and highlights the complex mechanisms underlying the modulation the GnRH pulse generator by gonadal steroids.

Here, using mathematical modelling along with optogenetic stimulation of ARC kisspeptin neurons, we embark to understand how the dynamics of the pulse generator are modulated across the ovarian cycle. Our mathematical model suggests that the level of excitability within the ARC kisspeptin network—the propensity of kisspeptin neurons to signal and activate each other—is one of the key parameters modulated in different stages of the cycle by gonadal steroids. Previous studies have shown that ARC kisspeptin neurons synapse on each other (Yip et al., 2015; Qiu et al., 2016) and are glutamatergic (Cravo et al., 2011; Qiu et al., 2011; Kelly et al., 2013; Nestor et al., 2016; Qiu et al., 2016; Wang et al., 2018). Based on these findings, we hypothesise that population excitability should be enabled primarily via glutamate signalling. We test our predictions in vivo and show that glutamatergic transmission is an important factor for the pulsatile behaviour of the KNDy network.

Results

The dynamic response of the KNDy network to sustained, low-frequency optic stimulation is estrous cycle dependent

Using optogenetics we perturbed the KNDy network to test whether and how sex steroids modulate the system’s dynamical response. ARC kisspeptin-expressing neurons were transduced with a Cre-dependent adeno-associated virus (AAV9-EF1-dflox-hChR2-(H134R)-mCherryWPRE-hGH) to express ChR2 (Figure 1; see Materials and methods) and were optogenetically stimulated at the estrous and the diestrous phase of the cycle, measuring LH pulse frequency as a readout. Sustained, low-frequency optic stimulation was used to emulate elevated basal activity in ARC kisspeptin neurons or persistent stimulatory signals to the KNDy population from other neuronal populations.

Figure 1

Download asset Open asset

Expression of arcuate nucleus (ARC) kisspeptin neurones with ChR2-mCherry in Kiss-Cre mouse.

Coronal section showing red mCherry fluorescence positive neurons in the ARC which indicates ChR2 receptor expressing kisspeptin neurones, following unilateral injection of AAV9.EF1.dflox.hChR2(H134R)-mCherry.WPRE.hGH into the ARC of Kiss-Cre mouse. Note the absence of mCherry fluorescence in the other side of ARC.

3 V, Third ventricle.

In estrous mice, we find that sustained optogenetic stimulation of ARC kisspeptin neurons at 5 Hz immediately triggers robust LH pulses at a frequency of 2.10 ± 0.24 pulses/hour (Figure 2A and C & E), which is in agreement with our previous findings (Voliotis et al., 2019) and highlights how pulsatile dynamics can emerge as a population phenomenon without the need of a pulsatile activation signal (Strogatz, 2018). In diestrous mice, on the other hand, optogenetic stimulation of ARC kisspeptin neurons at 5 Hz has a subtle slowdown effect on LH pulse frequency over the 1.5 hr stimulation period (Figure 2—figure supplement 1). To investigate this response in greater detail, we revised our experimental protocol, removing the control period and extending the stimulation period to 2.5 hr. With the extended protocol we measure 0.64 ± 0.09 and 0.40 ± 0.13 LH pulses/hr under sustained optic stimulation at 5 and 15 Hz, respectively; these frequencies are significantly lower than the LH pulse frequency we observe in control animals, which receive no optic stimulation (Figure 2B, D, F). We note that we observe normal LH pulse frequencies in WT animals receiving sustained optic stimulation for 2.5 hr (Figure 2—figure supplement 2).

Figure 2 with 2 supplements see all

Download asset Open asset

Differential effect of optic stimulation of ARC kisspeptin neurons in estrous and diestrous Kiss-Cre mice.

(**A–B**) Representative examples showing LH secretion in response to no stimulation (grey bar) or sustained blue light (473 nm, 5 ms pulse width, black bar) activation of kisspeptin neurons at 5 Hz in estrous (C) and diestrous (D) mice. (E) Summary showing mean ± SEM LH pulse frequency over the 60 min control period (white bars) and over the subsequent stimulation period (black bar) in estrous mice. (F) Summary showing mean ± SEM LH pulse frequency in the control (grey bar) and stimulated (black bars) diestrous mice. *Denote LH pulses. ^#p < 0.05 vs control; ^†p < 0.05 vs pre-stimulation; n = 5–6 per group.

Our data illustrate how natural variation of ovarian steroids across the ovarian cycle leads to qualitative changes in the dynamical response of ARC kisspeptin neurons to optical stimulation. These changes are most probably driven by the effect that gonadal steroids have on the intrinsic electrophysiological properties of ARC kisspeptin neurons (Ruka et al., 2016) and the neuromodulator signalling capacity within the KNDy network (Vanacker et al., 2017).

A mathematical model predicts key mechanisms modulating the behaviour of the KNDy pulse generator across the estrous cycle

Interrogating the KNDy network at different stages of the estous cycle via optic stimulation and measuring the effect on LH pulse frequency allows use of our mathematical model (Voliotis et al., 2019) to understand how key system parameters change under gonadal steroid control. The model describes the dynamical behaviour of ARC kisspeptin neurons using three dynamical variables: representing the levels of Dyn, NKB and neuronal activity (Figure 3A). Furthermore, rather than focusing on the biophyscal details of regulation, the model postulates that gonadal steroids could potentially modulate the behaviour of the KNDy system across the cycle via acting on four system-level parameters: (i) level of Dyn signalling, (ii) level of NKB signalling, (iii) network excitability (i.e. propensity of neurons in the population to transmit signals to one another), and (iv) basal neuronal activity.

Figure 3 with 1 supplement see all

Download asset Open asset

Model predictions on the key mechanisms modulating the behaviour of the KNDy pulse generator across the ovarian cycle.

(A) Schematic illustration of the coarse-grained model of the ARC KNDy population. The model comprises three dynamical variables representing the average levels of Dyn and NKB secreted by the population, and its average firing activity. We hypothesise that four key parameters modulate the behaviour of the system across the ovarian cycle: (i) Dyn signalling strength; (ii) NKB signalling strength; (iii) network excitability; and (iv) basal neuronal activity. Estimates for the four parameters in estrus and diestrus are inferred from LH pulse frequency data in estrus and diestrus animals; with or without 5 Hz optic stimulation (Figure 2E&F) (B) System response to low frequency stimulation during estrus and diestrus, using the maximum a-posteriori estimate of the parameter values inferred from the frequency data. (C) Analysis of parameter changes across the cycle. For each of the four parameter $(θ^{i}; i = 1, 2, 3, 4)$ , the diestrus-to-estrus change is defined as the log-ratio between the corresponding parameter values, that is, $l o g_{10} (θ_{e s t r u s}^{i} / θ_{d i e s t r u s}^{i}) .$ Normalised covariance (correlation) matrix of parameter changes reveals negative correlation between changes in NKB signalling strength and network excitability, and positive correlation between Dyn signalling strength and both NKB signalling. Eigen-parameters are visualised as pie charts. The eigen-parameter explaining the least of the variance in the posterior distribution corresponds to the stiffest parameter combination to which the system is most sensitive.

Employing Bayesian inference techniques (see Materials and methods), we sample values for these four parameters, which allow the model to replicate the mean LH frequency we observe experimentally in estrus and diestrus mice with and without 5 Hz optic stimulation (Figure 2E&F). Inspection of the dynamical behaviour of the model, using the identified diestrous parameter values, reveals that in response to optic stimulation in diestrus pulsatile dynamics could die out gradually (i.e. there is a transient period before activity shuts down; see Figure 3B for an illustrative example), which is confirmed by the delayed inhibition of LH pulses we observed experimentally in diestrous mice.

Next, we focus on how the four key parameters change between diestrus and estrus. We measure the change in each parameter using the log-ratio of its estrous to diestrous value and calculate the covariance matrix of these log-ratios from our set of inferred parameter values. We find a positive (linear) correlation between changes in Dyn and NKB signalling strength, and negative (linear) correlation between changes in NKB signalling strength and network excitability (Figure 3C). That is, the model predicts that NKB signalling strength and network excitability are characterised by opposite (in direction) correlations during the transition from diestrus to estrus (one decreasing the other increasing; we note the model predicts that both combinations are possible), whereas NKB and Dyn signalling remain correlated in the same direction (ether increasing or decreasing; we note the model predicts that both combinations are possible). Finally, we apply Principal Component Analysis to study the sensitivity of the system with respect to changes in the four parameters (see Materials and methods). We calculate the principal components in dataset with the inferred parameter changes. Principal components explaining small portions of the variance in the dataset (i.e. principal component with the smallest eigenvalue) correspond to parameter combinations to which the system dynamics are most sensitive (stiff parameter combinations). These combinations are the most critical in terms of regulation as small deviations in how these parameters co-vary result in significant shifts in the system’s dynamics. Interestingly, the principal component capturing the smallest share of the variance is comprised of the parameters controlling NKB signalling, Dyn signalling and network excitability in approximately equal portions, and therefore the model predicts that co-ordinated changes in these three parameters should be critical for the observed changes in system dynamics between diestrus and estrus. Interestingly, the second smallest principal component is largely determined by change in the basal activity parameter, suggesting that basal activity is another independent handle for modulating the system’s dynamics. Taken together our theoretical findings suggest that co-ordinated changes in KNDy signalling as well as changes in KNDy basal activity may be crucial pathways of regulation across the reproductive cycle.

Disrupting glutamatergic transmission in the KNDy population blocks lh pulses

Since KNDy neurons are primarily glutamatergic (Cravo et al., 2011; Nestor et al., 2016; Qiu et al., 2016; Qiu et al., 2018) and synapse to one another (Yip et al., 2015; Qiu et al., 2016), we hypothesise that glutamate transmission should directly affect the levels of excitability within the KNDy network. Hence, we disrupt signalling via glutamate receptors to explore in vivo how network excitability affects the ability of the system to generate and sustain LH pulses across the estous cycle.

First, using Kiss-Cre estrous mice we test whether glutamatergic transmission is necessary for the optogenetic induction of LH pulses. We drive the ARC kisspeptin population using sustained, low-frequency optic stimulation (5 Hz) in the presence of the combined NMDA and AMPA receptor antagonists (AP5 and CNQX, respectively). We find that blocking signalling via glutamate receptors inhibits the capacity of optic stimulation to generate and sustain pulsatile LH secretion (Figure 4A, B, D). This is in agreement with the model prediction that network excitability (ability of KNDy neurons to communicate and synchronise) is critical for sustained pulse generation. The combined AP5 and CNQX in the absence of optic stimulation had no effect (Figure 4C, F).

Figure 4

Download asset Open asset

Effect of NMDA + AMPA receptor antagonists on pulsatile LH secretion in estrus.

Representative examples showing LH secretion in estrous mice in response to optic stimulation (5 Hz blue light, 473 nm, 5 ms pulse width) (A) and optic stimulation combined with the NMDA + AMPA receptor antagonist (bolus ICV injection [12 nmol AP5 +5 nmol CNQX] over 5 min, followed by a continuous infusion [20 nmol AP5 and 10 nmol CNQX] for the remaining 90 min) treatment (B). NMDA + AMPA receptor antagonist alone had no effect (C). (D) Summary showing mean ± SEM LH pulse frequency over the 60 min non-stimulatory period (white bars) and over the subsequent 90 min stimulation period or appropriate non-stimulatory period in presence of. NMDA + AMPA receptor antagonist alone (black bar) in diestrous mice. *Denote LH pulses. ^†p < 0.05 vs pre-stimulation. ^#p < 0.05 compared to antagonist treatment groups; n = 5–6 per group.

Next, we test whether glutamatergic transmission is critical for the endogenous LH pulses observed in diestrus. Treatment of diestrous mice with the combined NMDA and AMPA receptor antagonists resulted in a significant reduction of LH pulse frequency from 2.50 ± 0.29–0.45 ± 0.15 pulses/hour (Figure 5B & D), confirming that the glutamatergic transmission is indeed critical for sustained pulsatility. Moreover, combining NMDA and AMPA receptor antagonist treatment with low frequency optic stimulation (5 Hz) partially restored LH pulsatility to 1.58 ± 0.17 pulses/hour (Figure 5C & D), suggesting low glutamatergic transmission within the KNDy population or from upstream neuronal populations could be offset by other exogenous inputs or elevated basal activity. This finding is in agreement with the model prediction that basal activity and signalling between KNDy neurons are independent pathways of modulating the system’s dynamical behaviour.

Figure 5

Download asset Open asset

Effect of NMDA + AMPA receptor antagonists on pulsatile LH secretion in diestrus.

Representative examples showing pulsatile LH secretion in response to ICV administration of aCSF as control (A), treatment with NMDA + AMPA receptor antagonists (AP5+ CNQX: bolus ICV injection [12 nmol AP5 +5 nmol CNQX] over 5 min, followed by a continuous infusion [20 nmol AP5 and 10 nmol CNQX] for the remaining 90 min) (B) and combined NMDA/AMPA receptor antagonist treatment and sustained optic stimulation (blue light 473 nm, 5 ms pulse width) at 5 Hz (C). (D) Summary showing mean ± SEM LH pulse frequency over the 60 min non-stimulatory period (white bars) and over the subsequent 90 min stimulation period in control mice (grey bar) and mice receiving treatment (black bar). *Denote LH pulses. ^†p < 0.05 vs pre-stimulation. ^#p < 0.05 compared to 5 Hz stimulation plus antagonist treatment and aCSF control groups; n = 5–6 per group.

Discussion

Using optogenetics, we perturbed the GnRH pulse generator at different stages of the ovarian cycle aiming to understand how gonadal steroids modulate key properties of the system. Previous studies have shown how the pulsatile activity generated by the kisspeptin neuronal network is modulated across the estrous cycle (Han et al., 2015; McQuillan et al., 2019). Our data show that the stage of the cycle also has a profound effect on the dynamical response of the kisspeptin population to sustained, low-frequency optic stimulation. Such stimulation triggers acceleration of LH pulses during estrus and deceleration during diestrus. Previously, our mathematical model of the KNDy network has predicted an upper and a lower bifurcation point that determine the system’s range of pulsatile behaviour as the system is driven externally (Voliotis et al., 2019). Our data suggest that the gonadal state plays a critical role in shifting these bifurcation points by modulating key parameters of the system. In particular, during estrous the system is positioned below the lower bifurcation point and optogenetic stimulation of ARC kisspeptin neurons at 5 Hz moves the system across the lower bifurcation point leading to the sudden emergence of pulsatile behaviour. In contrast, during diestrus the system is within the pulsatile regime and optogenetic stimulation of ARC kisspeptin neurons at frequencies greater than 5 Hz moves the system across the upper bifurcation point and its dynamics relax progressively from pulsatile to quiescent. Our data therefore highlight the critical role of gonadal steroids in modulating the dynamical response of the KNDy network to small changes in basal activity of ARC kisspeptin neurons or in how the population processes external perturbations and afferent inputs.

Using our mathematical model of the system, we gained insight into possible mechanisms via which gonadal steroids modulate the dynamic behaviour of the GnRH pulse generator. Based on the differential effect that optic stimulation had on LH pulse frequency in estrous versus diestrous animals, the model predicted that network excitability is an important parameter, which is actively regulated throughout the ovarian cycle. Importantly, KNDy network excitability is most probably co-regulated with parameters controlling the strength of Dyn and NKB signalling as the system transitions between the different phases of the ovarian cycle. In particular, our analysis predicts (i) a strong negative correlation between changes in NKB signalling strength and changes in network excitability, and (ii) a strong positive correlation in changes between NKB and Dyn signalling. We propose that these regulatory relationships ensure robust control of LH frequency over the estrous cycle (Figure 6, Figure 6—figure supplement 1). For example, the positive correlation in the regulation of NKB and Dyn signalling enables robust transition between pulsatile and quiescent dynamics, in contrast to negative correlation that would make the system far more sensitive to the magnitude of the change, that is, changes that are too small or too large could fail to trigger LH pulses (see Figure 6—figure supplement 1).

Figure 6 with 1 supplement see all

Download asset Open asset

Differential effect of perturbations on the dynamics of the pulse generator.

Two-parameter (NKB signalling and network excitability) bifurcation diagram showing the region in the parameter space for which the system exhibits pulsatile dynamics (grey area). Two points (denoted by A and B) illustrate how an increase in NKB signalling or network excitability could have a differential effect on the dynamics of the system. For point A an increase in network excitability or NKB signalling could lead to an increase in the frequency and width of pulses. However, for point B a similar increase leads to pulse inhibition and steady state system dynamics. Furthermore, a negative correlation in how NKB signalling and network excitability co-vary (i.e. aligned with the direction of the pulsatile regime) make the system dynamics less sensitive to small perturbations and enable more robust control over the cyle.

Recent transcriptomic data have revealed that treatment of ovariectomised mice with estradiol reduces expression of NKB and Dyn in KNDy neurons, but increases expression of glutamate transporters (vGlut2; leading to increased glutamate neurotransmission and neuronal excitability in the population) (Qiu et al., 2018). These findings are in line with the regulatory relationships predicted by the model and support the hypotheses that (i) correlated changes in the NKB and Dyn signalling strength should be reflected mainly in the expression levels of the two neuropeptides and consequently in their release availability, and (ii) a key mechanism for regulating network excitability could be through the expression of glutamate transporters. In normally cycling animals, these changes should be driven by the combined action of sex steroids. For example, although postovulatory increase in progesterone is linked to deceleration of LH pulses in the luteal-phase, this inhibitory effect of progesterone is conditional on prior exposure to high estradiol levels (Skinner et al., 1998) Moreover, data from the rhesus monkey support that estradiol and progesterone could play distinct roles in the deceleration of LH pulses from mid-cycle and throughout the luteal-phase (O’Byrne et al., 1991). Our model provides a novel, systems-level understanding of how the genomic changes in the KNDy population link to the dynamic behaviour of the pulse generator. Further experiments will be needed to validate model prediction and link cyclic changes in the sex steroid milieu to genomic pathways dynamically regulating the pulse generator.

The effect of cycle stage on the LH response to sustained optogenetic stimulation is reminiscent of the well documented effect that gonadal steroids have on LH response to various excitatory neurotransmitters and neuropeptides (e.g. NMDA). For instance, investigations in the female monkey revealed an unexplained inhibition of LH in OVX animals following treatment with NMDA, in contrast with the excitatory action of NMDA on LH secretion in the presence of ovarian steroids (Reyes et al., 1990; Reyes et al., 1991). Similar reversal of action on LH dynamics depending on the underlying ovarian steroid milieu has been documented for various other neurotransmitters and NKB receptor agonists (Kalra and Kalra, 1983; Scorticati et al., 2004). Our mathematical model supports that ovarian steroids change key parameters of the KNDy network, which control the dynamic behaviour of the system and its response to perturbations. As an illustration, Figure 6 shows how the dynamic behaviour of the model depends on network excitability and NKB signalling. Since these parameters are governed by gonadal steroids (Qiu et al., 2018; Wang et al., 2018), it is expected that the underlying steroid milieu will also modulate the effect of perturbations on the dynamics of the system. For instance, the effect of stimulating NKB signalling (e.g. via administration of NK3 receptor agonists) or network excitability (e.g. via NMDA administration) could result in inhibition of the pulse generator if the system is already located closer to the right boundary of the pulsatile dynamics region (e.g. point B in Figure 6). Such points correspond to states with high pulse generator activity similar to pulse generator dynamics observed in many animal models after gonadectomy (Reyes et al., 1990; Kinsey-Jones et al., 2012). In contrast, similar perturbations but from a different point in the parameter space, lying closer to the left edge of the pulsatile region (e.g. point A in Figure 6), could result in stimulation of the pulse generator (higher frequency). This illustrative example also highlights that the effect of gonadal steroids on the response of the pulse generator to perturbations is continuous rather than binary, that is, the behaviour of the pulse generator is modulated by the actual continuous levels of gonadal steroids rather than their mere presence of absence. Therefore, the underlying steroid levels could explain seemingly incompatible findings regarding the effect of NKB receptor agonism on LH secretion: ranging from stimulation (Navarro et al., 2011) or inhibition (Sandoval-Guzmán and Rance, 2004; Kinsey-Jones et al., 2012) to no effects (Navarro et al., 2009) in rodents.

Based on our model predictions regarding the importance of network excitability, we set out to uncover the role of this parameter on the dynamic response of the network in vivo. We hypothesised that network excitability should depend, partly at least, on the levels of glutamate signalling as ARC kisspeptin neurons are known to be interconnected (Yip et al., 2015; Qiu et al., 2016) and communicate via glutamate (Cravo et al., 2011; Qiu et al., 2011; Nestor et al., 2016; Qiu et al., 2016). Further evidence that estradiol regulates KNDy neuronal excitability (Qiu et al., 2018) and that cycle stage regulates spontaneous glutamatergic activity of KNDy neurons (Wang and Moenter, 2020) supports the model prediction that network excitability is a critical network property regulated by gonadal steroids. Therefore, to further test this prediction in vivo, we used glutamate receptor (NMDA and AMPA) antagonists to inhibit excitability in the KNDy network. In diestrus animals blocking glutamate receptors (NMDA and AMPA) inhibited LH pulses that were then rescued via low frequency optogenetic stimulation of kisspeptin neurons. Furthermore, in estrus animals, NMDA and AMPA receptor antagonism inhibited the induction of LH pulses via optic stimulation. These experimental findings highlight the complex fine-balanced mechanisms underlying pulse generation by the KNDy network. In particular, limited network excitability within the KNDy population blocks LH pulsatility but this can be mitigated by elevated basal neuronal activity. Similarly, increased basal neuronal activity can induce pulse generation but this effect can be negated by decreased excitability within the neuronal population. There is a caveat, however, as the glutamate receptor antagonists were given by intracerebroventricular injection, and this raises the possibility of having interferred with additional glutamatergic transmission from afferent populations. To the best of our knowledge, there is no evidence of such afferent populations in rodents, although supporting evidence can be found the sheep literature (Merkley et al., 2015). Nevertheless, the possibility of having blocked exogenous glutamatergic inputs does not invalidate our findings, but further supports a key model prediction that basal activation of KNDy neurons (either intrinsic or exogenous) is a critical pathway for modulating the dynamics of the pulse generator. Overall, our model predicts that pulse generation is an emergent property of the KNDy network depending both on single neuron properties such as basal activity and exogeneous activation but also on how the neurons signal and communicate with each other. Our results support this idea, highlighting the critical role of inter-neuronal communication in enabling the population to pulse in synchrony.

Materials and methods

Animals

Adult Kiss-Cre heterozygous transgenic female mice aged between 8 and 14 weeks, 25–30 g, were used for experiments (Yeo et al., 2016). Breeding pairs were obtained from the Department of Physiology, Development and Neuroscience, University of Cambridge, UK and mated in house at King’s College London. Genotyping was performed using a multiplex PCR protocol for detection of heterozygosity for the Kiss-Cre or wild-type allele as previously described (Lass et al., 2020). Only mice with normal estrous cycles were used. Daily vaginal smears were performed for the detection of the estrous and diestrous stages of the ovarian cycle. Mice were singularly housed and provided with food (standard maintenance diet; Special Dietary Services, Wittam, UK) and water ad libitum while being kept under a 12:12 hr light/dark cycle (lights on 0700 h) at 23°C ± 2°C. All animal procedures performed were approved by the Animal Welfare and Ethical Review Body Committee at King’s College London and conducted in accordance with the UK Home Office Regulations.

Surgical procedures

Request a detailed protocol

Stereotaxic injection of AAV9-EF1-dflox-hChR2-(H134R)-mCherry-WPRE-hGH (4.35 × 10¹³ GC/ml; Penn Vector Core; University of Pennsylvania, PA, USA) for targeted expression of channelrhodopsin (ChR2) in ARC kisspeptin neurons was done under aseptic conditions. The mice were anaesthetised using ketamine (Vetalar, 100 mg/kg, i.p.; Pfizer, New York City, NY, USA) and xylazine (Rompun, 10 mg/kg, i.p.; Bayer, Leverkusen, Germany). Kiss-Cre female mice (n = 12) or wilt-type (n = 3) were secured in a Kopf Instruments motorised stereotaxic frame (Kopf Instruments, Tujunga, CA, USA) and surgical procedures on the brain were performed using a Robot Stereotaxy system (Neurostar, Tubingen, Germany). Stereotaxic injection coordinates used to target the ARC were obtained from the mouse brain atlas of Paxinos and Franklin, 2004 (0.25 mm lateral, 1.94 mm posterior to bregma and at a depth of 5.8 mm). A skin incision was made and a small hole was drilled in the skull above the location of the ARC. A 2 μl Hamilton micro-syringe (Esslab, Essex, UK) was attached to the robot stereotaxy and used to inject 0.3 μl of the AAV-construct into the ARC, unilaterally, at a rate of 100 nl/min. After the injection, the needle was left in position for 5 min and then slowly lifted over 1 min. The same coordinates as the injection site were then used to insert a fibre-optic cannula (200 μm, 0.39 NA, 1.25 mm ceramic ferrule; Thorlabs, LTD, Ely, UK); however, a depth of 5.78 mm was reached to ensure the fibre-optic cannula was situated immediately above the injection site. Additionally, an intracerebroventricular (ICV) fluid guide cannulae (26 gauge; Plastics One) targeting the lateral ventricle (coordinates: 1.1 mm lateral, 1.0 mm posterior to bregma and at a depth of 3.0 mm) was chronically implanted. Dental cement (Superbond C&B kit Prestige Dental Products, Bradford UK) was used to fix the cannulae in place and the skin incision was sutured. A one week recovery period was given post-surgery. After this period, the mice were handled daily to acclimatise them to the tail-tip blood sampling procedure (Steyn et al., 2013). Mice were left for 4 weeks to achieve effective opsin expression before experimentation.

Validation of AAV injection site and fibre optic and ICV cannula position

Request a detailed protocol

Once experiments were completed, mice were given a lethal dose of ketamine and transcardially perfused for 5 min with heparinised saline, followed by 10 min of ice-cold 4 % paraformaldehyde (PFA) in phosphate buffer, pH 7.4, for 15 min using a pump (Minipuls; Gilson). Brains were collected immediately and post fixed at 4 °C in 15 % sucrose in 4 % PFA and left to sink. They were then transferred to 30 % sucrose in PBS until they sank. The brains were then snap-frozen on dry ice and stored at –80 °C. Using a cryostat, every third coronal brain section (30 µm) was collected between –1.34 mm and –2.70 mm from bregma and sections were mounted on microscope slides, left to air-dry and cover slipped with ProLong Antifade mounting medium (Molecular Probes, Inc, OR, USA). Verification and evaluation of the injection site was performed using an Axioskop 2 Plus microscope equipped with axiovision 4.7 (Zeiss). One of 12 Kiss-Cre mice failed to show mCherry fluorescence in the ARC and was excluded from the analysis.

Experimental design and blood sampling for LH measurement

Request a detailed protocol

For measurement of LH pulsatility during optogenetic stimulation, the tip of the mouse’s tail was removed with a sterile scalpel for tail-tip blood sampling (Czieselsky et al., 2016). The chronically implanted fibre-optic cannula was attached to a multimode fibre-optic rotary joint patch cables (Thorlabs) via a ceramic mating sleeve. This allows for freedom of movement and blue light delivery (473 nm wavelength) using a Grass SD9B stimulator controlled DPSS laser (Laserglow Technologies) during optogenetic stimulation.

The experimental protocol involved an hour long acclimatisation period, followed by 2.5 hr of blood sampling, where 5 μl of blood was collected every 5 min. For estrous and diestrous mice, optic stimulation was initiated after 1 hr of control blood sampling and was sustained for 1.5 hr. Optic stimulation was delivered as 5 ms pulses of light at 5 Hz with the laser intensity measured at the tip of the fibre-optic patch cable set to 5 mW (Voliotis et al., 2019). Additionally, in separate experiments, diestrous mice were optically stimulated at five or 15 Hz for 2.5 hr, that is entire blood sampling period. Control mice (in estrus or diestrus) received no optic stimulation. Wild-type mice (estrus and disetrus) received 5 Hz optic stimulation to verify that our optic stimulation protocol had no undesirable effects on LH secretion.

Neuropharmacological manipulation of glutamatergic signalling was performed using a combination of NMDA (AP5, Tocris, Abingdon, UK) and AMPA (CNQX, Alpha Aesar, Heysham, UK) receptor antagonist treatment with or without simultaneous optogenetic stimulation. The animals were prepared for optogenetic experimentation as described above with additional preparation of the ICV injection cannula. Immediately after connection of the fire-optic cannula, the ICV injection cannula with extension tubing, preloaded with drug solution (AP5 and CNQX dissolved in artificial CSF) or artificial CSF alone as control, was inserted into the guide cannula. The extension tubing, reaching outside of the cage, was connected to a 10 μl Hamilton syringe mounted in an automated pump (Harvard Apparatus) to allow for remote micro-infusion without disturbing the animals during experimentation. After a 55 min control blood sampling period, as described above, and 5 min before the onset of optic stimulation, a bolus ICV injection of drug solution (12 nmol AP5 and 5 nmol CNQX in 2.3 μl) was given over 5 min, followed by a continuous infusion (20 nmol AP5 and 10 nmol CNQX in 5.6 μl) for the remaining 90 min of experimentation. Artificial CSF controls, with or without optic stimulation, received the same ICV fluid regime. When no optic stimulation was applied the same ICV administration and blood sampling regimen described was applied. Stimulation and non-stimulation protocols were implemented in random order for Kiss-Cre mice.

The blood samples were snap-frozen on dry ice and stored at –80 °C until processed. In-house LH enzyme-linked immunosorbent assay (LH ELISA) similar to that described by Steyn et al. was used for processing of the mouse blood samples (Steyn et al., 2013). The mouse LH standard (AFP- 5,306 A; NIDDK-NHPP) was purchased from Harbor-UCLA along with the primary antibody (polyclonal antibody, rabbit LH antiserum, AFP240580Rb; NIDDK-NHPP). The secondary antibody (donkey anti-rabbit IgG polyclonal antibody [horseradish peroxidase]; NA934) was from VWR International. Validation of the LH ELISA was done in accordance with the procedure described in Steyn et al., 2013 derived from protocols defined by the International Union of Pure and Applied Chemistry. Serially diluted mLH standard replicates were used to determine the linear detection range. Nonlinear regression analysis was performed using serially diluted mLH standards of known concentration to create a standard curve for interpolating the LH concentration in whole blood samples, as described previously (Voliotis et al., 2019). The assay sensitivity was 0.031 ng/mL, with intra- and inter-assay coefficients of variation of 4.6% and 10.2%, respectively.

LH pulse detection and statistical analysis

Request a detailed protocol

Dynpeak algorithm was used for the detection of LH pulses (Vidal et al., 2012). The differential effect of optogenetic stimulation on LH pulsatility in estrus and diestrus was determined by looking at the frequency of LH pulses. For mice in estrus and for the neuropharmacological experiments, the mean ± SEM of LH pulses per hour were compared between the 60 min pre-stimulation/drug delivery control period and subsequent 90 min stimulation period. For mice in diestrus, the mean ± SEM of LH pulses per hour were compared between controls, 5 Hz and 15 Hz treatment groups, as optic stimulation was applied from the beginning of blood sample period. No optic stimulation was applied to control animals, however the same time points were compared. The frequency of LH pulses in the 90 min optic stimulation/drug delivery period was also compared between treatment groups. Mann-Whitney Rank Sum test was used to access LH frequency differences between groups and determine statistical significance (p < 0.05). LH data publicly available from http://doi.org/doi:10.18742/RDM01-750.

Mathematical model of the KNDy network

Request a detailed protocol

We used a modified version of our previously published mathematical model of the KNDy network (Voliotis et al., 2019). The model offers a high-level overview of the system, wrapping many biophysical details into a coarse grain description for the sake of simplicity and brevity. Importantly such a parsimonious model fits best to the high-level, holistic in vivo approach we use to study the system. Briefly, the model describes the ARC kisspeptin population in terms of three variables: $D$ , the average concentration of Dyn secreted by the population; $\overset{´}{N}$ , the average concentration of NKB secreted by the population; and $\overset{´}{v}$ , the average firing activity of the population, measured in spikes/min. The variables obey the following set of coupled ordinary differential equations (ODEs):

\frac{d D}{d t} = f_{D} (v) - d_{D} D;

\frac{d N}{d t} = f_{N} (v, D) - d_{N} N;

\frac{d v}{d t} = f_{v} (v, N) - d_{v} v .

Parameters $d_{D}$ , $d_{N}$ and $d_{v}$ control the characteristic timescale of each variable. The model describes Dyn and NKB secretion as independent processes based on the observation that Dyn and NKB are packaged in separate vesicles (Murakawa et al., 2016). The secretion rates of the two neuropeptides are given by:

f_{D} (\overset{´}{v}) = k_{D} \frac{{\overset{´}{v}}^{2}}{{\overset{´}{v}}^{2} + K_{v}^{2}};

f_{N} (\overset{´}{v}, \overset{´}{D}) = k_{N} \frac{{\overset{´}{v}}^{2}}{{\overset{´}{v}}^{2} + K_{v}^{2}} \frac{K_{D}^{2}}{{\overset{´}{D}}^{2} + K_{D}^{2}} .

In the equations above neuronal activity ( $\overset{´}{v}$ ) stimulates secretion of both neuropeptides, and Dyn represses NKB secretion. The maximum secretion rate for the two neuropeptides is controlled by parameters $k_{D}$ and $k_{N}$ and we refer to these parameters as the strength of Dyn and NKB singalling respectively. Furthermore, we assume that distinct modes of Dyn and NKB regulation (e.g. in terms of their synthesis and depletion rate, intracellular transport, packaging dynamics) are reflected in their secretion rate and therefore model Dyn and NKB regulation throughout the estrous cycle as changes of parameters $k_{D}$ and $k_{N}$ . The effector levels at which saturation occurs are controlled via parameters $K_{v}$ and $K_{D} .$ Here, we are interested in investigating the effect of network excitability on the dynamics therefore we modify the equation for the neuronal activity, $\overset{´}{v}$ , by setting:

f_{v} (\overset{´}{v}, \overset{´}{N}) = v_{0} \frac{1 - e x p (- I)}{1 + e x p (- I)}; I = - l o g \frac{1 - b}{1 + b} + k_{v} (e + \frac{N^{2}}{N^{2} + K_{N}^{2}}) \overset{´}{v} .

Here, we have introduced parameter $k_{v}$ capturing the intrinsic network excitability, that relates to the strength of the synaptic connections between KNDy neurons that are essential for synchronising neural activity and enabling pulse generation. We note that this parameter will be directly affected by any neuromodulator (including Glutamate) that affects KNDy activity as well as by processes that affect KNDy neurons’ synaptic density. Furthermore, parameter $b$ controls the basal neuronal activitation of the population, which could stem from synaptic noise or afferent inputs (extrinsic to the network). We assume that both $k_{v}$ and $b$ are regulated throughout the estrous cycle. Finally, $v_{0}$ is the maximum rate at which the firing rate increases in response to synaptic inputs $I$ . Note the stimulatory effect of NKB (which is secreted at the presynaptic terminal) on neuronal activity (Qiu et al., 2016). The full list of model parameters is given in Table 1.

Table 1

Model parameters values.

No	Parameter	Description	Value	Ref.
1	$d_{D}$	Dyn degradation rate	0.25 min^–1	Voliotis et al., 2019
2	$d_{N}$	NKB degradation rate	0.25 min^–1	Voliotis et al., 2019
3	$d_{v}$	Firing rate reset rate	10 min^–1	Qiu et al., 2016
4	$k_{D}$	Dyn singalling strength	inferred
5	$k_{N}$	NKB signalling strength	inferred	Ruka et al., 2016
6	$k_{v}$	Network excitability	inferred
7	$v_{0}$	Maximum rate of neuronal activity increase	30000 spikes min^–2	Qiu et al., 2016
8	$K_{D}$	Dyn IC₅₀	0.3 nM	Yasuda et al., 1993
9	$K_{N}$	NKB EC₅₀	32 nM	Seabrook et al., 1995
10	$K_{v}$	Firing rate for half-maximal NKB and Dyn secretion	1200 spikes min^–1	Dutton and Dyball, 1979
11	$b$	Basal activity	inferred
12	e	NKB-independent contribution to network excitability	inferred

Parameter inference

Request a detailed protocol

We used Approximate Bayesian Computation (ABC) based on sequential Monte Carlo (SMC) (Toni et al., 2009) to infer four key model parameters (Dyn signaling strength, $k_{D}$ ; NKB signalling strength, $k_{N}$ ; network excitability $k_{v}$ ; and basal activity, $b$ ) in the estrous and diestrous phase of the ovarian cycle. For inference we used the average LH inter-pulse interval observed in four different settings: estrous animals without optic stimulation ( $I_{E}$ ) and with 5 Hz optic stimulation ( $I_{E + 5 H z}$ ); diestrus animals without optic stimulation ( $I_{D}$ ) and with 5 Hz optic stimulation ( $I_{D + 5 H z}$ ). Model simulations were generated in Matlab using function ode45 under the four different settings for 6000 min and by calculating the frequencies after discarding the initial 1000 min. The following discrepancy function was used to compare simulated, $D^{} = (I_{E}^{*}, I_{E + 5 H z}^{*}, I_{D}^{*}, I_{D + 5 H z}^{*})$ , and experimental, $D = (I_{E}, I_{E + 5 H z}, I_{D}, I_{D + 5 z})$ , data:

d (D, D^{*}) = \sum_{i = 1}^{4} | D_{i} - D_{i}^{*} |

Furthermore, for the ABC SMC algorithm the size of the particle population was set to 500 and the algorithm was run for $T = 4$ populations with corresponding tolerance levels $ε_{i} = 54 - 2 i, i = 0, \dots, 26$ . Log-uniform prior distributions were used to explore the behaviour of the model under wide parameter ranges: $\log_{10} (k_{D}) \sim U n i f o r m (- 3, 3)$ ; $l o g_{10} (k_{N}) \sim U n i f o r m (- 3, 3)$ ; $l o g_{10} (k_{v}) \sim U n i f o r m (- 3, 3)$ ; and $\log_{10} (b) \sim U n i f o r m (- 3, 0)$ . All remaining parameters were fixed to values found in the literature (see Table 1). For each parameter, an independent log₁₀-normal perturbation kernel with variance 0.05 was used. Matlab code can be found at https://git.exeter.ac.uk/mv286/kndy-parameter-inference (Voliotis, 2021, copy archived at swh:1:rev:f5585f5f09d62a2569d21def31dc779774bc2476)

Sensitivity and principal component analysis

Request a detailed protocol

We used principal component analysis to study the sensitivity of the system with respect to changes in the four inferred parameters (Toni et al., 2009). We calculate the principal components in the dataset (sampled posterior distribution) of the inferred parameter changes. Principal component analysis produces a set of linearly uncorrelated eigen-parameters explaining the variance of the inferred changes (in the sampled posterior distribution.) These eigen-parameters are linear weighted combinations of the initial parameters. The eigen-parameter explaining the least of the variance in the posterior distribution corresponds to the stiffest parameter combination. That is small deviations from the inferred way these parameters co-vary would lead to changes in the model behaviour that make it incompatible with the data.

Data availability

The data and the code are publicly available via the following open access repositories: http://doi.org/doi:10.18742/RDM01-750 https://git.exeter.ac.uk/mv286/kndy-parameter-inference.git (copy archived at https://archive.softwareheritage.org/swh:1:rev:f5585f5f09d62a2569d21def31dc779774bc2476).

The following data sets were generated

1. de Burgh R
(2021) figshare
Modulation of pulsatile GnRH dynamics along the reproductive cycle - the role of excitability within the arcuate kisspeptin network.

https://doi.org/10.18742/RDM01-750

References

(1993) Estradiol modulates the LH release response to N-methyl-D-aspartate in adult female rats: studies on hypothalamic luteinizing hormone-releasing hormone and neurotransmitter release
Neuroendocrinology 57:710–715.

https://doi.org/10.1159/000126429
- PubMed
- Google Scholar
1. Bonavera JJ
2. Sahu A
3. Kalra SP
4. Kalra PS
(1994) The hypothalamic peptides, beta-endorphin, neuropeptide K and interleukin-1 beta, and the opiate morphine, enhance the excitatory amino acid-induced LH release under the influence of gonadal steroids
Journal of Neuroendocrinology 6:557–564.

https://doi.org/10.1111/j.1365-2826.1994.tb00619.x
- PubMed
- Google Scholar
1. Brann DW
2. Mahesh VB
(1992) Excitatory amino acid regulation of gonadotropin secretion: modulation by steroid hormones
The Journal of Steroid Biochemistry and Molecular Biology 41:847–850.

https://doi.org/10.1016/0960-0760(92)90436-m
- PubMed
- Google Scholar
1. Christian CA
2. Moenter SM
(2010) The neurobiology of preovulatory and estradiol-induced gonadotropin-releasing hormone surges
Endocrine Reviews 31:544–577.

https://doi.org/10.1210/er.2009-0023
- PubMed
- Google Scholar
1. Cravo RM
2. Margatho LO
3. Osborne-Lawrence S
4. Donato J Jr
5. Atkin S
6. Bookout AL
7. Rovinsky S
8. Frazão R
9. Lee CE
10. Gautron L
11. Zigman JM
12. Elias CF
(2011) Characterization of Kiss1 neurons using transgenic mouse models
Neuroscience 173:37–56.

https://doi.org/10.1016/j.neuroscience.2010.11.022
- PubMed
- Google Scholar
(2016) Pulse and Surge Profiles of Luteinizing Hormone Secretion in the Mouse
Endocrinology 157:4794–4802.

https://doi.org/10.1210/en.2016-1351
- PubMed
- Google Scholar
1. Dutton A
2. Dyball RE
(1979) Phasic firing enhances vasopressin release from the rat neurohypophysis
The Journal of Physiology 290:433–440.

https://doi.org/10.1113/jphysiol.1979.sp012781
- PubMed
- Google Scholar
1. Goodman RL
2. Holaskova I
3. Nestor CC
4. Connors JM
5. Billings HJ
6. Valent M
7. Lehman MN
8. Hileman SM
(2011) Evidence that the arcuate nucleus is an important site of progesterone negative feedback in the ewe
Endocrinology 152:3451–3460.

https://doi.org/10.1210/en.2011-0195
- PubMed
- Google Scholar
(2015) Selective optogenetic activation of arcuate kisspeptin neurons generates pulsatile luteinizing hormone secretion
PNAS 112:13109–13114.

https://doi.org/10.1073/pnas.1512243112
- PubMed
- Google Scholar
1. Kalra SP
2. Kalra PS
(1983) Neural regulation of luteinizing hormone secretion in the rat
Endocrine Reviews 4:311–351.

https://doi.org/10.1210/edrv-4-4-311
- PubMed
- Google Scholar
1. Kelly MJ
2. Zhang C
3. Qiu J
4. Rønnekleiv OK
(2013) Pacemaking kisspeptin neurons
Experimental Physiology 98:1535–1543.

https://doi.org/10.1113/expphysiol.2013.074559
- PubMed
- Google Scholar
(2012) The inhibitory effects of neurokinin B on GnRH pulse generator frequency in the female rat
Endocrinology 153:307–315.

https://doi.org/10.1210/en.2011-1641
- PubMed
- Google Scholar
1. Lass G
2. Li XF
3. de Burgh RA
4. He W
5. Kang Y
6. Hwa-Yeo S
7. Sinnett-Smith LC
8. Manchishi SM
9. Colledge WH
10. Lightman SL
11. O’Byrne KT
(2020) Optogenetic stimulation of kisspeptin neurones within the posterodorsal medial amygdala increases luteinising hormone pulse frequency in female mice
Journal of Neuroendocrinology 32:e12823.

https://doi.org/10.1111/jne.12823
- PubMed
- Google Scholar
(2019) GnRH Pulse Generator Activity Across the Estrous Cycle of Female Mice
Endocrinology 160:1480–1491.

https://doi.org/10.1210/en.2019-00193
- PubMed
- Google Scholar
(2015) Evidence for Changes in Numbers of Synaptic Inputs onto KNDy and GnRH Neurones during the Preovulatory LH Surge in the Ewe
Journal of Neuroendocrinology 27:624–635.

https://doi.org/10.1111/jne.12293
- PubMed
- Google Scholar
1. Moore AM
2. Coolen LM
3. Porter DT
4. Goodman RL
5. Lehman MN
(2018) KNDy Cells Revisited
Endocrinology 159:3219–3234.

https://doi.org/10.1210/en.2018-00389
- PubMed
- Google Scholar
(2016) Immunoelectron microscopic observation of the subcellular localization of kisspeptin, neurokinin B and dynorphin A in KNDy neurons in the arcuate nucleus of the female rat
Neuroscience Letters 612:161–166.

https://doi.org/10.1016/j.neulet.2015.12.008
- PubMed
- Google Scholar
(2009) Regulation of gonadotropin-releasing hormone secretion by kisspeptin/dynorphin/neurokinin B neurons in the arcuate nucleus of the mouse
The Journal of Neuroscience 29:11859–11866.

https://doi.org/10.1523/JNEUROSCI.1569-09.2009
- PubMed
- Google Scholar
(2011) Interactions between kisspeptin and neurokinin B in the control of GnRH secretion in the female rat
American Journal of Physiology. Endocrinology and Metabolism 300:202–210.

https://doi.org/10.1152/ajpendo.00517.2010
- PubMed
- Google Scholar
1. Nestor CC
2. Qiu J
3. Padilla SL
4. Zhang C
5. Bosch MA
6. Fan W
7. Aicher SA
8. Palmiter RD
9. Rønnekleiv OK
10. Kelly MJ
(2016) Optogenetic Stimulation of Arcuate Nucleus Kiss1 Neurons Reveals a Steroid-Dependent Glutamatergic Input to POMC and AgRP Neurons in Male Mice
Molecular Endocrinology 30:630–644.

https://doi.org/10.1210/me.2016-1026
- PubMed
- Google Scholar
(1991) Radiotelemetric monitoring of hypothalamic gonadotropin-releasing hormone pulse generator activity throughout the menstrual cycle of the rhesus monkey
Endocrinology 129:1207–1214.

https://doi.org/10.1210/endo-129-3-1207
- PubMed
- Google Scholar
Book
1. Paxinos G
2. Franklin KBJ
(2004)
The Mouse Brain in Stereotaxic Coordinates

Elsevier Academic Press.
- Google Scholar
1. Qiu J
2. Fang Y
3. Bosch MA
4. Rønnekleiv OK
5. Kelly MJ
(2011) Guinea pig kisspeptin neurons are depolarized by leptin via activation of TRPC channels
Endocrinology 152:1503–1514.

https://doi.org/10.1210/en.2010-1285
- PubMed
- Google Scholar
1. Qiu J
2. Nestor CC
3. Zhang C
4. Padilla SL
5. Palmiter RD
6. Kelly MJ
7. Rønnekleiv OK
(2016) High-frequency stimulation-induced peptide release synchronizes arcuate kisspeptin neurons and excites GnRH neurons
eLife 5:e16246.

https://doi.org/10.7554/eLife.16246
- PubMed
- Google Scholar
1. Qiu J
2. Rivera HM
3. Bosch MA
4. Padilla SL
5. Stincic TL
6. Palmiter RD
7. Kelly MJ
8. Rønnekleiv OK
(2018) Estrogenic-dependent glutamatergic neurotransmission from kisspeptin neurons governs feeding circuits in females
eLife 7:e35656.

https://doi.org/10.7554/eLife.35656
- PubMed
- Google Scholar
(1990) Unexpected inhibitory action of N-methyl-D,L-aspartate or luteinizing hormone release in adult ovariectomized rhesus monkeys: a role of the hypothalamic-adrenal axis
Endocrinology 127:724–729.

https://doi.org/10.1210/endo-127-2-724
- PubMed
- Google Scholar
1. Reyes A
2. Xia LN
3. Ferin M
(1991) Modulation of the effects of N-methyl-D,L-aspartate on luteinizing hormone by the ovarian steroids in the adult rhesus monkey
Neuroendocrinology 54:405–411.

https://doi.org/10.1159/000125921
- PubMed
- Google Scholar
(2016) Both Estrogen and Androgen Modify the Response to Activation of Neurokinin-3 and κ-Opioid Receptors in Arcuate Kisspeptin Neurons From Male Mice
Endocrinology 157:752–763.

https://doi.org/10.1210/en.2015-1688
- PubMed
- Google Scholar
1. Sandoval-Guzmán T
2. Rance NE
(2004) Central injection of senktide, an NK3 receptor agonist, or neuropeptide Y inhibits LH secretion and induces different patterns of Fos expression in the rat hypothalamus
Brain Research 1026:307–312.

https://doi.org/10.1016/j.brainres.2004.08.026
- PubMed
- Google Scholar
(2004) The inhibitory effect of anandamide on luteinizing hormone-releasing hormone secretion is reversed by estrogen
PNAS 101:11891–11896.

https://doi.org/10.1073/pnas.0404366101
- PubMed
- Google Scholar
(1995) Pharmacology of tachykinin receptors on neurones in the ventral tegmental area of rat brain slices
European Journal of Pharmacology 273:113–119.

https://doi.org/10.1016/0014-2999(94)00681-v
- PubMed
- Google Scholar
1. Skinner DC
2. Evans NP
3. Delaleu B
4. Goodman RL
5. Bouchard P
6. Caraty A
(1998) The negative feedback actions of progesterone on gonadotropin-releasing hormone secretion are transduced by the classical progesterone receptor
PNAS 95:10978–10983.

https://doi.org/10.1073/pnas.95.18.10978
- PubMed
- Google Scholar
1. Steyn FJ
2. Wan Y
3. Clarkson J
4. Veldhuis JD
5. Herbison AE
6. Chen C
(2013) Development of a methodology for and assessment of pulsatile luteinizing hormone secretion in juvenile and adult male mice
Endocrinology 154:4939–4945.

https://doi.org/10.1210/en.2013-1502
- PubMed
- Google Scholar
Book
1. Strogatz SH
(2018) Nonlinear Dynamics and Chaos with Student Solutions Manual: With Applications to Physics, Biology, Chemistry, and Engineering
CRC press.

https://doi.org/10.1201/9780429399640
- Google Scholar
1. Toni T
2. Welch D
3. Strelkowa N
4. Ipsen A
5. Stumpf MPH
(2009) Approximate Bayesian computation scheme for parameter inference and model selection in dynamical systems
Journal of the Royal Society, Interface 6:187–202.

https://doi.org/10.1098/rsif.2008.0172
- PubMed
- Google Scholar
(2017) Long-Term Recordings of Arcuate Nucleus Kisspeptin Neurons Reveal Patterned Activity That Is Modulated by Gonadal Steroids in Male Mice
Endocrinology 158:3553–3564.

https://doi.org/10.1210/en.2017-00382
- PubMed
- Google Scholar
1. Vidal A
2. Zhang Q
3. Médigue C
4. Fabre S
5. Clément F
(2012) DynPeak: an algorithm for pulse detection and frequency analysis in hormonal time series
PLOS ONE 7:e39001.

https://doi.org/10.1371/journal.pone.0039001
- PubMed
- Google Scholar
(2019) The Origin of GnRH Pulse Generation: An Integrative Mathematical-Experimental Approach
The Journal of Neuroscience 39:9738–9747.

https://doi.org/10.1523/JNEUROSCI.0828-19.2019
- PubMed
- Google Scholar
Software
1. Voliotis M
(2021) KNDY parameter inference, version swh:1:rev:f5585f5f09d62a2569d21def31dc779774bc2476
Software Heritage.

https://archive.softwareheritage.org/swh:1:dir:bd6fd08a3e4a03e2517bb45ec6312ce6e374d84f;origin=https://git.exeter.ac.uk/mv286/kndy-parameter-inference;visit=swh:1:snp:4bccc258033b8838a7e52d65dbcfded60dc5092a;anchor=swh:1:rev:f5585f5f09d62a2569d21def31dc779774bc2476
(2018) Glutamatergic Transmission to Hypothalamic Kisspeptin Neurons Is Differentially Regulated by Estradiol through Estrogen Receptor α in Adult Female Mice
The Journal of Neuroscience 38:1061–1072.

https://doi.org/10.1523/JNEUROSCI.2428-17.2017
- PubMed
- Google Scholar
1. Wang L
2. Moenter SM
(2020) Differential Roles of Hypothalamic AVPV and Arcuate Kisspeptin Neurons in Estradiol Feedback Regulation of Female Reproduction
Neuroendocrinology 110:172–184.

https://doi.org/10.1159/000503006
- PubMed
- Google Scholar
1. Yasuda K
2. Raynor K
3. Kong H
4. Breder CD
5. Takeda J
6. Reisine T
7. Bell GI
(1993) Cloning and functional comparison of kappa and delta opioid receptors from mouse brain
PNAS 90:6736–6740.

https://doi.org/10.1073/pnas.90.14.6736
- PubMed
- Google Scholar
1. Yeo SH
2. Kyle V
3. Morris PG
4. Jackman S
5. Sinnett-Smith LC
6. Schacker M
7. Chen C
8. Colledge WH
(2016) Visualisation of Kiss1 Neurone Distribution Using a Kiss1-CRE Transgenic Mouse
Journal of Neuroendocrinology 28:12435.

https://doi.org/10.1111/jne.12435
- PubMed
- Google Scholar
(2015) Conditional Viral Tract Tracing Delineates the Projections of the Distinct Kisspeptin Neuron Populations to Gonadotropin-Releasing Hormone (GnRH) Neurons in the Mouse
Endocrinology 156:2582–2594.

https://doi.org/10.1210/en.2015-1131
- PubMed
- Google Scholar

Article and author information

Author details

Margaritis Voliotis

Department of Mathematics and Living Systems Institute, College of Engineering, Mathematics and Physical Sciences, University of Exeter, Exeter, United Kingdom

Contribution
Conceptualization, Investigation, Methodology, Software, Writing – original draft

Contributed equally with
Xiao Feng Li and Ross Alexander De Burgh

For correspondence
m.voliotis@exeter.ac.uk

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-6488-7198
Xiao Feng Li

Department of Women and Children’s Health, School of Life Course Sciences, King’s College London, London, United Kingdom

Contribution
Conceptualization, Data curation, Investigation, Writing – review and editing

Contributed equally with
Margaritis Voliotis and Ross Alexander De Burgh

Competing interests
No competing interests declared
Ross Alexander De Burgh

Department of Women and Children’s Health, School of Life Course Sciences, King’s College London, London, United Kingdom

Contribution
Data curation, Investigation

Contributed equally with
Margaritis Voliotis and Xiao Feng Li

Competing interests
No competing interests declared
Geffen Lass

Department of Women and Children’s Health, School of Life Course Sciences, King’s College London, London, United Kingdom

Contribution
Investigation

Competing interests
No competing interests declared
Deyana Ivanova

Department of Women and Children’s Health, School of Life Course Sciences, King’s College London, London, United Kingdom

Contribution
Data curation, Investigation

Competing interests
No competing interests declared
Caitlin McIntyre

Department of Women and Children’s Health, School of Life Course Sciences, King’s College London, London, United Kingdom

Contribution
Data curation, Investigation

Competing interests
No competing interests declared
Kevin O'Byrne

Department of Women and Children’s Health, School of Life Course Sciences, King’s College London, London, United Kingdom

Contribution
Conceptualization, Investigation, Writing – original draft

Competing interests
No competing interests declared
Krasimira Tsaneva-Atanasova

Department of Mathematics and Living Systems Institute, College of Engineering, Mathematics and Physical Sciences, University of Exeter, Exeter, United Kingdom

Contribution
Conceptualization, Investigation, Writing – original draft, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-6294-7051

Funding

Engineering and Physical Sciences Research Council (EP/N014391/1)

Margaritis Voliotis
Krasimira Tsaneva-Atanasova

Biotechnology and Biological Sciences Research Council (BB/S000550/1)

Margaritis Voliotis
Xiao Feng Li
Kevin O'Byrne
Krasimira Tsaneva-Atanasova

Biotechnology and Biological Sciences Research Council (BB/S001255/1)

Margaritis Voliotis
Xiao Feng Li
Kevin O'Byrne
Krasimira Tsaneva-Atanasova

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

The authors gratefully acknowledge the financial support of the EPSRC via grant EP/N014391/1 (KTA and MV), and BBSRC via grants BB/S000550/1 and BB/S001255/1 (KTA, KOB, MV, XFL).

Ethics

All animal procedures performed were approved by the Animal Welfare and Ethical Review Body Committee at King's College London (PP4006193 ) and conducted in accordance with the UK Home Office Regulations.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

1,021

views
246

downloads
25

citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Margaritis Voliotis
Xiao Feng Li
Ross Alexander De Burgh
Geffen Lass
Deyana Ivanova
Caitlin McIntyre
Kevin O'Byrne
Krasimira Tsaneva-Atanasova

(2021)

Modulation of pulsatile GnRH dynamics across the ovarian cycle via changes in the network excitability and basal activity of the arcuate kisspeptin network

eLife 10:e71252.

https://doi.org/10.7554/eLife.71252