Dynamic predictive accuracy of electrocardiographic biomarkers of sudden cardiac death within a survival framework: the Atherosclerosis Risk in Communities (ARIC) study

Background

The risk of sudden cardiac death (SCD) is known to be dynamic. However, the accuracy of a dynamic SCD prediction is unknown. We aimed to measure the dynamic predictive accuracy of ECG biomarkers of SCD and competing non-sudden cardiac death (non-SCD).

Methods

Atherosclerosis Risk In Community study participants with analyzable ECGs in sinus rhythm were included (n = 15,716; 55% female, 73% white, age 54.2 ± 5.8 y). ECGs of 5 follow-up visits were analyzed. Global electrical heterogeneity and traditional ECG metrics (heart rate, QRS, QTc) were measured. Adjudicated SCD was the primary outcome; non-SCD was the competing outcome. Time-dependent area under the receiver operating characteristic curve (ROC(t) AUC) analysis was performed to assess the prediction accuracy of a continuous biomarker in a period of 3,6,9 months, and 1,2,3,5,10, and 15 years using a survival analysis framework. Reclassification improvement as compared to clinical risk factors (age, sex, race, diabetes, hypertension, coronary heart disease, stroke) was measured.

Results

Over a median 24.4 y follow-up, there were 577 SCDs (incidence 1.76 (95%CI 1.63–1.91)/1000 person-years), and 829 non-SCDs [2.55 (95%CI 2.37–2.71)]. No ECG biomarkers predicted SCD within 3 months after ECG recording. Within 6 months, spatial ventricular gradient (SVG) elevation predicted SCD (AUC 0.706; 95%CI 0.526–0.886), but not a non-SCD (AUC 0.527; 95%CI 0.303–0.75). SVG elevation more accurately predicted SCD if the ECG was recorded 6 months before SCD (AUC 0.706; 95%CI 0.526–0.886) than 2 years before SCD (AUC 0.608; 95%CI 0.515–0.701). Within the first 3 months after ECG recording, only SVG azimuth improved reclassification of the risk beyond clinical risk factors: 18% of SCD events were reclassified from low or intermediate risk to a high-risk category. QRS-T angle was the strongest long-term predictor of SCD (AUC 0.710; 95%CI 0.668–0.753 for ECG recorded within 10 years before SCD).

Conclusion

Short-term and long-term predictive accuracy of ECG biomarkers of SCD differed, reflecting differences in transient vs. persistent SCD substrates. The dynamic predictive accuracy of ECG biomarkers should be considered for competing SCD risk scores. The distinction between markers predicting short-term and long-term events may represent the difference between markers heralding SCD (triggers or transient substrates) versus markers identifying persistent substrate.

Sudden cardiac death (SCD) is a major contributor to cardiovascular mortality, accounting for 40–50% of the years of potential life lost from all cardiovascular diseases (CVD) [1, 2]. In the United States (US), more than 350,000 emergency medical services-assessed out-of-hospital sudden cardiac arrests occur annually [1]. There remains a lack of reliable, dynamic predictors of SCD [3]. An electrocardiogram (ECG) can characterize the presence and properties of the electrophysiological substrate of SCD. Our group recently showed that global electrical heterogeneity (GEH), as measured by five metrics [spatial QRS-T angle, spatial ventricular gradient (SVG) azimuth, elevation, and magnitude, and sum absolute QRST integral (SAI QRST)] is independently (after comprehensive adjustment for time-updated CVD events and their risk factors) associated with SCD, representing an underlying substrate of SCD [4]. The subsequent discovery of 10 genetic loci, associated with GEH at a genome-wide significance level, confirmed the presence of several underlying mechanisms behind the GEH ECG phenotype [5]. We also developed a competing risk score of SCD and showed that the addition of GEH measures to clinical risk factors significantly improves the reclassification of SCD risk [4]. The risk of SCD is known to be dynamic. However, current risk models predict SCD using baseline risk factors measured at a single point in time. Therefore, evaluating the accuracy of a dynamic prediction is, therefore, necessary to better understand the temporal relationship between substrate and events. The goal of this study was to investigate the dynamic predictive accuracy of GEH and traditional ECG biomarkers of SCD within a survival framework in comparison with competing non-sudden cardiac death (non-SCD) in the Atherosclerosis Risk in Community (ARIC) study participants.

Study population

