Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Blog
    • Collections
    • Podcast
  • TOPICS
    • Cognition and Behavior
    • Development
    • Disorders of the Nervous System
    • History, Teaching and Public Awareness
    • Integrative Systems
    • Neuronal Excitability
    • Novel Tools and Methods
    • Sensory and Motor Systems
  • ALERTS
  • FOR AUTHORS
  • ABOUT
    • Overview
    • Editorial Board
    • For the Media
    • Privacy Policy
    • Contact Us
    • Feedback
  • SUBMIT

User menu

Search

  • Advanced search
eNeuro

eNeuro

Advanced Search

 

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Blog
    • Collections
    • Podcast
  • TOPICS
    • Cognition and Behavior
    • Development
    • Disorders of the Nervous System
    • History, Teaching and Public Awareness
    • Integrative Systems
    • Neuronal Excitability
    • Novel Tools and Methods
    • Sensory and Motor Systems
  • ALERTS
  • FOR AUTHORS
  • ABOUT
    • Overview
    • Editorial Board
    • For the Media
    • Privacy Policy
    • Contact Us
    • Feedback
  • SUBMIT
PreviousNext
Research ArticleResearch Article: New Research, Neuronal Excitability

Reciprocal Changes in Voltage-Gated Potassium and Subthreshold Inward Currents Help Maintain Firing Dynamics of AVPV Kisspeptin Neurons during the Estrous Cycle

J. Rudolph Starrett, R. Anthony DeFazio and Suzanne M. Moenter
eNeuro 12 August 2021, 8 (5) ENEURO.0324-21.2021; DOI: https://doi.org/10.1523/ENEURO.0324-21.2021
J. Rudolph Starrett
1Department of Molecular and Integrative Physiology, University of Michigan, Ann Arbor, MI 48109
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for J. Rudolph Starrett
R. Anthony DeFazio
1Department of Molecular and Integrative Physiology, University of Michigan, Ann Arbor, MI 48109
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for R. Anthony DeFazio
Suzanne M. Moenter
1Department of Molecular and Integrative Physiology, University of Michigan, Ann Arbor, MI 48109
2Department of Obstetrics and Gynecology, University of Michigan, Ann Arbor, MI 48109
3Department of Internal Medicine, University of Michigan, Ann Arbor, MI 48109
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Suzanne M. Moenter
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Kisspeptin-expressing neurons in the anteroventral-periventricular nucleus (AVPV) are part of a neural circuit generating the gonadotropin-releasing hormone (GnRH) surge. This process is estradiol-dependent and occurs on the afternoon of proestrus in female mice. On proestrus, AVPV kisspeptin neurons express more kisspeptin and exhibit higher frequency action potentials and burst firing compared with diestrus, which is characterized by a pulsatile rather than a prolonged surge of GnRH secretion. We hypothesized changes in voltage-gated potassium conductances shape activity profiles of these cells in a cycle-dependent manner. Whole-cell voltage-clamp recordings of GFP-identified AVPV kisspeptin neurons in brain slices from diestrous and proestrous mice revealed three subcomponents of the voltage-sensitive K+ current: fast-transient slow-transient, and residual. During proestrus, the V50 of inactivation of the fast-transient current was depolarized and the amplitude of the slow-transient component was reduced compared with diestrus; the residual component was consistent across both stages. Computational models were fit to experimental data, including published estrous-cycle effects on other voltage-gated currents. Computer simulations suggest proestrus-typical K+ currents are suppressive compared with diestrus. Interestingly, larger T-type, persistent-sodium, and hyperpolarization-activated currents during proestrus compensate for this suppressive effect while also enabling postinhibitory rebound bursting. These findings suggest modulation of voltage-gated K+ and multiple subthreshold depolarizing currents across the negative to positive feedback transition maintain AVPV kisspeptin neuron excitability in response to depolarizing stimuli. These changes also enable firing in response to hyperpolarization, providing a net increase in neuronal excitability, which may contribute to activation of this population leading up to the preovulatory GnRH surge.

  • AVPV
  • estrous cycle
  • excitability
  • kisspeptin
  • positive feedback
  • potassium

Significance Statement

Gonadotropin-releasing hormone (GnRH) neurons provide the central signal to initiate ovulation by releasing a surge of hormone. GnRH neurons are regulated by other cells including those expressing kisspeptin, a potent stimulator of GnRH secretion. Kisspeptin neurons in the anteroventral-periventricular nucleus (AVPV) express more kisspeptin and become more active during the afternoon of proestrus, the phase of the rodent estrous (reproductive) cycle when the GnRH surge occurs. We found voltage-dependent potassium currents in AVPV kisspeptin neurons change with phase of the estrous cycle. Firing simulations indicated these changes are suppressive if occurring in isolation. But proestrous-typical increases in subthreshold depolarizing currents overcome this suppression and promote greater excitability by increasing rebound firing, possibly contributing to the preovulatory activation of this system.

Introduction

The anteroventral-periventricular (AVPV) nucleus is a critical site involved in the regulation of female fertility, specifically the control of the estradiol-dependent process of ovulation (Kalra and McCann, 1975; Goodman, 1978; Wiegand and Terasawa, 1982; Petersen and Barraclough, 1989; Petersen et al., 1989). The AVPV contains neurons that express both estrogen receptor (ER)α and kisspeptin (Simerly, 1998; Smith et al., 2005), a neuropeptide that acts at gonadotropin-releasing hormone (GnRH) neurons to potently induce GnRH neuron activity and hormone secretion (Han et al., 2005; Messager et al., 2005; Pielecka-Fortuna et al., 2008; Zhang et al., 2008; Caraty et al., 2013). Ovulation is initiated by positive feedback actions of estradiol at the central and pituitary levels (Döcke and Dörner, 1965; Sarkar et al., 1976; Moenter et al., 1991). In a prevailing theory, sustained elevation of circulating estradiol during the late follicular phase (the cycle day of proestrus in rodents) exerts positive feedback effects via ERα that drive increased kisspeptin expression in AVPV neurons (Smith et al., 2005; Gottsch et al., 2006; Oakley et al., 2009). Action potential firing by AVPV kisspeptin neurons is also increased on proestrus versus diestrus (a time of homeostatic negative feedback; Wang et al., 2016). Combined, these effects are thought to increase kisspeptin secretion onto the GnRH neuron, as increases in firing frequency are typically associated with increased neuropeptide release (Dutton and Dyball, 1979; Cropper et al., 2018). Estradiol positive feedback is likely conveyed to GnRH neurons by increased kisspeptin signaling, causing a surge of GnRH release which stimulates a luteinizing hormone (LH) surge from the pituitary. The LH surge stimulates ovulation (Greep et al., 1942).

If the above theory is correct, the firing activity of AVPV kisspeptin neurons is a critical piece of this surge induction process, serving as a gating mechanism controlling kisspeptin release and thus the physiologic cascade culminating in ovulation. Several studies have focused on the electrical properties of these cells and how they change across the reproductive cycle, measured via recordings performed in acutely prepared mouse brain slices during periods of estradiol negative (typically diestrus) versus positive (proestrus) feedback. AVPV kisspeptin neurons recorded during proestrus maintain higher firing rates than during diestrus even when fast synaptic transmission was blocked, suggesting a shift in intrinsic properties across the diestrus-proestrus transition enables increased firing during proestrus (Wang et al., 2016). Voltage-gated conductances may provide a basis for these shifts in activity since they help control membrane potential and shape firing dynamics in electrically-excitable cells. T-type calcium, persistent sodium, and hyperpolarization-activated currents were all found to be larger in AVPV kisspeptin neurons during proestrus (Piet et al., 2013; Zhang et al., 2015; Wang et al., 2016). Increases in these depolarizing currents, which can be active at subthreshold membrane potentials, could promote greater firing activity. Potassium currents, which strongly influence firing onset of many neurons and are vital for proper repolarization of the membrane (Hille, 2001), remained uncharacterized in AVPV kisspeptin neurons. K+ currents change across the estrous cycle and in response to estradiol in GnRH neurons and kisspeptin neurons in the arcuate nucleus (DeFazio et al., 2002, 2019; Pielecka-Fortuna et al., 2011).

We hypothesized voltage-gated potassium currents in AVPV kisspeptin neurons are regulated by estrous cycle stage. Classic voltage-clamp approaches revealed three conductances contribute to the total voltage-gated potassium current in these cells, and each component was characterized during the afternoon of diestrus (negative feedback) and proestrus (positive feedback). It can be difficult to predict how multiple voltage-gated conductances interact to control the membrane voltage. Computational modeling studies were thus conducted to better understand how estrous cycle modulation of these and other currents influence action potential firing output.

Materials and Methods

Animals

The University of Michigan Institutional Animal Care and Use Committee approved all procedures. Adult female Kiss1-hrGFP mice (age 55–154 d) (Cravo et al., 2011) which express humanized Renilla GFP under control of the kisspeptin promoter, were used for these studies. Mice were provided with Harlan 2916 chow and water ad libitum and were held on a 14/10 h light/dark cycle with lights on at 3:00 A.M. Eastern Standard Time. Estrous cycle stage was monitored by vaginal cytology for at least one week before experiments. Uterine mass was measured after brain slice preparation to confirm cycle stage. Uterine mass >100 mg indicated in vivo exposure to a high concentration of estradiol, typical of proestrus, whereas mass <60 mg indicated exposure to low estradiol, typical of diestrus (Shim et al., 2000).

Experimental design

Brain slices were prepared from cycling adult female mice during cycle stages corresponding to estradiol negative feedback (afternoon of diestrus) or positive feedback (afternoon of proestrus). Whole-cell voltage-clamp recordings of hrGFP-identified AVPV kisspeptin neurons were used to characterize macroscopic voltage-gated potassium currents; three currents (fast transient, slow transient, residual) were identified and separated using pharmacological and/or voltage-based subtraction methods for characterization. Current-clamp recordings were used to monitor action potential firing. Computational modeling was used to predict how each type of potassium current regulates AVPV neuron excitability, and how changes in multiple voltage-gated currents across the estrous cycle may impact excitability.

Brain slice preparation

