Ablation-induced left atrial mechanical dysfunction recovers in weeks after ablation.

BACKGROUND: The immediate impact of catheter ablation on left atrial mechanical function and the timeline for its recovery in patients undergoing ablation for atrial fibrillation (AF) remain uncertain. The mechanical function response to catheter ablation in patients with different AF types is poorly understood. METHODS: A total of 113 AF patients were included in this retrospective study. Each patient had three magnetic resonance imaging (MRI) studies in sinus rhythm: one pre-ablation, one immediate post-ablation (within 2 days after ablation), and one post-ablation follow-up MRI (≤ 3 months). We used feature tracking in the MRI cine images to determine peak longitudinal atrial strain (PLAS). We evaluated the change in strain from pre-ablation, immediately after ablation to post-ablation follow-up in a short-term study (< 50 days) and a 3-month study (3 months after ablation). RESULTS: The PLAS exhibited a notable reduction immediately after ablation, compared to both pre-ablation levels and those observed in follow-up studies conducted at short-term (11.1 ± 9.0 days) and 3-month (69.6 ± 39.6 days) intervals. However, there was no difference between follow-up and pre-ablation PLAS. The PLAS returned to 95% pre-ablation level within 10 days. Paroxysmal AF patients had significantly higher pre-ablation PLAS than persistent AF patients in pre-ablation MRIs. Both type AF patients had significantly lower immediate post-ablation PLAS compared with pre-ablation and post-ablation PLAS. CONCLUSION: The present study suggested a significant drop in PLAS immediately after ablation. Left atrial mechanical function recovered within 10 days after ablation. The drop in PLAS did not show a substantial difference between paroxysmal and persistent AF patients.

Influence of material parameter variability on the predicted coronary artery biomechanical environment via uncertainty quantification.

Central to the clinical adoption of patient-specific modeling strategies is demonstrating that simulation results are reliable and safe. Indeed, simulation frameworks must be robust to uncertainty in model input(s), and levels of confidence should accompany results. In this study, we applied a coupled uncertainty quantification-finite element (FE) framework to understand the impact of uncertainty in vascular material properties on variability in predicted stresses. Univariate probability distributions were fit to material parameters derived from layer-specific mechanical behavior testing of human coronary tissue. Parameters were assumed to be probabilistically independent, allowing for efficient parameter ensemble sampling. In an idealized coronary artery geometry, a forward FE model for each parameter ensemble was created to predict tissue stresses under physiologic loading. An emulator was constructed within the UncertainSCI software using polynomial chaos techniques, and statistics and sensitivities were directly computed. Results demonstrated that material parameter uncertainty propagates to variability in predicted stresses across the vessel wall, with the largest dispersions in stress within the adventitial layer. Variability in stress was most sensitive to uncertainties in the anisotropic component of the strain energy function. Moreover, unary and binary interactions within the adventitial layer were the main contributors to stress variance, and the leading factor in stress variability was uncertainty in the stress-like material parameter that describes the contribution of the embedded fibers to the overall artery stiffness. Results from a patient-specific coronary model confirmed many of these findings. Collectively, these data highlight the impact of material property variation on uncertainty in predicted artery stresses and present a pipeline to explore and characterize forward model uncertainty in computational biomechanics.

Influence of Material Parameter Variability on the Predicted Coronary Artery Biomechanical Environment via Uncertainty Quantification.

Central to the clinical adoption of patient-specific modeling strategies is demonstrating that simulation results are reliable and safe. Indeed, simulation frameworks must be robust to uncertainty in model input(s), and levels of confidence should accompany results. In this study, we applied a coupled uncertainty quantification-finite element (FE) framework to understand the impact of uncertainty in vascular material properties on variability in predicted stresses. Univariate probability distributions were fit to material parameters derived from layer-specific mechanical behavior testing of human coronary tissue. Parameters were assumed to be probabilistically independent, allowing for efficient parameter ensemble sampling. In an idealized coronary artery geometry, a forward FE model for each parameter ensemble was created to predict tissue stresses under physiologic loading. An emulator was constructed within the UncertainSCI software using polynomial chaos techniques, and statistics and sensitivities were directly computed. Results demonstrated that material parameter uncertainty propagates to variability in predicted stresses across the vessel wall, with the largest dispersions in stress within the adventitial layer. Variability in stress was most sensitive to uncertainties in the anisotropic component of the strain energy function. Moreover, unary and binary interactions within the adventitial layer were the main contributors to stress variance, and the leading factor in stress variability was uncertainty in the stress-like material parameter that describes the contribution of the embedded fibers to the overall artery stiffness. Results from a patient-specific coronary model confirmed many of these findings. Collectively, these data highlight the impact of material property variation on uncertainty in predicted artery stresses and present a pipeline to explore and characterize forward model uncertainty in computational biomechanics.

