The Relationship between Heart Rate Variability and Adiposity Differs for Central and Overall Adiposity.

While frank obesity is associated with reduced HRV, indicative of poorer autonomic nervous system (ANS) function, the association between body mass index (BMI) and HRV is less clear. We hypothesized that effects of adiposity on ANS are mostly mediated by visceral fat and less by subcutaneous fat; therefore, centrally distributed adipose tissue, that is, waist circumference (WC), should be more strongly associated with HRV than overall adiposity (BMI). To examine this hypothesis, we used data collected in a subset of the Baltimore Longitudinal Study of Aging to compare strength of association between HRV and WC to that of HRV and BMI. Time domain HRV variables SDNN (standard deviation of successive differences in normal-to-normal (N-N) intervals) and RMSSD (root mean square of successive differences in N-N intervals) were calculated from 24-hour Holter recordings in 159 participants (29-96 years). Increasing WC was associated with decreasing SDNN and RMSSD in younger but not older participants (P value for WC-by-age interaction = 0.003). BMI was not associated with either SDNN or RMSSD at any age. In conclusion, central adiposity may contribute to sympathetic and parasympathetic ANS declines early in life.

1. Introduction

The autonomic nervous system (ANS) directs voluntary and involuntary physiologic processes such as digestion, blood pressure, hormonal regulation, energy metabolism, and heart rate and is therefore considered an important regulator of homeostasis [1-3]. Heart rate variability (HRV) is widely considered a standard noninvasive method for assessing ANS function, due to ANS regulation of heart rate in a continuous, beat-to-beat manner [4], with lower HRV generally considered an indicator of poorer autonomic function. Lower HRV or autonomic dysfunction has been associated with several conditions of significant public health concern, including diabetes and diseases of the cardiovascular system [5-11].

Obesity, an important risk factor for diabetes and cardiovascular disease, is associated with dysregulation of autonomic function. It has been suggested that ANS dysfunction is an important mediator in the development of obesity-associated disease and insulin resistance although the nature of the link between adiposity and insulin sensitivity is still unclear [12-15]. Studies of the relationship between body mass index (BMI) and HRV have reported conflicting results with stronger associations observed in younger populations [16-19]. Understanding the mechanism that connects obesity and ANS function is important because of the increasing obesity prevalence documented among men and women of all ages [20]. One explanation for controversial data on the associations between BMI and HRV is that the distribution of body fat may be more important than a measure of overall obesity, such as BMI, due to the metabolic activity of visceral or abdominal adipose tissue [21-23]. This could be especially problematic in older persons for whom BMI might be a poorer indicator of adiposity [24]. We hypothesized that central (abdominal) adiposity, not overall adiposity, accounts for obesity-related autonomic dysfunction. To examine this hypothesis, we compared the strength of association between HRV and WC to that of HRV and BMI, anticipating stronger, negative associations between HRV and WC compared to that of HRV and BMI. Furthermore, since obesity tends to be a stronger risk factor for cardiovascular morbidity and mortality in young and middle-aged persons than in older persons, we examined whether WC and BMI relationships with ANS differ by age, independent of specific health conditions known to impact HRV.

2. Methods

2.1. Study Population

The Baltimore Longitudinal Study of Aging (BLSA) is an observational study sponsored by the National Institutes of Health and the National Institute on Aging (NIA) that began in 1958 to study normative aging in a volunteer cohort of healthy persons 18 years of age and older at study entry [25]. Most BLSA participants reside in the eastern USA, predominantly the Baltimore-Washington DC area, and visit the NIA Advanced Studies in Translational Research on Aging Unit every one to four years depending on their age. In 2004, 24-hour ambulatory electrocardiograms were implemented in the BLSA to assess HRV.

The first two hundred fourteen participants with Holter data analyzed consecutively in the BLSA were evaluated. These participants underwent a two-to-three day assessment between March 2004 and March 2006, which included a battery of tests including a standardized medical examination and interview, sociodemographic interview, anthropometric assessment, fasting blood tests, and an oral glucose tolerance test (OGTT) using a standard 75 gram glucose load. All participants signed informed consent. This study complied with the ethical rules for human experimentation as stated in the Declaration of Helsinki and was approved by the MedStar Research Institute institutional review board.

