How Accurate Is Inverse Electrocardiographic Mapping? A Systematic In Vivo Evaluation

Systematic review of computational techniques, dataset utilization, and feature extraction in electrocardiographic imaging

This study aimed to analyze computational techniques in ECG imaging (ECGI) reconstruction, focusing on dataset identification, problem-solving, and feature extraction. We employed a PRISMA approach to review studies from Scopus and Web of Science, applying Cochrane principles to assess risk of bias. The selection was limited to English peer-reviewed papers published from 2010 to 2023, excluding studies that lacked computational technique descriptions. From 99 reviewed papers, trends show a preference for traditional methods like the boundary element and Tikhonov methods, alongside a rising use of advanced technologies including hybrid techniques and deep learning. These advancements have enhanced cardiac diagnosis and treatment precision. Our findings underscore the need for robust data utilization and innovative computational integration in ECGI, highlighting promising areas for future research and advances. This shift toward tailored cardiac care suggests significant progress in diagnostic and treatment methods.

Solving the Inverse Problem of Electrocardiography for Cardiac Digital Twins: A Survey

Cardiac digital twins (CDTs) are personalized virtual representations used to understand complex cardiac mechanisms. A critical component of CDT development is solving the ECG inverse problem, which enables the reconstruction of cardiac sources and the estimation of patient-specific electrophysiology (EP) parameters from surface ECG data. Despite challenges from complex cardiac anatomy, noisy ECG data, and the ill-posed nature of the inverse problem, recent advances in computational methods have greatly improved the accuracy and efficiency of ECG inverse inference, strengthening the fidelity of CDTs. This paper aims to provide a comprehensive review of the methods for solving ECG inverse problems, their validation strategies, their clinical applications, and their future perspectives. For the methodologies, we broadly classify state-of-the-art approaches into two categories: deterministic and probabilistic methods, including both conventional and deep learning-based techniques. Integrating physics laws with deep learning models holds promise, but challenges such as capturing dynamic electrophysiology accurately, accessing accurate domain knowledge, and quantifying prediction uncertainty persist. Integrating models into clinical workflows while ensuring interpretability and usability for healthcare professionals is essential. Overcoming these challenges will drive further research in CDTs.

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.

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.

A reduced complexity ECG imaging model for regularized inversion optimization

The resolution of the inverse problem of electrocardiography represents a major interest in the diagnosis and catheter-based therapy of cardiac arrhythmia. In this context, the ability to simulate several cardiac electrical behaviors was crucial for evaluating and comparing the performance of inversion methods. For this application, existing models are either too complex or do not produce realistic cardiac patterns. In this work, a low-resolution heart-torso model generating realistic whole heart cardiac mappings and electrocardiograms in healthy and pathological cases is designed. This model was built upon a simplified heart-torso geometry and implements the monodomain formalism by using the finite element method. In addition, a model reduction step through a sensitivity analysis was proposed where parameters were identified using an evolutionary optimization approach. Finally, the study illustrates the usefulness of the proposed model by comparing the performance of different variants of Tikhonov-based inversion methods for the determination of the regularization parameter in healthy, ischemic and ventricular tachycardia scenarios. First, results of the sensitivity analysis show that among 58 parameters only 25 are influent. Note also that the level of influence of the parameters depends on the heart region. Besides, the synthesized electrocardiograms globally present the same characteristic shape compared to the reference once with a correlation value that reaches 88%. Regarding inverse problem, results highlight that only Robust Generalized Cross Validation and Discrepancy Principle provide best performance, with a quasi-perfect success rate for both, and a respective relative error, between the generated electrocardiograms to the reference one, of 0.75 and 0.62.

Non-invasive localization of the ventricular excitation origin without patient-specific geometries using deep learning

Cardiovascular diseases account for 17 million deaths per year worldwide. Of these, 25% are categorized as sudden cardiac death, which can be related to ventricular tachycardia (VT). This type of arrhythmia can be caused by focal activation sources outside the sinus node. Catheter ablation of these foci is a curative treatment in order to inactivate the abnormal triggering activity. However, the localization procedure is usually time-consuming and requires an invasive procedure in the catheter lab. To facilitate and expedite the treatment, we present two novel localization support techniques based on convolutional neural networks (CNNs) that address these clinical needs. In contrast to existing methods, our approaches were designed to be independent of the patient-specific geometry and directly applicable to surface ECG signals, while also delivering a binary transmural position. Moreover, one of the method's outputs can be interpreted as several ranked solutions. The CNNs were trained on a dataset containing only simulated data and evaluated both on simulated test data and clinical data. On a novel large and open simulated dataset, the median test error was below 3 mm. The median localization error on the unseen clinical data ranged from 32 mm to 41 mm without optimizing the pre-processing and CNN to the clinical data. Interpreting the output of one of the approaches as ranked solutions, the best median error of the top-3 solutions decreased to 20 mm on the clinical data. The transmural position was correctly detected in up to 82% of all clinical cases. These results demonstrate a proof of principle to utilize CNNs to localize the activation source without the intrinsic need for patient-specific geometrical information. Furthermore, providing multiple solutions can assist physicians in identifying the true activation source amongst more than one possible location. With further optimization to clinical data, these methods have high potential to accelerate clinical interventions, replace certain steps within these procedures and consequently reduce procedural risk and improve VT patient outcomes.

Localization of the ventricular pacing site from BSPM and standard 12-lead ECG: a comparison study

Inverse ECG imaging methods typically require 32-250 leads to create body surface potential maps (BSPM), limiting their routine clinical use. This study evaluated the accuracy of PaceView inverse ECG method to localize the left or right ventricular (LV and RV, respectively) pacing leads using either a 99-lead BSPM or the 12-lead ECG. A 99-lead BSPM was recorded in patients with cardiac resynchronization therapy (CRT) during sinus rhythm and sequential LV/RV pacing. The non-contrast CT was performed to localize precisely both ECG electrodes and CRT leads. From a BSPM, nine signals were selected to obtain the 12-lead ECG. Both BSPM and 12-lead ECG were used to localize the RV and LV lead, and the localization error was calculated. Consecutive patients with dilated cardiomyopathy, previously implanted with a CRT device, were enrolled (n = 19). The localization error for the RV/LV lead was 9.0 [IQR 4.8-13.6] / 7.7 [IQR 0.0-10.3] mm using the 12-lead ECG and 9.1 [IQR 5.4-15.7] / 9.8 [IQR 8.6-13.1] mm for the BSPM. Thus, the noninvasive lead localization using the 12-lead ECG was accurate enough and comparable to 99-lead BSPM, potentially increasing the capability of 12-lead ECG for the optimization of the LV/RV pacing sites during CRT implant or for the most favorable programming.

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.

A novel framework for the removal of pacing artifacts from bio-electrical recordings

BACKGROUND

Electroceuticals provide clinical solutions for a range of disorders including Parkinson's disease, cardiac arrythmias and are emerging as a potential treatment option for gastrointestinal disorders. However, pre-clinical investigations are challenged by the large stimulation artifacts registered in bio-electrical recordings.

METHOD

A generalized framework capable of isolating and suppressing stimulation artifacts with minimal intervention was developed. Stimulation artifacts with different pulse-parameters in synthetic and experimental cardiac and gastrointestinal signals were detected using a Hampel filter and reconstructed using 3 methods: i) autoregression, ii) weighted mean, and iii) linear interpolation.

RESULTS

Synthetic stimulation artifacts with amplitudes of 2 mV and 4 mV and pulse-widths of 50 ms, 100 ms, and 200 ms were successfully isolated and the artifact window size remained uninfluenced by the pulse-amplitude, but was influenced by pulse-width (e.g., the autoregression method resulted in an identical Root Mean Square Error (RMSE) of 1.64 mV for artifacts with 200 ms pulse-width and both 2 mV and 4 mV amplitudes). The performance of autoregression (RMSE = 1.45 ± 0.16 mV) and linear interpolation (RMSE = 1.22 ± 0.14 mV) methods were comparable and better than weighted mean (RMSE = 5.54 ± 0.56 mV) for synthetic data. However, for experimental recordings, artifact removal by autoregression was superior to both linear interpolation and weighted mean approaches in gastric, small intestinal and cardiac recordings.

CONCLUSIONS

A novel signal processing framework enabled efficient analysis of bio-electrical recordings with stimulation artifacts. This will allow the bio-electrical events induced by stimulation protocols to be efficiently and systematically evaluated, resulting in improved stimulation therapies.

Tipping the scales of understanding: An engineering approach to design and implement whole-body cardiac electrophysiology experimental models

The study of cardiac electrophysiology is built on experimental models that span all scales, from ion channels to whole-body preparations. Novel discoveries made at each scale have contributed to our fundamental understanding of human cardiac electrophysiology, which informs clinicians as they detect, diagnose, and treat complex cardiac pathologies. This expert review describes an engineering approach to developing experimental models that is applicable across scales. The review also outlines how we applied the approach to create a set of multiscale whole-body experimental models of cardiac electrophysiology, models that are driving new insights into the response of the myocardium to acute ischemia. Specifically, we propose that researchers must address three critical requirements to develop an effective experimental model: 1) how the experimental model replicates and maintains human physiological conditions, 2) how the interventions possible with the experimental model capture human pathophysiology, and 3) what signals need to be measured, at which levels of resolution and fidelity, and what are the resulting requirements of the measurement system and the access to the organs of interest. We will discuss these requirements in the context of two examples of whole-body experimental models, a closed chest in situ model of cardiac ischemia and an isolated-heart, torso-tank preparation, both of which we have developed over decades and used to gather valuable insights from hundreds of experiments.

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.

A generalized framework for pacing artifact removal using a Hampel filter

