Abstract

Context

Studies suggest that female reproductive hormones are under circadian regulation, although methodological differences have led to inconsistent findings.

Objective

To determine whether circulating levels of reproductive hormones exhibit circadian rhythms.

Design

Blood samples were collected across ∼90 consecutive hours, including 2 baseline days under a standard sleep-wake schedule and ∼50 hours of extended wake under constant routine (CR) conditions.

Setting

Intensive Physiological Monitoring Unit, Brigham and Women’s Hospital.

Participants

Seventeen healthy premenopausal women (22.8 ± 2.6 years; nine follicular; eight luteal).

Interventions

Fifty-hour CR.

Main Outcome Measures

Plasma estradiol (E2), progesterone (P4), LH, FSH, SHBG, melatonin, and core body temperature.

Results

All hormones exhibited significant 24-hour rhythms under both standard sleep-wake and CR conditions during the follicular phase (P < 0.05). In contrast, only FSH and SHBG were significantly rhythmic during the luteal phase. Rhythm acrophases and amplitudes were similar between standard sleep-wake and CR conditions. The acrophase occurred in the morning for P4; in the afternoon for FSH, LH, and SHBG; and during the night for E2.

Conclusions

Our results confirm previous reports of ∼24-hour rhythms in many female reproductive hormones in humans under ambulatory conditions but demonstrate that these hormones are under endogenous circadian regulation, defined as persisting in the absence of external time cues. These results may have important implications for the effects of circadian disruption on reproductive function.

Circadian regulation of female reproductive physiology is well characterized in rodent models [reviewed in Boden and Kennaway (1)] but is mainly supported by indirect evidence in humans. A recent meta-analysis of 16 female cohorts (123,403 women) showed that women who are engaged in shift work, which can disrupt circadian rhythms, have increased odds of menstrual cycle disruption, early spontaneous pregnancy loss, and infertility by ∼22%, ∼30%, and ∼80%, respectively, compared with women not engaged in shift work (2). Hormones related to these functional outcomes have been reported to exhibit 24-hour rhythms under normal sleep-wake, meal, and lighting conditions. Diurnal rhythms in estradiol (E2), progesterone (P4), LH, and FSH have been reported in premenopausal women (310). Findings of diurnal variability in these hormones have not been consistent, however (5, 1113). Differences in the study population (e.g., the menstrual phase at which women were studied) and analytic techniques [e.g., the application of nonlinear cosinor regression in some (46, 912) but not others (3, 7, 8, 13)] likely contributed to the inconsistent findings between studies in the identification of rhythmic hormones.

Moreover, these studies cannot confirm whether the daily rhythm represents the output of an endogenous circadian pacemaker or is a result of evoked responses, such as sleep-wake and feeding-fasting cycles, or a combination of both (14). The presence of an endogenous circadian rhythm can only be confirmed when external time cues are absent or distributed uniformly, such as during a constant routine (CR) protocol, where activity, posture, light, food intake, and wakefulness remain constant (14). When healthy, premenopausal women were studied under such CR conditions, LH, FSH, and E2 levels were reported to be arrhythmic (15) and therefore, not under endogenous circadian regulation. This study included only women in the early follicular phase of their menstrual cycle (approximately cycle-day 4 in relation to menses), however, which precludes examining whether circadian regulation of these hormones differs between the follicular and luteal phase. Moreover, these studies did not directly compare rhythms under constant and ambulatory conditions, which is important to elucidate the masking effects of behavioral events, such as sleep and feeding-fasting on the underlying circadian regulation of physiologic rhythms (16). Therefore, in the current study, we assessed the endogenous circadian regulation of E2, P4, FSH, LH, and SHBG in the follicular and luteal phases of the menstrual cycle and compared the rhythms in these hormones between CR and ambulatory conditions.

Methods

Participants

We studied 17 healthy premenopausal women (mean age ± SD = 22.8 ± 2.6 years; range 19 to 29 years; nine women in the follicular phase). All participants underwent comprehensive physical, psychiatric, and medical screening. Participants self-reported a regular menstrual cycle lasting 27 to 35 days (mean cycle length ± SD = 28.9 ± 1.84 days) and were not using oral contraception at the time of the study or within the previous 3 months. Menstrual phase was determined based on each individual participant’s daily average P4 levels. Women with daily mean P4 concentrations >3 ng/mL were considered to be postovulatory (17, 18) and assigned to the luteal phase, and women with P4 <3 ng/mL were assigned to the follicular phase. Given the protocol duration, one participant (2251V) changed from the follicular phase during baseline to the luteal phase during CR based on levels of P4, with concomitant changes in FSH and LH, indicating ovulation. Demographic information, P4 concentrations, and allocation of menstrual phase for each of the participants on the different study days are presented in Table 1.

Table 1.

Participant Demographic Information

IDAge, yMenst. Length, dBed Time, h:minDLMO, h:minBaselineCR
P4, ng/mL M (SD)Menst. PhaseP4, ng/mL M (SD)Menst. Phase
22JSV202800:3023:272.10 (0.54)F2.13 (0.31)F
22K7V193023:5823:062.06 (0.42)F2.04 (0.26)F
25N6V232822:4020:512.64 (0.77)F2.61 (0.48)F
2614V222800:0123:292.33 (0.44)F2.33 (0.42)F
2692V232823:3322:581.90 (0.33)F2.18 (0.46)F
26F2V223300:040:071.75 (0.64)F2.07 (0.45)F
26G3V242923:0322:362.48 (0.60)F2.22 (0.35)F
26G6V232921:0021:462.52 (0.83)F1.96 (0.37)F
26H6V1931.523:0423:342.35 (0.42)F
2251V232801:320:572.69 (0.77)F4.82 (1.17)L
21B8V272823:1221:0239.44 (10.58)L16.32 (6.51)L
22A1V242922:0121:554.70 (1.46)L10.74 (3.13)L
22K3V242722:0312.11 (3.68)L4.10 (1.28)L
25Q2V232723:0023:1338.92 (13.79)L28.52 (4.79)L
2622V222722:4121:0614.51 (5.01)L5.20 (3.14)L
26P1V292822:5121:3216.61 (4.68)L5.70 (3.06)L
26R1V2032.521:5619:1013.44 (3.51)L9.67 (2.34)L
Follicular mean ± SEM21.67 ± 0.6229.39 ± 0.6023:19 ± 00:2122:53 ± 00:202.27 ± 0.11n = 92.21 ± 0.07n = 9
Luteal mean ± SEM24.00 ± 1.0028.31 ± 0.6522:54 ± 00:2421:51 ± 00:4119.96 ± 5.16n = 710.63 ± 2.94n = 8
Group mean ± SEM22.76 ± 0.6328.88 ± 0.4523:07 ± 00:1522:26 ± 00:2210.01 ± 3.13n = 166.17 ± 1.70n = 17
IDAge, yMenst. Length, dBed Time, h:minDLMO, h:minBaselineCR
P4, ng/mL M (SD)Menst. PhaseP4, ng/mL M (SD)Menst. Phase
22JSV202800:3023:272.10 (0.54)F2.13 (0.31)F
22K7V193023:5823:062.06 (0.42)F2.04 (0.26)F
25N6V232822:4020:512.64 (0.77)F2.61 (0.48)F
2614V222800:0123:292.33 (0.44)F2.33 (0.42)F
2692V232823:3322:581.90 (0.33)F2.18 (0.46)F
26F2V223300:040:071.75 (0.64)F2.07 (0.45)F
26G3V242923:0322:362.48 (0.60)F2.22 (0.35)F
26G6V232921:0021:462.52 (0.83)F1.96 (0.37)F
26H6V1931.523:0423:342.35 (0.42)F
2251V232801:320:572.69 (0.77)F4.82 (1.17)L
21B8V272823:1221:0239.44 (10.58)L16.32 (6.51)L
22A1V242922:0121:554.70 (1.46)L10.74 (3.13)L
22K3V242722:0312.11 (3.68)L4.10 (1.28)L
25Q2V232723:0023:1338.92 (13.79)L28.52 (4.79)L
2622V222722:4121:0614.51 (5.01)L5.20 (3.14)L
26P1V292822:5121:3216.61 (4.68)L5.70 (3.06)L
26R1V2032.521:5619:1013.44 (3.51)L9.67 (2.34)L
Follicular mean ± SEM21.67 ± 0.6229.39 ± 0.6023:19 ± 00:2122:53 ± 00:202.27 ± 0.11n = 92.21 ± 0.07n = 9
Luteal mean ± SEM24.00 ± 1.0028.31 ± 0.6522:54 ± 00:2421:51 ± 00:4119.96 ± 5.16n = 710.63 ± 2.94n = 8
Group mean ± SEM22.76 ± 0.6328.88 ± 0.4523:07 ± 00:1522:26 ± 00:2210.01 ± 3.13n = 166.17 ± 1.70n = 17

Menstrual cycle length is self-reported. Where a range was given, the average length is reported. The means ± SEM of age, menstrual cycle length, bedtime, and dim light melatonin onset (DLMO) for each menstrual phase is based on the menstrual-phase classifications during CR. The daily mean and SD of P4 are shown for each participant on the different study days. Allocation to the luteal phase was based on a cutoff of P4 > 3 ng/mL.

Abbreviations: F, follicular phase; L, luteal phase; M (SD), mean (SD); Menst., menstrual.

Table 1.

Participant Demographic Information