2.2. Outcomes: HRV Assessment

HRV was assessed from ambulatory 24-hour Holter monitors that use digitalized technology with compact flashcards and HRV software (DelMar Reynolds Lifecard CF, 1024 samples per second, Impresario software). An initial analysis was planned for the first 200 recordings processed, and 214 had been completed at the time of this analysis. All recordings were visually reviewed and manually edited in a beat-by-beat analysis using a standardized protocol to identify artifact and nonsinus beats.

Exclusion criteria for HRV analysis on Holter data were (1) implanted pacemakers, (2) transient or persistent nonsinus rhythm during the recording, such as atrial fibrillation, and (3) wandering atrial pacemakers (WAPs). Although WAP is generally considered a benign condition, it may produce invalid readings and constitutes a typical exclusion [26-28]. HRV measures can be affected by frequent premature beats, artifact, and total hours of recordings. Thus, HRV data were excluded when recordings were less than 18 hours, more than 10% of beats were premature, or artifact time exceeded 5% of recorded time [4].

This analysis focused on standard time domain HRV measurements, SDNN (standard deviation of successive differences in normal-to-normal (N-N) intervals), and RMSSD (root mean square of successive differences between N-N intervals). SDNN represents a measure of overall ANS function (sympathetic and some parasympathetic function) while RMSSD represents a measure of parasympathetic function [4]. Time domain measures have been suggested for ambulatory studies as they are less susceptible to respiratory and movement artifacts. Because HRV can be affected by activity and sleep, SDNN and RMSSD were obtained separately for day (6:00 AM–10:30 PM unless otherwise reported in diary and was of similar hours) and night (1:00 AM–5:00 AM). Nighttime hours were also selected to provide a constant length of time for HRV analysis in which participants were all sleeping.

The autonomic nervous system contributes to the circadian patterns of HRV with sympathetic nerve discharges associated with activity, such as during daytime wakefulness [29]. The circadian patterns of parasympathetic activity also demonstrate a peak which occurs at night and plateaus during the day [30]. Because of this, the HRV outcomes of interest will focus on daytime SDNN, to describe overall function, sympathetic more so than parasympathetic function, and nighttime RMSSD, to describe parasympathetic function.

2.3. Predictors

Height, weight, and waist circumference, (WC, centimeter, one inch above the anterior superior iliac crests) were measured with participants disrobed or in light clothes and barefoot. Height was measured using a stadiometer. WC and height were measured three times and averaged. BMI was calculated as weight in kilograms/height in meters2 as a measure of overall adiposity whereas WC was considered a measure of central or visceral adiposity. This was based on findings from the Rancho Bernardo Heart and Chronic Disease Study, which showed WC to be strongly correlated with most measures of central obesity in men and women aged 65–96 years [31].

2.4. Covariates

Conditions known to be associated with HRV and body composition were assessed as covariates, including glucose intolerance, hypertension, and physical activity [10, 12, 32]. Glucose metabolism was assessed with a 2-hour oral glucose tolerance test (OGTT). Blood pressure was taken using a manual sphygmomanometer with the participant resting supine for 5 minutes in a quiet, climate-controlled room with low lighting, using the first and fifth Korotkoff sounds. The second and third systolic blood pressure (SBP) and diastolic blood pressure (DBP) measurements were averaged. Hypertension was defined as use of antihypertensive medication, self-reported diagnosis confirmed with medical record review, or SBP > 140 mm Hg or DBP > 90 mm Hg during the physical examination. Diabetes was defined as use of hypoglycemic medications, self-reported diagnosis of diabetes confirmed by medical records, fasting glucose of 126 mg/dL or greater, or a 2-hour OGTT value ≥200 mg/dL [33]. Angina and myocardial infarction were assessed by self-report, review of medical records, and by major Q-waves on resting electrocardiograms using standard Minnesota Coding system criteria. Due to the infrequent occurrence of diabetes, angina, and myocardial infarction, analyses were conducted including and then excluding participants with these conditions. Self-reported physical activity was categorized by metabolic equivalents per kg per week [34]. Age and race were obtained from a standardized interview. Race was categorized as black (self-reported race as “black” or “African American”) and nonblack.

2.5. Statistical Analysis

