Bioenergetics of the VO2 slow component between exercise intensity domains

During heavy and severe constant-load exercise, VO2 displays a slow component (VO2sc) typically interpreted as a loss of efficiency of locomotion. In the ongoing debate on the underpinnings of the VO2sc, recent studies suggested that VO2sc could be attributed to a prolonged shift in energetic sources rather than loss of efficiency. We tested the hypothesis that the total cost of cycling, accounting for aerobic and anaerobic energy sources, is affected by time during metabolic transitions in different intensity domains. Eight active men performed 3 constant load trials of 3, 6, and 9 min in the moderate, heavy, and severe domains (i.e., respectively below, between, and above the two ventilatory thresholds). VO2, VO2 of ventilation and lactate accumulation ([La−]) were quantified to calculate the adjusted oxygen cost of exercise (AdjO2Eq, i.e., measured VO2 − VO2 of ventilation + VO2 equivalent of [La−]) for the 0–3, 3–6, and 6–9 time segments at each intensity, and compared by a two-way RM-ANOVA (time × intensity). After the transient phase, AdjO2Eq was unaffected by time in moderate (ml*3 min−1 at 0–3, 0–6, 0–9 min: 2126 ± 939 < 2687 ± 1036, 2731 ± 1035) and heavy (4278 ± 1074 < 5121 ± 1268, 5225 ± 1123) while a significant effect of time was detected in the severe only (5863 ± 1413 < 7061 ± 1516 < 7372 ± 1443). The emergence of the VO2sc was explained by a prolonged shift between aerobic and anaerobic energy sources in heavy (VO2 − VO2 of ventilation: ml*3 min−1 at 0–3, 0–6, 0–9 min: 3769 ± 1128 < 4938 ± 1256, 5091 ± 1123, [La−]: 452 ± 254 < 128 ± 169, 79 ± 135), while a prolonged metabolic shift and a true loss of efficiency explained the emergence of the VO2sc in severe. Electronic supplementary material The online version of this article (10.1007/s00424-020-02437-7) contains supplementary material, which is available to authorized users.