Uncertainty quantification of the effect of cardiac position variability in the inverse problem of electrocardiographic imaging.

Objective.Electrocardiographic imaging (ECGI) is a functional imaging modality that consists of two related problems, the forward problem of reconstructing body surface electrical signals given cardiac bioelectric activity, and the inverse problem of reconstructing cardiac bioelectric activity given measured body surface signals. ECGI relies on a model for how the heart generates bioelectric signals which is subject to variability in inputs. The study of how uncertainty in model inputs affects the model output is known as uncertainty quantification (UQ). This study establishes develops, and characterizes the application of UQ to ECGI.Approach.We establish two formulations for applying UQ to ECGI: a polynomial chaos expansion (PCE) based parametric UQ formulation (PCE-UQ formulation), and a novel UQ-aware inverse formulation which leverages our previously established 'joint-inverse' formulation (UQ joint-inverse formulation). We apply these to evaluate the effect of uncertainty in the heart position on the ECGI solutions across a range of ECGI datasets.Main results.We demonstrated the ability of our UQ-ECGI formulations to characterize the effect of parameter uncertainty on the ECGI inverse problem. We found that while the PCE-UQ inverse solution provided more complex outputs such as sensitivities and standard deviation, the UQ joint-inverse solution provided a more interpretable output in the form of a single ECGI solution. We find that between these two methods we are able to assess a wide range of effects that heart position variability has on the ECGI solution.Significance.This study, for the first time, characterizes in detail the application of UQ to the ECGI inverse problem. We demonstrated how UQ can provide insight into the behavior of ECGI using variability in cardiac position as a test case. This study lays the groundwork for future development of UQ-ECGI studies, as well as future development of ECGI formulations which are robust to input parameter variability.

Performance of Off-the-Shelf Machine Learning Architectures and Biases in Detection of Low Left Ventricular Ejection Fraction.

Artificial intelligence - machine learning (AI-ML) is a computational technique that has been demonstrated to be able to extract meaningful clinical information from diagnostic data that are not available using either human interpretation or more simple analysis methods. Recent developments have shown that AI-ML approaches applied to ECGs can accurately predict different patient characteristics and pathologies not detectable by expert physician readers. There is an extensive body of literature surrounding the use of AI-ML in other fields, which has given rise to an array of predefined open-source AI-ML architectures which can be translated to new problems in an "off-the-shelf" manner. Applying "off-the-shelf" AI-ML architectures to ECG-based datasets opens the door for rapid development and identification of previously unknown disease biomarkers. Despite the excellent opportunity, the ideal open-source AI-ML architecture for ECG related problems is not known. Furthermore, there has been limited investigation on how and when these AI-ML approaches fail and possible bias or disparities associated with particular network architectures. In this study, we aimed to: (1) determine if open-source, "off-the-shelf" AI-ML architectures could be trained to classify low LVEF from ECGs, (2) assess the accuracy of different AI-ML architectures compared to each other, and (3) to identify which, if any, patient characteristics are associated with poor AI-ML performance.

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.

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.

UncertainSCI: Uncertainty quantification for computational models in biomedicine and bioengineering.

BACKGROUND: Computational biomedical simulations frequently contain parameters that model physical features, material coefficients, and physiological effects, whose values are typically assumed known a priori. Understanding the effect of variability in those assumed values is currently a topic of great interest. A general-purpose software tool that quantifies how variation in these parameters affects model outputs is not broadly available in biomedicine. For this reason, we developed the 'UncertainSCI' uncertainty quantification software suite to facilitate analysis of uncertainty due to parametric variability. METHODS: We developed and distributed a new open-source Python-based software tool, UncertainSCI, which employs advanced parameter sampling techniques to build polynomial chaos (PC) emulators that can be used to predict model outputs for general parameter values. Uncertainty of model outputs is studied by modeling parameters as random variables, and model output statistics and sensitivities are then easily computed from the emulator. Our approaches utilize modern, near-optimal techniques for sampling and PC construction based on weighted Fekete points constructed by subsampling from a suitably randomized candidate set. RESULTS: Concentrating on two test cases-modeling bioelectric potentials in the heart and electric stimulation in the brain-we illustrate the use of UncertainSCI to estimate variability, statistics, and sensitivities associated with multiple parameters in these models. CONCLUSION: UncertainSCI is a powerful yet lightweight tool enabling sophisticated probing of parametric variability and uncertainty in biomedical simulations. Its non-intrusive pipeline allows users to leverage existing software libraries and suites to accurately ascertain parametric uncertainty in a variety of applications.

