Is Baseline Cardiac Autonomic Modulation Related to Performance and Physiological Responses Following a Supramaximal Judo Test?

Little research exists concerning Heart Rate (HR) Variability (HRV) following supramaximal efforts focused on upper-body explosive strength-endurance. Since they may be very demanding, it seems of interest to analyse the relationship among performance, lactate and HR dynamics (i.e. HR, HRV and complexity) following them; as well as to know how baseline cardiac autonomic modulation mediates these relationships. The present study aimed to analyse associations between baseline and post-exercise HR dynamics following a supramaximal Judo test, and their relationship with lactate, in a sample of 22 highly-trained male judoists (20.70±4.56 years). A large association between the increase in HR from resting to exercise condition and performance suggests that individuals exerted a greater sympathetic response to achieve a better performance (Rating of Perceived Exertion: 20; post-exercise peak lactate: 11.57±2.24 mmol/L; 95.76±4.13 % of age-predicted HRmax). Athletes with higher vagal modulation and lower sympathetic modulation at rest achieved both a significant larger ∆HR and a faster post-exercise lactate removal. A enhanced resting parasympathetic modulation might be therefore related to a further usage of autonomic resources and a better immediate metabolic recovery during supramaximal exertions. Furthermore, analyses of variance displayed a persistent increase in α1 and a decrease in lnRMSSD along the 15 min of recovery, which are indicative of a diminished vagal modulation together with a sympathovagal balance leaning to sympathetic domination. Eventually, time-domain indices (lnRMSSD) showed no lactate correlations, while nonlinear indices (α1 and lnSaEn) appeared to be moderate to strongly correlated with it, thus pointing to shared mechanisms between neuroautonomic and metabolic regulation.


Introduction
During the last decade, Heart Rate (HR) Variability (HRV) has been widely used as a non-invasive marker of the Autonomic Nervous System (ANS) regulation of HR dynamics following exercises of different duration and intensities. Recently, new studies have been conducted to analyse cardiac autonomic regulation following exercise intensities described as supramaximal, since interval protocols (i.e. repeated bouts of brief supramaximal exercise,~15-30s, interspersed with recovery periods of similar duration) appear to be extremely effective in already highly-trained endurance athletes [1,2]; also in combat sports training [3].
Buchheit et al. [4] formerly showed a significant slower parasympathetic reactivation after a Repeated Sprint exercise compared with a continuous one of equivalent net caloric expenditure, concluding that anaerobic metabolism participation, rather than energy expenditure, modulates the level of post-exercise parasympathetic reactivation. A similar depressed post-exercise cardiac autonomic response has been reported by other authors following the anaerobic Wingate Test [5][6][7][8], Repeated Shuttle Sprints [9], or varying protocols of interval running exercises [10], suggesting that the depression of parasympathetic HRV indices (i.e., high frequency power, HF; and root-mean-square difference of successive normal R-R intervals, rMSSD) appears to be a common response to exercises with a large anaerobic component. Since all of the aforementioned studies have employed efforts mainly focused on lower limb, little research exists concerning cardiac autonomic regulation following supramaximal efforts, thus anaerobic, requiring high both neuromuscular and cardiovascular demands on the upper body; which it is a characteristic profile of some fighting sports like Judo [11][12][13]. Although these sports require the concomitant metabolic involvement of both anaerobic and aerobic pathways, the anaerobic capacity, focused in the upper body, points to be the prevalent and determining capacity due to the explosive throwing actions, as well as the gripping and grappling continuous actions in the contest [12,13].
On the other hand, it has been described higher lactate peaks in top athletes as an adaptive response to high intensity training [14,15]; but also a better lactate removal ability; this last due to an enhanced maximal muscle oxidative capacity and a greater content of monocarboxylate transporters [16,17]. This higher lactate clearance capacity, related in turn to a delayed appearance of muscle fatigue [18,19], might become a confounding factor in the assessment of the anaerobic capacity in the above cited multi-acceleration sport tasks [20,21]. However, despite this clearance capacity, the resulting stress metabolite accumulation importantly influences post-exercise sympathetic inhibition timeline through chemoreflex control of HR [22,23], thus interrelating metabolic and cardiac autonomic recovery processes.
To our knowledge, no previous study had undertaken a comprehensive analysis integrating: baseline cardiac autonomic modulation, HR indices (i.e., maximal HR increase during exercise and HR decay during recovery), performance, post-exercise blood lactate (i.e., peak accumulation and clearance) and nonlinear HR dynamics during the immediate recovery from an intervallic supramaximal exertion, mainly focused in upper limb strength. Given this scenario, the Blasco Specific Judo Test [BSJT;12], a Judo test which includes upper limb strength exercises, offers a validated [24] and suitable proposal to assess the relationship between cardiac autonomic control and lactate dynamics following intervallic supramaximal exercises primarily aimed at upper body, which have not yet been addressed in HRV studies.
Therefore, the present study aims: (1) to examine possible associations between: short-term baseline measures of HRV, performance, and both HR and blood lactate responses during and following a supramaximal judo test; (2) and to compare HRV and lactate dynamics during its immediate recovery. Following Nunan [25], we hypothesize that post-exercise HR dynamics might be related to baseline cardiac autonomic regulation. We also hypothesize an important mediating role for anaerobic metabolism handling (i.e., lactate accumulation and clearance) in post-exercise nonlinear HR dynamics.