IDAge, yMenst. Length, dBed Time, h:minDLMO, h:minBaselineCR
P4, ng/mL M (SD)Menst. PhaseP4, ng/mL M (SD)Menst. Phase
22JSV202800:3023:272.10 (0.54)F2.13 (0.31)F
22K7V193023:5823:062.06 (0.42)F2.04 (0.26)F
25N6V232822:4020:512.64 (0.77)F2.61 (0.48)F
2614V222800:0123:292.33 (0.44)F2.33 (0.42)F
2692V232823:3322:581.90 (0.33)F2.18 (0.46)F
26F2V223300:040:071.75 (0.64)F2.07 (0.45)F
26G3V242923:0322:362.48 (0.60)F2.22 (0.35)F
26G6V232921:0021:462.52 (0.83)F1.96 (0.37)F
26H6V1931.523:0423:342.35 (0.42)F
2251V232801:320:572.69 (0.77)F4.82 (1.17)L
21B8V272823:1221:0239.44 (10.58)L16.32 (6.51)L
22A1V242922:0121:554.70 (1.46)L10.74 (3.13)L
22K3V242722:0312.11 (3.68)L4.10 (1.28)L
25Q2V232723:0023:1338.92 (13.79)L28.52 (4.79)L
2622V222722:4121:0614.51 (5.01)L5.20 (3.14)L
26P1V292822:5121:3216.61 (4.68)L5.70 (3.06)L
26R1V2032.521:5619:1013.44 (3.51)L9.67 (2.34)L
Follicular mean ± SEM21.67 ± 0.6229.39 ± 0.6023:19 ± 00:2122:53 ± 00:202.27 ± 0.11n = 92.21 ± 0.07n = 9
Luteal mean ± SEM24.00 ± 1.0028.31 ± 0.6522:54 ± 00:2421:51 ± 00:4119.96 ± 5.16n = 710.63 ± 2.94n = 8
Group mean ± SEM22.76 ± 0.6328.88 ± 0.4523:07 ± 00:1522:26 ± 00:2210.01 ± 3.13n = 166.17 ± 1.70n = 17
IDAge, yMenst. Length, dBed Time, h:minDLMO, h:minBaselineCR
P4, ng/mL M (SD)Menst. PhaseP4, ng/mL M (SD)Menst. Phase
22JSV202800:3023:272.10 (0.54)F2.13 (0.31)F
22K7V193023:5823:062.06 (0.42)F2.04 (0.26)F
25N6V232822:4020:512.64 (0.77)F2.61 (0.48)F
2614V222800:0123:292.33 (0.44)F2.33 (0.42)F
2692V232823:3322:581.90 (0.33)F2.18 (0.46)F
26F2V223300:040:071.75 (0.64)F2.07 (0.45)F
26G3V242923:0322:362.48 (0.60)F2.22 (0.35)F
26G6V232921:0021:462.52 (0.83)F1.96 (0.37)F
26H6V1931.523:0423:342.35 (0.42)F
2251V232801:320:572.69 (0.77)F4.82 (1.17)L
21B8V272823:1221:0239.44 (10.58)L16.32 (6.51)L
22A1V242922:0121:554.70 (1.46)L10.74 (3.13)L
22K3V242722:0312.11 (3.68)L4.10 (1.28)L
25Q2V232723:0023:1338.92 (13.79)L28.52 (4.79)L
2622V222722:4121:0614.51 (5.01)L5.20 (3.14)L
26P1V292822:5121:3216.61 (4.68)L5.70 (3.06)L
26R1V2032.521:5619:1013.44 (3.51)L9.67 (2.34)L
Follicular mean ± SEM21.67 ± 0.6229.39 ± 0.6023:19 ± 00:2122:53 ± 00:202.27 ± 0.11n = 92.21 ± 0.07n = 9
Luteal mean ± SEM24.00 ± 1.0028.31 ± 0.6522:54 ± 00:2421:51 ± 00:4119.96 ± 5.16n = 710.63 ± 2.94n = 8
Group mean ± SEM22.76 ± 0.6328.88 ± 0.4523:07 ± 00:1522:26 ± 00:2210.01 ± 3.13n = 166.17 ± 1.70n = 17

Menstrual cycle length is self-reported. Where a range was given, the average length is reported. The means ± SEM of age, menstrual cycle length, bedtime, and dim light melatonin onset (DLMO) for each menstrual phase is based on the menstrual-phase classifications during CR. The daily mean and SD of P4 are shown for each participant on the different study days. Allocation to the luteal phase was based on a cutoff of P4 > 3 ng/mL.

Abbreviations: F, follicular phase; L, luteal phase; M (SD), mean (SD); Menst., menstrual.

For ∼3 weeks before the laboratory study, all participants maintained a self-selected, constant 8-hour sleep/rest/dark schedule, confirmed with calls to a time- and date-stamped voicemail at bedtime and wake time and with actigraphy (Actiwatch-L; Philips Respironics, Bend, OR) and sleep diaries for at least 7 days before entering the laboratory. Participants were asked to refrain from using any prescription and nonprescription medications, supplements, recreational drugs, caffeine, alcohol, or nicotine for the duration of the study. Compliance was confirmed with a urine toxicology test during the 3-week screening phase and upon admission to the laboratory. The study was approved by the Partners Human Research Committee (2007P000566), and written, informed consent was given by participants before starting the study.

Study protocol

Participants were studied individually for 9 days in an environment free of time cues (no access to windows, clocks, live television, radio, or internet and continually supervised by staff trained not to reveal the time) per a standard protocol and as described in Lockley et al. (19) and Gooley et al. (20). Only data from the first 6 days were analyzed and reported herein (Fig. 1); data from a separate analysis of the relationship among reproductive hormones, temperature, and performance are reported elsewhere (21). The first 6 days of the protocol consisted of 3 baseline days with an 8-hour sleep opportunity timed according to the average for 7 days before admission. The baseline study days were followed immediately by a 49-hour, 40-minute CR, followed by an 8-hour recovery sleep (Fig. 1). During the CRs, participants remained awake in a semirecumbent posture in dim light (<3 lux) and were fed hourly isocaloric snacks [150 mEq Na+/100 mEq K+ (±20%); 1.3 × basal energy expenditure; 2000 mL fluids/24 h/d].

Laboratory protocol. White bars represent wake episodes in <190 lux, black bars represent scheduled sleep episodes with lights off (0 lux), gray bars represent wake episodes in dim light (<3 lux) not under CR conditions, and hatched bars represent the CR in dim light (<3 lux).
Figure 1.

Laboratory protocol. White bars represent wake episodes in <190 lux, black bars represent scheduled sleep episodes with lights off (0 lux), gray bars represent wake episodes in dim light (<3 lux) not under CR conditions, and hatched bars represent the CR in dim light (<3 lux).

Lighting

Study lighting conditions have been described in detail previously (19, 20). During baseline days, maximum ambient light (ceiling mounted 4100K fluorescent lamps F96T12/41U/HO/EW, 95 W; F32T8/ADV841/A, 32 W; F25T8/TL841, 25 W; Philips Lighting, The Netherlands), during scheduled wake episodes, was ∼48 µW/cm2 (∼190 lux) when measured vertically and ∼23 µW/cm2 (∼88 lux) when measured horizontally at a height of 187 cm and 137 cm, respectively. Midway through day 3, maximum ambient light was reduced to <0.4 µW/cm2 (<3 lux) when measured vertically and ∼0.6 lux when measured horizontally. This level of light was maintained for the remainder of the study except during scheduled sleep episodes that occurred in darkness.

Hormone measurement

Plasma was collected from an indwelling IV cannula inserted into a forearm vein and kept patent with a heparinized saline infusion (5 IU heparin/mL 0.45% NaCl infused at 40 to 42 mL/h). The IV was inserted on the second baseline day. Blood samples were transferred to EDTA tubes and kept on ice before centrifugation. The plasma fraction was transferred into plastic tubes and stored at −20°C. During the baseline days and the first CR, blood samples were collected every 30 to 60 minutes, and two hourly samples were assayed. Plasma melatonin was assayed using radioimmunoassay (Alpco Diagnostics, Salem NH). Plasma intra- and interassay coefficients of variation were <9% and <11%, respectively, at 1.94 and 16.59 pg/mL. Plasma total E2, FSH, LH, total P4, and SHBG were assayed using the Access Chemiluminescent Immunoassay (Beckman Coulter, Fullerton CA). Intra- and interassay coefficients of variation were 12% to 20% for E2, 3.1% to 5.6% for FSH, 4.3% to 6.4% for LH, 6.11% to 11.19% for P4, and 4.5% to 5.5% for SHBG.

Data analysis

Data are expressed as means ± SEM unless otherwise specified. The first 5 hours of hormone data were excluded from the analysis of CR data to remove any masking effects from the prior sleep episode and changes in posture (22). One participant (26H6V) had no blood collected during baseline days 2 and 3. In one participant (26F2V), LH was below the assay limit of detection, and for another participant (22A1V), FSH was excluded as a result of abnormally high levels (means ± SD, 44.6 ± 7.5 mIU/mL, 4 SDs above the mean), although other reproductive hormone data were retained, as they were within normal ranges.

