Abstract
Autonomic imbalance, a feature of both diabetes and hypertension, may contribute to adverse cardiovascular outcomes. In animal models, sympathetic nerve activity contributes to renal damage but the extent to which autonomic dysfunction precedes the development of CKD and ESRD in humans is unknown. We measured resting heart rate and heart rate variability in 13,241 adults (45- to 64-years old) followed for a median of 16 years in the Atherosclerosis Risk in Communities (ARIC) Study. We examined heart rate parameters by quartiles, defining those in the lowest quartile (by time and frequency domain measures separately) as the risk group of interest. We identified 199 cases of incident ESRD and 541 patients with CKD-related hospitalizations; higher resting heart rate and lower heart rate variability associated with both outcomes. The fully adjusted hazard ratios for ESRD were 1.98 (95% confidence interval [CI] 1.45 to 2.70) among those in the highest heart rate quartile and 1.56 (95% CI 1.14 to 2.14) for high-frequency power. Other time and frequency domain measures were similarly and significantly associated with ESRD and CKD-related hospitalizations. These results suggest that autonomic dysfunction may be an important risk factor for ESRD and CKD-related hospitalizations and call for further studies to define the mechanisms that underlie these associations.
Even in the absence of hypertension and diabetes mellitus, 12% of adults aged 65 years and older have chronic kidney disease (CKD), suggesting that other mechanisms of renal injury may be important in the general population.1 An infrequently explored potential cause of renal injury is autonomic imbalance (high sympathetic tone and/or low parasympathetic tone). Histologic studies demonstrate that renal sympathetic nerve terminals are in direct contact not only with renal vasculature but also with tubules and juxtaglomerular cells.2 In animal models, alterations in the activity of these nerves modulate renal hemodynamics, tubular transport, and renin secretion.2–4 In animals, sympathetic denervation results in a natriuretic and diuretic response and attenuates the progression of renal failure.2–4 In humans, several phenomena associated with autonomic imbalance, such as impaired diurnal BP variation, have been linked to chronic kidney disease and its progression, independent of mean BP and diabetes.5–7 Whether the autonomic imbalance associated with abnormal diurnal BP variation is the physiologic driver of this association remains unclear.8,9
Although there is no gold-standard measure of autonomic balance,9 cardiac sympathetic and parasympathetic tone can be assessed using heart rate variability (HRV) measurements—a practical way to measure autonomic balance using standard electrocardiographic monitoring. By examination of the average heart rate, and quantification of beat-to-beat variability (the SD of a patient's RR intervals, for example), HRV can be quantified. Additionally, computer algorithms that utilize Fourier (spectral) analysis have been developed to quantify cyclic changes in heart rate that occur during the course of a few seconds (high-frequency [HF] power) or during longer time periods (low-frequency [LF] and very low frequency power). For example, in most healthy young adults, resting heart rate will predictably accelerate and decelerate with the respiratory cycle (every few seconds—a high frequency). The more the heart rate varies with respiration, the higher the HRV, and in this case, the higher the “high-frequency power.” A failure to exhibit this variation would be characterized as low HRV and is a marker of autonomic dysfunction. In general, lower resting heart rates and greater beat-to-beat variability in heart rate (high HRV) are associated with a healthy autonomic nervous system and good cardiovascular health. High resting heart rate and low HRV are associated with a host of adverse cardiovascular outcomes, as well as with precursors of cardiovascular disease, including features of the metabolic syndrome.10–19
Several small cross-sectional studies have found that patients with CKD have decreased HRV relative to those without CKD.20–28 However, it is not known to what extent autonomic imbalance precedes the development of CKD and may contribute to its development because reverse causality—nephrosclerosis causing autonomic imbalance—may be the primary mechanism for the association. If autonomic imbalance precedes the development of chronic kidney disease, it may serve as a marker to identify patients at higher risk of developing ESRD. More important, it may serve to provide evidence that the autonomic nervous system plays a role in the pathophysiology of nephrosclerosis. Therefore, using data from the Atherosclerosis Risk in Communities (ARIC) study, we sought to determine the relationship between low HRV and the subsequent development of adverse renal outcomes.
Results
Baseline Characteristics of Patients and Ultimate Renal Outcomes
Of 13,241 patients included, 199 developed ESRD and 541 had a CKD-related hospitalization over a median of 16-year follow-up. Baseline characteristics are shown in Table 1. Baseline estimated GFR (eGFR) was moderately lower in those who ultimately developed ESRD (77 ml/min per 1.73 m2) than in the total cohort (93 ml/min per 1.73 m2). More than half of those who developed ESRD were black, although blacks comprised only about a quarter of the total cohort. Diabetes was present in 62% of those who developed ESRD. Systolic and diastolic BP was higher in those who went on to develop ESRD, and treatment with antihypertensives was more prevalent as well. Features of the metabolic syndrome, including obesity, high triglycerides and low HDL cholesterol were also more prevalent in those who ultimately developed renal failure.
Baseline characteristics of the cohort of 13,241 by eventual renal outcome
Baseline characteristics of those with low HRV (quartile 1) versus those with higher HRV are shown in Table 2a (time domain measures) and Table 2b (frequency domain measures). Low HRV (quartile 1) was more prevalent in those with diabetes and hypertension, but the mean eGFR was 90 to 95 ml/min per 1.73 m2 for all HRV categories. Higher body mass index (BMI) was modestly but significantly associated with time domain measures, but not with frequency domain measures. Triglycerides were consistently highest in those with lower HRV.
Characteristics of the cohort by baseline HRV: time domain and heart rate measures
Characteristics of the cohort by baseline heart rate variability: frequency domain measures
As shown in Table 3, in univariate analysis, each of the HRV metrics was strongly associated with incidence rates of ESRD and CKD-related hospitalizations. There appeared to be a threshold effect such that those in the lowest quartile (Q1) of HRV measures were at increased risk for adverse renal outcomes, but the event rates among those in the remaining three quartiles (Q2 to Q4) were comparable. This is illustrated in Figure 1 and Figure 2. There was a suggestion of a dose-response relationship: among the 745 patients who were simultaneously in the worst quartile for heart rate and LF and HF power, the unadjusted hazard ratio (HR) for ESRD was 10.15 (6.95 to 14.87) and the fully adjusted hazard ratio was 3.39 (2.12 to 5.42), compared with the 6979 patients who were not in the worst quartile for any of these three measures.
Events and incidence rates of ESRD and CKD hospitalization events by quartiles of HRV measures
Those in the lowest HRV quartile are at increased risk of renal failure. Time to ESRD (Kaplan-Meier estimates) by quartiles of heart rate measures. Time to ESRD by quartiles of HRV measures: (A through C) time domain measures; (D through F) frequency domain measures. RMS, root-mean squared.
Those in the lowest HRV quartile are at increased risk of hospitalization. Time to kidney-related hospitalization (Kaplan-Meier estimates) by quartiles of heart rate measures. Time to kidney-related hospitalization by quartiles of HRV measures: (A through C) time domain measures; (D through F) frequency domain measures.
Unadjusted and adjusted hazard ratios associating low HRV with adverse renal outcomes are presented in Table 4. There was a moderately strong association between low HRV and both ESRD and CKD-related hospitalizations that remained significant after adjustment for other factors associated with adverse renal outcomes, including age, gender, race, baseline eGFR, systolic and diastolic BP, diabetes, use of antihypertensive drugs, baseline coronary disease, smoking, lipid variables, fasting insulin and glucose, BMI, and waist circumference.
HR associated with lowest quartiles of time and frequency domain measuresa
Although baseline albuminuria values were not available, they were available at the time of the third follow-up visit (approximately 9 years after study inception). Using these measures, we performed sensitivity analyses that excluded all patients who ultimately developed albuminuria,29 assuming that this would also exclude the vast majority of patients who had albuminuria at baseline. Even in this small subset, which included only 16 ESRD cases and 134 CKD-related hospitalization cases, there remained significant adjusted associations between some of the HRV measures and ultimate renal outcomes. Specifically, the fully adjusted HR for ESRD in patients with heart rates in the highest quartile was 3.17 (1.13 to 8.91) and was 4.84 (1.61 to 14.6) for those with low root mean square of successive differences in normal-to-normal RR intervals (rMSSD). The other HRV measures were not significantly associated with adverse renal outcomes in the fully adjusted models in this subset of patients.
Additional sensitivity analyses varied the baseline patients, modeling of baseline kidney function (continuous eGFR), and outcome definitions. Associations between HRV and risk of ESRD in sensitivity analyses that excluded all individuals with prevalent CKD (baseline eGFR <60 ml/min per 1.72 m2) and modeled baseline eGFR continuously (spline, with knots at eGFR values of 60 and 90 ml/min per 1.72 m2) were similar to the presented results. Varying the outcome to include less severe kidney dysfunction also yielded similar results. We also examined whether there was an interaction between the presence or absence of diabetes mellitus and the association between HRV measures and renal outcomes, but no significant interactions were found. When we excluded patients on antihypertensive medications at baseline, the relationship between impaired HRV and ESRD was modestly accentuated, as shown in Table 5.
HR for ESRD associated with HRV measures among those patients NOT on antihypertensive medications at baseline (n = 9947)a
Discussion
In this longitudinal study, we found that low HRV was associated with the subsequent development of renal impairment. This association remained significant after adjustment for other factors known to contribute to the development of renal failure, including diabetes, hypertension, baseline renal function, lipid variables, and parameters associated with insulin resistance and obesity. Although this does not demonstrate a cause-and-effect relationship,8 the temporal pattern is intriguing in that abnormal HRV preceded the development of clinical renal failure by many years. Furthermore, there are plausible mechanisms by which abnormal autonomic balance may lead to kidney damage.
Sympathetic nerve terminals innervate the kidneys directly, potentially affecting tubular function by enhancing solute and fluid resorption and modifying renal microvascular function by enhancing the effects of angiotensin.2,3,30–32 Additionally, the global health of the vasculature is affected by the autonomic nervous system's role in regulating hemodynamics, vascular tone, metabolism, and inflammation.9,33–37 This is relevant because the glomerulus is adversely affected by the same pathophysiological factors that lead to atherosclerosis, including hypertension, endothelial dysfunction, dyslipidemia (particularly high triglycerides and low HDL cholesterol), insulin resistance, inflammation, and oxidative stress.4,38–42 Given that these risk factors for cardiovascular disease are also associated with high sympathetic and low parasympathetic tone,19,43–45 it is not surprising that HRV and other measures of autonomic tone are related to the progression of renal dysfunction.
Our study has several limitations. First, baseline urinary albumin excretion was not measured. Therefore, it is difficult to exclude the possibility of reverse causality—early renal disease leading to dysregulated autonomic tone—because some patients may have had albuminuria at baseline despite having preserved eGFR. However, our sensitivity analyses suggested that even among those with normal eGFR at baseline and those who did not have albuminuria midway through follow-up, there was still a relationship between HRV and the development of clinical renal complications. We lacked serum creatinine measurements during the 3-year follow-up visits, which impaired our ability to assess CKD progression. Additionally, about 2000 patients enrolled in ARIC did not have baseline HRV measurements, and the HRV measurements themselves were limited by the absence of postural stimulation and longer recording intervals. Although continuous hospitalization surveillance was used to define incident ESRD and vital data, only events occurring in acute care hospitals were investigated, and events occurring in other institutions, such as nursing homes, might be missed. In addition, there is a possibility that hospitalizations occurred outside of the study area. We cannot exclude the possibility of residual confounding, particularly given that adjustment for other cardiovascular risk factors attenuated the relationship between HRV and renal outcomes; variables that we did not measure and analyze include markers of baseline inflammation, nutrition, and physical fitness. Finally, we acknowledge that heart rate and heart rate variability are affected by factors other than endogenous autonomic tone (such as medications that affect heart rate) and that there are other ways to examine autonomic tone (such as direct measurement of circulating norepinephrine levels, heart rate recovery after exercise, and direct measurement of sympathetic nerve electrical output).9 We hope that other researchers will use different cohorts to determine whether non-HRV markers of autonomic dysfunction predict subsequent renal dysfunction.
The clinical implications of our findings are uncertain. Although there is clarity that reducing BP by any number of means can retard the progression to renal failure,46,47 drugs that block the renin-angiotensin-aldosterone system remain the mainstay for preventing the progression to nephropathy in patients with diabetes or early renal damage.48 Some of the renoprotective effects of these agents may be mediated in part by cross-talk between the sympathetic nervous system and the renin-angiotensin-aldosterone system.49 Although animal models indicate that sympathetic denervation of the kidneys can provide renoprotective effects,50 there is little evidence that centrally or peripherally acting sympatholytics have unique renoprotective effects in humans.51 Rather than suggesting a change in clinical practice, we hope our findings will encourage further research to better define the putative role of the autonomic nervous system in precipitating and exacerbating renal disease in humans; this, in turn, may ultimately lead to novel therapeutic approaches once the mechanisms for our findings are better characterized.
Concise Methods
Study Cohort
The ARIC Study is a prospective observational cohort of 15,792 individuals aged 45 to 64 years, drawn from four U.S. communities. Examinations included a baseline visit between 1987 and 1989, follow-up examinations approximately every 3 years, and annual telephone interviews. Hospitalization events were ascertained through December 31, 2004. Further details of the ARIC cohort have been previously published.52
For this study, we excluded participants lacking baseline serum creatinine values (n = 149), reported a race other than white or black (n = 49), and excluded blacks from the Minnesota and Washington County study centers (n = 55) based on low numbers of blacks at these sites. We also excluded participants missing important exposure and covariate information (n = 2196 missing baseline HRV measures and 102 missing baseline hypertension or diabetes information). Analyses were based on the remaining 13,241 study participants.
Data Collection
Demographic data, medical history, and measurements of height, weight, and BP were obtained during each clinical examination. Laboratory measures were performed as described previously.52 Diabetes mellitus was defined as a fasting glucose of ≥126 mg/dl, nonfasting glucose of ≥200 mg/dl, self-reported physician diagnosis of diabetes, or use of diabetic medications.
Seated BP measurements were taken after 5 minutes of rest. Standard methods were used to determine lipid levels.53 Cigarette smoking was determined by self-report (current, former, or never). Coronary heart disease was defined as physician-diagnosed myocardial infarction, evidence of prior myocardial infarction by electrocardiogram, or prior coronary revascularization. eGFR was calculated using the four-variable Modification of Diet in Renal Disease (MDRD) Study equation.54 Urinary albumin was measured from a spot urine sample at visit 4 (1996 through 1998) as the ratio of albumin to creatinine (ACR). Albuminuria was defined as ACR ≥30 mg/g.
HRV was measured from 2-minute beat-to-beat heart rate recordings taken in a supine position at baseline. Heart rate data were analyzed by time domain analysis and frequency domain analysis as described previously.55 Time domain and heart rate measures included mean normal-to-normal RR interval length (which we converted to heart rate in beats per minute for ease of interpretation), the SD of all normal-to-normal RR intervals, and the rMSSD. For frequency domain measures, a fast Fourier transform algorithm was used.19 Spectral analysis was used to calculate power (area under the power spectral density curve), although it was not technically feasible to do so for 642 patients. Given the relatively short duration of these recordings, we report only HF and LF power (or their ratio) and do not report very low frequency power. HF power was defined as the total area between 0.15 and 0.4 Hz (cycles between 2.5 and 6.7 seconds), which is thought to correspond to the heart rate variation that results with normal respiration, mediated predominantly by the vagus nerve.56 LF power was defined as the total area between 0.04 and 0.15 Hz, which is thought to relate to changes in sympathetic neural output.56 Reproducibility of each of these measures has been shown to be high.57
Outcome Assessment—Hospitalization and ESRD
Deaths and hospitalizations through 2004 were identified via annual participant interviews, local hospital discharge lists, and county death certificates and included all those coded (International Classification of Diseases, Ninth Revision [ICD-9]) for chronic renal disease (581 to 583.91, 585 to 588.91), hypertensive renal disease (403 to 403.91), hypertensive heart and renal disease (404 to 404.93), unspecified disorder of kidney and ureter (593.9), diabetes with renal manifestations (250.40 to 250.43), kidney transplant, renal dialysis or adjustment/fitting of catheter (V42.0, V45.1, or V56), or either hemodialysis (39.95) or peritoneal dialysis (54.98) without simultaneous acute renal failure (584, 586, 788.9, and 958.5). Also included as cases were participants with acute renal failure (ARF) as a cause of death, if they had an earlier diagnosis of CKD.58,59 Excluded were individuals that had a transplant or dialysis code on the same date as another ARF code (586, 584, 788.9) without previous CKD and those with an ARF code of 958.5 (traumatic anuria). Corresponding ICD-10 codes were used for deaths after ICD-10 implementation.
Cases of ESRD were assessed through continuous surveillance of hospitalizations through December 31, 2004, and included ICD codes specified for kidney transplant, dialysis, or a procedural code indicating dialysis. Participants with an underlying cause of death of ARF, conditional upon having an earlier diagnosis of CKD (as defined by the ICD-9 codes above), and excluding individuals that have had a transplant or dialysis code on the same date as another ARF code (586, 584, 788.9) without previous CKD (and excluding anyone with an ARF code of 958.5) were defined as ESRD cases.
Statistical Analysis
Baseline characteristics of the cohort were compared across ESRD status, CKD hospitalization status, and HRV quartiles using χ2 and t tests. Follow-up time was calculated from baseline to the earliest date of ESRD (or CKD hospitalization); participants were censored at the earliest of the following: time of death, withdrawal, or December 31, 2004. Crude incidence rates and their 95% confidence intervals for the time to development of ESRD, or hospitalization, were computed using Poisson regression models.
On the basis of a prospectively defined analysis, HRV measures were divided into two groups with one group at or below the 25th percentile and the other group above the 25th percentile.17 Because we converted the RR interval to heart rate for ease of interpretation, the corresponding first quartile for the RR interval corresponds to the highest quartile of heart rate. We conducted Cox proportional hazards regression to examine the relationship between the lowest quartile of HRV and incident ESRD or CKD hospitalization. Multivariable models included age, gender, a race/study center combined variable, baseline eGFR category (defined as eGFR <45, 45 to 60, 60 to 75, 75 to 90, 90 to 120, or >120 ml/min per 1.73 m2), systolic BP, diastolic BP, use of antihypertensive agents, prevalent coronary heart disease, smoking, diabetes, total cholesterol, HDL cholesterol, and triglyceride concentrations. Additional models included fasting insulin and glucose, BMI, and waist circumference. Models including LDL cholesterol in place of total cholesterol were examined, but did not substantially change results so were dropped from final models. Sensitivity analyses, which included eGFR as a continuous variable, were performed.
Disclosures
None.
Acknowledgments
The Atherosclerosis Risk in Communities Study is carried out as a collaborative study supported by National Heart, Lung, and Blood Institute contracts N01-HC-55015, N01-HC-55016, N01-HC-55018, N01-HC-55019, N01-HC-55020, N01-HC-55021, and N01-HC-55022. The authors thank the staff and participants of the ARIC study for their important contributions.
We acknowledge methodological consultation from Hsin-Chieh Yeh, Ph.D., Director of the Johns Hopkins General Internal Medicine Methods Core.
This work was supported in part by grants T32-HL-007024 (L.D.B.), 1KL2RR025006 (R.Q.), and 1KL2RR025006-01 (D.C.) and Johns Hopkins Hospitalist Scholars Program (D.J.B.)
The Atherosclerosis Risk in Communities Study is carried out as a collaborative study supported by the National Heart, Lung, and Blood Institute contracts N01-HC-55015, N01-HC-55016, N01-HC-55018, N01-HC-55019, N01-HC-55020, N01-HC-55021, and N01-HC-55022.
Footnotes
Published online ahead of print. Publication date available at www.jasn.org.
- Copyright © 2010 by the American Society of Nephrology