Chemicals were purchased from Sigma-Aldrich unless noted. All solutions were bubbled with 95% O2/5% CO2 for at least 15 min before exposure to tissue. All mice were euthanized at 3–4 P.M. Eastern Standard Time, and the brain was rapidly removed and placed in ice-cold sucrose saline solution containing the following: 250 mm sucrose, 3.5 mm KCl, 25 mm NaHCO3, 10 mm D-glucose, 1.25 mm Na2HPO4, 1.2 mm MgSO4, and 3.8 mm MgCl2, at pH 7.6 and 345 mOsm. Coronal (300 μm) slices were cut with a VT1200S Microtome (Leica Biosystems). Slices were incubated in a 1:1 mixture of sucrose saline and artificial CSF (ACSF) containing the following 135 mm NaCl, 3.5 mm KCl, 26 mm NaHCO3, 10 mm D-glucose, 1.25 mm Na2HPO4, 1.2 mm MgSO4, and 2.5 mm CaCl2 at pH 7.4 and 305 mOsm for 30 min at room temperature (∼21–23°C) and then were transferred to 100% ACSF for an additional 30–180 min at room temperature before recording. For recording, slices were placed into a chamber and perfused (3 ml/min) with carboxygenated ACSF kept at 31°C with an inline heating unit (Warner Instruments). GFP-positive AVPV kisspeptin neurons were identified by brief illumination at 488 nm on an Olympus BX51WI microscope. Recordings were performed 1–4 h after brain slice preparation. No more than four cells were recorded per mouse; data values from cells from the same animal were not clustered in a manner that would typically contribute to reduced variability.

Voltage-clamp recordings

Recording micropipettes were pulled from borosilicate capillary glass using a Flaming/Brown P-97 puller (Sutter Instruments) to obtain pipettes with a resistance of 1.5–3.5 MΩ when filled with pipette solution, which consisted of 135 mm K-gluconate, 10 mm KCl, 10 mm HEPES, 5 mm EGTA, 0.1 mm CaCl2, 4 mm MgATP, and 0.4 mm NaGTP, 305 mOsm and pH 7.2 with NaOH. Pipettes were wrapped in Parafilm to reduce capacitive transients. All potentials reported were corrected on-line for a liquid junction potential of −15.7 mV (Barry, 1994). Recordings were performed with one channel of an EPC-10 dual patch-clamp amplifier and PatchMaster software (HEKA Elektronik). After achieving the on-cell configuration with seal resistance >2.0 GΩ, fast capacitive transients were minimized and the whole-cell configuration was achieved by rupturing the cell membrane with brief suction. The membrane potential was held at −70 mV between voltage-clamp protocols. Passive membrane properties and voltage-clamp quality were calculated from the averaged current response (after on-cell capacitive current subtraction) to sixteen −5 mV, 20-ms test pulses from a holding potential of −70 mV performed immediately before and after each protocol. To ensure adequate recording quality, the following criteria were required for inclusion for analysis: uncompensated series resistance (Rs) <20 MΩ, input resistance (Rinput) >500 MΩ, capacitance (Cm) between 8 and 30 pF, holding current (Ihold) between −60 and 10 pA (Table 1). Rs was compensated 50–85% for all protocols and recordings were excluded from analysis if >20% change in Rs occurred during the experiment.

View this table:
  • View inline
  • View popup
Table 1

Passive properties and statistical comparisons

To isolate voltage-gated potassium currents, TTX (2.5 μm, Tocris Bioscience), CdCl2 (100 μm), and NiCl2 (300 μm) were included in the ACSF to block Na+ and Ca2+ currents. Picrotoxin (100 μm) was used to block GABAA-receptor-mediated currents and D-(-)-2-amino-5-phosphonopentanoic acid (APV; 20 μMm, Tocris Bioscience), and 6-cyano-7-nitroquinoxaline-2,3-dione (CNQX; 10 μm) were used to block ionotropic glutamate receptor-mediated currents.

Total K+ current

In initial recordings of the total voltage-gated potassium current, two temporally-distinct peaks were visible in voltage steps from −100 to −10 or 0 mV. Further, a persistent negative slope in the outward current suggested a component with slow inactivation was present at potentials depolarized to −10 mV. Lengthening the voltage-clamp protocols allowed for more complete removal or induction of inactivation before the test pulse. To measure the voltage dependence of activation of the total voltage-gated K+ current, the membrane potential was held at −100 mV for 5 s to remove inactivation then stepped to test pulses of −70 to +40 mV (10-mV increments) for 1 s. To measure the voltage dependence of inactivation, the membrane was held at potentials ranging from −100 to −10 mV (10-mV increments) for 10 s and then stepped to a test pulse of 0 mV for 1 s. The long duration of these protocols makes P/N subtraction (Armstrong and Bezanilla, 1977) unfeasible. Instead, slow Cm and leak current subtraction were performed offline by subtracting the appropriately scaled average current response to one hundred −5 mV test pulses with Cm and Rs compensation activated (Kimm et al., 2015). Averages were scaled for each voltage step used in the protocol.

Slow-transient K+ current

To characterize the slow transient K+ current, the same activation/inactivation protocols as above were repeated with 5 mm 4-aminopyridine (4-AP) included in the ACSF to reduce primarily the fast-transient current. To measure the time course of inactivation of the slow-transient K+ current, the membrane was held at −100 mV for 5 s to remove inactivation, then at −10 mV for durations varying from 0 to 8 s, followed by a test pulse at 0 mV for 1 s. To measure the time dependence of recovery from inactivation, the membrane was held at −10 mV for 10s to inactivate the current, then at −100 mV for durations varying from 0 to 8 s, followed by a test pulse at 0 mV for 1 s. Offline leak subtraction was used as for recordings of the total K+ current (Kimm et al., 2015).

Fast-transient K+ current

The fast-transient K+ current was isolated and quantified using a voltage-based subtraction method after reducing primarily the slow-inactivating and residual components with ACSF in which 20 mm tetraethylammonium chloride (TEA) replaced 20 mm NaCl. To measure the voltage dependence of activation, two protocols (A and B) were run in series. For protocol A, the membrane potential was held at −100 mV for 200 ms, then current was measured during 150-ms test pulses from −70 to +40 mV (10-mV increments). Protocol B was identical to protocol A, except that the membrane potential was held at −30 mV rather than −100 mV during the first 200 ms of the protocol. This fully inactivates the fast-transient current, isolating any remaining residual component during the test pulse. The current response during the test pulse of protocol B was subtracted from that of protocol A to yield the fast-transient current. To measure the voltage dependence of inactivation, the membrane was held at −100 to −20 mV (10-mV increments) for 200 ms, before stepping to a test pulse at −10 mV (150 ms). Current occurring after a step from −30-mV holding potential was defined as residual current and was subtracted from earlier sweeps to yield isolated fast transient current. To measure the rate of inactivation of the fast-transient K+ current, inactivation was fully removed by holding the membrane at −100 mV for 300 ms. The membrane was then stepped to −30 mV for 0–250 ms followed by a test pulse of 0 mV for 100 ms. The component of the current that did not inactivate (residual current) was subtracted from the prior series to isolate the fast-transient current. To measure the rate of the recovery from inactivation, the fast-transient current was fully inactivated by holding the membrane at 0 mV for 300 ms, then recovery prepulses at −100 mV were applied for 0–250 ms before measuring current at 0 mV (100 ms). The current after the 0-ms recovery prepulse was subtracted from each test pulse to isolate the fast-transient K+ current. Slow Cm and leak currents were corrected using P/−5 online leak subtraction with Vhold = −70 mV (Armstrong and Bezanilla, 1977).

It was not possible to quantify all parameters from every cell because not all cells remained within our quality control standards long enough for inclusion in the analysis of the time dependence of activation and inactivation. The number of cells and animals per measurement are shown in Tables 2, 3.

View this table:
  • View inline
  • View popup
Table 2

Two-way ANOVA analyses of K+ current properties

View this table:
  • View inline
  • View popup
Table 3

Two-sample analyses of K+ current properties

Total, fast-transient, slow-transient and residual K+ current densities were calculated by dividing by Cm. Peak values were divided by the driving force according to the Goldman–Hodgkin–Katz (GHK) current equation using the calculated reversal potential (−94 mV; Clay, 2000, 2009). Conductance values were then normalized to the maximum conductance (gmax) to generate steady state activation/inactivation curves, which were fit with a Boltzmann function to calculate voltage of half (V50) activation or inactivation, and steepness (k). Inactivation and recovery from inactivation curves were fit to single exponentials to measure time constants (τ).

Current-clamp recordings

Depolarizing (2–28 pA, 2-pA increments) and hyperpolarizing (−25 to −5 pA, 5-pA increments) current pulses (500 ms) were applied in whole-cell configuration. The maximum amplitude of hyperpolarizing stimuli was limited to avoid hyperpolarizing the membrane beyond −100 mV, as further hyperpolarization was detrimental to recording stability. The targeted initial membrane potential was within 2 mV of −70 mV, which is close to the baseline membrane potential of these cells (DeFazio et al., 2014). Spikes were counted as action potentials only if their peak value was depolarized relative to −10 mV; spike number was plotted as a function of current injection.

Modeling

For physiologic data used, mean ± SEM normalized activation, inactivation, recovery, and rate of inactivation curves were calculated for the total, fast-transient (fast), slow-transient (slow), and residual (resid) K+ current densities in each group. Leak-subtracted current traces from activation protocols were averaged for each current at each voltage step to obtain mean ± SEM current traces used for fitting. T-type calcium and persistent sodium conductance data were collected from published measurements, which indicated larger amplitudes during proestrus (Zhang et al., 2015; Wang et al., 2016). The hyperpolarization-activated conductance was estimated based on the amplitude and time course of the sag potential, which has larger amplitude during proestrus (Piet et al., 2013; Wang et al., 2016). Neurons were modeled as a single compartment with whole-cell Cm of 17.62 pF based on the mean of experimentally-observed values (Table 1).

Equations

Na+, Ca2+, and hyperpolarization-activated currents were modeled based on Ohmic driving forces (Eq. 1–2). K+ currents were modeled based on nonlinear driving forces described by the GHK current equation (Eq. 3; Clay, 2000, 2009). I=CmdVdt + INaT + INaP + ICaT + Ih + Iresid + Ifast + Islow + ILeak (1) Ia=g¯amrh(V−Ea),a=Nap,CaT,h,Leak (2) Ib=g¯bmhqVkTeq(V−Eb)kT−1eqVkT−1,b=resid,fast,slow. (3)