Improving Generalization by Learning Geometry-Dependent and Physics-Based Reconstruction of Image Sequences.

Deep neural networks have shown promise in image reconstruction tasks, although often on the premise of large amounts of training data. In this paper, we present a new approach to exploit the geometry and physics underlying electrocardiographic imaging (ECGI) to learn efficiently with a relatively small dataset. We first introduce a non-Euclidean encoding-decoding network that allows us to describe the unknown and measurement variables over their respective geometrical domains. We then explicitly model the geometry-dependent physics in between the two domains via a bipartite graph over their graphical embeddings. We applied the resulting network to reconstruct electrical activity on the heart surface from body-surface potentials. In a series of generalization tasks with increasing difficulty, we demonstrated the improved ability of the network to generalize across geometrical changes underlying the data using less than 10% of training data and fewer variations of training geometry in comparison to its Euclidean alternatives. In both simulation and real-data experiments, we further demonstrated its ability to be quickly fine-tuned to new geometry using a modest amount of data.

Effects of Biventricular Pacing Locations on Anti-Tachycardia Pacing Success in a Patient-Specific Model.

Patients with drug-refractory ventricular tachycardia (VT) often undergo implantation of a cardiac defibrillator (ICD). While life-saving, shock from an ICD can be traumatic. To combat the need for defibrillation, ICDs come equipped with low-energy pacing protocols. These anti-tachycardia pacing (ATP) methods are conventionally delivered from a lead inserted at the apex of the right ventricle (RV) with limited success. Recent studies have shown the promise of biventricular leads placed in the left ventricle (LV) for ATP delivery. This study tested the hypothesis that stimulating ATP from multiple biventricular locations will improve termination rates in a patient-specific computational model. VT was first induced in the model, followed by ATP delivery from 1-4 biventricular stimulus sites. We found that combining stimulation sites does not alter termination success so long as a critical stimulus site is included. Combining the RV stimulus site with any combination of LV sites did not affect ATP success except for one case. Including the RV site may allow biventricular ATP to be a robust approach across different scar distributions without affecting the efficacy of other stimulation sites. Combining sites may increase the likelihood of including a critical stimulus site when such information cannot be ascertained.

Transfer Learning for Improved Classification of Drivers in Atrial Fibrillation.

"Drivers" are theorized mechanisms for persistent atrial fibrillation. Machine learning algorithms have been used to identify drivers, but the small size of current driver datasets limits their performance. We hypothesized that pretraining with unsupervised learning on a large dataset of unlabeled electrograms would improve classifier accuracy on a smaller driver dataset. In this study, we used a SimCLR-based framework to pretrain a residual neural network on a dataset of 113K unlabeled 64-electrode measurements and found weighted testing accuracy to improve over a non-pretrained network (78.6±3.9% vs 71.9±3.3%). This lays ground for development of superior driver detection algorithms and supports use of transfer learning for other datasets of endocardial electrograms.

All Roads Lead to Rome: Diverse Etiologies of Tricuspid Regurgitation Create a Predictable Constellation of Right Ventricular Shape Changes.

Introduction: Myriad disorders cause right ventricular (RV) dilation and lead to tricuspid regurgitation (TR). Because the thin-walled, flexible RV is mechanically coupled to the pulmonary circulation and the left ventricular septum, it distorts with any disturbance in the cardiopulmonary system. TR, therefore, can result from pulmonary hypertension, left heart failure, or intrinsic RV dysfunction; but once it occurs, TR initiates a cycle of worsening RV volume overload, potentially progressing to right heart failure. Characteristic three-dimensional RV shape-changes from this process, and changes particular to individual TR causes, have not been defined in detail. Methods: Cardiac MRI was obtained in 6 healthy volunteers, 41 patients with ≥ moderate TR, and 31 control patients with cardiac disease without TR. The mean shape of each group was constructed using a three-dimensional statistical shape model via the particle-based shape modeling approach. Changes in shape were examined across pulmonary hypertension and congestive heart failure subgroups using principal component analysis (PCA). A logistic regression approach based on these PCA modes identified patients with TR using RV shape alone. Results: Mean RV shape in patients with TR exhibited free wall bulging, narrowing of the base, and blunting of the RV apex compared to controls (p < 0.05). Using four primary PCA modes, a logistic regression algorithm identified patients with TR correctly with 82% recall and 87% precision. In patients with pulmonary hypertension without TR, RV shape was narrower and more streamlined than in healthy volunteers. However, in RVs with TR and pulmonary hypertension, overall RV shape continued to demonstrate the free wall bulging characteristic of TR. In the subgroup of patients with congestive heart failure without TR, this intermediate state of RV muscular hypertrophy was not present. Conclusion: The multiple causes of TR examined in this study changed RV shape in similar ways. Logistic regression classification based on these shape changes reliably identified patients with TR regardless of etiology. Furthermore, pulmonary hypertension without TR had unique shape features, described here as the "well compensated" RV. These results suggest shape modeling as a promising tool for defining severity of RV disease and risk of decompensation, particularly in patients with pulmonary hypertension.

