Article Text

Download PDFPDF

Heritability of resting heart rate and association with mortality in middle-aged and elderly twins
  1. Magnus T Jensen1,
  2. Mette Wod2,3,
  3. Søren Galatius4,
  4. Jacob B Hjelmborg2,
  5. Gorm B Jensen5,6,
  6. Kaare Christensen2,3,7
  1. 1 Department of Cardiology, Rigshospitalet, Copenhagen, Denmark
  2. 2 The Danish Twin Registry, Department of Public Health, University of Southern Denmark, Odense, Denmark
  3. 3 Department of Public Health, The Danish Aging Research Center, University of Southern Denmark, Odense C, Denmark
  4. 4 Department of Cardiology, Copenhagen University Hospital Bispebjerg, Copenhagen, Denmark
  5. 5 The Copenhagen City Heart Study, Copenhagen University Hospital Frederiksberg, Frederiksberg, Denmark
  6. 6 National Institute of Public Health, University of Southern Denmark, Copenhagen, Denmark
  7. 7 Department of Clinical Genetics, Clinical Biochemistry and Pharmacology, Odense University Hospital, Odense, Denmark
  1. Correspondence to Dr Magnus T Jensen, Department of Cardiology, Rigshospitalet, Blegdamsvej 9, 2100 Copenhagen, Denmark; magnustjensen{at}gmail.com

Abstract

Objective Resting heart rate (RHR) possibly has a hereditary component and is associated with longevity. We used the classical biometric twin study design to investigate the heritability of RHR in a population of middle-aged and elderly twins and, furthermore, studied the association between RHR and mortality.

Methods In total, 4282 twins without cardiovascular disease were included from the Danish Twin Registry, hereof 1233 twin pairs and 1816 ‘single twins’ (twins with a non-participating co-twin); mean age 61.7 (SD 11.1) years; 1334 (31.2%) twins died during median 16.3 (IQR 13.8–16.5) years of follow-up assessed through Danish national registers. RHR was assessed by palpating radial pulse.

