Regional conduction velocities determined by noninvasive mapping are associated with arrhythmia-free survival after atrial fibrillation ablation.

BACKGROUND: Atrial arrhythmogenic substrate is a key determinant of atrial fibrillation (AF) recurrence after pulmonary vein isolation (PVI), and reduced conduction velocities have been linked to adverse outcome. However, a noninvasive method to assess such electrophysiologic substrate is not available to date. OBJECTIVE: This study aimed to noninvasively assess regional conduction velocities and their association with arrhythmia-free survival after PVI. METHODS: A consecutive 52 patients scheduled for AF ablation (PVI only) and 19 healthy controls were prospectively included and received electrocardiographic imaging (ECGi) to noninvasively determine regional atrial conduction velocities in sinus rhythm. A novel ECGi technology obviating the need of additional computed tomography or cardiac magnetic resonance imaging was applied and validated by invasive mapping. RESULTS: Mean ECGi-determined atrial conduction velocities were significantly lower in AF patients than in healthy controls (1.45 ± 0.15 m/s vs 1.64 ± 0.15 m/s; P < .0001). Differences were particularly pronounced in a regional analysis considering only the segment with the lowest average conduction velocity in each patient (0.8 ± 0.22 m/s vs 1.08 ± 0.26 m/s; P < .0001). This average conduction velocity of the "slowest" segment was independently associated with arrhythmia recurrence and better discriminated between PVI responders and nonresponders than previously proposed predictors, including left atrial size and late gadolinium enhancement (magnetic resonance imaging). Patients without slow-conduction areas (mean conduction velocity <0.78 m/s) showed significantly higher 12-month arrhythmia-free survival than those with 1 or more slow-conduction areas (88.9% vs 48.0%; P = .002). CONCLUSION: This is the first study to investigate regional atrial conduction velocities noninvasively. The absence of ECGi-determined slow-conduction areas well discriminates PVI responders from nonresponders. Such noninvasive assessment of electrical arrhythmogenic substrate may guide treatment strategies and be a step toward personalized AF therapy.

Non-invasive estimation of atrial fibrillation driver position using long-short term memory neural networks and body surface potentials.

BACKGROUND AND OBJECTIVE: Atrial Fibrillation (AF) is a supraventricular tachyarrhythmia that can lead to thromboembolism, hearlt failure, ischemic stroke, and a decreased quality of life. Characterizing the locations where the mechanisms of AF are initialized and maintained is key to accomplishing an effective ablation of the targets, hence restoring sinus rhythm. Many methods have been investigated to locate such targets in a non-invasive way, such as Electrocardiographic Imaging, which enables an on-invasive and panoramic characterization of cardiac electrical activity using recording Body Surface Potentials (BSP) and a torso model of the patient. Nonetheless, this technique entails some major issues stemming from solving the inverse problem, which is known to be severely ill-posed. In this context, many machine learning and deep learning approaches aim to tackle the characterization and classification of AF targets to improve AF diagnosis and treatment. METHODS: In this work, we propose a method to locate AF drivers as a supervised classification problem. We employed a hybrid form of the convolutional-recurrent network which enables feature extraction and sequential data modeling utilizing labeled realistic computerized AF models. Thus, we used 16 AF electrograms, 1 atrium, and 10 torso geometries to compute the forward problem. Previously, the AF models were labeled by assigning each sample of the signals a region from the atria from 0 (no driver) to 7, according to the spatial location of the AF driver. The resulting 160 BSP signals, which resemble a 64-lead vest recording, are preprocessed and then introduced into the network following a 4-fold cross-validation in batches of 50 samples. RESULTS: The results show a mean accuracy of 74.75% among the 4 folds, with a better performance in detecting sinus rhythm, and drivers near the left superior pulmonary vein (R1), and right superior pulmonary vein (R3) whose mean sensitivity bounds around 84%-87%. Significantly good results are obtained in mean sensitivity (87%) and specificity (83%) in R1. CONCLUSIONS: Good results in R1 are highly convenient since AF drivers are commonly found in this area: the left atrial appendage, as suggested in some previous studies. These promising results indicate that using CNN-LSTM networks could lead to new strategies exploiting temporal correlations to address this challenge effectively.

Standardized 2D atrial mapping and its clinical applications.

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

AF driver detection in pulmonary vein area by electropcardiographic imaging: Relation with a favorable outcome of pulmonary vein isolation.

