Abstract
Hypothalamic kisspeptin (Kiss1) neurons are vital for pubertal development and reproduction. Arcuate nucleus Kiss1 (Kiss1ARH) neurons are responsible for the pulsatile release of Gonadotropin-releasing Hormone (GnRH). In females, the behavior of Kiss1ARH neurons, expressing Kiss1, Neurokinin B (NKB), and Dynorphin (Dyn), varies throughout the ovarian cycle. Studies indicate that 17β-estradiol (E2) reduces peptide expression but increases Vglut2 mRNA and glutamate neurotransmission in these neurons, suggesting a shift from peptidergic to glutamatergic signaling. To investigate this shift, we combined transcriptomics, electrophysiology, and mathematical modeling. Our results demonstrate that E2 treatment upregulates the mRNA expression of voltage-activated calcium channels, elevating the whole-cell calcium current and that contribute to high-frequency burst firing. Additionally, E2 treatment decreased the mRNA levels of Canonical Transient Receptor Potential (TPRC) 5 and G protein-coupled K+ (GIRK) channels. When TRPC5 channels in Kiss1ARH neurons were deleted using CRISPR, the slow excitatory postsynaptic potential (sEPSP) was eliminated. Our data enabled us to formulate a biophysically realistic mathematical model of the Kiss1ARH neuron, suggesting that E2 modifies ionic conductances in Kiss1ARH neurons, enabling the transition from high frequency synchronous firing through NKB-driven activation of TRPC5 channels to a short bursting mode facilitating glutamate release. In a low E2 milieu, synchronous firing of Kiss1ARH neurons drives pulsatile release of GnRH, while the transition to burst firing with high, preovulatory levels of E2 would facilitate the GnRH surge through its glutamatergic synaptic connection to preoptic Kiss1 neurons.
Introduction
Hypothalamic kisspeptin (Kiss1) neurons and its cognate receptor (GPR 54 or Kiss1 R) are essential for pubertal development and reproduction, and may also be involved in the control of energy homeostasis (Kotani et al., 2001) (De Roux et al., 2003) (Seminara et al., 2003) (Messager et al., 2005) (Shahab et al., 2005) (d’Anglemont de Tassigny et al., 2008) (Qiu J. et al., 2018) (Rønnekleiv et al., 2022). Kisspeptin neurons within the arcuate nucleus of the hypothalamus (Kiss1ARH) co-express Kiss1, Neurokinin B (NKB) and Dynorphin (Dyn), which are all down-regulated by 17β-estradiol (E2) (Goodman et al., 2007; Navarro V. M. et al., 2009). However, Kiss1ARH neurons also express vesicular glutamate transporter 2 (vGlut2) and release glutamate, and both vGlut2 expression and glutamate release are upregulated by E2 in females (Qiu J. et al., 2018). This indicates that peptides and glutamate in Kiss1ARH neurons are differently modulated by E2 and suggests that there is a complex E2 regulation in these neurons such that they transition from predominantly peptidergic to glutamatergic neurotransmission and hence from a “pulsatile” to “surge” mode. It has been known for decades that neurons located within the arcuate nucleus are responsible for the pulsatile release of GnRH and subsequently pulsatile release of LH from the pituitary gland (O’Byrne et al., 1991) (Moenter et al., 1993). In this respect, it is now generally accepted that Kiss1ARH neurons are the main neurons responsible for the generation of pulsatile LH release, but the underlying cellular conductances generating this activity have not been elucidated.
Morphological studies have provided evidence that Kiss1ARH neurons can communicate directly with each other (Lehman et al., 2010) (Navarro V. M. et al., 2009) (Navarro V.M. et al., 2011). Furthermore, Kiss1ARH neurons express the NKB receptor, Tacr3, as well as the kappa (κ) opioid receptor (KOR), whereas GPR54 is not expressed in Kiss1ARH neurons, rendering them unresponsive to kisspeptin (d’Anglemont de Tassigny et al., 2008; Navarro V. M. et al., 2009; Wakabayashi et al., 2010). Using optogenetics and whole-cell recordings we demonstrated that high-frequency photoactivation of Kiss1ARH neurons induces an NKB-mediated slow excitatory postsynaptic potential (EPSP), which is mediated by the recruitment of canonical transient receptor potential (TRCP5) channels. The release of NKB is limited by co-released dynorphin, which acts presynaptically to inhibit further release. Together the two peptides cause synchronized firing of Kiss1ARH neurons (Kelly et al., 2018; Qiu J. et al., 2016; Qiu Jian et al., 2021), whereas kisspeptin and glutamate appear to be the main output signals from Kiss1ARH neurons (Qiu J. et al., 2016; 3XXX) (Voliotis et al., 2021) (Liu et al., 2021).
Although single action potential-generated calcium influx is sufficient to trigger the release of classical neurotransmitters such as glutamate, high frequency (10-20 Hz) is required for the release of neuropeptides such as kisspeptin, NKB and dynorphin (Qiu J. et al., 2016). Indeed, the slow EPSP, which underlies the synchronization, is similar to the “plateau potential” that has been described in hippocampal and cortical neurons (Arboit et al., 2020; Zhang Z. et al., 2011). Many neurons, including Kiss1ARH neurons, express the biophysical properties that allow them to continue to persistently fire even after a triggering synaptic event has subsided (Zylberberg and Strowbridge, 2017) (Qiu J. et al., 2016). Moreover, the intrinsic bi-stability of neurons that generates persistent firing activity has been linked to a calcium-activated, non-selective cation current (ICAN) (Zylberberg and Strowbridge, 2017), and TRPC channels, specifically TRPC5 channels, are thought to be responsible for the ICAN in cortical neurons (Zhang Z. et al., 2011). Therefore, we postulated that TacR3 activation via NKB drives influx of Ca+2 through TRPC5 channels leading to greater build-up of [Ca2+]i that facilitates the opening of more TRPC5 channels in a self-sustaining manner. Indeed, using the fast intracellular calcium chelator BAPTA, which has been shown to robustly inhibit TRPC5 channel activation in heterologous cells (Blair et al., 2009), we have been able to abolish the slow EPSP and persistent firing in Kiss1ARH neurons following optogenetic stimulation in female mice (Qiu Jian et al., 2021).
Although the expression of peptides in Kiss1ARH neurons are downregulated by high circulating levels (late follicular levels) of E2, the intrinsic excitability of and the glutamate release by Kiss1ARH neurons are increased by Cacna1g (Cav3.1, T-type calcium channel), Hcn1 and Hcn2 (Hyperpolarization-activated, Cyclic Nucleotide Gated channels) mRNA expression and Vglut2 mRNA expression, respectively (Qiu J. et al., 2018). Burst firing in CNS neurons, which efficiently releases fast amino acid transmitters like glutamate, is generated primarily by the T-type calcium channel current (IT) (e.g., in thalamic relay neurons), and the rhythmicity of this burst firing is dependent on the h-current (Ih) (for review, see (Lüthi and McCormick, 1998; Zagotta and Siegelbaum, 1996)). Ih is mediated by the HCN channel family, which includes channel subtypes 1-4, of which Hcn1 and Hcn2 are the main channels in Kiss1ARH neurons. Ih depolarizes neurons from hyperpolarized states, raising the membrane potential into the range of IT activation (Erickson et al., 1993a, 1993b; Kelly and Rønnekleiv, 1994; Lüthi and McCormick, 1998; Zhang C. et al., 2009). IT is mediated by the low-threshold voltage-gated calcium channels, CaV3.1-3.3 (for review see (Perez-Reyes, 2003)). IT initiates a transient Ca2+-driven depolarization above the threshold for action potential initiation (i.e., a low threshold spike) (Llinás, 1988; Tsien et al., 1987). This depolarization then drives neurons to fire an ensemble (burst) of Na+-driven action potentials.
Based on the above compelling evidence we postulated that Kiss1ARH neurons transition from peptidergic neurotransmission, driving the pulsatile release of GnRH via kisspeptin release into the median eminence, to glutamatergic transmission that facilitates in the preovulatory surge of GnRH (Lin et al., 2021) (Shen et al., 2022) through their projection to the anteroventral periventricular/periventricular nucleus Kiss1 (Kiss1AVPV/PeN) neurons (Qiu J. et al., 2016). Therefore, we initiated studies to thoroughly characterize effects of high circulating (late follicular) levels of E2 on the expression of the full complement of voltage-activated calcium channels (and currents) and the opposing K+ channels, involved in the repolarization, on the excitability of Kiss1ARH neurons. We performed whole-cell recordings and single cell RT-PCR analysis of Kiss1ARH neurons to determine which channels are involved in the physiological transition from peptidergic to glutamatergic neurotransmission. Our physiological findings were incorporated into a mathematical model that accounts for the E2 effects on the firing activity of Kiss1ARH neurons and validates our hypothesis that high levels of E2 facilitate the transition from peptidergic to glutamatergic neurotransmission.
Results
Voltage-activated Ca2+ channels contribute to burst firing of Kiss1ARH neurons
Our whole-cell current clamp recordings of Kiss1ARH neurons revealed that there is a significant transition of these Kiss1 neurons from a silent/tonic firing mode in the ovariectomized state to bursting/irregular firing with E2 treatment (Figure 1). Therefore, we sought to elucidate the cationic channels contributing to this physiological transition. Previously, we have shown that an increase in the intracellular calcium concentrations can potentiate TRPC5 channel current in POMC neurons (Qiu J. et al., 2010). Additionally, chelating intracellular calcium with BAPTA abolishes the slow excitatory postsynaptic potential (EPSP) and persistent firing in Kiss1ARH neurons (Qiu Jian et al., 2021). Hence, to investigate the contributions of voltage-activated calcium channels (VGCCs) to the increase in intracellular calcium, we measured the peak calcium current contributed by both the low and high voltage-activated calcium channels. The inward currents evoked by the voltage pulses (150-ms pulses starting from a holding potential of -80 mV with 10-mV increments) were identified as calcium (Ca2+) currents, as they were blocked by the universal VGCC inhibitor Cd2+ (200 μM) (McNally et al., 2020) (Figures 2A-E). The maximum total inward and Cd2+-sensitive currents reached their peak amplitudes at -10 mV. To differentiate between various calcium channel subtypes present in Kiss1ARH neurons, we applied selective antagonists individually, allowing us to isolate the drug-sensitive current for each cell. When we individually applied specific antagonists, we observed partial inhibition of whole-cell Ca2+ currents. Treatment with 10 μM nifedipine (Hiraizumi et al., 2008; Kato et al., 2003; Lee et al., 2002), an L-type Ca2+ channel inhibitor, resulted in a partial inhibition (26.1%) of the whole-cell calcium current. Additionally, application of 2 μM ω-conotoxin GVIA (conoGVIA, targeting N-type channels) (Lee et al., 2002), 200 nM ω-agatoxin IVA (AgaIVA, targeting P/Q-type channels) (Kato et al., 2003; McNally et al., 2020), 100 nM SNX-482 (targeting R-type channels) (Hiraizumi et al., 2008), or 1 μM TTA-P2 (TTAP2, targeting T-type channels) (McNally et al., 2020) also led to partial inhibition of the Ca2+ current (Figure 2F). The observed reduction in current with each inhibitor indicates the presence of all the major subtypes of Ca2+ currents in Kiss1ARH neurons. Among the inhibitors used, the largest components of the whole-cell calcium current were sensitive to nifedipine (26.1%), conoGVIA (25.1%), and SNX-482 (31.1%) (Figures 2A, B, D and 2F). Subsequently, we documented contributions from TTA-P2-sensitive channels, accounting for approximately 6.7% of the total Ca2+ current, and AgaIVA-sensitive channels, which constituted approximately 3.9% (Figures 2E, C). These findings indicate that high voltage-activated L-, N-, and R-type channels constitute the largest component of the voltage-activated Ca2+ current in Kiss1ARH neurons that not only increases the overall excitability but also greatly facilitate TRPC5 channel opening (Blair et al., 2009), which is the major downstream target of TACR3 activation by NKB (Qiu J. et al., 2016).
Voltage-activated Ca2+ channels contribute to generation of slow-EPSP in Kiss1ARH neurons
Our previous study utilizing optogenetics demonstrated that high frequency photostimulation of Kiss1ARH neurons releases NKB. This release of NKB induces slow excitatory postsynaptic potentials (EPSPs) that facilitates the recruitment of other Kiss1ARH neurons, resulting in synchronous firing of the Kiss1ARH neuronal population (Qiu J. et al., 2016). Additionally, chelating intracellular calcium with the fast chelator BAPTA abolishes the slow EPSP and persistent firing in Kiss1ARH neurons, highlighting the role of calcium signaling in these processes (Qiu Jian et al., 2021). To access the involvement of HVA channels in the generation of the slow EPSP, we conducted experiments where we blocked the L-type Ca2+ channels with nifedipine (10 μM) and the N- and P/Q-type Ca2+ channels with ω-conotoxin MVIIC (1 μM) (Nunemaker et al., 2003). We then measured the slow EPSP. Indeed, both nifedipine and ω-conotoxin MVIIC significantly inhibited the slow EPSP by 35.2% and 60.4%, respectively (Figure 3A, B, E). In addition, we used SNX (100 nM), which selectively blocks R-type Ca2+ channels, and the slow EPSP was reduced to 28.7% of its control value (Figure 3C, E). The selective T-channel blocker TTA-P2 (5 µM) inhibited the slow EPSP by 28.6% (Figure 3D, E). Since the pharmacological blockade of the various calcium channels could also affect the presynaptic release of neurotransmitters (i.e., NKB), we measured the direct response of Kiss1ARH neurons to the TACR3 agonist senktide in the presence or absence of cadmium, a universal calcium channel blocker. In synaptically isolated whole-cell recordings of Kiss1ARH neurons, Cd2+ significantly inhibited the inward current induced by senktide (Figures 3F-H), which supports our posit that plasma membrane voltage-activated calcium channels contribute to the TRPC5 channel activation. Therefore, it appears that calcium channels contribute to maintaining the sustained depolarization underlying the slow EPSP.
E2 increases the mRNA expression and the whole-cell current of voltage-activated Ca2+ channels
The neuropeptides NKB (tachykinin2, Tac2) and kisspeptin (Kiss1), which are expressed in Kiss1ARH neurons, are crucial for the pulsatile release of gonadotropin-releasing hormone (GnRH) and reproductive processes. E2 decreases the expression of Kiss1 and Tac2 mRNA in Kiss1ARH neurons but enhances the excitability of Kiss1ARH neurons by amplifying the expression of Cacna1g, Hcn1, and Hcn2 mRNA, as well as increasing T-type calcium currents and h-currents, respectively (Qiu, 2018, 20613}. Moreover, E2 drives Slc17a6 mRNA expression and enhances glutamatergic synaptic input to arcuate neurons and Kiss1AVPV/PeN neurons (Qiu J. et al., 2018). As a result, the E2-driven increase in Kiss1ARH neuronal excitability and glutamate neurotransmission may play a crucial role in triggering the surge of GnRH, ultimately leading to the LH surge (Lin et al., 2021) (Shen et al., 2022).
To assess the impact of E2 on the modulation of voltage-activated calcium channels and Kiss1ARH neuronal excitability, we employed real-time PCR (qPCR) to measure the relative expression levels of ion channel subtypes in Kiss1ARH neurons. We compared the expression in E2-treated females to those treated with oil, using the specific primers listed in Table 1. The quantification was conducted on pools of 10 neurons, as indicated in the Methods and Materials. In both oil- and E2-treated females, we quantified the expression of Cav1.2 (L-type), Cav2.1 (P/Q-type), Cav2.2 (N-type) and Cav2.3 (R-type) mRNAs. Remarkably, all of these mRNA transcripts exhibited increased expression levels in response to E2 treatment (Figure 4). Congruent with our previous findings (Qiu J. et al., 2018), we observed that the mRNA expression of Cav3.1 and Hcn1 was also upregulated in response to E2 treatment (Figure 4). These results suggested that E2 has a regulatory effect on the expression and function of all of these ion channels in Kiss1ARH neurons. To determine whether the increased mRNA expression translated into functional changes at the cellular level, we measured the whole-cell calcium current in Kiss1ARH neurons obtained from ovariectomized mice treated with either vehicle or E2. We discovered that E2 treatment led to a significant increase in the peak calcium current density in Kiss1ARH neurons (Figure 5A-D), which was recapitulated by our computational modeling (Figure 5E). These findings indicate that the upregulation of the mRNA expression of calcium channels by E2 translated to an augmented peak calcium current in Kiss1ARH neurons.
The largest components of the calcium currents in Kiss1ARH neurons from the E2-treated, ovariectomized females were found to be sensitive to nifedipine, conoGVIA, and SNX-482, accounting for approximately 24.9%, 24.6%, and 27.0% of the total current across cells, respectively, which is very similar to their contributions to the whole-cell current from vehicle-treated, ovariectomized females (Figure 3). In addition, contributions from TTA-P2-sensitive channels accounted for approximately 11.1% of the total Ca2+ current, while agaIVA-sensitive channels contributed to approximately 11.0% (Figure 5D). These results highlight the prevalence of L-, N-, and R-type calcium channels as the major contributors to the whole-cell calcium current in Kiss1ARH neurons from E2-treated, ovariectomized females, but also the involvement of T-type and P/Q-type channels.
E2 does not alter the kinetics of calcium channel activation or de-inactivation
In order to determine whether the increase in peak voltage-activated calcium current density was the result of E2 regulating calcium channel kinetics, mRNA expression or both, we examined the voltage dependence of activation and inactivation. By measuring the voltage dependence of activation, we assessed how E2 affects the ability of calcium channels to open in response to membrane potential changes. Similarly, by examining the voltage dependence of inactivation, we determined how E2 influences the inactivation kinetics of calcium channels. Based on our results, there was no difference in the voltage dependence of activation between cells from the vehicle-treated control group (V1/2 = -32.3 ± 2.1 mV; n = 13) and cells from estrogen-treated females (V1/2 = -33.6 ± 2.5 mV; n = 11). Similarly, there was not a significant difference in the voltage dependence of inactivation between control cells (V1/2 = -48.9 ± 4.8 mV; n = 6) and estrogen-treated cells (V1/2 = -44.1 ± 1.9 mV; n = 5) (Figure 6). Therefore, although long-term E2 treatment increased the mRNA expression of HVA calcium channels (Figure 4), it did not affect the channel kinetics in Kiss1ARH neurons (Figure 6). Also, E2 downregulated the expression of Kcnd2 mRNA encoding Kv4.2, which is expressed in Kiss1ARH neurons (Mendonça et al., 2018) and has similar kinetics of activation as the T-type calcium channels (Oil-treated, ovariectomized females relative mRNA expression: 1.053, n = 5 animals versus E2-treated expression: 0.5643, n=5; t-test p = 0.0061). This opposing K+ current would dampen the inward calcium current. Therefore, it appears that long-term E2 treatment does not modulate calcium channel kinetics but rather alters the mRNA expression to increase the calcium channel conductance.
BK, SK and KCNQ channels are involved in modulating excitability of Kiss1ARH neurons
In the brain, calcium plays a crucial role in sculpting neuronal firing by activating potassium channels, which subsequently influence neuronal behavior (Nicoll, 1988; Storm, 1990). Since HVA and LVA calcium channels were expressed in Kiss1ARH neurons, all of which contribute to the elevation of intracellular calcium concentration ([Ca2+]i) that facilitates TRPC5 channel opening (Blair et al., 2009), our next step involved measuring the changes in Ca2+-activated K+ channel conductances and assessing their mRNA expression. In various cell types increases in cytosolic calcium levels, whether resulting from extracellular influx or intracellular release, lead to the activation of plasma membrane calcium-dependent potassium channels (Sah Pankaj and Louise Faber, 2002). Similarly, in Kiss1ARH neurons, these channels would be activated by calcium influx through all four types of high voltage-gated calcium channels, as well as the low voltage-activated calcium channel, which are all active during action potential firing. The activity of Ca2+-activated K+ channels play a crucial role in numerous physiological processes, including secretion and the regulation of neuronal firing properties. Two main families of Ca2+-activated K+ channel channels have been characterized, distinguished by their biophysical and pharmacological properties. These families are known as BK (Big Conductance K+) and SK (Small Conductance K+) channels in the CNS (Kshatri et al., 2018). BK channels are known for their high potassium selectivity and large single channel conductance, typically ranging from 100 to 300 pS. Activation of BK channels requires both calcium binding and membrane depolarization (Blatz and Magleby, 1987; Marty, 1989; Sah P., 1996; Storm, 1990). On the other hand, SK channels are simply activated by increases in cytosolic calcium levels, with their half-maximal activation at 0.3 μM (Bond et al., 1999).
To investigate K+ currents, the cells were maintained at a holding potential of -70 mV while being exposed to blockers CNQX, AP5, picrotoxin and TTX. Subsequently, the membrane potential was stepped by depolarizing voltages, ranging from -60 mV to +40 mV in 10 mV increments, for a duration of 500 ms (Brereton et al., 2013). This protocol was employed to activate K+ currents (Figure 7A). First, we examined SK currents. The mean current density was determined at the end of the voltage pulses. In vehicle-treated, OVX females the application of the SK channel blocker apamin (100 nM) (Spergel, 2007) led to a significant reduction in whole-cell currents in the +10 to +40 mV range (Figures 7A, B). The mean outward current density at +40 mV in the control group was 125.5 ± 7.2 pA/pF (n = 5) with the apamin-sensitive component contributing 60.2 ± 3.3 pA/pF (n = 5) (Figures 7A, C). In contrast, in the E2-treated females, the overall mean outward current density at +40 mV was 195.4 ± 14.6 pA/pF (n = 6), which was significantly greater than the vehicle control group (Figures 7D, E). However, there was no significant difference in the apamin-sensitive component between the vehicle-treated and E2-treated females, 60.2 ± 3.3 pA/pF (n = 5) versus 75.9 ± 12.3 pA/pF (n = 6), respectively (Figure 7F). Furthermore, to investigate the expression of the mRNAs encoding SK channel subunits in Kiss1ARH neurons from vehicle-treated and E2-treated OVX females, qPCR experiments were performed on 10-cell Kiss1ARH neuronal pools (Figure 7G). We focused on SK3 channels because these channels exhibit the highest expression in the hypothalamus and E2 regulates their expression (Bosch et al., 2002). E2 treatment had no effect on the mRNA expression of the SK3 subunit. These molecular findings support our electrophysiology results. Our computational model was calibrated so that SK channels contributed ∼50 pA/pF to the whole-cell outward K+ current in E2-treated females (Figure 7H).
Additionally, following the same protocol and in the presence of the same cocktail of blockers (CNQX, AP5, picrotoxin and TTX), we investigated the contribution of BK channels to Kiss1ARH neuronal excitability. In the OVX females, the application of the BK channel blocker iberiotoxin (ibTx, 200 nM) (Niday and Bean, 2021) resulted in only a slight attenuation of the outward current (n = 5) (Figures 8A, B). The ibTx-sensitive current density measured at +40 mV was 32.8 ± 12.0 pA/pF (Figures 8A, C). However, in the E2-treated females, the application of ibTx significantly attenuated the whole-cell K+ current from +20 to +40 mV, (Figures 8D, E). Additionally, the ibTx-sensitive current was significantly larger in the +0 to +40 mV range in the E2-treated females compared to the OVX females (Figure 8F). These findings indicate that E2 treatment modulates the activity of ibTx-sensitive BK current in Kiss1ARH neurons, resulting in increased current density (90.0 ± 10.1 pA/pF (n = 6) versus 31.1 ± 8.4 pA/pF (n = 5) at +40 mV). To investigate the expression of mRNA encoding BK channel subunits in Kiss1ARH neurons from vehicle-treated and E2-treated OVX females, qPCR experiments were performed on 10-cell Kiss1ARH neuronal pools (Figure 8G). E2 treatment significantly increased the mRNA expression of the BKα1 (Kcnma1) subunit. These findings support our electrophysiological findings that there is a significant increase in BK channel activity in Kiss1ARH neurons with E2 treatment (Figure 8F). In addition, E2 increased the mRNA expression of Kcnb1 encoding delayed rectifier Kv2.1 (E2-treated relative mRNA expression: 1.672, n = 5 versus oil-treated mRNA expression: 1.086, n = 5; t-test p-value = 0.0024). The up-regulation of the two of these K+ channels (BKα1 and Kv2.1) would facilitate rapid repolarization of Kiss1ARH following an action potential. Our computational model was calibrated so that BK channels contributed ∼100 pA/pF to the whole cell outward K+ current in the E2-treated females (Figures 8H).
Traditionally, the after hyperpolarization is divided into three distinct phases: fast (fAHP), medium (mAHP), and slow after hyperpolarization (sAHP) (Storm, 1990; Vogalis et al., 2003). The fast after hyperpolarization (fAHP) is primarily mediated by the BK family of potassium channels (Storm, 1987). The medium after hyperpolarization (mAHP) is predominantly mediated by apamin-sensitive SK2 channels (Bond et al., 2004; Peters et al., 2005). However, KCNQ family members contribute to both the mAHP and sAHP (Tzingounis A. V. et al., 2010; Tzingounis A. V. and Nicoll, 2008). Therefore, to investigate the contribution of KCNQ channels to Kiss1ARH neuronal excitability, voltage clamp experiments were conducted in the presence of TTX, CNQX, AP5, and picrotoxin, and a standard M-current protocol was run using the M-channel blocker XE-991 to isolate the M-current (Greene et al., 2017; Roepke et al., 2011) (Figure 9 A, B). The application of XE-991 resulted in the inhibition of M-current within a physiologically relevant voltage range of -60 to -30 mV in E2-treated OVX females but exhibited minimal impact in OVX females (Figure 9 C, D). Although the XE991-sensitive current was relatively small compared to other voltage-activated K+ conductances, it demonstrated a significant increase in E2-treated, OVX females (Figures 9 E). The maximum peak current density sensitive to XE-991 at -30 mV was found to be four times higher in E2-treated OVX females when compared to OVX females. This would contribute to the repolarization following burst firing. Furthermore, E2 increased the mRNA expression of Kcnq2, (Figure 9F), which suggests that KCNQ channels play a role in repolarizing Kiss1ARH neurons following burst firing. Indeed, our modeling predicted that M-current contributed to the repolarization following burst firing (Figure 9G).
E2 increases Vglut2 but down regulates Tac2, Trpc5 and Girk2 mRNA expression in Kiss1ARH neurons
Based on our electrophysiological results, Kiss1ARH neurons appear to transition from peptidergic to glutamatergic neurotransmission through E2-mediated changes in the expression of voltage-activated Ca2+ channels and K+ channels and their respective conductances. Therefore, we asked the question is there a difference in peptide and glutamate mRNA expression mediating this transition? Therefore, we ran a comparison between Tac 2 (NKB) and Vglut2 (surrogate for glutamate) expression. The cycle threshold (CT) was compared between Tac2 and Vglut2 as well as Kiss1, TRPC5 and GIRK2 in Kiss neuronal cell pools from OVX oil-treated and OVX E2-treated animals (Figures 10A, B). It is worth noting that lower number of cycles illustrate a higher quantity of mRNA expression because the fluorescence is detected earlier, and one cycle difference represents a doubling in expression. As expected, the reference gene Gapdh did not change with E2-treatment. However, quantitative PCR results revealed that E2 treatment of OVX females significantly reduced Tac2 expression (Figure 10C), whereas Vglut2 mRNA was significantly increased in Kiss1ARH neurons (Figure 10D). Moreover, both Trpc5 and Girk2 expression were significantly reduced in E2-treated, OVX females (Figures 10E, F).
CRISPR mutagenesis of Trpc5 attenuates slow EPSP and reduces excitability of Kiss1ARH neurons
Our previous (Qiu Jian et al., 2021) and current findings suggested that TRPC5 channels play a dominant role in regulating cell excitability in ovariectomized females. Therefore, as proof of principle, we utilized a CRISPR approach to mutate TRPC5 channels in Kiss1ARH neurons similar to our previous studies (Hunker et al., 2020; Stincic et al., 2021). Hunker et al. developed a single viral vector for conditional expression of the smaller Staphylococcus aureus (SaCas9) and sgRNA that yields high-efficiency mutagenesis in specific cell types (Hunker et al., 2020). To selectively mutate Trpc5 in Kiss1ARH neurons, we generated two guide RNA’s, one targeting exon 2, which is conserved across all splice variants, and the other targeting exon 7, the pore forming domain (Figures 11A, B). A cohort of Kiss1ARH mice were given bilateral stereotaxic injections into the ARH of the two AAV1-FLEX-SaCas9-sgTrpc5’s or a control virus containing with AAV1-FLEX-SaCas9-U6-sgRosa26.
An additional Cre-dependent virus of the same serotype (AAV1) that drove expression of a fluorophore (YFP or mCherry) was co-administered in order to visualize injection quality and facilitate harvesting of cells (Figure 11C). After three weeks, mice underwent ovariectomy since OVX mice express the maximum slow EPSP amplitude (Qiu J. et al., 2016). Brain slices were prepared, and cells harvested as previously described (Qiu J. et al., 2018) and analyzed with qPCR. We found that the Trpc5 mutagenesis group displayed a reduction in relative expression of Trpc5 in Kiss1ARH neurons compared to the control group (Figure 11D). Hence, the qPCR data verified that in the sgTrpc5-targeted mice we can selectively reduce Trpc5 gene expression in targeted cells.
As predicted, mutagenesis of Trpc5 channels in Kiss1ARH neurons significantly attenuated the slow EPSP (Figures 12A, B, C) such that the postsynaptic excitation was reduced to a “trickle” of action potential firing. What we would not have predicted is that the double sgRNA mutagenesis of Trpc5 channels in Kiss1ARH neurons significantly hyperpolarized the resting membrane potential by 7 mV (Figure 12D). Moreover, the rheobase (minimum current required to induce firing) significantly increased by ∼13 % in females bearing the sgTrpc5 double mutagenesis (Figure 12E). The firing frequency versus injected current (F-I) curve for sgTrpc5 double mutagenesis Kiss1ARH neurons was also significantly attenuated (Figure 12F). In agreement with these experimental findings, current ramp simulations of our mathematical model confirmed that a decrease of the TPRC5 conductance increased the rheobase and reduce the firing activity of Kiss1ARH neurons (Figures 12G, H). In our model, this contribution of TRPC5 channels to cell excitability is attributed to intracellular calcium levels which activate TRPC5 channels independently of NKB stimulation. Finally, we employed our mathematical model to further investigate the transition from synchronous firing driven by NKB release and TRPC5 channel activation to burst firing generated by E2-mediated upregulation of endogenous conductances. As we varied the conductance of TRCP5 channels and GIRK channels under saturating extracellular signals (e.g., NKB and Dyn), we find that balance between the two conductances controls the neuronal excitability in the OVX parameterized model, with synchronous firing eliminated when TRPC5 conductance is low (red triangle) relative to the GIRK conductance (Figure 13A). Furthermore, the burst firing of the OVX+E2 parameterized model was supported by elevated h- and Ca 2+-currents (Figure 13B) as well as by the high conductance of Ca2+ channels relative to the conductance of TRPC5 channels (Figure 13C). Importantly, the simulated burst firing was analogous to the phasic burst firing recorded ex vivo in Kiss1ARH neurons (Figure 1A).
Discussion
E2 appears to play a critical role in transitioning the glutamatergic/peptidergic Kiss1ARH neurons from a high frequency firing mode during synchronization, which is dependent on NKB-driven activation of TRPC5 channels, to a short bursting mode that would facilitate glutamate release. E2 decreased the expression of the peptide neurotransmitters NKB (kisspeptin and dynorphin) and TRPC5 channels but increased the mRNA expression of Vglut2, voltage-activated calcium channels and Hcn channels that contribute to burst firing and glutamate release from the Kiss1ARH neurons. We determined that the increase in mRNA expression of the HVA calcium channels translated into a significant increase in whole-cell current with all of the calcium channels contributing proportionally. Most importantly the kinetics of activation and inactivation were unaltered with E2 treatment, which indicates that other post-translation modifications were not affecting channel activity. Surprisingly and somewhat counter intuitive, the BK α1 subunit was also upregulated, but based on our modeling the rapid repolarization of the Kiss1ARH neurons (i.e., the fast AHP) would facilitate higher frequency of action potential firing. Moreover, our modeling confirmed that TRPC5 channels, which generate the slow EPSP (a.k.a., plateau potential in other CNS neurons), are vital for initiating and sustaining synchronous firing of Kiss1ARH neurons, while subsequent activation of GIRK channels repolarizes Kiss1ARH neurons. E2 treatment of ovariectomized females decreased both Trpc5 and Girk2 channel mRNA expression, which in our model correlated with the reduction in sustained high frequency firing of Kiss1ARH neurons. Therefore, the synchronous high frequency firing of Kiss1ARH neurons in a low E2 milieu correlates with the pulsatile release of GnRH (LH from the pituitary gland), whereas the transition to burst firing in the presence of high circulating levels of E2 (e.g., proestrus) would facilitate the GnRH (LH) surge through its glutamatergic synaptic connection with Kiss1AVPV/PeN neurons (Lin et al., 2021; Qiu J. et al., 2016).
Core calcium conductances underlying synchronous and burst firing of Kiss1ARH neurons
TRPC5 channels are highly expressed in Kiss1ARH neurons (Figure 11), and TRPC5 channels are essentially ligand-activated calcium channels with a high permeability to calcium (PCa/PNa = 9:1) (Venkatachalam and Montell, 2007). In general, mammalian TRPC channels are activated by both G protein-coupled receptors and receptor tyrosine kinases (Ambudkar and Ong, 2007; Clapham, 2003), and are one of the major downstream effectors activated by glutamate binding to group I metabotropic glutamate receptors (mGluR1 and mGluR5) in CNS neurons (Bengtson et al., 2004; Berg et al., 2007; Faber et al., 2006; Tozzi et al., 2003). In substantia nigra dopamine neurons mGluR1 agonists induce a current that exhibits the tell-tale double-rectifying current-voltage plot of TRPC channel activation (Tozzi et al., 2003), similar to what we see with the effects of the NKB agonist senktide in Kiss1ARH neurons (Qiu Jian et al., 2021). Both mGluR1 and TacR3 are Gq-coupled to phospholipase C (PLC) activation which leads to hydrolysis of phosphatidylinositol 4,5-bisphosphate (PIP2) to diacylglycerol (DAG) and inositol 1,4,5 triphosphate (IP3), which is involved in channel activation (Birnbaumer, 2009). TacR3 (NKB) signaling has additional consequences since many K+ (e.g., GIRK, KCNQ) channels are dependent on PIP2 for channel opening, and depletion of PIP2 by PLC leads to channel closure (Brown and Passmore, 2009) (Zhang C. et al., 2013) (Whorton and MacKinnon, 2011) (Zheng et al., 2022). Therefore, depletion of PIP2 by NKB signaling would further facilitate the sustained firing of Kiss1ARH neurons during synchronization.
A plateau potential has been characterized in hippocampal and cortical neurons (Arboit et al., 2020; Zhang Z. et al., 2011) as such these neurons express biophysical properties that allow them to continue to persistently fire even after a triggering synaptic event has subsided (Zylberberg and Strowbridge, 2017). The persistent firing activity of these neurons is linked to ICAN (Zylberberg and Strowbridge, 2017), and TRPC5 channels appear to be responsible for the ICAN (Zhang Z. et al., 2011). With TacR3 activation in Kiss1ARH neurons there is an influx of Ca+2 through TRPC5 channels leading to greater build-up of [Ca2+]i that facilitates the opening of more TRPC5 channels in a self-sustaining (autocatalytic) manner (Qiu J. et al., 2016). Using the fast intracellular calcium chelator BAPTA, which has been shown to robustly inhibit TRPC5 channel activation in heterologous cells (Blair et al., 2009), we abolished the slow EPSP and persistent firing of Kiss1ARH neurons following optogenetic stimulation (Qiu Jian et al., 2021). Moreover, HVA calcium channel blockers attenuated the generation of the slow EPSP (Figure 3) so it appears that they also contribute to the ICAN since calcium influx via both LVA and HVA calcium channels can also facilitate TRPC5 channel opening in Kiss1ARH neurons. In the ovariectomized female, treatment with E2 upregulated Cav1.2, Cav2.1, Cav2.2, and Cav2.3 mRNA by 1.5-to 2-fold and Cav3.1 mRNA expression by ∼3-fold. Hence, E2 significantly increased whole-cell calcium currents Kiss1ARH neurons, which greatly enhanced the excitability and contributed to the burst firing of Kiss1ARH neurons (Figure 1). However, the amplitude of the slow EPSP with E2 treatment is only ∼25% of the amplitude in the ovariectomized state (Qiu J. et al., 2016). Therefore, there appears to be a physiologic transition of Kiss1ARH neurons from the slow EPSP firing mode in the OVX state to the burst firing mode in the presence of E2, which has important physiologic ramifications as discussed below.
TRPC5 and GIRK channels are vital for synchronization of Kiss1ARH neurons
Recently, Tian et al. demonstrated that TRPC4, a close homolog of TRPC5 sharing ∼64 percent homology, is a “coincidence detector” of neurotransmission by both Gq/11 and Gi/o-coupled receptors in lateral septal (LS) neurons (Tian et al., 2022). In whole-cell recordings of LS neurons, TRPC 4 channels mediate a strong depolarizing plateau potential that in contrast to TRPC5 channel activation in Kiss1ARH neurons, abrogates action potential firing as a result of a depolarization block. In many instances the plateau potential in LS neurons is followed by an AHP, which is dependent on the activation of Gi/o-coupled receptors. In contrast, we have not observed an AHP in Kiss1ARH neurons following the slow EPSP (Qiu J. et al., 2016). Tian and colleagues showed that the depolarizing plateau in LS neurons is codependent on activation of both Gq/11-coupled mGluR1 glutamate receptors and Gi/o-coupled γ-aminobutyric acid type B receptors, the latter activating GIRK channels. Moreover, the firing patterns in LS neurons encodes information about the relative strengths of these contrasting inputs (i.e., Gq/11 versus Gi/o) such that only mGluR1 produces weak depolarization accompanied by increased firing of LS neurons, whereas pure GABAB receptor activation hyperpolarizes the cells and abrogates firing activity. Coincident input of both mGluR1 and GABAB receptors results in a brief burst of action potentials followed by a pause in firing, and both the pause duration and firing recovery patterns reflect the relative strengths of Gq/11 versus Gi/o inputs. Importantly, Tian and colleagues computationally simulated these various scenarios with computational modeling, and similar to our modeling, utilized only TRPC4 and GIRK channels. A notable difference between the Kiss1ARH neurons and the LS neurons is that dynorphin binds to kappa-opioid receptors to open GIRK channels in the nerve terminal of Kiss1ARH neurons (presynaptic effect) (Qiu J. et al., 2016), whereas GABA binds to some GABAB receptors to hyperpolarize LS neurons (Tian et al., 2022). Therefore, although the high frequency firing activity of both LS and Kiss1ARH neurons can be modeled around TRPC and GIRK channels, the generated firing patterns are dramatically different based on the timing (co-incident activation of TRPC4 and GIRK channels in LS) and localization of the kappa opioid-coupled GIRK channels in the axon terminal of Kiss1ARH neurons. Moreover, E2-treated, ovariectomized females show a significant down-regulation of both Trpc5 and Girk2 mRNA expression in Kiss1ARH neurons (Figure 10), which is important for the physiological transitioning as described below.
Interestingly, the hypothalamic A12 (ARH) dopamine neurons show a rhythmic “oscillatory” firing behavior that transitions to a tonic firing mode with synaptic input from Thyrotropin-releasing hormone (TRH) neurons (Lyons D.J. et al., 2010) or feedback by circulating prolactin, released by pituitary lactotrophs (Lyons D. J. et al., 2012). A more recent paper has revealed that TRPC5 channels mediate the plateau potential and tonic firing in A12 dopamine neurons in response to prolactin (Blum et al., 2019). Similar to our findings with CRISPR deletion of Trpc5 in Kiss1ARH neurons (Figures 11 and 12), conditional knockout of Trpc5 in dopamine neurons abrogated the prolactin-induced plateau potential and tonic firing. Although Blum and colleagues did not model the oscillatory firing or tonic firing of the A12 dopamine neurons, their findings are consistent with our results showing that the activation of TRPC5 channels underlies the slow EPSP (plateau potential) and sustained firing.
Contribution of endogenous K+ channels to synchronized and burst firing
Beyond the ligand-gated (e.g., baclofen) GIRK channels, there are endogenous K+ channels that help sculpt the firing activity of kisspeptin neurons. We focused on the calcium-activated K+ channel family: the large-conductance, calcium-activated potassium (BK, also called BKCa, KCa1.1, MaxiK, Slo), small conductance, calcium-activated K+ (SK1, SK2, SK3) (Bond et al., 2005), and the K+ channels underlying the M-current (KCNQ, Kv7.1-7.5) (Brown and Passmore, 2009), which mediate the fast afterhyperpolarization (AHP), the intermediate AHP/slow AHP, respectively (Andrade et al., 2012) BK channels are gated by both voltage and cytoplasmic calcium and sculpt action potential firing in CNS neurons (Blatz and Magleby, 1987; Marty, 1989; Sah P., 1996; Storm, 1990). Indeed, BK channels have been shown to mediate rapid spike repolarization—i.e., the fast AHP in hippocampal CA1 pyramidal neurons (Lancaster and Nicoll, 1987; Storm, 1987). Blockade of BK channels in CA1 neurons attenuates the initial discharge frequency in response to current injection, which is attributable to suppression of the BK channel-dependent rapid spike repolarization (Lancaster and Nicoll, 1987; Storm, 1987). Blockade of BK channels is thought to increase inactivation of the spike-generating transient Na+ current and activate more of the slower K+ currents, thereby enhancing refractoriness and reducing excitability during the immediate aftermath of the first action potential (Shao et al., 1999). Thus, BK channels facilitate high-frequency burst firing of CA1 neurons. Furthermore, extracellular field recordings confirmed that BK channels contribute to high-frequency burst firing in response to excitatory synaptic input to distal dendrites in CA1 neurons (Gu et al., 2007). Therefore, BK channels appear to play an important role for early high-frequency, rapidly adapting firing in hippocampal CA1 pyramidal neurons, thus promoting the type of bursting that is characteristic of these cells in vivo during behavior. Based on our in vitro electrophysiological recordings and computational modeling we see a similar physiological phenomenon in Kiss1ARH neurons (Figure 8). Not only does E2 increase the mRNA expression of Kcnma1 (Figure 8G), but also the maximum BK (IbTx-sensitive) current by 3-fold, which would facilitate rapid repolarization during burst firing and promote glutamate release similar to hippocampal CA1 neurons.
In contrast to BK channel expression, E2 did not affect the mRNA expression of SK3 channel mRNA. SK channels underlie the apamin-sensitive component of the medium duration AHP and are responsible for repolarization following a burst of action potentials (Andrade et al., 2012; Bond et al., 2005). The activation of SK channels is voltage-independent, but SK channels have a higher affinity for Ca2+ than BK channels (Andrade et al., 2012). SK channels are tightly coupled (within 100 nm) to L-type Ca2+ channels, and BK channels (within 30 nm) are tightly coupled to N-type Ca2+ channels in hippocampal CA1 pyramidal neurons (Marrion and Tavalin, 1998). The determination of the proximity of SK and BK channels to HVA calcium channels in kisspeptin neurons will require cell-attached patch recordings. However, in Kiss1ARH neurons the SK channels may come into play during a short burst of action potentials but would become overwhelmed with the higher frequency synchronized, sustained firing as a result of NKB stimulation (Qiu J. et al., 2016). As discussed above what limits the synchronized firing of Kiss1ARH neurons is the activation of GIRK channels.
Finally, the calcium-activated slow AHP probably plays a critical role in the repolarization of Kiss1ARH after burst firing. The molecular identification of the channels mediating the slow AHP has long been an area of intense investigation (Andrade et al., 2012; Vogalis et al., 2003). A critical feature of the slow AHP is that it activates very slowly (hundreds of milliseconds) long after the rise in cytoplasmic Ca2+ (Sah P. and Clements, 1999) so an intermediate Ca2+ signaling molecule has long been thought to be involved. Indeed, Tzingounis and colleagues (Tzingounis A.V. et al., 2007) have provided compelling evidence that the diffusible calcium sensor hippocalcin is the critical intermediate molecule involved in Ca2+ sensing. The slow AHP is abrogated in hippocalcin KO mice, and transfection of hippocalcin into cultured hippocampal neurons generates a pronounced slow AHP in response to a depolarizing stimulus (Tzingounis A.V. et al., 2007). Importantly, the slow AHP is activated by Ca2+ with an EC50 ≈ 300 nM, which is well within the operational range of hippocalcin but well below that of calmodulin (Andrade et al., 2012). Finally, two seminal papers from Tzingounis and colleagues demonstrate that KCNQ 2, 3 channels are responsible for the slow AHP in hippocampal dentate neurons (Tzingounis A. V. and Nicoll, 2008), and KCNQ 5 channels are responsible for the slow AHP in CA3 neurons (Tzingounis A. V. et al., 2010). Moreover, the KCNQ channel blocker XE991 attenuates the slow AHP in CA3 neurons (Tzingounis A. V. and Nicoll, 2008). Based on these seminal findings we investigated the role of the KCNQ channels, which “classically” underlie the M-current in Kiss1ARH. The M-current was first identified in Kiss1ARH neurons by Conde and Roepke (Conde and Roepke, 2019), and presently we found that Kcnq2 mRNA is expressed in Kiss1ARH neurons and up-regulated by E2, which translated into a greater M-current in Kiss1ARH neurons (Figure 9). Incorporating the M-current into our computational model indeed supports our hypothesis that M-current, along with ISK and IBK, contributes to membrane repolarization after burst firing (Figure 9G). Importantly the slow AHP, as opposed to the fast AHP (BK) and medium AHP (SK) is highly regulated by multiple neurotransmitters (Andrade et al., 2012), which sets the stage for further modulation of the slow EPSP in Kiss1ARH neurons.
Importance of E2-driven physiological transitioning
Since the expression of the peptide neurotransmitters in Kiss1ARH neurons are down-regulated by E2, the Kiss1ARH neurons are believed to be under “inhibitory” control by E2 and are important for “negative-feedback” regulation of GnRH release (Lehman et al., 2013; Navarro V. M. et al., 2009; Smith et al., 2005) (Rance Naomi E. and Young, 1991) (Rance N.E., 2009). However, our past (Gottsch et al., 2011) and current findings document that these Kiss1ARH neurons express HVA and LVA calcium and HCN (pacemaker) channels and are excited by co-released glutamate from neighboring Kiss1ARH neurons, which indicates that these neurons have pacemaker electrophysiological properties similar to other CNS neurons (Bal and McCormick, 1993; Lüthi and McCormick, 1998). Additionally, in contrast to the neuropeptides, E2 increases Slc17a6 (Vglut2) mRNA expression, in addition to mRNA for the voltage-activated calcium and HCN channels, and increases glutamate release onto Kiss1AVPV/PeN neurons (Qiu J. et al., 2018). Interestingly, Slc17a6 mRNA expression in Kiss1ARH neurons and the probability of glutamate release are decreased along with the neuropeptides in intact versus castrated males (Nestor et al., 2016), which indicates a profound sex difference in the glutamate signaling by Kiss1ARH neurons (Nestor et al., 2016; Qiu J. et al., 2018). Obviously, in the male there is no preovulatory LH surge so there is no need for excitatory glutamatergic input to the few Kiss1AVPV neurons in the male.
In females, conditional knockout of Slc17a6 in Kiss1 neurons abrogates glutamate release from Kiss1ARH neurons (Qiu J. et al., 2018). Kiss1AVPV/PeN neurons do not express Slc17a6 and do not release glutamate. Within the Kiss1ARH neurocircuitry the lack of glutamate transmission does not diminish the slow EPSP in ovariectomized females (Qiu J. et al., 2018). Indeed, a recent publication from the Herbison lab demonstrates that glutamate generates small “synchronizing” events that are dependent on the ionotropic receptors (Han et al., 2023), but the fast neurotransmitter is unable to support the sustained firing (i.e., slow EPSP) that is necessary for peptide release and synchronization of the KNDy network. Rather, we postulate that glutamate neurotransmission is more important for excitation of Kiss1AVPV/PeN neurons and facilitating the GnRH (LH) surge with high circulating levels of E2 when peptide neurotransmitters are at a nadir and glutamate levels are high in female Kiss1ARH neurons. Indeed, low frequency (5 Hz) optogenetic stimulation of Kiss1ARH neurons, which only releases glutamate in E2-treated, ovariectomized females (Qiu J. et al., 2016), generates a surge-like increase in LH release during periods of optical stimulation (Lin et al., 2021; Voliotis et al., 2021). In a subsequent study optical stimulation of Kiss1ARH neuron terminals in the AVPV at 20 Hz, a frequency commonly used for terminal stimulation in vivo, generated a similar surge of LH (Shen et al., 2022). Additionally, intra-AVPV infusion of glutamate antagonists, AP5+CNQX, completely blocked the LH surge induced by Kiss1ARH terminal photostimulation in the AVPV (Shen et al., 2022). Therefore, there appears to be a clear role for glutamatergic transmission from the Kiss1ARH to Kiss1AVPV/PeN neurons in amplifying the LH surge in the female mouse. Finally it is important to keep in mind that even in the presence of high physiological levels of E2, the mRNA expression of Tac2 is many-fold higher than Kiss1 (Figure 10), which is essential for NKB maintaining synchronous firing of Kiss1ARH neurons, albeit at a lower frequency, across all physiological states (Qiu J. et al., 2016). Indeed, there is a progressive change from a strictly pulsatile pattern of GnRH in the hypophyseal portal circulation to one containing both pulsatile and non-pulsatile components during the development of the GnRH surge in the ewe (Evans, Dahl, Mauger, & Karsch, 1995; Evans, Dahl, Mauger, Padmanabhan, et al., 1995), and a pulsatile mode of LH secretion during the preovulatory LH surge is also evident in other species including humans (Rossmanith et al., 1990). Therefore, we believe that our cellular molecular and electrophysiological findings in combination with our computational modelling provide a foundation for understanding the complex role of Kiss1ARH neurons in controlling fertility in the mammal. Finally, our model provides the first comprehensive biophysical description of the conductances underlying the neuronal activity of KissARH neurons, which can serve as a basis for future computational modelling of the KissARH neuronal network and its interactions with other brain regions involved in the complex regulation of mammalian female reproduction.
Methods and Materials
Animals
All the animal procedures described in this study were performed in accordance with institutional guidelines based on National Institutes of Health standards and approved by the Institutional Animal Care and Use Committee at Oregon Health and Science University (OHSU).
Mice
Kiss1Cre transgenic female mice version 2 (Padilla et al., 2018) were selectively bred at OHSU. They also were crossed with heterozygous Ai32 mice (RRID:IMSR_JAX:024109, C57BL/6 background) which carry ChR2 (H134R)–EYFP gene in their Gt(ROSA)26Sor locus (Madisen et al., 2012). All animals were maintained under controlled temperature and photoperiod (lights on at 0600h and off at 1800h) and given free access to food (Lab Diets 5L0D) and water. Where specified, Kiss1Cre mice received viral injections to express channelrhodopsin 2 (ChR2) in Kiss1ARH neurons, fourteen to twenty-one days prior to each experiment as described (Qiu J. et al., 2016). Some of the females were ovariectomized seven days prior to an experiment.
Each animal was injected on day 5 following OVX with 0.25 μg E2 or vehicle, followed on day 6 with 1.50 μg E2 or vehicle and used for experiments on day 7 (Bosch et al., 2013).
AAV delivery to Kiss1Cre mice
Fourteen to twenty-one days prior to each experiment, the Kiss1Cre mice (>60 d old) received bilateral ARH injections of a Cre-dependent adeno-associated viral (AAV; serotype 1) vector encoding mCherry (AAV1-Ef1α-DIO-mCherry) or AAV1 vectors designed to encode SaCas9 and single-guide RNAs (sgRNAs) (See the SaCas9 section for specifics on the sgRNA design). Using aseptic techniques, anesthetized female mice (1.5% isoflurane/O2) received a medial skin incision to expose the surface of the skull. The glass pipette with a beveled tip (diameter = 45 μm) was filled with mineral oil, loaded with an aliquot of AAV using a Nanoject II (Drummond Scientific). ARH injection coordinates were anteroposterior (AP): −1.20 mm, mediolateral (ML): ± 0.30 mm, dorsoventral (DV): −5.80 mm (surface of brain z = 0.0 mm); 500 nl of the AAV (2.0×1012 particles/ ml) was injected (100 nl/min) into each position, and the pipette left in place for 10 min post-injection, then slowly retracted from the brain. The skin incision was closed using Vetbond (3M) and each mouse received analgesia (Rimadyl, 4-5 mg/kg, s.c.).
Generation of AAV1-FLEX-SaCas9-U6-sgTrpc5
The generation of AAV1-FLEX-SaCas9-U6-sgTrpc5 viruses were done at the University of Washington using published methods (Gore et al., 2013; Hunker et al., 2020). The constructs of sgRNAs for Trpc5 (AAV1-FLEX-SaCas9-U6sgTrpc5) were designed to target exon2 and exon7, respectively (Figure 11A and B). To achieve Trpc5 mutagenesis in Kiss1ARH neurons, Kiss1Cre mice were co-injected with AAV1-DIO-mCherry, AAV1-FLEX-SaCas9-U6-sgTrpc5-exon2, and AAV1-FLEX-SaCas9-U6-sgTrpc5-exon7 at a ratio of 10%, 45%, and 45%, respectively. Control animals were co-injected with AAV1-FLEX-SaCas9-U6-sgRosa26 and AAV1-DIO-mCherry at a ratio of 90% and 10%, respectively. AAV1-DIO-mCherry was co-injected with Cas9 vectors to confirm the targeting of injections and visualize the infected Kiss1ARH neurons. Each animal was used for experiments three weeks after the viral injection, with OVX performed one week prior to the experiments.
Visualized whole-cell patch recording
Electrophysiological and optogenetic studies were made in coronal brain slices (250 μm) containing the ARH from AAV1-EF1α-DIO-mCherry injected Kiss1Cre:GFP or Kiss1-Cre:EYFP::AI32 mice, which were vehicle-treated OVX, and E2-treated OVX females 10 weeks and older as previously described (Qiu J. et al., 2016; Qiu J. et al., 2018). Whole-cell patch recordings were performed in voltage-clamp and current-clamp as previously described (Qiu J. et al., 2018) using an Olympus BX51 W1 fixed stage scope out-fitted with epifluorescence and IR-DIC video microscopy. Patch pipettes (A-M Systems; 1.5 μm outer diameter borosilicate glass) were pulled on a Brown/Flaming puller (Sutter Instrument, model P-97) and filled with the following solution (unless otherwise specified): 128 mM potassium gluconate, 10 mM NaCl, 1 mM MgCl2, 11 mM EGTA, 10 mM HEPES, 2 mM ATP, and 0.25 mM GTP adjusted to pH 7.3 with KOH; 295 mOsm. Pipette resistances ranged from 3.5–4 MΩ. In whole-cell configuration, access resistance was less than 30 MΩ; the access resistance was 80% compensated. The input resistance was calculated by measuring the slope of the I-V relationship curve between −70 and −50 mV. Standard whole-cell patch recording procedures and pharmacological testing were performed as previously described (Qiu J. et al., 2003; Qiu J. et al., 2014). Electrophysiological signals were digitized with a Digidata 1322A (Axon Instruments) and the data were analyzed using p-Clamp software (Molecular Devices, Foster City, CA). The liquid junction potential was corrected for all data analysis.
For optogenetic stimulation, a light-induced response was evoked using a light-emitting diode (LED) 470 nm blue light source controlled by a variable 2A driver (ThorLabs, Newton, NJ) with the light path directly delivered through an Olympus 40x water-immersion lens. For high-frequency (20 Hz) stimulation the length of stimulation was 10 seconds (Qiu J. et al., 2016).
For studying the firing pattern of Kiss1ARH neurons from OVX and E2-treated OVX mice, the electrodes were filled with an internal solution consisting of the following (in mM): 120 KMeSO4, 5 Na2 phosphocreatine, 10 KCl, 2 MgSO4, 10 HEPES, 4 K2ATP, and 0.25 GTP adjusted to pH 7.3 with KOH; 290 mOsm. The bath solution as described (Ashhad and Feldman, 2020) consisted of (in mM): 124 NaCl, 3.5 KCl, 25 NaHCO3, 0.5 NaH2PO4 and 30 D-glucose,1 MgSO4, 1.5 CaCl2, gassed with 95% O2 /5% CO2 (pH 7.4, 310 mOsm).
For studying the activation/inactivation characteristics of the Ca2+ current, the electrodes were filled with an internal solution as described (Qiu J. et al., 2003) consisting of the following (in mM): 100 Cs+ gluconate, 20 TEA-Cl, 10 NaCl, 1 MgCl2, 10 HEPES, 11 EGTA, 1 ATP, 0.25 GTP, the pH was adjusted to 7.3 with CsOH at 300 mOsm. The bath solution as described (Zhang X. B. and Spergel, 2012) consisted of (in mM) 117.5 NaCl, 25 NaHCO3, 1.25 NaH2PO4, 10 TEA-Cl, 2 CaCl2, 1 MgCl2, 20 sucrose and 5 glucose, gassed with 95% O2 /5% CO2 (pH 7.4, 309 mOsm) and supplemented with 1 μM TTX, 10 μM CNQX, 50 μM AP5 and 100 μM picrotoxin. The effects of estrogen treatment on the peak current, peak current density, and activation/inactivation characteristics of calcium current were measured. Activation curves were fitted by the Boltzmann equation: I/Imax=1/{1+exp[V1/2-Vs)/k]}, where I is the peak current at the step potential Vs, Imax is the peak current amplitude, V1/2 is the step potential yielding half-maximum current, and k is the slope factor. Inactivation curves were fit with the Boltzmann equation: I/Imax = 1 - 1/ {1 + exp [(VH - V1/2)/k]}, where I is the peak current at the step potential VH, Imax is the peak current amplitude, V1/2 is the step potential at which half the current is inactivated, and k is the slope factor.
To record M-currents, pipettes were filled with an internal solution consisting of 10 mM NaCl, 128 mM K-gluconate, 1 mM MgCl, 10 mM HEPES, 1 mM ATP, 1.1 mM EGTA, and 0.25 mM GTP (pH 7.3; 290 mOsm). During voltage-clamp, we employed a standard deactivation protocol (Conde and Roepke, 2019; Roepke et al., 2011) to measure potassium currents. This involved 500-ms voltage steps ranging from –30 to –75 mV in 5-mV increments, following a 300-millisecond prepulse to –20 mV. The amplitude of the M-current relaxation or deactivation was quantified as the difference between the initial (<10 ms) and sustained current (>475 ms) of the current trace.
The bath solution for whole-cell recording of BK, SK and M currents was aCSF supplemented with 1 μM TTX, 10 μM CNQX, 50 μM AP5 and 100 μM picrotoxin.
Electrophysiological solutions/drugs
A standard artificial cerebrospinal fluid (aCSF) was used (Qiu J. et al., 2003; Qiu J. et al., 2010). All drugs were purchased from Tocris Bioscience unless otherwise specified. 1 mM TTX (Alomone Labs), 50 mM DL-2-Amino-5-phosphonopentanoic acid sodium salt (AP5), 10 mM 6-Cyano-7-nitroquinoxaline-2,3-dione disodium (CNQX), 100 mM picrotoxin, 10 mM Nifedipine (Sigma), 2 mM ω-conotoxin GVIA (ConoGVIA; Alomone Labs), 0.5 mM ω-conotoxin MVIIC (ConoMVIIC; Alomone Labs, 50 μM ω-agatoxin IVA (AgaIVA; Alomone Labs), 50 μM SNX-482 (Alomone Labs), 10 mM TTA-P2 (TTAP2; Alomone Labs), 200 mM CdCl2 (Cd2+; Sigma), 40 mM XE991 (Alomone Labs), 100 μM Iberiotoxin (Alomone Labs), 500 μM Apamin (Alomone Labs), 100 mM glutamic acid and 1 mM senktide. Stocks (1000×) were prepared in dimethylsulfoxide DMSO (picrotoxin, TTAP2 and senktide) or water (TTX, AP5, CNQX, ConoGVIA, ConoMVIIC, AgaIVA, SNX-482, CdCl2 and glutamic acid) and stored at −20°C. Aliquots of the stock solutions were stored as appropriate until needed.
Cell harvesting of dispersed Kiss1Cre neurons and real-time quantitative PCR (qPCR)
Cell harvesting and qPCR was conducted as previously described (Bosch et al., 2013). The ARH was microdissected from basal hypothalamic coronal slices obtained from female Kiss1Cre version 2 mice (Padilla et al., 2018) (n = 4-6 animals/group). The dispersed cells were visualized, patched, and then harvested (10 cells/tube) as described previously (Bosch et al., 2013). Briefly, ARH tissue was incubated in papain (7 mg/ml in oxygenated aCSF) for 50 min at 37° C then washed 4 times in low Ca2+ aCSF and two times in aCSF. For cell dispersion, Pasteur pipettes were flame polished to decreasing tip sizes and gentle trituration used to disperse the neurons onto a glass bottom dish. The plated cells were bathed in oxygenated aCSF using a peristaltic pump to keep the cells viable and clear of debris. Healthy cells with processes and a smooth cell membrane were harvested. Pipettes (World Precision Instruments; 1.5 μm outer diameter borosilicate glass) were pulled on a Brown/Flaming puller (Sutter Instrument, model P-87) to a 10 µm diameter tip. The cells were harvested using the XenoWorks Microinjector System (Sutter Instruments, Navato, CA) which provides negative pressure in the pipette and fine control to draw the cell up into the pipette. Cell pools were harvested and stored at -80°C. All cell pools were DNAse-treated using DNase1. cDNA synthesis was performed as previously described (Bosch et al., 2013).
Primers for the genes that encode for low and high voltage-gated calcium channels, TRPC5, Vglut2, large conductance calcium-activated K+ (BKα) channels, small conductance calcium-activated K+ (SK3) channels and GAPDH were designed using Clone Manager software (Sci Ed Software) to cross at least one intron-exon boundary and optimized as previously described using Power Sybr Green method (Bosch et al., 2013). Real-time qPCR controls included neuronal pools without reverse transcriptase (-RT), hypothalamic RNA with RT (+) and without RT (-), as well as water blanks. Standard curves using ARH cDNA were utilized to determine the real-time PCR efficiency (E = 10(−1/m) – 1) (Biosystems, 2006; Pfaffl, 2001). Only primers resulting in efficiencies of 90-100% were used for analysis. Primer sequences, qPCR parameters and efficiency calculations are provided in Table 1.
mRNA expression analysis
qPCR was performed on a Quantstudio 7 Flex Real-Time PCR System (Applied Biosystems) using Power SYBR Green PCR Master Mix (Applied Biosystems) according to established protocols (Bosch et al., 2013). The comparative ΔΔCT method (Livak and Schmittgen, 2001; Pfaffl, 2001; Schmittgen and Livak, 2008) was used to determine values from duplicate samples of 4 µl for the target genes and 2 µl for the reference gene GAPDH in a 20 μl reaction volume containing 1x Power SYBR Green PCR Master Mix and 0.5 μM forward and reverse primers. Three to six 10-cell pools per animal were analyzed from 4-6 animals per group. The relative linear quantity was determined using the 2-ΔΔCT equation (Livak and Schmittgen, 2001; Pfaffl, 2001; Schmittgen and Livak, 2008). Relative mRNA expression level of target genes in Kiss1Cre neurons was obtained by comparing OVX Oil-treated controls to OVX E2-treated animals. The mean Δ CT for the target genes from the OVX Oil-treated control samples was used as the calibrator. The data were expressed as n-fold change in gene expression normalized to the reference gene GAPDH and relative to the calibrator.
Experimental design and Statistical analysis
For the whole-cell patch recording, pharmacological experiments, only one cell was recorded per slice. Two to three slices were analyzed from each Kiss1Cre mouse, with at least 3-5 mice contributing to each group. For cell harvesting of dispersed Kiss1Cre-YFP neurons and qPCR measurements, 10 cells per pool and 3-6 pools from each animal were used, unless otherwise specified. Statistical comparisons between two groups were performed using an unpaired two-tailed Student’s t-test. Comparisons between more than two groups were performed using the repeated measures, multifactorial ANOVA. If a significant interaction was encountered, we then moved to the one-way ANOVA, followed by the multiple range tests as specified in the appropriate figure legends. All data were analyzed using GraphPad Prism version 6. All data are presented as mean ± standard error of the mean (SEM). Differences were considered statistically significant if the probability of error was less than 5%.
Mathematical model and simulation of neuronal behavior
A mathematical model of the Kiss1ARH neuron was developed and calibrated based on our physiological findings. Here, we employed the Hodgkin-Huxley modelling approach (Hodgkin and Huxley, 1952). Accordingly, the equation describing the membrane potential Vm of the Kiss1ARH neuron is given by
where Cm is the membrane capacitance and I is the sum of 12 ionic currents:
INaT and INaP are the transient and persistent sodium currents, respectively; represents the A current; represents the M current; ISK and IBK are the potassium currents through the SK and BK channels respectively; Ih the HCN current; IT the T-type calcium current; ICa represents other calcium currents (L-, N-, P/Q-, and R-type); ITRPCS represents Calcium current throught the TRPC5 channel; IGIRK potassium current through the GIRK channels. Finally, Ileak represents the contribution of leak currents. We use the Hodgkin-Huxley formalism to model current dynamics and their dependence on the membrane potential. A complete specification of all currents along with the complete table of model parameters can be found in the supplementary material. Model simulations were conducted in Matlab using the built-in numerical solver (ode45; based on an explicit Runge-Kutta (4, 5) formula). The Matlab code is available at https://git.exeter.ac.uk/mv286/kiss1-arcuate-neuron-model.
Acknowledgements
Funding
All electrophysiology and molecular biological studies were funded by the National Institutes of Health Grant R01-DK68098 (OKR and MJK, multi-PI). The generation of the sgRNA’s was funded by National Institutes of Health Grants P30-MH048736 and R01104450 (LSZ). The computational modeling was funded by the BBSRC via grant BB/W0058831/1 (KTA and MV) and the EPSRC via grant EP/T017856/1 (KTA). We thank Zoe Plain for her contributions to the modeling. The BBSRC also provided an International Partnership Award, BB/3019978/1, to facilitate collaboration between the UK partners (KOB, XFLi, KTA and MV) and the USA partners (MJK, OKR, JQ and MAB).
Supplemental Information
A mathematical model of the arcuate nucleus kisspeptin neuron
A schematic diagram of the Arcuate nucleus Kiss1 (Kiss1ARH) neuron model is presented in Fig. S1 and parameter values used in the simulations are given in Table S1.
The equation describing the membrane potential, Vm, of Kiss1ARH neurons is given by
where Cm is the membrane capacitance and is the sum of 12 ionic currents:
INaT and INaP are the transient and persistent sodium currents, respectively;IA represents the A current;IM represents the M current;ISK and IBK are the potassium currents through the SK and BK channels respectively;IT the T-type calcium current; represents other calcium currents (L-, N-, P/Q-, and R-type);ITRPC5 represents Calcium current through the TRPC5 channel; potassium current through the GIRK channels. Finally,Ileak represents the contribution of leak currents.
We use the Hodgkin-Huxley formalism to model current dynamics and their dependence on the membrane potential. Below we detail are the equations governing the currents.
Transient sodium current
where gNaT is the maximum conductance;ENa is the sodium reversal potential; hNaT is the inactivation gating variable that obeys the following equation:
Parameter τh,NaT dictates the timescale of inactivation and hNaT,∞ (Vm) is the steady-state inactivation function:
Parameter Vh,NaT describes the voltage achieving half-maximal inactivation and parameter kh,NaT, is the associated scaling function.
Finally, in the current formulation mNaT,∞ (Vm) is the steady-state activation function given by:
The transient sodium channel is modelled using parameter values from the Purkinje neuron [1]. This neuron was chosen as a baseline as it contains the same subunits, i.e., NaV1.1-α, NaV1.2-α, and NaV1.6-α [1], as the transient sodium channel in Kiss1 neuron [2].
Persistent sodium current
gNaP is the maximum conductance; hNaP is the corresponding inactivation gating variable that obeys the following equation:
and the steady-state activation and inactivation functions are given by:
The above description of the persistent sodium current was taken from a model of the GnRH neuron [3].
A current
gNaP denotes the maximum conductance;EK is the potassium reversal potential; and mA and hA are the corresponding activation and inactivation gating variables, which are described by the following equations:
The steady-state activation and inactivation functions are given by:
The model of the A-current was based on Mendonca’s model of Kv4 channels [4], as these channels are also found in Kiss1ARH neurons [5].
BK current
Here, gBK is the maximum conductance; and bBK,∞ (Vm,C) is the steady-state activation function that depends on the membrane potential, Vm, as well as on the cytosolic calcium concertation C :
The model of the BK-current was based on the model presented in [6], with the conductance parameter fitted to the current-voltage relationships recorded from Kiss1ARH neurons in the absence and presence of the specific BK blocker, iberiotoxin.
SK current
gSK denotes the maximum conductance; and bSK,∞ (C) is the steady-state activation function, which depends on the cytosolic calcium concertation, C:
The model of the BK-current was based on the model presented in [7], with the conductance fitted to the current-voltage relationships recorded from Kiss1ARH neurons in the absence and presence of the specific SK blocker, apamin.
M current
gM denotes the maximum conductance, and mM is the corresponding activation gating variable:
with the steady-state activation function, mM,∞(Vm) taking the form:
The model of the M-current was parameterised using the steady-state voltage-clamp measurements from actuate Kiss1 neurons [8], while for the activation timescale we used the timescale used in a model of the CA1/3 pyramidal cells [9].
h current
gh denotes the maximum conductance; and mh,1 and mh,2 are separate activation gating variables operating on different timescales τm,h,1 and τm,h,2 respectively):
The corresponding steady-state activation functions are:
Finally, parameter ph dictates the relative contribution of and mh,1 and mh,2 to the total current. This model of the h-current is based on the hippocampal CA1 pyramidal neuron [7].
T-type calcium current
gT is the maximum conductance; and hT,1 and hT,2 are separate inactivation gating variables operating on different timescales (τh,T,1 and τh,T,1 respectively):
The corresponding steady-state inactivation functions are:
The steady-state activation function is given by:
Finally, parameter pT dictates the relative contribution of and hT,1 and hT,2 to the total current. To model of the T-current was based on AVPV kisspeptin neurons data presented in [2, 10].
L-, N-, P/Q-, R-type calcium currents
gCa denotes the maximum conductance; and mCa and hCa are the corresponding activation and inactivation gating variables, which are described by the following equations:
The steady-state activation and inactivation functions are given by:
Parameters of the model for the high voltage activated calcium channels were calibrated from the current-voltage relationships obtained from Kiss1ARH neurons (see Figure 2 and 6 of the main text).
TRPC5 current
gTRPC5, denotes the maximum conductance; and bTRPC5 (C, RTRPC5,act)the activating gating variable that depends both on cytosolic calcium concertation (c) and on NKB-mediated activation of an intermediary effector, RTRPC5,act [11]:
The dynamics of RTRPC5,act (activated form of RTRPC5) are described by:
where NKB is the extracellular NKB concertation; kR,0 is the basal rate of RTRPC5 activation; kR is the maximal rate of RTRPC5 activation in the presence of NKB; k−R is the rate of RTRPC5 inactivation; and RTRPC5,T is the total concentration of the effector.
GIRK current
gGIRK denotes the maximum conductance; and bGIRK (Vm, RGIRK,actB)the activating gating variable that depends on membrane potential, Vm, and on external activation of an intermediary effector, RGIRK,act:
The dynamics of mGIRK are described by:
where the steady state activation function and timescale function are given by:
The dynamics of RGIRK,act are described by:
where s is the extracellular concertation of the activation signal.
The model and parameters of the GIRK current is taken from [12].
Leak currents
The leak current parameters were calibrated to current-voltage relationships recorded from Kiss1ARH neurons in the absence/presence of iberiotoxin (BK blocker) and apamin (SK blocker).
Intracellular calcium dynamics
Finally, the intracellular calcium dynamics are described via the following equation:
where parameter γ converts the currents to molecule fluxes and parameter dca dictates the linear rate at which calcium is depleted or pumped out of the cell.
Modelling the effect of E2
The effect of E2 on ionic currents is modelled as a change in the maximum conductance parameter. For currents IM, IT, ICa and ITRPC5 this change is inferred from the qPCR data assuming that the conductance is directly proportional to the mRNA expression. For ISK, IBK,Ileak,the OVX and OVX+E2 conductances are obtained from current-voltage relationships recorded from Kiss1 neurons in the absence/presence of iberiotoxin (BK blocker) and apamin (SK blocker). All other currents were assumed to be unaffected by E2. Parameter Cm in calibrated using direct measurements of the membrane conductance in the OVX and OVX+E2 state. For the simulation presented in the Figure 13 of the main text, conductances were varied within ranges that captures the physiological effect of E2.
Computer simulations
Integration of the differential equations describing the model was carried out in MATLAB R2023b using a standard 4th order Runge-Kutta method. Parameter fitting was also conducted in MATLAB R2023b using the least squares curve fitting method.
References
- Organization and function of TRPC channelsomesPflügers Archiv: European Journal of Physiology 455:187–200https://doi.org/10.1007/s00424-007-0252-0
- The calcium-activated slow AHP: cutting through the Gordian knotFrontiers in Cellular Neuroscience 6https://doi.org/10.3389/fncel.2012.00047
- Involvement of TRPC4 and 5 Channels in Persistent Firing in Hippocampal CA1 Pyramidal CellsCells 9https://doi.org/10.3390/cells9020365
- Emergent elements of inspiratory rhythmogenesis: network synchronization and synchrony propagationNeuron 106:482–497
- Mechanisms of oscillatory activity in guinea-pig nucleus reticularis thalami in vitro: a mammalian pacemakerThe Journal of Physiology 468:669–691https://doi.org/10.1113/jphysiol.1993.sp019794
- Transient receptor potential-like channels mediate metabotropic glutamate receptor EPSCs in rat dopamine neuronesThe Journal of Physiology 555:323–330https://doi.org/10.1113/jphysiol.2003.060061
- TrpC3/C7 and Slo2.1 are molecular targets for metabotropic glutamate receptor signaling in rat striatal cholinergic interneuronsThe Journal of Neuroscience 27:8845–8856https://doi.org/10.1523/JNEUROSCI.0551-07.2007
- Applied Biosystems Application Note TaqMan Gene Expression AssaysUnited States: Applied Biosystems
- The TRPC class of ion channels: a critical review of their roles in slow, sustained increases in intracellular Ca2+ concentrationsAnnual Review of Pharmacology and Toxicology 49:395–426https://doi.org/10.1146/annurev.pharmtox.48.113006.094928
- Intracellular calcium strongly potentiates agonist-activated TRPC5 channelsThe Journal of General Physiology 133:525–546https://doi.org/10.1085/jgp.200810153
- Calcium-activated potassium channelsTrends in Neurosciences 10:463–467
- Trpc5 deficiency causes hypoprolactinemia and altered function of oscillatory dopamine neurons in the arcuate nucleusProceedings of the National Academy of Sciences 116:15236–15243https://doi.org/10.1073/pnas.1905705116
- Small conductance Ca2+-activated K+ channel knock-out mice reveal the identity of calcium-dependent afterhyperpolarization currentsThe Journal of Neuroscience 24:5301–5306https://doi.org/10.1523/JNEUROSCI.0182-04.2004
- Small-conductance calcium-activated potassium channelsAnnals of the New York Academy of Sciences 868:370–378
- SK channels in excitability, pacemaking and synaptic integrationCurrent Opinion in Neurobiology 15:305–311https://doi.org/10.1016/j.conb.2005.05.001
- Distribution, neuronal co-localization and 17β-E2 modulation of small conductance calcium-activated K+ channel (SK3) mRNA in the guinea pig brainEndocrinology 143:1097–1107
- mRNA expression of ion channels in GnRH neurons: subtype-specific regulation by 17β-EstradiolMolecular and Cellular Endocrinology 367:85–97https://doi.org/10.1016/j.mce.2012.12.021
- Characterisation of K+ Channels in Human Fetoplacental Vascular Smooth Muscle CellsPLOS One 8https://doi.org/10.1371/journal.pone.0057451
- Neural KCNQ (Kv7) channelsBritish Journal of Pharmacology 156:1185–1195https://doi.org/10.1111/j.1476-5381.2009.00111.x
- TRP channels as cellular sensorsNature 426:517–524https://doi.org/10.1038/nature02196
- 17β-estradiol increases arcuate KNDy neuronal sensitivity to ghrelin inhibition of the M-current in female miceNeuroendocrinology https://doi.org/10.1159/000503146
- Kisspeptin can stimulate gonadotropin-releasing hormone (GnRH) release by a direct action at GnRH nerve terminalsEndocrinology 149:3926–3932
- Hypogonadotropic hypogonadism due to loss of function of the KiSS 1-derived peptide receptor GPR54Proceedings of the National Academy of Sciences of the United States of America 100:10972–10976https://doi.org/10.1073/pnas.1834399100
- Electrophysiology of guinea-pig supraoptic neurones: role of a hyperpolarization-activated cation current in phasic firingThe Journal of Physiology 460:407–425
- Role of a T-type calcium current in supporting a depolarizing potential, damped oscillations, and phasic firing in vasopressinergic guinea pig supraoptic neuronsNeuroendocrinology 57:789–800https://doi.org/10.1159/000126438
- Estradiol induces both qualitative and quantitative changes in the pattern of gonadotropin-releasing hormone secretion during the presurge period in the eweEndocrinology 136:1603–1609
- Does estradiol induce the preovulatory gonadotropin-releasing hormone (GnRH) surge in the ewe by inducing a progressive change in the mode of operation of the GnRH neurosecretory systemEndocrinology 136:5511–5519
- Synaptic activation of transient receptor potential channels by metabotropic glutamate receptors in the lateral amygdalaNeuroscience 137:781–794https://doi.org/10.1016/j.neuroscience.2005.09.027
- Kisspeptin neurons in the arcuate nucleus of the ewe express both dynorphin A and neurokinin BEndocrinology 148:5752–5760https://doi.org/10.1210/en.2007-0961
- Manipulating gene expression in projection-specific neuronal populations using combinatorial viral approachesCurrent Protocols in Neuroscience 65:1–20https://doi.org/10.1002/0471142301.ns0435s65
- Molecular properties of Kiss1 neurons in the arcuate nucleus of the mouseEndocrinology 152:4298–4309https://doi.org/10.1210/en.2011-1521
- XE991 and linopirdine are state-dependent inhibitors for Kv7/KCNQ channels that favor activated single subunitsJournal of Pharmacology and Experimental Therapeutics 362:177–185https://doi.org/10.1124/jpet.117.241679
- BK potassium channels facilitate high-frequency firing and cause early spike frequency adaptation in rat CA1 hippocampal pyramidal cellsThe Journal of Physiology 580:859–882https://doi.org/10.1113/jphysiol.2006.126367
- Mechanism of kisspeptin neuron synchronization for pulsatile hormone secretion in male miceCell Reports 42https://doi.org/10.1016/j.celrep.2022.111914
- Rat GnRH neurons exhibit large conductance voltage- and Ca2+-Activated K+ (BK) currents and express BK channel mRNAsThe Journal of Physiological Sciences 58:21–29https://doi.org/10.2170/physiolsci.RP013207
- A quantitative description of membrane current and its application to conduction and excitation in nerveThe Journal of Physiology 117:500–544
- Conditional single vector CRISPR/SaCas9 viruses for efficient mutagenesis in the adult mouse nervous systemCell Reports 30:4303–4316https://doi.org/10.1016/j.celrep.2020.02.092
- Characterization of voltage-gated calcium currents in gonadotropin-releasing hormone neurons tagged with green fluorescent protein in ratsEndocrinology 144:5118–5125https://doi.org/10.1210/en.2003-0213
- TRPCing around the hypothalamusFrontiers in Neuroendocrinology 51:116–124https://doi.org/10.1016/j.yfrne.2018.05.004
- Electrophysiological analysis of neuroendocrine neuronal activity in hypothalamic slicesMethods in Neurosciences: Pulsatility in Neuroendocrine Systems San Diego: Academic Press, Inc :47–67
- The metastasis suppressor gene KiSS-1 encodes kisspeptins, the natural ligands of the orphan G protein-coupled receptor GPR54Journal of Biological Chemistry 276:34631–34636https://doi.org/10.1074/jbc.M104847200
- Physiological Roles and Therapeutic Potential of Ca2+ Activated Potassium Channels in the Nervous SystemFrontiers in Molecular Neuroscience 11https://doi.org/10.3389/fnmol.2018.00258
- Properties of two calcium-activated hyperpolarizations in rat hippocampal neuronesThe Journal of Physiology 389:187–203
- Molecular basis of R-type calcium channels in central amygdala neurons of the mouseProceedings of the National Academy of Sciences 99:3276–3281https://doi.org/10.1073/pnas.052697799
- Neuroanatomy of the kisspeptin signaling system in mammals: comparative and developmental aspectsAdvances in Experimental Medicine and Biology 784:27–62https://doi.org/10.1007/978-1-4614-6199-9_3
- Anatomy of the kisspeptin neural network in mammalsBrain Research 1364:90–102https://doi.org/10.1016/j.brainres.2010.09.020
- Optogenetic Activation of Arcuate Kisspeptin Neurons Generates a Luteinizing Hormone Surge-Like Secretion in an Estradiol-Dependent MannerFront Endocrinol (Lausanne) 12https://doi.org/10.3389/fendo.2021.775233
- Highly redundant neuropeptide volume co-transmission underlying episodic activation of the GnRH neuron dendroneLife 10https://doi.org/10.7554/eLife.62455
- Analysis of relative gene expression data using real-time quantitative PCR and the 2™ΔΔCT methodMethods 25:402–408https://doi.org/10.1006/meth.2001.1262
- The intrinsic electrophysiological properties of mammalian neurons: insights into central nervous system functionScience 242:1654–1664
- H-current: properties of a neuronal and network pacemakerNeuron 21:9–12https://doi.org/10.1016/S0896-6273(00)80509-7
- Prolactin regulates tuberoinfundibular dopamine neuron discharge pattern: novel feedback control mechanisms in the lactotrophic axisThe Journal of Neuroscience 32:8074–8083https://doi.org/10.1523/JNEUROSCI.0129-12.2012
- Synchronized network oscillations in rat tuberoinfundibular dopamine neurons: swith to tonic discharge by thyrotropin-releasing hormoneNeuron 65:217–229
- A toolbox of Cre-dependent optogenetic transgenic mice for light-induced activation and silencingNature Neuroscience 15:793–802https://doi.org/10.1038/nn.3078
- Selective activation of Ca2+-activated K+ channels by co-localize Ca2+ channels in hippocampal neuronsNature 395:900–905
- The physiological role of calcium-dependent channelsTrends in Neurosciences 12:420–424
- Diurnal properties of voltage-gated Ca2+ currents in suprachiasmatic nucleus and roles in action potential firingThe Journal of Physiology 598:1775–1790https://doi.org/10.1113/JP278327
- Kv4.2 channel activity controls intrinsic firing dynamics of arcuate kisspeptin neuronsThe Journal of Physiology 596:885–899https://doi.org/10.1113/JP274474
- Kisspeptin directly stimulates gonadotropin-releasing hormone release via G protein-coupled receptor 54Proceedings of the National Academy of Sciences of the United States of America 102:1761–1766https://doi.org/10.1073/pnas.0409330102
- Characterization and regulation of pre-ovulatory secretion of gonadotrophin-releasing hormoneHuman Reproduction 8:51–56
- Regulation of gonadotropin-releasing hormone secretion by kisspeptin/dynorphin/neurokinin B neurons in the arcuate nucleus of the mouseThe Journal of Neuroscience 29:11859–11866https://doi.org/10.1523/JNEUROSCI.1569-09.2009
- Regulation of NKB pathways and their roles in the control of Kiss1 neurons in the arcuate nucleus of the male mouseEndocrinology 152:4265–4275https://doi.org/10.1210/en.2011-1143
- Optogenetic stimulation of arcuate nucleus Kiss1 neurons reveals a steroid-dependent glutamatergic input to POMC and AgRP neurons in male miceMolecular Endocrinology 30:630–644https://doi.org/10.1210/me.2016-1026
- The coupling of neurotransmitter receptors to ion channels in the brainScience 241:545–551
- BK Channel Regulation of Afterpotentials and Burst Firing in Cerebellar Purkinje NeuronsThe Journal of Neuroscience 41:2854–2869https://doi.org/10.1523/jneurosci.0192-20.2021
- Calcium current subtypes in GnRH neuronsBiology of Reproduction 69:1914–1922https://doi.org/10.1095/biolreprod.103.019265
- Radiotelemetric monitoring of hypothalamic gonadotropin-releasing hormone pulse generator activity throughout the menstrual cycle of the rhesus monkeyEndocrinology 129:1207–1214https://doi.org/10.1210/endo-129-3-1207
- A neural circuit underlying the generation of hot flushesCell Reports 24:271–277https://doi.org/10.1016/j.celrep.2018.06.037
- Molecular physiology of low-voltage-activated t-type calcium channelsPhysiological Reviews 83:117–161https://doi.org/10.1152/physrev.00018.2002
- Conditional transgenic suppression of M channels in mouse brain reveals functions in neuronal excitability, resonance and behaviorNature Neuroscience 8:51–60https://doi.org/10.1038/nn1375
- A new mathematical model for relative quantification in real-time RT-PCRNucleic Acids Research 29https://doi.org/10.1093/nar/29.9.e45
- Rapid signaling of estrogen in hypothalamic neurons involves a novel G-protein-coupled estrogen receptor that activates protein kinase CThe Journal of Neuroscience 23:9529–9540https://doi.org/10.1523/JNEUROSCI.23-29-09529.2003
- Leptin excites proopiomelanocortin neurons via activation of TRPC channelsThe Journal of Neuroscience 30:1560–1565https://doi.org/10.1523/JNEUROSCI.4816-09.2010
- High-frequency stimulation-induced peptide release synchronizes arcuate kisspeptin neurons and excited GnRH neuronseLife 5https://doi.org/10.7554/eLife.16246
- Estrogenic-dependent glutamatergic neurotransmission from kisspeptin neurons governs feeding circuits in femaleseLife 7https://doi.org/10.7554/eLife.35656
- Deletion of Stim1 in hypothalamic arcuate nucleus Kiss1 neurons potentiates synchronous GCaMP activity and protects against diet-induced obesityJ Neurosc 41:9688–9701https://doi.org/10.1523/JNEUROSCI.0622-21.2021
- Insulin excites anorexigenic proopiomelanocortin neurons via activation of canonical transient receptor potential channelsCell Metabolism 19:682–693https://doi.org/10.1016/j.cmet.2014.03.004
- Menopause and the human hypothalamus: Evidence for the role of kisspeptin/neurokinin B neurons in the regulation of esstrogen negative feedbackPeptides 30:111–122
- Hypertrophy and Increased Gene Expression of Neurons Containing Neurokinin-B and Substance-P Messenger Ribonucleic Acids in the Hypothalami of Postmenopausal Women*Endocrinology 128:2239–2247https://doi.org/10.1210/endo-128-5-2239
- Fasting and 17β-estradiol differentially modulate the M-current in neuropeptide Y neuronsThe Journal of Neuroscience 17:11825–11835https://doi.org/10.1523/JNEUROSCI.1395-11.2011
- Hypothalamic Kisspeptin Neurons and the Control of HomeostasisEndocrinology 163:1–11https://doi.org/10.1210/endocr/bqab253
- Relative changes in LH pulsatility during the menstrual cycle: Using data from hypogonadal women as a reference pointClinical Endocrinology 32:647–660https://doi.org/10.1111/j.1365-2265.1990.tb00909.x
- Ca2+-activated K+ currents in neurones: Types, physiological roles and modulationTrends in Neurosciences 19:150–154
- Photolytic manipulation of [Ca2+]i reveals slow kinetics of potassium channels underlying the afterhyperpolarization in hippocampal pyramidal neuronsThe Journal of Neuroscience 19:3657–3664
- Channels underlying neuronal calcium-activated potassium currentsProgress in Neurobiology 66:345–353https://doi.org/10.1016/S0301-0082(02)00004-7
- Analyzing real-time PCR data by the comparative CT methodNature Protocols 3:1101–1108
- The GPR54 gene as a regulator of pubertyThe New England Journal of Medicine 349:1614–1627https://doi.org/10.1056/NEJMoa035322
- Increased hypothalamic GPR54 signaling: a potential mechanism for initiation of puberty in primatesProceedings of the National Academy of Sciences of the United States of America 102:2129–2134https://doi.org/10.1073/pnas.0409822102
- The role of BK-type Ca2+-dependent K+ channels in spike broadening during repetitive firing in rat hippocampal pyramidal cellsThe Journal of Physiology 521:135–146
- Optogenetic stimulation of Kiss1ARC terminals in the AVPV induces surge-like luteinizing hormone secretion via glutamate release in miceFrontiers in Endocrinology 13https://doi.org/10.3389/fendo.2022.1036235
- Regulation of Kiss1 gene expression in the brain of the female mouseEndocrinology 146:3686–3692https://doi.org/10.1210/en.2005-0488
- Calcium and small-conductance calcium-activated potassium channels in gonadotropin-releasing hormone neurons before, during, and after pubertyEndocrinology 148:2383–2390
- CRISPR knockdown of Kcnq3 attenuates the M current and increases excitability of NPY/AgRP neurons to alter energy balanceMolecular Metabolism https://doi.org/10.1016/j.molmet.2021.101218
- Action potential repolarization and a fast after-hyperpolarization in rat hippocampal pyramidal cellsThe Journal of Physiology 385:733–759
- Potassium currents in hippocampal pyramidal cellsProgress in Brain Research 83:161–187
- TRPC4 and GIRK channels underlie neuronal coding of firing patterns that reflect Gq/11–Gi/o coincidence signals of variable strengthsProceedings of the National Academy of Sciences 119https://doi.org/10.1073/pnas.2120870119
- Involvement of transient receptor potential-like channels in responses to mGluR-I activation in midbrain dopamine neuronsEuropean Journal of Neuroscience 18:2133–2145https://doi.org/10.1046/j.1460-9568.2003.02936.x
- Calcium channels: mechanisms of selectivity, permeation, and blockAnnual Review of Biophysics and Biophysical Chemistry 16:265–290https://doi.org/10.1146/annurev.bb.16.060187.001405
- The KCNQ5 potassium channel mediates a component of the afterhyperpolarization current in mouse hippocampusProceedings of the National Academy of Sciences of the United States of America 107:10232–10237https://doi.org/10.1073/pnas.1004644107
- Hippocalcin gates the calcium activation of the slow afterhyperpolarization in hippocampal pyramidal celllsNeuron 53:487–493
- Contribution of KCNQ2 and KCNQ3 to the medium and slow afterhyperpolarization currentsProceedings of the National Academy of Sciences of the United States of America 105:19974–19979https://doi.org/10.1073/pnas.0810535105
- TRP channelsAnnual Review of Biochemistry 76:387–417https://doi.org/10.1146/annurev.biochem.75.103004.142819
- SK channels and the varieties of slow after-hyperpolarizations in neuronsEuropean Journal of Neuroscience 18:3155–3166
- Modulation of pulsatile GnRH dynamics across the ovarian cycle via changes in the network excitability and basal activity of the arcuate kisspeptin networkeLife 10https://doi.org/10.7554/eLife.71252
- Neurokinin B and dynorphin A in kisspeptin neurons of the arcuate nucleus participate in generation of periodic oscillation of neural activity driving pulsatile gonadotropin-releasing hormone secretion in the goatThe Journal of Neuroscience 30:3124–3132https://doi.org/10.1523/JNEUROSCI.5848-09.2010
- Crystal structure of the mammalian GIRK2 K+ channel and gating regulation by G proteins, PIP2, and sodiumCell 147:199–208https://doi.org/10.1016/j.cell.2011.07.046
- Structure and function of cyclic nucleotide-gated channelsAnnual Review of Neuroscience 19:235–263https://doi.org/10.1146/annurev.ne.19.030196.001315
- 17β-estradiol regulation of T-type calcium channels in gonadotropin-releasing hormone neuronsThe Journal of Neuroscience 29:10552–10562https://doi.org/10.1523/JNEUROSCI.2962-09.2009
- Kisspeptin inhibits a slow afterhyperpolarization current via protein kinase C and reduces spike frequency adaptation in GnRH neuronsAmerican Journal of Physiology: Endocrinology and Metabolism 304:E1237–1244https://doi.org/10.1152/ajpendo.00058.2013
- Kisspeptin inhibits high-voltage activated Ca2+ channels in GnRH neurons via multiple Ca2+ influx and release pathwaysNeuroendocrinology 96:68–80https://doi.org/10.1159/000335985
- TRPC channels underlie cholinergic plateau potentials and persistent activity in entorhinal cortexHippocampus 21:386–397https://doi.org/10.1002/hipo.20755
- Structural insights into the lipid and ligand regulation of a human neuronal KCNQ channelNeuron 110:237–247https://doi.org/10.1016/j.neuron.2021.10.029
- Mechanisms of Persistent Activity in Cortical Circuits: Possible Neural Substrates for Working MemoryAnnual Review of Neuroscience 40:603–627https://doi.org/10.1146/annurev-neuro-070815-014006
- 1.Differentiated pattern of sodium channel expression in dissociated Purkinje neurons maintained in long-term cultureJournal of Neurochemistry 101:737–748
- 2.17β-Estradiol increases persistent Na current and excitability of AVPV/PeN Kiss1 neurons in female miceMolecular Endocrinology 29:518–527
- 3.A unified model for two modes of bursting in GnRH neuronsJournal of Computational Neuroscience 40:297–315
- 4.Stochastic and deterministic dynamics of intrinsically irregular firing in cortical inhibitory interneuronseLife 5
- 5.Kv4. 2 channel activity controls intrinsic firing dynamics of arcuate kisspeptin neuronsThe Journal of Physiology 596:885–899
- 6.Mechanism of spontaneous and receptor-controlled electrical activity in pituitary somatotrophs: experiments and theoryJournal of Neurophysiology 98:131–144
- 7.Altered intrinsic pyramidal neuron properties and pathway-specific synaptic dysfunction underlie aberrant hippocampal network function in a mouse model of tauopathyJournal of Neuroscience 13:350–63
- 8.17β-estradiol increases arcuate KNDy neuronal sensitivity to ghrelin inhibition of the M-current in female miceNeuroendocrinology 110:582–594
- 9.A unified model of CA1/3 pyramidal cells: An investigation into excitabilityProgress in Biophysics and Molecular Biology 105:34–48
- 10.Excitability and burst generation of AVPV kisspeptin neurons are regulated by the estrous cycle via multiple conductances modulated by estradiol actioneNeuro 3
- 11.Deletion of Stim1 in hypothalamic arcuate nucleus Kiss1 neurons potentiates synchronous GCaMP activity and protects against diet-induced obesityJournal of Neuroscience 41:9688–9701
- 12.TRPC4 and GIRK channels underlie neuronal coding of firing patterns that reflect Gq/11 - Gi/o coincidence signals of variable strengthsProceedings of the National Academy of Sciences 119
- 13.Small-conductance calcium-activated potassium channelsAnnals of New York Acadamy of Sciences 868:370–8
- 14.Molecular properties of Kiss1 neurons in the arcuate nucleus of the mouseEndocrinology 152:4298–4309
- 15.Estrogenic-dependent glutamatergic neurotransmission from kisspeptin neurons governs feeding circuits in femaleseLife 7
Article and author information
Author information
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Reviewed Preprint version 3:
- Version of Record published:
Copyright
© 2024, Qiu et al.
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
- views
- 390
- downloads
- 13
- citations
- 0
Views, downloads and citations are aggregated across all versions of this paper published by eLife.