I is current (nA) across the neuronal membrane, Cm is the cell Cm (pF), Vm is the membrane potential (mV), t is time (ms), g¯a is the specific conductance (nS), Ea is the reversal potential (mV) of a given ion, q is the elementary charge, k is the Boltzmann constant, T is the absolute temperature, and m and h represent activation and inactivation variables, respectively. Exponents r represents the number of activation particles. These variables are governed by the following differential equations (Eqs. 4, 5) dmdt=m∞(V)−mτm(V) (4) dhdt=h∞(V)−hτh(V), (5)where m∞(V) and h∞(V) are the steady-state activation/inactivation functions for each variable and τm and τh are functions determining the voltage-dependent time constants (in ms) of activation/deactivation and inactivation/recovery, respectively. These functions take the form (Eq 6–9). A-G are free parameters adjusted during fitting. m∞,h∞=11 + e(V50m,h−V)Km,h (6) τma,τha=A + B1 + e(V−CD)*E1 + e(V−FG),a=slow,resid (7) τma,τha=A + B1 + e(V−CD),a=fast,H,CaT (8) τNaP=EeA + VB + eC+VD + F. (9)

The transient sodium current (NaT) was modeled used a Markov chain with three states, open (O), closed (C), or inactivated (I). Transitions between states were either voltage sensitive or constants. The model used was modified from Adams et al., 2018. (Eq. 10–14). INaT=g¯NaTO3(V−ENa) (10) α(V),β(V),r1(V),r3(V)=r1 + e(V+s)k (11) dCdt=β(V)O + r3(V)I− α(V)C (12) dOdt=α(V)C−(β(V) + r1(V))O (13) I=1−C−O. (14)

Parameter optimization

Simulated voltage-clamp experiments were performed using xolotl (Gorur-Shandilya et al., 2018) in MATLAB (MathWorks) running on a Windows 10 PC. Voltage-step protocols matched those used experimentally for each conductance. Current traces generated by models were analyzed in the same manner as the slice recording data. For each voltage step, peak currents were converted to conductance and normalized to generate activation/inactivation curves. For K+ models, permeability as calculated by the GHK current equation was used (Clay, 2000, 2009). m∞ and h∞ function parameters were adjusted manually until the normalized activation/inactivation curves were well fit to the experimental data. τm, τh, and g¯ parameters were also adjusted so that the rise and decay (where applicable) phases of the current trace and peak amplitudes were well-fit to the data. This process was iterative and was greatly facilitated by xolotl’s “puppeteer” function, which allowed for rapid-recalculation and visualization of simulation results following parameter adjustment. Separate models were generated for each experimental group to reproduce estrous cycle effects on each of the voltage-gated currents.

The total K+ current was fit by simulating voltage-clamp of a cell containing the three K+ conductance subcomponents (fast, slow, resid) as they were fit to their respective drug-separated recordings. To better fit the total K+ recording, g¯ parameters were optimized using particle swarm optimization (Global Optimization Toolbox, MathWorks), which sought to minimize the sum of squares error between the current trace generated by the model and the mean data for the respective conductance.

Simulation of excitability

To study action potential firing, diestrous and proestrous neuron models were built with an identical sodium channel model. This model was made the same between groups because, although there are no experimental data characterizing this conductance in AVPV kisspeptin neurons, action potential threshold (defined as when rate of rise exceeds 2 V/s) and amplitude were the same between cycle stages (Wang et al., 2016), and these characteristics are highly influenced by transient sodium current properties (Kress et al., 2010).

To examine how conductances alter firing, a suite of conductance models corresponding to the respective experimental group were added to single-compartment neuron models, and the voltage of the model was allowed to vary, i.e., “current clamp.” Using the same protocols from in vitro recordings, firing was stimulated by applying increasing current steps (−30–0 pA at 10-pA intervals, 2–28 pA at 2-pA intervals) to the model cell at a baseline membrane potential of −70 mV. The number of spikes fired during the stimulation is plotted against the stimulation amplitude to generate a firing-current (F-I) curve. To achieve this baseline potential, Ihold and leak conductance were adjusted so that the membrane potential before the start of the current pulse was within 50 μV of −70 mV and Rinput was 1 GΩ (passive properties in Table 1). For diestrous and proestrous models, the simulated rheobase action potential waveform was compared with the peak-aligned mean rheobase action potential of the corresponding experimental group. This fit was improved by manually adjusting sodium conductance model parameters; adjustments were kept the same for the diestrous and proestrous models. To test predictions regarding estrous cycle effects or the role and/or interaction of various conductances, individual or multiple conductances fit to diestrous data were substituted into the proestrous model (and vice versa) for a particular simulation.

Code accessibility

The code described in the paper is freely available online at https://github.com/JRS92/Starrett_eNeuro_2021. The C++ files containing xolotl conductance code is available as Extended Data 1.

Extended Data 1

CaTcustom: hpp file containing code for T type calcium current; HCurrentcustom: hpp file containing code for hyperpolarization-activated current; markovNA: hpp file containing code for transient Na current; NaP_V2: hpp file containing code for persistent Na current; optifastK_GHK: hpp file containing code for fast K current; optiSlowKV5: hpp file containing code for slow K current; testcond: hpp file containing code for residual K current. These conductance files are written in C++ and can be added to xolotl neuron objects. Parameter values used for models are listed in Tables 4-Tables 6. Download Extended Data 1, ZIP file.

Statistics

Data were analyzed using Prism 9 (GraphPad). Shapiro-Wilk was used to test distribution of data. Parameters (mean ± SEM) were compared between cycle stages with unpaired statistical tests appropriate for experimental design and data distribution as indicated in the results section. The number of cells per group is indicated by n. For two-way ANOVAs and Mann–Whitney, difference of the means (Diff) was defined for cycle stage (diestrus–proestrus). For Kruskal–Wallis tests (Table 1), mean rank differences are reported as (total – slow), (total – fast), (slow – fast). Significance was set at p < 0.05.

Results

Passive properties

Rs, Rinput, Cm, and Ihold measurements were made before and after each protocol and values were averaged for each cell. There were no statistical differences between cells recorded during diestrus versus proestrus for any of the measured properties (Table 1). Comparing across drug treatments that are discussed below, Rinput was higher during 4-AP or TEA treatment (Table 1). Dunn’s multiple comparison post hoc tests of Rinput revealed increases consistent with antagonism of channels active at baseline membrane potentials by the respective drugs (Dunn’s, control vs 4-AP p = 0.002, control vs TEA p = 0.001, 4-AP vs TEA p >0.99).

Total voltage-dependent K+ current in AVPV kisspeptin neurons has multiple components

In nocturnal rodents, positive feedback needed to generate the LH surge occurs during the afternoon of proestrus (Brown-Grant et al., 1970; Sarkar et al., 1976). During this time, both overall firing rate and burst firing of AVPV kisspeptin neurons increase relative to diestrus (Wang et al., 2016). Given the importance of K+ currents in determining the excitability and firing patterns of neurons (Coetzee et al., 1999; Hille, 2001), we hypothesized voltage-gated potassium currents are targets of estradiol feedback in AVPV kisspeptin neurons; we tested this using whole-cell voltage-clamp. Figure 1A shows a representative recording of total voltage-dependent K+ current activation. Examining this family of traces from the most hyperpolarized to the most depolarized command potentials revealed induction of one component with faster inactivation at more hyperpolarized potentials that merged with a component that activated at more depolarized potentials and only slowly inactivated. The slow kinetics of the macroscopic current resembled a K+ current previously observed in dorsal root ganglion neurons (Gold et al., 1996), and led us to extend the duration of the prepulses beyond what is used in a typical inactivation protocol to 10 s (Fig. 1B). As described above, a fast component can be seen to activate and inactivate within ∼50 ms at more hyperpolarized potentials during the prepulse (expanded in Fig. 1C, left); this component is prominent as a spike during the test pulse following hyperpolarized prepulses (expanded in Fig. 1C, right). A second component required greater depolarization to begin activating and inactivated slowly only at more depolarized potentials. The peak total current density (derived from the voltage-clamp protocol in Fig. 1A) was not different between diestrus and proestrus at any of the potentials tested (Fig. 1D, two-way repeated-measures ANOVA for effect of cycle stage; Table 2). These observations, however, do not preclude latent changes attributable to individual components being differentially regulated.

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Total voltage-dependent K+ current in AVPV kisspeptin neurons has three distinct components. A, Representative total K+ current in response to the voltage-clamp protocol shown. B, Representative recording of total K+ current from a different cell in response to the voltage protocol shown; both a slowly-inactivating and residual (sustained) component are evident during the extended prepulse. C, Expansions of the areas within the dashed boxes in B, showing the change in inactivation rates as more depolarizing prepulses are applied (left), and how this affects both activation and inactivation during the test pulse (right). D, Mean ± SEM peak current density in cells from mice in diestrus (black symbols) and proestrus (magenta symbols). Error bars are smaller than symbols for some values.

Pharmacological treatments to isolate current components

To attempt to isolate the components of the total potassium current, recordings were made before and during application of the potassium channel blockers TEA (20 mM) or 4-AP (5 mM). TEA reduced the amplitude of both the fast and slow components, but the reduction in the slow component was proportionally greater, making the fast peak appear more prominent (Fig. 2A). Lower concentrations of TEA (500 μm to 10 mm) did not block most of the slow peak amplitude and were thus less effective in revealing the fast component. The fast-transient component was markedly reduced by 5 mm 4-AP (Fig. 2B, note right-shift of peak in red trace). Lower concentrations of 4-AP (100 μm to 1 mm) did not block the fast initial peak at −10 mV. These changes were attributable to drug action rather than rundown as the current remained stable in the absence of treatments over the 10 min required for these recordings (Fig. 2C).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Potassium current recorded from three different cells before and during application of potassium channel blockers. A, 20 mm TEA. B, 5 mm 4-AP. C, Time control.

Individual components of the total current may have distinct roles in determining firing output and may also be independently regulated by cycle stage. The differential effects of TEA and 4-AP were used to perform more detailed characterization of the sub-components of the total potassium current in AVPV neurons in slices from diestrous and proestrous mice. To increase recording yield, we switched from a within-cell design to a population (i.e., among-cell) design, with drug present from the time the slice was placed in the recording chamber. In these population studies, the longer exposure to TEA caused a larger reduction in peak current density. Specifically, the peak current density at +40 mV during TEA treatment was 76.6 ± 5.1% of pretreatment values in within-cell experiments, but only 53.0 ± 3.7% of untreated cells in population studies. This greater reduction was not observed with 4-AP [the peak current density during 4-AP treatment was 76.6 ± 4.9% of the control period for within-cell studies (n = 3), and 76.1 ± 9.9% of untreated controls in population studies (n = 4)].