Pulmonary vein isolation (PVI) is the most successful treatment for atrial fibrillation (AF) nowadays. However, not all AF patients benefit from PVI. In this study, we evaluate the use of ECGI to identify reentries and relate rotor density in the pulmonary vein (PV) area as an indicator of PVI outcome. Rotor maps were computed in a set of 29 AF patients using a new rotor detection algorithm. The relationship between the distribution of reentrant activity and the clinical outcome after PVI was studied. The number of rotors and proportion of PSs in different atrial regions were computed and compared retrospectively in two groups of patients: patients that remained in sinus rhythm 6 months after PVI and patients with arrhythmia recurrence. The total number of rotors obtained was higher in patients returning to arrhythmia after the ablation (4.31 ± 2.77 vs. 3.58 ± 2.67%, p = 0.018). However, a significantly higher concentration of PSs in the pulmonary veins was found in patients that remained in sinus rhythm (10.20 ± 12.40% vs. 5.19 ± 9.13%, p = 0.011) 6 months after PVI. The results obtained show a direct relationship between the expected AF mechanism and the electrophysiological parameters provided by ECGI, suggesting that this technology offers relevant information to predict the clinical outcome after PVI in AF patients.

Electrocardiographic imaging in the atria.

The inverse problem of electrocardiography or electrocardiographic imaging (ECGI) is a technique for reconstructing electrical information about cardiac surfaces from noninvasive or non-contact recordings. ECGI has been used to characterize atrial and ventricular arrhythmias. Although it is a technology with years of progress, its development to characterize atrial arrhythmias is challenging. Complications can arise when trying to describe the atrial mechanisms that lead to abnormal propagation patterns, premature or tachycardic beats, and reentrant arrhythmias. This review addresses the various ECGI methodologies, regularization methods, and post-processing techniques used in the atria, as well as the context in which they are used. The current advantages and limitations of ECGI in the fields of research and clinical diagnosis of atrial arrhythmias are outlined. In addition, areas where ECGI efforts should be concentrated to address the associated unsatisfied needs from the atrial perspective are discussed.

Effects of torso mesh density and electrode distribution on the accuracy of electrocardiographic imaging during atrial fibrillation.

Introduction: Electrocardiographic Imaging (ECGI) allows computing the electrical activity in the heart non-invasively using geometrical information of the patient and multiple body surface signals. In the present study we investigate the influence of the number of nodes of geometrical meshes and recording ECG electrodes distribution to compute ECGI during atrial fibrillation (AF). Methods: Torso meshes from 100 to 2000 nodes heterogeneously and homogeneously distributed were compared. Signals from nine AF realistic mathematical simulations were used for computing the ECGI. Results for each torso mesh were compared with the ECGI computed with a 4,000 nodes reference torso. In addition, real AF recordings from 25 AF patients were used to compute ECGI in torso meshes from 100 to 1,000 nodes. Results were compared with a reference torso of 2000 nodes. Torsos were remeshed either by reducing the number of nodes while maximizing the overall shape preservation and then assigning the location of the electrodes as the closest node in the new mesh or by forcing the remesher to place a node at each electrode location. Correlation coefficients, relative difference measurements and relative difference of dominant frequencies were computed to evaluate the impact on signal morphology of each torso mesh. Results: For remeshed torsos where electrodes match with a geometrical node in the mesh, all mesh densities presented similar results. On the other hand, in torsos with electrodes assigned to closest nodes in remeshed geometries performance metrics were dependent on mesh densities, with correlation coefficients ranging from 0.53 ± 0.06 to 0.92 ± 0.04 in simulations or from 0.42 ± 0.38 to 0.89 ± 0.2 in patients. Dominant frequency relative errors showed the same trend with values from 1.14 ± 0.26 to 0.55 ± 0.21 Hz in simulations and from 0.91 ± 0.56 to 0.45 ± 0.41 Hz in patients. Conclusion: The effect of mesh density in ECGI is minimal when the location of the electrode is preserved as a node in the mesh. Torso meshes constructed without imposing electrodes to constitute nodes in the torso geometry should contain at least 400 nodes homogeneously distributed so that a distance between nodes is below 4 cm.

Non-invasive Estimation of Atrial Fibrillation Driver Position With Convolutional Neural Networks and Body Surface Potentials.