The ARIC study is an ongoing prospective cohort study evaluating risk factors, progression, and outcomes of atherosclerosis in 15,792 participants (45% male, 74% white) enrolled in four US communities in 1987–1989. The ARIC study protocol and design have been previously described [6]. We excluded ARIC participants with absent or poor-quality ECGs due to noise, artifacts, or missing leads (n = 24), atrial fibrillation (AF) (n = 36), and ventricular pacing (n = 16) at the baseline study visit. Only participants in normal sinus rhythm were included in this study (n = 15,716). If AF or ventricular pacing were diagnosed at any time during follow-up, such participants were included for a period until AF or ventricular pacing were diagnosed on 12-lead ECG. Prevalent CVD was defined as the presence of at least one baseline prevalent condition: coronary heart disease (CHD), heart failure (HF), stroke, peripheral artery disease (PAD), atrioventricular (AV) block II-III, atrial or ventricular pacing, or Wolff-Parkinson-White ECG phenotype.

Clinical characteristics of participants

Prevalent CHD was defined as a history of myocardial infarction (MI), or coronary revascularization via coronary artery bypass surgery or percutaneous coronary intervention. Prevalent MI was defined as a self-reported history of MI and/or ECG evidence of MI as defined by the Minnesota code [7]. Prevalent AF was defined as either a self-reported and validated history of AF or a diagnosis of AF on the baseline ECG. Prevalent HF was defined as self-reported use of HF medication or evidence of symptomatic HF as defined by stage 3 of the Gothenburg criteria [8], which required the presence of specific cardiac and pulmonary symptoms in addition to medical treatment of HF. Prevalent stroke in ARIC was defined by a stroke and transient ischemic attack diagnostic algorithm, as previously described [9]. PAD was defined as the ankle-brachial index ≤0.90. Details of ankle-brachial index measurement in the ARIC study have been previously described [10].

Definition of a primary outcome: sudden cardiac death

A follow-up of ARIC participants was previously reported [11]. SCD was defined as a sudden pulseless condition in a previously stable individual without evidence of a non-cardiac cause of cardiac arrest if the cardiac arrest occurred out of the hospital or in the emergency room. To identify cases of SCD in ARIC, cases of fatal CHD that occurred by December 31, 2012 were reviewed and adjudicated by a committee of physicians in two phases, as previously described [12]. CHD deaths occurring on or before December 31, 2001 were adjudicated in the first phase. CHD deaths occurring between January 1, 2002 and December 31, 2012 were adjudicated in the second phase. Available data from death certificates, informant interviews, physician questionnaires, coroner reports, prior medical history, and hospital discharge summaries were reviewed, in addition to circumstances surrounding the event. Each event was adjudicated independently by two physicians. In cases of disagreement, a third reviewer independently reviewed the event to provide final classification. Definite, probable, or possible SCD was included in this study as a primary outcome: the strength of available evidence determined this stratification. Definite SCD included cases of witnessed SCD with definite evidence such as an available rhythm strip of life-threatening cardiac arrhythmia, or primary emergency medical services impression of cardiac arrest. Probable SCD was defined as SCD with uncertainty either due to concomitant clinical conditions that can muddle the exact cause of demise, or limited information to adjudicate an event. Possible SCDs were adjudicated only in the second phase of reviews. The strength of evidence for probable SCD was greater than the strength of evidence for possible SCD. Possible SCD included cases of death that were unwitnessed but had specified SCD on a death certificate and did not document another cause of death. Participants were censored at the time of loss to follow-up, incident AF or ventricular pacing on 12-lead ECG, or death if the cause of death was any other than SCD.

Competing mortality outcome: non-sudden cardiac death

Cases of fatal CHD were adjudicated by the ARIC Morbidity and Mortality Classification Committee, as previously described [11]. Fatal CHD that did not meet the criteria of SCD comprised the non-SCD outcome.

12 lead ECG recording

