Abstract
Aims Resting heart rate (RHR) is inversely associated with cardiorespiratory fitness (CRF) but few studies have investigated the nature of this relationship in large population samples. We examined the association between RHR and CRF in UK adults and explored factors that may influence this relationship.
Methods and Results In a population-based sample of 5,143 men and 5,722 women (aged 29-65 years), mean (SD) RHR while seated, supine, and during sleep was 67.6 (9.8), 63.5 (8.9), and 56.9 (6.9) bpm, respectively. The age- and sex-adjusted association with CRF as assessed by submaximal treadmill testing was −0.26 (95%CI −0.27; −0.24), −0.31 (95%CI −0.33; −0.29), and −0.31 (95%CI −0.34; −0.29) ml O2 kg-1 beat-1. Sequential adjustment for objectively measured obesity and physical activity attenuated the RHR coefficient by 10% and 50%, respectively. In longitudinal analyses of 6,589 participants re-examined after 6 years, each 1 bpm increase in supine RHR was associated with 0.23 (95%CI 0.20; 0.25) ml O2 min-1 kg-1 decrease in CRF.
Conclusions Across all measures, RHR is inversely associated with CRF; half of this association is explained by obesity and physical activity, suggesting CRF changes achieved through altered behaviour could be tracked through changes in RHR, a notion supported by longitudinal results. As well as its utility as a biomarker of CRF at population-level, serial measurements of RHR may facilitate personal goal setting/evaluation and remote patient monitoring.
Introduction
Cardiorespiratory fitness (CRF) is associated with the incidence of type 2 diabetes, cardiovascular disease, various types of cancer (1–3), as well as all-cause and cause-specific mortality (4,5). In spite of strong evidence of its association with clinical outcomes, it is relatively uncommon for CRF to be measured in clinical care, mostly because its assessment using graded exercise testing is time-consuming, costly and potentially unsafe for some patient groups without medical supervision. Measurement of resting heart rate (RHR) could serve as a viable alternative to clinical CRF testing as it is easy to undertake and is scalable to large populations through the use of wearable sensors. It also has a pattern of association with future health endpoints that is similar to that of CRF (6–12).
The potential clinical utility of using RHR as a proxy for CRF has not been realised in part due to unresolved methodological challenges. RHR values are known to be dependent on the physiological state at the time of measurement (13). In clinical practice, the most common states in which RHR measurement is taken are sitting during blood pressure measurement or lying supine during brief multi-lead ECG measurement. In free-living individuals, however, wearable sensors offer the potential for convenient assessment of heart rate in other states of rest, particularly sleep, and arousal. It is unknown whether differences in RHR measured in these different ways alters its relationship with CRF. It is also unknown to what extent the relationship between RHR and CRF is affected by factors such as adiposity and physical activity, which have established impact on CRF (14,15). Quantifying the degree to which these modifiable factors may explain the relationship between RHR and CRF would strengthen the argument for using RHR as a proxy measure of exercise capacity. Finally, although some studies have described the longitudinal relationship between RHR and CRF (8,16), there is uncertainty about how individual changes in RHR may reflect longer-term changes in CRF in the general population.
Here, we aimed to assess the cross-sectional associations between different measures of RHR and CRF in a large population-based study of UK adults, and how body composition and physical activity may alter the relationship between RHR and CRF. Secondly, we examined whether within-person change in RHR is associated with within-person change in CRF in longitudinal analysis.
Methods
Participants
Participants were recruited by letter from general practice lists in and around Cambridgeshire in the East of England to the Fenland Study, a population-based cohort study designed to investigate interactions between environmental and genetic factors in determining obesity, type 2 diabetes, and related metabolic disorders in young to middle-aged adults (17). Exclusion criteria included prevalent diabetes, pregnancy or lactation, inability to walk unaided, psychosis or terminal illness (life expectancy of ≤ 1 year at the time of recruitment). The study complies with the Declaration of Helsinki and was approved by the Health Research Authority NRES Committee East of England-Cambridge Central. All participants gave written informed consent.
Participants were invited to attend one of three testing facilities (Ely, Wisbech or Cambridge). Recruitment from general practice surgeries began in 2005, collecting data on a total of 12,435 participants (response rate 27%) as previously described (17). For the present analysis, data from 10,865 individuals was included after excluding participants with key missing variables such as RHR, VO2max, physical activity and body composition, as well as participants prescribed beta-blockers (n=315). Compared to the main cohort sample, included participants were 0.4 years younger, 2% more physically active and had a 0.2 bpm lower RHR. They were also 1 cm taller, had 0.5 kg lower body mass and 0.5% less bodyfat but were otherwise similar across other demographic characteristics (Supplementary Table 1).
RHR and CRF were assessed again after a median (interquartile range) of 6 (5-8) years in data that was available at the time of analysis from a subsample of 6,589 participants allowing the examination of the relationship between within-person change in RHR and CRF over time.
Experimental procedures
Participants arrived at the testing facility after an overnight fast to complete baseline clinical assessment and questionnaires. At least half an hour after arrival, resting pulse rate was measured in a seated position while resting blood pressure was being assessed using an Omron 705CP-II (OMRON Healthcare Europe, Hoofddorp, Netherlands) three times, at 1 minute intervals. The mean of these three pulse rate values was used as seated RHR. At least one hour after arrival, RHR was measured with the participant at rest in a supine position using a combined heart rate (HR) and movement sensor (Actiheart, CamNtech, Papworth, UK) attached to the chest at the base of the sternum by two standard ECG electrodes (18,19). HR was recorded for 6 minutes and RHR was calculated as the mean HR measured during the last 3 minutes. Following the visit to the testing facility, participants wore the combined HR and movement sensor continuously for 6 days and nights during free-living, recording at 60-second intervals. The daily minimum HR, determined robustly as the thirtieth lowest minute-by-minute HR reading during each 24-hour period, was averaged across the days to estimate habitual sleeping HR (20).
VO2max was assessed using HR response to a submaximal treadmill test. Participants exercised while treadmill speed and grade were progressively increased across several stages of level walking, inclined walking, and level running. The test was terminated if one of the following criteria were met: 1) the participant wanted to stop, 2) the participant reached 90% of age-predicted maximal HR (208 − 0.7 age) (21); or 3) the participant had exercised above 80% of age-predicted maximal HR for >2 minutes. Predicted workload (physical activity intensity) during the treadmill protocol (20) was regressed against HR (expressed in bpm above sleeping HR) to define the individual’s submaximal relationship between HR and physical activity intensity. The HR-to-activity intensity relationship was then extrapolated to age-predicted maximal HR to predict maximal work capacity. The resulting work capacity was converted to VO2max by adding an estimate of resting metabolic rate (RMR) (22) and dividing by the energetic equivalent of oxygen (23). VO2max was expressed in both ml O2 min-1 kg-1 body mass and as ml O2 min-1 kg-1 fat-free mass. The VO2max estimation procedure has been validated against directly measured VO2max, demonstrating low bias and high correlation (24).
Ethnicity (White, non-White), smoking status (never, former, current), and alcohol intake (units/week) were determined using a self-administered questionnaire. Anthropometric measures were collected by trained personnel. Weight was measured with a calibrated electronic scale (TANITA model BC-418 MA; Tanita, Tokyo, Japan) and height was assessed with a wall-mounted, calibrated stadiometer (SECA 240; Seca, Birmingham, United Kingdom). Body composition was assessed using dual-energy X-ray absorptiometry (DEXA) (GE Lunar Prodigy Advanced fan beam scanner, GE Healthcare, Bedford, United Kingdom) deriving fat, lean and bone mass measurements across all body regions (25).
Physical activity was objectively measured over 6 days using a combined HR and movement sensor (Actiheart, CamNtech, Papworth, UK) with individual calibration of the HR-to-physical activity intensity relationship performed using data from the treadmill test as described above (17). Free-living data was pre-processed (26) and modelled using a branched equation framework (27) to estimate intensity time-series, which was then summarised over time as daily physical activity energy expenditure (PAEE) (kJ/kg per day) (17). Intensity was also expressed in standard metabolic equivalents (METs), using the conversion 1 MET = 71 J/min/kg (∼3.5 ml O2 min-1 kg-1), and summarised as moderate-to-vigorous physical activity for intensity greater than 3.0 METs. In the present analyses, the proportional PAEE accumulated in moderate and vigorous physical activity was used alongside total PAEE to account for both total physical activity and intensity in the same model (17).
For the longitudinal analyses in the subsample, supine RHR was measured again at follow-up using a 15-min rest protocol but the mean of minutes 4 to 6 was used in the present analysis to match the baseline design for derivation of change in supine RHR. CRF was assessed in the same way as baseline.
Statistical analyses
We used sex-stratified regression models to examine cross-sectional associations between RHR and CRF while adjusting for groups of confounding or explanatory factors. The first model (model 1) only adjusted for age, whilst model 2 additionally adjusted for demographic and lifestyle factors (ethnicity, smoking status, alcohol intake). The third model (model 3) additionally adjusted for BMI. The fourth model (model 4) added adjustment for PAEE, and the fifth model (model 5) further adjusted for moderate to vigorous intensity activity (as a fraction of total PAEE). When cardiorespiratory fitness was expressed relative to fat-free mass (FFM) instead of body mass, fat mass index (fat mass divided by height squared) was used as the adjustment for adiposity instead of BMI. Subgroup analyses were performed as follows: Analyses were stratified by groups of age (less than 50y; 50-59y; 60y and greater), BMI (normal weight, BMI 18.5-25; overweight, BMI 25-30; obese, BMI above 30 kg m-2), and PAEE level (<40, 40-60, >60 kJ kg-1 day-1). Descriptive statistics were reported as means and standard deviations, unless specified otherwise. Interrelationships of RHR measures were examined by simple regression. For the longitudinal analysis of the subsample with repeat measures of RHR and CRF, the association between within-person change in RHR and CRF was adjusted for baseline age, sex, RHR, and fitness, as well as age at follow-up. Binscatter plots (5% bins) were used to visualise associations between RHR and VO2max. Statistical analyses were performed with STATA (Version 14.2; StataCorp, College Station, TX); a p-value of 0.05 or less was considered statistically significant for all analyses.
Results
Mean (SD) baseline RHR in the seated position, supine position, and during sleep was 67.6 (9.8), 63.5 (8.9), and 56.9 (6.9) bpm, respectively (Figure 1), and correlations (Pearson r) between these modalities ranged from 0.65 to 0.81 (Supplementary Table 2). On average, RHR was 3 bpm higher and VO2max was 7.7 ml O2 min-1 kg-1 body mass lower in women compared to men. Characteristics of participants included for analyses are summarised in Table 1, stratified by sex and supine RHR categories. Those with higher RHR had higher BMI and body fat levels, lower VO2max values, and were less physically active.
We used several regression models to examine the association between RHR and VO2max per kg body mass (Table 2, Figure 2) and per kg fat-free mass (Supplementary Table 3, supplementary Figure 1). In models only adjusted for age, RHR was significantly associated with VO2max in women and men, irrespective of how RHR was measured. These associations remained almost identical after additional adjustment for ethnicity, smoking, and alcohol use (Model 2). Further adjustment by BMI (Model 3) attenuated the strength of associations by about 10% for VO2max per kg body mass. The equivalent FMI adjustment for RHR associations with VO2max per kg fat-free mass resulted in stronger associations, particularly for sleeping RHR for which beta coefficients were 60% larger in magnitude.
Adjustment for PAEE (Model 4) attenuated the RHR-VO2max relationship by 30-40% irrespective of how VO2max was expressed, with 5-15% additional attenuation by also accounting for the proportion of PAEE expended at moderate and vigorous intensity (Model 5). Associations between RHR and VO2max were relatively similar across all RHR measurement modalities, especially in the maximally adjusted models, but BMI and physical activity attenuated more of the relationship for sleeping HR in women (30% plus another 40%, respectively). We also analysed the association between RHR and VO2max across age, BMI, and PAEE strata separately. In general, associations in these subgroups were similar in strength to pooled associations and all remained statistically significant.
In the subsample of participants (n = 6,589) with 6-yr repeat measures of supine RHR and VO2max, the overall group mean levels of RHR and fitness were relatively similar to baseline values but with a wide range of individual changes over time; the 5th to the 95th percentiles of changes being −11.4 to 9.2 bpm and −10.1 to 10.8 ml O2 min-1 kg-1, respectively. In the longitudinal association analysis, each 1 bpm increase in supine RHR was associated with a 0.23 (95%CI 0.20; 0.25) ml O2 min-1 kg-1 decline in VO2max, adjusted for follow-up time and baseline values of age, sex, RHR, and VO2max. Sex-stratified coefficients were −0.21 (95%CI −0.24; −0.17) and −0.25 (95%CI −0.28; −0.21) ml O2/min/kg per 1 bpm increase in RHR in women and men, respectively (Figure 3).
Discussion
We have documented the strong interrelationships between common measures of resting heart rate (RHR, when seated, lying supine, and during sleep) and investigated their relationship with cardiorespiratory fitness (CRF, estimated maximal oxygen consumption; VO2max) in a large population-based study of UK adults. Cross-sectional analyses showed a consistent inverse association between RHR and VO2max that was apparent across the different RHR measurement approaches and VO2max normalisation conventions (body mass, fat-free mass). A relatively small part of the association between RHR and VO2max was explained by adiposity but a much greater proportion was explained by the more readily modifiable behaviour physical activity. Concordantly in longitudinal analyses, within-person change in RHR was associated with within-person change in VO2max, similar in magnitude to the relationship that was observed cross-sectionally, thus supporting the notion that changes in factors determining CRF is paralleled by changes in RHR measures and hence RHR can be used as a biomarker for CRF.
This is the first study to examine these associations between multiple measures of RHR and CRF and interrogate them together with objective measures of adiposity and physical activity in a large cohort of men and women. We have previously shown that activity levels in this cohort are similar to those reported in national UK surveys (17,28), suggesting that CRF, as the capacity to undertake physical activity, may also be similar to national levels.
Previous detailed studies of the relationship between RHR and CRF have typically been small and focused on identifying electrophysiological predictors of exercise performance (29–32). There are, however, few large-scale studies that have explored this relationship with findings similar to those reported here. The Copenhagen Male Study found an inverse association between VO2max assessed submaximally in 1970 and supine RHR measured by 12-lead ECG about 15 years later in 2798 men (8); the RHR-to-CRF relationship (□ coefficient ≅ −0.19 ml O2 kg-1 beat-1) was shallower than values reported in the present study at around −0.30 ml O2 kg-1 beat-1. In the Danish Health Examination Survey, the relationship between seated RHR and VO2max was less pronounced (□ coefficient ≅ −0.12 ml O2 kg-1 beat-1) when assessed with maximal cycle ergometry in over 10 thousand men and women from 2007-8 (33,34). A weak prospective inverse association between RHR at baseline and VO2max at 23-year follow-up was reported in the Norwegian Nord-Trøndelag Health Study: −0.9 and −0.4 ml O2 kg-1 beat-1 in 807 men and 810 women, respectively (16). In the same study, within-person change in RHR between baseline and follow-up was inversely associated with VO2max at follow-up; within-person change in VO2max was not assessed. Among 56 thousand American patients referred by their physician for exercise testing and excluding those with known heart disease, the age-adjusted coefficient from meta-regression across RHR categories was −0.22 ml O2 kg-1 beat-1 (35), again similar to results reported in our present study, despite the difference in population sampling.
Our study is among the few to examine the influence of factors underpinning the RHR-to-CRF relationship. For all measures of RHR, we report significant inverse associations between RHR and CRF that are independent of age, sex, obesity and physical activity. The age, sex, BMI-, and physical activity-adjusted coefficient for seated RHR was about −0.13 ml O2 kg-1 beat-1. By comparison, a pooled cohort analysis of almost 50 thousand American and British individuals found a similarly adjusted RHR coefficient of about −0.17 ml O2 kg-1 beat-1 (36). However, physical activity was self-reported in those studies which may have inflated the value of the observed coefficients for RHR because of only partial adjustment for physical activity. In parallel, the Tromsø study compared seated RHR and CRF levels in 5017 men and 5607 women when stratified by self-reported activity levels (37), demonstrating significant inverse associations within sex and across physical activity strata. Together, these findings support the notion that RHR and habitual physical activity levels are intrinsically linked to exercise capacity. Thus, changes in CRF achieved through altered physical activity levels could be feasibly monitored with periodic RHR measurements. Additional work is needed to better elucidate the physical activity-CRF relationship, both cross-sectionally and longitudinally, so that the reliability of inferences from such an approach can be clarified.
It is well-recognised that RHR is associated with heart disease (18,38), diabetes (39,40), cancer (11,12), and all-cause, cardiovascular- and cancer-specific mortality (10,41) but the mechanisms underlying these are not fully understood. Knowing that both higher CRF and higher habitual physical activity levels are associated with lower incidence of such diseases and related mortality (5,42), the significant association of RHR with CRF and the degree to which that association is influenced by physical activity and body mass index offers an explanation for the association of RHR with these endpoints. In clinical practice, such as acute hospital settings, the use of RHR is primarily limited to the detection and management of bradycardia, tachycardia and related conditions. Our findings suggest that RHR may offer a relatively low-cost way of assessing an important risk factor for chronic disease among young-to-middle-aged adults and evaluating the effects of interventions targeting fitness in research, clinical practice and remote patient monitoring systems. We have facilitated this by documenting the interrelationships between common measures of RHR in the general population. However, for this to be directly clinically applicable, future work is needed to examine associations between RHR measures and their relation with such diseases in patient groups predisposed to reduced exercise capacity.
Our study has some limitations. We used heart response to a submaximal treadmill test to estimate rather than directly measure VO2max. Even though we have validated this approach (24), associations between RHR and CRF reported here may be influenced by residual error from the VO2max estimation process, which is largely dependent on reaching a percentage of age-predicted maximal HR. Additionally, VO2max estimated from HR response to submaximal exercise is unreliable in those taking medications such as beta-blockers (43). We therefore excluded participants on beta-blockers from analyses, as well as participants not passing the medical screening for treadmill testing and our results are unlikely to generalise to such individuals.
Conclusions
We have addressed several methodological challenges associated with the estimation of CRF from RHR. In a population sample of UK adults we have shown that RHR is inversely associated with CRF across different RHR measurement approaches. Half of this association is explained by modifiable factors such as body size and habitual physical activity, suggesting that changes in RHR may be used to track changes in CRF over time. We further showed that within-person change in RHR was associated with within-person change in CRF. These findings position RHR as a feasible and responsive biomarker of CRF in the general population and in clinical care and research, facilitating personal goal setting, evaluation and remote patient monitoring.
Data Availability
The data analysed during the current study are not publicly available because we have not obtained consent for public data sharing from the study participants.
Funding
The authors were supported by the UK Medical Research Council (MC_UU_12015/1, MC_UU_12015/3, MC_UU_12015/4, MC_UU_12015/5) and the NIHR Cambridge Biomedical Research Centre in Cambridge (IS-BRC-1215-20014).
Conflict of interest
No author declares any competing interest. No support from any additional organisations for the submitted work; no financial relationships with any organisations that might have an interest in the submitted work in the previous three years; no other relationships or activities that could appear to have influenced the submitted work.
Author’s contribution
TG, JYJ, TL, and SB conceptualised the design of the present analysis. NF, SG, and SB are principal investigators and NW is chief investigator of the Fenland Study. KW and SH contributed to data collection and cleaning. TG, TL, and SB analysed the data. TG, JYJ, TL, IP, and SB drafted the manuscript. All authors provided critical feedback and helped shape the research, analysis, and approved the final version of the manuscript. TG, TL, and SB had full access to the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. The corresponding author attests that all listed authors meet authorship criteria.
Acknowledgements
We are grateful to all Fenland Study participants who gave their time and effort. We also thank the functional teams of the MRC Epidemiology Unit at Cambridge (Field Epidemiology, Study Coordination, Data management and IT) for supporting this study.