In Vivo Validation of Electrocardiographic Imaging

Evaluation of five methods for the interpolation of bad leads in the solution of the inverse electrocardiography problem

Objective.This study aims to assess the sensitivity of epicardial potential-based electrocardiographic imaging (ECGI) to the removal or interpolation of bad leads.Approach.We utilized experimental data from two distinct centers. Langendorff-perfused pig (n= 2) and dog (n= 2) hearts were suspended in a human torso-shaped tank and paced from the ventricles. Six different bad lead configurations were designed based on clinical experience. Five interpolation methods were applied to estimate the missing data. Zero-order Tikhonov regularization was used to solve the inverse problem for complete data, data with removed bad leads, and interpolated data. We assessed the quality of interpolated ECG signals and ECGI reconstructions using several metrics, comparing the performance of interpolation methods and the impact of bad lead removal versus interpolation on ECGI.Main results.The performance of ECG interpolation strongly correlated with ECGI reconstruction. The hybrid method exhibited the best performance among interpolation techniques, followed closely by the inverse-forward and Kriging methods. Bad leads located over high amplitude/high gradient areas on the torso significantly impacted ECGI reconstructions, even with minor interpolation errors. The choice between removing or interpolating bad leads depends on the location of missing leads and confidence in interpolation performance. If uncertainty exists, removing bad leads is the safer option, particularly when they are positioned in high amplitude/high gradient regions. In instances where interpolation is necessary, the inverse-forward and Kriging methods, which do not require training, are recommended.Significance.This study represents the first comprehensive evaluation of the advantages and drawbacks of interpolating versus removing bad leads in the context of ECGI, providing valuable insights into ECGI performance.

The arrhythmic substrate of hypertrophic cardiomyopathy using ECG imaging

Introduction: Patients with hypertrophic cardiomyopathy (HCM) are at risk for lethal ventricular arrhythmia, but the electrophysiological substrate behind this is not well-understood. We used non-invasive electrocardiographic imaging to characterize patients with HCM, including cardiac arrest survivors. Methods: HCM patients surviving ventricular fibrillation or hemodynamically unstable ventricular tachycardia (n = 17) were compared to HCM patients without a personal history of potentially lethal arrhythmia (n = 20) and a pooled control group with structurally normal hearts. Subjects underwent exercise testing by non-invasive electrocardiographic imaging to estimate epicardial electrophysiology. Results: Visual inspection of reconstructed epicardial HCM maps revealed isolated patches of late activation time (AT), prolonged activation-recovery intervals (ARIs), as well as reversal of apico-basal trends in T-wave inversion and ARI compared to controls (p < 0.005 for all). AT and ARI were compared between groups. The pooled HCM group had longer mean AT (60.1 ms vs. 52.2 ms, p < 0.001), activation dispersion (55.2 ms vs. 48.6 ms, p = 0.026), and mean ARI (227 ms vs. 217 ms, p = 0.016) than structurally normal heart controls. HCM ventricular arrhythmia survivors could be differentiated from HCM patients without a personal history of life-threatening arrhythmia by longer mean AT (63.2 ms vs. 57.4 ms, p = 0.007), steeper activation gradients (0.45 ms/mm vs. 0.36 ms/mm, p = 0.011), and longer mean ARI (234.0 ms vs. 221.4 ms, p = 0.026). A logistic regression model including whole heart mean activation time and activation recovery interval could identify ventricular arrhythmia survivors from the HCM cohort, producing a C statistic of 0.76 (95% confidence interval 0.72-0.81), with an optimal sensitivity of 78.6% and a specificity of 79.8%. Discussion: The HCM epicardial electrotype is characterized by delayed, dispersed conduction and prolonged, dispersed activation-recovery intervals. Combination of electrophysiologic measures with logistic regression can improve differentiation over single variables. Future studies could test such models prospectively for risk stratification of sudden death due to HCM.

Hybrid Neural State-Space Modeling for Supervised and Unsupervised Electrocardiographic Imaging

State-space modeling (SSM) provides a general framework for many image reconstruction tasks. Error in a priori physiological knowledge of the imaging physics, can bring incorrectness to solutions. Modern deep-learning approaches show great promise but lack interpretability and rely on large amounts of labeled data. In this paper, we present a novel hybrid SSM framework for electrocardiographic imaging (ECGI) to leverage the advantage of state-space formulations in data-driven learning. We first leverage the physics-based forward operator to supervise the learning. We then introduce neural modeling of the transition function and the associated Bayesian filtering strategy. We applied the hybrid SSM framework to reconstruct electrical activity on the heart surface from body-surface potentials. In unsupervised settings of both in-silico and in-vivo data without cardiac electrical activity as the ground truth to supervise the learning, we demonstrated improved ECGI performances of the hybrid SSM framework trained from a small number of ECG observations in comparison to the fixed SSM. We further demonstrated that, when in-silico simulation data becomes available, mixed supervised and unsupervised training of the hybrid SSM achieved a further 40.6% and 45.6% improvements, respectively, in comparison to traditional ECGI baselines and supervised data-driven ECGI baselines for localizing the origin of ventricular activations in real data.

Variant patterns of electrical activation and recovery in normal human hearts revealed by noninvasive electrocardiographic imaging

AIMS

Although electrical activity of the normal human heart is well characterized by the electrocardiogram, detailed insights into within-subject and between-subject variations of ventricular activation and recovery by noninvasive electroanatomic mapping are lacking. We characterized human epicardial activation and recovery within and between normal subjects using non-invasive electrocardiographic imaging (ECGI) as a basis to better understand pathology.

METHODS AND RESULTS

Epicardial activation and recovery were assessed by ECGI in 22 normal subjects, 4 subjects with bundle branch block (BBB) and 4 with long-QT syndrome (LQTS). We compared characteristics between the ventricles [left ventricle (LV) and right ventricle (RV)], sexes, and age groups (<50/≥50years). Pearson's correlation coefficient (CC) was used for within-subject and between-subject comparisons. Age of normal subjects averaged 49 ± 14 years, 6/22 were male, and no structural/electrical heart disease was present. The average activation time was longer in LV than in RV, but not different by sex or age. Electrical recovery was similar for the ventricles, but started earlier and was on average shorter in males. Median CCs of between-subject comparisons of the ECG signals, activation, and recovery patterns were 0.61, 0.32, and 0.19, respectively. Within-subject beat-to-beat comparisons yielded higher CCs (0.98, 0.89, and 0.82, respectively). Activation and/or recovery patterns of patients with BBB or LQTS contrasted significantly with those found in the normal population.

CONCLUSION

Activation and recovery patterns vary profoundly between normal subjects, but are stable individually beat to beat, with a male preponderance to shorter recovery. Individual characterization by ECGI at baseline serves as reference to better understand the emergence, progression, and treatment of electrical heart disease.

High Density Body Surface Potential Mapping with Conducting Polymer-Eutectogel Electrode Arrays for ECG imaging

Electrocardiography imaging (ECGi) is a non-invasive inverse reconstruction procedure which employs body surface potential maps (BSPM) obtained from surface electrode array measurements to improve the spatial resolution and interpretability of conventional electrocardiography (ECG) for the diagnosis of cardiac dysfunction. ECGi currently lacks precision, which has prevented its adoption in clinical setups. The introduction of high-density electrode arrays could increase ECGi reconstruction accuracy but is not attempted before due to manufacturing and processing limitations. Advances in multiple fields have now enabled the implementation of such arrays which poses questions on optimal array design parameters for ECGi. In this work, a novel conducting polymer electrode manufacturing process on flexible substrates is proposed to achieve high-density, mm-sized, conformable, long-term, and easily attachable electrode arrays for BSPM with parameters optimally selected for ECGi applications. Temporal, spectral, and correlation analysis are performed on a prototype array demonstrating the validity of the chosen parameters and the feasibility of high-density BSPM, paving the way for ECGi devices fit for clinical application.

The Dutch Idiopathic Ventricular Fibrillation Registry: progress report on the quest to identify the unidentifiable

BACKGROUND

Idiopathic ventricular fibrillation (iVF) is a rare cause of sudden cardiac arrest and, by definition, a diagnosis of exclusion. Due to the rarity of the disease, previous and current studies are limited by their retrospective design and small patient numbers. Even though the incidence of iVF has declined owing to the identification of new disease entities, an important subgroup of patients remains.

AIM

To expand the existing Dutch iVF Registry into a large nationwide cohort of patients initially diagnosed with iVF, to reveal the underlying cause of iVF in these patients, and to improve arrhythmia management.

METHODS

The Dutch iVF Registry includes sudden cardiac arrest survivors with an initial diagnosis of iVF. Clinical data and outcomes are collected. Outcomes include subsequent detection of a diagnosis other than 'idiopathic', arrhythmia recurrence and death. Non-invasive electrocardiographic imaging is used to investigate electropathological substrates and triggers of VF.