12-lead ECG was recorded according to the ARIC study protocol and manual (version 1.0; August 1987). The method and procedure for 12-lead ECG recording, as described in the ARIC manual, are outlined below. During the baseline examination, a standard supine 12-lead ECG was recorded after a 12-h fast followed by a light snack and at least 1 h after smoking or ingestion of caffeine. The standard electrocardiograph for the ARIC study was the MAC PC by Marquette Electronics, Inc. A 12-lead resting ECG was obtained consisting of 10 s of each of the leads (I, II, III, aVR, aVL, aVF, Vl-V6) simultaneously recorded. In an effort to enable longitudinal comparisons of ECGs, ARIC investigators developed and implemented a uniform procedure for electrode placement and skin preparation. The participant, stripped to the waist, was instructed to lie on the recording bed with arms relaxed at the sides. The individual was asked to avoid movements that may cause errors in marking the electrode locations. For optimal electrode/skin interface, the electrodes were placed on the skin at least 2–3 min before recording the ECG. A pen was used to mark the six chest electrode positions. The chest was wiped with a sterile alcohol prep. Left leg electrode was placed on the medial surface of the left ankle. Right leg electrode was placed on the medial surface of the right ankle. Left arm electrode was placed on the medial surface of the left wrist. Right arm electrode was placed on the medial surface of the right wrist. Electrode V1 was located in the 4th intercostal space at the right sternal border, immediately to the right of the sternum. Electrode V2 was located in the 4th intercostal space, immediately to the left of the sternal border. Next, E-point was located by finding the 5th intercostal space, and following it horizontally to the midsternal line. Location of V6 electrode was found using the chest square. V6 was located at the same level as the E point in the midaxillary line (straight down from the center of the armpit). If breast tissue was over the V6 area, V6 was placed on top of the breast. No attempt was made to move the breast. Electrode V4 was located using a flexible ruler, as a midway between E and V6. Electrode V3 was located using a flexible ruler, midway between V2 and V4. Using a flexible ruler, electrode V5 was located midway between the locations of V4 and V6. After placing the electrodes on the skin, the participant’s information was input to the MAC PC. It was required that electrodes must be on the skin for at least 3 min before taking the ECG. During the ECG recording, special attention was paid to the quality of recording. Quality control and technical troubleshooting procedures were in place to minimize errors (lead switch), noise, and artifacts. Recorded 12-lead ECGs were originally saved in the memory of the ECG machine and transmitted to a MUSE database (GE Marquette, Milwaukee, WI) within the Halifax ECG Computing Center via the phone line. Later, the MUSE database was transferred from the Halifax ECG Computing Center to the Epidemiological Cardiology Research Center (EPICARE, Wake Forest University, NC), and then to the Tereshchenko laboratory at Oregon Health & Science University.

Electrocardiogram analyses

