Time-Warping Analysis of the T-Wave Peak-to-End Interval to Quantify Ventricular Repolarization Dispersion During Ischemia

Variations in the dispersion of ventricular repolarization can be quantified by T-wave time-warping based index, dw. However, the early phase of the T-wave can be affected by ST-segment changes during ischemia. We hypothesized that restricting dw to the T-wave peak-to-end ( Tpe) would circumvent this limitation while still quantifying variations in repolarization dispersion. A total of 101 ECG recordings from patients undergoing coronary occlusion, together with their control recordings, were analyzed. A series of dw values was calculated by quantifying the Tpe morphological variations between the T-waves at different occlusion stages and a baseline T-wave. We introduced a normalized version of dw, Rd, reflecting variations of dw during occlusion relative to control recordings ( Rd = 1 corresponds to the same level of variation). The dw series followed a gradually increasing trend with occlusion time, reaching median [range] Rd values of 9.44 [1.01, 80.74] at the occlusion end. Rd at occlusion end was significantly higher than threshold values of 1, 2, 5, and 10 in 94.1%, 85.11%, 64.4% and 48.5% of patients, respectively. The spatial lead-wise analysis of dw showed distinct distributions depending on the occluded artery, suggesting a relation with the ischemia location. The relative variation R with ischemia of index dw (9.4) is greater than that of the T-wave amplitude (7.7), Tpe interval (2.7) and T-wave width (3.0). In conclusion, dw tracks ischemic-induced variations in repolarization dispersion in a more robust manner than classical indexes, avoiding the impact of ST segment elevation/depression or early T-wave distortions, thus warranting further clinical studies.

Differences in ventricular wall composition may explain inter-patient variability in the ECG response to variations in serum potassium and calcium

Objective: Chronic kidney disease patients have a decreased ability to maintain normal electrolyte concentrations in their blood, which increases the risk for ventricular arrhythmias and sudden cardiac death. Non-invasive monitoring of serum potassium and calcium concentration, [K+] and [Ca2+], can help to prevent arrhythmias in these patients. Electrocardiogram (ECG) markers that significantly correlate with [K+] and [Ca2+] have been proposed, but these relations are highly variable between patients. We hypothesized that inter-individual differences in cell type distribution across the ventricular wall can help to explain this variability. Methods: A population of human heart-torso models were built with different proportions of endocardial, midmyocardial and epicardial cells. Propagation of ventricular electrical activity was described by a reaction-diffusion model, with modified Ten Tusscher-Panfilov dynamics. [K+] and [Ca2+] were varied individually and in combination. Twelve-lead ECGs were simulated and the width, amplitude and morphological variability of T waves and QRS complexes were quantified. Results were compared to measurements from 29 end-stage renal disease (ESRD) patients undergoing hemodialysis (HD). Results: Both simulations and patients data showed that most of the analyzed T wave and QRS complex markers correlated strongly with [K+] (absolute median Pearson correlation coefficients, r, ranging from 0.68 to 0.98) and [Ca2+] (ranging from 0.70 to 0.98). The same sign and similar magnitude of median r was observed in the simulations and the patients. Different cell type distributions in the ventricular wall led to variability in ECG markers that was accentuated at high [K+] and low [Ca2+], in agreement with the larger variability between patients measured at the onset of HD. The simulated ECG variability explained part of the measured inter-patient variability. Conclusion: Changes in ECG markers were similarly related to [K+] and [Ca2+] variations in our models and in the ESRD patients. The high inter-patient ECG variability may be explained by variations in cell type distribution across the ventricular wall, with high sensitivity to variations in the proportion of epicardial cells. Significance: Differences in ventricular wall composition help to explain inter-patient variability in ECG response to [K+] and [Ca2+]. This finding can be used to improve serum electrolyte monitoring in ESRD patients.

ECG T-Wave Morphologic Variations Predict Ventricular Arrhythmic Risk in Low- and Moderate-Risk Populations