Cosinor analysis was performed to assess the endogenous circadian regulation of sex hormones under standard sleep/wake baseline and CR conditions. For this analysis, all individual values for each participant were assigned a time relative to their scheduled wake time. For the group analysis, values from each individual were binned in 2-hour intervals, and values from matching bins were averaged across individuals. Hormone data were z-transformed for each individual using the means and SD across the entire data collection interval. Transformed data were fit for the baseline and CR intervals separately using cosinor regression analysis with a single harmonic model that included one 24-hour fundamental component and a linear term: y=A cos {(2πτ)(xϕ)}+(mx+c), where A = amplitude, τ = period (24 hours), ϕ = acrophase, m = slope of linear term, and c = vertical intercept of linear term (16, 23). The linear component was included in the model to estimate the influence of secular changes associated with the CR protocol (e.g., effects of accumulating time awake, repetitive meals, time in a semirecumbent posture, etc.). Melatonin data were z-transformed (24) before fitting with a single harmonic cosinor regression model with one 24-hour fundamental component (16): y=A cos {(2πτ)+(xϕ)}. For both regression analyses, the regressions were considered significant if the amplitude of the fitted rhythm was significantly different from 0, with a one- and two-sided type I error threshold of 0.05 used for individual and group analyses, respectively. Regression analyses were conducted in SAS 9.4. (SAS Inc., Cary, NC). The single harmonic regression analysis of melatonin is included in figures for illustrative purposes. Individual phase estimates (Table 1 and Fig. 2) for melatonin are reported as the time at which melatonin secretion began under dim light (dim light melatonin onset), defined as the time at which the melatonin rhythm crossed the 25% threshold of the three-harmonic, peak-to-trough fitted amplitude (25).

Circadian rhythms in female reproductive hormones in naturally cycling women in the follicular and luteal phases of the menstrual cycle under standard sleep/wake and CR conditions. Group-mean fitted (z-scored data) cosinor functions for women in the (left) follicular and (right) luteal phases of the menstrual cycle under baseline (●) and CR (○) conditions for (A and B) E2, (C and D) FSH, (E and F) LH, (G and H) P4, (I and J) SHBG, and (K and L) melatonin. (B, F, and H) Nonsignificant regression lines are denoted by dashed lines. Corresponding clock times are reported relative to scheduled wake. Time = 0 relative to scheduled wake, defined as 0700, based on the group mean (mean ± SD: 0710 ± 0111) for illustrative purposes. Open, vertical, dashed-lined bars represent group-mean scheduled sleep times.
Figure 2.

Circadian rhythms in female reproductive hormones in naturally cycling women in the follicular and luteal phases of the menstrual cycle under standard sleep/wake and CR conditions. Group-mean fitted (z-scored data) cosinor functions for women in the (left) follicular and (right) luteal phases of the menstrual cycle under baseline (●) and CR (○) conditions for (A and B) E2, (C and D) FSH, (E and F) LH, (G and H) P4, (I and J) SHBG, and (K and L) melatonin. (B, F, and H) Nonsignificant regression lines are denoted by dashed lines. Corresponding clock times are reported relative to scheduled wake. Time = 0 relative to scheduled wake, defined as 0700, based on the group mean (mean ± SD: 0710 ± 0111) for illustrative purposes. Open, vertical, dashed-lined bars represent group-mean scheduled sleep times.

Differences in rhythm acrophases and amplitudes between menstrual phases and between baseline and CR conditions were tested using the extra sum-of-squares F-tests (Prism 7; GraphPad Software, La Jolla CA). Differences in the proportion of participants displaying significant 24-hour rhythms between menstrual phases and baseline and CR conditions were tested using χ2 tests.

Results

Circadian rhythms in female reproductive hormones during the follicular and luteal phases of the menstrual cycle

In group-fitted data, including only women in the follicular phase, there were significant 24-hour rhythms observed in all of the reproductive hormones examined under both baseline and CR conditions (Fig. 2), including E2 (baseline: P < 0.05; CR: P < 0.05), FSH (baseline: P < 0.001; CR: P < 0.01), LH (baseline: P < 0.001; CR: P < 0.05), P4 (baseline: P < 0.001; CR: P < 0.001), and SHBG (baseline: P < 0.01; CR: P < 0.001). In contrast to women in the follicular phase, significant 24-hour rhythms were only observed in group-fitted FSH (baseline: P < 0.05; CR: P < 0.01) and SHBG (baseline: P < 0.01; CR: P < 0.001) in women in the luteal phase under both baseline and CR conditions (Fig. 2). The acrophase and amplitude of the FSH and SHBG rhythms (the only hormones significant in both the follicular and luteal phases) were not different between menstrual phases under either condition (Fig. 2). The unadjusted hormone data are shown in the data repository (26).

At the individual level, the proportions of women in the follicular phase (n = 9) who exhibited a significant 24-hour rhythm under baseline conditions were 44% for E2, 100% for FSH, 75% for LH, 67% for P4, and 44% for SHBG (Table 2). Compared with women in the follicular phase, a similar proportion of women in the luteal phase exhibited significant 24-hour rhythms in E2 (57%), P4 (42%), and SHBG (57%; P > 0.34 for all). In contrast to these hormones, FSH exhibited significant 24-hour rhythms in a substantially smaller proportion of women (33%; χ2 = 8.2, P < 0.01) in the luteal phase, and LH trended toward exhibiting significant 24-hour rhythms in a small proportion of women (29%; χ2 = 3.2, P = 0.07) compared with the follicular phase under baseline conditions (Table 2).

Table 2.

Acrophase and Amplitude of the Individual and Group-Fitted Curves of Female Reproductive Hormones Under Standard Baseline Sleep-Wake Conditions

IDE2FSHLHP4SHBG
Menst. PhasePh, h:minAMPPPh, h:minAMPPPh, h:minAMPPPh, h:minAMPPPh, h:minAMPP
2251VF02:510.360.2514:571.06<0.0113:111.11<0.0111:241.25<0.0122:030.610.06
22JSVF04:160.530.1421:060.700.0418:041.16<0.0114:071.14<0.0101:551.140.17
22K7VF18:530.680.0220:150.98<0.0119:461.03<0.0112:210.850.0113:460.350.14
25N6VF23:570.73<0.0118:161.07<0.0116:180.72<0.0108:010.88<0.0114:520.91<0.01
2614VF03:220.250.4918:190.81<0.0117:460.220.5211:140.500.1814:030.380.15
2692VF02:150.410.1122:520.540.0319:050.78<0.0113:280.440.2216:220.440.18
26F2VF21:410.580.0320:570.600.0512:060.50<0.010.930.92<0.01
26G3VF00:050.790.0220:130.550.1012:390.650.0608:000.79<0.0119:010.590.08
26G6VF03:440.440.2114:490.58<0.0115:020.210.4517:540.530.1011:060.700.10
26H6VF
21B8VL19:140.400.3206:270.110.7706:470.640.1015:070.350.3115:180.68<0.01
22A1VL20:050.75<0.0116:150.62<0.0110:030.780.160.440.210.06
22K3VL06:160.640.0920:170.610.1209:450.530.2013:220.450.1317:010.190.65
25Q2VL03:420.95<0.0108:020.240.3601:530.510.1605:241.01<0.0102:530.440.31
2622VL08:221.00<0.0113:340.580.2016:400.130.7410:321.07<0.0115:400.930.01
26P1VL23:440.090.8013:520.840.0213:150.850.0213:590.480.1414:450.320.40
26R1VL21:060.260.4415:390.83<0.0111:110.530.1218:381.02<0.0116:290.79<0.01
Group mean ± SEMF01:00 ± 01:020.53 ± 0.0619:05 ± 00:550.76 ± 0.0716:29 ± 00:560.74 ± 0.1310:44 ± 00:550.81 ± 0.2817:18 ± 01:320.67 ± 0.09
L00:56 ± 01:580.58 ± 0.1312:58 ± 02:040.53 ± 0.1210:49 ± 01:590.54 ± 0.0810:43 ± 01:590.69 ± 0.3317:37 ± 01:340.51 ± 0.11
Group fitted ± SEMF00:45 ± 01:260.36 ± 0.140.0218:27 ± 00:400.59 ± 0.10<0.0116:20 ± 00:340.56 ± 0.08<0.0110:35 ± 00:450.74 ± 0.16<0.0117:00 ± 01:060.42 ± 0.11<0.01
L01:40 ± 04:000.14 ± 0.150.3516:05 ± 01:470.36 ± 0.160.0412:57 ± 02:220.23 ± 0.150.1313:57 ± 05:020.13 ± 0.170.4517:03 ± 01:220.38 ± 0.13<0.01
IDE2FSHLHP4SHBG
Menst. PhasePh, h:minAMPPPh, h:minAMPPPh, h:minAMPPPh, h:minAMPPPh, h:minAMPP
2251VF02:510.360.2514:571.06<0.0113:111.11<0.0111:241.25<0.0122:030.610.06
22JSVF04:160.530.1421:060.700.0418:041.16<0.0114:071.14<0.0101:551.140.17
22K7VF18:530.680.0220:150.98<0.0119:461.03<0.0112:210.850.0113:460.350.14
25N6VF23:570.73<0.0118:161.07<0.0116:180.72<0.0108:010.88<0.0114:520.91<0.01
2614VF03:220.250.4918:190.81<0.0117:460.220.5211:140.500.1814:030.380.15
2692VF02:150.410.1122:520.540.0319:050.78<0.0113:280.440.2216:220.440.18
26F2VF21:410.580.0320:570.600.0512:060.50<0.010.930.92<0.01
26G3VF00:050.790.0220:130.550.1012:390.650.0608:000.79<0.0119:010.590.08
26G6VF03:440.440.2114:490.58<0.0115:020.210.4517:540.530.1011:060.700.10
26H6VF
21B8VL19:140.400.3206:270.110.7706:470.640.1015:070.350.3115:180.68<0.01
22A1VL20:050.75<0.0116:150.62<0.0110:030.780.160.440.210.06
22K3VL06:160.640.0920:170.610.1209:450.530.2013:220.450.1317:010.190.65
25Q2VL03:420.95<0.0108:020.240.3601:530.510.1605:241.01<0.0102:530.440.31
2622VL08:221.00<0.0113:340.580.2016:400.130.7410:321.07<0.0115:400.930.01
26P1VL23:440.090.8013:520.840.0213:150.850.0213:590.480.1414:450.320.40
26R1VL21:060.260.4415:390.83<0.0111:110.530.1218:381.02<0.0116:290.79<0.01
Group mean ± SEMF01:00 ± 01:020.53 ± 0.0619:05 ± 00:550.76 ± 0.0716:29 ± 00:560.74 ± 0.1310:44 ± 00:550.81 ± 0.2817:18 ± 01:320.67 ± 0.09
L00:56 ± 01:580.58 ± 0.1312:58 ± 02:040.53 ± 0.1210:49 ± 01:590.54 ± 0.0810:43 ± 01:590.69 ± 0.3317:37 ± 01:340.51 ± 0.11
Group fitted ± SEMF00:45 ± 01:260.36 ± 0.140.0218:27 ± 00:400.59 ± 0.10<0.0116:20 ± 00:340.56 ± 0.08<0.0110:35 ± 00:450.74 ± 0.16<0.0117:00 ± 01:060.42 ± 0.11<0.01
L01:40 ± 04:000.14 ± 0.150.3516:05 ± 01:470.36 ± 0.160.0412:57 ± 02:220.23 ± 0.150.1313:57 ± 05:020.13 ± 0.170.4517:03 ± 01:220.38 ± 0.13<0.01

