Modelling inter-individual variability in acute and adaptive responses to interval training: insights into exercise intensity normalisation

Purpose To investigate the influence of exercise intensity normalisation on intra- and inter-individual acute and adaptive responses to an interval training programme. Methods Nineteen cyclists were split in two groups differing (only) in how exercise intensity was normalised: 80% of the maximal work rate achieved in an incremental test (%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\dot{\text{W}}$$\end{document}W˙max) vs. maximal sustainable work rate in a self-paced interval training session (%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\dot{\text{W}}$$\end{document}W˙max-SP). Testing duplicates were conducted before and after an initial control phase, during the training intervention, and at the end, enabling the estimation of inter-individual variability in adaptive responses devoid of intra-individual variability. Results Due to premature exhaustion, the median training completion rate was 88.8% for the %\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\dot{\text{W}}$$\end{document}W˙max group, but 100% for the %\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\dot{\text{W}}$$\end{document}W˙max-SP the group. Ratings of perceived exertion and heart rates were not sensitive to how intensity was normalised, manifesting similar inter-individual variability, although intra-individual variability was minimised for the %\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\dot{\text{W}}$$\end{document}W˙max-SP group. Amongst six adaptive response variables, there was evidence of individual response for only maximal oxygen uptake (standard deviation: 0.027 L·min−1·week−1) and self-paced interval training performance (standard deviation: 1.451 W·week−1). However, inter-individual variability magnitudes were similar between groups. Average adaptive responses were also similar between groups across all variables. Conclusions To normalise completion rates of interval training, %\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\dot{\text{W}}$$\end{document}W˙max-SP should be used to prescribe relative intensity. However, the variability in adaptive responses to training may not reflect how exercise intensity is normalised, underlining the complexity of the exercise dose–adaptation relationship. True inter-individual variability in adaptive responses cannot always be identified when intra-individual variability is accounted for. Supplementary Information The online version contains supplementary material available at 10.1007/s00421-023-05340-y.


Introduction
It is commonly believed that the magnitude of physiological adaptations varies between individuals undertaking the same exercise training programme (Bouchard et al. 1999;Vollaard et al. 2009;Coakley and Passfield 2018;Montero and Lundby 2017;Astorino et al. 2018;Bonafiglia et al. 2019;McLellan and Skinner 1981;Preobrazenski et al. 2019;Weatherwax et al. 2019;Del Giudice et al. 2020;Hecksteden et al. 2018b).A major factor behind this phenomenon is suggested to be genetics (Mann et al. 2014;Meyler et al. 2021), estimated to account for approximately 50% of the changes in maximal oxygen uptake ( VO 2max ) (Bouchard et al. 1999(Bouchard et al. , 2011)).However, it has been proposed that methods of exercise intensity normalisation in experimental studies do not provide comparable metabolic stress across participants (Mann et al. 2013;Iannetta et al. 2020;Jamnick et al. 2020;Vollaard et al. 2009;Meyler et al. 2021Meyler et al. , 2023)), contributing to variability in the extent to which training adaptations occur (Mann et al. 2014;Meyler et al. 2021Meyler et al. , 2023)).For example, when exercise intensity is normalised as a percentage of maximal heart rate, a method with known limitations (Katch et al. 1978;Mann et al. 2013;Iannetta et al. 2020;Jamnick et al. 2020), VO 2max changes following identical training interventions separated by a washout are only moderately correlated (r = 0.31) (Del Giudice et al. 2020).Refining the scientific basis of exercise intensity prescription is, therefore, crucial to understand adaptive response heterogeneity.
The optimal method for exercise intensity normalisation may vary depending on target population, intensity domain of training (i.e.moderate, heavy, very heavy, or severe; see Rossiter (2011) for review), and exercise pattern (i.e.continuous or intermittent) (Mann et al. 2013;Jamnick et al. 2020;Meyler et al. 2023).Moreover, conflicting evidence exists as to whether certain intensity prescription methods could minimise adaptive variability (Weatherwax et al. 2019;McLellan and Skinner 1981;Preobrazenski et al. 2019).Typically, researchers compare two groups undertaking the same training programme but using different normalisation methods to set individual work rates (Weatherwax et al. 2019;McLellan and Skinner 1981;Preobrazenski et al. 2019).In this respect, Weatherwax et al. (2019) reported a reduced interindividual variability in VO 2max adaptive responses when exercise intensity domains were individually accounted for, compared with when they were not.McLellan and Skinner (1981), however, reported no differences.Preobrazenski et al. (2019) showed no differences in the magnitude of inter-individual variability of several adaptive responses, including VO 2max , when prescriptions based on the maximal work rate from an incremental test (% Ẇmax ) and the talk test were compared (see Reed and Pipe (2014) for talk test details).Nevertheless, they also found that the mean blood lactate concentration ([La − ]) of the first training session was positively associated with VO 2max changes within the % Ẇmax group (Preobrazenski et al. 2019), providing some evidence that individuals experiencing greater metabolic stress may also express larger adaptive response (and vice versa) (Mann et al. 2014(Mann et al. , 2013)).These inconsistent findings underscore the need for further research to ascertain the extent to which exercise intensity normalisation affects adaptive response variability, and in which contexts.
Interestingly, Montero and Lundby (2017) have demonstrated that a maximised training dose is essential for enhancing the maximal work rate achieved in an incremental test ( Ẇmax ) and VO 2max across all individuals within a study cohort.This suggests that McLellan and Skinner (1981) and Preobrazenski et al. (2019) may not have provided their participants with sufficient training stress, making it difficult to untangle potential between-group differences in adaptive variability (Joyner and Lundby 2018).The only study investigating exercise intensity normalisation that assessed adaptive response heterogeneity of a more intense, interval training intervention, compared results with the available literature as opposed to a comparative group, hampering interpretation of their findings (Astorino et al. 2018).Hence, filling this gap is important to elucidate this issue.
For intensive training, maximal self-paced intervals have been employed as a method of exercise intensity normalisation, both in cross-sectional (Brosnan et al. 2000;Villerius et al. 2008;Nicolò et al. 2014) and longitudinal interventions (Seiler and Sylta 2017;Seiler et al. 2013;Rønnestad et al. 2020).Frequently described as "how elite athletes train" (Brosnan et al. 2000;Villerius et al. 2008;Rønnestad et al. 2020;Seiler et al. 2013;Seiler and Sylta 2017;Nicolò et al. 2014), this approach is based on the premise that there exists an individualised maximal sustainable work rate for a given interval training format.Provided that work intervals are performed within the very heavy-intensity domain, this concept has precedents in the hyperbolic relationship between work rate and time to exhaustion (Ferguson et al. 2013;Jones and Vanhatalo 2017;Meyler et al. 2023).Surprisingly, there have been no attempts to assess the effectiveness of this method of intensity normalisation in comparison with other approaches such as % Ẇmax .
In the present study, we investigated inter-individual variability in acute and chronic (i.e.adaptive) responses to a training programme in two groups of cyclists.It was hypothesised that the group in which training intensity was prescribed relative to the maximal sustainable work rate in a self-paced interval training session (% Ẇmax-SP ) would exhibit less inter-individual variability in acute exercise responses, leading to less variability in adaptive responses, compared with the group in which training intensity was prescribed as % Ẇmax .We also hypothesised that the % Ẇmax group would demonstrate a higher proportion of unfinished training sessions, due to miscalculated work rate targets leading to premature exhaustion, potentially compromising the group's average adaptive responses.

