Applicability of Dmax Method on Heart Rate Variability to Estimate the Lactate Thresholds in Male Runners

Nascimento, Eduardo Marcel Fernandes; Antunes, Diego; do Nascimento Salvador, Paulo Cesar; Borszcz, Fernando Klitzke; de Lucas, Ricardo Dantas

doi:https://doi.org/10.1155/2019/2075371

Journal of Sports Medicine

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion Disclosure Abbreviations Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2019 | Article ID 2075371 | https://doi.org/10.1155/2019/2075371

Applicability of Dmax Method on Heart Rate Variability to Estimate the Lactate Thresholds in Male Runners

Eduardo Marcel Fernandes Nascimento,¹Diego Antunes,¹Paulo Cesar do Nascimento Salvador,¹Fernando Klitzke Borszcz,¹and Ricardo Dantas de Lucas¹

Academic Editor: Mark Willems

Received01 Apr 2019

Revised20 Jun 2019

Accepted03 Jul 2019

Published19 Sept 2019

Abstract

Introduction. The purpose of this study was to evaluate the application of the Dmax method on heart rate variability (HRV) to estimate the lactate thresholds (LT), during a maximal incremental running test (MIRT). Methods. Nineteen male runners performed two MIRTs, with the initial speed at 8 km·h⁻¹ and increments of 1 km·h⁻¹ every 3 minutes, until exhaustion. Measures of HRV and blood lactate concentrations were obtained, and lactate (LT₁ and LT₂) and HRV (HRVT_DMAX1 and HRVT_DMAX2) thresholds were identified. ANOVA with Scheffe’s post hoc test, effect sizes (d), the bias ± 95% limits of agreement (LoA), standard error of the estimate (SEE), Pearson’s (r), and intraclass correlation coefficient (ICC) were calculated to assess validity. Results. No significant differences were observed between HRVT_DMAX1 and LT₁ when expressed for speed (12.1 ± 1.4 km·h⁻¹ and 11.2 ± 2.1 km·h⁻¹; ; d = 0.45; r = 0.46; bias ± LoA = 0.8 ± 3.7 km·h⁻¹; SEE = 1.2 km·h⁻¹ (95% CI, 0.9–1.9)). Significant differences were observed between HRVT_DMAX2 and LT₂ when expressed for speed (12.0 ± 1.2 km·h⁻¹ and 14.1 ± 2.5 km·h⁻¹; ; d = 1.21; r = 0.48; bias ± LoA = −1.0 ± 1.8 km·h⁻¹; SEE = 1.1 km·h⁻¹ (95% CI, 0.8–1.6)), respectively. Reproducibility values were found for the LT₁ (ICC = 0.90; bias ± LoA = −0.7 ± 2.0 km·h⁻¹), LT₂ (ICC = 0.97; bias ± LoA = −0.1 ± 1.1 km·h⁻¹), HRVT_DMAX1 (ICC = 0.48; bias ± LoA = −0.2 ± 3.4 km·h⁻¹), and HRVT_DMAX2 (ICC = 0.30; bias ± LoA = 0.3 ± 3.5 km·h⁻¹). Conclusions. The Dmax method applied over a HRV dataset allowed the identification of LT₁ that is close to aerobic threshold, during a MIRT.

1. Introduction

The autonomic cardiac drive can be investigated by the heart rate variability (HRV), which is characterized as a variation quantified in milliseconds between RR intervals [1]. A predominance of parasympathetic nervous system (PNS) activity is observed at rest and low effort intensities. In approximately 50–60% of the maximum oxygen uptake (VO_2MAX), a significant vagal withdrawal occurs [2]. The aerobic threshold (AeT) has been related to that intensity [3–5], i.e., the exercise intensity which lactate concentrations [La] initiate to increase beyond resting values and are frequently called “lactate threshold” (LT) [6, 7]. On the other hand, in intensities above the AeT, there is a gradual and constant increase in activation in the sympathetic nervous system (SNS), and a marked increase in the physiological responses related to the anaerobic threshold (AnT) can be observed [3, 4]. That intensity is corresponding to maximal lactate steady state (MLSS), i.e., the highest constant exercise intensity output that can be maintained over time without continual [La] accumulation [3, 6].

Since AeT and AnT are good indexes of ideal training intensity and determinants of endurance performance [3, 4, 6], there is an obvious and growing interest in proposing different methods to estimate that intensities, mainly in relation to MSSL, considered gold standard endurance performance marker [3, 6, 8]. Among the different methods, the HRV thresholds (HRVTs) have been highlighted [5, 9–15]. The HRVT method is accessible and of low cost to use compared to traditional methods, such as [La] and gas exchange analysis, as well as a real noninvasive alternative to routine applications. The main methods for identifying the HRVTs are frequency domain [10, 12, 13], time domain [5, 9, 15], and nonlinear [10, 14] analysis. Success in identifying HRVTs and consequent estimation of AeT and AnT have been confirmed in different situations, such as in running [10–12], cycling [5, 9, 13–15], and ski mountaineering [16, 17], and also in high-level swimmers [18], trained boys [10], obese adolescents [19], and individuals with type 2 diabetes [20].

