The physiological, perceptual and neuromuscular responses of team sport athletes to a running and cycling high intensity interval training session

Purpose The acute physiological, perceptual and neuromuscular responses to volume-matched running and cycling high intensity interval training (HIIT) were studied in team sport athletes. Methods In a randomized cross-over design, 11 male team sport players completed 3 × 6 min (with 5 min between sets) repeated efforts of 15 s exercising at 120% speed (s\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\dot{\text{V}}$$\end{document}V˙O2max) or power (p\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\dot{\text{V}}$$\end{document}V˙O2max) at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\dot{\text{V}}$$\end{document}V˙O2max followed by 15 s passive recovery on a treadmill or cycle ergometer, respectively. Results Absolute mean \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\dot{\text{V}}$$\end{document}V˙O2 (ES [95% CI] = 1.46 [0.47–2.34], p < 0.001) and heart rate (ES [95% CI] = 1.53 [0.53–2.41], p = 0.001) were higher in running than cycling HIIT. Total time at > 90% \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\dot{\text{V}}$$\end{document}V˙O2max during the HIIT was higher for running compared to cycling (ES [95% CI] = 1.21 [0.26–2.07], p = 0.015). Overall differential RPE (dRPE) (ES [95% CI] = 0.55 [− 0.32–1.38], p = 0.094) and legs dRPE (ES [95% CI] = − 0.65 [− 1.48–0.23], p = 0.111) were similar, whereas breathing dRPE (ES [95% CI] = 1.01 [0.08–1.85], p = 0.012) was higher for running. Maximal isometric knee extension force was unchanged after running (ES [95% CI] = − 0.04 [− 0.80–0.8], p = 0.726) compared to a moderate reduction after cycling (ES [95% CI] = − 1.17 [− 2.02–0.22], p = 0.001). Conclusion Cycling HIIT in team sport athletes is unlikely to meet the requirements for improving run-specific metabolic adaptation but might offer a greater lower limb neuromuscular load.


Introduction
As an appropriate stimulus for improving V O 2max , high intensity interval training (HIIT) offers a training approach that is easily accommodated into an individual's training schedule to acutely or progressively overload the cardiopulmonary, metabolic and neuromuscular systems (Buchheit and Laursen 2013a, b;Dolci et al. 2020). Despite many studies examining the acute and chronic adaptive responses to HIIT in athletes (e.g., Dupont et al. 2002;Buchheit et al. 2009;Wong et al. 2010;Jones et al. 2015;Beard et al. 2019), little is known about how the selected training modality (i.e., running vs. cycling) influences an individual's response to HIIT. In many team sports, practitioners and rehabilitation staff will choose training practices that simultaneously reduce 1 3 musculoskeletal load while promoting appropriate central and peripheral stimuli. The use of cycle-based HIIT has been used with team sport athletes over periods of 2-6 weeks to improve intermittent running and cycling performance (Jones et al. 2015;Hamlin et al. 2017;Beard et al. 2019; Thom et al. 2020), while others have reported no effect of 5 weeks of cycle-based HIIT on running performance (Goods et al. 2015). Cycling might be adopted for athletes who require more careful load management, e.g., during rehabilitation after injury or in the days after match play, or to provide alternative training stimuli (Mallol et al. 2020; Thom et al. 2020). However, the application of cycle-based training with team sport athletes seems to have been applied without any direct comparison of the physiological and neuromuscular responses when compared to that from similar running-based HIIT approaches.
Running and cycling at the same relative intensity evoke distinct physiological and perceptual responses that suggest different stimuli could be applied when adopted by the same athlete (Carter et al. 2000;Hill et al. 2003;Millet et al. 2009;Mclaren et al. 2016). For example, compared to cycling, running is performed at a higher absolute metabolic state with a much faster oxygen uptake kinetic response (Carter et al. 2000;Hill et al. 2003). Exercise modality also alters localized perceived effort, with running eliciting higher central and cycling higher peripheral sensations (Mclaren et al. 2016;Rampinini et al. 2016). Differences in loading have also been used to suggest task-dependent neuromuscular responses for HIIT protocols using repeated sprint running and cycling (Rampinini et al. 2016;Tomazin et al. 2017). However, the extent to which these muscle responses are observed after lower intensity HIIT remains unclear. Describing how those athletes who regularly perform running respond to cycling activity of the same format would offer valuable insight for those seeking to better understand the application of HIIT in team sport athletes. Therefore, using a perimaximal HIIT format (~ 120% V O 2max ) that has been reported previously (15 s active: 15 s rest; Dupont et al. 2002;Buchheit et al. 2009), we sought to compare the acute metabolic, cardiovascular, perceptual and neuromuscular responses of team sport athletes to volume-matched running and cycling HIIT sessions. It was hypothesized that the relative metabolic and cardiovascular demand would be higher for running HIIT compared to cycling, despite a smaller reduction in muscle force. It was also hypothesized that perceptual responses using the differential RPE would be sensitive to the physiological and neuromuscular inputs of the exercise mode.