Background Early identification of individuals at risk of sudden cardiac death (SCD) remains a major challenge. The ECG is a simple, common test, with potential for large-scale application. We developed and tested the predictive value of a novel index quantifying T-wave morphologic variations with respect to a normal reference (TMV), which only requires one beat and a single-lead ECG. Methods and Results We obtained reference T-wave morphologies from 23 962 participants in the UK Biobank study. With Cox models, we determined the association between TMV and life-threatening ventricular arrhythmia in an independent data set from UK Biobank study without a history of cardiovascular events (N=51 794; median follow-up of 122 months) and SCD in patients with coronary artery disease from ARTEMIS (N=1872; median follow-up of 60 months). In UK Biobank study, 220 (0.4%) individuals developed life-threatening ventricular arrhythmias. TMV was significantly associated with life-threatening ventricular arrhythmias (hazard ratio [HR] of 1.13 per SD increase [95% CI, 1.03-1.24]; P=0.009). In ARTEMIS, 34 (1.8%) individuals reached the primary end point. Patients with TMV ≥5 had an HR for SCD of 2.86 (95% CI, 1.40-5.84; P=0.004) with respect to those with TMV <5, independently from QRS duration, corrected QT interval, and left ventricular ejection fraction. TMV was not significantly associated with death from a cause other than SCD. Conclusions TMV identifies individuals at life-threatening ventricular arrhythmia and SCD risk using a single-beat single-lead ECG, enabling inexpensive, quick, and safe risk assessment in large populations.

Monitoring of Serum Potassium and Calcium Levels in End-Stage Renal Disease Patients by ECG Depolarization Morphology Analysis

Objective: Non-invasive estimation of serum potassium, [K+], and calcium, [Ca2+], can help to prevent life-threatening ventricular arrhythmias in patients with advanced renal disease, but current methods for estimation of electrolyte levels have limitations. We aimed to develop new markers based on the morphology of the QRS complex of the electrocardiogram (ECG). Methods: ECG recordings from 29 patients undergoing hemodialysis (HD) were processed. Mean warped QRS complexes were computed in two-minute windows at the start of an HD session, at the end of each HD hour and 48 h after it. We quantified QRS width, amplitude and the proposed QRS morphology-based markers that were computed by warping techniques. Reference [K+] and [Ca2+] were determined from blood samples acquired at the time points where the markers were estimated. Linear regression models were used to estimate electrolyte levels from the QRS markers individually and in combination with T wave morphology markers. Leave-one-out cross-validation was used to assess the performance of the estimators. Results: All markers, except for QRS width, strongly correlated with [K+] (median Pearson correlation coefficients, r, ranging from 0.81 to 0.87) and with [Ca2+] (r ranging from 0.61 to 0.76). QRS morphology markers showed very low sensitivity to heart rate (HR). Actual and estimated serum electrolyte levels differed, on average, by less than 0.035 mM (relative error of 0.018) for [K+] and 0.010 mM (relative error of 0.004) for [Ca2+] when patient-specific multivariable estimators combining QRS and T wave markers were used. Conclusion: QRS morphological markers allow non-invasive estimation of [K+] and [Ca2+] with low sensitivity to HR. The estimation performance is improved when multivariable models, including T wave markers, are considered. Significance: Markers based on the QRS complex of the ECG could contribute to non-invasive monitoring of serum electrolyte levels and arrhythmia risk prediction in patients with renal disease.

Estimation of potassium levels in hemodialysis patients by T wave nonlinear dynamics and morphology markers

OBJECTIVE

Noninvasive screening of hypo- and hyperkalemia can prevent fatal arrhythmia in end-stage renal disease (ESRD) patients, but current methods for monitoring of serum potassium (K+) have important limitations. We investigated changes in nonlinear dynamics and morphology of the T wave in the electrocardiogram (ECG) of ESRD patients during hemodialysis (HD), assessing their relationship with K+ and designing a K+ estimator.

METHODS

ECG recordings from twenty-nine ESRD patients undergoing HD were processed. T waves in 2-min windows were extracted at each hour during an HD session as well as at 48 h after HD start. T wave nonlinear dynamics were characterized by two indices related to the maximum Lyapunov exponent (λt, λwt) and a divergence-related index (η). Morphological variability in the T wave was evaluated by three time warping-based indices (dw, reflecting morphological variability in the time domain, and da and daNL, in the amplitude domain). K+was measured from blood samples extracted during and after HD. Stage-specific and patient-specific K+ estimators were built based on the quantified indices and leave-one-out cross-validation was performed separately for each of the estimators.

RESULTS

The analyzed indices showed high inter-individual variability in their relationship with K+. Nevertheless, all of them had higher values at the HD start and 48 h after it, corresponding to the highest K+. The indices η and dw were the most strongly correlated with K+ (median Pearson correlation coefficient of 0.78 and 0.83, respectively) and were used in univariable and multivariable linear K+ estimators. Agreement between actual and estimated K+ was confirmed, with averaged errors over patients and time points being 0.000 ± 0.875 mM and 0.046 ± 0.690 mM for stage-specific and patient-specific multivariable K+ estimators, respectively.

CONCLUSION

ECG descriptors of T wave nonlinear dynamics and morphological variability allow noninvasive monitoring of K+ in ESRD patients.