The Poincaré plot is a nonlinear HRV analysis method that uses time domain markers [21] and is an important research area, since it allows its use in nonstationary data, a characteristic inherent to HRV, especially during the increase of effort intensity [22]. The Poincaréplot analysis provides the calculation of the standard deviation of instantaneous (SD1) and continuous long-term RR intervals (SD2) [1]. The SD1 marker has been shown to correlate strongly with vagal tone (PNS), and previous studies have pointed to an abrupt point of change in their behavior in intensities related to AeT [2, 5, 9, 11]. SD2 marker has been shown to correlate with the PNS and SNS, and this variable shows a nonlinear pattern in intensities close to heavy and severe domains [2, 11]. The applicability and efficacy of the Poincaré plot in the estimation of the AeT and AnT have been confirmed in previous studies in running [11] and cycling [9, 15, 23].

However, in addition to the heart rate (HR) which presents theoretical support for a nonlinear pattern, especially in intensity close to AnT [24], some aspects need to be better elucidated when using HRV markers for the estimation of AeT and AnT. Firstly, it would be the validity of the method since the majority of studies used visual analysis for HRVT identification [9–12, 16], which is influenced by the subjective aspect and experience of the evaluator. In order to remedy this limitation, Cheng et al. [25] proposed the Dmax method to identify the lactate and ventilatory thresholds. Previous studies demonstrated greater reliability of the Dmax method than visual analysis or the use of fixed [La] [26]. The Dmax method presents an important advantage which a breakpoint can always be detected [25], although a maximal test is needed. Only one study of our knowledge used the Dmax method to identify HRVTs [23]. Their results surprisingly on the contrary, as reported by the authors, pointed to the visual analysis as better indicators of reliability in the SD1 and RMSSD (square root of the mean squared differences of successive RR intervals) markers than the Dmax method. This way, it is doubted if the Dmax method is better or not than the visual analysis for the identification of AeT or AnT, when using HRV dataset. Nevertheless, the results of the aforementioned study [23] were not compared with traditional methods to estimate AeT and AnT, as [La] or gaseous exchanges; therefore, greater inferences are limited. Secondly, no study to date has analyzed the possibility to identify the HRVTs by the Dmax method in the different situations and conditions compared to lactate and ventilatory thresholds. Finally, the reproducibility of the method must be verified in relation to different situations and conditions, since it has only been tested on the cycle [15, 23].

The identification of HRVTs and consequently the estimation of MLSS can be a framework very important to control and monitor training workloads, as well as to assess the improvement in performance during an endurance training program [21]. The applicability of HRV thresholds in a single-day test perhaps can be very attractive for research studies of sports science, trainers, and practitioners users. Therefore, the aim of the study is to investigate the application of the Dmax method originally proposed by Cheng et al. [25], on HRV dataset to estimate the AeT and AnT, during a maximal incremental running test (MIRT) in male runners. Firstly, the hypothesis is that the HRVT identified by SD1 marker (HRVT_DMAX1) could be used to estimate the AeT, since this marker has been shown to correlate with the PNS [2, 9, 11, 15]. Secondly, the hypothesis is that the HRVT identified by SD2 marker (HRVT_DMAX2) could be used to estimate AnT, since this marker has been shown to correlate with a significant PNS and mainly with the SNS [11]. The reproducibility of the HRVT_DMAX1 and HRVT_DMAX2, as well as AeT and AnT, will be verified.

2. Materials and Methods

2.1. Participants

Nineteen male recreational long-distance runners (30.4 ± 4.0 years; body mass of 74.3 ± 8.4 kg; height of 176 ± 6.3 cm; body fat of 13.8 ± 4.5%) volunteered to participate in this study. All participants were healthy, without cardiovascular or orthopedic problems, nonsmokers, and not taking any medication. The study protocol complied with the Declaration of Helsinki for human experimentation [27] and was approved by the Institutional Ethics Committee of the University of São Paulo.

2.2. Study Design

All the participants performed two MIRTs interspersed by a washout period of 3–7 days. The tests were performed at the same time of the day and in standard laboratory conditions (humidity of ≈50% and temperature of ≈22°C). Participants were instructed to avoid intense exercises, alcohol, and caffeine beverages 24 hours before each test and to consume a light meal 3 hours before the tests.

2.3. Maximal Incremental Running Test Protocol

Before MIRT, participants used a cardio belt for beat-to-beat heart rate (HR) measures (S810 Polar®, Kempele, Finland), during rest for 20 min (10 min supine + 10 min sitting) for baseline measures of HR, HRV, and [La]. The [La] was obtained from a 25 μL blood, drawn from the tip of the forefinger, and blood samples were then stored in Eppendorf tubes containing 50 μL of 1% NaF in a −30°C environment, according to the recommendations of the manufacturer (YSI 1500 Sport, Yellow Springs, OH, USA). Later, the samples were analyzed using enzyme electrode technology (YSI 1500 Sport, Yellow Springs, OH, USA). Then, the participants were directed on the treadmill (CEFISE TK35, Nova Odessa, Brazil) and warmed up for 3 min at 5 km·h⁻¹ and 1% gradient. The test started at 8 km·h⁻¹, with 1 km·h⁻¹ increases every 3 min, until exhaustion, being a protocol adapted by Heck et al. [28]. The HR dataset was recorded continuously throughout the tests. Blood samples of 25 μL were collected during the last 30 s of every stage, while the participant was running.

