The Structural and Electrophysiological Properties of Progesterone Receptor-Expressing Neurons Vary along the Anterior-Posterior Axis of the Ventromedial Hypothalamus and Undergo Local Changes across the Reproductive Cycle

Visual Abstract

hypothalamus (VMHvl) contains neurons expressing receptors for sex hormones and its function is intimately linked to female sexual receptivity. However, recent findings suggest that the VMHvl is functionally heterogeneous. Here, we used whole recordings and intracellular labeling to characterize the electrophysiological and morphologic properties of individual VMHvl neurons in naturally cycling females and report the existence of multiple electrophysiological phenotypes within the VMHvl. We found that the properties of progesterone receptor expressing (PR1) neurons, but not PR-neurons, depended systematically on the neuron's location along the anterior-posterior (AP) axis of the VMHvl and the phase within the reproductive cycle. Prominent among this, the resting membrane potential of anterior PR1 neurons decreased during the receptive phase, while the excitability of medial PR1 neurons increased during the non-receptive phase. During the receptive phase of the cycle, posterior PR1 neurons simultaneously showed an increase in dendritic complexity and a decrease in spine density. These findings reveal an extensive diversity of local rules driving structural and physiological changes in response to fluctuating levels of sex hormones, supporting the anatomic and func-

Introduction
Female reproductive physiology and behavior are under the control of the ovarian sex hormones estrogen and progesterone, whose fluctuating levels act reversibly in the female brain, organizing the activity of neural circuits to synchronize sexual behavior with reproductive capacity (Snoeren, 2018;Jennings and de Lecea, 2020). The ventrolateral division of the ventromedial hypothalamus (VMHvl) is crucial for female reproductive behavior, in particular for the display of lordosis, the female acceptance posture. Non-specific electrolytic lesions (Pfaff and Sakuma, 1979a) or ablation of genetically delineated neuronal populations (Rissman et al., 1997;Yang et al., 2013) of the VMHvl virtually abolish lordosis, while electrical stimulation at the same location enhances the probability of females displaying it (Pfaff and Sakuma, 1979b). Importantly, VMHvl neurons have rich expression of receptors for estrogen (ER) and progesterone (PR) and therefore are sensitive to the fluctuating levels of these sex hormones across the reproductive cycle (Snoeren, 2018;Jennings and de Lecea, 2020). In accordance, local infusion of sex hormones in the VMHvl increases female receptivity (Rubin and Barfield, 1983a), while male-evoked activity in this hypothalamic nucleus is enhanced when females are sexually receptive (Nomoto and Lima, 2015). Moreover, the output of PR-expressing neurons of the VMHvl is altered across the reproductive cycle and this cyclic remodeling is fundamental for the expression of acceptance when females are sexual receptive (Inoue et al., 2019). In summary, vast evidence highlights the importance of the VMHvl for coordinating female's reproductive state with sexual receptivity.
Several studies have shown a large diversity of neuronal types with a wide variety of cellular identities within the VMHvl (McClellan et al., 2006;Flanagan-Cato, 2011;Kim et al., 2019). In addition, the generation of Cre lines providing access to ER-expressing and PR-expressing neurons uncovered further anatomic and functional subdivisions across the nucleus' anterior-posterior axis (AP axis;McClellan et al., 2006;Flanagan-Cato, 2011;Hashikawa et al., 2017;Inoue et al., 2019;Kim et al., 2019;Lo et al., 2019;Wang et al., 2019), revealing a much broader role of this hypothalamic region in the control of different aspects of female socio-sexual behavior. While regulation of female sexual receptivity seems to be localized to the most posterior-lateral part of the VMHvl, its posterior-medial division is involved in aggressive behavior toward intruders in mothers (Hashikawa et al., 2017). Moreover, it was recently shown that ER-expressing neurons in the anterior portion of the VMHvl are important for self-defense in males (Wang et al., 2019). The connectivity of ER-expressing neurons also varies across the AP axis, further supporting the idea of topographic heterogeneity within this small nucleus .
How fluctuating levels of sex hormones affect the properties and function of VMHvl neurons remains poorly understood. Most studies aimed at elucidating the impact of sex hormones on neuronal activity in this region have focused on the structural properties of VMHvl neurons. The majority of neurons in the VMHvl have a long primary dendrite (LPD) that extends outside the nucleus, receiving inputs from other forebrain regions, a short primary dendrite (SPD), and secondary dendrites that may integrate local inputs (Millhouse, 1973;Calizo and Flanagan-Cato, 2000). Previous studies in rodents have shown that externally primed estradiol treatment shortens LPDs and reduces dendritic density in the VMHvl, effects which are reversed by progesterone (Griffin and Flanagan-Cato, 2008;Griffin et al., 2010). In addition, externally applied ovarian hormones increase the density of dendritic spines on the SPD of VMHvl neurons (Calizo and Flanagan-Cato, 2000), and decrease the spine density on the LPD of VMHvl neurons (Calizo and Flanagan-Cato, 2002). Moreover, treatment with estrogen was found to lead to an increase in the size of neuronal somata in the VMHvl (Carrer and Aoki, 1982;Griffin and Flanagan-Cato, 2011). More recently, it was shown that elevation in the concentration of circulating estrogen increases the number of axon terminals of PR-expressing VMHvl neurons that project to the anteroventral periventricular nucleus of the hypothalamus in female mice and that this plasticity is necessary for sexual receptivity (Inoue et al., 2019). Interestingly, aggressive behavior in males correlates with androgen levels in circulation and is dependent on the physiological and structural plasticity of ER-expressing neurons in the VMHvl (Stagkourakis et al., 2020). However, the electrophysiological properties of VMHvl neurons across the reproductive cycle have received very little attention (but see Kow and Pfaff, 1985;Booth et al., 2010).
A noteworthy common factor among the vast majority of previous studies on the effects of sex hormones on VMHvl neurons and behavior is that they have been performed in ovariectomized females in which sex hormones were systemically administered at concentrations that do not always match the physiological levels that are observed in naturally cycling females (Liu and Shi, 2015). Moreover, despite the growing evidence suggesting anatomic and functional diversity of the VMHvl, most studies have been largely focused on its most posterior levels (Yang et al., 2013;Hashikawa et al., 2017;Inoue et al., 2019). However, the VMHvl extends slightly more than a millimeter across the base of the hypothalamus Lo et al., 2019), and sex hormone receptor expressing neurons are homogeneously distributed throughout (Sá and Fonseca, 2017;Kim et al., 2019). If and how the structural and physiological properties of VMHvl neurons along its AP axis vary across the female reproductive cycle remains elusive.
In order to gain a more complete understanding of the function of the VMHvl, in the present study we investigated the structural and electrophysiological properties of genetically delineated neurons along the VMHvl AP axis and across the reproductive cycle of naturally cycling females. We focused our efforts on neurons expressing the progesterone receptor (PR), as they are fundamental for female sexual behavior (Yang et al., 2013). To do so, we obtained slices from receptive and non-receptive females and performed in vitro whole-cell recordings of PR expressing (PR1) and PR non-expressing (PR-) neurons, with subsequent reconstruction of the recorded neurons to characterize their electrophysiological and structural properties.
Here, we report the existence of multiple electrophysiological phenotypes within the VMHvl and a wide variety of structural and physiological properties across the AP axis of the VMHvl, which are specific for PR1 neurons and that vary locally across the reproductive cycle. For instance, the membrane resting potential of anterior PR1 neurons decreases during the receptive phase, while the excitability of medial PR1 increases when females are non-receptive. During the receptive phase of the cycle, posterior PR1 neurons simultaneously undergo an increase in dendritic complexity and a decrease in spine density. These findings reveal an extensive diversity of local rules driving structural and physiological changes in response to fluctuating levels of sex hormones, supporting the anatomic and functional subdivision of the VMHvl and its possible role in the orchestration of different aspects of female socio-sexual behavior.
Animals were kept under controlled temperature of 23 6 1°C and photoperiod of reversed 12/12 h light/ dark cycle (light available from 8 P.M. to 8 A.M.) conditions and group-housed in standard cages with environmental enrichment elements. Food and water were provided ad libitum. Females were weaned at 20-21 d of age and group-housed with two to five animals.
After reaching six weeks of age, females were exposed to adult C57BL/6 male soiled bedding once per week to stimulate the natural reproductive cycle. All procedures were conducted in accordance with the animal protocols approved by the Portuguese National Authority for Animal Health (Direcção Geral de Alimentação e Veterinária; DGAV) and the Commission for Experimentation and Animal Welfare of the Champalimaud Centre for the Unknown (Órgão para o Bem Estar Animal; ORBEA).

