Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Daily fatigue-recovery balance monitoring with heart rate variability in well-trained female cyclists on the Tour de France circuit

  • Anna Barrero ,

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    annabarrero@gmail.com

    Affiliations M2S Laboratory, University of Rennes 2, Rennes, France, CIC-CIT INSERM UMR 1099, Rennes, France

  • Frédéric Schnell,

    Roles Conceptualization, Supervision

    Affiliation University of Rennes 1, Department of Sports Medicine University Hospital of Rennes, INSERM, LTSI-UMR1099, Rennes, France

  • Guy Carrault,

    Roles Conceptualization, Formal analysis

    Affiliations CIC-CIT INSERM UMR 1099, Rennes, France, University of Rennes 1, Department of Sports Medicine University Hospital of Rennes, INSERM, LTSI-UMR1099, Rennes, France

  • Gaelle Kervio,

    Roles Conceptualization

    Affiliation CIC-CIT INSERM UMR 1099, Rennes, France

  • David Matelot,

    Roles Conceptualization

    Affiliation CIC-CIT INSERM UMR 1099, Rennes, France

  • François Carré,

    Roles Conceptualization, Funding acquisition, Methodology, Supervision, Validation, Writing – review & editing

    Affiliation University of Rennes 1, Department of Sports Medicine University Hospital of Rennes, INSERM, LTSI-UMR1099, Rennes, France

  • Solène Le Douairon Lahaye

    Roles Conceptualization, Data curation, Funding acquisition, Methodology, Project administration, Resources, Supervision, Validation

    Affiliations M2S Laboratory, University of Rennes 2, Rennes, France, CIC-CIT INSERM UMR 1099, Rennes, France

Abstract

Objectives

This study aimed to analyze the daily heart rate variability (HRV) in well-trained female cyclists during the 2017 Tour de France circuit and to relate it to the load and perceived exertion response.

Methods

Ten female cyclists volunteered to participate in the study. HRV was recorded with a portable heart rate monitor each morning at rest in supine (7 min.) and upright (7 min.) positions, as well as throughout each day’s stage. Pre-Tour baseline HRV recordings were made, as well as during the four weeks following completion of the Tour. Exercise daily load was assessed using the training impulse score (TRIMPS). Post-exercise rate of perceived exertion (RPE) was assessed daily using the Borg CR-10 scale.

Results

The results show a HRV imbalance, increase of sympathetic and decrease of vagal activities respectively, along the event that correlated with rate of perceived exertion (r = 0.46), training impulse score (r = 0.60), and kilometers (r = 046). The greatest change in HRV balance was observed the days after the greatest relative physical load. Mean heart rate and heart rate variability values returned to their baseline values one week after completion of the event.

Conclusions

Despite incomplete recovery from day-to-day, fatigue is not summative or augmented with each successive stage and its physical load. Just one week is sufficient to restore baseline values. Heart rate and HRV can be used as a tool to strategically plan the effort of female cyclists that participate in multi-stage events.

Introduction

Intensive endurance training is very demanding for the human organism and its regulatory systems [1]. Among them, cardiovascular control by the autonomic nervous system (ANS) presents alterations varying with the training versus recovery balance. Heart rate variability (HRV) represents sinus node modulation by the sympathetic and parasympathetic branches of the ANS [2]. HRV analysis is validated as a noninvasive method to study individual functional adaptations occurring to a given training stimulus in athletes [1,3].

HRV analysis is proposed as a valuable tool to study the athlete’s training versus recovery equilibrium and to detect early overreaching state that can decrease the athlete’s performance level [2, 3]. To our knowledge, most studies performed on this topic have focused on male athletes or mixed populations, and very few have studied female athletes. This omission is significant because the neural control of circulation differs with sex, especially before menopause [4]. HRV markers of sympathetic activity response after an orthostatic challenge test have been reported to be higher in male than in female athletes presented with a similar training load [5].

Moreover, HRV analysis in athletes has been focused mainly on the changes observed after acute post-exercise or throughout a training season in order to prevent fatigue and overtraining [1, 68]. HRV parameters changes induced by repeated days of endurance exercise have been scarcely studied, and only in male athletes [9]. Yet, it seems that the daily monitoring of HRV throughout a multi-day sporting event would be of interest for coaches to anticipate fatigue and to guide the athlete for his best final performance.

