Original research

Patterns of endogenous and exogenous ovarian hormone modulation on recovery metrics across the menstrual cycle

Abstract

Introduction As the number of female athletes competing rises globally, training methodologies should reflect sex differences across critical metrics of adaptation to training. Surrogate markers of the autonomic nervous system (ANS) used for monitoring training load are heart rate variability (HRV) and resting heart rate (RHR). The aim was to investigate ovarian hormone effects on standard recovery metrics (HRV, RHR, respiratory rate (RR) and sleep duration) across a large population of female athletes.

Methods A retrospective study analysed 362 852 days of data representing 13 535 menstrual cycles (MC) from 4594 respondents (natural MC n=3870, BC n=455, progestin-only n=269) for relationships and/or differences between endogenous and exogenous ovarian hormones on ANS.

Results HRV and return to baseline (recovery) decreased as resting HR and RR increased (p<0.001) from the early follicular to the late luteal phase of the MC. Patterning was paradoxical across phases for users of combined hormonal contraception (BC) as compared with the patterning of the MC. HRV and recovery start elevated and drop off quickly during the withdrawal bleed, rising through the active pill weeks (p<0.001). Progestin-only users had similar patterning as the MC. The relationship between normalised recovery and previous day strain is modulated by birth control type. BC exhibited steeper declines in recovery with additional strain-normalised recovery decreases by an additional 0.0055±0.00135 (p<0.001) per unit of strain; with no significant difference between MC and progestin-only (p=0.19).

Conclusion The patterning of ANS modulation from ovarian hormones is significantly different between naturally cycling women and those on BC, with the patterning dependent on the type of contraception used.

Key messages

What are the new findings

  • Ovarian hormones modulate standard recovery metrics (respiratory rate, sleep measures) used to determine adaptation to training load and stress across the menstrual cycle.

  • Indices of lower recovery and higher strain are associated with the luteal phase of the natural menstrual cycle due to progesterone stimulatory influence on the autonomic nervous system (ANS).

  • Exogenous hormones of birth control (BC) modulate the recovery metrics in a different pattern than the natural menstrual cycle, with reduced indices of adaptation to stress across all pill phases.

  • Exogenous hormone use significantly impacts day to day recovery to increased cardiovascular strain as compared with the natural menstrual cycle.

  • Progestin-only BC exhibit a unique modulation of recovery indices, with a lessened effect on the ANS in the late luteal phase.

  • Progestin-only BC exhibit similar patterns to the natural menstrual cycle with regard to cardiovascular strain and next day recovery.

How might it impact on clinical practice in the future

  • Current practices in sport and exercise training programmes do not account for the influence of ovarian hormones on the ANS, nor are the differences between endogenous and exogenous hormone profiles acknowledged in coaching and exercise prescription. By understanding the modulatory effect exerted by the different hormone profiles, improved training methodologies may be employed to improve the adaptations to training loads of female athletes.

Introduction

The female hormonal cycle and its systemic effects on the body is a crucial area for research as it has been found that women frequently experience different adaptations and stress responses to their male counterparts, at times as a result of female hormonal fluctuations.1 The menstrual cycle (MC) is characterised by fluctuations in several hormones, most notably the gonadal steroids, oestrogen and progesterone, and is divided into four main hormonal phases, the early follicular, mid-follicular, early luteal and late-luteal phases. Throughout each phase, fluctuations in hormones trigger not only changes in the reproductive system, but also in all the tissues of the body,1 modulating sleep,2 3 circadian rhythms4 and the autonomic nervous system (ANS).5–7 For instance, there is evidence of a decrease of the vagal dominance on the heart from the follicular to the luteal phase.6–8 Moreover, baroreflex sensitivity (BRS) is higher in the luteal phase of the cycle, presumably due to the elevation of progesterone and its stimulation of the phrenic nerve and overall sympathetic drive.6 8