Cardiac pacing is a clinical therapy widely used for treating irregular heart rhythms. Equivalent techniques for the treatment of gastric functional motility disorders hold great potential. Accurate analysis of pacing studies is often hindered by the stimulus artifacts which are superimposed on the recorded signals. This paper presents a semi-automated artifact detection method using a Hampel filter accompanied by 2 separate artifact reconstruction methods: (i) an auto-regressive model, and (ii) weighted mean to estimate the underlying signal. The developed framework was validated on synthetic and experimental signals containing large periodic pacing artifacts alongside evoked bioelectrical events. The performance of the proposed algorithms was quantified for gastric and cardiac pacing data collected in vivo. A lower mean RMS difference was achieved by the artifact segment reconstructed using the auto-regression ([Formula: see text]), method compared to the weighted mean ([Formula: see text]) method. Therefore, a more accurate artifact reconstruction was provided by the auto-regression approach. Clinical Relevance- The ability to efficiently and accurately isolate evoked bioelectrical events by eliminating large artifacts is a critical advancement for the analysis of paced recordings. The developed framework allows more efficient analysis of preclinical pacing data and thereby contributes to the advancement of pacing as a clinical therapy.

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.

Modeling the His-Purkinje Effect in Non-invasive Estimation of Endocardial and Epicardial Ventricular Activation

Inverse electrocardiography (iECG) estimates epi- and endocardial electrical activity from body surface potentials maps (BSPM). In individuals at risk for cardiomyopathy, non-invasive estimation of normal ventricular activation may provide valuable information to aid risk stratification to prevent sudden cardiac death. However, multiple simultaneous activation wavefronts initiated by the His-Purkinje system, severely complicate iECG. To improve the estimation of normal ventricular activation, the iECG method should accurately mimic the effect of the His-Purkinje system, which is not taken into account in the previously published multi-focal iECG. Therefore, we introduce the novel multi-wave iECG method and report on its performance. Multi-wave iECG and multi-focal iECG were tested in four patients undergoing invasive electro-anatomical mapping during normal ventricular activation. In each subject, 67-electrode BSPM were recorded and used as input for both iECG methods. The iECG and invasive local activation timing (LAT) maps were compared. Median epicardial inter-map correlation coefficient (CC) between invasive LAT maps and estimated multi-wave iECG versus multi-focal iECG was 0.61 versus 0.31. Endocardial inter-map CC was 0.54 respectively 0.22. Modeling the His-Purkinje system resulted in a physiologically realistic and robust non-invasive estimation of normal ventricular activation, which might enable the early detection of cardiac disease during normal sinus rhythm.

Critical appraisal of technologies to assess electrical activity during atrial fibrillation: a position paper from the European Heart Rhythm Association and European Society of Cardiology Working Group on eCardiology in collaboration with the Heart Rhythm Society, Asia Pacific Heart Rhythm Society, Latin American Heart Rhythm Society and Computing in Cardiology

We aim to provide a critical appraisal of basic concepts underlying signal recording and processing technologies applied for (i) atrial fibrillation (AF) mapping to unravel AF mechanisms and/or identifying target sites for AF therapy and (ii) AF detection, to optimize usage of technologies, stimulate research aimed at closing knowledge gaps, and developing ideal AF recording and processing technologies. Recording and processing techniques for assessment of electrical activity during AF essential for diagnosis and guiding ablative therapy including body surface electrocardiograms (ECG) and endo- or epicardial electrograms (EGM) are evaluated. Discussion of (i) differences in uni-, bi-, and multi-polar (omnipolar/Laplacian) recording modes, (ii) impact of recording technologies on EGM morphology, (iii) global or local mapping using various types of EGM involving signal processing techniques including isochronal-, voltage- fractionation-, dipole density-, and rotor mapping, enabling derivation of parameters like atrial rate, entropy, conduction velocity/direction, (iv) value of epicardial and optical mapping, (v) AF detection by cardiac implantable electronic devices containing various detection algorithms applicable to stored EGMs, (vi) contribution of machine learning (ML) to further improvement of signals processing technologies. Recording and processing of EGM (or ECG) are the cornerstones of (body surface) mapping of AF. Currently available AF recording and processing technologies are mainly restricted to specific applications or have technological limitations. Improvements in AF mapping by obtaining highest fidelity source signals (e.g. catheter-electrode combinations) for signal processing (e.g. filtering, digitization, and noise elimination) is of utmost importance. Novel acquisition instruments (multi-polar catheters combined with improved physical modelling and ML techniques) will enable enhanced and automated interpretation of EGM recordings in the near future.

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.

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.

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.

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.

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.

Strategies to Refine Gastric Stimulation and Pacing Protocols: Experimental and Modeling Approaches

