Periodic repolarization dynamics as predictor of risk for sudden cardiac death in chronic heart failure patients

The two most common modes of death among chronic heart failure (CHF) patients are sudden cardiac death (SCD) and pump failure death (PFD). Periodic repolarization dynamics (PRD) quantifies low-frequency oscillations in the T wave vector of the electrocardiogram (ECG) and has been postulated to reflect sympathetic modulation of ventricular repolarization. This study aims to evaluate the prognostic value of PRD to predict SCD and PFD in a population of CHF patients. 20-min high-resolution (1000 Hz) ECG recordings from 569 CHF patients were analyzed. Patients were divided into two groups, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {PRD}^+$$\end{document}PRD+ and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {PRD}^-$$\end{document}PRD-, corresponding to PRD values above and below the optimum cutoff point of PRD in the study population. Univariate Cox regression analysis showed that SCD risk in the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {PRD}^+$$\end{document}PRD+ group was double the risk in the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {PRD}^-$$\end{document}PRD- group [hazard ratio (95% CI) 2.001 (1.127–3.554), \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {p}<0.05$$\end{document}p<0.05]. The combination of PRD with other Holter-based ECG indices, such as turbulence slope (TS) and index of average alternans (IAA), improved SCD prediction by identifying groups of patients at high SCD risk. PFD could be predicted by PRD only when combined with TS [hazard ratio 2.758 (1.572–4.838), \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {p}<0.001$$\end{document}p<0.001]. In conclusion, the combination of PRD with IAA and TS can be used to stratify the risk for SCD and PFD, respectively, in CHF patients.

Association of PRD with cardiac events. PRD was dichotomized according to the optimum cutoff point obtained from ROC analysis for each of the investigated endpoints. Optimum cutoff values were 1.33 • for SCD, 1.31 • for PFD (approximately coincident with the median PRD over the study population) and 1.32 • for cardiac death (CD), comprising SCD and PFD. In each case, PRD + and PRD − subgroups were defined to contain patients with PRD values above and below the cutoff point, respectively. Additional information on the derivation of the cutoff values is provided in "Methods" section.
The percentages of SCD, PFD and CD victims in PRD + and PRD − groups are presented in Fig. 1. There were 106 CD victims in the study population (18.6%), of which 53 died of SCD and 53 of PFD. The SCD mortality rate was significantly higher (p = 0.018) in the PRD + group (33 victims, 11.7%) than in the PRD − group (20 victims, 7%). Regarding PFD, no significant differences in mortality rates between PRD + and PRD − groups were found (28 victims, 10% in PRD + vs 25 victims, 8.7% in PRD − ). For CD, significant differences were found, with mortality being remarkably higher in the PRD + group as compared to the PRD − group (22.3% vs 15%).

SCD and PFD risk prediction based on PRD and other individual variables. Kaplan-Meier sur-
vival analysis showed that PRD + patients had significantly lower SCD survival probability than PRD − patients ( p = 0.024 ), as illustrated in Fig. 2. When accounting for PFD as a competing risk, Fine and Gray analysis 29 showed that the cumulative SCD incidence curves for PRD + and PRD − patients were statistically significantly different ( p = 0.048 ). When PFD was considered as endpoint, no significant differences in the survival rates were found between PRD + and PRD − groups.
Univariate Cox analysis results for SCD risk prediction by the PRD index as well as by demographical, clinical and other ECG variables are shown in Table 2. The variables significantly associated with SCD were New York Heart Association (NYHA) class III, LVEF ≤ 35% , NSVT and VPB > 240 , N-terminal prohormone of brain natriuretic peptide (NT-proBNP) > 1000 pg/mL , IAA + , TS + and PRD + . None of the HRV variables (LFn, HFn www.nature.com/scientificreports/ or LF/HF) was associated with SCD risk. Regarding PFD, PRD + was not able to predict this endpoint, whereas age, ischemic etiology, prior myocardial infarction, NYHA class, NT-proBNP > 1000 pg/mL and TS + were significant predictors (Table 3). Multivariate Cox proportional hazards regression model for SCD prediction including clinical variables such as the combination of NSVT and VPB > 240 , LVEF ≤ 35% , NYHA class and NT-proBNP > 1000 pg/mL together with PRD + led to the results presented in Table 4. PRD + and some of the clinical variables were independent predictors of SCD, being associated with similar hazard ratio (HR) values. When considering PFD as a competing risk event, the regression model for SCD as endpoint included clinical variables and PRD, the latter with a HR of 1.65, even if being only marginally predictive ( p = 0.095).
When tested in the subpopulation of patients with CHF of non-ischemic etiology, PRD was able to predict SCD in the multivariate model ( HR = 2.497 , p = 0.05 ). In the subpopulation of patients with CHF of ischemic etiology, PRD was not predictive of SCD in the univariate model and therefore it was not included in a multivariate model.
We finally tested the capacity of PRD + for CD prediction. The HRs for univariate and multivariate Cox analyses were 1.64 ( p = 0.014 ) and 1.63 ( p = 0.016 ), respectively. Other demographical, clinical and ECG variables were also able to predict CD in a univariate Cox model (Table 5).