2.4. Determination of Aerobic and Anaerobic Thresholds

2.4.1. Heart Rate Variability Thresholds

The Dmax method was used to analyze the behavior of SD1 and SD2 markers, to identify the HRVT_DMAX1 and HRVT_DMAX2, respectively (Figure 1). The Dmax method was determined according to Cheng et al. [25], thereby providing individualized lactate threshold values according to the following equation:where Y represents the predicted values of SD1 or SD2 at a given workload (km·h⁻¹); a₃, a₂, a₁, and a₀ are the intercepts; and x is the speed. Briefly, the Dmax method reflects the longest perpendicular distance between SD1 and SD2 values predicted by a third-order polynomial function over actual (SD1 and SD2) values and values derived from a linear regression calculated with the first and last values of each curve, respectively. The SD1 marker was used because it has been shown to correlate strongly with vagal tone (PNS), and previous studies have pointed to an abrupt point of change in their behavior in intensities related to AeT [2, 5, 9, 11, 15]. SD2 marker has been shown to correlate with the PNS and SNS, and this variable shows a nonlinear pattern in intensities close to heavy and severe domains, being these related to AnT [2, 11]. Thereafter, raw RR intervals were recorded during the last 60 s of each stage of the exercise, and then the Dmax method was applied on the measured values.

(a)

(b)

(c)

2.4.2. Lactate Thresholds

The first lactate threshold (LT₁) (i.e., AeT) was determined as the lowest value of the ratio [La]/speed [29]. After, the second lactate threshold (LT₂) (i.e. AnT) was determined as the running speed at 1.5 mmol·L⁻¹ above LT₁ (Figure 1) [29]. The LT₁ and LT₂ derived from the Dickhuth et al.’s [29] methods were used as criterium measures, because LT₁ has a high correlation with MLSS [7], and LT₂ (+1.5 mmol·L⁻¹) was used because it showed a high concordance with MLSS in runners during the MIRT with stages of 3 min [30].

All the thresholds, HRVT_DMAX1, HRVT_DMAX2, LT₁, and LT₂, were expressed as absolute and relative values for speed (km·h⁻¹), [La] (mmol·L⁻¹), milliseconds beat-to-beat RR intervals (ms), and HR (bpm).

2.5. Statistical Procedures

Values were expressed as mean and standard deviation (±SD). After ensuring Gaussian data distribution (normality and homoscedasticity), a spreadsheet was used for the analysis of concurrent validity [31] and statistical standards were followed [32]. Cohen’s [33] (d) effect sizes and ANOVA with Scheffe’s post hoc test were used to compare the magnitude of the differences between the thresholds LT₁ and LT₂ with HRVT_DMAX1 and HRVT_DMAX2, respectively. Additionally, the standard error of the estimate (SEE), the bias ± 95% of limits of agreement [LoA] of the Bland and Altman analysis [34], and the Pearson product-moment correlation were used to evaluate the association between the different methods for identifying thresholds. For measures, reliability determination, the intraclass correlation coefficient (ICC), and the typical error of measurement (TEM) were performed using a Hopkins spreadsheet [31]. The d values were interpreted using the following scale: <0.20 (trivial), 0.2–0.6 (small), 0.6–1.2 (moderate), 1.2–2.0 (large), 2.0–4.0 (very large), and >4.0 (extremely large) [33]. Additionally, the ICC and the Pearson correlation coefficients (r) were interpreted as follows: <0.10 (trivial), 0.30 (small), 0.50 (moderate), 0.70 (large), 0.90 (very large), 0.99 (nearly perfect), and 1 (perfect) [31]. The data analysis was performed using the SPSS (19.0). The significance adopted was set at .

3. Results

3.1. Identification of HRVT_DMAX1, HRVT_DMAX2, LT₁, and LT₂

No significant differences were observed, neither for absolute nor for relative values, between HRVT_DMAX1 and LT₁ when expressed for speed (; d = 0.45 and ; d = 0.76), lactate (; d = 0.79 and ; d = 0.06), RR (; d = 0.48 and ; d = 0.49), and HR (; d = 0.44 and ; d = 0.51), respectively. In the same way, significant differences were not observed between HRVT_DMAX2 and LT₂ when expressed for lactate (; d = 0.14 and ; d = 0.00) and RR (; d = 0.68 and ; d = 0.68), but on the other hand, significant differences were observed when expressed for speed (; d = 1.21 and ; d = 1.9) and HR (; d = 1.15 and ; d = 1.24), respectively. Further, no significant differences were observed, neither for absolute nor for relative values, between HRVT_DMAX2 and LT₁ when expressed for speed (; d = 0.41 and ; d = 0.75), RR (; d = 0.35 and ; d = 0.70), and HR (; d = 0.33 and ; d = 0.39). Table 1 shows the results of all the methods.