Therefore, the two objectives of this study were to describe a comprehensive characterization of resting heart rate (HR) and HRV changes in well-trained female cyclists during and following a multi-stage cycling event, and to propose an adapted method to follow fatigue and recovery in these athletes throughout this kind of sports event. We hypothesized that cyclists performing a Tour de France would incur a substantial HR increase and HRV decrease at rest along the days which would be worsened throughout the stages.

Materials and methods

Population

Ten healthy and well-trained (regional or national level) female cyclists coming from 4 countries (Belgium, France, Spain, Ukraine) completed the 21 stages of the men’s 2017 Tour de France, but one day before each stage of the official race. The event was performed without competition spirit. The cyclists completed the flat stages as a group, and each performed the time-trial and mountain stages at her own pace.

All participants gave informed written consent to participate in this study, which received the approval of the Rennes University Hospital Ethics Committee (Number 2013-A01524-41) and was conducted in accordance with the Declaration of Helsinki.

Pre-participation medical evaluation

All athletes had a medical examination before the event, including a physical exam, a resting electrocardiogram (ECG), a transthoracic echocardiogram, and an incremental maximal cardiopulmonary exercise test on an electronically braked cycle ergometer (Excalibur Sport, Lode, The Netherlands). The exercise protocol started with a warm-up period (100 W for 5 min and 150 W for 1 min) followed by a step load-increase of 25 W/min until exhaustion. None of the athletes ingested any contraceptive and their menstrual cycle was not controlled.

Description of the cycling event

The cycling event consisted of 21 stages and two rest days (Table 1). Due to the length of the stages, cyclists fed and hydrated continuously on the bike, according their needs. So individual food and fluid intake could not be controlled nor recorded during the stage nor between stages.

thumbnail
Table 1. Characteristics of the 2017 Tour de France stages.

https://doi.org/10.1371/journal.pone.0213472.t001

HRV study

Baseline HRV values were established from the average of four days before the first stage. HRV was monitored each day during the event and on the first day following completion of the Tour, then once per week during the four weeks post-Tour. Cyclists refrained from any intense efforts during both the pre- and post-Tour periods.

Recording RR samples.

The RR interval samples were recorded with a sampling rate of 1000 Hz with a HR monitor (Polar V800, Kempele, Finland) [10, 11]. Recordings were performed right after the cyclist woke up in the morning following an overnight fast, in a quiet, semi-darkened room, temperature range of 22–25°C. [2]. The RR samples were collected during two successive 7- minute periods, in supine and standing positions [1]. All athletes were familiarized with the monitor use.

HRV analysis.

The RR data recorded was downloaded via Polar FlowSync software for Mac version 2.6.4 (Polar, Kempele, Finland) and exported for analysis with the Kubios HRV Standard software v3.0.0 2 (Biosignal Analysis and Medical Imaging Group at the Department of Applied Physics, Kuopio, Finland) [10]. For the analysis, the last 5 minutes window for each position was used. All the ectopic beats were filtered with the artifact correction option of the software. A very low threshold was applied when needed. Both time and frequency domain analyses were performed [1]. The root-mean-square difference of successive normal RR intervals (RMSSD), which reflects HR vagal modulation, was calculated. The high (HF:0.15–0.40 Hz) and low (LF:0.04–0.15Hz) frequency domains were analyzed. The HF band reflects vagal modulation while the LF band indicates both sympathetic and parasympathetic influences []. RMSSD, HFnu, LFnu (normal units) absolute values and their difference between supine and standing positions were calculated. The normalized (or normalized unit) spectral indices are defined by the developers of the Kubios HRV Standard software v3.0.0 2 as HFnu = HF / (LF + HF) and LFnu = LF / (LF + HF) (Biosignal Analysis and Medical Imaging Group at the Department of Applied Physics, Kuopio, Finland) in accordance with the recommendations [2].

Load of exercise analysis

Each individual exercise daily load was calculated using the training impulse score (TRIMPS) method [12]. HR and GPS data were monitored continuously during each stage with the HR monitor. HR was divided into five zones, i.e. 50–60%, 60–70%, 70–80%, 80–90% and 90–100% of the individual maximal HR. Work load quantification was derived from the duration spent within the five HR zones [12].

