Front. Physiol. Frontiers in Physiology Front. Physiol. 1664-042X Frontiers Media S.A. 10.3389/fphys.2020.00629 Physiology Original Research Oral Contraceptive Use Influences On-Kinetic Adaptations to Sprint Interval Training in Recreationally-Active Women Schaumberg Mia Annalies 1 2 3 * Stanley Jamie 4 5 6 Jenkins David G. 1 3 Hume Emily A. 7 Janse de Jonge Xanne A. K. 8 Emmerton Lynne M. 9 Skinner Tina L. 3 1School of Health and Sport Sciences, University of the Sunshine Coast, Sippy Downs, QLD, Australia 2Sunshine Coast Health Institute, Birtinya, QLD, Australia 3School of Human Movement and Nutrition Sciences, The University of Queensland, Brisbane, QLD, Australia 4South Australian Sports Institute, Kidman Park, SA, Australia 5Australian Cycling Team, Gepps Cross, SA, Australia 6Allied Health and Human Performance, University of South Australia, Adelaide, SA, Australia 7School of Clinical Medicine, The University of Queensland, Brisbane, QLD, Australia 8School of Environmental and Life Sciences, The University of Newcastle, Ourimbah, NSW, Australia 9School of Pharmacy and Biomedical Sciences, Curtin University, Perth, WA, Australia

Edited by: Trine Moholdt, Norwegian University of Science and Technology, Norway

Reviewed by: Sarah J. Willis, Université de Lausanne, Switzerland; Zhaowei Kong, University of Macau, Macau

 *Correspondence: Mia Annalies Schaumberg, mschaum1@usc.edu.au

This article was submitted to Exercise Physiology, a section of the journal Frontiers in Physiology

 12 06 2020 2020 11 629 03 02 2020 18 05 2020 Copyright © 2020 Schaumberg, Stanley, Jenkins, Hume, Janse de Jonge, Emmerton and Skinner. 2020 Schaumberg, Stanley, Jenkins, Hume, Janse de Jonge, Emmerton and Skinner

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

 Introduction

Oral contraceptive (OC) use influences peak exercise responses to training, however, the influence of OC on central and peripheral adaptations to exercise training are unknown. This study investigated the influence of OC use on changes in time-to-fatigue, pulmonary oxygen uptake, cardiac output, and heart rate on-kinetics, as well as tissue saturation index to 4 weeks of sprint interval training in recreationally active women.

 Methods

Women taking an oral contraceptive (OC; n = 25) or experiencing natural menstrual cycles (MC; n = 22) completed an incremental exercise test to volitional exhaustion followed by a square-wave step-transition protocol to moderate (90% of power output at ventilatory threshold) and high intensity (Δ50% of power output at ventilatory threshold) exercise on two separate occasions. Time-to-fatigue, pulmonary oxygen uptake on-kinetics, cardiac output, and heart rate on-kinetics, and tissue saturation index responses were assessed prior to, and following 12 sessions of sprint interval training (10 min × 1 min efforts at 100–120% PPO in a 1:2 work:rest ratio) completed over 4 weeks.

 Results

Time-to-fatigue increased in both groups following training (p < 0.001), with no difference between groups. All cardiovascular on-kinetic parameters improved to the same extent following training in both groups. Greater improvements in pulmonary oxygen up-take kinetics were seen at both intensities in the MC group (p < 0.05 from pre-training) but were blunted in the OC group (p > 0.05 from pre-training). In contrast, changes in tissue saturation index were greater in the OC group at both intensities (p < 0.05); with the MC group showing no changes at either intensity.

 Discussion

Oral contraceptive use may reduce central adaptations to sprint interval training in women without influencing improvements in exercise performance - potentially due to greater peripheral adaptation. This may be due to the influence of exogenous oestradiol and progestogen on cardiovascular function and skeletal muscle blood flow. Further investigation into female-specific influences on training adaptation and exercise performance is warranted.

 female training adaptation ovarian hormones oral contraceptive (OC) exogenous hormones cardiorespiratory time-to-fatigue Introduction

Exogenous hormones introduced through oral contraceptive (OC) use can reduce maximal exercise capacity (Notelovitz, 1987; Casazza et al., 2002; Lebrun et al., 2003), increase fat-mass (Berenson and Rahman, 2009) and change the metabolic (Isacco et al., 2012), thermoregulatory (Stachenfeld et al., 2000), cardiovascular (Coney et al., 2001), and ventilatory (Charkoudian and Joyner, 2004) responses to exercise. Although OC use has been shown to impede peak exercise adaptations to training (Schaumberg et al., 2017a), the physiological mechanisms and whether performance adaptations are influenced by OC use is unclear (Casazza et al., 2002; Lebrun et al., 2003; Zierke, 2010).

One potential mechanism mediating training adaptation is skeletal muscle blood flow. During incremental exercise, skeletal muscle blood flow is determined by locally induced vasodilation and sympathetically mediated vasoconstriction, both of which are influenced by ovarian hormones (Charkoudian and Johnson, 1997; Charkoudian and Joyner, 2004; Vanheest et al., 2005). Indeed, chronic estrogen exposure has known vasodilatory responses (Stathokostas et al., 2009), and there is some evidence to suggest that estrogen and progestogen supplementation and the menstrual cycle may influence blood flow during exercise (Gurd et al., 2007), the directionality of which appears to depend on multiple factors, including the phase of the menstrual cycle (i.e., the ratio of estrogen and progesterone) and the type and concentration of exogenous ovarian hormone concentrations. Altered muscle blood flow due to endogenous and exogenous ovarian hormones may influence the ability of the muscle to meet oxidative demands during exercise. However, this has not been investigated, nor has OC use been considered as a potential mediator of skeletal muscle blood flow adaptations to exercise training in women.

A further potential mechanism is the integration of pulmonary and cardiovascular systems to deliver oxygenated blood to skeletal muscle during exercise. Pulmonary oxygen uptake on-kinetics (τV̇O2p) provide insight into how the cardiovascular system and mitochondria integrate to increase aerobic energy production in response to exercise (Murias et al., 2011). The speed of τV̇O2p is a good indicator of endurance performance (Jones and Burnley, 2009); a faster τV̇O2p indicates earlier achievement of physiological steady state (indicated by the time constant, τon, or the mean response time, MRT) resulting in reduced oxygen deficit (Xu and Rhodes, 1999). Additionally, a fast τV̇O2p has also been associated with reduced lactate accumulation and muscle glycogen depletion compared to a slow τV̇O2p (Berger et al., 2006). There is increasing interest in optimal training methods to elicit physiological adaptations within pulmonary, cardiovascular, and muscular systems to improve exercise performance.

Sprint interval training (SIT) can elicit adaptations traditionally associated with endurance training in a shorter period (Gibala et al., 2012; Shiraev and Barclay, 2012). These adaptations appear to be independent of sex and include improved oxidative enzyme activity (Carter et al., 2001; Astorino et al., 2011) coupled with increased capillarization (Shiraev and Barclay, 2012) and more efficient blood distribution (Murias et al., 2011), which can lead to improvements in V̇O2peak. These adaptations also improve the rate at which oxygen is extracted in the lungs (τV̇O2p) (Shoemaker et al., 1996; Tschakovsky and Hughson, 1999; Grassi, 2001; Hughson et al., 2001; Gibala et al., 2006). Despite these known physiological adaptations, the rate of oxygen extraction in muscle (represented by change in deoxyhaemoglobin: Δ[HHb]) or tissue saturation index (TSI), is not usually influenced by short training interventions (i.e., <6 weeks) (Overend et al., 1992; Berger et al., 2006; McKay et al., 2009); whether SIT can elicit adaptations at the muscular level following a shorter period of training is inconclusive.

The direct relationship between pulmonary and muscle oxygen extraction is important in understanding central (i.e., adaptations to cardiorespiratory function rather than adaptations to the peripheral vasculature and trained muscle) and peripheral (i.e., adaptations within skeletal muscle such as capillarization and/or mitochondrial biogenesis) adaptations to SIT; both parameters provide a measure of endurance capacity in recreationally active individuals (Murias et al., 2011; Spencer et al., 2013). Research in men suggests that an improvement in τV̇O2p with no concurrent change in TSI following a period of exercise training indicates that increased muscle oxygen utilization is accompanied by faster muscle oxygen extraction (MacPhee et al., 2005; McKay et al., 2009; Murias et al., 2011; Spencer et al., 2013). While research in women has found similar effects, hormone status or OC use has not been previously considered (Talanian et al., 2007; Astorino et al., 2011; Murias et al., 2011). Due to the potential influence of OC use on cardiovascular function and skeletal muscle blood flow due to chronic exogenous oestradiol and progestogen exposure, as well as our previous finding that OC use dampens V̇O2peak adaptations to SIT in women (Schaumberg et al., 2017a), the investigation of whether OC use influences pulmonary (τV̇O2p) and muscular oxygen extraction (TSI), as well as associated cardiovascular adaptation to SIT is warranted.