Reproductive cycle monitoring
To assess the reproductive state of naturally cycling female mice, vaginal cytology samples were collected every morning for at least one entire cycle (four/five consecutive days). This collection was done by vaginal lavage using a pipette (T-210-Y, AXYGEN, with its tip cut): 10 ml of 0.01 M PBS was gently flushed into the vagina and back out three or four times without touching the vaginal wall to avoid cervical stimulation and pseudopregnancy. The flush containing vaginal fluid was transferred to a glass slide and dried (Caligioni, 2009). Papanicolaou staining was used to differentiate cells (Bio-Optica protocol) and they were observed under a Zeiss AxioScope A1 brightfield microscope with a 10Â objective. The identification of the reproductive cycle phase was done based on the proportion of each cell type in the smear. Proestrus is characterized by nucleated epithelial cells. In estrus the main cells present in the vaginal secretion are cornified squamous epithelial cells, clustered with irregular shape. These start to be replaced by leukocytes during metestrus, being these predominant cells in diestrus (Pfaus, 1999;Caligioni, 2009;Snoeren, 2018). Experiments were performed when females were in proestrus/estrus (sexually receptive) or in diestrus (sexually nonreceptive).

Ex vivo electrophysiological recordings
PR-Cre-EYFP mice were deeply anesthetized with isoflurane and decapitated 1-2 h after the phase of the cycle identification by pap smearing. This protocol is used by the laboratory for behavioral experiments (Nomoto and Lima, 2015). After decapitation, brains were quickly removed and placed into "ice-cold" slicing solution containing the following: 0.66 mM kynurenic acid, 3.63 mM pyruvate, 2.5 mM KCl, 1.25 mM NaH 2 PO 4 , 26 mM NaHCO 3 , 10 mM D-glucose, 230 mM sucrose, 0.5 mM CaCl 2 , 10 mM MgSO 4 , and bubbled with 5% CO 2 and 95% O 2 . Coronal sections with 300 mm containing the VMH were cut using a vibratome (Leica VT1200) in the same ice-cold slicing solution. Slices were recovered in oxygenated artificial CSF (ACSF) containing the following: 127 mM NaCl, 2.5 mM KCl, 25 mM NaHCO 3 , 1.25 mM NaH 2 PO 4 , 25 mM D-glucose, 2 mM CaCl 2 , and 1 mM MgCl 2 , at 34°C for 30 min and stored in the same solution at room temperature before being transferred to the recording chamber.
Whole-cell recordings were made under a SliceScope Pro (Scientifica) microscope with the slices submerged in ACSF. A beam of blue light generated with a LED illumination system (CoolLED, pE300) was passed through an immersion objective (Olympus, 40Â) to achieve wide field illumination of the slice. PR1 neurons were identified by their natural fluorescence under the microscope and PRby the absence of fluorescence. A standard GFP filter cube (Thorlabs) was used to filter EYFP emission. EYFP expression, and therefore fluorescence levels, are independent of the reproductive cycle, as after Cre recombination, the expression of the fluorescent protein is dependent on the rosa locus. Patch recording pipettes (resistance 3-5 MV) were filled with internal solution containing the following: 135 mM K-gluconate, 10 mM HEPES, 10 mM Na-phosphocreatine, 3 mM Na-L-ascorbate, 4 mM MgCl 2 , 4 mM Na-ATP, 0.4 mM Na-GTP (pH 7.2 adjusted with NaOH and osmolarity ;292 mOsm), and 0.1% of biocytin. Internal solution was filtered with a 0.2-mm pore size cellulose acetate filter tip. A Multiclamp 700B amplifier and digitized at 10 kHz with a Digidata 1440a digitizer (both from Molecular Devices) were used and the data were filtered online with a 10-kHz low-pass filter.
We applied a test-pulse of -10 mV for 100 ms in voltageclamp mode at -70 mV to monitor membrane capacitance (Cmem) and series resistance. To determine excitability and firing properties, current-clamp experiments were performed at resting potential by applying series of 1000-ms current pulses between 0 and 240 pA with an increasing interval of 30 pA at a frequency of 0.1 Hz. Hyperpolarization-activated inward current (Ih current) and rebound firing were measured with a 1000-ms current pulse at -90 pA. The bridge balance was semi-automatically compensated using Axon pCLAMP 10 (Molecular Devices).
After recording, the pipette was carefully withdrawn from the neurons to allow the resealing of the membrane and the diffusion of biocytin inside the neuron. Slices were then transferred to a multi-well plate containing 4% paraformaldehyde (PFA) for ;1 h for fixation and then transferred to 0.01 M PBS for storage before immunostaining.

Histology
To visualize biocytin filling, slices were incubated in streptavidin-Alexa Fluor 488 conjugate (Invitrogen; concentration, 1:200) with 0.3% Triton X-100 in 0.01 M PBS at room temperature for 3 h in dark. After washing three times with 0.01 M PBS for 30 min each, sections were mounted in glass slides (Menzel-Glazer), coverslipped (Marienfeld) with Mowiol mounting medium (Sigma-Aldrich) and the edges sealed with clear nail polish. For tissue collection, VGat-Cre-tdtomato mice were deeply anesthetized with a lethal amount of a mixture of 12% of ketamine (Imalgene 1000, Merial) and 8% of xylazine (Rompun 2%, Bayer) in saline solution and perfused transcardially with 0.01 M PBS followed by 4% PFA in PBS. The brains were removed, fixated in the 4% PFA solution for ;1 h and transferred to a 30% sucrose (Sigma-Aldrich) in 0.01 M phosphate buffer and 0.1% sodium azide (ACROS Organics) to allow cryopreservation. Brains were then frozen and coronal sections with 45-mm thickness were obtained from the VMH using a freezing sliding microtome (SM2000R, Leica).
Histologic sections and fixed ex vivo brain slices were imaged with a Zeiss LSM 710 confocal laser scanning microscope with a 10Â and a 25Â magnification objectives.

