Assessing the cardiac autonomic response to bicycle exercise in Olympic athletes with different loads of endurance training: new insights from statistical indicators based on multilevel exploratory factor analysis

Aim: The mechanisms governing the organism’s response to exercise are complex and difficult to study. Spectral analysis of heart rate variability (HRV) could represent a convenient methodology for studying humans’ autonomic nervous system (ANS). However, difficulties in interpreting the multitude of correlated HRV-derived indices, mainly when computed over different time segments, may represent a barrier to its usage. This preliminary investigation addressed to elite athletes proposes a novel method describing the cardiac autonomic response to exercise based on multilevel exploratory factor analysis (MEFA), which reduces the multitude of HRV-derived indices to fewer uncorrelated ANS indicators capable of accounting for their interrelationships and overcoming the above difficulties. Methods: The study involved 30 Italian Olympic athletes, divided into 15 cyclists (prevalent high-intensity endurance training) and 15 shooters (prevalent technical training with low-intensity endurance component). All athletes underwent a complete test of a dynamic protocol, constituted by a rest-stand test followed by a stepwise bicycle stress test subdivided into a single bout of progressive endurance (from aerobic to anaerobic) exercise and recovery. Then, by spectral analysis, values of 12 ANS proxies were computed at each time segment (9 epochs in all) of the complete test. Results: We obtained two global ANS indicators (amplitude and frequency), expressing the athletes’ overall autonomic response to the complete test, and three dynamic ANS indicators (amplitude, signal self-similarity, and oscillatory), describing the principal dynamics over time of the variability of RR interval (RRV). Globally, cyclists have significantly higher amplitude levels (median ± MAD: cyclists 69.9 ± 20.5; shooters 37.2 ± 19.4) and lower frequency levels (median ± MAD: cyclists 37.4 ± 14.8; shooters 78.2 ± 10.2) than shooters, i.e., a parasympathetic predominance compared to shooters. Regarding the RRV dynamics, the signal self-similarity and oscillatory indicators have the strongest sensitivity in detecting the rest-stand change; the amplitude indicator is highly effective in detecting the athletes’ autonomic changes in the exercise fraction; the amplitude and oscillatory indicators present significant differences between cyclists and shooters in specific test epochs. Conclusion: This MEFA application permits a more straightforward representation of the complexity characterizing ANS modulation during exercise, simplifying the interpretation of the HRV-derived indices and facilitating the possible real-life use of this non-invasive methodology.