Table 2 shows in detail the results of the Pearson correlation coefficient between the methods expressed as absolute and relative values for speed, lactate, RR, and HR.

Figure 2 shows the magnitude of differences between HRVT_DMAX1 and LT₁, and HRVT_DMAX2 and LT₂. The Bland–Altman and regression analysis showed between HRVT_DMAX1 and LT₁ the bias ± LoA = 0.84 ± 3.7 km·h⁻¹ and SEE = 1.2 km·h⁻¹ (95% CI, 0.9–1.9), and between HRVT_DMAX2 and LT₂ the bias ± LoA = −1.07 ± 1.8 km·h⁻¹ and SEE = 1.1 km·h⁻¹ (95% CI, 0.8–1.6).

(a)

(b)

(c)

(d)

3.2. Reliability of the HRV and [La] Thresholds

In relation to baseline HRV values (423.0 ± 28 ms vs. 425.7 ± 25 ms; ; d = 0.09) and baseline lactate values (1.34 ± 0.4 mmol·L⁻¹ vs. 1.25 ± 0.3 mmol·L⁻¹. ; d = 0.41), no significant differences were observed between test and retest. With regard to values recorded at the exhaustion in MIRT, no significant differences were observed between test and retest for speed (16.4 ± 1.7 km·h⁻¹ vs. 16.6 ± 1.7 km·h⁻¹; ; d = 0.09), HR (192 ± 5 bpm vs. 191 ± 5 bpm; ; d = 0.12), and [La] (9.2 ± 1.8 mmol·L⁻¹ vs. 8.2 ± 1.9 mmol·L⁻¹; ; d = 0.52), respectively.

Significant differences were not observed between test and retest for all thresholds (HRVT_DMAX1 = 12.1 ± 1.4 km·h⁻¹ vs. 12.3 ± 1.5 km·h⁻¹, HRVT_DMAX2 = 12.0 ± 1.2 km·h⁻¹ vs. 11.7 ± 1.4 km·h⁻¹, LT₁ = 11.2 ± 2.1 km·h⁻¹ vs. 12.0 ± 1.5 km·h⁻¹, and LT₂ = 14.1 ± 2.1 km·h⁻¹ vs. 14.2 ± 1.7 km·h⁻¹). Large ICCs were found to HRVT_DMAX1 when expressed in relation to RR and HR (ICC = 0.80-0.81; TEM = 5.1–4.9%, respectively). The Bland–Altman and regression analysis showed the bias ± LoA = 0.84 ± 3.7 km·h⁻¹. In the same way, large ICCs were found to HRVT_DMAX2 when expressed in relation to RR and HR (ICC = 0.80–0.82; TEM = 5.5–5.6%, respectively). Results further suggested a consistent reproducibility for LT₁ and LT₂, since the large to nearly perfect ICCs were showed for speed and HR measures (ICC = 0.90–0.97 and 0.80–0.94; TEM = 6.4–4.9% and 2.8–1.7%, respectively). Table 3 shows all results of ICC and TEM for all the methods (Figure 3).

(a)

(b)

4. Discussion

The main findings of the present study were that the application of the Dmax method on HRV dataset (SD1 and SD2 markers by the Poincaré plot, being HRVT_DMAX1 and HRVT_DMAX2, respectively) enabled the estimation of the LT₁ during a MIRT in male runners. However, the results refute one of the hypotheses of the study, which was that the HRVT_DMAX2 method would estimate the AnT. The relative values found in HRVT_DMAX2 and HRVT_DMAX1 showed values with a greater approximation to the AeT. Consequently, it seems to suggest that when the Dmax method is applied to a HRV dataset extracted by the Poincaré plot, it is possible to identify a transition zone, with an approximation to the AeT, being an important intensity to improve cardiorespiratory and neuromuscular responses of runners. Furthermore, the results also suggested a consistent reproducibility for LT₁ and LT₂ (ICC =0.90 and 0.97, respectively), as well as, moderate to HRVT_DMAX1 (ICC = 0.48) with low bias (−0.18 ± 3.4 km·h⁻¹).

The results of the SD1 marker, which was used to identify the HRVT_DMAX1, showed values of approximately ≈73.4% of the peak speed value, being these intensities related to the AeT [3, 4]. However, the values were slightly above the values reported in previous studies such as Garcia-Tabar et al. [9] analyzing a homogeneous group of professional male world-class road cyclists (≈36–52% W_PEAK), Candido et al. [23] analyzing healthy individuals (≈50–60% W_PEAK), Sales et al. [20] analyzing individuals with type 2 diabetes (≈64–66% VO_2PEAK), and Tulppo et al. [2] investigating complete or not parasympathetic blockade. The HRVT_DMAX1 has elicited signiﬁcant correlation when compared to LT₁ (r = 0.46). The results of the present study in relation to HRVT_DMAX1 are slightly below those found by Garcia-Tabar et al. [9], which used the same marker of PNS (SD1) to estimate the LT (r = 0.66–0.88), although different protocols and methodologies were applied. However, it is important to note that only a study of our knowledge by Nascimento et al. [11] used the same ergometer when analyzing HRV indices by the Poincaré plot and [La], in case the treadmill, and all other studies used a cycle ergometer. Consequently, greater comparisons are limits due mainly to the specificity and differences in the movement gesture as well as in the recruitment of motor units [35, 36].