A surrogate marker for ANS activity is heart rate variability (HRV), the recording of in-between heartbeat variations that are modulated by parasympathetic and sympathetic inputs.9–11 HRV is sensitive to homoeostatic perturbations like fatigue, physiological and psychological stress.8 9 Across sports, HRV is widely used as an appropriate criterion for training load quantification,9 as the combined intensity and duration of a training session impacts the magnitude and time required for recovery from the induced stress response, shifting the autonomic balance. Periods of higher cumulative stress are typically highlighted by a reduction in HRV and in situations of positive adaptations to training, HRV should increase or remain stable, with increasing training loads. Thus, resting heart rate (RHR) and HRV have been used to adjust training loads and identify positive and negative adaptations to aerobic training.10 11

As an important marker for quantifying adaptation to training load, there is a lack of attention to how female sex hormones modulate ANS activity within exercise training/prescription, coaching and applied sport science,12 13 impacting quantification and monitoring training loads specific to female athletes. In addition, regardless of sport, there is a high percentage of female athletes who use hormonal contraception .14 Combined oral contraceptives (COC) contain derivatives of oestrogens and progestogens, which downregulate the regular cyclic hormone activity of the hypothalamic–pituitary–gonadal axis, resulting in physiological responses which differ to those of endogenous ovarian hormones.15 For example, Wilczak et al16 compared indices of the ANS between women using COC and naturally cycling women, finding a downregulation of BRS with COC use. Moreover, Teixeira et al17 determined COC do not affect HRV across the active pill phases of use.

The aim of this project was to investigate ovarian hormone effects on standard recovery metrics of HRV, RHR, respiratory rate (RR) and sleep duration across a large population of female athletes, and to determine any differences between endogenous and exogenous hormones on these metrics.

Methods

Data were collected retrospectively from subscribers to the WHOOP, Inc. platform. Specifically, RHR, HRV, RR and the percentage of time in the four sleep stages (wake, light sleep, REM and slow wave sleep) were captured using a wrist‐worn multisensor (triaxial accelerometer, optical sensor, capacitive touch sensor and ambient temperature sensor) device (WHOOP, Inc., Boston, Massachusetts USA). Participant characteristics, including age, weight, height and athletic status (elite or recreational), were determined through mandatory in-app questions that must be completed during the sign-up process. Data for the current study were analysed from an optional survey of women who recorded menstruation on the app from 29 October 2020 to 9 December 2020. Inclusion criteria included regular, natural MCs lasting between 25 and 35 days, or the use of hormonal contraceptives (including the inactive pill week). Exclusion criteria included: irregular cycles (as defined as missing two or more cycles in a 6-month period or shorter than 25 days or longer than 35 days), uncertainty of the type of birth control (BC) used, and/or having stopped using BC within the past 9 months. All participants in the study had consented at registration to the use of their data for the purposes of scientific research and could withdraw consent at any time.

Measures of RHR and HRV were evaluated during the last slow wave cycle of sleep, which was shown to be the best recording condition because of a high level of environmental and respiratory standardisation.18 19 RR was calculated as the median RR over the entire sleep period. RHR, HRV and RR are derived via photoplethysmography recordings at the wrist, which has been demonstrated to be a valid method for measurement of R‐R intervals via ECG while at rest.20–22 Recovery is a proprietary metric of the WHOOP, Inc. platform that combines nighttime measures of heart rate, HRV, RR and sleep duration into a once-daily metric used to guide recommended training load.23 Recovery ranges from 0% to 100% with higher values indicating greater readiness to adapt to a training load. All measures have been previously validated in female athletes as a means to quantify training load and training stress.12 23–26

The WHOOP data are the property of WHOOP, Inc. The authors (LW and ERC) are employees of WHOOP, Inc and are permitted to access the WHOOP data server.

Statistical analysis

This study’s primary interest was the effect of menstrual phase on physiological measures of recovery and the modulation of endogenous and exogenous hormones on these effects. Daily measures of HRV, resting HR, RR and recovery were analysed for the study population. The menstrual phase for each individual and day was estimated by taking the corresponding MC and dividing it into four equally spaced intervals with the first day of self-reported menstruation defining the beginning of the early follicular phase. Days that fell in the first 25% of the cycle were designated as early follicular phase. The next 25% interval was designated as late follicular phase. The next 25% was designated as early luteal phase, and the last 25% was designated as late luteal phase. Days corresponding to the withdrawal bleed of the BC were designated as the early follicular phase, with the subsequent active pill weeks designated as ‘week 2, week 3, week 4’ to correspond with the natural cycle ‘late follicular, early luteal, late luteal’ phases. All statistical analyses were performed in R (V.4.0.3).