Rate of perceived exertion

The individual rate of perceived exertion (RPE) was evaluated with the Borg CR-10 scale [13] within 30 minutes after the end of each stage.

Statistical analysis

Data are presented as mean ± SD. All analyses were performed using SPSS v.21 for Mac and STATISTICA v.7.1 for Windows. Normal Gaussian distribution of the data was verified by the Shapiro-Wilk test. Analysis of variance for repeated measurements was used. A Tukey’s post-hoc test was used to identify where the differences lie. In addition to the day-to-day effects on HR and HRV, the Tour was divided into 3 periods: period 1 (stages 1–9), period 2 (stages 10–15) and period 3 (stages 16–21), with rest days following stages 9 and 15.

Pearson’s product moment coefficient was calculated to assess the relationships among RPE, TRIMPS, the mileage of the stages, and the HR and HRV parameters. Significance was set at P<0.05.

Results

All cyclists successfully finished the Tour de France.

Table 2 shows the demographic characteristics and the result of the cardiopulmonary exercise test of the subjects.

HR, HRV, and workload values are presented for each stage (including resting days) of the cycling event in Table 3; and for each period in Table 4.

thumbnail
Table 3. Mean values (±SD) by stage for 10 female cyclists during the 2017 Tour de France.

https://doi.org/10.1371/journal.pone.0213472.t003

thumbnail
Table 4. Mean values (±SD) by period for 10 female cyclists during the 2017 Tour de France.

https://doi.org/10.1371/journal.pone.0213472.t004

In order to see the evolution of the subjects along the cycling event, Fig 1 presents HR, HRV, and workload values for each stage (including resting days).

thumbnail
Fig 1. Evolution of Heart rate and HRV parameters stage by stage.

Data are presented as mean group values. Absolute HR values different from baseline value* (P < 0.05); different from previous stage # (P < 0.05). HR standing- HR supine values different from baseline value $ (P < 0.05). When reading the differences, note that all the values were recorded on the morning of each stage, which reflects the load effect of the previous stage. HR: heart rate; RMSSD: root mean square of successive differences; HF: high frequency; LF: low frequency; ms: milliseconds, n.u.: normalized units; TT: time trial; FT: flat stage; MM: medium mountain stage; HM: high mountain stage; Post: post-cycling event period.

https://doi.org/10.1371/journal.pone.0213472.g001

To have a greater vision of the subjects’ evolution, Fig 2 illustrates HR and HRV parameters by periods.

thumbnail
Fig 2. Evolution of heart rate and HRV parameters by period.

Data are presented as mean group values for each period. Period 1: from stage 1 until stage 9 inclusive, period 2 from stage 10 until stage 15 inclusive, period 3 from stage 16 until stage 21 inclusive. Post: post-Tour period. Statistical differences: Different from baseline * (P < 0.05); different from the previous period # (P < 0.05). Standing-supine values different from baseline $ (P < 0.05); different from previous period £ (P < 0.05). HR: heart rate; RMSSD: root mean square of the successive differences; HF: high frequency; LF: low frequency; ms: milliseconds; n.u.: normalized units; TRIMPS: training impulse.

https://doi.org/10.1371/journal.pone.0213472.g002

Heart rate evolution

Supine HR during period 1 increased in comparison with its basal value after stage 2 (flat). This increase persisted until stage 9 with no difference among the stages. During period 2, HR increased after the first medium mountain stage (12th stage) and then decreased to its basal value. During period 3, HR was higher than its basal value only after the 17th (high-mountain) and 19th (flat) stages. For each of the three periods the global supine HR was higher than the basal value. HR was higher during period 1 than during periods 2 and 3, but no difference was observed between periods 2 and 3. The HR value returned to its basal values after each rest day and during the recovery period.

The standing HR values presented no difference during the cycling event, except after stage 2 (Fig 1). No differences were observed among the three periods. Positive correlations were observed between supine HR and the distance of the stage (r = 0.46, P = 0.037), TRIMPS (r = 0.60, P = 0.004), and RPE (r = 0.46, P = 0.038). The higher the workload, the higher supine HR was after a recovery night.

HFnu and LFnu evolution