SIGNIFICANCE

ECG markers have the potential to be used for hypo- and hyperkalemia screening in ESRD patients.

Weighted Time Warping Improves T-wave Morphology Markers Clinical Significance

Background: T-wave (TW) morphology indices based on time-warping (d_w) have shown significant cardiovascular risk stratification value. However, errors in the location of TW boundaries may impact their prognostic power. Our aim was to test the hypothesis that a weighted time-warping function (WF) would reduce the sensitivity of d_w to these errors and improve their clinical significance. Methods: The WFs were proportional to (i) the reference TW (T), and (ii) the absolute value of its derivative (D). The index d_w was recalculated using these WFs, and its performance was compared to the unweighted control case (C) in four different scenarios: 1) robustness against simulated TW boundaries location errors; 2) ability to retain physiological information in an electrophysiological cardiac model; 3) ability to monitor blood potassium concentration changes ([K⁺]) in 29 hemodialysis (HD) patients; 4) and the sudden cardiac death (SCD) risk stratification value of the TW morphology restitution (TMR) index, derived from d_w, in 651 chronic heart failure (CHF) patients. Results and Discussion: The WFs led to a reduced sensitivity (R) of d_w to TW boundary location errors as compared to C (median R=0.19 and 0.22 and 0.35 for T, D and C, respectively). They also preserved the physiological relationship between d_w and repolarization dispersion changes at ventricular level. No improvements in [K⁺] tracking were observed for the HD patients (Pearsons median correlation [r] between [K⁺] and d_w was 0.86r0.90 for T, D and C). In CHF patients, the SCD risk stratification value of TMR was improved by applying T (hazard ratio, HAR, of 2.80), followed by D (HAR=2.32) and C (HAR=2.23). Conclusions and Significance: The proposed WFs, with T showing the best performance, increased the robustness of time-warping based markers against TW location errors preserving their physiological information content and boosting their SCD risk stratification value. Results from this work support the use of T when deriving d_w for future clinical applications.

Inference of ventricular activation properties from non-invasive electrocardiography

The realisation of precision cardiology requires novel techniques for the non-invasive characterisation of individual patients’ cardiac function to inform therapeutic and diagnostic decision-making. Both electrocardiography and imaging are used for the clinical diagnosis of cardiac disease. The integration of multi-modal datasets through advanced computational methods could enable the development of the cardiac ‘digital twin’, a comprehensive virtual tool that mechanistically reveals a patient's heart condition from clinical data and simulates treatment outcomes. The adoption of cardiac digital twins requires the non-invasive efficient personalisation of the electrophysiological properties in cardiac models. This study develops new computational techniques to estimate key ventricular activation properties for individual subjects by exploiting the synergy between non-invasive electrocardiography, cardiac magnetic resonance (CMR) imaging and modelling and simulation. More precisely, we present an efficient sequential Monte Carlo approximate Bayesian computation-based inference method, integrated with Eikonal simulations and torso-biventricular models constructed based on clinical CMR imaging. The method also includes a novel strategy to treat combined continuous (conduction speeds) and discrete (earliest activation sites) parameter spaces and an efficient dynamic time warping-based ECG comparison algorithm. We demonstrate results from our inference method on a cohort of twenty virtual subjects with cardiac ventricular myocardial-mass volumes ranging from 74 cm³ to 171 cm³ and considering low versus high resolution for the endocardial discretisation (which determines possible locations of the earliest activation sites). Results show that our method can successfully infer the ventricular activation properties in sinus rhythm from non-invasive epicardial activation time maps and ECG recordings, achieving higher accuracy for the endocardial speed and sheet (transmural) speed than for the fibre or sheet-normal directed speeds.

•

Estimation of the ventricular speeds and earliest activation sites from ECG and CMR.

•

Evaluation with twenty virtual subjects shows the effect of anatomical variability.

•

Bayesian-inspired simultaneous estimation of continuous and discrete parameters.

•

Efficient dynamic time warping-based comparison of electrocardiograms (ECG).

•

Changing fibre and sheet-normal speed does not affect healthy activation sequence.

Analysing electrocardiographic traits and predicting cardiac risk in UK biobank

The electrocardiogram (ECG) is a commonly used clinical tool that reflects cardiac excitability and disease. Many parameters are can be measured and with the improvement of methodology can now be quantified in an automated fashion, with accuracy and at scale. Furthermore, these measurements can be heritable and thus genome wide association studies inform the underpinning biological mechanisms. In this review we describe how we have used the resources in UK Biobank to undertake such work. In particular, we focus on a substudy uniquely describing the response to exercise performed at scale with accompanying genetic information.