SCD and PFD risk prediction based on the combination of PRD with other variables.
To improve the predictive power for SCD risk, the combination of PRD with other Holter-based variables, such as IAA and TS, was assessed. As shown in Fig. 3, patients with PRD + and TS + had significantly higher SCD mortality ( p < 0.05 ) and PFD mortality ( p < 0.001 ) than the rest of patients. CD mortality was significantly higher too. Results for the combination of PRD and IAA are presented in Fig. 4. Patients with PRD + and IAA + had increased   www.nature.com/scientificreports/ SCD mortality ( p = 0.007 ) and CD mortality ( p = 0.023 ) than the rest of the population. No significant differences were found with respect to PFD mortality ( p = 0.816). Figures 5 and 6 show Kaplan-Meier probability curves for SCD survival using the combined variables PRD + and TS + and PRD + and IAA + , respectively. For the two variables, significantly lower survival was found for patients with positive PRD and TS or IAA. When PFD was considered as endpoint, differences in survival probabilities were only statistically significant when using PRD + and TS + . Considering PFD as a competing risk event for SCD, Fine and Gray analysis showed significantly higher cumulative SCD incidence for PRD + and TS + ( PRD + and IAA + , respectively) patients than the rest of the population ( p < 0.001 in all cases). When accounting for SCD as a competing risk event for PFD, the cumulative PFD incidence for patients with PRD + and TS + was significantly higher than for the rest of patients in the population ( p < 0.001).
In univariate Cox analysis (Table 6) with SCD as endpoint, PRD + and TS + as well as PRD + and IAA + were risk predictors, with associated HRs of 3.1 and 2.8, respectively. For PFD as endpoint, PRD + and TS + patients presented more than two and a half times higher risk than the rest of the population. For CD as endpoint, PRD + and TS + was associated with a HR of 3.1, while PRD + and IAA + with a HR of 1.9.
Results from multivariate Cox proportional hazard regression analysis including NYHA class, LVEF ≤ 35% , NSVT and VPB > 240 and NT-proBNP > 1000 pg/mL as well as one of the two combined ECG variables at a    Table 7 for SCD as endpoint. Both PRD + and IAA + and PRD + and TS + independently predicted SCD risk with HRs of 3.0 and 2.1, respectively. Multivariate regression analysis for SCD with PFD as competing risk indicated that the combined variable PRD + and IAA + ( PRD + and TS + , respectively) predicted SCD independently of other clinical variables, with an associated HR of 2.9 (2.4, respectively, being p < 0.01 in all cases).
In the subpopulation of patients with CHF of non-ischemic etiology, PRD + and TS + predicted SCD with a HR 3.2 (p = 0.012) and PFD with a HR of 3.43 ( p = 0.011 ) independently of other clinical variables. PRD + and IAA + predicted SCD risk in both ischemic and non-ischemic etiology subpopulations, with associated HR of 3.3 and 2.9, respectively.
To complete the study, we tested the capacity of the combined variables for CD risk prediction. Table 8 shows the results for the two tested multivariable proportional hazard models. Both PRD + and IAA + and PRD + and TS + predicted CD risk independently of other demographic and clinical variables with HRs of 2.4 and 2.0, respectively.