Methods
With Department of Sport and Exercise Sciences ethics approval, 11 male university standard team sport players (age 20.0 ± 0.8 y, stature 181 ± 5 cm, body mass 82.3 ± 12.4 kg) participated in this study after providing written informed consent. Participants represented a range of team sports, including soccer, rugby and basketball. An a priori sample size calculation using G*Power 3.1.9.6 (Faul et al. 2007) informed our sample recruitment. A sample size of nine was estimated to detect a one-tailed effect of d = 0.93 with a power of 80% and error rate of 5% using a paired-samples t test. The effect size of interest (d = 0.93) was estimated based on differences in maximal oxygen uptake for running and cycling in trained students (McArdle and Magel 1970), with the one-tailed option selected because of the systematically higher V O 2 values observed in running compared to cycling (Millet et al, 2009). Our choice of power and error rate were based on common practice (Lakens 2022), but were ultimately arbitrary. Twelve participants were initially recruited to account for participant drop out. Additional participants began testing before 9 complete sets of data had been collected-it was decided it was appropriate to complete the data collection for these additional participants. Participants first attended the laboratory on two separate occasions (temperature: 20.5 ± 1.1 compared to (cf.) 20.3 ± 1.1 °C, p = 0.518; humidity: 52.8 ± 3.3 cf. 52.1 ± 4.0%, p = 0.623; Pbar: 763 ± 11 cf. 764 ± 11 mmHg, p = 0.902; all ES < 0.2) completing incremental tests to exhaustion to establish speed (sVO 2max ) and power (pVO 2max ) corresponding to maximal oxygen uptake ( V O 2max ) during running (H/P Cosmos, Pulsar, Nussdorf-Traunstein, Germany) and cycling (Lode Excalibur Sport, Lode Medical Technology, Groningen, The Netherlands), respectively. The protocols started at 100 W (cycling) or 8 km h −1 with a 1% incline (running) and increased by 20 W min −1 (cycling) and speed by 0.5 km h −1 min −1 (running) until volitional exhaustion. Volitional exhaustion was defined as either the point at which participants could no longer maintain a cycling cadence of 50 rev min −1 , or the speed of the treadmill. Expired air was collected continuously throughout each exhaustive trial using a pre-calibrated metabolic cart (Quark RMR, Cosmed, Cosmed.S.R.L., Italy). Oxygen uptake ( V O 2 ), was recorded breath-by-breath and later averaged over 30 s, with heart rate (HR) collected via telemetry (Garmin Premium HR, Garmin Ltd, Kansas, USA). V O 2max was accepted as the highest V O 2max averaged over 30 s.
Participants completed two HIIT trials using either running (H/P Cosmos, Pulsar, Nussdorf-Traunstein, Germany) or cycling (Lode Excalibur Sport, Lode Medical Technology, Groningen, The Netherlands) in a randomized crossover design, with 5-7 days between trials (temperature: 20.6 ± 0.8 cf. 20.4 ± 0.7 ℃, p = 0.465; humidity 53.3 ± 5.0 cf. 54.7 ± 4.7%, p = 0.606; Pbar 763 ± 1.7 cf. 763 ± 4.4 mmHg, p = 0.872; all ES < 0.2). Each HIIT session comprised 15 s at 120% s V O 2max (running; 15.9 ± 1.7 km h −1 ) or p V O 2max (cycling; 301 ± 28 W) followed by 15 s passive recovery, repeated for 6 min. Participants completed 3 sets with a 5 min passive recovery between each 6 min bout. Passive recovery during running was achieved by the participants placing their hands on the handrails and straddling the treadmill belt, while during cycling, the participants remained seated and legs were stationary. During treadmill running, participants wore a safety harness. Oxygen uptake and heart rate were measured throughout, with values for mean V O 2 (absolute and relative to mode-specific maximum), time > 90% maximum values and energy expenditure (kcal min −1 ) calculated (Weir 1949). Blood lactate concentration (Lactate Pro II, Arkray, Japan) was recorded immediately after with differential rating of perceived exertion (dRPE) for overall exertion (dRPE-O), breathlessness (dRPE-B) and leg-muscle exertion (dRPE-L) recorded 30 min after each HIIT trial using the centiMax scale (CR100; Borg and Borg 2002). Maximal voluntary isometric contraction of the knee extensors (MVC) in the dominant limb was measured immediately before and after each HIIT trial (S Beam Load Cell, Richmond Industries, Reading, UK) with the participant seated and the knee angle fixed at 90 degrees.