Supine HFnu was decreased compared to its basal value from stage 2 until the first rest day, after which it went back to the basal value (Fig 1). Then, the HFnu value decreased from the 19th stage and returned to its basal value after the 21st stage. Lastly, the HFnu value was very close to its basal value after one week of recovery. As expected, LFnu mirrored the evolution of HFnu, with an increase of LFnu in compared to basal and recovery values. The changes in daily LFnu were more frequently different than the HFnu ones and overall LFnu was higher than HFnu during the entire cycling event. The opposite was observed only before the event start and during the recovery. No differences were observed between mean supine HFnu and LFnu values by period (Fig 2). The daily standing-supine differences for HFnu and LFnu were lower than their basal values during the entire Tour, except after the second rest day. No differences were observed after one week of recovery (Fig 1). The mean period analysis of standing-supine difference of these parameters showed the same significant differences as the daily values compared to baseline during all three periods, except for the second rest day (Fig 2). The standing HFnu values were lower and LFnu values were higher than supine, with no variation during the cycling event (Figs 1 and 2).

RMSSD evolution

The supine RMSSD global evolution showed the same trend as HFnu, even though more stages showed RMSSD declines (Fig 1). The same was observed with the periods analysis, with a more marked RMSSD decrease during all the periods in comparison to the basal and recovery periods (Fig 2). No significant variation was noted for the standing RMSSD.

The daily RMSSD standing-supine difference was lower than the basal value during the whole cycling period, except after the second rest day and after one week of recovery (Fig 2). There was a significant correlation between supine RMSSD and distance of the stages (r = -0.45, P = 0.0001). The mean period analysis of RMSSD standing-supine show a reduction compared to baseline but with no differences among periods (Fig 2).

Rate of perceived exertion evolution

There was a more marked decrease after the first rest day than after the second one (Table 3). The RPE evolution was positively related to TRIMPS (r = 0.61, P = 0.003) and to the distance of the stages (r = 0.91, P = 0.0001).

Discussion

To our knowledge, this study is the first to investigate the day-to-day resting HR and HRV responses in well-trained female cyclists during a multi-stage cycling event. We observed variable responses of HR and HRV parameters in relation to the distance of the stage throughout the Tour which partially confirms our hypothesis. The difference between standing and supine HR and HRV has a practical application for monitoring fatigue in a female cyclist.

Evolution of HR and HRV indices during the cycling event

As proposed, we studied the HR standing/supine difference and HRV indices responses to the orthostatic stress that reflects the adaptation of the sinus node [1]. Because significant changes concerned mainly resting values with no change for standing values (Figs 1 and 2), the discussion will concern resting HRV indices. Globally we observed an inversion of the ratio LFnu/HFnu with a value higher than 1 throughout the event, except for one value after the second rest day, in comparison to a LFnu/HFnu lower than 1 during the pre-race and post-race periods. During the first period (stages 1–9) of the event, both HR and HRV changes show a biphasic pattern. Indeed, in comparison to the basal value, period 1 showed three parts. No change was observed after stages 1 and 2, then after stages 3 to 5 a statistically significant decrease of HR response was noted. Concerning HRV indices, we observed a significant increase in the LFnu/HFnu ratio explained by a significant increase in sympathetic (LFnu) input associated with a significant decrease in parasympathetic (HFnu) input. Lastly, in stages 6 to 9 there was a less marked decrease of HR associated with a progressive increase in HFnu and a decrease in LFnu (Figs 1 and 2). The decrease of HR standing/supine difference observed during stages 3 to 9 was due to both an increase of resting HR and a decrease of standing HR. Throughout the second and third periods the HR as HRV responses to orthostatic stress were attenuated, with similar changes as seen during stages 6–9. The decrease of HR responses during these two periods was due to a lesser increase in standing HR than during the basal period, with fewer changes in supine HR. It must be noted that despite a very low level of mean TRIMPS during the last two stages of the Tour (20–21), 124 au and 276 au, respectively, we observed a marked decrease of HR response to orthostatic stress due to a low increase in HR and a marked increase in LFnu/HFnu when standing. This discrepancy between TRIMPS and HR response could be interpreted as a fatigue sign.