Ethics approval
The research protocols were submitted to and approved by the Research Ethics Committee at the University of Kent (Prop 18_2018_19), in compliance with the Declaration of Helsinki, except for registration in a database.All participants provided written informed consent prior to participating in this study.

Study design
Participants were involved for 16 weeks (see Table 1), with weeks designated for testing (4 weeks), control (6 weeks), and training intervention (6 weeks).While distinct methods of exercise intensity prescription were used for each group during the training intervention, testing and control phases consisted of identical requirements for all participants.Testing before and after the control phase served as a control against which to gauge the effects of the training interventions (Voisin et al. 2019).Moreover, the testing phase at week 11 enabled the estimation of inter-individual variability in adaptive responses without the need for repeating the training intervention (Hecksteden et al. 2018b).This experimental design makes it possible to estimate inter-individual variability in adaptive responses devoid of intra-individual variability (Voisin et al. 2019;Hecksteden et al. 2018b).

Testing phase
At consistent times of the day, participants visited the laboratory thrice, at least 48 h apart.In both the first and second visits, participants completed a lactate accumulation test and an incremental test to exhaustion (i.e.duplicate measures were averaged; see page 1 of supplementary material for reliability estimates).In the third visit, participants performed a self-paced interval training session.They were instructed to refrain from intense exercise before testing and to prepare as for competition.Participants were also requested to standardise meals 24 h prior and to consume their last large meal at least 2 h before arrival.The consumption of caffeine was not allowed in ] was ≥ 2.5 mmol•L −1 ), and terminating when [La − ] reached ≥ 4 mmol•L −1 .Blood samples taken from a fingertip (at the last 30 s of each 5-min bout) were immediately analysed for [La − ].Power output associated with 4 mmol•L −1 [La − ], sometimes referred to as the onset of blood lactate accumulation (Sjödin and Jacobs 1981), was calculated for each cyclist from the relationship between [La − ] and power output in the last two stages.Before the start, participants chose their preferred cadence for the entire test (91 ± 4 rev•min −1 ).Both the work rates and cadence of the first lactate accumulation test were held constant throughout the study.Breath-by-breath gas exchanges were monitored throughout the test and subsequently smoothed to 30-s averages.Gross efficiency, measured as the ratio between power output and energy expenditure (Hopker et al. 2009), was calculated at 150 W from the mean gas exchanges in the last 2.5 min of the stage.Energy expenditure was estimated assuming negligible protein oxidation according with the equations of Péronnet and Massicotte (1991).All participants met the criterion of a respiratory exchange ratio ≤ 1.0 in all tests.
After the lactate accumulation test, participants cycled for 10 min at a power output between 50 and 100 W. Subsequently, participants completed an incremental test in which work rate increased continuously at 25 W•min −1 until voluntary exhaustion, or participants' inability to maintain cadence above 70 rev•min −1 .Breath-by-breath gas exchanges were monitored throughout the test and subsequently smoothed to 15-s averages.VO 2max was identified as the highest 60-s mean oxygen uptake, and Ẇmax as the mean power output of the last 60 s.Immediately after the incremental test, a blood sample was taken from a fingertip to establish [La − ], and peak rating of perceived exertion (RPE) was noted.
The self-paced interval training session consisted of six 4-min work intervals interspersed with 2-min active recovery.Participants started immediately after a 10-min warmup at power outputs between 100 and 150 W. They were required to produce the highest possible amount of work to establish Ẇmax-SP (i.e.highest possible mean power output across all six work intervals) and received instructions to pace themselves by keeping power reasonably constant between and within work intervals.Recovery intervals had to be performed at power outputs ≤ 70 W. Heart rate was measured as the last-minute average of each work interval.RPE was noted immediately after each work interval.