Reconstruction of cardiac position using body surface potentials.

Electrocardiographic imaging (ECGI) is a noninvasive technique to assess the bioelectric activity of the heart which has been applied to aid in clinical diagnosis and management of cardiac dysfunction. ECGI is built on mathematical models that take into account several patient specific factors including the position of the heart within the torso. Errors in the localization of the heart within the torso, as might arise due to natural changes in heart position from respiration or changes in body position, contribute to errors in ECGI reconstructions of the cardiac activity, thereby reducing the clinical utility of ECGI. In this study we present a novel method for the reconstruction of cardiac geometry utilizing noninvasively acquired body surface potential measurements. Our geometric correction method simultaneously estimates the cardiac position over a series of heartbeats by leveraging an iterative approach which alternates between estimating the cardiac bioelectric source across all heartbeats and then estimating cardiac positions for each heartbeat. We demonstrate that our geometric correction method is able to reduce geometric error and improve ECGI accuracy in a wide range of testing scenarios. We examine the performance of our geometric correction method using different activation sequences, ranges of cardiac motion, and body surface electrode configurations. We find that after geometric correction resulting ECGI solution accuracy is improved and variability of the ECGI solutions between heartbeats is substantially reduced.

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

OBJECTIVE: To investigatecardiac 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 anegligible 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.

Combining endocardial mapping and electrocardiographic imaging (ECGI) for improving PVC localization: A feasibility study.

INTRODUCTION: Accurate reconstruction of cardiac activation wavefronts is crucial for clinical diagnosis, management, and treatment of cardiac arrhythmias. Furthermore, reconstruction of activation profiles within the intramural myocardium has long been impossible because electrical mapping was only performed on the endocardial surface. Recent advancements in electrocardiographic imaging (ECGI) have made endocardial and epicardial activation mapping possible. We propose a novel approach to use both endocardial and epicardial mapping in a combined approach to reconstruct intramural activation times. OBJECTIVE: To implement and validate a combined epicardial/endocardial intramural activation time reconstruction technique. METHODS: We used 11 simulations of ventricular activation paced from sites throughout myocardial wall and extracted endocardial and epicardial activation maps at approximate clinical resolution. From these maps, we interpolated the activation times through the myocardium using thin-plate-spline radial basis functions. We evaluated activation time reconstruction accuracy using root-mean-squared error (RMSE) of activation times and the percent of nodes within 1 ms of the ground truth. RESULTS: Reconstructed intramural activation times showed an RMSE and percentage of nodes within 1 ms of the ground truth simulations of 3 ms and 70%, respectively. In the worst case, the RMSE and percentage of nodes were 4 ms and 60%, respectively. CONCLUSION: We showed that a simple, yet effective combination of clinical endocardial and epicardial activation maps can accurately reconstruct intramural wavefronts. Furthermore, we showed that this approach provided robust reconstructions across multiple intramural stimulation sites.

Transient recovery of epicardial and torso ST-segment ischemic signals during cardiac stress tests: A possible physiological mechanism.