Nonlinear T-Wave Time Warping-Based Sensing Model for Non-Invasive Personalised Blood Potassium Monitoring in Hemodialysis Patients: A Pilot Study

Background: End-stage renal disease patients undergoing hemodialysis (ESRD-HD) therapy are highly susceptible to malignant ventricular arrhythmias caused by undetected potassium concentration ([K+]) variations (Δ[K+]) out of normal ranges. Therefore, a reliable method for continuous, noninvasive monitoring of [K+] is crucial. The morphology of the T-wave in the electrocardiogram (ECG) reflects Δ[K+] and two time-warping-based T-wave morphological parameters, dw and its heart-rate corrected version dw,c, have been shown to reliably track Δ[K+] from the ECG. The aim of this study is to derive polynomial models relating dw and dw,c with Δ[K+], and to test their ability to reliably sense and quantify Δ[K+] values. Methods: 48-hour Holter ECGs and [K+] values from six blood samples were collected from 29 ESRD-HD patients. For every patient, dw and dw,c were computed, and linear, quadratic, and cubic fitting models were derived from them. Then, Spearman’s (ρ) and Pearson’s (r) correlation coefficients, and the estimation error (ed) between Δ[K+] and the corresponding model-estimated values (Δ^[K+]) were calculated. Results and Discussions: Nonlinear models were the most suitable for Δ[K+] estimation, rendering higher Pearson’s correlation (median 0.77 ≤r≤ 0.92) and smaller estimation error (median 0.20 ≤ed≤ 0.43) than the linear model (median 0.76 ≤r≤ 0.86 and 0.30 ≤ed≤ 0.40), even if similar Spearman’s ρ were found across models (median 0.77 ≤ρ≤ 0.83). Conclusion: Results support the use of nonlinear T-wave-based models as Δ[K+] sensors in ESRD-HD patients.

Monitoring blood potassium concentration in hemodialysis patients by quantifying T-wave morphology dynamics

We investigated the ability of time-warping-based ECG-derived markers of T-wave morphology changes in time (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d_{w}$$\end{document}dw) and amplitude (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d_a$$\end{document}da), as well as their non-linear components (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${d_w^{{\mathrm{NL}}}}$$\end{document}dwNL and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${d_a^{\mathrm{NL}}}$$\end{document}daNL), and the heart rate corrected counterpart (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d_{w,c}$$\end{document}dw,c), to monitor potassium concentration (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$[K^{+}]$$\end{document}[K+]) changes (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta [K^+]$$\end{document}Δ[K+]) in end-stage renal disease (ESRD) patients undergoing hemodialysis (HD). We compared the performance of the proposed time-warping markers, together with other previously proposed \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$[K^{+}]$$\end{document}[K+] markers, such as T-wave width (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_w$$\end{document}Tw) and T-wave slope-to-amplitude ratio (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_{S/A}$$\end{document}TS/A), when computed from standard ECG leads as well as from principal component analysis (PCA)-based leads. 48-hour ECG recordings and a set of hourly-collected blood samples from 29 ESRD-HD patients were acquired. Values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d_w$$\end{document}dw, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d_a$$\end{document}da, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${d_w^{\mathrm{NL}}}$$\end{document}dwNL, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${d_a^{\mathrm{NL}}}$$\end{document}daNL and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d_{w,c}$$\end{document}dw,c were calculated by comparing the morphology of the mean warped T-waves (MWTWs) derived at each hour along the HD with that from a reference MWTW, measured at the end of the HD. From the same MWTWs \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_w$$\end{document}Tw and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_{S/A}$$\end{document}TS/A were also extracted. Similarly, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta [K^+]$$\end{document}Δ[K+] was calculated as the difference between the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$[K^{+}]$$\end{document}[K+] values at each hour and the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$[K^{+}]$$\end{document}[K+] reference level at the end of the HD session. We found that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d_{w}$$\end{document}dw and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d_{w,c}$$\end{document}dw,c showed higher correlation coefficients with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta [K^+]$$\end{document}Δ[K+] than \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_{S/A}$$\end{document}TS/A—Spearman’s (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho$$\end{document}ρ) and Pearson’s (r)—and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_w$$\end{document}Tw—Spearman’s (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho$$\end{document}ρ)—in both SL and PCA approaches being the intra-patient median \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho \ge 0.82$$\end{document}ρ≥0.82 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r \ge 0.87$$\end{document}r≥0.87 in SL and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho \ge 0.82$$\end{document}ρ≥0.82 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r \ge 0.89$$\end{document}r≥0.89 in PCA respectively. Our findings would point at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d_{w}$$\end{document}dw and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d_{w,c}$$\end{document}dw,c as the most suitable surrogate of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta [K^+]$$\end{document}Δ[K+], suggesting that they could be potentially useful for non-invasive monitoring of ESRD-HD patients in hospital, as well as in ambulatory settings. Therefore, the tracking of T-wave morphology variations by means of time-warping analysis could improve continuous and remote \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$[K^{+}]$$\end{document}[K+] monitoring of ESRD-HD patients and flagging risk of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$[K^{+}]$$\end{document}[K+]-related cardiovascular events.

