Abstract
The female reproductive cycle is gated by the circadian timing system and may be vulnerable to disruptions in the circadian system. Prior work suggests that vasoactive intestinal peptide (VIP)-expressing neurons in the suprachiasmatic nucleus (SCN) are one pathway by which the circadian clock can influence the estrous cycle, but the impact of the loss of this peptide on reproduction has not been assessed. In the present study, we first examine the impact of the genetic loss of the neuropeptide VIP on the reproductive success of female mice. Significantly, mutant females produce about half the offspring of their wild-type sisters even when mated to the same males. We also find that VIP-deficient females exhibit a disrupted estrous cycle; that is, ovulation occurs less frequently and results in the release of fewer oocytes compared with controls. Circadian rhythms of wheel-running activity are disrupted in the female mutant mice, as is the spontaneous electrical activity of dorsal SCN neurons. On a molecular level, the VIP-deficient SCN tissue exhibits lower amplitude oscillations with altered phase relationships between the SCN and peripheral oscillators as measured by PER2-driven bioluminescence. The simplest explanation of our data is that the loss of VIP results in a weakened SCN oscillator, which reduces the synchronization of the female circadian system. These results clarify one of the mechanisms by which disruption of the circadian system reduces female reproductive success.
Many people in our society experience chronic disruptions of their sleep-wake cycle. While not explicitly tested, there is every reason to suspect that many of the people suffering from disrupted sleep are also experiencing the disruption and misalignment of their circadian system. A number of lines of evidence link circadian disruption and declining reproductive success. For example, several studies have demonstrated reduced fecundity in mice with mutations in their core clock genes (Kennaway et al., 2004; Miller et al., 2004; Dolatshad et al., 2006; Alvarez et al., 2008; Ratajczak et al., 2009). Likewise, in Drosophila, clock mutants show reduced reproductive fitness (Beaver et al., 2002). Environmental disruptions of circadian timing also reduce reproductive success (Endo and Wanatabe, 1989; Summa et al., 2012), although this phenomenon is understudied. Finally, epidemiological studies have linked shift work with decreases in human reproductive success (Mahoney, 2010; Mong et al., 2011).
The female reproductive cycle is controlled by a finely tuned multiorgan feedback loop, involving hypothalamic control of pituitary hormone secretion and corresponding steroid release from the ovaries. In female rodents, ovulation occurs spontaneously in regular 4- to 5-day cycles and is gated by the circadian system (Everett and Sawyer, 1950; Turek et al., 1984; Barbacka-Surowiak et al., 2003). With the isolation or ablation of the master pacemaker of the circadian timing system, the suprachiasmatic nuclei (SCN), preovulatory hormone surges are lost in female rats (Nunez and Stephan, 1977; Gray et al., 1978; Mosko and Moore, 1979; Wiegand et al., 1980). These experimental paradigms in rodents, along with the observed deficits in clock mutants, models of deliberate circadian disruption, and epidemiological studies in humans, demonstrate that circadian disruption is negatively correlated with fertility and inspire examination of the links between the female circadian and reproductive cycles.
Previous anatomical studies have established that the connectivity from the SCN to gonadotropin-releasing hormone (GnRH) neurons in the medial preoptic area consists of both a direct connection with neurons expressing vasoactive intestinal peptide (VIP) (Van Der Beek et al., 1997; Kriegsfeld et al, 2002; Ward et al., 2009) and an indirect connection with kisspeptin- and arginine vasopressin (AVP)-positive neurons via the periventricular nucleus (De la Iglesia and Schwartz, 2006; Robertson et al., 2009; Vida et al., 2010; Williams et al., 2011). Physiologically, exogenous VIP increases the electrical activity of GnRH neurons during proestrus (Christian et al., 2005; Christian and Moenter, 2008), suggesting that VIP provides an excitatory signal from the circadian clock that helps time the GnRH surge. Transient reduction of VIP by infusing anti-sense to Vip mRNA or antibodies against VIP leads to a dampening of the magnitude of the proestrus surge in luteinizing hormone (LH) in female rats (Harney et al., 1996; Van Der Beek et al., 1999; Gerhold et al., 2005). Nonetheless, the impact of the loss of VIP on the female circadian and reproductive cycles has not been explored.
We postulate that VIP deficiency may be particularly detrimental to the SCN-GnRH circuit and may affect downstream measures of reproductive function in female mice, including fertility, regularity of the estrous cycle, and LH-induced ovulation. Additionally, the loss of VIP (VIP KO) in male mice results in severe disruption in behavioral rhythms (Colwell et al., 2003), which stems from the loss of neuronal synchrony (Aton et al., 2005; Ciarleglio et al., 2009), and also results in altered phasing of peripheral oscillators (Loh et al., 2011). The impact of VIP deficiency on female rhythms has been understudied, and we address this shortfall in our present study.
Materials and Methods
Mice
The experimental protocols used in this study were approved by the UCLA Animal Research Committee, and all recommendations for animal use and welfare, as dictated by the UCLA Division of Laboratory Animal Medicine and the guidelines from the National Institutes of Health, were followed. Vip–/– mutants (VIP–/–; also deficient in histidine isoleucine [PHI] encoded by the same locus) backcrossed to the C57BL/6 strain (Colwell et al., 2003) were crossed to PERIOD2::LUCIFERASE (PER2::LUC, Per2Luc) mice also on a C57BL/6 background (Yoo et al., 2004) as previously described (Loh et al., 2011). All wild-type (WT) and VIP–/– mice described in this study were maintained as homozygotes for the PER2::LUC knock-in fusion gene. All mice were housed under 12:12 light-dark (LD) conditions for a minimum of 2 weeks prior to any experimentation.
Measurement of Fertility
Female WT (n = 8) and VIP–/– (n = 8) littermates at 2 months of age were paired with male WT mice (2 months of age) to determine fecundity in terms of number of litters, number of pups born, and number of pups weaned. Breeders were checked daily for pregnancy and births, and dams were left undisturbed to avoid maternal stress leading to litter loss. As such, the number of pups born could only be accurately determined from P7.
Determination of Estrous Cycles by Vaginal Lavage
WT (n = 24) and VIP–/– (n = 24) female mice (2 months of age) housed under a 12:12 LD cycle were monitored using vaginal lavage at zeitgeber time (ZT) 0 to 1 (within 1 h of lights on). Care was taken to avoid insertion of the lavage pipet, and the lavage samples were spotted in small (~100 µL) drops onto glass microscope slides. Cell types contained within each vaginal smear were determined using light microscopy (10× power) and were scored for estrus state (Caligioni, 2009), where proestrus is indicated by the predominance of nucleated epithelial cells, estrus is denoted by the absence of nucleated cells and the presence of cornified squamous epithelial cells, and lavages from mice in metestrus and diestrus contain leukocytes. Mice were scored as having regular estrous cycles if they exhibited 4- to 5-day cycles throughout the monitoring period (14-21 days). Conversely, if mice exhibited more than 3 days in metestrus or diestrus, more than 2 days in proestrus, or more than 1 day in estrus or did not follow the order of progression of estrus states, they were scored as having irregular estrous cycles. The percentages of proestrus or estrus states were determined over the number of days sampled, and the percentages of proestrus smears followed by estrous were determined by visually scoring the number of proestrus smears that were immediately followed by an estrus-positive smear and dividing that by the total number of proestrus smears within each mouse.
Measurement of Serum Hormones
Serum was collected from 3-month-old WT (n = 10) and VIP–/– (n = 10) female mice at ZT 6 (midday) on the day of a proestrus smear using methods previously described (Loh et al., 2008). Serum was assayed for estradiol using radioimmunoassay by the National Peptide and Hormone Program laboratory at Harbor-UCLA (Torrance, CA). For measurement of LH concentrations, serum was collected from a separate cohort of age-matched WT (n = 4) and VIP–/– (n = 4) females (which were housed under a 12:12 LD cycle) at hourly intervals between ZT 7 (7 h after lights on) and ZT 15 (3 h after lights off) on the day of a proestrus smear. The serum samples were assayed for LH using the Milliplex MAP mouse pituitary magnetic bead assay (MPTMAG-49K, EMD Millipore, Billerica, MA). Assay sensitivity was reported as 1.7 pg/mL, with intra- and interassay %CVs reported as <15 and <20, respectively, for LH. No other analytes were run in this panel, thus ruling out cross-reactivity with other antibodies. Pituitary standard samples run in this assay contained LH ranging from 4.9 to 20 ng/mL.
Determination of Spontaneous Oocyte Release during Estrus
Vaginal smears were monitored for 2 estrous cycles for WT (n = 12) females or at least 7 days until a proestrus smear for VIP–/– (n = 12) females. During the estrous stage, mice were sacrificed at ZT 0 to 1, and the ovaries with oviducts attached were removed by cutting the ovary fat pad and the uterus. The ovary and oviduct were further dissected in Hank’s balanced salt solution (HBSS; Gibco, Invitrogen, Carlsbad, CA) to remove the fat and uterus and were transferred to fresh solution in an indented microscopic slide for analysis under a dissecting microscope. Each oviduct was carefully examined for released oocytes by making a small incision in the oviduct and teasing out the contents of the entire oviduct. The cumulus mass was dissociated by treatment in HBSS containing 300 µg/mL hyaluronidase (Sigma-Aldrich, St Louis, MO), and the oocytes were distinguished from other detritus by their spherical nature. Sampling at earlier time points (ZT 12-23 on the day of a proestrus smear) revealed an absence of oocytes in both WT (n = 5) and VIP–/– (n = 3) oviducts.
Oocyte Release in Response to Ovarian Stimulation
In a separate cohort, ovulation was induced using a combination of pregnant mare serum gonadotropin (PMSG) and human chorionic gonadotropin (hCG). Female WT (n = 6) and VIP–/– (n = 6) mice (2 months of age) were hormone-primed with an intraperitoneal (ip) injection of 5 IU of PMSG (G4877, Sigma-Aldrich) applied at ZT 11 on the day of metestrus, and ovulation was induced 48 h later by ip injection of 5 IU of hCG (Novarel, Ferring Pharmaceuticals, Parsippany, NJ). On the morning following the hCG injection, mice were sacrificed between ZT 0 and 1, and the number of oocytes released into the oviducts was counted.
Wheel-Running Activity
Adult female mice (2-4 months of age; WT n = 10, VIP–/– n = 10) were singly housed in cages containing wheels, and wheel-running activity was recorded as previously described (Colwell et al., 2003). Mice were entrained to a 12:12 h LD cycle for a minimum of 2 weeks prior to collection of 10 to 14 days of data under LD conditions, followed by 10 to 14 days in constant darkness (DD) to obtain free-running activity. Data were recorded using Mini Mitter (Bend, OR) data loggers in 3-min bins, and 10 days of data under each condition were averaged for analysis. Free-running period (tau) was determined using the χ2 periodogram, and the power of the rhythm was calculated by multiplying the Qp by 100/n, where n = number of data points examined using the El Temps program (Barcelona, Spain). Activity duration (alpha) was determined by the duration of activity over the threshold of the mean using an average waveform of 10 days of activity. The phase of activity onset was also defined as the time at which activity crossed the threshold of the mean of total activity. To calculate the phase angle of entrainment revealed by release into DD, the difference on the day of lighting change was determined between best-fit regression lines drawn through activity onset prior to and after release into DD. Precision was determined by calculating the daily variation in onset from a best-fit regression line drawn through 10 days of activity in both LD and DD conditions using the Clocklab (Actimetrics, Wilmette, IL) program. Fragmentation was defined by bouts per day with a maximum gap of 21 min.
Measurement of Spontaneous Firing Rate in SCN
Whole-cell patch-clamp electrophysiology methods used are similar to those described previously (Kudo et al., 2013). Young female mice (2-3 months old) were anesthetized in an isoflurane chamber in the day (ZT 2.5) or in the night (ZT 12.5). Then 300-µm coronal sections were prepared and superfused continuously with room temperature artificial cerebral spinal fluid (ACSF) aerated with 95% O2/5% CO2 in a recording chamber (PH-1, Warner Instruments, Hamden, CT) attached to the stage of a fixed-stage upright DIC microscope (OLYMPUS, Tokyo, Japan). Whole-cell patch-clamp SCN recordings were made using electrode micropipettes (3-7 MΩ), and recording electrodes were filled with standard internal solution (in mM): K-gluconate, 112.5; EGTA, 1; HEPES, 10; MgATP, 5; and MgCl2, 1. Internal solution pH was adjusted to 7.25 to 7.3, and osmolarity was adjusted to 290 to 300 mOsm. Recordings were obtained with the AXOPATCH 200B amplifier (Molecular Devices, Sunnyvale, CA) and monitored online with pCLAMP (Version 10, Molecular Devices).
For these experiments, recordings were made from cells located in the dorsal SCN (dSCN) as defined by being dorsal to the tip of the third ventricle. The access resistance of these cells ranged from 15 to 35 MΩ, and cell capacitance typically ranged from 6 to 18 pF. Data were not collected if access resistance values changed significantly (>20%) during the course of the experiment and/or if a neuron was not able to decrease firing rate following excitatory treatment with N-methyl-
Monitoring of PER2::LUC Bioluminescence
Female WT (n = 7) and VIP–/– (n = 10) mice (2-3 months old) were housed under a 12:12 LD cycle for at least 2 weeks prior to sampling. Vaginal smears were performed during the second week to determine estrous state, and mice in proestrus (nucleated cells) were sampled at ZT 10 to 11. SCN explants were dissected from 300-µm coronal sections, dissected ovaries were further halved, and 1- to 2-mm3 explants were dissected from the pituitary and midpoint of the uterus. PER2::LUC bioluminescence recordings using a Lumicycle photometer (Actimetrics, Wilmette, IL) were conducted and analyzed as previously described (Loh et al., 2011; Kuljis et al., 2013). Briefly, bioluminescence readings were taken every 10 min and, for analysis, were baseline-subtracted, detrended (24-h window), and smoothed (2-h window). Due to this process, we do not report the first recorded peak but instead report the first calculated peak for phase of PER2::LUC expression relative to the prior LD cycle (ZT). Amplitude of each cycle was determined by the sum of the peak and subsequent trough values, and period was determined from a minimum of 4 consecutive cycles. Damping rate was determined by the gradient of the best-fit line through a minimum of 4 consecutive cycles, and R2 is the coefficient of determination of the best-fit line, indicating the goodness of fit.
Statistical Analysis
For comparison of female WT and VIP–/– measures that passed normality and equal variance tests, we applied Student’s t tests and deemed differences as significant if p < 0.05. For parameters that failed either normality or equal variance tests, rank sum t tests were applied. Two-factor repeated measures analysis of variance (ANOVA) was applied to examine the relative contributions of genotype and time of day of proestrus serum LH concentration. Post hoc Bonferroni comparisons were performed to compare the effect of genotype at each sampling time. Values are reported as mean ± standard error of the mean (SEM).
Results
Female Reproduction Is Greatly Reduced in VIP-Deficient Mice
To determine whether VIP is critical for reproduction, we paired naïve WT and VIP–/– females littermates with WT males. VIP–/– females produced dramatically fewer pups (10.6 ± 1.5) than their sisters (21.5 ± 3.0, p < 0.01) during the 6 months of breeding, demonstrating reduced fecundity. The mean number of pups per litter was not significantly different between the genotypes (WT 4.8 ± 1.7, VIP–/– 4.0 ± 1.6 pups per litter; p = 0.32); however, reduced fecundity is accounted for by the fewer litters born to VIP–/– dams (2.9 ± 0.3) compared with WT females (4.5 ± 0.3, p < 0.01). Furthermore, VIP KO females spent a smaller percentage of time either pregnant or nursing (46.8% ± 3.1%) versus WT females (76.9% ± 5.2%, p < 0.001). We can conclude from these findings that although VIP–/– females are capable of reproduction, they are sub-fertile, spending less time pregnant or nursing and producing half the number of offspring compared with their sisters held under the same conditions.
Regularity of the Estrous Cycle Is Greatly Reduced in VIP-Deficient Mice
To determine whether the basic reproductive biology underlying ovulation is affected by VIP deficiency, we examined the estrous cycle by daily vaginal lavage in WT and VIP–/– female littermates (representative examples shown in Fig. 1A). We found irregular estrous cycles in the majority of VIP–/– females (20/24), while the majority of WT females exhibited regular cycling through proestrus, estrus, metestrus, and diestrus (19/24). Cycle duration as determined by the time between estrus was found to be longer in VIP–/– females (8.18 ± 1.05 days) compared with typical WT cycles (4.64 ± 0.17 days; p < 0.001; Fig. 1B). A common finding in VIP–/– females was that not every proestrus smear was followed by an estrus smear (51.6% ± 6.5%; Fig. 1A, C), whereas almost every proestrus smear in WT females was immediately followed by estrus (84.4% ± 5.5%, p < 0.001; typical example depicted in Fig. 1A, top). Additionally, VIP–/– females exhibited elongated duration in diestrus/metestrus as indicated by the presence of leukocytic cells (57.69% ± 3.88%) compared with WT females (44.68% ± 2.06%, p = 0.008; Fig. 1D). Over the entire period during which estrous cycles were assessed, VIP–/– females had significantly fewer percentages of days in proestrus (15.9% ± 1.8%) compared with WT females (20.7% ± 1.2%, p < 0.01) and also had fewer days in estrus (14.7% ± 2.6% vs. WT 20.5% ± 1.6%, p < 0.05). Our data indicate that VIP is essential for regular estrous cycles and that VIP deficiency is characterized by long spells in the diestrus and metestrus states. Furthermore, estrus is not guaranteed after proestrus, which further adds to the longer duration between estrous states in mutant females.

