Next Article in Journal
Prognostic Value of Choline and Other Metabolites Measured Using 1H-Magnetic Resonance Spectroscopy in Gliomas: A Meta-Analysis and Systemic Review
Next Article in Special Issue
Predictive Equations Overestimate Resting Metabolic Rate in Young Chilean Women with Excess Body Fat
Previous Article in Journal
The Importance of HDL-Cholesterol and Fat-Free Percentage as Protective Markers in Risk Factor Hierarchy for Patients with Metabolic Syndrome
Previous Article in Special Issue
Body Weight and Metabolic Rate Changes in Narcolepsy: Current Knowledge and Future Directions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Associations between Intra-Assessment Resting Metabolic Rate Variability and Health-Related Factors

by
Juan M. A. Alcantara
1,2,3,*,
Francisco J. Osuna-Prieto
1,4 and
Abel Plaza-Florido
1,5
1
PROFITH “PROmoting FITness and Health through Physical Activity” Research Group, Department of Physical and Sports Education, Faculty of Sport Sciences, Sport and Health University Research Institute (iMUDS), University of Granada, 18011 Granada, Spain
2
Department of Health Sciences, Institute for Innovation & Sustainable Food Chain Development, Public University of Navarra, Campus Arrosadía, s/n., 31006 Pamplona, Spain
3
Navarra Institute for Health Research, IdiSNA, 31008 Pamplona, Spain
4
Department of Analytical Chemistry, University of Granada, 18071 Granada, Spain
5
Pediatric Exercise and Genomics Research Center, Department of Pediatrics, School of Medicine, University of California at Irvine, Irvine, 92617 CA, USA
*
Author to whom correspondence should be addressed.
Submission received: 16 November 2022 / Revised: 28 November 2022 / Accepted: 1 December 2022 / Published: 4 December 2022
(This article belongs to the Special Issue Resting Metabolic Rate of Individuals)

Abstract

:
In humans, the variation in resting metabolic rate (RMR) might be associated with health-related factors, as suggested by previous studies. This study explored whether the intra-assessment RMR variability (expressed as a coefficient of variation (CV; %)) is similar in men and women and if it is similarly associated with diverse health-related factors. The RMR of 107 young, and relatively healthy adults, was assessed using indirect calorimetry. Then, the CV for volumes of oxygen consumption (VO2) and carbon dioxide production (VCO2), respiratory exchange ratio (RER), and resting energy expenditure (REE) were computed as indicators of intra-assessment RMR variability. Body composition, cardiorespiratory fitness (peak VO2 uptake), circulating cardiometabolic risk factors, and heart rate and its variability (HR and HRV) were assessed. Men presented higher CVs for VO2, VCO2, and REE (all p ≤ 0.001) compared to women. Furthermore, in men, the intra-assessment RER variability was associated with vagal-related HRV parameters and with mean HR (standardized β = −0.36, −0.38, and 0.41, respectively; all p < 0.04). In contrast, no associations were observed in women. In conclusion, men exhibited higher variability (CVs for VO2, VCO2, and REE) compared to women. The CV for RER could be a potential marker of cardiometabolic risk in young men.

1. Introduction

It is well known that the energy cost of self-maintenance varies within species, within and between days, and between sexes [1,2,3,4]. This self-maintenance component of energy cost is commonly referred to as the basal or resting metabolic rate (RMR). Nevertheless, whether this is a “basal” or a “resting” assessment depends on the methodology followed in the experiment [5,6]. The RMR, which is widely defined as the minimum energy needed for maintaining body homeostasis and normal body functions (organ functions, thermoregulation, etc.), may account for 60–70% of the 24-h energy expenditure of sedentary individuals [7].
The RMR component can be assessed in a relatively easy manner using indirect calorimetry and metabolic carts [8], which are considered the reference tool for its assessment [5,6]. Thus, metabolic carts allow us to assess and study the RMR of individuals. The variation in RMR between species is undoubtedly and mostly explained by differences in body mass, although other factors may also influence these differences [1,2,3,9]. In mammals (e.g., humans), the variation in RMR appears to influence behavioral traits and fitness (e.g., peak volume of oxygen consumption [VO2] uptake), among others [10,11,12,13]. Importantly, the RMR seems to be a repeatable component over time, showing within- and between-day reproducibility of 3–8%—that reproducibility percentage range corresponds to humans’ RMR assessments using a metabolic cart system [9]. Nowadays, in humans, studying the differences between sexes, and, more concretely, the variability in energy expenditure (EE), is a matter of interest. In a study by Halsey et al. [14], they observed a greater men vs. women variability in total EE, activity EE (estimated as: 0.9 × total EE − basal EE, aiming to determine EE cost from physical and/or exercise), and basal EE components. Interestingly, even after comparing men and women of the same age, height, and fat and fat-free mass-characteristics that directly influence EE—men exhibited more variability than women [14].
Concerning RMR assessments using metabolic carts, the variability (expressed as the coefficient of variation (CV) in percentage) of the measured parameters, such as the VO2, the volume of carbon dioxide production (VCO2), and the respiratory exchange ratio (RER), have also been of interest. In fact, the CV for these parameters is widely employed as criteria for determining the “gas exchange stability” and, thus, as a cut-off point for selecting the RMR data [15]. Of note, the CV for RMR itself (i.e., the CV for resting EE) is not commonly used as a cut-off point for this data selection. A previous study by Irving et al. [16], conducted in healthy participants, showed that 12 individuals (16% of the sample) exhibited elevated intra-assessment variability (CV for VO2 and VCO2 > 10%) during the entire 45-min RMR assessment. Intriguingly, they observed that 17% of them presented extreme low and high body mass index (BMI) values (BMIs < 17.5 and > 48, respectively) [16]. Thus, they suggested that an extreme BMI could influence the intra-assessment RMR variability, i.e., the higher the BMI, the higher the CVs for VO2 and VCO2. In addition, Reeves et al. [17] observed that 55% out of a sample of 39 participants did not accomplish the gas exchange stability criteria, or, in other words, they presented an elevated intra-assessment variability. It is important to acknowledge that more than half of the study cohort were cancer patients, which may influence that intra-assessment variability to an unknown extent [17]. In this line, the more ill the participant, the greater the variability in the RMR [17,18,19,20,21,22,23], an observation supported by studies performed in different populations with varying health status (e.g., patients suffering from cancer, hemodialysis, traumatic brain injuries, eating disorders, and mechanically ventilated patients). In addition, critically ill patients reduced their RMR daily variability later during their hospital course and stabilization [24]. However, daily RMR variability does not inherently correlate with lower or higher resting EE values during the assessment (i.e., lower or higher kilocalories per day), and hence, all daily variability cannot be completely explained by confounding factors such as nursing care procedures or surgery [24]. Therefore, it is plausible that the individuals’ intra-assessment variability could be influenced by their gender [14], age [21,25], and health status [16,17,18,19,20,21,22,23]. Unfortunately, whether a relationship exists between the intra-assessment variability (i.e., intra-assessment CV for VO2, VCO2, RER, and resting EE) and health-related markers (e.g., body composition, circulating cardiometabolic risk factors, and cardiac autonomic function) has not been deeply explored. Considering all of these factors together, we hypothesized that men will present a higher intra-assessment RMR variability compared to women and that this intra-assessment RMR variability will be associated with classical health-related factors.
Thus, in the present exploratory study, we aimed to examine: (i) the intra-assessment RMR variability (expressed as CV (in percentage) for VO2, VCO2, RER, and resting EE) exhibited by men and women, separately; and (ii) whether the intra-assessment RMR variability (expressed as CV for VO2, VCO2, RER, and resting EE) is associated with health-related factors such as body composition, cardiorespiratory fitness (i.e., peak VO2 uptake), circulating cardiometabolic risk factors, heart rate (HR), and heart rate variability (HRV) parameters in a cohort of relatively healthy young adults. To the best of the authors’ knowledge, our present work is the first study exploring the relationship between the intra-assessment RMR variability and objectively determined health-related markers.

2. Materials and Methods

2.1. Study Subjects

The present cross-sectional study used pre-intervention (i.e., baseline) data from the ACTIBATE randomized control trial (RCT; ClinicalTrials.gov ID: NCT02365129) study [26,27]. A total of 107 relatively healthy young adults were included in the present study. All subjects provided both oral and written informed consent (see the Institutional Review Board Statement Section presented below for extended information). In brief, the inclusion criteria were: (i) being sedentary; (ii) maintaining a stable body weight (change lower than 3 kg over the last months); (iii) not being on a weight loss program; (iv) presenting a normal electrocardiogram; (v) not suffering from chronic (or acute) illness; (vi) not being a smoker; and (vii) not being pregnant or lactating. Extended and detailed information concerning the ACTIBATE study can be found elsewhere [26,27].

2.2. Resting Indirect Calorimetry Assessments