RESULTS

To date, 432 patients have been included in the registry (median age at event 40 years (interquartile range 28-52)), 61% male. During a median follow-up of 6 (2-12) years, 38 patients (9%) received a diagnosis other than 'idiopathic'. Eleven iVF patients were characterised with electrocardiographic imaging.

CONCLUSION

The Dutch iVF Registry is currently the largest of its kind worldwide. In this heterogeneous population of index patients, we aim to identify common functional denominators associated with iVF. With the implementation of non-invasive electrocardiographic imaging and other diagnostic modalities (e.g. echocardiographic deformation, cardiac magnetic resonance), we advance the possibilities to reveal pro-fibrillatory substrates.

Standardized 2D atrial mapping and its clinical applications

The visualization and comparison of electrophysiological information in the atrium among different patients could be facilitated by a standardized 2D atrial mapping. However, due to the complexity of the atrial anatomy, unfolding the 3D geometry into a 2D atrial mapping is challenging. In this study, we aim to develop a standardized approach to achieve a 2D atrial mapping that connects the left and right atria, while maintaining fixed positions and sizes of atrial segments across individuals. Atrial segmentation is a prerequisite for the process. Segmentation includes 19 different segments with 12 segments from the left atrium, 5 segments from the right atrium, and two segments for the atrial septum. To ensure consistent and physiologically meaningful segment connections, an automated procedure is applied to open up the atrial surfaces and project the 3D information into 2D. The corresponding 2D atrial mapping can then be utilized to visualize different electrophysiological information of a patient, such as activation time patterns or phase maps. This can in turn provide useful information for guiding catheter ablation. The proposed standardized 2D maps can also be used to compare more easily structural information like fibrosis distribution with rotor presence and location. We show several examples of visualization of different electrophysiological properties for both healthy subjects and patients affected by atrial fibrillation. These examples show that the proposed maps provide an easy way to visualize and interpret intra-subject information and perform inter-subject comparison, which may provide a reference framework for the analysis of the atrial fibrillation substrate before treatment, and during a catheter ablation procedure.

Basis and applicability of noninvasive inverse electrocardiography: a comparison between cardiac source models

The body surface electrocardiogram (ECG) is a direct result of electrical activity generated by the myocardium. Using the body surface ECGs to reconstruct cardiac electrical activity is called the inverse problem of electrocardiography. The method to solve the inverse problem depends on the chosen cardiac source model to describe cardiac electrical activity. In this paper, we describe the theoretical basis of two inverse methods based on the most commonly used cardiac source models: the epicardial potential model and the equivalent dipole layer model. We discuss similarities and differences in applicability, strengths and weaknesses and sketch a road towards improved inverse solutions by targeted use, sequential application or a combination of the two methods.

Technical development and feasibility of a reusable vest to integrate cardiovascular magnetic resonance with electrocardiographic imaging

BACKGROUND

Electrocardiographic imaging (ECGI) generates electrophysiological (EP) biomarkers while cardiovascular magnetic resonance (CMR) imaging provides data about myocardial structure, function and tissue substrate. Combining this information in one examination is desirable but requires an affordable, reusable, and high-throughput solution. We therefore developed the CMR-ECGI vest and carried out this technical development study to assess its feasibility and repeatability in vivo.

METHODS

CMR was prospectively performed at 3T on participants after collecting surface potentials using the locally designed and fabricated 256-lead ECGI vest. Epicardial maps were reconstructed to generate local EP parameters such as activation time (AT), repolarization time (RT) and activation recovery intervals (ARI). 20 intra- and inter-observer and 8 scan re-scan repeatability tests.

RESULTS

77 participants were recruited: 27 young healthy volunteers (HV, 38.9 ± 8.5 years, 35% male) and 50 older persons (77.0 ± 0.1 years, 52% male). CMR-ECGI was achieved in all participants using the same reusable, washable vest without complications. Intra- and inter-observer variability was low (correlation coefficients [rs] across unipolar electrograms = 0.99 and 0.98 respectively) and scan re-scan repeatability was high (rs between 0.81 and 0.93). Compared to young HV, older persons had significantly longer RT (296.8 vs 289.3 ms, p = 0.002), ARI (249.8 vs 235.1 ms, p = 0.002) and local gradients of AT, RT and ARI (0.40 vs 0.34 ms/mm, p = 0,01; 0.92 vs 0.77 ms/mm, p = 0.03; and 1.12 vs 0.92 ms/mm, p = 0.01 respectively).

CONCLUSION

Our high-throughput CMR-ECGI solution is feasible and shows good reproducibility in younger and older participants. This new technology is now scalable for high throughput research to provide novel insights into arrhythmogenesis and potentially pave the way for more personalised risk stratification.

CLINICAL TRIAL REGISTRATION

Title: Multimorbidity Life-Course Approach to Myocardial Health-A Cardiac Sub-Study of the MRC National Survey of Health and Development (NSHD) (MyoFit46). National Clinical Trials (NCT) number: NCT05455125. URL: https://clinicaltrials.gov/ct2/show/NCT05455125?term=MyoFit&draw=2&rank=1.

Treating Stubborn Cardiac Arrhythmias-Looking Toward the Future

As animals can develop significant side effects or remain refractory while on antiarrhythmic medical therapy for tachyarrhythmias, interventional therapies are progressively being explored. This review will highlight the principles and utilities of implantable cardioverter-defibrillators, electrophysiological mapping and catheter ablation, three-dimensional electroanatomical mapping, and stereotactic arrhythmia radiotherapy. In particular, three-dimensional electroanatomical mapping is emerging as an adjunct electrophysiology tool to facilitate activation, substrate, and pace mapping for intuitive analysis of complex tachyarrhythmias. Unlike antiarrhythmic medications, these modalities offer potential for decreasing risk of sudden death and even permanent termination of tachyarrhythmias.

Innovations in ventricular tachycardia ablation

Catheter ablation of ventricular arrhythmias (VAs) has evolved significantly over the past decade and is currently a well-established therapeutic option. Technological advances and improved understanding of VA mechanisms have led to tremendous innovations in VA ablation. The purpose of this review article is to provide an overview of current innovations in VA ablation. Mapping techniques, such as ultra-high density mapping, isochronal late activation mapping, and ripple mapping, have provided improved arrhythmogenic substrate delineation and potential procedural success while limiting duration of ablation procedure and potential hemodynamic compromise. Besides, more advanced mapping and ablation techniques such as epicardial and intramyocardial ablation approaches have allowed operators to more precisely target arrhythmogenic substrate. Moreover, advances in alternate energy sources, such as electroporation, as well as stereotactic radiation therapy have been proposed to be effective and safe. New catheters, such as the lattice and the saline-enhanced radiofrequency catheters, have been designed to provide deeper and more durable tissue ablation lesions compared to conventional catheters. Contact force optimization and baseline impedance modulation are important tools to optimize VT radiofrequency ablation and improve procedural success. Furthermore, advances in cardiac imaging, specifically cardiac MRI, have great potential in identifying arrhythmogenic substrate and evaluating ablation success. Overall, VA ablation has undergone significant advances over the past years. Innovations in VA mapping techniques, alternate energy source, new catheters, and utilization of cardiac imaging have great potential to improve overall procedural safety, hemodynamic stability, and procedural success.

Comparison of dipole-based and potential-based ECGI methods for premature ventricular contraction beat localization with clinical data

Introduction: Localization of premature ventricular contraction (PVC) origin to guide the radiofrequency ablation (RFA) procedure is one of the prominent clinical goals of non-invasive electrocardiographic imaging. However, the results reported in the literature vary significantly depending on the source model and the level of complexity in the forward model. This study aims to compare the paced and spontaneous PVC localization performances of dipole-based and potential-based source models and corresponding inverse methods using the same clinical data and to evaluate the effects of torso inhomogeneities on these performances. Methods: The publicly available EP solution data from the EDGAR data repository (BSPs from a maximum of 240 electrodes) with known pacing locations and the Bratislava data (BSPs in 128 leads) with spontaneous PVCs from patients who underwent successful RFA procedures were used. Homogeneous and inhomogeneous torso models and corresponding forward problem solutions were used to relate sources on the closed epicardial and epicardial-endocardial surfaces. The localization error (LE) between the true and estimated pacing site/PVC origin was evaluated. Results: For paced data, the median LE values were 25.2 and 13.9 mm for the dipole-based and potential-based models, respectively. These median LE values were higher for the spontaneous PVC data: 30.2-33.0 mm for the dipole-based model and 28.9-39.2 mm for the potential-based model. The assumption of inhomogeneities in the torso model did not change the dipole-based solutions much, but using an inhomogeneous model improved the potential-based solutions on the epicardial-endocardial ventricular surface. Conclusion: For the specific task of localization of pacing site/PVC origin, the dipole-based source model is more stable and robust than the potential-based source model. The torso inhomogeneities affect the performances of PVC origin localization in each source model differently. Hence, care must be taken in generating patient-specific geometric and forward models depending on the source model representation used in electrocardiographic imaging (ECGI).