It is important to note that the SD1 marker, which presents the prevalence of PNS activity, correlates with indices representing high-frequency bands (HF), such as in the Fourier Transform, when analyzing the behavior of HRV by frequency domain [1]. In addition, studies suggest that the respiratory pattern has a strong effect on the HF-HRV bands, both at rest and at exercise [37–39]. During exercise in heavy domain occurs an increase in respiratory frequency with a constant final volume, triggering a mechanical effect on the sinus node, concomitant with an increase in HF. This can be demonstrated in previous studies which identified changes in HF behavior below and above of ventilatory threshold [37].

The SD2 marker used in the present study to identify the HRVT_DMAX2 showed significant differences in relation to LT₂ when expressed for absolute and relative values of speed and HR, but not when expressed to lactate and RR, respectively. On the other hand, it is important to note that HRVT_DMAX2 has elicited moderate coefficient values (r = 0.48) and signiﬁcant coefficient values (r = 0.71) when compared to LT₂, with the values expressed for speed and HR, respectively. This variation in the correlations may be explained in part by the size and heterogeneous characteristic of the sample (body fat coefficient of variation (CV) is 33%). The values found in relation to the HRVT_DMAX2 were approximately ≈72.9% of the peak speed value, being these intensities related to the AeT [3, 4]. Previous studies showed a breakpoint in SD2 to intensities of approximately ≈80%–90% VO_2MAX during maximal progressive cycling test [2] and intensities of approximately ≈86.1% of the peak speed value during MIRT [11]. However, as previously mentioned in relation to the possibility of identifying a single breakpoint by the Dmax method, these results suggest that the SD2 marker shows a nonlinearity behavior during a MIRT. Therefore, in addition to a breakpoint in intensity close to the AnT, as already demonstrated in previous studies [2, 11], a significant change point in SD2 also occurs at near intensities related to the AeT.

In the present study, the approximation of HRVT_DMAX2 with the AeT, possibly, is due to the use of the Dmax method. This method allows the identification of only one breakpoint, although it is a nonsubjective method when compared to the visual method, for example. However, there are questions concerning the determination of the breakpoint by the Dmax method, in relation to the amount of values used in the model construction, being that the initial and final values can influence and compromise greater inferences when comparing the identification of the AeT or AnT [40].

Perhaps the use of different mathematical models [40] on HRV dataset, even with the possibility of submaximal tests [9], could allow greater accuracy in the estimation of AeT and AnT in different situations and populations involved. Moreover, both HRVT_DMAX1 and HRVT_DMAX2 are methods relatively simple to analyze and do not require a fixed number of RR intervals nor long recording periods [41], which facilitates the evaluation during incremental exercise. The usefulness of the HRVT_DMAX1 and HRVT_DMAX2 to estimate the AeT and determine aerobic capacity to prescribe exercise training intensities in sports performance and rehabilitation programs, from SD1 and SD2 values, is of great interest. The HRVTs may be objectively and noninvasively determined and are of lower cost than lactate or ventilatory threshold assessments, where blood lactate or gas analysis equipment is required as well as specialized professionals.

The LT₁ and LT₂ methods demonstrated a high level of reliability (ICC = 0.90 and 0.97; ; TEM = 0.75 and 0.40 km·h⁻¹; TEM = 6.4 and 2.8%, respectively). These results presented lower values than previous studies using similar markers for determination of LT₁ and LT₂ (TEM = 2.8 and 3.6%, respectively) [42]. Bland and Altman presented low values for difference between the test and retest as a function of their mean (bias ± LoA = −0.73 ± 2.0 km·h⁻¹; bias ± LoA = −0.08 ± 1.1 km·h⁻¹, respectively). Regarding the HRVT_DMAX1, significant values were found in relation to RR and HR (ICC = 0.80 and 0.81; , respectively), but on the other hand, moderate values were found when expressed in relation to speed (ICC = 0.48), being below those found in a previous study (ICC = 0.73) using the same marker for identification of AeT [15], but through visual identification instead of the Dmax method. Regarding the TEM, results were slightly above the values found in the aforementioned study (TEM = 8.5%). Bland and Altman presented low values for difference between the test and retest as a function of their mean (bias ± LoA = −0.18 ± 3.4 km·h⁻¹). Already to the HRVT_DMAX2, significant values were found when expressed in relation to lactate (ICC = 0.60; ) and RR and HR (ICC = 0.80 and 0.82; , respectively), with low values only when expressed in relation to speed (ICC = 0.30; ), corroborating previous studies that used HRV markers by the Dmax method [23].

Due to its low cost, noninvasive nature, and high applicability, HRVTs is an important framework for researchers, trainers, and race practitioners. In the present study, its simplified and nonsubjective identification by the Dmax method suggests the possibility of planning a training program with a safe intensity in a metabolic transition zone, being slightly above the values found for AeT.