Breathing rate and heart rate as confounding factors in measuring T wave alternans and morphological variability in ECG

OBJECTIVE

High morphological variability magnitude (MVM) and microvolt T wave alternans (TWA) within an electrocardiogram (ECG) signifies increased electrical instability and risk of sudden cardiac death. However, the influence of breathing rate (BR), heart rate (HR), and signal-to-noise ratio (SNR) is unknown and may inflate measured values.

APPROACH

We synthesize ECGs with morphologies derived from the Physikalisch-Technische Bundesanstalt Database. We calculate MVM and TWA at varying BRs, HRs and SNRs. We compare the MVM and TWA of signal with versus without breathing at varying HRs and SNRs. We then quantify the percentage of MVM and TWA estimates affected by BR and HR in a healthy population and assess the effect of removing these affected estimates on a method for classifying individuals with and without post-traumatic stress disorder (PTSD).

MAIN RESULTS

For signals with high SNR (>15 dB), MVM is significantly increased when BRs are > 9 respirations/minute (rpm) and HRs are < 100 beats/minute (bpm). Increased TWAs are detected for HR/BR pairs of 60/15, 60/30 and 120/30 bpm/rpm. For 18 healthy participants, 8.33% of TWA windows and 66.76% of MVM windows are affected by BR and HR. On average, the number of windows with TWA elevations > 47 μV decreases by 23% after excluding regions with significant BR and HR effect. Adding HR and BR to a morphological variability feature increases the classification performance by 6% for individuals with and without PTSD.

SIGNIFICANCE

Physiological BR and HR significantly increase MVM and TWA , indicating that BR and HR should be considered separately as confounders. The code for this work has been released as part of an open-source toolbox.

A Method to Minimise the Impact of ECG Marker Inaccuracies on the Spatial QRS-T angle: Evaluation on 1,512 Manually Annotated ECGs

•

Inaccuracies of QRS and T-wave markers significantly impact QRS-Ta estimation.

•

These errors influence the classification of clinically relevant abnormal values.

•

Our algorithm provides robust measurements in the presence of inaccurate VCG markers.

•

We present for the first time, the distribution of the QRS-Ta in a large cohort.

The spatial QRS-T angle (QRS-Ta) derived from the vectorcardiogram (VCG) is a strong risk predictor for ventricular arrhythmia and sudden cardiac death with potential use for mass screening. Accurate QRS-Ta estimation in the presence of ECG delineation errors is crucial for its deployment as a prognostic test. Our study assessed the effect of inaccurate QRS and T-wave marker placement on QRS-Ta estimation and proposes a robust method for its calculation.

Reference QRS-Ta measurements were derived from 1,512 VCGs manually annotated by three expert reviewers. We systematically changed onset and offset timings of QRS and T-wave markers to simulate inaccurate placement. The QRS-Ta was recalculated using a standard approach and our proposed algorithm, which limits the impact of VCG marker inaccuracies by defining the vector origin as an interval preceding QRS-onset and redefines the beginning and end of QRS and T-wave loops.

Using the standard approach, mean absolute errors (MAE) in peak QRS-Ta were >40% and sensitivity and precision in the detection of abnormality (>105°) were <80% and <65% respectively, when QRS-onset was delayed or QRS-offset anticipated >15 ms. Using our proposed algorithm, MAE for peak QRS-Ta were reduced to <4% and sensitivity and precision of abnormality were >94% for inaccuracies up to ±15 ms. Similar results were obtained for mean QRS-Ta.

In conclusion, inaccuracies of QRS and T-wave markers can significantly influence the QRS-Ta. Our proposed algorithm provides robust QRS-Ta measurements in the presence of inaccurate VCG annotation, enabling its use in large datasets.

Characterization of T Wave Amplitude, Duration and Morphology Changes During Hemodialysis: Relationship With Serum Electrolyte Levels and Heart Rate

OBJECTIVE