Investigation into the importance of using natural PVCs and pathological models for potential-based ECGI validation

Introduction: Premature ventricular contractions (PVCs) are one of the most commonly targeted pathologies for ECGI validation, often through ventricular stimulation to mimic the ectopic beat. However, it remains unclear if such stimulated beats faithfully reproduce spontaneously occurring PVCs, particularly in the case of the R-on-T phenomenon. The objective of this study was to determine the differences in ECGI accuracy when reconstructing spontaneous PVCs as compared to ventricular-stimulated beats and to explore the impact of pathophysiological perturbation on this reconstruction accuracy. Methods: Langendorff-perfused pig hearts (n = 3) were suspended in a human torso-shaped tank, and local hyperkalemia was induced through perfusion of a high-K+ solution (8 mM) into the LAD. Recordings were taken simultaneously from the heart and tank surfaces during ventricular pacing and during spontaneous PVCs (including R-on-T), both at baseline and high K+. Epicardial potentials were reconstructed from torso potentials using ECGI. Results: Spontaneously occurring PVCs were better reconstructed than stimulated beats at baseline in terms of electrogram morphology [correlation coefficient (CC) = 0.74 ± 0.05 vs. CC = 0.60 ± 0.10], potential maps (CC = 0.61 ± 0.06 vs. CC = 0.51 ± 0.12), and activation time maps (CC = 0.86 ± 0.07 vs. 0.76 ± 0.10), though there was no difference in the localization error (LE) of epicardial origin (LE = 14 ± 6 vs. 15 ± 11 mm). High K+ perfusion reduced the accuracy of ECGI reconstructions in terms of electrogram morphology (CC = 0.68 ± 0.10) and AT maps (CC = 0.70 ± 0.12 and 0.59 ± 0.23) for isolated PVCs and paced beats, respectively. LE trended worse, but the change was not significant (LE = 17 ± 9 and 20 ± 12 mm). Spontaneous PVCs were less well when the R-on-T phenomenon occurred and the activation wavefronts encountered a line of block. Conclusion: This study demonstrates the differences in ECGI accuracy between spontaneous PVCs and ventricular-paced beats. We also observed a reduction in this accuracy near regions of electrically inactive tissue. These results highlight the need for more physiologically realistic experimental models when evaluating the accuracy of ECGI methods. In particular, reconstruction accuracy needs to be further evaluated in the presence of R-on-T or isolated PVCs, particularly when encountering obstacles (functional or anatomical) which cause line of block and re-entry.

Influence of the Tikhonov Regularization Parameter on the Accuracy of the Inverse Problem in Electrocardiography

The electrocardiogram (ECG) is the standard method in clinical practice to non-invasively analyze the electrical activity of the heart, from electrodes placed on the body's surface. The ECG can provide a cardiologist with relevant information to assess the condition of the heart and the possible presence of cardiac pathology. Nonetheless, the global view of the heart's electrical activity given by the ECG cannot provide fully detailed and localized information about abnormal electrical propagation patterns and corresponding substrates on the surface of the heart. Electrocardiographic imaging, also known as the inverse problem in electrocardiography, tries to overcome these limitations by non-invasively reconstructing the heart surface potentials, starting from the corresponding body surface potentials, and the geometry of the torso and the heart. This problem is ill-posed, and regularization techniques are needed to achieve a stable and accurate solution. The standard approach is to use zero-order Tikhonov regularization and the L-curve approach to choose the optimal value for the regularization parameter. However, different methods have been proposed for computing the optimal value of the regularization parameter. Moreover, regardless of the estimation method used, this may still lead to over-regularization or under-regularization. In order to gain a better understanding of the effects of the choice of regularization parameter value, in this study, we first focused on the regularization parameter itself, and investigated its influence on the accuracy of the reconstruction of heart surface potentials, by assessing the reconstruction accuracy with high-precision simultaneous heart and torso recordings from four dogs. For this, we analyzed a sufficiently large range of parameter values. Secondly, we evaluated the performance of five different methods for the estimation of the regularization parameter, also in view of the results of the first analysis. Thirdly, we investigated the effect of using a fixed value of the regularization parameter across all reconstructed beats. Accuracy was measured in terms of the quality of reconstruction of the heart surface potentials and estimation of the activation and recovery times, when compared with ground truth recordings from the experimental dog data. Results show that values of the regularization parameter in the range (0.01-0.03) provide the best accuracy, and that the three best-performing estimation methods (L-Curve, Zero-Crossing, and CRESO) give values in this range. Moreover, a fixed value of the regularization parameter could achieve very similar performance to the beat-specific parameter values calculated by the different estimation methods. These findings are relevant as they suggest that regularization parameter estimation methods may provide the accurate reconstruction of heart surface potentials only for specific ranges of regularization parameter values, and that using a fixed value of the regularization parameter may represent a valid alternative, especially when computational efficiency or consistency across time is required.

Understanding repolarization in the intracardiac unipolar electrogram: A long-lasting controversy revisited

Background: The optimal way to determine repolarization time (RT) from the intracardiac unipolar electrogram (UEG) has been a topic of debate for decades. RT is typically determined by either the Wyatt method or the "alternative method," which both consider UEG T-wave slope, but differently. Objective: To determine the optimal method to measure RT on the UEG. Methods: Seven pig hearts surrounded by an epicardial sock with 100 electrodes were Langendorff-perfused with selective cannulation of the left anterior descending (LAD) coronary artery and submersed in a torso-shaped tank containing 256 electrodes on the torso surface. Repolarization was prolonged in the non-LAD-regions by infusing dofetilide and shortened in the LAD-region using pinacidil. RT was determined by the Wyatt (t_Wyatt) and alternative (t_Alt) methods, in both invasive (recorded with epicardial electrodes) and in non-invasive UEGs (reconstructed with electrocardiographic imaging). t_Wyatt and t_Alt were compared to local effective refractory period (ERP). Results: With contact mapping, mean absolute error (MAE) of t_Wyatt and t_Alt vs. ERP were 21 ms and 71 ms, respectively. Positive T-waves typically had an earlier ERP than negative T-waves, in line with theory. t_Wyatt -but not t_Alt-shortened by local infusion of pinacidil. Similar results were found for the non-invasive UEGs (MAE of t_Wyatt and t_Alt vs. ERP were 30 ms and 92 ms, respectively). Conclusion: The Wyatt method is the most accurate to determine RT from (non) invasive UEGs, based on novel and historical analyses. Using it to determine RT could unify and facilitate repolarization assessment and amplify its role in cardiac electrophysiology.

High-resolution structural-functional substrate-trigger characterization: Future roadmap for catheter ablation of ventricular tachycardia

Introduction

Patients with ventricular tachyarrhythmias (VT) are at high risk of sudden cardiac death. When appropriate, catheter ablation is modestly effective, with relatively high VT recurrence and complication rates. Personalized models that incorporate imaging and computational approaches have advanced VT management. However, 3D patient-specific functional electrical information is typically not considered. We hypothesize that incorporating non-invasive 3D electrical and structural characterization in a patient-specific model improves VT-substrate recognition and ablation targeting.

Materials and methods

In a 53-year-old male with ischemic cardiomyopathy and recurrent monomorphic VT, we built a structural-functional model based on high-resolution 3D late-gadolinium enhancement (LGE) cardiac magnetic resonance imaging (3D-LGE CMR), multi-detector computed tomography (CT), and electrocardiographic imaging (ECGI). Invasive data from high-density contact and pace mapping obtained during endocardial VT-substrate modification were also incorporated. The integrated 3D electro-anatomic model was analyzed off-line.

Results

Merging the invasive voltage maps and 3D-LGE CMR endocardial geometry led to a mean Euclidean node-to-node distance of 5 ± 2 mm. Inferolateral and apical areas of low bipolar voltage (<1.5 mV) were associated with high 3D-LGE CMR signal intensity (>0.4) and with higher transmurality of fibrosis. Areas of functional conduction delay or block (evoked delayed potentials, EDPs) were in close proximity to 3D-LGE CMR-derived heterogeneous tissue corridors. ECGI pinpointed the epicardial VT exit at ∼10 mm from the endocardial site of origin, both juxtaposed to the distal ends of two heterogeneous tissue corridors in the inferobasal left ventricle. Radiofrequency ablation at the entrances of these corridors, eliminating all EDPs, and at the VT site of origin rendered the patient non-inducible and arrhythmia-free until the present day (20 months follow-up). Off-line analysis in our model uncovered dynamic electrical instability of the LV inferolateral heterogeneous scar region which set the stage for an evolving VT circuit.