The present study certainly has some limitations that must be considered in the analysis of the results and their applicability. Aspects such as heterogeneous characteristic of the sample may try to explain in part the variation of correlation values between the different methods, since previous studies report that the greater the heterogeneity of a group, the greater the magnitudes of correlation [32]. Another fact may be just the recruitment of male runners and the need to perform a MIRT. In addition, MLSS was not used as a gold standard, but LT₂ had a good approximation of MLSS in runners [30]. Therefore, it is suggested to carry out studies with different characteristics of samples, such as gender, age, and levels of training.

5. Conclusion

To the best of our knowledge, this is the first report on the application of the Dmax method on a HRV dataset to identify the lactate thresholds, during a MIRT. The results of this study showed that the Dmax method applied over a set of HRV during a MIRT in male recreational long-distance runners allowed the identification of HRVTs approaching the AeT. The HRVTs are of low cost, noninvasive nature, and high applicability. A limiting factor for the interpretation of the data is the recruitment of only males and trained, which do not allow the generalization of results to different populations. Thus, further studies are needed to confirm the reproducibility of HRVT as well as its use in different protocols, genders, age groups, and levels of training.

6. Disclosure

The level of evidence is Level II Study of Diagnostic Test.

Abbreviations

HRV:	Heart rate variability
VO_2MAX:	Maximum oxygen uptake
AeT:	Aerobic threshold
PNS:	Parasympathetic nervous system
SNS:	Sympathetic nervous system
AnT:	Anaerobic threshold
HRVT:	Heart rate variability threshold
SD1:	The standard deviation of instantaneous RR intervals
SD2:	The continuous long-term RR intervals
HR:	Heart rate
MIRT:	Maximal incremental running test
LA:	Lactate concentrations.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this manuscript.

Acknowledgments

This study was financed in part by the Coordenacão de Aperfeiçoamento de Pessoal de Nível Superior—Brasil (CAPES)—Finance Code 001. Dr. Paulo Cesar do Nascimento Salvador also acknowledges the support by grants from CNPq, Conselho Nacional de Desenvolvimento Científico e Tecnológico—Brasil (154191/2018-3).

References