Chronic kidney disease affects more than 10% of the world population. Changes in serum ion concentrations increase the risk for ventricular arrhythmias and sudden cardiac death, particularly in end-stage renal disease (ESRD) patients. We characterized how T wave amplitude, duration and morphology descriptors change with variations in serum levels of potassium and calcium and in heart rate, both in ESRD patients and in simulated ventricular fibers.

METHODS

Electrocardiogram (ECG) recordings from twenty ESRD patients undergoing hemodialysis (HD) and pseudo-ECGs (pECGs) calculated from twenty-two simulated ventricular fibers at varying transmural heterogeneity levels were processed to quantify T wave width ( T_w), T wave slope-to-amplitude ratio ([Formula: see text]) and four indices of T wave morphological variability based on time warping ( d_w, [Formula: see text], d_a and [Formula: see text]). Serum potassium and calcium levels and heart rate were measured along HD.

RESULTS

[Formula: see text] was the marker most strongly correlated with serum potassium, dw with calcium and da with heart rate, after correction for covariates. Median values of partial correlation coefficients were 0.75, -0.74 and -0.90, respectively. For all analyzed T wave descriptors, high inter-patient variability was observed in the pattern of such relationships. This variability, accentuated during the first HD time points, was reproduced in the simulations and shown to be influenced by differences in transmural heterogeneity.

CONCLUSION

Changes in serum potassium and calcium levels and in heart rate strongly affect T wave descriptors, particularly those quantifying morphological variability.

SIGNIFICANCE

ECG markers have the potential to be used for monitoring serum ion concentrations in ESRD patients.

Hypoxia-induced Effects on ECG Depolarization by Time Warping Analysis during Recurrent Obstructive Apnea

In this work, we evaluated a non-linear approach to estimate morphological variations in ECG depolarization, in the context of intermittent hypoxia (IH). Obstructive apnea sequences were provoked for 15 minutes in anesthetized Sprague-Dawley rats, alternating with equal periods of normal breathing, in a recurrent obstructive sleep apnea (OSA) model. Each apnea episode lasted 15 s, while the frequency used for each sequence was randomly selected. Average heartbeats obtained before the start and at the end of each episode, were delineated to extract only the QRS wave. Then, the segmented QRS waves were non-linearly aligned using the dynamic time warping (DWT) algorithm. Morphological QRS changes in both the amplitude and temporal domains were estimated from this alignment procedure. The hypoxic and basal segments were analyzed using ECG (lead I) recordings acquired during the experiment. To assess the effects of IH over time, the changes relative to the basal QRS wave were determined, in the intervals prior to each successive events until the end of the experiment. The results showed a progressive increase in the amplitude and time-domain morphological markers of the QRS wave along the experiment, which were strongly correlated with the changes in traditional QRS markers (r ≈ 0.9). Significant changes were found between pre-apnea and hypoxic measures only for the time-domain analysis (p<0.001), probably due to the short duration of the simulated apnea episodes.Clinical relevance Increased variability in ECG depolarization morphology during recurrent hypoxic episodes would be closely related to the expression of cardiovascular dysfunction in OSA patients.

Long-Term Microgravity Exposure Increases ECG Repolarization Instability Manifested by Low-Frequency Oscillations of T-Wave Vector

Ventricular arrhythmias and sudden cardiac death during long-term space missions are a major concern for space agencies. Long-duration spaceflight and its ground-based analog head-down bed rest (HDBR) have been reported to markedly alter autonomic and cardiac functioning, particularly affecting ventricular repolarization of the electrocardiogram (ECG). In this study, novel methods are developed, departing from previously published methodologies, to quantify the index of Periodic Repolarization Dynamics (PRD), an arrhythmic risk marker that characterizes sympathetically-mediated low-frequency oscillations in the T-wave vector. PRD is evaluated in ECGs from 42 volunteers at rest and during an orthostatic tilt table test recorded before and after 60-day –6° HDBR. Our results indicate that tilt test, on top of enhancing sympathetic regulation of heart rate, notably increases PRD, both before and after HDBR, thus supporting previous evidence on PRD being an indicator of sympathetic modulation of ventricular repolarization. Importantly, long-term microgravity exposure is shown to lead to significant increases in PRD, both when evaluated at rest and, even more notably, in response to tilt test. The extent of microgravity-induced changes in PRD has been associated with arrhythmic risk in prior studies. An exercise-based, but not a nutrition-based, countermeasure is able to partially reverse microgravity-induced effects on PRD. In conclusion, long-term exposure to microgravity conditions leads to elevated low-frequency oscillations of ventricular repolarization, which are potentiated following sympathetic stimulation and are related to increased risk for repolarization instabilities and arrhythmias. Tested countermeasures are only partially effective in counteracting microgravity effects.