Therefore, the aim of the present study was to investigate the influence of OC use on pulmonary, cardiovascular and muscular oxygen uptake kinetics adaptation at moderate and heavy exercise intensities following 4 weeks of SIT.

 Materials and Methods Overview

Physically active women with either regular menstrual cycles (MC; no current hormonal contraception) or using an OC completed two exercise tests – an incremental exercise test and a square-wave step-transition protocol (separated by a minimum of 48 h), prior to, and following a 4-week SIT program.

 Participants

Healthy, recreationally active women (regularly completing at least 150 min of self-reported moderate to vigorous physical activity per week, but not currently training for, or competing at state or national level sport competition), who were either long-term (minimum 6 months uninterrupted) monophasic combined OC users (n = 25) or experiencing regular natural menstrual cycles (MC; n = 22) participated in the study. All experimental procedures were approved by an ethics committee of The University of Queensland and participants provided written informed consent.

 Control Measures

The procedures relating to hormone verification and analysis and body composition assessment are described in detail by Schaumberg et al. (2017b). To summarize, all OC users completed testing in the ‘active pill’ phase of the oral contraceptive cycle, and all naturally menstruating women completed testing in the mid-luteal phase of the menstrual cycle with serum hormone verification conducted at each timepoint. Nutrition, hydration, and exercise control measures have also been previously described (Schaumberg et al., 2017a). In brief, prior to each trial participants completed a 24-h food diary, fasted overnight, consumed a standardized moderate carbohydrate pretrial meal 1 h before arrival at the laboratory, abstained from caffeine, alcohol, and other stimulants and depressants for 24 h, recorded any additional medications or supplements, and maintain an euhydrated state. Participants were encouraged to maintain their normal physical activity levels throughout the study; however, were asked to refrain from strenuous physical activity for 24 h before each trial and arrive at the laboratory in a rested state. A pretrial preparation checklist was completed to confirm compliance to pretesting requirements.

 Experimental Protocol

In each of the experimental trials, participants completed 3, 4-min step transitions to a moderate exercise intensity [90% of power output at ventilatory threshold (POVT); calculated as 90% of the power output (PO) in Watts (W) achieved at ventilatory threshold (VT)] and three, 3-min transitions to a high exercise intensity [Δ50% POVT; calculated as PO at VT plus 50% of the difference between PO at VT and peak power output (PPO)], as determined from the incremental exercise test (Joyce et al., 2013; Stanley et al., 2014) (Figure 1). The first transition was preceded by 4 min of baseline cycling (20 W) and each subsequent transition was separated by 4 min of baseline cycling (20 W). Participants were instructed to maintain a consistent cadence (70 ± 10 RPM) throughout the baseline, moderate, and high intensity cycling. Heart rate (HR) (Suunto®, United States), expired air (Parvo Medics’ TrueOne® 2400 Indirect Calorimetry System, Utah, United States), cardiovascular parameters (PhysioFlow, Manatec Biomedical, France), and muscle oxygenation via near infra-red spectroscopy (NIRS; Portalite, Artinis Medical Systems BV, Netherlands) were continuously monitored throughout the trial, and RPE (Borg, 1973) was recorded at the end of every step-transition.

Square-wave step transition exercise protocol to moderate and heavy exercise intensities. W, watts; POVT, power output at ventilatory threshold; min, minute; TTF, time to fatigue.

 Measures Respiratory Measures

A familiarization session was completed prior to the first experimental trial. The V̇O2peak protocol involved continuous incremental (25 W⋅min–1) exercise test on an electronically braked cycle ergometer (Lode Excalibur Sport, Lode BV, Groningen, Netherlands) as previously described (Schaumberg et al., 2017a). Respiratory gas exchange was continuously recorded via automated indirect calorimetry (Parvo Medics’ TrueOne® 2400 Indirect Calorimetry System, Utah, United States) for calculation of ventilatory parameters. Before each test, the analyzers were calibrated in accordance with the manufacturers’ recommendations. From the incremental test, data were averaged in 15-s epochs; V̇O2peak was defined as the highest V̇O2 value attained during a 15-s period (Rossiter et al., 2006; Stanley et al., 2014). During the experimental trials, average V̇O2 was determined from 5-s interval data.

 Cardiovascular Measures

During exercise, heart rate (HR), stroke volume (SV) and cardiac output (Q̇) were measured continuously using impedance cardiology (PhysioFlow®, Manatec Biomedical, France) (Charloux et al., 2000; Richard et al., 2001). The theoretical basis for determining cardiac output from this method and its validity during rest and exercise has been described previously (Charloux et al., 2000; Lepretre et al., 2004), and has been validated against the direct Fick method (Richard et al., 2001). Two sets of electrodes (Skintact FS-50, Leonhard Lang Gmbh, Austria) – one transmitting, one sensing – were applied above the supra-clavicular fossa at the left base of the neck, and along the xiphoid process. Another two electrodes were used to monitor a single electrocardiographic signal (ECG; CM5 position). Blood pressure was assessed (Digital blood pressure monitor, UA-767, A&D Instruments Ltd., United Kingdom) as part of standard calibration process for the PhysioFlow® prior to the exercise test. During the experimental trials (Figure 1), HR, SV, and Q̇ data were sampled at 5-s intervals, with the average for each interval determined during moderate and high loads. The coefficient of variation of cardiac output measures during peak exercise using this method has been reported as 3.4–3.6% (Hsu et al., 2006).

 Near-Infrared Spectroscopy Measurements

Near-infrared spectroscopy (NIRS) (Portalite, Artinis Medical Systems BV, Netherlands) estimated the oxygenation of the right vastus lateralis muscle during the performance trial. A three-wavelength continuous wave system was used, which simultaneously used the modified Beer–Lambert and spatially resolved spectroscopy methods (Ihsan et al., 2013; Stanley et al., 2013). Changes in total hemoglobin (tHb) oxyhaemoglobin (O2Hb) and deoxyhaemoglobin (HHb) were measured using the differences in absorption characteristics of light at 775, 810, and 850 mm (Stanley et al., 2013). An arbitrary value for the differential path length of 3.83 mm was used due to the uncertainty of proton path length at rest and during exercise (Stanley et al., 2013).

The NIRS device was connected via Bluetooth to a computer acquiring data at 10 Hz for later analysis. The probe was positioned one-third of the way along the vertical length of the thigh (from the quadriceps tendon on the patella). The NIRS device was wrapped in a zip-locked bag for water-proofing and covered in black material to prevent slipping and interference from ambient light and strapped to the leg securely with a bandage. Changes in tHb, O2Hb, and HHb were reported as a change from baseline measures after first recording the tissue saturation index (TSI). The balance between oxygen consumption of the muscle tissue and oxygen supply (Equation 1) represented TSI (as a percentage).

TSI = [ O 2 ⁢ Hb ⁢ consumed ] / [ O 2 ⁢ Hb ⁢ supplied ] × [ HHb ] × 100

Equation 1: TSI expressed as changes in O2Hb and HHb (Stanley et al., 2013). TSI, tissue saturation index; O2Hb, oxygenated hemoglobin; HHb, deoxygenated hemoglobin.

 Sprint Interval Training Protocol

As previously described (Schaumberg et al., 2017a), the SIT protocol required participants to complete three supervised SIT sessions per week for 4 weeks, with a minimum of 36 h between sessions. Each session involved 1 min of work followed by 2 min of passive recovery in a 1:2 work:rest ratio (Ready et al., 1981; Gosselin et al., 2012). The work interval intensity was self-selected at the maximal sustainable effort between 100 and 120% of PPO determined in the V̇O2peak test. Participants completed 10 1-min repetitions per session and peak heart rate, RPE, average power output and PPO were recorded. All exercise sessions were completed on an air- and magnetically braked cycle ergometer (Wattbike Ltd., Nottingham, United Kingdom).

 Data Analysis Assessment of Pulmonary Oxygen Uptake Kinetics

Recorded data for V˙O2 during each 4-min moderate (90% POVT) and 3-min high (Δ50% POVT) (plus an additional 1 min prior to each interval to determine baseline V˙O2) loads, were interpolated into 5 s intervals; aberrant data points (caused by swallowing and coughing) were filtered out (Ozyener et al., 2001; Berger et al., 2006). To decrease the signal-to-noise ratio caused by high variability between breaths, the common practize of including multiple exercise transitions into the same protocol was employed (Whipp et al., 2005; Stanley et al., 2014). The first interval of each intensity was excluded from analysis due to the lack of a priming effect (Whipp et al., 2005; DeLorey et al., 2007; Stanley et al., 2013). Therefore, V˙O2 recorded during the second and third intervals for each intensity of each experimental trial was time synchronized and ensemble averaged to yield a single response for each participant for each trial. A repeated iterative technique (Sigmaplot 10, SPSS Science; Chicago, IL, United States) and a mono-exponential function (Equation 2) (Ozyener et al., 2001; Dorado et al., 2004; Whipp et al., 2005; Stanley et al., 2014) were used to model τV˙O2p over moderate and high loads.

