# Modeling Stress-Recovery Status Through Heart Rate Changes Along a Cycling Grand Tour

Anna Barrero, Anne Le Cunuder, Guy Carrault, François Carré, Frédéric Schnell, Solène Le Douairon Lahaye

► To cite this version:

Anna Barrero, Anne Le Cunuder, Guy Carrault, François Carré, Frédéric Schnell, et al., Modeling Stress-Recovery Status Through Heart Rate Changes Along a Cycling Grand Tour. Frontiers in Neuroscience, 2020, 14, 10.3389/fnins.2020.576308. [https://hal.science/hal-03095621v1](https://hal.science/hal-03095621v1)

HAL Id: hal-03095621

Submitted on 4 Jan 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Modeling Stress-Recovery Status Through Heart Rate Changes Along a Cycling Grand Tour

Anna Barrero¹,²*†, Anne Le Cunuder³†, Guy Carrault²,⁴, François Carré⁵, Frédéric Schnell⁵ and Solène Le Douairon Lahaye¹

¹ University of Rennes 2, M2S Laboratory, Rennes, France, ² CHU Rennes, Inserm CIC 1414, Rennes, France, ³ University of Lyon, Ecole Normale Supérieure, CNRS UMR 5672, Lyon, France, ⁴ University of Rennes 1, Inserm, LTSI – UMR 1099, Rennes, France, ⁵ University of Rennes, CHU Rennes, Inserm, LTSI – UMR 1099, Rennes, France

**Background:** Heart rate (HR) and HR variability (HRV) indices are established tools to detect abnormal recovery status in athletes. A low HR and vagally mediated HRV index change between supine and standing positions reflected a maladaptive training stress-recovery status.

**Objectives:** Our study was focused on a female multistage cycling event. Its overall aim was twofold: (1) quantify the correlation between (a) the change in HR and HRV indices during an active orthostatic test and (b) subjective/objective fatigue, physical load, and training level indicators; and (2) formulate a model predicting the stress-recovery status as indexed by $\Delta\overline{RR}$ and $\Delta\lnRMSSD$ (defined as the difference between standing and supine mean RR intervals and $\ln$RMSSD, respectively), based on subjective/objective fatigue indicators, physical load, and training levels.

**Methods:** Ten female cyclists traveled the route of the 2017 Tour de France, comprising 21 stages of 200 km on average. From 4 days before the beginning of the event itself, and until 1 day after its completion, every morning, each cyclist was subjected to HR and HRV measurements, first at rest in a supine position and then in a standing position. The correlation between HR and HRV indices and subjective/objective fatigue, physical load, and training level indicators was then computed. Finally, several multivariable linear models were tested to analyze the relationships between HR and HRV indices, fatigue, workload, and training level indicators.

**Results:** HR changes appeared as a reliable indicator of stress-recovery status. Fatigue, training level, and $\Delta\overline{RR}$ displayed a linear relationship. Among a large number of linear models tested, the best one to predict stress-recovery status was the following: $\Delta\overline{RR} = 1,249.37 + 12.32\dot{V}O_{2max} + 0.36 \text{ km}\cdot\text{week}^{-1} - 8.83 \text{ HR}_{max} - 5.8 \text{ RPE} - 28.41$ perceived fatigue with an adjusted $R^2 = 0.322$.

**Conclusion:** The proposed model can help to directly assess the adaptation status of an athlete from RR measurements and thus to anticipate a decrease in performance due to fatigue, particularly during a multistage endurance event.

**Keywords:** females, cycling, endurance, mathematical model, performance, heart rate variability, stress-recovery status

**OPEN ACCESS**

**Edited by:**
Sylvain Laborde,
German Sport University Cologne,
Germany

**Reviewed by:**
Dorota Zyśko,
Wroclaw Medical University, Poland
Antonio Roberto Zamunér,
Catholic University of Maule, Chile

*Correspondence:
Anna Barrero
annabarrero@gmail.com

† These authors have contributed
equally to this work

**Specialty section:**
This article was submitted to Autonomic Neuroscience, a section of the journal Frontiers in Neuroscience

**Received:** 25 June 2020
**Accepted:** 26 October 2020

**Published:** 02 December 2020

**Citation:**
Barrero A, Le Cunuder A, Carrault G, Carré F, Schnell F and Le Douairon Lahaye S (2020) Modeling Stress-Recovery Status Through Heart Rate Changes Along a Cycling Grand Tour.
Front. Neurosci. 14:576308.
doi: 10.3389/fnins.2020.576308
# INTRODUCTION

Physical training must combine workloads and recovery periods (Bishop et al., 2008). An optimal match between these two parameters is requested to avoid fatigue accumulation and thus reach the best individual physical performance level. Conversely, in the event of an imbalance between the two parameters, a state of overreaching or overtraining with marked drop in performance may occur (Kuipers, 1998). The individual optimal balance between work and recovery is difficult to achieve, especially for highly trained athletes. Thus, having relevant indicators of the exercise stress-recovery status could be a real advantage for individual athlete's monitoring. The aim of this study was to provide a new tool for coaches to help them to assess the athlete's individual fatigue level.

Several tools have been proposed to assess the fatigue level in daily routine of athletes (Buchheit, 2014; Halson, 2014). The most used are those providing an indirect evaluation of heart rate (HR) control by the autonomic nervous system (ANS) as analysis of spontaneous RR interval duration or HR variability (HRV) temporal and spectral indices (Bellenger et al., 2016), although the spectral ones present some limits (Buchheit, 2014). The HRV analysis can be performed either statically or dynamically, reflecting the heart adaptation to a physiological stress (Buchheit, 2014). Training-induced fatigue causes a more or less weak response of the ANS to external stimuli. Specifically, under normal training conditions (stress/recovery balance), intense strenuous exercise results in a decrease in vagally mediated HRV indices, followed by cardiac parasympathetic reactivation, which takes place during the recovery period (24-48 h) during which the cardiovascular system plays an important role in restoring body's homeostasis (Stanley et al., 2013). In training with a stress/recovery imbalance, a more marked and prolonged alteration of only vagally mediated HRV indices is reported (Bosquet et al., 2008; Plews et al., 2012). These results seem to support a link between fatigue and mainly the parasympathetic nervous system. Thus, HRV analysis, and especially focusing on parasympathetic' activity, represents a noninvasive method to track and record training status, "exercise readiness," and post-exercise fatigue in athletes (Buchheit, 2014; Bellenger et al., 2016). Indeed, in the parasympathetic system presenting the more marked impact on the post-exercise recovery, the vagally mediated HRV index variation is able to evaluate the cardiac response (Michael et al., 2017): the root-mean-square difference of successive normal RR intervals (RMSSD or its natural log, LnRMSSD) in the time domain is the most recommended (Plews et al., 2012; Buchheit, 2014; Bellenger et al., 2016), since it represents pure parasympathetic HR modulation (Malik and Camm, 1993).

However, it seems that the correct interpretation of HR or HRV fluctuations during the training period requires the comparison of these markers with other objective signs of fatigue to detect the risk of overreaching or overtraining (Stanley et al., 2013; Buchheit, 2014). From nocturnal ANS activity recordings in swimmers, Chalencon et al. (2012) found strong relationships between resting high-frequency (HF) power of HRV and 400-m freestyle time-trial performance of the next morning (Chalencon

et al., 2012). The higher the HF power values, the better the performance. Moreover, using the same protocol of ANS activity and performance recordings, the authors demonstrated that intensive training periods have a negative impact on both performance and HF due to fatigue (negative influence). In addition, modeling the effect of a 30-week training on swimmers' performance allowed for an accurate prediction of individual performance (Chalencon et al., 2015), supporting the relevance of a mathematical modeling of HRV in order to predict responses to training. To our knowledge, such mathematical analysis was only conducted over long training periods. Therefore, it appears relevant to also assess daily athlete's response on a shorter timescale, such as during a competitive event spreading over a few weeks (e.g., Grand Tours in cycling and World Football Cup) or during training camps, in which stress and fatigue affect day-to-day performances.

In a previous study on well-trained female cyclists, we have demonstrated that during a multistage event, HR and HRV indices evolved along the event in correlation with the daily physical load (Barrero et al., 2019). Briefly, we have observed a progressive increase of resting HR during the event and a progressive imbalance in the autonomic balance with an increase in the low-frequency (LF) power value that partly reflects the effects of sympathetic tone and a decrease in the HF and RMSSD values that reflects the parasympathetic effects. Our results also highlighted that variation in HR and HRV indices when changing from supine-to-standing position during an active orthostatic test is strongly correlated with the fatigue status (Barrero et al., 2019). The use of active orthostatic test has been recommended to detect fatigue in athletes, because of its ability to detect autonomic responses not observed with isolated supine or standing measures (Bosquet et al., 2008; Schmitt et al., 2015). This test explores the reactivities of sympathetic (excitation) and parasympathetic systems (withdrawal) in response to the position change (Taylor, 1994) and is the most widely used physiological maneuver for assessing neuro-vegetative and cardiovascular responsiveness. When standing from a supine position, the normal response is an increase in HR to maintain blood pressure (Tse et al., 2005). In well-trained athletes, with respect to supine rest values, after active stand-up, a marked increase of HR associated with a decrease in global HRV, HF, and RMSSD and with an increase of LF/HF ratio has been reported (Hynynen et al., 2008). These observations suggest that standing up normally induced mainly parasympathetic withdrawal. In overtrained athletes, an attenuation of parasympathetic and sympathetic activity during both supine and standing positions has been observed. Moreover, in these athletes, in response to stand-up position, the HR increases, and the decrease of total power, HF, LF, and RMSSD seemed lower than in non-overtrained athletes (Hynynen et al., 2008). Although a complex phenomenon (Schmitt et al., 2015), it is proposed that training induced-fatigue attenuates baroreflex response to change in position (Uusitalo et al., 2000). Thus, a low HR change between supine and standing positions reflects a maladaptive training stress-recovery status.

The multidimensional monitoring of recovery status has been underlined (Heidari et al., 2018), and it was
recommended to associate other markers with HR/HRV indices to more accurately detect the state of overtraining (Stanley et al., 2013; Buchheit, 2014). Because progression to overtraining syndrome appears to be associated with psychological states, the use of self-administered questionnaires on perceived physical and psychological well-being (WB) by athletes has been recommended (Saw et al., 2016). In this context, monitoring the evolution of several indices such as rate of perceived exertion (RPE), the WB, the sleep quality, and the delayed-onset muscle soreness (DOMS) during intensive periods of training has been proposed (Ouergui et al., 2020).

In our previous study, we observed a progressive decrease of HR and HRV responses induced by the active orthostatic test. These alterations were positively associated with the athlete's daily physical load and thus with fatigue (Barrero et al., 2019). However, the descriptive method we used to analyze the day-to-day HR and HRV indices did not allow to accurately predict the cyclist's individual adaptation ability in response to the exercise performed.

Thus, based on these previous data, the aim of the present study was first to investigate the individual correlation between the supine-standing difference of HR and HRV index values presented above, fatigue, physical load, and training level and second to formulate models relating two dependent HR variables ($\Delta RR$ and $\Delta LnRMSSD$ defined as the difference between standing and supine mean RR intervals and LnRMSSD, respectively) to subjective/objective fatigue, physical load, and training level indicators in order to investigate the parameters better predicting $\Delta RR$ and $\Delta LnRMSSD$. The added value of the current study is to propose a model allowing coaches to understand the adaptation ability of athletes to a multistage endurance event and helping them to anticipate a decrease in performance due to fatigue.

The design of our study has been based on two assumptions.

First, a relation exists between ANS alterations, mainly the cardiac vagal tone and early fatigue detection in athletes. Indeed, the development of fatigue in athletes is considered as a continuum process with, one side, the voluntary and controlled fatigue and, the opposite side, the uncontrolled fatigue, so-called overtraining (Meeusen et al., 2013). All the states of fatigue in athletes classically associate the same more or less marked symptoms: decrease in physical performance, neuro-endocrinal abnormalities, and psychological trouble such as irritability. Thus, the early detection of fatigue in athletes by use of HR and HRV analyses could be in line with the biological behavioral model (Grossman and Taylor, 2007). In accordance with previously published data concerning this first hypothesis, we assumed that due to the repetition of the successive endurance stages, we would observe a lower change of HR and of HR vagally mediated HRV between supine and standing positions during the orthostatic test. Moreover, we also hypothesized that HR and time-domain HRV vagal indices would present the higher relation with fatigue parameter levels.

Second, based on previous research (Ouergui et al., 2020), we hypothesized that the most suitable prediction models for both

$\Delta RR$ and $\Delta LnRMSSD$ should include parameters related on the one hand with physical performance, such as amount of training and level of performance, and on the other hand with levels of WB and of exertion feeling. We hypothesized that parameters of $\dot{VO}_{2max}$, WB, RPE, sleep quality, and DOMS could be included in the prediction model.

## MATERIALS AND METHODS

The first results of this prospective study have been previously published (Barrero et al., 2019).

This scientific project took place as part of the sports project "Donnons des elles au vélo J-1" aimed to promote women's cycling. This sport project included only 11 well-trained female cyclists, and of them, 10 participated in the whole study.

### Population

All were healthy, with no medical history of cardiovascular disease, and were not currently taking medication. At the time of this study, their weekly training mileage ranged from 100 to 250 km per week. All cyclists had a minimum of 2 years of competitive cycling experience. The characteristics of the athletes already presented (Barrero et al., 2019) are recalled in Table 1.

This study received the approval of our hospital ethics committee. After information was given, all participants provided written informed consent. The study was conducted in accordance with the "Good Clinical Practice" guidelines as laid down in the Declaration of Helsinki.

### Cycling Event

The characteristics of the multistage cycling event were previously described (Barrero et al., 2019). In brief, cyclists performed the 21 stages of men's 2017 Tour de France (TdT) 1 day before each stage of the official race. All event stages were performed without any spirit of competition or performance goal. The unique cyclists' objective was to complete 3,540 km of the TdT to promote women's cycling.

<table><caption>**TABLE 1** | Anthropometric and physical performance characteristics of the cyclists (mean ± SD).</caption><thead><tr><td>Characteristic</td><td>Value</td></tr></thead><tbody><tr><td>Age (years)</td><td>31.7 ± 4.7</td></tr><tr><td>Weight (kg)</td><td>57.0 ± 6.3</td></tr><tr><td>Height (m)</td><td>1.60 ± 0.1</td></tr><tr><td>BMI (kg·m<sup>-2</sup>)</td><td>21.4 ± 1.9</td></tr><tr><td>VO<sub>2max</sub> (ml·min·kg<sup>-1</sup>)</td><td>53.6 ± 5.2</td></tr><tr><td>Maximal power output (W)</td><td>285.5 ± 19.2</td></tr><tr><td>Relative maximal power output (W kg<sup>-1</sup>)</td><td>5.0 ± 0.5</td></tr><tr><td>Maximal HR (bpm)</td><td>185.7 ± 7.2</td></tr><tr><td>Previous cycling experience (years)</td><td>14.0 ± 8.9</td></tr><tr><td>Current training level (km·week<sup>-1</sup>)</td><td>187.5 ± 51.5</td></tr></tbody></table>

BMI, body mass index; $\dot{VO}_{2max}$, maximal oxygen uptake; HR, heart rate.
## Preliminary Testing

Each cyclist had a preparticipation medical evaluation with a clinical exam, a resting ECG (Mac 1600, GE Healthcare, Chicago, IL, United States) and an incremental maximal cardiopulmonary exercise test performed on an electronically braked cycle ergometer (Excalibur Sport, Lode, Netherlands) with continuous ECG and blood pressure monitoring and gas exchange analysis (Case system-Power cube, GE Healthcare, Chicago, IL, United States). The French Cycling Federation incremental exercise protocol was used. It 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. This preliminary testing took place 1 week before the first stage.

## RR Interval Recording and Analysis

The RR interval recording and HRV analysis protocols used were previously described (Barrero et al., 2019). Baseline pre-TdF RR intervals were collected daily during 4 days before the event to obtain a basal HRV state. Then RR recordings were performed every day of the multistage event. All resting recordings were made in the morning fasting, right after awakening, before the cyclist gets up. In order to avoid mental activity and stress and thus to place the cyclists in an optimal physiological rest state, RR recordings were performed in spontaneous breathing (Bernardi et al., 2000).

Briefly, all RR interval samples were recorded with a portable HR monitor (Polar V800, Kempele, Finland) during the two successive phases of the test: 7 min in a supine position followed by 7 min in standing position as recommended (Bourdillon et al., 2017). Individual RR recorded data were downloaded via Polar FlowSync software for mac version 2.6.4 (Polar, Kempele, Finland) and exported for later analysis. The Kubios HRV Standard software version 3.0.0 2 (Biosignal Analysis and Medical Imaging Group at the Department of Applied Physics, University of Kuopio, Kuopio, Finland) was used. For the HRV analysis, the last 5-min 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 (<5% of corrected beats). Both time and frequency domain HRV analysis were performed. The HRV spectrum is calculated with fast Fourier transform-based Welch's periodogram for spectral analysis. The RMSSD, which reflects cardiac vagal tone, was calculated. The HF (0.15-0.40 Hz) and LF (0.04-0.15 Hz) domains were analyzed. The HF band reflects cardiac vagal tone, 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 difference of the natural logarithm LnRMSSD between supine and standing positions was also studied. 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 (No authors listed, 1996).

## Daily Collection and Analysis of Heart Rate and Workload

### Heart Rate and GPS Recording

Both HR and GPS data were continuously registered with the Polar V800 (Polar, Kempele, Finland) portable monitor during each stage.

### Workload Evaluation

The objective load of the daily individual exercise performed was estimated from the HR collected during each stage. For this, the recorded HR values were divided into five zones according to the percentages (i.e., 50-60, 60-70, 70-80, 80-90, and 90-100%) of the individual maximum HR obtained during the maximum effort test. Then, the individual daily internal workload was estimated using the training pulse score (TRIMP) method (Edwards, 1993). In this model proposed by Edwards, the quantification of the workload is derived from the duration of the exercise maintained in the five HR zones described above (Edwards, 1993).

The individual RPE of the 21 stages was evaluated with Borg CR-10 scale within 30 min after each stage (Borg, 1990). RPE was indicative of the subjective load of each stage.

### Perceived Fatigue Evaluation (Well-being Questionnaire)

Cyclists answered every morning, during breakfast, a questionnaire with four questions focusing on perceived general fatigue, sleep quality, DOMS, and stress level. Each question scored on a 7-point (with 1 and 7 representing poor and very good WB ratings, respectively) scale. Overall WB was determined by summing the four questions' scores (Hooper et al., 1995).

### Statistical Analysis

The first step of the statistical analysis concerned the first hypothesis of the study. We studied the evolution of HR and HRV indices throughout the cycling event in our previously published study (Barrero et al., 2019). Then, we investigated the correlation between individual HR (RR interval duration) and HRV indices, and subjective/objective fatigue, physical load, and training level indicators, in order to establish the best fatigue marker. We analyzed more precisely the RMSSD of RR interval duration in supine (RMSSDsup) and in standing (RMSSDsta) positions, the difference ΔRMSSD = RMSSDsup−RMSSDsta and its natural logarithm (ΔLnRMSSD), the LF and HF indices of RR time series, and the ratio LF/HF. We also computed the RR duration mean values in supine (MeanRRsup) and standing (MeanRRsta) positions. Finally, the difference of the mean time interval between two successive heart beats (ΔRR̄) between supine and standing positions was computed (White, 1980; Tse et al., 2005): ΔRR̄ = RRsupine−RRstanding.

We then calculated the correlation coefficient between each of these indices and subjective/objective fatigue, physical load, and training level indicators (see section “Materials and Methods”).

The second step of the analysis concerned the second hypothesis of the study. To test this second hypothesis, we first had to characterize the influence of the cyclist's training level
on her adaptation all along the multistage event. To do this, the impact of each cycling stage was evaluated through the daily change observed on WB, with $\Delta$WB = $WB_{(staget)}$ - $WB_{(staget-1)}$. A low daily $\Delta$WB indicates that the cyclist was slightly impacted by the stage. Then, all daily $\Delta$WB was averaged and plotted for each cyclist as a function of her weekly pre-TdF training load (km·week⁻¹).

The final step of the statistical analysis concerned the second hypothesis of the study. We aimed to test models to investigate which objective/subjective fatigue, physical load, and training level indicators best predict $\Delta RR$ and $\Delta LnRMSSD$. For this purpose, different multivariable linear models (LMMs) were tested in order to estimate which set of variables best explains the $\Delta RR$ and $\Delta LnRMSSD$.

For both, each model was a linear combination of up to eight parameters related to training and physical condition (km·week⁻¹, $\dot{VO}_{2max}$, $HR_{max}$) on the one hand and to the perceived difficulty of the cycling stage and its impact on the cyclist (RPE, DOMS, perceived fatigue, quality of sleep, and stress) on the other hand:

$\Delta RR$ = f (km·week⁻¹, $\dot{VO}_{2max}$, $HR_{max}$, RPE, perceived fatigue, DOMS, quality of sleep, stress)
where f is a linear function.

$\Delta LnRMSSD$ = f (km·week⁻¹, $\dot{VO}_{2max}$, $HR_{max}$, RPE, perceived fatigue, DOMS, quality of sleep, stress)
where f is a linear function.

The training level, assessed trough $\dot{VO}_{2max}$ and km·week⁻¹, clearly influences the fatigue accumulated during the event.

For each parameter, we thus tested all $2^8 = 256$ possible models, covering all possible combinations of the eight parameters. The predictive power of each model was estimated through its Akaike information criterion (AIC) and the highest adjusted R², which penalize for including extra fitting parameters (Akaike, 1974). The best model corresponds to the lowest AIC value. In particular, if two models have the same R², the one having the less parameters has the lowest AIC.

The statistical significance of each parameter was estimated through its *p*-value with a significant value stated at *p* < 0.05. The dispersion of our data did not allow us to evaluate nonlinear effects beyond LMMs.

All statistical analyses were performed with Python 3.7 software version.

## RESULTS

All participants cycled for 21 consecutive stages, including two resting days. They all successfully completed the whole circuit of the TdF 2017 (3,540 km) (Barrero et al., 2019).

### Effects of Cycling Event on Heart Rate and Heart Rate Variability Indices

The evolution of HR and HRV indices throughout the cycling event has been previously published (Barrero et al., 2019). To summarize, resting supine HR increased progressively in comparison with its basal value during the multistage event, and the HR value returned to its basal values after each rest day.

On the other hand, standing HR values showed no significant evolution during the cycling event. All along the multistage event, we observed a progressive decrease in $\Delta RR$. Indeed, compared with its initial value (day 0), a progressive decrease in $\Delta RR$ was observed through successive stages. A small increase in $\Delta RR$ was noted after each resting day (days 10, 17, and 24) (Figure 1A).

Regarding HRV, a progressive imbalance in the autonomic balance marked with a decrease in cardiac vagal activity, evaluated through RMSSD and HF, was noted all along the cycling event. The daily RMSSD standing-supine difference was lower than the basal value during the multistage event.

### Correlations Between Heart Rate and Heart Rate Variability Indices, Fatigue, Physical Load, and Training Levels Indicators

These correlations concern our first hypothesis. Table 2 shows correlations between HRV and subjective/objective fatigue, physical load, and training level indicators. As expected, MeanRRsup, MeanRRsta, $\Delta$RMSSD, $\Delta$LnRMSSD, and $\Delta RR$ were negatively correlated with workload. We also observed significant negative correlations between RMSSDsup and workload indicators.

No correlation was observed between RMSSDsta and RPE ($r$ = 0.19, $p$ = 0.41), RMSSDsta and TRIMPS ($r$ = 0.047, $p$ = 0.84), nor RMSSDsta and distance of the stages ($r$ = 0.21, $p$ = 0.35). Lastly, the LF, HF, and LF/HF indices were not correlated with workload indicators ($r$ < 0.15, $p$ > 0.46). The $\Delta$RMSSD index showed significant negative correlation with TRIMPS ($r$ = -0.20, $p$ = 0.02), RPE ($r$ = -0.17, $p$ = 0.03), and KMS ($r$ = -0.43, $p$ = 0.05). We also computed the correlations of $\Delta$LnRMSSD with TRIMPS ($r$ = -0.23, $p$ = 0.009), RPE ($r$ = -0.04, $p$ = 0.12), KMS ($r$ = -0.41, $p$ = 0.06), perceived fatigue ($r$ = -0.24, $p$ = 0.002), DOMS ($r$ = -0.22, $p$ = 0.004), quality of sleep ($r$ = -0.06, $p$ = 0.447), and stress ($r$ = -0.22, $p$ = 0.004).

The $\Delta RR$ index showed significant correlation with workload markers with negative correlations with TRIMPS ($r$ = -0.22, $p$ = 0.009), RPE ($r$ = -0.26, $p$ < 0.05), daily stage distance ($r$ = -0.41, $p$ = 0.07), perceived fatigue ($r$ = -0.39, $p$ < 0.001), DOMS ($r$ = -0.40, $p$ < 0.001), quality of sleep ($r$ = -0.13, $p$ = 0.107), and stress ($r$ = -0.24, $p$ = 0.001) (Table 2). Otherwise, the overall $\Delta RR$ decrease observed all along the multistage cycling event was correlated ($r$ = 0.41, $p$ < 0.05) with the increase of RPE (Figure 1B).

Correlations between $\Delta RR$, $\Delta$RMSSD or $\Delta$LnRMSSD, and subjective/objective fatigue, physical load, and training level indicators are statically comparable (the *p*-value of the William's test we used to compare correlations is > 0.1 for each indicator).

### Impact of Training Level

To test the second hypothesis of the study, we first had to characterize the influence of the cyclist's training level on her adaptation all along the multistage event, and we considered the
<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Evolution of ΔRR stage by stage (A) and with daily workload assessed through rate of perceived exertion (RPE) score on CR-10 Borg scale (B). (A) ΔRR of stage n reflect stress/fatigue induced by stage n-1. Statistical differences with baseline: \*p < 0.05, \*\*p < 0.01, and \*\*\*p < 0.001.</figcaption></figure>

change in mean WB on all the stages as a function of individual weekly training load. The results are shown in **Figure 2**. Overall, consequent to the accumulated fatigue, the feeling of WB gradually decreased all along the event. This decrease was less marked in the most trained cyclists. However, the relationship observed was not linear but curvilinear. Finally, we must emphasize the great individual variability of the evolution of WB for the same weekly training load (i.e., 100 or 200 km·week⁻¹ in **Figure 2**).

## Multivariable Linear Models

To test our second hypothesis, we had lastly to consider 256 models for both ΔRR and ΔLnRMSSD. Among them, the four presenting the lowest AIC and the highest adjusted R² for ΔRR and ΔLnRMSSD are shown in **Tables 3, 4**, respectively. The fitting coefficients are reported with their *p*-value.

It must be underlined that the VO₂max value and the pre-TdF training load (km·week⁻¹) that were included in all models appeared to be the most relevant parameters to explain the
<table><caption>TABLE 2 | Correlation coefficients (and their p-values) between heart rate (RR) and heart rate variability indices and subjective/objective fatigue, physical load, and training levels indicators.</caption><thead><tr><th></th><th>MeanRRsup</th><th>MeanRRsta</th><th>RMSSDsup</th><th>RMSSDsta</th><th>&lt;ΔRR&gt;</th><th>ΔRMSSD</th><th>ΔLn RMSSD</th><th>LF(nu)</th><th>HF(nu)</th><th>LF(ms²)</th><th>HF(ms²)</th><th>LF/HF</th></tr></thead><tbody><tr><td>KMS</td><td>-0.44<br>(p = 0.04)</td><td>-0.34<br>(p = 0.04)</td><td>-0.35<br>(p = 0.02)</td><td>0.21<br>(p = 0.35)</td><td>-0.41<br>(p = 0.07)</td><td>-0.43<br>(p = 0.05)</td><td>-0.41<br>(p = 0.06)</td><td>-0.07<br>(p = 0.77)</td><td>0.07<br>(p = 0.77)</td><td>0.15<br>(p = 0.51)</td><td>0.13<br>(p = 0.58)</td><td>0.09<br>(p = 0.70)</td></tr><tr><td>TRIMPS</td><td>-0.27<br>(p = 0.001)</td><td>-0.12<br>(p = 0.175)</td><td>-0.20<br>(p = 0.019)</td><td>-0.05<br>(p = 0.532)</td><td>-0.22<br>(p = 0.009)</td><td>-0.20<br>(p = 0.020)</td><td>-0.23<br>(p = 0.005)</td><td>0.07<br>(p = 0.415)</td><td>-0.07<br>(p = 0.408)</td><td>-0.13<br>(p = 0.112)</td><td>-0.14<br>(p = 0.100)</td><td>0.06<br>(p = 0.465)</td></tr><tr><td>RPE</td><td>-0.17<br>(p = 0.033)</td><td>-0.04<br>(p = 0.574)</td><td>-0.26<br>(p = 0.001)</td><td>-0.06<br>(p = 0.444)</td><td>-0.26<br>(p = 0.001)</td><td>-0.17<br>(p = 0.029)</td><td>-0.16<br>(p = 0.042)</td><td>-0.10<br>(p = 0.200)</td><td>-0.11<br>(p = 0.150)</td><td>0.05<br>(p = 0.545)</td><td>-0.05<br>(p = 0.538)</td><td>-0.01<br>(p = 0.900)</td></tr><tr><td>Perceived fatigue</td><td>0.10<br>(p = 0.184)</td><td>0.20<br>(p = 0.011)</td><td>-0.17<br>(p = 0.023)</td><td>0.21<br>(p = 0.006)</td><td>-0.39<br>(p &lt; 0.001)</td><td>-0.01<br>(p = 0.903)</td><td>-0.24<br>(p = 0.002)</td><td>0.11<br>(p = 0.169)</td><td>0.13<br>(p = 0.091)</td><td>-0.02<br>(p = 0.756)</td><td>0.02<br>(p = 0.759)</td><td>0.03<br>(p = 0.701)</td></tr><tr><td>DOMS</td><td>0.16<br>(p = 0.032)</td><td>0.23<br>(p = 0.003)</td><td>-0.26<br>(p = 0.001)</td><td>0.10<br>(p = 0.195)</td><td>-0.40<br>(p &lt; 0.001)</td><td>0.05<br>(p = 0.500)</td><td>-0.22<br>(p = 0.004)</td><td>0.21<br>(p = 0.005)</td><td>0.20<br>(p = 0.009)</td><td>-0.01<br>(p = 0.906)</td><td>0.01<br>(p = 0.904)</td><td>-0.02<br>(p = 0.797)</td></tr><tr><td>Quality of sleep</td><td>0.24<br>(p = 0.002)</td><td>0.21<br>(p = 0.007)</td><td>-0.16<br>(p = 0.042)</td><td>-0.09<br>(p = 0.259)</td><td>-0.13<br>(p = 0.107)</td><td>0.14<br>(p = 0.080)</td><td>-0.06<br>(p = 0.447)</td><td>0.29<br>(p = 0.000)</td><td>0.22<br>(p = 0.005)</td><td>-0.00<br>(p = 0.988)</td><td>0.00<br>(p = 0.993)</td><td>-0.09<br>(p = 0.248)</td></tr><tr><td>Stress</td><td>0.15<br>(p = 0.056)</td><td>-0.14<br>(p = 0.069)</td><td>-0.41<br>(p = 0.000)</td><td>-0.28<br>(p = 0.000)</td><td>-0.24<br>(p = 0.001)</td><td>-0.07<br>(p = 0.392)</td><td>0.03<br>(p = 0.696)</td><td>-0.14<br>(p = 0.063)</td><td>-0.13<br>(p = 0.097)</td><td>-0.01<br>(p = 0.872)</td><td>0.01<br>(p = 0.879)</td><td>-0.11<br>(p = 0.154)</td></tr></tbody></table>

KMS, km of the daily stage; TRIMPS, training impulse; RPE, rate of perceived exertion; DOMS, delayed-onset muscle soreness.

<figure><img src="image_4.png" /><figcaption>FIGURE 2 | Averaged change in well-being (ΔWB) as a function of weekly training load (km week⁻¹). Blue circles represent change in well-being averaged on cyclists covering the same distance per week. The small blue points and the dotted lines represent individual change in WB.</figcaption></figure>

ΔR̄R and the ΔLnRMSSD observed. To further quantify the importance of these indicators of the physical condition, we also included in Tables 3, 4 the best model, which does not include the V̇O₂max (model 5), the number of km·week⁻¹ (model 6), or none of them (model 7). Clearly, beyond the multistage event-induced fatigue, the physical condition is the key to account the ΔR̄R and ΔLnRMSSD variations.

Specifically, the most relevant model to explain ΔR̄R can be expressed as a multi-linear function:

ΔR̄R = 1,249.37 + 12.32 V̇O₂max + 0.36 km·week⁻¹ - 8.83 HRmax - 5.8 RPE - 28.41 perceived fatigue where the p-values of each coefficient are given in Table 3.

And the best model to predict ΔLnRMSSD is:

ΔLnRMSSD = 1.647 + 0.042 V̇O₂max + 0.002 km·week⁻¹ - 0.017 HRmax - 0.019 RPE - 0.100 perceived fatigue + 0.092 stress

where the p-values of each coefficient are given in Table 4.

Finally, to formulate a more accessible model for coaches and team managers, we looked for a linear relationship between ΔR̄R, a single indicator of fatigue and physiologic constants excluding the V̇O₂max, which is not easily measurable. This model can be expressed as ΔR̄R = -14.36 WB - 2.06 HRmax + 0.35 km·week⁻¹ + 689.62. The R² adjusted and AIC observed for this model were, respectively, 0.194 and 2033.

Concerning ΔLnRMSSD, as we can see from adjusted R², models predict less variance than for ΔR̄R index (Table 4).

DISCUSSION

Intense physical training exposes the athlete to the risk of overreaching or of overtraining, partly due to an imbalance between training and recovery (Bishop et al., 2008). The information from overreaching and overtraining markers, especially based on HR and HRV analyses, seems reinforced when it is associated with other fatigue parameters (Stanley et al., 2013; Buchheit, 2014). Our study aimed at modeling the relationship between the level of fatigue reported by well-trained female cyclists during a multistage cycling event, their physical load,
<table><caption>**TABLE 3** | Linear coefficients, their corresponding *p*-values (in parentheses), adjusted R² coefficients, and Akaike information criterion (AIC) of the different linear multivariable mixed models tested to explain the difference between standing and supine RR duration of cyclists.</caption><thead><tr><th></th><th>R1</th><th>R2</th><th>R3</th><th>R4</th><th>R5</th><th>R6</th><th>R7</th></tr></thead><tbody><tr><td>V̇O<sub>2max</sub></td><td>12.32** (&lt;0.001)</td><td>11.39** (&lt;0.001)</td><td>12.59** (&lt;0.001)</td><td>12.40** (&lt;0.001)</td><td></td><td>11.78** (&lt;0.001)</td><td></td></tr><tr><td>Training load (km·week<sup>-1</sup>)</td><td>0.36 (0.021)</td><td>0.36 (0.019)</td><td>0.35 (0.023)</td><td>0.36 (0.023)</td><td>0.35 (0.033)</td><td></td><td></td></tr><tr><td>HR<sub>max</sub></td><td>-8.83** (&lt;0.001)</td><td>-8.13** (&lt;0.001)</td><td>-9.02** (&lt;0.001)</td><td>-8.79** (&lt;0.001)</td><td></td><td>-8.97** (&lt;0.001)</td><td>-2.52 (0.075)</td></tr><tr><td>RPE</td><td>-5.90 (0.096)</td><td>-5.12 (0.16)</td><td>-5.99 (0.093)</td><td>-5.78 (0.156)</td><td>-4.71 (0.222)</td><td>-6.02 (0.094)</td><td>-5.94 (0.129)</td></tr><tr><td>Perceived fatigue</td><td>-28.41** (&lt;0.001)</td><td>-24.55* (0.006)</td><td>-28.80** (&lt;0.001)</td><td>-28.88 (&lt;0.001)</td><td>-24.58* (0.009)</td><td>-29.62** (&lt;0.001)</td><td>-19.32 (0.058)</td></tr><tr><td>DOMS</td><td></td><td>-6.20 (0.365)</td><td></td><td></td><td>-19.13* (0.004)</td><td></td><td>-19.33* (0.005)</td></tr><tr><td>Sleep quality</td><td></td><td></td><td></td><td>1.15 (0.85)</td><td></td><td></td><td>-10.81 (p = 0.10)</td></tr><tr><td>Stress</td><td></td><td></td><td>2.23 (0.748)</td><td></td><td></td><td></td><td></td></tr><tr><td>Constant</td><td>1,249.37** (&lt;0.001)</td><td>1,166.96** (&lt;0.001)</td><td>1,266.30** (&lt;0.001)</td><td>1,234.21** (&lt;0.001)</td><td>314.60** (&lt;0.001)</td><td>1,375.93** (&lt;0.001)</td><td>860.54** (&lt;0.001)</td></tr><tr><td>Adjusted R<sup>2</sup></td><td>0.322</td><td>0.321</td><td>0.318</td><td>0.319</td><td>0.234</td><td>0.302</td><td>0.223</td></tr><tr><td>AIC</td><td>1,805</td><td>1,807</td><td>1,807</td><td>1,807</td><td>1,823</td><td>1,824</td><td>1,825</td></tr></tbody></table>

Adjusted R² and AIC are used to estimate the predictive power of each model. These criteria are used instead of R² in order to penalize including extra fitting parameters. Training load corresponds to the weekly pre-TdF training load of the cyclists.

VO₂max, maximal oxygen uptake; HRmax, maximal heart rate; RPE, rate of perceived exertion; WB, well-being; DOMS, delayed-onset muscle soreness.

<table><caption>**TABLE 4** | Linear coefficients, their corresponding *p*-values (in parentheses), adjusted R² coefficients, and Akaike information criterion (AIC) of the different linear multivariable mixed models tested to explain the difference between standing and supine LnRMSSD of cyclists.</caption><thead><tr><th></th><th>R1"</th><th>R2"</th><th>R3"</th><th>R4"</th><th>R5"</th><th>R6"</th><th>R7"</th></tr></thead><tbody><tr><td>V̇O<sub>2max</sub></td><td>0.042* (0.002)</td><td>0.026* (0.004)</td><td>0.024* (0.008)</td><td>0.038* (0.009)</td><td></td><td>0.040* (0.003)</td><td></td></tr><tr><td>Training load (km·week<sup>-1</sup>)</td><td>0.002 (0.010)</td><td>0.002 (0.013)</td><td>0.002 (0.011)</td><td>0.002 (0.014)</td><td>0.0017 (0.036)</td><td></td><td></td></tr><tr><td>HR<sub>max</sub></td><td>-0.017 (0.094)</td><td></td><td></td><td>-0.014 (0.217)</td><td></td><td>-0.018 (0.081)</td><td></td></tr><tr><td>RPE</td><td>-0.019 (0.308)</td><td>-0.019 (0.303)</td><td>-0.013 (0.487)</td><td>-0.015 (0.421)</td><td>-0.016 (0.407)</td><td>-0.020 (0.297)</td><td>-0.017 (0.392)</td></tr><tr><td>Perceived fatigue</td><td>-0.100 (0.013)</td><td>-0.114* (0.005)</td><td>-0.081 (0.083)</td><td>-0.084 (0.072)</td><td>-0.075 (0.118)</td><td>-0.108* (0.009)</td><td>0.085 (0.078)</td></tr><tr><td>DOMS</td><td></td><td></td><td>-0.044 (0.018)</td><td>-0.026 (0.464)</td><td>-0.059 (0.079)</td><td></td><td>-0.054 (0.109)</td></tr><tr><td>Sleep quality</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Stress</td><td>0.092 (0.012)</td><td>0.074 (0.034)</td><td>0.074 (0.033)</td><td>0.089 (0.016)</td><td>0.059 (0.092)</td><td>0.095 (0.011)</td><td>0.063 (0.078)</td></tr><tr><td>Constant</td><td>1.647 (0.238)</td><td>-0.495 (0.388)</td><td>-0.427 (0.457)</td><td>1.271 (0.393)</td><td>0.982** (&lt;0.001)</td><td>2.365 (0.090)</td><td>1.331** (&lt;0.001)</td></tr><tr><td>Adjusted R<sup>2</sup></td><td>0.141</td><td>0.130</td><td>0.135</td><td>0.139</td><td>0.098</td><td>0.111</td><td>0.076</td></tr><tr><td>AIC</td><td>229.9</td><td>230.9</td><td>231.0</td><td>231.4</td><td>236.4</td><td>234.1</td><td>239.1</td></tr></tbody></table>

Adjusted R² and AIC are used to estimate the predictive power of each model. These criteria are used instead of R² in order to penalize including extra fitting parameters. Training load corresponds to the weekly pre-TdF training load of the cyclists.

VO₂max, maximal oxygen uptake; HRmax, maximal heart rate; RPE, rate of perceived exertion; WB, well-being; DOMS, delayed-onset muscle soreness.

their training level, and the variations in HRV and HR rate indices in response to an active orthostatic test.

Regarding our first hypothesis, we noted a lower change of HR and of HR vagally mediated HRV between supine and standing positions during the orthostatic test. However, only HR and RMSSD, a time domain HRV index, were correlated with subjective/objective fatigue, physical load, and training level indicators. Lastly, the index ΔRR, defined as the difference between the average RR intervals measured in a supine position and then in a standing position, appeared a new indicator of stress/recovery status. These results confirmed our first hypothesis.

Regarding our second hypothesis, we have then demonstrated that ΔRR and ΔLnRMSSD could be modeled as a linear function of training volume, V̇O₂max and fatigue level, assessed through the RPE and the WB questionnaire. Thus, the results observed

confirmed only partly our second hypothesis, because DOMS and sleep quality previously proposed (Ouergui et al., 2020) did not provide major information to specify the stress-recovery status in the multistage of endurance event studied (Table 3). Regarding ΔLnRMSSD, models predict less variance than for ΔRR index. This result underlines the interest of ΔRR in stress-recovery status prediction.

## Respective Values of Heart Rate Variability Indices

From all the HRV indices we used, time domain's is the one that was the best correlated with workload and fatigue parameters (Table 2). Limits of spectral indices have been previously reported (Buchheit, 2014). Our results confirm also that the time-domain markers of parasympathetic effects seemed to be better adapted

<footnote>*p-value < 10⁻², **p-value < 10⁻³.</footnote>

<footnote>*p-value < 10⁻², **p-value < 10⁻³.</footnote>
to explore fatigue level (Buchheit, 2014; Fazackerley et al., 2019). This can be explained by the fact that parasympathetic nervous system is implied in self-regulation mechanisms, which are critical for adaptation (Laborde et al., 2017). Lastly, among these indices, the ΔRR index was significantly correlated with subjective/objective fatigue, physical load, and training level indicators, as ΔRMSSD and ΔLnRMSSD. Previous reviews were focused on interests and limits of HR and HRV measures on monitoring training status (Buchheit, 2014; Bellenger et al., 2016). If Buchheit considers resting HRV (more precisely RMSSD) as the HR measure more sensitive to fatigue (Buchheit, 2014), Bellenger et al. (2016) underlined some limits of this ANS status analysis. Indeed, in their meta-analysis, Bellenger et al. (2016) have observed that overload training had little effect on resting HRV due to various effects on vagally mediated HRV indices. The authors explained that the disagreement between studies may be the result of methodological issues. In addition to these methodological aspects, HRV analysis appears as a complex process due to different fatigue-induced alterations of HRV pattern (Schmitt et al., 2015). Given these limits of resting HRV, it is interesting to bring out new tools in ANS status evaluation. Our study therefore outlines a new indicator of stress/recovery status, the ΔRR. Based on our results, a low ΔRR, translating a low HR change between supine and standing positions, could mean a lack of post-exercise recovery. This observation is in accordance with the increase cardiac sympathetic modulation during supine rest and attenuated baroreflex response to change position observed by Uusitalo et al. (2000) in heavily trained females. The decrease observed here seems explained by an increase in supine HR without modification of standing HR (Barrero et al., 2019). This observation may be due to a decrease in parasympathetic and/or an increase in sympathetic HR influence (Schmitt et al., 2015). To our knowledge, the precise cause of this observation, decrease in sensitivity, and/or density of sinus cell membrane receptors or other one has not been formally demonstrated (Bellenger et al., 2016).

## Mathematical Model

As specified in the *Introduction*, we tested eight parameters that quantify physical load, training level, and fatigue indicators. The $\dot{V}O_{2max}$ is widely used to assess both physical capacity and training level in endurance sports, but it is typically not repeatedly measured during a sporting season. Therefore, we included the training volume, summarized by the mean number of kilometers per week during the training period prior to the event. The chosen fatigue indicators reflect the internal load: RPE is commonly used to evaluate the perceived difficulty of an exercise and appeared strongly correlated to ΔRR, and the DOMS is specifically targeted at muscular fatigue; the perceived fatigue reflects general tiredness; quality of sleep and stress are associated with physical and mental fatigue.

We first studied the influence of the individual training level on the adaptation throughout the multistage cycling event. This training level depends on two main factors: the number of years of practice and the quality of training carried out during the weeks preceding the event. We noted no correlation between the number of years of practice and WB during the event. On the

other hand (see **Figure 2**), except for the least trained cyclist, we observed a positive correlation between the training volume per week before the event and the adaptation all along the event. We therefore observed a predominant influence of recent training on the level of exercise tolerance during the multistage event.

Then, several multivariable models were tested. We showed that the indices ΔRR and ΔLnRMSSD could be modeled linearly as a function of three main individual variables: training volume, $\dot{V}O_{2max}$, and fatigue assessed through RPE and WB questionnaire. These indices appear to be relevant indicators of the adaptation ability along multistage events. However, we have noticed that linear models based on ΔRR have more predictive power than those based on ΔLnRMSSD. Moreover, the individual ΔRR and ΔLnRMSSD appear to be reliable indicators of both athlete's training level and fatigue level. Indeed, we observed a positive correlation between ΔRR (respectively, ΔLnRMSSD) and $\dot{V}O_{2max}$, which reflects the training level and a negative correlation between ΔRR (respectively, ΔLnRMSSD) and fatigue, RPE, WB, and maximal HR. The positive correlation we observed with the individual fitness and recovery of altered autonomic regulation after prolonged exercise confirms previous observations (Hautala et al., 2001; Fazackerley et al., 2019). Finally, the negative correlation with maximal HR confirms that the latter decreases with chronic endurance training (Bailey and Davies, 1999).

## Practical Applications

A decrease in the value of HRV indices is a marker of weak adaptability of the cardiovascular system to stress conditions that it faces (Michael et al., 2017); for example, the supine vagally mediated HRV parameters (RMSSD, total spectral power, and HF but not LF/HF) were lower in athletes identified in the fatigue state compared with the nonfatigue one (Plews et al., 2012; Schmitt et al., 2015). The active orthostatic test, a well-described marked physiological stress, is recommended to study HRV in athletes (Uusitalo et al., 2000; Hynynen et al., 2008; Schmitt et al., 2015).

From our results, it seems that the higher is ΔRR, HR difference between supine and standing positions, the better is the cardiovascular adaptability to orthostatic stress. Conversely, a low ΔRR means stress-recovery imbalance as described in the paragraph *Respective Values of Heart Rate Variability Indices*. The negative impact of a Grand Tour on physical performance, mood, and WB of competitive cyclists is well reported, and a study performed with professional male cyclists during the Vuelta a España has noted that changes in supine HRV were inversely correlated to the exercise level (Earnest et al., 2004; Lastella et al., 2015; Rodríguez-Marroyo et al., 2017). Our results confirmed this observation. The impact of the preintervention physiological status on HRV alteration has also been proposed (Lastella et al., 2015; Rodríguez-Marroyo et al., 2017). Our data also underlined the importance of the pre-event training level to explain the HR adaptability change during a cyclist's multistage event.

To our knowledge, it is the first time that the HR changes observed during an easy physiological test are modeled as a function of two physical parameters, training level and $\dot{V}O_{2max}$, and one psychological parameter, fatigue. The model
proposed allowed us to understand the ability to adapt to a repeated endurance exercise measuring the mean RR interval changes observed between supine and standing positions in well-trained female cyclists. If the measured $\Delta RR$ is lower than predicted by the model, we could conclude that the imbalance of stress-recovery status is higher than perceived by the athlete in this context of a cyclist multistage event.

Thus, our findings support the use of $\Delta RR$ monitoring to quantify training load, as $\Delta RR$ can be directly predicted from fatigue and training level indicators. The use of $\Delta RR$ monitoring can help coaches and athletes to make strategic decisions during a multistage long-duration event. It should be noted that the proposed model can be used with unsophisticated HR monitors (i.e., those recording only RR average values). This therefore makes it accessible to a majority of coaches and athletes, interested in the scientific approach of the training and performance monitoring. Finally, we also presented a model more accessible for coaches and team managers, connecting HR changes to a single indicator of fatigue. However, in accordance with its $R^2$ adjusted and AIC values, the robustness of this model appears low, and it should be used with caution.

## Study Limitations

This study presents three main limitations. First, the small population sample studied can reduce the predictive power of the proposed model. However, the daily evolution of HRV in female athletes has been scarcely studied, and the model we proposed seems to be easy to use to follow the individual training and performance level. Second, our study focuses on well-trained female cyclists, and the validity of our model deserves to be confirmed in other sports. Specific studies are also needed in endurance male athletes, because HRV gender's difference has been reported (Schäfer et al., 2015). Third, as some HRV indices are closely related and dependent on the individual's breathing frequency during recording (No authors listed, 1996), it could have been relevant to impose controlled breathing (Gregoire et al., 1996), although this induces a mental activity and stress (Bernardi et al., 2000). At least, monitoring of respiratory rate could have been considered.

## REFERENCES

Akaike, H. (1974). A new look at the statistical model identification. *IEEE Trans. Auto. Control* 19, 716-723. doi: 10.1109/tac.1974.1100705

Bailey, D. M., and Davies, B. (1999). Decreased chronotropic drive as an adaptation to chronic exercise; possible mechanisms. *Int. J. Sports Med.* 20, 219-221. doi: 10.1055/s-2007-971120

Barrero, A., Schnell, F., Carrault, G., Kervio, G., Matelot, D., Carre, F., et al. (2019). Daily fatigue-recovery balance monitoring with heart rate variability in well-trained female cyclists on the Tour de France circuit. *PLoS One* 14:e0213472. doi: 10.1371/journal.pone.0213472

Bellenger, C. R., Fuller, J. T., Thomson, R. L., Davison, K., Robertson, E. Y., and Buckley, J. D. (2016). Monitoring athletic training status through autonomic heart rate regulation: a systematic review and meta-analysis. *Sports Med.* 10, 1461-1486. doi: 10.1007/s40279-016-0484-2

## CONCLUSION

From our data on well-trained female cyclists, we introduced a new indicator of post-endurance exercise recovery, the $\Delta RR$ based on the change of mean HR observed during an orthostatic active test. This index is influenced by the training level and by the $VO_{2max}$ of the athlete. The proposed quantitative model can help to assess the adaptation ability of an athlete and thus to anticipate a decrease in endurance performance due to fatigue, particularly during a long-duration multistage cycling event. Investigating larger populations of athletes, included in other sports than cycling, represents an exciting perspective for future studies.

## DATA AVAILABILITY STATEMENT

The original contributions presented in this study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

## ETHICS STATEMENT

The studies involving human participants were reviewed and approved by CHU Rennes, France. The patients/participants provided their written informed consent to participate in this study.

## AUTHOR CONTRIBUTIONS

AB, SL, GC, FC, and FS performed the experimental conception and design. AB and SL performed the experiments. AB and ALC analyzed the data. AB, ALC, SL, FC, and FS written the manuscript. All authors read and approved the final manuscript.

## 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.

Bernardi, L., Wdowczyk-Szulc, J., Valenti, C., Castoldi, S., Passino, C., Spadacini, G., et al. (2000). Effects of controlled breathing, mental activity and mental stress with or without verbalization on heart rate variability. *J Am Coll. Cardiol.* 35, 1462-1469. doi: 10.1016/s0735-1097(00)00595-7

Bishop, P. A., Jones, E., and Woods, A. K. (2008). Recovery from training: a brief review. *J. Strength Cond. Res.* 22, 1015-1024. doi: 10.1519/jsc.0b013e31816eb518

Borg, G. (1990). Psychophysical scaling with applications in physical work and the perception of exertion. *Scand. J. Work Environ. Health* 16(Suppl. 1), 55-58. doi: 10.5271/sjweh.1815

Bosquet, L., Merkari, S., Arvisais, D., and Aubert, A. E. (2008). Is heart rate a convenient tool to monitor over-reaching? A systematic review of the literature. *Br. J. Sports Med.* 42, 709-714. doi: 10.1136/bjsm.2007.042200

Bourdillon, N., Schmitt, L., Yazdani, S., Vesin, J. M., and Millet, G. P. (2017). Minimal window duration for accurate HRV recording in athletes. *Front. Neurosci.* 11:456. doi: 10.3389/fnins.2017.00456
Buchheit, M. (2014). Monitoring training status with HR measures: do all roads lead to Rome? *Front. Physiol.* 5:73. doi: 10.3389/fphys.2014.00073

Chalencon, S., Busso, T., Lacour, J. R., Garet, M., Pichot, V., Connes, P., et al. (2012). A model for the training effects in swimming demonstrates a strong relationship between parasympathetic activity, performance and index of fatigue. *PLoS One* 7:e52636. doi: 10.1371/journal.pone.005263

Chalencon, S., Pichot, V., Roche, F., Lacour, J. R., Garet, M., Connes, P., et al. (2015). Modeling of performance and ANS activity for predicting future responses to training. *Eur. J. Appl. Physiol.* 115, 589–596. doi: 10.1007/s00421-014-3035-2

Earnest, C. P., Jurca, R., Church, T. S., Chicharro, J. L., Hoyos, J., and Lucia, A. (2004). Relation between physical exertion and heart rate variability characteristics in professional cyclists during the Tour of Spain. *Br. J. Sports Med.* 38, 568–575. doi: 10.1136/bjsm.2003.005140

Edwards, S. (1993). "High performance training and racing," in *The Heart Rate Monitor Book*, ed. S. Edwards (Sacramento: Feet Fleet Press), 113–123.

Fazackerley, L. A., Fell, J. W., and Kitic, C. M. (2019). The effect of an ultra-endurance running race on heart rate variability. *Eur. J. Appl. Physiol.* 119, 2001–2009. doi: 10.1007/s00421-019-04187-6

Gregoire, J., Tuck, S., Yamamoto, Y., and Hughson, R. L. (1996). Heart rate variability at rest and exercise: influence of age, gender, and physical training. *Can. J. Appl. Physiol.* 21, 455–470. doi: 10.1139/h96-040

Grossman, P., and Taylor, E. W. (2007). Toward understanding respiratory sinus arrhythmia: relations to cardiac vagal tone, evolution and biobehavioral functions. *Biol. Psychol.* 74, 263–285. doi: 10.1016/j.biopsycho.2005.11.014

Halson, S. L. (2014). Monitoring training load to understand fatigue in athletes. *Sports Med.* 44(Suppl 2), S139-S147. doi: 10.1007/s40279-014-0253-z

Hautala, A., Tulppo, M. P., Mäkikallio, T. H., Laukkanen, R., Nissilä, S., and Huikuri, H. V. (2001). Changes in cardiac autonomic regulation after prolonged maximal exercise. *Clin. Physiol.* 21, 238–245. doi: 10.1046/j.1365-2281.2001.00309.x

Heidari, J., Beckmann, J., Bertollo, M., Brink, M., Kallus, W., Robazza, C., et al. (2018). Multidimensional monitoring of recovery status and implications for performance. *Int. J. Sports Physiol. Perform.* 15, 1–24. doi: 10.1123/ijspp.2017-0669

Hooper, S. L., Mackinnon, L. T., Howard, A., Gordon, R. D., and Bachmann, A. W. (1995). Markers for monitoring overtraining and recovery. *Med. Sci. Sports Exerc.* 27, 106–112.

Hynynen, E., Uusitalo, A., Konttinen, N., and Rusko, H. (2008). Cardiac autonomic responses to standing up and cognitive task in overtrained athletes. *Int. J. Sports Med.* 29, 552–558. doi: 10.1055/s-2007-989286

Kuipers, H. (1998). Training and overtraining: an introduction. *Med. Sci. Sports Exerc.* 30, 1137–1139. doi: 10.1097/0005768-199807000-00018

Laborde, S., Mosley, E., and Thayer, J. F. (2017). Heart rate variability and cardiac vagal tone in psychophysiological research - recommendations for experiment planning, data analysis, and data reporting. *Front. Physiol.* 8:213. doi: 10.3389/fpsyg.2017.00213

Lastella, M., Roach, G. D., Halson, S. L., Martin, D. T., West, N. P., and Sargent, C. (2015). The impact of a simulated grand tour on sleep, mood, and well-being of competitive cyclists. *J. Sports Med. Phys. Fitness* 55, 1555–1564.

Malik, M., and Camm, A. J. (1993). Components of heart rate variability-what they really mean and what we really measure. *Am. J. Cardiol.* 72, 821–822. doi: 10.1016/0002-9149(93)91070-x

Meeusen, R., Duclos, M., Foster, C., Fry, A., Gleeson, M., Nieman, D., et al. (2013). Prevention, diagnosis, and treatment of the overtraining syndrome: joint consensus statement of the European College of Sport Science and the American College of Sports Medicine. *Med. Sci. Sports Exer.* 45, 186–205. doi: 10.1249/MSS.0b013e318279a10a

Michael, S., Graham, K. S., and Davis, G. M. (2017). Cardiac autonomic responses during exercise and post-exercise recovery using heart rate variability and systolic time intervals-a review. *Front. Physiol.* 8:301. doi: 10.3389/fphys.2017.00301

No authors listed (1996). Heart rate variability: standards of measurement physiological interpretation and clinical use. task Force of the European society of Cardiology and the North American Society of Pacing and Electrophysiology. *Circulation* 93, 1043–1065. doi: 10.1161/01.cir.93.5.1043

Ouergui, I., Ardigò, L. P., Selmi, O., Levitt, D. E., Chtourou, H., Bouassida, A., et al. (2020). Changes in perceived exertion, well-being, and recovery during specific judo training: impact of training period and exercise modality. *Front. Physiol.* 11:931. doi: 10.3389/fphys10.3389/fphys.2020.00931

Plews, D. J., Laursen, P. B., Kilding, E. A., and Buchheit, M. (2012). Heart rate variability in elite triathletes, is variation in variability the key to effective training? A case comparison. *Eur. J. Appl. Physiol.* 112, 3729–3741. doi: 10.1007/s00421-012-2354-4

Rodríguez-Marroyo, J. A., Villa, J. G., Pernía, R., and Foster, C. (2017). Decrement in professional cyclists' performance after a grand tour. *Int. J. Sports Physiol. Perform.* 12, 1348–1355. doi: 10.1123/ijspp.2016-0294

Saw, A. E., Main, L. C., and Gastin, P. B. (2016). Monitoring the athlete training response: subjective self-reported measures trump commonly used objective measures: a systematic review. *Br. J. Sports Med.* 50, 281–291. doi: 10.1136/bjsports-2015-094758

Schäfer, D., Gjerdalen, G. F., Solberg, E. E., Khokhlova, M., Badtieva, V., Herzig, D., et al. (2015). Sex differences in heart rate variability: a longitudinal study in international elite crosscountry skiers. *Eur. J. Appl. Physiol.* 115, 2107–2114. doi: 10.1007/s00421-015-3190-0

Schmitt, L., Regnard, J., Parmentier, A. L., Mauny, F., Mourot, L., Coulmy, N., et al. (2015). Typology of "Fatigue" by heart rate variability analysis in elite nordic-skiers. *Int. J. Sports Med.* 36, 999–1007. doi: 10.1055/s-0035-1548885

Stanley, J., Peake, J. M., and Buchheit, M. (2013). Cardiac parasympathetic reactivation following exercise: implications for training prescription. *Sports Med.* 43, 1259–1277. doi: 10.1007/s40279-013-0083-4

Taylor, A. A. (1994). Autonomic control of cardiovascular function: clinical evaluation in health and disease. *J. Clin. Pharmacol.* 34, 363–374. doi: 10.1002/j.1552-4604.1994.tb04976.x

Tse, H. F., Siu, C. W., Tsang, V., Yu, C., Park, E., Bornzin, G. A., et al. (2005). Blood pressure response to transition from supine to standing posture using an orthostatic response algorithm. *Pacing Clin. Electrophysiol.* 28, S242–S245. doi: 10.1111/j.1540-8159.2005.00054.x

Uusitalo, A. L., Uusitalo, A. J., and Rusko, H. K. (2000). Heart rate and blood pressure vari- ability during heavy training and overtraining in the female athlete. *Int. J. Sports Med.* 21, 45–53. doi: 10.1055/s-2000-8853

White, N. J. (1980). Heart-rate changes on standing in elderly patients with orthostatic hypotension. *Clin. Sci.* 58, 411–413. doi: 10.1042/cs0580411

**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.

Copyright © 2020 Barrero, Le Cunuder, Carrault, Carré, Schnell and Le Douairon Lahaye. 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.