Discussion and conclusion

We developed a personalized 3D model that integrates high-resolution structural and electrical information and allows the investigation of their dynamic interaction during arrhythmia formation. This model enhances our mechanistic understanding of scar-related VT and provides an advanced, non-invasive roadmap for catheter ablation.

The circle of reentry: Characteristics of trigger-substrate interaction leading to sudden cardiac arrest

Sudden cardiac death is often caused by ventricular arrhythmias driven by reentry. Comprehensive characterization of the potential triggers and substrate in survivors of sudden cardiac arrest has provided insights into the trigger-substrate interaction leading to reentry. Previously, a "Triangle of Arrhythmogenesis", reflecting interactions between substrate, trigger and modulating factors, has been proposed to reason about arrhythmia initiation. Here, we expand upon this concept by separating the trigger and substrate characteristics in their spatial and temporal components. This yields four key elements that are required for the initiation of reentry: local dispersion of excitability (e.g., the presence of steep repolarization time gradients), a critical relative size of the region of excitability and the region of inexcitability (e.g., a sufficiently large region with early repolarization), a trigger that originates at a time when some tissue is excitable and other tissue is inexcitable (e.g., an early premature complex), and which occurs from an excitable region (e.g., from a region with early repolarization). We discuss how these findings yield a new mechanistic framework for reasoning about reentry initiation, the "Circle of Reentry." In a patient case of unexplained ventricular fibrillation, we then illustrate how a comprehensive clinical investigation of these trigger-substrate characteristics may help to understand the associated arrhythmia mechanism. We will also discuss how this reentry initiation concept may help to identify patients at risk, and how similar reasoning may apply to other reentrant arrhythmias.

Forward-Solution Noninvasive Computational Arrhythmia Mapping: The VMAP Study

BACKGROUND

The accuracy of noninvasive arrhythmia source localization using a forward-solution computational mapping system has not yet been evaluated in blinded, multicenter analysis. This study tested the hypothesis that a computational mapping system incorporating a comprehensive arrhythmia simulation library would provide accurate localization of the site-of-origin for atrial and ventricular arrhythmias and pacing using 12-lead ECG data when compared with the gold standard of invasive electrophysiology study and ablation.

METHODS

The VMAP study (Vectorcardiographic Mapping of Arrhythmogenic Probability) was a blinded, multicenter evaluation with final data analysis performed by an independent core laboratory. Eligible episodes included atrial and ventricular: tachycardia, fibrillation, pacing, premature atrial and ventricular complexes, and orthodromic atrioventricular reentrant tachycardia. Mapping system results were compared with the gold standard site of successful ablation or pacing during electrophysiology study and ablation. Mapping time was assessed from time-stamped logs. Prespecified performance goals were used for statistical comparisons.

RESULTS

A total of 255 episodes from 225 patients were enrolled from 4 centers. Regional accuracy for ventricular tachycardia and premature ventricular complexes in patients without significant structural heart disease (n=75, primary end point) was 98.7% (95% CI, 96.0%-100%; P<0.001 to reject predefined H0 <0.80). Regional accuracy for all episodes (secondary end point 1) was 96.9% (95% CI, 94.7%-99.0%; P<0.001 to reject predefined H0 <0.75). Accuracy for the exact or neighboring segment for all episodes (secondary end point 2) was 97.3% (95% CI, 95.2%-99.3%; P<0.001 to reject predefined H0 <0.70). Median spatial accuracy was 15 mm (n=255, interquartile range, 7-25 mm). The mapping process was completed in a median of 0.8 minutes (interquartile range, 0.4-1.4 minutes).

CONCLUSIONS

Computational ECG mapping using a forward-solution approach exceeded prespecified accuracy goals for arrhythmia and pacing localization. Spatial accuracy analysis demonstrated clinically actionable results. This rapid, noninvasive mapping technology may facilitate catheter-based and noninvasive targeted arrhythmia therapies.

REGISTRATION

URL: https://www.

CLINICALTRIALS

gov; Unique identifier: NCT04559061.

State-of-the-Art Multimodality Imaging in Sudden Cardiac Arrest with Focus on Idiopathic Ventricular Fibrillation: A Review

Idiopathic ventricular fibrillation is a rare cause of sudden cardiac arrest and a diagnosis by exclusion. Unraveling the mechanism of ventricular fibrillation is important for targeted management, and potentially for initiating family screening. Sudden cardiac arrest survivors undergo extensive clinical testing, with a growing role for multimodality imaging, before diagnosing "idiopathic" ventricular fibrillation. Multimodality imaging, considered as using multiple imaging modalities as diagnostics, is important for revealing structural myocardial abnormalities in patients with cardiac arrest. This review focuses on combining imaging modalities (echocardiography, cardiac magnetic resonance and computed tomography) and the electrocardiographic characterization of sudden cardiac arrest survivors and discusses the surplus value of multimodality imaging in the diagnostic routing of these patients. We focus on novel insights obtained through electrostructural and/or electromechanical imaging in apparently idiopathic ventricular fibrillation patients, with special attention to non-invasive electrocardiographic imaging.

A fast algorithm for spatiotemporal signals recovery using arbitrary dictionaries with application to electrocardiographic imaging

This paper presents a method to solve a linear regression problem subject to group lasso and ridge penalisation when the model has a Kronecker structure. This model was developed to solve the inverse problem of electrocardiography using sparse signal representation over a redundant dictionary or frame. The optimisation algorithm was performed using the block coordinate descent and proximal gradient descent methods. The explicit computation of the underlying Kronecker structure in the regression was avoided, reducing space and temporal complexity. We developed an algorithm that supports the use of arbitrary dictionaries to obtain solutions and allows a flexible group distribution.

Non-Contact Intracardiac Potential Mapping Using Mesh-Based and Meshless Inverse Solvers

Atrial fibrillation (AF) is the most common cardiac dysrhythmia and percutaneous catheter ablation is widely used to treat it. Panoramic mapping with multi-electrode catheters has been used to identify ablation targets in persistent AF but is limited by poor contact and inadequate coverage of the left atrial cavity. In this paper, we investigate the accuracy with which atrial endocardial surface potentials can be reconstructed from electrograms recorded with non-contact catheters. An in-silico approach was employed in which “ground-truth” surface potentials from experimental contact mapping studies and computer models were compared with inverse potential maps constructed by sampling the corresponding intracardiac field using virtual basket catheters. We demonstrate that it is possible to 1) specify the mixed boundary conditions required for mesh-based formulations of the potential inverse problem fully, and 2) reconstruct accurate inverse potential maps from recordings made with appropriately designed catheters. Accuracy improved when catheter dimensions were increased but was relatively stable when the catheter occupied >30% of atrial cavity volume. Independent of this, the capacity of non-contact catheters to resolve the complex atrial potential fields seen in reentrant atrial arrhythmia depended on the spatial distribution of electrodes on the surface bounding the catheter. Finally, we have shown that reliable inverse potential mapping is possible in near real-time with meshless methods that use the Method of Fundamental Solutions.

Intracardiac Inverse Potential Mapping Using the Method of Fundamental Solutions

Introduction: Atrial fibrillation (AF) is the most prevalent cardiac dysrhythmia and percutaneous catheter ablation is widely used to treat it. Panoramic mapping with multi-electrode catheters can identify ablation targets in persistent AF, but is limited by poor contact and inadequate coverage. Objective: To investigate the accuracy of inverse mapping of endocardial surface potentials from electrograms sampled with noncontact basket catheters. Methods: Our group has developed a computationally efficient inverse 3D mapping technique using a meshless method that employs the Method of Fundamental Solutions (MFS). An in-silico test bed was used to compare ground-truth surface potentials with corresponding inverse maps reconstructed from noncontact potentials sampled with virtual catheters. Ground-truth surface potentials were derived from high-density clinical contact mapping data and computer models. Results: Solutions of the intracardiac potential inverse problem with the MFS are robust, fast and accurate. Endocardial surface potentials can be faithfully reconstructed from noncontact recordings in real-time if the geometry of cardiac surface and the location of electrodes relative to it are known. Larger catheters with appropriate electrode density are needed to resolve complex reentrant atrial rhythms. Conclusion: Real-time panoramic potential mapping is feasible with noncontact intracardiac catheters using the MFS. Significance: Accurate endocardial potential maps can be reconstructed in AF with appropriately designed noncontact multi-electrode catheters.

Solving Inverse Electrocardiographic Mapping Using Machine Learning and Deep Learning Frameworks