Introduction
After the cardio-dynamic phase, oxygen consumption (VO 2 ) during constant heavy and severe intensity exercise, is better fitted by a two-component rather than a single-component model. Based on this observation, it has been suggested, that a slow component of VO 2 exists (VO 2sc ), that does not start at exercise onset but rather appears later in exercise (time delaỹ 120-180 s) [9]. Furthermore, it has been assumed, but not proven, that in association with VO 2sc , the muscle displays an increasing energy demand as a function of time [10]. VO 2sc represents an increased O 2 cost of locomotion when exercise is protracted more than 3 min at a constant workload above the lactate threshold (LT) [9]. Typically, when exercise is performed between LT and the critical power, (i.e. heavy intensity domain) VO 2sc tends to a delayed steady-state (i.e. compared to moderate). On the contrary, when effort rises above critical power (i.e. severe exercise domain) a steadystate is not achievable and VO 2 increases tending to the maximal oxygen consumption (VO 2max ) [9]. The magnitude of VO 2sc is considered to be associated with exercise intolerance and fatigue [7], and the study of this phenomenon is particularly important in diseased/frail populations (e.g. [19,24]). Nevertheless, despite the efforts of the past 40 years to clarify the exact physiological underpinnings of the VO 2sc , these remain debated even in healthy subjects. [9].
A key finding from the past studies was that roughly 85% of the VO 2sc originates from the contracting muscles, while the remaining 15% corresponds to the increased VO 2 cost of Electronic supplementary material The online version of this article (https://doi.org/10.1007/s00424-020-02437-7) contains supplementary material, which is available to authorized users. ventilation [20]. Focusing on the muscular component of the VO 2sc , subsequent investigations proposed that the recruitment of less efficient type II fibres necessary to maintain a specific power output (PO) [9] or metabolic instability due to metabolites accumulation occurring within the working fibres [26,28] could elicit an increased cost of locomotion. Nevertheless, a satisfactory theory explaining the mechanisms underpinning the VO 2sc is still missing.
In this context, a recent study [17] refuted that the energy demand of a constant, high-intensity exercise changes over time. By subtracting the VO 2 cost of ventilation and accounting for the contribution of the glycolytic energy sources [17,21,22] O'Connell et al. quantified the total, adjusted metabolic cost of muscle exercise over time during a constant load trial in the severe domain. The authors concluded that the overall cost of locomotion does not increase over time, other than what required by the augmenting cost of ventilation. The apparent VO 2 increases over time would be the result of the prolonged shift in metabolic sources beyond the first 3 min of exercise (i.e. an increased contribution of the aerobic metabolism to ATP resynthesis, mirrored by a decreased contribution of anaerobic ATP resynthesis over time) . In other words, these findings suggest that VO 2sc may in fact not represent a true loss of efficiency as a function of time but rather a prolonged adjustment of the oxidative metabolism. However, while this new hypothesis could explain why a satisfying explanatory theory of VO 2sc is still missing, the amount of literature that reported VO 2sc in different exercise modalities (i.e. whole-body vs isolated muscles) and models (i.e. humans vs animals) calls for further investigation [9]. In particular, the intensity domain in which O'Connell measures were obtained is not fully clear; while most studies in the field of VO 2 kinetics define intensity domains based on gas exchange thresholds, O'Connell used lactate-based measures. Yet, the correspondence between the above methods remains a controversial issue [3,12,13]. In addition, loss of efficiency, the metabolic shift among substrates and excess ventilation manifest and increase differently over time in the three domains of exercise (e.g. no loss of efficiency is present in the moderate domain while an elevated steady state is reached in the heavy yet not in the severe domain) [10]. Therefore, O'Connell's findings need to be confirmed and broadened with particular consideration to the three intensity domains of exercise.
Accordingly, by calculating the energy cost of ventilation, the glycolytic contribution to exercise, and directly measuring the aerobic cost of locomotion over time, we tested if and to what extent a true loss of efficiency during cycling explains the emergence of the slow component of VO 2 in different intensity domains. Specifically, we hypothesised that (i) the overall cost of locomotion would not be affected by the time during metabolic transitions in the moderate and heavy exercise domains (i.e. no loss of efficiency over time exists at these intensities but rather a prolonged metabolic shift is responsible for the observed VO 2sc ); (ii) the overall cost of locomotion would be affected by the time during metabolic transitions in the severe exercise domain (i.e. a prolonged metabolic shift is not sufficient to explain the emerge loss the observed VO 2sc and a true loss of efficiency over time will manifest at this intensity).

Ethical approval
The study was conducted according to the Declaration of Helsinki and all procedures were approved by the University of Verona Ethics Committee for Research on Human Subjects (CARU no 15/2019). Procedures and risks were explained to each subject, and all participants volunteered and gave informed written consent to participate before the start of the study.

Participants
Eight active men were recruited in the study (age 25 ± 2 years, body mass 74 ± 10 kg, height 181 ± 5 cm, VO 2max 49.3 ± 3.4 ml*kg −1 *min −1 , performing recreational aerobic activity 3-4 times × week). Inclusion criteria were male sex and age between 20 and 35 years; exclusion criteria were smoking and any condition that could influence the physiological responses during testing. All participants were instructed to avoid caffeine consumption and physical activity respectively for at least 8 h and 24 h before each testing session. Moreover, they received a standard and individualized food intake prescription before all the testing sessions to minimize variability of glycogen stores and glucose oxidation (i.e. 2 g of low glycemic index carbohydrates per kg of body weight, 2 h before testing; 0.5 L of water in the 90 min before testing; restriction from caffeine during the 8 h before testing).

Experimental protocol
After medical clearance, participants visited the laboratory on ten occasions within a maximum of 4 weeks. All subjects completed: (i) a preliminary maximal ramp incremental exercise test to exhaustion for the determination of gas exchange threshold (GET), respiratory compensation point (RCP) and the peak of oxygen uptake (VO 2peak ); (ii) three constant load trials (CLT) respectively of 3, 6 and 9 min in the "moderate" exercise intensity domain (iii) three CLT respectively of 3, 6 and 9 min in the "heavy" exercise intensity domain (iv) three CLT trials respectively of 3, 6 and 9 min in the "severe" exercise intensity domain. Tests of (ii), (iii) and (iv) aimed at determining VO 2 response and blood lactate ([La − ]) accumulation as a function of time in the three exercise intensity domains in order to evaluate the overall energetic contribution of the aerobic and anaerobic metabolisms [22] to ATP turnover. Moreover, these tests were executed in randomized order with the only exception of the longest CLT in the "severe" exercise domain, that was completed as first to assure that subjects were able to sustain the PO for the required time. All exercise tests were conducted on an electromagnetically braked cycle ergometer (Sport Excalibur, Lode, Groningen, Netherlands), at a similar time of the day in an environmentally controlled laboratory (18°C, 55-65% relative humidity).

Ramp incremental test
The ramp incremental test consisted of a 3-min baseline cycling at 50 W, followed by a 30-W*min −1 increase in PO until volitional exhaustion. Participants were asked to pick a selfselected cadence in the range of 70-90 rpm and to maintain it throughout all subsequent tests. Failure to maintain the indicated cadence within 5 rpm (for longer than 5 s) during testing despite strong verbal encouragement was considered as the criterion for exhaustion. Breath-by-breath pulmonary gas exchange and ventilation were continuously measured using a metabolic cart (Jaeger Oxycon Pro, Viasys Healtcare GmbH, Höchberg, Germany) as previously described [23]. Heart rate (HR) was monitored continuously (H7 Sensor, Polar, Kempele, Finland).

Constant load trials
After the preliminary ramp incremental test, subjects completed 3 CLT within each exercise intensity domain (i.e. moderate, heavy, and severe [18]) in a randomized order:
Each CLT was preceded by a 3-min warm-up at 20 W. Throughout the test, subjects kept the same, constant rpm and bike position as selected during the ramp incremental test. VO 2 and HR data were measured with the same method described for the ramp incremental test. Moreover, capillary blood samples (65 μl) were drawn from the fingertip in the last 30 s of warm-up and at the 1st, 3rd, 5th, and 7th min after each test and were immediately analysed (Radiometer ABL90 FLEX, Radiometer Medical ApS, Brønshøj, Denmark) to measure [La − ]. The highest value was considered as the peak of blood lactate concentration and used for further analysis to calculate anaerobic energetic contribution.

Ramp incremental test
For the gas exchange variables, aberrant data-points that lay 3 SD from the local mean were removed, and trials were linearly interpolated on a 1-s basis and then averaged every 5 s. VO 2peak was determined as the highest VO 2 obtained over a 10-s interval [6]. GET and RCP were determined with the standard technique from gas exchange variables by three blinded expert reviewers as detailed elsewhere [6]. Briefly, GET was determined by visual inspection as the VO 2 at which CO 2 output began to increase out of proportion in relation to VO 2 , with a systematic rise in the ventilation (VE)-to-VO 2 relation and end-tidal PO 2 whereas the ventilatory equivalent of VCO 2 (VE/VCO 2 ) and end-tidal PCO 2 is stable [2]. RCP was determined as the point where end-tidal PCO 2 began to fall after a period of isocapnic buffering [27]. This point was confirmed by examining VE/VCO 2 plotted against VO 2 and by identifying the sec breakpoint in the VE-to-VO 2 relation. Finally, we determined the constant workload equivalent to the specific moderate (80% of GET), heavy (50% Δ between GET and RCP), and severe (60%Δ between GET and VO 2peak ) VO 2 targets. To this aim, the VO 2 /W relationship identified with the incremental test was left-shifted to account for the mean response time [6].

Constant load trials
VO 2 , VCO 2 , and VE during CLT were sampled breath-bybreath and interpolated using the same procedure described for the ramp incremental test. Interpolated data from different CLT performed at the same exercise intensity were mediated in order to reduce breath-by-breath signals variability (data from 3 tests were mediated in the time segment from 0 to 3 min and data from 2 tests were mediated between 3 and 6 min) [11]. The sum of the oxygen consumed during each 3-min time segment was then considered as the aerobic energetic contribution to exercise. VO 2sc was calculated as the sum of the amount of oxygen exciding the VO 2 reached at the end of the 0-to-3-min time segment [25].
Work of breathing (WB) was calculated based on VE using the equation by Coast The so obtained values were utilized to estimate the energy contribution from anaerobic glycolysis in each time segment, based on the oxygen equivalent for lactate of Di Prampero (i.e. 1 mmol*L −1 [La − ] accumulation = 3.0 ml*kg −1 VO 2 ) [22].
Overall energetic cost of the activity (expressed as ml of oxygen) was calculated as described by O'Connell et al. and defined as "adjusted oxygen equivalent" (AdjO 2Eq ): o.
the amount of oxygen equivalent utilized to sustain each W during cycling) was calculated as the ratio between the amount of AdjO 2Eq required to sustain exercise during a specific time segment and the number of W imposed by the test:

Statistics
After assumptions verification (i.e., normality, homogeneity of variance), the within-subject coefficient of variation and a two-way repeated-measure ANOVA (trial × intensity domain) were used to evaluate VO 2 data repeatability measured at the end of the third min of exercise of each CLT.
Two-way repeated-measure ANOVA was also performed to assess differences over time between different intensities domains (time segment × intensity domain) for VO 2 , VO 2VM , O 2 equivalent of [La − ], AdjO 2Eq , VO 2gain . The amplitude of the VO 2sc in different intensity domains was compared using a one-way repeated measures ANOVA (intensity domain). Post-hoc analyses were performed using the Holm-Sidak method. Moreover, partial eta squares (η p 2 ) were calculated to quantify the effects sizes of different independent variables during the constant load trials [14]. Based on an expected standard deviation of breath-by-breath VO 2 measurements for steady-state exercise equal to 2.5%, and a minimum detectable change in VO 2 of 100-170 ml·min −1 at a VO 2 of 2.1 to 3.5 L·min −1 [12], the minimum sample size to obtain a power of 0.8 was 6 individuals.
Data are presented as means ± SD. All statistical analyses were performed using Sigmaplot version 12 and α was set in advance at the 0.05 level. Statistical significance was accepted when p < α.

Results
Subjects' anthropometrical and functional characteristics obtained during the ramp incremental test are reported in Table 1.
An overview of the energetic contributors to exercise (i.e., measured VO 2 , VO 2VM , [La − ] equivalent of O 2 , and AdjO 2Eq ) is reported as 10-s averages and over 3-min time segments in Fig. 2 Table 3. For all time segments, the contribution of anaerobic glycolysis (as represented by [La − ] equivalent) was significantly increased going from moderate to heavy to severe intensity trails (significant main effect of intensity, p ≤ 0.001, η p 2 : 0.91). Furthermore, the contribution of anaerobic glycolysis in the three domains was significantly affected by time (significant main effect of time segment p ≤ 0.001, η p 2 : 0.67). In particular, the post-hoc analysis revealed a limited and unvaried contribution of anaerobic glycolysis for the moderate-intensity domain (0-3 vs 3-6 time segments p = 0.772, 3-6 vs 6-9 p = 0.947); in heavy, a decreased contribution of anaerobic glycolysis from 0 to 3 to 3-6 time segments, yet no further after that (0-3 vs 3-6 p = 0.014, 3-6 vs 6-9 p = 0.429); finally, a continuous decrease in the contribution of anaerobic glycolysis over time in the severe exercise trials (0-3 vs 3-6 p ≤ 0.001, 3-6 vs 6-9 p = 0.024; Fig. 3).

