Abstract
The single leading cause of mortality on hemodialysis is sudden cardiac death. Whether measures of electrophysiologic substrate independently associate with mortality is unknown. We examined measures of electrophysiologic substrate in a prospective cohort of 571 patients on incident hemodialysis enrolled in the Predictors of Arrhythmic and Cardiovascular Risk in End Stage Renal Disease Study. A total of 358 participants completed both baseline 5-minute and 12-lead electrocardiogram recordings on a nondialysis day. Measures of electrophysiologic substrate included ventricular late potentials by the signal-averaged electrocardiogram and spatial mean QRS-T angle measured on the averaged beat recorded within a median of 106 days (interquartile range, 78–151 days) from dialysis initiation. The cohort was 59% men, and 73% were black, with a mean±SD age of 55±13 years. Transthoracic echocardiography revealed a mean±SD ejection fraction of 65.5%±12.0% and a mean±SD left ventricular mass index of 66.6±22.3 g/m2.7. During 864.6 person-years of follow-up, 77 patients died; 35 died from cardiovascular causes, of which 15 were sudden cardiac deaths. By Cox regression analysis, QRS-T angle ≥75° significantly associated with increased risk of cardiovascular mortality (hazard ratio, 2.99; 95% confidence interval, 1.31 to 6.82) and sudden cardiac death (hazard ratio, 4.52; 95% confidence interval, 1.17 to 17.40) after multivariable adjustment for demographic, cardiovascular, and dialysis factors. Abnormal signal–averaged electrocardiogram measures did not associate with mortality. In conclusion, spatial QRS-T angle but not abnormal signal–averaged electrocardiogram significantly associates with cardiovascular mortality and sudden cardiac death independent of traditional risk factors in patients starting hemodialysis.
CKD has an estimated prevalence in the general population of 12.1%. Approximately 20% of these patients progress to ESRD, with a 3-year mortality of nearly 50%.1–3 The most common cause of death is cardiovascular disease, with sudden cardiac death (SCD) being the largest proportion of these events.4,5 Understanding the mechanisms and risk factors of cardiovascular causes of mortality in ESRD is critical to improving survival.
Despite the increased risk of mortality from SCD, there are few systematic investigations of the electrophysiologic substrate in ESRD with longitudinal follow-up. Most studies focus on traditional cardiovascular risk factors, such as cholesterol and hypertension, which are not consistently associated with sudden death in the context of hemodialysis unlike in the general population, thus highlighting the need to better understand the pathogenesis of arrhythmias and SCD among patients on dialysis.6,7
A number of parameters probe the electrophysiologic substrate and have been used to risk stratify those at highest risk for SCD. The presence of late potentials on a signal-averaged electrocardiogram (SAECG) characterizes a scar-related area of slow conduction in the ventricles of the heart and is associated with re-entrant ventricular tachycardia.8 Wide QRS-T angle characterizes abnormal ventricular conduction and global electrical heterogeneity in patients at risk of SCD.9 Often, the presence of these abnormal electrocardiogram (ECG) metrics is amplified in those with left ventricular hypertrophy (LVH), which is also independently linked to SCD.10
The goal of this study is to determine and compare the association of spatial QRS-T angle and late potentials on SAECG with mortality in a large prospective cohort of patients on incident hemodialysis.11 Evaluations of electrophysiologic markers specifically in an incident rather than prevalent dialysis population may yield more promising targets for interventions to decrease risk of mortality, which is highest in the first few months on dialysis.
Results
In total, 571 patients on incident hemodialysis were enrolled in the Predictors of Arrhythmic and Cardiovascular Risk in End Stage Renal Disease (PACE) Study, and 400 completed the initial cardiovascular study visit. After exclusion of participants without complete digital SAECG measurement, 358 patients comprised the final study population. Demographic and clinical characteristics of patients are presented in Table 1. The cohort had a mean age of 55±13 years, 73% were black, 59% were men, 58% were patients with diabetes, and only 37% had prevalent coronary artery disease.
Comparison of clinical and electrocardiographic characteristics of participants on hemodialysis with QRS-T angle ≥75° and participants on hemodialysis with QRS-T angle <75°
During 864.6 person-years of follow-up, there were 15 SCDs (incidence rate, 17.3; 95% confidence interval [95% CI], 10.5 to 28.8 per 1000 person-years), 35 cardiovascular deaths including SCD (incidence rate, 40.5; 95% CI, 29.1 to 56.4 per 1000 person-years), and 77 all-cause deaths (incidence rate, 89.1; 95% CI, 71.2 to 111.3 per 1000 person-years). The median time to death was 2.4 (interquartile range, 1.4–3.4) years. Among participants with QRS-T angle ≥75°, there were 27 cardiovascular deaths (incidence rate, 55.4; 95% CI, 38.0 to 80.8 per 1000 person-years) and 12 SCDs (incidence rate, 24.6; 95% CI, 14.0 to 43.4 per 1000 person-years); among participants with QRS-T angle <75°, there were eight cardiovascular deaths (incidence rate, 21.2; 95% CI, 10.6 to 42.4 per 1000 person-years) and three SCDs (incidence rate, 7.9; 95% CI, 2.6 to 24.6 per 1000 person-years).
Clinical characteristics of participants stratified by QRS-T angle are compared in Table 1. Overall, the mean left ventricular (LV) ejection fraction was normal. Participants with QRS-T angle ≥75° were more likely to have diabetes, LVH, bundle branch block (BBB), greater change in intradialytic weight, higher systolic and diastolic BPs, and wider unfiltered and filtered QRS. Late potentials were found in 88 (25%) participants who were predominantly older, white men with BBB as well as longer QRS duration (Supplemental Table 1). We also compared the characteristics by primary outcome status (Supplemental Table 2). Those who died were slightly older; had higher body mass index; had lower serum calcium, albumin, and creatinine levels; and had greater mean QRS-T angle.
The unadjusted cumulative incidence of all-cause and cardiovascular mortality was higher in participants with abnormal QRS-T angle (Figure 1, A and B). The cumulative probability of SCD was higher in participants with abnormal QRS-T angle, but this was not statistically significant (Figure 1C). In the unadjusted Cox proportional hazards analysis (Table 2), each 10° increase in spatial QRS-T angle was associated with a 5% increased risk of all-cause death. Even after adjusting for demographic, clinical, and cardiovascular risk factors, every 10° increase in spatial QRS-T was consistently associated with a 5% increase in the risk of all-cause mortality. Higher QRS-T angle was also significantly associated with cardiovascular mortality. The association between continuous QRS-T angle and SCD was similar in direction and magnitude but was not statistically significant. When examined as a categorical variable, spatial QRS-T angle ≥75° was associated with a greater than twofold increase in the risk of all-cause mortality compared with spatial QRS-T angle <75° in both the unadjusted and fully adjusted models (Table 2). Spatial QRS-T angle ≥75° was also associated with a significantly increased risk of cardiovascular mortality and SCD. To confirm that there was a threshold effect, spatial QRS-T angle was divided into tertiles. Compared with the lowest tertile (<64.7°), both second (64.7°–111.5°; adjusted hazard ratio, 2.27; 95% CI, 0.90 to 5.68) and third tertiles (>115°; adjusted hazard ratio, 2.37; 95% CI, 0.93 to 6.09) had greater risk of cardiovascular mortality with similar magnitude of effect, but the associations were not statistically significant. The associations between QRS-T angle as a continuous variable and the outcomes were also consistent when deaths within 1 year after study enrollment were examined (Supplemental Table 3). Late potentials on SAECG were not associated with all-cause, cardiovascular, or sudden cardiac mortality in both unadjusted and adjusted Cox regression analyses (Supplemental Table 4).
Mortality among participants on incident hemodialysis with QRS-T angle ≥75° or <75°. Cumulative (A) all–cause mortality, (B) cardiovascular mortality, and (C) sudden cardiac death.
Risk for all-cause mortality, cardiovascular mortality, and SCD associated with QRS-T angle among 358 patients on incident hemodialysis
The associations of spatial QRS-T angle with all outcomes were not significantly modified by BBB. Associations with the risk of all-cause and cardiovascular mortality and SCD were robust, even when BBB was excluded from the analysis (Supplemental Table 5). Similarly, when patients with QRS duration ≥120 ms were excluded, spatial QRS-T angle remained significantly associated with the risk of all-cause mortality, cardiovascular mortality, and SCD (Supplemental Table 6). Adjusting for each of these additional factors, such as QRS duration, serum electrolytes (calcium, potassium, and magnesium), dialysate electrolytes (calcium and potassium), intradialytic weight change, and albumin, in the main model did not attenuate the association of spatial QRS-T angle with mortality (Supplemental Table 7).
Other traditional ECG measures had modest association with mortality (Supplemental Table 8). Every 10-ms increase in QRS interval was associated with a lower risk of all-cause and cardiovascular mortality. Conversely, heart rate and LVH as detected by ECG were not associated with all-cause mortality, cardiovascular mortality, or SCD.
Discussion
In this prospective cohort study of participants on incident hemodialysis, spatial QRS-T angle, a measure of the electrophysiologic substrate, is an independent risk factor of all-cause mortality, cardiovascular mortality, and importantly, SCD. Spatial QRS-T angle characterizes abnormal ventricular conduction and global electrical heterogeneity, which ultimately predisposes to ventricular arrhythmias and SCD.9,12,13 The independent and robust association of QRS-T angle with cardiovascular mortality and SCD supports the hypothesis that arrhythmias are an important mechanism leading to death after hemodialysis initiation.5,14,15
This compelling association occurred in a young dialysis cohort with a low baseline prevalence of coronary artery disease and also, preserved LV function. This suggests that underlying defects in depolarization and secondary abnormalities in repolarization occur independently of aging, coronary, or LV disease and place patients on dialysis at risk for lethal arrhythmias. This ECG-derived metric, spatial QRS-T angle, may provide evidence of subclinical myocardial fibrosis independent of structural and functional abnormalities typically documented by ECG.16–18 Conversely, other electrophysiologic markers, such as ventricular late potentials by SAECG, did not associate with mortality; however, increased risk may be only evident in those with subepicardial lesions in the late activated regions of the myocardium at the very end of QRS.18
Spatial QRS-T angle is a well established measure of global electrical heterogeneity. In the general population, it is associated with risk for total mortality, cardiovascular death, and SCD in those with myocardial infarction and cardiomyopathy, whereas among patients on prevalent dialysis, it is associated with heart failure.9,12,19,20 We now extend the findings to an incident hemodialysis population with SCD. A previous report studied patients on hemodialysis with ECGs performed up to 18 months after dialysis initiation. Notwithstanding the potential survivor bias, the study showed similar magnitude and direction of risk for arrhythmic-related events.12 Unlike our study, this prior analysis could not account for clinical and cardiovascular characteristics known to significantly affect mortality in dialysis, especially LV function. Moreover, we were able to show, in those with preserved ejection fraction at baseline, that the QRS-T angle is independently associated with mortality and not just heart failure. The results also show a consistently higher risk of spatial QRS-T angle with SCD and cardiovascular and all-cause mortality, even after adjustment of strong risk factors for mortality in dialysis, such as albumin.
Interestingly, in our study, the best threshold of QRS-T angle is lower than previously used in other studies but similar to that in the Multi-Ethnic Study of Atherosclerosis, which yielded a similar graded risk above the 95th percentile of the QRS-T angle.9,17 In a general community population, both QRS-T angle >105° and CKD were associated with increased risk of 1-year mortality by 1.5–2 times, indicating a possible interaction with kidney disease and common activated pathways leading to subclinical endocardial and kidney interstitial fibrosis.16 The specific thresholds for spatial QRS-T angle in a normal healthy population across ethnicities and all ages have not been derived; thus, clinically relevant cut points cannot be provided yet for routine clinical screening. Prior studies used ranges from 40° to >100°, depending on the population and method of derivation. The largest study population of >40,000 individuals free of cardiac disease and with age >55 years suggested that spatial QRS-T angles from 0° to 50° are normal, that spatial QRS-T angles from 50° to 100° are borderline, and that spatial QRS-T angles >100° are abnormal.9,21 Similarly, we examined QRS-T angle in tertiles given the limited power and showed a threshold effect around 75°, suggesting that there is a greater risk of mortality in individuals with QRS-T angle ≥75°.
This study did not find an association of late potentials on SAECG with mortality. Several prior studies reported intradialytic differences in SAECG, such as worsening in low-amplitude signals (<40 mV) in the terminal portion (LAS40) and prolongation of filtered QRS after dialysis.22,23 In our study, late potentials were found in 24% of participants, similar to previous reports of smaller populations.23 Although ventricular late potentials have a strong mechanistic association with ventricular arrhythmias, the SAECG method used may not detect areas of slow conduction occurring early on within the QRS complex.24,25 In addition, the delayed LV epicardial activation caused by interventricular conduction abnormality (which could manifest by wide QRS-T angle) can mask late potentials by SAECG.24 Because abnormal SAECG (late potentials on SAECG) characterizes the presence of a subepicardial scar, the lack of association between late potentials on SAECG and SCD suggests that persons with ESRD may have other arrhythmogenic properties of the myocardium that predispose them to SCD. Even after exclusion of participants with BBB and QRS≥120 ms, however, late potentials remained nonsignificant. The association of late potentials with myocardial fibrosis perhaps indicates late–stage cardiac disease, which was not evident in our incident dialysis population.
Consistent with previous reports, LVH was frequently observed in this study at initiation of dialysis.26 After adjusting for LV mass, there was no attenuation in risk of spatial QRS-T angle and mortality. This suggests that LV structural changes, such as increased LV mass, may be uncoupled from depolarization activity.27 Currently, an ECG is suggested at the time of dialysis by Kidney Disease: Improving Global Outcomes, but it may be warranted to include ECG metrics, such as the spatial QRS-T angle derived from a routine 12–lead ECG.28 Other ECG metrics derived from a 12-lead ECG have also been shown to be predictive of cardiovascular mortality in those with CKD before starting hemodialysis; however, we report no significant associations of the standard ECG parameters with mortality, and the magnitude of effect is not as great as with spatial QRS-T angle.29 The role of ECG monitoring has not been a central focus in the clinical management of patients on dialysis, and therefore, the inclusion of the spatial QRS-T angle and its effect on clinical decision making need to be considered.
There are several limitations to this study. Not all participants completed the cardiovascular study visit, because they were hospitalized, had events, or were too overwhelmed with initiating dialysis. Although the statistical power of the study was sufficient to determine independent association or tested exposure variables, it was not sufficient for many subgroup analyses. We did observe that there was a borderline association of ECG LVH and QRS duration with SCD with similar magnitude and direction as in past studies, but it was not statistically significant, most likely because of limited sample size. Independent validation of our findings is also required in other incident hemodialysis populations. Despite the limitations, important strengths need to be acknowledged. This is a large prospective cohort study of patients on incident hemodialysis with cardiovascular and clinical characterization and a large proportion of blacks. The cohort also has comprehensive baseline cardiovascular phenotyping, including a detailed assessment of electrophysiologic substrate. Importantly, we were able to investigate the cause of death by adjudication and address potential confounding by important factors, such as dialysis–related treatment parameters.
A straightforward measure of the electrophysiologic substrate using a 12-lead ECG, the spatial QRS-T angle, is associated with an increased risk of SCD within a few years after starting hemodialysis and independent of traditional, dialysis–related, and clinical risk factors. Future studies are needed to determine if interventions can mitigate this risk. Spatial QRS-T angle could be included in a future risk scores to guide management of patients on hemodialysis.
Concise Methods
Study Population and Data Collection
The PACE Study is a prospective cohort study of patients on incident hemodialysis in the Baltimore metropolitan area in Maryland. The study protocol was approved by the Johns Hopkins Institutional Review Board. All participants provided written informed consent. Details of the study protocol have been described previously.11 Briefly, ESRD patients were eligible for enrollment if they (1) were 18 years of age or older, (2) had started hemodialysis within 6 months of enrollment, and (3) were English speakers. Exclusion criteria were (1) diagnosed active cancer other than nonmelanoma skin cancer, (2) implanted pacemaker or implantable cardioverter defibrillator, and (3) pregnancy or nursing. In addition, in this analysis, we excluded study participants without available baseline SAECG. Participants underwent a baseline cardiovascular study visit on a nondialysis day within 6 months of initiation of hemodialysis. Detailed medical history, family history, and dialysis laboratories were collected as previously described.11 Average levels of serum albumin, creatinine, calcium, magnesium, and potassium were calculated using values from routine dialysis laboratory blood work during 90 days before the study clinic visit. Intradialytic weight change was calculated as the difference between postdialysis weight and predialysis weight and averaged using values during 90 days before the study visit. Prevalent coronary artery disease, hypertension, and diabetes were determined by physician abstraction from medical records and self-report questionnaires. Antihypertensive medications were assessed at each study visit. Dialysate calcium and potassium values closest to the study visit were recorded.
Echocardiographic examination was performed using a Vivid 7 Cardiac Ultrasound Machine (GE Healthcare, Waukesha, WI). Conventional measurements included LV end diastolic and end systolic diameters and LV ejection fraction by the biplane Simpson method.30 LV mass was estimated as recommended by the American Society of Echocardiography.30 LV mass index normalized by body surface area and height2.7 was calculated. Body surface area was calculated with the Mosteller equation using reported estimated dry weight. Echocardiographic LVH was defined as LV mass index >51 g/m2.7 in men and >47 g/m2.7 in women.31,32
SAECG Recording
SAECG (MAC 5000 with Hi-Res Module; GE Healthcare) signal was recorded using standard x, y, and z orthogonal leads placement at rest with 1000-Hz sampling rate, band pass filtering at 40–250 Hz, and average of 350 QRS complexes as an averaging target. The recorded SAECG signal was analyzed by two different approaches. First, SAECG analysis of late potentials was performed by a commercial Hi-Res SAECG Module (GE Healthcare). Second, the raw unfiltered averaged xyz ECG signal was exported via the ECG Research Workstation, HiRes Module (GE Healthcare), and mean spatial QRS-T angle was measured by the software developed in MATLAB (MathWorks, Inc., Natick, MA).
SAECG Analyses of Late Potentials
The following SAECG characteristics of the filtered QRS were evaluated: (1) total filtered QRS duration, (2) LAS40, and (3) root-mean-square voltage of the last 40 ms.33 Ventricular late potentials were considered positive when two or more of the following criteria were fulfilled: (1) filtered QRS duration>114 ms, (2) LAS40>38 ms, and (3) root-mean-square voltage of the last 40 ms <20 μV.
Measurement of Spatial QRS-T Angle
The unfiltered averaged xyz ECG signal (sampling rate of 1000 Hz; amplitude resolution of 1 μV) was used for the measurement of spatial QRS-T angle. In this study, spatial QRS-T angle was measured as an angle between spatial mean QRS vector and spatial peak T vector (Figure 1). The mean (but not peak) value of the spatial QRS vector over time during the QRS loop was used to characterize the spatial QRS vector as previously described.9 Peak spatial T vector was defined as a vector between the origin point of QRS and T loops and the farthest point in the T loop.34 We chose to use peak (but not mean) T vector, because the end of the T wave (and T loop) on the raw SAECG signal was truncated in patients with very long cardiac cycles. We chose to use mean (but not peak) QRS vector, because QRS loops in the study participants were asymmetric. Spatial QRS-T angle was calculated using the definition of the normalized inner product of spatial mean QRS vector and spatial peak T vector. Accuracy of an automated placement of fiducial points was reviewed by two investigators (E.G. and M.K.).
Definition of Complete BBB
Orthogonal xyz ECGs and Dower–derived 12–lead ECGs were reviewed by two investigators (L.G.T. and A.O.), who evaluated presence of complete left or right BBB according to the established criteria.35
Primary Outcome
The PACE Study participants were followed prospectively with yearly clinic visits. In addition, the research coordinator would conduct phone interviews every 6 months and visit the dialysis units frequently to check with patients and clinic staff.11 Mortality was confirmed by chart review, death certificate, or Centers for Medicare & Medicaid Services (CMS) Death Notification Form 2746. All records from family and dialysis unit interviews; hospital, emergency room, and ambulance records; death certificates; and the CMS-2746 Death Notification forms were reviewed by two trained physician abstractors, and the final determination of cause of death was conducted by the chair of the end point committee or delegate. The chair decided any discordant reviews after committee discussion. We defined SCD as a sudden collapse presumed to be caused by an arrhythmia occurring out of the hospital or in the emergency room witnessed or unwitnessed since the last dialysis treatment. Cardiovascular mortality was defined as the composite of SCD and death caused by arrhythmias, cardiac arrest, ischemic cardiovascular disease, and ischemic cerebrovascular disease either in or outside the hospital.
Statistical Analyses
Continuous normally distributed variables were expressed as means±SDs and compared using t tests. Categorical variables were expressed as frequency (%) and compared by Pearson chi–squared tests. Receiver operating characteristic curves were generated for spatial QRS-T angle. Sensitivity and specificity of each value of spatial QRS-T angle were evaluated. The optimal cut point value of spatial QRS-T angle was determined by several methods: (1) the method by Youden,36 which maximizes the sum of sensitivity and specificity; (2) the method by Liu,37 which maximizes the product of sensitivity and specificity; and (3) the nearest to (0,1) method, which finds the cut point on the receiver operating characteristic curve with perfect sensitivity and specificity. The cut point was bootstrapped (500 replications) to estimate 95% CIs. All three methods (the methods by Youden36 and Liu37 and the nearest to [0,1] method) identified a single optimal threshold value for the spatial QRS-T angle at 74.4° (Supplemental Figure 1), with a bootstrapped 95% CI from 62.6° to 85.8°. Sensitivity at the identified optimal threshold was satisfactory (73%), whereas specificity was low (47%). We rounded to the nearest integer value and dichotomized the exposure into those with abnormal wide spatial QRS-T angle (≥75°) and those with normal spatial QRS-T angle (<75°). Past studies have used other QRS-T angle thresholds on the basis of different methods.12,19,38–40 Additionally, we examined the association of QRS-T angle in tertiles with cardiovascular mortality.
Unadjusted Kaplan–Meier survival analysis was used to compare time to event probabilities in patients with spatial QRS-T angle above the threshold. The log rank statistic was computed to test the equality of survival distributions.
Multivariable Cox proportional hazards regression analyses were performed for all evaluated predictors. A forward model–building approach involving variable selection on the basis of changes in effect size, known associations with mortality in patients on dialysis, and significant P values was used to build the final multivariable model. Cox regression analyses were stratified by left or right BBB. In addition, sensitivity analyses were conducted excluding participants with interventricular conduction delay (QRS≥120 ms) and BBB. Finally, we tested for an interaction between QRS-T angle and BBB using the likelihood ratio test.41,42
Proportional hazards assumption was tested on the basis of Schoenfield residuals for each Cox model. To avoid overfitting, several Cox models were generated with adjustment of demographic, cardiovascular, dialysis-related, and inflammatory factors.5 We checked for nonlinear association between mean spatial QRS-T angle as a continuous variable and all-cause mortality using restricted cubic splines. All missing data were imputed using the multiple imputation by chained equations method.43 All statistical analyses were performed in Stata 14.1 (StataCorp., College Station, TX). A two–tailed P value <0.05 was considered statistically significant.
Disclosures
None.
Acknowledgments
We dedicate the Predictors of Arrhythmic and Cardiovascular Risk in End Stage Renal Disease (PACE) Study to our friend and colleague W.H.L.K. We thank the PACE Study staff and participants and the nephrologists and staff of the DaVita and MedStar dialysis units in the Baltimore area. We also thank Dr. Joel Xue (GE Healthcare, Waukesha, WI) for providing GE ECG Research Utility software for analysis.
The PACE Study was supported by National Institute of Diabetes and Digestive and Kidney Diseases grant R01DK072367 (to R.S.P.). This study was partially supported by National Heart, Lung, and Blood Institute grant R01HL118277 (to L.G.T.). This study was partially supported by Boston Scientific as an investigator-initiated research project (L.G.T.).
Some of the data in this article were presented at the American Health Association Scientific Sessions held November 7–11, 2015 in Orlando, FL.
Footnotes
Published online ahead of print. Publication date available at www.jasn.org.
This article contains supplemental material online at http://jasn.asnjournals.org/lookup/suppl/doi:10.1681/ASN.2015080916/-/DCSupplemental.
- Copyright © 2016 by the American Society of Nephrology