The VO2 and VCO2 gas exchanges were measured using either a CPX Ultima CardiO2 or a CCM Express metabolic cart (Medical Graphics Corp., St. Paul, MN, USA) during a 30-min period, while subjects were at rest (i.e., laying on bed in the supine position), in the morning (~9 AM), and following a 12 h overnight fast. Subjects were instructed not to perform moderate (24 h) and/or vigorous (48 h) exercise or physical activity before the assessment. Subjects were instructed to come to the lab by public transportation or motorized vehicle to avoid physical activity after they woke up. In addition, an acclimation period of 20–30 min was performed before the RMR assessment, as recommended by current guidelines [15]. During the entire gas exchange measurement, subjects were instructed to stay awake, remain silent, breathe normally, and avoid fidgeting [15].
Both metabolic carts mentioned above require a neoprene face mask equipped with a Directconnect™ low-flow sensor (Medical Graphics Corp., St. Paul, MN, USA), measured the VO2 and VCO2 using the same galvanic fuel cell and non-dispersive infrared analyzers [28], and require exactly the same calibration procedures. Before each RMR assessment, flow (using a 3 L syringe) and gas analyzers (using 2 gas bottles of standard gas concentrations) were calibrated accordingly to the manufacturers’ recommendations and instructions.
After each RMR assessment, the resulting gas exchange data was downloaded using the metabolic carts’ specific software (MGCDiagnostic® Breeze Suite, v. 8.1.0.54 SP7; Medical Graphics Corp., St. Paul, MN, USA) at a sampling frequency of 1 min, as extensively detailed elsewhere [29]. As recommended by current guidelines, the first 5 min of data were discarded [15], and thus, the remaining 25 min of data were processed and used for further calculations detailed below. The RMR (i.e., resting energy expenditure [REE]) was estimated using the equation proposed by Weir [30] and assuming no urinary nitrogen excretion (see Figure 1A). Then, the RER was calculated as the VCO2-to-VO2 ratio (see Figure 1A). In addition, for descriptive purposes, we computed the RMR relative to body weight (i.e., RMR/kg of body weight; RMRBW) and the RMR relative to fat-free mass (RMR/kg of fat-free mass; RMRFFM). Finally, as surrogate parameters of intra-assessment RMR variability, the CV for VO2, VCO2, RER, and RMR (thereinafter CV for REE) was calculated and expressed as a percentage for each subject (see Figure 1B). As an example, the CV for VO2 was computed for each participant as: (VO2 standard deviation/VO2 average) × 100. A summary of this entire process is depicted in Figure 1.

2.3. Anthropometry and Body Composition Assessments

The subject’s body weight and height were determined using a scale and a stadiometer (SECA model 799, Hamburg, Germany), respectively. Then, we computed BMI as body weight (kg)/height (m) squared. In addition, waist circumference was assessed twice, using a plastic tape, and the mean of both assessments was used for analyses. In this regard, we calculated a Z-score for waist circumference and used it to compute two different cardiometabolic risk Z-scores, which are detailed below.
Finally, body composition (fat mass, lean mass, and fat mass percentage) was assessed by a whole-body Dual Energy X-ray Absorptiometry scanner (Discovery Wi, Hologic Inc., Bedford, MA, USA).

2.4. Cardiorespiratory Fitness Assessment

The peak VO2 uptake (i.e., the cardiorespiratory fitness; CRF) was determined by indirect calorimetry (CPX Ultima CardiO2, Medical Graphics Corp., St. Paul, MN, USA; metabolic cart information was detailed above), while subjects elicited a maximum-effort graded exercise protocol [31]. As for the resting assessments, we calibrated both the volume (using a 3 L syringe and a high-flow sensor [Medical Graphics Corp., St. Paul, MN, USA]) and gas analyzers (using 2 gas bottles of standard gas concentrations) prior to every test, following the manufacturers’ recommendations. The CPX Ultima CardiO2 is a breath-by-breath metabolic cart that equipped a galvanic fuel cell for measuring VO2 and a non-dispersive infrared analyzer for measuring VCO2 (resolution of both gas analyzers ± 0.1%). Regarding the exercise protocol, briefly, every 1-min, the slope of the treadmill (H/P/cosmos pulsar; H/P/cosmos sports & medical GmbH, Nussdorf-Traunstein, Germany) increased by 1%, while subjects walked at 5.3 km/h until volitional exhaustion (self-reported by the subject). In addition, during the entire exercise protocol, the heart rate was monitored using an electrocardiogram. We considered that subjects performed their maximum effort when they met the following criteria [32]: (i) the RER was ≥1.1; (ii) the self-reported perceived exertion was >9 using the rating of perceived exertion category-ratio scale (RPE-CR10) [33]; and (iii) the heart rate was ≥90% of the age-predicted maximum heart rate (i.e., 209 − 0.73 × age). Lastly, all subjects were instructed not to consume any stimulant substance (e.g., caffeine) before the test (24 h), not to eat (3–5 h), and perform neither intense/vigorous nor moderate physical activity during the previous 48 h and/or the previous 24 h, respectively. Extended information can be found elsewhere [34].
To compute the peak VO2 uptake, we first downloaded the data at a sample frequency of 5 s. Subsequently, from the entire set of recorded data, we looked for the highest VO2 uptake value and then averaged the peak (i.e., highest) VO2 uptake value and the immediately 5-s values prior and after the peak VO2 uptake value (i.e., we used a 15-s data average to be used as that peak VO2 uptake value) [34]. Of note, as an approach to detect possible artefacts in the assessment (e.g., artefacts that may unmask the real VO2 value), we checked the entire data set from the 2nd to the 10th subsequent highest VO2 uptake values (i.e., the 2nd peak VO2 uptake value, the 3rd, the 4th, etc.) [34]. In addition, we computed the CRF relative to body weight (i.e., peak VO2 uptake/kg of body weight; CRFBW) and the CRF relative to fat-free mass (peak VO2 uptake/kg of fat-free mass; CRFFFM).

2.5. Circulating Cardiometabolic Risk Factors and Blood Pressure Assessments

For determining circulating cardiometabolic risk factors, subjects rested (sat) for at least 10 min and underwent an overnight fast of at least 12 h. Then, blood samples were collected from the antecubital vein in Vacutainer® SST™ II Advance tubes (Becton Dickinson, Plymouth, UK) to obtain serum. After collection, tubes were centrifuged, aliquoted, and stored (−80 °C) for later analyses.
Total cholesterol, high- and low-density lipoprotein cholesterol (HDL-C and LDL-C, respectively), triglycerides, and glucose were determined by spectrophotometry (model AU5800; Beckman Coulter, Brea, CA, USA), while insulin was determined by chemiluminescence immunoassay involving UniCel DxI 800 paramagnetic particles (Beckman Coulter, Brea, CA, USA). Finally, using the determined blood glucose and insulin values, the homeostatic model assessment of insulin resistance index (HOMA index) was computed:
HOMA   index = Insulin   ×   Glucose 22.5  
An automatic monitor (HEM 705 CP; Omron Healthcare Co., Kyoto, Japan) was used to assess blood pressure (BP). Two BP measurements were performed (in the right arm) while subjects were resting, and mean values were used for further analyses.
Finally, we computed a cardiometabolic risk Z-score (hereinafter cardiometabolic risk Z-score 1) using classical metabolic syndrome markers, i.e., glucose, triglycerides, HDL-C, BP, and waist circumference. To compute that cardiometabolic risk Z-score, we computed the individual Z-score for each of these markers as: (value − mean)/standard deviation. Regarding the Z-score for the HDL-C, only for this marker were its values inverted (i.e., multiplied by −1); thus, higher values in the Z-score can be interpreted as higher cardiometabolic risk. After computing all individual Z-scores, the cardiometabolic risk Z-score (i.e., the Z-score including all the individual Z-scores) was computed as detailed below:
Cardiometabolic   risk   Z - score   1 = ( Z - score   glucose   +   Z - score   triglycerides   +   Z - score   HDL     C   +   Z - score   mean   blood   pressure   +   Z - score   waist   circumference ) 5  
Moreover, we computed a second cardiometabolic risk Z-score (hereinafter cardiometabolic risk Z-score 2), adding to the metabolic syndrome markers mentioned above in the cardiometabolic risk Z-score 1 the TC, the LDL-C, the insulin, and the HOMA index Z-score values.

2.6. Heart Rate and Heart Rate Variability Assessment