Discussion
By comprehensively quantifying aerobic and glycolytic energy sources used during exercise, deprived of the VO 2 cost of ventilation, this study tested the hypothesis that the overall cost of cycling (i.e. AdjO 2Eq ) is affected by the time during metabolic transitions in the moderate, heavy and severe intensity domains. In agreement with our hypothesis, our findings suggest that (i) the overall cost of locomotion does not increase over time in the moderate domain after the 3rd min of Fig. 1 Repeatability of VO 2 : group mean VO 2 data are displayed as 10s-averages respectively for the "severe" (top panel), "heavy" (medium panel), and "moderate" (bottom panel) exercise intensity domain. Symbols represent the three duration trials: black dots = 9-min CLT, white dots = 6-min CLT, black triangles = 3-min CLT Values are means ±SD. PPO, peak PO; VO 2peak : peak of oxygen consumption; MRT, mean response time; GET, gas exchange threshold; RCP, respiratory compensation point exercise; (ii) in the heavy intensity domain, the emergence of a VO 2 slow component is ascribable to a "metabolic shift" between aerobic and anaerobic metabolisms that protracted beyond the 3rd min of exercise and to a small yet significant increase in the VO 2 of ventilation, rather than to an actual increased cost of locomotion after the 3rd min; (iii) finally, only in the severe domain the raise of the VO 2 slow component deprived of its ventilatory part was not completely explained by a prolonged metabolic shift between energy sources, and the further increase in AdjO 2Eq was possibly indicating that a true loss of efficiency over time occurred. This investigation provides new insight of the mechanisms underpinning VO 2sc and exercise tolerance. The fitness level of our sample measured during the ramp incremental tests indicated that the participants of this investigation were active subjects. (VO 2peak : 49 ± 3 ml*min −1 *kg, see Table 1 for the other characteristics) [1], and presented an overall metabolic response to constant load exercise consistent with previous studies (e.g. no VO 2sc in the moderate domain; VO 2sc in heavy +4.2 ± 3.4% tending to steady-state; VO 2sc in severe domain +11.4 ± 4.9% tending to VO 2max ) [17,25,29].
By applying an approach previously proposed by O'Connell et al. [17], our study quantified the individual contributors to the overall net VO 2 cost of exercise (i.e. AdjO 2Eq ) and their interplay during the metabolic transition and steady-state. While O'Connell et al. studied the VO 2 of active muscles, the VO 2 associated with ventilation and the VO 2 equivalent of [La − ] in the severe domain only, we extended the analysis to the three exercise intensity domains.
As expected, based on the changes in ventilation described above, the VO 2 associated with ventilation increased going from moderate (3 ± 3% of the VO 2sc ), heavy (7 ± 3%), to severe (14 ± 7%) intensity. In agreement with previous work [20], this finding confirms that~85% of the VO 2sc occurs within the working muscles. In addition, our study demonstrated that the VO 2 associated with ventilation increased over time in the severe domain only. Furthermore, for all the exercise domains, the O 2 equivalent of [La − ] accumulation decreased after the transition phase (i.e. 0-to-3-time segment), compatible with a progressive decrease of the contribution of glycolysis to ATP resynthesis over time [22]. Interestingly, in the domains below RCP, this decreased contribution of glycolysis to the energy provision was "mirrored" by an increased contribution of oxidative metabolism, to satisfy an invariant energetic demand (Table 3 and Fig. 3). In other words, the VO 2sc that is observed in the heavy domain may be related to a prolonged "metabolic shift" between anaerobic and aerobic energy sources of ATP resynthesis rather than to an augmented cost of locomotion over time. Such a view is in agreement with the findings and interpretations of previous  human studies, performed in the heavy/severe domain of exercise [16,17] as well as with studies based on computer modelling of the skeletal muscle bioenergetics [15]. For the exercise in the severe domain, however, we found that accounting for the total energetic cost of locomotion did not completely explain VO 2sc . We were the first to demonstrate that, in this domain, the AdjO 2Eq continued to increase showing an augment energy required to sustain the same PO, compatible with a true loss of efficiency over time.
Relative to the severe intensity domain, the partial discrepancy between our data and O'Connell's may be explained by methodological differences between both studies: (i) as we did, most studies in the field of VO 2 kinetics define intensity domains based on gas exchange thresholds [10,23]; on the contrary, O'Connell used lactate-based measures. Yet, the correspondence between the above methods remains controversial [3,12,13] and the difference may have led to unmatched intensity domains among our studies. (ii) O'Connell incremental protocol for VO 2peak detection was performed after 10 min of recovery from a~20-min protocol with 3-min steps used for LT determination. This approach could be responsible for an underestimation of peak-PO and consequently of Δ60% and may not have guaranteed an exercise intensity corresponding to the severe domain for all participants. In summary, while a direct comparison between different studies may be difficult, our data provides the first, comprehensive, domain-specific characterization of the contributors to the observed VO 2sc .
Traditionally, VO 2sc has been attributed to an increased cost of locomotion when exercise is protracted more than 3 min at a constant workload in the heavy and severe exercise intensity domains, in relation to either fatigue or recruitment of higher-order motor units or both [9]. Interestingly, while a clear distinction between exercises performed above and below RCP is a common concept in exercise physiology [10], VO 2sc measured in the heavy and in the severe exercise domains is usually considered as the expression of a single phenomenon [7,9]. In this context, a recent paper based on mathematical modelling of the muscle bioenergetics, proposed that a metabolic shift between the aerobic and the anaerobic energy systems, caused by a progressive inhibition of the glycolytic ATP supply by cytosol acidification, may also contribute Fig. 2 An overview of the energetic contributors to exercise is reported in 10-s averages for the "severe" (top panel), "heavy" (medium panel), and "moderate" (bottom panel) exercise intensity domain. White columns represent directly measured VO 2 . Grey columns indicate the O 2 cost requested by ventilation. The black dashed line displays the energy provided by glycolytic sources over 3 min segments. Finally, the black solid line represents the adjusted cost of exercise accounting for both aerobic and glycolytic energy sources. Please note that during all the first 3 min segments, the contribution of immediate energy sources (O 2 and phosphocreatine) was not accounted and was probably the cause of the lower adjusted VO 2 between the first and the successive time segments [8,22] to VO 2sc [15]. The authors also suggest that the size of the VO 2sc can increase when the contribution of glycolysis to ATP resynthesis (and in turn proton accumulation) is higher; a lower VO 2 in the initial stages of exercise rather than an increased VO 2 after 3-6 min into exercise may explain the larger VO 2sc at the higher exercise intensities [15]. This is the first experimental study to determine the contributors of VO 2sc in the three intensity domains of exercise. Our findings (i.e. a prolonged metabolic shift in the heavy domain and a true loss of efficiency over time in the severe intensity domain only) appear to support the theory proposed by Korzeniewski & Zoladz in 2015. We speculate that, with increasing exercise intensity, the recruitment of bigger, preferentially glycolytic muscle fibres could explain a higher contribution of glycolysis to ATP resynthesis at exercise onset and therefore a delayed VO 2 steady state. Furthermore, the recruitment of these intrinsically less efficient and fatigable type II fibres could explain the true loss of efficiency that appears over time in the severe domain of exercise [5,7,9]. The novelty and the practical implications of this, if confirmed by further studies, would be: (i) in the heavy domain, the VO 2sc is the expression of a delayed adjustment of VO 2 rather than of a loss of efficiency developing over time; (ii) contrary to what is currently accepted, the adjustment of VO 2 in the heavy domain may be described by a slow primary component rather than by the summation of a primary plus a slow component; (iii) in the severe domain the VO 2sc may be explained by both a prolonged metabolic shift and a true loss of efficiency over time; (iv) VO 2sc likely has different physiological underpinnings in the heavy and severe domain of exercise; therefore, we should be mindful of this when the two intensity paradigms are used in the context of VO 2sc as different interventions may produce different effects in the two domains.
In conclusion, the innovative methodological approach applied in this study allowed to discriminate three contributors to the VO 2sc : increased VO 2 cost of ventilation, prolonged shifting between aerobic and glycolytic metabolism, and loss of efficiency over time. How these mechanisms contribute to the VO 2sc depends on relative exercise intensity, with a true loss of efficiency over time occurring only in the severe domain.