Group-mean values are the average acrophase and amplitude values of each reproductive hormone across the number of participants studied. Individual participant acrophases are reported in clock time. Group-fitted values are the estimated acrophase and amplitude of each reproductive hormone derived from cosinor analysis of group-averaged data. Acrophases of group-fitted data are based on the group-averaged wake time, 0700, derived from the group mean (mean ± SD: 0710 ± 0111).

Abbreviations: AMP, amplitude; Ph, phase.

Table 2.

Acrophase and Amplitude of the Individual and Group-Fitted Curves of Female Reproductive Hormones Under Standard Baseline Sleep-Wake Conditions

IDE2FSHLHP4SHBG
Menst. PhasePh, h:minAMPPPh, h:minAMPPPh, h:minAMPPPh, h:minAMPPPh, h:minAMPP
2251VF02:510.360.2514:571.06<0.0113:111.11<0.0111:241.25<0.0122:030.610.06
22JSVF04:160.530.1421:060.700.0418:041.16<0.0114:071.14<0.0101:551.140.17
22K7VF18:530.680.0220:150.98<0.0119:461.03<0.0112:210.850.0113:460.350.14
25N6VF23:570.73<0.0118:161.07<0.0116:180.72<0.0108:010.88<0.0114:520.91<0.01
2614VF03:220.250.4918:190.81<0.0117:460.220.5211:140.500.1814:030.380.15
2692VF02:150.410.1122:520.540.0319:050.78<0.0113:280.440.2216:220.440.18
26F2VF21:410.580.0320:570.600.0512:060.50<0.010.930.92<0.01
26G3VF00:050.790.0220:130.550.1012:390.650.0608:000.79<0.0119:010.590.08
26G6VF03:440.440.2114:490.58<0.0115:020.210.4517:540.530.1011:060.700.10
26H6VF
21B8VL19:140.400.3206:270.110.7706:470.640.1015:070.350.3115:180.68<0.01
22A1VL20:050.75<0.0116:150.62<0.0110:030.780.160.440.210.06
22K3VL06:160.640.0920:170.610.1209:450.530.2013:220.450.1317:010.190.65
25Q2VL03:420.95<0.0108:020.240.3601:530.510.1605:241.01<0.0102:530.440.31
2622VL08:221.00<0.0113:340.580.2016:400.130.7410:321.07<0.0115:400.930.01
26P1VL23:440.090.8013:520.840.0213:150.850.0213:590.480.1414:450.320.40
26R1VL21:060.260.4415:390.83<0.0111:110.530.1218:381.02<0.0116:290.79<0.01
Group mean ± SEMF01:00 ± 01:020.53 ± 0.0619:05 ± 00:550.76 ± 0.0716:29 ± 00:560.74 ± 0.1310:44 ± 00:550.81 ± 0.2817:18 ± 01:320.67 ± 0.09
L00:56 ± 01:580.58 ± 0.1312:58 ± 02:040.53 ± 0.1210:49 ± 01:590.54 ± 0.0810:43 ± 01:590.69 ± 0.3317:37 ± 01:340.51 ± 0.11
Group fitted ± SEMF00:45 ± 01:260.36 ± 0.140.0218:27 ± 00:400.59 ± 0.10<0.0116:20 ± 00:340.56 ± 0.08<0.0110:35 ± 00:450.74 ± 0.16<0.0117:00 ± 01:060.42 ± 0.11<0.01
L01:40 ± 04:000.14 ± 0.150.3516:05 ± 01:470.36 ± 0.160.0412:57 ± 02:220.23 ± 0.150.1313:57 ± 05:020.13 ± 0.170.4517:03 ± 01:220.38 ± 0.13<0.01
IDE2FSHLHP4SHBG
Menst. PhasePh, h:minAMPPPh, h:minAMPPPh, h:minAMPPPh, h:minAMPPPh, h:minAMPP
2251VF02:510.360.2514:571.06<0.0113:111.11<0.0111:241.25<0.0122:030.610.06
22JSVF04:160.530.1421:060.700.0418:041.16<0.0114:071.14<0.0101:551.140.17
22K7VF18:530.680.0220:150.98<0.0119:461.03<0.0112:210.850.0113:460.350.14
25N6VF23:570.73<0.0118:161.07<0.0116:180.72<0.0108:010.88<0.0114:520.91<0.01
2614VF03:220.250.4918:190.81<0.0117:460.220.5211:140.500.1814:030.380.15
2692VF02:150.410.1122:520.540.0319:050.78<0.0113:280.440.2216:220.440.18
26F2VF21:410.580.0320:570.600.0512:060.50<0.010.930.92<0.01
26G3VF00:050.790.0220:130.550.1012:390.650.0608:000.79<0.0119:010.590.08
26G6VF03:440.440.2114:490.58<0.0115:020.210.4517:540.530.1011:060.700.10
26H6VF
21B8VL19:140.400.3206:270.110.7706:470.640.1015:070.350.3115:180.68<0.01
22A1VL20:050.75<0.0116:150.62<0.0110:030.780.160.440.210.06
22K3VL06:160.640.0920:170.610.1209:450.530.2013:220.450.1317:010.190.65
25Q2VL03:420.95<0.0108:020.240.3601:530.510.1605:241.01<0.0102:530.440.31
2622VL08:221.00<0.0113:340.580.2016:400.130.7410:321.07<0.0115:400.930.01
26P1VL23:440.090.8013:520.840.0213:150.850.0213:590.480.1414:450.320.40
26R1VL21:060.260.4415:390.83<0.0111:110.530.1218:381.02<0.0116:290.79<0.01
Group mean ± SEMF01:00 ± 01:020.53 ± 0.0619:05 ± 00:550.76 ± 0.0716:29 ± 00:560.74 ± 0.1310:44 ± 00:550.81 ± 0.2817:18 ± 01:320.67 ± 0.09
L00:56 ± 01:580.58 ± 0.1312:58 ± 02:040.53 ± 0.1210:49 ± 01:590.54 ± 0.0810:43 ± 01:590.69 ± 0.3317:37 ± 01:340.51 ± 0.11
Group fitted ± SEMF00:45 ± 01:260.36 ± 0.140.0218:27 ± 00:400.59 ± 0.10<0.0116:20 ± 00:340.56 ± 0.08<0.0110:35 ± 00:450.74 ± 0.16<0.0117:00 ± 01:060.42 ± 0.11<0.01
L01:40 ± 04:000.14 ± 0.150.3516:05 ± 01:470.36 ± 0.160.0412:57 ± 02:220.23 ± 0.150.1313:57 ± 05:020.13 ± 0.170.4517:03 ± 01:220.38 ± 0.13<0.01

Group-mean values are the average acrophase and amplitude values of each reproductive hormone across the number of participants studied. Individual participant acrophases are reported in clock time. Group-fitted values are the estimated acrophase and amplitude of each reproductive hormone derived from cosinor analysis of group-averaged data. Acrophases of group-fitted data are based on the group-averaged wake time, 0700, derived from the group mean (mean ± SD: 0710 ± 0111).

Abbreviations: AMP, amplitude; Ph, phase.

Under CR conditions, the proportion of women in the follicular phase who exhibited significant 24-hour rhythms was 89% for E2, 56% for FSH, 63% for LH, 89% for P4, and 78% for SHBG (Table 3). Compared with women in the follicular phase, a similar proportion of women in the luteal phase exhibited significant 24-hour rhythms in E2 (75%), FSH (43%), LH (25%), and SHBG (50%; P > 0.13 for all), but a smaller proportion exhibited significant 24-hour rhythms in P4 (38%; χ2 = 4.9, P < 0.05; Table 3). The phase distribution and time course of significant individual rhythms are shown in the data repository (26).