Discussion
This work investigates the performance of PRD, characterizing the oscillatory behavior of the T wave, for cardiac risk stratification in a population of CHF patients. Although PRD has been previously used in various contexts 19,27,[30][31][32][33][34] , in this study, we computed PRD in a large CHF cohort to predict the risk of the two most common modes of death, namely SCD and PFD, both considering PRD on its own and in combination with other ECG-based indices.
Our results show that CHF patients presenting high PRD values have nearly twice higher risk of suffering SCD than CHF patients with low PRD. These findings are in agreement with previous clinical studies showing a relationship between increased PRD and enhanced mortality, particularly due to ventricular arrhythmias, in different patient populations 19,35 . Specifically in post-myocardial infarction patients with impaired LVEF, elevated PRD has been related to enhanced risk for SCD 20 . Also, PRD has been suggested to be potentially useful in guiding clinical decisions on prophylactic implantation of ICDs in patients with ischemic and non-ischemic cardiomyopathy 36 .
PRD measures the magnitude of LF oscillations in the T wave vector of the ECG and has been postulated to reflect sympathetic modulation of ventricular repolarization 19 . The same range of LF oscillations in the T wave vector accounted for by PRD has been reported for the action potential duration (APD) in in vivo studies on HF patients 37 . Through in silico simulations, synergistic β-adrenergic and mechanical effects induced by sympathetic activation have been shown to contribute to the generation of these LF oscillations of APD 38,39 . Those results have been supported by subsequent experimental studies showing that increased sympathetic activity potentiates such oscillations and β-adrenergic blockade attenuates them 37,40 . In the presence of calcium overload and/ or reduced repolarization reserve, both being commonly associated with HF, the amplitude of LF oscillations of APD have been theoretically shown to be magnified, facilitating the occurrence of arrhythmogenic events 39,41 . A recent in vivo study in a canine model of ventricular remodeling caused by chronic atrioventricular block has provided evidence that dogs inducible for ventricular arrhythmias present higher LF oscillations of repolarization than non-inducible dogs 42 . Based on all these observations, the results of our study possibly suggest that CHF patients with high PRD present augmented repolarization variability that can lead to destabilization of repolarization and promote arrhythmogenesis.
Since CHF is a complex clinical syndrome, some studies have reported the benefit of using combined ECG risk markers and/or a risk score integrating information from various clinical and ECG variables to improve clinical decision making [43][44][45][46][47] . The use of markers providing information about different pathophysiological processes Table 8. Multivariable CD risk prediction including the following variables: age, ischemic etiology, prior myocardial infarction, NYHA class, LVEF ≤ 35% , combined NSVT and VPB > 240 , NT-proBNP > 1000 pg/mL and a combined ECG variable that can be either PRD + and TS + or PRD + and IAA + . * p < 0.05 , * * p < 0.01 , * * * p < 0.001. www.nature.com/scientificreports/ associated with CHF has been shown to be very useful 48,49 . Here, a combination of PRD and other ECG-based indices are tested for SCD and PFD risk prediction. For PFD prediction, the combination of PRD with TS, measuring heart rate turbulence, stratifies the study population into two groups according to the risk of PFD. The risk of dying from PFD is two and a half times higher in patients with high PRD ( PRD + ) and low TS ( TS + ) than in the rest of patients. Nevertheless, the capacity for PFD prediction of TS on its own is already very good, as reported in previous studies 50 , and PRD only slightly improves it. TS represents a vagally-mediated response of heart rate to ventricular premature beats with involvement of baroreflex sensitivity. TS has been considered an excellent marker of neurohormonal activation in CHF patients 50,51 . Since autonomic dysfunction and neurohormonal activation have a relevant role in the progression of CHF, our results on the capacity for PFD risk stratification of combined PRD and TS support assessment of such CHF landmarks noninvasively from short-term ECGs (for PRD evaluation) and Holter recordings (for TS evaluation). In terms of SCD prediction, the stratifying ability of PRD is enhanced when combined with the marker IAA evaluating T-wave alternans amplitude over 24 h.
While T-wave alternans, measured by IAA in this work, is related to abnormal cardiac function 17 , PRD measures changes in sympathetic modulation of ventricular repolarization 19 . The synergistic information provided by the two variables serves to improve SCD prediction. In all the above described findings, the risk stratification capacity of PRD is maintained when accounting for competing risk events in survival analyses. NT-proBNP concentration, dichotomized using the threshold (1000 pg/mL) proposed in the original MUSIC study 28 , was not a significant predictor of SCD in multivariable regression models comprising PRD or combined ECG markers. When the endpoint was CD, NT-proBNP > 1000 pg/mL could significantly predict CD in a multivariable model. Other studies have proposed alternative thresholds for NT-proBNP, like NT-proBNP > 5180 pg/mL 52 , and found the dichotomized variable to be a predictor of both SCD and PFD. However, as our CHF population includes only patients in NYHA classes II and III, only 7% of the patients had NT-proBNP values above 5180 pg/mL.
On the basis of previous studies showing that PRD is unrelated to HRV and respiratory activity 19 , we hypothesized that PRD and HRV markers could add complementary information. However, our results show that HRV markers do not have a relationship either with SCD or PFD and, thus, we have not combined them with PRD for risk stratification purposes. On the other hand, of the clinical variables assessed in our study, NYHA class and LVEF ≤ 35% are the two with the highest HR for both PFD and SCD. In any case, the HR of these two clinical variables are lower than those obtained for the combined ECG variables PRD and TS in the prediction of PFD and PRD and IAA in the prediction of SCD. LVEF is widely used in the clinics to identify high-risk patients but its accuracy is low 53 . Regarding NYHA class, as it reflects a subjective assessment and can change frequently over short periods of time 28 , its interpretation is more critical than other variables. The analysis conducted in this work suggests that the adjunct use of clinical variables and ECG-based markers can improve the prognostic value of all of them and render a risk score with remarkably superior performance to predict the two most common modes of death in CHF patients.
An important aspect regarding the use of PRD is the fact that it can be measured from 5-min ECGs by using the phase-rectified signal averaging (PRSA)-based method employed in this study, first described by Rizas et al. 35 and subsequently updated by Palacios et al. 27 Here, we use 20-min ECGs, from which the 5-min segment associated with minimum PRD is selected. This is an advantage as compared to a number of previously analyzed markers, such as IAA 17 , QTVi 12 , TS 54 , �α QT55 , �α Tpe56 or TMR 57 , which require longer ECG recordings or specific protocols, like stress tests or steady heart rate, for their evaluation. While the combination of PRD with other Holter-based markers improves risk stratification, it requires long-duration recordings. On this basis, we propose that PRD could be used in an initial step to select patients with high CD risk. Over those selected patients, the use of longer duration signals and associated variables could be useful to specifically predict each of the death modes. Several works have proposed other markers related to repolarization variability measured from short-term ECG recordings, such as the variance normalized by the mean of QT end, QT peak or T-peak-to-T-end (Te) intervals, and have evaluated them in CHF populations. In recent studies, the mean and/or standard deviation of Te have been shown to predict 30-day mortality 58,59 and mortality in hospital 60 among decompensated CHF patients.
To sum up, our study documents the prognostic value of PRD in combination with other Holter-derived markers to predict PFD and SCD in a large cohort of CHF patients. Future studies on CHF cohorts including a larger number of SCD and PFD victims would allow conducting more robust statistical analysis to confirm the findings here reported. Also, studies on larger populations would facilitate assessing the capacity of PRD for risk stratification in specific CHF patient subpopulations, such as those with reduced or preserved LVEF or at different NYHA classes, among others.

