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[5].
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). Only participants in normal sinus rhythm were included in this study (n=15,716).
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), angina pectoris, 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[6]. Prevalent AF was defined as either a self-reported and validated history of AF or 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,[7] 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[8]. PAD was defined as the ankle-brachial index ≤0.90. Details of ankle-brachial index measurement in the ARIC study have been previously described.[9]
Definition of a primary outcome: sudden cardiac death
Follow-up of ARIC participants was previously reported[10]. 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.[11] CHD deaths occurring on or before December 31, 2001 were adjudicated by 5 physicians. CHD deaths occurring between January 1, 2002 and December 31, 2012 were adjudicated by a committee of 11 physicians. 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. Witnessed SCD with an available rhythm strip of life-threatening arrhythmia, or primary emergency medical services impression of cardiac arrest constituted definite SCD events. Probable SCD was defined as SCD with greater uncertainty either due to (i) other clinical conditions that can cloud the issue of exact cause of demise; or (ii) limited information to adjudicate an event with greater confidence. Possible SCDs were adjudicated only at the second phase of reviews (2002-2012). Unwitnessed events with limited available data but specified SCD on a death certificate, without documentation of another cause of death, were defined as possible SCD. Inter‐reviewer agreement for the phase 2 effort was 83.2% and agreement across phases was 92.5%. Participants were censored at the time of loss to follow‐up 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 as described in the Supplement. Fatal CHD that did not meet 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 12-hour fast followed by a light snack and at least one hour 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 seconds of each of the leads simultaneously (I, II, III, aVR, aVL, aVF, Vl-V6).
Because it is essential for the study to be able to compare baseline ECG data with subsequent records, a uniform procedure for electrode placement and skin preparation was required. 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 which may cause errors in marking the electrode locations. For best electrode/skin interface, the electrodes were placed on the skin at least 2-3 minutes before taking the ECG. A pen was used to mark six chest electrode positions. The chest was wiped with a sterile alcohol prep. Left leg electrode was placed on the left ankle (inside). Right leg electrode was placed on the right ankle (inside). Left arm electrode was placed on the left wrist (inside). Right arm electrode was placed on the right wrist (inside). 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, participant’s information was entered into MAC PC. It was required that electrodes must be on the skin for at least 3 minutes 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 artefacts.
Recorded 12-lead ECGs were saved in the memory of the ECG machine, and transmitted to Halifax ECG Computing Center via the phone line. Recorded 12-lead ECGs were stored in the MUSE database (GE Marquette, Milwaukee, WI). Later, MUSE database was transferred from the Halifax ECG Computing Center to Epidemiological Cardiology Research Center (EPICARE, Wake Forest University, NC). Then, stored in a MUSE database (GE Marquette, Milwaukee, WI) ECGs were exported and transferred from EPICARE to Tereshchenko laboratory at the Oregon Health & Science University.
Electrocardiogram analyses
ECG data of all five follow-up visits were analyzed. Traditional ECG intervals were reported by the 12 SL algorithm (GE Marquette Electronics, Milwaukee, WI). QT interval was corrected for heart rate according to Bazett's formula.
We analyzed raw, digital, 10-second, 12-lead ECGs (sampling rate of 500Hz 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 normal sinus median beat was included in the analysis.
GEH was measured as previousl described,[3, 12] in a time-coherent median meat with identified origin of the heart vector.[13] We have provided the software code at Physionet (https://physionet.org/physiotools/geh/). In addition to previously reported “mean” GEH measures,[3] in this study we measured the spatial peak vectors (Figure 1).[12-14]
First, we transformed 12-lead ECG in orthogonal XYZ ECG, using Kors transformation.[15] Next, we constructed time-coherent median beat, and detected origin of the heart vector, using our novel approach, as recently described.[13] Then, we performed calculations of GEH metrics, using the following equations.
Due to technical limitations, the equations have been placed in the supplementary files section.
Statistical analysis
Time-dependent area under the (receiver operating characteristic, ROC) curve (AUC) analysis was performed to assess 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[16, 17]. 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 equation 0.25* ), where n is the number of observations. All available five visits’ ECG data were included in time-dependent AUC analysis.[3] To satisfy requirement for ROC analysis, and to perform internal validation of study findings, we divided the dataset in 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 a time from ECG recording to the time of outcome (or censoring). If 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 95% confidence interval (CI) of ROC AUC in each partition, separately. Then, ROC AUC point estimates, and lower and upper boundaries of 95%CI were averaged across 5 partitions and presented as a final summary result.
We summarized clinical characteristics of study participants with a 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[3] (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.[18, 19] 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.