Modeling the diving bradycardia: Toward an "oxygen-conserving breaking point"?

Purpose: Although it has been demonstrated that the exponential decay model fits the heart rate (HR) kinetics in short static breath holding (BH), this model might be inaccurate when BH is maintained for several minutes. The aim of this study was to build a new meaningful model to quantify HR kinetics during prolonged static BH.
Methods: Nonlinear regression analysis was used to build a model able to quantify the beat-to-beat HR reduction kinetics observed in prolonged static BH performed both in air and in immersed condition by 11 trained breath-hold divers. Dynamic changes in cardiac autonomic regulation through heart rate variability indices [root mean square of successive difference of R-R intervals (RMSSD); short-term fractal scaling exponent: (DFAα1)] and peripheral oxygen saturation (SpO2) were also analyzed to strengthen the model.
Results: The tri-phasic model showed a sharp exponential drop in HR immediately followed by a slight linear rise up until a breaking point preceding a linear drop in HR. The breaking points had similar level of SpO2 whether in air or in immersed condition (95.1 ± 2.1 vs. 95.2 ± 3.0 %, respectively; P = 0.49), and the subsequent linear drop in HR was concomitant with a shift in cardiac autonomic regulation in air (RMSSD: +109.0 ± 47.8 %; P < 0.001; DFAα1: -18.0 ± 17.4 %; P < 0.05) and in immersion (RMSSD: +112.6 ± 55.8 %; P < 0.001; DFAα1: -26.0 ± 12 %; P < 0.001).
Conclusion: In addition to accurately fitting the HR kinetics, the most striking finding is an "oxygen-conserving breaking point" highlighted by the model, which might be interpreted as unique adaptive feature against hypoxic damages in the human diving bradycardia.