ECG data from all five follow-up visits were analyzed. Traditional ECG intervals were reported by the 12 SL algorithm using Magellan ECG Research Workstation V2 (GE Marquette Electronics, Milwaukee, WI). QT interval was corrected for heart rate according to Bazett’s formula. We analyzed raw, digital, 10-s, 12-lead ECGs (sampling rate of 500 Hz and amplitude resolution of 1 μV). Origin and conduction path of each cardiac beat was adjudicated by the team of physicians (DG, AB, SVM, LGT), and each beat was manually labeled by investigators (CH, JAT) for subsequent automated analyses. A representative normal sinus median beat was constructed. For the development of a normal sinus median beat, sinus beats before and after premature ventricular complexes, and noisy or distorted beats were excluded. In this study, only the normal sinus median beat was used for our analysis. GEH was measured as previously described, [4, 13] in a time-coherent median beat with the identified origin of the heart vector [14]. We have provided the open-source software code at Physionet (https://physionet.org/physiotools/geh/). In addition to previously reported “mean” GEH measures, [4] in this study we measured the spatial peak vectors (Fig. 1) [13,14,15]. First, we transformed the 12-lead ECG into an orthogonal XYZ ECG, using Kors transformation [16]. Next, we constructed a time-coherent median beat, and detected the origin of the heart vector using our novel approach, as recently described [14]. Then, we performed calculations of GEH metrics using the following equations.

GEH measurement example of peak and area vectors, and vector magnitude (VM). a (i) VM plotted over time, and (ii) corresponding X, Y, and Z leads. b VM plotted in a three-dimensional space. Color-coded progression from QRS onset (red) to the end of T (purple) is shown. c Measurement of (i) QT integral on VM, and (ii) SAI QRST. d. Measurement of peak SVG vector magnitude, azimuth, and elevation

Full size image

Spatial QRS-T angles

Spatial peak QRS-T angle was calculated as the 3-dimensional angle between the spatial peak QRS vector and the spatial peak T vector:

Spatial area QRS-T angle was calculated as the 3-dimensional angle between the spatial area QRS vector and the spatial area T vector:

Spatial ventricular gradient vectors:

Magnitude and direction of spatial area (Wilson’s) and peak SVG vectors were measured.

Wilson’s (area) SVG magnitude was also calculated:

Statistical analysis

Time-dependent area under the receiver operating characteristic curve (ROC(t)AUC) analysis was performed to assess the predictive accuracy of a continuous biomarker in a period of 3, 6, 9 months, and 1,2,3,5,10, and 15 years, using an unadjusted survival analysis framework approach [17, 18]. We used the nearest neighbor estimator, which allows the censoring to depend on the marker and is therefore realistic. The percentage of observations included in each neighborhood was defined by the eq. 0.25* \( \Big(\sqrt[3]{n} \)), where n is the number of observations. All available five visits’ ECG data were included in time-dependent AUC analysis [4]. To satisfy the requirement for ROC analysis, and to perform internal validation of study findings, we divided the dataset into five unique (non-overlapping) partitions. In each partition, a study participant was presented not more than once, with a unique time to event defined as time from ECG recording to the time of outcome (or censoring). If a participant had five ECGs recorded at five study visits, each ECG contributed to a different partition. Those participants who had less than 5 visits/ECGs were randomly distributed across 5 partitions. Bootstrapping with 500 replications was performed to determine a 95% confidence interval (CI) of ROC(t) AUC in each partition, separately. Then, ROC(t) AUC point estimates, and lower and upper boundaries of 95%CI were averaged across 5 partitions, then presented as a final summary result. To assess the statistical power of ROC(t) AUC analysis at each time period, we compared observed 95% CI width with expected 95% CI width. Expected 95% CI width was calculated for observed sample size and observed AUC values at each time period. We summarized clinical characteristics of study participants with an SCD outcome within the first 3 months, 3–6 months, 6 months-1 year, 1–2 years, 2–5 years, and more than 5 years after ECG recording in a longitudinal dataset, reporting between-participant standard deviation (SD) for continuous variables, and between-participant frequencies for categorical variables. We then assessed whether the addition of traditional ECG metrics (heart rate, QRS, QTc) and GEH metrics to our previously identified clinical risk factors of SCD [4] (age, sex, race, diabetes, hypertension, CHD, and stroke) resulted in better predictive accuracy for SCD and non-SCD within the first 3 months, 3–6 months, 6 months-1 year, 1–2 years, 2–5 years, and more than 5 years after ECG recording. We calculated absolute integrated discrimination improvement (IDI), and net reclassification improvement (NRI) using multivariate logistic regression [19, 20]. IDI estimates improvement in average sensitivity and specificity. We estimated category-free NRI and two-category NRI for events, defining the high-risk category as a ≥ 25% risk of SCD/non-SCD within the first 3 months, 3–6 months, and 6 months-1 year after ECG recording. The high-risk category for events occurring 1–2 years, 2–5 years, and more than 5 years after ECG recording was defined as ≥10% risk of SCD/non-SCD. Statistical analysis was performed using STATA MP 15.1 (StataCorp LP, College Station, TX, USA). A P-value of < 0.05 was considered significant. PASS 2019 Power Analysis and Sample Size Software (NCSS, LLC. Kaysville, Utah, USA) was used for the calculations of the expected 95% CI width.

Study population

The clinical characteristics of the study population are shown in Table 1. Approximately half of the study participants were female, and 73% were white. Average traditional ECG parameters were normal. During the SCD adjudication, the inter-reviewer agreement was 83.2%, and agreement across phases was 92.5%. Over a median follow-up of 24.4 years, there were 577 SCDs (incidence 1.76 (95%CI 1.63–1.91) per 1000 person-years), and 829 non-SCDs (incidence 2.54 (95%CI 2.37–2.71) per 1000 person-years). SCD victims who died within the first 3 months after ECG recording were more likely to be CVD-free white males with fewer prevalent CVD risk factors. In contrast, SCD victims who died more than 5 years after ECG recording had nearly equal probabilities of being male or female, white or non-white (Table 1).

Time-dependent AUC analyses of SCD and competing outcome

Detailed results for all 5 partitions are provided in Additional file 1 Table S1 and Figs. 2, 3, 4. Considering the results of robust validation analysis, no ECG biomarkers predicted SCD within 3 months after ECG recording. Forward-directed SVG azimuth was significant predictor in some, but not all partitions: peak SVG azimuth partition #1 AUC 0.828 (95%CI 0.709–0.947), and partition #2 AUC 0.727 (95%CI 0.712–0.742). Predicted 95%CI width (Table 2) indicated insufficient statistical power for prediction of SCD within 3 months after ECG recording for most ECG biomarkers. At 6 months, SVG elevation was the only short-term statistically significant biomarker of SCD (area SVG elevation AUC 0.706; 95%CI 0.526–0.886), which specifically predicted SCD, but did not predict non-SCD (Fig. 6c). Figure 6 illustrates trends in variations in SVG vector direction (azimuth and elevation) predicting SCD versus non-SCD within 9 months – 1 year after ECG recording, when upward, and more likely forward-directed SVG vector predicted SCD, whereas leaning backward-directed SVG predicted non-SCD. In contrast, long-term (2 years and beyond) predictors of SCD and non-SCD had many similarities. Wide spatial QRS-T angle and backward [towards left ventricle (LV)] – directed SVG vector predicted both SCD and non-SCD. The most accurate prediction of SCD was provided by ECG biomarkers recorded within 2 years before the outcome (Figs. 5, 6). Most ECG biomarkers predicted both mortality outcomes, with few exceptions. SVG magnitude did not predict any outcome. At any time, heart rate and QTc were stronger predictors of non-SCD than SCD (Fig. 5). Neither heart rate nor QTc improved risk reclassification beyond traditional clinical risk factors of SCD (Table 3), whereas QRS duration improved reclassification of SCD risk for events occurring 2–5 years after ECG recording. Within the first 3 months after ECG recording, only SVG azimuth improved reclassification of the risk beyond traditional clinical risk factors (18% SCD events were reclassified from a low or intermediate-risk category to a high-risk category). All GEH metrics significantly improved reclassification of both SCD and non-SCD beyond clinical risk factors for events occurring at least 1 year after ECG recording (Table 3).

Internal validation of time-dependent AUC with 95% CI for prediction of SCD, for windows of prediction 3, 6, 9 months, 1, 2, 3, 5, 10, 15 years for heart rate, QTc, QRS duration, and SAI QRST in five partitions of the study dataset

Full size image

Internal validation of time-dependent AUC with 95% CI for prediction of SCD, for windows of prediction 3, 6, 9 months, 1, 2, 3, 5, 10, 15 years for spatial peak and area QRS-T angle and SVG magnitude in five partitions of the study dataset

Full size image

Internal validation of time-dependent AUC with 95% CI for prediction of SCD, for windows of prediction 3, 6, 9 months, 1, 2, 3, 5, 10, 15 years for peak and area SVG azimuth and elevation in five partitions of the study dataset

Full size image

Time-dependent AUC with 95% CI for prediction of SCD (red diamonds), and non-SCD (blue circles) for windows of prediction 3, 6, 9 months, 1, 2, 3, 5, 10, 15 years for (a) heart rate, (b) QTc interval, (c) QRS duration, (d) SAI QRST measured at visits 1, 2, 3, 4, and 5

Full size image

Time-dependent AUC with 95% CI for prediction of SCD (red diamonds), and non-SCD (blue circles) for windows of prediction 3, 6, 9 months, 1, 2, 3, 5, 10, 15 years for (a) spatial peak QRS-T angle, (b) area SVG azimuth, (c) area SVG elevation, (d) area SVG magnitude measured at visits 1, 2, 3, 4, and 5

Full size image

Comparison of area-based versus peak-based SVG and QRS-T angle measurements

SCD predictive accuracy of peak-vector-based versus area-based spatial QRS-T angle, SVG direction, and SVG magnitude was similar (Fig. 7), with a trend towards slightly better prediction by peak-vector-based QRS-T angle and SVG azimuth, as compared to area-based QRS-T angle and SVG azimuth metrics.

Comparison of time-dependent AUC (with 95% CI) for windows of SCD prediction 3, 6, 9 months, 1, 2, 3, 5, 10, 15 years, for peak-based SVG vector measurements (orange squares) vs. area-based SVG vector measurements (green circles): a QRS-T angle, b SVG azimuth, c SVG elevation, d SVG magnitude

Full size image

Internal validation

Results of internal validation are shown in Additional file 1: Table S1 and Figs. 2, 3, 4. As expected, long-term prediction of SCD was consistent across all 5 study partitions, whereas short-term (within 3–9 months) prediction was unreliable for most ECG variables (except SVG elevation). Consistently, statistical power was reliably sufficient for ROC(t) AUC analyses for outcomes occurring at least 1 year after ECG recording, and beyond (Table 2).

In this study, we described the dynamic predictive accuracy of ECG and VCG biomarkers of two competing mortality outcomes: SCD and non-SCD within a survival framework. Within the identified dynamic predictors of SCD, there was a distinction between markers predicting short-term events (within 6 months) and markers predicting more intermediate- and long-term events. This may represent the difference between markers heralding SCD (triggers or transient substrates) versus markers identifying persistent substrate. As expected, transient substrate of non-SCD (describing structural heart disease substrate) was characterized by wide QRS-T angle, SVG vector pointing backward (towards LV), wide QRS, prolonged QTc, and increased heart rate. A transient substrate of SCD was characterized by SVG vector pointing upward (towards the outflow tract). Dynamic predictive accuracy of ECG and VCG biomarkers of SCD should be taken into account for development of dynamic and life-long prediction of SCD and non-SCD. Importantly, the addition of GEH metrics (but not QTc) to known demographic and clinical risk factors (age, sex race, CHD, stroke, diabetes, and hypertension) significantly improved reclassification of risk, supporting inclusion of GEH metrics into dynamic risk scores for SCD.

Triggers, or transient substrate of SCD event

An SCD event represents a “perfect storm,” requiring both susceptible anatomical/functional substrate and a trigger/transient initiating event [21]. Short-term predictors of SCD in our study reflect possible SCD triggers. The SVG vector direction predicting short-term SCD risk differed from SVG vector direction of the intermediate- and long-term risk of SCD. The short-term risk was uniquely predicted by an SVG vector pointing upward (toward the outflow tract), suggesting that short total recovery time in the outflow tracts (as the SVG vector points towards an area with the shortest total refractory time) [22] may represent an SCD trigger. SVG azimuth was the only ECG metric which improved reclassification of the risk beyond known clinical and demographic risk factors. Indeed, it is known that “malignant” idiopathic ventricular fibrillation and polymorphic ventricular tachycardia can be triggered by ventricular ectopy arising from the outflow tracts [23]. Early cardiac development affects the generation of electrophysiological heterogeneities in the adult heart [24]. There may be a genetic basis for this phenomenon as GEH-associated genetic polymorphisms indicated the involvement of HAND1 and TBX3 genes [5], both of which play a role in outflow tract development.

Intermediate and long-term substrates of SCD

All ECG measurements (except SVG magnitude) predicted SCD long-term. The long-term substrate of SCD was characterized by an SVG vector pointing backward (toward the LV), a wide spatial QRS-T angle, and a large SAI QRST. Reliable long-term prediction of SCD offers an opportunity for early preventive intervention. Many GEH-associated genetic loci are implicated in cardiac development [5]. Further studies of the underlying biology behind GEH-associated loci will help to uncover novel mechanisms of SCD and develop primary prevention strategies. A recent case-control genome-wide association study of sudden cardiac arrest [25] did not identify any variants at genome-wide statistical significance. An ideal case-control study of paroxysmal life-threatening arrhythmias (e.g. SCD) would require evidence of freedom from arrhythmogenic substrate in controls, which is difficult to achieve. As both trigger and substrate are required for the development of sudden cardiac arrest, a low yield from a case-control genome-wide association study of sudden cardiac arrest is to be expected. In contrast, genomic studies of electrophysiological substrates have the advantage of a more accurate measurement of phenotype and larger statistical power (as an outcome is a continuous variable), providing higher yield.

Dynamic predictive accuracy of biomarkers within a survival framework

The dynamic nature of SCD risk is well-recognized. However, the dynamic predictive accuracy of SCD risk markers has not been previously studied. Our large prospective epidemiological study used repeated ECG measures, obtained at five follow-up visits, which ensured stable estimates of the dynamic predictive accuracy of ECG biomarkers within a survival framework. An analytical framework for the assessment of a dynamic predictive accuracy of biomarkers for censored survival data was developed fairly recently [17]. Heagerty et al. [17] showed that a simple estimator based on Kaplan-Meier method has serious shortcomings for characterization of accuracy for censored survival outcomes, and developed the nearest neighbor estimator as a valid ROC solution for prediction accuracy assessment, allowing the censoring process to depend on the marker. In this study, we used an analytical approach to answer an agnostic predictive accuracy question. To mimic the real-life clinical scenarios, we intentionally did not adjust for confounders and therefore did not comment on the independence of association of ECG biomarkers with SCD at any given time point. There were noticeable differences in the clinical characteristics of study participants who died suddenly within 3 months after ECG recording, as compared to those who experienced SCD 5 years after ECG recording. Nevertheless, observed dissimilarities in a dynamic predictive accuracy of ECG biomarkers suggested different mechanisms behind short-term SCD triggers (or transient substrates) and long-term SCD substrates. A study of SCD triggers is objectively difficult to conduct. The methodological approach of the dynamic predictive accuracy of ECG biomarkers within the survival framework can provide unique perspective on transient substrates and triggers of SCD, which prompts further investigation.

Global electrical heterogeneity

The GEH concept is based on Wilson’s ventricular gradient [22, 26]. SVG defines a vector that characterizes the magnitude and direction of the steepest gradient between the longest and the shortest total recovery time across the entire heart, both left and right ventricles [27]. SVG vector is directed towards a zone with the shortest total recovery time. Therefore, GEH is a global measure of the dispersion of total recovery time across the heart, a marker of an underlying arrhythmogenic substrate, encompassing dispersion in both activation and refractoriness. We comprehensively characterize GEH by measuring five features of the SVG vector (Fig. 1): SVG magnitude, SVG direction (azimuth and elevation), QRS-T angle, and SVG’s scalar SAI QRST. The spatial QRS-T angle [28] is a well-known marker of the risk of SCD and cardiovascular mortality. SAI QRST is a scalar analog of the SVG [15, 29,30,31], associated with ventricular tachyarrhythmias in heart failure [29, 30, 32]. We previously showed that GEH is associated with SCD after adjustment for traditional ECG metrics and a comprehensive list of cardiovascular risk factors [4]. Five GEH metrics only weekly correlate with each other and traditional ECG measures, [5] and, therefore, are well-suited for inclusion in SCD risk scores. Importantly, GEH genome-wide association study [5] provided first clues about underlying biology behind GEH, which can lead to the future development of novel methods of SCD prevention. In vectorcardiography, there are two major approaches to define spatial vectors: either measuring spatial peak or area vectors [28]. In our study, some peak-based GEH metrics outperformed area-based GEH metrics. This finding may be at least partially explained by the fact that we used a physiologically sound definition of the heart vector origin point and time-coherent global median beat [14], which permitted accurate measurement of peak vectors. GEH metrics are highly reproducible [33], which is an essential factor for their future implementation in clinical practice. Of note, GEH can be calculated using 12-lead ECG that was recorded by any manufacturer ECG machine. The open-source software code for GEH calculation is provided at https://github.com/Tereshchenkolab/Origin and https://www.physionet.org/content/geh.

Clinical implications

Our study is the first step towards the development of dynamic SCD risk prediction. Recently, we developed an SCD risk score [4], available at www.ecgpredictscd.org. However, our SCD risk model, similarly to other risk models [34, 35], includes risk factors that were measured once at baseline. Further improvement of SCD risk stratification requires knowledge of the dynamic predictive accuracy of ECG biomarkers of SCD. After validation of our study findings in another cohort, recommendations regarding an optimal frequency of ECG recordings can be developed.

Strengths and limitations

The strength of our study derives from the large prospective cohort design, with five longitudinal ECG recordings, long-term (median 24 years) follow-up, and a well-adjudicated SCD outcome. However, limitations of the study should be taken into account. The small number of events within 3 and 6 months after ECG recording limited the statistical power of SCD trigger analyses. It is worth noting that 9 out of 11 ARIC participants who succumbed to SCD within 3 months after ECG recording were men. We cannot rule out the possibility that the observed transient substrate of SCD is sex-specific. Further studies of SCD triggers in women are needed. While we performed robust internal validation of our findings using bootstrapping in five partitions of the dataset, replication of the SCD trigger analyses in another prospective cohort is needed. Nevertheless, this is the largest prospective study of SCD triggers and substrates, suggesting differences between long-term and transient SCD-specific and non-SCD-specific substrates [36]. In this study, correlation between GEH metrics and heart rate was weak (r values between 0.1–0.2), and we did not normalize GEH metrics by heart rate. However, further studies are needed to determine whether normalization by heart rate can further improve the predictive value of GEH.

Dynamic predictive accuracy of ECG and VCG biomarkers should be taken into account for the development of dynamic risk scores of competing SCD risk. The distinction between markers predicting short-term and long-term events may represent the difference between markers heralding SCD (triggers or transient substrates) versus markers identifying persistent substrate.

After the ethical approval at a requestor institution, the de-identified ARIC data can be obtained from the National Heart, Lung, and Blood Institute – maintained BioLINCC repository [37], and via ARIC Coordinating Center at the University of North Carolina—Chapel Hill [38]. The BioLINCC website [37] includes detailed information about the process to obtain such data. Details about the procedures for data request from the ARIC Coordinating Center at the University of North Carolina—Chapel Hill can be found online [38]. GEH results were reported to the ARIC Coordinating Center by the ARIC ancillary study “Novel ECG measures and risk of sudden cardiac death” (principal investigator Tereshchenko). We provided open-source MATLAB (MathWorks, Natick, MA, USA) code at: https://github.com/Tereshchenkolab/Global-Electrical-Heterogeneity, and https://github.com/Tereshchenkolab/Origin. SCD risk calculator is available at http://www.ecgpredictscd.org/.

AF:

Atrial fibrillation

ARIC:

Atherosclerosis risk in community study

AUC:

Area under the curve

AV:

Atrioventricular

CHD:

Coronary heart disease

CI:

Confidence interval

CVD:

Cardiovascular disease

ECG:

Electrocardiogram

GEH:

Global electrical heterogeneity

HF:

Heart failure

IDI:

Integrated discrimination improvement

LV:

Left ventricle

MI:

Myocardial infarction

Non-SCD:

Non-sudden cardiac death

NRI:

Net reclassification improvement

PAD:

Peripheral artery disease

ROC:

Receiver operating characteristic

RVOT:

Right ventricular outflow tract

SAI QRST:

Sum absolute QRST integral

SCD:

Sudden cardiac death

SD:

Standard deviation

SVG:

Spatial ventricular gradient

US:

United States

VM:

Vector magnitude

The authors thank the staff and participants of the ARIC study for their important contributions. We would like to acknowledge the SCD mortality classification committee members: Nona Sotoodehnia (lead), Selcuk Adabag, Sunil Agarwal, Lin Chen, Rajat Deo, Leonard Ilkhanoff, Liviu Klein, Saman Nazarian, Ashleigh Owen, Kris Patton, and Larisa Tereshchenko.

Funding

The ARIC study has been funded in whole or in part with Federal funds from the National Heart, Lung, and Blood Institute, National Institutes of Health, Department of Health and Human Services, under Contract nos. (HHSN268201700001I, HHSN268201700002I, HHSN268201700003I, HHSN268201700004I, HHSN268201700005I). This work was supported by HL118277 (LGT). The funding body had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Affiliations

1. Knight Cardiovascular Institute, Department of Medicine, Oregon Health & Science University School of Medicine, 3181 SW Sam Jackson Park Rd; UHN62, Portland, OR, 97239, USA

   • Erick A. Perez-Alday
   • , Aron Bender
   • , David German
   • , Srini V. Mukundan
   • , Christopher Hamilton
   • , Jason A. Thomas
   • , Yin Li-Pershing
   •  & Larisa G. Tereshchenko

2. Department of Biomedical Informatics, Emory University, Atlanta, GA, USA

   • Erick A. Perez-Alday

3. UCLA Cardiac Arrhythmia Center, University of California Los Angeles, Los Angeles, CA, USA

   • Aron Bender

4. Rush University, Chicago, IL, USA

   • Srini V. Mukundan

5. Rosalind Franklin University of Medicine and Science, North Chicago, IL, USA

   • Christopher Hamilton

6. University of Washington, Seattle, WA, USA

   • Jason A. Thomas

7. Cardiovascular Division, Department of Medicine, Johns Hopkins University School of Medicine, Baltimore, MD, USA

   • Larisa G. Tereshchenko

Contributions

EAPA and LGT designed the study. AB, DG, SVM and LGT reviewed 12-lead ECGs and adjudicated origin and conduction of each cardiac beat. EAPA and YLP developed MATLAB software code for automated ECG analyses and analyzed digital ECG signal. CH and JT checked quality of automated ECG analysis and verified accuracy of fiducial points detection. LGT performed the statistical analyses. EAPA and LGT wrote the first draft of the manuscript. All authors participated in data interpretation, discussion and preparation of the final manuscript, and read and approved the final manuscript.

Corresponding author

Correspondence to Larisa G. Tereshchenko.

Ethics approval and consent to participate

This study was approved by the Oregon Health & Science University and the Johns Hopkins University Institutional Review boards (IRB). The ARIC study protocol was approved by the IRB at all participating institutions: University of North Carolina at Chapel Hill IRB, Johns Hopkins University IRB, University of Mississippi Medical Center IRB, and University of Minnesota IRB. All participants signed informed consent before entering the study. The study protocol conforms to the ethical guidelines of the 1975 Declaration of Helsinki.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Reprints and Permissions