Control phase
During this phase, participants did not attend the laboratory.However, they were required to keep their weekly training duration similar to the last two weeks before joining the study, and to avoid structured interval training.

Training intervention phase
Due to the relatively small number of participants recruited for this study, the first participant was truly randomised, with subsequent participants allocated to one of the two training interventions to keep groups closely matched with regard to dependent variables; i.e. minimisation approach (Hecksteden et al. 2018a).Participants were blinded to their group assignment and unaware of the methods of intensity normalisation used.Both groups attended the laboratory twice per week, at least 72 h apart, to perform interval training sessions consisting of 4-min work intervals interspersed with 2-min active recovery, at predefined work rates.Six training sessions were performed from weeks 8 to 10, and another six from weeks 12 to 14.While in one training intervention (% Ẇmax ), the work intervals were performed at 80% Ẇmax meas- ured on the first incremental test (i.e.visit one of testing; see Table 1); in the other (% Ẇmax-SP ), the work intervals were performed at 100% Ẇmax-SP .Recovery intervals were performed at 20% of the work rate prescribed for the work intervals, irrespective of group allocation; i.e. 0.2•(mean [80% Ẇmax, 100% Ẇmax-SP] ).Participants of both groups were prescribed six work intervals in each training session, except for weeks 8 and 12, in which five work intervals were prescribed to boost their confidence that sessions could be completed.Despite strong encouragement, voluntary exhaustion or inability to maintain cadence above 70 rev•min −1 were utilised as criteria to establish individual completion rates in the event of premature termination.Cadence was recorded as the average of each work interval (or the average of completed duration in case of exhaustion), and heart rate as the last-minute average of each work interval (or the average of completed duration if shorter than one minute).RPE was noted immediately after each work interval or at exhaustion.All interval training sessions commenced with a 15-min warm-up and finish with a 3-min cool-down, at, respectively, 60% and 40% of the work rate prescribed for the work intervals, irrespective of group allocation; i.e. 0.6•(mean [80% Ẇmax, 100% Ẇmax−SP] ), and 0.4•(mean [80% Ẇmax, 100% Ẇmax−SP] ).As VO 2max gains have been shown to plateau after 3 weeks of high-intensity training at the same work rates (Hickson et al. 1981), training targets were re-adjusted following the testing phase of week 11, no matter if participants exhibited an increase or a decrease in performance.Participants were instructed to perform their remaining training sessions (i.e.outside the laboratory) at work rates below the power output associated with 4 mmol•L −1 [La − ], and to keep weekly training duration similar to the control phase.

Equipment
Cyclists used their own bikes mounted on a cycle ergometer (Cyclus 2, RBM Elektronik-Automation, Leipzig, Germany).For the lactate accumulation tests, incremental tests, and predefined interval training sessions, the ergometer was set at power mode (i.e.cadence independent).For the selfpaced interval training sessions, the ergometer was set at inclination mode (i.e.0% gradient; cadence dependent), and participants were required to change gears, as if they were riding outdoors.Heart rate was continuously monitored during all sessions through an ANT + belt transmitter (Cyclus 2, RBM Elektronik-Automation, Leipzig, Germany).Elapsed time, power output, heart rate, and cadence were not concealed from participants.
Breath-by-breath gas exchanges were monitored through a metabolic cart (MetaLyzer 3B, Cortex Biophysik, Leipzig, Germany).Prior to every test, calibration was performed according to the manufacturer's instructions.[La − ] was assessed using an automatic analyser (Biosen C-Line, EKF Diagnostics, Penarth, UK).RPE was assessed based on the 6-20 Borg's scale (Borg 1982).The same trained experimenter conducted all testing and training sessions to minimise procedural variability.

Data analysis
Data were assessed for normality using Shapiro-Wilk's test and normal quantile plots.To investigate between-group differences in target work rates for the training sessions, independent samples t tests were used.Training intervention completion rates were assessed for between-group difference using a Mann-Whitney test.Training RPE, heart rate, and cadence were investigated via linear mixed models with participant as a random effect, and group, training session, and work interval as fixed effects.To identify evidence of between-group differences in the magnitude of inter-and intra-individual variability, models were fitted with homogeneous and heterogeneous inter-and intra-individual variance structures for group.
To investigate between-group differences in the adaptive response variables (i.e.VO 2max , Ẇmax , power output associ- ated with 4 mmol•L −1 [La − ], gross efficiency, Ẇmax-SP , and body mass) prior to the intervention, independent samples t tests were used.Differences in adaptive response variables between testing weeks 0, 7, 11, and 15 were assessed using repeated-measures analysis of variance, with Bonferroni pairwise comparisons used to identify where significant differences existed within the data.Linear mixed models, with participant as a random effect, and group and testing occasion as fixed effects, were used to test for a group effect on adaptive response variables' change from week 0 while controlling for their absolute baseline scores.To investigate inter-individual variability in adaptive responses to training, piecewise linear mixed models were used with participant and participant-by-intervention week interaction as random effects, and control week and intervention week as fixed effects.The standard error of the participant-by-intervention week interaction was used to calculate confidence intervals associated with individual adaptive responses.Individuals whose confidence intervals overlapped '0' were considered non-responders, whereas those whose confidence intervals did not overlap '0' were considered responders or adverse responders based on a positive or negative response, respectively.Optimal models were selected using likelihood ratio tests.Pearson's correlation was employed to examine the relationship between modelled adaptive responses.
Data were analysed using Prism 8 (GraphPad, San Diego, USA), with model fitting performed in R 4.0.4(R Foundation for Statistical Computing, Vienna, Austria).Significance level was set at P ≤ 0.05, and confidence level was set at 95%. Results are presented as mean ± SD unless otherwise stated.The reader unfamiliarised with linear mixed models is referred to Brown (2021), Naumova et al. (2001), and Pinheiro and Bates (2020).

Training intervention
All participants attended all sessions of the training intervention.Target work rates for training are presented in Table 2.No between-group differences were detected for any of the target work rates, at either the first or second half of the intervention (all P ≥ 0.220).
Rates of training completion are presented in Table 3. Due to premature exhaustion, participants of the % Ẇmax group did not complete the entire sessions as often as participants of the % Ẇmax-SP group, resulting in lower overall completion rates.
There were no between-group differences for RPE, heart rate, or cadence (see % Ẇmax-SP estimates on page 2 of sup- plementary material), and there was no evidence of betweengroup differences in the magnitude of inter-individual variability for these variables (Table 4).In contrast, there was evidence of lower intra-individual variability in acute training responses for the % Ẇmax-SP group (Table 4).

Training outcomes
Participants attended all testing sessions, except for one participant of the % Ẇmax-SP group that did not attend the second visit of week 0. No between-group differences were detected for any of the adaptive response variables, at either week 0 or week 7 (all P ≥ 0.233), suggesting the participant allocation into groups was successful.During the self-paced session, RPEs associated with each work interval were not different between weeks (all P ≥ 0.102).This is despite an increased heart rate and power output after the start of the training intervention (see Fig. 1 for details), suggesting that participants consistently adhered to instructions.
When considering all participants together, changes over the 16 weeks of the study were evident for all adaptive response variables, except gross efficiency (Fig. 2).During the control phase (from week 0 to 7), Ẇmax and power output associated with 4 mmol•L −1 [La − ] increased by 11 W (P < 0.001) and 8 W (P = 0.027), respectively, but there was no change for any other adaptive response (all P ≥ 0.414).During the training intervention (from week 7 to 15), VO 2max increased by 0.215 L•min −1 (P = 0.038), Ẇmax increased by 14 W (P < 0.001), and body mass increased by 1.1 kg (P = 0.009).While there was also an increasing trend for Ẇmax-SP from week 7 to 15 (8 W, P = 0.085), it reached statistical significance only compared with week 0 (12 W, P = 0.014).Power output associated with 4 mmol•L −1 [La − ] did not increase further from week 7 (P = 0.636).When changes in adaptive response variables were modelled, a group difference was evident only for gross efficiency (% Ẇmax-SP group: − 0.8%, P = 0.044).However, adding a testing occasion-group interaction did not further improve the gross efficiency model (P = 0.119).
After accounting for intra-individual variability associated with control and intervention phases (see Table 5 for fixed effects), there was evidence of inter-individual variability in adaptive responses for VO 2max (P = 0.003 -see page 3 of supplementary material) and Ẇmax-SP (P = 0.001 -see page 4 of supplementary material).However, adding an intervention week-group interaction as a fixed or random effect did not improve the models (P ≥ 0.197 for all model comparisons), indicating that there was no evidence of between-group differences in the magnitude of inter-individual variability for either variable.Accordingly, the confidence intervals for the SD of the individual intervention-week coefficients overlapped substantially ( VO 2max : 0.017-0.045L•min −1 •week −1 for the % Ẇmax group, and 0.014-0.040L•min −1 •week −1 for the % Ẇmax-SP group; Ẇmax-SP : 0.948-2.517W•week −1 for the % Ẇmax group, and 0.762-2.161W•week −1 for the % Ẇmax-SP group).Unlike VO 2max and Ẇmax-SP , there was no evidence of inter-individual variability in adaptive responses for Ẇmax (P = 0.207), power output associated with 4 mmol•L −1 [La − ] (P = 0.466), gross efficiency (P = 0.348), or body mass (P = 0.173).
Both VO 2max and Ẇmax-SP models yielded large resid- ual errors relative to the variability in intervention-week slopes (see pages 3 and 4 of supplementary material), resulting in wide confidence intervals for individual responses (Fig. 3-panels A and B), and making it difficult to categorise most participants.There were three and two responders for VO 2max in the % Ẇmax and % Ẇmax-SP groups, respectively, with the remaining participants being categorised as non-responders.There were three responders, six non-responders, and one adverse responder for Ẇ max-SP in the % Ẇmax group; and one responder, seven non- responders, and one adverse responder in the % Ẇmax-SP group.However, modelled VO 2max and Ẇmax-SP responses were not correlated (Fig. 3-panel C).