Simple comparisons of continuous variables by sex were performed with nonparametric Wilcoxon tests and categorical variables with chi-square tests. The correlation between continuous variables was assessed using Spearman's correlation coefficient. Because BMI and WC have different units, we used standardized coefficients of each to assess whether one was associated more with the HRV outcomes than the other. Akaike's Information Criterion (AIC) was used to compare models (WC versus BMI), with lower AIC indicating a better fitted model. To account for nonnormal distributions of HRV, standard generalized linear models with gamma distribution and log link were used to assess relationships between adiposity measures (BMI and WC) and HRV (nighttime RMSSD and daytime SDNN). To determine whether these relationships differ by age, age-by-BMI and age-by-WC interaction terms were included. All models were adjusted for age, sex, race, physical activity, and hyperglycemia using the OGTT 2-hour glucose. Previously published studies suggest that HRV may differ across age by sex and race, and some studies suggest that glucose metabolism differs by sex [35, 36]. Thus, we also [34] evaluated these interactions by including age-by-sex, age-by-race, and sex-by-OGTT glucose terms. Significance level was defined as P < .05.

Relationships were reported as relative ratios (RR), to describe the multiplicative increases and decreases in HRV associated with a unit change in the explanatory variable, similar to a relative risk for binary data. For example, an RR of 1.2 for male sex suggests a 20% higher HRV in men than in women.

Lastly, because of the potential for beta blockers to improve HRV [37], we repeated all analyses after excluding participants on beta blockers. All analyses were conducted using SAS 9.1 statistical software (SAS Institute Inc., Cary, NC).

3. Results

Of 214 participants with Holter data, we excluded two participants with recordings less than 18 hours, one with nonsinus rhythm who also had a pacemaker, and 33 with wandering atrial pacemakers, leaving 178 with HRV data. Among these, nine were excluded from the OGTT (two were taking insulin, six used corticosteroids in the previous three months, and one refused). Of the 169 remaining participants, three were missing data on WC and seven were missing OGTT glucose due to difficult intravenous access, leaving 159 participants for the current analysis. Of these, one participant lacked a nighttime assessment. Only 15 participants (8%) were diabetic and 8 (4.5%) had angina or myocardial infarction; 116 (65%) were hypertensive. Sixty-seven participants (37%) were 70 years or older.

Characteristics of the 178 participants with valid HRV data are presented in Table 1. Figure 1 shows histograms of the HRV variables. The similarity between the lines representing the gamma distribution and the kernel density smoother of the actual data suggests a good fit of the gamma model.

Table 1

Sample characteristics stratified by sex.

	Women (n = 90)	Men (n = 88)	P value
	Mean (SD) or N (%)
SDNN (ms2)	103 (25.6)	110 (30.0)	0.096
RMSSD (ms2)	30 (14.6)	30 (16.7)	0.815
Age years	64 (14.3)	66 (11.8)	0.198
Systolic blood pressure (mm Hg)	119 (16.6)	123 (14.2)	0.107
Body mass index (kg/m2)	27 (4.9)	27 (3.9)	0.226
Waist circumference (cm)	88 (14.0)	98 (12.9)	<0.001
Fasting glucose (mg/dL)	93 (13.1)	98 (17.7)	0.026
OGTT 2-hour glucose (mg/dL)	118 (43.5)	133 (54.7)	0.067
Black race†	26 (29%)	24 (26%)	0.681
Antihypertensive use	40 (44%)	46 (52%)	0.371
Beta blocker use	16 (18%)	21 (24%)	0.317
Hypertension	50 (56%)	56 (63%)	0.353

†Self-reported race/ethnicity (N) included: Caucasian (118), black or African American (50), Asian or Pacific islanders (4), American Indian or native Alaskan (1), Hispanic (2), Chinese (3), and other or unknown (1).

ms: milliseconds; OGTT: oral glucose tolerance test; RMSSD: root mean square of successive differences in normal-to-normal intervals; SDNN: standard deviation of normal-to-normal intervals.

Figure 1

Histograms of heart rate variability variables with a smooth kernel density estimate (dashed line) and gamma distribution (solid line). RMSSD: root mean square of successive differences; SDNN: standard deviation of N-N intervals.

WC and BMI were significantly correlated (r = 0.73, P < .0001). Despite this high correlation, a wide range of WC was observed for any given BMI, suggesting that WC describes potentially different information than BMI in this population.