Electrocardiographic imaging (ECGi) reconstructs electrograms at the heart's surface using the potentials recorded at the body's surface. This is called the inverse problem of electrocardiography. This study aimed to improve on the current solution methods using machine learning and deep learning frameworks. Electrocardiograms were simultaneously recorded from pigs' ventricles and their body surfaces. The Fully Connected Neural network (FCN), Long Short-term Memory (LSTM), Convolutional Neural Network (CNN) methods were used for constructing the model. A method is developed to align the data across different pigs. We evaluated the method using leave-one-out cross-validation. For the best result, the overall median of the correlation coefficient of the predicted ECG wave was 0.74. This study demonstrated that a neural network can be used to solve the inverse problem of ECGi with relatively small datasets, with an accuracy compatible with current standard methods.

Generalization and Regularization for Inverse Cardiac Estimators

Electrocardiographic Imaging (ECGI) aims to estimate the intracardiac potentials noninvasively, hence allowing the clinicians to better visualize and understand many arrhythmia mechanisms. Most of the estimators of epicardial potentials use a signal model based on an estimated spatial transfer matrix together with Tikhonov regularization techniques, which works well specially in simulations, but it can give limited accuracy in some real data. Based on the quasielectrostatic potential superposition principle, we propose a simple signal model that supports the implementation of principled out-of-sample algorithms for several of the most widely used regularization criteria in ECGI problems, hence improving the generalization capabilities of several of the current estimation methods. Experiments on simple cases (cylindrical and Gaussian shapes scrutinizing fast and slow changes, respectively) and on real data (examples of torso tank measurements available from Utah University, and an animal torso and epicardium measurements available from Maastricht University, both in the EDGAR public repository) show that the superposition-based out-of-sample tuning of regularization parameters promotes stabilized estimation errors of the unknown source potentials, while slightly increasing the re-estimation error on the measured data, as natural in non-overfitted solutions. The superposition signal model can be used for designing adequate out-of-sample tuning of Tikhonov regularization techniques, and it can be taken into account when using other regularization techniques in current commercial systems and research toolboxes on ECGI.

Spatiotemporal approximation of cardiac activation and recovery isochrones

BACKGROUND

The sequence of myocardial activation and recovery can be studied in detail by invasive catheter recordings of cardiac electrograms (EGMs), or noninvasive inverse reconstructions thereof with electrocardiographic imaging (ECGI). Local activation and recovery times are obtained from a unipolar EGM by the moment of maximum downslope of the QRS complex or maximum upslope of the T wave, respectively. However, both invasive and noninvasive recordings of intracardiac EGMs may suffer from noise and fractionation, making reliable detection of these deflections nontrivial.

METHODS

Here, we introduce a novel method that benefits from the spatial coupling of these processes, and incorporate not only the temporal EGM deflection, but also the spatial gradients. We validated this approach in computer simulations, in animal data with ECGI and invasive electrode recordings, and illustrated its use in a clinical case.

RESULTS

In the simulated data, the spatiotemporal approach was able to incorporate spatial information to better select the correct deflection in artificially fractionated EGMs and outperformed the traditional temporal-only method. In experimental data, the accuracy of time estimation from ECGI compared to invasive recordings significantly increased from R = 0.73 (activation) and R = 0.58 (recovery) with the temporal-only method to R = 0.79 (activation) and R = 0.72 (recovery) with the novel approach. Localization of the pacing origin of paced beats improved significantly from 36 mm mean error with the temporal-only approach to 23 mm with the spatiotemporal approach.

CONCLUSION

The spatiotemporal method to compute activation and recovery times from EGMs outperformed the traditional temporal-only approach in which spatial information was not taken into account.

Reducing Line-of-block Artifacts in Cardiac Activation Maps Estimated Using ECG Imaging: A Comparison of Source Models and Estimation Methods

OBJECTIVE

To investigate cardiac activation maps estimated using electrocardiographic imaging and to find methods reducing line-of-block (LoB) artifacts, while preserving real LoBs.

METHODS

Body surface potentials were computed for 137 simulated ventricular excitations. Subsequently, the inverse problem was solved to obtain extracellular potentials (EP) and transmembrane voltages (TMV). From these, activation times (AT) were estimated using four methods and compared to the ground truth. This process was evaluated with two cardiac mesh resolutions. Factors contributing to LoB artifacts were identified by analyzing the impact of spatial and temporal smoothing on the morphology of source signals.

RESULTS

AT estimation using a spatiotemporal derivative performed better than using a temporal derivative. Compared to deflection-based AT estimation, correlation-based methods were less prone to LoB artifacts but performed worse in identifying real LoBs. Temporal smoothing could eliminate artifacts for TMVs but not for EPs, which could be linked to their temporal morphology. TMVs led to more accurate ATs on the septum than EPs. Mesh resolution had a negligible effect on inverse reconstructions, but small distances were important for cross-correlation-based estimation of AT delays.

CONCLUSION

LoB artifacts are mainly caused by the inherent spatial smoothing effect of the inverse reconstruction. Among the configurations evaluated, only deflection-based AT estimation in combination with TMVs and strong temporal smoothing can prevent LoB artifacts, while preserving real LoBs.

SIGNIFICANCE

Regions of slow conduction are of considerable clinical interest and LoB artifacts observed in non-invasive ATs can lead to misinterpretations. We addressed this problem by identifying factors causing such artifacts and methods to reduce them.

Body Surface Potential Mapping: Contemporary Applications and Future Perspectives

Body surface potential mapping (BSPM) is a noninvasive modality to assess cardiac bioelectric activity with a rich history of practical applications for both research and clinical investigation. BSPM provides comprehensive acquisition of bioelectric signals across the entire thorax, allowing for more complex and extensive analysis than the standard electrocardiogram (ECG). Despite its advantages, BSPM is not a common clinical tool. BSPM does, however, serve as a valuable research tool and as an input for other modes of analysis such as electrocardiographic imaging and, more recently, machine learning and artificial intelligence. In this report, we examine contemporary uses of BSPM, and provide an assessment of its future prospects in both clinical and research environments. We assess the state of the art of BSPM implementations and explore modern applications of advanced modeling and statistical analysis of BSPM data. We predict that BSPM will continue to be a valuable research tool, and will find clinical utility at the intersection of computational modeling approaches and artificial intelligence.

Noninvasive detection of spatiotemporal activation-repolarization interactions that prime idiopathic ventricular fibrillation

[Figure: see text].

The Use of Electrocardiographic Imaging in Localising the Origin of Arrhythmias During Catheter Ablation of Ventricular Tachycardia

Non-invasive electrocardiographic imaging (ECGI) is a novel clinical tool for mapping ventricular arrhythmia. Using multiple body surface electrodes to collect unipolar electrograms and conventional medical imaging of the heart, an epicardial shell can be created to display calculated electrograms. This calculation is achieved by solving the inverse problem and allows activation times to be calculated from a single beat. The technology was initially pioneered in the US using an experimental torso-shaped tank. Accuracy from studies in humans has varied. Early data was promising, with more recent work suggesting only moderate accuracy when reproducing cardiac activation. Despite these limitations, the system has been successfully used in pioneering work with non-invasive cardiac radioablation to treat ventricular arrhythmia. This suggests that the resolution may be sufficient for treatment of large target areas. Although untested in a well conducted clinical study it is likely that it would not be accurate enough to guide more discreet radiofrequency ablation.

Ex vivo Validation of Noninvasive Epicardial and Endocardial Repolarization Mapping

Background: The detection and localization of electrophysiological substrates currently involve invasive cardiac mapping. Electrocardiographic imaging (ECGI) using the equivalent dipole layer (EDL) method allows the noninvasive estimation of endocardial and epicardial activation and repolarization times (AT and RT), but the RT validation is limited to in silico studies. We aimed to assess the temporal and spatial accuracy of the EDL method in reconstructing the RTs from the surface ECG under physiological circumstances and situations with artificially induced increased repolarization heterogeneity.

Methods: In four Langendorff-perfused pig hearts, we simultaneously recorded unipolar electrograms from plunge needles and pseudo-ECGs from a volume-conducting container equipped with 61 electrodes. The RTs were computed from the ECGs during atrial and ventricular pacing and compared with those measured from the local unipolar electrograms. Regional RT prolongation (cooling) or shortening (pinacidil) was achieved by selective perfusion of the left anterior descending artery (LAD) region.

Results: The differences between the computed and measured RTs were 19.0 ± 17.8 and 18.6 ± 13.7 ms for atrial and ventricular paced beats, respectively. The region of artificially delayed or shortened repolarization was correctly identified, with minimum/maximum RT roughly in the center of the region in three hearts. In one heart, the reconstructed region was shifted by ~2.5 cm. The total absolute difference between the measured and calculated RTs for all analyzed patterns in selectively perfused hearts (n = 5) was 39.6 ± 27.1 ms.

Conclusion: The noninvasive ECG repolarization imaging using the EDL method of atrial and ventricular paced beats allows adequate quantitative reconstruction of regions of altered repolarization.