Discussion
We investigated acute and chronic (i.e.adaptive) responses to a training programme in which recreationally trained cyclists were split into two groups differing in exercise intensity normalisation methods, but with identical prescriptions otherwise.The main findings are: a) performance in a maximal self-paced interval training session (i.e. % Ẇmax-SP ) may be used to normalise the exercise intensity of interval training performed at predefined work rates, particularly if premature exhaustion is to be avoided; b) after accounting for sources of intra-individual variability, there was evidence of adaptive response heterogeneity for VO 2max and Ẇmax-SP only, but no between-group dif- ferences in magnitude; and c) average adaptive responses were not different between groups, meaning that the higher prevalence of incomplete training sessions in the % Ẇmax group, due to premature exhaustion, did not compromise participants' training effect.

Self-paced performance to normalise interval training intensity
Research on exercise intensity normalisation has been ongoing since the late 1970's, with Katch et al. (1978) and McLellan and Skinner (1985) amongst the first to challenge the efficacy of percentages of maximal heart rate and VO 2max .Subsequent studies have also identified shortcomings in other traditional methods, namely % Ẇmax and VO 2max /heart rate reserve (Iannetta et al. 2020;Jamnick et al. 2020;Mann et al. 2013;Marini et al. 2021).Yet, all these methods continue to be used (Vollaard et al. 2009;Bouchard et al. 1999Bouchard et al. , 2011;;Hecksteden et al. 2018b;Coakley and Passfield 2018;Bonafiglia et al. 2019;Montero and Lundby 2017;Del Giudice et al. 2020), most likely due to the limited empirical support for alternative approaches, such as the delta concept, which considers different physiological anchors (Lansley et al. 2011;McLellan and Skinner 1985;Meyler et al. 2023), critical power modelling (Ferguson et al. 2013;Jones and Vanhatalo 2017;Meyler et al. 2023), and maximal self-paced intervals (Villerius et al. 2008;Seiler and Sylta 2017;Nicolò et al. 2014;Brosnan et al. 2000).Accordingly, the present investigation reveals that prescribing exercise intensity of interval training as 100% Ẇmax-SP minimises performance variability between individuals compared with 80% Ẇmax .Only occasionally (8.3% of the sessions), did participants of the % Ẇmax-SP group experience premature exhaustion throughout the training intervention, with a median completion rate of 100%.In contrast, premature exhaustion was very common (44.2% of the sessions) amongst participants of the % Ẇmax group, with a median completion rate of 88.8% (see Table 3 for individual data).This was despite a similar exercise intensity between groups on average (see Table 2).These data thus reinforce previous critiques of % Ẇmax (Iannetta et al. 2020; Jamnick et al.  2020) and substantiate the use of 100% Ẇmax-SP for interval training intensity normalisation.
The relationship between work rate and sustainable duration is largely individual, particularly for intermittent exercise (Ferguson et al. 2013;Jones and Vanhatalo 2017;Meyler et al. 2023).Accordingly, it makes sense to establish a common duration, and allow individuals to select the maximal sustainable work rate, instead of presuming that a single variable (e.g. % Ẇmax ) is able to predict their exercise capacity.Prescribing exercise at 100% Ẇmax-SP neverthe- less assumes that a) individuals can pace maximal efforts to deliver performances consistent with their capacity; and b) self-and ergometer-paced performances are equivalent when the mean work rate is the same, which may not be universally true (Black et al. 2015;Thomas et al. 2013).Crucially, our data ease concerns about both assumptions.In line with other studies (Villerius et al. 2008;Seiler and Sylta 2017;Nicolò et al. 2014;Brosnan et al. 2000), RPE increased quasi-linearly during the self-paced interval training sessions, approaching 20 in the last work interval (see Fig. 1).Moreover, no differences between testing occasions were detected for the RPEs associated with each work interval, despite an increased heart rate and performance after the start of the training intervention.These observations suggest that performance gains most likely reflected an improved exercise capacity rather than a different approach to the task.As for the self-vs.ergometer-paced performances, While it is tempting to conclude that maximal self-paced intervals should be incorporated into training, replacing intervals at predefined work rates (Rønnestad et al. 2020;Seiler and Sylta 2017;Seiler et al. 2013), some physiological responses such as oxygen uptake are sensitive to large power output variations within (Bossi et al. 2020) and between work intervals (Ferguson et al. 2013).Whether variability in power output distribution would contribute to increased inter-individual variability in acute and adaptive responses to training is unclear.We, therefore, preliminarily recommend that maximal self-paced intervals are used only to determine % Ẇmax-SP .