BACKGROUND: Acute myocardial ischemia has several characteristic ECG findings, including clinically detectable ST-segment deviations. However, the sensitivity and specificity of diagnosis based on ST-segment changes are low. Furthermore, ST-segment deviations have been shown to be transient and spontaneously recover without any indication the ischemic event has subsided. OBJECTIVE: Assess the transient recovery of ST-segment deviations on remote recording electrodes during a partial occlusion cardiac stress test and compare them to intramyocardial ST-segment deviations. METHODS: We used a previously validated porcine experimental model of acute myocardial ischemia with controllable ischemic load and simultaneous electrical measurements within the heart wall, on the epicardial surface, and on the torso surface. Simulated cardiac stress tests were induced by occluding a coronary artery while simultaneously pacing rapidly or infusing dobutamine to stimulate cardiac function. Postexperimental imaging created anatomical models for data visualization and quantification. Markers of ischemia were identified as deviations in the potentials measured at 40% of the ST-segment. Intramural cardiac conduction speed was also determined using the inverse gradient method. We assessed changes in intramyocardial ischemic volume proportion, conduction speed, clinical presence of ischemia on remote recording arrays, and regional changes to intramyocardial ischemia. We defined the peak deviation response time as the time interval after onset of ischemia at which maximum ST-segment deviation was achieved, and ST-recovery time was the interval when ST deviation returned to below thresholded of ST elevation. RESULTS: In both epicardial and torso recordings, the peak ST-segment deviation response time was 4.9±1.1 min and the ST-recovery time was approximately 7.9±2.5 min, both well before the termination of the ischemic stress. At peak response time, conduction speed was reduced by 50% and returned to near baseline at ST-recovery. The overall ischemic volume proportion initially increased, on average, to 37% at peak response time; however, it recovered to only 30% at the ST-recovery time. By contrast, the subepicardial region of the myocardial wall showed 40% ischemic volume at peak response time and recovered much more strongly to 25% as epicardial ST-segment deviations returned to baseline. CONCLUSION: Our data show that remote ischemic signal recovery correlates with a recovery of the subepicardial myocardium, whereas subendocardial ischemic development persists.

Pharmacological and simulated exercise cardiac stress tests produce different ischemic signatures in high-resolution experimental mapping studies.

OBJECTIVE: Test the hypothesis that exercise and pharmacological cardiac stressors create different electrical ischemic signatures. INTRODUCTION: Current clinical stress tests for detecting ischemia lack sensitivity and specificity. One unexplored source of the poor detection is whether pharmacological stimulation and regulated exercise produce identical cardiac stress. METHODS: We used a porcine model of acute myocardial ischemia in which animals were instrumented with transmural plunge-needle electrodes, an epicardial sock array, and torso arrays to simultaneously measure cardiac electrical signals within the heart wall, the epicardial surface, and the torso surface, respectively. Ischemic stress via simulated exercise and pharmacological stimulation were created with rapid electrical pacing and dobutamine infusion, respectively, and mimicked clinical stress tests of five 3-minute stages. Perfusion to the myocardium was regulated by a hydraulic occluder around the left anterior descending coronary artery. Ischemia was measured as deflections to the ST-segment on ECGs and electrograms. RESULTS: Across eight experiments with 30 (14 simulated exercise and 16 dobutamine) ischemic interventions, the spatial correlations between exercise and pharmacological stress diverged at stage three or four during interventions (p<0.05). We found more detectable ST-segment changes on the epicardial surface during simulated exercise than with dobutamine (p<0.05). The intramyocardial ischemia formed during simulated exercise had larger ST40 potential gradient magnitudes (p<0.05). CONCLUSION: We found significant differences on the epicardium between cardiac stress types using our experimental model, which became more pronounced at the end stages of each test. A possible mechanism for these differences was the larger ST40 potential gradient magnitudes within the myocardium during exercise. The presence of microvascular dysfunction during exercise and its absence during dobutamine stress may explain these differences.

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

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

The effect of interpolating low amplitude leads on the inverse reconstruction of cardiac electrical activity.

Electrocardiographic imaging is an imaging modality that has been introduced recently to help in visualizing the electrical activity of the heart and consequently guide the ablation therapy for ventricular arrhythmias. One of the main challenges of this modality is that the electrocardiographic signals recorded at the torso surface are contaminated with noise from different sources. Low amplitude leads are more affected by noise due to their low peak-to-peak amplitude. In this paper, we have studied 6 datasets from two torso tank experiments (Bordeaux and Utah experiments) to investigate the impact of removing or interpolating these low amplitude leads on the inverse reconstruction of cardiac electrical activity. Body surface potential maps used were calculated by using the full set of recorded leads, removing 1, 6, 11, 16, or 21 low amplitude leads, or interpolating 1, 6, 11, 16, or 21 low amplitude leads using one of the three interpolation methods - Laplacian interpolation, hybrid interpolation, or the inverse-forward interpolation. The epicardial potential maps and activation time maps were computed from these body surface potential maps and compared with those recorded directly from the heart surface in the torso tank experiments. There was no significant change in the potential maps and activation time maps after the removal of up to 11 low amplitude leads. Laplacian interpolation and hybrid interpolation improved the inverse reconstruction in some datasets and worsened it in the rest. The inverse forward interpolation of low amplitude leads improved it in two out of 6 datasets and at least remained the same in the other datasets. It was noticed that after doing the inverse-forward interpolation, the selected lambda value was closer to the optimum lambda value that gives the inverse solution best correlated with the recorded one.

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.