Electrophysiological data analysis
A total of 144 recorded neurons were used in the present study. Please see Table 1 and Figure 1 for the numbers of recorded PR1 and PR-neurons, across the reproductive cycle and AP axis.
The action potentials recorded in current-clamp mode were analyzed with custom-written MATLAB code. Each spike within the spike trains obtained on current injection were detected and used to determine the firing frequency, interspike interval (ISI), threshold to spike, latency to spike and coefficient of variation calculated at a fixed current injection of 150 pA). The coefficient of variation was calculated using the following formula (Holt et al., 1996): The rebound action potentials observed after injecting hyperpolarizing currents were obtained using the same method.
The first action potential obtained of each spiking sweep was used to calculate the spike amplitude and spike half-width to prevent frequency-dependent artifacts in the spike shape caused by sodium and potassium channel inactivation.
The resting potential was determined by averaging the membrane voltage values before the current injection and the Ih current amplitude was calculated by the difference between the early and late phase of a hyperpolarizing current pulse.
Voltage clamp test-pulses were analyzed using Clampfit 10.7 Software (Molecular Devices, LLC) and access resistance (Ra), membrane resistance (Rm), time constant (t ), and capacitance were calculated. Only neurons with Ra 1/ 10th lower than the Rm were used for analysis. Ra was comparable across location, genotype and the phase of the cycle in the present study.
The passive membrane properties Rm, Cmem, and membrane t were obtained immediately after membrane rupture, using a square voltage step (À10 mV, 100 ms). The Ra was determined by measuring the amplitude of the current response to the command voltage step and the Rm as the difference between the baseline and the holding current in the steady state after the capacitive decay, by applying the Ohm's law. The membrane t was determined by a single exponential fit of the decay phase in response to the square pulse. An approximation of the capacitance was made by using the formula: Figure 1. The predominant firing pattern of PR1 neurons varies across the AP axis and differs from the PR-population. A, Low-dimensional embedding of VMHvl neurons input current to output voltage responses (I-V curves) using tSNE. Four representative examples of the different firing responses found correspond to tonic neurons with high firing rate (A1), tonic neurons with low firing rate (A2), phasic neurons (A3), and tonic neurons with depolarization block (A4). B, Sum of squared distances between each data point and its assigned cluster as a function of increasing number of possible clusters. C, D, Tonic neurons with depolarization block, which initially were largely clustered with the phasic neurons, can be distinguished by obtaining the ratio between the maximum firing rate and the firing rate at the maximum current injected (240 pA). E, The majority of PR1 neurons in the anterior and medial VMHvl display tonic firing, while the posterior PR1 neurons have mainly phasic firing (p , 0.001), with no changes across the reproductive cycle. F, The PR-neurons display mainly tonic firing across the whole AP axis and do not change across the reproductive cycle. x 2 test; **p , 0.01 t ¼ membranecapacitance p accessresistance: The principal component analysis (PCA)-initialized tdistributed stochastic neighbor embedding (tSNE) on the firing rates and its subsequent clustering were performed using the tsne and kmeans functions of MATLAB 2020a.

Neuronal reconstruction and image data analysis
A total of 180 neurons VMHvl neurons were filled with biocytin in acute slices obtained from PR-EYFP female mice using whole-cell recording pipettes. These slices containing biocytin filled neurons were fixed and stained with streptavidin-Alexa Fluor 488 conjugated for the subsequent reconstruction of neuronal structure. In our hands, the native EYFP signal systematically bleached after slice fixation, being several fold dimmer than the fluorescence obtained with streptavidin-Alexa Fluor 488. Therefore, with our confocal settings, no somatic EYFP signal was obtained despite their spectral overlap. Please refer to Table  2 for the numbers of filled PR1 and PR-neurons across the reproductive cycle and AP axis.
To determine the bregma coordinate for each reconstructed neuron and ensure that its location fell within the VMHvl, in a separate experiment we first assessed the extension of the VMHvl, taking advantage of the VGAT-tdTomato. The VMH is organized in a core of glutamatergic neurons (Ziegler et al., 2002;Yamamoto et al., 2018) and a surrounding shell with higher density of GABAergic neurons (Yamamoto et al., 2018). Therefore, we assessed the AP extension of the VMH by counting consecutive slices in which the VMH was detectable, as shown by a decrease of the VGAT-tdTomato signal (Fig. 6A) and multiplying by the thickness of the histologic slices. A sharp transition from low to high tdTomato fluorescence was observed both in its anterior limit (near the anterior hypothalamus) as well as with its posterior limit (near the premammilary hypothalamus; data not shown). The extension obtained matched the extension illustrated by the Paxinos Brain Atlas: anterior limit = 1.05 to posterior limit = 2.06 from bregma (that is, the VMHvl spans for ;1 mm). For each mouse, we systematically sampled three consecutive levels of the VMH in slices of 300 mm of thickness. We thus categorized each sampling level as anterior (bregma À1.05 to À1.30 mm approximately), medial (bregma À1.30 to À1.65 mm approximately), and posterior (bregma À1.65 to À2.00 mm approximately) levels of the VMH (Fig. 6B,C). Localization of the reconstructed neurons in the brain was assessed by matching in Adobe Illustrator CS6 (Adobe Systems Incorporated) the confocal images of the neurons with adapted sections from the Paxinos brain atlas (Franklin and Paxinos, 2008). Only the neurons inside the defined boundaries of the VMHvl were considered for quantification. We used the simple neurite tracer (Longair et al., 2011) package from Fiji/ImageJ software (Schindelin et al., 2012) to analyze the morphologic properties and Sholl profiles of the neurons filled with biocytin. The Cell Counter Fiji/ImageJ plugin was used to manually count spines.
Cut and resealed dendrites (identified as a globular thickening in the extreme of a dendrite) were quantified for every neuron and did not reveal significant differences across groups (Number of cut dendrites/total number of dendrites quantified per group for PR1 neurons: 21/126 in anteriorNon-Rec, 19/120 in anteriorRec, 20/ 116 in medialNon-Rec, 31/151 in medialRec, 26/131 in posteriorNon-Rec and 26/121 in posteriorRec. And for PR-neurons: 26/130 in anteriorNon-Rec, 16/70 in anteriorRec, 19/95 in medialNon-Rec, 15/85 in medialRec, 26/100 in posteriorNon-Rec and 17/100 in posteriorRec). When tested with a Fisher's exact test, neither PR1 nor PRrevealed statistical difference in the proportion of cut dendrites. We thus assume that this unavoidable artifact of our method of neuronal recording and reconstruction did not add any group-specific bias in the structural quantification.
In our hands, as well as in previous studies (Calizo and Flanagan-Cato, 2000), the intracellular labeling of VMHvl neurons did not always make their thin axonal process visible, and therefore, the axons identified were not considered for analysis.
Graphs were made using custom made MATLAB code or GraphPad Prism 8 Software.