Participants and experimental protocol
Twenty two male judokas from under-60 to under-81 kg categories (20.70 ± 4.56 years, 1.73 ± 0.06 m, 68.25 ± 7.17 kg, 22.84 ± 1.64 kg*m -1 ), took part in the study. They all had 3 to 7 years of competitive experience at senior national and international level and trained for 13 ± 6 hours a week. The study conformed to the Declaration of Helsinki and participants gave voluntary written consent to participate in the investigation, approved by the ethics committee of the University of Valencia. In the event of underage athletes, coaches signed a consent statement indicating that they have been advised of the potential dangers of participation and that they had understood them. The underage written consent also included the judokas' signature to show that they also acknowledged of these terms [26].
Special care was taken to standardize tapering prior to the test. Vigorous exercise was not allowed for 48 h before the testing day, and no training was permitted for 24 h before the testing day. Participants were summoned at 8:00 am on a fasting condition. On arrival to the gymnasium, they rested quietly for 10 min in the seated position to obtain baseline heart-rate derived data. They were allowed to breathe spontaneously but reminded to avoid irregular respiration. Breakfast was instructed to be appropriate for the subsequent exertion, and 2 h later a standardized warm-up was conducted to prepare participants for BSJT. Immediately after the test (<5 s), the judokas were assisted to sat passively on a chair placed adjacent to the tatami and remained still for the following 15 min. Sitting position and spontaneous breathing were chosen in order to match baseline and post-exercise conditions, as lying down position and controlled breathing would have been impractical during recovery. Moreover, normal respiratory rate does not result in significantly different heart rate-derived indices compared with controlled breathing [27].

Performance assessment: The BSJT
The assessment of sport performance and the associated physiological profile remains under debate in most complex anaerobic sports [28], since activity patterns are very changing and considerably influenced by the opponent features [21]. Besides, there is the issue of the specific nature of the sport and the ecological validity of some classical laboratory variables such as the maximum oxygen consumption attained in a treadmill or a cycloergometer maximal graded test [29]. In this context, since Judo contest cannot always be considered as a reference of maximum exertion, BSJT is aimed at reproducing both neuromuscular and energy demands of judo competition within a time structure similar to that of a contest, in order to assess the judoist's maximum performance and the use of specific skills under fatigue conditions. Time structure and the drills employed in the test (pull-ups, gripping up and down a rope, and four Judo specific tasks) have been previously described in detail, and the test have already shown to be valid (with significant performance differences between elite and national level judoist; n=53, 29 male and 24 female; p<0.05 for both, male and female athletes) [12]; and reliable (Test-Retest Intraclass Correlation 0.875; p = 0.003) [24].
Briefly, BSJT comprises six 15-s exercise periods, separated by either 10-s or 15-s recovery periods, repeated for three consecutive bouts (i.e., eighteen periods of 15-s maximum intensity); then a larger 30-s recovery period and a final "timeto-exhaustion" isometric pull-up ( TE PU, see Table 1 for a glossary of abbreviations). This final exercise assesses judoist's upper limb strength reserve (i.e., residual isometric strength of pulling muscles) following a series of strenuous explosive-strength efforts; a strength reserve that may become crucial in the Golden Score extra-time of Judo contests. Accordingly, total duration of the test varies from 8 to 9 minutes, depending on each judoist's capacity to maintain the TE PU. Because the BSJT assesses some specific judo drills, three judokas of similar body weights are needed to perform the test: 1 participant is evaluated and 2 other individuals receive throws. Finally, BSJT offers seven different performance assessment sub-scores: whole number of repetitions performed in each of the six drills and time achieved in the final TE PU, which after being converted into t-scores, are added up to generate a global performance score (BSJT score ).