The slow-transient component is increased on diestrus

The extended duration voltage-clamp protocols were used to characterize the slow-transient potassium current in the presence of 4-AP; representative traces are shown in Figure 3. Peak current density was greater in cells from diestrous (n = 12) than proestrous (n = 13) mice at more depolarized potentials (Fig. 3C, two-way repeated-measures ANOVA for effect of cycle stage; Table 2). The slow current was activated at command voltages at and depolarized relative to −20 mV. Normalized activation and inactivation curves were not different between groups (Fig. 3D, two-tailed Student’s t test with Welch’s correction for V50 activation, inactivation; Table 3). Time course of inactivation was not different between diestrus and proestrus (Fig. 3F,H, two-way repeated-measures ANOVA for effect of cycle stage; Table 2). In contrast, recovery from inactivation was more rapid during proestrus (Fig. 3G,H, two-way repeated-measures ANOVA for effect of cycle stage; Table 2). In recovery from inactivation experiments, the recovery of two currents was evident. The fast transient was visible when the hyperpolarizing recovery prepulse was shorter than 200 ms (Fig. 3G, arrow). This may be attributable to the voltage dependence of 4-AP action, in which the channel becomes unblocked during depolarization (Choquet and Korn, 1992; Kehl, 2017). The slow current was visible 70–150 ms into the test pulse. To avoid contamination with the fast current, for this protocol only the peak of the slow transient current was taken as the maximum current 70–150 ms after the start of the test pulse.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

4-AP-resistant slow-transient voltage-dependent K+ current is larger on diestrus. A, Representative traces (top) in response to the activation voltage-clamp protocol (bottom) in the presence of 5 mm 4-AP. B, Representative traces (top) in response to the inactivation voltage-clamp protocol (bottom). Offline leak subtraction was applied to test pulses only, hence the Cm transient is visible at the start of the recording. C, Mean ± SEM peak current density from cells in mice in diestrus (black symbols) and proestrus (magenta symbols). D, Mean ± SEM normalized conductance. Solid lines are Boltzmann fits to the mean data. E, Individual values and mean ± SEM parameters obtained from Boltzmann fits to normalized conductance curves for each cell. F, G, Representative traces (top) in response to the voltage-clamp protocol (bottom) used to measure the time dependence of inactivation (F) and recovery (G). Arrow denotes peak of fast transient current. H, Mean ± SEM normalized peak current versus prepulse duration for time dependence of inactivation (inact) and recovery from inactivation. Error bars are smaller than symbols for some values; *p < 0.05.

The inactivation curve of the fast-transient component is depolarized during proestrus

Bath-applied 20 mm TEA preferentially reduced the amplitude of the slow-transient current, allowing better characterization of the fast-transient current. Because this fast-transient component activated and inactivated quickly, voltage-clamp protocols with more traditional durations were used (Fig. 4A). Under these recording conditions, a fast-transient current began to activate between −50 and −40 mV (Fig. 4C); current reached a non-zero steady state by the end of the test pulse, suggesting a contribution by non-inactivating or slowly-inactivating channels. To isolate the fast-transient current from this non-inactivating residual component, a voltage-based subtraction method was used. Subtraction of the residual current following a depolarizing prepulse (Fig. 4A, middle) yielded a transient current that was almost fully inactivated by the end of the test-pulse (Fig. 4A, right). The V50 activation of the residual current was more depolarized than that of the fast-transient current at both cycle stages (Fig. 4B–D, diestrus n = 11, proestrus n = 13, one-way ANOVA F(3,44) =10.8, p < 0.001; Tukey residual vs fast p = 0.001 for proestrus, residual vs fast p = 0.002 for diestrus). There was no effect of cycle stage on the peak current density of either the fast-transient or residual component (Fig. 4B, two-way repeated-measures ANOVA for effect of cycle stage; Table 2). In contrast, cycle stage affected voltage dependence of the fast-transient component. Specifically, inactivation of the fast transient occurred at more depolarized potentials on proestrus than diestrus (Fig. 4C,E, two-way repeated-measures ANOVA for effect of cycle stage on inactivation; Table 2; unpaired two-tailed Student’s t test with Welch’s correction of V50 inactivation; Table 3), Cycle stage did not, however, affect the voltage dependence of activation (Fig. 4C, two-way repeated-measures ANOVA for effect of cycle stage on activation curve; Table 2) or V50 activation (Fig. 4D, unpaired two-tailed Student’s t test; Table 3). Full inactivation and full recovery of the fast transient occurred in <100 ms (Fig. 4F,G). There was no group difference in either of these properties (two-way repeated-measures ANOVAs for effect of cycle stage; Table 2).

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

TEA-resistant fast-transient voltage-dependent K+ current has a depolarized inactivation curve during proestrus. A, left, Representative unsubtracted K+ current (top) in response to the voltage-clamp protocol (bottom) in presence of 20 mm TEA. Middle, Residual current in the same cell with a −30-mV prepulse. Right, Fast transient current yielded by subtracting residual current from raw current. B, Mean ± SEM peak current density in cells from mice in diestrus (black symbols) and proestrus (magenta symbols). C, Mean ± SEM normalized conductance. Solid lines represent Boltzmann sigmoidal fits to the mean data for the fast transient, dashed lines are the fits for the residual current. D, Individual values and mean ± SEM parameters obtained from Boltzmann fits to normalized conductance curves for each cell. E, F, Representative traces (top) in response to the voltage-clamp protocol (bottom) used to measure the time dependence of inactivation (E) and recovery (F). G, Mean ± SEM normalized peak current versus prepulse duration for time dependence of inactivation (inact) and recovery from inactivation. Error bars are smaller than symbols for some values; *p < 0.05.

Modeling conductances

Computational modeling was used to better understand the role of estrous cycle modulation of multiple conductance types in the control of AVPV kisspeptin neuron firing. Specifically, we examined whether estrous cycle-related shifts in K+ current components and/or subthreshold depolarizing currents alter the firing dynamics of these cells. To model subcomponents in silico, Hodgkin–Huxley conductance models were fit to mean ± SEM recordings of drug-separated and voltage-separated subcomponents of the K+ current: slow, fast, and residual. This involved adding individual subcomponents to single-compartment neuron models and running such models through the same voltage-clamp protocols used experimentally, generating membrane current over time plots. Model output was compared with experimental data (Fig. 5), and m∞(V), h∞(V), τm(V), and τh(V) function parameters were changed iteratively to improve the quality of fit. Steady-state activation and inactivation curves and I(density)-V plots were calculated from the model output and compared with experimental data for K+ currents (Fig. 6A,B,D) and for CaT (Fig. 6C,F) and NaP (Fig. 6E) currents reported (Zhang et al., 2015; Wang et al., 2016).

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

Potassium conductance model output versus experimental data for voltage steps. Rainbow colors indicate mean ± SEM current traces recorded at different test potentials; black lines show model simulations. In each panel, diestrus is on the left and proestrus on the right. A, Slow (voltage protocol as in Fig. 3A). B, Fast (voltage protocol as in Fig. 4A). C, Residual (voltage protocol as in Fig. 4A).

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

Conductance models versus data for steady state activation/inactivation and current density experimental data. A–C, Steady state activation/inactivation curves calculated from voltage-clamp simulations (dashed lines) compared with corresponding mean ± SEM experimental data (symbols) from diestrous (black) and proestrous (magenta) groups. D–F, Peak current density for various conductances. Experimental data for K+ current subcomponents are the same as in Figures 3, 4 and replotted here for ease of comparison. Data points for NaP (E) and CaT (C, F) are adapted from Wang et al. (2016).

To reconstruct the total K+ current from the sum of the subcomponents, the set of Islow, Ifast, and Iresid currents for each group were included in single neuron models, and the models were run through the same voltage-clamp protocol used for recordings in Figure 1A. Initial simulations indicated the sum of fast, slow, and residual conductance models as fit to their respective isolated recordings was not enough to account for the full amplitude of the total K+ current, suggesting subcomponents were blocked to some degree by the drug used to isolate them (i.e., 4-AP blocked a portion of the slow current, TEA blocked some of the fast and residual currents). To overcome this problem, particle-swarm optimization was used to minimize the error between the model-generated current trace and the mean ± SEM trace by optimizing g¯ parameters of each subcomponent, with equal weighting across the trace. This preserved the voltage-sensitivity and kinetics of the subcomponents as fit in the drug-separated and voltage-separated recordings (Fig. 5) while providing a better fit to the experimental data. Parameter values are shown in Tables 4, 5. This approach provided a good fit to peak current density (Fig. 7A) and the first ∼250 ms of the total K+ current (Fig. 7B). The error between the model and data increased as the simulation progressed, however. This error could be reduced by decreasing the inactivation rate of Islow, but this modification had no effect on firing behavior as measured in response to 500-ms current pulses, suggesting this shift in kinetics does not substantially influence firing output at these timescales (data not shown).

View this table:
  • View inline
  • View popup
Table 4

Model parameters for diestrus

View this table:
  • View inline
  • View popup
Table 5

Model parameters for proestrus

Figure 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 7.

Reconstruction of the total K+ current from the sum of the three subcomponents. A, Peak current density when Ifast, Islow, and Iresid are simulated together in the same model (dashed lines) after g¯ optimization to correct for suppression by TEA/4-AP. Mean ± SEM symbols (black: diestrus, magenta: proestrus) and voltage-clamp protocols are the same as shown in Figure 1A,D and are reproduced for ease of comparison. B, Mean ± SEM experimental current traces at different test pulses (rainbow colors) and model simulation (black).

Simulations of firing