Gastric pacing and stimulation strategies were first proposed in the 1960s to treat motility disorders. However, there has been relatively limited clinical translation of these techniques. Experimental investigations have been critical in advancing our understanding of the control mechanisms that innervate gut function. In this review, we will discuss the use of pacing to modulate the rhythmic slow wave conduction patterns generated by interstitial cells of Cajal in the gastric musculature. In addition, the use of gastric high-frequency stimulation methods that target nerves in the stomach to either inhibit or enhance stomach function will be discussed. Pacing and stimulation protocols to modulate gastric activity, effective parameters and limitations in the existing studies are summarized. Mathematical models are useful to understand complex and dynamic systems. A review of existing mathematical models and techniques that aim to help refine pacing and stimulation protocols are provided. Finally, some future directions and challenges that should be investigated are discussed.

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.

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.

A case report of successful elimination of recurrent ventricular tachycardia by repeated stereotactic radiotherapy: the importance of accurate target volume delineation

Background

Stereotactic body radiotherapy (SBRT) has emerged recently as a novel therapeutic alternative for patients with ventricular tachycardias (VTs) resistant to convetional treatment. Nevertheless, many aspects related to SBRT are currently unknown.

Case summary

A 66-year-old man with ischaemic heart disease, a history of coronary artery bypass graft surgery and left ventricular dysfunction was referred for recurrent symptomatic episodes of slow VT (108 b.p.m.). The arrhythmia was resistant to antiarrhythmic drug therapy with amiodarone and repeated catheter ablation. The patient was scheduled to SBRT, however, the first session failed to suppress VT recurrences. After 20 months, the patient underwent re-do ablation procedure that revealed a newly developed scar with its core adjacent to the presumed critical part of the VT substrate. Catheter ablation again failed to eliminate VT and the second session of SBRT was scheduled. To improve targeting of the VT substrate for SBRT, we applied our recently developed original method for integration of data from the electroanatomical mapping system with computer tomography images. The second session of SBRT with precise targeting using the novel strategy led within 3 months to the successful elimination of VT.

Discussion

This case report describes a patient in whom the recurrent VT was abolished only by properly targeted SBRT. Above all, the case highlights the importance of precise identification and targeting for SBRT. Our case also documents in vivo, by electroanatomical voltage mapping, the development of SBRT-related myocardial lesion. This represents an important mechanistic proof of the concept of SBRT.

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.

Spatial-Temporal Signals and Clinical Indices in Electrocardiographic Imaging (I): Preprocessing and Bipolar Potentials

During the last years, Electrocardiographic Imaging (ECGI) has emerged as a powerful and promising clinical tool to support cardiologists. Starting from a plurality of potential measurements on the torso, ECGI yields a noninvasive estimation of their causing potentials on the epicardium. This unprecedented amount of measured cardiac signals needs to be conditioned and adapted to current knowledge and methods in cardiac electrophysiology in order to maximize its support to the clinical practice. In this setting, many cardiac indices are defined in terms of the so-called bipolar electrograms, which correspond with differential potentials between two spatially close potential measurements. Our aim was to contribute to the usefulness of ECGI recordings in the current knowledge and methods of cardiac electrophysiology. For this purpose, we first analyzed the basic stages of conventional cardiac signal processing and scrutinized the implications of the spatial-temporal nature of signals in ECGI scenarios. Specifically, the stages of baseline wander removal, low-pass filtering, and beat segmentation and synchronization were considered. We also aimed to establish a mathematical operator to provide suitable bipolar electrograms from the ECGI-estimated epicardium potentials. Results were obtained on data from an infarction patient and from a healthy subject. First, the low-frequency and high-frequency noises are shown to be non-independently distributed in the ECGI-estimated recordings due to their spatial dimension. Second, bipolar electrograms are better estimated when using the criterion of the maximum-amplitude difference between spatial neighbors, but also a temporal delay in discrete time of about 40 samples has to be included to obtain the usual morphology in clinical bipolar electrograms from catheters. We conclude that spatial-temporal digital signal processing and bipolar electrograms can pave the way towards the usefulness of ECGI recordings in the cardiological clinical practice. The companion paper is devoted to analyzing clinical indices obtained from ECGI epicardial electrograms measuring waveform variability and repolarization tissue properties.

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.

Granger Causality-Based Analysis for Classification of Fibrillation Mechanisms and Localization of Rotational Drivers

BACKGROUND

The mechanisms sustaining myocardial fibrillation remain disputed, partly due to a lack of mapping tools that can accurately identify the mechanism with low spatial resolution clinical recordings. Granger causality (GC) analysis, an econometric tool for quantifying causal relationships between complex time-series, was developed as a novel fibrillation mapping tool and adapted to low spatial resolution sequentially acquired data.

METHODS

Ventricular fibrillation (VF) optical mapping was performed in Langendorff-perfused Sprague-Dawley rat hearts (n=18), where novel algorithms were developed using GC-based analysis to (1) quantify causal dependence of neighboring signals and plot GC vectors, (2) quantify global organization with the causality pairing index, a measure of neighboring causal signal pairs, and (3) localize rotational drivers (RDs) by quantifying the circular interdependence of neighboring signals with the circular interdependence value. GC-based mapping tools were optimized for low spatial resolution from downsampled optical mapping data, validated against high-resolution phase analysis and further tested in previous VF optical mapping recordings of coronary perfused donor heart left ventricular wedge preparations (n=12), and adapted for sequentially acquired intracardiac electrograms during human persistent atrial fibrillation mapping (n=16).