The biphasic response observed during period 1 was not linked to a workload difference (mean TRIMPS 1205 au and 1397 au, respectively, for stages 3–5 and 6–9; NS). The HR and HRV response we observed could be explained by the change in the cardiac preload reported after a few days of intense training, due to delayed hormonal responses [14, 15]. A similar biphasic response observed in left ventricular function has been reported after a four-day simulated multi-stage cycling [16]. The higher workload recorded during the first period compared to the second and third periods could explain the differences observed in HR and HRV responses.

The beneficial effect of a rest day for the responses to orthostatic stress appears clearly in Figs 1 and 2, for one day as well as during the recovery period post-Tour. Concerning the effects of the rest days proposed during the race, to our knowledge only one study investigated the effects on HRV indices of rest days (days 10 and 17) during the Vuelta a España performed by male cyclists. Similar to the present study, no difference was observed between rest days and pre-race HRV indices [9]. Our results show that after one week of post-race recovery, all HR parameters and HRV indices were similar to the pre-event basal values.

Supine RMSSD values, another parasympathetic parameter, mimics the globally responses of HFnu (Figs 1 and 2). After the two hardest stages of the event (HM stages 9 and 12, with 181 and 214 km), an acute decrease of supine and standing RMSSD values was observed (Fig 1).

To summarize, in comparison to the basal value, the day-to-day HR and HRV responses to an orthostatic active test vary with the duration of the cycling event. The temporary marked alterations of HR responses and of autonomic balance (LFnu/HFnu) observed during the first stages do not seem to be due to a real fatigue state but to acute stress stimuli (subjects not used to these distances). During the next stages both a stable and modest decrease of HR responses and an increase of LFnu/HFnu were observed. We did not observe a pattern of fatigue accumulation with the repetition of stages. However, at the end of the event (two last stages) marked alteration of the parameters studied was again observed. These results are in accordance with the most common ‘fatigue’ pattern described in athletes during intensive training [17]. This is confirmed by the quick and complete recovery after one week. Another observation is that supine HR response correlated with the distance of the stages, TRIMPS, and RPE.

Practical applications from the study

Several studies, reviews, and meta-analyses are in favor of the value of the HRV follow-up in male and female athletes to guide their training to prevent overreaching and overtraining [1, 3, 18, 19]. Our study corroborates that HRV is a useful non-invasive tool that can help to program training and identify fatigue in endurance cyclists. From a practical point of view, the results of this study can help to propose the more useful validated parameters for a daily practice during a multistage cycling event. Globally, given their lack of sensitivity the contribution of the isolated variations of the HR and HRV standing parameters used in this study appear as the less contributive (Figs 1 and 2). Concerning the HR survey, the difference between standing and supine HR after an active orthostatic test showed a good value and the variations of resting supine HR appear as a very simple tool. Concerning the HRV parameters, the use of two indices (RMSSD and HFnu) reflecting the response of the sinus node to the parasympathetic stimuli does not seem useful. The association of the two spectral indices LFnu and HFnu in supine position appears to be the most informative.

Limitations of the study

This study presents some limitations. First, the number of athletes studied was small. This limitation was mainly due to the specific and heavy daily logistic associated with this project that was not an official competition. Second, the cycling event was not a competitive one and our results need to be confirmed in competition because of the specific somatic stress linked to competition. Third, the daily recovery included only one night’s rest without a scientific adapted nutrition or massage session. A more scientific recovery might alter the observed results.

Conclusions

The characterization of the HRV response during the entire Tour de France adds new information concerning the fluctuating sympathovagal response of an ultra-endurance event. Despite incomplete recovery, the extent of cardiac suppression with each successive stage, and its physical load, is not summative or augmented. Just one week is enough to restore baseline values. These results suggest that well-trained female cyclists need only one week to recover from an effort such as the Tour de France circuit.

Acknowledgments

The authors gratefully acknowledge the enthusiastic participation of the subjects in this study and the generous technical support of Johan Cassirame and Armel Cretual. We also thank Dave Tanner for his help in English proofreading.