Statistical analysis: MC phase and HC

In the analysis, linear models were used to relate each physiological measure (normalised HRV, normalised RR, normalised resting HR, normalised recovery score) to menstrual phase and BC type. Interaction effects were included in the models. Extra sum-of-squares F-tests were performed to test for the significance of the main effects and interaction effects. For hypothesis tests with multiple comparisons, a Bonferroni correction was applied to the alpha level of 0.05 to generate familywise alpha levels.

Statistical analysis: cardiovascular strain and recovery

Previous day cardiovascular strain was added to the existing model, including interactions with menstrual phase and with BC category. Extra sum-of-squares F-tests were performed to test for the significance of the main effect of cardiovascular strain and the interaction terms.

Results

The final dataset included 362 852 days of data representing 13 535 MCs (2.9±2.0 cycles per participant; 28.4±2.4 days per cycle) from 4594 respondents (natural MC n=3870, BC n=455, progestin-only n=269). The participants were 33.5 (±7.3) years, with a body mass index of 24.5 (±4.1), exercising 5.6 (4.1) hours per week, with a daily strain of 10.9 (±2.6) (table 1).

Table 1
|
Participant characteristics

Natural MC versus BC: HRV and recovery decreased across the MC from the mid-follicular phase whereas resting HR and RR increased. The effect of menstrual phase was found to be statistically significant at an alpha level of 0.0125 (p<0.001) for all measures. The patterning of these physiological measures is paradoxical across phases for those individuals using BC. HRV and recovery start elevated and drop off quickly during the withdrawal bleed (inactive pill week), rising through the active pill weeks (corresponding ‘week 2: late follicular’, ‘week 3: early luteal’ and ‘week 4: late luteal’ phases). Resting HR starts low during the first few days of the withdrawal bleed (inactive pill week) and quickly rises in the later days of the bleed week, then slowly decreases across the 3 weeks of active pill use. RR dips to its lowest point in the later days of the withdrawal bleed (corresponding early follicular), then increases across the active pill weeks, reaching its peak at the end of the fourth week of the active pill (corresponding late luteal phase). The effect of menstrual phase was found to be statistically significant at an alpha level of 0.0125 (p<0.001) for all measures. The interaction between BC category (Natural MC vs BC) and menstrual phase was also found to be statistically significant for all measures at an alpha of 0.0125 (p<0.001), suggesting that the use of BC modulates the patterning of these physiological measures over the cycle (figures 1 and 2).

Figure 1
Figure 1

Effects of natural versus BC menstrual phases on physiological measures. BC, birth control; HRV, heart rate variability. *denotes significance with α set at 0.003125; p<0.001.

Figure 2
Figure 2

Effects of natural versus BC menstrual phases on unnormalised physiological measures. BC, birth control; HRV, heart rate variability. *denotes significance with α set at 0.003125; p<0.001.

Natural MC versus Progestin only: For individuals using progestin-only BC, the effect of menstrual phase was similarly found to be statistically significant at an alpha level of 0.0125 (p<0.001) for all measures. Again, the patterning of these measures was found to be distinct from the patterning for the natural MC group. The interaction between hormonal BC category (natural MC vs progestin-only BC) and menstrual phase was found to be statistically significant at an alpha of 0.0125 (p<0.001) for all measures (figures 3 and 4).

Figure 3
Figure 3

Effects of natural versus progestin-only menstrual phases on physiological measures. BC, birth control; HRV, heart rate variability. *denotes significance with α set at 0.003125; p<0.001.

Figure 4
Figure 4

Effects of natural versus progestin-only BC menstrual phases on unnormalised physiological measures. BC, birth control; HRV, heart rate variability. *denotes significance with α set at 0.003125; p<0.001.

Cardiovascular strain and next-day recovery