Statistical analysis
All comparisons are reported as effect sizes (Cohen's d; mean difference between trials/pooled standard deviation) and 95% confidence intervals (ES [95% CI]), with threshold values of 0.0-0.2, trivial; 0.21-0.6, small; 0.61-1.2, moderate; 1.21-2.0, large; > 2.0, very large. These arbitrary thresholds were used in the absence of accepted minimum thresholds for changes in the measurements of interest. Effects with confidence intervals that crossed a small positive or negative change were classified as unclear. For those wishing to interpret the analysis using a more traditional approach, we provide p values based on appropriate null hypothesis tests, although any ES confidence interval that includes zero can be considered as p > 0.05. Data were checked for assumptions of normality using the Shapiro-Wilk test and were found to be normally distributed (p > 0.05). Differences in physiological and perceptual responses were analyzed using separate paired-samples t tests, whereas differences in time spent > 90%VO 2max and % V O 2max during each bout, and changes in MVC were examined using separate repeatedmeasures analysis of variance. All data were analyzed using SPSS (version 27, Chicago, Illinois, USA) or a custom-made spreadsheet (https:// www. cem. org/ effect-size-calcu lator).   (Fig. 2).

Discussion
Running elicited a higher V O 2max and HRmax for team sport athletes than cycling. A ~ 10% higher V O 2max in running compared to cycling is attributed to the greater active muscle mass and the mode-specific adaptations that occur when running is a large proportion of the participant's habitual training (Hill et al. 2003;Millet et al. 2009). While differences in HR max are of a smaller magnitude (~ 2%) and less certain (i.e., with confidence intervals spanning a small decrease and increase), our observations are consistent with previous studies that report slightly lower HR values in cycling compared to running exercise (Roecker et al. 2003;Millet et al. 2009).
Given the greater capacity and requirement for oxygen consumption during running, absolute mean V O 2 during the HIIT running session was understandably higher than cycling. However, when mean V O 2 during the HIIT session was expressed relative to mode-specific V O 2max , the difference between exercise modes was less clear-we observed a moderate difference, but our data were equally compatible with a large increase in running to a decrease with cycling. The relatively low overall oxygen demand for running (~ 76% V O 2max ) and cycling (~ 70% V O 2max ) during the 3 × 6 min bouts was noticeable, and similar to that reported before for studies adopting HIIT sessions comprising 15 s  (Rozenek et al. 2007). Rozenek and colleagues (2007) had participants run at intensities lower than those used in our study (i.e., 100% s V O 2max ) but also incorporated an active recovery at 50% s V O 2max . The use of a passive recovery interval is likely to have contributed to the low mean % V O 2max observed in our study (Buchheit and Laursen 2013a) and is therefore an important consideration when planning HIIT sessions, more so when using stationary cycling as the mode of exercise. The proportion of V O 2max used during the running is also lower than that reported by Buchheit et al. (2009), who showed that during intermittent shuttle runs comprising 15 s work at 120% s V O 2max : 15 s passive recovery, mean V O 2 was ~ 88% of peak value. The use of treadmill running in the current study, where participants jumped on and off the treadmill moving at the set speed, is in contrast to the shuttle running on an indoor surface used by Buchheit et al. (2009). Participants on the treadmill would not be subject to the accelerations and decelerations inherent in shuttle running that would increase the metabolic cost of exercise (Stevens et al. 2015). The mode of how running is performed, i.e., outdoor cf. treadmill running, is therefore an important consideration for practitioners. The mean time at > 90% sVO 2max during running in our study (288 ± 132 s) was similar to values reported by Dupont et al. (2002) using the same running protocol (323 ± 272 s), albeit the within group variability was lower in our study. Differences are probably explained by Dupont and colleagues' participants running over ground from a stationary start to cover a fixed distance in 15 s, meaning the need to accelerate and time at the required running speed would have fluctuated more, compared to our participants who ran on a treadmill at the fixed speed during the 15 s (Dupont et al. 2002).
Time spent at > 90% V O 2max is a key parameter for adaptation to HIIT, with target times of ~ 5-7 min of total exercise time proposed for team sport athletes (Buchheit and Laursen 2013b;Paquette et al. 2019;Dolci et al. 2020). Running elicited a greater time > 90% V O 2max than cycling HIIT, which equated to more than twice the total exercise time (27 ± 12% cf. 12 ± 12% total training session for running and cycling, respectively). Our confidence intervals for this observation provide some certainty that the effect is to increase time > 90%VO 2max , albeit this effect could range from small to very large. More time above the pre-defined threshold for adaptation in running compared to cycling might be explained by running possessing a higher metabolic demand that leads to a larger and faster primary phase of the V O 2 response (Hill et al. 2003). Millet et al. (2003) also reported a positive association between the time constant (τ) of the primary phase of the V O 2 kinetics and time above > 90% V O 2max . Faster oxygen uptake kinetics during running would mean that during the repeated 15 s efforts, the oxygen demands were met much sooner than in cycling and over the 18 min exercise period running elicited an improved opportunity to increase V O 2 above the threshold. Notably, some participants recorded no or very limited time above the defined threshold during cycling (Fig. 1a). The use of short intervals (i.e., ≤ 15 s) using cycling in team sport athletes might therefore offer little value in targeting adaptations in V O 2max . Future studies exploring the responses to longer intervals and other HIIT types (e.g., sprint interval, repeated sprints) during off-feet training are needed.
While understanding the V O 2 response is an important measure of determining adaptation to HIIT, other measures of exercise intensity are typically used in practice. Mean HR was higher for running compared to cycling (Millet et al. 2009) and the time spent at > 90%HRmax was consistent with a moderate-to-very large difference between exercise modes. Albeit our selection of 90% of maximum values is arbitrary and might not represent the same physiological intensity (Achten and Jeukendrup 2003), the data suggest that in short intervals, such as those used here, heart rate is slower to respond in cycling compared to V O 2 (Midgley et al. 2007). This 'lag' in heart rate during short intervals reaffirms the challenges of using this measure to monitor exercise intensity and load during short-duration HIIT in team sport athletes (Buchheit et al. 2009).
The higher reported breathlessness (dRPE-B) for running probably reflected the higher metabolic demand and central responses (e.g., oxygen uptake, breathing rate, etc.) of this exercise modality, albeit our data were consistent with this effect being trivial to very large. Differences in favour of a higher rating of leg-exertion (dRPE-L) for cycling were more certain (small to large effects), and accompanied a greater reduction in MVC after this mode of exercise with similar certainty (small to large). Our use of dRPE therefore offered a potentially sensitive measure capable of differentiating between the specific central and peripheral inputs during HIIT (Mclaren et al. 2016).
Understanding the neuromuscular response to short-duration HIIT is necessary because of the potential impact on subsequent training sessions (Leveritt and Abernethy 1999) and that very little data examining muscle force after HIIT exercise exists. The cycling HIIT session caused moderate reductions (~ 16%) in MVC that were not observed after running HIIT (~ 1%). Indeed, reductions in MVC after running were unclear, in that our data were consistent with a moderate reduction and a moderate increase. If our observation of a trivial difference is correct, there are several potential mechanisms that could work either in isolation or combination to explain a greater force loss in cycling. A greater eccentric loading would be anticipated in running compared to cycling due to activation of the stretch-shortening cycle. Therefore, in running, force production would possibly have been enhanced for a given neural input (de Haan et al. 1991) delaying the onset of peripheral fatigue with a lower recruitment of type II motor units during running compared with cycling HIIT. Muscle activation across a number of muscle groups (i.e., knee extensors, knee flexors and plantar flexors) is likely to have occurred in running, whereas concentric actions of the knee extensors would predominate in cycling leading to a lower efficiency (Bijker et al. 2002). A greater contribution of the upper body musculature to overall V O 2 during running means the metabolic cost of upper body exercise during cycling makes a smaller contribution to the total exercise V O 2 . Therefore, in cycling HIIT, the lower body is more likely to be closer to its individual maximal oxygen consumption and its maximal voluntary contraction (Carter et al. 2000). This might require a progressive recruitment of the less efficient type II muscle fibres as the initially recruited type I fibres become fatigued, particularly given the exercising muscle is the principal origin of the slow component (Burnley and Jones 2007). During heavy cycling, there is high intramuscular tension and the recruitment of type II motor units is closely related to the requirements for muscle force generation (Carter et al. 2000). Higher intramuscular pressures might also cause partial occlusion of femoral arterial blood flow, reducing oxygen delivery that increases type II motor unit recruitment (Carter et al. 2000). While mechanisms remain speculative, we identify that cycling HIIT in team sport players resulted in a greater decrease in muscle force immediately after exercise than did running. This has clear implications for team sports that might use cycling as part of concurrent training practices. Further work exploring the recovery of muscle force after cycling HIIT in team sport athletes and impact on subsequent training is needed.
The study is not without limitations. First, the use of male university standard team sport athletes means that our findings might not translate directly to those athletes of a higher or lower standard or to female participants. Given the task-dependent nature of fatigue (Enoka and Duchateau 2008), we were unable to ascertain the central and peripheral components that contributed to the changes muscle force after cycling and if these differed to running. We also did not establish the cellular stress and molecular responses to HIIT for cycling and running. Finally, we acknowledge that our approach to sample size estimation has several limitations. Our anticipated effect size was estimated from a single study, which might have overestimated the population effect size. We selected a one-tailed test owing the typically higher V O 2 associated with running, but accept that we could not be certain of the direction of all effects, and therefore, a two-tailed option might have been more appropriate. Studies with sufficient resources might also wish to base their sample on a greater power than the commonly used 80% that we adopted. Indeed, a larger and more powerful sample would likely have improved the precision of our population estimates, which often had wide confidence intervals.

Conclusion
This study examines the responses to a specific HIIT training session using cycle ergometry and running in team sport athletes, offering valuable insight to those team sport practitioners using cycle-based training with their athletes. These data highlight cycling elicited lower responses compared to running during short-duration (15 s) high intensity interval training. The time above the threshold for adaptation (i.e., > 90% V O 2max ) is less likely to be met when using off-feet cycle training in team sport athletes. Short-interval cycling might offer a greater lower limb neuromuscular load without the need for accelerations, decelerations and changes of direction that are observed in running. When using cycling-based HIIT in team sport athletes where a high fractional utilization of maximum oxygen uptake is required, practitioners should also consider a prior warm-up that speeds the oxygen response (e.g., Jones et al. 2003) and the appropriateness of the work and rest interval intensity to elicit the required physiological stimulus.
Author contributions CT and JH conceived and designed research. RB, CT and JH conducted experiments. CT and JH analyzed data. CT wrote the manuscript. All authors read and approved the manuscript.
Funding This project did not receive any funding.

Conflict of Interest
The authors declare that they have no conflict of interest.
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:// creat iveco mmons. org/ licen ses/ by/4. 0/.