Statistical analyses
Statistical analyses were performed using GraphPad Prism 8 Software. Normality of the residuals was tested with the D'Agostino-Pearson omnibus K2 test. When normally distributed, three-way ANOVA tests were performed to compare groups in different phases of the cycle (D vs PE) versus location in the AP axis (anterior vs medial vs posterior) versus genotype (PR1 and PR-neurons), using the Sidak test to correct for multiple comparisons. In Figure 5E, a two-way ANOVA, with Sidak test for multiple comparisons, was used to investigate the origin of the differences in the location factor that were specific for PR-neurons. In the properties whose residuals did not pass the normality test, a logarithmic or square root transformation was applied and normality was consequently reassessed. As specified in the legend, in the properties that the transformations did not make the residuals become normally distributed, we performed a Kruskal-Wallis test followed by a Dunn's post hoc test for multiple comparisons.
For Figure 2A-F, a mixed-effects test with repeated measures, and Figure 7C-H a two-way ANOVA with repeated measures were performed. For Figures 1E-F, 5B, C, a x 2 test was performed, the pairwise comparisons we tested followed the same structure as the one obtained for the three-way ANOVA multiple comparisons. Whisker plots represent median with interquartile range. Error bars represent mean 6 SEM p values were rounded to two decimal places. Significance was noted as *p , 0.05, **p , 0.01.

Results
The proportion of tonic and phasic firing PR1 neurons varies across the anteroposterior axis In order to investigate the neuronal excitability of PR1 and PR-neurons, we recorded current to voltage inputoutput response curves (I-V curves) at resting potential in current-clamp mode in acute slices of naturally cycling PR-Cre x EYFP adult female mice. Briefly, the female reproductive stage was determined by vaginal lavage and slices were obtained from females in the least receptive state (diestrus) and the most receptive state (proestrus/ estrus). PR1 neurons were identified by their natural fluorescence under the microscope and PR-by the absence of fluorescence (for details, see Materials and Methods).
We performed a PCA-initialized tSNE analysis on the I-V curve firing rates obtained for our full sample of recorded PR1 and PR-neurons. To understand how many clusters can be optimally distinguished in the resulting embedded data (Fig. 1A), we performed k-means clustering with a possible number of clusters varying between one and eight. When quantifying the sum of squared distances between each datapoint and its assigned cluster centroid over an increasing number of clusters, we observed a characteristic "elbow" shape when there were three possible clusters, suggesting this is the optimal number of clusters in our dataset (Fig. 1B). The first two clusters corresponded to tonic neurons with high and low firing rates (Fig. 1A1,A2, respectively). Tonic firing neurons exhibited a linear increase in firing rate in response to the increases in the magnitude of current injected (Fig.  1C, top). The third cluster roughly corresponded to phasic firing neurons (Fig. 1A3), neurons in which the action potentials did not linearly increase in response to increased current injection (Fig. 1C, bottom). These neurons exhibited a step function-like behavior, where after reaching the firing threshold, the number of action potentials remained constant with increased amounts of current. However, when visually inspecting the firing profiles of neurons in this cluster (Fig. 1C), we observed that this cluster included phasic neurons and neurons in which the firing rate initially increased linearly but then entered depolarization block and drastically reduced their firing rate (Fig.  1A4). When dividing the peak firing rate of each neuron by their firing rate at the maximum current injected (240 pA), a subset of neurons with particularly high ratio values was clearly distinguishable from the rest (Fig. 1D). These neurons were tonic firing neurons with depolarization block.
We therefore considered two major categories of neurons. The first category being tonic firing neurons, comprising neurons that did not reach the depolarization block within the range of currents we injected (Fig. 1A1, A2, tonic neurons), and neurons that did (Fig. 1A4, tonic neurons with depolarization block), independently of their firing rate. And the second category being phasic firing neurons (Fig. 1A3).
The majority of the anterior and medial PR1 neurons displayed tonic firing, while most of the posterior ones displayed phasic firing (x 2 test, p , 0.001; Fig. 1E), with no differences across the reproductive cycle. In contrast, the PR-population is homogenously composed of a higher proportion of tonic firing neurons across the AP axis (Fig. 1F).
To summarize, the VMHvl is composed of multiple electrophysiological phenotypes and posterior PR1 neurons have a higher density of phasic neurons in comparison to the anterior and medial levels of the VMHvl.
Local changes in the intrinsic excitability and threshold to spike of PR1 neurons across the reproductive cycle The analysis of the I-V curves revealed that in the anterior and posterior VMHvl, the excitability of PR1 neurons does not to change across the cycle ( Fig. 2A,C). In contrast, medial PR1 neurons of non-receptive females showed significantly higher excitability when compared with the excitability of medial PR1 neurons originating from receptive females, particularly in response to higher input current (Fig. 2B, top, strong interaction between the phase of the cycle and the amount of current injected; mixed-effects test, F (8,160) = 6.87, p , 0.0001). This modulation is caused by a change in the excitability of the tonic firing medial neurons across the reproductive cycle while the behavior of medial phasic neurons remained unchanged (Fig. 2B, bottom). In addition, this change was specific to the PR1 population since no changes were observed across the cycle in the PR-neurons at any of the AP levels ( Fig. 2D-F).
We also observed that the PR1 population was less excitable than the PR-in the anterior and posterior levels of the VMHvl (mixed-effects test, anterior PR1 vs anterior Insets, Representative recordings of voltage response from PR1 neurons in response to the same current injection. Anterior n = 11 NR, n = 13 R, medial n = 11 NR, n = 12 R, posterior n = 10 NR, n = 13 R. D-F, IV curves across the reproductive cycle reporting no changes in the firing frequency of PR-neurons after injecting 1000-ms current pulses between 0 and 240 pA (30-pA steps). Insets, Representative recordings of voltage response from PR-neurons in response to the same current injection. Anterior n = 15 NR, n = 11 R, medial n = 10 NR, n = 13 R, posterior n = 12 NR, n = 13 R. A, C, D, F, PR1 neurons are less excitable than PR-in the anterior and posterior VMHvl (F (8,360) = 3.53, p , 0.001 and F (8,361) = 2.83, p , 0.01, respectively) but not in (B, E) the medial. Mixed-effects test with repeated measures. Mean (6SEM); *p , 0.05 PR-F (8,360) = 3.53, p , 0.001 and posterior PR1 vs posterior PR-F (8,361) = 2.83, p , 0.01) but not in the medial subdivision, probably in part because of the increased excitability of medial PR1 neurons from non-receptive females.
It is worth mentioning that the Rm did not vary across experimental groups (Table 1) and only a small decrease was observed in the capacitance of the anterior and posterior PR1 of non-receptive females (three-way ANOVA, location/phase interaction F (2,132) = 3.20, p = 0.04). The membrane t also did not vary across the reproductive cycle, however posterior PR1 neurons exhibited a smaller t when compared with the anterior PR1 population (three-way ANOVA, location F (2,132) = 3.23, p = 0.04, multiple comparison anterior vs posterior p = 0.04).
The resting potential of anterior PR1 neurons of receptive females was lower when compared with non-receptive females (Fig. 3A,B; three-way ANOVA, phase F (1,132) = 6.39, p = 0.01, multiple comparison antNon-Rec vs antRec p , 0.01). However, this difference was not large enough to produce changes in the threshold to spike (Fig. 3C). The latency to spike of PR1 and PRneurons did not vary across the cycle (Fig. 3F,H). In addition, we observed a mild, yet significant, lower latency to spike of the posterior neurons when compared with anterior (Fig. 3F,H; three-way ANOVA, location F (2,132) = 5.48, p , 0.01, multiple comparison anterior vs posterior p = 0.01) and medial neurons (multiple comparison medial vs posterior p = 0.01) that was independent of the phase of the cycle and the genotype.
The regularity of the recorded spike trains was similar across experimental groups, as shown by comparable CV2 values across phase and location for both the PR1 and PR-populations (regularity was determined at 150 pA of injected current; Fig. 3G,I).
In focusing on the action potential shape, we observed a comparable amplitude and half-width of the action potentials across the reproductive cycle in PR1 and PRneurons ( Fig. 4A-E). The spike amplitude of the posterior neurons of both PR1 and PR-populations is moderately but consistently smaller compared with anterior neurons (Fig. 4B,D; three-way ANOVA, location F (2,132) = 5.74, p , 0.01, multiple comparison anterior vs posterior p = 0.03) and medial (multiple comparison medial vs posterior p , 0.01).
Together, these results suggest that the electrophysiological properties of VMHvl neurons vary across the reproductive cycle, but the observed changes depend on their location in the AP axis. In the sexually receptive phase, anterior PR1 neurons are more hyperpolarized in their resting state and the medial PR1 population exhibits lower excitability.
The anterior VMHvl presents a higher proportion of rebound firing neurons during the non-sexually receptive phase of the cycle Finally, and given that it is well established that hyperpolarizing input can trigger rebound depolarizing responses increasing the firing rate of neurons in a wide variety of brain areas including the hypothalamus (Burdakov et al., 2004;Israel et al., 2008), we sought to characterize the response of VMHvl neurons to hyperpolarizing input. To do so, we quantified the proportion of VMHvl neurons which exhibited rebound firing after hyperpolarization ( Fig. 5A; -90 pA for 1 s). We observed that 6 out of 32 PR1 neurons of non-receptive females and only two out of 38 PR1 neurons of receptive females fired rebound action potentials (Fig. 5B). This difference originates mainly from the anterior PR1 neurons (x 2 test, antNon-Rec1 vs antRec1 p = 0.01), where four out of 11 neurons of non-receptive females displayed rebound firing while in receptive females none of the 13 recorded neurons produced rebound action potentials. Furthermore, no significant differences were observed between PR1 neurons across the AP axis nor between the PR1 and the PR-populations (Fig. 5B,C).
Rebound depolarization is caused by the hyperpolarization activated inward current (Ih) which represents a powerful modulator of neuronal firing frequency and timing (Momin et al., 2008). Thus, we quantified the Ih current amplitude in VMHvl neurons (Fig. 5D,E). A significant difference was obtained in the interaction between location and genotype (three-way ANOVA, location/genotype interaction F (2,132) = 4.36 p = 0.01), that however, was not large enough to yield statistical differences across groups in the multiple comparisons test. To pinpoint the origin of such changes, we analyzed separately PR1 and PRpopulations. Our results show that while the Ih magnitude was unchanged across the AP axis for the PR1 population, the anterior PR-neurons have modestly higher Ih magnitude compared with medial and posterior PR-neurons (two-way ANOVA, location F (2,68) = 6.30, p , 0.01, multiple comparison anterior vs medial p = 0.04, anterior vs posterior p , 0.01), however, these differences are not reflected in a significantly different fraction of rebound anterior PR-neurons compared medial and posterior PRneurons (Fig. 5C).
To summarize, the proportion of neurons with rebound firing is fairly homogenous in the VMHvl across the reproductive cycle and independently of genotype, with the exception of its most anterior subdivision, where PR1 neurons exhibit cyclical alterations. Given that these changes are not accompanied by a modulation of Ih current, they likely reflect the change in resting membrane potential that drives receptive neurons away from their firing threshold, causing Ih-mediated depolarization to fail to produce rebound firing (Figs. 3B, 5A).