Methods
Study population. The MUSIC study is a prospective, multicenter, longitudinal study designed to assess risk predictors of cardiac mortality and SCD in ambulatory patients with CHF. Patients were consecutively enrolled from the specialized CHF clinics of eight University Hospitals between April 2003 and December 2004. All had symptomatic CHF (NYHA class II-III) and were treated according to institutional guidelines. The study included patients with either depressed or preserved LVEF. Patients with preserved LVEF were included if they had CHF symptoms and a prior hospitalization for HF or some objective signs of HF confirmed by chest X-ray (findings of pulmonary congestion) and/or echocardiography (abnormal LV filling pattern and LV hypertrophy). Patients were excluded if they had recent acute coronary syndrome or severe valvular disease amenable for surgical repair. Patients with other concomitant diseases expected to reduce life-expectancy were also excluded. All patients gave written informed consent and the study protocol was approved by all the institutional investigation and ethics committees from the following participant hospitals: Valme Hospital, Santiago de Compostela Hospital, Son Dureta Hospital, Arrixaca Hospital, Gregorio Marañon Hospital, Joan XXIII Hospital, Insular Las Study protocol. For each patient, two ECG recordings were available: a 24-h ambulatory ECG sampled at 200 Hz in 3 orthogonal X, Y and Z leads using SpiderView recorders (ELA Medical, SorinGroup, Paris, France) and a prior 20-min ECG sampled at 1000 Hz while patients were resting in supine position. In this study, PRD was measured from the 20-min ECG recording. Patients were followed every 6 months on an outpatient basis for an average of 48 months. Death was defined as SCD if it was: (1) a witnessed death occurring within 60 min from the onset of new symptoms, unless a cause other than cardiac was obvious; (2) an unwitnessed death ( < 24 h ) in the absence of preexisting progressive circulatory failure or other causes of death; or (3) a death during attempted resuscitation. Deaths occurring in hospital as a result of refractory progressive end-stage CHF, or CHF patients undergoing heart transplantation, were defined as PFD. Endpoints were reviewed and classified by the MUSIC Study Endpoint Committee. A description of clinical data for the overall population and more detailed information about the study protocol can be found in Vazquez et al. 28 .
ECG pre-processing. High-resolution ECG signals were preprocessed by a 50 Hz notch filter to remove powerline interference. QRS complexes were detected using Aristotle software 61 . The heart timing method described in Mateo and Laguna 62 was applied to detect irregular sinus beats. As the main analysis was focused on the T-wave, a 40-Hz low-pass filter was subsequently applied onto the ECG signals to remove noise without distorting the T-wave shape. Finally, baseline wander was estimated by cubic splines interpolation and cancelled. For each beat i, a window including the T-wave was defined from the QRS position, denoted by QRS i , and the associated RR interval preceding it, denoted by RR i . The T-wave window onset, denoted by T on i , was set at 90 ms after the QRS i mark: T on i = QRS i + 90 ms . The T-wave window end, denoted by T end i , was defined as T end i = QRS i + min(360 ms, 2 3 RR i ) for RR i below 720 ms. For RR i equal or higher than 720 ms, T end i = QRS i + 360 ms.
The noise level was estimated in each lead and beat by measuring the root mean square value of the highfrequency components (above 15 Hz) in a window around the T-wave, similarly to previous works 56 . A global measure of noise for each beat i was obtained by summing the noise levels of all leads and denoted by V RMS noise (i) . If a beat was too noisy ( V RMS noise >140 µV), the associated T-wave marks ( T on and T end ) were discarded for further analysis, but the associated QRS mark was retained. PRD computation. PRD was measured in the preprocessed signal by using the method described in Palacios et al. 27 , a modified version of the original method described in Rizas et al. 19 . This method included the following steps: 1. T waves were selected using the T wave windows defined in section ECG pre-processing. 2. A constant value was subtracted from each T wave in each of the analyzed leads so that the amplitude at T end was set to 0 mV. 3. The average electrical vector was calculated for each T-wave window. The angle dT • between two consecutive T-wave windows was calculated by the dot product of the corresponding average vectors. 4. A 10th-order median filter was used to attenuate outliers and artifacts in the dT • time series.
Steps (i) to (iii) are illustrated in Fig. 7.
Over the obtained dT • time series, a method based on PRSA 63 was applied to evaluate the oscillations measured by the PRD index 35  The value of M = 9 was established because it allows detecting frequencies in the range of interest (from 0.025 to 0.1 Hz), as fully described in Bauer et al. 63 . 6. Windows of 2L values were defined around each anchor point. Anchor points in the last L samples of the dT • series were discarded, as windows of length 2L could not be defined around them. In this study, L = 20 was chosen, as it allowed detection of frequencies in the range of interest. 7. The PRSA series was obtained by averaging the dT • values over all 2L-sample windows contained in each 5-min segment.
For each 5-min segment, a PRD value was defined as the difference between the maximum and minimum values of the PRSA series. For the 20-min recording of each subject, a unique PRD value was calculated as the minimum PRD over the analyzed 5-min segments with 4-min overlap. Illustrative examples of dT • and PRSA time series from two subjects, a SCD victim and a survivor, are shown in Fig. 8.
By setting a cutoff point of 1.33 • , an optimum thresholds was identified as that maximizing the geometric mean of sensitivity and specificity for SCD as endpoint, two groups were defined: PRD + , containing those patients     www.nature.com/scientificreports/ Heart rate variability analysis. For the 5-min segment with minimum PRD, HRV indices were computed as in Bailón et al. 64 . Power spectral density (PSD) was estimated based on the periodogram. LF and highfrequency (HF) powers were calculated as the areas under the PSD within the 0.04-0.15 Hz and 0.15-0.5 Hz frequency bands, respectively. Normalized powers, denoted by LFn and HFn, were obtained by dividing LF and HF powers by the sum of the two. The ratio between LF and HF powers, denoted by LF/HF, was additionally computed 65 .