Cardiovascular Predictive Value and Genetic Basis of Ventricular Repolarization Dynamics

BACKGROUND

Early prediction of cardiovascular risk in the general population remains an important issue. The T-wave morphology restitution (TMR), an ECG marker quantifying ventricular repolarization dynamics, is strongly associated with cardiovascular mortality in patients with heart failure. Our aim was to evaluate the cardiovascular prognostic value of TMR in a UK middle-aged population and identify any genetic contribution.

METHODS

We analyzed ECG recordings from 55 222 individuals from a UK middle-aged population undergoing an exercise stress test in UK Biobank (UKB). TMR was used to measure ventricular repolarization dynamics, exposed in this cohort by exercise (TMR during exercise, TMR^ex) and recovery from exercise (TMR during recovery, TMR^rec). The primary end point was cardiovascular events; secondary end points were all-cause mortality, ventricular arrhythmias, and atrial fibrillation with median follow-up of 7 years. Genome-wide association studies for TMR^ex and TMR^rec were performed, and genetic risk scores were derived and tested for association in independent samples from the full UKB cohort (N=360 631).

RESULTS

A total of 1743 (3.2%) individuals in UKB who underwent the exercise stress test had a cardiovascular event, and TMR^rec was significantly associated with cardiovascular events (hazard ratio, 1.11; P=5×10^-7), independent of clinical variables and other ECG markers. TMR^rec was also associated with all-cause mortality (hazard ratio, 1.10) and ventricular arrhythmias (hazard ratio, 1.16). We identified 12 genetic loci in total for TMR^ex and TMR^rec, of which 9 are associated with another ECG marker. Individuals in the top 20% of the TMR^rec genetic risk score were significantly more likely to have a cardiovascular event in the full UKB cohort (18 997, 5.3%) than individuals in the bottom 20% (hazard ratio, 1.07; P=6×10^-3).

CONCLUSIONS

TMR and TMR genetic risk scores are significantly associated with cardiovascular risk in a UK middle-aged population, supporting the hypothesis that increased spatio-temporal heterogeneity of ventricular repolarization is a substrate for cardiovascular risk and the validity of TMR as a cardiovascular risk predictor.

Mechanisms Underlying Interactions Between Low-Frequency Oscillations and Beat-to-Beat Variability of Celullar Ventricular Repolarization in Response to Sympathetic Stimulation: Implications for Arrhythmogenesis

Background and Objectives: Enhanced beat-to-beat variability of ventricular repolarization (BVR) has been linked to arrhythmias and sudden cardiac death. Recent experimental studies on human left ventricular epicardial electrograms have shown that BVR closely interacts with low-frequency (LF) oscillations of activation recovery interval during sympathetic provocation. In this work human ventricular computational cell models are developed to reproduce the experimentally observed interactions between BVR and its LF oscillations, to assess underlying mechanisms and to establish a relationship with arrhythmic risk.

Materials and Methods: A set of human ventricular action potential (AP) models covering a range of experimental electrophysiological characteristics was constructed. These models incorporated stochasticity in major ionic currents as well as descriptions of β-adrenergic stimulation and mechanical effects to investigate the AP response to enhanced sympathetic activity. Statistical methods based on Automatic Relevance Determination and Canonical Correlation Analysis were developed to unravel individual and common factors contributing to BVR and LF patterning of APD in response to sympathetic provocation.

Results: Simulated results reproduced experimental evidences on the interactions between BVR and LF oscillations of AP duration (APD), with replication of the high inter-individual variability observed in both phenomena. I_CaL, I_Kr and I_K1 currents were identified as common ionic modulators of the inter-individual differences in BVR and LF oscillatory behavior and were shown to be crucial in determining susceptibility to arrhythmogenic events.

Conclusions: The calibrated family of human ventricular cell models proposed in this study allows reproducing experimentally reported interactions between BVR and LF oscillations of APD. Ionic factors involving I_CaL, I_Kr and I_K1 currents are found to underlie correlated increments in both phenomena in response to sympathetic provocation. A link to arrhythmogenesis is established for concomitantly elevated levels of BVR and its LF oscillations.

Sudden cardiac death and pump failure death prediction in chronic heart failure by combining ECG and clinical markers in an integrated risk model

Background

Sudden cardiac death (SCD) and pump failure death (PFD) are common endpoints in chronic heart failure (CHF) patients, but prevention strategies are different. Currently used tools to specifically predict these endpoints are limited. We developed risk models to specifically assess SCD and PFD risk in CHF by combining ECG markers and clinical variables.