Ultrafast four-dimensional imaging of cardiac mechanical wave propagation with sparse optoacoustic sensing

Propagation of electromechanical waves in excitable heart muscles follows complex spatiotemporal patterns holding the key to understanding life-threatening arrhythmias and other cardiac conditions. Accurate volumetric mapping of cardiac wave propagation is currently hampered by fast heart motion, particularly in small model organisms. Here we demonstrate that ultrafast four-dimensional imaging of cardiac mechanical wave propagation in entire beating murine heart can be accomplished by sparse optoacoustic sensing with high contrast, ∼115-µm spatial and submillisecond temporal resolution. We extract accurate dispersion and phase velocity maps of the cardiac waves and reveal vortex-like patterns associated with mechanical phase singularities that occur during arrhythmic events induced via burst ventricular electric stimulation. The newly introduced cardiac mapping approach is a bold step toward deciphering the complex mechanisms underlying cardiac arrhythmias and enabling precise therapeutic interventions.

Comparing Non-invasive Inverse Electrocardiography With Invasive Endocardial and Epicardial Electroanatomical Mapping During Sinus Rhythm

This study presents a novel non-invasive equivalent dipole layer (EDL) based inverse electrocardiography (iECG) technique which estimates both endocardial and epicardial ventricular activation sequences. We aimed to quantitatively compare our iECG approach with invasive electro-anatomical mapping (EAM) during sinus rhythm with the objective of enabling functional substrate imaging and sudden cardiac death risk stratification in patients with cardiomyopathy. Thirteen patients (77% males, 48 ± 20 years old) referred for endocardial and epicardial EAM underwent 67-electrode body surface potential mapping and CT imaging. The EDL-based iECG approach was improved by mimicking the effects of the His-Purkinje system on ventricular activation. EAM local activation timing (LAT) maps were compared with iECG-LAT maps using absolute differences and Pearson’s correlation coefficient, reported as mean ± standard deviation [95% confidence interval]. The correlation coefficient between iECG-LAT maps and EAM was 0.54 ± 0.19 [0.49–0.59] for epicardial activation, 0.50 ± 0.27 [0.41–0.58] for right ventricular endocardial activation and 0.44 ± 0.29 [0.32–0.56] for left ventricular endocardial activation. The absolute difference in timing between iECG maps and EAM was 17.4 ± 7.2 ms for epicardial maps, 19.5 ± 7.7 ms for right ventricular endocardial maps, 27.9 ± 8.7 ms for left ventricular endocardial maps. The absolute distance between right ventricular endocardial breakthrough sites was 30 ± 16 mm and 31 ± 17 mm for the left ventricle. The absolute distance for latest epicardial activation was median 12.8 [IQR: 2.9–29.3] mm. This first in-human quantitative comparison of iECG and invasive LAT-maps on both the endocardial and epicardial surface during sinus rhythm showed improved agreement, although with considerable absolute difference and moderate correlation coefficient. Non-invasive iECG requires further refinements to facilitate clinical implementation and risk stratification.

A Novel Data-Adaptive Regression Framework Based on Multivariate Adaptive Regression Splines for Electrocardiographic Imaging

OBJECTIVE

Noninvasive electrocardiographic imaging (ECGI) is a promising tool for revealing crucial cardiac electrical events with diagnostic potential. We propose a novel nonparametric regression framework based on multivariate adaptive regression splines (MARS) for ECGI.

METHODS

The inverse problem was solved by using the regression model trained with body surface potentials (BSP) and corresponding electrograms (EGM). Simulated data as well as experimental data from torso-tank experiments were used as to assess the performance of the proposed method. The robustness of the method to measurement noise and geometric errors were assessed in terms of electrogram reconstruction quality, activation time accuracy, and localization error metrics. The methods were compared with Tikhonov regularization and neural network (NN)-based methods. The resulting mapping functions between the BSPs and EGMs were also used to evaluate the most influential measurement leads.

RESULTS

MARS-based method outperformed Tikhonov regularization in terms of reconstruction accuracy and robustness to measurement noise. The effects of geometric errors were remedied to some extent by enriching the training set composition including model errors. The MARS-based method had a comparable performance with NN-based methods, which require the adjustment of many parameters.

CONCLUSION

MARS-based method successfully discovers the inverse mapping functions between the BSPs and EGMs yielding accurate reconstructions, and quantifies the contribution of each BSP lead.

SIGNIFICANCE

MARS-based method is adaptive, requires fewer parameter adjustments than NN-based methods, and robust to errors. Thus, it can be a feasible data-driven approach for accurately solving inverse imaging problems.

Cardiac Activation Maps Reconstruction: A Comparative Study Between Data-Driven and Physics-Based Methods

One of the essential diagnostic tools of cardiac arrhythmia is activation mapping. Noninvasive current mapping procedures include electrocardiographic imaging. It allows reconstructing heart surface potentials from measured body surface potentials. Then, activation maps are generated using the heart surface potentials. Recently, a study suggests to deploy artificial neural networks to estimate activation maps directly from body surface potential measurements. Here we carry out a comparative study between the data-driven approach DirectMap and noninvasive classic technique based on reconstructed heart surface potentials using both Finite element method combined with L1-norm regularization (FEM-L1) and the spatial adaptation of Time-delay neural networks (SATDNN-AT). In this work, we assess the performance of the three approaches using a synthetic single paced-rhythm dataset generated on the atria surface. The results show that data-driven approach DirectMap quantitatively outperforms the two other methods. In fact, we observe an absolute activation time error and a correlation coefficient, respectively, equal to 7.20 ms, 93.2% using DirectMap, 14.60 ms, 76.2% using FEM-L1 and 13.58 ms, 79.6% using SATDNN-AT. In addition, results show that data-driven approaches (DirectMap and SATDNN-AT) are strongly robust against additive gaussian noise compared to FEM-L1.

Fusion Pacing with Biventricular, Left Ventricular-only and Multipoint Pacing in Cardiac Resynchronisation Therapy: Latest Evidence and Strategies for Use

Despite advances in the field of cardiac resynchronisation therapy (CRT), response rates and durability of therapy remain relatively static. Optimising device timing intervals may be the most common modifiable factor influencing CRT efficacy after implantation. This review addresses the concept of fusion pacing as a method for improving patient outcomes with CRT. Fusion pacing describes the delivery of CRT pacing with a programming strategy to preserve intrinsic atrioventricular (AV) conduction and ventricular activation via the right bundle branch. Several methods have been assessed to achieve fusion pacing. QRS complex duration (QRSd) shortening with CRT is associated with improved clinical response. Dynamic algorithm-based optimisation targeting narrowest QRSd in patients with intact AV conduction has shown promise in people with heart failure with left bundle branch block. Individualised dynamic programming achieving fusion may achieve the greatest magnitude of electrical synchrony, measured by QRSd narrowing.

Reconstruction of three-dimensional biventricular activation based on the 12-lead electrocardiogram via patient-specific modelling

Aims

Non-invasive imaging of electrical activation requires high-density body surface potential mapping. The nine electrodes of the 12-lead electrocardiogram (ECG) are insufficient for a reliable reconstruction with standard inverse methods. Patient-specific modelling may offer an alternative route to physiologically constraint the reconstruction. The aim of the study was to assess the feasibility of reconstructing the fully 3D electrical activation map of the ventricles from the 12-lead ECG and cardiovascular magnetic resonance (CMR).

Methods and results

Ventricular activation was estimated by iteratively optimizing the parameters (conduction velocity and sites of earliest activation) of a patient-specific model to fit the simulated to the recorded ECG. Chest and cardiac anatomy of 11 patients (QRS duration 126–180 ms, documented scar in two) were segmented from CMR images. Scar presence was assessed by magnetic resonance (MR) contrast enhancement. Activation sequences were modelled with a physiologically based propagation model and ECGs with lead field theory. Validation was performed by comparing reconstructed activation maps with those acquired by invasive electroanatomical mapping of coronary sinus/veins (CS) and right ventricular (RV) and left ventricular (LV) endocardium. The QRS complex was correctly reproduced by the model (Pearson’s correlation r = 0.923). Reconstructions accurately located the earliest and latest activated LV regions (median barycentre distance 8.2 mm, IQR 8.8 mm). Correlation of simulated with recorded activation time was very good at LV endocardium (r = 0.83) and good at CS (r = 0.68) and RV endocardium (r = 0.58).

Conclusion

Non-invasive assessment of biventricular 3D activation using the 12-lead ECG and MR imaging is feasible. Potential applications include patient-specific modelling and pre-/per-procedural evaluation of ventricular activation.

ESC working group on cardiac cellular electrophysiology position paper: relevance, opportunities, and limitations of experimental models for cardiac electrophysiology research