The heart rhythm was recorded over a 15-min period using a Polar RS800CX (Polar Electro, Kempele, Finland; sample frequency 1000 Hz), early in the morning (between 8 and 9 AM), while the subjects were lying on a bed (awake, in the supine position). Of note, this measurement was performed immediately before the indirect calorimetry assessment (i.e., in the same room and under the same ambient temperature and humidity conditions). In addition, subjects were instructed to stay awake, not move too much (i.e., be motionless), and remain silent while the heart rhythm was recorded.
The heart rhythm data was processed using the Kubios software (free version, v.3.0.0, HRV analysis, University of Eastern Finland) [35]. In brief, we excluded the first 5 min of data and manually selected the 5-min period [36,37]. Subsequently, in the selected 5-min period, we applied the medium Kubios threshold-based artefact correction level following current recommendations [38]. Finally, the R-R interval series were detrended using the smoothness prior method with alpha set at 500 ms and a cubic interpolation at the default rate of 4 Hz [36,38].
From the heart rhythm recorded, we derived the HR (in beats per minute). In addition, and following the Guidelines of The European Society of Cardiology task force and The North American Society of Pacing and Electrophysiology [39], we derived the vagal-related HRV variables from time and frequency domains. In the first domain, we derived: (i) the squared root of the mean of the sum of the squares of successive R-R interval differences (RMSSD); (ii) the standard deviation of all normal R-R intervals (SDNN); and (iii) the percentage of pairs of adjacent R-R intervals differing by more than 50 milliseconds (pNN50). In the frequency-domain, we derived the power of the high frequency band (HF; 0.15–0.4 Hz) using the Fast Fourier transformation algorithm.
Finally, we also computed the vagal-related HRV score as proposed in our previous study [40]. In brief, we calculated an individual Z-score for all the aforementioned time and frequency domain parameters (i.e., the RMSSD, the SDNN, the pNN50, and the HF), and subsequently, we ran the following equation [40]:
HRV   score = ( Z - score   RMSSD   +   Z - score   SDNN   +   Z - score   pNN 50 +   Z - score   HF ) 4  

2.7. Statistical Analyses

The normal distribution of the variables was examined using the Kolmogorov-Smirnov test and the visual inspection of histograms. As part of our descriptive analyses, for all health-related outcomes and for both RMR ratios, we compared men vs. women using non-paired t-tests. Then, for analytical purposes, the variables that presented a skewed distribution were transformed using a natural logarithm (ln). Analyses were performed separately for men and women as suggested by previous literature (e.g., [14]). Although the sample sizes were different (n = 35 and n = 72 for men and women, respectively), we performed an analysis of covariance (ANCOVA) with Bonferroni comparisons to examine sex-adjusted mean differences (age as a confounder factor) on indirect calorimetry parameters (CVs for VO2, VCO2, RER, and REE). Multiple linear regressions analyses, adjusting for age, were conducted to test the associations between the CVs obtained for the assessed indirect calorimetry parameters (i.e., the CVs for VO2, VCO2, RER, and REE) and: (i) anthropometric and body composition parameters; (ii) CRF (expressed as absolute and relative values); (iii) circulating cardiometabolic risk factors and BP; and (iv) HR and HRV derived parameters. In addition, the associations between the abovementioned CVs and both cardiometabolic risk Z-scores (i.e., cardiometabolic risk Z-scores 1 and 2) and the vagal-related HRV score, adjusting for age, were also conducted.
Analyses were performed using the Statistical Package for the Social Sciences v.22.0 (IBM SPSS Statistics, IBM Corporation, Chicago, IL, USA). The significance level was set at p < 0.05. The results are presented as mean ± standard deviation (SD), and as standardized β unless otherwise stated. Graphs were created using the software Graph Pad Prism (GraphPad, v. 8.0.2, San Diego, CA, USA).

3. Results

A total of 35 men (age = 23 ± 2 years old) and 72 women (age = 22 ± 2 years old) were included in the study. Table 1 provides descriptive data for the participants in both groups (i.e., men and women). In brief, mean differences were observed in all anthropometry and body composition parameters except fat mass (expressed in kilograms; Table 1). In addition, mean CRF and CRFBW were higher in men compared to women (Table 1). Concerning circulating cardiometabolic risk factors and blood pressure parameters, men showed lower mean HDL-C, higher mean systolic BP, and higher cardiometabolic risk Z-score values (Table 1). Finally, no mean differences were observed for HR and HRV-related outcomes (Table 1). Regarding RMR ratios, we observed no mean differences in RMRBW between men and women (21.2 ± 4.5 kcal/kg/day vs. 22.1 ± 3.4 kcal/kg/day; p = 0.299), but mean differences were observed for the RMRFFM ratio (30.1 ± 5.0 kcal/kg/day vs. 36.4 ± 4.7 kcal/kg/day for men and women, respectively; p < 0.001).
The intra-assessment variability, the individual CVs for VO2, VCO2, RER, and REE for men and women, are depicted in Figure 2. In general terms, men had a higher intra-assessment variability for all parameters. Men presented higher CVs for VO2 (adjusted mean difference of 5.9 ± 1.5%; p < 0.001; Figure 2A), for VCO2 (adjusted mean difference of 4.8 ± 1.3%; p < 0.001; Figure 2B), and for REE (adjusted mean difference of 3.4 ± 1.0%; p = 0.001; Figure 2D), compared to women. No differences were observed for CV for RER (adjusted mean difference of 0.59 ± 0.71%; p = 0.409; Figure 2C). Of note, some men presented CVs for VO2, VCO2, and REE (only one participant for this parameter) > 35% (Figure 2, red circles).
Table 2 shows the results from multiple linear regression analyses between the CVs obtained for the indirect calorimetry parameters and the anthropometric and body composition parameters, the CRF (expressed as absolute and relative values), circulating cardiometabolic risk factors and BP, HR and HRV derived parameters, both cardiometabolic risk Z-scores, and the vagal-related HRV score (all associations were performed separately for men and women and adjusted for age). For men, we observed a negative association between the CV for RER and vagal-related HRV parameters, specifically the RMSSD and the pNN50 (β = −0.36 and −0.38, respectively; both p < 0.04; Table 2). Moreover, a positive association was observed between the CV for RER and the mean HR (β = 0.41; p < 0.03) for men. Concerning the women group, we observed no associations between the CVs for VO2, VCO2, RER, and REE and the health-related parameters (all p > 0.05; Table 2).

4. Discussion