RESULTS

Global VF organization quantified by causality pairing index showed a negative correlation at progressively lower resolutions (50% resolution: P=0.006, R2=0.38, 12.5% resolution, P=0.004, R2=0.41) with a phase analysis derived measure of disorganization, locations occupied by phase singularities. In organized VF with high causality pairing index values, GC vector mapping characterized dominant propagating patterns and localized stable RDs, with the circular interdependence value showing a significant difference in driver versus nondriver regions (0.91±0.05 versus 0.35±0.06, P=0.0002). These findings were further confirmed in human VF. In persistent atrial fibrillation, a positive correlation was found between the causality pairing index and presence of stable RDs (P=0.0005,R2=0.56). Fifty percent of patients had RDs, with a low incidence of 0.9±0.3 RDs per patient.

CONCLUSIONS

GC-based fibrillation analysis can measure global fibrillation organization, characterize dominant propagating patterns, and map RDs using low spatial resolution sequentially acquired data.

Evaluation of ECG Imaging to Map Hemodynamically Stable and Unstable Ventricular Arrhythmias

BACKGROUND

ECG imaging (ECGI) has been used to guide treatment of ventricular ectopy and arrhythmias. However, the accuracy of ECGI in localizing the origin of arrhythmias during catheter ablation of ventricular tachycardia (VT) in structurally abnormal hearts remains to be fully validated.

METHODS

During catheter ablation of VT, simultaneous mapping was performed using electroanatomical mapping (CARTO, Biosense-Webster) and ECGI (CardioInsight, Medtronic) in 18 patients. Sites of entrainment, pace-mapping, and termination during ablation were used to define the VT site of origin (SoO). Distance between SoO and the site of earliest activation on ECGI were measured using co-registered geometries from both systems. The accuracy of ECGI versus a 12-lead surface ECG algorithm was compared.

RESULTS

A total of 29 VTs were available for comparison. Distance between SoO and sites of earliest activation in ECGI was 22.6, 13.9 to 36.2 mm (median, first to third quartile). ECGI mapped VT sites of origin onto the correct AHA segment with higher accuracy than a validated 12-lead ECG algorithm (83.3% versus 38.9%; P=0.015).

CONCLUSIONS

This simultaneous assessment demonstrates that CardioInsight localizes VT circuits with sufficient accuracy to provide a region of interest for targeting mapping for ablation. Resolution is not sufficient to guide discrete radiofrequency lesion delivery via catheter ablation without concomitant use of an electroanatomical mapping system but may be sufficient for segmental ablation with radiotherapy.

Validation of Dipole Density Mapping During Atrial Fibrillation and Sinus Rhythm in Human Left Atrium

OBJECTIVES

This study sought to validate the accuracy of noncontact electrograms against contact electrograms in the left atrium during sinus rhythm (SR) and atrial fibrillation (AF).

BACKGROUND

Noncontact mapping offers the opportunity to assess global cardiac activation in the chamber of interest. A novel noncontact mapping system, which records intracardiac voltage to derive cellular charge sources (dipole density), allows real-time mapping of AF to guide ablation.

METHODS

Noncontact and contact unipolar electrogram pairs were recorded simultaneously from multiple locations. Morphology correlation and timing difference of reconstructed electrograms obtained from a noncontact catheter were compared with those from contact electrograms obtained from a contact catheter at the same endocardial locations.

RESULTS

A total of 796 electrogram pairs in SR and 969 electrogram pairs in AF were compared from 20 patients with persistent AF. The median morphology correlation and timing difference (ms) in SR was 0.85 (interquartile range [IQR]: 0.71 to 0.94) and 6.4 ms (IQR: 2.6 to 17.1 ms); in AF was 0.79 (IQR: 0.69 to 0.88) and 14.4 ms (IQR: 6.7 to 26.2 ms), respectively. The correlation was stronger and the timing difference was less when the radial distance (r) from the noncontact catheter center to the endocardium was ≤ 40 versus > 40 mm; 0.87 (IQR: 0.72 to 0.94) versus 0.73 (IQR: 0.56 to 0.88) and 5.7 ms (IQR: 2.6 to 15.4 ms) versus 15.1 ms (IQR: 4.1 to 27.7 ms); p < 0.01 when in SR; 0.81 (IQR: 0.69 to 0.89) versus 0.67 (IQR: 0.45 to 0.82) and 12.3 ms (IQR: 5.9 to 21.8 ms) versus 28.3 ms (IQR: 16.2 to 36.0 ms); p < 0.01 when in AF.

CONCLUSIONS

This novel noncontact dipole density mapping system provides comparable reconstructed atrial electrogram measurements in SR or AF in human left atrium when the anatomical site of interest is ≤40 mm from the mapping catheter.

Understanding unipolar electrograms and global activation from noninvasive mapping for diagnosing arrhythmias

