Physiological responses during simulated 16 km recumbent handcycling time trial and determinants of performance in trained handcyclists

To characterise the physiological profiles of trained handcyclists, during recumbent handcycling, to describe the physiological responses during a 16 km time trial (TT) and to identify the determinants of this TT performance. Eleven male handcyclists performed a sub-maximal and maximal incremental exercise test in their recumbent handbike, attached to a Cyclus II ergometer. A physiological profile, including peak aerobic power output (POPeak), peak rate of oxygen uptake (V˙\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\dot{V}$$\end{document}O2Peak), aerobic lactate threshold (AeLT) and PO at 4 mmol L−1 (PO4), were determined. Participants also completed a 16 km simulated TT using the same experimental set-up. Determinants of TT performance were identified using stepwise multiple linear regression analysis. Mean values of POPeak = 252 ± 9 W, V˙\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\dot{V}$$\end{document}O2Peak = 3.30 ± 0.36 L min−1 (47.0 ± 6.8 mL kg−1 min−1), AeLT = 87 ± 13 W and PO4 = 154 ± 14 W were recorded. The TT was completed in 29:21 ± 0:59 min:s at an intensity equivalent to 69 ± 4% POPeak and 87 ± 5% V˙\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\dot{V}$$\end{document}O2Peak. POPeak (r = − 0.77, P = 0.006), PO4 (r = − 0.77, P = 0.006) and AeLT (r = − 0.68, P = 0.022) were significantly correlated with TT performance. PO4 and POPeak were identified as the best predictors of TT performance (r = 0.89, P < 0.001). POPeak, PO4 and AeLT are important physiological TT performance determinants in trained handcyclists, differentiating between superior and inferior performance, whereas V˙\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\dot{V}$$\end{document}O2peak was not. The TT took place at an intensity corresponding to 69% POPeak and 87% V˙\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\dot{V}$$\end{document}O2peak.