RPE and heart rate as indicators of exercise response variability
Even though a lower performance variability between individuals was detected for the % Ẇmax-SP compared with the % Ẇmax group, RPE and heart rate data only partially corroborate this finding.Within the % Ẇmax-SP group, the magnitude of intra-, but not inter-individual variability, was lower for both RPE and heart rate (see Table 4).From an intra-individual perspective, this outcome likely stems from the fact that participants of the % Ẇmax-SP group con- sistently completed their training sessions, stopping at the same timepoint, whereas premature exhaustion occurred at different timepoints when participants of the % Ẇmax group struggled.This is expected, due to normal day-today performance variability (Midgley et al. 2007) plus the combined effects of gradual training adaptation and work rate adjustment at week 11 (see RPE on page 2 of supplementary material for evidence of the latter effects).From an inter-individual perspective, our findings align with those of Meyler et al. (2023), who showed that heart rate, oxygen uptake, and [La − ] do not always reflect between-group differences in inter-individual variability in performance.Therefore, physiological and perceptual responses to highintensity training may not be as sensitive as performance to quantify variability and inform the normalisation of exercise intensity.Alternatively, an effective normalisation of exercise intensity based on performance may not ensure uniform physiological and perceptual responses across individuals.While more studies are required to elucidate these hypotheses, it is important to underscore that RPE and heart rate data were modelled to factor in the fixed effects of group, training session, and work interval, meaning that our estimates are conservative compared with other studies (Lansley  . 2011;Scharhag-Rosenberger et al. 2010;Vollaard et al. 2009;Katch et al. 1978;Meyler et al. 2023), and certainly closer to the true inter-individual variability (Voisin et al. 2019;Hecksteden et al. 2018b).Unless the pitfalls inherent to the analysis of raw variability are avoided (Voisin et al. 2019;Hecksteden et al. 2018b;Williamson et al. 2017;Atkinson et al. 2019), future investigations are unlikely to clarify our findings.

Adaptive response heterogeneity and the impact of intra-individual variability
Since the influential work of Bouchard et al. (1999) investigating the heritability of VO 2max responses to training, several authors have claimed that the extent to which each individual adapts to a standardised programme is fairly unique (Vollaard et al. 2009;Coakley and Passfield 2018;Bonafiglia et al. 2019;Preobrazenski et al. 2019;Astorino et al. 2018;Weatherwax et al. 2019;Montero and Lundby 2017;Del Giudice et al. 2020;Hecksteden et al. 2018b).However, apart from Hecksteden et al. (2018b), they did not account for all sources of variability affecting the observed inter-individual variability (Voisin et al. 2019;Hecksteden et al. 2018b;Williamson et al. 2017;Atkinson et al. 2019), prompting questions as to the existence of true adaptive response heterogeneity (Williamson et al. 2017).By following the best design and analytical practices (Voisin et al. 2019;Hecksteden et al. 2018b), we demonstrate that variability between individuals in Ẇmax , power output associ- ated with 4 mmol•L −1 [La − ], gross efficiency, and body mass responses to an interval training programme is likely a manifestation of intra-individual variability associated with the control phase and/or the intervention phase itself.This interpretation is strengthened by the use of averaged duplicate measures to minimise day-to-day biological and technical fluctuations, facilitating the identification of a true inter-individual variability (if present) (Voisin et al. 2019).Adaptive response heterogeneity was nevertheless detected for VO 2max and Ẇmax-SP , even though the latter variable was not analysed in duplicates.Together, these distinct outcomes indicate that inter-individual variability in training adaptations can occur, although it may be difficult to demonstrate statistically when all confounding sources of variability are accounted for and/ or the magnitude of changes associated with an intervention is relatively small.
Upon re-analysis of the HERITAGE Family Study data (Bouchard et al. 1999), Shephard et al. (2004) have demonstrated that the true inter-individual variability in VO 2max adaptive responses was much smaller than originally estimated.The raw SD of 0.010 L•min −1 •week −1 represented in reality 0.007 or 0.006 L•min −1 •week −1 , whether a 2-day or a 2-week test-retest coefficient of variation for VO 2max was considered, respectively, to factor in the intra-individual variability expected for assessments conducted 20 weeks apart.Given that Shephard et al. (2004) were not able to account for intra-individual variability associated with identical training programmes, either through repeated testing or repeated interventions (Voisin et al. 2019;Hecksteden et al. 2018b), the 0.007-0.006L•min −1 •week −1 figure likely still overestimates the true inter-individual variability.Accordingly, the question that arises is whether recreationally trained cyclists, as employed herein, are more susceptible to adaptive response heterogeneity than sedentary individuals, as employed in the HERITAGE Family Study (Bouchard et al. 1999) and elsewhere (Hecksteden et al. 2018b).The SD for VO 2max responses reached 0.027 L•min −1 •week −1 in the current study, notably higher than those estimated by Shephard et al. (2004) (see above), and Hecksteden et al. (2018b) as 0.042 ml•kg −1 •min −1 •week −1 .In theory, sedentary individuals have untapped genetic potential for VO 2max improvements, unlike recreationally trained cyclists, resulting in inconsistent adaptive gains within the latter cohort.Further studies are necessary to test this hypothesis.
Interestingly, we also found evidence of adaptive response heterogeneity for Ẇmax-SP , implying that the extent to which participants improved intermittent self-paced performance varied (with a SD of 1.451 W•week −1 ).While no comparable studies exist, this may suggest that Ẇmax-SP is char- acterised by a high signal-to-noise ratio, being sensitive to small changes in exercise capacity, and thus suitable as an intensity prescription benchmark.Conversely, adaptive response heterogeneity was not detected for Ẇmax .These findings are consistent with the fact that 100% Ẇmax-SP suc- cessfully normalised the completion rates of interval training, while 80% Ẇmax led to premature exhaustion in 44% of the sessions.