In the past several decades Electrocardiographic Imaging has been developed, validated, and recently commercialized for clinical use (Medtronic CardioInsight™). The technology noninvasively maps cardiac electrical potentials by combining multi-channel body surface electrocardiograms and patient specific heart-torso geometry. Unlike contact catheter mapping where bipolar electrograms are commonly used, noninvasive mapping is unipolar in nature. While bipolar electrograms reflect local activation time, noninvasive mapping reflects both local and global information. Characterizing arrhythmias using noninvasive mapping requires a different approach than what is used with bipolar mapping during electrophysiological studies. With multiple example cases, this article aims to demonstrate how in-depth analysis of noninvasive electrogram features and global activation patterns can reveal clinically important characteristics of arrhythmias, which may not be evident on the automatically generated maps. These characteristics include transmural location of the ectopic focus, differential locations within structural proximity, engagement of bundle or specific tissue, and substrate properties. The analysis approach involves applying unipolar potential fundamentals and relating electrogram features, global activation with underlying electrophysiology and cardiac anatomy to explain specific mapping phenomenon. If validated, easy-to-use diagnostic criteria or algorithms can be built and automated as product features in the commercial system to facilitate future noninvasive mapping in patients.

Spatial Accuracy of a Clinically Established Noninvasive Electrocardiographic Imaging System for the Detection of Focal Activation in an Intact Porcine Model

Background:

Noninvasive electrocardiographic imaging (ECGi) is used clinically to map arrhythmias before ablation. Despite its clinical use, validation data regarding the accuracy of the system for the identification of arrhythmia foci is limited.

Methods:

Nine pigs underwent closed-chest placement of endocardial fiducial markers, computed tomography, and pacing in all cardiac chambers with ECGi acquisition. Pacing location was reconstructed from biplane fluoroscopy and registered to the computed tomography using the fiducials. A blinded investigator predicted the pacing location from the ECGi data, and the distance to the true pacing catheter tip location was calculated.

Results:

A total of 109 endocardial and 9 epicardial locations were paced in 9 pigs. ECGi predicted the correct chamber of origin in 85% of atrial and 92% of ventricular sites. Lateral locations were predicted in the correct chamber more often than septal locations (97% versus 79%, P =0.01). Absolute distances in space between the true and predicted pacing locations were 20.7 (13.8–25.6) mm (median and [first–third] quartile). Distances were not significantly different across cardiac chambers.

Conclusions:

The ECGi system is able to correctly identify the chamber of origin for focal activation in the vast majority of cases. Determination of the true site of origin is possible with sufficient accuracy with consideration of these error estimates.

Noninvasive Mapping and Electrocardiographic Imaging in Atrial and Ventricular Arrhythmias (CardioInsight)

Electrocardiographic imaging is a mapping technique aiming to noninvasively characterize cardiac electrical activity using signals collected from the torso to reconstruct epicardial potentials. Its efficacy has been demonstrated clinically, from mapping premature ventricular complexes and accessory pathways to of complex arrhythmias. Electrocardiographic imaging uses a standardized workflow. Signals should be checked manually to avoid automatic processing errors. Reentry is confirmed in the presence of local activation covering the arrhythmia cycle length. Focal breakthroughs demonstrate a QS pattern associated with centrifugal activation. Electrocardiographic imaging offers a unique opportunity to better understand the mechanism of cardiac arrhythmias and guide ablation.

Synchronization-based reconstruction of electromechanical wave dynamics in elastic excitable media

The heart is an elastic excitable medium, in which mechanical contraction is triggered by nonlinear waves of electrical excitation, which diffuse rapidly through the heart tissue and subsequently activate the cardiac muscle cells to contract. These highly dynamic excitation wave phenomena have yet to be fully observed within the depths of the heart muscle, as imaging technology is unable to penetrate the tissue and provide panoramic, three-dimensional visualizations necessary for adequate study. As a result, the electrophysiological mechanisms that are associated with the onset and progression of severe heart rhythm disorders such as atrial or ventricular fibrillation remain insufficiently understood. Here, we present a novel synchronization-based data assimilation approach with which it is possible to reconstruct excitation wave dynamics within the volume of elastic excitable media by observing spatiotemporal deformation patterns, which occur in response to excitation. The mechanical data are assimilated in a numerical replication of the measured elastic excitable system, and within this replication, the data drive the intrinsic excitable dynamics, which then coevolve and correspond to a reconstruction of the original dynamics. We provide a numerical proof-of-principle and demonstrate the performance of the approach by recovering even complicated three-dimensional scroll wave patterns, including vortex filaments of electrical excitation from within a deformable bulk tissue with fiber anisotropy. In the future, the reconstruction approach could be combined with high-speed imaging of the heart's mechanical contractions to estimate its electrophysiological activity for diagnostic purposes.

Non-invasive electrocardiographic imaging in patients with idiopathic premature ventricular contractions from the right ventricular outflow tract: New insights into arrhythmia substrate