Impact of VIP deficiency on estrous cycles. (A) VIP deficiency results in disrupted estrous cycles (middle and bottom) compared with the 4- to 5-day estrous cycles exhibited by WT females (top). P = proestrus, majority nucleated epithelial cells. E = estrus, cornified cells. D/M = diestrus or metestrus, presence of leukocytic cells. (B) The average length of an estrous-to-estrous cycle is longer in VIP-deficient mice than in WT females (*p < 0.05). (C) VIP deficiency reduces the number of proestrus to estrus transitions. (D) The number of consecutive days in diestrus and metestrus increases with the loss of VIP (*p < 0.05).
Serum Estradiol Levels Are Not Altered by the Loss of VIP
To determine whether the steroidal signals that define the ovulatory process are affected by the VIP mutation, we measured serum concentrations of estradiol during proestrus. The increase in circulating estradiol during the day of proestrus precedes the activation of GnRH neurons, which initiates the LH surge and in turn stimulates ovulation (Legan et al., 1975). The midday (ZT 6) serum concentration of estradiol was no different in VIP–/– females (52.9 ± 4.2 pg/µL) and WT females (54.4 ± 3.8 pg/µL, p = 0.77; Fig. 2A). We also performed hourly measurements of LH from ZT 7 to ZT 15 to determine whether surge time and concentration were altered in VIP–/– females. The maximum LH measured in VIP–/– females (0.2 ± 0.04 pg/mL) was not significantly different from that found in WT females (1.1 ± 0.8 ng/mL; rank sum test, p = 0.11; Suppl. Fig. S1). The maximum LH concentration in WT females was found to be lower than surge-like concentrations of LH (>5 ng/mL) despite the hourly sampling on the day of proestrus. Given this finding, coupled with the variability in LH concentrations within WT females, we are unable to draw conclusions about genotypic differences in LH surge concentration. We ascertained that the temporal window of peak LH concentration within animals was not significantly altered by the loss of VIP (WT ZT 11.3 ± 1.4, VIP–/– ZT 11.4 ± 1.3; rank sum test p = 0.73; Suppl. Fig. S1). Our findings indicate that the midday increase in estradiol concentration during proestrus is not altered in the VIP–/– mice, which do not exhibit any measurable changes in LH concentration or peak timing.