There is evidence that the relationship between normalised recovery and previous day strain is modulated by BC type. Specifically, those in the HC groups exhibited steeper declines in recovery with additional strain-normalised recovery decreases by an additional 0.0055±0.00135 (p<0.001, significant at an alpha of 0.025) per unit of strain for individuals using combination BC, and an additional 0.0025±0.0019 (p=0.19, not significant at an alpha of 0.025) for progestin-only BC individuals (figure 5). The main effect of previous day strain was found to be statistically significant (p<0.001, alpha of 0.0125), the inclusion of the interaction between previous day strain and menstrual phase was found to be statistically significant (p<0.001, alpha of 0.0125), and the inclusion of the interaction between previous day strain and BC category was found to be statistically significant (p<0.001, alpha of 0.0125).

Figure 5
Figure 5

Modulation of cardiovascular strain and next day recovery by hormonal contraception. BC, birth control.

Discussion

The ANS is essential for regulating heart rate and arterial pressure via the arterial baroreflex. The activity of the autonomic nerves is controlled centrally via nuclei which relay excitatory or inhibitory information to the preganglionic autonomic fibres in the sympathetic and/or parasympathetic pathways.8 16 Sex hormones directly affect the activity of these pathways, thus altering peripheral sympathetic and parasympathetic neural activity, and subsequently, cardiovascular function.6 26

HRV analysis is a valuable, objective marker of ANS used to study an athlete’s training vs recovery equilibrium and detect early overreaching state that can decrease the athlete’s performance.11–13 19 24 To our knowledge, existing research and its application to female athletes has been generalised from male data, despite limited research indicating significant sex differences and responses of HRV and the ANS.5 12 13 This generalisation is a significant misstep for female athletes as the neural control of circulation differs with sex, mainly attributed to the effects of ovarian hormones on ANS.8 16 23 25 26 Moreover, endogenous and exogenous hormones exert different modulatory effects; specific to the progestogen of the exogenous hormone formulation.16 17 26 27 The current study is the first to investigate HRV, RHR, RR and sleep metrics in a large population of female athletes; to include naturally cycling women, women using COC and women using progestin-only contraception.

Endogenous versus exogenous hormone modulatory effects

In this robust dataset, it was observed that the patterning of ANS modulation from ovarian hormones is significantly different between naturally cycling women and those on BC; with the patterning dependent on the type of BC used.

Specifically, in naturally cycling women, HRV and recovery are elevated in the early and mid-follicular phases with a significant decline of HRV and recovery metrics into the late luteal phase. As a brief review, the first half of the MC comprised by the menstrual and follicular phases during which time oestrogen levels are low (menstrual phase) and rise (follicular phase) and ends with the periovulatory phase in which follicular-stimulating hormone and luteinising hormones peak. The second half of the cycle is comprised by the luteal (during which time oestrogen level rises with a progesterone peak) and the premenstrual phases during which time oestrogen and progesterone levels fall (figure 6). It is known that oestrogens act centrally to modulate the autonomic nervous system, increasing vagal and decreasing sympathetic activity, whereas progesterone appears to have an opposing effect, elevating central norepinephrine release inducing a sympathetic drive.24 27 Thus, due to the unopposed oestrogen in the early and mid-follicular phases, followed by the elevation of progesterone in the luteal phase which antagonises oestrogenic effects, it is unsurprising to note this pattern of HRV and recovery in naturally cycling women.

Figure 6
Figure 6

Diagram of the menstrual cycle. FSH, follicular-stimulating hormone; LH, luteinising hormone.

However, in women using combined BC, a different, paradoxical pattern emerges. We observed an elevated HRV and recovery in the withdrawal phase (inactive pill phase), with a significant decline at the onset of active pill use. This patterning may be representative of the perturbation of exogenous hormones. In the withdrawal phase of BC use, exogenous hormones are not consumed, thus their effect on the downregulation of the HPA axis is removed; given the different half-lives of the exogenous steroids and variable impact on the endogenous hormones, the withdrawal week should be considered a transient hormonal profile phase.1 28 With resuming exogenous hormone consumption, a stable hormone profile emerges; the inhibitory effects of both the oestrogen and the progestin components are established through synergistic interactions at the hypothalamic–pituitary level. Recently, Assadpour et al29 suggested the chronic upregulation of β-adrenergic receptors of BC users may play a role in the upregulation of ventilatory rate, and autonomic reflex function, which is reflected in the patterning of the women on BC of this study of decreased HRV, increased RR, increased RHR across the active pill phases.