K+ conductance models from fits to the total K+ current, CaT, NaP, and HCN conductances for the diestrous and proestrous groups were combined into single-compartment neuron models, thus creating neuron models with seven voltage-gated conductances and one leak conductance each. A transient sodium (NaT) conductance is needed to enable action-potential firing and was thus added to the models at this stage (Table 6). Transient sodium currents of AVPV kisspeptin neurons have not been measured. We thus estimated NaT parameters using action potential properties such as threshold and peak amplitude. NaT was kept identical between diestrus and proestrus as these action potential properties did not change with cycle stage previously (Wang et al., 2016) or in this study (threshold: diestrus −51.8 ± 1.0 mV, n = 12 cells from 5 animals, proestrus −51.6 ± 0.6 mV, n = 17 cells from 6 animals, two-tailed, unpaired Student’s t test, p = 0.85, t = 0.1906, df = 27, amplitude diestrus 82.1 ± 1.9 mV, n = 12, proestrus 80.5 ± 1.1 mV, n = 17, two-tailed unpaired Student’s t test, p = 0.4501, t = 0.72, df = 22), suggesting estrous cycle modulation of this conductance is minimal. Using the same NaT conductance in the diestrous and proestrous models generated action potential shapes similar to experimental action potentials (Fig. 8C). The diestrous and proestrous models were tested with the same current injection protocols used in recordings and demonstrated F-I curves similar to the mean ± SEM data, indicating that the models could reproduce the firing behavior of AVPV kisspeptin neurons recorded in brain slices (Fig. 8D).

View this table:
  • View inline
  • View popup
Table 6

NaT Parameters

Figure 8.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 8.

Simulations of firing from a baseline of −70 mV. A, Performance of diestrous (black) and proestrous (magenta) models in response to −30 pA (top), +6 pA (middle), and +28 pA (bottom) applied current. Boxed regions are shown with expanded axes in the indicated subfigure. B, Posthyperpolarization rebound of diestrous and proestrous models from panel A. C, Rheobase action potentials for both models (dashed lines) from A (middle) and experimental counterparts (solid lines, mean ± SEM). D, F-I curves for diestrous and proestrous models (lines) and experimental counterparts (circles are means, shading is SEM). E, F-I performance of hybrid models in which one or multiple K+ conductances in the diestrous model was substituted for a proestrous counterpart. Non-hybrid models from part D (di full and pro full) are reproduced here to facilitate comparison in this panel as well as F, G. F, F-I performance of hybrid models in which one or multiple subthreshold depolarizing currents in the diestrous model were replaced with proestrous counterparts. G, Rebound bursting performance of hybrid models in response to −30-pA current injection.

These models also faithfully reproduced cycle-dependent differences in posthyperpolarization rebound firing behavior observed (Wang et al., 2016). Specifically, cessation of hyperpolarization initiated rebound spikes in the proestrous but not the diestrous model (Fig. 8B). To determine the basis for this difference as well as to further understand how each conductance influences firing behavior, we performed a series of simulations involving swapping individual or multiple conductances from one model into the other. In the diestrous model, replacing all K+ conductances with their proestrus counterparts had a suppressive effect, depolarizing the rheobase and reducing the number of spikes fired in response to all depolarizing stimuli (Fig. 8E). Of the three potassium current subcomponents, proestrous Islow was the most suppressive when individually substituted into the diestrous model (Fig. 8E). The diestrous model with only Ifast substituted for the proestrous counterpart had a depolarized rheobase but the F-I curve (Fig. 8E) had a similar slope. Changing from diestrous to proestrus Iresidual did not change firing behavior (Fig. 8E), consistent with this current being similar between groups (Fig. 4). Of the subthreshold depolarizing conductances, proestrous CaT caused the greatest increase in excitability when substituted into the diestrous model, generating more spikes at every depolarizing stimulus after rheobase was achieved, but not changing rheobase itself (Fig. 8F). Rebound firing only occurred in the diestrous model when all three subthreshold depolarizing currents were replaced with their proestrous counterparts (Fig. 8F,G). This model, which had diestrous K+ conductances and proestrus CaT, NaP, and h-current, had the greatest excitability of any model tested. It is important to point out that changes in excitability were not attributable to differences in either baseline membrane potential or Rinput as these were the same among models. Comparing this hybrid model to the proestrus-only model thus further indicates the suppressive effect of proestrous K+ currents.

Discussion

AVPV kisspeptin-expressing neurons are key mediators of estradiol positive feedback, a process critical for ovulation. We studied intrinsic properties of these neurons to determine how they become more active during proestrus, the phase of the rodent estrous cycle when positive feedback culminates in the ovulation-inducing GnRH/LH surges. AVPV kisspeptin neurons exhibit three types of voltage-gated potassium currents; two were modulated in an estrous-cycle dependent manner. Unexpectedly, effects of the proestrous K+ currents were suppressive but countered by cycle-dependent shifts in inward currents active in the subthreshold range. The proestrous shift in K+ currents did not, however, block increased rebound firing enabled by the proestrous shift in inward currents.

This finding that the K+ currents on proestrus were suppressive was initially surprising given substantial evidence indicating AVPV kisspeptin neurons are activated during estradiol positive feedback (Zhang et al., 2015; Wang et al., 2016, 2019). Much of this evidence was gathered using experimental models of estradiol feedback involving ovariectomy plus estradiol replacement (OVX+E) as a constant-release implant producing daily surges (Christian et al., 2005), or an implant followed several days later by an estradiol benzoate injection that triggers an LH surge the next day (Bronson and Vom Saal, 1979). These models allow study of estradiol feedback in isolation of other ovarian factors. In both models, AVPV kisspeptin neurons have increased c-fos expression during the LH surge (Adachi et al., 2007; Porteous and Herbison, 2019). Increases in spontaneous firing rate during positive feedback in the daily surge model are similar to those observed in naturally cycling animals on proestrus (Wang et al., 2016, 2019). Interestingly, given the changes in spontaneous firing, no studies have reported that feedback or estrous cycle stage affect the F-I curve of AVPV kisspeptin neurons. Our data and simulations suggest this lack of effect on the F-I curve is attributable to more suppressive K+ currents during proestrus counteracting larger inward currents during depolarization, without inhibiting rebound burst firing.

Estradiol feedback also affects voltage-gated K+ currents in other central neurons controlling reproduction. In GnRH neurons, inactivation of the transient K+ current is depolarized by positive feedback in the daily surge model to a similar degree observed in the present study (DeFazio et al., 2002). In contrast to AVPV neurons, however, peak K+ current densities in GnRH neurons are reduced during positive feedback, and excitability is increased (Adams et al., 2018). Interestingly, positive feedback effects in GnRH neurons may be at least in part attributable to increased kisspeptin receptor activation, as kisspeptin treatment has similar effects on transient K+ currents in these cells (Pielecka-Fortuna et al., 2011). In arcuate kisspeptin neurons, K+ currents were compared between open-loop (OVX) and estradiol negative feedback (DeFazio et al., 2019), with negative feedback reducing peak transient and sustained K+ currents compared with open-loop but not affecting voltage sensitivity of activation/inactivation. Together these observations suggest the effects of estradiol and cycle stage on voltage-gated K+ currents in the reproductive neuroendocrine circuitry depend on both cell type and type of feedback involved.

The subcellular mechanisms linking estradiol feedback to these effects have not been extensively studied but may include changes in channel subunit expression/composition, subcellular location of channels, binding by partner proteins, and/or subunit phosphorylation state (Coetzee et al., 1999; Levitan, 2006). In GnRH neurons, estradiol or cycle stage regulates expression of Ca2+, HCN, and SK channel genes (Zhang et al., 2009; Bosch et al., 2013; Rønnekleiv et al., 2015; Vastagh et al., 2019). K+ channel genes are also estradiol-sensitive in the paraventricular and arcuate nuclei (Qiu et al., 2006; Roepke et al., 2007; Lee et al., 2013; Yang et al., 2016). Estradiol affects transcription of ion channels outside the central nervous system in neural (Du et al., 2014) and non-neural tissue (Marni et al., 2009; Banciu et al., 2018). In addition to transcription, estradiol action through membrane-associated signaling pathways can rapidly modulate L-type and R-type Ca2+ and K-ATP currents (Sun et al., 2010; Zhang et al., 2010). Estradiol positive feedback effects on K+ currents in GnRH neurons are attenuated by a broad-spectrum kinase inhibitor, suggesting a role of phosphorylation (DeFazio and Moenter, 2002). Estradiol may also interact directly with channels, as with voltage-gated BK channels in vascular smooth muscle (Granados et al., 2019). Effects of estradiol on CaT currents in AVPV kisspeptin neurons, however, are dependent on actions via ERα as cre-lox or CRISPR mediated ERα knock-down in these cells eliminates estradiol effects (Wang et al., 2019).

The present results add important information to a growing list of voltage-gated conductances in AVPV kisspeptin neurons that are modified by the estrous cycle (Piet et al., 2013; Zhang et al., 2015; Wang et al., 2016). We used in silico approaches to amalgamate multiple experimental findings into a more comprehensive understanding of how hormonal manipulations/estrous cycle shape the membrane response of these cells as in other reproductive neuroendocrine neurons (Moran et al., 2016; Adams et al., 2018; Mendonça et al., 2018; DeFazio et al., 2019). Simulation results suggest modulation of K+ currents across the negative to positive feedback transition do not increase AVPV neuron excitability or enable rebound burst firing, as had been suggested by blocking the fast-transient current with 4-AP in current-clamp studies (Wang et al., 2016). Rather, proestrous shifts in subthreshold inward currents CaT, NaP, and HCN counteract suppressive K+ currents during proestrus and enable rebound firing. Involvement of the former currents in rebound firing is consistent with their role in shaping firing behavior (Lüthi and McCormick, 1998; Bevan and Wilson, 1999; Hille, 2001; Perez-Reyes, 2003).

The finding that cycle-dependent changes in intrinsic properties lead to no net change in excitability begs the question of how increased firing rates during proestrus are maintained when receptors for fast glutamate and GABA transmission are blocked (Wang et al., 2016). First, neuromodulation via metabotropic receptors could occur, and may alter non-voltage-sensitive mechanisms, such as G-protein-coupled inward rectifying potassium (GIRK) channels (Kelly and Qiu, 2010). Second, proestrous increases in spontaneous firing rate may be because of the emergence of rebound bursting (Wang et al., 2016). The suppressive K+ current changes observed on proestrus did not inhibit this process and may bestow properties to the cell that are otherwise beneficial. For example, these more suppressive K+ currents may serve a protective/homeostatic function by limiting the increase in excitability enabled by increases in depolarizing currents. An optogenetic study examining the relationship between AVPV neuron stimulation and GnRH neuron firing found delayed activation of GnRH neurons characteristic of kisspeptin receptor activation plateaus at 10 Hz (Piet et al., 2018), suggesting firing frequencies above this are unnecessary. Third, excitation-release coupling in AVPV neurons may be regulated by the estrous cycle. Small increases in release probability or the readily releasable pool of kisspeptin could have substantial downstream effects on GnRH neuron activation with only modest differences in AVPV neuron firing frequency. Of note, kisspeptin expression by these cells is increased during positive feedback (Smith et al., 2005; Gottsch et al., 2006). An important consideration for brain slice preparations is the disruption of neuronal projections that likely play a role in cycle-dependent changes. In rodents, input from the circadian clock is critical for surge induction (van der Beek, 1996); these and other afferents may be regulators of AVPV neuron activity in vivo and would not be present in the slices used (Watson et al., 1995).