3.1. RMSSD (Parasympathetic ANS Function)

The relationship between RMSSD and WC differed by age, with RMSSD decreasing as WC increased in younger but not older participants (P for interaction = 0.008) (Table 2 and Figure 2, left panel). Figure 2 illustrates a comparison of the expected nighttime RMSSD, based on model results, for an average participant at two different ages: 45 years and 85 years if the standardized WC increased. With increasing WC, the RMSSD decreases in the 45-year-old but not in the 85-year-old participant. In contrast, neither the main effect of BMI nor the BMI-by-age interaction was statistically significant, suggesting that BMI was not associated with RMSSD at young or old ages (P = 0.403) (Figure 2, right panels). Furthermore, the AIC for the RMSSD models suggested that WC provided a better fitted model than did BMI (AIC = 819 versus 842).

Table 2

Relative risks (RR) percent change in SDNN and RMSSD for each standard deviation increase in the obesity variables: waist circumference (WC) and body mass index (BMI). Smaller Akaike's Information Criterion (AIC) indicates better fitted models.

A. WC model	Night RMSSD	Day SDNN
	RR	RR
	(95% CI) P value	(95% CI) P value
Age (per 10 year)	0.41	0.64
	(0.23, 0.75) 0.004	(0.48, 0.85), 0.002
Waist circumference (per 5 cm)	0.73	0.85
	(0.58, 0.92) 0.008	( 0.76, 0.95), 0.005
Age-by-waist circumference	1.05	1.02
	(1.01, 1.08) 0.008	(1.01, 1.04), 0.004
AIC	819	953
B. BMI model	Night RMSSD	Day SDNN
	RR	RR
	(95% CI) P value	(95% CI) P value
Age (per 10 year)	1.10	0.89
	(0.71, 1.68), 0.673	(0.72, 1.10), 0.300
Body mass index (BMI)	1.04	0.98
	(0.93, 1.17), 0.496	(0.92, 1.03), 0.404
Age-by-BMI	0.99	1.00
	(0.98, 1.01), 0.403	(0.999, 1.01), 0.428
AIC	842	982

Models adjusted for sex, race, blood pressure, oral glucose tolerance test (OGTT) 2-hour glucose, age-by-race, age-by-sex, OGTT 2-hour glucose-by-sex, and physical activity.

RR:relative risk;SDNN:standard deviation of normal-to-normal intervals.RMSSD:root mean square of successive differences in normal-to-normal intervals.

Figure 2

Data points are observed data; lines represent the predicted value of nighttime RMSSD (top row) and daytime SDNN (bottom row) (if the adiposity measure increased for an average participant at two ages: 45 years (solid lines) and 85 years (dashed lines)). Standardized waist circumference and standardized BMI are shown on x-axes to scale the measures similarly. The figures show differing relationships between waist circumference and HRV measures at young and old ages but not between BMI and HRV measures. Negative relationships were observed between waist circumference and RMSSD and SDNN in the 45 yo but no significant relationship in the 85 yo. Comparatively, the SDNN-BMI relationship is attenuated, compared to WC, for the 45-year olds and the RMSSD-BMI relationship is relatively flat at 45 and 85 years, with neither being statistically significant. *x-axes are standard deviations from mean waist circumference (women: 88 cm, men: 98 cm) or BMI (27 kg/m2). *Graded gray-scale data points from population illustrate increasing age (lighter is younger, darker is older, 26–96 y). Models were adjusted for sex, race, OGTT glucose, blood pressure, and the interactions age-by-sex, age-by-race, the sex-by-OGTT glucose. OGTT glucose = 2-hour oral glucose tolerance test glucose result; RMSSD:root mean square of successive differences; SDNN:standard deviation of N-N intervals.

3.2. SDNN (Overall ANS Function)

Increasing WC was also associated with decreasing SDNN in younger but not older participants (P for interaction = 0.003) (Figure 2, Table 2). There was no statistically significant association between BMI and SDNN (Figure 2, Table 2). In addition, the model fit for WC was superior to the model fit for BMI (AIC = 953 versus 982).

When participants with coronary artery disease or diabetes or those taking beta blockers were excluded, we found similar results for both SDNN and RMSSD. Specifically, the parameter estimates remained similar, and statistical significance was unchanged despite a smaller sample (data not shown).