Impact of VIP deficiency on circulating estradiol and ovulation in adult females. (A) The proestrus midday increase in serum estradiol is not affected by the loss of VIP in female mice. (B) The number of oocytes spontaneously released on the morning of estrus is reduced in VIP–/– oviducts compared with WT (*p = 0.002). (C) The number of oocytes released in response to PMSG-hCG ovarian stimulation technique in hormone-primed female mice is not different between the genotypes.
Spontaneous Ovulation Is Dependent on VIP
Mice ovulate spontaneously following estrus, and we examined the number of oocytes released in both WT and VIP–/– mice on the day following a proestrus smear. We found that 6 of 12 VIP–/– females did not release any oocytes after a proestrus smear compared with 2 of 12 WT females. Of the mice that released oocytes, significantly fewer oocytes were found in VIP–/– oviducts (5.3 ± 0.7 vs. WT 10.7 ± 0.7 oocytes/pair of oviducts, p < 0.001; Fig. 2B). To determine whether the reduced number of oocytes in VIP–/– females during ovulation is due to an inability of the ovaries to release the oocytes, we hormone-primed female mice using the PMSG-hCG ovarian stimulation technique commonly used to increase fertility in the production of transgenic mice. Following ovarian stimulation, we found that VIP–/– mice release 26.7 ± 1.7 oocytes, a number that is not significantly different from the 24.7 ± 2.8 oocytes released by WT females (p = 0.51; Fig. 2C). In contrast to our finding that 50% of VIP–/– females do not spontaneously ovulate on the day after proestrus, the entire cohort of hormone-primed VIP–/– females released oocytes, demonstrating that the VIP–/– ovaries are able to respond to gonadotropin stimulation.
Circadian Behavior Is Disrupted in Female VIP-Deficient Mice
Several studies have characterized the circadian phenotype in male VIP-deficient mice (Colwell et al., 2003; Aton et al., 2005; Ciarleglio et al., 2009) while surprisingly neglecting to describe the behavior of the female mutants. We used wheel-running activity to determine the impact of the loss of VIP on diurnal and circadian rhythms of behavior in female mice (Fig. 3). In LD conditions, the diurnal rhythms are largely similar between the genotypes, with the only significant change seen in the extent of fragmentation of the nocturnal activity bout (Table 1). It is in DD that the full extent of the circadian phenotype in these mutant animals emerges. We observed significant decreases in the power and precision of the wheel-running activity. Additionally, there is a significant decrease in the free-running period and a dramatic decrease in the magnitude of light-induced phase shifts to light exposure at Circadian Time (CT) 16 (Table 1). Overall activity levels are not altered, highlighting the general conclusion that the loss of VIP does not produce a gross motor phenotype. A characteristic feature of the behavior of VIP–/– mice is the vastly advanced angle of entrainment to the LD cycle, as revealed when the mice are placed in DD. Like their male counterparts, VIP–/– females are abnormally entrained to the LD cycle, but at a significantly smaller magnitude, suggesting a sex-dependent effect of the mutation (Table 2). It is worth noting that there was no problem measuring key circadian parameters in WT female mice despite their likely estrous cycling. Overall, the female VIP-deficient mice showed clear and dramatic disruption of their circadian rhythms in wheel-running activity.