Table 3.

Acrophase and Amplitude of the Individual and Group-Fitted Curves of Female Reproductive Hormones Under CR Conditions

E2FSHLHP4SHBG
IDMenst. PhasePh, h:minAMPPPh, h:minAMPPPh, h:minAMPPPh, h:minAMPPPh, h:minAMPP
22JSVF03:050.690.0200:350.180.5418:490.64<0.0114:080.370.1916:040.660.05
22K7VF20:080.32<0.0103:390.130.3223:490.78<0.0111:160.650.0406:240.200.32
25N6VF12:300.88<0.0102:310.24<0.0116:060.320.1908:550.78<0.0110:220.72<0.01
2614VF03:300.390.1623:070.250.2705:440.220.4310:480.96<0.0112:570.650.03
2692VF20:280.56<0.0112:290.240.0621:440.480.0210:450.81<0.0109:330.79<0.01
26F2VF20:060.610.0114:360.350.0209:320.79<0.0115:540.72<0.01
26G3VF00:570.710.0516:090.390.0921:330.440.1107:460.540.0316:580.310.05
26G6VF15:520.370.0114:300.260.2417:440.660.04408:540.78<0.0112:220.750.05
26H6VF22:540.760.0216:360.79<0.0113:571.09<0.0111:160.620.0312:540.420.19
21B8VL15:010.300.1705:400.87<0.0123:520.030.9414:330.140.4513:480.79<0.01
2251VL16:030.420.0223:120.070.6318:180.110.6913:590.170.2204:170.160.41
22A1VL23:140.570.1210:340.140.7002:340.200.3908:270.280.21
22K3VL03:350.560.0909:540.470.2400:290.800.0207:270.510.0210:510.270.56
25Q2VL22:300.470.0319:270.420.1119:530.410.1520:220.94<0.0113:441.01<0.01
2622VL18:410.270.0417:001.03<0.0116:330.65<0.0119:270.100.4513:170.76<0.01
26P1VL11:030.210.0623:370.150.0616:120.230.4209:170.230.0206:020.190.29
26R1VL18:070.550.0820:200.250.3219:050.330.3114:500.100.7405:320.720.023
Group mean ±SEMF21:22 ± 01:410.59 ± 0.0619:35 ± 01:560.31 ± 0.0620:26 ± 01:450.58 ± 0.1010:22 ± 00:370.70 ± 0.0612:36 ± 01:080.08 ± 0.04
L19:02 ± 01:580.42 ± 0.0523:53 ± 02:260.47 ± 0.1418:37 ± 01:400.34 ± 0.1012:49 ± 02:170.30 ± 0.1009:30 ± 01:290.19 ± 0.08
Group fitted ±SEMF21:56 ± 01:380.20 ± 0.080.0316:59 ± 01:120.18 ± 0.05<0.0118:17 ± 01:260.29 ± 0.100.0210:11 ± 00:280.64 ± 0.08<0.0112:57 ± 00:550.44 ± 0.10<0.01
L19:26 ± 02:070.12 ± 0.070.0919:08 ± 01:290.23 ± 0.090.0218:50 ± 02:160.15 ± 0.090.1120:28 ± 02:450.11 ± 0.080.1812:05 ± 00:460.43 ± 0.08<0.01
E2FSHLHP4SHBG
IDMenst. PhasePh, h:minAMPPPh, h:minAMPPPh, h:minAMPPPh, h:minAMPPPh, h:minAMPP
22JSVF03:050.690.0200:350.180.5418:490.64<0.0114:080.370.1916:040.660.05
22K7VF20:080.32<0.0103:390.130.3223:490.78<0.0111:160.650.0406:240.200.32
25N6VF12:300.88<0.0102:310.24<0.0116:060.320.1908:550.78<0.0110:220.72<0.01
2614VF03:300.390.1623:070.250.2705:440.220.4310:480.96<0.0112:570.650.03
2692VF20:280.56<0.0112:290.240.0621:440.480.0210:450.81<0.0109:330.79<0.01
26F2VF20:060.610.0114:360.350.0209:320.79<0.0115:540.72<0.01
26G3VF00:570.710.0516:090.390.0921:330.440.1107:460.540.0316:580.310.05
26G6VF15:520.370.0114:300.260.2417:440.660.04408:540.78<0.0112:220.750.05
26H6VF22:540.760.0216:360.79<0.0113:571.09<0.0111:160.620.0312:540.420.19
21B8VL15:010.300.1705:400.87<0.0123:520.030.9414:330.140.4513:480.79<0.01
2251VL16:030.420.0223:120.070.6318:180.110.6913:590.170.2204:170.160.41
22A1VL23:140.570.1210:340.140.7002:340.200.3908:270.280.21
22K3VL03:350.560.0909:540.470.2400:290.800.0207:270.510.0210:510.270.56
25Q2VL22:300.470.0319:270.420.1119:530.410.1520:220.94<0.0113:441.01<0.01
2622VL18:410.270.0417:001.03<0.0116:330.65<0.0119:270.100.4513:170.76<0.01
26P1VL11:030.210.0623:370.150.0616:120.230.4209:170.230.0206:020.190.29
26R1VL18:070.550.0820:200.250.3219:050.330.3114:500.100.7405:320.720.023
Group mean ±SEMF21:22 ± 01:410.59 ± 0.0619:35 ± 01:560.31 ± 0.0620:26 ± 01:450.58 ± 0.1010:22 ± 00:370.70 ± 0.0612:36 ± 01:080.08 ± 0.04
L19:02 ± 01:580.42 ± 0.0523:53 ± 02:260.47 ± 0.1418:37 ± 01:400.34 ± 0.1012:49 ± 02:170.30 ± 0.1009:30 ± 01:290.19 ± 0.08
Group fitted ±SEMF21:56 ± 01:380.20 ± 0.080.0316:59 ± 01:120.18 ± 0.05<0.0118:17 ± 01:260.29 ± 0.100.0210:11 ± 00:280.64 ± 0.08<0.0112:57 ± 00:550.44 ± 0.10<0.01
L19:26 ± 02:070.12 ± 0.070.0919:08 ± 01:290.23 ± 0.090.0218:50 ± 02:160.15 ± 0.090.1120:28 ± 02:450.11 ± 0.080.1812:05 ± 00:460.43 ± 0.08<0.01

Group-mean values are the average acrophase and amplitude values of each reproductive hormone across the number of participants studied. Individual participant acrophases are reported in clock time. Group-fitted values are the estimated acrophase and amplitude of each reproductive hormone derived from cosinor analysis of group-averaged data. Acrophases of group-fitted data are reported in clock time based on the group-averaged wake time, 0700, derived from the group mean (mean ± SD: 0710 ± 0111).

Abbreviations: AMP, amplitude; Ph, phase.

Table 3.

Acrophase and Amplitude of the Individual and Group-Fitted Curves of Female Reproductive Hormones Under CR Conditions

