- 1PROFITH “PROmoting FITness and Health Through Physical Activity” Research Group, Department of Physical and Sports Education, Faculty of Sport Sciences, University of Granada, Granada, Spain
- 2Department of Rehabilitation Sciences, KU Leuven – University of Leuven, Leuven, Belgium
- 3Center for Cognitive and Brain Health, Department of Psychology, Northeastern University, Boston, MA, United States
- 4IRyS Group, School of Physical Education, Pontificia Universidad Católica de Valparaíso, Valparaíso, Chile
- 5Andalusian Centre of Sport Medicine (CAMD), Junta de Andalucía, Granada, Spain
- 6Department of Medical Physiology, School of Medicine, University of Granada, Granada, Spain
- 7Department of Public Health, Faculty of Medicine and Health Sciences, Ghent University, Ghent, Belgium
- 8Faculty of Physical Education and Physiotherapy, Opole University of Technology, Opole, Poland
- 9Department of Cardiology, University Hospital in Opole, University of Opole, Opole, Poland
- 10Department of Biosciences and Nutrition, Karolinska Institutet, Huddinge, Sweden
Cardiac autonomic function can be quantified through mean heart rate (HR) or heart rate variability (HRV). Numerous studies have supported the utility of different HRV parameters as indicators of cardiorespiratory fitness (CRF). However, HR has recently shown to be a stronger predictor of CRF than HRV in healthy young adults, yet these findings need to be replicated, in other age groups such as children. Therefore, this study aimed: (1) to study the associations between indicators of cardiac autonomic function (HR, standard and corrected HRV parameters) and CRF in overweight/obese children; and (2) to test which of the two indicators (i.e., HR or HRV) is a stronger predictor of CRF. This study used cross-sectional baseline data of 107 overweight/obese children (10.03 ± 1.13 years, 58% boys) from the ActiveBrains project. Cardiac autonomic indicators were measured with Polar RS800CX®. CRF was assessed using a gas analyzer while performing a maximal incremental treadmill test. Correlations and stepwise linear regressions were performed. Mean HR and standard HRV parameters (i.e., pNN50, RMSSD, and SDNN) were associated with CRF (r coefficients ranging from -0.333 to 0.268; all p ≤ 0.05). The association of HR with CRF persisted after adjusting for sex, peak height velocity (PHV), adiposity moderate-to-vigorous physical activity, energy intake and circadian-related variable intradaily variability of activity patterns whilst for HRV parameters (i.e., pNN50, RMSSD, and SDNN) disappeared. Stepwise linear regression models entering HR and all HRV parameters showed that mean HR was the strongest predictor of CRF (β = -0.333, R2 = 0.111, p < 0.001). Standard and corrected HRV parameters did not provide additional value to the coefficient of determination (all p > 0.05). Our findings suggest that HR is the strongest indicator of CRF. It seems that quantification of HRV parameters in time and frequency domain do not add relevant clinical information about the cardiovascular health status (as measured by CRF) in overweight/obese children beyond the information already provided by the simple measure of HR.
Introduction
A number of comorbidities linked to childhood obesity could be explained by their altered cardiac autonomic functions (Thayer and Lane, 2007; Zhou et al., 2012). Particularly, overweight/obese children and adults usually present a decreased parasympathetic activity (Gutin et al., 2000; Nagai et al., 2003; Rodríguez-Colón et al., 2011; Triggiani et al., 2017) which, in turn, is related to a higher risk of cardiovascular disease (CVD) (Lahiri et al., 2008). Therefore, the quantification of cardiac autonomic function is a matter of interest in overweight/obese populations.
Heart rate (HR) and heart rate variability (HRV) have been proposed as indicators of the cardiac autonomic nervous system functioning (Koizumi et al., 1985; Task Force, 1996; Lahiri et al., 2008). Specifically, HR reflects the frequency of heart beats per minute. HRV refers to the variability in time intervals between consecutive heart beats, namely between the R peaks registered in the electrocardiogram (ECG) (Task Force, 1996). Both HR and HRV have been used as predictors of CVD, showing that higher HR and lower HRV at rest are related to a higher CVD risk (Kannel et al., 1987; Thayer and Lane, 2007; Trimmel et al., 2015; Dilaveris and Tousoulis, 2018).
A good indicator of the CVD risk in adults is cardiorespiratory fitness (CRF) (Ortega et al., 2016; Ross et al., 2016; Tikkanen et al., 2017), which is also known as a powerful marker of health in children and adolescents (Ortega et al., 2008). In this context, HR is inversely associated with CRF (Kang et al., 2017), but based on a recent systematic review, associations between HRV and CRF are inconsistent in children and adolescents (Oliveira et al., 2017). For example, some studies have found positive associations between HRV and CRF in obese youth (Gutin et al., 2005; Michels et al., 2013; Da Silva et al., 2014), whereas others did not find differences in HRV parameters between 3 groups of healthy adolescents with different levels of CRF (Brunetto et al., 2005). In fact, Oliveira et al. concluded that several methodological differences between studies could contribute to these differences and that relevant studies in pediatric population are currently lacking.
Importantly, HRV reveals a non-linear inverse association with average HR (or direct association with average R-R interval), and consequently some clinical value of HRV must originate from average HR, – to elucidate this issue one should liberate HRV from the influence of HR (Sacha, 2013, 2014a; Sacha et al., 2013a; Sacha and Pluta, 2008). To overcome this HRV dependence on HR, a correction procedure has been proposed which is based on calculating the ratios between HRV parameters and different powers of their corresponding mean R-R interval (Sacha et al., 2013b). Of note, by employing this procedure, the HRV prediction capacity after myocardial infarction has been improved (Sacha, 2014b) and it was also possible to show that HR is a better predictor of CRF than HRV in young adults (Grant et al., 2013), although little is known on how it works in overweight/obese children.
Thus, the aims of this study were: (1) to explore the associations between indicators of cardiac autonomic function (HR, standard and corrected HRV parameters) and CRF; and (2) to test which of the two indicators (i.e., HR or HRV) is a stronger predictor of CRF in overweight/obese children.
Materials and Methods
Participants and Study Design
This cross-sectional study used baseline data from the ActiveBrains project1 (Clinical Trial: NCT02295072). A total of 107 overweight/obese children (10.03 ± 1.13 years, 58% boys) were included in the present study. Information about the methodology and protocol can be found elsewhere (Cadenas-Sánchez et al., 2016). This study was conducted according to the Declaration of Helsinki. The protocol was approved by the Committee for Research Involving Human Subjects at the University of Granada (Reference: 848, February 2014). All parents had received information about the study and gave written informed consent in accordance with the Declaration of Helsinki.
Anthropometry and Weight Status
Body weight and height were measured to the nearest 0.1 kg and cm with an electronic scale and a stadiometer, respectively (SECA instruments, Germany, Ltd). Body mass index (BMI) was calculated as kg/m2 and the participants were classified as overweight, mild obesity or severe/morbid obesity according to the sex-and-age specific international BMI standards (World Obesity Federation, formerly named International Obesity Task Force) (Cole and Lobstein, 2012; Bervoets and Massa, 2014). Body fat percentage was assessed with dual-energy X-ray Absorptiometry (DXA, discovery densitometer from Hologic). All DXA scans and analyses were performed using the GE encore software (version 4.0.2) and were completed following the same protocol by the same evaluator. The positioning of the participants and the analyses of the results were undertaken following recommendations from the International Society of Clinical Densitometry (Crabtree et al., 2014).
Physical Activity
Physical activity (PA) levels were estimated from the triaxial accelerometers ActiGraph GT3X+ (ActiGraph, Pensacola, FL, United States). Accelerometers were initialized to capture and store accelerations at 100 Hz and placed on non-dominant wrist for seven consecutive days. We processed the raw accelerations in the GGIR R package (Van Hees et al., 2013). We considered every participant wearing the accelerometer a minimum of 16 h per day during at least 3 weekdays and 1 weekend day as recommended elsewhere (Migueles et al., 2017, 2019). The daily time spent in moderate-to-vigorous PA intensity (i.e., >3 metabolic equivalents) was estimated using previously proposed cut-points for children (Hildebrand et al., 2014).
Energy Intake
Energy intake was obtained by means of two non-consecutive 24-h recalls during a time span of a week by trained dietitians-nutritionists. The presence of parents or legal guardians was obligatory for the collection of dietary data due to the difficulty for children to remember recipes or amounts of foods. A book with pictures of different food servings and sizes was used to help the participants to estimate the amount of food consumed. The Diet software (Xyris Software, Brisbane, Australia), supported by the Spanish Association of Dietetics and Nutritionists, was used to calculated dietary data.
Circadian-Related Variable Based on Fragmentation of Daytime Activity Patterns: Intradaily Variability
We calculated non-parametric indexes described by Van Someren et al. (1999) to characterize activity patterns. The circadian-related variable selected was fragmentation of daytime activity, measured by intradaily variability (IV), due to has been related to cardiorespiratory fitness and adiposity in adolescents (Garaulet et al., 2016). This circadian-related variable indicates the frequency of changes between high and low activity. Its values oscillated between 0, when the wave was perfectly sinusoidal, and 2, when the wave was as Gaussian noise.
Heart Rate and Heart Rate Variability
Participants were placed supine for 10 min in a quiet and comfortable room between 9 a.m. and 12 p.m. Supine position has shown a higher reliability for the HRV measurement in children than sitting or standing (Silva et al., 2016). The POLAR RS800CX (Polar Electro Oy Inc., Kempele, Finland) recorded HRV during 10 min at a sampling frequency of 1000 Hz. This HR monitor is valid (Gamelin et al., 2008) and reliable when compared with electrocardiography (Vasconcellos et al., 2015) for the HRV assessment in children and adolescents. Participants were encouraged to keep relaxed, breathe normally and do not speak or move during the evaluation.
The indicators of cardiac autonomic function were HR and HRV parameters. For HRV analyses, we derived time and frequency-domain HRV parameters from the normal R-R intervals after excluding the extreme values using an automatized “low filter” in the Kubios software (HRV analysis, University of Eastern Finland) (Niskanen et al., 2004; Tarvainen et al., 2014). The R-R intervals series were detrended using the smoothness prior method with alpha set at 500 and a cubic interpolation at the default rate of 4 Hz. Out of the 10 min recorded, the middle 5 min (i.e., from 3 to 8 min) were visually checked for quality (i.e., normal distribution of the R-R intervals, no large R-R interval outliers and R-R intervals equidistance and minimal variation), and a different period of 5-min was selected when necessary. Time- and frequency-domain HRV parameters were selected based on the Guidelines of Task Force of The European Society of Cardiology and The North American Society of Pacing and Electrophysiology (Task Force, 1996). In the time-domain, we computed the squared root of the mean of the sum of the squares of successive normal R-R interval differences (RMSSD) and the percentage number of pairs of adjacent normal R-R intervals differing by more than 50 ms in the entire recording (pNN50), both as indexes of parasympathetic activity. We also computed the standard deviation of all normal R-R intervals (SDNN). In the frequency-domain, we performed spectral analyses using the non-parametric fast fourier transformation algorithm (FFT), with Welch’s periodogram method (i.e., 50% overlap Hanning window as pre-processing technique and calculating area under the curve with an integration). We derived total power of HRV spectrum (TP), the power in the high (HF) and the low frequency (LF) (TP: 0–0.4 Hz; HF: 0.15–0.4 Hz; LF: 0.04–0.15 Hz) in absolute units (ms2). HF and the LF/HF ratio were used in the analyses as indicators of parasympathetic activity and sympatho-vagal balance, respectively (Task Force, 1996), although physiological significance of LF and LF/HF ratio is not clear (Goldstein et al., 2011; Billman, 2013).
To remove the HRV dependence on HR, we calculated the corrected HRV parameters proposed by Sacha et al. (2013b) based in two assumptions: (i) If HRV parameters were negatively correlated with HR, correction procedure consisted in calculating ratios between HRV parameters and different powers of their corresponding mean normal R-R interval; (ii) if HRV parameters were positively correlated with HR, the correction procedure was performed by multiplying HRV parameters by the adequate powers of mean normal R-R interval as follows: RMMSDc = RMSSD/meanRR3.7, pNN50c = pNN50/meanRR5.3, SDNNc = SDNN/meanRR2.8, FFT TPc = TP/meanRR6.3, FFT HFc = HF/meanRR6.3, FFT LFc = LF/meanRR5.0, FFT LF/HFc = LF/HF × meanRR1.8.
Cardiorespiratory Fitness
Cardiorespiratory fitness was assessed using a gas analyzer (General Electric Corporation) while performing a maximal incremental treadmill test (HP-Cosmos ergometer) adapted for low-fit children (Davis et al., 2012). The incremental test consisted of walking on a treadmill at a constant speed (4.8 Km/h) starting at a 6% slope with grade increments of 1% every minute until volitional exhaustion. Children were highly encouraged to walk as long as possible. Maximal peak oxygen consumption (VO2peak, ml/kg/min), HR (beats/min) and respiratory exchange ratio (RER) were continuously recorded each 10 s, whilst the rating of perceived exertion (RPE) scale was reported at the end of each 1 min stage using children’s OMNI scale ranging from 0 to 10 (Suminski et al., 2008). VO2max was confirmed when meeting 3 out of 4 following criteria: achieving >85% of aged-predicted HRmax, a RER of ≥1.10, volitional fatigue (i.e., >8 points in the OMNI scale) and a plateau in the oxygen consumption during the last two exercise work rates (<2.0 ml/kg/min). We performed analysis with the complete sample that performed the treadmill test (VO2peak; n = 107) and with the sub-sample of children that met at least 3 out of 4 criteria described above (VO2max; n = 75). As conclusions were the same in both occasions, we presented here the results analyzing the complete sample with VO2peak (n = 107).
Demographic Variables
Peak height velocity was assessed as an accurate and discriminant measure of maturational status (Malina et al., 2015). PHV was obtained from anthropometric variables (weight, height, and seated height) using Moore equations (Moore et al., 2015). Years from PHV offset were calculated by subtracting the age of PHV from the chronological age.
Socioeconomic status was assessed by the parental education level. Both father and mother self-reported their maximum educational level (i.e., elementary school, middle school, high school and university degree completed). Parents responses were then combined as follows: none of the parents had a university degree, one of the parents had a university degree or both parents had a university degree.
Statistical Analyses
Descriptive data are presented as frequency and percentage for categorical variables. Kolmogorov-Smirnov test and a visual inspection of histograms were performed to assess the normality of data distribution. Variables that exhibited normal distribution were presented as mean and standard deviation (mean ± SD). Standard HRV parameters did not exhibit normal distribution and were presented as median and interquartile range. For analytical purpose, normal scores of standard and corrected HRV parameters were calculated using Blom formula (Blom, 1958).
In order to test the dependence between HRV and HR parameters, we performed Spearman correlations between standard and corrected HRV parameters and HR. First, the associations of HR, standard and corrected HRV parameters with CRF were studied using Bivariate Pearson correlations. Second, to determine the strongest predictors of CRF, we performed stepwise linear regressions of HR, standard and corrected HRV parameters with CRF. Predictors were included in two different stepwise models. The first model included HR and standard HRV parameters. The second model included HR and corrected HRV parameters. Lastly, we conducted exploratory analyses to test whether the observed associations were independent of potential confounders, such as sex, PHV, adiposity (body fat percentage), moderate-to-vigorous PA, energy intake and circadian-related variable. Analyses were performed using SPSS version 21.0 (IBM Corp., NY, United States). The significance level was set at 0.05.
Results
Table 1 shows the descriptive characteristics of the participants. Figures 1, 2 present the Spearman correlations between HRV parameters and HR. Significant correlations found between the standard HRV parameters and HR (r ranging from -0.802 to 0.250; all p ≤ 0.009) disappeared after the correction procedure (r ranging from -0.120 to 0.033; p ≥ 0.220).
Figure 1. Scatter plots showing the Spearman correlations of time domain standard (left panels) and corrected (right panels) HRV parameters with mean heart rate.
Figure 2. Scatter plots showing the Spearman correlations of frequency domain standard (left panels) and corrected (right panels) HRV parameters with mean heart rate.
Table 2 presents Pearson correlations of HR, and time and frequency normal scores (according to Blom, 1958) of standard and corrected HRV parameters with CRF. HR was inversely associated with CRF (r = -0.333, p < 0.001). The standard HRV parameters in time-domain (RMSSD, pNN50, and SDNN) were positively associated with CRF (r ranging from 0.209 to 0.268, all p < 0.050), whereas frequency-domain parameters were not significantly associated with CRF (p < 0.05). After correction for the prevailing HR, the RMSSD, pNN50 and SDNN completely lost their association with CRF, i.e., none of the corrected HRV parameter correlated with CRF (all p > 0.05).
Table 2. Pearson bivariate correlations of HR, standard and corrected HRV parameters with cardiorespiratory fitness (CRF) in overweight/obese children.
Table 3 shows stepwise linear regressions of HR and normal scores of standard and corrected HRV parameters with CRF. Out of all variables entered in the model, only HR was selected in the stepwise model as a significant predictor of CRF (R2 = 0.111, p < 0.001). The associations of HR with CRF remained significant (β = -0.169, p = 0.030), after adjustment for potential confounders (i.e., sex, PHV, adiposity [body fat percentage], moderate-to-vigorous PA, energy intake and circadian-related variable); whereas the associations of pNN50 and the rest of the HRV parameters with CRF became/remained non-significant (p > 0.05) (Supplementary Table S1).
Table 3. Stepwise linear regressions of HR, standard and corrected HRV parameters with cardiorespiratory fitness (CRF) in overweight/obese children.
Discussion
The main findings of this study were that: (1) HR showed stronger associations with CRF compared with standard and corrected HRV parameters in either time or frequency domains, regardless of sex, maturation, adiposity, moderate-to-vigorous PA, energy intake and circadian-related variable; (2) HR, was the only independent parameter associated with CRF when all the HRV parameters were included in the model. These results suggest that the complex concept of HRV does not seem to add clinically relevant information about the cardiovascular health status, as measured by CRF, in overweight/obese children beyond the information already provided by the “simple” measure of HR.
We found positive associations between standard HRV parameters in time domain (i.e., RMSSD, pNN50, and SDNN) and CRF, with pNN50 being the HRV parameter more strongly associated. In several studies, parasympathetic activity quantified through standard HRV parameters in time and frequency domains (i.e., RMSSD, pNN50, and HF) has been shown to be positively associated with CRF in children and obese adolescents (Gutin et al., 2005; Michels et al., 2013; Da Silva et al., 2014; Redón et al., 2016). However, a recent systematic review showed controversial associations between standard HRV parameters, used as indicators of parasympathetic activity, and CRF in children and adolescents (Oliveira et al., 2017). The authors concluded that the main limitations of most of the studies were the lack of a detailed explanation of HRV protocol, the use of different protocols to measure CRF and the lack of adjustment for relevant confounders (Oliveira et al., 2017).
In this context, it has been reported that HRV was significantly associated with CRF independently of age, race, sex and adiposity in adolescents (Gutin et al., 2005). However, in our study we show that the associations between HRV and CRF was not significant after controlling for potential confounders (sex, maturation, adiposity, moderate-to-vigorous physical activity, energy intake, circadian-related variable and HR). These contradictory findings can be explained by different factors: participants of different ranges of age (adolescents 14–18 years vs. overweight/obese children 8–11 years); use of different HRV protocol and instrument (ECG vs. HR monitor) and different methodology of CRF assessment (VO2 at HR of 170 bpm vs. maximal incremental treadmill test using a gas analyzer). Also, lack of studies in pediatric population studying associations between HRV and CRF should be considered (Oliveira et al., 2017), specifically in overweight/obese children.
To date, the majority of studies have often omitted the dependence of HRV on HR in overweight/obese children and adolescents (Rodríguez-Colón et al., 2011; Baum et al., 2013; Da Silva et al., 2014; Redón et al., 2016). This dependency might have significant influence on HRV since resting HR is different per each individual. In this sense, a correction procedure has been proposed which is based on ratios between the HRV parameters and different powers of their corresponding mean R-R intervals (Sacha et al., 2013b). In fact, the correction of HRV for prevailing HR is a kind of adjustment for HR at the individual level and it provides HRV parameters which are independent on HR. This method allows us to investigate the HR impact on the prognostic power of HRV in various clinical circumstances, e.g., this impact turned out to be positive for the prediction of cardiac death but negative for the prediction of non-cardiac death among patients after myocardial infarction (Sacha et al., 2013a). Moreover, in women, who usually present high HR, the exclusion of the HR impact on HRV improved the HRV prediction capability for all modes of death (Sacha et al., 2014). However, only one study has used the corrected parameters in relation to CRF in young adults and found that HR was a better predictor of CRF than the standard and corrected HRV parameters (Grant et al., 2013). In addition, they found that from all HRV indices, pNN50 was the most strongly associated with CRF, but HR was a better indicator of CRF than all the standard and corrected HRV parameters in young adults (Grant et al., 2013).
Our results in overweight/obese children are in line with the previous study in young adults (Grant et al., 2013). We found that some standard HRV parameters were associated with CRF, yet they lost this association after correction for HR, i.e., none of the corresponding corrected HRV parameter correlated with CRF (Table 2) – thus, HR seems to be a driving force in the association between standard HRV and CRF. The multiple regression analysis also confirmed that HR, but not HRV, is the only independent indicator of CRF. However, comparing with the study by Grant et al., our multiple stepwise regression revealed that only HR was a significant determinant (explained 11.1% of changes on CRF), but Grant et al. showed that both HR and HF were relevant (explained 17 and 3.1% respectively). Beyond the age range and weight status, the only methodological remarkable difference between our study and the previous one was that we measured CRF directly from gas analyzer (i.e., lab test) and Grant et al. (Grant et al., 2013) estimated CRF from the Cooper 12 min test (i.e., field test). Collectively, previous findings in young adults and present findings in overweight/obese children suggest that the associations between HRV and CRF exit mainly due to relationship between HR and CRF.
To our knowledge this is the first study showing that HR, as indicator of cardiac autonomic function, is a stronger predictor of CRF than HRV parameters in time and frequency domain in overweight/obese children. Therefore, it seems that the more complex concept of HRV do not provide additional information compared with the HR measurement in children with weight disturbances. This may have an important practical consequence, i.e., therapeutic interventions in such children should be directed at the reduction of HR which is an easy parameter and provides valuable information on CRF. However, future studies should perform these analyses in other populations of different age ranges (i.e., older adults) and health status (i.e., healthy children, obese adults, older adults with cardiovascular diseases, etc.).
Strengths and Limitations
Several limitations need to be acknowledged: (1) since our study design is cross-sectional, we cannot assume a causal relationship; (2) we did not use a gold standard for the HRV measurement, however, the RS800CX has demonstrated to be valid (Gamelin et al., 2008) and reliable (Vasconcellos et al., 2015) for this measurement; (3) some studies have found HRV parameters to be affected by breathing (Wessel et al., 2009; Sidorenko et al., 2016), but this effect is not clear and some parameters can be less affected by changes in respiration (Hill and Siebenbrock, 2009), so we decided not to control the breathing to avoid disturbing the resting status of the participants; (4) HRV analyses were performed only in time and frequency domain, maybe other complex approaches to analyze HRV (i.e., approximate and sample entropies, detrended fluctuation analyses short term exponent, etc.,) (Sassi et al., 2015) could add additional information in relation to the prediction of CRF. The main strengths of our study were: (1) CRF was objectively measured with a gas analyzer (gold standard), using an adapted protocol recommended to low-fit children and (2) to the best of our knowledge this is the first study confirming that HR is a better predictor of CRF than complex HRV parameters in overweight/obese children.
Conclusion
In conclusion, HR is a better predictor of CRF than HRV parameters, in either time or frequency domains, in overweight/obese children. Further, HRV parameters did not add any significant information to the models in our overweight/obese participants. Thus, the complex concept of HRV might not provide additional information to the prediction of CRF. Clinicians and sport scientists that use cardiac autonomic activity as indicator of cardiovascular health represented by CRF should take into consideration that the “simple” HR metric gives more information about CRF levels of overweight/obese children than any HRV parameter. Furthermore, the HR assessment and its interpretation are markedly easier than the ones of HRV parameters.
Ethics Statement
This study was conducted according to the Committee for Research Involving Human Subjects at the University of Granada (Reference: 848, February 2014).
Author Contributions
JM, J M-G, CC-S, and FO contributed to design the experiment and study design. AP-F, JM, JM-G, PM-G, MR-A, CC-S, IE-C, JS, and FO contributed to acquisition, analysis, or interpretation of data for the manuscript. All authors drafted or critically revised the manuscript for important intellectual content and approved the final version of the manuscript to be published.
Funding
The project was funded by the Spanish Ministry of Economy and Competitiveness (Reference DEP2013-47540 and DEP2016-79512-R). CC-S is supported by a grant from the Spanish Ministry of Economy and Competitiveness (BES-2014-068829). AP-F, JM, and JM-G are supported by the Spanish Ministry of Education, Culture and Sport (FPU 16/02760, FPU15/02645 and FPU14/06837; respectively). IE-C is supported by a grant from the Alicia Koplowitz Foundation. MR-A is supported by the National Operational Programme on Youth Employment. PS-U was supported by a grant from CONICYT/BECAS Chile/72180543. Additional support was obtained from Unit of Excellence on Exercise and Health (UCEES) and EXERNET Research Network on Exercise and Health in Special Populations (DEP2005-00046/ACTI). This study has been partially funded by the University of Granada, Plan Propio de Investigación 2016, Excellence actions: Units of Excellence; Unit of Excellence on Exercise and Health (UCEES), and by the Junta de Andalucía, Consejería de Conocimiento, Investigación y Universidades and European Regional Development Fund (ERDF), ref. SOMM17/6107/UGR.
Conflict of Interest Statement
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.
Acknowledgments
The authors would like to thank all the participants who volunteered for this investigation. This work is part of a Ph.D. thesis conducted in the Biomedicine Doctoral Studies of the University of Granada, Spain.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2019.00510/full#supplementary-material
Footnotes
References
Baum, P., Petroff, D., Classen, J., Kiess, W., and Blüher, S. (2013). Dysfunction of autonomic nervous system in childhood obesity: a cross-sectional study. PLoS One 8:e54546. doi: 10.1371/journal.pone.0054546
Bervoets, L., and Massa, G. (2014). Defining morbid obesity in children based on BMI 40 at age 18 using the extended international (IOTF) cut-offs. Pediatr. Obes. 9, e94–e98. doi: 10.1111/j.2047-6310.2014.00217.x
Billman, G. E. (2013). The LF/HF ratio does not accurately measure cardiac sympatho-vagal balance. Front. Physiol. 4:26. doi: 10.3389/fphys.2013.00026
Blom, G. (1958). Statistical Estimates and Transformed Beta Variables. New York, NY: John Wiley & Sons.
Brunetto, A. F., Roseguini, B. T., Silva, B. M., Hirai, D. M., and Guedes, D. P. (2005). Effects of gender and aerobic fitness on cardiac autonomic responses to head-up tilt in healthy adolescents. Pediatr. Cardiol. 26, 418–424. doi: 10.1007/s00246-004-0808-0
Cadenas-Sánchez, C., Mora-González, J., Migueles, J. H., Martín-Matillas, M., Gómez-Vida, J., Escolano-Margarit, M. V., et al. (2016). An exercise-based randomized controlled trial on brain, cognition, physical health and mental health in overweight/obese children (ActiveBrains project): rationale, design and methods. Contemp. Clin. Trials 47, 315–324. doi: 10.1016/j.cct.2016.02.007
Cole, T. J., and Lobstein, T. (2012). Extended international (IOTF) body mass index cut-offs for thinness, overweight and obesity. Pediatr. Obes. 7, 284–294. doi: 10.1111/j.2047-6310.2012.Y00064.x
Crabtree, N. J., Arabi, A., Bachrach, L. K., Fewtrell, M., El-Hajj Fuleihan, G., Kecskemethy, H. H., et al. (2014). Dual-Energy X-Ray absorptiometry interpretation and reporting in children and adolescents: the revised 2013 ISCD pediatric official positions. J. Clin. Desnsitom. 17, 225–242. doi: 10.1016/j.jocd.2014.01.003
Da Silva, D. F., Bianchini, J. A. A., Antonini, V. D. S., Hermoso, D. A. M., Lopera, C. A., Pagan, B. G. M., et al. (2014). Parasympathetic cardiac activity is associated with cardiorespiratory fitness in overweight and obese adolescents. Pediatr. Cardiol. 35, 684–690. doi: 10.1007/s00246-013-0838-6
Davis, C. L., Pollock, N. K., Waller, J. L., Allison, J. D., Dennis, B. A., Bassali, R., et al. (2012). Exercise dose and diabetes risk in overweight and obese children: a randomized controlled trial. JAMA 308, 1103–1112. doi: 10.1001/2012.jama.10762
Dilaveris, P., and Tousoulis, D. (2018). Big data analysis from sweden confirms that resting heart rate in adolescencents is significantly associated with incident heart failure and all-cause mortality. Int. J. Cardiol. 259, 220–221. doi: 10.1016/j.ijcard.2018.02.067
Gamelin, F. X., Baquet, G., Berthoin, S., and Bosquet, L. (2008). Validity of the polar S810 to measure R-R intervals in children. Int. J. Sports Med. 29, 134–138. doi: 10.1055/s-2007-964995
Garaulet, M., Martinez-Nicolas, A., Ruiz, J. R., Konstabel, K., Labayen, I., González-Gross, M., et al. (2016). Fragmentation of daily rhythms associates with obesity and cardiorespiratory fitness in adolescents: the HELENA study. Clin. Nutr. 36, 1558–1566. doi: 10.1016/j.clnu.2016.09.026
Goldstein, D. S., Bentho, O., Park, M. Y., and Sharabi, Y. (2011). Low-frequency power of heart rate variability is not a measure of cardiac sympathetic tone but may be a measure of modulation of cardiac autonomic outflows by baroreflexes. Exp. Physiol. 96, 1255–1261. doi: 10.1113/expphysiol.2010.056259
Grant, C. C., Murray, C., Janse van Rensburg, D. C., and Fletcher, L. (2013). A comparison between heart rate and heart rate variability as indicators of cardiac health and fitness. Front. Physiol. 4:337. doi: 10.3389/fphys.2013.00337
Gutin, B., Barbeau, P., Litaker, M. S., Ferguson, M., and Owens, S. (2000). Heart rate variability in obese children: relations to total body and visceral adiposity, and changes with physical training and detraining. Obes. Res. 8, 12–19. doi: 10.1038/oby.2000.3
Gutin, B., Howe, C. A., Johnson, M. H., Humphries, M. C., Snieder, H., and Barbeau, P. (2005). Heart rate variability in adolescents: relations to physical activity, fitness, and adiposity. Med. Sci. Sports Exerc. 37, 1856–1863. doi: 10.1249/01.mss.0000175867.98628.27
Hildebrand, M., Van Hees, V. T., Hansen, B. H., and Ekelund, U. (2014). Age group comparability of raw accelerometer output from wrist-and hip-worn monitors. Med. Sci. Sports Exerc. 46, 1816–1824. doi: 10.1249/MSS.0000000000000289
Hill, L. K., and Siebenbrock, A. (2009). Are all measures created equal? Heart rate variability and respiration - biomed 2009. Biomed. Sci. Instrum. 45, 71–76.
Kang, S., Ha, G., and Ko, K. (2017). Journal of exercise science & fitness association between resting heart rate, metabolic syndrome and cardiorespiratory fitness in korean male adults. J. Exerc. Sci. Fit. 15, 27–31. doi: 10.1016/j.jesf.2017.06.001
Kannel, W. B., Kannel, C., Paffenbarger, R. S., and Cupples, L. A. (1987). Heart rate and cardiovascular mortality: the framingham study. Am. Heart J. 113, 1489–1494. doi: 10.1016/0002-8703(87)90666-1
Koizumi, K., Terui, N., and Kollai, M. (1985). Effect of cardiac vagal and sympathetic nerve activity on heart rate in rhythmic fluctuations. J. Auton. Nerv. Syst. 12, 251–259. doi: 10.1016/0165-1838(85)90065-7
Lahiri, M. K., Kannankeril, P. J., and Goldberger, J. J. (2008). Assessment of autonomic function in cardiovascular disease. physiological basis and prognostic implications. J. Am. Coll. Cardiol. 51, 1725–1733. doi: 10.1016/j.jacc.2008.01.038
Malina, R. M., Rogol, A. D., Cumming, S. P., Coelho e Silva, M. J., and Figueiredo, A. J. (2015). Biological maturation of youth athletes: assessment and implications. Br. J. Sports Med. 49, 852–859. doi: 10.1136/bjsports-2015-094623
Michels, N., Clays, E., De Buyzere, M., Huybrechts, I., Marild, S., Vanaelst, B., et al. (2013). Determinants and reference values of short-term heart rate variability in children. Eur. J. Appl. Physiol. 113, 1477–1488. doi: 10.1007/s00421-012-2572-9
Migueles, J. H., Cadenas-Sanchez, C., Ekelund, U., Delisle Nyström, C., Mora-Gonzalez, J., Löf, M., et al. (2017). Accelerometer data collection and processing criteria to assess physical activity and other outcomes: a systematic review and practical considerations. Sport Med. 47, 1821–1845. doi: 10.1007/s40279-017-0716-0
Migueles, J. H., Cadenas-Sanchez, C., Tudor-Locke, C., Löf, M., Esteban-Cornejo, I., Molina-Garcia, P., et al. (2019). Comparability of published cut - points for the assessment of physical activity: implications for data harmonization. Scand. J. Med. Sci. Sports 29, 566–574. doi: 10.1111/sms.13356
Moore, S. A., McKay, H. A., Macdonald, H., Nettlefold, L., Baxter-Jones, A. D. G., Cameron, N., et al. (2015). Enhancing a somatic maturity prediction model. Med. Sci. Sports Exerc. 47, 1755–1764. doi: 10.1249/MSS.0000000000000588
Nagai, N., Matsumoto, T., Kita, H., and Moritani, T. (2003). Autonomic nervous system activity and the state and development of obesity in Japanese school children. Obesity 11, 25–32. doi: 10.1038/oby.2003.6
Niskanen, J. P., Tarvainen, M. P., Ranta-Aho, P. O., and Karjalainen, P. A. (2004). Software for advanced HRV analysis. Comput. Methods Prog. Biomed. 76, 73–81. doi: 10.1016/j.cmpb.2004.03.004
Oliveira, R. S., Barker, A. R., Wilkinson, K. M., Abbott, R. A., and Williams, C. A. (2017). Is cardiac autonomic function associated with cardiorespiratory fitness and physical activity in children and adolescents? A systematic review of cross-sectional studies. Int. J. Cardiol. 236, 113–122. doi: 10.1016/j.ijcard.2017.02.022
Ortega, F. B., Lavie, C. J., and Blair, S. N. (2016). Obesity and cardiovascular disease. Circ. Res. 118, 1752–1770. doi: 10.1161/CIRCRESAHA.115.306883
Ortega, F. B., Ruiz, J. R., Castillo, M. J., and Sjöström, M. (2008). Physical fitness in childhood and adolescence: a powerful marker of health. Int. J. Obes. 32, 1–11. doi: 10.1038/sj.ijo.0803774
Redón, P., Grassi, G., Redon, J., Álvarez-Pitti, J., and Lurbe, E. (2016). Sympathetic neural activity, metabolic parameters and cardiorespiratory fitness in obese youths. J. Hypertens. 35, 571–577. doi: 10.1097/HJH.0000000000001200
Rodríguez-Colón, S. M., Bixler, E. O., Li, X., and Vgontzas, A. N. (2011). Obesity is associated with impaired cardiac autonomic modulation in children. Int. J. Pediatr. Obes. 6, 128–134. doi: 10.1038/jid.2014.371
Ross, R., Blair, S. N., Arena, R., Church, T. S., Després, J. P., Franklin, B. A., et al. (2016). Importance of assessing cardiorespiratory fitness in clinical practice: a case for fitness as a clinical vital sign: a scientific statement from the american heart association. Circulation 134, e653–e699. doi: 10.1161/CIR.0000000000000461
Sacha, J. (2013). Why should one normalize heart rate variability with respect to average heart rate. Front. Physiol. 4:90. doi: 10.1186/1475-925X-10-90
Sacha, J. (2014a). Interaction between heart rate and heart rate variability. Ann. Noninvasive Electrocardiol. 19, 207–216. doi: 10.1111/anec.12148
Sacha, J. (2014b). Interplay between heart rate and its variability: a prognostic game. Front. Physiol. 5:347. doi: 10.3389/fphys.2014.00347
Sacha, J., Barabach, S., Statkiewicz-Barabach, G., Sacha, K., Müller, A., Piskorski, J., et al. (2013a). How to select patients who will not benefit from ICD therapy by using heart rate and its variability? Int. J. Cardiol. 168, 1655–1658. doi: 10.1016/j.ijcard.2013.03.040
Sacha, J., Barabach, S., Statkiewicz-Barabach, G., Sacha, K., Müller, A., Piskorski, J., et al. (2013b). How to strengthen or weaken the HRV dependence on heart rate - description of the method and its perspectives. Int. J. Cardiol. 168, 1660–1663. doi: 10.1016/j.ijcard.2013.03.038
Sacha, J., Barabach, S., Statkiewicz-Barabach, G., Sacha, K., Müller, A., Piskorski, J., et al. (2014). Gender differences in the interaction between heart rate and its variability - how to use it to improve the prognostic power of heart rate variability. Int. J. Cardiol. 171, e42–e45. doi: 10.1016/j.ijcard.2013.11.116
Sacha, J., and Pluta, W. (2008). Alterations of an average heart rate change heart rate variability due to mathematical reasons. Int. J. Cardiol. 128, 444–447. doi: 10.1016/j.ijcard.2007.06.047
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
Sidorenko, L., Kraemer, J. F., and Wessel, N. (2016). Standard heart rate variability spectral analysis: does it purely assess cardiac autonomic function? Europace 18:1085. doi: 10.1093/europace/euw078
Silva, C. C., Bertollo, M., Reichert, F. F., Boullosa, D. A., and Nakamura, F. Y. (2016). Reliability of heart rate variability in children: influence of sex and body position during data collection. Pediatr. Exerc. Sci. 29, 228–236. doi: 10.1123/pes.2016-0085
Suminski, R. R., Robertson, R. J., Goss, F. L., and Olvera, N. (2008). Validation O F the omni scale of perceived exertion in a sample of spanish-skeaking youth from the usa. Percept. Mot. Skills 107, 181–188. doi: 10.2466/PMS.107.5.181-188
Tarvainen, M. P., Niskanen, J. P., Lipponen, J. A., Ranta-aho, P. O., and Karjalainen, P. A. (2014). Kubios HRV - heart rate variability analysis software. Comput. Methods Prog. Biomed. 113, 210–220. doi: 10.1016/j.cmpb.2013.07.024
Task Force (1996). Heart rate variability: standards of measurement,physiological interpretation and clinical use. Task force of the european societyof cardiology and the north american society of pacing and electrophysiology. Eur. Hear. J. 17, 354–381. doi: 10.1093/oxfordjournals.eurheartj.a014868
Thayer, J. F., and Lane, R. D. (2007). The role of vagal function in the risk for cardiovascular disease and mortality. Biol. Psychol. 74, 224–242. doi: 10.1016/j.biopsycho.2005.11.013
Tikkanen, E., Gustafsson, S., and Ingelsson, E. (2017). Fitness, physical activity, and cardiovascular disease: longitudinal and genetic analyses in the UK biobank study. bioRxiv [Preprint]. doi: 10.1101/136192
Triggiani, A. I., Valenzano, A., Ciliberti, M. A. P., Moscatelli, F., Villani, S., Monda, M., et al. (2017). Heart rate variability is reduced in underweight and overweight healthy adult women. Clin. Physiol. Funct. Imaging 37, 162–167. doi: 10.1111/cpf.12281
Trimmel, K., Sacha, J., and Huikuri, H. V. (eds) (2015). Heart Rate Variability: Clinical Applications and Interaction Between HRV and Heart Rate. Lausanne: Frontiers Media, doi: 10.3389/978-2-88919-652-4
Van Hees, V. T., Gorzelniak, L., Dean León, E. C., Eder, M., Pias, M., Taherian, S., et al. (2013). Separating movement and gravity components in an acceleration signal and implications for the assessment of human daily physical activity. PLoS One 8:e61691. doi: 10.1371/journal.pone.0061691
Van Someren, E. J. W., Swaab, D. F., Colenda, C. C., Cohen, W., McCall, W. V., and Rosenquist, P. B. (1999). Bright light therapy: improved sensitivity to its effects on rest- activity rhythms in Alzheimer patients by application of nonparametric methods. Chronobiol. Int. 16, 505–518. doi: 10.3109/07420529908998724
Vasconcellos, F. V. A., Seabra, A., Cunha, F. A., Montenegro, R. A., Bouskela, E., and Farinatti, P. (2015). Heart rate variability assessment with fingertip photoplethysmography and polar RS800cx as compared with electrocardiography in obese adolescents. Blood Press. Monit. 20, 351–360. doi: 10.1097/MBP.0000000000000143
Wessel, N., Riedl, M., and Kurths, J. (2009). Is the normal heart rate “chaotic” due to respiration? CHAOS 19:028508. doi: 10.1063/1.3133128
Keywords: parasympathetic, sympathetic, heart rate variability, treadmill, adiposity, youth
Citation: Plaza-Florido A, Migueles JH, Mora-Gonzalez J, Molina-Garcia P, Rodriguez-Ayllon M, Cadenas-Sanchez C, Esteban-Cornejo I, Solis-Urra P, de Teresa C, Gutiérrez Á, Michels N, Sacha J and Ortega FB (2019) Heart Rate Is a Better Predictor of Cardiorespiratory Fitness Than Heart Rate Variability in Overweight/Obese Children: The ActiveBrains Project. Front. Physiol. 10:510. doi: 10.3389/fphys.2019.00510
Received: 23 October 2018; Accepted: 11 April 2019;
Published: 07 May 2019.
Edited by:
Gary Iwamoto, University of Illinois at Urbana–Champaign, United StatesReviewed by:
Fiorenzo Moscatelli, University of Foggia, ItalyZhaowei Kong, University of Macau, China
J. A. Taylor, Harvard Medical School, United States
Copyright © 2019 Plaza-Florido, Migueles, Mora-Gonzalez, Molina-Garcia, Rodriguez-Ayllon, Cadenas-Sanchez, Esteban-Cornejo, Solis-Urra, de Teresa, Gutiérrez, Michels, Sacha and Ortega. 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: Abel Plaza-Florido, YWJlbGFkcmlhbkB1Z3IuZXM=