V . ⁢ O 2 ⁢ p ⁢ ( t ) = V . ⁢ O 2 ⁢ baseline + Ampl × [ 1 - e - ( τ - TD / τ on ) ] × U 1

Equation 2: Mono-exponential function (Sigmaplot 10, SPSS Science; Chicago, IL, United States) where u = if {t < TD1, VO2b, VO2b + A1∗[1−exp(−(t-TD1/tau1)]}. Initial parameters V˙O2b = 0.5 TD1 (time delay) = −0.5; A1 (amplitude) = 1.5; Tau1 = 30. Constraints VO2 > 0; A1 > 0; Tau1 > 0; TD1 ≥ 0. V˙O2p, pulmonary oxygen kinetics; V˙O2, oxygen consumption; Ampl, amplitude; TD, time delay.

In Equation 2 V˙O2 baseline is the average V˙O2 during the 60 s prior to onset of the rest- (or active recovery) to-exercise transition, Ampl is the asymptotic amplitude for the exponential term and τon is the time constant of the exponential (seconds) (Stanley et al., 2014). Use of the mono-exponential function was most appropriate due to its simplicity when considering the rare occurrence of τV˙O2p phase III (due to the priming effect from the first interval, the short duration of the interval and the submaximal intensity of exercise, therefore if Phase III were to occur the amplitude would be minor) (Stanley et al., 2014). Phase I to Phase II transition occurred approximately 15 s after the onset of exercise for all participants (determined by visual examination), therefore this initial cardiodynamic component was excluded by deleting the first 20 s of data (Ozyener et al., 2001; Dorado et al., 2004; Stanley et al., 2014). Overall τV˙O2p, tau1 and mean response time (MRT = time delay + tau1) were calculated and very low residuals (r2 > 0.98) obtained (Stanley et al., 2014), providing an overall description of on-transient oxygen uptake kinetics. While tau1 and MRT are closely related, MRT has been shown to be more reliable that tau1 in recreationally active women (Schaumberg et al., unpublished data), and therefore MRT was used as the primary measure of interest.

 NIRS Data and Assessment of De-Oxygenation Rates

Data from the NIRS device (tHb, O2Hb, and [HHb]) were sampled down from 10 to 1 Hz and then averaged into 5 s intervals. Since tissue saturation index has been shown to provide a more accurate indication of muscle oxygenation status than Δ[HHb] (Wolf et al., 2007), TSI data was modeled. Using the Sigmaplot 10 program (SPSS Science; Chicago, IL, United States), a linear model was used to calculate the results from the average of the second and third intervals for both moderate and high intensity. Similar to τV˙O2p, the first interval was excluded from analysis due to the priming effect (Stanley et al., 2013). Moderate exercise prior to heavy workloads has been shown to influence Δ [HHb] (Whipp et al., 2005; Spencer et al., 2013), therefore, both workloads were analyzed using a linear model (Equation 3) (Bae et al., 2000):

TSI = a × t + b

Equation 3: Linear equation (Sigmaplot 10, SPSS Science; Chicago, IL, United States). TSI, tissue saturation index; t, time; a, slope; b, y-intercept.

TSI was modeled without time delay during the first 20 s of moderate intensity and 30 s of high intensity. However, the slope (a) was retained as an index of deoxygenation rate (Stanley et al., 2013). This particular model yielded very low residuals (r2 > 0.98) and provided an overall descriptor of the muscular deoxygenation rate, and therefore muscle deoxygenation kinetic response to exercise.

 Assessment of Cardiovascular Kinetics

Heart rate and Q̇ on-transient kinetics were modeled using the same iterative technique adopted for V˙O2 on-transient kinetics. HR and Q̇ data were fitted with a mono-exponential function consistent with Equation 2 (using the same 4 or 5 min window), with the HR and Q̇ data substituted for V˙O2 (Stanley et al., 2014). Unlike the V˙O2 on-transient kinetics analysis, the initial 20 s of data was not deleted due to lack of a cardiodynamic (Phase I to Phase II) transition. The mean response time (MRT = time delay + tau1) was calculated to provide an overall description of on-transient cardiovascular kinetics.

 Statistical Analysis

Data were analyzed using Microsoft Excel® 2007 and SPSS® (version 22.0, SPSS, Inc., Chicago, IL, United States). Normality of distribution was tested using the Kolmogorov–Smirnov test; when not normally distributed, data were log-transformed and re-checked for normality of distribution. Analyses included standard descriptive statistics, paired t-test, and two-way repeated measures analysis of variance (ANOVA) (with a main effect for training × group). To locate the source of significant differences, the Bonferroni post hoc test was used. Homogeneity of variance was confirmed using Mauchly’s test of sphericity. When the assumption of sphericity was violated (p < 0.05), the F-statistic was adjusted using the Greenhouse–Geisser correction. Where Mauchly’s test of sphericity was not found to be significant, post hoc analyses assumed sphericity (Vincent, 1999). Magnitude-based inferences (Hopkins et al., 2009; Batterham and Hopkins, 2015) calculated the between-trial standardized differences or effect sizes [ES, 95% confidence interval (CI)] using the pooled standard deviation (Cohen, 1988) and standard threshold values (Batterham and Hopkins, 2005). All tests were two-tailed and statistical significance was set at p < 0.05. Parametric results are given as the mean, standard deviation and 95% confidence interval (CI), [mean ± SD (95% CI)]; non-parametric results are given as the median and interquartile range and 95% CI, [median (IQR) (95% CI)] unless stated otherwise.

 Results Participant Characteristics, Control Measures, and Training Protocol

Participant recruitment and retention has previously been described (Schaumberg et al., 2017a). Due to the nature of the outcome measures, we included all naturally menstruating participants within the main data set, including six of the 22 MC participants who had ovulatory, regular menstrual cycles but exhibited potential luteal phase deficiency (LPD) based on failing to meet the serum progesterone criterion of >6 ng.mL–1 (Schaumberg et al., 2017b). Supplementary Table S1 comparing the normal MC (n = 16) versus the potential LPD MC (n = 8) participants has been included; no significant differences between outcome measures pre- or post-training were found, though the MC LPD group showed dampened pulmonary oxygen up-take kinetic responses to training, discussed below. All 25 participants recruited to the OC group were taking a low-dose, monophasic combined oestradiol and progestin formulation. There were variations in androgenic (n = 5), antiandrogenic (n = 5) and non-androgenic (n = 15) formulations [calculated using the method of Greer et al. (2005)]; subsequent analyses confirmed androgenicity of OC type did not influence baseline characteristics or outcome measures. Physical activity, energy intake and body composition parameters were not different within or between groups at any timepoint, however, due to the inclusion of the potential LPD participants (n = 8) there were some differences in hormone concentrations between groups that have been previously described (Schaumberg et al., 2017b), with the MC group having higher oestradiol, progesterone and free androgen index (p < 0.001), and lower sex-hormone binding globulin (SHBG; p < 0.001). As such, participant characteristics are re-presented in Table 1.

Participant characteristics.

Oral contraceptive group (n = 25)
 Menstrual cycle group (n = 22)
Pre-training Post-training Pre-training Post-training
Participant demographics and control measures
Age (years) 25.5 ± 5.4 (23.1–27.8) X 26.4 ± 5.2 (24.0–28.8) X
MC length (days) 28.0 ± 0.0 (28.0–28.0) X 30.5[28.8–33.3]* (29.9–32.6) X
Physical activity (min.wk–1) 247 ± 64 (222–272) 235 ± 61 (211–258) 235 ± 58 (211–260) 222 ± 51 (201–243)
Energy intake (kJ.kg–1.d–1) 8461 ± 3194 (6896–10026) 8490 ± 2452 (7103–9877) 8373 ± 2360 (7179–9567) 8499 ± 1971 (7269–9712)
Body composition
Body mass (kg) 63.6 ± 7.8 (60.3–66.8) 63.4 ± 7.2 (60.4–66.3) 64.6 ± 9.2 (60.5–68.8) 64.2 ± 8.8 (60.2–68.2)
Body mass index (kg.m–2) 22.6 ± 2.1 (21.7–23.4) 22.6 ± 2.1 (21.7–23.4) 22.7 ± 2.3 (21.7–23.7) 22.7 ± 2.3 (21.7–23.7)
Fat mass (kg) 20.7 ± 4.8 (18.7–22.7) 20.5 ± 4.6 (18.6–22.4) 21.0 ± 5.2 (18.8–23.3) 20.8 ± 5.6 (18.2–23.3)
Lean body mass (kg) 40.6 ± 4.4 (38.8–42.5) 40.6 ± 4.1 (38.9–42.3) 41.3 ± 5.9 (38.7–43.9) 41.2 ± 5.2 (38.8–43.5)
Lean body mass – legs (kg) 13.6 ± 1.8 (12.9–14.4) 13.8 ± 1.9 (13.0–14.5) 13.8 ± 2.5 (12.7–14.9) 13.8 ± 2.3 (12.7–14.8)
Body fat (%) 32.3 ± 4.8 (30.3–34.3) 32.1 ± 4.8 (30.1–34.0) 32.3 ± 5.3 (30.0–34.7) 32.0 ± 5.9 (29.3–34.7)
Hormone measures
Estradiol (pg.mL–1) 5.6 [5.0–10.3] (5.7–13.5) 5.1[5.0–10.2] (11.0–63.8) 124.4 ± 67.5Λ (92.8–156.1) 93.2 ± 67.5#Λ (62.4–123.9)
Progesterone (ng.mL–1) 0.6 ± 0.3 (0.4–0.7) 0.5 ± 0.3 (0.4–0.6) 10.0 ± 7.9Λ (6.3–13.7) 1.0[0.7–1.4]# (0.2–7.0)
Total testosterone (ng.mL–1) 0.20 ± 0.10 (0.10–0.20) 0.13 ± 0.07 (0.10–0.16) 0.2[0.1–0.4] (0.2–0.3) 0.3 ± 0.2Λ (0.2–0.4)
SHBG (pg.mL–1) 209.0 ± 87.1 (170.3–247.6) 189.6 ± 99.3 (145.6–233.7) 62.6[40.8–90.0]Λ (50.7–86.3) 58.5[35.6–81.9]Λ (47.5–76.9)
Free androgen index (%) 7.8[4.2–8.5] (6.9–14.7) 8.6 ± 5.6 (6.1–11.0) 33.2[14.7–77.7]Λ (22.3–99.4) 55.3[19.6–96.5]Λ (35.3–106.0)
Parametric data are presented as mean ± SD (95%CI); non-parametric data are presented as median [IQR] (95% CI). #p < 0.05 vs. pre-training; *p < 0.001 vs. pre-training; Λp < 0.05 vs. OC group. Peak Exercise Responses and Time to Fatigue

Peak exercise adaptations have previously been reported (Schaumberg et al., 2017a), with the OC group showing dampened peak exercise adaptations to SIT compared to the MC group. There was no significant difference in time-to-fatigue (TTF) between groups at any timepoint (p > 0.05) (Figure 2). Following training, TTF was increased from pre-training in both groups (p < 0.001). Standardized between-group differences for within-group changes (Cohen’s d) demonstrated that there was a likely higher TTF adaptation to training in the OC-group compared to the MC-group (0.96 ± 1.04).

Time to fatigue (seconds) at a heavy intensity (Δ50%POVT) for the oral contraceptive group (n = 25) and menstrual cycle group (n = 22) pre- and post-training. Data are mean ± SD. *p < 0.001 vs. pre-training; paired samples t-test.

 Pulmonary Oxygen Uptake Kinetics

There was no significant difference between groups for τV˙O2p [expressed as the mean response time (MRT) in seconds] at any time point at moderate or high intensity (p > 0.05). Following training, τV˙O2p was improved in the MC-group at both moderate (p = 0.021) and high (p = 0.015) intensities but the OC-group showed no change from baseline at both intensities (p > 0.05). When the MC-group was sub-grouped for potential LPD, the LPD group showed dampened τV˙O2p at both intensities (p < 0.05) as measured by MRT; but there was no significant group × time interaction between normal MC and LPD MC (see Supplementary Table S1). There was a significant group x time interaction for τV˙O2p at moderate intensity (p = 0.020); the interaction was approaching significance at high intensity (p = 0.086). Post hoc analyses suggested that the MC-group showed greater improvement in τV˙O2p at high intensity following training compared to the OC-group (OC-group −1.4 s vs. MC-group −6.1 s; p = 0.021); but this was not significant between groups at moderate intensity (OC-group +0.5 s vs. MC-group −4.4 s; p = 0.097). Data are presented in Table 2 and the τV˙O2p profile at both intensities from a representative participant is presented in Figure 3. Standardized between-group differences for within-group changes (Cohen’s d) demonstrated that the OC-group had a likely lower τV˙O2p adaptation to training compared to the MC-group at moderate (−0.54 ± 0.67) and high (−0.57 ± 0.49) intensities.

Pulmonary oxygen uptake, cardiac output, and heart rate on-kinetic responses to moderate and heavy intensity exercise, pre- and post-training in oral contraceptive users (n = 25) and naturally menstruating women (n = 22).

Oral contraceptive group (n = 25)
 Menstrual cycle group (n = 22)
Pre-training Post-training Pre-training Post-training
Moderate intensity exercise (90% POVT)
Pulmonary oxygen uptake on-kinetic response
Baseline V˙O2 (L.min–1) 0.8 ± 0.1(0.7–0.8) 0.7 ± 0.1(0.6–0.7)# 0.7 ± 0.1(0.7–0.8) 0.7 ± 0.1(0.6–0.7)
Amplitude (L.min–1) 0.8 ± 0.2(0.7–0.9) 0.8 ± 0.2(0.7–0.9) 1.0 ± 0.5(0.7–1.2) 0.7 ± 0.4(0.6–0.9)#
Time delay (sec) 23.2 ± 7.3(20.1–26.3) 27.1 ± 8.9(23.4–30.9) 22.0 ± 8.0(18.3–25.6) 26.0 ± 9.3(21.8–30.2)
Time constant (tau1; sec) 32.2 ± 10.5(27.7–36.6) 26.9 ± 9.1(23.0–30.7) 36.0 ± 12.0(30.5–41.4) 25.9 ± 9.6(21.5–30.3)#
Mean response time (sec) 55.4 ± 6.7(52.5–58.2) 54.0 ± 7.8(50.7–57.3) 58.0 ± 10.1(53.4–62.6) 51.9 ± 5.5(49.3–54.4)#
Cardiac output on-kinetic response
Baseline Q̇ (L.min–1) 11.0 ± 1.3(10.5–11.6) 12.3 ± 2.2(11.3–13.3)* 11.0 ± 1.6(10.2–11.7) 10.6 ± 1.0(10.2–11.1)
Amplitude (L.min–1) 6.6 ± 3.1(5.2–8.0) 5.2 ± 1.2(4.7–5.8)# 4.0 ± 1.3(3.4–4.6) 4.6 ± 1.5(3.9–5.3)#
Time delay (sec) 7.2 ± 4.0(5.4–9.0) 4.0 ± 1.9(3.1–4.8) 6.0 ± 2.5(4.8–7.2) 8.8 ± 6.6(5.7–11.9)#
Time constant (tau1; sec) 45.9 ± 10.1(41.5–50.4) 38.4 ± 10.2(33.9–42.9)* 50.0 ± 10.7(45.0–55.0) 34.4 ± 10.5(29.5–39.3)*
Mean response time (sec) 53.1 ± 10.3(48.6–57.7) 38.1 ± 10.2(33.9–42.9)* 56.0 ± 10.0(51.3–60.7) 43.2 ± 5.8(40.5–45.9)*
Heart rate on-kinetic response
Baseline heart rate (bpm) 105.5 ± 10.4(100.9–110.1) 106.8 ± 12.1(101.4–112.2) 113.8 ± 8.5(109.8–117.8) 107.3 ± 10.0(102.6–112.0)*
Amplitude (bpm) 44.6 ± 8.1(41.0–48.3) 39.9 ± 9.4(35.8–44.1)* 32.0 ± 6.8(28.8–35.2) 32.4 ± 9.6(27.9–36.9)
Time delay (sec) 4.1 ± 4.3(2.2–6.0) 2.8 ± 1.0(2.3–3.2) 4.2 ± 3.4(2.7–5.8) 6.9 ± 6.3(4.0–9.9)
Time constant (tau1; sec) 55.4 ± 9.1(51.4–59.5) 43.2 ± 9.9(38.8–47.6)* 53.1 ± 15.7(45.8–60.5) 38.2 ± 10.4(33.4–43.1)*
Mean response time (sec) 59.5 ± 10.4(54.9–64.1) 46.0 ± 10.2(41.5–50.5)* 57.3 ± 14.0(50.8–63.9) 45.1 ± 6.1(42.3–48.0)*
Heavy intensity exercise (Δ50% POVT)
Pulmonary oxygen uptake on-kinetic response
Baseline V˙O2 (L.min–1) 0.9 ± 0.2(0.8–0.9) 0.8 ± 0.1(0.7–0.8) 0.8 ± 0.2(0.8–0.9) 0.8 ± 0.1(0.8–0.9)
Amplitude (L.min–1) 1.3 ± 0.3(1.2–1.4) 1.3 ± 0.3(1.2–1.4) 1.3 ± 0.4(1.1–1.5) 1.3 ± 0.4(1.1–1.4)
Time delay (sec) 18.2 ± 7.2(15.2–21.3) 22.7 ± 5.9(20.2–25.3)# 20.6 ± 5.9(17.9–23.3) 23.1 ± 5.4(20.7–25.6)
Time constant (tau1; sec) 32.5 ± 10.0(28.3–36.7) 28.5 ± 10.0(24.3–32.7) 33.5 ± 8.5(29.7–37.4) 26.6 ± 9.0(22.5–30.7)#
Mean response time (sec) 50.7 ± 7.6(47.5–53.9) 51.2 ± 8.0(47.9–54.6) 54.1 ± 8.9(50.1–58.2) 49.7 ± 6.9(46.6–52.9)#
Cardiac output on-kinetic response
Baseline Q̇ (L.min–1) 13.0 ± 1.7(12.3–13.8) 14.2 ± 1.7(13.5–15.0)# 12.7 ± 1.0(12.3–13.1) 12.2 ± 1.0(11.8–12.7)#
Amplitude (L.min–1) 6.7 ± 2.0(5.8–7.6) 6.7 ± 2.1(5.7–7.6) 5.4 ± 1.6(4.7–6.2) 6.1 ± 1.5(5.4–6.8)#
Time delay (sec) 5.8 ± 4.4(3.8–7.7) 9.7 ± 3.8(8.0–11.4)# 8.7 ± 3.7(7.0–10.4) 9.9 ± 5.0(7.6–12.2)
Time constant (tau1; sec) 45.4 ± 9.4(41.3–49.6) 31.8 ± 6.6(28.9–34.8)* 44.3 ± 13.8(37.9–50.8) 30.0 ± 10.1(25.2–34.7)*
Mean response time (sec) 51.2 ± 9.9(46.8–56.0) 41.5 ± 6.8(38.5–44.5)* 53.0 ± 13.8(46.6–59.5) 39.9 ± 6.3(36.9–42.8)*
Heart rate
Baseline heart rate (bpm) 122.8 ± 7.8(119.4–126.3) 120.3 ± 10.1(115.8–124.8) 125.9 ± 10.5(121.0–130.8) 120.7 ± 11.0(115.6–125.8)#
Amplitude (bpm) 56.2 ± 8.2(52.6–59.9) 52.9 ± 8.0(49.4–56.5)* 49.0 ± 10.9(44.0–54.1) 48.1 ± 10.5(43.2–53.0)
Time delay (sec) 7.4 ± 5.3(5.0–9.7) 5.1 ± 3.8(3.4–6.8)# 10.3 ± 3.1(8.8–11.7) 6.8 ± 1.9(5.9–7.7)*
Time constant (tau1; sec) 45.2 ± 6.1(42.5–47.9) 43.6 ± 10.9(38.8–48.4) 43.4 ± 9.8(38.8–48.0) 41.3 ± 4.9(39.0–43.6)
Mean response time (sec) 52.6 ± 10.1(48.1–57.0) 48.7 ± 12.2(43.3–54.1)* 53.7 ± 9.9(49.1–58.3) 48.1 ± 5.0(45.8–50.4)*
POVT, power output at ventilatory threshold; sec, seconds; V˙O2, volume of oxygen consumed. Q̇, cardiac output; sec, seconds. Parametric data are presented as mean ± SD (95%CI); non-parametric data are presented as median [IQR] (95% CI). #p < 0.05 vs. pre-training; *p < 0.001 vs. pre-training.

Pulmonary oxygen uptake profile (with model best fit line and residuals) for a representative naturally menstruating participant (A) and oral contraceptive using participant (B) during the transition to moderate (left) and heavy (right) intensity exercise, pre-training (open circles and gray line), post-training (closed circles and black line). Mean response time (MRT) is indicated by a matched vertical line. Shading represents exercise onset.

 Cardiac Output Kinetic Response

There were no differences between groups for cardiac output kinetic response represented as mean response time (Q̇MRT) pre- or post-training at either intensity (p > 0.05). Following training, Q̇MRT improved in both groups at both intensities (p < 0.001). There was a significant group x time interaction for Q̇MRT at moderate (p = 0.003) and high (p = 0.039) intensity. Post hoc analyses suggested there was no difference between groups for pre- to post-training improvements in Q̇MRT at moderate [OC-group −15.0 s vs. MC-group −12.8 s; p = 0.405], and high intensities [OC-group −9.7 s vs. MC-group −13.1 s; p = 0.282]. Data are presented in Table 2 and the cardiac output profile at both intensities from a representative participant is presented in Supplementary Figure S1. Standardized between-group differences for within-group changes (Cohen’s d) demonstrated that there was no clear between-group difference in Q̇MRT adaptation to training at both moderate (0.19 ± 0.48) and heavy (0.29 ± 0.53) exercise intensities (Table 3).

Standardized between-group differences for within-group changes for the oral contraceptive group versus the menstrual cycle group following training and de-training.

Cohen’s d ± SD 95% CI % chance for OC change to be higher/trivial/lower than MC Descriptive of difference
Time to fatigue (sec) 0.96 ± 1.04 −0.08–2.01 93/6/1 Likely higher
Moderate intensity (90% POVT)
Mean response time (sec) 0.54 ± 0.67 −0.13–1.21 85/14/2 Likely higher
Tissue saturation index (slope) −1.03 ± 0.71 −1.74–−0.33 0/1/99 V. likely lower
Mean response time (Q̇) 0.19 ± 0.48 −0.29–0.68 49/46/5 Unclear
Mean response time (HR) −0.10 ± 0.46 −0.56–0.36 10/57/33 Unclear
Heavy intensity (Δ50% POVT)
Mean response time (sec) 0.57 ± 0.49 0.08–1.07 93/6/0 V. likely higher
Tissue saturation index (slope) −0.66 ± 0.41 −1.07–−0.25 0/1/99 V. likely lower
Mean response time (Q̇) 0.29 ± 0.53 −0.25–0.82 63/34/4 Unclear
Mean response time (HR) 0.17 ± 0.50 −0.32–0.67 45/48/7 Unclear
SD, standard deviation; CI, confidence interval; OC, oral contraceptive group; MC, naturally menstruating group; sec, seconds; Q, cardiac output; HR, heart rate; POVT, power output at ventilatory threshold; V. likely, very likely. Heart Rate Kinetic Response

There was no difference between groups for HRMRT at any time point at either intensity (p > 0.05). Following training, HRMRT improved in both groups and at each intensity (p < 0.05). There was no significant group × time interaction for HRMRT at either moderate (p = 0.269) or high intensity (p = 0.712). Data are presented in Table 2 and the heart rate profile at both intensities from a representative participant is presented in Supplementary Figure S2. Standardized between-group differences for within-group changes (Cohen’s d) demonstrated no clear differences between groups (Table 3).

 Muscle Deoxygenation Rates

There was no difference between groups for TSI (presented here as an index/percentage) at baseline for either intensity (moderate: MC −0.18 ± 0.08; OC −0.16 ± 0.05; high: MC −0.24 ± 0.07; OC −0.26 ± 0.12; all p > 0.05). Post-training, the OC-group had improved TSI compared to the MC-group at both moderate (OC −0.24 ± 0.08; MC −0.19 ± 0.07; p = 0.027) and high (OC −0.34 ± 0.14; MC −0.25 ± 0.07; p = 0.018) intensities. Following training, TSI was improved in the OC-group at both intensities [moderate: Δ−0.08 (50%); high: Δ−0.08 (30.8%); both p < 0.001]; but did not change in the MC-group at either intensity [moderate: Δ−0.01 (6.0%), p = 0.295; high: Δ−0.01 (4.1%) p = 0.422]. There was a significant group × time interaction for TSI at moderate (p = 0.001), and high (p < 0.001) intensities. Post hoc analyses suggested the OC-group showed greater improvement in TSI at both moderate (p = 0.004) and high (p = 0.002) intensities following training compared to the MC-group. Data are presented in Figure 4. Standardized between-group differences for within-group changes (Cohen’s d) demonstrated that the OC-group had a likely greater TSI adaptation to training compared to the MC-group at moderate (−1.03 ± 0.71) and heavy (−0.66 ± 0.41) intensities (Table 3).

Tissue saturation index at a moderate exercise intensity (A) and a high exercise intensity (B) for the oral contraceptive group (n = 25) and menstrual cycle group (n = 22) pre- and post-training. Data are mean ± SD. *p < 0.001 vs. pre- training; #p < 0.05 group × time interaction.

 Discussion

The present study investigated the influence of OC use on time to fatigue and changes in pulmonary, cardiovascular and muscle oxygen uptake kinetics to sprint interval training (SIT). Pulmonary oxygen uptake kinetics (τV˙O2p) as determined by MRT, improved following training in the MC-group only. However, tissue saturation index (TSI) improved in the OC-group only. Despite these differences, improvements in time to fatigue in response to SIT did not differ between groups.

Significant improvements (30%) in the mean response time for τV˙O2p at both moderate and high intensity exercise were observed with training in the MC-group only; the OC-group showed no significant change from baseline. When the MC-group was sub-grouped for potential LPD, the MC LPD sub-group also demonstrated dampened τV˙O2p. The improvement in τV˙O2p observed in the MC-group is consistent with previous research, where Murias et al. (2011) demonstrated that τV˙O2p was improved in response to endurance training in young women. However, these authors did not include women using an OC and our results suggest that exogenous hormones may dampen the τV˙O2p adaptations to training. It is also possible that the absence of endogenous ovarian hormones (within the probable LPD group) may also dampen τV˙O2p adaptations. This is consistent with findings from our previous study (Schaumberg et al., 2017a), where V˙O2peak improvement was lower in OC users compared to naturally menstruating women following training.

Improvement in τV˙O2p can increase exercise capacity at both moderate and higher intensities of exercise (Poole and Richardson, 1997). Oxygen delivery to skeletal muscle is challenged during high intensity exercise training, thereby eliciting structural adaptations. This leads to an increased blood flow and an increase in pulmonary and/or muscle oxygen kinetic responses at the onset of constant-load exercise (Poole and Richardson, 1997; Jones and Poole, 2005). Therefore, improvements in τV˙O2p may be limited by either central (e.g., cardiac output) or peripheral (e.g., skeletal muscle blood flow) adaptations (Barstow et al., 1990; Grassi, 2006). OC use has been shown to decrease collagen synthesis, and therefore capillarization, in response to exercise (Hansen et al., 2011). Therefore, reduced capillarization adaptation with exercise may dampen further central adaptations such as τV˙O2p. In addition, OC use increases the release of hormones (e.g., growth hormone and other glucoregulatory hormones) (62), which can influence carbohydrate metabolism during exercise (Davidson and Holzman, 1973; Bunt, 1990; Bonen et al., 1991). As τV˙O2p is controlled intracellularly via oxidative phosphorylation, glucose is required to fuel the turnover of ATP for energy production (McKay et al., 2009). If there is limited blood glucose available, this could limit the improvement of τV˙O2p. Luteal phase deficiency is a menstrual disturbance that may reflect early stages of low energy availability in women (De Souza et al., 2003). Therefore, the dampened τV˙O2p adaptation seen within the probable LPD group further support this potential mechanism.

Cardiac output on-kinetics improved in both groups at both moderate and high intensity loads following training. It is possible that this measure was not sensitive enough to detect any between-group differences. However, it is important to note that the adaptation of the women in this study was comparable to men, where it has previously been demonstrated that the heart rate time constant to moderate intensity exercise improves with just eight sessions of HIIT, and that end heart rate is also improved (McKay et al., 2009). Recent work by Howden et al. (2015) demonstrated that women had a markedly blunted cardiovascular response to 1 year of endurance training, compared to males. Therefore, further investigation into cardiovascular on-kinetic responses is warranted to explain the lack of adaptation in the present study.

The present study investigated TSI as an overall index of the muscle deoxygenation kinetics, representative of the dynamic balance between oxygen consumption of the muscle tissue and supply. While most literature reports total oxy- or deoxy-hemoglobin, TSI arguably provides a more effective measure of the overall muscular response (because of its incorporation of both oxy- and deoxy-hemoglobin and consumption and supply of oxygen), and is closely related to [HHB] (Ihsan et al., 2013).

The first, and most obvious, conclusion for the significant TSI improvement in the OC-group and no change in the MC-group is that the OC-group showed significant peripheral adaptation with training compared to the MC-group, yet, when analyzed in conjunction with the τV˙O2p adaptations, this assumption must be considered with caution, as TSI is a rate relative to both supply and utilization of oxygen at the muscle. TSI and τV˙O2p are often analyzed and reported together (Gurd et al., 2007; McKay et al., 2009; Murias et al., 2011), as TSI is indicative of peripheral adaptation whilst τV˙O2p is indicative of central adaptations (Murias et al., 2011). An increase in TSI (also demonstrated with the increase relative to exercise intensity) may indicate a mismatch between oxygen delivery and utilization during exercise (Spencer et al., 2013). As TSI adjusts more rapidly in response to exercise compared to τV˙O2p there is insufficient oxygen delivery for the working muscles and temporal dissociation occurs, which may result in greater reliance on anaerobic energy production at the working muscles (MacPhee et al., 2005; Spencer et al., 2013), and be potentially detrimental to exercise capacity.

The concomitant changes in TSI and τV˙O2p in the MC- and OC-groups can be further considered through the concept of a physiological phenomenon known as the ‘transient overshoot’ (Murias et al., 2011), whereby training instigates more efficient oxygen utilization as well as blood distribution, thereby causing a speeding of τV˙O2p with no change in TSI. An increase in muscle oxygen extraction causes insufficient oxygen delivery, instigating a temporal dissociation period between [HHB] and τV˙O2p (MacPhee et al., 2005; Spencer et al., 2013). This overshoot may be present with OC use, with greater utilization of oxygen at the muscle occurring after training due to peripheral adaptations, with concurrent blunting of oxygen supply indicated by τV˙O2p.

Dynamic changes in NIRS-derived muscle deoxygenation (represented by TSI) provide insights into the balance between local muscle oxygen availability and utilization during exercise (Wolf et al., 2007; Spencer et al., 2013). During exercise muscle deoxygenation adjusts to an increased workload more rapidly than τV˙O2p, resulting in a transient period characterized by an increased relative reliance on oxygen extraction to support a given metabolic rate. This temporary dissociation between adjustments of [HHB] and τV˙O2p suggest transient oxygen delivery insufficiency (Murias et al., 2013; Spencer et al., 2013) for the rate of oxygen consumption.

Therefore, we observed an abnormal response of τV˙O2p and TSI to exercise training in OC users, demonstrated by an apparent acceleration of TSI, which may suggest a mismatch of oxygen delivery and utilization at the exercising muscles. In our participant sample, oxygen delivery did not improve with training in OC users compared to naturally menstruating women, despite the apparent increase in local muscular oxygen utilization. Therefore, the next consideration based on these results is that exogenous hormones do not appear to influence peripheral adaptation to training and may indeed be beneficial. Indeed, the mechanism for the change in TSI in OC users warrants further investigation. Previous research has speculated the mismatch between oxygen utilization and delivery with OC use is due to its direct effect on the female sex hormones (estrogen and progesterone, and their exogenous forms). However, several other hormones (e.g., growth hormone, inflammatory factors, and free androgens) are also influenced by OC use and should therefore be considered. For example, growth factors, including growth hormone and free androgens are significantly different from naturally menstruating women; therefore it is necessary to further investigate this possible factor, due to their influence on capillarization and subsequent oxygen utilization at the muscle (Hansen et al., 2011; Hansen and Kjaer, 2014). Further, the influence of OC use on mitochondrial/oxidative enzymes is unknown, and further investigation into the potential influence of OC use on mitochondrial oxidative capacity is warranted.

 Conclusion

Exogenous ovarian hormones found in the oral contraceptive pill may be, at least in part, responsible for the dampened physiological adaptations to training in OC users. Although both OC users and naturally menstruating women improved TTF, we observed a dampened response of central physiological adaptation, demonstrated by pulmonary oxygen uptake kinetics in the OC group. This may have been offset by the greater improvements in muscle oxygen utilization in OC users, compared to the MC group. These results provide insight into potential mechanisms related to training adaptation in women. Based on these results, potential mechanisms may include the lack of endogenous ovarian hormones, as well as the influence of exogenous hormones on the overall endocrinological profile, including growth hormone and free androgens, may be implicated in dampening the physiological adaptations to training with OC use. These potential mechanisms warrant investigation to further elucidate the influence of OC use on adaptations to training in women.

 Data Availability Statement

The datasets generated for this study are available on request to the corresponding author.

 Ethics Statement

This study involving human participants was reviewed and approved by The University of Queensland Human Research Ethics Committee. The patients/participants provided their written informed consent to participate in this study.

 Author Contributions

MS, JS, DJ, XJ, LE, and TS designed the protocol and methods. MS and EH collected the data. MS, JS, and EH analyzed the data. All authors contributed to the writing and drafting of the manuscript and read and approved the manuscript.

 Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The authors wish to thank Mr. Gary Wilson for his biochemical analysis guidance and technical support, and the dedicated international and Honors research students who assisted with the data collection.

 Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2020.00629/full#supplementary-material

Cardiac output profile (with model best fit line and residuals; time in seconds on the x-axis) for a representative naturally menstruating participant (A) and oral contraceptive using participant (B) during the transition to moderate (left) and heavy (right) intensity exercise, pre-training (open circles and gray line) and post-training (closed circles and black line). Mean response time (MRT) is indicated by a matched vertical line. Shading represents exercise onset.

Heart rate profile (with model best fit line and residuals) for a representative naturally menstruating participant (A) and oral contraceptive using participant (B) during the transition to moderate (left) and heavy (right) intensity exercise, pre-training (open circles and gray line) and post-training (closed circles and black line). Mean response time (MRT) is indicated by a matched vertical line. Shading represents exercise onset.

Pulmonary oxygen uptake, cardiac output and heart rate on-kinetic responses to moderate and heavy intensity exercise, pre- and post-training in naturally-menstruating women, sub-grouped for normal luteal phase characteristics (n = 16) and apparent luteal phase deficient women (n = 8); as determined by serum progesterone concentrations.

 References Astorino T. A. Allen R. P. Robertson D. W. Jurancich M. Lewis R. McCarthy K. (2011). Adaptations to high-intensity training are independent of gender. Eur. J. Appl. Physiol. 111 1279–1286. 10.1007/s00421-010-1741-y 21132441 Bae S. Y. Hamaoka T. Katsumura T. Shiga T. Ohno H. Haga S. (2000). Comparison of muscle oxygen consumption measured by near infrared continuous wave spectroscopy during supramaximal and intermittent pedalling exercise. Int. J. Sports Med. 21 168–174. 10.1055/s-2000-8880 10834347 Barstow T. J. Lamarra N. Whipp B. J. (1990). Modulation of muscle and pulmonary O2 uptakes by circulatory dynamics during exercise. J. Appl. Physiol. 68 979–989. 10.1152/jappl.1990.68.3.979 2341363 Batterham A. M. Hopkins W. G. (2005). Making meaningful inferences about magnitudes. Sport Sci. 9 6–13. Batterham A. M. Hopkins W. G. (2015). The case for magnitude-based inference. Med. Sci. Sports Exerc. 47:885. 10.1249/mss.0000000000000551 25783666 Berenson A. B. Rahman M. (2009). Changes in weight, total fat, percent body fat, and central-to-peripheral fat ratio associated with injectable and oral contraceptive use. Am. J. Obstet. Gynecol. 200 e321–e329. Berger N. J. Tolfrey K. Williams A. G. Jones A. M. (2006). Influence of continuous and interval training on oxygen uptake on-kinetics. Med. Sci. Sports Exerc. 38 504–512. 10.1249/01.mss.0000191418.37709.81 Bonen A. Haynes F. W. Graham T. E. (1991). Substrate and hormonal responses to exercise in women using oral contraceptives. J. Appl. Physiol. 70 1917–1927. 10.1152/jappl.1991.70.5.1917 1864771 Borg G. A. (1973). Perceived exertion: a note on “history” and methods. Med. Sci. Sports 5 90–93. Bunt J. C. (1990). Metabolic actions of estradiol: significance for acute and chronic exercise responses. Med. Sci. Sports Exerc. 22 286–290. Carter S. L. Rennie C. D. Hamilton S. J. Tarnopolsky H. (2001). Changes in skeletal muscle in males and females following endurance training. Can. J. Physiol. Pharmacol. 79 386–392. 10.1139/y01-008 Casazza G. A. Suh S. H. Miller B. F. Navazio F. M. Brooks G. A. (2002). Effects of oral contraceptives on peak exercise capacity. J. Appl. Physiol. 93 1698–1702. 10.1152/japplphysiol.00622.2002 12381756 Charkoudian N. Johnson J. M. (1997). Modification of active cutaneous vasodilation by oral contraceptive hormones. J. Appl. Physiol. 83 2012–2018. 10.1152/jappl.1997.83.6.2012 9390975 Charkoudian N. Joyner M. J. (2004). Physiologic considerations for exercise performance in women. Clin. Chest Med. 25 247–255. 10.1016/j.ccm.2004.01.001 15099886 Charloux A. Lonsdorfer-Wolf E. Richard R. Lampert E. Oswald-Mammosser M. Mettauer B. (2000). A new impedance cardiograph device for the non-invasive evaluation of cardiac output at rest and during exercise: comparison with the “direct”. Fick method. Eur. J. Appl. Physiol. 82 313–320. 10.1007/s004210000226 10958374 Cohen J. (1988). Statistical Power Analysis For The Behavioural Sciences. Hillsdale, NJ: Lawrence Erlbaum Associates. Coney P. J. Washenik K. Langley R. G. B. DiGiovanna J. J. Harrison D. D. (2001). Weight change and adverse event incidence with a low-dose oral contraceptive: two randomized, placebo-controlled trials. Contraception 63 297–302. 10.1016/s0010-7824(01)00208-6 Davidson M. B. Holzman G. B. (1973). Role of growth hormone in the alteration of carbohydrate metabolism induced by oral contraceptive agents. J. Clin. Endocrinol. Metab. 36 246–255. 10.1210/jcem-36-2-246 4683182 De Souza M. J. Van Heest J. Demers L. M. Lasley B. L. (2003). Luteal phase deficiency in recreational runners: evidence for a hypometabolic state. J. Clin. Endocrinol. Metab. 88 337–346. 10.1210/jc.2002-020958 12519874 DeLorey D. S. Kowalchuk J. M. Heenan A. P. Dumanoir G. R. Paterson D. H. (2007). Prior exercise speeds pulmonary O2 uptake kinetics by increases in both local muscle O2 availability and O2 utilization. J. Appl. Physiol. 103 771–778. 10.1152/japplphysiol.01061.2006 17495116 Dorado C. Sanchis-Moysi J. Calbet J. A. (2004). Effects of recovery mode on performance, O2 uptake, and O2 deficit during high-intensity intermittent exercise. Can. J. Appl. Physiol. 29 227–244. 10.1139/h04-016 15199225 Gibala M. J. Little J. P. MacDonald M. J. Hawley J. A. (2012). Physiological adaptations to low-volume, high-intensity interval training in health and disease. J. Physiol. 590 1077–1084. 10.1113/jphysiol.2011.224725 22289907 Gibala M. J. Little J. P. van Essen M. Wilkin G. P. Burgomaster K. A. Safdar A. (2006). Short-term sprint interval versus traditional endurance training: similar initial adaptations in human skeletal muscle and exercise performance. J. Physiol. 575 901–911. 10.1113/jphysiol.2006.112094 16825308 Gosselin L. E. Kozlowski K. F. DeVinney-Boymel L. Hambridge C. (2012). Metabolic response of different high-intensity aerobic interval exercise protocols. J. Strength Cond. Res. 26 2866–2871. 10.1519/jsc.0b013e318241e13d 22124355 Grassi B. (2001). Regulation of oxygen consumption at exercise onset: is it really controversial? Exerc. Sport Sci. Rev. 29 134–138. 10.1097/00003677-200107000-00009 11474962 Grassi B. (2006). Oxygen uptake kinetics: why are they so slow? And what do they tell us? J. Physiol. Pharmacol. 57(Suppl. 10), 53–65. Greer J. B. Modugno F. Allen G. O. Ness R. B. (2005). Androgenic progestins in oral contraceptives and the risk of epithelial ovarian cancer. Obstet. Gynecol. 105 731–740. 10.1097/01.aog.0000154152.12088.48 Gurd B. J. Scheid J. Paterson D. H. Kowalchuk J. M. (2007). O2 uptake and muscle deoxygenation kinetics during the transition to moderate-intensity exercise in different phases of the menstrual cycle in young adult females. Eur. J. Appl. Physiol. 101 321–330. 10.1007/s00421-007-0505-9 17618450 Hansen M. Kjaer M. (2014). Influence of sex and estrogen on musculotendinous protein turnover at rest and after exercise. Exerc. Sport Sci. Rev. 42 183–192. 10.1249/jes.0000000000000026 25062001 Hansen M. Langberg H. Holm L. Miller B. F. Petersen S. G. Doessing S. (2011). Effect of administration of oral contraceptives on the synthesis and breakdown of myofibrillar proteins in young women. Scand. J. Med. Sci. Sports 21 62–72. 10.1111/j.1600-0838.2009.01002.x 19883384 Hopkins W. G. Marshall S. W. Batterham A. M. Hanin J. (2009). Progressive statistics for studies in sports medicine and exercise science. Med. Sci. Sports Exerc. 41 3–13. 10.1249/MSS.0b013e31818cb278 19092709 Howden E. J. Perhonen M. Peshock R. M. Zhang R. Arbab-Zadeh A. Adams-Huet B. (2015). Females have a blunted cardiovascular response to one year of intensive supervised endurance training. J. Appl. Physiol. 119 37–46. 10.1152/japplphysiol.00092.2015 25930024 Hsu A. R. Barnholt K. E. Grundmann N. K. Lin J. H. McCallum S. W. Friedlander A. L. (2006). Sildenafil improves cardiac output and exercise performance during acute hypoxia, but not normoxia. Science 100 2031–2040. 10.1152/japplphysiol.00806.2005 16455814 Hughson R. L. Tschakovsky M. E. Houston M. E. (2001). Regulation of oxygen consumption at the onset of exercise. Exerc. Sport Sci. Rev. 29 129–133. 10.1097/00003677-200107000-00008 11474961 Ihsan M. Watson G. Lipski M. Abbiss C. R. (2013). Influence of postexercise cooling on muscle oxygenation and blood volume changes. Med. Sci. Sports Exerc. 45 876–882. 10.1249/MSS.0b013e31827e13a2 23247707 Isacco L. Duche P. Boisseau N. (2012). Influence of hormonal status on substrate utilisation at rest and during exercise in the female population. Sports Med. 42 327–342. 10.2165/11598900-000000000-00000 22380007 Jones A. M. Burnley M. (2009). Oxygen uptake kinetics: an underappreciated determinant of exercise performance. Int. J. Sports Physiol. Perform. 4 524–532. 10.1123/ijspp.4.4.524 20029103 Jones A. M. Poole D. C. (2005). Oxygen Uptake Kinetics In Sport, Exercise And Medicine. New York, NY: Routledge. Joyce S. Sabapathy S. Bulmer A. Minahan C. (2013). Effect of long-term oral contraceptive use on determinants of endurance performance. J. Strength Cond. Res. 27 1891–1896. 10.1519/JSC.0b013e3182736935 22996028 Lebrun C. M. Petit M. A. McKenzie D. C. Taunton J. E. Prior J. C. (2003). Decreased maximal aerobic capacity with use of a triphasic oral contraceptive in highly active women: a randomised controlled trial. Br. J. Sports Med. 37 315–320. 10.1136/bjsm.37.4.315 12893716 Lepretre P.-M. Koralsztein J.-P. Billat V. L. (2004). Effect of exercise intensity on relationship between V˙O2max and cardiac output. Cell 36 1357–1363. 10.1249/01.Mss.0000135977.12456.8f MacPhee S. L. Shoemaker J. K. Paterson D. H. Kowalchuk J. M. (2005). Kinetics of O2 uptake, leg blood flow, and muscle deoxygenation are slowed in the upper compared with lower region of the moderate-intensity exercise domain. J. Appl. Physiol. 99 1822–1834. 10.1152/japplphysiol.01183.2004 16037398 McKay B. R. Paterson D. H. Kowalchuk J. M. (2009). Effect of short-term high-intensity interval training vs. continuous training on O2 uptake kinetics, muscle deoxygenation and exercise performance. J. Appl. Physiol. 107 128–138. 10.1152/japplphysiol.90828.2008 19443744 Murias J. M. Kowalchuk J. M. Paterson D. H. (2011). Speeding of VO2 kinetics in response to endurance-training in older and young women. Eur. J. Appl. Physiol. 111 235–243. 10.1007/s00421-010-1649-6 20857137 Murias J. M. Spencer M. D. Keir D. A. Paterson D. H. (2013). Systemic and vastus lateralis muscle blood flow and O2 extraction during ramp incremental cycle exercise. Am. J. Physiol. Regul. Integr. Comp. Physiol. 304 R720–R725. 10.1152/ajpregu.00016.2013 23515617 Notelovitz M. (1987). The effect of low-dose oral contraceptives on cardiorespiratory function, coagulation, and lipids in exercising young women: a preliminary report. Am. J. Obstet. Gynecol. 156 591–598. 10.1016/0002-9378(87)90059-7 Overend T. J. Paterson D. H. Cunningham D. A. (1992). The effect of interval and continuous training on the aerobic parameters. Can. J. Sport Sci. 17 129–134. Ozyener F. Rossiter H. B. Ward S. A. Whipp B. J. (2001). Influence of exercise intensity on the on- and off-transient kinetics of pulmonary oxygen uptake in humans. J. Physiol. 533(Pt 3), 891–902. 10.1111/j.1469-7793.2001.t01-1-00891.x 11410644 Poole D. C. Richardson R. S. (1997). Determinants of oxygen uptake. Implications for exercise testing. Sports Med. 24 308–320. 10.2165/00007256-199724050-00003 9368277 Ready A. E. Eynon R. B. Cunningham D. A. (1981). Effect of interval training and detraining on anaerobic fitness in women. Can. J. Appl. Sport Sci. 6 114–118. Richard R. Lonsdorfer-Wolf E. Charloux A. Doutreleau S. Buchheit M. Oswald-Mammosser M. (2001). Non-invasive cardiac output evaluation during a maximal progressive exercise test, using a new impedance cardiograph device. Eur. J. Appl. Physiol. 85 202–207. 10.1007/s004210100458 11560071 Rossiter H. B. Kowalchuk J. M. Whipp B. J. (2006). A test to establish maximum O2 uptake despite no plateau in the O2 uptake response to ramp incremental exercise. J. Appl. Physiol. 100 764–770. 10.1152/japplphysiol.00932.2005 16282428 Schaumberg M. A. Jenkins D. G. Janse D. E. Emmerton L. M. Skinner T. L. (2017a). Oral contraceptive use dampens physiological adaptations to sprint interval training. Med. Sci. Sports Exerc. 49 717–727. 10.1249/mss.0000000000001171 27898641 Schaumberg M. A. Jenkins D. G. Jonge X. A. Emmerton L. M. Skinner T. L. (2017b). Three-step method for menstrual and oral contraceptive cycle verification. J. Sci. Med. Sport. 20 965–969. 10.1016/j.jsams.2016.08.013 28684053 Shiraev T. Barclay G. (2012). Evidence based exercise - clinical benefits of high intensity interval training. Aust. Fam. Phys. 41 960–962. Shoemaker J. K. Phillips S. M. Green H. J. Hughson R. L. (1996). Faster femoral artery blood velocity kinetics at the onset of exercise following short-term training. Cardiovasc. Res. 31 278–286. 10.1016/s0008-6363(95)00199-9 Spencer M. D. Murias J. M. Kowalchuk J. M. Paterson D. H. (2013). Effect of moderate-intensity work rate increment on phase II tauVO(2), functional gain and Delta[HHb]. Eur. J. Appl. Physiol. 113 545–557. 10.1007/s00421-012-2460-3 22829340 Stachenfeld N. S. Silva C. Keefe D. L. (2000). Estrogen modifies the temperature effects of progesterone. J. Appl. Physiol. 88 1643–1649. 10.1152/jappl.2000.88.5.1643 10797125 Stanley J. Peake J. M. Buchheit M. (2013). Consecutive days of cold water immersion: effects on cycling performance and heart rate variability. Eur. J. Appl. Physiol. 113 371–384. 10.1007/s00421-012-2445-2 22752345 Stanley J. Peake J. M. Coombes J. S. Buchheit M. (2014). Central and peripheral adjustments during high-intensity exercise following cold water immersion. Eur. J. Appl. Physiol. 114 147–163. 10.1007/s00421-013-2755-z 24158407 Stathokostas L. Kowalchuk J. M. Petrella R. J. Paterson D. H. (2009). Moderate and heavy oxygen uptake kinetics in postmenopausal women. Appl. Physiol. Nutr. Metab. 34 1065–1072. 10.1139/h09-107 20029515 Talanian J. L. Galloway S. D. R. Heigenhauser G. J. F. Bonen A. Spriet L. L. (2007). Two weeks of high-intensity aerobic interval training increases the capacity for fat oxidation during exercise in women. J. Appl. Physiol. 102 1439–1447. 10.1152/japplphysiol.01098.2006 17170203 Tschakovsky M. E. Hughson R. L. (1999). Interaction of factors determining oxygen uptake at the onset of exercise. J. Appl. Physiol. 86 1101–1113. 10.1152/jappl.1999.86.4.1101 10194190 Vanheest J. L. Mahoney C. E. Rodgers C. D. (2005). “Oral contraceptive use and physical performance,” in The Endocrine Sysytem in Sports and Exercise, eds Kraemer W. J. Rogol A. D. (Carlton: Blackwell Publishing), 250–260. 10.1002/9780470757826.ch19 Vincent J. (1999). Statistics in Kinesiology. Champaign, IL: Human Kinetics. Whipp B. J. Ward S. A. Rossiter H. B. (2005). Pulmonary O2 uptake during exercise: conflating muscular and cardiovascular responses. Med. Sci. Sports Exerc. 37 1574–1585. 10.1249/01.mss.0000177476.63356.22 Wolf M. Ferrari M. Quaresima V. (2007). Progress of near-infrared spectroscopy and topography for brain and muscle clinical applications. J. Biomed. Opt. 12:e062104. 10.1117/1.2804899 Xu F. Rhodes E. C. (1999). Oxygen uptake kinetics during exercise. Sports Med. 27 313–327. 10.2165/00007256-199927050-00003 10368878 Zierke J. M. (2010). The Effects Of High Intensity Interval Training On Sub-Maximal Exercise Responses and VO2max in Active Women Taking Oral Contraceptives. Masters thesis, California State University, Long Beach, CA.