E2FSHLHP4SHBG
IDMenst. PhasePh, h:minAMPPPh, h:minAMPPPh, h:minAMPPPh, h:minAMPPPh, h:minAMPP
22JSVF03:050.690.0200:350.180.5418:490.64<0.0114:080.370.1916:040.660.05
22K7VF20:080.32<0.0103:390.130.3223:490.78<0.0111:160.650.0406:240.200.32
25N6VF12:300.88<0.0102:310.24<0.0116:060.320.1908:550.78<0.0110:220.72<0.01
2614VF03:300.390.1623:070.250.2705:440.220.4310:480.96<0.0112:570.650.03
2692VF20:280.56<0.0112:290.240.0621:440.480.0210:450.81<0.0109:330.79<0.01
26F2VF20:060.610.0114:360.350.0209:320.79<0.0115:540.72<0.01
26G3VF00:570.710.0516:090.390.0921:330.440.1107:460.540.0316:580.310.05
26G6VF15:520.370.0114:300.260.2417:440.660.04408:540.78<0.0112:220.750.05
26H6VF22:540.760.0216:360.79<0.0113:571.09<0.0111:160.620.0312:540.420.19
21B8VL15:010.300.1705:400.87<0.0123:520.030.9414:330.140.4513:480.79<0.01
2251VL16:030.420.0223:120.070.6318:180.110.6913:590.170.2204:170.160.41
22A1VL23:140.570.1210:340.140.7002:340.200.3908:270.280.21
22K3VL03:350.560.0909:540.470.2400:290.800.0207:270.510.0210:510.270.56
25Q2VL22:300.470.0319:270.420.1119:530.410.1520:220.94<0.0113:441.01<0.01
2622VL18:410.270.0417:001.03<0.0116:330.65<0.0119:270.100.4513:170.76<0.01
26P1VL11:030.210.0623:370.150.0616:120.230.4209:170.230.0206:020.190.29
26R1VL18:070.550.0820:200.250.3219:050.330.3114:500.100.7405:320.720.023
Group mean ±SEMF21:22 ± 01:410.59 ± 0.0619:35 ± 01:560.31 ± 0.0620:26 ± 01:450.58 ± 0.1010:22 ± 00:370.70 ± 0.0612:36 ± 01:080.08 ± 0.04
L19:02 ± 01:580.42 ± 0.0523:53 ± 02:260.47 ± 0.1418:37 ± 01:400.34 ± 0.1012:49 ± 02:170.30 ± 0.1009:30 ± 01:290.19 ± 0.08
Group fitted ±SEMF21:56 ± 01:380.20 ± 0.080.0316:59 ± 01:120.18 ± 0.05<0.0118:17 ± 01:260.29 ± 0.100.0210:11 ± 00:280.64 ± 0.08<0.0112:57 ± 00:550.44 ± 0.10<0.01
L19:26 ± 02:070.12 ± 0.070.0919:08 ± 01:290.23 ± 0.090.0218:50 ± 02:160.15 ± 0.090.1120:28 ± 02:450.11 ± 0.080.1812:05 ± 00:460.43 ± 0.08<0.01
E2FSHLHP4SHBG
IDMenst. PhasePh, h:minAMPPPh, h:minAMPPPh, h:minAMPPPh, h:minAMPPPh, h:minAMPP
22JSVF03:050.690.0200:350.180.5418:490.64<0.0114:080.370.1916:040.660.05
22K7VF20:080.32<0.0103:390.130.3223:490.78<0.0111:160.650.0406:240.200.32
25N6VF12:300.88<0.0102:310.24<0.0116:060.320.1908:550.78<0.0110:220.72<0.01
2614VF03:300.390.1623:070.250.2705:440.220.4310:480.96<0.0112:570.650.03
2692VF20:280.56<0.0112:290.240.0621:440.480.0210:450.81<0.0109:330.79<0.01
26F2VF20:060.610.0114:360.350.0209:320.79<0.0115:540.72<0.01
26G3VF00:570.710.0516:090.390.0921:330.440.1107:460.540.0316:580.310.05
26G6VF15:520.370.0114:300.260.2417:440.660.04408:540.78<0.0112:220.750.05
26H6VF22:540.760.0216:360.79<0.0113:571.09<0.0111:160.620.0312:540.420.19
21B8VL15:010.300.1705:400.87<0.0123:520.030.9414:330.140.4513:480.79<0.01
2251VL16:030.420.0223:120.070.6318:180.110.6913:590.170.2204:170.160.41
22A1VL23:140.570.1210:340.140.7002:340.200.3908:270.280.21
22K3VL03:350.560.0909:540.470.2400:290.800.0207:270.510.0210:510.270.56
25Q2VL22:300.470.0319:270.420.1119:530.410.1520:220.94<0.0113:441.01<0.01
2622VL18:410.270.0417:001.03<0.0116:330.65<0.0119:270.100.4513:170.76<0.01
26P1VL11:030.210.0623:370.150.0616:120.230.4209:170.230.0206:020.190.29
26R1VL18:070.550.0820:200.250.3219:050.330.3114:500.100.7405:320.720.023
Group mean ±SEMF21:22 ± 01:410.59 ± 0.0619:35 ± 01:560.31 ± 0.0620:26 ± 01:450.58 ± 0.1010:22 ± 00:370.70 ± 0.0612:36 ± 01:080.08 ± 0.04
L19:02 ± 01:580.42 ± 0.0523:53 ± 02:260.47 ± 0.1418:37 ± 01:400.34 ± 0.1012:49 ± 02:170.30 ± 0.1009:30 ± 01:290.19 ± 0.08
Group fitted ±SEMF21:56 ± 01:380.20 ± 0.080.0316:59 ± 01:120.18 ± 0.05<0.0118:17 ± 01:260.29 ± 0.100.0210:11 ± 00:280.64 ± 0.08<0.0112:57 ± 00:550.44 ± 0.10<0.01
L19:26 ± 02:070.12 ± 0.070.0919:08 ± 01:290.23 ± 0.090.0218:50 ± 02:160.15 ± 0.090.1120:28 ± 02:450.11 ± 0.080.1812:05 ± 00:460.43 ± 0.08<0.01

Group-mean values are the average acrophase and amplitude values of each reproductive hormone across the number of participants studied. Individual participant acrophases are reported in clock time. Group-fitted values are the estimated acrophase and amplitude of each reproductive hormone derived from cosinor analysis of group-averaged data. Acrophases of group-fitted data are reported in clock time based on the group-averaged wake time, 0700, derived from the group mean (mean ± SD: 0710 ± 0111).

Abbreviations: AMP, amplitude; Ph, phase.

Comparison of circadian rhythms in female reproductive hormones under standard sleep-wake and CR conditions

For women in the follicular phase, the amplitude of the group-fitted FSH rhythm was significantly smaller (P < 0.001; Fig 2C) under CR compared with standard sleep-wake conditions (baseline), and a similar trend was observed for LH (P = 0.051; Fig. 2E). The acrophase of the SHBG rhythm was significantly earlier (∼4 hours) under CR conditions compared with baseline (P < 0.05; Fig. 2I). There were no differences in the phase or amplitude of the other hormones between baseline and CR in women in the follicular phase. At the individual level, compared with baseline, a larger proportion of women in the follicular phase exhibited 24-hour rhythms in E2 (baseline = 44%, CR = 89%; χ2 = 4.0, P < 0.05) under CR conditions, and a smaller proportion exhibited 24-hour rhythms in FSH (baseline = 100%, CR = 56%; χ2 = 5.1, P < 0.05). The proportion of participants in the follicular phase exhibiting significant rhythms at the individual level was similar for baseline and CR for LH, P4, and SHBG (P > 0.15 for all).

For women in the luteal phase, the acrophase of the SHBG rhythm was significantly earlier (∼5 hours) under CR compared with baseline (P < 0.01; Fig. 2J). At the individual level, the number of women in the luteal phase who exhibited 24-hour rhythms during CR was similar to under baseline conditions for all hormones (P > 0.46 for all).

Discussion

Our findings demonstrate that key endocrine regulators of female reproductive function, including sex-steroid hormones, gonadotropins, and SHBG, are under endogenous circadian regulation. Moreover, this regulation appears to depend on menstrual phase, such that most variables were robustly rhythmic during the follicular phase but not the luteal phase of the menstrual cycle. These results have important implications for understanding the regulation of reproductive function.

Our results show that the circadian clock differentially regulates female reproductive hormones in humans, depending on the menstrual phase. Although the data from our study preclude directly testing neuroendocrine mechanisms that can explain this differential circadian regulation of female reproductive hormones across the menstrual cycle, there are a number of possibilities. The loss of rhythmicity in P4 may be a result of a shift in the primary source of P4 from the adrenal cortex during the follicular phase to the corpus luteum during the luteal phase, which may not be rhythmic. Although the expression of steroidogenic acute regulatory protein, which catalyzes the rate-limiting step of steroidogenesis, is regulated by the circadian core clock gene brain and muscle aryl hydrocarbon receptor nuclear translocator-like l (BMAL1), P4 secretion by the corpus luteum is arrhythmic in both bmal1−/− and heterozygote controls (27). Additionally, the arrhythmic and elevated P4 levels during the luteal phase may explain the loss in LH rhythmicity during the luteal phase, given that P4 suppresses the pulsatile secretion of LH and mean circulating levels of LH (28). The significant circadian regulation of sex steroids, gonadotropins, and SHBG during the follicular phase, preceding ovulation but not after ovulation, suggests that complex reproductive outcomes, such as ovulation, may be temporally coordinated by the circadian system.

Although our study did not directly assess the effect of circadian disruption of these rhythms, several epidemiologic studies have reported the adverse effects of shift work, which disrupts circadian rhythms, on female reproductive health. Data from the 1993 cohort of the Women’s Health Study that followed 71,077 women showed a 23% increased risk of irregular and abnormally long (>40 days) or short (<21 days) menstrual cycles with increasing months of rotating shift work (29). Other studies have reported that female shift workers self-report changes to cycle length, menstrual flow, menstrual pain, and duration of menses (30) and also report more problems conceiving [i.e., increased time to pregnancy (31)] and maintaining pregnancy [i.e., increased odds of miscarriage (32)]. Furthermore, two prospective studies have shown a higher rate of monophasic (i.e., no basal body temperature increase; anovulatory cycles) and irregular ovulatory cycles in rotating shift workers compared with day workers (33, 34). We hypothesize, therefore, that these circadian rhythms in reproductive hormones are disrupted during shift work and that their disruption contributes to the impaired reproductive function observed in female shift workers. Future studies should examine this hypothesis explicitly. Diurnal changes in gonadotropins and sex steroids have been reported previously [for example, see Veldhuis et al. (3), Rossmanith and Lauritzen (4), Mortola et al. (5), and Bao et al. (6)], although the rhythmic oscillations have been examined under constant environmental and behavioral conditions only once previously (15). We examined both ambulatory and constant conditions in the same study to compare directly evoked responses and the endogenous circadian regulation of these hormones. Our results corroborate previous reports of significant 24-hour oscillations under ambulatory conditions. In contrast, our results are not consistent with the previous study in which a CR was used, wherein no significant circadian oscillations in gonadotropins and E2 were detected in healthy young women in the early follicular phase of the menstrual cycle (15). Subtle differences in study protocols and possible differences in menstrual-phase estimation may have contributed to these different conclusions. For example, the ambient light levels were dimmer in our protocol (<3 lux) than in the previous protocol (50 to 150 lux). Whereas the effect of ambient light exposure on reproductive hormones in women is not well characterized, melatonin has been shown to modulate reproductive hormone levels (35), and melatonin is suppressed by light exposure. Therefore, it is possible that the higher intensity ambient light used in the prior study may have altered melatonin levels, which indirectly then modulated the levels of the reproductive hormones, masking the underlying circadian regulation. Future studies are necessary to characterize systematically the influence of environmental factors (e.g., light) on these reproductive rhythms.