AIMS

The aim of this study was to use non-invasive electrocardiographic imaging (ECGI) to study the electrophysiological properties of right ventricular outflow tract (RVOT) in patients with frequent premature ventricular contractions (PVCs) from the RVOT and in controls.

METHODS

ECGI is a combined application of body surface electrocardiograms and computed tomography or magnetic resonance imaging data. Unipolar electrograms are reconstructed on the epicardial and endocardial surfaces. Activation time (AT) was defined as the time of maximal negative slope of the electrogram (EGM) during QRS, recovery time (RT) as the time of maximal positive slope of the EGM during T wave, Activation recovery interval (ARI) was defined as the difference between RT and AT. ARI dispersion (Δ ARI) and RT dispersion (Δ RT) were calculated as the difference between maximal and minimal ARI and RT respectively. We evaluated those parameters in patients with frequent PVCs from the RVOT, defined as >10.000 per 24 h, and in a control group.

RESULTS

We studied 7 patients with frequent RVOT PVCs and 17 controls. Patients with PVCs from the RVOT had shorter median RT than controls, in the endocardium and in the epicardium, respectively 380 (239-397) vs 414 (372-448) ms, p = 0.047 and 275 (236-301) vs 330 (263-418) ms, p = 0.047. The dispersion of ARI and of RT in the epicardium was higher than in controls, Δ ARI of 145 (68-216) vs 17 (3-48) ms, p = 0.001 and Δ RT of 201 (160-235) vs 115 (65-177), p = 0.019.

CONCLUSION

In this group of patients we found a shorter median RT in the endocardium and in the epicardium of the RVOT and a higher dispersion of the ARI and RT across the epicardium in patients with PVCs from the RVOT when comparing to controls.

Effects of Heart Rate and Ventricular Wall Thickness on Non-invasive Mapping: An in silico Study

Background: Non-invasive cardiac mapping—also known as Electrocardiographic imaging (ECGi)—is a novel, painless and relatively economic method to map the electrical activation and repolarization patterns of the heart, providing a valuable tool for early identification and diagnosis of conduction abnormalities and arrhythmias. Moreover, the ability to obtain information on cardiac electrical activity non-invasively using ECGi provides the potential for a priori information to guide invasive surgical procedures, improving success rates, and reducing procedure time.

Previous studies have shown the influence of clinical variables, such as heart rate, heart size, endocardial wall, and body composition on surface electrocardiogram (ECG) measurements. The influence of clinical variables on the ECG variability has provided information on cardiovascular control and its abnormalities in various pathologies. However, the effects of such clinical variables on the Body Surface Potential (BSP) and ECGi maps have yet to be systematically investigated.

Methods: In this study we investigated the effects of heart size, intracardiac thickness, and heart rate on BSP and ECGi maps using a previously-developed 3D electrophysiologically-detailed ventricles-torso model. The inverse solution was solved using the three different Tikhonov regularization methods.

Results: Through comparison of multiple measures of error/accuracy on the ECGi reconstructions, our results showed that using different heart geometries to solve the forward and inverse problems produced a larger estimated focal excitation location. An increase of ~2 mm in the Euclidean distance error was observed for an increase in the heart size. However, the estimation of the location of focal activity was still able to be obtained. Similarly, a Euclidean distance increase was observed when the order of regularization was reduced.

For the case of activation maps reconstructed at the same ectopic focus location but different heart rates, an increase in the errors and Euclidean distance was observed when the heart rate was increased.

Conclusions: Non-invasive cardiac mapping can still provide useful information about cardiac activation patterns for the cases when a different geometry is used for the inverse problem compared to the one used for the forward solution; rapid pacing rates can induce order-dependent errors in the accuracy of reconstruction.

Effects of ECG Signal Processing on the Inverse Problem of Electrocardiography

The inverse problem of electrocardiography is ill-posed. Errors in the model such as signal noise can impact the accuracy of reconstructed cardiac electrical activity. It is currently not known how sensitive the inverse problem is to signal processing techniques. To evaluate this, experimental data from a Langendorff-perfused pig heart (n=1) suspended in a human-shaped torso-tank was used. Different signal processing methods were applied to torso potentials recorded from 128 electrodes embedded in the tank surface. Processing methods were divided into three categories i) high-frequency noise removal ii) baseline drift removal and iii) signal averaging, culminating in n=72 different signal sets. For each signal set, the inverse problem was solved and reconstructed signals were compared to those directly recorded by the sock around the heart. ECG signal processing methods had a dramatic effect on reconstruction accuracy. In particular, removal of baseline drift significantly impacts the magnitude of reconstructed electrograms, while the presence of high-frequency noise impacts the activation time derived from these signals (p<0.05).

High-resolution noncontact charge-density mapping of endocardial activation

BACKGROUND

Spatial resolution in cardiac activation maps based on voltage measurement is limited by far-field interference. Precise characterization of electrical sources would resolve this limitation; however, practical charge-based cardiac mapping has not been achieved.