References

  1. 1. Buchheit M. Monitoring training status with HR measures: do all roads lead to Rome? Front Physiol. 2014;5:73. pmid:24578692
  2. 2. Heart rate variability: standards of measurement physiologicla interpretatioin and clinicla use. task Force of the European society of Cardiology and the North American Society of Pacing Circulation. 1996;93(5):1043–65.
  3. 3. Bellenger CR, Fuller JT, Thomson RL, Davison K, Robertson EY, Buckley JD. Monitoring Athletic Training Status Through Autonomic Heart Rate Regulation: A Systematic Review and Meta-Analysis. Sports medicine (Auckland, NZ). 2016;46(10):1461–86
  4. 4. Joyner MJ, Barnes JN, Hart EC, Wallin BG, Charkoudian N. Neural control of the circulation: how sex and age differences interact in humans. Comprehensive Physiology. 2015;5(1):193–215. pmid:25589269
  5. 5. Schafer D, Gjerdalen GF, Solberg EE, Khokhlova M, Badtieva V, Herzig D, et al. Sex differences in heart rate variability: a longitudinal study in international elite cross-country skiers. European journal of applied physiology. 2015;115(10):2107–14. pmid:26002403
  6. 6. da Silva VP, de Oliveira NA, Silveira H, Mello RG, Deslandes AC. Heart rate variability indexes as a marker of chronic adaptation in athletes: a systematic review. Annals of noninvasive electrocardiology: the official journal of the International Society for Holter and Noninvasive Electrocardiology, Inc. 2015;20(2):108–18.
  7. 7. Hautala AJ, Kiviniemi AM, Tulppo MP. Individual responses to aerobic exercise: the role of the autonomic nervous system. Neurosci Biobehav Rev. 2009;33(2):107–15. pmid:18514313
  8. 8. Perkins S, Jelinek H, Al-Aubaidy H, de Jong B. Immediate and long term effects of endurance and high intensity interval exercise on linear and nonlinear heart rate variability. Journal of science and medicine in sport. 2016;S1440-2440(16):30152–9.
  9. 9. Earnest CP, Jurca R, Church TS, Chicharro JL, Hoyos J, Lucia A. Relation between physical exertion and heart rate variability characteristics in professional cyclists during the Tour of Spain. Br J Sports Med. 2004;38:568–75 pmid:15388541
  10. 10. Caminal P, Sola F, Gomis P, Guasch E, Perera A, Soriano N, et al. Validity of the Polar V800 monitor for measuring heart rate variability in mountain running route conditions. European journal of applied physiology. 2018;118(3):669–77. pmid:29356949
  11. 11. Vanderlei LC, Silva RA, Pastre CM, Azevedo FM, Godoy MF. Comparison of the Polar S810i monitor and the ECG for the analysis of heart rate variability in the time and frequency domains. Braz J Med Biol Res. 2008;41(10):854–9. pmid:18853042
  12. 12. Edwards S. The heart rate monitor book. Sacramento: Fleet Feet Press; 1993.
  13. 13. Borg G. Borg’s perceived exertion and pain scales. Champaign, IL, US: Human Kinetics; 1998. viii, 104–viii, p.
  14. 14. Buchheit M, Laursen PB, Al Haddad H, Ahmaidi S. Exercise-induced plasma volume expansion and post-exercise parasympathetic reactivation. European journal of applied physiology. 2009;105(3):471–81. pmid:19009307
  15. 15. Convertino VA. Blood volume: its adaptation to endurance training. Medicine and science in sports and exercise. 1991;23(12):1338–48. pmid:1798375
  16. 16. Oosthuyse T, Avidon I, Likuwa I, Woodiwiss AJ. Progression of changes in left ventricular function during four days of simulated multi-stage cycling. European journal of applied physiology. 2012;112(6):2243–55. pmid:21997679
  17. 17. Schmitt L, Regnard J, Desmarets M, Mauny F, Mourot L, Fouillot JP, et al. Fatigue shifts and scatters heart rate variability in elite endurance athletes. PLoS One. 2013;8(8):e71588. pmid:23951198
  18. 18. Plews DJ, Laursen PB, Stanley J, Kilding AE, Buchheit M. Training adaptation and heart rate variability in elite endurance athletes: opening the door to effective monitoring. Sports medicine (Auckland, NZ). 2013;43(9):773–81.
  19. 19. Vesterinen V, Nummela A, Heikura I, Laine T, Hynynen E, Botella J, et al. Individual Endurance Training Prescription with Heart Rate Variability. Medicine and science in sports and exercise. 2016;48(7):1347–54. pmid:26909534