This study represents a first exploration of the intra-assessment RMR variability (expressed as CV for VO2, VCO2, RER, and REE) exhibited by young and relatively healthy men and women and its relationship with health-related parameters. Our results show a remarkable intra-assessment RMR variability, which was significantly greater in men than in women (Figure 2A,B,D). We also found that in men, the intra-assessment RMR variability for the assessed RER variability (i.e., CV for RER) was negatively associated with RMSSD and pNN50 and positively associated with mean HR (Table 2). However, in women, no associations were observed between the intra-assessment RMR variability and the different parameters included in the study.
Previous literature showed that certain anthropometric characteristics such as body weight and height and body composition (e.g., fat-free mass) are directly related to energy expenditure and RMR [41,42,43]. We should highlight that our groups were different in all anthropometric and body composition parameters, except fat mass. In brief, men were heavier, taller, and presented a higher FFM compared to women (Table 1). However, these anthropometric differences may be explained by the fact that men presented an obese phenotype (based on BMI and waist circumference values; Table 1). After accounting for “body size”, the mean RMR of men was similar to that yielded by women, as suggested by the RMRBW ratio (p = 0.299). However, after considering the FFM (i.e., RMRFFM), we observed that men expended less energy during resting compared to women 30.1 ± 5.0 vs. 36.4 ± 4.7 kcal/kg/day (p < 0.001). Interestingly, this issue should be acknowledged, as low RMR values could predispose to weight (re)gain [44], since the RMR component could account for up to 70% of the 24-h energy expenditure of sedentary individuals [7]. On the other hand, and regarding variation in EE components, a recent study by Halsey et al. [14] observed that after comparing men and women of similar characteristics (i.e., same age, height, and fat and fat-free masses), men presented a larger variation in different EE components, including total EE, activity EE (the authors estimated this component as: 0.9 × total EE − basal EE), and basal EE. Thus, based on their results, it is reasonable to hypothesize that morphometric and body composition parameters would have an impact on the intra-assessment RMR variability (i.e., CV for VO2, VCO2, RER, and REE). However, we did not observe a relationship between intra-assessment RMR variability and morphometric and body composition parameters (Table 2). Notably, this lack of association was consistent in both men and women.
Concerning CRF, no association was observed between the intra-assessment RMR variability and the CRF, regardless of whether it was expressed as absolute or relative (CRFBW or CRFFFM) values (Table 2). As occurred previously for the anthropometric parameters, the same results were observed for men and women. This absence of association could be explained, at least partially, by the fact that our sample was composed of sedentary individuals, thus this issue may have influenced our results to an unknown extent, as regular exercise and/or physical activity directly influence several systems and organs [45]. Halsey et al. [14] suggested that the dimorphism (i.e., greater basal energy expenditure variability) observed in men compared to women could be influenced by adaptations produced by exercise and training. Congruently, it is known that physical activity levels are more variable in men than in women [46], and thus Halsey et al. [14] proposed that these adaptations may partially drive the observed differences in variability between men and women.
Contrary to our expectations, no associations were observed between the intra-assessment RMR variability and all the circulating cardiometabolic risk factors included in our study. In this case, the absence of associations was consistent in men and women, as were the systolic and diastolic BP assessments (Table 2). This could be partially explained by the fact that participants were young and metabolically healthy according to their body composition and cardiometabolic risk factor levels (Table 1). In fact, these results concur with previous literature in which an exercise intervention carried out in relatively healthy individuals did not induce changes in the assessed cardiometabolic risk factors [47,48], as well as with the results observed in the ACTIBATE study [27], an RCT from which the data of the present study arose. Trying to overcome this possible limitation, we computed two different cardiometabolic risk Z-scores (see Materials and Methods Section) to further study the association between intra-assessment RMR variability and circulating cardiometabolic risk factors. However, the results remained unaltered even after computing both cardiometabolic risk Z-scores (Table 2).
Finally, we observed that certain vagal-related HRV-derived parameters, which are widely employed as surrogate markers of cardiometabolic health [40], were inversely associated with CV for RER in men (Table 2). In a previous study, we also observed that HRV was differently associated with circulating cardiometabolic risk factors in men and women [49]. In brief, we computed three different HRV ratios and studied whether they were similarly associated with circulating cardiometabolic risk factors in young men and women [49]. Interestingly, we observed that the HRV ratios were associated with the circulating cardiometabolic risk factors in the group of young women but not in men. Here, we observed a positive association between the intra-assessment RER variability and HR and a negative association with vagal-related HRV parameters such as RMSSD and pNN50 in men (Table 2). Previous literature showed that, while the subjects are in the resting state (e.g., sitting or lying), those individuals exhibiting lower values of vagal-related HRV parameters (e.g., low RMSSD values) present a higher risk of suffering from cardiovascular diseases and an increased morality risk [39,50,51,52]. Moreover, literature have also suggested that these individuals presenting higher HR values in the resting state present a worse health status compared to these individuals presenting lower HR values [34,40,53,54]. This is in line with a recent meta-analysis of prospective studies, in which Aune et al. [55] observed an increased risk of suffering from cardiovascular diseases and all-cause mortality in those subjects exhibiting greater resting HR values. Therefore, our results show that men exhibiting a higher CV for RER presented a higher cardiometabolic risk (as suggested by the directions of the associations mentioned above; Table 2). To the best of our knowledge, this is the first study investigating whether the intra-measurement RMR variability is related to HR and HRV in young men and women, and therefore, we cannot perform comparisons between studies. Nevertheless, this relationship is of interest as it could be related to metabolic flexibility, a parameter of interest due to its possible relationship with cardiometabolic health [56]. In brief, metabolic flexibility is considered the ability to shift from one substrate (e.g., carbohydrates) to another (e.g., fat) based on fuel availability [56], an ability that is widely considered a cardiometabolic health marker. In this regard, during prolonged fasting, the endogenous fat concentration increases, stimulating fat oxidation [57]. Based on the fasting state of our participants (12 h), they should theoretically be predominantly oxidizing fat, and thus, the expected CV for RER should be low. However, we observed a high CV for RER in men (Figure 2), which may indicate that they were alternating from one substrate to another, thereby suggesting a putatively impaired metabolic flexibility in comparison to women. This issue suggests that exhibiting a high CV for RER in resting conditions could be a potential marker of cardiometabolic risk. Of note, about 90% of VO2 is coupled to adenosine triphosphate production by the mitochondria, with nearly 19−28% being used by the sodium–potassium pump [3], and approximately 20−30% being coupled to the basal mitochondrial proton and electron leak [58]. Thus, variation in RMR between individuals, populations, and/or species could be mostly attributed to differences in mitochondrial function [59], although individual differences in humans are mostly explained by a different mitochondrial O2 affinity [60]. Unfortunately, our study design precludes further study of this issue and whether differences in mitochondrial O2 affinity might be mediating our observations. Future studies involving different study populations (e.g., unhealthy vs. healthy participants), study designs (e.g., muscle biopsies and mitochondrial respiration analyses), and larger cohorts are needed to confirm our results.
Growing evidence also highlights the potential effect of the subjects’ dietary habits on the RMR and the energy balance of the individuals [61,62,63]. In fact, diet seems not only to modulate gut microbiota composition but also energy balance [63,64]. In line with this, a recent study that compared the effects of the Mediterranean diet—the dietary pattern that likely resembles the diet followed by our study participants—vs. a vegan diet found that the Mediterranean diet positively shaped salivary microbiota composition (higher abundance of Subflava and Prevotella species) [61]. Specifically, Prevotella abundance was inversely associated with RER, whereas Subflavan abundance was positively associated with RMR [61]. Unfortunately, in their study, the authors did not analyze whether salivary microbiota were further associated with any of the intra-assessment RMR variability parameters of our study [61]. Considering all these findings, future studies evaluating the association between intra-assessment RMR variability and health-related factors should also consider not only the impact of the diet followed by the individuals but also the individuals’ gut and/or salivary microbiota composition as possible modulators of human metabolism and energy balance.
Considering all of this together, and although no associations were observed between the intra-assessment RMR variability and certain health-related factors included in our present study (only associations with HR, RMSSD, and pNN50 were observed), we recommend its inclusion in future studies. In this regard, computing the CV for the gas exchange outcomes is an easy and feasible procedure that does not require additional measurements and, furthermore, can be retrospectively calculated. Thus, considering our results and these advantages, the intra-assessment RMR variability expressed as CV should be considered in the future as a potential marker for evaluating cardiometabolic health. There were a few limitations to our study that deserve attention. The assessment of RMR (thus, intra-assessment RMR variability) was performed using two different metabolic carts equipped with a face mask. Therefore, the use of other metabolic carts or other gas exchange collection systems (e.g., a canopy hood system) may influence the results. The cross-sectional design of our study does not allow us to establish any cause-and-effect relationships, thus, longitudinal studies in which the intra-assessment RMR variability as well as the health-related parameters are measured at different time points within the same participant are needed. Furthermore, it is important to consider that our study was performed in young and relatively healthy adults; thus, studies carried out in other populations (e.g., older populations, ill patients, etc.) and in larger cohorts to increase the statistical power of the analyses are warranted. Finally, we assessed the heart rhythm using a heart rate monitor instead of an electrocardiograph. However, we used the Polar RS800CX heart rate monitor, which has been validated against the gold-standard technique (i.e., electrocardiography); in this regard, the heart rate monitor used in our study is considered a valid instrument for the heart rhythm recording [65,66].

5. Conclusions

Young and relatively healthy men exhibit a remarkable intra-assessment RMR variability in terms of CV for VO2, VCO2, and REE that is greater than that observed for women. Conversely, no differences in CV for RER were observed between men and women. Moreover, in men, the intra-assessment RER variability was negatively associated with RMSSD and pNN50 and positively associated with mean HR. Interestingly, in women, no associations were observed between the intra-assessment RMR variability and the diverse health-related parameters included. Our results suggest that intra-assessment RMR variability, and specifically the CV for RER, could be a potential marker of cardiometabolic risk in young and relatively healthy men. The present study provides novel, preliminary findings suggesting a sexual dimorphism in the association between intra-assessment RMR variability and health-related factors in young, relatively healthy adults.

Author Contributions

Conceptualization, J.M.A.A. and A.P.-F.; methodology, J.M.A.A. and A.P.-F.; formal analysis, J.M.A.A. and A.P.-F.; data curation, J.M.A.A. and A.P.-F.; writing—original draft preparation, J.M.A.A.; writing—review and editing, J.M.A.A., F.J.O.-P. and A.P.-F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Spanish Ministry of Economy and Competitiveness (DEP2016-79512-R and PTA 12264-I). J.M.A.A., F.J.O.P., and A.P.F. are supported by the Spanish Ministry of Education, Culture and Sport (FPU15/04059, FPU16/02828, and FPU 16/02760, respectively). J.M.A.A. is supported by the University of Granada, Plan Propio de Investigación 2020 Programa de Contratos Puente, and by the Grant FJC2020-044453-I funded by MCIN/AEI/10.13039/501100011033 and by “European Union NextGenerationEU/PRTR”. F.J.O.P. is supported by the University of Granada, Plan Propio de Investigación 2022 Programa de Contratos Puente. A.P.F. is supported in part by NIH grant #: U01 TR002004 (REACH project). Additional support was obtained from the Unit of Excellence on Exercise and Health (UCEES) and the EXERNET Research Network on Exercise and Health in Special Populations (DEP2005-00046/ACTI). This study was additionally supported by the Unit of Excellence in Sport and Health (UCEES), granted by the University of Granada and Junta de Andalucía, Consejería de Conocimiento, Investigación y Universidades and European Regional Development Funds (ref. SOMM17/6107/UGR). The funding agencies had no role in study design, data collection and analysis, or manuscript writing.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Human Research Ethics Committee of both the University of Granada (#924) and the Servicio Andaluz de Salud (Centro de Granada, CEI-Granada). In addition, as previously mentioned, the study was registered in ClinicalTrials.gov (ID: NCT02365129).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author, J.M.A.A., upon reasonable request. The data are not publicly available.