Similar group-level adaptive responses despite contrasting rates of training completion
Despite the evidence in favour of Ẇmax-SP as a bench- mark for interval training prescription, there were no between-group differences in the magnitude of adaptive variability, either for VO 2max or Ẇmax-SP .Contrary to our hypothesis, adaptive responses may be too complex to reflect the manipulation of a single element of training prescription (i.e.how exercise intensity is normalised).(Bouchard et al. 1999(Bouchard et al. , 2011)), the isolated or combined impact of training prescription, nutrition, and recovery remains unclear.Given that the average changes in VO 2max , % Ẇmax , power output associated with 4 mmol•L −1 [La − ], gross efficiency, Ẇmax-SP , and body mass were also not different between groups, it may be speculated that fine-tuning exercise intensity is irrelevant from an adaptive point of view.To shed light on this possibility, a literature overview is instructive.McLellan and Skinner (1981) compared the inter-individual variability in VO 2max responses between groups; one in which exercise intensity was normalised as % VO 2max , and another in which exercise intensity was normalised relative to the first ventilatory threshold (%VT 1 ).No between-group differences were detected for the magnitude of inter-individual variability.Likewise, the dataset of Preobrazenski et al. (2019), which included VO 2max , Ẇmax , and power output associated with 4 mmol•L −1 [La − ], displays a similar magnitude of inter-individual variability between groups (i.e.65% Ẇmax vs. first negative stage of the talk test).Current results, therefore, corroborate these previous findings.In contrast, by comparing the heart rate reserve method with an individualised approach using the heart rate associated with each ventilatory threshold, Weatherwax et al. (2019) concluded that how exercise intensity is normalised affects the inter-individual variability in VO 2max responses.However, a detailed inspection reveals that the individualised approach group trained at a higher intensity on average and made a larger VO 2max gain.As Weatherwax et al. (2019) relied on the responder counting approach, which has been shown to reflect the magnitude of mean differences rather than inter-individual differences (Atkinson et al. 2019), their inference could be questioned.
Even though the evidence mostly indicates that adaptive response heterogeneity is not directly influenced by how exercise intensity is normalised, we cannot discard a small contribution.For example, Preobrazenski et al. (2019)  revealed that the mean [La − ] of the first training session in a series was positively associated with VO 2max changes within their 65% Ẇmax group.Similarly, Gaskill et al. (2001)  showed that training intensity in the HERITAGE Family Study, originally normalised as % VO 2max (Bouchard et al. 1999), accounted for 26% of the gains in the oxygen uptake associated with VT 1 when expressed as %VT 1 .In other words, the higher the intensity relative to VT 1 , the greater the VT 1 gain (Gaskill et al. 2001).These two studies provide some evidence that the metabolic stress experienced by each individual is associated with their adaptive response (Mann et al. 2014(Mann et al. , 2013)).Thus, sample sizes of less than twenty participants per group, as employed herein and elsewhere (McLellan and Skinner 1981;Preobrazenski et al. 2019;Weatherwax et al. 2019), may not be enough to investigate adaptive response heterogeneity from an intensity normalisation perspective.This possibility requires careful consideration by those designing future studies.

Limited utility of the responder counting approach
The responder counting approach has been frequently adopted to investigate inter-individual variability in adaptive responses to a training programme (Bouchard et al. 1999;Vollaard et al. 2009;Coakley and Passfield 2018;Montero and Lundby 2017;Astorino et al. 2018;Bonafiglia et al. 2019;Weatherwax et al. 2019;Del Giudice et al. 2020;Hecksteden et al. 2018b).While definitions for responders, non-responders, and adverse responders vary between studies, constituting a problem in itself (see Hecksteden et al. (2018b) and Voisin et al. (2019) for overview), it has been argued that this approach is flawed for two main reasons: a) observed responses may simply reflect intra-individual variability of different sorts (Voisin et al. 2019;Hecksteden et al. 2018b;Williamson et al. 2017;Atkinson et al. 2019); and b) the number of responders, non-responders, and adverse responders of a sample is expected to conform with a normal distribution, reflecting deviations of the mean, rather than the true magnitude of adaptive response heterogeneity (Atkinson et al. 2019).In light of these criticisms, we used the standard error of the participant-by-intervention week interaction to calculate confidence intervals associated with individual adaptive responses for VO 2max and Ẇmax-SP .As results demonstrate, most participants were clas- sified as non-responders due to the uncertainty with which individual responses are estimated.For both VO 2max and Ẇmax-SP models, there was a large residual error compared with the inter-individual variability in intervention-week slopes, suggesting a great level of intra-individual variability.These findings, therefore, corroborate previous demonstrations that the responder counting approach may be untenable (Atkinson et al. 2019;Hecksteden et al. 2018b).

