Fractal Correlation Properties of Heart Rate Variability as a Biomarker for Intensity Distribution and Training Prescription in Endurance Exercise: An Update

Background The non-linear index alpha 1 of Detrended Fluctuation Analysis (DFA a1) of heart rate variability, has been shown to be a marker of fatigue during endurance exercise. This report aims to explore its ability to assess the physiological status as a surrogate metric for “readiness to train” while performing simulated warm-up sessions the day after two different exercise sessions. Methods 11 triathletes were recruited to determine the first ventilatory threshold (VT1) during a baseline assessment and to perform 10-min of cycling at 90% of VT1 (simulating a warm-up bout) before (PRE) and within 36 h after (POST) light and heavy running exercise. RR intervals were recorded for DFA a1 analysis along with neuromuscular testing to verify the effects of the performed exercise sessions. In addition to common statistical methods, magnitude-based inferences (MBI) were applied to assess the changes in true score and thus also the practical relevance of the magnitude. Results Rating of perceived exertion for the heavy exercise session showed a significant higher rating as opposed to the light exercise session (p < 0.001, d = 0.89). In regard of MBIs, PRE versus POST comparisons revealed a significant reduced DFA a1 with large effect size after the heavy exercise session (p = 0.001, d = − 1.44) and a 99% chance that this negative change was clinically relevant. Conclusions Despite inter-individual differences, DFA a1 offers potential to assess physiological status and guide athletes in their training as an easy-to-apply monitoring procedure during a standardized warm-up. A regular assessment including individual data history and statistical references for identification of response is recommended. Further data are necessary to confirm the results in a larger and more homogeneous population.

Monitoring of the physiologic metric, respiratory frequency (RF), has been shown to be of value in health, disease, and exercise science. Both heart rate (HR) and variability (HRV), as represented by variation in RR interval timing, as well as analysis of ECG waveform variability, have shown potential in its measurement. Validation of RF accuracy using newer consumer hardware and software applications have been sparse. The intent of this report is to assess the precision of the RF derived using Kubios HRV Premium software version 3.5 with the Movesense Medical sensor single-channel ECG (MS ECG) and the Polar H10 (H10) HR monitor. Gas exchange data (GE), RR intervals (H10), and continuous ECG (MS ECG) were recorded from 21 participants performing an incremental cycling ramp to failure. Results showed high correlations between the reference GE and both the H10 (r = 0.85, SEE = 4.2) and MS ECG (r = 0.95, SEE = 2.6). Although median values were statistically different via Wilcoxon testing, adjusted median differences were clinically small for the H10 (RF about 1 breaths/min) and trivial for the MS ECG (RF about 0.1 breaths/min). ECG based measurement with the MS ECG showed reduced bias, limits of agreement (maximal bias, −2.0 breaths/min, maximal LoA, 6.1 to −10.0 breaths/min) compared to the H10 (maximal bias, −3.9 breaths/min, maximal LoA, 8.2 to −16.0 breaths/min). In conclusion, RF derived from the combination of the MS ECG sensor with Kubios HRV Premium software, tracked closely to the reference device through an exercise ramp, illustrates the potential for this system to be of practical usage during endurance exercise.

Identifying exercise intensity boundaries has been shown to be important during endurance training for performance enhancement and rehabilitation. Unfortunately, even though surrogate markers show promise when assessed on a group level, substantial deviation from gold standards can be present in each individual. The aim of this study was to evaluate whether combining two surrogate intensity markers improved this agreement. Electrocardiogram (ECG) and gas exchange data were obtained from 21 participants who performed an incremental cycling ramp to exhaustion and evaluated for first (VT1) and second (VT2) ventilatory thresholds, heart rate (HR) variability (HRV), and ECG derived respiratory frequency (EDR). HRV thresholds (HRVT) were based on the non-linear index a1 of a Detrended Fluctuation Analysis (DFA a1) and EDR thresholds (EDRT) upon the second derivative of the sixth-order polynomial of EDR over time. The average of HRVT and EDRT HR was set as the combined threshold (Combo). Mean VT1 was reached at a HR of 141 ± 15, HRVT1 at 152 ± 14 (p < 0.001), EDRT1 at 133 ± 12 (p < 0.001), and Combo1 at 140 ± 13 (p = 0.36) bpm with Pearson’s r of 0.83, 0.78, and 0.84, respectively, for comparisons to VT1. A Bland–Altman analysis showed mean biases of 8.3 ± 7.9, −8.3 ± 9.5, and −1.7 ± 8.3 bpm, respectively. A mean VT2 was reached at a HR of 165 ± 13, HRVT2 at 167 ± 10 (p = 0.89), EDRT2 at 164 ± 14 (p = 0.36), and Combo2 at 164 ± 13 (p = 0.59) bpm with Pearson’s r of 0.58, 0.95, and 0.94, respectively, for comparisons to VT2. A Bland–Altman analysis showed mean biases of −0.3 ± 8.9, −1.0 ± 4.6, and −0.6 ± 4.6 bpm, respectively. Both the DFA a1 and EDR intensity thresholds based on HR taken individually had moderate agreement to targets derived through gas exchange measurements. By combining both non-invasive approaches, there was improved correlation, reduced bias, and limits of agreement to the respective corresponding HRs at VT1 and VT2.