HR, HRV, RPE and blood lactate concentrations measurement
HR and Rating of Perceived Exertion (RPE, scale 6-20) [30] were obtained immediately after the BSJT as measures of the relative intensity achieved by each subject. Final HR was retained for analysis as an absolute value (HR final ) and as a percentage of the age-predicted maximum HR (%HR final ). Recordings were performed using a Polar RS800 HR monitor . Baseline data were obtained from the last 5min segment of the resting recording, whereas recovery data were extracted on four 180-s windows (Post1, Post2, Post3 and Post4). The initial 180 s were not evaluated to minimize the influence of pronounced downward trends [31]. 5-min segments during resting recordings were employed to equate as much as possible the number of data points analyzed in resting and recovery conditions. Immediately before the warm-up, and 1, 3, 5, 10 and 15 minutes following BSJT, an earlobe capillary blood sample was collected for lactate concentration analysis ( 15 ) using a Lactate Pro LT-1710 analyser (Arkray Inc, Kyoto, Japan) [32]. The maximum value reached during the recovery period was defined as the peak concentration ([La] peak ), whereas [La] cle was calculated by subtracting lactate concentration at the 15 th minute of recovery from [La] peak , see Figure 1.

HR and HRV data analysis
The absolute change between baseline HR (HR baseline ) and maximum HR achieved during the test was defined as ∆HR; whereas HRR 1 and HRR 2 were calculated by subtracting HR recorded following 1 and 2 min of recovery from HR final . In these latter cases, HR was obtained as a mean value of a 5-s window.
The whole analysis process was carried out by the same researcher to ensure consistency. Artifacts were identified and corrected, and resting recordings were also detrended (Smooth priors, ƛ = 500) according to manufacturer's recommendations. Regarding linear analysis during resting condition, the standard deviation of normal RR intervals (SDNN) and rMSSD were calculated for the time domain. Low frequency power (LF, 0.04-0.15 Hz) and High Frequency (HF 0.15-0.4 Hz) were calculated as integrals of the respective power spectral density curve, using a Fast Fourier Transform algorithm; and LF/HF ratio was also retained for statistical analysis.
Examination of nonlinear HR dynamics was completed using measures of fractal scaling properties and complexity. Detrended Fluctuation Analysis (DFA) technique was applied to quantify self-similarity correlations [33]. Based on previous research and because of our relatively short recording time, only the short-term (4 to 11 beats) scaling exponent (α 1 ) was retained for statistical analysis. Eventually, Sample Entropy (SaEn) was calculated to provide a general indication of predictability of the time-series [34]. To characterize the stringency of match recognition, the length (m) of the subseries and the tolerance (r) of the matches were fixed at m = 2 and r = 20% of the SD of the datasets.
Later on, during recovery condition, only rMSSD was considered within linear indices, in order to measure the vagal reactivation [35], while frequency-domain indices were discarded given that its reliability is worsened following supramaximal efforts [36]. Nonlinear indices completed the analysis: α 1 is considered to estimate sympathovagal balance [37], and SaEn to provide an indication of the complexity of the time-series under these circumstances [34].