Representative actograms and periodograms of female WT (A, C) and VIP–/– (B, D) mice under LD and DD conditions. (A, C) Wheel-running activity records are double plotted, and gray shading on the actograms denotes lights-off. (B, D) Plot of power of activity as a function of period (%V) under DD. The diagonal line represents 0.1% significance level.
Parameters of wheel-running activity in female WT and VIP–/– mice.
LD = light-dark cycle; DD = constant darkness; rev = wheel revolutions; %V = normalized Qp; CT = Circadian Time. Values are reported as mean ± SEM. n = 10 for WT and n = 10 for VIP–/– females.
p < 0.05, **p < 0.01, and ***p < 0.001 indicate significant differences between WT and VIP–/– mice.
Sex-dependent effects of VIP deficiency on wheel-running activity.
LD = light-dark cycle; DD = constant darkness. Wheel-running activity was monitored in male VIP–/– littermates from the same colony of VIP;PER2::LUC double transgenics. As previously described (Colwell et al., 2003), male VIP–/–;PER2::LUC mice have a shorter free-running period, vastly altered angle of entrainment as indicated by the large phase advance upon release to constant darkness, greatly reduced power of rhythm, and poor precision of daily activity onset. The same differences are observed in female VIP–/–;PER2::LUC mice. However, the angle of entrainment is less affected by the loss of VIP in females than males (2-way ANOVA sex × genotype interaction F = 11.22, p = 0.002; post hoc Holm-Sidak t = 7.69, p < 0.001). In contrast, female VIP–/– mice have more imprecise activity onset than males (2-way ANOVA sex × genotype interaction F = 4.11, p = 0.05; post hoc Holm-Sidak t = 3.19, p = 0.003).
Significant differences between WT and VIP–/– mice within sex as determined by Student’s t test.
Significant differences between male and female VIP–/– mice as determined by ANOVA post hoc tests.
p < 0.05, **p < 0.01, and ***p < 0.001 indicate significant differences between WT and VIP -/ - mice.
SCN Neural Activity Is Compromised in Female VIP-Deficient Mice
SCN neurons are spontaneously electrically- active cells that generate rhythms in action potential frequency with peak activity during the day (Colwell, 2011). We examined daytime (ZT 4-6) and nighttime (ZT 14-16) neural activity of dSCN neurons from WT and VIP–/– mice using whole-cell patch-clamp electrophysiology. As expected, female WT mice exhibited a day-night difference in SFR, with higher frequency of firing during the day (3.57 ± 0.56 Hz) than during the night (1.25 ± 0.44, p = 0.011; Fig. 4A). In comparison, female VIP–/– dSCN neurons exhibited lower firing rates during the day (1.89 ± 0.33 Hz, vs. WT, p = 0.019), and there was no longer a significant day-night difference in SFR (night SFR 3.74 ± 1.08 Hz, vs. day SFR, p = 0.054; Fig. 4B). When the cells were synaptically isolated using a GABA-blocker (GabZ), VIP–/– dSCN neurons were again found to have lost the day-night differences in SFR (day 2.61 ± 0.56 Hz vs. night 3.29 ± 1.10 Hz; p = 0.54 Fig. 4C). The reduced daytime SFR in the dSCN in female VIP –/– mice demonstrates that the mutants exhibit a weakened output from the SCN circuit in vitro and suggests that the rest of the circadian system receives a low-amplitude signal from the SCN in vivo.