4. Discussion

To our knowledge, this is the first study to assess the roles of overall and abdominal adiposity on ANS function in young and old persons. The findings showed that central adiposity, measured by waist circumference, was associated with poorer ANS function at younger but not older ages. This relationship was observed for a measure of overall ANS function that includes sympathetic and parasympathetic activity and for a measure of primarily parasympathetic activity. Overall adiposity, measured using BMI, was not associated with either measure of ANS function. Thus, abdominal adiposity, as opposed to overall adiposity tissue, could adversely affect ANS function at younger ages with potentially unfavorable effects in later life.

Our results are in accordance with a small study by Christou et al. [38], which reported higher abdominal-to-peripheral body fat distributions measured by DXA were strongly correlated with lower sympathetic and parasympathetic function in young (mean age 25 ± 1 year) and old (mean age 65 ± 1 year) healthy men. Other studies have also shown that body fat distribution accounts for obesity-associated risk more so than total fat, with visceral fat carrying the greatest risk for insulin resistance, hyperglycemia, hypertriglyceridemia, cardiovascular disease, metabolic syndrome, and mortality [22, 23, 39–41]. Our findings suggest that abdominal adiposity may also contribute to poorer sympathetic and parasympathetic function as well.

The relationship between central adiposity and ANS function is not a surprising finding. The ANS innervates fat depots, which are positively associated with catecholamine production by the ANS [42, 43]. The biologic mechanisms whereby visceral fat could contribute to ANS dysfunction could involve adipocytokines secreted by fat cells [44, 45]. Although metabolic activity of central adipose tissue could account for the associations between ANS function and obesity, lack of association with BMI at older ages could be explained by the poor ability of BMI to accurately reflect body fat at older age [24]. This does not explain, however, the lack of associations observed in young persons. Future studies should examine the role of adipocytokines in mediating the effects of obesity on ANS function.

The lack of association for either measure of obesity at older ages supports recent arguments that “a little extra” around the middle may be good for older adults, or at least not as harmful as in young and middle age [46, 47]. Our findings lend some support for this or, at least, fail to demonstrate that central adiposity, within the range observed in this population, was associated with poorer ANS function in the oldest participants. In contrast, younger participants with more central adiposity appeared to have ANS function similar to older participants, that is, premature aging of the ANS.

The long-term sequelae of impaired ANS function manifested as lower HRV are unknown but could contribute to cardiovascular, lipid, and endocrine abnormalities as these mechanisms are regulated by the ANS. Evidence for this includes studies that demonstrate associations between lower HRV and coronary artery disease, heart failure, hyperlipidemia, hypertension, and prediabetic and diabetic states [9–11, 13, 48, 49]. Additionally, the relationship between central adiposity and ANS function observed in this study independent of blood pressure and glucose could explain how three of the most important components of the metabolic syndrome (i.e., hypertension, diabetes, and obesity) act synergistically on the cardiovascular system—that is, through the derangement of ANS function.

Regarding limitations, this study included a relatively small sample in a cross-sectional analysis. However, the BLSA is an ongoing study so future analyses will assess the robustness of these findings both cross-sectionally and longitudinally. Study participants were in relatively good health. Consequently, findings may not be generalizable to persons in poorer health. Nevertheless, observing the association between WC and HRV in healthy persons suggests that ANS dysfunction can be detected early and at relatively young ages, therefore providing an opportunity for early intervention. The cross-sectional design limits our ability to establish causal relationships, and, in fact, one hypothesis proposes that central effects of adipocytokines could alter sympathetic and parasympathetic balance, leading to changes in the distribution of visceral and subcutaneous fat [50]. Nevertheless, our findings provide important clues to guide future studies. For instance, longitudinal analyses will be conducted to determine how changes in WC and ANS function at younger ages contribute to health risk in late life.

In summary, we observed an inverse association between central adiposity measured by WC and ANS function at young and middle but not old-age. This was true for a measure of primarily parasympathetic function (RMSSD) as well as a measure of mixed sympathetic and parasympathetic function (SDNN). BMI was not associated with ANS function in this study. Future studies should examine whether metabolic activity of central adipose tissue explains these associations.