Localized structural plasticity of VMHvl PR1 neurons across the reproductive cycle
To study the impact of the phase of the reproductive cycle on the morphologic properties of VMHvl neurons, we reconstructed and performed morphometric quantifications of the soma area, number of dendrites, primary dendrites and branch points per neuron of PR1 and PRneurons in different locations across the AP axis of the VMHvl. As previously reported, the VMHvl is primarily composed of glutamatergic neurons, that are surrounded by GABAergic neurons (Yamamoto et al., 2018). In fact, the absence of GABAergic somas has been used to delineate the boundaries of the VMH (Jarvie and Hentges, 2012). Taking advantage of mice expressing the fluorescent reporter tdTomato under the promoter of the vesicular transporter of GABA (VGat-tdTomato; Kaneko et al., 2018) we observed, as expected, lower signal intensity inside the boundaries of the VMH (Fig. 6A) that was consistent across the AP regions of the VMH. The limits obtained with the inspection of GABAergic markers were used to delineate the boundaries of the VMHvl, and ensure that the neurons characterized in this study were indeed within the nucleus (Fig. 6B,C).
The somatic area was unaltered across the reproductive cycle and was not different between the PR1 and PR-populations (Table 2). We observed that PR1 neurons of receptive females exhibit a robust increase in the number of dendrites per neuron across the AP axis (Table 2; three-way ANOVA, genotype/phase interaction F (1,173) = 7.85, p = 0.01), and a higher number of branching points per neuron compared with neurons originating from non-receptive females (Table 2; three-way ANOVA, genotype/phase interaction F (1,173) = 5.81, p = 0.02). The changes across the reproductive cycle are specific for the PR1 population, as they were not observed in PR-neurons (Table 2). Both in medial PR1 and PR-neurons, independent of the phase of the cycle, we observed a moderately lower number of primary dendrites per neuron, compared with the anterior and posterior neurons which yielded small yet significant differences across the AP axis (three-way ANOVA, location F (2,173) = 3.14, p = 0.04).
Neurons in the VMHvl are characterized by having a LPD that can be several fold longer than its SPDs (Calizo and Flanagan-Cato, 2002). Therefore, we analyzed LPD and SPD separately (Table 2). While the length of the LPD of PR1 was unchanged across the reproductive cycle, PR-neurons exhibited longer LPDs in non-receptive females (three-way ANOVA, genotype/phase interaction F (1,168) = 6.67, p = 0.01). No overall significant differences were found in the SPD length of PR1 and PR-neurons across the cycle and AP axis. Interestingly, the posterior neurons were shown to have shorter LPDs than the anterior and medial neurons regardless of the phase of the cycle or the genotype (three-way ANOVA, location F (2,168) = 9.73, p , 0.0001). No differences in the length of LPDs and SPDs were found between the PR1 and PRpopulations.
In addition to the analysis of dendritic length, we sought to investigate whether the complexity of the dendritic  (HW). B, D, The amplitude of the first action potential for each current injection step of the PR1 and PRneurons does not change across the reproductive cycle; however, it differs mildly across the AP axis (location F (2,132) = 5.74, p , 0.01, multiple comparison L*, anterior vs posterior p = 0.03 and L**, medial vs posterior p , 0.01). C, E, The spike half-width of PR1 and PR-neurons is consistent across the reproductive cycle and the AP axis. PR1 anterior n = 11 NR, n = 13 R, medial n = 11 NR, n = 12 R, posterior n = 10 NR, n = 13 R. PR-anterior n = 15 NR, n = 11 R, medial n = 10 NR, n = 13 R, posterior n = 12 NR, n = 13 R. Three-way ANOVA, Sidak correction for multiple comparison. Medians (625th-75th percentile); *p , 0.05, **p , 0.01. Individual data points are shown per group color coded according their tonic or phasic firing type. trees was modified across the reproductive cycle. To do so, we used the Sholl method that measures the number of dendritic processes intersections as a function of the radial distance from the soma. No changes were observed in the maximum number of dendritic intersections of PR1 and PR-neurons across the phase of the cycle (Fig. 7A-H) nor across their location in the VMHvl, indicating that these neurons reach comparable maximum complexities in their dendritic trees. Nevertheless, the analysis of the Sholl profiles of these neurons revealed that the posterior PR1 neurons of non-receptive females had a significantly lower complexity of their dendritic tree (Fig.  7E, repeated measures ANOVA p , 0.0001) compared with those from females in the receptive phase. In addition, the medial PR-neurons of non-receptive females had a higher complexity than their counterparts in receptive females (Fig. 7G, repeated measures ANOVA F (40,136) = 3.10, p , 0.0001).
The increased complexity observed in the posterior PR1 neurons of females in the receptive phase was large enough to provide differences when tested with the same method specifically in the proximal branching ( Fig. 7E; ,250 mm from the soma), while the increased complexity of the medial PR-neurons of non-receptive females were observed at more distal parts of the soma (240-340 mm) and were large enough to yield multiple comparisons significant changes (Fig. 7G).
We also observed that medial PR1 neurons have lower dendritic branching than medial PR-neurons ( Fig. 7G; repeated measures ANOVA F (35,2030) = 2.84, p , 0.0001), particularly in the distal parts from the soma (320-380 mm). The posterior PR1 neurons have lower proximal Figure 5. The proportion of anterior PR1 neurons with rebound firing changes across the reproductive cycle. A, Example of recordings of voltage response at -90 pA. Proportion of neurons presenting rebound firing after hyperpolarization of (B) PR1 neurons from non-receptive and receptive females, and from (C) PR-neurons from non-receptive and receptive females. The anterior PR1 neurons of non-receptive females display a significantly higher fraction of neurons with rebound spikes compared with receptive females (p = 0.01) x 2 test. D, The Ih amplitude of the PR1 neurons does not change across the reproductive cycle nor the AP axis, while the (E) Ih amplitude of PR-neurons varies across the AP axis (three-way ANOVA location/genotype interaction F (2,132) = 4.36 p = 0.01, two-way ANOVA for PR-, location F (2,68) = 6.30, p , 0.01, multiple comparison L*, anterior vs medial p = 0.04, L**, anterior vs posterior p , 0.01). Three-way ANOVA and (E) two-way, Sidak correction for multiple comparison; *p , 0.05. Individual data points are shown per group color coded according their tonic or phasic firing type. dendritic branching than the posterior PR-neurons ( Fig.  7H; repeated measures ANOVA F (40,2520) = 1.77, p , 0.01). These observations indicate that PR1 neurons have a different dendritic complexity proximal to the cell body compared with that of the PR-population, independent of the phase of the reproductive cycle.
To summarize, similarly to what we report for the electrophysiological properties, the structural properties of PR1 and PR-neurons are diverse, with some varying across the AP axis and the reproductive cycle.