Statistical Analyses
Statistical analyses were carried out using the Statistical Package for the Social Sciences software (SPSS version 15.0, SPSS Inc., Chicago, USA). After testing for normal distribution (Kolmogorov-Smirnov test, with Lilliefor's correction), SDNN, RMSDD, LF, HF, LF/HF and SaEn were logarithmically transformed to allow parametric comparisons. Those variables will henceforth be referred as lnSDNN, lnRMSSD, lnHF, lnLF, lnLFHF and lnSaEn respectively.
Possible relationships between HR baseline , HRV indices and BSJT score , heart rate measures (∆HR, HRR 1 , HRR 2 ), and postexercise lactate concentrations ([La] peak, [La] cle ), were initially examined using a Pearson correlation. Besides, considering the previously described strong dependence of resting HRV indices on age [38] and intrinsic HR [39], and also on total spectral power in the case of linear measures [31], this analysis was re-conducted using a partial correlation controlling for either these three or two variables. A second Pearson correlation analysis was later performed to assess possible associations between BSJT score , heart rate measures (∆HR, HRR 1 , HRR 2 ) and post-exercise lactate concentrations ([La] peak, [La] cle ). Eventually, a third partial correlation analysis, now controlling for age, was performed to examine possible relationships between lnRMSSD, α 1  Finally, a repeated measures multivariate analysis of variance, controlling for age and followed by Bonferroni posthoc test, was used to assess the effect of the effort on lnRMSSD, lnSaEn and α 1 (i.e. Baseline, Post1, Post2, Post3, and Post4) and lactate (i.e. Baseline, 3, 5, 10 and 15 minutes after the test). Whenever Mauchly's Sphericity test was violated, necessary technical corrections were performed using the Greenhouse-Geisser test. Partial estimated effect size (partial η2) was reported as a measure of the magnitude of the effect. The significance level was set at p<0.05. However, trends towards signification (p<0.1) were also named [41]. Data are presented as means and standard deviations (±SD).

Results
Artifact beats accounted for less than 1% in the 22 participants presented in the study. Descriptive data of baseline HRV, performance, RPE at the end of the test, HR indices and lactate concentration measures are presented in Table 2.

General Discussion
Previous to the discussion of HRV and lactate results, it was important to ensure that the test was indeed a supramaximal upper limb effort. The values of RPE and HR at the end of the test confirmed this statement.
[LA] peak (11.57 ± 2.24 mmol/L) was close to the higher lactate productions described in previous studies based on judo combat simulation, i.e. 10.3 ± 2.6 and 9.2 ± 2.0 mmol/L for male and female Italian Olympic judoist [13], 10 mmol/L [42], or 12.3 ± 0.8 mmol/L [43]; and just a bit under lactate productions reached in real competition (13.4 ± 5.2 mmol/L in the first contest in an international judo tournament) [44]. However, although [LA] peak during recovery from the BSJT test was also similar to those previously reported following usual running supramaximal intermittent protocols [4,9], it seemed slightly lower compared to those measured after repeated Wingate test [45,46]. Considering that the abovementioned studies were carried out with recreational athletes, which are known to reach lower [LA] peak levels after supramaximal exercises [15,46], our lactate results seem even lower. A plausible explanation may lie in the fact that, even among specifically trained upper-limb athletes (i.e., judokas), it has been reported that the ability to perform anaerobic work with their upper body is only 81% of their lower limbs homologous capacity [47]. Additionally, it is plausible that during the 10 and 15-s recovery periods interspersed in the test, our athletes might have removed more lactate through aerobic pathways [17], given their higher fitness level (i.e., compared to recreational athletes).

Performance, lactate, HR response and baseline cardiac autonomic modulation
In this anaerobic context, one of the major outcomes of the present study was that those athletes with higher vagal modulation (i.e., higher lnHF) and lower sympathetic modulation (i.e., lower HR) before the test, were capable of increasing their HR to a greater extent during the exercise protocol, which it is considered a measure of maximal sympathetic response [48]. This result is in contrast with those obtained by Nunan et al. [25], who found no association between baseline HRV and ∆HR. Nunan et al. [25] argued that resting HRV was only related to those HR responses that are vagally mediated (i.e., HR recovery indices). However, our findings are in accordance with Leeper et al. [48], who postulated that the absolute change from baseline HR to peak HR during a maximal exercise was connected to the level of basal vagal modulation.
The association found between BSJT score and ∆HR may be justified by the need to exert at the maximum individual capacity, thus the maximum sympathetic response, to optimize the performance in such a high-intensity intervallic test. Moreover, according to Halliwill et al. [49], the final isometric pull-up until exhaustion would have increased "the pressor response mediated by progressive, parallel increases in muscle sympathetic nerve activity and vascular resistance". Similarly, resting lnRMSSD trend to correlation with BSJT score may reinforce the idea that a higher baseline vagal modulation would allow further usage of autonomic resources for judokas' responses to a supramaximal exercise, probably enabling a better performance during this kind of efforts. This postulate would be in accordance with the "autonomic resource hypothesis" proposed by Hynynen et al. [50]. These authors, by examining overtrained and non-overtrained athletes HRV response to a cognitive task, found that overtrained athletes were incapable of decreasing their HRV during a mental challenge (i.e., in relation to his individual resting values), unlike their non-overtrained counterparts. This absence of change (~response) was related to a worse performance in the task. Furthermore, since α 1 matched lnRMSSD associations, it may be suggested that a sympathovagal balance leaning to vagal dominance (i.e., a lower α 1 ) during resting conditions may contribute to a higher capacity to manage stress responses, and thus a better performance.
A relevant issue may be the differences in time-domain (i.e. lnRMSSD) and frequency-domain indices (i.e. lnHF) when controlling for Age, RRi and lnTP. lnHF kept its significance in the partial correlation with ∆HR, which is also a cardiac index, while lost it regarding BSJT score ., a measure of performance which is also mediated by factors unrelated to the cardiovascular functioning. Instead, lnRMSSD kept the trend to correlate with Performance, but lost its correlation regarding ∆HR, thus suggesting that the variables of control could be responsible for the pearson correlation with no controlls. Both, lnHF and lnRMSSD are vagal-mediated and intend to to measure the respiratory sinus arrhythmia [51] , but according to our results, the fast inter-beat changes in cardiac cycle measured by lnRMSSD might be better related to the overall factors that integrate performance. Although not exclusively, recent investigations point that power spectral analysis might provide a means to evaluate the ability to modulate autonomic outflows via baroreflexes [51,52] rather than the autonomic tone per se [51], what could explain some of these differences.
At the same time, and opposite to previous studies which found no relationship between baseline HRV and HRR indices [53][54][55] or even a positive association between resting parasympathetic activity and HR recovery post-exercise [25], the results of the present study have shown an inverse relationship between overall HRV (i.e., lnSDNN) at rest and HR recovery (i.e., HRR 2 ) and a positive association among baseline and recovery HR indices following a supramaximal exertion. This discordance may be attributed to differences in the exercise protocols employed, since the individualized graded exercise protocols used in previous studies require a similar relative workload for each subject [48], whereas during supramaximal interval exercises, a better/worse handling of peripheral fatigue importantly mediates the timing of performance decay [18,19]. Namely, the better the capacity to maintain performance at high intensities (i.e., greater explosive strength-endurance), the greater the relative workload accumulated during a supramaximal interval exercise. Thereby, it seems that our inverse correlation between overall HRV at rest and HR recovery indices together with the positive association among baseline and recovery HR indices were highly mediated by the different relative workload exerted by each judoka, as BSJT score was inversely correlated with HRR 2 . Large lactate productions during the test might explain the delay in recovery processes and thus the absence of correlations in regard to HRR 1 Another novel finding of our study was the relationship found between baseline HRV and post-exercise lactate removal after controlling for the proper confounding variables. Those judokas who presented greater baseline vagal activity (i.e., higher lnHF) and a sympathovagal balance leaning to vagal dominance (i.e., lower α 1 and lnLFHF) achieve a larger [LA] cle during recovery. Therefore, it may be suggested that an enhanced vagal modulation and HR dynamics unpredictability during resting conditions could be related to improved lactate removal ability following high-intensity exercises, since Thomas et al. [18,19] demonstrated that lactate removal ability during recovery from supramaximal efforts keeps relationship with both the maximal oxidative capacity and the content of monocarboxylate transporters. Notwithstanding, further studies including biochemical measures are required to clarify whether this apparently improved metabolic capacity in those participants with an enhanced vagal modulation is related or not to better lactate removal mechanisms.

Post-exercise HR dynamics and lactate
Our results displayed a persistent significant increase in α 1 and a decrease in lnRMSSD following BSJT (i.e., compared to baseline condition), which are indicative of diminished vagal modulation together with increased sympathetic modulation. These results are similar to previous studies assessing recovery from Wingate test [7], cycle ergometer incremental maximal exercise [56] and isometric handgrip [57]. The first two coincide with us in the anaerobic character of the effort, whereas we share with the latter the isometric handgrip ability. The increase in lnRMSSD from Post2 to Post3 condition shows a vagal modulation reactivation which, however, does not achieve returning cardiac parasympathetic activity to resting values. An unexpected absence of change in lnSaEn, from resting to post-exercise condition, suggest that HR complexity is affected to a different extent by a supramaximal effort, at least among highly-trained athletes, since previous studies involving non-athletes did show a decrease in lnSaEn during recovery [5,7]. Notwithstanding, these studies were conducted on very short lower-limb anaerobic efforts (i.e. wingate test), so further investigation is needed to clarify whether HR complexity following high-intensity exercise behaves differently among athletes and non-athletes, or it is related to differences in the protocol employed..
On the other hand, previous studies have also shown significant relationships between parasympathetic reactivation (both by means of HR recovery and HRV indices) and postexercise blood lactate concentration [58,59]. Besides, recently Simoes et al. [60] have shown a concordance between HRV and lactate regarding anaerobic threshold determination. Our results corroborate these previous investigations [4,58,59] in relation to the fact that initial parasympathetic reactivation (i.e. HRR 1 and HRR 2 ) following exercises with a large anaerobic component is hindered by a delay in the inhibition of sympathetic activity, due to lactate effect on chemoreflex control of HR [22,23]. Actually, both HRR 1 and HRR 2 were strongly correlated with post-exercise [LA] peak . Furthermore, some significant, moderate to strong correlations found between lactate concentration and both SaEn and α 1 during recovery (mainly considering the time course of the assessments; see the shadowed areas in Table 5) suggest that sympathovagal balance and HR complexity following a supramaximal effort may be also influenced by lactate concentration. According to these correlations, those judokas with higher lactate concentrations in our study, exhibited a sympathovagal balance leaning to vagal dominance, jointly with more complex HR dynamics. This is in agreement with Tulppo et al. [61], who found that subjects who have a higher sympathetic activity following an exercise also exhibited an augmented vagal activity (i.e., post-exercise muscle sympathetic nervous activity was inversely associated with LF/HF ratio). Eventually, such greater interplay between vagal and sympathetic branches of the ANS would result in more complex HR dynamics.
Finally, unlike Obmiński et al. [44], who described a lactate clearance parabola following official Judo matches (2 to 5 min), our results displayed a lactate nearly steady-state, similarly to Bastos et al. [62]. As in this latter study, our athletes remained seated during the whole recovery period in order to match baseline conditions in regard to body position because this factor is known to significantly affect HR dynamics [63]. This seating position probably prevented our athletes from achieving a larger lactate removal [64], although the significant difference between [LA] 15 and [LA] 5 seems to indicate that athletes began to effectively remove lactate from blood from this moment on; unfortunately, we could not confirm this hypothesis as no further lactate samples were collected.

Limitations
The short follow-up and the absence of a control exercise with whom to compare HRV response obtained following BSJT makes necessary further studies that tackle with these issues. Another limitation of this study may be the employment of Polar HR monitors instead of ECG recordings to estimate RR intervals, as possible physiological artifacts, such as ectopic beats or arrhythmic events, are more difficult to detect. Notwithstanding, this instrument has been previously validated for the measurement of RR intervals and subsequent HR dynamics analysis [65,66]. Actually, Wallén et al. [67] recently reported that the agreement between RS800 HR monitor and electrocardiographic recordings range from good to excellent among men, thus reinforcing previous results from Weippert et al. [68]. Finally, spontaneous breathing may be another limitation, although this is a persistent limitation in those fieldtest studies following supramaximal efforts. Moreover, recent studies [69] have shown that there are no differences between spontaneous breathing and controlled breathing in the longitudinal follow-up of athletes, at least for RMSSD, which is an important index in our study.

Conclusions
In summary, the major outcomes of this investigation were: (1) In highly-trained athletes, higher baseline vagal modulation and lower sympathetic activity prior to a supramaximal exertion is related to a greater response capacity during exercise, as evidenced by a larger ∆HR. Additionally, lower baseline sympathetic activity is also related to a better performance, and a clear positive tendency towards significance is present between resting vagal modulation and performance. These results endorse previous assumptions from Leeper et al. [48], and extend Hynynen et al. [50] related Autonomic Resource Hypothesis to supramaximal interval efforts (2). Following supramaximal exertions, initial HR decay is mediated by preceding individual relative workload limiting their use in interindividual comparisons regarding parasympathetic reactivation capacity (3). An enhanced resting vagal modulation and a sympathovagal balance leaning to vagal dominance (i.e. lower α 1 ) before performing a high-intensity exercise is related to an improved post-exercise lactate removal capacity. Future studies including biochemical measures should clarify whether this apparently improved metabolic capacity is related or not to better lactate removal mechanisms (4). Larger lactate concentrations after supramaximal upper limb explosive strength exercises are related, in our study, to a slower initial HR decay followed by a delayed sympathovagal balance leaning to vagal dominance and a more complex HR dynamics in highly trained judokas. Following Tulppo [61], long-lasting sympathetic hyperactivation due to lactate accumulation might have hindered initial HR recovery but subsequently enhanced cardiac sympathovagal outflow and complexity of HR dynamics.