In silico approaches provide the advantage of being able to test multiple hypotheses quickly but also have several caveats. First, while we sought to make the models reasonably comprehensive by including characterizations of currents reported in the literature, there are likely uncharacterized currents expressed in AVPV neurons that may change with cycle stage and/or influence firing behavior that were not included. Second, and related to the first point, no studies have characterized the transient sodium current in AVPV neurons. We adapted a Markov model of NaT used in a GnRH neuron model to enable firing in these models (Adams et al., 2018). Estimating this current was guided by fitting the rising phase of the action potential to experimental data, the phase of the action potential most highly influenced by sodium channel properties (Kress et al., 2010). Third, no computational model is ever a perfect representation, and it is possible that other combinations of properties for the components of this model could produce similar results in firing output. Despite these limitations, our voltage-clamp models are well fit, and the firing simulations reasonably approximated the firing behavior of these cells.

Together, these results point to changes in multiple voltage-dependent currents in these cells through the reproductive cycle. Our data and firing simulations suggest estrous cycle effects on these currents are unlikely to explain the increase in AVPV kisspeptin neuron spontaneous activity during the afternoon of proestrus, as changes in typically hyperpolarizing and depolarizing currents exert reciprocal effects on excitability, though the latter contribute to the emergence of postinhibitory rebound burst firing. These findings motivate further study of the inputs controlling firing behavior of these cells as well as the mechanisms regulating their neurochemical and peptidergic output to GnRH neurons, and how these contribute to GnRH/LH surge generation.

Acknowledgments

Acknowledgements: We thank Elizabeth Wagenmaker and Laura Burger for expert technical assistance and editorial comments. We also thank James L. Kenyon, University of Nevada, Reno, for the Excel spreadsheet used to calculate junction potentials and Srinivas Gorur-Shandilya for help with xolotl software.

Footnotes

  • The authors declare no competing financial interests.

  • This work was supported by the National Institute of Health/Eunice Kennedy Shriver National Institute of Child Health and Human Development Grant R01HD41469 (to S.M.M.). J.R.S. was supported by the National Institute of Health/Eunice Kennedy Shriver National Institute of Child Health and Human Development Grant F31HD097830.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.