Additionally, in our study, menstrual phase was determined using a threshold for circulating P4, which does not allow further categorization into the early or late follicular phase. A trend toward an increase in amplitude of the diurnal rhythm in free E2 between the early to late follicular phase has been previously reported, although the difference did not reach statistical significance in that study with a limited sample size (6). The study did, however, observe a significant difference in the acrophase of the rhythm between the early and late follicular phase. Additionally, the study by Bao et al. (6) did not characterize the rhythms under constant conditions, which may have masked the subtle differences in rhythm characteristics between the early and late follicular phase. It is possible that there was a wider distribution in the follicular phase in our study compared with the study by Klingman et al. (15), with more individuals in our study in the late follicular phase than in the early follicular phase. Therefore, there may be differences in circadian rhythms in reproductive hormones, not only between the follicular and luteal phases but also perhaps within a menstrual phase (e.g., early vs late follicular phase) that warrants further examination. Furthermore, we cannot rule out a possible order effect or the time-into-protocol effect that may have contributed to the differential hormone levels between CR and baseline conditions.

Another pertinent finding in the study was the remarkable variability in the proportion of women who were rhythmic for a given endocrine variable. The proportion ranged between 100% for FSH in the follicular phase under baseline conditions to 29% for LH in the luteal phase under baseline conditions. We and others have previously reported similarly large variability in the number of significant rhythms observed at the individual level for proinflammatory cytokines and chemokines (16), lipids (24), and metabolites (36). It is possible that the cosinor-based approach of the determination of rhythms is not adequately robust for individual-level data, which may not fit a sinusoid profile as well as group-averaged data. Additionally, this variability may also represent varying levels of susceptibility of each individual to rhythm perturbations caused by external or internal factors, which may, in turn, explain the differences in outcomes among individuals in response to external influences, such as shift work.

This study has several limitations. Whereas we collected blood samples every 2 hours longitudinally over ∼96 hours, most of these hormones are known to be pulsatile. Whereas a less frequent sampling over a long duration of time may be sufficient when trying to detect gross 24-hour rhythmicity, future studies with more a frequent sampling are needed to examine whether amplitude and frequency of gonadotropin pulsatility also exhibit circadian rhythms. Additionally, future studies with higher sensitivity assays (e.g., HPLC) are warranted to confirm the precise amplitude of the rhythms. The limited sample size and the age range in the current study also limit the generalizability of the results. Additionally, results from individual-based analyses should be interpreted with caution because of the limited sample size and limited statistical power. Additional studies with a greater sample size and statistical power are warranted, however, to test the interesting individual-level variability observed in the current study.

Overall, our study demonstrates that many female reproductive hormones are under endogenous circadian regulation. Our findings provide a potential mechanism by which shift work adversely affects reproductive function through disruption of the circadian regulation of reproductive hormones.

Acknowledgments

We thank the technical, dietary, and laboratory staff, nurses, and physicians; participant recruiters; and the study participants at the Center for Clinical Investigation and Division of Sleep and Circadian Disorders, Brigham and Women’s Hospital.

Financial Support: This work was supported by the National Institute of Mental Health (2R01 MH45130-11A1; to C.A.C. and S.W.L.), National Center for Complementary and Alternative Medicine (R01 AT002129; to C.A.C. and S.W.L.), and National Institute of Environmental Health Sciences (R21 ES017112-01A1; to S.W.L.). C.A.C. and S.W.L. were supported, in part, by the National Space Biomedical Research Institute through NASA NCC 9-58. The project was supported by Brigham and Women’s Hospital General Clinical Research Center Grant M01-RR02635.

Additional Information

Disclosure Summary: L.K.G. and J.J.G. have nothing to disclose. S.A.R. owns equity in Melcort Inc.; has provided paid consulting services to Sultan & Knight Limited and Bambu Vault LLC; and has received honoraria as an invited speaker and travel funds from Starry Skies Lake Superior, University of Minnesota Medical School, PennWell Corp., and Seoul Semiconductor Co. Ltd. S.M.W.R. is a program leader and serves as a consultant to the CRC for Alertness, Safety and Productivity, Australia; reports receiving research grants from the CRC for Alertness, Safety and Productivity, Philips Respironics, Rio Tinto, Shell, Linfox Australia, and Teva Pharma Australia; and has received equipment support and consultancy fees through his institution from Vanda Pharmaceuticals, Optalert, Tyco Healthcare, Compumedics, BHP, and Teva Pharmaceuticals, which are not related to this paper. C.A.C. reports grants from Cephalon Inc.; Jazz Pharmaceuticals Plc., Inc.; National Football League Charities; Optum; Philips Respironics, Inc.; Regeneron Pharmaceuticals; ResMed Foundation; San Francisco Bar Pilots; Sanofi S.A.; Sanofi-Aventis, Inc.; Schneider Inc.; Sepracor, Inc.; Mary Ann & Stanley Snider via Combined Jewish Philanthropies; Sysco; Takeda Pharmaceuticals; Teva Pharmaceuticals Industries, Ltd.; and Wake Up Narcolepsy and personal fees from Bose Corporation; Boston Celtics; Boston Red Sox; Cephalon, Inc.; Columbia River Bar Pilots; Ganésco Inc.; Institute of Digital Media and Child Development; Klarman Family Foundation; Samsung Electronics; Quest Diagnostics, Inc.; Teva Pharma Australia; Vanda Pharmaceuticals; Washington State Board of Pilotage Commissioners; and Zurich Insurance Company, Ltd. In addition, C.A.C. holds a number of process patents in the field of sleep/circadian rhythms (e.g., photic resetting of the human circadian pacemaker) and holds an equity interest in Vanda Pharmaceuticals, Inc. Since 1985, C.A.C. has also served as an expert on various legal and technical cases related to sleep and/or circadian rhythms, including those involving the following commercial entities: Casper Sleep Inc., Comair/Delta Airlines, Complete General Construction Company, Federal Express, Greyhound, HG Energy LLC, Purdue Pharma L.P., South Carolina Central Railroad Co., Steel Warehouse Inc., Stric-Lan Companies LLC, Texas Premier Resource LLC, and United Parcel Service. C.A.C. receives royalties from the New England Journal of Medicine; McGraw Hill; Houghton Mifflin Harcourt/Penguin; and Philips Respironics, Inc., for the Actiwatch 2 and Actiwatch Spectrum devices. C.A.C.’s interests were reviewed and managed by Brigham and Women’s Hospital and Partners HealthCare in accordance with their conflict of interest policies. S.W.L. has had a number of commercial interests in the last 24 months (2017 to 2019). No other interests are directly related to the research or topic reported in this paper but in the interests of full disclosure, are mentioned. S.W.L. has received consulting fees from BHP Billiton, EyeJust Inc., Noble Insights, and Team C Racing; honoraria and/or paid travel from BHP Billiton, DIN, IES, Ineos, SLTBR, and Teague; has current consulting contracts with Akili Interactive, Apex 2100 Ltd., Consumer Sleep Solutions, Headwaters Inc., Hintsa Performance AG, Light Cognitive, Lighting Science Group Corporation, Mental Workout, PlanLED, Six Senses, Stantec, and Wyle Integrated Science and Engineering; has received unrestricted equipment gifts from Bionetics Corporation and F. Lux Software LLC and royalties from Oxford University Press; and has served as a paid expert in legal proceedings related to light, sleep, and health.

Data Availability: Restrictions apply to the availability of data generated or analyzed during this study to preserve patient confidentiality or because they were used under license. The corresponding author will on request detail the restrictions and any conditions under which access to some data may be provided.

Abbreviations:

    Abbreviations:
     
  • CR

    constant routine

  •  
  • E2

    estradiol

  •  
  • P4

    progesterone

References and Notes

1.

Boden
MJ
,
Kennaway
DJ
.
Circadian rhythms and reproduction
.
Reproduction
.
2006
;
132
(
3
):
379
392
.

2.

Stocker
LJ
,
Macklon
NS
,
Cheong
YC
,
Bewley
SJ
.
Influence of shift work on early reproductive outcomes: a systematic review and meta-analysis
.
Obstet Gynecol
.
2014
;
124
(
1
):
99
110
.

3.

Veldhuis
JD
,
Christiansen
E
,
Evans
WS
,
Kolp
LA
,
Rogol
AD
,
Johnson
ML
.
Physiological profiles of episodic progesterone release during the midluteal phase of the human menstrual cycle: analysis of circadian and ultradian rhythms, discrete pulse properties, and correlations with simultaneous luteinizing hormone release
.
J Clin Endocrinol Metab
.
1988
;
66
(
2
):
414
421
.

4.

Rossmanith
WG
,
Lauritzen
C
.
The luteinizing hormone pulsatile secretion: diurnal excursions in normally cycling and postmenopausal women
.
Gynecol Endocrinol
.
1991
;
5
(
4
):
249
265
.

5.

Mortola
JF
,
Laughlin
GA
,
Yen
SS
.
A circadian rhythm of serum follicle-stimulating hormone in women
.
J Clin Endocrinol Metab
.
1992
;
75
(
3
):
861
864
.

6.

Bao
AM
,
Liu
RY
,
van Someren
EJ
,
Hofman
MA
,
Cao
YX
,
Zhou
JN
.
Diurnal rhythm of free estradiol during the menstrual cycle
.
Eur J Endocrinol
.
2003
;
148
(
2
):
227
232
.

