Abstract
Background: Focal ventricular tachycardia (VT) is a life-threating arrhythmia, responsible for high morbidity rates and sudden cardiac death (SCD). Radiofrequency ablation is the only curative therapy against incessant VT; however, its success is dependent on accurate localization of its source, which is highly invasive and time-consuming.
Objective: The goal of our study is, as a proof of concept, to demonstrate the possibility of utilizing electrogram (EGM) recordings from cardiac implantable electronic devices (CIEDs). To achieve this, we utilize fast and accurate whole torso electrophysiological (EP) simulations in conjunction with convolutional neural networks (CNNs) to automate the localization of focal VTs using simulated EGMs.
Materials and Methods: A highly detailed 3D torso model was used to simulate ∼4000 focal VTs, evenly distributed across the left ventricle (LV), utilizing a rapid reaction-eikonal environment. Solutions were subsequently combined with lead field computations on the torso to derive accurate electrocardiograms (ECGs) and EGM traces, which were used as inputs to CNNs to localize focal sources. We compared the localization performance of a previously developed CNN architecture (Cartesian probability-based) with our novel CNN algorithm utilizing universal ventricular coordinates (UVCs).
Results: Implanted device EGMs successfully localized VT sources with localization error (8.74 mm) comparable to ECG-based localization (6.69 mm). Our novel UVC CNN architecture outperformed the existing Cartesian probability-based algorithm (errors = 4.06 mm and 8.07 mm for ECGs and EGMs, respectively). Overall, localization was relatively insensitive to noise and changes in body compositions; however, displacements in ECG electrodes and CIED leads caused performance to decrease (errors 16–25 mm).
Conclusion: EGM recordings from implanted devices may be used to successfully, and robustly, localize focal VT sources, and aid ablation planning.
Introduction
Ventricular tachycardia (VT) is a serious cardiac arrhythmia that represents an important source of morbidity and, upon degeneration into more lethal arrhythmias such as ventricular fibrillation (VF) (), sudden cardiac death (SCD) (; ). Hence, the prevention of VT, and its degeneration into VF, is of primary clinical importance to improve morbidity and reduce mortality.
In structurally healthy hearts, VT occurs primarily as a consequence of abnormal ectopic firing in the ventricles, overtaking sino-atrial activation and leading to premature ventricular contractions (PVCs). An effective treatment against ectopic VT is radiofrequency catheter ablation, which aims to target the tachycardia by first locating, and then electrically isolating the region causing the episode. However, procedure success is heavily dependent on an accurate localization of the VT source. Often, recordings of the focal VT, in the form of an electrocardiogram (ECG) or implanted device electrograms (EGM), exist prior to an ablation procedure, which inherently contain important information related to the focal origin of the VT source. Integration of computational studies and deep learning approaches provides an exciting opportunity to utilize the information contained within these recordings to potentially facilitate automated VT localization into clinical practice.
In recent decades, computational studies (; ; ; ; ) have enhanced greatly our knowledge of VT mechanisms and have strengthened diagnostic, therapeutic, and prognostic VT clinical strategies (; ; ; ), helping in the growth of personalized modeling (; ; ; ; ; ). One limitation of the majority of these studies is the dependence on monodomain formulations to represent electric sources in the form of transmembrane voltages. These models are time-consuming, and thus to achieve clinical translation, fast reaction-eikonal (RE) simulations (; ) have been the preferred choice. More recently, realistic simulations of full extracellular potentials at specific locations (e.g., ECG electrodes) have been obtained from the combination of lead field (LF) methods () with fast RE models (), achieving accuracy comparable to pseudodomain or bidomain formulations, but within a fraction of the computational time.
Using computational simulations of electrophysiological (EP) behavior has also been exploited to provide training datasets for machine and deep learning algorithms (; ); however, these studies did not utilize rapid RE models (; ), or LF methods (). were among the first to utilize convolutional neural networks (CNNs) to localize focal VT sources from simulated ECGs. The novelty of the study was in the integration of computational simulated data with CNN architectures; previous studies had in fact attempted to localize focal VTs from either simulated ECGs—utilizing myocardial activation imaging techniques—with no use of artificial intelligence ()—or clinical ECGs utilizing machine learning algorithms ().
One important limitation of was the restriction of the method to the use of ECGs. Although ECGs are widely used as a routine modality for VT management, they are not always available for VT patients, particularly focal VT patients in which the clinical VT is not inducible. Utilizing cardiac implantable electronic device (CIED) EGMs, which the majority of pre-ablation patients have in situ (; ), and which continuously record and store any abnormal arrhythmic activity, could bring great improvements to the automated localization of focal VT. Recent clinical studies have demonstrated that stored EGM recordings of re-entrant VT episodes from implanted devices can be successfully used to guide the construction of pace-maps during an ablation procedure, with similar accuracy to the use of ECGs, but with the advantage of ensuring that the clinical VT is targeted (; ). In our own recent work, we demonstrated how such EGM recordings might be utilized to perform patient-specific in silico pace-mapping (), improving pre-procedural ablation planning for complex scar-related VTs. However, literature lacks further investigation on the power of EGM recordings for the localization of focal VTs, for which computational models can address and answer a variety of different questions, and their use in AI-based algorithms.
In this study, we demonstrate the utility of leveraging the information contained within simulated implanted device EGM recordings for the automated localization of focal VTs in the LV. This could benefit clinical procedures by providing pre-procedural ablation information of the VT episodes without the necessity of acquiring ECG recordings of the focal VT, which represents the long-term aim of our study. Although the majority of idiopathic VTs originate in the right ventricular outflow track (RVOT) (), focal VTs can also originate from a variety of different locations in the LV (; ; ), and automating their localization could be beneficial to clinical procedures. To achieve our goal, we extend the previous work of and utilize fast computational simulations (RE combined with LF) on a realistic image-based torso model to generate ECG and EGM traces, which serve as inputs to a CNN architecture. We show the possibility of obtaining comparable localization in Cartesian coordinates between ECG-based and EGM-based trained CNNs. Moreover, we show improvement in the overall localization by introducing a novel CNN algorithm, utilizing a local ventricular-specific coordinate system ().
Materials and Methods
The workflow of this study is summarized in Figure 1. Briefly, this involved using a previously generated 3D torso model () (step 1) to rapidly simulate focal pacing across the LV within a RE environment (step 2). These simulated paced beats were combined with LF matrices computed on the standard 12-lead ECG electrodes and manufacturer-guided CIED right (RV) and left (LV) ventricular leads (step 3) to reproduce accurate 12-lead ECG and EGM recordings of the paced beats (step 4). The data were then processed and used as input to an adapted version of a previously developed CNN architecture by (step 5) and also to a novel network, consisting of a two-output regression and a classification CNN (step 6), and utilizing UVCs, to localize the paced beats (step 7).
FIGURE 1
Model Preparation
As in , a 3D torso model was generated from a computed tomography (CT) trans-catheter aortic valve implantation (TAVI) panning scan. The torso model included all major organs, with conductivities reported in Table 1, and a detailed four-chamber heart, extracted from a separate cardiac CT scan. The patient did not present any visible structural heart disease and consented to the use of their data in ethically approved research: UK Research Ethics Committee reference number 19/HRA/0502 and 15/LO/1803. To decrease computational time without a loss of physiological electric signals, the average ventricular edge length of the biventricular mesh was kept at 738 μm. Realistic fiber orientation was incorporated into the ventricular myocardium using a well-established rule-based approach ().
TABLE 1
| Organs | Tissue conductivities (S/m) |
| Lungs | 0.0714 |
| Bones | 0.05 |
| Skin | 0.05 |
| Fat/Muscle | 0.24725 |
| Liver | 0.1667 |
| Spleen | 0.1 |
| Kidneys | 0.1667 |
| Aorta | 0.6667 |
| Ventricular blood pools | 0.6667 |
| Atrial blood pools and walls | 0.6667 |
| Pericardium | 0.2 |
Organ conductivities of our torso model.
See also .
The LV was geometrically divided into 17 segments, according to the American Heart Association (AHA) guidelines (), as shown in Figures 2A–F. In addition, each of the segments was subsequently divided into four, for a total of 68 (Figures 2G,H). These models were used as guidance for the collection of pacing locations, for the generation of training and testing labels for the existing CNNs, and the visualization of the localized VT sources.
FIGURE 2
To replicate focal ectopic VTs across the LV segments, ∼3767 randomly chosen paced beats—single stimuli, with a basic cycle length (BCL) of 400 ms—were simulated using a computationally efficient RE formulation () within the Cardiac Arrhythmia Research Package (CARP) (), utilizing the 10 Tusscher ventricular cell model (). Intra- and extracellular tissue conductivities were tuned to achieve physiological QRSs (), comparable to equivalent pseudo bidomain simulations on a higher-resolution mesh (). Intra- and extracellular conductivities were 0.1845S/m and 0.6628S/m along the fiber direction, respectively, and 0.0493S/m and 0.1769S/m transverse to it. The corresponding RE conduction velocities (CVs) were 0.5455m/s and 0.1802m/s, along and transverse the fiber direction, respectively.
To allow the computation of extracellular potential signals from specific locations within the torso (Figure 3), the simulated cardiac potentials of each paced beat were combined with the LF Method (). Specifically, LF matrices were calculated within CARP () on the standard ECG lead locations and on the RV and LV lead sensing parts of a standard Boston Scientific implanted device (). This virtual device had a non-septal RV lead, with a superior vena cava (SVC) coil in the right atrium (RA), and a straight LV lead through the coronary sinus, with four sensing LV tips distanced equally at 7.5mm. Configurations of both 12-lead ECG and CIED are shown in Figures 4A–D. All sensing electrodes, including the can of the implanted device (CAN), were approximated to single points, to increase the speed of LF computations and subsequent simulations. The computation of these matrices was only performed once for each torso configuration and took ∼8 min (128 cores). Their combination with the RE solutions produced high-fidelity 12-lead ECGs and EGM traces (Figures 4E,F) in ∼20 s (256 cores) for each paced location. Eight EGM vectors were chosen as the main EGM signals (), and included far-field CAN-SVC, CAN-RV, and SVC-RV, and near-field RV tip-RV ring and each LV tip-RV tip. However, importantly, additional vector combinations (four for ECGs and eight for EGMs) were added to the standard signals to facilitate integration into the CNN algorithms (see section “CNNs Training and Testing”).
FIGURE 3
FIGURE 4
Finally, a standardized universal ventricular coordinate (UVC) system was computed on the biventricular mesh () to facilitate the development of a novel CNN specific to the ventricles, which should be advantageous as it identifies and constrains the localization of the paced beats inherently within the myocardium. As shown in Figure 5, UVCs describe the biventricular mesh using three parameters: z—normalized distance between apex (0) and base (1) along the long axis, ρ —normalized distance between endocardial (0) and epicardial (1) surfaces along the short axis, and φ —rotational distance from LV septum.
FIGURE 5
CNN Architectures
In this study, we developed two separate 2D CNN architectures, which used the same ECG and EGM traces as inputs to identify the location of a simulated paced beat (representing an ectopic VT). The first architecture, based on , locates the origin of the paced beat in Cartesian coordinates, after converting the outputs of the CNNs. The second utilizes a regression and a classification CNN to locate the VT in UVC space, naturally constraining the final localization of the focal VT source to the myocardium.
The existing Cartesian probability-based architecture was reproduced from and is composed of two classification CNNs named Segment and EpiEndo CNNs. Segment CNN classified in which LV segment the pacing beat originated, whereas EpiEndo CNN determined whether the pacing was endocardial or epicardial (binary decision). In our study, we developed two separate Segment CNNs: one classified between 17 LV segments (CNN with 17 output neurons) and the other between 68 LV segments (68 output neurons). Briefly, the structure of both Segment and EpiEndo CNNs consisted of two hidden layers alternating with two pooling layers and terminating with a fully connected (FC) layer. The output of both Segment and EpiEndo CNNs was a probability distribution (likelihood of each output neuron being the correct class). These distributions were obtained utilizing a softmax function on the output of the final FC layer. As performed in , the probability distributions and for Segment and EpiEndo, respectively) of each output segment (largest probability) and its adjacent segments were then combined with the centers of gravity CoG(x, y, z)ij of the corresponding endocardial and epicardial surfaces, as shown in Equation 1, to localize a paced beat in Cartesian coordinates.
The distance of the localized sources to the ground truths (simulated sources) was expressed in terms of localization errors (computed as Euclidean distance in millimeters).
Our novel UVC-based algorithm is composed of one regression CNN, outputting z and ρ, and one 68-feature classification CNN, predicting the rotational coordinate φ. The structure of both CNNs was similar to the Cartesian probability-based network (hidden layer–hidden layer–pooling layer–FC layer), as shown in Figure 6. Because of the cyclic nature of φ, a three-output regression would have not returned satisfying and accurate results; hence, we used φ to divide the LV into 68 “wedges” (φ was grouped into intervals of 0.09 radians with each class assigned a label from 1 to 68, spanning φ = −π to φ=π). Using a higher number of LV “wedges” would have not returned desirable accuracies; thus, we decided to use a number of features (68) that had worked for Segment CNN and that was still suitable for achieving precise localization along φ. For the final localization of the paced beats, the outputs of the 68-feature classification (“wedges” with the highest probabilities) were converted back to φ, and combined with z and ρ regression predictions.
FIGURE 6
Both Cartesian probability-based and UVC-based algorithms were implemented in Python using Scikit-learn () and TensorFlow ().
CNN Localization Performance
Localization performance of the Cartesian probability-based algorithm was expressed in terms of localization error in millimeters, as described above. The same metric was used for our novel UVC-based architecture; however, the predicted values were first transformed from UVC space to Cartesian space [by locating the closest node in the mesh, with appropriate scaling of UVC coordinates ()] and then the distance with the ground truths (simulated pace beats) was evaluated (Euclidian distance, in millimeters).
For the Segment CNN of the Cartesian probability-based architecture, testing performance was evaluated in terms of accuracy, to allow comparison with results from . Accuracy is defined as the percentage of paced beats correctly classified within each segment.
CNNs Training and Testing
Training and testing inputs of both Cartesian probability-based and UVC-based CNNs were ECG and EGM traces computed from 3767 pacing excitations across the LV. To facilitate the execution of the 2D CNNs, the ECG and EGM signals had to be placed in square matrices; hence, we added four additional leads to the standard 12-lead ECGs [LL-RA-LA, LA-RA-LL, RA-LA-LL, and (RA+LA+LL)/2], as performed in , and eight more EGM vectors to the standard eight EGMs (CAN-each LV tip and SVC-each LV tip), to achieve a total of 16 combinations of ECG and EGM vectors. QRSs were then extracted and sampled in time (16 time points), and stacked in 16 × 16 matrices, as shown in Figures 4E,F. A total of 2767 sets of these ECG and EGM matrices were used for training and were uniformly distributed across the myocardium (∼36% intramural/mid-wall, ∼32% epicardial, and ∼32% endocardial), with the exception of EpiEndo training data, which were epicardial and endocardial only (). White Gaussian noise with a signal-to-noise ratio (SNR) of 25 dB was then added 10 consecutive times to all 16 ECG and EGM leads of the training set to augment the data by 10-fold (∼27,670) and increase robustness of the CNN training. A 10-fold cross-validation was performed in the existing Cartesian probability-based CNNs as part of the training (), with a 90% (training)–10% (validation/testing) split. On the other hand, the cross-validation was used for hyper-parameter tuning in the UVC-based networks. After training, the localization performance of both Cartesian probability-based and UVC-based networks was tested by feeding the retained 1000 sets of ECG and EGM QRSs, with a SNR of 25 dB. Parameters of both Segment and EpiEndo CNNs were taken from ; batch size was set to 23, number of epochs was set to 10, learning rate was set to 0.001, and cross-entropy was used as loss function. A ReLU function was used as the activation function for feed-forward propagation, and a gradient-descent-projection method was used as the back propagation algorithm. In our UVC-based networks, we used similar parameters, except for the regression where we set the batch size to 23 and the number of epochs to 15, and we used mean absolute error as loss function.
Investigation of Model Uncertainties
Localization performance of both Cartesian probability-based and UVC-based CNN architectures, trained on the data described above, was also investigated by introducing different noise levels to the retained 1000 sets (SNR = 5, 10, 15, 20, and 30 dB). Moreover, we investigated the localization performance of both architectures as body compositions of the torso model were also varied, shown in Table 2, as well as different ECG electrode configurations (Figure 7) and different CIED configurations (Figure 8). For all these variations, LF matrices were recomputed (according to the new organ conductivities or electrodes positions) and combined with the retained 1000 intramural excitations to obtain new ECG and EGM matrices. These traces were then used to test both previously trained CNN architectures. Some of the major organ conductivities were varied according to physiological variations (; ); however, we chose to pair specific changes (for instance, liver and lungs, fat/muscle, named “bath” and liver, and different blood pools, etc.) to challenge CNN localization performance. ECG electrodes were displaced by 5 cm in all major orthogonal directions, and across all leads. Specifically, we shifted all ECGs leads upward (Figure 7A) and downward (Figure 7B)—RA and LA were always shifted downward, and LL upward—toward the left (Figure 7C) and the right (Figure 7D). Moreover, in one configuration (Figure 7E), the distance between ECG leads was increased by ∼10 cm. Finally, we simulated variations in electrode location and diameter of the virtual implanted device, as reported in for different cardiac-resynchronization therapy (CRT-D) devices available in the market. Specifically, we changed the spacing between the sensing electrodes of the straight LV lead, to account for shorter or longer inter-electrode distance; in addition, we increased the diameter of RV and LV tips to ∼2 mm. In the latter scenario, instead of considering the EGM signals from single point electrodes, we averaged the signals obtained from a cloud of points within a 2-mm radius, to simulate more realistic conditions, and investigate whether our single point approximation of the CIED leads could affect the final localization. Finally, we considered the case of a septal RV coil configuration, which has been tested in CRT-D () and cardioverter defibrillators (ICD) ().
TABLE 2
| (1) Liver: 0.023 (S/m) and lungs: 0.039 (S/m) |
| (2) Bath: 0.45 (S/m) (pure muscle) and lungs: 0.039 (S/m) |
| (3) Bath: 0.05 (S/m) (pure fat) and lungs: 0.203 |
| (4) Liver: 0.2 (S/m) and lungs: 0.039 (S/m) |
| (5) Bath: 0.05 (S/m) and lungs: 0.039 (S/m) |
| (6) Bath: 0.45 (S/m) and lungs: 0.203 (S/m) |
| (7) Bath (all organs except lungs): 0.24 (S/m) and lungs: 0.07 (S/m) |
| (8) Skin: 0.117 (S/m) |
| (9) Atria, ventricles, and aorta: 0.84 (S/m) |
Variations in body compositions.
Combinations of different organ conductivities within physiological changes are here reported, used in our CNNs sensitivity analysis.
FIGURE 7
FIGURE 8
Results
Utility of EGMs in Existing Segment/EpiEndo-Based CNNs
We successfully reproduced the existing classification CNNs, namely, Segment and EpiEndo, introduced by , to be trained and tested not only on ECG traces, but also on 16 different combinations of EGM vectors from a standard CIED with RV and LV leads. Testing performance of Segment CNN was similar for both ECG-based and EGM-based testing. As shown in Figure 9A, 86.76% accuracy was achieved for ECG-based testing and 79.70% was achieved for EGM-based testing (SNR = 25 dB).
FIGURE 9
The precision of Segment CNN in each segment, which defined how correct the CNN is at classifying one segment, is shown in Figure 9B. Here, we see that ECGs and EGMs have a similar influence on the network in almost every segment, with only few exceptions. The three highest precisions are in segments 1, 10, and 14 for ECG-based testing, and 1, 10, and 16 for EGM-based testing. The three lowest are in 2, 3, and 15 for ECGs, and 2, 7, and 9 for EGMs.
Utility of EGMs in Cartesian Probability-Based Localization
Localization in Cartesian space of each paced beat, from either ECG or EGM signals, was possible by combining probability distributions of Segment and EpiEndo CNNs (as shown in Equation 1). The localization performance, defined as the Euclidean distance in millimeters between an estimated source and the real location of the simulated paced beat, for the testing dataset of 1000 cases, is reported in Figure 10 for ECG-based and EGM-based testing. ECG-based localization and EGM-based localization produced a mean localization error of 11.76 ± 5.32 mm and 13.25 ± 6.79 mm, respectively.
FIGURE 10
Application of the 68-feature Segment CNN, based on the 68-segment LV AHA model shown in Figure 2G, was able to reduce localization errors of both ECGs and EGMs to 6.69 ± 3.19 mm and 8.74 ± 6.41 mm, respectively, as shown in Figure 10.
UVC-Based Localization
Further improvements in localization performance were made by developing two CNNs, which returned the position of a paced beat in a reference frame specific to the ventricles (UVCs). This UVC-based localization outperformed the Cartesian probability-based localization, as shown in Figure 11A, reducing localization errors to 4.06 ± 2.47 mm and 8.07 ± 8.26 mm for ECG and EGM-based testing, respectively.
FIGURE 11
UVC-based localized sources are visualized in 2D in a patient-specific bull’s-eye diagram, shown in Figure 11B for 30 beats, as previously illustrated in Figure 5. Here, a paced beat can be visualized using its UVC coordinates and can be compared to the ground truth, revealing a close match between all pairs. The radius of the diagram describes the distance of a paced beat to the LV apex (center of the diagram), relatable to UVC z, and its circumferential direction (φ) facilitates the positioning of the beat within a specific segment. The intramural location (ρ) of the beat (how far from the endocardium and/or epicardium) is color coded.
Sensitivity to Noise
Overall, localization was only slightly affected by noise, as seen in Figure 12A (ECG-based localization) and Figure 12B (EGM-based localization). As SNR decreased (increased noise), localization errors increased only slightly, with one exception (SNR = 5 dB), where the performance of both UVC-based and Cartesian probability-based localization was reduced. However, all localization errors were < 12.5 mm for ECG-based localization. Moreover, noise seemed to affect EGM-based localization more than ECG-based localization.
FIGURE 12
Sensitivity to Electrode Locations
Displacements of ECG leads and different CIED configurations did affect the localization performance of UVC-based and Cartesian probability-based algorithms (errors > ∼15 mm) for ECG- and EGM-based testing, as shown in Figures 12C,D, respectively. ECG-based localization was more affected by displacements away from the heart (20 mm)—right and downward shifts. Errors in EGM-based localization were higher (20 mm) when considering longer inter-electrode distance (20 mm) and increased lead surface diameter (2 mm). For UVC-based localization, a septal RV coil configuration caused errors to increase > 20 mm as well.
Sensitivity to Tissue Conductivities in Torso Model
A comparison between ECG-based and EGM-based localization for different body compositions is shown in Figure 13A (for UVC-based localization) and Figure 13B (for Cartesian probability-based localization). ECG-based localization was only affected by a high increase of fat in the torso bath (scenarios 3 and 5) and when the whole torso was simplified to bath and lungs (scenario 7). In those three scenarios, mean localization errors increased to 17.75 ± 9.88 mm, 20.72 ± 10.99 mm, and 14.08 ± 7.38 mm for UVC-based testing, respectively, and to 13.01 ± 8.89 mm, 15.07 ± 11.20 mm, and 13.88 ± 8.40 mm for 68-segment Cartesian probability-based testing. Other variations of tissue conductivity did not affect the performance of either algorithm (localization errors < ∼8 mm).
FIGURE 13
Finally, EGM-based testing was less affected by changes in body compositions during Cartesian probability-based localization. Similar to above for ECG-based localization, simplification of the torso to bath and lungs (scenario 7) caused the highest mean error (14.83 ± 11.24 mm), but in all other scenarios, localization ranged between 9.04 and 10.66 mm. On the other hand, UVC-based EGM localization had a similar trend to ECG-based localizations, with errors > 15 mm for scenarios 3, 5, and 7.
Discussion
In this study, we successfully utilized simulated implanted device EGMs to localize virtual focal VT sources using CNNs, achieving accuracies that could be useful in clinical settings. A previous algorithm () utilized 12-lead ECGs for a similar purpose; here, we managed to replicate the structure of the CNN architecture for EGM traces and improve the overall localization by introducing a higher number of segments in the AHA LV model. Moreover, we also improved the overall localization precision by introducing a novel architecture composed of regression and classification algorithms, which was able to identify the source in a framework specific to the ventricles, easily interpretable by clinicians. Finally, we investigated the robustness of both CNN algorithms to the introduction of uncertainties, such as different noise levels in the data, and possible inter-patient variabilities (different body compositions, ECG lead positions, and CIED configurations).
Successful Application of EGMs in Existing Cartesian Probability-Based Algorithm
Simulated focal VT sources were successfully identified from 16 combinations of implanted device EGM vectors. In the previous study, achieved localization precision in the range of 10–11 mm when utilizing a combination of 16 ECG vectors; here, we reduced localization errors to 6.69 mm and 8.74 mm for ECG and EGM traces, respectively, by incrementing the number of segments in the LV to 68. In clinical practice, the average diameters of catheter tips are between 4 mm and 8 mm (), limiting the average lesion size to a minimum of ∼8.5 mm (). Hence, we achieved localization precisions in a range suitable for improving ablation planning. Especially in patients with a non-stable condition, pre-planning of these procedures could be expedited and aided if the acquisition of ECG data during VT would not be required, which can be achieved by utilizing information stored in implanted devices. Our algorithm thus proposes a first level of investigation that could direct clinicians to the region of interest with high precision. Moreover, we achieved ECG-based and EGM-based testing accuracies in ranges comparable to (77%). Similar patterns to the previous study were seen when investigating how noise affects the localization; only a loss in accuracy and localization precision is seen with SNR = 5 dB. Interestingly, noise seems to have a greater impact on EGM-based localization than on ECGs. This could be explained by the fact that implanted device sensing vectors are closer to one another and to the cardiac electrical activity, amplifying variations caused by noise, thus affecting EGM-based localization to a greater degree.
Novel UVC-Based Algorithm Improves Localization
Our novel UVC-based algorithm improves localization to 4.06 mm and 8.07 mm for ECGs and EGMs, respectively, outperforming the existing study. Whereas the Cartesian probability-based algorithm relies on combining probabilities of two networks with the geometrical centers of gravity of each LV segment to locate a VT origin, our architecture predicts the actual location of the source in terms of its (normalized) distance from the apex, the LV septum and, most importantly, the endocardium. Furthermore, it intrinsically bounds the localization to the myocardium. Knowing the exact intramural (mid-wall) location, a VT source could help in the choice of power, tip diameter, and lesion size to apply, as well as access direction (epicardial or endocardial), in pre-procedural planning. Finally, our novel localization facilitates the visualization of focal estimates, by plotting a patient-specific bull’s-eye diagram, where the radius represents the distance from the apex and the circumferential direction relates to septal, anterior, inferior, and lateral LV segments.
Automated Localization Is Only Affected by Extreme Changes in ECG Lead Positions and Implanted Device Lead Configurations
Focal VT localization is only marginally affected by differences in body compositions. However, to increase the accuracy of the results, a torso model constructed for algorithm training should at least include all major thoracic and abdominal organs with realistic conductivity values; our findings suggest that representing bath and lungs (as used in Yang et al) only produce signals that differ substantially from more complete torso models, importantly affecting localization accuracy. Moreover, EGM-based algorithms seem to be more robust to tissue variations, possibly due to the closer positioning of the device leads to the electrical cardiac source, with extracellular potentials being less affected by the surrounding tissue/organ conductivities. Displacements in ECG leads and differences in common CIED configurations do seem to have an impact on the final localization; this suggests the necessity of integrating a higher variability in the training data, or extrapolating ECG or CIED patient-specific information from imaging data to strengthen future automated algorithms and allow clinical validation and translation.
Limitations
A notable limitation of this study is the absence of clinical validation. However, our main goal was to strengthen the automated localization of focal VTs and investigate the properties of our 3D pipeline that need improvement for future clinical studies. For future validation of our in silico EGM model and corresponding CNN localization, we will need to generate patient-specific 3D models that have been registered and tuned to the clinical framework used during EP mapping and ablation, collect simulated data on such models for CNN training, and test the latter on clinical EGM recordings of the focal VT(s), and/or paced beats, that have been collected from CIEDs directly or from the latter recording during the mapping. When attempting clinical translation in the future, we will also investigate other aspects of our work regarding patient-specific EP properties that were not taken into consideration in this study. Our model required certain simplifications, such as rule-based fibers and lack of Purkinje activation, which we believe would not make an impact in the final performance of our algorithms when dealing with focal beats, but that could be useful to take into account for more complex patient-specific approaches. Although our cardiac model was static, we do not believe that the absence of electro-mechanical feedback significantly influenced the final ECG or EGM signals, when considering only QRSs (ventricular activation); many studies have validated static simulated EGM signals against clinical data (; ; ), showing that it is not necessary to couple mechanical simulations with EP for these types of problems. Moreover, we only considered single beats originating in the LV. In future studies, it will be worth including focal VTs in the RV, which is a common region of VT especially around the outflow track (RVOT). This could be easily achieved by using the UVC system, which covers the RV, to generate labels and prepare simulations, facilitating both modeling and localization pipelines. Although we believe that simulating multiple paced beats would not have an impact in the final CNN performance and localization, it will be necessary to achieve more realistic scenarios, as it can influence the waveforms of ECG and EGM traces. Furthermore, extending the automated localization of VT to more complex episodes (for instance, in presence of micro re-entries and/or infarction) represents an interest of ours that will be addressed in future studies. The investigation on how different signal uncertainties influence the performance of our CNNs could also be extended to include more complex and realistic ways of adding noise to customize computational models to patient-specific settings (; ). Another aspect of this study that could be refined is the overall structure of our novel UVC-based architecture; both regression and classification networks were implemented following the structure proposed by , although some parameters were optimized to fit the new tasks. In future studies, deeper networks could be developed, and different input data shapes could be investigated (e.g., 2D vs. 1D). Moreover, to tackle the problem of computational efficiency and decrease even further our simulation time when dealing with more complex arrhythmias, we may investigate the possibility of GPU-based models, which have recently opened new perspectives in terms of real-time, physiologically detailed simulations ().
Conclusion
By integrating fast EP simulations with deep-learning algorithms, we have demonstrated the utility of our in silico pipeline for the simulations of EGMs stored in implanted devices, which, in addition to 12-lead ECGs, can accurately localize focal VTs. Our novel in silico automated algorithm, which utilizes a coordinate frame specific to the ventricles, increased localization precision above previous segment-classification approaches, facilitating clinical interpretation. Moreover, we showed the necessity of including more variability in the training data regarding lead positions, and the stability, on the other hand, of the localization to changes in body compositions.
Disclaimer
The views expressed are those of the author(s) and not necessarily those of the National Health Service, the National Institute for Health Research, or the Department of Health.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
SM, KG, and EP-A were primarily led the work and consulted in the initial coding of lead field computations and CNN structure, respectively. RR provided the CT data used for torso model creation. GP was responsible for developing and providing the reaction-eikonal and lead field packages within CARP. Finally, AK and MB supervised the project, guided the study design, and provided useful insights into and feedback on every aspect of the work (from the modeling side to the localization pipeline and the CNN algorithms). All authors contributed to the article and approved the submitted version.
Funding
SM was funded by the Engineering and Physical Sciences Research Council (EPSRC; EP/L015226/1). This work was supported by the National Institute for Health Research Biomedical Research Centre at Guy’s and St. Thomas’ Trust and King’s College, and the Centre of Excellence in Medical Engineering funded by the Wellcome Trust and Engineering and Physical Sciences Research Council (EPSRC; WT 088641/Z/09/Z). MB was also supported by a Medical Research Council New Investigator Grant (MR/N011007/1).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AbadiM.AgarwalA.BarhamP.BrevdoE.ChenZ.CitroC.et al (2015). TensorFlow: Large-Scale Machine Learning on Heterogeneous Systems. Available online at: tensorflow.org(accessed Febraury 17, 2021).
2
AntoniadisA. P.BeharJ. M.SieniewiczB.GouldJ.NiedererS.RinaldiC. A. (2017). A comparison of the different features of quadripolar left ventricular pacing leads to deliver cardiac resynchronization therapy.Exp. Rev. Med. Dev.14697–706. 10.1080/17434440.2017.1369404
3
BaroneA.GizziA.FentonF.FilippiS.VenezianiA. (2020). Experimental validation of a variational data assimilation procedure for estimating space-dependent cardiac conductivities.Comput. Methods Appl. Mech. Eng.358:112615. 10.1016/j.cma.2019.112615
4
BayerJ. D.BlakeR. C.PlankG.TrayanovaN. A. (2012). A novel rule-based algorithm for assigning myocardial fiber orientation to computational heart models.Annal. Biomed. Eng.402243–2254. 10.1007/s10439-012-0593-5
5
BayerJ.PrasslA. J.PashaeiA.GomezJ. F.FronteraA.NeicA.et al (2018). Universal ventricular coordinates: a generic framework for describing position within the heart and transferring data.Med. Image Anal.4583–93. 10.1016/j.media.2018.01.005
6
Cardone-NoottL.Bueno-OrovioA.MincholéA.ZemzemiN.RodriguezB. (2016). Human ventricular activation sequence and the simulation of the electrocardiographic qrs complex and its variability in healthy and intraventricular block conditions.Europace18 (Suppl. 4)iv4–iv15. 10.1093/europace/euw346
7
CedilnikN.DuchateauJ.DuboisR.SacherF.JaïsP.CochetH.et al (2018). Fast personalized electrophysiological models from computed tomography images for ventricular tachycardia ablation planning.Europace20(Suppl. 3)III94–III101. 10.1093/europace/euy228
8
CedilnikN.SermesantM.EikonalM. S. (2019). Model Personalisation Using Invasive Data to Predict Cardiac Resynchronisation Therapy Electrophysiological Response. Available online at: https://hal.inria.fr/hal-02368288(accessed October 20, 2020).
9
ClaytonR. H. H.BishopM. J. J. (2014). Computational models of ventricular arrhythmia mechanisms: recent developments and future prospects.Drug Discov. Today Dis. Models4:1016. 10.1016/j.ddmod.2014.04.002
10
CostaC. M.HoetzlE.RochaB. M.PrasslA. J.PlankG. (2013). Automatic parameterization strategy for cardiac electrophysiology simulations.Comput. Cardiol.40373–376.
11
GilletteK.MatthiasA. F. G.PrasslA. J.KarabelasE.ReiterU.ReiterG.et al (2021). A framework for the generation of digital twins of cardiac electrophysiology from clinical 12-Leads ECGs.Med. Image Anal.71:102080. 10.1016/j.media.2021.102080
12
HarrisP.LysitsasD. (2016). Ventricular arrhythmias and sudden cardiac death.BJA Educ.16221–229. 10.1016/j.ccep.2017.08.004
13
HenriquezC. S. (2014). A brief history of tissue models for cardiac electrophysiology.IEEE Transact. Biomed. Eng.611457–1465. 10.1109/TBME.2014.2310515
14
IlgK.BamanT. S.GuptaS. K.SwansonS.GoodE.ChughA.et al (2010). Assessment of radiofrequency ablation lesions by CMR imaging after ablation of idiopathic ventricular arrhythmias.JACC Cardiovas. Imaging3278–285. 10.1016/j.jcmg.2009.09.028
15
ItoS.TadaH.NaitoS.KurosakiK.UedaM.HoshizakiH.et al (2003). Development and validation of an ECG algorithm for identifying the optimal ablation site for idiopathic ventricular outflow tract tachycardia.J. Cardiovasc. Electrophysiol.141280–1286. 10.1046/j.1540-8167.2003.03211.x
16
Le BrasA. (2018). 3D virtual heart to guide VT ablation.Nat. Rev. Cardiol.15:654. 10.1038/s41569-018-0102-0
17
LeclercqC.SadoulN.MontL.DefayeP.OscaJ.MoutonE.et al (2016). Comparison of right ventricular septal pacing and right ventricular apical pacing in patients receiving cardiac resynchronization therapy defibrillators: the SEPTAL CRT study.Eur. Heart J.37473–483. 10.1093/eurheartj/ehv422
18
MarcotteC. D.FentonF. H.HoffmanM. J.CherryE. M. (2021). Robust data assimilation with noise: applications to cardiac dynamics.Chaos31:013118. 10.1063/5.0033539
19
Medtronic (2016). Cardioinsight Mapping Vest. Available online at: https://www.medtronic.com/us-en/healthcare-professionals/products/cardiac-rhythm/cardiac-mapping/cardioinsight-mapping-vest.html(accessed November 23, 2020).
20
Mendonca CostaC.NeicA.KerfootE.PorterB.SieniewiczB.GouldJ.et al (2019). Pacing in proximity to scar during cardiac resynchronization therapy increases local dispersion of repolarization and susceptibility to ventricular arrhythmogenesis.Heart Rhythm161475–1483. 10.1016/j.hrthm.2019.03.027
21
MonaciS.StrocchiM.RoderoC.GilletteK.WhitakerJ.RajaniR.et al (2020). In-silico pace-mapping using a detailed whole torso model and implanted electronic device electrograms for more efficient ablation planning.Comput. Biol. Med.125:104005. 10.1016/j.compbiomed.2020.104005
22
NeicA.CamposF. O.PrasslA. J.NiedererS. A.BishopM. J.VigmondE. J.et al (2017). Efficient computation of electrograms and ECGs in human whole heart simulations using a reaction-eikonal model.J. Comput. Phys.346191–211. 10.1016/j.jcp.2017.06.020
23
NiedererS. A.LumensJ.TrayanovaN. A. (2019). Computational models in cardiology.Nat. Rev. Cardiol.16100–111. 10.1038/s41569-018-0104-y
24
PedregosaF.MichelV.GriselO.BlondelM.PrettenhoferP.WeissR.et al (2011). Scikit-learn: machine learning in python gaël varoquaux bertrand thirion vincent dubourg alexandre passos pedregosa, varoquaux, gramfort et al. matthieu perrot.J. Mach. Learn. Res.12:54474.
25
Pekka RaatikainenM. J.ArnarD. O.ZeppenfeldK.MerinoJ. L.LevyaF.HindriksG.et al (2014). Statistics on the use of cardiac electronic devices and electrophysiological procedures in the european society of cardiology countries: 2014 report from the european heart rhythm association.Europace17(Suppl. 1)i1–i75. 10.1093/europace/euu300
26
PlanckeA. M.AdamC. P. M.NeicG. A.McSpaddenL. C.WhitakerJ.O’NeillM.et al (2019). Generation of a cohort of whole-torso cardiac models for assessing the utility of a novel computed shock vector efficiency metric for ICD optimisation.Comput. Biol. Med.112:103368. 10.1016/j.compbiomed.2019.103368
27
PotseM. (2018). Scalable and accurate ECG simulation for reaction-diffusion models of the human heart.Front. Physiol.9:370. 10.3389/fphys.2018.00370
28
PrasslA. J.KickingerF.AhammerH.GrauV.SchneiderJ. E.HoferE.et al (2009). Automatically generated, anatomically accurate meshes for cardiac electrophysiology problems.IEEE Transact. Biomed. Eng.561318–1330. 10.1109/TBME.2009.2014243
29
RantnerL. J.VadakkumpadanF.SpevakP. J.CrossonJ. E.TrayanovaN. A. (2013). Placement of implantable cardioverter-defibrillators in paediatric and congenital heart defect patients: a pipeline for model generation and simulation prediction of optimal configurations.J. Physiol.5914321–4334. 10.1113/jphysiol.2013.255109
30
RelanJ.ChinchapatnamP.SermesantM.RhodeK.GinksM.DelingetteH.et al (2011). Coupled personalization of cardiac electrophysiology models for prediction of ischaemic ventricular tachycardia.Interface Focus1396–407. 10.1098/rsfs.2010.0041
31
RitchieH.RoserM. (2018). “Causes of Death.” Our World in Data. Available online at: https://ourworldindata.org/causes-of-death(accessed November 12, 2020).
32
SelvaduraiB. S. N.PuntmannV. O.BluemkeD. A.FerrariV. A.FriedrichM. G.KramerC. M.et al (2018). Definition of left ventricular segments for cardiac magnetic resonance imaging.JACC Cardiovas. Imag.11926–928. 10.1016/j.jcmg.2017.09.010
33
ShadeJ. K.AliR. L.BasileD.PopescuD.AkhtarT.MarineJ. E.et al (2020). Preprocedure application of machine learning and mechanistic simulations predicts likelihood of paroxysmal atrial fibrillation recurrence following pulmonary vein isolation.Circulation Arrhyth. Electrophysiol.13617–627. 10.1161/CIRCEP.119.008213
34
SoviljS.MagjarevićR.AbedA. A.LovellN. H.DokosS. (2014). Simplified 2D bidomain model of whole heart electrical activity and ECG generation.Measur. Sci. Rev.14136–143. 10.2478/msr-2014-0018
35
SrinivasanN. T.SchillingR. J. (2018). Sudden cardiac death and arrhythmias.Arrhyth. Electrophysiol. Rev.7:111. 10.15420/aer.2018:15:2
36
SrivathsanK.LesterS. J.AppletonC. P.ScottL. R. P.MungerT. M. (2005). Ventricular Tachycardia in the Absence of Structural Heart Disease: Indian Pacing and Electrophysiology Journal.Amsterdam: Elsevier. 10.1016/b978-1-4377-2029-7.00012-7
37
ten TusscherK. H. W. J. D.NobleP. J.PanfilovA. V. (2004). A model for human ventricular tissue.Am. J. Physiol. Heart Circulat. Physiol.286H1573–H1589. 10.1152/ajpheart.00794.2003
38
TrakicA.AkhandM.WangH.MasonD.LiuF.WilsonS.et al (2010). Computational modelling of blood-flow-induced changes in blood electrical conductivity and its contribution to the impedance cardiogram.Physiol. Measur.3113–33. 10.1088/0967-3334/31/1/002
39
TrayanovaN. A. (2011). Whole-heart modeling: applications to cardiac electrophysiology and electromechanics.Circulat. Res.108113–128. 10.1161/CIRCRESAHA.110.223610
40
TrayanovaN. A.PashakhanlooF.WuK. C.HalperinH. R. (2017). Imaging-based simulations for predicting sudden death and guiding ventricular tachycardia ablation.Circulat. Arrhyth. Electrophysiol.10:4743. 10.1161/CIRCEP.117.004743
41
VasconcellosE. C.CluaE. W. G.FentonF. H.ZamithM. (2020). Accelerating simulations of cardiac electrical dynamics through a multi-GPU platform and an optimized data structure.Concur. Comput.32:e5528. 10.1002/cpe.5528
42
VigmondE. J.HughesM.PlankG.LeonL. J. (2003). Computational tools for modeling electrical activity in cardiac tissue.J. Electrocardiol.3669–74. 10.1016/j.jelectrocard.2003.09.017
43
van DamP. M.OostendorpT. F.van OosteromA. (2009). Application of the fastest route algorithm in the interactive simulation of the effect of local ischemia on the ECG.Med. Biol. Eng. Comput.4711–20. 10.1007/s11517-008-0391-2
44
WinterJ.HeilJ. E.SchumannC.LinY.SchannwellC. M.MichelU.et al (1998). Effect of implantable cardioverter/defibrillator lead placement in the right ventricle on defibrillation energy requirements. a combined experimental and clinical study1.Eur. J. Cardio Thoracic Surg.14419–425. 10.1016/S1010-7940(98)00215-2
45
WinterfieldJ. R.KentA. R.KarstE.DalalN.MahapatraS.BunchT. J.et al (2018). Impact of ventricular tachycardia ablation on health care utilization.Heart Rhythm15355–362. 10.1016/j.hrthm.2017.10.009
46
WittkampfF. H.HauerR. N.Robles de MedinaE. O. (1989). Control of radiofrequency lesion size by power regulation.Circulation80962–968. 10.1161/01.CIR.80.4.962
47
YamadaT.McElderryH. T.DoppalapudiH.MurakamiY.YoshidaY.YoshidaN.et al (2008). Idiopathic ventricular arrhythmias originating from the aortic root prevalence, electrocardiographic and electrophysiologic characteristics, and results of radiofrequency catheter ablation.J. Am. Coll. Cardiol.52:40. 10.1016/j.jacc.2008.03.040
48
YangT.YuL.JinQ.WuL.HeB. (2018). Localization of origins of premature ventricular contraction by means of convolutional neural network from 12-lead ECG.IEEE Transact. Biomed. Eng.651662–1671. 10.1109/TBME.2017.2756869
49
YokokawaM.KimH. M.DabbaghG. S.SiontisK. C.Lathkar-PradhanS.JongnarangsinK.et al (2019). Targeting noninducible clinical ventricular tachycardias in patients with prior myocardial infarctions based on stored electrograms.Circulat. Arrhythm. Electrophysiol.12:6978.
50
YoshidaK.LiuT. Y.ScottC.HeroA.YokokawaM.GuptaS.et al (2010). The value of defibrillator electrograms for recognition of clinical ventricular tachycardias and for pace mapping of post-infarction ventricular tachycardia.J. Am. Coll. Cardiol.56969–979. 10.1016/j.jacc.2010.04.043
51
YuJ. K.FranceschiW.HuangQ.PashakhanlooF.BoyleP. M.TrayanovaN. A. (2019). A comprehensive, multiscale framework for evaluation of arrhythmias arising from cell therapy in the whole post-myocardial infarcted heart.Sci. Rep.91–16. 10.1038/s41598-019-45684-0
52
ZhouS.AbdelWahabA.SappJ. L.WarrenJ. W.HoráčekB. M. (2019). Localization of ventricular activation origin from the 12-Lead ECG: a comparison of linear regression with non-linear methods of machine learning.Annal. Biomed. Eng.47403–412. 10.1007/s10439-018-02168-y
Summary
Keywords
ventricular tachycardia, implanted devices, electrograms, automated localization, torso modeling, deep learning
Citation
Monaci S, Gillette K, Puyol-Antón E, Rajani R, Plank G, King A and Bishop M (2021) Automated Localization of Focal Ventricular Tachycardia From Simulated Implanted Device Electrograms: A Combined Physics–AI Approach. Front. Physiol. 12:682446. doi: 10.3389/fphys.2021.682446
Received
18 March 2021
Accepted
31 May 2021
Published
01 July 2021
Volume
12 - 2021
Edited by
Rafael Sebastian, University of Valencia, Spain
Reviewed by
Juan Pablo Martínez, University of Zaragoza, Spain; Alessio Gizzi, Campus Bio-Medico University, Italy; Rodrigo Weber Dos Santos, Juiz de Fora Federal University, Brazil
Updates
Copyright
© 2021 Monaci, Gillette, Puyol-Antón, Rajani, Plank, King and Bishop.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sofia Monaci, sofia.monaci@kcl.ac.uk
This article was submitted to Computational Physiology and Medicine, a section of the journal Frontiers in Physiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.