M. Malik, “Task force of the European society of cardiology and the north American society of pacing and electrophysiology. Heart rate variability: standards of measurement, physiological interpretation and clinical use,” Circulation, vol. 93, no. 5, pp. 1043–1065, 1996.
View at: Google Scholar
M. P. Tulppo, T. H. Makikallio, T. E. Takala, T. Seppanen, and H. V. Huikuri, “Quantitative beat-to-beat analysis of heart rate dynamics during exercise,” American Journal of Physiology-Heart and Circulatory Physiology, vol. 271, no. 1, pp. H244–H252, 1996.
View at: Publisher Site | Google Scholar
O. Faude, W. Kindermann, and T. Meyer, “Lactate threshold concepts: how valid are they?” Sports Medicine, vol. 39, no. 6, pp. 469–490, 2009.
View at: Publisher Site | Google Scholar
J. S. Skinner and T. M. McLellan, “The transition from aerobic to anaerobic metabolism,” Research Quarterly for Exercise and Sport, vol. 51, no. 1, pp. 234–248, 1980.
View at: Publisher Site | Google Scholar
G. K. Karapetian, H. J. Engels, and R. J. Gretebeck, “Use of heart rate variability to estimate LT and VT,” International Journal of Sports Medicine, vol. 29, no. 8, pp. 652–657, 2008.
View at: Publisher Site | Google Scholar
N. L. Jones and R. E. Ehrsam, “The anaerobic threshold,” Exercise and Sport Sciences Reviews, vol. 10, no. 1, pp. 49–83, 1982.
View at: Publisher Site | Google Scholar
I. Garcia-Tabar and E. M. Gorostiaga, “A “blood relationship” between the overlooked minimum lactate equivalent and maximal lactate steady state in trained runners. Back to the old days?” Frontiers in Physiology, vol. 9, p. 1034, 2018.
View at: Publisher Site | Google Scholar
I. Llodio, E. M. Gorostiaga, I. Garcia-Tabar, C. Granados, and L. Sánchez-Medina, “Estimation of the maximal lactate steady state in endurance runners,” International Journal of Sports Medicine, vol. 37, no. 7, pp. 539–546, 2016.
View at: Publisher Site | Google Scholar
I. Garcia-Tabar, L. Sánchez-Medina, J. F. Aramendi, M. Ruesta, J. Ibañez, and E. M. Gorostiaga, “Heart rate variability thresholds predict lactate thresholds in professional world-class road cyclists,” Journal of Exercise Physiology Online, vol. 16, no. 5, pp. 38–50, 2013.
View at: Google Scholar
M. Buchheit, R. Solano, and G. P. Millet, “Heart-rate deflection point and the second heart-rate variability threshold during running exercise in trained boys,” Pediatric Exercise Science, vol. 19, no. 2, pp. 192–204, 2007.
View at: Publisher Site | Google Scholar
E. M. F. Nascimento, M. A. P. D. M. Kiss, T. M. Santos, M. Lambert, and F. O. Pires, “Determination of lactate thresholds in maximal running test by heart rate variability data set,” Asian Journal of Sports Medicine, vol. 8, no. 3, Article ID e58480, 2017.
View at: Publisher Site | Google Scholar
F. Cottin, C. Medigue, P. Lopes, P. M. Lepretre, R. Heubert, and V. Billat, “Ventilatory thresholds assessment from heart rate variability during an incremental exhaustive running test,” International Journal of Sports Medicine, vol. 28, no. 4, pp. 287–294, 2007.
View at: Publisher Site | Google Scholar
F. Cottin, P. M. Lepretre, P. Lopes, Y. Papelier, C. Medigue, and V. Billat, “Assessment of ventilatory thresholds from heart rate variability in welltrained subjects during cycling,” International Journal of Sports Medicine, vol. 27, no. 12, pp. 959–967, 2006.
View at: Publisher Site | Google Scholar
C. Liu, C. Liu, P Shao et al., “Comparison of different threshold values r for approximate entropy: application to investigate the heart rate variability between heart failure and healthy control groups,” Physiological Measurement, vol. 32, no. 2, pp. 167–180, 2011.
View at: Publisher Site | Google Scholar
F. Novelli, J. de Araújo, G. Tolazzi, G. Tricot, G. Arsa, and L. Cambri, “Reproducibility of heart rate variability threshold in untrained individuals,” International Journal of Sports Medicine, vol. 40, no. 2, pp. 95–99, 2018.
View at: Publisher Site | Google Scholar
J. Cassirame, N. Tordi, N. Fabre, S. Duc, F. Durand, and L. Mourot, “Heart rate variability to assess ventilatory thresholds in ski-mountaineering,” European Journal of Sport Science, vol. 15, no. 7, pp. 615–622, 2015.
View at: Publisher Site | Google Scholar
L. Mourot, N. Fabre, A. Savoldelli, and F. Schena, “Second ventilatory threshold from heart rate variability: valid when the upper body is involved?” International Journal of Sports Physiology and Performance, vol. 9, no. 4, pp. 695–701, 2014.
View at: Publisher Site | Google Scholar
R. Di Michele, G. Gatta, A. Di Leo et al., “Estimation of the anaerobic threshold from heart rate variability in an incremental swimming test,” Journal of Strength and Conditioning Association, vol. 26, no. 11, pp. 3059–3066, 2012.
View at: Publisher Site | Google Scholar
S. Quinart, L. Mourot, V. Negre et al., “Ventilatory thresholds determined from HRV: comparison of 2 methods in obese adolescents,” International Journal of Sports Medicine, vol. 35, no. 3, pp. 203–208, 2014.
View at: Publisher Site | Google Scholar
M. M. Sales, C. S. Campbell, P. K Morais et al., “Noninvasive method to estimate anaerobic threshold in individuals with type 2 diabetes,” Diabetology & Metabolic Syndrome, vol. 3, no. 1, pp. 1–8, 2011.
View at: Publisher Site | Google Scholar
L. Mourot, M. Bouhaddi, S. Perrey, J. D. Rouillon, and J. Regnard, “Quantitative Poincare plot analysis of heart rate variability: effect of endurance training,” European Journal of Applied Physiology, vol. 91, no. 1, pp. 79–87, 2003.
View at: Publisher Site | Google Scholar
B. Casadei, S. Cochrane, J. Johnston, J. Conway, and P. Sleight, “Pitfalls in the interpretation of spectral analysis of the heart rate variability during exercise in humans,” Acta Physiologica Scandinavica, vol. 153, no. 2, pp. 125–131, 1995.
View at: Publisher Site | Google Scholar
N. Candido, N. M. Okuno, C. C. da Silva, F. A. Machado, and F. Y. Nakamura, “Reliability of the heart rate variability threshold using visual inspection and Dmax methods,” International Journal of Sports Medicine, vol. 36, no. 13, pp. 1076–1080, 2015.
View at: Publisher Site | Google Scholar
F. Conconi, M. Ferrari, P. G. Ziglio, P. Droghetti, and L. Codeca, “Determination of the anaerobic threshold by a noninvasive field test in runners,” Journal of Applied Physiology, vol. 52, no. 4, pp. 869–873, 1982.
View at: Publisher Site | Google Scholar
B. Cheng, H. Kuipers, A. C. Snyder, H. A. Keizer, A. Jeukendrup, and M. Hesselink, “A new approach for the determination of ventilatory and lactate thresholds,” International Journal of Sports Medicine, vol. 13, no. 7, pp. 518–522, 1992.
View at: Publisher Site | Google Scholar
R. H. Morton, S. R. Stannard, and B. Kay, “Low reproducibility of many lactate markers during incremental cycle exercise,” British Journal of Sports Medicine, vol. 46, no. 1, pp. 64–69, 2012.
View at: Publisher Site | Google Scholar
D. Harriss and G. Atkinson, “Ethical standards in sport and exercise science research: 2016 update,” International Journal of Sports Medicine, vol. 36, no. 14, pp. 1121–2112, 2015.
View at: Publisher Site | Google Scholar
H. Heck, A. Mader, G. Hess, S. Mucke, R. Muller, and W. Hollmann, “Justification of the 4-mmol/l lactate threshold,” International Journal of Sports Medicine, vol. 6, no. 3, pp. 117–130, 1985.
View at: Publisher Site | Google Scholar
H. H. Dickhuth, L. Yin, A. Niess et al., “Ventilatory, lactate derived and catecholamine thresholds during incremental treadmill running: relationship and reproducibility,” International Journal of Sports Medicine, vol. 20, no. 2, pp. 122–127, 1999.
View at: Publisher Site | Google Scholar
K. M. Souza, T. Grossl, R. J. Babel Junior, R. D. de Lucas, V. P. Costa, and L. G. A. Guglielmo, “Maximal lactate steady state estimated by different methods of anaerobic threshold,” Revista Brasileira de Cineantropometria e Desempenho Humano, vol. 14, no. 3, pp. 264–275, 2012.
View at: Publisher Site | Google Scholar
W. G. Hopkins, “Spreadsheets for analisys of validity and reliability,” Sportscience, vol. 19, pp. 36–42, 2015.
View at: Google Scholar
W. G. Hopkins, S. W. Marshall, A. M. Batterham, and J. Hanin, “Progressive statistics for studies in sports medicine and exercise science,” Medicine & Science in Sports & Exercise, vol. 41, no. 1, pp. 3–13, 2009.
View at: Publisher Site | Google Scholar
J. Cohen, Statiscal Power Analysis for the Behavioral Sciences, Lawrance Erlbaum, Mahwah, NJ, USA, 1986.
J. M. Bland and D. G. Altman, “Statistical methods for assessing agreement between two methods of clinical measurement,” The Lancet, vol. 327, no. 8476, pp. 307–310, 1986.
View at: Publisher Site | Google Scholar
I. E. Klein, J. B. White, and S. R. Rana, “Comparison of physiological variables between the elliptical bicycle and run training in experienced runners,” Journal of Strength and Conditioning Research, vol. 30, no. 11, pp. 2998–3006, 2016.
View at: Publisher Site | Google Scholar
G. P. Millet, V. E. Vleck, and D. J. Bentley, “Physiological differences between cycling and running: lessons from triathletes,” Sports Medicine, vol. 39, no. 3, pp. 179–206, 2009.
View at: Publisher Site | Google Scholar
F. Cottin, C. Medigue, P. M. Lepretre, Y. Papelier, J. P. Koralsztein, and V. Billat, “Heart rate variability during exercise performed below and above ventilatory threshold,” Medicine & Science in Sports & Exercise, vol. 36, no. 4, pp. 594–600, 2004.
View at: Publisher Site | Google Scholar
S. Laborde, E. Mosley, and J. F. Thayer, “Heart rate variability and cardiac vagal tone in psychophysiological research—recommendations for experiment planning, data analysis, and data reporting physiol,” Frontiers in Psychology, vol. 8, p. 213, 2017.
View at: Publisher Site | Google Scholar
T. E. Brown, L. A. Beightol, J. Koh, and D. L. Eckberg, “Important influence of respiration on human R-R interval power spectra is largely ignored,” Journal of Applied Physiology, vol. 75, no. 5, pp. 2310–2317, 1993.
View at: Publisher Site | Google Scholar
S. Chalmers, A. Esterman, R. Eston, and K. Norton, “Standardization of the Dmax method for calculating the second lactate threshold,” International Journal of Sports Physiology and Performance, vol. 10, no. 7, pp. 921–926, 2015.
View at: Publisher Site | Google Scholar
G. R. H. Sandercock and D. A. Brodie, “The use of heart rate variability measures to assess autonomic control during exercise,” Scandinavian Journal of Medicine and Science in Sports, vol. 16, no. 5, pp. 302–313, 2006.
View at: Publisher Site | Google Scholar
J. A. A. C. Heuberger, P. Gal, F. E. Stuurman, W. A. S. de Muinck Keizer, Y. Mejia Miranda, and A. F. Cohen, “Repeatability and predictive value of lactate threshold concepts in endurance sports,” PLoS One, vol. 13, no. 11, Article ID e0206846, 2018.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2019 Eduardo Marcel Fernandes Nascimento et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

2375

Downloads

1385

Citations

Journal of Sports Medicine

Applicability of Dmax Method on Heart Rate Variability to Estimate the Lactate Thresholds in Male Runners

Abstract

1. Introduction

2. Materials and Methods

2.1. Participants

2.2. Study Design

2.3. Maximal Incremental Running Test Protocol

2.4. Determination of Aerobic and Anaerobic Thresholds

2.4.1. Heart Rate Variability Thresholds

2.4.2. Lactate Thresholds

2.5. Statistical Procedures

3. Results

3.1. Identification of HRVTDMAX1, HRVTDMAX2, LT1, and LT2

3.2. Reliability of the HRV and [La] Thresholds

4. Discussion

5. Conclusion

6. Disclosure

Abbreviations

Data Availability

Conflicts of Interest

Acknowledgments

References

Copyright

3.1. Identification of HRVT_DMAX1, HRVT_DMAX2, LT₁, and LT₂