Methods

The relation of clinical and ECG markers with SCD and PFD risk was assessed in 597 patients enrolled in the MUSIC (MUerte Súbita en Insuficiencia Cardiaca) study. ECG indices included: turbulence slope (TS), reflecting autonomic dysfunction; T-wave alternans (TWA), reflecting ventricular repolarization instability; and T-peak-to-end restitution (Δα^Tpe) and T-wave morphology restitution (TMR), both reflecting changes in dispersion of repolarization due to heart rate changes. Standard clinical indices were also included.

Results

The indices with the greatest SCD prognostic impact were gender, New York Heart Association (NYHA) class, left ventricular ejection fraction, TWA, Δα^Tpe and TMR. For PFD, the indices were diabetes, NYHA class, Δα^Tpe and TS. Using a model with only clinical variables, the hazard ratios (HRs) for SCD and PFD for patients in the high-risk group (fifth quintile of risk score) with respect to patients in the low-risk group (first and second quintiles of risk score) were both greater than 4. HRs for SCD and PFD increased to 9 and 11 when using a model including only ECG markers, and to 14 and 13, when combining clinical and ECG markers.

Conclusion

The inclusion of ECG markers capturing complementary pro-arrhythmic and pump failure mechanisms into risk models based only on standard clinical variables substantially improves prediction of SCD and PFD in CHF patients.

Long-term intra-individual reproducibility of heart rate dynamics during exercise and recovery in the UK Biobank cohort

Background

The heart rate (HR) response to exercise provides useful information about the autonomic function and has prognostic value, but its reproducibility over a long period of time, a critical requirement for using it as a clinical biomarker, is undetermined.

Aim

To determine the intra-individual reproducibility of HR dynamics during sub-maximum exercise and one minute recovery.

Methods

1187 individuals from the Cardio physical fitness assessment test of the UK Biobank repeated a standard exercise stress test twice (recall time 34.2 ± 2.8 months) and were prospectively studied.

Results

821 individuals complied with inclusion criteria for reproducibility analysis, including peak workload differences between assessments ≤10 W. Intra-individual correlation between HR profile during the first and the second assessment was very high and higher than inter-individual correlation (0.92±0.08 vs 0.87±0.11, p<0.01). Intra-individual correlation of indices describing HR dynamics was: ρ = 0.81 for maximum HR during exercise; ρ = 0.71 for minimum HR during recovery; ρ = 0.70 for HR changes during both exercise and recovery; Intra-individual correlation was higher for these indices of HR dynamics than for resting HR (ρ = 0.64). Bland-Altman plots demonstrated good agreement between HR indices estimated during the first and second assessment. A small but consistent bias was registered for all repeated measurements. The intra-individual consistency of abnormal values was about 60–70%.

Conclusions

The HR dynamics during exercise and recovery are reproducible over a period of 3 years, with moderate to strong intra-individual reproducibility of abnormal values.

T-Wave Morphology Restitution Predicts Sudden Cardiac Death in Patients With Chronic Heart Failure

BACKGROUND

Patients with chronic heart failure are at high risk of sudden cardiac death (SCD). Increased dispersion of repolarization restitution has been associated with SCD, and we hypothesize that this should be reflected in the morphology of the T-wave and its variations with heart rate. The aim of this study is to propose an electrocardiogram (ECG)-based index characterizing T-wave morphology restitution (TMR), and to assess its association with SCD risk in a population of chronic heart failure patients.

METHODS AND RESULTS

Holter ECGs from 651 ambulatory patients with chronic heart failure from the MUSIC (MUerte Súbita en Insuficiencia Cardiaca) study were available for the analysis. TMR was quantified by measuring the morphological variation of the T-wave per RR increment using time-warping metrics, and its predictive power was compared to that of clinical variables such as the left ventricular ejection fraction and other ECG-derived indices, such as T-wave alternans and heart rate variability. TMR was significantly higher in SCD victims than in the rest of patients (median 0.046 versus 0.039, P<0.001). When TMR was dichotomized at TMR=0.040, the SCD rate was significantly higher in the TMR≥0.040 group (P<0.001). Cox analysis revealed that TMR≥0.040 was strongly associated with SCD, with a hazard ratio of 3.27 (P<0.001), independently of clinical and ECG-derived variables. No association was found between TMR and pump failure death.

CONCLUSIONS

This study shows that TMR is specifically associated with SCD in a population of chronic heart failure patients, and it is a better predictor than clinical and ECG-derived variables.