Investigating inter-individual variability in performance as a meaningful target
Interestingly, modelled gains in VO 2max and Ẇmax-SP resulting from the training intervention did not correlate, despite VO 2max being generally considered the main endurance performance determinant (Joyner and Coyle 2008).Björklund et al. (2007) also found no correlation between VO 2max and time to exhaustion during an interval training session, whilst Daniels et al. (1978) and Vollaard et al. (2009) reported no association between changes in VO 2max and time-trial performance following a training intervention.Together, these findings suggest that, within athletic populations, the scientific interest for inter-individual variability in adaptive responses should perhaps shift from V O 2max to performance.

Methodological considerations for future studies
Research on exercise intensity normalisation predominantly falls into three categories: a) those that demonstrate the variability in work rate targets based on percentages of a maximal benchmark (e.g.70% VO 2max , 60% Ẇmax ) in relation to the intensity domains of exercise (Katch et al. 1978;Iannetta et al. 2020); b) those that assess inter-individual variability in performance, physiological, and/or perceptual responses to acute bouts of exercise (McLellan and Skinner 1985;Lansley et al. 2011;Scharhag-Rosenberger et al. 2010;Vollaard et al. 2009;Meyler et al. 2023); and c) those that compare inter-individual variability in training adaptations between groups differing in how exercise intensity is normalised (McLellan and Skinner 1981;Preobrazenski et al. 2019;Weatherwax et al. 2019).The rationale linking these rather distinct experimental designs is that a large inter-individual variability in acute exercise responses is likely to manifest as a large variability in adaptive responses to a training programme (Mann et al. 2014;Meyler et al. 2021).This theory draws upon molecular biology evidence that chronic adaptations to training originate from the cumulative effects of transient homeostatic perturbations associated with each exercise session (Perry et al. 2010;Egan and Zierath 2013).The findings of the present study and others (McLellan and Skinner 1981;Preobrazenski et al. 2019) suggest that future investigations should look at acute exercise responses to different methods of exercise intensity normalisation rather than adaptive responses, thus avoiding waste of time and resources.

Limitations
One relevant characteristic of our study is that the total training load was not only dependent on laboratory-based training sessions, unlike comparable investigations (Weatherwax et al. 2019;McLellan and Skinner 1981;Preobrazenski et al. 2019).While the experimental design (with a control phase prior to the intervention) and analytical procedures permitted that this source of variability was accounted for, we cannot discard the possibility that participants changed their habitual training routine after the start of the intervention, affecting the magnitude of inter-individual variability in adaptive responses.Despite our efforts to persuade participants of the importance of training load documentation, very few did so with sufficient detail to provide insights in this respect.

Conclusions
In summary, this study suggests that Ẇmax-SP may be used to normalise the intensity of interval training performed at predefined work rates.This approach prevents premature exhaustion, although without necessarily minimising interindividual variability in RPE and heart rate.The inter-individual variability in adaptive responses to training, while only detected for VO 2max and Ẇmax-SP amongst six variables, was similar in magnitude between groups differing (only) in how exercise intensity was normalised (i.e.100% Ẇmax-SP vs. 80% Ẇmax ).Furthermore, no between-group differences in the magnitude of average responses to training were observed across all variables.These results underline the complexity of the relationship between acute training dose and chronic adaptations.From a methodological point of view, true inter-individual variability in adaptive responses cannot always be identified when intra-individual variability is accounted for.

Fig. 3
Fig. 3 Individual estimates with confidence intervals for weekly changes in maximal oxygen uptake ( Δ .VO 2max -panel A) and maximal sustainable work rate in a self-paced interval training session (Δ Ẇmax-SP -panel B) beyond the increase associated with the control phase, and related scatterplot (panel C)

Table 1
Timeline of the study % Ẇmax , group in which training intensity was prescribed relative to the maximal work rate achieved in an incremental test; % Ẇmax-SP , group in which training intensity was prescribed relative to the maximal sustainable work rate in a self-paced interval training session; N/A, not applicable

Table 2
Target work rates for training (W•kg −1 ) % Ẇmax , group in which training intensity was prescribed relative to the maximal work rate achieved in an incremental test; % Ẇmax-SP , group in which training intensity was prescribed relative to the maximal sustainable work rate in a self-paced interval training session

Table 4
Variability in acute exercise responses (SD) % Ẇmax , group in which training intensity was prescribed relative to the maximal work rate achieved in an incremental test; % Ẇmax-SP , group in which training intensity was prescribed relative to the maximal sustainable work rate in a self-paced interval training session; RPE, ratings of perceived exertion

Table 5
Fixed effects upon adaptive responses to training Formula: dependent variable = intercept + control week coefficient • x + intervention week coefficient • x.For control week, x = 0 to 15; for intervention week, x = 0 to 8 (where intervention week 1 corresponds to control week 8).VO 2max , maximal oxygen uptake;