Results Within pair correlations for RHR adjusted for sex and age were 0.23 (95% CI 0.14 to 0.32) and 0.10 (0.03 to 0.17) for RHR in monozygotic (MZ) and dizygotic (DZ) twin pairs, respectively. Overall, heritability estimates were 0.23 (95% CI 0.15 to 0.30); 0.27 (0.15 to 0.38) for males and 0.17 (0.06 to 0.28) for females. In multivariable models adjusting for age, gender, body mass index, diabetes, hypertension, pulmonary function, smoking, physical activity and zygosity, RHR was significantly associated with mortality (eg, RHR >90 vs 61–70 beats per min: all-cause HR 1.56 (95% CI 1.21 to 2.03); cardiovascular 2.19 (1.30 to 3.67). Intrapair twin comparison revealed that the twin with the higher RHR was significantly more likely to die first and the probability increased with increase in intrapair difference in RHR.

Conclusions RHR is a trait with a genetic influence in middle-aged and elderly twins free of cardiovascular disease. RHR is independently associated with longevity even when familial factors are controlled for in a twin design.

  • Resting heart rate
  • longevity
  • cardiovascular
  • twin study
  • heritability
  • mortality

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Introduction

Resting heart rate (RHR) is a fundamental physiological parameter associated with longevity and may be a heritable trait. Genome-wide association studies (GWAS) have so far found several genetic loci that are significantly associated with RHR. However, when combining the GWAS identified heart rate-increasing alleles, only a few percentages of the variance in heart rate is explained.1 Using the classical twin design, the relative contribution of genetic and environmental factors can be determined. Examination of twins from early infancy to the age of 7 years has shown that RHR in monozygotic (MZ) twins are significantly more correlated than in dizygotic (DZ) twins indicating that variation in RHR in children has a significant genetic component of magnitudes between 0.22 and 0.50.2 In adolescents, twins’ RHR also showed a substantial heritability of 0.68, which was of similar magnitude in twins with a mean age of 44–47 years.3 4 The heritability of RHR in older age groups, however, is not well characterised.

A high heart rate has been shown to be associated with cardiovascular and all-cause mortality in general populations5–7 and in various patient populations from patients with chronic obstructive pulmonary disease,8 cardiovascular disease,9 to patients with diabetes.10 Studies specifically investigating elderly subjects have shown RHR to be associated with short-term mortality in subjects aged 65 years or older followed for up to 6 years,11 but the long-term relationship between RHR and mortality in this age-group is unknown.

In the present study, we investigated the heritability of RHR and its association with cardiovascular and all-cause mortality in 4282 middle-aged and elderly twins without known cardiovascular disease followed for a median of 16 years. The twin design enables us to make intrapair comparisons of RHR and subsequent survival, hereby controlling for familial, genetic and environmental confounding.

Methods

Study population

The individuals in this study were identified among 7023 twins who participated in one of two population-based nationwide cohorts: The Longitudinal Study of Aging Danish Twins (LSADT; third wave) and the Middle Age Danish Twin study (MADT). Both studies were ascertained through the nationwide Danish Twin Registry (DTR).12 LSADT is a cohort sequential study of Danish same-sex twins aged 70 years and older and was initiated in 1995 with follow-up every second year through 2005.13 Data in this study are for the third wave in 1999, as this was the only wave in which RHR was measured. The MADT study included twins born between 1931 and 1952 and were first assessed in 1998.14 The participation rates were 69.9% in the third wave of LSADT (n=2709) and 83.1% in MADT (n=4314). Each survey comprises multidimensional interviews conducted by trained interviewers. All participants provided informed consent and the Danish Scientific Ethical Committees approved both studies which complied with the Helsinki Declaration.

Follow-up and data linkage

Since 1968, all Danish residents have been assigned a unique central registration number (CPR number). Population-based registries in Denmark offer unique possibilities for studies on long-term occurrence of cardiovascular outcomes and mortality as well as discharge records. Exact dates of death were ascertained from the Danish Civil Registration System and causes of death through the Danish Cause of Death Registry. Causes of death have been coded according to the International Classification of Diseases Tenth Edition (ICD-10) codes.

Exclusion criteria

Individual participants with known cardiovascular disease before the survey were excluded. Information about cardiovascular disease and pacemakers was collected from the Danish National Patient Registry. Cardiovascular disease was defined by ICD-8 and ICD-10 codes (400-438 and I20-I99, DZ950, respectively). Self-reported medicine use (beta-blockers, adrenergic stimulators and calcium-channel inhibitors) also leads to exclusion as well as self-reported coronary thrombosis, angina pectoris, arrhythmia, prescription-treated hypertension and ‘general heart trouble’. To exclude the presence of undiagnosed bradycardia or tachyarrhythmia, participants with very low (<40 beats per min (bpm)) or very high (>120 bpm) RHR were excluded.

Resting heart rate

RHR was measured, as recommended,15 by palpating the radial pulse for 30 s and multiplying by 2 to obtain the number of bpm. RHR was measured with the participant in a sitting position after a written questionnaire had been completed by the participant (>100 items in the questionnaire).

Covariates/confounders/mediators in the survival analyses

Diabetes was defined by ICD-8 as 249-250 and by ICD-10 as E10-11 or self-reported during interview. Hypertension was defined as use of medication with Anatomical Therapeutic Chemical Classification (ATC) C02, C03, C08 and C09 obtained from the interview.

Forced expiratory volume in 1 s (FEV1) was measured during the study visit. Self-reported data included medicine use, smoking, physical activity, height (cm) and weight (kg). Smoking was categorised as never, previous or current. The questions on physical activity were not identical in LSADT and MADT, thus physical activity was categorised as light, moderate or high based on self-reported walking, biking and other physical activity.

Outcomes

The main outcome, heritability, was determined as described below. Endpoints for the survival analyses included: (1) death from any cause and (2) cardiovascular death (ICD-10 I20-99).

Statistical analysis

Descriptive analyses

Median RHR, status (alive/dead at end of follow-up) and age at death are tabulated by sex. Unpaired t-test was used to test for differences between genders.

Twin similarity

Intra-class correlations were calculated for RHR and adjusted for gender and age, and heritability analyses were performed using standard biometric methods.16 It is assumed that twins are exposed to shared and unique environmental factors, which are independent of zygosity. Correlations in the genetic component differ between MZ twin pairs (one for both additive and dominant genetic components) and DZ twin pairs (0.5 for additive and 0.25 for dominant). Hence, higher correlation between MZ twin pairs as compared with DZ twin pairs reflects genetic influences. The classical twin model decomposes total phenotypic variance into additive genetic (A), shared family environment (C) dominant genetic (D) and unique environment (E) components.16 We used the free R package mets to fit and compare five models (ACE, ADE, AE, CE and E) to see which model was the best fit of data. The model with the lowest Akaike Information Criterion represents the best balance of goodness of fit and parsimony. Considering the large age span of the participants, we stratified participants into two age groups: 45–65 and >65 years, a cut-off which has previously been used to study RHR in the elderly.11 All models were adjusted for sex and age when fitted.

Survival analyses

Survival analyses using Cox proportional hazard models were performed to study the association between RHR and mortality from all-cause and cardiovascular death, computing crude and adjusted HRs for intrapair comparison using the stratified Cox model with pair-specific baseline hazard functions, thereby exploiting the paired structure of the data which is in analogy to conditional logistic regression in case of complete follow-up.17 Participants were categorised according to RHR (40–50, 51–60, 61–70, 71–80, 81–90, 91–120), with the largest group as reference (61–70 bpm). Analyses included the following covariates: (1) age at interview and gender; (2) age at interview, gender and body mass index (BMI) and (3) age at interview, gender, BMI, diabetes, hypertension, FEV1, smoking, physical activity and zygosity. Follow-up extended from date of interview until 1 June 2015 or death, whatever came first. The proportional hazards assumption was assessed graphically and found to be met. Estimated HR curves of RHR relative to the hazard of RHR with best mortality were performed using a smoothing method based on splines to introduce flexibility into the Cox model18 (see online supplementary data).

Supplementary Material

Survival analysis—intrapair difference in RHR

Finally, among pairs in which at least one member had died during follow-up, we calculated the proportion of pairs in which the co-twin with the higher RHR died first (null hypothesis was 50%). Participants were divided according to intrapair difference in RHR with increasing intrapair difference (≥1 bpm, ≥5 bpm, ≥10 bpm, ≥15 bpm, ≥20 bpm).

Data were analysed using STATA V.13.1 and R V.3.1.3 with R package mets V.1.1.0 and R package smooth HR V.1.0.2.

Results

Descriptive analyses

There were 7023 participants of which 6374 (90.8%) had data on RHR, 31 were excluded due to missing information on zygosity and 2014 due to known heart disease or use of cardiac medication, leaving 4282 individuals for further analyses (males: n=2092; females: n=2190), figure 1. There were 1233 twin pairs and 1816 ‘single twins’ (twins with a non-participating co-twin). Of these 479 were MZ pairs, 431 same-sex DZ pairs (ssDZ) and 323 opposite-sex DZ pairs (osDZ). A total of 1334 (31.2%) of the 4282 participants had died by 1 June 2015. Median RHR was virtually constant over age. Overall, males had a lower RHR than females (70 bpm (IQR 64 to 78) and 72 bpm (IQR 66 to 78), respectively; p<0.001), and females had a higher age of death than males (82.7±10.8 years vs 78.6±10.5 years; p<0.001). Online supplementary figure 1 shows RHR for males and females with 95% confidence bands. An RHR below 60 bpm was found in 16.2% of males and 12.5% of females. Table 1 displays the median RHR, status, age of death of males, females and combined.

Figure 1

Cardiovascular disease pacemaker, etc,: include previous cardiovascular disease, pacemaker, self-reported medicine (beta-blockers, adrenergic stimulators, and calcium channel inhibitors) and self-reported conditions: coronary thrombosis, angina pectoris, arrhythmia, prescription-treated hypertension and ‘general heart trouble’. LSADT: The Longitudinal Study of Aging Danish Twins (third wave); MADT, The Middle Age Danish Twin study; RHR, resting heart rate.

Table 1

Median resting heart rate (median (IQR)), alive/dead at end of follow-up and mean age at death—males, females and both

Twin similarity

Intraclass correlation coefficients and heritability

The intraclass correlation coefficients (ie, within-pair correlation) adjusted for sex and age for RHR were 0.23 (95% CI 0.14 to 0.32) for MZ twin pairs and 0.10 (0.03 to 0.17) in DZ twin pairs (table 2). For ssDZ and osDZ twin pairs, the intraclass correlation coefficients were 0.07 (95% CI 0.00 to 0.16) and 0.14 (95% CI 0.03 to 0.25), respectively. When stratifying into age groups, no significant changes are observed. Table 2 gives the results for the genetic modelling for males, females and both genders. Among both males and females, a combined AE model showed the best fit in terms of AIC. All the models showed moderate genetic influence on RHR with an overall heritability of 0.23 (95% CI 0.15 to 0.30) for both genders, 0.27 (0.15 to 0.38) for males and 0.17 (0.06 to 0.28) for females. Heritability in age groups ranged from 0.17 (0.01 to 0.33; >65 years) to 0.24 (0.15 to 0.32; 45–65 years). The highest heritability estimate was found in males >65 years (0.42 (95% CI 0.14 to 0.69)) and the lowest in females >65 years (0.07 (0 to 0.27)).

Table 2

Intrapair ICC and heritability (from AE model) of RHR

Survival analyses

Among participants, 31.2% died during follow-up and 6.6% died of cardiovascular disease. Median follow-up time was 16.3 (IQR 13.8 to 16.5) years.

Risk increased with increasing RHR (table 3 and online supplementary table 1, stratified by gender). Multivariable adjustments left risk estimates virtually unchanged. A Wald test showed a significant interaction between RHR and gender (p=0.018). The HR curves by RHR confirmed that males and females may differ in RHR optimal for survival (shown in online supplementary figure 2). The curves showed a reference value of 40 bpm optimal for survival for males, adjusted by age at measurement, with increase in risk as RHR increased. For females, the corresponding reference value was 66 bpm and indicated a ‘J-shaped’ risk pattern of increased HRs at lower RHR and increased HRs at higher RHR.

Table 3

Cox regression* adjusted for within-pair difference. HR of mortality according to RHR. Time extends from time of RHR measurement† to death or 1 June 2015

Cox regression stratified by twin pair, that is, controlling for intrapair similarity (see online supplementary table 2), did not change the overall association. When performing Cox regression stratified by twin pair including only MZ twins, a significantly higher mortality was observed in groups with RHR of 81–90 and 90–120 bpm compared with the reference group (see online supplementary table 3).

Intrapair difference in RHR and mortality in twins

Out of the total number of 1233 complete twin pairs, there were 391 pairs in which at least one twin had died during follow-up, and 302 of these pairs were same-sex twins. Fifteen pairs were excluded as they had the exact same RHR, leaving 287 pairs for analysis. The analysis is shown in figure 2, displaying a clear ‘dose–response’ relationship and suggests that an increased intrapair difference in RHR was associated with a higher probability that the co-twin with the higher RHR would die first. In 160 of the 287 (55.7%) same-sex pairs with a difference in RHR of ≥1 bpm, the twin with the higher RHR died first, and in 38 of the 51 (74.5%) same-sex pairs with a difference of ≥20 bpm, the twin with higher RHR died first.

Figure 2

Proportion of pairs in which the co-twin with the higher RHR died first, according to the intrapair difference in RHR with increasing intrapair difference (≥0, ≥5, ≥10, ≥15 and ≥20 bpm). Numbers in brackets show the number of same-sex pairs were the twin with higher RHR died first/number of available pairs. bpm, beats per min; RHR, resting heart rate.

Discussion

In the present study of 4282 middle-aged and elderly twins free from cardiovascular disease at inclusion and followed for a median of 16 years, the main findings were as follows: first, in the middle-aged and elderly persons, RHR is a trait with a significant heritability; second, the higher the RHR the greater the risk of all-cause and cardiovascular mortality; third, as an extension of this result, we found a clear dose–response relationship which showed that the greater the intrapair difference in RHR the more likely that the twin with the higher RHR died first.

RHR as a possible risk factor for adverse events has received increasing attention in recent years and has been investigated in both observational studies5 6 and randomised clinical trials.19 Levine demonstrated that there is a linear, inverse semilogarithmic relation between RHR and life expectancy across the animal kingdom, and that the number of heart beats per lifetime is nearly constant across species.20 It would be reasonable to think that if a certain RHR characterises an individual species, then there must be a certain degree of heritability in RHR, also in humans.

In a large GWAS, combining 46 loci associated with RHR, only 2.5% of the variance of RHR was explained.1 Here, a genetically predicted RHR increase of 5 bpm increased the risk of mortality by 20%. This finding is within the same magnitude as the findings from The Copenhagen City Heart Study,6 where crude analyses showed an increased risk of 15% per 10 bpm. RHR is a modifiable trait influenced by environmental factors such as fitness,7 diabetes,10 chronic obstructive pulmonary disease8 and medication. For instance, in a population of patient with diabetes, 3 months of exercise lowered RHR, in average, by 7 bpm.21 Also, heart rate-modifying drugs play a central role in clinical medicine. However, although RHR is subject to individual differences in health behaviour and life circumstances, a significant degree of the variance of RHR, ranging from weak to moderate, is explained by genetic factors, as shown in the present study, even in the elderly individuals.

Previous studies have found a significant heritability of 0.40–0.66 in twins aged 15 to 20 years,3 22 and a heritability of similar magnitude in subjects around age 40 years.3 23 In the present study, we studied individuals above 45 years with a mean age of 62 years with an equal distribution of males and females. Here, there was a slightly lower, but highly significant, overall heritability of 0.23 for the whole population, 0.27 for males and 0.17 for females with no clear age pattern. Our study demonstrates that genetic influences, which can vary or be stable by age, are present. The findings points towards a stable heritability, which is also seen for a number of ageing-related traits in studies of twins24 in the latter half of the lifespan.

Our study indicates that RHR is a trait with a life-long twin similarity, which is in line with other phenotypes, demonstrating that when twins grow old, they do not grow apart on central ageing phenotypes.24

RHR was significantly associated with both cardiovascular and all-cause mortality. In males, a low heart rate was associated with a decreased risk of adverse events without a lower limit and with a stepwise increase in risk with increase in RHR. For females, our results indicated a ‘J-shaped’ relationship with a trend towards increased risk of mortality at very low RHR. This could be due to clinically undiagnosed conduction delay or block. The present findings are generally in line with other population studies finding an increased risk of mortality with increase in RHR.25 Studies on RHR in the elderly, however, are not all consistent; Perk et al found RHR to be associated with mortality in elderly females but not in males26; in contrast, other studies have found RHR to be associated with mortality only in elderly males27 and elderly females.28 The present findings are possibly due to longer follow-up and greater number of endpoints compared with these studies. We extended these findings by showing that the twin with the higher heart rate was more likely to die first. This intrapair comparison design provides an opportunity to take into account the unobserved familial confounding exploiting that twins share their childhood environment and are matched partly (DZ) or fully (MZ) on genetic makeup. Hence, intrapair comparisons of exposure discordant twin pairs (difference in RHR) will per design be controlling for these familial factors.29 To our knowledge, this methodology has not previously been used in studies of RHR.

Possible study limitations should be addressed. First, in the present study, RHR was assessed by palpating the pulse. While this is in line with recommendations,15 the presence of an undiagnosed arrhythmia could not be detected. Second, there may have been a healthy survivor bias since all subjects survived into their middle age or older and were all without known heart disease. The ‘true’ risk estimate of RHR may therefore have been underestimated. In terms of the heritability estimates, this may have skewed the population towards twins with lower heart rates. Lastly, the present population was primarily of Scandinavian descent, and the results may therefore not be applicable to other populations.

In conclusion, we found that in middle-aged and elderly individuals free of cardiovascular disease, RHR is a trait with a life-long heritability. Also, elevated RHR is associated with increased risk of all-cause and cardiovascular mortality, and the greater the intrapair difference in RHR the more likely the twin with the higher RHR is to die first, demonstrating that the association persists even after controlling for familial factors. RHR is a clinical variable that is very easy to measure. Although there is currently not yet a specific cut-off for what is considered to be a high RHR, we recommend that RHR should be part of the routine assessment in primary prevention, and that individuals with an elevated heart rate should receive particular attention from the primary care physician, initially by way of risk factor assessment and through lifestyle modification.15 30

Key messages

What is already known on this subject?

Resting heart has been shown to be associated with longevity and may have a hereditary determinant.

What might this study add?

In the present study, we show that resting heart rate (RHR) is a trait with a life-long heritability. Also, elevated RHR is associated with increased risk of all-cause and cardiovascular mortality. The greater the intrapair difference in RHR the more likely that the twin with the higher RHR dies first, demonstrating that the association between RHR and mortality persists even after controlling for familial factors.

How might this impact on clinical practice?

RHR is a clinical variable that is very easy to measure, can be obtained with only minimal training and without any advanced equipment. Although there is currently not yet a specific cut-off for what is considered to be a high RHR, we recommend that RHR should be part of the routine assessment in primary prevention, and that individuals with an elevated heart rate should receive particular attention from the primary care physician, initially by way of risk factor assessment and through lifestyle modification.

Acknowledgments

We thank the participants of the Danish Twin Registry.

References

Footnotes

  • Contributors MTJ made substantial contributions to the conception, analysis, interpretation of the data, drafting the work and revising it critically for important intellectual content, performed final approval of the version to be published and agreed to be accountable for all aspects of the work in ensuring that questions related to the accuracy and integrity of any part of the work are appropriately investigated and resolved. MW made substantial contributions to the analysis, interpretation of the data, drafting the work and revising it critically for important intellectual content, performed final approval of the version to be published and agreed to be accountable for all aspects of the work in ensuring that questions related to the accuracy and integrity of any part of the work are appropriately investigated and resolved. SG and JH made substantial contributions to the analysis, interpretation of the data and revising it critically for important intellectual content, performed final approval of the version to be published and agreed to be accountable for all aspects of the work in ensuring that questions related to the accuracy and integrity of any part of the work are appropriately investigated and resolved. GBJ made substantial contributions to the conception, interpretation of the data and revising it critically for important intellectual content, performed final approval of the version to be published and agreed to be accountable for all aspects of the work in ensuring that questions related to the accuracy and integrity of any part of the work are appropriately investigated and resolved. KC made substantial contributions to the conception, acquisition, analysis, interpretation of the data and revising it critically for important intellectual content, performed final approval of the version to be published and agreed to be accountable for all aspects of the work in ensuring that questions related to the accuracy and integrity of any part of the work are appropriately investigated and resolved.

  • Funding The Danish Twin Registry and the Danish Aging Research Center have received grants from the National Program for Research Infrastructure, the Danish Agency for Science Technology and Innovation, Copenhagen, Denmark, 2007 (grant no. 09-063256), the VELUX Foundation, Copenhagen, Denmark and the US National Institutes of Health, Bethesda, Maryland (grant no. P01 AG08761).

  • Competing interests None declared.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data sharing statement Enquiries about statistical code can be directed to the corresponding author.

  • Correction notice This paper has been amended since it was published Online First. Owing to a scripting error, some of the publisher names in the references were replaced with ’BMJ Publishing Group'. This only affected the full text version, not the PDF. We have since corrected these errors and the correct publishers have been inserted into the references.

Linked Articles