Clinical variables and other Holter-based ECG indices. Clinical variables such as LVEF ≤ 35% ,
NT-proBNP > 1000 pg/mL , non-sustained ventricular tachycardia (NSVT) and a number of ventricular premature beats (VPB) > 240 , whose capacity for SCD and PFD risk prediction has been previously shown 28 , were included in our analysis. Additionally, other Holter-based indices calculated in previous studies of our group were considered based on their risk prediction power. On the one hand, the index of average alternans, IAA, quantifying the average amplitude of T-wave alternans over a 24-h period, was computed by automatic ECG analysis 17 . IAA was shown to risk stratify for SCD in the study population here analyzed when dichotomized at 3.7 µ V to define IAA + and IAA − groups 17 . On the other hand, turbulence slope, TS, describing the initial phase of sinus rhythm deceleration after a VPB, was determined as the maximum positive slope of a regression line assessed over any of 5 consecutive RR intervals within the first 20 sinus RR intervals after a VPB during the 24-h ECG recording 51 . TS was shown to stratify risk for both SCD and PFD when dichotomized at 2.5 ms/RR, with TS + , denoting the group of patients with TS below the threshold, being associated with higher risk than the group TS − containing the remaining patients 51 .
Statistical analysis. Continuous variables are presented as median [interquartile range (IQR)]. The Mann-Whitney U test (or Wilcoxon rank-sum test) was used for univariate comparisons of continuous variables between patient groups. Survival probability was estimated by Kaplan-Meier analysis, with the log-rank test used to assess group differences. Univariate and multivariate Cox 66 proportional hazards regression models and Fine and Gray analysis for competing risk 29 were used for prediction of endpoints. Hazard ratios (HR) and 95% confident intervals (CI), expressed as HR [95% CI] were quantified. P values < 0.05 were considered for statistical significance. Variables being significantly different between groups in the univariate analysis were input to the multivariate Cox regression model. Data were analyzed using MATLAB R2017a (9.2), SPSS (version 24.0) and R software (version 4.1).

Conclusions
This study tests PRD, a non-invasive marker of repolarization instability associated with low-frequency oscillations in sympathetic activity, as a predictor of SCD and PFD risk in CHF patients. The combination of PRD with an index of T wave alternans further enhances the capacity of PRD for SCD risk stratification. Additionally, the combination of PRD with heart rate turbulence slope allows predicting PFD risk. Altogether, the value of ECG markers, either derived from short-term ECGs or ambulatory Holter recordings, is highlighted as a means to improve prognosis in CHF patients beyond commonly used clinical variables.

Data availability
The current study analyzed datasets which are not publicly available due to restrictions in the ethical permission but the data can be accessed through the corresponding author upon reasonable request and with permission of the MUSIC Study Committee.