Limitations
It should be acknowledged that data interpretation in this investigation depends upon estimates of the energetic yield of lactate accumulation and of the VO 2 cost of ventilation and ignores the contribution of oxygen stores and anaerobic alactic mechanism of ATP resynthesis. These estimates rely on some assumptions (e.g. fingertip capillary lactate reflects wholebody lactate accumulation) and systematic biases in applying the same values/equations to different subjects. In the belief that these limitations should not have altered the interpretation of our results; we provided a deep discussion on these issues for the interested reader as supplementary material of this article. [La − ] Equivalent (ml O 2 *3min −1 ) influenced the present work. The results of the study are presented clearly, honestly and without fabrication, falsification or inappropriate data manipulation.
Authors contributions ALC and SP concepted the study methodology. ALC and KC collected the data. All the authors contributed to the analysis of the results. ALC, JB and SP drafted and revised the manuscript. All the authors approved the final version of the manuscript.

Compliance with ethical standards
Conflict of interest The authors declare that they have no conflicts of interest.
Ethical approval All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee (CARU no 15/2019) and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.
Informed consent Written informed consent was obtained from all individual participants included in the study. Fig. 3 Three-min mean oxygen equivalents of aerobic (grey, i.e. VO 2 subtracted by VO 2VM ) and glycolytic (white, i.e. [La − ] equivalent) cost of exercise are represented respectively for the "severe" (top panel), "heavy" (medium panel) and "moderate" (bottom panel) exercise intensity domain. ANOVA revealed a significant "intensity" × "time" interaction both for aerobic (p ≤ 0.001, η p 2 : 0.81) and anaerobic (p ≤ 0.001, η p 2 : 0.58) cost of exercise. Multiple comparisons are displayed as: # significant statistical difference with 0 to 3 min segment; * significant statistical difference with 3 to 6 min segment.°and § represent respectively statistical difference with the "moderate" and "heavy" exercise intensity domain. Please note that values measured during warm-up were subtracted Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.