Dendritic spine density largely differs between the PR1 and PR-populations
Previous studies have shown that externally primed estrogen exerts effects on the dendritic spine density of VMHvl neurons that are different depending on the type of dendrite, specifically, an increase in spine density on SPDs of VMHvl neurons (Calizo and Flanagan-Cato, 2000) and a decrease of the spine density on LPDs of VMHvl neurons expressing the receptor for estrogen (Calizo and Flanagan-Cato, 2002). In order to determine whether such modulation of the spine density of VMHvl neurons is present in the physiological range of sex hormone fluctuation, the dendritic spine density of PR1 and PR-neurons was analyzed depending not only on their location in the VMHvl, but also on the primary dendrite category across the reproductive cycle of naturally cycling females. Although overall PR1 and PR-neurons have comparable spine densities on SPDs across the reproductive cycle (Fig. 8A,C) and location, posterior VMHvl neurons in non-receptive females have higher spine density on SPDs compared with posterior VMHvl neurons in the receptive phase resulting in a significant interaction between the phase of the cycle and location within the VMHvl (three-way ANOVA, phase/location interaction F (2,153) = 3.88, p = 0.02). Furthermore, the PR1 population has significantly lower spine density on the SPDs than PR-neurons ( Fig. 8E; three-way ANOVA, genotype F (1,153) = 26.21, p , 0.01).
Similar to what was observed for the spine densities of SPDs, for LPDs we observed that the posterior VMHvl neurons of receptive females present a decreased spine density compared with posterior neurons from females that were non-receptive, which in the case of LPDs the change is specific to PR1 neurons (Fig. 8B,D; three-way In addition, by comparing PR1 with PR-neurons, we observed that PR1 neurons have lower spine density on the LPDs (Fig. 8F; three-way ANOVA, genotype F (1,150) = 7.25, p = 0.01).
Altogether, we observed that neurons expressing PR have an overall lower spine density when compared with their neighboring PR-neurons. Across the reproductive cycle, posterior PR1 neurons undergo more pronounced structural changes, exhibiting even lower spine density in the receptive phase.

Discussion
We combined whole-cell recordings with labeling of individual neurons to investigate the structural and electrophysiological properties of PR1 and PR-neurons along the AP axis of the VMHvl and across the reproductive cycle of naturally cycling female mice. First, we report the existence of multiple electrophysiological phenotypes within the VMHvl, including tonic firing neurons (with and without depolarization block) and phasic firing neurons. We show that PR1 cells are distinct from PR-in three major aspects: first, because of structural properties (PR1 neurons have lower spine densities in general); second, because of the great extent to which the properties of PR1 neurons vary across the AP axis; and third because of local changes in the properties of PR1 across the reproductive cycle (which are minimal for the PRpopulation). Our results further support the existence of subdivisions in the VMHvl and its possible role in coordinating female behavior with the internal reproductive state.
The impact of the reproductive cycle on VMHvl function has been interrogated in vivo and in vitro, but most studies were performed in females whose ovaries were removed and then supplemented with estrogen and progesterone (Millhouse, 1973;Barfield, 1980, 1983a,b;Figure 7. The dendritic properties of PR1 neurons are differently modulated across the reproductive cycle depending on the location in the AP axis and differ from the PR-population. A, Schematic representation of the Sholl analysis of a posterior PR1 neuron of a non-receptive female and a posterior PR1 neuron of a receptive female (the axon, visible in this case, depicted in lighter shade and not considered for analysis). B, Confocal image of the posterior PR1 neuron of a receptive female depicted in A filled with biocytin and stained with streptavidin-Alexa Fluor 488 conjugated. Scale bar: 50 mm. Sholl profiles of (C-E) PR1 anterior n = 19 NR, n = 15 R, medial n = 17 NR, n = 19 R, posterior n = 20 NR, n = 16 R and (F-H) PR-neurons anterior n = 16 NR, n = 10 R, medial n = 12 NR, n = 12 R, posterior n = 15 NR, n = 14 R across the reproductive cycle. E, The posterior PR1 neurons of receptive females show higher proximal complexity than the ones from non-receptive females (distance from soma Â phase of the cycle F (40,136) = 3.10, p , 0.0001, multiple comparison p , 0.001 at 120 mm, p , 0.0001 at 140-160 mm, p , 0.01 at 180-200 mm, and p = 0.04 at 220 mm). G, Medial PR-neurons from receptive females show lower complexity than the ones from non-receptive females (distance from soma Â phase of the cycle F (33,726) = 2.27, p , 0.0001, multiple comparison p = 0.03 in 320-340 mm). The medial PR1 neurons have lower distal branching than the medial PR-neurons (distance from soma Â genotype F (35,2030) = 2.84, p , 0.0001, G*, multiple comparison p = 0.04 at 320 mm, p , 0.01 at 340-360 mm p = 0.03 at 380 mm) and (H) the posterior PR1 neurons have lower proximal branching than the PR-neurons (distance from soma Â genotype F (40,2520) = 1.77, p , 0.01). Two-way ANOVA with repeated measures. Mean (6SEM). *p,0.05. Flanagan-Cato, 2000, 2002;Griffin and Flanagan-Cato, 2008;Griffin et al., 2010; for examples of naturally cycling studies, see Nomoto and Lima, 2015;Inoue et al., 2019). While these manipulations are convenient, easing experimental planning and decreasing interindividual variability (with the extra benefit of allowing a reduction in the number of females used), ovariectomized females with hormonal replacement differ substantially from intact females: first, before hormonal replacement the levels of sex hormones are extremely low, as the main source of estrogen and progesterone is absent and second because the hormonal treatment exposes females to concentrations of sex hormones that differ from physiological levels (Liu and Shi, 2015). The expression of hormone receptors, such as the PR, are under the control of sex hormone levels as well (MacLusky and McEwen, 1978), meaning that that the neuronal response is affected at the level of the receiver (receptor) and message (sex hormone). Persistent hormonal replacement of ovariectomized females leads to tumor development (Kordon et al., 1993) further suggesting that the hormonal treatment leads to undesired physiological effects. Finally, the sexual behavior that hormonally treated females exhibits differs from the behavior of naturally cycling females (Zipse et al., 2000), suggesting that the treatment fails to fully recapitulate the effects of natural sex hormones levels. To the best of our knowledge, this is the first study interrogating the intrinsic properties of neurons across the reproductive cycle of naturally cycling females and therefore the results observed reflect endogenous changes. However, we cannot claim that the changes we observe across the reproductive cycle are an effect of sex hormone levels because those were not directly manipulated. Manipulations with more naturalistic levels of sex hormones or local manipulations in the expression levels of sex hormone receptors should be employed to establish a direct causal link between hormonal levels, the structural and physiological properties of VMHvl neurons and the changes observed.
In the present study, we report, for the first time, the existence of different neuronal types in the VMHvl according to their firing pattern. We found tonic and phasic neurons, in both PR1 and PR-populations, intermingled across the AP axis of the VMHvl. However, the presence of phasic neurons was particularly abundant in the posterior VMHvl, where it comprised more than half of the recorded PR1 neurons.
This AP gradient may have strong implications not only for the role of different portions of the VMHvl in controlling behavioral output, but also for the interpretation of activity manipulations, such as the optogenetic control of neural activity. While anterior VMHvl neurons can reach high Figure 8. The spine density of the LPD of posterior PR1 neurons changes across the reproductive cycle. A, Spine density of the SPDs anterior n = 15 NR, n = 13 R, medial n = 15 NR, n = 17 R, posterior n = 15 NR, n = 14 R and the (B) LPDs anterior n = 15 NR, n = 13 R, medial n = 15 NR, n = 17 R, posterior n = 13 NR, n = 14 R of PR1 neurons across the reproductive cycle. The spine density of the LPDs of posterior PR1 neurons is reduced in receptive females (multiple comparison p = 0.0170, causing a change in the three-way interaction between the phase of the reproductive cycle, the location in the AP axis, and the genotype on the LPDs, F (2,150) = 4.81, p = 0.01). C, Spine density of the SPDs anterior n = 15 NR, n = 10 R, medial n = 12 NR, n = 17 R, posterior n = 15 NR, n = 12 R and the (D) LPDs anterior n = 15 NR, n = 10 R, medial n = 12 NR, n = 11 R, posterior n = 15 NR, n = 12 R of PR-neurons across the reproductive cycle. E, PR1 neurons have lower spine density on the SPDs (G**, genotype F (1,153) = 26.21, p , 0.01) and also (F) lower spine density on the LPDs (G*, genotype F (1,150) = 7.25, p = 0.01) than the PR-neurons. Three-way ANOVA, Sidak correction for multiple comparison. Medians (625th-75th percentile); *p , 0.05, **p , 0.01 firing frequencies, posterior phasic neurons present a robust stop of action potential generation after reaching the spiking threshold that does not depend on the current injected in the cell. Importantly, this may be particularly relevant for the interpretation of results where different laser intensities produce diverse behavioral outputs (Lee et al., 2014;Kunwar et al., 2015), as optogenetic activation at particular AP levels may bias the activation of neuronal populations with different intrinsic properties.
Other hypothalamic neuronal populations, such as the tuberoinfundibular dopamine neurons (Lyons et al., 2012) or the gonadotropin-releasing hormone neurons (Jarry et al., 1991) present in the tuberal and preoptic hypothalamus, are known to display phasic firing behavior. Furthermore, the presence of intermingled tonic and phasic firing neuronal population has been previously reported in the magnocellular nucleus of the hypothalamus (Israel et al., 2016), where oxytocin neurons show tonic and vasopressin neurons phasic firing behavior. The differential role of tonic and phasic firing neurons in the control of goal directed behaviors (Grace et al., 2007;Budygin et al., 2020) but also in their contributions for encoding stimulus-environment features (Pilkiw et al., 2017) has been widely studied throughout the brain. In addition, it is known that changes in voltage-sensitive channels (Suzuki and Rogawskit, 1989), that can be triggered by neuromodulation (Lyons et al., 2012;Valentino and Van Bockstaele, 2015), are able to mediate the dynamic transition between tonic to phasic firing regimes. In the VMHvl, whether tonic and phasic firing neuronal populations, or the dynamic transition between these firing regimes, contribute to different aspects of socio-sexual behavior remains elusive and should be the focus of future studies.
Our results indicate that PR1 cells are distinct from their neighbors, a difference that is probably established early in development, and later on by the fact that their intrinsic properties of this population can be directly affected by estrogen and progesterone, as PR1 neurons co-express ER (Blaustein and Turcotte, 1985;Hashikawa et al., 2017;Sá and Fonseca, 2017). Interestingly, the properties of PR-neurons were similar along the AP axis and across the reproductive cycle, with the exception of changes in the dendritic complexity of neurons in the medial VMHvl. This observation may be explained by the fact that not all ER1 neurons co-express PR. ER1/PR-neurons were not labeled with our genetic strategy, but their properties are sensitive to fluctuating levels of estrogen Flanagan-Cato, 2000, 2002). Intersectional strategies are needed to address this question.
Interestingly, although anterior PR1 neurons had a comparable excitability profile across the reproductive cycle, the resting potential of neurons in the receptive phase was significantly hyperpolarized. This change in resting potential likely underlies a moderate increase in the threshold to fire in these neurons, which required more current to start generating action potentials. Even if this trend did not reach significant differences with our current statistical power, the visual difference between anterior non-receptive and receptive PR1 neurons can hardly be overlooked (Fig. 3C). In addition, the anterior PR1 neurons of receptive females do not show hyperpolarized-induced rebound firing. The fact that these neurons have a more hyperpolarized resting potential, thus more distant from the firing threshold, might explain the fact that the same magnitude of Ih current fails to evoke rebound firing in anterior PR1 at the receptive phase (Fig.  5B). We also report a robust increase of the firing rate under same current injections in medial VMHvl neurons obtained from females in the non-receptive phase compared with those in the receptive phase. These changes were not accompanied by changes in capacitance or Rm across the reproductive cycle of medial neurons, thus it is unlikely that the changes in firing rate are caused by differences in the passive propagation of current into the neuron and instead points toward a modulation of voltage-sensitive channels in the medial PR1 neurons, that may be mediated by progesterone (Scharfman and MacLusky, 2006). Overall, these results suggest that the anterior and medial VMHvl have the potential to generate more action potentials in non-receptive females, which is in seeming disagreement with the enhanced-male evoked responses that we previously observed, at the population level, in the VMHvl of sexually receptive female mice (Nomoto and Lima, 2015). However, the results of this study were obtained with extracellular recordings of nonidentified neurons and therefore we do not know whether they were PR1 or PR-. Also, in the present study, we investigated the intrinsic properties of VMHvl neurons and cannot make any claim regarding the driving input they receive (synaptic and/or neuromodulatory). The VMHvl receives indirect input from the vomeronasal organ, which includes some neurons whose activity is increased in response to male stimuli when females are sexually receptive (Dey et al., 2015). It is conceivable that the enhanced VMHvl activity that we previously observed reflects modifications at the sensory level.
It has been previously reported that hormonal treatment in ovariectomized rats reduced the amount and length of secondary dendrites of a non-defined VMHvl neuronal population (Griffin and Flanagan-Cato, 2009). Even if in seeming contradiction with our findings, since we report the opposite effect in dendritic complexity and unchanged dendritic lengths, we would like to point out several possible reasons for such difference. First, the changes we observed are specific to the posterior PR1 population (and absent in our sample of PR-neurons), therefore, sampling VMHvl neurons independently of their genotype could have masked a specific plasticity process in posterior PR1 neurons. Second, in our dataset we observed that medial PR-neurons undergo a reduction in dendritic complexity that goes in line with previous reports. Thus, pooling samples from different levels of the VMHvl may mask changes, or alternatively sampling exclusively in a single AP region of the VMHvl could lead to findings that are limited to that level alone, and thus not generalizable to the whole VMHvl. Interestingly, posterior PR1 neurons did not only undergo changes at the level of their dendritic complexity, but also at the level of spine density, which was reduced in the receptive phase specifically in their LPDs. The fact that such dendritic spine changes are specific to LPDs may suggest that some pathway specific plasticity is happening during the reproductive cycle. Although it has been previously hypothesized that SPDs may preferentially integrate local excitatory inputs, while LPDs may preferentially interact with long range inputs, inhibitory neurons of the shell and neuropeptides (Griffin and Flanagan-Cato, 2011;Yamamoto et al., 2018), as of yet, no circuit mapping study has quantitatively explored what is the anatomic origin of the inputs. Thus, both the anatomic origin of those inputs and their possible relevance for the dendritic coding presented here needs to be further investigated. Our findings suggest a modification of the overall synaptic weights that the LPDs and SPDs will have in the dendritic integration for output generation, with a bias toward the inputs contacting SPDs.
It is worth noting that while the spine densities we observed are comparable to the ones reported in previous studies (Calizo and Flanagan-Cato, 2002), the dendritic lengths we report both for PR1 and for PR-neurons are several fold longer than those previously reported (Calizo and Flanagan-Cato, 2002). These changes may be explained by the differences in the technical approaches used to investigate neuronal structure. While previous reports have used fixed histologic slices (100-150 mm) filled a posteriori with lucifer yellow, in the present study we have used acute slices for ex vivo electrophysiological recordings (300 mm), in which neurons were filled with biocytin during their electrophysiological monitoring. Thus, several key differences may explain the contradictory results: first, the thickness of the slice, which could allow us to track dendrites for longer extensions; second, the fact that neurons in acute slices for ex vivo recordings reseal their membranes, producing a visible thickening in the extremity of their cut dendrites, thus allowing us to recognize incomplete dendrites with our technical approach.
In summary, we have observed multiple electrophysiological phenotypes within the VMHvl and a surprising diversity of structural and physiological plasticity processes both along the AP axis of the VMHvl, and across the reproductive cycle, within the genetically defined population of PR1 cells. These findings highlight the repertoire of local plasticity rules across the VMHvl that is probably explained by the different specialized transcriptomic clusters previously characterized across the AP extension of the VMHvl (McClellan et al., 2006;Kim et al., 2019), supporting the existence of subdivisions in the VMHvl and a more complex role of this hypothalamic structure in socio-sexual behavior, which has been recently acknowledged (Yang et al., 2013;Sakurai et al., 2016;Hashikawa et al., 2017;Inoue et al., 2019;Wang et al., 2019). The existence of multiple subdivisions that undergo local changes in response to sex hormone levels hint to the hypothesis that this hypothalamic structure might be able to control a complex and coordinated set of responses in reaction to a ubiquitous message, the systemic levels of sex hormones in circulation, from metabolism to sexual receptivity. This idea is supported by the connectivity of ER1 neurons, which present a high degree of output divergence, with anterior ER1 neurons projecting to premotor areas, while posterior ER1 are engaged in amygdalohypothalamic loops . Similar circuit motifs have been described in other hypothalamic nuclei, such as the organum vasculosum of the lamina terminalis (OVLT), where different subnuclei can coordinate a multipronged response in reaction to thirst (Graebner et al., 2015). Altogether, it becomes evident that acknowledging the different neuronal properties within the AP axis of the VMHvl is crucial to have a holistic understanding of its involvement in behavior and further studies with focal manipulations are needed to understand the functional role of observed local properties and their modulation by the reproductive cycle.