Reduced amplitude in rhythm of spontaneous firing in VIP–/– female SCN neurons. (A, B) Representative traces of spontaneous action potentials in day (left) and night (right) in WT (A) and VIP–/– (B) dorsal SCN neurons. (C) Population averages of firing rate (Hz) across WT (day n = 11, night n = 9) and VIP–/– (day n = 14, night n = 9) SCN neurons synaptically isolated by treatment with GabZ. Student’s t tests determined a significant difference between daytime and nighttime firing rates in WT (*p < 0.05) but not in VIP–/– SCN. VIP–/– daytime firing rates were also found to be significantly different from WT (^p < 0.05).
PER2 Rhythms in SCN, Pituitary, and Reproductive Organs Are Differentially Affected in Female VIP-Deficient Mice
The circadian timing system consists of a central clock located in the SCN and peripheral clocks found throughout the body (Albrecht, 2012; Mohawk et al., 2012). To determine the impact of the loss of VIP on the molecular clockwork, we measured the phase of the expression of a critical circadian gene, Period2, using the PER2::LUC reporter, which drives expression of luciferase fused to the PERIOD2 protein in double transgenics for the VIP mutation and the reporter (Loh et al., 2011). To control for the effects of estrous cycle-driven hormones, we only used proestrus female mice. While the period and phase of peak PER2::LUC bioluminescence in the SCN were not affected by loss of VIP (Fig. 5A, Table 3), the phasing of the pituitary and uterus were advanced in VIP–/– PER2::LUC explants (Fig. 5B). Conversely, the amplitude of the oscillations in the peripheral explants were unaffected by loss of VIP (Fig. 5C), whereas the amplitude of PER2::LUC bioluminescence in the SCN and pituitary were blunted in the mutant mice (Fig. 5A, C). The loss of VIP did not affect the period of PER2::LUC rhythms in any of the explants (Table 3). We can thus conclude that VIP is not needed for intrinsic rhythmicity of peripheral oscillators in the ovary or uterus, since the period and amplitude of these rhythms are not different from WT controls. In contrast, VIP is important for the amplitude of PER2::LUC bioluminescence in the SCN and the pituitary and is also critical for the normal phasing of the pituitary and uterus.