References

  1. ↵
    Adachi S, Yamada S, Takatsu Y, Matsui H, Kinoshita M, Takase K, Sugiura H, Ohtaki T, Matsumoto H, Uenoyama Y, Tsukamura H, Inoue K, Maeda K-I (2007) Involvement of anteroventral periventricular metastin/kisspeptin neurons in estrogen positive feedback action on luteinizing hormone release in female rats. J Reprod Dev 53:367–378. doi:10.1262/jrd.18146 pmid:17213691
    OpenUrlCrossRefPubMed
  2. ↵
    Adams C, Stroberg W, DeFazio RA, Schnell S, Moenter SM (2018) Gonadotropin-releasing hormone (GnRH) neuron excitability is regulated by estradiol feedback and kisspeptin. J Neurosci 38:1249–1263. doi:10.1523/JNEUROSCI.2988-17.2017 pmid:29263236
    OpenUrlAbstract/FREE Full Text
  3. ↵
    Armstrong CM, Bezanilla F (1977) Inactivation of the sodium channel. II. Gating current experiments. J Gen Physiol 70:567–590. doi:10.1085/jgp.70.5.567 pmid:591912
    OpenUrlAbstract/FREE Full Text
  4. ↵
    Banciu A, Banciu DD, Mustaciosu CC, Radu M, Cretoiu D, Xiao J, Cretoiu SM, Suciu N, Radu BM (2018) Beta-estradiol regulates voltage-gated calcium channels and estrogen receptors in telocytes from human myometrium. Int J Mol Sci 19:1413. doi:10.3390/ijms19051413
    OpenUrlCrossRef
  5. ↵
    Barry PH (1994) JPCalc, a software package for calculating liquid junction potential corrections in patch-clamp, intracellular, epithelial and bilayer measurements and for correcting junction potential measurements. J Neurosci Methods 51:107–116. doi:10.1016/0165-0270(94)90031-0 pmid:8189746
    OpenUrlCrossRefPubMed
  6. ↵
    Bevan MD, Wilson CJ (1999) Mechanisms underlying spontaneous oscillation and rhythmic firing in rat subthalamic neurons. J Neurosci 19:7617–7628. doi:10.1523/JNEUROSCI.19-17-07617.1999
    OpenUrlAbstract/FREE Full Text
  7. ↵
    Bosch MA, Tonsfeldt KJ, Rønnekleiv OK (2013) mRNA expression of ion channels in GnRH neurons: subtype-specific regulation by 17β-estradiol. Mol Cell Endocrinol 367:85–97. doi:10.1016/j.mce.2012.12.021 pmid:23305677
    OpenUrlCrossRefPubMed
  8. ↵
    Bronson FH, Vom Saal FS (1979) Control of the preovulatory release of luteinizing hormone by steroids in the mouse. Endocrinology 104:1247–1255. doi:10.1210/endo-104-5-1247 pmid:571329
    OpenUrlCrossRefPubMed
  9. ↵
    Brown-Grant K, Exley D, Naftolin F (1970) Peripheral plasma oestradiol and luteinizing hormone concentrations during the oestrous cycle of the rat. J Endocrinol 48:295–296. doi:10.1677/joe.0.0480295 pmid:5528614
    OpenUrlAbstract/FREE Full Text
  10. ↵
    Caraty A, Lomet D, Sébert ME, Guillaume D, Beltramo M, Evans NP (2013) Gonadotrophin-releasing hormone release into the hypophyseal portal blood of the ewe mirrors both pulsatile and continuous intravenous infusion of kisspeptin: an insight into kisspeptin’s mechanism of action. J Neuroendocrinol 25:537–546. doi:10.1111/jne.12030 pmid:23387514
    OpenUrlCrossRefPubMed
  11. ↵
    Choquet D, Korn H (1992) Mechanism of 4-aminopyridine action on voltage-gated potassium channels in lymphocytes. J Gen Physiol 99:217–240. doi:10.1085/jgp.99.2.217 pmid:1613484
    OpenUrlAbstract/FREE Full Text
  12. ↵
    Christian CA, Mobley JL, Moenter SM (2005) Diurnal and estradiol-dependent changes in gonadotropin-releasing hormone neuron firing activity. Proc Natl Acad Sci USA 102:15682–15687. doi:10.1073/pnas.0504270102 pmid:16230634
    OpenUrlAbstract/FREE Full Text
  13. ↵
    Clay JR (2000) Determining K+ channel activation curves from K+ channel currents. Eur Biophys J 29:555–557. doi:10.1007/s002490000091 pmid:11156298
    OpenUrlCrossRefPubMed
  14. ↵
    Clay JR (2009) Determining k channel activation curves from k channel currents often requires the Goldman-Hodgkin-Katz equation. Front Cell Neurosci 3:20.
    OpenUrlCrossRefPubMed
  15. ↵
    Coetzee WA, Amarillo Y, Chiu J, Chow A, Lau D, McCormack T, Moreno H, Nadal MS, Ozaita A, Pountney D, Saganich M, Vega-Saenz de Miera E, Rudy B (1999) Molecular diversity of K+ channels. Ann NY Acad Sci 868:233–285. doi:10.1111/j.1749-6632.1999.tb11293.x pmid:10414301
    OpenUrlCrossRefPubMed
  16. ↵
    Cravo RM, Margatho LO, Osborne-Lawrence S, Donato J, Atkin S, Bookout AL, Rovinsky S, Frazão R, Lee CE, Gautron L, Zigman JM, Elias CF (2011) Characterization of Kiss1 neurons using transgenic mouse models. Neuroscience 173:37–56. doi:10.1016/j.neuroscience.2010.11.022 pmid:21093546
    OpenUrlCrossRefPubMed
  17. ↵
    Cropper EC, Jing J, Vilim FS, Weiss KR (2018) Peptide cotransmitters as dynamic, intrinsic modulators of network activity. Front Neural Circuits 12:78. doi:10.3389/fncir.2018.00078 pmid:30333732
    OpenUrlCrossRefPubMed
  18. ↵
    DeFazio RA, Moenter SM (2002) Estradiol feedback alters potassium currents and firing properties of gonadotropin-releasing hormone neurons. Mol Endocrinol 16:2255–2265. doi:10.1210/me.2002-0155 pmid:12351691
    OpenUrlCrossRefPubMed
  19. ↵
    DeFazio RA, Heger S, Ojeda SR, Moenter SM (2002) Activation of A-type gamma-aminobutyric acid receptors excites gonadotropin-releasing hormone neurons. Mol Endocrinol 16:2872–2891. doi:10.1210/me.2002-0163 pmid:12456806
    OpenUrlCrossRefPubMed
  20. ↵
    DeFazio RA, Elias CF, Moenter SM (2014) GABAergic transmission to kisspeptin neurons is differentially regulated by time of day and estradiol in female mice. J Neurosci 34:16296–16308. doi:10.1523/JNEUROSCI.3057-14.2014 pmid:25471569
    OpenUrlAbstract/FREE Full Text
  21. ↵
    DeFazio RA, Navarro MA, Adams CE, Milescu LS, Moenter SM (2019) Estradiol enhances the depolarizing response to GABA and AMPA synaptic conductances in arcuate kisspeptin neurons by diminishing voltage-gated potassium currents. J Neurosci 39:9532–9545. doi:10.1523/JNEUROSCI.0378-19.2019 pmid:31628184
    OpenUrlAbstract/FREE Full Text
  22. ↵
    Döcke F, Dörner G (1965) The mechanism of the induction of ovulation by oestrogens. J Endocrinol 33:491–499. doi:10.1677/joe.0.0330491 pmid:5891806
    OpenUrlAbstract/FREE Full Text
  23. ↵
    Du J, Wang Q, Hu F, Wang J, Ding H, Gao R, Xiao H, Wang L (2014) Effects of estradiol on voltage-gated potassium channels in mouse dorsal root ganglion neurons. J Membr Biol 247:541–548. doi:10.1007/s00232-014-9670-z pmid:24838692
    OpenUrlCrossRefPubMed
  24. ↵
    Dutton A, Dyball REJ (1979) Phasic firing enhances vasopressin release from the rat neurohypophysis. J Physiol 290:433–440. doi:10.1113/jphysiol.1979.sp012781
    OpenUrlCrossRefPubMed
  25. ↵
    Gold MS, Shuster MJ, Levine JD (1996) Characterization of six voltage-gated K+ currents in adult rat sensory neurons. J Neurophysiol 75:2629–2646. doi:10.1152/jn.1996.75.6.2629 pmid:8793767
    OpenUrlCrossRefPubMed
  26. ↵
    Goodman RL (1978) The site of the positive feedback action of estradiol in the rat. Endocrinology 102:151–159. doi:10.1210/endo-102-1-151 pmid:570478
    OpenUrlCrossRefPubMed
  27. ↵
    Gorur-Shandilya S, Hoyland A, Marder E (2018) Xolotl: an intuitive and approachable neuron and network simulator for research and teaching. Front Neuroinform 12:87. doi:10.3389/fninf.2018.00087 pmid:30534067
    OpenUrlCrossRefPubMed
  28. ↵
    Gottsch ML, Clifton DK, Steiner RA (2006) Kisspepeptin-GPR54 signaling in the neuroendocrine reproductive axis. Mol Cell Endocrinol 254–255:91–96. doi:10.1016/j.mce.2006.04.030 pmid:16762492
    OpenUrlCrossRefPubMed
  29. ↵
    Granados ST, Castillo K, Bravo-Moraga F, Sepúlveda RV, Carrasquel-Ursulaez W, Rojas M, Carmona E, Lorenzo-Ceballos Y, González-Nilo F, González C, Latorre R, Torres YP (2019) The molecular nature of the 17β-Estradiol binding site in the voltage- and Ca 2+ -activated K + (BK) channel β1 subunit. Sci Rep 9:9965. doi:10.1038/s41598-019-45942-1 pmid:31292456
    OpenUrlCrossRefPubMed
  30. ↵
    Greep RO, Van Dyke HB, Chow BF (1942) Gonadotropins of the swine pituitary I. Various biological effects of purified Thylakentrin (FSH) and pure metakentrin (ICSH). Endocrinology 30:635–649. doi:10.1210/endo-30-5-635
    OpenUrlCrossRef
  31. ↵
    Han SK, Gottsch ML, Lee KJ, Popa SM, Smith JT, Jakawich SK, Clifton DK, Steiner RA, Herbison AE (2005) Activation of gonadotropin-releasing hormone neurons by kisspeptin as a neuroendocrine switch for the onset of puberty. J Neurosci 25:11349–11356. doi:10.1523/JNEUROSCI.3328-05.2005
    OpenUrlAbstract/FREE Full Text
  32. ↵
    Hille B (2001) Ion channels of excitable membranes, Ed 3. Sunderland: Sinauer Associates.
  33. ↵
    Kalra PS, McCann SM (1975) The stimulatory effect on gonadotropin release of implants of estradiol or progesterone in certain sites in the central nervous system. Neuroendocrinology 19:289–302. doi:10.1159/000122450 pmid:1226226
    OpenUrlCrossRefPubMed
  34. ↵
    Kehl SJ (2017) A Model of the block of voltage-gated potassium Kv4.2 ionic currents by 4-aminopyridine. J Pharmacol Exp Ther 363:184–195. doi:10.1124/jpet.117.243097 pmid:28864468
    OpenUrlAbstract/FREE Full Text
  35. ↵
    Kelly MJ, Qiu J (2010) Estrogen signaling in hypothalamic circuits controlling reproduction. Brain Res 1364:44–52. doi:10.1016/j.brainres.2010.08.082 pmid:20807512
    OpenUrlCrossRefPubMed
  36. ↵
    Kimm T, Khaliq ZM, Bean BP (2015) Differential regulation of action potential shape and burst-frequency firing by BK and Kv2 channels in substantia nigra dopaminergic neurons. J Neurosci 35:16404–16417. doi:10.1523/JNEUROSCI.5291-14.2015 pmid:26674866
    OpenUrlAbstract/FREE Full Text
  37. ↵
    Kress GJ, Dowling MJ, Eisenman LN, Mennerick S (2010) Axonal sodium channel distribution shapes the depolarized action potential threshold of dentate granule neurons. Hippocampus 20:558–571. doi:10.1002/hipo.20667 pmid:19603521
    OpenUrlCrossRefPubMed
  38. ↵
    Lee SK, Ryu PD, Lee SY (2013) Estrogen replacement modulates voltage-gated potassium channels in rat presympathetic paraventricular nucleus neurons. BMC Neurosci 14:134. doi:10.1186/1471-2202-14-134 pmid:24180323
    OpenUrlCrossRefPubMed
  39. ↵
    Levitan IB (2006) Signaling protein complexes associated with neuronal ion channels. Nat Neurosci 9:305–310. doi:10.1038/nn1647 pmid:16498425
    OpenUrlCrossRefPubMed
  40. ↵
    Lüthi A, McCormick DA (1998) H-current: properties of a neuronal and network pacemaker. Neuron 21:9–12. doi:10.1016/s0896-6273(00)80509-7 pmid:9697847
    OpenUrlCrossRefPubMed
  41. ↵
    Marni F, Wang Y, Morishima M, Shimaoka T, Uchino T, Zheng M, Kaku T, Ono K (2009) 17Beta-estradiol modulates expression of low-voltage-activated Ca(V)3.2 T-type calcium channel via extracellularly regulated kinase pathway in cardiomyocytes. Endocrinology 150:879–888. doi:10.1210/en.2008-0645 pmid:18832095
    OpenUrlCrossRefPubMed
  42. ↵
    Mendonça PRF, Kyle V, Yeo S-H, Colledge WH, Robinson HPC (2018) Kv4.2 channel activity controls intrinsic firing dynamics of arcuate kisspeptin neurons. J Physiol 596:885–899. doi:10.1113/JP274474 pmid:29214635
    OpenUrlCrossRefPubMed
  43. ↵
    Messager S, Chatzidaki EE, Ma D, Hendrick AG, Zahn D, Dixon J, Thresher RR, Malinge I, Lomet D, Carlton MBL, Colledge WH, Caraty A, Aparicio SAJR (2005) Kisspeptin directly stimulates gonadotropin-releasing hormone release via G protein-coupled receptor 54. Proc Natl Acad Sci USA 102:1761–1766. doi:10.1073/pnas.0409330102 pmid:15665093
    OpenUrlAbstract/FREE Full Text
  44. ↵
    Moenter SM, Caraty A, Locatelli A, Karsch FJ (1991) Pattern of gonadotropin-releasing hormone (GnRH) secretion leading up to ovulation in the ewe: existence of a preovulatory GnRH surge. Endocrinology 129:1175–1182. doi:10.1210/endo-129-3-1175 pmid:1874164
    OpenUrlCrossRefPubMed
  45. ↵
    Moran S, Moenter SM, Khadra A (2016) A unified model for two modes of bursting in GnRH neurons. J Comput Neurosci 40:297–315. doi:10.1007/s10827-016-0598-4 pmid:26975615
    OpenUrlCrossRefPubMed
  46. ↵
    Oakley AE, Clifton DK, Steiner RA (2009) Kisspeptin signaling in the brain. Endocr Rev 30:713–743. doi:10.1210/er.2009-0005 pmid:19770291
    OpenUrlCrossRefPubMed
  47. ↵
    Perez-Reyes E (2003) Molecular physiology of low-voltage-activated t-type calcium channels. Physiol Rev 83:117–161. doi:10.1152/physrev.00018.2002 pmid:12506128
    OpenUrlCrossRefPubMed
  48. ↵
    Petersen SL, Barraclough CA (1989) Suppression of spontaneous LH surges in estrogen-treated ovariectomized rats by microimplants of antiestrogens into the preoptic brain. Brain Res 484:279–289. doi:10.1016/0006-8993(89)90371-5 pmid:2713688
    OpenUrlCrossRefPubMed
  49. ↵
    Petersen SL, Cheuk C, Hartman RD, Barraclough CA (1989) Medial preoptic microimplants of the antiestrogen, keoxifene, affect luteinizing hormone-releasing hormone mRNA levels, median eminence luteinizing hormone-releasing hormone concentrations and luteinizing hormone release in ovariectomized, estrogen-treated rats. J Neuroendocrinol 1:279–283. doi:10.1111/j.1365-2826.1989.tb00116.x pmid:19210441
    OpenUrlCrossRefPubMed
  50. ↵
    Pielecka-Fortuna J, Chu Z, Moenter SM (2008) Kisspeptin acts directly and indirectly to increase gonadotropin-releasing hormone neuron activity and its effects are modulated by estradiol. Endocrinology 149:1979–1986. doi:10.1210/en.2007-1365 pmid:18162521
    OpenUrlCrossRefPubMed
  51. ↵
    Pielecka-Fortuna J, DeFazio RA, Moenter SM (2011) Voltage-gated potassium currents are targets of diurnal changes in estradiol feedback regulation and kisspeptin action on gonadotropin-releasing hormone neurons in mice. Biol Reprod 85:987–995. doi:10.1095/biolreprod.111.093492 pmid:21778142
    OpenUrlCrossRefPubMed
  52. ↵
    Piet R, Boehm U, Herbison AE (2013) Estrous cycle plasticity in the hyperpolarization-activated current ih is mediated by circulating 17β-estradiol in preoptic area kisspeptin neurons. J Neurosci 33:10828–10839. doi:10.1523/JNEUROSCI.1021-13.2013 pmid:23804103
    OpenUrlAbstract/FREE Full Text
  53. ↵
    Piet R, Kalil B, McLennan T, Porteous R, Czieselsky K, Herbison AE (2018) Dominant neuropeptide cotransmission in kisspeptin-GABA regulation of GnRH neuron firing driving ovulation. J Neurosci 38:6310–6322. doi:10.1523/JNEUROSCI.0658-18.2018 pmid:29899026
    OpenUrlAbstract/FREE Full Text
  54. ↵
    Porteous R, Herbison AE (2019) Genetic deletion of Esr1 in the mouse preoptic area disrupts the LH surge and estrous cyclicity. Endocrinology 160:1821–1829. doi:10.1210/en.2019-00284
    OpenUrlCrossRef
  55. ↵
    Qiu J, Bosch MA, Jamali K, Xue C, Kelly MJ, Rønnekleiv OK (2006) Estrogen upregulates T-type calcium channels in the hypothalamus and pituitary. J Neurosci 26:11072–11082. doi:10.1523/JNEUROSCI.3229-06.2006 pmid:17065449
    OpenUrlAbstract/FREE Full Text
  56. ↵
    Roepke TA, Malyala A, Bosch MA, Kelly MJ, Rønnekleiv OK (2007) Estrogen regulation of genes important for K+ channel signaling in the arcuate nucleus. Endocrinology 148:4937–4951. doi:10.1210/en.2007-0605 pmid:17595223
    OpenUrlCrossRefPubMed
  57. ↵
    Rønnekleiv OK, Zhang C, Bosch MA, Kelly MJ (2015) Kisspeptin and gonadotropin-releasing hormone neuronal excitability: molecular mechanisms driven by 17β-estradiol. Neuroendocrinology 102:184–193. doi:10.1159/000370311 pmid:25612870
    OpenUrlCrossRefPubMed
  58. ↵
    Sarkar DK, Chiappa SA, Fink G, Sherwood NM (1976) Gonadotropin-releasing hormone surge in pro-oestrous rats. Nature 264:461–463. doi:10.1038/264461a0 pmid:794737
    OpenUrlCrossRefPubMed
  59. ↵
    Shim WS, Conaway M, Masamura S, Yue W, Wang JP, Kumar R, Santen RJ (2000) Estradiol hypersensitivity and mitogen-activated protein kinase expression in long-term estrogen deprived human breast cancer cells in vivo. Endocrinology 141:396–405. doi:10.1210/endo.141.1.7270
    OpenUrlCrossRefPubMed
  60. ↵
    Simerly RB (1998) Organization and regulation of sexually dimorphic neuroendocrine pathways. Behav Brain Res 92:195–203. doi:10.1016/s0166-4328(97)00191-5 pmid:9638961
    OpenUrlCrossRefPubMed
  61. ↵
    Smith JT, Cunningham MJ, Rissman EF, Clifton DK, Steiner RA (2005) Regulation of Kiss1 gene expression in the brain of the female mouse. Endocrinology 146:3686–3692. doi:10.1210/en.2005-0488 pmid:15919741
    OpenUrlCrossRefPubMed
  62. ↵
    Sun J, Chu Z, Moenter SM (2010) Diurnal in vivo and rapid in vitro effects of estradiol on voltage-gated calcium channels in gonadotropin-releasing hormone neurons. J Neurosci 30:3912–3923. doi:10.1523/JNEUROSCI.6256-09.2010 pmid:20237262
    OpenUrlAbstract/FREE Full Text
  63. ↵
    van der Beek EM (1996) Circadian control of reproduction in the female rat. Prog Brain Res 111:295–320. doi:10.1016/s0079-6123(08)60415-x pmid:8990922
    OpenUrlCrossRefPubMed
  64. ↵
    Vastagh C, Solymosi N, Farkas I, Liposits Z (2019) Proestrus differentially regulates expression of ion channel and calcium homeostasis genes in GnRH neurons of mice. Front Mol Neurosci 12:137. doi:10.3389/fnmol.2019.00137 pmid:31213979
    OpenUrlCrossRefPubMed
  65. ↵
    Wang L, DeFazio RA, Moenter SM (2016) Excitability and burst generation of AVPV kisspeptin neurons are regulated by the estrous cycle via multiple conductances modulated by estradiol action. eNeuro 3:ENEURO.0094-16.2016. doi:10.1523/ENEURO.0094-16.2016
    OpenUrlAbstract/FREE Full Text
  66. ↵
    Wang L, Vanacker C, Burger LL, Barnes T, Shah YM, Myers MG, Moenter SM (2019) Genetic dissection of the different roles of hypothalamic kisspeptin neurons in regulating female reproduction. Elife 8:e43999. doi:10.7554/eLife.43999
    OpenUrlCrossRefPubMed
  67. ↵
    Watson RE, Langub MC, Engle MG, Maley BE (1995) Estrogen-receptive neurons in the anteroventral periventricular nucleus are synaptic targets of the suprachiasmatic nucleus and peri-suprachiasmatic region. Brain Res 689:254–264. doi:10.1016/0006-8993(95)00548-5 pmid:7583329
    OpenUrlCrossRefPubMed
  68. ↵
    Wiegand SJ, Terasawa E (1982) Discrete lesions reveal functional heterogeneity of suprachiasmatic structures in regulation of gonadotropin secretion in the female rat. Neuroendocrinology 34:395–404. doi:10.1159/000123335 pmid:6808412
    OpenUrlCrossRefPubMed
  69. ↵
    Yang JA, Mamounis KJ, Yasrebi A, Roepke TA (2016) Regulation of gene expression by 17β-estradiol in the arcuate nucleus of the mouse through ERE-dependent and ERE-independent mechanisms. Steroids 107:128–138. doi:10.1016/j.steroids.2016.01.003 pmid:26768413
    OpenUrlCrossRefPubMed
  70. ↵
    Zhang C, Roepke TA, Kelly MJ, Rønnekleiv OK (2008) Kisspeptin depolarizes gonadotropin-releasing hormone neurons through activation of TRPC-like cationic channels. J Neurosci 28:4423–4434. doi:10.1523/JNEUROSCI.5352-07.2008 pmid:18434521
    OpenUrlAbstract/FREE Full Text
  71. ↵
    Zhang C, Bosch MA, Rick EA, Kelly MJ, Rønnekleiv OK (2009) 17Beta-estradiol regulation of T-type calcium channels in gonadotropin-releasing hormone neurons. J Neurosci 29:10552–10562. doi:10.1523/JNEUROSCI.2962-09.2009 pmid:19710308
    OpenUrlAbstract/FREE Full Text
  72. ↵
    Zhang C, Kelly MJ, Rønnekleiv OK (2010) 17 β-estradiol rapidly increases ATP-sensitive potassium channel activity in gonadotropin-releasing hormone neurons [corrected] via a protein kinase signaling pathway. Endocrinology 151:4477–4484. doi:10.1210/en.2010-0177 pmid:20660067
    OpenUrlCrossRefPubMed
  73. ↵
    Zhang C, Bosch MA, Qiu J, Rønnekleiv OK, Kelly MJ (2015) 17β-Estradiol increases persistent Na(+) current and excitability of AVPV/PeN Kiss1 neurons in female mice. Mol Endocrinol 29:518–527. doi:10.1210/me.2014-1392 pmid:25734516
    OpenUrlCrossRefPubMed