7.

Filicori
M
,
Santoro
N
,
Merriam
GR
,
Crowley
WF
Jr
.
Characterization of the physiological pattern of episodic gonadotropin secretion throughout the human menstrual cycle
.
J Clin Endocrinol Metab
.
1986
;
62
(
6
):
1136
1144
.

8.

Hall
JE
,
Schoenfeld
DA
,
Martin
KA
,
Crowley
WF
Jr
.
Hypothalamic gonadotropin-releasing hormone secretion and follicle-stimulating hormone dynamics during the luteal-follicular transition
.
J Clin Endocrinol Metab
.
1992
;
74
(
3
):
600
607
.

9.

Fujimoto
VY
,
Clifton
DK
,
Cohen
NL
,
Soules
MR
.
Variability of serum prolactin and progesterone levels in normal women: the relevance of single hormone measurements in the clinical setting
.
Obstet Gynecol
.
1990
;
76
(
1
):
71
78
.

10.

Venturoli
S
,
Porcu
E
,
Fabbri
R
,
Magrini
O
,
Gammi
L
,
Paradisi
R
,
Forcacci
M
,
Bolzani
R
,
Flamigni
C
.
Episodic pulsatile secretion of FSH, LH, prolactin, oestradiol, oestrone, and LH circadian variations in polycystic ovary syndrome
.
Clin Endocrinol (Oxf)
.
1988
;
28
(
1
):
93
107
.

11.

Soules
MR
,
Steiner
RA
,
Cohen
NL
,
Bremner
WJ
,
Clifton
DK
.
Nocturnal slowing of pulsatile luteinizing hormone secretion in women during the follicular phase of the menstrual cycle
.
J Clin Endocrinol Metab
.
1985
;
61
(
1
):
43
49
.

12.

Soules
MR
,
Clifton
DK
,
Steiner
RA
,
Cohen
NL
,
Bremner
WJ
.
The corpus luteum: determinants of progesterone secretion in the normal menstrual cycle
.
Obstet Gynecol
.
1988
;
71
(
5
):
659
666
.

13.

Kivelä
A
,
Kauppila
A
,
Ylöstalo
P
,
Vakkuri
O
,
Leppäluoto
J
.
Seasonal, menstrual and circadian secretions of melatonin, gonadotropins and prolactin in women
.
Acta Physiol Scand
.
1988
;
132
(
3
):
321
327
.

14.

Czeisler
CA
,
Klerman
EB
.
Circadian and sleep-dependent regulation of hormone release in humans
.
Recent Prog Horm Res
.
1999
;
54
:
97
130
, discussion 130–132.

15.

Klingman
KM
,
Marsh
EE
,
Klerman
EB
,
Anderson
EJ
,
Hall
JE
.
Absence of circadian rhythms of gonadotropin secretion in women
.
J Clin Endocrinol Metab
.
2011
;
96
(
5
):
1456
1461
.

16.

Rahman
SA
,
Castanon-Cervantes
O
,
Scheer
FA
,
Shea
SA
,
Czeisler
CA
,
Davidson
AJ
,
Lockley
SW
.
Endogenous circadian regulation of pro-inflammatory cytokines and chemokines in the presence of bacterial lipopolysaccharide in humans
.
Brain Behav Immun
.
2015
;
47
:
4
13
.

17.

Israel
R
,
Mishell
DR
Jr
,
Stone
SC
,
Thorneycroft
IH
,
Moyer
DL
.
Single luteal phase serum progesterone assay as an indicator of ovulation
.
Am J Obstet Gynecol
.
1972
;
112
(
8
):
1043
1046
.

18.

Stricker
R
,
Eberhart
R
,
Chevailler
MC
,
Quinn
FA
,
Bischof
P
,
Stricker
R
.
Establishment of detailed reference values for luteinizing hormone, follicle stimulating hormone, estradiol, and progesterone during different phases of the menstrual cycle on the Abbott ARCHITECT analyzer
.
Clin Chem Lab Med
.
2006
;
44
(
7
):
883
887
.

19.

Lockley
SW
,
Brainard
GC
,
Czeisler
CA
.
High sensitivity of the human circadian melatonin rhythm to resetting by short wavelength light
.
J Clin Endocrinol Metab
.
2003
;
88
(
9
):
4502
4505
.

20.

Gooley
JJ
,
Rajaratnam
SMW
,
Brainard
GC
,
Kronauer
RE
,
Czeisler
CA
,
Lockley
SW
.
Spectral responses of the human circadian system depend on the irradiance and duration of exposure to light
.
Sci Transl Med
.
2010
;
2
(
31
):
31ra33
.

21.

Grant
LK
,
Gooley
JJ
,
St Hilaire
MA
,
Rajaratnam
SMW
,
Brainard
GC
,
Czeisler
CA
,
Lockley
SW
,
Rahman
SA
.
Menstrual phase-dependent differences in neurobehavioral performance: the role of temperature and the progesterone/estradiol ratio
.
Sleep
2019
(in press).

22.

Brown
EN
,
Czeisler
CA
.
The statistical analysis of circadian phase and amplitude in constant-routine core-temperature data
.
J Biol Rhythms
.
1992
;
7
(
3
):
177
202
.

23.

Chua
ECP
,
Shui
G
,
Cazenave-Gassiot
A
,
Wenk
MR
,
Gooley
JJ
.
Changes in plasma lipids during exposure to total sleep deprivation
.
Sleep
.
2015
;
38
(
11
):
1683
1691
.

24.

Chua
ECP
,
Shui
G
,
Lee
ITG
,
Lau
P
,
Tan
LC
,
Yeo
SC
,
Lam
BD
,
Bulchand
S
,
Summers
SA
,
Puvanendran
K
,
Rozen
SG
,
Wenk
MR
,
Gooley
JJ
.
Extensive diversity in circadian regulation of plasma lipids and evidence for different circadian metabolic phenotypes in humans
.
Proc Natl Acad Sci USA
.
2013
;
110
(
35
):
14468
14473
.

25.

Klerman
EB
,
Gershengorn
HB
,
Duffy
JF
,
Kronauer
RE
.
Comparisons of the variability of three markers of the human circadian pacemaker
.
J Biol Rhythms
.
2002
;
17
(
2
):
181
193
.

26.

Rahman
SA
,
Grant
LK
,
Gooley
JJ
,
Rajaratnam
SMW
,
Czeisler
CA
,
Lockley
SW
.
Data from: Endogenous circadian regulation of female reproductive hormones. figshare 2019. Deposited 31 May 2019. https://dx.doi.org/10.6084/m9.figshare.8214779.v1
.

27.

Ratajczak
CK
,
Boehle
KL
,
Muglia
LJ
.
Impaired steroidogenesis and implantation failure in Bmal1-/- mice
.
Endocrinology
.
2009
;
150
(
4
):
1879
1885
.

28.

Soules
MR
,
Steiner
RA
,
Clifton
DK
,
Cohen
NL
,
Aksel
S
,
Bremner
WJ
.
Progesterone modulation of pulsatile luteinizing hormone secretion in normal women
.
J Clin Endocrinol Metab
.
1984
;
58
(
2
):
378
383
.

29.

Lawson
CC
,
Whelan
EA
,
Lividoti Hibert
EN
,
Spiegelman
D
,
Schernhammer
ES
,
Rich-Edwards
JW
.
Rotating shift work and menstrual cycle characteristics
.
Epidemiology
.
2011
;
22
(
3
):
305
312
.

30.

Labyak
S
,
Lava
S
,
Turek
F
,
Zee
P
.
Effects of shiftwork on sleep and menstrual function in nurses
.
Health Care Women Int
.
2002
;
23
(
6-7
):
703
714
.

31.

Ahlborg
G
Jr,
Axelsson
G
,
Bodin
L
.
Shift work, nitrous oxide exposure and subfertility among Swedish midwives
.
Int J Epidemiol
.
1996
;
25
(
4
):
783
790
.

32.

Axelsson
G
,
Ahlborg
G
Jr
,
Bodin
L
.
Shift work, nitrous oxide exposure, and spontaneous abortion among Swedish midwives
.
Occup Environ Med
.
1996
;
53
(
6
):
374
378
.

33.

Hatch
MC
,
Figa-Talamanca
I
,
Salerno
S
.
Work stress and menstrual patterns among American and Italian nurses
.
Scand J Work Environ Health
.
1999
;
25
(
2
):
144
150
.

34.

Wan
GH
,
Chung
FF
.
Working conditions associated with ovarian cycle in a medical center nurses: a Taiwan study
.
Jpn J Nurs Sci
.
2012
;
9
(
1
):
112
118
.

35.

Cagnacci
A
,
Paoletti
AM
,
Soldani
R
,
Orrù
M
,
Maschio
E
,
Melis
GB
.
Melatonin enhances the luteinizing hormone and follicle-stimulating hormone responses to gonadotropin-releasing hormone in the follicular, but not in the luteal, menstrual phase
.
J Clin Endocrinol Metab
.
1995
;
80
(
4
):
1095
1099
.

36.

Grant
LK
,
Ftouni
S
,
Nijagal
B
,
De Souza
DP
,
Tull
D
,
McConville
MJ
,
Rajaratnam
SMW
,
Lockley
SW
,
Anderson
C
.
Circadian and wake-dependent changes in human plasma polar metabolites during prolonged wakefulness: a preliminary analysis
.
Sci Rep
.
2019
;
9
(
1
):
4428
.