METHODS

A prototype algorithm, developed from first principles of electrostatic field theory, derives charge density (CD) as a spatial representation of the true sources of the cardiac field. The algorithm processes multiple, simultaneous, noncontact voltage measurements within the cardiac chamber to inversely derive the global distribution of CD sources across the endocardial surface.

RESULTS

Comparison of CD to an established computer-simulated model of atrial conduction demonstrated feasibility in terms of spatial, temporal, and morphologic metrics. Inverse reconstruction matched simulation with median spatial errors of 1.73 mm and 2.41 mm for CD and voltage, respectively. Median temporal error was less than 0.96 ms and morphologic correlation was greater than 0.90 for both CD and voltage. Activation patterns observed in human atrial flutter reproduced those established through contact maps, with a 4-fold improvement in resolution noted for CD over voltage. Global activation maps (charge density-based) are reported in atrial fibrillation with confirmed reduction of far-field interference. Arrhythmia cycle-length slowing and termination achieved through ablation of critical points demonstrated in the maps indicates both mechanistic and pathophysiological relevance.

CONCLUSION

Global maps of cardiac activation based on CD enable classification of conduction patterns and localized nonpulmonary vein therapeutic targets in atrial fibrillation. The measurement capabilities of the approach have roles spanning deep phenotyping to therapeutic application.

TRIAL REGISTRATION

ClinicalTrials.gov NCT01875614.

FUNDING

The National Institute for Health Research (NIHR) Translational Research Program at Royal Papworth Hospital and Acutus Medical.

Optical Imaging of Ventricular Action Potentials in a Torso Tank: A New Platform for Non-Invasive Electrocardiographic Imaging Validation

Background: Non-invasive electrocardiographic imaging (ECGI) is a promising tool to provide high-resolution panoramic imaging of cardiac electrical activity noninvasively from body surface potential measurements. Current experimental methods for ECGI validation are limited to comparison with unipolar electrograms and the relatively low spatial resolution of cardiac mapping arrays. We aim to develop a novel experimental set up combining a human shaped torso tank with high-resolution optical mapping allowing the validation of ECGI reconstructions.

Methods: Langendorff-perfused pig hearts (n = 3) were suspended in a human torso-shaped tank, with the left anterior descending artery (LAD) cannulated on a separate perfusion. Electrical signals were recorded from an 108-electrode epicardial sock and 128 electrodes embedded in the tank surface. Simultaneously, optical mapping of the heart was performed through the anterior surface of the tank. Recordings were made in sinus rhythm and ventricular pacing (n = 55), with activation and repolarization heterogeneities induced by perfusion of hot and cold solutions as well as Sotalol through the LAD. Fluoroscopy provided 3D cardiac and electrode geometries in the tank that were transformed to the 2D optical mapping window using an optimization algorithm. Epicardial unipolar electrograms were reconstructed from torso potentials using ECGI and validated using optical activation and repolarization maps.

Results: The transformation and alignment of the 3D geometries onto the 2D optical mapping window was good with an average correlation of 0.87 ± 0.10 and error of 7.7 ± 3.1 ms with activation derived from the sock. The difference in repolarization times were more substantial (error = 17.4 ± 3.7 ms) although the sock and optical repolarization patterns themselves were very similar (correlation = 0.83 ± 0.13). Validation of ECGI reconstructions revealed ECGI accurately captures the pattern of activation (correlation = 0.79 ± 0.11) and identified regions of late and/or early repolarization during different perfusions through LAD. ECGI also correctly demonstrated gradients in both activation and repolarization, although in some cases these were under or over-estimated or shifted slightly in space.

Conclusion: A novel experimental setup has been developed, combining a human-shaped torso tank with optical mapping, which can be effectively used in the validation of ECGI techniques; including the reconstruction of activation and repolarization patterns and gradients.

Impact of the Endocardium in a Parameter Optimization to Solve the Inverse Problem of Electrocardiography

Electrocardiographic imaging aims at reconstructing cardiac electrical events from electrical signals measured on the body surface. The most common approach relies on the inverse solution of the Laplace equation in the torso to reconstruct epicardial potential maps from body surface potential maps. Here we apply a method based on a parameter identification problem to reconstruct both activation and repolarization times. From an ansatz of action potential, based on the Mitchell-Schaeffer ionic model, we compute body surface potential signals. The inverse problem is reduced to the identification of the parameters of the Mitchell-Schaeffer model. We investigate whether solving the inverse problem with the endocardium improves the results or not. We solved the parameter identification problem on two different meshes: one with only the epicardium, and one with both the epicardium and the endocardium. We compared the results on both the heart (activation and repolarization times) and the torso. The comparison was done on validation data of sinus rhythm and ventricular pacing. We found similar results with both meshes in 6 cases out of 7: the presence of the endocardium slightly improved the activation times. This was the most visible on a sinus beat, leading to the conclusion that inclusion of the endocardium would be useful in situations where endo-epicardial gradients in activation or repolarization times play an important role.