Synthesis

Reviewing Editor: Anna Menini, SISSA (International School of Advances Studies)

Decisions are customarily a result of the Reviewing Editor and the peer reviewers coming together and discussing their recommendations until a consensus is reached. When revisions are invited, a fact-based synthesis statement explaining their decision and outlining what is needed to prepare a revision will be listed below. Note: If this manuscript was transferred from JNeurosci and a decision was made to accept the manuscript without peer review, a brief statement to this effect will instead be what is listed below.

The authors answered reviewers’ comments and made the suggested revisions.

Back to top

In this issue

eneuro: 8 (5)
eNeuro
Vol. 8, Issue 5
September/October 2021
  • Table of Contents
  • Index by author
  • Ed Board (PDF)
Email

Thank you for sharing this eNeuro article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Reciprocal Changes in Voltage-Gated Potassium and Subthreshold Inward Currents Help Maintain Firing Dynamics of AVPV Kisspeptin Neurons during the Estrous Cycle
(Your Name) has forwarded a page to you from eNeuro
(Your Name) thought you would be interested in this article in eNeuro.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Reciprocal Changes in Voltage-Gated Potassium and Subthreshold Inward Currents Help Maintain Firing Dynamics of AVPV Kisspeptin Neurons during the Estrous Cycle
J. Rudolph Starrett, R. Anthony DeFazio, Suzanne M. Moenter
eNeuro 12 August 2021, 8 (5) ENEURO.0324-21.2021; DOI: 10.1523/ENEURO.0324-21.2021

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Share
Reciprocal Changes in Voltage-Gated Potassium and Subthreshold Inward Currents Help Maintain Firing Dynamics of AVPV Kisspeptin Neurons during the Estrous Cycle
J. Rudolph Starrett, R. Anthony DeFazio, Suzanne M. Moenter
eNeuro 12 August 2021, 8 (5) ENEURO.0324-21.2021; DOI: 10.1523/ENEURO.0324-21.2021
Reddit logo Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Significance Statement
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Acknowledgments
    • Footnotes
    • References
    • Synthesis
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • AVPV
  • estrous cycle
  • excitability
  • kisspeptin
  • positive feedback
  • potassium

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Research Article: New Research

  • Effects of Cortical FoxP1 Knockdowns on Learned Song Preference in Female Zebra Finches
  • Impaired AMPARs Translocation into Dendritic Spines with Motor Skill Learning in the Fragile X Mouse Model
  • Taste-Odor Association Learning Alters the Dynamics of Intraoral Odor Responses in the Posterior Piriform Cortex of Awake Rats
Show more Research Article: New Research

Neuronal Excitability

  • SK and Kv4 channels limit spike timing perturbations in pacemaking dopamine neurons
  • Gas7 is a novel dendritic spine initiation factor
  • CaMKIIα promoter-controlled circuit manipulations target both pyramidal cells and inhibitory interneurons in cortical networks
Show more Neuronal Excitability

Subjects

  • Neuronal Excitability

  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Latest Articles
  • Issue Archive
  • Blog
  • Browse by Topic

Information

  • For Authors
  • For the Media

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
  • Feedback
(eNeuro logo)
(SfN logo)

Copyright © 2023 by the Society for Neuroscience.
eNeuro eISSN: 2373-2822

The ideas and opinions expressed in eNeuro do not necessarily reflect those of SfN or the eNeuro Editorial Board. Publication of an advertisement or other product mention in eNeuro should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in eNeuro.