Cardiac arrhythmias are a major cause of death and disability. A large number of experimental cell and animal models have been developed to study arrhythmogenic diseases. These models have provided important insights into the underlying arrhythmia mechanisms and translational options for their therapeutic management. This position paper from the ESC Working Group on Cardiac Cellular Electrophysiology provides an overview of (i) currently available in vitro, ex vivo, and in vivo electrophysiological research methodologies, (ii) the most commonly used experimental (cellular and animal) models for cardiac arrhythmias including relevant species differences, (iii) the use of human cardiac tissue, induced pluripotent stem cell (hiPSC)-derived and in silico models to study cardiac arrhythmias, and (iv) the availability, relevance, limitations, and opportunities of these cellular and animal models to recapitulate specific acquired and inherited arrhythmogenic diseases, including atrial fibrillation, heart failure, cardiomyopathy, myocarditis, sinus node, and conduction disorders and channelopathies. By promoting a better understanding of these models and their limitations, this position paper aims to improve the quality of basic research in cardiac electrophysiology, with the ultimate goal to facilitate the clinical translation and application of basic electrophysiological research findings on arrhythmia mechanisms and therapies.

Mapping Technologies for Catheter Ablation of Atrial Fibrillation Beyond Pulmonary Vein Isolation

Catheter ablation remains the most effective and relatively minimally invasive therapy for rhythm control in patients with AF. Ablation has consistently shown a reduction of arrhythmia-related symptoms and significant improvement in patients’ quality of life compared with medical treatment. The ablation strategy relies on a well-established anatomical approach of effective pulmonary vein isolation. Additional anatomical targets have been reported with the aim of increasing procedure success in complex substrates. However, larger ablated areas with uncertainty of targeting relevant regions for AF initiation or maintenance are not exempt from the potential risk of complications and pro-arrhythmia. Recent developments in mapping tools and computational methods for advanced signal processing during AF have reported novel strategies to identify atrial regions associated with AF maintenance. These novel tools – although mainly limited to research series – represent a significant step forward towards the understanding of complex patterns of propagation during AF and the potential achievement of patient-tailored AF ablation strategies for the near future.

Simultaneous Multi-Heartbeat ECGI Solution with a Time-Varying Forward Model: a Joint Inverse Formulation

Electrocardiographic imaging (ECGI) is an effective tool for noninvasive diagnosis of a range of cardiac dysfunctions. ECGI leverages a model of how cardiac bioelectric sources appear on the torso surface (the forward problem) and uses recorded body surface potential signals to reconstruct the bioelectric source (the inverse problem). Solutions to the inverse problem are sensitive to noise and variations in the body surface potential (BSP) recordings such as those caused by changes or errors in cardiac position. Techniques such as signal averaging seek to improve ECGI solutions by incorporating BSP signals from multiple heartbeats into an averaged BSP with a higher SNR to use when estimating the cardiac bioelectric source. However, signal averaging is limited when it comes to addressing sources of BSP variability such as beat to beat differences in the forward solution. We present a novel joint inverse formulation to solve for the cardiac source given multiple BSP recordings and known changes in the forward solution, here changes in the heart position. We report improved ECGI accuracy over signal averaging and averaged individual inverse solutions using this joint inverse formulation across multiple activation sequence types and regularization techniques with measured canine data and simulated heart motion. Our joint inverse formulation builds upon established techniques and consequently can easily be applied with many existing regularization techniques, source models, and forward problem formulations.

3-Dimensional ventricular electrical activation pattern assessed from a novel high-frequency electrocardiographic imaging technique: principles and clinical importance

The study introduces and validates a novel high-frequency (100–400 Hz bandwidth, 2 kHz sampling frequency) electrocardiographic imaging (HFECGI) technique that measures intramural ventricular electrical activation. Ex-vivo experiments and clinical measurements were employed. Ex-vivo, two pig hearts were suspended in a human-torso shaped tank using surface tank electrodes, epicardial electrode sock, and plunge electrodes. We compared conventional epicardial electrocardiographic imaging (ECGI) with intramural activation by HFECGI and verified with sock and plunge electrodes. Clinical importance of HFECGI measurements was performed on 14 patients with variable conduction abnormalities. From 3 × 4 needle and 108 sock electrodes, 256 torso or 184 body surface electrodes records, transmural activation times, sock epicardial activation times, ECGI-derived activation times, and high-frequency activation times were computed. The ex-vivo transmural measurements showed that HFECGI measures intramural electrical activation, and ECGI-HFECGI activation times differences indicate endo-to-epi or epi-to-endo conduction direction. HFECGI-derived volumetric dyssynchrony was significantly lower than epicardial ECGI dyssynchrony. HFECGI dyssynchrony was able to distinguish between intraventricular conduction disturbance and bundle branch block patients.

The electrocardiographic forward problem: A benchmark study

BACKGROUND

Electrocardiographic forward problems are crucial components for noninvasive electrocardiographic imaging (ECGI) that compute torso potentials from cardiac source measurements. Forward problems have few sources of error as they are physically well posed and supported by mature numerical and computational techniques. However, the residual errors reported from experimental validation studies between forward computed and measured torso signals remain surprisingly high.

OBJECTIVE

To test the hypothesis that incomplete cardiac source sampling, especially above the atrioventricular (AV) plane is a major contributor to forward solution errors.

METHODS

We used a modified Langendorff preparation suspended in a human-shaped electrolytic torso-tank and a novel pericardiac-cage recording array to thoroughly sample the cardiac potentials. With this carefully controlled experimental preparation, we minimized possible sources of error, including geometric error and torso inhomogeneities. We progressively removed recorded signals from above the atrioventricular plane to determine how the forward-computed torso-tank potentials were affected by incomplete source sampling.

RESULTS

We studied 240 beats total recorded from three different activation sequence types (sinus, and posterior and anterior left-ventricular free-wall pacing) in each of two experiments. With complete sampling by the cage electrodes, all correlation metrics between computed and measured torso-tank potentials were above 0.93 (maximum 0.99). The mean root-mean-squared error across all beat types was also low, less than or equal to 0.10 mV. A precipitous drop in forward solution accuracy was observed when we included only cage measurements below the AV plane.

CONCLUSION

First, our forward computed potentials using complete cardiac source measurements set a benchmark for similar studies. Second, this study validates the importance of complete cardiac source sampling above the AV plane to produce accurate forward computed torso potentials. Testing ECGI systems and techniques with these more complete and highly accurate datasets will improve inverse techniques and noninvasive detection of cardiac electrical abnormalities.

Electrocardiographic Imaging for Atrial Fibrillation: A Perspective From Computer Models and Animal Experiments to Clinical Value

Electrocardiographic imaging (ECGI) is a technique to reconstruct non-invasively the electrical activity on the heart surface from body-surface potential recordings and geometric information of the torso and the heart. ECGI has shown scientific and clinical value when used to characterize and treat both atrial and ventricular arrhythmias. Regarding atrial fibrillation (AF), the characterization of the electrical propagation and the underlying substrate favoring AF is inherently more challenging than for ventricular arrhythmias, due to the progressive and heterogeneous nature of the disease and its manifestation, the small volume and wall thickness of the atria, and the relatively large role of microstructural abnormalities in AF. At the same time, ECGI has the advantage over other mapping technologies of allowing a global characterization of atrial electrical activity at every atrial beat and non-invasively. However, since ECGI is time-consuming and costly and the use of electrical mapping to guide AF ablation is still not fully established, the clinical value of ECGI for AF is still under assessment. Nonetheless, AF is known to be the manifestation of a complex interaction between electrical and structural abnormalities and therefore, true electro-anatomical-structural imaging may elucidate important key factors of AF development, progression, and treatment. Therefore, it is paramount to identify which clinical questions could be successfully addressed by ECGI when it comes to AF characterization and treatment, and which questions may be beyond its technical limitations. In this manuscript we review the questions that researchers have tried to address on the use of ECGI for AF characterization and treatment guidance (for example, localization of AF triggers and sustaining mechanisms), and we discuss the technological requirements and validation. We address experimental and clinical results, limitations, and future challenges for fruitful application of ECGI for AF understanding and management. We pay attention to existing techniques and clinical application, to computer models and (animal or human) experiments, to challenges of methodological and clinical validation. The overall objective of the study is to provide a consensus on valuable directions that ECGI research may take to provide future improvements in AF characterization and treatment guidance.

Electrocardiographic Imaging of Repolarization Abnormalities

Background

Dispersion and gradients in repolarization have been associated with life‐threatening arrhythmias, but are difficult to quantify precisely from surface electrocardiography. The objective of this study was to evaluate electrocardiographic imaging (ECGI) to noninvasively detect repolarization‐based abnormalities.

Methods and Results