Interestingly, women taking the progestin-only contraception exhibited a similar pattern as naturally cycling women except for the late luteal phase in which there was an observed attenuation of sympathetic drive. This may be due to the generation of the progestogens used in the progestin-only formulations. The progestin-only formulations use only norethisterone, a first-generation progestin, or desogestrel, a third-generation progestin, which are 1.5 times more potent than endogenous progesterone,27 thus may attenuate the magnitude of the phase-based increase of sympathetic drive in the later stages of the cycle.

Cardiovascular strain and next day recovery

The application of appropriate training load is one of the fundamental factors for positive physiological adaptations to occur with improvements in performance. Monitoring training loads and recovery should reduce the incidence of excessive loading, which results in negative adaptations (including non-functional overtraining).30 31 A major tenet in determining load application is based on cardiovascular strain; periods of higher cumulative stress are typically highlighted by a reduction in HRV and in situations of positive adaptations to training, HRV should increase, or remain stable, with increasing training loads. The data from the current study indicates a reduced recovery in women using combined BC with every unit of increased strain, as compared with naturally cycling or progestin-only pill women across regardless of cycle phase. In addition to the exogenous effects described above, recovery may have been attenuated due to greater inflammatory responses associated with BC use32–34 and to the difference in stress responsiveness of women on BC. Larsen et al34 reported a higher incidence of IL-6 and C reactive protein (CRP) in Olympic-calibre female athletes using BC as compared with the same calibre naturally cycling female athlete. Moreover, Cauci et al33 observed a marked elevation of low-grade elevation (as measured by CRP) in female athletes using BC, which may exacerbate inflammatory responses to physical stress, increasing the time course for recovery. Of greater interest as it pertains to the ANS is the difference between naturally cycling and BC users free cortisol. It has been repeatedly shown that BC increase free circulating cortisol in the blood35 36 which is associated with chronically elevated cortisol levels and a reduced diurnal rhythm.37 As cortisol is a known stimulatory hormone, increasing the threshold for parasympathetic drive, this may be a mechanism responsible for decreased recovery metrics per unit of strain.

Limitations and strengths

Strengths of this project include the large dataset of physically active premenopausal women; to include normally cycling, COC users, and progestin-only pill users; allowing for robust comparisons which have not been identified in previous literature to the authors knowledge.

This project has some limitations inherent to the data retrieved. One limitation is the classification of the MC phase. Determination of the MC phase was somewhat uncertain as the model assumed fertility in a regular 28-day cycle, which may not correspond to the phases determined by modern hormonal documentation or the presence of anovulatory cycles. As this was a large retrospective study, MC phases were not confirmed with blood tests (as is the gold standard,1 which does not allow for identification of extended follicular phases or luteal phase deficiencies. Second, identification of the generation of the progestin in the oral contraceptive pills was not feasible except for the progestin-only (as there are only two progestins available to market for the mini-pill). The mechanism of progestins is a multimodal one, and each subsequent generation of progestogen has a different binding affinity to steroid receptors and serum binding proteins27 which may affect the threshold of parasympathetic downregulation. Third, due to privacy and deidentification of data, specific identifiers such as ethnicity, socioeconomic status, medication usage, cardiovascular and metabolic disease risk factors and other dietary-lifestyle choices were not included in the data aggregation and analysis. Further research lends to the investigation of these different progestin generations on parasympathetic thresholds and vagal tone. To note, a unique limitation, in this era of COVID-19, is the psychological stress of the pandemic which may have influenced HRV and recovery metrics.

Conclusion

As the number of female athletes competing across all levels rises, training methodologies should evolve to reflect sex differences across critical metrics of adaptation to training. Using a robust data set (13 535 MCs from 4594 respondents) we presented novel data with regards to recovery metrics across three different hormone profiles of female athletes: naturally cycling, COC and progestin-only users. Understanding how endogenous and exogenous hormones modulate the autonomous nervous system and subsequent recovery metrics is critical for monitoring training loads and the application of appropriate training stress.