Impact of the loss of VIP on the amplitude and phasing of PER2::LUC bioluminescence in female SCN and peripheral organs. (A) Representative bioluminescence rhythms from WT and VIP–/– SCN ex vivo culture. (B) Phases of peak expression of PER2::LUC of the pituitary and the uterus are altered in the VIP-deficient mice when compared with WT controls (*p < 0.05). (C
PER2::LUC bioluminescence parameters of WT and VIP–/– ex vivo explants.
Values are reported as mean ± SEM. n = 7 for WT and n = 10 for VIP–/–.
p < 0.05, **p < 0.01, and ***p < 0.001 indicate significant differences between WT and VIP–/– mice as determined by Student’s t tests.
Discussion
We wished to address the impact of VIP deficiency on the reproductive success of female mice, due to the importance of both VIP and the circadian timing system on reproduction. Anecdotal evidence from our colony indicated difficulty maintaining the homozygote VIP-deficient line. Careful measurements over 6 months confirmed that the VIP–/– females are subfertile, spending less time pregnant or nursing and producing half the number of offspring compared with their sisters held under the same conditions. Low fecundity has also been reported in VIPR2 knock-out mice (Dolatshad et al., 2006). These studies indicate that the loss of VIP signaling can result in the decline of reproduction, and they raise questions about the underlying mechanisms.
Previous studies have established that VIP-expressing neurons provide one of the connections between the SCN and GnRH neurons in the medial preoptic area (Van Der Beek et al., 1997; Horvath et al., 1998; Kriegsfeld et al., 2002; Ward et al., 2009; Vida et al., 2010). One of the receptors for VIP (VIPR2) is expressed on GnRH neurons (Smith et al., 2000), and, physiologically, the exogenous treatment of VIP increases the electrical activity of GnRH neurons whereas blocking the receptor decreases firing (Christian and Moenter, 2008). Therefore, we examined the impact of the loss of VIP on the estrous cycle. We found that VIP-deficient female mice have irregular and lengthened estrous cycles, due to a failure to progress from proestrus to estrus (Fig. 1). The VIPR2 KO mice also show a disrupted estrous cycle, but these mice spent more time in estrus (Dolatshad et al., 2006). We cannot explain the differences between the mutant lines but point out that VIPR2 also binds PACAP. Both lines of mutants are able to reproduce, and the AVP and kisspeptin pathway most likely allows for continued, albeit lower, fertility. Current thinking places greater importance on the indirect multisynaptic connection between the SCN and GnRH neurons via AVP and kisspeptin (Robertson et al., 2009), which may be timed or gated by the direct VIPergic connection (Williams et al., 2011). In rats, a sophisticated model of internal desynchrony within the SCN provides evidence that the AVP- and VIP-expressing subdivisions need to have coupled cross-talk to achieve maximal GnRH-dependent LH surge activity (Smarr et al., 2012). Ovulation involves the coincidence between high estrogen levels sensitizing the gonadotrophs with a gating signal from the SCN controlling the time of day of ovulation (Legan and Karsch, 1975; Christian and Moenter, 2010; Williams and Krigsfeld, 2012). We predicted a blunted LH surge in the mutant females as has been seen in rats treated with anti-sense against Vip (Harney et al., 1996; Gerhold et al., 2005) but were unable to successfully measure the surge with hourly sampling. We were able to determine that the necessary increase in serum estradiol is not affected by the loss of VIP (Fig. 2). The end measure of the estrous cycle in females is ultimately ovulation, which we found to be deficient in VIP–/– females, in which oocytes were released less frequently and in fewer numbers than in WT controls. VIP has been shown to mediate the development and survival of female rat ovarian follicles (Flaws et al., 1995; Cecconi et al., 2004; Chen et al., 2013) and has a steroidogenic role in the ovaries (Ahmed et al, 1986; Kowalewski et al., 2010). In our study, ovarian stimulation in response to exogenous treatment with gonadotropins did not vary between the genotypes (Fig. 2), demonstrating that the deficit in the mutants is not in the oocyte content of the ovaries but rather in the luteinizing signal. Additionally, the comparable proestrus concentration of serum estradiol in VIP-competent and VIP–/– females suggests that VIP is not critical for ovarian steroidogenesis. It is worth noting that prior studies examining the effect of VIP on ovarian development, survival, and steroidogenesis used in vitro manipulations, supplementing cultures with exogenous VIP. In contrast, development in our VIP–/– model may have used redundant pathways (e.g., PACAP) (Cecconi et al., 2004), allowing for reproduction to continue, albeit at suboptimal levels. The disruption in the estrous cycle provides a proximate explanation for the observed decline in reproductive success and is likely a common consequence of circadian disruption.
The circadian timing system of female VIP-deficient mice is compromised at several levels. Behaviorally, female mutant mice have severe deficits in circadian rhythms of locomotor activity, characterized by a shortened free-running period, low power of rhythms, and poor precision of cycle-to-cycle activity onset (Fig. 3, Table 1). A further characteristic observed in both male and female VIP–/– mice is the abnormally advanced angle of entrainment revealed when mice are released from LD to free-running DD conditions. This suggests that the mice are abnormally entrained to the prior LD cycle, and we and others have reported phase-advanced physiology and behavior in VIP–/– and VIPR2–/– males under an LD cycle (Colwell et al., 2003; Aton et al., 2005; Ciarleglio et al., 2009). Although the VIP–/– female mice also exhibit an advanced phase angle of entrainment, its magnitude is half that of their male counterparts (Table 2), suggesting there is a sex-dependent effect of the mutation. Additional sex-dependent effects include the reduced amount of activity specific to the VIP–/– males, which may be mediated by testosterone, also reduced in VIP–/– males (Lacombe et al., 2007). Physiologically, our data indicate that the VIP–/– females have severe deficits in the expression of daily rhythms in spontaneous electrical activity (Fig. 4). VIP is critical for coupling SCN neurons, as mice deficient for the peptide or its receptor lose coherence and synchrony at the level of electrical firing (Aton et al., 2005; Brown et al., 2007) and gene expression (Harmar et al., 2002; Maywood et al., 2006; Hughes et al., 2008; Loh et al., 2008; Ciarleglio et al., 2009; Dragich et al., 2010; Loh et al., 2011). Our finding that daytime SFR of SCN neurons is reduced in VIP-deficient females even in the presence of the GABA-blocker GabZ suggests that even when GABA-dependent neuronal cross-talk is disabled acutely, the deficits remain. This is consistent with earlier work indicating that VIP signaling is a critical regulator of intrinsic membrane events in SCN neurons (Pakhotin et al., 2006; Kudo et al., 2013). This in vitro electrophysiological data coupled with the behavioral analysis indicate a weakly rhythmic SCN in the female mutant mice.
Finally, we assessed the molecular circadian clockwork, measuring PER2-driven bioluminescence in female VIP-deficient mice (Fig. 5, Table 3). At the level of the SCN, the PER2::LUC rhythms showed reduced amplitude but a normal phase relationship to the prior LD cycle. With the exception of the pituitary, the peripheral clocks mostly showed no change in amplitude but exhibited changes in phasing relative to each other and to the prior LD cycle. These alterations in the timing of gene expression in the peripheral oscillators were similar to those previously found in male VIP-deficient mice (Loh et al., 2011). Curiously, although the loss of VIP led to a phase advance of most of the peripheral oscillators (Table 3), the phasing of VIP–/– ovaries was not different from WT (Fig. 5), suggesting that VIP is not necessary for the circadian rhythm of PER2 expression in the ovary (Sellix et al., 2010) and that the temporal gating of ovarian response to gonadotropins is unlikely to underlie the ovulation deficits. The circadian phase of ovaries can be altered by gonadotropins (Yoshikawa et al., 2009) and is also dependent on estrous state (Nakamura et al., 2008). By fixing our collection of PER2::LUC organs to proestrus, we may have picked an estrous state in which the phase of the circadian oscillator in the ovary is governed by circulating steroids. It is also important to note that the phase advance of the other peripheral oscillators is not due to a change in the intrinsic period of PER2 expression, which remains unaltered between the WT and VIP–/– tissues. Again, the gene expression data are best explained by a weak SCN clock whose output is less able to synchronize the network of circadian oscillators. These 3 lines of evidence for a weakened SCN clock suggest a severe impact on the strength of SCN output signals. These include a dampening of the indirect AVP-kisspeptin route to the GnRH neurons and the loss of the direct VIP-GnRH relay, both of which ultimately result in a gross decrease in reproductive success.
One limitation of our study is that using a whole-body knock-out of VIP does not allow us to determine which of these deficits has the strongest impact on the observed sub-fertility of VIP KO females. The tripartite functions of VIP as a synchronizer between SCN neurons, its role as an SCN output signal to GnRH neurons, and its local role in ovarian follicular development cannot be easily distinguished in our model. Ideally, to distinguish the degree to which these varying functions of VIP affect reproduction, one would use gene ablation targeted to specific organs or cell types, but a floxed-VIP mutation does not exist at present. Alternatively, one could use targeted rescue of VIP loss, for instance, using viral constructs to restore VIP expression in SCN neurons. Nevertheless, our study reinforces the importance of VIP in female reproduction.
For the purposes of this study, the VIP-deficient female mice served as a model for circadian disruption caused by deficits in cellular communication rather than direct disruption of the molecular clockwork. In this sense, the present study brings the females into the story being generated by a number of laboratories working with VIP- and VIPR2-deficient mice. Although in general, the circadian phenotype observed with female mutants was similar to that gathered from males, there were also some interesting differences that highlight the importance of including both sexes in our preclinical work (Clayton and Collins, 2014). In addition, there is good reason to be specifically interested in the role of VIP in the circadian regulation of female reproduction. VIP-expressing neurons couple the SCN with GnRH neurons and alter their physiology. In particular, as has been previously described (Duncan, 2006; Downs and Wise, 2009; Christian and Moenter, 2010), a decline in VIP has particular relevance to aging in women, where the issue of menopause and its corresponding index of symptoms have negative consequences on the quality of life for middle-aged women. In female rats, the beginnings of reproductive senescence are heralded by a decline in Vip expression in middle-age (Krajnak et al., 1998) that is correlated with a decline in the sensitivity of VIP-innervated GnRH neurons (Krajnak et al., 2001). When middle-aged female rats are treated with intracerebroventricular injections of VIP, the LH surge is restored to levels similar to those of young females (Sun et al., 2012). These findings suggest that reproductive aging may be due in part to hypothalamic changes as well as gonadal decline. Our finding that the absence of VIP has a greater impact on the SCN than the ovaries supports the view that reproductive aging has a hypothalamic cause, and this finding provides at least one mechanism (see Sellix, 2013, for a discussion of other mechanisms) by which circadian disruption can cause a decline in reproductive function.
Footnotes
Acknowledgements
The authors acknowledge the contributions of the following people. Dr. Analyne Schroeder, Tamara Cutler, Dr. Takashi Kudo, and Dr. Takahiro Nakamura helped with dissections for bioluminescence organotypic cultures. Donna Crandall prepared the figures for submission. Dr. A. Parlow’s laboratory at UCLA Harbor, Long Beach, CA, performed the radioimmunoassay for estradiol. DHL and CSC received funding support from the Iris-Cantor Women’s Health Center at UCLA. DAK was supported by a NIH training grant from the Laboratory of Neuroendocrinology (T32 HD 07228-26).
Author Contributions
Designed study: DHL, DAK, and CSC. Draft of manuscript: DHL, DAK, and CSC. Behavior: DHL, LA, and DT. Electrophysiology: DAK. PER2::LUC bioluminescence: DHL. Estrous and oocyte counts: DHL, LA, and YW. Hormones: HBW and DHL.
Conflict of Interest Statement
The author(s) have no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Notes
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