Atrial fibrillation (AF) is characterized by complex and irregular propagation patterns, and AF onset locations and drivers responsible for its perpetuation are the main targets for ablation procedures. ECG imaging (ECGI) has been demonstrated as a promising tool to identify AF drivers and guide ablation procedures, being able to reconstruct the electrophysiological activity on the heart surface by using a non-invasive recording of body surface potentials (BSP). However, the inverse problem of ECGI is ill-posed, and it requires accurate mathematical modeling of both atria and torso, mainly from CT or MR images. Several deep learning-based methods have been proposed to detect AF, but most of the AF-based studies do not include the estimation of ablation targets. In this study, we propose to model the location of AF drivers from BSP as a supervised classification problem using convolutional neural networks (CNN). Accuracy in the test set ranged between 0.75 (SNR = 5 dB) and 0.93 (SNR = 20 dB upward) when assuming time independence, but it worsened to 0.52 or lower when dividing AF models into blocks. Therefore, CNN could be a robust method that could help to non-invasively identify target regions for ablation in AF by using body surface potential mapping, avoiding the use of ECGI.

Analysis of the response of human iPSC-derived cardiomyocyte tissue to ICaL block. A combined in vitro and in silico approach.

The high incidence of cardiac arrythmias underlines the need for the assessment of pharmacological therapies. In this field of drug efficacy, as in the field of drug safety highlighted by the Comprehensive in Vitro Proarrhythmia Assay initiative, new pillars for research have become crucial: firstly, the integration of in-silico experiments, and secondly the evaluation of fully integrated biological systems, such as human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs). In this study, we therefore aimed to combine in-vitro experiments and in-silico simulations to evaluate the antiarrhythmic effect of L-type calcium current (ICaL) block in hiPSC-CMs. For this, hiPSC-CM preparations were cultured and an equivalent virtual tissue was modeled. Re-entry patterns of electrical activation were induced and several biomarkers were obtained before and after ICaL block. The virtual hiPSC-CM simulations were also reproduced using a tissue composed of adult ventricular cardiomyocytes (hAdultV-CMs). The analysis of phases, currents and safety factor for propagation showed an increased size of the re-entry core when ICaL was blocked as a result of depressed cellular excitability. The bigger wavefront curvature yielded reductions of 12.2%, 6.9%, and 4.2% in the frequency of the re-entry for hiPSC-CM cultures, virtual hiPSC-CM, and hAdultV-CM tissues, respectively. Furthermore, ICaL block led to a 47.8% shortening of the vulnerable window for re-entry in the virtual hiPSC-CM tissue and to re-entry vanishment in hAdultV-CM tissue. The consistent behavior between in-vitro and in-silico hiPSC-CMs and between in-silico hiPSC-CMs and hAdultV-CMs evidences that virtual hiPSC-CM tissues are suitable for assessing cardiac efficacy, as done in the present study through the analysis of ICaL block.

Ranolazine-Mediated Attenuation of Mechanoelectric Feedback in Atrial Myocyte Monolayers.

BACKGROUND: Mechanical stretch increases Na+ inflow into myocytes, related to mechanisms including stretch-activated channels or Na+/H+ exchanger activation, involving Ca2+ increase that leads to changes in electrophysiological properties favoring arrhythmia induction. Ranolazine is an antianginal drug with confirmed beneficial effects against cardiac arrhythmias associated with the augmentation of I NaL current and Ca2+ overload. OBJECTIVE: This study investigates the effects of mechanical stretch on activation patterns in atrial cell monolayers and its pharmacological response to ranolazine. METHODS: Confluent HL-1 cells were cultured in silicone membrane plates and were stretched to 110% of original length. The characteristics of in vitro fibrillation (dominant frequency, regularity index, density of phase singularities, rotor meandering, and rotor curvature) were analyzed using optical mapping in order to study the mechanoelectric response to stretch under control conditions and ranolazine action. RESULTS: HL-1 cell stretch increased fibrillatory dominant frequency (3.65 ± 0.69 vs. 4.35 ± 0.74 Hz, p < 0.01) and activation complexity (1.97 ± 0.45 vs. 2.66 ± 0.58 PS/cm2, p < 0.01) under control conditions. These effects were related to stretch-induced changes affecting the reentrant patterns, comprising a decrease in rotor meandering (0.72 ± 0.12 vs. 0.62 ± 0.12 cm/s, p < 0.001) and an increase in wavefront curvature (4.90 ± 0.42 vs. 5.68 ± 0.40 rad/cm, p < 0.001). Ranolazine reduced stretch-induced effects, attenuating the activation rate increment (12.8% vs. 19.7%, p < 0.01) and maintaining activation complexity-both parameters being lower during stretch than under control conditions. Moreover, under baseline conditions, ranolazine slowed and regularized the activation patterns (3.04 ± 0.61 vs. 3.65 ± 0.69 Hz, p < 0.01). CONCLUSION: Ranolazine attenuates the modifications of activation patterns induced by mechanical stretch in atrial myocyte monolayers.

Automatic quality electrogram assessment improves phase-based reentrant activity identification in atrial fibrillation.

Identification of reentrant activity driving atrial fibrillation (AF) is increasingly important to ablative therapies. The goal of this work is to study how the automatically-classified quality of the electrograms (EGMs) affects reentrant AF driver localization. EGMs from 259 AF episodes obtained from 29 AF patients were recorded using 64-poles basket catheters and were manually classified according to their quality. An algorithm capable of identifying signal quality was developed using time and spectral domain parameters. Electrical reentries were identified in 3D phase maps using phase transform and were compared with those obtained with a 2D activation-based method. Effect of EGM quality was studied by discarding 3D phase reentries detected in regions with low-quality EGMs. Removal of reentries identified by 3D phase analysis in regions with low-quality EGMs improved its performance, increasing the area under the ROC curve (AUC) from 0.69 to 0.80. The EGMs quality classification algorithm showed an accurate performance for EGM classification (AUC 0.94) and reentry detection (AUC 0.80). Automatic classification of EGM quality based on time and spectral signal parameters is feasible and accurate, avoiding the manual labelling. Discard of reentries identified in regions with automatically-detected poor-quality EGMs improved the specificity of the 3D phase-based method for AF driver identification.

Noninvasive Assessment of Complexity of Atrial Fibrillation: Correlation With Contact Mapping and Impact of Ablation.

BACKGROUND: It is difficult to noninvasively phenotype atrial fibrillation (AF) in a way that reflects clinical end points such as response to therapy. We set out to map electrical patterns of disorganization and regions of reentrant activity in AF from the body surface using electrocardiographic imaging, calibrated to panoramic intracardiac recordings and referenced to AF termination by ablation. METHODS: Bi-atrial intracardiac electrograms of 47 patients with AF at ablation (30 persistent, 29 male, 63±9 years) were recorded with 64-pole basket catheters and simultaneous 57-lead body surface ECGs. Atrial epicardial electrical activity was reconstructed and organized sites were invasively and noninvasively tracked in 3-dimension using phase singularity. In a subset of 17 patients, sites of AF organization were targeted for ablation. RESULTS: Body surface mapping showed greater AF organization near intracardially detected drivers than elsewhere, both in phase singularity density (2.3±2.1 versus 1.9±1.6; P=0.02) and number of drivers (3.2±2.3 versus 2.7±1.7; P=0.02). Complexity, defined as the number of stable AF reentrant sites, was concordant between noninvasive and invasive methods (r2=0.5; CC=0.71). In the subset receiving targeted ablation, AF complexity showed lower values in those in whom AF terminated than those in whom AF did not terminate (P<0.01). CONCLUSIONS: AF complexity tracked noninvasively correlates well with organized and disorganized regions detected by panoramic intracardiac mapping and correlates with the acute outcome by ablation. This approach may assist in bedside monitoring of therapy or in improving the efficacy of ongoing ablation procedures.

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

Introduction: Regional differences in activation rates may contribute to the electrical substrates that maintain atrial fibrillation (AF), and estimating them non-invasively may help guide ablation or select anti-arrhythmic medications. We tested whether non-invasive assessment of regional AF rate accurately represents intracardiac recordings. Methods: In 47 patients with AF (27 persistent, age 63 ± 13 years) we performed 57-lead non-invasive Electrocardiographic Imaging (ECGI) in AF, simultaneously with 64-pole intracardiac signals of both atria. ECGI was reconstructed by Tikhonov regularization. We constructed personalized 3D AF rate distribution maps by Dominant Frequency (DF) analysis from intracardiac and non-invasive recordings. Results: Raw intracardiac and non-invasive DF differed substantially, by 0.54 Hz [0.13 - 1.37] across bi-atrial regions (R 2 = 0.11). Filtering by high spectral organization reduced this difference to 0.10 Hz (cycle length difference of 1 - 11 ms) [0.03 - 0.42] for patient-level comparisons (R 2 = 0.62), and 0.19 Hz [0.03 - 0.59] and 0.20 Hz [0.04 - 0.61] for median and highest DF, respectively. Non-invasive and highest DF predicted acute ablation success (p = 0.04). Conclusion: Non-invasive estimation of atrial activation rates is feasible and, when filtered by high spectral organization, provide a moderate estimate of intracardiac recording rates in AF. Non-invasive technology could be an effective tool to identify patients who may respond to AF ablation for personalized therapy.

Optical imaging of voltage and calcium in isolated hearts: Linking spatiotemporal heterogeneities and ventricular fibrillation initiation.

BACKGROUND: Alternans have been associated with the development of ventricular fibrillation and its control has been proposed as antiarrhythmic strategy. However, cardiac arrhythmias are a spatiotemporal phenomenon in which multiple factors are involved (e.g. calcium and voltage spatial alternans or heterogeneous conduction velocity) and how an antiarrhythmic drug modifies these factors is poorly understood. OBJECTIVE: The objective of the present study is to evaluate the relation between spatial electrophysiological properties (i.e. spatial discordant alternans and conduction velocity) and the induction of ventricular fibrillation (VF) when a calcium blocker is applied. METHODS: The mechanisms of initiation of VF were studied by simultaneous epicardial voltage and calcium optical mapping in isolated rabbit hearts using an incremental fast pacing protocol. The additional value of analyzing spatial phenomena in the generation of unidirectional blocks and reentries as precursors of VF was depicted. Specifically, the role of action potential duration (APD), calcium transients (CaT), spatial alternans and conduction velocity in the initiation of VF was evaluated during basal conditions and after the administration of verapamil. RESULTS: Our results enhance the relation between (1) calcium spatial alternans and (2) slow conduction velocities with the dynamic creation of unidirectional blocks that allowed the induction of VF. In fact, the administration of verapamil demonstrated that calcium and not voltage spatial alternans were the main responsible for VF induction. CONCLUSIONS: VF induction at high activation rates was linked with the concurrence of a low conduction velocity and high magnitude of calcium alternans, but not necessarily related with increases of APD. Verapamil can postpone the development of cardiac alternans and the apparition of ventricular arrhythmias.

Clinical Characteristics and Electrophysiological Mechanisms Underlying Brugada ECG in Patients With Severe Hyperkalemia.

Background Several metabolic conditions can cause the Brugada ECG pattern, also called Brugada phenotype (BrPh). We aimed to define the clinical characteristics and outcome of BrPh patients and elucidate the mechanisms underlying BrPh attributed to hyperkalemia. Methods and Results We prospectively identified patients hospitalized with severe hyperkalemia and ECG diagnosis of BrPh and compared their clinical characteristics and outcome with patients with hyperkalemia but no BrPh ECG. Computer simulations investigated the roles of extracellular potassium increase, fibrosis at the right ventricular outflow tract, and epicardial/endocardial gradients in transient outward current. Over a 6-year period, 15 patients presented severe hyperkalemia with BrPh ECG that was transient and disappeared after normalization of their serum potassium. Most patients were admitted because of various severe medical conditions causing hyperkalemia. Six (40%) patients presented malignant arrhythmias and 6 died during admission. Multiple logistic regression analysis revealed that higher serum potassium levels (odds ratio, 15.8; 95% CI, 3.1-79; P=0.001) and male sex (odds ratio, 17; 95% CI, 1.05-286; P=0.045) were risk factors for developing BrPh ECG in patients with severe hyperkalemia. In simulations, hyperkalemia yielded BrPh by promoting delayed and heterogeneous right ventricular outflow tract activation attributed to elevation of resting potential, reduced availability of inward sodium channel conductance, and increased right ventricular outflow tract fibrosis. An elevated transient outward current gradient contributed to, but was not essential for, the BrPh phenotype. Conclusions In patients with severe hyperkalemia, a BrPh ECG is associated with malignant arrhythmias and all-cause mortality secondary to resting potential depolarization, reduced sodium current availability, and fibrosis at the right ventricular outflow tract.

Solving Inaccuracies in Anatomical Models for Electrocardiographic Inverse Problem Resolution by Maximizing Reconstruction Quality.

Electrocardiographic Imaging has become an increasingly used technique for non-invasive diagnosis of cardiac arrhythmias, although the need for medical imaging technology to determine the anatomy hinders its introduction in the clinical practice. This paper explores the ability of a new metric based on the inverse reconstruction quality for the location and orientation of the atrial surface inside the torso. Body surface electrical signals from 31 realistic mathematical models and four AF patients were used to estimate the optimal position of the atria inside the torso. The curvature of the L-curve from the Tikhonov method, which was found to be related to the inverse reconstruction quality, was measured after application of deviations in atrial position and orientation. Independent deviations in the atrial position were solved by finding the maximal L-curve curvature with an error of 1.7 ± 2.4 mm in mathematical models and 9.1 ± 11.5 mm in patients. For the case of independent angular deviations, the error in location by using the L-curve was 5.8±7.1° in mathematical models and 12.4° ± 13.2° in patients. The ability of the L-curve curvature was tested also under superimposed uncertainties in the three axis of translation and in the three axis of rotation, and the error in location was of 2.3 ± 3.2 mm and 6.4° ± 7.1° in mathematical models, and 7.9±10.7 mm and 12.1°±15.5° in patients. The curvature of L-curve is a useful marker for the atrial position and would allow emending the inaccuracies in its location.

Corrigendum: Regularization Techniques for ECG Imaging during Atrial Fibrillation: A Computational Study.

[This corrects the article on p. 466 in vol. 7, PMID: 27790158.].

Regularization Techniques for ECG Imaging during Atrial Fibrillation: A Computational Study.

The inverse problem of electrocardiography is usually analyzed during stationary rhythms. However, the performance of the regularization methods under fibrillatory conditions has not been fully studied. In this work, we assessed different regularization techniques during atrial fibrillation (AF) for estimating four target parameters, namely, epicardial potentials, dominant frequency (DF), phase maps, and singularity point (SP) location. We use a realistic mathematical model of atria and torso anatomy with three different electrical activity patterns (i.e., sinus rhythm, simple AF, and complex AF). Body surface potentials (BSP) were simulated using Boundary Element Method and corrupted with white Gaussian noise of different powers. Noisy BSPs were used to obtain the epicardial potentials on the atrial surface, using 14 different regularization techniques. DF, phase maps, and SP location were computed from estimated epicardial potentials. Inverse solutions were evaluated using a set of performance metrics adapted to each clinical target. For the case of SP location, an assessment methodology based on the spatial mass function of the SP location, and four spatial error metrics was proposed. The role of the regularization parameter for Tikhonov-based methods, and the effect of noise level and imperfections in the knowledge of the transfer matrix were also addressed. Results showed that the Bayes maximum-a-posteriori method clearly outperforms the rest of the techniques but requires a priori information about the epicardial potentials. Among the purely non-invasive techniques, Tikhonov-based methods performed as well as more complex techniques in realistic fibrillatory conditions, with a slight gain between 0.02 and 0.2 in terms of the correlation coefficient. Also, the use of a constant regularization parameter may be advisable since the performance was similar to that obtained with a variable parameter (indeed there was no difference for the zero-order Tikhonov method in complex fibrillatory conditions). Regarding the different targets, DF and SP location estimation were more robust with respect to pattern complexity and noise, and most algorithms provided a reasonable estimation of these parameters, even when the epicardial potentials estimation was inaccurate. Finally, the proposed evaluation procedure and metrics represent a suitable framework for techniques benchmarking and provide useful insights for the clinical practice.

Balance between sodium and calcium currents underlying chronic atrial fibrillation termination: An in silico intersubject variability study.

BACKGROUND: Atrial remodeling as a result of long-standing persistent atrial fibrillation (AF) induces substrate modifications that lead to different perpetuation mechanisms than in paroxysmal AF and a reduction in the efficacy of antiarrhythmic treatments. OBJECTIVE: The purpose of this study was to identify the ionic current modifications that could destabilize reentries during chronic AF and serve to personalize antiarrhythmic strategies. METHODS: A population of 173 mathematical models of remodeled human atrial tissue with realistic intersubject variability was developed based on action potential recordings of 149 patients diagnosed with AF. The relationship of each ionic current with AF maintenance and the dynamics of functional reentries (rotor meandering, dominant frequency) were evaluated by means of 3-dimensional simulations. RESULTS: Self-sustained reentries were maintained in 126 (73%) of the simulations. AF perpetuation was associated with higher expressions of INa and ICaL (P <.01), with no significant differences in the remaining currents. ICaL blockade promoted AF extinction in 30% of these 126 models. The mechanism of AF termination was related with collisions between rotors because of an increase in rotor meandering (1.71 ± 2.01cm2) and presented an increased efficacy in models with a depressed INa (P <.01). CONCLUSION: Mathematical simulations based on a population of models representing intersubject variability allow the identification of ionic mechanisms underlying rotor dynamics and the definition of new personalized pharmacologic strategies. Our results suggest that the underlying mechanism of the diverging success of ICaL block as an antiarrhythmic strategy is dependent on the basal availability of sodium and calcium ion channel conductivities.