Acknowledgments

The authors thank all the participants who took part in the study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Burton, T.; Killen, S.S.; Armstrong, J.D.; Metcalfe, N.B. What causes intraspecific variation in resting metabolic rate and what are its ecological consequences? Proc. R. Soc. B Biol. Sci. 2011, 278, 3465–3473. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Speakman, J.R.; Król, E.; Johnson, M.S. The functional significance of individual variation in basal metabolic rate. Physiol. Biochem. Zool. 2004, 77, 900–915. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. White, C.R.; Kearney, M.R. Determinants of inter-specific variation in basal metabolic rate. J. Comp. Physiol. B 2013, 183, 1–26. [Google Scholar] [CrossRef] [PubMed]
  4. Konarzewski, M.; Książek, A. Determinants of intra-specific variation in basal metabolic rate. J. Comp. Physiol. B Biochem. Syst. Environ. Physiol. 2013, 183, 27–41. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Chen, K.Y.; Smith, S.; Ravussin, E.; Krakoff, J.; Plasqui, G.; Tanaka, S.; Murgatroyd, P.; Brychta, R.; Bock, C.; Carnero, E.; et al. Room Indirect Calorimetry Operating and Reporting Standards (RICORS 1.0): A Guide to Conducting and Reporting Human Whole-Room Calorimeter Studies. Obesity 2020, 28, 1613–1625. [Google Scholar] [CrossRef]
  6. Schoffelen, P.F.M.; Plasqui, G. Classical experiments in whole-body metabolism: Open-circuit respirometry—Diluted flow chamber, hood, or facemask systems. Eur. J. Appl. Physiol. 2018, 118, 33–49. [Google Scholar] [CrossRef] [Green Version]
  7. Lam, Y.Y.; Ravussin, E. Indirect calorimetry: An indispensable tool to understand and predict obesity. Eur. J. Clin. Nutr. 2017, 71, 318–322. [Google Scholar] [CrossRef]
  8. Alcantara, J.M.A.; Galgani, J.E.; Jurado-Fasoli, L.; Dote-Montero, M.; Merchan-Ramirez, E.; Ravussin, E.; Ruiz, J.R.; Sanchez-Delgado, G. Validity of four commercially available metabolic carts for assessing resting metabolic rate and respiratory exchange ratio in non-ventilated humans. Clin. Nutr. 2022, 41, 746–754. [Google Scholar] [CrossRef]
  9. Donahoo, W.T.; Levine, J.A.; Melanson, E.L. Variability in energy expenditure and its components. Curr. Opin. Clin. Nutr. Metab. Care 2004, 7, 599–605. [Google Scholar] [CrossRef]
  10. Nevill, A.; Holder, R.; Markovic, G. Scaling maximum oxygen uptake using lower leg muscle volume provides further insight into the pitfalls of whole body-mass power laws. J. Appl. Physiol. 2006, 101, 1006–1007. [Google Scholar] [CrossRef] [PubMed]
  11. Nevill, A.M.; Markovic, G.; Vucetic, V.; Holder, R. Can greater muscularity in larger individuals resolve the 3/4 power-law controversy when modelling maximum oxygen uptake? Ann. Hum. Biol. 2004, 31, 436–445. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Tolfrey, K.; Barker, A.; Thom, J.M.; Morse, C.I.; Narici, M.V.; Batterham, A.M. Scaling of maximal oxygen uptake by lower leg muscle volume in boys and men. J. Appl. Physiol. 2006, 100, 1851–1856. [Google Scholar] [CrossRef] [Green Version]
  13. Eliakim, A.; Barstow, T.J.; Brasel, J.A.; Ajie, H.; Lee, W.N.P.; Renslo, R.; Berman, N.; Cooper, D.M. Effect of exercise training on energy expenditure, muscle volume, and maximal oxygen uptake in female adolescents. J. Pediatrics 1996, 129, 537–543. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Halsey, L.G.; Careau, V.; Pontzer, H.; Ainslie, P.N.; Andersen, L.F.; Anderson, L.J.; Arab, L.; Baddou, I.; Bedu-Addo, K.; Blaak, E.E.; et al. Variability in energy expenditure is much greater in males than females. J. Hum. Evol. 2022, 171, 103229. [Google Scholar] [CrossRef]
  15. Fullmer, S.; Benson-Davies, S.; Earthman, C.P.; Frankenfield, D.C.; Gradwell, E.; Lee, P.S.P.; Piemonte, T.; Trabulsi, J. Evidence analysis library review of best practices for performing indirect calorimetry in healthy and non-critically ill individuals. J. Acad. Nutr. Diet. 2015, 115, 1417–1446.e2. [Google Scholar] [CrossRef] [PubMed]
  16. Irving, C.J.; Eggett, D.L.; Fullmer, S. Comparing Steady State to Time Interval and Non-Steady State Measurements of Resting Metabolic Rate. Nutr. Clin. Pract. 2017, 32, 77–83. [Google Scholar] [CrossRef]
  17. Reeves, M.M.; Davies, P.S.W.; Bauer, J.; Battistutta, D. Reducing the time period of steady state does not affect the accuracy of energy expenditure measurements by indirect calorimetry. J. Appl. Physiol. 2004, 97, 130–134. [Google Scholar] [CrossRef] [Green Version]
  18. Olejnik, L.A.; Peters, E.N.; Parrott, J.S.; Marcus, A.F.; Brody, R.A.; Hand, R.K.; Fiutem, J.J.; Byham-Gray, L.D. Abbreviated Steady State Intervals for Measuring Resting Energy Expenditure in Patients on Maintenance Hemodialysis. JPEN J. Parenter. Enter. Nutr. 2017, 41, 1348–1355. [Google Scholar] [CrossRef] [Green Version]
  19. McEvoy, C.; Cooke, S.R.; Young, I.S. A reduced abbreviated indirect calorimetry protocol is clinically acceptable for use in spontaneously breathing patients with traumatic brain injury. Nutr. Clin. Pract. 2009, 24, 513–519. [Google Scholar] [CrossRef] [PubMed]
  20. Mathisen, T.F.; Engen, K.M.; Sundgot-Borgen, J.; Stensrud, T. Evaluation of a short protocol for indirect calorimetry in females with eating disorders and healthy controls. Clin. Nutr. ESPEN 2017, 22, 28–35. [Google Scholar] [CrossRef]
  21. Frankenfield, D.C.; Sarson, G.Y.; Blosser, S.A.; Cooney, R.N.; Smith, J.S. Validation of a 5-minute steady state indirect calorimetry protocol for resting energy expenditure in critically ill patients. J. Am. Coll. Nutr. 1996, 15, 397–402. [Google Scholar] [CrossRef]
  22. Popp, C.; Butler, M.; Curran, M.; Illiano, P.; Sevick, M.A.; St-Jules, D. Evaluating steady-state resting energy expenditure using indirect calorimetry in adults with overweight and obesity. Clin. Nutr. 2019, 39, 2220–2226. [Google Scholar] [CrossRef]
  23. McClave, S.A.; Spain, D.A.; Skolnick, J.L.; Lowen, C.C.; Kleber, M.J.; Wickerham, P.S.; Vogt, J.R.; Looney, S.W. Achievement of steady state optimizes results when performing indirect calorimetry. J. Parenter. Enter. Nutr. 2003, 27, 16–20. [Google Scholar] [CrossRef] [PubMed]
  24. Weissman, C.; Kemper, M.; Hyman, A.I. Variation in the resting metabolic rate of mechanically ventilated critically ill patients. Anesth. Analg. 1989, 68, 457–461. [Google Scholar] [CrossRef] [PubMed]
  25. Sanchez-Delgado, G.; Alcantara, J.; Ortiz-Alvarez, L.; Xu, H.; Martinez-Tellez, B.; Labayen, I.; Ruiz, J. Reliability of resting metabolic rate measurements in young adults: Impact of methods for data analysis. Clin. Nutr. 2018, 37, 1618–1624. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Sanchez-Delgado, G.; Martinez-Tellez, B.; Olza, J.; Aguilera, C.M.; Labayen, I.; Ortega, F.B.; Chillon, P.; Fernandez-Reguera, C.; Alcantara, J.M.A.; Martinez-Avila, W.D.; et al. Activating brown adipose tissue through exercise (ACTIBATE) in young adults: Rationale, design and methodology. Contemp. Clin. Trials 2015, 45, 416–425. [Google Scholar] [CrossRef] [PubMed]
  27. Martinez-Tellez, B.; Sanchez-Delgado, G.; Acosta, F.M.; Alcantara, J.M.A.; Amaro-Gahete, F.J.; Martinez-Avila, W.D.; Merchan-Ramirez, E.; Muñoz-Hernandez, V.; Osuna-Prieto, F.J.; Jurado-Fasoli, L.; et al. No evidence of brown adipose tissue activation after 24 weeks of supervised exercise training in young sedentary adults in the ACTIBATE randomized controlled trial. Nat. Commun. 2022, 13, 5259. [Google Scholar] [CrossRef]
  28. Alcantara, J.M.A.; Sanchez-Delgado, G.; Martinez-Tellez, B.; Merchan-Ramirez, E.; Labayen, I.; Ruiz, J.R. Congruent validity and inter-day reliability of two breath by breath metabolic carts to measure resting metabolic rate in young adults. Nutr. Metab. Cardiovasc. Dis. 2018, 28, 929–936. [Google Scholar] [CrossRef]
  29. Alcantara, J.M.A.; Delgado, G.S.; Gahete, F.J.A.; Galgani, J.E.; Ruiz, J.R. Impact of the Method Used to Select Gas Exchange Data for Estimating the Resting Metabolic Rate, as Supplied by Breath-by-Breath Metabolic Carts. Nutrients 2020, 12, 487. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. De Weir, J.B.V. New methods for calculating metabolic rate with special reference to protein metabolism. J. Physiol. 1949, 109, 1–9. [Google Scholar] [CrossRef]
  31. Balke, B.; Ware, R. An experimental study of physical fitness of Air Force personnel. U. S. Armed Forces Med. J. 1959, 10, 675–688. [Google Scholar]
  32. Midgley, A.W.; McNaughton, L.R.; Polman, R.; Marchant, D. Criteria for determination of maximal oxygen uptake: A brief critique and recommendations for future research. Sports Med. 2007, 37, 1019–1028. [Google Scholar] [CrossRef] [PubMed]
  33. Borg, E.; Kaijser, L. A comparison between three rating scales for perceived exertion and two different work tests. Scand. J. Med. Sci. Sports 2006, 16, 57–69. [Google Scholar] [CrossRef]
  34. Plaza-Florido, A.; Amaro-Gahete, F.J.; Acosta, F.M.; Sacha, J.; Alcantara, J.M.A. Heart rate rather than heart rate variability is better associated with cardiorespiratory fitness in adults. Eur. J. Sport Sci. 2022, 22, 836–845. [Google Scholar] [CrossRef]
  35. Tarvainen, M.P.; Niskanen, J.P.; Lipponen, J.A.; Ranta-aho, P.O.; Karjalainen, P.A. Kubios HRV—Heart rate variability analysis software. Comput. Methods Programs Biomed. 2014, 113, 210–220. [Google Scholar] [CrossRef]
  36. Plaza-Florido, A.; Alcantara, J.M.A.; Migueles, J.H.; Amaro-Gahete, F.J.; Acosta, F.M.; Mora-Gonzalez, J.; Sacha, J.; Ortega, F.B. Inter- and intra-researcher reproducibility of heart rate variability parameters in three human cohorts. Sci. Rep. 2020, 10, 11399. [Google Scholar] [CrossRef]
  37. Plaza-Florido, A.; Sacha, J.; Alcantara, J.M.A. Short-term heart rate variability in resting conditions: Methodological considerations. Kardiol. Pol. 2021, 79, 745–755. [Google Scholar] [CrossRef]
  38. Alcantara, J.M.A.; Plaza-florido, A.; Amaro-gahete, F.J.; Acosta, F.M.; Migueles, J.H.; Molina-garcia, P.; Sacha, J.; Sanchez-delgado, G.; Martinez-tellez, B. Impact of Using Different Levels of Threshold-Based Artefact Correction on the Quantification of Heart Rate Variability in Three Independent Human Cohorts. J. Clin. Med. 2020, 9, 325. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Force, T. Heart rate variability. Standards of measurement, physiological interpretation, and clinical use. Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. Eur. Heart J. 1996, 17, 354–381. [Google Scholar]
  40. Plaza-Florido, A.; Alcantara, J.M.A.; Amaro-Gahete, F.J.; Sacha, J.; Ortega, F.B. Cardiovascular Risk Factors and Heart Rate Variability: Impact of the Level of the Threshold-Based Artefact Correction Used to Process the Heart Rate Variability Signal. J. Med. Syst. 2020, 45, 2. [Google Scholar] [CrossRef] [PubMed]
  41. Heymsfield, S.B.; Childers, D.; Beetsch, J.; Allison, D.B.; Pietrobelli, A. Body size and human energy requirements: Reduced mass-specific resting energy expenditure in tall adults. J. Appl. Physiol. 2007, 103, 1543–1550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Pontzer, H.; Yamada, Y.; Sagayama, H.; Ainslie, P.N.; Andersen, L.F.; Anderson, L.J.; Arab, L.; Baddou, I.; Bedu-Addo, K.; Blaak, E.E.; et al. Daily energy expenditure through the human life course. Science 2021, 373, 808–812. [Google Scholar] [CrossRef] [PubMed]
  43. Johnstone, A.M.; Murison, S.D.; Duncan, J.S.; Rance, K.A.; Speakman, J.R. Factors influencing variation in basal metabolic rate include fat-free mass, fat mass, age, and circulating thyroxine but not sex, circulating leptin, or triiodothyronine. Am. J. Clin. Nutr. 2005, 82, 941–948. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Fernández-Verdejo, R.; Aguirre, C.; Galgani, J.E. Issues in Measuring and Interpreting Energy Balance and Its Contribution to Obesity. Curr. Obes. Rep. 2019, 8, 88–97. [Google Scholar] [CrossRef]
  45. Fiuza-Luces, C.; Garatachea, N.; Berger, N.A.; Lucia, A. Exercise is the real polypill. Physiology 2013, 28, 330–358. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Althoff, T.; Sosič, R.; Hicks, J.L.; King, A.C.; Delp, S.L.; Leskovec, J. Large-scale physical activity data reveal worldwide activity inequality. Nature 2017, 547, 336–339. [Google Scholar] [CrossRef]
  47. Sarzynski, M.A.; Burton, J.; Rankinen, T.; Blair, S.N.; Church, T.S.; Després, J.P.; Hagberg, J.M.; Landers-Ramos, R.; Leon, A.S.; Mikus, C.R.; et al. The effects of exercise on the lipoprotein subclass profile: A meta-analysis of 10 interventions. Atherosclerosis 2015, 243, 364–372. [Google Scholar] [CrossRef] [Green Version]
  48. Harrington, D.M.; Champagne, C.M.; Broyles, S.T.; Johnson, W.D.; Tudor-Locke, C.; Katzmarzyk, P.T. Cardiometabolic risk factor response to a lifestyle intervention: A randomized trial. Metab. Syndr. Relat. Disord. 2015, 13, 125–131. [Google Scholar] [CrossRef] [Green Version]
  49. Plaza-Florido, A.; Sacha, J.; Alcantara, J.M.A. Associations between different heart rate variability ratios and cardiometabolic risk factors in young adults. Am. J. Hum. Biol. 2022, e23797, (Online ahead of print). [Google Scholar] [CrossRef]
  50. La Rovere, M.T.; Pinna, G.D.; Maestri, R.; Mortara, A.; Capomolla, S.; Febo, O.; Ferrari, R.; Franchini, M.; Gnemmi, M.; Opasich, C.; et al. Short-term heart rate variability strongly predicts sudden cardiac death in chronic heart failure patients. Circulation 2003, 107, 565–570. [Google Scholar] [CrossRef] [Green Version]
  51. Bigger, J.T.; Fleiss, J.L.; Steinman, R.C.; Rolnitzky, L.M.; Kleiger, R.E.; Rottman, J.N. Frequency domain measures of heart period variability and mortality after myocardial infarction. Circulation 1992, 85, 164–171. [Google Scholar] [CrossRef] [Green Version]
  52. Kleiger, R.E.; Miller, J.P.; Bigger, J.T.; Moss, A.J. Decreased heart rate variability and its association with increased mortality after acute myocardial infarction. Am. J. Cardiol. 1987, 59, 256–262. [Google Scholar] [CrossRef] [PubMed]
  53. Grant, C.C.; Murray, C.; Janse van Rensburg, D.C.; Fletcher, L. A comparison between heart rate and heart rate variability as indicators of cardiac health and fitness. Front. Physiol. 2013, 4, 337. [Google Scholar] [CrossRef] [Green Version]
  54. González-Acedo, A.; Plaza-Florido, A.; Amaro-Gahete, F.J.; Sacha, J.; Alcantara, J.M.A. Associations between heart rate variability and maximal fat oxidation in two different cohorts of healthy sedentary adults. Nutr. Metab. Cardiovasc. Dis. 2022, 32, 2338–2347. [Google Scholar] [CrossRef]
  55. Aune, D.; Sen, A.; ó’Hartaigh, B.; Janszky, I.; Romundstad, P.R.; Tonstad, S.; Vatten, L.J. Resting heart rate and the risk of cardiovascular disease, total cancer, and all-cause mortality—A systematic review and dose–response meta-analysis of prospective studies. Nutr. Metab. Cardiovasc. Dis. 2017, 27, 504–517. [Google Scholar] [CrossRef]
  56. Galgani, J.E.; Fernández-Verdejo, R. Pathophysiological role of metabolic flexibility on metabolic health. Obes. Rev. 2021, 22, e13131. [Google Scholar] [CrossRef] [PubMed]
  57. Samra, J.S.; Clark, M.L.; Humphreys, S.M.; Macdonald, I.A.; Frayn, K.N. Regulation of lipid metabolism in adipose tissue during early starvation. Am. J. Physiol. 1996, 271, E541–E554. [Google Scholar] [CrossRef]
  58. Jastroch, M.; Divakaruni, A.S.; Mookerjee, S.; Treberg, J.R.; Brand, M.D. Mitochondrial proton and electron leaks. Essays Biochem. 2010, 47, 53–67. [Google Scholar] [PubMed] [Green Version]
  59. Speakman, J.R.; Talbot, D.A.; Selman, C.; Snart, S.; McLaren, J.S.; Redman, P.; Krol, E.; Jackson, D.M.; Johnson, M.S.; Brand, M.D. Uncoupled and surviving: Individual mice with high metabolism have greater mitochondrial uncoupling and live longer. Aging Cell 2004, 3, 87–95. [Google Scholar] [CrossRef]
  60. Larsen, F.J.; Schiffer, T.A.; Sahlin, K.; Ekblom, B.; Weitzberg, E.; Lundberg, J.O. Mitochondrial oxygen affinity predicts basal metabolic rate in humans. FASEB J. 2011, 25, 2843–2852. [Google Scholar] [CrossRef]
  61. Daniele, S.; Scarfò, G.; Ceccarelli, L.; Fusi, J.; Zappelli, E.; Biagini, D.; Lomonaco, T.; Di Francesco, F.; Franzoni, F.; Martini, C. The Mediterranean Diet Positively Affects Resting Metabolic Rate and Salivary Microbiota in Human Subjects: A Comparison with the Vegan Regimen. Biology 2021, 10, 1292. [Google Scholar] [CrossRef] [PubMed]
  62. Belstrøm, D. The salivary microbiota in health and disease. J. Oral Microbiol. 2020, 12, 1723975. [Google Scholar] [CrossRef] [PubMed]
  63. Riedl, R.A.; Atkinson, S.N.; Burnett, C.M.L.; Grobe, J.L.; Kirby, J.R. The Gut Microbiome, Energy Homeostasis, and Implications for Hypertension. Curr. Hypertens. Rep. 2017, 19, 27. [Google Scholar] [CrossRef] [Green Version]
  64. Scheithauer, T.P.M.; Dallinga-Thie, G.M.; de Vos, W.M.; Nieuwdorp, M.; van Raalte, D.H. Causality of small and large intestinal microbiota in weight regulation and insulin resistance. Mol. Metab. 2016, 5, 759–770. [Google Scholar] [CrossRef] [PubMed]
  65. Williams, D.W.P.; Jarczok, M.N.; Ellis, R.J.; Hillecke, T.K.; Thayer, J.F.; Koenig, J. Two-week test-retest reliability of the Polar® RS800CXTM to record heart rate variability. Clin. Physiol. Funct. Imaging 2017, 37, 776–781. [Google Scholar] [CrossRef] [PubMed]
  66. Tsitoglou, K.I.; Koutedakis, Y.; Dinas, P.C.; Hernando, A.; Nanomedicina, B.; Khadka, R.; Koirala, B. Validation of the Polar RS800CX for assessing heart rate variability during rest, moderate cycling and post-exercise recovery. F1000Research 2018, 7, 1501. [Google Scholar] [CrossRef]