Ex vivo data were obtained from Langendorff‐perfused pig hearts (n=8) and a human donor heart. Unipolar electrograms were recorded simultaneously during sinus rhythm from an epicardial sock and the torso‐shaped tank within which the heart was suspended. Regional repolarization heterogeneities were introduced through perfusion of dofetilide and pinacidil into separate perfusion beds. In vivo data included torso and epicardial potentials recorded simultaneously in anesthetized, closed‐chest pigs (n=5), during sinus rhythm, and ventricular pacing. For both data sets, ECGI accurately reconstructed T‐wave electrogram morphologies when compared with those recorded by the sock (ex vivo: correlation coefficient, 0.85 [0.52–0.96], in vivo: correlation coefficient, 0.86 [0.52–0.96]) and repolarization time maps (ex‐vivo: correlation coefficient, 0.73 [0.63–0.83], in vivo: correlation coefficient, 0.76 [0.67–0.82]). ECGI‐reconstructed repolarization time distributions were strongly correlated to those measured by the sock (both data sets, R² ≥0.92). Although the position of the gradient was slightly shifted by 8.3 (0–13.9) mm, the mean, max, and SD between ECGI and recorded gradient values were highly correlated (R²=0.87, 0.75, and 0.86 respectively). There was no significant difference in ECGI accuracy between ex vivo and in vivo data.

Conclusions

ECGI reliably and accurately maps potentially critical repolarization abnormalities. This noninvasive approach allows imaging and quantifying individual parameters of abnormal repolarization‐based substrates in patients with arrhythmogenesis, to improve diagnosis and risk stratification.

Non-invasive Spatial Mapping of Frequencies in Atrial Fibrillation: Correlation With Contact Mapping

Introduction: Regional differences in activation rates may contribute to the electrical substrates that maintain atrial fibrillation (AF), and estimating them non-invasively may help guide ablation or select anti-arrhythmic medications. We tested whether non-invasive assessment of regional AF rate accurately represents intracardiac recordings.

Methods: In 47 patients with AF (27 persistent, age 63 ± 13 years) we performed 57-lead non-invasive Electrocardiographic Imaging (ECGI) in AF, simultaneously with 64-pole intracardiac signals of both atria. ECGI was reconstructed by Tikhonov regularization. We constructed personalized 3D AF rate distribution maps by Dominant Frequency (DF) analysis from intracardiac and non-invasive recordings.

Results: Raw intracardiac and non-invasive DF differed substantially, by 0.54 Hz [0.13 – 1.37] across bi-atrial regions (R² = 0.11). Filtering by high spectral organization reduced this difference to 0.10 Hz (cycle length difference of 1 – 11 ms) [0.03 – 0.42] for patient-level comparisons (R² = 0.62), and 0.19 Hz [0.03 – 0.59] and 0.20 Hz [0.04 – 0.61] for median and highest DF, respectively. Non-invasive and highest DF predicted acute ablation success (p = 0.04).

Conclusion: Non-invasive estimation of atrial activation rates is feasible and, when filtered by high spectral organization, provide a moderate estimate of intracardiac recording rates in AF. Non-invasive technology could be an effective tool to identify patients who may respond to AF ablation for personalized therapy.

The Impact of Torso Signal Processing on Noninvasive Electrocardiographic Imaging Reconstructions

Goal:

To evaluate state-of-the-art signal processing methods for epicardial potential-based noninvasive electrocardiographic imaging reconstructions of single-site pacing data.

Methods:

Experimental data were obtained from two torso-tank setups in which Langendorff-perfused hearts (n = 4) were suspended and potentials recorded simultaneously from torso and epicardial surfaces. 49 different signal processing methods were applied to torso potentials, grouped as i) high-frequency noise removal (HFR) methods ii) baseline drift removal (BDR) methods and iii) combined HFR+BDR. The inverse problem was solved and reconstructed electrograms and activation maps compared to those directly recorded.

Results:

HFR showed no difference compared to not filtering in terms of absolute differences in reconstructed electrogram amplitudes nor median correlation in QRS waveforms (p > 0.05). However, correlation and mean absolute error of activation times and pacing site localization were improved with all methods except a notch filter. HFR applied post-reconstruction produced no differences compared to pre-reconstruction. BDR and BDR+HFR significantly improved absolute and relative difference, and correlation in electrograms (p < 0.05). While BDR+HFR combined improved activation time and pacing site detection, BDR alone produced significantly lower correlation and higher localization errors (p < 0.05).

Conclusion:

BDR improves reconstructed electrogram morphologies and amplitudes due to a reduction in lambda value selected for the inverse problem. The simplest method (resetting the isoelectric point) is sufficient to see these improvements. HFR does not impact electrogram accuracy, but does impact post-processing to extract features such as activation times. Removal of line noise is insufficient to see these changes. HFR should be applied post-reconstruction to ensure over-filtering does not occur.

IRN-MLSQR: An improved iterative reweight norm approach to the inverse problem of electrocardiography incorporating factorization-free preconditioned LSQR

The inverse problem of electrocardiography (ECG) of computing epicardial potentials from body surface potentials, is an ill-posed problem and needs to be solved by regularization techniques. The L2-norm regularization can cause considerable smoothing of the solution, while the L1-norm scheme promotes a solution with sharp boundaries/gradients between piecewise smooth regions, so L1-norm is widely used in the ECG inverse problem. However, large amount of computation and long computation time are needed in the L1-norm scheme. In this paper, by combining iterative reweight norm (IRN) with a factorization-free preconditioned LSQR algorithm (MLSQR), a new IRN-MLSQR method was proposed to accelerate the convergence speed of the L1-norm scheme. We validated the IRN-MLSQR method using experimental data from isolated canine hearts and clinical procedures in the electrophysiology laboratory. The results showed that the IRN-MLSQR method can significantly reduce the number of iterations and operation time while ensuring the calculation accuracy. The number of iterations of the IRN-MLSQR method is about 60%-70% that of the conventional IRN method, and at the same time, the accuracy of the solution is almost the same as that of the conventional IRN method. The proposed IRN-MLSQR method may be used as a new approach to the inverse problem of ECG.

Electrocardiographic imaging for cardiac arrhythmias and resynchronization therapy

Use of the 12-lead electrocardiogram (ECG) is fundamental for the assessment of heart disease, including arrhythmias, but cannot always reveal the underlying mechanism or the location of the arrhythmia origin. Electrocardiographic imaging (ECGi) is a non-invasive multi-lead ECG-type imaging tool that enhances conventional 12-lead ECG. Although it is an established technology, its continuous development has been shown to assist in arrhythmic activation mapping and provide insights into the mechanism of cardiac resynchronization therapy (CRT). This review addresses the validity, reliability, and overall feasibility of ECGi for use in a diverse range of arrhythmias. A systematic search limited to full-text human studies published in peer-reviewed journals was performed through Medline via PubMed, using various combinations of three key concepts: ECGi, arrhythmia, and CRT. A total of 456 studies were screened through titles and abstracts. Ultimately, 42 studies were included for literature review. Evidence to date suggests that ECGi can be used to provide diagnostic insights regarding the mechanistic basis of arrhythmias and the location of arrhythmia origin. Furthermore, ECGi can yield valuable information to guide therapeutic decision-making, including during CRT. Several studies have used ECGi as a diagnostic tool for atrial and ventricular arrhythmias. More recently, studies have tested the value of this technique in predicting outcomes of CRT. As a non-invasive method for assessing cardiovascular disease, particularly arrhythmias, ECGi represents a significant advancement over standard procedures in contemporary cardiology. Its full potential has yet to be fully explored.

In silico validation of electrocardiographic imaging to reconstruct the endocardial and epicardial repolarization pattern using the equivalent dipole layer source model

The solution of the inverse problem of electrocardiology allows the reconstruction of the spatial distribution of the electrical activity of the heart from the body surface electrocardiogram (electrocardiographic imaging, ECGI). ECGI using the equivalent dipole layer (EDL) model has shown to be accurate for cardiac activation times. However, validation of this method to determine repolarization times is lacking. In the present study, we determined the accuracy of the EDL model in reconstructing cardiac repolarization times, and assessed the robustness of the method under less ideal conditions (addition of noise and errors in tissue conductivity). A monodomain model was used to determine the transmembrane potentials in three different excitation-repolarization patterns (sinus beat and ventricular ectopic beats) as the gold standard. These were used to calculate the body surface ECGs using a finite element model. The resulting body surface electrograms (ECGs) were used as input for the EDL-based inverse reconstruction of repolarization times. The reconstructed repolarization times correlated well (COR > 0.85) with the gold standard, with almost no decrease in correlation after adding errors in tissue conductivity of the model or noise to the body surface ECG. Therefore, ECGI using the EDL model allows adequate reconstruction of cardiac repolarization times. Graphical abstract Validation of electrocardiographic imaging for repolarization using forward calculated body surface ECGs from simulated activation-repolarization sequences.