Introduction
Recumbent handcycling is an endurance sport for athletes with lower limb impairments, such as spinal cord injuries, lower limb amputations and congenital conditions (Abel et al. 2006). Since 2004 (Athens Paralympic Games), handcycling has been an integral part of the paracycling road programme, with 65 handcyclists contesting 13 events at the most recent 2016 Paralympic Games in Rio (Paralympics. org 2016). At national and international events (e.g. British Championships, Paralympic Games), handcyclists compete 1 3 in individual time trials (TT) (10-30 km, lasting 20-40 min) and/or road races (40-80 km, lasting 60-150 min) with many handcycling events being scheduled on consecutive days (Zeller et al. 2017). To date, most studies examining physiological responses during handcycling have gathered and described heart rate (HR), capillary blood lactate concentration (BLa), power output (PO) or rate of oxygen uptake ( V O 2 ) during laboratory-based exercise tests (Janssen et al. 2001;Abel et al. 2010;de Groot et al. 2014;Kouwijzer et al. 2018b;Quittmann et al. 2018). These studies have been limited to recreational or touring bike configurations which are disparate to that of elite handcyclists. It is only recently that recumbent handcycling studies have been conducted during competition (West et al. 2015) or protocols replicating recumbent handcycling or competitive race intensities (Fischer et al. 2015(Fischer et al. , 2020Graham-Paulson et al. 2016;Stangier et al. 2019;Stone et al. 2019b).
Physiological variables, such as peak rate of oxygen uptake ( V O 2Peak ) and ventilatory thresholds are the best predictors of handcycling TT performance. Moreover, peak power output achieved during incremental tests to exhaustion (PO Peak ) has been related to performance (Janssen et al. 2001;Lovell et al. 2012;de Groot et al. 2014;Fischer et al. 2015). However, large variance in PO Peak (129-247 W) and V O 2Peak (2.0-3.45 L min −1 equivalent to 26.5-42.3 mL kg −1 min −1 ) have been reported in the literature (Janssen et al. 2001;Abel et al. 2006;Lovell et al. 2012;Fischer et al. 2015;Graham-Paulson et al. 2018;Kouwijzer et al. 2018a;Stangier et al. 2019;Stone et al. 2019a). Confounding factors, such as training status, athlete classification and the number of years involved in the sport and indeed the protocols adopted (e.g. ramp vs. step incremental test), all contribute to these reported physiological profiles (Miller et al. 2004;Goosey-Tolfrey et al. 2008;van Drongelen et al. 2009;Lovell et al. 2012). Moreover, since the experimental designs of these studies varied considerably, using attachable-units, touring handbikes or bespoke ergometers, this makes the findings difficult to transfer to the recumbent handbikes used in the present day. These handbikes are considerably lighter, but more importantly, are bespoke to the individual, with respect to crank position, crank width and backrest inclination (Stone et al. 2018). Therefore, the purpose of this study was to characterise the physiological profiles of trained handcyclists, during recumbent handcycling and to identify which physiological variables are related to 16 km TT performance.

Experimental protocol
Participants refrained from exercise in the 24 h preceding the testing. The experiment protocol was performed on two consecutive days (Fig. 1). On the first day, participants completed a sub-maximal exercise test, to determine-the aerobic threshold (AeLT), power output (PO) at 4 mmol L −1 (PO 4 ), and maximal incremental exercise test, to assess PO Peak , V O 2Peak , and peak heart rate (HR Peak ). Ventilatory thresholds were not determined as there has been mixed results to their suitability in this sample (Kouwijzer et al. 2019;Leicht et al. 2014). Both tests were conducted in the participants' bespoke recumbent handbike (5 'Carbonbike', 4 'Top End' and 2 'Schmicking'), which were attached to a Cyclus II ergometer (RBM electronic automation GmbH, Leipzig, Germany). Following a 10-min warm-up at a selfselected cadence and intensity, the sub-maximal test commenced with an initial load of 20 W, increasing by 20 W every 4 min until BLa exceeded 4.0 mmol L −1 when the test was terminated. BLa was determined from 20 μL earlobe capillary blood samples, collected in the last minute of each stage and analysed using a Biosen C-Line (EFK Diagnostics GmbH, Germany). Breath-by-breath gas analysis (Cortex Metalyzer 3B, Cortex, Leipzig, Germany)-for calculation of V O 2 , carbon dioxide production and respiratory exchange ratio (RER)-HR (Polar RS400, Kempele, Finland) and rating of perceived exertion (RPE; Borg 1998), were also collected in the last minute of each stage. AeLT was calculated using the log-log transformation method (Beaver et al. 1985). Moreover, this test enabled the identification of the PO corresponding to a fixed BLa of 4 mmol L −1 (PO 4 ), by using linear interpolation methods as used in the handcycling literature (Stangier et al. 2019;Zeller et al. 2017).
Following 30 min' passive rest, participants performed a maximal test to exhaustion. The starting PO was equivalent to their AeLT and was maintained for 2 min. PO was then increased by 5 W every 15 s until the athlete reached volitional exhaustion (failure to maintain cadence ≥ 50 rpm and an overall RPE ≥ 19) (Graham-Paulson et al. 2016). Verbal encouragement was provided during this test. Breathby-breath gas analysis and HR were recorded continuously throughout the maximal test and BLa reported 2 min postcompletion of the maximal test. The ergometer was set in a power control mode, ensuring the pre-set PO was controlled independently of cadence or gear selection. Cadence was self-selected in all trials as reported elsewhere (Graham-Paulson et al. 2016).
On the second day, participants completed a 16-km TT in the shortest time possible. Participants selected the gear ratio to commence the TT, which could then be changed virtually by the investigators throughout the time trial, as instructed by participants. No motivation was provided during the TT and the feedback provided was PO, cadence and cumulative distance displayed on the ergometer, to maximise ecological validity. Breath-by-breath gas analysis and HR were recorded continuously throughout the TT and BLa collected every 4 km.

Data and statistical analysis
In the sub-maximal exercise test V O 2 , HR, RER and cadence were averaged across the last minute of each stage. Mechanical efficiency (ME%), calculated as the ratio of external work to energy expended (Powers et al. 1984), in 1 min of exercise, was determined at AeLT and PO 4 (RER < 1.00 for all participants). Energy expenditure was calculated from the product of V O 2 and RER and the standard conversion table (Péronnet and Massicotte 1991 Shapiro-Wilks tests were used to determine the distribution of the data. Repeated measures analysis of variance was used to determine differences in V O 2 , HR, PO and cadence across each km of the TT. Sphericity was assessed using the Mauchly's test of sphericity. If the data were aspherical, a correction, using the Greenhouse-Geisser epsilon, was applied to the calculated P value (Girden 1992). If a significant difference was identified, post hoc paired t tests, with Bonferroni corrections, were applied to determine the differences in cadence or PO.
Pearson's product-moment correlation was used to analyse the relationship between variables from the sub-maximal (AeLT, ME% at AeLT, V O 2 at AeLT, PO 4 , ME% at PO 4 and V O 2 at PO 4 ) and maximal test [ V O 2Peak (absolute and relative), HR Peak and PO Peak ] to TT performance (IBM SPSS Statistics 24, Chicago, IL). Variables that were significantly correlated with 16 km TT performance were inputted into a stepwise multiple linear regression to establish the most important determinants of 16 km TT performance. Durbin-Watson test was used to indicate the independence of variables entered into the regression. Tolerance and variance inflation factors were also used to assess multicollinearity within the data.

Determinants of TT performance
The highest correlations were observed between PO Peak (r = − 0.77, P = 0.006), PO 4 (r = − 0.77, P = 0.006), AeLT (r = − 0.68, P = 0.022) and TT performance (Table 2). Moderate to low correlations were observed for all other physiological parameters ( Table 2). The multiple linear regression indicated that PO 4 was the strongest predictor accounting for 59.3% of the variance in 16 km TT performance (P = 0.006). The addition of PO Peak (19.0% of overall variance) to PO 4 provided a stronger prediction model, accounting for 78.3% of the variance in performance (P = 0.002). A Durbin-Watson value of 1.501 indicated that variables within this model were sufficiently independent and not autocorrelated. Tolerance values of 0.741 and a variance inflation factor of 1.350 further indicate that multicollinearity was not an issue for this regression model.

Discussion
The current study investigated the physiological responses during a 16-km TT, in a population of national and international recumbent handcyclists. Extending the recent work of Stangier et al. (2019), this study focused on recumbent trained handcyclists yet was designed to measure the submaximal exercise, maximal exercise and simulated TT performance in their own bespoke recumbent handbikes. A key finding was that better 16 km TT performances were achieved by handcyclists with a higher PO Peak , PO 4 and AeLT. Conversely, V O 2Peak , both absolute and relative, were not correlated with TT performance.
The results from the maximal exercise test revealed that the handcyclists in the current study had higher average  *Significant difference to 1st km (P < 0.05), # Significant difference to 2nd km (P < 0.05), § Significant difference to 15th km (P < 0.05), ¥ Significant difference to 16th km (P < 0.05) ≥ 175 km week −1 ) to the participants from previous studies (Lovell et al. 2012;Fischer et al. 2015;Stangier et al. 2019). Furthermore, the handcyclists that achieved the greatest V O 2Peak (4.02 L min −1 and 57.6 mL kg min −1 ) and PO Peak (281 W) values averaged 5 handcycling sessions a week totalling 8 h week −1 with a self-reported distance of 170 km week −1 . Therefore, the higher physiological values observed in the current study are likely to be due to advances in endurance training regimes (Zeller et al. 2017), strength training (Nevin et al. 2018) and the use of bespoke recumbent handbikes in comparison to the modified ergometers used previously (Lovell et al. 2012;Fischer et al. 2015). Participants were likely to be more stable and comfortable in their handbike, which is perceived to be essential for the application of power in handcycling, potentially facilitating the increased V O 2Peak and PO Peak . The trained handcyclists completed the simulated 16 km TT at a similar intensity (87 ± 5% V O 2Peak ) to a 20 and 22 km TT, reported by Graham-Paulson et al. (2018) and Fischer et al. (2015), respectively, yet with much higher end BLa (> 6 mmol L −1 ). During the TT participants maintained a constant PO, similar to Fischer et al. (2015), although ME% reduced from 19.3% in the first 4 km to 16.3% in the final 4 km. Based on this data, future constant work studies replicating handcycling TTs, should select a PO equivalent to ~ 70% PO Peak or ~ 25% greater than PO 4 (using methods similar to the current study). This high exercise intensity indicates the contribution of the anaerobic pathways to supply energy during a typical handcycling TT. Therefore, a high anaerobic capacity may play an important role in handcycling TT success, like able-bodied cycling (Støren et al. 2013), although this requires future investigation.
TT performance was significantly and highly correlated with PO Peak (r = − 0.77, P = 0.006), which supports previous research (Janssen et al. 2001;Lovell et al. 2012;de Groot et al. 2014;Fischer et al. 2015). The present study found, for the first time, that recumbent handcycling 16 km TT performance was strongly correlated with AeLT (r = − 0.68, P = 0.022) and PO 4 (r = − 0.77, P = 0.006). Contrary to the previous handcycling literature (Janssen et al. 2001;Lovell et al. 2012;de Groot et al. 2014;Fischer et al. 2015), absolute V O 2Peak (r = − 0.06) and relative V O 2Peak (r = − 0.21) were not correlated with TT performance. Therefore, these results suggest that, in highly trained handcyclists with a heterogenous injury description, PO Peak and blood lactate (AeLT and PO 4 ) are better predictors of TT performance than V O 2Peak . In able-bodied athletes, fractional utilisation of V O 2Peak , i.e. intensity at lactate threshold, has long been understood to better predict endurance performance in highly trained athletes whilst PO peak has consistently been shown to predict cycling performance, regardless of task duration or test profile (Bentley et al. 2001;Joyner and Coyle 2008). This has also been shown to be true when assessing physiological correlates to cycling performance within a paratriathlon race (Stephenson et al. 2020). As such, this finding is not unique in endurance sport but is presently confirmed in elite handcyclists.
The current study extends previous work, which was limited to handcyclists with a spinal cord injury (Lovell et al. 2012;Fischer et al. 2015;West et al. 2015;Graham-Paulson et al. 2018), by including participants with lower limb impairments and cerebral palsy, which is more representative of athletes competing at the Paralympic Games. The increased recruitable muscle mass of athletes with nonparalysed lower limbs (e.g. amputation or cerebral palsy), results in greater rates of oxygen uptake relative to athletes with complete spinal lesions (Baumgart et al. 2018). For example, in amputee athletes this is due to the capability to use lower limb musculature to brace within the handbike for greater stability and power transfer. The increased available muscle mass and the potential loss of lower limbs makes variables such as absolute V O 2Peak , relative V O 2Peak and power to weight ratio unrepresentative or misleading without scaling parameters (Goosey-Tolfrey et al. 2003). Therefore, in a population of handcyclists competing in the H3 and H4 classes [spinal lesions, lower limb amputations and cerebral palsy (UCI 2019)] variables such as PO Peak and blood lactate (AeLT and PO 4 ) are better indicators of TT performance than V O 2Peak . From this data, it may be recommended that testing for PO Peak or PO 4 is conducted with handcyclists to infer performance potential. Although assessing the latter requires specialist equipment for BLa measurement, PO 4 is also commonly used for training intensity prescription.

Limitations
An electronically braked ergometer was used to simulate the TT. The electronic braking meant that handcyclists had to pedal continually to apply power, which differs from a road TT. The corners and gradients throughout a road TT course would allow the handcyclists to recover for short periods (0-5 s). Furthermore, as the ergometer was fixed to the floor, the influence of aerodynamic drag, steering and braking were removed. Future research should aim to collect TT data on the road, in a competition if possible, and compare with laboratory-based findings. Additionally, recruiting participants from the H1, H2 and H5 classes, along with a greater number within H3 (n = 5) and H4 (n = 6) would allow a comparison between classes. Finally, the inter-individual variability of physiological parameters in the study's cohort is acknowledged. Heterogeneity is inherent in Paralympic populations due to athletes' spectra of impairments. Furthermore, variation in physiological parameters are likely largely mediated by training history and performance standard.

Conclusions
The current study revealed that the best predictors of 16 km TT performance, in a population of trained recumbent handcyclists, were PO 4 and PO Peak . It is suggested that PO 4 and PO Peak are used to infer performance level/potential rather than V O 2Peak within a population of H3 and H4 recumbent handcyclists (spinal cord lesion vs. lower limb amputation vs. cerebral palsy). A protocol equivalent to 70% PO Peak is recommended, in future studies, to replicate 16 km TT intensity.