Figure 1. Representation of the procedure for computing the respiratory exchange ratio (RER) and the resting metabolic rate (RMR; i.e., resting energy expenditure [REE] estimated using the equation proposed by Weir [30] and assuming no urinary nitrogen excretion) using the measured volume of oxygen consumption (VO2) and volume of carbon dioxide production (VCO2) gas exchange (Panel (A)), and the procedure for computing each coefficient of variation (CV, in percentage; Panel B). These are actual data from a participant included in the present study. The blue line represents the VO2, the yellow line represents the VCO2, the black line represents the RER, and the red line represents the RMR (i.e., REE). In both panels, the data highlighted in gray represents the 5 min data discarded following current recommendations [15]. In Panel (B), the dashed line represents the data used to compute the CV for each gas exchange parameter (i.e., CV for VO2, VCO2, RER, and REE). In Panel B, the left y-axis represents the gas exchange in milliliters per minute (mL/min) and RER, and the right y-axis represents the REE in kilocalories per day (kcal/day).
Figure 1. Representation of the procedure for computing the respiratory exchange ratio (RER) and the resting metabolic rate (RMR; i.e., resting energy expenditure [REE] estimated using the equation proposed by Weir [30] and assuming no urinary nitrogen excretion) using the measured volume of oxygen consumption (VO2) and volume of carbon dioxide production (VCO2) gas exchange (Panel (A)), and the procedure for computing each coefficient of variation (CV, in percentage; Panel B). These are actual data from a participant included in the present study. The blue line represents the VO2, the yellow line represents the VCO2, the black line represents the RER, and the red line represents the RMR (i.e., REE). In both panels, the data highlighted in gray represents the 5 min data discarded following current recommendations [15]. In Panel (B), the dashed line represents the data used to compute the CV for each gas exchange parameter (i.e., CV for VO2, VCO2, RER, and REE). In Panel B, the left y-axis represents the gas exchange in milliliters per minute (mL/min) and RER, and the right y-axis represents the REE in kilocalories per day (kcal/day).
Metabolites 12 01218 g001
Figure 2. Column plots for each coefficient of variation (CV, expressed as a percentage) of the data separated by men (red circles) and women (blue circles). Results are presented as mean and standard deviation, and all data points (i.e., individual values). CV were calculated for each individual and parameter, i.e., for the volume of oxygen consumption (VO2; Panel (A)), for the volume of carbon dioxide production (VCO2; Panel (B)), for the respiratory exchange ratio (RER; Panel (C)), and for the resting metabolic rate (RMR; i.e., resting energy expenditure [REE] estimated using the equation proposed by Weir [30] and assuming no urinary nitrogen excretion; Panel (D)). * represents p-values < 0.05 derived from analysis of covariance (ANCOVA) to examine between sex (i.e., men vs. women) adjusted mean differences while taking participants’ age into account as a confounder (i.e., ANCOVA models adjusted for age). N = 35 and 72 for men and women, respectively.
Figure 2. Column plots for each coefficient of variation (CV, expressed as a percentage) of the data separated by men (red circles) and women (blue circles). Results are presented as mean and standard deviation, and all data points (i.e., individual values). CV were calculated for each individual and parameter, i.e., for the volume of oxygen consumption (VO2; Panel (A)), for the volume of carbon dioxide production (VCO2; Panel (B)), for the respiratory exchange ratio (RER; Panel (C)), and for the resting metabolic rate (RMR; i.e., resting energy expenditure [REE] estimated using the equation proposed by Weir [30] and assuming no urinary nitrogen excretion; Panel (D)). * represents p-values < 0.05 derived from analysis of covariance (ANCOVA) to examine between sex (i.e., men vs. women) adjusted mean differences while taking participants’ age into account as a confounder (i.e., ANCOVA models adjusted for age). N = 35 and 72 for men and women, respectively.
Metabolites 12 01218 g002
Table 1. Participants’ descriptive characteristics.
Table 1. Participants’ descriptive characteristics.
Men (n = 35)Women (n = 72)
MeanSDMeanSDp
Anthropometry and body composition parameters
Weight (kg)82166312<0.001
Height (cm)17571647<0.001
BMI (kg/m2)275234<0.001
Fat mass (kg)25112480.558
Fat mass (%)307386<0.001
Fat free mass (kg)557385<0.001
Waist circumference (cm)89147611<0.001
Cardiorespiratory fitness parameters
CRF (ml/min)37457102528433<0.001
CRFBW (ml/[kg/BW]/min)46.28.940.75.90.002
CRFFFM (ml/[kg/FFM]/min)67.79.567.08.10.708
Circulating cardiometabolic risk factors and blood pressure parameters
Glucose (mg/dl)8978650.052
Insulin (UI/ml)96730.257
HOMA index22210.185
Total cholesterol (mg/dl)16132167320.431
HDL-C (mg/dl)4585511<0.001
LDL-C (mg/dl)982896240.707
Triglycerides (mg/dl)884577420.256
Systolic BP (mm Hg)1271211211<0.001
Diastolic BP (mm Hg)70126780.122
Cardiometabolic risk Z-score 10.600.70−0.300.50<0.001
Cardiometabolic risk Z-score 20.400.80−0.200.500.002
Heart rate and heart rate variability parameters
Mean HR (bpm)67116990.537
RMSSD (ms)593462320.694
SDNN (ms)542553230.888
pNN50 (%)332336210.513
HF (ms2)71710.380
HRV Z-score0.011.000.101.000.700
Results are presented as mean and standard deviation (SD), and p-values from students t-test analyses. BMI: body mass index; CRF: cardiorespiratory fitness (i.e., peak VO2 uptake); CRFBW: CRF relative to body weight (peak VO2 uptake/kg of body weight); CRFFFM: CRF relative to fat-free mass (peak VO2 uptake/kg of fat-free mass); HOMA: homeostatic model assessment of insulin resistance index; HDL-C, high-density lipoprotein cholesterol; LDL-C, low-density lipoprotein cholesterol; BP, blood pressure; Cardiometabolic risk Z-scores 1 and 2: the first Z-score included classical metabolic syndrome markers (glucose, triglycerides, HDL-C, BP, and waist circumference), while the second included these classical metabolic syndrome markers and insulin related markers (insulin and HOMA index) and total cholesterol and LDL-C Z-score values; HR, heart rate (in beats per minute [bpm]); RMSSD, the square root of the mean of the sum of the squares of the R-R interval differences; SDNN, standard deviation of normal R-R intervals; pNN50, the percentage of R-R intervals that show a difference higher than 50 milliseconds (ms); HF, the power of the high-frequency band (HF: 0.15–0.4 Hertz); HRV Z-score: heart rate variability Z-score including the RMSSD, the SDNN, the pNN50, and the HF time and frequency domain parameters’ Z-score values.
Table 2. Association between the intra-assessment resting metabolic rate (RMR) variability expressed as a coefficient of variation (CV) for each of the assessed indirect calorimetry parameters, i.e., the CVs for volume of oxygen consumption (VO2), volume of carbon dioxide production (VCO2), respiratory exchange ratio (RER), and resting energy expenditure (REE), and different health-related parameters.
Table 2. Association between the intra-assessment resting metabolic rate (RMR) variability expressed as a coefficient of variation (CV) for each of the assessed indirect calorimetry parameters, i.e., the CVs for volume of oxygen consumption (VO2), volume of carbon dioxide production (VCO2), respiratory exchange ratio (RER), and resting energy expenditure (REE), and different health-related parameters.
CV for VO2CV for VCO2CV for RERCV for REE
MenWomenMenWomenMenWomenMenWomen
βpβpβpβpβpβpβpβp
Anthropometry and body composition parameters
Weight (kg)−0.1000.5580.1030.851−0.1720.3080.0300.800−0.1580.351−0.0660.598−0.1540.362−0.0250.838
Height (cm)0.0700.6950.0670.5800.0670.709−0.0430.7160.0790.656−0.1230.3170.1340.453−0.0450.709
BMI (kg/m2)−0.1430.3930.0540.661−0.2150.1930.0400.740−0.1950.240−0.0090.942−0.2140.197−0.0130.916
Fat mass (kg)−0.1900.2720.1170.336−0.2440.1550.0640.594−0.1090.530−0.0650.603−0.2180.205−0.0220.855
Fat mass (%)−0.2100.2310.0810.508−0.2350.1770.0830.493−0.0610.729−0.0370.769−0.2280.192−0.0210.865
Fat free mass (kg)0.0930.5850.0400.7450.0100.953−0.0330.784−0.1770.291−0.0400.7510.0120.943−0.0330.783
Waist circumference (cm)−0.0860.6130.1260.301−0.1670.3200.0320.794−0.2570.122−0.1610.195−0.2560.125−0.0910.453
Cardiorespiratory fitness parameters
CRF (mL/min)0.0420.821−0.0170.894−0.0540.769−0.0480.698−0.2830.116−0.0060.964−0.1040.576−0.0600.623
CRFBW (mL/[kg/BW]/min)0.0820.650−0.1350.2730.0680.706−0.0850.489−0.1120.5300.0840.5070.0630.728−0.0380.758
CRFFFM (mL/[kg/FFM]/min)0.0050.978−0.0750.543−0.0380.835−0.0280.818−0.1770.3250.0490.696−0.0190.918−0.0480.693
Circulating cardiometabolic risk factors and blood pressure parameters
Glucose (mg/dL)−0.1580.3900.0710.560−0.0370.839−0.0140.9070.2600.153−0.0720.564−0.0250.893−0.0440.713
Insulin (UI/mL)−0.1850.2990.1990.096−0.2530.1500.0980.408−0.0170.925−0.0200.872−0.2770.1230.0810.499
HOMA index −0.1620.3670.2180.072−0.1400.4360.1590.1860.0220.904−0.0200.876−0.1570.3910.1270.289
Total cholesterol (mg/dL)0.0990.5890.0160.8960.1490.4120.1150.330−0.0560.7610.0040.9740.0300.8730.0740.528
HDL-C (mg/dl)0.2060.253−0.1480.2240.1820.313−0.1580.187−0.1190.514−0.0820.5090.1080.559−0.1620.177
LDL-C (mg/dl)0.0860.638−0.0140.9060.1320.4680.0920.435−0.0650.7230.0090.9390.0210.9120.0650.584
Triglycerides (mg/dl)−0.0020.990−0.0570.640−0.0510.781−0.0670.5750.0910.623−0.0430.728−0.0830.659−0.1400.240
Systolic BP (mm Hg)0.0810.6560.0460.701−0.0130.9440.0690.562−0.2310.1960.0490.692−0.0460.8000.0730.542
Diastolic BP (mm Hg)−0.0230.8980.0430.726−0.1480.395−0.0900.454−0.2350.170−0.1330.286−0.2640.123−0.1140.345
Cardiometabolic risk Z-score 1−0.0970.6000.0660.588−0.1520.4090.0460.701−0.0900.6280.0100.934−0.2290.218−0.0460.702
Cardiometabolic risk Z-score 2−0.0620.7400.0620.614−0.1290.4870.0700.562−0.0960.6060.0090.940−0.1950.296−0.0150.900
Heart rate and heart rate variability parameters
Mean HR (bpm)−0.1450.443−0.0390.771−0.0690.7170.0540.6800.4070.0250.0880.5210.0860.6510.0660.616
RMSSD (ms)0.0850.6430.0300.823−0.0150.933−0.0120.928−0.3640.039−0.0710.602−0.0400.828−0.0210.873
SDNN (ms)0.0700.702−0.0460.730−0.0230.900−0.0280.834−0.3300.0620.0710.606−0.0740.686−0.0280.833
pNN50 (%)0.1210.5030.1200.3620.0100.9550.0230.862−0.3760.031−0.0820.550−0.0470.7950.0360.787
HF (ms2)0.0740.6830.0180.891−0.0130.943−0.0760.563−0.2880.101−0.0740.586−0.0330.855−0.0250.848
HRV Z-score0.0470.800−0.0080.9490.0020.989−0.0440.737−0.2310.201−0.0700.611−0.0030.988−0.0260.844
Results are presented as standardized β and p-values from multiple linear regression analyses (adjusted for age). BMI: body mass index; CRF: cardiorespiratory fitness (i.e., peak VO2 uptake); CRFBW: CRF relative to body weight (peak VO2 uptake/kg of body weight); CRFFFM: CRF relative to fat free mass (peak VO2 uptake/kg of fat-free mass); HOMA: homeostatic model assessment of insulin resistance index; HDL-C, high-density lipoprotein cholesterol; LDL-C, low-density lipoprotein cholesterol; BP, blood pressure; Cardiometabolic risk Z-scores 1 and 2: the first Z-score included classical metabolic syndrome markers (glucose, triglycerides, HDL-C, BP, and waist circumference), while the second included these classical metabolic syndrome markers and insulin-related markers (insulin and HOMA index) and total cholesterol and LDL-C Z-score values; HR: heart rate (in beats per minute [bpm]); RMSSD: the square root of the mean of the sum of the squares of the R-R interval differences; SDNN: standard deviation of normal R-R intervals; pNN50: the percentage of R-R intervals that shows a difference higher than 50 milliseconds (ms); HF: the power of the high frequency band (HF: 0.15–0.4 Hertz); HRV Z-score: heart rate variability Z-score including the RMSSD, the SDNN, the pNN50, and the HF time and frequency domains parameters’ Z-score values.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alcantara, J.M.A.; Osuna-Prieto, F.J.; Plaza-Florido, A. Associations between Intra-Assessment Resting Metabolic Rate Variability and Health-Related Factors. Metabolites 2022, 12, 1218. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo12121218

AMA Style

Alcantara JMA, Osuna-Prieto FJ, Plaza-Florido A. Associations between Intra-Assessment Resting Metabolic Rate Variability and Health-Related Factors. Metabolites. 2022; 12(12):1218. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo12121218

Chicago/Turabian Style

Alcantara, Juan M. A., Francisco J. Osuna-Prieto, and Abel Plaza-Florido. 2022. "Associations between Intra-Assessment Resting Metabolic Rate Variability and Health-Related Factors" Metabolites 12, no. 12: 1218. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo12121218

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop