Jakub S. Gąsior1*
Maciej Rosoł2
Marcel Młyńczak2
Andrew A. Flatt3
Bartosz Hoffmann4
Rafał Baranowski5
Bożena Werner1- 1Department of Pediatric Cardiology and General Pediatrics, Medical University of Warsaw, Warsaw, Poland
- 2Faculty of Mechatronics, Institute of Metrology and Biomedical Engineering, Warsaw University of Technology, Warsaw, Poland
- 3Biodynamics and Human Performance Center, Department of Health Sciences and Kinesiology, Georgia Southern University (Armstrong Campus), Savannah, GA, United States
- 4Physiotherapy Division, Faculty of Medical Sciences, Medical University of Warsaw, Warsaw, Poland
- 5Department of Heart Rhythm Disorders, National Institute of Cardiology, Warsaw, Poland
Background and Purpose: Most studies on heart rate variability (HRV) in professional athletes concerned linear, time-, and frequency-domain indices, and there is lack of studies on non-linear parameters in this group. The study aimed to determine the inter-day reliability, and group-related and individual changes of short-term symbolic dynamics (SymDyn) measures during sympathetic nervous system activity (SNSa) stimulation among elite modern pentathletes.
Methods: Short-term electrocardiographic recordings were performed in stable measurement conditions with a 7-day interval between tests. SNSa stimulation via isometric handgrip strength test was conducted on the second day of study. The occurrence rate of patterns without variations (0V), with one variation (1V), two like (2LV), and two unlike variations (2UV) obtained using three approaches (the Max–min, the σ, and the Equal-probability methods) were analyzed. Relative and absolute reliability were evaluated.
Results: All SymDyn indices obtained using the Max–min method, 0V, and 2UV obtained using the σ method, 2UV obtained using the Equal-probability method presented acceptable inter-day reliability (the intraclass correlation coefficient between.91 and.99, Cohen’s d between −.08 and.10, the within-subject coefficient of variation between 4% and 22%). 2LV, 2UV, and 0V obtained using the Max–min and σ methods significantly decreased and increased, respectively, during SNSa stimulation—such changes were noted for all athletes. There was no significant association between differences in SymDyn parameters and respiratory rate in stable conditions and while comparing stable conditions and SNSa stimulation.
Conclusion: SymDyn indices may be used as reliable non-respiratory-associated parameters in laboratory settings to detect autonomic nervous system (ANS) activity modulations in elite endurance athletes. These findings provide a potential solution for addressing the confounding influence of respiration frequency on HRV-derived inferences of cardiac autonomic function. For this reason, SymDyn may prove to be preferable for field-based monitoring where measurements are unsupervised.
Introduction
Comprehensive monitoring and identification of the physiological state of athletes by coaches, exercise scientists, and/or sports physicians can be accomplished using sensitive, non-invasive, time-efficient, and cost-effective testing methods and biomarkers. Knowledge of status is important for optimizing training adaptations and improving performance. Status markers may also aid in the diagnosis and risk prediction of medical conditions, such as sport-related concussion in elite athletes (Meeusen et al., 2013; Bourdon et al., 2017; Gabbett et al., 2017; Bishop et al., 2018; Heidari et al., 2018; Kellmann et al., 2018; Schneider et al., 2018). To distinguish between intended (e.g., training-related, lifestyle-related, or injury, e.g., concussion-induced) and unintended (measurement error) changes in selected physiological parameters, it is crucial to use objective and reliable measurements with validated tools (Lachin, 2004; Matheson, 2019).
Heart rate variability (HRV) parameters are becoming increasingly popular as non-invasive, reliable and sensitive biomarkers reflecting changes in autonomic nervous system (ANS) activity in athletes (Baumert et al., 2006; Buchheit, 2014; Flatt and Esco, 2016; Pereira et al., 2016; Esco et al., 2018; Schneider et al., 2018; Vescovi, 2019; Perrone et al., 2021). A vast majority of studies on HRV in professional athletes concerned only linear (time- and frequency-domain) parameters, which may not provide comprehensive description of ANS activity within this populations. Importantly, selected linear HRV parameters are inadequate to assess autonomic control and activity during short-time series commonly observed in exercise physiology or sports medicine (Gronwald et al., 2020; Rogers et al., 2021a,b). The signal stationarity and “controlled” respiratory rate are required for the short-term frequency-domain HRV analysis (Magagnin et al., 2011; Quintana et al., 2016). Therefore, in studies where the immediate effect is measured or specific exercise task influences breathing pattern, these measures are inappropriate. In recent years, measures of the non-linear dynamics of HRV have provided new opportunities to monitor cardiac autonomic regulation (Voss et al., 2009; Sassi et al., 2015; Henriques et al., 2020). It is known that cardiovascular control involves non-linear interactions between physiological systems, with complex dynamics (Schiecke et al., 2018; Saul and Valenza, 2021). Considering that many biomarkers are non-stationary signals (Goldberger, 1991; Faust et al., 2004) and vary in a complex and non-linear way (Voss et al., 2009), some authors highlight that non-linear HRV indices are more suitable than linear for the evaluation of individual changes in test–retest studies (Maestri et al., 2007) and underline the need to perform non-linear analysis to provide holistic information on HRV (Acharya et al., 2006; Huikuri et al., 2009; Voss et al., 2009; Hoshi et al., 2013; Sassi et al., 2015).
Recently (2020), we found that short-term (5-min) non-linear indices of entropy measures (Approximate Entropy and Sample Entropy), Poincaré plot (SD2/SD1), and index of HRV based on fractal correlation properties (short-term scaling exponent of detrended fluctuation analysis) presented large relative or high absolute test–retest reliability among elite endurance athletes (Hoffmann et al., 2020). Apart from these indices, symbolic dynamics have also gained wide acceptance for assessing various complex systems, but their applications in HRV studies are substantially less popular than the aforementioned non-linear parameters (Henriques et al., 2020).
The symbolic analysis method consists mainly of the transformation of a time series (RR intervals) into short patterns (three beats long), their classification, and the evaluation of their rates of occurrence (Porta et al., 2001). Porta et al. (2007a) showed that two non-linear symbolic indexes, that is, patterns with no variation (0V), and patterns with two unlike variations (2UV), represent a valid alternative to linear spectral indexes for assessment of the cardiac autonomic modulation from short-term heart period variability. The enhancement of cardiac sympathetic modulation and reduction of vagal modulation during sympathetic activation was associated with an increase of 0V and a decrease of 2UV (Porta et al., 2007b; Cysarz et al., 2013).
Very recently (2021), Matsumura et al. (2021) showed that sympathetic activity is dominant prior to competition attempts in elite athletes in a real-world sport event and concluded that sympathetic predominance before and/or during professional competition would be advantageous for performance in risky sports, such as snowboard jumping. Recognizing and modifying cardiac sympathetic and parasympathetic activity could be useful for coaches, exercise scientists and/or sports physicians in helping athletes to improve their physical performance and achieve better sport results.
Summarizing, little is known about the reliability of symbolic analysis when applied to short-term measurements in athletic populations (Gronwald and Hoos, 2020; Uhlig et al., 2020). HRV analysis has been mostly performed on the basis of group changes in endurance athletes (Baumert et al., 2006; Schmitt et al., 2018; Deus et al., 2019), which is ineffective in detecting individual athletes’ responses (Muñoz-López and Naranjo-Orellana, 2020). Greater attention to individual analysis seems to be necessary when assessing physiological responses to a specific training or interventions in athletes (Gąsior et al., 2020; Hoffmann et al., 2020). Cardiac autonomic modulation is linked to spontaneous respiration (Eckberg and Eckberg, 1982; Eckberg et al., 1985; Eckberg, 2003; Narkiewicz et al., 2006; Bari et al., 2016), hence, HRV (mostly frequency-domain) parameters are affected by respiratory depth and frequency (Hirsch and Bishop, 1981; Brown et al., 1993; De Souza et al., 2018). Sports field practitioners should interpret changes in HRV parameters with concomitant respiratory rate alterations. Therefore, the presented study has the following aims: (i) to assess the inter-day reliability of symbolic dynamics (SymDyn) indices obtained using three different transformations; (ii) to determine group and individual changes during SNSa stimulation in short-term SymDyn measures; and (iii) to quantify associations between differences in respiratory rate and SymDyn measures among elite modern pentathletes.
Materials and Methods
Study protocol with details is presented in Figure 1. Details of the study population (participants), study protocol, ECG acquisition, respiratory rate recordings, and SNSa stimulation have been presented elsewhere (Gąsior et al., 2020; Hoffmann et al., 2020). The study was approved by the University Ethical Committee (SKE 01-01/2017, 7 March 2017, Warsaw, Poland) and followed the rules and principles of the Helsinki Declaration. All athletes were informed of the aims and risks involved with the protocol and subsequently provided written informed consent prior to data collection.
Figure 1. Study protocol.
The ECGs were visually inspected for potential non-sinus or aberrant beats, and such erroneous beats were corrected from the cardiac interval series (RR series) before HRV analysis. The erroneous beats were manually corrected, that is, one R-R interval before, and one after each non-sinus beat were eliminated and replaced by R-R intervals computed by interpolation of degree zero based on the surrounding normal beats (Peltola, 2012).
Screened and corrected 5-min RR series were imported for analysis into PyBiOS software from ASCII text files (Silva et al., 2020). Imported RR series were not corrected for artifacts, segmented, filtered, or detrended using software tools.
Analysis of Symbolic Dynamics
In general, symbolic analysis consist of three steps: (1) conversion of the original series into a series of symbols; (2) definition of words (sequences of symbols) of interest; and (3) quantification of the different words, quantities that may be related to the phenomena of interest. The second and third steps in three transformation methods (Max–min method, σ method, and Equal-probability method) are similar. They differ in the way the signal of symbols is generated (Cysarz et al., 2018). Details concerning time series transformation into the symbol sequences using the appropriate alphabet for each method were precisely described by Wessel et al. (2007), Porta et al., (2007b) and Cysarz et al., (2013, 2018).
In this study, SymDyn indexes were obtained using three transformation methods: Max–min method, σ method, and Equal-probability method (Cysarz et al., 2013, 2018).
In Max–min method, the series of RR intervals were converted into a series of symbols through a uniform quantization of six levels (quantization level: 6). This means that six equal ranges were defined from the minimum to the maximum value within the series [l = (max(xi)—min(xi)]/ξ, quantization level ξ = 6), and each value in the original series was converted into a symbol (0–5) according to the range it belongs to Porta et al. (2001, 2007b).
In the σ method, three levels were defined using the following quantization lines: the signal average (μ), the signal average shifted up by a factor a, that is, (1 + a) μ, and the signal average shifted down by a factor a, that is, (1 − a) μ. The parameter a (sigma rate) was set to.05 (Kurths et al., 1995).
The Equal-probability method divides the full range of the signal into quantization levels ensuring that all levels will contain the same number of points. Thus, if the signal has length L, each level will have L/quantization level samples. If L is not a multiple of quantization level, the number of points can vary by one within the levels (Cysarz et al., 2018). The transformation was used with two different quantization levels: 4 (q = 4) and 6 (q = 6) enabling a direct comparison with the σ method and Max–min method, respectively.
For all the three methods, all sequences of three consecutive symbols (words) are classified into one of four families: 0V (zero variation), if the three symbols are equal. Examples of words from this family are {1,1,1} and {5,5,5}; 1V (one variation), which represent sequences with only one variation. The words {1,1,2} and {3,3,0} are examples from this family; 2LV (two like variation), representing sequences with two variations in the same direction, that is, the symbols are all different and form an increasing or decreasing ramp. Examples from this family are {0,3,5} and {2,1,0}; 2UV (two unlike variation), where symbols vary two times, in opposite directions, composing a peak or a valley. The words {1,2,0} and {3,0,3} are examples from this family. Finally, the percentages of words classified into each family are taken for analysis of the series dynamics (Porta et al., 2001, 2007a,b; Guzzetti et al., 2005; Cysarz et al., 2013, 2015, 2018).
Evidently, words from 0V family represent the slowest oscillations in the series, whereas words from 2UV represent the fastest. The 1V and 2LV families are intermediate levels. In case of HRV series, higher 0V% and 2UV% were associated with the higher the sympathetic and vagal modulation to the heart, respectively (Porta et al., 2001, 2007a,b, 2015; Guzzetti et al., 2005; Baumert et al., 2009; Tobaldini et al., 2009; Cysarz et al., 2013; Silva et al., 2017).
Statistical Analysis
For reliability assessment, the intraclass correlation coefficient (ICC), the within-subject coefficient of variation (WSCV) and Cohen’s d were calculated. ICC (Shrout and Fleiss, 1979) and WSCV (Bland and Altman, 1996; Atkinson and Nevill, 1998; Hopkins, 2000a,b) were used to analyze the relative reliability and absolute reliability respectively, while Cohen’s d was used to assess the effect size of the mean differences between repeated measurements. The interpretation of ICC values was determined a priori as follows: from 0 to.30—small, from.31 to.49—moderate, from.50 to.69—large, from.70 to.89 very large and from.90 to 1.00 the value was considered nearly perfect (Hopkins et al., 2009). Smaller WSCV means better reproducibility (Shoukri et al., 2008). The interpretation of Cohen’s d (Cohen, 1988) was also determined a priori, considering the value less than.20 as trivial, from.21 to.60 as small, from.61 to 1.20 as moderate, from 1.21 to 2.00 as large and greater than 2.01 as very large (Hopkins et al., 2009). Overall, if ICC was greater than.90 and Cohen’s d value trivial, this was interpreted as acceptable reliability. For visualization of the data, Bland–Altman plots were created with marked 95% limits of agreement (LoA) and maximum allowed difference (Abu-Arafeh et al., 2016), which was equal to the smallest worthwhile change (SWC) calculated as a.2 times SD of the Test measurements values (Buchheit, 2018).
In order to compare the results obtained from test and retest, Wilcoxon signed-rank tests were used (due to the small sample size). To assess the correlation between repeated measurements, Spearman’s rank correlation coefficient, Lin’s concordance coefficient, and bias correction factor (from Lin’s analysis), were calculated. Also, the dependency between the obtained results was evaluated by intercept from linear regression models. The ratio of the intercept from the model, divided by the mean value of the parameter from both measurements, enables the relative comparison of the results for individual parameters. The obtained value is associated with the bias for the consecutive measurements.
To compare values from Retest and SNSa stimulation, and hence the impact of the intervention on the parameter values, the mean and SD of the difference between measurements were calculated. Then, those values were divided by the mean value of the given parameter in both groups and presented as a percentage rate of change to enable the comparison of the influence of the intervention on the results. To compare the results from Retest and SNSa stimulation, Wilcoxon signed-rank tests were performed. Bland–Altman plots were also created for each comparison.
The changes in HRV parameter values between two measurements (test–retest and retest-SNSa) were studied in association with the changes in respiratory rates (identically like presented elsewhere Hoffmann et al., 2020; Gąsior et al., 2020). Correlation plots, along with linear regression model fitting were created for consecutive measurements with the results of Pearson correlation coefficients (r) for each of those associations. In this and other analyses, the change or difference of the parameters was studied by subtracting the values for Test from Retest, and values for Retest from SNSa.
Analysis based on group changes is ineffective in detecting individual athletes’ responses (Muñoz-López and Naranjo-Orellana, 2020). It was suggested that HRV be assessed in athletes on an individual basis (Plews et al., 2013). Therefore, individual athletes’ responses have been graphically presented in this study.
Statistical analysis and figures were, respectively, performed and prepared using Python 3.7.10 (including primarily modules: NumPy version 1.18.1, Pandas 1.0.1, SciPy 1.4.1 and scikit-learn version 0.22.1) and R 4.1.0 (including primarily packages: irr version 0.84.1 and effsize version 0.8.1). The significance level is.05.
Results
Results of four participants out of 12 were excluded due to diagnosis of prolonged QTc interval >450 ms (n = 2) or left bundle branch block (n = 2)—athletes were referred for further medical investigation. Results of eight male active elite modern Caucasian pentathletes living in Warsaw (Poland), medallists of the World Championships (n = 3), European Championships (n = 5), and National Championships (n = 8) were included in the statistical analysis. The mean (±SD) age, weight, height, body mass index (BMI), and duration of professional athletic career were: 21.7 years (±3.1), 75.9 kg (±9.5), 182.6 cm (±6.1), 22.7 kg/m2 (±2.3), and 10.8 years (±2.9). Athletes declared participating in 19 training sessions (±2) per week during the normal in-season time.
Reliability Analysis and Differences Between Test and Retest
Results for reliability analysis are presented in Table 1. For all SymDyn indices obtained using the Max–min method, 0V and 2UV obtained using the σ method, and 2UV obtained using the Equal-probability methods, ICC and Cohen’s d were nearly perfect and trivial, respectively. WSCV was between 2% and 12% for all SymDyn indices obtained using the Max–min method, between 8% and 25% using the σ-method, between 7% and 28% using the Equal-probability methods.
Table 1. Results of reliability analysis.
The Bland–Altman plots for SymDyn indices obtained using the Max–min method, the σ method, and the Equal-probability method q = 6, q = 4 presented in Figures 2, 3, respectively. The smallest mean difference between Test and Retest was observed for SymDyn indices obtained using the Max–min method. The maximum allowed differences (i.e., smallest worthwhile changes) for 0V, 1V, 2LV, and 2UV were: 1.9,.5, 1.5, 2.2 (the Max–min method); 2.6, 1.1, 1.0, 2.0 (the σ method); 1.3, 1.0, 1.5, 2.4 (the Equal-probability method q = 6); 2.2,.8, 1.2, 2.2 (the Equal-probability method q = 4), respectively. For all analyzed parameters, regardless of the method used, LoA exceeded the defined maximum allowed difference.
Figure 2. Bland–Altman plots for SymDyn indices: 0V, 1V, 2LV, and 2UV (Test and Retest): (A) the Max–min method, (B) the σ-method. The black dashed line indicates the mean difference, the grey dashed lines are the limits of agreement (±1.96 SD), and the green lines are the maximum allowed difference.
Figure 3. Bland–Altman plots for SymDyn indices: 0V, 1V, 2LV, and 2UV (Test and Retest): the Equal-probability method (A) q = 6, (B) q = 4. The black dashed line indicates the mean difference, the grey dashed lines are the limits of agreement (±1.96 SD), and the green lines are the maximum allowed difference.
There were no significant differences in all SymDyn indices obtained using all methods between Test and Retest. There was a significant correlation between values from Test and Retest for all indices obtained using the Max–min method, for 0V and 2UV obtained using the σ method and the Equal-probability method (q = 6), for 1V and 2UV obtained using the Equal-probability method (q = 4). Lin’s concordance coefficient was the highest for 2UV regardless of the method used. For the rest of SymDyn indices obtained using the Max–min method, Lin’s concordance coefficient was between.89 and.94. The smallest value of the intercept of the linear model scaled by the mean value of the parameter was observed for 2UV obtained using the Max–min method, while the highest for 1V obtained using the σ method (Table 2).
Table 2. Results of the comparison between Test and Retest in analyzed symbolic dynamics indices.
SNSa Stimulation (Impact of Intervention)
Statistically significant increases in 0V and decreases in both 2LV and 2UV with no changes in 1V were observed during SNSa stimulation, independently on the method used to obtain indices (Table 3). The largest percentage change was observed for 2LV. The Bland–Altman plots for SymDyn indices obtained using the Max–min method, the σ method and the Equal-probability method q = 6, q = 4 are presented in Figures 4, 5, respectively. The smallest mean difference between Retest and SNSa was observed for 1V obtained using all methods.
Table 3. Results of the comparison between Retest and SNSa stimulation in analyzed symbolic dynamics indices.
Figure 4. Bland–Altman plots for SymDyn indices 0V, 1V, 2LV, and 2UV (Retest and sympathetic nervous system activity—SNSa stimulation): (A) the Max–min method, (B) the σ-method. The black dashed line indicates the mean difference, the grey dashed lines are the limits of agreement (±1.96 SD), and the green lines are the maximum allowed difference.
Figure 5. Bland–Altman plots for SymDyn indices 0V, 1V, 2LV, and 2UV (Retest and SNSa stimulation): the Equal-probability method (A) q = 6, (B) q = 4. The black dashed line indicates the mean difference, the grey dashed lines are the limits of agreement (±1.96 SD), and the green lines are the maximum allowed difference.
Correlation Between Changes in Non-linear HRV Parameters and Changes in RespRate
The correlations for retest–test and SNSa stimulation–retest differences between SymDyn indices and RespRate are shown in Figures 6–9, respectively. The correlation coefficients were not significant for all indices for all methods.
Figure 6. Correlation between retest–test differences in SymDyn indices and differences in respiratory rate (RespRate): (A) the Max–min method, (B) the σ-method.
Figure 7. Correlation between retest–test differences in SymDyn indices and differences in respiratory rate (RespRate): the Equal-probability method (A) q = 6, (B) q = 4.
Figure 8. Correlation between SNSa stimulation–retest differences in SymDyn indices and differences in respiratory rate (RespRate): (A) the Max–min method, (B) the σ-method.
Figure 9. Correlation between SNSa stimulation–retest differences in SymDyn indices and differences in respiratory rate (RespRate): the Equal-probability method (A) q = 6, (B) q = 4.
Individual Changes
By analyzing individual HRV changes, professionals can identify athletes who show different responsiveness to the selected stressors (Mann et al., 2014; Vescovi, 2019; Muñoz-López and Naranjo-Orellana, 2020). SymDyn indices values for individuals during subsequent measurements are shown in Figure 10 for the Max–min method and σ method and in Figure 11 for the Equal-probability methods. SymDyn indices obtained using the Max–min method and the σ method presented a similar pattern of changes, that is, increase in 0V and decrease in 2LV, 2UV during SNSa stimulation were observed for all athletes (Figure 10). For indices obtained using the Equal-probability methods, there was one athlete (#6) who presented a lower value of 0V during SNSa stimulation than in Retest and one athlete (#5) who had a higher value of 2UV during SNSa than in Retest (but only for the Equal-probability method q = 6; Figure 11). In all other athletes, the same pattern was observed as described for the Max–min and the σ methods. The highest variations in individual changes were observed for 1V.
Figure 10. Individual changes in SymDyn indices between Test, Retest, and SNSa stimulation: (A) the Max–min method, (B) σ-method.
Figure 11. Individual changes in SymDyn indices between Test, Retest, and SNSa stimulation: the Equal-probability method (A) q = 6, (B) q = 4.
Discussion
All symbolic dynamics measures obtained using the Max–min method, 0V and 2UV obtained using the σ method, and 2UV obtained using the Equal-probability methods, calculated based on short-term (5-min) RR series derived from ECG recordings performed in stable measurement conditions with 1-week time interval between examinations presented acceptable inter-day reliability in elite modern pentathletes. We found that the rate of patterns with two variations (2LV and 2UV) and the frequency of patterns with no variations (0V) obtained using the Max–min and σ methods significantly decreased and increased, respectively, during SNSa stimulation. Such changes were noted for all athletes. The rate of patterns with one variation (1V) remained unchanged with different responses among athletes. There was no significant association between differences in SymDyn parameters and respiratory rate in stable conditions and while comparing stable conditions and SNSa stimulation.
To summarize findings from our research project, short-term (5-min) time-domain log transformed RMSSD (ICC = .98, Cohen’s d = .11, WSCV = 2%), frequency-domain log transformed HF (ICC = .86, Cohen’s d = .17, WSCV = 6%), descriptors of the Poincaré plot—SD2/SD1 (ICC = .87, Cohen’s d = −.07, WSCV = 9%; Hoffmann et al., 2020), 0V obtained using the Max–min and the σ method (ICC between.91 and.95, Cohen’s d between.03 and.06, WSCV between 12% and 22%) and 2UV obtained using all methods (ICC between.94 and.99, Cohen’s d between.02 and.10, WSCV between 4% and 18%) may be used by sport practitioners and researchers as reliable parameters in elite modern pentathletes to detect changes during baseline examinations in laboratory settings. Good absolute and substantial relative reliability for SymDyn variable (2UV) obtained from laboratory recordings was confirmed for healthy volunteers by Maestri et al. (2007).
Sport professionals should keep in mind that some athletes may present values during re-examination in the same stable conditions that are higher or lower than maximum acceptable difference (smallest worthwhile change). Recent reviews have underlined that important data from the Bland–Altman method are often omitted (Abu-Arafeh et al., 2016; Gerke, 2020). The definition of the a priori acceptable LoA, to define the minimal agreement between measurement was set as the first key item from comprehensive list of items published in Abu-Arafeh et al. (2016). In our study, all analyzed SymDyn indices showed a LoA that exceeded the defined a priori maximum acceptable difference (SWC).
Significant correlations between test–retest differences in respiratory rate and HF in elite endurance athletes were shown (Hoffmann et al., 2020). In presented results, there was no significant association between test–retest differences in respiratory rate and SymDyn indices. Non-linear HRV measures do not require a stationary signal (Michael et al., 2017). Lack of significant dependency on respiratory rate may be considered as another advantage of SymDyn indices. However, alterations in SymDyn parameters may be associated with respiratory depth. In some applications depth of breathing appeared more important for the relation with HRV parameters that its rate (Młyńczak and Krysztofiak, 2019). There is no state-of-the-art solution on how to record respiratory mechanics along with cardiac activity in HRV studies (Quintana and Heathers, 2014; De Souza et al., 2018). One possible solution would be to use Pneumonitors—portable devices measuring both depth and rate of breathing using impedance pneumography, and single-lead ECG (designed specifically for environmental physiology and sports medicine analysis; Młyńczak et al., 2014, 2017).
Porta et al. (2007b) and Cysarz et al. (2013) performed a head-up tilt test to assess cardiac autonomic regulation by symbolic indexes obtained from short-term recordings (250 cardiac beats; range 220–260) in healthy subjects. It was shown that SymDyn indexes are linearly associated with tilt angles during the test (the correlation was observed in more participants for 0V and 2UV than for any spectral parameters) and represent different aspects of the autonomic response to head-up tilt test. Authors suggested that symbolic indexes might be considered as an alternative to frequency-domain HRV parameters in scenarios where short-term recordings are analyzed (Porta et al., 2007b). One of the key differences between these two approaches is that spectral analysis is bounded to long sinusoid components, with different frequencies (in Hz), whereas the SymDyn decomposes the series into small sequences (three values), and the oscillatory patterns are not restricted to any frequency in Hz.
Very recently (2021) Storniolo et al. (2021) analyzed results of symbolic markers pre, post, and during a plateau phase—very short epoch (22 ± 8 cardiac beats, corresponding to 9 ± 3 s.) following maximal sprint exercise—in physically active subjects. Authors observed increased 2UV% and unvaried 0V% during plateau compared to pre sequences, suggesting a remarkable vagal modulation and persistency of the sympathetic control, respectively; 0V% significantly decreased and 2UV% was significantly higher in post compared to pre sequences, suggesting a less active sympathetic and more active vagal control in post than pre sequences, respectively (Storniolo et al., 2021). Authors underlined that such conclusions cannot be achieved from the LF/HF ratio parameter. Moreover, they suggested that future studies should evaluate the eventual association of the symbolic measures with the athlete’s performance (Storniolo et al., 2021).
In our study, to stimulate and assess cardiovascular autonomic function and responsiveness in elite endurance athletes, the isometric handgrip strength test was used (Khurana and Setty, 1996). The effect and procedure of squeezing the dynamometer’s handle may be, to some extent, considered as pre-competition stress and comparable in sports like modern pentathlon to shooting performance. During SNSa stimulation, not unexpectedly, 0V and 2UV were, respectively, higher and lower compared to baseline examination, suggesting more active sympathetic and less active vagal control. Importantly, such changes of these SymDyn indexes (calculated using the Max–min and σ methods) were observed for all athletes. Successful detecting and recognizing individual athletes’ cardiac autonomic activity could be helpful in achieving better sport results as it was shown that, for example, sympathetic predominance prior and/or during competition could be advantageous for performance in extreme sports professional athletes (Matsumura et al., 2021). Analysis of inter-individual responses of modifiable biomarkers to a specified stimulus may help in the identification of athletes that will benefit from practical techniques aimed at avoiding pre-performance stress and improve performance in sports.
To monitor athletes in the field using HRV indices Buchheit discouraged practitioners using spectral indices and recommended more appropriate linear time-domain parameters (RMSSD or SD1 reflected parasympathetic modulation) due to their calculation simplicity (using, e.g., an Excel spreadsheet) based on short-term or ultra-short-term recordings and very low sensitivity to breathing pattern (Buchheit, 2014). SymDyn indices may be a good alternative for frequency-domain parameters: 0V and 2UV have been demonstrated to be correlated with the sympathetic and parasympathetic autonomic modulation to the heart (Porta et al., 2001, 2007a,b; Guzzetti et al., 2005; Baumert et al., 2009; Tobaldini et al., 2009; Cysarz et al., 2013, 2018; Silva et al., 2017), can be easily calculated based on very short-term recordings (Storniolo et al., 2021) using recently presented, available free of charge, software tool—PyBioS (Silva et al., 2020).
Small sample size of only male participants, cross-sectional study design, data obtained during supine ECG recordings performed in controlled laboratory settings and lack of additional collection of RR intervals using HR monitor to provide more practical aspect of sports field should be recognized as limitations of this pilot study.
Conclusion
Data from the current pilot investigation indicate that short-term SymDyn indices may be used as reliable non-respiratory-associated parameters in laboratory settings to detect ANS activity modulations in elite modern male pentathlon athletes. These findings provide a potential solution for addressing the confounding influence of respiration frequency on HRV-derived inferences of cardiac autonomic function. For this reason, SymDyn may prove to be preferable for field-based monitoring where measurements are unsupervised.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by University Ethical Committee (SKE 01-01/2017, 7 March 2017, Józef Piłsudski University of Physical Education in Warsaw, Warsaw, Poland). The patients/participants provided their written informed consent to participate in this study.
Author Contributions
JG and BH: conceptualization and methodology. JG, MR, MM, BH, and RB: formal analysis and investigation. JG, MR, and MM: writing—original draft preparation. JG, MR, MM, AF, RB, and BW: writing—review and editing. All authors contributed to the article and approved the submitted version.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We gratefully acknowledge Luiz Eduardo Virgilio da Silva for his contributions to the analysis of the data and for support in developing this manuscript. We thank the athletes for participation in the study. We thank reviewers for their valuable suggestions.
References
Abu-Arafeh, A., Jordan, H., and Drummond, G. (2016). Reporting of method comparison studies: a review of advice, an assessment of current practice, and specific suggestions for future reports. Br. J. Anaesth. 117, 569–575. doi: 10.1093/bja/aew320
Acharya, U. R., Joseph, K. P., Kannathal, N., Lim, C. M., and Suri, J. S. (2006). Heart rate variability: a review. Med. Biol. Eng. Comput. 44, 1031–1051. doi: 10.1007/s11517-006-0119-0
Al Haddad, H., Laursen, P. B., Chollet, D., Ahmaidi, S., and Buchheit, M. (2011). Reliability of resting and postexercise heart rate measures. Int. J. Sports Med. 32, 598–605. doi: 10.1055/s-0031-1275356
Araújo, C. G., and Scharhag, J. (2016). Athlete: a working definition for medical and health sciences research. Scand. J. Med. Sci. Sports 26, 4–7. doi: 10.1111/sms.12632
Atkinson, G., and Nevill, A. M. (1998). Statistical methods for assessing measurement error (reliability) in variables relevant to sports medicine. Sports Med. 26, 217–238. doi: 10.2165/00007256-199826040-00002
Bari, V., Marchi, A., De Maria, B., Rossato, G., Nollo, G., Faes, L., et al. (2016). Nonlinear effects of respiration on the crosstalk between cardiovascular and cerebrovascular control systems. Philos. Trans. A Math. Phys. Eng. Sci. 374:20150179. doi: 10.1098/rsta.2015.0179
Baumert, M., Brechtel, L., Lock, J., and Voss, A. (2006). Changes in heart rate variability of athletes during a training camp. Biomed. Tech. 51, 201–204. doi: 10.1515/BMT.2006.037
Baumert, M., Lambert, G. W., Dawood, T., Lambert, E. A., Esler, M. D., McGrane, M., et al. (2009). Short-term heart rate variability and cardiac norepinephrine spillover in patients with depression and panic disorder. Am. J. Physiol. Heart Circ. Physiol. 297, H674–H679. doi: 10.1152/ajpheart.00236.2009
Bishop, S. A., Dech, R. T., Guzik, P., and Neary, J. P. (2018). Heart rate variability and implication for sport concussion. Clin. Physiol. Funct. Imaging 38, 733–742. doi: 10.1111/cpf.12487
Bland, J. M., and Altman, D. G. (1996). Statistics notes: measurement error proportional to the mean. Br. Med. J. 313:106. doi: 10.1136/bmj.313.7049.106
Bourdon, P. C., Cardinale, M., Murray, A., Gastin, P., Kellmann, M., Varley, M. C., et al. (2017). Monitoring athlete training loads: consensus statement. Int. J. Sports Physiol. Perform. 12, S2-161–S2-170. doi: 10.1123/IJSPP.2017-0208
Brown, T. E., Beightol, L. A., Koh, J., and Eckberg, D. L. (1993). Important influence of respiration on human R-R interval power spectra is largely ignored. J. Appl. Physiol. 75, 2310–2317. doi: 10.1152/jappl.1993.75.5.2310
Buchheit, M. (2014). Monitoring training status with HR measures: do all roads lead to Rome? Front. Physiol. 5:73. doi: 10.3389/fphys.2014.00073
Cohen, J. (1988). Statistical Power Analysis for the Behavioral Sciences, 2nd ed. London, UK, Routledge, pages: 179–206.
Cysarz, D., Edelhauser, F., Javorka, M., Montano, N., and Porta, A. (2015). A percentile-based coarse graining approach is helpful in symbolizing heart rate variability during graded head-up tilt. Annu. Int. Conf. IEEE Eng. Med. Biol. Soc. 2015, 286–289. doi: 10.1109/EMBC.2015.7318356
Cysarz, D., Edelhäuser, F., Javorka, M., Montano, N., and Porta, A. (2018). On the relevance of symbolizing heart rate variability by means of a percentile-based coarse graining approach. Physiol. Meas. 39:105010. doi: 10.1088/1361-6579/aae302
Cysarz, D., Porta, A., Montano, N., Leeuwen, P. V., Kurths, J., and Wessel, N. (2013). Quantifying heart rate dynamics using different approaches of symbolic dynamics. Eur. Phys. J. Spec. Top. 222, 487–500. doi: 10.1140/epjst/e2013-01854-7
De Souza, L. A., Ferreira, J. B., Schein, A. S. O., Dartora, D. R., Casali, A. G., Scassola, C. M. C., et al. (2018). Optimization of vagal stimulation protocol based on spontaneous breathing rate. Front. Physiol. 9:1341. doi: 10.3389/fphys.2018.01341
Deus, L. A., Sousa, C. V., Rosa, T. S., Filho, J. M. S., Santos, P. A., Barbosa, L. D., et al. (2019). Heart rate variability in middle-aged sprint and endurance athletes. Physiol. Behav. 205, 39–43. doi: 10.1016/j.physbeh.2018.10.018
Eckberg, D. L. (2003). The human respiratory gate. J. Physiol. 548, 339–352. doi: 10.1113/jphysiol.2002.037192
Eckberg, D. L., and Eckberg, M. J. (1982). Human sinus node responses to repetitive, ramped carotid baroreceptor stimuli. Am. J. Phys. 242, H638–H644. doi: 10.1152/ajpheart.1982.242.4.H638
Eckberg, D. L., Nerhed, C., and Wallin, B. G. (1985). Respiratory modulation of muscle sympathetic and vagal cardiac outflow in man. J. Physiol. 365, 181–196. doi: 10.1113/jphysiol.1985.sp015766
Esco, M. R., Williford, H. N., Flatt, A. A., Freeborn, T. J., and Nakamura, F. Y. (2018). Ultra-shortened time-domain HRV parameters at rest and following exercise in athletes: an alternative to frequency computation of sympathovagal balance. Eur. J. Appl. Physiol. 118, 175–184. doi: 10.1007/s00421-017-3759-x
Faust, O., Acharya, U. R., Krishnan, S. M., and Min, L. C. (2004). Analysis of cardiac signals using spatial filling index and time-frequency domain. Biomed. Eng. Online 3:30. doi: 10.1186/1475-925X-3-30
Flatt, A. A., and Esco, M. R. (2016). Heart rate variability stabilization in athletes: towards more convenient data acquisition. Clin. Physiol. Funct. Imaging 36, 331–336. doi: 10.1111/cpf.12233
Gabbett, T. J., Nassis, G. P., Oetter, E., Pretorius, J., Johnston, N., Medina, D., et al. (2017). The athlete monitoring cycle. A practical guide to interpreting and applying training monitoring data. Br. J. Sports Med. 51, 1451–1452. doi: 10.1136/bjsports-2016-097298
Gąsior, J. S., Hoffmann, B., Silva, L. E. V., Małek, Ł., Flatt, A. A., Baranowski, R., et al. (2020). Changes in short-term and ultra-short term heart rate, respiratory rate, and time-domain heart rate variability parameters during sympathetic nervous system activity stimulation in elite modern pentathlonists: a pilot study. Diagnostics 10:1104. doi: 10.3390/diagnostics10121104
Gerke, O. (2020). Reporting standards for a Bland-Altman agreement analysis: a review of methodological reviews. Diagnostics 10:334. doi: 10.3390/diagnostics10050334
Goldberger, A. L. (1991). Is the normal heartbeat chaotic or homeostatic? News Physiol. Sci. 6, 87–91. doi: 10.1152/physiologyonline.1991.6.2.87
Gronwald, T., and Hoos, O. (2020). Correlation properties of heart rate variability during endurance exercise: a systematic review. Ann. Noninvasive Electrocardiol. 25:e12697. doi: 10.1111/anec.12697
Gronwald, T., Rogers, B., and Hoos, O. (2020). Fractal correlation properties of heart rate variability: a new biomarker for intensity distribution in endurance exercise and training prescription? Front. Physiol. 11:550572. doi: 10.3389/fphys.2020.550572
Guzzetti, S., Borroni, E., Garbelli, P. E., Ceriani, E., Della Bella, P., Montano, N., et al. (2005). Symbolic dynamics of heart rate variability: a probe to investigate cardiac autonomic modulation. Circulation 112, 465–470. doi: 10.1161/CIRCULATIONAHA.104.518449
Heidari, J., Beckmann, J., Bertollo, M., Brink, M., Kallus, W., Robazza, C., et al. (2018). Multidimensional monitoring of recovery status and implications for performance. Int. J. Sports Physiol. Perform. 1–24. doi: 10.1123/ijspp.2017-0669
Henriques, T., Ribeiro, M., Teixeira, A., Castro, L., Antunes, L., and Costa-Santos, C. (2020). Nonlinear methods Most applied to heart-rate time series: a review. Entropy 22:309. doi: 10.3390/e22030309
Hirsch, J. A., and Bishop, B. (1981). Respiratory sinus arrhythmia in humans: how breathing pattern modulates heart rate. Am. J. Physiol. Heart Circ. Physiol. 241, H620–H629. doi: 10.1152/ajpheart.1981.241.4.H620
Hoffmann, B., Flatt, A. A., Silva, L. E. V., Młyńczak, M., Baranowski, R., Dziedzic, E., et al. (2020). A pilot study of the reliability and agreement of heart rate, respiratory rate and short-term heart rate variability in elite modern pentathlon athletes. Diagnostics 10:833. doi: 10.3390/diagnostics10100833
Holmes, C. J., Fedewa, M. V., Dobbs, W. C., Liu, Y., Flatt, A. A., Nakamura, F. B., et al. (2020). The effects of different body positions on the accuracy of ultra-short-term heart rate variability indexes. J. High Technol. Manag. Res. 31:100375. doi: 10.1016/j.hitech.2020.100375
Hopkins, W. G. (2000a). A New View of Statistics. Available at: http://sportsci.org/resource/stats (Accessed August 3, 2020).
Hopkins, W. G. (2000b). Measures of reliability in sports medicine and science. Sports Med. 30, 1–15. doi: 10.2165/00007256-200030010-00001
Hopkins, W. G., Marshall, S. W., Batterham, A. M., and Hanin, J. (2009). Progressive statistics for studies in sports medicine and exercise science. Med. Sci. Sports Exerc. 41, 3–12. doi: 10.1249/MSS.0b013e31818cb278
Hoshi, R. A., Pastre, C. M., Vanderlei, L. C., and Godoy, M. F. (2013). Poincaré plot indexes of heart rate variability: relationships with other nonlinear variables. Auton. Neurosci. 177, 271–274. doi: 10.1016/j.autneu.2013.05.004
Huikuri, H. V., Perkiomaki, J. S., Maestri, R., and Pinna, G. D. (2009). Clinical impact of evaluation of cardiovascular control by novel methods of heart rate dynamics. Philos. Trans. R. Soc-Math. Phys. Eng. Sci. 367, 1223–1238. doi: 10.1098/rsta.2008.0294
Kellmann, M., Bertollo, M., Bosquet, L., Brink, M., Coutts, A. J., Duffeld, R., et al. (2018). Recovery and performance in sport: consensus statement. Int. J. Sports Physiol. Perform. 13, 240–245. doi: 10.1123/ijspp.2017-0759
Khurana, R. K., and Setty, A. (1996). The value of the isometric hand-grip test-studies in various autonomic disorders. Clin. Auton. Res. 6, 211–218. doi: 10.1007/BF02291136
Kurths, J., Voss, A., Saparin, P., Witt, A., Kleiner, H. J., and Wessel, N. (1995). Quantitative analysis of heart rate variability. Chaos 5, 88–94. doi: 10.1063/1.166090
Laborde, S., Mosley, E., and Thayer, J. F. (2017). Heart rate variability and cardiac vagal tone in psychophysiological research - recommendations for experiment planning, data analysis, and data reporting. Front. Psychol. 8:213. doi: 10.3389/fpsyg.2017.00213
Lachin, J. M. (2004). The role of measurement reliability in clinical trials. Clin. Trials 1, 553–566. doi: 10.1191/1740774504cn057oa
Lucini, D., Marchetti, I., Spataro, A., Malacarne, M., Benzi, M., Tamorri, S., et al. (2017). Heart rate variability to monitor performance in elite athletes: criticalities and avoidable pitfalls. Int. J. Cardiol. 240, 307–312. doi: 10.1016/j.ijcard.2017.05.001
Maestri, R., Pinna, G. D., Porta, A., Balocchi, R., Sassi, R., Signorini, M. G., et al. (2007). Assessing nonlinear properties of heart rate variability from short-term recordings: are these measurements reliable? Physiol. Meas. 28, 1067–1077. doi: 10.1088/0967-3334/28/9/008
Magagnin, V., Bassani, T., and Bari, V. (2011). Non-stationarities significantly distort short-term spectral, symbolic and entropy heart rate variability indices. Physiol. Meas. 32, 1775–1786. doi: 10.1088/0967-3334/32/11/S05
Mann, T. N., Lamberts, R. P., and Lambert, M. I. (2014). High responders and low responders: factors associated with individual variation in response to standardized training. Sports Med. 44, 1113–1124. doi: 10.1007/s40279-014-0197-3
Matheson, G. J. (2019). We need to talk about reliability: making better use of test-retest studies for study design and interpretation. Peer J. 7:e6918. doi: 10.7717/peerj.6918
Matsumura, S., Watanabe, K., Saijo, N., Ooishi, Y., Kimura, T., and Kashino, M. (2021). Positive relationship between precompetitive sympathetic predominance and competitive performance in elite extreme sports athletes. Front. Sports Act. Living 3:712439. doi: 10.3389/fspor.2021.712439
Meeusen, R., Duclos, M., Foster, C., Fry, A., Gleeson, M., Nieman, D., et al. (2013). Prevention, diagnosis, and treatment of the overtraining syndrome: joint consensus statement of the European College of Sport Science and the American College of Sports Medicine. Med. Sci. Sports Exerc. 45, 186–205. doi: 10.1249/MSS.0b013e318279a10a
Michael, S., Graham, K. S., and Davis Oam, G. M. (2017). Cardiac autonomic responses during exercise and post-exercise recovery using heart rate variability and systolic time intervals—a review. Front. Physiol. 8:301. doi: 10.3389/fphys.2017.00301
Młyńczak, M., and Krysztofiak, H. (2019). Cardiorespiratory temporal causal links and the differences by sport or lack thereof. Front. Physiol. 10:45. doi: 10.3389/fphys.2019.00045
Młyńczak, M., Niewiadomski, W., Żyliński, M., and Cybulski, G. P. (2014). “Ambulatory impedance pneumography device for quantitative monitoring of volumetric parameters in respiratory and cardiac applications.” in Proceedings of the Computing in Cardiology Conference (CinC), Cambridge, MA, USA, September 7–10, 2014.
Młyńczak, M., Niewiadomski, W., Żyliński, M., and Cybulski, G. (2017). “Ambulatory devices measuring cardiorespiratory activity with motion.” in Proceedings of the 10th International Joint Conference on Biomedical Engineering Systems and Technologies, Porto, Portugal. February 21–23, 2017; 91–97.
Muñoz-López, A., and Naranjo-Orellana, J. (2020). Individual versus team heart rate variability responsiveness analyses in a national soccer team during training camps. Sci. Rep. 10:11726. doi: 10.1038/s41598-020-68698-5
Narkiewicz, K., van de Borne, P., Montano, N., Hering, D., Kara, T., and Somers, V. K. (2006). Sympathetic neural outflow and chemoreflex sensitivity are related to spontaneous breathing rate in normal men. Hypertension 47, 51–55. doi: 10.1161/01.HYP.0000197613.47649.02
Peltola, M. A. (2012). Role of editing of R-R intervals in the analysis of heart rate variability. Front. Physiol. 3:148. doi: 10.3389/fphys.2012.00148
Pereira, L. A., Flatt, A. A., Ramirez-Campillo, R., Loturco, I., and Nakamura, F. Y. (2016). Assessing shortened field-based heart-rate-variability-data acquisition in team-sport athletes. Int. J. Sports Physiol. Perform. 11, 154–158. doi: 10.1123/ijspp.2015-0038
Perrone, M. A., Volterrani, M., Manzi, V., Barchiesi, F., and Iellamo, F. (2021). Heart rate variability modifications in response to different types of exercise training in athletes. J. Sports Med. Phys. Fitness 61, 1411–1415. doi: 10.23736/S0022-4707.21.12480-6
Plews, D. J., Laursen, P. B., Stanley, J., Kilding, A. E., and Buchheit, M. (2013). Training adaptation and heart rate variability in elite endurance athletes: opening the door to effective monitoring. Sports Med. 43, 773–781. doi: 10.1007/s40279-013-0071-8
Porta, A., Faes, L., Masé, M., D’Addio, G., Pinna, G. D., Maestri, R., et al. (2007a). An integrated approach based on uniform quantization for the evaluation of complexity of short-term heart period variability: application to 24h Holter recordings in healthy and heart failure humans. Chaos 17:015117. doi: 10.1063/1.2404630
Porta, A., Guzzetti, S., Montano, N., Furlan, R., Pagani, M., Malliani, A., et al. (2001). Entropy, entropy rate, and pattern classification as tools to typify complexity in short heart period variability series. IEEE Trans. Biomed. Eng. 48, 1282–1291. doi: 10.1109/10.959324
Porta, A., Marchi, A., Bari, V., Heusser, K., Tank, J., Jordan, J., et al. (2015). Conditional symbolic analysis detects nonlinear influences of respiration on cardiovascular control in humans. Philos. Trans. A Math. Phys. Eng. Sci. 373:20140096. doi: 10.1098/rsta.2014.0096
Porta, A., Tobaldini, E., Guzzetti, S., Furlan, R., Montano, N., and Gnecchi-Ruscone, T. (2007b). Assessment of cardiac autonomic modulation during graded head-up tilt by symbolic analysis of heart rate variability. Am. J. Physiol. Heart Circ. Physiol. 293, H702–H708. doi: 10.1152/ajpheart.00006.2007
Quintana, D. S., Alvares, G. A., and Heathers, J. A. J. (2016). Guidelines for reporting articles on psychiatry and heart rate variability (GRAPH): recommendations to advance research communication. Transl. Psychiatry 6:e803. doi: 10.1038/tp.2016.73
Quintana, D. S., and Heathers, J. A. J. (2014). Considerations in the assessment of heart rate variability in biobehavioral research. Front. Psychol. 5:805. doi: 10.3389/fpsyg.2014.00805
Rogers, B., Giles, D., Draper, N., Hoos, O., and Gronwald, T. (2021a). A new detection method defining the aerobic threshold for endurance exercise and training prescription based on fractal correlation properties of heart rate variability. Front. Physiol. 11:596567. doi: 10.3389/fphys.2020.596567
Rogers, B., Mourot, L., Doucende, G., and Gronwald, T. (2021b). Fractal correlation properties of heart rate variability as a biomarker of endurance exercise fatigue in ultramarathon runners. Phys. Rep. 9:e14956. doi: 10.14814/phy2.14956
Sassi, R., Cerutti, S., Lombardi, F., Malik, M., Huikuri, H. V., Peng, C. K., et al. (2015). Advances in heart rate variability signal analysis: joint position statement by the e-cardiology ESC working group and the European heart rhythm association co-endorsed by the Asia Pacific Heart Rhythm Society. Europace 17, 1341–1353. doi: 10.1093/europace/euv015
Saul, J. P., and Valenza, G. (2021). Heart rate variability and the dawn of complex physiological signal analysis: methodological and clinical perspectives. Philos. Trans. A Math. Phys. Eng. Sci. 379:20200255. doi: 10.1098/rsta.2020.0255
Schiecke, K., Pester, B., Schumann, A., and Baer, K. J. (2018). Nonlinear interaction analysis of cardiovascular-respiratory data by means of convergent cross mapping. Annu. Int. Conf. IEEE Eng. Med. Biol. Soc. 2018, 267–270. doi: 10.1109/EMBC.2018.8512213
Schmitt, L., Willis, S. J., Fardel, A., Coulmy, N., and Millet, G. P. (2018). Live high-train low guided by daily heart rate variability in elite nordic-skiers. Eur. J. Appl. Physiol. 118, 419–428. doi: 10.1007/s00421-017-3784-9
Schneider, C., Hanakam, F., Wiewelhove, T., Döweling, A., Kellmann, M., Meyer, T., et al. (2018). Heart rate monitoring in team sports: a conceptual framework for contextualizing heart rate measures for training and recovery prescription. Front. Physiol. 9:639. doi: 10.3389/fphys.2018.00639
Shoukri, M. M., Colak, D., Kaya, N., and Donner, A. (2008). Comparison of two dependent within subject coefficients of variation to evaluate the reproducibility of measurement devices. BMC Med. Res. Methodol. 8:24. doi: 10.1186/1471-2288-8-24
Shrout, P. E., and Fleiss, J. L. (1979). Intraclass correlations: uses in assessing rater reliability. Psychol. Bull. 86, 420–428. doi: 10.1037//0033-2909.86.2.420
Silva, L. E. V., Fazan, R. Jr., and Marin-Neto, J. A. (2020). PyBioS: a freeware computer software for analysis of cardiovascular signals. Comput. Methods Prog. Biomed. 197:105718. doi: 10.1016/j.cmpb.2020.105718
Silva, L. E. V., Rezende Geraldini, V., Potratz de Oliveira, B., Aguiar Silva, C. A., Porta, A., and Fazan, R. (2017). Comparison between spectral analysis and symbolic dynamics for heart rate variability analysis in the rat. Sci. Rep. 7:8428. doi: 10.1038/s41598-017-08888-w
Storniolo, J. L., Cairo, B., Porta, A., and Cavallari, P. (2021). Symbolic analysis of the heart rate variability during the plateau phase following maximal sprint exercise. Front. Physiol. 12:632883. doi: 10.3389/fphys.2021.632883
Tobaldini, E., Porta, A., Wei, S.-G., Zhang, Z.-H., Francis, J., Casali, K. R., et al. (2009). Symbolic analysis detects alterations of cardiac autonomic modulation in congestive heart failure rats. Auton. Neurosci. 150, 21–26. doi: 10.1016/j.autneu.2009.03.009
Uhlig, S., Meylan, A., and Rudolph, U. (2020). Reliability of short-term measurements of heart rate variability: findings from a longitudinal study. Biol. Psychol. 154:107905. doi: 10.1016/j.biopsycho.2020.107905
Vescovi, J. D. (2019). Intra-individual variation of HRV during orthostatic challenge in elite male field hockey players. J. Med. Syst. 43:328. doi: 10.1007/s10916-019-1478-z
Voss, A., Schulz, S., Schroeder, R., Baumert, M., and Caminal, P. (2009). Methods derived from nonlinear dynamics for analysing heart rate variability. Phil. Trans. R. Soc. A 367, 277–296. doi: 10.1098/rsta.2008.0232
Keywords: symbolic dynamics, heart rate variability, non-linear, athletes, modern pentathlon athletes
Citation: Gąsior JS, Rosoł M, Młyńczak M, Flatt AA, Hoffmann B, Baranowski R and Werner B (2022) Reliability of Symbolic Analysis of Heart Rate Variability and Its Changes During Sympathetic Stimulation in Elite Modern Pentathlon Athletes: A Pilot Study. Front. Physiol. 13:829887. doi: 10.3389/fphys.2022.829887
Edited by:
Julien Louis, Liverpool John Moores University, United KingdomReviewed by:
Robin Pla, French Institute of Sport (INSEP), FranceMathias Baumert, University of Adelaide, Australia
Copyright © 2022 Gąsior, Rosoł, Młyńczak, Flatt, Hoffmann, Baranowski and Werner. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jakub S. Gąsior, jgasior@wum.edu.pl