| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
(Circulation. 2004;110:2792-2796.)
© 2004 American Heart Association, Inc.
Arrhythmia/Electrophysiology |
From the Vrije Universiteit, Amsterdam (N.H.M.K., G.W., D.d.B., D.P., D.I.B., E.J.C.d.G.), and Trimbos Institute, Utrecht (M.v.d.B.), the Netherlands.
Correspondence to Nina H.M. Kupper, Department of Biological Psychology, Van der Boechorststraat 1, 1081 BT Amsterdam, The Netherlands. E-mail hm.kupper{at}psy.vu.nl
Received December 29, 2003; de novo received April 27, 2004; revision received June 9, 2004; accepted June 10, 2004.
| Abstract |
|---|
|
|
|---|
Methods and Results In 772 healthy twins and singleton siblings, ambulatory ECG was recorded during 24 hours. Two time domain measures of HRV were used: the standard deviations of all normal-to-normal intervals across 5-minute segments (SDNN index) and the root mean square of successive differences between adjacent normal RR intervals (RMSSD). Multivariate genetic analyses across 4 periods of day (morning, afternoon, evening, night) yielded significant estimates for genetic contribution to the mean ambulatory SDNN index (ranging from 35% to 47%) and the mean ambulatory RMSSD (ranging from 40% to 48%).
Conclusions Ambulatory HRV measures are highly heritable traits that can be used to support genetic association and linkage studies in their search for genetic variation influencing cardiovascular disease risk.
Key Words: genetics heart rate nervous system, autonomic
| Introduction |
|---|
|
|
|---|
| Methods |
|---|
|
|
|---|
Of the 1332 offspring who returned a DNA sample (buccal swabs) for the linkage study, 1008 were successfully contacted for a cardiovascular ambulatory monitoring study, of which 192 refused or were excluded. Reasons for exclusion were pregnancy, heart transplantation, pacemaker and known ischemic heart disease, congestive heart failure, or diabetic neuropathy. In 14 of the remaining 816 subjects, no data were available because of equipment failures. Ten subjects showed a very noisy ECG signal and were excluded from the analyses. Subjects (20 in total) using rhythm-altering medication (n=2) and HRV-reducing antidepressants (tricyclic antidepressants [n=2] and benzodiazepines [n=10]) and antihypertensive medication (ß-blockers [n=13]) or a combination of these were excluded from the analysis. The final sample consisted of 218 monozygotic (MZ) twins (79 men), 301 DZ twins (107 men), and 253 singleton siblings (97 men) from 339 families. For the majority of the twin pairs, zygosity was determined by DNA typing; in a small part (8%), zygosity questionnaires were used. The mean age was 31.3 years (SD=10.6) for men and 30.8 years (SD=10.9) for women. The Ethics Committee of Vrije Universiteit approved the study protocol, and all subjects gave written consent before entering the study. No payment was made for participation.
Study Design
Subjects were visited at home on a weekday, before starting their normal daily activities. They were subjected to an interview on health status and current medication use. The Vrije Universiteit Ambulatory Monitoring System 46 (VU-AMS device15,16) was attached, and its operation was explained. Subjects wore the VU-AMS device the entire day and night until awakening the next morning. Every 30 (±10) minutes, the ambulatory device produced an audible alarm beep to prompt them to complete a detailed diary. Subjects wrote a chronological account of activity, posture, location, presence of other persons, and amount of perceived stress during each past 30 minutes. On the following day, the research assistant collected the device at home.
Heart Rate Variability
The VU-AMS device continuously recorded the ECG from a 6-electrode configuration. Two HRV measures were extracted from the interbeat interval time series: the standard deviations of all normal-to-normal intervals (SDNN) and the root mean squares of the successive differences between adjacent normal-to-normal intervals (RMSSD). In addition to cardiac measures, the device also recorded vertical acceleration as a proxy for gross body movement. The vertical accelerometer information was combined with the diary information to divide the entire recording into smaller fragments that were stationary with regard to physical activity and posture, eg, within each fragment no shifts in activity/posture occurred. The fragments were never <5 minutes or >1 hour. They were coded for posture (lying, sitting, standing, walking, and bicycling), activity (eg, desk work, housekeeping, watching television), and location (eg, at home, at work, at a public place). SDNN was computed across all 5-minute periods that fitted in the coded fragment, effectively yielding the SDNN index. SDNN index and RMSSD were averaged over the entire fragment. On the basis of the reported times of dinner and lunch, awakening, and bedtime, mean RMSSD and SDNN index were computed across all fragments in the morning, afternoon, evening, and nighttime sleep periods. In 8% of the subjects, the exact time of dinner, lunch, awakening, or bedtime could not be extracted from either diary or body movement. For these subjects, the missing time was imputed with the use of the mean times of these events in the rest of the sample.
Statistical Analysis
Confounding
The individual differences in ambulatory HRV were expected to be sensitive to 3 main confounders: differences in sex and age17 and differences in physical activity patterns on the measurement day.18 All analyses below were adjusted for age and sex to control for the first 2 of these confounders. Because ambulatory recorded subjects may differ in their activity patterns, the potential influence of physical activity and postural changes on interindividual variance in HRV measures must be taken into account. This was done by calculating the twin correlations twice: once including the entire recording, containing data during all postures, and once including those fragments of the recording during which a subject was either sitting or lying.
Twin Correlations
MZ twins share all their genetic material, whereas DZ twins and siblings share on average 50% of their segregating genes. A larger resemblance of MZ than DZ twins or other first-degree relatives thus indicates that their larger genetic resemblance is associated with a larger phenotypic resemblance.19 To determine the extent to which MZ twin pairs are more similar than DZ or sibling pairs, Pearson correlation coefficients were calculated per zygosity with the use of SPSS-11 (SPSS Inc). All possible MZ and DZ/sib pairs were used.
Structural Equation Modeling
To answer the question of to what extent genes and shared and unshared environment contribute to the variance of SDNN index and RMSSD, biometric genetic models were fitted to the data with the use of the structural equation program Mx.20 First, nested univariate unconstrained models were fitted to test assumptions of the (extended) twin model. For each period of day, we tested the equality of means and variances for MZ twins, DZ twins, and singleton siblings. Likewise, we examined the presence of sex and age effects on the means and variances. In a final step, we tested for heterogeneity of correlations of men versus women and of DZ twins versus singletons.
The resulting most parsimonious unconstrained models were the ones against which the variance decomposition models were tested. The observed variance was decomposed into 3 sources: additive genetic influences (A), shared environment (C), and unshared environment (E), following Neale and Cardon.21 For DZ twins and sibling pairs, similarity in shared environmental influences was fixed at 100%, and similarity of additive genetic influences was fixed at 50%. For MZ twins, similarities of additive genetic and shared environmental influences were fixed at 100%. Unshared environmental influences are uncorrelated in all twin and sibling pairs. After establishing the most parsimonious variance components model (ACE, AE, CE, or E) for each period of day, we used a full 4-variate Cholesky decomposition to test whether the same or different genetic and environmental factors influenced HRV at each of the 4 periods of the day. A priori, we expected a single genetic factor to underlie the variance across all 4 periods for both SDNN index and RMSSD. This was tested by contrasting a full Cholesky decomposition against a genetic factor model, which allows for a common genetic factor and specific additive genetic influences at each period. It was further tested whether unique environmental influences could also be better described by such a factor structure or whether a Cholesky decomposition should be preferred.
Nested models were compared by likelihood ratio test, with the use of twice the difference between the log-likelihoods of 2 models, which is asymptotically distributed as
2. A high
2 against a low gain in degrees of freedom will generate a significant probability value and denotes a worsening of the fit (related to the more parsimonious model).
| Results |
|---|
|
|
|---|
|
Twin and Sibling Correlations
Table 2 shows the resemblance between MZ and DZ/sibling pairs for SDNN index and RMSSD. All correlations were calculated twice: once on all available data and once with the use of fragments in which subjects were either sitting or lying. Despite the potentially large effects of differences in posture and physical activity on HRV similarity in twins and siblings, the correlations based on sitting/lying-only fragments differed only marginally from those that included the entire recording. Although the potential confounding by mixing data across different postures had little actual impact, we decided to restrict further model fitting to the most "pure" data, ie, fragments in which subjects had been sitting or lying (sleep).
|
Structural Equation Modeling
First, we fitted a series of univariate unconstrained models. The means and variances of both SDNN index and RMSSD were equal for MZ and DZ twins and singleton siblings. Importantly, equating male or female correlations or DZ correlations to correlations across any of the other sib/sib pairings (MZ twin/singleton sibling, DZ twin/singleton sibling, singleton sibling/singleton sibling) yielded no significant worsening in the fit of the model. This allowed us to reduce the number of parameters to be estimated but also implies that the results obtained in twins can be generalized to singletons. Both RMSSD and SDNN index decreased with age, in accordance with previous findings,17 and were higher in men at all periods of day. The sex difference for SDNN index and RMSSD repeats previous findings,22 although the opposite has been found for RMSSD.17 In view of their effects, sex and age were retained as covariates in the final variance components analyses.
The resulting most parsimonious unconstrained models were contrasted against different multivariate variance components models (ACE, AE, CE, E). Only additive genetic (A) and unshared environmental (E) sources contributed significantly to individual variation in SDNN index and RMSSD. For the additive genetic variance, the genetic factor model showed the best fit. Unique environmental variance had to be left in full Cholesky decomposition. In this final model, the common genetic factor explained between 28% and 45% of the variance in SDNN index and between 32% and 48% of the variance in RMSSD (Table 3). Specific genetic influences on SDNN index were always present except for the afternoon and added between 2% and 12% to total heritability. Specific genetic influences on RMSSD were present only in the afternoon and during nighttime sleep and added 2% and 8%, respectively, to total heritability.
|
| Discussion |
|---|
|
|
|---|
Our heritability estimates correspond well with those in a previous study in which much shorter ambulatory recording periods (<4 hours) were used. Using segregation analysis, the Kibbutzim family study12 found genetic influences to account for 45% of age- and sex-corrected RMSSD. However, our estimate for SDNN index is substantially larger than the estimate from a family study based on the Framingham Heart Study and the Framingham Offspring Study.11 In this latter study, SDNN was obtained rather than SDNN index. SDNN was averaged over a 2-hour fragment obtained during a routine, scheduled examination at the Framingham Heart Study clinic. Genetic factors accounted for only 19% of the interindividual variation in SDNN. A major difference in the genetic study design might account for these diverging findings. The Framingham studies used spouse and sibling correlations to produce synthetic estimates of variance components. Because a significant spouse correlation was found, the resemblance between siblings was attributed in part to a shared household. Spouse correlation, however, may also reflect assortative mating for exercise behavior, a variable known to be associated with SDNN index.26 It is of note that the age-corrected sibling correlations in the Framingham Study (0.23 to 0.26) correspond very closely to our age-corrected sibling correlations, suggesting that, when the studies can be compared directly, they are actually very consistent.
Our study made use of an extended twin design that strongly increases statistical power to distinguish between components A, C, and E compared with a design including only MZ and DZ twins.27 Although there was sufficient power (at ß=0.80,
=0.05) to detect effects of
23%, no significant common environmental effect was found. The extended twin design further allowed us to test the possibility that results obtained on singleton sibling pairs were identical to those obtained in twin pairs. This is important because the much lower birth weight in twins might be considered to reflect an impaired fetal environment, which according to the "Barker hypothesis" may influence autonomic function.28,29 MZ or DZ twins did not differ from singleton siblings in means, variances, and covariances on any of the measures. The absence of any twin-singleton difference repeats previous findings in other cardiovascular risk factors30 and indicates that our results can be generalized safely to the population at large.
We separated the entire ambulatory recording into 4 periods of day to allow for the possibility that different genetic factors would affect heart rate regulation during awake and sleeping periods or during leisure (evening) and work (morning, afternoon) periods. Although evidence was found for separate genetic factors influencing HRV at different daily periods, their contribution to the total genetic variance was marginal in comparison to the common genetic factor that influenced HRV at all times of day. From a gene-finding point of view, the common genetic factor structure is advantageous on at least 2 accounts. First, using highly genetically correlated multivariate phenotypes can yield higher statistical power to find genes in linkage analysis.31 Second, these genes, by virtue of having a pervasive influence on HRV across all situations, will also have the largest clinical relevance. This assumes that the genes causing low HRV in this relatively young population remain of importance in later life. Note that we cannot exclude expression of different HRV genes throughout the life span.
In conclusion, this study provides a strong confirmation that genes are important in the regulation of ambulatory HRV, a clinically relevant phenotype for a wide range of cardiovascular diseases.3,5,6,24 The next step is to trace the actual genetic polymorphisms that influence ambulatory HRV to provide new angles for preventive therapy. The first whole-genome screens and candidate gene studies have been initiated already,13,32 and these initiatives should be rapidly extended. Because the power to detect genes increases with the availability of genetically correlated repeated measurements, we believe ambulatory HRV to be an important asset in the search for genetic variation influencing cardiovascular disease risk.
| Acknowledgments |
|---|
| References |
|---|
|
|
|---|
2. Nolan J, Batin P, Andrews R, et al. Prospective study of heart rate variability and mortality in chronic heart failure: results of the United Kingdom failure evaluation and assessment of risk trial (UK-Heart). Circulation. 1998; 98: 15101516.
3. Huikuri HV, Tapanainen JM, Lindgren K, et al. Prediction of sudden cardiac death after myocardial infarction in the beta-blocking era. J Am Coll Cardiol. 2003; 42: 652658.
4. Dekker JM, Schouten EG, Klootwijk P, et al. Heart rate variability from short electrocardiographic recordings predicts mortality from all causes in middle-aged and elderly men: the Zutphen Study. Am J Epidemiol. 1997; 145: 899908.
5. Tsuji H, Larson MG, Venditti FJ Jr, et al. Impact of reduced heart rate variability on risk for cardiac events: the Framingham Heart Study. Circulation. 1996; 94: 28502855.
6. Bigger JT, Fleiss JL, Rolnitzky LM, et al. The ability of several short-term measures of RR variability to predict mortality after myocardial infarction. Circulation. 1993; 88: 927934.
7. Schwartz PJ, Vanoli E, Stramba-Badiale M, et al. Autonomic mechanisms and sudden death: new insights from analysis of baroreceptor reflexes in conscious dogs with and without a myocardial infarction. Circulation. 1988; 78: 969979.
8. La Rovere MT, Bersano C, Gnemmi M, et al. Exercise-induced increase in baroreflex sensitivity predicts improved prognosis after myocardial infarction. Circulation. 2002; 106: 945949.
9. Boomsma DI, van Baal GC, Orlebeke JF. Genetic influences on respiratory sinus arrhythmia across different task conditions. Acta Genet Med Gemellol (Roma). 1990; 39: 181191.[Medline] [Order article via Infotrieve]
10. Snieder H, Boomsma DI, Van Doornen LJ, et al. Heritability of respiratory sinus arrhythmia: dependency on task and respiration rate. Psychophysiology. 1997; 34: 317328.[Medline] [Order article via Infotrieve]
11. Singh JP, Larson MG, ODonnell CJ, et al. Heritability of heart rate variability: the Framingham Heart Study. Circulation. 1999; 99: 22512254.
12. Sinnreich R, Friedlander Y, Luria MH, et al. Inheritance of heart rate variability: the Kibbutzim Family Study. Hum Genet. 1999; 105: 654661.[CrossRef][Medline] [Order article via Infotrieve]
13. Busjahn A, Voss A, Knoblauch H, et al. Angiotensin-converting enzyme and angiotensinogen gene polymorphisms and heart rate variability in twins. Am J Cardiol. 1998; 81: 755760.[CrossRef][Medline] [Order article via Infotrieve]
14. Boomsma DI, Beem AL, Dolan CV, et al. Netherlands twin family study of anxious depression (NETSAD). Twin Res. 2000; 3: 323334.[CrossRef][Medline] [Order article via Infotrieve]
15. de Geus EJC, Willemsen AHM, Klaver CHAM, et al. Ambulatory measurement of respiratory sinus arrhythmia and respiration rate. Biol Psychol. 1995; 42: 205227.
16. de Geus EJC, van Doornen LJP. Ambulatory assessment of parasympathic/sympathic balance by impedance cardiography. In: Fahrenberg J, Myrtek M, eds. Ambulatory Assessment: Computer-Assisted Psychological and Psychophysiological Methods in Monitoring and Field Studies. Göttingen, Germany: Hogrefe & Huber Publishers; 1996: 141163.
17. Antelmi I, De Paula RS, Shinzato AR, et al. Influence of age, gender, body mass index, and functional capacity on heart rate variability in a cohort of subjects without heart disease. Am J Cardiol. 2004; 93: 381385.[CrossRef][Medline] [Order article via Infotrieve]
18. Osterhues HH, Hanzel SR, Kochs M, et al. Influence of physical activity on 24-hour measurements of heart rate variability in patients with coronary artery disease. Am J Cardiol. 1997; 80: 14341437.[CrossRef][Medline] [Order article via Infotrieve]
19. Boomsma DI, Busjahn A, Peltonen L. Classical twin studies and beyond. Nat Rev Genet. 2002; 3: 872882.[CrossRef][Medline] [Order article via Infotrieve]
20. Neale MC, Boker SM, Xie G, et al. Mx: Statistical Modeling. 6th ed. Richmond, Va: VCU Department of Psychiatry; 2003.
21. Neale M, Cardon L. Methodology for Genetic Studies of Twins and Families. Dordrecht, Netherlands: Kluwer Academic Publishers; 1992.
22. Umetani K, Singer DH, McCraty R, et al. Twenty-four hour time domain heart rate variability and heart rate: relations to age and gender over nine decades. J Am Coll Cardiol. 1998; 31: 593601.
23. La Rovere MT, Bigger JT Jr, Marcus FI, et al. Baroreflex sensitivity and heart-rate variability in prediction of total cardiac mortality after myocardial infarction. Lancet. 1998; 351: 478484.[CrossRef][Medline] [Order article via Infotrieve]
24. Janszky I, Ericson M, Mittleman MA, et al. Heart rate variability in long-term risk assessment in middle-aged women with coronary heart disease: the Stockholm Female Coronary Risk Study. J Intern Med. 2004; 255: 1321.[CrossRef][Medline] [Order article via Infotrieve]
25. Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. Heart rate variability: standards of measurement, physiological interpretation, and clinical use. Circulation. 1996; 93: 10431065.
26. Pardo Y, Merz CN, Velasquez I, et al. Exercise conditioning and heart rate variability: evidence of a threshold effect. Clin Cardiol. 2000; 23: 615620.[Medline] [Order article via Infotrieve]
27. Posthuma D, Boomsma DI. A note on the statistical power in extended twin designs. Behav Genet. 2000; 30: 147158.[CrossRef][Medline] [Order article via Infotrieve]
28. Phillips DIW, Walker BR, Reynolds RM, et al. Low birth weight predicts elevated plasma cortisol concentrations in adults from 3 populations. Hypertension. 2000; 35: 13011306.
29. Ijzerman RG, Stehouwer CD, De Geus EJ, et al. Low birth weight is associated with increased sympathetic activity: dependence on genetic factors. Circulation. 2003; 108: 566571.
30. de Geus EJC, Posthuma D, Ijzerman RG, et al. Comparing blood pressure of twins and their singleton siblings: being a twin does not affect adult blood pressure. Twin Res. 2001; 4: 385391.[CrossRef][Medline] [Order article via Infotrieve]
31. Allison DB, Thiel B, St Jean P, et al. Multiple phenotype modeling in gene-mapping studies of quantitative traits: power advantages. Am J Hum Genet. 1998; 63: 11901201.[CrossRef][Medline] [Order article via Infotrieve]
32. Singh JP, Larson MG, ODonnell CJ, et al. Genome scan linkage results for heart rate variability (the Framingham Heart Study). Am J Cardiol. 2002; 90: 12901293.[CrossRef][Medline] [Order article via Infotrieve]
This article has been cited by other articles:
![]() |
C. M. M. Licht, E. J. C. de Geus, R. van Dyck, and B. W. J. H. Penninx Association between Anxiety Disorders and Heart Rate Variability in The Netherlands Study of Depression and Anxiety (NESDA) Psychosom Med, June 1, 2009; 71(5): 508 - 518. [Abstract] [Full Text] [PDF] |
||||
![]() |
K.-J. Bar, S. Berger, M. Metzner, M. K. Boettger, S. Schulz, C. T. Ramachandraiah, J. Terhaar, A. Voss, V. K. Yeragani, and H. Sauer Autonomic Dysfunction in Unaffected First-Degree Relatives of Patients Suffering From Schizophrenia Schizophr Bull, April 14, 2009; (2009) sbp024v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. U. Viola, L. M. James, S. N. Archer, and D.-J. Dijk PER3 polymorphism and cardiac autonomic control: effects of sleep debt and circadian phase Am J Physiol Heart Circ Physiol, November 1, 2008; 295(5): H2156 - H2163. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. Vaccarino, R. Lampert, J. D. Bremner, F. Lee, S. Su, C. Maisano, N. V. Murrah, L. Jones, F. Jawed, N. Afzal, et al. Depressive Symptoms and Heart Rate Variability: Evidence for a Shared Genetic Substrate in a Study of Twins Psychosom Med, July 1, 2008; 70(6): 628 - 636. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. K. Lahiri, P. J. Kannankeril, and J. J. Goldberger Assessment of Autonomic Function in Cardiovascular Disease: Physiological Basis and Prognostic Implications J. Am. Coll. Cardiol., May 6, 2008; 51(18): 1725 - 1733. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Kupper, J. Denollet, E. J. C. de Geus, D. I. Boomsma, and G. Willemsen Heritability of Type-D Personality Psychosom Med, September 1, 2007; 69(7): 675 - 681. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. L. T. Uusitalo, E. Vanninen, E. Levalahti, M. C. Battie, T. Videman, and J. Kaprio Role of genetic and environmental influences on heart rate variability in middle-aged men Am J Physiol Heart Circ Physiol, August 1, 2007; 293(2): H1013 - H1022. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. J. C. De Geus, N. Kupper, D. I. Boomsma, and H. Snieder Bivariate Genetic Modeling of Cardiovascular Stress Reactivity: Does Stress Uncover Genetic Variance? Psychosom Med, May 1, 2007; 69(4): 356 - 364. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Circulation Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 2004 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |